Source author record

Krzysztof Nalewajko

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

29works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

29 published item(s)

preprint2022arXiv

A simple analytical model of magnetic jets

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

preprint2022arXiv

Kinetic Simulations of Instabilities and Particle Acceleration in Cylindrical Magnetized Relativistic Jets

Relativistic magnetized jets, such as those from AGN, GRBs and XRBs, are susceptible to current- and pressure-driven MHD instabilities that can lead to particle acceleration and non-thermal radiation. Here we investigate the development of these instabilities through 3D kinetic simulations of cylindrically symmetric equilibria involving toroidal magnetic fields with electron-positron pair plasma. Generalizing recent treatments by Alves et al. (2018) and Davelaar et al. (2020), we consider a range of initial structures in which the force due to toroidal magnetic field is balanced by a combination of forces due to axial magnetic field and gas pressure. We argue that the particle energy limit identified by Alves et al. (2018) is due to the finite duration of the fast magnetic dissipation phase. We find a rather minor role of electric fields parallel to the local magnetic fields in particle acceleration. In all investigated cases a kink mode arises in the central core region with a growth timescale consistent with the predictions of linearized MHD models. In the case of a gas-pressure-balanced (Z-pinch) profile, we identify a weak local pinch mode well outside the jet core. We argue that pressure-driven modes are important for relativistic jets, in regions where sufficient gas pressure is produced by other dissipation mechanisms.

preprint2022arXiv

Multiwavelength variability power spectrum analysis of the blazars 3C 279 and PKS 1510-089 on multiple timescales

