Accreting Neutron Stars Research Articles

Partially accreted crusts of neutron stars

Context. Neutron stars in low-mass binary systems are subject to accretion when material originating from the companion star accumulates on the surface. In most cases, the justified and common assumption in studying the properties of the neutron star crust is the fully accreted crust approximation. However, observations of some X-ray transient sources indicate that the original crust has not been completely replaced by accreted material, but is partly composed of the compressed original crust. Aims. The crust of an accreting neutron star beyond the fully accreted crust approximation was studied; a two-part (or hybrid) crust made of the original crust that is compressed and of the accreted material crashing onto it was reconstructed as a function of the accretion stage. The differences in the composition and energy sources for the fully accreted and hybrid crusts influence the cooling and transport properties. Methods. A simple semi-empirical formula of a compressible liquid drop was used to compute the equation of state and composition of the hybrid crust. Calculations were based on the single-nucleus model, with a more accurate treatment of the neutron drip point. We compared the nuclear reactions triggered by compression in the original crust and in the accreted matter part of the hybrid crust. We discuss another crust compression astrophysical phenomenon related to spinning neutron stars. Results. The compression of the originally catalyzed outer crust triggers exothermic reactions (electron captures and pycnonuclear fusions) that deposit heat in the crust. The heat sources are cataloged as a function of the compression until the fully accreted crust approximation is reached. The pressure at which neutron drip occurs is a nonmonotonic function of the depth, leading to a temporary neutron drip anomaly. The additional potential source of energy for partially accreted crusts is the occurrence of a density inversion phenomenon between some compressed layers. Conclusions. The original crust of a neutron star cannot be neglected for the initial period of accretion, when the original crust is not fully replaced by the accreted matter. The amount of heat associated with the compression of the original crust is on the same order of magnitude as that from the sources acting in the accreted part of the hybrid crust.

Accreting neutron stars: heating of the upper layers of the inner crust

ABSTRACT Neutron stars in low-mass X-ray binaries are thought to be heated up by accretion-induced exothermic nuclear reactions in the crust. The energy release and the location of the heating sources are important ingredients of the thermal evolution models. Here, we present thermodynamically consistent calculations of the energy release in three zones of the stellar crust: at the outer–inner crust interface, in the upper layers of the inner crust (up to the density ρ ≤ 2 × 1012 g cm−3), and in the underlying crustal layers. We consider three representative models of thermonuclear ashes (superburst, extreme rp, and Kepler ashes). The energy release in each zone is parametrized by the pressure at the outer–inner crust interface, Poi, which encodes all uncertainties related to the physics of the deepest inner-crust layers. Our calculations allow us to set new theoretical lower limits on the net energy release (per accreted baryon): Q ≳ 0.28 MeV for extreme rp ashes and Q ≳ 0.43–0.51 MeV for superburst and Kepler ashes. Our results can be directly incorporated into numerical codes and provide an opportunity to constrain Poi by comparing thermal evolution models of accreting neutron stars with observations.

Localized thermonuclear bursts from accreting magnetic white dwarfs

Nova explosions are caused by global thermonuclear runaways triggered in the surface layers of accreting white dwarfs1-3. It has been predicted4-6 that localized thermonuclear bursts on white dwarfs can also take place, similar to type-I X-ray bursts observed in accreting neutron stars. Unexplained rapid bursts from the binary system TV Columbae, in which mass is accreted onto a moderately strong magnetized white dwarf from a low-mass companion, have been observed on several occasions in the past 40 years7-11. During these bursts, the optical/ultraviolet luminosity increases by a factor of more than three in less than an hour and fades in around ten hours. Fast outflows have been observed in ultraviolet spectral lines7, with velocities of more than 3,500 kilometres per second, comparable to the escape velocity from the white dwarf surface. Here we report on optical bursts observed in TV Columbae and in two additional accreting systems, EI Ursae Majoris and ASASSN-19bh. The bursts have a total energy of approximately 10-6 times thanthose of classical nova explosions (micronovae) and bear a strong resemblance to type-I X-ray bursts12-14. We exclude accretion or stellar magnetic reconnection events as their origin and suggest thermonuclear runaway events in magnetically confined accretion columns as a viable explanation.

