The Cassini and Juno probes have revealed large coherent cyclonic vortices in the polar regions of Saturn and Jupiter, a dramatic contrast from the east-west banded jet structure seen at lower latitudes. Debate has centered on whether the jets are shallow, or extend to greater depths in the planetary envelope. Recent experiments and observations have demonstrated the relevance of deep convection models to a successful explanation of jet structure and cyclonic coherent vortices away from the polar regions have been simulated recently including an additional stratified shallow layer. Here we present new convective models able to produce long-lived polar vortices. Using simulation parameters relevant for giant planet atmospheres we find flow regimes that are in agreement with geostrophic turbulence (GT) theory in rotating convection for the formation of large scale coherent structures via an upscale energy transfer fully three-dimensional. Our simulations generate polar characteristics qualitatively similar to those seen by Juno and Cassini: they match the structure of cyclonic vortices seen on Jupiter; or can account for the existence of a strong polar vortex extending downwards to lower latitudes with a marked spiral morphology and the hexagonal pattern seen on Saturn. Our findings indicate that these vortices can be generated deep in the planetary interior. A transition differentiating these two polar flows regimes is described, interpreted in terms of different force balances and compared with previous shallow atmospheric models which characterised polar vortex dynamics in giant planets. In addition, the heat transport properties are investigated confirming recent scaling laws obtained in the context of reduced models of GT.

Waves in thin fluid layers are important in various stellar and planetary problems. Due to rapid rotation such systems will become oblate, with a latitudinal variation in the gravitational acceleration across the surface of the object. In the case of accreting neutron stars, rapid rotation could lead to a polar radius smaller than the equatorial radius by a factor $\sim 0.8$. We investigate how the oblateness and a changing gravitational acceleration affect different hydrodynamic modes that exist in such fluid layers through analytic approximations and numerical calculations. The wave vectors of $g$-modes and Yanai modes increase for more oblate systems compared to spherical counterparts, although the impact of variations in the changing gravitational acceleration is effectively negligible. We find that for increased oblateness, Kelvin modes show less equatorial confinement and little change in their wave vector. For $r$-modes, we find that for more oblate systems the wave vector decreases. The exact manner of these changes for the $r$-modes depends on the model for the gravitational acceleration across the surface.

We present the largest sample of type-I (thermonuclear) X-ray bursts yet assembled, comprising 7083 bursts from 85 bursting sources. The sample is drawn from observations with Xenon-filled proportional counters on the long-duration satellites RXTE, BeppoSAX and INTEGRAL, between 1996 February 8, and 2012 May 3. The burst sources were drawn from a comprehensive catalog of 115 burst sources, assembled from earlier catalogs and the literature. We carried out a consistent analysis for each burst lightcurve (normalised to the relative instrumental effective area), and provide measurements of rise time, peak intensity, burst timescale, and fluence. For bursts observed with the RXTE/PCA and BeppoSAX/WFC we also provide time-resolved spectroscopy, including estimates of bolometric peak flux and fluence, and spectral parameters at the peak of the burst. For 950 bursts observed with the PCA from sources with previously detected burst oscillations, we include an analysis of the high-time resolution data, providing information on the detectability and amplitude of the oscillations, as well as where in the burst they are found. We also present analysis of 117883 observations of the burst sources within the sample timeframe. We extracted 3-25 keV X-ray spectra from most observations, and (for observations meeting our signal-to-noise criterion), we provide measurements of the flux, spectral colours, and for selected sources, the position on the colour-colour diagram, for the best-fit spectral model. We present a description of the sample, a summary of the science investigations completed to date, and suggestions for further studies.

The NICER collaboration recently published a joint estimate of the mass and the radius of PSR J0030+0451, derived via X-ray pulse-profile modeling. Raaijmakers et al. (2019) explored the implications of this measurement for the dense matter equation of state (EOS) using two parameterizations of the high-density EOS: a piecewise-polytropic model, and a model based on the speed of sound in neutron stars. In this work we obtain further constraints on the EOS following this approach, but we also include information about the tidal deformability of neutron stars from the gravitational wave signal of the compact binary merger GW170817. We compare the constraints on the EOS to those set by the recent measurement of a 2.14 solar mass pulsar, included as a likelihood function approximated by a Gaussian, and find a small increase in information gain. To show the flexibility of our method, we also explore the possibility that GW170817 was a neutron star-black hole merger, which yields weaker constraints on the EOS.

We present an analysis of asymmetries in global HI spectra from the extended GALEX Arecibo SDSS Survey (xGASS), a stellar mass-selected and gas fraction-limited survey which is representative of the HI properties of galaxies in the local Universe. We demonstrate that the asymmetry in a HI spectrum is strongly linked to its signal-to-noise meaning that, contrary to what was done in previous works, asymmetry distributions for different samples cannot be compared at face value. We develop a method to account for noise-induced asymmetry and find that the typical galaxy detected by xGASS exhibits higher asymmetry than what can be attributed to noise alone, with 37% of the sample showing asymmetry greater than 10% at an 80% confidence level. We find that asymmetric galaxies contain, on average, 29% less HI mass compared to their symmetric counterparts matched in both stellar mass and signal-to-noise. We also present clear evidence that satellite galaxies, as a population, exhibit more asymmetric HI spectra than centrals and that group central galaxies show a slightly higher rate of HI asymmetries compared to isolated centrals. All these results support a scenario in which environmental processes, in particular those responsible for gas removal, are the dominant driver of asymmetry in xGASS.

We report on Bayesian parameter estimation of the mass and equatorial radius of the millisecond pulsar PSR J0030$+$0451, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer (NICER) X-ray spectral-timing event data. We perform relativistic ray-tracing of thermal emission from hot regions of the pulsar's surface. We assume two distinct hot regions based on two clear pulsed components in the phase-folded pulse-profile data; we explore a number of forms (morphologies and topologies) for each hot region, inferring their parameters in addition to the stellar mass and radius. For the family of models considered, the evidence (prior predictive probability of the data) strongly favors a model that permits both hot regions to be located in the same rotational hemisphere. Models wherein both hot regions are assumed to be simply-connected circular single-temperature spots, in particular those where the spots are assumed to be reflection-symmetric with respect to the stellar origin, are strongly disfavored. For the inferred configuration, one hot region subtends an angular extent of only a few degrees (in spherical coordinates with origin at the stellar center) and we are insensitive to other structural details; the second hot region is far more azimuthally extended in the form of a narrow arc, thus requiring a larger number of parameters to describe. The inferred mass $M$ and equatorial radius $R_\mathrm{eq}$ are, respectively, $1.34_{-0.16}^{+0.15}$ M$_{\odot}$ and $12.71_{-1.19}^{+1.14}$ km, whilst the compactness $GM/R_\mathrm{eq}c^2 = 0.156_{-0.010}^{+0.008}$ is more tightly constrained; the credible interval bounds reported here are approximately the $16\%$ and $84\%$ quantiles in marginal posterior mass.

Both the mass and radius of the millisecond pulsar PSR J0030+0451 have been inferred via pulse-profile modeling of X-ray data obtained by NASA's NICER mission. In this Letter we study the implications of the mass-radius inference reported for this source by Riley et al. (2019) for the dense matter equation of state (EOS), in the context of prior information from nuclear physics at low densities. Using a Bayesian framework we infer central densities and EOS properties for two choices of high-density extensions: a piecewise-polytropic model and a model based on assumptions of the speed of sound in dense matter. Around nuclear saturation density these extensions are matched to an EOS uncertainty band obtained from calculations based on chiral effective field theory interactions, which provide a realistic description of atomic nuclei as well as empirical nuclear matter properties within uncertainties. We further constrain EOS expectations with input from the current highest measured pulsar mass; together, these constraints offer a narrow Bayesian prior informed by theory as well as laboratory and astrophysical measurements. The NICER mass-radius likelihood function derived by Riley et al. (2019) using pulse-profile modeling is consistent with the highest-density region of this prior. The present relatively large uncertainties on mass and radius for PSR J0030+0451 offer, however, only a weak posterior information gain over the prior. We explore the sensitivity to the inferred geometry of the heated regions that give rise to the pulsed emission, and find a small increase in posterior gain for an alternative (but less preferred) model. Lastly, we investigate the hypothetical scenario of increasing the NICER exposure time for PSR J0030+0451.

We describe the model of surface emission from a rapidly rotating neutron star that is applied to Neutron Star Interior Composition Explorer X-ray data of millisecond pulsars in order to statistically constrain the neutron star mass-radius relation and dense matter equation of state. To ensure that the associated calculations are both accurate and precise, we conduct an extensive suite of verification tests between our numerical codes for both the Schwarzschild + Doppler and Oblate Schwarzschild approximations, and compare both approximations against exact numerical calculations. We find superb agreement between the code outputs, as well as in comparison against a set of analytical and semi-analytical calculations, which combined with their speed, demonstrates that the codes are well-suited for large-scale statistical sampling applications. A set of verified, high-precision reference synthetic pulse profiles is provided to the community to facilitate testing of other independently developed codes.

We present the set of deep Neutron Star Interior Composition Explorer (NICER) X-ray timing observations of the nearby rotation-powered millisecond pulsars PSRs J0437-4715, J0030+0451, J1231-1411, and J2124-3358, selected as targets for constraining the mass-radius relation of neutron stars and the dense matter equation of state via modeling of their pulsed thermal X-ray emission. We describe the instrument, observations, and data processing/reduction procedures, as well as the series of investigations conducted to ensure that the properties of the data sets are suitable for parameter estimation analyses to produce reliable constraints on the neutron star mass-radius relation and the dense matter equation of state. We find that the long-term timing and flux behavior and the Fourier-domain properties of the event data do not exhibit any anomalies that could adversely affect the intended measurements. From phase-selected spectroscopy, we find that emission from the individual pulse peaks is well described by a single-temperature hydrogen atmosphere spectrum, with the exception of PSR J0437-4715, for which multiple temperatures are required.

Recent modeling of NICER observations of thermal X-ray pulsations from the surface of the isolated millisecond pulsar PSR J0030+0451 suggests that the hot emitting regions on the pulsar's surface are far from antipodal, which is at odds with the classical assumption that the magnetic field in the pulsar magnetosphere is predominantly that of a centered dipole. Here, we review these results and examine previous attempts to constrain the magnetospheric configuration of PSR J0030+0451. To the best of our knowledge, there is in fact no direct observational evidence that PSR J0030+0451's magnetic field is a centered dipole. Developing models of physically motivated, non-canonical magnetic field configurations and the currents that they can support poses a challenging task. However, such models may have profound implications for many aspects of pulsar research, including pulsar braking, estimates of birth velocities, and interpretations of multi-wavelength magnetospheric emission.

Accreting neutron stars (NS) can exhibit high frequency modulations in their lightcurves during thermonuclear X-ray bursts, known as burst oscillations. These frequencies can be offset from the NS spin frequency by several Hz (where known independently) and can drift by 1-3 Hz. One plausible explanation is that a wave is present in the bursting ocean, the rotating frame frequency of which is the offset. The frequency of the wave should decrease (in the rotating frame) as the burst cools hence explaining the drift. A strong candidate is a buoyant $r$-mode. To date, models that calculated the frequency of this mode taking into account the radial structure neglected relativistic effects and predicted rotating frame frequencies of $\sim$ 4 Hz and frequency drifts of > 5 Hz; too large to be consistent with observations. We present a calculation that includes frame-dragging and gravitational redshift that reduces the rotating frame frequency by up to 30 % and frequency drift by up to 20 %. Updating previous models for the ocean cooling in the aftermath of the burst to a model more representative of detailed calculations of thermonuclear X-ray bursts reduces the frequency of the mode still further. This model, combined with relativistic effects, can reduce the rotating frequency of the mode to $\sim$ 2 Hz and frequency drift to $\sim$ 2 Hz, which is closer to the observed values.

