An Energetic Blast Wave from the 2004 December 27 Giant Flare of the Soft Gamma-Ray Repeater SGR 1806-20

This article-based dissertation provides a review on the broad subject of magnetars-their characteristics, giant flares (GFs) and associated observations of X-ray, gamma-ray, and radio emissions and their proposed physical mechanisms. The primary purpose of this dissertation is to provide an extensive description of the two research projects I undertook during my tenure as a Graduate Research Assistant, under the guidance of my advisor. Broadly, my research was focused on building analytical models and running three-dimensional (3-D), high-resolution magnetohydrodynamic (MHD) simulations using the astrophysical PLUTO code to investigate the physical mechanisms behind high-energy (X-ray and gamma-ray) and radio emissions associated with magnetar GFs using observational constraints. This, in turn, aided in either validating or disfavoring existing theories behind such energetic explosions.Chapter 1 provides a review on magnetars, their GFs and associated high-energy and radio emissions, largely based on excellent reviews by [1]–[5]. I summarize interesting observational features of magnetars, specifically those of soft gamma-ray repeaters (SGRs) and anomalous X-ray pulsars (AXPs), along with known aspects of their X-ray and gamma-ray activity. I focus on the December 27, 2004 GF emitted by SGR 1806-20, the most energetic GF out of the three that occurred to date, describe its energetics and summarize existing theories behind the physical mechanisms that give rise to two emission characteristics associated with the GF - (i) quasi-periodic oscillations (QPOs) seen in the tail, and (ii) a radio afterglow detected a week after the GF. Lastly, I describe the methods I used to hypothesize the physical mechanisms behind QPOs and the radio emission and compare and contrast them with those suggested previously.In chapter 2, I present a version of the research article in preparation and pending publication in the Monthly Notices of the Royal Astronomical Society. The work titled “Radio afterglow of magnetars’ giant flares”, undertaken under the supervision of Dr. Maxim Lyutikov and in collaboration with Dr. Maxim Barkov, explores the possible physical mechanisms behind the radio afterglow associated with the SGR 1806-20 GF using high-resolution 3-D MHD simulations.In chapter 3, I present a version of the research article previously published by the Journal of Plasma Physics. The work titled “Tilting instability of magnetically confined spheromaks”, undertaken under the supervision of Dr. Maxim Lyutikov, in collaboration with Dr. Lorenzo Sironi and Dr. Maxim Barkov, investigates the tilting instability of a magnetically confined spheromak using 3-D MHD and relativistic particle-in-cell (PIC) simulations with an application to astrophysical plasmas, specifically to explain the QPOs arising in the tail of the SGR 1806-20 GF.I summarize the main results and conclusions of the two research projects and describe future prospects in chapter 4, followed by appendices A and B which describe additional theoretical concepts and simulation results for a better understanding of the nature of radio afterglows associated with GFs, and structure of spheromaks. References are compiled after the appendices in order that they are first cited, followed by a brief autobiographical sketch, and a list of publications.

The knowledge of the rate of soft gamma-ray repeater (SGR) giant flares is important for understanding the giant flare mechanism and the SGR energy budget in the framework of the magnetar model. We estimate the upper limit to the rate using the results of an extensive search for extragalactic soft gamma-repeater giant flares (GFs) among 140 short gamma-ray bursts detected between 1994 and 2010 by Konus-Wind using InterPlanetary Network (IPN) localizations and temporal parameters. We show that Konus-Wind and the IPN are capable of detecting GFs with energies of 2.3x10^46 erg (which is the energy of the GF from SGR 1806-20 assuming a distance of 15 kpc) at distances of up to about 30 Mpc and GFs with energies of <10^45 erg (which is the energy of the GF from SGR 0526-66) at distances of up to about 6 Mpc. Using a sample of 1896 nearby galaxies we found that only two bursts, GRB 051103 and GRB 070201, have a low chance coincidence probability between an IPN localization and a nearby galaxy. We found the upper limit to the fraction of GFs among short GRBs with fluence above ~5x10^-7 erg cm^-2 to be <8% (95% confidence level). Assuming that the number of active SGRs in nearby galaxies is proportional to their core-collapse supernova rate, we derived the one-sided 95% upper limit to the rate of GFs with energy output similar to the GF from SGR 1806-20 to be (0.6--1.2)x10^-4 Q_46^-1.5 yr^-1 per SGR, where Q_46 is the GF energy output in 10^46 erg.

