Gamma-ray burst (GRB) 150910A was detected by {\it Swift}/BAT, and then rapidly observed by {\it Swift}/XRT, {\it Swift}/UVOT, and ground-based telescopes. We report Lick Observatory spectroscopic and photometric observations of GRB~150910A, and we investigate the physical origins of both the optical and X-ray afterglows, incorporating data obtained with BAT and XRT. The light curves show that the jet emission episode lasts $\sim 360$~s with a sharp pulse from BAT to XRT (Episode I). In Episode II, the optical emission has a smooth onset bump followed by a normal decay ($\alpha_{\rm R,2} \approx -1.36$), as predicted in the standard external shock model, while the X-ray emission exhibits a plateau ($\alpha_{\rm X,1} \approx -0.36$) followed by a steep decay ($\alpha_{\rm X,2} \approx -2.12$). The light curves show obvious chromatic behavior with an excess in the X-ray flux. Our results suggest that GRB 150910A is an unusual GRB driven by a newly-born magnetar with its extremely energetic magnetic dipole (MD) wind in Episode II, which overwhelmingly dominates the observed early X-ray plateau. The radiative efficiency of the jet prompt emission is $\eta_{\gamma} \approx 11\%$. The MD wind emission was detected in both the BAT and XRT bands, making it the brightest among the current sample of MD winds seen by XRT. We infer the initial spin period ($P_0$) and the surface polar cap magnetic field strength ($B_p$) of the magnetar as $1.02 \times 10^{15}~{\rm G} \leq B_{p} \leq 1.80 \times 10^{15}~{\rm G}$ and 1~ms $\leq P_{0}v\leq 1.77$~ms, and the radiative efficiency of the wind is $\eta_w \geq 32\%$.

Ultra-long-duration gamma-ray burst GRB 111209A was found to be associated with a very luminous supernovae (SNe) SN 2011kl. The physics of GRB 111209A/SN 2011kl has been extensively studied in the literures, but does not settle down yet. By investigating in detail the characteristics of the X-ray light curve of GRB 111209A, coupled with the temporal and spectral features observed in SN 2011kl, we argue that a short-living supramassive magnetar can be responsible for the initial shallow X-ray emission. Then the electromagnetic extraction of spin energy from a black hole results in the steeply declining X-ray flux when the magnetar collapses into a black hole (BH). A fraction of the envelope materials falls back and activates the accretion onto the newborn BH, which produces the X-ray rebrightening bump at late times. During this process, a centrifugally driven baryon-rich quasi-isotropic Blandford \& Payne outflow from the revived accretion disk deposits its kinetic energy on the SN ejecta, which powers luminous SN 2011kl. Finally, we place a limitation on the magnetar's physical parameters based on the observations.

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

Violent multi-wavelength variabilities are observed in gamma-ray-selected blazars. We present an analysis of long-term light curves for eight bright blazars to explore the co-variation pattern in the gamma-ray and radio bands. We extract their gamma-ray light curves and spectra with data observed by the Fermi/LAT since 2008. We find diverse co-variation patterns between the gamma-ray and radio (at 43 GHz) fluxes in these sources. The gamma-ray and radio fluxes of 3C 454.3 and PKS 1633+382 are correlated without any time-lag, suggesting that they are from the same radiation region. Similar correlation is also observed in 3C 273 and PKS 1222+216, but the radio flux is lag behind the gamma-ray flux approximately ~160 days and ~290 days, respectively. This likely suggests that their gamma-ray emission regions are located at the upstream of their radio cores at 43 GHz. The gamma-ray and radio fluxes of the other four blazars are not correlated, implying that the gamma-ray and radio emission may be from different regions in their jets. The gamma-ray light curves of the eight blazars can be decomposed into some long timescale variability components and fast spike flares. We propose that they may be attributed to the central engine activity and the magnetic reconnection process or turbulence in the local emission region, respectively.

Short GRBs with extended emission (EE) that are composed initially of a short-hard spike and followed by a long-lasting EE, are thought to be classified as a subsection of short GRBs. The narrow energy band available during the {\em Swift} era combined with a lack of spectral information prevented discovery of the intrinsic properties of those events. In this paper, we performed a systematic search of short GRBs with EE by using all available {\em Fermi}/GBM data. The search identified 26 GBM-detected short GRBs with EE that are similar to GRB 060614 observed by {\em Swift}/BAT. We focus on investigating the spectral and temporal properties for both the hard spike and the EE components of all 26 GRBs, and explore differences and possible correlations between them. We find that while the peak energy ($E_{\rm p}$) of the hard spikes is a little bit harder than that of the EE, but their fluences are comparable. The harder $E_{\rm p}$ seems to correspond to a larger fluence and peak flux with a large scatter for both the hard spike and EE components. Moreover, the $E_{\rm p}$ of both the hard spikes and EE are compared to other short GRBs. Finally, we also compare the properties of GRB 170817A with those short GRBs with EE and find no significant statistical differences between them. We find that GRB 170817A has the lowest $E_{\rm p}$, likely because it was off-axis.

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

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

The Lorentz factor of a relativistic jet and its evolution during the jet expansion are difficult to estimate, especially for the jets in gamma-ray bursts (GRBs). However, it is related to the understanding of jet physics. Owing to the absorption of two-photon pair production ($\gamma\gamma{\leftrightarrow}e^+e^-$), a high-energy spectral cutoff may appear in the radiation spectrum of GRBs. We search such kind of high-energy cutoff in GRB 160625B, which is one of the brightest bursts in recent years. It is found that the high-energy spectral cutoff is obvious for the first pulse in the second emission episode of GRB 160625B (i.e., $\sim$186-192 s after the burst first trigger), which is smooth and well-shaped. Then, we estimate the Lorentz factor and radiation location of the jet shell associated with the first pulse in the second emission episode of GRB 160625B. It is found that the radiation location increases with time. In addition, the Lorentz factor remains almost constant during the expansion of the jet shell. This reveals that the magnetization of the jet is low or intermediate in the emission region, event though the jet could be still Poynting flux dominated at smaller radii to avoid a bright thermal component in the emission episode.

The Galactic Center supermassive black hole Sagittarius A* (Sgr A*) is one of the most promising targets to study the dynamics of black hole accretion and outflow via direct imaging with very long baseline interferometry (VLBI). At 3.5 mm (86 GHz), the emission from Sgr A* is resolvable with the Global Millimeter VLBI Array (GMVA). We present the first observations of Sgr A* with the phased Atacama Large Millimeter/submillimeter Array (ALMA) joining the GMVA. Our observations achieve an angular resolution of ~87{\mu}as, improving upon previous experiments by a factor of two. We reconstruct a first image of the unscattered source structure of Sgr A* at 3.5 mm, mitigating effects of interstellar scattering. The unscattered source has a major axis size of 120 $\pm$ 34{\mu}as (12 $\pm$ 3.4 Schwarzschild radii), and a symmetrical morphology (axial ratio of 1.2$^{+0.3}_{-0.2}$), which is further supported by closure phases consistent with zero within 3{\sigma}. We show that multiple disk-dominated models of Sgr A* match our observational constraints, while the two jet-dominated models considered are constrained to small viewing angles. Our long-baseline detections to ALMA also provide new constraints on the scattering of Sgr A*, and we show that refractive scattering effects are likely to be weak for images of Sgr A* at 1.3 mm with the Event Horizon Telescope. Our results provide the most stringent constraints to date for the intrinsic morphology and refractive scattering of Sgr A*, demonstrating the exceptional contribution of ALMA to millimeter VLBI.

We report results from a deep polarization imaging of the nearby radio galaxy 3C$\,$84 (NGC$\,$1275). The source was observed with the Global Millimeter VLBI Array (GMVA) at 86$\,$GHz at an ultra-high angular resolution of $50\mu$as (corresponding to 250$R_{s}$). We also add complementary multi-wavelength data from the Very Long Baseline Array (VLBA; 15 & 43$\,$GHz) and from the Atacama Large Millimeter/submillimeter Array (ALMA; 97.5, 233.0, and 343.5$\,$GHz). At 86$\,$GHz, we measure a fractional linear polarization of $\sim2$% in the VLBI core region. The polarization morphology suggests that the emission is associated with an underlying limb-brightened jet. The fractional linear polarization is lower at 43 and 15$\,$GHz ($\sim0.3-0.7$% and $<0.1$%, respectively). This suggests an increasing linear polarization degree towards shorter wavelengths on VLBI scales. We also obtain a large rotation measure (RM) of $\sim10^{5-6}~{\rm rad/m^{2}}$ in the core at $\gtrsim$43$\,$GHz. Moreover, the VLBA 43$\,$GHz observations show a variable RM in the VLBI core region during a small flare in 2015. Faraday depolarization and Faraday conversion in an inhomogeneous and mildly relativistic plasma could explain the observed linear polarization characteristics and the previously measured frequency dependence of the circular polarization. Our Faraday depolarization modeling suggests that the RM most likely originates from an external screen with a highly uniform RM distribution. To explain the large RM value, the uniform RM distribution, and the RM variability, we suggest that the Faraday rotation is caused by a boundary layer in a transversely stratified jet. Based on the RM and the synchrotron spectrum of the core, we provide an estimate for the magnetic field strength and the electron density of the jet plasma.

4C +49.22 is a gamma-ray flat spectrum radio quasar with a bright and knotty jet. We investigate the properties of the core and large-scale knots by using their spectral energy distributions (SEDs). Analyzing its Fermi/LAT data in the past 8 years, a long-term steady gamma-ray emission component is found besides bright outbursts. For the core region, the gamma-ray emission together with the simultaneous emission in the low-energy bands at different epochs is explained with the single-zone leptonic model. The derived magnetization parameters and radiation efficiencies of the radio-core jet decrease as gamma-ray flux decays, likely indicating that a large part of the magnetic energy is converted to the kinetic energy of particles in pc-scale. For the large-scale knots, their radio-optical-X-ray SEDs can be reproduced with the leptonic model by considering the inverse Compton scattering of cosmic microwave background photons. The sum of the predicted gamma-ray fluxes of these knots is comparable to that observed with LAT at 10^{24} Hz of the steady gamma-ray component, indicating that the steady gamma-ray emission may be partially contributed by these large-scale knots. This may conceal the flux variations of the low-level gamma-ray emission from the radio-core. The derived bulk Lorentz factors of the knots decrease along the distance to the core, illustrating as deceleration of jet in large-scale. The powers of the core and knots are roughly in the same order, but the jet changes from highly magnetized at the core region into particle-dominated at the large-scale knots.

The radiation mechanism of the prompt emission of gamma-ray bursts (GRBs) remains an open question. Although their spectra are usually well fitted with the empirical Band function, which is widely believed to be fully non-thermal and interpreted as an optically thin synchrotron emission, accumulating evidence shows that a thermal component actually exists. In this paper, a multi-color blackbody (mBB) model is proposed for the time-integrated spectrum of GRB 081221 by assuming a power-law distribution of the thermal luminosities with temperature, which manifests photospheric emissions from a different radius and/or angle. The effects of the minimum temperature $kT_{\min }$, the maximum temperature $kT_{\max }$ and the power law index $m$ of the luminosity distribution of an mBB are discussed. The fitting to the time-integrated spectrum during the bright phase (from 20s to 30s since the trigger) of GRB 081221 by the mBB model yields $kT_{\min }$ = $4.4\pm 0.3$ keV, $kT_{\max }$ = 57.0 $_{-1.4}^{+1.6}$ keV, and $m=-0.46_{-0.06}^{+0.05}$. When the time bin is small enough, the time-resolved spectra of GRB 081221 are well fitted with a series of single-temperature blackbodies. Our results imply the prompt emission of GRB 081221 is dominated by the photosphere emission and its time-integrated spectrum is a superposition of pure blackbody components at different times, indicating that some empirical Band spectra may be interpreted as mBB if the temperature is widely distributed.

The spectral evolution and spectral lag behavior of 92 bright pulses from 84 gamma-ray bursts (GRBs) observed by the Fermi GBM telescope are studied. These pulses can be classified into hard-to-soft pulses (H2S, 64/92), H2S-dominated-tracking pulses (21/92), and other tracking pulses (7/92). We focus on the relationship between spectral evolution and spectral lags of H2S and H2S-dominated-tracking pulses. %in hard-to-soft pulses (H2S, 64/92) and H2S-dominating-tracking (21/92) pulses. The main trend of spectral evolution (lag behavior) is estimated with $\log E_p\propto k_E\log(t+t_0)$ (${\hat{\tau}} \propto k_{\hat{\tau}}\log E$), where $E_p$ is the peak photon energy in the radiation spectrum, $t+t_0$ is the observer time relative to the beginning of pulse $-t_0$, and ${\hat{\tau}}$ is the spectral lag of photons with energy $E$ with respect to the energy band $8$-$25$ keV. For H2S and H2S-dominated-tracking pulses, a weak correlation between $k_{{\hat{\tau}}}/W$ and $k_E$ is found, where $W$ is the pulse width. We also study the spectral lag behavior with peak time $t_{\rm p_E}$ of pulses for 30 well-shaped pulses and estimate the main trend of the spectral lag behavior with $\log t_{\rm p_E}\propto k_{t_p}\log E$. It is found that $k_{t_p}$ is correlated with $k_E$. We perform simulations under a phenomenological model of spectral evolution, and find that these correlations are reproduced. We then conclude that spectral lags are closely related to spectral evolution within the pulse. The most natural explanation of these observations is that the emission is from the electrons in the same fluid unit at an emission site moving away from the central engine, as expected in the models invoking magnetic dissipation in a moderately-high-$\sigma$ outflow.

