Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

Hot spots, or plasmoids, forming due to magnetic reconnection in thin current sheets, are conjectured to power frequent bright X-ray and near-infrared flares from Sgr A*, the supermassive black hole in the center of our Galaxy. It is of yet unclear how, where, and when such thin current sheets form in black-hole accretion disks or astrophysical jets. In this work we perform axisymmetric general-relativistic resistive magnetohydrodynamics simulations to model magnetic reconnection and associated plasmoid formation in a wide range of accretion flows. We show that current sheets and plasmoids are ubiquitous features which form regardless of the initial size of the disk and the magnetization in the quasi-steady-state phase of accretion. In all cases we observe plasmoids forming in current sheets close to the event horizon within 5 to 10 Schwarzschild radii. These plasmoids can then merge, grow to macroscopic scales of the order of a Schwarzschild radius, and are ultimately advected along the jet's sheath or into the disk. The largest plasmoids are efficiently energized to relativistic temperatures via magnetic reconnection and contribute to the heating of the jet's sheath. In all cases we find reconnection rates between 0.01c and 0.03c, consistent with studies of reconnection in isolated Harris-type current sheets. We quantify magnetic dissipation and strong non-ideal electric fields which can efficiently inject non-thermal particles. We also show that an explicit resistivity allows for converged numerical solutions, such that the electromagnetic energy density evolution and dissipation become independent of the grid scale for the extreme resolutions considered here.

We review and discuss insights on ideal magnetohydrodynamic (MHD) instabilities that can play a role in destabilizing solar coronal flux rope structures. For single flux ropes, failed or actual eruptions may result from internal or external kink evolutions, or from torus unstable configurations. We highlight recent findings from 3D magnetic field reconstructions and simulations where kink and torus instabilities play a prominent role. For interacting current systems, we critically discuss different routes to coronal dynamics and global eruptions, due to current channel coalescence or to tilt-kink scenarios. These scenarios involve the presence of two nearby current channels and are clearly distinct from the popular kink or torus instability. Since the solar corona is pervaded with myriads of magnetic loops -- creating interacting flux ropes typified by parallel or antiparallel current channels as exemplified in various recent observational studies -- coalescence or tilt-kink evolutions must be very common for destabilizing adjacent flux rope systems. Since they also evolve on ideal MHD timescales, they may well drive many sympathetic eruptions witnessed in the solar corona. Moreover, they necessarily lead to thin current sheets that are liable to reconnection. We review findings from 2D and 3D MHD simulations for tilt and coalescence evolutions, as well as on particle acceleration aspects derived from computed charged particle motions in tilt-kink disruptions and coalescing flux ropes. The latter were recently studied in two-way coupled kinetic-fluid models, where the complete phase-space information of reconnection is incorporated.

Recent advances in black hole astrophysics, particularly the first visual evidence of a supermassive black hole at the center of the galaxy M87 by the Event Horizon Telescope (EHT), and the detection of an orbiting "hot spot" nearby the event horizon of Sgr A* in the Galactic center by the Gravity Collaboration, require the development of novel numerical methods to understand the underlying plasma microphysics. Non-thermal emission related to such hot spots is conjectured to originate from plasmoids that form due to magnetic reconnection in thin current layers in the innermost accretion zone. Resistivity plays a crucial role in current sheet formation, magnetic reconnection, and plasmoid growth in black hole accretion disks and jets. We included resistivity in the three-dimensional general-relativistic magnetohydrodynamics (GRMHD) code BHAC and present the implementation of an Implicit-Explicit scheme to treat the stiff resistive source terms of the GRMHD equations. The algorithm is tested in combination with adaptive mesh refinement to resolve the resistive scales and a constrained transport method to keep the magnetic field solenoidal. Several novel methods for primitive variable recovery, a key part in relativistic magnetohydrodynamics codes, are presented and compared for accuracy, robustness, and efficiency. We propose a new inversion strategy that allows for resistive-GRMHD simulations of low gas-to-magnetic pressure ratio and highly magnetized regimes as applicable for black hole accretion disks, jets, and neutron star magnetospheres. We apply the new scheme to study the effect of resistivity on accreting black holes, accounting for dissipative effects as reconnection.

The Kelvin-Helmholtz instability (KHI) is a known mechanism for penetration of solar wind matter into the magnetosphere. Using three-dimensional, resistive magnetohydrodynamic simulations, the double mid-latitude reconnection (DMLR) process was shown to efficiently exchange solar wind matter into the magnetosphere, through mixing and reconnection. Here, we compute test particle orbits through DMLR configurations. In the instantaneous electromagnetic fields, charged particle trajectories are integrated using the guiding centre approximation. The mechanisms involved in the electron particle orbits and their kinetic energy evolutions are studied in detail, to identify specific signatures of the DMLR through particle characteristics. The charged particle orbits are influenced mainly by magnetic curvature drifts. We identify complex, temporarily trapped, trajectories where the combined electric field and (reconnected) magnetic field variations realize local cavities where particles gain energy before escaping. By comparing the orbits in strongly deformed fields due to the KHI development, with the textbook mirror-drift orbits resulting from our initial configuration, we identify effects due to current sheets formed in the DMLR process. We do this in various representative stages during the DMLR development.

