Einstein's general relativity predicts that orbital motion of accreted gas approaching a neutron star (NS) in a NS low-mass X-ray binary (LMXB) system occurs on a time scale proportional to the NS mass. Radiation of the gas accounts for most of the observed LMXBs variability. In more than a dozen of sources twin-peak quasi-periodic oscillations (QPOs) have been observed. Inspired by the expected proportionality between periods of orbital motion and NS mass we present a straightforward comparison among these sources. We investigate relations between QPO periods and their ratios and identify characteristic time scales of QPOs associated to individual sources. These timescales are likely determined by the relative mass of each NS. We show that the characteristic time scale of the millisecond pulsar XTE J1807.4-294 is longer than for most other NS LMXBs. Consequently, models of QPOs that consider geodesic orbital frequencies imply that the X-ray pulsars' mass has to be about 50 % higher than the average mass of other sources. Consideration of other X-ray pulsars indicates that the exceptionality of XTE J1807.4-294 cannot be related to NS magnetic field in any simple manner. We suggest that QPOs observed in this source can help to discriminate between the proposed versions of the NS equation of state.

The external Hartle-Thorne geometry, which describes the space-time outside a slowly-rotating compact star, is characterized by the gravitational mass $M$, angular momentum $J$ and quadrupole moment $Q$ of the star and gives a convenient description which, for the rotation frequencies of more than 95 % of known pulsars, is sufficiently accurate for most purposes. We focus here on the motion of particles in these space-times, presenting a detailed systematic analysis of the frequency properties of radial and vertical epicyclic motion and of orbital motion. Our investigation is motivated by X-ray observations of binary systems containing a rotating neutron star which is accreting matter from its binary companion. In these systems, twin high-frequency quasi-periodic oscillations are sometimes observed with a frequency ratio approaching $3:2$ or $5:4$ and these may be explained by models involving the orbital and epicyclic frequencies of quasi-circular geodesic motion. In our analysis, we use realistic equations of state for the stellar matter and proceed in a self-consistent way, following the Hartle-Thorne approach in calculating both the corresponding values of $Q$, $M$ and $J$ for the stellar model and the properties of the surrounding spacetime. Our results are then applied to a range of geodetical models for QPOs. A key feature of our study is that it implements the recently-discovered universal relations among neutron star parameters so that the results can be directly used for models with different masses $M$, radii $R$ and rotational frequencies $f_\mathrm{rot}$.

In this paper we describe the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to accretion flows in the strong field gravity regime around both stellar-mass and supermassive black-holes. eXTP has the unique capability of using advanced 'spectral-timing-polarimetry' techniques to analyze the rapid variations with three orthogonal diagnostics of the flow and its geometry, yielding unprecedented insight into the inner accreting regions, the effects of strong field gravity on the material within them and the powerful outflows which are driven by the accretion process.

In this White Paper we present the potential of the Enhanced X-ray Timing and Polarimetry (eXTP) mission for determining the nature of dense matter; neutron star cores host an extreme density regime which cannot be replicated in a terrestrial laboratory. The tightest statistical constraints on the dense matter equation of state will come from pulse profile modelling of accretion-powered pulsars, burst oscillation sources, and rotation-powered pulsars. Additional constraints will derive from spin measurements, burst spectra, and properties of the accretion flows in the vicinity of the neutron star. Under development by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

Twin-peak quasi-periodic oscillations (QPOs) are observed in several low-mass X-ray binary systems containing neutron stars (NSs). Timing analysis of X-ray fluxes of more than dozen of such systems reveals remarkable correlations between the frequencies of two characteristic peaks present in the power density spectra. The individual correlations clearly differ, but they roughly follow a common individual pattern. High values of measured QPO frequencies and strong modulation of the X-ray flux both suggest that the observed correlations are connected to orbital motion in the innermost part of an accretion disc. Several attempts to model these correlations with simple geodesic orbital models or phenomenological relations have failed in the past. We find and explore a surprisingly simple analytic relation that reproduces individual correlations for a group of several sources through a single parameter. When an additional free parameter is considered within our relation, it well reproduces the data of a large group of 14 sources. The very existence and form of this simple relation supports the hypothesis of the orbital origin of QPOs and provides the key for further development of QPO models. We discuss a possible physical interpretation of our relation's parameters and their links to concrete QPO models.

