Can all breaks in gamma-ray burst afterglows be explained by jet effects?

Whether gamma-ray bursts are highly beamed or not is a very important question, as it has been pointed out that the beaming will lead to a sharp break in the afterglow light curves during the ultrarelativistic phase, with the breaking point determined by  Γ∼1/ θ ₀  , where Γ is the bulk Lorentz factor and θ ₀ is the initial half opening angle of the ejecta, and such a break is claimed to be present in the light curves of some GRBs. In this paper we will examine whether all the observed breaks in GRB afterglow light curves can be explained by jet effects. Here we present a detailed calculation of the jet evolution and emission, and have obtained a simple formula of bulk Lorentz factor evolution. We show that the light curves are very smoothly steepened by jet effect, and the shape of the light curve is determined by only one parameter –     , where E and n are the fireball energy and surrounding medium density, respectively. We find that for GRB 990123 and GRB 991216, the jet model can approximately fit their light curves, and the values of     are about 0.17 and 0.22, respectively. On the other hand, the light curves of GRB 990510, GRB 000301c, GRB 000926 and GRB 010222 cannot be fitted by the jet model, which suggests that the breaks may be caused by some other reasons, and the jet effect should be not the unique reason.     

Prompt optical emission and synchrotron self-absorption constraints on emission site of GRBs

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/γ-ray spectral energy distribution, assuming that the optical and the γ-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 Γ and the magnetic field strength B in the emission region (with a weak dependence). We identify a small sample of four bursts that satisfy the following three criteria: (1) they all have simultaneous optical and γ-ray detections in multiple observational time intervals, (2) they all show temporal correlations between the optical and γ-ray light curves and (3) the optical emission is consistent with belonging to the same spectral component as the γ-ray emission. For all the time intervals of these four bursts, it is inferred that   R ≥ 10¹⁴  (Γ/300)^3/4 ( B /10⁵ G)^1/4  cm. For a small fraction of the sample, the constraint can be pinned down to   R ≈ 10¹⁴–10¹⁵ cm  for  Γ∼ 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.     

High-energy afterglow emission from giant flares of soft gamma-ray repeaters: the case of the 2004 December 27 event from SGR 1806–20

We discuss the high-energy afterglow emission (including high-energy photons, neutrinos and cosmic rays) following the 2004 December 27 giant flare from the soft gamma-ray repeater (SGR) 1806−20. If the initial outflow is relativistic with a bulk Lorentz factor  Γ₀∼  tens, the high-energy tail of the synchrotron emission from electrons in the forward shock region gives rise to a prominent sub-GeV emission, if the electron spectrum is hard enough and if the initial Lorentz factor is high enough. This signal could serve as a diagnosis of the initial Lorentz factor of the giant flare outflow. This component is potentially detectable by the Gamma-Ray Large Area Telescope ( GLAST ) if a similar giant flare occurs in the GLAST era. With the available 10-MeV data, we constrain that  Γ₀ < 50  if the electron distribution is a single power law. For a broken power-law distribution of electrons, a higher Γ₀ is allowed. At energies higher than 1 GeV, the flux is lower because of a high-energy cut-off of the synchrotron emission component. The synchrotron self-Compton emission component and the inverse Compton scattering component off the photons in the giant flare oscillation tail are also considered, but they are found not significant given a moderate Γ₀ (e.g. ≤ 10). The forward shock also accelerates cosmic rays to the maximum energy 10¹⁷ eV, and generates neutrinos with a typical energy 10¹⁴ eV through photomeson interaction with the X-ray tail photons. However, they are too weak to be detectable.     

Particle Acceleration and High-Frequency (X-ray and γ-ray) Emission in the Jets of Active Galactic Nuclei

