On 2019 August 14, the LIGO and Virgo interferometers detected a high-significance event labelled S190814bv. Preliminary analysis of the GW data suggests that the event was likely due to the merger of a compact binary system formed by a BH and a NS. ElectromagNetic counterparts of GRAvitational wave sources at the VEry Large Telescope (ENGRAVE) collaboration members carried out an intensive multi-epoch, multi-instrument observational campaign to identify the possible optical/near infrared counterpart of the event. In addition, the ATLAS, GOTO, GRAWITA-VST, Pan-STARRS and VINROUGE projects also carried out a search on this event. Our observations allow us to place limits on the presence of any counterpart and discuss the implications for the kilonova (KN) possibly generated by this NS-BH merger, and for the strategy of future searches. Altogether, our observations allow us to exclude a KN with large ejecta mass $M\gtrsim 0.1\,\mathrm{M_\odot}$ to a high ($>90\%$) confidence, and we can exclude much smaller masses in a subsample of our observations. This disfavours the tidal disruption of the neutron star during the merger. Despite the sensitive instruments involved in the campaign, given the distance of S190814bv we could not reach sufficiently deep limits to constrain a KN comparable in luminosity to AT 2017gfo on a large fraction of the localisation probability. This suggests that future (likely common) events at a few hundreds Mpc will be detected only by large facilities with both high sensitivity and large field of view. Galaxy-targeted observations can reach the needed depth over a relevant portion of the localisation probability with a smaller investment of resources, but the number of galaxies to be targeted in order to get a fairly complete coverage is large, even in the case of a localisation as good as that of this event.

Gamma-ray bursts (GRBs) are known to be the most violent explosions in the universe, and a variety of correlations between observable GRB properties have been proposed in literature, but none of these correlations is valid for both long GRBs and short GRBs. In this paper we report the discovery of a universal correlation which is suitable for both long and short GRBs using three prompt emission properties of GRBs, i.e. the isotropic peak luminosity $L_{\rm iso}$, the peak energy of the time-integtated prompt emission spectrum $E_{\rm peak}$, and the "high signal" timescale $T_{\rm 0.45}$, $L_{\rm iso} \propto E_{\rm peak}^{1.94} T_{0.45}^{0.37}$. This universal correlation just involves properties of GRB prompt emission and does not require any information of afterglow phase, which can be used as a relatively unbiased redshift estimator. Here we use this correlation to estimate the pseudo-redshifts for short Gamma Ray Bursts and then use Lynden-Bell method to obtain a non-parametric estimate of their luminosity function and formation rate. The luminosity function is $\psi(L_0)\propto{L_0^{-0.63\pm{0.07}}}$ for dim SGRBs and $\psi(L_0)\propto{L_0^{-1.96\pm{0.28}}}$ for bright SGRBs, with the break point $6.95_{-0.76}^{+0.84}\times10^{50} erg/s$. The local formation rate of SGRBs is about 15 events $\rm Gpc^{-3}yr^{-1}$ . This universal correlation may have important implications for GRB physics, implying that the long and short GRBs should share similar radiation processes.

The post-merger gravitational wave radiation of the remnant formed in the double neutron star coalescence has not been directly measured, yet. We show in this work that the properties of the binary neutron stars involved in GW170817, additionally constrained by PSR J0030+0451, the lower limit on the maximum gravitational mass of non-rotating neutron star and some nuclear data, are in favor of strong post-merger gravitational wave radiation. This conclusion applies to the mergers of Galactic binary neutron star systems as well. Significant post-merger gravitational wave radiation is also preferred to improve the consistency between the maximum gravitational mass of the non-rotating neutron star inferred from GW170817/GRB170817A/AT2017gfo and the latest mass measurements of pulsars. The prominent post-merger gravitational radiation of GW170817-like events are expected to be detectable by advanced LIGO/Virgo detectors in the next decade and then shed valuable lights on the properties of the matter in the extremely high density.

For on-axis typical short gamma-ray bursts (sGRBs), the forward shock emission is usually so bright that renders the identification of kilonovae (also known as macronovae) in the early afterglow ($t<0.5$ d) phase rather challenging. This is why previously no thermal-like kilonova component has been identified at such early time except in the off-axis dim GRB 170817A associated with GW170817. Here we report the identification of an unusual optical radiation component in GRB 070809 at $t\sim 0.47$ d, thanks plausibly to the very-weak/subdominant forward shock emission. The optical emission with a very red spectrum is well in excess of the extrapolation of the X-ray emission that is distinguished by an unusually hard spectrum, which is at odds with the forward shock afterglow prediction but can be naturally interpreted as a kilonova. Our finding supports the speculation that kilonovae are ubiquitous , and demonstrates the possibility of revealing the neutron star merger origin with the early afterglow data of some typical sGRBs that take place well beyond the sensitive radius of the advanced gravitational wave detectors and hence the opportunity of organizing dedicated follow-up observations for events of interest.

Assuming that the differential rotation of the massive neutron star (NS) formed in the double NS (DNS) mergers has been effectively terminated by the magnetic braking and a uniform rotation has been subsequently established (i.e., a supramassive NS is formed), we analytically derive in this work an approximated expression for the critical total gravitational mass ($M_{\rm tot,c}$) to form supramassive NS (SMNS) in the DNS mergers, benefited from some equation of state (EoS) insensitive relationships. The maximum gravitational mass of the nonrotating NSs ($M_{\rm TOV}$) as well as the dimensionless angular momentum of the remnant ($j$) play the dominant roles in modifying $M_{\rm tot,c}$, while the radius and mass differences of the premerger NSs do not. The GW170817/GRB 170817A/AT2017gfo observations have provided so far the best opportunity to quantitatively evaluate $M_{\rm TOV}$. Supposing the central engine for GRB 170817A is a black hole quickly formed in the collapse of an SMNS, we find $M_{\rm TOV}=2.13^{+0.09}_{-0.08}M_\odot$ (68.3% credibility interval, including also the uncertainties of the EoS insensitive relationships), which is consistent with the constraints set by current NS mass measurements.

GW190425 is the second neutron star merger event detected by the Advanced LIGO/Virgo detectors. If interpreted as a double neutron star merger, the total gravitational mass is substantially larger than that of the binary systems identified in the Galaxy. In this work we analyze the gravitational-wave data within the neutron star$-$black hole merger scenario. For the black hole, we yield a mass of $2.40^{+0.36}_{-0.32}M_\odot$ and an aligned spin of $0.141^{+0.067}_{-0.064}$. As for the neutron star we find a mass of $1.15^{+0.15}_{-0.13}M_\odot$ and the dimensionless tidal deformability of $1.4^{+3.8}_{-1.2}\times 10^{3}$. These parameter ranges are for 90\% credibility. The inferred masses of the neutron star and the black hole are not in tension with current observations and we suggest that GW190425 is a viable candidate of a neutron star$-$black hole merger event. Benefitting from the continual enhancement of the sensitivities of the advanced gravitational detectors and the increase of the number of the observatories, similar events are anticipated to be much more precisely measured in the future and the presence of black holes below the so-called mass gap will be unambiguously clarified. If confirmed, the mergers of neutron stars with (quickly rotating) low-mass black holes are likely important production sites of the heaviest r-process elements.

