Recent observations of VER J2227+608 reveal a broken power $\gamma$-ray spectrum with the spectral index increasing from $\sim 1.8$ in the GeV energy range to $\sim 2.3$ in the TeV range. Such a spectral break can be attributed to radiative energy loss of energetic electrons in the leptonic scenario for the $\gamma$-ray emission, which, in combination with characteristic age of the nearby pulsar, can be used to constrain magnetic field in the emission region. We show that the radio and X-ray observations can also be explained in such a scenario. In the hadronic scenario, the spectral break can be attributed to diffusion of energetic ions in a turbulent medium and detailed spectral measurement can be used to constrain the diffusion coefficient. These two models, however, predict drastically different spectra above 100 TeV, which will be uncovered with future high-resolution observations, such as LHAASO and/or CTA.

We installed two sets of Astronomical Site Monitoring System(ASMS) at Lijiang Observatory(GMG), for the running of the 2.4-meter Lijiang optical telescope(LJT) and the 1.6-meter Multi-channel Photometric Survey Telescope (Mephisto). The Mephistro is under construction. ASMS has been running on robotic mode since 2017. The core instruments: Cloud Sensor, All-Sky Camera and Autonomous-DIMM that are developed by our group, together with the commercial Meteorological Station and Sky Quality Meter, are combined into the astronomical optical site monitoring system. The new Cloud Sensor's Cloud-Clear Relationship is presented for the first time, which is used to calculate the All-Sky cloud cover. We designed the Autonomous-DIMM located on a tower, with the same height as LJT. The seeing data have been observed for a full year. ASMS's data for the year 2019 are also analysed in detail, which are valuable to observers.

Observations by HAWC and Milagro have detected spatially extended TeV sources surrounding middle-aged ($t\sim100-400 \,{\rm kyr}$) pulsars like Geminga and PSR B0656+14, which have been named "TeV Halos", representing very extended TeV pulsar wind nebulae (PWNe) powered by relatively old pulsars. A few more HAWC-detected sources have been suggested to be TeV halo candidates. In this paper, we search for possible GeV counterparts of three TeV halo candidates with Fermi Large Area Telescopes. We detect a new spatially extended GeV source in the vicinity of the TeV halo candidate 2HWC J1912+099, which is also detected by HESS (HESS J1912+101). We find that the size of the GeV source is significantly larger than that of the TeV emission measured by HESS, and a spatial template characteristic of a PWN can fit the GeV data. We suggest that the GeV source is an extended PWN powered by the central middle-aged pulsar PSR J1913+1011. This discovery favors the TeV halo scenario for the TeV source 2HWC J1912+099 (HESS J1912+101), although the possible shell-like morphology measured by HESS challenges this interpretation. Alternatively, the TeV emission could be dominated by a supernova remnant (SNR) via the hadronic process. Future multi-wavelength observations of the source and more precise measurements of the spatial profile of the TeV emission will be useful to distinguish between the two scenarios.

In order to study short timescale optical variability of $\gamma$-ray blazar S5 0716+714, quasi-simultaneous spectroscopic and multi-band photometric observations were performed from 2018 November to 2019 March with the 2.4 m optical telescope located at Lijiang Observatory of Yunnan Observatories. The observed spectra are well fitted with a power-law $F_{\lambda}=A\lambda ^{-\alpha}$ (spectral index $\alpha >0$). Correlations found between $\dot{\alpha}$, $\dot{A}$, $\dot{A}/A$, $\dot{F_{\rm{\lambda}}}$, and $\dot{F_{\rm{\lambda}}}/F_{\rm{\lambda}}$ are consistent with the trend of bluer-when-brighter (BWB). \textbf{The same case is for colors, magnitudes, color variation rates, and magnitude variation rates of photometric observations.} The variations of $\alpha$ lead those of $F_{\rm{\lambda}}$. Also, the color variations lead the magnitude variations. The observational data are mostly distributed in the I(+,+) and III(-,-) quadrants of coordinate system. Both of spectroscopic and photometric observations show BWB behaviors in S5 0716+714. The observed BWB may be explained by the shock-jet model, and its appearance may depend on the relative position of the observational frequency ranges with respect to the synchrotron peak frequencies, e.g., at the left of the peak frequencies. \textbf{Fractional variability amplitudes are $F_{\rm{var}}\sim 40\%$ for both of spectroscopic and photometric observations. Variations of $\alpha$ indicate variations of relativistic electron distribution producing the optical spectra. }

The spectra fits to a sample of 34 supernova remnants (Zeng et al., 2019) are updated. $\gamma$-ray spectra of 20 supernova remnants (SNRs) with a soft TeV spectrum are further analyzed. We found that 17 of them can be fitted in the hadronic scenario with a single power-law ion distribution with an index of $\sim$ 2.6, which is significantly softer than the ion distribution inferred from $\gamma$-ray observations of star-burst galaxies. If Galactic cosmic rays are mostly produced by SNRs, this result suggests that SNRs in star-burst galaxies may never reach the phase with a soft $\gamma$-ray spectra or escape of high-energy particles from SNRs before they reach this phase with a soft $\gamma$-ray spectrum dominates the contribution of SNRs to Galactic cosmic rays.

We developed a spectroscopic monitoring project to investigate the kinematics of the broad-line region (BLR) in active galactic nuclei (AGN) with ultra-fast outflows (UFOs). Mrk~79 is a radio-quiet AGN with UFOs and warm absorbers, had been monitored by three reverberation mapping (RM) campaigns, but its BLR kinematics is not understood yet. In this paper, we report the results from a new RM-campaign of Mrk~79, which was undertaken by Lijiang 2.4-m telescope. Mrk~79 is seeming to come out the faint state, the mean flux approximates a magnitude fainter than historical record. We successfully measured the lags of the broad emission lines including H$\beta~\lambda4861$, H$\gamma~\lambda4340$, He II $\lambda4686$ and He I $\lambda5876$ with respect to the varying AGN continuum. Based on the broad H$\beta~\lambda4861$ line, we measured black hole (BH) mass of $M_{\bullet}=5.13^{+1.57}_{-1.55}\times10^{7}M_{\odot}$, estimated accretion rates of ${\dot{M}_{\bullet}}=(0.05\pm0.02)~L_{\rm Edd}~c^{-2}$, indicating that Mrk~79 is a sub-Eddington accretor. We found that Mrk~79 deviates from the canonical Radius$-$Luminosity relationship. The marginal blueshift of the broad He II $\lambda4686$ line detected from rms spectrum indicates outflow of high-ionization gas. The velocity-resolved lag profiles of the broad H$\gamma~\lambda4340$, H$\beta~\lambda4861$, and He I $\lambda5876$ lines show similar signatures that the largest lag occurs in the red wing of the lines then the lag decreases to both sides. These signatures should suggest that the BLR of Keplerian motion probably exists the outflow gas motion. All findings including UFOs, warm absorbers, and the kinematics of high- and low-ionization BLR, may provide an indirect evidence that the BLR of Mrk~79 probably originates from disk wind.

