Showing 1–50 of 123 results for author: Sikora, M

Observations of a black hole X-ray binary indicate formation of a magnetically arrested disk

Authors: Bei You, Xinwu Cao, Zhen Yan, Jean-Marie Hameury, Bozena Czerny, Yue Wu, Tianyu Xia, Marek Sikora, Shuang-Nan Zhang, Pu Du, Piotr T. Zycki

Abstract: Accretion of material onto a black hole drags any magnetic fields present inwards, increasing their strength. Theory predicts that sufficiently strong magnetic fields can halt the accretion flow, producing a magnetically arrested disk (MAD). We analyze archival multi-wavelength observations of an outburst from the black hole X-ray binary MAXI J1820+070 in 2018. The radio and optical fluxes are del… ▽ More Accretion of material onto a black hole drags any magnetic fields present inwards, increasing their strength. Theory predicts that sufficiently strong magnetic fields can halt the accretion flow, producing a magnetically arrested disk (MAD). We analyze archival multi-wavelength observations of an outburst from the black hole X-ray binary MAXI J1820+070 in 2018. The radio and optical fluxes are delayed by about 8 and 17 days respectively, compared to the X-ray flux. We interpret this as evidence for the formation of a MAD. In this scenario, the magnetic field is amplified by an expanding corona, forming a MAD around the time of the radio peak. The optical delay is then due to thermal viscous instability in the outer disk. △ Less

Submitted 11 October, 2023; v1 submitted 31 August, 2023; originally announced September 2023.

Comments: The author's version of the article which will appear in Science on 31 August 2023, 49 pages including the extended data. The online publication version can be found at https://doi.org/10.1126/science.abo4504

Journal ref: Science 381, 961-964 (2023)

Formation and evolution of transient jets and their cavities in black-hole X-ray binaries

Authors: Marek Sikora, Andrzej Zdziarski

Abstract: We propose a model explaining the origin of transient/episodic jets in black-hole X-ray binaries, in which they are caused by transitions from a collimated, strongly magnetized, jet to a wide, un-collimated, outflow. The change occurs when the accretion flow leaves the magnetically-choked state due to an increase of the accretion rate for a weakly varying magnetic flux. The formed powerful jet the… ▽ More We propose a model explaining the origin of transient/episodic jets in black-hole X-ray binaries, in which they are caused by transitions from a collimated, strongly magnetized, jet to a wide, un-collimated, outflow. The change occurs when the accretion flow leaves the magnetically-choked state due to an increase of the accretion rate for a weakly varying magnetic flux. The formed powerful jet then detaches from its base, and propagates as a discrete ejection. The uncollimated outflow then produces a relativistic plasma that fills surrounding of the black hole, contributing to the formation of a low-density cavity. While the pressure in the cavity is in equilibrium with the surrounding interstellar medium (ISM), its inertia is orders of magnitude lower than that of the ISM. This implies that the plasma cannot efficiently decelerate the ejecta, explaining most of the observations. The modest deceleration within the cavities observed in some cases can be then due to the presence of clouds and/or filaments, forming a wide transition zone between the cavity and the ISM. △ Less

Submitted 12 August, 2023; v1 submitted 6 July, 2023; originally announced July 2023.

Comments: ApJL in press, 6 pages

No Need for an Extreme Jet Energy in the Black-Hole X-Ray Binary MAXI J1348--630

Authors: Andrzej A. Zdziarski, Marek Sikora, Michal Szanecki, Markus Boettcher

Abstract: We model interaction with the surrounding medium of the main discrete jet ejection in the accreting black-hole binary MAXI J1348--630. The kinetic energy in the ejection of that jet was estimated before to be $>10^{46}$ erg. That energy requires that the jet power was about two orders of magnitude above the limit corresponding to a magnetically arrested accretion onto a maximally rotating black ho… ▽ More We model interaction with the surrounding medium of the main discrete jet ejection in the accreting black-hole binary MAXI J1348--630. The kinetic energy in the ejection of that jet was estimated before to be $>10^{46}$ erg. That energy requires that the jet power was about two orders of magnitude above the limit corresponding to a magnetically arrested accretion onto a maximally rotating black hole. That large estimate was obtained by considering the initial ballistic jet propagation in a surrounding cavity followed by a sudden deceleration in interstellar medium under the assumption of its standard density of $\sim$1 cm$^{-3}$. Such densities are likely in the surrounding of this source given its location in the Galactic Plane. Here, we show that the estimate of the kinetic energy can be reduced to realistic values of $\sim\! 10^{44}$ erg by considering the presence of a transition layer with an exponential density growth separating the cavity and the interstellar medium. In that case, the jet is found to decelerate mostly in the transition layer, in regions with the densities $\ll$1 cm$^{-3}$, which strongly reduces the energy requirement. Still, the required jet masses are large, ruling out the presence of a significant number of electron-positron pairs. △ Less

Submitted 7 April, 2023; v1 submitted 2 March, 2023; originally announced March 2023.

Comments: ApJL, in press

A simple analytical model of magnetic jets

Authors: Andrzej A. Zdziarski, Lukasz Stawarz, Marek Sikora, Krzysztof Nalewajko

Abstract: 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… ▽ More 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. △ Less

Submitted 31 May, 2022; v1 submitted 25 April, 2022; originally announced April 2022.

Comments: MNRAS, Letters, in press

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

Authors: Andrzej A. Zdziarski, Dakalo G. Phuravhathu, Marek Sikora, Markus Boettcher, James O. Chibueze

Abstract: 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… ▽ More 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. △ Less

Submitted 8 March, 2022; v1 submitted 22 February, 2022; originally announced February 2022.

Comments: ApJL, in press

The dependence of the fraction of radio luminous quasars on redshift and its theoretical implications

Authors: Katarzyna Rusinek-Abarca, Marek Sikora

Abstract: While radio emission in quasars can be contributed to by a variety of processes (involving star forming regions, accretion disk coronas and winds, and jets), the powering of the radio loudest quasars must involve very strong jets, presumably launched by the Blandford-Znajek mechanism incorporating the magnetically arrested disk (MAD) scenario. We focus on the latter and investigate the dependence… ▽ More While radio emission in quasars can be contributed to by a variety of processes (involving star forming regions, accretion disk coronas and winds, and jets), the powering of the radio loudest quasars must involve very strong jets, presumably launched by the Blandford-Znajek mechanism incorporating the magnetically arrested disk (MAD) scenario. We focus on the latter and investigate the dependence of their fraction on redshift. We also examine the dependence of the radio-loud fraction (RLF) on BH mass ($M_{\rm BH}$) and Eddington ratio ($λ_{\rm Edd}$) while excluding the redshift bias by narrowing its range. In both these investigations we remove the bias associated with: (1) the diversity of source selection by constructing two well-defined, homogeneous samples of quasars (first within $0.7 \leq z < 1.9$, second within $0.5 \leq z < 0.7$); (2) a strong drop in the RLF of quasars at smaller BH masses by choosing those with BH masses larger than $10^{8.5} M_{\odot}$. We confirm some previous results showing the increase in the fraction of radio-loud quasars with cosmic time and that this trend can be even steeper if we account for the bias introduced by the dependence of the RLF on BH mass whereas the bias introduced by the dependence of the RLF on Eddington ratio is shown to be negligible. Assuming that quasar activities are triggered by galaxy mergers we argue that such an increase can result from the slower drop with cosmic time of mixed mergers than of wet mergers. △ Less

