Accreting Neutron Stars Research Articles

INTERFERENCE AS AN ORIGIN OF THE PEAKED NOISE IN ACCRETING X-RAY BINARIES

ABSTRACT We propose a physical model for the peaked noise in the X-ray power density spectra of accreting X-ray binaries. We interpret its appearance as an interference of two Comptonization continua: one coming from the upscattering of seed photons from the cold thin disk and the other fed by the synchrotron emission of the hot flow. Variations of both X-ray components are caused by fluctuations in mass accretion rate, but there is a delay between them corresponding to the propagation timescale from the disk Comptonization radius to the region of synchrotron Comptonization. If the disk and synchrotron Comptonization are correlated, the humps in the power spectra are harmonically related and the dips between them appear at frequencies related as odd numbers 1:3:5. If they are anti-correlated, the humps are related as 1:3:5, but the dips are harmonically related. Similar structures are expected to be observed in accreting neutron star binaries and supermassive black holes. The delay can be easily recovered from the frequency of peaked noise and further used to constrain the combination of the viscosity parameter and disk height-to-radius ratio α(H/R)2 of the accretion flow. We model multi-peak power spectra of black hole X-ray binaries GX 339–4 and XTE J1748–288 to constrain these parameters.

OCEAN g-MODES ON TRANSIENT NEUTRON STARS

ABSTRACT The neutron star ocean is a plasma of ions and electrons that extends from the base of the neutron star’s envelope to a depth where the plasma crystallizes into a solid crust. During an accretion outburst in an X-ray transient, material accumulates in the envelope of the neutron star primary. This accumulation compresses the neutron star’s outer layers and induces nuclear reactions in the ocean and crust. Accretion-driven heating raises the ocean’s temperature and increases the frequencies of -modes in the ocean; when accretion halts, the ocean cools and ocean -mode frequencies decrease. If the observed low-frequency quasi-periodic oscillations on accreting neutron stars are -modes in the ocean, the observed quasi-periodic oscillation frequencies will increase during the outburst—reaching a maximum when the ocean temperature reaches steady state—and subsequently decrease during quiescence. For time-averaged accretion rates during outbursts between the predicted -mode fundamental n = 1 l = 2 frequency is between ≈3–7 Hz for slowly rotating neutron stars. Accreting neutron stars that require extra shallow heating, such as the Z-sources MAXI J0556-332, MXB 1659-29, and XTE J1701-462, have predicted -mode fundamental frequencies between ≈3–16 HZ. Therefore, observations of low-frequency quasi-periodic oscillations between in these sources, or in other transients that require shallow heating, will support a -mode origin for the observed quasi-periodic oscillations.

Gravitational waves from surface inhomogeneities of neutron stars

Surface asymmetries of accreting neutron stars are investigated for their mass quadrupole moment content. Though the amplitude of the gravitational waves from such asymmetries seem to be beyond the limit of detectability of the present generation of detectors, it appears that rapidly rotating neutron stars with strong magnetic fields residing in HMXBs would be worth considering for targeted search for continuous gravitational waves with the next generation of instruments.

DISCOVERY OF COHERENT PULSATIONS FROM THE ULTRALUMINOUS X-RAY SOURCE NGC 7793 P13

ABSTRACT We report the detection of coherent pulsations from the ultraluminous X-ray source (ULX) NGC 7793 P13. The ≈0.42 s nearly sinusoidal pulsations were initially discovered in broadband X-ray observations using XMM-Newton and NuSTAR taken in 2016. We subsequently also found pulsations in archival XMM-Newton data taken in 2013 and 2014. The significant (≫5σ) detection of coherent pulsations demonstrates that the compact object in P13 is a neutron star, and given the observed peak luminosity of ≈1040 erg s − 1 (assuming isotropy), it is well above the Eddington limit for a 1.4 M ⊙ accretor. This makes P13 the second ULX known to be powered by an accreting neutron star. The pulse period varies between epochs, with a slow but persistent spin-up over the 2013–2016 period. This spin-up indicates a magnetic field of B ≈ 1.5 × 1012 G, typical of many Galactic accreting pulsars. The most likely explanation for the extreme luminosity is a high degree of beaming; however, this is difficult to reconcile with the sinusoidal pulse profile.

A NEW γ-RAY LOUD, ECLIPSING LOW-MASS X-RAY BINARY

