Using an explicitly isospin-dependent parametric Equation of State (EOS) for the core of neutron stars (NSs) within the Bayesian statistical approach, we infer the EOS parameters of super-dense neutron-rich nuclear matter from three sets of imagined mass-radius correlation data representing typical predictions by various nuclear many-body theories, i.e, the radius stays the same, decreases or increases with increasing NS mass within $\pm 15\%$ between 1.4 M$_{\odot}$ and 2.0 M$_{\odot}$. The corresponding average density increases quickly, slowly or slightly decreases as the NS mass increases from 1.4 M$_{\odot}$ to 2.0 M$_{\odot}$. Using the posterior probability distribution functions (PDFs) of EOS parameters inferred from GW170817 and NICER radius data for canonical NSs as references, we investigate how future radius measurements of massive NS will improve our knowledge about the EOS of super-dense neutron-rich nuclear matter, especially its symmetry energy term, compared to what people have already learned from analyzing the GW170817 and NICER data. While the EOS of symmetric nuclear matter (SNM) inferred from the three data sets are approximately the same, the corresponding high-density symmetry energies at densities above about $2\rho_0$ are very different, indicating that the radii of massive NSs carry reliable information about the high-density behavior of nuclear symmetry energy with little influence from the remaining uncertainties of the SNM EOS.

Within the minimum model of neutron stars (NS) consisting of neutrons, protons and electrons, a new approach is proposed for inferring the symmetry energy of super-dense neutron-rich nucleonic matter above twice the saturation density $\rho_0$ of nuclear matter directly from integrating iteratively barotropic pressures in both neutron stars and heavy-ion reactions. Simultaneously, the proton fraction of NSs at $\beta$ equilibrium is extracted as a function of baryon density from the same procedure. An application of this approach using the NS pressure from GW170817 and the pressure in cold symmetric nuclear matter (SNM) extracted earlier by analyzing nuclear collective flow data in relativistic heavy-ion collisions provides a useful constraining band for the symmetry energy above $2\rho_0$.

Evidence for the central engine of gamma-ray bursts (GRBs) has been collected in the Neil Gehrels \emph{Swift} data. For instance, some GRBs show an internal X-ray plateau followed by very steep decay, which is difficult to be interpreted within the framework of a black hole (BH) central engine, but are consistent within a rapidly spinning magnetar engine picture. The very steep decay at the end of the plateau suggests a sudden cessation of the central engine, which is explained as the collapse of a supra-massive magnetar into a black hole when it spins down. Here we propose that some additional evidence, such as a second X-ray plateau feature would show up, if the fall-back accretion could activate the newborn BH and sufficient energy could be transferred from the newborn BH to the GRB blast wave. With a systematic data analysis for all long GRBs, we find three candidates in \emph{Swift} sample, i.e., GRBs 070802, 090111, and 120213A, whose X-ray afterglow light curves contain two plateaus, with the first one being an internal plateau. We find that in a fairly loose and reasonable parameter space, the second X-ray plateau data for all 3 GRBs could be well interpreted with our proposed model. Future observations are likely to discover more similar events, which could offer more information of the properties of the magnetar as well as the newborn BH.

