Source author record

Markus Boettcher

Markus Boettcher 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

42works
7topics
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

42 published item(s)

preprint2022arXiv

Multiwavelength and Multimessenger Observations of Blazars and Theoretical Modeling: Blazars as Astrophysical Neutrino Sources

This contribution reviews recent advances in the possible identification of blazars as potential sources of at least some of the very-high-energy neutrinos detected by the IceCube neutrino detector at the South Pole. The basic physical requirements for neutrino production and physics constraints that may be drawn from neutrino - blazar associations are reviewed. Several individual cases of possible associations will be discussed in more detail. It is emphasized that due to $γγ$ opacity constraints in efficiently neutrino-producing blazars, an association between X-ray -- soft $γ$-ray activity and very-high-energy neutrino production is more naturally expected than a connection between neutrino and high-energy/very-high-energy $γ$-ray activity.

preprint2022arXiv

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

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

preprint2020arXiv

gSeaGen: the KM3NeT GENIE-based code for neutrino telescopes

The gSeaGen code is a GENIE-based application developed to efficiently generate high statistics samples of events, induced by neutrino interactions, detectable in a neutrino telescope. The gSeaGen code is able to generate events induced by all neutrino flavours, considering topological differences between track-type and shower-like events. Neutrino interactions are simulated taking into account the density and the composition of the media surrounding the detector. The main features of gSeaGen are presented together with some examples of its application within the KM3NeT project.

preprint2019arXiv

Multiwavelength Variability Signatures of Relativistic Shocks in Blazar Jets

Mildly-relativistic shocks that are embedded in colliding magnetohydrodynamic flows are prime sites for relativistic particle acceleration and production of strongly variable, polarized multi-wavelength emission from relativistic jet sources such as blazars and gamma-ray bursts. The principal energization mechanisms at these shocks are diffusive shock acceleration and shock drift acceleration. In recent work, we had self-consistently coupled shock acceleration and radiation transfer simulations in blazar jets in a basic one-zone scenario. These one-zone models revealed that the observed spectral energy distributions (SEDs) of blazars strongly constrain the nature of the shock layer hydromagnetic turbulence. In this paper, we expand our previous work by including full time dependence and treating two zones, one being the site of acceleration, and the second being a larger emission zone. This construction is applied to multiwavelength flares of the flat spectrum radio quasar 3C 279, fitting snap-shot SEDs and generating light curves that are consistent with observed variability timescales. We also present a generic study for the typical flaring behavior of the BL Lac object Mrk 501. The model predicts correlated variability across all wavebands, but cross-band time lags depending on the type of blazar (FSRQ vs. BL Lac), as well as distinctive spectral hysteresis patterns in all wavelength bands, from mm radio waves to gamma-rays. These evolutionary signatures serve to provide diagnostics on the competition between acceleration and radiative cooling.

preprint2019arXiv

Time-Dependent, Multi-Wavelength Shock Acceleration Models for Active Flares of 3C 279

Jets in blazars are an excellent forum for studying acceleration at relativistic shocks using the highly-variable emission seen across the electromagnetic spectrum. Our recent work on combining multi-wavelength leptonic emission models with simulated thermal+non-thermal distributions from shock acceleration theory has resulted in new insights into plasma conditions in blazars. This has demonstrated the ability to infer the cyclotron frequency, the plasma density and thus also the Alfven speed, thereby determining the rapidity of particle energization. An important inference was that MHD turbulence levels decline with remoteness from jet shocks. This paper outlines new results from our recent extension of this program to a two-zone, time-evolving construction, modeling together both extended, enhanced emission states from larger radiative regions, and prompt flare events from compact acceleration zones. These are applied to flares in the FSRQ blazar 3C 279 monitored by Fermi-LAT in gamma-rays in late 2013. With impulsive injection episodes from the shock zone, as the acceleration first proceeds and then abates, the radiative simulations obtain a pronounced spectral hardening in the optical and gamma-ray bands as the flare grows, followed by a softening during the decay phase. For 3C 279, while model radio and X-ray synchrotron flares are temporally correlated, there is a lag in both bands relative to GeV gamma rays and optical emission on timescales of a number of hours. This delay is governed by the short cooling time associated with the bright external Compton signal.

preprint2016arXiv

3C 279 in Outburst in June 2015: A Broadband-SED Study Based on the INTEGRAL Detection

Blazars radiate from radio through gamma-ray frequencies thereby being ideal targets for multifrequency studies. Such studies allow constraining the properties of the emitting jet. 3C 279 is among the most notable blazars and therefore subject to extensive multifrequency campaigns. We report the results of a campaign ranging from near-IR to gamma-ray energies of an outburst of 3C 279 in June 2015. The overall spectral energy distribution from near-IR to gamma rays can be well represented by either a leptonic or a lepto-hadronic radiation transfer model. Even though the data are equally well represented by the two models, their inferred parameters challenge the physical conditions in the jet. In fact, the leptonic model requires parameters with a magnetic field far below equipartition with the relativistic particle energy density. On the contrary, equipartition may be achieved with the lepto-hadronic model, which however implies an extreme total jet power close to Eddington luminosity.

preprint2016arXiv

Gamma-Gamma Absorption in the Broad Line Region Radiation Fields of Gamma-Ray Blazars

