Adiabatic expansion, early X-ray data and the central engine in GRBs

The Swift satellite early X-ray data show a very steep decay in most of the gamma-ray bursts light curves. This decay is either produced by the rapidly declining continuation of the central engine activity or by some leftover radiation starting right after the central engine shuts off. The latter scenario consists of the emission from an 'ember' that cools via adiabatic expansion and, if the jet angle is larger than the inverse of the source Lorentz factor, the large angle emission. In this work, we calculate the temporal and spectral properties of the emission from such a cooling ember, providing a new treatment for the microphysics of the adiabatic expansion. We use the adiabatic invariance of   p ²_⊥/ B ( p _⊥  is the component of the electrons' momentum normal to the magnetic field, B ) to calculate the electrons' Lorentz factor during the adiabatic expansion; the electron momentum becomes more and more aligned with the local magnetic field as the expansion develops. We compare the theoretical expectations of the adiabatic expansion (and the large angle emission) with the current observations of the early X-ray data and find that only ∼20 per cent of our sample of 107 bursts are potentially consistent with this model. This leads us to believe that, for most bursts, the central engine does not turn off completely during the steep decay of the X-ray light curve; therefore, this phase is produced by the continued rapidly declining activity of the central engine.     

A Rigid-Field Hydrodynamics approach to modelling the magnetospheres of massive stars

We introduce a new Rigid-Field Hydrodynamics approach to modelling the magnetospheres of massive stars in the limit of very strong magnetic fields. Treating the field lines as effectively rigid, we develop hydrodynamical equations describing the one-dimensional flow along each, subject to pressure, radiative, gravitational and centrifugal forces. We solve these equations numerically for a large ensemble of field lines to build up a three-dimensional time-dependent simulation of a model star with parameters similar to the archetypal Bp star σ Ori E. Since the flow along each field line can be solved independently of other field lines, the computational cost of this approach is a fraction of an equivalent magnetohydrodynamical treatment.
The simulations confirm many of the predictions of previous analytical and numerical studies. Collisions between wind streams from opposing magnetic hemispheres lead to strong shock heating. The post-shock plasma cools initially via X-ray emission, and eventually accumulates into a warped, rigidly rotating disc defined by the locus of minima of the effective (gravitational plus centrifugal) potential. However, a number of novel results also emerge. For field lines extending far from the star, the rapid area divergence enhances the radiative acceleration of the wind, resulting in high shock velocities (up to  ∼3000 km s⁻¹  ) and hard X-rays. Moreover, the release of centrifugal potential energy continues to heat the wind plasma after the shocks, up to temperatures around twice those achieved at the shocks themselves. Finally, in some circumstances the cool plasma in the accumulating disc can oscillate about its equilibrium position, possibly due to radiative cooling instabilities in the adjacent post-shock regions.     

Jet breaks at the end of the slow decline phase of Swift GRB light curves

The Swift mission has discovered an intriguing feature of gamma-ray burst (GRBs) afterglows, a phase of shallow decline of the flux in the X-ray and optical light curves. This behaviour is typically attributed to energy injection into the burst ejecta. At some point this phase ends, resulting in a break in the light curve, which is commonly interpreted as the cessation of the energy injection. In a few cases, however, while breaks in the X-ray light curve are observed, optical emission continues its slow flux decline. This behaviour suggests a more complex scenario. In this paper, we present a model that invokes a double component outflow, in which narrowly collimated ejecta are responsible for the X-ray emission while a broad outflow is responsible for the optical emission. The narrow component can produce a jet break in the X-ray light curve at relatively early times, while the optical emission does not break due to its lower degree of collimation. In our model both components are subject to energy injection for the whole duration of the follow-up observations. We apply this model to GRBs with chromatic breaks, and we show how it might change the interpretation of the GRBs canonical light curve. We also study our model from a theoretical point of view, investigating the possible configurations of frequencies and the values of GRB physical parameters allowed in our model.     

Magnetic deformation of magnetars for the giant flares of the soft gamma-ray repeaters

We present one possible mechanism for the giant flares of the soft gamma-ray repeaters (SGRs) within the framework of the magnetar (superstrongly magnetized neutron star) model, motivated by the positive period increase associated with the August 27 event from SGR 1900+14. From second-order perturbation analysis of the equilibrium of the magnetic polytrope, we find that there exist different equilibrium states separated by the energy of the giant flares and the shift in the moment of inertia to cause the period increase. This suggests that, if we assume that global reconfiguration of the internal magnetic field of     suddenly occurs, the positive period increase     as well as the energy ≳10⁴⁴ erg of the giant flares may be explained. The moment of inertia can increase with a release of energy, because the star shape deformed by the magnetic field can be prolate rather than oblate. In this mechanism, since oscillation of the neutron star will be excited, a ∼ ms-period pulsation of the burst profile and an emission of gravitational waves are expected. The gravitational waves could be detected by planned interferometers such as LIGO, VIRGO and LCGT.     

Polarization evolution accompanying the very early sharp decline of gamma-ray burst X-ray afterglows

In the synchrotron radiation model, the polarization property depends on both the configuration of the magnetic field and the geometry of the visible emitting region. Some peculiar behaviours in the X-ray afterglows of gamma-ray bursts (GRBs) observed with Swift , such as energetic flares and a plateau followed by a sharp drop, might be highly linearly polarized because the outflows powering these behaviours may be dominated by Poynting flux. The breakdown of the symmetry of the visible emitting region may also be well hidden in the peculiar X-ray data and may give rise to interesting polarization signatures. In this paper, we focus on the polarization accompanying the very early sharp decline of GRB X-ray afterglows. We show that strong polarization evolution is possible in both the high latitude emission model and the dying central engine model, which are used to interpret this sharp X-ray decline. It is thus not easy to efficiently probe the physical origin of the very early X-ray sharp decline with future polarimetry. Strong polarization evolution is also possible in the decline phase of X-ray flares and in the shallow decline phase of X-ray light curves characterized by chromatic X-ray versus optical breaks. A detector such as the X-ray Telescope (XRT), but with polarization capability, on board a satellite like Swift would be suitable for testing our predictions.     

Spectral features in gamma-rays expected from millisecond pulsars

In the advent of next generation gamma-ray missions, we present general properties of spectral features of high-energy emission above 1 MeV expected for a class of millisecond, low magnetic field (∼10⁹ G) pulsars. We extend polar-cap model calculations of Rudak & Dyks by including inverse Compton scattering events in an ambient field of thermal X-ray photons and by allowing for two models of particle acceleration. In the range between 1 MeV and a few hundred GeV, the main spectral component is the result of curvature radiation of primary particles. The synchrotron component arising from secondary pairs becomes dominant only below 1 MeV. The slope of the curvature radiation spectrum in the energy range from 100 MeV to 10 GeV strongly depends on the model of longitudinal acceleration, whereas below ∼100 MeV all slopes converge to a unique value of 4/3 (in a ν ℱ _ν convention). The thermal soft X-ray photons, which come either from the polar cap or from the surface, are Compton upscattered to a very high energy domain and form a separate spectral component peaking at ∼1 TeV. We discuss the observability of millisecond pulsars by future high‐energy instruments and present two rankings relevant for GLAST and MAGIC. We point to the pulsar J0437−4715 as a promising candidate for observations.     

On the global structure of self-gravitating discs for softened gravity

The effects of gravitational softening on the global structure of self-gravitating discs in centrifugal equilibrium are examined in relation to hydrodynamical/gravitational simulations. The one-parameter spline softening proposed by Hernquist & Katz is used.
It is found that if the characteristic size of a disc, r , is comparable to or less than the gravitational softening length, ε, then the cross-section of the simulated disc is significantly larger than that of a no-softening (Newtonian) disc with the same mass and angular momentum.
We demonstrate, furthermore, that if r ≲ε/2 then the scaling relation r ∝ε^3/4 holds for a given mass and specific angular momentum distribution with mass. Finally, we compare some of the theoretical results obtained in this paper and a previous one with the results of numerical Tree-SPH simulations and find qualitative agreement.     

Supranova Events from Spun-up Neutron Stars: An Explosion in Search of an Observation

We consider a formation scenario for supramassive neutron stars (SMNSs) that takes place through mass and angular momentum transfer from a close companion during a low-mass X-ray binary phase, with the ensuing suppression of the magnetic field. After the end of the mass transfer phase, SMNSs will lose, through magnetic dipole radiation, most of their angular momentum, triggering the star's collapse to a black hole. We discuss the rate of occurrence of these collapses and propose that these stars, because of the baryon-clear environment in which the implosion/explosion takes place, are the originators of gamma-ray bursts.     

Tail emission of prompt gamma-ray burst jets

Tail emission of the prompt gamma-ray burst (GRB) is discussed using a multiple emitting sub-shell (inhomogeneous jet, sub-jets or mini-jets) model, where the whole GRB jet consists of many emitting sub-shells. One may expect that such a jet with angular inhomogeneity should produce spiky tail emission. However, we found that the tail is not spiky but is decaying roughly monotonically. The global decay slope of the tail is not so much affected by the local angular inhomogeneity but affected by the global sub-shell energy distribution. The fact that steepening GRB tail breaks appeared in some events prefers the structured jets. If the angular size of the emitting sub-shell is around 0.01–0.02 rad, some bumps or fluctuations appear in the tail emission observed frequently in long GRBs. If the parameter differences of sub-shell properties are large, the tail has frequent changes of the temporal slope observed in a few bursts. Therefore, the multiple emitting sub-shell model has the advantage of explaining the small-scale structure in the observed rapid decay phase.     

A unified treatment of the gamma-ray burst 021211 and its afterglow

The gamma-ray burst (GRB) 021211 had a simple light curve, containing only one peak and the expected Poisson fluctuations. Such a burst may be attributed to an external shock, offering the best chance for a unified understanding of the gamma-ray burst and afterglow emissions. We analyse the properties of the prompt (burst) and delayed (afterglow) emissions of GRB 021211 within the fireball model. Consistency between the optical emission during the first 11 min (which, presumably, comes from the reverse shock heating of the ejecta) and the later afterglow emission (arising from the forward shock) requires that, at the onset of deceleration (∼2 s), the energy density in the magnetic field in the ejecta, expressed as a fraction of the equipartition value  (ɛ _B )  , is larger than in the forward shock at 11 min by a factor of approximately 10³. We find that synchrotron radiation from the forward shock can account for the gamma-ray emission of GRB 021211; to explain the observed GRB peak flux requires that, at 2 s,  ɛ _B   in the forward shock is larger by a factor 100 than at 11 min. These results suggest that the magnetic field in the reverse shock and early forward shock is a frozen-in field originating in the explosion and that most of the energy in the explosion was initially stored in the magnetic field. We can rule out the possibility that the ejecta from the burst for GRB 021211 contained more than 10 electron–positron pairs per proton.     

Swift observations of GRB 050712

We present the results of X-ray and optical observations of GRB 050712 performed by Swift . The X-ray light curve of this burst exhibits episodes of flares in the first 1000 s, the same epoch at which the UVOT detected an optical counterpart. A shallow X-ray decay, with a decay slope of  α=−0.73  , followed and lasted ∼70 ks. This behaviour can be explained in terms of activity of the gamma-ray burst 'inner engine', with the possibility that the last flare is caused by the interaction of the ejecta with the surrounding medium.
We also find interesting spectral parameters for the X-ray emission. In particular, data suggest the presence of an intrinsic absorption in the first 1000 s, which can be explained if circumburst medium clouds lie along the line of sight.     

Formation of stars and planets: the role of magnetic fields

Star formation is thought to be triggered by gravitational collapse of the dense cores of molecular clouds. Angular momentum conservation during the collapse results in the progressive increase of the centrifugal force, which eventually halts the inflow of material and leads to the development of a central mass surrounded by a disc. In the presence of an angular momentum transport mechanism, mass accretion onto the central object proceeds through this disc, and it is believed that this is how stars typically gain most of their mass. However, the mechanisms responsible for this transport of angular momentum are not well understood. Although the gravitational field of a companion star or even gravitational instabilities (particularly in massive discs) may play a role, the most general mechanisms are turbulence viscosity driven by the magnetorotational instability (MRI), and outflows accelerated centrifugally from the surfaces of the disc. Both processes are powered by the action of magnetic fields and are, in turn, likely to strongly affect the structure, dynamics, evolutionary path and planet-forming capabilities of their host discs. The weak ionisation of protostellar discs, however, may prevent the magnetic field from effectively coupling to the gas and shear and driving these processes. Here I examine the viability and properties of these magnetically-driven processes in protostellar discs. The results indicate that, despite the weak ionisation, the magnetic field is able to couple to the gas and shear for fluid conditions thought to be satisfied over a wide range of radii in these discs.     

Precessing jets interacting with interstellar material as the origin for the light curves of gamma-ray bursts

We present an internal shock model with external characteristics for explaining the complicated light curves of gamma-ray bursts. Shocks produce gamma-rays in the interaction between a precessing beam of relativistic particles and the interstellar medium. Each time the particle beam passes the same line of sight with the observer the interstellar medium is pushed outward. Subsequent interactions between the medium and the beam are delayed by the extra distance to be travelled for the particles before the shock can form. This results in a natural retardation and leads to an intrinsic asymmetry in the light curves produced for gamma-ray bursts. In addition, we account for the cooling of the electron–proton plasma in the shocked region, which gives rise to an exponential decay in the gamma-ray flux. The combination of these effects and the precessing jet of ultrarelativistic particles produces light curves that can be directly compared with observed gamma-ray burst light curves. We illustrate the model by fitting a number of observed gamma-ray bursts that are difficult to explain with only a precessing jet. We develop a genetic algorithm to fit several observed gamma-ray bursts with remarkable accuracy. We find that for different bursts the observed fluence, assuming isotropic emission, easily varies over four orders of magnitude from the energy generated intrinsically.     

A possible energy mechanism for cosmological γ-ray bursts

We suggest that an extreme Kerr black hole with a mass ∼10⁶ M_⊙, a dimensionless angular momentum     and a marginally stable orbital radius     located in a normal galaxy may produce a γ -ray burst (GRB) by capturing and disrupting a star. During the capture period, a transient accretion disc is formed and a strong transient magnetic field ∼     lasting for     may be produced at the inner boundary of the accretion disc. A large amount of rotational energy of the black hole is extracted and released in an ultrarelativistic jet with a bulk Lorentz factor Γ larger than 10³ via the Blandford–Znajek process. The relativistic jet energy can be converted into γ -radiation via an internal shock mechanism. The GRB duration should be the same as the lifetime of the strong transient magnetic field. The maximum number of sub-bursts is estimated to be     because the disc material is likely to break into pieces with a size about the thickness of the disc h at the cusp     The shortest risetime of the burst estimated from this model is ∼     The model GRB density rate is also estimated.     

The mass and angular momentum losses from spinars

Mestel's stellar wind theory is applied to estimate the mass and angular momentum losses for the recently proposed spinar model of quasars. If a spinar is uniformly rotating and has temperatures in the corotating regions of over a billion degrees, then it is found that all the rotational energy will be lost in 10³ yr. This time scale is much shorter than the total life time of the spinar model which is about 10⁶ yr. Hence it is suggested that the temperatures of the corotating regions of spinars must be less than a billion degrees.     

Activation of the Blandford–Znajek mechanism in collapsing stars

The collapse of massive stars may result in the formation of accreting black holes in their interiors. The accreting stellar matter may advect substantial magnetic flux on to the black hole and promote the release of its rotational energy via magnetic stresses (the Blandford–Znajek mechanism). In this paper we explore whether this process can explain the stellar explosions and relativistic jets associated with long gamma-ray bursts. In particular, we show that the Blandford–Znajek mechanism is activated when the rest mass–energy density of matter drops below the energy density of the magnetic field in the near vicinity of the black hole (within its ergosphere). We also discuss whether such a strong magnetic field is in conflict with the rapid rotation of the stellar core required in the collapsar model, and suggest that the conflict can be avoided if the progenitor star is a component of a close binary. In this case the stellar rotation can be sustained via spin-orbital interaction. In an alternative scenario the magnetic field is generated in the accretion disc, but in this case the magnetic flux through the black hole ergosphere is not expected to be sufficiently high to explain the energetics of hypernovae by the BZ mechanism alone. However, this energy deficit can be recovered via the additional power provided by the disc.     

Evidence for energy injection and a fine-tuned central engine at optical wavelengths in GRB 070419A

We present a comprehensive multiwavelength temporal and spectral analysis of the 'fast rise exponential decay' GRB 070419A. The early-time emission in the γ-ray and X-ray bands can be explained by a central engine active for at least 250 s, while at late times the X-ray light curve displays a simple power-law decay. In contrast, the observed behaviour in the optical band is complex (from 10² up to 10⁶ s). We investigate the light-curve behaviour in the context of the standard forward/reverse shock model; associating the peak in the optical light curve at ∼450 s with the fireball deceleration time results in a Lorenz factor  Γ≈ 350  at this time. In contrast, the shallow optical decay between 450 and 1500 s remains problematic, requiring a reverse shock component whose typical frequency is above the optical band at the optical peak time for it to be explained within the standard model. This predicts an increasing flux density for the forward shock component until   t ∼ 4 × 10⁶ s  , inconsistent with the observed decay of the optical emission from   t ∼ 10⁴ s  . A highly magnetized fireball is also ruled out due to unrealistic microphysic parameters and predicted light-curve behaviour that is not observed. We conclude that a long-lived central engine with a finely tuned energy injection rate and a sudden cessation of the injection is required to create the observed light curves, consistent with the same conditions that are invoked to explain the plateau phase of canonical X-ray light curves of γ-ray bursts.     

The possible role of r-modes in post-glitch relaxation of the Crab pulsar

The loss of angular momentum through gravitational radiation, driven by the excitation of r-modes, is considered for neutron stars that have rotation frequencies lower than the associated critical frequency. We find that for reasonable values of the initial amplitudes of such pulsation modes of the star, being excited at the event of a glitch in a pulsar, the total post-glitch losses correspond to a negligible fraction of the initial rise of the spin frequency in the case of Vela and older pulsars. However, for the Crab pulsar the same effect would result, within a few months, in a decrease in its spin frequency by an amount larger than its glitch-induced frequency increase. This could provide an explanation for the peculiar behaviour observed in the post-glitch relaxations of the Crab pulsar.     

A method for detecting gravitational waves coincident with gamma-ray bursts

The mechanism for gamma-ray bursters and the detection of gravitational waves (GWs) are two outstanding problems facing modern physics. Many models of gamma-ray bursters predict copious GW emission, so the assumption of an association between GWs and gamma-ray bursts (GRBs) may be testable with existing bar GW detector data. We consider Weber bar data streams in the vicinity of known GRB times and present calculations of the expected signal after co-addition of 1000 GW/GRBs that have been shifted to a common zero time. Our calculations are based on assumptions concerning the GW spectrum and the redshift distribution of GW/GRB sources that are consistent with current GW/GRB models. We discuss further possibilities of GW detection associated with GRBs in light of future bar detector improvements and suggest that co-addition of data from several improved bar detectors may result in detection of GWs (if the GW/GRB assumption is correct) on a time-scale comparable to the LIGO projects.     

A Test on Different Types of the Time Curve of Hardness Ratio of Gamma-Ray Bursts based on the Curvature Effect

We analyzed a sample of 66 gamma-ray bursts (GRBs) and statistically confirmed the prediction on the time curve of the hardness ratio of GRBs made by Qin et al. based on the curvature effect. In their analysis, GRB pulses are divided into three types according to the shape of their raw hardness ratio (RHR) time curves, defined as to include the background counts to the signal counts, so as to make use of counts within small time intervals. Of the three types, very hard sources exhibit a perfect pulse-like profile (type 1), hard bursts possess a pulse-like profile with a dip in the decay phase (type 2), and soft bursts show no pulse-like profile but have only a dipped profile (type 3). In terms of the conventional hardness ratio, type 3 sources are indeed generally softer than those of type 1 and type 2, in agreement with the prediction. We found that the minimum value of RHR is sensitive in distinguishing the different types. We propose that GRB pulses can be classified according to the minimum value of RHR and that the different type sources may be connected with different strengths of the shock or/and the magnetic field.