Triggering micronovae through magnetically confined accretion flows in accreting white dwarfs

ABSTRACT Rapid bursts at optical wavelengths have been reported for several accreting white dwarfs. In these bursts, the optical luminosity can increase by up to a factor of 30 in less than an hour, before fading on time-scales of several hours, and the energy release can reach ~1039 erg (‘micronovae’). Several systems have also shown these bursts to be semirecurrent on time-scales of days to months, and the temporal profiles of these bursts strongly resemble those observed in Type-I X-ray bursts in accreting neutron stars. It has been suggested that the observed micronovae may be the result of localized thermonuclear runaways in the surface layers of accreting white dwarfs. Here we propose a model in which the magnetic confinement of accretion streams on to the accreting magnetic white dwarf may trigger localized thermonuclear runaways. The model proposed to trigger micronovae appears to favour magnetic systems with both a high white dwarf mass and a high mass-transfer rate.

A Common Envelope Jets Supernova (CEJSN) Impostor Scenario for Fast Blue Optical Transients

Abstract I propose a new scenario, the polar common envelope jets supernova (CEJSN) impostor scenario, to account for AT2018cow-like fast blue optical transients (FBOTs). The polar CEJSN impostor scenario evolves through four main phases. (1) A red supergiant (RSG) star expands to tidally interact with a neutron star (NS) companion (or a black hole). The interaction increases the RSG mass loss rate to form a circumstellar matter (CSM) halo to r ≃ 0.1 pc. (2) Shortly before the onset of a common envelope evolution (CEE) and about a year before explosion the NS accretes mass from the RSG envelope and launches jets that inflate two opposite lobes in the CSM within ≈100 au. (3) The NS-RSG system enters a CEE phase during which the system ejects most of the envelope mass in a dense equatorial outflow. (4) At the termination of the CEE the leftover envelope forms a circumbinary disk around the NS-core system. The NS accretes mass from the circumbinary disk and launches energetic jets that, when colliding with the fronts of the CSM lobes, power an FBOT event. The low mass of the jets-lobes interaction zones and their large distance, of about 100 au, from the center account for the fast transient. In the future the core collapses to form a second NS. In the far future the two NSs might merge. I suggest that FBOTs and similar fast transients are CEJSN impostors which compose a large fraction of the progenitors of NS-NS merger binaries.

Impact of disc-outflows in accreting neutron stars

Impact of disc-outflows in accreting neutron stars

Probing spectral and temporal evolution of the neutron star low-mass X-ray binary 4U 1724–30 with AstroSat

ABSTRACT We report the broad-band spectro-temporal study of the poorly studied accreting neutron star (NS) low mass X-ray binary (LMXB) 4U 1724–30 using data from Soft X-ray Telescope (SXT) and Large Area X-ray Proportional Counters (LAXPC) instruments on board AstroSat. The dim persistent LMXB source was observed with AstroSat over 4 epochs in 2017, all of which corresponded to a low-luminosity non-thermal emission dominated (hard/island) emission state with modest spectral evolution. All the X-ray broad-band spectra can be modelled by a combination of thermal emission from the NS boundary layer (BL) or NS surface and a non-thermal emission component possibly originating from the inverse Comptonization of the disc seed photons. We investigate the presence of frequency and energy-dependent variabilities to probe the origin of the disc/coronal fluctuations. We also report the detection of a Type-I X-ray burst displaying a photospheric radius expansion (PRE). During the burst, a hard X-ray shortage in the 30–80 keV energy band and the enhancement of the persistent emission reveal the burst feedback on the overall accretion process. Using the touch-down burst flux ∼4.25 × 10−8 erg s−1 cm−2, the distance of the source is estimated as ∼8.4 kpc.

A persistent ultraviolet outflow from an accreting neutron star binary transient.