Linking a coronagraph instrument to a spectrograph via a single mode optical fiber is a pathway towards detailed characterization of exoplanet atmospheres with current and future ground- and space-based telescopes. However, given the extreme brightness ratio and small angular separation between planets and their host stars, the planet signal-to-noise ratio will likely be limited by the unwanted coupling of starlight into the fiber. To address this issue, we utilize a wavefront control loop and a deformable mirror to systematically reject starlight from the fiber by measuring what is transmitted through the fiber. The wavefront control algorithm is based on the formalism of electric field conjugation (EFC), which in our case accounts for the spatial mode selectivity of the fiber. This is achieved by using a control output that is the overlap integral of the electric field with the fundamental mode of a single mode fiber. This quantity can be estimated by pair-wise image plane probes injected using a deformable mirror. We present simulation and laboratory results that demonstrate our approach offers a significant improvement in starlight suppression through the fiber relative to a conventional EFC controller. With our experimental setup, which provides an initial normalized intensity of $3\times10^{-4}$ in the fiber at an angular separation of $4\lambda/D$, we obtain a final normalized intensity of $3\times 10^{-6}$ in monochromatic light at $\lambda=635$~nm through the fiber (100x suppression factor) and $2\times 10^{-5}$ in $\Delta\lambda/\lambda=8\%$ broadband light about $\lambda=625$~nm (10x suppression factor). The fiber-based approach improves the sensitivity of spectral measurements at high contrast and may serve as an integral part of future space-based exoplanet imaging missions as well as ground-based instruments.

We report the detection of GeV $\gamma$-ray emission from the very-high-energy (VHE) $\gamma$-ray source VER J2227+608 associated with the "tail" region of SNR G106.3+2.7. The GeV $\gamma$-ray emission is extended and spatially coincident with molecular clouds traced by CO emission. The broadband GeV to TeV emission of VER J2227+608 can be well fitted by a single power-law function with an index of 1.90$\pm$0.04, without obvious indication of spectral cutoff toward high energies. The pure leptonic model for the $\gamma$-ray emission can be marginally ruled out by the X-ray and TeV data. In the hadronic model, the low energy content of CRs and the hard $\gamma$-ray spectrum, in combination with the center-bright source structure, suggest that VER J2227+608 may be powered by the PWN instead of shocks of the SNR. And the cutoff energy of the proton distribution needs to be higher than $\sim$ 400 TeV, which makes it an attractive PeVatron candidate. Future observations by the upcoming Large High Altitude Air Shower Observatory (LHAASO) and the Cherenkov Telescope Array in the north (CTA-North) could distinguish these models and constrain the maximum energy of cosmic rays in supernova remnants.

Circinus galaxy is a nearby composite starburst/AGN system. In this work we re-analyze the GeV emission from Circinus with 10 years of {\it Fermi}-LAT Pass 8 data. In the energy range of 1-500 GeV, the spectrum can be well fitted by a power-law model with a photon index of $\Gamma$ = $2.20\pm0.14$, and its photon flux is $(5.90\pm1.04) \times 10^{-10}$ photons cm$^{-2}$ s$^{-1}$. Our 0.1-500 GeV flux is several times lower than that reported in the previous literature, which is roughly in compliance with the empirical relation for star-forming and local group galaxies and might be reproduced by the interaction between cosmic rays and the interstellar medium. The ratio between the $\gamma$-ray luminosity and the total infrared luminosity is near the proton calorimetric limit, indicating that Circinus may be a proton calorimeter. However, marginal evidence for variability of the $\gamma$-ray emission is found in the timing analysis, which may indicate the activity of AGN jet. More {\it Fermi}-LAT data and future observation of CTA are required to fully reveal the origin of its $\gamma$-ray emission.

We present the photometric observations of blazars S5 0716+714 and 3C 273 with high temporal resolution (30--60s) in the $I$ or $R$ bands. The observations were performed with a 1.02 m optical telescope from 2007 March 07 to 2012 May 16. The $F$-test, one-way analysis of variance (ANOVA) test, and $z$-transformed discrete correlation function (ZDCF) cross-correlation zero lag test are used to search for intra-day variability (IDV). Four and five reliable IDVs survive three tests for S5 0716+714 and 3C 273, respectively. IDVs are found for S5 0716+714 and 3C 273. A flare on 2008 May 08 has $\Delta I \approx$ 0.06$\pm$0.01 mag in a duration of 0.54 hr for S5 0716+714. A flare on 2011 May 10 shows $\Delta R \approx$ 0.05$\pm$0.01 mag in a duration of 0.40 hr for 3C 273. Sharp dips appear on 2011 May 9 for 3C 273, and show $\Delta R \approx$ 0.05$\pm$0.01 mag. Under the assumptions that the IDV is tightly connected to black hole mass, $M_{\bullet}$, and that the flare durations are representative of the minimum characteristic timescales, we can estimate upper bounds to $M_{\bullet}$. In the case of the Kerr black holes, $M_{\bullet} \la 10^{8.91} M_{\odot}$ and $M_{\bullet} \la 10^{9.02} M_{\odot}$ are given for S5 0716+714 and 3C 273, respectively. These mass measurements are consistent with those measurements reported in the literatures. Also, we discuss the origins of optical variations found in this work.

The DArk Matter Particle Explorer (DAMPE), a satellite-based cosmic ray and gamma-ray detector, was launched on December 17, 2015, and began its on-orbit operation on December 24, 2015. In this work we document the on-orbit calibration procedures used by DAMPE and report the calibration results of the Plastic Scintillator strip Detector (PSD), the Silicon-Tungsten tracKer-converter (STK), the BGO imaging calorimeter (BGO), and the Neutron Detector (NUD). The results are obtained using Galactic cosmic rays, bright known GeV gamma-ray sources, and charge injection into the front-end electronics of each sub-detector. The determination of the boundary of the South Atlantic Anomaly (SAA), the measurement of the live time, and the alignments of the detectors are also introduced. The calibration results demonstrate the stability of the detectors in almost two years of the on-orbit operation.

Lijiang 2.4-meter Telescope(LJT), the largest common-purpose optical telescope in China, has been applied to the world-wide astronomers since 2008. It is located at Gaomeigu site, Lijiang Observatory(LJO), the southwest of China. The site has very good observational conditions. Since 10-year operation, several instruments have been equipped on the LJT. Astronomers can perform both photometric and spectral observations. The main scientific goals of LJT include photometric and spectral evolution of supernova, reverberation mapping of active galactic nucleus, physical properties of binary star and near-earth object(comet and asteroid), identification of exoplanet, and all kinds of transients. Until now, the masses of 41 high accretion rate black holes have been measured, and more than 168 supernova have been identified by the LJT. More than 190 papers related to the LJT have been published. In this paper, the general observation condition of the Gaomeigu site is introduced at first. Then, the LJT structure is described in detail, including the optical, mechanical, motion and control system. The specification of all the instruments, and some detailed parameters of the YFOSC is also presented. Finally, some important scientific results and future expectations are summarized.

We investigate the prompt emission and the afterglow properties of short duration gamma-ray burst (sGRB) 130603B and another eight sGRB events during 2012-2015, observed by several multi-wavelength facilities including the GTC 10.4m telescope. Prompt emission high energy data of the events were obtained by INTEGRAL/SPI/ACS, Swift/BAT and Fermi/GBM satellites. The prompt emission data by INTEGRAL in the energy range of 0.1-10 MeV for sGRB 130603B, sGRB 140606A, sGRB 140930B, sGRB 141212A and sGRB 151228A do not show any signature of the extended emission or precursor activity and their spectral and temporal properties are similar to those seen in case of other short bursts. For sGRB130603B, our new afterglow photometric data constraints the pre jet-break temporal decay due to denser temporal coverage. For sGRB 130603B, the afterglow light curve, containing both our new as well as previously published photometric data is broadly consistent with the ISM afterglow model. Modeling of the host galaxies of sGRB 130603B and sGRB 141212A using the LePHARE software supports a scenario in which the environment of the burst is undergoing moderate star formation activity. From the inclusion of our late-time data for 8 other sGRBs we are able to; place tight constraints on the non-detection of the afterglow, host galaxy or any underlying kilonova emission. Our late-time afterglow observations of the sGRB 170817A/GW170817 are also discussed and compared with the sub-set of sGRBs.