It is suggested that the outflowing plasma in the jets of active galactic nuclei (AGNs) is inhomogeneous and consists of separate clouds. These clouds are strongly magnetized and move away from the central engine at relativistic speeds. The clouds interact with an ambient medium which is assumed to be at rest. In the process of this interaction, particles of the ambient medium are accelerated to high energies at the cloud front and flow ahead of the front. It is shown that the radiation of the accelerated particles may be responsible for the X-ray and γ-ray emission from AGN jets. TeV γ-ray emission is generated in the inner parts of AGN jets where the Lorentz factor of the cloud fronts is Γ0≥ 30, while GeV γ-ray emission emanates from the outer parts of AGN jets where Γ0 is ∼ 10. This revised version was published online in July 2006 with corrections to the Cover Date.     

Combined Analysis of Ultraviolet and Radio Observations of the 7 May 2004 CME/Shock Event

We report results from the combined analysis of UV and radio observations of a CME-driven shock observed on 7 May 2004 above the southeast limb of the Sun at 1.86 R _⊙ with the Ultraviolet Coronagraph Spectrometer (UVCS) on board the Solar and Heliospheric Observatory (SOHO). The coronal mass ejection (CME) was first detected in white-light by the SOHO’s Large Angle and Spectrometric Coronagraph (LASCO) C2 telescope and shock-associated type II metric emission was recorded simultaneously by ground-based radio spectrographs. The shock speed (∼ 690 km s⁻¹), as deduced from the analysis of the type II emission drift in the radio spectra and the pre-shock local electron density estimated with the diagnostics provided by UVCS observations of the O vi λλ 1031.9, 1037.6 doublet line intensities, is just a factor ∼ 0.1 higher than the CME speed inferred by means of the white-light (and EUV) data in the middle corona. The local magnetosonic speed, computed from a standard magnetic field model, was estimated as high as ∼ 600 km s⁻¹, implying that the CME speed was probably just sufficient to drive a weak fast-mode MHD shock ahead of the front. Simultaneously with the type II radio emission, significant changes in the O vi doublet line intensities and profiles were recorded in the UVCS spectra and found compatible with abrupt post-shock plasma acceleration and modest ion heating. This work provides further evidence for the CME-driven origin of the shocks observed in the middle corona.     

Fitting X-ray Afterglow Light curves of Gamma-ray Bursts by Using the Magnetar Energy Injection Model

The central compact object for some gamma-ray bursts (GRBs) may be a strongly magnetized millisecond pulsar. It can inject energy to the outer shock of the GRB by through the magnetic dipole radiation, and therefore causes the shallow decay of the early afterglow. Recently, from a large number of GRB X-ray afterglows observed by Swift/XRT(X-ray telescope), it is revealed that many of them exhibit the shallow decay about 10²∼10⁴ s after the burst prompt emission. We have fitted the X-ray afterglow light curves of 11 GRBs by using the energy injection model of a magnetar with the rotation period in the millisecond order of magnitude. The obtained result shows the validity and universality of the magnetar energy injection model in explaining the shallow decay of afterglows, and simultaneously provides some constraints on the magnetic field strength and rotation period of the central magnetar.     

Kinematical diagrams for conical relativistic jets

We present diagrams depicting the expected inter-dependences of two key kinematical parameters of radio knots in the parsec-scale jets of blazars, deduced from VLBI observations. The two parameters are the apparent speed (υ _app = cβ _app) and the effective Doppler boosting factor (δ _eff) of the relativistically moving radio knot. A novel aspect of these analytical computations of β-δ diagrams is that they are made for parsec-scale jets having a conical shape, with modest opening angles (ω up to 10°), in accord with the VLBI observations of the nuclei of the nearest radio galaxies. Another motivating factor is the recent finding that consideration of a conical geometry can have important implications for the interpretation of a variety of radio observations of blazar jets. In addition to uniform jet flows (i.e., those having a uniform bulk Lorentz factor, Γ), computational results are also presented for stratified jets where an ultra-relativistic central spine along the jet axis is surrounded by a slower moving sheath, possibly arising from a velocity shear.     

Gamma-Ray Bursts and magnetar-forming Supernovae

