Nuclear reactions on accreting neutron star surface

The production of X-rays and gamma-rays in bursts is believed to be due to the rapid burning of matter accreted onto a neutron star surface from its companion, most probably a giant star. The accreted matter consists mainly of hydrogen and helium and a very small amount of heavy elements. Due to the infall of matter, the temperature at the bottom layers is raised to a value of the order of 10⁸ K. The neutron star surface density is>10⁷ g cm^–3. As hydrogen burning is a slow process under any temperature and density conditions, we consider the helium-burning reactions as the source of gamma-rays in the neutron star surface. Under high-density conditions the ordinary laboratory reaction rates should become modified. At high-density conditions, the strong screening effect due to the polarising cloud of electrons around the ions become important and enhances the reaction rates considerably. The helium-burning reactions are calculated under such conditions. The abundances of helium-burning products such as¹²C, 1¹⁶O, and²⁰Ne, etc., are computed. Under high-density and temperature conditions carbon is found to be more abundant than oxygen. Neon is completely absent in almost all the relevant physical conditions in which a strong screening effect is operative. It is suggested that explosive burning of accreted helium of 10^–13 M will account for the observed energy of gamm-ray burst.     

Boundary layer on the surface of a neutron star

In an attempt to model the accretion on to a neutron star in low-mass X-ray binaries, we present 2D hydrodynamical models of the gas flow in close vicinity of the stellar surface. First, we consider a gas pressure-dominated case, assuming that the star is non-rotating. For the stellar mass we take   M _star= 1.4 × 10⁻² M_⊙  and for the gas temperature   T = 5 × 10⁶ K  . Our results are qualitatively different in the case of a realistic neutron star mass and a realistic gas temperature of T ≃ 10⁸ K, when the radiation pressure dominates. We show that to get the stationary solution in a latter case, the star most probably has to rotate with the considerable velocity.     

Dynamics of the flows accreting onto a magnetized neutron star

We investigate the unsteady column accretion of material at a rate \(10^{15} g s^{ - 1} \leqslant \dot M \leqslant 10^{16} g s^{ - 1}\) onto the surface of a magnetized neutron star using a modified first-order Godunov method with splitting. We study the dynamics of the formation and evolution of a shock in an accretion column near the surface of a star with a magnetic field 5×10¹¹≤B≤10¹³ G. An effective transformation of the accretion flow energy into cyclotron radiation is shown to be possible for unsteady accretion with a collisionless shock whose front executes damped oscillations. The collisionless deceleration of the accreting material admits the conservation of a fraction of the heavy nuclei that have not been destroyed in spallation reactions. The fraction of the CNO nuclei that reach the stellar atmosphere is shown to depend on the magnetic field strength of the star.     

Spectra of the spreading layers on the neutron star surface and constraints on the neutron star equation of state

Spectra of the spreading layers on the neutron star surface are calculated on the basis of the Inogamov–Sunyaev model taking into account general relativity correction to the surface gravity and considering various chemical composition of the accreting matter. Local (at a given latitude) spectra are similar to the X-ray burst spectra and are described by a diluted blackbody. Total spreading layer spectra are integrated accounting for the light bending, gravitational redshift and the relativistic Doppler effect and aberration. They depend slightly on the inclination angle and on the luminosity. These spectra also can be fitted by a diluted blackbody with the colour temperature depending mainly on a neutron star compactness. Owing to the fact that the flux from the spreading layer is close to the critical Eddington, we can put constraints on a neutron star radius without the need to know precisely the emitting region area or the distance to the source. The boundary layer spectra observed in the luminous low-mass X-ray binaries, and described by a blackbody of colour temperature   T _c= 2.4 ± 0.1 keV  , restrict the neutron star radii to   R = 14.8 ± 1.5 km  (for a  1.4-M_⊙  star and solar composition of the accreting matter), which corresponds to the hard equation of state.     

Comptonization and QPO origins in accreting neutron star systems

We develop a simple, time-dependent Comptonization model to probe the origins of spectral variability in accreting neutron star systems. In the model, soft 'seed photons' are injected into a corona of hot electrons, where they are Compton upscattered before escaping as hard X-rays. The model describes how the hard X-ray spectrum varies when the properties of either the soft photon source or the Comptonizing medium undergo small oscillations. Observations of the resulting spectral modulations can determine whether the variability is due to (i) oscillations in the injection of seed photons, (ii) oscillations in the coronal electron density, or (iii) oscillations in the coronal energy dissipation rate. Identifying the origin of spectral variability should help clarify how the corona operates and its relation to the accretion disc. It will also help in finding the mechanisms underlying the various quasi-periodic oscillations (QPOs) observed in the X-ray outputs of many accreting neutron star and black hole systems. As a sample application of our model, we analyse a kilohertz QPO observed in the atoll source 4U 1608–52. We find that the QPO is driven predominantly by an oscillation in the electron density of the Comptonizing gas.     

