Source author record

Juri Poutanen

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

ResearcherUnclaimed source record

astro-ph.HE astro-ph.GA astro-ph.SR astro-ph astro-ph.CO nucl-th astro-ph.IM gr-qc hep-th physics.plasm-ph

Catalog footprint

What is connected

77works

10topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2026arXiv

Discovery of High X-Ray Polarization from the Neutron Star Low-Mass X-Ray Binary Cyg X-2 in the Horizontal Branch

We present results from simultaneous X-ray polarimetric and spectroscopic observations of the bright neutron star low-mass X-ray binary Cyg X-2, performed by the Imaging X-ray Polarimetry Explorer (IXPE) and the Nuclear Spectroscopic Telescope Array (NuSTAR). IXPE detected significant polarization (15 sigma) from the source in the 2-8 keV energy band with an average polarization degree (PD) of 4.5% +/- 0.3% and a polarization angle (PA) of 128 +/- 2 degrees as the source moved along the horizontal branch of its Z-track. The PD increases with energy reaching 9.9% +/- 2.8% in the 7-8 keV band, with no evidence for energy-dependent variation in the PA. The PA is roughly consistent with previous measurements obtained during the normal and flaring branches and also with the known radio jet axis. From spectropolarimetric analysis, the main contribution to the polarized radiation is due to Comptonized photons, but the polarization is higher than predicted in typical spreading layer geometries. The observed high polarization may be due to a combination of a highly polarized reflected component and a moderately polarized spreading layer on the neutron star surface or produced by electron scattering in an equatorial wind.

preprint2023arXiv

Determination of X-ray pulsar geometry with IXPE polarimetry

Using observations of X-ray pulsar Her X-1 by the Imaging X-ray Polarimetry Explorer, we report on a highly significant ($>17σ$) detection of the polarization signal from an accreting neutron star. The observed degree of the linear polarization of $\sim$10\% is found to be far below theoretical expectations for this object, and stays low throughout the spin cycle of the pulsar. Both the polarization degree and the angle exhibit variability with pulse phase, which allowed us to measure the pulsar spin position angle 57(2) deg and the magnetic obliquity 12(4) deg, which is an essential step towards detailed modelling of the intrinsic emission of X-ray pulsars. Combining our results with the optical polarimetric data, we find that the spin axis of the neutron star and the angular momentum of the binary orbit are misaligned by at least $\sim$20 deg, which is a strong argument in support of the models explaining stability of the observed super-orbital variability with the precession of the neutron star.

preprint2022arXiv

Black hole spin-orbit misalignment in the X-ray binary MAXI J1820+070

The observational signatures of black holes in x-ray binary systems depend on their masses, spins, accretion rate and the misalignment angle between the black hole spin and the orbital angular momentum. We present optical polarimetric observations of the black hole x-ray binary MAXI J1820+070, from which we constrain the position angle of the binary orbital axis. Combining this with previous determinations of the relativistic jet orientation axis, which traces the black hole spin, and the inclination of the orbit, we determine a lower limit of 40 deg on the spin-orbit misalignment angle. The misalignment has to originate from either the binary or black hole formation stage. If other x-ray binaries have similarly large misalignments, these would bias measurements of black hole masses and spins from x-ray observations.

preprint2022arXiv

Constraints on the magnetic field structure in accreting compact objects from aperiodic variability

We investigate the aperiodic variability for a relatively large sample of accreting neutron stars and intermediate polars, focusing on the properties of the characteristic break commonly observed in power spectra of accreting objects. In particular, we investigate the relation of the break frequency and the magnetic field strength, both of which are connected to the size of the magnetosphere. We find that for the majority of objects in our sample the measured break frequency values indeed agree with estimated inner radii of the accretion disc, which allows to use observed break frequencies to independently assess the magnetic field strength and structure in accreting compact objects. As a special case, we focus on Hercules X-1 which is a persistent, medium-luminosity X-ray pulsar accreting from its low-mass companion. In the literature, it has been suggested that the complex pulse profiles, the spin-up behaviour and the luminosity-correlation of the cyclotron energy seen in Her X-1 can be explained with a complex magnetic field structure of the neutron star. Here, we connect the measured break frequency to the magnetospheric radius and show that the magnetic field strength derived assuming a dipole configuration is nearly an order of magnitude smaller than the magnetic field strength corresponding to the cyclotron energy. Accordingly, this discrepancy can be explained with the magnetic field having strong multipole components. The multipolar structure would also increase the accreting area on the neutron star surface, explaining why the critical luminosity for accretion column formation is puzzlingly high in this source.

preprint2022arXiv

Limits on X-ray Polarization at the Core of Centaurus A as Observed with the Imaging X-ray Polarimetry Explorer

We present measurements of the polarization of X-rays in the $2-8 \thinspace \mathrm{keV}$ band from the nucleus of the radio galaxy Centaurus A (Cen A), using a 100ks observation from the Imaging X-ray Polarimetry Explorer (IXPE). Nearly simultaneous observations of Cen A were also taken with the Swift, NuSTAR, and INTEGRAL observatories. No statistically significant degree of polarization is detected with IXPE. These observations have a minimum detectable polarization at $99 \%$ percent confidence (MDP$_{99}$) of $6.5 \%$ using a weighted, spectral model-independent calculation in the $2-8 \thinspace \mathrm{keV}$ band. The polarization angle $ψ$ is consequently unconstrained. Spectral fitting across three orders of magnitude in X-ray energy ($0.3-400 \thinspace \mathrm{keV}$) demonstrates that the SED of Cen A is well described by a simple power law with moderate intrinsic absorption ($N_H \sim 10^{23} \thinspace \mathrm{cm}^{-2}$) and a Fe K$α$ emission line, although a second unabsorbed power law is required to account for the observed spectrum at energies below $2 \thinspace \mathrm{keV}$. This spectrum suggests that the reprocessing material responsible for this emission line is optically thin and distant from the central black hole. Our upper limits on the X-ray polarization are consistent with the predictions of Compton scattering, although the specific seed photon population responsible for production of the X-rays cannot be identified. The low polarization degree, variability in the core emission, and the relative lack of variability in the Fe K$α$ emission line support a picture where electrons are accelerated in a region of highly disordered magnetic fields surrounding the innermost jet.

preprint2022arXiv

Optical polarization signatures of black hole X-ray binaries

Polarimetry provides an avenue for probing the geometry and physical mechanisms producing optical radiation in many astrophysical objects, including stellar binary systems. We present the results of multiwavelength (BVR) polarimetric studies of a sample of historical black hole X-ray binaries, observed during their outbursts or in the quiescent (or near-quiescent) state. We surveyed both long- and short-period systems, located at different Galactic latitudes. We performed careful analysis of the interstellar polarization in the direction on the sources to reliably estimate the intrinsic source polarization. Intrinsic polarization was found to be small (< 0.2 per cent) in sources observed in bright soft states (MAXI J0637-430 and 4U 1957+115). It was found to be significant in the rising hard state of MAXI J1820+070 at the level of 0.5 per cent and negligible in the decaying hard state and during its failed outbursts, while Swift J1357.2-0933 showed its absence in the rising hard state. Three (XTE J1118+480, V4641 Sgr, V404 Cyg) sources observed during quiescence show no evidence of significant intrinsic polarization, while MAXI J1820+070 is the only black hole X-ray binary which showed substantial (> 5 per cent) intrinsic quiescent-state polarization with a blue spectrum. The absence of intrinsic polarization at the optical wavelengths puts constraints on the potential contribution of non-stellar (jet, hot flow, accretion disc) components to the total spectra of quiescent black hole X-ray binaries.

preprint2022arXiv

Simultaneous space and phase resolved X-ray polarimetry of the Crab Pulsar and Nebula

The Crab pulsar and its nebula are among the most studied astrophysical systems, and constitute one of the most promising environments where high energy processes and particle acceleration can be investigated. They are the only objects for which previous X-ray polarisation has been reported. We present here the first Imaging X-ray Polarimetry Explorer (IXPE) observation of the Crab pulsar and nebula. The total pulsar pulsed emission in the [2-8] keV energy range is unpolarised. Significant polarisation up to 15% is detected only in the core of the main peak. The nebula has a total space integrated polarised degree of 20% and polarisation angle of 145deg. The polarised maps show a large variation in the local polarisation, and regions with polarised degree up to 45-50%. The polarisation pattern suggests a predominantly toroidal magnetic field.

preprint2021arXiv

Analytical techniques for polarimetric imaging of accretion flows in Schwarzschild metric

Emission from an accretion disc around compact objects, such as neutron stars and black holes, is expected to be significantly polarized. The polarization can be used to put constraints on geometrical and physical parameters of the compact sources -- their radii, masses and spins -- as well as to determine the orbital parameters. The radiation escaping from the innermost parts of the disc is strongly affected by the gravitational field of the compact object and relativistic velocities of the matter. The straightforward calculation of the observed polarization signatures involves computationally expensive ray-tracing technique. At the same time, having fast computational routines for direct data fitting becomes increasingly important in light of the currently observed images of the accretion flow around supermassive black hole in M87 by the Event Horizon Telescope, infrared polarization signatures coming from Sgr A*, as well as for the upcoming X-ray polarization measurements by the Imaging X-ray Polarimetry Explorer and enhanced X-ray Timing and Polarimetry mission. In this work, we obtain an exact analytical expression for the rotation angle of polarization plane in Schwarzschild metric accounting for the effects of light bending and relativistic aberration. We show that the calculation of the observed flux, polarization degree and polarization angle as a function of energy can be performed analytically with high accuracy using approximate light-bending formula, lifting the need for the pre-computed tabular models in fitting routines.

preprint2021arXiv

Broad-band analysis of X-ray pulsar 2S 1845$-$024