Two-episode emission components separated by quiescent gaps in the prompt emission of gamma-ray bursts (GRBs) have been observed in the {\em Swift} era, but there is a lack of spectral information due to the narrow energy band of the {\em Swift}/Burst Alert Telescope. In this paper, a systematic analysis of the spectral and temporal properties of the prompt emission of 101 {\em Fermi}/Gamma-ray Burst Monitor detected long GRBs show the existence of two-episode emission components in the light curves, with quiescent times of up to hundreds of seconds. We focus on investigating the differences of those two emission episodes. We find that the light curves of the two emission components exhibit different behavior, e.g., a soft emission component that either precedes or follows the main prompt emission or that the intensity of the two emission episodes are comparable with each other. No statistically significant correlation in the duration of the two emission episodes can be claimed. We define a new parameter $\varepsilon$ as the ratio of the peak flux of the first and second emission episodes and find that a higher $\varepsilon$ corresponds to a larger fluence. The preferred spectral model in our analysis is a cutoff power-law model for most GRBs. The distribution of $E_p$ for episodes I and II range from tens of keV to 1000 keV with a lognormal fit and there are no significant differences between them. Moreover, we do not find significant relationships between $\varepsilon$ and $E_p$ for the two emission episodes. Those results suggest that these two-episode emission components likely share the same physical origin.

We report results from very long baseline interferometric (VLBI) observations of the supermassive black hole in the Galactic center, Sgr A*, at 1.3 mm (230 GHz). The observations were performed in 2013 March using six VLBI stations in Hawaii, California, Arizona, and Chile. Compared to earlier observations, the addition of the APEX telescope in Chile almost doubles the longest baseline length in the array, provides additional {\it uv} coverage in the N-S direction, and leads to a spatial resolution of $\sim$30 $\mu$as ($\sim$3 Schwarzschild radii) for Sgr A*. The source is detected even at the longest baselines with visibility amplitudes of $\sim$4-13% of the total flux density. We argue that such flux densities cannot result from interstellar refractive scattering alone, but indicate the presence of compact intrinsic source structure on scales of $\sim$3 Schwarzschild radii. The measured nonzero closure phases rule out point-symmetric emission. We discuss our results in the context of simple geometric models that capture the basic characteristics and brightness distributions of disk- and jet-dominated models and show that both can reproduce the observed data. Common to these models are the brightness asymmetry, the orientation, and characteristic sizes, which are comparable to the expected size of the black hole shadow. Future 1.3 mm VLBI observations with an expanded array and better sensitivity will allow a more detailed imaging of the horizon-scale structure and bear the potential for a deep insight into the physical processes at the black hole boundary.

Centaurus A (Cen A) is a bright radio source associated with the nearby galaxy NGC 5128 where high-resolution radio observations can probe the jet at scales of less than a light-day. The South Pole Telescope (SPT) and the Atacama Pathfinder Experiment (APEX) performed a single-baseline very-long-baseline interferometry (VLBI) observation of Cen A in January 2015 as part of VLBI receiver deployment for the SPT. We measure the correlated flux density of Cen A at a wavelength of 1.4 mm on a $\sim$7000 km (5 G$\lambda$) baseline. Ascribing this correlated flux density to the core, and with the use of a contemporaneous short-baseline flux density from a Submillimeter Array observation, we infer a core brightness temperature of $1.4 \times 10^{11}$ K. This is close to the equipartition brightness temperature, where the magnetic and relativistic particle energy densities are equal. Under the assumption of a circular Gaussian core component, we derive an upper limit to the core size $\phi = 34.0 \pm 1.8~\mu\textrm{as}$, corresponding to 120 Schwarzschild radii for a black hole mass of $5.5 \times 10^7 M_{\odot}$.

M87 is one of the nearest radio galaxies with a prominent jet extending from sub-pc to kpc-scales. Because of its proximity and large mass of the central black hole, it is one of the best radio sources to study jet formation. We aim at studying the physical conditions near the jet base at projected separations from the BH of $\sim7-100$ Schwarzschild radii ($R_{\rm sch}$). Global mm-VLBI Array (GMVA) observations at 86 GHz ($\lambda=3.5\,$mm) provide an angular resolution of $\sim50\mu$as, which corresponds to a spatial resolution of only $7~R_{\rm sch}$ and reach the small spatial scale. We use five GMVA data sets of M87 obtained during 2004--2015 and present new high angular resolution VLBI maps at 86GHz. In particular, we focus on the analysis of the brightness temperature, the jet ridge lines, and the jet to counter-jet ratio. The imaging reveals a parabolically expanding limb-brightened jet which emanates from a resolved VLBI core of $\sim(8-13) R_{\rm sch}$ size. The observed brightness temperature of the core at any epoch is $\sim(1-3)\times10^{10}\,$K, which is below the equipartition brightness temperature and suggests magnetic energy dominance at the jet base. We estimate the diameter of the jet at its base to be $\sim5 R_{\rm sch}$ assuming a self-similar jet structure. This suggests that the sheath of the jet may be anchored in the very inner portion of the accretion disk. The image stacking reveals faint emission at the center of the edge-brightened jet on sub-pc scales. We discuss its physical implication within the context of the spine-sheath structure of the jet.

Gamma-ray Burst (GRB) collimation has been inferred with the observations of achromatic steepening in GRB light curves, known as jet breaks. Identifying a jet break from a GRB afterglow lightcurve allows a measurement of the jet opening angle and true energetics of GRBs. In this paper, we reinvestigate this problem using a large sample of GRBs that have an optical jet break which is consistent with being achromatic in the X-ray band. Our sample includes 99 GRBs from February 1997 to March 2015 that have optical and, for Swift GRBs, X-ray lightcurves that are consistent with the jet break interpretation. Out of 99 GRBs we have studied, 55 GRBs are found to have temporal and spectral behaviors both before and after the break consistent with the theoretical predictions of the jet break models, respectively. These include 53 long/soft (Type II) and 2 short/hard (Type I) GRBs. Only 1 GRB is classified as the candidate of a jet break with energy injection. Another 41 and 3 GRBs are classified as the candidates with the lower and upper limits of the jet break time, respectively. The typical beaming correction factor $f_b^{-1} \sim 1000$ for Type II GRBs, suggesting an even higher total GRB event rate density in the universe. Both isotropic and jet-corrected energies have a wide span in their distributions. We also investigate several empirical correlations (Amati, Frail, Ghirlanda and Liang-Zhang) previously discussed in the literature. We find that in general most of these relations are less tight than before. The existence of early jet breaks and hence small opening angle jets, which were detected in the {\em Swfit era}, is most likely the source of scatter. If one limits the sample to jet breaks later than $10^4$ s, the Liang-Zhang relation remains tight and the Ghirlanda relation still exists. These relations are derived from Type II GRBs, and Type I GRBs usually deviate from them.

We analyze the long-term lightcurve of 3C 454.3 observed with Fermi/LAT and investigate its relation to the flux in the radio, optical, and X-ray bands. By fitting the 1-day binned GeV lightcurve with multiple Gaussian function (MGF), we propose that the typical variability timescale in the GeV band is 1--10 days. The GeV flux variation is accompanied by the spectral variation characterized as flux-tracking, i.e., "harder when brighter". The GeV flux is correlated with the optical and X-ray fluxes, and a weak correlation between gamma-ray flux and radio flux is also observed. The gamma-ray flux is not correlated with the optical linear polarization degree for the global lightcurves, but they show a correlation for the lightcurves before MJD 56000. The power density spectrum of the global lightcurve shows an obvious turnover at ~7.7 days, which may indicate a typical variability timescale of 3C 454.3 in the gamma-ray band. This is also consistent with the derived timescales by fitting the global lightcurve with MGF. The spectral evolution and an increase of the optical linear polarization degree along with the increase of the gamma-ray flux may indicate that the radiation particles are accelerated and the magnetic field is ordered by the shock processes during the outbursts. In addition, the nature of 3C 454.3 may be consistent with the self-organized criticality system, similar to Sagittarius A, and thus the outbursts are from plasmoid ejections driven by magnetic reconnection. This may further support the idea that the jet radiation regions are magnetized.

Motivating by the discovery of association between GW 170817 and sGRB 170817A, we present a comprehensive analysis for sGRBs observed with Fermi/GBM in 9 operation years and study the properties of sGRB 170817A -like events. We derive a catalog of 275 typical sGRBs and 48 sGRB 170817A-like weak events from the GBM data of 2217 GRBs. We visibly identify two patterns of their light curve, single episode (Pattern I, 61\% of the SGRBs) and multiple episodes (Pattern II, 39\% of the SGRBs). Their duration distribution shows a tentative bimodal feature. Their spectra can be fitted with a cutoff power-law model, except for 4 sGRBs, and the spectral indices normally distribute at $\Gamma=0.69\pm 0.40$. Their $E_p$ values show a tentative bimodal distribution with peaks at 145 keV and 500 keV. No correlation among $T_{90}$, $E_p$, and $\Gamma$ is found. GRB 170817A is a soft, weak sGRB with $ E_{p}=124\pm 106$ keV, $L_{\rm iso}=(5.67\pm4.65)\times10^{46}\rm ~erg~s^{-1}$, and $E_{\rm iso}=(3.23\pm2.65)\times10^{46}\rm ~erg$. It follows the $E_{\rm iso}-E_{\rm p}$ relation of typical short GRBs. Its lightcurve is of Pattern II. Two lightcurve patterns, together with the potential two components in the $E_{\rm p}$ and $T_{90}$ distributions, we suspect that the current sample may include two distinct types of sGRBs from different progenitors. sGRB 170817A-like events may be from NS-NS mergers and those sGRBs with a Pattern I lightcurve may be from another distinct type of compact binary.

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.

We present a comprehensive analysis of multiple wavelength observational data of the first GeV-selected narrow-line Seyfert 1 galaxy PMN J0948+0022. We derive its lightcurves in the \gamma-ray and X-ray bands from the data observed with Fermi/LAT and Swift/XRT, and make the optical and radio lightcurves by collecting the data from the literature. These lightcurves show significant flux variations. With the LAT data we show that this source is analogue to typical flat spectrum radio quasars in the $L_\gamma-\Gamma_\gamma$ plane, where $L_\gamma$ and $\Gamma_\gamma$ are the luminosity and spectral index in the LAT energy band. The gamma-ray flux is correlated with the V-band flux with a lag of ~44 days, and a moderate quasi-periodic oscillation (QPO) with a periodicity of ~490 days observed in the LAT lightcurve. A similar QPO signature is also found in the V-band lightcurve. The gamma-ray flux is not correlated with the radio flux in 15 GHz, and no similar QPO signature is found in a confidence level of 95%. Possible mechanisms of the QPO are discussed. We propose that gravitational wave observations in the future may clarify the current plausible models for the QPO.

The jet compositions, central engines, and progenitors of gamma-ray bursts (GRBs) remain open questions in GRB physics. Applying broadband observations, including GRB prompt emission and afterglow properties derived from {\em Fermi} and {\em Swift} data, as well as from Keck host-galaxy observations, we address these questions for the peculiar, bright GRB 110731A. By using the pair-opacity method, we derive $\Gamma_{0}>190$ during the prompt emission phase. Alternatively, we derive $\Gamma_{0} \approx 580$ and $\Gamma_{0} \approx 154$ by invoking the early-afterglow phase within the homogeneous density and wind cases, respectively. On the other hand, nondetection of a thermal component in the spectra suggests that the prompt emission is likely powered by dissipation of a Poynting-flux-dominated jet leading to synchrotron radiation in an optically thin region. The nondetection of a jet break in the X-ray and optical bands allows us to place a lower limit on the jet opening angle $\theta_{\rm j}>5.5^{\circ}$. Within a millisecond magnetar central engine scenario, we derive the period $P_{0}$ and polar magnetic field strength $B_{\rm p}$, which have extreme (but still allowed) values. The moderately short observed duration (7.3\,s) and relatively large redshift ($z=2.83$) places the burst as a "rest-frame short" GRB, so the progenitor of the burst is subject to debate. Its relatively large $f_{{\rm eff}, z}$ parameter (ratio of the 1\,s peak flux of a pseudo-GRB and the background flux) and a large physical offset from a potential host galaxy suggest that the progenitor of GRB 110731A may be a compact-star merger.