Submitted 15 September, 2021; originally announced September 2021.

Comments: 15 pages, 4 figures. Accepted for publication in ApJ

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

Authors: Andrzej A. Zdziarski, Alexandra J. Tetarenko, Marek Sikora

Abstract: 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… ▽ More 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. △ Less

Submitted 10 November, 2021; v1 submitted 24 August, 2021; originally announced August 2021.

Comments: ApJ, in press

Identifying radio active galactic nuclei among radio-emitting galaxies

Authors: Dorota Kozieł-Wierzbowska, Natalia Vale Asari, Grażyna Stasińska, Fabio R. Herpich, Marek Sikora, Natalia Żywucka, Arti Goyal

Abstract: Basing our analysis on ROGUE I, a catalog of over 32,000 radio sources associated with optical galaxies, we provide two diagnostics to select the galaxies where the radio emission is due to an active galactic nucleus (AGN). Each of these diagnostics can be applied independently. The first one, dubbed MIRAD, compares the flux $F_{W3}$ in the $W3$ mid-infrared band of the WISE telescope, with the ra… ▽ More Basing our analysis on ROGUE I, a catalog of over 32,000 radio sources associated with optical galaxies, we provide two diagnostics to select the galaxies where the radio emission is due to an active galactic nucleus (AGN). Each of these diagnostics can be applied independently. The first one, dubbed MIRAD, compares the flux $F_{W3}$ in the $W3$ mid-infrared band of the WISE telescope, with the radio flux at 1.4 GHz, $\Frad$. MIRAD requires no optical spectra. The second diagnostic, dubbed DLM, relates the 4000 Å break strength, $D_{\rm n}(4000)$, with the radio luminosity per unit stellar mass. The DLM diagram has already been used in the past, but not as standalone. For these two diagrams, we propose simple, empirical dividing lines that result in the same classification for the objects in common. These lines correctly classify as radio-AGN 99.5 percent of the extended radio sources in the ROGUE~I catalog, and as star-forming (SF) galaxies 98--99 percent of the galaxies identified as such by their emission line ratios. Both diagrams clearly show that radio AGNs are preferentially found among elliptical galaxies and among galaxies hosting the most massive black holes. Most of the radio sources classified as radio-AGNs in the MIRAD or DLM diagrams are either optically weak AGNs or retired galaxies. △ Less

Submitted 1 December, 2020; originally announced December 2020.

Comments: 14 pages, 14 figures, submitted to ApJ

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

Authors: Katarzyna Rusinek, Marek Sikora, Dorota Kozieł-Wierzbowska, Maitrayee Gupta

Abstract: 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),… ▽ More 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. △ Less

Submitted 30 December, 2020; v1 submitted 19 June, 2020; originally announced June 2020.

Comments: 18 pages, 4 figures. Published in ApJ on 2020 September 8

Journal ref: The Astrophysical Journal, 2020, Volume 900, Issue 2, id.125, 13 pp

On the significance of relativistically hot pairs in the jets of FR II radio galaxies

Authors: Marek Sikora, Krzysztof Nalewajko, Greg M. Madejski

Abstract: The energetic composition of radio lobes in the FR II galaxies $-$ estimated by comparing their radio luminosities with the powers required to inflate cavities in the external medium $-$ seems to exclude the possibility of their energetic domination by protons. Furthermore, if the jets were dominated by the kinetic energy of cold protons, it would be difficult to efficiently accelerate leptons in… ▽ More The energetic composition of radio lobes in the FR II galaxies $-$ estimated by comparing their radio luminosities with the powers required to inflate cavities in the external medium $-$ seems to exclude the possibility of their energetic domination by protons. Furthermore, if the jets were dominated by the kinetic energy of cold protons, it would be difficult to efficiently accelerate leptons in the jets' terminal shocks. Assuming that the relative energy contents of leptons, protons and magnetic fields are preserved across the shocks, the above implies that the large-scale jets should also be energetically dominated by leptons: $P_{\rm e,j} \gtrsim P_{\rm p,j}$. On the other hand, previous studies of small-scale jets in blazars and radio cores suggest a pair content (number of electrons and positrons per proton) of the order of $n_{\rm e}/n_{\rm p} \sim 20$. Assuming further that the particle composition of jets does not evolve beyond the blazar scales, we show that this implies an average random Lorentz factor of leptons in large-scale jets of $\barγ_{\rm e,j} \gtrsim 70(1+χ_{\rm p})(20n_{\rm p}/n_{\rm e})$, and that the protons should be mildly relativistic with $χ_{\rm p} \equiv (ε_{\rm p} + p_{\rm p})/ρ_{\rm p} c^2 \lesssim 2$, $p_{\rm p}$ the pressure of protons, $ε_{\rm p}$ the internal energy density of protons, and $ρ_{\rm p} c^2$ the rest-mass energy density of protons. We derive the necessary conditions for loading the inner jets by electron-positron pairs and proton-electron plasma, and provide arguments that heating of leptons in jets is dominated by magnetic reconnection. △ Less

Submitted 25 October, 2020; v1 submitted 7 April, 2020; originally announced April 2020.

Comments: 7 pages, accepted for publication in the MNRAS, updated list of references

Orientation of the crescent image of M 87*

Authors: Krzysztof Nalewajko, Marek Sikora, Agata Różańska

Abstract: 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 f… ▽ More 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. △ Less

Submitted 26 February, 2021; v1 submitted 27 August, 2019; originally announced August 2019.

Comments: 7 pages, 7 figures, published in Astronomy & Astrophysics, matching the publication title

Journal ref: A&A 634, A38 (2020)

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

Authors: Maitrayee Gupta, Marek Sikora, Katarzyna Rusinek

Abstract: 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}$.… ▽ More 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. △ Less

Submitted 14 November, 2019; v1 submitted 9 April, 2019; originally announced April 2019.

Comments: Submitted to MNRAS, 11 pages, 14 figures, Comments are welcome

The effect of adiabatic losses on spectra of stationary jets and the origin of soft radio spectra of accreting black-hole sources

Authors: Andrzej A. Zdziarski, Lukasz Stawarz, Marek Sikora

Abstract: It has been suggested that adiabatic energy losses are not effective in stationary jets, where the jet expansion is not associated with net work. Here, we study jet solutions without them, assuming that adiabatic losses are balanced by electron reacceleration. The absence of effective adiabatic losses makes electron advection along the jet an important process, and we solve the electron kinetic eq… ▽ More It has been suggested that adiabatic energy losses are not effective in stationary jets, where the jet expansion is not associated with net work. Here, we study jet solutions without them, assuming that adiabatic losses are balanced by electron reacceleration. The absence of effective adiabatic losses makes electron advection along the jet an important process, and we solve the electron kinetic equation including that process. We find analytical solutions for the case of conical jets with advection and synchrotron losses. We show that accounting for adiabatic losses in the case of sources showing soft partially self-absorbed spectra with the spectral index of $α<0$ in the radio-to-IR regime requires deposition of large amounts of energy at large distances in the jet. On the other hand, such spectra can be accounted for by advection of electrons in the jet. We compare our results to the quiescent spectrum of the blazar Mrk 421. We find its soft radio-IR spectrum can be fitted either by a model without adiabatic losses and advection of electrons or by one with adiabatic losses, but the latter requires injection of a very large power at large distances. △ Less

