Since Sep. 2018, LAMOST started the medium-resolution (R $\sim 7500$) spectral survey (MRS). We proposed the spectral survey of Galactic nebulae, including \hii\ regions, HH objects, supernova remnants, planetary nebulae and the special stars with MRS (LAMOST MRS-N). LAMOST MRS-N covers about 1700 square degrees of the northern Galactic plane within $40^\circ < l < 215^\circ$ and $-5^\circ < b < 5^\circ$. In this 5-year survey, we plan to observe about 500 thousand nebulae spectra. According to the commissioning observations, the nebulae spectra can provide precise radial velocity with uncertainty less than 1 km s$^{-1}$. These high-precision spectral data are of great significance to our understanding of star formation and evolution.

Providing a possible connection between neutrino emission and gravitational-wave (GW) bursts is important to our understanding of the physical processes that occur when black holes or neutron stars merge. In the Daya Bay experiment, using data collected from December 2011 to August 2017, a search has been performed for electron-antineutrino signals coinciding with detected GW events, including GW150914, GW151012, GW151226, GW170104, GW170608, GW170814, and GW170817. We used three time windows of $\mathrm{\pm 10~s}$, $\mathrm{\pm 500~s}$, and $\mathrm{\pm 1000~s}$ relative to the occurrence of the GW events, and a neutrino energy range of 1.8 to 100 MeV to search for correlated neutrino candidates. The detected electron-antineutrino candidates are consistent with the expected background rates for all the three time windows. Assuming monochromatic spectra, we found upper limits (90% confidence level) on electron-antineutrino fluence of $(1.13~-~2.44) \times 10^{11}~\rm{cm^{-2}}$ at 5 MeV to $8.0 \times 10^{7}~\rm{cm^{-2}}$ at 100 MeV for the three time windows. Under the assumption of a Fermi-Dirac spectrum, the upper limits were found to be $(5.4~-~7.0)\times 10^{9}~\rm{cm^{-2}}$ for the three time windows.

We present panchromatic observations and modeling of the Calcium-rich supernova 2019ehk in the star-forming galaxy M100 (d$\approx$16.2 Mpc) starting 10 hours after explosion and continuing for ~300 days. SN 2019ehk shows a double-peaked optical light curve peaking at $t = 3$ and $15$ days. The first peak is coincident with luminous, rapidly decaying $\textit{Swift}$-XRT discovered X-ray emission ($L_x\approx10^{41}~\rm{erg~s^{-1}}$ at 3 days; $L_x \propto t^{-3}$), and a Shane/Kast spectral detection of narrow H$\alpha$ and He II emission lines ($v \approx 500$ km/s) originating from pre-existent circumstellar material. We attribute this phenomenology to radiation from shock interaction with extended, dense material surrounding the progenitor star at $r<10^{15}$ cm and the resulting cooling emission. We calculate a total CSM mass of $\sim$ $7\times10^{-3}$ $\rm{M_{\odot}}$ with particle density $n\approx10^{9}\,\rm{cm^{-3}}$. Radio observations indicate a significantly lower density $n < 10^{4}\,\rm{cm^{-3}}$ at larger radii. The photometric and spectroscopic properties during the second light curve peak are consistent with those of Ca-rich transients (rise-time of $t_r =13.4\pm0.210$ days and a peak B-band magnitude of $M_B =-15.1\pm0.200$ mag). We find that SN 2019ehk synthesized $(3.1\pm0.11)\times10^{-2} ~ \rm{M_{\odot}}$ of ${}^{56}\textrm{Ni}$ and ejected $M_{\rm ej} = (0.72\pm 0.040)~\rm{M_{\odot}}$ total with a kinetic energy $E_{\rm k}=(1.8\pm0.10)\times10^{50}~\rm{erg}$. Finally, deep $\textit{HST}$ pre-explosion imaging at the SN site constrains the parameter space of viable stellar progenitors to massive stars in the lowest mass bin (~10 $\rm{M_{\odot}}$) in binaries that lost most of their He envelope or white dwarfs. The explosion and environment properties of SN 2019ehk further restrict the potential WD progenitor systems to low-mass hybrid HeCO WD + CO WD binaries.

It is generally believed that the variability of photospheric emission in gamma-ray bursts (GRBs) traces that of the jet power. This work further investigates the variability of photospheric emission in a variable jet. By setting a constant $\eta$ (dimensionless entropy of the jet), we find that the light curve of the photospheric emission shows a ``tracking'' pattern on the time profile of jet power. However, the relative variability is significantly low in the photospheric emission compared with that in the jet power. If the $\eta$ is genetic variable, the variability of the photospheric emission is not only limited by the jet power but also affected by $\eta$ strongly. It becomes complex and is generally different from that of the jet power. Moreover, the opposite phase may stand in the variabilities of the photospheric emission at different photon energies. We also find that the relative variability does not remain constant over the photon energies with an obvious reduction at a certain energy. This is consistent with the analysis of GRB 090902B in which an appreciable thermal component has been detected in a wide energy range. For several other GRBs coupling with the thermal component, we conservatively evaluate the variability of the thermal and non-thermal emission, respectively. Our results show that the relative variability of the thermal emission is likely comparable to that of the non-thermal emission for these bursts. In addition, the analysis of GRB~120323A reveals that the variability of the photospheric emission may be of the opposite phase from that of the non-thermal emission.

We present a spectroscopic survey of 230 white dwarf candidates within 40 pc of the Sun from the William Herschel Telescope and Gran Telescopio Canarias. All candidates were selected from $\textit{Gaia}$ Data Release 2 (DR2) and in almost all cases had no prior spectroscopic classifications. We find a total of 191 confirmed white dwarfs and 39 main-sequence star contaminants. The majority of stellar remnants in the sample are relatively cool ($\langle T_{\rm eff} \rangle$ = 6200 K), showing either hydrogen Balmer lines or a featureless spectrum, corresponding to 89 DA and 76 DC white dwarfs, respectively. We also recover two DBA white dwarfs and 9--10 magnetic remnants. We find two carbon-bearing DQ stars and 14 new metal-rich white dwarfs. This includes the possible detection of the first ultra-cool white dwarf with metal lines. We describe three DZ stars for which we find at least four different metal species, including one which is strongly Fe- and Ni-rich, indicative of the accretion of a planetesimal with core-Earth composition. We find one extremely massive (1.31 $\pm$ 0.01 M$_{\odot}$) DA white dwarf showing weak Balmer lines, possibly indicating stellar magnetism. Another white dwarf shows strong Balmer line emission but no infrared excess, suggesting a low-mass sub-stellar companion. High spectroscopic completeness ($>$99%) has now been reached for $\textit{Gaia}$ DR2 sources within 40 pc sample, in the northern hemisphere ($\delta >$ 0 deg) and located on the white dwarf cooling track in the Hertzsprung-Russell diagram. A statistical study of the full northern sample is presented in a companion paper.

Quark matter with only $u$ and $d$ quarks ($ud$QM) might be the ground state of baryonic matter at large baryon number $A>A_{\rm min}$. With $A_{\rm min}\gtrsim 300$, this has no direct conflict with the stability of ordinary nuclei. An intriguing test of this scenario is to look for quantum nucleation of $ud$QM inside neutron stars due to their large baryon densities. In this paper, we study the transition rate of cold neutron stars to $ud$ quark stars ($ud$QSs) and the astrophysical implications, considering the relevant theoretical uncertainties and observational constraints. It turns out that a large portion of parameter space predicts an instantaneous transition, and so the observed neutron stars are mostly $ud$QSs. We find this possibility still viable under the recent gravitational wave and pulsar observations, although there are debates on its compatibility with some observations that involve complicated structure of quark matter. The tension could be partially relieved in the two-families scenario, where the high-mass stars ($M\gtrsim2 M_{\odot}$) are all $ud$QSs and the low-mass ones ($M\sim1.4\, M_{\odot}$) are mostly hadronic stars. In this case, the slow transition of the low-mass hadronic stars points to a very specific class of hadronic models with moderately stiff EOSs, and $ud$QM properties are also strongly constrained.

We compare the infrared excess (IRX) and Balmer decrement (${\rm H\alpha/H\beta }$) as dust attenuation indicators in relation to other galaxy parameters using a sample of $\sim$32 000 local star-forming galaxies (SFGs) carefully selected from SDSS, GALEX and WISE. While at fixed ${\rm H\alpha/H\beta }$, IRX turns out to be independent on galaxy stellar mass, the Balmer decrement does show a strong mass dependence at fixed IRX. We find the discrepancy, parameterized by the color excess ratio $R_{\rm EBV} \equiv E(B-V)_{\rm IRX}/E(B-V)_{\rm H\alpha/H\beta }$, is not dependent on the gas-phase metallicity and axial ratio but on the specific star formation rate (SSFR) and galaxy size ($R_{\rm e}$) following $R_{\rm EBV}=0.79+0.15\log({\rm SSFR}/R_{\rm e}^{2})$. This finding reveals that the nebular attenuation as probed by the Balmer decrement becomes increasingly larger than the global (stellar) attenuation of SFGs with decreasing SSFR surface density. This can be understood in the context of an enhanced fraction of intermediate-age stellar populations that are less attenuated by dust than the HII region-traced young population, in conjunction with a decreasing dust opacity of the diffuse ISM when spreading over a larger spatial extent. Once the SSFR surface density of an SFG is known, the conversion between attenuation of nebular and stellar emission can be well estimated using our scaling relation.

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.

