Relativistic Magnetohydrodynamics Research Articles

Relativistic dissipative magnetohydrodynamics from the Boltzmann equation for a two-component gas

Relativistic dissipative magnetohydrodynamics from the Boltzmann equation for a two-component gas

Evidence of a toroidal magnetic field in the core of 3C 84

The spatial scales of relativistic radio jets, probed by relativistic magneto-hydrodynamic (RMHD) jet launching simulations and by most very long baseline interferometry (VLBI) observations differ by an order of magnitude. Bridging the gap between these RMHD simulations and VLBI observations requires selecting nearby active galactic nuclei (AGN), the parsec-scale region of which can be resolved. The radio source is a nearby bright AGN fulfilling the necessary requirements: it is launching a powerful, relativistic jet powered by a central supermassive black hole, while also being very bright. Using 22\,GHz globe-spanning VLBI measurements of we studied its sub-parsec region in both total intensity and linear polarisation to explore the properties of this jet, with a linear resolution of $ We tested different simulation set-ups by altering the bulk Lorentz factor Gamma of the jet, as well as the magnetic field configuration (toroidal, poloidal helical). We confirm the persistence of a limb brightened structure, which reaches deep into the sub-parsec region. The corresponding electric vector position angles (EVPAs) follow the bulk jet flow inside but tend to be orthogonal to it near the edges. Our state-of-the-art RMHD simulations show that this geometry is consistent with a spine-sheath model, associated with a mildly relativistic flow and a toroidal magnetic field configuration.

Dynamical Properties of Magnetized Low-angular-momentum Accretion Flows around a Kerr Black Hole

An essential factor for determining the characteristics of an accretion flow is its angular momentum. According to the angular momentum of the flow, semi-analytical analysis suggests various types of accretion solutions. It is critical to test these with numerical simulations, using the most advanced framework available (general relativistic magnetohydrodynamics), to understand how the flow changes with different angular momentum. By changing the initial condition of the accretion torus minimally, we can simulate a steady, low-angular-momentum accretion flow around a Kerr black hole. We focus primarily on the lower limits of angular momentum and find that an accretion flow with an intermediate range of angular momentum differs significantly from high- or very-low-angular-momentum flows. The intermediate-angular-momentum accretion flow has the highest density, pressure, and temperature near the black hole, making it easier to observe. We find that the density and pressure have power-law scalings ρ ∝ r n−3/2 and p g ∝ r n−5/2, which only hold for very-low-angular-momentum cases. With the increase in flow angular momentum, it develops a nonaxisymmetric nature. In this case, simple self-similarity does not hold. We also find that the sonic surface moves away from the innermost stable circular orbit as the angular momentum decreases. Finally, we emphasize that an intermediate-angular-momentum flow could provide a possible solution to explaining the complex observation features of the supermassive black hole Sgr A* at our galactic center.

Magnetic Flux Plays an Important Role during a Black Hole X-Ray Binary Outburst in Radiative Two-temperature General Relativistic Magnetohydrodynamic Simulations

Black hole (Bh) X-ray binaries cycle through different spectral states of accretion over the course of months to years. Although persistent changes in the Bh mass accretion rate are generally recognized as the most important component of state transitions, it is becoming increasingly evident that magnetic fields play a similarly important role. In this article, we present the first radiative two-temperature general relativistic magnetohydrodynamics simulations in which an accretion disk transitions from a quiescent state at an accretion rate of Ṁ∼10−10ṀEdd to a hard-intermediate state at an accretion rate of Ṁ∼10−2ṀEdd . This huge parameter space in mass accretion rate is bridged by artificially rescaling the gas density scale of the simulations. We present two jetted BH models with varying degrees of magnetic flux saturation. We demonstrate that in “standard and normal evolution” models, which are unsaturated with magnetic flux, the hot torus collapses into a thin and cold accretion disk when Ṁ≳5×10−3ṀEdd . On the other hand, in “magnetically arrested disk” models, which are fully saturated with vertical magnetic flux, the plasma remains mostly hot with substructures that condense into cold clumps of gas when Ṁ≳1×10−2ṀEdd . This suggests that the spectral signatures observed during state transitions are closely tied to the level of magnetic flux saturation.

An operator-splitting numerical scheme for relativistic magnetohydrodynamics