Using Hubble Space Telescope ACS/WFC data we present the photometry and spatial distribution of resolved stellar populations of four fields within the extended ultraviolet disk (XUV disk) of M83. These observations show a clumpy distribution of main-sequence stars and a mostly smooth distribution of red giant branch stars. We constrain the upper-end of the initial mass function (IMF) in the outer disk using the detected population of main-sequence stars and an assumed constant star formation rate (SFR) over the last 300 Myr. By comparing the observed main-sequence luminosity function to simulations, we determine the best-fitting IMF to have a power law slope $\alpha=-2.35 \pm 0.3$ and an upper-mass limit $\rm M_{u}=25_{-3}^{+17} \, M_\odot$. This IMF is consistent with the observed H$\alpha$ emission, which we use to provide additional constraints on the IMF. We explore the influence of deviations from the constant SFR assumption, finding that our IMF conclusions are robust against all but strong recent variations in SFR, but these are excluded by causality arguments. These results, along with our similar studies of other nearby galaxies, indicate that some XUV disks are deficient in high-mass stars compared to a Kroupa IMF. There are over one hundred galaxies within 5 Mpc, many already observed with HST, thus allowing a more comprehensive investigation of the IMF, and how it varies, using the techniques developed here.

Convection in rotating spherical geometries is an important physical process in planetary and stellar systems. Using continuation methods at low Prandtl number, we find both strong equatorially asymmetric and symmetric polar nonlinear rotating waves in a model of thermal convection in thin rotating spherical shells with stress-free boundary conditions. For the symmetric waves convection is confined to high latitude in both hemispheres but is only restricted to one hemisphere close to the pole in the case of asymmetric waves. This is in contrast to what is previously known from studies in the field. These periodic flows, in which the pattern is rotating steadily in the azimuthal direction, develop a strong axisymmetric component very close to onset. Using stability analysis of periodic orbits the regions of stability are determined and the topology of the stable/unstable oscillatory flows bifurcated from the branches of rotating waves is described. By means of direct numerical simulations of these oscillatory chaotic flows, we show that these three-dimensional convective polar flows exhibit characteristics, such as force balance or mean physical properties, which are similar to flows occuring in planetary atmospheres. We show that these results may open a route to understanding unexplained features of gas giant atmospheres, in particular for the case of Jupiter. These include the observed equatorial asymmetry with a pronounced decrease at the equator (the so-called dimple), and the coherent vortices surrounding the poles recently observed by the Juno mission.

We describe a blind uniform search for thermonuclear burst oscillations (TBOs) in the majority of Type-I bursts observed by RXTE (2118 bursts from 57 neutron stars). We examined 2-2002 Hz power spectra from the Fourier transform in sliding 0.5-2 s windows, using fine-binned light curves in 2-60 keV energy range. The significance of the oscillation candidates was assessed by simulations which took into account light curve variations, dead time and sliding time windows. Some of our sources exhibited multi-frequency variability below approximately 15 Hz that cannot be readily removed with light-curve modeling and may have an astrophysical (non-TBO) nature. Overall, we found that the number and strength of potential candidates depends strongly on the parameters of the search. We found candidates from all previously known RXTE TBO sources, with pulsations that had been detected at similar frequencies in multiple independent time windows, and discovered TBOs from SAX J1810.8-2658. We could not confirm most previously-reported tentative TBO detections or identify any obvious candidates just below the detection threshold at similar frequencies in multiple bursts. We computed fractional amplitudes of all TBO candidates and placed upper limits on non-detections. Finally, for a few sources we noted small excess of candidates with powers comparable to fainter TBOs, but appearing in single independent time bins at random frequencies. At least some of these candidates may be noise spikes that appear interesting due to selection effects. The potential presence of such candidates calls for extra caution if claiming single-bin TBO detections.

One very promising technique for measuring the dense matter Equation of State exploits hotspots that form on the neutron star surface due to the pulsar mechanism, accretion streams, or during thermonuclear explosions in the neutron star ocean. This article explains how Pulse Profile Modeling of hotspots is being used by the Neutron Star Interior Composition Explorer (NICER), an X-ray telescope installed on the International Space Station in 2017 - and why the technique is a mission driver for the next, larger-area generation of telescopes including the enhanced X-ray Timing and Polarimetry (eXTP) mission and the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X).

Burst oscillations are brightness asymmetries that develop in the burning ocean during thermonuclear bursts on accreting neutron stars. They have been observed during H/He-triggered (Type I) bursts and Carbon-triggered superbursts. The mechanism responsible is not unknown, but the dominant burst oscillation frequency is typically within a few Hz of the spin frequency, where this is independently known. One of the best-studied burst oscillation sources, 4U 1636-536, has oscillations at $581\,\text{Hz}$ in both its regular Type I bursts and in one superburst. Recently however, Strohmayer \& Mahmoodifar reported the discovery of an additional signal at a higher frequency, $835\,\text{Hz}$, during the superburst. This higher frequency is consistent with the predictions for several types of global ocean mode, one of the possible burst oscillation mechanisms. If this is the case then the same physical mechanism may operate in the normal Type I bursts of this source. In this paper we report a stacked search for periodic signals in the regular Type I bursts: we found no significant signal at the higher frequency, with upper limits for the single trial root mean square (rms) fractional amplitude of 0.57(6)\%. Our analysis did however reveal that the dominant $581\,\text{Hz}$ burst oscillation signal is present at a weak level even in the sample of bursts where it cannot be detected in individual bursts. This indicates that any cutoff in the burst oscillation mechanism occurs below the detection threshold of existing X-ray telescopes.

The unknown state of matter at ultra-high density, large proton/neutron number asymmetry, and low temperature is a major long-standing problem in modern physics. Neutron stars provide the only known setting in the Universe where matter in this regime can stably exist. Valuable information about the interior structure of neutron stars can be extracted via sensitive observations of their exteriors. There are several complementary techniques that require different combinations of high time resolution, superb spectral resolution, and high spatial resolution. In the upcoming decade and beyond, measurements of the masses and radii of an ensemble of neutron stars using these techniques, based on data from multiple proposed next-generation X-ray telescopes, can produce definitive empirical constraints on the allowed dense matter equation of state.

We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science.

Understanding the dense matter equation of state at extreme conditions is an important open problem. Astrophysical observations of neutron stars promise to solve this, with NICER poised to make precision measurements of mass and radius for several stars using the waveform modelling technique. What has been less clear, however, is how these mass-radius measurements might translate into equation of state constraints and what are the associated equation of state sensitivities. We use Bayesian inference to explore and contrast the constraints that would result from different choices for the equation of state parametrization; comparing the well-established piecewise polytropic parametrization to one based on physically motivated assumptions for the speed of sound in dense matter. We also compare the constraints resulting from Bayesian inference to those from simple compatibility cuts. We find that the choice of equation of state parametrization and particularly its prior assumptions can have a significant effect on the inferred global mass-radius relation and the equation of state constraints. Our results point to important sensitivities when inferring neutron star and dense matter properties. This applies also to inferences from gravitational wave observations.

We present an extensive study of Boussinesq thermal convection including a temperature-dependent internal heating source, based on numerical three-dimensional simulations. The temperature dependence mimics triple-$\alpha$ nuclear reactions and the fluid geometry is a rotating spherical shell. These are key ingredients for the study of convective accreting neutron star oceans. A dimensionless parameter ${{\rm Ra}}_n$, measuring the relevance of nuclear heating, is defined. We explore how flow characteristics change with increasing ${{\rm Ra}}_n$ and give an astrophysical motivation. The onset of convection is investigated with respect to this parameter and periodic, quasiperiodic, chaotic flows with coherent structures, and fully turbulent flows are exhibited as ${{\rm Ra}}_n$ is varied. Several regime transitions are identified and compared with previous results on differentially heated convection. Finally, we explore (tentatively) the potential applicability of our results to the evolution of thermonuclear bursts in accreting neutron star oceans.

In this paper we present the science potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies of strongly magnetized objects. We will focus on the physics and astrophysics of strongly magnetized objects, namely magnetars, accreting X-ray pulsars, and rotation powered pulsars. We also discuss the science potential of eXTP for QED studies. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Sciences, the eXTP mission is expected to be launched in the mid 2020s.

In this White Paper we present the potential of the Enhanced X-ray Timing and Polarimetry (eXTP) mission for determining the nature of dense matter; neutron star cores host an extreme density regime which cannot be replicated in a terrestrial laboratory. The tightest statistical constraints on the dense matter equation of state will come from pulse profile modelling of accretion-powered pulsars, burst oscillation sources, and rotation-powered pulsars. Additional constraints will derive from spin measurements, burst spectra, and properties of the accretion flows in the vicinity of the neutron star. Under development by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

In this paper we present the enhanced X-ray Timing and Polarimetry mission - eXTP. eXTP is a space science mission designed to study fundamental physics under extreme conditions of density, gravity and magnetism. The mission aims at determining the equation of state of matter at supra-nuclear density, measuring effects of QED, and understanding the dynamics of matter in strong-field gravity. In addition to investigating fundamental physics, eXTP will be a very powerful observatory for astrophysics that will provide observations of unprecedented quality on a variety of galactic and extragalactic objects. In particular, its wide field monitoring capabilities will be highly instrumental to detect the electro-magnetic counterparts of gravitational wave sources. The paper provides a detailed description of: (1) the technological and technical aspects, and the expected performance of the instruments of the scientific payload; (2) the elements and functions of the mission, from the spacecraft to the ground segment.

In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to study one common aspect of these objects: their often transient nature. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

Accreting neutron stars (NS) can exhibit high-frequency modulations, known as burst oscillations, in their lightcurves during thermonuclear X-ray bursts. Their frequencies can be offset from the spin frequency of the NS (known independently) by several Hz, and can drift by 1-3 Hz. One plausible explanation for this phenomenon is that a wave is present in the bursting ocean that decreases in frequency (in the rotating frame) as the burst cools. The strongest candidate is the buoyant $r$-mode; however, models for the burning ocean background used in previous studies over-predict frequency drifts by several Hz. Using new background models (which include shallow heating, and burning in the tail of the burst) the evolution of the buoyant $r$-mode is calculated. The resulting frequency drifts are smaller, in line with observations. This illustrates the importance of accounting for the detailed nuclear physics in these bursts.

We describe the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept that will provide an unprecedented view of the X-ray sky, performing timing and spectroscopy over both a broad energy band (0.2-30 keV) and a wide range of timescales from microseconds to years. STROBE-X comprises two narrow-field instruments and a wide field monitor. The soft or low-energy band (0.2-12 keV) is covered by an array of lightweight optics (3-m focal length) that concentrate incident photons onto small solid-state detectors with CCD-level (85-175 eV) energy resolution, 100 ns time resolution, and low background rates. This technology has been fully developed for NICER and will be scaled up to take advantage of the longer focal length of STROBE-X. The higher-energy band (2-30 keV) is covered by large-area, collimated silicon drift detectors that were developed for the European LOFT mission concept. Each instrument will provide an order of magnitude improvement in effective area over its predecessor (NICER in the soft band and RXTE in the hard band). Finally, STROBE-X offers a sensitive wide-field monitor (WFM), both to act as a trigger for pointed observations of X-ray transients and also to provide high duty-cycle, high time-resolution, and high spectral-resolution monitoring of the variable X-ray sky. The WFM will boast approximately 20 times the sensitivity of the RXTE All-Sky Monitor, enabling multi-wavelength and multi-messenger investigations with a large instantaneous field of view. This mission concept will be presented to the 2020 Decadal Survey for consideration.