The expected level of gamma-gamma absorption in the Broad Line Region (BLR) radiation field of gamma-ray loud Flat Spectrum Radio Quasars (FSRQs)is evaluated as a function of the location of the gamma-ray emission region. This is done self-consistently with parameters inferred from the shape of the spectral energy distribution (SED) in a single-zone leptonic EC-BLR model scenario. We take into account all geometrical effects both in the calculation of the gamma-gamma opacity and the normalization of the BLR radiation energy density. As specific examples, we study the FSRQs 3C279 and PKS 1510-089, keeping the BLR radiation energy density at the location of the emission region fixed at the values inferred from the SED. We confirm previous findings that the optical depth due to $γγ$ absorption in the BLR radiation field exceeds unity for both 3C279 and PKS 1510-089 for locations of the gamma-ray emission region inside the inner boundary of the BLR. It decreases monotonically, with distance from the central engine and drops below unity for locations within the BLR. For locations outside the BLR, the BLR radiation energy density required for the production of GeV gamma-rays rapidly increases beyond observational constraints, thus making the EC-BLR mechanism implausible. Therefore, in order to avoid significant gamma-gamma absorption by the BLR radiation field, the gamma-ray emission region must therefore be located near the outer boundary of the BLR.

preprint2016arXiv

Gravitational Light Bending Prevents $γγ$ Absorption in Gravitational Lenses

The magnification effect due to gravitational lensing enhances the chances of detecting moderate-redshift ($z \sim 1$) sources in very-high-energy (VHE; $E > 100$ GeV) $γ$-rays by ground-based Atmospheric Cherenkov Telescope facilities. It has been shown in previous work that this prospect is not hampered by potential $γ-γ$ absorption effects by the intervening (lensing) galaxy, nor by any individual star within the intervening galaxy. In this paper, we expand this study to simulate the light-bending effect of a realistic ensemble of stars. We first demonstrate that, for realistic parameters of the galaxy's star field, it is extremely unlikely (probability $\lesssim 10^{-6}$) that the direct line of sight between the $γ$-ray source and the observer passes by any star in the field close enough to be subject to significant $γγ$ absorption. Our simulations then focus on the rare cases where $γγ$ absorption by (at least) one individual star might be non-negligible. We show that gravitational light bending will have the effect of avoiding the $γ-γ$ absorption spheres around massive stars in the intervening galaxy. This confirms previous results by Barnacka et al. and re-inforces arguments in favour of VHE $γ$-ray observations of lensed moderate-redshift blazars to extend the redshift range of objects detected in VHE $γ$-rays, and to probe the location of the $γ$-ray emission region in those blazars.

preprint2016arXiv

Leptonic and Lepto-Hadronic Modeling of the Nov. 2010 flare from 3C 454.3

In this study, we use a one-zone leptonic and a lepto-hadronic model to investigate the multi-wavelength emission and the prominent flare of the flat spectrum radio quasar 3C 454.3 in Nov 2010. We perform a parameter study with both models to obtain broadband fits to the spectral energy distribution of 3C 454.3. Starting with the baseline parameters obtained from the fits, we then investigate different flaring scenarios for both models to explain an extreme outburst and spectral hardening of 3C 454.3 that occurred in Nov 2010. We find that the one zone lepto-hadronic model can successfully explain both the broadband multi-wavelength spectral energy distribution and light curves in the optical R, Swift XRT and Fermi gamma-ray bandpasses for 3C 454.3 during quiescence and the peak of the Nov. 2010 flare. We also find that the one-zone leptonic model produces poor fits to the broadband spectra in the X-ray and high energy gamma-ray band passes for the Nov. 2010 flare.

preprint2016arXiv

Spectral and Polarization Signatures of Relativistic Shocks in Blazars

Relativistic shocks are one of the most plausible sites of the emission of strongly variable, polarized multi-wavelength emission from relativistic jet sources such as blazars, via diffusive shock acceleration (DSA) of relativistic particles. This paper summarizes recent results on a self-consistent coupling of diffusive shock acceleration and radiation transfer in blazar jets. We demonstrate that the observed spectral energy distributions (SEDs) of blazars strongly constrain the nature of hydromagnetic turbulence responsible for pitch-angle scattering by requiring a strongly energy-dependent pitch-angle mean free path. The prominent soft X-ray excess ("Big Blue Bump") in the SED of the BL Lac object AO 0235+164 can be modelled as the signature of bulk Compton scattering of external radiation fields by the thermal electron population, which places additional constraints on the level of hydromagnetic turbulence. It has further been demonstrated that internal shocks propagating in a jet pervaded by a helical magnetic field naturally produce polarization-angle swings by 180$^o$, in tandem with multi-wavelength flaring activity, without requiring any helical motion paths or other asymmetric jet structures. The specific application of this model to 3C279 presents the first consistent, simultaneous modeling of snap-shot SEDs, multi-wavelength light curves and time-dependent polarization signatures of a blazar during a polarization-angle (PA) rotation. This model has recently been generalized to a lepto-hadronic model, in which the high-energy emission is dominated by proton synchrotron radiation. It is shown that in this case, the high-energy (X-ray and $γ$-ray) polarization signatures are expected to be significantly more stable (not showing PA rotations) than the low-energy (electron-synchrotron) signatures.

preprint2016arXiv

The Aid of Optical Studies in Understanding Millisecond Pulsar Binaries

A large number of new "black widow" and "redback" energetic millisecond pulsars with irradiated stellar companions have been discovered through radio and optical searches of unidentified \textit{Fermi} sources. Synchrotron emission, from particles accelerated up to several TeV in the intrabinary shock, exhibits modulation at the binary orbital period. Our simulated double-peaked X-ray light curves modulated at the orbital period, produced by relativistic Doppler-boosting along the intrabinary shock, are found to qualitatively match those observed in many sources. In this model, redbacks and transitional pulsar systems where the double-peaked X-ray light curve is observed at inferior conjunction have intrinsically different shock geometry than other millisecond pulsar binaries where the light curve is centered at superior conjunction. We discuss, and advocate, how current and future optical observations may aid in constraining the emission geometry, intrabinary shock and the unknown physics of pulsar winds.

preprint2016arXiv

Theoretical Study Of The Effects Of Magnetic Field Geometry On The High-Energy Emission Of Blazars

