Source author record

Marek Sikora

Marek Sikora 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

36works
8topics
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

36 published item(s)

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

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.

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.

preprint2021arXiv

Orientation of the crescent image of M 87*

The first image of the black hole (BH) M 87* obtained by the Event Horizon Telescope (EHT) has the shape of a crescent extending from the E to WSW position angles, while the observed direction of the large-scale jet is WNW. Images based on numerical simulations of BH accretion flows suggest that on average the projected BH spin axis should be oriented SSW. We explore highly simplified toy models for geometric distribution and kinematics of emitting regions in the Kerr metric, perform ray tracing to calculate the corresponding images, and simulate their observation by the EHT to calculate the corresponding visibilities and closure phases. We strictly assume that (1) the BH spin vector is fixed to the jet axis, (2) the emitting regions are stationary and symmetric with respect to the BH spin, and that (3) the emissivities are isotropic in the local rest frames. Emission from the crescent sector between SSE and WSW can be readily explained in terms of an equatorial ring with either circular or plunging geodesic flows, regardless of the value of BH spin. In the case of plane-symmetric polar caps with plunging geodesic flows, the dominant image of the cap located behind the BH is sensitive to the angular momentum of the emitter. Within the constraints of our model, we have not found a viable explanation for the observed brightness of the ESE sector. Most likely, the ESE "hotspot" has been produced by a non-stationary localised perturbation in the inner accretion flow. Alternatively, it could result from locally anisotropic synchrotron emissivities. Multi-epoch and polarimetric results from the EHT will be essential to verify the theoretically expected alignment of the BH spin with the large-scale jet.

preprint2020arXiv

On the diversity of the jet production efficiency in Swift/BAT AGNs

We address the very large diversity of the jet production efficiency in active galactic nuclei (AGNs) by using data on low redshift AGNs selected from the Swift/BAT catalog and having black hole (BH) masses larger than $10^{8.5}\,M_{\odot}$. Most of these AGNs accrete at intermediate rates and have bolometric luminosities dominated by mid-IR radiation. Our sample contains $14\%$ radio-loud (RL), $6\%$ radio-intermediate (RI), and $80\%$ radio-quiet (RQ) AGNs. All RL objects are found to have extended radio structures and most of them have classical FR II morphology. Converting their radio loudness to the jet production efficiency, we find that the median of this efficiency is on the order of $(ε_d/0.1)\%$, where $ε_d=L_{\rm bol}/\dot{M}c^2$ is the radiation efficiency of the accretion disk. Without knowing the contribution of jets to the radio emission in the RQ AGNs, we are only able to estimate their efficiencies using upper limits. Their median is found to be $0.002(ε_d/0.1)\%$. Our results suggest that some threshold conditions must be satisfied to allow production of strong, relativistic jets in RL AGNs. We discuss several possible scenarios and argue that the production of collimated, relativistic jets must involve the Blandford-Znajek mechanism and can be activated only in those AGNs whose lifetime is longer than the time required to enter the magnetically arrested disk (MAD). Presuming that MAD is required to collimate relativistic jets, we expect that the weak nonrelativistic jets observed in some RQ AGNs are produced by accretion disks rather than by rotating BHs.

preprint2019arXiv

Comparison of SEDs of very massive radio-loud and radio-quiet AGN

The main objective of this work is to establish and interpret the dominant spectral components and their differences in radio-loud (RL) and radio-quiet (RQ) AGN with very massive black holes, and accreting at moderate rates. Such a sample is selected from the \textit{Swift}/BAT catalogue of AGN having determined optical spectra types and hosting black holes with masses $\rm > 10^{8.5}\,M_{\odot}$. We confirm our previous results, that radio loudness distribution of \textit{Swift}/BAT AGN is bimodal and that radio galaxies are about two times X-ray louder than their radio-quiet counterparts. We show that the average X-ray loudness (defined as a ratio of luminosity in the $14-195$ keV band to that at 12 $μ$m) of Type 1 and Type 2 AGN is very similar. This similarity holds for both RL and RQ subsamples and indicates negligible dependence of the observed X-ray luminosities on the inclination angle in both populations. In both the radiative output is dominated by mid-IR and hard X-ray components, and relatively weak UV luminosities indicate large amounts of dust in polar regions.