Twin-peak quasi-periodic oscillations (QPOs) are observed in the X-ray power-density spectra of several accreting low-mass neutron star (NS) binaries. In our previous work we have considered several QPO models. We have identified and explored mass-angular-momentum relations implied by individual QPO models for the atoll source 4U 1636-53. In this paper we extend our study and confront QPO models with various NS equations of state (EoS). We start with simplified calculations assuming Kerr background geometry and then present results of detailed calculations considering the influence of NS quadrupole moment (related to rotationally induced NS oblateness) assuming Hartle-Thorne spacetimes. We show that the application of concrete EoS together with a particular QPO model yields a specific mass-angular-momentum relation. However, we demonstrate that the degeneracy in mass and angular momentum can be removed when the NS spin frequency inferred from the X-ray burst observations is considered. We inspect a large set of EoS and discuss their compatibility with the considered QPO models. We conclude that when the NS spin frequency in 4U 1636-53 is close to 580Hz we can exclude 51 from 90 of the considered combinations of EoS and QPO models. We also discuss additional restrictions that may exclude even more combinations. Namely, there are 13 EOS compatible with the observed twin peak QPOs and the relativistic precession model. However, when considering the low frequency QPOs and Lense-Thirring precession, only 5 EOS are compatible with the model.

Isolated pulsars are rotating neutron stars with accurately measured angular velocities $\Omega$, and their time derivatives that show unambiguously that the pulsars are slowing down. The commonly accepted view is that it arises through emission of magnetic dipole radiation (MDR) from a rotating magnetized body. The calculated energy loss by a rotating pulsar with a constant moment of inertia is assumed proportional to a model dependent power of $\Omega$. This relation leads to the power law $\dot{\Omega}$ = -K $\Omega^{\rm n}$ where $n$ is called the braking index. The MDR model predicts $n$ exactly equal to 3. Selected observations of isolated pulsars provide rather precise values of $n$, individually accurate to a few percent or better, in the range 1$ <$ n $ < $ 2.8, which is consistently less than the predictions of the MDR model. In spite of an extensive investigation of various modifications of the MDR model, no satisfactory explanation of observation has been found yet. The aim of this work is to determine the deviation of the value of $n$ from the canonical $n = 3$ for a star with a frequency dependent moment of inertia in the region of frequencies from zero to the Kepler velocity (onset of mass shedding by a rotating deformed star), in the macroscopic MDR model. For the first time, we use microscopic realistic equations of state (EoS) of the star to determine its behavior and structure. In addition, we examine the effects of the baryonic mass M$_{\rm B}$ of the star, and possible core superfluidity, on the value of the braking index within the MDR model.

Serious theoretical effort has been devoted to explain the observed frequencies of twin-peak quasi-periodic oscillations (HF QPOs) observed in low-mass X-ray neutron star binaries. Here we propose a new model of HF QPOs. Within its framework we consider an oscillating torus with cusp that changes location $r_0$ of its centre around radii very close to innermost stable circular orbit. The observed variability is assigned to global modes of accreted fluid motion that may give strong modulation of both accretion disc radiation and the accretion rate. For a given spacetime geometry, the model predicts that QPO frequencies are function of single parameter $r_0$. We illustrate that the model can provide fits of data comparable to those reached by other models, or even better. In particular it is compared to relativistic precession model. Moreover, we also illustrate that the model consideration is compatible with consideration of models of a rotating neutron star in the atoll source 4U~1636-53.

