Search

We investigate the external reverse shock region of relativistic jets as the origin of X-ray afterglows of jetted tidal disruption events (TDEs) that exhibit luminous jets accompanied by fast-declining non-thermal X-ray emissions. We model the dynamics of jet propagating within an external density medium, accounting for continuous energy injection driven by accretion activities. We compute the time-dependent synchrotron and inverse Compton emissions from the reverse shock region. Our analysis demonstrates that the reverse shock scenario can potentially explain the X-ray light curves and spectra of four jetted TDEs, AT 2022cmc, Swift J1644, Swift J2058, and Swift J1112. Notably, the rapid steepening of the late-stage X-ray light curves can be attributed jointly to the jet break and cessation of the central engine as the accretion rate drops below the Eddington limit. Using parameters obtained from X-ray data fitting, we also discuss the prospects for $\gamma$-ray and neutrino detection., Comment: 10 pages, 3 figures, 1 table. Comments are welcome

Tidal disruption events (TDEs), where stars are captured or tidally disrupted by supermassive black holes, are potential sources of high-energy neutrinos. We report the discovery of a potential neutrino flare that is spatially and temporally associated with X-ray emission from TDE ATLAS17jrp. The best-fit spectrum of the neutrino flare follows a power-law with an index of $\rm{\gamma=2.7\pm0.4}$ and a flux normalization of $\rm{\Phi_0 =1.7^{+6.3}_{-1.5}\times 10^{-18}\;GeV^{-1} cm^{-2} s^{-1}}$ at 100 TeV based on an analysis of 10-year track data from IceCube, and the flare duration is 61 days. We calculate that the probability of this association occurring by chance is $0.17\%$. Therefore, ATLAS17jrp is the second TDE (not including candidates) associated with high-energy neutrinos, following TDE AT2019dsg associated with an IceCube neutrino alert. This association can be attributed to the interaction of X-ray photons produced by the hot corona with high-energy particles accelerated by disk winds or outflows, resulting in the production of neutrinos., Comment: 8 pages, 2 figures, submitted to ApJL

We revisit the Ultra-High-Energy Cosmic Ray (UHECRs) production in Tidal Disruption Events (TDEs) in the light of recent neutrino-TDE associations. We use an isotropically emitting source-propagation model, which has been developed to describe the neutrino production in AT2019dsg, AT2019fdr, and AT2019aalc. These TDEs have strong dust echoes in the infrared range, which are potentially linked with the neutrino production. A mechanism where neutrinos originate from cosmic ray scattering on infrared photons implies cosmic rays in the ultra-high energy range, thus suggesting a natural connection with the observed UHECR. We extrapolate the three TDE associations to a population of neutrino- and UHECR-emitting TDEs, and postulate that these TDEs power the UHECRs. We then infer the source composition, population parameters, and local rates that are needed to describe UHECR data. We find that UHECR data point towards a mix of light to mid-heavy injection isotopes, which could be found, e.g., in oxygen-neon-magnesium white dwarfs, and to a contribution of at least two groups of TDEs with different characteristics, dominated by AT2019aalc-type events. The required local TDE rates of ${\mathcal O}(10^2)~\mathrm{Gpc^{-3} \, yr^{-1}}$, however, are more indicative of the disruption of main sequence stars. We propose an enhanced efficiency in the acceleration of heavier nuclei that could address this discrepancy. The predicted diffuse neutrino fluxes suggest a population of astrophysical neutrino sources that can be observed by future radio neutrino detection experiments., Comment: 27 pages, 15 figures