preprint2019arXiv

RuleKit: A Comprehensive Suite for Rule-Based Learning

Rule-based models are often used for data analysis as they combine interpretability with predictive power. We present RuleKit, a versatile tool for rule learning. Based on a sequential covering induction algorithm, it is suitable for classification, regression, and survival problems. The presence of a user-guided induction facilitates verifying hypotheses concerning data dependencies which are expected or of interest. The powerful and flexible experimental environment allows straightforward investigation of different induction schemes. The analysis can be performed in batch mode, through RapidMiner plug-in, or R package. A documented Java API is also provided for convenience. The software is publicly available at GitHub under GNU AGPL-3.0 license.

preprint2016arXiv

A spine-sheath model for strong-line blazars

We have developed a quasi-analytical model for the production of radiation in strong-line blazars, assuming a spine-sheath jet structure. The model allows us to study how the spine and sheath spectral components depend on parameters describing the geometrical and physical structure of "the blazar zone". We show that typical broad-band spectra of strong-line blazars can be reproduced by assuming the magnetization parameter to be of order unity and reconnection to be the dominant dissipation mechanism. Furthermore, we demonstrate that the spine-sheath model can explain why gamma-ray variations are often observed to have much larger amplitudes than the corresponding optical variations. The model is also less demanding of jet power than one-zone models, and can reproduce the basic features of extreme gamma-ray events.

preprint2016arXiv

Covering factors of the dusty obscurers in radio-loud and radio-quiet quasars

We compare covering factors of circumnuclear dusty obscurers in radio-loud and radio-quiet quasars. The radio-loud quasars are represented by a sample of FR II quasars obtained by cross-matching a catalog of the FR II radio sources selected by van Velzen et al. with the SDSS DR7 catalog of quasars. Covering factors of FR II quasars are compared with covering factors of the radio-quiet quasars matched with them in redshift, black hole mass, and Eddington-ratio. We found that covering factors, proxied by the infrared-to-bolometric luminosity ratio, are on average slightly smaller in FR II quasars than in radio-quiet quasars, however, this difference is statistically significant only for the highest Eddington ratios. For both samples, no statistically significant dependence of a median covering factor on Eddington ratio, black hole mass, nor redshift can be claimed.

preprint2016arXiv

Powers and Magnetization of Blazar Jets

In this work I review the observational constraints imposed on the energetics and magnetisation of quasar jets, in the context of theoretical expectations. The discussion is focused on issues regarding the jet production efficiency, matter content, and particle acceleration. I show that if the ratio of electron-positron-pairs to protons is of order $15$, as is required to achieve agreement between jet powers computed using blazar spectral fits and those computed using radio-lobe calorimetry, the magnetization of blazar jets in flat-spectrum-radio-quasars (FSRQ) must be significant. This result favors the reconnection mechanism for particle acceleration and explains the large Compton-dominance of blazar spectra that is often observed, without the need to postulate very low jet magnetization.

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.

preprint2015arXiv

Application of the spine-layer jet radiation model to outbursts in the broad-line radio galaxy 3C 120