Li et al. (2018) (hereafter L18) detected two apparently distinct populations of blue straggler stars (BSSs) in the young globular cluster NGC 2173, a similar feature as observed in numerous Galactic globular clusters (GCs). Recently, Dalessandro et al. (2018) (D18) compared the observed color--magnitude diagrams (CMDs) of both the cluster and a nearby reference field (although observed with the ACS/WFC instrument). They conclude that the bifurcated pattern of BSSs in NGC 2173 observed by L18 is a field contamination artefact. In this note, we explore the central concentration properties of the removed `field stars' identified by D18. Our purpose is to examine if these `field stars' are spatially homogeneously distributed. Employing a Monte Carlo-based approach, we have carefully studied the probability that any such central concentration may be caused by small number statistics. We find that, in most cases of, the `field stars' removed by D18 exhibit a clear central concentration, which cannot be explained on the basis of small number statistics alone. Therefore, we suggest that D18 may well have overestimated the field contamination level, implying that the bifurcated BSS pattern in NGC 2173 cannot, in fact, be explained by field contamination.

Supernova remnants (SNRs) have been considered as the dominant contributors to Galactic cosmic rays. However, the relation between high-energy particles trapped in SNRs and cosmic rays observed at the Earth remains obscure. In this paper, we fit the spectral energy distributions of 35 SNRs with a simple one-zone emission model and analyze correlations of model parameters to uncover the evolution of high-energy particle distribution in SNRs. We find that (1) the particle distribution in general can be described by a broken power-law function with a high-energy cutoff for all SNRs; (2) the low-energy spectrum becomes harder and the break energy decreases with aging of SNRs, (3) for most middle-age SNRs, the energy loss timescale of electrons at the high-energy cutoff is approximately equal to the age of the corresponding remnant implying quenching of very high-energy electron acceleration; for young SNRs, this energy loss timescale is shorter than the age of SNRs implying continuous electon acceleration at the cutoff energy; and for a few old age SNRs, the energy loss timescale is longer than the corresponding age which may suggest escaping of higher energy particles from SNRs. Finally, we comment on the implications of these results on the SNR origin of Galactic cosmic rays.

HESS J1640-465 is an extended TeV $\gamma$-ray source and its $\gamma$-ray emission whether from the shell of a supernova remnant (SNR) or a pulsar wind nebula (PWN) is still under debate. We reanalyze the GeV $\gamma$-ray data in the field of HESS J1640-465 using eight years of Pass 8 data recorded by the Fermi Large Area Telescope. An extended GeV $\gamma$-ray source positionally coincident with HESS J1640-465 is found. Its photon spectrum can be described by a power-law with an index of $1.42\pm0.19$ in the energy range of 10-500 GeV, and smoothly connects with the TeV spectrum of HESS J1640-465. The broadband spectrum of HESS J1640-465 can be well fit by a leptonic model with a broken power-law spectrum of electrons with an exponential cut-off at $\sim$ 300 TeV. The spectral properties of HESS J1640-465 are broadly consistent with the characteristics of other sources identified as PWNe, such as the correlations between high-energy luminosity ratios and the physical parameters of pulsar, including spin-down luminosity $\dot{E}$ and characteristic age $\tau_c$. All these pieces of evidence support that the $\gamma$-ray emission of HESS J1640-465 may originate from the PWN powered by PSR J1640-4631 rather than the shell of the SNR G338.3-0.0.

Recently the IceCube collaboration and 15 other collaborations reported the spatial and temporal coincidence between the neutrino event IceCube-170922A and the radio-TeV activity of the blazar TXS 0506+056. Their further analysis on 9.5 years of IceCube data discovered neutrino flare between September 2014 and March 2015, when TXS 0506+056 is however in "quiescent" state. We analyze the Fermi-LAT data in that direction, and find another bright GeV source PKS 0502+049, which is at an angle of $1.2^{\circ}$ from TXS 0506+056, with strong activties during the neutrino flare. No other bright GeV source was detected in the region of interest. Though PKS 0502+049 is $1.2^\circ$ separated from TXS 0506+056, it locates within the directional reconstruction uncertainties of 7 neutrinos, out of the 13 neutrino events during the neutrino flare. Together with the observed high flux of the $\gamma$-ray flare, it may be unreasonable to fully discard the (partial) contribution of PKS 0502+049 to the neutrino flare. The single source assumption used in the neutrino data analysis might need to be revisited.

We present optical spectroscopic and photometric observations of the nearby type Ic supernova (SN Ic) SN 2014L. This SN was discovered by the Tsinghua-NAOC Transient Survey (TNTS) in the nearby type-Sc spiral galaxy M99 (NGC 4254). Fitting to the early-time light curve indicates that SN 2014L was detected at only a few hours after the shock breakout, and it reached a peak brightness of $M_{\rm V} = -17.73 \pm 0.28$ mag ($L= [2.06\pm 0.50] \times 10^{42}$ erg s$^{-1}$) approximately 13 days later. SN 2014L shows a close resemblance to SN 2007gr in the photometric evolution, while it shows stronger absorption features of intermediate-mass elements (especially \CaII) in the early-time spectra. Based on simple modeling of the observed light curves, we derived the mass of synthesized $^{56}$Ni as $M_{\rm Ni} = 0.075 \pm 0.025 M_{\odot}$, and the mass and total energy of the ejecta as $M_{\rm ej} = 1.00 \pm 0.20 M_\odot$ and $E_{\rm ej}=1.45\pm0.25$ foe, respectively. Given these typical explosion parameters, the early detection, and the extensive observations, we suggest that SN 2014L could be a template sample for the investigation of SNe Ic.

On 2017 September 22 the IceCube neutrino observatory detected a track-like, very-high-energy event (IceCube-170922A) that is spatially associated with TXS 0506+056, a quasar at a redshift of $z=0.3365$. This source is characterized by the increased acitivies in a very wide energy range (from radio to TeV) during these days. To investigate the possible connection of the PeV neutrino emission with the GeV activity of blazars, in this work we select 116 bright sources and analyze their lightcurves and spectra. We focus on the sources displaying GeV activities. Among these blazars, TXS 0506+056 seems to be typical in many aspects but is distinguished by the very strong GeV activties. We suggest to search for neutrino outburst in the historical data of IceCube, as recently done for TXS 0506+056, from the directions of these more energetic and harder blazars with strong GeV activities.

