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.

We propose a revision of the system developed by L'epine et al. (2007) for spectroscopic M subdwarf classification. Based on an analysis of subdwarf spectra and templates from Savcheva et al. (2014), we show thatthe CaH1 feature originally proposed by Gizis (1997) is important in selecting reliable cool subdwarf spectra. This index should be used in combination with the [TiO5, CaH2+CaH3] relation provided by L\'epine et al. (2007) to avoid misclassification results. In the new system, the dwarf-subdwarf separators are first derived from a sample of more than 80,000 M dwarfs and a "labeled" subdwarf subsample, these objects being all visually identified from their optical spectra. Based on these two samples, we re-fit the initial [TiO5, CaH1] relation, and propose a new [CaOH, CaH1] relation supplementing the [TiO5, CaH1] relation to reduce the impact of uncertainty in flux calibration on classification accuracy. In addition, we recalibrate the $\zeta_{TiO/CaH}$ parameter defined in L'epine et al. (2007) to enable its successful application to LAMOST spectra. Using this new system, we select candidates from LAMOST Data Release 4 and finally identify a set of 2791 new M subdwarf stars, covering the spectral sequence from type M0 to M7. This sample contains a large number of objects located at low Galactic latitudes, especially in the Galactic anti-center direction, expanding beyond previously published halo- and thick disk-dominated samples. Besides, we detect magnetic activity in 141 objects. We present a catalog for this M subdwarf sample, including radial velocities, spectral indices and errors, activity flags, with a compilation of external data (photometric and GAIA DR2 astrometric parameters). The catalog is provided on-line, and the spectra can be retrieved from the LAMOST Data Release web portal.

A sample of Cataclysmic Variables (CVs) is presented including spectroscopically identified 380 spectra of 245 objects, of which 58 CV candidates are new discoveries. The BaggingTopPush and the Random Forest algorithms are applied to the Fifth Data Release (DR5) of LAMOST to retrieve CVs with strong emission lines and with broad absorption lines respectively. Based on spectroscopic classification, 134 dwarf novae, 41 nova-like variables and 19 magnetic CVs are identified from the sample. In addition, 89 high--inclination systems and 33 CVs showing companion stars are recognized and discussed for their distinct spectral characteristics. Comparisons between CVs from LAMOST and from published catalogs are made in spatial and magnitude distribution, and the difference of their locus in Gaia color--absolute magnitude diagram (CaMD) are also investigated. More interestingly, for two dwarf novae observed through LAMOST and SDSS in different epoch, their spectra both in quiescence phase and during outburst are exhibited.

LAMOST DR5 released more than 200,000 low resolution spectra of early-type stars with S/N>50. Searching for metallic-line (Am) stars in such a large database and study of their statistical properties are presented in this paper. Six machine learning algorithms were experimented with using known Am spectra, and both the empirical criteria method(Hou et al. 2015) and the MKCLASS package(Gray et al. 2016) were also investigated. Comparing their performance, the random forest (RF) algorithm won, not only because RF has high successful rate but also it can derives and ranks features. Then the RF was applied to the early type stars of DR5, and 15,269 Am candidates were picked out. Manual identification was conducted based on the spectral features derived from the RF algorithm and verified by experts. After manual identification, 9,372 Am stars and 1,131 Ap candidates were compiled into a catalog. Statistical studies were conducted including temperature distribution, space distribution, and infrared photometry. The spectral types of Am stars are mainly between F0 and A4 with a peak around A7, which is similar to previous works. With the Gaia distances, we calculated the vertical height Z from the Galactic plane for each Am star. The distribution of Z suggests that the incidence rate of Am stars shows a descending gradient with increasing jZj. On the other hand, Am stars do not show a noteworthy pattern in the infrared band. As wavelength gets longer, the infrared excess of Am stars decreases, until little or no excess in W1 and W2 bands.

