It has been recently discovered that the transition of X-ray pulsars to the low luminosity state ($L\lesssim 10^{35}\,{\rm erg\,s^{-1}}$) is accompanied by a dramatic spectral changes. Namely, the typical power-law-like spectrum with high energy cutoff transforms into a two-component structure with a possible cyclotron absorption feature on top of it. It was proposed that these spectral characteristics can be explained qualitatively by the emission of cyclotron photons in the atmosphere of the neutron star caused by collisional excitation of electrons to upper Landau levels and further comptonization of the photons by electron gas. The latter is expected to be overheated in a thin top layer of the atmosphere. In this paper, we perform Monte Carlo simulations of the radiative transfer in the atmosphere of an accreting neutron star while accounting for a resonant scattering of polarized X-ray photons by thermally distributed electrons. The spectral shape is shown to be strongly polarization-dependent in soft X-rays ($\lesssim 10\,{\rm keV}$) and near the cyclotron scattering feature. The results of our numerical simulations are tested against the observational data of X-ray pulsar A 0535+262 in the low luminosity state. We show that the spectral shape of the pulsar can be reproduced by the proposed theoretical model. The applications of the discovery to the observational studies of accreting neutron stars are discussed.

Neutron stars (NSs) in low-mass X-ray binaries rotate at frequencies high enough to significantly deviate from sphericity ($\nu_* \sim$ 200--600 Hz). We investigate the effects of rapid rotation on the observational appearance of a NS. We propose analytical formulae relating gravitational mass and equatorial radius of the rapidly rotating NS to the mass $M$ and radius $R$ of a non-rotating NS of the same baryonic mass using accurate fully relativistic computations. We compute spectra from an oblate rotating NS observed at different inclination angles using the modified oblate Schwarzschild (MOS) approximation, where light bending is computed in Schwarzschild metric, but frame dragging and quadrupole moment of a NS are approximately accounted for in the photon redshift calculations. We generalize the cooling tail method to the case of a rapidly rotating NS to obtain the most probable values of $M$ and $R$ of the corresponding non-rotating NS with the same baryonic mass. We approximate the local spectra from the NS surface by a diluted blackbody using previously computed NS atmosphere models. We show that the NS radius could be overestimated by 3--3.5 km for face-on stars of $R\approx 11$ km rotating at $\nu_* =$ 700 Hz if the version of the cooling tail method for a non-rotating NS is used. We apply the method to an X-ray burst observed from the NS rotating at $\nu_* \approx$ 532 Hz in SAX J1810.8$-$2609. The resulting radius of the non-rotating NS (assuming $M=1.5 M_\odot$) becomes $11.8\pm0.5$ km if it is viewed at inclination i=60 deg and $R=11.2\pm0.5$ km for a face-on view, which are smaller by 0.6 and 1.2 km than the radius obtained using standard cooling tail method ignoring rotation. The corresponding equatorial radii of these rapidly rotating NSs are 12.3$\pm 0.6$ km (for i=60 deg) and 11.6$\pm 0.6$\,km (for i=0 deg).

Black hole X-ray binaries show signs of non-thermal emission in the optical/near-infrared range. We analyze the optical/near-infrared SMARTS data on GX339$-$4 over the 2002--2011 period. Using the soft state data, we estimate the interstellar extinction towards the source and characteristic color temperatures of the accretion disk. We show that various spectral states of regular outbursts occupy similar regions on the color-magnitude diagrams, and that transitions between the states proceed along the same tracks despite substantial differences in the observed light curves morphology. We determine the typical duration of the hard-to-soft and soft-to-hard state transitions and the hard state at the decaying stage of the outburst to be one, two and four weeks, respectively. We find that the failed outbursts cannot be easily distinguished from the regular ones at their early stages, but if the source reaches 16 mag in $V$-band, it will transit to the soft state. By subtracting the contribution of the accretion disk, we obtain the spectra of the non-thermal component, which have constant, nearly flat shape during the transitions between the hard and soft states. In contrast to the slowly evolving non-thermal component seen at optical and near-infrared wavelengths, the mid-infrared spectrum is strongly variable on short timescales and sometimes shows a prominent excess with a cutoff below $10^{14}$ Hz. We show that the radio to optical spectrum can be modeled using three components corresponding to the jet, hot flow and irradiated accretion disk.

We compute accurate atmosphere models of rotation-powered millisecond pulsars in which the polar caps of a neutron star (NS) are externally heated by magnetospheric return currents. The external ram pressure, energy losses, and stopping depth of the penetrating charged particles are computed self-consistently with the atmosphere model, instead of assuming a simplified deep-heated atmosphere in radiative equilibrium. We use exact Compton scattering formalism to model the properties of the emergent X-ray radiation. The deep-heating approximation is found to be valid only if most of the heat originates from ultra-relativistic bombarding particles with Lorentz factors $\gamma \gtrsim 100$. In the opposite regime, the atmosphere attains a distinct two-layer structure with an overheated optically thin skin on top of an optically thick cool plasma. The overheated skin strongly modifies the emergent radiation producing a Compton-upscattered high-energy tail in the spectrum and changes the radiation beaming pattern from limb darkening to limb brightening for emitted hard X-rays. This kind of drastic change in the emission properties can have a significant impact on the inferred NS pulse profile parameters as performed, for example, by Neutron star Interior Composition ExploreR (NICER). Finally, the connection between the energy distribution of the return current particles and the atmosphere emission properties offers a new tool for probing the exact physics of pulsar magnetospheres.

We report on our analysis of the 2019 outburst of the X-ray accreting pulsar 4U 1901+03 observed with Insight-HXMT and NICER. Both spectra and pulse profiles evolve significantly in the decaying phase of the outburst. Dozens of flares are observed throughout the outburst. They are more frequent and brighter at the outburst peak. We find that the flares, which have a duration from tens to hundreds of seconds, are generally brighter than the persistent emission by a factor of $\sim$ 1.5. The pulse profile shape during the flares can be significantly different than that of the persistent emission. In particular, a phase shift is clearly observed in many cases. We interpret these findings as direct evidence of changes of the pulsed beam pattern, due to transitions between the sub- and super-critical accretion regimes on a short time scale. We also observe that at comparable luminosities the flares' pulse profiles are rather similar to those of the persistent emission. This indicates that the accretion on the polar cap of the neutron star is mainly determined by the luminosity, i.e., the mass accretion rate.

Type I X-ray bursts in GS 1826-24, and in several other systems, may induce cooling of the hot inner accretion flow that surrounds the bursting neutron star. Given that GS 1826-24 remained persistently in the hard state over the period 2003-2008 and presented regular bursting properties, we stacked the spectra of the X-ray bursts detected by INTEGRAL (JEM-X and ISGRI) and XMM-Newton (RGS) during that period to study the effect of the burst photons on the properties of the Comptonizing medium. The extended energy range provided by these instruments allows the simultaneous observation of the burst and persistent emission spectra. We detect an overall change in the shape of the persistent emission spectrum in response to the burst photon shower. For the first time, we observe simultaneously a drop in the hard X-ray emission, together with a soft X-ray excess with respect to the burst blackbody emission. The hard X-ray drop can be explained by burst-induced coronal cooling, while the bulk of the soft X-ray excess can be described by fitting the burst emission with an atmosphere model, instead of a simple blackbody model. Traditionally, the persistent emission was assumed to be invariant during X-ray bursts, and more recently to change only in normalization but not in spectral shape; the observed change in the persistent emission level during X-ray bursts may thus trigger the revision of existing neutron star mass-radius constraints, as the derived values rely on the assumption that the persistent emission does not change during X-ray bursts. The traditional burst fitting technique leads to up to a 10% overestimation of the bolometric burst flux in GS 1826-24, which significantly hampers the comparisons of the KEPLER and MESA model against this 'textbook burster'.

Magnetars are believed to host the strongest magnetic fields in the present universe ($B\gtrsim10^{14}$ G) and the study of their persistent emission in the X-ray band offers an unprecendented opportunity to gain insight into physical processes in the presence of ultra-strong magnetic fields. Up to now, most of our knowledge about magnetar sources came from spectral analysis, which allowed to test the resonant Compton scattering scenario and to probe the structure of the star magnetosphere. On the other hand, radiation emitted from magnetar surface is expected to be strongly polarized and its observed polarization pattern bears the imprint of both scatterings onto magnetospheric charges and QED effects as it propagates in the magnetized vacuum around the star. X-ray polarimeters scheduled to fly in the next years will finally allow to exploit the wealth of information stored in the polarization observables. Here we revisit the problem of assessing the spectro-polarimetric properties of magnetar persistent emission. At variance with previous investigations, proper account for more physical surface emission models is made by considering either a condensed surface or a magnetized atmosphere. Results are used to simulate polarimetric observations with the forthcoming Imaging X-ray Polarimetry Explorer (IXPE). We find that X-ray polarimetry will allow to detect QED vacuum effects for all the emission models we considered and to discriminate among them.

The X-ray Polarization Probe (XPP) is a second generation X-ray polarimeter following up on the Imaging X-ray Polarimetry Explorer (IXPE). The XPP will offer true broadband polarimetery over the wide 0.2-60 keV bandpass in addition to imaging polarimetry from 2-8 keV. The extended energy bandpass and improvements in sensitivity will enable the simultaneous measurement of the polarization of several emission components. These measurements will give qualitatively new information about how compact objects work, and will probe fundamental physics, i.e. strong-field quantum electrodynamics and strong gravity.

The aim of this work is to study the possible effects and biases on the radius constraints for rotation-powered millisecond pulsars when using Thomson approximation to describe electron scattering in the atmosphere models, instead of using exact formulation for Compton scattering. We compare the differences between the two models in the energy spectrum and angular distribution of the emitted radiation. We also analyse a self-generated synthetic phase-resolved energy spectrum, based on Compton atmosphere and the most X-ray luminous rotation-powered millisecond pulsars observed by the Neutron star Interior Composition ExploreR (NICER). We derive constraints for the neutron star parameters using both the Compton and Thomson models. The results show that the method works by reproducing the correct parameters with the Compton model. However, biases are found in size and the temperature of the emitting hot spot, when using the Thomson model. The constraints on the radius are still not significantly changed, and therefore the Thomson model seems to be adequate if we are interested only in the radius measurements using NICER.

We present the results of an XMM-Newton observation of the slowly rotating ($P = 3.4$ s), highly magnetized ($B \approx 3 \times 10^{13}$ G) radio pulsar PSR J0726-2612. A previous X-ray observation with the Chandra satellite showed that some of the properties of PSR J0726-2612 are similar to those of the X-ray Dim Isolated Neutron Stars (XDINSs), a small class of nearby slow pulsars characterized by purely thermal X-ray spectra and undetected in the radio band. We confirm the thermal nature of the X-ray emission of PSR J0726-2612, which can be fit by the sum of two blackbodies with temperatures $kT_1 = 0.074_{-0.011}^{+0.006}$ keV and $kT_2 = 0.14_{-0.02}^{+0.04}$ keV and emitting radii $R_1 = 10.4_{-2.8}^{+10.8}$ km and $R_2 = 0.5_{-0.3}^{+0.9}$ km, respectively (assuming a distance of 1 kpc). A broad absorption line modeled with a Gaussian profile centred at $0.39_{-0.03}^{+0.02}$ keV is required in the fit. The pulse profile of PSR J0726-2612 is characterized by two peaks with similar intensity separated by two unequal minima, a shape and pulsed fraction that cannot be reproduced without invoking magnetic beaming of the X-ray emission. The presence of a single radio pulse suggests that in PSR J0726-2612 the angles that the dipole axis and the line of sight make with the rotation axis, $\xi$ and $\chi$ respectively, are similar. This geometry differs from that of the two radio-silent XDINSs with a double peaked pulse profile similar to that of PSR J0726-2612, for which $\xi \sim 90^\circ$ and $\chi \sim 45^\circ$ have been recently estimated. These results strengthen the similarity between PSR J0726-2612 and the XDINSs and support the possibility that the lack of radio emission from the latter might simply be due to an unfavourable viewing geometry.

Deep NuSTAR observation of X-ray pulsar A 0535+262, performed at a very low luminosity of $\sim7\times10^{34}$ erg s$^{-1}$, revealed the presence of two spectral components. We argue that the high-energy component is associated with cyclotron emission from recombination of electrons collisionally excited to the upper Landau levels. The cyclotron line energy of $E_{\rm cyc}=47.7\pm0.8$ keV was measured at the luminosity of almost an order of magnitude lower than what was achieved before. The data firmly exclude a positive correlation of the cyclotron energy with the mass accretion rate in this source.

The X-ray pulsar GRO J1744-28 is a unique source which shows both pulsations and type-II X-ray bursts, allowing studies of the interaction of the accretion disk with the magnetosphere at huge mass accretion rates exceeding $10^{19}$ g s$^{-1}$ during its super-Eddington outbursts. The magnetic field strength in the source, $B\approx 5\times 10^{11}$ G, is known from the cyclotron absorption feature discovered in the energy spectrum around 4.5 keV. Here, we explore the flux variability of the source in context of interaction of its magnetosphere with the radiation-pressure dominated accretion disk. Particularly, we present the results of the analysis of noise power density spectra (PDS) using the observations of the source in 1996-1997 by RXTE. Accreting compact objects commonly exhibit a broken power-law shape of the PDS with a break corresponding to the Keplerian orbital frequency of matter at the innermost disk radius. The observed frequency of the break can thus be used to estimate the size of the magnetosphere. We found, however, that the observed PDS of GRO J1744-28 differs dramatically from the canonical shape. Furthermore, the observed break frequency appears to be significantly higher than what is expected based on the magnetic field estimated from the cyclotron line energy. We argue that these observational facts can be attributed to the existence of the radiation-pressure dominated region in the accretion disk at luminosities above $\sim$2$\times 10^{37}$ erg s$^{-1}$. We discuss a qualitative model for the PDS formation in such disks, and show that its predictions are consistent with our observational findings. The presence of the radiation-pressure dominated region can also explain the observed weak luminosity-dependence of the inner radius, and we argue that the small inner radius can be explained by a quadrupole component dominating the magnetic field of the neutron star.

