Source author record

O. Porth

O. Porth appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.HE astro-ph.GA astro-ph.IM physics.plasm-ph

Catalog footprint

What is connected

10works

4topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2026arXiv

Magnetic reconnection with a 0.1 rate: Effective resistivity in general relativistic magnetohydrodynamics

Relativistic magnetic reconnection is thought to power various multi-wavelength emission signatures from neutron stars and black holes. Relativistic resistive magnetohydrodynamics (RRMHD) offers the simplest model of reconnection. However, a small uniform resistivity underestimates the reconnection rate compared to first-principles kinetic models. By employing an effective resistivity based on kinetic models - which connects the reconnection electric field to the charge-starved current density - we show that RRMHD can reproduce the increased reconnection rate of kinetic models, both in local current sheets and in global black hole magnetospheres.

preprint2025arXiv

On the universality of the split monopole black hole magnetosphere

Black holes can acquire magnetic flux from their magnetized progenitor or via prolonged accretion. We study the evolution of black hole magnetospheres by means of axisymmetric general relativistic magnetohydrodynamic simulations. We show that all simulated initial magnetic field geometries of varying complexity ultimately evolve into a split monopole magnetosphere. The magnetospheric evolution consists of two phases. In the first phase, the magnetosphere evolves toward pressure equilibrium accompanied by a large magnetic flux decrease on the event horizon on a fast Alfvénic timescale of $\sim 60$ light-crossing times of the gravitational radius. The second phase proceeds in a pressure balance in which the magnetic flux decays and current sheets shift in polar angle over the event horizon on slower resistive timescales. We present an analytic model for the second phase. Furthermore, we show that in a split monopole magnetosphere the magnetic flux on the event horizon decays exponentially with a timescale that depends on the black hole spin, where higher spin results in slower decay. Our results can have an implications for the timescales of reconnection-powered flares and for multimessenger counterparts to gravitational wave events.

preprint2022arXiv

Formation of Magnetically Truncated Accretion Disks in 3D Radiation-Transport Two-Temperature GRMHD Simulations

Multi-wavelength observations suggest that the accretion disk in the hard and intermediate states of X-ray binaries (XRBs) and active galactic nuclei (AGN) transitions from a cold, thin disk at large distances into a hot, thick flow close to the black hole. However, the formation, structure and dynamics of such truncated disks are poorly constrained due to the complexity of the thermodynamic, magnetic, and radiative processes involved. We present the first radiation-transport two-temperature general relativistic magnetohydrodynamic (GRMHD) simulations of truncated disks radiating at ~35% of the Eddington luminosity with and without large-scale poloidal magnetic flux. We demonstrate that when a geometrically-thin accretion disk is threaded by large-scale net poloidal magnetic flux, it self-consistently transitions at small radii into a two-phase medium of cold gas clumps floating through a hot, magnetically dominated corona. This transition occurs at a well-defined truncation radius determined by the distance out to which the disk is saturated with magnetic flux. The average ion and electron temperatures in the semi-opaque corona reach, respectively, T_i ~ 10^10K and T_e ~ 5 10^8K. The system produces radiation, powerful collimated jets and broader winds at the total energy efficiency exceeding ~90%, the highest ever energy extraction efficiency from a spinning black hole by a radiatively efficient flow in a GRMHD simulation. This is consistent with jetted ejections observed during XRB outbursts. The two-phase medium may naturally lead to broadened iron line emission observed in the hard state.

preprint2021arXiv

Flares in the Galactic center I: orbiting flux tubes in Magnetically Arrested Black Hole Accretion Disks

Recent observations of SgrA* by the GRAVITY instrument have astrometrically tracked infrared flares (IR) at distances of $\sim 10$ gravitational radii ($r_g$). In this paper, we study a model for the flares based on 3D general relativistic magnetohydrodynamic (GRMHD) simulations of magnetically arrested accretion disks (MADs) which exhibit violent episodes of flux escape from the black hole magnetosphere. These events are attractive for flare modeling for several reasons: i) the magnetically dominant regions can resist being disrupted via magneto-rotational turbulence and shear, ii) the orientation of the magnetic field is predominantly vertical as suggested by the GRAVITY data, iii) magnetic reconnection associated with the flux eruptions could yield a self-consistent means of particle heating/acceleration during the flare events. In this analysis we track erupted flux bundles and provide distributions of sizes, energies and plasma parameter. In our simulations, the orbits tend to circularize at a range of radii from $\sim 5-40 r_g$. The magnetic energy contained within the flux bundles ranges up to $\sim10^{40}$ erg, enough to power IR and X-ray flares. We find that the motion within the magnetically supported flow is substantially sub-Keplerian, in tension with the inferred period-radius relation of the three GRAVITY flares.

