We investigate the three-dimensional motion of a test particle in the gravitational field generated by a non-spherical compact object endowed with a mass quadrupole moment, described by the Erez-Rosen metric, and a radiation field, including the general relativistic Poynting-Robertson effect, coming from a rigidly rotating spherical emitting source located outside of the compact object. We derive the equations of motion for test particles influenced by such radiation field, recovering the two-dimensional description, and the weak-field approximation. This dynamical system admits the existence of a critical hypersurface, region where gravitational and radiation forces balance. Selected test particle orbits for different set of input parameters are displayed. The possible configurations on the critical hypersurfaces can be either latitudinal drift towards the equatorial ring or suspended orbits. We discuss about the existence of multiple hypersurface solutions through a simple method to perform the calculations. We graphically prove also that the critical hypersurfaces are stable configurations within the Lyapunov theory.

We present a new method to measure the radius-to-mass ratio (R/M) of weakly magnetic, disc-accreting neutron stars by exploiting the occultation of parts of the inner disc by the star itself. This occultation imprints characteristic features on the X-ray line profile that are unique and are expected to be present in low mass X-ray binary systems seen under inclinations higher than ~65 degrees. We analyse a NuSTAR observation of a good candidate system, 4U 1636-53, and find that X-ray spectra from current instrumentation are unlikely to single out the occultation features owing to insufficient signal-to-noise. Based on an extensive set of simulations we show that large-area X-ray detectors of the future generation could measure R/M to ~2{\div}3% precision over a range of inclinations. Such is the precision in radius determination required to derive tight constraints on the equation of state of ultradense matter and it represents the goal that other methods too aim to achieve in the future.

We analyzed the thermodynamics of hypothetical exoplanets at very low Keplerian circular orbits in close vicinity of rapidly spinning supermassive black holes. Such black hole exoplanets are heated by strongly blueshifted and focused flux of the incoming cosmic microwave background (CMB) and cooled by the cold part of the local sky containing the black hole shadow. This gives rise to a temperature difference, which can drive processes far from thermodynamic equilibrium in a hypothetical life form inhabiting black hole exoplanets, similar to the case of a planet heated by the radiation of the parent star and cooled by the night sky. We found that for a narrow range of radii of very low Keplerian circular orbits and for very high spin of a supermassive black hole, the temperature regime of the black hole exoplanets corresponds to the habitable zone around standard stars. The thermodynamics of black hole exoplanets therefore, in principle, does not exclude the existence of life based on known biology. The peak of the multiblackbody spectral profile of the CMB heating the exoplanet is located in the ultraviolet band, but a significant fraction of the flux comes also in the visible and infrared bands. The minimum mass of a black hole ensuring the resistance to tidal disruption of an Earth-like exoplanet orbiting in the habitable zone is estimated to $1.63 \cdot 10^8 \, m_{\odot}$.

We prove the stability of the critical hypersurfaces associated with the three-dimensional general relativistic Poynting-Robertson effect. The equatorial ring configures to be as a stable attractor and the whole critical hypersurface as a basin of attraction for this dynamical system. We introduce a new, simpler (in terms of calculations), and more physical approach within the Lyapunov theory. We propose three different Lyapunov functions, each one carrying important information and very useful for understanding such phenomenon under different aspects.

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.

Life is dependent on the income of energy with low entropy and the disposal of energy with high entropy. On Earth, the low-entropy energy is provided by solar radiation and the high-entropy energy is disposed as infrared radiation emitted into the cold space. Here we turn the situation around and assume cosmic background radiation as the low-entropy source of energy for a planet orbiting a black hole into which the high-entropy energy is disposed. We estimate the power that can be produced by thermodynamic processes on such a planet, with a particular interest in planets orbiting a fast rotating Kerr black hole as in the science fiction movie {\em Interstellar}. We also briefly discuss a reverse Dyson sphere absorbing cosmic background radiation from the outside and dumping waste energy to a black hole inside.

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.

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.

