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}$.

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.