RadioAstron is a 10 m orbiting radio telescope mounted on the Spektr-R satellite, launched in 2011, performing Space Very Long Baseline Interferometry (SVLBI) observations supported by a global ground array of radio telescopes. With an apogee of about 350 000 km, it is offering for the first time the possibility to perform {\mu}as-resolution imaging in the cm-band. We present observations at 22 GHz of 3C 273, performed in 2014, designed to reach a maximum baseline of approximately nine Earth diameters. Reaching an angular resolution of 0.3 mas, we study a particularly low-activity state of the source, and estimate the nuclear region brightness temperature, comparing with the extreme one detected one year before during the RadioAstron early science period. We also make use of the VLBA-BU-BLAZAR survey data, at 43 GHz, to study the kinematics of the jet in a 1.5-year time window. We find that the nuclear brightness temperature is two orders of magnitude lower than the exceptionally high value detected in 2013 with RadioAstron at the same frequency (1.4x10^13 K, source-frame), and even one order of magnitude lower than the equipartition value. The kinematics analysis at 43 GHz shows that a new component was ejected 2 months after the 2013 epoch, visible also in our 22 GHz map presented here. Consequently this was located upstream of the core during the brightness temperature peak. These observations confirm that the previously detected extreme brightness temperature in 3C 273, exceeding the inverse Compton limit, is a short-lived phenomenon caused by a temporary departure from equipartition. Thus, the availability of interferometric baselines capable of providing {\mu}as angular resolution does not systematically imply measured brightness temperatures over the known physical limits for astrophysical sources.

This White Paper details the intentions and plans of the East Asian Very Long Baseline Interferometry (VLBI) community for pushing the frontiers of millimeter/submillimeter VLBI. To this end, we shall endeavor to actively promote coordinated efforts in the East Asia region. Our goal is to establish firm collaborations among the East Asia VLBI community in partnership with related institutes in North America and Europe and to expand existing global mm/submm VLBI arrays for (a) exploring the vicinity of black holes with an ultimate angular resolution down to 1~R$_{S}$ (Schwarzschild radius) and (b) investigating the dynamics of circumstellar gas in star-forming regions and late-type stars, and circumnuclear gas around active galactic nuclei (AGNs). In the first half of this White Paper, we highlight scientific accomplishments of the East Asia (EA) VLBI community. Various VLBI research results on M87, Sgr A*, blazars, narrow-line Seyfert~1 galaxies, and compact symmetric objects are described, and future visions of our VLBI science are briefly presented. Maser science of star formation, stellar evolution, and physics of accretion disks around AGNs are also discussed. A new vision for conducting multi-transition maser studies using mm/submm VLBI together with the Atacama Large Millimeter/Submillimeter Array (ALMA) is described. In the second half of this White Paper, we describe the EA community's vision for using mm/submm VLBI arrays in the framework or extended version of the Event Horizon Telescope (EHT) and the Global Millimeter VLBI Array (GMVA). The accomplishment of the aforementioned goal will maximize the overall scientific outcomes of mm/submm VLBI in the world.

We derive a simple analytical formula to describe the evolution of spectral index $\beta$ in the steep decay phase shaped by the curvature effect with assumption that the spectral parameters and Lorentz factor of jet shell is the same for different latitude. Here, the value of $\beta$ is estimated in 0.3$-$10keV energy band. For a spherical thin shell with a cutoff power law (CPL) intrinsic radiation spectrum, the spectral evolution can be read as a linear function of observer time. For the situation with Band function intrinsic radiation spectrum, the spectral evolution may be complex. If the observed break energy of radiation spectrum is larger than 10keV, the spectral evolution is the same as that shaped by jet shells with a CPL spectrum. If the observed break energy is less than 0.3keV, the value of $\beta$ would be a constant. Others, the spectral evolution can be approximated as a logarithmal function of the observer time in generally.

The curvature effect may be responsible for the steep decay phase observed in gamma-ray bursts. For testing the curvature effect with observations, the zero time point $t_0$ adopted to plot observer time and flux on a logarithmic scale should be appropriately selected. In practice, however, the true $t_0$ cannot be directly constrained from the data. Then, we move $t_0$ to a certain time in the steep decay phase, which can be easily identified. In this situation, we derive an analytical formula to describe the flux evolution of the steep decay phase. The analytical formula is read as $F_\nu\propto (1+\tilde t_{\rm obs}/{\tilde t_c})^{-\alpha}$ with $\alpha(\tilde{t}_{\rm obs})=2+{\int_{0}^{\log (1+\tilde{t}_{\rm obs}/{\tilde{t}_c})} {\beta(\tau)d[\log(1+\tau/{\tilde{t}_c})]}}/{\log (1 + {\tilde t}_{\rm obs}/{{\tilde t}_c})}$, where $F_\nu$ is the flux observed at frequency $\nu$, $\tilde t_{\rm obs}$ is the observer time by setting zero time point $t_0$ at a certain time in the steep decay phase, $\beta$ is the spectral index estimated around $\nu$, and ${\tilde t}_c$ is the decay timescale of the phase with $\tilde{t}_{\rm obs}{\geqslant}0$. We test the analytical formula with the data from numerical calculations. It is found that the analytical formula presents a well estimation about the evolution of flux shaped by the curvature effect. Our analytical formula can be used to confront the curvature effect with observations and estimate the decay timescale of the steep decay phase.

We aim to test if a blazar candidate of uncertain-type (BCU) in the third Fermi active galactic nuclei catalog (3LAC) can be potentially classified as a BL Lac object or a flat spectrum radio quasar (FSRQ) by performing a statistical analysis of its broadband spectral properties. We find that 34% of the radio-selected BCUs (583 BCUs) are BL Lac-like and 20% of them are FSRQ-like, which maybe within 90% level of confidence. Similarly, 77.3% of the X-ray selected BCUs (176 BCUs) are evaluated as BL Lac-like and 6.8% of them may be FSRQ-like sources. And 88.7% of the BL Lac-like BCUs that have synchrotron peak frequencies available are high synchrotron peaked BL Lacs in the X-ray selected BCUs. The percentages are accordingly 62% and 7.3% in the sample of 124 optical-selected BCUs. The high ratio of source numbers of the BL Lac-like to the FSRQ-like BCUs in the X-ray and optically selected BCU samples is due to the selection effect. Examining the consistency between our evaluation and spectroscopic identification case by case with a sample of 78 radio-selected BCUs, it is found that the statistical analysis and its resulting classifications agree with the results of the optical follow-up spectroscopic observations. Our observation campaign for high-|\rho_{s}| BCUs , i.e., |\rho_{s}|>0.8, selected with our method is ongoing.

Einstein's General Theory of Relativity (GR) successfully describes gravity. The most fundamental predictions of GR are black holes (BHs), but in spite of many convincing BH candidates in the Universe, there is no conclusive experimental proof of their existence using astronomical observations in the electromagnetic spectrum. Are BHs real astrophysical objects? Does GR hold in its most extreme limit or are alternatives needed? The prime target to address these fundamental questions is in the center of our own Galaxy, which hosts the closest and best-constrained supermassive BH candidate in the Universe, Sagittarius A* (Sgr A*). Three different types of experiments hold the promise to test GR in a strong-field regime using observations of Sgr A* with new-generation instruments. The first experiment aims to image the relativistic plasma emission which surrounds the event horizon and forms a "shadow" cast against the background, whose predicted size (~50 microarcseconds) can now be resolved by upcoming VLBI experiments at mm-waves such as the Event Horizon Telescope (EHT). The second experiment aims to monitor stars orbiting Sgr A* with the upcoming near-infrared interferometer GRAVITY at the Very Large Telescope (VLT). The third experiment aims to time a radio pulsar in tight orbit about Sgr A* using radio telescopes (including the Atacama Large Millimeter Array or ALMA). The BlackHoleCam project exploits the synergy between these three different techniques and aims to measure the main BH parameters with sufficient precision to provide fundamental tests of GR and probe the spacetime around a BH in any metric theory of gravity. Here, we review our current knowledge of the physical properties of Sgr A* as well as the current status of such experimental efforts towards imaging the event horizon, measuring stellar orbits, and timing pulsars around Sgr A*.

M87 is one of the nearest radio galaxies with a central SMBH and a prominent relativistic jet. Due to its close distance to the observer and the large SMBH mass, the source is one of the best laboratories to obtain strong observational constraints on the theoretical models for the formation and evolution of the AGN jets. In this article, we present preliminary results from our ongoing observational study about the innermost jet of M87 at an ultra-high resolution of $\sim$50${\mu}$as achieved by the Global Millimeter-VLBI Array (GMVA). The data obtained between 2004 and 2015 clearly show limb-brightened jets at extreme resolution and sensitivity. Our preliminary analysis reveals that the innermost jet expands in an edge-brightened cone structure (parabolic shape) but with the jet expansion profile slightly different from the outer regions of the jet. Brightness temperatures of the VLBI core obtained from cm- to mm-wavelengths show a systematic evolution, which can be interpreted as the evolution as a function of distance from the BH. We also adopt an alternative imaging algorithm, BSMEM, to test reliable imaging at higher angular resolution than provided by the standard CLEAN method (i.e. super-resolution). A demonstration with a VLBA 7mm example data set shows consistent results with a near-in-time 3mm VLBI image. Application of the method to the 2009 GMVA data yields an image with remarkable fine-scale structures that have been never imaged before. We present a brief interpretation of the complexity in the structure.

General Relativity predicts that the emission close to a black hole must be lensed by its strong gravitational field, illuminating the last photon orbit. This results in a dark circular area known as the black hole 'shadow'. The Event Horizon Telescope (EHT) is a (sub)mm VLBI network capable of Schwarzschild-radius resolution on Sagittarius A* (or Sgr A*), the 4 million solar mass black hole at the Galactic Center. The goals of the Sgr A* observations include resolving and measuring the details of its morphology. However, EHT data are sparse in the visibility domain, complicating reliable detailed image reconstruction. Therefore, direct pixel imaging should be complemented by other approaches. Using simulated EHT data from a black hole emission model we consider an approach to Sgr A* image reconstruction based on a simple and computationally efficient analytical model that produces images similar to the synthetic ones. The model consists of an eccentric ring with a brightness gradient and a two-dimensional Gaussian. These elemental forms have closed functional representations in the visibility domain, which lowers the computational overhead of fitting the model to the EHT observations. For model fitting we use a version of the Markov chain Monte-Carlo (MCMC) algorithm based on the Metropolis-Hastings sampler with replica exchange. Over a series of simulations we demonstrate that our model can be used for determining geometric measures of a black hole, thus providing information on the shadow size, linking General Relativity with accretion theory.

M87 is arguably the best supermassive black hole (BH) to explore the jet and/or accretion physics due to its proximity and fruitful high-resolution multi-waveband observations. We model the multi-wavelength spectral energy distribution (SED) of M87 core that observed at a scale of 0.4 arcsec ($\sim 10^5R_{\rm g}$, $R_{\rm g}$ is gravitational radius) as recently presented by Prieto et al. Similar to Sgr A*, we find that the millimeter bump as observed by Atacama Large Millimeter/submillimeter Array (ALMA) can be modeled by the synchrotron emission of the thermal electrons in advection dominated accretion flow (ADAF), while the low-frequency radio emission and X-ray emission may dominantly come from the jet. The millimeter radiation from ADAF dominantly come from the region within $10R_{\rm g}$, which is roughly consistent with the recent very long baseline interferometry observations at 230\,GHz. We further calculate the Faraday rotation measure (RM) from both ADAF and jet models, and find that the RM predicted from the ADAF is roughly consistent with the measured value while the RM predicted from the jet is much higher if jet velocity close to the BH is low or moderate (e.g., $v_{\rm jet}\lesssim0.6\,c$). With the constraints from the SED modeling and RM, we find that the accretion rate close to the BH horizon is $\sim (0.2-1)\times10^{-3}{M}_{\odot} \rm yr^{-1}\ll\dot{\it M}_{\rm B}\sim 0.2\it {M}_{\odot} \rm yr^{-1}$ ($\dot{M}_{\rm B}$ is Bondi accretion rate), where the electron density profile, $n_{\rm e}\propto r^{\sim -1}$, in the accretion flow is consistent with that determined from X-ray observation inside the Bondi radius and recent numerical simulations.

The radio emission from Sgr A$^\ast$ is thought to be powered by accretion onto a supermassive black hole of $\sim\! 4\times10^6~ \rm{M}_\odot$ at the Galactic Center. At millimeter wavelengths, Very Long Baseline Interferometry (VLBI) observations can directly resolve the bright innermost accretion region of Sgr A$^\ast$. Motivated by the addition of many sensitive, long baselines in the north-south direction, we developed a full VLBI capability at the Large Millimeter Telescope Alfonso Serrano (LMT). We successfully detected Sgr A$^\ast$ at 3.5~mm with an array consisting of 6 Very Long Baseline Array telescopes and the LMT. We model the source as an elliptical Gaussian brightness distribution and estimate the scattered size and orientation of the source from closure amplitude and self-calibration analysis, obtaining consistent results between methods and epochs. We then use the known scattering kernel to determine the intrinsic two dimensional source size at 3.5 mm: $(147\pm7~\mu\rm{as}) \times (120\pm12~\mu\rm{as})$, at position angle $88^\circ\pm7^\circ$ east of north. Finally, we detect non-zero closure phases on some baseline triangles, but we show that these are consistent with being introduced by refractive scattering in the interstellar medium and do not require intrinsic source asymmetry to explain.

