SAX J1808 Research Articles

Lack of emission lines in the optical spectra of SAX J1808.4–3658 during reflaring of the 2019 outburst

Aims. We present spectroscopy of the accreting X-ray binary and millisecond pulsar SAX J1808.4–3658. These observations are the first to be obtained during a reflaring phase. We collected spectroscopic data during the beginning of the reflaring of the 2019 outburst and compared them to previous datasets taken at different epochs, both of the same outburst and across the years. To this end, we also present spectra of the source taken during quiescence in 2007, one year before the next outburst. Methods. We made use of data taken by the Very Large Telescope (VLT) X-shooter spectrograph on August 31, 2019, three weeks after the outburst peak. For flux calibration, we used photometric data taken during the same night by the 1m telescopes from the Las Cumbres Observatory network that are located in Chile. We compare our spectra to the quiescent data taken by the VLT-FORS1 spectrograph in September, 2007. We inspected the spectral energy distribution by fitting our data with a multicolored accretion-disk model and sampled the posterior probability density function for the model parameters with a Markov chain Monte Carlo (MCMC) algorithm. Results. We find the optical spectra of the 2019 outburst to be unusually featureless, with no emission lines present despite the high resolution of the instrument. Fitting the UV-optical spectral energy distribution with a disk plus irradiated star model results in a very large value for the inner disk radius of ∼5130 ± 240 km, which could suggest that the disk was emptied of material during the outburst, possibly accounting for the emission-less spectra. Alternatively, the absence of emission lines could be due to a significant contribution of the jet emission at optical wavelengths.

Inferring system parameters from the bursts of the accretion-powered pulsar IGR J17498–2921

ABSTRACT Thermonuclear (type-I) bursts exhibit properties that depend both on the local surface conditions of the neutron stars on which they ignite, as well as the physical parameters of the host binary system. However, constraining the system parameters requires a comprehensive method to compare the observed bursts to simulations. We have further developed the beansp code for this purpose and analysed the bursts observed from IGR J17498$-$2921, a 401-Hz accretion-powered pulsar, discovered during it’s 2011 outburst. We find good agreement with a model having H-deficient fuel with $X=0.15\pm 0.4$, and CNO metallicity $Z=0.0014^{+0.0004}_{-0.0003}$, about a tenth of the solar value. The model has the system at a distance of $5.7^{+0.6}_{-0.5}$ kpc, with a massive (${\approx} 2\ \mathrm{ M}_\odot$) neutron star and a likely inclination of $60^\circ$. We also re-analysed the data from the 2002 outburst of the accretion-powered millisecond pulsar SAX J1808.4$-$3658. For that system we find a substantially closer distance than previously inferred, at $2.7\pm 0.3$ kpc, likely driven by a larger degree of burst emission anisotropy. The other system parameters are largely consistent with the previous analysis. We briefly discuss the implications for the evolution of these two systems.

NuSTAR and Swift observations of two supergiant fast X-ray transients: AX J1841.0−0536 and SAX J1818.6−1703

ABSTRACT Supergiant fast X-ray transients are wind-fed binaries hosting neutron star accretors, which display a peculiar variability in the X-ray domain. Different models have been proposed to explain this variability and the strength of the compact object magnetic field is generally considered a key parameter to discriminate among possible scenarios. We present here the analysis of two simultaneous observational campaigns carried out with Swift and NuSTAR targeting the supergiant fast X-ray transient sources AX J1841.0−0536 and SAX J1818.6−1703. A detailed spectral analysis is presented for both sources, with the main goal of hunting for cyclotron resonant scattering features that can provide a direct measurement of the neutron star magnetic field intensity. AX J1841.0−0536 was caught during the observational campaign at a relatively low flux. The source broad-band spectrum was featureless and could be well-described by using a combination of a hot blackbody and a power-law component with no measurable cut-off energy. In the case of SAX J1818.6−1703, the broad-band spectrum presented a relatively complex curvature which could be described by an absorbed cut-off power law (including both a cut-off and a folding energy) and featured a prominent edge at ∼7 keV, compatible with being associated to the presence of a ‘screen’ of neutral material partly obscuring the X-ray source. The fit to the broad-band spectrum also required the addition of a moderately broad (∼1.6 keV) feature centred at ∼14 keV. If interpreted as a cyclotron resonant scattering feature, our results would indicate for SAX J1818.6−1703 a relatively low-magnetized neutron star (∼1.2 × 1012 G).