ABSTRACT We report the discovery of an eclipsing low-mass X-ray binary at the center of the 3FGL error ellipse of the unassociated Fermi/Large Area Telescope γ-ray source 3FGL J0427.9–6704. Photometry from OGLE and the SMARTS 1.3 m telescope and spectroscopy from the SOAR telescope have allowed us to classify the system as an eclipsing low-mass X-ray binary (P = 8.8 hr) with a main-sequence donor and a neutron-star accretor. Broad double-peaked H and He emission lines suggest the ongoing presence of an accretion disk. Remarkably, the system shows separate sets of absorption lines associated with the accretion disk and the secondary, and we use their radial velocities to find evidence for a massive (∼1.8–1.9 M ⊙) neutron-star primary. In addition to a total X-ray eclipse with a duration of ∼2200 s observed with NuSTAR, the X-ray light curve also shows properties similar to those observed among known transitional millisecond pulsars: short-term variability, a hard power-law spectrum ( ), and a comparable 0.5–10 keV luminosity ( erg s−1). We find tentative evidence for a partial ( ) γ-ray eclipse at the same phase as the X-ray eclipse, suggesting the γ-ray emission may not be confined to the immediate region of the compact object. The favorable inclination of this binary is promising for future efforts to determine the origin of γ-rays among accreting neutron stars.

The spin period formation of millisecond pulsar by the torques of accretion and gravitational wave emission

We investigated the spin period evolution of the accreting neutron star (NS) in binary systems, based on the accretion-induced magnetic decay model, while both the accretion spin-up torque and the gravitational wave (GW) emission induced spin-down torque are taken into account. We found that the spin period of millisecond pulsar (MSP) can stop at the value of several milliseconds, if the accretion spin-up torque balances the spin-down torque by the GW radiation. Furthermore, we obtained the minimum spin period of MSP and its relation to the deformation ellipticity of NS that accounts for the NS mass quadrupole moment. The comparisons between the pulsar observations and the model results are discussed in the diagram of magnetic field versus spin period, and the consistency can be obtained.

DEPENDENCE OF X-RAY BURST MODELS ON NUCLEAR REACTION RATES

ABSTRACT X-ray bursts are thermonuclear flashes on the surface of accreting neutron stars, and reliable burst models are needed to interpret observations in terms of properties of the neutron star and the binary system. We investigate the dependence of X-ray burst models on uncertainties in (p, γ), (α, γ), and (α, p) nuclear reaction rates using fully self-consistent burst models that account for the feedbacks between changes in nuclear energy generation and changes in astrophysical conditions. A two-step approach first identified sensitive nuclear reaction rates in a single-zone model with ignition conditions chosen to match calculations with a state-of-the-art 1D multi-zone model based on the Kepler stellar evolution code. All relevant reaction rates on neutron-deficient isotopes up to mass 106 were individually varied by a factor of 100 up and down. Calculations of the 84 changes in reaction rate with the highest impact were then repeated in the 1D multi-zone model. We find a number of uncertain reaction rates that affect predictions of light curves and burst ashes significantly. The results provide insights into the nuclear processes that shape observables from X-ray bursts, and guidance for future nuclear physics work to reduce nuclear uncertainties in X-ray burst models.

Ultraluminous X-ray bursts in two ultracompact companions to nearby elliptical galaxies

A flaring X-ray source was found near the galaxy NGC 4697 (ref. 1). Two brief flares were seen, separated by four years. During each flare, the flux increased by a factor of 90 on a timescale of about one minute. There is no associated optical source at the position of the flares, but if the source was at the distance of NGC 4697, then the luminosities of the flares were greater than 1039 erg per second. Here we report the results of a search of archival X-ray data for 70 nearby galaxies looking for similar flares. We found two ultraluminous flaring sources in globular clusters or ultracompact dwarf companions of parent elliptical galaxies. One source flared once to a peak luminosity of 9 × 1040 erg per second; the other flared five times to 1040 erg per second. The rise times of all of the flares were less than one minute, and the flares then decayed over about an hour. When not flaring, the sources appear to be normal accreting neutron-star or black-hole X-ray binaries, but they are located in old stellar populations, unlike the magnetars, anomalous X-ray pulsars or soft γ repeaters that have repetitive flares of similar luminosities.

Constraining the dipolar magnetic field of M82 X-2 by the accretion model