The Galactic Center black hole Sagittarius A* (Sgr A*) is a prime observing target for the Event Horizon Telescope (EHT), which can resolve the 1.3 mm emission from this source on angular scales comparable to that of the general relativistic shadow. Previous EHT observations have used visibility amplitudes to infer the morphology of the millimeter-wavelength emission. Potentially much richer source information is contained in the phases. We report on 1.3 mm phase information on Sgr A* obtained with the EHT on a total of 13 observing nights over 4 years. Closure phases, the sum of visibility phases along a closed triangle of interferometer baselines, are used because they are robust against phase corruptions introduced by instrumentation and the rapidly variable atmosphere. The median closure phase on a triangle including telescopes in California, Hawaii, and Arizona is nonzero. This result conclusively demonstrates that the millimeter emission is asymmetric on scales of a few Schwarzschild radii and can be used to break 180-degree rotational ambiguities inherent from amplitude data alone. The stability of the sign of the closure phase over most observing nights indicates persistent asymmetry in the image of Sgr A* that is not obscured by refraction due to interstellar electrons along the line of sight.

We study the contribution of Galactic sources to the flux of astrophysical neutrinos recently observed by the IceCube Collaboration. We show that in the simplest model of homogeneous and isotropic cosmic ray diffusion in the Milky Way the Galactic diffuse neutrino emission consistent with $\gamma$-ray (Fermi-LAT) and cosmic ray data (KASCADE, KASCADE-Grande and CREAM) is expected to account for only $4\%-8\%$ of the IceCube flux above 60 TeV. Direct neutrino emission from cosmic ray-gas ($pp$) interactions in the sources would require an unusually large average opacity above 0.01. On the other hand, we find that the IceCube events already probe Galactic neutrino scenarios via the distribution of event arrival directions. Based on the latter, we show that most Galactic scenarios can only have a limited contribution to the astrophysical signal: diffuse Galactic emission ($\lesssim50\%$), quasi-diffuse emission of neutrino sources ($\lesssim65\%$), extended diffuse emission from the Fermi Bubbles ($\lesssim25\%$) or unidentified TeV $\gamma$-ray sources ($\lesssim25\%$). The arguments discussed here leave, at present, dark matter decay unconstrained.

The cosmic evolution of gamma-ray burst (GRB) luminosity is essential for revealing the GRB physics and for using GRBs as cosmological probes. We investigate the luminosity evolution of long GRBs with a large sample of 258 {\em Swift}/BAT GRBs. Parameterized the peak luminosity of individual GRBs evolves as $L_{\rm p}\propto{\rm }(1+z)^{k}$, we get $k=1.49\pm0.19$ using the non-parametric $\tau$ statistics method without considering observational biases of GRB trigger and redshift measurement. By modeling these biases with the observed peak flux and characterizing the peak luminosity function of long GRBs as a smoothly broken power-law with a break that evolves as $L_{\rm b}\propto (1+z)^{k_{\rm b}}$, we obtain $k_{\rm b}=1.14^{+0.99}_{-0.47}$ through simulations based on assumption that the long GRB rate follows the star formation rate (SFR) incorporating with cosmic metallicity history. The derived $k$ and $k_b$ values are systematically smaller than that reported in previous papers. By removing the observational biases of the GRB trigger and redshift measurement based on our simulation analysis, we generate mock {\em complete} samples of 258 and 1000 GRBs to examine how these biases affects on the $\tau$ statistics method. We get $k=0.94\pm 0.14$ and $k=0.80\pm 0.09$ for the two samples, indicating that these observational biases may lead to overestimate the $k$ value. With the large uncertain of $k_b$ derived from our simulation analysis, one even cannot convincingly argue a robust evolution feature of the GRB luminosity.

The black hole in the center of the Galaxy, associated with the compact source Sagittarius A* (Sgr A*), is predicted to cast a shadow upon the emission of the surrounding plasma flow, which encodes the influence of general relativity in the strong-field regime. The Event Horizon Telescope (EHT) is a Very Long Baseline Interferometry (VLBI) network with a goal of imaging nearby supermassive black holes (in particular Sgr A* and M87) with angular resolution sufficient to observe strong gravity effects near the event horizon. General relativistic magnetohydrodynamic (GRMHD) simulations show that radio emission from Sgr A* exhibits vari- ability on timescales of minutes, much shorter than the duration of a typical VLBI imaging experiment, which usually takes several hours. A changing source structure during the observations, however, violates one of the basic assumptions needed for aperture synthesis in radio interferometry imaging to work. By simulating realistic EHT observations of a model movie of Sgr A*, we demonstrate that an image of the average quiescent emission, featuring the characteristic black hole shadow and photon ring predicted by general relativity, can nonetheless be obtained by observing over multiple days and subsequent processing of the visibilities (scaling, averaging, and smoothing) before imaging. Moreover, it is shown that this procedure can be combined with an existing method to mitigate the effects of interstellar scattering. Taken together, these techniques allow the black hole shadow in the Galactic center to be recovered on the reconstructed image.

Near a black hole, differential rotation of a magnetized accretion disk is thought to produce an instability that amplifies weak magnetic fields, driving accretion and outflow. These magnetic fields would naturally give rise to the observed synchrotron emission in galaxy cores and to the formation of relativistic jets, but no observations to date have been able to resolve the expected horizon-scale magnetic-field structure. We report interferometric observations at 1.3-millimeter wavelength that spatially resolve the linearly polarized emission from the Galactic Center supermassive black hole, Sagittarius A*. We have found evidence for partially ordered fields near the event horizon, on scales of ~6 Schwarzschild radii, and we have detected and localized the intra-hour variability associated with these fields.

We select a sample of 34 gamma-ray bursts (GRBs) whose $\Gamma_0$ values are derived with the onset peaks observed in the afterglow lightcurves (except for GRB 060218 whose $\Gamma_0$ is estimated with its radio data), and investigate the correlations among $\Gamma_0$, the isotropic peak luminosity ($L_{\rm iso}$), and the peak energy ($E_{\rm p,z}$) of the $\nu f_\nu$ spectrum in the cosmological rest frame. An analysis of pair correlations among these observables well confirms the results reported by the previous papers. More interestingly, a tight correlation among $L_{\rm iso}$, $E_{\rm p,z}$, and $\Gamma_0$ is found from a multiple regression analysis, which takes the form of $L_{\rm iso} \propto E_{\rm p,z}^{1.34\pm 0.14} \Gamma_0^{1.32\pm 0.19}$ or $E_{\rm p,z} \propto L_{\rm iso}^{0.55\pm 0.06}\Gamma_0^{-0.50\pm 0.17}$. Nine other GRBs whose $\Gamma_0$ are derived via the pair production opacity constraint also follow such a correlation. Excluding GRB 060218, the $L_{\rm iso}-E_{\rm p,z}-\Gamma_0$ correlation is valid, and it even holds in the jet co-moving frame. However, GRB 060218 deviates the $L^{'}_{\rm iso}-E^{'}_{\rm p}$ relation of typical GRBs in the jet co-moving frame with $3\sigma$. We argue that the $L_{\rm iso} - E_{\rm p, z} - \Gamma_0$ correlation may be more physical than the $L_{\rm iso} - E_{\rm p,z}$ correlation, since physically the relationship between the observed $L_{\rm iso}$ and $E_{\rm p,z}$ not only depends on radiation physics, but also depends on the bulk motion of the jet. We explore the possible origins of this correlation and discuss its physical implications for understanding GRB jet composition and radiation mechanism.

We report on 230 GHz (1.3 mm) VLBI observations of M87 with the Event Horizon Telescope using antennas on Mauna Kea in Hawaii, Mt. Graham in Arizona and Cedar Flat in California. For the first time, we have acquired 230 GHz VLBI interferometric phase information on M87 through measurement of closure phase on the triangle of long baselines. Most of the measured closure phases are consistent with 0$^{\circ}$ as expected by physically-motivated models for 230 GHz structure such as jet models and accretion disk models. The brightness temperature of the event-horizon-scale structure is $\sim 1 \times 10^{10}$ K derived from the compact flux density of $\sim 1$ Jy and the angular size of $\sim 40 $ $\rm \mu$as $\sim$ 5.5 $R_{{\rm s}}$, which is broadly consistent with the peak brightness of the radio cores at 1-86 GHz located within $\sim 10^2$ $R_{{\rm s}}$. Our observations occurred in the middle of an enhancement in very-high-energy (VHE) $\rm \gamma$-ray flux, presumably originating in the vicinity of the central black hole. Our measurements, combined with results of multi-wavelength observations, favor a scenario in which the VHE region has an extended size of $\sim$20-60 $R_{{\rm s}}$.

We report about a 230 GHz very long baseline interferometry (VLBI) fringe finder observation of blazar 3C 279 with the APEX telescope in Chile, the phased submillimeter array (SMA), and the SMT of the Arizona Radio Observatory (ARO). We installed VLBI equipment and measured the APEX station position to 1 cm accuracy (1 sigma). We then observed 3C 279 on 2012 May 7 in a 5 hour 230 GHz VLBI track with baseline lengths of 2800 M$\lambda$ to 7200 M$\lambda$ and a finest fringe spacing of 28.6 micro-arcseconds. Fringes were detected on all baselines with SNRs of 12 to 55 in 420 s. The correlated flux density on the longest baseline was ~0.3 Jy/beam, out of a total flux density of 19.8 Jy. Visibility data suggest an emission region <38 uas in size, and at least two components, possibly polarized. We find a lower limit of the brightness temperature of the inner jet region of about 10^10 K. Lastly, we find an upper limit of 20% on the linear polarization fraction at a fringe spacing of ~38 uas. With APEX the angular resolution of 230 GHz VLBI improves to 28.6 uas. This allows one to resolve the last-photon ring around the Galactic Center black hole event horizon, expected to be 40 uas in diameter, and probe radio jet launching at unprecedented resolution, down to a few gravitational radii in galaxies like M 87. To probe the structure in the inner parsecs of 3C 279 in detail, follow-up observations with APEX and five other mm-VLBI stations have been conducted (March 2013) and are being analyzed.

X-ray flares of gamma-ray bursts (GRBs) are usually observed in the soft X-ray range and the spectral coverage is limited. In this paper, we present an analysis of 32 GRB X-ray flares that are simultaneously observed by both BAT and XRT on board the Swift mission, so a joint spectral analysis with a wider spectral coverage is possible. Our results show that the joint spectra of 19 flares are fitted with the absorbed single power-law or the Band function models. More interestingly, the joint spectra of the other 13 X-ray flares are fitted with the absorbed single power-law model plus a black body (BB) component. Phenomenally, the observed spectra of these 13 flares are analogous to several GRBs with a thermal component, but only with a much lower temperature of $kT=1\sim 3$ keV. Assuming that the thermal emission is the photosphere emission of the GRB fireball, we derive the fireball properties of the 13 flares that have redshift measurements, such as the bulk Lorentz factor $\Gamma_{\rm ph}$ of the outflow. The derived $\Gamma_{\rm ph}$ range from $50$ to $150$ and a relation of $\Gamma_{\rm ph}$ to the thermal emission luminosity is found. It is consistent with the $\Gamma_0-L_{\rm iso}$ relation that are derived for the prompt gamma-ray emission. We discuss the physical implications of these results within the content of jet composition and radiation mechanism of GRBs and X-ray flares.

The image of the emission surrounding the black hole in the center of the Milky Way is predicted to exhibit the imprint of general relativistic (GR) effects, including the existence of a shadow feature and a photon ring of diameter ~50 microarcseconds. Structure on these scales can be resolved by millimeter-wavelength very long baseline interferometry (VLBI). However, strong-field GR features of interest will be blurred at lambda >= 1.3 mm due to scattering by interstellar electrons. The scattering properties are well understood over most of the relevant range of baseline lengths, suggesting that the scattering may be (mostly) invertible. We simulate observations of a model image of Sgr A* and demonstrate that the effects of scattering can indeed be mitigated by correcting the visibilities before reconstructing the image. This technique is also applicable to Sgr A* at longer wavelengths.