The Imprint of Convection on Type I X-Ray Bursts: Pauses in Photospheric Radius Expansion Lightcurves

Motivated by the recent observation by NICER of a type I X-ray burst from SAX J1808.4–3658 with a distinct “pause” feature during its rise, we show that bursts which ignite in a helium layer underneath a hydrogen-rich shell naturally give rise to such pauses, as long as enough energy is produced to eject the outer layers of the envelope by super-Eddington winds. The length of the pause is determined by the extent of the convection generated after ignition, while the rate of change of luminosity following the pause is set by the hydrogen gradient left behind by convection. Using the MESA stellar evolution code, we simulate the accumulation, nuclear burning, and convective mixing prior to and throughout the ignition of the burst, followed by the hydrodynamic wind. We show that the results are sensitive to the treatment of convection adopted within the code. In particular, the efficiency of mixing at the H/He interface plays a key role in determining the shape of the lightcurve. The data from SAX J1808.4–3658 favor strong mixing scenarios. Multidimensional simulations will be needed to properly model the interaction between convection and nuclear burning during these bursts, which will then enable a new way to use X-ray burst lightcurves to study neutron star surfaces.

Hydrogen-triggered X-Ray Bursts from SAX J1808.4−3658? The Onset of Nuclear Burning

We present a study of weak, thermonuclear X-ray bursts from the accreting millisecond X-ray pulsar SAX J1808.4−3658. We focus on a burst observed with the Neutron Star Interior Composition Explorer on 2019 August 9, and describe a similar burst observed with the Rossi X-ray Timing Explorer in 2005 June. These bursts occurred soon after outburst onset, 2.9 and 1.1 days, after the first indications of fresh accretion. We measure peak burst bolometric fluxes of 6.98 ± 0.50 × 10−9 and 1.54 ± 0.10 × 10−8 erg cm−2 s−1, respectively, which are factors of ≈30 and 15 less than the peak flux of the brightest, helium-powered bursts observed from this source. From spectral modeling we estimate the accretion rates and accreted columns at the time of each burst. For the 2019 burst we estimate an accretion rate of M ⊙ yr−1, and a column in the range 3.9–5.1 × 107 g cm−2. For the 2005 event the accretion rate was similar, but the accreted column was half of that estimated for the 2019 burst. The low accretion rates, modest columns, and evidence for a cool neutron star in quiescence, suggest these bursts are triggered by thermally unstable CNO cycle hydrogen burning. The post-burst flux level in the 2019 event appears offset from the pre-burst level by an amount consistent with quasi-stable hydrogen burning due to the temperature-insensitive, hot-CNO cycle, further suggesting hydrogen burning as the primary fuel source. This provides strong observational evidence for hydrogen-triggered bursts. We discuss our results in the context of previous theoretical modeling.

Study of a Minimally Deformed Anisotropic Solution for Compact Objects with Massive Scalar Field in Brans–Dicke Gravity Admitting the Karmarkar Condition

In the present paper, we focused on exploring the possibility of providing a new class of exact solutions for viable anisotropic stellar systems by means of the massive Brans–Dicke (BD) theory of gravity. In this respect, we used the decoupling of gravitational sources by minimal geometric deformation (MGD) (e−η=Ψ+βh) for compact stellar objects in the realm of embedding class-one space-time to study anisotropic solutions for matter sources through the modified Einstein field equations. For this purpose, we used the ansatz for Ψ relating to the prominent, well-known and well-behaved Finch–Skea model via Karmarkar condition, and the determination scheme for deformation function h(r) was proposed via mimic requirement on radial pressure component: θ11(r)=pr(r) and matter density: θ00(r)=ρ(r) for the anisotropic sector. Moreover, we analyzed the main physical highlights of the anisotropic celestial object by executing several physical tests for the case θ11(r)=pr(r). We have clearly shown how the parameters α, β and ωBD introduced by massive BD gravity via the MGD approach incorporating the anisotropic profile of the matter distribution have an immense effect on many physical parameters of compact bodies such as LMC X-4, LMC X-4, Her X-1, 4U 1820-30, 4U 1608-52, SAX J1808.4–658 and many others that can be fitted.

