Current-driven Kink Research Articles

The Detection of Possible Transient Quasiperiodic Oscillations in the γ-Ray Light Curve of PKS 0244-470 and 4C+38.41

The continuous monitoring capability of Fermi-LAT has enabled the exploration of quasiperiodic oscillations (QPOs) in the γ-ray light curve of blazars that has given a new perspective to probe these sources over a wide range of timescales. We report the presence of transient QPOs in the long-term γ-ray light curve of blazars PKS 0244-470 and 4C +38.41. We first identified different flux states using the Bayesian block algorithm and then explored the possible transient QPOs in the segments of each flux phase where the flux level changes over fairly regular intervals. Combining this with the source’s intrinsic variance, we identified two flux phases for PKS 0244-470: one activity (AP-1) and one quiescent phase (QP-1). For 4C+38.41, we similarly identified four activity (AP-1, AP-2, AP-3, and AP-4) and two quiescent (QP-1 and QP-2) phases. The AP-1 phase of PKS 0244-470 shows QPO of ∼225 days persisting for eight cycles (∼4.1σ). In 4C+38.41, AP-1 and AP-2 phases show QPO-like behavior of ∼110 days and ∼60 days, respectively, persisting for five cycles. In AP-3, we identified three subphases, and all show a ∼7 day scale possible recurrent rise with five complete cycles, while in QP-1, we could identify two subphases (Q1 and Q2). The Q1 phase shows a period of ∼104 days with six complete cycles. The Q2 phase also shows QPO but with only ∼3.7 cycles. We discuss the possible origin and argue that the current-driven kink instability and curved jet model seem to be the most likely causes for shorter and longer QPOs.

Rapid quasi-periodic oscillations in the relativistic jet of BL Lacertae.

Blazars are active galactic nuclei (AGN) with relativistic jets whose non-thermal radiation is extremely variable on various timescales1-3. This variability seems mostly random, although some quasi-periodic oscillations (QPOs), implying systematic processes, have been reported in blazars and other AGN. QPOs with timescales of days or hours are especially rare4 in AGN and their nature is highly debated, explained by emitting plasma moving helically inside the jet5, plasma instabilities6,7 or orbital motion in an accretion disc7,8. Here we report results of intense optical and γ-ray flux monitoring of BL Lacertae (BL Lac) during a dramatic outburst in 2020 (ref. 9). BL Lac, the prototype of a subclass of blazars10, is powered by a 1.7 × 108 MSun (ref. 11) black hole in an elliptical galaxy (distance = 313 megaparsecs (ref. 12)). Our observations show QPOs of optical flux and linear polarization, and γ-ray flux, with cycles as short as approximately 13 h during the highest state of the outburst. The QPO properties match the expectations of current-driven kink instabilities6 near a recollimation shock about 5 parsecs (pc) from the black hole in the wake of an apparent superluminal feature moving down the jet. Such a kink is apparent in a microwave Very Long Baseline Array (VLBA) image.

A simple analytical model of magnetic jets

ABSTRACT We propose a simple analytical jet model of magnetic jets, in which radially averaged profiles of main physical quantities are obtained based on conservation laws and some results of published general relativistic magnetohydrodynamic jet simulations. We take into account conversion of the magnetic energy flux to bulk acceleration in jets formed around rotating black holes assuming the mass continuity equation and constant jet power, which leads to the Bernoulli equation. For assumed profiles of the bulk Lorentz factor and the radius, this gives us the profile of the toroidal magnetic field component along the jet. We then consider the case where the poloidal field component is connected to a rotating black hole surrounded by an accretion disc. Our formalism then recovers the standard formula for the power extracted from a rotating black hole. We find that the poloidal field strength dominates over the toroidal one in the comoving frame up to large distances, which means that jets should be more stable to current-driven kink modes. The resulting magnetic field profiles can then be used to calculate the jet synchrotron emission.

Current-driven kink instabilities in relativistic jets: dissipation properties

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.

Deep learning based surrogate models for first-principles global simulations of fusion plasmas

