No Arabic abstract
Across black hole (BH) and neutron star (NS) low-mass X-ray binaries (LMXBs), there appears to be some correlation between certain high- and low-frequency quasi-periodic oscillations (QPOs). In a previous paper, we showed that for BH LMXBs, this could be explained by the simultaneous oscillation and precession of a hot, thick, torus-like corona. In the current work, we extend this idea to NS LMXBs by associating the horizontal branch oscillations (HBO) with precession and the upper-kiloHertz (ukHz) QPO with vertical epicyclic motion. For the Atoll source 4U 1608-52, the model can match many distinct, simultaneous observations of the HBO and ukHz QPO by varying the inner and outer radius of the torus, while maintaining fixed values for the mass (M_{NS}) and spin (a_*) of the neutron star. The best fit values are M_{NS} = 1.38 pm 0.03 M_odot and a_* = 0.325 pm 0.005. By combining these constraints with the measured spin frequency, we are able to obtain an estimate for the moment of inertia of I_{NS} = 1.40 pm 0.02 times 10^{45} g cm^2, which places constraints on the equation of state. The model is unable to fit the lower-kHz QPO, but evidence suggests that QPO may be associated with the boundary layer between the accretion flow and the neutron star surface, which is not treated in this work.
The high-frequency quasi-periodic oscillations (HF QPOs) that appear in the X-ray fluxes of low-mass X-ray binaries remain an unexplained phenomenon. Among other ideas, it has been suggested that a non-linear resonance between two oscillation modes in an accretion disc orbiting either a black hole or a neutron star plays a role in exciting the observed modulation. Several possible resonances have been discussed. A particular model assumes resonances in which the disc-oscillation modes have the eigenfrequencies equal to the radial and vertical epicyclic frequencies of geodesic orbital motion. This model has been discussed for black hole microquasar sources as well as for a group of neutron star sources. Assuming several neutron (strange) star equations of state and Hartle-Thorne geometry of rotating stars, we briefly compare the frequencies expected from the model to those observed. Our comparison implies that the inferred neutron star radius RNS is larger than the related radius of the marginally stable circular orbit rms for nuclear matter equations of state and spin frequencies up to 800Hz. For the same range of spin and a strange star (MIT) equation of state, the inferrred radius RNS is roughly equal to rms. The Paczynski modulation mechanism considered within the model requires that RNS < rms. However, we find this condition to be fulfilled only for the strange matter equation of state, masses below one solar mass, and spin frequencies above 800Hz. This result most likely falsifies the postulation of the neutron star 3:2 resonant eigenfrequencies being equal to the frequencies of geodesic radial and vertical epicyclic modes. We suggest that the 3:2 epicyclic modes could stay among the possible choices only if a fairly non-geodesic accretion flow is assumed, or if a different modulation mechanism operates.
It is shown that depletion of the magnetic field pressure in a quaking neutron star undergoing Lorentz-force-driven torsional seismic vibrations about axis of its dipole magnetic moment is accompanied by the loss of vibration energy of the star that causes its vibration period to lengthen at a rate proportional to the rate of magnetic field decay. Highlighted is the magnetic-field-decay induced conversion of the energy of differentially rotational Alfven vibrations into the energy of oscillating magneto-dipole radiation. A set of representative examples of magnetic field decay illustrating the vibration energy powered emission with elongating periods produced by quaking neutron star are considered and discussed in the context of theory of magnetars.
3D MHD simulation of accretion onto neutron stars have shown in the last few years that the footprint (hotspot) of the accretion flow changes with time. Two different kinds of accretion, namely the funnel flow and the equatorial accretion produced by instabilities at the inner disk, produce different kinds of motion of the hotspot. The funnel flow produces hotspots that move around the magnetic pole, while instabilities produce other hotspots that appear randomly and move along the equator or slightly above. The angular velocities of the two hotspots are different, the equatorial one being higher and both close to the Keplerian velocity in the inner region. Modeling of the lightcurves of these hotspots with Monte Carlo simulations show that the signatures produced in power specra by them, if observed, are QPOs plus low frequency components. Their frequencies, general behavior and features describe correctly most of the properties of kHz QPOs, if we assume the funnel flow hotspots as the origin of the lower kHz QPO and instabilities as the origin of the upper kHz QPO.
We describe the model of surface emission from a rapidly rotating neutron star that is applied to Neutron Star Interior Composition Explorer X-ray data of millisecond pulsars in order to statistically constrain the neutron star mass-radius relation and dense matter equation of state. To ensure that the associated calculations are both accurate and precise, we conduct an extensive suite of verification tests between our numerical codes for both the Schwarzschild + Doppler and Oblate Schwarzschild approximations, and compare both approximations against exact numerical calculations. We find superb agreement between the code outputs, as well as in comparison against a set of analytical and semi-analytical calculations, which combined with their speed, demonstrates that the codes are well-suited for large-scale statistical sampling applications. A set of verified, high-precision reference synthetic pulse profiles is provided to the community to facilitate testing of other independently developed codes.
We use two-dimensional, general relativistic, viscous, radiation hydrodynamic simulations to study the impact of a Type I X-ray burst on a hot and geometrically thick accretion disk surrounding an unmagnetized, non-rotating neutron star. The disk is initially consistent with a system in its low/hard spectral state, and is subject to a burst which rises to a peak luminosity of $10^{38}$ erg s$^{-1}$ in $2.05$ s. At the peak of the burst, the temperature of the disk has dropped by more than three orders of magnitude and its scale height has gone down by more than one order of magnitude. The simulations show that these effects predominantly happen due to Compton cooling of the hot plasma, and clearly illustrate the potential cooling effects of bursts on accretion disk coronae. In addition, we demonstrate the presence of Poynting-Robertson drag, though it only enhances the mass accretion rate onto the neutron star by a factor of $sim 3$-$4$ compared to a simulation with no burst. Simulations such as these are important for building a general understanding of the response of an accretion disk to an intense X-ray impulse, which, in turn, will be crucial for deciphering burst spectra. Detailed analysis of such spectra offers the potential to measure neutron star radii, and hence constrain the neutron star equation of state, but only if the contributions coming from the impacted disk and its associated corona can be understood.