Source author record

Andrzej A. Zdziarski

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

50works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

50 published item(s)

preprint2024arXiv

Black hole spin measurements in LMC X-1 and Cyg X-1 are highly model-dependent

The black-hole spin parameter, $a_*$, was measured to be close to its maximum value of 1 in many accreting X-ray binaries. In particular, $a_*\gtrsim 0.9$ was found in a number of studies of LMC X-1. These measurements were claimed to take into account both statistical and systematic uncertainties. We perform new measurements using a recent simultaneous observation by NICER and NuSTAR, providing a data set of high quality. We use the disk continuum method together with improved models for coronal Comptonization. With the standard relativistic disk model and optically thin Comptonization, we obtain values of $a_*$ similar to those obtained before. We then consider modifications to the standard model. Using a color correction of 2, we find $a_*\approx 0.64$--0.84. We then consider disks with dissipation in surface layers. To account for that, we assume the standard disk is covered by a warm and optically thick Comptonizing layer. Our model with the lowest $χ^2$ yields then $a_*\approx 0.40^{+0.41}_{-0.32}$. In order to test the presence of such effects in other sources, we also study an X-ray observation of Cyg X-1 by Suzaku in the soft state. We confirm the previous findings of $a_*>0.99$ using the standard model, but then we find a weakly constrained $a_*\approx 0.82^{+0.16}_{-0.74}$ when including an optically thick Comptonizing layer. We conclude that determinations of the spin using the continuum method can be highly sensitive to the assumptions about the disk structure.

preprint2022arXiv

A simple analytical model of magnetic jets

We propose a simple analytical jet model of magnetic jets, in which radially-averaged profiles of main physical quantities are obtained based on conservation laws and some results of published GRMHD jet simulations. We take into account conversion of the magnetic energy flux to bulk acceleration in jets formed around rotating black holes assuming the mass continuity equation and constant jet power, which leads to the Bernoulli equation. For assumed profiles of the bulk Lorentz factor and the radius, this gives us the profile of the toroidal magnetic field component along the jet. We then consider the case where the poloidal field component is connected to a rotating black hole surrounded by an accretion disc. Our formalism then recovers the standard formula for the power extracted from a rotating black hole. We find that the poloidal field strength dominates over the toroidal one in the comoving frame up to large distances, which means that jets should be more stable to current-driven kink modes. The resulting magnetic field profiles can then be used to calculate the jet synchrotron emission.

preprint2022arXiv

Improved model of X-ray emission from hot accretion flows

We have developed an improved model of X-ray emission from optically thin, two-temperature accretion flows, \texttt{kerrflow}, using an exact Monte Carlo treatment of global Comptonization as well as with a fully general relativistic description of both the radiative and hydrodynamic processes. It also includes pion-decay electrons, whose synchrotron emission dominates the seed photons yield at high accretion rates in flows around supermassive black holes. We consider in detail the dependence of the model spectra on the black hole spin, the electron heating efficiency, the plasma magnetization and the accretion rate, and we discuss feasibility of constraining these parameters by analyzing X-ray spectra of nearby low-luminosity active galactic nuclei. We note some degeneracies which hinder precise estimations of these parameters when individual X-ray spectra are analyzed. These degeneracies are eliminated when several spectra from a given source are fitted jointly, which then allows us to reliably measure the model parameters. We find significant differences with previous spectral models of hot-flow emission, related with the computational methods for Comptonization. Finally, we briefly consider and discuss the dependence on the viscosity parameter and on the outflow strength.

preprint2022arXiv

Insight-HXMT, NuSTAR and INTEGRAL Data Show Disk Truncation in the Hard State of the Black-Hole X-Ray Binary MAXI J1820+070

We study X-ray and soft gamma-ray spectra from the hard state of the accreting black-hole binary MAXI J1820+070. We perform analysis of joint spectra from HXMT, NuSTAR and INTEGRAL. We find an overall agreement between the spectra from all three satellites. Satisfactory fits to the data require substantial spectral complexity, with our models including two Comptonization regions and their associated disk reflection, a disk blackbody and a narrow Fe K$α$ line. Our fits confirm the presence of the truncation of the reflecting optically-thick disk at least at $>$10 gravitational radii. However, we find that the HXMT data alone cannot significantly constrain the disk inner radii.

preprint2022arXiv

The Composition and Power of the Jet of the Broad-Line Radio Galaxy 3C 120

We calculated the electron-positron pair production rate at the base of the jet of 3C 120 by collisions of photons from the hot accretion flow using the measurement of its average soft gamma-ray spectrum by Compton Gamma Ray Observatory. We found that this rate approximately equals the flow rate of the leptons emitting the observed synchrotron radio-to-IR spectrum of the jet core, calculated using the extended jet model following Blandford \& Konigl. This coincidence shows the jet composition is likely to be pair-dominated. We then calculated the jet power in the bulk motion of ions, and found it greatly exceeds that achievable by the magnetically arrested disk scenario for the maximum black hole spin unless the jet contains mostly pairs. Next, we found that the magnetic flux through the synchrotron-emitting jet equals the maximum poloidal flux that can thread the black hole. Finally, we compared two estimates of the magnetization parameter at the onset of the synchrotron emission and found they are in agreement only if pairs dominate the jet content.

preprint2022arXiv

The Donor of the Black-Hole X-Ray Binary MAXI J1820+070

We estimate the parameters of the donor of the accreting black-hole binary MAXI J1820+070. The measured values of the binary period, rotational and radial velocities and constraints on the orbital inclination imply the donor is a subgiant with the mass of $M_2\approx 0.49^{+0.10}_{-0.10}M_\odot$ and the radius of $R_2\approx 1.19^{+0.08}_{-0.08}R_\odot$. We re-analyze the previously obtained optical spectrum from the Gran Telescopio Canarias, and found it yields a strict lower limit on the effective temperature of $T>4200$ K. We compile optical and infrared fluxes observed during the quiescence of this system. From the minima $r$ and $i$-band fluxes found in Pan-STARSS1 Data Release 2 pre-discovery imaging and for a distance of $D\approx3$ kpc, reddening of $E(B$--$V)=0.23$ and $R_2\approx{1.11R_\odot}$, we find $T\lesssim4230$ K, very close to the above lower limit. For a larger distance, the temperature can be higher, up to about 4500 K (corresponding to a K5 spectral type, preferred by previous studies) at $D=3.5$ kpc, allowed by the Gaia parallax. We perform evolutionary calculations for the binary system and compare them to the observational constraints. Our model fitting the above temperature and radius constraints at $D\approx 3$ kpc has the mass of $0.4M_\odot$, $T\approx4200$ K and solar metallicity. Two alternative models require $D\gtrsim 3.3$--3.4 kpc at $0.4 M_\odot$, $T\approx4500$ K and half solar metallicity, and $0.5M_\odot$, $T\approx4300$ K and solar metallicity. These models yield mass transfer rates of $\sim\!\!10^{-10}M_\odot$/yr, compatible with those based on the estimated accreted mass of $\approx\!2\times 10^{25}$ g and the time between the 2018 discovery and the 1934 historical outburst.