We present results of detailed investigation of the poorly studied X-ray pulsar 2S 1845$-$024 based on the data obtained with $NuSTAR$ observatory during the type I outburst in 2017. Neither pulse phase-averaged, nor phase-resolved spectra of the source show evidence for a cyclotron absorption feature. We also used the data obtained from other X-ray observatories ($Swift$, $XMM-Newton$ and $Chandra$) to study the spectral properties as a function of orbital phase. The analysis revealed a high hydrogen column density for the source reaching $\sim$10$^{24}$ cm$^{-2}$ around the periastron. Using high-quality $Chandra$ data we were able to obtain an accurate localization of 2S 1845$-$024 at R.A. = 18$^{h}$48$^{m}$16$^{s}$.8 and Dec. = $-$2$^{\circ}$25'25".1 (J2000) that allowed us to use infrared (IR) data to roughly classify the optical counterpart of the source as an OB supergiant at the distance of $\gtrsim$15 kpc.

preprint2021arXiv

Mechanical model of a boundary layer for the parallel tracks of kilohertz quasi-periodic oscillations in accreting neutron stars

Kilohertz-scale quasi-periodic oscillations (kHz QPOs) are a distinct feature of the variability of neutron star low-mass X-ray binaries. Among all the variability modes, they are especially interesting as a probe for the innermost parts of the accretion flow, including the accretion boundary layer (BL) on the surface of the neutron star. All the existing models of kHz QPOs explain only part of their rich phenomenology. Here, we show that some of their properties may be explained by a very simple model of the BL that is spun up by accreting rapidly rotating matter from the disk and spun down by the interaction with the neutron star. In particular, if the characteristic time scales for the mass and the angular momentum transfer from the BL to the star are of the same order of magnitude, our model naturally reproduces the so-called parallel tracks effect, when the QPO frequency is correlated with luminosity at time scales of hours but becomes uncorrelated at time scales of days. The closeness of the two time scales responsible for mass and angular momentum exchange between the BL and the star is an expected outcome of the radial structure of the BL.

preprint2021arXiv

Neutron star parameter constraints for accretion-powered millisecond pulsars from the simulated IXPE data

We have simulated the X-ray polarization data that can be obtained with the Imaging X-ray Polarimetry Explorer, when observing accretion-powered millisecond pulsars. We estimated the necessary exposure times for SAX J1808.4$-$3658 in order to obtain different accuracies in the measured time-dependent Stokes profiles integrated over all energy channels. We found that the measured relative errors depend strongly on the relative configuration of the observer and the emitting hotspot. The improvement in the minimum relative error in Stokes $Q$ and $U$ parameters as a function of observing time $t$ scales as $1/\sqrt{t}$, and spans the range from 30-90% with 200 ks exposure time to 20-60% with 500 ks exposure time (in case of data binned in 19 phase bins). The simulated data were also used to predict how accurate measurements of the geometrical parameters of the neutron star can be made when modelling only $Q$ and $U$ parameters, but not the flux. We found that the observer inclination and the hotspot co-latitude could be determined with better than 10 deg accuracy for most of the cases we considered. These measurements can be used to further constrain neutron star mass and radius when combined with modelling of the X-ray pulse profile.

preprint2020arXiv

Accurate analytic formula for light bending in Schwarzschild metric

We propose new analytic formulae describing light bending in Schwarzschild metric. For emission radii above the photon orbit at 1.5 Schwarzschild radius, the formulae have an accuracy of better than 0.2% for the bending angle and 3% for the lensing factor for any trajectories that turn around a compact object by less than about 160 deg. In principle, they can be applied to any emission point above the horizon of the black hole. The proposed approximation can be useful for problems involving emission from neutron stars and accretion discs around compact objects when fast accurate calculations of light bending are required. It can also be used to test the codes that compute light bending using exact expressions via elliptical integrals.

preprint2020arXiv

Colors and patterns of black hole X-ray binary GX 339-4

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.

preprint2020arXiv

Disc and wind in black hole X-ray binary MAXIJ1820+070 observed through polarized light during its 2018 outburst

We describe the first complete polarimetric dataset of the entire outburst of a low-mass black hole X-ray binary system and discuss the constraints for geometry and radiative mechanisms it imposes. During the decaying hard state, when the optical flux is dominated by the non-thermal component, the observed polarization is consistent with the interstellar values in all filters. During the soft state, the intrinsic polarization of the source is small, $\sim 0.15$ per cent in $B$ and $V$ filters, and is likely produced in the irradiated disc. A much higher polarization, reaching $\sim 0.5$ per cent in $V$ and $R$ filters, at position angle of $\sim 25^\circ$ observed in the rising hard state coincides in time with the detection of winds in the system. This angle coincides with the position angle of the jet. The detected optical polarization is best explained by scattering of the non-thermal (hot flow or jet base) radiation in an equatorial wind.

preprint2020arXiv

Discovery of a retrogradely rotating neutron star in the X-ray pulsar GX 301$-$2

We report on the analysis of the spin evolution of a slow X-ray pulsar GX 301$-$2 along the orbit using long-term monitoring by Fermi/GBM. Based on the observationally confirmed accretion scenario and an analytical model for the accretion of angular momentum we demonstrate that in this system, the neutron star spins retrogradely, that is, in a direction opposite to the orbital motion. This first-of-a-kind discovery of such a system proves the principal possibility of retrograde rotation in accreting systems with suitable accretion torque, and might have profound consequences for our understanding of the spin evolution of X-ray pulsars, estimates of their initial spin periods, and the ultimate result of their evolution.

preprint2020arXiv

Kilohertz quasi-periodic oscillations from neutron star spreading layers

When the accretion disc around a weakly magnetised neutron star (NS) meets the stellar surface, it should brake down to match the rotation of the NS, forming a boundary layer. As the mechanisms potentially responsible for this braking are apparently inefficient, it is reasonable to consider this layer as a spreading layer (SL) with negligible radial extent and structure. We perform hydrodynamical 2D spectral simulations of an SL, considering the disc as a source of matter and angular momentum. Interaction of new, rapidly rotating matter with the pre-existing, relatively slow material co-rotating with the star leads to instabilities capable of transferring angular momentum and creating variability on dynamical timescales. For small accretion rates, we find that the SL is unstable for heating instability that disrupts the initial latitudinal symmetry and produces large deviations between the two hemispheres. This instability also results in breaking of the axial symmetry as coherent flow structures are formed and escape from the SL intermittently. At enhanced accretion rates, the SL is prone to shearing instability and acts as a source of oblique waves that propagate towards the poles, leading to patterns that again break the axial symmetry. We compute artificial light curves of an SL viewed at different inclination angles. Most of the simulated light curves show oscillations at frequencies close to 1kHz. We interpret these oscillations as inertial modes excited by shear instabilities near the boundary of the SL. Their frequencies, dependence on flux, and amplitude variations can explain the high-frequency pair quasi-periodic oscillations observed in many low-mass X-ray binaries.

preprint2020arXiv

Magnetospheric return-current-heated atmospheres of rotation-powered millisecond pulsars

We computed 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 were computed self-consistently with the atmosphere model, instead of assuming a simplified deep-heated atmosphere in radiative equilibrium. We used exact Compton scattering formalism to model the properties of the emergent X-ray radiation. The deep-heating approximation was found to be valid only if most of the heat originates from ultra-relativistic bombarding particles with Lorentz factors of $γ\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: it produces a Compton-upscattered high-energy tail in the spectrum and alters 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. Finally, the connection between the energy distribution of the return current particles and the atmosphere emission properties offers a new tool to probe the exact physics of pulsar magnetospheres.

preprint2020arXiv

Observational appearance of rapidly rotating neutron stars: X-ray bursts, cooling tail method, and radius determination

Neutron stars (NSs) in low-mass X-ray binaries rotate at frequencies high enough to significantly deviate from sphericity ($ν_* \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 $ν_* =$ 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 $ν_* \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).

preprint2020arXiv

Spectral and temporal properties of Compton scattering by mildly relativistic thermal electrons

We have obtained new solutions and methods for the process of thermal Comptonization. We modify the solution to the kinetic equation of Sunyaev \& Titarchuk to allow its application up to mildly relativistic electron temperatures and optical depths $\gtrsim{1}$. The solution can be used for spectral fitting of X-ray spectra from astrophysical sources. We also have developed an accurate Monte Carlo method for calculating spectra and timing properties of thermal-Comptonization sources. The accuracy of our kinetic-equation solution is verified by comparison with the Monte Carlo results. We also compare our results with those of other publicly available methods. Furthermore, based on our Monte Carlo code, we present distributions of the photon emission times and the evolution of the average photon energy for both up and down-scattering.

preprint2020arXiv

The unusual behaviour of the young X-ray pulsar SXP 1062 during the 2019 outburst

We present results of the first dedicated observation of the young X-ray pulsar SXP 1062 in the broad X-ray energy band obtained during its 2019 outburst with the NuSTAR and XMM-Newton observatories. The analysis of the pulse-phase averaged and phase-resolved spectra in the energy band from 0.5 to 70 keV did not reveal any evidence for the presence of a cyclotron line. The spin period of the pulsar was found to have decreased to 979.48+/-0.06 s implying a ~10% reduction compared to the last measured period during the monitoring campaign conducted about five years ago, and is puzzling considering that the system apparently did not show major outbursts ever since. The switch of the pulsar to the spin-up regime supports the common assumption that torques acting on the accreting neutron star are nearly balanced and thus SXP 1062 likely also spins with a period close to the equilibrium value for this system. The current monitoring of the source revealed also a sharp drop of its soft X-ray flux right after the outburst, which is in drastic contrast to the behaviour during the previous outburst when the pulsar remained observable for years with only a minor flux decrease after the end of the outburst. This unexpected off state of the source lasted for at most 20 days after which SXP 1062 returned to the level observed during previous campaigns. We discuss this and other findings in context of the modern models of accretion onto strongly magnetized neutron stars.

preprint2016arXiv

Compton scattering S-matrix and cross section in strong magnetic field