The optical observations of Ic-4 supernova (SN) 2016coi/ASASSN-16fp, from $\sim 2$ to $\sim450$ days after explosion, are presented along with analysis of its physical properties. The SN shows the broad lines associated with SNe Ic-3/4 but with a key difference. The early spectra display a strong absorption feature at $\sim 5400$ \AA\ which is not seen in other SNe~Ic-3/4 at this epoch. This feature has been attributed to He I in the literature. Spectral modelling of the SN in the early photospheric phase suggests the presence of residual He in a C/O dominated shell. However, the behaviour of the He I lines are unusual when compared with He-rich SNe, showing relatively low velocities and weakening rather than strengthening over time. The SN is found to rise to peak $\sim 16$ d after core-collapse reaching a bolometric luminosity of Lp $\sim 3\times10^{42}$ \ergs. Spectral models, including the nebular epoch, show that the SN ejected $2.5-4$ \msun\ of material, with $\sim 1.5$ \msun\ below 5000 \kms, and with a kinetic energy of $(4.5-7)\times10^{51}$ erg. The explosion synthesised $\sim 0.14$ \msun\ of 56Ni. There are significant uncertainties in E(B-V)host and the distance however, which will affect Lp and MNi. SN 2016coi exploded in a host similar to the Large Magellanic Cloud (LMC) and away from star-forming regions. The properties of the SN and the host-galaxy suggest that the progenitor had $M_\mathrm{ZAMS}$ of $23-28$ \msun\ and was stripped almost entirely down to its C/O core at explosion.

Bifurcated patterns of blue straggler stars in their color--magnitude diagrams have atracted significant attention. This type of special (but rare) pattern of two distinct blue straggler sequences is commonly interpreted as evidence of cluster core-collapse-driven stellar collisions as an efficient formation mechanism. Here, we report the detection of a bifurcated blue straggler distribution in a young Large MagellanicCloud cluster, NGC 2173. Because of the cluster's low central stellar number density and its young age, dynamical analysis shows that stellar collisions alone cannot explain the observed blue straggler stars. Therefore, binary evolution is instead the most viable explanation of the origin of these blue straggler stars. However, the reason why binary evolution would render the color--magnitude distribution of blue straggler stars bifurcated remains unclear.

The remnant of supernova explosion is widely believed to be the acceleration site of high-energy cosmic ray particles. The acceleration timescale is, however, typically very long. Here we report the detection of a variable $\gamma$-ray source with the Fermi Large Area Telescope, which is positionally and temporally consistent with a peculiar supernova, iPTF14hls. A quasi-stellar object SDSS J092054.04+504251.5, which is probably a blazar according to the infrared data, is found in the error circle of the $\gamma$-ray source. More data about the $\gamma$-ray source and SDSS J092054.04+504251.5 are needed to confirm their association. On the other hand, if the association between the $\gamma$-ray source and the supernova is confirmed, this would be the first time to detect high-energy $\gamma$-ray emission from a supernova, suggesting very fast particle acceleration by supernova explosions.

We report the detection of GeV $\gamma$-ray emission from supernova remnant HESS J1731-347 using 9 years of {\it Fermi} Large Area Telescope data. We find a slightly extended GeV source in the direction of HESS J1731-347. The spectrum above 1 GeV can be fitted by a power-law with an index of $\Gamma = 1.77\pm0.14$, and the GeV spectrum connects smoothly with the TeV spectrum of HESS J1731-347. Either a hadronic-leptonic or a pure leptonic model can fit the multi-wavelength spectral energy distribution of the source. However, the hard GeV $\gamma$-ray spectrum is more naturally produced in a leptonic (inverse Compton scattering) scenario, under the framework of diffusive shock acceleration. We also searched for the GeV $\gamma$-ray emission from the nearby TeV source HESS J1729-345. No significant GeV $\gamma$-ray emission is found, and upper limits are derived.

High energy cosmic ray electrons plus positrons (CREs), which lose energy quickly during their propagation, provide an ideal probe of Galactic high-energy processes and may enable the observation of phenomena such as dark-matter particle annihilation or decay. The CRE spectrum has been directly measured up to $\sim 2$ TeV in previous balloon- or space-borne experiments, and indirectly up to $\sim 5$ TeV by ground-based Cherenkov $\gamma$-ray telescope arrays. Evidence for a spectral break in the TeV energy range has been provided by indirect measurements of H.E.S.S., although the results were qualified by sizeable systematic uncertainties. Here we report a direct measurement of CREs in the energy range $25~{\rm GeV}-4.6~{\rm TeV}$ by the DArk Matter Particle Explorer (DAMPE) with unprecedentedly high energy resolution and low background. The majority of the spectrum can be properly fitted by a smoothly broken power-law model rather than a single power-law model. The direct detection of a spectral break at $E \sim0.9$ TeV confirms the evidence found by H.E.S.S., clarifies the behavior of the CRE spectrum at energies above 1 TeV and sheds light on the physical origin of the sub-TeV CREs.

The DArk Matter Particle Explorer (DAMPE), one of the four scientific space science missions within the framework of the Strategic Pioneer Program on Space Science of the Chinese Academy of Sciences, is a general purpose high energy cosmic-ray and gamma-ray observatory, which was successfully launched on December 17th, 2015 from the Jiuquan Satellite Launch Center. The DAMPE scientific objectives include the study of galactic cosmic rays up to $\sim 10$ TeV and hundreds of TeV for electrons/gammas and nuclei respectively, and the search for dark matter signatures in their spectra. In this paper we illustrate the layout of the DAMPE instrument, and discuss the results of beam tests and calibrations performed on ground. Finally we present the expected performance in space and give an overview of the mission key scientific goals.

Puppis A is a very famous and extensively studied supernova remnant (SNR) that shows strong evidence of shock-cloud interaction. We re-analyze the GeV $\gamma$-ray emission of it using seven years Pass 8 data recorded by the Fermi Large Area Telescope (Fermi-LAT). The morphology of the $\gamma$-ray emission is more compatible with that of the thermal X-ray and IR emissions than the radio image, which suggests a possible correlation between the gamma-ray emitting region and dense clouds. The $\gamma$-ray spectrum in the energy range of 1-500 GeV shows a break at $7.92\pm1.91$ GeV with the photon indices of $1.81\pm0.08$ below the break and $2.53\pm0.12$ above the break, which can naturally explain the lack of TeV $\gamma$-ray emission from Puppis A. The multi-wavelength observations favor a hadronic origin for the $\gamma$-ray emission.

Very few low-ionization broad absorption line (LoBAL) QSOs have been found at high redshifts to date. One high-redshift LoBAL QSO, J0122+1216, was recently discovered at the Lijiang 2.4-m Telescope with an initial redshift determination of 4.76. Aiming to investigate its physical properties, we carried out follow-up observations in the optical and near-IR spectroscopy. Near-IR spectra from UKIRT and P200 confirms that it is a LoBAL, with a new redshift determination of $4.82\pm0.01$ based on the \mgii~ emission-line. The new \mgii~ redshift determination reveals strong blueshifts and asymmetry of the high-ionization emission lines. We estimated a black hole mass of $\sim 2.3\times 10^9 M_\odot$ and Eddington ratio of $\sim 1.0$ according to the empirical \mgii-based single-epoch relation and bolometric correction factor. It is possible that strong outflows are the result of an extreme quasar environment driven by the high Eddington ratio. A lower limit on the outflowing kinetic power ($>0.9\% L_{Edd}$) was derived from both emission and absorption lines, indicating these outflows play a significant role in the feedback process to regulate the growth of its black hole as well as host galaxy evolution.