Abstract Recently, ultraluminous X-ray source (ULX) M82 X-2 has been identified to be an accreting neutron star, which has a P = 1.37 s spin period, and is spinning up at a rate $\dot{P}=-2.0\times 10^{-10}\,\rm s\,s^{-1}$. Interestingly, its isotropic X-ray luminosity Liso = 1.8 × 1040 erg s− 1 during outbursts is 100 times the Eddington limit for a 1.4 M⊙ neutron star. In this Letter, based on the standard accretion model we attempt to constrain the dipolar magnetic field of the pulsar in ULX M82 X-2. Our calculations indicate that the accretion rate at the magnetospheric radius must be super-Eddington during outbursts. To support such a super-Eddington accretion, a relatively high multipole field ( ≳ 1013 G) near the surface of the accretor is invoked to produce an accreting gas column. However, our constraint shows that the surface dipolar magnetic field of the pulsar should be in the range of 1.0−3.5 × 1012 G. Therefore, our model supports that the neutron star in ULX M82 X-2 could be a low-magnetic-field magnetar (proposed by Tong) with a normal dipolar field (∼1012 G) and relatively strong multipole field. For the large luminosity variations of this source, our scenario can also present a self-consistency interpretation.

EVIDENCE OF SPREADING LAYER EMISSION IN A THERMONUCLEAR SUPERBURST

ABSTRACT When a neutron star (NS) accretes matter from a companion star in a low-mass X-ray binary, the accreted gas settles onto the stellar surface through a boundary/spreading layer. On rare occasions the accumulated gas undergoes a powerful thermonuclear superburst powered by carbon burning deep below the NS atmosphere. In this paper, we apply the non-negative matrix factorization spectral decomposition technique to show that the spectral variations during a superburst from 4U 1636–536 can be explained by two distinct components: (1) the superburst emission characterized by a variable temperature blackbody radiation component and (2) a quasi-Planckian component with a constant, ∼2.5 keV, temperature varying by a factor of ∼15 in flux. The spectrum of the quasi-Planckian component is identical in shape and characteristics to the frequency-resolved spectra observed in the accretion/persistent spectrum of NS low-mass X-ray binaries and agrees well with the predictions of the spreading layer model by Inogamov & Sunyaev. Our results provide yet more observational evidence that superbursts—and possibly also normal X-ray bursts—induce changes in the disc–star boundary.

The nature of the X-ray pulsar in M 31: An intermediate-mass X-ray binary?

Abstract The first finding of the spin period of an accreting neutron star in M 31 was recently reported. The observed spin period is 1.2 s, and it shows 1.27 d modulations due to orbital motion. From the orbital information, the mass donor could not be a giant massive star. On the other hand, its observed properties are very odd as those of typical low-mass X-ray binaries. In this study, we compare the observed binary parameters with theoretical models given by a stellar evolution track, and give a restriction on the possible mass range of the donor. According to the standard stellar evolution model, the donor star should be larger than 1.5 M⊙, which suggests that this system is a new member of a rare category, an intermediate-mass X-ray binary. The magnetic field strength of the neutron star suggested by the spin-up/down tendency in this system supports the possibility of an intermediate-mass donor.

Searching for supergiant fast X-ray transients withSwift

Supergiant fast X-ray transients (SFXTs) are high mass X-ray binaries (HMXBs) hosting a neutron star and an OB supergiant companion. We examine the available Swift data, as well as other new or archival/serendipitous data, on three sources: IGR J17407-2808, 2XMM J185114.3-000004, and IGR J18175-2419, whose X-ray characteristics qualify them as candidate SFXT, in order to explore their properties and test whether they are consistent with an SFXT nature. As IGR J17407-2808 and 2XMM J185114.3-000004 triggered the Burst Alert Telescope on board Swift, the Swift data allow us to provide their first arcsecond localisations, leading to an unequivocal identification of the source CXOU J174042.0-280724 as the soft X-ray counterpart of IGR J17407-2808, as well as their first broadband spectra, which can be fit with models generally describing accreting neutron stars in HMXBs. While still lacking optical spectroscopy to assess the spectral type of the companion, we propose 2XMM J185114.3-000004 as a very strong SFXT candidate. The nature of IGR J17407-2808 remains, instead, more uncertain. Its broad band properties cannot exclude that the emission originates from either a HMXB (and in that case, a SFXT) or, more likely, a low mass X-ray binary. Finally, based on the deep non-detection in our XRT monitoring campaign and a careful reanalysis of the original Integral data in which the discovery of the source was first reported, we show that IGR J18175-2419 is likely a spurious detection.

TheChandraACIS Timing Survey Project: glimpsing a sample of faint X-ray pulsators