We present a spectral analysis for a sample of redshift known GRBs observed with {\em Fermi}/GBM. Together with the results derived from our systematical spectral energy distribution modeling with the leptonic models for a {\em Fermi}/LAT blazar sample, we compare the distributions of the GRBs and the blazars by plotting the synchrotron peak luminosity ($L_{\rm s}$) and the corresponding peak photon energy $E_{\rm s}$ of blazars in the $L_{\rm p}-E_{\rm p}$ plane of GRBs, where $L_{\rm p}$ and $E_{\rm p}$ are the peak luminosity and peak photon energy of the GRB time-integrated $\nu f_\nu$ spectrum, respectively. The GRBs are in the high-$L_{\rm p}$, high-$E_{\rm p}$ corner of the plane and a tight $L_{\rm p}-E_{\rm p}$ relation is found, i.e., $L_{\rm p}\propto E_{\rm p}^{2.13^{+0.54}_{-0.46}}$. Both FSRQs and LBLs are clustered in the low-$E_{\rm p}$, low-$L_{\rm p}$ corner. IBLs and HBLs have $E_{\rm s}\sim 2\times 10^{-3} - 10^{2}$ keV and $L_{\rm s} \sim 10^{44} - 10^{47}$ erg s$^{-1}$, but no dependence of $L_{\rm s}$ on $E_{\rm s}$ is found. We show that the tight $L_p-E_p$ relation of GRBs is potentially explained with the synchrotron radiation of fast-cooling electrons in a highly magnetized ejecta, and the weak anti-correlation of $L_{\rm s}-E_{\rm s}$ for FSRQs and LBLs may be attributed to synchrotron radiation of slow-cooling electrons in a moderately magnetized ejecta. The distributions of IBLs and HBLs in the $L_{\rm p}-E_{\rm p}$ plane may be interpreted with synchrotron radiation of fast-cooling electrons in a matter-dominated ejecta. These results may present a unified picture for the radiation physics of relativistic jets in GRBs and blazars within the framework of the leptonic synchrotron radiation models.

We investigate whether a subset of high-energy events observed by IceCube may be due to neutrinos from Sagittarius A*. We check both spatial and temporal coincidences of IceCube events with other transient activities of Sagittarius A*. Among the seven IceCube shower events nearest to the galactic center, we have found that event 25 has a time very close to (around three hours after) the brightest X-ray flare of Sagittarius A* observed by the Chandra X-ray Observatory with a p-value of 0.9%. Furthermore, two of the seven events occurred within one day of each other (there is a 1.6% probability that this would occur for a random distribution in time). Thus, the determination that some IceCube events occur at similar times as X-ray flares and others occur in a burst could be the smoking gun that Sagittarius A* is a point source of very high energy neutrinos. We point out that if IceCube galactic center neutrino events originate from charged pion decays, then TeV gamma rays should come from neutral pion decays at a similar rate. We show that the CTA, HAWC, H.E.S.S. and VERITAS experiments should be sensitive enough to test this.

We jointly analyze the gamma-ray burst (GRB) data observed with BAT and XRT on board the Swift mission to present a global view on the internal energy dissipation processes in GRBs, including precursors, prompt gamma-ray emission, extended soft gamma-ray emission, and late X-ray flares. The Bayesian block method is utilized to analyze the BAT lightcurves to identify various emission episodes. Our results suggest that these emission components likely share a same physical origin, which is repeated activation of the GRB central engine. What we observe in the gamma-ray band may be the tip-of-iceberg of more extended underlying activities. The precursor emission, which is detected in about 10% of {\em Swift} GRBs, is preferably detected in those GRBs that have a massive star core-collapse origin. The soft extended emission (EE) tail, on the other hand, is preferably detected in those GRBs that have a compact star merger origin. Bright X-ray emission is detected during the BAT quiescent phases prior to subsequent gamma-ray peaks, implying that X-ray emission may be detectable prior the BAT trigger time. Future GRB alert instruments with soft X-ray capability would be essential to reveal the early stage of GRB central engine activities, sheding light into jet composition and jet launching mechanism in GRBs.

The Event Horizon Telescope (EHT) is a project to assemble a Very Long Baseline Interferometry (VLBI) network of mm wavelength dishes that can resolve strong field General Relativistic signatures near a supermassive black hole. As planned, the EHT will include enough dishes to enable imaging of the predicted black hole "shadow", a feature caused by severe light bending at the black hole boundary. The center of M87, a giant elliptical galaxy, presents one of the most interesting EHT targets as it exhibits a relativistic jet, offering the additional possibility of studying jet genesis on Schwarzschild radius scales. Fully relativistic models of the M87 jet that fit all existing observational constraints now allow horizon-scale images to be generated. We perform realistic VLBI simulations of M87 model images to examine detectability of the black shadow with the EHT, focusing on a sequence of model images with a changing jet mass load radius. When the jet is launched close to the black hole, the shadow is clearly visible both at 230 and 345 GHz. The EHT array with a resolution of 20-30$\mu$as resolution ($\sim$2-4 Schwarzschild radii) is able to image this feature independent of any theoretical models and we show that imaging methods used to process data from optical interferometers are applicable and effective for EHT data sets. We demonstrate that the EHT is also capable of tracing real-time structural changes on a few Schwarzschild radii scales, such as those implicated by VHE flaring activity of M87. While inclusion of ALMA in the EHT is critical for shadow imaging, generally the array is robust against loss of a station.

We fit the spectral energy distributions (SEDs) of a GeV-TeV FSRQ sample with the leptonic model. Their gamma_min of the relativistic electron distributions, which significantly affect the estimates of the jet properties, are constrained, with a typical value of 48. Their jet power, magnetized parameter, radiation efficiency, and jet production/radiation rates per central black hole (BH) mass are derived and compared to that of BL Lacs. We show that he FSRQ jets may be dominated by the Poynting flux and have a high radiation efficiency, whereas the BL Lac jets are likely dominated by particles and have a lower radiation efficiency than FSRQs. Being different from BL Lacs, the jet powers of FSRQs are proportional to their central BH masses. The jet production and radiation rates of the FSRQs distribute in narrow ranges and are correlated with each other, whereas no similar feature is found for the BL Lacs. We also show that the jet power is correlated with the cavity kinetic power, and the magnetic field energy in the jets may provide the cavity kinetic energy of FSRQs and the kinetic energy of cold protons in the jets may be crucial for cavity kinetic energy of BL Lacs. We suggest that the dominating formation mechanism of FSRQ jets may be the BZ process, but BL Lac jets may be produced via the BP and/or BZ processes, depending on the structures and accretion rates of accretion disks.

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.

We perform a geometric analysis for the sky map of the IceCube TeV-PeV neutrino excess and test its compatibility with the sky map of decaying dark matter signals in our galaxy. Using both Kolmogorov-Smirnov and the likelihood-ratio tests, we have found that the observed event sky map prefers to have a combination of the galactic dark matter and a homogeneous background contributions, compared to a purely galactic dark matter origin. For the assumption that the galactic dark matter is responsible for all neutrino excess, the current data can also exclude a wide range of dark matter profiles except flatter profiles such as the isothermal one. We also consider several representative decaying dark matter spectra, which can provide a good fit to the observed spectrum at IceCube with a dark matter lifetime of around 12 orders of magnitude longer than the age of the universe.

An international consortium is presently constructing a beamformer for the Atacama Large Millimeter/submillimeter Array (ALMA) in Chile that will be available as a facility instrument. The beamformer will aggregate the entire collecting area of the array into a single, very large aperture. The extraordinary sensitivity of phased ALMA, combined with the extremely fine angular resolution available on baselines to the Northern Hemisphere, will enable transformational new very long baseline interferometry (VLBI) observations in Bands 6 and 7 (1.3 and 0.8 mm) and provide substantial improvements to existing VLBI arrays in Bands 1 and 3 (7 and 3 mm). The ALMA beamformer will have impact on a variety of scientific topics, including accretion and outflow processes around black holes in active galactic nuclei (AGN), tests of general relativity near black holes, jet launch and collimation from AGN and microquasars, pulsar and magnetar emission processes, the chemical history of the universe and the evolution of fundamental constants across cosmic time, maser science, and astrometry.

We continue our systematic statistical study on optical afterglow data of gamma-ray bursts (GRBs). We present the apparent magnitude distributions of early optical afterglows at different epochs (t= 10^2 s, t = 10^3 s, and 1 hour) for the optical lightcurves of a sample of 93 GRBs (the global sample), and for sub-samples with an afterglow onset bump or a shallow decay segment. For the onset sample and shallow decay sample we also present the brightness distribution at the peak time t_{p} and break time t_{b}, respectively. All the distributions can be fit with Gaussian functions. We further perform Monte Carlo simulations to infer the luminosity function of GRB optical emission at the rest-frame time 10^3 seconds, t_{p}, and t_{b}, respectively. Our results show that a single power-law luminosity function is adequate to model the data, with indices -1.40+/-0.10, -1.06+/- 0.16, and -1.54\+/- 0.22, respectively. Based on the derived rest-frame 10^3 s luminosity function, we generate the intrinsic distribution of the R-band apparent magnitude M_{R} at the observed time 10^{3} seconds post trigger, which peaks at M_{R}=22.5 mag. The fraction of GRBs whose R-band magnitude is fainter than 22 mag, and 25 mag and at the observer time 10^3 seconds are ~63% and ~25%, respectively. The detection probabilities of the optical afterglows with ground-based robotic telescopes and UVOT onboard {Swift} are roughly consistent with that inferred from this intrinsic M_{R} distribution, indicating that the variations of the dark GRB fraction among the samples with different telescopes may be due to the observational selection effect, although the existence of an intrinsically dark GRB population cannot be ruled out.

We continue our systematic statistical study of various components in gamma-ray burst (GRB) optical lightcurves. We decompose the early onset bump and the late re-brightening bump with empirical fits and analyze their statistical properties. Among the 146 GRBs that have well-sampled optical lightcurves, the onset and re-brightening bumps are observed in 38 and 26 GRBs, respectively. It is found that the typical rising and decaying slopes for both the onset and re-brightening bumps are ~1.5 and -1.15, respectively. No early onset bumps in the X-ray band are detected to be associated with the optical onset bumps, while an X-ray re-brightening bump is detected for half of the re-brightening optical bumps. The peak luminosity is anti-correlated with the peak time, L_p\propto t_{p}^{-1.81+/-0.32} and L_p\propto t_{p}^{-0.83+/-0.17} for the onset and re-brightening bumps, respectively. Both L_p and the isotropic energy release of the onset bumps are correlated with E_{\gamma, iso}, whereas no similar correlation is found for the re-brightening bumps. Taking the onset bumps as probes for the properties of the fireballs and their ambient medium, we find that the typical power-law index of the relativistic electrons is 2.5 and the medium density profile behaves as n\propto r^{-1} within the framework of the synchrotron external shock models. With the medium density profile obtained from our analysis, we also confirm the correlation between initial Lorentz factor (\Gamma_0) and E_{\gamma, iso} in our previous work. The jet component that produces the re-brightening bump seems to be on-axis and independent of the prompt emission jet component. Its typical kinetic energy budget would be about one order of magnitude larger than the prompt emission component, but with a lower \Gamma_0, typically several tens.

We report results from 5-day VLBI observations of the well-known quasar 3C 279 at 1.3 mm (230 GHz) in 2011. The measured nonzero closure phases on triangles including stations in Arizona, California and Hawaii indicate that the source structure is spatially resolved. We find an unusual inner jet direction at scales of $\sim$1 parsec extending along the northwest-southeast direction (PA = $127^{\circ}\pm3^{\circ}$), as opposed to other (previously) reported measurements on scales of a few parsecs showing inner jet direction extending to the southwest. The 1.3 mm structure corresponds closely with that observed in the central region of quasi-simultaneous super-resolution VLBA images at 7 mm. The closure phase changed significantly on the last day when compared with the rest of observations, indicating that the inner jet structure may be variable on daily timescales. The observed new direction of the inner jet shows inconsistency with the prediction of a class of jet precession models. Our observations indicate a brightness temperature of $\sim 8\times10^{10}$ K in the 1.3 mm core, much lower than that at centimeter wavelengths. Observations with better uv coverage and sensitivity in the coming years will allow the discrimination between different structure models and will provide direct images of the inner regions of the jet with 20--30 $\mu$as (5--7 light months) resolution.

GRB 121027A is un-usual with its extremely long-lasting, energetic X-ray flares. The total energy release in X-ray flares is about one order of magnitude higher than prompt gamma-rays, making it special from most long GRBs. We show that while the prompt gamma-ray emission satisfies the empirical E_{iso}-E_{p} relation of typical long GRBs, the X-ray flares, whose spectra can be fit with a cutoff-power-law model with well-constrained E_p, significantly deviate from such a relation. Nonetheless, a time-resolved spectral analysis of X-ray flares suggest that the X-ray emission is consistent with the L_{iso}-E_{p} relation of long GRBs. We constrain the minimum Lorentz factor of the X-ray flares to be ~14, which is consistent with the Gamma-L_{iso} relation. Our results imply that prompt gamma-ray emission and late X-ray flares share the similar radiation mechanism, but originate from the outflows with different Lorentz factors. We search for similar GRBs from the Swift GRB archives, and find that the z=6.29 GRB 050904 is a carbon copy of GRB 121027A. The long-lasting, energetic X-ray flares in these GRBs demand significant accretion at late times, which point towards a large-radius progenitor star.