We analyse new optical spectroscopic, direct-image and X-ray observations of the recently discovered a high proper motion cataclysmic variable V1838 Aql. The data were obtained during its 2013 superoutburst and its subsequent quiescent state. An extended emission around the source was observed up to 30 days after the peak of the superoutburst, interpreted it as a bow--shock formed by a quasi-continuous outflow from the source in quiescence. The head of the bow--shock is coincident with the high--proper motion vector of the source ($v_{\perp}=123\pm5$ km s$^{-1}$) at a distance of $d=202\pm7$ pc. The object was detected as a weak X-ray source ($0.015\pm0.002$ counts s$^{-1}$) in the plateau of the superoutburst, and its flux lowered by two times in quiescence (0.007$\pm$0.002 counts s$^{-1}$). Spectroscopic observations in quiescence we confirmed the orbital period value $P_{\rm{orb}}=0.0545\pm 0.0026$ days, consistent with early-superhump estimates, and the following orbital parameters: $\gamma= -21\pm3$ km s$^{-1}$ and $K_1 = 53\pm3$ km s$^{-1}$. The white dwarf is revealed as the system approaches quiescence, which enables us to infer the effective temperature of the primary $T_{eff}=11,600\pm400$K. The donor temperature is estimated $\lesssim 2200$K and suggestive of a system approaching the period minimum. Doppler maps in quiescence show the presence of the hot spot in HeI line at the expected accretion disc-stream shock position and an unusual structure of the accretion disc in H$\alpha$.

PSR B0943+10 is a mode-switching radio pulsar characterized by two emission modes with different radio and X-ray properties. Previous studies, based on simple combinations of blackbody and power law models, showed that its X-ray flux can be decomposed in a pulsed thermal plus an unpulsed non-thermal components. However, if PSR B0943+10 is a nearly aligned rotator seen pole-on, as suggested by the radio data, it is difficult to reproduce the high observed pulsed fraction unless magnetic beaming is included. In this work we reanalyze all the available X-ray observations of PSR B0943+10 with simultaneous radio coverage, modeling its thermal emission with polar caps covered by a magnetized hydrogen atmosphere or with a condensed iron surface. The condensed surface model provides good fits to the spectra of both pulsar modes, but, similarly to the blackbody, it can not reproduce the observed pulse profiles, unless an additional power law with an ad hoc modulation is added. Instead, the pulse profiles and phase-resolved spectra are well described using the hydrogen atmosphere model to describe the polar cap emission, plus an unpulsed power law. For the X-ray brighter state (Q-mode) we obtain a best fit with a temperature kT~0.09 keV, an emitting radius R~260 m, a magnetic field consistent with the value of the dipole field of 4x10^12 G inferred from the timing parameters, and a small angle between the magnetic and spin axis, $\xi$=5. The corresponding parameters for the X-ray fainter state (B-mode) are kT~0.08 keV and R~170 m.

We report on the discovery of a dramatic change in the energy spectrum of the X-ray pulsar GX 304-1 appearing at low luminosity. Particularly, we found that the cutoff power-law spectrum typical for accreting pulsars, including GX 304-1 at higher luminosities of $L_{\rm X}\sim 10^{36} - 10^{37}$ erg s$^{-1}$, transformed at lower luminosity of $L_{\rm X}\sim 10^{34}$ erg s$^{-1}$ to a two-component spectrum peaking around 5 and 40 keV. We suggest that the observed transition corresponds to a change of the dominant mechanism responsible for the deceleration of the accretion flow. We argue that the accretion flow energy at low accretion rates is released in the atmosphere of the neutron star, and the low-energy component in the source spectrum corresponds to the thermal emission of the optically thick, heated atmospheric layers. The most plausible explanations for the high-energy component are either the cyclotron emission reprocessed by the magnetic Compton scattering or the thermal radiation of deep atmospheric layers partly Comptonized in the overheated upper layers. Alternative scenarios are also discussed.

Some thermonuclear (type I) X-ray bursts at the neutron star surfaces in low-mass X-ray binaries take place during hard persistent states of the systems. Spectral evolution of these bursts is well described by the atmosphere model of a passively cooling neutron star when the burst luminosity is high enough. The observed spectral evolution deviates from the model predictions when the burst luminosity drops below a critical value of 20-70% of the maximum luminosity. We suggest that these deviations are induced by the additional heating of the accreted particles. We present a method for computation of the neutron star atmosphere models heated by accreted particles assuming that their energy is released via Coulomb interactions with electrons. We compute the temperature structures and the emergent spectra of the atmospheres of various chemical compositions and investigate the dependence of the results on the other model parameters. We show that the heated atmosphere develops the hot (20--100 keV) corona-like surface layer cooled by Compton scattering, and the deeper, almost isothermal optically thick region with a temperature of a few keV. The emergent spectra deviate strongly from those of undisturbed neutron star atmospheres, with the main differences being the presence of a high-energy tail and a strong excess in the low-energy part of the spectrum. They also lack the iron absorption edge, which is visible in the spectra of undisturbed low-luminosity atmospheres with solar chemical composition. Using the computed spectra, we obtained the dependences of the dilution and color-correction factors as functions of relative luminosities for pure helium and solar abundance atmospheres. We show that the helium model atmosphere heated by accretion corresponding to 5% of the Eddington luminosity describes well the late stages of the X-ray bursts in 4U 1820-30.

In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to study one common aspect of these objects: their often transient nature. Developed 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.

Thermonuclear X-ray bursts on the surface of neutron stars (NSs) can enrich photosphere with metals, which may imprint the photoionization edges on the burst spectra. We report here the discovery of absorption edges in the spectra of the type I X-ray burst from the NS low-mass X-ray binary GRS 1747$-$312 in Terzan 6 during observations by the Rossi X-ray Timing Explorer. We find that the edge energy evolves from $9.45\pm0.51$ keV to $\sim6$ keV and then back to $9.44\pm0.40$ keV during the photospheric radius expansion phase and remains at $8.06\pm0.66$ keV in the cooling tail. The photoionization absorption edges of hydrogen-like Ni, Fe or Fe/Ni mixture and the bound-bound transitions of metals may be responsible for the observed spectral features. The ratio of the measured absorption edge energy in the cooling tail to the laboratory value of hydrogen-like Ni(Fe) edge energy allows us to estimate the gravitational redshift factor $1+z=1.34\pm0.11$($1+z=1.15\pm0.09$). The evolution of the spectral parameters during the cooling tail are well-described by the metal-rich atmosphere models. The combined constraints on the NS mass and radius from the direct cooling method and the tidal deformability strongly suggest very high atmospheric abundance of the iron group elements and limit the distance to the source to $11\pm1$ kpc.

We analyzed RXTE/PCA observations of the bursting pulsar GRO J1744-28 during its two outbursts in 1995-1997. We have found a significant transition, i.e., a sharp change, of both spectral and timing properties around a flux of $\rm 10^{-8}\ ergs\ cm^{-2}\ s^{-1}$ (3-30 keV), which corresponds to a luminosity of 2--8$\rm \times 10^{37}\ ergs\ s^{-1}$ for a distance of 4--8 kpc. We define the faint (bright) state when the flux is smaller (larger) than this threshold. In the faint state, the spectral hardness increases significantly with the increasing flux, while remaining almost constant in the bright state. In addition, we find that the pulsed fraction is positively related to both the flux and the energy (<30keV) in both states. Thanks to the very stable pulse profile shape in all energy bands, a hard X-ray lag could be measured. This lag is only significant in the faint state (reaching $\sim$ 20\,ms between 3-4 keV and 16-20 keV), and much smaller in the bright state ($\lesssim$2ms). We speculate that the discovered transition might be caused by changes of the accretion structure on the surface of the neutron star, and the hard X-ray lags are due to a geometry effect.

The hardness of the X-ray spectra of intermediate polars (IPs) is determined mainly by the white dwarf (WD) compactness (mass-radius ratio, M/R) and, thus, hard X-ray spectra can be used to constrain the WD mass. An accurate mass estimate requires the finite size of the WD magnetosphere R_m to be taken into the account. We suggested to derive it either directly from the observed break frequency in power spectrum of X-ray or optical lightcurves of a polar, or assuming the corotation. Here we apply this method to all IPs observed by NuSTAR (10 objects) and Swift/BAT (35 objects). For the dwarf nova GK Per we also observe a change of the break frequency with flux, which allows to constrain the dependence of the magnetosphere radius on the mass-accretion rate. For our analysis we calculated an additional grid of two-parameter (M and R_m/R) model spectra assuming a fixed, tall height of the accretion column H_sh/R=0.25, which is appropriate to determine WD masses in low mass-accretion IPs like EX\,Hya. Using the Gaia Data Release 2 we obtain for the first time reliable estimates of the mass-accretion rate and the magnetic field strength at the WD surface for a large fraction of objects in our sample. We find that most IPs accrete at rate of ~10^{-9} M_Sun/yr, and have magnetic fields in the range 1--10 MG. The resulting WD mass average of our sample is 0.79 +/- 0.16 M_Sun, which is consistent with earlier estimates.

G15.9+0.2 is a Galactic shell-type supernova remnant (SNR), which was detected in radio and has been confirmed in X-rays based on Chandra observations. An X-ray point source CXOUJ181852.0-150213 has been detected and suggested to be an associated neutron star. In a recent study, we have confirmed the source to be a central compact object (CCO). We have studied the SNR using high-resolution X-ray data taken with Chandra in combination with infrared (IR) data in order to understand its emission and to derive its physical parameters. This will also help to constrain, e.g., the age of the CCO and the environment in which it was born. The spectral analysis of the X-ray emission using the new Chandra data and the comparison to the IR data have shown that the SNR is relatively young with an age of a few thousand years and that its emission is dominated by that of shocked interstellar medium (ISM). However, the analysis of the spectrum of the bright eastern shell shows that there is an additional emission component with enhanced abundances of {\alpha} elements and Fe, suggesting ejecta emission. The multi-wavelength emission is consistent with SNR G15.9+0.2 expanding in an ISM with a density gradient, while there is also colder material located in front of the SNR, which absorbs its thermal X-ray emission in the softer bands.

We report on the analysis of XMM-Newton observations of the central compact object CXOU J160103.1-513353 located in the center of the non-thermally emitting supernova remnant (SNR) G330.2+1.0. The X-ray spectrum of the source is well described with either single-component carbon or two-component hydrogen atmosphere models. In the latter case, the observed spectrum is dominated by the emission from a hot component with a temperature ~3.9MK, corresponding to the emission from a hotspot occupying ~1% of the stellar surface (assuming a neutron star with mass M = 1.5M$_\odot$, radius of 12 km, and distance of ~5 kpc as determined for the SNR). The statistics of the spectra and obtained upper limits on the pulsation amplitude expected for a rotating neutron star with hot spots do not allow us to unambiguously distinguish between these two scenarios. We discuss, however, that while the non-detection of the pulsations can be explained by the unfortunate orientation in CXOU J160103.1-513353, this is not the case when the entire sample of similar objects is considered. We therefore conclude that the carbon atmosphere scenario is more plausible.

Neutron stars are the densest, directly observable stellar objects in the universe and serve as unique astrophysical laboratories to study the behavior of matter under extreme physical conditions. This book chapter is devoted to describing how electromagnetic observations, particularly at X-ray, optical and radio wavelengths, can be used to measure the mass and radius of neutron stars and how this leads to constraints on the equation of state of ultra-dense matter. Having accurate theoretical models to describe the astrophysical data is essential in this effort. We will review different methods to constrain neutron star masses and radii, discuss the main observational results and theoretical developments achieved over the past decade, and provide an outlook of how further progress can be made with new and upcoming ground-based and space-based observatories.

The classical limit on the accretion luminosity of a neutron star is given by the Eddington luminosity. The advanced models of accretion onto magnetized neutron stars account for the appearance of magnetically confined accretion columns and allow the accretion luminosity to be higher than the Eddington value by a factor of tens. However, the recent discovery of pulsations from ultraluminous X-ray source (ULX) in NGC 5907 demonstrates that the accretion luminosity can exceed the Eddington value up to by a factor of 500. We propose a model explaining observational properties of ULX-1 in NGC 5907 without any ad hoc assumptions. We show that the accretion column at extreme luminosity becomes advective. Enormous energy release within a small geometrical volume and advection result in very high temperatures at the bottom of accretion column, which demand to account for the energy losses due to neutrino emission which can be even more effective than the radiation energy losses. We show that the total luminosity at the mass accretion rates above $10^{21}\,{\rm g\,s^{-1}}$ is dominated by the neutrino emission similarly to the case of core-collapse supernovae. We argue that the accretion rate measurements based on detected photon luminosity in case of bright ULXs powered by neutron stars can be largely underestimated due to intense neutrino emission. The recently discovered pulsating ULX-1 in galaxy NGC 5907 with photon luminosity of $\sim 10^{41}\,{\rm erg\,s^{-1}}$ is expected to be even brighter in neutrinos and is thus the first known Neutrino Pulsar.

[Abridged] We present more than 4 years of Swift X-ray observations of the 2013 superoutburst, subsequent decline and quiescence of the WZ Sge-type dwarf nova SSS J122221.7-311525 (SSS122222) from 6 days after discovery. Only a handful of WZ Sge-type dwarf novae have been observed in X-rays, and until recently GW Lib was the only binary of this type with complete coverage of an X-ray light curve throughout a superoutburst. We collected extensive X-ray data of a second such system to understand the extent to which the unexpected properties of GW Lib are common to the WZ Sge class. We analysed the X-ray light curve and compared it with the behaviour of superhumps which were detected in the optical light curve. We also performed spectral analysis of the data. The results were compared with the properties of GW Lib, for which new X-ray observations were also obtained. SSS122222 was variable and around five times brighter in 0.3-10 keV X-rays during the superoutburst than in quiescence, mainly because of a significant strengthening of a high-energy component of the X-ray spectrum. The post-outburst decline of the X-ray flux lasted at least 500 d. The data show no evidence of the expected optically thick boundary layer in the system during the outburst. SSS122222 also exhibited a sudden X-ray flux change in the middle of the superoutburst, which occurred exactly at the time of the superhump stage transition. A similar X-ray behaviour was also detected in GW Lib. This result demonstrates a relationship between the outer disc and the white dwarf boundary layer for the first time, and suggests that models for accretion discs in high mass ratio accreting binaries are currently incomplete. The very long decline to X-ray quiescence is also in strong contrast to the expectation of low viscosity in the disc after outburst.

The aim of this paper is to investigate the transition of a strongly magnetized neutron star into the accretion regime with very low accretion rate. For this purpose we monitored the Be-transient X-ray pulsar GRO J1008-57 throughout a full orbital cycle. The current observational campaign was performed with the Swift/XRT telescope in the soft X-ray band (0.5-10 keV) between two subsequent Type I outbursts in January and September 2016. The expected transition to the propeller regime was not observed. However, the transitions between different regimes of accretion were detected. In particular, after an outburst the source entered a stable accretion state characterised by the accretion rate of ~10^14-10^15 g/s. We associate this state with accretion from a cold (low-ionised) disc of temperature below ~6500 K. We argue that a transition to such accretion regime should be observed in all X-ray pulsars with certain combination of the rotation frequency and magnetic field strength. The proposed model of accretion from a cold disc is able to explain several puzzling observational properties of X-ray pulsars.