preprint2022arXiv

What are the Composition and Power of the Jet in Cyg X-1?

We calculate the electron-positron pair production rate at the base of the jet of Cyg X-1 by collisions of photons from its hot accretion flow using the measurement of its average soft gamma-ray spectra by the Compton Gamma Ray Observatory and INTEGRAL satellites. We have found that this rate approximately equals the flow rate of the leptons emitting the observed synchrotron radio-to-IR spectrum of the jet, calculated using an extended jet model following that of Blandford \& Koenigl. This coincidence shows the jet composition is likely to be pair-dominated. The same coincidences were found before in the microquasar MAXI J1820+070 and in the radio galaxy 3C 120, which shows that the considered mechanism can be universal for at least some classes of relativistic jets. Furthermore, we recalculate the jet power of Cyg X-1. The presence of pairs can strongly reduces the power in the bulk motion of ions, which then limits the parameter space at which the jet can power the $\sim$5-pc nebular structure present in its vicinity.

preprint2021arXiv

Accretion Geometry in the Hard State of the Black-Hole X-Ray Binary MAXI J1820+070

We study X-ray spectra from the outburst rise of the accreting black-hole binary MAXI J1820+070. We find that models having the disk inclinations within those of either the binary or the jet imply significant changes of the accretion disk inner radius during the luminous part of the hard spectral state, with that radius changing from $>$100 to $\sim$10 gravitational radii. The main trend is a decrease with the decreasing spectral hardness. Our analysis requires the accretion flow to be structured, with at least two components with different spectral slopes. The harder component dominates the bolometric luminosity and produces strong, narrow, X-ray reflection features. The softer component is responsible for the underlying broader reflection features. The data are compatible with the harder component having a large scale height, located downstream the disk truncation radius, and reflecting mostly from remote parts of the disk. The softer component forms a corona above the disk up to some transition radius. Our findings can explain the changes of the characteristic variability time scales, found in other works, as being driven by the changes of the disk characteristic radii.

preprint2021arXiv

Jet Parameters in the Black-Hole X-Ray Binary MAXI J1820+070

We study the jet in the hard state of the accreting black-hole binary MAXI J1820+070. From the available radio-to-optical spectral and variability data, we put strong constraints on the jet parameters. We find while it is not possible to uniquely determine the jet Lorentz factor from the spectral and variability properties alone, we can estimate the jet opening angle ($\approx 1.5\pm1^\circ$), the distance at which the jet starts emitting synchrotron radiation ($\sim 3 \times10^{10}$cm), and the magnetic field strength there ($\sim$10$^4$G), with relatively low uncertainty, as they depend weakly on the bulk Lorentz factor. We find the breaks in the variability power spectra from radio to sub-mm are consistent with variability damping over the time scale equal to the travel time along the jet at any Lorentz factor. This factor can still be constrained by the electron-positron pair production rate within the jet base, which we calculate based on the observed X-ray/soft gamma-ray spectrum, and the jet power, required to be less than the accretion power. The minimum ($\sim$1.5) and maximum ($\sim$4.5) Lorentz factors correspond to the dominance of pairs and ions, and the minimum and maximum jet power, respectively. We estimate the magnetic flux threading the black hole and find the jet can be powered by the Blandford-Znajek mechanism in a magnetically-arrested flow accretion flow. We point out the similarity of our derived formalism to that of core shifts, observed in extragalactic radio sources.

preprint2020arXiv

Correlations and lags between X-ray and radio emission of Cyg X-1

Cygnus X-1 is a Galactic black-hole binary and a microquasar. Its X-ray emission originates most likely in the accretion flow, while the radio emission is from a magnetized jet. We study long-term cross-correlations between its X-ray and radio fluxes, searching for lags of the peak of the correlation. We use soft X-ray data from the All Sky Monitor on board of Rossi X-ray Timing Explorer, hard X-ray data from the Burst Alert Telescope on board of Neil Gehrels Swift, and 15 GHz radio monitoring data from the Ryle and AMI telescopes, which light curves cover all of the spectral states of Cyg X-1. Previous studies of radio/X-ray correlations in black-hole binaries concentrated on the soft X-ray/radio correlation in the hard spectral state. In contrast, our strongest correlation is between the hard X-rays and radio in the entire data at a $\lesssim 1$-day lag. This indicates the jet formation from a hard X-ray emitting corona/hot accretion flow through all states of Cyg X-1. We also find evidence for a long lag of the radio with respect to the soft X-rays in the soft state. This indicates long-term advection of the magnetic field through the accretion disc from its high-mass donor, which is similar to the case of another microquasar with a high-mass donor, Cyg X-3.

preprint2020arXiv

Distinct accretion modes of Cygnus X-1 revealed from hard X-rays

Thanks to recurrent observations of the black hole binary Cyg X-1 carried out over 15 years the INTEGRAL satellite has collected the largest data set in the hard X-ray band for this source. We have analyzed these data, complemented by data collected by other X-ray satellites and radio flux at 15 GHz. To characterize the spectral and variability properties of the system we have examined parameters such as the hard X-ray flux, photon index and fractional variability. Our main result is that the 2D distribution of the photon index and flux determined for the 22-100 keV band forms six clusters. This result, interpreted within the Comptonization scenario as the dominant process responsible for the hard X-ray emission, leads to a conclusion that the hot plasma in Cyg X-1 takes the form of six specific geometries. The distinct character of each of these plasma states is reinforced by their different X-ray and radio variability patterns. In particular, the hardest and softest plasma states show no short-term flux - photon index correlation typical for the four other states, implying a lack of interaction between the plasma and accretion disk. The system evolves between these two extreme states, with the spectral slope regulated by a variable cooling of the plasma by the disk photons.

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 persistent radio jet coupled to hard X-rays in the soft state of Cyg X-1

We study long-term radio/X-ray correlations in Cyg X-1. We find the persistent existence of a compact radio jet in its soft state. This represents a new phenomenon in black-hole binaries, in addition to compact jets in the hard state and episodic ejections of ballistic blobs in the intermediate state. While the radio emission in the hard state is strongly correlated with both the soft and hard X-rays, the radio flux in the soft state is not directly correlated with the flux of the dominant disk blackbody in soft X-rays, but instead it is lagged by about a hundred days. We interpret the lag as occurring in the process of advection of the magnetic flux from the donor through the accretion disk. On the other hand, the soft-state radio flux is very tightly correlated with the hard X-ray, 15--50 keV, flux without a measurable lag and at the same rms. This implies that the X-ray emitting disk corona and the soft-state jet are powered by the same process, probably magnetically.

preprint2020arXiv

Two major constraints on the inner radii of accretion disks