We describe a novel operator-splitting approach to numerical relativistic magnetohydrodynamics designed to expand its applicability to the domain of ultra-high magnetisation. In this approach, the electromagnetic field is split into the force-free component, governed by the equations of force-free degenerate electrodynamics (FFDE), and the perturbation component, governed by the perturbation equations derived from the full system of relativistic magnetohydrodynamics (RMHD). The combined system of the FFDE and perturbation equations is integrated simultaneously, for which various numerical techniques developed for hyperbolic conservation laws can be used. At the end of every time-step of numerical integration, the force-free and the perturbation components of the electromagnetic field are recombined and the result is regarded as the initial value of the force-free component for the next time-step, whereas the initial value of the perturbation component is set to zero. To explore the potential of this approach, we build a 3rd-order WENO code, which was used to carry out 1D and 2D test simulations. Their results show that this operator-splitting approach allows us to bypass the stiffness of RMHD in the ultra-high-magnetisation regime where the perturbation component becomes very small. At the same time, the code performs very well in problems with moderate and low magnetisation too.

Preface

The Center for Space Plasma and Aeronomic Research (CSPAR) at the University of Alabama in Huntsville (UAH) and Maison de la Simulation, a joint laboratory at the French Alternative Energies and Atomic Energy Commission (CEA), National Center for Scientific Research (CNRS), University of Paris-Saclay, and University of Versailles-St Quentin, organized the 15th International Conference on Numerical Modeling of Space Plasma Flows (ASTRONUM-2023) on June 26—30, 2023 in Pasadena, California, USA.The Program Committee consisted of Nikolai Pogorelov (University of Alabama in Huntsville/CSPAR, USA, chair), Edouard Audit (CEA Maison de la Simulation, Gif-sur-Yvette, France, co-chair), Wes Bethel (San Francisco State University, USA), Amitava Bhattacharjee (Princeton Plasma Physics Laboratory, USA), Phillip Colella (Lawrence Berkeley National Laboratory, USA), Tomoyuki Hanawa (Chiba University, Japan), Maria Elena Innocenti (Ruhr University Bochum, Germany), Kanya Kusano (Nagoya University, Japan), Dongwook Lee (University of California, Santa Cruz), Jon Linker (Predictive Science Inc., USA), Anthony Mezzacappa (University of Tennessee, Knoxville, USA), Ewald Müller (Max-Planck-Institute for Astrophysics, Garching, Germany), Dongsu Ryu (Ulsan National Institute of Science and Technology, Korea), James Stone (Princeton University), and Gary P. Zank (University of Alabama in Huntsville, USA).The conference attracted over 90 scientists from 15 countries representing different branches of the plasma simulation community. The distinctive feature of this conference is a combination of diverse research topics, all of which are essential for performing high-resolution, continuum mechanics and particle, simulations of physical phenomena in space physics and astrophysics. Among such topics were software packages for modeling and analyzing plasma flows; advanced numerical methods for space and astrophysical flows; large-scale fluid-based, kinetic, and hybrid simulations; turbulence and cosmic ray transport; and magnetohydrodynamics. The discussed applications included cosmology and galaxy formation, supernova explosions, physics of the Sun-heliosphere-magnetosphere environments, space weather, the interstellar medium and star formation, stellar physics, experimental plasma physics, astrophysical accretion, numerical methods for ideal and non-ideal, relativistic and nonrelativistic MHD, etc. The proceedings volume is structured so that it covers all of these topics.The structure of the ASTRONUM conference series is based on the idea that modelers working in seemingly distant fields should have an opportunity to share their scientific achievements with the broad community of computational scientists performing numerical experiments. As in previous ASTRONUM meetings, we were interested in physical systems that are coupled across a multiplicity of spatial and temporal that incorporate diverse physical processes.The contributors to this volume are both young researchers and renowned experts in space physics and astrophysics, applied mathematics, and computational physics. This volume describes the application of numerical methods and the algorithms themselves, allowing us to discuss the challenges that theory imposes on numerical schemes for solving partial differential equations describing collisional and collisionless processes in space and astrophysical plasmas.We would like to thank the participants who submitted their papers to Proceedings of ASTRONUM-2023 and especially to those who reviewed manuscripts thus ensuring the high quality of this publication. We also are grateful to Adele Corona (ICNSMeetings.com) for the excellent management of the conference.This publication will be useful to graduate and postgraduate students majoring in space physics, astrophysics, numerical, engineering, and applied mathematics. It is also aimed at specialists in applied mathematics and various fields of physics that involve flows of partially ionized plasmas at both the collisional and collisionless levels.Nikolai V. Pogorelov, Edouard Audit, and Gary P. ZankMarch 3, 2024

Radiation Process in Relativistic MHD Waves: The Case of Circularly Polarized Alfvén Waves