The 2019 outburst of AMXP SAX J1808.4–3658 and radio follow up of MAXI J0911–655 and XTE J1701–462

ABSTRACTWe present radio coverage of the 2019 outburst of the accreting millisecond X-ray pulsar (AMXP) SAX J1808.4–3658, obtained with MeerKAT. We compare these data to contemporaneous X-ray and optical measurements in order to investigate the coupling between accretion and jet formation in this system, while the optical light curve provides greater detail of the outburst. The reflaring activity following the main outburst peak was associated with a radio re-brightening, indicating a strengthening of the jet in this phase of the outburst. We place quasi-simultaneous radio and X-ray measurements on the global radio:X-ray plane for X-ray binaries, and show they reside in the same region of luminosity space as previous outburst measurements, but significantly refine the correlation for this source. We also present upper limits on the radio emission from the AMXP MAXI J0911–655 and the transitional Z/Atoll-type transient XTE J1701–462. In the latter source, we also confirm that nearby large-scale structures reported in previous radio observations of the source are persistent over a period of ∼15 yr, and so are almost certainly background radio galaxies and not associated with the X-ray transient.

Photometric Observations of Nine High Mass X-Ray Binaries and Analysis of Potential Periodicities and Variations

Beginning with American Association of Variable Star Observers Alert Notice 348 (2007 April 6) and being continually extended and updated through 2020 April 17 with Alert Notice 377, Dr. Gordon Sarty requested continuing assistance with his ongoing radial velocity observations of High Mass X-Ray Binaries. As of Alert Notice 377, nine targets were identified. These are: 1H 1936+541, RX J2030.5+4751, V2175 Cyg, V1008 Cep, V831 Cas, RX J0440.9+4431/LS V+44 17, V420 Aur, HD 249179, and SAX J0635+0533. While these Alert Notices prioritized radial velocity observations with the objective of determining orbital periods of these binary systems, optical photometry can also reveal other periods and variations due to mass transfer, pulsations, rotations, spin rates, and other effects. This Research Note will address in more detail the observation and analysis of one of these nine High Mass X-Ray Binaries, specifically 1H 1936+541.

Timing Analysis of the 2022 Outburst of the Accreting Millisecond X-Ray Pulsar SAX J1808.4-3658: Hints of an Orbital Shrinking

We present a pulse timing analysis of NICER observations of the accreting millisecond X-ray pulsar SAX J1808.4−3658 during the outburst that started on 2022 August 19. Similar to previous outbursts, after decaying from a peak luminosity of ≃1 × 1036 erg s−1 in about a week, the pulsar entered a ∼1 month long reflaring stage. Comparison of the average pulsar spin frequency during the outburst with those previously measured confirmed the long-term spin derivative of Hz s−1, compatible with the spin-down torque of a ≈1026 G cm3 rotating magnetic dipole. For the first time in the last twenty years, the orbital phase evolution shows evidence for a decrease of the orbital period. The long-term behavior of the orbit is dominated by an ∼11 s modulation of the orbital phase epoch consistent with a ∼21 yr period. We discuss the observed evolution in terms of a coupling between the orbit and variations in the mass quadrupole of the companion star.

AstroSat observation of the accreting millisecond X-ray pulsar SAX J1808.4–3658 during its 2019 outburst

ABSTRACT We report on the analysis of the AstroSat data set of the accreting millisecond X-ray pulsar SAX J1808.4–3658, obtained during its 2019 outburst. We found coherent pulsations at ∼401 Hz and an orbital solution consistent with previous studies. The 3–20 keV pulse profile can be well fitted with three harmonically related sinusoidal components with background-corrected fractional amplitudes of $\sim 3.5 {{\ \rm per\ cent}}$, $\sim 1.2 {{\ \rm per\ cent}}$ and $\sim 0.37 {{\ \rm per\ cent}}$ for the fundamental, second and third harmonics, respectively. Our energy-resolved pulse profile evolution study indicates a strong energy dependence. We also observed a soft lag in the fundamental and hard lags during its harmonic. The broad-band spectrum of SAX J1808.4–3658 can be described well using a combination of the thermal emission component with kT ∼ 1 keV, a thermal Comptonization (Γ ∼ 1.67) from the hot corona and broad emission lines due to Fe.