The Stefan-Boltzmann law yields a fundamental constraint on the geometry of inner accretion disks in black-hole X-ray binaries. It follows from considering the irradiating flux and the effective temperature of the inner parts of the disk, which implies that a strong quasi-thermal component with the average energy higher than that of a blackbody at the effective temperature has to be present whenever relativistic Fe K fluorescence and reflection features are observed. The apparent absence of such quasi-thermal component with the color temperature of $\sim$1 keV in high-luminosity hard states is not compatible with a strongly irradiated disk extending close to the innermost stable circular orbit. Instead, the disk should be either truncated at a relatively large radius or irradiated by a corona at a large height, which would reduce the effective temperature and bring it to an agreement with the data. We also study constraints on disk/corona models following from comparing the disk densities fitted in literature using variable-density reflection codes with those calculated by us from the ionization parameter, the luminosity and the disk inner radius. We find that the fitted densities are much higher/lower in the hard/soft state of binaries, implying significant problems with the used assumptions and methods.

preprint2019arXiv

Jets in the soft state in Cyg X-3 caused by advection of the donor magnetic field and unification with low-mass X-ray binaries

The high-mass accreting binary Cyg X-3 is distinctly different from low-mass X-ray binaries (LMXBs) in having powerful radio and $γ$-ray emitting jets in its soft spectral state. However, the transition from the hard state to the soft one is first associated with quenching of the hard-state radio emission, as in LMXBs. The powerful soft-state jets in Cyg X-3 form, on average, $\sim$50 d later. We interpret the initial jet quenching as due to the hard-state vertical magnetic field quickly diffusing away in the thin disc extending to the innermost stable circular orbit in the soft state, or, if that field is produced in situ, also cessation of its generation. The subsequent formation of the powerful jets occurs due to advection of the magnetic field from the donor. We find this happens only above certain threshold accretion rate associated with appearance of magnetically driven outflows. The $\sim$50 d lag is of the order of the viscous timescale in the outer disc, while the field advection is much faster. This process does not happen in LMXBs due to the magnetic fluxes available from their donors being lower than that for the wind accretion from the Wolf-Rayet donor of Cyg X-3. In our model, the vertical magnetic field in the hard state, required to form the jets both in Cyg X-3 and LMXBs, is formed in situ rather than advected from the donor. Our results provide a unified scenario of the soft and hard states in both Cyg X-3 and LMXBs.

preprint2019arXiv

Radiative Properties of Magnetically-Arrested Disks

Magnetically-arrested disks (MADs) appear when accretion flows are supplied with a sufficient amount of magnetic flux. In this work, we use results of magnetohydrodynamic simulations to set the configuration of the magnetic field and investigate the dynamics and radiative properties of the resulting accretion flow (i.e., without that of the jet) of MAD. The method developed here is applied to both the MAD and the standard and normal evolution (SANE) accretion flow with or without large scale magnetic fields. For the radiative processes, we include synchrotron, bremsstrahlung and Compton scattering. We find that, in general, MAD accretion flows have similar spectra to those of the SANE, which complicates the task of distinguishing MADs from SANEs. At the same accretion rates, MADs are systematically brighter than SANEs. However, the critical accretion rate above which the hot solution ceases to exist is lower in MAD. Consequently, the maximum luminosity the MAD can reach is comparable but slightly lower than that of SANE, and the dependence on the magnetic flux is weak. We then discuss implications of our results for active galactic nuclei and accreting black-hole binaries.

preprint2016arXiv

A luminous hot accretion flow in the low-luminosity active galactic nucleus NGC 7213

The active galactic nucleus (AGN) NGC 7213 shows a complex correlation between the monochromatic radio luminosity $L_R$ and the 2--10 keV X-ray luminosity $L_X$, i.e. the correlation is unusually weak with $p\sim 0$ (in the form $L_R\propto L_X^p$) when $L_X$ is below a critical luminosity, and steep with $p>1$ when $L_X$ is above that luminosity. Such a hybrid correlation in individual AGNs is unexpected as it deviates from the Fundamental Plane of AGN activity. Interestingly, a similar correlation pattern is observed in the black hole X-ray binary H1743--322, where it has been modelled by switching between different modes of accretion. We propose that the flat $L_R$--$L_X$ correlation of NGC 7213 is due to the presence of a luminous hot accretion flow, an accretion model whose radiative efficiency is sensitive to the accretion rate. Given the low luminosity of the source, $L_X\sim 10^{-4}$ of the Eddington luminosity, the viscosity parameter is determined to be small, $α\approx 0.01$. We also modelled the broad-band spectrum from radio to $γ$-rays, the time lag between the radio and X-ray light curves, and the implied size and the Lorentz factor of the radio jet. We predict that NGC 7213 will enter into a two-phase accretion regime when $L_X > 1.5 \times 10^{42}\, {\rm erg\,s^{-1}}$. When this happens, we predict a softening of the X-ray spectrum with the increasing flux and a steep radio/X-ray correlation.

preprint2016arXiv

Anisotropy of partially self-absorbed jets and the jet of Cyg X-1

We study the angular dependence of the flux from partially synchrotron self-absorbed conical jets (proposed by Blandford \& K{ö}nigl). We consider the jet viewed from either a side or close to on axis, and in the latter case, either from the jet top or bottom. We derive analytical formulae for the flux in each of these cases, and find the exact solution for an arbitrary angle numerically. We find that the maximum of the emission occurs when the jet is viewed from top on-axis, which is contrast to a previous result, which found the maximum at some intermediate angle and null emission on-axis. We then calculate the ratio of the jet-to-counterjet emission for this model, which depends on the viewing angle and the index of power-law electrons. We apply our results to the black-hole binary Cyg X-1. Given the jet-to-counterjet flux ratio of $\gtrsim$50 found observationally and the current estimates of the inclination, we find the jet velocity to be $\gtrsim 0.8c$. We also point out that when the projection effect is taken into account, the radio observations imply the jet half-opening angle of $\lesssim 1^{\rm o}$, a half of the value given before. When combined with the existing estimates of $Γ_{\rm j}$, the jet half-opening angle is low, $\ll 1/Γ_{\rm j}$, and much lower than values observed in blazars, unless $Γ_{\rm j}$ is much higher than currently estimated.

preprint2016arXiv

IGR J17451--3022: constraints on the nature of the donor star

We constrain the binary parameters of the eclipsing accreting X-ray binary IGR J17451--3022 and the nature of its donor star. The donor mass, its radius, and the system inclination angle are computed based on the system orbital period and eclipse duration recently reported by Bozzo et al.\ (2016). We find that the donor is most likely a main-sequence star with the mass comprised within the range $\sim$(0.5--$0.8){\rm M}_{\odot}$ and the radius of $\sim\! 0.7{\rm R}_{\odot}$. Assuming that the accreting compact object in IGR J17451--3022 is a neutron star, the duration of the nearly total rectangular eclipses yields the inclination angle of the system of $71^\circ\lesssim i\lesssim 76^\circ$, compatible with the presence of dips in this system. We rule out the presence of either brown or white dwarf. However, we find an alternative possibility that the donor star in IGR J17451--3022 could be a partially stripped giant with a very low mass, $\sim\! 0.2{\rm M}_{\odot}$. This case requires a substantial mass loss prior to the formation of the giant-star He core. According to that solution, the radius would be $\sim\! 0.4{\rm R}_{\odot}$, at $i\sim 80^\circ$. We additionally show that the well-known approximate dependence of the giant-star radius exclusively on its core mass breaks down below $\sim\!0.3{\rm M}_{\odot}$.