The knowledge of the structure of the magnetic field inside a blazar jet, as deduced from polarization observations at radio to optical wavelengths, is closely related to the formation and propagation of relativistic jets that result from accretion onto supermassive black holes. However, a largely unexplored aspect of the theoretical understanding of radiation transfer physics in blazar jets has been the magnetic field geometry as revealed by the polarized emission and the connection between the variability in polarization and flux across the spectrum. Here, we explore the effects of various magnetic geometries that can exist inside a blazar jet: parallel, oblique, toroidal, and tangled. We investigate the effects of changing the orientation of the magnetic field, according to the above-mentioned geometries, on the resulting high-energy spectral energy distributions (SEDs) and spectral variability patterns (SVPs) of a typical blazar. We use the MUlti-ZOne Radiation Feedback (MUZORF) model of Joshi et al. (2014) to carry out this study and to relate the geometry of the field to the observed SEDs at X-ray and gamma-ray energies. One of the goals of the study is to understand the relationship between synchrotron and inverse Compton peaks in blazar SEDs and the reason for the appearance of gamma-ray "orphan flares" observed in some blazars. This can be associated with the directionality of the magnetic field, which creates a difference in the radiation field as seen by an observer versus that seen by the electrons in the emission region.

preprint2015arXiv

Hadronic models of blazars require a change of the accretion paradigm

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

preprint2015arXiv

Polarization Swings Reveal Magnetic Energy Dissipation in Blazars

The polarization signatures of the blazar emissions are known to be highly variable. In addition to small fluctuations of the polarization angle around a mean value, sometimes large (> 180^o) polarization angle swings are observed. We suggest that such p henomena can be interpreted as arising from light-travel-time effects within an underlying axisymmetric emission region. We present the first simultaneous fitting of the multi-wavelength spectrum, variability and time-dependent polarization features of a correlated optical and gamma-ray flaring event of the prominent blazar 3C279, which was accompanied by a drastic change of its polarization signatures. This unprecedented combination of spectral, variability, and polarization information in a coherent physical model allows us to place stringent constraints on the particle acceleration and magnetic-field topology in the relativistic jet of a blazar, strongly favoring a scenario in which magnetic energy dissipation is the primary driver of the flare event.

preprint2015arXiv

Studies of Relativistic Jets in Active Galactic Nuclei with SKA

Relativistic jets in active galactic nuclei (AGN) are among the most powerful astrophysical objects discovered to date. Indeed, jetted AGN studies have been considered a prominent science case for SKA, and were included in several different chapters of the previous SKA Science Book (Carilli & Rawlings 2004). Most of the fundamental questions about the physics of relativistic jets still remain unanswered, and await high-sensitivity radio instruments such as SKA to solve them. These questions will be addressed specially through analysis of the massive data sets arising from the deep, all-sky surveys (both total and polarimetric flux) from SKA1. Wide-field very-long-baseline-interferometric survey observations involving SKA1 will serve as a unique tool for distinguishing between extragalactic relativistic jets and star forming galaxies via brightness temperature measurements. Subsequent SKA1 studies of relativistic jets at different resolutions will allow for unprecedented cosmological studies of AGN jets up to the epoch of re-ionization, enabling detailed characterization of the jet composition, magnetic field, particle populations, and plasma properties on all scales. SKA will enable us to study the dependence of jet power and star formation on other properties of the AGN system. SKA1 will enable such studies for large samples of jets, while VLBI observations involving SKA1 will provide the sensitivity for pc-scale imaging, and SKA2 (with its extraordinary sensitivity and dynamic range) will allow us for the first time to resolve and model the weakest radio structures in the most powerful radio-loud AGN.

preprint2015arXiv

Time Dependent Hadronic Modeling of Flat Spectrum Radio Quasars

We introduce a new time-dependent lepto-hadronic model for blazar emission that takes into account the radiation emitted by secondary particles, such as pions and muons, from photo hadronic interactions. Starting from a baseline parameter set guided by a fit to the spectral energy distribution of the blazar 3C 279, we perform a parameter study to investigate the effects of perturbations of the input parameters to mimic different flaring events to study the resulting lightcurves in the optical, X-ray, high energy (HE: E > 100 MeV) and very-high-energy (VHE: E > 100 GeV) gamma-rays as well as the neutrino emission associated with charged-pion and muon decay. We find that flaring events from an increase in the efficiency of Fermi II acceleration will produce a positive correlation between all bandpasses and a marked plateau in the HE gamma-ray lightcurve. We also predict a distinctive dip in the HE lightcurve for perturbations caused by a change in the proton injection spectral index. These plateaus / dips could be a tell tale signature of hadronic models for perturbations that lead to more efficient acceleration of high energy protons in parameter regimes where pion and muon synchrotron emission is non-negligible.

preprint2015arXiv

X-Ray Polarimetry with the Polarization Spectroscopic Telescope Array (PolSTAR)

This paper describes the Polarization Spectroscopic Telescope Array (PolSTAR), a mission proposed to NASA's 2014 Small Explorer (SMEX) announcement of opportunity. PolSTAR measures the linear polarization of 3-50 keV (requirement; goal: 2.5-70 keV) X-rays probing the behavior of matter, radiation and the very fabric of spacetime under the extreme conditions close to the event horizons of black holes, as well as in and around magnetars and neutron stars. The PolSTAR design is based on the technology developed for the Nuclear Spectroscopic Telescope Array (NuSTAR) mission launched in June 2012. In particular, it uses the same X-ray optics, extendable telescope boom, optical bench, and CdZnTe detectors as NuSTAR. The mission has the sensitivity to measure ~1% linear polarization fractions for X-ray sources with fluxes down to ~5 mCrab. This paper describes the PolSTAR design as well as the science drivers and the potential science return.

preprint2014arXiv

Magnetic Field Amplification and Flat Spectrum Radio Quasars

