Source author record

V. F. Suleimanov

V. F. Suleimanov appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

21works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

21 published item(s)

preprint2022arXiv

Likely optical counterpart of the cool middle-aged pulsar J1957+5033

The 840 kyr old pulsar PSR J1957+5033, detected so far only in $γ$- and X-rays, is a nearby and rather cool neutron star with a temperature of 0.2--0.3 MK, a distance of $\la$1 kpc, and a small colour reddening excess $E(B-V) \approx 0.03$. These properties make it an ideal candidate to detect in the optical to get additional constraints on its parameters. We thus performed the first deep optical observations of the pulsar with the 10.4-meter Gran Telescopio Canarias in the $g'$ band and found its possible counterpart with $g'=27.63\pm 0.26$. The counterpart candidate position is consistent with the X-ray coordinates of the pulsar within the 0.5 arcsec accuracy. Assuming that this is the real counterpart, we analysed the pulsar X-ray spectrum together with the derived optical flux density. As a result, we found that the thermal emission from the bulk surface of the cooling neutron star can significantly contribute to its optical flux. Our multi-wavelength spectral analysis favours the pulsar nature of the detected optical source, since it provides physically adequate parameters of the pulsar emission. We show that the optical data can provide new constraints on the pulsar temperature and distance.

preprint2021arXiv

Middle aged $γ$-ray pulsar J1957+5033 in X-rays: pulsations, thermal emission and nebula

We analyze new XMM-Newton and archival Chandra observations of the middle-aged $γ$-ray radio-quiet pulsar J1957+5033. We detect, for the first time, X-ray pulsations with the pulsar spin period of the point-like source coinciding by position with the pulsar. This confirms the pulsar nature of the source. In the 0.15--0.5 keV band, there is a single pulse per period and the pulsed fraction is $\approx18\pm6$ per cent. In this band, the pulsar spectrum is dominated by a thermal emission component that likely comes from the entire surface of the neutron star, while at higher energies ($\gtrsim0.7$ keV) it is described by a power law with the photon index $Γ\approx 1.6$. We construct new hydrogen atmosphere models for neutron stars with dipole magnetic fields and non-uniform surface temperature distributions with relatively low effective temperatures. We use them in the spectral analysis and derive the pulsar average effective temperature of $\approx(2-3)\times10^5$ K. This makes J1957+5033 the coldest among all known thermally emitting neutron stars with ages below 1 Myr. Using the interstellar extinction--distance relation, we constrain the distance to the pulsar in the range of 0.1--1 kpc. We compare the obtained X-ray thermal luminosity with those for other neutron stars and various neutron star cooling models and set some constraints on latter. We observe a faint trail-like feature, elongated $\sim 8$ arcmin from J1957+5033. Its spectrum can be described by a power law with a photon index $Γ=1.9\pm0.5$ suggesting that it is likely a pulsar wind nebula powered by J1957+5033.

preprint2020arXiv

Burst-induced coronal cooling in GS 1826-24 The clock wagging its tail

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'.

preprint2020arXiv

Measuring the masses of magnetic white dwarfs: A NuSTAR Legacy Survey

The hard X-ray spectrum of magnetic cataclysmic variables can be modelled to provide a measurement of white dwarf mass. This method is complementary to radial velocity measurements, which depend on the (typically rather uncertain) binary inclination. Here we present results from a Legacy Survey of 19 magnetic cataclysmic variables with NuSTAR. We fit accretion column models to their 20-78 keV spectra and derive the white dwarf masses, finding a weighted average $\bar{M}_{\rm WD}=0.77\pm0.02$ $M_{\odot}$, with a standard deviation $σ=0.10$ $M_{\odot}$, when we include the masses derived from previous NuSTAR observations of seven additional magnetic cataclysmic variables. We find that the mass distribution of accreting magnetic white dwarfs is consistent with that of white dwarfs in non-magnetic cataclysmic variables. Both peak at a higher mass than the distributions of isolated white dwarfs and post-common-envelope binaries. We speculate as to why this might be the case, proposing that consequential angular momentum losses may play a role in accreting magnetic white dwarfs and/or that our knowledge of how the white dwarf mass changes over accretion-nova cycles may also be incomplete.

preprint2016arXiv

Equation of state constraints for the cold dense matter inside neutron stars using the cooling tail method

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.

preprint2016arXiv

Radiative properties of magnetic neutron stars with metallic surfaces and thin atmospheres

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.

preprint2015arXiv

Measuring the basic parameters of neutron stars using model atmospheres

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.

preprint2015arXiv

The imprint of carbon combustion on a superburst from the accreting neutron star 4U 1636-536

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.

preprint2015arXiv

Transient X-ray pulsar V0332+53: pulse phase-resolved spectroscopy and the reflection model

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.

preprint2014arXiv

Relation between the X-ray and Optical Luminosities in Binary Systems with Accreting Nonmagnetic White Dwarfs

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

preprint2013arXiv

Carbon neutron star atmospheres

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.

preprint2012arXiv

Fitting the Chandra LETG spectrum of SS Cygni in outburst with model atmosphere spectra

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.

preprint2012arXiv

Investigation of the new cataclysmic variable 1RXS J180834.7+101041

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.

preprint2012arXiv

Investigations of a new eclipsing cataclysmic variable HBHA 4705-03

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

preprint2012arXiv

Magnetized neutron star atmospheres: beyond cold plasma approximation

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.

preprint2010arXiv

1RXS J180834.7+101041 is a new cataclysmic variable with non-uniform disc

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.

preprint2010arXiv

Dark spot, Spiral waves and the SW Sex behaviour: it is all about UX Ursae Majoris

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.

preprint2010arXiv

Modeling of the spectral energy distribution of the cataclysmic variable TT Ari and evaluation of the system parameters

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 α-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 ε(\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.

preprint2010arXiv

Population synthesis of DA white dwarfs: constraints on soft X-ray spectra evolution

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.

preprint2010arXiv

Quantum nature of cyclotron harmonics in thermal spectra of neutron stars

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.