preprint2016arXiv

On the lamppost model of accreting black holes

We study the lamppost model, in which the X-ray source in accreting black-hole systems is located on the rotation axis close to the horizon. We point out a number of inconsistencies in the widely used lamppost model relxilllp, e.g., the neglect of the redshift of the photons emitted by the lamppost and directly observed. They appear to invalidate those model fitting results for which the source distances from the horizon are within several gravitational radii. Furthermore, if those results were correct, most of the photons produced in the lamppost would be trapped by the black hole, and the luminosity generated in the source as measured at infinity would be much larger than that observed. This appears to be in conflict with the observed smooth state transitions between the hard and soft states of X-ray binaries. The required increase of the accretion rate and the associated efficiency reduction present also a problem for AGNs. Then, those models imply the luminosity measured in the local frame to be much higher than that produced in the source and measured at infinity, due to the additional effects of time dilation and redshift, and the electron temperature to be higher significantly than that observed. We show that these conditions imply that the fitted sources would be out of the pair equilibrium. On the other hand, the above issues pose relatively minor problems for sources at large distances from the black hole, where relxilllp can still be used.

preprint2016arXiv

Spectral analysis of the XMM-Newton data of GX 339-4 in the low/hard state: disc truncation and reflection

We analyse all available observations of GX 339--4 by XMM-Newton in the hard spectral state. We jointly fit the spectral data by Comptonization and the currently best reflection code, relxill. We consider in detail a contribution from a standard blackbody accretion disc, testing whether its inner radius can be set equal to that of the reflector. However, this leads to an unphysical behaviour of the disc truncation radius, implying the soft X-ray component is not a standard blackbody disc. This appears to be due to irradiation by the hard X-rays, which strongly dominate the total emission. We consider a large array of models, testing, e.g., the effects of the chosen energy range, of adding unblurred reflection, and assuming a lamppost geometry. We find the effects of relativistic broadening to be relatively weak in all cases. In the coronal models, we find the inner radius to be large. In the lamppost model, the inner radius is unconstrained, but when fixed to the innermost stable orbit, the height of the source is large, which also implies a weak relativistic broadening. In the former models, the inner radius correlates with the X-ray hardness ratio, which is consistent with the presence of a truncated disc turning into a complete disc in the soft state. We also find the degree of the disc ionization to anti-correlate with the hardness, leading to strong spectral broadening due to scattering of reflected photons in the reflector in the softest studied states.

preprint2016arXiv

The power and production efficiency of blazar jets

We use published data on the power and production efficiency of jets in blazars with double radio lobes in order to compare results obtained using different methods. In order to eliminate selection effects, we use cross-matched sub-samples containing only luminous blazars. We compare the three main existing methods, namely those based on the emission of radio lobes, on spectral fitting, and on radio core shift. We find the average jet power obtained for identical samples with the radio-lobe method is $\sim10$ times lower than that from the spectral fitting. In turn, the power from spectral fitting is compatible with that from core-shift modelling for plausible parameters of the latter. We also consider a phenomenological estimator based on the γ-ray luminosity. We examine uncertainties of those methods and discuss two alternative hypotheses. In one, the blazar-fit and core-shift methods are assumed to be correct, and the lower power from radio lobes is caused by intermittency of accretion. Long periods of quiescence cause the energy in the radio lobes, accumulated over the lifetime of the blazar, to be much less than that estimated based on the present luminous state. In addition, the power calculated using the radio lobes can be underestimated for intrinsically compact jets, in which the radio core flux can be over-subtracted. In our second hypothesis, the radio-lobe method is assumed to be correct, and the blazar-fit and core-shift powers are reduced due to the presence of $\sim15$ pairs per proton, and a larger magnetization than usually assumed, respectively.

preprint2016arXiv

The radio/X-ray correlation in Cyg X-3 and the nature of its hard spectral state

We study the radio/X-ray correlation in Cyg X-3. It has been known that the soft and hard X-ray fluxes in the hard spectral state are correlated positively and negatively, respectively, with the radio flux. We show that this implies that the observed $\sim$1--100 keV flux (which is a fair approximation to the bolometric flux) is completely uncorrelated with the radio flux. We can recover a positive correlation (seen in other sources and expected theoretically) if the soft X-rays are strongly absorbed by a local medium. Then, however, the intrinsic X-ray spectrum of Cyg X-3 in its hard state becomes relatively soft, similar to that of an intermediate spectral state of black-hole binaries, but not to their true hard state. We also find the radio spectra in the hard state of Cyg X-3 are hard on average, and the flux distributions of the radio emission and soft X-rays can be described by sums of two log-normal functions. We compare Cyg X-3 with other X-ray binaries using colour-colour, colour-Eddington ratio and Eddington ratio-radio flux diagrams. We find Cyg X-3 to be spectrally most similar to GRS 1915+105, except that Cyg X-3 is substantially more radio loud, which appears to be due to its jet emission enhanced by interaction with the powerful stellar wind from the Wolf-Rayet donor.

preprint2015arXiv

Core shifts, magnetic fields and magnetization of extragalactic jets

We study the effect of radio-jet core shift, which is a dependence of the position of the jet radio core on the observational frequency. We derive a new method of measuring the jet magnetic field based on both the value of the shift and the observed radio flux, which complements the standard method that assumes equipartition. Using both methods, we re-analyse the blazar sample of Zamaninasab et al. We find that equipartition is satisfied only if the jet opening angle in the radio core region is close to the values found observationally, $\simeq$0.1--0.2 divided by the bulk Lorentz factor, $Γ_{\rm j}$. Larger values, e.g., $1/Γ_{\rm j}$, would imply magnetic fields much above equipartition. A small jet opening angle implies in turn the magnetization parameter of $\ll 1$. We determine the jet magnetic flux taking into account this effect. We find that the transverse-averaged jet magnetic flux is fully compatible with the model of jet formation due to BH spin energy extraction and the accretion being a magnetically arrested disc (MAD). We calculate the jet average mass-flow rate corresponding to this model and find it consists of a substantial fraction of the mass accretion rate. This suggests the jet composition with a large fraction of baryons. We also calculate the average jet power, and find it moderately exceeds the accretion power, $\dot M c^2$, reflecting BH spin energy extraction. We find our results for radio galaxies at low Eddington ratios are compatible with MADs but require a low radiative efficiency, as predicted by standard accretion models.

preprint2015arXiv

Formation of recollimation shocks in jets of high-mass X-ray binaries

