Search for: All records

Magnetic reconnection—a fundamental plasma physics process, where magnetic field lines of opposite polarity annihilate—is invoked in astrophysical plasmas as a powerful mechanism of nonthermal particle acceleration, able to explain fast-evolving, bright high-energy flares. Near black holes and neutron stars, reconnection occurs in the relativistic regime, in which the mean magnetic energy per particle exceeds the rest mass energy. This review reports recent advances in our understanding of the kinetic physics of relativistic reconnection:▪Kinetic simulations have elucidated the physics of plasma heating and nonthermal particle acceleration in relativistic reconnection (RR).▪The physics of radiative RR, with its self-consistent interplay between photons and reconnection-accelerated particles—a peculiarity of luminous, high-energy astrophysical sources—is the new frontier of research.▪RR plays a key role in global models of high-energy sources, in terms of both global-scale layers as well as reconnection sites generated as a by-product of local magnetohydrodynamic instabilities. We summarize themes of active investigation and future directions, emphasizing the role of upcoming observational capabilities, laboratory experiments, and new computational tools. 

Abstract A wide range of astrophysical sources exhibit extreme and rapidly varying electromagnetic emission indicative of efficient nonthermal particle acceleration. Understanding these sources often involves comparing data with a broad range of theoretical scenarios. To this end, it is beneficial to have tools that enable not only fast and efficient parametric investigation of the predictions of a specific scenario but also the flexibility to explore different theoretical ideas. In this paper, we introduceTleco, a versatile and lightweight toolkit for developing numerical models of relativistic outflows, including their particle acceleration mechanisms and resultant electromagnetic signature. Built on the Rust programming language and wrapped into a Python library,Tlecooffers efficient algorithms for evolving relativistic particle distributions and for solving the resulting emissions in a customizable fashion.Tlecouses a fully implicit discretization algorithm to solve the Fokker–Planck equation with user-defined diffusion, advection, cooling, injection, and escape and offers prescriptions for radiative emission and cooling. These include, but are not limited to, synchrotron, inverse-Compton, and self-synchrotron absorption.Tlecois designed to be user friendly and adaptable to model particle acceleration and the resulting electromagnetic spectrum and temporal variability in a wide variety of astrophysical scenarios, including, but not limited to, gamma-ray bursts, pulsar wind nebulae, and jets from active galactic nuclei. In this work, we outline the core algorithms and proceed to evaluate and demonstrate their effectiveness. The code is open source and available in the GitHub repository:https://github.com/zkdavis/Tleco. 

Abstract In 2023, the Pulsar Timing Array Collaborations announced the discovery of a gravitational wave background (GWB), predominantly attributed to supermassive black hole binary (SMBHB) mergers. However, the detected GWB is several times stronger than the default value expected from galactic observations at low and moderate redshifts. Recent findings by the James Webb Space Telescope have unveiled a substantial number of massive, high-redshift galaxies, suggesting more massive SMBHB mergers at these early epochs. Motivated by these findings, we propose an “early merger” model that complements the standard merger statistics by incorporating these early, massive galaxies. We compare the early and standard “late merger” models, which assume peak merger rates in the local Universe, and match both merger models to the currently detected GWB. Our analysis shows that the early merger model has a significantly lower detection probability for single binaries and predicts a ∼30% likelihood that the first detectable single source will be highly redshifted and remarkably massive with rapid frequency evolution. In contrast, the late merger model predicts a nearly monochromatic first source at low redshift. The future confirmation of an enhanced population of massive high-redshift galaxies and the detection of fast-evolving binaries would strongly support the early merger model, offering significant insights into the evolution of galaxies and SMBHs. 