Primordial black holes that survive until the present have been considered as a dark matter candidate. In this paper we argue that primordial 2-2-hole remnants provide a more promising and testable option. 2-2-holes arise in quadratic gravity as a new family of classical solutions for ultracompact matter distributions and they possess the black hole exterior without an event horizon. They may serve as the endpoint of gravitational collapse, providing a resolution for the information loss problem. Intriguing thermodynamic behavior is found for these objects when sourced by a thermal gas. A large 2-2-hole radiates with a Hawking-like temperature and exhibits an entropy-area law. At a late stage, the evaporation slows down and essentially stops as the mass asymptotically approaches a minimal value. This remnant mass is determined by a fundamental scale in quadratic gravity. We study the cosmological and astrophysical implications of having these remnants as dark matter and derive the corresponding constraints. A distinctive phenomenon associated with remnant mergers occurs, predicting fluxes of high-energy astrophysical particles due to the spectacular evaporation of the merger product. Measurements of high-energy photon and neutrino fluxes could possibly bound the remnant mass to be not far above the Planck mass. Early-universe physics, on the other hand, requires that 2-2-holes quickly evolve into the remnant state after formation, putting an upper bound on the formation mass.

We are entering a new era to test the strong gravity regime around astrophysical black holes. The possibility that they are actually horizonless ultracompact objects and then free from the information loss paradox can be examined more closely with observational data. In this paper, we systematically develop a thermal gas model of the 2-2-hole in quadratic gravity, as one step further to look for more tractable models of black hole mimickers. Concrete predictions for departures from black holes are made all the way down to the high curvature interior. The simple form of matter further enables an explicit study of the relation between geometry and thermodynamics. Within this unified framework, we identify notably different behaviors at two limits. On one side is the astrophysically large 2-2-hole, as characterized by a minuscule deviation outside the would-be horizon and a highly squeezed interior along the radial direction. Anomalous features of black hole thermodynamics emerge from the ordinary gas. On the other side is the minimal 2-2-hole with an isotropic and shrinking interior, which behaves more like a normal thermodynamic system. This brings a new perspective to the related theoretical questions as well as phenomenological implications.

The first detected gravitational wave GW170817 from a binary neutron star merger is associated with an important optical transient AT 2017gfo, which is a direct observation of kilonova. Recent observations suggest that the remnant compact object of the binary neutron star merger associated with GW170817/GRB 170817A may be a stable long-lived magnetized neutron star. In this situation, there would be a pulsar wind nebula (PWN) embedded inside the dynamic ejecta. The PWN emission may be absorbed by the ejecta or leak out of the system. We study the effect of the PWN emission on the observed light curves and radiation spectra. Different from previous works, the absorption and leakage of the PWN emission are all involved in our model, where the absorption of the PWN emission heats up the ejecta and alters its radiation. It is found that the characteristic emission of the embedded PWN quickly evolves. For the multiband and long-term observations of AT 2017gfo, we find that the dynamic ejecta with a PWN emission can fit the observational data very well, especially for the light curves at $t\sim 5$ days and those in the late phase. In addition, our model can naturally generate the thermal to nonthermal spectrum evolution of AT 2017gfo. Our fitting result suggests that a PWN is embedded in the AT 2017gfo.

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}$.

We report our photometric and spectroscopic follow-up observations of the superoutbursts of three dwarf novae (GWAC\,180415A, GWAC\,181017A and GWAC\,181211A) identified independently by the Ground Wide-angle Cameras system, one of the ground-based instruments of the China-France SVOM mission. Based on a combination of our photometry and that taken from the AAVSO, our period analysis of the superhumps enables us to determine the mass ratios to be 0.0967-0.1163, 0.1879-0.1883 and 0.0981-0.1173 for GWAC\,180415A, GWAC\,181017A and GWAC\,181211A, respectively. GWAC\,180415A can be firmly identified as a WZ sge-type dwarf novae due to its long duration ($\sim2$ weeks) multiple rebrightenings with amplitudes of 3-4 magnitudes, the early superhump associated with a double-wave modulation and the low mass ratio. The inferred low mass ratio and location in the $\varepsilon-P_{\mathrm{orb}}$ diagram suggest that GWAC\,181211A is a WZ sge-type dwarf novae candidate. The measured Balmer decrements suggest the Balmer line emission is produced from an optical thick region in GWAC\,180415A and GWAC\,181017A, and from an optical thin region in GWAC\,181211A.

We obtained 7 spectra of the Be star V423 Aur on Dec. 5th, 2017 using the LAMOST Medium-Resolution Spectrograph with exposures from 600 to 1200 seconds. These spectra show the irregular H{\alpha} emission line profile variations (LPVs). In the seven spectra, from the 4th to 7th, the left part of H{\alpha} profile even show excess. However, no variation can be seen from the follow-up observation of photometry by 1.26-m telescope and High-Resolution spectra by 2.16-m telescope. According to the High-Resolution spectra, we conclude that it is a B7V type star with E(B-V)= 0.709{\pm}0.036 and its vsini is {\sim} 221.8km/s. The short-term H{\alpha} LPVs could be explained as a result of the transient ejection of matter from rotating disk or shell around V423 Aur.

Recent research shows that the spectral lag is closely related to the spectral evolution in GRBs. In this paper, we study the spectral lag for a radiating jet shell with a high-energy cut-off radiation spectrum. For the jet shell with a cut-off power-law spectrum, the spectral lag monotonically increases with the photon energy and levels off at a certain photon energy. It is the same for the jet shell with a Band cut-off spectrum (Bandcut). However, a turn-over from the positive lags to negative lags appears in the high-energy range for the jet shell with Bandcut, which is very similar to that observed in GRB~160625B. The dependence of the spectral lags on the spectral shape/evolution are studied in details. In addition, the spectral lag behavior observed in GRB~160625B is naturally reproduced based on our theoretical outcome.

We present the revised and high quality CCD time series observations in 2017 and 2018 of the W UMa-type contact binary BH~Cas. The combination of our new determinations of minimum times and those collected from the literatures reveal a long-term period increase with a rate of $\textit{dP/dt}\,=\,+3.27 \times 10^{-7} \textrm{d}\,\textrm{yr}^{-1}$, which can be accounted for by angular momentum transfer as the less-massive component loses its mass to the more massive one. One cyclical oscillation ($A'=0.00300$~days and $P'=20.09$~yrs) superimposed on the long-term period increase tendency is attributed to one additional object that orbits BH Cas. The low-resolution spectra of this source are observed at phase $\sim$0.0 and $\sim$0.5 in 2018 and 2019, and the spectral types are identified as K3V$\,\pm\,$1 and G8V$\,\pm\,$2. According to the luminosity contributions and spectra at these phases, the spectral type of the primary star is close to K3V$\,\pm\,$1. We infer a cool spot on the hotter and less-massive component to fit the asymmetry of light curves. In addition, changes in the location, temperature, and area of spots on the secondary star in different observations may indicate that the magnetic field activity of the secondary is more active than that of the primary. Such strong spot activities and optical flare are both detected at the same time from a W UMa-type star for the first time. With the \emph{XMM-Newton} observations taken in 2003 and 2012, the X-ray light curve is obtained and the X-ray spectra can be well described with a black-body or a bremsstrahlung model. The soft X-ray luminosity of BH Cas is estimated as 8.8\,--\,9.7$\,\times\,10^{30}$\,erg\,s$^{-1}$.

The structured jet is involved to explain the afterglows and even the prompt emission of GRB 170817A. In this paper, we stress that for a precessing jet, the jet structure in the prompt emission phase and that in the afterglow phase may be different. The jet structure in the afterglow phase can be non-uniform even if a narrow-uniform jet is presented in the prompt emission phase. We estimate the jet structure in the afterglow phase under the situation that a narrow-uniform-precessing jet is launched from the central engine of gamma-ray burst. With different precession angles, it is found that the structured jet can be roughly described as follows: a narrow uniform core with power-law wings and sharp cut-off edges, a Gaussian profile, a ring shape, or other complex profile in energy per solid angle. Correspondingly, the afterglows for our obtained structured jets are also estimated. We find that the estimates of the intrinsic kinetic energy, the electron index, and the jet opening angle based on the afterglows formed in a precessing system may be incorrect. Our obtained structured jet is likely to be revealed by future observations for a fraction of gravitational wave detected merging compact binary systems (e.g., black hole-neutron star mergers).

We investigate dust obscuration as parameterised by the infrared excess IRX$\equiv$$L_{\rm IR}/L_{\rm UV}$ in relation to global galaxy properties, using a sample of $\sim$32$\,$000 local star-forming galaxies (SFGs) selected from SDSS, GALEX and WISE. We show that IRX generally correlates with stellar mass ($M_\ast$), star formation rate (SFR), gas-phase metallicity ($Z$), infrared luminosity ($L_{\rm IR}$) and the half-light radius ($R_{\rm e}$). A weak correlation of IRX with axial ratio (b/a) is driven by the inclination and thus seen as a projection effect. By examining the tightness and the scatter of these correlations, we find that SFGs obey an empirical relation of the form $IRX$=$10^\alpha\,(L_{\rm IR})^{\beta}\,R_{\rm e}^{-\gamma}\,(b/a)^{-\delta}$ where the power-law indices all increase with metallicity. The best-fitting relation yields a scatter of $\sim$0.17$\,$dex and no dependence on stellar mass. Moreover, this empirical relation also holds for distant SFGs out to $z=3$ in a population-averaged sense, suggesting it to be universal over cosmic time. Our findings reveal that IRX approximately increases with $L_{\rm IR}/R_{\rm e}^{[1.3 - 1.5]}$ instead of $L_{\rm IR}/R_{\rm e}^{2}$ (i.e., surface density). We speculate this may be due to differences in the spatial extent of stars versus star formation and/or complex star-dust geometries. We conclude that not stellar mass but IR luminosity, metallicity and galaxy size are the key parameters jointly determining dust obscuration in SFGs.