We study conditions for formation of recollimation shocks in jets interacting with stellar winds in high-mass X-ray binaries. We show the existence of a critical jet power, dependent on the wind rate and velocity and the jet velocity, above which a recollimation shock is not formed. For jet powers below critical, we derive the location of the shock. We show that surface shocks may still exist above the critical power, but only occupy a small volume of the jet and do not significantly alter the jet opening angle. We test these prediction by 3-D numerical simulations, which confirm the existence and the value of the critical power. We apply our results to Cyg X-1 and Cyg X-3.

preprint2015arXiv

Hadronic models of blazars require a change of the accretion paradigm

We study hadronic models of broad-band emission of jets in radio-loud active galactic nuclei, and their implications for the accretion in those sources. We show that the models that account for broad-band spectra of blazars emitting in the GeV range in the sample of Boettcher et al. have highly super-Eddington jet powers. Furthermore, the ratio of the jet power to the radiative luminosity of the accretion disc is $\sim 3000$ on average and can be as high as $\sim 10^5$. We then show that the measurements of the radio core shift for the sample imply low magnetic fluxes threading the black hole, which rules out the Blandford-Znajek mechanism to produce powerful jets. These results require that the accretion rate necessary to power the modelled jets is extremely high, and the average radiative accretion efficiency is $\sim 4 \times 10^{-5}$. Thus, if the hadronic model is correct, the currently prevailing picture of accretion in AGNs needs to be significantly revised. Also, the obtained accretion mode cannot be dominant during the lifetimes of the sources, as the modelled very high accretion rates would result in too rapid growth of the central supermassive black holes. Finally, the extreme jet powers in the hadronic model are in conflict with the estimates of the jet power by other methods.

preprint2015arXiv

Long-term quasi-periodicity of 4U 1636-536 resulting from accretion disc instability

We present the results of a study of the low-mass X-ray binary 4U 1636-536. We have performed temporal analysis of all available RXTE/ASM, Swift/BAT and MAXI data. We have confirmed the previously discovered quasi-periodicity of ~45 d present during ~2004, however we found it continued to 2006. At other epochs, the quasi-periodicity is only transient, and the quasi-period, if present, drifts. We have then applied a time-dependent accretion disc model to the interval with the significant X-ray quasi-periodicity. For our best model, the period and the amplitude of the theoretical light curve agree well with that observed. The modelled quasi-periodicity is due to the hydrogen thermal-ionization instability occurring in outer regions of the accretion disc. The model parameters are the average mass accretion rate (estimated from the light curves), and the accretion disc viscosity parameters, for the hot and cold phases. Our best model gives relatively low values of viscosity parameter for cold phase 0.01 and for hot phase 0.03.

preprint2015arXiv

The high-energy gamma-ray detection of G73.9+0.9, a supernova remnant interacting with a molecular cloud

We have analysed the Fermi LAT data on the SNR G73.9+0.9. We have confirmed a previous detection of high-energy gamma-rays from this source at a high significance of $\simeq 12σ$. The observed spectrum shows a significant curvature, peaking in $E F_E$ at $\sim$1 GeV. We have also calculated the flux upper limits in the mm-wavelength and X-ray ranges from Planck and XMM-Newton, respectively. We have inspected the intensity of the CO (1$\rightarrow $0) emission line and found a large peak at a velocity range corresponding to the previously estimated source distance of $\sim$4 kpc, which may indicate an association between a molecular cloud and the SNR. The gamma-ray emission appears due to interaction of accelerated particles within the SNR with the matter of the cloud. The most likely radiative process responsible for the gamma-ray emission is decay of neutral pions produced in ion-ion collisions. While a dominant leptonic origin of this emission can be ruled out, the relativistic electron population related to the observed radio flux will necessarily lead to a certain level of bremsstrahlung gamma-ray emission. Based on this broad-band modelling, we have developed a method to estimate the magnetic field, yielding $B\geq 80\,μ$G at our best estimate of the molecular cloud density (or less at a lower density). G73.9+0.9 appears similar, though somewhat weaker, to other SNRs interacting with a local dense medium detected by the LAT.

preprint2014arXiv

Correlation between the photon index and X-ray luminosity of black hole X-ray binaries and active galactic nuclei: observations and interpretation

We investigate the observed correlation between the 2--10 keV X-ray luminosity (in unit of the Eddington luminosity; $l_X \equiv L_X/L_{Edd}$) and the photon index ($Γ$) of the X-ray spectrum for both black hole X-ray binaries (BHBs) and active galactic nuclei (AGNs). We construct a large sample, with $10^{-9} < l_X < 10^{-1}$. We find that $Γ$ is positively and negatively correlated with $l_X$ when $l_X > 10^{-3}$ and $10^{-6.5} < l_X < 10^{-3}$ respectively, while $Γ$ is nearly a constant when $l_X < 10^{-6.5}$. We explain the above correlation in the framework of a coupled hot accretion flow -- jet model. The radio emission always come from the jet while the X-ray emission comes from the accretion flow and jet when $l_X$ is above and below $10^{-6.5}$, respectively. More specifically, we assume that with the increase of mass accretion rate, the hot accretion flow develops into a clumpy and further a disc -- corona two-phase structure because of thermal instability. We argue that such kind of two-phase accretion flow can explain the observed positive correlation.

preprint2014arXiv

Jet contributions to the broad-band spectrum of Cyg X-1 in the hard state

We apply the jet model developed in the preceding paper of Zdziarski et al. to the hard-state emission spectra of Cyg X-1. We augment the model for the analytical treatment of the particle evolution beyond the energy dissipation region, and allow for various forms of the acceleration rate. We calculate the resulting electron and emission spectra as functions of the jet height, along with the emission spectra integrated over the outflow. The model accounts well for the observed radio, infrared, and GeV fluxes of the source, although the available data do not provide unique constraints on the model free parameters. The contribution of the jet emission in the UV--to--X-ray range turns out to be in all the cases negligible compared to the radiative output of the accretion component. Nevertheless, we find out that it is possible to account for the observed flux of Cyg X-1 at MeV energies by synchrotron jet emission, in accord with the recent claims of the detection of strong linear polarization of the source in that range. However, this is possible only assuming a very efficient particle acceleration leading to the formation of flat electron spectra, and jet magnetic fields much above the equipartition level.

preprint2014arXiv

Jet models for black-hole binaries in the hard spectral state

This is part one of our study of models of jets with distributed electron acceleration. We present here our assumptions, basic equations, and their solutions for the steady-state electron distribution. We assume the shape of the rate of electron acceleration and the dependencies of its normalization and the magnetic field strength on the height along the jet. Our focus is on the hard spectral state of black-hole binaries, for which we take into account that their typical radio spectra are flat. This appears to require a constant dissipation rate per unit logarithmic length and conservation of the magnetic energy flux. Our electron kinetic equation includes adiabatic and radiative losses and advection, and our photon radiative transfer equation includes synchrotron absorption and emission and Compton emission. Apart from the self-Compton process, we take into account Compton scattering of stellar and accretion photons and absorption of very-high energy gamma-rays by pair production on soft photons. We present a general solution of the kinetic equation with advection and radiative and adiabatic losses and an analytic solution in the case of dominant synchrotron losses in conical jets. In the following paper, we present detailed spectra resulting from our equations as applied to Cyg X-1.

