We report the detection of an infrared burst lagging a thermonuclear Type I X-ray burst from the accreting neutron star 4U 1728-34 (GX 354-0). Observations were performed simultaneously with XMM-Newton (0.7-12 keV), NuSTAR (3-79 keV) and HAWK-I@VLT (2.2$\mu$m). We measure a lag of $4.75 \pm 0.5$ s between the peaks of the emission in the two bands. Due to the length of the lag and the shape of the IR burst, we found that the most plausible cause for such a large delay is reprocessing of the Type I burst X-rays by the companion star. The inferred distance between the neutron star and the companion can be used to constrain the orbital period of the system, which we find to be larger than $\sim$ 66 minutes (or even $\gtrsim$ 2 hours, for a realistic inclination $< 75^\circ$). This is much larger than the current tentatively estimated period of $\sim 11$ minutes. We discuss the physical implications on the nature of the binary and conclude that most likely the companion of 4U 1728-34 is a helium star.

We investigate the structure of laterally-propagating flames through the highly-stratified burning layer in an X-ray burst. Two-dimensional hydrodynamics simulations of flame propagation are performed through a rotating plane-parallel atmosphere, exploring the structure of the flame. We discuss the approximations needed to capture the length and time scales at play in an X-ray burst and describe the flame acceleration observed. Our studies complement other multidimensional studies of burning in X-ray bursts.

We present the first realistic 3D simulations of flame front instabilities during type I X-ray bursts. The unperturbed front is characterised by the balance between the pressure gradient and the Coriolis force of a spinning neutron star ({\nu} = 450 Hz in our case). This balance leads to a fast horizontal velocity field parallel to the flame front. This flow is strongly sheared in the vertical direction. When we perturb the front an instability quickly corrugates the front. We identify this instability as the baroclinic instability. Most importantly, the flame is not disrupted by the instability and there are two major consequences: the overall flame propagation speed is {\sim} 10 times faster than in the unperturbed case and distinct flame vortices appear. The speedup is due to the corrugation of the front and the dynamics of the vortices. These vortices may also be linked to the oscillations observed in the lightcurves of the bursts.

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.

We discuss the challenges of modeling X-ray bursts in multi-dimensions, review the different calculations done to date, and discuss our new set of ongoing simulations. We also describe algorithmic improvements that may help in the future to offset some of the expense of these simulations, and describe what may be possible with exascale computing.

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.

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.

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.

EXO 1745-248 is a transient neutron-star low-mass X-ray binary that resides in the globular cluster Terzan 5. We studied the transient during its quiescent state using 18 Chandra observations of the cluster acquired between 2003 and 2016. We found an extremely variable source, with a luminosity variation in the 0.5-10 keV energy range of $\sim3$ orders of magnitude (between $3\times10^{31}$ erg s$^{-1}$ and $2\times10^{34}$ erg s$^{-1}$) on timescales from years down to only a few days. Using an absorbed power-law model to fit its quiescent spectra, we obtained a typical photon index of $\sim1.4$, indicating that the source is even harder than during outburst and much harder than typical quiescent neutron stars if their quiescent X-ray spectra are also described by a single power-law model. This indicates that EXO 1745-248 is very hard throughout the entire observed X-ray luminosity range. At the highest luminosity, the spectrum fits better when an additional (soft) component is added to the model. All these quiescent properties are likely related to strong variability in the low-level accretion rate in the system. However, its extreme variable behavior is strikingly different from the one observed for other neutron star transients that are thought to still accrete in quiescence. We compare our results to these systems. We also discuss similarities and differences between our target and the transitional millisecond pulsar IGR J18245-2452, which also has hard spectra and strong variability during quiescence.

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

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.

The paradigm in which magnetic fields play a crucial role in launching/collimating outflows in many astrophysical objects continues to gain support. However, semi-analytical models including the effect of magnetic fields on the dynamics and morphology of jets are still missing due to the intrinsic difficulties in integrating the equations describing a collimated, relativistic flow in the presence of gravity. Only few solutions have been found so far, due to the highly nonlinear character of the equations together with the need to blindly search for singularities. These numerical problems prevented a full exploration of the parameter space. We present a new integration scheme to solve r-self-similar, stationary, axisymmetric magnetohydrodynamics equations describing collimated, relativistic outflows crossing smoothly all the singular points (Alfven point and modified slow/fast points). For the first time, we are able to integrate from the disk mid-plane to downstream of the modified fast point. We discuss an ensemble of jet solutions, emphasising trends and features that can be compared to observables. We present, for the first time with a semi-analytical MHD model, solutions showing counter-rotation of the jet for a substantial fraction of its extent. We find diverse jet configurations with bulk Lorentz factors up to 10 and potential sites for recollimation between $10^3-10^7$ gravitational radii. Such extended coverage of the intervals of quantities, such as magnetic-to-thermal energy ratios at the base or the heights/widths of the recollimation region, makes our solutions suitable for application to many different systems where jets are launched.

