Model atmospheres of X-ray bursting neutron stars

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.Comment: 2 pages, 1 figure, conference "Astrophysics of Neutron Stars - 2010" in honor of M. Ali Alpar, Izmir, Turke

Modeling the EUV spectra of optically thick boundary layers of dwarf novae in outburst

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.Comment: Accepted for publication in A&A, 14 pages, 13 figures, 4 table

On the compactness of the isolated neutron star RX J0720.4-3125

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$.Comment: Comments: 19 pages, 15 figures and 5 tables, Astronomy and Astrophysics accepted. arXiv admin note: text overlap with arXiv:1108.389

Absorption Features in Spectra of Magnetized Neutron Stars

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.Comment: 4 pages, 3 figures, conference "Astrophysics of Neutron Stars - 2010" in honor of M. Ali Alpar, Izmir, Turke

Vertical Structure of the Outer Accretion Disk in Persistent Low-Mass X-Ray Binaries

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.Comment: 19 pages, 8 figures, published in Astronomy Letter

The boundary layer in compact binaries

Context. 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. Although 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. Aims. It is the aim of this work to obtain more insight into the BL 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. Methods. We perform numerical hydrodynamical simulations, where the problem is treated in a one-dimensional, radial approximation (slim disk). The turbulence is described by the α parameter viscosity. We include both cooling from the disk surfaces and radial radiation transport. The radiation energy is treated in a one-temperature approximation. Results. For a given M our results show a strong dependence on the stellar mass and rotation rate. The midplane and the effective temperature rise considerably with increasing stellar mass or decreasing stellar rotation rate. Our simulations also 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. Conclusions. The combination of a high mass and a small radius, characteristic of 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*), 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. © ESO, 2013