Cool subdwarfs are metal-poor low-mass stars that formed during the early stages of the evolution of our Galaxy. Because they are relatively rare in the vicinity of the Sun, we know of few cool subdwarfs in the solar neighbourhood, and none with both the mass and the radius accurately determined. This hampers our understanding of stars at the low-mass end of the main-sequence. Here we report the discovery of SDSSJ235524.29+044855.7 as an eclipsing binary containing a cool subdwarf star, with a white dwarf companion. From the light-curve and the radial-velocity curve of the binary we determine the mass and the radius of the cool subdwarf and we derive its effective temperature and luminosity by analysing its spectral energy distribution. Our results validate the theoretical mass-radius-effective temperature-luminosity relations for low-mass low-metallicity stars.

We have investigated the metallicity distributions of mono-age stellar populations across the disc of $6 \lesssim R \lesssim 12$\,kpc and $|Z| \lesssim 2$\,kpc using samples selected from the main-sequence turn-off and sub-giant (MSTO-SG) stars targeted by the LAMOST Galactic Spectroscopic surveys. Both the mean values and the profiles of the distributions exhibit significant variations with age and position. We confirm that the oldest ($>11$\,Gyr) stars have nearly flat radial [Fe/H] gradients at all heights above the disc but show negative vertical [Fe/H] gradients. For stars younger than 11\,Gyr, the radial [Fe/H] gradients flatten with $|Z|$, while the vertical [Fe/H] gradients flatten with $R$. Stars of 4--6\,Gyr exhibit steeper negative radial [Fe/H] gradients than those of either younger or older ages. Values of [$\alpha$/Fe] of mono-age stellar populations also show significant radial and vertical gradients, with patterns varying with age. The [Fe/H] distribution profiles of old ($>8$\,Gyr) stars vary little with $R$, while those of younger stars exhibit strong radial variations, probably a consequence of significant radial migration. The [$\alpha$/Fe] radial distribution profiles show opposite patterns of variations with age compared to those of [Fe/H]. We have also explored the impacts of stellar mixing by epicycle motions (blurring) on the [Fe/H] and [$\alpha$/Fe] distributions, and found that blurring mainly change the widths of the distribution profiles. Our results suggest that the disc may have experienced a complex assemblage history, in which both the "inside-out" and "upside-down" formation processes may have played an important role.

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.

In the minimal supersymmetric standard model (MSSM) the lightest superparticle (LSP) can be a TeV neutralino (mainly dominated by higgsino or wino) which serves as a dark matter candidate with correct thermal relic density. In this work we confront the 1-2 TeV neutralino dark matter with the latest direct and indirect detections from PandaX and AMS-02/DAMPE. Considering various scenarios with decoupled sfermions, with \textit{A}-mediated annihilation, with squark or stop coannihilation, we find that the parameter space is stringently constrained by the direct detection limits. In the allowed parameter space, the TeV neutralino dark matter annihilation contribution to the anti-proton flux is found to agree with the AMS-02 data while its contribution to eletron/positron flux is too small to cause any visible excess. The current survived parameter space can be mostly covered by the future direct detection experiment LZ7.2T.

The Beijing Faint Object Spectrograph and Camera (BFOSC) is one of the most important instruments of the 2.16-m telescope of the Xinglong Observatory. Every year there are ~ 20 SCI-papers published based on the observational data of this telescope. In this work, we have systemically measured the total efficiency of the BFOSC of the 2.16-m reflector, based on the observations of two ESO flux standard stars. We have obtained the total efficiencies of the BFOSC instrument of different grisms with various slit widths in almost all ranges, and analysed the factors which effect the efficiency of telescope and spectrograph. For the astronomical observers, the result will be useful for them to select a suitable slit width, depending on their scientific goals and weather conditions during the observation; For the technicians, the result will help them systemically find out the real efficiency of telescope and spectrograph, and further to improve the total efficiency and observing capacity of the telescope technically.

We present extensive spectroscopic observations of supernova remnant (SNR) S147 collected with the Large sky Area Multi-Object fiber Spectroscopic Telescope (LAMOST). The spectra were carefully sky-subtracted taking into account the complex filamentary structure of S147. We have utilized all available LAMOST spectra toward S147, including sky and stellar spectra. By measuring the prominent optical emission lines including H{\alpha}, [NII]{\lambda}6584, and [SII]{\lambda}{\lambda}6717,6731, we present maps of radial velocity and line intensity ratio covering the whole nebula of S147 with unprecedented detail. The maps spatially correlated well with the complex filamentary structure of S147. For the central 2 deg of S147, the radial velocity varies from -100 to 100 km/s and peaks between ~ 0 and 10 km/s. The intensity ratios of H{\alpha}/[SII]{\lambda}{\lambda}6717,6731, [SII]{\lambda}6717/{\lambda}6731 and H{\alpha}/[NII]{\lambda}6584 peak at about 0.77, 1.35 and 1.48, respectively, with a scatter of 0.17, 0.19 and 0.37, respectively. The intensity ratios are consistent with the literature values. However, the range of variations of line intensity ratios estimated here and representative of the whole nebula, are larger than previously estimated.

We present the data release (DR) 5 catalogue of white dwarf-main sequence (WDMS) binaries from the Large Area Multi-Object fiber Spectroscopic Telescope (LAMOST). The catalogue contains 876 WDMS binaries, of which 757 are additions to our previous LAMOST DR1 sample and 357 are systems that have not been published before. We also describe a LAMOST-dedicated survey that aims at obtaining spectra of photometrically-selected WDMS binaries from the Sloan Digital Sky Survey (SDSS) that are expected to contain cool white dwarfs and/or early type M dwarf companions. This is a population under-represented in previous SDSS WDMS binary catalogues. We determine the stellar parameters (white dwarf effective temperatures, surface gravities and masses, and M dwarf spectral types) of the LAMOST DR5 WDMS binaries and make use of the parameter distributions to analyse the properties of the sample. We find that, despite our efforts, systems containing cool white dwarfs remain under-represented. Moreover, we make use of LAMOST DR5 and SDSS DR14 (when available) spectra to measure the Na I {\lambda}{\lambda} 8183.27, 8194.81 absorption doublet and/or H{\alpha} emission radial velocities of our systems. This allows identifying 128 binaries displaying significant radial velocity variations, 76 of which are new. Finally, we cross-match our catalogue with the Catalina Surveys and identify 57 systems displaying light curve variations. These include 16 eclipsing systems, two of which are new, and nine binaries that are new eclipsing candidates. We calculate periodograms from the photometric data and measure (estimate) the orbital periods of 30 (15) WDMS binaries.

We present the observed H$\alpha$ flux and derived star formation rates (SFRs) for a fall sample of low$-$surface$-$brightness galaxies (LSBGs). The sample is selected from the fall sky region of the 40$\%$ ALFALFA {\ion{H}{1}} survey $-$ SDSS DR7 photometric data, and all the $H\alpha$ images were obtained using the 2.16 m telescope, operated by the National Astronomy Observatories, Chinese Academy of Sciences. A total of 111 LSBGs were observed and $H\alpha$ flux was measured in 92 of them. Though almost all the LSBGs in our sample are {\ion{H}{1}}$-$rich, their SFRs derived from the extinction and filter$-$transmission$-$corrected $H\alpha$ flux, are less than 1$M_{\sun}$$yr^{-1}$. LSBGs and star forming galaxies have similar {\ion{H}{1}} surface densities, but LSBGs have much lower SFRs and SFR surface densities than star$-$forming galaxies. Our results show that LSBGs deviate from the Kennicutt-Schmidt law significantly, which indicate that they have low star formation efficiency. The SFRs of LSBGs are close to average SFRs in Hubble time and support the previous arguments that most of the LSBGs are stable systems and they tend to seldom contain strong interactions or major mergers during their star formation histories.

Stellar systems composed of single, double, triple or high-order systems are rightfully regarded as the fundamental building blocks of the Milky Way. Binary stars play an important role in formation and evolution of the Galaxy. Through comparing the radial velocity variations from multi-epoch observations, we analyze the binary fraction of dwarf stars observed with the LAMOST. Effects of different model assumptions such as orbital period distributions on the estimate of binary fractions, are investigated. The results based on log-normal distribution of orbital periods reproduce the previous complete analyses better than the power-law distribution. We find that the binary fraction increases with $T_{\rm eff}$ and decreases with [Fe/H]. We first investigate the relation between $\alpha$-elements and binary fraction in such a large sample as the LAMOST. The old stars with high [$\alpha$/Fe] dominate higher binary fraction than young stars with low [$\alpha$/Fe]. At the same mass, former forming stars possess a higher binary fraction than newly forming ones, which may be related with the evolution of the Galaxy.

We have studied the variability of \object{S5 0716+714} at radio 15 GHz and $\gamma$-ray band using three different methods. A possible periodicity of $P_{15 GHz}=266.0\pm11.5$ and $P_{\gamma}=344.0\pm16.4$ days are obtained for radio 15 GHz and $\gamma$-ray light curves, respectively. The variability may be related to the intrinsically emission mechanism. The difference between the variability timescales of radio 15 GHz and $\gamma$-ray may be due to that the emission of radio 15 GHz is produced via the synchrotron process, while the $\gamma$-ray is produced by both the SSC and EC processes.