Observations of thermonuclear X-ray bursts from accreting neutron stars (NSs) in low-mass X-ray binary systems can be used to constrain NS masses and radii. Most previous work of this type has set these constraints using Planck function fits as a proxy: both the models and the data are fit with diluted blackbody functions to yield normalizations and temperatures which are then compared against each other. Here, for the first time, we fit atmosphere models of X-ray bursting NSs directly to the observed spectra. We present a hierarchical Bayesian fitting framework that uses state-of-the-art X-ray bursting NS atmosphere models with realistic opacities and relativistic exact Compton scattering kernels as a model for the surface emission. We test our approach against synthetic data, and find that for data that are well-described by our model we can obtain robust radius, mass, distance, and composition measurements. We then apply our technique to Rossi X-ray Timing Explorer observations of five hard-state X-ray bursts from 4U 1702-429. Our joint fit to all five bursts shows that the theoretical atmosphere models describe the data well but there are still some unmodeled features in the spectrum corresponding to a relative error of 1-5% of the energy flux. After marginalizing over this intrinsic scatter, we find that at 68% credibility the circumferential radius of the NS in 4U 1702-429 is R = 12.4+-0.4 km, the gravitational mass is M=1.9+-0.3 Msun, the distance is 5.1 < D/kpc < 6.2, and the hydrogen mass fraction is X < 0.09.

The ultracompact low-mass X-ray binary 4U 1820-30 situated in the globular cluster NGC 6624 has an orbital period of only $\approx$11.4 min which likely implies a white dwarf companion. The observed X-ray bursts demonstrate a photospheric radius expansion phase and therefore are believed to reach the Eddington luminosity allowing us to estimate the mass and the radius of the neutron star (NS) in this binary. Here we re-analyse all Rossi X-ray Timing Explorer observations of the system and confirm that almost all the bursts took place during the hard persistent state of the system. This allows us to use the recently developed direct cooling tail method to estimate the NS mass and radius. However, because of the very short, about a second, duration of the cooling tail phases that can be described by the theoretical atmosphere models, the obtained constraints on the NS radius are not very strict. Assuming a pure helium NS atmosphere we found that the NS radius is in the range 10-12 km, if the NS mass is below 1.7 $M_\odot$, and in a wider range of 8-12 km for a higher 1.7-2.0 $M_\odot$ NS mass. The method also constrains the distance to the system to be 6.5$\pm$0.5 kpc, which is consistent with the distance to the cluster. For the solar composition atmosphere, the NS parameters are in strong contradiction with the generally accepted range of possible NS masses and radii.

Thermonuclear (type-I) X-ray bursts, observed from neutron star (NS) low-mass X-ray binaries (LMXB), provide constraints on NS masses and radii and consequently the equation of state of NS cores. In such analyses various assumptions are made without knowing if they are justified. We have analyzed X-ray burst spectra from the LMXB 4U 1608-52, with the aim of studying how the different persistent emission components react to the bursts. During some bursts in the soft spectral state we find that there are two variable components; one corresponding to the burst blackbody component and another optically thick Comptonized component. We interpret the latter as the spreading layer between the NS surface and the accretion disc, which is not present during the hard state bursts. We propose that the spectral changes during the soft state bursts are driven by the spreading layer that could cover almost the entire NS in the brightest phases due to the enhanced radiation pressure support provided by the burst, and that the layer subsequently returns to its original state during the burst decay. When deriving the NS mass and radius using the soft state bursts two assumptions are therefore not met: the NS is not entirely visible and the burst emission is reprocessed in the spreading layer, causing distortions of the emitted spectrum. For these reasons the NS mass and radius constraints using the soft state bursts are different compared to the ones derived using the hard state bursts.

The data from all observations of RX J0720.4-3125 conducted by XMM-Newton EPIC-pn with the same instrumental setup in 2000-2012 were reprocessed to form a homogenous data set of solar barycenter corrected photon arrival times registered from RX J0720.4-3125. A Bayesian method for the search, detection, and estimation of the parameters of an unknown-shaped periodic signal was employed as developed by Gregory & Loredo (1992). A number of complex models (single and double peaked) of light curves from pulsating neutron stars were statistically analyzed. The distribution of phases for the registered photons was calculated by folding the arrival times with the derived spin-period and the resulting distribution of phases approximated with a mixed von Mises distribution, and its parameters were estimated by using the Expected Maximization method. Spin phase-resolved spectra were extracted, and a number of highly magnetized atmosphere models of an INS were used to fit simultaneously, the results were verified via an MCMC approach. The phase-folded light curves in different energy bands with high S/N ratio show a high complexity and variations depending on time and energy. They can be parameterized with a mixed von Mises distribution, i.e. with double-peaked light curve profile showing a dependence of the estimated parameters (mean directions, concentrations, and proportion) upon the energy band, indicating that radiation emerges from at least two emitting areas. The genuine spin-period of the isolated neutron star RX J0720-3125 derived as more likely is twice of that reported in the literature (16.78s instead of 8.39s). The gravitational redshift of RX J0720.4-3125 was determined to $z=0.205_{-0.003}^{+0.006}$ and the compactness was estimated to $(M/M_{Sun})/R(km)=0.105 \pm 0.002 $.

We report extensive 3-yr multiwavelength observations of the WZ Sge-type dwarf nova SSS J122221.7-311525 during its unusual double superoutburst, the following decline and in quiescence. The second segment of the superoutburst had a long duration of 33 d and a very gentle decline with a rate of 0.02 mag/d, and it displayed an extended post-outburst decline lasting at least 500 d. Simultaneously with the start of the rapid fading from the superoutburst plateau, the system showed the appearance of a strong near-infrared excess resulting in very red colours, which reached extreme values (B-I~1.4) about 20 d later. The colours then became bluer again, but it took at least 250 d to acquire a stable level. Superhumps were clearly visible in the light curve from our very first time-resolved observations until at least 420 d after the rapid fading from the superoutburst. The spectroscopic and photometric data revealed an orbital period of 109.80 min and a fractional superhump period excess <0.8 per cent, indicating a very low mass ratio q<0.045. With such a small mass ratio the donor mass should be below the hydrogen-burning minimum mass limit. The observed infrared flux in quiescence is indeed much lower than is expected from a cataclysmic variable with a near-main-sequence donor star. This strongly suggests a brown-dwarf-like nature for the donor and that SSS J122221.7-311525 has already evolved away from the period minimum towards longer periods, with the donor now extremely dim.

We report on the analysis of NuSTAR observations of the Be-transient X-ray pulsar V 0332+53 during the giant outburst in 2015 and another minor outburst in 2016. We confirm the cyclotron-line energy-luminosity correlation previously reported in the source and the line energy decrease during the giant outburst. Based on 2016 observations, we find that a year later the line energy has increased again essentially reaching the pre-outburst values. We discuss this behaviour and conclude that it is likely caused by a change of the emission region geometry rather than previously suggested accretion-induced decay of the neutron stars magnetic field. At lower luminosities, we find for the first time a hint of departure from the anticorrelation of line energy with flux, which we interpret as a transition from super- to sub-critical accretion associated with the disappearance of the accretion column. Finally, we confirm and briefly discuss the orbital modulation observed in the outburst light curve of the source.

The X-ray spectra of the neutron stars located in the centers of supernova remnants Cas A and HESS J1731-347 are well fit with carbon atmosphere models. These fits yield plausible neutron star sizes for the known or estimated distances to these supernova remnants. The evidence in favor of the presence of a pure carbon envelope at the neutron star surface is rather indirect and is based on the assumption that the emission is generated uniformly by the entire stellar surface. Although this assumption is supported by the absence of pulsations, the observational upper limit on the pulsed fraction is not very stringent. In an attempt to quantify this evidence, we investigate the possibility that the observed spectrum of the neutron star in HESS J1731-347 is a combination of the spectra produced in a hydrogen atmosphere of the hotspots and of the cooler remaining part of the neutron star surface. The lack of pulsations in this case has to be explained either by a sufficiently small angle between the neutron star spin axis and the line of sight, or by a sufficiently small angular distance between the hotspots and the neutron star rotation poles. As the observed flux from a non-uniformly emitting neutron star depends on the angular distribution of the radiation emerging from the atmosphere, we have computed two new grids of pure carbon and pure hydrogen atmosphere model spectra accounting for Compton scattering. Using new hydrogen models, we have evaluated the probability of a geometry that leads to a pulsed fraction below the observed upper limit to be about 8.2 %. Such a geometry thus seems to be rather improbable but cannot be excluded at this stage.

Magnetized neutron stars power at least some ultra-luminous X-ray sources. The accretion flow in these cases is interrupted at the magnetospheric radius and then reaches the surface of a neutron star following magnetic field lines. Accreting matter moving along magnetic field lines forms the accretion envelope around the central object. We show that, in case of high mass accretion rates $\gtrsim 10^{19}\,{\rm g\,s^{-1}}$ the envelope becomes closed and optically thick, which influences the dynamics of the accretion flow and the observational manifestation of the neutron star hidden behind the envelope. Particularly, the optically thick accretion envelope results in a multi-color black-body spectrum originating from the magnetospheric surface. The spectrum and photon energy flux vary with the viewing angle, which gives rise to pulsations characterized by high pulsed fraction and typically smooth pulse profiles. The reprocessing of radiation due to interaction with the envelope leads to the disappearance of cyclotron scattering features from the spectrum. We speculate that the super-orbital variability of ultra-luminous X-ray sources powered by accreting neutron stars can be attributed to precession of the neutron star due to interaction of magnetic dipole with the accretion disc.

Determining neutron star (NS) radii and masses can help to understand the properties of matter at supra-nuclear densities. Thermal emission during thermonuclear X-ray bursts from NSs in low-mass X-ray binaries provides a unique opportunity to study NS parameters, because of the high fluxes, large luminosity variations and the related changes in the spectral properties. The standard cooling tail method uses hot NS atmosphere models to convert the observed spectral evolution during cooling stages of X-ray bursts to the Eddington flux F_Edd and the stellar angular size \Omega. These are then translated to the constraints on the NS mass M and radius R. Here we present the improved, direct cooling tail method that generalises the standard approach. First, we adjust the cooling tail method to account for the bolometric correction to the flux. Then, we fit the observed dependence of the blackbody normalization on flux with a theoretical model directly on the M-R plane by interpolating theoretical dependences to a given gravity, hence ensuring only weakly informative priors for M and R instead of F_Edd and \Omega. The direct cooling method is demonstrated using a photospheric radius expansion burst from SAX J1810.8--2609, which has happened when the system was in the hard state. Comparing to the standard cooling tail method, the confidence regions are shifted by 1 sigma towards larger radii, giving R=11.5-13.0 km at M=1.3-1.8 M_sun for this NS.

We present a study of the evolution of the optical and X-ray fluxes during outbursts of two short-period cataclysmic variables, the confirmed intermediate polar CC Scl and the intermediate polar candidate FS Aur. We found that the X-ray and optical light curves are well correlated in both objects, although the amplitudes of outbursts in X-rays are smaller than those in the optical. The ratio of the outburst amplitudes in X-rays and the optical in both objects is close to ~0.6. This is significantly higher than was observed during the outburst of the non-magnetic dwarf nova U Gem, in which this ratio was only ~0.03. The obtained data also suggest that the dependence between the X-ray and optical fluxes must steepen significantly toward very low accretion rates and very low fluxes. Similarities in the behaviour of CC Scl and FS Aur indicate strongly the magnetic nature of the white dwarf in FS Aur.

When neutron stars (NS) accrete gas from low-mass binary companions, explosive nuclear burning reactions in the NS envelope fuse hydrogen and helium into heavier elements. The resulting thermonuclear (type-I) X-ray bursts produce energy spectra that are fit well with black bodies, but a significant number of burst observations show deviations from Planck spectra. Here we present our analysis of RXTE/PCA observations of X-ray bursts from the NS low-mass X-ray binary HETE J1900.1-2455. We have discovered that the non-Planckian spectra are caused by photo-ionization edges. The anti-correlation between the strength of the edges and the colour temperature suggests that the edges are produced by the nuclear burning ashes that have been transported upwards by convection and become exposed at the photosphere. The atmosphere model fits show that occasionally the photosphere can consist entirely of metals, and that the peculiar changes in black body temperature and radius can be attributed to the emergence and disappearance of metals in the photosphere. As the metals are detected already in the Eddington-limited phase, it is possible that a radiatively driven wind ejects some of the burning ashes into the interstellar space.

We present our study of the central point source CXOU J181852.0-150213 in the young Galactic supernova remnant (SNR) G15.9+0.2 based on the recent ~90 ks Chandra observations. The point source was discovered in 2005 in shorter Chandra observations and was hypothesized to be a neutron star associated with the SNR. Our X-ray spectral analysis strongly supports the hypothesis of a thermally emitting neutron star associated with G15.9+0.2. We conclude that the object belongs to the class of young cooling low-magnetized neutron stars referred to as central compact objects (CCOs). We modeled the spectrum of the neutron star with a blackbody spectral function and with our hydrogen and carbon neutron star atmosphere models, assuming that the radiation is uniformly emitted by the entire stellar surface. Under this assumption, only the carbon atmosphere models yield a distance that is compatible with a source located in the Galaxy. In this respect, CXOU J181852.0-150213 is similar to two other well-studied CCOs, the neutron stars in Cas A and in HESS J1731-347, for which carbon atmosphere models were used to reconcile their emission with the known or estimated distances.

We present the results of the monitoring programmes performed with the Swift/XRT telescope and aimed specifically to detect an abrupt decrease of the observed flux associated with a transition to the propeller regime in two well known X-ray pulsars 4U 0115+63 and V 0332+53 during their giant outbursts in 2015. Such transitions were detected at the threshold luminosities of $(1.4\pm0.4)\times10^{36}$ erg s$^{-1}$ and $(2.0\pm0.4)\times10^{36}$ erg s$^{-1}$ for 4U 0115+63 and V 0332+53, respectively. Spectra of the sources are shown to be significantly softer during the low state. In both sources, the accretion at rates close to the aforementioned threshold values briefly resumes during the periastron passage following the transition into propeller regime. The strength of the dipole component of the magnetic field required to inhibit the accretion agrees well with estimates based on the position of the cyclotron lines in their spectra, thus excluding presence of a strong multipole component of the magnetic field in the vicinity of the neutron star.

Symbiotic stars are a heterogeneous class of interacting binaries. Among them, RT Cru has been classified as prototype of a subclass that is characterised by hard X-ray spectra extending past ~20 keV. We analyse ~8.6 Ms of archival INTEGRAL data collected in the period 2003-2014, ~140 ks of Swift/XRT data, and a Suzaku observation of 39 ks, to study the spectral X-ray emission and investigate the nature of the compact object. Based on the 2MASS photometry, we estimate the distance to the source of 1.2-2.4 kpc. The X-ray spectrum obtained with Swift/XRT, JEM-X, IBIS/ISGRI, and Suzaku data is well fitted by a cooling flow model modified by an absorber that fully covers the source and two partial covering absorbers. Assuming that the hard X-ray emission of RT Cru originates from an optically thin boundary layer around a non-magnetic white dwarf, we estimated a mass of the WD of about 1.2 M_Sun. The mass accretion rate obtained for this source might be too high for the optically thin boundary layer scenario. Therefore we investigate other plausible scenarios to model its hard X-ray emission. We show that, alternatively, the observed X-ray spectrum can be explained with the X-ray emission from the post-shock region above the polar caps of a magnetised white dwarf with mass ~0.9-1.1 M_Sun.

