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.     

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.     

Continuous frequency spectrum of the global hydromagnetic oscillations of a magnetically confined mountain on an accreting neutron star

We compute the continuous part of the ideal-magnetohydrodynamic (ideal-MHD) frequency spectrum of a polar mountain produced by magnetic burial on an accreting neutron star. Applying the formalism developed by Hellsten & Spies, extended to include gravity, we solve the singular eigenvalue problem subject to line-tying boundary conditions. This spectrum divides into an Alfvén part and a cusp part. The eigenfunctions are chirped and anharmonic with an exponential envelope, and the eigenfrequencies cover the whole spectrum above a minimum ω_low. For equilibria with accreted mass  1.2 × 10⁻⁶≲ M _a/M_⊙≲ 1.7 × 10⁻⁴  and surface magnetic fields  10¹¹≲ B _*/G ≲ 10¹³, ω_low  is approximately independent of   B _*  , and increases with M _a. The results are consistent with the Alfvén spectrum excited in numerical simulations with the zeus-mp solver. The spectrum is modified substantially by the Coriolis force in neutron stars spinning faster than ∼100 Hz. The implications for gravitational-wave searches for low-mass X-ray binaries are considered briefly.     

The radial disc structure around a magnetic neutron star: analytic and semi-analytic solutions

The radial structure of a thin accretion disc is calculated in the presence of a central dipole magnetic field aligned with the rotation axis. The problem is treated using a modified expression for the turbulent magnetic diffusion, which allows the angular momentum equation to be integrated analytically. The governing algebraic equations are solved iteratively between 1 and 10⁴ stellar radii. An analytic approximation is provided that is valid near the disruption radius at about 100 stellar radii. At that point, which is approximately 60 per cent of the Alfvén radius and typically about 30 per cent of the corotation radius, the disc becomes viscously unstable. This instability results from the fact that both radiation pressure and opacity caused by electron scattering become important. This in turn is a consequence of the magnetic field which leads to an enhanced temperature in the inner parts. This is because the magnetic field gives rise to a strongly enhanced vertically integrated viscosity, so that the viscous torque can balance the magnetic torque.     

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.     

Burial of the polar magnetic field of an accreting neutron star – II. Hydromagnetic stability of axisymmetric equilibria

The theory of polar magnetic burial in accreting neutron stars predicts that a mountain of accreted material accumulates at the magnetic poles of the star, and that, as the mountain spreads equatorward, it is confined by, and compresses, the equatorial magnetic field. Here, we extend previous, axisymmetric, Grad–Shafranov calculations of the hydromagnetic structure of a magnetic mountain up to accreted masses as high as   M _a= 6 × 10⁻⁴ M_⊙  , by importing the output from previous calculations (which were limited by numerical problems and the formation of closed bubbles to   M _a < 10⁻⁴ M_⊙  ) into the time-dependent, ideal-magnetohydrodynamic code zeus-3d and loading additional mass on to the star dynamically. The rise of buoyant magnetic bubbles through the accreted layer is observed in these experiments. We also investigate the stability of the resulting hydromagnetic equilibria by perturbing them in zeus-3d . Surprisingly, it is observed that the equilibria are marginally stable for all   M _a≤ 6 × 10⁻⁴ M_⊙  ; the mountain oscillates persistently when perturbed, in a combination of Alfvén and acoustic modes, without appreciable damping or growth, and is therefore not disrupted (apart from a transient Parker instability initially, which expels <1 per cent of the mass and magnetic flux).     

Constraining the neutron star equation of state using XMM‐Newton