In this paper we report the results of the first $\sim$four years of spectroscopic and photometric monitoring of the Type IIn supernova SN 2015da (also known as PSN J13522411+3941286, or iPTF16tu). The supernova exploded in the nearby spiral galaxy NGC 5337 in a relatively highly extinguished environment. The transient showed prominent narrow Balmer lines in emission at all times and a slow rise to maximum in all bands. In addition, early observations performed by amateur astronomers give a very well-constrained explosion epoch. The observables are consistent with continuous interaction between the supernova ejecta and a dense and extended H-rich circumstellar medium. The presence of such an extended and dense medium is difficult to reconcile with standard stellar evolution models, since the metallicity at the position of SN 2015da seems to be slightly subsolar. Interaction is likely the mechanism powering the light curve, as confirmed by the analysis of the pseudo bolometric light curve, which gives a total radiated energy $\gtrsim10^{51}\,\rm{erg}$. Modeling the light curve in the context of a supernova shock breakout through a dense circumstellar medium allowed us to infer the mass of the prexisting gas to be $\simeq8\,\rm{M_{\odot}}$, with an extreme mass-loss rate for the progenitor star $\simeq0.6\,\rm{M_{\odot}}\,\rm{yr^{-1}}$, suggesting that most of the circumstellar gas was produced during multiple eruptive events. Near- and mid-infrared observations reveal a flux excess in these domains, similar to those observed in SN 2010jl and other interacting transients, likely due to preexisting radiatively heated dust surrounding the supernova. By modeling the infrared excess, we infer a mass $\gtrsim0.4\times10^{-3}\,\rm{M_{\odot}}$ for the dust.

In this work we parameterize the Equation of State of dense neutron star (NS) matter with four pressure parameters of $\{\hat{p}_1, \hat{p}_2, \hat{p}_3, \hat{p}_4\}$ and then set the combined constraints with the data of GW 170817 and the data of 6 Low Mass X-ray Binaries (LMXBs) with thermonuclear burst or alternatively the symmetry energy of the nuclear interaction. We find that the nuclear data effectively narrow down the possible range of $\hat{p}_1$, the gravitational wave data plays the leading role in bounding $\hat{p}_2$, and the LMXB data as well as the lower bound on maximal gravitational mass of non-rotating NSs govern the constraints on $\hat{p}_3$ and $\hat{p}_4$. Using posterior samples of pressure parameters and some universal relations, we further investigate how the current data sets can advance our understanding of tidal deformability ($\Lambda$), moment of inertia ($I$) and binding energy ($BE$) of NSs. For a canonical mass of $1.4M_\odot$, we have $I_{1.4} = {1.43}^{+0.30}_{-0.13} \times 10^{38}~{\rm kg \cdot m^2}$, $\Lambda_{1.4} = 390_{-210}^{+280}$ , $R_{1.4} = 11.8_{-0.7}^{+1.2}~{\rm km}$ and $BE_{1.4} = {0.16}^{+0.01}_{-0.02} M_{\odot}$ if the constraints from the nuclear data and the gravitational wave data have been jointly applied. For the joint analysis of gravitational wave data and the LMXB data, we have $I_{1.4} = {1.28}^{+0.15}_{-0.08} \times 10^{38}~{\rm kg \cdot m^2}$, $\Lambda_{1.4} = 220_{-90}^{+90}$, $R_{1.4} = 11.1_{-0.6}^{+0.7}~{\rm km}$ and $BE_{1.4} = {0.18}^{+0.01}_{-0.01} M_{\odot}$. These results suggest that the current constraints on $\Lambda$ and $R$ still suffer from significant systematic uncertainties, while $I_{1.4}$ and $BE_{1.4}$ are better constrained.

This paper is devoted to the primary spectro-polarimetric observation performed at the New Vacuum Solar Telescope of China since 2017, and our aim is to precisely evaluate the real polarimetric accuracy and sensitivity of this polarimetry by using full Stokes spectro-polarimetric observations of the photospheric line Fe I 532.4 nm. In the work, we briefly describe the salient characteristic of the NVST as a polarimeter in technology and then characterize its instrumental polarization based on the operation in 2017 and 2019. It is verified that the calibration method making use of the instrumental polarization calibration unit (ICU) is stable and credible. The calibration accuracy can reach up to 3$\times 10^{-3}$ . Based on the scientific observation of the NOAA 12645 on April 5th, 2017, we estimate that the residual cross-talk from Stokes $I$ to Stokes $Q$, $U$ and $V$, after the instrumental polarization calibration, is about 4$\times10^{-3}$ on average, which is consistent with the calibration accuracy and close to the photon noise. The polarimetric sensitivity (i.e., the detection limit) for polarized light is of the order of $10^{-3}$ with an integration time over 20 seconds. Slow modulation rate is indeed an issue for the present system. The present NVST polarimeter is expected to be integrated with an high-order adaptive optics system and a field scanner to realize 2D magnetic field vector measurements in the following instrumentation update.

As revealed recently by the modeling of the multi-wavelength data of the emission following GW170817/GRB 170817A, there was an off-axis energetic relativistic outflow component launched by this historic double neutron star merger event. In this work we use the results of these modeling to examine the energy extraction process of the central engine. We show that the magnetic process (i.e., the Blandford-Znajek mechanism) is favored, while the neutrino process usually requires a too massive accretion disk if the duration of the central engine activity is comparable to the observed $T_{90}$ of GRB 170817A, unless the timescale of the central engine activity is less than $\sim$ 0.2s. We propose that the GRB observations are helpful to constrain the combined tidal parameter $\tilde{\Lambda}$, and by adopting the accretion disk mass distribution estimated in BZ mechanism, the $90\%$ credible interval of $\tilde{\Lambda}$ for the progenitor of GW170817 is inferred as $309-954$.

GW170817, the first neutron star merger event detected by advanced LIGO/Virgo detectors, was associated with an underluminous short duration GRB 170817A. In this work we compare the forward shock afterglow emission of GW170817/GRB 170817A to other luminous short GRBs (sGRBs) with both a known redshift and an afterglow emission lasting at least one day after the burst. In the rapid decay phase, the afterglow emission of the bright sGRBs and GW170817/GRB 170817A form a natural and continuous sequence, though separated by an observation time gap. If viewed on-axis, the forward shock afterglow emission of GW170817/GRB 170817A would be among the brightest ones detected so far. This provides a strong evidence for the GW170817-like merger origin of bright sGRBs, and suggests that the detection of the forward shock afterglow emission of most neutron star merger events are more challenging than the case of GW170817 since usually the mergers will be more distant and the viewing angles are plausibly higher.

The Gravitational Wave (GW) event GW 170817 was generated by the coalescence of two neutron stars (NS) and produced an electromagnetic transient, labelled AT 2017gfo, that was target of a massive observational campaign. Polarimetry, a powerful diagnostic tool for probing the geometry and emission processes of unresolved sources, was obtained for this event. The observed linear polarization was consistent with being mostly induced by intervening dust, suggesting that the intrinsic emission was weakly polarized ($P < 0.4-0.5$ %). In this paper, we present and discuss a detailed analysis of the linear polarization expected from a merging NS binary system by means of 3D Monte Carlo radiative transfer simulations assuming a range of possible configurations, wavelengths, epochs and viewing angles. We find that polarization originates from the non-homogeneous opacity distribution within the ejecta and can reach levels of $P\sim1$ % at early times (1$-$2 days after the merger) and in the optical R band. Smaller polarization signals are expected at later epochs and/or different wavelengths. From the viewing-angle dependence of the polarimetric signal, we constrain the observer orientation of AT 2017gfo within $\sim$65$^\circ$ from the polar direction. The detection of non-zero polarization in future events will unambiguously reveal the presence of a lanthanide-free ejecta component and unveil its spatial and angular distribution.