The Resonant Switch (RS) model of twin high-frequency quasi-periodic oscillations (HF QPOs) observed in neutron star binary systems, based on switch of the twin oscillations at a resonant point, has been applied to the atoll source 4U 1636-53 under assumption that the neutron star exterior can be approximated by the Kerr geometry. Strong restrictions of the neutron star parameters M (mass) and a (spin) arise due to fitting the frequency pairs admitted by the RS model to the observed data in the regions related to the resonant points. The most precise variants of the RS model are those combining the relativistic precession frequency relations with their modifications. Here, the neutron star mass and spin estimates given by the RS model are confronted with a variety of equations of state (EoS) governing structure of neutron stars in the framework of the Hartle-Thorne theory of rotating neutron stars applied for the observationally given rotation frequency f_rot~580 Hz (or alternatively f_rot~290 Hz) of the neutron star at 4U 1636-53. It is shown that only two variants of the RS model based on the Kerr approximation are compatible with two EoS applied in the Hartle-Thorne theory for f_rot~580 Hz, while no variant of the RS model is compatible for f_rot~290 Hz. The two compatible variants of the RS model are those giving the best fits of the observational data. However, a self-consistency test by fitting the observational data to the RS model with oscillation frequencies governed by the Hartle-Thorne geometry described by three spacetime parameters M, a and (quadrupole moment) q related by the two available EoS puts strong restrictions. The test admits only one variant of the RS model of twin HF QPOs for the Hartle-Thorne theory with the Gandolfi et al. (2010) EoS predicting the parameters of the neutron star $M \sim 2.10 \mathrm{M}_{\odot}$, $a \sim 0.208$, and $q/a^2 \sim 1.77$.

Isolated pulsars are rotating neutron stars with accurately measured angular velocities $\Omega$, and their time derivatives which show unambiguously that the pulsars are slowing down. Although the exact mechanism of the spin-down is a question of debate in detail, the commonly accepted view is that it arises through emission of magnetic dipole radiation (MDR) from a rotating magnetized body. Other processes, including the emission of gravitational radiation, and of relativistic particles (pulsar wind), are also being considered. The calculated energy loss by a rotating pulsar with a constant moment of inertia is assumed proportional to a model dependent power of $\Omega$. This relation leads to the power law $\dot{\Omega}$ = -K $\Omega^{\rm n}$ where $n$ is called the braking index. The MDR model predicts $n$ exactly equal to 3. Selected observations of isolated pulsars provide rather precise values of $n$, individually accurate to a few percent or better, in the range 1$ <$ n $ < $ 2.8, which is consistently less than the predictions of the MDR model. In spite of an extensive investigation of various modifications of the MDR model, no satisfactory explanation of observation has been found as yet. We employ four physically realistic equations of state, and two computational codes, to model the dynamical effects of rotation on the braking index in the MDR model. In addition to this we simulate an effect on moment of inertia where we assume a certain amount of superfluid matter has manifest between the crust and core region of the star thus essentially eliminating momentum transfer between the two regions. We find that the effects of rotation on braking index are significant at high frequencies, but have little effect at frequencies consistent with the most accurately measured pulsars to-date.

Twin peak quasi-periodic oscillations (QPOs) appear in the X-ray power-density spectra of several accreting low-mass neutron star (NS) binaries. Observations of the peculiar Z-source Circinus X-1 display unusually low QPO frequencies. Using these observations, we have previously considered the relativistic precession (RP) twin peak QPO model to estimate the mass of central NS in Circinus X-1. We have shown that such an estimate results in a specific mass-angular-momentum (M-j) relation rather than a single preferred combination of M and j. Here we confront our previous results with another binary, the atoll source 4U 1636-53 that displays the twin peak QPOs at very high frequencies, and extend the consideration to various twin peak QPO models. In analogy to the RP model, we find that these imply their own specific M-j relations. We explore these relations for both sources and note differences in the chi-square behavior that represent a dichotomy between high- and low-frequency sources. Based on the RP model, we demonstrate that this dichotomy is related to a strong variability of the model predictive power across the frequency plane. This variability naturally comes from the radial dependence of characteristic frequencies of orbital motion. As a consequence, the restrictions on the models resulting from observations of low-frequency sources are weaker than those in the case of high-frequency sources. Finally we also discuss the need for a correction to the RP model and consider the removing of M-j degeneracies, based on the twin peak QPO-independent angular momentum estimates.