SAX J1810.8-2609 is a faint X-ray transient, mostly known for its abnormally low quiescent thermal luminosity, which disagrees with standard cooling models. It is also one of a small sample of stars whose mass and radius have been estimated using spectral modeling of one of its thermonuclear bursts. Here we report the discovery of millisecond oscillations in a type I thermonuclear X-ray burst from SAX J1810.8-2609 observed by RXTE during the 2007 outburst. A strong signal (Leahy-normalized power of 71.5, 4.5e-9 chance of coincidence with a conservative estimate for the number of trials) was present at 531.8 Hz during the decay of one out of six bursts observed. Oscillations were detected for about 6 seconds, during which their frequency increased from 531.4 to 531.9 Hz in a manner similar to other burst oscillation sources. The millisecond oscillations discovered pinpoint the spin frequency of the neutron star, which is important for the spectral modeling, associated mass-radius inference, and understanding the evolutionary status and cooling behavior of the star. As of April 2018 the source is in outburst again, providing a fleeting opportunity to acquire new material for the burst oscillation searches.

This study considers the onset of stress-free Boussinesq thermal convection in rotating spherical shells with aspect ratio $\eta=r_i/r_o=0.9$ ($r_i$ and $r_o$ being the inner and outer radius), Prandtl numbers ${\rm Pr} \in[10^{-4},10^{-1}]$, and Taylor numbers ${\rm Ta}\in[10^{4},10^{12}]$. We are particularly interested in the form of the convective cell pattern that develops, and in its time scales, since this may have observational consequences. For a fixed ${\rm Ta}<10^{9}$ and by decreasing ${\rm Pr}$ from 0.1 to $10^{-4}$ a transition between spiralling columnar (SC) and equatorially-attached (EA) modes, and a transition between EA and equatorially antisymmetric or symmetric polar (AP/SP) weakly multicellular modes are found. The latter modes are preferred at very low ${\rm Pr}$. Surprisingly, for ${\rm Ta}>3\times 10^{9}$ the unicellular polar modes become also preferred at moderate ${\rm Pr}\sim10^{-2}$ because two new transition curves between EA and AP/SP and between AP/SP and SC modes are born at a triple-point bifurcation. The dependence on ${\rm Pr}$ and ${\rm Ta}$ of the transitions is studied to estimate the type of modes, and their critical parameters, preferred at different stellar regimes.

The only messenger radiation in the Universe which one can use to statistically probe the Equation of State (EOS) of cold dense matter is that originating from the near-field vicinities of compact stars. Constraining gravitational masses and equatorial radii of rotating compact stars is a major goal for current and future telescope missions, with a primary purpose of constraining the EOS. From a Bayesian perspective it is necessary to carefully discuss prior definition; in this context a complicating issue is that in practice there exist pathologies in the general relativistic mapping between spaces of local (interior source matter) and global (exterior spacetime) parameters. In a companion paper, these issues were raised on a theoretical basis. In this study we reproduce a probability transformation procedure from the literature in order to map a joint posterior distribution of Schwarzschild gravitational masses and radii into a joint posterior distribution of EOS parameters. We demonstrate computationally that EOS parameter inferences are sensitive to the choice to define a prior on a joint space of these masses and radii, instead of on a joint space interior source matter parameters. We focus on the piecewise-polytropic EOS model, which is currently standard in the field of astrophysical dense matter study. We discuss the implications of this issue for the field.

Constraining the equation of state of cold dense matter in compact stars is a major science goal for observing programmes being conducted using X-ray, radio, and gravitational wave telescopes. We discuss Bayesian hierarchical inference of parametrised dense matter equations of state. In particular we generalise and examine two inference paradigms from the literature: (i) direct posterior equation of state parameter estimation, conditioned on observations of a set of rotating compact stars; and (ii) indirect parameter estimation, via transformation of an intermediary joint posterior distribution of exterior spacetime parameters (such as gravitational masses and coordinate equatorial radii). We conclude that the former paradigm is not only tractable for large-scale analyses, but is principled and flexible from a Bayesian perspective whilst the latter paradigm is not. The thematic problem of Bayesian prior definition emerges as the crux of the difference between these paradigms. The second paradigm should in general only be considered as an ill-defined approach to the problem of utilising archival posterior constraints on exterior spacetime parameters; we advocate for an alternative approach whereby such information is repurposed as an approximative likelihood function. We also discuss why conditioning on a piecewise-polytropic equation of state model - currently standard in the field of dense matter study - can easily violate conditions required for transformation of a probability density distribution between spaces of exterior (spacetime) and interior (source matter) parameters.

We present optical photometry of Hubble Space Telescope (HST) ACS/WFC data of the resolved stellar populations in the outer disc of the dwarf irregular galaxy DDO 154. The photometry reveals that young main sequence stars are almost absent from the outermost HI disc. Instead, most are clustered near the main stellar component of the galaxy. We constrain the stellar initial mass function (IMF) by comparing the luminosity function of the main sequence stars to simulated stellar populations assuming a constant star formation rate over the dynamical timescale. The best-fitting IMF is deficient in high mass stars compared to a canonical Kroupa IMF, with a best-fit slope $\alpha = -2.45$ and upper mass limit $M_U = 16\ M_{\odot}$. This top-light IMF is consistent with predictions of the Integrated Galaxy-wide IMF theory. Combining the HST images with HI data from The HI Nearby Galaxy Survey Treasury (THINGS) we determine the star formation law (SFL) in the outer disc. The fit has a power law exponent $N = 2.92 \pm0.22$ and zero point $A=4.47 \pm 0.65 \times 10^{-7} \ M_{\odot} \ \text{yr}^{-1} \ \text{kpc}^{-2}$. This is depressed compared to the Kennicutt-Schmidt Star Formation Law, but consistent with weak star formation observed in diffuse HI environments. Extrapolating the SFL over the outer disc implies that there could be significant star formation occurring that is not detectable in H$\alpha$. Last, we determine the Toomre stability parameter $Q$ of the outer disc of DDO 154 using the THINGS HI rotation curve and velocity dispersion map. 72% of the HI in our field has $Q\leq 4$ and this incorporates 96% of the observed MS stars. Hence 28% of the HI in the field is largely dormant.

Runaway thermonuclear burning of a layer of accumulated fuel on the surface of a compact star provides a brief but intense display of stellar nuclear processes. For neutron stars accreting from a binary companion, these events manifest as thermonuclear (type-I) X-ray bursts, and recur on typical timescales of hours to days. We measured the burst rate as a function of accretion rate, from seven neutron stars with known spin rates, using a burst sample accumulated over several decades. At the highest accretion rates, the burst rate is lower for faster spinning stars. The observations imply that fast (> 400 Hz) rotation encourages stabilization of nuclear burning, suggesting a dynamical dependence of nuclear ignition on the spin rate. This dependence is unexpected, because faster rotation entails less shear between the surrounding accretion disk and the star. Large-scale circulation in the fuel layer, leading to enhanced mixing of the burst ashes into the fuel layer, may explain this behavior; further numerical simulations are required to confirm.

Accreting neutron stars (NS) can exhibit high frequency modulations in their lightcurves during thermonuclear X-ray bursts, known as burst oscillations. The frequencies can be offset from the spin frequency of the NS by several Hz, and can drift by 1-3 Hz. One possible explanation is a mode in the bursting ocean, the frequency of which would decrease (in the rotating frame) as the burst cools, hence explaining the drifts. Most burst oscillations have been observed during H/He triggered bursts, however there has been one observation of oscillations during a superburst; hours' long Type I X-ray bursts caused by unstable carbon burning deeper in the ocean. This paper calculates the frequency evolution of an oceanic r-mode during a superburst. The rotating frame frequency varies during the burst from 4-14 Hz, and is sensitive to the background parameters, in particular the temperature of the ocean and ignition depth. This calculation is compared to the superburst oscillations observed on 4U-1636-536. The predicted mode frequencies ($\sim$ 10 Hz) would require a spin frequency of $\sim$ 592 Hz to match observations; 6 Hz higher than the spin inferred from an oceanic r-mode model for the H/He triggered burst oscillations. This model also over-predicts the frequency drift during the superburst by 90 %.

The trigger for the short bursts observed in $\gamma$-rays from many magnetar sources remains unknown. One particular open question in this context is the localization of burst emission to a singular active region or a larger area across the neutron star. While several observational studies have attempted to investigate this question by looking at the phase dependence of burst properties, results have been mixed. At the same time, it is not obvious a priori that bursts from a localized active region would actually give rise to a detectable phase-dependence, taking into account issues such as geometry, relativistic effects, and intrinsic burst properties such brightness and duration. In this paper, we build a simple theoretical model to investigate the circumstances under which the latter effects could affect detectability of a dependence of burst emission on rotational phase. We find that even for strongly phase-dependent emission, inferred burst properties may not show a rotational phase dependence depending on the geometry of the system and the observer. Furthermore, the observed properties of bursts with durations short as 10-20% of the spin period can vary strongly depending on the rotational phase at which the burst was emitted. We also show that detectability of a rotational phase dependence depends strongly on the minimum number of bursts observed, and find that existing burst samples may simply be too small to rule out a phase dependence.

Nuclear burning and its dependence on the mass accretion rate are fundamental ingredients for describing the complicated observational phenomenology of neutron stars in binary systems. Motivated by high quality burst rate data emerging from large statistical studies, we report general calculations relating bursting rate to mass accretion rate and neutron star rotation frequency. In this first work we neglect general relativistic effects and accretion topology, though we discuss where their inclusion should play a role. The relations we derive are suitable for different burning regimes and provide a direct link between parameters predicted by theory and what is to be expected in observations. We illustrate this for analytical relations of different unstable burning regimes that operate on the surface of an accreting neutron star. We also use the observed behaviour of burst rate to suggest new constraints on burning parameters. We are able to provide an explanation for the long standing problem of the observed decrease of burst rate with increasing mass accretion that follows naturally from these calculations: when accretion rate crosses a certain threshold, ignition moves away from its initially preferential site and this can cause a net reduction of the burst rate due to the effects of local conditions that set local differences in both burst rate and stabilization criteria. We show under which conditions this can happen even if locally the burst rate keeps increasing with accretion.

Matter in neutron star cores reaches extremely high densities, forming states of matter that cannot be generated in the laboratory. The Equation of State (EOS) of the matter links to macroscopic observables, such as mass M and radius R, via the stellar structure equations. A promising technique for measuring M and R exploits hotspots (burst oscillations) that form on the stellar surface when material accreted from a companion star undergoes a thermonuclear explosion. As the star rotates, the hotspot gives rise to a pulsation, and relativistic effects encode information about M and R into the pulse profile. However the burst oscillation mechanism remains unknown, introducing uncertainty when inferring the EOS. I review the progress that we are making towards cracking this long-standing problem, and establishing burst oscillations as a robust tool for measuring M and R. This is a major goal for future large area X-ray telescopes.