Turbulence in highly magnetized plasma can be relativistic and induce an electric field comparable to the background magnetic field. Such a strong electric field can affect the emission process of nonthermal electrons. As the first step toward elucidating the emission process in relativistic turbulence, we study the radiation process of electrons in relativistic circularly polarized Alfvén waves. While the induced electric field boosts the average energy of low-energy electrons with a Larmor radius smaller than the wavelength, the emissivity for such electrons is suppressed because of the elongated gyromotion trajectory. The trajectory of high-energy electrons is shaken by the small-scale electric field, which enhances the emissivity. Since the effective Lorentz factor of E × B drift is ≃2 in the circularly polarized Alfvén waves, the deviation from the standard synchrotron emission is not so prominent. However, a power-law energy injection in the waves can produce a concave photon spectrum, which is similar to the GeV extra component seen in GRB spectra. If the turbulence electric field is responsible for the GeV extra component in GRBs, the estimates of the typical electron energy and magnetic field should be largely altered.

Magnetic Field Strength Effects on Nucleosynthesis from Neutron Star Merger Outflows

Magnetohydrodynamic turbulence drives the central engine of post-merger remnants, potentially powering both a nucleosynthetically active disk wind and the relativistic jet behind a short gamma-ray burst. We explore the impact of the magnetic field on this engine by simulating three post-merger black hole accretion disks using general relativistic magnetohydrodynamics with Monte Carlo neutrino transport, in each case varying the initial magnetic field strength. We find increasing ejecta masses associated with increasing magnetic field strength. We find that a fairly robust main r-process pattern is produced in all three cases, scaled by the ejected mass. Changing the initial magnetic field strength has a considerable effect on the geometry of the outflow and hints at complex central engine dynamics influencing lanthanide outflows. We find that actinide production is especially sensitive to magnetic field strength, with the overall actinide mass fraction calculated at 1 Gyr post-merger increasing by more than a factor of 6 with a tenfold increase in magnetic field strength. This hints at a possible connection to the variability in actinide enhancements exhibited by metal-poor, r-process-enhanced stars.

Imaging the event horizon of M87* from space on different timescales

Imaging the event horizon of M87* from space on different timescales

The 230 GHz Variability of Numerical Models of Sagittarius A*. I. Parameter Surveys on Varying the Ion-to-electron Temperature Ratio Under Strongly Magnetized Conditions

The 230 GHz lightcurves of Sagittarius A* (Sgr A*) predicted by general relativistic magnetohydrodynamics and general relativistic ray-tracing (GRRT) models by the Event Horizon Telescope Collaboration have higher variability M ΔT compared to observations. In this series of papers, we explore the origin of such large brightness variability. In this first paper, we performed large GRRT parameter surveys that span from the optically thin to the optically thick regimes, covering the ion-to-electron temperature ratio under strongly magnetized conditions, R Low, from 1 to 60. We find that increasing R Low can lead to either an increase or a reduction in M ΔT depending on the other model parameters, making it consistent with the observed variability of Sgr A* in some cases. Our analysis of GRRT image snapshots finds that the major contribution to the large M ΔT for the R Low = 1 models comes from the photon ring. However, secondary contributions from the accretion flow are also visible depending on the spin parameter. Our work demonstrates the importance of the electron temperature used for modeling radiatively inefficient accretion flows and places new constraints on the ion-to-electron temperature ratio. A more in-depth analysis for understanding the dependencies of M ΔT on R Low will be performed in subsequent papers.

Emergent Nucleosynthesis from a 1.2 s Long Simulation of a Black Hole Accretion Disk

We simulate a black hole accretion disk system with full-transport general relativistic neutrino radiation magnetohydrodynamics for 1.2 s. This system is likely to form after the merger of two compact objects and is thought to be a robust site of r-process nucleosynthesis. We consider the case of a black hole accretion disk arising from the merger of two neutron stars. Our simulation time coincides with the nucleosynthesis timescale of the r-process (∼1 s). Because these simulations are time-consuming, it is common practice to run for a “short” duration of approximately 0.1–0.3 s. We analyze the nucleosynthetic outflow from this system and compare the results of stopping at 0.12 and 1.2 s. We find that the addition of mass ejected in the longer simulation as well as more favorable thermodynamic conditions from emergent viscous ejecta greatly impacts the nucleosynthetic outcome. We quantify the error in nucleosynthetic outcomes between short and long cuts.

Relativistic shocks in conductive media