Compton scattering of polarized radiation in a strong magnetic field is considered. The recipe for calculation of the scattering matrix elements, the differential and total cross sections based on quantum electrodynamic (QED) second order perturbation theory is presented for the case of arbitrary initial and final Landau level, electron momentum along the field and photon momentum. Photon polarization and electron spin state are taken into account. The correct dependence of natural Landau level width on the electron spin state is taken into account in general case of arbitrary initial photon momentum for the first time. A number of steps in calculations were simplified analytically making the presented recipe easy-to-use. The redistribution functions over the photon energy, momentum and polarization states are presented and discussed. The paper generalizes already known results and offers a basis for accurate calculation of radiation transfer in strong $B$-field, for example, in strongly magnetized neutron stars.

preprint2016arXiv

Gamma-ray opacity of the anisotropic stratified broad-line regions in blazars

The GeV-range spectra of blazars are shaped not only by non-thermal emission processes internal to the relativistic jet but also by external pair-production absorption on the thermal emission of the accretion disc and the broad-line region (BLR). For the first time, we compute here the pair-production opacities in the GeV range produced by a realistic BLR accounting for the radial stratification and radiation anisotropy. Using photoionization modelling with the CLOUDY code, we calculate a series of BLR models of different sizes, geometries, cloud densities, column densities and metallicities. The strongest emission features in the model BLR are Ly$α$ and HeII Ly$α$. Contribution of recombination continua is smaller, especially for hydrogen, because Ly continuum is efficiently trapped inside the large optical depth BLR clouds and converted to Lyman emission lines and higher-order recombination continua. The largest effects on the gamma-ray opacity are produced by the BLR geometry and localization of the gamma-ray source. We show that when the gamma-ray source moves further from the central source, all the absorption details move to higher energies and the overall level of absorption drops because of decreasing incidence angles between the gamma-rays and BLR photons. The observed positions of the spectral breaks can be used to measure the geometry and the location of the gamma-ray emitting region relative to the BLR. Strong dependence on geometry means that the soft photons dominating the pair-production opacity may be actually produced by a different population of BLR clouds than the bulk of the observed broad line emission.

preprint2016arXiv

Measuring the neutron star equation of state using X-ray timing

One of the primary science goals of the next generation of hard X-ray timing instruments is to determine the equation of state of the matter at supranuclear densities inside neutron stars, by measuring the radius of neutron stars with different masses to accuracies of a few percent. Three main techniques can be used to achieve this goal. The first involves waveform modelling. The flux we observe from a hotspot on the neutron star surface offset from the rotational pole will be modulated by the star's rotation, giving rise to a pulsation. Information about mass and radius is encoded into the pulse profile via relativistic effects, and tight constraints on mass and radius can be obtained. The second technique involves characterising the spin distribution of accreting neutron stars. The most rapidly rotating stars provide a very clean constraint, since the mass-shedding limit is a function of mass and radius. However the overall spin distribution also provides a guide to the torque mechanisms in operation and the moment of inertia, both of which can depend sensitively on dense matter physics. The third technique is to search for quasi-periodic oscillations in X-ray flux associated with global seismic vibrations of magnetars (the most highly magnetized neutron stars), triggered by magnetic explosions. The vibrational frequencies depend on stellar parameters including the dense matter equation of state. We illustrate how these complementary X-ray timing techniques can be used to constrain the dense matter equation of state, and discuss the results that might be expected from a 10m$^2$ instrument. We also discuss how the results from such a facility would compare to other astronomical investigations of neutron star properties. [Modified for arXiv]

preprint2016arXiv

Propeller effect in action in the ultraluminous accreting magnetar M82 X-2

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.

preprint2015arXiv

Discovery of correlated optical/X-ray quasi-periodic oscillations in black hole binary SWIFT J1753.5-0127

We report the discovery of the correlated optical/X-ray low-frequency quasi-periodic oscillations (QPOs) in black hole binary SWIFT J1753.5-0127. The phase lag between two light-curves at the QPO frequency is close to zero. This result puts strong constraints on the nature of the optical emission in this object and on the origin of the QPOs in general. We demonstrate that the QPO signal and the broadband variability can be explained in terms of the hot accretion flow radiating in both optical and X-ray bands. In this model, the QPO appears due to the Lense-Thirring precession of entire flow, while the broadband variability in the optical is produced by two components: the hot flow and the irradiated disc. Using the phase-lag spectra, we put a lower limit on the orbital inclination i>50 deg, which can be used to constrain the mass of the compact object.

preprint2015arXiv

Models of neutron star atmospheres enriched with nuclear burning ashes

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.

preprint2015arXiv

NuSTAR discovery of a cyclotron absorption line in the transient X-ray pulsar 2S 1553-542

We report results of a spectral and timing analysis of the poorly studied transient X-ray pulsar 2S 1553-542 using data collected with the NuSTAR and Chandra observatories and the Fermi/GBM instrument during an outburst in 2015. Properties of the source at high energies (>30 keV) are studied for the first time and the sky position had been essentially improved. The source broadband spectrum has a quite complicated shape and can be reasonably described by a composite model with two continuum components - a black body emission with the temperature about 1 keV at low energies and a power law with an exponential cutoff at high energies. Additionally an absorption feature at $\sim23.5$ keV is discovered both in phase-averaged and phase-resolved spectra and interpreted as the cyclotron resonance scattering feature corresponding to the magnetic field strength of the neutron star $B\sim3\times10^{12}$ G. Based on the Fermi/GBM data the orbital parameters of the system were substantially improved, that allowed us to determine the spin period of the neutron star P = 9.27880(3) s and a local spin-up $\dot P \simeq -7.5\times10^{-10}$ s s$^{-1}$ due to the mass accretion during the NuSTAR observations. Assuming accretion from the disk and using standard torque models we have estimated the distance to the system $d=20\pm4$ kpc.

preprint2015arXiv

On the maximum accretion luminosity of magnetized neutron stars: connecting X-ray pulsars and ultraluminous X-ray sources

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.

preprint2015arXiv

Polarization modulation from Lense-Thirring precession in X-ray binaries

It has long been recognised that quasi-periodic oscillations (QPOs) in the X-ray light curves of accreting black hole and neutron star binaries have the potential to be powerful diagnostics of strong field gravity. However, this potential cannot be fulfilled without a working theoretical model, which has remained elusive. Perhaps the most promising model associates the QPO with Lense-Thirring precession of the inner accretion flow, with the changes in viewing angle and Doppler boosting modulating the flux over the course of a precession cycle. Here, we consider the polarization signature of a precessing inner accretion flow. We use simple assumptions about the Comptonization process generating the emitted spectrum and take all relativistic effects into account, parallel transporting polarization vectors towards the observer along null geodesics in the Kerr metric. We find that both the degree of linear polarization and the polarization angle should be modulated on the QPO frequency. We calculate the predicted absolute rms variability amplitude of the polarization degree and angle for a specific model geometry. We find that it should be possible to detect these modulations for a reasonable fraction of parameter space with a future X-ray polarimeter such as NASA's Polarization Spectroscopic Telescope Array (PolSTAR: the satellite incarnation of the recent balloon experiment X-Calibur).

preprint2015arXiv

Positive correlation between the cyclotron line energy and luminosity in sub-critical X-ray pulsars: Doppler effect in the accretion channel

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.

preprint2015arXiv

Reprocessing model for the optical quasi-periodic oscillations in black hole binaries

A number of black hole X-ray transients show quasi-periodic oscillations (QPOs) in the optical (ultraviolet) and X-ray bands at the same frequency, which challenge models for production of radiation at these wavelengths. We propose a model where the optical radiation is modulated by the oscillating X-ray flux resulting in varying irradiation of the outer parts of the accretion disc. The proposed QPO mechanism inevitably takes place in the systems with sufficiently small ratio of the outer disc radius to the QPO period. We show that, unlike in the case of the aperiodic variability, it is not possible to obtain the optical QPO profiles from those observed in the X-rays through the transfer function, because of different X-ray signals seen by the disc and by the observer. We demonstrate that with the increasing QPO frequency, occurring at the rising phase of the X-ray outburst, the rms should be constant for sufficiently low frequencies, then to increase reaching the peak and finally to drop substantially when the QPO period becomes comparable to the light-crossing time to the outer disc. We predict that the QPO rms in this model should increase towards shorter wavelengths and this fact can be used to distinguish it from other QPO mechanisms.

preprint2014arXiv

Colours of black holes: infrared flares from the hot accretion disc in XTE J1550-564

Outbursts of the black hole (BH) X-ray binaries are dramatic events occurring in our Galaxy approximately once a year. They are detected by the X-ray telescopes and often monitored at longer wavelengths. We analyse the X-ray and optical/infrared (OIR) light curves of the BH binary XTE J1550-564 during the 2000 outburst. By using the observed extreme colours as well as the characteristic decay time-scales of the OIR and X-ray light curves, we put strong constraints on the extinction towards the source. We accurately separate the contributions to the OIR flux of the irradiated accretion disc and a non-thermal component. We show that the OIR non-thermal component appears during the X-ray state transitions both during the rising and the decaying part of the outburst at nearly the same X-ray hardness but at luminosities differing by a factor of 3. The line marking the quenching/recovery of the non-thermal component at the X-ray hardness-flux diagram seems to coincide with the "jet line" that marks the presence of the compact radio jet. The inferred spectral shape and the evolution of the non-thermal component during the outburst, however, are not consistent with the jet origin, but are naturally explained in terms of the hybrid hot flow scenario, where non-thermal electrons emit synchrotron radiation in the OIR band. This implies a close, possibly causal connection between the presence of the hot flow and the compact jet. We find that the non-thermal component is hardening during the hard state at the decaying stage of the outburst, which indicates that the acceleration efficiency is a steep function of radius at low accretion rate.

preprint2014arXiv

Spectroscopic evidence for a low-mass black hole in SWIFT J1753.5-0127

