Magnetar giant flares show oscillatory modulations in the tails of their light curves, which can only be explained via some form of beaming. The fireball model for magnetar bursts has been used successfully to fit the phase-averaged light curves of the tails of giant flares, but so far no attempts have been made to fit the pulsations. We present a relatively simple numerical model to simulate beaming of magnetar flare emission. In our simulations, radiation escapes from the base of a fireball trapped in a dipolar magnetic field, and is scattered through the optically thick magnetosphere of the magnetar until it escapes. Beaming is provided by the presence of a relativistic outflow, as well as by the geometry of the system. We find that a simple picture for the relativistic outflow is enough to create the pulse fraction and sharp peaks observed in pulse profiles of magnetar flares, while without a relativistic outflow the beaming is insufficient to explain giant flare rotational modulations.

We report on time-resolved spectroscopy of the 63 brightest bursts of SGR J1550-5418, detected with Fermi/Gamma-ray Burst Monitor during its 2008-2009 intense bursting episode. We performed spectral analysis down to 4 ms time-scales, to characterize the spectral evolution of the bursts. Using a Comptonized model, we find that the peak energy, E_peak, anti-correlates with flux, while the low-energy photon index remains constant at -0.8 up to a flux limit F~10^-5 erg s-1 cm-2. Above this flux value the E_peak-flux correlation changes sign, and the index positively correlates with flux reaching 1 at the highest fluxes. Using a two black-body model, we find that the areas and fluxes of the two emitting regions correlate positively. Further, we study here for the first time, the evolution of the temperatures and areas as a function of flux. We find that the area-kT relation follows lines of constant luminosity at the lowest fluxes, R^2 \propto kT^-4, with a break at higher fluxes ($F>10^-5.5 erg s-1 cm-2). The area of the high-kT component increases with flux while its temperature decreases, which we interpret as due to an adiabatic cooling process. The area of the low-kT component, on the other hand, appears to saturate at the highest fluxes, towards R_max~30 km. Assuming that crust quakes are responsible for SGR bursts and considering R_max as the maximum radius of the emitting photon-pair plasma fireball, we relate this saturation radius to a minimum excitation radius of the magnetosphere, and put a lower limit on the internal magnetic field of SGR J1550-5418, B_int>~4.5x10^15 G.

We investigate the possibility of Photospheric Radius Expansion (PRE) during magnetar bursts. Identification of PRE would enable a determination of the magnetic Eddington limit (which depends on field strength and neutron star mass and radius), and shed light on the burst mechanism. To do this we model hydrostatic atmospheres in a strong radial magnetic field, determining both their maximum extent and photospheric temperatures. We find that spatially-extended atmospheres cannot exist in such a field configuration: typical maximum extent for magnetar-strength fields is ~10 m (as compared to 200 km in the non-magnetic case). Achieving balance of gravitational and radiative forces over a large range of radii, which is critical to the existence of extended atmospheres, is rendered impossible in strong fields due to the dependence of opacities on temperature and field strength. We conclude that high luminosity bursts in magnetars do not lead to expansion and cooling of the photosphere, as in the non-magnetic case. We also find the maximum luminosity that can propagate through a hydrostatic magnetar atmosphere to be lower than previous estimates. The proximity and small extent of the photospheres associated with the two different polarization modes also calls into question the interpretation of two blackbody fits to magnetar burst spectra as being due to extended photospheres.