Showing 1–49 of 49 results for author: Ripperda, B

Resolution analysis of magnetically arrested disk simulations

Authors: Leon Sosapanta Salas, Gibwa Musoke, Koushik Chatterjee, Sera Markoff, Oliver Porth, Matthew Liska, Bart Ripperda

Abstract: Polarisation measurements by the Event Horizon Telescope from M87$^{\ast}$ and Sgr A$^\ast$ suggest that there is a dynamically strong, ordered magnetic field, typical of what is expected of a magnetically arrested accretion disk (MAD). In such disks the strong poloidal magnetic field can suppress the accretion flow and cause episodic flux eruptions. Recent work shows that General Relativistic Mag… ▽ More Polarisation measurements by the Event Horizon Telescope from M87$^{\ast}$ and Sgr A$^\ast$ suggest that there is a dynamically strong, ordered magnetic field, typical of what is expected of a magnetically arrested accretion disk (MAD). In such disks the strong poloidal magnetic field can suppress the accretion flow and cause episodic flux eruptions. Recent work shows that General Relativistic Magnetohydrodynamic (GRMHD) MAD simulations feature dynamics of turbulence and mixing instabilities that are becoming resolved at higher resolutions. We perform a convergence study of MADs exceeding the status quo by an order of magnitude in resolution. We use existing 3D simulations performed with the H-AMR code, up to resolution of 5376 x 2304 x 2304 in a logarithmic spherical-polar grid. We find consistent time-averaged disk properties across all resolutions. However, higher resolutions reveal signs of inward angular momentum transport attributed to turbulent convection, particularly evident when mixing instabilities occur at the surfaces of flux tubes during flux eruptions. Additionally, we see wave-like features in the jet sheath, which become more prominent at higher resolutions, that may induce mixing between jet and disk. At higher resolutions, we observe the sheath to be thinner, resulting in increased temperature, reduced magnetisation, and greater variability. Those differences could affect the dissipation of energy, that would eventually result in distinct observable radiative emission from high-resolution simulations. With higher resolutions, we can delve into crucial questions about horizon-scale physics and its impact on the dynamics and emission properties of larger-scale jets. △ Less

Submitted 1 May, 2024; originally announced May 2024.

Comments: 14 pages, 11 figures, submitted to the monthly notices of the royal astronomical society

Black Hole-Disk Interactions in Magnetically Arrested Active Galactic Nuclei: General Relativistic Magnetohydrodynamic Simulations Using A Time-Dependent, Binary Metric

Authors: Sean M. Ressler, Luciano Combi, Xinyu Li, Bart Ripperda, Huan Yang

Abstract: Perturber objects interacting with supermassive black hole accretion disks are often invoked to explain observed quasi-periodic behavior in active galactic nuclei (AGN). We present global, 3D general relativistic magnetohydrodynamic (GRMHD) simulations of black holes on inclined orbits colliding with magnetically arrested thick AGN disks using a binary black hole spacetime with mass ratio $0.1$. W… ▽ More Perturber objects interacting with supermassive black hole accretion disks are often invoked to explain observed quasi-periodic behavior in active galactic nuclei (AGN). We present global, 3D general relativistic magnetohydrodynamic (GRMHD) simulations of black holes on inclined orbits colliding with magnetically arrested thick AGN disks using a binary black hole spacetime with mass ratio $0.1$. We do this by implementing an approximate time-dependent binary black hole metric into the GRMHD code Athena++. The secondary enhances the unbound mass outflow rate 2-4 times above that provided by the disk in quasi-periodic outbursts, eventually merging into a more continuous outflow at larger distances. We present a simple analytic model that qualitatively agrees well with this result and can be used to extrapolate to unexplored regions of parameter space. We show self-consistently for the first time that spin-orbit coupling between the primary black hole spin and the binary orbital angular momentum causes the accretion disk and jet directions to precess significantly (by 60$^\circ$-80$^\circ$) on long time-scales (e.g., $\sim$ 20 times the binary orbital period). Because this effect may be the only way for thick AGN disks to consistently precess, it could provide strong evidence of a secondary black hole companion if observed in such a system. Besides this new phenomenology, the time-average properties of the disk and accretion rates onto the primary are only marginally altered by the presence of the secondary, consistent with our estimate for a perturbed thick disk. This situation might drastically change in cooled thin disks. △ Less

Submitted 2 April, 2024; originally announced April 2024.

Comments: Accepted by ApJ. Animations are available at https://www.youtube.com/playlist?list=PL3pLmTeUPcqSBQM0OvfkR47opm8ucjA0_

Monster shocks, gamma-ray bursts and black hole quasi-normal modes from neutron-star collapse

Authors: Elias R. Most, Andrei M. Beloborodov, Bart Ripperda

Abstract: We perform the first magnetohydrodynamic simulation tracking the magnetosphere of a collapsing magnetar. The collapse is expected for massive rotating magnetars formed in merger events, and it may occur many hours after the merger. The results suggest a novel gamma-ray burst (GRB) scenario, which creates a delayed high-energy counterpart of the merger gravitational waves. The collapse launches an… ▽ More We perform the first magnetohydrodynamic simulation tracking the magnetosphere of a collapsing magnetar. The collapse is expected for massive rotating magnetars formed in merger events, and it may occur many hours after the merger. The results suggest a novel gamma-ray burst (GRB) scenario, which creates a delayed high-energy counterpart of the merger gravitational waves. The collapse launches an outgoing magnetospheric shock, and a hot magnetized outflow forms behind the shock. The outflow is modulated by the ring-down of the nascent black hole, imprinting its kilohertz quasi-normal modes on the GRB tail. △ Less

Submitted 1 April, 2024; originally announced April 2024.

Comments: 7 pages, 4 figures

Current sheet alignment in oblique black hole magnetospheres -- a black hole pulsar?

Authors: Sebastiaan Selvi, Oliver Porth, Bart Ripperda, Lorenzo Sironi

Abstract: We study the magnetospheric evolution of a non-accreting spinning black hole (BH) with an initially inclined split monopole magnetic field by means of three-dimensional general relativistic magnetohydrodynamics simulations. This serves as a model for a neutron star (NS) collapse or a BH-NS merger remnant after the inherited magnetosphere has settled into a split monopole field creating a striped w… ▽ More We study the magnetospheric evolution of a non-accreting spinning black hole (BH) with an initially inclined split monopole magnetic field by means of three-dimensional general relativistic magnetohydrodynamics simulations. This serves as a model for a neutron star (NS) collapse or a BH-NS merger remnant after the inherited magnetosphere has settled into a split monopole field creating a striped wind. We show that the initially inclined split monopolar current sheet aligns over time with the BH equatorial plane. The inclination angle evolves exponentially towards alignment, with an alignment timescale that is inversely proportional to the square of the BH angular velocity, where higher spin results in faster alignment. Furthermore, magnetic reconnection in the current sheet leads to exponential decay of event horizon penetrating magnetic flux with nearly the same timescale for all considered BH spins. In addition, we present relations for the BH mass and spin in terms of the period and alignment timescale of the striped wind. The explored scenario of a rotating, aligning and reconnecting current sheet can potentially lead to multimessenger electromagnetic counterparts to a gravitational wave event due to the acceleration of particles powering high-energy radiation, plasmoid mergers resulting in coherent radio signals, and pulsating emission due to the initial misalignment of the BH magnetosphere. △ Less

Submitted 25 February, 2024; originally announced February 2024.

Comments: 16 pages, 9 figures

A Case for a Binary Black Hole System Revealed via Quasi-Periodic Outflows

Authors: Dheeraj R. Pasham, Francesco Tombesi, Petra Sukova, Michal Zajacek, Suvendu Rakshit, Eric Coughlin, Peter Kosec, Vladimir Karas, Megan Masterson, Andrew Mummery, Thomas W. -S. Holoien, Muryel Guolo, Jason Hinkle, Bart Ripperda, Vojtech Witzany, Ben Shappee, Erin Kara, Assaf Horesh, Sjoert van Velzen, Itai Sfaradi, David L. Kaplan, Noam Burger, Tara Murphy, Ronald Remillard, James F. Steiner , et al. (11 additional authors not shown)

Abstract: Binaries containing a compact object orbiting a supermassive black hole are thought to be precursors of gravitational wave events, but their identification has been extremely challenging. Here, we report quasi-periodic variability in X-ray absorption which we interpret as quasi-periodic outflows (QPOuts) from a previously low-luminosity active galactic nucleus after an outburst, likely caused by a… ▽ More Binaries containing a compact object orbiting a supermassive black hole are thought to be precursors of gravitational wave events, but their identification has been extremely challenging. Here, we report quasi-periodic variability in X-ray absorption which we interpret as quasi-periodic outflows (QPOuts) from a previously low-luminosity active galactic nucleus after an outburst, likely caused by a stellar tidal disruption. We rule out several models based on observed properties and instead show using general relativistic magnetohydrodynamic simulations that QPOuts, separated by roughly 8.3 days, can be explained with an intermediate-mass black hole secondary on a mildly eccentric orbit at a mean distance of about 100 gravitational radii from the primary. Our work suggests that QPOuts could be a new way to identify intermediate/extreme-mass ratio binary candidates. △ Less

Submitted 15 February, 2024; originally announced February 2024.