ABSTRACT Magnetic reconnection is often considered as the primary particle acceleration mechanism in a magnetized blazar zone environment. The majority of radiation in the reconnection layer comes from plasmoids and their mergers. In particular, plasmoid mergers can produce strong multiwavelength flares and major variations in synchrotron polarization signatures. However, radiative properties of plasmoid mergers have not been well explored due to difficulties in tracking the merging processes. Here we use an image processing method that combines the magnetic vector potential and density to identify isolated and merging plasmoids. We find that this method can clearly distinguish radiation contributions from isolated plasmoids, merging plasmoids, and the primary current sheet of reconnection. This new method enables us to study the radiative properties of plasmoids and mergers statistically. Our results show that isolated plasmoids have similar emissivity regardless of their sizes, and they generally have non-zero polarization degree (PD) due to their quasi-circular shape. Flares due to plasmoid mergers have relative amplitudes that are antiproportional to the size ratio of the plasmoids participating in the mergers. Finally, only mergers between plasmoids of comparable sizes (width ratio ≲5) can lead to significant spectral hardening and polarization angle (PA) variations; the amplitude of the PA variations is between 0 and 180° and has a mean value of 90°. Our analyses on 2D simulations can pave the way for future analyses and machine learning techniques on radiative properties of 3D magnetic reconnection simulations. 

Abstract Gamma-ray burst (GRB) afterglows are emissions from ultrarelativistic blast waves produced by a narrow jet interacting with surrounding matter. Since the first multimessenger observation of a neutron star merger, hydrodynamic modeling of GRB afterglows for structured jets with smoothly varying angular energy distributions has gained increased interest. While the evolution of a jet is well described by self-similar solutions in both ultrarelativistic and Newtonian limits, modeling the transitional phase remains challenging. This is due to the nonlinear spreading of a narrow jet to a spherical configuration and the breakdown of self-similar solutions. Analytical models are limited in capturing these nonlinear effects, while relativistic hydrodynamic simulations are computationally expensive, which restricts the exploration of various initial conditions. In this work, we introduce a reduced hydrodynamic model that approximates the blast wave as an infinitely thin two-dimensional surface. Further assuming axial symmetry, this model simplifies the simulation to one dimension and drastically reduces the computational costs. We have compared our modeling to relativistic hydrodynamic simulations and semianalytic methods, and applied it to fit the light curve and flux centroid motion of GRB 170817A. These comparisons demonstrate good agreement and validate our approach. We have developed this method into a numerical tool,jetsimpy, which models the synchrotron GRB afterglow emission from a blast wave with arbitrary angular energy and Lorentz factor distribution. Although the code is built with GRB afterglow in mind, it applies to any relativistic jet. This tool is particularly useful in Markov Chain Monte Carlo studies and is provided to the community. 

Abstract Kinetic simulations of relativistic turbulence have significantly advanced our understanding of turbulent particle acceleration. Recent progress has highlighted the need for an updated acceleration theory that can account for particle acceleration within the plasma’s coherent structures. Here, we investigate how intermittency modeling connects statistical fluctuations in turbulence to regions of high-energy dissipation. This connection is established by employing a generalized She–Leveque model to characterize the exponentsζ_pfor the structure functions $S^p\propto l^{{\zeta }_p}$. The fitting of the scaling exponents provides us with a measure of the codimension of the dissipative structures, for which we subsequently determine the filling fraction. We perform our analysis for a range of magnetizationsσand relative fluctuation amplitudesδB₀/B₀. We find that increasing values ofσandδB₀/B₀allow the turbulent cascade to break sheetlike structures into smaller regions of dissipation that resemble chains of flux ropes. However, as their dissipation measure increases, the dissipative regions become less volume filling. With this work, we aim to inform future turbulent acceleration theories that incorporate particle energization from interactions with coherent structures within relativistic turbulence. 

ABSTRACT Kilonovae are optical transients following the merger of neutron star binaries, which are powered by the r-process heating of merger ejecta. However, if a merger remnant is a long-lived supramassive neutron star supported by its uniform rotation, it will inject energy into the ejecta through spin-down power. The energy injection can boost the peak luminosity of a kilonova by many orders of magnitudes, thus significantly increasing the detectable volume. Therefore, even if such events are only a small fraction of the kilonova population, they could dominate the detection rates. However, after many years of optical sky surveys, no such event has been confirmed. In this work, we build a boosted kilonova model with rich physical details, including the description of the evolution and stability of a proto neutron star, and the energy absorption through X-ray photoionization. We simulate the observation prospects and find the only way to match the absence of detection is to limit the energy injection by the newly born magnetar to only a small fraction of the neutron star rotational energy, thus they should collapse soon after the merger. Our result indicates that most supramassive neutron stars resulting from binary neutron star mergers are short lived and they are likely to be rare in the Universe. 

