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.

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. This suggests that the jets powering the lobes are dominated by pairs not only in number but also in enthalpy. Such a possibility is confronted against the constraints imposed by the jet formation scenarios, and by the properties of jets on parsec and kiloparsec scales. We find a lower limit on the pair content from the energetic dominance of leptons over protons in the radio lobes, which exceeds the estimate of $n_e/n_p \sim 20$ indicated by models of the blazar zones and of the radio cores. This allows us to establish that the average energy of protons in the jet co-moving frame is $<2m_pc^2$, and to estimate the number fluxes of protons and $e^{+}e^{-}$ pairs. The required proton flux is achievable, if the jet carries away about $1\%$ of the accreting matter. Loading the jets by pairs can be achieved by pair creation using soft gamma-ray photons from the accretion disk corona, provided that about $2-3\%$ of the disk radiation extends beyond $1\;{\rm MeV}$. Finally, we discuss possible dissipative mechanisms which may keep the jets dominated by relativistically hot pairs.

The first image of the black hole (BH) M87* obtained by the Event Horizon Telescope (EHT) has the shape of a crescent extending from the E to SWW position angles, while the observed direction of the large-scale jet is NWW. Images based on numerical simulations of BH accretion flows suggest that on average the projected BH spin axis would 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, (3) the emissivities are isotropic in the local rest frames. Emission from the crescent sector SSE - SWW 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 SEE sector. Most likely, the SEE '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.

Characterisation of the long-term variations of the broad line region in a luminous blazar, where Comptonization 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 analyse 10 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. The optical spectra are dominated by a highly variable non-thermal synchrotron continuum with a prominent Mg II broad emission line. The line flux is obtained by spectral decomposition including significant contribution from the Fe II pseudo-continuum. Three methods are used to characterise variations of the line flux: (1) stacking of the continuum-subtracted spectra, (2) subtracting the running mean light curves calculated for different time scales, (3) Discrete Correlation Function (DCF) for evaluating potential time delays. Despite very large variations of 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 of 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 increase in the line luminosity is delayed by ~600 days. If this time delay results from reverberation of poorly constrained accretion disk 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 have not identified additional examples of short-term 'flares' of the line flux, besides the case of MJD 55512 reported by Le\'on-Tavares et al. (2013).

PKS 1510-089 is one of the most variable blazars in the third Fermi-LAT source catalog. During 2015, this source has shown four flares identified as flare A, B, C, and D in between three quiescent states Q1, Q2, and Q3. The multi-wavelength data from Fermi-LAT, Swift-XRT/UVOT, OVRO, and SMA observatory are used in our work to model these states. Different flux doubling times have been observed in different energy bands which indicate there could be multiple emission zones. The flux doubling time from the gamma-ray and X-ray light curves are found to be 10.6 hr, 2.5 days, and the average flux doubling time in the optical/UV band is 1 day. It is possible that the gamma-ray and optical/UV emission are produced in the same region whereas X-ray emission is coming from a different region along the jet axis. We have also estimated the discrete correlations functions (DCFs) among the light curves of different energy bands to infer about their emission regions. However, our DCF analysis does not show significant correlation in different energy bands though it shows peaks in some cases at small time lags. We perform a two-zone multi-wavelength time-dependent SED modeling with one emission zone located near the outer edge of the broad line region (BLR) and another further away in the dusty/molecular torus (DT/MT) region to study this high state.

Active galactic nuclei (AGN) with relativistic jets, powered by gas accretion onto their central supermassive black hole (SMBH), are unique laboratories for studying the physics of matter and elementary particles in extreme conditions that cannot be realized on Earth. For a long time since the discovery of AGN, photons were the only way to probe the underlying physical processes. The recent discovery of a very high energy neutrino, IceCube-170922A, coincident with a flaring blazar, TXS 0506+056, provides the first evidence that AGN jets are multi-messenger sources; they are capable of accelerating hadrons to very high energies, while producing non-thermal EM radiation and high-energy neutrinos. This new era of multi-messenger astronomy, which will mature in the next decade, offers us the unprecedented opportunity to combine more than one messenger to solve some long-standing puzzles of AGN jet physics: How do jets dissipate their energy to accelerate particles? What is the jet total kinetic power? Where and how do jets produce high-energy emission and neutrinos? What physical mechanisms drive the particle acceleration?