Context. Relativistic shocks are present in all high-energy astrophysical processes involving relativistic plasma outflows interacting with their ambient medium. While they are well understood in the context of relativistic hydrodynamics and ideal magnetohydrodynamics (MHD), there is a limited understanding of the properties related to their propagation in media characterized by finite electrical conductivity. Aims. This work presents a systematic method for the derivation and solution of the jump conditions for relativistic shocks propagating in MHD media with finite electrical conductivity. This method is applied to the numerical solution of the Riemann problem and the determination of the conditions inside the blastwave that is formed when ultrarelativistic magnetized ejecta interact with the circumburst medium during a gamma-ray burst. Methods. We derived the covariant relations expressing the jump conditions in a frame-independent manner. The resulting algebraic equations expressing the Rankine-Hugoniot conditions in the propagation medium’s frame were then solved numerically. A variable adiabatic index equation of state was used in order to obtain a realistic description of the post-shock fluid’s thermodynamics. This method was then employed for the solution of the Riemann problem for the case of a forward and a reverse shock, both of which form during the interaction of a gamma-ray burst ejecta with the circumburst medium. This allowed us to determine the kinematics of the resulting blastwave and the dynamical conditions in its interior. Results. Our solutions clearly depict the impact of the plasma’s electrical conductivity in the properties of the post-shock medium. Two characteristic regimes are identified with respect to the value of a dimensionless parameter that has a linear dependence on the conductivity. For small values of this parameter, the shock affects only the hydrodynamic properties of the propagation medium and leaves its electromagnetic field unaffected. No current layer forms in the shock front; thus, we refer to this as the current-free regime. For large values of this parameter, the ideal MHD regime has been retrieved. We also show that the assumption of a finite electrical conductivity can lead to higher efficiencies in the conversion of the ejecta energy into thermal energy of the blastwave through the reverse shock. The theory developed in this work can be applied to the construction of Riemann solvers for resistive relativistic MHD (RR4MHD).

PeV proton acceleration in gamma-ray binaries

Current generation of ground based gamma-ray telescopes observed dozens of sources of photons above 100 TeV. Supernova remnants, pulsar wind nebulae, young stellar clusters and superbubbles are considered as possible sites of PeV-regime particles producing the radiation. Another possible source of PeV particles could be gamma-ray binary systems. In these systems, a strong relativistic outflow from a compact object (neutron star or black hole) collides with the dense wind from a massive companion early-type star. Gamma-ray binaries are observed from radio to high energy gamma-rays as luminous non-thermal sources. Apart from acceleration of very high energy leptons producing most of the non-thermal radiation, these systems may also efficiently accelerate protons. We present here the results of numerical simulation of the PeV-regime proton acceleration in gamma-ray binaries. The simulation is based on relativistic MHD modeling of local flows of magnetized plasma in the region of interaction of two colliding winds. We then inject 0.1 PeV protons into the system and directly follow their trajectories to demonstrate that they are accelerated to energies above PeV. High magnetization of the wind of the young massive star providing a Gauss range field in the winds interaction region is of paramount importance for the acceleration of protons above PeV. The maximum energies of protons accelerated by colliding winds in gamma-ray binaries can significantly exceed the energy of the pulsar potential’s drop, which limits from above the energy of particles accelerated by an isolated pulsar.

A 4th-order accurate finite volume method for ideal classical and special relativistic MHD based on pointwise reconstructions

We present a novel implementation of a genuinely 4th-order accurate finite volume scheme for multidimensional classical and special relativistic magnetohydrodynamics (MHD) based on the constrained transport (CT) formalism. The scheme introduces several novel aspects when compared to its predecessors yielding a more efficient computational tool. Among the most relevant ones, our scheme exploits pointwise to pointwise reconstructions (rather than one-dimensional finite volume ones), employs the generic upwind constrained transport averaging and sophisticated limiting strategies that include both a discontinuity detector and an order reduction procedure. Selected numerical benchmarks demonstrate the accuracy and robustness of the method.

Charged participants and their electromagnetic fields in an expanding fluid

We investigate the space-time dependence of electromagnetic fields produced by charged participants in an expanding fluid. To address this problem, we need to solve the Maxwell's equations coupled to the hydrodynamics conservation equation, specifically the relativistic magnetohydrodynamics (RMHD) equations, since the charged participants move with the flow. To gain analytical insight, we approximate the problem by solving the equations in a fixed background Bjorken flow, onto which we solve Maxwell's equations. The dynamical electromagnetic fields interact with the fluid's kinematic quantities such as the shear tensor and the expansion scalar, leading to additional non-trivial coupling. We use mode decomposition of Green's function to solve the resulting non-linear coupled wave equations. We then use this function to calculate the electromagnetic field for two test cases: a point source and a transverse charge distribution. The results show that the resulting magnetic field vanishes at very early times, grows, and eventually falls at later times.