We report the detection of a GeV $\gamma$-ray source which is spatially overlapping and thus very likely associated with the unidentified very-high-energy (VHE) $\gamma$-ray source HESS J1427-608 with the Pass 8 data recorded by the Fermi Large Area Telescope. The photon spectrum of this source is best described by a power-law with an index of $1.85\pm0.17$ in the energy range of $3-500$ GeV, and the measured flux connects smoothly with that of HESS J1427-608 at a few hundred GeV. This source shows no significant extension and time variation. The broadband GeV-TeV emission over four decades of energies can be well fitted by a single power-law function with an index of 2.0, without obvious indication of spectral cutoff toward high energies. Such a result implies that HESS J1427-608 may be a PeV particle accelerator. We discuss possible nature of HESS J1427-608 according to the multi-wavelength spectral fittings. Given the relatively large errors, either a leptonic or a hadronic model can explain the multi-wavelength data from radio to VHE $\gamma$-rays. The inferred magnetic field strength is a few $\mu$G, which is smaller than typical values of supernova remnants (SNRs), and is consistent with some pulsar wind nebulae (PWNe). On the other hand, the flat $\gamma$-ray spectrum is slightly different from typical PWNe but similar to that of some known SNRs.

Multi-wavelength observations of mature supernova remnants (SNRs), especially with recent advances in gamma-ray astronomy, make it possible to constrain energy distribution of energetic particles within these remnants. In consideration of the SNR origin of Galactic cosmic rays and physics related to particle acceleration and radiative processes, we use a simple one-zone model to fit the nonthermal emission spectra of three shell-type SNRs located within 2 degrees on the sky: RX J1713.7-3946, CTB 37B, and CTB 37A. Although radio images of these three sources all show a shell (or half-shell) structure, their radio, X-ray, and gamma-ray spectra are quite different, offering an ideal case to explore evolution of energetic particle distribution in SNRs. Our spectral fitting shows that 1) the particle distribution becomes harder with aging of these SNRs, implying a continuous acceleration process, and the particle distributions of CTB 37A and CTB 37B in the GeV range are harder than the hardest distribution that can be produced at a shock via the linear diffusive shock particle acceleration process, so spatial transport may play a role; 2) the energy loss timescale of electrons at the high-energy cutoff due to synchrotron radiation appears to be always a bit (within a factor of a few) shorter than the age of the corresponding remnant, which also requires continuous particle acceleration; 3) double power-law distributions are needed to fit the spectra of CTB 37B and CTB 37A, which may be attributed to shock interaction with molecular clouds.

III Zw 2 is the prototype of radio-intermediate quasars. Although there is the evidence of possessing strong jet, significant $\gamma$-ray emission has not been reported before. In this work, we carry out a detailed analysis of the latest {\it Fermi}-LAT {\it Pass} 8 data. No significant $\gamma$-ray signal has been detected in the time-averaged 7-year {\it Fermi}-LAT data of III Zw 2, we however have identified two distinct $\gamma$-ray flares with isotropic luminosities of $\sim 10^{45}$ erg $\rm s^{-1}$. Multiwavelength data analysis (including also the optical photometric observations from Yunnan Observatories) are presented and the main finding is the simultaneous optical and $\gamma$-ray flares of III Zw 2 appearing in Nov. 2009. Violent $\gamma$-ray variability with doubling timescale of 2.5 hours was detected in another $\gamma$-ray flare in May 2010, for which the 3-hour $\gamma$-ray peak flux is $\sim 250$ times of the average flux in 7 years. Rather similar behaviors are observed in blazars and the blazar model can reasonably reproduce the spectral energy distribution of III Zw 2 in a wide energy range, strongly suggesting that its central engine resembles that of blazars. In view of its core which shares radio similarities with young radio sources together with the weak extended radio lobe emission, we suggest that III Zw 2 harbors a recurrent activity core and serves as a valuable target for investigating the fuelling and triggering of the activity in radio loud AGNs.

We used the newly commissioned 50 cm Binocular Network (50BiN) telescope at Qinghai Station of Purple Mountain Observatory (Chinese Academy of Sciences) to observe the old open cluster NGC 188 in V and R as part of a search for variable objects. Our time-series data span a total of 36 days. Radial velocity and proper-motion selection resulted in a sample of 532 genuine cluster members. Isochrone fitting was applied to the cleaned cluster sequence, yielding a distance modulus of (m - M)0V = 11.35 \pm 0.10 mag and a total foreground reddening of E(V - R) = 0.062 \pm 0.002 mag. Light-curve solutions were obtained for eight W Ursae Majoris eclipsing-binary systems (W UMas) and their orbital parameters were estimated. Using the latter parameters, we estimate a distance to the W UMas which is independent of the host cluster's physical properties. Based on combined fits to six of the W UMas (EP Cep, EQ Cep, ES Cep, V369 Cep, and--for the first time--V370 Cep and V782 Cep), we obtain an average distance modulus of (m - M)0V = 11.31 \pm 0.08 mag, which is comparable with that resulting from our isochrone fits. These six W UMas exhibit an obvious period-luminosity relation. We derive more accurate physical parameters for the W UMa systems and discuss their initial masses and ages. The former show that these W UMa systems have likely undergone angular-momentum evolution within a convective envelope (W-type evolution). The ages of the W UMa systems agree well with the cluster's age.

Extensive and independent observations of Type Ia supernova (SN Ia) SN 2013dy are presented, including a larger set of $UBVRI$ photometry and optical spectra from a few days before the peak brightness to $\sim$ 200 days after explosion, and ultraviolet (UV) photometry spanning from $t \approx -10$ days to $t \approx\ +15$ days referring to the $B$ band maximum. The peak brightness (i.e., $M_{\rm B} = -19.65 \pm 0.40$ mag, $L_{\rm max} = [1.95 \pm 0.55] \times 10^{43}$ erg s$^{-1}$) and the mass of synthesised $^{56}$Ni (i.e., $M$($^{56}$Ni) = 0.90 $\pm$ 0.26 M$_{\odot}$) are calculated, and they conform to the expectation for a SN Ia with a slow decline rate (i.e., $\Delta m_{15}(B)$ = 0.90 $\pm$ 0.03 mag, Phillips 1993). However, the near infrared (NIR) brightness of this SN (i.e., $M_{\rm H} = -17.33 \pm 0.30$ mag) is at least 1.0 mag fainter than usual. Besides, spectroscopy classification reveals that SN 2013dy resides on the border of "core normal" and "shallow silicon" subclasses in the Branch et al. (2009) classification scheme, or on the border of the "normal velocity" SNe Ia and 91T/99aa-like events in the Wang et al. (2009a) system. These suggest that SN 2013dy is a slow-declining SN Ia located on the transitional region of nominal spectroscopic subclasses and might not be a typical normal sample of SNe Ia.

In the third Fermi catalogue (3FGL) there are sixteen gamma-ray globular clusters. After analyzing the recent released Pass 8 data of Fermi Large Area Telescope (LAT), we report the discovery of significant gamma-ray emission from M 15 and NGC 6397, confirm that NGC 5904 is a gamma-ray emitter and find evidence of gamma-ray emission from NGC 6218 and NGC 6139. Inter- estingly, in the globular clusters M 15, NGC 6397 and NGC 5904, millisecond pulsars (MSPs) have been found in radio or X-rays, which are strongly in support of the MSP origin of the gamma-ray emission. However, due to the relative low luminosity of the gamma-ray emission we do not find any evidence for the gamma-ray pulsation or flux variability of these sources.