Soft gamma-ray repeaters (SGRs) are a peculiar family of bursting neutron stars that, occasionally, have been observed to emit extremely energetic giant flares (GFs), with energy release up to approximately 10(47) ergs(-1). These are exceptional and rare events. It has been recently proposed that GFs, if emitted by extragalactic SGRs, may appear at Earth as short gamma-ray bursts. Here, I will discuss the properties of the GFs observed in SGRs, with particular emphasis on the spectacular event registered from SGR 1806-20 in December 2004. I will review the current scenario for the production of the flare, within the magnetar model, and the observational implications.

It was suggested that some of the short-duration gamma-ray bursts (GRBs) are giant flares of soft gamma-ray repeaters (SGRs) in nearby galaxies. To test this hypothesis, I have constructed a sample of 47 short GRBs, detected by the Interplanetary Network (IPN), for which the position is constrained by at least one annulus on the celestial sphere. For each burst, I have checked whether its IPN 3 σ error region coincides with the apparent disk of one of 316 bright, star-forming galaxies found within 20 Mpc. I find a single match of GRB 000420B with M74, which could, however, be due to a chance coincidence. I estimate the IPN efficiency as a function of fluence and derive the galaxy sample completeness. I find that assuming there is a cutoff in the observed energy distribution of SGR flares at ≤10^(47) ergs, the fraction of SGRs among short GRBs with fluence above ~10^(-5) ergs cm^(-2) is <16% (95% confidence). I estimate the number of active SGRs in each one of the galaxies in the sample, and combine it with the distances to these galaxies, the IPN efficiency, and the SGR flare energy distribution, to derive the rate of giant flares with energy above 4 × 10^(46) ergs. I find that the rate of such giant flares is about (0.4-5) × 10^(-4) yr^(-1) per SGR. This rate is marginally consistent with the observed Galactic rate. Comparison of the Galactic rate with the inferred extragalactic rate implies a gradual cutoff (or steepening) of the flare energy distribution at ≾3 × 10^(46) ergs (95% confidence). Using the Galactic SGR flare rate, I set a lower limit of 1% on the fraction of SGR flares among short GRBs.

Soft gamma-ray repeaters (SGRs) are 'magnetars', a small class of slowly spinning neutron stars with extreme surface magnetic fields, B approximately 10(15) gauss (refs 1 , 2 -3). On 27 December 2004, a giant flare was detected from the magnetar SGR 1806-20 (ref. 2), only the third such event recorded. This burst of energy was detected by a variety of instruments and even caused an ionospheric disturbance in the Earth's upper atmosphere that was recorded around the globe. Here we report the detection of a fading radio afterglow produced by this outburst, with a luminosity 500 times larger than the only other detection of a similar source. From day 6 to day 19 after the flare from SGR 1806-20, a resolved, linearly polarized, radio nebula was seen, expanding at approximately a quarter of the speed of light. To create this nebula, at least 4 x 10(43) ergs of energy must have been emitted by the giant flare in the form of magnetic fields and relativistic particles.

The giant flare observed on 2004 December 27 from SGR 1806−20 has revived the idea that a fraction of short (&amp;lt;2 s) gamma-ray bursts (GRBs) are due to giant flares from soft gamma-ray repeaters (SGRs) located in nearby galaxies. One of the distinguishing characteristics of these events is the thermal (blackbody) spectrum with temperatures ranging from ∼50 to ∼180 keV, with the highest temperature observed for the initial 0.2-s spike of the 2004 December 27 event. We have analysed the spectra of a complete sample of short GRBs with peak fluxes greater than 4 photon s−1 cm−2 detected by BATSE. Of the 115 short GRBs so selected, only 76 had sufficient signal-to-noise ratio to allow the spectral analysis. We find only three short GRBs with a spectrum well fitted by a blackbody, with 60 ≲kT≲ 90 keV, albeit with a considerably longer duration (i.e. ≳1 s) and a more complex light curve than the 2004 December 27 event. This implies a stringent limit on the rate of extragalactic SGR giant flares with spectral properties analogous to the December 27 flare. We conclude that up to 4 per cent of the short GRBs could be associated with giant flares (2σ confidence). This implies that either the distance to SGR 1806−20 is smaller than 15 kpc or the rate of Galactic giant flares is lower than the estimated 0.033 yr−1.

