Search for: All records

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

Abstract Models invoking magnetic reconnection as the particle acceleration mechanism within relativistic jets often adopt a gradual energy dissipation profile within the jet. However, such a profile has yet to be reproduced in first-principles simulations. Here we perform a suite of 3D general relativistic magnetohydrodynamic simulations of post–neutron star merger disks with an initially purely toroidal magnetic field. We explore the variations in both the microphysics (e.g., nuclear recombination, neutrino emission) and system parameters (e.g, disk mass). In all of our simulations, we find the formation of magnetically striped jets. The stripes result from the reversals in the poloidal magnetic flux polarity generated in the accretion disk. The simulations display large variations in the distributions of stripe duration,τ, and power, 〈P_Φ〉. We find that more massive disks produce more powerful stripes, the most powerful of which reaches 〈P_Φ〉 ∼ 10⁴⁹erg s⁻¹atτ∼ 20 ms. The power and variability that result from the magnetic reconnection of the stripes agree with those inferred in short-duration gamma-ray bursts. We find that the dissipation profile of the cumulative energy is roughly a power law in both radial distance,z, andτ, with a slope in the range of ∼1.7–3; more massive disks display larger slopes. 

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 Blazar emission is dominated by nonthermal radiation processes that are highly variable across the entire electromagnetic spectrum. Turbulence, which can be a major source of nonthermal particle acceleration, can widely exist in the blazar emission region. The Turbulent Extreme Multi-Zone (TEMZ) model has been used to describe turbulent radiation signatures. Recent particle-in-cell (PIC) simulations have also revealed the stochastic nature of the turbulent emission region and particle acceleration therein. However, radiation signatures have not been systematically studied via first-principles-integrated simulations. In this paper, we perform combined PIC and polarized radiative transfer simulations to study synchrotron emission from magnetic turbulence in the blazar emission region. We find that the multiwavelength flux and polarization are generally characterized by stochastic patterns. Specifically, the variability timescale and average polarization degree (PD) are governed by the correlation length of the turbulence. Interestingly, magnetic turbulence can result in polarization angle swings with arbitrary amplitudes and duration, in either direction, that are not associated with changes in flux or PD. Surprisingly, these swings, which are stochastic in nature, can appear either bumpy or smooth, although large-amplitude swings (>180°) are very rare, as expected. Our radiation and polarization signatures from first-principles-integrated simulations are consistent with the TEMZ model, except that in the latter, there is a weak correlation, with zero lag, between flux and degree of polarization. 

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