We perform time-dependent, spatially-resolved simulations of blazar emission to evaluate several flaring scenarios related to magnetic-field amplification and enhanced particle acceleration. The code explicitly accounts for light-travel-time effects and is applied to flares observed in the flat spectrum radio quasar (FSRQ) PKS 0208-512, which show optical/γ-ray correlation at some times, but orphan optical flares at other times. Changes in both the magnetic field and the particle acceleration efficiency are explored as causes of flares. Generally, external Compton emission appears to describe the available data better than a synchrotron self-Compton scenario, and in particular orphan optical flares are difficult to produce in the SSC framework. X-ray soft-excesses, γ-ray spectral hardening, and the detections at very high energies of certain FSRQs during flares find natural explanations in the EC scenario with particle acceleration change. Likewise, optical flares with/without γ-ray counterparts can be explained by different allocations of energy between the magnetization and particle acceleration, which may be related to the orientation of the magnetic field relative to the jet flow. We also calculate the degree of linear polarization and polarization angle as a function of time for a jet with helical magnetic field. Tightening of the magnetic helix immediately downstream of the jet perturbations, where flares occur, can be sufficient to explain the increases in the degree of polarization and a rotation by >= 180 degree of the observed polarization angle, if light-travel-time effects are properly considered.

preprint2014arXiv

Particle diffusion and localized acceleration in inhomogeneous AGN jets - Part I: Steady-state spectra

We study the acceleration, transport, and emission of particles in relativistic jets. Localized stochastic particle acceleration, spatial diffusion, and synchrotron as well as synchrotron self-Compton emission are considered in a leptonic model. To account for inhomogeneity, we use a 2D axi-symmetric cylindrical geometry for both relativistic electrons and magnetic field. In this first phase of our work, we focus on steady-state spectra that develop from a time-dependent model. We demonstrate that small isolated acceleration region in a much larger emission volume are sufficient to accelerate particles to high energy. Diffusive escape from these small regions provides a natural explanation for the spectral form of the jet emission. The location of the acceleration regions within the jet is found to affect the cooling break of the spectrum in this diffusive model. Diffusion-caused energy-dependent inhomogeneity in the jets predicts that the SSC spectrum is harder than the synchrotron spectrum. There can also be a spectral hardening towards the high-energy section of the synchrotron spectrum, if particle escape is relatively slow. These two spectral hardening effects indicate that the jet inhomogeneity might be a natural explanation for the unexpected hard γ-ray spectra observed in some blazars.

preprint2014arXiv

Seed Photon Fields of Blazars in the Internal Shock Scenario

We extend our approach of modeling spectral energy distribution (SED) and lightcurves of blazars to include external Compton (EC) emission due to inverse Compton scattering of an external anisotropic target radiation field. We describe the time-dependent impact of such seed photon fields on the evolution of multifrequency emission and spectral variability of blazars using a multi-zone time-dependent leptonic jet model, with radiation feedback, in the internal shock model scenario. We calculate accurate EC-scattered high-energy spectra produced by relativistic electrons throughout the Thomson and Klein-Nishina regimes. We explore the effects of varying the contribution of (1) a thermal Shakura-Sunyaev accretion disk, (2) a spherically symmetric shell of broad-line clouds, the broad line region (BLR), and (3) a hot infrared emitting dusty torus (DT), on the resultant seed photon fields. We let the system evolve to beyond the BLR and within the DT and study the manifestation of the varying target photon fields on the simulated SED and lightcurves of a typical blazar. The calculations of broadband spectra include effects of gamma-gamma absorption as gamma-rays propagate through the photon pool present inside the jet due to synchrotron and inverse Compton processes, but neglect gamma-gamma absorption by the BLR and DT photon fields outside the jet. Thus, our account of gamma-gamma absorption is a lower limit to this effect. Here, we focus on studying the impact of parameters relevant for EC processes on high-energy (HE) emission of blazars.

preprint2014arXiv

Synchrotron Polarization in Blazars

We present a detailed analysis of time- and energy-dependent synchrotron polarization signatures in a shock-in-jet model for gamma-ray blazars. Our calculations employ a full 3D radiation transfer code, assuming a helical magnetic field throughout the jet. The code considers synchrotron emission from an ordered magnetic field, and takes into account all light-travel-time and other relevant geometric effects, while the relevant synchrotron self-Compton and external Compton effects are taken care of with the 2D MCFP code. We consider several possible mechanisms through which a relativistic shock propagating through the jet may affect the jet plasma to produce a synchrotron and high-energy flare. Most plausibly, the shock is expected to lead to a compression of the magnetic field, increasing the toroidal field component and thereby changing the direction of the magnetic field in the region affected by the shock. We find that such a scenario leads to correlated synchrotron + SSC flaring, associated with substantial variability in the synchrotron polarization percentage and position angle. Most importantly, this scenario naturally explains large PA rotations by > 180 deg., as observed in connection with gamma-ray flares in several blazars, without the need for bent or helical jet trajectories or other non-axisymmetric jet features.

preprint2014arXiv

The Polarization Dependence of Gamma-Gamma Absorption - Implications for Gamma-Ray Bursts and Blazars

This paper presents an analysis of the dependence of the opacity for high-energy gamma-rays to gamma-gamma absorption by low-energy photons, on the polarization of the gamma-ray and target photons. This process has so far only been considered using the polarization-averaged gamma-gamma absorption cross section. It is demonstrated that in the case of polarized gamma-ray emission, subject to source-intrinsic gamma-gamma absorption by polarized target photons, this may lead to a slight over-estimation of the gamma-gamma opacity by up to ~ 10 % in the case of a perfectly ordered magnetic field. Thus, for realistic astrophysical scenarios with partially ordered magnetic fields, the use of the polarization-averaged gamma-gamma cross section is justified for practical purposes, such as estimates of minimum Doppler factors inferred for gamma-ray bursts and blazars, based on gamma-gamma transparency arguments, and this paper quantifies the small error incurred by the unpolarized-radiation approximation. Furthermore, it is shown that polarization-dependent gamma-gamma absorption of initially polarized gamma-rays can lead to a slight increase in the polarization beyond the spectral break caused by gamma-gamma absorption, to an amount distinctly different from the change in polarization expected if the same spectral break was produced by a break in the underlying electron distribution. This may serve as a diagnostic of whether gamma-gamma absorption is relevant in sources such as gamma-ray bursts and blazars where the gamma-ray emission may be intrinsically highly polarized.