A positivity-preserving adaptive-order finite-difference scheme for GRMHD

We present an adaptive-order positivity-preserving conservative finite-difference scheme that allows a high-order solution away from shocks and discontinuities while guaranteeing positivity and robustness at discontinuities. This is achieved by monitoring the relative power in the highest mode of the reconstructed polynomial and reducing the order when the polynomial series no longer converges. Our approach is similar to the multidimensional optimal order detection strategy, but differs in several ways. The approach is a priori and so does not require retaking a time step. It can also readily be combined with positivity-preserving flux limiters that have gained significant traction in computational astrophysics and numerical relativity. This combination ultimately guarantees a physical solution both during reconstruction and time stepping. We demonstrate the capabilities of the method using a standard suite of very challenging 1d, 2d, and 3d general relativistic magnetohydrodynamics test problems.

Ameliorating the Courant-Friedrichs-Lewy condition in spherical coordinates: A double FFT filter method for general relativistic MHD in dynamical spacetimes

Ameliorating the Courant-Friedrichs-Lewy condition in spherical coordinates: A double FFT filter method for general relativistic MHD in dynamical spacetimes

On the Escape of Low-frequency Waves from Magnetospheres of Neutron Stars

We study the nonlinear decay of the fast magnetosonic (fms) into the Alfvén waves in relativistic force-free magnetohydrodynamics. The work has been motivated by models of pulsar radio emission and fast radio bursts (FRBs), in which the emission is generated in neutron star magnetospheres at conditions when not only the Larmor but also the plasma frequencies significantly exceed the radiation frequency. The decay process places limits on the source luminosity in these models. We estimated the decay rate and showed that the phase volume of Alfvén waves available for the decay of an fms wave is infinite. Therefore, the energy of fms waves could be completely transferred to the small-scale Alfvén waves not via a cascade, as in the Kolmogorov turbulence, but directly. Our results explain the anomalously low radio efficiency of the Crab pulsar and show that FRBs could not be produced well within magnetar magnetospheres.

Resistive relativistic MHD simulations of astrophysical jets

Aims. The main goal of the present paper is to provide the first systematic numerical study of the propagation of astrophysical relativistic jets, in the context of high-resolution, shock-capturing Resistive Relativistic MagnetoHydroDynamic (RRMHD) simulations. We aim to investigate different values and models for the plasma resistivity coefficient, and to assess their impact on the level of turbulence, the formation of current sheets and reconnection plasmoids, the electromagnetic energy content, and the dissipated power. Methods. We used the PLUTO code for simulations and we assumed an axisymmetric setup for the jets, endowed with both poloidal and toroidal magnetic fields, and propagating in a uniform magnetized medium. The gas was assumed to be characterized by a realistic, Synge-like equation of state (the Taub equation), appropriate for such astrophysical jets. The Taub equation was combined here for the first time with the implicit-explicit Runge-Kutta time-stepping procedure, as required in RRMHD simulations. Results. The main result is that turbulence is clearly suppressed for the highest values of resistivity (low Lundquist numbers), current sheets are broader, and plasmoids are barely present, while for low values of resistivity the results are very similar to ideal runs, in which dissipation is purely numerical. We find that recipes employing a variable resistivity based on the advection of a jet tracer or on the assumption of a uniform Lundquist number improve on the use of a constant coefficient and are probably more realistic possible sites for the acceleration of the nonthermal particles that produce the observed high-energy emission, preserving as they do the development of turbulence and of sharp current sheets.

κmonty: a Monte Carlo Compton scattering code including non-thermal electrons

ABSTRACT Low-luminosity active galactic nuclei are strong sources of X-ray emission produced by Compton scattering originating from the accretion flows surrounding their supermassive black holes. The shape and energy of the resulting spectrum depend on the shape of the underlying electron distribution function (DF). In this work, we present an extended version of the grmonty code, called κmonty. The grmonty code previously only included a thermal Maxwell–Jütner electron DF. We extend the grmonty code with non-thermal electron DFs, namely the κ and power-law DFs, implement Cartesian Kerr–Schild coordinates, accelerate the code with mpi, and couple the code to the non-uniform adaptive mesh refinement grid data from the general relativistic magnetohydrodynamics code bhac. For the Compton scattering process, we derive two sampling kernels for both DFs. Finally, we present a series of code tests to verify the accuracy of our schemes. The implementation of non-thermal DFs opens the possibility of studying the effect of non-thermal emission on previously developed black hole accretion models.