The connection between long Gamma Ray Bursts (GRBs) and Supernovae (SNe) have been established through the well observed cases. These events can be explained as the prompt collapse to a black hole (BH) of the core of a massive star (M≳40M _⊙). The energies of these GRB-SNe were much larger than those of typical SNe, thus these SNe are called Hypernovae (HNe). The case of SN 2006aj/GRB060218 appears different: the GRB was weak and soft, being called an X-Ray Flash (XRF); the SN is dimmer and has very weak oxygen lines. The explosion energy of SN 2006aj was smaller, as was the ejected mass. In our model, the progenitor star had a smaller mass than other GRB-SNe (M∼20M _⊙), suggesting that a neutron star (NS) rather than a BH was formed. If the nascent NS was strongly magnetized as a magnetar and rapidly spinning, it may launch a weak GRB or an XRF. The peculiar light curve of Type Ib SN 2005bf may also be powered by a magnetar. The blue-shifted nebular emission lines of 2005bf indicate the unipolar explosion possibly related to standing accretion shock instability (SASI) associated with a newly born NS.     

A massive, dusty toroid with large grains in the pre-planetary nebula IRAS22036+5306

Using the Submillimeter Array (SMA), we have obtained high angular-resolution (∼1″) interferometric maps of the submillimeter (0.88 mm) continuum and CO J=3–2 line from IRAS 22036+5306 (I 22036), a bipolar pre-planetary nebula (PPN) with knotty jets discovered in our HST SNAPshot survey of young PPNe. In addition, we have obtained supporting lower-resolution (∼10″) 2.6 mm continuum and CO, ¹³CO J=1–0 observations with the Owens Valley Radio Observatory (OVRO) interferometer. We find an unresolved source of submillimeter (and millimeter-wave) continuum emission in I 22036, implying a very substantial mass (0.02–0.04M _⊙) of large (i.e., radius ≳1 mm), cold (≲50 K) dust grains associated with I 22036’s toroidal waist. The CO J=3–2 observations show the presence of a very fast (∼220 km s⁻¹), highly collimated, massive (0.03M _⊙) bipolar outflow with a very large scalar momentum (about 10³⁹ g cm s⁻¹), and the characteristic spatio-kinematic structure of bow-shocks at the tips of this outflow. The fast outflow in I 22036, as in most PPNe, cannot be driven by radiation pressure. The large mass of the torus suggests that it has most likely resulted from common-envelope evolution in a binary, however it remains to be seen whether or not the time-scales required for the growth of grains to millimeter sizes in the torus are commensurate with such a formation scenario. The presence of the torus should facilitate the formation of the accretion disk needed to launch the jet. We also find that the ¹³C/¹²C ratio in I 22036 is very high (0.16), close to the maximum value achieved in equilibrium CNO-nucleosynthesis (0.33). The combination of the high circumstellar mass (i.e., in the torus and an extended dust shell inferred from ISO far-infrared spectra) and the high ¹³C/¹²C ratio in I 22036 provides strong support for this object having evolved from a massive (≳4M _⊙) progenitor in which hot-bottom-burning has occurred.     

Acceleration of particles in the vicinity of a massive black hole

