ABSTRACT
We investigate with relativistic MHD (magnetohydrodynamic) simulations the dissipation physics of BL Lac jets, by studying the synchrotron polarization signatures of particles accelerated by ...the kink instability in a magnetically dominated plasma column. The non-linear stage of the kink instability generates current sheets, where particles can be efficiently accelerated via magnetic reconnection. We identify current sheets as regions where s = Jδ/B is above some predefined threshold (where B is the field strength, J is the current density, and δ is the grid scale), and assume that the particle injection efficiency scales as ∝J2. X-ray emitting particles have short cooling times, so they only probe the field geometry of their injection sites. In contrast, particles emitting in the optical band, which we follow self-consistently as they propagate away from their injection sites while cooling, sample a larger volume, and so they may be expected to produce different polarimetric signatures. We find that the degree of polarization is roughly the same between X-ray and optical bands, because even the optical emitting particles do not travel far from the current sheet where they were injected, due to lack of sufficient kink-generated turbulence. The polarization angle shows a different temporal evolution between the two bands, due to the different regions probed by X-ray and optical emitting particles. In view of the upcoming IXPE satellite, our results can help constrain whether kink-induced reconnection (as opposed to shocks) can be the source of multiwavelength emission from BL Lacs.
We describe an implementation of a particle physics module available for the PLUTO code appropriate for the dynamical evolution of a plasma consisting of a thermal fluid and a nonthermal component ...represented by relativistic charged particles or cosmic rays (CRs). While the fluid is approached using standard numerical schemes for magnetohydrodynamics, CR particles are treated kinetically using conventional Particle-In-Cell (PIC) techniques. The module can be used either to describe test-particle motion in the fluid electromagnetic field or to solve the fully coupled magnetohydrodynamics (MHD)-PIC system of equations with particle backreaction on the fluid as originally introduced by Bai et al. Particle backreaction on the fluid is included in the form of momentum-energy feedback and by introducing the CR-induced Hall term in Ohm's law. The hybrid MHD-PIC module can be employed to study CR kinetic effects on scales larger than the (ion) skin depth provided that the Larmor gyration scale is properly resolved. When applicable, this formulation avoids resolving microscopic scales, offering substantial computational savings with respect to PIC simulations. We present a fully conservative formulation that is second-order accurate in time and space, and extends to either the Runge-Kutta (RK) or the corner transport upwind time-stepping schemes (for the fluid), while a standard Boris integrator is employed for the particles. For highly energetic relativistic CRs and in order to overcome the time-step restriction, a novel subcycling strategy that retains second-order accuracy in time is presented. Numerical benchmarks and applications including Bell instability, diffusive shock acceleration, and test-particle acceleration in reconnecting layers are discussed.
We present a description of the adaptive mesh refinement (AMR) implementation of the PLUTO code for solving the equations of classical and special relativistic magnetohydrodynamics (MHD and RMHD). ...The current release exploits, in addition to the static grid version of the code, the distributed infrastructure of the CHOMBO library for multidimensional parallel computations over block-structured, adaptively refined grids. We employ a conservative finite-volume approach where primary flow quantities are discretized at the cell center in a dimensionally unsplit fashion using the Corner Transport Upwind method. Time stepping relies on a characteristic tracing step where piecewise parabolic method, weighted essentially non-oscillatory, or slope-limited linear interpolation schemes can be handily adopted. A characteristic decomposition-free version of the scheme is also illustrated. The solenoidal condition of the magnetic field is enforced by augmenting the equations with a generalized Lagrange multiplier providing propagation and damping of divergence errors through a mixed hyperbolic/parabolic explicit cleaning step. Among the novel features, we describe an extension of the scheme to include non-ideal dissipative processes, such as viscosity, resistivity, and anisotropic thermal conduction without operator splitting. Finally, we illustrate an efficient treatment of point-local, potentially stiff source terms over hierarchical nested grids by taking advantage of the adaptivity in time. Several multidimensional benchmarks and applications to problems of astrophysical relevance assess the potentiality of the AMR version of PLUTO in resolving flow features separated by large spatial and temporal disparities.
We present the results of state-of-the-art simulations of recollimation shocks induced by the interaction of a relativistic jet with an external medium, including the effect of radiative losses of ...the shocked gas. Our simulations confirm that – as suggested by earlier semi-analytical models – the post-shock pressure loss induced by radiative losses may lead to a stationary equilibrium state characterized by a very strong focusing of the flow, with the formation of quite narrow nozzles, with cross-sectional radii as small as 10-3 times the length scale of the jet. We also study the time-dependent evolution of the jet structure induced by a density perturbation injected at the flow base. The set-up and the results of the simulations are particularly relevant for the interpretation of the observed rapid variability of the γ-ray emission associated to flat spectrum radio quasars. In particular, the combined effects of jet focusing and Doppler beaming on the observed radiation make it possible to explain the sub-hour flaring events such as that observed in the flat specrum radio quasar PKS 1222+216 by MAGIC.
ABSTRACT
We carry out a linear stability analysis of a magnetized relativistic rotating cylindrical jet flow using the approximation of zero thermal pressure. We identify several modes of instability ...in the jet: Kelvin–Helmholtz, current-driven and two kinds of centrifugal-buoyancy modes – toroidal and poloidal. The Kelvin–Helmholtz mode is found at low magnetization and its growth rate depends very weakly on the pitch parameter of the background magnetic field and on rotation. The current-driven mode is found at high magnetization, the values of its growth rate and the wavenumber, corresponding to the maximum growth, increase as we decrease the pitch parameter of the background magnetic field. This mode is stabilized by rotation, especially, at high magnetization. The centrifugal-buoyancy modes, arising due to rotation, tend also to be more stable when magnetization is increased. Overall, relativistic jet flows appear to be more stable with respect to their non-relativistic counterpart.
ABSTRACT
We analyse the evolution of current-driven kink instabilities of a highly magnetized relativistic plasma column, focusing in particular on its dissipation properties. The instability ...evolution leads to the formation of thin current sheets where the magnetic energy is dissipated. We find that the total amount of dissipated magnetic energy is independent of the dissipation properties. Dissipation occurs in two stages: a peak when the instability saturates that is characterized by the formation of a helicoidal current sheet at the boundary of the deformed plasma column, followed by a weaker almost flat phase, in which turbulence develops. The detailed properties of these two phases depend on the equilibrium configuration and other parameters, in particular on the steepness of the pitch radial profile, on the presence of an external axial magnetic field and on the amount of magnetization. These results are relevant for high energy astrophysical sources, since current sheets can be the sites of magnetic reconnection where particles can be accelerated to relativistic energies and give rise to the observed radiation.
ABSTRACT
We describe a novel Godunov-type numerical method for solving the equations of resistive relativistic magnetohydrodynamics. In the proposed approach, the spatial components of both magnetic ...and electric fields are located at zone interfaces and are evolved using the constrained transport formalism. Direct application of Stokes’ theorem to Faraday’s and Ampere’s laws ensures that the resulting discretization is divergence-free for the magnetic field and charge-conserving for the electric field. Hydrodynamic variables retain, instead, the usual zone-centred representation commonly adopted in finite-volume schemes. Temporal discretization is based on Runge–Kutta implicit–explicit schemes in order to resolve the temporal scale disparity introduced by the stiff source term in Ampere’s law. The implicit step is accomplished by means of an improved and more efficient Newton–Broyden multidimensional root-finding algorithm. The explicit step relies on a multidimensional Riemann solver to compute the line-averaged electric and magnetic fields at zone edges and it employs a 1D Riemann solver at zone interfaces to update zone-centred hydrodynamic quantities. For the latter, we introduce a five-wave solver based on the frozen limit of the relaxation system whereby the solution to the Riemann problem can be decomposed into an outer Maxwell solver and an inner hydrodynamic solver. A number of numerical benchmarks demonstrate that our method is superior in stability and robustness to the more popular charge-conserving divergence cleaning approach where both primary electric and magnetic fields are zone-centred. In addition, the employment of a less diffusive Riemann solver noticeably improves the accuracy of the computations.
Context. A significant fraction of extended radio sources presents a peculiar X-shaped radio morphology: in addition to the classical double lobed structure, radio emission is also observed along a ...second axis of symmetry in the form of diffuse wings or tails. In a previous investigation we showed the existence of a connection between the radio morphology and the properties of the host galaxies. Motivated by this connection we performed two-dimensional numerical simulations showing that X-shaped radio sources may naturally form as a jet propagates along the major axis a highly elliptical density distribution, because of the fast expansion of the cocoon along the minor axis of the distribution. Aims. We intend to extend our analysis by performing three-dimensional numerical simulations and investigating the role of different parameters in determining the formation of the X-shaped morphology. Methods. The problem is addressed by numerical means, carrying out three-dimensional relativistic magnetohydrodynamic simulations of bidirectional jets propagating in a triaxial density distribution. Results. We show that only jets with power ≲ 1044 erg s-1 can give origin to an X-shaped morphology and that a misalignment of 30° between the jet axis and the major axis of the density distribution is still favourable to the formation of this kind of morphology. In addition we compute synthetic radio emission maps and polarization maps. Conclusions. In our scenario for the formation of X-shaped radio sources only low power FRII can give origin to such kind of morphology. Our synthetic emission maps show that the different observed morphologies of X-shaped sources can be the result of similar structures viewed under different perspectives.
Context. Extragalactic radio sources have been classified into two classes, Fanaroff-Riley I and II, which differ in morphology and radio power. Strongly emitting sources belong to the ...edge-brightened FR II class, and weakly emitting sources to the edge-darkened FR I class. The origin of this dichotomy is not yet fully understood. Numerical simulations are successful in generating FR II morphologies, but they fail to reproduce the diffuse structure of FR Is. Aims. By means of hydro-dynamical 3D simulations of supersonic jets, we investigate how the displayed morphologies depend on the jet parameters. Bow shocks and Mach disks at the jet head, which are probably responsible for the hot spots in the FR II sources, disappear for a jet kinetic power L sub(kin)< or = 10 super(43) ergs super(-1). This threshold compares favorably with the luminosity at which the FR I/FR II transition is observed. Methods. The problem is addressed by numerical means carrying out 3D HD simulations of supersonic jets that propagate in a non-homogeneous medium with the ambient temperature that increases with distance from the jet origin, which maintains constant pressure. Results. The jet energy in the lower power sources, instead of being deposited at the terminal shock, is gradually dissipated by the turbulence. The jets spread out while propagating, and they smoothly decelerate while mixing with the ambient medium and produce the plumes characteristic of FR I objects. Conclusions. Three-dimensionality is an essential ingredient to explore the FR I evolution becausethe properties of turbulence in two and three dimensions are very different, since there is no energy cascade to small scales in two dimensions, and two-dimensional simulations with the same parameters lead to FRII-like behavior.
We perform three-dimensional numerical simulations of relativistic (with a Lorentz factor of 10), non-magnetized jets propagating in a uniform density environment in order to study the effect of the ...entrainment and the consequent deceleration. Our simulations investigate the jet propagation inside the galaxy core, where the deceleration most likely occurs more efficiently. We compare cases with different density and pressure ratios with respect to the ambient medium and find that low density jets are efficiently decelerated and reach a quasi-steady state in which, over a length of 600 jet radii, they slow down from highly relativistic to sub-relativistic velocities. Conversely, denser jets keep highly relativistic velocities over the same length. We discuss these results in relation to the Faranoff Riley (FR) radio source classification. We infer that lower density jets can give rise to FR 0 and FR I radio sources, while higher density jets may be connected to FR II radio sources.