Source author record

Jonathan Granot

Jonathan Granot appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

41works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

41 published item(s)

preprint2023arXiv

Prompt GRB Polarization from Non-Axisymmetric Jets

Time-resolved linear polarization ($Π$) measurements of the prompt gamma-ray burst emission can reveal its dominant radiation mechanism. A widely considered mechanism is synchrotron radiation, for which linear polarization can be used to probe the jet's magnetic-field structure, and in turn its composition. In axisymmetric jet models the polarization angle (PA) can only change by $90^\circ$, as $Π$ temporarily vanishes. However, some time-resolved measurements find a continuously changing PA, which requires the flow to be non-axisymmetric in at least one out of its emissivity, bulk Lorentz factor or magnetic field. Here we consider synchrotron emission in non-axisymmetric jets, from an ultrarelativistic thin shell, comprising multiple radially-expanding mini-jets (MJs) or emissivity patches within the global jet, that yield a continuously changing PA. We explore a wide variety of possibilities with emission consisting of a single pulse or multiple overlapping pulses, presenting time-resolved and integrated polarization from different magnetic field configurations and jet angular structures. We find that emission from multiple incoherent MJs/patches reduces the net polarization due to partial cancellation in the Stokes plane. When these contain a large-scale ordered field in the plane transverse to the radial direction, $Π$ always starts near maximal and then declines over the single pulse or shows multiple highly polarized peaks due to multiple pulses. Observing $Π\lesssim40\%$ ($15\%$) integrated over one (several) pulse(s) will instead favor a shock-produced small-scale field either ordered in the radial direction or tangled in the plane transverse to it.

preprint2022arXiv

Gamma-Ray Bursts at TeV Energies: Theoretical Considerations

Gamma-ray bursts (GRBs) are the most luminous explosions in the Universe and are powered by ultra-relativistic jets. Their prompt $γ$-ray emission briefly outshines the rest of the $γ$-ray sky making them detectable from cosmological distances. It is followed by, and sometimes partially overlaps with, a similarly energetic but very broadband and longer-lasting afterglow emission. While most GRBs are detected below a few MeV, over a hundred were detected at high ($\gtrsim0.1\;$GeV) energies and several have now been observed up to tens of GeV with the \textit{Fermi} Large Area Telescope (LAT). A new electromagnetic window in the very high energy (VHE) domain ($\gtrsim0.1\;$TeV) was recently opened with the detection of afterglow emission in the $(0.1$\textendash$1)\,$TeV energy band by ground-based imaging atmospheric Cherenkov telescopes. The emission mechanism for the VHE spectral component is not fully understood, and its detection offers important constraints for GRB physics. This review provides a brief overview of the different leptonic and hadronic mechanisms capable of producing VHE emission in GRBs. The same mechanisms possibly give rise to the high-energy spectral component seen during the prompt emission of many \textit{Fermi}-LAT GRBs. Possible origins of its delayed onset and long duration, well into the afterglow phase, with implications for the emission region and relativistic collisionless shock physics are discussed. Key results for using GRBs as ideal probes for constraining models of extra-galactic background light and intergalactic magnetic fields, as well as for testing Lorentz invariance violation, are presented.

preprint2022arXiv

Robust Features of Off-Axis Gamma-Ray Burst Afterglow Lightcurves

The ultra-relativistic outflows powering gamma-ray bursts (GRBs) acquire angular structure through their interaction with external material. They are often characterized by a compact, nearly uniform narrow core (with half-opening angle $θ_{c,\{ε,Γ\}}$) surrounded by material with energy per unit solid angle ($ε=ε_cΘ_ε^{-a}$, where $Θ_{\{ε,Γ\}}=[1+θ^2/θ_{c,\{ε,Γ\}}^2]^{1/2}$) and initial specific kinetic energy ($Γ_0-1=[Γ_c-1]Θ_Γ^{-b}$) declining as power laws. Multi-wavelength afterglow lightcurves of off-axis jets (with viewing angle $θ_{\rm obs} > θ_c$) offer robust ways to constrain $a$, $b$ and the external density radial profile ($ρ\propto R^{-k}$), even while other burst parameters may remain highly degenerate. We extend our previous work on such afterglows to include more realistic angular structure profiles derived from three-dimensional hydrodynamic simulations of both long and short GRBs (addressing also jets with shallow angular energy profiles, whose emission exhibits unique evolution). We present afterglow lightcurves based on our parameterized power-law jet angular profiles for different viewing angles $θ_{\rm obs}$ and $k=\{0,1,2\}$. We identify a unique evolutionary power-law phase of the characteristic synchrotron frequencies ($ν_m$ and $ν_c$) that manifests when the lightcurve is dominated by emission sensitive to the angular structure of the outflow. We calculate the criterion for obtaining single or double peaked light-curves in the general case when $θ_{c,Γ}\neqθ_{c,ε}$. We emphasize how the shape of the lightcurve and the temporal evolution of $ν_m$ and $ν_c$ can be used to constrain the outflow structure and potentially distinguish between magnetic and hydrodynamic jets.

preprint2020arXiv

Afterglow Lightcurves from Misaligned Structured Jets

GRB 170817A/GW 170817 is the first GRB clearly viewed far from the GRB jet's symmetry axis. Its afterglow was densely monitored over a wide range of frequencies and times. It has been modeled extensively, primarily numerically, and although this endeavour was very fruitful, many of the underlying model parameters remain undetermined. We provide analytic modelling of GRB afterglows observed off-axis, considering jets with a narrow core (of half-opening angle $θ_{\rm c}$) and power-law wings in energy per unit solid angle ($ε=ε_cΘ^{-a}$ where $Θ=[1+(θ/θ_{\rm c})^2]^{1/2}$) and initial specific kinetic energy ($Γ_0-1=[Γ_{\rm c,0}-1]Θ^{-b}$), as well as briefly discuss Gaussian jets. Our study reveals qualitatively different types of lightcurves that can be viewed in future off-axis GRBs, with either single or double peaks, depending on the jet structure and the viewing angle. Considering the lightcurve shape rather than the absolute normalizations of times and / or fluxes, removes the dependence of the lightcurve on many of the highly degenerate burst parameters. This study can be easily used to determine the underlying jet structure, significantly reduce the effective parameter space for numerical fitting attempts and provide physical insights. As an illustration, we show that for GRB 170817A, there is a strong correlation between the allowed values of $Γ_{\rm c,0}$ and $b$, leading to a narrow strip of allowed solutions in the $Γ_{\rm c,0}$-$b$ plane above some minimal values $Γ_{\rm c,0}\gtrsim 40, b\gtrsim1.2$. Furthermore, the Lorentz factor of the material dominating the early lightcurve can be constrained by three independent techniques to be $Γ_{0}(θ_{\rm min,0})\approx5-7$.