Direct observations of compact objects, in the form of radiation spectra, gravitational waves from VIRGO/LIGO, and forthcoming direct imaging, are currently one of the primary source of information on the physics of plasmas in extreme astrophysical environments. The modeling of such physical phenomena requires numerical methods that allow for the simulation of microscopic plasma dynamics in presence of both strong gravity and electromagnetic fields. In Bacchini et al. (2018) we presented a detailed study on numerical techniques for the integration of free geodesic motion. Here we extend the study by introducing electromagnetic forces in the simulation of charged particles in curved spacetimes. We extend the Hamiltonian energy-conserving method presented in Bacchini et al. (2018) to include the Lorentz force and we test its performance compared to that of standard explicit Runge-Kutta and implicit midpoint rule schemes against analytic solutions. Then, we show the application of the numerical schemes to the integration of test particle trajectories in general relativistic magnetohydrodynamic (GRMHD) simulations, by modifying the algorithms to handle grid-based electromagnetic fields. We test this approach by simulating ensembles of charged particles in a static GRMHD configuration obtained with the Black Hole Accretion Code (BHAC).

We analyze particle acceleration in explosive reconnection events in magnetically dominated proton-electron plasmas. Reconnection is driven by large-scale magnetic stresses in interacting current-carrying flux tubes. Our model relies on development of current-driven instabilities on macroscopic scales. These tilt-kink instabilities develop in an initially force-free equilibrium of repelling current channels. Using MHD methods we study a 3D model of repelling and interacting flux tubes in which we simultaneously evolve test particles, guided by electromagnetic fields obtained from MHD. We identify two stages of particle acceleration; Initially particles accelerate in the current channels, after which the flux ropes start tilting and kinking and particles accelerate due to reconnection processes in the plasma. The explosive stage of reconnection produces non-thermal energy distributions with slopes that depend on plasma resistivity and the initial particle velocity. We also discuss the influence of the length of the flux ropes on particle acceleration and energy distributions. This study extends previous 2.5D results to 3D setups, providing all ingredients needed to model realistic scenarios like solar flares, black hole flares and particle acceleration in pulsar wind nebulae: formation of strong resistive electric fields, explosive reconnection and non-thermal particle distributions. By assuming initial energy equipartition between electrons and protons, applying low resistivity in accordance with solar corona conditions and limiting the flux rope length to a fraction of a solar radius we obtain realistic energy distributions for solar flares with non-thermal power law tails and maximum electron energies up to 11 MeV and maximum proton energies up to 1 GeV.

We compare relativistic particle integrators commonly used in plasma physics showing several test cases relevant for astrophysics. Three explicit particle pushers are considered, namely the Boris, Vay, and Higuera-Cary schemes. We also present a new relativistic fully implicit particle integrator that is energy conserving. Furthermore, a method based on the relativistic guiding center approximation is included. The algorithms are described such that they can be readily implemented in magnetohydrodynamics codes or Particle-in-Cell codes. Our comparison focuses on the strengths and key features of the particle integrators. We test the conservation of invariants of motion, and the accuracy of particle drift dynamics in highly relativistic, mildly relativistic, and non-relativistic settings. The methods are compared in idealized test cases, i.e., without considering feedback on the electrodynamic fields, collisions, pair creation, or radiation. The test cases include uniform electric and magnetic fields, $\mathbf{E}\times\mathbf{B}$-fields, force-free fields, and setups relevant for high-energy astrophysics, e.g., a magnetic mirror, a magnetic dipole, and a magnetic null. These tests have direct relevance for particle acceleration in shocks and in magnetic reconnection.

We apply the general relativistic resistive magnetohydrodynamics code {\tt BHAC} to perform a 2D study of the formation and evolution of a reconnection layer in between two merging magnetic flux tubes in Minkowski spacetime. Small-scale effects in the regime of low resistivity most relevant for dilute astrophysical plasmas are resolved with very high accuracy due to the extreme resolutions obtained with adaptive mesh refinement. Numerical convergence in the highly nonlinear plasmoid-dominated regime is confirmed for a sweep of resolutions. We employ both uniform resistivity and non-uniform resistivity based on the local, instantaneous current density. For uniform resistivity we find Sweet-Parker reconnection, from $\eta = 10^{-2}$ down to $\eta = 10^{-4}$, for a reference case of magnetisation $\sigma = 3.33$ and plasma-$\beta = 0.1$. {For uniform resistivity $\eta=5\times10^{-5}$ the tearing mode is recovered, resulting in the formation of secondary plasmoids. The plasmoid instability enhances the reconnection rate to $v_{\rm rec} \sim 0.03c$ compared to $v_{\rm rec} \sim 0.01c$ for $\eta=10^{-4}$.} For non-uniform resistivity with a base level $\eta_0 = 10^{-4}$ and an enhanced current-dependent resistivity in the current sheet, we find an increased reconnection rate of $v_{\rm rec} \sim 0.1c$. The influence of the magnetisation $\sigma$ and the plasma-$\beta$ is analysed for cases with uniform resistivity $\eta=5\times10^{-5}$ and $\eta=10^{-4}$ in a range $0.5 \leq \sigma \leq 10$ and $0.01 \leq \beta \leq 1$ in regimes that are applicable for black hole accretion disks and jets. The plasmoid instability is triggered for Lundquist numbers larger than a critical value of $S_{\rm c} \approx 8000$.