The accurate identification and control of plasma instabilities is important for successful fusion experiments. First-principle simulations that can provide physics-based instability information such as the mode structure are generally not fast enough for real-time applications. In this work, a workflow has been presented to develop deep-learning based surrogate models for the first-principle simulations using the gyrokinetic toroidal code (GTC). The trained surrogate models of GTC (SGTC) can be used as physics-based fast instability simulators that run on the order of milliseconds, which fits the requirement of the real-time plasma control system. We demonstrate the feasibility of this workflow by first creating a big database from GTC systematic linear global electromagnetic simulations of the current-driven kink instabilities in DIII-D plasmas, and then developing SGTC linear internal kink instability simulators through supervised training. SGTC linear internal kink simulators demonstrate predictive capabilities for the mode instability properties including the growth rate and mode structure.

Fast Magnetic Reconnection Structures in Poynting Flux-dominated Jets

Abstract The ubiquitous relativistic jet phenomena associated with black holes play a major role in high and very-high-energy (VHE) astrophysics. In particular, observations have demonstrated that blazars show VHE emission with time variability from days to minutes (in the gigaelectronvolt and teraelectronvolt bands), implying very compact emission regions. The real mechanism of the particle acceleration process responsible for this emission is still being debated, but magnetic reconnection has lately been discussed as a strong potential candidate. In this work, we present the results of three-dimensional special relativistic magnetohydrodynamic simulations of the development of reconnection events driven by turbulence induced by current-driven kink instability along a relativistic jet. We have performed a systematic identification of all reconnection regions in the system, characterizing their local magnetic field topology and quantifying the reconnection rates. We obtained average rates of 0.051 ± 0.026 (in units of Alfvén speed), which are comparable to the predictions of the theory of turbulence-induced fast reconnection. A detailed statistical analysis also demonstrated that the fast reconnection events follow a log-normal distribution, which is a signature of its turbulent origin. To probe the robustness of our method, we have applied our results to the blazar Mrk 421. Building a synthetic light curve from the integrated magnetic reconnection power, we evaluated the time variability from a power spectral density analysis, obtaining good agreement with observations in the gigaelectronvolt band. This suggests that turbulent fast magnetic reconnection can be a possible process behind the high-energy emission variability phenomena observed in blazars.

Particle Acceleration by Relativistic Magnetic Reconnection Driven by Kink Instability Turbulence in Poynting Flux–Dominated Jets

Abstract Particle acceleration in magnetized relativistic jets still puzzles theorists. In this work, we investigate the acceleration of particles injected into a three-dimensional relativistic magnetohydrodynamical jet subject to current-driven kink (CDK) instability. We find that, once turbulence driven by CDK fully develops, achieving a nearly stationary state, the amplitude of excited wiggles along the jet spine attains maximum growth, causing disruption of the magnetic field lines and the formation of several sites of fast reconnection. Low-energy protons injected into the jet at this state experience exponential acceleration, mostly in directions parallel to the local magnetic field, up to maximum energies eV for G and eV for G. The Larmor radius of the particles attaining these energies corresponds to the size of the acceleration region (∼the diameter of the perturbed jet). There is a clear association of the accelerated particles with regions of fast reconnection. In the early nonlinear growth stage of the CDK, when there are no sites of fast reconnection yet, injected particles with initially much larger energy are accelerated by magnetic curvature drift. We have also obtained the acceleration time due to reconnection with a dependence on the particles’ energy, . The energy spectrum of the accelerated particles develops a power-law index ∼ −1.2 in the beginning, in agreement with earlier works. Our results provide a multidimensional framework for exploring this process in real systems and explain their emission patterns, especially at very high energies, and associated neutrino emission recently detected in some blazars.

Particle acceleration and the origin of the very high energy emission around black holes and relativistic jets

AbstractParticle acceleration induced by fast magnetic reconnection may help to solve current puzzles related to the interpretation of the very high energy (VHE) and neutrino emissions from AGNs and compact sources in general. Our general relativistic-MHD simulations of accretion disk-corona systems reveal the growth of turbulence driven by MHD instabilities that lead to the development of fast magnetic reconnection in the corona. In addition, our simulations of relativistic MHD jets reveal the formation of several sites of fast reconnection induced by current-driven kink turbulence. The injection of thousands of test particles in these regions causes acceleration up to energies of several PeVs, thus demonstrating the ability of this process to accelerate particles and produce VHE and neutrino emission, specially in blazars. Finally, we discuss how reconnection can also explain the observed VHE luminosity-black hole mass correlation, involving hundreds of non-blazar sources like Perseus A, and black hole binaries.