We present a detailed Fermi/LAT data analysis for the broad-line radio galaxy 3C 120. This source has recently entered into a state of increased gamma-ray activity which manifested itself in two major flares detected by Fermi/LAT in September 2014 and April 2015 with no significant flux changes reported in other wavelengths. We analyse available data focusing our attention on aforementioned outbursts. We find very fast variability timescale during flares (of the order of hours) together with a significant gamma-ray flux increase. We show that the 6.8 years averaged gamma-ray emission of 3C 120 is likely a sum of the external radiation Compton and the synchrotron self-Compton radiative components. To address the problem of "orphan" gamma-ray flares and fast variability we model the jet radiation dividing the jet structure into two components: the wide and relatively slow outer layer and the fast, narrow spine. We show that with the addition of the fast spine occasionally bent towards the observer we are able to explain observed spectral energy distribution of 3C 120 during flares with the Compton up-scattered broad-line region and dusty torus photons as main gamma-rays emission mechanism.

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

Magnetization of jets in luminous blazars

Luminosities of many powerful blazars are strongly dominated by $γ$-rays which most likely result from Comptonization of radiation produced outside a jet. This observation sets certain constraints on composition and energetics of the jet as well as the surrounding quasar environment. We study the dependence of a Compton dominance on a jet magnetization (the magnetic-to-matter energy flux) and on the location of the 'blazar zone'. Calculations are performed for two geometries of broad-emission-line and hot-dust regions: spherical and planar. The jet magnetization corresponding to the large observed Compton dominance is found to be $\sim 0.1 (θ_{\rm j} Γ)^2$ for spherical geometries and $\sim 0.01 (θ_{\rm j} Γ)^2$ for planar geometries, where $θ_{\rm j}$ is the jet half-opening angle and $Γ$ is the jet Lorentz factor. This implies that jets in luminous blazars are matter dominated and that this domination is particularly strong for the flattened geometry of external radiation sources.

preprint2015arXiv

Time-Dependent Modeling of Gamma-ray Flares in Blazar PKS1510-089

Here we present a new approach for constraining luminous blazars, incorporating fully time-dependent and self-consistent modeling of bright gamma-ray flares of PKS1510-089 resolved with Fermi-LAT, in the framework of the internal shock scenario. The results of our modeling imply the location of the gamma-ray flaring zone outside of the broad-line region, namely around 0.3pc from the core for a free-expanding jet with the opening angle Gamma, θ_\mathrm{jet} \simeq 1 (where Gamma is the jet bulk Lorentz factor), up to \simeq 3pc for a collimated outflow with Gamma, θ_\mathrm{jet} \simeq 0.1. Moreover, under the Gamma, θ_\mathrm{jet} \simeq 1 condition, our modeling indicates the maximum efficiency of the jet production during the flares, with the total jet energy flux strongly dominated by protons and exceeding the available accretion power in the source. This is in contrast to the quiescence states of the blazar, characterized by lower jet kinetic power and an approximate energy equipartition between different plasma constituents. We demostrate how strictly simultaneous observations of flaring PKS1510-089 at optical, X-ray, and GeV photon energies on hourly timescales, augmented by extensive simulations as presented in this paper, may help to impose further precise constraints on the magnetization and opening angle of the emitting region. Our detailed modeling implies in addition that a non-uniformity of the Doppler factor across the jet, caused by the radial expansion of the outflow, may lead to a pronounced time distortion in the observed gamma-ray light curves, resulting in particular in asymmetric flux profiles with substantially extended decay phases.

preprint2014arXiv

Constraining the Location of Gamma-Ray Flares in Luminous Blazars