Non-singular T–K axion stars with/without the dynamical bosonic field in the presence of negative Λ term

In the present work authors derive exact solutions for the relativistic compact stars in the presence of two fields axion (Dante’s Inferno model) and with/without the complex scalar field (with the quartic self-interaction) coupled to gravity. The matter source is assumed as the perfect fluid one, and we also use barotropic/MIT bag EoS to derive energy density and anisotropic/isotropic pressures from Einstein Field Equations (EFE’s). For simplicity, as the metric potentials we use Tolman–Kuchowicz metric potentials, which are non-singular and physically acceptable. Unknown variables for the T–K metric potentials were derived from the junction conditions. To examine the redshift function, adiabatic index and energy, causality conditions we used such compact stars as PSR J1416-223, PSR J1903+32, 4U 1820-30, Cen X-3, EXO 1785-248, SAX J1808.4365.

Search for continuous gravitational waves from 20 accreting millisecond x-ray pulsars in O3 LIGO data

Results are presented of searches for continuous gravitational waves from 20 accreting millisecond x-ray pulsars with accurately measured spin frequencies and orbital parameters, using data from the third observing run of the Advanced LIGO and Advanced Virgo detectors. The search algorithm uses a hidden Markov model, where the transition probabilities allow the frequency to wander according to an unbiased random walk, while the J-statistic maximum-likelihood matched filter tracks the binary orbital phase. Three narrow subbands are searched for each target, centered on harmonics of the measured spin frequency. The search yields 16 candidates, consistent with a false alarm probability of 30% per subband and target searched. These candidates, along with one candidate from an additional target-of-opportunity search done for SAX J1808.4−3658, which was in outburst during one month of the observing run, cannot be confidently associated with a known noise source. Additional follow-up does not provide convincing evidence that any are a true astrophysical signal. When all candidates are assumed nonastrophysical, upper limits are set on the maximum wave strain detectable at 95% confidence, h095%. The strictest constraint is h095%=4.7×10−26 from IGR J17062−6143. Constraints on the detectable wave strain from each target lead to constraints on neutron star ellipticity and r-mode amplitude, the strictest of which are ε95%=3.1×10−7 and α95%=1.8×10−5 respectively. This analysis is the most comprehensive and sensitive search of continuous gravitational waves from accreting millisecond x-ray pulsars to date.Received 19 September 2021Accepted 16 December 2021DOI:https://doi.org/10.1103/PhysRevD.105.022002© 2022 American Physical SocietyPhysics Subject Headings (PhySH)Research AreasGravitational wavesPhysical SystemsBinary starsNeutron stars & pulsarsGravitation, Cosmology & Astrophysics

Sensitivity Study of Nuclear Reactions Influencing Photospheric Radius Expansion X-Ray Bursts

Up to now, more than 62 of the 115 X-ray sources of low-mass-X-ray binaries have been identified as photospheric radius expansion (PRE) bursters [1]. Galloway and collaborators expect more PRE bursters in their near future analysis [2]. Although more than half of the discovered X-ray sources are PRE bursters, the bursting mechanism of PRE burster is still not adequately understood. This is because of the complicated hydrodynamics and variable accretion rates. An example is the accretion-powered millisecond pulsar SAX J1808.4–3658 [3, 4] that powered up the brightest Type-I X-ray burst (XRB) recorded by NICER in recent history [5]. The first 1D multi-zone model of SAX J1808.4–3658 was recently constructed [6, 7]. The pioneering model offers a first concurrent and direct comparison with the observed light curves, fluences, and recurrence times. With the three observables, a comparison between theory and observations could be more sensitive than the previous studies of the clocked burster and post-processing models. We perform a sensitivity study on (α,p), (α,γ), (p, α), and (p,γ) reactions with a total up to ~1,500 reactions. Our current result indicates that the observables are more sensitive to the competition between the reactions involving alpha-capture, e.g., the 22Mg(α, p) and 22Mg(p,γ) reactions competing at the 22Mg branch point [8].

Modelling spin-up episodes in accreting millisecond X-ray pulsars