The black hole (BH) candidate SWIFT J1753.5-0127 has remained active since the onset of its 2005 outburst. Emission lines in the optical spectrum were observed at the very beginning of the outburst, but since then the spectrum has been featureless making a precise BH mass estimation impossible. Here we present results from our optical and UV observations of SWIFT J1753.5-0127 taken in 2012-2013. Our new observations show extremely broad, double-peaked emission lines in the optical and UV spectra. The optical data also show narrow absorption and emission features with nearly synchronous and significant Doppler motions. A radial velocity study of these lines which we associate with the secondary star, yields a semi-amplitude of K_2=382 km/s. A time-series analysis of the spectral and photometric data revealed a possible orbital periodicity of 2.85 h, significantly shorter than the reported 3.2 h periodic signal by Zurita et al. The observed variability properties argue against a low orbital inclination angle and we present several observational arguments in favour of the BH interpretation. However, the measured radial velocity semi-amplitude of the donor star and the short orbital period imply that SWIFT J1753.5-0127 has one of the lowest measured mass function for a BH in a low-mass X-ray binary. We show that the compact object mass in excess of 5 Msun is highly improbable. Thus, SWIFT J1753.5-0127 is a BH binary that has one of the shortest orbital period and hosts probably one of the smallest stellar-mass BH found to date.

preprint2014arXiv

The critical accretion luminosity for magnetized neutron stars

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.

preprint2014arXiv

The effect of accretion on the measurement of neutron star mass and radius in the low-mass X-ray binary 4U 1608-52

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.

preprint2014arXiv

The influence of accretion geometry on the spectral evolution during thermonuclear (type-I) X-ray bursts

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.

preprint2014arXiv

The mystery of spectral breaks: Lyman continuum absorption by photon-photon pair production in the Fermi GeV spectra of bright blazars

We reanalyze Fermi/LAT gamma-ray spectra of bright blazars with a higher photon statistics than in previous works and with new Pass 7 data representation. In the spectra of the brightest blazar 3C 454.3 and possibly of 4C +21.35 we detect breaks at 5 GeV (in the rest frame) associated with the photon-photon pair production absorption by He II Lyman continuum (LyC). We also detect confident breaks at 20 GeV associated with hydrogen LyC both in the individual spectra and in the stacked redshift-corrected spectrum of several bright blazars. The detected breaks in the stacked spectra univocally prove that they are associated with atomic ultraviolet emission features of the quasar broad-line region (BLR). The dominance of the absorption by hydrogen Ly complex over He II, rather small detected optical depth, and the break energy consistent with the head-on collisions with LyC photons imply that the gamma-ray emission site is located within the BLR, but most of the BLR emission comes from a flat disk-like structure producing little opacity. Alternatively, the LyC emission region size might be larger than the BLR size measured from reverberation mapping, and/or the gamma-ray emitting region is extended. These solutions would resolve a long-standing issue how the multi-hundred GeV photons can escape from the emission zone without being absorbed by softer photons.

preprint2013arXiv

A reflection model for the cyclotron lines in the spectra of X-ray pulsars

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.

preprint2013arXiv

A unified Lense-Thirring precession model for optical and X-ray quasi-periodic oscillations in black hole binaries

Recent observations of accreting black holes reveal the presence of quasi-periodic oscillations (QPO) in the optical power density spectra. The corresponding oscillation periods match those found in the X-rays, implying a common origin. Among the numerous suggested X-ray QPO mechanisms, some may also work in the optical. However, their relevance to the broadband -- optical through X-ray -- spectral properties have not been investigated. For the first time, we discuss the QPO mechanism in the context of the self-consistent spectral model. We propose that the QPOs are produced by Lense-Thirring precession of the hot accretion flow, whose outer parts radiate in the optical wavelengths. At the same time, its innermost parts are emitting the X-rays, explaining the observed connection of QPO periods. We predict that the X-ray and optical QPOs should be either in phase or shifted by half a period, depending on the observer position. We investigate the QPO harmonic content and find that the variability amplitudes at the fundamental frequency are larger in the optical, while the X-rays are expected to have strong harmonics. We then discuss the QPO spectral dependence and compare the expectations to the existing data.

preprint2013arXiv

Hot accretion flow in black hole binaries: a link connecting X-rays to the infrared

Multiwavelength observations of Galactic black hole transients have opened a new path to understanding the physics of the innermost parts of the accretion flows. While the processes giving rise to their X-ray continuum have been studied extensively, the emission in the optical and infrared (OIR) energy bands was less investigated and remains poorly understood. The standard accretion disc, which may contribute to the flux at these wavelengths, is not capable of explaining a number of observables: the infrared excesses, fast OIR variability and a complicated correlation with the X-rays. It was suggested that these energy bands are dominated by the jet emission, however, this scenario does not work in a number of cases. We propose here an alternative, namely that most of the OIR emission is produced by the extended hot accretion flow. In this scenario, the OIR bands are dominated by the synchrotron radiation from the non-thermal electrons. An additional contribution is expected from the outer irradiated part of the accretion disc heated by the X-rays. We discuss properties of the model and compare them to the data. We show that the hot-flow scenario is consistent with many of the observed spectral data, at the same time naturally explaining X-ray timing properties, fast OIR variability and its correlation with the X-rays.

preprint2013arXiv

Modelling spectral and timing properties of accreting black holes: the hybrid hot flow paradigm

The general picture that emerged by the end of 1990s from a large set of optical and X-ray, spectral and timing data was that the X-rays are produced in the innermost hot part of the accretion flow, while the optical/infrared (OIR) emission is mainly produced by the irradiated outer thin accretion disc. Recent multiwavelength observations of Galactic black hole transients show that the situation is not so simple. Fast variability in the OIR band, OIR excesses above the thermal emission and a complicated interplay between the X-ray and the OIR light curves imply that the OIR emitting region is much more compact. One of the popular hypotheses is that the jet contributes to the OIR emission and even is responsible for the bulk of the X-rays. However, this scenario is largely ad hoc and is in contradiction with many previously established facts. Alternatively, the hot accretion flow, known to be consistent with the X-ray spectral and timing data, is also a viable candidate to produce the OIR radiation. The hot-flow scenario naturally explains the power-law like OIR spectra, fast OIR variability and its complex relation to the X-rays if the hot flow contains non-thermal electrons (even in energetically negligible quantities), which are required by the presence of the MeV tail in Cyg X-1. The presence of non-thermal electrons also lowers the equilibrium electron temperature in the hot flow model to <100 keV, making it more consistent with observations. Here we argue that any viable model should simultaneously explain a large set of spectral and timing data and show that the hybrid (thermal/non-thermal) hot flow model satisfies most of the constraints.

preprint2013arXiv

On the association of the ultraluminous X-ray sources in the Antennae galaxies with young stellar clusters

The nature of the ultra-luminous X-ray sources (ULXs) in the nearby galaxies is a matter of debates. One of the popular hypothesis associates them with accretion at a sub-Eddington rate on to intermediate mass black holes. Another possibility is a stellar-mass black hole in a high-mass X-ray binary accreting at super-Eddington rates. In this paper we find a highly significant association between brightest X-ray sources in the Antennae galaxies and stellar clusters. On the other hand, we show that most of the X-ray sources are located outside of these clusters. We study clusters associated with the ULXs using the ESO Very Large Telescope spectra and the Hubble Space Telescope data together with the theoretical evolutionary tracks and determine their ages to be below 6 Myr. This implies that the ULX progenitor masses certainly exceed 30 and for some objects are closer to 100 solar masses. We also estimate the ages of clusters situated close to the less luminous X-ray sources (with luminosity in the range 3x10^38 < L_X < 10^39 erg/s) and find that most of them are younger than 10 Myr, because they are surrounded by strong H$α$ emission. These findings are consistent with the idea that majority of ULXs are massive X-ray binaries that have been ejected in the process of formation of stellar clusters by a few-body encounters and at the same time rules out the proposal that most of the ULXs are intermediate mass black holes.

preprint2013arXiv

On the origin of cyclotron lines in the spectra of X-ray pulsars

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.

preprint2013arXiv

On the spreading layer emission in luminous accreting neutron stars

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.

preprint2013arXiv

White Paper for Blazar Observations with a GEMS-like X-ray Polarimetry Mission

In this document, we describe the scientific potential of blazar observations with a X-ray polarimetry mission like GEMS (Gravity and Extreme Magnetism SMEX). We describe five blazar science investigations that such a mission would enable: (i) the structure and the role of magnetic fields in AGN jets, (ii) analysis of the polarization of the synchrotron X-ray emission from AGN jets, (iii) discrimination between synchrotron self-Compton and external Compton models for blazars with inverse Compton emission in the X-ray band, (iv) a precision study of the polarization properties of the X-ray emission from Cen-A, (v) tests of Lorentz Invariance based on X-ray polarimetric observations of blazars. We conclude with a discussion of a straw man observation program and recommended accompanying multiwavelength observations.

preprint2013arXiv

White Paper on GEMS Study of Polarized X-rays from Neutron Stars

We examine the expected X-ray polarization properties of neutron-star X-ray sources of various types, e.g., accretion and rotation powered pulsars, magnetars, and low-mass X-ray binaries. We summarize the model calculations leading to these expected properties. We describe how a comparison of these with their observed properties, as inferred from GEMS data, will probe the essential dynamical, electromagnetic, plasma, and emission processes in neutron-star binaries, discriminate between models of these processes, and constrain model parameters. An exciting goal is the first observational demonstration in this context of the existence of vacuum resonance, a fundamental quantum electrodynamical phenomenon first described in the 1930s.

preprint2013arXiv

X-ray Polarization from Black Holes: GEMS Scientific White Paper

We present here a summary of the scientific goals behind the Gravity and Extreme Magnetism SMEX (GEMS) X-ray polarimetry mission's black hole (BH) observing program. The primary targets can be divided into two classes: stellar-mass galactic BHs in accreting binaries, and super-massive BHs in the centers of active galactic nuclei (AGN). The stellar-mass BHs can in turn be divided into various X-ray spectral states: thermal-dominant (disk), hard (radio jet), and steep power-law (hot corona). These different spectral states are thought to be generated by different accretion geometries and emission mechanisms. X-ray polarization is an ideal tool for probing the geometry around these BHs and revealing the specific properties of the accreting gas.

preprint2013arXiv

XIPE: the X-ray Imaging Polarimetry Explorer