Locating the gamma-ray emission sites in blazar jets is a long-standing and highly controversial issue. We investigate jointly several constraints on the distance scale r and Lorentz factor Gamma of the gamma-ray emitting regions in luminous blazars (primarily flat spectrum radio quasars, FSRQs). Working in the framework of one-zone external radiation Comptonization (ERC) models, we perform a parameter space study for several representative cases of actual gamma-ray flares in their multiwavelength context. We find a particularly useful combination of three constraints: from an upper limit on the collimation parameter Gamma*theta <~ 1, from an upper limit on the synchrotron self-Compton (SSC) luminosity L_SSC <~ L_X, and from an upper limit on the efficient cooling photon energy E_cool,obs <~ 100 MeV. These three constraints are particularly strong for sources with low accretion disk luminosity L_d. The commonly used intrinsic pair-production opacity constraint on Gamma is usually much weaker than the SSC constraint. The SSC and cooling constraints provide a robust lower limit on the collimation parameter Gamma*theta >~ 0.1 - 0.7. Typical values of r corresponding to moderate values of Gamma ~ 20 are in the range 0.1 - 1 pc, and are determined primarily by the observed variability time scale t_var,obs. Alternative scenarios motivated by the observed gamma-ray/mm connection, in which gamma-ray flares of t_var,obs ~ a few days are located at r ~ 10 pc, are in conflict with both the SSC and cooling constraints. Moreover, we use a simple light travel time argument to point out that the gamma-ray/mm connection does not provide a significant constraint on the location of gamma-ray flares. We argue that spine-sheath models of the jet structure do not offer a plausible alternative to external radiation fields at large distances, however, an extended broad-line region is an idea worth exploring.

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

Reconciling models of luminous blazars with magnetic fluxes determined by radio core shift measurements

Estimates of magnetic field strength in relativistic jets of active galactic nuclei (AGN), obtained by measuring the frequency-dependent radio core location, imply that the total magnetic fluxes in those jets are consistent with the predictions of the magnetically-arrested disk (MAD) scenario of jet formation. On the other hand, the magnetic field strength determines the luminosity of the synchrotron radiation, which forms the low-energy bump of the observed blazar spectral energy distribution (SED). The SEDs of the most powerful blazars are strongly dominated by the high-energy bump, which is most likely due to the external radiation Compton (ERC) mechanism. This high Compton dominance may be difficult to reconcile with the MAD scenario, unless 1) the geometry of external radiation sources (broad-line region, hot-dust torus) is quasi-spherical rather than flat, or 2) most gamma-ray radiation is produced in jet regions of low magnetization, e.g., in magnetic reconnection layers or in fast jet spines.

preprint2013arXiv

Constraining jet production scenarios by studies of Narrow-Line-Radio-Galaxies

We study a large sample of narrow-line radio galaxies (NLRGs) with extended radio structures. Using 1.4 GHz radio luminosities, $L_{1.4}$, narrow optical emission line luminosities, $L_{\oiii}$ and $L_{H_α}$, as well as black hole masses $M_{BH}$ derived from stellar velocity dispersions measured from the optical spectra obtained with the Sloan Digital Sky Survey, we find that: (i) NLRGs cover about 4 decades of the Eddington ratio, $λ\equiv L_{bol}/L_{Edd} \propto L_{line}/M_{BH}$; (ii) $L_{1.4}/M_{BH}$ strongly correlates with $λ$; (iii) radio-loudness, ${\cal R} \equiv L_{1.4}/L_{line}$, strongly anti-correlates with $λ$. A very broad range of the Eddington ratio indicates that the parent population of NLRGs includes both radio-loud quasars (RLQs) and broad-line radio galaxies (BLRGs). The correlations they obey and their high jet production efficiencies favor a jet production model which involves the so-called 'magnetically choked' accretion scenario. In this model, production of the jet is dominated by the Blandford-Znajek mechanism, and the magnetic fields in the vicinity of the central black hole are confined by the ram pressure of the accretion flow. Since large net magnetic flux accumulated in central regions of the accretion flow required by the model can take place only via geometrically thick accretion, we speculate that the massive, 'cold' accretion events associated with luminous emission-line AGN can be accompanied by an efficient jet production only if preceded by a hot, very sub-Eddington accretion phase.

preprint2013arXiv

Magnetic flux paradigm for radio-loudness of AGN