AT 2022cmc is a recently documented tidal disruption event (TDE) that exhibits a luminous jet, accompanied by fast-declining X-ray and long-lasting radio/millimeter emission. Motivated by the distinct spectral and temporal signatures between X-ray and radio observations, we propose a multizone model involving relativistic jets with different Lorentz factors. We systematically study the evolution of the faster and slower jets in an external density profile, considering the continuous energy injection rate associated with the time-dependent accretion rates before and after the mass fallback time. We investigate time-dependent multiwavelength emission from both the forward shock and reverse shock regions of the fast and slow jets, in a self-consistent manner. Our analysis demonstrates that the energy injection rate can significantly impact the jet evolution and subsequently influence the lightcurves. We find that the X-ray spectra and lightcurves can be described by the electron synchrotron emission from the reverse shock of the faster jet, in which the late-time X-ray upper limits, extending to 400 days after the disruption, could be interpreted as the jet break steepening. Meanwhile, the radio observations can be interpreted as a result of synchrotron emissions from the forward shock region of the slower jet. We also discuss prospects for testing the model with current and future observations., Comment: 17 pages, 5 figures and 1 table. Accepted for publication in ApJ

Recent multi-wavelength observations of gamma-ray burst afterglows observed in the TeV energy range challenge the simplest Synchrotron Self-Compton (SSC) interpretation of this emission and are consistent with a single power-law component spanning over eight orders of magnitude in energy. To interpret this generic behaviour in the single-zone approximation without adding further free parameters, we perform an exhaustive parameter space study using the public, time-dependent, multi-messenger transport software AM3. This description accounts for the radiation from non-thermal protons and the lepto-hadronic cascade induced by pp- and p{\gamma}-interactions. We summarise the main scenarios which we have found (SSC, Extended-syn, Proton-syn, pp-cascade, and p{\gamma}-cascade), and discuss their advantages and limitations. We find that possible high-density environments, as may be typical for surrounding molecular cloud material, offer an alternative explanation for producing flat hard (source) spectra up to and beyond energies of 10 TeV.

We discuss the possible association of an astrophysical neutrino (IC220405B) with the recently reported, extremely energetic tidal disruption event (TDE) candidate AT2021lwx (ZTF20abrbeie, aka ``Scary Barbie'') at redshift $z=0.995$. Although the TDE is about $2.6^\circ$ off the direction of the reconstructed neutrino event (outside the 90% confidence level localization region), the TDE candidate shares some important characteristics with so far reported neutrino-TDE associations: a strong infrared dust echo, high bolometric luminosity, a neutrino time delay with respect to the peak mass accretion rate of the order of hundred days, and a high observed X-ray luminosity. We interpret this new association using an isotropic emission model, where neutrinos are produced by the collision of accelerated protons with infrared photons. After accounting for the high redshift of AT2021lwx (by interpreting the data in the SMBH frame), we find that the expected neutrino fluences and neutrino time delays are qualitatively comparable to the other TDEs. Since data are only available up to 300 days post-peak in the SMBH frame, significant uncertainties exist in the dust echo interpretation, and therefore in the predicted number of neutrinos detected, $\mathcal N_{\nu}\simeq3.0\times10^{-3}-0.012$. We recommend further follow-up on this object for an extended period and suggest refining the reconstruction of the neutrino arrival direction in this particular case., Comment: 10 pages, 5 figures, 1 table. Accepted for publication in ApJ

We present the AM$^3$ ("Astrophysical Multi-Messenger Modeling") software. AM$^3$ is a documented open source software that efficiently solves the coupled integro-differential equations describing the temporal evolution of the spectral densities of particles interacting in astrophysical environments, in-cluding photons, electrons, positrons, protons, neutrons, pions, muons, and neutrinos. The software has been extensively used to simulate the multi-wavelength and neutrino emission from active galactic nuclei (including blazars), gamma-ray bursts, and tidal disruption events. The simulations include all relevant non-thermal processes, namely synchrotron emission, inverse Compton scattering, photon-photon annihilation, proton-proton and proton-photon pion production, and photo-pair production. The software self-consistently calculates the full cascade of primary and secondary particles, including non-linear feedback processes and predictions in the time domain. It also allows to track separately the particle densities produced by means of each distinct interaction processes, including the different hadronic channels. With its efficient hybrid solver combining analytical and numerical techniques, AM$^3$ combines efficiency and accuracy at a user-adjustable level. We describe the technical details of the numerical framework and present three examples of applications to different astrophysical environments., Comment: The BSD-licensed source code can be found under https://gitlab.desy.de/am3/am3, for the user guide and documentation see https://am3.readthedocs.io/en/latest/