Submitted 13 February, 2019; v1 submitted 29 December, 2018; originally announced December 2018.

Comments: MNRAS in press, 10 pages

Comparison of hard X-ray spectra of luminous radio galaxies and their radio-quiet counterparts

Authors: Maitrayee Gupta, Marek Sikora, Katarzyna Rusinek, Greg M. Madejski

Abstract: We study the differences in X-ray properties of luminous radio galaxies and their radio-quiet counterparts. In order to avoid biases associated with the dependence of X-ray properties on the black hole mass and Eddington ratio, the radio-loud and radio-quiet objects are selected to cover similar ranges of these parameters. Our studies are based on the X-ray data from the Swift/BAT catalogue. We co… ▽ More We study the differences in X-ray properties of luminous radio galaxies and their radio-quiet counterparts. In order to avoid biases associated with the dependence of X-ray properties on the black hole mass and Eddington ratio, the radio-loud and radio-quiet objects are selected to cover similar ranges of these parameters. Our studies are based on the X-ray data from the Swift/BAT catalogue. We confirm previous results that radio galaxies are on average X-ray-louder than radio-quiet AGNs, but find that their spectral slopes are very similar. This suggests that in radio-loud and radio-quiet AGNs the hard X-rays are produced in the same region and by the same mechanism. We argue that this region can be associated with the hot, geometrically thick, central portion of the accretion flow, where production of hard X-rays is likely to be dominated by Comptonization of the optical/UV radiation of the truncated 'cold' accretion disc by hot electrons. The larger X-ray luminosities in radio-loud AGNs may result from larger radiative efficiencies of the innermost portions of the accretion flows around faster rotating black holes. △ Less

Submitted 21 August, 2018; originally announced August 2018.

Comments: 12 pages, 10 figures, 3 tables; MNRAS, in press

Pair-matching of radio-loud and radio-quiet AGNs

Authors: D. Kozieł-Wierzbowska, G. Stasińska, N. Vale Asari, M. Sikora, E. Goettems, A. Wójtowicz

Abstract: Active galactic nuclei (AGNs) are known to cover an extremely broad range of radio luminosities and the spread of their radio-loudness is very large at any value of the Eddington ratio. This implies very diverse jet production efficiencies which can result from the spread of the black hole spins and magnetic fluxes. Magnetic fluxes can be developed stochastically in the innermost zones of accretio… ▽ More Active galactic nuclei (AGNs) are known to cover an extremely broad range of radio luminosities and the spread of their radio-loudness is very large at any value of the Eddington ratio. This implies very diverse jet production efficiencies which can result from the spread of the black hole spins and magnetic fluxes. Magnetic fluxes can be developed stochastically in the innermost zones of accretion discs, or can be advected to the central regions prior to the AGN phase. In the latter case there could be systematic differences between the properties of galaxies hosting radio-loud (RL) and radio-quiet (RQ) AGNs. In the former case the differences should be negligible for objects having the same Eddington ratio. To study the problem we decided to conduct a comparison study of host galaxy properties of RL and RQ AGNs. In this study we selected type II AGNs from SDSS spectroscopic catalogues. Our RL AGN sample consists of the AGNs appearing in the Best & Heckman (2012) catalogue of radio galaxies. To compare RL and RQ galaxies that have the same AGN parameters we matched the galaxies in black hole mass, Eddington ratio and redshift. We compared several properties of the host galaxies in these two groups of objects like galaxy mass, colour, concentration index, line widths, morphological type and interaction signatures. We found that in the studied group RL AGNs are preferentially hosted by elliptical galaxies while RQ ones are hosted by galaxies of later type. We also found that the fraction of interacting galaxies is the same in both groups of AGNs. These results suggest that the magnetic flux in RL AGNs is advected to the nucleus prior to the AGN phase. △ Less

Submitted 5 November, 2017; v1 submitted 31 October, 2017; originally announced November 2017.

Comments: Fixed small typo in abstract. 7 pages, 4 figures. Front. Astron. Space Sci. - Milky Way and Galaxies in press. Proceedings for "Quasars at all cosmic epochs" meeting, held in Padua, Italy, 2-7 April 2017

What distinguishes the host galaxies of radio-loud and radio-quiet AGNs?

Authors: D. Kozieł-Wierzbowska, N. Vale Asari, G. Stasińska, M. Sikora, E. I. Goettems, A. Wójtowicz

Abstract: We compare the optical properties of the host galaxies of radio-quiet (RQ) and radio-loud (RL) Type 2 active galactic nuclei (AGNs) to infer whether the jet production efficiency depends on the host properties or is determined just by intrinsic properties of the accretion flows. We carefully select galaxies from SDSS, FIRST, and NVSS catalogs. We confirm previous findings that the fraction of RL A… ▽ More We compare the optical properties of the host galaxies of radio-quiet (RQ) and radio-loud (RL) Type 2 active galactic nuclei (AGNs) to infer whether the jet production efficiency depends on the host properties or is determined just by intrinsic properties of the accretion flows. We carefully select galaxies from SDSS, FIRST, and NVSS catalogs. We confirm previous findings that the fraction of RL AGNs depends on the black-hole (BH) masses, and on the Eddington ratio. The comparison of the nature of the hosts of RL and RQ AGNs, therefore, requires pair-matching techniques. Matching in BH mass and Eddington ratio allows us to study the differences between galaxies hosting RL and RQ AGNs that have the same basic accretion parameters. We show that these two samples differ predominantly in the host-galaxy concentration index, morphological type (in the RL sample the frequency of elliptical galaxies becoming larger with increasing radio loudness), and nebular extinction (galaxies with highest radio loudness showing only low nebular extinction). Contrary to some previous studies, we find no significant difference between our radio-loud and radio-quiet samples regarding merger/interaction features. △ Less

Submitted 28 September, 2017; originally announced September 2017.

Comments: 16 pages, 13 figures, published in ApJ

Journal ref: ApJ, 2017, 846, 42

Disk-Jet Connection in Active Supermassive Black Holes in the Standard Accretion Disk Regime

Authors: Yoshiyuki Inoue, Akihiro Doi, Yasuyuki T. Tanaka, Marek Sikora, Gregorz M. Madejski

