We present a summary of the long-term evolution of various properties of the five non-transient Anomalous X-ray Pulsars (AXPs) 1E 1841-045, RXS J170849.0-400910, 1E 2259+586, 4U 0142+61, and 1E 1048.1-5937, regularly monitored with RXTE from 1996 to 2012. We focus on three properties of these sources: the evolution of the timing, pulsed flux, and pulse profile. We report several new timing anomalies and radiative events, including a putative anti-glitch seen in 1E 2259+586 in 2009, and a second epoch of very large spin-down rate fluctuations in 1E 1048.1-5937 following a large flux outburst. We compile the properties of the 11 glitches and 4 glitch candidates observed from these 5 AXPs between 1996 and 2012. Overall, these monitoring observations reveal several apparent patterns in the behavior of this sample of AXPs: large radiative changes in AXPs (including long-lived flux enhancements, short bursts, and pulse profile changes) are rare, occurring typically only every few years per source; large radiative changes are almost always accompanied by some form of timing anomaly, usually a spin-up glitch; only 20-30% of timing anomalies are accompanied by any form of radiative change. We find that AXP radiative behavior at the times of radiatively loud glitches is sufficiently similar to suggest common physical origins. The similarity in glitch properties when comparing radiatively loud and radiatively silent glitches in AXPs suggests a common physical origin in the stellar interior. Finally, the overall similarity of AXP and radio pulsar glitches suggests a common physical origin for both phenomena.

We present the results of Rossi X-ray Timing Explorer (RXTE) and Swift monitoring observations of the magnetar 1E 1547.0-5408 following the pulsar's radiative outbursts in 2008 October and 2009 January. We report on a study of the evolution of the timing properties and the pulsed flux from 2008 October 4 through 2009 December 26. We show that the pulsed flux decrease which followed an initial rise in the 2008 outburst was interrupted by a spike ~9 days after the initial outburst. In our timing study, a phase-coherent analysis shows that for the first 29 days following the 2008 outburst, there was a very fast increase in the magnitude of the rotational frequency derivative nudot, such that the second derivative was a factor of ~60 larger than that reported in data from 2007. This nudot magnitude increase occurred in concert with the decay of the pulsed flux following the start of the 2008 event. Following the 2009 outburst, for the first 23 days, the second derivative was consistent with zero, and nudot had returned to close to its 2007 value. In contrast to the 2008 event, the 2009 outburst showed a major increase in persistent flux, relatively little change in the pulsed flux, and sudden significant spectral hardening ~15 days after the outburst. We show that, excluding the month following each of the outbursts, and because of the noise and the sparsity in the data, multiple plausible timing solutions fit the pulsar's frequency behavior. We note similarities in the behavior of 1E 1547.0-5408 following the 2008 outburst to that seen in the AXP 1E 1048.1-5937 following its 2001-2002 outburst and discuss this in terms of the magnetar model.

We present an analysis of the extended emission around the anomalous X-ray pulsar 1E 1547.0-5408 using four XMM-Newton observations taken with the source in varying states of outburst as well as in quiescence. We find that the extended emission flux is highly variable and strongly correlated with the flux of the magnetar. Based on this result, as well as on spectral and energetic considerations, we conclude that the extended emission is dominated by a dust-scattering halo and not a pulsar wind nebula (PWN), as has been previously argued. We obtain an upper limit on the 2-10 keV flux of a possible PWN of 4.7e-14 erg/s/cm^2, three times less than the previously claimed value, implying an efficiency for conversion of spin-down energy into nebular luminosity of <9e-4 (assuming a distance of 4 kpc). We do, however, find strong evidence for X-ray emission from the supernova remnant shell surrounding the pulsar, as previously reported.

We report on radio observations of five magnetars and two magnetar candidates carried out at 1950 MHz with the Green Bank Telescope in 2006-2007. The data from these observations were searched for periodic emission and bright single pulses. Also, monitoring observations of magnetar 4U0142+61 following its 2006 X-ray bursts were obtained. No radio emission was detected was detected for any of our targets. The non-detections allow us to place luminosity upper limits (at 1950 MHz) of approximately L < 1.60 mJy kpc^2 for periodic emission and L < 7.6 Jy kpc^2 for single pulse emission. These are the most stringent limits yet for the magnetars observed. The resulting luminosity upper limits together with previous results are discussed, as is the importance of further radio observations of radio-loud and radio-quiet magnetars.