On 27 August 1998, the Soft Gamma Ray Repeater (SGR) 1900+14, which is an exotic neutron star located near the Galactic Center, produced a giant flare. Gamma rays from the giant flare unusually ionized the lower ionosphere, and the ionospheric disturbance was detected as a large‐amplitude change of the VLF signal whose propagation distance is relatively short (870 km). The peak flux of the flare was so huge that it saturated all the gamma ray detectors on the spacecraft; consequently, the flux and the spectrum during the most intense period was poorly determined. Tanaka et al. (2007) have recently derived the accurate peak flux of the flare from the Geotail data. By means of model calculations based on this accurate estimation and their comparisons with the short‐distance VLF data, we have found that the spectrum during the most intense period was one temperature (kT = 240 keV) optically thin thermal Bremsstrahlung (OTTB). This result provides us with a clue to reveal the emission and triggering mechanisms of the giant flare.

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 Gamma(0) similar to 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 Gamma(0) < 50 if the electron distribution is a single power law. For a broken power-law distribution of electrons, a higher Gamma(0) is allowed. At energies higher than I 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 Gamma(0) (e.g. < 10). The forward shock also accelerates cosmic rays to the maximum energy 10(17) eV, and generates neutrinos with a typical energy 10(14) eV through photomeson interaction with the X-ray tail photons. However, they are too weak to be detectable.

Anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs) are magnetar candidates, i.e., neutron stars powered by a strong magnetic field. If they are indeed magnetars, they will emit high-energy gamma rays that are detectable by the Fermi Large Area Telescope (LAT), according to the outer gap model. However, no significant detection is reported in recent Fermi-LAT observations of all known AXPs and SGRs. Considering the discrepancy between theory and observations, we calculate the theoretical spectra for all AXPs and SGRs with sufficient observational parameters. Our results show that most AXPs and SGRs are high-energy gamma-ray emitters if they are really magnetars. The four AXPs 1E 1547.0−5408, XTE J1810−197, 1E 1048.1−5937, and 4U 0142+61 should have been detected by Fermi-LAT. There is therefore a conflict between the outer gap model in the case of magnetars and Fermi observations. Possible explanations in the magnetar model are discussed. On the other hand, if AXPs and SGRs are fallback disk systems, i.e., accretion-powered for the persistent emissions, most of them are not high-energy gamma-ray emitters. Future deep Fermi-LAT observations of AXPs and SGRs will help us make clear whether they are magnetars or fallback disk systems.

We study magneto-elastic oscillations of highly magnetized neutron stars (magnetars) which have been proposed as an explanation for the quasi-periodic oscillations (QPOs) appearing in the decaying tail of the giant flares of soft gamma-ray repeaters (SGRs). We extend previous studies by investigating various magnetic field configurations, computing the Alfv\'en spectrum in each case and performing magneto-elastic simulations for a selected number of models. By identifying the observed frequencies of 28 Hz (SGR 1900+14) and 30 Hz (SGR 1806-20) with the fundamental Alfv\'en QPOs, we estimate the required surface magnetic field strength. For the magnetic field configurations investigated (dipole-like poloidal, mixed toroidal-poloidal with a dipole-like poloidal component and a toroidal field confined to the region of field lines closing inside the star, and for poloidal fields with an additional quadrupole-like component) the estimated dipole spin-down magnetic fields are between 8x10^14 G and 4x10^15 G, in broad agreement with spin-down estimates for the SGR sources producing giant flares. A number of these models exhibit a rich Alfv\'en continuum revealing new turning points which can produce QPOs. This allows one to explain most of the observed QPO frequencies as associated with magneto-elastic QPOs. In particular, we construct a possible configuration with two turning points in the spectrum which can explain all observed QPOs of SGR 1900+14. Finally, we find that magnetic field configurations which are entirely confined in the crust (if the core is assumed to be a type I superconductor) are not favoured, due to difficulties in explaining the lowest observed QPO frequencies (f<30 Hz).

