For some neutron stars (NSs) in the binary systems, the masses have been accurately measured. While for the isolated neutron stars (INSs), no mass measurement has been reported yet. The situation will change soon thanks to the successful performance of the Neutron Star Interior Composition Explorer (NICER), with which the radius and mass of the isolated PSR J0030+0451 can be simultaneously measured. For most INSs, no mass measurements are possible for NICER because of observational limitations. Benefiting from recent significant progress made on constraining the equation of state of NSs, in this work we propose a way to estimate the masses of the INSs with the measured gravitational redshifts. We apply our method to RX J1856.5-3754, RX J0720.4-3125, and RBS 1223, three members of "The Magnificent Seven" (M7), and estimate their masses to be $1.24_{-0.29}^{+0.29}M_{\odot}$, $1.23_{-0.05}^{+0.10}M_{\odot}$, and $1.08_{-0.11}^{+0.20}M_{\odot}$, respectively. These masses are consistent with that of binary NS systems, suggesting no evidence for experiencing significant accretion of these isolated objects.

The hadronic interaction of cosmic rays with solar atmosphere can produce high energy gamma rays. The gamma-ray luminosity is correlated both with the flux of primary cosmic rays and the intensity of the solar magnetic field. The gamma rays below 200 GeV have been observed by $Fermi$ without any evident energy cutoff. The bright gamma-ray flux above 100 GeV has been detected only during solar minimum. The only available data in TeV range come from the HAWC observations, however outside the solar minimum. The ARGO-YBJ dataset has been used to search for sub-TeV/TeV gamma rays from the Sun during the solar minimum from 2008 to 2010, the same time period covered by the Fermi data. A suitable model containing the Sun shadow, solar disk emission and inverse-Compton emission has been developed, and the chi-square minimization method was used to quantitatively estimate the disk gamma-ray signal. The result shows that no significant gamma-ray signal is detected and upper limits to the gamma-ray flux at 0.3$-$7 TeV are set at 95\% confidence level. In the low energy range these limits are consistent with the extrapolation of the Fermi-LAT measurements taken during solar minimum and are compatible with a softening of the gamma-ray spectrum below 1 TeV. They provide also an experimental upper bound to any solar disk emission at TeV energies. Models of dark matter annihilation via long-lived mediators predicting gamma-ray fluxes > $10^{-7}$ GeV $cm^{-2}$ $s^{-1}$ below 1 TeV are ruled out by the ARGO-YBJ limits.

Using photometric and spectroscopic data of Supernova (SN) 2012cu, a fairly reddened type Ia supernova, we derived its color excess curves and probed the dust extinction in its host galaxy, NGC 4772. In order to derive the extinction as a function of wavelength (i.e., $A_\lambda$), we model the color excess curves of SN 2012cu in terms of dust models consisting of silicate and carbonaceous (graphite or amorphous carbon) dust. The modeled extinction law towards SN 2012cu extends flatly to the far-ultraviolet (UV) bands, which is much flatter than those of the Milky Way and Magellanic Clouds, and the 2175 $\AA$ feature is very weak or absent. The flatness of the modeled extinction curve in the UV bands suggests a "grey" extinction law of the active galactic nucleus in the vicinity of the SN 2012cu-Earth line of sight. Our results indicate that the sizes of the dust in the ISM towards SN 2012cu in NGC 4772 are larger than those of the Milky Way and the Large Magellanic Cloud, and much larger than that of the Small Magellanic Cloud. The best fitting gives an observed visual extinction towards SN 2012cu of $A_V \approx 2.6$ mag, a reddening of $E(B-V) \approx 1.0$ mag, with a total-to-selective extinction ratio $R_V \approx 2.7$, consistent with previous results.

This paper reports on the observation of the sidereal large-scale anisotropy of cosmic rays using data collected by the ARGO-YBJ experiment over 5 years (2008$-$2012). This analysis extends previous work limited to the period from 2008 January to 2009 December,near the minimum of solar activity between cycles 23 and 24.With the new data sample the period of solar cycle 24 from near minimum to maximum is investigated. A new method is used to improve the energy reconstruction, allowing us to cover a much wider energy range, from 4 to 520 TeV. Below 100 TeV, the anisotropy is dominated by two wide regions, the so-called "tail-in" and "loss-cone" features. At higher energies, a dramatic change of the morphology is confirmed. The yearly time dependence of the anisotropy is investigated. Finally, no noticeable variation of cosmic-ray anisotropy with solar activity is observed for a median energy of 7 TeV.