Recent results of the gamma-ray Cherenkov astronomy definitely prove the existence of fast variability in the very high energy (V.H.E.) gamma-ray flux of some active galactic nuclei. The BL Lac PKS 2155-304 for instance showed variations down to a few minutes time scale. From standard light travel time argument, these variations put extremely strong constraints on the size of the TeV emitting zone, which has to be of the order of a few Schwarzschild radius, even for high values of the relativistic Doppler factor of the emitting jets. Such discovery is a challenge for particle acceleration scenarios, which have to imagine efficient acceleration processes at work in a very compact zone. Eventually, the immediate vicinity of the central black hole appears as the most conservative choice for the location of the TeV emission region of active galactic nuclei. In this paper, we propose a two-step mechanism for charged particle acceleration in the magnetosphere of a massive black hole surrounded by an accretion disk. Particles first gain energy by a stochastic process during the accretion phase. It is shown that effective proton acceleration up to energies 10¹⁷–10¹⁹ eV is possible in a low-luminosity magnetized accretion disk with 2D turbulent motion. The distribution function of energetic protons over energies is a power law function with typical index ≃−1. Here electrons are not very efficiently accelerated because of their drastic losses by synchrotron radiation. In a second time, part of the fast particles escape from the disk and are then entrained by the magnetic structure above the disk, in the rotating black hole magnetosphere. They thus gain additional energy by direct centrifugal mechanism, up to about 10²⁰ eV for the protons and to 10–100 TeV for the electrons when they cross the light cylinder surface. Such energetic particles can further radiate in the TeV spectral range observed by Cherenkov experiments as HESS, MAGIC and VERITAS. Energetic protons can produce γ-radiation in the energy band 1 GeV–100 TeV and above mainly by nuclei collisions with the disk matter, clouds, or ambient low energy photons. Energetic electrons can also reach the required spectral range by inverse Compton emission. However their acceleration is less efficient due to heavy radiation losses, and only gained by centrifugal process during the second phase of the whole mechanism we describe. Our present analysis would therefore favor hadronic scenarios for TeV emission of active galactic nuclei. It is tempting to relate long term variability over years of TeV active galactic nuclei to the first stochastic acceleration phase, which also provides the needed power law particle distributions, while short term variability over minutes is more likely due to perturbations of the second fast direct acceleration phase.     

Gamma-ray bursts in the Swift era

GRBs are the most energetic combination of jets and disks in the Universe. Observations made using Swift reveal a complex temporal and spectral behaviour. We propose that this behaviour can be used to refine the GRB classification scheme and align it better with progenitor types. The early (prompt) X-ray light curve can be well described by an exponential which relaxes into a power law. The transition time between the exponential and the power law gives a well-defined timescale, T _p , for the burst duration which we use with the spectral index of the prompt emission, β _p , and the prompt power law decay index, α _p to define four classes of burst: short, slow, fast and soft. Short bursts tend to decline more gradually than long bursts. Most GRBs display a second “afterglow” component which can be fitted in a similar way to the early emission. During the decay of this second component, few GRBs show jet breaks in accord with pre-Swift predictions. However, the start time of the final afterglow decay, T _a , correlates with the peak of the prompt γ-ray emission spectrum, E _peak, in an analogous way to the Ghirlanda relation found between optical “jet-break” times, t _j , and E _peak. These data are inconsistent with simple achromatic jet-break models casting doubt on the reliability of using late temporal breaks to determine the jet collimation.     

Damped Oscillations of Coronal Loops

A mechanism of damped oscillations of a coronal loop is investigated. The loop is treated as a thin toroidal flux rope with two stationary photospheric footpoints, carrying both toroidal and poloidal currents. The forces and the flux-rope dynamics are described within the framework of ideal magnetohydrodynamics (MHD). The main features of the theory are the following: i) Oscillatory motions are determined by the Lorentz force that acts on curved current-carrying plasma structures and ii) damping is caused by drag that provides the momentum coupling between the flux rope and the ambient coronal plasma. The oscillation is restricted to the vertical plane of the flux rope. The initial equilibrium flux rope is set into oscillation by a pulse of upflow of the ambient plasma. The theory is applied to two events of oscillating loops observed by the Transition Region and Coronal Explorer (TRACE). It is shown that the Lorentz force and drag with a reasonable value of the coupling coefficient (c _d ) and without anomalous dissipation are able to accurately account for the observed damped oscillations. The analysis shows that the variations in the observed intensity can be explained by the minor radial expansion and contraction. For the two events, the values of the drag coefficient consistent with the observed damping times are in the range c _d ≈2 – 5, with specific values being dependent on parameters such as the loop density, ambient magnetic field, and the loop geometry. This range is consistent with a previous MHD simulation study and with values used to reproduce the observed trajectories of coronal mass ejections (CMEs).     

On-axis orphan afterglows