All disc-accreting astrophysical objects produce powerful disc winds. In compact binaries containing neutron stars or black holes, accretion often takes place during violent outbursts. The main disc wind signatures during these eruptions are blue-shifted X-ray absorption lines, which are preferentially seen in disc-dominated 'soft states'1,2. By contrast, optical wind-formed lines have recently been detected in 'hard states', when a hot corona dominates the luminosity3. The relationship between these signatures is unknown, and no erupting system has as yet revealed wind-formed lines between the X-ray and optical bands, despite the many strong resonance transitions in this ultraviolet (UV) region4. Here we report that the transient neutron star binary Swift J1858.6-0814 exhibits wind-formed, blue-shifted absorption linesassociated with C IV, N Vand He IIin time-resolved UV spectroscopy during a luminous hard state, which we interpret as a warm, moderately ionized outflow component in this state. Simultaneously observed optical lines also display transient blue-shifted absorption. Decomposing the UV data into constant and variable components, the blue-shifted absorption is associated with the former. This implies that the outflow is not associated with the luminous flares in the data. The joint presence of UV and optical wind features reveals a multi-phase and/or spatially stratified evaporative outflow from the outer disc5. This type of persistent mass loss across all accretion states has been predicted by radiation-hydrodynamic simulations6 and helps to explain the shorter-than-expected duration of outbursts7.

MAXI and NuSTAR Observations of the Faint X-Ray Transient MAXI J1848-015 in the GLIMPSE-C01 Cluster

Abstract We present the results of Monitor of All-sky X-ray Image (MAXI) monitoring and two Nuclear Spectroscopic Telescope Array (NuSTAR) observations of the recently discovered faint X-ray transient MAXI J1848015. Analysis of the MAXI light curve shows that the source underwent a rapid flux increase beginning on 2020 December 20, followed by a rapid decrease in flux after only ∼5 days. NuSTAR observations reveal that the source transitioned from a bright soft state with unabsorbed, bolometric (0.1–100 keV) flux F = 6.9 ± 0.1 × 10−10 erg cm−2 s−1, to a low hard state with flux F = 2.85 ± 0.04 × 10−10 erg cm−2 s−1. Given a distance of 3.3 kpc, inferred via association of the source with the GLIMPSE-C01 cluster, these fluxes correspond to an Eddington fraction of the order of 10−3 for an accreting neutron star (NS) of mass M = 1.4M ⊙, or even lower for a more massive accretor. However, the source spectra exhibit strong relativistic reflection features, indicating the presence of an accretion disk that extends close to the accretor, for which we measure a high spin, a = 0.967 ± 0.013. In addition to a change in flux and spectral shape, we find evidence for other changes between the soft and hard states, including moderate disk truncation with the inner disk radius increasing from R in ≈ 3 R g to R in ≈ 8 R g, narrow Fe emission whose centroid decreases from 6.8 ± 0.1 keV to 6.3 ± 0.1 keV, and an increase in low-frequency (10−3–10−1 Hz) variability. Due to the high spin, we conclude that the source is likely to be a black hole rather than an NS, and we discuss physical interpretations of the low apparent luminosity as well as the narrow Fe emission.

Radius Constraints from Reflection Modeling of Cygnus X-2 with NuSTAR and NICER

Abstract We present a spectral analysis of NuSTAR and NICER observations of the luminous, persistently accreting neutron star (NS) low-mass X-ray binary Cygnus X-2. The data were divided into different branches that the source traces out on the Z-track of the X-ray color–color diagram; namely, the horizontal branch, the normal branch, and the vertex between the two. The X-ray continuum spectrum was modeled in two different ways that produced comparable quality fits. The spectra showed clear evidence of a reflection component in the form of a broadened Fe K line, as well as a lower-energy emission feature near 1 keV likely due to an ionized plasma located far from the innermost accretion disk. We account for the reflection spectrum with two independent models (relxillns and rdblur*rfxconv). The inferred inclination is in agreement with earlier estimates from optical observations of ellipsoidal lightcurve modeling (relxillns: i = 67° ± 4°; rdblur*rfxconv: i = 60° ± 10°). The inner disk radius remains close to the NS (R in ≤ 1.15 R ISCO) regardless of the source position along the Z-track or how the 1 keV feature is modeled. Given the optically determined NS mass of 1.71 ± 0.21 M ⊙, this corresponds to a conservative upper limit of R in ≤ 19.5 km for M = 1.92 M ⊙ or R in ≤ 15.3 km for M = 1.5 M ⊙. We compare these radius constraints to those obtained from NS gravitational wave merger events and recent NICER pulsar lightcurve modeling measurements.