We argue that the magnetic flux threading the black hole, rather than black hole spin or Eddington ratio, is the dominant factor in launching powerful jets and thus determining the radio loudness of active galactic nuclei (AGN). Most AGN are radio quiet because the thin accretion disks that feed them are inefficient in depositing magnetic flux close to the black hole. Flux accumulation is more likely to occur during a hot accretion (or thick disk) phase, and we argue that radio-loud quasars and strong emission-line radio galaxies occur only when a massive, cold accretion event follows an episode of hot accretion. Such an event might be triggered by the merger of a giant elliptical galaxy with a disk galaxy. This picture supports the idea that flux accumulation can lead to the formation of a so-called magnetically choked accretion flow (MCAF). The large observed range in radio loudness reflects not only the magnitude of the flux pressed against the black hole, but also the decrease in UV flux from the disk, due to its disruption by the "magnetosphere" associated with the accumulated flux. While the strongest jets result from the secular accumulation of flux, moderate jet activity can also be triggered by fluctuations in the magnetic flux deposited by turbulent, hot inner regions of otherwise thin accretion disks, or by the dissipation of turbulent fields in accretion disk coronae. These processes could be responsible for jet production in Seyferts and low-luminosity AGN, as well as jets associated with X-ray binaries.

preprint2013arXiv

On the origin of X-ray spectra in luminous blazars

Gamma-ray luminosities of some quasar-associated blazars imply jet powers reaching values comparable to the accretion power even if assuming very strong Doppler boosting and very high efficiency of gamma-ray production. With much lower radiative efficiencies of protons than of electrons, and the recent reports of very strong coupling of electrons with shock-heated protons indicated by Particle-in-Cell (PIC) simulations, the leptonic models seem to be strongly favored over the hadronic ones. However, the electron-proton coupling combined with the ERC (External-Radiation-Compton) models of gamma-ray production in leptonic models predict extremely hard X-ray spectra, with energy indices about 0. This is inconsistent with the observed 2-10 keV slopes of blazars, which cluster around an index value of 0.6. This problem can be resolved by assuming that electrons can be cooled down radiatively to non-relativistic energies, or that blazar spectra are entirely dominated by the SSC (Synchrotron-Self Compton) component up to at least 10 keV. Here, we show that the required cooling can be sufficiently efficient only at distances r < 0.03pc. SSC spectra, on the other hand, can be produced roughly co-spatially with the observed synchrotron and ERC components, which are most likely located roughly at a parsec scale. We show that the dominant SSC component can also be produced much further than the dominant synchrotron and ERC components, at distances larger than 10 parsecs. Hence, depending on the spatial distribution of the energy dissipation along the jet, one may expect to see gamma-ray/optical events with either correlated or uncorrelated X-rays. In all cases the number of electron-positron pairs per proton is predicted to be very low. The direct verification of the proposed SSC scenario requires sensitive observations in the hard X-ray band which is now possible with the NuSTAR satellite.

preprint2013arXiv

The evolution of active galactic nuclei and their spins

Massive black holes (MBHs) in contrast to stellar mass black holes are expected to substantially change their properties over their lifetime. MBH masses increase by several order of magnitude over the Hubble time, as illustrated by Soltan's argument. MBH spins also must evolve through the series of accretion and mergers events that grow the MBH's masses. We present a simple model that traces the joint evolution of MBH masses and spins across cosmic time. Our model includes MBH-MBH mergers, merger-driven gas accretion, stochastic fueling of MBHs through molecular cloud capture, and a basic implementation of accretion of recycled gas. This approach aims at improving the modeling of low-redshift MBHs and AGN, whose properties can be more easily estimated observationally. Despite the simplicity of the model, it captures well the global evolution of the MBH population from z\sim6 to today. Under our assumptions, we find that the typical spin and radiative efficiency of MBHs decrease with cosmic time because of the higher incidence of stochastic processes in gas-rich galaxies and MBH-MBH mergers in gas-poor galaxies. At z=0 the spin distribution in gas-poor galaxies peaks at spins 0.4-0.8, and it is not strongly mass dependent. MBHs in gas-rich galaxies have a more complex evolution, with low-mass MBHs at low redshift having low spins, and spins increasing at larger masses and redshifts. We also find that at z>1 MBH spins are on average highest in high luminosity AGN, while at lower redshifts these differences disappear.

