-
Tleco: A Toolkit for Modeling Radiative Signatures from Relativistic Outflows
Authors:
Zachary Davis,
Jesús M. Rueda-Becerril,
Dimitrios Giannios
Abstract:
A wide range of astrophysical sources exhibit extreme and rapidly varying electromagnetic emission indicative of efficient non-thermal 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 specif…
▽ More
A wide range of astrophysical sources exhibit extreme and rapidly varying electromagnetic emission indicative of efficient non-thermal 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 introduce \texttt{Tleco}, 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, \texttt{Tleco} offers efficient algorithms for evolving relativistic particle distributions and for solving the resulting emissions in a customizable fashion. \texttt{Tleco} uses a fully implicit discretization algorithm to solve the Fokker-Planck (FP) 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. \texttt{Tleco} is 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: \href{https://github.com/zkdavis/Tleco
△ Less
Submitted 27 May, 2024;
originally announced May 2024.
-
Balancing Turbulent Heating with Radiative Cooling in Blazars
Authors:
Zachary Davis,
Jesús M. Rueda-Becerril,
Dimitrios Giannios
Abstract:
Recently, particle in cell (PIC) simulations have shown that relativistic turbulence in collisionless plasmas can result in an equilibrium particle distribution function where turbulent heating is balanced by radiative cooling of electrons. Strongly magnetized plasmas are characterized by higher energy peaks and broader particle distributions. In relativistically moving astrophysical jets, it is b…
▽ More
Recently, particle in cell (PIC) simulations have shown that relativistic turbulence in collisionless plasmas can result in an equilibrium particle distribution function where turbulent heating is balanced by radiative cooling of electrons. Strongly magnetized plasmas are characterized by higher energy peaks and broader particle distributions. In relativistically moving astrophysical jets, it is believed that the flow is launched Poynting flux dominated and that the resulting magnetic instabilities may create a turbulent environment inside the jet, i.e., the regime of relativistic turbulence. In this paper, we extend previous PIC simulation results to larger values of plasma magnetization by linearly extrapolating the diffusion and advection coefficients relevant for the turbulent plasmas under consideration. We use these results to build a single zone turbulent jet model that is based on the global parameters of blazar emission region, and consistently calculate the particle distribution and resulting synchrotron and inverse Compton emission spectra. We then test our model by comparing its predictions with the broad-band quiescent emission spectra of a dozen blazars. Our results show good agreement with observations of low-synchrotron peaked (LSP) sources and find that LSPs are moderately Poynting flux dominated with magnetization $1\lesssim σ\lesssim 5$, have bulk Lorentz factor $Γ\sim 10-30$, and that the turbulent region is located at the edge, or just beyond, the broad line region (BLR). The turbulence is found to be driven at an area comparable to the jet cross section.
△ Less
Submitted 19 January, 2022;
originally announced January 2022.
-
Handing off the outcome of binary neutron star mergers for accurate and long-term post-merger simulations
Authors:
Federico G. Lopez Armengol,
Zachariah B. Etienne,
Scott C. Noble,
Bernard J. Kelly,
Leonardo R. Werneck,
Brendan Drachler,
Manuela Campanelli,
Federico Cipolletta,
Yosef Zlochower,
Ariadna Murguia-Berthier,
Lorenzo Ennoggi,
Mark Avara,
Riccardo Ciolfi,
Joshua Faber,
Grace Fiacco,
Bruno Giacomazzo,
Tanmayee Gupte,
Trung Ha,
Julian H. Krolik,
Vassilios Mewes,
Richard O'Shaughnessy,
Jesús M. Rueda-Becerril,
Jeremy Schnittman
Abstract:
We perform binary neutron star (BNS) merger simulations in full dynamical general relativity with IllinoisGRMHD, on a Cartesian grid with adaptive-mesh refinement. After the remnant black hole has become nearly stationary, the evolution of the surrounding accretion disk on Cartesian grids over long timescales (1s) is suboptimal, as Cartesian coordinates over-resolve the angular coordinates at larg…
▽ More
We perform binary neutron star (BNS) merger simulations in full dynamical general relativity with IllinoisGRMHD, on a Cartesian grid with adaptive-mesh refinement. After the remnant black hole has become nearly stationary, the evolution of the surrounding accretion disk on Cartesian grids over long timescales (1s) is suboptimal, as Cartesian coordinates over-resolve the angular coordinates at large distances, and the accreting plasma flows obliquely across coordinate lines dissipating angular momentum artificially from the disk. To address this, we present the Handoff, a set of computational tools that enables the transfer of general relativistic magnetohydrodynamic (GRMHD) and spacetime data from IllinoisGRMHD to HARM3D, a GRMHD code that specializes in modeling black hole accretion disks in static spacetimes over long timescales, making use of general coordinate systems with spherical topology. We demonstrate that the Handoff allows for a smooth and reliable transition of GRMHD fields and spacetime data, enabling us to efficiently and reliably evolve BNS dynamics well beyond merger. We also discuss future plans, which involve incorporating advanced equations of state and neutrino physics into BNS simulations using the \handoff approach.
△ Less
Submitted 31 October, 2022; v1 submitted 17 December, 2021;
originally announced December 2021.
-
HARM3D+NUC: A new method for simulating the post-merger phase of binary neutron star mergers with GRMHD, tabulated EOS and neutrino leakage
Authors:
Ariadna Murguia-Berthier,
Scott C. Noble,
Luke F. Roberts,
Enrico Ramirez-Ruiz,
Leonardo R. Werneck,
Michael Kolacki,
Zachariah B. Etienne,
Mark Avara,
Manuela Campanelli,
Riccardo Ciolfi,
Federico Cipolletta,
Brendan Drachler,
Lorenzo Ennoggi,
Joshua Faber,
Grace Fiacco,
Bruno Giacomazzo,
Tanmayee Gupte,
Trung Ha,
Bernard J. Kelly,
Julian H. Krolik,
Federico G. Lopez Armengol,
Ben Margalit,
Tim Moon,
Richard O'Shaughnessy,
Jesús M. Rueda-Becerril
, et al. (3 additional authors not shown)
Abstract:
The first binary neutron star merger has already been detected in gravitational waves. The signal was accompanied by an electromagnetic counterpart including a kilonova component powered by the decay of radioactive nuclei, as well as a short $γ$-ray burst. In order to understand the radioactively-powered signal, it is necessary to simulate the outflows and their nucleosynthesis from the post-merge…
▽ More
The first binary neutron star merger has already been detected in gravitational waves. The signal was accompanied by an electromagnetic counterpart including a kilonova component powered by the decay of radioactive nuclei, as well as a short $γ$-ray burst. In order to understand the radioactively-powered signal, it is necessary to simulate the outflows and their nucleosynthesis from the post-merger disk. Simulating the disk and predicting the composition of the outflows requires general relativistic magnetohydrodynamical (GRMHD) simulations that include a realistic, finite-temperature equation of state (EOS) and self-consistently calculating the impact of neutrinos. In this work, we detail the implementation of a finite-temperature EOS and the treatment of neutrinos in the GRMHD code HARM3D+NUC, based on HARM3D. We include formal tests of both the finite-temperature EOS and the neutrino leakage scheme. We further test the code by showing that, given conditions similar to those of published remnant disks following neutron star mergers, it reproduces both recombination of free nucleons to a neutron-rich composition and excitation of a thermal wind.
△ Less
Submitted 31 January, 2022; v1 submitted 9 June, 2021;
originally announced June 2021.
-
The blazar sequence revised
Authors:
Jesús M. Rueda-Becerril,
Amanda O. Harrison,
Dimitrios Giannios
Abstract:
We propose and test a fairly simple idea that could account for the blazar sequence: all jets are launched with similar energy per baryon, independently of their power. For instance, flat-spectrum radio quasars (FSRQs), the most powerful jets, manage to accelerate to high bulk Lorentz factor, as observed in the radio. As a result, the emission region will have a rather modest magnetization which w…
▽ More
We propose and test a fairly simple idea that could account for the blazar sequence: all jets are launched with similar energy per baryon, independently of their power. For instance, flat-spectrum radio quasars (FSRQs), the most powerful jets, manage to accelerate to high bulk Lorentz factor, as observed in the radio. As a result, the emission region will have a rather modest magnetization which will induce a steep particle spectra therein and a rather steep emission spectra in the gamma-rays; particularly in the \textit{Fermi}-LAT band. For the weaker jets, namely BL Lacertae objects (BL Lacs), the opposite holds true; i.e., the jet does not achieve a very high bulk Lorentz factor, leading to more magnetic energy available for non-thermal particle acceleration and harder emission spectra. Moreover, this model requires but a handful of parameters. By means of numerical simulations we have accomplished to reproduce the spectral energy distributions and light-curves from fiducial sources following the aforementioned model. With the a complete evolution of the broadband spectra we were able to study in detail the spectral features at any particular frequency band at any given stage. Finally numerical results are compared and contrasted with observations.
△ Less
Submitted 27 November, 2020;
originally announced November 2020.
-
A numerical approach for radiative cooling in relativistic outflows
Authors:
Jesús M. Rueda-Becerril
Abstract:
In high energy astrophysics scenarios such as blazars, GRBs or PWNe, it is highly probable that ultra-relativistic particles interact with photons in their environment through scattering. As long as the energy of the particle is greater than the energy of the interacting photon, the (classical) scattering is known to be in the Thomson regime. Otherwise, quantum effects will affect the scattering c…
▽ More
In high energy astrophysics scenarios such as blazars, GRBs or PWNe, it is highly probable that ultra-relativistic particles interact with photons in their environment through scattering. As long as the energy of the particle is greater than the energy of the interacting photon, the (classical) scattering is known to be in the Thomson regime. Otherwise, quantum effects will affect the scattering cross section, and we enter into the so-called Klein-Nishina regime. It is well known that radiative cooling in the Thomson regime is very efficient, leading to soft high-energy spectra. However, observations have shown that, in many cases, the high energy spectrum of some objects is rather hard. This has led to think that maybe particles are not being cooled down efficiently. Asymptotic approximations of the Klein-Nishina regime have been formulated in the last decades in order to account for these corrections in the distribution of particles responsible for the observed spectrum of high energy sources. In this work we presenta a numerical approach of the Klein-Nishina corrections to the radiative cooling. It has been developed to simulate the evolution of a distribution of particles interacting with photons in their surroundings via inverse Compton scattering.
△ Less
Submitted 27 November, 2020;
originally announced November 2020.
-
Blazar jets launched with similar energy per baryon, independently of their power
Authors:
Jesus M. Rueda-Becerril,
Amanda O. Harrison,
Dimitrios Giannios
Abstract:
The most extreme active galactic nuclei (AGN) are the radio active ones whose relativistic jet propagates close to our line of sight. These objects were first classified according to their emission line features into flat-spectrum radio quasars (FSRQs) and BL Lacertae objects (BL Lacs). More recently, observations revealed a trend between these objects known as the \emph{blazar sequence}, along wi…
▽ More
The most extreme active galactic nuclei (AGN) are the radio active ones whose relativistic jet propagates close to our line of sight. These objects were first classified according to their emission line features into flat-spectrum radio quasars (FSRQs) and BL Lacertae objects (BL Lacs). More recently, observations revealed a trend between these objects known as the \emph{blazar sequence}, along with an anti-correlation between the observed power and the frequency of the synchrotron peak. In the present work, we propose a fairly simple idea that could account for the whole blazar population: all jets are launched with similar energy per baryon, independently of their power. In the case of FSRQs, the most powerful jets, manage to accelerate to high bulk Lorentz factors, as observed in the radio. As a result, they have a rather modest magnetization in the emission region, resulting in magnetic reconnection injecting a steep particle energy distribution and, consequently, steep emission spectra in the $γ$-rays. For the weaker jets, namely BL Lacs, the opposite holds true; i.e., the jet does not achieve a very high bulk Lorentz factor, leading to more magnetic energy available for non-thermal particle acceleration, and harder emission spectra at frequencies $\gtrsim$ GeV. In this scenario, we recover all observable properties of blazars with our simulations, including the \emph{blazar sequence} for models with mild baryon loading ($50 \lesssim μ\lesssim 80$). This interpretation of the blazar population, therefore, tightly constrains the energy per baryon of blazar jets regardless of their accretion rate.
△ Less
Submitted 6 January, 2021; v1 submitted 4 September, 2020;
originally announced September 2020.
-
Inverse Compton Signatures of Gamma-Ray Burst Afterglows
Authors:
Hao Zhang,
Ian Christie,
Maria Petropoulou,
Jesus M. Rueda-Becerril,
Dimitrios Giannios
Abstract:
The afterglow emission from gamma-ray bursts (GRBs) is believed to originate from a relativistic blast wave driven into the circumburst medium. Although the afterglow emission from radio up to X-ray frequencies is thought to originate from synchrotron radiation emitted by relativistic, non-thermal electrons accelerated by the blast wave, the origin of the emission at high energies (HE; $\gtrsim$~G…
▽ More
The afterglow emission from gamma-ray bursts (GRBs) is believed to originate from a relativistic blast wave driven into the circumburst medium. Although the afterglow emission from radio up to X-ray frequencies is thought to originate from synchrotron radiation emitted by relativistic, non-thermal electrons accelerated by the blast wave, the origin of the emission at high energies (HE; $\gtrsim$~GeV) remains uncertain. The recent detection of sub-TeV emission from GRB~190114C by MAGIC raises further debate on what powers the very high-energy (VHE; $\gtrsim 300$GeV) emission. Here, we explore the inverse Compton scenario as a candidate for the HE and VHE emissions, considering two sources of seed photons for scattering: synchrotron photons from the blast wave (synchrotron self-Compton or SSC) and isotropic photon fields external to the blast wave (external Compton). For each case, we compute the multi-wavelength afterglow spectra and light curves. We find that SSC will dominate particle cooling and the GeV emission, unless a dense ambient infrared photon field, typical of star-forming regions, is present. Additionally, considering the extragalactic background light attenuation, we discuss the detectability of VHE afterglows by existing and future gamma-ray instruments for a wide range of model parameters. Studying GRB~190114C, we find that its afterglow emission in the \fermi-LAT band is synchrotron-dominated.The late-time \fermi-LAT measurement (i.e., $t\sim 10^4$~s), and the MAGIC observation also set an upper limit on the energy density of a putative external infrared photon field (i.e. $\lesssim 3\times 10^{-9}\,{\rm erg\,cm^{-3}}$), making the inverse Compton dominant in the sub-TeV energies.
△ Less
Submitted 17 June, 2020; v1 submitted 30 October, 2019;
originally announced October 2019.
-
On the influence of a Hybrid Thermal-Non thermal distribution in the Internal Shocks model for blazars
Authors:
Jesus M. Rueda-Becerril,
Petar Mimica,
Miguel-A. Aloy
Abstract:
Internal shocks occurring in blazars may accelerate both thermal and non-thermal electrons. In this paper we examine the consequences that such a hybrid (thermal/non-thermal) EED has on the spectrum of blazars. Since the thermal component of the EED may extend to very low energies. We replace the standard synchrotron process by the more general magneto-bremsstrahlung (MBS). Significant differences…
▽ More
Internal shocks occurring in blazars may accelerate both thermal and non-thermal electrons. In this paper we examine the consequences that such a hybrid (thermal/non-thermal) EED has on the spectrum of blazars. Since the thermal component of the EED may extend to very low energies. We replace the standard synchrotron process by the more general magneto-bremsstrahlung (MBS). Significant differences in the energy flux appear at low radio frequencies when considering MBS instead of the standard synchrotron emission. A drop in the spectrum appears in the all the radio band and a prominent valley between the infrared and soft X-rays bands when a hybrid EED is considered, instead of a power-law EED. In the $γ$-ray band an EED of mostly thermal particles displays significant differences with respect to the one dominated by non-thermal particles. A thermally-dominated EED produces a synchrotron self-Compton (SSC) peak extending only up to a few MeV, and the valley separating the MBS and the SSC peaks is much deeper than if the EED is dominated by non-thermal particles. The combination of these effects modifies the Compton dominance of a blazar, suggesting that the vertical scatter in the distribution of FSRQs and BL Lac objects in the peak synchrotron frequency - Compton dominance parameter space could be attributed to different proportions of thermal/non-thermal particles in the EED of blazars. Finally, the temperature of the electrons in the shocked plasma is shown to be a degenerated quantity for different magentizations of the ejected material.
△ Less
Submitted 21 February, 2017; v1 submitted 19 December, 2016;
originally announced December 2016.
-
Numerical simulations of the internal shock model in magnetized relativistic jets of blazars
Authors:
Jesus M. Rueda-Becerril,
Petar Mimica,
Miguel A. Aloy
Abstract:
The internal shocks scenario in relativistic jets is used to explain the variability of the blazar emission. Recent studies have shown that the magnetic field significantly alters the shell collision dynamics, producing a variety of spectral energy distributions and light-curves patterns. However, the role played by magnetization in such emission processes is still not entirely understood. In this…
▽ More
The internal shocks scenario in relativistic jets is used to explain the variability of the blazar emission. Recent studies have shown that the magnetic field significantly alters the shell collision dynamics, producing a variety of spectral energy distributions and light-curves patterns. However, the role played by magnetization in such emission processes is still not entirely understood. In this work we numerically solve the magnetohydodynamic evolution of the magnetized shells collision, and determine the influence of the magnetization on the observed radiation. Our procedure consists in systematically varying the shell Lorentz factor, relative velocity, and viewing angle. The calculations needed to produce the whole broadband spectral energy distributions and light-curves are computationally expensive, and are achieved using a high-performance parallel code.
△ Less
Submitted 27 February, 2015;
originally announced February 2015.
-
The influence of the magnetic field on the spectral properties of blazars
Authors:
J. M. Rueda-Becerril,
P. Mimica,
M. A. Aloy
Abstract:
We explore the signature imprinted by dynamically relevant magnetic fields on the spectral energy distribution (SED) of blazars. It is assumed that the emission from these sources originates from the collision of cold plasma shells, whose magnetohydrodynamic evolution we compute by numerically solving Riemann problems. We compute the SEDs including the most relevant radiative processes and scan a…
▽ More
We explore the signature imprinted by dynamically relevant magnetic fields on the spectral energy distribution (SED) of blazars. It is assumed that the emission from these sources originates from the collision of cold plasma shells, whose magnetohydrodynamic evolution we compute by numerically solving Riemann problems. We compute the SEDs including the most relevant radiative processes and scan a broad parameter space that encompasses a significant fraction of the commonly accepted values of not directly measurable physical properties. We reproduce the standard double hump SED found in blazar observations for unmagnetized shells, but show that the prototype double hump structure of blazars can also be reproduced if the dynamical source of the radiation field is very ultrarelativistic both, in a kinematically sense (namely, if it has Lorentz factors \gtrsim 50) and regarding its magnetization (e.g., with flow magnetizations σ\simeq 0.1). A fair fraction of the {\em blazar sequence} could be explained as a consequence of shell magnetization: negligible magnetization in FSRQs, and moderate or large (and uniform) magnetization in BL Lacs. The predicted photon spectral indices (\gph) in the γ-ray band are above the observed values (Γ_{\rm ph, obs} \lesssim 2.6 for sources with redshifts 0.4\le z \le 0.6) if the magnetization of the sources is moderate (σ\simeq 10^{-2}).
△ Less
Submitted 2 December, 2013; v1 submitted 21 October, 2013;
originally announced October 2013.
-
Numerical study of broadband spectra caused by internal shocks in magnetized relativistic jets of blazars
Authors:
Jesus M. Rueda-Becerril,
Petar Mimica,
Miguel A. Aloy,
Carmen Aloy
Abstract:
The internal-shocks scenario in relativistic jets has been used to explain the variability of blazars' outflow emission. Recent simulations have shown that the magnetic field alters the dynamics of these shocks producing a whole zoo of spectral energy density patterns. However, the role played by magnetization in such high-energy emission is still not entirely understood. With the aid of \emph{Fer…
▽ More
The internal-shocks scenario in relativistic jets has been used to explain the variability of blazars' outflow emission. Recent simulations have shown that the magnetic field alters the dynamics of these shocks producing a whole zoo of spectral energy density patterns. However, the role played by magnetization in such high-energy emission is still not entirely understood. With the aid of \emph{Fermi}'s second LAT AGN catalog, a comparison with observations in the $γ$-ray band was performed, in order to identify the effects of the magnetic field.
△ Less
Submitted 18 September, 2013;
originally announced September 2013.
-
Numerical simulations of dynamics and emission from relativistic astrophysical jets
Authors:
Petar Mimica,
Miguel Angel Aloy,
Jesus M. Rueda-Becerril,
Siham Tabik,
Carmen Aloy
Abstract:
Broadband emission from relativistic outflows (jets) of active galactic nuclei (AGN) and gamma-ray bursts (GRBs) contains valuable information about the nature of the jet itself, and about the central engine which launches it. Using special relativistic hydrodynamics and magnetohydronamics simulations we study the dynamics of the jet and its interaction with the surrounding medium. The observation…
▽ More
Broadband emission from relativistic outflows (jets) of active galactic nuclei (AGN) and gamma-ray bursts (GRBs) contains valuable information about the nature of the jet itself, and about the central engine which launches it. Using special relativistic hydrodynamics and magnetohydronamics simulations we study the dynamics of the jet and its interaction with the surrounding medium. The observational signature of the simulated jets is computed using a radiative transfer code developed specifically for the purpose of computing multi-wavelength, time-dependent, non-thermal emission from astrophysical plasmas. We present results of a series of long-term projects devoted to understanding the dynamics and emission of jets in parsec-scale AGN jets, blazars and the afterglow phase of the GRBs.
△ Less
Submitted 20 January, 2013; v1 submitted 8 November, 2012;
originally announced November 2012.
-
Spherical Boson Stars as Black Hole mimickers
Authors:
F. S. Guzman,
J. M. Rueda-Becerril
Abstract:
We present spherically symmetric boson stars as black hole mimickers based on the power spectrum of a simple accretion disk model. The free parameters of the boson star are the mass of the boson and the fourth order self-interaction coefficient in the scalar field potential. We show that even if the mass of the boson is the only free parameter it is possible to find a configuration that mimics the…
▽ More
We present spherically symmetric boson stars as black hole mimickers based on the power spectrum of a simple accretion disk model. The free parameters of the boson star are the mass of the boson and the fourth order self-interaction coefficient in the scalar field potential. We show that even if the mass of the boson is the only free parameter it is possible to find a configuration that mimics the power spectrum of the disk due to a black hole of the same mass. We also show that for each value of the self-interaction a single boson star configuration can mimic a black hole at very different astrophysical scales in terms of the mass of the object and the accretion rate. In order to show that it is possible to distinguish one of our mimickers from a black hole we also study the deflection of light.
△ Less
Submitted 7 September, 2010;
originally announced September 2010.