The Black Hole Accretion Code (BHAC) has recently been extended with the ability to evolve charged test particles according to the Lorentz force within resistive relativistic magnetohydrodynamics simulations. We apply this method to evolve particles in a reconnecting current sheet that forms due to the coalescence of two magnetic flux tubes in 2D Minkowski spacetime. This is the first analysis of charged test particle evolution in resistive relativistic magnetohydrodynamics simulations. The energy distributions of an ensemble of 100.000 electrons are analyzed, as well as the acceleration of particles in the plasmoids that form in the reconnection layer. The effect of the Lundquist number, magnetization, and plasma-$\beta$ on the particle energy distribution is explored for a range of astrophysically relevant parameters. We find that electrons accelerate to non-thermal energies in the thin current sheets in all cases. We find two separate acceleration regimes: An exponential increase of the Lorentz factor during the island coalescence where the acceleration depends linearly on the resistivity and a nonlinear phase with high variability. These results are relevant for determining energy distributions and acceleration sites obtaining radiation maps in large-scale magnetohydrodynamics simulations of black hole accretion disks and jets.

We present recent developments on numerical algorithms for computing photon and particle trajectories in the surrounding of compact objects. Strong gravity around neutron stars or black holes causes relativistic effects on the motion of massive particles and distorts light rays due to gravitational lensing. Efficient numerical methods are required for solving the equations of motion and compute i) the black hole shadow obtained by tracing light rays from the object to a distant observer, and ii) obtain information on the dynamics of the plasma at the microscopic scale. Here, we present generalized algorithms capable of simulating ensembles of photons or massive particles in any spacetime, with the option of including external forces. The coupling of these tools with GRMHD simulations is the key point for obtaining insight on the complex dynamics of accretion disks and jets and for comparing simulations with upcoming observational results from the Event Horizon Telescope.

The numerical integration of particle trajectories in curved spacetimes is fundamental for obtaining realistic models of the particle dynamics around massive compact objects such as black holes and neutron stars. Generalized algorithms capable of handling generic metrics are required for studies of both standard spacetimes (Schwarzschild and Kerr metrics) and non-standard spacetimes (e.g. Schwarzschild metric plus non-classical perturbations or multiple black hole metrics). The most commonly employed explicit numerical schemes (e.g. Runge-Kutta) are incapable of producing highly accurate results at critical points, e.g. in the regions close to the event horizon where gravity causes extreme curvature of the spacetime, at an acceptable computational cost. Here, we describe a generalized algorithm for the numerical integration of time-like (massive particles) and null (photons) geodesics in any given 3+1 split spacetime. We introduce a new, exactly energy-conserving implicit integration scheme based on the preservation of the underlying Hamiltonian, and we compare its properties with a standard fourth-order Runge-Kutta explicit scheme and an implicit midpoint scheme. We test the numerical performance of the three schemes against analytic solutions of particle and photon orbits in Schwarzschild and Kerr spacetimes. We also prove the versatility of our framework in handling more exotic metrics such as Morris-Thorne wormholes and quantum-perturbed Schwarzschild black holes. The generalized approach is also discussed in the perspective of future extensions to more complex particle dynamics, e.g. the addition of the Lorentz force acting on charged particles, which allows for test particle diagnostics in GRMHD simulations.

Magnetic reconnection and non-thermal particle distributions associated with current-driven instabilities are investigated by means of resistive magnetohydrodynamics (MHD) simulations combined with relativistic test particle methods. We propose a system with two parallel, repelling current channels in an initially force-free equilibrium, as a simplified representation of flux ropes in a stellar magnetosphere. The current channels undergo a rotation and separation on Alfv\'enic timescales, forming secondary islands and (up to tearing unstable) current sheets in which non-thermal energy distributions are expected to develop. Using the recently developed particle module of our open-source grid-adaptive MPI-AMRVAC software, we simulate MHD evolution combined with test particle treatments in MHD snapshots. We explore under which plasma-$\beta$ conditions the fastest reconnection occurs in two-and-a-half dimensional (2.5D) scenarios and in these settings test particles are evolved. We quantify energy distributions, acceleration mechanisms, relativistic corrections to the particle equations of motion and effects of resistivity in magnetically dominated proton-electron plasmas. Due to large resistive electric fields and indefinite acceleration of particles in the infinitely long current channels, hard energy spectra are found in 2.5D configurations. Solutions to these numerical artifacts are proposed for both 2.5D setups and future 3D work. We discuss the magnetohydrodynamics of an additional kink instability in 3D setups and the expected effects on energy distributions. The obtained results hold as a proof-of-principle for test particle approaches in MHD simulations, relevant to explore less idealised scenarios like solar flares and more exotic astrophysical phenomena, like black hole flares, magnetar magnetospheres and pulsar wind nebulae.