A Model of Polarisation Rotations in Blazars from Kink Instabilities in Relativistic Jets

This paper presents a simple model of polarisation rotation in optically thin relativistic jets of blazars. The model is based on the development of helical (kink) mode of current-driven instability. A possible explanation is suggested for the observational connection between polarisation rotations and optical/gamma-ray flares in blazars, if the current-driven modes are triggered by secular increases of the total jet power. The importance of intrinsic depolarisation in limiting the amplitude of coherent polarisation rotations is demonstrated. The polarisation rotation amplitude is thus very sensitive to the viewing angle, which appears to be inconsistent with the observational estimates of viewing angles in blazars showing polarisation rotations. Overall, there are serious obstacles to explaining large-amplitude polarisation rotations in blazars in terms of current-driven kink modes.

Magnetic Dissipation in Relativistic Jets

The most promising mechanisms for producing and accelerating relativistic jets, and maintaining collimated structure of relativistic jets involve magnetohydrodynamical (MHD) processes. We have investigated the magnetic dissipation mechanism in relativistic jets via relativistic MHD simulations. We found that the relativistic jets involving a helical magnetic field are unstable for the current-driven kink instability, which leads to helically distorted structure in relativistic jets. We identified the regions of high current density in filamentary current sheets, indicative of magnetic reconnection, which are associated to the kink unstable regions and correlated to the converted regions of magnetic to kinetic energies of the jets. We also found that an over-pressured relativistic jet leads to the generation of a series of stationary recollimation shocks and rarefaction structures by the nonlinear interaction of shocks and rarefaction waves. The differences in the recollimation shock structure due to the difference of the magnetic field topologies and strengths may be observable through mm-VLBI observations and space-VLBI mission.

The collimation of magnetic jets by disc winds

The collimation of a Poynting-flux dominated jet by a wind emanating from the surface of an accretion flow is computed using a semi-analytic model. The injection of the disk wind is treated as a boundary condition in the equatorial plane, and its evolution is followed by invoking a prescribed geometry of streamlines. Solutions are obtained for a wide range of disk wind parameters. It is found that jet collimation generally occurs when the total wind power exceeds about ten percents of the jet power. For moderate wind powers we find gradual collimation. For strong winds we find rapid collimation followed by focusing of the jet, after which it remains narrow over many Alfv\'en crossing times before becoming conical. We estimate that in the later case the jet's magnetic field may be dissipated by the current-driven kink instability over a distance of a few hundreds gravitational radii. We apply the model to M87 and show that the observed parabolic shape of the radio jet within the Bondi radius can be reproduced provided that the wind injection zone extends to several hundreds gravitational radii, and that its total power is about one third of the jet power. The radio spectrum can be produced by synchrotron radiation of relativistically hot, thermal electrons in the sheath flow surrounding the inner jet.

SPATIAL GROWTH OF CURRENT-DRIVEN INSTABILITY IN RELATIVISTIC ROTATING JETS AND THE SEARCH FOR MAGNETIC RECONNECTION

ABSTRACT Using the three-dimensional relativistic magnetohydrodynamic code RAISHIN, we investigated the influence of the radial density profile on the spatial development of the current-driven kink instability along magnetized rotating, relativistic jets. For the purposes of our study, we used a nonperiodic computational box, the jet flow is initially established across the computational grid, and a precessional perturbation at the inlet triggers the growth of the kink instability. We studied light and heavy jets with respect to the environment depending on the density profile. Different angular velocity amplitudes have been also tested. The results show the propagation of a helically kinked structure along the jet and a relatively stable configuration for the lighter jets. The jets appear to be collimated by the magnetic field, and the flow is accelerated owing to conversion of electromagnetic into kinetic energy. We also identify regions of high current density in filamentary current sheets, indicative of magnetic reconnection, which are associated with the kink-unstable regions and correlated with the decrease of the sigma parameter of the flow. We discuss the implications of our findings for Poynting-flux-dominated jets in connection with magnetic reconnection processes. We find that fast magnetic reconnection may be driven by the kink-instability turbulence and govern the transformation of magnetic into kinetic energy, thus providing an efficient way to power and accelerate particles in active galactic nucleus and gamma-ray-burst relativistic jets.