A correlation between the secondary cosmic ray flux and the near-earth electric field intensity, measured during thunderstorms, has been found by analyzing the data of the ARGO-YBJ experiment, a full coverage air shower array located at the Yangbajing Cosmic Ray Laboratory (4300 m a. s. l., Tibet, China). The counting rates of showers with different particle multiplicities, have been found to be strongly dependent upon the intensity and polarity of the electric field measured during the course of 15 thunderstorms. In negative electric fields (i.e. accelerating negative charges downwards), the counting rates increase with increasing electric field strength. In positive fields, the rates decrease with field intensity until a certain value of the field EFmin (whose value depends on the event multiplicity), above which the rates begin increasing. By using Monte Carlo simulations, we found that this peculiar behavior can be well described by the presence of an electric field in a layer of thickness of a few hundred meters in the atmosphere above the detector, which accelerates/decelerates the secondary shower particles of opposite charge, modifying the number of particles with energy exceeding the detector threshold. These results, for the first time, give a consistent explanation for the origin of the variation of the electron/positron flux observed for decades by high altitude cosmic ray detectors during thunderstorms.

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.

The merger of two dense stellar remnants including at least one neutron star (NS) is predicted to produce gravitational waves (GWs) and short duration gamma ray bursts (GRBs). In the process, neutron-rich material is ejected from the system and heavy elements are synthesized by r-process nucleosynthesis. The radioactive decay of these heavy elements produces additional transient radiation termed "kilonova" or "macronova". We report the detection of linear optical polarization P = (0.50 +/- 0.07)% at 1.46 days after detection of the GWs from GW170817, a double neutron star merger associated with an optical macronova counterpart and a short GRB. The optical emission from a macronova is expected to be characterized by a blue, rapidly decaying, component and a red, more slowly evolving, component due to material rich of heavy elements, the lanthanides. The polarization measurement was made when the macronova was still in its blue phase, during which there is an important contribution from a lanthanide-free outflow. The low degree of polarization is consistent with intrinsically unpolarized emission scattered by Galactic dust, suggesting a symmetric geometry of the emitting region and low inclination of the merger system. Stringent upper limits to the polarization degree from 2.45 - 9.48 days post-burst are consistent with the lanthanides-rich macronova interpretation.

The ARGO-YBJ detector, located at the Yangbajing Cosmic Ray Laboratory (4300 m a. s. l., Tibet, China), was a full coverage air shower array dedicated to gamma ray astronomy and cosmic ray studies. The wide field of view (~ 2 sr) and high duty cycle (> 86%), made ARGO-YBJ suitable to search for short and unexpected gamma ray emissions like gamma ray bursts (GRBs). Between 2007 November 6 and 2013 February 7, 156 satellite-triggered GRBs (24 of them with known redshift) occurred within the ARGO-YBJ field of view. A search for possible emission associated to these GRBs has been made in the two energy ranges 10-100 GeV and 10-1000 GeV. No significant excess has been found in time coincidence with the satellite detections nor in a time window of one hour after the bursts. Taking into account the EBL absorption, upper limits to the energy fluence at 99% of confidence level have been evaluated,with values ranging from ~ 10-5 erg cm-2 to ~10-1 erg cm-2.

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.

A toy detector array is designed to detect a shower generated by the interaction between a TeV cosmic ray and the atmosphere. In the present paper, the primary energies of showers detected by the detector array are reconstructed with the algorithm of Bayesian neural networks (BNNs) and a standard method like the LHAASO experiment \cite{lhaaso-ma}, respectively. Compared to the standard method, the energy resolutions are significantly improved using the BNNs. And the improvement is more obvious for the high energy showers than the low energy ones.

Part of a review paper entitled "Gravitational wave astronomy: the current status.", appeared in " Science China Physics, Mechanics & Astronomy 58.12 (2015): 1-41.

A direct WIMP (Weakly Interacting Massive Particle) detector with a neutron veto system is designed to better reject neutrons. An experimental configuration is studied in the present paper: 984 Ge modules are placed inside a reactor neutrino detector. In order to discriminate between nuclear and electron recoil, both ionization and heat signatures are measured using cryogenic germanium detectors in this detection. The neutrino detector is used as a neutron veto device.The neutron background for the experimental design has been estimated using the Geant4 simulation. The results show that the neutron background can decrease to O(0.01) events per year per tonne of high purity Germanium. We calculate the sensitivity to spin-independent WIMP-nucleon elastic scattering. An exposure of one tonne $\times$ year could reach a cross-section of about 2$\times$$10^{-11}$ pb.