Following previous work, we further confirm that the cosmic evolution of steep-spectrum radio-loud AGNs (active galactic nuclei) can be reproduced by a simple combination of density evolution (DE) and luminosity evolution (LE). This mixture evolution scenario can naturally explain the luminosity-dependent evolution of radio-loud AGNs. Our models successfully fitted a large amount of data on radio luminosity functions (RLFs) of steep-spectrum sources and multi-frequency source counts. The modeling indicates that the DE slowly increase as $(1+z)^{0.3 \thicksim 1.3}$ out to $z \thicksim 0.8$, and then rapidly decreases as $(1+z)^{-6.8 \thicksim -5.7}$, while the LE rapidly increase as $(1+z)^{4.8}$ out to a higher redshift (at least $z>3.5$). We find a high-redshift decline (i.e. redshift cutoff) in the number density of steep-spectrum radio sources, but we cannot conclude whether such decline is sharp or shallow. We believe that whether a redshift cutoff occurs or not depends mainly on DE, while its steepness is decided by LE, which, however, cannot be well constrained due to the lack of high-redshift samples. Most intriguingly, according to our mixture evolution scenario, there appears to be no need for different evolution for the low- and high-power radio-loud AGNs. Both types of sources experience the same combined evolution of DE and LE.

We spectroscopically identify 101 Galactic HII regions using spectra from the Large Sky Area Multi- Object Fiber Spectroscopic Telescope (LAMOST) survey, cross-matched with an HII region catalog derived from the all-sky Wide-Field Infrared Survey Explorer(WISE) data. Among all HII regions in our sample, 47 sources are newly confirmed. Spatially, most of our identified HII regions are located in the anti-center area of the Galaxy. For each of the HII regions, we accurately extract and measure the nebular emission lines of the spectra, and estimate the oxygen abundances using the strong-line method. We focus on the abundance distribution of HII regions in the Galactic anti-center area. Accordingly, we derive the oxygen abundance gradient with a slope of -0.036 +/-0.004 dex/kpc, covering a range of RG from 8.1 to 19.3 kpc. In particular, we also fit the outer disk objects with a slope of -0.039 +/- 0.012 dex /kpc, which indicates that there is no flattening of the radial oxygen gradient in the outer Galactic disk.

In this work, we present a catalog of 2651 carbon stars from the fourth Data Release (DR4) of the Large Sky Area Multi-Object Fiber Spectroscopy Telescope (LAMOST). Using an efficient machine-learning algorithm, we find out these stars from more than seven million spectra. As a by-product, 17 carbon-enhanced metal-poor (CEMP) turnoff star candidates are also reported in this paper, and they are preliminarily identified by their atmospheric parameters. Except for 176 stars that could not be given spectral types, we classify the other 2475 carbon stars into five subtypes including 864 C-H, 226 C-R, 400 C-J, 266 C-N, and 719 barium stars based on a series of spectral features. Furthermore, we divide the C-J stars into three subtypes of CJ( H), C-J(R), C-J(N), and about 90% of them are cool N-type stars as expected from previous literature. Beside spectroscopic classification, we also match these carbon stars to multiple broadband photometries. Using ultraviolet photometry data, we find that 25 carbon stars have FUV detections and they are likely to be in binary systems with compact white dwarf companions.