Observed hard X-ray spectra of intermediate polars are determined mainly by the accretion flow velocity at the white dwarf surface, which is normally close to the free-fall velocity. This allows to estimate the white dwarf masses as the white dwarf mass-radius relation M-R and the expected free-fall velocities at the surface are well known. This method is widely used, however, derived white dwarf masses M can be systematically underestimated because the accretion flow is stopped at and re-accelerates from the magnetospheric boundary R_m, and therefore, its velocity at the surface will be lower than free-fall.To avoid this problem we computed a two-parameter set of model hard X-ray spectra, which allows to constrain a degenerate M - R_m dependence. On the other hand, previous works showed that power spectra of accreting X-ray pulsars and intermediate polars exhibit breaks at the frequencies corresponding to the Keplerian frequencies at the magnetospheric boundary. Therefore, the break frequency \nu_b in an intermediate polar power spectrum gives another relation in the M - R_m plane. The intersection of the two dependences allows, therefore, to determine simultaneously the white dwarf mass and the magnetospheric radius. To verify the method we analyzed the archival Suzaku observation of EX Hya obtaining M /M_sun= 0.73 \pm 0.06 and R_ m / R = 2.6 \pm 0.4 consistent with the values determined by other authors. Subsequently, we applied the same method to a recent NuSTAR observation of another intermediate polar GK~Per performed during an outburst and found M/M_sun = 0.86 \pm 0.02 and R_ m / R = 2.8 \pm 0.2. The long duration observations of GK Per in quiescence performed by Swift/BAT and INTEGRAL observatories indicate increase of magnetosphere radius R_m at lower accretion rates.

The cooling phase of thermonuclear (type-I) X-ray bursts can be used to constrain the neutron star (NS) compactness by comparing the observed cooling tracks of bursts to accurate theoretical atmosphere model calculations. By applying the so-called cooling tail method, where the information from the whole cooling track is used, we constrain the mass, radius, and distance for three different NSs in low-mass X-ray binaries 4U 1702-429, 4U 1724-307, and SAX J1810.8-260. Care is taken to only use the hard state bursts where it is thought that only the NS surface alone is emitting. We then utilize a Markov chain Monte Carlo algorithm within a Bayesian framework to obtain a parameterized equation of state (EoS) of cold dense matter from our initial mass and radius constraints. This allows us to set limits on various nuclear parameters and to constrain an empirical pressure-density relation for the dense matter. Our predicted EoS results in NS radius between 10.5-12.8 km (95% confidence limits) for a mass of 1.4 $M_{\odot}$. Due to systematic errors and uncertainty in the composition these results should be interpreted as lower limits for the radius.

Central compact objects (CCOs) are thought to be young thermally emitting isolated neutron stars that were born during the preceding core-collapse supernova explosion. Here we present evidence that at least in one case the CCO could have been formed within a binary system. We show that the highly reddened optical source IRAS~17287$-$3443, located $25^{\prime \prime}$ away from the CCO candidate XMMUJ173203.3$-$344518 and classified previously as a post asymptotic giant branch star, is indeed surrounded by a dust shell. This shell is heated by the central star to temperatures of $\sim90$\,K and observed as extended infrared emission in 8-160\,$\mu$m band. The dust temperature also increases in the vicinity of the CCO which implies that it likely resides within the shell. We estimate the total dust mass to be $\sim0.4-1.5\,M_\odot$ which significantly exceeds expected dust yields by normal stars and thus likely condensed from supernova ejecta. Taking into account that both the age of the supernova remnant and the duration of active mass loss phase by the optical star are much shorter than the total lifetime of either object, the supernova and the onset of the active mass loss phase of the companion have likely occurred approximately simultaneously. This is most easily explained if the evolution of both objects is interconnected. We conclude, therefore, that both stars were likely members of the same binary system disrupted by a supernova.

The goal of this work is to develop a simple analytic description of the emission properties (spectrum and polarization) of the condensed, strongly magnetized surface of neutron stars. We have improved the method of van Adelsberg et al. (2005) (arXiv:astro-ph/0406001) for calculating the spectral properties of condensed magnetized surfaces. Using the improved method, we calculate the reflectivity of an iron surface at magnetic field strengths B \sim (10^{12} - 10^{14}) G, with various inclinations of the magnetic field lines and radiation beam with respect to the surface and each other. We construct analytic expressions for the emissivity of this surface as functions of the photon energy, magnetic field strength, and the three angles that determine the geometry of the local problem. Using these expressions, we calculate X-ray spectra for neutron stars with condensed iron surfaces covered by thin partially ionized hydrogen atmospheres. We develop simple analytic descriptions of the intensity and polarization of radiation emitted or reflected by condensed iron surfaces of neutron stars with strong magnetic fields typical for isolated neutron stars. This description provides boundary conditions at the bottom of a thin atmosphere, which are more accurate than previously used approximations. The spectra calculated with this improvement show absorption features different from those in simplified models. The approach developed in this paper yields results that can facilitate modeling and interpretation of the X-ray spectra of isolated, strongly magnetized, thermally emitting neutron stars.

Magnetic fields are important for accretion disc structure. Magnetic fields in a disc system may be transported with the accreted matter. They can be associated with either the central body and/or jet, and be fossil or dynamo excited in situ. We consider dynamo excitation of magnetic fields in accretion discs of accreting binary systems in an attempt to clarify possible configurations of dynamo generated magnetic fields. We first model the entire disc with realistic radial extent and thickness using an alpha-quenching non-linearity. We then study the simultaneous effect of feedback from the Lorentz force from the dynamo-generated field. We perform numerical simulations in the framework of a relatively simple mean-field model which allows the generation of global magnetic configurations. We explore a range of possibilities for the dynamo number, and find quadrupolar-type solutions with irregular temporal oscillations that might be compared to observed rapid luminosity fluctuations. The dipolar symmetry models with $R_\alpha<0$ have lobes of strong toroidal field adjacent to the rotation axis that could be relevant to jet launching phenomena. We have explored and extended the solutions known for thin accretion discs.

We present here the first convincing observational manifestation of a magnetar-like magnetic field in an accreting neutron star in binary system - the first pulsating ultra-luminous X-ray source X-2 in the galaxy M82. Using the Chandra X-ray observatory data we show that the source exhibit the bimodal distribution of the luminosity with two well-defined peaks separated by a factor of 40. This behaviour can be interpreted as the action of the "propeller regime" of accretion. The onset of the propeller in a 1.37 s pulsar at luminosity of ~$10^{40}$ erg/s implies the dipole component of the neutron star magnetic field of ~$10^{14}$ G.

We study properties of luminous X-ray pulsars using a simplified model of the accretion column. The maximal possible luminosity is calculated as a function of the neutron star (NS) magnetic field and spin period. It is shown that the luminosity can reach values of the order of $10^{40}\,{\rm erg/s}$ for the magnetar-like magnetic field ($B\gtrsim 10^{14}\,{\rm G}$) and long spin periods ($P\gtrsim 1.5\,{\rm s}$). The relative narrowness of an area of feasible NS parameters which are able to provide higher luminosities leads to the conclusion that $L\simeq 10^{40}\,\,{\rm erg/s}$ is a good estimate for the limiting accretion luminosity of a NS. Because this luminosity coincides with the cut-off observed in the high mass X-ray binaries luminosity function which otherwise does not show any features at lower luminosities, we can conclude that a substantial part of ultra-luminous X-ray sources are accreting neutron stars in binary systems.

Cyclotron resonance scattering features observed in the spectra of some X-ray pulsars show significant changes of the line centroid energy with the pulsar luminosity. Whereas for bright sources above the so called critical luminosity these variations are established to be connected with the appearance of the high accretion column above the neutron star surface, at low, sub-critical luminosities the nature of the variations (but with the opposite sign) has not been discussed widely. We argue here that the cyclotron line is formed when the radiation from a hotspot propagates through the plasma falling with a mildly relativistic velocity onto the neutron star surface. The position of the cyclotron resonance is determined by the Doppler effect. The change of the cyclotron line position in the spectrum with luminosity is caused by variations of the velocity profile in the line-forming region affected by the radiation pressure force. The presented model has several characteristic features: (i) the line centroid energy is positively correlated with the luminosity; (ii) the line width is positively correlated with the luminosity as well; (iii) the position and the width of the cyclotron absorption line are variable over the pulse phase; (iv) the line has a more complicated shape than widely used Lorentzian or Gaussian profiles; (v) the phase-resolved cyclotron line centroid energy and the width are negatively and positively correlated with the pulse intensity, respectively. The predictions of the proposed theory are compared with the variations of the cyclotron line parameters in the X-ray pulsar GX 304-1 over a wide range of sub-critical luminosities as seen by the INTEGRAL observatory.

Superbursts are hours-long X-ray flares attributed to the thermonuclear runaway burning of carbon-rich material in the envelope of accreting neutron stars. By studying the details of the X-ray light curve, properties of carbon combustion can be determined. In particular, we show that the shape of the rise of the light curve is set by the the slope of the temperature profile left behind by the carbon flame. We analyse RXTE/PCA observations of 4U 1636-536 and separate the direct neutron star emission from evolving photoionized reflection and persistent spectral components. This procedure results in the highest quality light curve ever produced for the superburst rise and peak, and interesting behaviour is found in the tail. The rising light curve between 100 and 1000 seconds is inconsistent with the idea that the fuel burned locally and instantaneously everywhere, as assumed in some previous models. By fitting improved cooling models, we measure for the first time the radial temperature profile of the superbursting layer. We find $d\ln T/d\ln P=1/4$. Furthermore, 20% of the fuel may be left unburned. This gives a new constraint on models of carbon burning and propagation in superbursts.

Model spectra of neutron star atmospheres are nowadays widely used to fit the observed thermal X-ray spectra of neutron stars. This fitting is the key element in the method of the neutronstar radius determination. Here, we present the basic assumptions used for the neutron star atmosphere modeling as well as the main qualitative features of the stellar atmospheres leading to the deviations of the emergent model spectrum from blackbody. We describe the properties of two of our model atmosphere grids: (i) pure carbon atmospheres for relatively cool neutron stars (1--4 MK) and (ii) hot atmospheres with Compton scattering taken into account. The results obtained by applying these grids to model the X-ray spectra of the central compact object in supernova remnant HESS 1731-347, and two X-ray bursting neutron stars in low-mass X-ray binaries, 4U 1724-307 and 4U 1608-52, are presented. Possible systematic uncertainties associated with the obtained neutron star radii are discussed.

The central compact object in the supernova remnant HESS J1731-347 appears to be the hottest observed isolated cooling neutron star. The cooling theory of neutron stars enables one to explain observations of this star by assuming the presence of strong proton superfluidity in the stellar core and the existence of the surface heat blanketing envelope which almost fully consists of carbon. The cooling model of this star is elaborated to take proper account of the neutrino emission due to neutron-neutron collisions which is not suppressed by proton superfluidity. Using the results of spectral fits of observed thermal spectra for the distance of 3.2 kpc and the cooling theory for the neutron star of age 27 kyr, new constraints on the stellar mass and radius are obtained which are more stringent than those derived from the spectral fits alone.

Low-mass X-ray binaries hosting neutron stars (NS) exhibit thermonuclear (type-I) X-ray bursts, which are powered by unstable nuclear burning of helium and/or hydrogen into heavier elements deep in the NS "ocean". In some cases the burning ashes may rise from the burning depths up to the NS photosphere by convection, leading to the appearance of the metal absorption edges in the spectra, which then force the emergent X-ray burst spectra to shift toward lower energies. These effects may have a substantial impact on the color correction factor $f_c$ and the dilution factor $w$, the parameters of the diluted blackbody model $F_E \approx w B_E(f_c T_{eff})$ that is commonly used to describe the emergent spectra from NSs. The aim of this paper is to quantify how much the metal enrichment can change these factors. We have developed a new NS atmosphere modeling code, which has a few important improvements compared to our previous code required by inclusion of the metals. The opacities and the internal partition functions (used in the ionization fraction calculations) are now taken into account for all atomic species. In addition, the code is now parallelized to counter the increased computational load. We compute a detailed grid of atmosphere models with different exotic chemical compositions that mimic the presence of the burning ashes. From the emerging model spectra we compute the color correction factors $f_c$ and the dilution factors $w$ that can then be compared to the observations. We find that the metals may change $f_c$ by up to about 40%, which is enough to explain the scatter seen in the blackbody radius measurements. The presented models open up the possibility for determining NS mass and radii more accurately, and may also act as a tool to probe the nuclear burning mechanisms of X-ray bursts.

We use the best available X-ray data from the intermediate polar EX Hydrae to study the cooling-flow model often applied to interpret the X-ray spectra of these accreting magnetic white dwarf binaries. First, we resolve a long-standing discrepancy between the X-ray and optical determinations of the mass of the white dwarf in EX Hya by applying new models of the inner disk truncation radius. Our fits to the X-ray spectrum now agree with the white dwarf mass of 0.79 M$_{\odot}$sun determined using dynamical methods through spectroscopic observations of the secondary. We use a simple isobaric cooling flow model to derive the emission line fluxes, emission measure distribution, and H-like to He-like line ratios for comparison with the 496 ks Chandra High Energy Transmission Grating observation of EX Hydrae. We find that the H/He ratios are not well reproduced by this simple isobaric cooling flow model and show that while H-like line fluxes can be accurately predicted, fluxes of lower-Z He-like lines are significantly underestimated. This discrepancy suggests that some extra heating mechanism plays an important role at the base of the accretion column, where cooler ions form. We thus explored more complex cooling models including the change of gravitational potential with height in the accretion column and a magnetic dipole geometry. None of these modifications to the standard cooling flow model are able to reproduce the observed line ratios. While a cooling flow model with subsolar (0.1 $\odot$) abundances is able to reproduce the line ratios by reducing the cooling rate at temperatures lower than $\sim 10^{7.3}$ K, the predicted line-to-continuum ratios are much lower than observed. We discuss and discard mechanisms such as photoionization, departures from constant pressure, resonant scattering, different electron-ion temperatures, and Compton cooling. [Abridged]