We present a three dimensional (3D) extinction analysis in the region toward the supernova remnant (SNR) S147 (G180.0-1.7) using multi-band photometric data from the Xuyi Schmidt Telescope Photometric Survey of the Galactic Anticentre (XSTPS-GAC), 2MASS and WISE. We isolate a previously unrecognised dust structure likely to be associated with SNR S147. The structure, which we term as "S147 dust cloud", is estimated to have a distance $d$ = 1.22 $\pm$ 0.21 kpc, consistent with the conjecture that S147 is associated with pulsar PSR J0538 + 2817. The cloud includes several dense clumps of relatively high extinction that locate on the radio shell of S147 and coincide spatially with the CO and gamma-ray emission features. We conclude that the usage of CO measurements to trace the SNR associated MCs is unavoidably limited by the detection threshold, dust depletion, and the difficulty of distance estimates in the outer Galaxy. 3D dust extinction mapping may provide a better way to identify and study SNR-MC interactions.

We present the second paper of a series of publications aiming at obtaining a better understanding regarding the nature of type Ia supernovae (SNIa) progenitors by studying a large sample of detached F, G and K main sequence stars in close orbits with white dwarf companions (i.e. WD+FGK binaries). We employ the LAMOST (Large Sky Area Multi-Object Fibre Spectroscopic Telescope) data release 4 spectroscopic data base together with GALEX (Galaxy Evolution Explorer) ultraviolet fluxes to identify 1,549 WD+FGK binary candidates (1,057 of which are new), thus doubling the number of known sources. We measure the radial velocities of 1,453 of these binaries from the available LAMOST spectra and/or from spectra obtained by us at a wide variety of different telescopes around the globe. The analysis of the radial velocity data allows us to identify 24 systems displaying more than 3sigma radial velocity variation that we classify as close binaries. We also discuss the fraction of close binaries among WD+FGK systems, which we find to be ~10 per cent, and demonstrate that high-resolution spectroscopy is required to efficiently identify double-degenerate SNIa progenitor candidates.

We present estimates of stellar age and mass for 0.93 million Galactic disk main sequence turn-off and sub-giant stars from the LAMOST Galactic Spectroscopic Surveys. The ages and masses are determined by matching with stellar isochrones using Bayesian algorithm, utilizing effective temperature $T_{\rm eff}$, absolute magnitude ${\rm M}_V$, metallicity [Fe/H] and $\alpha$-element to iron abundance ratio [$\alpha$/Fe] deduced from the LAMOST spectra. Extensive examinations suggest the age and mass estimates are robust. The overall sample stars have a median error of 34 per cent for the age estimates, and half of the stars older than 2\,Gyr have age uncertainties of only 20--30 per cent. Median error for the mass estimates of the whole sample stars is $\sim8$ per cent. The huge dataset demonstrates good correlations among stellar age, [Fe/H] ([$\alpha$/H]) and [$\alpha$/Fe]. Particularly, double sequence features are revealed in the both the age--[$\alpha$/Fe] and age--[Fe/H]([$\alpha$/H]) spaces. In the [Fe/H]--[$\alpha$/Fe] space, stars of 8--10\,Gyr exhibit both the thin and thick disk sequences, while younger (older) stars show only the thin (thick) disk sequence, indicating that the thin disk became prominent 8--10\,Gyr ago, while the thick disk formed earlier and almost quenched 8\,Gyr ago. Stellar ages exhibit positive vertical and negative radial gradients across the disk, and the outer disk of $R\gtrsim$\,9\,kpc exhibits a strong flare in stellar age distribution.

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.

Accurate determination of stellar atmospheric parameters and elemental abundances is crucial for Galactic archeology via large-scale spectroscopic surveys. In this paper, we estimate stellar atmospheric parameters -- effective temperature T_{\rm eff}, surface gravity log g and metallicity [Fe/H], absolute magnitudes M_V and M_{Ks}, {\alpha}-element to metal (and iron) abundance ratio [{\alpha}/M] (and [{\alpha}/Fe]), as well as carbon and nitrogen abundances [C/H] and [N/H] from the LAMOST spectra with amultivariate regressionmethod based on kernel-based principal component analysis, using stars in common with other surveys (Hipparcos, Kepler, APOGEE) as training data sets. Both internal and external examinations indicate that given a spectral signal-to-noise ratio (SNR) better than 50, our method is capable of delivering stellar parameters with a precision of ~100K for Teff, ~0.1 dex for log g, 0.3 -- 0.4mag for M_V and M_{Ks}, 0.1 dex for [Fe/H], [C/H] and [N/H], and better than 0.05 dex for [{\alpha}/M] ([{\alpha}/Fe]). The results are satisfactory even for a spectral SNR of 20. The work presents first determinations of [C/H] and [N/H] abundances from a vast data set of LAMOST, and, to our knowledge, the first reported implementation of absolute magnitude estimation directly based on the observed spectra. The derived stellar parameters for millions of stars from the LAMOST surveys will be publicly available in the form of value-added catalogues.

Photometric data from the Xuyi Schmidt Telescope Photometric Survey of the Galactic Anticentre (XSTPS-GAC) and the Sloan Digital Sky Survey (SDSS) are used to derive the global structure parameters of the smooth components of the Milky Way. The data, which cover nearly 11,000 deg$^2$ sky area and the full range of Galactic latitude, allow us to construct a globally representative Galactic model. The number density distribution of Galactic halo stars is fitted with an oblate spheroid that decays by power law. The best-fit yields an axis ratio and a power law index $\kappa=0.65$ and $p=2.79$, respectively. The $r$-band differential star counts of three dwarf samples are then fitted with a Galactic model. The best-fit model yielded by a Markov Chain Monte Carlo analysis has thin and thick disk scale heights and lengths of $H_{1}=$ 322\,pc and $L_{1}=$2343\,pc, $H_{2}=$794\,pc and $L_{2}=$3638\,pc, a local thick-to-thin disk density ratio of $f_2=$11\,per\,cent, and a local density ratio of the oblate halo to the thin disk of $f_h=$0.16\,per\,cent. The measured star count distribution, which is in good agreement with the above model for most of the sky area, shows a number of statistically significant large scale overdensities, including some of the previously known substructures, such as the Virgo overdensity and the so-called "north near structure", and a new feature between 150\degr $< l < $ 240\degr~and $-1$5\degr $< b < $ $-$5\degr, at an estimated distance between 1.0 and 1.5\,kpc. The Galactic North-South asymmetry in the anticentre is even stronger than previously thought.

The number of spatially unresolved white dwarf plus main-sequence star binaries has increased rapidly in the last decade, jumping from only ~30 in 2003 to over 3000. However, in the majority of known systems the companion to the white dwarf is a low mass M dwarf, since these are relatively easy to identify from optical colours and spectra. White dwarfs with more massive FGK type companions have remained elusive due to the large difference in optical brightness between the two stars. In this paper we identify 934 main-sequence FGK stars from the Radial Velocity Experiment (RAVE) survey in the southern hemisphere and the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) survey in the northern hemisphere, that show excess flux at ultraviolet wavelengths which we interpret as the likely presence of a white dwarf companion. We obtained Hubble Space Telescope ultraviolet spectra for nine systems which confirmed that the excess is indeed caused, in all cases, by a hot compact companion, eight being white dwarfs and one a hot subdwarf or pre-helium white dwarf, demonstrating that this sample is very clean. We also address the potential of this sample to test binary evolution models and type Ia supernovae formation channels.

The rotation curve (RC) of the Milky Way out to $\sim$ 100 kpc has been constructed using $\sim$ 16,000 primary red clump giants (PRCGs) in the outer disk selected from the LSS-GAC and the SDSS-III/APOGEE survey, combined with $\sim$ 5700 halo K giants (HKGs) selected from the SDSS/SEGUE survey. To derive the RC, the PRCG sample of the warm disc population and the HKG sample of halo stellar population are respectively analyzed using a kinematical model allowing for the asymmetric drift corrections and re-analyzed using the spherical Jeans equation along with measurements of the anisotropic parameter $\beta$ currently available. The typical uncertainties of RC derived from the PRCG and HKG samples are respectively 5-7 km/s and several tens km/s. We determine a circular velocity at the solar position, $V_c (R_0)$ = 240 $\pm$ 6 km/s and an azimuthal peculiar speed of the Sun, $V_{\odot}$ = 12.1 $\pm$ 7.6 km/s, both in good agreement with the previous determinations. The newly constructed RC has a generally flat value of 240 km/s within a Galactocentric distance $r$ of 25 kpc and then decreases steadily to 150 km/s at $r$ $\sim$ 100 kpc. On top of this overall trend, the RC exhibits two prominent localized dips, one at $r$ $\sim$ 11 kpc and another at $r$ $\sim$ 19 kpc. From the newly constructed RC, combined with other constraints, we have built a parametrized mass model for the Galaxy, yielding a virial mass of the Milky Way's dark matter halo of $0.90^{+0.07}_{-0.08} \times 10^{12}$ ${\rm M}_{\odot}$ and a local dark matter density, $\rho_{\rm \odot, dm} = 0.32^{+0.02}_{-0.02}$ GeV cm$^{-3}$.

Constrained flavor violation is a recent proposal for predicting the down-quark Yukawa matrix in terms of those for up quarks and charged leptons. We study the viability of CFV with respect to its predictions for the lepton mass ratios, showing that this remains a challenge, and suggest some possible means for improving this shortcoming. We then extend CFV to include neutrinos, and show that it leads to interesting predictions for hierachical heavy neutrinos, and leptogenesis dominated by decays of the second heaviest one ("N2 leptogenesis"), as well as the possibility of low-scale leptoquark-mediated exotic decays.