Abstract Magnetic reconnection is often invoked as a source of high-energy particles, and in relativistic astrophysical systems it is regarded as a prime candidate for powering fast and bright flares. We present a novel analytical model—supported and benchmarked with large-scale three-dimensional kinetic particle-in-cell simulations in electron–positron plasmas—that elucidates the physics governing the generation of power-law energy spectra in relativistic reconnection. Particles with Lorentz factorγ≳ 3σ(here,σis the magnetization) gain most of their energy in the inflow region, while meandering between the two sides of the reconnection layer. Their acceleration time is $t_{\mathrm{acc}}\sim \gamma {\eta }_{\mathrm{rec}}^{-1}{\omega }_{\mathrm{c}}^{-1}\simeq 20\gamma {\omega }_{\mathrm{c}}^{-1}$, whereη_rec≃ 0.06 is the inflow speed in units of the speed of light andω_c=eB₀/mcis the gyrofrequency in the upstream magnetic field. They leave the region of active energization aftert_esc, when they get captured by one of the outflowing flux ropes of reconnected plasma. We directly measuret_escin our simulations and find thatt_esc∼t_accforσ≳ few. This leads to a universal (i.e.,σ-independent) power-law spectrum ${\mathit{dN}}_{\mathrm{free}}/d\gamma \propto {\gamma }^{-1}$for the particles undergoing active acceleration, and $\mathit{dN}/d\gamma \propto {\gamma }^{-2}$for the overall particle population. Our results help to shed light on the ubiquitous presence of power-law particle and photon spectra in astrophysical nonthermal sources. 

ABSTRACT Despite a generally accepted framework for describing the gamma-ray burst (GRB) afterglows, the nature of the compact object at the central engine and the mechanism behind the prompt emission remain debated. The striped jet model is a promising venue to connect the various GRB stages since it gives a robust prediction for the relation of jet bulk acceleration, magnetization, and dissipation profile as a function of distance. Here, we use the constraints of the magnetization and bulk Lorentz of the jet flow at the large scales, where the jet starts interacting with the ambient gas in a large sample of bursts to (i) test the striped jet model for the GRB flow and (ii) study its predictions for the prompt emission and the constraints on the nature of the central engine. We find that the peak of the photospheric component of the emission predicted by the model is in agreement with the observed prompt emission spectra in the majority of the bursts in our sample, with a radiative efficiency of about 10 per cent. Furthermore, we adopt two different approaches to correlate the peak energies of the bursts with the type of central engine to find that more bursts are compatible with a neutron star central engine compared to a black hole one. Lastly, we conclude that the model favours broader distribution of stripe length-scales which results in a more gradual dissipation profile in comparison to the case, where the jet stripes are characterized by a single length-scale. 

ABSTRACT Blazars are a rare class of active galactic nuclei (AGNs) with relativistic jets pointing towards the observer. Jets are thought to be launched as Poynting-flux dominated outflows that accelerate to relativistic speeds at the expense of the available magnetic energy. In this work, we consider electron–proton jets and assume that particles are energized via magnetic reconnection in parts of the jet where the magnetization is still high (σ ≥ 1). The magnetization and bulk Lorentz factor Γ are related to the available jet energy per baryon as μ = Γ(1 + σ). We adopt an observationally motivated relation between Γ and the mass accretion rate into the black hole $$\dot{m}$$, which also controls the luminosity of external radiation fields. We numerically compute the photon and neutrino jet emission as a function of μ and σ. We find that the blazar SED is produced by synchrotron and inverse Compton radiation of accelerated electrons, while the emission of hadronic-related processes is subdominant except for the highest magnetization considered. We show that low-luminosity blazars (Lγ ≲ 1045 erg s−1) are associated with less powerful, slower jets with higher magnetizations in the jet dissipation region. Their broad-band photon spectra resemble those of BL Lac objects, and the expected neutrino luminosity is $$L_{\nu +\bar{\nu }}\sim (0.3-1)\, L_{\gamma }$$. High-luminosity blazars (Lγ ≫ 1045 erg s−1) are associated with more powerful, faster jets with lower magnetizations. Their broad-band photon spectra resemble those of flat spectrum radio quasars, and they are expected to be dim neutrino sources with $$L_{\nu +\bar{\nu }}\ll L_{\gamma }$$.