preprint2020arXiv

Relativistic AGN jets III. Synthesis of synchrotron emission from Double-Double Radio Galaxies

The class of Double-Double Radio Galaxies (DDRGs) relates to episodic jet outbursts. How various regions and components add to the total intensity in radio images is less well known. In this paper we synthesize synchrotron images for DDRGs based on special relativistic hydrodynamic simulations, making advanced approximations for the magnetic fields. We study the synchrotron images for: Three different radial jet profiles; Ordered, entangled or mixed magnetic fields; Spectral ageing from synchrotron cooling; The contribution from different jet components; The viewing angle and Doppler (de-)boosting; The various epochs of the evolution of the DDRG. To link our results to observational data, we adopt to J1835+6204 as a reference source. In all cases the synthesized synchrotron images show two clear pairs of hotspots, in the inner and outer lobes. The best resemblance is obtained for the piecewise isochoric jet model, for a viewing angle of approximately $\vartheta \sim -71^{\circ}$, i.e. inclined with the lower jet towards the observer, with predominantly entangled ($\gtrsim 70$ per cent of the magnetic pressure) in turbulent, rather than ordered fields. The effects of spectral ageing become significant when the ratio of observation frequencies and cut-off frequency $ν_{\rm obs}/ν_{\infty,0} \gtrsim 10^{-3}$, corresponding to $\sim 3 \cdot 10^2$ MHz. For viewing angles $\vartheta \lesssim -30^{\circ}$, a DDRG morphology can no longer be recognized. The second jets must be injected within $\lesssim$ 4 per cent of the lifetime of the first jets for a DDRG structure to emerge, which is relevant for Active Galactic Nuclei feedback constraints.

preprint2016arXiv

Diffusion in Pulsar Wind Nebulae: An Investigation using Magnetohydrodynamic and Particle Transport Models

We study the transport of high-energy particles in pulsar wind nebulae (PWN) using three-dimensional MHD (see Porth et al. (2014) for details) and test-particle simulations, as well as a Fokker-Planck particle transport model. The latter includes radiative and adiabatic losses, diffusion, and advection on the background flow of the simulated MHD nebula. By combining the models, the spatial evolution of flux and photon index of the X-ray synchrotron emission is modelled for the three nebulae G21.5-0.9, the inner regions of Vela, and 3C 58, thereby allowing us to derive governing parameters: the magnetic field strength, average flow velocity and spatial diffusion coefficient. For comparison, the nebulae are also modelled with the semi-analytic Kennel & Coroniti (1984) model but the Porth et al. (2014) model generally yields better fits to the observational data. We find that high velocity fluctuations in the turbulent nebula (downstream of the termination shock) give rise to efficient diffusive transport of particles, with average Péclet number close to unity, indicating that both advection and diffusion play an important role in particle transport. We find that the diffusive transport coefficient of the order of $\sim2\times 10^{27} (L_{\rm s}/0.42\rm Ly) cm^{2}s^{-1}$ ($L_{\rm s}$ is the size of the termination shock) is independent of energy up to extreme particle Lorentz factors of $γ_{p}\sim10^{10}$.

preprint2016arXiv

Stationary Relativistic Jets

In this paper we describe a simple numerical approach which allows to study the structure of steady-state axisymmetric relativistic jets using one-dimensional time-dependent simulations. It is based on the fact that for narrow jets with v~c the steady-state equations of relativistic magnetohydrodynamics can be accurately approximated by the one-dimensional time-dependent equations after the substitution z=ct. Since only the time-dependent codes are now publicly available this is a valuable and efficient alternative to the development of a high-specialized code for the time-independent equations. The approach is also much cheaper and more robust compared to the relaxation method. We tested this technique against numerical and analytical solutions found in literature as well as solutions we obtained using the relaxation method and found it sufficiently accurate. In the process, we discovered the reason for the failure of the self-similar analytical model of the jet reconfinement in relatively flat atmospheres and elucidated the nature of radial oscillations of steady-state jets.