Biases for neutron star mass, radius and distance measurements from Eddington-limited X-ray bursts

Eddington-limited X-ray bursts from neutron stars can be used in conjunction with other spectroscopic observations to measure neutron star masses, radii and distances. In order to quantify some of the uncertainties in the determination of the Eddington limit, we analysed a large sample of photospheric radius-expansion thermonuclear bursts observed with the Rossi X-ray Timing Explorer . We identified the instant at which the expanded photosphere 'touches down' back on to the surface of the neutron star and compared the corresponding touchdown flux to the peak flux of each burst. We found that for the majority of sources, the ratio of these fluxes is smaller than ≃1.6, which is the maximum value expected from the changing gravitational redshift during the radius expansion episodes (for a  2 M_⊙  neutron star). The only sources for which this ratio is larger than ≃1.6 are high-inclination sources that include dippers and Cyg X-2. We discuss two possible geometric interpretations of this effect and show that the inferred masses and radii of neutron stars are not affected by this bias. On the other hand, systematic uncertainties as large as ∼50 per cent may be introduced to the distance determination.     

On hard X-ray spectra of accreting neutron stars

Formation of the spectra of X-ray pulsars and gamma bursters is investigated. Interpretation of a hard X-ray spectrum of pulsars containing cyclotron lines is feasible on the basis of an isothermal model of a polar spot heated due to accretion to a neutron star. It has been ascertained that in the regions responsible for the formation of continuum radiation and lines the mode polarization is determined by a magnetized vacuum rather than by a plasma. Bearing this in mind, the influence of the magnetic field of a star on the wide wings of the cyclotron line and on its depth is discussed. The part played by the accreting column in the case of strong accretion (10¹⁹ el cm^–3) needed for long sustaining of the high level of X-rays from a neutron star-pulsar is studied. There occur the gaps in spectrum at frequencies close to the electron gyro-frequency and its harmonics due to the screening of the hot spot by the opaque gyro-resonant layers located within the accreting column. These gaps ensure the formation of cyclotron lines in absorption irrespective of the presence of such lines in the X-ray spectrum of a polar hot spot.The spectra of gamma-bursters recorded by Venus 11 and Venus 12 are interpreted in terms of a two-layer model of a polar hot spot. The estimates are given of the distance to some of the bursters, of the emission measure from a high-temperature layer responsible for continuum radiation and of the dispersion measure of a colder layer forming cyclotron lines in absorption. It is noted that the action of an accreting column leads generally to the radiation depression at frequencies below cyclotron lines. By the observed depression for one of the bursters the electron density of near-star accreting plasma during the burst has been directly estimated (4×10^–14 el cm^–3). Possible appearance of false cyclotron lines associated with cyclotron scattering in accreting column has been revealed.The problem of measuring the magnetic fields of neutron stars taking account of the gravitational redshift and the quantum recoil effect in emission and in absorption is discussed. Possibilty for a more precise measurement of the magnetic fields of those bursters whose spectrum contains both a cyclotron and an annihilation lines is noted.     

Resistive relaxation of a magnetically confined mountain on an accreting neutron star

Three-dimensional numerical magnetohydrodynamic (MHD) simulations are performed to investigate how a magnetically confined mountain on an accreting neutron star relaxes resistively. No evidence is found for non-ideal MHD instabilities on a short time-scale, such as the resistive ballooning mode or the tearing mode. Instead, the mountain relaxes gradually as matter is transported across magnetic surfaces on the diffusion time-scale, which evaluates to  τ_I∼ 10⁵–10⁸ yr  (depending on the conductivity of the neutron star crust) for an accreted mass of   M _a= 1.2 × 10⁻⁴ M_⊙  . The magnetic dipole moment simultaneously re-emerges as the screening currents dissipate over  τ_I  . For non-axisymmetric mountains, ohmic dissipation tends to restore axisymmetry by magnetic reconnection at a filamentary neutral sheet in the equatorial plane. Ideal-MHD oscillations on the Alfvén time-scale, which can be excited by external influences, such as variations in the accretion torque, compress the magnetic field and hence decrease  τ_I  by one order of magnitude relative to its standard value (as computed for the static configuration). The implications of long-lived mountains for gravitational wave emission from low-mass X-ray binaries are briefly explored.     

Mechanism of thermonuclear burning propagation in a helium layer on a neutron star surface: A simplified adiabatic model