preprint2014arXiv

Time Dependent Leptonic Modeling of Fermi II Processes in the Jets of Flat Spectrum Radio Quasars

In this paper, we discuss the light-curve features of various flaring scenarios in a time-dependent leptonic model for low-frequency-peaked blazars. The quasar 3C273 is used as an illustrative example. Our code takes into account Fermi-II acceleration and all relevant electron cooling terms, including the external radiation fields generally found to be important in the modeling of the SEDs of FSRQs, as well as synchrotron self absorption and gamma-gamma pair-production. General parameters are constrained through a fit to the average spectral energy distribution (SED) of the blazar by numerically solving the time-dependent Fokker-Planck equation for the electron evolution in a steady-state situation. We then apply perturbations to several input parameters (magnetic field, particle injection luminosity, acceleration time scale) to simulate flaring events and compute time-dependent SEDs and light curves in representative energy bands (radio, optical, X-rays, gamma-rays). Time lags between different bands are evaluated using a discrete cross correlation analysis. We find that Fermi-II acceleration has a significant effect on the distributions and that flaring events caused by increased acceleration efficiency of the Fermi II process will produce a correlation between the radio, optical and gamma-ray bandpasses, but an anti-correlation between these three bandpasses and the X-ray band, with the X-rays lagging behind the variations in other bands by up to several hours.

preprint2013arXiv

Magnetic Field Generation and Particle Energization at Relativistic Shear Boundaries in Collisionless Electron-Positron Plasmas

Using 2.5-dimensional Particle-in-Cell simulations, we study the kinetic physics of relativistic shear flow boundary in collisionless electron-positron (e+e-) plasmas. We find efficient magnetic field generation and particle energization at the shear boundary, driven by streaming instabilities across the shear interface and sustained by the shear flow. Nonthermal, anisotropic high-energy particles are accelerated across field lines to produce a power-law tail, truncated at energies below the shear Lorentz factor. These results have important implications for the dissipation and radiation of jets in blazars, gamma-ray bursts and other relativistic outflows.

preprint2013arXiv

Multiwavelength Probes of the Environs of Relativistic Shocks in Blazar Jets

Diffusive shock acceleration (DSA) at relativistic shocks is likely to be an important acceleration mechanism in various astrophysical jet sources, including radio-loud AGN. An important recent development for blazar science is the ability of Fermi-LAT data to pin down the power-law index of the high energy portion of emission in these sources, and therefore also the index of the underlying non-thermal particle population. This diagnostic potential was not possible prior to Fermi launch, when gamma-ray information was dominated by the highly-absorbed TeV band. This paper highlights how multiwavelength spectra including X-ray band and Fermi data can be used to probe diffusive acceleration in relativistic, oblique, MHD shocks in blazar jets. The spectral index of the non-thermal particle distributions resulting from Monte Carlo simulations of DSA, and the fraction of thermal particles accelerated to non-thermal energies, depend sensitively on the particles' mean free path scale, and also on the magnetic field obliquity to the shock normal. We investigate self-consistently the radiative synchrotron/Compton signatures of the resulting thermal and non-thermal particle distributions. Important constraints on the frequency of particle scattering and the level of field turbulence are identified for the blazar AO 0235+164. The possible interpretation that turbulence levels decline with remoteness from jet shocks, and a significant role for non-gyroresonant diffusion, are discussed.

preprint2013arXiv

Relativistic Positron-Electron-Ion Shear Flows and Application to Gamma-Ray Bursts

We present Particle-in-Cell simulation results of relativistic shear flows for hybrid positron-electron-ion plasmas and compare to those for pure e+e- and pure e-ion plasmas. Among the three types of relativistic shear flows, we find that only hybrid shear flow is able to energize the electrons to form a high-energy spectral peak plus a hard power-law tail. Such electron spectra are needed to model the observational properties of gamma-ray bursts.

preprint2013arXiv

X-Ray and Gamma-Ray Polarization in Leptonic and Hadronic Jet Models of Blazars

We present a theoretical analysis of the expected X-ray and gamma-ray polarization signatures resulting from synchrotron self-Compton emission in leptonic models, compared to the polarization signatures from proton synchrotron and cascade synchrotron emission in hadronic models for blazars. Source parameters resulting from detailed spectral-energy-distribution modeling are used to calculate photon-energy-dependent upper limits on the degree of polarization, assuming a perfectly organized, mono-directional magnetic field. In low-synchrotron-peaked blazars, hadronic models exhibit substantially higher maximum degrees of X-ray and gamma-ray polarization than leptonic models, which may be within reach for existing X-ray and gamma-ray polarimeters. In high-synchrotron-peaked blazars (with electron-synchrotron-dominated X-ray emission), leptonic and hadronic models predict the same degree of X-ray polarization, but substantially higher maximum gamma-ray polarization in hadronic models than leptonic ones. These predictions are particularly relevant in view of the new generation of balloon-borne X-ray polarimeters (and possibly GEMS, if revived), and the ability of Fermi-LAT to measure gamma-ray polarization at < 200 MeV. We suggest observational strategies combining optical, X-ray, gamma-ray polarimetry to determine the degree of ordering of the magnetic field and to distinguish between leptonic and hadronic high-energy emission.

preprint2012arXiv