We present a LIGO search for short-duration gravitational waves (GWs) associated with soft gamma ray repeater (SGR) bursts. This is the first search sensitive to neutron star f modes, usually considered the most efficient GW emitting modes. We find no evidence of GWs associated with any SGR burst in a sample consisting of the 27 Dec. 2004 giant flare from SGR 1806-20 and 190 lesser events from SGR 1806-20 and SGR 1900+14. The unprecedented sensitivity of the detectors allows us to set the most stringent limits on transient GW amplitudes published to date. We find upper limit estimates on the model-dependent isotropic GW emission energies (at a nominal distance of 10 kpc) between 3x10;{45} and 9x10;{52} erg depending on waveform type, detector antenna factors and noise characteristics at the time of the burst. These upper limits are within the theoretically predicted range of some SGR models.

Giant flares from soft gamma-ray repeaters (SGRs) are one of the most violent phenomena in neutron stars. Quasi-periodic oscillations (QPOs) with frequencies ranging from 18 to 1840 Hz have been discovered in the tails of giant flares from two SGRs, and were ascribed to be seismic vibrations or torsional oscillations of magnetars. Here we propose an alternative explanation for the QPOs in terms of standing sausage mode oscillations of flux tubes in the magnetar coronae. We show that most of the QPOs observed in SGR giant flares could be well accounted for except for those with very high frequencies (625 and 1840 Hz).

Radio afterglows have been detected following two giant flares from soft gamma repeaters (i.e. SGR1900+14, SGR 1806‐20). Recent follow‐up observations of the December 27 giant flare of SGR 1806‐20 have detected a multi‐frequency radio afterglow from 240 MHz to 8.46 GHz, extending in time from one week to about one month after the flare. The angular size of the source was also measured for the first time. Here we show that this radio afterglow gives the first piece of clear evidence that an energetic blast wave sweeps up its surrounding medium and produces a synchrotron afterglow, the same mechanism as established for GRB afterglows.

We examine two trigger mechanisms, one internal and the other external to the neutron star, that give rise to the intense soft gamma-ray repeater (SGR) giant flares. So far, three giant flares have been observed from the three out of the seven confirmed SGRs on March 5, 1979, August 27, 1998, and December 27, 2004. The last two events were found to be much more powerful than the first, and both showcased the existence of a precursor, that we show to have had initiated the main flare. In the internal mechanism, we propose that the strongly wound up poloidal magnetic field develops tangential discontinuities and dissipates its torsional energy in heating the crust. The timescale for the instability to develop coincides with the duration of the quiescent state that followed the precursor. Alternatively, we develop a reconnection model based on the hypothesis that shearing motion of the footpoints causes the materialization of a Sweet-Parker current layer in the magnetosphere. The thinning of this macroscopic layer due to the development of an embedded super-hot turbulent current layer switches on the impulsive Hall reconnection, which powers the giant flare. Again, we show that the thinning time is on the order of the preflare quiescent time. This model naturally explains the origin of the observed nonthermal radiation during the flares, as well as the post flare radio afterglows.

In January 2009, multiple short bursts of soft gamma-rays were detected from the direction of the anomalous X-ray pulsar 1E 1547.0-5408 by different satellites. Here we report on the observations obtained with the INTEGRAL SPI-ACS detector during the period with the strongest bursting activity. More than 200 bursts were detected at energies above 80 keV in a few hours on January 22. Among these, two remarkably bright events showed pulsating tails lasting several seconds and modulated at the 2.1 s spin period of 1E 1547.0-5408. The energy released in the brightest of these bursts was of a few 10^43 erg, for an assumed distance of 10 kpc. This is smaller than that of the three giant flares seen from soft gamma-ray repeaters, but higher than that of typical bursts from soft gamma-ray repeaters and anomalous X-ray pulsars.