We present the results of variability power spectral density (PSD) analysis using multiwavelength radio to GeV\,$γ$-ray light curves covering decades/years to days/minutes timescales for the blazars 3C 279 and PKS 1510-089. The PSDs are modeled as single power-laws, and the best-fit spectral shape is derived using the `power spectral response' method. With more than ten years of data obtained with weekly/daily sampling intervals, most of the PSDs cover ~2-4 decades in temporal frequency; moreover, in the optical band, the PSDs cover ~6 decades for 3C 279 due to the availability of intranight light curves. Our main results are the following: (1) on timescales ranging from decades to days, the synchrotron and the inverse Compton spectral components, in general, exhibit red-noise (slope ~2) and flicker-noise (slope ~1) type variability, respectively; (2) the slopes of $γ$-ray variability PSDs obtained using a 3-hr integration bin and a 3-weeks total duration exhibit a range between ~1.4 and ~2.0 (mean slope = 1.60$\pm$0.70), consistent within errors with the slope on longer timescales; (3) comparisons of fractional variability indicate more power on timescales $\leq$100\, days at $γ$-ray frequencies as compared to longer wavelengths, in general (except between $γ$-ray and optical frequencies for PKS 1510$-$089); (4) the normalization of intranight optical PSDs for 3C\,279 appears to be a simple extrapolation from longer timescales, indicating a continuous (single) process driving the variability at optical wavelengths; (5) the emission at optical/infrared wavelengths may involve a combination of disk and jet processes for PKS\,1510$-$089.

preprint2022arXiv

Rapid X-ray Variability in Mkn 421 during a Multiwavelength Campaign

The study of short-term variability properties in AGN jets has the potential to shed light on their particle acceleration and emission mechanisms. We report results from a four-day coordinated multi-wavelength campaign on the highly-peaked blazar (HBL) Mkn 421 in 2019 January. We obtained X-ray data from AstroSAT, BVRI photometry with the Whole Earth Blazar Telescope (WEBT), and TeV data from FACT to explore short-term multi-wavelength variability in this HBL. The X-ray continuum is rapidly variable on time-scales of tens of ks. Fractional variability amplitude increases with energy across the synchrotron hump, consistent with previous studies; we interpret this observation in the context of a model with multiple cells whose emission spectra contain cutoffs that follow a power-law distribution. We also performed time-averaged and time-resolved (time-scales of 6 ks) spectral fits; a broken power-law model fits all spectra well; time-resolved spectral fitting reveals the usual hardening when brightening behaviour. Intra-X-ray cross correlations yield evidence for the 0.6-0.8 keV band to likely lead the other bands by an average of 4.6 +- 2.6 ks, but only during the first half of the observation. The source displayed minimal night-to-night variability at all wavebands thus precluding significant interband correlations during our campaign. The broadband SED is modeled well with a standard one-zone leptonic model, yielding jet parameters consistent with those obtained from previous SEDs of this source.

preprint2021arXiv

Orientation of the crescent image of M 87*

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

preprint2020arXiv

Long-term optical spectroscopic variations in blazar 3C 454.3

Characterisation of the long-term variations in the broad line region in a luminous blazar, where Comptonisation of broad-line emission within a relativistic jet is the standard scenario for production of gamma-ray emission that dominates the spectral energy distribution. We analysed ten years of optical spectroscopic data from the Steward Observatory for the blazar 3C 454.3, as well as gamma-ray data from the Fermi Large Area Telescope (LAT). The optical spectra are dominated by a highly variable non-thermal synchrotron continuum with a prominent Mg II broad emission line. The line flux was obtained by spectral decomposition including significant contribution from the Fe II pseudo-continuum. Three methods were used to characterise variations in the line flux: (1) stacking of the continuum-subtracted spectra, (2) subtracting the running mean light curves calculated for different timescales, and (3) evaluating potential time delays via the discrete correlation function (DCF). Despite very large variations in the gamma-ray and optical continua, the line flux changes only moderately (< 0.1 dex). The data suggest that the line flux responds to a dramatic change in the blazar activity from a very high state in 2010 to a deep low state in 2012. Two interpretations are possible: either the line flux is anti-correlated with the continuum or the increase in the line luminosity is delayed by ~600 days. If this time delay results from the reverberation of poorly constrained accretion disc emission in both the broad-line region (BLR) and the synchrotron emitting blazar zone within a relativistic jet, we would obtain natural estimates for the BLR radius [R_{BLR,MgII} >~ 0.28 pc] and for the supermassive black hole mass [M_SMBH ~ 8.5x10^8 M_sun]. We did not identify additional examples of short-term 'flares' of the line flux, in addition to the previously reported case observed in 2010.

preprint2020arXiv

Radiative kinetic simulations of steady-state relativistic plasmoid magnetic reconnection

We present the results of 2D particle-in-cell (PIC) simulations of relativistic magnetic reconnection (RMR) in electron-positron plasma, including the dynamical influence of the synchrotron radiation process, and integrating the observable emission signatures. The simulations are initiated with a single Harris current layer with a central gap that triggers the RMR process. We achieve a steady-state reconnection with unrestricted outflows by means of open boundary conditions. The radiative cooling efficiency is regulated by the choice of initial plasma temperature Theta. We explore different values of Theta and of the background magnetisation sigma_0. Throughout the simulations, plasmoids are generated in the central region of the layer, and they evolve at different rates, achieving a wide range of sizes. The gaps between plasmoids are filled by smooth relativistic outflows called minijets, whose contribution to the observed radiation is very limited due to their low particle densities. Small-sized plasmoids are rapidly accelerated, however, they have lower contributions to the observed emission, despite stronger relativistic beaming. Large-sized plasmoids are slow, but produce most of the observed synchrotron emission, with major part of their radiation produced within the central cores, the density of which is enhanced by radiative cooling. Synchrotron lightcurves show rapid bright flares that can be identified as originating from mergers between small/fast plasmoids and large/slow targets moving in the same direction. In the high-magnetisation case, the accelerated particles form a broken power-law energy distribution with a soft tail produced by particles accelerated in the minijets.

preprint2016arXiv

Applying Relativistic Reconnection to Blazar Jets

Rapid and luminous flares of non-thermal radiation observed in blazars require an efficient mechanism of energy dissipation and particle acceleration in relativistic active galactic nuclei (AGN) jets. Particle acceleration in relativistic magnetic reconnection is being actively studied by kinetic numerical simulations. Relativistic reconnection produces hard power-law electron energy distributions N(gamma) = N_0 gamma^(-p) exp(-gamma/gamma_max) with index p -> 1 and exponential cut-off Lorentz factor gamma_max ~ sigma in the limit of magnetization sigma = B^2/(4 pi w) >> 1 (where w is the relativistic enthalpy density). Reconnection in electron-proton plasma can additionally boost gamma_max by the mass ratio m_p/m_e. Hence, in order to accelerate particles to gamma_max ~ 10^6 in the case of BL Lacs, reconnection should proceed in plasma of very high magnetization sigma_max >~ 10^3. On the other hand, moderate mean jet magnetization values are required for magnetic bulk acceleration of relativistic jets, sigma_mean ~ Gamma_j <~ 20 (where Gamma_j is the jet bulk Lorentz factor). I propose that the systematic dependence of gamma_max on blazar luminosity class -- the blazar sequence -- may result from a systematic trend in sigma_max due to homogeneous loading of leptons by pair creation regulated by the energy density of high-energy external radiation fields. At the same time, relativistic AGN jets should be highly inhomogeneous due to filamentary loading of protons, which should determine the value of sigma_mean roughly independently of the blazar class.

preprint2016arXiv

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

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

preprint2016arXiv

Kinetic simulations of the lowest-order unstable mode of relativistic magnetostatic equilibria

We present the results of particle-in-cell numerical pair plasma simulations of relativistic 2D magnetostatic equilibria known as the 'ABC' fields. In particular, we focus on the lowest-order unstable configuration consisting of two minima and two maxima of the magnetic vector potential. Breaking of the initial symmetry leads to exponential growth of the electric energy and to the formation of two current layers, which is consistent with the picture of 'X-point collapse' first described by Syrovatskii. Magnetic reconnection within the layers heats a fraction of particles to very high energies. After the saturation of the linear instability, the current layers are disrupted and the system evolves chaotically, diffusing the particle energies in a stochastic second-order Fermi process leading to the formation of power-law energy distributions. The power-law slopes harden with the increasing mean magnetization, but they are significantly softer than those produced in simulations initiated from Harris-type layers. The maximum particle energy is proportional to the mean magnetization, which is attributed partly to the increase of the effective electric field and partly to the increase of the acceleration time scale. We describe in detail the evolving structure of the dynamical current layers, and report on the conservation of magnetic helicity. These results can be applied to highly magnetized astrophysical environments, where ideal plasma instabilities trigger rapid magnetic dissipation with efficient particle acceleration and flares of high-energy radiation.

preprint2016arXiv

Kinetic study of radiation-reaction-limited particle acceleration during the relaxation of unstable force-free equilibria

Many powerful and variable gamma-ray sources, including pulsar wind nebulae, active galactic nuclei and gamma-ray bursts, seem capable of accelerating particles to gamma-ray emitting energies efficiently over very short time scales. These are likely due to rapid dissipation of electromagnetic energy in a highly magnetized, relativistic plasma. In order to understand the generic features of such processes, we have investigated simple models based on relaxation of unstable force-free magnetostatic equilibria. In this work, we make the connection between the corresponding plasma dynamics and the expected radiation signal, using 2D particle-in-cell simulations that self-consistently include synchrotron radiation reaction. We focus on the lowest order unstable force-free equilibrium in a 2D periodic box. We find that rapid variability, with modest apparent radiation efficiency as perceived by a fixed observer, can be produced during the evolution of the instability. The "flares" are accompanied by an increased polarization degree in the high energy band, with rapid variation in the polarization angle. Furthermore, the separation between the acceleration sites and the synchrotron radiation sites for the highest energy particles facilitates acceleration beyond the synchrotron radiation reaction limit. We also discuss the dynamical consequences of radiation reaction, and some astrophysical applications of this model. Our current simulations with numerically tractable parameters are not yet able to reproduce the most dramatic gamma-ray flares, e.g., from Crab Nebula. Higher magnetization studies are promising and will be carried out in the future.

preprint2016arXiv

Systematic Study of Gamma-ray bright Blazars with Optical Polarization and Gamma-ray Variability

Blazars are highly variable active galactic nuclei which emit radiation at all wavelengths from radio to gamma-rays. Polarized radiation from blazars is one key piece of evidence for synchrotron radiation at low energies and it also varies dramatically. The polarization of blazars is of interest for understanding the origin, confinement, and propagation of jets. However, even though numerous measurements have been performed, the mechanisms behind jet creation, composition and variability are still debated. We performed simultaneous gamma-ray and optical photopolarimetry observations of 45 blazars between Jul. 2008 and Dec. 2014 to investigate the mechanisms of variability and search for a basic relation between the several subclasses of blazars. We identify a correlation between the maximum degree of optical linear polarization and the gamma-ray luminosity or the ratio of gamma-ray to optical fluxes. Since the maximum polarization degree depends on the condition of the magnetic field (chaotic or ordered), this result implies a systematic difference in the intrinsic alignment of magnetic fields in pc-scale relativistic jets between different types blazars (FSRQs vs. BL Lacs), and consequently between different types of radio galaxies (FR Is vs. FR IIs).

preprint2015arXiv

Turbulent spectra of the brightest gamma-ray flares of blazars

We investigate the spectral properties of the brightest gamma-ray flares of blazars detected by the Fermi Large Area Telescope. We search for the presence of spectral breaks and measure the spectral curvature on typical time scales of a few days. We identify significant spectral breaks in fewer than half of the analyzed flares, but their parameters do not show any discernible regularities, and in particular there is no indication for gamma-ray absorption at any fixed source-frame photon energy. More interestingly, we find that the studied blazars are characterized by significant spectral variability. Gamma-ray flares of short duration are often characterized by strong spectral curvature, with the spectral peak located above 100 MeV. Since these spectral variations are observed despite excellent photon statistics, they must reflect temporal fluctuations in the energy distributions of the emitting particles. We suggest that highly regular gamma-ray spectra of blazars integrated over long time scales emerge from a superposition of many short-lived irregular components with relatively narrow spectra. This would imply that the emitting particles are accelerated in strongly turbulent environments.

preprint2014arXiv

Constraining the Location of Gamma-Ray Flares in Luminous Blazars

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

preprint2014arXiv

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

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

preprint2014arXiv

Stochastic Modeling of the Fermi/LAT Gamma-ray Blazar Variability

We study the gamma-ray variability of 13 blazars observed with the Fermi Large Area Telescope (LAT). These blazars have the most complete light curves collected during the first 4 years of the Fermi sky survey. We model them with the Ornstein-Uhlenbeck (OU) process or a mixture of the OU processes. The OU process has power spectral density (PSD) proportional to 1/f^alpha with alpha changing at a characteristic time scale, tau_0, from 0 (tau>>tau_0) to 2 (tau<<tau_0). The PSD of the mixed OU process has two characteristic time scales and an additional intermediate region with 0<alpha<2. We show that the OU model provides a good description of the Fermi/LAT light curves of three blazars in our sample. For the first time we constrain a characteristic gamma-ray time scale of variability in two BL Lac sources, 3C 66A and PKS 2155-304 (tau_0=25 day and tau_0=43 day, respectively, in the observer's frame), which are longer than the soft X-ray time scales detected in blazars and Seyfert galaxies. We find that the mixed OU process approximates the light curves of the remaining 10 blazars better than the OU process. We derive limits on their long and short characteristic time scales, and infer that their Fermi/LAT PSDs resemble a power-law function. We constrain the PSD slopes for all but one source in the sample. We find hints for sub-hour Fermi/LAT variability in four flat spectrum radio quasars. We discuss the implications of our results for theoretical models of blazar variability.

preprint2013arXiv

Implications of the Anomalous Outburst in the Blazar PKS 0208-512

The flat spectrum radio quasar (FSRQ) PKS 0208-512 underwent three outbursts at the optical-near-infrared (OIR) wavelengths during 2008-2011. The second OIR outburst did not have a gamma-ray counterpart despite being comparable in brightness and temporal extent to the other two. We model the time variable spectral energy distribution of PKS 0208-512 during those three flaring episodes with leptonic models to investigate the physical mechanism that can produce this anomalous flare. We show that the redder-when-brighter spectral trend in the OIR bands can be explained by the superposition of a fixed thermal component from the accretion disk and a synchrotron component of fixed shape and variable normalization. We estimate the accretion disk luminosity at L_d ~8 X 10^45 erg/s. Using the observed variability timescale in the OIR band t_{var,obs} ~2 d and the X-ray luminosity L_X ~3.5 X 10^45 erg/s, we constrain the location of the emitting region to distance scales that are broadly comparable with the dusty torus. We show that variations in the Compton dominance parameter by a factor of ~4 --- which may result in the anomalous outburst --- can be relatively easily accounted for by moderate variations in the magnetic field strength or the location of the emission region. Since such variations appear to be rare among FSRQs, we propose that most gamma-ray/OIR flares in these objects are produced in jet regions where the magnetic field and external photon fields vary similarly with distance along the jet, e.g., u_B ~u_ext ~r^{-2}.

preprint2013arXiv

On the origin of X-ray spectra in luminous blazars

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

preprint2012arXiv

A Fast Flare and Direct Redshift Constraint in Far-UV Spectra of the Blazar S50716+714

The BL Lacertae object S50716+714 is one of the most studied blazars on the sky due to its active variability and brightness in many bands, including VHE gamma rays. We present here two serendipitous results from recent far-ultraviolet spectroscopic observations by the Cosmic Origins Spectrograph onboard the Hubble Space Telescope. First, during the course of our 7.3 hour HST observations, the blazar increased in flux rapidly by ~40% (-0.45 mag/h) followed by a slower decline (+0.36 mag/h) to previous far-UV flux levels. We model this flare using asymmetric flare templates and constrain the physical size and energetics of the emitting region. Furthermore, the spectral index of the object softens considerably during the course of the flare from alpha(nu)=-1.0 to alpha(nu)=-1.4. Second, we constrain the source redshift directly using the ~30 intervening absorption systems. A system at z=0.2315 is detected in Lya, Lyb, OVI, and NV and defines the lower bound on the source redshift. No absorbers are seen in the remaining spectral coverage (0.2315<z(Lya)<0.47) and we set a statistical upper bound of z<0.322 (95% confidence) on the blazar. This is the first direct redshift limit for this object and is consistent with literature estimates of z=0.31+-0.08 based on the detection of a host galaxy.

preprint2012arXiv

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

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

preprint2012arXiv

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

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

preprint2012arXiv

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

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

preprint2012arXiv

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

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

preprint2012arXiv

Polarization swings from curved trajectories of the emitting regions

We present a model of polarization swings in blazars from axially symmetric blobs propagating on curved trajectories. If the minimum inclination of the velocity vector to the line of sight is smaller than Gamma^{-1}, the polarization angle maximum rotation rate is simultaneous with the polarization degree minimum and a spike in the total flux. By measuring the maximum rotation rate and the moment of the polarization maximum, it is possible to estimate the distance covered by the blob and thus its approximate position. We apply this model to the recent polarization event in blazar 3C 279.

preprint2012arXiv

The brightest gamma-ray flares of blazars

I present a systematic study of gamma-ray flares in blazars. For this purpose, I propose a very simple and practical definition of a flare as a period of time, associated with a given flux peak, during which the flux is above half of the peak flux. I select a sample of 40 brightest gamma-ray flares observed by Fermi/LAT during the first 4 years of its mission. The sample is dominated by 4 blazars: 3C 454.3, PKS 1510-089, PKS 1222+216 and 3C 273. For each flare, I calculate a light curve and variations of the photon index. For the whole sample, I study the distributions of the peak flux, peak luminosity, duration, time asymmetry, average photon index and photon index scatter. I find that: 1) flares produced by 3C 454.3 are longer and have more complex light curves than those produced by other blazars; 2) flares shorter than 1.5 days in the source frame tend to be time-asymmetric with the flux peak preceding the flare midpoint. These differences can be largely attributed to a smaller viewing angle of 3C 454.3 as compared to other blazars. Intrinsically, the gamma-ray emitting regions in blazar jets may be structured and consist of several domains. I find no regularity in the spectral gamma-ray variations of flaring blazars.

preprint2012arXiv

The effect of poloidal velocity shear on the local development of current-driven instabilities

We perform a local (short-wavelength) linear stability analysis of an axisymmetric column of magnetized plasma with a nearly toroidal magnetic field and a smooth poloidal velocity shear by perturbing the equations of relativistic magnetohydrodynamics. We identify two types of unstable modes, which we call 'exponential' and 'overstable', respectively. The exponential modes are present in the static equilibria and their growth rates decrease with increasing velocity shear. The overstable modes are driven by the effects of velocity shear and dominate the exponential modes for sufficiently high shear values. We argue that these local instabilities can provide an important energy dissipation mechanism in astrophysical relativistic jets. Strong co-moving velocity shear arises naturally in the magnetic acceleration mechanism, therefore it may play a crucial role in converting Poynting-flux-dominated jets into matter-dominated jets, regulating the global acceleration and collimation processes, and producing the observed emission of blazars and gamma-ray bursts.

preprint2011arXiv

Dissipation Efficiency of Reconfinement Shocks in Relativistic Jets

We calculate the dissipation efficiency of relativistic reconfinement shocks. Building on previous work (Nalewajko & Sikora 2009), we consider different distributions of the external pressure. The average dissipation efficiency epsilon_diss is a function of the product of two parameters - the jet Lorentz factor Gamma_j and the reconfinement angle Theta_r, which is related to the opening angle Theta_j and the external pressure index eta. The spatial distribution of the dissipation rate strongly depends on eta. We discuss the significance of these results for the properties of relativistic jets in gamma-ray bursts and active galactic nuclei and propose that reconfinement shocks may explain a very high dissipation efficiency of the former and a moderate dissipation efficiency of the latter. Finally, we estimate the dissipation efficiency of the reconfinement shock associated with the quasi-stationary knot HST-1 in the jet of radio galaxy M87 and show that it is roughly consistent with the observational constraints.

preprint2010arXiv

Radiative Properties of Reconnection-Powered Minijets in Blazars

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

preprint2009arXiv

Constraining Emission Models of Luminous Blazar Sources

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