Comments: Accepted for publication in Science Advances. We report a new supermassive black hole phenomenon that we call quasi-periodic outflows (QPOuts)

A Beginner's Guide to Black Hole Imaging and Associated Tests of General Relativity

Authors: Alexandru Lupsasca, Daniel R. Mayerson, Bart Ripperda, Seppe Staelens

Abstract: Following the 2019 release by the Event Horizon Telescope Collaboration of the first pictures of a supermassive black hole, there has been an explosion of interest in black hole images, their theoretical interpretation, and their potential use in tests of general relativity. The literature on the subject has now become so vast that an introductory guide for newcomers would appear welcome. Here, we… ▽ More Following the 2019 release by the Event Horizon Telescope Collaboration of the first pictures of a supermassive black hole, there has been an explosion of interest in black hole images, their theoretical interpretation, and their potential use in tests of general relativity. The literature on the subject has now become so vast that an introductory guide for newcomers would appear welcome. Here, we aim to provide an accessible entry point to this growing field, with a particular focus on the black hole "photon ring": the bright, narrow ring of light that dominates images of a black hole and belongs to the black hole itself, rather than to its surrounding plasma. Far from an exhaustive review, this beginner's guide offers a pedagogical review of the key basic concepts and a brief summary of some results at the research frontier. △ Less

Submitted 2 February, 2024; originally announced February 2024.

Comments: 45 pages (+ appendix and references). Invited chapter for the forthcoming book "Recent Progress on Gravity Tests" (eds. C. Bambi and A. Cárdenas-Avendaño, Springer Singapore)

Ordered magnetic fields around the 3C 84 central black hole

Authors: G. F. Paraschos, J. -Y. Kim, M. Wielgus, J. Röder, T. P. Krichbaum, E. Ros, I. Agudo, I. Myserlis, M. Moscibrodzka, E. Traianou, J. A. Zensus, L. Blackburn, C. -K. Chan, S. Issaoun, M. Janssen, M. D. Johnson, V. L. Fish, K. Akiyama, A. Alberdi, W. Alef, J. C. Algaba, R. Anantua, K. Asada, R. Azulay, U. Bach , et al. (258 additional authors not shown)

Abstract: 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures a… ▽ More 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures and understand the physical conditions in the compact region of 3C84. We used EHT 228GHz observations and, given the limited (u,v)-coverage, applied geometric model fitting to the data. We also employed quasi-simultaneously observed, multi-frequency VLBI data for the source in order to carry out a comprehensive analysis of the core structure. We report the detection of a highly ordered, strong magnetic field around the central, SMBH of 3C84. The brightness temperature analysis suggests that the system is in equipartition. We determined a turnover frequency of $ν_m=(113\pm4)$GHz, a corresponding synchrotron self-absorbed magnetic field of $B_{SSA}=(2.9\pm1.6)$G, and an equipartition magnetic field of $B_{eq}=(5.2\pm0.6)$G. Three components are resolved with the highest fractional polarisation detected for this object ($m_\textrm{net}=(17.0\pm3.9)$%). The positions of the components are compatible with those seen in low-frequency VLBI observations since 2017-2018. We report a steeply negative slope of the spectrum at 228GHz. We used these findings to test models of jet formation, propagation, and Faraday rotation in 3C84. The findings of our investigation into different flow geometries and black hole spins support an advection-dominated accretion flow in a magnetically arrested state around a rapidly rotating supermassive black hole as a model of the jet-launching system in the core of 3C84. However, systematic uncertainties due to the limited (u,v)-coverage, however, cannot be ignored. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: 15 pages, 6 figures, published in A&A

Journal ref: Issue: A&A Volume 682, February 2024; Article number: L3; Number of pages: 15

Charged Particle Motion in Neutron Star Magnetic Fields: A Comparison Between the Boris Algorithm and the Guiding Center Approximation

Authors: Minghao Zou, Bart Ripperda

Abstract: Neutron star emission originates typically from its magnetosphere due to radiating electrons. Trajectories of relativistic charged particles under uniform electromagnetic fields can be calculated analytically. However, under more complex and realistic fields, numerical solutions are required. Two common schemes are the Boris method, which solves the full equations of motion, and the guiding center… ▽ More Neutron star emission originates typically from its magnetosphere due to radiating electrons. Trajectories of relativistic charged particles under uniform electromagnetic fields can be calculated analytically. However, under more complex and realistic fields, numerical solutions are required. Two common schemes are the Boris method, which solves the full equations of motion, and the guiding center approximation (GCA), which only evolves the orbital center. We compare both methods in a series of tests, discuss their characteristics and quantify their accuracy. We apply the methods to dipolar, quadrupolar, and quadrudipolar magnetic fields, as applicable for magnetospheres. It is essential to consider such realistic neutron star magnetic field geometries to model the emission from magnetars and pulsars. Our work can assist the Neutron Star Interior Composition ExploreR (NICER) to understand emission from the surface and magnetosphere of neutron stars and to study their composition. △ Less

Submitted 24 November, 2023; originally announced November 2023.

Comments: 5 pages, 4 figures, to be published in IEEE Xplore

Synchrotron polarization signatures of surface waves in supermassive black hole jets

Authors: Jordy Davelaar, Bart Ripperda, Lorenzo Sironi, Alexander A. Philippov, Hector Olivares, Oliver Porth, Bram van den Berg, Thomas Bronzwaer, Koushik Chatterjee, Matthew Liska

Abstract: Supermassive black holes in active galactic nuclei (AGN) are known to launch relativistic jets, which are observed across the entire electromagnetic spectrum and are thought to be efficient particle accelerators. Their primary radiation mechanism for radio emission is polarized synchrotron emission produced by a population of non-thermal electrons. In this Letter, we present a global general relat… ▽ More Supermassive black holes in active galactic nuclei (AGN) are known to launch relativistic jets, which are observed across the entire electromagnetic spectrum and are thought to be efficient particle accelerators. Their primary radiation mechanism for radio emission is polarized synchrotron emission produced by a population of non-thermal electrons. In this Letter, we present a global general relativistic magnetohydrodynamical (GRMHD) simulation of a magnetically arrested disk (MAD). After the simulation reaches the MAD state, we show that waves are continuously launched from the vicinity of the black hole and propagate along the interface between the jet and the wind. At this interface, a steep gradient in velocity is present between the mildly relativistic wind and the highly relativistic jet. The interface is, therefore, a shear layer, and due to the shear, the waves generate roll-ups that alter the magnetic field configuration and the shear layer geometry. We then perform polarized radiation transfer calculations of our GRMHD simulation and find signatures of the waves in both total intensity and linear polarization, effectively lowering the fully resolved polarization fraction. The tell-tale polarization signatures of the waves could be observable by future Very Long Baseline Interferometric observations, e.g., by the next-generation Event Horizon Telescope. △ Less

Submitted 24 October, 2023; v1 submitted 14 September, 2023; originally announced September 2023.

Comments: 20 pages, 17 figures, accepted for publication in ApJL

A search for pulsars around Sgr A* in the first Event Horizon Telescope dataset

Authors: Pablo Torne, Kuo Liu, Ralph P. Eatough, Jompoj Wongphechauxsorn, James M. Cordes, Gregory Desvignes, Mariafelicia De Laurentis, Michael Kramer, Scott M. Ransom, Shami Chatterjee, Robert Wharton, Ramesh Karuppusamy, Lindy Blackburn, Michael Janssen, Chi-kwan Chan, Geoffrey B. Crew, Lynn D. Matthews, Ciriaco Goddi, Helge Rottmann, Jan Wagner, Salvador Sanchez, Ignacio Ruiz, Federico Abbate, Geoffrey C. Bower, Juan J. Salamanca , et al. (261 additional authors not shown)

Abstract: The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission… ▽ More The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission spectra - are expected to be very faint. However, it also negates pulse scattering, an effect that could hinder pulsar detections in the Galactic Center. Additionally, magnetars or a secondary inverse Compton emission could be stronger at millimeter wavelengths than at lower frequencies. We present a search for pulsars close to Sgr A* using the data from the three most-sensitive stations in the EHT 2017 campaign: the Atacama Large Millimeter/submillimeter Array, the Large Millimeter Telescope and the IRAM 30 m Telescope. We apply three detection methods based on Fourier-domain analysis, the Fast-Folding-Algorithm and single pulse search targeting both pulsars and burst-like transient emission; using the simultaneity of the observations to confirm potential candidates. No new pulsars or significant bursts were found. Being the first pulsar search ever carried out at such high radio frequencies, we detail our analysis methods and give a detailed estimation of the sensitivity of the search. We conclude that the EHT 2017 observations are only sensitive to a small fraction ($\lesssim$2.2%) of the pulsars that may exist close to Sgr A*, motivating further searches for fainter pulsars in the region. △ Less

Submitted 29 August, 2023; originally announced August 2023.

Comments: 33 pages, 7 figures, 6 Tables. Accepted for publication in ApJ

Neutron Star Atmosphere-Ocean Dynamics

Authors: Joonas Nättilä, James Y-K. Cho, Jack W. Skinner, Elias R. Most, Bart Ripperda