High energy ($> 250$ keV) emission has been detected persisting for several tens of seconds after the initial spike of magnetar giant flares. It has been conjectured that this emission might arise via inverse Compton scattering in a highly extended corona generated by super-Eddington outflows high up in the magnetosphere. In this paper we undertake a detailed examination of this model. We investigate the properties of the required scatterers, and whether the mechanism is consistent with the degree of pulsed emission observed in the tail of the giant flare. We conclude that the mechanism is consistent with current data, although the origin of the scattering population remains an open question. We propose an alternative picture in which the emission is closer to that star and is dominated by synchrotron radiation. The $RHESSI$ observations of the December 2004 flare modestly favor this latter picture. We assess the prospects for the Fermi Gamma-Ray Space Telescope to detect and characterize a similar high energy component in a future giant flare. Such a detection should help to resolve some of the outstanding issues.

We observed the Rapid Burster with Chandra when it was in the 'banana' state that usually precedes the type-II X-ray bursting 'island' state for which the source is particularly known. We employed the High-Energy Transmission Grating Spectrometer in combination with the ACIS-S detector in continuous clocking mode. The observation yielded 20 thermonuclear type-I X-ray bursts emitted from the neutron star surface with recurrence times between 0.9 and 1.2 hr, and an e-folding decay time scale of 1 min. We searched for narrow spectral features in the burst emission that could constrain the composition of the ashes of the nuclear burning and the compactness of the neutron star, but found none. The upper limit on the equivalent width of narrow absorption lines between 2 and 6 keV is between 5 and 20 eV (single trial 3 sigma confidence level) and on those of absorption edges between 150 and 400 eV. The latter numbers are comparable to the levels predicted by Weinberg, Bildsten & Schatz (2006) for Eddington-limited thermonuclear bursts.

Even as our measurements of cosmological parameters improve, the physical nature of the dark sector of the universe largely remains a mystery. Many effects of dark sector models are most prominent at very large scales and will rely on future galaxy surveys to elucidate. In this paper we compare the topological properties of the large scale dark matter distribution in a number of cosmological models using hydrodynamical simulations and the cosmological genus statistic. Genus curves are computed from z = 11 to z = 0 for {\Lambda}CDM, Quintessence and Warm Dark Matter models, over a scale range of 1 to 20 Mpc/h. The curves are analysed in terms of their Hermite spectra to describe the power contained in non-Gaussian deformations to the cosmological density field. We find that the {\Lambda}CDM and {\Lambda}WDM models produce nearly identical genus curves indicating no topological differences in structure formation. The Quintessence model, which differs solely in its expansion history, produces significant differences in the strength and redshift evolution of non-Gaussian modes associated with higher cluster abundances and lower void abundances. These effects are robust to cosmic variance and are characteristically different from those produced by tweaking the parameters of a {\Lambda}CDM model. Given the simplicity and similarity of the models, detecting these discrepancies represents a promising avenue for understanding the effect of non-standard cosmologies on large-scale structure.

We have compiled the most comprehensive burst sample from magnetar 4U 0142+61, comprising 27 bursts from its three burst-active episodes in 2011, 2012 and the latest one in 2015 observed with Swift/BAT and Fermi/GBM. Bursts from 4U 0142+61 morphologically resemble typical short bursts from other magnetars. However, 4U 0142+61 bursts are less energetic compared to the bulk of magnetar bursts. We uncovered an extended tail emission following a burst on 2015 February 28, with a thermal nature, cooling over a time-scale of several minutes. During this tail emission, we also uncovered pulse peak phase aligned X-ray bursts , which could originate from the same underlying mechanism as that of the extended burst tail, or an associated and spatially coincident but different mechanism.

Neutron stars in low mass X-ray binaries exhibit oscillations during thermonuclear bursts, attributed to asymmetric brightness patterns on the burning surfaces. All models that have been proposed to explain the origin of these asymmetries (spreading hotspots, surface waves, and cooling wakes) depend on the accretion rate. By analysis of archival RXTE data of six oscillation sources, we investigate the accretion rate dependence of the amplitude of burst oscillations. This more than doubles the size of the sample analysed previously by Muno et al. (2004), who found indications for a relationship between accretion rate and oscillation amplitudes. We find that burst oscillation signals can be detected at all observed accretion rates. Moreover, oscillations at low accretion rates are found to have relatively small amplitudes ($A_\text{rms}\leq0.10$) while oscillations detected in bursts observed at high accretion rates cover a broad spread in amplitudes ($0.05\leq A_\text{rms}\leq0.20$). In this paper we present the results of our analysis and discuss these in the light of current burst oscillation models. Additionally, we investigate the bursts of two sources without previously detected oscillations. Despite that these sources have been observed at accretion rates where burst oscillations might be expected, we find their behaviour to be not anomalous compared to oscillation sources.

Two energetic hard X-ray bursts have recently triggered the Fermi and Swift space observatories from the rotation powered pulsar, PSR J1119-6127. We have performed in depth spectral and temporal analyses of these two events. Our extensive searches in both observatory data for lower luminosity bursts uncovered 10 additional events from the source. We report here on the timing and energetics of the 12 bursts from PSR J1119-6127 during its burst active phase of 2016 July 26 and 28. We also found a spectral softer X-ray flux enhancement in a post burst episode, which shows evidence of cooling. We discuss here the implications of these results on the nature of this unusual high-field radio pulsar, which firmly place it within the typical magnetar population.

The origin of the observed steep rotation curves of blue compact dwarf galaxies (BCDs) remains largely unexplained by theoretical models of BCD formation. We therefore investigate the rotation curves in BCDs formed from mergers between gas-rich dwarf irregular galaxies based on the results of numerical simulations for BCD formation. The principal results are as follows. The dark matter of merging dwarf irregulars undergoes a central concentration so that the central density can become up to six times higher than those of the initial dwarf irregulars. However, the more compact dark matter halo alone can not reproduce the gradient differences observed between dwarf irregulars and BCDs. We provide further support that the central concentration of gas due to rapid gas-transfer to the central regions of dwarf--dwarf mergers is responsible for the observed difference in rotation curve gradients. The BCDs with central gas concentration formed from merging can thus show steeply rising rotation curves in their central regions. Such gas concentration is also responsible for central starbursts of BCDs and the high central surface brightness and is consistent with previous BCD studies. We discuss the relationship between rotational velocity gradient and surface brightness, the dependence of BCD rotation curves on star formation threshold density, progenitor initial profile, interaction type and merger mass ratio, as well as potential evolutionary links between dwarf irregulars, BCDs and compact dwarf irregulars.