The historic first joint detection of both gravitational wave and electromagnetic emission from a binary neutron star merger cemented the association between short gamma-ray bursts (SGRBs) and compact object mergers, as well as providing a well sampled multi-wavelength light curve of a radioactive kilonova (KN) for the first time. Here we compare the optical and near-infrared light curves of this KN, AT2017gfo, to the counterparts of a sample of nearby (z < 0.5) SGRBs to characterize their diversity in terms of their brightness distribution. Although at similar epochs AT2017gfo appears fainter than every SGRB-associated KN claimed so far, we find three bursts (GRBs 050509B, 061201 and 080905A) where, if the reported redshifts are correct, deep upper limits rule out the presence of a KN similar to AT2017gfo by several magnitudes. Combined with the properties of previously claimed KNe in SGRBs this suggests considerable diversity in the properties of KN drawn from compact object mergers, despite the similar physical conditions that are expected in many NS-NS mergers. We find that observer angle alone is not able to explain this diversity, which is likely a product of the merger type (NS-NS versus NS-BH) and the detailed properties of the binary (mass ratio, spins etc). Ultimately disentangling these properties should be possible through observations of SGRBs and gravitational wave sources, providing direct measurements of heavy element enrichment throughout the Universe.

So far, in total 15 double neutron star systems (DNSs) with a reliable measurement of the total gravitational mass ($M_{\rm T}$) have been detected in the Galaxy. In this work we study the distribution of $M_{\rm T}$. The data prefer the double Gaussian distribution over a single Gaussian distribution and the low and high mass populations center at $M_{\rm T}\sim 2.58M_\odot$ and $\sim 2.72M_\odot$, respectively. The progenitor stars of GW170817 have a $M_{\rm T}=2.74^{+0.04}_{-0.01}M_\odot$, falling into the high mass population. With a local neutron star merger rate of $\sim 10^{3}~{\rm Gpc^{-3}~yr^{-1}}$, supposing the $M_{\rm T}$ of those merging neutron stars also follow the double Gaussian distribution, the upcoming runs of the advanced LIGO/Virgo will soon detect some events with a $M_{\rm T}\lesssim 2.6M_\odot$ that can effectively probe the equation of state of the neutron stars and the distribution function is expected to be reliably reconstructed in the next decade.

The gravitational wave data of GW170817 favor the equation of state (EoS) models that predict compact neutron stars (NSs), consistent with the radius constraints from X-ray observations. Motivated by such a remarkable progress, we examine the fate of the remnants formed in NS mergers and focus on the roles of the angular momentum and the mass distribution of the binary NSs. In the mass shedding limit (for which the dimensionless angular momentum equals to the Keplerian value, i.e., $j=j_{\rm Kep}$), the adopted { seven EoS models, except H4 and ALF2,} yield supramassive NSs in more than half of the mergers. However, for $j\lesssim 0.7j_{\rm Kep}$, the presence or absence of a non-negligible fraction of supramassive NSs formed in the mergers depends sensitively on both the EoS and the mass distribution of the binary systems. The NS mergers with a total gravitational mass $\leq 2.6M_\odot$ are found to be able to shed valuable light on both the EoS model and the angular momentum of the remnants if supramassive NSs are still absent. We have also discussed the uncertainty on estimating the maximum gravitational mass of non-rotating NSs ($M_{\rm max}$) due to the unknown $j$ of the pre-collapse remnants. With the data of GW170817 and the assumption of the mass loss of $0.03M_\odot$, we have $M_{\rm max}<(2.19,~2.32)M_\odot$ (90\% confidence level) for $j=(1.0,~0.8)j_{\rm Kep}$, respectively.

The jet breaks in the afterglow lightcurves of short gamma-ray bursts (SGRBs), rarely detected so far, are crucial for estimating the half-opening angles of the ejecta ($\theta_{\rm j}$) and hence the neutron star merger rate. In this work we report the detection of jet decline behaviors in GRB 150424A and GRB 160821B and find $\theta_{\rm j}\sim 0.1$ rad. Together with five events reported before 2015 and other three "identified" recently (GRB 050709, GRB 060614 and GRB 140903A), we have a sample consisting of nine SGRBs and one long-short GRB with reasonably estimated $\theta_{\rm j}$. In particular, three {\it Swift} bursts in the sample have redshifts $z\leq 0.2$, with which we estimate the local neutron star merger rate density {to be $\sim 1109^{+1432}_{-657}~{\rm Gpc^{-3}~yr^{-1}}$ or $162^{+140}_{-83} {\rm Gpc^{-3}yr^{-1}}$ if the narrowly-beamed GRB 061201 is excluded}. Inspired by the typical $\theta_{\rm j}\sim 0.1$ rad found currently, we further investigate whether the off-beam GRBs (in the uniform jet model) or the off-axis events (in the structured jet model) can significantly enhance the GRB/GW association or not. For the former the enhancement is at most moderate, while for the latter the enhancement can be much greater and a high GRB/GW association probability of $\sim 10\%$ is possible. We also show that the data of GRB 160821B may contain a macronova/kilonova emission component with a temperature of $\sim 3100$ K at $\sim 3.6$ days after the burst and more data are needed to ultimately clarify.

GRB 170817A is the first short gamma-ray burst (GRB) with direct detection of the gravitational-wave radiation and also the spectroscopically identified macronova emission (i.e., AT 2017gfo). The prompt emission of this burst, however, is underluminous in comparison with the other short GRBs with known redshift. In this work, we examine whether GRB 170817A is indeed unique. We firstly show that GRB 130603B/macronova may be the on-axis "analogs" of GRB 170817A/AT 2017gfo, and the extremely dim { but long-lasting} afterglow emission of GRB 170817A may suggest a low number density ($\sim 10^{-5}~{\rm cm^{-3}}$) of its circumburst medium { and a structured outflow}. We then discuss whether GRB 070923, GRB 080121, GRB 090417A, GRB 111005A, and GRB 170817A form a new group of very nearby underluminous GRBs originated from neutron star mergers. If the short events GRB 070923, GRB 080121, and GRB 090417A are indeed at a redshift of $\sim 0.076,~0.046,~0.088$, respectively, their isotropic energies of the prompt emission are $\sim 10^{47}$ erg and thus comparable to the other two events. The non-detection of optical counterparts of GRB 070923, GRB 080121, GRB 090417A, and GRB 111005A, however, strongly suggests that the macronovae from neutron star mergers are significantly diverse in luminosities or, alternatively, there is the other origin channel (for instance, the white dwarf and black hole mergers). We finally suggest that GW170817/GRB 170817A are likely not alone and similar events will be detected by the upgraded/upcoming gravitational-wave detectors and the electromagnetic monitors.

GRB 111005A, one long duration gamma-ray burst (GRB) occurred within a metal-rich environment that lacks massive stars with $M_{\rm ZAMS}\geq 15M_\odot$, is not coincident with supernova emission down to stringent limit and thus should be classified as a "long-short" GRB (lsGRB; also known as SN-less long GRB or hybrid GRB), like GRB 060505 and GRB 060614. In this work we show that in the neutron star merger model, the non-detection of the optical/infrared emission of GRB 111005A requires a sub-relativistic neutron-rich ejecta with the mass of $\leq 0.01~M_\odot$, (significantly) less massive than that of GRB 130603B, GRB 060614 and GRB 050709. The lsGRBs are found to have a high rate density and the neutron star merger origin model can be unambiguously tested by the joint observations of the second generation gravitational wave (GW) detectors and the full-sky gamma-ray monitors such as Fermi-GBM and the proposing GECAM. If no lsGRB/GW association is observed in 2020s, alternative scenarios have to be systematically investigated. With the detailed environmental information achievable for the very-nearby events, a novel kind of merger or explosion origin may be identified.