Abstract: We analyze the structure and dynamics of the plasma atmospheres and Coulomb-liquid oceans on neutron stars. Salient dynamical parameters are identified and their values estimated for the governing set of magnetohydrodynamics equations. Neutron star atmospheres and oceans are strongly stratified and, depending on the rotation period, contain a multitude of long-lived vortices (spots) and/or narrow… ▽ More We analyze the structure and dynamics of the plasma atmospheres and Coulomb-liquid oceans on neutron stars. Salient dynamical parameters are identified and their values estimated for the governing set of magnetohydrodynamics equations. Neutron star atmospheres and oceans are strongly stratified and, depending on the rotation period, contain a multitude of long-lived vortices (spots) and/or narrow zonal jets (free-shear zones) in the large plasma-beta regime - i.e., $β_p \gg 1$ (hydrodynamic regime). In contrast, when $β_p \lesssim 1$ (magnetohydrodynamic regime), the flow is dominated by a global lattice of effectively fixed magnetic islands (plasmoids) - without any jets. Understanding the spatio-temporal variability of dynamic atmospheres and oceans on neutron stars is crucial for interpreting observations of their X-ray emissions. △ Less

Submitted 5 June, 2024; v1 submitted 13 June, 2023; originally announced June 2023.

Comments: 15 pages, 6 figures, 2 appendices. Accepted to ApJ

Cosmic ray transport in large-amplitude turbulence with small-scale field reversals

Authors: Philipp Kempski, Drummond B. Fielding, Eliot Quataert, Alisa K. Galishnikova, Matthew W. Kunz, Alexander A. Philippov, Bart Ripperda

Abstract: The nature of cosmic ray (CR) transport in the Milky Way remains elusive. The predictions of current micro-physical CR transport models in magneto-hydrodynamic (MHD) turbulence are drastically different from what is observed. These models usually focus on MHD turbulence with a strong guide field and ignore the impact of turbulent intermittency on particle propagation. This motivates our studying t… ▽ More The nature of cosmic ray (CR) transport in the Milky Way remains elusive. The predictions of current micro-physical CR transport models in magneto-hydrodynamic (MHD) turbulence are drastically different from what is observed. These models usually focus on MHD turbulence with a strong guide field and ignore the impact of turbulent intermittency on particle propagation. This motivates our studying the alternative regime of large-amplitude turbulence with $δB/B_0 \gg 1$, in which intermittent small-scale magnetic field reversals are ubiquitous. We study particle transport in such turbulence by integrating trajectories in stationary snapshots. To quantify spatial diffusion, we use a setup with continuous particle injection and escape, which we term the turbulent leaky box. We find that particle transport is very different from the strong-guide-field case. Low-energy particles are better confined than high-energy particles, despite less efficient pitch-angle isotropization at small energies. In the limit of weak guide field, energy-dependent confinement is driven by the energy-dependent (in)ability to follow reversing magnetic field lines exactly and by the scattering in regions of ``resonant curvature", where the field line bends on a scale that is of order the local particle gyro-radius. We derive a heuristic model of particle transport in magnetic folds that approximately reproduces the energy dependence of transport found numerically. We speculate that CR propagation in the Galaxy is regulated by the intermittent field reversals highlighted here and discuss the implications of our findings for CR transport in the Milky Way. △ Less

Submitted 19 December, 2023; v1 submitted 24 April, 2023; originally announced April 2023.

Comments: Published in MNRAS, 16 pages, 9 Figures

Comparison of Polarized Radiative Transfer Codes used by the EHT Collaboration

Authors: Ben S. Prather, Jason Dexter, Monika Moscibrodzka, Hung-Yi Pu, Thomas Bronzwaer, Jordy Davelaar, Ziri Younsi, Charles F. Gammie, Roman Gold, George N. Wong, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Uwe Bach, Anne-Kathrin Baczko, David Ball, Mislav Baloković, John Barrett, Michi Bauböck, Bradford A. Benson, Dan Bintley , et al. (248 additional authors not shown)

Abstract: Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curve… ▽ More Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I, Q, U , and V respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties. △ Less

Submitted 21 March, 2023; originally announced March 2023.

Comments: Accepted for publication in ApJ

Black Hole Photon Rings Beyond General Relativity

Authors: Seppe Staelens, Daniel R. Mayerson, Fabio Bacchini, Bart Ripperda, Lorenzo Küchler

Abstract: We investigate whether photon ring observations in black hole imaging are able to distinguish between the Kerr black hole in general relativity and alternative black holes that deviate from Kerr. Certain aspects of photon rings have been argued to be robust observables in Very-Long-Baseline Interferometry (VLBI) black hole observations which carry imprints of the underlying spacetime. The photon r… ▽ More We investigate whether photon ring observations in black hole imaging are able to distinguish between the Kerr black hole in general relativity and alternative black holes that deviate from Kerr. Certain aspects of photon rings have been argued to be robust observables in Very-Long-Baseline Interferometry (VLBI) black hole observations which carry imprints of the underlying spacetime. The photon ring shape, as well as its Lyapunov exponent (which encodes the narrowing of successive photon subrings), are detailed probes of the underlying geometry; measurements thereof have been argued to provide a strong null test of general relativity and the Kerr metric. However, a more complicated question is whether such observations of the photon ring properties can distinguish between Kerr and alternative black holes. We provide a first answer to this question by calculating photon rings of the Johannsen, Rasheed-Larsen, and Manko-Novikov black holes. We find that large deviations from Kerr and large observer inclinations are needed to obtain measurable differences in the photon ring shape. In other words, the Kerr photon ring shape appears to be the universal shape even for deviating black holes at low inclinations. On the other hand, the Lyapunov exponent shows more marked variations for deviations from the Kerr metric. Our analysis lays out the groundwork to determine deviations from the Kerr spacetime in photon rings that are potentially detectable by future observing missions. △ Less

Submitted 3 March, 2023; originally announced March 2023.

Comments: 31 pages + appendices; 20 figures

Three-dimensional dynamics of strongly twisted magnetar magnetospheres: Kinking flux tubes and global eruptions

Authors: J. F. Mahlmann, A. A. Philippov, V. Mewes, B. Ripperda, E. R. Most, L. Sironi

Abstract: The origins of the various outbursts of hard X-rays from magnetars (highly magnetized neutron stars) are still unknown. We identify instabilities in relativistic magnetospheres that can explain a range of X-ray flare luminosities. Crustal surface motions can twist the magnetar magnetosphere by shifting the frozen-in footpoints of magnetic field lines in current-carrying flux bundles. Axisymmetric… ▽ More The origins of the various outbursts of hard X-rays from magnetars (highly magnetized neutron stars) are still unknown. We identify instabilities in relativistic magnetospheres that can explain a range of X-ray flare luminosities. Crustal surface motions can twist the magnetar magnetosphere by shifting the frozen-in footpoints of magnetic field lines in current-carrying flux bundles. Axisymmetric (2D) magnetospheres exhibit strong eruptive dynamics, as to say, catastrophic lateral instabilities triggered by a critical footpoint displacement of $ψ_{\rm crit}\gtrsimπ$. In contrast, our new three-dimensional (3D) twist models with finite surface extension capture important non-axisymmetric dynamics of twisted force-free flux bundles in dipolar magnetospheres. Besides the well-established global eruption resulting (as in 2D) from lateral instabilities, such 3D structures can develop helical, kink-like dynamics, and dissipate energy locally (confined eruptions). Up to $25\%$ of the induced twist energy is dissipated and available to power X-ray flares in powerful global eruptions, with most of our models showing an energy release in the range of the most common X-ray outbursts, $\lesssim 10^{43}$erg. Such events occur when significant energy builds up deeply buried in the dipole magnetosphere. Less energetic outbursts likely precede powerful flares, due to intermittent instabilities and confined eruptions of a continuously twisting flux tube. Upon reaching a critical state, global eruptions produce the necessary Poynting-flux-dominated outflows required by models prescribing the fast radio burst production in the magnetar wind -- for example, via relativistic magnetic reconnection or shocks. △ Less

Submitted 25 April, 2023; v1 submitted 14 February, 2023; originally announced February 2023.

Comments: 21 pages, 11 figures, Accepted for publication by ApJL

Particle acceleration by magnetic Rayleigh-Taylor instability: mechanism for flares in black-hole accretion flows

Authors: Vladimir Zhdankin, Bart Ripperda, Alexander A. Philippov

Abstract: We study the magnetic Rayleigh-Taylor instability in relativistic collisionless plasma, as an astrophysical process for nonthermal particle acceleration. We consider dense plasma on top of a highly magnetized cavity with sheared magnetic field. Using particle-in-cell simulations, we show that small plumes grow and merge progressively to form a large-scale plume, which broadens to drive rapid magne… ▽ More We study the magnetic Rayleigh-Taylor instability in relativistic collisionless plasma, as an astrophysical process for nonthermal particle acceleration. We consider dense plasma on top of a highly magnetized cavity with sheared magnetic field. Using particle-in-cell simulations, we show that small plumes grow and merge progressively to form a large-scale plume, which broadens to drive rapid magnetic reconnection in the cavity. We find that this leads to efficient particle acceleration capable of explaining flares from the inner accretion flow onto the black hole Sgr A*. △ Less

Submitted 10 October, 2023; v1 submitted 10 February, 2023; originally announced February 2023.