High frequency (millisecond) quasi-periodic oscillations (HF QPOs) are observed in the X-ray power-density spectra of several microquasars and low mass X-ray binaries. Two distinct QPO peaks, so-called twin peak QPOs, are often detected simultaneously exhibiting their frequency ratio close or equal to 3/2. Following the analytic theory and previous studies of observable spectral signatures, we aim to model the twin peak QPOs as a spectral imprint of specific dual oscillation regime defined by a combination of the lowest radial and vertical oscillation mode of optically thick slender tori with constant specific angular momentum. We examined power spectra and fluorescent K$\alpha$ iron line profiles for two different simulation setups with the mode frequency relations corresponding to the epicyclic resonance HF QPOs model and modified relativistic precession QPOs model. We use relativistic ray-tracing implemented in parallel simulation code LSDplus. In the background of the Kerr spacetime geometry, we analyze the influence of the distant observer inclination and the spin of the central compact object. Relativistic optical projection of the oscillating slender torus is illustrated by images in false colours related to the frequency shift. We show that performed simulations yield power spectra with the pair of dominant peaks corresponding to the frequencies of radial and vertical oscillation modes with the proper ratio equal to 3/2 on a wide range of inclinations and spin values. We also discuss exceptional cases of a very small and very high inclination as well as unstable high spin relativistic precession-like configuration predicting constant frequency ratio equal to 1/2. We demonstrate signifiant dependency of broadened K$\alpha$ iron line profiles on the inclination of the distant observer.

Estimations of black hole spin in the three Galactic microquasars GRS 1915+105, GRO J1655-40, and XTE J1550-564 have been carried out based on spectral and timing X-ray measurements and various theoretical concepts. Among others, a non-linear resonance between axisymmetric epicyclic oscillation modes of an accretion disc around a Kerr black hole has been considered as a model for the observed high-frequency quasi-periodic oscillations (HF QPOs). Estimates of spin predicted by this model have been derived based on the geodesic approximation of the accreted fluid motion. Here we assume accretion flow described by the model of a pressure-supported torus and carry out related corrections to the mass-spin estimates. We find that for dimensionless black hole spin a<0.9, the resonant eigenfrequencies are very close to those calculated for the geodesic motion. Their values slightly grow with increasing torus thickness. These findings agree well with results of a previous study carried out in the pseudo-Newtonian approximation. The situation becomes different for a>0.9, in which case the resonant eigenfrequencies rapidly decrease as the torus thickness increases. We conclude that the assumed non-geodesic effects shift the lower limit of the spin, implied for the three microquasars by the epicyclic model and independently measured masses, from a~0.7 to a~0.6. Their consideration furthermore confirms compatibility of the model with the rapid spin of GRS 1915+105 and provides highly testable predictions of the QPO frequencies. Individual sources with a moderate spin (a<0.9) should exhibit a smaller spread of the measured 3:2 QPO frequencies than sources with a near-extreme spin (a~1). This should be further examined using the large amount of high-resolution data expected to become available with the next generation of X-ray instruments, such as the proposed Large Observatory for X-ray Timing (LOFT).

High frequency quasi-periodic oscillations (HF QPOs) are observed in the X-ray power-density spectra (PDS) of several microquasars and low mass X-ray binaries. Many proposed QPO models are based on oscillations of accretion toroidal fluid structures orbiting in the vicinity of a compact object. We study oscillating accretion tori orbiting in the vicinity of a Kerr black hole. We demonstrate that significant variation of the observed flux can be caused by the combination of radial and vertical oscillation modes of a slender, polytropic, perfect fluid, non-self-graviting torus with constant specific angular momentum. We investigate two combinations of the oscillating modes corresponding to the direct resonance QPO model and the modified relativistic precession QPO model.

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

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.

In the context of high-frequency quasi-periodic oscillation (HF QPOs) we further explore the appearance of an observable signal generated by hot spots moving along quasi-elliptic trajectories close to the innermost stable circular orbit in the Schwarzschild spacetime. The aim of our investigation is to reveal whether observable characteristics of the Fourier power-spectral density can help us to distinguish between the two competing models, namely, the idea of bright spots orbiting on the surface of an accretion torus versus the scenario of intrinsic oscillations of the torus itself. We take the capabilities of the present observatories (represented by the Rossi X-ray Timing Explorer, RXTE) into account, and we also consider the proposed future instruments (represented here by the Large Observatory for X-ray Timing, LOFT).

Context. Some microquasars exhibit millisecond quasi-periodic oscillations (QPO) that are likely related to phenomena occuring in the immediate vicinity of the central black hole. Oscillations of accretion tori have been proposed to model these QPOs. Aims. Here, we aim at determining the observable spectral signature of slender accretion tori surrounding Kerr black holes. We analyze the impact of the inclination and spin parameters on the power spectra. Methods. Ray-traced power spectra of slender tori oscillation modes are computed in the Kerr metric. Results. We show that the power spectral densities of oscillating tori are very sensitive to the inclination and spin parameters. This strong dependency of the temporal spectra on inclination and spin may lead to observable constraints of these parameters. Conclusions. This work goes a step further in the analysis of the oscillating torus QPO model. It is part of a long-term study that will ultimately lead to comparison with observed data.

