We report the discovery of a Tidal Disruption Event (TDE) candidate occurring in the Active Galactic Nucleus SDSS J022700.77-042020.6. A sudden increase in flux of J0227-0420 during the second half of 2009 is clearly shown in the long-term optical, UV, and NIR light curves. A plateau phase, following an initial decline, is seen in the NUV and optical u, g, r, i light curves. Moreover, we find possible evidence that the plateau phase in the NUV band may lag behind the optical ones by approximately 70-80 days with also a much shorter duration, i.e. $\sim$7-15 days against $\sim$40-50 days. The long-term NUV/optical (after the plateau phase), NIR and MIR light curves can be well characterized with a form of $L(t)\propto t^{-\beta}$, consistent with the expectation of a TDE. The plateaus can be explained if the stellar streams collide with the pre-existing AGN disk at different radii. Though the overall fallback rate decreases, the material in the outer disk gradually drifts inward and increases the local accretion rate at the inner region, producing the optical and UV plateaus. The possible lag between the optical and NUV plateaus can then be attributed to viscosity delay. The index $\beta$ of the NIR $J, H, K_s$ bands ($\sim1.4-3.3$) is steeper than that of the UV/optical ($\sim0.7-1.3$) and MIR bands ($\sim0.9-1.8$), which may suggest that a certain fraction of the dust in the inner region of the dusty torus may be sublimated during the TDE phase. Our results indicate that, due to collisions between stellar debris and pre-existing disk, the light curves of TDEs occurring in AGN may show distinctive features, which may shed new light on the accretion process.

