-
Radiative Properties of Plasmoids and Plasmoid Mergers in Magnetic Reconnection
Authors:
Haocheng Zhang,
Lingyi Dong,
Dimitrios Giannios
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 multi-wavelength flares and major variations in synchrotron polarization signatures. However, radiative properties of plasmoid…
▽ More
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 multi-wavelength 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 nonzero polarization degree (PD) due to their quasi-circular shape. Flares due to plasmoid mergers have relative amplitudes that are anti-proportional to the size ratio of the plasmoids participating in the mergers. Finally, only mergers between plasmoids of comparable sizes (width ratio $\lesssim 5$) can lead to significant spectral hardening and polarization angle (PA) variations; the amplitude of the PA variations is between 0 and $180^{\circ}$ and has a mean value of $90^{\circ}$. Our analyses on 2D simulations can pave the way for future analyses and machine learning techniques on radiative properties of 3D magnetic reconnection simulations.
△ Less
Submitted 14 June, 2024;
originally announced June 2024.
-
The role of magnetospheric current sheets in pair enrichment and ultra-high energy proton acceleration in M87*
Authors:
S. I. Stathopoulos,
M. Petropoulou,
L. Sironi,
D. Giannios
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 $σ_{\rm e}$ in the upstream region and the mass accretion rate $\dot{m}$ (in units of the Eddington mass accretion rate) as free parameters we estim…
▽ More
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 $σ_{\rm e}$ in the upstream region and the mass accretion rate $\dot{m}$ (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 electron 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 $\sim 10^{41}$~erg s$^{-1}$ -- independent of $σ_e$ -- for a black hole accreting at $\dot{m}=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 \sim 10^3-10^4$ for $\dot{m}=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.
△ Less
Submitted 3 June, 2024;
originally announced June 2024.
-
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.
-
High-Redshift Merger Model for Low-Frequency Gravitational Wave Background
Authors:
Zhao-Feng Wu,
Dimitrios Giannios
Abstract:
In 2023, the Pulsar Timing Array (PTA) 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 (JWST) have…
▽ More
In 2023, the Pulsar Timing Array (PTA) 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 (JWST) 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.
△ Less
Submitted 3 August, 2024; v1 submitted 27 May, 2024;
originally announced May 2024.
-
jetsimpy: A Highly Efficient Hydrodynamic Code for Gamma-ray Burst Afterglow
Authors:
Hao Wang,
Ranadeep G. Dastidar,
Dimitrios Giannios,
Paul C. Duffell
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 describe…
▽ More
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, \texttt{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.
△ Less
Submitted 11 July, 2024; v1 submitted 29 February, 2024;
originally announced February 2024.
-
Intermittency and Dissipative Structures Arising from Relativistic Magnetized Turbulence
Authors:
Zachary Davis,
Luca Comisso,
Dimitrios Giannios
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 acceleration within the plasma's coherent structures. Here, we investigate how turbulent intermittency models connect statistical fluctuations in turbulence to regions of hig…
▽ More
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 acceleration within the plasma's coherent structures. Here, we investigate how turbulent intermittency models connect statistical fluctuations in turbulence to regions of high dissipation. This connection is established by employing a generalized She-Leveque model to describe the exponents $ζ_p$ for the structure functions $S^p \propto l^{ζ_p}$. The fitting of the scaling exponents provide us with a measure of the co-dimension of the dissipative structures, and we subsequently measure their filling fraction. We perform our analysis for a range of magnetizations $σ$ and magnetic field fluctuations ${δB_0}/{B_0}$. We find that increasing the values of $σ$ and ${δB_0}/{B_0}$ allows the cascade to break sheets into smaller regions of dissipation that resemble chains of plasmoids. However, as their dissipation 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.
△ Less
Submitted 15 March, 2024; v1 submitted 14 November, 2023;
originally announced November 2023.
-
Constraining the long-lived supramassive neutron stars by magnetar boosted kilonovae
Authors:
Hao Wang,
Paz Beniamini,
Dimitrios Giannios
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 spindown power. The energy injection can boost the peak luminosity of a kilonova by many orders of magnitud…
▽ More
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 spindown 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 kilonovae 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.
△ Less
Submitted 24 November, 2023; v1 submitted 17 August, 2023;
originally announced August 2023.
-
The origin of power-law spectra in relativistic magnetic reconnection
Authors:
Hao Zhang,
Lorenzo Sironi,
Dimitrios Giannios,
Maria Petropoulou
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 particle-in-cell simulations - that elucidates the physics governing the generation of power-law energy spect…
▽ More
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 particle-in-cell simulations - that elucidates the physics governing the generation of power-law energy spectra in relativistic reconnection. Particles with Lorentz factor $γ\gtrsim 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_{\rm acc}\sim γ\,η_{\rm rec}^{-1}ω_{\rm c}^{-1}\simeq 20\,γ\,ω_{\rm c}^{-1}$, where $η_{\rm rec}\simeq0.06$ is the inflow speed in units of the speed of light and $ω_{\rm c}=eB_0/mc$ is the gyrofrequency in the upstream magnetic field. They leave the region of active energization after $t_{\rm esc}$, when they get captured by one of the outflowing flux ropes of reconnected plasma. We directly measure $t_{\rm esc}$ in our simulations and find that $t_{\rm esc}\sim t_{\rm acc}$ for $σ\gtrsim {\rm few}$. This leads to a universal (i.e., $σ$-independent) power-law spectrum $dN_{\rm free}/dγ\propto γ^{-1}$ for the particles undergoing active acceleration, and $dN/dγ\propto γ^{-2}$ for the overall particle population. Our results help shedding light on the ubiquitous presence of power-law particle and photon spectra in astrophysical non-thermal sources.
△ Less
Submitted 23 February, 2023;
originally announced February 2023.
-
Scary Barbie: An Extremely Energetic, Long-Duration Tidal Disruption Event Candidate Without a Detected Host Galaxy at z = 0.995
Authors:
Bhagya M. Subrayan,
Dan Milisavljevic,
Ryan Chornock,
Raffaella Margutti,
Kate D. Alexander,
Vandana Ramakrishnan,
Paul C. Duffell,
Danielle A. Dickinson,
Kyoung-Soo Lee,
Dimitrios Giannios,
Geoffery Lentner,
Mark Linvill,
Braden Garretson,
Matthew J. Graham,
Daniel Stern,
Daniel Brethauer,
Tien Duong,
Wynn Jacobson-Galán,
Natalie LeBaron,
David Matthews,
Huei Sears,
Padma Venkatraman
Abstract:
We report multi-wavelength observations and characterization of the ultraluminous transient AT 2021lwx (ZTF20abrbeie; aka ``Barbie'') identified in the alert stream of the Zwicky Transient Facility (ZTF) using a Recommender Engine For Intelligent Transient Tracking (REFITT) filter on the ANTARES alert broker. From a spectroscopically measured redshift of 0.995, we estimate a peak observed pseudo-b…
▽ More
We report multi-wavelength observations and characterization of the ultraluminous transient AT 2021lwx (ZTF20abrbeie; aka ``Barbie'') identified in the alert stream of the Zwicky Transient Facility (ZTF) using a Recommender Engine For Intelligent Transient Tracking (REFITT) filter on the ANTARES alert broker. From a spectroscopically measured redshift of 0.995, we estimate a peak observed pseudo-bolometric luminosity of log (L$_{\text{max}} / [\text{erg}/\text{s}]$) = 45.7 from slowly fading ztf-$\it{g}$ and ztf-$r$ light curves spanning over 1000 observer-frame days. The host galaxy is not detected in archival Pan-STARRS observations ($g > 23.3$ mag), implying a lower limit to the outburst amplitude of more than 5 mag relative to the quiescent host galaxy. Optical spectra from Lick and Keck Observatories exhibit strong emission lines with narrow cores from the H Balmer series and ultraviolet semi-forbidden lines of Si III] $λ$1892, C III] $λ$1909, and C II] $λ$2325. Typical nebular lines in AGN spectra from ions such as [O II] and [O III] are not detected. These spectral features, along with the smooth light curve that is unlike most AGN flaring activity, and the luminosity that exceeds any observed or theorized supernova, lead us to conclude that AT 2021lwx is most likely an extreme tidal disruption event (TDE). Modeling of ZTF photometry with MOSFiT suggests that the TDE was between a $\approx 14 M_{\odot}$ star and a supermassive black hole of mass $M_{\text{BH}} \sim$ $10^{8} M_{\odot}$. Continued monitoring of the still-evolving light curve along with deep imaging of the field once AT 2021lwx has faded can test this hypothesis and potentially detect the host galaxy.
△ Less
Submitted 8 June, 2023; v1 submitted 21 February, 2023;
originally announced February 2023.
-
First-Principle-Integrated Study of Blazar Synchrotron Radiation and Polarization Signatures from Magnetic Turbulence
Authors:
Haocheng Zhang,
Alan Marscher,
Fan Guo,
Dimitrios Giannios,
Xiaocan Li,
Michela Negro
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 widely used to describe turbulent radiation signatures. Recent particle-in-cell (PIC) si…
▽ More
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 widely 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-principle-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 multi-wavelength flux and polarization are generally characterized by stochastic patterns. Specifically, the variability time scale and average polarization degree (PD) are governed by the correlation length of the turbulence. Interestingly, magnetic turbulence can result in polarization angle (PA) swings with arbitrary amplitudes and duration, in either directions, that are not associated to changes in flux or PD. Surprisingly, these swings, which are of stochastic nature, can appear either bumpy or smooth, although large amplitude swings ($>180^{\circ}$) are very rare as expected. Our radiation and polarization signatures from first-principle-integrated simulations are consistent with the TEMZ model.
△ Less
Submitted 20 March, 2023; v1 submitted 30 January, 2023;
originally announced January 2023.
-
Connecting the early afterglow to the prompt GRB and the central engine in the striped jet model
Authors:
Michail Damoulakis,
Rodolfo Barniol Duran,
Dimitrios Giannios
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…
▽ More
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 favors 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.
△ Less
Submitted 23 June, 2023; v1 submitted 21 December, 2022;
originally announced December 2022.
-
Striped Jets in Post Neutron Star Merger Systems
Authors:
Emma Kaufman,
I. M. Christie,
A. Lalakos,
A. Tchekhovskoy,
D. Giannios
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…
▽ More
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 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, $\langle P_Φ \rangle$. We find that more massive disks produce more powerful stripes, the most powerful of which reaches $\langle P_Φ \rangle \sim 10^{49}$~erg~s$^{-1}$ at $τ\sim 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 the slope in the range, $\sim 1.7-3$; more massive disks display larger slopes.
△ Less
Submitted 12 December, 2022;
originally announced December 2022.
-
Hadronic signatures from magnetically dominated baryon-loaded AGN jets
Authors:
Maria Petropoulou,
Filippos Psarras,
Dimitrios Giannios
Abstract:
Blazars are a rare class of active galactic nuclei (AGN) 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…
▽ More
Blazars are a rare class of active galactic nuclei (AGN) 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 ($σ\ge 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_γ \lesssim 10^{45}$ erg s$^{-1}$) are associated with less powerful, slower jets with higher magnetizations in the jet dissipation region. Their broadband photon spectra resemble those of BL Lac objects, and the expected neutrino luminosity is $L_{ν+\barν}\sim (0.3-1)\, L_γ$. High-luminosity blazars ($L_γ \gg 10^{45}$ erg s$^{-1}$) are associated with more powerful, faster jets with lower magnetizations. Their broadband photon spectra resemble those of flat spectrum radio quasars, and they are expected to be dim neutrino sources with $L_{ν+\barν}\ll L_γ$.
△ Less
Submitted 3 November, 2022;
originally announced November 2022.
-
Black hole to photosphere: 3D GRMHD simulations of collapsars reveal wobbling and hybrid composition jets
Authors:
Ore Gottlieb,
Matthew Liska,
Alexander Tchekhovskoy,
Omer Bromberg,
Aretaios Lalakos,
Dimitrios Giannios,
Philipp Mösta
Abstract:
Long-duration $γ$-ray bursts (GRBs) accompany the collapse of massive stars and carry information about the central engine. However, no 3D models have been able to follow these jets from their birth by a black-hole (BH) to the photosphere. We present the first such 3D general-relativity magnetohydrodynamic simulations, which span over 6 orders of magnitude in space and time. The collapsing stellar…
▽ More
Long-duration $γ$-ray bursts (GRBs) accompany the collapse of massive stars and carry information about the central engine. However, no 3D models have been able to follow these jets from their birth by a black-hole (BH) to the photosphere. We present the first such 3D general-relativity magnetohydrodynamic simulations, which span over 6 orders of magnitude in space and time. The collapsing stellar envelope forms an accretion disk, which drags inwardly the magnetic flux that accumulates around the BH, becomes dynamically important and launches bipolar jets. The jets reach the photosphere at $\sim10^{12}$ cm with an opening angle $θ_j\sim6^\circ$ and a Lorentz factor $Γ_j\lesssim 30$, unbinding $\gtrsim90\%$ of the star. We find that (i) the disk-jet system spontaneously develops misalignment relative to the BH rotational axis. As a result, the jet wobbles with an angle $θ_t\sim12^\circ$, which can naturally explain quiescent times in GRB lightcurves. The effective opening angle for detection $θ_j+θ_t$ suggests that the intrinsic GRB rate is lower by an order of magnitude than standard estimates. This suggests that successful GRBs may be rarer than currently thought and emerge in only $\sim 0.1\%$ of supernovae Ib/c, implying that jets are either not launched or choked inside most supernova Ib/c progenitors. (ii) The magnetic energy in the jet decreases due to mixing with the star, resulting in jets with a hybrid composition of magnetic and thermal components at the photosphere, where $\sim 10\%$ of the gas maintains magnetization $σ\gtrsim 0.1$. This indicates that both a photospheric component and reconnection may play a role in the prompt emission.
△ Less
Submitted 20 June, 2022; v1 submitted 26 April, 2022;
originally announced April 2022.
-
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.
-
Radiation and Polarization Signatures from Magnetic Reconnection in Relativistic Jets--II. Connection with $γ$-rays
Authors:
Haocheng Zhang,
Xiaocan Li,
Dimitrios Giannios,
Fan Guo,
Hannes Thiersen,
Markus Böttcher,
Tiffany Lewis,
Tonia Venters
Abstract:
It is commonly believed that blazar jets are relativistic magnetized plasma outflows from supermassive black holes. One key question is how the jets dissipate magnetic energy to accelerate particles and drive powerful multi-wavelength flares. Relativistic magnetic reconnection has been proposed as the primary plasma physical process in the blazar emission region. Recent numerical simulations have…
▽ More
It is commonly believed that blazar jets are relativistic magnetized plasma outflows from supermassive black holes. One key question is how the jets dissipate magnetic energy to accelerate particles and drive powerful multi-wavelength flares. Relativistic magnetic reconnection has been proposed as the primary plasma physical process in the blazar emission region. Recent numerical simulations have shown strong acceleration of nonthermal particles that may lead to multi-wavelength flares. Nevertheless, previous works have not directly evaluated $γ$-ray signatures from first-principle simulations. In this paper, we employ combined particle-in-cell and polarized radiation transfer simulations to study multi-wavelength radiation and optical polarization signatures under the leptonic scenario from relativistic magnetic reconnection. We find harder-when-brighter trends in optical and {\it Fermi-LAT} $γ$-ray bands as well as closely correlated optical and $γ$-ray flares. The optical polarization angle swings are also accompanied by $γ$-ray flares with trivial time delays. Intriguingly, we find highly variable synchrotron self Compton signatures due to inhomogeneous particle distributions during plasmoid mergers. This feature may result in fast $γ$-ray flares or orphan $γ$-ray flares under the leptonic scenario, complementary to the frequently considered mini-jet scenario. It may also infer neutrino emission with low secondary synchrotron flux under the hadronic scenario, if plasmoid mergers can accelerate protons to very high energy.
△ Less
Submitted 3 November, 2021;
originally announced November 2021.
-
Fast particle acceleration in three-dimensional relativistic reconnection
Authors:
Hao Zhang,
Lorenzo Sironi,
Dimitrios Giannios
Abstract:
Magnetic reconnection is invoked as one of the primary mechanisms to produce energetic particles. We employ large-scale three-dimensional (3D) particle-in-cell simulations of reconnection in magnetically-dominated ($σ=10$) pair plasmas to study the energization physics of high-energy particles. We identify a novel acceleration mechanism that only operates in 3D. For weak guide fields, 3D plasmoids…
▽ More
Magnetic reconnection is invoked as one of the primary mechanisms to produce energetic particles. We employ large-scale three-dimensional (3D) particle-in-cell simulations of reconnection in magnetically-dominated ($σ=10$) pair plasmas to study the energization physics of high-energy particles. We identify a novel acceleration mechanism that only operates in 3D. For weak guide fields, 3D plasmoids / flux ropes extend along the $z$ direction of the electric current for a length comparable to their cross-sectional radius. Unlike in 2D simulations, where particles are buried in plasmoids, in 3D we find that a fraction of particles with $γ\gtrsim 3σ$ can escape from plasmoids by moving along $z$, and so they can experience the large-scale fields in the upstream region. These "free" particles preferentially move in $z$ along Speiser-like orbits sampling both sides of the layer, and are accelerated linearly in time -- their Lorentz factor scales as $γ\propto t$, in contrast to $γ\propto \sqrt{t}$ in 2D. The energy gain rate approaches $\sim eE_{\rm rec}c$, where $E_{\rm rec}\simeq 0.1 B_0$ is the reconnection electric field and $B_0$ the upstream magnetic field. The spectrum of free particles is hard, $dN_{\rm free}/dγ\propto γ^{-1.5}$, contains $\sim 20\%$ of the dissipated magnetic energy independently of domain size, and extends up to a cutoff energy scaling linearly with box size. Our results demonstrate that relativistic reconnection in GRB and AGN jets may be a promising mechanism for generating ultra-high-energy cosmic rays.
△ Less
Submitted 30 April, 2021;
originally announced May 2021.
-
Evidence for X-ray Emission in Excess to the Jet Afterglow Decay 3.5 yrs After the Binary Neutron Star Merger GW 170817: A New Emission Component
Authors:
A. Hajela,
R. Margutti,
J. S. Bright,
K. D. Alexander,
B. D. Metzger,
V. Nedora,
A. Kathirgamaraju,
B. Margalit,
D. Radice,
C. Guidorzi,
E. Berger,
A. MacFadyen,
D. Giannios,
R. Chornock,
I. Heywood,
L. Sironi,
O. Gottlieb,
D. Coppejans,
T. Laskar,
Y. Cendes,
R. Barniol Duran,
T. Eftekhari,
W. Fong,
A. McDowell,
M. Nicholl
, et al. (12 additional authors not shown)
Abstract:
For the first $\sim3$ years after the binary neutron star merger event GW 170817 the radio and X-ray radiation has been dominated by emission from a structured relativistic off-axis jet propagating into a low-density medium with n $< 0.01\,\rm{cm^{-3}}$. We report on observational evidence for an excess of X-ray emission at $δt>900$ days after the merger. With…
▽ More
For the first $\sim3$ years after the binary neutron star merger event GW 170817 the radio and X-ray radiation has been dominated by emission from a structured relativistic off-axis jet propagating into a low-density medium with n $< 0.01\,\rm{cm^{-3}}$. We report on observational evidence for an excess of X-ray emission at $δt>900$ days after the merger. With $L_x\approx5\times 10^{38}\,\rm{erg\,s^{-1}}$ at 1234 days, the recently detected X-ray emission represents a $\ge 3.2\,σ$ (Gaussian equivalent) deviation from the universal post jet-break model that best fits the multi-wavelength afterglow at earlier times. In the context of JetFit afterglow models, current data represent a departure with statistical significance $\ge 3.1\,σ$, depending on the fireball collimation, with the most realistic models showing excesses at the level of $\ge 3.7\,σ$. A lack of detectable 3 GHz radio emission suggests a harder broad-band spectrum than the jet afterglow. These properties are consistent with the emergence of a new emission component such as synchrotron radiation from a mildly relativistic shock generated by the expanding merger ejecta, i.e. a kilonova afterglow. In this context, we present a set of ab-initio numerical-relativity BNS merger simulations that show that an X-ray excess supports the presence of a high-velocity tail in the merger ejecta, and argues against the prompt collapse of the merger remnant into a black hole. Radiation from accretion processes on the compact-object remnant represents a viable alternative. Neither a kilonova afterglow nor accretion-powered emission have been observed before, as detections of BNS mergers at this phase of evolution are unprecedented.
△ Less
Submitted 5 March, 2022; v1 submitted 5 April, 2021;
originally announced April 2021.
-
First-Principle Prediction of X-ray Polarization from Magnetic Reconnection in High-Frequency BL Lacs
Authors:
Haocheng Zhang,
Xiaocan Li,
Dimitrios Giannios,
Fan Guo
Abstract:
Relativistic magnetic reconnection is a potential particle acceleration mechanism for high-frequency BL Lacs (HBLs). The {\it Imaging X-ray Polarimetry Explorer} ({\it IXPE}) scheduled to launch in 2021 has the capability to probe the magnetic field evolution in HBLs, examining the magnetic reconnection scenario for the HBL flares. In this paper, we make the first attempt to self-consistently pred…
▽ More
Relativistic magnetic reconnection is a potential particle acceleration mechanism for high-frequency BL Lacs (HBLs). The {\it Imaging X-ray Polarimetry Explorer} ({\it IXPE}) scheduled to launch in 2021 has the capability to probe the magnetic field evolution in HBLs, examining the magnetic reconnection scenario for the HBL flares. In this paper, we make the first attempt to self-consistently predict HBL X-ray polarization signatures arising from relativistic magnetic reconnection via combined particle-in-cell (PIC) and polarized radiation transfer simulations. We find that although the intrinsic optical and X-ray polarization degrees are similar on average, the X-ray polarization is much more variable in both polarization degree and angle (PD and PA). Given the sensitivity of the {\it IXPE}, it may obtain one to a few polarization data points for one flaring event of nearby bright HBLs Mrk~421 and 501. However, it may not fully resolve the highly variable X-ray polarization. Due to the temporal depolarization, where the integration of photons with variable polarization states over a finite period of time can lower the detected PD, the measured X-ray PD can be considerably lower than the optical counterpart or even undetectable. The lower X-ray PD than the optical thus can be a characteristic signature of relativistic magnetic reconnection. For very bright flares where the X-ray polarization is well resolved, relativistic magnetic reconnection predicts smooth X-ray PA swings, which originate from large plasmoid mergers in the reconnection region.
△ Less
Submitted 27 March, 2021;
originally announced March 2021.
-
Radiation Signatures From Striped Blazar Jet
Authors:
Haocheng Zhang,
Dimitrios Giannios
Abstract:
Relativistic jets from supermassive black holes are among the most powerful and luminous astrophysical systems in Universe. We propose that the open magnetic field lines through the black hole, which drive a strongly magnetized jet, may have their polarity reversing over time scales related to the growth of the magneto-rotational dynamo in the disc, resulting in dissipative structures in the jet c…
▽ More
Relativistic jets from supermassive black holes are among the most powerful and luminous astrophysical systems in Universe. We propose that the open magnetic field lines through the black hole, which drive a strongly magnetized jet, may have their polarity reversing over time scales related to the growth of the magneto-rotational dynamo in the disc, resulting in dissipative structures in the jet characterized by reversing toroidal field polarities, referred to as ``stripes''. Magnetic reconnection between the stripes dissipates the magnetic energy and powers jet acceleration. The striped jet model can explain the jet acceleration, large-scale jet emission, and blazar emission signatures consistently in a unified physical picture. Specifically, we find that the jet accelerates to the bulk Lorentz factor $Γ\gtrsim 10$ within one parsec distance from the central engine. The acceleration slows down but continues at larger distances, with intrinsic acceleration rate $\dotΓ/Γ$ between $0.0005~\rm{yr^{-1}}$ and $0.005~\rm{yr^{-1}}$ at tens of parsecs, which is in very good agreement with recent radio observations. Magnetic reconnection continuously accelerates nonthermal particles over large distances from the central engine, resulting in the core-shift effect and overall flat-to-inverted synchrotron spectrum. The large-scale spectral luminosity peak $ν_{\rm peak}$ is anti-proportional to the location of the peak of the dissipation, which is set by the minimal stripe width $l_{\rm min}$. The blazar zone is approximately at the same location. At this distance, the jet is moderately magnetized, with the comoving magnetic field strength and dissipation power consistent with typical leptonic blazar model parameters.
△ Less
Submitted 30 December, 2020;
originally announced December 2020.
-
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.
-
Shock-powered radio precursors of neutron star mergers from accelerating relativistic binary winds
Authors:
Navin Sridhar,
Jonathan Zrake,
Brian D. Metzger,
Lorenzo Sironi,
Dimitrios Giannios
Abstract:
During the final stages of a compact object merger, if at least one of the binary components is a magnetized neutron star (NS), then its orbital motion substantially expands the NS's open magnetic flux -- and hence increases its wind luminosity -- relative to that of an isolated pulsar. As the binary orbit shrinks due to gravitational radiation, the power and speed of this binary-induced inspiral…
▽ More
During the final stages of a compact object merger, if at least one of the binary components is a magnetized neutron star (NS), then its orbital motion substantially expands the NS's open magnetic flux -- and hence increases its wind luminosity -- relative to that of an isolated pulsar. As the binary orbit shrinks due to gravitational radiation, the power and speed of this binary-induced inspiral wind may (depending on pair loading) secularly increase, leading to self-interaction and internal shocks in the outflow beyond the binary orbit. The magnetized forward shock can generate coherent radio emission via the synchrotron maser process, resulting in an observable radio precursor to binary NS merger. We perform 1D relativistic hydrodynamical simulations of shock interaction in the accelerating binary NS wind, assuming that the inspiral wind efficiently converts its Poynting flux into bulk kinetic energy prior to the shock radius. This is combined with the shock maser spectrum from particle-in-cell simulations, to generate synthetic radio light curves. The precursor burst with a fluence of $\sim1$ Jy$\cdot$ms at $\sim$GHz frequencies lasts $\sim 1-500$ ms following the merger for a source at $\sim3$ Gpc ($B_{\rm d}/10^{12}$ G)$^{8/9}$, where $B_{\rm d}$ is the dipole field strength of the more strongly-magnetized star. Given an outflow geometry concentrated along the binary equatorial, the signal may be preferentially observable for high-inclination systems, i.e. those least likely to produce a detectable gamma-ray burst.
△ Less
Submitted 18 January, 2021; v1 submitted 19 October, 2020;
originally announced October 2020.
-
Multimessenger parameter estimation of GW170817: from jet structure to the Hubble constant
Authors:
Hao Wang,
Dimitrios Giannios
Abstract:
The electromagnetic radiation that followed the neutron star merger event GW170817 revealed that gamma-ray burst afterglows from jets misaligned with our line of sight exhibit a light curve with slowly rising flux. The slope of the rising light curve depends sensitively on the angle of the observer with respect to the jet axis, which is likely to be perpendicular to the merger plane of the neutron…
▽ More
The electromagnetic radiation that followed the neutron star merger event GW170817 revealed that gamma-ray burst afterglows from jets misaligned with our line of sight exhibit a light curve with slowly rising flux. The slope of the rising light curve depends sensitively on the angle of the observer with respect to the jet axis, which is likely to be perpendicular to the merger plane of the neutron star binary. Therefore, the afterglow emission can be used to constrain the inclination of the merging system. Here, we calculate the gamma-ray burst afterglow emission based on the realistic jet structure derived from general-relativistic magnetohydrodynamical simulations of a black hole torus system for the central engine of the gamma-ray burst. Combined with gravitational wave parameter estimation, we fit the multi-epoch afterglow emission of GW170817. We show that with such a jet model, the observing angle can be tightly constrained by multi messenger observations. The best fit observing angle of GW170817 is $θ_{\rm v} = 0.38\pm 0.02$ rad. With such a constraint, we can break the degeneracy between inclination angle and luminosity distance in gravitational wave parameter estimation, and substantially increase the precision with which the Hubble constant is constrained by the standard siren method. Our estimation of the distance is $D_{\rm L}=43.4\pm 1\ \rm Mpc$ and the Hubble constant constraint is $69.5\pm 4\ \mathrm{km\ s^{-1}\ Mpc^{-1}}$. As a result, multimessenger observations of short-duration gamma-ray bursts, combined with a good theoretical understanding of the jet structure, can be powerful probes of cosmological parameters.
△ Less
Submitted 31 January, 2021; v1 submitted 9 September, 2020;
originally announced September 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.
-
Radiation and Polarization Signatures from Magnetic Reconnection in Relativistic Jets--I. A Systematic Study
Authors:
Haocheng Zhang,
Xiaocan Li,
Dimitrios Giannios,
Fan Guo,
Yi-Hsin Liu,
Lingyi Dong
Abstract:
Blazars are relativistic magnetized plasma outflows from supermassive black holes that point very close to our line of sight. Their emission is nonthermal dominated and highly variable across the entire electromagnetic spectrum. Relativistic magnetic reconnection has been proposed as the driver of particle acceleration during blazar flares. While recent particle-in-cell simulations have self-consi…
▽ More
Blazars are relativistic magnetized plasma outflows from supermassive black holes that point very close to our line of sight. Their emission is nonthermal dominated and highly variable across the entire electromagnetic spectrum. Relativistic magnetic reconnection has been proposed as the driver of particle acceleration during blazar flares. While recent particle-in-cell simulations have self-consistently studied the evolution of magnetic reconnection and particle acceleration therein, the resulting radiation signatures have not been systematically explored. In particular, the polarization signatures, which directly reflect the characteristic strongly dynamical magnetic field evolution during reconnection, have not been carefully investigated. In this paper, we present a systematic study of radiation and polarization signatures arising from magnetic reconnection in blazars, based on combined PIC and polarized radiation transfer simulations with various physical parameters. We identify a harder-when-brighter trend in the spectral evolution. Moreover, higher-frequency bands tend to flare earlier than lower-frequency bands in the synchrotron spectral component. Most importantly, polarization signatures appear more variable with higher frequencies. We find that the temporal polarization variations strongly depends on the guide field strength. Specifically, reconnection with significant guide field component leads to very high polarization degree that contradict to typical blazar observations, while large polarization angle rotations are unique signatures of magnetic reconnection between nearly anti-parallel magnetic field lines. These rotations are at least $90^o$ and can extend to $>180^o$, and they may rotate in both directions. These results imply that blazars that have shown large polarization angle rotations intrinsically have more nearly anti-parallel magnetic field morphology.
△ Less
Submitted 21 August, 2020;
originally announced August 2020.
-
Deciphering the properties of the central engine in GRB collapsars
Authors:
M. Petropoulou,
P. Beniamini,
G. Vasilopoulos,
D. Giannios,
R. Barniol Duran
Abstract:
The central engine in long gamma-ray bursts (GRBs) is thought to be a compact object produced by the core collapse of massive stars, but its exact nature (black hole or millisecond magnetar) is still debatable. Although the central engine of GRB collapsars is hidden to direct observation, its properties may be imprinted on the accompanying electromagnetic signals. We aim to decipher the generic pr…
▽ More
The central engine in long gamma-ray bursts (GRBs) is thought to be a compact object produced by the core collapse of massive stars, but its exact nature (black hole or millisecond magnetar) is still debatable. Although the central engine of GRB collapsars is hidden to direct observation, its properties may be imprinted on the accompanying electromagnetic signals. We aim to decipher the generic properties of central engines that are consistent with prompt observations of long GRBs detected by the Burst Alert Telescope (BAT) on board the Neil Gehrels Swift Observatory. Adopting a generic model for the central engine, in which the engine power and activity timescale are independent of each other, we perform Monte Carlo simulations of long GRBs produced by jets that successfully breakout from the star. Our simulations consider the dependence of the jet breakout timescale on the engine luminosity and the effects of the detector's flux threshold. The two-dimensional (2D) distribution of simulated detectable bursts in the gamma-ray luminosity versus gamma-ray duration plane is consistent with the observed one for a range of parameter values describing the central engine. The intrinsic 2D distribution of simulated collapsar GRBs peaks at lower gamma-ray luminosities and longer durations than the observed one, a prediction that can be tested in the future with more sensitive detectors. Black-hole accretors, whose power and activity time are set by the large-scale magnetic flux through the progenitor star and stellar structure, respectively, are compatible with the properties of the central engine inferred by our model.
△ Less
Submitted 12 June, 2020;
originally announced June 2020.
-
Kink Instabilities In Relativistic Jets Can Drive Quasi-Periodic Radiation Signatures
Authors:
Lingyi Dong,
Haocheng Zhang,
Dimitrios Giannios
Abstract:
Relativistic jets are highly collimated plasma outflows emerging from accreting black holes. They are launched with a significant amount of magnetic energy, which can be dissipated to accelerate nonthermal particles and give rise to electromagnetic radiation at larger scales. Kink instabilities can be an efficient mechanism to trigger dissipation of jet magnetic energy. While previous works have s…
▽ More
Relativistic jets are highly collimated plasma outflows emerging from accreting black holes. They are launched with a significant amount of magnetic energy, which can be dissipated to accelerate nonthermal particles and give rise to electromagnetic radiation at larger scales. Kink instabilities can be an efficient mechanism to trigger dissipation of jet magnetic energy. While previous works have studied the conditions required for the growth of kink instabilities in relativistic jets, the radiation signatures of these instabilities have not been investigated in detail. In this paper, we aim to self-consistently study radiation and polarization signatures from kink instabilities in relativistic jets. We combine large-scale relativistic magnetohydrodynamic (RMHD) simulations with polarized radiation transfer of a magnetized jet, which emerges from the central engine and propagates through the surrounding medium. We observe that a localized region at the central spine of the jet exhibits the strongest kink instabilities, which we identify as the jet emission region. Very interestingly, we find quasi-periodic oscillation (QPO) signatures in the light curve from the emission region. Additionally, the polarization degree appears to be anti-correlated to flares in the light curves. Our analyses show that these QPO signatures are intrinsically driven by kink instabilities, where the period of the QPOs is associated to the kink growth time scale. The latter corresponds to weeks to months QPOs in blazars. The polarization signatures offer unique diagnostics for QPOs driven by kink instabilities.
△ Less
Submitted 17 March, 2020;
originally announced March 2020.
-
FRB-periodicity: mild pulsars in tight O/B-star binaries
Authors:
Maxim Lyutikov,
Maxim Barkov,
Dimitrios Giannios
Abstract:
Periodicities observed in two Fast Radio Burst (FRB) sources (16 days in FRB 180916.J0158+65 and 160 days in FRB 121102) are consistent with that of tight, stellar mass binary systems. In the case of FRB 180916.J0158+65 the primary is an early OB-type star with mass loss rate $\dot{M} \sim 10^{-8}- 10^{-7} M_\odot$ yr$^{-1}$, and the secondary a neutron star. The observed periodicity is not intrin…
▽ More
Periodicities observed in two Fast Radio Burst (FRB) sources (16 days in FRB 180916.J0158+65 and 160 days in FRB 121102) are consistent with that of tight, stellar mass binary systems. In the case of FRB 180916.J0158+65 the primary is an early OB-type star with mass loss rate $\dot{M} \sim 10^{-8}- 10^{-7} M_\odot$ yr$^{-1}$, and the secondary a neutron star. The observed periodicity is not intrinsic to the FRB's source, but is due to the orbital phase-dependent modulation of the absorption conditions in the massive star's wind. The observed relatively narrow FRB activity window implies that the primary's wind dynamically dominates that of the pulsar, $η= L_{sd}/(\dot{M} v_w c) \leq 1$, where $L_{sd} $ is pulsar spin-down, $\dot{M}$ is the primary's wind mass loss rate and $v_w$ is its velocity. The condition $η\leq 1$ requires mildly powerful pulsar with $L_{sd} \lesssim 10^{37}$ erg $s^{-1}$. The observations are consistent with magnetically-powered radio emission originating in the magnetospheres of strongly magnetized neutron stars, the classical magnetars.
△ Less
Submitted 7 August, 2022; v1 submitted 5 February, 2020;
originally announced February 2020.
-
Ready, set, launch: time interval between BNS merger and short GRB jet formation
Authors:
Paz Beniamini,
Rodolfo Barniol Duran,
Maria Petropoulou,
Dimitrios Giannios
Abstract:
The joint detection of GW~170817/GRB 170817 confirmed the long-standing theory that binary neutron star mergers produce short gamma-ray burst (sGRB) jets that can successfully break out of the surrounding ejecta. At the same time, the association with a kilonova provided unprecedented information regarding the physical properties (such as masses and velocities) of the different ejecta constituents…
▽ More
The joint detection of GW~170817/GRB 170817 confirmed the long-standing theory that binary neutron star mergers produce short gamma-ray burst (sGRB) jets that can successfully break out of the surrounding ejecta. At the same time, the association with a kilonova provided unprecedented information regarding the physical properties (such as masses and velocities) of the different ejecta constituents. Combining this knowledge with the observed luminosities and durations of cosmological sGRBs detected by the Burst Alert Telescope (BAT) onboard the Neil Gehrels Swift Observatory, we revisit the breakout conditions of sGRB jets. Assuming self-collimation of sGRB jets does not play a critical role, we find that the time interval between the binary merger and the launching of a typical sGRB jet is $\lesssim0.1$~s. We also show that for a fraction of at least $\sim 30\%$ of sGRBs, the usually adopted assumption of static ejecta is inconsistent with observations, even if the polar ejecta mass is an order of magnitude smaller than the one in GRB 170817. Our results disfavour magnetar central engines for powering cosmological sGRBs, limit the amount of energy deposited in the cocoon prior to breakout, and suggest that the observed delay of $\sim 1.$7~s in GW 170817 /GRB 170817 between the gravitational wave and $γ$-ray signals is likely dominated by the propagation time of the jet to the $γ$-ray production site.
△ Less
Submitted 18 May, 2020; v1 submitted 3 January, 2020;
originally announced January 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.
-
Two years of non-thermal emission from the binary neutron star merger GW170817: rapid fading of the jet afterglow and first constraints on the kilonova fastest ejecta
Authors:
A. Hajela,
R. Margutti,
K. D. Alexander,
A. Kathirgamaraju,
A. Baldeschi,
C. Guidorzi,
D. Giannios,
W. Fong,
Y. Wu,
A. MacFadyen,
A. Paggi,
E. Berger,
P. K. Blanchard,
R. Chornock,
D. L. Coppejans,
P. S. Cowperthwaite,
T. Eftekhari,
S. Gomez,
G. Hosseinzadeh,
T. Laskar,
B. D. Metzger,
M. Nicholl,
K. Paterson,
D. Radice,
L. Sironi
, et al. (5 additional authors not shown)
Abstract:
We present Chandra and VLA observations of GW170817 at ~521-743 days post merger, and a homogeneous analysis of the entire Chandra data set. We find that the late-time non-thermal emission follows the expected evolution from an off-axis relativistic jet, with a steep temporal decay $F_ν\propto t^{-1.95\pm0.15}$ and a simple power-law spectrum $F_ν\propto ν^{-0.575\pm0.007}$. We present a new metho…
▽ More
We present Chandra and VLA observations of GW170817 at ~521-743 days post merger, and a homogeneous analysis of the entire Chandra data set. We find that the late-time non-thermal emission follows the expected evolution from an off-axis relativistic jet, with a steep temporal decay $F_ν\propto t^{-1.95\pm0.15}$ and a simple power-law spectrum $F_ν\propto ν^{-0.575\pm0.007}$. We present a new method to constrain the merger environment density based on diffuse X-ray emission from hot plasma in the host galaxy and we find $n\le 9.6 \times 10^{-3}\,\rm{cm^{-3}}$. This measurement is independent from inferences based on the jet afterglow modeling and allows us to partially solve for model degeneracies. The updated best-fitting model parameters with this density constraint are a fireball kinetic energy $E_0 = 1.5_{-1.1}^{+3.6}\times 10^{49}\,\rm{erg}$ ($E_{iso}= 2.1_{-1.5}^{+6.4}\times10^{52}\, \rm{erg}$), jet opening angle $θ_{0}= 5.9^{+1.0}_{-0.7}\,\rm{deg}$ with characteristic Lorentz factor $Γ_j = 163_{-43}^{+23}$, expanding in a low-density medium with $n_0 = 2.5_{-1.9}^{+4.1} \times 10^{-3}\, \rm{cm^{-3}}$ and viewed $θ_{obs} = 30.4^{+4.0}_{-3.4}\, \rm{deg}$ off-axis. The synchrotron emission originates from a power-law distribution of electrons with $p=2.15^{+0.01}_{-0.02}$. The shock microphysics parameters are constrained to $ε_{\rm{e}} = 0.18_{-0.13}^{+0.30}$ and $ε_{\rm{B}}=2.3_{-2.2}^{+16.0} \times 10^{-3}$. We investigate the presence of X-ray flares and find no statistically significant evidence of $\ge2.5σ$ of temporal variability at any time. Finally, we use our observations to constrain the properties of synchrotron emission from the deceleration of the fastest kilonova ejecta with energy $E_k^{KN}\propto (Γβ)^{-α}$ into the environment, finding that shallow stratification indexes $α\le6$ are disfavored.
△ Less
Submitted 8 November, 2019; v1 submitted 13 September, 2019;
originally announced September 2019.
-
The Optical Afterglow of GW170817: An Off-axis Structured Jet and Deep Constraints on a Globular Cluster Origin
Authors:
Wen-fai Fong,
P. K. Blanchard,
K. D. Alexander,
J. Strader,
R. Margutti,
A. Hajela,
V. A. Villar,
Y. Wu,
C. S. Ye,
E. Berger,
R. Chornock,
D. Coppejans,
P. S. Cowperthwaite,
T. Eftekhari,
D. Giannios,
C. Guidorzi,
A. Kathirgamaraju,
T. Laskar,
A. MacFadyen,
B. D. Metzger,
M. Nicholl,
K. Paterson,
G. Terreran,
D. J. Sand,
L. Sironi
, et al. (3 additional authors not shown)
Abstract:
We present a revised and complete optical afterglow light curve of the binary neutron star merger GW170817, enabled by deep Hubble Space Telescope (HST) F606W observations at $\approx\!584$ days post-merger, which provide a robust optical template. The light curve spans $\approx 110-362$ days, and is fully consistent with emission from a relativistic structured jet viewed off-axis, as previously i…
▽ More
We present a revised and complete optical afterglow light curve of the binary neutron star merger GW170817, enabled by deep Hubble Space Telescope (HST) F606W observations at $\approx\!584$ days post-merger, which provide a robust optical template. The light curve spans $\approx 110-362$ days, and is fully consistent with emission from a relativistic structured jet viewed off-axis, as previously indicated by radio and X-ray data. Combined with contemporaneous radio and X-ray observations, we find no spectral evolution, with a weighted average spectral index of $\langle β\rangle = -0.583 \pm 0.013$, demonstrating that no synchrotron break frequencies evolve between the radio and X-ray bands over these timescales. We find that an extrapolation of the post-peak temporal slope of GW170817 to the luminosities of cosmological short GRBs matches their observed jet break times, suggesting that their explosion properties are similar, and that the primary difference in GW170817 is viewing angle. Additionally, we place a deep limit on the luminosity and mass of an underlying globular cluster of $L \lesssim 6.7 \times 10^{3}\,L_{\odot}$, or $M \lesssim 1.3 \times 10^{4}\,M_{\odot}$, at least 4 standard deviations below the peak of the globular cluster mass function of the host galaxy, NGC4993. This limit provides a direct and strong constraint that GW170817 did not form and merge in a globular cluster. As highlighted here, HST (and soon JWST) enables critical observations of the optical emission from neutron star merger jets and outflows.
△ Less
Submitted 21 August, 2019;
originally announced August 2019.
-
Inter-Plasmoid Compton Scattering and the Compton Dominance of BL Lacs
Authors:
I. M. Christie,
M. Petropoulou,
L. Sironi,
D. Giannios
Abstract:
Blazar emission models based on magnetic reconnection succeed in reproducing many observed spectral and temporal features, including the short-duration luminous flaring events. Plasmoids, a self-consistent by-product of the tearing instability in the reconnection layer, can be the main source of blazar emission. Kinetic simulations of relativistic reconnection have demonstrated that plasmoids are…
▽ More
Blazar emission models based on magnetic reconnection succeed in reproducing many observed spectral and temporal features, including the short-duration luminous flaring events. Plasmoids, a self-consistent by-product of the tearing instability in the reconnection layer, can be the main source of blazar emission. Kinetic simulations of relativistic reconnection have demonstrated that plasmoids are characterized by rough energy equipartition between their radiating particles and magnetic fields. This is the main reason behind the apparent shortcoming of plasmoid-dominated emission models to explain the observed Compton ratios of BL Lac objects. Here, we demonstrate that the radiative interactions among plasmoids, which have been neglected so far, can assist in alleviating this contradiction. We show that photons emitted by large, slow-moving plasmoids can be a potentially important source of soft photons to be then up-scattered, via inverse Compton, by small fast-moving, neighboring plasmoids. This inter-plasmoid Compton scattering process can naturally occur throughout the reconnection layer, imprinting itself as an increase in the observed Compton ratios from those short and luminous plasmoid-powered flares within BL Lac sources, while maintaining energy equipartition between radiating particles and magnetic fields.
△ Less
Submitted 6 November, 2019; v1 submitted 7 August, 2019;
originally announced August 2019.
-
Relativistic Magnetic Reconnection in Electron-Positron-Proton Plasmas: Implications for Jets of Active Galactic Nuclei
Authors:
Maria Petropoulou,
Lorenzo Sironi,
Anatoly Spitkovsky,
Dimitrios Giannios
Abstract:
Magnetic reconnection is often invoked to explain the non-thermal radiation of relativistic outflows, including jets of active galactic nuclei (AGN). Motivated by the largely unknown plasma composition of AGN jets, we study reconnection in the unexplored regime of electron-positron-proton (pair-proton) plasmas with large-scale two-dimensional particle-in-cell simulations. We cover a wide range of…
▽ More
Magnetic reconnection is often invoked to explain the non-thermal radiation of relativistic outflows, including jets of active galactic nuclei (AGN). Motivated by the largely unknown plasma composition of AGN jets, we study reconnection in the unexplored regime of electron-positron-proton (pair-proton) plasmas with large-scale two-dimensional particle-in-cell simulations. We cover a wide range of pair multiplicities (lepton-to-proton number ratio $κ=1-199$) for different values of the all-species plasma magnetization ($σ=1,3$ and 10) and electron temperature ($Θ_e\equiv kT_e/m_ec^2=0.1-100$). We focus on the dependence of the post-reconnection energy partition and lepton energy spectra on the hot pair plasma magnetization $σ_{e,h}$ (i.e., the ratio of magnetic to pair enthalpy densities). We find that the post-reconnection energy is shared roughly equally between magnetic fields, pairs, and protons for $σ_{e,h}\gtrsim 3$. We empirically find that the mean lepton Lorentz factor in the post-reconnection region depends on $σ, Θ_e$, and $σ_{e,h}$ as $\langle γ_e-1\rangle \approx \sqrtσ(1+4Θ_e)\left(1+σ_{e,h}/30\right)$, for $σ\ge1$. The high-energy part of the post-reconnection lepton energy distributions can be described by a power law, whose slope is mainly controlled by $σ_{e,h}$ for $κ\gtrsim 3-6$, with harder power laws obtained for higher magnetizations. We finally show that reconnection in pair-proton plasmas with multiplicities $κ\sim 1-20$, magnetizations $σ\sim 1-10$, and temperatures $Θ_e \sim 1-10$ results in particle power law slopes and average electron Lorentz factors that are consistent with those inferred in leptonic models of AGN jet emission.
△ Less
Submitted 7 June, 2019;
originally announced June 2019.
-
Astro2020 Science White Paper: Using X-Ray Polarimetry to Probe the Physics of Black Holes and Neutron Stars
Authors:
Henric Krawczynski,
Giorgio Matt,
Adam R. Ingram,
Roberto Taverna,
Roberto Turolla,
Fabian Kislat,
C. C. Teddy Cheung,
Andrei Bykov,
Kuver Sinha,
Haocheng Zhang,
Jeremy Heyl,
Niccolo Bucciantini,
Greg Madejski,
Tim Kallman,
Keith M. Jahoda,
Quin Abarr,
Matthew G. Baring,
Luca Baldini,
Mitchell Begelman,
Markus Boettcher,
Edward Cackett,
Ilaria Caiazzo,
Paolo Coppi,
Enrico Costa,
Jason Dexter
, et al. (32 additional authors not shown)
Abstract:
This white paper highlights compact object and fundamental physics science opportunities afforded by high-throughput broadband (0.1-60 keV) X-ray polarization observations. X-ray polarimetry gives new observables with geometric information about stellar remnants which are many orders of magnitude too small for direct imaging. The X-ray polarimetric data also reveal details about the emission mecha…
▽ More
This white paper highlights compact object and fundamental physics science opportunities afforded by high-throughput broadband (0.1-60 keV) X-ray polarization observations. X-ray polarimetry gives new observables with geometric information about stellar remnants which are many orders of magnitude too small for direct imaging. The X-ray polarimetric data also reveal details about the emission mechanisms and the structure of the magnetic fields in and around the most extreme objects in the Universe. Whereas the Imaging X-ray Polarimetry Explorer (IXPE) to be launched in 2021 will obtain first results for bright objects, a follow-up mission could be one order of magnitude more sensitive and would be able to use a broader bandpass to perform physics type experiments for representative samples of sources.
△ Less
Submitted 19 April, 2019;
originally announced April 2019.
-
Neutrinos, Cosmic Rays and the MeV Band
Authors:
R. Ojha,
H. Zhang,
M. Kadler,
N. K. Neilson,
M. Kreter,
J. McEnery,
S. Buson,
R. Caputo,
P. Coppi,
F. D'Ammando,
A. De Angelis,
K. Fang,
D. Giannios,
S. Guiriec,
F. Guo,
J. Kopp,
F. Krauss,
H. Li,
M. Meyer,
A. Moiseev,
M. Petropoulou,
C. Prescod-Weinstein,
B. Rani,
C. Shrader,
T. Venters
, et al. (1 additional authors not shown)
Abstract:
The possible association of the blazar TXS 0506+056 with a high-energy neutrino detected by IceCube holds the tantalizing potential to answer three astrophysical questions: 1. Where do high-energy neutrinos originate? 2. Where are cosmic rays produced and accelerated? 3. What radiation mechanisms produce the high-energy γ-rays in blazars? The MeV gamma-ray band holds the key to these questions, be…
▽ More
The possible association of the blazar TXS 0506+056 with a high-energy neutrino detected by IceCube holds the tantalizing potential to answer three astrophysical questions: 1. Where do high-energy neutrinos originate? 2. Where are cosmic rays produced and accelerated? 3. What radiation mechanisms produce the high-energy γ-rays in blazars? The MeV gamma-ray band holds the key to these questions, because it is an excellent proxy for photo-hadronic processes in blazar jets, which also produce neutrino counterparts. Variability in MeV gamma-rays sheds light on the physical conditions and mechanisms that take place in the particle acceleration sites in blazar jets. In addition, hadronic blazar models also predict a high level of polarization fraction in the MeV band, which can unambiguously distinguish the radiation mechanism. Future MeV missions with a large field of view, high sensitivity, and polarization capabilities will play a central role in multi-messenger astronomy, since pointed, high-resolution telescopes will follow neutrino alerts only when triggered by an all-sky instrument.
△ Less
Submitted 13 March, 2019;
originally announced March 2019.
-
Extreme Plasma Astrophysics
Authors:
D. Uzdensky,
M. Begelman,
A. Beloborodov,
R. Blandford,
S. Boldyrev,
B. Cerutti,
F. Fiuza,
D. Giannios,
T. Grismayer,
M. Kunz,
N. Loureiro,
M. Lyutikov,
M. Medvedev,
M. Petropoulou,
A. Philippov,
E. Quataert,
A. Schekochihin,
K. Schoeffler,
L. Silva,
L. Sironi,
A. Spitkovsky,
G. Werner,
V. Zhdankin,
J. Zrake,
E. Zweibel
Abstract:
This is a science white paper submitted to the Astro-2020 and Plasma-2020 Decadal Surveys. The paper describes the present status and emerging opportunities in Extreme Plasma Astrophysics -- a study of astrophysically-relevant plasma processes taking place under extreme conditions that necessitate taking into account relativistic, radiation, and QED effects.
This is a science white paper submitted to the Astro-2020 and Plasma-2020 Decadal Surveys. The paper describes the present status and emerging opportunities in Extreme Plasma Astrophysics -- a study of astrophysically-relevant plasma processes taking place under extreme conditions that necessitate taking into account relativistic, radiation, and QED effects.
△ Less
Submitted 20 July, 2019; v1 submitted 13 March, 2019;
originally announced March 2019.
-
Multi-Messenger Astronomy with Extremely Large Telescopes
Authors:
Ryan Chornock,
Philip S. Cowperthwaite,
Raffaella Margutti,
Dan Milisavljevic,
Kate D. Alexander,
Igor Andreoni,
Iair Arcavi,
Adriano Baldeschi,
Jennifer Barnes,
Eric Bellm,
Paz Beniamini,
Edo Berger,
Christopher P. L. Berry,
Federica Bianco,
Peter K. Blanchard,
Joshua S. Bloom,
Sarah Burke-Spolaor,
Eric Burns,
Dario Carbone,
S. Bradley Cenko,
Deanne Coppejans,
Alessandra Corsi,
Michael Coughlin,
Maria R. Drout,
Tarraneh Eftekhari
, et al. (60 additional authors not shown)
Abstract:
The field of time-domain astrophysics has entered the era of Multi-messenger Astronomy (MMA). One key science goal for the next decade (and beyond) will be to characterize gravitational wave (GW) and neutrino sources using the next generation of Extremely Large Telescopes (ELTs). These studies will have a broad impact across astrophysics, informing our knowledge of the production and enrichment hi…
▽ More
The field of time-domain astrophysics has entered the era of Multi-messenger Astronomy (MMA). One key science goal for the next decade (and beyond) will be to characterize gravitational wave (GW) and neutrino sources using the next generation of Extremely Large Telescopes (ELTs). These studies will have a broad impact across astrophysics, informing our knowledge of the production and enrichment history of the heaviest chemical elements, constrain the dense matter equation of state, provide independent constraints on cosmology, increase our understanding of particle acceleration in shocks and jets, and study the lives of black holes in the universe. Future GW detectors will greatly improve their sensitivity during the coming decade, as will near-infrared telescopes capable of independently finding kilonovae from neutron star mergers. However, the electromagnetic counterparts to high-frequency (LIGO/Virgo band) GW sources will be distant and faint and thus demand ELT capabilities for characterization. ELTs will be important and necessary contributors to an advanced and complete multi-messenger network.
△ Less
Submitted 11 March, 2019;
originally announced March 2019.
-
High-Energy Polarimetry - a new window to probe extreme physics in AGN jets
Authors:
B. Rani,
H. Zhang,
S. D. Hunter,
F. Kislat,
M. Böttcher,
J. E. McEnery,
D. J. Thompson,
D. Giannios,
F. Guo,
H. Li,
M. Baring,
I. Agudo,
S. Buson,
M. Petropoulou,
V. Pavlidou,
E. Angelakis,
I. Myserlis,
Z. Wadiasingh,
R. M. Curado da Silva,
P. Kilian,
S. Guiriec,
V. V. Bozhilov,
J. Hodgson,
S. Antón,
D. Kazanas
, et al. (9 additional authors not shown)
Abstract:
The constantly improving sensitivity of ground-based and space-borne observatories has made possible the detection of high-energy emission (X-rays and gamma-rays) from several thousands of extragalactic sources. Enormous progress has been made in measuring the continuum flux enabling us to perform imaging, spectral and timing studies. An important remaining challenge for high-energy astronomy is m…
▽ More
The constantly improving sensitivity of ground-based and space-borne observatories has made possible the detection of high-energy emission (X-rays and gamma-rays) from several thousands of extragalactic sources. Enormous progress has been made in measuring the continuum flux enabling us to perform imaging, spectral and timing studies. An important remaining challenge for high-energy astronomy is measuring polarization. The capability to measure polarization is being realized currently at X-ray energies (e.g. with IXPE), and sensitive gamma-ray telescopes capable of measuring polarization, such as AMEGO, AdEPT, e-ASTROGAM, etc., are being developed. These future gamma-ray telescopes will probe the radiation mechanisms and magnetic fields of relativistic jets from active galactic nuclei at spatial scales much smaller than the angular resolution achieved with continuum observations of the instrument. In this white paper, we discuss the scientific potentials of high-energy polarimetry, especially gamma-ray polarimetry, including the theoretical implications, and observational technology advances being made. In particular, we will explore the primary scientific opportunities and wealth of information expected from synergy of multi-wavelength polarimetry that will be brought to multi-messenger astronomy.
△ Less
Submitted 11 March, 2019;
originally announced March 2019.
-
Multi-Physics of AGN Jets in the Multi-Messenger Era
Authors:
B. Rani,
M. Petropoulou,
H. Zhang,
F. D'Ammando,
J. Finke,
M. Baring,
M. Böttcher,
S. Dimitrakoudis,
Z. Gan,
D. Giannios,
D. H. Hartmann,
T. P. Krichbaum,
A. P. Marscher,
A. Mastichiadis,
K. Nalewajko,
R. Ojha,
D. Paneque,
C. Shrader,
L. Sironi,
A. Tchekhovskoy,
D. J. Thompson,
N. Vlahakis,
T. M. Venters
Abstract:
Active galactic nuclei (AGN) with relativistic jets, powered by gas accretion onto their central supermassive black hole (SMBH), are unique laboratories for studying the physics of matter and elementary particles in extreme conditions that cannot be realized on Earth. For a long time since the discovery of AGN, photons were the only way to probe the underlying physical processes. The recent discov…
▽ More
Active galactic nuclei (AGN) with relativistic jets, powered by gas accretion onto their central supermassive black hole (SMBH), are unique laboratories for studying the physics of matter and elementary particles in extreme conditions that cannot be realized on Earth. For a long time since the discovery of AGN, photons were the only way to probe the underlying physical processes. The recent discovery of a very high energy neutrino, IceCube-170922A, coincident with a flaring blazar, TXS 0506+056, provides the first evidence that AGN jets are multi-messenger sources; they are capable of accelerating hadrons to very high energies, while producing non-thermal EM radiation and high-energy neutrinos. This new era of multi-messenger astronomy, which will mature in the next decade, offers us the unprecedented opportunity to combine more than one messenger to solve some long-standing puzzles of AGN jet physics: How do jets dissipate their energy to accelerate particles? What is the jet total kinetic power? Where and how do jets produce high-energy emission and neutrinos? What physical mechanisms drive the particle acceleration?
△ Less
Submitted 11 March, 2019;
originally announced March 2019.
-
Probing the Emission Mechanism and Magnetic Field of Neutrino Blazars with Multi-Wavelength Polarization Signatures
Authors:
Haocheng Zhang,
Ke Fang,
Hui Li,
Dimitrios Giannios,
Markus Böttcher,
Sara Buson
Abstract:
The characteristic two-component blazar spectral energy distribution (SED) can be of either leptonic and/or hadronic origins. The potential association of the high-energy neutrino event IceCube-170922A with the flaring blazar TXS~0506+056 indicates that hadronic processes may operate in a blazar jet. Despite multi-wavelength follow-ups of the event and extensive theoretical modelings, the radiatio…
▽ More
The characteristic two-component blazar spectral energy distribution (SED) can be of either leptonic and/or hadronic origins. The potential association of the high-energy neutrino event IceCube-170922A with the flaring blazar TXS~0506+056 indicates that hadronic processes may operate in a blazar jet. Despite multi-wavelength follow-ups of the event and extensive theoretical modelings, the radiation mechanisms and the underlying magnetic field strength and configuration remain poorly understood. In this paper, we consider generic leptonic and hadronic blazar spectral models with distinct magnetic field strengths and radiation mechanisms. We analytically reproduce the SEDs and the neutrino flux of hadronic models, and predict their X-ray to $γ$-ray polarization degrees. Furthermore, by performing relativistic magnetohydrodynamic (RMHD) simulations taking into account the polarization-dependent radiation transfer, we study the time-dependent multi-wavelength polarization variability of the proton synchrotron model under a shock scenario. Our results suggest that the high-energy polarization degree and the neutrino flux can be jointly used to pinpoint the leptonic and/or hadronic blazar radiation mechanisms in the X-ray and $γ$-ray bands, and to infer the magnetic field strength in the emission region. Additionally, the temporal multi-wavelength polarization signatures in the proton synchrotron model shed light on the jet energy composition and the dynamical importance of magnetic fields in the blazar emission region. Future multi-wavelength polarimetry facilities such as {\it IXPE} and {\it AMEGO} together with neutrino telescopes such as {\it IceCube} can provide unprecedented observational constraints to probe the blazar radiation mechanisms and jet dynamics.
△ Less
Submitted 5 March, 2019;
originally announced March 2019.
-
Observational Constraints on Late-Time Radio Rebrightening of GRB-Supernovae
Authors:
Charee Peters,
Alexander J. van der Horst,
Laura Chomiuk,
Adithan Kathirgamaraju,
Rodolfo Barniol Duran,
Dimitrios Giannios,
Cormac Reynolds,
Zsolt Paragi,
Eric Wilcots
Abstract:
We present a search for late-time rebrightening of radio emission from three supernovae (SNe) with associated gamma-ray bursts (GRBs). It has been previously proposed that the unusually energetic SNe associated with GRBs should enter the Sedov-Taylor phase decades after the stellar explosion, and this SN "remnant" emission will outshine the GRB radio afterglow and be detectable at significant dist…
▽ More
We present a search for late-time rebrightening of radio emission from three supernovae (SNe) with associated gamma-ray bursts (GRBs). It has been previously proposed that the unusually energetic SNe associated with GRBs should enter the Sedov-Taylor phase decades after the stellar explosion, and this SN "remnant" emission will outshine the GRB radio afterglow and be detectable at significant distances. We place deep limits on the radio luminosity of GRB 980425/SN 1998bw, GRB 030329/SN 2003dh and GRB 060218/SN 2006aj, 10-18 years after explosion, with our deepest limit being $L_ν$ $< 4 \times 10^{26}$ erg s$^{-1}$ Hz$^{-1}$ for GRB 980425/SN 1998bw. We put constraints on the density of the surrounding medium for various assumed values of the microphysical parameters related to the magnetic field and synchrotron-emitting electrons. For GRB 060218/SN 2006aj and GRB 980425/SN 1998bw, these density limits have implications for the density profile of the surrounding medium, while the non-detection of GRB 030329/SN 2003dh implies that its afterglow will not be detectable anymore at GHz frequencies.
△ Less
Submitted 4 January, 2019;
originally announced January 2019.
-
Observable Features of GW170817 Kilonova Afterglow
Authors:
Adithan Kathirgamaraju,
Dimitrios Giannios,
Paz Beniamini
Abstract:
The neutron star merger, GW170817, was followed by an optical-infrared transient (a kilonova) which indicated that a substantial ejection of mass at trans-relativistic velocities occurred during the merger. Modeling of the kilonova is able to constrain the kinetic energy of the ejecta and its characteristic velocity, but not the high-velocity distribution of the ejecta. Yet, this distribution cont…
▽ More
The neutron star merger, GW170817, was followed by an optical-infrared transient (a kilonova) which indicated that a substantial ejection of mass at trans-relativistic velocities occurred during the merger. Modeling of the kilonova is able to constrain the kinetic energy of the ejecta and its characteristic velocity, but not the high-velocity distribution of the ejecta. Yet, this distribution contains crucial information on the merger dynamics. In this work, we assume a power-law distribution of the form $E(>βΓ)\propto(βΓ)^{-α}$ for the energy of the kilonova ejecta and calculate the non-thermal signatures produced by the interaction of the ejecta with the ambient gas. We find that ejecta with minimum velocity $β_0\simeq 0.3$ and energy $E\sim 10^{51}$ erg, as inferred from kilonova modeling, has a detectable radio, and possibly X-ray, afterglow for a broad range of parameter space. This afterglow component is expected to dominate the observed emission on a timescale of a few years post merger and peak around a decade later. Its light curve can be used to determine properties of the kilonova ejecta and in particular the ejecta velocity distribution $α$, the minimum velocity $β_0$ and its total kinetic energy $E$. We also predict that an afterglow rebrightening, that is associated with the kilonova component, will be accompanied by a shift of the centroid of the radio source towards the initial position of the explosion.
△ Less
Submitted 3 January, 2019;
originally announced January 2019.
-
Gamma-ray bursts as cool synchrotron sources
Authors:
J. Michael Burgess,
Damien Bégué,
Ana Bacelj,
Dimitrios Giannios,
Francesco Berlato,
Jochen Greiner
Abstract:
Gamma-ray bursts are the most energetic electromagnetic sources in the Universe. Their prompt gamma-ray radiation corresponds to an energy release of 1E42-1E47J. Fifty years after their discovery and several dedicated space-based instruments, the physical origin of this emission is still unknown. Synchrotron emission has been one of the early contenders but was criticized because spectral fits of…
▽ More
Gamma-ray bursts are the most energetic electromagnetic sources in the Universe. Their prompt gamma-ray radiation corresponds to an energy release of 1E42-1E47J. Fifty years after their discovery and several dedicated space-based instruments, the physical origin of this emission is still unknown. Synchrotron emission has been one of the early contenders but was criticized because spectral fits of empirical models (such as a smoothly-connected broken power law or a cut-off power law) suggest too hard a slope of the low-energy power law, violating the so-called synchrotron line-of-death. We perform time-resolved gamma-ray spectroscopy of single-peaked GRBs as measured with Fermi/GBM. We demonstrate that idealized synchrotron emission, when properly incorporating time-dependent cooling of the electrons, is capable of fitting ~95% of all these GBM spectra. The comparison with spectral fit results based on previous empirical models demonstrates that the past exclusion of synchrotron radiation as an emission mechanism derived via the line-of-death was misleading. Our analysis probes the physics of these ultra-relativistic outflows and the related microphysical processes, and for the first time provides estimates of magnetic field strength and Lorentz factors of the emitting region directly from spectral fits. Our modeling of the Fermi/GBM observations provides evidence that GRBs are produced by moderately magnetized jets in which relativistic mini-jets emit optically-thin synchrotron radiation at large emission radii.
△ Less
Submitted 16 October, 2018;
originally announced October 2018.
-
EM counterparts of structured jets from 3D GRMHD simulations
Authors:
Adithan Kathirgamaraju,
Alexander Tchekhovskoy,
Dimitrios Giannios,
Rodolfo Barniol Duran
Abstract:
GW170817/GRB170817A has offered unprecedented insight into binary neutron star post-merger systems. Its Prompt and afterglow emission imply the presence of a tightly collimated relativistic jet with a smooth transverse structure. However, it remains unclear whether and how the central engine can produce such structured jets. Here, we utilize 3D GRMHD simulations starting with a black hole surround…
▽ More
GW170817/GRB170817A has offered unprecedented insight into binary neutron star post-merger systems. Its Prompt and afterglow emission imply the presence of a tightly collimated relativistic jet with a smooth transverse structure. However, it remains unclear whether and how the central engine can produce such structured jets. Here, we utilize 3D GRMHD simulations starting with a black hole surrounded by a magnetized torus with properties typically expected of a post-merger system. We follow the jet, as it is self-consistently launched, from the scale of the compact object out to more than 3 orders of magnitude in distance. We find that this naturally results in a structured jet, which is collimated by the disk wind into a half-opening angle of roughly 10 degrees, its emission can explain features of both the prompt and afterglow emission of GRB170817A for a 30 degree observing angle. Our work is the first to compute the afterglow, in the context of a binary merger, from a relativistic magnetized jet self-consistently generated by an accreting black hole, with the jet's transverse structure determined by the accretion physics and not prescribed at any point.
△ Less
Submitted 2 December, 2018; v1 submitted 13 September, 2018;
originally announced September 2018.
-
Extreme scattering events from axisymmetric plasma lenses
Authors:
Lingyi Dong,
Maria Petropoulou,
Dimitrios Giannios
Abstract:
Frequency-dependent brightness fluctuations of radio sources, the so-called extreme scattering events (ESEs), have been observed over the last three decades. They are caused by Galactic plasma structures whose geometry and origin are still poorly understood. In this paper, we construct axisymmentric two-dimensional (2D) column density profiles for the plasma lens and explore the resulting ESEs for…
▽ More
Frequency-dependent brightness fluctuations of radio sources, the so-called extreme scattering events (ESEs), have been observed over the last three decades. They are caused by Galactic plasma structures whose geometry and origin are still poorly understood. In this paper, we construct axisymmentric two-dimensional (2D) column density profiles for the plasma lens and explore the resulting ESEs for both point-like and extended sources. A quantity that becomes relevant is the impact parameter $b$, namely the distance of the observer's path from the lens symmetry axis. We demonstrate its effects on the shape of ESE light curves and use it for a phenomenological classification of ESEs into four main types. Three of them are unique outcomes of the 2D model and do not show a characteristic U-shaped dip in the light curve, which has been traditionally used as an identification means of ESEs. We apply our model to five well-studied ESEs and show that elongated plasma tubes or quasi-spherical clouds are favoured over plasma sheets for four of them, while the remaining one is compatible with both lens geometries.
△ Less
Submitted 31 August, 2018;
originally announced September 2018.
-
A lesson from GW170817: most neutron star mergers result in tightly collimated successful GRB jets
Authors:
Paz Beniamini,
Maria Petropoulou,
Rodolfo Barniol Duran,
Dimitrios Giannios
Abstract:
The joint detection of gravitational waves (GWs) and $γ$-rays from a binary neutron-star (NS) merger provided a unique view of off-axis GRBs and an independent measurement of the NS merger rate. Comparing the observations of GRB170817 with those of the regular population of short GRBs (sGRBs), we show that an order unity fraction of NS mergers result in sGRB jets that breakout of the surrounding e…
▽ More
The joint detection of gravitational waves (GWs) and $γ$-rays from a binary neutron-star (NS) merger provided a unique view of off-axis GRBs and an independent measurement of the NS merger rate. Comparing the observations of GRB170817 with those of the regular population of short GRBs (sGRBs), we show that an order unity fraction of NS mergers result in sGRB jets that breakout of the surrounding ejecta. We argue that the luminosity function of sGRBs, peaking at $\approx 2\times 10^{52}\, \mbox{erg s}^{-1}$, is likely an intrinsic property of the sGRB central engine and that sGRB jets are typically narrow with opening angles $θ_0 \approx 0.1$. We perform Monte Carlo simulations to examine models for the structure and efficiency of the prompt emission in off axis sGRBs. We find that only a small fraction ($\sim 0.01-0.1$) of NS mergers detectable by LIGO/VIRGO in GWs is expected to be also detected in prompt $γ$-rays and that GW170817-like events are very rare. For a NS merger rate of $\sim 1500$ Gpc$^{-3}$ yr$^{-1}$, as inferred from GW170817, we expect within the next decade up to $\sim 12$ joint detections with off-axis GRBs for structured jet models and just $\sim 1$ for quasi-spherical cocoon models where $γ$-rays are the result of shock breakout. Given several joint detections and the rates of their discoveries, the different structure models can be distinguished. In addition the existence of a cocoon with a reservoir of thermal energy may be observed directly in the UV, given a sufficiently rapid localisation of the GW source.
△ Less
Submitted 10 December, 2018; v1 submitted 14 August, 2018;
originally announced August 2018.
-
Large-Amplitude Blazar Polarization Angle Swing as a Signature of Magnetic Reconnection
Authors:
Haocheng Zhang,
Xiaocan Li,
Fan Guo,
Dimitrios Giannios
Abstract:
Relativistic magnetic reconnection events can widely exist in magnetized plasmas in astrophysical systems. During this process, oppositely directed magnetic field lines reconnect and release magnetic energy, efficiently accelerating nonthermal particles. However, so far there is little clear observational signatures of relativistic magnetic reconnection events in astrophysical systems. Blazars are…
▽ More
Relativistic magnetic reconnection events can widely exist in magnetized plasmas in astrophysical systems. During this process, oppositely directed magnetic field lines reconnect and release magnetic energy, efficiently accelerating nonthermal particles. However, so far there is little clear observational signatures of relativistic magnetic reconnection events in astrophysical systems. Blazars are relativistic magnetized plasma outflows from supermassive black holes. Their multi-wavelength flares may be powered by relativistic magnetic reconnection. The highly variable radiation and polarization signatures are well covered by multi-wavelength observation campaigns, making them ideal targets to examine the magnetic reconnection model. Recent observations have found that several blazar flares are accompanied by optical polarization angle swings which may be of as large amplitude as $> 180^{\circ}$, challenging existing theoretical models. In this paper, we present integrated particle-in-cell (PIC) and polarized radiation transfer simulations of magnetic reconnection events. We find that plasmoid coalescences in the reconnection layer can give rise to highly variable light curves, low and fluctuating polarization degree, and rotating polarization angle. In particular, large-amplitude polarization angle swings, similar to those observed during blazar flares, can be a unique signature of relativistic magnetic reconnection events.
△ Less
Submitted 22 July, 2018;
originally announced July 2018.
-
Radiative Signatures of Plasmoid-Dominated Reconnection in Blazar Jets
Authors:
I. M. Christie,
M. Petropoulou,
L. Sironi,
D. Giannios
Abstract:
The multi-wavelength spectral and temporal variability observed in blazars set tight constraints on current theoretical emission models. Here, we investigate the relativistic magnetic reconnection process as a source of blazar emission in which quasi-spherical plasmoids, containing relativistic particles and magnetic fields, are associated with the emission sites in blazar jets. By coupling recent…
▽ More
The multi-wavelength spectral and temporal variability observed in blazars set tight constraints on current theoretical emission models. Here, we investigate the relativistic magnetic reconnection process as a source of blazar emission in which quasi-spherical plasmoids, containing relativistic particles and magnetic fields, are associated with the emission sites in blazar jets. By coupling recent two-dimensional particle-in-cell simulations of relativistic reconnection with a time-dependent radiative transfer code, we compute the non-thermal emission from a chain of plasmoids formed during a reconnection event. The derived photon spectra display characteristic features observed in both BL Lac sources and flat spectrum radio quasars, with the distinction made by varying the strength of the external photon fields, the jet magnetization, and the number of pairs per proton contained within. Light curves produced from reconnection events are composed of many fast and powerful flares that appear on excess of a slower evolving envelope produced by the cumulative emission of medium-sized plasmoids. The observed variability is highly dependent upon the orientation of the reconnection layer with respect to the blazar jet axis and to the observer. Our model provides a physically motivated framework for explaining the multi-timescale blazar variability across the entire electromagnetic spectrum.
△ Less
Submitted 3 December, 2018; v1 submitted 20 July, 2018;
originally announced July 2018.
-
GRB and blazar jets shining through their stripes
Authors:
Dimitrios Giannios,
Dmitri A. Uzdensky
Abstract:
Black-hole driven relativistic astrophysical jets, such as blazars and gamma-ray bursts (GRBs), are powerful sources of electromagnetic radiation. Their emission is powered by some energy dissipation and particle acceleration mechanism operating in-situ at large distances from the black hole. We propose that the formation of the dissipative structures in the jet is controlled by the time variabili…
▽ More
Black-hole driven relativistic astrophysical jets, such as blazars and gamma-ray bursts (GRBs), are powerful sources of electromagnetic radiation. Their emission is powered by some energy dissipation and particle acceleration mechanism operating in-situ at large distances from the black hole. We propose that the formation of the dissipative structures in the jet is controlled by the time variability of the accretion disc. We argue that the open magnetic field lines through the black hole, which drive a strongly magnetized jet, may have their polarity reversing over time-scales related to the growth of the magneto-rotational dynamo in the disc. Consequently, the jet is characterized by an alternating toroidal field polarity along its propagation axis, i.e., it is a "striped jet". Magnetic reconnection in the current sheets that form between the stripes dissipates the alternating-field energy and thus powers further jet acceleration. Here, we consider a jet with a broad distribution of stripe widths $l>l_{\rm min}$, above a dominant scale $l_{\rm min}$. We find that the bulk acceleration of the jet, driven by the annihilation of the stripes, is very gradual. The dissipation rate peaks at a distance $z_{\rm peak}\sim 10^6 R_g\, (Γ_\infty/30)^2\, (l_{\rm min}/1000R_g)$, where $R_{g}$ is the black-hole's gravititional radius and $Γ_{\infty}$ the jet's asymptotic Lorentz factor, and exhibits a very broad plateau extending by $\sim 4-5$ orders of magnitude in distance. The prolonged energy dissipation accounts for the flat-to-inverted long-wavelength spectra commonly observed in jets. The model can also account for the broad range of flaring timescales of blazars and the fact that their bulk acceleration appears to continue out to ~100 pc scales. In GRB jets, the model predicts comparable power for the photopsheric and Thomson-thin emission components.
△ Less
Submitted 23 May, 2018;
originally announced May 2018.