We present the results of three-dimensional kinetic particle-in-cell (PIC) simulations of isotropic periodic relativistically magnetized pair-plasma equilibria known as the ABC fields. We performed several simulations for initial wavenumbers k_ini = 2 or k_ini = 4, different efficiencies of radiative cooling (including radiation reaction from synchrotron and inverse Compton processes), and different mean magnetization values. These equilibria evolve by means of ideal coalescence instability, the saturation of which generates ab initio localized kinetically-thin current layers -- sites of magnetic reconnection and non-thermal particle acceleration -- eventually relaxing to a state of lower magnetic energy at conserved total magnetic helicity. We demonstrate that magnetic relaxation involves in addition localized collapses of magnetic minima and bulk mergers of current layer pairs, which represents a novel scenario of spontaneous magnetic dissipation with application to the rapid gamma-ray flares of blazars and of the Crab Nebula. Particle acceleration under strong radiative losses leads to formation of power-law indices N(gamma) ~ gamma^(-p) up to p ~= -2.3 at mean hot magnetization values of <sigma_hot> ~ 6. Individual energetic particles can be accelerated within one light-crossing time by electric fields that are largely perpendicular to the local magnetic fields. The energetic particles are highly anisotropic due to the kinetic beaming effect, implying complex patterns of rapid variability. A significant fraction of the initial total energy can be radiated away in the overall process of magnetoluminescence.

Cosmic sources of gamma-ray radiation in the GeV range are often characterized by violent variability, in particular this concerns blazars, gamma-ray bursts, and the pulsar wind nebula Crab. Such gamma-ray emission requires a very efficient particle acceleration mechanism. If the environment, in which such emission is produced, is relativistically magnetized (i.e., that magnetic energy density dominates even the rest-mass energy density of matter), then the most natural mechanism of energy dissipation and particle acceleration is relativistic magnetic reconnection. Basic research into this mechanism is performed by means of kinetic numerical simulations of various configurations of collisionless relativistic plasma with the use of the particle-in-cell algorithm. Such technique allows to investigate the details of particle acceleration mechanism, including radiative energy losses, and to calculate the temporal, spatial, spectral and angular distributions of synchrotron and inverse Compton radiation. The results of these simulations indicate that the effective variability time scale of the observed radiation can be much shorter than the light-crossing time scale of the simulated domain.

The BL Lac object S5~0716+714, a highly variable blazar, underwent an impressive outburst in January 2015 (Phase A), followed by minor activity in February (Phase B). The MAGIC observations were triggered by the optical flux observed in Phase A, corresponding to the brightest ever reported state of the source in the R-band. The comprehensive dataset collected is investigated in order to shed light on the mechanism of the broadband emission. Multi-wavelength light curves have been studied together with the broadband Spectral Energy Distributions (SEDs). The data set collected spans from radio, optical photometry and polarimetry, X-ray, high-energy (HE, 0.1 GeV < E < 100 GeV) with \textit{Fermi}-LAT to the very-high-energy (VHE, E>100 GeV) with MAGIC. The flaring state of Phase A was detected in all the energy bands, providing for the first time a multi-wavelength sample of simultaneous data from the radio band to the VHE. In the constructed SED the \textit{Swift}-XRT+\textit{NuSTAR} data constrain the transition between the synchrotron and inverse Compton components very accurately, while the second peak is constrained from 0.1~GeV to 600~GeV by \textit{Fermi}+MAGIC data. The broadband SED cannot be described with a one-zone synchrotron self-Compton model as it severely underestimates the optical flux in order to reproduce the X-ray to $\gamma$-ray data. Instead we use a two-zone model. The EVPA shows an unprecedented fast rotation. An estimation of the redshift of the source by combined HE and VHE data provides a value of $z = 0.31 \pm 0.02_{stats} \pm 0.05_{sys}$, confirming the literature value. The data show the VHE emission originating in the entrance and exit of a superluminal knot in and out a recollimation shock in the inner jet. A shock-shock interaction in the jet seems responsible for the observed flares and EVPA swing. This scenario is also consistent with the SED modelling.

First results are presented from kinetic numerical simulations of relativistic collisionless magnetic reconnection in pair plasma that include radiation reaction from both synchrotron and inverse Compton (IC) processes, motivated by non-thermal high-energy astrophysical sources, including in particular blazars. These simulations are initiated from a configuration known as 'ABC fields' that evolves due to coalescence instability and generates thin current layers in its linear phase. Global radiative efficiencies, instability growth rates, time-dependent radiation spectra, lightcurves, variability statistics and the structure of current layers are investigated for a broad range of initial parameters. We find that the IC radiative signatures are generally similar to the synchrotron signatures. The luminosity ratio of IC to synchrotron spectral components, the Compton dominance, can be modified by more than one order of magnitude with respect to its nominal value. For very short cooling lengths, we find evidence for modification of the temperature profile across the current layers, no systematic compression of plasma density, and very consistent profiles of E.B. We decompose the profiles of E.B with the use of the Vlasov momentum equation, demonstrating a contribution from radiation reaction at the thickness scale consistent with the temperature profile.