Comments: 10 pages, 7 figures; accepted for publication in Physical Review Research

Journal ref: Phys. Rev. Research 5, 043023 (2023)

Millimeter Observational Signatures of Flares in Magnetically Arrested Black Hole Accretion Models

Authors: He Jia, Bart Ripperda, Eliot Quataert, Christopher J. White, Koushik Chatterjee, Alexander Philippov, Matthew Liska

Abstract: In general relativistic magneto-hydrodynamic (GRMHD) simulations, accreted magnetic flux on the black hole horizon episodically decays, during which magnetic reconnection heats up the plasma near the horizon, potentially powering high-energy flares like those observed in M87* and Sgr A*. We study the mm observational counterparts of such flaring episodes. The change in 230 GHz flux during the expe… ▽ More In general relativistic magneto-hydrodynamic (GRMHD) simulations, accreted magnetic flux on the black hole horizon episodically decays, during which magnetic reconnection heats up the plasma near the horizon, potentially powering high-energy flares like those observed in M87* and Sgr A*. We study the mm observational counterparts of such flaring episodes. The change in 230 GHz flux during the expected high energy flares depends primarily on the efficiency of accelerating $γ\gtrsim 100$ ($T_e \gtrsim 10^{11}$ K) electrons. For models in which the electrons are heated to $T_e \sim 10^{11}$ K during flares, the hot plasma produced by reconnection significantly enhances 230 GHz emission and increases the size of the 230 GHz image. By contrast, for models in which the electrons are heated to higher temperatures (which we argue are better motivated), the reconnection-heated plasma is too hot to produce significant 230 GHz synchrotron emission, and the 230 GHz flux decreases during high energy flares. We do not find a significant change in the mm polarization during flares as long as the emission is Faraday thin. We also present expectations for the ring-shaped image as observed by the Event Horizon Telescope during flares, as well as multi-wavelength synchrotron spectra. Our results highlight several limitations of standard post-processing prescriptions for the electron temperature in GRMHD simulations. We also discuss the implications of our results for current and future observations of flares in Sgr A*, M87*, and related systems. Appendices contain detailed convergence studies with respect to resolution and plasma magnetization. △ Less

Submitted 23 September, 2023; v1 submitted 21 January, 2023; originally announced January 2023.

Comments: 11+7 pages, 9+7 figures, 1 table, accepted by MNRAS

Collisionless accretion onto black holes: dynamics and flares

Authors: Alisa Galishnikova, Alexander Philippov, Eliot Quataert, Fabio Bacchini, Kyle Parfrey, Bart Ripperda

Abstract: We study the accretion of collisionless plasma onto a rotating black hole from first principles using axisymmetric general-relativistic particle-in-cell simulations. We carry out a side-by-side comparison of these results to analogous general-relativistic magnetohydrodynamic simulations. Although there are many similarities in the overall flow dynamics, three key differences between the kinetic an… ▽ More We study the accretion of collisionless plasma onto a rotating black hole from first principles using axisymmetric general-relativistic particle-in-cell simulations. We carry out a side-by-side comparison of these results to analogous general-relativistic magnetohydrodynamic simulations. Although there are many similarities in the overall flow dynamics, three key differences between the kinetic and fluid simulations are identified. Magnetic reconnection is more efficient, and rapidly accelerates a nonthermal particle population, in our kinetic approach. In addition, the plasma in the kinetic simulations develops significant departures from thermal equilibrium, including pressure anisotropy that excites kinetic-scale instabilities, and a large field-aligned heat flux near the horizon that approaches the free-streaming value. We discuss the implications of our results for modeling event-horizon scale observations of Sgr A* and M87 by GRAVITY and the Event Horizon Telescope. △ Less

Submitted 5 December, 2022; originally announced December 2022.

Comments: 10 pages, 8 figures, submitted

Tracing the hot spot motion using the next generation Event Horizon Telescope (ngEHT)

Authors: Razieh Emami, Paul Tiede, Sheperd S. Doeleman, Freek Roelofs, Maciek Wielgus, Lindy Blackburn, Matthew Liska, Koushik Chatterjee, Bart Ripperda, Antonio Fuentes, Avery Broderick, Lars Hernquist, Charles Alcock, Ramesh Narayan, Randall Smith, Grant Tremblay, Angelo Ricarte, He Sun, Richard Anantua, Yuri Y. Kovalev, Priyamvada Natarajan, Mark Vogelsberger

Abstract: We propose to trace the dynamical motion of a shearing hot spot near the SgrA* source through a dynamical image reconstruction algorithm, StarWarps. Such a hot spot may form as the exhaust of magnetic reconnection in a current sheet near the black hole horizon. A hot spot that is ejected from the current sheet into an orbit in the accretion disk may shear and diffuse due to instabilities at its bo… ▽ More We propose to trace the dynamical motion of a shearing hot spot near the SgrA* source through a dynamical image reconstruction algorithm, StarWarps. Such a hot spot may form as the exhaust of magnetic reconnection in a current sheet near the black hole horizon. A hot spot that is ejected from the current sheet into an orbit in the accretion disk may shear and diffuse due to instabilities at its boundary during its orbit, resulting in a distinct signature. We subdivide the motion to two distinct phases; the first phase refers to the appearance of the hot spot modelled as a bright blob, followed by a subsequent shearing phase simulated as a stretched ellipse. We employ different observational arrays, including EHT(2017,2022) and the next generation event horizon telescope (ngEHTp1, ngEHT) arrays, in which few new additional sites are added to the observational array. We make dynamical image reconstructions for each of these arrays. Subsequently, we infer the hot spot phase in the first phase followed by the axes ratio and the ellipse area in the second phase. We focus on the direct observability of the orbiting hot spot in the sub-mm wavelength. Our analysis demonstrates that newly added dishes may easily trace the first phase as well as part of the second phase, before the flux is reduced substantially. The algorithm used in this work can be extended to any other types of the dynamical motion. Consequently, we conclude that the ngEHT is a key to directly observe the dynamical motions near variable sources, such as SgrA*. △ Less

Submitted 15 November, 2022; v1 submitted 12 November, 2022; originally announced November 2022.

Comments: 16 pages, 8 pages

Plasmoid Instability in the Multiphase Interstellar Medium

Authors: Drummond B. Fielding, Bart Ripperda, Alexander A. Philippov

Abstract: The processes controlling the complex clump structure, phase distribution, and magnetic field geometry that develops across a broad range of scales in the turbulent interstellar medium remains unclear. Using unprecedentedly high resolution three-dimensional magnetohydrodynamic simulations of thermally unstable turbulent systems, we show that large current sheets unstable to plasmoid-mediated recon… ▽ More The processes controlling the complex clump structure, phase distribution, and magnetic field geometry that develops across a broad range of scales in the turbulent interstellar medium remains unclear. Using unprecedentedly high resolution three-dimensional magnetohydrodynamic simulations of thermally unstable turbulent systems, we show that large current sheets unstable to plasmoid-mediated reconnection form regularly throughout the volume. The plasmoids form in three distinct environments: (i) within cold clumps, (ii) at the asymmetric interface of the cold and warm phases, and (iii) within the warm, volume-filling phase. We then show that the complex magneto-thermal phase structure is characterized by a predominantly highly magnetized cold phase, but that regions of high magnetic curvature, which are the sites of reconnection, span a broad range in temperature. Furthermore, we show that thermal instabilities change the scale dependent anisotropy of the turbulent magnetic field, reducing the increase in eddy elongation at smaller scales. Finally, we show that most of the mass is contained in one contiguous cold structure surrounded by smaller clumps that follow a scale free mass distribution. These clumps tend to be highly elongated and exhibit a size versus internal velocity relation consistent with supersonic turbulence, and a relative clump distance-velocity scaling consistent with subsonic motion. We discuss the striking similarity of cold plasmoids to observed tiny scale atomic and ionized structures and HI fibers, and consider how the prevalence of plasmoids will modify the motion of charged particles thereby impacting cosmic ray transport and thermal conduction in the ISM and other similar systems. △ Less

Submitted 29 April, 2024; v1 submitted 11 November, 2022; originally announced November 2022.

Comments: 21 pages, 13 figures. For associated movies, see https://dfielding14.github.io/movies/

Effective resistivity in relativistic collisionless plasmoid-mediated reconnection

Authors: Sebastiaan Selvi, Oliver Porth, Bart Ripperda, Fabio Bacchini, Lorenzo Sironi, Rony Keppens

Abstract: Magnetic reconnection can power spectacular high-energy astrophysical phenomena by producing non-thermal energy distributions in highly magnetized regions around compact objects. By means of two-dimensional fully kinetic particle-in-cell (PIC) simulations we investigate relativistic collisionless plasmoid-mediated reconnection in magnetically dominated pair plasmas with and without guide field. In… ▽ More Magnetic reconnection can power spectacular high-energy astrophysical phenomena by producing non-thermal energy distributions in highly magnetized regions around compact objects. By means of two-dimensional fully kinetic particle-in-cell (PIC) simulations we investigate relativistic collisionless plasmoid-mediated reconnection in magnetically dominated pair plasmas with and without guide field. In X-points, where diverging flows result in a non-diagonal thermal pressure tensor, a finite residence time for particles gives rise to a localized collisionless effective resistivity. Here, for the first time for relativistic reconnection in a fully developed plasmoid chain we identify the mechanisms driving the non-ideal electric field using a full Ohm's law by means of a statistical analysis based on our PIC simulations. We show that the non-ideal electric field is predominantly driven by gradients of nongyrotropic thermal pressures. We propose a kinetic physics motivated non-uniform effective resistivity model, which is negligible on global scales and becomes significant only locally in X-points, that captures the properties of collisionless reconnection with the aim of mimicking its essentials in non-ideal magnetohydrodynamic descriptions. This effective resistivity model provides a viable opportunity to design physically grounded global models for reconnection-powered high-energy emission. △ Less