After at least 6 years of quiescence, Anomalous X-ray Pulsar (AXP) 4U 0142+61 entered an active phase in 2006 March that lasted several months and included six X-ray bursts as well as many changes in the persistent X-ray emission. The bursts, the first seen from this AXP in >11 years of Rossi X-ray Timing Explorer monitoring, all occurred in the interval between 2006 April 6 and 2007 February 7. The burst durations ranged from 0.4-1.8\times103 s. The first five burst spectra are well modeled by blackbodies, with temperatures kT ~ 2-9 keV. However, the sixth burst had a complicated spectrum that is well characterized by a blackbody plus two emission features whose amplitude varied throughout the burst. The most prominent feature was at 14.0 keV. Upon entry into the active phase the pulsar showed a significant change in pulse morphology and a likely timing glitch. The glitch had a total frequency jump of (1.9\pm0.4)\times10-7 Hz, which recovered with a decay time of 17\pm 2 days by more than the initial jump, implying a net spin-down of the pulsar. Within the framework of the magnetar model, the net spin-down of the star could be explained by regions of the superfluid that rotate slower than the rest. The bursts, flux enhancements, and pulse morphology changes can be explained as arising from crustal deformations due to stresses imposed by the highly twisted internal magnetic field. However, unlike other AXP outbursts, we cannot account for a major twist being implanted in the magnetosphere.

We present results from observations of the magnetar 1E 1547.0-5408 (SGR J1550-5418) taken with the Chandra X-ray Observatory and the Rossi X-ray Timing Explorer (RXTE) following the source's outbursts in 2008 October and 2009 January. During the time span of the Chandra observations, which covers days 4 through 23 and days 2 through 16 after the 2008 and 2009 events, respectively, the source spectral shape remained stable, while the pulsar's spin-down rate in the same span in 2008 increased by a factor of 2.2 as measured by RXTE. The lack of spectral variation suggests decoupling between magnetar spin-down and radiative changes, hence between the spin-down-inferred magnetic field strength and that inferred spectrally. We also found a strong anti-correlation between the phase-averaged flux and the pulsed fraction in the 2008 and 2009 Chandra data, but not in the pre-2008 measurements. We discuss these results in the context of the magnetar model.

After three years of no unusual activity, Anomalous X-ray Pulsar 1E 1048.1-5937 reactivated in 2007 March. We report on the detection of a large glitch (Delta(nu)/nu =1.63(2)X~10^{-5}) on 2007 March 26 (MJD 54185.9), contemporaneous with the onset of a pulsed-flux flare, the third flare observed from this source in 10 years of monitoring with the Rossi X-ray Timing Explorer. Additionally, we report on a detailed study of the evolution of the timing properties, the pulsed flux, and the pulse profile of this source as measured by RXTE from 1996 July to 2008 January. In our timing study, we attempted phase coherent timing of all available observations. We show that in 2001, a timing anomaly of uncertain nature occurred near the rise of the first pulsed flux flare; we show that a likely glitch (Delta(nu)/nu =2.91(9)X10^{-6}) occurred in 2002, near the rise of the second flare, and we present a detailed description of the variations in the spin-down. In our pulsed flux study, we compare the decays of the three flares and discuss changes in the hardness ratio. In our pulse profile study, we show that the profile exhibited large variations near the peak of the first two flares, and several small short-term profile variations during the most recent flare. Finally, we report on the discovery of a small burst 27 days after the peak of the last flare, the fourth burst discovered from this source. We discuss the relationships between the observed properties in the framework of the magnetar model.

We present high-speed optical photometry of the anomalous X-ray pulsar 1E 1048.1-5937 obtained with ULTRACAM on the 8.2-m Very Large Telescope in June 2007. We detect 1E 1048.1-5937 at a magnitude of i'=25.3+/-0.2, consistent with the values found by Wang et al. (2008) and hence confirming their conclusion that the source was approximately 1 mag brighter than in 2003-2006 due to an on-going X-ray flare that started in March 2007. The increased source brightness enabled us to detect optical pulsations with an identical period (6.458 s) to the X-ray pulsations. The rms pulsed fraction in our data is 21+/-7%, approximately the same as the 2-10 keV X-ray rms pulsed fraction. The optical and X-ray pulse profiles show similar morphologies and appear to be approximately in phase with each other, the latter lagging the former by only 0.06+/-0.02 cycles. The optical pulsations in 1E 1048.1-5937 are very similar in nature to those observed in 4U 0142+61. The implications of our observations for models of anomalous X-ray pulsars are discussed.