On 17 August 2017, a gravitational wave event (GW170817) and an associated short gamma-ray burst (GRB 170817A) from a binary neutron star merger had been detected. The followup optical/infrared observations also identified the macronova/kilonova emission (AT2017gfo). In this work we discuss some implications of the remarkable GW170817/GRB 170817A/AT2017gfo association. We show that the $\sim 1.7$s time delay between the gravitational wave (GW) and GRB signals imposes very tight constraint on the superluminal movement of gravitational waves (i.e., the relative departure of GW velocity from the speed of light is $\leq 4.3\times 10^{-16}$) or the possible violation of weak equivalence principle (i.e., the difference of the gamma-ray and GW trajectories in the gravitational field of the galaxy and the local universe should be within a factor of $\sim 3.4\times 10^{-9}$). The so-called Dark Matter Emulators and a class of contender models for cosmic acceleration ("Covariant Galileon") are ruled out, too. The successful identification of Lanthanide elements in the macronova/kilonova spectrum also excludes the possibility that the progenitors of GRB 170817A are a binary strange star system. The high neutron star merger rate (inferred from both the local sGRB data and the gravitational wave data) together with the significant ejected mass strongly suggest that such mergers are the prime sites of heavy r-process nucleosynthesis.

The merger of two neutron stars is predicted to give rise to three major detectable phenomena: a short burst of gamma-rays, a gravitational wave signal, and a transient optical/near-infrared source powered by the synthesis of large amounts of very heavy elements via rapid neutron capture (the r-process). Such transients, named "macronovae" or "kilonovae", are believed to be centres of production of rare elements such as gold and platinum. The most compelling evidence so far for a kilonova was a very faint near-infrared rebrightening in the afterglow of a short gamma-ray burst at z = 0.356, although findings indicating bluer events have been reported. Here we report the spectral identification and describe the physical properties of a bright kilonova associated with the gravitational wave source GW 170817 and gamma-ray burst GRB 170817A associated with a galaxy at a distance of 40 Mpc from Earth. Using a series of spectra from ground-based observatories covering the wavelength range from the ultraviolet to the near-infrared, we find that the kilonova is characterized by rapidly expanding ejecta with spectral features similar to those predicted by current models. The ejecta is optically thick early on, with a velocity of about 0.2 times light speed, and reaches a radius of about 50 astronomical units in only 1.5 days. As the ejecta expands, broad absorption-like lines appear on the spectral continuum indicating atomic species produced by nucleosynthesis that occurs in the post-merger fast-moving dynamical ejecta and in two slower (0.05 times light speed) wind regions. Comparison with spectral models suggests that the merger ejected 0.03-0.05 solar masses of material, including high-opacity lanthanides.

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.

Short Gamma-Ray Bursts(SGRBs) are widely believed to be from mergers of binary compact objects involving at least one neutron star and hence have a broad range of spatial offsets from their host galaxies. In this work we search for possible correlations between the emission properties of 18 SGRBs and their offsets from the host galaxies. The SGRBs with and without extended emission do not show significant difference between their offset distribution, in agreement with some previous works. There are however possible correlations between the optical and X-ray afterglow emission and the offsets. The underlying physical origins are examined.

Intermediate mass black holes (IMBHs) with a mass between $10^{2}$ and $10^{5}$ times that of the sun, which bridges the {mass gap between the} stellar-mass black holes and the supermassive black holes, are crucial in understanding the evolution of the black holes. Although they are widely believed to exist, decisive evidence has long been absent. Motivated by the successful detection of massive stellar-mass black holes by advanced LIGO, through the gravitational wave radiation during the binary merger, in this work we investigate the prospect of detecting/identifying the lightest IMBHs (LIMBHs; the black holes $\gtrsim 100M_\odot$) with the second generation gravitational wave detectors. In general, the chance of hearing the birth of the LIMBHs is significantly higher than that to identify pre-merger IMBHs. The other formation channel of LIMBHs, where stars with huge mass/low-metallicity directly collapse, is likely "silent", so the merger-driven birth of the LIMBHs may be the only "observable" scenario in the near future. Moreover, the prospect of establishing the presence of (lightest) intermediate mass black holes in the O3 run and beyond of advanced LIGO is found quite promising, implying that such an instrument could make another breakthrough on astronomy in the near future. The joining of other detectors like advanced Virgo would only increase the detection rate.

GRB 160625B is an extremely-bright outburst with well-monitored afterglow emission. The geometry-corrected energy is high up to $\sim 5.2\times10^{52}$ erg or even $\sim 8\times 10^{52}$ erg, rendering it the most energetic GRB prompt emission recorded so far. We analyzed the time-resolved spectra of the prompt emission and found that in some intervals there were likely thermal-radiation components and the high energy emission were characterized by significant cutoff. The bulk Lorentz factors of the outflow material are estimated accordingly. We found out that the Lorentz factors derived in the thermal-radiation model are consistent with the luminosity-Lorentz factor correlation found in other bursts as well as in GRB 090902B for the time-resolved thermal-radiation components. While the spectral cutoff model yields much lower Lorentz factors that are in tension with the constraints set by the electron pair Compoton scattering process. We then suggest that these spectral cutoffs are more likely related to the particle acceleration process and that one should be careful in estimating the Lorentz factors if the spectrum cuts at a rather low energy (e.g., $\sim$ tens MeV). The nature of the central engine has also been discussed and a stellar-mass black hole is likely favored.

The modeling of Li-Paczy\'{n}ski macronova/kilonova signals gives reasonable estimate on the neutron-rich material ejected during the neutron star mergers. Usually the accretion disk is more massive than the macronova ejecta, with which the efficiencies of converting the disk mass into prompt emission of three merger-driven GRBs can hence be directly constrained. Supposing the macronovae/kilonovae associated with GRB 050709, GRB 060614 and GRB 130603B arose from radioactive decay of the r-process material, the upper limit on energy conversion efficiencies are found to be as low as $\sim 10^{-6}-10^{-4}$. Moreover, for all three events, neutrino annihilation is likely powerful enough to account for the brief gamma-ray flashes. Neutrino annihilation can also explain the "extended" emission lasting $\sim 100$ s in GRB 050709, but does not work for the one in GRB 060614. These progresses demonstrate that the macronova can serve as a novel probe of the central engine activity.

All next generation ground-based and space-based solar telescopes require a good quality assessment metric in order to evaluate their imaging performance. In this paper, a new image quality metric, the median filter gradient similarity (MFGS) is proposed for photospheric images. MFGS is a no-reference/blind objective image quality metric (IQM) by a measurement result between 0 and 1 and has been performed on short-exposure photospheric images captured by the New Vacuum Solar Telescope (NVST) of the Fuxian Solar Observatory and by the Solar Optical Telescope (SOT) onboard the Hinode satellite, respectively. The results show that: (1)the measured value of MFGS changes monotonically from 1 to 0 with degradation of image quality; (2)there exists a linear correlation between the measured values of MFGS and root-mean-square-contrast (RMS-contrast) of granulation; (3)MFGS is less affected by the image contents than the granular RMS-contrast. Overall, MFGS is a good alternative for the quality assessment of photospheric images.

Many technologies have been developed to help improve spatial resolution of observational images for ground-based solar telescopes, such as adaptive optics (AO) systems and post-processing reconstruction. As any AO system correction is only partial, it is indispensable to use post-processing reconstruction techniques. In the New Vacuum Solar Telescope (NVST), speckle masking method is used to achieve the diffraction limited resolution of the telescope. Although the method is very promising, the computation is quite intensive, and the amount of data is tremendous, requiring several months to reconstruct observational data of one day on a high-end computer. To accelerate image reconstruction, we parallelize the program package on a high performance cluster. We describe parallel implementation details for several reconstruction procedures. The code is written in C language using Message Passing Interface (MPI) and optimized for parallel processing in a multi-processor environment. We show the excellent performance of parallel implementation, and the whole data processing speed is about 71 times faster than before. Finally, we analyze the scalability of the code to find possible bottlenecks, and propose several ways to further improve the parallel performance. We conclude that the presented program is capable of executing in real-time reconstruction applications at NVST.