We present an updated version of the spectroscopic catalogue of white dwarf-main sequence (WDMS) binaries from the Sloan Digital Sky Survey (SDSS). We identify 939 WDMS binaries within the data releases (DR) 9-12 of SDSS plus 40 objects from DR 1-8 that we missed in our previous works, 646 of which are new. The total number of spectroscopic SDSS WDMS binaries increases to 3294. This is by far the largest and most homogeneous sample of compact binaries currently available. We use a decomposition/fitting routine to derive the stellar parameters of all systems identified here (white dwarf effective temperatures, surface gravities and masses, and secondary star spectral types). The analysis of the corresponding stellar parameter distributions shows that the SDSS WDMS binary population is seriously affected by selection effects. We also measure the NaI 8183.27, 8194.81 absorption doublet and Halpha emission radial velocities (RV) from all SDSS WDMS binary spectra identified in this work. 98 objects are found to display RV variations, 62 of which are new. The RV data are sufficient enough to estimate the orbital periods of three close binaries.

In recent decades, inorganic crystals have been widely used in dark matter direct search experiments. To contribute to the understanding of the capabilities of $CsI(Na)$ and $CaF_{2}(Eu)$ crystals, a mono-energetic neutron beam is utilized to study the properties of nuclear recoils, which are expected to be similar to signals of dark matter direct detection. The quenching factor of nuclear recoils in $CsI(Na)$ and $CaF_{2}(Eu)$, as well as an improved discrimination factor between nuclear recoils and ${\gamma}$ backgrounds in $CsI(Na)$, are reported.

We introduce a method for identifying "twin" Type Ia supernovae, and using them to improve distance measurements. This novel approach to Type Ia supernova standardization is made possible by spectrophotometric time series observations from the Nearby Supernova Factory (SNfactory). We begin with a well-measured set of supernovae, find pairs whose spectra match well across the entire optical window, and then test whether this leads to a smaller dispersion in their absolute brightnesses. This analysis is completed in a blinded fashion, ensuring that decisions made in implementing the method do not inadvertently bias the result. We find that pairs of supernovae with more closely matched spectra indeed have reduced brightness dispersion. We are able to standardize this initial set of SNfactory supernovae to 0.083 +/- 0.012 magnitudes, implying a dispersion of 0.072 +/- 0.010 magnitudes in the absence of peculiar velocities. We estimate that with larger numbers of comparison SNe, e.g, using the final SNfactory spectrophotometric dataset as a reference, this method will be capable of standardizing high-redshift supernovae to within 0.06-0.07 magnitudes. These results imply that at least 3/4 of the variance in Hubble residuals in current supernova cosmology analyses is due to previously unaccounted-for astrophysical differences among the supernovae

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.

In Gentile Fusillo et al. (2015) we developed a selection method for white dwarf candidates which makes use of photometry, colours and proper motions to calculate a probability of being a white dwarf (Pwd). The application of our method to the Sloan Digital Sky Survey (SDSS) data release 10 resulted in nearly 66,000 photometrically selected objects with a derived Pwd, approximately 21000 of which are high confidence white dwarf candidates. Here we present an independent test of our selection method based on a sample of spectroscopically confirmed white dwarfs from the LAMOST (Large Sky Area Multi-Fiber Spectroscopic Telescope) survey. We do this by cross matching all our $\sim$66,000 SDSS photometric white dwarf candidates with the over 4 million spectra available in the third data release of LAMOST. This results in 1673 white dwarf candidates with no previous SDSS spectroscopy, but with available LAMOST spectra. Among these objects we identify 309 genuine white dwarfs. We find that our Pwd can efficiently discriminate between confirmed LAMOST white dwarfs and contaminants. Our white dwarf candidate selection method can be applied to any multi-band photometric survey and in this work we conclusively confirm its reliability in selecting white dwarfs without recourse to spectroscopy. We also discuss the spectroscopic completeness of white dwarfs in LAMOST, as well as deriving effective temperatures, surface gravities and masses for the hydrogen-rich atmosphere white dwarfs in the newly identified LAMOST sample.

We have compared the stellar parameters, temperature, gravity and metallicity, between the LAMOST-DR2 and SDSS-DR12/APOGEE database for stars in common. It is found that the LAMOST database provides a better red-clump feature than the APOGEE database in the Teff versus logg diagram. With this advantage, we have separated red clump stars from red giant stars, and attempt to establish the calibrations between the two datasets for the two groups of stars respectively. It shows that there is a good consistency in temperature with a calibration close to the one-to-one line, and we can establish a satisfied metallicity calibration of [Fe/H]_APO=1.18[Fe/H]_APO+0.11 with a scatter of 0.08 dex for both red clump and red giant branch samples. For gravity, there is no any correlation for red clump stars between the two databases, and scatters around the calibrations of red giant stars are substantial. We found two main sources of the scatters of logg for red giant stars. One is a group of stars with $0.00253*Teff-8.67<logg<2.6 locating at the forbidden region, and the other is the contaminated red clump stars, which could be picked out from the unmatched region where stellar metallicity is not consistent with the position in the Teff versus logg diagram. After excluding stars at the two regions, we have established two calibrations for red giant stars, logg_APO=0.000615*Teff_LAMO+0.697*logg_LAMO-2.208$ (sigma=0.150) for [Fe/H]>-1 and logg_APO=0.000874*Teff_LAMO+0.588*logg_LAMO-3.117 (sigma=0.167) for [Fe/H]<-1. The calibrations are valid for stars with Teff=3800-5400 K and logg=0-3.8 dex, and are useful in a work aiming to combine the LAMOST and APOGEE databases in the future study. In addition, we find that an SVM method based on seismic logg is a good way to greatly improve the accuracy of gravity for these two regions at least in the LAMOST database.

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 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 on the search for new eclipsing white dwarf plus main-sequence (WDMS) binaries in the light curves of the Catalina surveys. We use a colour selected list of almost 2000 candidate WDMS systems from the Sloan Digital Sky Survey, specifically designed to identify WDMS systems with cool white dwarfs and/or early M type main-sequence stars. We identify a total of 17 eclipsing systems, 14 of which are new discoveries. We also find 3 candidate eclipsing systems, 2 main-sequence eclipsing binaries and 22 non-eclipsing close binaries. Our newly discovered systems generally have optical fluxes dominated by the main-sequence components, which have earlier spectral types than the majority of previously discovered eclipsing systems. We find a large number of ellipsoidally variable binaries with similar periods, near 4 hours, and spectral types M2--3, which are very close to Roche-lobe filling. We also find that the fraction of eclipsing systems is lower than found in previous studies and likely reflects a lower close binary fraction among WDMS binaries with early M-type main-sequence stars due to their enhanced angular momentum loss compared to fully convective late M type stars, hence causing them to become cataclysmic variables quicker and disappear from the WDMS sample. Our systems bring the total number of known detached, eclipsing WDMS binaries to 71.

So far all evidences of dark matter (DM) come from astrophysical and cosmological observations, due to gravitational interactions of the DM. It is possible that the true DM particle in the universe joins gravitational interactions only, but nothing else. Such a Gravitational DM (GDM) acts as a weakly interacting massive particle (WIMP), which is conceptually simple and attractive. In this work, we explore this direction by constructing the simplest scalar GDM particle $\chi_s$. It is a $Z_2$ odd singlet under the standard model (SM) gauge group, and naturally joins the unique dimension-4 interaction with Ricci curvature, $\xi_s \chi_s^2 R$, where $\xi_s$ is the dimensionless nonminimal coupling. We demonstrate that this gravitational interaction $\xi_s \chi_s^2 R$, together with Higgs-curvature nonminimal coupling term $\xi_h H^\dag H R$, induces effective couplings between $\chi_s^2$ and SM fields which can account for the observed DM thermal relic abundance. We analyze the annihilation cross sections of GDM particles and derive the viable parameter space for realizing the DM thermal relic density. We further study the direct/indirect detections and the collider signatures of such a scalar GDM. These turn out to be highly predictive and testable.

In this work, we examine the possibility of realizing a strongly first-order electroweak phase transition within the minimal classically scale invariant extension of the standard model (SM), previously proposed and analyzed as a potential solution to the hierarchy problem. By introducing one complex singlet scalar and three right-handed Majorana neutrinos, the scenario was successfully capable of achieving a radiative breaking of the electroweak symmetry (Coleman-Weinberg Mechanism), inducing non-zero masses for the SM neutrinos (seesaw mechanism), presenting a pseudoscalar dark matter candidate, and predicting the existence of a second $CP$-even boson in addition to the 125 GeV scalar. We construct the full finite-temperature one-loop effective potential of the model, including the resummed thermal daisy loops, and demonstrate that finite-temperature effects induce a first-order electroweak phase transition. Requiring the thermally-driven first-order phase transition to be sufficiently strong further constrains the model's parameter space; in particular, an $\mathcal O(0.01)$ fraction of the dark matter in the universe may be simultaneously accommodated with a strongly first-order electroweak phase transition. Moreover, such a phase transition disfavors right-handed Majorana neutrino masses above several hundreds of GeV, confines the pseudoscalar dark matter masses to $\sim 1$-2 TeV, predicts the mass of the second $CP$-even scalar to be $\sim 100$-300 GeV, and requires the mixing angle between the $CP$-even components of the SM doublet and the complex singlet to lie within the range $0.2 \lesssim \sin\omega \lesssim 0.4$. The obtained results are displayed in comprehensive exclusion plots, identifying the viable regions of the parameter space. Many of these predictions lie within the reach of the next LHC run.