An Angle Dependent SSC Model for Relativistic Jet Sources

We report on the development of a numerical code to calculate the angle-dependent synchrotron + synchrotron self-Compton radiation from relativistic jet sources with partially ordered magnetic fields and anisotropic particle distributions. Using a multi-zone radiation transfer approach, we can simulate magnetic-field configurations ranging from perfectly ordered (unidirectional) to randomly oriented (tangled). We demonstrate that synchrotron self-Compton model fits to the spectral energy distributions (SEDs) of extragalactic jet sources may be possible with a wide range of magnetic-field values, depending on their orientation with respect to the jet axis and the observer. This is illustrated with the example of a spectral fit to the SED of Mrk~421 from multi-wavelength observations in 2006, where acceptable fits are possible with magnetic-field values varying within a range of an order of magnitude for different degrees of B-field alignment and orientation.

preprint2012arXiv

Modeling the Spectral Energy Distributions and Variability of Blazars

In this review, recent progress in theoretical models for the broadband (radio through gamma-ray) emission from blazars are summarized. The salient features of both leptonic and hadronic models are reviewed. I present sample modeling results of spectral energy distributions (SEDs) of different types of Fermi-detected blazars along the traditional blazar sequence, using both types of models. In many cases, the SEDs of high-frequency peaked blazars (HBLs) have been found to be well represented by simple synchrotron + synchrotron self-Compton (SSC) models. However, a few HBLs recently discovered as very-high-energy (VHE) gamma-ray emitters by VERITAS are actually better represented by either external-Compton or hadronic models. Often, spectral modeling with time-independent single-zone models alone is not sufficient to constrain models, as both leptonic and lepto-hadronic models are able to provide acceptable fits to the overall SED. This degeneracy can be lifted by considering further constraints from spectral variability. Recent developments of time-dependent and inhomogeneous blazar models will be discussed, including detailed numerical simulations as well as a semi-analytical approach to the time-dependent radiation signatures of shock-in-jet models.

preprint2012arXiv

Simultaneous B'V'R' Monitoring of BL Lacertae Object S5~0716+714 and Detection of Inter-Band Time Delay

We present the results of our optical monitoring of the BL Lac object S5 0716+714 on seven nights in 2006 December. The monitoring was carried out simultaneously at three optical wavelengths with a novel photometric system. The object did not show large-amplitude internight variations during this period. Intranight variations were observed on four nights and probably on one more. Strong bluer-when-brighter chromatism was detected on both intranight and internight timescales. The intranight variation amplitude decreases in the wavelength sequence of B', R', and V'. Cross correlation analyses revealed that the variability at the $B'$ and $V'$ bands lead that at the $R'$ band by about 30 minutes on one night.

preprint2012arXiv

Synchrotron Emission from Pair Cascades in AGN Environments

Recent detections of very-high-energy (VHE, E > 100 GeV) gamma-ray blazars which do not belong to the high frequency peaked BL Lac (HBL) class, suggest that gamma-gamma absorption and pair cascading might occur in those objects. In the presence of even weak magnetic fields, these Compton-supported pair cascades will be deflected and contribute to the Fermi gamma-ray flux of radio galaxies. We demonstrate that, in this scenario, the magnetic field can not be determined from a fit of the cascade emission to the gamma-ray spectrum alone, and the degeneracy can only be lifted if the synchrotron emission from the cascades is observed as well. We illustrate this fact with the example of NGC 1275. We point out that the cascade synchrotron emission may produce spectral features reminiscent of the big blue bump observed in the spectral energy distributions of several blazars, and illustrate this idea for 3C 279.

preprint2012arXiv

Synchrotron Emission from VHE Gamma-Ray Induced Pair Cascades in AGN Environments

The discovery of very-high-energy (VHE, E > 100 GeV) gamma-ray emission from intermediate- and low-frequency peaked blazars suggests that gamma-gamma absorption and pair cascading might occur in those objects. In previous papers, we investigated the Compton emission from VHE gamma-ray induced pair cascades, deflected by moderate magnetic fields, in a largely model-independent way, and demonstrated that this emission can explain the Fermi fluxes and spectra of the radio galaxies Cen A and NGC 1275. In this paper, we describe a generalization of our Monte-Carlo cascade code to include the angle-dependent synchrotron output from the cascades, allowing for the application to situations with non-negligible magnetic fields, leading to potentially observable synchrotron signatures, but still not dominating the radiative energy loss of cascade particles. We confirm that the synchrotron radiation from the cascades in NGC 1275 and Cen A are negligible for the parameters used in our previous works. We demonstrate that the magnetic field can not be determined from a fit of the cascade emission to the gamma-ray spectrum alone, and the degeneracy can only be lifted if the synchrotron emission from the cascades is observed as well. We illustrate this fact with the example of NGC 1275. We point out that the cascade synchrotron emission may produce spectral features in the same energy range in which the big blue bump is observed in the spectral energy distributions of several blazars, and may make a non-negligible contribution to this feature. We illustrate this idea with the example of 3C 279.

preprint2012arXiv

Time-dependent simulations of emission from FSRQ PKS1510-089: multiwavelength variability of external Compton and SSC models