GRB 140903A, a short duration $\gamma-$ray burst (SGRB) detected by {\it Swift}, is characterized by its long-lasting radio emission among SGRBs. In addition to the $\sim 10^{6}$ s radio afterglow emission, the afterglow of GRB 140903A displays a plateau from $10^3$ s to $7 \times 10^3 $ s in the X-rays. In this work, we attribute the X-ray plateau to the energy injection into the decelerating blast wave and then model the later radio/optical/X-ray afterglow emission within the standard fireball afterglow model. The afterglow emission has been well reproduced with reasonable physical parameters, including a jet half-opening angle $\sim 0.05$. We also compare the physical properties of SGRBs with and without radio afterglow detections and find that nearby SGRBs are brighter and easier to be detected in radio. Such a phenomena has been interpreted within the afterglow model.

Recently Jin et al. (2016) reanalyzed the optical observation data of GRB 050709 and reported a line-like spectral energy distribution (SED) component observed by the Very Large Telescope at $t\sim 2.5$ days after the trigger of the burst, which had been interpreted as a broadened line signal arising from a macronova dominated by iron group. In this work we show that an optical flare origin of such a peculiar optical SED is still possible. Interestingly, even in such a model, an "unusual" origin of the late-time long-lasting Hubble Space Telescope F814W-band emission is still needed and a macronova/kilonova is the natural interpretation.

Neutron star$-$black hole (NS$-$BH) coalescences are widely believed to be promising gravitational wave sources in the era of advanced detectors of LIGO/Virgo but such binaries have never been directly detected yet. Evidence for NS$-$BH coalescences have been suggested in short and hybrid GRB observations, which are examined critically. Based on the suggested connection between the observed macronovae/kilonovae events and NS$-$BH coalescences, we get a fiducial lower limit of NS$-$BH coalescence rate density ${\cal R}_{\rm nsbh} \approx 18.8^{+12.5}_{-8.6} ~{\rm Gpc^{-3}~ yr^{-1}~ (\theta_j/0.1~{\rm rad})^{-2}}$, where $\theta_{\rm j}$ is the typical half-opening angle of the GRB ejecta. The real value of ${\cal R}_{\rm nsbh}$ is likely at least $\sim {\rm a~few}$ times larger, depending upon the equation of state of NS material and the properties of the NS$-$BH system, such as the mass and spin distribution of the black hole. If the link between macronovae/kilonovae and NS$-$BH coalescence is valid, one can expect that at design sensitivity the aLIGO/AdVirgo network will detect NS$-$BH coalescence signals at a rate of at least a dozen per year, and to consequently place constraints on certain physical properties of NS$-$BH systems.

GRB 050709 was the first short Gamma-ray Burst (sGRB) with an identified optical afterglow. In this work we report a re-analysis of the publicly available data of this event and the discovery of a Li-Paczynski macronova/kilonova that dominates the optical/IR signal at t $>$ 2.5 days. Such a signal would arise from ~0.05 M$_\odot$ r-process material launched by a compact binary merger. The implied mass ejection supports the suggestion that compact binary mergers are significant and possibly main sites of heavy r-process nucleosynthesis. We have re-analyzed all afterglow data from nearby short and hybrid GRBs. A statistical study of sGRB/macronova connection reveals that macronova may have taken place in all these GRBs though the fraction as low as 0.18 cannot be ruled out. The identification of two of the three macronova candidates in the I-band implies a more promising detection prospect for ground-based surveys.

The merger-driven Gamma-ray Bursts (GRBs) and their associated gravitational wave (GW) radiation, if both successfully detected, have some far-reaching implications, including for instance: (i) The statistical comparison of the physical properties of the short/long-short GRBs with and without GW detection can test the general origin model; (ii) Revealing the physical processes taking place at the central engine; (iii) Measuring the velocity of the Gravitational wave directly/accurately. In this work we discuss these implications in the case of possible association of GW150914/ GBM transient 150914. We compared GBM transient 150914 with other SGRBs and found that such an event {may be} a distinct outlier in some statistical diagrams, possibly due to its specific binary-black-hole merger origin. However, the presence of a "new" group of SGRBs with "unusual" physical parameters is also possible. If the outflow of GBM transient 150914 was launched by the accretion onto the nascent black hole, the magnetic activity rather than the neutrino process is likely responsible for the energy extraction and the accretion disk mass is estimated to be $\sim 10^{-5}~M_\odot$. The GW150914/GBM transient 150914 association, {if confirmed, would} provide the first opportunity to directly measure the GW velocity and its departure from the speed of the light {should be within} a factor of $\sim 10^{-17}$.

After the jet break at $t\sim 1.4$ days, the optical afterglow emission of the long-short burst GRB 060614 can be divided into two components. One is the power-law decaying forward shock afterglow emission. The other is an excess of flux in several multi-band photometric observations, which emerges at $\sim$4 days after the burst, significantly earlier than that observed for a supernova associated with a long-duration GRB. At $t>13.6$ days, the F814W-band flux drops faster than $t^{-3.2}$. Moreover, the spectrum of the excess component is very soft and the luminosity is extremely low. These observed signals are incompatible with those from weak supernovae, but the ejection of $\sim 0.1~M_\odot$ of $r-$process material from a black hole-neutron star merger, as recently found in some numerical simulations, can produce it. If this interpretation is correct, it represents the first time that a multi-epoch/band lightcurve of a Li-Paczynski macronova (also known as kilonova) has been obtained and black hole-neutron star mergers are sites of significant production of $r-$process elements.

The {\it Swift}-detected GRB 060614 was a unique burst that straddles an imaginary divide between long- and short-duration gamma-ray bursts (GRBs), and its physical origin has been heavily debated over the years. Recently, a distinct very-soft F814W-band excess at $t\sim 13.6$ days after the burst was identified in a joint-analysis of VLT and HST optical afterglow data of GRB~060614, which has been interpreted as evidence for an accompanying Li-Paczynski macronova (also called a kilonova). Under the assumption that the afterglow data in the time interval of $1.7-3.0$ days after the burst are due to external forward shock emission, when this assumption is extrapolated to later times it is found that there is an excess of flux in several multi-band photometric observations. This component emerges at $\sim$4 days after the burst, and it may represent the first time that a multi-epoch/band lightcurve of a macronova has been obtained. The macronova associated with GRB 060614 peaked at $t\lesssim 4$ days after the burst, which is significantly earlier than that observed for a supernova associated with a long-duration GRB. Due to the limited data, no strong evidence for a temperature evolution is found. We derive a conservative estimate of the macronova rate of $\sim 16.3^{+16.3}_{-8.2}~{\rm Gpc^{-3}}{\rm yr^{-1}}$, implying a promising prospect for detecting the gravitational wave radiation from compact object mergers by upcoming Advanced LIGO/VIRGO/KAGRA detectors (i.e., the rate is ${\cal R}_{\rm GW} \sim 0.5^{+0.5}_{-0.25}(D/200~{\rm Mpc})^{3}~{\rm yr^{-1}}$).

Long-duration ($>2$ s) $\gamma$-ray bursts that are believed to originate from the death of massive stars are expected to be accompanied by supernovae. GRB 060614, that lasted 102 s, lacks a supernova-like emission down to very stringent limits and its physical origin is still debated. Here we report the discovery of near-infrared bump that is significantly above the regular decaying afterglow. This red bump is inconsistent with even the weakest known supernova. However, it can arise from a Li-Paczy\'{n}ski macronova $-$ the radioactive decay of debris following a compact binary merger. If this interpretation is correct GRB 060614 arose from a compact binary merger rather than from the death of a massive star and it was a site of a significant production of heavy r-process elements. The significant ejected mass favors a black hole-neutron star merger but a double neutron star merger cannot be ruled out.