The Spectral Galaxy Pairs (SGPs) are defined as the composite galaxy spectra which contain two independent redshift systems. These spectra are useful for studying dust properties of the foreground galaxies. In this paper, a total of 165 spectra of SGPs are mined out from Sloan Digital Sky Survey (SDSS) Data Release 9 (DR9) using the concept of membership degree from the fuzzy set theory particularly defined to be suitable for fuzzily identifying emission lines. The spectra and images of this sample are classified according to the membership degree and their image features, respectively. Many of these 2nd redshift systems are too small or too dim to select from the SDSS images alone, making the sample a potentially unique source of information on dust effects in low-luminosity or low-surface-brightness galaxies that are underrepresented in morphological pair samples. The dust extinction of the objects with high membership degree is also estimated by Balmer decrement. Additionally, analyses for a series of spectroscopic observations of one SGP from 165 systems indicate that a newly star-forming region of our Milky Way might occur.

Context. White dwarf-main sequence (WDMS) binaries are used to study several different important open problems in modern astrophysics. Aims. The Sloan Digital Sky Survey (SDSS) identified the largest catalogue of WDMS binaries currently known. However, this sample is seriously affected by selection effects and the population of systems containing cool white dwarfs and early-type companions is under-represented.Here we search for WDMS binaries within the spectroscopic data release 1 of the LAMOST (Large sky Area Multi-Object fiber Spectroscopic Telescope) survey. LAMOST and SDSS follow different target selection algorithms. Hence, LAMOST WDMS binaries may be drawn from a different parent population and thus help in overcoming the selection effects incorporated by SDSS on the current observed population. Methods. We develop a fast and efficient routine based on the wavelet transform to identify LAMOST WDMS binaries containing a DA white dwarf and a M dwarf companion, and apply a decomposition/fitting routine to their LAMOST spectra to estimate their distances and measure their stellar parameters, namely the white dwarf effective temperatures, surface gravities and masses, and the secondary star spectral types. Results. We identify 121 LAMOST WDMS binaries, 80 of which are new discoveries, and estimate the sample to be \sim90 per cent complete. The LAMOST and SDSS WDMS binaries are found to be statistically different. However, this result is not due to the different target selection criteria of both surveys, but likely a simple consequence of the different observing conditions. Thus, the LAMOST population is found at considerably shorter distances (\sim50-450 pc) and is dominated by systems containing early-type companions and hot white dwarfs. (abridged)

The China Dark Matter Experiment reports results on light WIMP dark matter searches at the China Jinping Underground Laboratory with a germanium detector array with a total mass of 20 g. The physics threshold achieved is 177 eVee ("ee" represents electron equivalent energy) at 50% signal efficiency. With 0.784 kg-days of data, exclusion region on spin-independent coupling with the nucleon is derived, improving over our earlier bounds at WIMP mass less than 4.6 GeV.

The China Dark Matter Experiment reports results on light WIMP dark matter searches at the China Jinping Underground Laboratory with a germanium detector array with a total mass of 20 g. The physics threshold achieved is 177 eVee ("ee" represents electron equivalent energy) at 50% signal efficiency. With 0.784 kg-days of data, exclusion region on spin-independent coupling with the nucleon is derived, improving over our earlier bounds at WIMP mass less than 4.6 GeV.

In a previous work, a classically scale invariant extension of the standard model was proposed, as a potential candidate for resolving the hierarchy problem, by minimally introducing a complex gauge singlet scalar, and generating radiative electroweak symmetry breaking by means of the Coleman- Weinberg Mechanism. Postulating the singlet sector to respect the CP-symmetry, the existence of a stable pseudoscalar dark matter candidate with a mass in the TeV range was demonstrated. More- over, the model predicted the presence of another physical CP-even Higgs boson (with suppressed tree-level couplings), in addition to the 125 GeV scalar discovered by the LHC. The viable region of the parameter space was determined by various theoretical and experimental considerations. In this work, we continue to examine the phenomenological implications of the proposed minimal sce- nario by considering the constraints from the dark matter relic density, as determined by the Planck collaboration, as well as the direct detection bounds from the LUX experiment. Furthermore, we investigate the implications of the collider Higgs searches for the additional Higgs boson. Our results are comprehensively demonstrated in unified exclusion plots, which analyze the viable region of the parameter space from all relevant angles, demonstrating the testability of the proposed scenario.

We study gravitational interaction of Higgs boson through the unique dimension-4 operator $\xi H^\dag H R$, with $H$ the Higgs doublet and $R$ the Ricci scalar curvature. We analyze the effect of this dimensionless nonminimal coupling $\xi$ on weak gauge boson scattering in both Jordan and Einstein frames. We explicitly establish the longitudinal-Goldstone boson equivalence theorem with nonzero $\xi$ coupling in both frames, and analyze the unitarity constraints. We study the $\xi$-induced weak boson scattering cross sections at O(1-30)TeV scales, and propose to probe the Higgs-gravity coupling via weak boson scattering experiments at the LHC(14TeV) and the next generation pp colliders (50-100TeV). We further extend our study to Higgs inflation, and quantitatively derive the perturbative unitarity bounds via coupled channel analysis, under large field background at the inflation scale. We analyze the unitarity constraints on the parameter space in both the conventional Higgs inflation and the improved models in light of the recent BICEP2 data.

We present a set of white dwarf-main sequence (WDMS) binaries identified spectroscopically from the Large sky Area Multi-Object fiber Spectroscopic Telescope (LAMOST, also called the Guo Shou Jing Telescope) pilot survey. We develop a color selection criteria based on what is so far the largest and most complete Sloan Digital Sky Survey (SDSS) DR7 WDMS binary catalog and identify 28 WDMS binaries within the LAMOST pilot survey. The primaries in our binary sample are mostly DA white dwarfs except for one DB white dwarf. We derive the stellar atmospheric parameters, masses, and radii for the two components of 10 of our binaries. We also provide cooling ages for the white dwarf primaries as well as the spectral types for the companion stars of these 10 WDMS binaries. These binaries tend to contain hot white dwarfs and early-type companions. Through cross-identification, we note that nine binaries in our sample have been published in the SDSS DR7 WDMS binary catalog. Nineteen spectroscopic WDMS binaries identified by the LAMOST pilot survey are new. Using the 3$\sigma$ radial velocity variation as a criterion, we find two post-common-envelope binary candidates from our WDMS binary sample.

The China Dark matter Experiment collaboration reports the first experimental limit on WIMP dark matter from 14.6 kg-day of data taken with a 994 g p-type point-contact germanium detector at the China Jinping underground Laboratory where the rock overburden is more than 2400 m. The energy threshold achieved was 400 eVee. According to the 14.6 kg-day live data, we placed the limit of N= 1.75 * 10^{-40} cm^{2} at 90% confidence level on the spin-independent cross-section at WIMP mass of 7 GeV before differentiating bulk signals from the surface backgrounds.

The China Dark matter Experiment collaboration reports the first experimental limit on WIMP dark matter from 14.6 kg-day of data taken with a 994 g p-type point-contact germanium detector at the China Jinping underground Laboratory where the rock overburden is more than 2400 m. The energy threshold achieved was 400 eVee. According to the 14.6 kg-day live data, we placed the limit of N= 1.75 * 10^{-40} cm^{2} at 90% confidence level on the spin-independent cross-section at WIMP mass of 7 GeV before differentiating bulk signals from the surface backgrounds.

We present a spectroscopic catalog of 67082 M dwarfs from the LAMOST pilot survey. For each spectrum of the catalog, spectral subtype, radial velocity, equivalent width of H${\alpha}$, a number of prominent molecular band indices and the metal sensitive parameter $\zeta$ are provided . Spectral subtype have been estimated by a remedied Hammer program (Original Hammer: Covey et al. 2007), in which indices are reselected to obtain more accurate auto-classified spectral subtypes. All spectra in this catalog have been visually inspected to confirm the spectral subtypes. Radial velocities have been well measured by our developed program which uses cross-correlation method and estimates uncertainty of radial velocity as well. We also examine the magnetic activity properties of M dwarfs traced by H${\alpha}$ emission line. The molecular band indices included in this catalog are temperature or metallicity sensitive and can be used for future analysis of the physical properties of M dwarfs. The catalog is available on the website \url{http://sciwiki.lamost.org/MCatalogPilot/}.

China JinPing underground Laboratory (CJPL) is the deepest underground laboratory presently running in the world. In such a deep underground laboratory, the cosmic ray flux is a very important and necessary parameter for rare event experiments. A plastic scintillator telescope system has been set up to measure the cosmic ray flux. The performance of the telescope system has been studied using the cosmic ray on the ground laboratory near CJPL. Based on the underground experimental data taken from November 2010 to December 2011 in CJPL, which has effective live time of 171 days, the cosmic ray muon flux in CJPL is measured to be (2.0+-0.4)*10^(-10)/(cm^2)/(s). The ultra-low cosmic ray background guarantees CJPL's ideal environment for dark matter experiment.

The CDEX Collaboration has been established for direct detection of light dark matter particles, using ultra-low energy threshold p-type point-contact germanium detectors, in China JinPing underground Laboratory (CJPL). The first 1 kg point-contact germanium detector with a sub-keV energy threshold has been tested in a passive shielding system located in CJPL. The outputs from both the point-contact p+ electrode and the outside n+ electrode make it possible to scan the lower energy range of less than 1 keV and at the same time to detect the higher energy range up to 3 MeV. The outputs from both p+ and n+ electrode may also provide a more powerful method for signal discrimination for dark matter experiment. Some key parameters, including energy resolution, dead time, decay times of internal X-rays, and system stability, have been tested and measured. The results show that the 1 kg point-contact germanium detector, together with its shielding system and electronics, can run smoothly with good performances. This detector system will be deployed for dark matter search experiments.

The CDEX Collaboration has been established for direct detection of light dark matter particles, using ultra-low energy threshold p-type point-contact germanium detectors, in China JinPing underground Laboratory (CJPL). The first 1 kg point-contact germanium detector with a sub-keV energy threshold has been tested in a passive shielding system located in CJPL. The outputs from both the point-contact p+ electrode and the outside n+ electrode make it possible to scan the lower energy range of less than 1 keV and at the same time to detect the higher energy range up to 3 MeV. The outputs from both p+ and n+ electrode may also provide a more powerful method for signal discrimination for dark matter experiment. Some key parameters, including energy resolution, dead time, decay times of internal X-rays, and system stability, have been tested and measured. The results show that the 1 kg point-contact germanium detector, together with its shielding system and electronics, can run smoothly with good performances. This detector system will be deployed for dark matter search experiments.

Weakly Interacting Massive Particles (WIMPs) are the candidates of dark matter in our universe. Up to now any direct interaction of WIMP with nuclei has not been observed yet. The exclusion limits of the spin-independent cross section of WIMP-nucleon which have been experimentally obtained is about 10^{-7}pb at high mass region and only 10^{-5}pb} at low mass region. China Jin-Ping underground laboratory CJPL is the deepest underground lab in the world and provides a very promising environment for direct observation of dark matter. The China Dark Matter Experiment (CDEX) experiment is going to directly detect the WIMP flux with high sensitivity in the low mass region. Both CJPL and CDEX have achieved a remarkable progress in recent two years. The CDEX employs a point-contact germanium semi-conductor detector PCGe whose detection threshold is less than 300 eV. We report the measurement results of Muon flux, monitoring of radioactivity and Radon concentration carried out in CJPL, as well describe the structure and performance of the 1 kg PCGe detector CDEX-1 and 10kg detector array CDEX-10 including the detectors, electronics, shielding and cooling systems. Finally we discuss the physics goals of the CDEX-1, CDEX-10 and the future CDEX-1T detectors.