The durations (T90) of 315 GRBs detected with Fermi/GBM (8-1000 keV) by 2011 September are calculated using the Bayesian Block method. We compare the T90 distributions between this sample and those derived from previous/current GRB missions. We show that the T90 distribution of this GRB sample is bimodal, with a statistical significance level being comparable to those derived from the BeppoSAX/GRBM sample and the Swift/BAT sample, but lower than that derived from the CGRO/BATSE sample. The short-to-long GRB number ratio is also much lower than that derived from the BATSE sample, i.e., 1:6.5 vs 1:3. We measure T90 in several bands, i.e., 8-15, 15-25, 25-50, 50-100, 100-350, and 350-1000 keV, to investigate the energy-dependence effect of the bimodal T90 distribution. It is found that the bimodal feature is well observed in the 50-100 and 100-350 keV bands, but is only marginally acceptable in the 25-50 keV and 350-1000 keV bands. The hypothesis of the bimodality is confidently rejected in the 8-15 and 15-25 keV bands. The T90 distributions in these bands are roughly consistent with those observed by missions with similar energy bands. The parameter T90 as a function of energy follows \bar T90 \propto E^{-0.20\pm 0.02} for long GRBs. Considering the erratic X-ray and optical flares, the duration of a burst would be even much longer for most GRBs. Our results, together with the observed extended emission of some short GRBs, indicate that the central engine activity time scale would be much longer than T90} for both long and short GRBs and the observed bimodal T90 distribution may be due to an instrumental selection effect.

The prompt emission spectra of gamma-ray bursts (GRBs) usually have a dominant component that is well described by a phenomenological "Band" function. The physical origin of this spectral component is debated. Although the traditional interpretation is synchrotron radiation of non-thermal electrons accelerated in internal shocks or magnetic dissipation regions, a growing trend in the community is to interpret this component as modified thermal emission from a dissipative photosphere of a GRB fireball. We analyze the time dependent spectrum of GRB 110721A detected by {\em Fermi} GBM and LAT, and pay special attention to the rapid evolution of the peak energy $E_p$. We define a "death line" of thermally-dominated dissipative photospheric emission in the $E_p - L$ plane, and show that $E_p$ of GRB 110721A at the earliest epoch has a very high $E_p \sim 15$ MeV that is beyond the "death line". Together with the finding that an additional "shoulder" component exists in this burst that is consistent with a photospheric origin, we suggest that at least for some bursts, the "Band" component is not from a dissipative photosphere, but must invoke a non-thermal origin (e.g. synchrotron or inverse Compton) in the optically thin region of a GRB outflow. We also suggest that the rapid "hard-to-soft" spectral evolution is consistent with the quick discharge of magnetic energy in a magnetically-dominated outflow in the optically thin region.

We present the first 1.3 mm (230 GHz) very long baseline interferometry model image of an AGN jet using closure phase techniques with a four-element array. The model image of the quasar 1924-292 was obtained with four telescopes at three observatories: the James Clerk Maxwell Telescope (JCMT) on Mauna Kea in Hawaii, the Arizona Radio Observatory's Submillimeter Telescope (SMT) in Arizona, and two telescopes of the Combined Array for Research in Millimeterwave Astronomy (CARMA) in California in April 2009. With the greatly improved resolution compared with previous observations and robust closure phase measurement, the inner jet structure of 1924-292 was spatially resolved. The inner jet extends to the northwest along a position angle of $-53^\circ$ at a distance of 0.38\,mas from the tentatively identified core, in agreement with the inner jet structure inferred from lower frequencies, and making a position angle difference of $\sim 80^{\circ}$ with respect to the cm-jet. The size of the compact core is 0.15\,pc with a brightness temperature of $1.2\times10^{11}$\,K. Compared with those measured at lower frequencies, the low brightness temperature may argue in favor of the decelerating jet model or particle-cascade models. The successful measurement of closure phase paves the way for imaging and time resolving Sgr A* and nearby AGN with the Event Horizon Telescope.

Well-sampled optical lightcurves of 146 gamma-ray bursts (GRBs) are compiled from the literature. By empirical fitting we identify eight possible emission components and summarize the results in a "synthetic" lightcurve. Both optical flare and early shallow-decay components are likely related to long-term central engine activities. We focus on their statistical properties in this paper. Twenty-four optical flares are obtained from 19 GRBs. The isotropic R-band energy is smaller than 1% of $E_{\gamma, \rm iso}$. The relation between isotropic luminosities of the flares and gamma-rays follows $L^{\rm F}_{\rm R, iso}\propto L_{{\gamma}, \rm iso}^{1.11\pm 0.27}$. Later flares tend to be wider and dimmer, i.e., $w^{\rm F}\sim t^{\rm F}_{\rm p}/2$ and $L^{\rm F}_{\rm R, iso}\propto [t^{\rm F}_{\rm p}/(1+z)]^{-1.15\pm0.15}$. The detection probability of the optical flares is much smaller than that of X-ray flares. An optical shallow decay segment is observed in 39 GRBs. The relation between the break time and break luminosity is a power-law, with an index of $-0.78\pm 0.08$, similar to that derived from X-ray flares. The X-ray and optical breaks are usually chromatic, but a tentative correlation is found. We suggest that similar to the prompt optical emission that tracks $\gamma$-rays, the optical flares are also related to the erratic behavior of the central engine. The shallow decay component is likely related to a long-lasting spinning-down central engine or piling up of flare materials onto the blastwave. Mixing of different emission components may be the reason of the diverse chromatic afterglow behaviors.

We present a time-resolved spectral analysis of 51 long and 11 short bright GRBs observed with the {\em Femri}/GBM, paying special attention to $E_{\rm p}$ evolution within a same burst. Among 8 single-pulse long GRBs, 5 show hard-to-soft evolution, while 3 show intensity-tracking. The multi-pulse long GRBs have more complicated patterns. Among the GRBs whose time-resolved spectrum is available for the first pulse, almost half (15/32 GRBs) show clear hard-to-soft evolution, and the other half (17/32 GRBs) show clear intensity-tracking. Later pulses typically show the tracking behavior, although a hard-to-soft evolution pattern was identified in the 2nd pulse of 2 GRBs whose pulses are well separated. Statistically, the hard-to-soft evolution pulses tend to be more asymmetric than the intensity-tracking ones, with a steeper rising wing than the falling wing. Short GRBs have $E_{\rm p}$ tracking intensity exclusively with the 16ms time resolution analysis. We performed a simulation analysis, and suggest that at least for some bursts, the late intensity-tracking pulses could be a consequence of overlapping hard-to-soft pulses. However, the fact that the intensity-tracking pattern exists in the first pulse of multi-pulse long GRBs and some single-pulse GRBs suggest that intensity tracking is an independent component, which may operate in some late pulses as well. For the GRBs with measured redshifts, we present a time-resolved $E_{\rm p}-L_{\gamma, \rm iso}$ correlation analysis and show that the scatter of the correlation is comparable to that of the global Amati/Yonetoku relation. We discuss the predictions of various radiation models regarding $E_{\rm p}$ evolution, as well as the possibility of a precession jet in GRBs. It seems that the data pose great challenge to all these models, and hold the key to unveil the physics of GRB prompt emission.

String/M theory compactifications with low energy supersymmetry tend to predict that dark matter has two components: axions and WIMPs \cite{1004.5138,1204.2795}. In accord with this, we show that the tentative 130 GeV gamma-line signal reported in \cite{1204.2797} can be interpreted as arising from the annihilation of 145 GeV mass, Wino-like WIMPs into a Z-boson and a photon. In this context, the signal implies a second component of dark matter which we interpret as being composed of axions - the relative Wino/Axion abundances being approximately equal. Further predictions are implied: signals in both diffuse and monochromatic photons from dwarf spheroidal galaxies; monochromatic photons with energy 145 GeV; for the LHC, the Higgs boson mass has been predicted in this framework \cite{1112.1059}, and the current Higgs limits provide interesting constraints on the mass of the Gluino.

PKS 1749+096 is a BL Lac object showing weak extended jet emission to the northeast of the compact VLBI core on parsec scales. We aim at better understanding the jet kinematics and variability of this source and finding clues that may applicable to other BL Lac objects. The jet was studied with multi-epoch multi-frequency high-resolution VLBI observations. The jet is characterized by a one-sided curved morphology at all epochs and all frequencies. The VLBI core, located at the southern end of the jet, was identified based on its spectral properties. The equipartition magnetic field of the core was investigated, through which we derived a Doppler factor of 5, largely consistent with that derived from kinematics (component C5). The study of the detailed jet kinematics at 22 and 15 GHz, spanning a period of more than 10 years, indicates the possible existence of a bimodal distribution of the jet apparent speed. Ballistic and non-ballistic components are found to coexist in the jet. Superluminal motions in the range of 5-21 c were measured in 11 distinct components. We estimated the physical jet parameters with the minimum Lorentz factor of 10.2 and Doppler factors in the range of 10.2-20.4 (component C5). The coincidence in time of the component's ejection and flares supports the idea that, at least in PKS 1749+096, ejection of new jet components is connected with major outbursts in flux density. For the best-traced component (C5) we found that the flux density decays rapidly as it travels downstream the jet, accompanied by a steepening of its spectra, which argues in favor of a contribution of inverse Compton cooling. These properties make PKS 1749+096 a suitable target for an intensive monitoring to decipher the variability phenomenon of BL Lac objects.

Apparent redshift dependence of the jet opening angles ($\theta_{\rm j}$) of gamma-ray bursts (GRBs) is observed from current GRB sample. We investigate whether this dependence can be explained with instrumental selection effects and observational biases by a bootstrapping method. Assuming that (1) the GRB rate follows the star formation history and the cosmic metallicity history and (2) the intrinsic distributions of the jet-corrected luminosity ($L_{\rm \gamma}$) and $\theta_{\rm j}$ are a Gaussian or a power-law function, we generate a mock {\em Swift}/BAT sample by considering various instrumental selection effects, including the flux threshold and the trigger probability of BAT, the probabilities of a GRB jet pointing to the instrument solid angle and the probability of redshift measurement. Our results well reproduce the observed $\theta_{\rm j}-z$ dependence. We find that in case of $L_{\gamma}\propto \theta_{\rm j}^2$ good consistency between the mock and observed samples can be obtained, indicating that both $L_{\rm \gamma}$ and $\theta_{\rm j}$ are degenerate for a flux-limited sample. The parameter set $(L_{\rm \gamma}, \theta_{\rm j})=(4.9\times 10^{49} \rm {erg\ s}^{-1},\ 0.054 {rad})$ gives the best consistency for the current {\em Swift} GRB sample. Considering the beaming effect, the derived intrinsic local GRB rate accordingly is $2.85\times 10^2$ Gpc$^{-3}$ yr$^{-1}$, inferring that $\sim 0.59%$ of Type Ib/c SNe may be accompanied by a GRB.