Magnetic reconnection in relativistic collisionless plasmas can accelerate particles and power high-energy emission in various astrophysical systems. Whereas most previous studies focused on relativistic reconnection in pair plasmas, less attention has been paid to electron-ion plasma reconnection, expected in black hole accretion flows and relativistic jets. We report a comprehensive particle-in-cell numerical investigation of reconnection in an electron-ion plasma, spanning a wide range of ambient ion magnetizations $\sigma_i$, from the semirelativistic regime (ultrarelativistic electrons but nonrelativistic ions, 0.001<<$\sigma_i$<<1) to the fully relativistic regime (both species are ultrarelativistic, $\sigma_i$>>1). We investigate how the reconnection rate, electron and ion plasma flows, electric and magnetic field structures, electron/ion energy partitioning, and nonthermal particle acceleration depend on $\sigma_i$. Our key findings are: (1) the reconnection rate is about 0.1 of the Alfvenic rate across all regimes; (2) electrons can form concentrated moderately relativistic outflows even in the semirelativistic, small-$\sigma_i$ regime; (3) while the released magnetic energy is partitioned equally between electrons and ions in the ultrarelativistic limit, the electron energy fraction declines gradually with decreased $\sigma_i$ and asymptotes to about 0.25 in the semirelativistic regime; (4) reconnection leads to efficient nonthermal electron acceleration with a $\sigma_i$-dependent power-law index, $p(\sigma_i) \simeq $const$+0.7 {\sigma_i}^{-1/2}$. These findings are important for understanding black hole systems and lend support to semirelativistic reconnection models for powering nonthermal emission in blazar jets, offering a natural explanation for the spectral indices observed in these systems.

Recent detection of suborbital gamma-ray variability of Flat Spectrum Radio Quasar (FSRQ) 3C 279 by Fermi Large Area Telescope (LAT) is in severe conflict with established models of blazar emission. This paper presents the results of suborbital analysis of the Fermi/LAT data for the brightest gamma-ray flare of another FSRQ blazar 3C 454.3 in November 2010 (MJD 55516-22). Gamma-ray light curves are calculated for characteristic time bin lengths as short as 3 min. The measured variations of the 0.1-10 GeV photon flux are tested against the hypothesis of steady intraorbit flux. In addition, the structure function is calculated for absolute photon flux differences and for their significances. Significant gamma-ray flux variations are measured only over time scales longer than ~5h, which is consistent with the standard blazar models.

This paper presents a simple model of polarisation rotation in optically thin relativistic jets of blazars. The model is based on the development of helical (kink) mode of current-driven instability. A possible explanation is suggested for the observational connection between polarisation rotations and optical/gamma-ray flares in blazars, if the current-driven modes are triggered by secular increases of the total jet power. The importance of intrinsic depolarisation in limiting the amplitude of coherent polarisation rotations is demonstrated. The polarisation rotation amplitude is thus very sensitive to the viewing angle, which appears to be inconsistent with the observational estimates of viewing angles in blazars showing polarisation rotations. Overall, there are serious obstacles to explaining large-amplitude polarisation rotations in blazars in terms of current-driven kink modes.