Abstract: We study the disk-jet connection in supermassive black holes by investigating the properties of their optical and radio emissions utilizing the SDSS-DR7 and the NVSS catalogs. Our sample contains 7017 radio-loud quasars with detection both at 1.4~GHz and SDSS optical spectrum. Using this radio-loud quasar sample, we investigate the correlation among the jet power ($P_{\rm jet}$), the bolometric di… ▽ More We study the disk-jet connection in supermassive black holes by investigating the properties of their optical and radio emissions utilizing the SDSS-DR7 and the NVSS catalogs. Our sample contains 7017 radio-loud quasars with detection both at 1.4~GHz and SDSS optical spectrum. Using this radio-loud quasar sample, we investigate the correlation among the jet power ($P_{\rm jet}$), the bolometric disk luminosity ($L_{\rm disk}$), and the black hole mass ($M_{\rm BH}$) in the standard accretion disk regime. We find that the jet powers correlate with the bolometric disk luminosities as $\log P_{\rm jet} = (0.96\pm0.012)\log L_{\rm disk} + (0.79 \pm 0.55)$. This suggests that the jet production efficiency of $η_{\rm jet}\simeq1.1_{-0.76}^{+2.6}\times10^{-2}$ assuming the disk radiative efficiency of $0.1$ implying low black hole spin parameters and/or low magnetic flux for radio-loud quasars. But it can be also due to dependence of the efficiency on geometrical thickness of the accretion flow which is expected to be small for quasars accreting at the disk Eddington ratios $0.01 \lesssim λ\lesssim 0.3$. This low jet production efficiency does not significantly increase even if we set the disk radiative efficiency of 0.3. We also investigate the fundamental plane in our samples among $P_{\rm jet}$, $L_{\rm disk}$, and $M_{\rm BH}$. We could not find a statistically significant fundamental plane for radio-loud quasars in the standard accretion regime. △ Less

Submitted 1 April, 2017; originally announced April 2017.

Comments: 14 pages, 14 figures, accepted for publication in ApJ

On the efficiency of jet production in FR II radio galaxies and quasars

Authors: Katarzyna Rusinek, Marek Sikora, Dorota Kozieł-Wierzbowska, Leith Godfrey

Abstract: Jet powers in many radio galaxies with extended radio structures appear to exceed their associated accretion luminosities. In systems with very low accretion rates, this is likely due to the very low accretion luminosities resulting from radiatively inefficient accretion flows. In systems with high accretion rates, the accretion flows are expected to be radiatively efficient, and the production of… ▽ More Jet powers in many radio galaxies with extended radio structures appear to exceed their associated accretion luminosities. In systems with very low accretion rates, this is likely due to the very low accretion luminosities resulting from radiatively inefficient accretion flows. In systems with high accretion rates, the accretion flows are expected to be radiatively efficient, and the production of such powerful jets may require an accretion scenario which involves magnetically arrested discs (MADs). However, numerical simulations of the MAD scenario indicate that jet production efficiency is large only for geometrically thick accretion flows and scales roughly with $(H/R)^2$, where $H$ is the disc height and $R$ is the distance from the BH. Using samples of FRII radio galaxies and quasars accreting at moderate accretion rates we show that their jets are much more powerful than predicted by the MAD scenario. We discuss possible origins of this discrepancy, suggesting that it can be related to approximations adopted in MHD simulations to treat optically thick accretion flow within the MAD-zone, or may indicate that accretion disks are geometrically thicker than the standard theory predicts. △ Less

Submitted 21 December, 2016; originally announced December 2016.

Comments: 9 pages, 5 figures, 5 tables. Accepted for publication in MNRAS

Powers and Magnetization of Blazar Jets

Authors: Marek Sikora

Abstract: 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… ▽ More 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. △ Less

Submitted 13 September, 2016; originally announced September 2016.

Comments: 8 pages, 1 figure, accepted for publication in the special issue of Galaxies, "Blazars through Sharp Multi-wavelength Eyes", edited by J. L. Gómez, A. P. Marscher, and S. G. Jorstad

First NuSTAR observations of the BL Lac - type blazar PKS~2155-304: constraints on the jet content and distribution of radiating particles

Authors: G. M. Madejski, K. Nalewajko, K. K. Madsen, J. Chiang, M. Baloković, D. Paneque, A. K. Furniss, M. Hayashida, C. M. Urry, M. Sikora, M. Ajello, R. D. Blandford, F. A. Harrison, D. Sanchez, B. Giebels, D. Stern, D. M. Alexander, D. Barret, S. E. Boggs, F. E. Christensen, W. W. Craig, K. Forster, P. Giommi, B. Grefenstette, C. Hailey , et al. (13 additional authors not shown)

Abstract: We report the first hard X-ray observations with NuSTAR of the BL Lac type blazar PKS 2155-304, augmented with soft X-ray data from XMM-Newton and gamma-ray data from the Fermi Large Area Telescope, obtained in April 2013 when the source was in a very low flux state. A joint NuSTAR and XMM spectrum, covering the energy range 0.5 - 60 keV, is best described by a model consisting of a log-parabola c… ▽ More We report the first hard X-ray observations with NuSTAR of the BL Lac type blazar PKS 2155-304, augmented with soft X-ray data from XMM-Newton and gamma-ray data from the Fermi Large Area Telescope, obtained in April 2013 when the source was in a very low flux state. A joint NuSTAR and XMM spectrum, covering the energy range 0.5 - 60 keV, is best described by a model consisting of a log-parabola component with curvature beta = 0.3(+0.2,-0.1) and a (local) photon index 3.04 +/- 0.15 at photon energy of 2 keV, and a hard power-law tail with photon index 2.2 +/- 0.4. The hard X-ray tail can be smoothly joined to the quasi-simultaneous gamma-ray spectrum by a synchrotron self-Compton component produced by an electron distribution with index p = 2.2. Assuming that the power-law electron distribution extends down to the minimum electron Lorentz factor gamma_min = 1 and that there is one proton per electron, an unrealistically high total jet power L_p of roughly 10^47 erg/s is inferred. This can be reduced by two orders of magnitude either by considering a significant presence of electron-positron pairs with lepton-to-proton ratio of at least 30, or by introducing an additional, low-energy break in the electron energy distribution at the electron Lorentz factor gamma_br1 of roughly 100. In either case, the jet composition is expected to be strongly matter-dominated. △ Less

Submitted 7 September, 2016; originally announced September 2016.

Comments: Astrophysical Journal, in press

The power and production efficiency of blazar jets

Authors: Patryk Pjanka, Andrzej A. Zdziarski, Marek Sikora

Abstract: 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 rad… ▽ More 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. △ Less

Submitted 14 November, 2016; v1 submitted 29 July, 2016; originally announced July 2016.

Comments: Accepted for publication in MNRAS

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

Authors: Maitrayee Gupta, Marek Sikora, Krzysztof Nalewajko

Abstract: 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 matc… ▽ More 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. △ Less

Submitted 13 August, 2016; v1 submitted 11 March, 2016; originally announced March 2016.

Comments: Published in MNRAS, comments are welcome

A spine-sheath model for strong-line blazars

Authors: Marek Sikora, Mieszko Rutkowski, Mitchell C. Begelman

Abstract: 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 th… ▽ More 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. △ Less

Submitted 26 January, 2016; v1 submitted 28 November, 2015; originally announced November 2015.

Comments: 7 pages, 9 figures. Minor changes are made to match the version accepted for publication in MNRAS

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

Authors: Mateusz Janiak, Marek Sikora, Rafal Moderski