We present an analysis of five X-ray Multi-Mirror Mission (XMM) observations of the anomalous X-ray pulsar (AXP) 1E 2259+586 taken in 2004 and 2005 during its relaxation following its 2002 outburst. We compare these data with those of five previous XMM observations taken in 2002 and 2003, and find the observed flux decay is well described by a power-law of index -0.69+/-0.03. As of mid-2005, the source may still have been brighter than preoutburst, and was certainly hotter. We find a strong correlation between hardness and flux, as seen in other AXP outbursts. We discuss the implications of these results for the magnetar model.

We present detailed spectral and temporal characteristics both in the hard X-ray (>10 keV) and soft X-ray (<10 keV) domains, obtained using data from INTEGRAL, XMM-Newton, ASCA and RXTE. The INTEGRAL time-averaged total spectrum shows a power-law like shape with photon index Gamma = 0.93 +/- 0.06. 4U 0142+61 is detected up to 229 keV and the flux between 20 keV and 229 keV is (15.01 +/- 0.82) x 10^(-11) erg/cm^2/s. Using simultaneously collected data with the spectrometer SPI of INTEGRAL the combined total spectrum yields the first evidence for a spectral break with a peak energy of 228 +65/-41 keV. There is no evidence for significant long-term time variability of the total emission. Pulsed emission is measured with ISGRI up to 160 keV. The 20-160 keV profile shows a broad double-peaked pulse with a 6.2 sigma detection significance. The total pulsed spectrum can be described with a very hard power-law shape with a photon index \Gamma = 0.40 +/- 0.15. We performed phase-resolved spectroscopy over the total high-energy band (2.8-300 keV) and identify at least three genuinely different pulse components with different spectra. The high level of consistency between the detailed results from the four missions is indicative for a remarkable stable geometry underlying the emission scenario.

After 6 years of quiescence, Anomalous X-ray Pulsar (AXP) 4U 0142+61 entered an active phase in 2006 March that lasted several months. During the active phase, several bursts were detected, and many aspects of the X-ray emission changed. We report on the discovery of six X-ray bursts, the first ever seen from this AXP in ~10 years of Rossi X-ray Timing Explorer (RXTE) monitoring. All the bursts occurred in the interval between 2006 April 6 and 2007 February 7. The bursts had the canonical fast rise slow decay profiles characteristic of SGR/AXP bursts. The burst durations ranged from 8-3x10^3 s as characterized by T90,these are very long durations even when compared to the broad T90 distributions of other bursts from SGRs and AXPs. The first five burst spectra are well modeled by simple blackbodies, with temperature kT ~2-6 keV. However, the sixth burst had a complicated spectrum consisting of at least three emission lines with possible additional emission and absorption lines. The most significant feature was at ~14 keV. Similar 14-keV spectral features were seen in bursts from AXPs 1E 1048.1-5937 and XTE J1810-197. If this feature is interpreted as a proton cyclotron line, then it supports the existence of a magnetar-strength field for these AXPs. Several of the bursts were accompanied by a short-term pulsed flux enhancement. We discuss these events in the context of the magnetar model.

Anomalous X-ray Pulsars (AXPs), thought to be magnetars, exhibit poorly understood deviations from a simple spin-down called "timing noise". AXP timing noise has strong low-frequency components which pose significant challenges for quantification. We describe a procedure for extracting two quantities of interest, the intensity and power spectral index of timing noise. We apply this procedure to timing data from three sources: a monitoring campaign of five AXPs, observations of five young pulsars, and the stable rotator PSR B1937+21.

We present the results of X-ray and near-IR observations of the anomalous X-ray pulsar 1E 1048.1-5937, believed to be a magnetar. This AXP underwent a period of extreme variability during 2001-2004, but subsequently entered an extended and unexpected quiescence in 2004-2006, during which we monitored it with RXTE, CXO, and HST. Its timing properties were stable for >3 years throughout the quiescent period. 1E 1048.1-5937 again went into outburst in March 2007, which saw a factor of >7 total X-ray flux increase which was anti-correlated with a pulsed fraction decrease, and correlated with spectral hardening, among other effects. The near-IR counterpart also brightened following the 2007 event. We discuss our findings in the context of the magnetar and other models.