In this paper, we report the detections of stellar variabilities from the first 2-year observations of sky area of about 1300 square degrees from the Tsinghua University-NAOC Transient Survey (TNTS). A total of 1237 variable stars (including 299 new ones) were detected with brightness < 18.0 mag and magnitude variation >= 0.1 mag on a timescale from a few hours to few hundred days. Among such detections, we tentatively identified 661 RR Lyrae stars, 431 binaries, 72 Semiregular pulsators, 29 Mira stars, 11 slow irregular variables, 11 RS Canum Venaticorum stars, 7 Gamma Doradus stars, 5 long period variables, 3 W Virginis stars, 3 Delta Scuti stars, 2 Anomalous Cepheids, 1 Cepheid, and 1 nove-like star based on their time-series variability index Js and their phased diagrams. Moreover, we found that 14 RR Lyrae stars show the Blazhko effect and 67 contact eclipsing binaries exhibit the O'Connell effect. Since the period and amplitude of light variations of RR Lyrae variables depend on their chemical compositions, their photometric observations can be used to investigate distribution of metallicity along the direction perpendicular to the Galactic disk. We find that the metallicity of RR Lyrae stars shows large scatter at regions closer to the Galactic plane (e.g., -3.0 < [Fe/H] < 0) but tends to converge at [Fe/H]~ -1.7 at larger Galactic latitudes. This variation may be related to that the RRAB Lyrae stars in the Galactic halo come from globular clusters with different metallicity and vertical distances, i.e. OoI and OoII populations, favoring for the dual-halo model.

The short burst GRB 130912A was detected by Swift, Fermi satellites and several ground-based optical telescopes. Its X-ray light curve decayed with time normally. The optical emission, however, displayed a long term plateau, which is the longest one in current short GRB observations. In this work we examine the physical origin of the X-ray and optical emission of this peculiar event. We find that the canonical forward shock afterglow emission model can account for the X-ray and optical data self-consistently and the energy injection model that has been widely adopted to interpret the shallowly-decaying afterglow emission is not needed. We also find that the burst was born in a very-low density interstellar medium, consistent with the compact object merger model. Significant fractions of the energy of the forward shock have been given to accelerate the non-thermal electrons and amplify the magnetic fields (i.e., $\epsilon_{\rm e}\sim 0.37$ and $\epsilon_{\rm B}\sim 0.16$, respectively), which are much larger than those inferred in most short burst afterglow modeling and can explain why the long-lasting optical afterglow plateau is rare in short GRBs.

Gamma-ray bursts (GRBs) are most probably powered by collimated relativistic outflows (jets) from accreting black holes at cosmological distances. Bright afterglows are produced when the outflow collides with the ambient medium. Afterglow polarization directly probes the magnetic properties of the jet, when measured minutes after the burst, and the geometric properties of the jet and the ambient medium when measured hours to days after the burst. High values of optical polarization detected minutes after burst in GRB 120308A indicate the presence of large-scale ordered magnetic fields originating from the central engine (the power source of the GRB). Theoretical models predict low degrees of linear polarization and negligable circular polarization at late times, when the energy in the original ejecta is quickly transferred to the ambient medium and propagates farther into the medium as a blastwave. Here we report the detection of circularly polarized optical light in the afterglow of GRB 121024A, measured 0.15 days after the burst. We show that the circular polarization is intrinsic to the afterglow and unlikely to be produced by dust scattering or plasma propagation effects. A possible explanation is to invoke anisotropic (rather than the commonly assumed isotropic) electron pitch angle distributions, and we suggest that new models are required to produce the complex microphysics of realistic shocks in relativistic jets.

GRB 120308A, a long duration $\gamma-$ray burst detected by {\it Swift}, was distinguished by a highly-polarized early optical afterglow emission that strongly suggests an ordered magnetic field component in the emitting region. In this work we model the optical and X-ray emission in the reverse and forward shock scenario and show that the strength of the magnetic field in reverse shock region is $\sim 10$ times stronger than that in the forward shock region. Consequently the outflow powering the highly-polarized reverse shock optical emission was mildly-magnetized at a degree $\sigma \sim$ a few percent. Considering the plausible magnetic energy dissipation in both the acceleration and prompt emission phases of the Gamma-ray Burst (GRB) outflow, the afterglow data of GRB 120308A provides us the compelling evidence that at least for some GRBs a non-ignorable fraction of the energy was released in the form of Poynting-flux, confirming the finding firstly made in the reverse-forward shock emission modeling of the optical afterglow of GRB 990123 (Fan et al. 2002; Zhang et al. 2003).

GRB 131231A was detected by the Large Area Telescope onboard Fermi Space Gamma-ray Telescope. The high energy gamma-ray ($> 100$ MeV) afterglow emission spectrum is $F_\nu \propto \nu^{-0.54\pm0.15}$ in the first $\sim 1300$ s after the trigger and the most energetic photon has an energy $\sim 62$ GeV arriving at $t\sim 520$ s. With reasonable parameters of the GRB outflow as well as the density of the circum-burst medium, the synchrotron radiation of electrons or protons accelerated at an external forward shock have difficulty accounting for the data. The synchrotron self-Compton radiation of the forward shock-accelerated electrons, instead, can account for both the spectrum and temporal behavior of the GeV afterglow emission. We also show that the prospect for detecting GRB 131231A$-$like GRBs with Cherenkov Telescope Array (CTA) is promising.

The New Vacuum Solar Telescope (NVST) is a 1 meter vacuum solar telescope that aims to observe the fine structures on the Sun. The main tasks of NVST are high resolution imaging and spectral observations, including the measurements of solar magnetic field. NVST is the primary ground-based facility of Chinese solar community in this solar cycle. It is located by the Fuxian Lake of southwest China, where the seeing is good enough to perform high resolution observations. In this paper, we first introduce the general conditions of Fuxian Solar Observatory and the primary science cases of NVST. Then, the basic structures of this telescope and instruments are described in detail. Finally, some typical high resolution data of solar photosphere and chromosphere are also shown.

In this paper we give a brief review of our recent studies on the long and short gamma-ray bursts (GRBs) detected Swift, in an effort to understand the puzzle of classifying GRBs. We consider that it is still an appealing conjecture that both long and short GRBs are drawn from the same parent sample by observational biases.

We show that the peculiar early optical and in particular X-ray afterglow emission of the short duration burst GRB 130603B can be explained by continuous energy injection into the blastwave from a supra-massive magnetar central engine. The observed energetics and temporal/spectral properties of the late infrared bump (i.e., the "kilonova") are also found consistent with emission from the ejecta launched during an NS-NS merger and powered by a magnetar central engine. The isotropic-equivalent kinetic energies of both the GRB blastwave and the kilonova are about $E_{\rm k}\sim 10^{51}$ erg, consistent with being powered by a near-isotropic magnetar wind. However, this relatively small value demands that most of the initial rotational energy of the magnetar $(\sim {\rm a~ few \times 10^{52}~ erg})$ is carried away by gravitational wave radiation. Our results suggest that (i) the progenitor of GRB 130603B would be a NS-NS binary system, whose merger product would be a supra-massive neutron star that lasted for about $\sim 1000$ seconds; (ii) the equation-of-state of nuclear matter would be stiff enough to allow survival of a long-lived supra-massive neutron star, so that it is promising to detect bright electromagnetic counterparts of gravitational wave triggers without short GRB associations in the upcoming Advanced LIGO/Virgo era.

