Relativistic Magnetohydrodynamics Research Articles

On a degenerate boundary value problem to relativistic magnetohydrodynamics with a general pressure law

On a degenerate boundary value problem to relativistic magnetohydrodynamics with a general pressure law

Binary neutron star mergers using a discontinuous Galerkin-finite difference hybrid method

Abstract We present a discontinuous Galerkin-finite difference hybrid scheme that allows high-order shock capturing with the discontinuous Galerkin method for general relativistic magnetohydrodynamics in dynamical spacetimes. We present several optimizations and stability improvements to our algorithm that allow the hybrid method to successfully simulate single, rotating, and binary neutron stars. The hybrid method achieves the efficiency of discontinuous Galerkin methods throughout almost the entire spacetime during the inspiral phase, while being able to robustly capture shocks and resolve the stellar surfaces. We also use Cauchy-Characteristic evolution to compute the first gravitational waveforms at future null infinity from binary neutron star mergers. The simulations presented here are the first successful binary neutron star inspiral and merger simulations using discontinuous Galerkin methods.

Bridging simulations of kink instability in relativistic magnetized jets with radio emission and polarisation

Relativistic outflows emanating from active galactic nuclei can extend up to kiloparsec scales in length, displaying a variety of complex morphologies. This study explores the intricate morphologies of such relativistic jets, mainly focusing on creating a bridge between magnetic instabilities in jets with observational signatures from complex radio galaxies. In particular, we aim to study the role of dynamical instabilities in forming distinctive morphological features by employing 3D relativistic magnetohydrodynamic (RMHD) simulations of rotating jets. Our simulations have further used the hybrid Eulerian-Lagrangian framework of the PLUTO code and generated the synthetic synchrotron emission and polarisation maps to compare with the observed signatures. Our analysis based on simulations of a continuously injected jet suggests that current-driven instabilities, notably the |m|=1 mode, generate ribs-like structures that are seen in some of the recent radio galaxies using MeerKat, e.g. MysTail. In our contrasting simulations of the restarted jet, the kink-instability driven ribs-like structures were formed relatively near the nozzle. In both cases, the jet dissipates its pre-existing magnetic energy through these instabilities, transitioning to a more kinetic energy dominant state. The turbulent structures resulting from this dissipation phase are filamentary and resemble the tethers as observed for the case of MysTail. This pilot study essentially provides a plausible qualitative explanation by bridging simulations of kink instability to produce synthetic radio features resembling the observed complex radio morphology of MysTail.

3D hybrid fluid-particle jet simulations and the importance of synchrotron radiative losses

Relativistic jets in active galactic nuclei are known for their exceptional energy output, and imaging the synthetic synchrotron emission of numerical jet simulations is essential for a comparison with observed jet polarization emission. Through the use of 3D hybrid fluid-particle jet simulations (with the PLUTO code), we overcome some of the commonly made assumptions in relativistic magnetohydrodynamic (RMHD) simulations by using non-thermal particle attributes to account for the resulting synchrotron radiation. Polarized radiative transfer and ray-tracing (via the RADMC-3D code) highlight the differences in total intensity maps when (i) the jet is simulated purely with the RMHD approach, (ii) a jet tracer is considered in the RMHD approach, and (iii) a hybrid fluid-particle approach is used. The resulting emission maps were compared to the example of the radio galaxy Centaurus A. We applied the Lagrangian particle module implemented in the latest version of the PLUTO code. This new module contains a state-of-the-art algorithm for modeling diffusive shock acceleration and for accounting for radiative losses in RMHD jet simulations. The module implements the physical postulates missing in RMHD jet simulations by accounting for a cooled ambient medium and strengthening the central jet emission. We find a distinction between the innermost structure of the jet and the back-flowing material by mimicking the radio emission of the Seyfert II radio galaxy Centaurus\,A when considering an edge-brightened jet with an underlying purely toroidal magnetic field. We demonstrate the necessity of synchrotron cooling as well as the improvements gained when directly accounting for non-thermal synchrotron radiation via non-thermal particles.

General Relativistic Radiation Magnetohydrodynamics Simulations of Precessing Tilted Super-Eddington Disks

We perform a three-dimensional general relativistic radiation magnetohydrodynamics simulation of a tilted super-Eddington accretion disk around a spinning black hole (BH). The disk, which tilts and twists as it approaches the BH, precesses while maintaining its shape. The gas is mainly ejected around the rotation axis of the outer part of the disk rather than around the spin axis of the BH. The disk precession changes the ejection direction of the gas with time. The radiation energy is also released in approximately the same direction as the outflow, so the precession is expected to cause a quasiperiodic time variation of the observed luminosity. The timescale of the precession is about 10 s for a 10 M ⊙ BH and for the radial extent of a disk of several tens of gravitational radii, where M ⊙ is the solar mass. This timescale is consistent with the frequency of the low-frequency quasiperiodic oscillation (0.01–1 Hz) observed in some ultraluminous X-ray sources.