[abridged] We present results of modeling the SED and multiwavelength variability of the bright FSRQ PKS1510-089 with our time-dependent multizone Monte Carlo/Fokker-Planck code (Chen et al. 2001). As primary source of seed photons for inverse Compton scattering, we consider radiation from the broad line region (BLR), from the molecular torus, and the local synchrotron radiation (SSC). Different scenarios are assessed by comparing simulated light curves and SEDs with one of the best flares by PKS1510-089, in March 2009. The time-dependence of our code and its correct handling of light travel time effects allow us to fully take into account the effect of the finite size of the active region, and in turn to fully exploit the information carried by time resolved observed SEDs, increasingly available since the launch of Fermi. We confirm that the spectrum adopted for the external radiation has an important impact on the modeling of the SED, in particular for the lower energy end of the Compton component, observed in the X-ray band, which in turn is one of the most critical bands to assess the differences between EC and SSC emission. In the context of the scenario presented here, where the flaring is caused by the increase of the number of relativistic electrons ascribed to the effect of the interaction of a portion of the jet (blob) with a shock, we can not firmly discriminate the three main scenarios for gamma-ray emission. However, results show clearly the differences produced by a more realistic treatment of the emitting source in the shape of SEDs and their time variability over relevant, observable time-scales, and demonstrate the crucial importance of time-dependent multi-zone models to advance our understanding of the physics of these sources, by taking full advantage of the wealth of information offered by the high quality data of current multiwavelength campaigns.

preprint2011arXiv

Multiwavelength Spectral Studies Of Fermi-LAT Blazars

We present multiwavelength spectral analyses of two Fermi-LAT blazars, OJ 287 and 3C 279, that are part of the Boston University multiwaveband polarization program. The data have been compiled from observations with Fermi, RXTE, the VLBA, and various ground-based optical and radio telescopes. We simulate the dynamic spectral energy distributions (SEDs) within the framework of a multi-slice, time-dependent leptonic jet model for blazars, with radiation feedback, in the internal shock scenario. We use the physical jet parameters obtained from the VLBA monitoring to guide our modeling efforts. We discuss the role of intrinsic parameters and the interplay between synchrotron and inverse Compton radiation processes responsible for producing the resultant SEDs.

preprint2011arXiv

Time dependent simulations of multiwavelength variability of the blazar Mrk 421 with a Monte Carlo multi-zone code