Short-hard gamma-ray bursts (GRBs) are widely believed to be produced by the merger of two binary compact objects, specifically by two neutron stars or by a neutron star orbiting a black hole. According to the Li-Paczynski kilonova model, the merger would launch sub-relativistic ejecta and a near-infrared/optical transient would then occur, lasting up to days, which is powered by the radioactive decay of heavy elements synthesized in the ejecta. The detection of a late bump using the {\em Hubble Space Telescope} ({\em HST}) in the near-infrared afterglow light curve of the short-hard GRB 130603B is indeed consistent with such a model. However, as shown in this Letter, the limited {\em HST} near-infrared lightcurve behavior can also be interpreted as the synchrotron radiation of the external shock driven by a wide mildly relativistic outflow. In such a scenario, the radio emission is expected to peak with a flux of $\sim 100 \mu$Jy, which is detectable for current radio arrays. Hence, the radio afterglow data can provide complementary evidence on the nature of the bump in GRB 130603B. It is worth noting that good spectroscopy during the bump phase in short-hard bursts can test validity of either model above, analogous to spectroscopy of broad-lined Type Ic supernova in long-soft GRBs.

We discuss the results of the analysis of multi-wavelength data for the afterglows of GRB 081007 and GRB 090424, two bursts detected by Swift. One of them, GRB 081007, also shows a spectroscopically confirmed supernova, SN 2008hw, which resembles SN 1998bw in its absorption features, while the maximum luminosity is only about half as large as that of SN 1998bw. Bright optical flashes have been detected in both events, which allows us to derive solid constraints on the circumburst-matter density profile. This is particularly interesting in the case of GRB 081007, whose afterglow is found to be propagating into a constant-density medium, yielding yet another example of a GRB clearly associated with a massive star progenitor which did not sculpt the surroundings with its stellar wind. There is no supernova component detected in the afterglow of GRB 090424, likely due to the brightness of the host galaxy, comparable to the Milky Way. We show that the afterglow data are consistent with the presence of both forward- and reverse-shock emission powered by relativistic outflows expanding into the interstellar medium. The absence of optical peaks due to the forward shock strongly suggests that the reverse shock regions should be mildly magnetized. The initial Lorentz factor of outflow of GRB 081007 is estimated to be \Gamma ~ 200, while for GRB 090424 a lower limit of \Gamma > 170 is derived. We also discuss the prompt emission of GRB 081007, which consists of just a single pulse. We argue that neither the external forward-shock model nor the shock-breakout model can account for the prompt emission data and suggest that the single-pulse-like prompt emission may be due to magnetic energy dissipation of a Poynting-flux dominated outflow or to a dissipative photosphere.

A nearby super-luminous burst GRB 130427A was simultaneously detected by six $\gamma$-ray space telescopes ({\it Swift}, Fermi-GBM/LAT, Konus-Wind, SPI-ACS/INTEGRAL, AGILE and RHESSI) and by three RAPTOR full-sky persistent monitors. The isotropic $\gamma-$ray energy release is of $\sim 10^{54}$ erg, rendering it the most powerful explosion among the GRBs with a redshift $z\leq 0.5$. The emission above 100 MeV lasted about one day and four photons are at energies greater than 40 GeV. We show that the count rate of 100 MeV-100 GeV emission may be mainly accounted for by the forward shock synchrotron radiation and the inverse Compton radiation likely dominates at GeV-TeV energies. In particular, an inverse Compton radiation origin is established for the $\sim (95.3,~47.3,~41.4,~38.5,~32)$ GeV photons arriving at $t\sim (243,~256.3,~610.6,~3409.8,~34366.2)$ s after the trigger of Fermi-GBM. Interestingly, the external-inverse-Compton-scattering of the prompt emission (the second episode, i.e., $t\sim 120-260$ s) by the forward-shock-accelerated electrons is expected to produce a few $\gamma-$rays at energies above 10 GeV, while five were detected in the same time interval. A possible unified model for the prompt soft $\gamma-$ray, optical and GeV emission of GRB 130427A, GRB 080319B and GRB 090902B is outlined. Implication of the null detection of $>1$ TeV neutrinos from GRB 130427A by IceCube is discussed.

In this paper we compute rest-frame extinctions for the afterglows of a sample of gamma-ray bursts complete in redshift. The selection criteria of the sample are based on observational high-energy parameters of the prompt emission and therefore our sample should not be biased against dusty sight-lines. It is therefore expected that our inferences hold for the general population of gamma-ray bursts. Our main result is that the optical/near-infrared extinction of gamma-ray burst afterglows in our sample does not follow a single distribution. 87% of the events are absorbed by less than 2 mag, and 50% suffer from less than 0.3-0.4 mag extinction. The remaining 13% of the afterglows are highly absorbed. The true percentage of gamma-ray burst afterglows showing high absorption could be even higher since a fair fraction of the events without reliable redshift measurement are probably part of this class. These events may be due to highly dusty molecular clouds/star forming regions associated with the gamma-ray burst progenitor or along the afterglow line of sight, and/or to massive dusty host galaxies. No clear evolution in the dust extinction properties is evident within the redshift range of our sample, although the largest extinctions are at z~1.5-2, close to the expected peak of the star formation rate. Those events classified as dark are characterized, on average, by a higher extinction than typical events in the sample. A correlation between optical/near-infrared extinction and hydrogen-equivalent column density based on X-ray studies is shown although the observed NH appears to be well in excess compared to those observed in the Local Group. Dust extinction does not seem to correlate with GRB energetics or luminosity.

We present an analysis of the unusual optical light curve of the gamma-ray burst GRB 081029, a long-soft burst with a redshift of z = 3.8479. We combine X-ray and optical observations from the Swift X-Ray Telescope and the Swift UltraViolet/Optical Telescope with ground-based optical and infrared data obtained using the REM, ROTSE, and CTIO 1.3-m telescopes to construct a detailed data set extending from 86 s to approximately 100,000 s after the BAT trigger. Our data cover a wide energy range, from 10 keV to 0.77 eV (1.24 to 16,000 Angstrom). The X-ray afterglow shows a shallow initial decay followed by a rapid decay starting at about 18,000 s. The optical and infrared afterglow, however, shows an uncharacteristic rise at about 3000 s that does not correspond to any feature in the X-ray light curve. Our data are not consistent with synchrotron radiation from a jet interacting with an external medium, a two-component jet, or continuous energy injection from the central engine. We find that the optical light curves can be broadly explained by a collision between two ejecta shells within a two-component jet. A growing number of gamma-ray burst afterglows are consistent with complex jets, which suggests that some (or all) gamma-ray burst jets are complex and will require detailed modelling to fully understand them.

In order to study the effect of dust extinction on the afterglow of gamma-ray bursts (GRBs), we carry out numerical calculations with high precision based on rigorous Mie theory and latest optical properties of interstellar dust grains, and analyze the different extinction curves produced by dust grains with different physical parameters. Our results indicate that the absolute extinction quantity is substantially determined by the medium density and metallicity. However, the shape of the extinction curve is mainly determined by the size distribution of the dust grains. If the dust grains aggregate to form larger ones, they will cause a flatter or grayer extinction curve with lower extinction quantity. On the contrary, if the dust grains are disassociated to smaller ones due to some uncertain processes, they will cause a steeper extinction curve with larger amount of extinction. These results might provide an important insight into understanding the origin of the optically dark GRBs.

In an effort to understand the puzzle of classifying gamma-ray bursts (GRBs), we perform a systematic study of {\it Swift} GRBs and investigate several short GRB issues. Though short GRBs have a short ($\lesssim2$ s) prompt duration as monitored by the Burst Alert Telescope, the composite light curves including both the prompt and afterglow emission suggest that most of the short GRBs have a similar radiative feature to long GRBs. Further, some well-studied short GRBs might also have an intrinsically long prompt duration, which renders them as a type of short GRB imposters. Genuine short GRBs detected by {\it Swift} might be rare that discriminating the observed short GRBs is, not surprisingly, troublesome. In particular, the observational biases in the host identification and redshift measurement of GRBs should be taken with great caution. The redshift distribution which has been found to be different for long and short GRBs might have been strongly affected by the measurement methods. We find that the redshifts measured from the presumed host galaxies of long and short GRBs appear to have a similar distribution.

