Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userFrom electrons to Janskys: Full stokes polarized radiative transfer in 3D relativistic particle-in-cell jet simulations
215
0
0.0
(
0
)
Added by
Nicholas MacDonald Dr.
Publication date
2021
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
The underlying plasma composition of relativistic extragalactic jets remains largely unknown. Relativistic magnetohydrodynamic (RMHD) models are able to reproduce many of the observed macroscopic features of these outflows. The nonthermal synchrotron emission detected by very long baseline interferometric (VLBI) arrays, however, is a by-product of the kinetic-scale physics occurring within the jet, physics that is not modeled directly in most RMHD codes. This paper attempts to discern the radiative differences between distinct plasma compositions within relativistic jets using small-scale 3D relativistic particle-in-cell (PIC) simulations. We generate full Stokes imaging of two PIC jet simulations, one in which the jet is composed of an electron-proton ($e^{-}$-$p^{+}$) plasma (i.e., a normal plasma jet), and the other in which the jet is composed of an electron-positron ($e^{-}$-$e^{+}$) plasma (i.e., a pair plasma jet). We examined the differences in the morphology and intensity of the linear polarization (LP) and circular polarization (CP) emanating from these two jet simulations. We find that the fractional level of CP emanating from the $e^{-}$-$p^{+}$ plasma jet is orders of magnitude larger than the level emanating from an $e^{-}$-$e^{+}$ plasma jet of a similar speed and magnetic field strength. In addition, we find that the morphology of both the linearly and circularly polarized synchrotron emission is distinct between the two jet compositions. We also demonstrate the importance of slow-light interpolation and we highlight the effect that a finite light-crossing time has on the resultant polarization when ray-tracing through relativistic plasma.
rate research
Read More
We compare maps of scattering polarization signals obtained from three-dimensional (3D) radiation transfer calculations in a magneto-convection model of the solar atmosphere using formal solvers based on the short characteristics (SC) and the long characteristics (LC) methods. The SC method requires less computational work, but it is known to introduce spatial blurring in the emergent radiation for inclined lines of sight. For polarized radiation this effect is generally more severe due to it being a signed quantity and to the sensitivity of the scattering polarization to the models inhomogeneities. We study the differences in the polarization signals of the emergent spectral line radiation calculated with such formal solvers. We take as a case study already published results of the scattering polarization in the Sr I $4607~unicode{xC5}$ line obtained with the SC method, demonstrating that in high-resolution grids it is accurate enough for that type of study. In general, the LC method is the preferred one for accurate calculations of the emergent radiation, reason why it is now one of the options in the public version of the 3D radiative transfer code PORTA.
We study the interaction of relativistic jets with their environment, using 3-dimensional relativistic particle-in-cell simulations for two cases of jet composition: (i) electron-proton ($e^{-}-p^{+}$) and (ii) electron-positron ($e^{pm}$) plasmas containing helical magnetic fields. We have performed simulations of global jets containing helical magnetic fields in order to examine how helical magnetic fields affect kinetic instabilities such as the Weibel instability, the kinetic Kelvin-Helmholtz instability and the Mushroom instability. We have found that these kinetic instabilities are suppressed and new types of instabilities can grow. For the $e^{-}-p^{+}$ jet, a recollimation-like instability occurs and jet electrons are strongly perturbed, whereas for the $e^{pm}$ jet, a recollimation-like instability occurs at early times followed by kinetic instability and the general structure is similar to a simulation without a helical magnetic field. We plan to perform further simulations using much larger systems to confirm these new findings.
Accreting supermassive black holes are sources of polarized radiation that propagates through highly curved spacetime before reaching the observer. In order to help interpret observations of such polarized emission, accurate and efficient numerical schemes for polarized radiative transfer in curved spacetime are needed. In this manuscript we extend our publicly available radiative transfer code RAPTOR to include polarization. We provide a brief review of different codes and methods for covariant polarized radiative transfer available in the literature and existing codes, and present an efficient new scheme. For the spacetime-propagation aspect of the computation, we develop a compact, Lorentz-invariant representation of a polarized ray. For the plasma-propagation aspect of the computation, we perform a formal analysis of the stiffness of the polarized radiative-transfer equation with respect to our explicit integrator, and develop a hybrid integration scheme that switches to an implicit integrator in case of stiffness, in order to solve the equation with optimal speed and accuracy for all possible values of the local optical/Faraday thickness of the plasma. We perform a comprehensive code verification by solving a number of well-known test problems using RAPTOR and comparing its output to exact solutions. We also demonstrate convergence with existing polarized radiative-transfer codes in the context of complex astrophysical problems. RAPTOR is capable of performing polarized radiative transfer in arbitrary, highly curved spacetimes. This capability is crucial for interpreting polarized observations of accreting black holes, which can yield information about the magnetic-field configuration in such accretion flows. The efficient formalism implemented in RAPTOR is computationally light and conceptually simple. The code is publicly available.
Numerical algorithms to load relativistic Maxwell distributions in particle-in-cell (PIC) and Monte-Carlo simulations are presented. For stationary relativistic Maxwellian, the inverse transform method and the Sobol algorithm are reviewed. To boost particles to obtain relativistic shifted-Maxwellian, two rejection methods are proposed in a physically transparent manner. Their acceptance efficiencies are ${approx}50%$ for generic cases and $100%$ for symmetric distributions. They can be combined with arbitrary base algorithms.
We investigate the temporal evolution of an axisymmetric magnetosphere around a rapidly rotating, stellar-mass black hole, applying a two-dimensional particle-in-cell simulation scheme. Adopting a homogeneous pair production, and assuming that the mass accretion rate is much less than the Eddington limit, we find that the black holes rotational energy is preferentially extracted from the middle latitudes, and that this outward energy flux exhibits an enhancement that lasts approximately 160 dynamical time scales. It is demonstrated that the magnetohydrodynamic approximations cannot be justified in such a magnetically-dominated magnetosphere, because the Ohms law completely breaks down, and because the charge-separated electron-positron plasmas are highly non-neutral. An implication is given regarding the collimation of relativistic jets.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