preprint2012arXiv

Energetic Constraints on a Rapid Gamma-Ray Flare in PKS 1222+216

We study theoretical implications of a rapid Very-High-Energy (VHE) flare detected by MAGIC in the Flat-Spectrum Radio Quasar PKS 1222+216. The minimum distance from the jet origin at which this flare could be produced is 0.5 pc. A moderate Doppler factor of the VHE source, D_{VHE} ~ 20, is allowed by all opacity constraints. The concurrent High-Energy (HE) emission observed by Fermi provides estimates of the total jet power and the jet magnetic field strength. Energetic constraints for the VHE flare are extremely tight: for an isotropic particle distribution they require a huge co-moving energy density in the emitting region and a very efficient radiative process. We disfavor hadronic processes due to their low radiative efficiency, as well as the synchrotron scenario recently proposed for the case of HE flares in the Crab Nebula, since the parameters needed to overcome the radiative losses are quite extreme. The VHE emission can be explained by the Synchrotron Self-Compton (SSC) mechanism for D_{VHE} ~ 20 or by the External Radiation Compton (ERC) mechanism involving the infrared radiation of the dusty torus for D_{VHE} ~ 50. After discussing several alternative scenarios, we propose that the extreme energy density constraint can be satisfied when the emission comes from highly anisotropic short-lived bunches of particles formed by the kinetic beaming mechanism in magnetic reconnection sites. By focusing the emitting particles into very narrow beams, this mechanism allows one to relax the causality constraint on the source size, decreasing the required energy density by 4 orders of magnitude.

preprint2012arXiv

Herschel PACS and SPIRE observations of blazar PKS 1510-089: a case for two blazar zones

We present the results of observations of blazar PKS 1510-089 with the Herschel Space Observatory PACS and SPIRE instruments, together with multiwavelength data from Fermi/LAT, Swift, SMARTS and SMA. The source was found in a quiet state, and its far-infrared spectrum is consistent with a power-law with a spectral index of alpha ~ 0.7. Our Herschel observations were preceded by two 'orphan' gamma-ray flares. The near-infrared data reveal the high-energy cut-off in the main synchrotron component, which cannot be associated with the main gamma-ray component in a one-zone leptonic model. This is because in such a model the luminosity ratio of the External-Compton and synchrotron components is tightly related to the frequency ratio of these components, and in this particular case an unrealistically high energy density of the external radiation would be implied. Therefore, we consider a well-constrained two-zone blazar model to interpret the entire dataset. In this framework, the observed infrared emission is associated with the synchrotron component produced in the hot-dust region at the supra-pc scale, while the gamma-ray emission is associated with the External-Compton component produced in the broad-line region at the sub-pc scale. In addition, the optical/UV emission is associated with the accretion disk thermal emission, with the accretion disk corona likely contributing to the X-ray emission.

preprint2012arXiv

On the origin of the gamma-ray/optical lags in luminous blazars

Blazars are strongly variable sources that occasionally show spectacular flares visible in various energy bands. These flares are often, but not always, correlated. In a number of cases the peaks of optical flares are found to be somewhat delayed with respect to the gamma-ray peaks. One notable example of such a delay was found in 3C 279 by Hayashida et al. and interpreted as a result of steeper drop with distance of the energy density of external radiation field than of the magnetic energy density. In this paper we demonstrate that in general, depending on the respective energy density profile along the jet, such lags can have both signs and that they can take place for any ratio of these energy densities. We study the dependence of such lags on the ratio of these energy densities at a distance of a maximal energy dissipation in a jet, on their gradients, as well as on the time profile of the relativistic electron injection within the moving source. We show how prominent such lags can be, and what are their expected time scales. We suggest that studies of such lags can provide a powerful tool to resolve the structure of relativistic jets and their radiative environment. As an example we model the lag observed in 3C 279, showing that in this object the flare is produced at a distance of a few parsecs from the central black hole, consistent with our previous inferences based on the spectra and optical polarization properties.