Some thermonuclear X-ray bursters exhibit a high-frequency (about 300 Hz or more) brightness modulation at the rising phase of some bursts. These oscillations are explained by inhomogeneous heating of the surface layer on a rapidly rotating neutron star due to the finite propagation speed of thermonuclear burning. We suggest and substantiate a mechanism of this propagation that is consistent with experimental data. Initially, thermonuclear ignition occurs in a small region of the neutron star surface layer. The burning products rapidly rise and spread in the upper atmospheric layers due to turbulent convection. The accumulation of additional matter leads to matter compression and ignition at the bottom of the layer. This determines the propagation of the burning front. To substantiate this mechanism, we use the simplifying assumptions about a helium composition of the neutron star atmosphere and its initial adiabatic structure with a density of 1.75 × 10⁸ g cm⁻³ at the bottom. 2D numerical simulations have been performed using a modified particle method in the adiabatic approximation.     

The de-excited energy of electron capture in accreting neutron star crusts

When a daughter nucleus produced by electron capture takes part in a level transition from an excited state to its ground state in accreting neutron star crusts, thermal energy will be released and heat the crust, increasing crust temperature and changing subsequent carbon ignition conditions. Previous studies show that the theoretical carbon ignition depth is deeper than the value inferred from observations because the thermal energy is not sufficient. In this paper, we present the de-excited energy from electron capture of rp-process ash before carbon ignition, especially for the initial evolution stage of rp-process ash, by using a level-to-level transition method. We find the theoretical column density of carbon ignition in the resulting superbursts and compare it with observations. The calculation of the electron capture process is based on a more reliable level-to-level transition, adopting new data from experiments or theoretical models(e.g., large-scale shell model and proton-neutron quasi-particle random phase approximation). The new carbon ignition depth is estimated by fitting from previous results of a nuclear reaction network. Our results show the average de-excited energy from electron capture before carbon ignition is ~0.026 Me V/u, which is significantly larger than the previous results. This energy is beneficial for enhancing the crust's temperature and decreasing the carbon ignition depth of superbursts.     

Statified model of a neutron star

Evolutionary calculations based on realistic equations of state indicate the stratified nature of the distribution of hadron matter in the interiors of neutron stars. In the proposed model, the stratified structure of a neutron star is treated as a rigid inert core surrounded by a dynamical layer. The physical basis for the model is the concept of the stellar matter of the peripheral envelope as an elastic Fermi continuum, the motions of which are described by the equations of nuclear elastodynamics, proposed in the macroscopic theory of collective processes in laboratory nuclear physics. It is shown that the vibrational dynamics of a neutron star is characterized by two branches of gravitational—elastic, spheroidal (s-mode) and torsional (t-mode) nonradial eigenvibrations. Estimates obtained for the periods of global, gravitational nonradial modes suggest that variations in the intensity of micropulses observed in the millisecond range of the spectra of C-pulsars may be ascribed to these vibrations. The proposed two-component model of a neutron star enables one to consider a glitch in a pulsar’s radio emission as a starquake due to the passage of the companion through periastron of the binary system. Translated from Astrofizika, Vol. 42, No. 2, pp. 235–252, April–June, 1999     

The formation of X-ray radiation in a boundary layer during disk accretion onto a neutron star

We present computed radiation spectra for the boundary layer (BL) of the accretion disk that is formed near the surface of a neutron star. Both free-free processes and Comptonization were taken into account. Our computations are based on the hydrodynamic solution obtained by Popham and Sunyaev (2001) for the BL structure. The computed spectra are highly diluted compared to the Planck spectra of the same surface temperature. They are complex in shape; in particular, an intense Wien emission component is formed in their high-energy region at high accretion rates. In general, the computed spectra are harder than those observed in actual X-ray sources. This is the result of a very high temperature found by Popham and Sunyaev (2001) for the BL. We show that such temperatures could result from an oversimplified treatment of radiative transfer in their paper, which completely ignored the frequency dependence of the matter opacity and radiation intensity. Our computations indicate that at moderate accretion rates, a proper treatment of radiative transfer with allowance for Comptonization leads to appreciably lower plasma temperatures and to softer radiation spectra.     

Nuclear deflagration on the surface of accreting neutron stars

We analyze the structure of a nuclear deflagration front in the crust of accreting neutron stars. Models of quasi-stationary deflagration fronts are calculated and subsequently evolved in time and space in order to check their stability. Unlike white dwarfs, where the velocity of aninwards propagating combustion front is governed mainly by energy losses to the gravitational field, the structure of a deflagration front in neutron stars is determined essentially by the heat fluxes into the stellar core.     

Analysing the atolls: X-ray spectral transitions of accreting neutron stars