The two-component broad-band spectral energy distributions of blazars were suggested to form a sequence in which (1) the peak frequency of the low-energy (synchrotron) component nu_syn is anticorrelated with the synchrotron luminosity L_syn, and (2) the luminosity ratio of the high-energy (inverse Compton) to synchrotron components q = L_IC / L_syn (Compton dominance) increases with L_syn from the BL Lac objects (BL Lacs) to the flat spectrum radio quasars (FSRQs). The Compton dominance parameter is an important probe of plasma magnetisation in the blazar zones within relativistic jets. We investigate a sample of blazars detected by WISE in the mid-infrared (MIR) band and by Fermi/LAT in the GeV gamma-ray band, with the focus on the distribution of luminosities and photon indices. Our findings are the following: (1) the MIR photon index Gamma_W12 is a useful probe of the blazar sequence, with the exception of low-luminosity BL Lacs that are most likely contaminated by their host galaxies (L_W1 ~ 10^44 erg/s and Gamma_W12 < 1); (2) Gamma_W12 is correlated with the gamma-ray photon index Gamma_1-100GeV, with the MIR luminosity L_W1, and with the Fermi/WISE Compton dominance q_FW = L_1GeV/L_W1; (3) a clean separation between FSRQs and BL Lacs can be seen in the parameter space of Gamma_W12 and q_FW; (4) the observed distribution of MIR luminosity L_W1 vs. Compton dominance q_FW for the entire sample of blazars can be modeled as a sequence of lepto-magnetic jet powers in the range (log_10 P_eB) in [42:45] with the preference for sub-equipartition magnetic fields P_B/P_e in [0.05:1], assuming fixed bulk Lorentz factor Gamma_j = 15, fixed jet opening angle Gamma_j Theta_j = 0.3, fixed radiative efficiency of jet electrons epsilon_em = 50%, and that external radiation luminosity scales like L_ext propto (P_eB)^(1.6) (parameter degeneracies are discussed).

The flat spectrum radio quasar 3C 279 is a known $\gamma$-ray variable source that has recently exhibited minute-scale variability at energies $>100$ MeV. One-zone leptonic models for blazar emission are severely constrained by the short timescale variability that implies a very compact emission region at a distance of hundreds of Schwarzschild radii from the central black hole. Here, we investigate a hadronic scenario where GeV $\gamma$-rays are produced via proton synchrotron radiation. We also take into account the effects of the hadronically initiated electromagnetic cascades (EMC). For a $\gamma$-ray emitting region in rough equipartition between particles and kG magnetic fields, located within the broad-line region (BLR), the development of EMC redistributes the $\gamma$-ray luminosity to softer energy bands and eventually leads to broad-band spectra that differ from the observed ones. Suppression of EMC and energy equipartition are still possible, if the $\gamma$-ray emitting region is located beyond the BLR, is fast moving with Doppler factor ($>70$), and contains strong magnetic fields ($>100$ G). Yet, these conditions cannot be easily met in parsec-scale jets, thus disfavouring a proton synchrotron origin of the Fermi-LAT flare.

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).

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.

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.

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.

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.

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.

Using two-dimensional particle-in-cell simulations, we characterize the energy spectra of particles accelerated by relativistic magnetic reconnection (without guide field) in collisionless electron-positron plasmas, for a wide range of upstream magnetizations $\sigma$ and system sizes $L$. The particle spectra are well-represented by a power law $\gamma^{-\alpha}$, with a combination of exponential and super-exponential high-energy cutoffs, proportional to $\sigma$ and $L$, respectively. For large $L$ and $\sigma$, the power-law index $\alpha$ approaches about 1.2.

Jets associated with Active Galactic Nuclei (AGN) have been observed for almost a century, initially at optical and radio wavelengths. They are now widely accepted as "exhausts" produced electromagnetically by the central, spinning, massive black hole and its orbiting, accreting gas. Observations at X-ray and, especially, gamma-ray energies have transformed our understanding of how these jets evolve dynamically, accelerate electrons (and positrons) and radiate throughout the entire electromagnetic spectrum. Some new approaches to modeling the powerful and rapidly variable TeV emission observed from many blazars are sketched. Observations at the highest TeV energies, to which the High Altitude Water Cherenkov Gamma-Ray Observatory (HAWC) will contribute, promise crucial discrimination between rival models of AGN jets.

We investigate the distribution of particle acceleration sites, independently of the actual acceleration mechanism, during plasmoid-dominated, relativistic collisionless magnetic reconnection by analyzing the results of a particle-in-cell numerical simulation. The simulation is initiated with Harris-type current layers in pair plasma with no guide magnetic field, negligible radiative losses, no initial perturbation, and using periodic boundary conditions. We find that the plasmoids develop a robust internal structure, with colder dense cores and hotter outer shells, that is recovered after each plasmoid merger on a dynamical time scale. We use spacetime diagrams of the reconnection layers to probe the evolution of plasmoids, and in this context we investigate the individual particle histories for a representative sample of energetic electrons. We distinguish three classes of particle acceleration sites associated with (1) magnetic X-points, (2) regions between merging plasmoids, and (3) the trailing edges of accelerating plasmoids. We evaluate the contribution of each class of acceleration sites to the final energy distribution of energetic electrons -- magnetic X-points dominate at moderate energies, and the regions between merging plasmoids dominate at higher energies. We also identify the dominant acceleration scenarios, in order of decreasing importance -- (1) single acceleration between merging plasmoids, (2) single acceleration at a magnetic X-point, and (3) acceleration at a magnetic X-point followed by acceleration in a plasmoid. Particle acceleration is absent only in the vicinity of stationary plasmoids. The effect of magnetic mirrors due to plasmoid contraction does not appear to be significant in relativistic reconnection.

