LeHaMoC: a versatile time-dependent lepto-hadronic modeling code for high-energy astrophysical sources

Abstract

Recent associations of high-energy neutrinos with active galactic nuclei (AGN) have revived the interest in leptohadronic models of radiation from astrophysical sources. The rapid increase in the amount of acquired multi-messenger data will require soon fast numerical models that may be applied to large source samples. We develop a time-dependent leptohadronic code, LeHaMoC, that offers several notable benefits compared to other existing codes, such as versatility and speed. LeHaMoC solves the Fokker-Planck equations of photons and relativistic particles (i.e. electrons, positrons, protons, and neutrinos) produced in a homogeneous magnetized source that may also be expanding. The code utilizes a fully implicit difference scheme that allows fast computation of steady-state and dynamically evolving physical problems. We first present test cases where we compare the numerical results obtained with LeHaMoC against exact analytical solutions and numerical results computed with ATHE$\nu$A, a well-tested code of similar philosophy but different numerical implementation. We find a good agreement (within 10-30%) with the numerical results obtained with ATHE$\nu$A without evidence of systematic differences. We then demonstrate the capabilities of the code through illustrative examples. First, we fit the spectral energy distribution from a jetted AGN in the context of a synchrotron-self Compton model and a proton-synchrotron model using Bayesian inference. Second, we compute the high-energy neutrino signal and the electromagnetic cascade induced by hadronic interactions in the corona of NGC 1068. LeHaMoC is easily customized to model a variety of high-energy astrophysical sources and has the potential to become a widely utilized tool in multi-messenger astrophysics.

19 pages, version accepted in A&A. The code is available at: https://github.com/mariapetro/LeHaMoC.git

References in corpus (46)

• emcee: The MCMC Hammer
• Unified Schemes for Radio-Loud Active Galactic Nuclei
• The Nuclear Spectroscopic Telescope Array (NuSTAR) Mission
• Evidence for High-Energy Extraterrestrial Neutrinos at the IceCube Detector
• Fermi Large Area Telescope Fourth Source Catalog
• Multi-messenger observations of a flaring blazar coincident with high-energy neutrino IceCube-170922A
• Neutrino emission from the direction of the blazar TXS 0506+056 prior to the IceCube-170922A alert
• Energy spectra of gamma-rays, electrons and neutrinos produced at proton-proton interactions in the very high energy regime
• Leptonic and Hadronic Modeling of Fermi-Detected Blazars
• The HIPASS Catalogue - I. Data Presentation
• TeV Gamma Rays from BL Lac Objects due to Synchrotron Radiation of Extremely High Energy Protons
• Evidence for neutrino emission from the nearby active galaxy NGC 1068
• Time-integrated Neutrino Source Searches with 10 years of IceCube Data
• A Proton Synchrotron Blazar Model for Flaring in Markarian~501
• The Challenge of Machine Learning in Space Weather Nowcasting and Forecasting
• Monte-Carlo simulations of photohadronic processes in astrophysics
• On the X-ray, optical emission line and black hole mass properties of local Seyfert galaxies
• X-CIGALE: fitting AGN/galaxy SEDs from X-ray to infrared
• Energy spectra of gamma-rays, electrons and neutrinos produced at interactions of relativistic protons with low energy radiation
• Fast Automated Analysis of Strong Gravitational Lenses with Convolutional Neural Networks
• A hadronic origin for ultra-high-frequency-peaked BL Lac objects
• Propagation of Ultra High Energy Protons over Cosmological Distances and Implications for Topological Defect Models
• Hidden Cores of Active Galactic Nuclei as the Origin of Medium-Energy Neutrinos: Critical Tests with the MeV Gamma-Ray Connection
• Photohadronic origin of $\gamma$-ray BL Lac emission: implications for IceCube neutrinos
• NuSTAR Spectroscopy of Multi-Component X-ray Reflection from NGC 1068
• SEDfit: Software for Spectral Energy Distribution Fitting of Photometric Data
• The time-dependent one-zone hadronic model - First principles
• Hidden Hearts of Neutrino Active Galaxies
• Progress in unveiling extreme particle acceleration in persistent astrophysical jets
• On the direct correlation between gamma-rays and PeV neutrinos from blazars
• Active galactic nuclei at gamma-ray energies
• Constraints on gamma-ray and neutrino emission from NGC 1068 with the MAGIC telescopes
• Seyfert 2 galaxies in the GeV band: jets and starburst
• 3HSP J095507.9+355101: a flaring extreme blazar coincident in space and time with IceCube-200107A
• Spectral and temporal signatures of ultrarelativistic protons in compact sources
• Time-dependent lepto-hadronic modeling of the emission from blazar jets with SOPRANO: the case of TXS 0506+056, 3HSP J095507.9+355101 and 3C 279
• Neutrino emission in the hadronic Synchrotron Mirror Model: the "orphan" TeV flare from 1ES 1959+650
• Comprehensive Multimessenger Modeling of the Extreme Blazar 3HSP J095507.9+355101 and Predictions for IceCube
• Temporal signatures of leptohadronic feedback mechanisms in compact sources
• The TeV emission of Ap Librae: a hadronic interpretation and prospects for CTA
• Improved photomeson model for interactions of cosmic ray nuclei
• Multi-Frequency Observations of the Candidate Neutrino Emitting Blazar BZB J0955+3551
• An expanding one-zone model for studying blazars emission
• A Bayesian approach for torque modelling of BeXRB pulsars with application to super-Eddington accretors
• The redshift and the host galaxy of the neutrino candidate 4FGL J0955.1+3551 (3HSP J095507.9+355101)
• A roadmap to hadronic supercriticalities: a comprehensive study of the parameter space for high-energy astrophysical sources

Cited by in corpus (8)

• AM${}^3$: An Open-Source Tool for Time-Dependent Lepto-Hadronic Modeling of Astrophysical Sources
• Astroparticles from X-ray Binary Coronae
• A closer look at the electromagnetic signatures of Bethe-Heitler pair production process in blazars
• Camelidae on BOAT: observation of a second spectral component in GRB 221009A
• Stochastic acceleration in Extreme TeV BL Lacs through MCMC
• Insights from leptohadronic modelling of the brightest blazar flare
• Impact of model parameter degeneracy on leptonic radiation models. The case of blazar multi-wavelength spectra
• Constraining the synchrotron peak and estimating the VHE brightness of a sample of extreme high synchrotron peak blazars