The potential association between Tidal Disruption Events (TDEs) and high-energy astrophysical neutrinos implies the acceleration of cosmic rays. These accelerated particles will initiate electromagnetic (EM) cascades spanning from keV to GeV energies by the processes related to neutrino production. We model the EM cascade and neutrino emissions by numerically solving the time-dependent transport equations and discuss the implications for AT2019dsg and AT2019fdr in the X-ray and $\gamma$-ray bands. We show that the $\gamma$-ray constraints from \emph{Fermi} can constrain the size of the radiation zone and the maximum energy of injected protons, and that the corresponding expected neutrino event numbers in follow-up searches are limited to be less than about 0.1. Depending on the efficiency of $p\gamma$ interactions, the X-ray and $\gamma$-ray signals can be expected closer to the peak of the optical-ultraviolet (OUV) luminosity, or to the time of the neutrino production., Comment: 18 pages, 7 figures, 1 table. Revised version. Accepted for publication in ApJ

The detection of the hyper-bright gamma-ray burst (GRB) 221009A enables us to explore the nature of GRB emission and the origin of very-high-energy (VHE) gamma-rays. We analyze the ${\it Fermi}$-LAT data and investigate GeV-TeV emission in the framework of the external reverse shock model. We show that early $\sim1-10$ GeV emission can be explained by the external inverse-Compton mechanism via upscattering MeV gamma-rays by electrons accelerated at the reverse shock, in addition to the synchrotron self-Compton component. The predicted early optical flux could have been brighter than the naked-eye GRB 080319B. We also show that proton synchrotron emission from accelerated ultra-high-energy cosmic rays (UHECRs) is detectable, and could potentially explain $\gtrsim \rm TeV$ photons detected by LHAASO or UHECR acceleration can be constrained. Our model suggests that the detection of $\mathcal{O}(10\rm~TeV)$ photons with energy up to $\sim18$ TeV is possible for reasonable models of the extragalactic background light without invoking new physics, and predicts anti-correlations between MeV photons and TeV photons, which can be tested with the LHAASO data., Comment: 9 pages, 4 figures, accepted for publication in ApJL. Results and figures are updated, and the main conclusions remain unchanged