(abridged) We present a new time-dependent multi-zone radiative transfer code and its application to study the SSC emission of Mrk 421. The code couples Fokker-Planck and Monte Carlo methods, in a 2D geometry. For the first time all the light travel time effects (LCTE) are fully considered, along with a proper treatment of Compton cooling, which depends on them. We study a set of simple scenarios where the variability is produced by injection of relativistic electrons as a `shock front' crosses the emission region. We consider emission from two components, with the second one either being pre-existing and co-spatial and participating in the evolution of the active region, or spatially separated and independent, only diluting the observed variability. Temporal and spectral results of the simulation are compared to the multiwavelength observations of Mrk 421 in March 2001. We find parameters that can adequately fit the observed SEDs and multiwavelength light curves and correlations. There remain however a few open issues, most notably: i) systematic soft intra-band X-ray lags. ii) The quadratic correlation between the TeV and X-ray flux during the flare decay has not been reproduced. These features are among those affected by the spatial extent and geometry of the source. The difficulty of producing hard X-ray lags is exacerbated by a bias towards soft lags caused by the combination of energy dependent radiative cooling time-scales and LCTE. About the second emission component, our results strongly favor the scenario where it is co-spatial and it participates in the flare evolution, suggesting that different phases of activity may occur in the same region. The cases presented in this paper represent only an initial study, and despite their limited scope they make a strong case for the need of true time-dependent and multi-zone modeling.

preprint2010arXiv

Models for the Spectral Energy Distributions and Variability of Blazars

In this review, recent progress in theoretical models for blazar emission will be summarized. The salient features of both leptonic and lepto-hadronic approaches to modeling blazar spectral energy distributions will be reviewed. I will present sample modeling results of spectral energy distributions (SEDs) of different types of blazars along the blazar sequence, including Fermi high-energy gamma-ray data, using both types of models. Special emphasis will be placed on the implications of the recent very-high-energy (VHE) gamma-ray detections of non-traditional VHE gamma-ray blazars, including intermediate and low-frequency-peaked BL Lac objects and even flat-spectrum radio quasars. Due to the featureless optical spectra of BL Lac objects, the redshifts of several BL Lacs remain unknown. I will briefly discuss possible constraints on their redshift using spectral modeling of their SED including Fermi + ground-based VHE gamma-ray data. It will be shown that in some cases, spectral modeling with time-independent single-zone models alone is not sufficient to constrain models, as both leptonic and lepto-hadronic models are able to provide acceptable fits to the overall SED. Subsequently, recent developments of time-dependent and inhomogeneous blazar models will be discussed, including detailed numerical simulations as well as a semi-analytical approach to the time-dependent radiation signatures of shock-in-jet models.

preprint2010arXiv

Time-dependent Radiation Transfer in the Internal Shock Model Scenario for Blazar Jets

We describe the time-dependent radiation transfer in blazar jets, within the internal shock model. We assume that the central engine, which consists of a black hole and an accretion disk, spews out relativistic shells of plasma with different velocity, mass, and energy. We consider a single inelastic collision between a faster (inner) and a slower (outer) moving shell. We study the dynamics of the collision and evaluate the subsequent emission of radiation via the synchrotron and synchrotron self Compton (SSC) processes after the interaction between the two shells has begun. The collision results in the formation of a forward shock (FS) and a reverse shock (RS) that convert the ordered bulk kinetic energy of the shells into magnetic field energy and accelerate the particles, which then radiate. We assume a cylindrical geometry for the emission region of the jet. We treat the self-consistent radiative transfer by taking into account the inhomogeneity in the photon density throughout the region. In this paper, we focus on understanding the effects of varying relevant input parameters on the simulated spectral energy distribution (SED) and spectral variability patterns.

preprint2010arXiv

Timing Signatures of the Internal-Shock Model for Blazars

We investigate the spectral and timing signatures of the internal-shock model for blazars. For this purpose, we develop a semi-analytical model for the time-dependent radiative output from internal shocks arising from colliding relativistic shells in a blazar jet. The emission through synchrotron and synchrotron-self Compton (SSC) radiation as well as Comptonization of an isotropic external radiation field are taken into account. We evaluate the discrete correlation function (DCF) of the model light curves in order to evaluate features of photon-energy dependent time lags and the quality of the correlation, represented by the peak value of the DCF. The almost completely analytic nature of our approach allows us to study in detail the influence of various model parameters on the resulting spectral and timing features. This paper focuses on a range of parameters in which the gamma-ray production is dominated by Comptonization of external radiation, most likely appropriate for gamma-ray bright flat-spectrum radio quasars (FSRQs) or low-frequency peaked BL Lac objects (LBLs). In most cases relevant for FSRQs and LBLs, the variability of the optical emission is highly correlated with the X-ray and high-energy (HE: > 100 MeV) gamma-ray emission. Our baseline model predicts a lead of the optical variability with respect to the higher-energy bands by 1 - 2 hours and of the HE gamma-rays before the X-rays by about 1 hour. We show that variations of certain parameters may lead to changing signs of inter-band time lags, potentially explaining the lack of persistent trends of time lags in most blazars.

preprint2010arXiv

VHE Gamma-Ray Induced Pair Cascades in Blazars and Radio Galaxies: Application to NGC 1275

Recent blazar detections by HESS, MAGIC, and VERITAS suggest that very-high-energy (VHE, E > 100 GeV) gamma-rays may be produced in most, if not all, types of blazars, including those that possess intense circumnuclear radiation fields. In this paper, we investigate the interaction of nuclear VHE gamma-rays with the circumnuclear radiation fields through gamma-gamma absorption and pair production, and the subsequent Compton-supported pair cascades. We have developed a Monte-Carlo code to follow the spatial development of the cascade in full 3-dimensional geometry, and calculate the radiative output due to the cascade as a function of viewing angle with respect to the primary VHE gamma-ray beam (presumably the jet axis of the blazar). We show that even for relatively weak magnetic fields, the cascades can be efficiently isotropized, leading to substantial off-axis cascade emission peaking in the Fermi energy range at detectable levels for nearby radio galaxies. We demonstrate that this scenario can explain the Fermi flux and spectrum of the radio galaxy NGC 1275.

preprint2010arXiv

VHE Gamma-Ray Induced Pair Cascades in the Radiation Fields of Dust Tori of AGN: Application to Cen A

The growing number of extragalactic high-energy (HE, E > 100 MeV) and very-high-energy (VHE, E > 100 GeV) gamma-ray sources that do not belong to the blazar class suggests that VHE gamma-ray production may be a common property of most radio-loud Active Galactic Nuclei (AGN). In a previous paper, we have investigated the signatures of Compton-supported pair cascades initiated by VHE gamma-ray absorption in monochromatic radiation fields, dominated by Ly-alpha line emission from the Broad Line Region. In this paper, we investigate the interaction of nuclear VHE gamma-rays with the thermal infrared radiation field from a circumnuclear dust torus. Our code follows the spatial development of the cascade in full 3-dimensional geometry. We provide a model fit to the broadband SED of the dust-rich, gamma-ray loud radio galaxy Cen A and show that typical blazar-like jet parameters may be used to model the broadband SED, if one allows for an additional cascade contribution to the Fermi gamma-ray emission.

preprint2007arXiv

X-ray Spectral Signatures of the Photon Bubble Model for Ultraluminous X-ray Sources

The nature of ultraluminous X-ray sources in nearby galaxies is one of the major open questions in modern X-ray astrophysics. One possible explanation for these objects is an inhomogeneous, radiation dominated accretion disk around a $\sim 10 M_{\odot}$ black hole -- the so-called ``photon bubble'' model. While previous studies of this model have focused primarily on its radiation-hydrodynamics aspects, in this paper, we provide an analysis of its X-ray spectral (continuum and possible edge and line) characteristics. Compton reflection between high and low density regions in the disk may provide the key to distinguishing this model from others, such as accretion onto an intermediate mass black hole. We couple a Monte Carlo/Fokker-Planck radiation transport code with the XSTAR code for reflection to simulate the photon spectra produced in a photon bubble model for ULXs. We find that reflection components tend to be very weak and in most cases not observable, and make predictions for the shape of the high-energy Comptonizing spectra. In many cases the Comptonization dominates the spectra even down to $\sim$ a few keV. In one simulation, a $\sim 9 \kev$ feature was found, which may be considered a signature of photon bubbles in ULXs; furthermore, we make predictions of high energy power-laws which may be observed by future instruments.

preprint2006arXiv

A Time Dependent Leptonic Model for Microquasar Jets: Application to LSI 61 303

The Galactic high-mass X-ray binary and jet source (microquasar) LSI +61 303 has recently been detected at TeV gamma-ray energies by the MAGIC telescope. We have applied a time-dependent leptonic jet model to the broadband spectral energy distribution and suggested (though not unambiguously detected) orbital modulation of the very high energy gamma-ray emission of this source. Our model takes into account time dependent electron injection and acceleration, and the adiabatic and radiative cooling of non-thermal electrons. It includes synchrotron, synchrotron self-Compton and external inverse Compton (with seed photons from the companion star), as well as gamma-gamma absorption of gamma-rays by starlight photons. The model can successfully reproduce the available multiwavelength observational data. Our best fit to the SED indicates that a magnetic field of B_0 ~ 5 X 10^3 G at ~ 10^3 R_g is required, and electrons need to be accelerated out to TeV energies (gamma_2 = 10^6) with a nonthermal injection spectrum with a spectral index of q = 1.7, indicating the operation of acceleration mechanisms beyond the standard first-order Fermi mechanism at relativistic or non-relativistic shocks. The orbital modulation of the VHE gamma-ray emission can be explained solely by the geometrical effect of changes in the relative orientation of the stellar companion with respect to the compact object and jet as it impacts the position and depth of the gamma-gamma absorption trough. Such a scenario predicts a trend of spectral hardening during VHE gamma-ray low orbital phases.