We report observations of a transient source \fermi\ from radio to \grs. \fermi\ was discovered by the {\it Fermi-LAT} in May 2017. Follow-up {\it Swift-XRT} observations revealed three flaring episodes through March 2018, and the peak X-ray flux is about $10^3$ higher than the {\it ROSAT all-sky survey (RASS)} flux upper limit. Optical spectral measurements taken by the {\it Magellan 6.5-m telescope} and the {\it Lick-Shane telescope} both show a largely featureless spectrum, strengthening the BL Lac interpretation first proposed by \citet{Bruni18}. The optical and mid-infrared (MIR) emission goes to a higher state in 2018, when the flux in high energies goes down to a lower level. Our {\it RATAN-600m} measurements at 4.8~GHz and 8.2~GHz do not indicate any significant radio flux variation over the monitoring seasons in 2017 and 2018, nor deviate from the archival {\it NVSS} flux level. During GeV flaring times, the spectrum is very hard ($\Gamma_\gamma\sim$1.7) in the GeV band and at times also very hard (($\Gamma_{\rm X}\lesssim2$) in the X-rays, similar to a high-synchrotron-peak (or even an extreme) BL Lac object, making \fermi\ a good target for ground-based {\it Cherenkov telescopes}.

We present and analyze the optical/UV and X-ray observations of a nearby tidal disruption event (TDE) candidate AT2019azh, spanning from 30 d before to ~ 250 d after its early optical peak. The X-rays show a late brightening by a factor of ~ 30-100 around 250 days after discovery, while the UV/opticals continuously decayed. The early X-rays show two flaring episodes of variation, temporally uncorrelated with the early UV/opticals. We found a clear sign of X-ray hardness evolution, i.e., the source is harder at early times, and becomes softer as it brightens later. The drastically different temporal behaviors in X-rays and UV/opticals suggest that the two bands are physically distinct emission components, and probably arise from different locations. These properties argue against the reprocessing of X-rays by any outflow as the origin of the UV/optical peak. The full data are best explained by a two-process scenario, in which the UV/optical peak is produced by the debris stream-stream collisions during the circularization phase; some low angular momentum, shocked gas forms an early, low-mass accretion disk which emits the early X-rays. The major body of the disk is formed after the circularization finishes, whose enhanced accretion rate produces the late X-ray brightening. AT2019azh is a strong case of TDE whose emission signatures of stream-stream collision and delayed accretion are both identified.

At about 70 solar masses, the recently-discovered dark object orbited by a B-type star in the system LB-1 is difficult to understand as the end point of standard stellar evolution, except as a binary black hole (BBH). LB-1 shows a strong, broad H-alpha emission line that is best attributed to a gaseous disk surrounding the dark mass. We use the observed H-alpha line shape, particularly its wing extension, to constrain the inner radius of the disk and thereby the separation of a putative BBH. The hypothesis of a current BBH is effectively ruled out on the grounds that its merger time must be a small fraction of the current age of the B star. The hypothesis of a previous BBH that merged to create the current dark mass is also effectively ruled out by the low orbital eccentricity, due to the combination of mass loss and kick resulted from gravitational wave emission in any past merger. We conclude that the current dark mass is a single black hole produced by the highly mass-conserving, monolithic collapse of a massive star.

All stellar mass black holes have hitherto been identified by X-rays emitted by gas that is accreting onto the black hole from a companion star. These systems are all binaries with black holes below 30 M$_{\odot}$$^{1-4}$. Theory predicts, however, that X-ray emitting systems form a minority of the total population of star-black hole binaries$^{5,6}$. When the black hole is not accreting gas, it can be found through radial velocity measurements of the motion of the companion star. Here we report radial velocity measurements of a Galactic star, LB-1, which is a B-type star, taken over two years. We find that the motion of the B-star and an accompanying H$\alpha$ emission line require the presence of a dark companion with a mass of $68^{+11}_{-13}$ M$_{\odot}$, which can only be a black hole. The long orbital period of 78.9 days shows that this is a wide binary system. The gravitational wave experiments have detected similarly massive black holes$^{7,8}$, but forming such massive ones in a high-metallicity environment would be extremely challenging to current stellar evolution theories$^{9-11}$.

Recently two fast X-ray transients (XT1 and XT2) have been reported from the search in the Chandra Deep Field (CDF) data. Each transient shows an initial plateau lasting around hundreds to thousands seconds followed by a rapid decay in the light curve. In particular, CDF-S XT2 is found to be associated with a galaxy at redshift z= 0.738 and was explained as a counterpart of a binary neutron-star merger event. In this paper, motivated by the short duration and decay slopes of the two transients, we consider an alternative interpretation in which both events are accretion-driven flares from tidal disruption of white dwarfs by intermediate-mass black holes. We derive a theoretical model of the accretion rate history, and find that it fits the observed X-ray light curves well. The extremely super-Eddington peak luminosity of XT2 can be explained by the beaming effect of the system, likely in the form of a jet.

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

A multi-messenger, black hole (BH) - neutron star (NS) merger event still remains to be detected. The tidal (dynamical) ejecta from such an event, thought to produce a kinonova, is concentrated in the equatorial plane and occupies only part of the whole azimuthal angle. In addition, recent simulations suggest that the outflow or wind from the post-merger remnant disk, presumably anisotropic, can be a major ejecta component responsible for a kilonova. For any ejecta whose photosphere shape deviates from the spherical symmetry, the electron scattering at the photosphere causes a net polarization in the kilonova light. Recent observational and theoretical polarization studies have been focused to the NS-NS merger kilonova AT2017gfo. We extend those work to the case of a BH-NS merger kilonova. We show that the degree of polarization at the first $\sim 1$ hr can be up to $\sim$ 3\% if a small amount ($10^{-4} M_{\odot}$) of free neutrons have survived in the fastest component of the dynamical ejecta, whose beta-decay causes a precursor in the kilonova light. The polarization degree can be $\sim$ 0.6\% if free neutrons survived in the fastest component of the disk wind. Future polarization detection of a kilonova will constrain the morphology and composition of the dominant ejecta component, therefore help to identify the nature of the merger.

Two X-ray sources were recently discovered by Irwin et al. in compact companions to elliptical galaxies to show ultra-luminous flares with fast rise (~ minute) and decay (~ hour), and with a peak luminosity ~10^{40-41} erg/s. Together with two other sources found earlier, they constitute a new type of fast transients which cannot be attributed to neutron stars but might be due to intermediate-mass black holes (IMBHs; 10^{2-4} M_sun). The flaring behavior is recurrent for at least two sources. If the flare represents a short period of accretion onto an IMBH during the periastron passage of a donor star on an eccentric (i.e., repeating) or parabolic (non-repeating) orbit, we argue that the flare's rise time corresponds to the duration during which the donor's tidally stripped mass joins a residual disk at the pericenter. This duration is in turn equal to three other time scales: the duration of stripping, the sound crossing time of the donor, and the circular orbit time at the pericenter radius. Only a white dwarf can have a sound crossing time as short as one minute. Therefore, the donor must be a white dwarf and it was stripped of ~10^{-10} M_sun upon each passage at several to tens of Schwarzschild radii from the IMBH. The flux decay corresponds to the viscous drainage of the supplied mass toward the hole. Aided with long-term X-ray monitoring, this type of fast transients would be an ideal target for next-generation gravitational wave detectors.

There has been suggestive evidence of intermediate-mass black holes (IMBHs; 10^{3-5} M_sun) existing in some globular clusters (GCs) and dwarf galaxies, but IMBHs as a population remain elusive. As a main-sequence star passes too close by an IMBH it might be tidally captured and disrupted. We study the long-term accretion and observational consequence of such tidal disruption events. The disruption radius is hundreds to thousands of the BH's Schwarzschild radius, so the circularization of the falling-back debris stream is very inefficient due to weak general relativity effects. Due to this and a high mass fallback rate, the bound debris initially goes through a ~10 yr long super-Eddington accretion phase. The photospheric emission of the outflow ejected during this phase dominates the observable radiation and peaks in the UV/optical bands with a luminosity of 10^42 erg/s. After the accretion rate drops below the Eddington rate, the bolometric luminosity follows the conventional t^{-5/3} power-law decay, and X-rays from the inner accretion disk start to be seen. Modeling the newly reported IMBH tidal disruption event candidate 3XMM J2150-0551, we find a general consistency between the data and predictions. The search for these luminous, long-term events in GCs and nearby dwarf galaxies could unveil the IMBH population.

According to the traditional scenario for core-collapse supernovae, the core of the collapsing star forms a neutron star and its gravitational energy release sends out a shockwave into the stellar envelope. However, in a significant number of numerical simulations, the shock stalls and the star cannot be exploded successfully, especially for a massive, compact star. We consider an alternative scenario that with mass fallback, the collapsing star forms a black hole in the center, surrounded by a dense, hot accretion disk, which blows out an intense outflow (wind). The kinetic energy of the wind may result in a successful stellar explosion. With an improved version of the formulism in Kohri et al. (2005) who studied neutron star accretion of minor fallback, we study this disk wind-driven explosion by calculating the accretion history for a suite of pre-SN stellar models with different initial surface rotational velocities, masses and metallicities, and by comparing the disk wind energy with the binding energy of the infalling stellar envelope. We show that the most promising models to be exploded successfully by this new channel are those relatively compact pre-SN stars with relatively low metallicities and not too low specific angular momenta. The total energies of the explosions are $\sim 10^{51-52}$ergs, and a more massive progenitor may produce a more energetic explosion.

The angular distribution of GRB jets is not yet clear. The observed luminosity of GRB 170817A is the lowest among all known short GRBs, which is best explained by the fact that our line of sight is outside of the jet opening angle, $\theta_{obs}$ $>$ $\theta_j$, where $\theta_{obs}$ is the angle between our line of sight and jet axis. Inferred by gravitational wave observations, as well as radio and X-ray afterglow modeling of GRB 170817A, it is likely that $\theta_{obs}$ $\sim$ 20$^{\circ}$ -- 28$^{\circ}$. In this work, we quantitatively consider two scenarios of angular energy distribution of GRB ejecta: top-hat jet, and structured jet with a power law index $s$. For top-hat jet model, we get a large $\theta_j$ (e.g., $\theta_{j}$$>$$10^{\circ}$), a rather high local (i.e., z$<$0.01) short GRB rate $\sim$8--15$\times10^3$ $Gpc^{-3}$$yr^{-1}$ (estimated 90$\sim$1850 $ Gpc^{-3}$$yr^{-1}$ in Fong et al. 2015) and an extremely high $E_{peak,0}$(on-axis,rest-frame)$>$7.5$\times10^4$keV($\sim$500keV for typical short GRB). For structured jet model, we use $\theta_{obs}$ to give limits on $s$ and $\theta_{j}$ for a typical on-axis luminosity short GRB (e.g., $10^{49}$ erg s$^{-1}$$\sim$$10^{51}$ erg s$^{-1}$), and a low on-axis luminosity case (e.g., $10^{49}$ erg s$^{-1}$) gives more reasonable values of $s$. A structured jet model is more feasible for GRB 170817A than a top-hat jet model, due to the rather high local short GRB rate and the extremely high on-axis $E_{peak,0}$ almost rules out the top-hat jet model. GRB 170817A is likely a low on-axis luminosity GRB ($10^{49}$ erg s$^{-1}$) with a structured jet.

Recent observation of some luminous transient sources with low color temperatures suggests that the emission is dominated by optically thick winds driven by super-Eddington accretion. We present a general analytical theory of the dynamics of radiation pressure-driven, optically thick winds. Unlike the classical adiabatic stellar wind solution whose dynamics are solely determined by the sonic radius, here the loss of the radiation pressure due to photon diffusion also plays an important role. We identify two high mass loss rate regimes ($\dot{M} > L_{\rm Edd\,}/c^2$). In the large total luminosity regime the solution resembles an adiabatic wind solution. Both the radiative luminosity, $L$, and the kinetic luminosity, $L_k$, are super-Eddington with $L < L_k$ and $L \propto L_k^{1/3}$. In the lower total luminosity regime most of the energy is carried out by the radiation with $L_k < L \approx L_{\rm Edd\,}$. In a third, low mass loss regime ($\dot{M} < L_{\rm Edd\,}/c^2$), the wind becomes optically thin early on and, unless gas pressure is important at this stage, the solution is very different from the adiabatic one. The results are independent from the energy generation mechanism at the foot of the wind, therefore they are applicable to a wide range of mass ejection systems, from black hole accretion, to planetary nebulae, and to classical novae.

The nature of ultra-luminous X-ray sources (ULXs) has long been plagued by an ambiguity about whether the central compact objects are intermediate-mass (IMBH, >~ 10^3 M_sun) or stellar-mass (a few tens M_sun) black holes (BHs). The high luminosity (~ 10^39 erg/s) and super-soft spectrum (T ~ 0.1 keV) during the high state of the ULX source X-1 in the galaxy M101 suggest a large emission radius (>~ 10^9 cm), consistent with being an IMBH accreting at a sub-Eddington rate. However, recent kinematic measurement of the binary orbit of this source and identification of the secondary as a Wolf-Rayet star suggest a stellar-mass BH primary with a super-Eddington accretion. If that is the case, a hot, optically thick outflow from the BH can account for the large emission radius and the soft spectrum. By considering the interplay of photons' absorption and scattering opacities, we determine the radius and mass density of the emission region of the outflow and constrain the outflow mass loss rate. The analysis presented here can be potentially applied to other ULXs with thermally dominated spectra, and to other super-Eddington accreting sources.

During a stellar tidal disruption event (TDE), an accretion disk forms as stellar debris returns to the disruption site and circularizes. Rather than being confined within the circularizing radius, the disk can spread to larger radii to conserve angular momentum. A spreading disk is a source of matter for re-accretion at rates which can exceed the later stellar fallback rate, although a disk wind can suppress its contribution to the central black hole accretion rate. A spreading disk is detectible through a break in the central accretion rate history, or, at longer wavelengths, by its own emission. We model the evolution of TDE disk size and accretion rate, by accounting for the time-dependent fallback rate, for the influence of wind losses in the early, advective stage, and for the possibility of thermal instability for accretion rates intermediate between the advection-dominated and gas-pressure dominated states. The model provides a dynamic basis for modeling TDE light curves. All or part of a young TDE disk will precess as a solid body due to Lense-Thirring effect, and precession may manifest itself as quasi-periodic modulation of light curve. The precession period increases with time. Applying our results to the jetted TDE candidate Swift J1644+57, whose X-ray light curve shows numerous quasi-periodic dips, we argue that the data best fit a scenario in which a main-sequence star was fully disrupted by an intermediate mass black hole on an orbit significantly inclined from the black hole equator, with the apparent jet shutoff at t= 500 d corresponding to a disk transition from the advective state to the gas-pressure dominated state.

GRB 120422A is a low-luminosity Gamma-ray burst (GRB) associated with a bright supernova, which distinguishes itself by its relatively short T90 ~ 5 s and an energetic X-ray tail. We analyze the Swift BAT and XRT data and discuss the physical implications. We show that the early steep decline in the X-ray light curve can be interpreted as the curvature tail of a late emission episode around 58-86 s, with a curved instantaneous spectrum at the end of the emission episode. Together with the main activity in the first ~ 20 s and the weak emission from 40 s to 60 s, the prompt emission is variable, which points towards a central engine origin, in contrast to the shock breakout origin as invoked to interpret some other nearby low-luminosity supernova GRBs. The curvature effect interpretation and interpreting the early shallow decay as the coasting external forward shock emission in a wind medium both give a constraint on the bulk Lorentz factor \Gamma to be around several. Comparing the properties of GRB 120422A and other supernova GRBs, we found that the main criterion to distinguish engine-driven GRBs from the shock breakout GRBs is the time-averaged luminosity, with a separation line of about ~ 10^48 erg s^-1.

We present a comprehensive analysis of a bright, long duration (T90 ~ 257 s) GRB 110205A at redshift z= 2.22. The optical prompt emission was detected by Swift/UVOT, ROTSE-IIIb and BOOTES telescopes when the GRB was still radiating in the gamma-ray band. Nearly 200 s of observations were obtained simultaneously from optical, X-ray to gamma-ray, which makes it one of the exceptional cases to study the broadband spectral energy distribution across 6 orders of magnitude in energy during the prompt emission phase. By fitting the time resolved prompt spectra, we clearly identify, for the first time, an interesting two-break energy spectrum, roughly consistent with the standard GRB synchrotron emission model in the fast cooling regime. Although the prompt optical emission is brighter than the extrapolation of the best fit X/gamma-ray spectra, it traces the gamma-ray light curve shape, suggesting a relation to the prompt high energy emission. The synchrotron + SSC scenario is disfavored by the data, but the models invoking a pair of internal shocks or having two emission regions can interpret the data well. Shortly after prompt emission (~ 1100 s), a bright (R = 14.0) optical emission hump with very steep rise (alpha ~ 5.5) was observed which we interpret as the emission from the reverse shock. It is the first time that the rising phase of a reverse shock component has been closely observed. The full optical and X-ray afterglow lightcurves can be interpreted within the standard reverse shock (RS) + forward shock (FS) model. In general, the high quality prompt emission and afterglow data allow us to apply the standard fireball shock model to extract valuable information about the GRB including the radiation mechanism, radius of prompt emission R, initial Lorentz factor of the outflow, the composition of the ejecta, as well as the collimation angle and the total energy budget.

Observational evidence suggests a link between long duration gamma ray bursts (LGRBs) and Type Ic supernovae. Here, we propose a potential mechanism for Type Ic supernovae in LGRB progenitors powered solely by accretion energy. We present spherically-symmetric hydrodynamic simulations of the long-term accretion of a rotating gamma-ray burst progenitor star, a "collapsar," onto the central compact object, which we take to be a black hole. The simulations were carried out with the adaptive mesh refinement code FLASH in one spatial dimension and with rotation, an explicit shear viscosity, and convection in the mixing length theory approximation. Once the accretion flow becomes rotationally supported outside of the black hole, an accretion shock forms and traverses the stellar envelope. Energy is carried from the central geometrically thick accretion disk to the stellar envelope by convection. Energy losses through neutrino emission and nuclear photodisintegration are calculated but do not seem important following the rapid early drop of the accretion rate following circularization. We find that the shock velocity, energy, and unbound mass are sensitive to convective efficiency, effective viscosity, and initial stellar angular momentum. Our simulations show that given the appropriate combinations of stellar and physical parameters, explosions with energies ~5x10^50 ergs, velocities 3000 km/s, and unbound material masses >6 solar masses are possible in a rapidly rotating 16 solar mass main sequence progenitor star. Further work is needed to constrain the values of these parameters, to identify the likely outcomes in more plausible and massive LRGB progenitors, and to explore nucleosynthetic implications.

The plateaus observed in about one half of the early X-ray afterglows are the most puzzling feature in gamma-ray bursts (GRBs) detected by Swift. By analyzing the temporal and spectral indices of a large X-ray plateau sample, we find that 55% can be explained by external, forward shock synchrotron emission produced by a relativistic ejecta coasting in a \rho ~ r^{-2}, wind-like medium; no energy injection into the shock is needed. After the ejecta collects enough medium and transitions to the adiabatic, decelerating blastwave phase, it produces the post-plateau decay. For those bursts consistent with this model, we find an upper limit for the initial Lorentz factor of the ejecta, \Gamma_0 \leq 46 (\epsilon_e/0.1)^{-0.24} (\epsilon_B/0.01)^{0.17}; the isotropic equivalent total ejecta energy is E_{iso} ~ 10^{53} (\epsilon_e/0.1)^{-1.3} (\epsilon_B/0.01)^{-0.09} (t_b/10^4 s) erg, where \epsilon_e and \epsilon_B are the fractions of the total energy at the shock downstream that are carried by electrons and the magnetic field, respectively, and t_b is the end of the plateau. Our finding supports Wolf-Rayet stars as the progenitor stars of some GRBs. It raises intriguing questions about the origin of an intermediate-\Gamma_0 ejecta, which we speculate is connected to the GRB jet emergence from its host star. For the remaining 45% of the sample, the post-plateau decline is too rapid to be explained in the coasting-in-wind model, and energy injection appears to be required.

The association of long-duration gamma-ray bursts (LGRBs) with Type Ic supernovae presents a challenge to supernova explosion models. In the collapsar model for LGRBs, gamma rays are produced in an ultrarelativistic jet launching from the magnetosphere of the black hole that forms in the aftermath of the collapse of a rotating progenitor star. The jet is collimated along the star's rotation axis, but the concomitant luminous supernova should be relatively--though certainly not entirely--spherical, and should synthesize a substantial mass of 56Ni. Our goal is to provide a qualitative assessment of the possibility that accretion of the progenitor envelope onto the black hole, which powers the LGRB, could also deposit sufficient energy and nickel mass in the envelope to produce a luminous supernova. For this, the energy dissipated near the black hole during accretion must be transported outward, where it can drive a supernova-like shockwave. Here we suggest that the energy is transported by convection and develop an analytical toy model, relying on global mass and energy conservation, for the dynamics of stellar collapse. The model suggests that a ~10,000 km/s shock can be driven into the envelope and that ~10^51 erg explosions are possible. The efficiency with which the accretion energy is being transferred to the envelope is governed by the competition of advection and convection at distances ~100-1,000 km from the black hole and is sensitive to the values of the convective mixing length, the magnitude of the effective viscous stress, and the specific angular momentum of the infalling envelope. Substantial masses of 56Ni may be synthesized in the convective accretion flow over the course of tens of seconds from the initial circularization of the infalling envelope around the black hole. The synthesized nickel is convectively mixed with a much larger mass of unburned ejecta.

Gamma-ray burst X-ray flares are believed to mark the late time activity of the central engine. We compute the temporal evolution of the average flare luminosity $< L >$ in the common rest frame energy band of 44 GRBs taken from the large \emph{Swift} 5-years data base. Our work highlights the importance of a proper consideration of the threshold of detection of flares against the contemporaneous continuous X-ray emission. In the time interval $30 \rm{s}<t<1000\,\rm{s}$ we find $< L >\propto t^{-2.7\pm 0.1}$; this implies that the flare isotropic energy scaling is $E_{\rm{iso,flare}}\propto t^{-1.7}$. The decay of the continuum underlying the flare emission closely tracks the average flare luminosity evolution, with a typical flare to steep-decay luminosity ratio which is $L_{\rm{flare}}/L_{\rm{steep}}=4.7$: this suggests that flares and continuum emission are deeply related to one another. We infer on the progenitor properties considering different models. According to the hyper-accreting black hole scenario, the average flare luminosity scaling can be obtained in the case of rapid accretion ($t_{\rm{acc}}\ll t$) or when the last $\sim 0.5 M_{\sun}$ of the original $14 M_{\sun}$ progenitor star are accreted. Alternatively, the steep $\propto t^{-2.7}$ behaviour could be triggered by a rapid outward expansion of an accretion shock in the material feeding a convective disk. If instead we assume the engine to be a rapidly spinning magnetar, then its rotational energy can be extracted to power a jet whose luminosity is likely to be between the monopole ($L\propto e^{-2t}$) and dipole ($L\propto t^{-2}$) cases. In both scenarios we suggest the variability, which is the main signature of the flaring activity, to be established as a consequence of different kinds of instabilities.

GRB 090417B was an unusually long burst with a T_90 duration of at least 2130 s and a multi-peaked light curve at energies of 15-150 keV. It was optically dark and has been associated with a bright star-forming galaxy at a redshift of 0.345 that is broadly similar to the Milky Way. This is one of the few cases where a host galaxy has been clearly identified for a dark gamma-ray burst and thus an ideal candidate for studying the origin of dark bursts. We find that the dark nature of GRB 090417B cannot be explained by high redshift, incomplete observations, or unusual physics in the production of the afterglow. Assuming the standard relativistic fireball model for the afterglow we find that the optical flux is at least 2.5 mag fainter than predicted by the X-ray flux. The Swift/XRT X -ray data are consistent with the afterglow being obscured by a dense, localized sheet of dust approximately 30-80 pc from the burst along the line of sight. Our results suggest that this dust sheet imparts an extinction of A_V >~ 12 mag, which is sufficient to explain the missing optical flux. GRB 090417B is an example of a gamma-ray burst that is dark due to the localized dust structure in its host galaxy.

Late X-ray flares observed in X-ray afterglows of gamma-ray bursts (GRBs) suggest late central engine activities at a few minuets to hours after the burst. A few unambiguously confirmed cases of supernova associations with nearby long GRBs imply that an accompanying supernova-like component might be a common feature in all long GRB events. These motivate us to study the interactions of a late jet, responsible for a x-ray flare, with various components in a stellar explosion, responsible for a GRB. These components include a supernova shell-like ejecta, and a cocoon that was produced when the main jet producing the GRB itself was propagating through the progenitor star. We find that the interaction between the late jet and the supernova ejecta may produce a luminous (up to 10^49 erg s^-1) thermal X-ray transient lasting for ~ 10 s. The interaction between the late jet and the cocoon produces synchrotron self-absorbed non-thermal emission, with the observed peak X-ray flux density from 0.001 micro Jy to 1 mJy at 1 keV and a peak optical flux density from 0.01 micro Jy to 0.1 Jy (for a redshift z= 2). The light curve due to the late jet - cocoon interaction has very small pulse-width-to-time ratio, \Delta t / t \approx 0.01 - 0.5, where t is the pulse peak time since the burst trigger. Identifying these features in current and future observations would open a new frontier in the study of GRB progenitor stars.

We constrain the distance of the Gamma-Ray Burst (GRB) prompt emission site from the explosion centre, R, by determining the location of the electron's self absorption frequency in the GRB prompt optical-to-X/gamma-ray spectral energy distribution, assuming that the optical and the gamma-ray emissions are among the same synchrotron radiation continuum of a group of hot electrons. All possible spectral regimes are considered in our analysis. The method has only two assumed parameters, namely, the bulk Lorentz factor of the emitting source Gamma, and the magnetic field strength B in the emission region (with a weak dependence). We identify a small sample of 4 bursts that satisfy the following three criteria: (1) they all have simultaneous optical and gamma-ray detections in multiple observational time intervals; (2) they all show temporal correlations between the optical and gamma-ray light curves; and (3) the optical emission is consistent with belonging to the same spectral component as the gamma-ray emission. For all the time intervals of these 4 bursts, it is inferred that R \geq 10^{14} (Gamma/300)^{3/4} (B/10^5 Gauss)^{1/4} cm. For a small fraction of the sample, the constraint can be pinned down to R \approx 10^{14} - 10^{15} cm for Gamma ~ 300. For a second sample of bursts with prompt optical non-detections, only upper limits on R can be obtained. We find no inconsistency between the R-constraints for this non-detection sample and those for the detection sample.

A dust scattering model was recently proposed to explain the shallow X-ray decay (plateau) observed prevalently in Gamma-Ray Burst (GRB) early afterglows. In this model the plateau is the scattered prompt X-ray emission by the dust located close (about 10 to a few hundred pc) to the GRB site. In this paper we carefully investigate the model and find that the scattered emission undergoes strong spectral softening with time, due to the model's essential ingredient that harder X-ray photons have smaller scattering angle thus arrive earlier, while softer photons suffer larger angle scattering and arrive later. The model predicts a significant change, i.e., $\Delta \b \sim 2 - 3$, in the X-ray spectral index from the beginning of the plateau toward the end of the plateau, while the observed data shows close to zero softening during the plateau and the plateau-to-normal transition phase. The scattering model predicts a big difference between the harder X-ray light curve and the softer X-ray light curve, i.e., the plateau in harder X-rays ends much earlier than in softer X-rays. This feature is not seen in the data. The large scattering optical depths of the dust required by the model imply strong extinction in optical, $A_V \gtrsim $ 10, which contradicts current findings of $A_V= 0.1 - 0.7$ from optical and X-ray afterglow observations. We conclude that the dust scattering model can not explain the X-ray plateaus.

We investigate a scenario of photons scattering by electrons within a relativistic outflow. The outflow is composed of discrete shells with different speeds. One shell emits radiation for a short duration. Some of this radiation is scattered by the shell(s) behind. We calculate in a simple two-shell model the observed scattered flux density as a function of the observed primary flux density, the normalized arrival time delay between the two emission components, the Lorentz factor ratio of the two shells and the scattering shell's optical depth. Thomson scattering in a cold shell and inverse Compton scattering in a hot shell are both considered. The results of our calculations are applied to the Gamma-Ray Bursts and the afterglows. We find that the scattered flux from a cold slower shell is small and likely to be detected only for those bursts with very weak afterglows. A hot scattering shell could give rise to a scattered emission as bright as the X-ray shallow decay component detected in many bursts, on a condition that the isotropically equivalent total energy carried by the hot electrons is large, $\sim 10^{52-56}$ erg. The scattered emission from a faster shell could appear as a late short $\gamma$-ray/MeV flash or become part of the prompt emission depending on the delay of the ejection of the shell.

The Gamma-Ray Burst (GRB) prompt emission is believed to be from highly relativistic electrons accelerated in relativistic shocks. From the GRB high-energy power-law spectral indices $\beta$ observed by the Burst and Transient Source Experiment (BATSE) Large Area Detectors (LAD), we determine the spectral index, $p$, of the electrons' energy distribution. Both the theoretical calculations and numerical simulations of the particle acceleration in relativistic shocks show that $p$ has a universal value $\approx 2.2-2.3$. We show that the observed distribution of $p$ during GRBs is not consistent with a $\delta$-function distribution or an universal $p$ value, with the width of the distribution $\ge$ 0.54. The distributions of $p$ during X-ray afterglows are also investigated and found to be inconsistent with a $\delta$-function distribution. The $p$-distributions in blazars and pulsar wind nebulae are also broad, inconsistent with a $\delta$-function distribution.

We compute the temporal profiles of the gamma-ray burst pulse in the four Burst and Transient Source Experiment (BATSE) Large Area Detector (LAD) discriminator energy channels, with the relativistic curvature effect of a expanding fireball being explicitly investigated. Assuming an intrinsic 'Band' shape spectrum and an intrinsic energy-independent emission profile, we show that merely the curvature effect can produce detectable spectral lags if the intrinsic pulse profile has a gradually decaying phase. We examine the spectral lag's dependences on some physical parameters, such as the Lorentz factor \Gamma, the low-energy spectral index, \alpha, of the intrinsic spectrum, the duration of the intrinsic radiation t_d' and the fireball radius R. It is shown that approximately the lag \propto \Gamma ^{-1} and \propto t_d', and a spectrum with a more extruded shape (a larger \alpha) causes a larger lag. We find no dependence of the lag on R. Quantitatively, the lags produced from the curvature effect are marginally close to the observed ones, while larger lags require extreme physical parameter values, e.g., \Gamma < 50, or \alpha > -0.5. The curvature effect causes an energy-dependent pulse width distribution but the energy dependence of the pulse width we obtained is much weaker than the observed W \propto E^{-0.4} one. This indicates that some intrinsic mechanism(s), other than the curvature effect, dominates the pulse narrowing of GRBs.

We determined the characteristic variability time scales (\Delta t_p) of 410 bright and long GRBs, by locating the peaks of their Power Density Spectra, defined and calculated in the time domain. We found that the averaged variability time scale decreases with the peak flux. This is consistent with the time-dilation effect expected for the cosmological origin of GRBs. We also found that the occurrence distribution of the characteristic variability time scale shows bimodality, which might be interpreted as that the long GRB sample is composed of two sub-classes with different variability time scales. However, we found no difference for some other characteristics of these two sub-classes.

We have studied the averaged power density spectra (PDSs) of two classes of long-duration gamma-ray bursts in the recent classification by Balastegui et al.(2001) based on neural network analysis. Both PDSs follow a power law over a wide frequency range with approximately the same slope, which indicates that a process with a self-similar temporal property may underlie the emission mechanisms of both. The two classes of bursts are divided into groups according to their brightness and spectral hardness respectively and each group's PDS was calculated; For both classes, the PDS is found to flatten both with increasing burst brightness and with increasing hardness.