Origin of young accreting neutron stars in high-mass X-ray binaries in supernova remnants

ABSTRACT Recently, several accreting neutron stars (NSs) in X-ray binary systems inside supernova remnants have been discovered. They represent a puzzle for the standard magnetorotational evolution of NSs, as their ages (≲105 yr) are much shorter than the expected duration of Ejector and Propeller stages preceding the onset of wind accretion. To explain appearance of such systems, we consider rotational evolution of NSs with early fallback accretion and asymmetry in forward/backward transitions between Ejector and Propeller stages (so-called hysteresis effect proposed by V. Shvartsman in 1970). It is shown that after a successful fallback episode with certain realistic values of the initial spin period, stellar wind properties, and magnetic field, a young NS may not enter the Ejector stage during its evolution which results in a relatively rapid initiation of accretion within the lifetime of a supernova remnant. For a standard magnetic field ∼1012 G and initial spin period ∼0.1–0.2 s accretion rate ≳1014–1015 g s−1 is enough to avoid the Ejector stage.

The Impact of Neutron Transfer Reactions on the Heating and Cooling of Accreted Neutron Star Crusts

Abstract Nuclear reactions heat and cool the crust of accreting neutron stars and need to be understood to interpret observations of X-ray bursts and long-term cooling in transiently accreting systems. It was recently suggested that previously ignored neutron transfer reactions may play a significant role in the nuclear processes. We present results from full nuclear network calculations that now include these reactions and determine their impact on crust composition, crust impurity, heating, and cooling. We find that a large number of neutron transfer reactions indeed occur and impact crust models. In particular, we identify a new type of reaction cycle that brings a pair of nuclei across the nuclear chart into equilibrium via alternating neutron capture and neutron release, interspersed with a neutron transfer. While neutron transfer reactions lead to changes in crust model predictions and need to be considered in future studies, previous conclusions concerning heating, cooling, and compositional evolution are remarkably robust.

Starquakes in millisecond pulsars and gravitational waves emission

ABSTRACT So far, only transient gravitational waves (GWs) produced by catastrophic events of extra-galactic origin have been detected. However, it is generally believed that there should be also continuous sources of GWs within our Galaxy, such as accreting neutron stars (NSs), that could in principle be detected in the next near future. In fact, in these objects, centrifugal forces can be so strong to break the NS crust (causing a starquake), thus producing a quadrupole moment responsible for the continuous emission of GWs. At equilibrium, the angular momentum gained by accretion and the one lost via GWs emission should balance each other, stopping the stellar spin-up. We hereinafter investigate the above physical picture within the framework of a Newtonian model describing compressible, non-magnetized, and self-gravitating NSs. In particular, we calculate the rotational frequency need to break the stellar crust of an accreting pulsar and we estimate the upper limit for the ellipticity due to this event. We find that the maximum starquake-induced ellipticity ranges from 10−9 to 10−5, depending on the stellar mass and its equation of state. The corresponding equilibrium frequency that we calculate is in good agreement with observations and, for all the scenarios, it is below the higher NS frequency observed of 716.36 Hz. Finally, we also discuss possible observational constraints on the ellipticity upper limit of accreting pulsars.

Sensitivity Study of Type-I X-ray Bursts to Nuclear Reaction Rates