(Abridged) We report on new and archival X-ray and near-infrared observations of the anomalous X-ray pulsar 1E 1048.1-5937 performed between 2001-2007 with RXTE, CXO, Swift, HST, and VLT. During its ~2001-2004 active period, 1E 1048.-5937 exhibited two large, long-term X-ray pulsed-flux flares as well as short bursts, and large (>10x) torque changes. Monitoring with RXTE revealed that the source entered a phase of timing stability in 2004; at the same time, a series of four simultaneous observations with CXO and HST in 2006 showed that its X-ray flux and spectrum and near-IR flux, all variable prior to 2005, stabilized. The near-IR flux, when detected by HST (H~22.7 mag) and VLT (K_S~21.0 mag), was considerably fainter than previously measured. Recently, in 2007 March, this newfound quiescence was interrupted by a sudden flux enhancement, X-ray spectral changes and a pulse morphology change, simultaneous with a large spin-up glitch and near-IR enhancement. Our RXTE observations revealed a sudden pulsed flux increase by a factor of ~3 in the 2-10 keV band. In observations with CXO and Swift, we found that the total X-ray flux increased much more than the pulsed flux, reaching a peak value of >7 times the quiescent value (2-10 keV). With these recent data, we find a strong anti-correlation between X-ray flux and pulsed fraction, and a correlation between X-ray spectral hardness and flux. Simultaneously with the radiative and timing changes, we observed a significant X-ray pulse morphology change such that the profile went from nearly sinusoidal to having multiple peaks. We compare these remarkable events with other AXP outbursts and discuss implications in the context of the magnetar model and other models of AXP emission.

(Abridged) We report on 8.7 and 7.6 yr of RXTE observations of the Anomalous X-ray Pulsars (AXPs) RXS J170849.0-400910 and 1E 1841-045, respectively. These observations, part of a larger RXTE AXP monitoring program, have allowed us to study the long-term timing, pulsed flux, and pulse profile evolution of these objects. We report on four new glitches, one from RXS J170849.0-400910 and three from 1E 1841-045. With nearly all known persistent AXPs now seen to glitch, such behavior is clearly generic to this source class. We show that in terms of fractional frequency change, AXPs are among the most actively glitching neutron stars. However, in terms of absolute glitch amplitude, AXP glitches are unremarkable. We show that the largest AXP glitches observed thus far have recoveries that are unusual among those of radio pulsar glitches, with the combination of recovery time scale and fraction yielding changes in spin-down rates following the glitch similar to, or larger than, the long-term average. We also observed a large long-term fractional increase in the magnitude of the spin-down rate of 1E 1841-045 following its largest glitch. These observations are challenging to interpret in standard glitch models, as is the frequent occurence of large glitches given AXPs' high measured temperatures. We speculate that the stellar core may be involved in the largest AXP glitches. Furthermore, we show that AXP glitches appear to fall in two classes: radiatively loud and radiatively quiet. The latter, of which the glitches of J170849.0-400910 and 1E 1841-045 are examples, show little evidence for an accompanying radiative event. We also show, however, that pulse profile and pulsed flux changes are common in these AXPs, but do not apprear closely correlated with any timing behavior.

We present results obtained from X-ray observations of the anomalous X-ray pulsar (AXP) 4U 0142+61 taken between 2000-2007 using XMM-Newton, Chandra and Swift. In observations taken before 2006, the pulse profile is observed to become more sinusoidal and the pulsed fraction increased with time. These results confirm those derived using the Rossi X-ray Timing Explorer and expand the observed evolution to energies below 2 keV. The XMM-Newton total flux in the 0.5-10 keV band is observed to be nearly constant in observations taken before 2006, while an increase of ~10% is seen afterwards and coincides with the burst activity detected from the source in 2006-2007. After these bursts, the evolution towards more sinusoidal pulse profiles ceased while the pulsed fraction showed a further increase. No evidence for large-scale, long-term changes in the emission as a result of the bursts is seen. The data also suggest a correlation between the flux and hardness of the spectrum, with brighter observations on average having a harder spectrum. As pointed out by other authors, we find that the standard blackbody plus power-law model does not provide the best spectral fit to the emission from 4U 0142+61. We also report on observations taken with the Gemini telescope after two bursts. These observations show source magnitudes consistent with previous measurements. Our results demonstrate the wide range of X-ray variability characteristics seen in AXPs and we discuss them in light of current emission models for these sources.

We report on 10 years of monitoring of the 8.7-s Anomalous X-ray Pulsar 4U 0142+61 using the Rossi X-Ray Timing Explorer (RXTE). This pulsar exhibited stable rotation from 2000 March until 2006 February: the RMS phase residual for a spin-down model which includes nu, nudot, and nuddot is 2.3%. We report a possible phase-coherent timing solution valid over a 10-yr span extending back to March 1996. A glitch may have occured between 1998 and 2000, but is not required by the existing timing data. The pulse profile has been evolving since 2000. In particular, the dip of emission between its two peaks got shallower between 2002 and 2006, as if the profile were evolving back to its pre-2000 morphology, following an earlier event, which possibly also included the glitch suggested by the timing data. These profile variations are seen in the 2-4 keV band but not in 6-8 keV. We also detect a slow increase in the pulsed flux between 2002 May and 2004 December, such that it has risen by 36+/-3% over 2.6 years in the 2-10 keV band. The pulsed flux variability and the narrow-band pulse profile changes present interesting challenges to aspects of the magnetar model.