LOFT-P is a concept for a NASA Astrophysics Probe-Class (<$1B) X-ray timing mission, based on the LOFT concept originally proposed to ESAs M3 and M4 calls. LOFT-P requires very large collecting area (>6 m^2, >10x RXTE), high time resolution, good spectral resolution, broad-band spectral coverage (2-30 keV), highly flexible scheduling, and an ability to detect and respond promptly to time-critical targets of opportunity. It addresses science questions such as: What is the equation of state of ultra dense matter? What are the effects of strong gravity on matter spiraling into black holes? It would be optimized for sub-millisecond timing to study phenomena at the natural timescales of neutron star surfaces and black hole event horizons and to measure mass and spin of black holes. These measurements are synergistic to imaging and high-resolution spectroscopy instruments, addressing much smaller distance scales than are possible without very long baseline X-ray interferometry, and using complementary techniques to address the geometry and dynamics of emission regions. A sky monitor (2-50 keV) acts as a trigger for pointed observations, providing high duty cycle, high time resolution monitoring of the X-ray sky with ~20 times the sensitivity of the RXTE All-Sky Monitor, enabling multi-wavelength and multi-messenger studies. A probe-class mission concept would employ lightweight collimator technology and large-area solid-state detectors, technologies which have been recently greatly advanced during the ESA M3 study. Given the large community interested in LOFT (>800 supporters, the scientific productivity of this mission is expected to be very high, similar to or greater than RXTE (~2000 refereed publications). We describe the results of a study, recently completed by the MSFC Advanced Concepts Office, that demonstrates that LOFT-P is feasible within a NASA probe-class mission budget.

eXTP is a science mission designed to study the state of matter under extreme conditions of density, gravity and magnetism. Primary targets include isolated and binary neutron stars, strong magnetic field systems like magnetars, and stellar-mass and supermassive black holes. The mission carries a unique and unprecedented suite of state-of-the-art scientific instruments enabling for the first time ever the simultaneous spectral-timing-polarimetry studies of cosmic sources in the energy range from 0.5-30 keV (and beyond). Key elements of the payload are: the Spectroscopic Focusing Array (SFA) - a set of 11 X-ray optics for a total effective area of about 0.9 m^2 and 0.6 m^2 at 2 keV and 6 keV respectively, equipped with Silicon Drift Detectors offering <180 eV spectral resolution; the Large Area Detector (LAD) - a deployable set of 640 Silicon Drift Detectors, for a total effective area of about 3.4 m^2, between 6 and 10 keV, and spectral resolution <250 eV; the Polarimetry Focusing Array (PFA) - a set of 2 X-ray telescope, for a total effective area of 250 cm^2 at 2 keV, equipped with imaging gas pixel photoelectric polarimeters; the Wide Field Monitor (WFM) - a set of 3 coded mask wide field units, equipped with position-sensitive Silicon Drift Detectors, each covering a 90 degrees x 90 degrees FoV. The eXTP international consortium includes mostly major institutions of the Chinese Academy of Sciences and Universities in China, as well as major institutions in several European countries and the United States. The predecessor of eXTP, the XTP mission concept, has been selected and funded as one of the so-called background missions in the Strategic Priority Space Science Program of the Chinese Academy of Sciences since 2011. The strong European participation has significantly enhanced the scientific capabilities of eXTP. The planned launch date of the mission is earlier than 2025.

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.

Type I X-ray bursts are produced by thermonuclear runaways that develop on accreting neutron stars. Once one location ignites, the flame propagates across the surface of the star. Flame propagation is fundamental in order to understand burst properties like rise time and burst oscillations. Previous work quantified the effects of rotation on the front, showing that the flame propagates as a deflagration and that the front strongly resembles a hurricane. However the effect of magnetic fields was not investigated, despite the fact that magnetic fields strong enough to have an effect on the propagating flame are expected to be present on many bursters. In this paper we show how the coupling between fluid layers introduced by an initially vertical magnetic field plays a decisive role in determining the character of the fronts that are responsible for the Type I bursts. In particular, on a star spinning at 450 Hz (typical among the bursters) we test seed magnetic fields of $10^{7} - 10^{10}$ G and find that for the medium fields the magnetic stresses that develop during the burst can speed up the velocity of the burning front, bringing the simulated burst rise time close to the observed values. By contrast, in a magnetic slow rotator like IGR J17480--2446, spinning at 11 Hz, a seed field $\gtrsim 10^9$ G is required to allow localized ignition and the magnetic field plays an integral role in generating the burst oscillations observed during the bursts.

We report on the analysis of two deep XMM-Newton observations of the magnetar Swift J1834.9-0846 and its surrounding extended emission taken in March 2014 and October 2014, 2.5 and 3.1 years after the source went into outburst. The magnetar is only weakly detected in the first observation with an absorption corrected flux $F_{\rm 0.5-10 keV}\approx4\times10^{-14}$ erg s$^{-1}$ cm$^{-2}$, and a $3\sigma$ upper limit during the second observation of about $3\times10^{-14}$ erg s$^{-1}$ cm$^{-2}$. This flux level is more than 3 orders of magnitude lower than the flux measured at the outburst onset on September 2011. The extended emission, centered at the magnetar position and elongated towards the south-west, is clearly seen in both observations; it is best fit by a highly absorbed power-law (PL), with a hydrogen column density of $N_{\rm H}=8.0\times10^{22}$ cm$^{-2}$ and PL photon index $\Gamma=2.2\pm0.2$. Its flux is constant between the two observations at $F_{\rm 0.5-10 keV}=1.3\times10^{-12}$ erg s$^{-1}$ cm$^{-2}$. We find no statistically significant changes in the spectral shape or the flux of this extended emission over a period of 9 years from 2005 to 2014. These new results strongly support the extended emission nature as a wind nebula and firmly establish Swift J1834.9-0846 as the first magnetar to show a surrounding wind nebula. Further, our results imply that such nebulae are no longer exclusive to rotation-powered pulsars and narrow the gap between these two sub-populations of isolated neutron stars. The size and spectrum of the nebula are compatible with those of pulsar-wind nebulae but its radiative efficiency $\eta_{\rm X}=L_{\rm X}/\dot{E}\approx0.1$ is markedly high, possibly pointing to an additional wind component in Swift J1834.9-0846.

One of the primary science goals of the next generation of hard X-ray timing instruments is to determine the equation of state of the matter at supranuclear densities inside neutron stars, by measuring the radius of neutron stars with different masses to accuracies of a few percent. Three main techniques can be used to achieve this goal. The first involves waveform modelling. The flux we observe from a hotspot on the neutron star surface offset from the rotational pole will be modulated by the star's rotation, giving rise to a pulsation. Information about mass and radius is encoded into the pulse profile via relativistic effects, and tight constraints on mass and radius can be obtained. The second technique involves characterising the spin distribution of accreting neutron stars. The most rapidly rotating stars provide a very clean constraint, since the mass-shedding limit is a function of mass and radius. However the overall spin distribution also provides a guide to the torque mechanisms in operation and the moment of inertia, both of which can depend sensitively on dense matter physics. The third technique is to search for quasi-periodic oscillations in X-ray flux associated with global seismic vibrations of magnetars (the most highly magnetized neutron stars), triggered by magnetic explosions. The vibrational frequencies depend on stellar parameters including the dense matter equation of state. We illustrate how these complementary X-ray timing techniques can be used to constrain the dense matter equation of state, and discuss the results that might be expected from a 10m$^2$ instrument. We also discuss how the results from such a facility would compare to other astronomical investigations of neutron star properties. [Modified for arXiv]

Giant $\gamma$-ray flares comprise the most extreme radiation events observed from magnetars. Developing on (sub)millisecond timescales and generating vast amounts of energy within a fraction of a second, the initial phase of these extraordinary bursts present a significant challenge for candidate trigger mechanisms. Here we assess and critically analyse the linear growth of the relativistic tearing instability in a globally twisted magnetosphere as the trigger mechanism for giant $\gamma$-ray flares. Our main constraints are given by the observed emission timescales, the energy output of the giant flare spike, and inferred dipolar magnetic field strengths. We find that the minimum growth time of the linear mode is comparable to the $e$-folding rise time, i.e. $\sim10^{-1}$ ms. With this result we constrain basic geometric parameters of the current sheet. We also discuss the validity of the presumption that the $e$-folding emission timescale may be equated with the growth time of an MHD instability.

Magnetars are the strongest magnets in the present universe and the combination of extreme magnetic field, gravity and density makes them unique laboratories to probe current physical theories (from quantum electrodynamics to general relativity) in the strong field limit. Magnetars are observed as peculiar, burst--active X-ray pulsars, the Anomalous X-ray Pulsars (AXPs) and the Soft Gamma Repeaters (SGRs); the latter emitted also three "giant flares," extremely powerful events during which luminosities can reach up to 10^47 erg/s for about one second. The last five years have witnessed an explosion in magnetar research which has led, among other things, to the discovery of transient, or "outbursting," and "low-field" magnetars. Substantial progress has been made also on the theoretical side. Quite detailed models for explaining the magnetars' persistent X-ray emission, the properties of the bursts, the flux evolution in transient sources have been developed and confronted with observations. New insight on neutron star asteroseismology has been gained through improved models of magnetar oscillations. The long-debated issue of magnetic field decay in neutron stars has been addressed, and its importance recognized in relation to the evolution of magnetars and to the links among magnetars and other families of isolated neutron stars. The aim of this paper is to present a comprehensive overview in which the observational results are discussed in the light of the most up-to-date theoretical models and their implications. This addresses not only the particular case of magnetar sources, but the more fundamental issue of how physics in strong magnetic fields can be constrained by the observations of these unique sources.

Neutron stars are a prime laboratory for testing physical processes under conditions of strong gravity, high density, and extreme magnetic fields. Among the zoo of neutron star phenomena, magnetars stand out for their bursting behaviour, ranging from extremely bright, rare giant flares to numerous, less energetic recurrent bursts. The exact trigger and emission mechanisms for these bursts are not known; favoured models involve either a crust fracture and subsequent energy release into the magnetosphere, or explosive reconnection of magnetic field lines. In the absence of a predictive model, understanding the physical processes responsible for magnetar burst variability is difficult. Here, we develop an empirical model that decomposes magnetar bursts into a superposition of small spike-like features with a simple functional form, where the number of model components is itself part of the inference problem. The cascades of spikes that we model might be formed by avalanches of reconnection, or crust rupture aftershocks. Using Markov Chain Monte Carlo (MCMC) sampling augmented with reversible jumps between models with different numbers of parameters, we characterise the posterior distributions of the model parameters and the number of components per burst. We relate these model parameters to physical quantities in the system, and show for the first time that the variability within a burst does not conform to predictions from ideas of self-organised criticality. We also examine how well the properties of the spikes fit the predictions of simplified cascade models for the different trigger mechanisms.

Since launch in 2008, the Fermi Gamma-ray Burst Monitor (GBM) has detected many hundreds of bursts from magnetar sources. While the vast majority of these bursts have been attributed to several known magnetars, there is also a small sample of magnetar-like bursts of unknown origin. Here we present the Fermi/GBM magnetar catalog, giving the results of the temporal and spectral analyses of 440 magnetar bursts with high temporal and spectral resolution. This catalog covers the first five years of GBM magnetar observations, from July 2008 to June 2013. We provide durations, spectral parameters for various models, fluences and peak fluxes for all the bursts, as well as a detailed temporal analysis for SGR J1550-5418 bursts. Finally, we suggest that some of the bursts of unknown origin are associated with the newly discovered magnetar 3XMM J185246.6+0033.7.

Crustquake events may be connected with both rapid spin-up `glitches' within the regular slowdown of neutron stars, and high-energy magnetar flares. We argue that magnetic field decay builds up stresses in a neutron star's crust, as the elastic shear force resists the Lorentz force's desire to rearrange the global magnetic-field equilibrium. We derive a criterion for crust-breaking induced by a changing magnetic-field configuration, and use this to investigate strain patterns in a neutron star's crust for a variety of different magnetic-field models. Universally, we find that the crust is most liable to break if the magnetic field has a strong toroidal component, in which case the epicentre of the crustquake is around the equator. We calculate the energy released in a crustquake as a function of the fracture depth, finding that it is independent of field strength. Crust-breaking is, however, associated with a characteristic local field strength of $2.4\times 10^{14}$ G for a breaking strain of $0.001$, or $2.4\times 10^{15}$ G at a breaking strain of $0.1$. We find that even the most luminous magnetar giant flare could have been powered by crustal energy release alone.

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 $2\Omega\cos\theta$ 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.

This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of very faint X-ray binaries, orbital period distribution of black hole X-ray binaries and neutron star spin up. For a summary, we refer to the paper.

This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of thermonuclear X-ray bursts on accreting neutron stars. For a summary, we refer to the paper.

With an average density higher than the nuclear density, neutron stars (NS) provide a unique test-ground for nuclear physics, quantum chromodynamics (QCD), and nuclear superfluidity. Determination of the fundamental interactions that govern matter under such extreme conditions is one of the major unsolved problems of modern physics, and -- since it is impossible to replicate these conditions on Earth -- a major scientific motivation for SKA. The most stringent observational constraints come from measurements of NS bulk properties: each model for the microscopic behaviour of matter predicts a specific density-pressure relation (its `Equation of state', EOS). This generates a unique mass-radius relation which predicts a characteristic radius for a large range of masses and a maximum mass above which NS collapse to black holes. It also uniquely predicts other bulk quantities, like maximum spin frequency and moment of inertia. The SKA, in Phase 1 and particularly in Phase 2 will, thanks to the exquisite timing precision enabled by its raw sensitivity, and surveys that dramatically increase the number of sources: 1) Provide many more precise NS mass measurements (high mass NS measurements are particularly important for ruling out EOS models); 2) Allow the measurement of the NS moment of inertia in highly relativistic binaries such as the Double Pulsar; 3) Greatly increase the number of fast-spinning NS, with the potential discovery of spin frequencies above those allowed by some EOS models; 4) Improve our knowledge of new classes of binary pulsars such as black widows and redbacks (which may be massive as a class) through sensitive broad-band radio observations; and 5) Improve our understanding of dense matter superfluidity and the state of matter in the interior through the study of rotational glitches, provided that an ad-hoc campaign is developed.

Providing a link between magnetars and radio pulsars, high magnetic field neutron stars are ideal targets to investigate how bursting/magnetospheric activity and braking torque variations are connected to rotational glitches. The last spin-up glitch of the highly magnetised pulsar J1119$-$6127 back in 2007 was the first glitch in a rotationally powered radio pulsar to be accompanied by radiative changes. Moreover, it was followed by an uncommon glitch relaxation that resulted in a smaller spin-down rate relative to the prediction of the pre-glitch timing model. Here, we present 4 years of new radio timing observations and analyse the total of 16 years of timing data for this source. The new data uncover an ongoing evolution of the spin-down rate, thereby allowing us to exclude permanent changes in the external or internal torque as a stand-alone cause of the peculiar features of the glitch recovery. Furthermore, no additional variations of the radio pulse profile are detected, strengthening the association of the previously observed transient emission features with the glitching activity. A self-consistent measurement of the braking index yields a value $n\simeq2.7$, indicating a trajectory in the $P-\dot{P}$ plane inclined towards the magnetars. Such a potential evolutionary link might be strengthened by a, possibly permanent, reduction of $\sim15\%$ in $n$ at the epoch of the 2007 glitch.

Quasi-periodic oscillations (QPOs) observed in the giant flares of magnetars are of particular interest due to their potential to open up a window into the neutron star interior via neutron star asteroseismology. However, only three giant flares have been observed. We therefore make use of the much larger data set of shorter, less energetic recurrent bursts. Here, we report on a search for QPOs in a large data set of bursts from the two most burst-active magnetars, SGR 1806-20 and SGR 1900+14, observed with the Rossi X-ray Timing Explorer (RXTE). We find a single detection in an averaged periodogram comprising 30 bursts from SGR 1806-20, with a frequency of 57 Hz and a width of 5 Hz, remarkably similar to a giant flare QPO observed from SGR 1900+14. This QPO fits naturally within the framework of global magneto-elastic torsional oscillations employed to explain the giant flare QPOs. Additionally, we uncover a limit on the applicability of Fourier analysis for light curves with low background count rates and strong variability on short timescales. In this regime, standard Fourier methodology and more sophisticated Fourier analyses fail in equal parts by yielding an unacceptably large number of false positive detections. This problem is not straightforward to solve in the Fourier domain. Instead, we show how simulations of light curves can offer a viable solution for QPO searches in these light curves.

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

Quasi-periodic oscillations (QPOs) detected in the 2004 giant flare from SGR 1806-20 are often interpreted as global magneto-elastic oscillations of the neutron star. There is, however, a large discrepancy between theoretical models, which predict that the highest frequency oscillations should die out rapidly, and the observations, which suggested that the highest-frequency signals persisted for ~100s in X-ray data from two different spacecraft. This discrepancy is particularly important for the high-frequency QPO at ~625 Hz. However, previous analyses did not systematically test whether the signal could also be there in much shorter data segments, more consistent with the theoretical predictions. Here, we test for the presence of the high-frequency QPO at 625 Hz in data from both the Rossi X-ray Timing Explorer (RXTE) and the Ramaty High Energy Solar Spectroscopic Imager (RHESSI) systematically both in individual rotational cycles of the neutron star, as well as averaged over multiple successive rotational cycles at the same phase. We find that the QPO in the RXTE data is consistent with being only present in a single cycle, for a short duration of ~0.5s, whereas the RHESSI data are as consistent with a short-lived signal that appears and disappears as with a long-lived QPO. Taken together, this data provides evidence for strong magnetic interaction between the crust and the core.

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.

Glitches are sudden spin-up events that punctuate the steady spin down of pulsars and are thought to be due to the presence of a superfluid component within neutron stars. The precise glitch mechanism and its trigger, however, remain unknown. The size of glitches is a key diagnostic for models of the underlying physics. While the largest glitches have long been taken into account by theoretical models, it has always been assumed that the minimum size lay below the detectability limit of the measurements. In this paper we define general glitch detectability limits and use them on 29 years of daily observations of the Crab pulsar, carried out at Jodrell Bank Observatory. We find that all glitches lie well above the detectability limits and by using an automated method to search for small events we are able to uncover the full glitch size distribution, with no biases. Contrary to the prediction of most models, the distribution presents a rapid decrease of the number of glitches below ~0.05 $\mu$Hz. This substantial minimum size indicates that a glitch must involve the motion of at least several billion superfluid vortices and provides an extra observable which can greatly help the identification of the trigger mechanism. Our study also shows that glitches are clearly separated from all the other rotation irregularities. This supports the idea that the origin of glitches is different to that of timing noise, which comprises the unmodelled random fluctuations in the rotation rates of pulsars.

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.

LOFT, the Large Observatory For X-ray Timing, is one of five ESA M3 candidate missions. It will address the Cosmic Vision theme: "Matter under Extreme Conditions". By coupling for the first time a huge collecting area for the detection of X-ray photons with CCD-quality spectral resolution (15 times bigger in area than any previously flown X-ray instrument and >100 times bigger for spectroscopy than any similar-resolution instrument), the instruments on-board LOFT have been designed to (i) determine the properties of ultradense matter by reconstructing its Equation of State through neutron star mass and radius measurements of unprecedented accuracy; (ii) measure General Relativity effects in the strong field regime in the stationary spacetimes of neutron stars and black holes of all masses down to a few gravitational radii. Besides the above two themes, LOFT's observations will be devoted to "observatory science", providing new insights in a number of research fields in high energy astrophysics (e.g. Gamma-ray Bursts). The assessment study phase of LOFT, which ended in September 2013, demonstrated that the mission is low risk and the required Technology Readiness Level can be easily reached in time for a launch by the end of 2022.

Accreting Millisecond X-Ray Pulsars (AMXPs) are astrophysical laboratories without parallel in the study of extreme physics. In this chapter we review the past fifteen years of discoveries in the field. We summarize the observations of the fifteen known AMXPs, with a particular emphasis on the multi-wavelength observations that have been carried out since the discovery of the first AMXP in 1998. We review accretion torque theory, the pulse formation process, and how AMXP observations have changed our view on the interaction of plasma and magnetic fields in strong gravity. We also explain how the AMXPs have deepened our understanding of the thermonuclear burst process, in particular the phenomenon of burst oscillations. We conclude with a discussion of the open problems that remain to be addressed in the future.

Long suspected to be due to unstable accretion events, the type II bursts exhibited by the Rapid Burster (RB, or MXB 1730-335) still lack an exhaustive explanation. Apart from type II bursts, the transient RB also shows the better-understood thermonuclear shell flashes known as type I bursts. In search of links between these two phenomena, we carried out a comprehensive analysis of all $\textit{Rossi X-ray Timing Explorer}$ observations of the RB and found six atypical type I bursts, featuring a double-peaked profile that is not due to photospheric radius expansion. The bursts appear in a phase of the outburst decay close to the onset of the type II bursts, when the source also switches from the high/soft to the low/hard state. We also report the discovery of a simultaneous low-frequency quasi-periodic oscillation present in the persistent emission as well as in the burst decaying emission. We discuss several scenarios to understand the nature of the peculiar bursts and of the accompanying oscillation, as well as their possible relation with each other and with the type II burst phenomenon. We favour a model consisting of two accretion channels, one polar and one equatorial, in a source viewed at low inclination. We also, for the first time, clearly identify the atoll nature of the RB.

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.

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.

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.

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.

We perform time-resolved spectroscopy of all the type I bursts from the Rapid Burster (MXB 1730-335) detected with the Rossi X-ray Timing Explorer. Type I bursts are detected at high accretion rates, up to \sim 45% of the Eddington luminosity. We find evidence that bursts lacking the canonical cooling in their time-resolved spectra are, none the less, thermonuclear in nature. The type I bursting rate keeps increasing with the persistent luminosity, well above the threshold at which it is known to abruptly drop in other bursting low-mass X-ray binaries. The only other known source in which the bursting rate keeps increasing over such a large range of mass accretion rates is the 11 Hz pulsar IGR J17480$-$2446. This may indicate a similarly slow spin for the neutron star in the Rapid Burster.

The LOFT mission concept is one of four candidates selected by ESA for the M3 launch opportunity as Medium Size missions of the Cosmic Vision programme. The launch window is currently planned for between 2022 and 2024. LOFT is designed to exploit the diagnostics of rapid X-ray flux and spectral variability that directly probe the motion of matter down to distances very close to black holes and neutron stars, as well as the physical state of ultra-dense matter. These primary science goals will be addressed by a payload composed of a Large Area Detector (LAD) and a Wide Field Monitor (WFM). The LAD is a collimated (<1 degree field of view) experiment operating in the energy range 2-50 keV, with a 10 m^2 peak effective area and an energy resolution of 260 eV at 6 keV. The WFM will operate in the same energy range as the LAD, enabling simultaneous monitoring of a few-steradian wide field of view, with an angular resolution of <5 arcmin. The LAD and WFM experiments will allow us to investigate variability from submillisecond QPO's to year-long transient outbursts. In this paper we report the current status of the project.

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.

Neutron stars enable us to study both the highest densities and the highest magnetic fields in the known Universe. In this article I review what can be learned about such fundamental physics using magnetar bursts. Both the instability mechanisms that trigger the bursts, and the subsequent dynamical and radiative response of the star, can be used to explore stellar and magnetospheric structure and composition.

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.

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.

Burst oscillations, a phenomenon observed in a significant fraction of Type I (thermonuclear) X-ray bursts, involve the development of highly asymmetric brightness patches in the burning surface layers of accreting neutron stars. Intrinsically interesting as nuclear phenomena, they are also important as probes of dense matter physics and the strong gravity, high magnetic field environment of the neutron star surface. Burst oscillation frequency is also used to measure stellar spin, and doubles the sample of rapidly rotating (above 10 Hz) accreting neutron stars with known spins. Although the mechanism remains mysterious, burst oscillation models must take into account thermonuclear flame spread, nuclear processes, rapid rotation, and the dynamical role of the magnetic field. This review provides a comprehensive summary of the observational properties of burst oscillations, an assessment of the status of the theoretical models that are being developed to explain them, and an overview of how they can be used to constrain neutron star properties such as spin, mass and radius.

We have performed detailed temporal and time-integrated spectral analysis of 286 bursts from SGR J1550-5418 detected with the Fermi Gamma-ray Burst Monitor (GBM) in January 2009, resulting in the largest uniform sample of temporal and spectral properties of SGR J1550-5418 bursts. We have used the combination of broadband and high time-resolution data provided with GBM to perform statistical studies for the source properties. We determine the durations, emission times, duty cycles and rise times for all bursts, and find that they are typical of SGR bursts. We explore various models in our spectral analysis, and conclude that the spectra of SGR J1550-5418 bursts in the 8-200 keV band are equally well described by optically thin thermal bremsstrahlung (OTTB), a power law with an exponential cutoff (Comptonized model), and two black-body functions (BB+BB). In the spectral fits with the Comptonized model we find a mean power-law index of -0.92, close to the OTTB index of -1. We show that there is an anti-correlation between the Comptonized Epeak and the burst fluence and average flux. For the BB+BB fits we find that the fluences and emission areas of the two blackbody functions are correlated. The low-temperature BB has an emission area comparable to the neutron star surface area, independent of the temperature, while the high-temperature blackbody has a much smaller area and shows an anti-correlation between emission area and temperature. We compare the properties of these bursts with bursts observed from other SGR sources during extreme activations, and discuss the implications of our results in the context of magnetar burst models.

The most strongly magnetized neutron stars, the magnetars, have spectacular outbursts of gamma-ray flares powered by decay of the magnetic field. The rapidly changing field is strong enough that it should be able to stress and rupture the star's crust, with potentially interesting seismic consequences. This was confirmed in dramatic fashion a few years ago with the discovery of long-lived seismic vibrations excited by rare giant flares. The starquakes associated with the giant flares are, it seems, so catastrophic that they leave the whole star ringing. These discoveries have opened up the possibility of using asteroseismology to study the extreme conditions of the neutron star interior - the dense matter equation of state in crust and core, and the core magnetic field. This chapter reviews the main observational properties of the oscillations seen in the giant flares, and the state of the art in terms of theoretical models of global seismic vibrations in magnetars. I discuss efforts to extend oscillation searches to smaller magnetar bursts (which are more frequent but less energetic), and the role that crust rupturing might play in triggering magnetar flares and exciting global seismic oscillations.

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.

The recently-discovered accreting X-ray pulsar IGR J17480--2446 spins at a frequency of ~11 Hz. We show that Type I X-ray bursts from this source display oscillations at the same frequency as the stellar spin. IGR J17480--2446 is the first secure case of a slowly rotating neutron star which shows Type I burst oscillations, all other sources featuring such oscillations spin at hundreds of Hertz. This means that we can test burst oscillation models in a completely different regime. We explore the origin of Type I burst oscillations in IGR J17480--2446 and conclude that they are not caused by global modes in the neutron star ocean. We also show that the Coriolis force is not able to confine an oscillation-producing hot-spot on the stellar surface. The most likely scenario is that the burst oscillations are produced by a hot-spot confined by hydromagnetic stresses.

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).