preprint2020arXiv

Modelling the prompt optical emission of GRB 180325A: the evolution of a spike from the optical to gamma-rays

The transition from prompt to the afterglow emission is one of the most exciting and least understood phases in gamma-ray bursts (GRBs). Correlations among optical, X-ray and gamma-ray emission in GRBs have been explored, to attempt to answer whether the earliest optical emission comes from internal and/or external shocks. We present optical photometric observations of GRB 180325A collected with the TAROT and RATIR ground-based telescopes. These observations show two strong optical flashes with separate peaks at $\sim50\;$s and $\sim120\;$s, followed by a temporally extended optical emission. We also present X-rays and gamma-ray observations of GRB 180325A, detected by the Burst Alert Telescope (BAT) and X-ray Telescope (XRT), on the Neil Gehrels Swift observatory, which both observed a narrow flash at $\sim80\;$s. We show that the prompt gamma- and X-ray early emission shares similar temporal and spectral features consistent with internal dissipation within the relativistic outflow (e.g. by internal shocks or magnetic reconnection), while the early optical flashes are likely generated by the reverse shock that decelerates the ejecta as it sweeps up the external medium.

preprint2020arXiv

Spectral Modification of Magnetar Flares by Resonant Cyclotron Scattering

Spectral modification of energetic magnetar flares by resonant cyclotron scattering (RCS) is considered. During energetic flares, photons emitted from the magnetically-trapped fireball near the stellar surface should resonantly interact with magnetospheric electrons or positrons. We show by a simple thought experiment that such scattering particles are expected to move at mildly relativistic speeds along closed magnetic field lines, which would slightly shift the incident photon energy due to the Doppler effect. We develop a toy model for the spectral modification by a single RCS that incorporates both a realistic seed photon spectrum from the trapped fireball and the velocity field of particles, which is unique to the flaring magnetosphere. We show that our spectral model can be effectively characterized by a single parameter; the effective temperature of the fireball, which enables us to fit observed spectra with low computational cost. We demonstrate that our single scattering model is in remarkable agreement with Swift/BAT data of intermediate flares from SGR 1900+14, corresponding to effective fireball temperatures of $T_{\rm eff}=6$-$7$ keV, whereas BeppoSAX/GRBM data of giant flares from the same source may need more elaborate models including the effect of multiple scatterings. Nevertheless, since there is no standard physically-motivated model for magnetar flare spectra, our model could be a useful tool to study magnetar bursts, shedding light on the hidden properties of the flaring magnetosphere.

preprint2019arXiv

Constraining the magnetic field structure in collisionless relativistic shocks with a radio afterglow polarization upper limit in GW170817

Gamma-ray burst afterglows arise from relativistic collisionless shocks in which the postshock tangled magnetic field $\vec{B}$ is produced by the two-stream and/or Weibel instabilities on plasma skin-depth scales $(c/ω_p)$. The field is expected to be oriented predominantly within the shock plane ($B_{\perp}$; transverse to the shock normal, $\hat{n}_{\rm{sh}}$), and is often approximated to be completely within it ($B_\parallel\equiv\hat{n}_{\rm{sh}}\,\cdot\,\vec{B}=0$). Current 2D/3D particle-in-cell simulations are limited to short timescales and box sizes $\lesssim10^4(c/ω_p)\ll R/Γ_{\rm{sh}}$ much smaller than the shocked region's comoving width, and cannot probe the asymptotic downstream $\vec{B}$ structure. We constrain the latter using the linear polarization upper limit, $|Π|<12\%$, on the radio afterglow of GW170817/GRB170817A. Afterglow polarization depends on the jet's angular structure, our viewing angle, and the $\vec{B}$ structure. In GW170817/GRB170817A the latter can be tightly constrained since the former two are constrained from observations. We model $\vec{B}$ as an isotropic field in 3D that is stretched along $\hat{n}_{\rm{sh}}$ by a factor $ξ\equiv{}B_{\parallel}/B_{\perp}$, whose initial value $ξ_f\equiv{}B_{\parallel,f}/B_{\perp,f}$ describes the field that survives downstream on plasma scales $\ll{}R/Γ_{\rm{sh}}$. We calculate $Π(ξ_f)$ by integrating over the shocked volume for core-dominated structured jets, with a local Blandford-McKee self-similar radial profile. We find that independent of the exact jet structure, $\vec{B}$ has a finite, but initially sub-dominant, parallel component: $0.57\lesssimξ_f\lesssim0.89$, making it less anisotropic. While this motivates numerical studies of the asymptotic $\vec{B}$ structure in relativistic collisionless shocks, it may be consistent with turbulence amplified magnetic field.

preprint2019arXiv

Linear polarization in gamma-ray burst prompt emission

Despite being hard to measure, GRB prompt $γ$-ray emission polarization is a valuable probe of the dominant emission mechanism and the outflow's composition and angular structure. During the prompt emission the outflow is ultra-relativistic with Lorentz factors $Γ\gg1$. We describe in detail the linear polarization properties of various emission mechanisms: synchrotron radiation from different magnetic field structures (ordered: toroidal $B_{\rm tor}$ or radial $B_\parallel$, and random: normal to the radial direction $B_\perp$), Compton drag, and photospheric emission. We calculate the polarization for different GRB jet angular structures (e.g. top-hat, Gaussian, power-law) and viewing angles $θ_{\rm obs}$. Synchrotron with $B_\perp$ can produce large polarizations, up to $25\%\lesssimΠ\lesssim45\%$, for a top-hat jet but only for lines of sight just outside the jet's sharp edge. The same also holds for Compton drag, albeit with a slightly higher overall $Π$. Moreover, we demonstrate how $Γ$-variations during the GRB or smoother jet edges would significantly reduce $Π$. We construct a semi-analytic model for non-dissipative photospheric emission from structured jets. Such emission can produce up to $Π\lesssim15\%$ with reasonably high fluences, but this requires steep gradients in $Γ(θ)$. A polarization of $50\%\lesssimΠ\lesssim65\%$ can robustly be produced only by synchrotron emission from a transverse magnetic field ordered on angles $\gtrsim\!1/Γ$ around our line of sight (like a global toroidal field). Therefore, such a model would be strongly favored even by a single secure measurement within this range. We find that such a model would also be favored if $Π\gtrsim20\%$ is measured in most GRBs within a large enough sample, by deriving the polarization distribution for our different emission and jet models.

preprint2016arXiv

GRBs from Magnetic Reconnection: Variability and Robustness of Lightcurves