Based on Kepler data, we present the results of a search for white-light flares on 1049 close binaries. We identify 234 flare binaries, on which 6818 flares are detected. We compare the flare-binary fraction in different binary morphologies ("detachedness"). The result shows that the fractions in over-contact and ellipsoidal binaries are approximately 10-20 percent lower than those in detached and semi-detached systems. We calculate the binary flares activity level (AL) of all the flare binaries, and discuss its variations along the orbital period (P_orb) and rotation period (P_rot, calculated for only detached binaries). We find that AL increases with decreasing P_orb or P_rot up to the critical values at P_orb near 3 days or P_rot near 1.5 days, thereafter, the AL starts decreasing no matter how fast the stars rotate. We examine the flaring rate as a function of orbital phase in 2 eclipsing binaries on which a large number of flares are detected. It appears that there is no correlation between flaring rate and orbital phase in these 2 binaries. In contrast, when we examine the function with 203 flares on 20 non-eclipse ellipsoidal binaries, bimodal distribution of amplitude weighted flare numbers shows up at orbital phase 0.25 and 0.75. Such variation could be larger than what is expected from the cross-section modification.

SONG is a global ground based network of 1 meter telescopes for stellar time-domain science, an international collaboration involving many countries across the world. In order to enable a favourable duty cycle, the SONG network plans to create a homogeneous distribution of 4 nodes in each of the northern and southern hemispheres. A natural possibility was building one of the northern nodes in East Asia, preferably on the Qinghai-Tibetan Plateau. During the last decade, a great deal of effort has been invested in searching for high a quality site for ground based astronomy in China, since this has been one of the major concerns for the development of Chinese astronomy. A number of sites on the plateau have been in operation for many years, but most of them are used only for radio astronomy, as well as small optical telescopes for applied astronomy. Several potential sites for large optical instruments have been identified by the plateau site survey, but as yet none of them have been adequately quantitatively characterised. Here we present results from a detailed multi-year study of the Delingha site, which was eventually selected for the SONG-China node. We also describe the site monitoring system that will allow an isolated SONG and 50BiN node to operate safely in an automated mode.

Stars in star clusters are thought to form in a single burst from a common progenitor cloud of molecular gas. However, massive, old globular clusters -- with ages greater than 10 billion years and masses of several hundred thousand solar masses -- often harbour multiple stellar populations, indicating that more than one star-forming event occurred during their lifetimes. Colliding stellar winds from late-stage, asymptotic-giant-branch stars are often invoked as second-generation star-formation trigger. The initial cluster masses should be at least 10 times more massive than they are today for this to work. However, large populations of clusters with masses greater than a few million solar masses are not found in the local Universe. Here we report on three 1-2 billion-year-old, massive star clusters in the Magellanic Clouds, which show clear evidence of burst-like star formation that occurred a few hundred million years after their initial formation era. We show that such clusters could accrete sufficient gas reservoirs to form new stars if the clusters orbited in their host galaxies' gaseous discs throughout the period between their initial formation and the more recent bursts of star formation. This may eventually give rise to the ubiquitous multiple stellar populations in globular clusters.

Supernova remnants (SNRs) are the most attractive candidates for the acceleration sites of Galactic cosmic rays. We report the detection of GeV $\gamma$-ray emission with the Pass 8 events recorded by Fermi Large Area Telescope (Fermi-LAT) in the vicinity of the shell type SNR CTB 37B that is likely associated with the TeV $\gamma-$ray source HESS J1713-381. The photon spectrum of CTB 37B is consistent with a power-law with an index of $1.89\pm0.08$ in the energy range of $0.5-500$ GeV, and the measured flux connects smoothly with that of HESS J1713-381 at a few hundred GeV. No significant spatial extension and time variation are detected. The multi-wavelength data can be well fitted with either a leptonic model or a hadronic one. However, parameters of both models suggest more efficient particle acceleration than typical SNRs. Meanwhile, the X-ray and $\gamma$-ray spectral properties of CTB 37B show that it is an interesting source bridging young SNRs dominated by non-thermal emission and old SNRs interacting with molecular clouds.

Photometric and spectroscopic observations of a slowly declining, luminous Type Ia supernova (SN Ia) SN 2011hr in the starburst galaxy NGC 2691 are presented. SN 2011hr is found to peak at $M_{B}=-19.84 \pm 0.40\,\rm{mag}$, with a post-maximum decline rate $\Delta$m$_{15}$(B) = 0.92 $\pm$ 0.03\,$\rm{mag}$. From the maximum-light bolometric luminosity, $L=(2.30 \pm 0.90) \times 10^{43}\,\rm{erg\,s^{-1}}$, we estimate the mass of synthesized \Nifs\ in SN 2011hr to be $M(\rm{^{56}Ni})=1.11 \pm 0.43\,M_{\sun}$. SN 2011hr appears more luminous than SN 1991T at around maximum light, and the absorption features from its intermediate-mass elements (IMEs) are noticeably weaker than the latter at similar phases. Spectral modeling suggests that SN 2011hr has the IMEs of $\sim$\,0.07 M$_{\sun}$ in the outer ejecta, which is much lower than the typical value of normal SNe Ia (i.e., 0.3 -- 0.4 M$_{\sun}$) and is also lower than the value of SN 1991T (i.e., $\sim$\,0.18 M$_{\sun}$). These results indicate that SN 2011hr may arise from a Chandrasekhar-mass white dwarf progenitor that experienced a more efficient burning process in the explosion. Nevertheless, it is still possible that SN 2011hr may serve as a transitional object connecting the SN 1991T-like SNe Ia with the superluminous subclass like SN 2007if given that the latter also shows very weak IMEs at all phases.

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

We systematically analyze the 6-year {\it Fermi}/LAT data of the lobe-dominated quasars (LDQs) in the complete LDQ sample from 3CRR survey and report the discovery of high-energy $\gamma$-ray emission from 3C 275.1. The $\gamma$-ray emission of 3C 207 is confirmed and significant variability of the lightcurve is identified. We do not find statistically significant $\gamma$-ray emission from other LDQs. 3C 275.1 is the known $\gamma$-ray quasar with the lowest core dominance parameter (i.e., $R=0.11$). We also show that both the northern radio hotspot and parsec jet models can reasonably reproduce the $\gamma$-ray data. The parsec jet model, however, is favored by the potential $\gamma$-ray variability at the timescale of months. We suggest that some dimmer $\gamma$-ray LDQs will be detected in the future and LDQs could contribute non-negligibly to the extragalactic $\gamma$-ray background.

We present optical observations of a SN 2002cx-like supernova SN 2013en in UGC 11369, spanning from a phase near maximum light (t= +1 d) to t= +60 d with respect to the R-band maximum. Adopting a distance modulus of mu=34.11 +/- 0.15 mag and a total extinction (host galaxy+Milky Way) of $A_V \sim1.5$ mag, we found that SN 2013en peaked at $M(R)\sim -18.6$ mag, which is underluminous compared to the normal SNe Ia. The near maximum spectra show lines of Si II, Fe II, Fe III, Cr II, Ca II and other intermediate-mass and iron group elements which all have lower expansion velocities (i.e., ~ 6000 km/s). The photometric and spectroscopic evolution of SN 2013en is remarkably similar to those of SN 2002cx and SN 2005hk, suggesting that they are likely to be generated from a similar progenitor scenario or explosion mechanism.