We have identified three possible ways in which future XMM‐Newton observations can provide significant constraints on the equation of state of neutron stars. First, using a long observation of the neutron star X‐ray transient Cen X‐4 in quiescence one can use the RGS spectrum to constrain the interstellar extinction to the source. This removes this parameter from the X‐ray spectral fitting of the pn and MOS spectra and allows us to investigate whether the variability observed in the quiescent X‐ray spectrum of this source is due to variations in the soft thermal spectral component or variations in the power law spectral component coupled with variations in N_H. This will test whether the soft thermal spectral component can indeed be due to the hot thermal glow of the neutron star. Potentially such an observation could also reveal redshifted spectral lines from the neutron star surface. Second, XMM‐Newton observations of radius expansion type I Xray bursts might reveal redshifted absorption lines from the surface of the neutron star. Third, XMM‐Newton observations of eclipsing quiescent low‐mass X‐ray binaries provide the eclipse duration. With this the system inclination can be determined accurately. The inclination determined from the X‐ray eclipse duration in quiescence, the rotational velocity of the companion star and the semi‐amplitude of the radial velocity curve determined through optical spectroscopy, yield the neutron star mass. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

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.     

Evidence for a jet contribution to the optical/infrared light of neutron star X-ray binaries

Optical/near-infrared (optical/NIR, OIR) light from low-mass neutron star X-ray binaries (NSXBs) in outburst is traditionally thought to be thermal emission from the accretion disc. Here we present a comprehensive collection of quasi-simultaneous OIR and X-ray data from 19 low magnetic field NSXBs, including new observations of three sources: 4U 0614+09, LMC X−2 and GX 349+2. The average radio–OIR spectrum for NSXBs is  α≈+ 0.2  (where   L _ν∝ν^α  ) at least at high luminosities when the radio jet is detected. This is comparable to, but slightly more inverted than the  α≈ 0.0  found for black hole X-ray binaries. The OIR spectra and relations between OIR and X-ray fluxes are compared to those expected if the OIR emission is dominated by thermal emission from an X-ray or viscously heated disc, or synchrotron emission from the inner regions of the jets. We find that thermal emission due to X-ray reprocessing can explain all the data except at high luminosities for some NSXBs, namely, the atolls and millisecond X-ray pulsars. Optically thin synchrotron emission from the jets (with an observed OIR spectral index of  α_thin < 0  ) dominate the NIR light above     and the optical above     in these systems. For NSXB Z-sources, the OIR observations can be explained by X-ray reprocessing alone, although synchrotron emission may make a low-level contribution to the NIR, and could dominate the OIR in one or two cases.     

Interaction and ablation of fall-back discs in isolated neutron stars

An analysis of ablation processes is made for a fall-back disc with inner and outer radii external to the neutron-star light cylinder. The calculated ablation rate leads, with certain other assumptions, to a simple expression relating the inner radius and mean mass per unit area of any long-lived fall-back disc. Expressions for the torque components generated by interaction with the pulsar wind are obtained. It is not impossible that these could be responsible for small observable variations in pulse shape and spin-down rate but they are unlikely to be the source of the periodic changes seen in several pulsars.     

Influence of the r-mode instability on hypercritically accreting neutron stars

We have investigated the influence of the r-mode instability on hypercritically accreting neutron stars in close binary systems during their common envelope phases, based on the scenario proposed by Brown et al. On the one hand, neutron stars are heated by the accreted matter at the stellar surface, but on the other hand they are also cooled down by the neutrino radiation. At the same time, the accreted matter transports its angular momentum and mass to the star. We have studied the evolution of the stellar mass, temperature and rotational frequency.
The gravitational-wave-driven instability of the r-mode oscillation strongly suppresses spinning up of the star, the final rotational frequency of which is well below the mass-shedding limit, in fact typically as low as 10 per cent of that of the mass-shedding state. On a very short time-scale the rotational frequency tends to approach a certain constant value and saturates there, as long as the amount of accreted mass does not exceed a certain limit to collapse to a black hole. This implies that a similar mechanism of gravitational radiation to that in the so-called 'Wagoner star' may work in this process. The star is spun up by accretion until the angular momentum loss by gravitational radiation balances the accretion torque. The time-integrated dimensionless strain of the radiated gravitational wave may be large enough to be detectable by gravitational wave detectors such as LIGO II.     

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.