Type-I X-ray bursts are frequently occurring thermonuclear runaways on the surface of an accreting neutron star. The different nuclear reactions that empower the cataclysmic event play a key role in accurately comparing theoretical models and observations. We investigate the effect of the uncertainties in the nuclear reactions using a ONEZONE model for a set of different compositions and accretion rates that are within the range of the standard observed burst sources. A combination of a full-reaction network and a semi-analytic model is used to obtain the conditions at the time of X-ray burst ignition via simulating the settling process of the accreted material. We then evaluate the sensitivity of the X-ray burst model by varying the proton and alpha-induced reaction rates in JINA REACLIBV2.2 within representative nuclear physics uncertainties.

Rapid Bursts of Magnetically Gated Accretion in the Intermediate Polar V1025 Cen

Abstract Magnetically gated accretion has emerged as a proposed mechanism for producing extremely short, repetitive bursts of accretion onto magnetized white dwarfs in intermediate polars (IPs), but this phenomenon has not been detected previously in a confirmed IP. We report the 27 day TESS light curve of V1025 Cen, an IP that shows a remarkable series of 12 bursts of accretion, each lasting for less than 6 hours. The extreme brevity of the bursts and their short recurrence times (∼1–3 days) are incompatible with the dwarf-nova instability, but they are natural consequences of the magnetic gating mechanism developed by Spruit and Taam to explain the Type II bursts of the accreting neutron star known as the Rapid Burster. In this model, the accretion flow piles up at the magnetospheric boundary and presses inward until it couples with the star’s magnetic field, producing an abrupt burst of accretion. After each burst, the reservoir of matter at the edge of the magnetosphere is replenished, leading to cyclical bursts of accretion. A pair of recent studies applied this instability to the suspected IPs MV Lyr and TW Pic, but the magnetic nature of these two systems has not been independently confirmed. In contrast, previous studies have unambiguously established the white dwarf in V1025 Cen to be significantly magnetized. The detection of magnetically gated bursts in a confirmed IP therefore validates the extension of the Spruit and Taam instability to magnetized white dwarfs.

Quiescent luminosities of transiently accreting neutron stars with neutrino heating due to charged pion decay

We study the quiescent luminosities of accreting neutron stars by a new mechanism as neutrino heating for additional deep crustal heating, where the neutrino heating is produced by charged pions decay from the nuclear collisions on the surface of neutron star during its active accretion. For low mass neutron star($\lesssim1.4~M_{\odot}$), as the neutrino heating is little($\lesssim1$ MeV per accreted nucleon) or there would be no neutrino heating, the quiescent luminsoity will be not affected or slightly affected. While for massive neutron star ($\gtrsim2~M_{\odot}$), the quiescent luminosity will be enhanced more obviously with neutrino heating in the range 2-6 MeV per accreted nucleon. The observations on cold neutron stars such as 1H 19605+00, SAX J1808.4-3658 can be explained with neutrino heating if a fast cooling and heavy elements surface are considered. The observations on a hot neutron star such as RX J0812.4-3114 can be explained with neutrino heating if the direct Urca process is forbidden for a massive star with light elements surface, which is different from the previous work that the hot observations should be explained with small mass neutron star and the effect of superfluidity.

The MAVERIC Survey: Variable Jet-accretion Coupling in Luminous Accreting Neutron Stars in Galactic Globular Clusters