ABSTRACT Accreting millisecond X-ray pulsars are known to provide a wealth of physical information during their successive states of outburst and quiescence. Based on the observed spin-up and spin-down rates of these objects, it is possible, among other things, to infer the stellar magnetic field strength and test models of accretion disc flow. In this paper, we consider the three accreting X-ray pulsars (XTE J1751–305, IGR J00291+5934 and SAX J1808.4–3658) with the best available timing data, and model their observed spin-up rates with the help of a collection of standard torque models that describe a magnetically threaded accretion disc truncated at the magnetospheric radius. Whilst none of these models is able to explain the observational data, we find that the inclusion of the physically motivated phenomenological parameter ξ, which controls the uncertainty in the location of the magnetospheric radius, leads to an enhanced disc-integrated accretion torque. These ‘new’ torque models are compatible with the observed spin-up rates as well as the inferred magnetic fields of these objects provided that ξ ≈ 0.1−0.5. Our results are supplemented with a discussion of the relevance of additional physics effects that include the presence of a multipolar magnetic field and general relativistic gravity.

A relativistic model of stellar objects with core-crust-envelope division

In this work, we present a cogent and physically well-behaved solution for neutron stars envisaged with a core layer having quark matter satisfying the MIT-bag equation of state (EoS), meso layer with Bose-Einstein condensate (BEC) matter satisfying modified BEC EoS and an envelope having neutron fluid and Coulomb liquids satisfying quadratic EoS. All the required physical and geometrical parameters like gravitational potentials, pressures, radial velocity, anisotropy, adiabatic index, mass function, compactification factor, and gravitational and surface redshift functions show a feasible trend and are continuous with smooth variation throughout the interior and across the regions of the star. Further, causality condition, energy conditions, static stability criterion (using Tolman-Oppenheimer-Volkoff equation) and Herrera cracking stability criterion are met throughout the star. The approach seems to be resulting in more realistic and accurate modeling of stellar objects, particularly realized by us for X-ray binary stars 4U 1608–52 (M = 1.7 M ⊙, R = 9.5 km) and SAX J1808.4–3658 (M = 1.2 M ⊙, R = 7.2 km). Furthermore, we have ascertained that the continuity of the stability factor in all three regions of the stars demand a smaller core. As the core region of the star increases, the stability factor becomes discontinuous at all the interfaces inside the star.

A new class of viable and exact solutions of EFEs with Karmarkar conditions: an application to cold star modeling

In this work we present a theoretical framework within Einstein’s classical general relativity which models stellar compact objects such as PSR J1614–2230 and SAX J1808.4–3658. The Einstein field equations are solved by assuming that the interior of the compact object is described by a class I spacetime. The so-called Karmarkar condition arising from this requirement is integrated to reduce the gravitational behaviour to a single generating function. By appealing to physics we adopt a form for the gravitational potential which is sufficiently robust to accurately describe compact objects. Our model satisfies all the requirements for physically realistic stellar structures.

Thermal evolution of relativistic hyperonic compact stars with calibrated equations of state

A set of unified relativistic mean-field equations of state for hyperonic compact stars recently built in [M. Fortin, Ad. R. Raduta, S. Avancini, and C. Providencia, Phys. Rev. D {\bf 101}, 034017 (2020)] is used to study the thermal evolution of non-magnetized and non-rotating spherically-symmetric isolated and accreting neutron stars under different hypothesis concerning proton $S$-wave superfluidity. These equations of state have been obtained in the following way: the slope of the symmetry energy is in agreement with experimental data; the coupling constants of $\Lambda$ and $\Xi$-hyperons are determined from experimental hypernuclear data; uncertainties in the nucleon-$\Sigma$ interaction potential are accounted for; current constraints on the lower bound of the maximum neutron star mass are satisfied. Within the considered set of equations of state, the presence of hyperons is essential for the description of the cooling/heating curves. One of the conclusions we reach is that the criterion of best agreement with observational data leads to different equations of states and proton $S$-wave superfluidity gaps when applied separately for isolated neutron stars and accreting neutron stars in quiescence. This means that at least in one situation the traditional simulation framework that we employ is not complete and/or the equations of state are inappropriate. Another result is that, considering equations of state which do not allow for nucleonic dUrca or allow for it only in very massive NS, the low luminosity of SAX J1808 requires a repulsive $\Sigma$-hyperon potential in symmetric nuclear matter in the range $U_\Sigma^{(N)}\approx 10-30$ MeV. This range of values for $U_\Sigma^{(N)} $ is also supported by the criterion of best agreement with all available data from INS and XRT.