High frequency quasi-periodic oscillations (HF QPOs) appear in the X-ray variability of several accreting low-mass binaries. In a series of works it was suggested that these QPOs may have connection to inhomogeneities orbiting close to an inner edge of the accretion disc. In this paper we explore the appearance of an observable signal generated by small radiating circular hot spots moving along quasi-elliptic trajectories close to the innermost stable circular orbit in the Schwarzschild spacetime. Our consideration takes into account the capabilities of observatories that have been operating in the past two decades represented by the Rossi X-ray Timing Explorer (RXTE) and the proposed future instruments represented by the Large Observatory for X-ray Timing (LOFT). For these purposes we choose such model parameters that lead to lightcurves comparable to those observed in Galactic black hole sources, in particular the microquasar GRS 1915+105. We find that when a weak signal corresponding to the hot-spot Keplerian frequency is around the limits of the RXTE detectability, the LOFT observations can clearly reveal its first and second harmonics. Moreover, in some specific situations the radial epicyclic frequency of the spot can be detected as well. Finally, we also compare the signal produced by the spots to the signal produced by axisymmetric epicyclic disc-oscillation modes and discuss the key differences that could be identified via the proposed future technology. We conclude that the ability to recognize the harmonic content of the signal can help to distinguish between the different proposed physical models.

Context. One of the often discussed models for X-ray binaries high-frequency quasi-periodic oscillations is the oscillating torus model that considers oscillation modes of slender accretion tori. Aims. Here, we aim at developing this model by considering the observable signature of an optically thick slender accretion torus subject to simple periodic deformations. Methods. We compute light curves and power spectra of a slender accretion torus subject to simple periodic deformations: vertical or radial translation, rotation, expansion and shear. Results. We show that different types of deformations lead to very different Fourier power spectra and therefore could be observationally distinguished. Conclusions. This work is a first step in a longer-term study of the observable characteristics of the oscillating torus model. It gives promising perspectives on the possibility to constrain this model by studying the observed power spectra of quasi-periodic oscillations.

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.

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.

A long discussion has been devoted to the issue of clustering of the kHz QPO frequency ratios in the neutron star sources. While the distribution of ratios inferred from an occurrence of a single QPO seems to be consistent with a random walk, the distribution based on simultaneous detections of both peaks indicates a preference of ratios of small integers. Based on the public RXTE data we further investigate this issue for the source 4U 1636-53. Quality factors and rms amplitudes of both the QPOs nearly equal to the points where the frequencies are commensurable, and where the twin QPO detections cluster. We discuss a connection of the clustering with the varying properties of the two QPO modes. Assuming approximative relations for the observed correlations of the QPO properties, we attempt to reproduce the frequency and ratio distributions using a simple model of a random-walk evolution along the observed frequency-frequency correlation. We obtain results which are in qualitative agreement with the observed distributions.

A recently published study on long term evolution of the frequencies of the kilohertz quasi-periodic oscillations (QPOs) in the atoll source 4U 1636-53 concluded that there is no preferred frequency ratio in a distribution of twin QPOs that was inferred from the distribution of a single frequency alone. However, we find that the distribution of the ratio of actually observed pairs of kHz QPO frequencies is peaked close to the 3/2 value, and possibly also close to the 5/4 ratio. To resolve the apparent contradiction between the two studies, we examine in detail the frequency distributions of the lower kHz QPO and the upper kHz QPO detected in our data set. We demonstrate that for each of the two kHz QPOs (the lower or the upper), the frequency distribution in all detections of a QPO differs from the distribution of frequency of the same QPO in the subset of observations where both the kHz QPOs are detected. We conclude that detections of individual QPOs alone should not be used for calculation of the distribution of the frequency ratios.

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.

We have developed a realistic, fully general relativistic computer code to simulate optical projection in a strong, spherically symmetric gravitational field. The standard theoretical analysis of optical projection for an observer in the vicinity of a Schwarzschild black hole is extended to black hole spacetimes with a repulsive cosmological constant, i.e, Schwarzschild-de Sitter spacetimes. Influence of the cosmological constant is investigated for static observers and observers radially free-falling from the static radius. Simulations include effects of the gravitational lensing, multiple images, Doppler and gravitational frequency shift, as well as the intensity amplification. The code generates images of the sky for the static observer and a movie simulations of the changing sky for the radially free-falling observer. Techniques of parallel programming are applied to get a high performance and a fast run of the BHC simulation code.

The results we presented were biased because of the typing error in the code transcription of the equation (A.12). We will replace by the corrected version soon. The brief of the corrected results can be found at http://www.physics.cz/research/doc/posters/1181038112_0.pdf