The dissipation mechanism that powers gamma-ray bursts (GRBs) remains uncertain almost half a century after their discovery. The two main competing mechanisms are the extensively studied internal shocks and the less studied magnetic reconnection. Here we consider GRB emission from magnetic reconnection accounting for the relativistic bulk motions that it produces in the jet's bulk rest frame. Far from the source the magnetic field is almost exactly normal to the radial direction, suggesting locally quasi-spherical thin reconnection layers between regions of oppositely directed magnetic field. We show that if the relativistic motions in the jet's frame are confined to such a quasi-spherical uniform layer, then the resulting GRB lightcurves are independent of their direction distribution within this layer. This renders previous results for a delta-function velocity-direction distribution (Beniamini & Granot 2016) applicable to a much more general class of reconnection models, which are suggested by numerical simulations. Such models that vary in their velocity-direction distribution differ mainly in the size of the bright region that contributes most of the observed flux at a given emission radius or observed time. The more sharply peaked this distribution, the smaller this bright region, and the stronger the lightcurve variability that may be induced by deviations from a uniform emission over the thin reconnection layer, which may be expected in a realistic GRB outflow. This is reflected both in the observed image at a given observed time and in the observer-frame emissivity map at a given emission radius, which are calculated here for three simple velocity-direction distributions.

preprint2016arXiv

Learning About the Magnetar Swift J1834.9-0846 from its Wind Nebula

The first wind nebula around a magnetar was recently discovered in X-rays around Swift~J1834.9$-$0846. We study this magnetar's global energetics and the properties of its particle wind or outflows. At a distance of $\sim4\;$kpc, Swift~J1834.9$-$0846 is located at the center of the supernova remnant (SNR) W41 whose radius is $\sim 19\;$pc, an order of magnitude larger than that of the X-ray nebula ($\sim2\;$pc). The association with SNR W41 suggests a common age of $\sim5-100\;$kyr, while its spin-down age is $4.9$~kyr. A small natal kick velocity may partly explain why a wind nebula was detected around this magnetar but not around other magnetars, most of which appear to have larger kick velocities and may have exited their birth SNR. We find that the GeV and TeV source detected by Fermi/LAT and H.E.S.S., respectively, of radius $\sim11\;$pc is most likely of hadronic origin. The dynamics and internal structure of the nebula are examined analytically to explain the nebula's current properties. Its size may naturally correspond to the diffusion-dominated cooling length of the X-ray emitting $e^+e^-$ pairs. This may also account for the spectral softening of the X-ray emission from the nebula's inner to outer parts. Analysis of the X-ray synchrotron nebula implies that (i) the nebular magnetic field is $\gtrsim 11\;μ$G (and likely $\lesssim30\;μ$G), and (ii) the nebula is not powered predominantly by the magnetar's quiescent spin-down-powered MHD wind, but by other outflows that contribute most of its energy. The latter are most likely associated with the magnetar's bursting activity, and possibly dominated by outflows associated with its past giant flares. The energy source for the required outflows cannot be the decay of the magnetar's dipole field alone, and is most likely the decay of its much stronger internal magnetic field.

preprint2016arXiv

Properties of GRB Lightcurves from Magnetic Reconnection

The energy dissipation mechanism within Gamma-Ray Burst (GRB) outflows, driving their extremely luminous prompt $γ$-ray emission is still uncertain. The leading candidates are internal shocks and magnetic reconnection. While the emission from internal shocks has been extensively studied, that from reconnection still has few quantitative predictions. We study the expected prompt-GRB emission from magnetic reconnection and compare its temporal and spectral properties to observations. The main difference from internal shocks is that for reconnection one expects relativistic bulk motions with Lorentz factors $Γ'\gtrsim\;$a few in the jet's bulk frame. We consider such motions of the emitting material in two anti-parallel directions (e.g. of the reconnecting magnetic-field lines) within an ultra-relativistic (with $Γ\gg1$) thin spherical reconnection layer. The emission's relativistic beaming in the jet's frame greatly affects the light-curves. For emission at radii $R_0<R<R_0+ΔR$ (with $Γ=\rm{const}$) the observed pulse width is $ΔT\sim(R_0/2cΓ^2)\max(1/Γ',\,ΔR/R_0)$, i.e. up to $\simΓ'$ times shorter than for isotropic emission in the jet's frame. We consider two possible magnetic reconnection modes: a quasi steady-state with continuous plasma flow into and out of the reconnection layer, and sporadic reconnection in relativistic turbulence that produces relativistic plasmoids. Both of these modes can account for many observed prompt-GRB properties: variability, pulse asymmetry, the very rapid declines at their end, and pulse evolutions that are either hard to soft (for $Γ'\lesssim2$) or intensity tracking (for $Γ'>2$). However, the relativistic turbulence mode is more likely to be relevant for the prompt sub-MeV emission and can naturally account also for the peak luminosity$\,$--$\,$peak frequency correlation.

preprint2015arXiv

Gamma-Ray Bursts as Sources of Strong Magnetic Fields

Gamma-Ray Bursts (GRBs) are the strongest explosions in the Universe, which due to their extreme character likely involve some of the strongest magnetic fields in nature. This review discusses the possible roles of magnetic fields in GRBs, from their central engines, through the launching, acceleration and collimation of their ultra-relativistic jets, to the dissipation and particle acceleration that power their $γ$-ray emission, and the powerful blast wave they drive into the surrounding medium that generates their long-lived afterglow emission. An emphasis is put on particular areas in which there have been interesting developments in recent years.

preprint2014arXiv

On the Composition of GRBs' Collapsar Jets

The duration distribution of long Gamma Ray Bursts reveals a plateau at durations shorter than ~20 s (in the observer frame) and a power-law decline at longer durations (Bromberg et al., 2012). Such a plateau arises naturally in the Collapsar model. In this model the engine has to operate long enough to push the jet out of the stellar envelope and the observed duration of the burst is the difference between the engine's operation time and the jet breakout time. We compare the jet breakout time inferred from the duration distribution (~10 s in the burst's frame) to the breakout time of a hydrodynamic jet (~10 s for typical parameters) and of a Poynting flux dominated jet with the same overall energy (<~1 s). As only the former is compatible with the duration of the plateau in the GRB duration distribution, we conclude that the jet is hydrodynamic during most of the time that its head is within the envelope of the progenitor star and around the time when it emerges from the star. This would naturally arise if the jet forms as a hydrodynamic jet in the first place or if it forms Poynting flux dominated but dissipates most of its magnetic energy early on within the progenitor star and emerges as a hydrodynamic jet.

preprint2014arXiv

Quasi-Periodic Oscillations in Short Recurring Bursts of the Soft-Gamma Repeater J1550-5418