The present work reports on a feasibility study commissioned by the Chinese Academy of Sciences of China to explore various possible mission options to detect gravitational waves in space alternative to that of the eLISA/LISA mission concept. Based on the relative merits assigned to science and technological viability, a few representative mission options descoped from the ALIA mission are considered. A semi-analytic Monte Carlo simulation is carried out to understand the cosmic black hole merger histories starting from intermediate mass black holes at high redshift as well as the possible scientific merits of the mission options considered in probing the light seed black holes and their coevolution with galaxies in early Universe. The study indicates that, by choosing the armlength of the interferometer to be three million kilometers and shifting the sensitivity floor to around one-hundredth Hz, together with a very moderate improvement on the position noise budget, there are certain mission options capable of exploring light seed, intermediate mass black hole binaries at high redshift that are not readily accessible to eLISA/LISA, and yet the technological requirements seem to within reach in the next few decades for China.

Swift J164449+573451 is a peculiar outburst which is most likely powered by the tidal disruption of a star by a massive black hole. Within the tidal disruption scenario, we show that the periastron distance is considerably smaller than the disruption radius and the outflow should be launched mainly via magnetic activities (e.g., Blandford-Znajek process) otherwise the observed long-lasting X-ray afterglow emission satisfying the relation $L_{X}\propto\dot{M}$ can not be reproduced, where $L_{X}$ is the X-ray luminosity and $\dot{M}$ is the accretion rate. We also suggest that $L_{X}\propto\dot{M}$ may hold in the quick decline phase of Gamma-ray Bursts.

The optical emission of GRB 110205A is distinguished by two flares. In this work we examine two possible scenarios for the optical afterglow emission. In the first scenario, the first optical flare is the reverse shock emission of the main outflow and the second one is powered by the prolonged activity of central engine. We however find out that it is rather hard to interpret the late ($t>0.1$ day) afterglow data reasonably unless the GRB efficiency is very high ($\sim 0.95$). In the second scenario, the first optical flare is the low energy prompt emission and the second one is the reverse shock of the initial outflow. Within this scenario we can interpret the late afterglow emission self-consistently. The reverse shock region may be weakly magnetized and the decline of the second optical flare may be dominated by the high latitude emission, for which strong polarization evolution accompanying the quick decline is possible, as suggested by Fan et al. in 2008. Time-resolved polarimetry by RINGO2-like polarimeters will test our prediction directly.

We constrain the physical composition of the outflows of GRBs 080916C and 090510 with the prompt emission data and find that the former is likely magnetic while the latter may be baryonic. The X-ray and optical afterglow emission of both GRBs can be reasonably fitted using the standard external shock model but the density profiles of the circum-burst medium are different. We also propose a simple method to estimate the number of the seed photons suppose the GeV afterglow photons are due to the inverse Compton radiation of external forward shock electrons. The seed photons needed in the modeling are too many to be realistic for both events. The synchrotron radiation of the forward shock seems able to account for the GeV afterglow data.

Strong optical and near-infrared (NIR) flares were discovered in the afterglow of GRB 080129. Their temporal behaviors, the sudden emergence and the quick disappearance, are rather similar to that of many X-ray flares (for instance, the giant flare of GRB 050502B). We argue that the optical/NIR flares following GRB 080129 are a low energy analogy of the X-ray flares and the most likely interpretation is the "late internal shock model". In this model, both the very sharp decline and the very small ratio between the duration and the occurrence time of the optical/NIR flares in GRB 080129 can be naturally interpreted. The initial Lorentz factor of the flare outflow is found to be $\sim 30$, consistent with the constraint $\leq 120$ set by the forward shock afterglow modeling. Other possibilities, like the reverse shock emission or the radiation from the continued but weaker and weaker collision between the initial GRB outflow material, are disfavored.

We suggest a short-living supermassive magnetar model to account for the X-ray flares following short $\gamma-$ray bursts. In this model, the central engine of the short $\gamma-$ray bursts is a supermassive millisecond magnetar. The X-ray flares are powered by the dipole radiation of the magnetar. When the magnetar has lost a significant part of its angular momentum, it collapses to a black hole and the X-ray flares disappear abruptly. Two important predictions of this model are (i) X-ray flares much more energetic than that detected in GRB 050724 may be detectable in the coming months and years by the XRT onboard {\it Swift}. (ii) The short GRBs with X-ray flares may occur outside of their host galaxy.