We report on simultaneous broadband observations of the TeV-emitting blazar Markarian 501 between 1 April and 10 August 2013, including the first detailed characterization of the synchrotron peak with Swift and NuSTAR. During the campaign, the nearby BL Lac object was observed in both a quiescent and an elevated state. The broadband campaign includes observations with NuSTAR, MAGIC, VERITAS, the Fermi Large Area Telescope (LAT), Swift X-ray Telescope and UV Optical Telescope, various ground-based optical instruments, including the GASP-WEBT program, as well as radio observations by OVRO, Mets\"ahovi and the F-Gamma consortium. Some of the MAGIC observations were affected by a sand layer from the Saharan desert, and had to be corrected using event-by-event corrections derived with a LIDAR (LIght Detection And Ranging) facility. This is the first time that LIDAR information is used to produce a physics result with Cherenkov Telescope data taken during adverse atmospheric conditions, and hence sets a precedent for the current and future ground-based gamma-ray instruments. The NuSTAR instrument provides unprecedented sensitivity in hard X-rays, showing the source to display a spectral energy distribution between 3 and 79 keV consistent with a log-parabolic spectrum and hard X-ray variability on hour timescales. None (of the four extended NuSTAR observations) shows evidence of the onset of inverse-Compton emission at hard X-ray energies. We apply a single-zone equilibrium synchrotron self-Compton model to five simultaneous broadband spectral energy distributions. We find that the synchrotron self-Compton model can reproduce the observed broadband states through a decrease in the magnetic field strength coinciding with an increase in the luminosity and hardness of the relativistic leptons responsible for the high-energy emission.

We present results from a 244\,ks \textit{NuSTAR} observation of 3C\,273 obtained during a cross-calibration campaign with the \textit{Chandra}, \textit{INTEGRAL}, \textit{Suzaku}, \textit{Swift}, and \textit{XMM-Newton} observatories. We show that the spectrum, when fit with a power-law model using data from all observatories except \textit{INTEGRAL} over the 1--78\,keV band, leaves significant residuals in the \textit{NuSTAR} data between 30--78\,keV. The \nustar\ 3--78\,keV spectrum is well-described by an exponentially cutoff power-law ($\Gamma = 1.646 \pm 0.006$, E$_\mathrm{cutoff} = 202_{-34}^{+51}$\,keV) with a weak reflection component from cold, dense material. There is also evidence for a weak ($EW = 23 \pm 11$ eV) neutral iron line. We interpret these features as arising from coronal emission plus reflection off an accretion disk or distant material. Beyond 80\,keV \textit{INTEGRAL} data show clear excess flux relative to an extrapolation of the AGN model fit to \nustar. This high-energy power-law is consistent with the presence of a beamed jet, which begins to dominate over emission from the inner accretion flow at 30-40 keV. Modeling the jet locally (in the \textit{NuSTAR} + \textit{INTEGRAL} band) as a power-law, we find the coronal component is fit by $\Gamma_\mathrm{AGN} = 1.638 \pm 0.045$, $E_\mathrm{cutoff} = 47 \pm 15$\,keV, and jet photon index by $\Gamma_\mathrm{jet} = 1.05 \pm 0.4$. We also consider \textit{Fermi}/LAT observations of 3C\,273 and here the broad-band spectrum of the jet can be described by a log-parabolic model, peaking at $\sim 2$\,MeV. Finally, we investigate the spectral variability in the \textit{NuSTAR} band and find an inverse correlation between flux and $\Gamma$.

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 $\gamma$-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 $\gamma$-ray spectrum during one of the flares was very hard, with an index of $\Gamma_\gamma = 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_\gamma/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 $\Delta\Gamma_{\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 $\gamma$-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 $\gamma$-ray and optical spectra cannot satisfactorily explain the production of X-ray emission.

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.

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.

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.

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.

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.

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}.

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.

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.

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.

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.

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.

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.

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.

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 {\gamma} -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 {\gamma} -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.

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.

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.

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.