Two separated sequences of blue straggler stars (BSSs) have been revealed by Ferraro et al. (2009) in the color-magnitude diagram (CMD) of the Milky Way globular cluster M30. Their presence has been suggested to be related to the two BSS formation channels (namely, collisions and mass-transfer in close binaries) operating within the same stellar system. The blue sequence was indeed found to be well reproduced by collisional BSS models. In contrast, no specific models for mass transfer BSSs were available for an old stellar system like M30. Here we present binary evolution models, including case-B mass transfer and binary merging, specifically calculated for this cluster. We discuss in detail the evolutionary track of a $0.9+0.5 M_\odot$ binary, which spends approximately 4 Gyr in the BSS region of the CMD of a 13 Gyr old cluster. We also run Monte-Carlo simulations to study the distribution of mass transfer BSSs in the CMD and to compare it with the observational data. Our results show that: (1) the color and magnitude distribution of synthetic mass transfer BSSs defines a strip in the CMD that nicely matches the observed red BSS sequence, thus providing strong support to the mass transfer origin for these stars; (2) the CMD distribution of synthetic BSSs never attains the observed location of the blue BSS sequence, thus reinforcing the hypothesis that the latter formed through a different channel (likely collisions); (3) most ($\sim 60\%$) of the synthetic BSSs are produced by mass-transfer models, while the remaining $< 40\%$ requires the contribution from merger models.

Only very few z>5 quasars discovered to date are radio-loud, with a radio-to-optical flux ratio (radio-loudness parameter) higher than 10. Here we report the discovery of an optically luminous radio-loud quasar, SDSS J013127.34-032100.1 (J0131-0321 in short), at z=5.18+-0.01 using the Lijiang 2.4m and Magellan telescopes. J0131-0321 has a spectral energy distribution consistent with that of radio-loud quasars. With an i-band magnitude of 18.47 and radio flux density of 33 mJy, its radio-loudness parameter is ~100. The optical and near-infrared spectra taken by Magellan enable us to estimate its bolometric luminosity to be L_bol ~ 1.1E48 erg/s, approximately 4.5 times greater than that of the most distant quasar known to date. The black hole mass of J0131-0321 is estimated to be 2.7E9 solar masses, with an uncertainty up to 0.4 dex. Detailed physical properties of this high-redshift, radio-loud, potentially super-Eddington quasar can be probed in the future with more dedicated and intensive follow-up observations using multi-wavelength facilities.

Optical and ultraviolet observations for the nearby type II-plateau supernova (SN IIP) 2013am in the nearby spiral galaxy M65 are presented in this paper. The early spectra are characterized by relatively narrow P-Cygni features, with ejecta velocities much lower than observed in normal SNe IIP (i.e., $\sim$2000 km s$^{-1}$ vs. $\sim$5000 km $^{-1}$ in the middle of the plateau phase). Moreover, prominent Ca II absorptions are also detected in SN 2013am at relatively early phases. These spectral features are reminiscent of those seen in the low-velocity and low-luminosity SN IIP 2005cs. However, SN 2013am exhibits different photometric properties, having shorter plateau phases and brighter light-curve tails if compared to SN 2005cs. Adopting $R_{V}$=3.1 and a mean value of total reddening derived from the photometric and spectroscopic methods(i.e., $E(B-V)=0.55\pm$0.19 mag), we find that SN 2013am may have reached an absolute $V$-band peak magnitude of $-15.83\pm0.71$ mag, and produced a $^{56}$Ni mass of $0.016^{+0.010}_{-0.006} solar masses in the explosion. These parameters are close to those derived for SN 2008in and SN 2009N which have been regarded as "gap-filler" objects linking the faint SNe IIP to the normal ones. This indicates that some low-velocity SNe IIP may not necessarily result from the low-energetic explosions, and the low expansion velocities could be due to a lower metallicity of the progenitor stars, a larger envelope mass ejected in the explosion, or that is was observed at an angle that is away from the polar direction.

Blue straggler stars are the most prominent bright objects in the colour-magnitude diagram of a star cluster that challenges the theory of stellar evolution. Star clusters are the closest counterparts of the theoretical concept of simple stellar populations (SSPs) in the Universe. SSPs are widely used as the basic building blocks to interpret stellar contents in galaxies. The concept of an SSP is a group of coeval stars which follows a given distribution in mass, and has the same chemical property and age. In practice, SSPs are more conveniently made by the latest stellar evolutionary models of single stars. In reality, however, stars can be more complicated than just single either at birth time or during the course of evolution in a typical environment. Observations of star clusters show that there are always exotic objects which do not follow the predictions of standard theory of stellar evolution. Blue straggler stars (BSSs), as discussed intensively in this book both observationally and theoretically, are very important in our context when considering the integrated spectral properties of a cluster, or a simple stellar population. In this chapter, we are going to describe how important the contribution of BSSs is to the total light of a cluster.

We present a support vector machine classifier to identify the K giant stars from the LAMOST survey directly using their spectral line features. The completeness of the identification is about 75% for tests based on LAMOST stellar parameters. The contamination in the identified K giant sample is lower than 2.5%. Applying the classification method to about 2 million LAMOST spectra observed during the pilot survey and the first year survey, we select 298,036 K giant candidates. The metallicities of the sample are also estimated with uncertainty of $0.13\sim0.29$\,dex based on the equivalent widths of Mg$_{\rm b}$ and iron lines. A Bayesian method is then developed to estimate the posterior probability of the distance for the K giant stars, based on the estimated metallicity and 2MASS photometry. The synthetic isochrone-based distance estimates have been calibrated using 7 globular clusters with a wide range of metallicities. The uncertainty of the estimated distance modulus at $K=11$\,mag, which is the median brightness of the K giant sample, is about 0.6\,mag, corresponding to $\sim30$% in distance. As a scientific verification case, the trailing arm of the Sagittarius stream is clearly identified with the selected K giant sample. Moreover, at about 80\,kpc from the Sun, we use our K giant stars to confirm a detection of stream members near the apo-center of the trailing tail. These rediscoveries of the features of the Sagittarius stream illustrate the potential of the LAMOST survey for detecting substructures in the halo of the Milky Way.

We present multiwavelength observations of the afterglow of GRB 130427A, the brightest (in total fluence) gamma-ray burst of the past 29 years. Optical spectroscopy from Gemini-North reveals the redshift of the GRB to be z=0.340, indicating that its unprecedented brightness is primarily the result of its relatively close proximity to Earth; the intrinsic luminosities of both the GRB and its afterglow are not extreme in comparison to other bright GRBs. We present a large suite of multiwavelength observations spanning from 300 s to 130 d after the burst and demonstrate that the afterglow shows relatively simple, smooth evolution at all frequencies with no significant late-time flaring or rebrightening activity. The entire dataset from 1 GHz to 10 GeV can be modeled as synchrotron emission from a combination of reverse and forward shocks in good agreement with the standard afterglow model, providing strong support to the applicability of the underlying theory and clarifying the nature of the GeV emission observed to last for minutes to hours following other very bright GRBs. A tenuous, wind-stratified circumburst density profile is required by the observations, suggesting a massive-star progenitor with a low mass-loss rate, perhaps due to low metallicity. GRBs similar in nature to GRB 130427A, inhabiting low-density media and exhibiting strong reverse shocks, are probably not uncommon but may have been difficult to recognize in the past due to their relatively faint late-time radio emission; more such events should be found in abundance by the new generation of sensitive radio and millimeter instruments.

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