Images of Kerr-MOG black holes surrounded by geometrically thick magnetized equilibrium tori

We adopt general relativistic ray-tracing (GRRT) schemes to study images of Kerr-MOG black holes surrounded by geometrically thick magnetized equilibrium tori, which belong to steady-state solutions of thick accretion disks within the framework of general relativistic magnetohydrodynamics (GRMHD). The black hole possesses an extra dimensionless MOG parameter described its deviation from usual Kerr one. Our results show that the presence of the MOG parameter leads to smaller disks in size, but enhances the total flux density and peak brightness in their images. Combining with observation data of black hole M87* from the Event Horizon Telescope (EHT), we make a constraint on parameters of the Kerr-MOG black hole and find that the presence of the MOG parameter broadens the allowable range of black hole spin.

GRMHD Simulations of Accretion Structures with Different Angular Momentum Profiles

In this article, we explore the dynamics of accretion structures encircling spherically symmetric black holes, comparing three accretion disk models with distinct angular momentum profiles: (i) the geometrically thin Keplerian disk, (ii) the Fishbone–Moncrief torus; and (iii) the Polish Doughnut. Employing general relativistic magnetohydrodynamics simulations with the High Accuracy Relativistic Magnetohydrodynamics code, we investigate these three models, considering the magnetic field’s influence on the accretion disk angular momentum redistribution. We show that the magnetic field is a key factor in accretion disk structures, especially in regions with lower mass density. Our investigation verifies the well-established fact that the presence of a magnetic field significantly influences the accretion rate and its temporal variability.

A Revisited Equilibrium Solution of the Fishbone and Moncrief Torus for Extended General Relativistic Magnetohydrodynamic Simulations

Accretion physics has become more important recently due to the detection of the first horizon-scale images of the supermassive black holes of M 87* and Sgr A* by the Event Horizon Telescope. General relativistic magnetohydrodynamic (GRMHD) simulations of magnetized accretion flows onto a Kerr black hole have been used to interpret them. However, further testing the theory of gravity by using horizon-scale images requires performing consistent GRMHD simulations in non-Kerr spacetime. In this paper, we revisited the hydrodynamical equilibrium solution of the Fishbone and Moncrief (FM) torus that can be used to study any stationary, axisymmetric, vacuum, or nonvacuum spacetime. Further, we check the stability of the FM torus in non-Kerr spacetime by general relativistic hydrodynamic simulations. We find that FM torus in non-Kerr spacetime is indeed stable under long-term evolution. We conclude that the generalized FM torus solution would be very useful for creating new GRMHD libraries in extended Kerr black holes.

A fourth-order accurate finite volume scheme for resistive relativistic MHD

A fourth-order accurate finite volume scheme for resistive relativistic MHD

Parametric decay instability of circularly polarized Alfvén waves in magnetically dominated plasma.

We investigate parametric decay instability (PDI) of circularly polarized Alfvén wave into daughter acoustic wave and backward Alfvén wave in magnetically dominated plasma, in which the magnetization parameter σ (energy density ratio of background magnetic field to matter) exceeds unity. We analyze relativistic magnetohydrodynamics (MHD), focusing on wave frequencies sufficiently lower than the plasma and cyclotron frequencies. We derive analytical formulas for the dispersion relation and growth rate of the instability as a function of the magnetization σ, wave amplitude η, and plasma temperature θ. We find that PDI persists even in high magnetization σ, albeit with a decreased growth rate up to σ→∞. Our formulas are useful for estimating the decay of Alfvén wave into acoustic wave and heat in high magnetization σ plasma, which is a ubiquitous phenomenon such as in pulsars, magnetars, and fast radio bursts.

Revisiting shear stress tensor evolution: Nonresistive magnetohydrodynamics with momentum-dependent relaxation time

This study aims to develop second-order relativistic viscous magnetohydrodynamics (MHD) derived from kinetic theory within an extended relaxation time approximation (momentum/energy dependent) for the collision kernel. The investigation involves a detailed examination of shear stress tensor evolution equations and associated transport coefficients. The Boltzmann equation is solved using a Chapman-Enskog-like gradient expansion for a charge-conserved conformal system, incorporating a momentum-dependent relaxation time. The derived relativistic nonresistive, viscous second-order MHD equations for the shear stress tensor reveal significant modifications in the coupling with dissipative charge current and magnetic field due to the momentum dependence of the relaxation time. By utilizing a power law parametrization to quantify the momentum dependence of the relaxation time, the anisotropic magnetic field-dependent shear coefficients in the Navier-Stokes limit have been investigated. The resulting viscous coefficients are seen to be sensitive to the momentum dependence of the relaxation time. Published by the American Physical Society 2024