Tracing the formation and evolution of all supermassive black holes, including the obscured ones, understanding how black holes influence their surroundings and how matter behaves under extreme conditions, are recognized as key science objectives to be addressed by the next generation of instruments. These are the main goals of the COSPIX proposal, made to ESA in December 2010 in the context of its call for selection of the M3 mission. In addition, COSPIX, will also provide key measurements on the non thermal Universe, particularly in relation to the question of the acceleration of particles, as well as on many other fundamental questions as for example the energetic particle content of clusters of galaxies. COSPIX is proposed as an observatory operating from 0.3 to more than 100 keV. The payload features a single long focal length focusing telescope offering an effective area close to ten times larger than any scheduled focusing mission at 30 keV, an angular resolution better than 20 arcseconds in hard X-rays, and polarimetric capabilities within the same focal plane instrumentation. In this paper, we describe the science objectives of the mission, its baseline design, and its performances, as proposed to ESA.

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.

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.

We study the polarization properties of relativistic reconfinement shocks with chaotic magnetic fields. Using our hydrodynamical model of their structure, we calculate synthetic polarization maps, longitudinal polarization profiles and discuss the spatially averaged polarization degree as a function of jet half-opening angle Theta_j, jet Lorentz factor Gamma_j and observer inclination angle to the jet axis theta_{obs}. We find, that for theta_{obs} <= Theta_j the wave electric vectors are parallel in the vicinity of the structure ends and perpendicular in between, while for theta_{obs} > Theta_j the polarization can only be perpendicular. The spatially averaged polarization degree does not exceed 30%. Parallel average polarization, with polarization degrees lower than 10%, have been found for theta_{obs} < Theta_j under the condition Gamma_j * Theta_j > 1. As earlier works predicted the parallel polarization from relativistic conical shocks, we explain our results by discussing conical shocks with divergent upstream flow.

We present a semi-analytical hydrodynamical model for the structure of reconfinement shocks formed in astrophysical relativistic jets interacting with external medium. We take into account exact conservation laws, both across the shock front and in the zone of the shocked matter, and exact angular relations. Our results confirm a good accuracy of the approximate formulae derived by Komissarov & Falle (1997). However, including the transverse pressure gradient in the shocked jet, we predict an absolute size of the shock to be about about twice larger. We calculate the efficiency of the kinetic energy dissipation in the shock and show a strong dependence on both the bulk Lorentz factor and opening angle of the jet.

We have obtained single-phase near-infrared magnitudes in the J and K bands for a sample of 78 RR Lyrae stars in the Sculptor dSph galaxy. Applying different theoretical and empirical calibrations of the period-luminosity-metallicity relation for RR Lyrae stars in the infrared, we find consistent results and obtain a true, reddening-corrected distance modulus of 19.67 $\pm$ 0.02 (statistical) $\pm$ 0.12 (systematic) mag for Sculptor from our data. This distance value is consistent with the value of 19.68 $\pm$ 0.08 mag which we obtain from earlier V-band data of RR Lyrae stars in Sculptor, and the V magnitude-metallicity calibration of Sandage (1993). It is also in a very good agreement with the results obtain by Rizzi (2002) based on tip of the red giant branch (TRGB, 19.64 $\pm$ 0.08 mag) and horizontal branch (HB, 19.66 $\pm$ 0.15 mag).

We have measured near-infrared magnitudes in the J and K bands for 56 Cepheid variables in the Local Group galaxy NGC 6822 with well-determined periods and optical light curves in the V and I bands. Using the template light curve approach of Soszynski, Gieren and Pietrzynski, accurate mean magnitudes were obtained from these data which allowed us to determine with unprecedented accuracy the distance to NGC 6822 from a multi-wavelength period-luminosity solution in the VIJK bands. From our data, we obtain a distance to NGC 6822 of (m-M)_{0} = 23.312 +- 0.021 (random error) mag, with an additional systematic uncertainty of about 3 %. This distance value is tied to an assumed LMC distance modulus of 18.50. From our multiwavelength approach, we find for the total (average) reddening to the NGC 6822 Cepheids E(B-V) = 0.356 +- 0.013 mag, which is in excellent agreement with a previous determination of McGonegal et al. from near-infrared photometry and implies significant internal reddening of the Cepheids in NGC 6822. Our present, definitive distance determination of NGC 6822 from Cepheids agrees within 2 % with the previous distance we had derived from optical photometry alone, but has significantly reduced error bars. Our Cepheid distance to NGC 6822 is in excellent agreement with the recent independent determination of Cioni and Habing from the I-band magnitude of the tip of the red giant branch. It also agrees well, within the errors, with the early determination of McGonegal et al. (1983) from random-phase H-band photometry of nine Cepheids.