Modelling the kinked jet of the Crab nebula

We investigate the dynamical propagation of the South-East jet from the Crab pulsar interacting with supernova ejecta by means of three-dimensional relativistic MHD numerical simulations with the PLUTO code. The initial jet structure is set up from the inner regions of the Crab Nebula. We study the evolution of hot, relativistic hollow outflows initially carrying a purely azimuthal magnetic field. Our jet models are characterized by different choices of the outflow magnetization ($\sigma$ parameter) and the bulk Lorentz factor ($\gamma_{j}$). We show that the jet is heavily affected by the growth of current-driven kink instabilities causing considerable deflection throughout its propagation length. This behavior is partially stabilized by the combined action of larger flow velocities and/or reduced magnetic field strengths. We find that our best jet models are characterized by relatively large values of $\sigma$ ($\gtrsim 1$) and small values of $\gamma_{j}\simeq 2$. Our results are in good agreement with the recent X-ray (\textit{Chandra}) data of the Crab Nebula South-East jet indicating that the jet changes direction of propagation on a time scale of the order of few years. The 3D models presented here may have important implications in the investigation of particle acceleration in relativistic outflows.

Current-Driven Kink Instability in Magnetically Dominated Rotating Relativistic Jet

We have investigated the influence of jet rotation and differential motion on the linear and nonlinear development of the current-driven (CD) kink instability of force-free helical magnetic equilibria via threedimensional relativistic magnetohydrodynamic simulations. In this study, we follow the temporal development within a periodic computational box. Displacement of the initial helical magnetic field leads to the growth of the CD kink instability. In rotating relativistic jet case, developed helical kink structure propagates along jet axis with continuous growth of kink amplitude. The growth rate of CD kink instability does not depend on the jet rotation. The coupling of multiple unstable wavelengths is crucial to determining whether the jet is eventually disrupted in nonlinear stage. The strongly deformed magnetic field via CD kink instability may trigger of magnetic reconnection in the jet.

Current-Driven Kink Instability in Relativistic Jets

AbstractWe have investigated the development of current-driven (CD) kink instability in relativistic jets via 3D RMHD simulations. In this investigation a static force-free equilibrium helical magnetic field configuration is considered in order to study the influence of the initial configuration on the linear and nonlinear evolution of the instability. We found that the initial configuration is strongly distorted but not disrupted by the CD kink instability. The linear growth and nonlinear evolution of the CD kink instability depends moderately on the radial density profile and strongly on the magnetic pitch profile. Kink amplitude growth in the nonlinear regime for decreasing magnetic pitch leads to a slender helically twisted column wrapped by magnetic field. On the other hand, kink amplitude growth in the nonlinear regime nearly ceases for increasing magnetic pitch.

High-resolution 3D relativistic MHD simulations of jets

Relativistic magnetized jets are key elements in active galactic nuclei and in other astrophysical environments. Their structure and evolution involve a complex non-linear physics that can be approached by numerical studies only. Still, owing to a number of challenging computational aspects, only a few numerical investigations have been undertaken so far. In this paper, we present high-resolution three-dimensional numerical simulations of relativistic magnetized jets carrying an initially toroidal magnetic field. The presence of a substantial toroidal component of the field is nowadays commonly invoked and held responsible for the process of jet acceleration and collimation. We find that the typical nose cone structures, commonly observed in axisymmetric two-dimensional simulations, are not produced in the three-dimensional case. Rather, the toroidal field gives rise to strong current-driven kink instabilities leading to jet wiggling. However, it appears to be able to maintain a highly relativistic spine along its full length. By comparing low- and high-resolution simulations, we emphasize the impact of resolution on the jet dynamical properties.

A Numerical Model of Hercules A by Magnetic Tower: Jet/Lobe Transition, Wiggling, and the Magnetic Field Distribution