NRAO 530 is an optically violent variable source and has been studied with multi-epoch multi-frequency high-resolution VLBI observations. NRAO 530 was monitored with the VLBA at three frequencies (22, 43 and 86 GHz) on 10 consecutive days in 2007 May during observations of the Galactic Center (Sgr A*). Furthermore, analysis of archival data of NRAO 530 at 15 GHz over the last ten years allows us to study its detailed jet kinematics. We identified the compact component located at the southern-end of the jet as the VLBI core, consistent with previous studies. The 10-d monitoring data at the 3 high frequencies were shown to produce high quality and self-consistent measurements of the component positions, from which we detected for the first time a two-dimensional frequency-dependent position shift. In addition, the repeated measurements also permit us to investigate the interday flux density and structure variability of NRAO 530. We find that it is more variable for the inner jet components than those further out. We obtained apparent velocities for eight jet components with $\beta_{\rm app} ranging from 2 to 26 c. Accordingly, we estimated physical jet parameters with the minimum Lorentz factor of 14 and Doppler factors in the range of 14--28 (component f). The changes in the morphology of NRAO 530 were related to the motion of separate jet components with the most pronounced changes occurring in the regions close to the core. For NRAO 530, we estimated a P.A. swing of $3^{\circ}.4$ per year for the entire inner jet (components d and e). The non-ballistic motion and change of jet orientation makes this source another prominent example of a helical and possibly `swinging' jet.

Gamma-ray bursts (GRBs) are divided into two classes according to their durations. We investigate if the softness of bursts plays a role in the conventional classification of the objects. We employ the BATSE (Burst and Transient Source Experiment) catalog and analyze the duration distributions of different groups of GRBs associated with distinct softness. Our analysis reveals that the conventional classification of GRBs with the duration of bursts is influenced by the softness of the objects. There exits a bimodality in the duration distribution of GRBs for each group of bursts and the time position of the dip in the bimodality histogram shifts with the softness parameter. Our findings suggest that the conventional classification scheme should be modified by separating the two well-known populations in different softness groups, which would be more reasonable than doing so with a single sample. According to the relation between the dip position and the softness parameter, we get an empirical function that can roughly set apart the short-hard and long-soft bursts: $SP = (0.100 \pm 0.028) T_{90}^{-(0.85 \pm 0.18)}$, where $SP$ is the softness parameter adopted in this paper.

The compact radio source Sgr\,A*, associated with the super massive black hole at the center of the Galaxy, has been studied with VLBA observations at 3 frequencies (22, 43, 86\,GHz) performed on 10 consecutive days in May 2007. The total VLBI flux density of Sgr\,A* varies from day to day. The variability is correlated at the 3 observing frequencies with higher variability amplitudes appearing at the higher frequencies. For the modulation indices, we find 8.4\,% at 22\,GHz, 9.3\,% at 43\,GHz, and 15.5\,% at 86\,GHz. The radio spectrum is inverted between 22 and 86\,GHz, suggesting inhomogeneous synchrotron self-absorption with a turnover frequency at or above 86\,GHz. The radio spectral index correlates with the flux density, which is harder (more inverted spectrum) when the source is brighter. The average source size does not appear to be variable over the 10-day observing interval. However, we see a tendency for the sizes of the minor axis to increase with increasing total flux, whereas the major axis remains constant. Towards higher frequencies, the position angle of the elliptical Gaussian increases, indicative of intrinsic structure, which begins to dominate the scatter broadening. At cm-wavelength, the source size varies with wavelength as $\lambda^{2.12\pm0.12}$, which is interpreted as the result of interstellar scatter broadening. After removal of this scatter broadening, the intrinsic source size varies as $\lambda^{1.4 ... 1.5}$. The VLBI closure phases at 22, 43, and 86\,GHz are zero within a few degrees, indicating a symmetric or point-like source structure. In the context of an expanding plasmon model, we obtain an upper limit of the expansion velocity of about 0.1\,c from the non-variable VLBI structure. This agrees with the velocity range derived from the radiation transport modeling of the flares from the radio to NIR wavelengths.}

We report on a successful, simultaneous observation and modelling of the millimeter (mm) to near-infrared (NIR) flare emission of the Sgr A* counterpart associated with the supermassive black hole at the Galactic centre (GC). We present a mm/sub-mm light curve of Sgr A* with one of the highest quality continuous time coverages and study and model the physical processes giving rise to the variable emission of Sgr A*.

A time-resolved spectral analysis for a sample of 22 intense, broad GRB pulses from the BATSE GRB sample is presented. We fit the spectra with the Band function and investigate the correlation between the observed flux (F) and the peak energy (E_p) of the $\nu f_\nu$ spectrum in the rising and decaying phases of these pulses. Two kinds of E_p evolution trends, i.e., hard-to-soft (the two-third pulses in our sample) and $E_{\rm p}$-tracing-$F$ (the one-third pulses in our sample) are observed in pulses from different GRBs and even from different pulses of the same burst. No dependence of spectral evolution feature on the pulse shape is found. A tight $F-E_{\rm p}$ positive correlation is observed in the decaying phases, with a power-law index $\sim 2.2$, which is much shallower than that expectation of the curvature effect. In the rising phase, the observed $F$ is either correlated or anti-correlated with $E_{\rm p}$, depending on the spectral evolution feature, and the power-law index of the correlation is dramatically different among pulses. More than $80%$ of the low energy photon indices in the time-resolved spectra whose $E_{\rm p}$ is anti-correlated with $F$ during the rising phase violate the death line of the synchrotron radiation, disfavoring the synchrotron radiation model for these gamma-rays. The $F-E_{\rm p}$ correlation, especially for those GRBs with $E_{\rm p}$-tracking-$F$ spectral evolution, may be due to the viewing angle and jet structure effects. In this cenario, the observed $F-E_{\rm p}$ correlation in the rising phase may be due to the line of sight from off-beam to on-beam toward a structured jet (or jitter), and the decaying phase is contributed by both the on-beam emission and the decayed photons from high latitude of the GRB fireball,resulting in a shallower slope of the observed $F-E_{\rm p}$ correlation than that predicted by the pure curvature effect.

Since the launch of Swift satellite, the detections of high-z (z>4) long gamma-ray bursts (LGRBs) have been rapidly growing, even approaching the very early Universe (the record holder currently is z=8.3). The observed high-z LGRB rate shows significant excess over that estimated from the star formation history. We investigate what may be responsible for this high productivity of GRBs at high-z through Monte Carlo simulations, with effective Swif/BAT trigger and redshift detection probabilities based on current Swift/BAT sample and CGRO/BATSE LGRB sample. We compare our simulations to the Swift observations via log N-log P, peak luminosity (L) and redshift distributions. In the case that LGRB rate is purely proportional to the star formation rate (SFR), our simulations poorly reproduce the LGRB rate at z>4, although the simulated log N-log P distribution is in good agreement with the observed one. Assuming that the excess of high-z GRB rate is due to the cosmic metallicity evolution or unknown LGRB rate increase parameterized as (1+z)^delta, we find that although the two scenarios alone can improve the consistency between our simulations and observations, incorporation of them gives much better consistency. We get 0.2<epsilon<0.6 and delta<0.6, where epsilon is the metallicity threshold for the production of LGRBs. The best consistency is obtained with a parameter set (epsilon, delta)=(~0.4, ~0.4), and BAT might trigger a few LGRBs at z~14. With increasing detections of GRBs at z>4 (~15% of GRBs in current Swift LGRB sample based on our simulations), a window for very early Universe is opening by Swift and up-coming SVOM missions.

Satellite data is accumulating that suggests and constrains dark matter physics. We argue there is a very well motivated theoretical preexisting framework consistent with dark matter annihilation being observed by the PAMELA satellite detector. The dark matter is (mainly) the neutral W boson superpartner, the wino with mass below 200 GeV. Using the program GALPROP we study the annihilation products and backgrounds together. Antimatter and gammas from annihilating winos contribute below this energy. We explain why PAMELA data does not imply no antiproton signal was observed by PAMELA or earlier experiments, and explain why the antiproton analysis was misunderstood by earlier papers. Wino annihilation does not describe the Fermi e+ + e- data (except partially below ~ 100 GeV). At higher energies we expect astrophysical mechanisms, and we simply parameterize them so the combination can describe all the data. We emphasize several predictions for satellite data to test the wino interpretation, particularly the turndown of the positron and antiproton spectra above 100 GeV. Most other interpretations require a large rise in the positron rates above 100 GeV. We focus on studying this well-motivated and long predicted wino interpretation, rather than comparisons with other interpretations. We emphasize that interpretations also depend very strongly on assumptions about the cosmological history of the universe, and on propagation of antiprotons and positrons in the galaxy. The winos PAMELA is observing arose from some non-thermal sources such as moduli decay rather than a universe that cooled in thermal equilibrium after the big bang. Then it is appropriate to normalize the wino density to the local relic density, and no "boost factors" are needed to obtain the reported PAMELA rates.

We report on a successful, simultaneous observation and modeling of the sub-millimeter to near-infrared flare emission of the Sgr A* counterpart associated with the super-massive black hole at the Galactic center. Our modeling is based on simultaneous observations that have been carried out on 03 June, 2008 using the NACO adaptive optics (AO) instrument at the ESO VLT and the LABOCA bolometer at the APEX telescope. Inspection and modeling of the light curves show that the sub-mm follows the NIR emission with a delay of 1.5+/-0.5 hours. We explain the flare emission delay by an adiabatic expansion of the source components.

The softening process observed in the steep decay phase of early X-ray afterglows of Swift bursts has remained a puzzle since its discovery. The softening process can also be observed in the later phase of the bursts and its cause has also been unknown. Recently, it was suggested that, influenced by the curvature effect, emission from high latitudes would shift the Band function spectrum from higher energy band to lower band, and this would give rise to the observed softening process accompanied by a steep decay of the flux density. The curvature effect scenario predicts that the terminating time of the softening process would be correlated with the duration of the process. In this paper, based on the data from the UNLV GRB group web-site, we found an obvious correlation between the two quantities. In addition, we found that the softening process can be divided into two classes: the early type softening ($t_{s,max}\leq "4000"s$) and the late type softening ($t_{s,max} > "4000"s$). The two types of softening show different behaviors in the duration vs. terminating time plot. In the relation between the variation rates of the flux density and spectral index during the softening process, a discrepancy between the two types of softening is also observed. According to their time scales and the discrepancy between them, we propose that the two types are of different origins: the early type is of internal shock origin and the late type is of external shock origin. The early softening is referred to the steep decay just following the prompt emission, whereas the late decay typically conceives the transition from flat decay to late afterglow decay. We suspect that there might be a great difference of the Lorentz factor in two classes which is responsible for the observed discrepancy.

Employing a sample presented by Kaneko et al. (2006) and Kocevski et al. (2003), we select 42 individual tracking pulses (here we defined tracking as the cases in which the hardness follows the same pattern as the flux or count rate time profile) within 36 Gamma-ray Bursts (GRBs) containing 527 time-resolved spectra and investigate the spectral hardness, $E_{peak}$ (where $E_{peak}$ is the maximum of the $\nu F_{\nu}$ spectrum), evolutionary characteristics. The evolution of these pulses follow soft-to-hard-to-soft (the phase of soft-to-hard and hard-to-soft are denoted by rise phase and decay phase, respectively) with time. It is found that the overall characteristics of $E_{peak}$ of our selected sample are: 1) the $E_{peak}$ evolution in the rise phase always start on the high state (the values of $E_{peak}$ are always higher than 50 keV); 2) the spectra of rise phase clearly start at higher energy (the median of $E_{peak}$ are about 300 keV), whereas the spectra of decay phase end at much lower energy (the median of $E_{peak}$ are about 200 keV); 3) the spectra of rise phase are harder than that of the decay phase and the duration of rise phase are much shorter than that of decay phase as well. In other words, for a complete pulse the initial $E_{peak}$ is higher than the final $E_{peak}$ and the duration of initial phase (rise phase) are much shorter than the final phase (decay phase). This results are in good agreement with the predictions of Lu et al. (2007) and current popular view on the production of GRBs. We argue that the spectral evolution of tracking pulses may be relate to both of kinematic and dynamic process even if we currently can not provide further evidences to distinguish which one is dominant. Moreover, our statistical results give some witnesses to constrain the current GRB model.

We report on recent near-infrared (NIR) and X-ray observations of Sagittarius A* (Sgr A*), the electromagnetic manifestation of the ~4x10^6 solar masses super-massive black hole (SMBH) at the Galactic Center. The goal of these coordinated multi-wavelength observations is to investigate the variable emission from Sgr A* in order to obtain a better understanding of the underlying physical processes in the accretion flow/outflow. The observations have been carried out using the NACO adaptive optics (AO) instrument at the European Southern Observatory's Very Large Telescope (July 2005, May 2007) and the ACIS-I instrument aboard the Chandra X-ray Observatory (July 2005). We report on a polarized NIR flare synchronous to a 8x1033 erg/s X-ray flare in July 2005, and a further flare in May 2007 that shows the highest sub-flare to flare contrast observed until now. The observations can be interpreted in the framework of a model involving a temporary disk with a short jet. In the disk component flux density variations can be explained due to hot spots on relativistic orbits around the central SMBH. The variations of the sub-structures of the May 2007 flare are interpreted as a variation of the hot spot structure due to differential rotation within the disk.

At the center of the Milky Way, with a distance of ~8 kpc, the compact source Sagittarius A* (SgrA*) can be associated with a super massive black hole of ~4x10^6 solar masses. SgrA* shows strong variability from the radio to the X-ray wavelength domains. Here we report on simultaneous NIR/sub-millimeter/X-ray observations from May 2007 that involved the NACO adaptive optics (AO) instrument at the European Southern Observatory's Very Large Telescope, the Australian Telescope Compact Array (ATCA), the US mm-array CARMA, the IRAM 30m mm-telescope, and other telescopes. We concentrate on the time series of mm/sub-mm data from CARMA, ATCA, and the MAMBO bolometer at the IRAM 30m telescope.

In May 2007 the compact radio source Sgr A* was observed in a global multi-frequency monitoring campaign, from radio to X-ray bands. Here we present and discuss first and preliminary results from polarization sensitive VLBA observations, which took place during May 14-25, 2007. Here, Sgr A* was observed in dual polarization on 10 consecutive days at 22, 43, and 86 GHz. We describe the VLBI experiments, our data analysis, monitoring program and show preliminary images obtained at the various frequencies. We discuss the data with special regard also to the short term variability.

We report on VLBA observations of a gamma-ray bright blazar NRAO 530 at multiple frequencies (5, 8, 15, 22, 39, 43, and 45 GHz) in 1997 and 1999. These multi-epoch multi-frequency high-resolution VLBI images exhibit a consistent core-dominated morphology with a bending jet to the north of the core. The quasi-simultaneous data observed at five frequencies (5, 8, 15, 22 and 43 GHz) in February 1997 enable us to estimate the spectra of compact VLBI components in this highly variable source. Flat spectra are seen in central two components (A and B), and the most compact component A with the flattest spectral index at the south end is identified as the core.

Using 64 ms count data for long gamma-ray bursts ($T_{90}>$2.6s), we analyze the quantity, relative spectral lag (RSL), which is defined as $\tau_{31}/FWHM_{(1)}$, where $\tau_{31}$ is the spectral lag between energy bands 1 and 3, and $FWHM_{(1)}$ denotes full width at half maximum of the pulse in energy channel 1. To get insights into features of the RSLs, we investigate in detail all the correlations between them and other parameters with a sub-sample including nine long bursts. The general cross-correlation technique is adopted to measure the lags between two different energy bands. We can derive the below conclusions. Firstly, the distribution of RSLs is normal and oncentrates on around the value of 0.1. Secondly, the RSLs are weakly correlated with $FWHM$, asymmetry, peak flux ($F_{p}$), peak energy ($E_{p}$) and spectral indexes ($\alpha$ and $\beta$), while they are uncorrelated with $\tau_{31}$, hardness-ratio ($HR_{31}$) and peak time ($t_m$). The final but important discovery is that redshift ($z$) and peak luminosity ($L_{p}$) are strongly correlated with the RSLs which can be measured easily and directly. We find that the RSL is a good redshift and peak luminosity estimator.

Recently, under the assumption that the Doppler effect of the relativistically expanding fireball surface is important, Qin et al. showed that in most cases the power law relationship between the pulse width and energy of gamma-ray bursts (GRBs)would exist in a certain energy range. We check this prediction with two GRB samples which contain well identified pulses. A power law anti-correlation between the full pulse width and energy and a power law correlation between the pulse width ratio and energy are seen in the light curves of the majority (around 65%) of bursts of the two samples within the energy range of BATSE, suggesting that these bursts are likely to arise from the emission associated with the shocks occurred on a relativistically expanding fireball surface. For the rest of the bursts, the relationships between these quantities were not predicted previously. We propose to consider other spectral evolutionary patterns or other radiation mechanisms such as a varying synchrotron or Comptonized spectrum to check if the observed relationships for these rest bursts can also be accounted for by the Doppler model. In addition, we find that the upper limits of the width ratio for the two samples do not exceed 0.9, in agrement with what predicted previously by the Doppler model. The plateau/power law/plateau and the peaked features predicted and detected previously by Qin et al. are generally observed, with the exceptions being noticed only in a few cases. According to the distinct values of two power law indices of FWHM and ratio and energy, we divide the bursts into three subsets which are located in different areas of the two indices plane. We suspect that different locations of the two indices might correspond to different mechanisms.

We present a detail analysis on the spectral lags of the short gamma-ray bursts (GRBs) and compare them with that of the long GRBs by using the CGRO/BATSE GRB Catalog. Our results are as follows. (1)The spectral lag distribution of the short GRBs is significantly different from that of the long GRBs. Excluding the statistical fluctuation effect, a proportion of ~17% of the short GRBs have a negative spectral lag, i.e., the hard photons being lag behind the soft photons. We do not find any peculiar features from their light curves to distinguish these bursts from those with a positive spectral lag. We argue that a more physical mechanism dominated the hard lag may be hid behind the morphological features of the light curves. This should be a great challenge to the current GRB models. We notice that this proportion is consistent with the proportion of short GRBs correlated with nearby galaxies newly discovered by Tanvir et al., although it is unclear if these short GRBs are indeed associated with the sources originated at low redshift. (2)While the spectral lags of the long GRBs are strongly correlated with the pulse durations, they are not for the short GRBs. However, the ratios of the spectral lag to the pulse duration for the short and long GRBs normally distributed at 0.023 and 0.046, respectively, with the sample width, indicating that the curvature effect alone could not explain the difference of the spectral lags between the two types of GRBs. The hydrodynamic timescales of the outflows and the radiative processes at work in GRBs might also play an important role as suggested by Daigne and Mochkovitch.

Recently, Shen et al. (2005) studied the contributions of curvature effect of fireballs to the spectral lag and showed that the observed lags could be accounted for by the effect. Here we check their results by performing a more precise calculation with both the formulas presented in Shen et al. (2005) and Qin et al. (2004). Several other aspects which were not considered in Shen et al. (2005) are investigated. We find that in the case of ultra-relativistic motions, both formulas are identical as long as the whole fireball surface is concerned. In our analysis, the previous conclusion that the detected spectral lags could be accounted for by the curvature effect is confirmed, while the conclusion that the lag has no dependence on the radius of fireballs is not true. We find that introducing extreme physical parameters is not the only outlet to explain those observed large lags. Even for the larger lags ($\sim5s$) (see Norris et al. 2005), a wider local pulse ($\Delta t_{\theta,FWHM} = 10^{7}s$) could account for it. Some conclusions not presented in Shen et al. (2005) or those modified in our analysis are listed below: a) lag $\propto \Gamma^{-\epsilon}$ with $\epsilon>$ 2; b) lag is proportional to the local pulse width and the $FWHM$ of the observed light curves; c) a large lag requires a large $\alpha_0$ and a small $\beta_0$ as well as a large $E_{0p}$; d) when the rest frame spectrum varies with time, the lag would become larger; e) lag decreases with the increasing of $ R_c$; f) lag $\propto$ E within the certain energy range for a given Lorentz factor; g) lag is proportional to the opening angle of uniform jets when $\theta_j < 0.6\Gamma^{-1}$.

In this paper, the effect of the intrinsic distribution of cosmological candles is investigated. We find that, in the case of a narrow distribution, the deviation of the observed modulus of sources from the expected central value could be estimated within a ceratin range. We thus introduce a lower and upper limits of $\chi ^{2}$, $\chi_{\min}^{2}$ and $ \chi_{\max}^{2}$, to estimate cosmological parameters by applying the conventional minimizing $\chi ^{2}$ method. We apply this method to a gamma-ray burst (GRB) sample as well as to a combined sample including this GRB sample and an SN Ia sample. Our analysis shows that: a) in the case of assuming an intrinsic distribution of candles of the GRB sample, the effect of the distribution is obvious and should not be neglected; b) taking into account this effect would lead to a poorer constraint of the cosmological parameter ranges. The analysis suggests that in the attempt of constraining the cosmological model with current GRB samples, the results tend to be worse than what previously thought if the mentioned intrinsic distribution does exist.

Kocevski et al. (2003) found that there is a linear relation between the rise width and the full width of gamma-ray burst pulses detected by the BATSE instrument based on their empirical functions. Motivated by this, we investigate the relationship based on Qin et al. (2004) model. Theoretical analysis shows that each of the two quantities, the rise width and the full width of observed pulse, is proportional to $\Gamma^{-2}\Delta\tau_{\theta,\rm FWHM}\frac{R_c}{c}$, where $\Gamma$ is the Lorentz factor for the bulk motion, $\Delta \tau_{\theta,\rm FWHM}$ is a local pulse's width, $R_{c}$ is the radius of fireballs and c is the velocity of light. We employ the observed pulses coming from four samples to study the relationship and find that: (1) Merely the curvature effect could produce the relationship with the same slope as those derived from Qin et al. (2004) model in the rise width vs. the full width panel. (2) Gamma-ray burst pulses, long or short ones (selected from the short and long GRBs), follow the same sequence in the rise width vs. the full width panel, with the shorter pulses at the end of this sequence. (3) All GRBs may intrinsically result from local Gaussian pulses. These features place constraints on the physical mechanism(s) for producing long and short GRBs.

Qin et al. (2004) have derived a formula of count rates based on a model of highly symmetric expanding fireballs, where the Doppler effect is the key factor to be concerned. In this paper, we employ the formula to both the Qin model and the uniform jet model to study how the rising timescale, $\Delta\tau_{\theta,r}$, and the decay timescale, $\Delta\tau_{\theta,d}$, of a local pulse affect the light curve of GRBs. Our analysis shows that they do make contributions to both rising and decay portions of the light curve of GRBs. Associated with a local pulse with both rising and decay portions, the light curve of GRBs in the rising portion is expected to undergo a concave phase and then a convex one, which we call a "concave-to-convex" character, whereas that in the decay portion is expected to evolve an opposite process which is called a "convex-to-concave" character, regardless of being in Qin model or in the uniform jet model. These characteristics are independent of local pulse shapes and their rest frame radiation forms. The diversity of light curves are caused by different forms of local pulses and different ratios of $\Delta\tau_{\theta,r}$ to $\Delta\tau_{\theta,d}$. We study a sample of 86 GRB pulses detected by the BATSE instrument on board the Compton Gamma Ray Observatory and find that the "concave-to-convex" (expected in the rising portion of the light curve) and "convex-to-concave" (expected in the decay portion) characters, which result from the curvature effect of spherical surface (or the Doppler effect), do exist in the light curve of some GRBs.

As revealed previously, under the assumption that some pulses of gamma-ray bursts are produced by shocks in spherical fireballs or uniform jets of large opening angles, there exists a standard decay form of the profile of pulses arising from very narrow or suddenly dimming local (or intrinsic) pulses due to the relativistic curvature effect (the Doppler effect over the spherical shell surface). Profiles of pulses arising from other local pulses were previously found to possess a reverse S-feature deviation from the standard decay form. We show in this paper that, in addition to the standard decay form shown in Qin et al. (2004), there exists a marginal decay curve associated with a local $\delta$ function pulse with a mono-color radiation. We employ the sample of Kocevski et al. (2003) to check this prediction and find that the phenomenon of the reverse S-feature is common, when compared with both the standard decay form and the marginal decay curve. We accordingly propose to take the marginal decay curve (whose function is simple) as a criteria to check if an observed pulse could be taken as a candidate suffered from the curvature effect. We introduce two quantities $A_1$ and $A_2$ to describe the mentioned deviations within and beyond the $FWHM$ position of the decay phase, respectively. The values of $A_1$ and $A_2$ of pulses of the sample are calculated, and the result suggests that for most of these pulses their corresponding local pulses might contain a long decay time relative to the time scale of the curvature effect.

We study in details how the pulse width of gamma-ray bursts is related with energy under the assumption that the sources concerned are in the stage of fireballs. Due to the Doppler effect of fireballs, there exists a power law relationship between the two quantities within a limited range of frequency. The power law range and the power law index depend strongly on the observed peak energy $E_p$ as well as the rest frame radiation form, and the upper and lower limits of the power law range can be determined by $E_p$. It is found that, within the same power law range, the ratio of the $FWHM$ of the rising portion to that of the decaying phase of the pulses is also related with energy in the form of power laws. A platform-power-law-platform feature could be observed in the two relationships. In the case of an obvious softening of the rest frame spectrum, the two power law relationships also exist, but the feature would evolve to a peaked one. Predictions on the relationships in the energy range covering both the BATSE and Swift bands for a typical hard burst and a typical soft one are made. A sample of FRED (fast rise and exponential decay) pulse bursts shows that 27 out of the 28 sources belong to either the platform-power-law-platform feature class or the peaked feature group, suggesting that the effect concerned is indeed important for most of the sources of the sample. Among these bursts, many might undergo an obvious softening evolution of the rest frame spectrum.

Assuming that the observed gamma-ray burst (GRB) rate as a function of redshift is proportional to a corrected star formation rate, we derive the empirical distribution of the viewing angles of long BATSE GRBs, $P^{\rm em}(\theta)$, and the distribution of these bursts in the plane of $\theta$ against redshift, $P^{\rm em}(\theta, z)$, by using a tight correlation between $E_{\gamma}$) and $E_{\rm p}^{'}$). Our results show that $P^{\rm em}(\theta)$ is well fitted by a log-normal distribution centering at $\log \theta/{\rm rad}=-0.76$ with a width of $\sigma_{\log \theta}=0.57$. We test different universal structured jet models by comparing model predictions with our empirical results. To make the comparisons reasonable, an "effective" threshold, which corresponds to the sample selection criteria of the long GRB sample, is used. We find that the predictions of a two-Gaussian jet model are roughly consistent with our empirical results. A brief discussion shows that cosmological effect on the $E_{\gamma}-E_{\rm p}^{'}$ relation does not significantly affect our results, but sample selection effects on this relationship might significantly influence our results.

The narrowness of the distribution of the peak energy of $\nu F_{\nu}$ spectrum of gamma-ray bursts (GRBs) and the unification of GRB population are great puzzles yet to be solved. We investigate the two puzzles based on the global spectral behaviors of different GRB population in the $HR-E_{\rm{p}}$ plane (HR the spectral hardness ratio) with BATSE and HETE-2 observations. It is found that long GRBs and XRFs observed by HETE-2 seem to follow the same sequence in the $HR-E_{\rm{p}}$ plane, with the XRFs at the low end of this sequence. The long and short GRBs observed by BATSE follow significantly different sequences in the $HR-E_{\rm p}$ plane, with most of the short GRBs having a larger hardness ratio than the long GRBs at a given $E_{\rm{p}}$. These results indicate that the global spectral behaviors of the long GRB sample and the XRF sample are similar, while that of short GRBs is different. The short GRBs seem to be a unique subclass of GRBs, and they are not the higher energy extension of the long GRBs (abridged).