Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new user3D PIC Simulations for Relativistic Jets with a Toroidal Magnetic Field
120
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The properties of relativistic jets, their interaction with the ambient environment, and particle acceleration due to kinetic instabilities are studied self-consistently with Particle-in-Cell simulations. An important key issue is how a toroidal magnetic field affects the evolution of an e$^{pm}$ and an e$^{-}$ - p$^{+}$ jet, how kinetic instabilities such as the Weibel instability (WI), the mushroom instability (MI) and the kinetic Kelvin-Helmholtz instability (kKHI) are excited, and how such instabilities contribute to particle acceleration. We show that WI, MI and kKHI excited at the linear stage, generate a quasi-steady $x$-component of electric field which accelerates and decelerates electrons. In this work, we use a new jet injection scheme where an electric current is self-consistently generated at the jet orifice by the jet particles. We inject both e$^{pm}$ and e$^{-}$ - p$^{+}$ jets with a toroidal magnetic field (with a top-hat jet density profile) and for a sufficiently long time in order to examine the non-linear effects of the jet evolution. Despite the weakness of the initial magnetic field, we observe significant differences in the structure of the strong electromagnetic fields that are driven by the kinetic instabilities. We find that different jet compositions present different strongly excited instability modes. The magnetic field in the non-linear stage generated by different instabilities becomes dissipated and reorganized into a new topology. The 3-dimensional magnetic field topology indicates possible reconnection sites and the accelerated particles are significantly accelerated in the non-linear stage by the dissipation of the magnetic field and/or reconnection. This study will shed further light on the nature of astrophysical relativistic magnetized jet phenomena.
rate research
Read More
Relativistic jets naturally occur in astrophysical systems that involve accretion onto compact objects, such as core collapse of massive stars in gamma-ray bursts (GRBs) and accretion onto supermassive black holes in active galactic nuclei (AGN). It is generally accepted that these jets are powered electromagnetically, by the magnetised rotation of a central compact object. However, how they produce the observed emission and survive the propagation for many orders of magnitude in distance without being disrupted by current-driven non-axisymmetric instabilities is the subject of active debate. We carry out time-dependent 3D relativistic magnetohydrodynamic simulations of relativistic, Poynting flux dominated jets. The jets are launched self-consistently by the rotation of a strongly magnetised central compact object. This determines the natural degree of azimuthal magnetic field winding, a crucial factor that controls jet stability. We find that the jets are susceptible to two types of instability: (i) a global, external kink mode that grows on long time scales and causes the jets to bodily bend sideways. Whereas this mode does not cause jet disruption over the simulated distances, it substantially reduces jet propagation speed. We show, via an analytic model, that the growth of the external kink mode depends on the slope of the ambient medium density profile. In flat density distributions characteristic of galactic cores, an AGN jet may stall, whereas in stellar envelopes the external kink weakens as the jet propagates outward; (ii) a local, internal kink mode that grows over short time scales and causes small-angle magnetic reconnection and conversion of about half of jet electromagnetic energy flux into heat. Based on the robustness and energetics of the internal kink mode, we suggest that this instability is the main dissipation mechanism responsible for powering GRB prompt emission.
Accreting black holes launch relativistic collimated jets, across many decades in luminosity and mass, suggesting the jet launching mechanism is universal, robust and scale-free. Theoretical models and general relativistic magnetohydrodynamic (GRMHD) simulations indicate that the key jet-making ingredient is large-scale poloidal magnetic flux. However, its origin is uncertain, and it is unknown if it can be generated in situ or dragged inward from the ambient medium. Here, we use the GPU-accelerated GRMHD code HAMR to study global 3D black hole accretion at unusually high resolutions more typical of local shearing box simulations. We demonstrate that accretion disc turbulence in a radially-extended accretion disc can generate large-scale poloidal magnetic flux in situ, even when starting from a purely toroidal magnetic field. The flux accumulates around the black hole till it becomes dynamically-important, leads to a magnetically arrested disc (MAD), and launches relativistic jets that are more powerful than the accretion flow. The jet power exceeds that of previous GRMHD toroidal field simulations by a factor of 10,000. The jets do not show significant kink or pinch instabilities, accelerate to $gamma sim 10$ over 3 decades in distance, and follow a collimation profile similar to the observed M87 jet.
The Particle-In-Cell (PIC) method has been developed by Oscar Buneman, Charles Birdsall, Roger W. Hockney, and John Dawson in the 1950s and, with the advances of computing power, has been further developed for several fields such as astrophysical, magnetospheric as well as solar plasmas and recently also for atmospheric and laser-plasma physics. Currently more than 15 semi-public PIC codes are available which we discuss in this review. Its applications have grown extensively with increasing computing power available on high performance computing facilities around the world. These systems allow the study of various topics of astrophysical plasmas, such as magnetic reconnection, pulsars and black hole magnetosphere, non-relativistic and relativistic shocks, relativistic jets, and laser-plasma physics. We review a plethora of astrophysical phenomena such as relativistic jets, instabilities, magnetic reconnection, pulsars, as well as PIC simulations of laser-plasma physics (until 2021) emphasizing the physics involved in the simulations. Finally, we give an outlook of the future simulations of jets associated to neutron stars, black holes and their merging and discuss the future of PIC simulations in the light of petascale and exascale computing.
We report the results of a 3D particle-in-cell (PIC) simulation carried out to study the early-stage evolution of the shock formed when an unmagnetized relativistic jet interacts with an ambient electron-ion plasma. Full-shock structures associated with the interaction are observed in the ambient frame. When open boundaries are employed in the direction of the jet; the forward shock is seen as a hybrid structure consisting of an electrostatic shock combined with a double layer, while the reverse shock is seen as a double layer. The ambient ions show two distinct features across the forward shock: a population penetrating into the shocked region from the precursor region and an accelerated population escaping from the shocked region into the precursor region. This behavior is a signature of a combination of an electrostatic shock and a double layer. Jet electrons are seen to be electrostatically trapped between the forward and reverse shock structures showing a ring-like distribution in a phase-space plot, while ambient electrons are thermalized and become essentially isotropic in the shocked region. The magnetic energy density grows to a few percent of the jet kinetic energy density at both the forward and the reverse shock transition layers in a rather short time scale. We see little disturbance of the jet ions over this time scale.
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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