Neutron star parameter constraints for accretion-powered millisecond pulsars from the simulated IXPE data

We have simulated the X-ray polarization data that can be obtained with the Imaging X-ray Polarimetry Explorer, when observing accretion-powered millisecond pulsars. We estimated the necessary exposure times for SAX J1808.4−3658 in order to obtain different accuracy in the measured time-dependent Stokes profiles integrated over all energy channels. We found that the measured relative errors strongly depend on the relative configuration of the observer and the emitting hotspot. The improvement in the minimum relative error in Stokes Q and U parameters as a function of observing time t scales as 1/√t, and it spans the range from 30–90% with a 200 ks exposure time to 20–60% with a 500 ks exposure time (in the case of data binned in 19 phase bins). The simulated data were also used to predict how accurate measurements of the geometrical parameters of the neutron star can be made when modelling only Q and U parameters, but not the flux. We found that the observer inclination and the hotspot co-latitude could be determined with better than 10° accuracy for most of the cases we considered. In addition, we show that the position of a secondary hotspot can also be constrained when the spot is not obscured by an accretion disc. These measurements can be used to further constrain the neutron star mass and radius when combined with modelling of the X-ray pulse profile.

Probing Jet Launching in Neutron Star X-Ray Binaries: The Variable and Polarized Jet of SAX J1808.4–3658

Abstract We report on an optical photometric and polarimetric campaign on the accreting millisecond X-ray pulsar (AMXP) SAX J1808.4–3658 during its 2019 outburst. The emergence of a low-frequency excess in the spectral energy distribution in the form of a red excess above the disk spectrum (seen most prominently in the z, i, and R bands) is observed as the outburst evolves. This is indicative of optically thin synchrotron emission due to a jet, as seen previously in this source and in other AMXPs during outburst. At the end of the outburst decay, the source enters a reflaring state. The low-frequency excess is still observed during the reflares. Our optical (BVRI) polarimetric campaign shows variable linear polarization (LP) throughout the outburst. We show that this is intrinsic to the source, with low-level but significant detections (0.2%–2%) in all bands. The LP spectrum is red during both the main outburst and the reflaring state, favoring a jet origin for this variable polarization over other interpretations, such as Thomson scattering with free electrons from the disk or the propelled matter. During the reflaring state, a few episodes with stronger LP levels (1%–2%) are observed. The low-level, variable LP is suggestive of strongly tangled magnetic fields near the base of the jet. These results clearly demonstrate how polarimetry is a powerful tool for probing the magnetic field structure in X-ray binary jets, as for active galactic nuclei jets.

Discovery of thermonuclear Type-I X-ray bursts from the X-ray binary MAXI J1807+132

ABSTRACT MAXI J1807+132 is a low-mass X-ray binary (LMXB) first detected in outburst in 2017. Observations during the 2017 outburst did not allow for an unambiguous identification of the nature of the compact object. MAXI J1807+132 that was detected in outburst again in 2019 and was monitored regularly with Neutron Star Interior Composition Explorer(NICER). In this paper, we report on 5 days of observations during which we detected three thermonuclear (Type-I) X-ray bursts, identifying the system as a neutron star LMXB. Time-resolved spectroscopy of the three Type-I bursts revealed typical characteristics expected for these phenomena. All three Type-I bursts show slow rises and long decays, indicative of mixed H/He fuel. We find no strong evidence that any of the Type-I bursts reached the Eddington Luminosity; however, under the assumption that the brightest X-ray burst underwent photospheric radius expansion, we estimate a <12.4 kpc upper limit for the distance. We searched for burst oscillations during the Type-I bursts from MAXI J1807+132 and found none (<10 per cent amplitude upper limit at 95 per cent confidence level). Finally, we found that the brightest Type-I burst shows a ∼1.6 s pause during the rise. This pause is similar to one recently found with NICER in a bright Type-I burst from the accreting millisecond X-ray pulsar SAX J1808.4–3658. The fact that Type-I bursts from both sources can show this type of pause suggests that the origin of the pauses is independent of the composition of the burning fuel, the peak luminosity of the Type-I bursts, or whether the NS is an X-ray pulsar.