We report on the discovery and the timing analysis of the first eclipsing accretion-powered millisecond X-ray pulsar (AMXP): SWIFT J1749.4-2807. The neutron star rotates at a frequency of ~517.9 Hz and is in a binary system with an orbital period of 8.8 hrs and a projected semi-major axis of ~1.90 lt-s. Assuming a neutron star between 0.8 and 2.2 M_o and using the mass function of the system and the eclipse half-angle, we constrain the mass of the companion and the inclination of the system to be in the ~0.46-0.81 M_o and $\sim74.4^o-77.3^o range, respectively. To date, this is the tightest constraint on the orbital inclination of any AMXP. As in other AMXPs, the pulse profile shows harmonic content up to the 3rd overtone. However, this is the first AMXP to show a 1st overtone with rms amplitudes between ~6% and ~23%, which is the strongest ever seen, and which can be more than two times stronger than the fundamental. The fact that SWIFT J1749.4-2807 is an eclipsing system which shows uncommonly strong harmonic content suggests that it might be the best source to date to set constraints on neutron star properties including compactness and geometry.

We perform phase-resolved spectroscopy of the accreting millisecond pulsar, SAX J1808.4-3658, during the slow-decay phase of the 2002 outburst. Simple phenomenological fits to RXTE PCA data reveal a pulsation in the iron line at the spin frequency of the neutron star. However, fitting more complex spectral models reveals a degeneracy between iron-line pulsations and changes in the underlying hotspot blackbody temperature with phase. By comparing with the variations in reflection continuum, which are much weaker than the iron line variations, we infer that the iron-line is not pulsed. The observed spectral variations can be explained by variations in blackbody temperature associated with rotational Doppler shifts at the neutron star surface. By allowing blackbody temperature to vary in this way, we also find a larger phase-shift between the pulsations in the Comptonised and blackbody components than has been seen in previous work. The phase-shift between the pulsation in the blackbody temperature and normalisation is consistent with a simple model where the Doppler shift is maximised at the limb of the neutron star, ~90 degrees prior to maximisation of the hot-spot projected area.