Weakly Interacting Massive Particles (WIMPs) are the candidates of dark matter in our universe. Up to now any direct interaction of WIMP with nuclei has not been observed yet. The exclusion limits of the spin-independent cross section of WIMP-nucleon which have been experimentally obtained is about 10^{-7}pb at high mass region and only 10^{-5}pb} at low mass region. China Jin-Ping underground laboratory CJPL is the deepest underground lab in the world and provides a very promising environment for direct observation of dark matter. The China Dark Matter Experiment (CDEX) experiment is going to directly detect the WIMP flux with high sensitivity in the low mass region. Both CJPL and CDEX have achieved a remarkable progress in recent two years. The CDEX employs a point-contact germanium semi-conductor detector PCGe whose detection threshold is less than 300 eV. We report the measurement results of Muon flux, monitoring of radioactivity and Radon concentration carried out in CJPL, as well describe the structure and performance of the 1 kg PCGe detector CDEX-1 and 10kg detector array CDEX-10 including the detectors, electronics, shielding and cooling systems. Finally we discuss the physics goals of the CDEX-1, CDEX-10 and the future CDEX-1T detectors.

We report on a pilot study on identifying metal-poor stars pre-enriched by Pair-Instability Supernovae (PISNe). Very massive, first generation (Population III) stars (140M\odot \leq M \leq 260M\odot) end their lives as PISNe, which have been predicted by theories, but no relics of PISNe have been observed yet. Among the distinct characteristics of the yields of PISN, as predicted by theoretical calculations, are a strong odd-even effect, and a strong overabundance of Ca with respect to iron and the Solar ratio. We use the latter characteristic to identify metal-poor stars in the Galactic halo that have been pre-enriched by PISN, by comparing metallicites derived from strong, co-added Fe lines detected in low-resolution (i.e., R \sim 2000) spectra of the Sloan Digital Sky Survey (SDSS), with metallicities determined by the SDSS Stellar Parameters Pipeline (SSPP). The latter are based on the strength of the Ca II K line and assumptions on the Ca/Fe abundance ratio. Stars are selected as candidates if their metallicity derived from Fe lines is significantly lower than the SSPP metallicities. In a sample of 12,300 stars for which SDSS spectroscopy is available, we have identified 18 candidate stars. Higher resolution and signal-to-noise ratio spectra of these candidates are being obtained with the Very Large Telescope of the European Southern Observatory and the XSHOOTER spectrograph, to determine their abundance patterns, and to verify our selection method. We plan to apply our method to the data base of several million stellar spectra to be acquired with the LAMOST telescope in the next five years.

We report thorium abundances for 77 metal-poor stars in the metallicity range of -3.5<[Fe/H]}<-1.0, based on "snapshot" spectra obtained with VLT-UT2/UVES during the HERES Survey. We were able to determine the thorium abundances with better than 1-sigma confidence for 17 stars, while for 60 stars we derived upper limits. For five stars common with previous studies, our results were in good agreement with the literature results. The thorium abundances span a wide range of about 4.0 dex, and scatter exists in the distribution of log(Th/Eu) ratios for lower metallicity stars, supporting previous studies suggesting the r-process is not universal. We derived ages from the log(Th/Eu) ratios for 12 stars, resulting in large scattered ages, and two stars with significant enhancement of Th relative to Eu are found, indicating the "actinide boost" does not seem to be a rare phenomenon and thus highlighting the risk in using log(Th/Eu) to derive stellar ages.

Hypervelocity stars are believed to be ejected out from the Galactic center through dynamical interactions of (binary) stars with the central massive black hole(s). In this letter, we report 13 metal-poor F-type hypervelocity star candidates selected from 370,000 stars of the data release 7 of the Sloan Digital Sky Survey. With a detailed analysis of the kinematics of these stars, we find that seven of them were likely ejected from the Galactic center (GC) or the Galactic disk, four neither originated from the GC nor the Galactic disk, and the other two were possibly ejected from either the Galactic disk or other regions. Those candidates which unlikely originated from the GC or the Galactic disk, may be explained by other mechanisms, like the tidal disruption of the Milky Way's dwarf galaxies in the Galactic potential, or the gravitational interactions with a massive black hole at the center of M31 or M32.

We present an observational study towards the young high-mass star forming region G23.44-0.18 using the Submillimeter Array. Two massive, radio-quiet dusty cores MM1 and MM2 are observed in 1.3 mm continuum emission and dense molecular gas tracers including thermal CH$_3$OH, CH$_3$CN, HNCO, SO, and OCS lines. The $^{12}$CO (2--1) line reveals a strong bipolar outflow originated from MM2. The outflow consists of a low-velocity component with wide-angle quasi-parabolic shape and a more compact and collimated high-velocity component. The overall geometry resembles the outflow system observed in the low-mass protostar which has a jet-driven fast flow and entrained gas shell. The outflow has a dynamical age of $6\times10^3$ years and a mass ejection rate $\sim10^{-3} M_{\odot}$ year$^{-1}$. A prominent shock emission in the outflow is observed in SO and OCS, and also detected in CH$_3$OH and HNCO. We investigated the chemistry of MM1, MM2 and the shocked region. The dense core MM2 have molecular abundances of 3 to 4 times higher than those in MM1. The abundance excess, we suggest, can be a net effect of the stellar evolution and embedded shocks in MM2 that calls for further inspection.

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.

We study Hawking-like radiation in a Friedmann-Robertson-Walker (FRW) universe using the quasi-classical WKB/tunneling method which pictures this process as a "tunneling" of particles from behind the apparent horizon. The correct temperature of the Hawking-like radiation from the FRW spacetime is obtained using a canonical invariant tunneling amplitude. In contrast to the usual quantum mechanical WKB/tunneling problem where the tunneling amplitude has only a spatial contribution, we find that the tunneling amplitude for FRW spacetime (i.e. the imaginary part of the action) has both spatial and temporal contributions. In addition we study back reaction and energy conservation of the radiated particles and find that the tunneling probability and change in entropy, ${\cal S}$ are related by the relationship: $\Gamma\propto\exp[-\Delta {\cal S}]$ which differs from the standard result $\Gamma\propto\exp[\Delta {\cal S}]$. By regarding the whole FRW universe as an isolated adiabatic system the change in the total entropy is zero. Then splitting the entropy between interior and exterior parts of the horizon ($\Delta {\cal S}_{total}=\Delta {\cal S}_{int} + \Delta {\cal S}_{ext}=0$), we can explain the origin of the minus sign difference with the usual result: our $\Delta {\cal S}$ is for the interior region while the standard result from black hole physics is for the exterior region.

By making use of Duan-Ge's decomposition theory of gauge potential and the Duan's topological current theory proposed by Prof. Duan Yi-Shi, we study a two component superfluid Bose condensed system, which is supposed being realized in the interior of neutron stars in the form of a coexistent neutron superfluid and protonic superconductor. We propose that this system possesses vortex lines. The topological charge of the vortex lines are characterized by the Hopf indices and the Brower degrees of $\phi$-mapping.

Recently, a Hamilton-Jacobi method beyond semiclassical approximation in black hole physics was developed by \emph{Banerjee} and \emph{Majhi}\cite{beyond0}. In this paper, we generalize their analysis of black holes to the case of Friedmann-Robertson-Walker (FRW) universe. It is shown that all the higher order quantum corrections in the single particle action are proportional to the usual semiclassical contribution. The corrections to the Hawking-like temperature and entropy of apparent horizon for FRW universe are also obtained. In the corrected entropy, the area law involves logarithmic area correction together with the standard inverse power of area term.

We explore the noncommutative effect on single field inflation and compare with WMAP five-year data. First, we calculate the noncommutative effect from the potential and dynamical terms, and construct the general form of modified power spectrum. Second, we consider the leading order modification of slow-roll, DBI and K-inflation and unite the modification, which means the modification is nearly model independent at this level. Finally, comparing with the WMAP5 data, we find that the modified can be well realized as the origin of the relative large spectral index and the quite small running.