X-ray polarimetry, sometimes alone, and sometimes coupled to spectral and temporal variability measurements and to imaging, allows a wealth of physical phenomena in astrophysics to be studied. X-ray polarimetry investigates the acceleration process, for example, including those typical of magnetic reconnection in solar flares, but also emission in the strong magnetic fields of neutron stars and white dwarfs. It detects scattering in asymmetric structures such as accretion disks and columns, and in the so-called molecular torus and ionization cones. In addition, it allows fundamental physics in regimes of gravity and of magnetic field intensity not accessible to experiments on the Earth to be probed. Finally, models that describe fundamental interactions (e.g. quantum gravity and the extension of the Standard Model) can be tested. We describe in this paper the X-ray Imaging Polarimetry Explorer (XIPE), proposed in June 2012 to the first ESA call for a small mission with a launch in 2017 but not selected. XIPE is composed of two out of the three existing JET-X telescopes with two Gas Pixel Detectors (GPD) filled with a He-DME mixture at their focus and two additional GPDs filled with pressurized Ar-DME facing the sun. The Minimum Detectable Polarization is 14 % at 1 mCrab in 10E5 s (2-10 keV) and 0.6 % for an X10 class flare. The Half Energy Width, measured at PANTER X-ray test facility (MPE, Germany) with JET-X optics is 24 arcsec. XIPE takes advantage of a low-earth equatorial orbit with Malindi as down-link station and of a Mission Operation Center (MOC) at INPE (Brazil).

preprint2013arXiv

Young rotation-powered pulsars as ultra-luminous X-ray sources

The aim of the present paper is to investigate a possible contribution of the rotation-powered pulsars and pulsar wind nebulae to the population of ultraluminous X-ray sources (ULXs). We first develop an analytical model for the evolution of the distribution function of pulsars over the spin period and find both the steady-state and the time-dependent solutions. Using the recent results on the X-ray efficiency dependence on pulsar characteristic age, we then compute the X-ray luminosity function (XLF) of rotation-powered pulsars. In a general case it has a broken power-law shape with a high luminosity cutoff, which depends on the distributions of the birth spin period and the magnetic field. Using the observed XLF of sources in the nearby galaxies and the condition that the pulsar XLF does not exceed that, we find the allowed region for the parameters describing the birth period distribution. We find that the mean pulsar period should be greater than 10-40 ms. These results are consistent with the constraints obtained from the X-ray luminosity of core-collapse supernovae. We estimate that the contribution of the rotation-powered pulsars to the ULX population is at a level exceeding 3 per cent. For a wide birth period distribution, this fraction grows with luminosity and above 10E40 erg/s pulsars can dominate the ULX population.

preprint2012arXiv

Joint spectral-timing modelling of the hard lags in GX 339-4: constraints on reflection models

The X-ray variations of hard state black hole X-ray binaries above 2 keV show 'hard lags', in that the variations at harder energies follow variations at softer energies, with a time-lag τdepending on frequency νapproximately as τ\propto ν^{-0.7}. Several models have so far been proposed to explain this time delay, including fluctuations propagating through an accretion flow, spectral variations during coronal flares, Comptonisation in the extended hot corona or a jet, or time-delays due to large-scale reflection from the accretion disc. In principle these models can be used to predict the shape of the energy spectrum as well as the frequency-dependence of the time-lags, through the construction of energy-dependent response functions which map the emission as a function of time-delay in the system. Here we use this approach to test a simple reflection model for the frequency-dependent lags seen in the hard state of GX 339-4, by simultaneously fitting the model to the frequency-dependent lags and energy spectrum measured by XMM-Newton in 2004 and 2009. Our model cannot simultaneously fit both the lag and spectral data, since the relatively large lags require an extremely flared disc which subtends a large solid angle to the continuum at large radii, in disagreement with the observed Fe Kαemission. Therefore, we consider it more likely that the lags > 2 keV are caused by propagation effects in the accretion flow, possibly related to the accretion disc fluctuations which have been observed previously.

preprint2012arXiv

Relativistic kinetic equation for Compton scattering of polarized radiation in strong magnetic field

We derive the relativistic kinetic equation for Compton scattering of polarized radiation in strong magnetic field using the Bogolyubov method. The induced scattering and the Pauli exclusion principle are taken into account. The electron polarization is also considered in the general form of the kinetic equation. The special forms of the equation for the cases of the non-polarized electrons, the rarefied electron gas and the two polarization mode description of radiation are found. The derived equations are valid for any photon and electron energies and the magnetic field strength below about 10^{16} G. These equations provide the basis for formulation of the equation for polarized radiation transport in atmospheres and magnetospheres of strongly magnetized neutron stars.

preprint2011arXiv

A Complexity-Brightness Correlation in Gamma Ray Bursts

We observe strong correlations between the temporal properties of gamma ray bursts (GRBs) and their apparent peak brightness. The strongest effect (with a significance level of 10^{-6}) is the difference between the brightness distributions of simple bursts (dominated by a single smooth pulse) and complex bursts (consisting of overlapping pulses). The latter has a break at a peak flux of 1.5 ph/cm^2/s, while the distribution of simple bursts is smooth down to the BATSE threshold. We also observe brightness dependent variations in the shape of the average peak aligned time profile (ATP) of GRBs. The decaying slope of the ATP shows time dilation when comparing bright and dim bursts while the rising slope hardly changes. Both slopes of the ATP are deformed for weak bursts as compared to strong bursts. The interpretation of these effects is simple: a complex burst where a number of independent pulses overlap in time appears intrinsically stronger than a simple burst. Then the BATSE sample of complex bursts covers larger redshifts where some cosmological factor causes the break in the peak brightness distribution. This break could correspond to the peak in the star formation rate that was recently shown to occur at a redshift of z~1.5.

preprint2011arXiv

A self-consistent hybrid Comptonization model for broad-band spectra of accreting supermassive black holes

The nature of the broad-band spectra of supermassive accreting black holes in active galactic nuclei (AGNs) is still unknown. The hard X-ray spectra of Seyferts as well as of Galactic stellar-mass black holes (GBHs) are well represented by thermal Comptonization, but the origin of the seed photons is less certain. The MeV tails observed in GBHs provide evidence in favour of non-thermal electron tails and it is possible that such electrons are also present in the X-ray emitting regions of AGNs. Using simulations with the kinetic code that self-consistently models electron and photon distributions, we find that the power-law-like X-ray spectra in AGNs can be explained in terms of the synchrotron self-Compton radiation of hybrid thermal/non-thermal electrons, similarly to the hard/low state of GBHs. Under a very broad range of parameters the model predicts a rather narrow distribution of photon spectral slopes consistent with that observed from LINERs and Seyferts at luminosities less than 3 per cent of the Eddington luminosity. The entire infrared to X-ray spectrum of these objects can be described in terms of our model, suggesting a tight correlation between the two energy bands. We show that the recently found correlation between slope and the Eddington ratio at higher luminosities can be described by the increasing fraction of disc photons in the emitting region, which may be associated with the decreasing inner radius of the optically thick accretion disc. The increasing flux of soft photons is also responsible for the transformation of the electron distribution from nearly thermal to almost completely non-thermal. The softer X-ray spectra observed in narrow-line Seyfert galaxies may correspond to non-thermal Comptonization of the disc photons, predicting that no cutoff should be observed up to MeV energies in these sources, similarly to the soft-state GBHs.

preprint2011arXiv

A synchrotron self-Compton -- disk reprocessing model for optical/X-ray correlation in black hole X-ray binaries

Physical picture of the emission mechanisms operating in the X-ray binaries was put under question by the simultaneous optical/X-ray observations with high time resolution. The light curves of the two energy bands appeared to be connected and the cross-correlation functions observed in three black hole binaries exhibited a complicated shape. They show a dip of the optical emission a few seconds before the X-ray peak and the optical flare just after the X-ray peak. This behavior could not be explained in terms of standard optical emission candidates (e.g., emission from the cold accretion disk or a jet). We propose a novel model, which explains the broadband optical to the X-ray spectra and the variability properties. We suggest that the optical emission consists of two components: synchrotron radiation from the non-thermal electrons in the hot accretion flow and the emission produced by reprocessing of the X-rays in the outer part of the accretion disk. The first component is anti-correlated with the X-rays, while the second one is correlated, but delayed and smeared relative to the X-rays. The interplay of the components explains the complex shape of the cross-correlation function, the features in the optical power spectral density as well as the time lags.

preprint2011arXiv

Fermi Observations of Blazars: Implications for Gamma-ray Production

The brightest blazars detected by the Fermi Gamma-ray Space Telescope Large Area Telescope (Fermi/LAT) show significant breaks in their spectra at a few GeV. The sharpness and the position of the breaks can be well reproduced by absorption of $γ$-rays via photon--photon pair production on He ii and H i Lyman recombination continua (LyC) produced in the broad-line region (BLR). Using 138 weeks of LAT observations of the brightest GeV blazar 3C 454.3 we find a power-law dependence of the peak energy on flux and discover anti-correlation between flux and the column density of the He ii LyC which is responsible for absorption of the >2.5 GeV photons in this object. The strength and the variability of the absorption implies the location of the gamma-ray emitting zone close to the boundary of the high-ionization part of the BLR and moving away from the black hole when the flux increases. A combination of the GeV breaks with the detection of a few powerful blazars in the TeV range puts strong constraints on the BLR size. Additional spectral breaks at ~100 and ~400 GeV due to absorption by the Balmer and Paschen lines could be detected by the Cherenkov Telescope Array.

preprint2011arXiv

Gamma-ray bursts from magnetized collisionally-heated jets