Submitted 8 November, 2022; originally announced November 2022.

Comments: 12 pages, 4 figures

Relativistic Signatures of Flux Eruption Events Near Black Holes

Authors: Zachary Gelles, Koushik Chatterjee, Michael Johnson, Bart Ripperda, Matthew Liska

Abstract: Images of supermassive black holes produced using very long baseline interferometry provide a pathway to directly observing effects of a highly curved spacetime, such as a bright ``photon ring'' that arises from strongly lensed emission. In addition, the emission near supermassive black holes is highly variable, with bright high-energy flares regularly observed. We demonstrate that intrinsic varia… ▽ More Images of supermassive black holes produced using very long baseline interferometry provide a pathway to directly observing effects of a highly curved spacetime, such as a bright ``photon ring'' that arises from strongly lensed emission. In addition, the emission near supermassive black holes is highly variable, with bright high-energy flares regularly observed. We demonstrate that intrinsic variability can introduce prominent associated changes in the relative brightness of the photon ring. We analyze both semianalytic toy models and GRMHD simulations with magnetic flux eruption events, showing that they each exhibit a characteristic ``loop'' in the space of relative photon ring brightness versus total flux density. For black holes viewed at high inclination, the relative photon ring brightness can change by an order of magnitude, even with variations in total flux density that are comparatively mild. We show that gravitational lensing, Doppler boosting, and magnetic field structure all significantly affect this feature, and we discuss the prospects for observing it in observations of M87* and Sgr A* with the next-generation Event Horizon Telescope. △ Less

Submitted 24 November, 2022; v1 submitted 13 October, 2022; originally announced October 2022.

Comments: Accepted in Galaxies

Radiative reconnection-powered TeV flares from the black hole magnetosphere in M87

Authors: Hayk Hakobyan, Bart Ripperda, Alexander Philippov

Abstract: Active Galactic Nuclei in general, and the supermassive black hole in M87 in particular, show bright and rapid gamma-ray flares up to energies of 100 GeV and above. For M87, the flares show multiwavelength components, and the variability timescale is comparable to the dynamical time of the event horizon, suggesting that the emission may come from a compact region nearby the nucleus. However, the e… ▽ More Active Galactic Nuclei in general, and the supermassive black hole in M87 in particular, show bright and rapid gamma-ray flares up to energies of 100 GeV and above. For M87, the flares show multiwavelength components, and the variability timescale is comparable to the dynamical time of the event horizon, suggesting that the emission may come from a compact region nearby the nucleus. However, the emission mechanism for these flares is not well understood. Recent high-resolution general-relativistic magnetohydrodynamic simulations show the occurrence of episodic magnetic reconnection events that can power flares nearby the black hole event horizon. In this work, we analyze the radiative properties of the reconnecting current layer under the extreme plasma conditions applicable to the black hole in M87 from the first principles. We show that abundant pair production is expected in the vicinity of the reconnection layer, to the extent that the produced secondary pair plasma dominates the reconnection dynamics. Using analytic estimates backed by two-dimensional particle-in-cell simulations we demonstrate that in the presence of strong synchrotron cooling, reconnection can produce a hard power-law distribution of pair plasma imprinted in the outgoing synchrotron (up to a few tens of MeV) and the inverse-Compton signal (up to TeV). We produce synthetic radiation spectra from our simulations, which can be directly compared with the results of future multiwavelength observations of M87* flares. △ Less

Submitted 17 March, 2023; v1 submitted 5 September, 2022; originally announced September 2022.

Comments: 16 pages, 5 figures, 1 table

Journal ref: ApJL 943 (2023) L29

The Photon Ring in M87*

Authors: Avery E. Broderick, Dominic W. Pesce, Paul Tiede, Hung-Yi Pu, Roman Gold, Richard Anantua, Silke Britzen, Chiara Ceccobello, Koushik Chatterjee, Yongjun Chen, Nicholas S. Conroy, Geoffrey B. Crew, Alejandro Cruz-Osorio, Yuzhu Cui, Sheperd S. Doeleman, Razieh Emami, Joseph Farah, Christian M. Fromm, Peter Galison, Boris Georgiev, Luis C. Ho, David J. James, Britton Jeter, Alejandra Jimenez-Rosales, Jun Yi Koay , et al. (26 additional authors not shown)

Abstract: We report measurements of the gravitationally lensed secondary image -- the first in an infinite series of so-called "photon rings" -- around the supermassive black hole M87* via simultaneous modeling and imaging of the 2017 Event Horizon Telescope (EHT) observations. The inferred ring size remains constant across the seven days of the 2017 EHT observing campaign and is consistent with theoretical… ▽ More We report measurements of the gravitationally lensed secondary image -- the first in an infinite series of so-called "photon rings" -- around the supermassive black hole M87* via simultaneous modeling and imaging of the 2017 Event Horizon Telescope (EHT) observations. The inferred ring size remains constant across the seven days of the 2017 EHT observing campaign and is consistent with theoretical expectations, providing clear evidence that such measurements probe spacetime and a striking confirmation of the models underlying the first set of EHT results. The residual diffuse emission evolves on timescales comparable to one week. We are able to detect with high significance a southwestern extension consistent with that expected from the base of a jet that is rapidly rotating in the clockwise direction. This result adds further support to the identification of the jet in M87* with a black hole spin-driven outflow, launched via the Blandford-Znajek process. We present three revised estimates for the mass of M87* based on identifying the modeled thin ring component with the bright ringlike features seen in simulated images, one of which is only weakly sensitive to the astrophysics of the emission region. All three estimates agree with each other and previously reported values. Our strongest mass constraint combines information from both the ring and the diffuse emission region, which together imply a mass-to-distance ratio of $4.20^{+0.12}_{-0.06}~μ{\rm as}$ and a corresponding black hole mass of $(7.13\pm0.39)\times10^9M_\odot$, where the error on the latter is now dominated by the systematic uncertainty arising from the uncertain distance to M87*. △ Less

Submitted 18 August, 2022; originally announced August 2022.

Journal ref: ApJ, 935, 61 (2022)

Millimeter light curves of Sagittarius A* observed during the 2017 Event Horizon Telescope campaign

Authors: Maciek Wielgus, Nicola Marchili, Ivan Marti-Vidal, Garrett K. Keating, Venkatessh Ramakrishnan, Paul Tiede, Ed Fomalont, Sara Issaoun, Joey Neilsen, Michael A. Nowak, Lindy Blackburn, Charles F. Gammie, Ciriaco Goddi, Daryl Haggard, Daeyoung Lee, Monika Moscibrodzka, Alexandra J. Tetarenko, Geoffrey C. Bower, Chi-Kwan Chan, Koushik Chatterjee, Paul M. Chesler, Jason Dexter, Sheperd S. Doeleman, Boris Georgiev, Mark Gurwell , et al. (6 additional authors not shown)

Abstract: The Event Horizon Telescope (EHT) observed the compact radio source, Sagittarius A* (Sgr A*), in the Galactic Center on 2017 April 5-11 in the 1.3 millimeter wavelength band. At the same time, interferometric array data from the Atacama Large Millimeter/submillimeter Array and the Submillimeter Array were collected, providing Sgr A* light curves simultaneous with the EHT observations. These data s… ▽ More The Event Horizon Telescope (EHT) observed the compact radio source, Sagittarius A* (Sgr A*), in the Galactic Center on 2017 April 5-11 in the 1.3 millimeter wavelength band. At the same time, interferometric array data from the Atacama Large Millimeter/submillimeter Array and the Submillimeter Array were collected, providing Sgr A* light curves simultaneous with the EHT observations. These data sets, complementing the EHT very-long-baseline interferometry, are characterized by a cadence and signal-to-noise ratio previously unattainable for Sgr A* at millimeter wavelengths, and they allow for the investigation of source variability on timescales as short as a minute. While most of the light curves correspond to a low variability state of Sgr A*, the April 11 observations follow an X-ray flare, and exhibit strongly enhanced variability. All of the light curves are consistent with a red noise process, with a power spectral density (PSD) slope measured to be between -2 and -3 on timescales between 1 min and several hours. Our results indicate a steepening of the PSD slope for timescales shorter than 0.3 h. The spectral energy distribution is flat at 220 GHz and there are no time-lags between the 213 and 229 GHz frequency bands, suggesting low optical depth for the event horizon scale source. We characterize Sgr A*'s variability, highlighting the different behavior observed just after the X-ray flare, and use Gaussian process modeling to extract a decorrelation timescale and a PSD slope. We also investigate the systematic calibration uncertainties by analyzing data from independent data reduction pipelines. △ Less

Submitted 14 July, 2022; originally announced July 2022.