preprint2014arXiv

The jet kinetic power, distance and inclination of GRS 1915+105

We apply a recently developed technique of calculating the minimum jet kinetic power to the major mass ejections of the black-hole binary GRS 1915+105 observed in radio wavelengths in 1994 and 1997. We derive for them the distance-dependent minimum power, and the corresponding mass flow rate and the total energy and mass content. We find that a fast increase of the jet power with the increasing distance combined with the jet power estimates based on the bolometric luminosity imply the source distance is <10 kpc. If the jet in GRS 1915 contains ions, their bulk motion dominates the jet power, which was either neglected or not properly taken into account earlier. We also reconsider the parameters of the binary, and derive the current best estimates of the distance-dependent black-hole mass and the inclination based on existing measurements combined with the kinematic constraints from the mass ejections. We also find the measurement of the donor radius of Steeghs et al. implies the distance to the system of <10 kpc, in agreement with the estimate from the jet power.

preprint2014arXiv

The minimum jet power and equipartition

We derive the minimum power of jets and their magnetic field strength based on their observed non-thermal synchrotron emission. The correct form of this method takes into account both the internal energy in the jet and the ion rest-mass energy associated with the bulk motion. The latter was neglected in a number of papers, which instead adopted the well-known energy-content minimization method. That method was developed for static sources, for which there is no bulk-motion component of the energy. In the case of electron power-law spectra with index >2 in ion-electron jets, the rest-mass component dominates. The minimization method for the jet power taking it into account was considered in some other work, but only based on either an assumption of a constant total synchrotron flux or a fixed range of the Lorentz factors. Instead, we base our method on an observed optically-thin synchrotron spectrum. We find the minimum jet power is independent of its radius when the ion rest-mass power dominates, which becomes the case below certain critical radius. This allows for robust minimum power estimates. We also present results for the case with observed turnover frequency, at which the source becomes optically thick. This method allows a determination of the source size, in addition to the power and the magnetic field. We also point out that when the ion rest-mass power dominates, the estimates of the minimum power lead to very different equipartition parameters than those based on minimization of the energy content. The former and latter lead to approximate equipartition between the internal energy in magnetic field and in particles including and excluding, respectively, their rest mass energy.

preprint2013arXiv

Compton scattering of blackbody photons by relativistic electrons

We present simple and accurate analytical formulae for the rates of Compton scattering by relativistic electrons integrated over the energy distribution of blackbody seed photons. Both anisotropic scattering, in which blackbody photons arriving from one direction are scattered by an anisotropic electron distribution into another direction, and scattering of isotropic seed photons are considered. Compton scattering by relativistic electrons off blackbody photons from either stars or CMB takes place, in particular, in microquasars, colliding-wind binaries, supernova remnants, interstellar medium, and the vicinity of the Sun.

preprint2013arXiv

High-energy gamma-ray emission from Cyg X-1 measured by Fermi and its theoretical implications

We have obtained measurements and upper limits on the emission of Cyg X-1 in the photon energy range of 0.03--300 GeV based on observations by Fermi. We present the results separately for the hard and soft spectral states, as well for all of the analysed data. In the hard state, we detect a weak steady emission in the 0.1--10 GeV range with a power-law photon index of Gamma=2.6+-0.2 at a 4 sigma statistical significance. This measurement, even if considered to be an upper limit, strongly constrains Compton emission of the steady radio jet, present in that state. The number of relativistic electrons in the jet has to be low enough for the spectral components due to Compton upscattering of the stellar blackbody and synchrotron radiation to be within the observed fluxes. If optically-thin synchrotron emission of the jet is to account for the MeV tail, as implied by the recently-claimed strong polarization in that energy range, the magnetic field in the jet has to be much above equipartition. The GeV-range measurements also strongly constrain models of hot accretion flows, most likely present in the hard state, in which gamma-rays are produced from decay of neutral pions produced in collisions of energetic ions in an inner part of the flow. In the soft state, the obtained upper limits constrain electron acceleration in a non-thermal corona, most likely present around a blackbody accretion disc. The coronal emission above 30 MeV has to be rather weak, which is most readily explained by absorption of gamma-rays in pair-producing photon-photon collisions. Then, the size of the bulk of the corona is less than a few tens of the gravitational radii.

preprint2012arXiv

Cyg X-3: a low-mass black hole or a neutron star

Cyg X-3 is a highly interesting accreting X-ray binary, emitting from the radio to high-energy gamma-rays. It consists of a compact object wind-fed by a Wolf-Rayet (WR) star, but the masses of the components and the mass-loss rate have been a subject of controversies. Here, we determine its masses, inclination, and the mass-loss rate using our derived relationship between the mass-loss rate and the mass for WR stars of the WN type, published infrared and X-ray data, and a relation between the mass-loss rate and the binary period derivative (observed to be >0 in Cyg X-3). Our obtained mass-loss rate is almost identical to that from two independent estimates and consistent with other ones, which strongly supports the validity of this solution. The found WR and compact object masses are 10.3_{-2.8}^{+3.9}, 2.4_{-1.1}^{+2.1} solar masses, respectively. Thus, our solution still allows for the presence of either a neutron star or a black hole, but the latter only with a low mass. However, the radio, infrared and X-ray properties of the system suggest that the compact object is a black hole. Such a low-mass black-hole could be formed via accretion-induced collapse or directly from a supernova.

preprint2012arXiv

Energy-dependent orbital modulation of X-rays and constraints on emission of the jet in Cyg X-3

We study orbital modulation of X-rays from Cyg X-3, using data from Swift, INTEGRAL and RXTE. Using the wealth of the presently available data and an improved averaging method, we obtain energy-dependent folded and averaged light curves with unprecedented accuracy. We find that above ~5 keV, the modulation depth decreases with the increasing energy, which is consistent with the modulation being caused by both bound-free absorption and Compton scattering in the stellar wind of the donor, with minima corresponding to the highest optical depth, which occurs around the superior conjunction. We find a decrease of the depth below ~3 keV, which appears to be due to re-emission of the absorbed continuum by the wind in soft X-ray lines. Based on the shape of the folded light curves, any X-ray contribution from the jet in Cyg X-3, which emits gamma-rays detected at energies >0.1 GeV in soft spectral states, is found to be minor up to ~100 keV. This implies the presence of a rather sharp low-energy break in the jet MeV-range spectrum. We also calculate phase-resolved RXTE X-ray spectra, and show the difference between the spectra corresponding to phases around the superior and inferior conjunctions can indeed be accounted for by a combined effect of bound-free absorption in an ionized medium and Compton scattering.

preprint2012arXiv