Observations show that nearly all star clusters and stellar populations contain blue straggler stars (BSs). BSs in a cluster can significantly enhance the integrated spectrum of the host population, preferentially at short wavelengths, and render it much bluer in photometric colours. Current theoretical simple stellar population (SSP) models constructed within the traditional framework of single and binary stellar evolution cannot fully account for the impact of these objects on the integrated spectral properties of stellar populations. Using conventional SSP models without taking into account BS contributions may significantly underestimate a cluster's age and/or metallicity, simply because one has to balance the observed bluer colours (or a bluer spectrum) with a younger age and/or a lower metallicity. Therefore, inclusion of BS contributions in SSP models is an important and necessary improvement for population synthesis and its applications. Here, we present a new set of SSP models, which include BS contributions based on our analysis of individual star clusters. The models cover the wavelength range from 91~{\AA} to 160~$\mu$m, ages from 0.1 to 20 Gyr and metallicities $Z=0.0004, 0.004, 0.008, 0.02$ (solar) and 0.05. We use the observed integrated spectra of several Magellanic Cloud star clusters to cross-check and validate our models. The results show that the age predictions from our models are closer to those from isochrone fitting in the clusters' colour-magnitude diagrams compared to age predictions based on standard SSP models.

(abridged) We present a comprehensive study of stellar stratification in young star clusters in the Large Magellanic Cloud (LMC). We apply our recently developed effective radius method for the assessment of stellar stratification on imaging data obtained with the Advanced Camera for Surveys of three young LMC clusters to characterize the phenomenon and develop a comparative scheme for its assessment in such clusters. The clusters of our sample, NGC 1983, NGC 2002 and NGC 2010, are selected on the basis of their youthfulness, and their variety in appearance, structure, stellar content, and surrounding stellar ambient. Our photometry is complete for magnitudes down to m_814 ~ 23 mag, allowing the calculation of the structural parameters of the clusters, the estimation of their ages and the determination of their stellar content. Our study shows that each cluster in our sample demonstrates stellar stratification in a quite different manner and at different degree from the others. Specifically, NGC 1983 shows to be partially segregated only for the faintest stars of the cluster, NGC 2002 shows evidence of strong stellar stratification for both bright and faint stars, and NGC 2010 is found not to be segregated. For the parametrization of the phenomenon of stellar stratification and its quantitative comparison among these clusters, we propose the slope derived from the change in the effective radius over the corresponding magnitude range as indicative parameter of the degree of stratification in the clusters. A positive value of this slope indicates mass segregation in the cluster, while a negative or zero value signifies the lack of the phenomenon.

We propose a new method for the characterization of stellar stratification in stellar systems. The method uses the mean-square radius (also called the Spitzer radius) of the system as a diagnostic tool. An estimate of the observable counterpart of this radius for stars of different magnitude ranges is used as the effective radius of each stellar species in a star cluster. We explore the dependence of these radii on magnitude as a possible indication of stellar stratification. This method is the first of its kind to use a dynamically stable radius, and though seemingly trivial it has never been applied before. We test the proposed method using model star clusters, which are constructed to be segregated on the basis of a Monte Carlo technique, and on Hubble Space Telescope observations of mass-segregated star clusters in order to explore the limitations of the method in relation to actual data. We conclude that the method performs efficiently in the detection of stellar stratification and its results do not depend on the data, provided that incompleteness has been accurately measured and the contamination by the field population has been thoroughly removed. Our diagnosis method is also independent of any model or theoretical prediction, in contrast to the `classical' methods used so far for the detection of mass segregation.

The presence of blue straggler stars (BSs) in star clusters has proven a challenge to conventional simple stellar population (SSP) models. Conventional SSP models are based on the evolution theory of single stars. Meanwhile, the typical locations of BSs in the colour-magnitude diagram of a cluster are brighter and bluer than the main sequence turn-off point. Such loci cannot be predicted by single-star evolution theory. However, stars with such properties contribute significantly to the integrated light of the cluster. In this paper, we reconstruct the integrated properties of the Large Magellanic Cloud cluster ESO 121-SC03, based on a detailed exploration of the individual cluster stars, and with particular emphasis on the cluster's BSs. We find that the integrated light properties of ESO 121-SC03 are dramatically modified by its BS component. The integrated spectral energy distribution (ISED) flux level is significantly enhanced toward shorter wavelengths, and all broad-band colours become bluer. When fitting the fully integrated ISED of this cluster based on conventional SSP models, the best-fitting values of age and metallicity are significantly underestimated compared to the true cluster parameters. The age underestimate is $\sim40$ per cent if we only include the BSs within the cluster's half-light radius and $\sim60$ per cent if all BSs are included. The corresponding underestimates of the cluster's metallicity are $\sim30$ and $\sim60$ per cent, respectively. The populous star clusters in the Magellanic Clouds are ideal objects to explore the potential importance of BSs for the integrated light properties of more distant unresolved star clusters in a statistically robust manner, since they cover a large range in age and metallicity.

Blue straggler stars present as secure members in the Galactic open clusters form a major challenge to the conventional picture of evolutionary population synthesis based on the stellar evolution theory of single stars, as illustrated in our previous work. Expansion of our sample in the current work to include younger age clusters provides a larger data base to expose the question raised for the simple stellar population model. The working sample now includes 97 Galactic open clusters of ages ranging from 0.1 to several Gyrs. The contributions of blue straggler stars to the integrated light of the host clusters are calculated on an individual cluster base. A data base of observational constrained simple stellar population model is made which has a larger age coverage than our previous work. It is shown in this work that the general existence of blue stragglers in star clusters of our sample dramatically altered the predictions of convectional stellar population model in terms of spectral energy distribution. The integrated spectral energy distributions of the synthetic spectra of the clusters are enhanced towards shorter wavelengths, therefore the results of the present work will cast new lights in understanding the properties of stellar populations.

Using a volume-limited sample of Main Galaxies from the SDSS Data Release 5, we investigate the dependence of galaxy properties on local environment. For each galaxy, the local three-dimensional density is calculated. We find that galaxy morphologies strongly depend on local environment: galaxies in dense environments have predominantly early type morphologies, but other galaxy properties do not present significant dependence on local environment. Clearly, this puts a important constraint on proposed physical mechanisms.

Synthetic integrated spectral properties of the old Galactic open clusters are studies in this work, where twenty-seven Galactic open clusters of ages >= 1Gyr are selected as the working sample. Based on the photometric observations of these open clusters, synthetic integrated spectrum has been made for the stellar population of each cluster. The effects of blue straggler stars (BSSs) on the conventional simple stellar population (SSP) model are analyzed on an individual cluster base. It is shown that the BSSs, whose holding positions in the color-magnitude diagrams (CMDs) cannot be predicted by the current single-star evolution theory, present significant modifications to the integrated properties of theoretical SSP model. The synthesized integrated spectral energy distributions (ISEDs) of our sample clusters are dramatically different from the SSPs based on isochrone only. The BSSs corrected ISEDs of stellar populations show systematic enhancements towards shorter wavelength in the spectra. When measured with wide-band colors in unresolvable conditions, the age of a stellar population can be seriously under-estimated by the conventional SSP model. Therefore, considering the common existence of BSS component in real stellar populations, a considerable amount of alternations on the conventional ISEDs should be expected when applying the technique of evolutionary population synthesis (EPS) to more complicated stellar systems.