Journal ref: The Astrophysical Journal Letters, 930:L19 (2022)

Dynamic alignment and plasmoid formation in relativistic magnetohydrodynamic turbulence

Authors: Alexander Chernoglazov, Bart Ripperda, Alexander A. Philippov

Abstract: We present high resolution 2D and 3D simulations of magnetized decaying turbulence in relativistic resistive magneto-hydrodynamics. The simulations show dynamic formation of large scale intermittent long-lived current sheets being disrupted by the tearing instability into plasmoid chains. These current sheets are locations of enhanced magnetic field dissipation and heating of the plasma. We find m… ▽ More We present high resolution 2D and 3D simulations of magnetized decaying turbulence in relativistic resistive magneto-hydrodynamics. The simulations show dynamic formation of large scale intermittent long-lived current sheets being disrupted by the tearing instability into plasmoid chains. These current sheets are locations of enhanced magnetic field dissipation and heating of the plasma. We find magnetic energy spectra $\propto k^{-3/2}$, together with strongly pronounced dynamic alignment of Elsasser fields and of velocity and magnetic fields, for strong guide-field turbulence, whereas we retrieve spectra $\propto k^{-5/3}$ for the case of a weak guide-field. △ Less

Submitted 15 November, 2021; originally announced November 2021.

Comments: accepted for publication in the ApJ Letters

Event Horizon Telescope observations of the jet launching and collimation in Centaurus A

Authors: Michael Janssen, Heino Falcke, Matthias Kadler, Eduardo Ros, Maciek Wielgus, Kazunori Akiyama, Mislav Baloković, Lindy Blackburn, Katherine L. Bouman, Andrew Chael, Chi-kwan Chan, Koushik Chatterjee, Jordy Davelaar, Philip G. Edwards, Christian M. Fromm, José L. Gómez, Ciriaco Goddi, Sara Issaoun, Michael D. Johnson, Junhan Kim, Jun Yi Koay, Thomas P. Krichbaum, Jun Liu, Elisabetta Liuzzo, Sera Markoff , et al. (215 additional authors not shown)

Abstract: Very-long-baseline interferometry (VLBI) observations of active galactic nuclei at millimeter wavelengths have the power to reveal the launching and initial collimation region of extragalactic radio jets, down to $10-100$ gravitational radii ($r_g=GM/c^2$) scales in nearby sources. Centaurus A is the closest radio-loud source to Earth. It bridges the gap in mass and accretion rate between the supe… ▽ More Very-long-baseline interferometry (VLBI) observations of active galactic nuclei at millimeter wavelengths have the power to reveal the launching and initial collimation region of extragalactic radio jets, down to $10-100$ gravitational radii ($r_g=GM/c^2$) scales in nearby sources. Centaurus A is the closest radio-loud source to Earth. It bridges the gap in mass and accretion rate between the supermassive black holes (SMBHs) in Messier 87 and our galactic center. A large southern declination of $-43^{\circ}$ has however prevented VLBI imaging of Centaurus A below $λ1$cm thus far. Here, we show the millimeter VLBI image of the source, which we obtained with the Event Horizon Telescope at $228$GHz. Compared to previous observations, we image Centaurus A's jet at a tenfold higher frequency and sixteen times sharper resolution and thereby probe sub-lightday structures. We reveal a highly-collimated, asymmetrically edge-brightened jet as well as the fainter counterjet. We find that Centaurus A's source structure resembles the jet in Messier 87 on ${\sim}500r_g$ scales remarkably well. Furthermore, we identify the location of Centaurus A's SMBH with respect to its resolved jet core at $λ1.3$mm and conclude that the source's event horizon shadow should be visible at THz frequencies. This location further supports the universal scale invariance of black holes over a wide range of masses. △ Less

Submitted 5 November, 2021; originally announced November 2021.

Comments: 27 pages, 9 figures. This is a post-peer-review, pre-copyedit version of an article published in Nature Astronomy. The final authenticated version is available online at: https://doi.org/10.1038/s41550-021-01417-w

Journal ref: Nature Astronomy, July 2021, Volume 5, p. 1017-1028

The Variability of the Black-Hole Image in M87 at the Dynamical Time Scale

Authors: Kaushik Satapathy, Dimitrios Psaltis, Feryal Ozel, Lia Medeiros, Sean T. Dougall, Chi-kwan Chan, Maciek Wielgus, Ben S. Prather, George N. Wong, Charles F. Gammie, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Anne-Kathrin Baczko, David R. Ball, Mislav Baloković, John Barrett, Bradford A. Benson, Dan Bintley, Lindy Blackburn, Raymond Blundell , et al. (213 additional authors not shown)

Abstract: The black-hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5-61 days) is comparable to the 6-day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expect… ▽ More The black-hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5-61 days) is comparable to the 6-day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expected structural changes of the images but are free of station-based atmospheric and instrumental errors. We explored the day-to-day variability in closure phase measurements on all six linearly independent non-trivial baseline triangles that can be formed from the 2017 observations. We showed that three triangles exhibit very low day-to-day variability, with a dispersion of $\sim3-5^\circ$. The only triangles that exhibit substantially higher variability ($\sim90-180^\circ$) are the ones with baselines that cross visibility amplitude minima on the $u-v$ plane, as expected from theoretical modeling. We used two sets of General Relativistic magnetohydrodynamic simulations to explore the dependence of the predicted variability on various black-hole and accretion-flow parameters. We found that changing the magnetic field configuration, electron temperature model, or black-hole spin has a marginal effect on the model consistency with the observed level of variability. On the other hand, the most discriminating image characteristic of models is the fractional width of the bright ring of emission. Models that best reproduce the observed small level of variability are characterized by thin ring-like images with structures dominated by gravitational lensing effects and thus least affected by turbulence in the accreting plasmas. △ Less

Submitted 1 November, 2021; originally announced November 2021.

Comments: Accepted for Publication in ApJ

Black hole flares: ejection of accreted magnetic flux through 3D plasmoid-mediated reconnection

Authors: Bart Ripperda, Matthew Liska, Koushik Chatterjee, Gibwa Musoke, Alexander A. Philippov, Sera B. Markoff, Alexander Tchekhovskoy, Ziri Younsi

Abstract: Magnetic reconnection can power bright, rapid flares originating from the inner magnetosphere of accreting black holes. We conduct extremely high resolution ($5376\times2304\times2304$ cells) general-relativistic magnetohydrodynamics simulations, capturing plasmoid-mediated reconnection in a 3D magnetically arrested disk for the first time. We show that an equatorial, plasmoid-unstable current she… ▽ More Magnetic reconnection can power bright, rapid flares originating from the inner magnetosphere of accreting black holes. We conduct extremely high resolution ($5376\times2304\times2304$ cells) general-relativistic magnetohydrodynamics simulations, capturing plasmoid-mediated reconnection in a 3D magnetically arrested disk for the first time. We show that an equatorial, plasmoid-unstable current sheet forms in a transient, non-axisymmetric, low-density magnetosphere within the inner few Schwarzschild radii. Magnetic flux bundles escape from the event horizon through reconnection at the universal plasmoid-mediated rate in this current sheet. The reconnection feeds on the highly-magnetized plasma in the jets and heats the plasma that ends up trapped in flux bundles to temperatures proportional to the jet's magnetization. The escaped flux bundles can complete a full orbit as low-density hot spots, consistent with Sgr A$^{*}$ observations by the GRAVITY interferometer. Reconnection near the horizon produces sufficiently energetic plasma to explain flares from accreting black holes, such as the TeV emission observed from M87. The drop in mass accretion rate during the flare, and the resulting low-density magnetosphere make it easier for very high energy photons produced by reconnection-accelerated particles to escape. The extreme resolution results in a converged plasmoid-mediated reconnection rate that directly determines the timescales and properties of the flare. △ Less

Submitted 6 February, 2022; v1 submitted 30 September, 2021; originally announced September 2021.

Comments: ApJ Letters (2022), 924, L32

Magnetic Hair and Reconnection in Black Hole Magnetospheres

Authors: Ashley Bransgrove, Bart Ripperda, Alexander Philippov

Abstract: The no-hair theorem of general relativity states that isolated black holes are characterized by three parameters: mass, spin, and charge. In this Letter we consider Kerr black holes endowed with highly magnetized plasma-filled magnetospheres. Using general relativistic kinetic plasma and resistive magnetohydrodynamics simulations, we show that a dipole magnetic field on the event horizon opens int… ▽ More The no-hair theorem of general relativity states that isolated black holes are characterized by three parameters: mass, spin, and charge. In this Letter we consider Kerr black holes endowed with highly magnetized plasma-filled magnetospheres. Using general relativistic kinetic plasma and resistive magnetohydrodynamics simulations, we show that a dipole magnetic field on the event horizon opens into a split-monopole and reconnects in a plasmoid-unstable current-sheet. The no-hair theorem is satisfied, in the sense that all components of the stress-energy tensor decay exponentially in time. We measure the decay time of magnetic flux on the event horizon for plasmoid-dominated reconnection in collisionless and collisional plasma. The reconnecting magnetosphere should be a powerful source of hard X-ray emission when the magnetic field is strong. △ Less

Submitted 29 September, 2021; originally announced September 2021.

Comments: 5 pages, 4 figures, accepted PRL