The gamma-ray emitting region of the jet in Cyg X-3

We study models of the gamma-ray emission of Cyg X-3 observed by Fermi. We calculate the average X-ray spectrum during the gamma-ray active periods. Then, we calculate spectra from Compton scattering of a photon beam into a given direction by isotropic relativistic electrons with a power-law distribution, both based on the Klein-Nishina cross section and in the Thomson limit. Applying the results to scattering of stellar blackbody radiation in the inner jet of Cyg X-3, we find that a low-energy break in the electron distribution at a Lorentz factor of ~ 300--1000 is required by the shape of the observed X-ray/gamma-ray spectrum in order to avoid overproducing the observed X-ray flux. The electrons giving rise to the observed \g-rays are efficiently cooled by Compton scattering, and the power-law index of the acceleration process is ~ 2.5--3. The bulk Lorentz factor of the jet and the kinetic power before the dissipation region depend on the fraction of the dissipation power supplied to the electrons; if it is ~ 1/2, the Lorentz factor is ~ 2.5, and the kinetic power is ~ 10^38 erg/s, which represents a firm lower limit on the jet power, and is comparable to the bolometric luminosity of Cyg X-3. Most of the power supplied to the electrons is radiated. The broad band spectrum constrains the synchrotron and self-Compton emission from the gamma-ray emitting electrons, which requires the magnetic field to be relatively weak, with the magnetic energy density < a few times 10^-3 of that in the electrons. The actual value of the magnetic field strength can be inferred from a future simultaneous measurement of the IR and gamma-ray fluxes.

preprint2012arXiv

The MeV spectral tail in Cyg X-1 and optically-thin emission of jets

We study the average X-ray and soft gamma-ray spectrum of Cyg X-1 in the hard spectral state, using data from INTEGRAL. We compare these results with those from CGRO, and find a good agreement. Confirming previous studies, we find the presence of a high-energy MeV tail beyond a thermal-Comptonization spectrum; however, the tail is much softer and weaker than that recently published by Laurent et al. In spite of this difference, the observed high-energy tail could still be due to the synchrotron emission of the jet of Cyg X-1, as claimed by Laurent et al. To test this possibility, we study optically-thin synchrotron and self-Compton emission from partially self-absorbed jets. We develop formalisms for calculating both emission of the jet base (which we define here as the region where the jet starts its emission) and emission of the entire jet. We require the emission to match that observed at the turnover energy. The optically thin emission is dominated by that from the jet base, and it has to become self-absorbed within it at the turnover frequency. We find this implies the magnetic field strength at the jet base of B_0 prop. to z_0^4, where z_0 is the distance of the base from the black-hole centre. The value of B_0 is then constrained from below by the condition that the self-Compton emission is below an upper limit in the GeV range, and from above by the condition that the Poynting flux does not exceed the jet kinetic power. This yields B_0 of the order of ~10^4 G and the location of the jet base at ~10^3 gravitational radii. Using our formalism, we find the MeV tail can be due to jet synchrotron emission, but this requires the electron acceleration at a rather hard power-law index, p~1.3-1.6. For acceleration indices of p> 2, the amplitude of the synchrotron component is much below that of MeV tail, and its origin is likely to be due to hybrid Comptonization in the accretion flow.

preprint2012arXiv

The structure of the jet in Cyg X-1 inferred from orbital modulation of the radio emission

We study free-free absorption of radio emission by winds of massive stars. We derive formulae for the optical depth through the wind measured from a point of emission along a jet, taking into account Compton and photoionization heating and Compton, recombination, line and advection cooling. We apply the developed formalism to radio monitoring data for Cyg X-1, which allows us to obtain strong constraints on the structure of its inner jet. With the data at 15 GHz, and taking into account an anisotropy of the stellar wind in Cyg X-1, we estimate the location of the peak of that emission along the jet at about one orbital separation, i.e., $\sim 3\times 10^{12}$ cm. Given a previous determination of the turnover frequency in Cyg X-1, this implies the location of the base of the jet at $\sim 10^3$ gravitational radii. We also obtain corresponding results at 8.3 GHz and 2.25 GHz, which roughly follow the standard conical partially self-absorbed jet model. Furthermore, we find that the level of the orbital modulation depends on the radio flux, with the modulation being substantially stronger when the radio flux is lower. This is explained by the height of the radio emission along the jet decreasing with the decreasing radio flux, as predicted by jet models. Based on the finding of the flux-dependent orbital modulation, we are able to estimate a range of the possible changes of the form of the radio/X-ray correlation in Cyg X-1 due to free-free absorption. We also derive predictions for the orbital modulation and flux attenuation at frequencies beyond the 2.25--15 GHz range.

preprint2011arXiv

X-ray variability patterns and radio/X-ray correlations in Cyg X-1

We have studied the X-ray variability patterns and correlations of the radio and X-ray fluxes in all spectral states of Cyg X-1 using X-ray data from RXTE/ASM, CGRO/BATSE, and Swift/BAT. In the hard state, the dominant spectral variability is a changing of normalisation with fixed spectral shape, while in the intermediate state the slope changes, with a pivot point around 10 keV. In the soft state, the low energy X-ray emission dominates the bolometric flux which is only loosely correlated with the high energy emission. In black hole binaries in the hard state, the radio flux is generally found to depend on a power of the X-ray flux, F_R ~ F_X^p. We confirm this for Cyg X-1. Our new finding is that this correlation extends to the intermediate and soft states provided the broad-band X-ray flux in the Comptonization part of the spectrum (excluding the blackbody component) is considered instead of a narrow-band medium-energy X-ray flux. We find an index p=1.7+-0.1 for 15 GHz radio emission, decreasing to p=1.5+-0. at 2.25 GHz. We conclude that the higher value at 15 GHz is due to the effect of free-free absorption in the wind from the companion. The intrinsic correlation index remains uncertain. However, based on a theoretical model of the wind in Cyg X-1, it appears to be close to 1.3, which, in the framework of accretion/jet models, implies that the accretion flow in Cyg X-1 is radiatively efficient. The correlation with the flux due to Comptonization emission indicates that the radio jet is launched by the hot electrons in the accretion flow in all spectral states of Cyg X-1. On the other hand, we rule out the X-ray jet model. Finally, we find that the index of the correlation, when measured using the X-ray flux in a narrow energy band, strongly depends on the band chosen and is, in general, different from that for either the bolometric flux or the flux in the hot-electron emission.

preprint2010arXiv

A study of spectral and timing properties of Cyg X-1 based on a large sample of pointed RXTE observations

We study a large sample of RXTE PCA/HEXTE observations of Cyg X-1. We characterize the spectra by soft and hard X-ray colours (which define the spectral states), and fit them with a physical model of hybrid, thermal/non-thermal Comptonization. We then fit the power spectra by a sum of Lorentzians. We show the resulting correlations between the spectral colour (defining the hard, intermediate and soft spectral states), the Comptonization amplification factor, the fractional Compton reflection strength, and the peak frequency of the Lorentzians. We also calculate the fractional variability (rms) as a function of photon energy. We fit the obtained rms dependencies by physical models of varying accretion rate, a varying soft-photon input from an outer disc to a hot inner plasma, and a varying dissipation rate in a hot corona above a disc.