Abstract: 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 outbur… ▽ More 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. △ Less

Submitted 26 August, 2015; originally announced August 2015.

Comments: 12 pages, 8 figures, 2 tables

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

Authors: Shinya Saito, Lukasz Stawarz, Yasuyuki Tanaka, Tadayuki Takahashi, Marek Sikora, Rafal Moderski

Abstract: 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 fre… ▽ More 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. △ Less

Submitted 9 July, 2015; originally announced July 2015.

Comments: 11 pages, 7 figures, accepted for publication in ApJ

Rapid Variability of Blazar 3C 279 during Flaring States in 2013-2014 with Joint Fermi-LAT, NuSTAR, Swift, and Ground-Based Multi-wavelength Observations

Authors: M. Hayashida, K. Nalewajko, G. M. Madejski, M. Sikora, R. Itoh, M. Ajello, R. D. Blandford, S. Buson, J. Chiang, Y. Fukazawa, A. K. Furniss, C. M. Urry, I. Hasan, F. A. Harrison, D. M. Alexander, M. Baloković, D. Barret, S. E. Boggs, F. E. Christensen, W. W. Craig, K. Forster, P. Giommi, B. Grefenstette, C. Hailey, A. Hornstrup , et al. (27 additional authors not shown)

Abstract: We report the results of a multi-band observing campaign on the famous blazar 3C 279 conducted during a phase of increased activity from 2013 December to 2014 April, including first observations of it with NuSTAR. The $γ$-ray emission of the source measured by Fermi-LAT showed multiple distinct flares reaching the highest flux level measured in this object since the beginning of the Fermi mission,… ▽ More We report the results of a multi-band observing campaign on the famous blazar 3C 279 conducted during a phase of increased activity from 2013 December to 2014 April, including first observations of it with NuSTAR. The $γ$-ray emission of the source measured by Fermi-LAT showed multiple distinct flares reaching the highest flux level measured in this object since the beginning of the Fermi mission, with $F(E > 100\,{\rm MeV})$ of $10^{-5}$ photons cm$^{-2}$ s$^{-1}$, and with a flux doubling time scale as short as 2 hours. The $γ$-ray spectrum during one of the flares was very hard, with an index of $Γ_γ= 1.7 \pm 0.1$, which is rarely seen in flat spectrum radio quasars. The lack of concurrent optical variability implies a very high Compton dominance parameter $L_γ/L_{\rm syn} > 300$. Two 1-day NuSTAR observations with accompanying Swift pointings were separated by 2 weeks, probing different levels of source activity. While the 0.5$-$70 keV X-ray spectrum obtained during the first pointing, and fitted jointly with Swift-XRT is well-described by a simple power law, the second joint observation showed an unusual spectral structure: the spectrum softens by $ΔΓ_{\rm X} \simeq 0.4$ at $\sim$4 keV. Modeling the broad-band SED during this flare with the standard synchrotron plus inverse Compton model requires: (1) the location of the $γ$-ray emitting region is comparable with the broad line region radius, (2) a very hard electron energy distribution index $p \simeq 1$, (3) total jet power significantly exceeding the accretion disk luminosity $L_{\rm j}/L_{\rm d} \gtrsim 10$, and (4) extremely low jet magnetization with $L_{\rm B}/L_{\rm j} \lesssim 10^{-4}$. We also find that single-zone models that match the observed $γ$-ray and optical spectra cannot satisfactorily explain the production of X-ray emission. △ Less

Submitted 2 April, 2015; v1 submitted 16 February, 2015; originally announced February 2015.

Comments: 17 pages, 10 figures, accepted for publication in The Astrophysical Journal

Magnetization of jets in luminous blazars

Authors: Mateusz Janiak, Marek Sikora, Rafał Moderski

Abstract: 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 t… ▽ More 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. △ Less

Submitted 18 March, 2015; v1 submitted 26 November, 2014; originally announced November 2014.

Comments: 19 pages, 6 figures, published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society 2015 449 (4): 431-439

Core shifts, magnetic fields and magnetization of extragalactic jets

Authors: Andrzej A. Zdziarski, Marek Sikora, Patryk Pjanka, Alexander Tchekhovskoy

Abstract: 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… ▽ More 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. △ Less

Submitted 2 May, 2015; v1 submitted 27 October, 2014; originally announced October 2014.

Comments: MNRAS, in press

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

Authors: Krzysztof Nalewajko, Marek Sikora, Mitchell C. Begelman

Abstract: 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 synchrot… ▽ More 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. △ Less

Submitted 16 October, 2014; originally announced October 2014.

Comments: 5 pages, accepted in ApJ Letters

Constraining the Location of Gamma-Ray Flares in Luminous Blazars

Authors: Krzysztof Nalewajko, Mitchell C. Begelman, Marek Sikora

Abstract: 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 para… ▽ More 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. △ Less

Submitted 29 May, 2014; originally announced May 2014.

Comments: 21 pages, 10 figures, accepted for publication in ApJ

Peculiar Radio Structures in the Central Regions of Galaxy Cluster Abell 585

Authors: M. Jamrozy, Ł Stawarz, V. Marchenko, A. Kuzmicz, M. Ostrowski, C. C. Cheung, M. Sikora

Abstract: In this paper we analyze the peculiar radio structure observed across the central region of the galaxy cluster Abell 585 (z=0.12). In the low-resolution radio maps, this structure appears uniform and diffuse on angular scales of ~3 arcmin, and is seemingly related to the distant (z=2.5) radio quasar B3 0727+409 rather than to the cluster itself. However, after a careful investigation of the unpubl… ▽ More In this paper we analyze the peculiar radio structure observed across the central region of the galaxy cluster Abell 585 (z=0.12). In the low-resolution radio maps, this structure appears uniform and diffuse on angular scales of ~3 arcmin, and is seemingly related to the distant (z=2.5) radio quasar B3 0727+409 rather than to the cluster itself. However, after a careful investigation of the unpublished archival radio data with better angular resolution, we resolve the structure into two distinct arcmin-scale features, which resemble typical lobes of cluster radio galaxies with no obvious connection to the background quasar. We support this conclusion by examining the spectral and polarization properties of the features, demonstrating in addition that the analyzed structure can hardly be associated with any sort of a radio mini-halo or relics of the cluster. Yet at the same time we are not able to identify host galaxies of the radio lobes in the available optical and infrared surveys. We consider some speculative explanations for our findings, including gravitational wave recoil kicks of SMBHs responsible for the lobes' formation in the process of merging massive ellipticals within the central parts of a rich cluster environment, but we do not reach any robust conclusions regarding the origin of the detected radio features. △ Less

Submitted 1 April, 2014; originally announced April 2014.

Comments: 11 pages, 3 figures, 2 tables, accepted for publication in MNRAS

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

Authors: Andrzej A. Zdziarski, Patryk Pjanka, Marek Sikora, Lukasz Stawarz

Abstract: 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 emiss… ▽ More 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. △ Less

Submitted 25 May, 2014; v1 submitted 19 March, 2014; originally announced March 2014.

Comments: MNRAS, in press

On the origin of X-ray spectra in luminous blazars