We present an analysis of BeppoSAX observations of the unique transient bursting X-ray pulsar GRO J1744-28. The observations took place in March 1997 during the decay phase of the outburst. We find that the persistent broadband X-ray continuum of the source is consistent with a cutoff power law typical for the accreting pulsars. We also detect the fluorescence iron line at 6.7 keV and an absorption feature at ~4.5 keV, which we interpret as a cyclotron line. The corresponding magnetic field strength in the line forming region is ~3.7 x 10^11 G. Neither line is detected in the spectra of the bursts. However, additional soft thermal component with kT ~2 keV was required to describe the burst spectrum. We briefly discuss the nature of this component and argue that among other possibilities it might be connected with thermonuclear flashes at the neutron star surface which accompany the accretion-powered bursts in the source.

We present the results of the pulse phase- and luminosity-resolved spectroscopy of the transient X-ray pulsar V0332+53, performed for the first time in a wide luminosity range (1-40)x10^{37} erg/s during a giant outburst observed by the RXTE observatory in Dec 2004 - Feb 2005. We characterize the spectra quantitatively and built the detailed "three-dimensional" picture of spectral variations with pulse phase and throughout the outburst. We show that all spectral parameters are strongly variable with the pulse phase, and the pattern of this variability significantly changes with luminosity directly reflecting the associated changes in the structure of emission regions and their beam patterns. Obtained results are qualitatively discussed in terms of the recently developed reflection model for the formation of cyclotron lines in the spectra of X-ray pulsars.

This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of thermonuclear X-ray bursts on accreting neutron stars. For a summary, we refer to the paper.

The accretion flow around X-ray pulsars with a strong magnetic field is funnelled by the field to relatively small regions close to the magnetic poles of the neutron star (NS), the hotspots. During strong outbursts regularly observed from some X-ray pulsars, the X-ray luminosity can be so high, that the emerging radiation is able to stop the accreting matter above the surface via radiation-dominated shock, and the accretion column begins to rise. This border luminosity is usually called the "critical luminosity". Here we calculate the critical luminosity as a function of the NS magnetic field strength $B$ using exact Compton scattering cross section in strong magnetic field. Influence of the resonant scattering and photon polarization is taken into account for the first time. We show that the critical luminosity is not a monotonic function of the B-field. It reaches a minimum of a few 10^{36} erg s^{-1} when the cyclotron energy is about 10 keV and a considerable amount of photons from a hotspot have energy close to the cyclotron resonance. For small B, this luminosity is about 10^{37} erg s^{-1}, nearly independent of the parameters. It grows for the B-field in excess of 10^{12} G because of the drop in the effective cross-section of interaction below the cyclotron energy. We investigate how different types of the accretion flow and geometries of the accretion channel affect the results and demonstrate that the general behaviour of the critical luminosity on B-field is very robust. The obtained results are shown to be in a good agreement with the available observational data and provide a necessary ground for the interpretation of upcoming high quality data from the currently operating and planned X-ray telescopes.

Context: Central Compact Objects (CCOs) in supernova remnants are isolated thermally emitting neutron stars (NSs). They are most probably characterized by a magnetic field strength that is roughly two orders of magnitude lower than that of most of the radio and accreting pulsars. The thermal emission of CCOs can be modeled to obtain constraints on the physical parameters of the star such as its mass, radius, effective temperature, and chemical composition. Aims: The CCO in HESS, J1731-347 is one of the brightest objects of this class. We present our analysis of two new XMM-Newton observations of the source which increase the total exposure time of the by a factor of five compared to the analyses presented before. Methods: We use our numerical spectral models for carbon and hydrogen atmospheres to fit the spectrum of the CCO. From our fits, we derive constraints on the physical parameters of the emitting star. We also use the new data to derive new upper limits on pulsations and to confirm the absence of a long-term variability. Results: The analysis shows that atmosphere spectral models are clearly preferred by the fit over the blackbody function. Under the assumption that the X-ray emission is uniformly produced by the entire star surface, hydrogen atmosphere models lead to uncomfortably large distances of the CCO, above 7-8 kpc. On the other hand, the carbon atmosphere model formally excludes distances above 5-6 kpc and is compatible with the source located in the Scutum-Crux (3 kpc) or Norma-Cygnus (4.5 kpc) Galactic spiral arm. We provide and discuss the corresponding confidence contours in the NS mass-radius plane. The measured effective temperature indicates that the NS is exceptionally hot for the estimated age of 30 kyr. We discuss possible cooling scenarios to explain this property, as well as possible additional constraints on the star mass and radius from cooling theory.

Neutron star (NS) masses and radii can be estimated from observations of photospheric radius-expansion X-ray bursts, provided the chemical composition of the photosphere, the spectral colour-correction factors in the observed luminosity range, and the emission area during the bursts are known. By analysing 246 X-ray bursts observed by the Rossi X-ray Timing Explorer from 11 low-mass X-ray binaries, we find a dependence between the persistent spectral properties and the time evolution of the black body normalisation during the bursts. All NS atmosphere models predict that the colour-correction factor decreases in the early cooling phase when the luminosity first drops below the limiting Eddington value, leading to a characteristic pattern of variability in the measured blackbody normalisation. However, the model predictions agree with the observations for most bursts occurring in hard, low-luminosity, 'island' spectral states, but rarely during soft, high-luminosity, 'banana' states. The observed behaviour may be attributed to the accretion flow, which influences cooling of the NS preferentially during the soft state bursts. This result implies that only the bursts occurring in the hard, low-luminosity spectral states can be reliably used for NS mass and radius determination.

Here we compute detailed model spectra of recently published optically thick one-dimensional radial baundary layer (BL) models in cataclysmic variables and compare them with observed soft X-ray/extreme ultraviolet (EUV) spectra of dwarf novae in outburst. Every considered BL model is divided into a number of rings, and for each ring, a structure model along the vertical direction is computed using the stellar-atmosphere method. The ring spectra are then combined into a BL spectrum taking Doppler broadening and limb darkening into account. Two sets of model BL spectra are computed, the first of them consists of BL models with fixed white dwarf (WD) mass (1 M_sun) and various relative WD angular velocities (0.2, 0.4, 0.6 and 0.8 break-up velocities), while the other deals with a fixed relative angular velocity (0.8 break-up velocity) and various WD masses (0.8, 1, and 1.2 M_sun). The model spectra show broad absorption features because of blending of numerous absorption lines, and emission-like features at spectral regions with only a few strong absorption lines. The model spectra are very similar to observed soft X-ray/EUV spectra of SS Cyg and U Gem in outburst. The observed SS Cyg spectrum could be fitted by BL model spectra with WD masses 0.8 - 1 M_sun and relative angular velocities 0.6 - 0.8 break up velocities. These BL models also reproduce the observed ratio of BL luminosity and disk luminosity. The difference between the observed and the BL model spectra is similar to a hot optically thin plasma spectrum and could be associated with the spectrum of outflowing plasma with a mass loss rate compatible with the BL mass accretion rate. The suggested method of computing BL spectra seems very promising and can be applied to other BL models for comparison with EUV spectra of dwarf novae in outburst.

Spectral measurements of thermonuclear (type-I) X-ray bursts from low mass X-ray binaries have been used to measure neutron star (NS) masses and radii. A number of systematic issues affect such measurements and have raised concerns as to the robustness of the methods. We present analysis of the X-ray emission from bursts observed from 4U 1608-52 at various persistent fluxes. We find a strong dependence of the burst properties on the flux and spectral hardness of the persistent emission before burst. Bursts occurring during the low-accretion rate (hard) state exhibit evolution of the black body normalisation consistent with the theoretical predictions of NS atmosphere models. However, bursts occurring during the high-accretion rate (soft) state show roughly constant normalisation, which is inconsistent with the NS atmosphere models and therefore these bursts cannot be easily used to determine NS parameters. We analyse the hard-state burst to put the lower limit on the neutron star radius in 4U 1608-52 of 13 km (for masses 1.2-2.4 M$_\odot$). The best agreement with the theoretical NS mass-radius relations is achieved for source distances in the range 3.1-3.7 kpc. We expect that the radius limit will be 10 per cent lower if spectral models including rapid rotation are used instead.

We investigate the relation between the optical (g-band) and X-ray (0.5-10 keV) luminosities of accreting nonmagnetic white dwarfs. According to the present-day counts of the populations of star systems in our Galaxy, these systems have the highest space density among the close binary systems with white dwarfs. We show that the dependence of the optical luminosity of accreting white dwarfs on their X-ray luminosity forms a fairly narrow one-parameter curve. The typical half-width of this curve does not exceed 0.2-0.3 dex in optical and X-ray luminosities, which is essentially consistent with the amplitude of the aperiodic flux variability for these objects. At X-ray luminosities Lx~1e32 erg/sec or lower, the optical g-band luminosity of the accretion flow is shown to be related to its X-ray luminosity by a factor ~2-3. At even lower X-ray luminosities (Lx~1e30 erg/sec), the contribution from the photosphere of the white dwarf begins to dominate in the optical spectrum of the binary system and its optical brightness does not drop below Mg~13-14. Using the latter fact, we show that in current and planned X-ray sky surveys, the family of accreting nonmagnetic white dwarfs can be completely identified to the distance determined by the sensitivity of an optical sky survey in this region. For the Sloan Digital Sky Survey (SDSS) with a limiting sensitivity m_g~22.5, this distance is ~400-600 pc

Disk accretion onto stars leads to the formation of a Boundary Layer (BL) near the stellar surface where the disk makes contact with the star. Albeit a large fraction of the total luminosity of the system originates from this tiny layer connecting the accretion disk and the accreting object, its structure has not been fully understood yet. It is the aim of this work, to obtain more insight into the Boundary Layer around the white dwarf in compact binary systems. There are still many uncertainties concerning the extent and temperature of the BL and the rotation rate of the white dwarf. We perform numerical hydrodynamical simulations, where the problem is treated in a one-dimensional, radial approximation (slim disk). The turbulence is described by the alpha parameter viscosity. We include both cooling from the disk surfaces and radial radiation transport. The radiation energy is treated in the one-temperature approximation. For a given M_dot our results show a strong dependence on the stellar mass and rotation rate. Both, the midplane and the effective temperature rise considerably with increasing stellar mass or decreasing stellar rotation rate. Our simulations further show, that the radiation energy and pressure are indeed important in the BL. However, some models show a low optical depth in the BL, making it necessary to find a better representation for optically thin regions. The combination of a high mass and a small radius, such as in white dwarfs, can lead to an enormous energy release in the BL, provided the WD rotates slowly. Since the radial extent of BLs is typically very small (about 0.02 to 0.05 R_star), this leads to surface temperatures of a few hundred thousand Kelvin. All of our models showed subsonic infall velocities with Mach numbers of < 0.4 at most.

The accuracy of measuring the basic parameters of neutron stars is limited in particular by uncertainties in chemical composition of their atmospheres. For example, atmospheres of thermally - emitting neutron stars in supernova remnants might have exotic chemical compositions, and for one of them, the neutron star in CasA, a pure carbon atmosphere has recently been suggested by Ho & Heinke (2009). To test such a composition for other similar sources, a publicly available detailed grid of carbon model atmosphere spectra is needed. We have computed such a grid using the standard LTE approximation and assuming that the magnetic field does not exceed 10^8 G. The opacities and pressure ionization effects are calculated using the Opacity Project approach. We describe the properties of our models and investigate the impact of the adopted assumptions and approximations on the emergent spectra.

To examine the previously claimed fast cooling of the Central Compact Object (CCO) in the Cas A supernova remnant (SNR), we analyzed two Chandra observations of this CCO, taken in a setup minimizing instrumental spectral distortions. We fit the two CCO X-ray spectra from 2006 and 2012 with hydrogen and carbon neutron star atmosphere models. The temperature and flux changes in the 5.5 years between the two epochs depend on the adopted constraints on the fitting parameters and the uncertainties of the effective area calibrations. If we allow a change of the equivalent emitting region size, R_Em, the effective temperature remains essentially the same. If R_Em is held constant, the best-fit temperature change is negative, but its statistical significance ranges from 0.8sigma to 2.5sigma, depending on the model. If we assume that the optical depth of the ACIS filter contaminant in 2012 was +/-10% different from its default calibration value, the significance of the temperature drop becomes 0.8sigma to 3.1sigma, for the carbon atmospheres with constant R_Em. Thus, we do not see a statistically significant temperature drop in our data, but the involved uncertainties are too large to firmly exclude the previously reported fast cooling. Our analysis indicate a decrease of 4%-6% (1.9-2.9sigma significance) for the absorbed flux in the energy range 0.6-6keV between 2006 and 2012, most prominent in the 1.4-1.8 keV energy range. It could be caused by unaccounted changes of the detector response or contributions from unresolved SNR material along the line of sight to the CCO.

Cyclotron resonance scattering features are observed in the spectra of some X-ray pulsars and show significant changes in the line energy with the pulsar luminosity. In a case of bright sources, the line centroid energy is anti-correlated with the luminosity. Such a behaviour is often associated with the onset and growth of the accretion column, which is believed to be the origin of the observed emission and the cyclotron lines. However, this scenario inevitably implies large gradient of the magnetic field strength within the line-forming region, and it makes the formation of the observed line-like features problematic. Moreover, the observed variation of the cyclotron line energy is much smaller than could be anticipated for the corresponding luminosity changes. We argue that a more physically realistic situation is that the cyclotron line forms when the radiation emitted by the accretion column is reflected from the neutron star surface. The idea is based on the facts that a substantial part of column luminosity is intercepted by the neutron star surface and the reflected radiation should contain absorption features. The reflection model is developed and applied to explain the observed variations of the cyclotron line energy in a bright X-ray pulsar V 0332+53 over a wide range of luminosities.

Cyclotron resonance scattering features observed in the spectra of some X-ray pulsars show significant changes of the line energy with the pulsar luminosity. At high luminosities, these variations are often associated with the onset and growth of the accretion column, which is believed to be the origin of the observed emission and of the cyclotron lines. However, this scenario inevitably implies large gradient of the magnetic field strength within the line-forming region, which makes the formation of the observed line-like features problematic. Moreover, the observed variation of the cyclotron line energy is much smaller than could be anticipated for the corresponding luminosity changes. We argue here that a more physically realistic situation is that the cyclotron line forms when the radiation emitted by the accretion column is reflected from the neutron star surface, where the gradient of the magnetic field strength is significantly smaller. We develop here the reflection model and apply it to explain the observed variations of the cyclotron line energy in a bright X-ray pulsar V 0332+53 over a wide range of luminosities.