We develop two methods to estimate the bulk Lorentz factor of X-ray flare outflow. In the first method the outflow is assumed to be baryonic and is accelerated by the thermal pressure, for which the final bulk Lorentz factor is limited by the outflow luminosity as well as the initial radius of the outflow getting accelerated. Such a method may be applied to a considerable fraction of flares. The second method, based on the curvature effect interpretation of the quick decline of the flare, can give a tightly constrained estimate of the bulk Lorentz factor but can only be applied to a few giant flares. The results obtained in these two different ways are consistent with each other. The obtained bulk Lorentz factor (or just upper limit) of the X-ray flare outflows, ranging from ten to a few hundred, is generally smaller than that of the Gamma-ray Burst (GRB) outflows.

VisIVO is an integrated suite of tools and services specifically designed for the Virtual Observatory. This suite constitutes a software framework for effective visual discovery in currently available (and next-generation) very large-scale astrophysical datasets. VisIVO consists of VisiVO Desktop - a stand alone application for interactive visualization on standard PCs, VisIVO Server - a grid-enabled platform for high performance visualization and VisIVO Web - a custom designed web portal supporting services based on the VisIVO Server functionality. The main characteristic of VisIVO is support for high-performance, multidimensional visualization of very large-scale astrophysical datasets. Users can obtain meaningful visualizations rapidly while preserving full and intuitive control of the relevant visualization parameters. This paper focuses on newly developed integrated tools in VisIVO Server allowing intuitive visual discovery with 3D views being created from data tables. VisIVO Server can be installed easily on any web server with a database repository. We discuss briefly aspects of our implementation of VisiVO Server on a computational grid and also outline the functionality of the services offered by VisIVO Web. Finally we conclude with a summary of our work and pointers to future developments.

The scientific community is presently witnessing an unprecedented growth in the quality and quantity of data sets coming from simulations and real-world experiments. To access effectively and extract the scientific content of such large-scale data sets (often sizes are measured in hundreds or even millions of Gigabytes) appropriate tools are needed. Visual data exploration and discovery is a robust approach for rapidly and intuitively inspecting large-scale data sets, e.g. for identifying new features and patterns or isolating small regions of interest within which to apply time-consuming algorithms. This paper presents a high performance parallelized implementation of Splotch, our previously developed visual data exploration and discovery algorithm for large-scale astrophysical data sets coming from particle-based simulations. Splotch has been improved in order to exploit modern massively parallel architectures, e.g. multicore CPUs and CUDA-enabled GPUs. We present performance and scalability benchmarks on a number of test cases, demonstrating the ability of our high performance parallelized Splotch to handle efficiently large-scale data sets, such as the outputs of the Millennium II simulation, the largest cosmological simulation ever performed.

We present the prompt BAT and afterglow XRT data of Swift-discovered GRB081109A up to ~ 5\times 10^5 sec after the trigger, and the early ground-based optical follow-ups. The temporal and spectral indices of the X-ray afterglow emission change remarkably. We interpret this as the GRB jet first traversing the freely expanding supersonic stellar wind of the progenitor with density varying as $\rho \propto r^{-2}$. Then after approximately 300 sec the jet traverses into a region of apparent constant density similar to that expected in the stalled-wind region of a stellar wind bubble or the interstellar medium (ISM). The optical afterglow data are generally consistent with such a scenario. Our best numerical model has a wind density parameter {$A_{*} \sim 0.02$, a density of the stalled wind $n\sim 0.12 {\rm cm}^{-3}$, and a transition radius $ \sim 4.5 \times 10^{17}$ cm}. Such a transition radius is smaller than that predicted by numerical simulations of the stellar wind bubbles and may be due to a rapidly evolving wind of the progenitor close to the time of its core-collapse.

Rich soft X-ray emission lines of highly charged silicon ions (Si VI--Si XII) were observed by irradiating an ultra-intense laser pulse with width of 200 fs and energy of $\sim$90 mJ on the solid silicon target. The high resolution spectra of highly charged silicon ions with full-width at half maximum (FWHM) of $\sim$0.3--0.4\AA is analyzed in wavelength range of 40--90 \AA . The wavelengths of 53 prominent lines are determined with statistical uncertainties being up to 0.005 \AA . Collisional-radiative models were constructed for Si VI -- Si XII ions, which satisfactorily reproduces the experimental spectra, and helps the line identification. Calculations at different electron densities reveal that the spectra of dense plasmas are more complicate than the spectra of thin plasmas. A comparison with the Kelly database reveals a good agreement for most peak intensities, and differences for a few emission lines.

We constrain the circum-burst medium profile with the rise behavior of the very early afterglow light curves of gamma-ray bursts (GRBs). Using this method, we find a constant and low-density medium profile for GRB 060418 and GRB 060607A, which is consistent with the inference from the late afterglow data. In addition, we show that the absence of the IR flashes in these two GRB afterglows is consistent with the standard external reverse shock model, which thus renders models like the highly magnetized GRB outflow being unnecessary to a certain extent.

In the double neutron star merger or neutron star-black hole merger model for short GRBs, the outflow launched might be mildly magnetized and neutron rich. The magnetized neutron-rich outflow will be accelerated by the magnetic and thermal pressure and may form a two component jet finally, as suggested by Vlahakis, Peng & K\"{o}nigl (2003). We show in this work that such a two component jet model could well reproduce the multi-wavelength afterglow lightcurves, in particular the X-ray flat segment, of short GRB 051221A. In this model, the central engine need not to be active much longer than the prompt $\gamma-$ray emission.

It is widely believed that outflows of gamma-ray bursts are jetted. Some people also suggest that the jets may have structures like $\epsilon(\theta)\propto\theta^{-k}$. We test the possibility of X-ray flashes coming from such jets in this paper. Both qualitative and quantitative analyses have shown that this model can reproduce most observational features for both X-ray flashes and gamma-ray bursts. Using common parameters of gamma-ray bursts, we have done both uniform and nonuniform jets' numerical calculations for their fluxes, spectra and peak energies. It seems that nonuniform jets are more appropriated to these observational properties than uniform ones. We also give a spectrum and flux fit to the most significant X-ray flash, XRF971019 by our model. We also have shown that in our model the observational ratio of gamma-ray bursts to X-ray flashes is about several.

Here we calculate the GRB afterglow light curves from a relativistic jet as seen by observers at a wide range of viewing angles from the jet axis, and the jet is uniform or non-uniform. We find that, for uniform jet the afterglow light curves for different viewing angles are somewhat different: in general, there are two breaks in the light curve, corresponding to the time $\gamma\sim (\theta_j-\theta_v)^{-1}$ and $\gamma\sim (\theta_j+\theta_v)^{-1}$ respectively. However, for non-uniform jet, the things become more complicated. For the case $\theta_v=0$, we can obtain the analytical results, for $k<8/(p+4)$ there should be two breaks in the light curve correspond to $\gamma\sim\theta_c^{-1}$ and $\gamma\sim\theta_j^{-1}$ respectively, while for $k>8/(p+4)$ there should be only one break corresponds to $\gamma\sim\theta_c^{-1}$, and this provides a possible explanation for some rapidly fading afterglows whose light curves have no breaks since the time at which $\gamma\sim\theta_c^{-1}$ is much earlier than our first observation time. For the case $\theta_v\neq 0$, our numerical results show that, the afterglow light curves are strongly affected by the values of $\theta_v$, $\theta_c$ and $k$. If the values of $\theta_v/\theta_c$ and $k$ are larger, there will be a prominent flattening in the afterglow light curve, which is quite different from the uniform jet, and after the flattening a very sharp break will be occurred at the time $\gamma\sim (\theta_v + \theta_c)^{-1}.