Recently, X-ray emission lines have been observed in X-ray afterglows of several $\gamma$-ray bursts. It is a major breakthrough for understanding the nature of the progenitors. It is proposed that the X-ray emission lines can be well explained by the Geometry-Dominated models, but in these models the illuminating angle is much larger than that of the collimated jet of the $\gamma$-ray bursts(GRBs). For GRB 011211, we obtain the illuminating angle is about $\theta\sim45^{\circ}$, while the angle of GRB jet is only $3.6^{\circ}$, so we propose that the outflow of the GRBs with emission lines should have two distinct components. The wide component illuminates the reprocessing material, and produces the emission lines, while the narrow one produces the $\gamma$-ray bursts. The observations show that the energy for producing the emission lines is higher than that of the GRBs. In this case, when the wide component dominates the afterglows, a bump will appear in the GRBs afterglows. For GRB 011211, the emergence time of the bump is less than 0.05 days after the GRB, it is obviously too early for the observation to catch it. With the presence of the X-ray emission lines there should also be a bright emission component between the UV and the soft X-rays. These features can be tested by the $Swift$ satellite in the near future.

Recently iron emission features have been observed in several X-ray afterglows of GRBs. It is found that the energy obtained from the illuminating continuum which produces the emission lines is much higher than that of the main burst.The observation of SN-GRB association indicates a fallback disk should be formed after the supernovae explosion. The disk is optically thick and advection-dominated and dense. We suggest that the delayed injection energy after the initial main burst, much higher than energy of the main burst, causes the re-burst appearance in GRB afterglow and illuminates the region of the disk surface with $\tau\approx1$($\tau$ is the optical depth for the Thomson scatter ) and produces the iron emission line whose luminosity can be up to $10^{45}$ erg$s^{-1}$. The duration of the iron line emission can be $10^{4}-10^{5}$ s. This model can explain the appearance of re-burst and emission lines in GRB afterglow and disappearance of the iron emission lines, and also can naturally solve the problem of higher energy of the illuminating continuum than that of the main burst. This scenario is different from the models put forward to explain the emission lines before, that can be tested by SWIFT satellite.

Extremely powerful emission lines are observed in the X-ray afterglow of several GRBs. The energy contained in the illuminating continuum which is responsible for the line production exceeds 10$^{51}$ erg, much higher than that of the collimated GRBs. It constrains the models which explain the production of X-ray emission lines. In this paper, We argue that this energy can come from a continuous postburst outflow. Focusing on a central engine of highly magnetized millisecond pulsar or magnetar we find that afterglow can be affected by the illuminating continuum, and therefore a distinct achromatic bump may be observed in the early afterglow lightcurves. With the luminosity of the continuous outflow which produces the line emission, we define the upper limit of the time when the bump feature appears. We argue that the reason why the achromatic bumps have not been detected so far is that the bumps should appear at the time too early to be observed.

The variability of gamma-ray burst (GRB) is thought to be correlated with its absolute peak luminosity, and this relation had been used to derive an estimate of the redshifts of GRBs. Recently Amati et al. presented the results of spectral and energetic properties of several GRBs with known redshifts. Here we analyse the properties of two group GRBs, one group with known redshift from afterglow observation, and another group with redshift derived from the luminosity - variability relation. We study the redshift dependence of various GRBs features in their cosmological rest frames, including the burst duration, the isotropic luminosity and radiated energy, and the peak energy $E_p$ of $\nu F_\nu$ spectra. We find that, for these two group GRBs, their properties are all redshift dependent, i.e. their intrinsic duration, luminosity, radiated energy and peak energy $E_p$, are all correlated with the redshift, which means that there are cosmological evolution effects on gamma-ray bursts features, and this can provide an interesting clue to the nature of GRBs. If this is true, then the results also imply that the redshift derived from the luminosity - variability relation may be reliable.

We report the photometry of a peculiar SU UMa-type dwarf nova - ER UMa for ten nights during 1998 December and 1999 March covering a complete rise to the supermaximum and a normal outburst cycle. Superhumps have been found during the rise to the superoutburst. A negative superhump appeared in Dec.22 light curve, while the superhump on the next night became positive and had large amplitude and distinct waveform from that of the previous night. In the normal outburst we captured, superhumps with larger or smaller amplitudes seem to always exist, although it is not necessarily true for every normal outburst. These results show great resemblance with V1159 Ori (Patterson et al. 1995). It is more likely that superhumps occasionally exist at essentially all phases of the eruption cycles of ER UMa stars, which should be considered in modeling.