We report the discovery of burst oscillations at the spin frequency in ten thermonuclear bursts from the accreting millisecond X-ray pulsar (AMXP) IGR J17511-3057. The burst oscillation properties are, like those from the persistent AMXPs SAX J1808.4-3658 and XTE J1814-338, anomalous compared to burst oscillations from intermittent pulsars or non-pulsing LMXBs. Like SAX J1808.4-3658 they show frequency drifts in the rising phase rather than the tail. There is also evidence for harmonic content. Where IGR J17511-3057 is unusual compared to the other two persistent pulsars is that oscillations are not detected throughout all bursts. As accretion rate drops the bursts get brighter and their rise/decay time scales become shorter, while the oscillation amplitude falls below the detection threshold: first in the burst peak and then also in the rise. None of the bursts from IGR J17511-3057 show evidence for photospheric radius expansion (which might be expected to suppress oscillation amplitude) which allow us to set an upper limit to the distance of 6.9 kpc. We discuss the implications of our results for models of the burst oscillation mechanism.

We report the detection of 15 X-ray bursts with RXTE and Swift observations of the peculiar X-ray binary Circinus X-1 during its May 2010 X-ray re-brightening. These are the first X-ray bursts observed from the source after the initial discovery by Tennant and collaborators, twenty-five years ago. By studying their spectral evolution, we firmly identify nine of the bursts as type I (thermonuclear) X-ray bursts. We obtain an arcsecond location of the bursts that confirms once and for all the identification of Cir X-1 as a type I X-ray burst source, and therefore as a low magnetic field accreting neutron star. The first five bursts observed by RXTE are weak and show approximately symmetric light curves, without detectable signs of cooling along the burst decay. We discuss their possible nature. Finally, we explore a scenario to explain why Cir X-1 shows thermonuclear bursts now but not in the past, when it was extensively observed and accreting at a similar rate.

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.

On 2009 June 5, the Gamma-ray Burst Monitor (GBM) onboard the Fermi Gamma-ray Space Telescope triggered on two short, and relatively dim bursts with spectral properties similar to Soft Gamma Repeater (SGR) bursts. Independent localizations of the bursts by triangulation with the Konus-RF and with the Swift satellite, confirmed their origin from the same, previously unknown, source. The subsequent discovery of X-ray pulsations with the Rossi X-ray Timing Explorer (RXTE), confirmed the magnetar nature of the new source, SGR J0418+5729. We describe here the Fermi/GBM observations, the discovery and the localization of this new SGR, and our infrared and Chandra X-ray observations. We also present a detailed temporal and spectral study of the two GBM bursts. SGR J0418+5729 is the second source discovered in the same region of the sky in the last year, the other one being SGR J0501+4516. Both sources lie in the direction of the galactic anti-center and presumably at the nearby distance of ~2 kpc (assuming they reside in the Perseus arm of our galaxy). The near-threshold GBM detection of bursts from SGR J0418+5729 suggests that there may be more such dim SGRs throughout our galaxy, possibly exceeding the population of bright SGRs. Finally, using sample statistics, we conclude that the implications of the new SGR discovery on the number of observable active magnetars in our galaxy at any given time is <10, in agreement with our earlier estimates.

SGR J1550-5418 (previously known as AXP 1E 1547.0-5408 or PSR J1550-5418) went into three active bursting episodes in 2008 October and in 2009 January and March, emitting hundreds of typical Soft Gamma Repeater (SGR) bursts in soft gamma-rays. The second episode was especially intense, and our untriggered burst search on Fermi/GBM data (8-1000 keV) revealed ~450 bursts emitted over 24 hours during the peak of this activity. Using the GBM data, we identified a ~150-s-long enhanced persistent emission during 2009 January 22 that exhibited intriguing timing and spectral properties: (i) clear pulsations up to ~110 keV at the spin period of the neutron star (P ~ 2.07 s, the fastest of all magnetars), (ii) an additional (to a power-law) blackbody component required for the enhanced emission spectra with kT ~ 17 keV, (iii) pulsed fraction that is strongly energy dependent and highest in the 50-74 keV energy band. A total isotropic-equivalent energy emitted during this enhanced emission is estimated to be 2.9 x 10^{40} (D/5 kpc)^2 erg. The estimated area of the blackbody emitting region of ~0.046(D/5 kpc)^2 km^2 (roughly a few x 10^{-5} of the neutron star area) is the smallest "hot spot" ever measured for a magnetar and most likely corresponds to the size of magnetically-confined plasma near the neutron star surface.

We present a simultaneous periodic and aperiodic timing study of the accreting millisecond X-ray pulsar SAX J1808.4-3658. We analyze five outbursts of the source and for the first time provide a full and systematic investigation of the enigmatic phenomenon of the 1 Hz flares observed during the final stages of some of the outbursts. We show that links between pulsations and 1 Hz flares might exist, and suggest they are related with hydrodynamic disk instabilities that are triggered close to the disk-magnetosphere boundary layer when the system is entering the propeller regime.

We show that the fundamental seismic shear mode, observed as a quasi-periodic oscillation in giant flares emitted by highly-magnetized neutron stars, is particularly sensitive to the nuclear physics of the crust. The identification of an oscillation at ~ 30 Hz as the fundamental crustal shear mode requires a nuclear symmetry energy that depends very weakly on density near saturation. If the nuclear symmetry energy varies more strongly with density, then lower frequency oscillations, previously identified as torsional Alfven modes of the fluid core, could instead be associated with the crust. If this is the case, then future observations of giant flares should detect oscillations at around 18 Hz. An accurate measurement of the neutron skin thickness of lead will also constrain the frequencies predicted by the model.

We report the discovery of burst oscillations from the intermittent accretion-powered millisecond pulsar (AMP) HETE J1900.1-2455, with a frequency approximately 1 Hz below the known spin frequency. The burst oscillation properties are far more similar to those of the non-AMPs and Aql X-1 (an intermittent AMP with a far lower duty cycle), than those of the AMPs SAX J1808.4-3658 and XTE J1814-338. We discuss the implications for models of the burst oscillation and intermittency mechanisms.

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 gravitational waves emitted by neutron stars carry unique information about their structure and composition. Direct detection of these gravitational waves, however, is a formidable technical challenge. In a recent study we quantified the hurdles facing searches for gravitational waves from the known accreting neutron stars, given the level of uncertainty that exists regarding spin and orbital parameters. In this paper we reflect on our conclusions, and issue an open challenge to the theoretical community to consider how searches should be designed to yield the most astrophysically interesting upper limits. With this in mind we examine some more optimistic emission scenarios involving spin-down, and show that there are technically feasible searches, particularly for the accreting millisecond pulsars, that might place meaningful constraints on torque mechanisms. We finish with a brief discussion of prospects for indirect detection.

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.

X-ray timing of the accretion-powered pulsations during the 2003 outburst of the accreting millisecond pulsar XTE J1814-338 has revealed variation in the pulse time of arrival residuals. These can be interpreted in several ways, including spin-down and wandering of the fuel impact point around the magnetic pole. In this Letter we show that the burst oscillations of this source are coherent with the persistent pulsations, to the level where they track all of the observed fluctuations. Only one burst, which occurs at the lowest accretion rates, shows a significant phase offset. We discuss what might lead to such rigid phase-locking between the modulations in the accretion and thermonuclear burst emission, and consider the implications for spin variation and the burst oscillation mechanism. Wandering of the fuel impact hot spot around a fixed magnetic pole seems the most likely cause for the accretion-powered pulse phase variations. This means that the burst asymmetry is coupled to the hot spot, not the magnetic pole. If premature ignition at this point (due to higher local temperatures) triggers a burning front that stalls before spreading over the entire surface, the resulting localized nuclear hot spot may explain the unusual burst and burst oscillation properties of this source.