Jets producing gamma-ray bursts (GRBs) are likely to carry a neutron component that drifts with respect to the proton component. The neutron-proton collisions strongly heat the jet and generate electron-positron pairs. We investigate radiation produced by this heating using a new numerical code. Our results confirm the recent claim that collisional heating generates the observed Band-type spectrum of GRBs. We extend the model to study the effects of magnetic fields on the emitted spectrum. We find that the spectrum peak remains near 1 MeV for the entire range of the magnetization parameter $0<ε_B<2$ that is explored in our simulations. The low-energy part of the spectrum softens with increasing $ε_B$, and a visible soft excess appears in the keV band. The high-energy part of the spectrum extends well above the GeV range and can contribute to the prompt emission observed by Fermi/LAT. Overall, the radiation spectrum created by the collisional mechanism appears to agree with observations, with no fine-tuning of parameters.

preprint2011arXiv

Neutron star stiff equation of state derived from cooling phases of the X-ray burster 4U 1724-307

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]

preprint2011arXiv

Variation of the gamma-gamma opacity by the He II Lyman continuum constrains the location of the gamma-ray emission region in the blazar 3C 454.3

We study spectral properties of the brightest gamma-ray blazar 3C454.3 using 138 weeks of observations by the Fermi Gamma-ray Space Telescope (Fermi). We probe the behaviour of the source as a function of time at different brightness levels. The Fermi spectra in the GeV range can be well described by a wide underlying lognormal distribution with the photon-photon absorption breaks produced by the He II and H I Lyman recombination continua (LyC). We find a power-law dependence of the peak energy on flux and discover anti-correlation between the column density of the He II LyC and flux. This implies that the gamma-ray emission zone lies close to the boundary of the high-ionization part of the broad-line region and moves away from the black hole when the flux increases. Identification of the gamma-ray production with the relativistic jet, implies that the jet is already accelerated at sub-parsec distances from the central black hole, which favours the Blandford-Znajek process as the jet launching mechanism.

preprint2011arXiv

Varying disc-magnetosphere coupling as the origin of pulse profile variability in SAX J1808.4-3658

Accreting millisecond pulsars show significant variability of their pulse profiles, especially at low accretion rates. On the other hand, their X-ray spectra are remarkably similar with not much variability over the course of the outbursts. For the first time, we have discovered that during the 2008 outburst of SAX J1808.4-3658 a major pulse profile change was accompanied by a dramatic variation of the disc luminosity at almost constant total luminosity. We argue that this phenomenon is related to a change in the coupling between the neutron star magnetic field and the accretion disc. The varying size of the pulsar magnetosphere can influence the accretion curtain geometry and affect the shape and the size of the hotspots. Using this physical picture, we develop a self-consistent model that successfully describes simultaneously the pulse profile variation as well as the spectral transition. Our findings are particularly important for testing the theories of accretion onto magnetized neutron stars, better understanding of the accretion geometry as well as the physics of disc-magnetosphere coupling. The identification that varying hotspot size can lead to pulse profile changes has profound implications for determination of the neutron star masses and radii.

preprint2010arXiv

Constraining compactness and magnetic field geometry of X-ray pulsars from the statistics of their pulse profiles

The light curves observed from X-ray pulsars and magnetars reflect the radiation emission pattern, the geometry of the magnetic field, and the neutron star compactness. We study the statistics of X-ray pulse profiles in order to constrain the neutron star compactness and the magnetic field geometry. We collect the data for 124 X-ray pulsars, which are mainly in high-mass X-ray binary systems, and classify their pulse profiles according to the number of observed peaks seen during one spin period, dividing them into two classes, single- and double-peaked. We find that the pulsars are distributed about equally between both groups. We also compute the probabilities predicted by the theoretical models of two antipodal point-like spots that emit radiation according to the pencil-like emission patterns. These are then compared to the observed fraction of pulsars in the two classes. Assuming a blackbody emission pattern, it is possible to constrain the neutron star compactness if the magnetic dipole has arbitrary inclinations to the pulsar rotational axis. More realistic pencil-beam patterns predict that 79% of the pulsars are double-peaked independently of their compactness. The theoretical predictions can be made consistent with the data if the magnetic dipole inclination to the rotational axis has an upper limit of 40+/-4 deg. We also discuss the effect of limited sensitivity of the X-ray instruments to detect weak pulses, which lowers the number of detected double-peaked profiles and makes the theoretical predictions to be consistent with the data even if the magnetic dipole does have random inclinations. This shows that the statistics of pulse profiles does not allow us to constrain the neutron star compactness. In contrast to the previous claims by Bulik et al. (2003), the data also do not require the magnetic inclination to be confined in a narrow interval.

preprint2010arXiv

GeV breaks in blazars as a result of gamma-ray absorption within the broad-line region

Spectra of the brightest blazars detected by the Fermi Gamma-ray Space Telescope Large Area Telescope cannot be described by a simple power law model. A much better description is obtained with a broken power law, with the break energies of a few GeV. We show here that the sharpness and the position of the breaks can be well reproduced by absorption of gamma-rays via photon--photon pair production on HeII Lyman recombination continuum and lines. This implies that the blazar zone lies inside the region of the highest ionization of the broad-line region (BLR) within a light-year from a super-massive black hole. The observations of gamma-ray spectral breaks open a way of studying the BLR photon field in the extreme-UV/soft X-rays, which are otherwise hidden from our view.

preprint2010arXiv

Theory of Compton scattering by anisotropic electrons

Compton scattering plays an important role in various astrophysical objects such as accreting black holes and neutron stars, pulsars, and relativistic jets, clusters of galaxies as well as the early Universe. In most of the calculations it is assumed that the electrons have isotropic angular distribution in some frame. However, there are situations where the anisotropy may be significant due to the bulk motions, or anisotropic cooling by synchrotron radiation, or anisotropic source of seed soft photons. We develop here an analytical theory of Compton scattering by anisotropic distribution of electrons that can simplify significantly the calculations. Assuming that the electron angular distribution can be represented by a second order polynomial over cosine of some angle (dipole and quadrupole anisotropy), we integrate the exact Klein-Nishina cross-section over the angles. Exact analytical and approximate formulae valid for any photon and electron energies are derived for the redistribution functions describing Compton scattering of photons with arbitrary angular distribution by anisotropic electrons. The analytical expressions for the corresponding photon scattering cross-section on such electrons as well as the mean energy of scattered photons, its dispersion and radiation pressure force are also derived. We applied the developed formalism to the accurate calculations of the thermal and kinematic Sunyaev-Zeldovich effects for arbitrary electron distributions.

preprint2009arXiv

Accreting millisecond pulsar SAX J1808.4-3658 during its 2002 outburst: evidence for a receding disc

An outburst of the accreting X-ray millisecond pulsar SAX J1808.4-3658 in October-November 2002 was followed by the Rossi X-ray Timing Explorer for more than a month. We demonstrate how the area covered by the hotspot at the neutron star surface is decreasing in the course of the outburst together with the reflection amplitude. These trends are in agreement with the natural scenario, where the disc inner edge is receding from the neutron star as the mass accretion rate drops. These findings are further supported by the variations of the pulse profiles, which clearly show the presence of the secondary maximum at the late stages of the outburst after October 29. This fact can be interpreted as the disc receding sufficiently far from the neutron star to open the view of the lower magnetic pole. In that case, the disc inner radius can be estimated. Assuming that disc is truncated at the Alfvén radius, we constrain the stellar magnetic moment to μ=(9\pm5) 10^{25} G cm^3, which corresponds to the surface field of 10^8 G. On the other hand, using the magnetic moment recently obtained from the observed pulsar spin-down rate we show that the disc edge has to be within factor of two of the Alfvén radius, putting interesting constraints on the models of the disc-magnetosphere interaction. We also demonstrate that the sharp changes in the phase of the fundamental are intimately related to the variations of the pulse profile, which we associate with the varying obscuration of the antipodal spot. The pulse profile amplitude allows us to estimate the colatitude of the hotspot centroid to be 4-10 deg.

preprint2009arXiv

On the nature of pulse profile variations and timing noise in accreting millisecond pulsars

Timing noise in the data on accretion-powered millisecond pulsars (AMP) appears as irregular pulse phase jumps on timescales from hours to weeks. A large systematic phase drift is also observed in the first discovered AMP SAX J1808.4-3658. To study the origin of these timing features, we use here the data of the well studied 2002 outburst of SAX J1808.4-3658. We develop first a model for pulse profile formation accounting for the screening of the antipodal emitting spot by the accretion disk. We demonstrate that the variations of the visibility of the antipodal spot associated with the receding accretion disk cause a systematic shift in Fourier phases, observed together with the changes in the pulse form. We show that a strong secondary maximum can be observed only in a narrow intervals of inner disk radii, which explains the very short appearance of the double-peaked profiles in SAX J1808.4-3658. By directly fitting the pulse profile shapes with our model, we find that the main parameters of the emitting spot such as its mean latitude and longitude as well as the emissivity pattern change irregularly causing small shifts in pulse phase, and the strong profile variations are caused by the increasing inner disk radius. We finally notice that significant variations in the pulse profiles in the 2002 and 2008 outbursts of SAX J1808.4-3658 happen at fluxes differing by a factor of 2, which can be explained if the inner disk radius is not a simple function of the accretion rate, but depends on the previous history.

preprint2009arXiv

Polarization properties of X-ray millisecond pulsars

Radiation of X-ray bursts and of accretion shocks in weakly magnetized neutron stars in low-mass X-ray binaries is produced in plane-parallel atmospheres dominated by electron scattering. We first discuss polarization produced by single (non-magnetic) Compton scattering, in particular the depolarizing effect of high electron temperature, and then the polarization due to multiply electron scattering in a slab. We further predict the X-ray pulse profiles and polarization properties of nuclear- and accretion-powered millisecond pulsars. We introduce a relativistic rotation vector model, which includes the effect of rotation of polarization plane due to the rapid motion of the hot spot as well as the light bending. Future observations of the X-ray polarization will provide a valuable tool to test the geometry of the emission region in pulsars and its physical characteristics.

preprint2009arXiv

Spectral variability of ultraluminous X-ray sources