The influence of the generalized uncertainty principle (GUP) and extended uncertainty principle (EUP) on the thermodynamics of the Friedmann-Robertson-Walker (FRW) universe has been investigated. It is shown that the entropy of the apparent horizon of the FRW universe gets a correction if one considers the effect of the GUP or EUP. Moreover, starting with the modified entropy and applying the first law of thermodynamics, $dE=TdS$, to the apparent horizon of the FRW universe, we obtain the modified Friedmann equations. The influence of the GUP or EUP on the thermodynamics of the FRW universe provides a deep insight into the understanding of the quantum gravity or large length scale corrections to the dynamics of the FRW universe.

The vacuum inflation with a boundary condition specified at a short-distance scale in a generally primeval non-flat Universe, and implications of the correspondingly modified primordial power spectrum via the WMAP data are investigated. We obtain a general form of the modified primordial power spectrum including the effects of both possibly new physics scale and the spatial curvature. The modulation of the primordial power spectrum due to new physics is of the order $H/\Lambda$, where $H$ is the Hubble parameter during inflation and $\Lambda$ is the new physics scale that is key to inflation, while the modulation from the curvature effect is of the order $K/k^2$, where $K$ is the spatial curvature before inflation and $k$ is the comoving wave number. We also find that a closed Universe before inflation is favored by WMAP five-year data.

Using the data taken from Tibet II High Density (HD) Array (1997 February-1999 September) and Tibet-III array (1999 November-2005 November), our previous northern sky survey for TeV $\gamma-$ray point sources has now been updated by a factor of 2.8 improved statistics. From $0.0^{\circ}$ to $60.0^{\circ}$ in declination (Dec) range, no new TeV $\gamma-$ray point sources with sufficiently high significance were identified while the well-known Crab Nebula and Mrk421 remain to be the brightest TeV $\gamma-$ray sources within the field of view of the Tibet air shower array. Based on the currently available data and at the 90% confidence level (C.L.), the flux upper limits for different power law index assumption are re-derived, which are approximately improved by 1.7 times as compared with our previous reported limits.

One of the simplest way to extend 4D cosmological models is to add another spatial dimension to make them 5D. In particular, it has been shown that the simplest of such 5D models, i.e., one in which the right hand side of the Einstein equation is empty, induces a 4D nonempty Universe. Accordingly, the origin of matter in the 4D real Universe might be mathematically attributed to the existence of one (fictitious) extra spatial dimension. Here we consider the case of an empty 5D Universe possessing a cosmological constant $\Lambda$ and obtain exact solutions for both positive and negative values of $\Lambda$. It is seen that such a model can naturally reduce to a power law $\Lambda$CDM model for the real Universe. Further, it can be seen that, the arbitrary constants and functions appearing in this model are endowed with definite physical meanings.

The tachyon field in cosmology is studied by applying the generating function method to obtain exact solutions. The equation of state parameter of the tachyon field is $w=-1+\epsilon\dot{\phi^2}$, which can be expressed as a function in terms of the redshift $z$. Based on these solutions, we propose some tachyon-inspired dark energy models to explore the properties of the corresponding cosmological evolution. The explicit relations between Hubble parameter and redshift enable us to test the models with SNe Ia data sets easily. In the current work we employ the SNe Ia data with the parameter $\mathcal{A}$ measured from the SDSS and the shift parameter $\mathcal{R}$ from WMAP observations to constrain the parameters in our models.

The amplitude of the Compton-Getting (CG) anisotropy contains the power-law index of the cosmic-ray energy spectrum. Based on this relation and using the Tibet air-shower array data, we measure the cosmic-ray spectral index to be $-3.03 \pm 0.55_{stat} \pm < 0.62_{syst}$ between 6 TeV and 40 TeV, consistent with $-$2.7 from direct energy spectrum measurements. Potentially, this CG anisotropy analysis can be utilized to confirm the astrophysical origin of the ``knee'' against models for non-standard hadronic interactions in the atmosphere.

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We have a plan to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV energy region by discriminating them from cosmic-ray hadrons. We have deduced the attainable sensitivity of the muon detector array using our Monte Carlo simulation. We report here on the detailed procedure of our Monte Carlo simulation.

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We are planning to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV (10 - 1000 TeV) energy region by discriminating them from cosmic-ray hadrons. We report on the possibility of detection of gamma rays in the 100 TeV energy region in our field of view, based on the improved sensitivity of our air shower array deduced from the full Monte Carlo simulation.

We generalize the $\Lambda$CDM model by introducing a unified EOS to describe the Universe contents modeled as dark viscous fluid, motivated by the fact that a single constant equation of state (EOS) $p=-p_0$ ($p_0>0$) reproduces the $\Lambda$CDM model exactly. This EOS describes the perfect fluid term, the dissipative effect, and the cosmological constant in a unique framework and the Friedmann equations can be analytically solved. Especially, we find a relation between the EOS parameter and the renormalizable condition of a scalar field. We develop a completely numerical method to perform a $\chi^2$ minimization to constrain the parameters in a cosmological model directly from the Friedmann equations, and employ the SNe data with the parameter $\mathcal{A}$ measured from the SDSS data to constrain our model. The result indicates that the dissipative effect is rather small in the late-time Universe.

With considerations of the recently released WMAP year three and supernova legacy survey (SNLS) data set analysis that favors models similar to the $% \Lambda CDM$ model by possibly mild fluctuations around the vacuum energy or the cosmological constant, we extend the original Chaplygin Gas model (ECG) via modifying the Chaplygin Gas equation of state by two parameters to describe an entangled mixture state from an available matter and the rest component (which can take the cosmological constant or dark energy as in the current cosmic stage, or `curvature-like' term, or radiation component in the early epoch, as various phases) coexistence. At low redshifts, the connection of the ECG model and the Born-infeld field is set up. As paradigms, we use the data coming from the recently released SNLS for the first year and also the famous 157 type Ia supernova (Ia SNe) gold dataset to constrain the model parameters. The restricted results demonstrate clearly how large the entangled degree or the ratio between the energy density parameters of the two entangled phases being. The fact that the ECG models are consistent with the observations of Ia SNe is obtained through the redshift-luminosity distance diagram, hence the ECG can be regarded possible candidates for mimicking the current speed-up expansion of our universe.

A generalized dynamical equation for the scale factor of the universe is proposed to describe the cosmological evolution, of which the $\Lambda$CDM model is a special case. It also provides a general example to show the equivalence of the modified equation of state (EOS) and a scalar field model. In the mathematical aspect, the EOS, the scalar field potential $V(\phi)$, and the scale factor $a(t)$ all have possessed analytical solutions. Such features are due to a simple form invariance of the equation inherited which determines the Hubble parameter. From the physical point of view, this dynamical equation can be regarded as the $\Lambda$CDM model with bulk viscosity, an existence content in the universe. We employ the SNe Ia data with the parameter $\mathcal{A}$ measured from the SDSS data and the shift parameter $\mathcal{R}$ measured from WMAP data to constrain the parameters in our model. The result is that the contribution of the bulk viscosity, accumulated as an effective dark energy responsible for the current cosmic accelerating expansion, is made approximately ten percent to that of the cosmological constant.

A generally parameterized equation of state (EOS) is investigated in the cosmological evolution with bulk viscosity media modelled as dark fluid, which can be regarded as a unification of dark energy and dark matter. Compared with the case of the perfect fluid, this EOS has possessed four additional parameters, which can be interpreted as the case of the non-perfect fluid with time-dependent viscosity or the model with variable cosmological constant. From this general EOS, a completely integrable dynamical equation to the scale factor is obtained with its solution explicitly given out. (i) In this parameterized model of cosmology, for a special choice of the parameters we can explain the late-time accelerating expansion universe in a new view. The early inflation, the median (relatively late time) deceleration, and the recently cosmic acceleration may be unified in a single equation. (ii) A generalized relation of the Hubble parameter scaling with the redshift is obtained for some cosmology interests. (iii) By using the SNe Ia data to fit the effective viscosity model we show that the case of matter described by $p=0$ plus with effective viscosity contributions can fit the observational gold data in an acceptable level

In this paper the accelerating expansion of our universe at the late cosmic evolution time in a generally modified (extended) \emph{Chaplygin gas} (Dark Fluid) model is detailed, which is characterized by two parameters ($m$, $\alpha$). Different choices to the parameters $m$ and $\alpha$ divide this model into two main kinds of situation by different properties with cosmological interests. With proper choices of parameters, we find that this extended model can realize the phantom divide $w=-1$ (Equation Of State parameter) crossing phenomenon with interesting ranges of the scale factor value $a(t)$ corresponding to $w=-1$ and present value of state parameter $w$. Additionally, through Taylor series expansion of the function $A(a)$ in the extended Chaplygin gas model, a specific Equation of State is gained. Under the framework of Friedman-Robertson-Walker cosmic model, it can successfully explain the accelerating expansion of our universe. However, the value of $w$ in this case is large than -1, that is, indicating it like a quintessence fluid and no $w=-1$ crossing occurs.

The universe media is considered as a non-perfect fluid with bulk viscosity and described by a more general equation of state. We assume the bulk viscosity is a linear combination of the two terms: one is constant, and the other is proportional to the scalar expansion $\theta=3\dot{a}/a$. The equation of state is described as $p=(\gamma-1)\rho+p_0$, where $p_0$ is a parameter. This model can be used to explain the dark energy dominated universe. Different choices of the parameters may lead to three kinds of fates of the cosmological evolution: no future singularity, big rip, or Type III singularity of Ref. [S. Nojiri, S.D. Odintsov, and S. Tsujikawa, Phys. Rev. D \textbf{71}, 063004 (2005)].