For the first time a quasi-simultaneous multi-wavelength campaign has been performed on an Anomalous X-ray Pulsar from the radio to the hard X-ray band. 4U 0142+61 was an INTEGRAL target for 1 Ms in July 2005. During these observations it was also observed in the X-ray band with Swift and RXTE, in the optical and NIR with Gemini North and in the radio with the WSRT. In this paper we present the source-energy distribution. The spectral results obtained in the individual wave bands do not connect smoothly; apparently components of different origin contribute to the total spectrum. Remarkable is that the INTEGRAL hard X-ray spectrum (power-law index 0.79 +/- 0.10) is now measured up to an energy of ~230 keV with no indication of a spectral break. Extrapolation of the INTEGRAL power-law spectrum to lower energies passes orders of magnitude underneath the NIR and optical fluxes, as well as the low ~30 microJy (2 sigma) upper limit in the radio band.

We report on 10 yr of monitoring of the 8.7-s Anomalous X-ray Pulsar 4U 0142+61 using the Rossi X-Ray Timing Explorer (RXTE). This pulsar exhibited stable rotation from 2000 until February 2006: the RMS phase residual for a spin-down model which includes nu, nudot, and nuddot is 2.3%. We report a possible phase-coherent timing solution valid over a 10-yr span extending back to March 1996. A glitch may have occured between 1998 and 2000, but it is not required by the existing data. We also report that the source's pulse profile has been evolving in the past 6 years, such that the dip of emission between its two peaks has been getting shallower since 2000, almost as if the profile is recovering to its pre-2000 morphology, in which there was no clear distinction between the peaks. These profile variations are seen in the 2-4 keV band but not in 6-8 keV. Finally, we present the pulsed flux time series of the source in 2-10 keV. There is evidence of a slow but steady increase in the source's pulsed flux since 2000. The pulsed flux variability and the narrow-band pulse profile changes present interesting challenges to aspects of the magnetar model.

As part of a large-scale search for coherent pulsations from LMXBs in the RXTE archive, we have completed a detailed series of searches for coherent pulsations of 4U 1820-30 -- an ultracompact LMXB with a binary period of 11.4 min, located in the globular cluster NGC6624. The short binary period implies any coherent signal would be highly accelerated, so we used phase modulation searches, orbital-parameter-fitting coherent searches, and standard acceleration searches to give significant sensitivity to millisecond pulsations. We searched, in four energy bands and at a range of luminosities, a total of 34 archival RXTE observations, 32 of which had on-source integration times longer than 10 ks, and some of which were made consecutively which allowed us to combine them. We found no pulsations. Using our phase modulation search technique, which we ran on all 34 observations, we have been able to place the first stringent (95% confidence) pulsed fraction limits of <~0.8% for all realistic spin frequencies (i.e. <~2kHz) and likely companion masses (0.02Msun <= M_c <= 0.3Msun). Using our orbital-parameter-fitting coherent search, which we ran on only 11 selected observations, we have placed a pulsed fraction limit of <~0.3% for spin frequencies <~1.25kHz and companion masses M_ <= 0.106Msun. By contrast, all five LMXBs known to emit coherent pulsations have intrinsic pulsed fractions in the range 3% to 7% when pulsations are observed. Hence, our searches rule out pulsations with significantly lower pulsed fractions than those already observed.

As part of a large-scale search for coherent pulsations from LMXBs in the RXTE archive, we have completed a detailed series of searches for coherent pulsations of 4U 1820-30 -- an ultracompact LMXB with a binary period of 11.4 minutes located in the globular cluster NGC6624. The small binary period leads to a very high acceleration, so we used phase modulation searches as well as acceleration searches to give significant sensitivity to millisecond pulsations. We searched a total of 34 archival RXTE observations, 32 of which had an on-source integration time longer than 10 ks, and some of which were made consecutively which allowed us to combine them. While we found no pulsations, we have been able to place the first stringent (95% confidence) pulsed fraction limits of < 0.8% for all realistic spin frequencies (i.e. < 1 kHz) and likely companion masses (0.02 Msun < Mc < 0.3 Msun). By contrast all five LMXBs known to emit coherent pulsations have intrinsic pulsed fractions in the range 3% to 7% when pulsations are observed.