The innermost stable cicular orbit (ISCO) of an accretion disc orbiting a neutron star (NS) is often assumed a unique prediction of general relativity. However, it has been argued that ISCO also appears around highly elliptic bodies described by Newtonian theory. In this sense, the behaviour of an ISCO around a rotating oblate neutron star is formed by the interplay between relativistic and Newtonian effects. Here we briefly explore the consequences of this interplay using a straightforward analytic approach as well as numerical models that involve modern NS equations of state. We examine the ratio K between the ISCO radius and the radius of the neutron star. We find that, with growing NS spin, the ratio K first decreases, but then starts to increase. This non-monotonic behaviour of K can give rise to a neutron star spin interval in which ISCO appears for two very different ranges of NS mass. This may strongly affect the distribution of neutron stars that have an ISCO (ISCO-NS). When (all) neutron stars are distributed around a high mass M0, the ISCO-NS spin distribution is roughly the same as the spin distribution corresponding to all neutron stars. In contrast, if M0 is low, the ISCO-NS distribution can only have a peak around a high value of spin. Finally, an intermediate value of M0 can imply an ISCO-NS distribution divided into two distinct groups of slow and fast rotators. Our findings have immediate astrophysical applications. They can be used for example to distinguish between different models of high-frequency quasiperiodic oscillations observed in low-mass NS X-ray binaries.

We present results for models of neutron stars and strange stars constructed using the Hartle-Thorne slow-rotation method with a wide range of equations of state, focusing on the values obtained for the angular momentum $J$ and the quadrupole moment $Q$, when the gravitational mass $M$ and the rotational frequency $\Omega$ are specified. Building on previous work, which showed surprising uniformity in the behaviour of the moment of inertia for neutron-star models constructed with widely-different equations of state, we find similar uniformity for the quadrupole moment. These two quantities, together with the mass, are fundamental for determining the vacuum space-time outside neutron stars. We study particularly the dimensionless combination of parameters $QM/J^2$ (using units for which $c=G=1$). This quantity goes to 1 in the case of a Kerr-metric black hole and deviations away from 1 then characterize the difference between neutron-star and black-hole space-times. It is found that $QM/J^2$ for both neutron stars and strange stars decreases with increasing mass, for a given equation of state, reaching a value of around 2 (or even less) for maximum-mass models, meaning that their external space-time is then rather well approximated by the Kerr metric. If $QM/J^2$ is plotter against compactness $R/2M$ (where $R$ is the radius), it is found that the relationship is nearly unique for neutron-star models, independent of the equation of state, while it is significantly different for strange stars. This gives a new way of possibly distinguishing between them.

We study non-geodesic corrections to the quasicircular motion of charged test particles in the field of magnetized slowly rotating neutron stars. The gravitational field is approximated by the Lense-Thirring geometry, the magnetic field is of the standard dipole character. Using a fully-relativistic approach we determine influence of the electromagnetic interaction (both attractive and repulsive) on the quasicircular motion. We focus on the behaviour of the orbital and epicyclic frequencies of the motion. Components of the four-velocity of the orbiting charged test particles are obtained by numerical solution of equations of motion, the epicyclic frequencies are obtained by using the standard perturbative method. The role of the combined effect of the neutron star magnetic field and its rotation in the character of the orbital and epicyclic frequencies is discussed.

We consider twin-peak quasi-periodic oscillations observed in the accreting low-mass neutron star binaries and explore restrictions to central compact object properties that are implied by various QPO models. For each model and each source, the consideration results in a specific relation between the compact object mass $M$ and the angular-momentum $j$ rather than in their single preferred combination. Moreover, restrictions on the models resulting from observations of the low-frequency sources are weaker than those in the case of the high-frequency sources.

Extremely Compact Stars (ECS) contain trapped null geodesics. When such objects enter the evolution period admitting geodetical motion of neutrinos, certain part of neutrinos produced in their interior will be trapped influencing their neutrino luminosity and thermal evolution. We study neutrino trapping in the braneworld ECS, assuming uniform distribution of neutrino emissivity and massless neutrinos. We give the efficiency of the neutrino trapping effects in the framework of the simple model of the internal spacetime with uniform distribution of energy density, and external spacetime described by the Reissner-Nordstr\"om geometry characterized by the braneworld "tidal" parameter $b$. For $b < 0$ the external spacetime is of the black-hole type, while for $b > 0$ the external spacetime can be of both black-hole and naked-singularity type. Then the ECS surface radius $R$ can be located also above the unstable (outer) photon circular orbit. Such basically new types of the spacetimes strongly alter the trapping phenomena as compared to the standard case of $b = 0$. It is shown that the neutrino trapping effects are slightly lowered by the presence of physically more plausible case of $b < 0$, as compared to the standard internal Schwarzschild spacetime, while they can be magnified by positive tidal charges if $b < 1$ and lowered for $b > 1$. However, potential astrophysical relevance of the trapping phenomena is strongly enhanced for negative tidal charges enabling a significant enlargement of the ECS surface radius to values coherent with recent observations.