Reconstruction of Fermi and eROSITA Bubbles from Magnetized Jet Eruption with Simulations

The Fermi bubbles and the eROSITA bubbles around the Milky Way Galaxy are speculated to be the aftermaths of past jet eruptions from a supermassive black hole in the galactic center. In this work, a 2.5D axisymmetric relativistic magnetohydrodynamic (RMHD) model is applied to simulate a jet eruption from our galactic center and to reconstruct the observed Fermi bubbles and eROSITA bubbles. High-energy non-thermal electrons are excited around forward shock and discontinuity transition regions in the simulated plasma distributions. The γ-ray and X-ray emissions from these electrons manifest patterns on the skymap that match the observed Fermi bubbles and eROSITA bubbles, respectively, in shape, size and radiation intensity. The influence of the background magnetic field, initial mass distribution in the Galaxy, and the jet parameters on the plasma distributions and hence these bubbles is analyzed. Subtle effects on the evolution of plasma distributions attributed to the adoption of a galactic disk model versus a spiral-arm model are also studied.

GQL-based bound-preserving and locally divergence-free central discontinuous Galerkin schemes for relativistic magnetohydrodynamics

This paper develops novel and robust central discontinuous Galerkin (CDG) schemes of arbitrarily high-order accuracy for special relativistic magnetohydrodynamics (RMHD) with a general equation of state (EOS). These schemes are provably bound-preserving (BP), meaning they consistently preserve the upper bound for subluminal fluid velocity and the positivity of density and pressure, while also (locally) maintaining the divergence-free (DF) constraint for the magnetic field. For 1D RMHD, the standard CDG method is exactly DF, and its BP property is proven under a condition achievable by the BP limiter. For 2D RMHD, we design provably BP and locally DF CDG schemes based on the suitable discretization of a modified RMHD system, which is the relativistic analogue of Godunov's symmetrizable form of the non-relativistic MHD system (Godunov, 1972 [19]). A key novelty in our schemes is the meticulous discretization of additional source terms in the modified RMHD equations, so as to precisely counteract the influence of divergence errors on the BP property across overlapping meshes. Notably, we provide rigorous proofs of the BP property for our CDG schemes and first establish the theoretical connection between BP and discrete DF properties on overlapping meshes for RMHD. Owing to the absence of explicit expressions for primitive variables in terms of conserved variables, the constraints of physical bounds are strongly nonlinear, making the BP proofs highly nontrivial. We overcome these challenges through technical estimates within the geometric quasilinearization (GQL) framework (Wu and Shu, 2023 [49]), which equivalently converts the nonlinear constraints into linear ones. Furthermore, we introduce a new 2D cell average decomposition on overlapping meshes, which relaxes the theoretical BP CFL constraint and reduces the number of internal nodes, thereby enhancing the efficiency of the 2D BP CDG method. Finally, we implement the proposed CDG schemes for extensive RMHD problems with various EOSs, demonstrating their robustness and effectiveness in challenging scenarios like ultra-relativistic blasts and jets in strongly magnetized environments.

A unified accretion disc model for supermassive black holes in galaxy formation simulations: method and implementation

ABSTRACT It is well established that supermassive black hole (SMBH) feedback is crucial for regulating the evolution of massive, if not all, galaxies. However, modelling the interplay between SMBHs and their host galaxies is challenging due to the vast dynamic range. Previous simulations have utilized simple subgrid models for SMBH accretion, while recent advancements track the properties of the unresolved accretion disc, usually based on the thin α-disc model. However, this neglects accretion in the radiatively inefficient regime, expected to occur through a thick disc for a significant portion of an SMBH’s lifetime. To address this, we present a novel ‘unified’ accretion disc model for SMBHs, harnessing results from the analytical advection-dominated inflow–outflow solution (ADIOS) model and state-of-the-art general relativistic (radiation-)magnetohydrodynamics (GR(R)MHD) simulations. Going from low to high Eddington ratios, our model transitions from an ADIOS flow to a thin α-disc via a truncated disc, incorporating self-consistently SMBH spin evolution due to Lense–Thirring precession. Utilizing the moving mesh code arepo, we perform simulations of single and binary SMBHs within gaseous discs to validate our model and assess its impact. The disc state significantly affects observable luminosities, and we predict markedly different electromagnetic counterparts in SMBH binaries. Crucially, the assumed disc model shapes SMBH spin magnitudes and orientations, parameters that gravitational wave observatories like LISA and IPTA are poised to constrain. Our simulations emphasize the importance of accurately modelling SMBH accretion discs and spin evolution, as they modulate the available accretion power, profoundly shaping the interaction between SMBHs and their host galaxies.

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

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

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 3C 84 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 3C 84 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 ∼0.1 parsec. We tested different simulation set-ups by altering the bulk Lorentz factor Γ 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.

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

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

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.