Early shallow-decaying X-ray afterglows of gamma-ray bursts (GRBs) may be attributed to the dipole radiations of newly-born magnetars. Assuming that the GRB jets powered by magnetars are quasi-universal, we find that the jet structure can be parameterized as a uniform jet with a luminosity of $\log L_{\rm j}/{\rm erg\ s^{-1}}=52.68^{+0.76}_{-0.33}$ (1$\sigma$) and an opening angle $\theta_{\rm j}=2.10_{-1.28}^{+1.90}$ (50\% confidence level), surrounding by a power-law decay component with an index of ${-4.00^{+0.27}_{-0.37}}$ (1$\sigma$). The inferred local GRB rate is $\rho=9.6$ Gpc$^{-3}$ yr$^{-1}$ by including both the typical GRBs and LL-GRBs as the same population. The typical viewing angle is $3.3^{o}$, and may be $20^{o}\sim30^{o}$ for LL-GRBs. The X-ray luminosity function of the dipole radiation wind can be empirically described by a broken power-law function with indices $\beta_1=0.78^{+0.16}_{-0.15}$ and $\beta_2>1.6$ broken at $\log L_{b, w}/{\rm erg\ s^{-1}}=48.51^{+0.53}_{-0.65}$. In case of that the wind outflow is collimated and co-axial with the GRB jet, we find that the wind structure is similar to the GRB jet, i.e., $\log L_{\rm c, w}/{\rm erg\ s^{-1}}=48.38^{+0.30}_{-0.48}$, $\theta_{\rm c, w}={2.65^{o}}_{-1.19^{o}}^{+0.1.73^{o}}$, and $k_{\rm w}=4.57^{+1.21}_{-0.75}$. The observed correlation between the prompt gamma-ray luminosity and X-ray luminosity of the wind may be resulted from the viewing angle effect in such a jet-wind system. Discussion on survey with the X-ray instruments on board the {\em Einstein\ Probe} mission in the soft X-ray band for the jet and wind emission is also presented.

Within the Bayesian framework using available constraining bands on the pressure in symmetric nuclear matter (SNM) derived earlier by others in the density range of 1.3$\rho_0$ to 4.5$\rho_0$ from kaon production and nuclear collective flow data in energetic heavy-ion collisions, we infer the posterior probability distribution functions (PDFs) of SNM incompressibility $K_0$ and skewness $J_0$ using uniform prior PDFs for them in the ranges of $220\leq K_0\leq 260$ MeV and $-800\leq J_0\leq 400$ MeV. The 68\% posterior credible boundaries around the most probable values of $K_0$ and $J_0$ are found to be 222$\pm$2 MeV and -215$\pm$20 MeV, respectively, much narrower than their prior ranges widely used currently in the literature and are consistent with the results of a recent Bayesian analysis of neutron star properties constrained by available X-ray and gravitational wave observations.

The radius $R_{1.4}$ of neutron stars (NSs) with a mass of 1.4 M$_{\odot}$ has been extracted consistently in many recent studies in the literature. Using representative $R_{1.4}$ data, we infer high-density nuclear symmetry energy $E_{\rm{sym}}(\rho)$ and the associated nucleon specific energy $E_0(\rho)$ in symmetric nuclear matter (SNM) within a Bayesian statistical approach using an explicitly isospin-dependent parametric Equation of State (EOS) for nucleonic matter. We found that: (1) The available astrophysical data can already improve significantly our current knowledge about the EOS in the density range of $\rho_0-2.5\rho_0$. In particular, the symmetry energy at twice the saturation density $\rho_0$ of nuclear matter is determined to be $E_{\mathrm{sym}}(2\rho_0)$ =39.2$_{-8.2}^{+12.1}$ MeV at 68\% confidence level. (2) A precise measurement of the $R_{1.4}$ alone with a 4\% 1$\sigma$ statistical error but no systematic error will not improve much the constraints on the EOS of dense neutron-rich nucleonic matter compared to what we extracted from using the available radius data. (3) The $R_{1.4}$ radius data and other general conditions, such as the observed NS maximum mass and causality condition introduce strong correlations for the high-order EOS parameters. Consequently, the high-density behavior of $E_{\rm{sym}}(\rho)$ inferred depends strongly on how the high-density SNM EOS $E_0(\rho)$ is parameterized, and vice versa. (4) The value of the observed maximum NS mass and whether it is used as a sharp cut-off for the minimum maximum mass or through a Gaussian distribution affect significantly the lower boundaries of both the $E_0(\rho)$ and $E_{\rm{sym}}(\rho)$ only at densities higher than about $2.5\rho_0$.