The discovery of quasi-periodic oscillations (QPOs) in magnetar giant flares has opened up prospects for neutron star asteroseismology. The scarcity of giant flares makes a search for QPOs in the shorter, far more numerous bursts from Soft Gamma Repeaters (SGRs) desirable. In Huppenkothen et al (2013), we developed a Bayesian method for searching for QPOs in short magnetar bursts, taking into account the effects of the complicated burst structure, and have shown its feasibility on a small sample of bursts. Here, we apply the same method to a much larger sample from a burst storm of 286 bursts from SGR J1550-5418. We report a candidate signal at 260 Hz in a search of the individual bursts, which is fairly broad. We also find two QPOs at 93 Hz and one at 127 Hz, when averaging periodograms from a number of bursts in individual triggers, at frequencies close to QPOs previously observed in magnetar giant flares. Finally, for the first time, we explore the overall burst variability in the sample, and report a weak anti-correlation between the power-law index of the broadband model characterising aperiodic burst variability, and the burst duration: shorter bursts have steeper power law indices than longer bursts. This indicates that longer bursts vary over a broader range of time scales, and are not simply longer versions of the short bursts.

preprint2014arXiv

The dynamics of a highly magnetized jet propagating inside a star

The collapsar model explains the association of long duration gamma-Ray Bursts (GRBs) with stellar collapse. It involves a relativistic jet that forms at the core of a collapsing massive star. The jet penetrates the stellar envelope and the prompt GRB emission is produced once the jet is well outside the star. Most current models for generation of relativistic jets involve Poynting flux dominated outflows. We explore here the propagation of such a jet through a stellar envelope. The jet forms a bow shock around it. Energy dissipation at the head of this shock supplies energy to a cocoon that surrounds the jet. This cocoon exerts pressure on the jet and collimates it. While this description resembles the propagation of a hydrodynamic jets there are significant qualitative differences. Two Strong shocks, the reverse shock that slows down the hydrodynamic jet and the collimation shock that collimates it, cannot form within the Poynting flux dominated jet. As a result this jet moves much faster and dissipates much less energy while it crosses the stellar envelope. We construct here a simple analytic model that explores, self consistently, the jet-cocoon interaction and dynamics. Using this model we determine the properties of the jet, including its velocity, propagation time and shape.

preprint2013arXiv

Constraints on Lorentz Invariance Violation with Fermi-LAT Observations of Gamma-Ray Bursts

Some Quantum Gravity (QG) theories allow for a violation of Lorentz invariance (LIV), manifesting as a dependence of the velocity of light in vacuum on its energy. If such a dependence exists, then photons of different energies emitted together by a distant source will arrive at the Earth at different times. High-energy (GeV) transient emissions from distant astrophysical sources such as Gamma-ray Bursts (GRBs) and Active Galaxy Nuclei can be used to search for and constrain LIV. The Fermi collaboration has previously analyzed two GRBs in order to put constraints on the dispersion parameter in vacuum, and on the energy scale at which QG effects causing LIV may arise. We used three different methods on four bright GRBs observed by the Fermi-LAT to get more stringent and robust constraints. No delays have been detected and strong limits on the QG energy scale are derived: for linear dispersion we set tight constraints placing the QG energy scale above the Planck mass; a quadratic leading LIV effect is also constrained.

preprint2013arXiv

New Methods for Timing Analysis of Transient Events, Applied to Fermi/GBM Magnetar Bursts

In order to discern the physical nature of many gamma-ray sources in the sky, we must look not only in spectral and spatial dimensions, but also understand their temporal variability. However, timing analysis of sources with a highly transient nature, such as magnetar bursts, is difficult: standard Fourier techniques developed for long-term variability generally observed, for example, from AGN often do not apply. Here, we present newly developed timing methods applicable to transient events of all kinds, and show their successful application to magnetar bursts observed with Fermi/GBM. Magnetars are a prime subject for timing studies, thanks to the detection of quasi-periodicities in magnetar Giant Flares and their potential to help shed light on the structure of neutron stars. Using state-of-the art statistical techniques, we search for quasi-periodicities (QPOs) in a sample of bursts from Soft Gamma Repeater SGR J0501+4516 observed with Fermi/GBM and provide upper limits for potential QPO detections. Additionally, for the first time, we characterise the broadband variability behaviour of magnetar bursts and highlight how this new information could provide us with another way to probe these mysterious objects.

preprint2012arXiv

Broadband Spectral Investigations of SGR J1550-5418 Bursts

We present the results of our broadband spectral analysis of 42 SGR J1550-5418 bursts simultaneously detected with the Swift/X-ray Telescope (XRT) and the Fermi/Gamma-ray Burst Monitor (GBM), during the 2009 January active episode of the source. The unique spectral and temporal capabilities of the XRT Windowed Timing mode have allowed us to extend the GBM spectral coverage for these events down to the X-ray domain (0.5-10 keV). Our earlier analysis of the GBM data found that the SGR J1550-5418 burst spectra were described equally well with a Comptonized model or with two blackbody functions; the two models were statistically indistinguishable. Our new broadband (0.5 - 200 keV) spectral fits show that, on average, the burst spectra are better described with two blackbody functions than with the Comptonized model. Thus, our joint XRT/GBM analysis clearly shows for the first time that the SGR J1550-5418 burst spectra might naturally be expected to exhibit a more truly thermalized character, such as a two-blackbody or even a multi-blackbody signal. Using the Swift and RXTE timing ephemeris for SGR J1550-5418 we construct the distribution of the XRT burst counts with spin phase and find that it is not correlated with the persistent X-ray emission pulse phase from SGR J1550-5418. These results indicate that the burst emitting sites on the neutron star need not be co-located with hot spots emitting the bulk of the persistent X-ray emission. Finally, we show that there is a significant pulse phase dependence of the XRT burst counts, likely demonstrating that the surface magnetic field of SGR J1550-5418 is not uniform over the emission zone, since it is anticipated that regions with stronger surface magnetic field could trigger bursts more efficiently.

preprint2012arXiv

Detection of spectral evolution in the bursts emitted during the 2008-2009 active episode of SGR J1550 - 5418

In early October 2008, the Soft Gamma Repeater SGRJ1550 - 5418 (1E 1547.0 - 5408, AXJ155052 - 5418, PSR J1550 - 5418) became active, emitting a series of bursts which triggered the Fermi Gamma-ray Burst Monitor (GBM) after which a second especially intense activity period commenced in 2009 January and a third, less active period was detected in 2009 March-April. Here we analyze the GBM data all the bursts from the first and last active episodes. We performed temporal and spectral analysis for all events and found that their temporal characteristics are very similar to the ones of other SGR bursts, as well the ones reported for the bursts of the main episode (average burst durations \sim 170 ms). In addition, we used our sample of bursts to quantify the systematic uncertainties of the GBM location algorithm for soft gamma-ray transients to < 8 deg. Our spectral analysis indicates significant spectral evolution between the first and last set of events. Although the 2008 October events are best fit with a single blackbody function, for the 2009 bursts an Optically Thin Thermal Bremsstrahlung (OTTB) is clearly preferred. We attribute this evolution to changes in the magnetic field topology of the source, possibly due to effects following the very energetic main bursting episode.