The Polarized Image of a Synchrotron Emitting Ring of Gas Orbiting a Black Hole

Authors: Ramesh Narayan, Daniel C. M. Palumbo, Michael D. Johnson, Zachary Gelles, Elizabeth Himwich, Dominic O. Chang, Angelo Ricarte, Jason Dexter, Charles F. Gammie, Andrew A. Chael, The Event Horizon Telescope Collaboration, :, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Anne-Kathrin Baczko, David Ball, Mislav Balokovic, John Barrett, Bradford A. Benson, Dan Bintley , et al. (215 additional authors not shown)

Abstract: Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equ… ▽ More Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equatorial accretion disk around a Schwarzschild black hole. By using an approximate expression for the null geodesics derived by Beloborodov (2002) and conservation of the Walker-Penrose constant, we provide analytic estimates for the image polarization. We test this model using currently favored general relativistic magnetohydrodynamic simulations of M87*, using ring parameters given by the simulations. For a subset of these with modest Faraday effects, we show that the ring model broadly reproduces the polarimetric image morphology. Our model also predicts the polarization evolution for compact flaring regions, such as those observed from Sgr A* with GRAVITY. With suitably chosen parameters, our simple model can reproduce the EVPA pattern and relative polarized intensity in Event Horizon Telescope images of M87*. Under the physically motivated assumption that the magnetic field trails the fluid velocity, this comparison is consistent with the clockwise rotation inferred from total intensity images. △ Less

Submitted 13 May, 2021; v1 submitted 4 May, 2021; originally announced May 2021.

Comments: 29 pages, 11 figures, published in ApJ on May 3

Journal ref: ApJ 912 35 (2021)

Weak Alfvénic turbulence in relativistic plasmas. Part 1. Dynamical equations and basic dynamics of interacting resonant triads

Authors: J. M. TenBarge, B. Ripperda, A. Chernoglazov, A. Bhattacharjee, J. F. Mahlmann, E. R. Most, J. Juno, Y. Yuan, A. A. Philippov

Abstract: Alfvén wave collisions are the primary building blocks of the non-relativistic turbulence that permeates the heliosphere and low-to-moderate energy astrophysical systems. However, many astrophysical systems such as gamma-ray bursts, pulsar and magnetar magnetospheres, and active galactic nuclei have relativistic flows or energy densities. To better understand these high energy systems, we derive r… ▽ More Alfvén wave collisions are the primary building blocks of the non-relativistic turbulence that permeates the heliosphere and low-to-moderate energy astrophysical systems. However, many astrophysical systems such as gamma-ray bursts, pulsar and magnetar magnetospheres, and active galactic nuclei have relativistic flows or energy densities. To better understand these high energy systems, we derive reduced relativistic MHD equations and employ them to examine weak Alfvénic turbulence, dominated by three-wave interactions, in reduced relativistic magnetohydrodynamics, including the force-free, infinitely magnetized limit. We compare both numerical and analytical solutions to demonstrate that many of the findings from non-relativistic weak turbulence are retained in the relativistic system. But, an important distinction in the relativistic limit is the inapplicability of a formally incompressible limit, i.e, there exists finite coupling to the compressible fast mode regardless of the strength of the magnetic field. Since fast modes can propagate across field lines, this mechanism provides a route for energy to escape strongly magnetized systems, e.g., magnetar magnetospheres. However, we find that the fast-Alfvén coupling is diminished in the limit of oblique propagation. △ Less

Submitted 23 February, 2022; v1 submitted 3 May, 2021; originally announced May 2021.

Comments: Accepted

Weak Alfvénic turbulence in relativistic plasmas II: Current sheets and dissipation

Authors: B. Ripperda, J. F. Mahlmann, A. Chernoglazov, J. M. TenBarge, E. R. Most, J. Juno, Y. Yuan, A. A. Philippov, A. Bhattacharjee

Abstract: Alfvén waves as excited in black hole accretion disks and neutron star magnetospheres are the building blocks of turbulence in relativistic, magnetized plasmas. A large reservoir of magnetic energy is available in these systems, such that the plasma can be heated significantly even in the weak turbulence regime. We perform high-resolution three-dimensional simulations of counter-propagating Alfvén… ▽ More Alfvén waves as excited in black hole accretion disks and neutron star magnetospheres are the building blocks of turbulence in relativistic, magnetized plasmas. A large reservoir of magnetic energy is available in these systems, such that the plasma can be heated significantly even in the weak turbulence regime. We perform high-resolution three-dimensional simulations of counter-propagating Alfvén waves, showing that an $E_{B_{\perp}}(k_{\perp}) \propto k_{\perp}^{-2}$ energy spectrum develops as a result of the weak turbulence cascade in relativistic magnetohydrodynamics and its infinitely magnetized (force-free) limit. The plasma turbulence ubiquitously generates current sheets, which act as locations where magnetic energy dissipates. We show that current sheets form as a natural result of nonlinear interactions between counter-propagating Alfvén waves. These current sheets form due to the compression of elongated eddies, driven by the shear induced by growing higher order modes, and undergo a thinning process until they break-up into small-scale turbulent structures. We explore the formation of {current sheets} both in overlapping waves and in localized wave packet collisions. The relativistic interaction of localized Alfvén waves induces both Alfvén waves and fast waves and efficiently mediates the conversion and dissipation of electromagnetic energy in astrophysical systems. Plasma energization through reconnection in current sheets emerging during the interaction of Alfvén waves can potentially explain X-ray emission in black hole accretion coronae and neutron star magnetospheres. △ Less

Submitted 25 October, 2021; v1 submitted 3 May, 2021; originally announced May 2021.

Comments: Accepted in JPP, accompanied by Paper I: Weak Alfvénic turbulence in relativistic plasmas I: dynamical equations and basic dynamics of interacting resonant triads

Fuzzball Shadows: Emergent Horizons from Microstructure

Authors: Fabio Bacchini, Daniel R. Mayerson, Bart Ripperda, Jordy Davelaar, Héctor Olivares, Thomas Hertog, Bert Vercnocke

Abstract: We study the physical properties of four-dimensional, string-theoretical, horizonless "fuzzball" geometries by imaging their shadows. Their microstructure traps light rays straying near the would-be horizon on long-lived, highly redshifted chaotic orbits. In fuzzballs sufficiently near the scaling limit this creates a shadow much like that of a black hole, while avoiding the paradoxes associated w… ▽ More We study the physical properties of four-dimensional, string-theoretical, horizonless "fuzzball" geometries by imaging their shadows. Their microstructure traps light rays straying near the would-be horizon on long-lived, highly redshifted chaotic orbits. In fuzzballs sufficiently near the scaling limit this creates a shadow much like that of a black hole, while avoiding the paradoxes associated with an event horizon. Observations of the shadow size and residual glow can potentially discriminate between fuzzballs away from the scaling limit and alternative models of black compact objects. △ Less

Submitted 20 October, 2021; v1 submitted 22 March, 2021; originally announced March 2021.

Comments: 6 pages, 3 figures

Gravitational Test Beyond the First Post-Newtonian Order with the Shadow of the M87 Black Hole

Authors: Dimitrios Psaltis, Lia Medeiros, Pierre Christian, Feryal Ozel, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Keiichi Asada, Rebecca Azulay, David Ball, Mislav Balokovic, John Barrett, Dan Bintley, Lindy Blackburn, Wilfred Boland, Geoffrey C. Bower, Michael Bremer, Christiaan D. Brinkerink, Roger Brissenden, Silke Britzen, Dominique Broguiere, Thomas Bronzwaer, Do-Young Byun, John E. Carlstrom, Andrew Chael , et al. (163 additional authors not shown)

Abstract: The 2017 Event Horizon Telescope (EHT) observations of the central source in M87 have led to the first measurement of the size of a black-hole shadow. This observation offers a new and clean gravitational test of the black-hole metric in the strong-field regime. We show analytically that spacetimes that deviate from the Kerr metric but satisfy weak-field tests can lead to large deviations in the p… ▽ More The 2017 Event Horizon Telescope (EHT) observations of the central source in M87 have led to the first measurement of the size of a black-hole shadow. This observation offers a new and clean gravitational test of the black-hole metric in the strong-field regime. We show analytically that spacetimes that deviate from the Kerr metric but satisfy weak-field tests can lead to large deviations in the predicted black-hole shadows that are inconsistent with even the current EHT measurements. We use numerical calculations of regular, parametric, non-Kerr metrics to identify the common characteristic among these different parametrizations that control the predicted shadow size. We show that the shadow-size measurements place significant constraints on deviation parameters that control the second post-Newtonian and higher orders of each metric and are, therefore, inaccessible to weak-field tests. The new constraints are complementary to those imposed by observations of gravitational waves from stellar-mass sources. △ Less

Submitted 2 October, 2020; originally announced October 2020.

Comments: Physical Review Letters

A coupled guiding center-Boris particle pusher for magnetized plasmas in compact-object magnetospheres

Authors: Fabio Bacchini, Bart Ripperda, Alexander A. Philippov, Kyle Parfrey