The Gamma-Ray Integrated Detectors (GRID) is a space mission concept dedicated to monitoring the transient gamma-ray sky in the energy range from 10 keV to 2 MeV using scintillation detectors onboard CubeSats in low Earth orbits. The primary targets of GRID are the gamma-ray bursts (GRBs) in the local universe. The scientific goal of GRID is, in synergy with ground-based gravitational wave (GW) detectors such as LIGO and VIRGO, to accumulate a sample of GRBs associated with the merger of two compact stars and study jets and related physics of those objects. It also involves observing and studying other gamma-ray transients such as long GRBs, soft gamma-ray repeaters, terrestrial gamma-ray flashes, and solar flares. With multiple CubeSats in various orbits, GRID is unaffected by the Earth occultation and serves as a full-time and all-sky monitor. Assuming a horizon of 200 Mpc for ground-based GW detectors, we expect to see a few associated GW-GRB events per year. With about 10 CubeSats in operation, GRID is capable of localizing a faint GRB like 170817A with a 90% error radius of about 10 degrees, through triangulation and flux modulation. GRID is proposed and developed by students, with considerable contribution from undergraduate students, and will remain operated as a student project in the future. The current GRID collaboration involves more than 20 institutes and keeps growing. On August 29th, the first GRID detector onboard a CubeSat was launched into a Sun-synchronous orbit and is currently under test.

We offer an explanation for the anti-correlation between the minimum variability timescale ($MTS$) in the prompt emission light curve of gamma-ray bursts (GRBs) and the estimated bulk Lorentz factor of these GRBs, in the context of a magnetically arrested disk (MAD) model. In particular, we show that previously derived limits on the maximum available energy per baryon in a Blandford-Znajek jet leads to a relationship between the characteristic MAD timescale, $t_{MAD}$, in GRBs and the maximum bulk Lorentz factor: $t_{MAD} \propto \Gamma^{-6}$, somewhat steeper than (although within the error bars of) the fitted relationship found in the GRB data. Similarly, the MAD model also naturally accounts for the observed anti-correlation between $MTS$ and gamma-ray luminosity $L$ in the GRB data, and we estimate the accretion rates of the GRB disks (given these luminosities) in the context of this model. Both of these correlations ($MTS-\Gamma$ and $MTS-L$) are also observed in the AGN data, and we discuss the implications of our results in the context of both GRB and blazar systems.

Mergers of neutron star-neutron star (NS-NS) or neutron star-black hole (NS-BH) binaries are candidate sources of gravitational waves (GWs). At least a fraction of the merger remnant should be a stellar mass BH with a sub-relativistic ejecta. A collimated jet is launched via Blandford-Znajek mechanism from the central BH to trigger a short gamma-ray burst (SGRB). At the same time, a near-isotropic wind can be driven by the Blandford-Payne mechanism (BP). In previous work, additional energy injection to the ejecta from the BP mechanism was ignored, and radioactive decay has long been thought as the main source of energy for "kilonova" in the BH scenario. In this letter, we propose that the wind driven by BP mechanism from the new-born BH-disk can heat up and push the ejecta during the prompt emission phase or even at late time when there is mass fall-back. Such a BP-powered merger-nova could be bright at optical band even for a low-luminosity SGRB. The detection of a GW event with a merger product of BH, and accompanied by a bright merger-nova, would be a robust test of our model.

