Search for: All records

Abstract Gamma-ray bursts (GRBs) are among the most energetic events in the Universe, driven by relativistic jets launched from black holes (BHs) formed during the collapse of massive stars or after the merger of two neutron stars. The jet power depends on the BH spin and the magnetic flux accreted onto it. In the standard thin disk model, jet power is limited by insufficient magnetic flux, even when the spin approaches maximum possible value. In contrast, the magnetically arrested disk (MAD) state limits jet energy by extracting significant angular momentum, braking BH rotation. We propose a unified model incorporating both standard thin disk and MAD states, identifying a universal curve for jet power per accretion rate as a function of the magnetic flux ratio, ${\Delta }_{\mathrm{eq}}={({\Phi }_{\mathrm{BH}}/{\Phi }_{\mathrm{MAD}})}_{\mathrm{eq}}$, at spin equilibrium. For long GRBs (lGRBs), the model predicts a maximum jet energy of ∼1.5% of the accretion energy, occurring at Δ_eq ∼ 0.4, where the BH equilibrium spin isa ∼ 0.5. Both long and short GRBs are unlikely to be produced by a MAD: for short GRBs, this requires an accreted mass orders of magnitude smaller than that available, while for lGRBs, the narrow progenitor mass distribution challenges the ability to produce the observed broad distribution of jet energies. This framework provides a consistent explanation for both standard and luminous GRBs, emphasizing the critical role of magnetic flux. Both long and short GRBs require magnetic flux distributions that peak around 10²⁷G cm². 

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 Recent advances in numerical simulations of magnetically arrested accretion onto supermassive black holes have shed light on the formation and dynamics of magnetospheric current sheets near the black hole horizon. By considering the pair magnetizationσ_ein the upstream region and the mass accretion rateṁ(in units of the Eddington mass accretion rate) as free parameters we estimate the strength of the magnetic field and develop analytical models, motivated by recent three-dimensional particle-in-cell simulations, to describe the populations of relativistic electrons and positrons (pairs) in the reconnection region.Applying our model to M87*, we numerically compute the non-thermal photon spectra for various values ofσ_e. We show that pairs that are accelerated up to the synchrotron radiation-limited energy while meandering across both sides of the current sheet, can produce MeV flares with luminosity of ∼ 10⁴¹ erg s^-1— independent ofσ_e— for a black hole accreting atṁ=10^-5. Pairs that are trapped in the transient current sheet can produce X-ray counterparts to the MeV flares, lasting about a day for current sheets with length of a few gravitational radii. We also show that the upstream plasma can be enriched due to photon-photon pair creation, and derive a new equilibrium magnetization ofσ_e∼ 10³-10⁴forṁ= 10^-6- 10^-5. Additionally, we explore the potential of magnetospheric current sheets to accelerate protons to ultra-high energies, finding that while acceleration to such energies is limited by various loss mechanisms, such as synchrotron and photopion losses from the non-thermal emission from pairs, maximal proton energies in the range of a few EeV are attainable in magnetospheric sheets forming around supermassive sub-Eddington accreting black holes. 

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 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 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 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. 

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.