We systematically analyse all the available X-ray spectra of disc accreting neutron stars (atolls and millisecond pulsars) from the RXTE data base. We show that while all these have similar spectral evolution as a function of mass accretion rate, there are also subtle differences. There are two different types of hard/soft transition, those where the spectrum softens at all energies, leading to a diagonal track on a colour–colour diagram, and those where only the higher energy spectrum softens, giving a vertical track. The luminosity at which the transition occurs is correlated with this spectral behaviour, with the vertical transition at   L / L _Edd∼ 0.02  while the diagonal one is at ∼0.1. Superimposed on this is the well-known hysteresis effect, but we show that classic, large-scale hysteresis occurs only in the outbursting sources, indicating that its origin is in the dramatic rate of change of mass accretion rate during the disc instability. We show that the long-term mass accretion rate correlates with the transition behaviour, and speculate that this is due to the magnetic field being able to emerge from the neutron star surface for low average mass accretion rates. While this is not strong enough to collimate the flow except in the millisecond pulsars, its presence may affect the inner accretion flow by changing the properties of the jet.     

Type-I X-ray bursts from voids: tracers of weakly accreting bursters?

We obtained long-term (10–20 years) light curves for seven X-ray bursters. These sources exhibited no prolonged episodes of luminosities exceeding several percent of the Eddington luminosity over the entire observing period. For four sources, we found upper limits for the luminosity of over 5 years. These limits proved to be below 10³⁶ erg s^?1. We estimated the total number of such sources in our Galaxy.     

Quasi-periodic oscillations in low-mass X-ray binaries and constraints on the equation of state of neutron star matter

Recently discovered quasi-periodic oscillations in the X-ray brightness of low-mass X-ray binaries are used to derive constraints on the mass of the neutron star component and the equation of state of neutron star matter. The observations are compared with models of rapidly rotating neutron stars which are calculated by means of an exact numerical method in full relativity. For the equations of state we select a broad collection of models representing different assumptions about the many-body structure and the complexity of the composition of superdense matter. The mass constraints differ from their values in the approximate treatment by ∼10 per cent. Under the assumption that the maximum frequency of the quasi-periodic oscillations originates from the innermost stable orbit, the mass of the neutron star is in the range M ∼1.92–2.25 M_⊙. The quasi-periodic oscillation in the Atoll-source 4U 1820−30 in particular is only consistent with equations of state that are rather stiff at high densities, which is explainable, so far, only with pure nucleonic/leptonic composition. This interpretation contradicts the hypothesis that the protoneutron star formed in SN 1987A collapsed to a black hole, since this would demand a maximum neutron star mass below 1.6 M_⊙. The recently suggested identification of quasi-periodic oscillations with frequencies of about 10 Hz with the Lense–Thirring precession of the accretion disc is found to be inconsistent with the models studied in this work, unless it is assumed that the first overtone of the precession is observed.     

Torus formation in neutron star mergers and well-localized short gamma-ray bursts

Merging neutron stars (NSs) are hot candidates for the still enigmatic sources of short gamma-ray bursts (GRBs). If the central engines of the huge energy release are accreting relic black holes (BHs) of such mergers, it is important to understand how the properties of the BH–torus systems, in particular disc masses and mass and rotation rate of the compact remnant, are linked to the characterizing parameters of the NS binaries. For this purpose, we present relativistic smoothed particle hydrodynamic simulations with conformally flat approximation of the Einstein field equations and a physical, non-zero temperature equation of state. Thick disc formation is highlighted as a dynamical process caused by angular momentum transfer through tidal torques during the merging process of asymmetric systems or in the rapidly spinning triaxial post-merger object. Our simulations support the possibility that the first well-localized short and hard GRBs 050509b, 050709, 050724, 050813 have originated from NS merger events and are powered by neutrino-antineutrino annihilation around a relic BH–torus system. Using model parameters based on this assumption, we show that the measured GRB energies and durations lead to estimates for the accreted masses and BH mass accretion rates which are compatible with theoretical expectations. In particular, the low-energy output and short duration of GRB 050509b set a very strict upper limit of less than 100 ms for the time interval after the merging until the merger remnant has collapsed to a BH, leaving an accretion torus with a small mass of only  ∼0.01 M_⊙  . This favours a (nearly) symmetric NS+NS binary with a typical mass as progenitor system.     

A study of type I X-ray bursts from an NS accreting pure helium

Using the Modules for Experiments in Stellar Astrophysics(MESA) code, we investigate Type I X-ray bursts(XRBs) produced by neutron stars(NSs) accreting pure helium, which are called intermediate XRBs in observations. We simulate 21 models for intermediate XRBs with various mass-accretion rates(■) from 2.5 × 10~(-8) to 5 × 10~(-10) M_⊙yr~(-1). Compared with normal XRBs, in which the NS accretes matter with solar metallicity, intermediate XRBs have higher luminosity and longer recurrence time,which are essentially consistent with observations. We find that the recurrence time of intermediate XRBs is proportional to ■~(-2.0).