GRB 160625B is an extremely bright GRB with three distinct emission episodes. By analyzing its data observed with the GBM and LAT on board the {\em Fermi} mission, we find that a multi-color black body (mBB) model can be used to fit the spectra of initial short episode (Episode I) very well within the hypothesis of photosphere emission of a fireball model. The time-resolved spectra of its main episode (Episode II), which was detected with both GBM and LAT after a long quiet stage ($\sim 180$ seconds) of the initial episode, can be fitted with a model composing of an mBB component plus a cutoff power-law (CPL) component. This GRB was detected again in the GBM and LAT bands with a long extended emission (Episode III) after a quiet period of $\sim 300$ seconds. The spectrum of Episode III is adequately fitted with a CPL plus a single power-law models, and no mBB component is required. These features may imply that the emission of three episodes are dominated by distinct physics process, i.e., Episode I is possible from cocoon emission surrounding the relativistic jet, Episode II may be from photosphere emission and internal shock of relativistic jet, and Episode III is contributed by internal and external shocks of relativistic jet. On the other hand, both X-ray and optical afterglows are consistent with standard external shocks model.

The X-ray plateau followed by a steep decay ("internal plateau") has been observed in both long and short gamma-ray burst (GRBs), implying a millisecond magnetar operating in some GRBs. The sharp decay at the end of plateau, marking the abrupt cessation of the magnetar central engine, has been considered as the collapse of a supra-massive magnetar to a black hole (BH) when it spins down. If "internal plateau" is indeed the evidence of a magnetar central engine, a natural expectation is a signature from the new-born BH in some candidates. In this work, we find that GRB 070110 is a particular case, which shows a small X-ray bump following its "internal plateau". We interpret the plateau with a spin-down supra-massive magnetar and the X-ray bump with a fall-back BH accretion. This indicates that the new-born BH is likely active in some GRBs. Therefore, GRB 070110-like events may provide a further support to the magnetar central engine model and enable us to investigate the properties of the magnetar as well as the new-born BH.

Recently, two empirical correlations related to the minimum variability timescale ($\rm MTS$) of the lightcures are discovered in gamma-ray bursts (GRBs). One is the anti-correlation between $\rm MTS$ and Lorentz factor $\Gamma$, the other is the anti-correlation between the $\rm MTS$ and gamma-ray luminosity $L_\gamma$. Both the two correlations might be used to explore the activity of the central engine of GRBs. In this paper we try to understand these empirical correlations by combining two popular black hole (BH) central engine models (namely, Blandford \& Znajek mechanism and neutrino-dominated accretion flow). By taking the $\rm MTS$ as the timescale of viscous instability of the neutrino-dominated accretion flow (NDAF), we find that these correlations favor the scenario in which the jet is driven by Blandford-Znajek (BZ) mechanism.

A stellar mass black hole (BH) surrounded by a neutrino-dominated accretion flow (NDAF) has been discussed in a number of works as the central engine of gamma-ray bursts (GRBs). It is widely believed that NDAF cannot liberate enough energy for bright GRBs. However, these works have been based on the assumption of "no torque" boundary condition, which is invalid when the disk is magnetized. In this paper, we present both numerical and analytical solutions for NDAFs with non-zero boundary stresses, and reexamine their properties. We find that NDAF with such boundary torque can be powerful enough to account for those bright short GRBs, energetic long GRBs and ultra-long GRBs. The disk becomes viscously unstable, which makes it possible to interpret the variability of GRB prompt emission and the steep decay phase in the early X-ray afterglow. Finally, we study the gravitational waves radiated from a processing BH-NDAF. We find that the effects of the boundary torque on the strength of the gravitational waves can be ignored.

Recently, the first association between an ultra-long gamma-ray burst (GRB) and a supernova is reported, i.e., GRB 111209A/SN 2011kl, which enables us to investigate the physics of central engines or even progenitors for ultra-long GRBs. In this paper, we inspect the broad-band data of GRB 111209A/SN 2011kl. The late-time X-ray lightcurve exhibits a GRB 121027A-like fall-back bump, suggesting a black hole central engine. We thus propose a collapsar model with fall-back accretion for GRB 111209A/SN 2011kl. The required model parameters, such as the total mass and radius of the progenitor star, suggest that the progenitor of GRB 111209A is more likely a Wolf-Rayet star instead of blue supergiant, and the central engine of this ultra-long burst is a black hole. The implications of our results is discussed.