preprint2012arXiv

Polarization of synchrotron emission from relativistic reconfinement shocks with ordered magnetic fields

We calculate the polarization of synchrotron radiation produced at the relativistic reconfinement shocks, taking into account globally ordered magnetic field components, in particular toroidal and helical fields. In these shocks, toroidal fields produce high parallel polarization (electric vectors parallel to the projected jet axis), while chaotic fields generate moderate perpendicular polarization. Helical fields result in a non-axisymmetric distribution of the total and polarized brightness. For a diverging downstream velocity field, the Stokes parameter U does not vanish and the average polarization is neither strictly parallel nor perpendicular. A distance at which the downstream flow is changing from diverging to converging can be easily identified on polarization maps as the turning point, at which polarization vectors switch, e.g., from clockwise to counterclockwise.

preprint2012arXiv

Quasi-star jets as unidentified gamma-ray sources

Gamma-ray catalogs contain a considerable amount of unidentified sources. Many of these are located out of the Galactic plane and therefore may have extragalactic origin. Here we assume that the formation of massive black holes in galactic nuclei proceeds through a quasi-star stage and consider the possibility of jet production by such objects. Those jets would be the sources of collimated synchrotron and Compton emission, extending from radio to gamma rays. The expected lifetimes of quasi-stars are of the order of million of years while the jet luminosities, somewhat smaller than that of quasar jets, are sufficient to account for the unidentified gamma-ray sources. The jet emission dominates over the thermal emission of a quasi-star in all energy bands, except when the jet is not directed towards an observer. The predicted synchrotron emission peaks in the IR band, with the flux close to the limits of the available IR all sky surveys. The ratio of the $γ$-ray flux to the IR flux is found to be very large ($\sim 60$), much larger than in BL Lac objects but reached by some radio-loud quasars. On the other hand, radio-loud quasars show broad emission lines while no such lines are expected from quasi-stars. Therefore the differentiation between various scenarios accounting for the unidentified gamma-ray sources will be possible at the basis of the photometry and spectroscopy of the IR/optical counterparts.

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.

preprint2010arXiv

Hadronic jet models today

The matter content of relativistic jets in AGNs is dominated by a mixture of protons, electrons, and positrons. During dissipative events these particles tap a significant portion of the internal and/or kinetic energy of the jet and convert it into electromagnetic radiation. While leptons - even those with only mildly relativistic energies - can radiate efficiently, protons need to be accelerated up to energies exceeding $10^{16-19}$ eV to dissipate radiatively a significant amount of energy via either trigerring pair cascades or direct synchrotron emission. Here I review various constraints imposed on the role of hadronic non-adiabatic cooling processes in shaping the high energy spectra of blazars. It will be argued that protons, despite being efficiently accelerated and presumably playing a crucial role in jet dynamics and dissipation of the jet kinetic energy to the internal energy of electrons and positrons, are more likely to remain radiatively passive in AGN jets.

preprint2010arXiv

Launching of jets by cold, magnetized disks in Kerr Metric

We confirm recent discovery by Cao that in the vicinity of fast rotating black holes jets can be launched centrifugally by cold, magnetized disks even for nearly vertically shaped magnetic flux surfaces. Outflows produced under such extreme conditions are investigated via studying kinematics of test particles in the force-free magnetosphere approximation. Implications of a possibility of magneto-centrifugal launching of very well collimated central outflows around the fast rotating black holes are discussed in the general context of the jet formation scenarios in AGNs.

preprint2010arXiv

Radiative Properties of Reconnection-Powered Minijets in Blazars

