Rotational effects in thermonuclear Type I Bursts: equatorial crossing and directionality of flame spreading

In a previous study on thermonuclear (type I) nursts on accreting neutron stars we addressed and demonstrated the importance of the effects of rotation, through the Coriolis force, on the propagation of the burning flame. However, that study only analysed cases of longitudinal propagation, where the Coriolis force coefficient $2Omegacostheta$ was constant. In this paper, we study the effects of rotation on propagation in the meridional (latitudinal) direction, where the Coriolis force changes from its maximum at the poles to zero at the equator. We find that the zero Coriolis force at the equator, while affecting the structure of the flame, does not prevent its propagation from one hemisphere to another. We also observe structural differences between the flame propagating towards the equator and that propagating towards the pole, the second being faster. In the light of the recent discovery of the low spin frequency of burster IGR~J17480-2446 rotating at 11 Hz (for which Coriolis effects should be negligible) we also extend our simulations to slow rotation.

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.

Flame Propagation on the Surfaces of Rapidly Rotating Neutron Stars during Type I X-ray Bursts

We present the first vertically resolved hydrodynamic simulations of a laterally propagating, deflagrating flame in the thin helium ocean of a rotating accreting neutron star. We use a new hydrodynamics solver tailored to deal with the large discrepancy in horizontal and vertical length scales typical of neutron star oceans, and which filters out sound waves that would otherwise limit our timesteps. We find that the flame moves horizontally with velocities of order $10^5$ cm s$^{-1}$, crossing the ocean in few seconds, broadly consistent with the rise times of Type I X-ray bursts. We address the open question of what drives flame propagation, and find that heat is transported from burning to unburnt fuel by a combination of top-to-bottom conduction and mixing driven by a baroclinic instability. The speed of the flame propagation is therefore a sensitive function of the ocean conductivity and spin: we explore this dependence for an astrophysically relevant range of parameters and find that in general flame propagation is faster for slower rotation and higher conductivity.

Quasi-Periodic Oscillations and broadband variability in short magnetar bursts

The discovery of quasi-periodic oscillations (QPOs) in magnetar giant flares has opened up prospects for neutron star asteroseismology. However, with only three giant flares ever recorded, and only two with data of sufficient quality to search for QPOs, such analysis is seriously data limited. We set out a procedure for doing QPO searches in the far more numerous, short, less energetic magnetar bursts. The short, transient nature of these bursts requires the implementation of sophisticated statistical techniques to make reliable inferences. Using Bayesian statistics, we model the periodogram as a combination of red noise at low frequencies and white noise at high frequencies, which we show is a conservative approach to the problem. We use empirical models to make inferences about the potential signature of periodic and quasi-periodic oscillations at these frequencies. We compare our method with previously used techniques and find that although it is on the whole more conservative, it is also more reliable in ruling out false positives. We illustrate our Bayesian method by applying it to a sample of 27 bursts from the magnetar SGR J0501+4516 observed by the Fermi Gamma-ray Burst Monitor, and we find no evidence for the presence of QPOs in any of the bursts in the unbinned spectra, but do find a candidate detection in the binned spectra of one burst. However, whether this signal is due to a genuine quasi-periodic process, or can be attributed to unmodeled effects in the noise is at this point a matter of interpretation.

Observable QPOs produced by steep pulse profiles in Magnetar Flares

Strong quasi-periodic oscillations in the tails of the giant gamma-ray flares seen in SGR 1806-20 and SGR 1900+14 are thought to be produced by starquakes in the flaring magnetar. However, the large fractional amplitudes (up to ~20%) observed are difficult to reconcile with predicted amplitudes of starquakes. Here we demonstrate that the steeply pulsed emission profile in the tail of the giant flare can enhance the observed amplitude of the underlying oscillation, analogously to a beam of light oscillating in and out of the line of sight. This mechanism will also broaden the feature in the power spectrum and introduce power at harmonics of the oscillation. The observed strength of the oscillation depends on the amplitude of the underlying starquake, the orientation and location of the emission on the surface of the star, and the gradient of the light curve profile. While the amplification of the signal can be significant, we demonstrate that even with uncertainties in the emission geometry, this effect is not sufficient to produce the observed QPOs. This result excludes the direct observation of a starquake, and suggests that the observed variations come from modulations in the intensity of the emission.

Photospheric radius expansion during magnetar bursts

On August 24th 2008 the new magnetar SGR 0501+4516 (discovered by SWIFT) emitted a bright burst with a pronounced double-peak structure in hard X-rays, reminiscent of the double-peak temporal structure seen in some bright thermonuclear bursts on accreting neutron stars. In the latter case this is due to Photospheric Radius Expansion (PRE): when the flux reaches the Eddington limit, the photosphere expands and cools so that emission becomes softer and drops temporarily out of the X-ray band, re-appearing as the photosphere settles back down. We consider the factors necessary to generate double-peaked PRE events, and show that such a mechanism could plausibly operate in magnetar bursts, despite the vastly different emission process. Identification of the magnetic Eddington limit in a magnetar would constrain magnetic field and distance and could, in principle, enable a measurement of gravitational redshift. It would also locate the emitting region at the neutron star surface, constraining the burst trigger mechanism. Conclusive confirmation of PRE events will require more detailed radiative models for bursts. However for SGR 0501+4516 the predicted critical flux (using the magnetic field strength inferred from timing and the distance suggested by its probable location in the Perseus arm of our Galaxy) is consistent with that observed in the August 24th burst.

The Luminosity and Energy Dependence of Pulse Phase Lags in the Accretion-powered Millisecond Pulsar SAX J1808.4-3658

Soft phase lags, in which X-ray pulses in lower energy bands arrive later than pulses in higher energy bands, have been observed in nearly all accretion-powered millisecond pulsars, but their origin remains an open question. In a study of the 2.5 ms accretion-powered pulsar SAX J1808.4-3658, we report that the magnitude of these lags is strongly dependent on the accretion rate. During the brightest stage of the outbursts from this source, the lags increase in magnitude as the accretion rate drops; when the outbursts enter their dimmer flaring-tail stage, the relationship reverses. We evaluate this complex dependence in the context of two theoretical models for the lags, one relying on the scattering of photons by the accretion disk and the other invoking a two-component model for the photon emission. In both cases, the turnover suggests that we are observing the source transitioning into the propeller accretion regime.

The Swift capture of a long X-ray burst from XTE J1701-407

XTE J1701-407 is a new transient X-ray source discovered on June 8th, 2008. More than one month later it showed a rare type of thermonuclear explosion: a long type I X-ray burst. We report herein the results of our study of the spectral and flux evolution during this burst, as well as the analysis of the outburst in which it took place. We find an upper limit on the distance to the source of 6.1 kpc by considering the maximum luminosity reached by the burst. We measure a total fluence of 3.5*10^{-6} erg/cm^2 throughout the ~20 minutes burst duration and a fluence of 2.6*10^{-3} erg/cm^2 during the first two months of the outburst. We show that the flux decay is best fitted by a power law (index ~1.6) along the tail of the burst. Finally, we discuss the implications of the long burst properties, and the presence of a second and shorter burst detected by Swift ten days later, for the composition of the accreted material and the heating of the burning layer.