preprint2012arXiv

Interaction of a highly magnetized impulsive relativistic flow with an external medium

Important astrophysical sources, such as gamma-ray bursts (GRBs) or tidal disruption events, are impulsive -- strongly varying with time. These outflows are likely highly magnetized near the central source, but their interaction with the external medium is not yet fully understood. Here I consider the combined impulsive magnetic acceleration of an initially highly magnetized shell of plasma and its deceleration by the external medium. I find four main dynamical regimes, that (for a given outflow) depend on the external density. (I) For small enough external densities the shell becomes kinetically dominated before it is significantly decelerated, thus reverting to the familiar unmagnetized "thin shell" case, which produces bright reverse shock emission that peaks well after the prompt GRB. (II) For larger external densities the shell remains highly magnetized and the reverse shock is strongly suppressed. It eventually transfers most of its energy through pdV work to the shocked external medium, whose afterglow emission peaks on a timescale similar to the prompt GRB duration. (III) For even larger external densities there is no initial impulsive acceleration phase. (IV) For the highest external densities the flow remains Newtonian.

preprint2012arXiv

On the lateral expansion of GRB jets

The dynamics of GRB jets during the afterglow phase have an important effect on the interpretation of their observations and for inferring key physical parameters such as their true energy and event rate. Semi-analytic models generally predict a fast lateral expansion, where the jet opening angle asymptotically grows exponentially with its radius. Numerical simulations, however, show a much more modest lateral expansion, where the jet retains memory of its initial opening angle for a very long time, and the flow remains non-spherical until it becomes sub-relativistic, and only then gradually approaches spherical symmetry. Here we suggest a new analytic model based on a new physically derived recipe for the lateral expansion. We also generalize the model by relaxing the common approximations of ultra-relativistic motion and a narrow jet opening angle. We find that the new analytic model fits much better the results of numerical simulations, mainly because it remains valid also in the mildly relativistic, quasi spherical regime. This model shows that for modest initial jet half-opening angles, θ_0, the outflow is not sufficiently ultra-relativistic when its Lorentz factor reaches Γ= 1/θ_0 and therefore the sideways expansion is rather slow, showing no rapid, exponential phase. On the other hand, we find that jets with an extremely narrow initial half-opening angle, of about θ_0 << 10^{-1.5} or so, which are still sufficiently ultra-relativistic at Γ= 1/θ_0, do show a phase of rapid, exponential lateral expansion. However, even such jets that expand sideways exponentially are still not spherical when they become sub-relativistic.

preprint2012arXiv

Scaling relations between numerical simulations and physical systems they represent

The dynamical equations describing the evolution of a physical system generally have a freedom in the choice of units, where different choices correspond to different physical systems that are described by the same equations. Since there are three basic physical units, of mass, length and time, there are up to three free parameters in such a rescaling of the units, $N_f \leq 3$. In Newtonian hydrodynamics, e.g., there are indeed usually three free parameters, $N_f = 3$. If, however, the dynamical equations contain a universal dimensional constant, such as the speed of light in vacuum $c$ or the gravitational constant $G$, then the requirement that its value remains the same imposes a constraint on the rescaling, which reduces its number of free parameters by one, to $N_f = 2$. This is the case, for example, in magneto-hydrodynamics (MHD) or special relativistic hydrodynamics, where $c$ appears in the dynamical equations and forces the length and time units to scale by the same factor, or in Newtonian gravity where the gravitational constant $G$ appears in the equations. More generally, when there are $N_{udc}$ independent (in terms of their units) universal dimensional constants, then the number of free parameters is $N_f = max(0,3-N_{udc})$. When both gravity and relativity are included, there is only one free parameter ($N_f = 1$, as both $G$ and $c$ appear in the equations so that $N_{udc} = 2$), and the units of mass, length and time must all scale by the same factor. The explicit rescalings for different types of systems are discussed and summarized here. Such rescalings of the units also hold for discrete particles, e.g. in N-body or particle in cell simulations. They are very useful when numerically investigating a large parameter space or when attempting to fit particular experimental results, by significantly reducing the required number of simulations.

preprint2012arXiv

The effects of sub-shells in highly magnetized relativistic flows

Astrophysical sources of relativistic jets or outflows, such as gamma-ray bursts (GRBs), active galactic nuclei (AGN) or micro-quasars, often show strong time variability. Despite such impulsive behavior, most models of these sources assume a steady state for simplicity. Here I consider a time-dependent outflow that is initially highly magnetized and divided into many well-separated sub-shells, as it experiences impulsive magnetic acceleration and interacts with the external medium. In AGN the deceleration by the external medium is usually unimportant and most of the initial magnetic energy is naturally converted into kinetic energy, leading to efficient dissipation in internal shocks as the sub-shells collide. Such efficient low-magnetization internal shocks can also naturally occur in GRBs, where the deceleration by the external medium can be important. A strong low-magnetization reverse shock can develop, and the initial division into sub-shells allows it to be relativistic and its emission to peak on the timescale of the prompt GRB duration (which is not possible for a single shell). Sub-shells also enable the outflow to reach much higher Lorentz factors that help satisfy existing constraints on GRBs from intrinsic pair opacity and from the afterglow onset time.

preprint2011arXiv

Burst and Persistent Emission Properties during the Recent Active Episode of the Anomalous X-ray Pulsar 1E 1841-045

Swift/BAT detected the first burst from 1E 1841-045 in May 2010 with intermittent burst activity recorded through at least July 2011. Here we present Swift and Fermi/GBM observations of this burst activity and search for correlated changes to the persistent X-ray emission of the source. The T90 durations of the bursts range between 18-140 ms, comparable to other magnetar burst durations, while the energy released in each burst ranges between (0.8 - 25)E38 erg, which is in the low side of SGR bursts. We find that the bursting activity did not have a significant effect on the persistent flux level of the source. We argue that the mechanism leading to this sporadic burst activity in 1E 1841-045 might not involve large scale restructuring (either crustal or magnetospheric) as seen in other magnetar sources.

preprint2011arXiv

Fermi/GBM Observations of SGR J0501+4516 Bursts