We study spectral variability of 11 ultraluminous X-ray sources (ULX) using archived XMM-Newton and Chandra observations. We use three models to describe the observed spectra: a power-law, a multi-colour disc (MCD) and a combination of these two models. We find that 7 ULXs show a correlation between the luminosity Lx and the photon index Gamma. Furthermore, 4 out of these 7 ULXs also show spectral pivoting in the observed energy band. We also find that two ULXs show an Lx-Gamma anti-correlation. The spectra of 4 ULXs in the sample can be adequately fitted with a MCD model. We compare these sources to known black hole binaries (BHB) and find that they follow similar paths in their luminosity-temperature diagrams. Finally we show that the `soft excess' reported for many of these ULXs at about 0.2 keV seems to roughly follow a trend Lsoft \propto T^{-3.5} when modelled with a power-law plus a `cool' MCD model. This is contrary to the L \propto T^4 relation that is expected from theory and what is seen for many accreting BHBs. The observed trend could instead arise from disc emission beamed by an outflowing wind around a about 10 solar mass black hole.

preprint2009arXiv

Time-dependent modeling of radiative processes in hot magnetized plasmas

Numerical simulations of radiative processes in magnetized compact sources such as hot accretion disks around black holes, relativistic jets in active galaxies and gamma-ray bursts are complicated because the particle and photon distributions span many orders of magnitude in energy, they also strongly depend on each other, the radiative processes behave significantly differently depending on the energy regime, and finally due to the enormous difference in the time-scales of the processes. We have developed a novel computer code for the time-dependent simulations that overcomes these problems. The processes taken into account are Compton scattering, electron-positron pair production and annihilation, Coulomb scattering as well as synchrotron emission and absorption. No approximation has been made on the corresponding rates. For the first time, we solve coupled integro-differential kinetic equations for photons and electrons/positrons without any limitations on the photon and lepton energies. A numerical scheme is proposed to guarantee energy conservation when dealing with synchrotron processes in electron and photon equations. We apply the code to model non-thermal pair cascades in the blackbody radiation field, to study the synchrotron self-absorption as particle thermalization mechanism, and to simulate time evolution of stochastically heated pairs and corresponding synchrotron self-Compton photon spectra which might be responsible for the prompt emission of gamma-ray bursts. Good agreement with previous works is found in the parameter regimes where comparison is feasible, with the differences attributable to our improved treatment of the microphysics.

preprint2008arXiv

Gamma-ray emission of relativistic jets as a supercritical process

Supercriticality of the same kind as that in a nuclear pile can take place in high-energy astrophysical objects producing a number of impressive effects. For example, it could cause an explosive release of the energy of a cloud of ultrarelativistic protons into radiation. More certainly, supercriticality should be responsible for energy dissipation of very energetic relativistic fluids such as ultrarelativistic shocks in gamma-ray bursts and jets in active galactic nuclei (AGNs). In this case, the photon breeding process operates. It is a kind of the converter mechanism with the high-energy photons and e^+ e^- pairs converting into each other via pair production and inverse Compton scattering. Under certain conditions, which should be satisfied in powerful AGNs, the photon breeding mechanism becomes supercritical: the high-energy photons breed exponentially until their feedback on the fluid changes its velocity pattern. Then the system comes to a self-adjusting near-critical steady state. Monte-Carlo simulations with the detailed treatment of particle propagation and interactions demonstrate that a jet with the Lorentz factor Gamma ~ 20 can radiate away up to a half of its total energy and for Gamma=40 the radiation efficiency can be up to 80 per cent. Outer layers of the jet decelerate down to a moderate Lorentz factor 2-4, while the spine of the jet has the final Lorentz factor in the range 10-20 independently on the initial Gamma. Such sharp deceleration under the impact of radiation must cause a number of interesting phenomena such as formation of internal shocks and an early generation of turbulence.

preprint2008arXiv

Photon breeding mechanism in relativistic jets: astrophysical implications

Photon breeding in relativistic jets involves multiplication of high-energy photons propagating from the jet to the external environment and back with the conversion into electron-positron pairs. The exponential growth of the energy density of these photons is a super-critical process powered by the bulk energy of the jet. The efficient deceleration of the jet outer layers creates a structured jet morphology with the fast spine and slow sheath. In initially fast and high-power jets even the spine can be decelerated efficiently leading to very high radiative efficiencies of conversion of the jet bulk energy into radiation. The decelerating, structured jets have angular distribution of radiation significantly broader than that predicted by a simple blob model with a constant Lorentz factor. This reconciles the discrepancy between the high Doppler factors determined by the fits to the spectra of TeV blazars and the low apparent velocities observed at VLBI scales as well as the low jet Lorentz factors required by the observed statistics and luminosity ratio of Fanaroff-Riley I radio galaxies and BL Lac objects. Photon breeding produces a population of high-energy leptons in agreement with the constraints on the electron injection function required by spectral fits of the TeV blazars. Relativistic pairs created outside the jet and emitting gamma-rays by inverse Compton process might explain the relatively high level of the TeV emission from the misaligned jet in the radio galaxies. The mechanism reproduces basic spectral features observed in blazars including the blazar sequence (shift of the spectral peaks towards lower energies with increasing luminosity). The mechanism is very robust and can operate in various environments characterised by the high photon density.

preprint2008arXiv

Spectral variability of ultraluminous X-ray sources

We study spectral variability of 11 ultraluminous X-ray sources (ULX) using archived XMM-Newton and Chandra observations. We use three models to describe the observed spectra; a power-law, a multi-colour disk (MCD) and a combination of these two models. We find that out of the 11 ULXs in our sample, 7 ULXs show a correlation between the luminosity and the photon index Gamma (hereafter L-Gamma correlation). Furthermore, out of the 7 ULXs that have the L-Gamma correlation, 4 ULXs also show spectral pivoting in the observed energy band. We also find that two ULXs show an L-Gamma anti-correlation. The spectra of 4 ULXs in the sample can be adequately fitted with a MCD model. We compare these sources to known black hole binaries (BHB) and find that they follow similar paths in their luminosity-temperature (hereafter L-T) diagrams. Finally we show that the 'soft excess' reported for many of these ULXs at 0.2 keV seem to follow a trend L \propto T^{-4} when modeled with a power-law plus a 'cool' MCD model. This is contrary to the expected L \propto T^4 relation that is expected from theory and what is seen for many accreting BHBs.

preprint2005arXiv

Gamma-ray burst spectra from continuously accelerated electrons

We discuss here constraints on the particle acceleration models from the observed gamma-ray bursts spectra. The standard synchrotron shock model assumes that some fraction of available energy is given instantaneously to the electrons which are injected at high Lorentz factor. The emitted spectrum in that case corresponds to the spectrum of cooling electrons, F_ν~ ν^{-1/2}, is much too soft to account for the majority of the observed spectral slopes. We show that continuous heating of electrons over the life-time of a source is needed to produce hard observed spectra. In this model, a prominent peak develops in the electron distribution at energy which is a strong function of Thomson optical depth τ_T of heated electrons (pairs). At τ_T>1, a typical electron Lorentz factor γ~ 1-2 and quasi-thermal Comptonization operates. It produces spectrum peaking at a too high energy. Optical depths below 10^{-4} would be difficult to imagine in any physical scenario. At τ_T =10^{-4}-10^{-2}, γ~ 30-100 and synchrotron self-Compton radiation is the main emission mechanism. The synchrotron peak should be observed at 10--100 eV, while the self-absorbed low-energy tail with F_ν~ ν^2 can produce the prompt optical emission (like in the case of GRB 990123). The first Compton scattering radiation by nearly monoenergetic electrons peaks in the BATSE energy band and can be as hard as F_ν~ ν^1 reproducing the hardness of most of the observed GRB spectra. The second Compton peak should be observed in the high-energy gamma-ray band, possibly being responsible for the 10-100 MeV emission detected in GRB 941017. A significant electron-positron pair production reduces the available energy per particle, moving spectral peaks to lower energies as the burst progresses.

preprint2002arXiv

Understanding the long-term spectral variability of Cygnus X-1 with BATSE and ASM observations

We present an analysis of all observations of Cygnus X-1 by the BATSE (20-300 keV) and by ASM (1.5-12 keV) until 2002 June, including 1200 days of simultaneous data. We find a number of correlations between fluxes and hardnesses in different energy bands. In the hard state, the variability can be explained by softening the overall spectrum with a pivot at 50 keV and another, independent variability pattern where the spectral shape does not change when the luminosity changes. In the soft state, the variability is caused by a variable hard tail of a constant shape superimposed on a constant soft component. These variability patterns are in agreement with the energy-dependent rms variability in the two states. We also study in detail recent soft states in 2000-02. The last of them has lasted so far for >200 days. Their spectra are harder in the 1.5-5 keV band but similar or in the 3-12 keV band than those of the 1996 soft state whereas the rms variability is stronger in all the ASM bands. On the other hand, the 1994 soft state transition observed by BATSE appears very similar to the 1996 one. We interpret the variability patterns by theoretical Comptonization models. In the hard state, the variability appears to be driven mostly by changing flux in seed photons Comptonized in a hot thermal plasma cloud with an approximately constant power supply. In the soft state, the variability is consistent with flares of hybrid, thermal/nonthermal, plasma with variable power above a stable cold disk. The spectral and timing differences between the 1996 and 2000-02 soft states are explained by a decrease of the color disk temperature. Also, based on broad-band pointed observations, we find the intrinsic bolometric luminosity increases by a factor of 3-4 from the hard state to the soft one.

preprint2000arXiv

Time Lags in Compact Objects: Constraints on the Emission Models

Accreting black holes and neutron stars in their hard (low) state show not only very similar X/gamma-ray spectra but also that the behaviour of their light curves is quite similar which can be quantified as having similar power-density spectra and Fourier-frequency-dependent time/phase lags. Taken together this argues for a common mechanism of the X/gamma-ray production in these objects. This mechanism is probably a property of the accretion flow only since it does not depend on the nature of the compact object. In this paper, I review the observational data paying most attention to the properties of the temporal variability such as the time/phase lags that hopefully can help us to discriminate between different theoretical models. I also discuss the models developed to account for the basic observational facts. Particularly, I show that the commonly used Compton cloud models with constant temperature cannot explain variable sources without violating the energy conservation law. Alternative models where time lags are related to the spectral evolution during X-ray flares are discussed and compared with observations. Compton reflection from the outer edge of the accretion disc is shown to markedly affect the time lag Fourier spectrum.