The key contribution of the discovery of nuclear-powered pulsations from the accretion-powered millisecond pulsars (AMPs) has been the establishment of burst oscillation frequency as a reliable proxy for stellar spin rate. This has doubled the sample of rapidly-rotating accreting neutron stars and revealed the unexpected absence of any stars rotating near the break-up limit. The resulting `braking problem' is now a major concern for theorists, particularly given the possible role of gravitational wave emission in limiting spin. This, however, is not the only area where burst oscillations from the AMPs are having an impact. Burst oscillation timing is developing into a promising technique for verifying the level of spin variability in the AMPs (a topic of considerable debate). These sources also provide unique input to our efforts to understand the still-elusive burst oscillation mechanism. This is because they are the only stars where we can reliably gauge the role of uneven fuel deposition and, of course, the magnetic field.

Detection of gravitational waves from accreting neutron stars (NSs) in our galaxy, due to ellipticity or internal oscillation, would be a breakthrough in our understanding of compact objects and explain the absence of NSs rotating near the break-up limit. Direct detection, however, poses a formidable challenge. Using the current data available on the properties of the accreting NSs in Low Mass X-Ray Binaries (LMXBs), we quantify the detectability for the known accreting NSs, considering various emission scenarios and taking into account the negative impact of parameter uncertainty on the data analysis process. Only a few of the persistently bright NSs accreting at rates near the Eddington limit are detectable by Advanced LIGO if they are emitting gravitational waves at a rate matching the torque from accretion. A larger fraction of the known population is detectable if the spin and orbital parameters are known in advance, especially with the narrow-band Advanced LIGO. We identify the most promising targets, and list specific actions that would lead to significant improvements in detection probability. These include astronomical observations (especially for unknown orbital periods), improvements in data analysis algorithms and capabilities, and further detector development.

The shape of the lightcurve during the rising phase of Type I X-ray bursts is determined by many factors including the ignition latitude, the accretion rate, and the rotation rate of the star. We develop a phenomenological model of the burst rise process and show that simple measures of the burst morphology can be robust diagnostics of ignition latitude and burning regime. We apply our results to the large sample of bursts from the Low Mass X-ray Binary 4U 1636-536, and find evidence for off-equatorial ignition for many of the bursts. We argue that such behaviour may be associated with the transition from hydrogen to helium ignition at accretion rates a few percent of Eddington. We show that this model can also explain variations in the detectability of burst oscillations, and discuss the implications for other burst sources.

Compact relativistic stars allow us to study the nature of matter under extreme conditions, probing regions of parameter space that are otherwise inaccessible. Nuclear theory in this regime is not well constrained: one key issue is whether neutron stars are in fact composed primarily of strange quark matter. Distinguishing the two possibilities, however, has been difficult. The recent detection of seismic vibrations in the aftermath of giant flares from two magnetars (highly magnetized compact stars) is a major breakthrough. The oscillations excited seem likely to involve the stellar crust, the properties of which differ dramatically for strange stars. We show that the resulting mode frequencies cannot be reconciled with the observations for reasonable magnetar parameters. Ruling out strange star models would place a strong constraint on models of dense quark matter.

Most Type I X-ray bursts from accreting neutron stars have a lightcurve with a single peak, but there is a rare population of faint bursts that are double or even triple peaked. Suggested mechanisms include polar ignition with equatorial stalling, or multi-step energy release; the latter being caused by hydrodynamic instabilities or waiting points in the nuclear reaction sequence. We present an analysis of the accretion rate dependence of the multi-peak bursts, and discuss the consequences for the various models. The observations pose particular challenges for the polar ignition mechanism given current models of ignition latitude dependence.

The high frequency oscillations discovered in the tails of giant flares from two magnetars are thought to be the first direct detections of seismic vibrations from neutron stars. The possibility of starquakes associated with the giant flares triggering global vibrations opens up the prospect of using seismology to study the interior structure and composition of neutron stars. This is a major breakthrough in the study of the nature of matter under conditions of extreme pressure. In this paper we provide an up to date summary of the observations and the theoretical framework, including a brief discussion of gravitational wave searches for the QPOs. We summarize the status of alternative non-seismic mechanisms, and give a critique of a recent paper by Levin that argued against seismic vibrations as a viable mechanism. We conclude with an overview of current results using the seismological technique that constrain parameters such as the equation of state and crust structure.

The nature of the asymmetry that gives rise to Type I X-ray burst oscillations on accreting neutron stars remains a matter of debate. Of particular interest is whether the burst oscillation mechanism differs between the bursting millisecond pulsars and the non-pulsing systems. One means to diagnose this is to study the energy dependence of the burst oscillations: here we present an analysis of oscillations from 28 bursts observed during the 2003 outburst of the accreting millisecond pulsar XTE J1814-338. We find that the fractional amplitude of the burst oscillations falls with energy, in contrast to the behaviour found by Muno et al. (2003) in the burst oscillations from a set of non-pulsing systems. The drop with energy mirrors that seen in the accretion-powered pulsations; in this respect XTE J1814-338 behaves like the other accreting millisecond pulsars. The burst oscillations show no evidence for either hard or soft lags, in contrast to the persistent pulsations, which show soft lags of up to 50 $\mu$s. The fall in amplitude with energy is inconsistent with current surface mode and simple hot spot models of burst oscillations. We discuss improvements to the models and uncertainties in the physics that might resolve these issues.

The recent discovery of high frequency oscillations during giant flares from the Soft Gamma Repeaters SGR 1806-20 and SGR 1900+14 may be the first direct detection of vibrations in a neutron star crust. If this interpretation is correct it offers a novel means of testing the neutron star equation of state, crustal breaking strain, and magnetic field configuration. We review the observational data on the magnetar oscillations, including new timing analysis of the SGR 1806-20 giant flare using data from the Ramaty High Energy Solar Spectroscopic Imager (RHESSI) and the Rossi X-ray Timing Explorer (RXTE). We discuss the implications for the study of neutron star structure and crust thickness, and outline areas for future investigation.

We report an analysis of the archival Rossi X-ray Timing Explorer (RXTE) data from the December 2004 hyperflare from SGR 1806-20. In addition to the 90 Hz QPO first discovered by Israel et al., we report the detection of higher frequency oscillations at 150, 625, and 1,840 Hz. In addition, we also find indications of oscillations at 720, and 2,384 Hz, but with lower significances. The 150 Hz QPO has a width (FWHM) of about 17 Hz, an average amplitude (rms) of 6.8 %, and is detected in average power spectra centered on the rotational phase of the strongest peak in the pulse profile. This is approximately half a cycle from the phase at which the 90 Hz QPO is detected. The 625 Hz oscillation was detected in an average power spectrum from nine successive cycles beginning approximately 180 s after the initial hard spike. It has a width (FWHM) of 2 Hz and an average amplitude (rms) during this interval of 8.5 %. We find a strong detection of the 625 Hz oscillation in a pair of successive rotation cycles begining about 230 s after the start of the flare. In these cycles we also detect the 1,840 Hz QPO. When the 625 Hz QPO is detected we also confirm the simultaneous presence of 30 and 92 Hz QPOs. The centroid frequency of the 625 Hz QPO detected with RXTE is within 1 Hz of the 626 Hz oscillation recently found in RHESSI data from this hyperflare by Watts & Strohmayer. We argue that these new findings provide further evidence for a connection of these oscillations with global oscillation modes of neutron stars, in particular, the high frequency signals may represent toroidal modes with at least one radial node in the crust. We discuss their implications in the context of this model, in particular for the depth of neutron star crusts

The recent discovery of high frequency oscillations in giant flares from SGR 1806-20 and SGR 1900+14 may be the first direct detection of vibrations in a neutron star crust. If this interpretation is correct it offers a novel means of testing the neutron star equation of state, crustal breaking strain, and magnetic field configuration. Using timing data from RHESSI, we have confirmed the detection of a 92.5 Hz Quasi-Periodic Oscillation (QPO) in the tail of the SGR 1806-20 giant flare. We also find another, stronger, QPO at higher energies, at 626.5 Hz. Both QPOs are visible only at particular (but different) rotational phases, implying an association with a specific area of the neutron star surface or magnetosphere. At lower frequencies we confirm the detection of an 18 Hz QPO, at the same rotational phase as the 92.5 Hz QPO, and report the additional presence of a broad 26 Hz QPO. We are however unable to make a robust confirmation of the presence of a 30 Hz QPO, despite higher countrates. We discuss our results in the light of neutron star vibration models.

In this paper we discuss the spin-equilibrium of accreting neutron stars in LMXBs. We demonstrate that, when combined with a naive spin-up torque, the observed data leads to inferred magnetic fields which are at variance with those of galactic millisecond radiopulsars. This indicates the need for either additional spin-down torques (eg. gravitational radiation) or an improved accretion model. We show that a simple consistent accretion model can be arrived at by accounting for radiation pressure in rapidly accreting systems (above a few percent of the Eddington accretion rate). In our model the inner disk region is thick and significantly sub-Keplerian, and the estimated equilibrium periods are such that the LMXB neutron stars have properties that accord well with the galactic millisecond radiopulsar sample. The implications for future gravitational-wave observations are also discussed briefly.

We report the discovery of complex high frequency variability during the August 27, 1998 giant flare from SGR 1900+14 using the Rossi X-ray Timing Explorer (RXTE). We detect an 84 Hz oscillation (QPO) during a 1 s interval beginning approximately 1 min after the initial hard spike. The modulation amplitude is energy dependent, reaching a maximum of 26% (rms) for photons above 30 keV, and is not detected below 11 keV, with a 90% confidence upper limit of 14% (rms). Remarkably, additional QPOs are detected in the average power spectrum of data segments centered on the rotational phase at which the 84 Hz signal was detected. Two signals, at 53.5 and 155.1 Hz, are strongly detected, while a third feature at 28 Hz is found with lower significance. These QPOs are not detected at other rotational phases. The phenomenology seen in the SGR 1900+14 flare is similar to that of QPOs recently reported by Israel et al. from the December 27, 2004 flare from SGR 1806-20, suggesting they may have a common origin, perhaps torsional vibrations of the neutron star crust. Indeed, an association of the four frequencies (in increasing order) found in SGR 1900+14 with l = 2, 4, 7, and 13 toroidal modes appears plausible. We discuss our findings in the context of this model and show that if the stars have similar masses then the magnetic field in SGR 1806-20 must be about twice as large as in SGR 1900+14, broadly consistent with magnetic field estimates from pulse timing.

The accreting millisecond pulsar XTE J1814-338 exhibits oscillations at the known spin frequency during Type I X-ray bursts. The properties of the burst oscillations reflect the nature of the thermal asymmetry on the stellar surface. We present an analysis of the variability of the burst oscillations of this source, focusing on three characteristics: fractional amplitude, harmonic content and frequency. Fractional amplitude and harmonic content constrain the size, shape and position of the emitting region, whilst variations in frequency indicate motion of the emitting region on the neutron star surface. We examine both long-term variability over the course of the outburst, and short-term variability during the bursts. For most of the bursts, fractional amplitude is consistent with that of the accretion pulsations, implying a low degree of fuel spread. There is however a population of bursts whose fractional amplitudes are substantially lower, implying a higher degree of fuel spread, possibly forced by the explosive burning front of a precursor burst. For the first harmonic, substantial differences between the burst and accretion pulsations suggest that hotspot geometry is not the only mechanism giving rise to harmonic content in the latter. Fractional amplitude variability during the bursts is low; we cannot rule out the hypothesis that the fractional amplitude remains constant for bursts that do not exhibit photospheric radius expansion (PRE). There are no significant variations in frequency in any of the bursts except for the one burst that exhibits PRE. This burst exhibits a highly significant but small ($\approx 0.1$Hz) drop in frequency in the burst rise. The timescale of the frequency shift is slower than simple burning layer expansion models predict, suggesting that other mechanisms may be at work.

Recent numerical simulations indicate the presence of dynamical instabilities of the f-mode in differentially rotating stars even at very low values of $T/|W|$, the ratio of kinetic to potential energy. In this Letter we argue that these may be shear instabilities which occur when the degree of differential rotation exceeds a critical value and the f-mode develops a corotation point associated with the presence of a continuous spectrum. Our explanation, which is supported by detailed studies of a simple shell model, offers a straightforward way of understanding all of the key features of these instabilities.

An understanding of the oscillations of differentially rotating systems is key to many areas of astrophysics. It is of particular relevance to the emission of gravitational waves from oscillating neutron stars, which are expected to possess significant differential rotation immediately after birth or binary merger. In a previous paper we analysed the normal modes of a simple system exhibiting differential rotation. In this complementary paper we address the initial value problem for the same simple model using both analytical methods and numerical time evolutions. We derive a necessary and sufficient condition for dynamical shear instability. We discuss the dynamical behaviour of the continuous spectrum in response to an initial perturbation, and show that certain singular solutions within the continuous spectrum appear physically indistinguishable from the discrete modes outside the continuous spectrum.

An understanding of the dynamics of differentially rotating systems is key to many areas of astrophysics. We investigate the oscillations of a simple system exhibiting differential rotation, and discuss issues concerning the role of corotation points and the emergence of dynamical instabilities. This problem is of particular relevance to the emission of gravitational waves from oscillating neutron stars, which are expected to possess significant differential rotation immediately after birth or binary merger.

The loss of angular momentum due to unstable r-modes in hot young neutron stars has been proposed as a mechanism for achieving the spin rates inferred for young pulsars. One factor that could have a significant effect on the action of the r-mode instability is fallback of supernova remnant material. The associated accretion torque could potentially counteract any gravitational-wave induced spin-down, and accretion heating could affect the viscous damping rates and hence the instability. We discuss the effects of various external agents on the r-mode instability scenario within a simple model of supernova fallback on to a hot young magnetized neutron star. We find that the outcome depends strongly on the strength of the star's magnetic field. Our model is capable of generating spin rates for young neutron stars that accord well with initial spin rates inferred from pulsar observations. The combined action of r-mode instability and fallback appears to cause the spin rates of neutron stars born with very different spin rates to converge, on a timescale of about a year. The results suggest that stars with magnetic fields $\le10^{13}$ G could emit a detectable gravitational wave signal for perhaps several years after the supernova event. Stars with higher fields (magnetars) are unlikely to emit a detectable gravitational wave signal via the r-mode instability. The model also suggests that the r-mode instability could be extremely effective in preventing young neutron stars from going dynamically unstable to the bar-mode.