Probing dissipation mechanisms in BL Lac jets through X-ray polarimetry

The dissipation of energy flux in blazar jets plays a key role in the acceleration of relativistic particles. Two possibilities are commonly considered for the dissipation processes, magnetic reconnection -- possibly triggered by instabilities in magnetically-dominated jets -- , or shocks -- for weakly magnetized flows. We consider the polarimetric features expected for the two scenarios analyzing the results of state-of-the-art simulations. For the magnetic reconnection scenario we conclude, using results from global relativistic MHD simulations, that the emission likely occurs in turbulent regions with unstructured magnetic fields, although the simulations do not allow us to draw firm conclusions. On the other hand, with local particle-in-cell simulations we show that, for shocks with a magnetic field geometry suitable for particle acceleration, the self-generated magnetic field at the shock front is predominantly orthogonal to the shock normal and becomes quasi-parallel downstream. Based on this result we develop a simplified model to calculate the frequency-dependent degree of polarization, assuming that high-energy particles are injected at the shock and cool downstream. We apply our results to HBLs, blazars with the maximum of their synchrotron output at UV-soft X-ray energies. While in the optical band the predicted degree of polarization is low, in the X-ray emission it can ideally reach 50\%, especially during active/flaring states. The comparison between measurements in the optical and in the X-ray band made during active states (feasible with the planned {\it IXPE} satellite) are expected to provide valuable constraints on the dissipation and acceleration processes.Comment: 9 pages, 6 figures, accepted for publication by MNRA

Reduced magnetohydrodynamic theory of oblique plasmoid instabilities

The three-dimensional nature of plasmoid instabilities is studied using the reduced magnetohydrodynamic equations. For a Harris equilibrium with guide field, represented by \vc{B}_o = B_{po} \tanh (x/\lambda) \hat{y} + B_{zo} \hat{z}, a spectrum of modes are unstable at multiple resonant surfaces in the current sheet, rather than just the null surface of the polodial field $B_{yo} (x) = B_{po} \tanh (x/\lambda)$, which is the only resonant surface in 2D or in the absence of a guide field. Here $B_{po}$ is the asymptotic value of the equilibrium poloidal field, $B_{zo}$ is the constant equilibrium guide field, and $\lambda$ is the current sheet width. Plasmoids on each resonant surface have a unique angle of obliquity $\theta \equiv \arctan(k_z/k_y)$. The resonant surface location for angle $\theta$ is x_s = - \lambda \arctanh (\tan \theta B_{zo}/B_{po}), and the existence of a resonant surface requires $|\theta| < \arctan (B_{po} / B_{zo})$. The most unstable angle is oblique, i.e. $\theta \neq 0$ and $x_s \neq 0$, in the constant-$\psi$ regime, but parallel, i.e. $\theta = 0$ and $x_s = 0$, in the nonconstant-$\psi$ regime. For a fixed angle of obliquity, the most unstable wavenumber lies at the intersection of the constant-$\psi$ and nonconstant-$\psi$ regimes. The growth rate of this mode is $\gamma_{\textrm{max}}/\Gamma_o \simeq S_L^{1/4} (1-\mu^4)^{1/2}$, in which $\Gamma_o = V_A/L$, $V_A$ is the Alfv\'{e}n speed, $L$ is the current sheet length, and $S_L$ is the Lundquist number. The number of plasmoids scales as $N \sim S_L^{3/8} (1-\mu^2)^{-1/4} (1 + \mu^2)^{3/4}$.Comment: 9 pages, 8 figures, to be published in Physics of Plasma

Interpretation of the I-Regime and transport associated with relevant heavy particle modes

The excitation of a novel kind of heavy particle [1, 2] mode at the edge of the plasma column is considered as the signature of the I-con nement Regime [3{7]. The outward transport of impurities produced by this mode is in fact consistent with the observed expulsion of them from the main body of the plasma column (a high degree of plasma purity is a necessary feature for fusion burning plasmas capable of approaching ignition). Moreover, the theoretically predicted mode phase velocity, in the direction of the electron diamagnetic velocity, has been con rmed by relevant experimental analyses [8] of the excited uctuations (around 200 kHz). The plasma \spontaneous rotation" in the direction of the ion diamagnetic velocity is also consistent, according to the Accretion Theory [9] of this phenomenon, with the direction of the mode phase velocity. Another feature of the mode that predicted by the theory is that the I-Regime exhibits a knee of the ion temperature at the edge of the plasma column but not one of the particle density as the mode excitation factor is the relative main ion temperature gradient exceeding the local relative density gradient. The net plasma current density appearing in the saturation stage of the relevant instability, where the induced particle and energy uxes are drastically reduced, is associated with the signi cant amplitudes of the poloidal magnetic eld uctuations [6, 7] observed to accompany the density uctuations. The theoretical implications of the signi cant electron temperature uctuations [10] observed are discussed.United States. Dept. of Energ

Time dependent numerical model for the emission of radiation from relativistic plasma

We describe a numerical model constructed for the study of the emission of radiation from relativistic plasma under conditions characteristic, e.g., to gamma-ray bursts (GRB's) and active galactic nuclei (AGN's). The model solves self consistently the kinetic equations for e^\pm and photons, describing cyclo-synchrotron emission, direct Compton and inverse Compton scattering, pair production and annihilation, including the evolution of high energy electromagnetic cascades. The code allows calculations over a wide range of particle energies, spanning more than 15 orders of magnitude in energy and time scales. Our unique algorithm, which enables to follow the particle distributions over a wide energy range, allows to accurately derive spectra at high energies, >100 \TeV. We present the kinetic equations that are being solved, detailed description of the equations describing the various physical processes, the solution method, and several examples of numerical results. Excellent agreement with analytical results of the synchrotron-SSC model is found for parameter space regions in which this approximation is valid, and several examples are presented of calculations for parameter space regions where analytic results are not available.Comment: Minor changes; References added, discussion on observational status added. Accepted for publication in Ap.

Iron K Lines from Gamma Ray Bursts

We present models for reprocessing of an intense flux of X-rays and gamma rays expected in the vicinity of gamma ray burst sources. We consider the transfer and reprocessing of the energetic photons into observable features in the X-ray band, notably the K lines of iron. Our models are based on the assumption that the gas is sufficiently dense to allow the microphysical processes to be in a steady state, thus allowing efficient line emission with modest reprocessing mass and elemental abundances ranging from solar to moderately enriched. We show that the reprocessing is enhanced by down-Comptonization of photons whose energy would otherwise be too high to absorb on iron, and that pair production can have an effect on enhancing the line production. Both "distant" reprocessors such as supernova or wind remnants and "nearby" reprocessors such as outer stellar envelopes can reproduce the observed line fluxes with Fe abundances 30-100 times above solar, depending on the incidence angle. The high incidence angles required arise naturally only in nearby models, which for plausible values can reach Fe line to continuum ratios close to the reported values.Comment: 37 pages, 10 figures. Ap. J in pres

Plasmas and Controlled Nuclear Fusion

Contains report on one research project.U.S. Atomic Energy Commission (Contract AT(30-1)-3980