Abstract Accreting neutron stars in low-mass X-ray binaries show outflows—and sometimes jets—in the general manner of accreting black holes. However, the quantitative link between the accretion flow (traced by X-rays) and outflows and/or jets (traced by radio emission) is much less well understood for neutron stars than for black holes, other than the general observation that neutron stars are fainter in the radio at a given X-ray luminosity. We use data from the deep MAVERIC radio continuum survey of Galactic globular clusters for a systematic radio and X-ray study of six luminous (L X > 1034 erg s−1) persistent neutron star X-ray binaries in our survey, as well as two other transient systems also captured by our data. We find that these neutron star X-ray binaries show an even larger range in radio luminosity than previously observed. In particular, in quiescence at L X ∼ 3 × 1034 erg s−1, the confirmed neutron star binary GRS 1747–312 in Terzan 6 sits near the upper envelope of the black hole radio/X-ray correlation, and the persistently accreting neutron star systems AC 211 (in M15) and X1850–087 (in NGC 6712) show unusual radio variability and luminous radio emission. We interpret AC 211 as an obscured “Z source” that is accreting at close to the Eddington limit, while the properties of X1850–087 are difficult to explain, and motivate future coordinated radio and X-ray observations. Overall, our results show that neutron stars do not follow a single relation between inflow and outflow, and confirm that their accretion dynamics are more complex than for black holes.

Urca Cooling in Neutron Star Crusts and Oceans: Effects of Nuclear Excitations.

The excited-state structure of atomic nuclei can modify nuclear processes in stellar environments. In this Letter, we study the influence of nuclear excitations on Urca cooling (repeated back-and-forth β decay and electron capture in a pair of nuclear isotopes) in the crust and ocean of neutron stars. We provide for the first time an expression for Urca process neutrino luminosity which accounts for a thermal Boltzmann distribution of excited states in both members of an Urca pair. We use our new formula with state-of-the-art nuclear structure inputs to compute neutrino luminosities of candidate Urca cooling pairs. Our nuclear inputs consist of the latest experimental data supplemented with calculations using the projected shell model. We show that, in contrast to previous results that only consider the ground states of both nuclei in the pair, our calculated neutrino luminosities for different Urca pairs vary sensitively with the environment temperature and can be radically different from those obtained in the one-transition approximation. We find that nuclear excitations can lead to an enhancement in total Urca neutrino luminosities in the accreted neutron star crust by about 5 times as compared with the previous Urca results, which is expected to cause significant observational effects.

Nonequilibrium thermodynamics of accreted neutron-star crust

We show that, in order to determine the equation of state of the inner crust of an accreting neutron star, one should minimize not the Gibbs free energy, as it is generally assumed in the literature, but a different thermodynamic potential $\Psi$, which tends to the minimum at fixed pressure and neutron chemical potential. Once this potential is specified, one can calculate the heat-release distribution in the stellar crust due to nonequilibrium nuclear reactions induced by accretion of matter onto the neutron-star surface. The results are important for adequate modeling of the accreted crust and interpretation of the observations of accreting neutron stars in low-mass X-ray binaries.

A Comprehensive X-Ray Report on AT2019wey

Abstract Here, we present MAXI, Swift, NICER, NuSTAR, and Chandra observations of the X-ray transient AT2019wey (SRGA J043520.9+552226, SRGE J043523.3+552234). From spectral and timing analyses we classify it as a Galactic low-mass X-ray binary (LMXB) with a black hole (BH) or neutron star (NS) accretor. AT2019wey stayed in the low/hard state (LHS) from 2019 December to 2020 August 21, and the hard-intermediate state (HIMS) from 2020 August 21 to 2020 November. For the first six months of the LHS, AT2019wey had a flux of ∼1 mCrab, and displayed a power-law X-ray spectrum with photon index Γ = 1.8. From 2020 June to August, it brightened to ∼20 mCrab. Spectral features characteristic of relativistic reflection became prominent. On 2020 August 21, the source left the “hard line” on the rms–intensity diagram, and transitioned from LHS to HIMS. The thermal disk component became comparable to the power-law component. A low-frequency quasi-periodic oscillation (QPO) was observed. The QPO central frequency increased as the spectrum softened. No evidence of pulsation was detected. We are not able to decisively determine the nature of the accretor (BH or NS). However, the BH option is favored by the position of this source on the Γ–L X, L radio–L X, and L opt–L X diagrams. We find the BH candidate XTE J1752−223 to be an analog of AT2019wey. Both systems display outbursts with long plateau phases in the hard states. We conclude by noting the potential of SRG in finding new members of this emerging class of low luminosity and long-duration LMXB outbursts.