Context: The CCO candidate in the center of the supernova remnant shell HESS J1731-347 / G353.6-0.7 shows no pulsations and exhibits a blackbody-like X-ray spectrum. If the absence of pulsations is interpreted as evidence for the emitting surface area being the entire neutron star surface, the assumption of the measured flux being due to a blackbody emission translates into a source distance that is inconsistent with current estimates of the remnant's distance. Aims: With the best available observational data, we extended the pulse period search down to a sub-millisecond time scale and used a carbon atmosphere model to describe the X-ray spectrum of the CCO and to estimate geometrical parameters of the neutron star. Methods: To search for pulsations we used data of an observation of the source with XMM-Newton performed in timing mode. For the spectral analysis, we used earlier XMM-Newton observations performed in imaging mode, which permits a more accurate treatment of the background. The carbon atmosphere models used to fit the CCO spectrum are computed assuming hydrostatic and radiative equilibria and take into account pressure ionization and the presence of spectral lines. Results: Our timing analysis did not reveal any pulsations with a pulsed fraction above ~8% down to 0.2 ms. This finding further supports the hypothesis that the emitting surface area is the entire neutron star surface. The carbon atmosphere model provides a good fit to the CCO spectrum and leads to a normalization consistent with the available distance estimates of the remnant. The derived constraints on the mass and radius of the source are consistent with reasonable values of the neutron star mass and radius. After the CCO in Cas A, the CCO in HESS J1731-347 / G353.6-0.7 is the second object of this class for which a carbon atmosphere model provides a consistent description of X-ray emission.

Emission of the neutron star surface potentially contains information about its size and thus of vital importance for high energy astrophysics. In spite of the wealth of data on the emission of luminous accreting neutron stars, the emission of their surfaces is hard to disentangle from their time averaged spectra. A recent X-ray transient source XTE J1701-462 has provided a unique dataset covering the largest ever observed luminosity range for a single source. In this paper, we extract the spectrum of the boundary layer between the inner part of the accretion disc and the neutron star surface with the help of maximally spectral model-independent method. We show compelling evidences that the energy spectrum of the boundary layer stays virtually the same over factor of 20 variations of the source luminosity. It is rather wide and cannot be described by a single temperature blackbody spectrum, probably because of the inhomogeneity of the boundary layer and a spread in the colour temperature. The observed maximum colour temperature of the boundary/spreading layer emission of kT~2.4-2.6 keV is very close to the maximum observed colour temperature in the photospheric radius expansion X-ray bursts, which is set by the limiting Eddington flux at the neutron star surface. Observed stability of the boundary layer spectrum and its maximum colour temperature strongly supports theoretical models of the boundary/spreading layers on surfaces of luminous accreting neutron stars, which assume the presence of a region emitting at the local Eddington limit. Variations in the luminosity in that case lead to changes in the size of this region, but affect less the spectral shape. Elaboration of this model will provide solid theoretical grounds for measurements of the neutron star sizes using the emission of the boundary/spreading layers of luminous accreting neutron stars.

A model with a condensed iron surface and partially ionized hydrogen-thin atmosphere allows us to fit simultaneously the observed general spectral shape and the broad absorption feature (observed at 0.3 keV) in different spin phases of the isolated neutron star RBS 1223. We constrain some physical properties of the X-ray emitting areas, i.e. the temperatures (Tpole1 ~ 105eV, Tpole2 ~99eV), magnetic field strengths Bpole1 ~ Bpole2 ~ 8.6x10^13G) at the poles, and their distribution parameters (a1 ~ 0.61, a2 ~ 0.29, indicating an absence of strong toroidal magnetic field component). In addition, we are able to place some constraints on the geometry of the emerging X-ray emission and the gravitational redshift (z ~0.16+0.03-0.01) of the isolated neutron star RBS 1223.

Results of photometric and spectroscopic investigations of the recently discovered eclipsing cataclysmic variable star HBHA 4705-03 are presented. The emission spectra of the system show broad hydrogen and helium emission lines. The bright spots with an approximately zero velocity components are found in the Doppler maps for the hydrogen and ionized helium lines. The disc structure is more prominent in the maps for the neutral helium lines. The masses of the components (M_WD = 0.54 \pm 0.10 M_sun and M_RD = 0.45 \pm 0.05 M_sun), and the orbit inclination (i = 71.8 \pm 0.^7 deg) were estimated using the radial velocity light curve and the eclipse width. The modeling of the light curve allows us to evaluate the bright spot parameters and the mass accretion rate (\dot M \approx 2 10^{17} g s^{-1}).

The Chandra / LETG spectrum of SS Cyg in outburst shows broad (\approx 5 A) spectral features that have been interpreted as a large number of absorption lines on a blackbody continuum with a temperature of 250 kK (Mauche 2004). It is most probable that this is the spectrum of the fast-rotating optically thick boundary layer on the white dwarf surface. Here we present the results of fitting this spectrum with high gravity hot stellar model atmospheres. An extended set of LTE model atmospheres with solar chemical composition was computed for this purpose. The best fit is obtained with the following parameters: T_eff=190 kK, log g=6.2, and N_H=8 10^{19} cm^{-2}. The spectrum of this model describes the observed spectrum in the 60--125 A range reasonably well, but at shorter wavelengths the observed spectrum has much higher flux. The reasons for this are discussed. The derived low surface gravity supports the hypothesis of the fast rotating boundary layer.