preprint2010arXiv

The doubling of the superorbital period of Cyg X-1

We study properties of the superorbital modulation of the X-ray emission of Cyg X-1. We find that it has had a stable period of about 300 d in soft and hard X-rays and in radio since 2005 until at least 2010, which is about double the previously seen period. This new period, seen in the hard spectral state only, is detected not only in the light curves but also in soft X-ray hardness ratios and in the amplitude of the orbital modulation. On the other hand, the spectral slope in hard X-rays, >20 keV, averaged over superorbital bins is constant, and the soft and hard X-rays and the radio emission change in phase. This shows that the superorbital variability consists of changing the normalization of an intrinsic spectrum of a constant shape and of changes of the absorbing column density with the phase. The maximum column density is achieved at the superorbital minimum. The amplitude changes are likely to be caused by a changing viewing angle of an anisotropic emitter, most likely a precessing accretion disc. The constant shape of the intrinsic spectrum shows that this modulation is not caused by a changing accretion rate. The modulated absorbing column density shows the presence of a bulge around the disc centre, as proposed previously. We also find the change of the superorbital period from about 150 d to about 300 d to be associated with almost unchanged average X-ray fluxes, making the period change difficult to explain in the framework of disc-irradiation models. Finally, we find no correlation of the X-ray and radio properties with the reported detections in the GeV and TeV gamma-ray range.

preprint2009arXiv

Monte Carlo simulations of global Compton cooling in inner regions of hot accretion flows

Hot accretion flows such as advection-dominated accretion flows are generally optically thin in the radial direction. Thus photons generated at some radii can cool or heat electrons at other radii via Compton scattering. Such global Compton scattering has previously been shown to be important for the dynamics of accretion flows. Here, we extend previous treatments of this problem by using accurate global general relativistic Monte Carlo simulations. We focus on an inner region of the accretion flow (R < 600R_g), for which we obtain a global self-consistent solution. As compared to the initial, not self-consistent solution, the final solution has both the cooling rate and the electron temperature significantly reduced at radii >=10 gravitational radii. On the other hand, the radiation spectrum of the self-consistent solution has the shape similar to that of the initial iteration, except for the high-energy cut-off being at an energy lower by a factor of ~2 and the bolometric luminosity decreased by a factor of ~2. We also compare the global Compton scattering model with local models in spherical and slab geometry. We find that the slab model approximates the global model significantly better than the spherical one. Still, neither local model gives a good approximation to the radial profile of the cooling rate, and the differences can be up to two orders of magnitude. The local slab model underestimates the cooling rate at outer regions whereas it overestimates that rate at inner regions. We compare our modelling results to observed hard-state spectra of black-hole binaries and find an overall good agreement provided any disc outflow is weak. We find that general-relativistic effects in flows which dynamics is modified by global Comptonization is crucial in approaching this agreement.

preprint2006arXiv

Modeling the hard states of XTE J1550--564 during its 2000 outburst

We study hard states of the black-hole binary XTE J1550--564 during its 2000 outburst. In order to explain those states at their highest luminosities, $L\sim 10%$ of the Eddington luminosity, $L_{\rm E}$, we propose a specific hot accretion flow model. We point out that the highest values of the hard-state $L$ are substantially above the $L$ an advection-dominated accretion flow (ADAF) can produce, $\sim 0.4α^2 L_{\rm E}$, which is only $\sim (3$--$4)%L_{\rm E}$ even for $α$ as high as 0.3. On the other hand, we successfully explain the hard states with $L\sim (4$--$10)%$ using the luminous hot accretion flow (LHAF) model. As $10%L_{\rm E}$ is also roughly the highest luminosity an LHAF can produce, such an agreement between the predicted and observed highest luminosities provides by itself strong support for this model. Then, we study multi-waveband spectral variability during the 2000 outburst. In addition to the primary maxima in the optical light curves, secondary maxima were detected after the transition from the very high state to the hard state. We show that the secondary maxima are well modeled by synchrotron emission from a jet formed during the state transition. We argue that the absence of the corresponding secondary peak in the X-ray light curve indicates that the X-ray jet emission, regardless of its radiative process, synchrotron or its Comptonization, is not important in the hard state compared to the emission from the accretion flow.

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.

preprint1998arXiv

Hot accretion discs with thermal Comptonisation and advection in luminous black hole sources

We solve for the structure of a hot accretion disc with unsaturated thermal Comptonisation of soft photons and with advection, generalizing the classical model of Shapiro et al. The upper limit on the accretion rate due to advection constrains the luminosity to $\simless 0.15 y^{3/5} α^{7/5}$ of the Eddington limit, where $y$ and $α$ are the Compton and viscosity parameters, respectively. The characteristic temperature and Thomson optical depth of the inner flow at accretion rates within an order of magnitude of that upper limit are $\sim 10^9$ K and $\sim 1$, respectively. The resulting spectra are then in close agreement with the X-ray and soft \g-ray spectra from black-hole binaries in the hard state and Seyferts. At low accretion rates, bremsstrahlung becomes the dominant radiative process.

preprint1997arXiv

Simultaneous X-ray and gamma-ray observations of Cyg X-1 in the hard state by Ginga and OSSE

We present four X-ray/gamma-ray spectra of Cyg X-1 observed in the hard ('low') state simultaneously by Ginga and GRO/OSSE on 1991 July 6. The 3-30 keV Ginga spectra are well represented by power laws with an energy spectral index of alpha~0.6 and a Compton reflection component including a fluorescent Fe K-alpha corresponding to the solid angle of the reflector of ~0.3 times 2 Pi. The overall Ginga/OSSE spectra can be modelled by repeated Compton scattering in a mildly-relativistic, tau ~1, plasma. However, the high-energy cutoff is steeper than that due to single-temperature thermal Comptonisation. It can be described by a superposition of dominant optically-thin, thermal emission at kT~140 keV and a Wien-like component from an optically-thick plasma at kT~50 keV. The X-ray spectra do not show the presence of an anisotropy break required if thermal Compton scattering takes place in a corona above a cold disc. Also, the flat spectral index shows that the plasma is soft-photon starved, i.e., the luminosity in incident soft X-ray seed photons is very much less than that in the hard X-rays. Furthermore, the observed solid angle of the reflector is significantly less than 2 pi. These facts taken together strongly rule out a disc-corona geometry. Rather, the observed spectra are consistent with a geometry in which the cold accretion disc (which both supplies the seed soft X-rays and reflects hard X-rays) only exists at large radii, while the Comptonising hot plasma is located in an inner region with no cold disc. This hot plasma consists of either pure electron-positron pairs if the source size is ~5 Schwarzschild radii or it contains also protons if the size is larger.

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.