Authors: Marek Sikora, Mateusz Janiak, Krzysztof Nalewajko, Greg M. Madejski, Rafał Moderski

Abstract: 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-… ▽ More 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. △ Less

Submitted 4 November, 2013; originally announced November 2013.

Comments: 19 pages, 1 figure, accepted for publication in ApJ

CTA contributions to the 33rd International Cosmic Ray Conference (ICRC2013)

Authors: The CTA Consortium, :, O. Abril, B. S. Acharya, M. Actis, G. Agnetta, J. A. Aguilar, F. Aharonian, M. Ajello, A. Akhperjanian, M. Alcubierre, J. Aleksic, R. Alfaro, E. Aliu, A. J. Allafort, D. Allan, I. Allekotte, R. Aloisio, E. Amato, G. Ambrosi, M. Ambrosio, J. Anderson, E. O. Angüner, L. A. Antonelli, V. Antonuccio , et al. (1082 additional authors not shown)

Abstract: Compilation of CTA contributions to the proceedings of the 33rd International Cosmic Ray Conference (ICRC2013), which took place in 2-9 July, 2013, in Rio de Janeiro, Brazil Compilation of CTA contributions to the proceedings of the 33rd International Cosmic Ray Conference (ICRC2013), which took place in 2-9 July, 2013, in Rio de Janeiro, Brazil △ Less

Submitted 29 July, 2013; v1 submitted 8 July, 2013; originally announced July 2013.

Comments: Index of CTA conference proceedings at the ICRC2013, Rio de Janeiro (Brazil). v1: placeholder with no arXiv links yet, to be replaced once individual contributions have been all submitted. v2: final with arXiv links to all CTA contributions and full author list

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

Authors: Andrzej A. Zdziarski, Lukasz Stawarz, Patryk Pjanka, Marek Sikora

Abstract: 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 stat… ▽ More 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. △ Less

Submitted 3 March, 2014; v1 submitted 4 July, 2013; originally announced July 2013.

Comments: MNRAS in press; the original submission has been split into two, this is part I

Magnetic flux paradigm for radio-loudness of AGN

Authors: Marek Sikora, Mitchell C. Begelman

Abstract: 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 lik… ▽ More 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. △ Less

Submitted 23 January, 2013; originally announced January 2013.

Comments: 5 pages, no figures, accepted for publication the The Astrophysical Journal Letters

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

Authors: Mateusz Janiak, Marek Sikora, Krzysztof Nalewajko, Rafał Moderski, Greg M. Madejski

Abstract: 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 d… ▽ More 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. △ Less

Submitted 18 October, 2012; originally announced October 2012.

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

Authors: Krzysztof Nalewajko, Marek Sikora, Greg M. Madejski, Katrina Exter, Anna Szostek, Ryszard Szczerba, Mark R. Kidger, Rosario Lorente

Abstract: 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-… ▽ More 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. △ Less

Submitted 16 October, 2012; originally announced October 2012.

Comments: 13 pages, 8 figures, 7 tables; accepted for publication in the Astrophysical Journal

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

Authors: Marek Sikora, Grażyna Stasińska, Dorota Kozieł-Wierzbowska, Greg M. Madejski, Natalia V. Asari

Abstract: 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 abo… ▽ More 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. △ Less

Submitted 18 February, 2013; v1 submitted 9 October, 2012; originally announced October 2012.

Comments: 24 pages, 6 figures, published in ApJ, moderate revisions to match the published version

Journal ref: Astrophysical Journal 765 (2013) 62-70

The evolution of active galactic nuclei and their spins

Authors: Marta Volonteri, Marek Sikora, Jean-Pierre Lasota, Andrea Merloni

Abstract: 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 j… ▽ More 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. △ Less

Submitted 15 August, 2013; v1 submitted 3 October, 2012; originally announced October 2012.

Comments: ApJ, in press

Multi-wavelength observations of blazar AO 0235+164 in the 2008-2009 flaring state

Authors: M. Ackermann, M. Ajello, J. Ballet, G. Barbiellini, D. Bastieri, R. Bellazzini, R. D. Blandford, E. D. Bloom, E. Bonamente, A. W. Borgland, E. Bottacini, J. Bregeon, M. Brigida, P. Bruel, R. Buehler, S. Buson, G. A. Caliandro, R. A. Cameron, P. A. Caraveo, J. M. Casandjian, E. Cavazzuti, C. Cecchi, E. Charles, A. Chekhtman, J. Chiang , et al. (186 additional authors not shown)

Abstract: The blazar AO 0235+164 (z = 0.94) has been one of the most active objects observed by Fermi Large Area Telescope (LAT) since its launch in Summer 2008. In addition to the continuous coverage by Fermi, contemporaneous observations were carried out from the radio to γ -ray bands between 2008 September and 2009 February. In this paper, we summarize the rich multi-wavelength data collected during the… ▽ More The blazar AO 0235+164 (z = 0.94) has been one of the most active objects observed by Fermi Large Area Telescope (LAT) since its launch in Summer 2008. In addition to the continuous coverage by Fermi, contemporaneous observations were carried out from the radio to γ -ray bands between 2008 September and 2009 February. In this paper, we summarize the rich multi-wavelength data collected during the campaign (including F-GAMMA, GASP- WEBT, Kanata, OVRO, RXTE, SMARTS, Swift, and other instruments), examine the cross-correlation between the light curves measured in the different energy bands, and interpret the resulting spectral energy distributions in the context of well-known blazar emission models. We find that the γ -ray activity is well correlated with a series of near-IR/optical flares, accompanied by an increase in the optical polarization degree. On the other hand, the X-ray light curve shows a distinct 20 day high state of unusually soft spectrum, which does not match the extrapolation of the optical/UV synchrotron spectrum. We tentatively interpret this feature as the bulk Compton emission by cold electrons contained in the jet, which requires an accretion disk corona with an effective covering factor of 19% at a distance of 100 Rg . We model the broadband spectra with a leptonic model with external radiation dominated by the infrared emission from the dusty torus. △ Less

Submitted 12 July, 2012; originally announced July 2012.

Comments: 48 pages, 8 figures

Journal ref: ApJ 751, 159 (2012)

Quasi-star jets as unidentified gamma-ray sources

Authors: Bozena Czerny, Agnieszka Janiuk, Marek Sikora, Jean-Pierre Lasota

Abstract: 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 synchro… ▽ More 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. △ Less

Submitted 6 July, 2012; originally announced July 2012.

Comments: 5 pages, 1 figure; accepted by ApJ Letters

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

Authors: Krzysztof Nalewajko, Marek Sikora

Abstract: 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 polarizatio… ▽ More 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. △ Less

Submitted 28 June, 2012; originally announced June 2012.

Comments: 10 pages, 6 figures, accepted for publication in A&A

The structure and emission model of the relativistic jet in the quasar 3C 279 inferred from radio to high-energy gamma-ray observations in 2008-2010