preprint2015arXiv

Causality and stability of cosmic jets

In stark contrast to their laboratory and terrestrial counterparts, the cosmic jets appear to be very stable. We propose that the reason behind this remarkable property is the loss of causal connectivity across these jets, caused by their rapid expansion in response to fast decline of external pressure with the distance from the "jet engine". In atmospheres with power-law pressure distribution, the total loss of causal connectivity occurs, when the power index k>2 - the steepness which is expected to be quite common for many astrophysical environments. This conclusion does not seem to depend on the physical nature of jets - it applies both to relativistic and non-relativistic flows, both magnetically-dominated and unmagnetized jets. In order to verify it, we have carried out numerical simulations of moderately magnetized and moderately relativistic jets. Their results give strong support to our hypothesis and provide with valuable insights. In particular, we find that the z-pinched inner cores of magnetic jets expand slower than their envelopes and become susceptible to instabilities even when the whole jet is stable. This may result in local dissipation and emission without global disintegration of the flow. Cosmic jets may become globally unstable when they enter flat sections of external atmospheres. We propose that the Fanaroff-Riley morphological division of extragalactic radio sources into two classes is related to this issue. In particular, we argue that the low power FR-I jets become re-confined, causally connected and globally unstable on the scale of galactic X-ray coronas, whereas more powerful FR-II jets re-confine much further out and remain largely intact until they terminate at hot spots. Using this idea, we derived the relationship between the critical jet power and the optical luminosity of the host galaxy, which is in a very good agreement with the observations.

preprint2014arXiv

MPI-AMRVAC for Solar and Astrophysics

In this paper we present an update on the open source MPI-AMRVAC simulation toolkit where we focus on solar- and non-relativistic astrophysical magneto-fluid dynamics. We highlight recent developments in terms of physics modules such as hydrodynamics with dust coupling and the conservative implementation of Hall magnetohydrodynamics. A simple conservative high-order finite difference scheme that works in combination with all available physics modules is introduced and demonstrated at the example of monotonicity preserving fifth order reconstruction. Strong stability preserving high order Runge-Kutta time steppers are used to obtain stable evolutions in multidimensional applications realizing up to fourth order accuracy in space and time. With the new distinction between active and passive grid cells, MPI-AMRVAC is ideally suited to simulate evolutions where parts of the solution are controlled analytically, or have a tendency to progress into or out of a stationary state. Typical test problems and representative applications are discussed, with an outlook to follow-up research. Finally, we discuss the parallel scaling of the code and demonstrate excellent weak scaling up to 30 000 processors allowing to exploit modern petascale infrastructure.

preprint2009arXiv

Acceleration and collimation of relativistic MHD disk winds

We perform axisymmetric relativistic magnetohydrodynamic (MHD) simulations to investigate the acceleration and collimation of jets and outflows from disks around compact objects. The fiducial disk surface (respectively a slow disk wind) is prescribed as boundary condition for the outflow. We apply this technique for the first time in the context of relativistic jets. The strength of this approach is that it allows us to run a parameter study in order to investigate how the accretion disk conditions govern the outflow formation. Our simulations using the PLUTO code run for 500 inner disk rotations and on a physical grid size of 100x200 inner disk radii. In general, we obtain collimated beams of mildly relativistic speed and mass-weighted half-opening angles of 3-7 degrees. When we increase the outflow Poynting flux by injecting an additional disk toroidal field into the inlet, Lorentz factors up to 6 are reached. These flows gain super-magnetosonic speed and remain Poynting flux dominated. The light surface of the outflow magnetosphere tends to align vertically - implying three relativistically distinct regimes in the flow - an inner sub-relativistic domain close to the jet axis, a (rather narrow) relativistic jet and a surrounding subrelativistic outflow launched from the outer disk surface - similar to the spine-sheath structure currently discussed for asymptotic jet propagation and stability. The outer subrelativistic disk wind is a promising candidate for the X-ray absorption winds that are observed in many radio-quiet AGN.