We construct a numerical model of emission from minijets, localized flows driven by magnetic reconnection inside Poynting-flux-dominated jets proposed to explain the ultrafast variability of blazars. The geometrical structure of the model consists of two wedge-like regions of relativistically flowing gas, separated by a stationary shock. The dynamics is based on solutions of relativistic magnetic reconnection with a guide field from Lyubarsky (2005). Electron distributions in each region are chosen to the match the pressure and density of the local plasma. Synchrotron emission from both regions is used to calculate Compton scattering, Compton drag and photon-photon opacity effects, with exact treatment of anisotropy and the Klein-Nishina regime. Radiative effects on plasma are taken into account, including the dependence of pressure on electron radiative losses and adiabatic heating of the flow decelerating under Compton drag. The results are applied to the July 2006 flare in the BL Lac object PKS 2155-304, with the aim of matching TeV flux measurements by H.E.S.S. with models that satisfy the variability constraints, while keeping X-ray emission below simultaneous Chandra observations. We find that models of isolated minijets with a significant guide field overproduce X-ray emission, and that we must take into account the radiative interaction of oppositely-oriented minijets in order to achieve a high enough dominance by Comptonized TeV radiation. We argue that such interactions are likely to occur in a jet where there is substantial internal reconnection, producing a large number of misaligned minijets. Finally, we show that large jet magnetizations are indeed required to satisfy all observational constraints and that the effective Lorentz factor of the minijet plasma has to be larger than 50, in agreement with earlier one-zone estimates.

preprint2009arXiv

Constraining Emission Models of Luminous Blazar Sources

Many luminous blazars which are associated with quasar-type active galactic nuclei display broad-band spectra characterized by a large luminosity ratio of their high-energy (gamma-ray) and low-energy (synchrotron) spectral components. This large ratio, reaching values up to 100, challenges the standard synchrotron self-Compton models by means of substantial departures from the minimum power condition. Luminous blazars have also typically very hard X-ray spectra, and those in turn seem to challenge hadronic scenarios for the high energy blazar emission. As shown in this paper, no such problems are faced by the models which involve Comptonization of radiation provided by a broad line-region, or dusty molecular torus. The lack or weakness of bulk Compton and Klein-Nishina features indicated by the presently available data favors production of gamma-rays via up-scattering of infrared photons from hot dust. This implies that the blazar emission zone is located at parsec-scale distances from the nucleus, and as such is possibly associated with the extended, quasi-stationary reconfinement shocks formed in relativistic outflows. This scenario predicts characteristic timescales for flux changes in luminous blazars to be days/weeks, consistent with the variability patterns observed in such systems at infrared, optical and gamma-ray frequencies. We also propose that the parsec-scale blazar activity can be occasionally accompanied by dissipative events taking place at sub-parsec distances and powered by internal shocks and/or reconnection of magnetic fields. These could account for the multiwavelength intra-day flares occasionally observed in powerful blazars sources.

preprint2008arXiv

Suzaku Observations of Extreme MeV Blazar SWIFT J0746.3+2548

We report the Suzaku observations of the high luminosity blazar SWIFT J0746.3+2548 (J0746) conducted in November 2005. This object, with z = 2.979, is the highest redshift source observed in the Suzaku Guaranteed Time Observer (GTO) period, is likely to show high gamma-ray flux peaking in the MeV range. As a result of the good photon statistics and high signal-to-noise ratio spectrum, the $Suzaku$ observation clearly confirms that J0746 has an extremely hard spectrum in the energy range of 0.3-24 keV, which is well represented by a single power-law with a photon index of 1.17 and Galactic absorption. The multiwavelength spectral energy distribution of J0746 shows two continuum components, and is well modeled assuming that the high-energy spectral component results from Comptonization of the broad-line region photons. In this paper we search for the bulk Compton spectral features predicted to be produced in the soft X-ray band by scattering external optical/UV photons by cold electrons in a relativistic jet. We discuss and provide constraints on the pair content resulting from the apparent absence of such features.

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.