In this paper we show that X-ray spectral observations of the ATHENA mission, which is planned to launch in 2031, can constrain the equation of state of superdense matter. We use our well-constrained continuum fitting method for mass and radius determination of the neutron star. Model spectra of the emission from a neutron star were calculated using the atmosphere code ATM24. In the next step, those models were fitted to a simulated spectra of the neutron star calculated for ATHENA's WFI detector, using the satellite calibration files. To simulate the spectra we assumed three different values of effective temperatures, surface gravities and gravitational redshifts. There cases are related to the three different neutron star masses and radii. This analysis allows us to demonstrate the precision of our method and demonstrate the need for a fast detector onboard of ATHENA. A large grid of theoretical spectra was calculated with various parameters and a hydrogen-helium-iron composition of solar proportion. These spectra were fitted to the simulated spectrum to estimate the precision of mass and radius determination. In each case, we obtained very precise mass and radius values with errors in the range 3--10% for mass and in the range 2--8% for radius within the 1-sigma confidence error. We show here that with the ATHENA WFI detector, such a determination could be used to constrain the equation of state of superdense neutron star matter.

In this paper we fit the NuSTAR and XMM-Newton data of three sources: NGC7793~P13, NGC5907~ULX1, and Circinus~ULX5. Our single model contains emission form non-spherical system: neutron star plus accretion disk directed towards observer. We obtained the very good fit with the reduced $\chi^2$ per degree of freedom equal 1.08 for P13, 1.01 for ULX1, and 1.14 for ULX5.The normalization of our model constrains the distance to the source. The resulting distances are $D=3.41^{+0.11}_{-0.10}$, $6.55_{-0.81}^{+0.69}$ and $2.60^{+0.05}_{-0.03}$~Mpc for P13, ULX1 and ULX5 respectively. The distances to P13 and ULX5 are in perfect agreement with previous distance measurements to their host galaxies. Our results confirm that P13, ULX1 and ULX5 may contain central hot neutron star. When the outgoing emission is computed by integration over the emitting surface and successfully fitted to the data, then the resulting model normalization is the direct distance indicator.

Thermal dominated X-ray spectra of neutron stars in quiescent transient X-ray binaries and neutron stars that undergo thermonuclear bursts are sensitive to mass and radius. The mass-radius relation of neutron stars depends on the equation of state that governs their interior. Constraining this relation accurately is thus of fundamental importance to understand the nature of dense matter. In this context we introduce a pipeline to calculate realistic model spectra of rotating neutron stars with hydrogen and helium atmospheres. An arbitrarily fast rotating neutron star with a given equation of state generates the spacetime in which the atmosphere emits radiation. We use the Lorene/nrotstar code to compute the spacetime numerically and the ATM24 code to solve the radiative transfer equations self-consistently. Emerging specific intensity spectra are then ray-traced through the neutron star's spacetime from the atmosphere to a distant observer with the Gyoto code. Here, we present and test our fully relativistic numerical pipeline. To discuss and illustrate the importance of realistic atmosphere models we compare our model spectra to simpler models like the commonly used isotropic color-corrected blackbody emission. We highlight the importance of considering realistic model-atmosphere spectra together with relativistic ray tracing to obtain accurate predictions. We also insist on the crucial impact of the star's rotation on the observables. Finally, we close a controversy that has been appearing in the literature in the recent years regarding the validity of the ATM24 code.

We study the appearance of the neutron star - accretion disk system as seen by a distant observer in the UV/X-ray domain. The observed intensity spectra are computed assuming non-spherical geometry of the whole system, in which outgoing spectrum is not represented by the flux spectrum, the latter being valid for spherically symmetric objects. Intensity spectra of our model display double bumps in UV/X-ray energy domains. Such structure is caused by the fact that the the source is not spherically symmetric, and the proper integration of intensity over emitted area is needed to reproduce observed spectral shape. Relative normalization of double bump is self consistently computed by our model. X-ray spectra of such a type were often observed in LMXB with accretion disk, ultra luminous X-ray sources, and accreting black hole systems with hot inner compact corona. Our model naturally explains high energy broadening of the disk spectrum observed in some binaries. We attempted to fit our model to X-ray data of XTE~J1709-267 from {\it XMM-Newton}. Unfortunately, the double intensity bump predicted by our model for LMXB is located in soft X-ray domain, uncovered by existing data for this source.

We have studied the future evolution of a recently discovered ULX source P13 in NGC7793. This source was shown to contain a 5-15$M_{\odot}$ black hole and a massive 18-23$M_{\odot}$ B9Ia companion on a 64 day orbit. For low black hole mass (5-10$M_{\odot}$) and high companion mass $\gtrsim 20 M_{\odot}$ the binary is predicted to initiate a common envelope evolution in near future and significantly decrease the orbital separation (6 hr orbit). This leads to a high probability ($\sim70\%$) of the system surviving a supernova explosion that will form a neutron star out of the companion. About one third of the surviving BH-NS systems will merge within Hubble time and be a source of high frequency gravitational radiation. We estimate that the chances of detection of BH-NS systems with advanced LIGO/Virgo that form via P13-evolutionary channel are at the level of 0.1 yr$^{-1}$ with wide range of allowed probability $0 \div 0.6$yr$^{-1}$. This is the fourth empirical estimate of BH-NS merger rate.