Theoretical spectra of X-ray bursting neutron star (NS) model atmospheres are widely used to determine the basic NS parameters such as their masses and radii. We construct accurate NS atmosphere models using for the first time an exact treatment of Compton scattering via the integral relativistic kinetic equation. We also compare the results with the previous calculations based on the Kompaneets operator. We solve the radiation transfer equation together with the hydrostatic equilibrium equation accounting exactly for the radiation pressure by electron scattering. We thus construct a new set of plane-parallel atmosphere models in LTE for hot NSs. The models were computed for six chemical compositions (pure H, pure He, solar H/He mix with various heavy elements abundances Z = 1, 0.3, 0.1, and 0.01 Z_sun, and three log g = 14.0, 14.3, and 14.6. For each chemical composition and log g, we compute more than 26 model atmospheres with various luminosities relative to the Eddington luminosity L_Edd computed for the Thomson cross-section. The maximum relative luminosities L/L_Edd reach values of up to 1.1 for high gravity models. The emergent spectra of all models are redshifted and fitted by diluted blackbody spectra in the 3--20 keV energy range appropriate for the RXTE/PCA. We also compute the color correction factors f_c. The radiative acceleration g_rad in our luminous, hot-atmosphere models is significantly smaller than in corresponding models based on the Kompaneets operator, because of the Klein-Nishina reduction of the electron scattering cross-section, and therefore formally "super-Eddington" model atmospheres do exist. The differences between the new and old model atmospheres are small for L / L_Edd < 0.8. For the same g_rad / g, the new f_c are slightly larger (by approximately 1%) than the old values.

We present the results of our photometric and spectroscopic studies of the new eclipsing cataclysmic variable star 1RXS J180834.7+101041. Its spectrum exhibits double-peaked hydrogen and helium emission lines. The Doppler maps constructed from hydrogen lines show a nonuniform distribution of emission in the disk similar to that observed in IP Peg. This suggests that the object can be a cataclysmic variable with tidal density waves in the disk. We have determined the component masses (M_WD =0.8 \pm 0.22 M_sun and M_RD =0.14 \pm 0.02 M_sun) and the binary inclination (i =78 \pm 1.5 deg) based on well-known relations between parameters for cataclysmic variable stars. We have modeled the binary light curves and showed that the model of a disk with two spots is capable of explaining the main observed features of the light curves.

All the neutron star (NS) atmosphere models published so far have been calculated in the "cold plasma approximation", which neglects the relativistic effects in the radiative processes, such as cyclotron emission/absorption at harmonics of cyclotron frequency. Here we present new NS atmosphere models which include such effects. We calculate a set of models for effective temperatures T_eff =1-3 MK and magnetic fields B \sim 10^{10}-10^{11} G, typical for the so-called central compact objects (CCOs) in supernova remnants, for which the electron cyclotron energy E_{c,e} and its first harmonics are in the observable soft X-ray range. Although the relativistic parameters, such as kT_eff /(m_e c^2) and E_{c,e} /(m_e c^2), are very small for CCOs, the relativistic effects substantially change the emergent spectra at the cyclotron resonances, E \approx sE_{c,e} (s=1, 2,...). Although the cyclotron absorption features can form in a cold plasma due to the quantum oscillations of the free-free opacity, the shape and depth of these features change substantially if the relativistic effects are included. In particular, the features acquire deep Doppler cores, in which the angular distribution of the emergent intensity is quite different from that in the cold plasma approximation. The relative contributions of the Doppler cores to the equivalent widths of the features grow with increasing the quantization parameter b_eff = E_{c,e}/kT_eff and harmonic number s. The total equivalent widths of the features can reach \sim 150-250 eV; they increase with growing b_eff and are smaller for higher harmonics.

Thermal emission during X-ray bursts is a powerful tool to determine neutron star masses and radii, if the Eddington flux and the apparent radius in the cooling tail can be measured accurately, and distances to the sources are known. We propose here an improved method of determining the basic stellar parameters using the data from the cooling phase of photospheric radius expansion bursts covering a large range of luminosities. Because at that phase the blackbody apparent radius depends only on the spectral hardening factor (color-correction), we suggest to fit the theoretical dependences of the color-correction versus flux in Eddington units to the observed variations of the inverse square root of the apparent blackbody radius with the flux. We show that spectral variations observed during a long photospheric radius expansion burst from 4U 1724-307 are entirely consistent with the theoretical expectations for the passively cooling neutron star atmospheres. Our method allows us to determine both the Eddington flux (which is found to be smaller than the touchdown flux by 15%) and the ratio of the stellar apparent radius to the distance much more reliably. We then find a lower limit on the neutron star radius of 14 km for masses below 2.2M_sun, independently of the chemical composition. These results suggest that the matter inside neutron stars is characterized by a stiff equation of state. We finally show that the apparent blackbody emitting area in the cooling tails of the short bursts from 4U 1724-307 is two times smaller than that for the long burst and their evolution does not follow the theory. This makes their usage for determination of the neutron star parameters questionable and casts serious doubts on the results of previous works that used for the analysis similar bursts from other sources. [abridged]

We have investigated the influence of X-ray irradiation on the vertical structure of the outer accretion disk in low-mass X-ray binaries by performing a self-consistent calculation of the vertical structure and X-ray radiation transfer in the disk. Penetrating deep into the disk, the field of scattered X-ray photons with energy $E\gtrsim10$\,keV exerts a significant influence on the vertical structure of the accretion disk at a distance $R\gtrsim10^{10}$\,cm from the neutron star. At a distance $R\sim10^{11}$\,cm, where the total surface density in the disk reaches $\Sigma_0\sim20$\,g\,cm$^{-2}$, X-ray heating affects all layers of an optically thick disk. The X-ray heating effect is enhanced significantly in the presence of an extended atmospheric layer with a temperature $T_{atm}\sim(2\div3)\times10^6$\,K above the accretion disk. We have derived simple analytic formulas for the disk heating by scattered X-ray photons using an approximate solution of the transfer equation by the Sobolev method. This approximation has a $\gtrsim10$\,% accuracy in the range of X-ray photon energies $E<20$\,keV.

The data from all observations of RBS 1223 (1RXS J130848.6+212708) conducted by XMM-Newton EPIC pn with the same instrumental setup in 2003-2007 were combined to form spin-phase resolved spectra. A number of complex models of neutron stars with strongly magnetized ($B_{pole} 10^{13}-10^{14}$ G) surface, with temperature and magnetic field distributions around magnetic poles, and partially ionized hydrogen thin atmosphere above it have been implemented into the X-ray spectral fitting package XSPEC for simultaneous fitting of phase-resolved spectra. A Markov-Chain-Monte-Carlo (MCMC) approach is also applied to verify results of fitting and estimating in multi parameter models. The spectra in different rotational phase intervals and light curves in different energy bands with high S/N ratio show a high complexity. The spectra can be parameterized with a Gaussian absorption line superimposed on a blackbody spectrum, while the light curves with double-humped shape show strong dependence of pulsed fraction upon the energy band (13%-42%), which indicates that radiation emerges from at least two emitting areas. A model with condensed iron surface and partially ionized hydrogen thin atmosphere above it allows us to fit simultaneously the observed general spectral shape and the broad absorption feature observed at 0.3 keV in different spin phases of RBS 1223. It allowed to constrain some physical properties of X-ray emitting areas, i.e. the temperatures ($T_{p1} ~ 105$ eV, $T_{p2} ~ 99$ eV), magnetic field strengths ($B_{p1}= B_{p2} ~ 8.6x10^{13}$G) at the poles, and their distributions parameters ($a_{1} 0.61, a_{2} 0.29$, indicating an absence of strong toroidal magnetic field component). In addition, it puts some constraints on the geometry of the emerging X-ray emission and gravitational redshift ($z=0.16^{+0.03}_{-0.01}$) of RBS 1223.

The eruption of the recurrent nova U Scorpii on 28 January 2010 is now the all-time best observed nova event. We report 36,776 magnitudes throughout its 67 day eruption, for an average of one measure every 2.6 minutes. This unique and unprecedented coverage is the first time that a nova has any substantial amount of fast photometry. With this, two new phenomena have been discovered: the fast flares in the early light curve seen from days 9-15 (which have no proposed explanation) and the optical dips seen out of eclipse from days 41-61 (likely caused by raised rims of the accretion disk occulting the bright inner regions of the disk as seen over specific orbital phases). The expanding shell and wind cleared enough from days 12-15 so that the inner binary system became visible, resulting in the sudden onset of eclipses and the turn-on of the supersoft X-ray source. On day 15, a strong asymmetry in the out-of-eclipse light points to the existence of the accretion stream. The normal optical flickering restarts on day 24.5. For days 15-26, eclipse mapping shows that the optical source is spherically symmetric with a radius of 4.1 R_sun. For days 26-41, the optical light is coming from a rim-bright disk of radius 3.4 R_sun. For days 41-67, the optical source is a center-bright disk of radius 2.2 R_sun. Throughout the eruption, the colors remain essentially constant. We present 12 eclipse times during eruption plus five just after the eruption.

We present an analysis of the Vela X-1's "off-states" based on Suzaku observations taken in June 2008. Defined as states in which the flux sudden decreases below the instrumental sensitivity, these "off-states" have been interpreted by several authors as the onset of the "propeller regime". For the first time ever, however, we find that the source does not turn off and, although the flux drops by a factor of 20 during the three recorded "off-states", pulsations are still observed. The spectrum and the pulse-profiles of the "off-states" are also presented. Eventually, we discuss our findings in framework of the "gated accretion" scenario and conclude that most likely the residual flux is due to the accretion of matter leaking through the magnetosphere by means of Kelvin-Helmholz instabilities (KHI).

We present the results of our optical identification of the X-ray source IGR J16547-1916 detected by the INTEGRAL observatory during a deep all-sky survey. Analysis of the spectroscopic data from the SWIFT and INTEGRAL observatories in the X-ray energy band and from the BTA (Special Astrophysical Observatory) telescope in the optical band has shown that the source is most likely an intermediate polar -- an accreting white dwarf with the mass of M~0.85 M_Sun in a low-mass binary system. Subsequent studies of the object's rapid variability with the RTT-150 telescope have confirmed this conclusion by revealing periodic pulsations of its optical emission with a period of ~550 s.

We present an extended set of model atmospheres and emergent spectra of X-ray bursting neutron stars in low mass X-ray binaries. Compton scattering is taken into account. The models were computed in LTE approximation for six different chemical compositions: pure hydrogen and pure helium atmospheres, and atmospheres with a solar mix of hydrogen and helium and various heavy elements abundances: Z = 1, 0.3, 0.1, and 0.01 Z_sun, for three values of gravity, log g =14.0, 14.3, and 14.6 and for 20 values of relative luminosity l = L/L_Edd in the range 0.001 - 0.98. The emergent spectra of all models are fitted by diluted blackbody spectra in the observed RXTE/PCA band 3 - 20 keV and the corresponding values of color correction factors f_c are presented. We also show how to use these dependencies to estimate the neutron star's basic parameters.

The X-ray spectra of some magnetized isolated neutron stars (NSs) show absorption features with equivalent widths (EWs) of 50 - 200 eV, whose nature is not yet well known. To explain the prominent absorption features in the soft X-ray spectra of the highly magnetized (B ~ 10^{14} G) X-ray dim isolated NSs (XDINSs), we theoretically investigate different NS local surface models, including naked condensed iron surfaces and partially ionized hydrogen model atmospheres, with semi-infinite and thin atmospheres above the condensed surface. We also developed a code for computing light curves and integral emergent spectra of magnetized neutron stars with various temperature and magnetic field distributions over the NS surface. We compare the general properties of the computed and observed light curves and integral spectra for XDINS RBS\,1223 and conclude that the observations can be explained by a thin hydrogen atmosphere above the condensed iron surface, while the presence of a strong toroidal magnetic field component on the XDINS surface is unlikely. We suggest that the harmonically spaced absorption features in the soft X-ray spectrum of the central compact object (CCO) 1E 1207.4-5209 (hereafter 1E 1207) correspond to peaks in the energy dependence of the free-free opacity in a quantizing magnetic field, known as quantum oscillations. To explore observable properties of these quantum oscillations, we calculate models of hydrogen NS atmospheres with B ~ 10^{10} - 10^{11} G (i.e., electron cyclotron energy E_{c,e} ~ 0.1 - 1 keV) and T_eff = 1 - 3 MK. Such conditions are thought to be typical for 1E 1207. We show that observable features at the electron cyclotron harmonics with EWs \approx 100 - 200 eV can arise due to these quantum oscillations.

Results of photometric and spectroscopic investigations of the recently discovered disc cataclysmic variable star 1RXS J180834.7+101041 are presented. Emission spectra of the system show broad double peaked hydrogen and helium emission lines. Doppler maps for the hydrogen lines demonstrate strongly non-uniform emissivity distribution in the disc, similar to that found in IP Peg. It means that the system is a new cataclysmic variable with a spiral density wave in the disc. Masses of the components (M_WD = 0.8 +/- 0.22 M_sun and M_RD = 0.14 +/- 0.02 M_sun), and the orbit inclination (i = 78 +/- 1.5 deg) were estimated using the various well-known relations for cataclysmic variables.

The spectral energy distribution (SED) of the TT Ari system, which is well known from published IUE and optical photometric observations, was modeled by a steady-state accretion \alpha-disc around a white dwarf. Parameters of the system were derived from time-resolved optical spectral observations in the bright state that we obtained in Sep. 1998. The radial velocity semi-amplitude of the white dwarf (33.8 +/- 2.5 km/s) and corresponding mass function (f(M) = 5.5 +/- 1.2 *10^{-4}~ M_sun) were derived from the motion of the emission components of Balmer lines. The mass ratio q (\approx 0.315) was evaluated from the fractional period excess of the superhump period over the orbital period \epsilon (\approx 0.085), and a secondary mass range (0.18 - 0.38 M_sun) was estimated from the orbital period. Therefore, the white dwarf mass range is 0.57 - 1.2 M_sun and the inclination angle of the system to the line of sight is 17 - 22.5 degrees. The adopted distance to the system is 335 +/- 50 pc. To fit the observed SED it is necessary to add a thermal spectrum with T \approx 11600 K and luminosity \approx 0.4 L_disk to the accretion disc spectrum. This combined spectrum successfully describes the observed Balmer lines absorption components. Formally the best fit of the HeI 4471 line gives minimum masses of the components (M_RD = 0.18 M_sun and M_WD = 0.57 M_sun), with the corresponding inclination angle i = 22.1 deg and mass-accretion rate \dot M = 2.6 * 10^{17} g/s.

X-ray bursting neutron stars in low mass X-ray binaries constitute an appropriate source class to constrain masses and radii of neutron stars, but a sufficiently extended set of corresponding model atmospheres is necessary for these investigations. We computed such a set of model atmospheres and emergent spectra in a plane-parallel, hydrostatic, and LTE approximation with Compton scattering taken into account. The models were calculated for six different chemical compositions: pure hydrogen and pure helium atmospheres, and atmospheres with solar mix of hydrogen and helium, and various heavy element abundances Z = 1, 0.3, 0.1, and 0.01 Z_sun. For each chemical composition the models are computed for three values of surface gravity, log g =14.0, 14.3, and 14.6, and for 20 values of the luminosity in units of the Eddington luminosity, L/L_Edd, in the range 0.001--0.98. The emergent spectra of all models are redshifted and fitted by a diluted blackbody in the RXTE/PCA 3--20 keV energy band, and corresponding values of the color correction (hardness factors) f_c are presented. Theoretical dependences f_c - L/L_Edd can fitted to the observed dependence K^{-1/4} - F of the blackbody normalization K on flux during cooling stages of X-ray bursts to determine the Eddington flux and the ratio of the apparent neutron star radius to the source distance. If the distance is known, these parameters can be transformed to the constraints on neutron star mass and radius. The theoretical atmosphere spectra can also be used for direct comparison with the observed X-ray burst spectra.

Extending the population synthesis method to isolated young cooling white dwarfs we are able to confront our model assumptions with observations made in ROSAT All-Sky Survey (Fleming et al., 1996). This allows us to check model parameters such as evolution of spectra and separation of heavy elements in DA WD envelopes. It seems like X-ray spectrum temperature of these objects is given by the formula T_{X-ray} = min(T_eff, T_max). We have obtained DA WD's birth rate and upper limit of the X-ray spectrum temperature: DA birth rate $= 0.61\times 10^{-12}$ in cubic parsec per year and T_max = 41000 K. These values are in good correspondence with values obtained by other authors (Liebert et al., 2004; Wolff et al., 1996). From this fact we also conclude that our population synthesis method is applicable to the population of close-by isolated cooling white dwarfs as well as to the population of the isolated cooling neutron stars.

We present an analysis of time-resolved, medium resolution optical spectroscopic observations of UX UMa in the blue (3920-5250 A) and red (6100-7200 A) wavelength ranges, that were obtained in April 1999 and March 2008 respectively. The observed characteristics of our spectra indicate that UX UMa has been in different states during those observations. The blue spectra are very complex. They are dominated by strong and broad single-peaked emission lines of hydrogen. The high-excitation lines of HeII 4686 and the Bowen blend are quite strong as well. All the lines consist of a mixture of absorption and emission components. Using Doppler tomography we have identified four distinct components of the system: the accretion disc, the secondary star, the bright spot from the gas stream/disc impact region, and the unique compact area of absorption in the accretion disc seen as a dark spot in the lower-left quadrant of the tomograms. In the red wavelength range, both the hydrogen (H_alpha) and neutral helium (HeI 6678 and HeI 7065) lines were observed in emission and both exhibited double-peaked profiles. Doppler tomography of these lines reveals spiral structure in the accretion disc, but in contrast to the blue wavelength range, there is no evidence for either the dark spot or the gas stream/disc impact region emission, while the emission from the secondary star is weak. During the observations in 1999, UX UMa showed many of the defining properties of the SW Sex stars. However, all these features almost completely disappeared in 2008. We have also estimated the radial velocity semi-amplitudes K_1 and K_2 and evaluated the system parameters of UX UMa. These estimates are inconsistent with previous values derived by means of analysis of WD eclipse features in the light curve in the different wavelength ranges.

In the X-ray spectra of most X-ray dim isolated neutron stars (XDINSs) absorption features with equivalent widths (EWs) of 50 -- 200 eV are observed. We theoretically investigate different models to explain absorption features and compare their properties with the observations. We consider various theoretical models for the magnetized neutron star surface: naked condensed iron surfaces and partially ionized hydrogen model atmospheres, including semi-infinite and thin atmospheres above a condensed surface. The properties of the absorption features (especially equivalent widths) and the angular distributions of the emergent radiation are described for all models. A code for computing light curves and integral emergent spectra of magnetized neutron stars is developed. We assume a dipole surface magnetic field distribution with a possible toroidal component and corresponding temperature distribution. A model with two uniform hot spots at the magnetic poles can also be employed. Light curves and spectra of highly magnetized neutron stars with parameters typical for XDINSs are computed using different surface temperature distributions and various local surface models. Spectra of magnetized model atmospheres are approximated by diluted blackbody spectra with one or two Gaussian lines having parameters, which allow us to describe the model absorption features. To explain the prominent absorption features in the soft X-ray spectra of XDINSs a thin atmosphere above the condensed surface can be invoked, whereas a strong toroidal magnetic field component on the XDINS surfaces can be excluded.

Some isolated neutron stars show harmonically spaced absorption features in their thermal soft X-ray spectra. The interpretation of the features as a cyclotron line and its harmonics has been suggested, but the usual explanation of the harmonics as caused by relativistic effects fails because the relativistic corrections are extremely small in this case. We suggest that the features correspond to the peaks in the energy dependence of the free-free opacity in a quantizing magnetic field, known as quantum oscillations. The peaks arise when the transitions to new Landau levels become allowed with increasing the photon energy; they are strongly enhanced by the square-root singularities in the phase-space density of quantum states in the case when the free (non-quantized) motion is effectively one-dimensional. To explore observable properties of these quantum oscillations, we calculate models of hydrogen neutron star atmospheres with B \sim 10^{10} - 10^{11} G (i.e., electron cyclotron energy E_{c,e} = 0.1 - 1 keV) and T_{eff} = 1 - 3 MK. Such conditions are thought to be typical for the so-called central compact objects in supernova remnants, such as 1E 1207.4-5209 in PKS 1209-51/52. We show that observable features at the electron cyclotron harmonics form at moderately large values of the quantization parameter, b_{eff} = E_{c,e}/kT_{eff} = 0.5 - 20. The equivalent widths of the features can reach 100 - 200 eV; they grow with increasing b_{eff} and are lower for higher harmonics.

We present the results of an in-depth study of the long-period X-ray pulsar GX 301-2. Using archival data of INTEGRAL, RXTE ASM, and CGRO BATSE, we study the spectral and timing properties of the source. Comparison of our timing results with previously published work reveals a secular decay of the orbital period at a rate of \simeq -3.25 \times 10^{-5} d yr^{-1}, which is an order of magnitude faster than for other known systems. We argue that this is probably result either of the apsidal motion or of gravitational coupling of the matter lost by the optical companion with the neutron star, although current observations do not allow us to distinguish between those possibilities. We also propose a model to explain the observed long pulse period. We find that a very strong magnetic field B \sim 10^{14} G can explain the observed pulse period in the framework of existing models for torques affecting the neutron star. We show that the apparent contradiction with the magnetic field strength B_{CRSF} \sim 4 \times 10^{12} G derived from the observed cyclotron line position may be resolved if the line formation region resides in a tall accretion column of height \sim 2.5 - 3 R_{NS}. The color temperature measured from the spectrum suggests that such a column may indeed be present, and our estimates show that its height is sufficient to explain the observed cyclotron line position.

We present a detailed investigation of SBS1150+599A, a close binary star hosted by the planetary nebula PN G135.9+55.9 (TS01, Stasinska et al, 2009). The nebula, located in the Galactic halo, is the most oxygen-poor one known to date and is the only one known to harbor a double degenerate core. We present XMM-Newton observations of this object, which allowed the detection of the previously invisible component of the binary core, whose existence was inferred so far only from radial velocity and photometric variations. The parameters of the binary system were deduced from a wealth of information via three independent routes using the spectral energy distribution (from the infrared to X-rays), the light and radial velocity curves, and a detailed model atmosphere fitting of the stellar absorption features of the optical/UV component. We find that the cool component must have a mass of 0.54+/-0.2 Msun, an average effective temperature, Teff, of 58000+/-3000 K, a mean radius of 0.43+/-0.3 Rsun, a gravity log g=5.0+/-0.3, and that it nearly fills its Roche lobe. Its surface elemental abundances are found to be: 12 + log He/H = 10.95+/-0.04 dex, 12 + log C/H = 7.20+/-0.3 dex, 12 + log N/H < 6.92 and 12 + log O/H < 6.80, in overall agreement with the chemical composition of the planetary nebula. The hot component has Teff = 160-180 kK, a luminosity of about ~10e4 Lsun and a radius slightly larger than that of a white dwarf. It is probably bloated and heated as a result of intense accretion and nuclear burning on its surface in the past. The total mass of the binary system is very close to Chandrasekhar limit. This makes TS01 one of the best type Ia supernova progenitor candidates. We propose two possible scenarios for the evolution of the system up to its present stage.

We present a new computer code for modeling magnetized neutron star atmospheres in a wide range of magnetic fields (10^{12} - 10^{15} G) and effective temperatures (3 \times 10^5 - 10^7 K). The atmosphere is assumed to consist either of fully ionized electron-ion plasmas or of partially ionized hydrogen. Vacuum resonance and partial mode conversion are taken into account. Any inclination of the magnetic field relative to the stellar surface is allowed. We use modern opacities of fully or partially ionized plasmas in strong magnetic fields and solve the coupled radiative transfer equations for the normal electromagnetic modes in the plasma. Using this code, we study the possibilities to explain the soft X-ray spectra of isolated neutron stars by different atmosphere models. In particular, the outgoing spectrum using the "sandwich" model (thin atmosphere with a hydrogen layer above a helium layer) is constructed. Thin partially ionized hydrogen atmospheres with vacuum polarization are shown to be able to improve our understanding of the observed spectrum of the nearby isolated neutron star RBS 1223 (RX J1308.8+2127).

Observed X-ray spectra of some isolated magnetized neutron stars display absorption features, sometimes interpreted as ion cyclotron lines. Modeling the observed spectra is necessary to check this hypothesis and to evaluate neutron star parameters.We develop a computer code for modeling magnetized neutron star atmospheres in a wide range of magnetic fields (10^{12} - 10^{15} G) and effective temperatures (3 \times 10^5 - 10^7 K). Using this code, we study the possibilities to explain the soft X-ray spectra of isolated neutron stars by different atmosphere models. The atmosphere is assumed to consist either of fully ionized electron-ion plasmas or of partially ionized hydrogen. Vacuum resonance and partial mode conversion are taken into account. Any inclination of the magnetic field relative to the stellar surface is allowed. We use modern opacities of fully or partially ionized plasmas in strong magnetic fields and solve the coupled radiative transfer equations for the normal electromagnetic modes in the plasma. Spectra of outgoing radiation are calculated for various atmosphere models: fully ionized semi-infinite atmosphere, thin atmosphere, partially ionized hydrogen atmosphere, or novel "sandwich" atmosphere (thin atmosphere with a hydrogen layer above a helium layer. Possibilities of applications of these results are discussed. In particular, the outgoing spectrum using the "sandwich" model is constructed. Thin partially ionized hydrogen atmospheres with vacuum polarization are shown to be able to improve the fit to the observed spectrum of the nearby isolated neutron star RBS 1223 (RX J1308.8+2127).

We address the task of modeling soft X-ray and optical light curves of X-ray novae in the high/soft state. The analytic model of viscous evolution of an externally truncated accretion \alpha-disk is used. Relativistic effects near a Kerr black hole and self-irradiation of an accretion disk are taken into account. The model is applied to the outbursts of X-ray nova Monocerotis 1975 (A 0620-00) and X-ray nova Muscae 1991 (GRS 1124-68). Comparison of observational data with the model yields constraints on the angular momentum (the Kerr parameter) of the black holes in A 0620-00 and GRS 1124-68: 0.3-0.6 and \leq 0.4, and on the viscosity parameter \alpha of the disks: 0.7-0.95 and 0.55-0.75. We also conclude that the accretion disks should have an effective geometrical thickness 1.5-2 times greater than the theoretical value of the distance between the photometric layers.

The discovery of photospheric absorption lines in XMM-Newton spectra of the X-ray bursting neutron star in EXO0748-676 by Cottam and collaborators allows us to constrain the neutron star mass-radius ratio from the measured gravitational redshift. A radius of R=9-12km for a plausible mass range of M=1.4-1.8Msun was derived by these authors. It has been claimed that the absorption features stem from gravitationally redshifted (z=0.35) n=2-3 lines of H- and He-like iron. We investigate this identification and search for alternatives. We compute LTE and non-LTE neutron-star model atmospheres and detailed synthetic spectra for a wide range of effective temperatures (effective temperatures of 1 - 20MK) and different chemical compositions. We are unable to confirm the identification of the absorption features in the X-ray spectrum of EXO0748-676 as n=2-3 lines of H- and He-like iron (Fe XXVI and Fe XXV). These are subordinate lines that are predicted by our models to be too weak at any effective temperature. It is more likely that the strongest feature is from the n=2-3 resonance transition in Fe XXIV with a redshift of z=0.24. Adopting this value yields a larger neutron star radius, namely R=12-15km for the mass range M=1.4-1.8Msun, favoring a stiff equation-of-state and excluding mass-radius relations based on exotic matter. Combined with an estimate of the stellar radius R>12.5km from the work of Oezel and collaborators, the z=0.24 value provides a minimum neutron-star mass of M>1.48Msun, instead of M>1.9Msun, when assuming z=0.35.

The majority of cataclysmic variables observed in the hard X-ray energy band are intermediate polars where the magnetic field is strong enough to channel the accreting matter to the magnetic poles of the white dwarf. A shock above the stellar surface heats the gas to fairly high temperatures (10--100 keV). The post-shock region cools mostly via optically thin bremsstrahlung. We investigate the influence of Compton scattering on the structure and the emergent spectrum of the post-shock region. We also study the effect it has on the mass of the white dwarfs obtained from fitting the observed X-ray spectrum of intermediate polars. We construct the model of the post-shock region taking Compton scattering into account. The radiation transfer equation is solved in the plane-parallel approximation. The feedback of Compton scattering on the structure of the post-shock region is also accounted for. A set of the post-shock region model spectra for various white dwarf masses is calculated. We find that Compton scattering does not change the emergent spectra significantly for low accretion rates or low white dwarf masses. However, it becomes important at high accretion rates and high white dwarf masses. The time-averaged, broad-band X-ray spectrum of intermediate polar V709 Cas obtained by the RXTE and INTEGRAL observatories is fitted using the set of computed spectral models. We obtained the white dwarf mass of 0.91 +/- 0.02 M_sun and 0.88 +/- 0.02 M_sun using models with Compton scattering taken into account and without it, respectively.

Emergent model spectra of neutron star atmospheres are widely used to fit the observed soft X-ray spectra of different types of isolated neutron stars. We investigate the effect of Compton scattering on the emergent spectra of hot (T_eff > 10^6 K) isolated neutron stars with weak magnetic fields. In order to compute model atmospheres in hydrostatic and radiative equilibrium we solve the radiation transfer equation with the Kompaneets operator. We calculate a set of models with effective temperatures in the range 1 - 5 * 10^6 K, with two values of surface gravity (log g = 13.9 and 14.3) and different chemical compositions. Radiation spectra computed with Compton scattering are softer than those computed without Compton scattering at high energies (E > 5 keV) for light elements (H or He) model atmospheres. The Compton effect is more significant in H model atmospheres and models with low surface gravity. The emergent spectra of the hottest (T_eff > 3 * 10^6 K) model atmospheres can be described by diluted blackbody spectra with hardness factors ~ 1.6 - 1.9. Compton scattering is less important in models with solar abundance of heavy elements.

Model atmospheres of isolated neutron stars with low magnetic field are calculated with Compton scattering taking into account. Models with effective temperatures 1, 3 and 5 MK, with two values of surface gravity log(g)g = 13.9 and 14.3), and different chemical compositions are calculated. Radiation spectra computed with Compton scattering are softer than the computed with Thomson scattering at high energies (E > 5 keV) for hot (T_eff > 1 MK) atmospheres with hydrogen-helium composition. Compton scattering is more significant to hydrogen models with low surface gravity. The emergent spectra of the hottest (T_eff > 3 MK) model atmospheres can be described by diluted blackbody spectra with hardness factors ~ 1.6 - 1.9. Compton scattering is less important for models with solar abundance of heavy elements.

We describe our recent progress in modeling supernova-fallback disks and neutron star (NS) atmospheres. We present a first detailed spectrum synthesis calculation of a SN-fallback disk composed of iron. We assume a thin disk with a radial structure described by the alpha-disk model. The vertical structure and emission spectrum are computed self-consistently by solving the structure equations simultaneously with the radiation transfer equations under non-LTE conditions. We describe the properties of a specific disk model and discuss various effects on the emergent UV/optical spectrum. We investigate Compton scattering effects on the thermal spectrum of NSs. In addition, we constructed a new generation of metal line-blanketed non-LTE model atmospheres for NSs. It is compared to X-ray burst spectra of EXO0748-676. It is possible that the gravitational redshift, deduced from absorption lines, is lower (z=0.24) than hitherto assumed (z=0.35). Accordingly, this would result in NS mass and radius lower limits of M>1.63Msun and R>13.8 km.

We present the results of spectroscopic-- and orbit--sampled photometric observations of the faint UV-excess object PG 2200+085. The optical CCD photometry observations of this object were performed by the Russian-Turkish 1.5-meter telescope RTT150 at the TUBITAK National Observatory (Turkey). The long-slit optical spectroscopy observations with 2.6 A resolution were carried out by 6-meter telescope BTA at the Special Astrophysical Observatory (Russia). The photometric variations over two nights are almost sinusoidal with an amplitude \Delta m_V = 0.04 m and a period of P = 0.3186 d. Such a light curve is typical of a detached close binary with an illumination effect or the ellipsoidal deformation of a secondary star. The observed spectrum clearly displays a featureless blue continuum of a hot component and a rich absorption--line and molecular band K--star spectrum. The CaII line profiles with strong emission cores are remarkably similar to those of V471 Tau. We tentatively classify PG 2200+085 as a pre-cataclysmic binary of the V471 Tau type.

Spectra of the spreading layers on the neutron star surface are calculated on the basis of the Inogamov-Sunyaev model taking into account general relativity correction to the surface gravity and considering various chemical composition of the accreting matter. Local (at a given latitude) spectra are similar to the X-ray burst spectra and are described by a diluted black body. Total spreading layer spectra are integrated accounting for the light bending, gravitational redshift, and the relativistic Doppler effect and aberration. They depend slightly on the inclination angle and on the luminosity. These spectra also can be fitted by a diluted black body with the color temperature depending mainly on a neutron star compactness. Owing to the fact that the flux from the spreading layer is close to the critical Eddington, we can put constraints on a neutron star radius without the need to know precisely the emitting region area or the distance to the source. The boundary layer spectra observed in the luminous low-mass X-ray binaries, and described by a black body of color temperature Tc=2.4+-0.1 keV, restrict the neutron star radii to R=14.8+- 1.5 km (for a 1.4-Msun star and solar composition of the accreting matter), which corresponds to the hard equation of state.

We re-examine the effects of Compton scattering on the emergent spectra of hot DA white dwarfs in the soft X-ray range. Earlier studies have implied that sensitive X-ray observations at wavelengths $\lambda < 50$ \AA might be capable of probing the flux deficits predicted by the redistribution of electron-scattered X-ray photons toward longer wavelengths. We adopt two independent numerical approaches to the inclusion of Compton scattering in the computation of pure hydrogen atmospheres in hydrostatic equilibrium. One employs the Kompaneets diffusion approximation formalism, while the other uses the cross-sections and redistribution functions of Guilbert. Models and emergent spectra are computed for stellar parameters representative of HZ 43 and Sirius B, and for models with an effective temperature $T_{\rm eff} = 100 000$ K. The differences between emergent spectra computed for Compton and Thomson scattering cases are completely negligible in the case of both HZ 43 and Sirius B models, and are also negligible for all practical purposes for models with temperatures as high as $T_{\rm eff} = 100 000$ K. Models of the soft X-ray flux from these stars are instead dominated by uncertainties in their fundamental parameters.

ROSAT spectra of 11 supersoft X-ray sources (RX J0439.8-6809, RX J0513.9-6951, RX J0527.8-6954, CAL 87, CAL 83, 1E 0035.4-7230, RX J0048.4-7332, 1E 0056.8-7154, RX J0019.8 +2156, RX J0925.7-4758, AG Draconis) are approximated with theoretical spectra obtained in LTE models for the atmospheres of hot white dwarfs with line blanketing. The confidence intervals of parameters derived from these approximations $ T_{\rm eff}$, $\log g$, $N_H$, and $R^2/d^2$ are determined. The results are compared with predictions for a model with stable/recurrent thermonuclear burning on the white-dwarf surface.

The parameters of ten supersoft X-ray sources (RX J0439.8-6809, RX J0513.9-6951, RX J0527.8-6954, CAL 87, CAL 83, 1E 0035.4-7230, RX J0048.4-7332, 1E 0056.8-7154, RX J0019.8 +2156, RX J0925.7-4758) observed by ROSAT obtained using blanketing LTE model atmospheres are analyzed. The consistency of the resulting parameters with a model with stable/recurrent burning on the surface of the white dwarf is studied. The luminosity and sizes of seven of the sources are in good agreement with this model. The masses of the white dwarfs in these sources are estimated. A formula that can be used to estimate the masses of white dwarfs in classical supersoft sources based on their effective temperatures is presented.

We present the hard X-ray time-averaged spectrum of the intermediate polar V709 Cas observed with INTEGRAL. We performed the observation using data from the IBIS/ISGRI instrument in the 20-100 keV energy band and from JEM-X at lower energy (5-20 keV). Using different multi-temperature and density X-ray post-shock models we measured a improved post-shock temperature of ~40 keV and estimated the V709 Cas mass to be 0.82(+0.12/-0.25} Msun. We compare the resulting spectral parameters with previously reported BeppoSAX and RXTE observations.

We present results of an analysis of broadband spectra of 14 intermediate polars obtained with the RXTE observatory (PCA and HEXTE spectrometers, 3-100 keV). By means of our calculations of the structure and the emergent spectrum of the post-shock region of intermediate polars we fitted the observed spectra and obtained estimates of the white dwarf masses. The estimated masses are compared with masses obtained by other authors.

We present the broadband phase averaged spectrum of one of the brightest intermediate polars V1223 Sgr, obtained with INTEGRAL and RXTE observatories (3-100 keV). Good statistical quality of the spectrum in a hard X-ray energy band (INTEGRAL/IBIS and RXTE/HEXTE) allowed us to disentangle contributions of a direct optically thin plasma emission and a reflected component to the spectrum of V1223 Sgr. The obtained measurement of the post-shock temperature of the accreting matter give us the information about the mass of the white dwarf and the inclination of the system.

Detailed analysis of the lightcurve of CAL 87 clearly has shown that the high optical luminosity comes from the accretion disc rim and can only be explained by a severe thickening of the disc rim near the location where the accretion stream impinges. This area is irradiated by the X-rays where it faces the white dwarf. Only if the reprocessing rate of X-rays to optical light is high a luminosity as high as observed can be understood. But a recent detailed study of the soft X-ray radiation reprocessing in supersoft X-ray sources has shown that the efficiency is not high enough. We here propose a solution for this problem. As already discussed in the earlier lightcurve analysis the impact of the accretion stream at the outer disc rim produces a ``spray'', consisting of a large number of individual gas blobs imbedded in a surrounding corona. For the high mass flow rate this constitutes an optically thick vertically extended screen at the rim of the accretion disc. We analyse the optical properties of this irradiated spray and find that the multiple scattering between these gas blobs leads to an effective reprocessing of soft X-rays to optical light as required by the observations.

A Chandra ACIS-S imaging observation of the nearby galaxy M81 (NGC 3031) reveals 9 luminous soft X-ray sources. The local environments, X-ray spectral properties, and X-ray light curves of the sources are presented and discussed in the context of prevailing physical models for supersoft sources. It is shown that the sample falls within expectations based on population synthesis models taken from the literature though the high observed luminosities (L~2.0e36 to \~3.0e38 erg/s in the 0.2-2.0 keV band) and equivalent blackbody temperatures T~40 to 80 eV) place the brightest detected M81 objects at the high luminosity end of the class of supersoft sources defined by previous ROSAT and Einstein studies of nearby galaxies. This is interpreted as a natural consequence of the higher sensitivity of Chandra to hotter and more luminous systems. Most of the sources can be explained as canonical \ssss, accreting white dwarfs powered by steady surface nuclear burning, with X-ray spectra well-fit by hot white dwarf local thermodynamic equilibrium atmosphere models. An exceptionally bright source is scrutinized in greater detail as its estimated bolometric luminosity, L~1.5e39 erg/s, greatly exceeds theoretical estimates for supersoft sources. This source may be beyond the stability limit and undergoing a phase of mass outflow under extreme conditions. Alternatively, a model in which the observed X-ray spectrum arises from an accretion disk around a blackhole of mass \~1200/sqrt(cos(i)) solar masses (viewed at an inclination angle i) cannot be excluded.