Boutloukos et al. (2006) discovered twin-peak quasi-periodic oscillations (QPOs) in 11 observations of the peculiar Z-source Circinus X-1. Among several other conjunctions the authors briefly discussed the related estimate of the compact object mass following from the geodesic relativistic precession model for kHz QPOs. Neglecting the neutron star rotation they reported the inferred mass M_0 = 2.2 +/- 0.3 M_\sun. We present a more detailed analysis of the estimate which involves the frame-dragging effects associated with rotating spacetimes. For a free mass we find acceptable fits of the model to data for (any) small dimensionless compact object angular momentum j=cJ/GM^2. Moreover, quality of the fit tends to increase very gently with rising j. Good fits are reached when M ~ M_0[1+0.55(j+j^2)]. It is therefore impossible to estimate the mass without the independent knowledge of the angular momentum and vice versa. Considering j up to 0.3 the range of the feasible values of mass extends up to 3M_\sun. We suggest that similar increase of estimated mass due to rotational effects can be relevant for several other sources.

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.

Set of neutron star observational results is used to test some selected equations of state of dense nuclear matter. The first observational result comes from the mass--baryon number relation for pulsar B of the double pulsar system J 0737--3039. The second one is based on the mass--radius relation coming from observation of the thermal radiation of the neutron star RX J 1856.35--3754. The third one follows the population analysis of isolated neutron star thermal radiation sources. The last one is the test of maximum mass. The equation of state of asymmetric nuclear matter is given by the parameterized form of the relativistic Brueckner-Hartree-Fock mean field, and we test selected parameterizations that represent fits of full relativistic mean field calculation. We show that only one of them is capable to pass the observational tests. This equation of state represents the first equation of state that is able to explain all the mentioned observational tests, especially the very accurate test given by the double pulsar even if no mass loss is assumed.

Assuming a resonant origin of the quasiperiodic oscillations observed in the X-ray neutron star binary systems, we apply a genetic algorithm method for selection of neutron star models. It was suggested that pairs of kilo-Hertz peaks in the X-ray Fourier power density spectra of some neutron stars reflect a non-linear resonance between two modes of accretion disk oscillations. In several specific models, the two modes are related to physically plausible combinations of Keplerian, vertical and radial frequencies of geodesic orbital motion. We investigate this concept for a specific neutron star source, a fixed pair of modes and various neutron star equations of state. Each neutron star model is characterized by the equation of state (EOS), rotation frequency ($\Omega$) and central energy density ($\rho_\mathrm c$). These determine the spacetime structure governing geodesic motion and position dependent radial and vertical epicyclic oscillations related to the stable circular geodesics. When the parameters of neutron star model are fixed, the two considered modes imply a frequency-frequency relation which can be compared to the observation in order to eliminate the unsatisfactory sets (KR,$\rho_\mathrm c, \Omega$, EOS). For the elimination we use the advanced genetic algorithm. Genetic algorithm comes out from the method of natural selection when subjects with the best adaptation to assigned conditions have best chances to survive. The chosen genetic algorithm with sexual reproduction contains one chromosome with restricted lifetime, uniform crossing and genes of type 3/3/5. For encryption of physical description (KR,$\rho_\mathrm c, \Omega$, EOS) into chromosome we use the Gray code. As a fitness function we use correspondence between the observed and calculated pairs of eigenfrequencies.

Trapping of neutrinos in extremely compact stars containing trapped null geodesics is studied. We calculated the ratio of produced to trapped neutrinos in the simplest model of uniform density stars. This gives the upper limit on trapping coefficients in real objects.