The joint detection of gravitational waves and the gamma-ray counterpart of a binary neutron star merger event, GW170817, unambiguously validates the connection between short gamma-ray bursts and compact binary object (CBO) mergers. We focus on a special scenario where short gamma-ray bursts produced by CBO mergers are embedded in disks of active galactic nuclei (AGN), and we investigate the $\gamma$-ray emission produced in the internal dissipation region via synchrotron, synchrotron self-Compton and external inverse-Compton (EIC) processes. In this scenario, isotropic thermal photons from the AGN disks contribute to the EIC component. We show that a low-density cavity can be formed in the migration traps, leading to the embedded mergers producing successful GRB jets. We find that the EIC component would dominate the GeV emission for typical CBO mergers with an isotropic-equivalent luminosity of $L_{j,\rm iso}=10^{48.5}~\rm erg~s^{-1}$ which are located close to the central supermassive black hole. Considering a long-lasting jet of duration $T_{\rm dur}\sim10^2-10^3$ s, we find that the future CTA will be able to detect its $25-100$ GeV emission out to a redshift $z=1.0$. {In the optimistic case, it is possible to detect the on-axis extended emission simultaneously with GWs within one decade using MAGIC, H.E.S.S., VERITAS, CTA, and LHAASO-WCDA. Early diagnosis of prompt emissions with Fermi-GBM and HAWC can provide valuable directional information for the follow-up observations., Comment: 11 pages, 5 figures. Accepted for publication in The Astrophysical Journal

As a powerful source of gravitational waves (GW), a supermassive black hole (SMBH) merger may be accompanied by a relativistic jet that leads to detectable electromagnetic (EM) emission. We model the propagation of post-merger jets inside a pre-merger circumnuclear environment formed by disk winds, and calculate multi-wavelength EM spectra from the forward shock region. We show that the non-thermal EM signals from SMBH mergers are detectable up to the detection horizon of future GW facilities such as the Laser Interferometer Space Antenna (LISA). Calculations based on our model predict slowly fading transients with time delays from days to months after the coalescence, leading to implications for EM follow-up observations after the GW detection., Comment: 6 pages, 3 figures, accepted for publication in ApJL. Schematic diagrams added

The recent detection of TeV photons from two gamma-ray bursts (GRBs), GRB 190114C and GRB 180720B, has opened a new window for multi-messenger and multi-wavelength astrophysics of high-energy transients. We study the origin of very-high-energy (VHE) $\gamma$-rays from the short GRB 160821B, for which the MAGIC Collaboration reported a $\sim 3 \sigma$ statistical significance. Short GRBs are often accompanied by extended and plateau emission, which is attributed to internal dissipation resulting from activities of a long-lasting central engine, and Murase et al. (2018) recently suggested the external inverse-Compton (EIC) scenario for VHE counterparts of short GRBs and neutron star mergers. Applying this scenario to GRB 160821B, we show that the EIC flux can reach $\sim 10^{-12}\rm~erg~cm^{-2}~s^{-1}$ within a time period of $\sim 10^3 - 10^4\rm~s$, which is consistent with the MAGIC observations. EIC $\gamma$-rays expected during the extended and plateau emission will be detectable with greater significance by future detectors such as the Cherenkov Telescope Array (CTA). The resulting light curve has a distinguishable feature, where the VHE emission is predicted to reach the peak around the end of the seed emission., Comment: 6 pages, 2 figures, accepted for publication in ApJL

Supermassive black hole (SMBH) coalescences are ubiquitous in the history of the Universe and often exhibit strong accretion activities and powerful jets. These SMBH mergers are also promising candidates for future gravitational wave detectors such as Laser Space Inteferometric Antenna (LISA). In this work, we consider neutrino counterpart emission originating from the jet-induced shocks. The physical picture is that relativistic jets launched after the merger will push forward inside the premerger disk wind material, and then they subsequently get collimated, leading to the formation of internal shocks, collimation shocks, forward shocks and reverse shocks. Cosmic rays can be accelerated in these sites and neutrinos are expected via the photomeson production process. We formulate the jet structures and relevant interactions therein, and then evaluate neutrino emission from each shock site. We find that month-to-year high-energy neutrino emission from the postmerger jet after the gravitational wave event is detectable by IceCube-Gen2 within approximately five to ten years of operation in optimistic cases where the cosmic-ray loading is sufficiently high and a mildly super-Eddington accretion is achieved. We also estimate the contribution of SMBH mergers to the diffuse neutrino intensity, and find that a significant fraction of the observed very high-energy ($E_\nu\gtrsim1$ PeV) IceCube neutrinos could originate from them in the optimistic cases. In the future, such neutrino counterparts together with gravitational wave observations can be used in a multimessenger approach to elucidate in greater detail the evolution and the physical mechanism of SMBH mergers., Comment: 17 pages, 9 figures, 1 table. Accepted for publication in PRD

We investigate the blazar contribution to the cumulative neutrino intensity assuming a generic relationship between neutrino and gamma-ray luminosities, $L_{\nu}\propto (L_{\rm ph})^{\gamma_{\rm lw}}$. Using the gamma-ray luminosity functions for blazars including flat spectrum radio quasars (FSRQs) and BL Lac objects, as well as the $Fermi$-LAT detection efficiency, we estimate contributions from {blazars resolved by $Fermi$-LAT as well as the unresolved counterpart.} Combining the existing upper limits from stacking analyses, the cumulative neutrino flux from all blazars ({including $Fermi$-LAT resolved and unresolved ones}) are constrained { in the range $0\lesssim\gamma_{\rm lw}\lesssim2.5$}. We also evaluate the effects of the redshift evolution and the effective local number densities for each class of FSRQs, BL Lacs, and all blazars, by which we place another type of constraints on the blazar contribution using the non-detection of high-energy neutrino multiplets. We demonstrate that these two upper limits are complementary, and that the joint consideration of the stacking and multiplet analyses { not only supports the argument that blazars are disfavored as the dominant sources of the 100-TeV neutrino background, but it extends this argument by including also $Fermi$-LAT-unresolved blazars as well, for a more generic luminosity correlation $L_{\nu}\propto (L_{\rm ph})^{\gamma_{\rm lw}}$}., Comment: Accepted for publication in The Astrophysical Journal. 8 pages, 4 figures

Shocks arising in galaxy mergers could accelerate cosmic-ray (CR) ions to TeV-PeV energies. While propagating in the intergalactic medium, these CRs can produce high-energy neutrinos, electron-positron pairs and gamma-rays. In the presence of intergalactic magnetic fields, the secondary pairs will radiate observable emissions through synchrotron radiation and inverse Compton scattering. In this paper, we demonstrate that these emissions can explain the radio and X-ray fluxes of merging galaxies such as NGC 660 and NGC 3256. Using our model in combination with the observations, we can constrain the gas mass, shock velocity, magnetic field and the CR spectral index $s$ of these systems. For NGC 660 a single-zone model with a spectral index $2.1\lesssim s\lesssim2.2$ is able to reproduce simultaneously the radio and X-ray observations, while a simple one-zone scenario with $s\sim2$ can describe the radio and a large fraction of X-ray observations of NGC 3256. Our work provides a useful approach for studying the dynamics and physical parameters of galaxy mergers, which can play an important part in future multi-messenger studies of similar and related extragalactic sources., Comment: 11 pages, 6 figures. Accepted by ApJ

The merger of dark matter halos and the gaseous structures embedded in them, such as proto-galaxies, galaxies, and groups and clusters of galaxies, results in strong shocks that are capable of accelerating cosmic rays (CRs) to $\sim10~\rm PeV$. These shocks will produce high-energy neutrinos and $\gamma$-rays through inelastic $pp$ collisions with ambient gaseous environments. In this work, we study the contributions of these halo mergers to the diffuse neutrino flux measured in IceCube and to the non-blazar portion of the extragalactic $\gamma$-ray background measured by $Fermi$. In order to calculate them, we formulate the redshift dependence of the shock velocity, galactic radius, halo gas content and galactic/intergalactic magnetic fields over the dark matter halo distribution up to a redshift $z=10$. We find that high-redshift mergers contribute a significant amount of the cosmic-ray energy luminosity density, and the resulting neutrino spectra could explain a large part of the observed diffuse neutrino flux above 0.1 PeV up to several PeV. We also show that our model can somewhat alleviate tensions with the extragalactic $\gamma$-ray background. First, since a larger fraction of the CR energy luminosity density comes from high redshifts, the accompanying $\gamma$-rays are more strongly suppressed through $\gamma\gamma$ annihilations with the cosmic microwave background (CMB) and the extragalactic background light (EBL). Second, mildly radiative-cooled shocks may lead to a harder CR spectrum with spectral indices of $1.5\lesssim s\lesssim2.0$. Our study suggests that halo mergers, a fraction of which may also induce starbursts in the merged galaxies, can be promising neutrino emitters without violating the existing $Fermi$ $\gamma$-ray constraints on the non-blazar component of the extragalactic $\gamma$-ray background., Comment: 13 pages, 7 figures, ApJ accepted