We fit the spectral energy distributions (SEDs) of 46 GeV - TeV BL Lac objects in the frame of leptonic one-zone synchrotron self-Compton (SSC) model and investigate the physical properties of these objects. We use the Markov Chain Monte Carlo (MCMC) method to obtain the basic parameters, such as magnetic field (B), the break energy of the relativistic electron distribution ($\gamma'_{\rm{b}}$) and the electron energy spectral index. Based on the modeling results, we support the following scenarios on GeV-TeV BL Lac objects: (1) Some sources have large Doppler factors, implying other radiation mechanism should be considered. (2) Comparing with FSRQs, GeV-TeV BL Lac objects have weaker magnetic field and larger Doppler factor, which cause the ineffective cooling and shift the SEDs to higher bands. Their jet powers are around $4.0\times 10^{45}~\rm{ erg\cdot s}^{-1}$, comparing with radiation power, $5.0\times 10^{42}~\rm{ erg\cdot s}^{-1}$, indicating that only a small fraction of jet power is transformed into the emission power. (3) For some BL Lacs with large Doppler factors, their jet components could have two substructures, e.g., the fast core and the slow sheath. For most GeV-TeV BL Lacs, Kelvin-Helmholtz instabilities are suppressed by their higher magnetic fields, leading few micro-variability or intro-day variability in the optical bands. (4) Combined with a sample of FSRQs, an anti-correlation between the peak luminosity $L_{\rm {pk}}$ and the peak frequency $\nu_{\rm {pk}}$ is obtained, favoring the blazar sequence scenario. In addition, an anti-correlation between the jet power $P_{\rm {jet}}$ and the break Lorentz factor $\gamma_{\rm {b}}$ also supports the blazar sequence.

We present a model-dependent method to estimate the redshifts of three TeV BL Lac objects (BL Lacs) through fltting their (quasi-) simultaneous multi-waveband spectral energy distributions (SEDs) by one-zone leptonic synchrotron self-Compton (SSC) model. Considering the impact of electron energy distributions (EEDs) on the results, we use three types of EEDs, such as the power-law EED with exponential cut-ofi (PLC), the log-parabola (PLLP) and the broken power-law (BPL) EEDs, to flt the SEDs. We also use a parameter fi to describe the uncertainties of the extragalactic background light (EBL) models, as in (Abdo et al. 2010b). We then use Markov Chain Monte Carlo (MCMC) method to explore multi-dimensional parameter space and obtain the uncertainties of the model parameters based on the observational data. We apply our method to obtain the redshifts of three TeV BL Lac objects in the marginalized 68% confldence, and flnd that the PLC EED does not flt the SEDs. For 3C66A, the redshift is 0.14 - 0.31 and 0.16 - 0.32 in the BPL and PLLP EEDs; for PKS1424+240, the redshift is 0.55 - 0.68 and 0.55 - 0.67 in the BPL and PLLP EEDs; for PG1553+113, the redshift is 0.22 - 0.48 and 0.22 - 0.39 in the BPL and PLLP EEDs. We also estimate the redshift of PKS1424+240 in the high stage to be 0.46 - 0.67 in the PLLP EED, roughly consistent with that in the low stage.

Systematic variations in 54Cr/52Cr ratios between meteorite classes (Qin et al., 2010a; Trinquier et al., 2007) point to large scale spatial and/or temporal isotopic heterogeneity in the solar protoplanetary disk. Two explanations for these variations have been proposed, with important implications for the formation of the Solar System: heterogeneous seeding of the disk with dust from a supernova, or energetic-particle irradiation of dust in the disk. The key to differentiating between them is identification of the carrier(s) of the 54Cr anomalies. Here we report the results of our recent NanoSIMS imaging search for the 54Cr-rich carrier in the acid-resistant residue of the CI chondrite Orgueil. A total of 10 regions with extreme 54Cr-excesses ({\delta}54Cr values up to 1500 %) were found. Comparison between SEM, Auger and NanoSIMS analyses showed that these 54Cr-rich regions are associated with one or more sub-micron (typically less than 200 nm) Cr oxide grains, most likely spinels. Because the size of the NanoSIMS primary O- ion beam is larger than the typical grain size on the sample mount, the measured anomalies are lower limits, and we estimate that the actual 54Cr enrichments in three grains are at least 11 times Solar and in one of these may be as high as 50 times Solar. Such compositions strongly favor a Type II supernova origin. The variability in bulk 54Cr/52Cr between meteorite classes argues for a heterogeneous distribution of the 54Cr carrier in the solar protoplanetary disk following a late supernova injection event. Such a scenario is also supported by the O-isotopic distribution and variable abundances in different planetary materials of other presolar oxide and silicate grains from supernovae.

The thermal evolution of a neutron star is studied by including the energy loss due to axion emission. Two axion models and three types of neutron-star matter equation of state are used with the effects of nucleon superfluidity properly taken into account. In comparison with the observational data of PSR0656+14 from ROSAT, the upper limits on the axion mass are found to be $m_a < 0.06 - 0.3$ eV and 0.08$ -$ 0.8 eV for the KSVZ and DFSZ axion models, respectively, with the soft equation of state giving the most stringent limits.