We report on the discovery of 41 new pulsating sources in the data of the Chandra Advanced CCD Imaging Spectrometer, which is sensitive to X-ray photons in the 0.3-10 keV band. The archival data of the first 15 years of Chandra observations were retrieved and analysed by means of fast Fourier transforms, employing a peak-detection algorithm able to screen candidate signals in an automatic fashion. We carried out the search for new X-ray pulsators in light curves with more than 50 photons, for a total of about 190,000 lightcurves out of about 430,000 extracted. With these numbers, the ChAndra Timing Survey at Brera And Roma astronomical observatories (CATS@BAR) - as we called the project - represents the largest ever systematic search for coherent signals in the classic X-ray band. More than 50 per cent of the signals were confirmed by further Chandra (for those sources with two or more pointings), XMM-Newton or ROSAT data. The period distribution of the new X-ray pulsators above about 2,000s resembles that of cataclysmic variables, while there is a paucity of sources with shorter period and low fluxes. Since there is not an obvious bias against these detections, a possible interpretation is in terms of a magnetic gating mechanism in accreting neutron stars. Finally, we note that CATS@BAR is a living project and the detection algorithm will continue to be routinely applied to the new Chandra data as they become public. Based on the results obtained so far, we expect to discover about three new pulsators every year.

ACCRETION DISK SIGNATURES IN TYPE I X-RAY BURSTS: PROSPECTS FOR FUTURE MISSIONS

ABSTRACT Type I X-ray bursts and superbursts from accreting neutron stars illuminate the accretion disk and produce a reflection signal that evolves as the burst fades. Examining the evolution of reflection features in the spectra will provide insight into the burst–disk interaction, a potentially powerful probe of accretion disk physics. At present, reflection has been observed during only two bursts of exceptional duration. We investigate the detectability of reflection signatures with four of the latest well-studied X-ray observatory concepts: Hitomi, Neutron Star Interior Composition Explorer (NICER), Athena, and Large Observatory For X-ray Timing (LOFT). Burst spectra are modeled for different values for the flux, temperature, and the disk ionization parameter, which are representative for most known bursts and sources. The effective area and throughput of a Hitomi-like telescope are insufficient for characterizing burst reflection features. NICER and Athena will detect reflection signatures in Type I bursts with peak fluxes ≳10−7.5 erg cm−2 s−1 and also effectively constrain the reflection parameters for bright bursts with fluxes of ∼10−7 erg cm−2 s−1 in exposures of several seconds. Thus, these observatories will provide crucial new insight into the interaction of accretion flows and X-ray bursts. For sources with low line-of-sight absorption, the wide bandpass of these instruments allows for the detection of soft X-ray reflection features, which are sensitive to the disk metallicity and density. The large collecting area that is part of the LOFT design would revolutionize the field by tracing the evolution of the accretion geometry in detail throughout short bursts.

Thermal conductivity and impurity scattering in the accreting neutron star crust

We calculate the thermal conductivity of electrons for the strongly correlated multi-component ion plasma expected in the outer layers of neutron star's crust employing a Path Integral Monte Carlo (PIMC) approach. This allows us to isolate the low energy response of the ions and use it to calculate the electron scattering rate and the electron thermal conductivity. We find that the scattering rate is enhanced by a factor 2-4 compared to earlier calculations based on the simpler electron-impurity scattering formalism. This findings directly impacts the interpretation of thermal relaxation observed in transiently accreting neutron stars and has implications for the composition and nuclear reactions in the crust that occur during accretion.

The magnetic field evolution of ULX NuSTAR J095551+6940.8 in M82 – a legacy of accreting magnetar

Ultra-luminous X-ray sources are usually believed to be black holes with mass about $10^{2-3}M_{\odot}$. However, the recent discovery of NuSTAR J095551+6940.8 in M82 by Bachetti et al. shows that it holds the spin period $P=1.37\rm\,s$ and period {\bf derivative} $\dot{P}\approx-2\times10^{-10}\rm\,s\,s^{-1}$, which provides a strong evidence that some ultra-luminous X-ray sources could be neutron stars. We obtain that the source may be an evolved magetar, according to our simulation by employing the model of accretion induced the polar magnetic field decay and standard spin-up torque of an accreting neutron star. The results show that NuSTAR J095551+6940.8 is still in the spin-up process, and the polar magnetic field decays to about $4.5\times10^{12}\rm\,G$ after accreting $\sim 10^{-2.5}$\ms, while the strong magnetic field exists in the out-polar region, which could be responsible for the observed low field magnetar. The ultra luminosity of the source can be explained by the beaming effort and two kinds of accretion--radial random accretion and disk accretion. Since the birth rate of magnetars is about ten percent of the normal neutron stars, we guess that several ultra-luminous X-ray sources should share the similar properties to that of NuSTAR J095551+6940.8.