Abstract: We present a novel numerical scheme for simulating the motion of relativistic charged particles in magnetospheres of compact objects, typically filled with highly magnetized collisionless plasmas. The new algorithm is based on a dynamic switch between the full system of equations of motion and a guiding center approximation. The switch between the two formulations is based on the magnetization of… ▽ More We present a novel numerical scheme for simulating the motion of relativistic charged particles in magnetospheres of compact objects, typically filled with highly magnetized collisionless plasmas. The new algorithm is based on a dynamic switch between the full system of equations of motion and a guiding center approximation. The switch between the two formulations is based on the magnetization of the plasma particles, such that the dynamics are accurately captured by the guiding center motion even when the gyro-frequency is under-resolved by the time step. For particles with a large gyro-radius, due to acceleration in, e.g., reconnecting current sheets, the algorithm adaptively switches to solve the full equations of motion instead. The new scheme is directly compatible with standard Particle-in-Cell codes, and is readily applicable in curved spacetimes via a dedicated covariant formulation. We test the performance of the coupled algorithm by evolving charged particles in electromagnetic configurations of reconnecting current sheets in magnetized plasma, obtained from special- and general-relativistic Particle-in-Cell simulations. The new coupled pusher is capable of producing highly accurate particle trajectories even when the time step is many orders of magnitude larger than the gyro-period, substantially reducing the restrictions of the temporal resolution. △ Less

Submitted 3 November, 2020; v1 submitted 28 July, 2020; originally announced July 2020.

SYMBA: An end-to-end VLBI synthetic data generation pipeline

Authors: F. Roelofs, M. Janssen, I. Natarajan, R. Deane, J. Davelaar, H. Olivares, O. Porth, S. N. Paine, K. L. Bouman, R. P. J. Tilanus, I. M. van Bemmel, H. Falcke, K. Akiyama, A. Alberdi, W. Alef, K. Asada, R. Azulay, A. Baczko, D. Ball, M. Baloković, J. Barrett, D. Bintley, L. Blackburn, W. Boland, G. C. Bower , et al. (183 additional authors not shown)

Abstract: 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 capabili… ▽ More 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. △ Less

Submitted 2 April, 2020; originally announced April 2020.

Comments: 20 pages, 15 figures, accepted for publication in A&A

Magnetic Reconnection and Hot Spot Formation in Black Hole Accretion Disks

Authors: Bart Ripperda, Fabio Bacchini, Alexander Philippov

Abstract: Hot spots, or plasmoids, forming due to magnetic reconnection in current sheets, are conjectured to power frequent X-ray and near-infrared flares from Sgr A*, the black hole in the center of our Galaxy. It is unclear how, where, and when current sheets form in black-hole accretion disks. We perform axisymmetric general-relativistic resistive magnetohydrodynamics simulations to model reconnection a… ▽ More Hot spots, or plasmoids, forming due to magnetic reconnection in current sheets, are conjectured to power frequent X-ray and near-infrared flares from Sgr A*, the black hole in the center of our Galaxy. It is unclear how, where, and when current sheets form in black-hole accretion disks. We perform axisymmetric general-relativistic resistive magnetohydrodynamics simulations to model reconnection and plasmoid formation in a range of accretion flows. Current sheets and plasmoids are ubiquitous features which form regardless of the initial magnetic field in the disk, the magnetization in the quasi-steady-state phase of accretion, and the spin of the black hole. Within 10 Schwarzschild radii from the event horizon, we observe plasmoids forming, after which they can merge, grow to macroscopic scales of the order of a few Schwarzschild radii, and are ultimately advected along the jet's sheath or into the disk. Large plasmoids are energized to relativistic temperatures via reconnection and contribute to the jet's limb-brightening. We find that only hot spots forming in magnetically arrested disks can potentially explain the energetics of Sgr A* flares. The flare period is determined by the reconnection rate, which we find to be between 0.01c and 0.03c in all cases, consistent with studies of reconnection in isolated Harris-type current sheets. We quantify magnetic dissipation and non-ideal electric fields which can efficiently inject non-thermal particles. We show that explicit resistivity allows for converged numerical solutions, such that the electromagnetic energy evolution and dissipation become independent of the grid scale for the extreme resolutions considered here. △ Less

Submitted 7 September, 2020; v1 submitted 9 March, 2020; originally announced March 2020.

Comments: Added Sgr A* flare analysis, flare energetics and statistics, retrograde and Schwarzschild cases, new figure on convergence of adaptive mesh refinement

Journal ref: Astrophysical Journal, 900, 2 (2020)

Ideal MHD instabilities for coronal mass ejections

Authors: Rony Keppens, Yang Guo, Kirit Makwana, Zhixing Mei, Bart Ripperda, Chun Xia, Xiaozhou Zhao

Abstract: 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 toru… ▽ More 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. △ Less

Submitted 28 October, 2019; originally announced October 2019.

Comments: Reviews of Modern Plasma Physics (accepted), 14 figures, 29 pages

General relativistic resistive magnetohydrodynamics with robust primitive variable recovery for accretion disk simulations

Authors: Bart Ripperda, Fabio Bacchini, Oliver Porth, Elias R. Most, Hector Olivares, Antonios Nathanail, Luciano Rezzolla, Jannis Teunissen, Rony Keppens

Abstract: 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… ▽ More 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. △ Less

Submitted 6 August, 2019; v1 submitted 16 July, 2019; originally announced July 2019.

Comments: Matching accepted version in ApJS

Relativistic resistive magnetohydrodynamic reconnection and plasmoid formation in merging flux tubes

Authors: Bart Ripperda, Oliver Porth, Lorenzo Sironi, Rony Keppens

Abstract: 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 obtaine… ▽ More 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 $η= 10^{-2}$ down to $η= 10^{-4}$, for a reference case of magnetisation $σ= 3.33$ and plasma-$β= 0.1$. {For uniform resistivity $η=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 $η=10^{-4}$.} For non-uniform resistivity with a base level $η_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 $σ$ and the plasma-$β$ is analysed for cases with uniform resistivity $η=5\times10^{-5}$ and $η=10^{-4}$ in a range $0.5 \leq σ\leq 10$ and $0.01 \leq β\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$. △ Less

Submitted 6 February, 2019; v1 submitted 23 October, 2018; originally announced October 2018.

Comments: Matching accepted version in MNRAS

Particle orbits at the magnetopause: Kelvin-Helmholtz induced trapping

Authors: Matthieu Leroy, Bart Ripperda, Rony Keppens

Abstract: 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… ▽ More 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. △ Less

Submitted 6 August, 2019; v1 submitted 9 October, 2018; originally announced October 2018.

Comments: Matching accepted version in AGU JGR: Space Physics

Test particles in relativistic resistive magnetohydrodynamics

Authors: Bart Ripperda, Oliver Porth, Rony Keppens

Abstract: 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… ▽ More 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-$β$ 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. △ Less

Submitted 6 February, 2019; v1 submitted 9 October, 2018; originally announced October 2018.

Comments: Matching accepted version in J. Phys.: Conf. Ser. Astronum 2018 Proceedings

Numerical methods for General Relativistic particles

Authors: Fabio Bacchini, Bart Ripperda, Lorenzo Sironi

Abstract: 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)… ▽ More 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. △ Less

Submitted 3 October, 2018; v1 submitted 2 October, 2018; originally announced October 2018.

Comments: Proceedings of the IAU meeting (Noto, Italy, May 2018)

Journal ref: Proc. IAU 14 (2018) 19-23

Generalized, energy-conserving numerical simulations of particles in general relativity. II. Test particles in electromagnetic fields and GRMHD

Authors: Fabio Bacchini, Bart Ripperda, Oliver Porth, Lorenzo Sironi

Abstract: 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 presen… ▽ More 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). △ Less

Submitted 25 March, 2019; v1 submitted 1 October, 2018; originally announced October 2018.

Generalized, energy-conserving numerical simulations of particles in general relativity. I. Time-like and null geodesics

Authors: Fabio Bacchini, Bart Ripperda, Alexander Yuran Chen, Lorenzo Sironi

Abstract: 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. Schwarzs… ▽ More 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. △ Less

Submitted 31 May, 2018; v1 submitted 8 January, 2018; originally announced January 2018.

A comprehensive comparison of relativistic particle integrators

Authors: Bart Ripperda, Fabio Bacchini, Jannis Teunissen, Chun Xia, Oliver Porth, Lorenzo Sironi, Giovanni Lapenta, Rony Keppens

Abstract: 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 approximatio… ▽ More 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. △ Less

Submitted 6 February, 2019; v1 submitted 25 October, 2017; originally announced October 2017.

Comments: Matching published version in ApJS

Reconnection and particle acceleration in interacting flux ropes -- II. 3D effects on test particles in magnetically dominated plasmas

Authors: B. Ripperda, O. Porth, C. Xia, R. Keppens

Abstract: 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 curren… ▽ More 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. △ Less

Submitted 6 February, 2019; v1 submitted 27 July, 2017; originally announced July 2017.

Comments: Matches published version in MNRAS, including Erratum

Reconnection and particle acceleration in interacting flux ropes I. Magnetohydrodynamics and test particles in 2.5D

Authors: Bart Ripperda, Oliver Porth, Chun Xia, Rony Keppens

Abstract: 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 ma… ▽ More 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énic 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-$β$ 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. △ Less

Submitted 10 February, 2017; v1 submitted 29 November, 2016; originally announced November 2016.

Comments: Accepted in MNRAS