preprint1997arXiv

The Nature of Spectral Transitions in Accreting Black Holes: The Case of Cyg X-1

Accreting black holes radiate in one of several spectral states, switching from one to another for reasons that are as yet not understood. Using the best studied example, Cyg X-1, we identify the geometry and physical conditions characterizing these states. In particular, we show that in the hard state most of the accretion energy is dissipated in a corona-like structure which fills the inner few tens of gravitational radii around the black hole and has Compton optical depth of order unity. In this state, an optically thick accretion disc extends out to greater distance, but penetrates only a short way into the coronal region. In the soft state, the optically thick disc moves inward and receives the majority of the dissipated energy, while the "corona" becomes optically thin and extends around much of the inner disc. The mass accretion rate in both states is $~10^{-8} M_{\odot}$ yr$^{-1}$.

preprint1997arXiv

The X/gamma-ray spectral properties of NGC 7172

We present a combined, non-simultaneous, ASCA GIS and CGRO OSSE spectrum of the Seyfert 2 galaxy, NGC 7172, and make broad band spectral fits. The only Seyfert 2 galaxy previously studied over such a broad band is NGC 4945. We find that the most probable model for the data is a power law with an exponential cut-off being affected by a neutral absorber. The best fit parameters are found to be $Γ= 1.47 \pm 0.15$ and $\NH = (7.8 \pm 0.6) 10^{22}$ cm$^{-2}$. The spectral index of the underlying power law of NGC 7172 has therefore varied from 1.8 to 1.5 since the Ginga observations in 1989. For this simple model the e-folding energy at $88^{+65}_{-28}$ keV is relatively well constrained. The observed flux in the 2-10 keV range is $F_{2-10}=4.7 10^{-11}$ erg cm$^{-2}$ s$^{-1}$, which corresponds to a small increase since the Ginga measurement in October 1989.

preprint1997arXiv

Thermal/non-thermal model of Cyg X-1 in the soft state

We analyze X-ray/gamma-ray observations of Cyg X-1 in the soft state. We find the hybrid thermal/non-thermal model fits the data very well. The significant contribution from the non-thermal electrons is required to account for the observed power law extending without a break to at least 800 keV. The presence of electron-positron pairs in the hot source is unlikely. We propose the geometry in which a cold accretion disk, surrounded by a hot corona, extends down to the last stable orbit. At the observed accretion rate the corona ensures the stability of the cold disk.

preprint1995arXiv

ON THE GEOMETRY OF THE X-RAY EMITTING REGION IN SEYFERT GALAXIES

For the first time, detailed radiative transfer calculations of Comptonized X-ray and gamma-ray radiation in a hot pair plasma above a cold accretion disk are performed using two independent codes and methods. The simulations include both energy and pair balance as well as reprocessing of the X- and gamma-rays by the cold disk. We study both plane-parallel coronae as well as active dissipation regions having shapes of hemispheres and pill boxes located on the disk surface. It is shown, contrary to earlier claims, that plane-parallel coronae in pair balance have difficulties in selfconsistently reproducing the ranges of 2-20 keV spectral slopes, high energy cutoffs, and compactnesses inferred from observations of type 1 Seyfert galaxies. Instead, the observations are consistent with the X-rays coming from a number of individual active regions located on the surface of the disk. A number of effects such as anisotropic Compton scattering, the reflection hump, feedback to the soft photon source by reprocessing, and an active region in pair equilibrium all conspire to produce the observed ranges of X-ray slopes, high energy cutoffs, and compactnesses. The spread in spectral X-ray slopes can be due to a spread in the properties of the active regions such as their compactnesses and their elevations above the disk surface. Simplified models invoking isotropic Comptonization in spherical clouds are no longer sufficient when interpreting the data.

Juri Poutanen

What is connected

Connect this record

See the researcher in context

Building this map preview

77 published item(s)

Discovery of High X-Ray Polarization from the Neutron Star Low-Mass X-Ray Binary Cyg X-2 in the Horizontal Branch

Determination of X-ray pulsar geometry with IXPE polarimetry

Black hole spin-orbit misalignment in the X-ray binary MAXI J1820+070

Constraints on the magnetic field structure in accreting compact objects from aperiodic variability

Limits on X-ray Polarization at the Core of Centaurus A as Observed with the Imaging X-ray Polarimetry Explorer

Optical polarization signatures of black hole X-ray binaries

Simultaneous space and phase resolved X-ray polarimetry of the Crab Pulsar and Nebula

Analytical techniques for polarimetric imaging of accretion flows in Schwarzschild metric

Broad-band analysis of X-ray pulsar 2S 1845$-$024

Mechanical model of a boundary layer for the parallel tracks of kilohertz quasi-periodic oscillations in accreting neutron stars

Neutron star parameter constraints for accretion-powered millisecond pulsars from the simulated IXPE data

Accurate analytic formula for light bending in Schwarzschild metric

Colors and patterns of black hole X-ray binary GX 339-4

Disc and wind in black hole X-ray binary MAXIJ1820+070 observed through polarized light during its 2018 outburst

Discovery of a retrogradely rotating neutron star in the X-ray pulsar GX 301$-$2

Kilohertz quasi-periodic oscillations from neutron star spreading layers

Magnetospheric return-current-heated atmospheres of rotation-powered millisecond pulsars

Observational appearance of rapidly rotating neutron stars: X-ray bursts, cooling tail method, and radius determination

Spectral and temporal properties of Compton scattering by mildly relativistic thermal electrons

The unusual behaviour of the young X-ray pulsar SXP 1062 during the 2019 outburst

Compton scattering S-matrix and cross section in strong magnetic field

Gamma-ray opacity of the anisotropic stratified broad-line regions in blazars

Measuring the neutron star equation of state using X-ray timing

Propeller effect in action in the ultraluminous accreting magnetar M82 X-2

Discovery of correlated optical/X-ray quasi-periodic oscillations in black hole binary SWIFT J1753.5-0127

Models of neutron star atmospheres enriched with nuclear burning ashes

NuSTAR discovery of a cyclotron absorption line in the transient X-ray pulsar 2S 1553-542

On the maximum accretion luminosity of magnetized neutron stars: connecting X-ray pulsars and ultraluminous X-ray sources

Polarization modulation from Lense-Thirring precession in X-ray binaries

Positive correlation between the cyclotron line energy and luminosity in sub-critical X-ray pulsars: Doppler effect in the accretion channel

Reprocessing model for the optical quasi-periodic oscillations in black hole binaries

Colours of black holes: infrared flares from the hot accretion disc in XTE J1550-564

Spectroscopic evidence for a low-mass black hole in SWIFT J1753.5-0127

The critical accretion luminosity for magnetized neutron stars

The effect of accretion on the measurement of neutron star mass and radius in the low-mass X-ray binary 4U 1608-52

The influence of accretion geometry on the spectral evolution during thermonuclear (type-I) X-ray bursts

The mystery of spectral breaks: Lyman continuum absorption by photon-photon pair production in the Fermi GeV spectra of bright blazars

A reflection model for the cyclotron lines in the spectra of X-ray pulsars

A unified Lense-Thirring precession model for optical and X-ray quasi-periodic oscillations in black hole binaries

Hot accretion flow in black hole binaries: a link connecting X-rays to the infrared

Modelling spectral and timing properties of accreting black holes: the hybrid hot flow paradigm

On the association of the ultraluminous X-ray sources in the Antennae galaxies with young stellar clusters

On the origin of cyclotron lines in the spectra of X-ray pulsars

On the spreading layer emission in luminous accreting neutron stars

White Paper for Blazar Observations with a GEMS-like X-ray Polarimetry Mission

White Paper on GEMS Study of Polarized X-rays from Neutron Stars

X-ray Polarization from Black Holes: GEMS Scientific White Paper

XIPE: the X-ray Imaging Polarimetry Explorer

Young rotation-powered pulsars as ultra-luminous X-ray sources

Joint spectral-timing modelling of the hard lags in GX 339-4: constraints on reflection models

Relativistic kinetic equation for Compton scattering of polarized radiation in strong magnetic field

A Complexity-Brightness Correlation in Gamma Ray Bursts

A self-consistent hybrid Comptonization model for broad-band spectra of accreting supermassive black holes

A synchrotron self-Compton -- disk reprocessing model for optical/X-ray correlation in black hole X-ray binaries

Fermi Observations of Blazars: Implications for Gamma-ray Production

Gamma-ray bursts from magnetized collisionally-heated jets

Neutron star stiff equation of state derived from cooling phases of the X-ray burster 4U 1724-307

Variation of the gamma-gamma opacity by the He II Lyman continuum constrains the location of the gamma-ray emission region in the blazar 3C 454.3

Varying disc-magnetosphere coupling as the origin of pulse profile variability in SAX J1808.4-3658

Constraining compactness and magnetic field geometry of X-ray pulsars from the statistics of their pulse profiles

GeV breaks in blazars as a result of gamma-ray absorption within the broad-line region

Theory of Compton scattering by anisotropic electrons

Accreting millisecond pulsar SAX J1808.4-3658 during its 2002 outburst: evidence for a receding disc

On the nature of pulse profile variations and timing noise in accreting millisecond pulsars

Polarization properties of X-ray millisecond pulsars

Spectral variability of ultraluminous X-ray sources

Time-dependent modeling of radiative processes in hot magnetized plasmas

Gamma-ray emission of relativistic jets as a supercritical process

Photon breeding mechanism in relativistic jets: astrophysical implications

Spectral variability of ultraluminous X-ray sources

Gamma-ray burst spectra from continuously accelerated electrons

Understanding the long-term spectral variability of Cygnus X-1 with BATSE and ASM observations

Time Lags in Compact Objects: Constraints on the Emission Models

The Nature of Spectral Transitions in Accreting Black Holes: The Case of Cyg X-1