In many GRB inner engine models the highly relativistic GRB jets are engulfed by slower moving matter. This could result in different beaming for the prompt γ-ray emission and for the lower energy afterglow. In this case we should expect that some observer will see on-axis orphan afterglows: X-ray, optical and radio afterglows within the initial relativistic ejecta with no preceding GRB; the prompt γ-ray emission is pointing elsewhere. We show that the observations of the WFC on BeppoSAX constrain with high certainty the prompt X-ray beaming factor to be less than twice the prompt γ-ray beaming. The results of Ariel 5 are consistent with this interpretation. The RASS from ROSAT and HEAO-1 constrain the X-ray beaming factor at 400 and 20 min after the burst, respectively, to be comparable and certainly not much larger than the γ-ray beaming factor. There is no direct limit on the optical beaming. However, we show that observations of several months with existing hardware could result in a useful limit on the optical beaming factor of GRB afterglows.     

Coronal Field Opens at Lower Height During the Solar Cycles 22 and 23 Minimum Periods: IMF Comparison Suggests the Source Surface Should Be Lowered

The solar cycle 23 minimum period has been characterized by a weaker solar and interplanetary magnetic field. This provides an ideal time to study how the strength of the photospheric field affects the interplanetary magnetic flux and, in particular, how much the observed interplanetary fields of different cycle minima can be understood simply from differences in the areas of the coronal holes, as opposed to differences in the surface fields within them. In this study, we invoke smaller source surface radii in the potential-field source-surface (PFSS) model to construct a consistent picture of the observed coronal holes and the near-Earth interplanetary field strength as well as polarity measurements for the cycles 23 and 22 minimum periods. Although the source surface value of 2.5 R _⊙ is typically used in PFSS applications, earlier studies have shown that using smaller source surface heights generates results that better match observations during low solar activity periods. We use photospheric field synoptic maps from Mount Wilson Observatory (MWO) and find that the values of ≈ 1.9 R _⊙ and ≈ 1.8 R _⊙ for the cycles 22 and 23 minimum periods, respectively, produce the best results. The larger coronal holes obtained for the smaller source surface radius of cycle 23 somewhat offsets the interplanetary consequences of the lower magnetic field at their photospheric footpoints. For comparison, we also use observations from the Michelson Doppler Imager (MDI) and find that the source surface radius of ≈ 1.5 R _⊙ produces better results for cycle 23, rather than ≈ 1.8 R _⊙ as suggested from MWO observations. Despite this difference, our results obtained from MWO and MDI observations show a qualitative consistency regarding the origins of the interplanetary field and suggest that users of PFSS models may want to consider using these smaller values for their source surface heights as long as the solar activity is low.     

Estimating black hole masses in quasars using broad optical and UV emission lines

We review past work using broad emission lines as virial estimators of black hole masses in quasars. Basically one requires estimates of the emitting region radius and virial velocity dispersion to obtain black hole masses. The three major ways to estimate the broad-line emitting region (BLR) radius involve: (1) direct reverberation mapping, (2) derivation of BLR radius for larger samples using the radius-luminosity correlation derived from reverberation measures, and (3) estimates of BLR radius using the definition of the ionization parameter solved for BLR radius (photoionization method). At low redshift (z ? 0.7) FWHM Hβ serves as the most widely used estimator of virial velocity dispersion. FWHM Hβ can provide estimates for tens of thousands of quasars out to z ≈ 3.8 (IR spectroscopy beyond z ≈ 1). A new photoionization method also shows promise for providing many reasonable estimates of BLR radius via high S/N IR spectroscopy of the UV region 1300–2000 Å. FWHM Mgiiλ2800 can serve as a surrogate for FWHM Hβ in the range 0.4 ? z ? 6.5 while civλ1549 is affected by broadening due to non-virial motions and best avoided (i.e. there is no clear conversion factor between FWHM Hβ and FWHM Civλ1549). Most quasars yield mass estimates in the range 7 ? log M_BH ? 9.7. There is no strong evidence for values above 10.0 and there may be evidence for a turnover in the maximum black hole mass near z ≈ 5.     