Authors: M. Hayashida, G. M. Madejski, K. Nalewajko, M. Sikora, A. E. Wehrle, P. Ogle, W. Collmar, S. Larsson, Y. Fukazawa, R. Itoh, J. Chiang, L. Stawarz, R. D. Blandford, J. L. Richards, W. Max-Moerbeck, A. Readhead, R. Buehler, E. Cavazzuti, S. Ciprini, N. Gehrels, A. Reimer, A. Szostek, T. Tanaka, G. Tosti, Y. Uchiyama , et al. (68 additional authors not shown)

Abstract: We present time-resolved broad-band observations of the quasar 3C 279 obtained from multi-wavelength campaigns conducted during the first two years of the Fermi Gamma-ray Space Telescope mission. While investigating the previously reported gamma-ray/optical flare accompanied by a change in optical polarization, we found that the optical emission appears delayed with respect to the gamma-ray emissi… ▽ More We present time-resolved broad-band observations of the quasar 3C 279 obtained from multi-wavelength campaigns conducted during the first two years of the Fermi Gamma-ray Space Telescope mission. While investigating the previously reported gamma-ray/optical flare accompanied by a change in optical polarization, we found that the optical emission appears delayed with respect to the gamma-ray emission by about 10 days. X-ray observations reveal a pair of `isolated' flares separated by ~90 days, with only weak gamma-ray/optical counterparts. The spectral structure measured by Spitzer reveals a synchrotron component peaking in the mid-infrared band with a sharp break at the far-infrared band during the gamma-ray flare, while the peak appears in the mm/sub-mm band in the low state. Selected spectral energy distributions are fitted with leptonic models including Comptonization of external radiation produced in a dusty torus or the broad-line region. Adopting the interpretation of the polarization swing involving propagation of the emitting region along a curved trajectory, we can explain the evolution of the broad-band spectra during the gamma-ray flaring event by a shift of its location from ~ 1 pc to ~ 4 pc from the central black hole. On the other hand, if the gamma-ray flare is generated instead at sub-pc distance from the central black hole, the far-infrared break can be explained by synchrotron self-absorption. We also model the low spectral state, dominated by the mm/sub-mm peaking synchrotron component, and suggest that the corresponding inverse-Compton component explains the steady X-ray emission. △ Less

Submitted 4 June, 2012; originally announced June 2012.

Comments: 23 pages, 18 figures 5 tables, Accepted for publication in The Astrophysical Journal

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

Authors: Andrzej A. Zdziarski, Piotr Lubinski, Marek Sikora

Abstract: 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 a… ▽ More 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. △ Less

Submitted 12 March, 2012; originally announced March 2012.

Comments: MNRAS, in press, 13 pages

Journal ref: Mon. Not. R. Astron. Soc. 423, 663-675 (2012)

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

Authors: Krzysztof Nalewajko, Mitchell C. Begelman, Benoit Cerutti, Dmitri A. Uzdensky, Marek Sikora

Abstract: 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 esti… ▽ More 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. △ Less

Submitted 13 July, 2012; v1 submitted 9 February, 2012; originally announced February 2012.

Comments: 12 pages, 2 figures, accepted for publication in MNRAS

A multiwavelength view of the flaring state of PKS 2155-304 in 2006

Authors: HESS Collaboration, A. Abramowski, F. Acero, F. Aharonian, A. G. Akhperjanian, G. Anton, A. Balzer, A. Barnacka, U. Barres de Almeida, Y. Becherini, J. Becker, B. Behera, W. Benbow, K. Bernlöhr, A. Bochow, C. Boisson, J. Bolmont, P. Bordas, T. Boutelier, J. Brucker, F. Brun, P. Brun, T. Bulik, I. Büsching, S. Carrigan , et al. (186 additional authors not shown)

Abstract: Multiwavelength (MWL) observations of the blazar PKS 2155-304 during two weeks in July and August 2006, the period when two exceptional flares at very high energies (VHE, E>= 100 GeV) occurred, provide a detailed picture of the evolution of its emission. The complete data set from this campaign is presented, including observations in VHE gamma-rays (H.E.S.S.), X-rays (RXTE, CHANDRA, SWIFT XRT), op… ▽ More Multiwavelength (MWL) observations of the blazar PKS 2155-304 during two weeks in July and August 2006, the period when two exceptional flares at very high energies (VHE, E>= 100 GeV) occurred, provide a detailed picture of the evolution of its emission. The complete data set from this campaign is presented, including observations in VHE gamma-rays (H.E.S.S.), X-rays (RXTE, CHANDRA, SWIFT XRT), optical (SWIFT UVOT, Bronberg, Watcher, ROTSE), and in the radio band (NRT, HartRAO, ATCA). Optical and radio light curves from 2004 to 2008 are compared to the available VHE data from this period, to put the 2006 campaign into the context of the long-term evolution of the source. The X-ray and VHE gamma-ray emission are correlated during the observed high state of the source, but show no direct connection with longer wavelengths. The long-term flux evolution in the optical and radio bands is found to be correlated and shows that the source reaches a high state at long wavelengths after the occurrence of the VHE flares. Spectral hardening is seen in the SWIFT XRT data. The nightly averaged high-energy spectra of the non-flaring nights can be reproduced by a stationary one-zone SSC model, with only small variations in the parameters. The spectral and flux evolution in the high-energy band during the night of the second VHE flare is modelled with multi-zone SSC models, which can provide relatively simple interpretations for the hour time-scale evolution of the high-energy emission, even for such a complex data set. For the first time in this type of source, a clear indication is found for a relation between high activity at high energies and a long-term increase in the low frequency fluxes. △ Less

Submitted 19 January, 2012; originally announced January 2012.

Comments: 23 pages, 13 figures, accepted for publication in Astronomy & Astrophysics

Discovery of hard-spectrum γ-ray emission from the BL Lac object 1ES 0414+009

Authors: The HESS Collaboration, A. Abramowski, F. Acero, F. Aharonian, A. G. Akhperjanian, G. Anton, A. Balzer, A. Barnacka, U. Barres de Almeida, Y. Becherini, J. Becker, B. Behera, K. Bernlöhr, E. Birsin, J. Biteau, A. Bochow, C. Boisson, J. Bolmont, P. Bordas, J. Brucker, F. Brun, P. Brun, T. Bulik, I. Büsching, S. Carrigan , et al. (177 additional authors not shown)

Abstract: 1ES 0414+009 (z = 0.287) is a distant high-frequency-peaked BL Lac object, and has long been considered a likely emitter of very-high energy (VHE, E>100 GeV) gamma-rays due to its high X-ray and radio flux. Observations in the VHE gamma-ray band and across the electromagnetic spectrum can provide insights into the origin of highly energetic particles present in the source and the radiation process… ▽ More 1ES 0414+009 (z = 0.287) is a distant high-frequency-peaked BL Lac object, and has long been considered a likely emitter of very-high energy (VHE, E>100 GeV) gamma-rays due to its high X-ray and radio flux. Observations in the VHE gamma-ray band and across the electromagnetic spectrum can provide insights into the origin of highly energetic particles present in the source and the radiation processes at work. Because of the distance of the source, the gamma-ray spectrum might provide further limits on the level of the Extragalactic Background Light (EBL). We report observations made between October 2005 and December 2009 with H.E.S.S., an array of four imaging atmospheric Cherenkov telescopes. Observations at high energies (HE, 100 MeV - 100 GeV) with the Fermi-LAT instrument in the first 20 months of its operation are also reported. To complete the multi-wavelength picture, archival UV and X-ray observations with the Swift satellite and optical observations with the ATOM telescope are also used. Based on the observations with H.E.S.S., 1ES 0414+009 is detected for the first time in the VHE band. An excess of 224 events is measured, corresponding to a significance of 7.8 sigma. The photon spectrum of the source is well described by a power law, with photon index of 3.45 \pm 0.25stat \pm 0.20syst. The integral flux above 200 GeV is (1.88 \pm 0.20stat \pm 0.38syst) \times10-12 cm-2 s-1. Observations with the Fermi-LAT in the first 20 months of operation show a flux between 200 MeV and 100 GeV of (2.3 \pm 0.2stat) \times 10-9 erg cm-2 s-1, and a spectrum well described by a power-law function with a photon index 1.85 \pm 0.18. Swift/XRT observations show an X-ray flux between 2 and 10 keV of (0.8 - 1) \times 10-11 erg cm-2 s-1, and a steep spectrum (2.2 - 2.3). Combining X-ray with optical-UV data, a fit with a log-parabolic function locates the synchrotron peak around 0.1 keV. ... △ Less

Submitted 10 January, 2012; originally announced January 2012.

Comments: 9 pages, 4 figures, accepted for publication in A&A

Simultaneous multi-wavelength campaign on PKS 2005-489 in a high state

Authors: The H. E. S. S. Collaboration, the Fermi LAT Collaboration, :, A. Abramowski, F. Acero, F. Aharonian, A. G. Akhperjanian, G. Anton, A. Barnacka, U. Barres de Almeida, A. R. Bazer-Bachi, Y. Becherini, J. Becker, B. Behera, K. Bernlöhr, A. Bochow, C. Boisson, J. Bolmont, P. Bordas, V. Borrel, J. Brucker, F. Brun, P. Brun, T. Bulik, I. Büsching , et al. (323 additional authors not shown)

Abstract: The high-frequency peaked BL Lac object PKS 2005-489 was the target of a multi-wavelength campaign with simultaneous observations in the TeV gamma-ray (H.E.S.S.), GeV gamma-ray (Fermi/LAT), X-ray (RXTE, Swift), UV (Swift) and optical (ATOM, Swift) bands. This campaign was carried out during a high flux state in the synchrotron regime. The flux in the optical and X-ray bands reached the level of th… ▽ More The high-frequency peaked BL Lac object PKS 2005-489 was the target of a multi-wavelength campaign with simultaneous observations in the TeV gamma-ray (H.E.S.S.), GeV gamma-ray (Fermi/LAT), X-ray (RXTE, Swift), UV (Swift) and optical (ATOM, Swift) bands. This campaign was carried out during a high flux state in the synchrotron regime. The flux in the optical and X-ray bands reached the level of the historical maxima. The hard GeV spectrum observed with Fermi/LAT connects well to the very high energy (VHE, E>100GeV) spectrum measured with H.E.S.S. with a peak energy between ~5 and 500 GeV. Compared to observations with contemporaneous coverage in the VHE and X-ray bands in 2004, the X-ray flux was ~50 times higher during the 2009 campaign while the TeV gamma-ray flux shows marginal variation over the years. The spectral energy distribution during this multi-wavelength campaign was fit by a one zone synchrotron self-Compton model with a well determined cutoff in X-rays. The parameters of a one zone SSC model are inconsistent with variability time scales. The variability behaviour over years with the large changes in synchrotron emission and small changes in the inverse Compton emission does not warrant an interpretation within a one-zone SSC model despite an apparently satisfying fit to the broadband data in 2009. △ Less

Submitted 14 November, 2011; originally announced November 2011.

Comments: published in A&A, 9pages, 4figures. Due to some style difficulties using the standard style file with the long author list, the abstract has been included manually. A version with the correct style can be found in the website http://www.lsw.uni-heidelberg.de/users/skaufman/publications/aa16170-10published.pdf

Journal ref: A&A 533, A110 (2011)

Discovery of extended VHE γ-ray emission from the vicinity of the young massive stellar cluster Westerlund 1

Authors: HESS Collaboration, A. Abramowski, F. Acero, F. Aharonian, A. G. Akhperjanian, G. Anton, A. Balzer, A. Barnacka, U. Barres de Almeida, Y. Becherini, J. Becker, B. Behera, K. Bernlöhr, E. Birsin, J. Biteau, A. Bochow, C. Boisson, J. Bolmont, P. Bordas, J. Brucker, F. Brun, P. Brun, T. Bulik, I. Büsching, S. Carrigan , et al. (176 additional authors not shown)

Abstract: Results obtained in very-high-energy (VHE; E > 100 GeV) γ-ray observations performed with the H.E.S.S. telescope array are used to investigate particle acceleration processes in the vicinity of the young massive stellar cluster Westerlund 1 (Wd 1). Imaging of Cherenkov light from γ-ray induced particle cascades in the Earth's atmosphere is used to search for VHE γ rays from the region around Wd 1.… ▽ More Results obtained in very-high-energy (VHE; E > 100 GeV) γ-ray observations performed with the H.E.S.S. telescope array are used to investigate particle acceleration processes in the vicinity of the young massive stellar cluster Westerlund 1 (Wd 1). Imaging of Cherenkov light from γ-ray induced particle cascades in the Earth's atmosphere is used to search for VHE γ rays from the region around Wd 1. Possible catalogued counterparts are searched for and discussed in terms of morphology and energetics of the H.E.S.S. source. The detection of the degree-scale extended VHE γ-ray source HESS J1646-458 is reported based on 45 hours of H.E.S.S. observations performed between 2004 and 2008. The VHE γ-ray source is centred on the nominal position of Wd 1 and detected with a total statistical significance of ~20σ. The emission region clearly extends beyond the H.E.S.S. point-spread function (PSF). The differential energy spectrum follows a power law in energy with an index of Γ=2.19 \pm 0.08_{stat} \pm 0.20_{sys} and a flux normalisation at 1 TeV of Φ_0 = (9.0 \pm 1.4_{stat} \pm 1.8_{sys}) x 10^{-12} TeV^{-1} cm^{-2} s^{-1}. The integral flux above 0.2 TeV amounts to (5.2 \pm 0.9) x 10^{-11} cm^{-2} s^{-1}. Four objects coincident with HESS J1646-458 are discussed in the search of a counterpart, namely the magnetar CXOU J164710.2-455216, the X-ray binary 4U 1642-45, the pulsar PSR J1648-4611 and the massive stellar cluster Wd 1. In a single-source scenario, Wd 1 is favoured as site of VHE particle acceleration. Here, a hadronic parent population would be accelerated within the stellar cluster. Beside this, there is evidence for a multi-source origin, where a scenario involving PSR J1648-4611 could be viable to explain parts of the VHE γ-ray emission of HESS J1646-458. △ Less

Submitted 9 November, 2011; originally announced November 2011.

Comments: 11 pages, 6 figures, accepted for publication in A&A

Journal ref: A&A 537, A114 (2012)