The continued observations of Sw J1644+57 in X-ray and radio bands accumulated a rich data set to study the relativistic jet launched in this tidal disruption event. We find that the re-brightening feature in the radio light curve can be naturally explained by the two-component jet model. The possible origin of this structured jet are the Blandford-Znajek and Blandford-Payne mechanisms. We also show that this two-component jet model can interpret the two kinds of quasi-periodic variations in the X-ray light curve: a 200 second quasi-periodic oscillation (QPO) and a 2.7-day quasi-periodic variation. The latter is interpreted by a precessing outer jet launched near the Bardeen-Petterson radius of a warped disk. The $\sim$ 200s QPO could be associated with a second, narrower jet sweeping the observer line-of-sight periodically, which is launched from a spinning black hole in the misaligned direction with respect to the black hole's angular momentum.

The variation of the solar diameter in time and in position angle has implications in astrophysics and in general relativity, as the long series of studies attest. The Transits of Venus in 2004 and 2012 have been carefully studied because of the rarity of the phenomenon and its historical importance due the AU measure and to the discovery of Venus atmosphere. The characterization of Venus atmosphere and the measure of the solar diameter to the milliarcsecond level of precision have been studied also from satellite images. The results of the solar diameter measurements made with the observations in Athens (2004) and at the Huairou Solar Observing Station in China (2012) are presented. The topic of the oblateness of the Sun at sunset and its intrinsic value is drafted to introduce the general public to the relativistic relevance of measuring the solar figure, in the occasion of the International Year of Light 2015.

An isolated and fast decayed active region (NOAA 9729) was observed when passing through solar disk. There is only one CME related with it that give us a good opportunity to investigate the whole process of the CME. Filament in this active region rises up rapidly and then hesitates and disintegrates into flare loops. The rising filament from EIT images separates into two parts just before eruption. A new filament reforms several hours later after CME and the axis of this new one rotates clockwise about 22 degrees comparing with that of the former one. We also observed a bright transient J-shaped X-ray sigmoid immediately appears after filament eruption. It quickly develops into a soft X-ray cusp and rises up firstly then drops down. Two magnetic cancelation regions have been observed clearly just before filament eruption. Moreover, the magnetic flux rope erupted as the magnetic helicity approach the maximum and the normalized helicity is -0.036 when the magnetic flux rope erupted, which is close to the prediction value based on the theoretical non-linear force-free model.

The possibility to measure the solar diameter using the transits of Mercury has been exploited to investigate the past three centuries of its evolution and to calibrate these measurements made with satellites. This measurement basically consists to compare the ephemerides of the internal contact timings with the observed timings. The transits of Venus of 2004 and 2012 gave the possibility to apply this method, involving a planet with atmosphere, with the refraction of solar light through it creating a luminous arc all around the disk of the planet. The observations of the 2012 transit made to measure the solar diameter participate to the project Venus Twilight Experiment to study the aureole appearing around it near the ingress/egress phases.

We propose a two-component jet model consistent with the observations of several gamma ray bursts (GRBs) and active galactic nuclei (AGNs). The jet consists of inner and outer components, and they are supposed to be driven by the Blandford-Znajek (BZ) and Blandford-Payne (BP) processes, respectively. The baryons in the BP jet is accelerated centrifugally via the magnetic field anchored in the accretion disk. The BZ jet is assumed to be entrained a fraction of accreting matter leaving the inner edge of the accretion disk, and the baryons are accelerated in the conversion from electromagnetic energy to the kinetic energy. By fitting the Lorentz factors of some GRBs (GRB 030329, GRB 051221A, GRB 080413B) and AGNs (Cen A, Mkn 501 and Mkn 421) with this model, we constrain the physical parameters related to the accretion and outflow of these two kind of objects. We conclude that the spine/sheath structure of the jet from these sources can be interpreted naturally by the BZ and BP processes.