We present our temporal and spectral analyses of 29 bursts from SGR J0501+4516, detected with the Gamma-ray Burst Monitor onboard the Fermi Gamma-ray Space Telescope during the 13 days of the source activation in 2008 (August 22 to September 3). We find that the T90 durations of the bursts can be fit with a log-normal distribution with a mean value of ~ 123 ms. We also estimate for the first time event durations of Soft Gamma Repeater (SGR) bursts in photon space (i.e., using their deconvolved spectra) and find that these are very similar to the T90s estimated in count space (following a log-normal distribution with a mean value of ~ 124 ms). We fit the time-integrated spectra for each burst and the time-resolved spectra of the five brightest bursts with several models. We find that a single power law with an exponential cutoff model fits all 29 bursts well, while 18 of the events can also be fit with two black body functions. We expand on the physical interpretation of these two models and we compare their parameters and discuss their evolution. We show that the time-integrated and time-resolved spectra reveal that Epeak decreases with energy flux (and fluence) to a minimum of ~30 keV at F=8.7e-6 erg/cm2/s, increasing steadily afterwards. Two more sources exhibit a similar trend: SGRs J1550-5418 and 1806-20. The isotropic luminosity corresponding to these flux values is roughly similar for all sources (0.4-1.5 e40 erg/s).

preprint2011arXiv

Gamma-Ray Burst Dynamics and Afterglow Radiation from Adaptive Mesh Refinement, Special Relativistic Hydrodynamic Simulations

We report on the development of Mezcal-SRHD, a new adaptive mesh refinement, special relativistic hydrodynamics (SRHD) code, developed with the aim of studying the highly relativistic flows in Gamma-Ray Burst sources. The SRHD equations are solved using finite volume conservative solvers. The correct implementation of the algorithms is verified by one-dimensional (1D) shock tube and multidimensional tests. The code is then applied to study the propagation of 1D spherical impulsive blast waves expanding in a stratified medium with $ρ\propto r^{-k}$, bridging between the relativistic and Newtonian phases, as well as to a two-dimensional (2D) cylindrically symmetric impulsive jet propagating in a constant density medium. It is shown that the deceleration to non-relativistic speeds in one-dimension occurs on scales significantly larger than the Sedov length. This transition is further delayed with respect to the Sedov length as the degree of stratification of the ambient medium is increased. This result, together with the scaling of position, Lorentz factor and the shock velocity as a function of time and shock radius, is explained here using a simple analytical model based on energy conservation. The method used for calculating the afterglow radiation by post-processing the results of the simulations is described in detail. The light curves computed using the results of 1D numerical simulations during the relativistic stage correctly reproduce those calculated assuming the self-similar Blandford-McKee solution for the evolution of the flow. The jet dynamics from our 2D simulations and the resulting afterglow lightcurves, including the jet break, are in good agreement with those presented in previous works. Finally, we show how the details of the dynamics critically depend on properly resolving the structure of the relativistic flow.

preprint2011arXiv

Simulations of GRB Jets in a Stratified External Medium: Dynamics, Afterglow Lightcurves, Jet Breaks and Radio Calorimetry

The dynamics of GRB jets during the afterglow phase is most reliably and accurately modelled using hydrodynamic simulations. All published simulations, however, have considered only a uniform external medium, while a stratified external medium is expected around long duration GRB progenitors. Here we present simulations of the dynamics of GRB jets and the resulting afterglow emission for both uniform and stratified external media with $ρ\propto r^{-k}$ for k = 0, 1, 2. The simulations are performed in 2D using the special relativistic version of the Mezcal code. The dynamics for stratified external media are broadly similar to those derived for expansion into a uniform external medium. The jet half-opening angle start increasing logarithmically with time once the Lorentz factor drops below 1/theta_0. For larger k values the lateral expansion is faster at early times and slower at late times with the jet expansion becoming Newtonian and slowly approaching spherical symmetry over progressively longer timescales. We find that contrary to analytic expectations, there is a reasonably sharp jet break in the lightcurve for k = 2 although the shape of the break is affected more by the viewing angle than by the slope of the external density profile. Steeper density profiles are found to produce more gradual jet breaks while larger viewing angles cause smoother and later appearing jet breaks. The counter-jet becomes visible as it becomes sub-relativistic, and for k=0 this results in a clear bump-like feature in the light curve. However, for larger k values the jet decelerates more gradually, causing only a mild flattening in the radio light curve that might be hard to discern when k=2. Late time radio calorimetry is likely to consistently over-estimate the true energy by up to a factor of a few for k=2, and either over-predict or under-predict it by a smaller factor for k = 0,1.

preprint2011arXiv

The missing link: Merging neutron stars naturally produce jet-like structures and can power short Gamma-Ray Bursts

Short Gamma-Ray Bursts (SGRBs) are among the most luminous explosions in the universe, releasing in less than one second the energy emitted by our Galaxy over one year. Despite decades of observations, the nature of their "central-engine" remains unknown. Considering a binary of magnetized neutron stars and solving Einstein equations, we show that their merger results in a rapidly spinning black hole surrounded by a hot and highly magnetized torus. Lasting over 35 ms and much longer than previous simulations, our study reveals that magnetohydrodynamical instabilities amplify an initially turbulent magnetic field of ~ 10^{12} G to produce an ordered poloidal field of ~ 10^{15} G along the black-hole spin-axis, within a half-opening angle of ~ 30 deg, which may naturally launch a relativistic jet. The broad consistency of our ab-initio calculations with SGRB observations shows that the merger of magnetized neutron stars can provide the basic physical conditions for the central-engine of SGRBs.

preprint2010arXiv

Impulsive Acceleration of Strongly Magnetized Relativistic Flows

The strong variability of magnetic central engines of AGN and GRBs may result in highly intermittent strongly magnetized relativistic outflows. We find a new magnetic acceleration mechanism for such impulsive flows that can be much more effective than the acceleration of steady-state flows. This impulsive acceleration results in kinetic-energy-dominated flows at astrophysically relevant distances from the central source. For a spherical flow, a discrete shell ejected from the source over a time t_0 with Lorentz factor Gamma~1 and initial magnetization sigma_0 = B_0^2/(4 pi rho_0 c^2) >> 1 quickly reaches a typical Lorentz factor Gamma ~ sigma_0^{1/3} and magnetization sigma ~ sigma_0^{2/3} at the distance R_0 ~ ct_0. At this point the magnetized shell of width Delta ~ R_0 in the lab frame loses causal contact with the source and continues to accelerate by spreading significantly in its own rest frame. The expansion is driven by the magnetic pressure gradient and leads to relativistic relative velocities between the front and back of the shell. While the expansion is roughly symmetric in the center of momentum frame, in the lab frame most of the energy and momentum remain in a region (or shell) of width Delta ~ R_0 at the head of the flow. This acceleration proceeds as Gamma ~ (sigma_0 R/R_0)^{1/3} and sigma ~ sigma_0^{2/3} (R/R_0)^{-1/3} until reaching a coasting radius R_c ~ R_0 sigma_0^2 where the kinetic energy becomes dominant: Gamma ~ sigma_0 and sigma ~ 1 at R_c. Then the shell starts coasting and spreading (radially) causing its magnetization to drop as sigma ~ R_c/R at R>R_c. Given the typical variability time-scales of AGN and GRBs, the magnetic acceleration in these sources is a combination of the quasi-steady-state collimation acceleration close to the source and the impulsive acceleration further out. The interaction with the external medium is briefly discussed.