We apply magnetohydrodynamic (MHD) modeling to the radio galaxy Hercules A for investigating the jetdriven shock, jet/lobe transition, wiggling, and magnetic field distribution associated with this source. The model consists of magnetic tower jets in a galaxy cluster environment, which has been discussed in a series of our papers. The profile of underlying ambient gas plays an imp ortant role in jet-lobe morphology. The balance between the magnetic pressure generated by axial current and the ambient gas pressure can determine the lobe radius. The jet body is confined jointly by the external press ure and gravity inside the cluster core radius Rc, while outside Rc it expands radially to form fat lobes in a steeply decreasing ambient thermal pressure gradient. The current-carrying jets are responsible for ge nerating a strong, tightly wound helical magnetic field. This magnetic configuration will be unstable against t he current-driven kink mode and it visibly grows beyond Rc where a separation between the jet forward and return currents occurs. The reversed pinch profile of global magnetic field associated with the jet and lobes pro duces projected B-vector distributions aligned with the jet flow and the lobe edge. AGN-driven shock powered b y the expanding magnetic tower jet surrounds the jet/lobe structure and heats the ambient ICM. The lobes expand subsonically; no obvious hot spots are produced at the heads of lobes. Several key features in our MHD modeling may be qualitatively supported by the observations of Hercules A. Subject headings:galaxies:individual: Hercules A — galaxies: active — galaxies: jets — methods: numerical — MHD

Relationship of the Multimode Kink Instability to δ‐Spot Formation

We study the current-driven kink instability with a three-dimensional MHD spectral code, under conditions appropriate for the convection zone of the Sun. Our goal is to determine whether the kink instability can explain the unusual morphology of certain active regions. The characteristics of these δ-spots are unusually large tilt upon emergence, subsequent rotation, compactness during their lifetime, magnetic shear that develops along the neutral line of their vertical magnetic field, and unusually high flare activity. We find that highly twisted tubes perturbed with a single unstable kink mode evolve to a new helical equilibrium with finite amplitude. This kink would produce an active region that emerged tilted and then rotated and so is suggestive of δ-spot behavior. It would not, however, produce other characteristics of δ-spots such as compactness and shear. We then investigate the kink instability in tubes perturbed simultaneously with several unstable modes at different wavelengths. These tubes develop a concentrated kink at the point at which the modes interact constructively. We show that such a concentrated kink would first emerge highly tilted, rotate during its subsequent emergence while remaining compact, and develop strong shear along its magnetic neutral line, in agreement with δ-spot observations. In addition we find that very strong concentrated kinks develop a current sheet at which the magnetic field reconnects, causing field lines near the center of the tube to become knotted. This reconnection could be related to the high flare activity of δ-spots.

Structure of the n = 1 mode responsible for relaxation and current drive during sustainment of the SPHEX spheromak

The structure of the n = 1 mode in the SPHEX spheromak, which plays a central role in relaxation during sustainment, is investigated by analysing the measured voltage fluctuations in the central plasma column. By combining these results with a suitably defined helical magnetic flux function, the mode is found to be due to a rotating helical distortion of the open linked flux. We propose that the distortion is due to a saturated current-driven kink mode of the open flux tube. The prolongation of this `helical column' on its return around the outside of the closed flux is found to be strongly asymmetric. Previously published measurements of the Poynting flux and -profile are re-analysed in the light of these results, and implications for the mechanism of relaxation and non-inductive current drive are discussed.

Present state of the theory of the MHD equilibrium and stability of stellarator plasmas

The MHD theory of the equilibrium and stability of the plasma in a stellarator is set forth. Various ways to develop magnetic configurations of the stellarator type are examined. The basic characteristics of the devices presently in operation are reported. The method of averaging over the rapidly varying spatial variable is described. This method is used to study the plasma equilibrium and also the stability with respect to current-driven kink modes and ballooning modes. The limiting plasma pressure, determined jointly by the conditions for plasma equilibrium and stability, is discussed. Methods for raising this limit are described. The present state of the theory of the MHD equilibrium and stability of plasmas in stellarators is summarized. The basic problems for further development of this theory are outlined.