The Be/X-ray transient 4U 0115+63 exhibited a giant, type-II outburst in October 2015. The source did not decay to its quiescent state but settled in a meta-stable plateau state (a factor ~10 brighter than quiescence) in which its luminosity slowly decayed. We used XMM-Newton to observe the system during this phase and we found that its spectrum can be well described using a black-body model with a small emitting radius. This suggests emission from hot spots on the surface, which is confirmed by the detection of pulsations. In addition, we obtained a relatively long (~7.9 ksec) Swift/XRT observation ~35 days after our XMM-Newton one. We found that the source luminosity was significantly higher and, although the spectrum could be fitted with a black-body model the temperature was higher and the emitting radius smaller. Several weeks later the system started a sequence of type-I accretion outbursts. In between those outbursts, the source was marginally detected with a luminosity consistent with its quiescent level. We discuss our results in the context of the three proposed scenarios (accretion down to the magnestospheric boundary, direct accretion onto neutron star magnetic poles or cooling of the neutron star crust) to explain the plateau phase.

In this work we discuss the recent criticism by Zdziarski of the maximal jet model derived in Falcke & Biermann (1995). We agree with Zdziarski that in general a jet's internal energy is not bounded by its rest-mass energy density. We describe the effects of the mistake on conclusions that have been made using the maximal jet model and show when a maximal jet is an appropriate assumption. The maximal jet model was used to derive a 1-D hydrodynamical model of jets in agnjet, a model that does multiwavelength fitting of quiescent/hard state X-ray binaries and low-luminosity active galactic nuclei. We correct algebraic mistakes made in the derivation of the 1-D Euler equation and relax the maximal jet assumption. We show that the corrections cause minor differences as long as the jet has a small opening angle and a small terminal Lorentz factor. We find that the major conclusion from the maximal jet model, the jet-disk symbiosis, can be generally applied to astrophysical jets. We also show that isothermal jets are required to match the flat radio spectra seen in low-luminosity X-ray binaries and active galactic nuclei, in agreement with other works.

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 search the literature for reports on the spectral properties of neutron-star low-mass X-ray binaries when they have accretion luminosities between 1E34 and 1E36 ergs/s. We found that in this luminosity range the photon index (obtained from fitting a simple absorbed power-law in the 0.5-10 keV range) increases with decreasing 0.5-10 keV X-ray luminosity (i.e., the spectrum softens). Such behaviour has been reported before for individual sources, but here we demonstrate that very likely most (if not all) neutron-star systems behave in a similar manner and possibly even follow a universal relation. When comparing the neutron-star systems with black-hole systems, it is clear that most black-hole binaries have significantly harder spectra at luminosities of 1E34 - 1E35 erg/s. Despite a limited number of data points, there are indications that these spectral differences also extend to the 1E35 - 1E36 erg/s range. This observed difference between the neutron-star binaries and black-hole ones suggests that the spectral properties (between 0.5-10 keV) at 1E34 - 1E35 erg/s can be used to tentatively determine the nature of the accretor in unclassified X-ray binaries. We discuss our results in the context of properties of the accretion flow at low luminosities and we suggest that the observed spectral differences likely arise from the neutron-star surface becoming dominantly visible in the X-ray spectra. We also suggest that both the thermal component and the non-thermal component might be caused by low-level accretion onto the neutron-star surface for luminosities below a few times 1E34 erg/s.

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 thermonuclear X-ray bursts on accreting neutron stars. For a summary, we refer to the paper.

We study the exceptionally short (32-41 ms) precursors of two intermediate-duration thermonuclear X-ray bursts observed with RXTE from the neutron stars in 4U 0614+09 and 2S 0918-549. They exhibit photon fluxes that surpass those at the Eddington limit later in the burst by factors of 2.6 to 3.1. We are able to explain both the short duration and the super-Eddington flux by mildly relativistic outflow velocities of 0.1$c$ to 0.3$c$ subsequent to the thermonuclear shell flashes on the neutron stars. These are the highest velocities ever measured from any thermonuclear flash. The precursor rise times are also exceptionally short: about 1 ms. This is inconsistent with predictions for nuclear flames spreading laterally as deflagrations and suggests detonations instead. This is the first time that a detonation is suggested for such a shallow ignition column depth ($y_{\rm ign}$ = 10$^{10}$ g cm$^{-2}$). The detonation would possibly require a faster nuclear reaction chain, such as bypassing the alpha-capture on $^{12}$C with the much faster $^{12}$C(p,$\gamma$)$^{13}$N($\alpha$,p)$^{16}$O process previously proposed. We confirm the possibility of a detonation, albeit only in the radial direction, through the simulation of the nuclear burning with a large nuclear network and at the appropriate ignition depth, although it remains to be seen whether the Zel'dovich criterion is met. A detonation would also provide the fast flame spreading over the surface of the neutron star to allow for the short rise times. (...) As an alternative to the detonation scenario, we speculate on the possibility that the whole neutron star surface burns almost instantly in the auto-ignition regime. This is motivated by the presence of 150 ms precursors with 30 ms rise times in some superexpansion bursts from 4U 1820-30 at low ignition column depths of ~10$^8$ g cm$^{-2}$.

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.

In gravitationally stratified fluids, length scales are normally much greater in the horizontal direction than in the vertical one. When modelling these fluids it can be advantageous to use the hydrostatic approximation, which filters out vertically propagating sound waves and thus allows a greater timestep. We briefly review this approximation, which is commonplace in atmospheric physics, and compare it to other approximations used in astrophysics such as Boussinesq and anelastic, finding that it should be the best approximation to use in context such as radiative stellar zones, compact objects, stellar or planetary atmospheres and other contexts. We describe a finite-difference numerical scheme which uses this approximation, which includes magnetic fields.

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

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.