preprint2010arXiv

Jets and Gamma-Ray Burst Unification Schemes

There are several lines of evidence indicating that the ultra-relativistic outflows powering gamma-ray bursts (GRBs) are collimated into narrow jets. However, these are indirect, and the jet structure is rather poorly constrained. What is more, the jet dynamics have still not been investigated in detail. It has been suggested that the observed variety between different long duration events, ranging from bright spectrally hard GRBs, to dimmer and spectrally softer X-ray flashes (XRFs) may be largely due to different viewing angles (or lines of sight) relative to rather similar relativistic jets. Here we describe the current state of knowledge on these topics, explain some of the most relevant physics behind some of the basic principles, and discuss prospects for the future.

preprint2010arXiv

Late time observations of GRB080319B: jet break, host galaxy and accompanying supernova

The Swift discovered GRB080319B was by far the most distant source ever observed at naked eye brightness, reaching a peak apparent magnitude of 5.3 at a redshift of z=0.937. We present our late-time optical (HST, Gemini & VLT) and X-ray (Chandra) observations, which confirm that an achromatic break occurred in the power-law afterglow light curve at ~11 days post-burst. This most likely indicates that the gamma-ray burst (GRB) outflow was collimated, which for a uniform jet would imply a total energy in the jet E_{jet} \gsim 10^{52} erg. Our observations also show a late-time excess of red light, which is well explained if the GRB was accompanied by a supernova (SN), similar to those seen in some other long-duration GRBs. The latest observations are dominated by light from the host and show that the GRB took place in a faint dwarf galaxy (r(AB)\approx27.0, rest-frame M_B\approx-17.2). This galaxy is small even by the standards of other GRB hosts, which is suggestive of a low metallicity environment. Intriguingly, the properties of this extreme event - a small host and bright supernova - are entirely typical of the very low-luminosity bursts such as GRB980425 and GRB060218.

preprint2010arXiv

The long rapid decay phase of the extended emission from the short GRB 080503

GRB080503 was classified as a short GRB with extended emission (Perley et al. 2009). The origin of such extended emission (found in about a quarter of Swift short GRBs) is still unclear and may provide some clues to the identity of the elusive progenitors of short GRBs. The extended emission from GRB 080503 is followed by a rapid decay phase (RDP) that is detected over an unusually large dynamical range (one decade in time and ~3.5 decades in flux). We model the broad envelope of extended emission and the subsequent RDP using a physical model (Genet & Granot 2009), in which the prompt emission (and its tail) is the sum of its individual pulses (and their tails). For GRB 080503, a single pulse fit is found to be unacceptable. The RDP displays very strong spectral evolution and shows some evidence for the presence of two spectral components with different temporal behaviour, likely arising from distinct physical regions. A two pulse fit provides a much better fit to the data. The shallow gamma-ray and steep hard X-ray decays are hard to account for simultaneously, and require the second pulse to deviate from the simplest version of the model we use. Therefore, while high latitude emission is a viable explanation for the RDP in GRB080503, it is quite plausible that another mechanism is at work here. Finally, we note that the properties of the RDP following the extended emission of short GRBs appear to have different properties than that following the prompt emission of long GRBs. However, a larger sample of short GRBs with extended emission is required before any strong conclusion can be drawn.

preprint2009arXiv

Limits on Radioactive-Powered Emission Associated With a Short-Hard GRB 070724A in a Star-Forming Galaxy

We present results of an extensive observing campaign of the short duration, hard spectrum gamma-ray burst (GRB) 070724A, aimed at detecting the radioactively-powered emission that might follow from a binary merger or collapse involving compact objects. Our multi-band observations span the range in time over which this so-called Li-Paczynski mini-supernova could be active, beginning within 3 hours of the GRB trigger, and represent some of the deepest and most comprehensive searches for such emission. We find no evidence for such activity and place limits on the abundances and the lifetimes of the possible radioactive nuclides that could form in the rapid decompression of nuclear-density matter. Furthermore, our limits are significantly fainter than the peak magnitude of any previously detected broad-lined Type Ic supernova (SN) associated with other GRBs, effectively ruling out a long GRB-like SN for with this event. Given the unambiguous redshift of the host galaxy (z=0.456), GRB 070724A represents one of a small, but growing, number of short-hard GRBs for which firm physical/restframe quantities currently exist. The host of GRB 070724A is a moderately star-forming galaxy with an older stellar population component and a relatively high metallicity of 12+log(O/H)_KD02=9.1. We find no significant evidence for large amounts of extinction along the line of sight that could mask the presence of a SN explosion and estimate a small probability for chance alignment with the putative host. We discuss how our derived constraints fit into the evolving picture of short-hard GRBs, their potential progenitors, and the host environments in which they are thought to be produced.

preprint2007arXiv

Opacity Build-up in Impulsive Relativistic Sources

Opacity effects in relativistic sources of high-energy gamma-rays, such as gamma-ray bursts (GRBs) or Blazars, can probe the Lorentz factor of the outflow and the distance of the emission site from the source, and thus help constrain the composition of the outflow (protons, pairs, magnetic field) and the emission mechanism. Most previous works consider the opacity in steady state. Here we study time dependent effects of the opacity to pair production ($γγ\to e^+e^-$) in impulsive relativistic sources. We present a simple, yet rich, semi-analytic model for the time and energy dependence of the optical depth, $τ$, where an ultra-relativistic thin spherical shell emits isotropically in its own rest frame over a finite range of radii, $R_0 < R < R_0 + ΔR$. This is very relevant for GRB internal shocks. We find that in an impulsive source ($ΔR <~ R_0$), while the instantaneous spectrum has an exponential cutoff above the photon energy $ε_1(T)$ where $τ= 1$, the time integrated spectrum has a power-law high-energy tail above a photon energy $ε_{1*} \sim ε_1(ΔT)$ where $ΔT$ is the duration of the emission episode. Furthermore, photons with $ε> ε_{1*}$ will arrive mainly near the onset of the spike or flare from to the short emission episode, as in impulsive sources it takes time to build-up the (target) photon field, and thus $τ(ε)$ initially increases with time and $ε_1(T)$ decreases with time, so that photons of energy $ε> ε_{1*}$ are able to escape the source mainly very early on while $ε_1(T) > ε$. As the source approaches a quasi-steady state ($ΔR >> R_0$), the time integrated spectrum develops an exponential cutoff, while the power-law tail becomes increasingly suppressed.