Discovery of a Giant Radio Halo in a Massive Merging Cluster at z?=?0.443

We have discovered a giant radio halo in the massive merging cluster MACSJ0417.5-1154. This cluster, at a redshift of 0.443, is one of the most X-ray luminous galaxy cluster in the MAssive Cluster Survey (MACS) with an X-ray luminosity in the 0.1–2.4 keV band of 2.9×10⁴⁵ erg s^− 1. Recent observations from GMRT at 230 and 610 MHz have revealed a radio halo of ∼ 1.2 × 0.3 Mpc² in extent. This halo is elongated along the North-West, similar to the morphology of the X-ray emission from Chandra. The 1400 MHz radio luminosity (L _r) of the halo is ∼2 × 10²⁵ W Hz^− 1, in good agreement with the value expected from the L _x − L _r correlation for cluster halos.     

Electron-positron annihilation lines and decaying sterile neutrinos

If massive sterile neutrinos exist, their decays into photons and/or electron-positron pairs may give rise to observable consequences. We consider the possibility that MeV sterile neutrino decays lead to the diffuse positron annihilation line in the Milky Way center, and we thus obtain bounds on the sterile neutrino decay rate Γ _e ≥10⁻²⁸ s⁻¹ from relevant astrophysical/cosmological data. Also, we expect a soft gamma flux of 1.2×10⁻⁴–9.7×10⁻⁴ ph cm⁻² s⁻¹ from the Milky Way center which shows up as a small MeV bump in the background photon spectrum. Furthermore, we estimate the flux of active neutrinos produced by sterile neutrino decays to be 0.02–0.1 cm⁻² s⁻¹ passing through the earth.     

A systematic description of shocks in gamma-ray bursts – II. Simulation

In Paper I, we presented a detailed formulation of the relativistic shocks and synchrotron emission in the context of gamma-ray burst (GRB) physics. To see how well this model reproduces the observed characteristics of the GRBs and their afterglows, here we present the results of some simulations based on this model. They are meant to reproduce the prompt and afterglow emissions in some intervals of time during a burst. We show that this goal is achieved for both short and long GRBs and their afterglows, at least for part of the parameter space. Moreover, these results are evidence of the physical relevance of the two phenomenological models we have suggested in Paper I for the evolution of the active region – synchrotron emitting region in a shock. The dynamical active region model seems to reproduce the observed characteristics of prompt emissions and late afterglow better than the quasi-steady model which is more suitable for the onset of afterglows. Therefore, these simulations confirm the arguments presented in Paper I about the behaviour of these models based on their physical properties.     

Doppler boosting, superluminal motion, and the kinematics of AGN jets

We discuss results from a decade long program to study the fine-scale structure and the kinematics of relativistic AGN jets with the aim of better understanding the acceleration and collimation of the relativistic plasma forming AGN jets. From the observed distribution of brightness temperature, apparent velocity, flux density, time variability, and apparent luminosity, the intrinsic properties of the jets including Lorentz factor, luminosity, orientation, and brightness temperature are discussed. Special attention is given to the jet in M87, which has been studied over a wide range of wavelengths and which, due to its proximity, is observed with excellent spatial resolution. Most radio jets appear quite linear, but we also observe curved non-linear jets and non-radial motions. Sometimes, different features in a given jet appear to follow the same curved path but there is evidence for ballistic trajectories as well. The data are best fit with a distribution of Lorentz factors extending up to γ∼30 and intrinsic luminosity up to ∼10²⁶ W Hz⁻¹. In general, gamma-ray quasars may have somewhat larger Lorentz factors than non gamma-ray quasars. Initially the observed brightness temperature near the base of the jet extend up to ∼5×10¹³ K which is well in excess of the inverse Compton limit and corresponds to a large excess of particle energy over magnetic energy. However, more typically, the observed brightness temperatures are ∼2×10¹¹ K, i.e., closer to equipartition.