Swift J174540.7−290015: a new accreting binary in the Galactic Centre

We report on the identification of the new Galactic Center (GC) transient Swift J174540.7-290015 as a likely low mass X-ray binary (LMXB) located at only 16 arcsec from Sgr A*. This transient was detected on 2016 February 6th during the Swift GC monitoring, and it showed long-term spectral variations compatible with a hard to soft state transition. We observed the field with XMM-Newton on February 26th for 35 ks, detecting the source in the soft state, characterised by a low level of variability and a soft X-ray thermal spectrum with a high energy tail (detected by INTEGRAL up to ~50 keV), typical of either accreting neutron stars or black holes. We observed: i) a high column density of neutral absorbing material, suggesting that Swift J174540.7-290015 is located near or beyond the GC and; ii) a sub-Solar Iron abundance, therefore we argue that Iron is depleted into dust grains. The lack of detection of FeK absorption lines, eclipses or dipping suggests that the accretion disc is observed at a low inclination angle. Radio (VLA) observations did not detect any radio counterpart to Swift J174540.7-290015. No evidence for X-ray or radio periodicity is found. The location of the transient was observed also in the near-IR with GROND at MPG/ESO La Silla 2.2m telescope and VLT/NaCo pre- and post-outburst. Within the Chandra error region we find multiple objects that display no significant variations.

A new model for the X-ray continuum of the magnetized accreting pulsars

Accreting highly magnetized pulsars in binary systems are among the brightest X-ray emitters in our Galaxy. Although a number of high statistical quality broad-band (0.1-100 keV) X-ray observations are available, the spectral energy distribution of these sources is usually investigated by adopting pure phenomenological models, rather than models linked to the physics of accretion. In this paper, a detailed spectral study of the X-ray emission recorded from the high-mass X-ray binary pulsars Cen X-3, 4U 0115+63, and Her X-1 is carried out by using BeppoSAX and joined Suzaku+NuStar data, together with an advanced version of the compmag model. The latter provides a physical description of the high energy emission from accreting pulsars, including the thermal and bulk Comptonization of cyclotron and bremsstrahlung seed photons along the neutron star accretion column. The compmag model is based on an iterative method for solving second-order partial differential equations, whose convergence algorithm has been improved and consolidated during the preparation of this paper. Our analysis shows that the broad-band X-ray continuum of all considered sources can be self-consistently described by the compmag model. The cyclotron absorption features, not included in the model, can be accounted for by using Gaussian components. From the fits of the compmag model to the data we inferred the physical properties of the accretion columns in all sources, finding values reasonably close to those theoretically expected according to our current understanding of accretion in highly magnetized neutron stars. The updated version of the compmag model has been tailored to the physical processes that are known to occur in the columns of highly magnetized accreting neutron stars and it can thus provide a better understanding of the high energy radiation from these sources.

A SURVEY OF CHEMICAL SEPARATION IN ACCRETING NEUTRON STARS

ABSTRACT The heavy element ashes of rp-process hydrogen and helium burning in accreting neutron stars are compressed to high density where they freeze, forming the outer crust of the star. We calculate the chemical separation on freezing for a number of different nuclear mixtures resulting from a range of burning conditions for the rp-process. We confirm the generic result that light nuclei are preferentially retained in the liquid and heavy nuclei in the solid. This is in agreement with the previous study of a 17-component mixture of rp-process ashes by Horowitz et al., but extends that result to a much larger range of compositions. We also find an alternative phase separation regime for the lightest ash mixtures which does not demonstrate this generic behavior. With a few exceptions, we find that chemical separation reduces the expected in the outer crust compared to the initial rp-process ash, where measures the mean-square dispersion in atomic number Z of the nuclei in the mixture. We find that the fractional spread of Z plays a role in setting the amount of chemical separation and is strongly correlated to the divergence between the two/three-component approximations and the full component model. The contrast in Y e between the initial rp-process ashes and the equilibrium liquid composition is similar to that assumed in earlier two-component models of compositionally driven convection, except for very light compositions which produce nearly negligible convective driving. We discuss the implications of these results for observations of accreting neutron stars.

Colloquium: Measuring the neutron star equation of state using x-ray timing

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]