preprint2006arXiv

The flat decay phase in the early X-ray afterglows of Swift GRBs

Many Swift GRBs show an early phase of shallow decay in their X-ray afterglows, lasting from $t \sim 10^{2.5} $s to $\sim 10^4 $s after the GRB, where the flux decays as $\sim t^{-0.2}-t^{-0.8}$. This is perhaps the most mysterious of the new features discovered by Swift in the early X-ray afterglow, since it is still not clear what causes it. I discuss different possible explanations for this surprising new discovery, as well as their potential implications for the gamma-ray efficiency, the afterglow kinetic energy, and perhaps even for the physics of collisionless relativistic shocks.

preprint2006arXiv

The Structure and Dynamics of GRB Jets

There are several lines of evidence which suggest that the relativistic outflows in gamma-ray bursts (GRBs) are collimated into narrow jets. The jet structure has important implications for the true energy release and the event rate of GRBs, and can constrain the mechanism responsible for the acceleration and collimation of the jet. Nevertheless, the jet structure and its dynamics as it sweeps up the external medium and decelerates, are not well understood. In this review I discuss our current understanding of GRB jets, stressing their structure and dynamics.

preprint2004arXiv

Implications of the Measured Image Size for the Radio Afterglow of GRB 030329

We use data on the image size of the radio afterglow of GRB 030329 (Taylor et al. 2004) to constrain the physical parameters of this explosion. Together with the observed broad band spectrum, this data over-constrains the physical parameters, thus enabling to test different GRB jet models for consistency. We consider two extreme models for the lateral spreading of the jet: model 1 with relativistic expansion in the local rest frame, and model 2 with little lateral expansion as long as the jet is highly relativistic. We find that both models are consistent with the data for a uniform external medium, while for a stellar wind environment model 1 is consistent with the data but model 2 is disfavored by the data. Our derivations can be used to place tighter constraints on the dynamics and structure of GRB jets in future afterglows, following a denser monitoring campaign for the temporal evolution of their image size.

preprint2002arXiv

Constraining the Structure of GRB Jets Through the Afterglow Light Curves

We investigate the effect that the structure of GRB jets has on the afterglow light curves for observers located at different viewing angles, $θ_{obs}$, from the jet symmetry axis. The largest uncertainty in the jet dynamics is the degree of lateral energy transfer. Thus, we use two simple models, that make opposite and extreme assumptions for this point, and calculate the light curves for an external density that is either homogeneous, or decreases as the square of the distance from the source. The Lorentz factor and kinetic energy per unit solid angle are initially taken to be power laws of the angle from the jet axis. We perform a qualitative comparison between the resulting light curves and afterglow observations. This constrains the jet structure, and poses problems for a `universal' jet model, where all GRB jets are assumed to be intrinsically identical, and differ only by our viewing angle, $θ_{obs}$.

preprint2002arXiv

The Mean Number of Extra Micro-Image Pairs for Macro-Lensed Quasars

When a gravitationally lensed source crosses a caustic, a pair of images is created or destroyed. We calculate the mean number of such pairs of micro-images $<n>$ for a given macro-image of a gravitationally lensed point source, due to microlensing by the stars of the lensing galaxy. This quantity was calculated by Wambsganss, Witt & Schneider (1992) for the case of zero external shear, $γ=0$, at the location of the macro-image. Since in realistic lens models a non-zero shear is expected to be induced by the lensing galaxy, we extend this calculation to a general value of $γ$. We find a complex behavior of $<n>$ as a function of $γ$ and the normalized surface mass density in stars $κ_*$. Specifically, we find that at high magnifications, where the average total magnification of the macro-image is $<μ>=|(1-κ_*)^2-γ^2|^{-1}\gg 1$, $<n>$ becomes correspondingly large, and is proportional to $<μ>$. The ratio $<n>/<μ>$ is largest near the line $γ=1-κ_*$ where the magnification $<μ>$ becomes infinite, and its maximal value is 0.306. We compare our semi-analytic results for $<n>$ to the results of numerical simulations and find good agreement. We find that the probability distribution for the number of extra image pairs is reasonably described by a Poisson distribution with a mean value of $<n>$, and that the width of the macro-image magnification distribution tends to be largest for $<n>\sim 1$.

preprint2001arXiv

The Shape of Spectral Breaks in GRB Afterglows

Gamma-Ray Burst (GRB) afterglows are well described by synchrotron emission from relativistic blast waves expanding into an external medium. The blast wave is believed to amplify the magnetic field and accelerate the electrons into a power law distribution of energies promptly behind the shock. These electrons then cool both adiabatically and by emitting synchrotron and inverse Compton radiation. The resulting spectra is known to consist several power law segments, which smoothly join at certain break frequencies. Here, we give a complete description of all possible spectra under those assumptions, and find that there are 5 possible regimes, depending on the ordering of the break frequencies. The flux density is calculated by integrating over the contributions from all the shocked region, using the Blandford McKee solution. This allows us to calculate more accurate expressions for the value of these break frequencies, and describe the shape of the spectral breaks around them. This also provides the shape of breaks in the light curves caused by the passage of a break frequency through the observed band. These new, more exact, estimates are different from more simple calculations by up to a factor of about 70, and describe some new regimes which where previously ignored.

preprint1998arXiv

Images and Spectra From the Interior of a Relativistic Fireball

The detection of an afterglow, following a gamma-ray burst (GRB), can be described reasonably well by synchrotron emission from a relativistic spherical expanding blast wave, driven by an expanding fireball. We perform detailed calculations considering the emission from the whole region behind the shock front. We use the Blandford McKee self similar adiabatic solution to describe the fluid behind the shock. Using this detailed model, we derive expressions for the peak flux, and the peak frequency at a given observed time. These expressions provide important numerical corrections to previous, more simplified models. We calculate the observed light curve and spectra for several magnetic field models. We show that both the light curve and the spectra are flat near the peak. This rules out the interpretation of the optical peak of GRB970508 as the peak of the light curve, predicted by the existing fireball models. We calculate the observed image of a GRB afterglow. The observed image is bright near the edge and dimmer at the center, thus creating a ring. The contrast between the edge and the center is larger at high frequencies and the width of the ring is smaller.