Garret Cotter

Abstract:

Context. Composite galaxies that contain both Seyfert and starburst components may produce very high-energy (VHE; > 100 GeV) γ -ray emission at a wide range of spatial scales, from a few Schwarzschild radii of a supermassive black hole (SMBH; R S = 10 −6 pc for M SMBH = 10 7 M ⊙ ) to dimensions of kiloparsec-size jet-driven outflows. In addition to supernova remnants, various sources have been suggested to explain data collected on composite galaxies, including multi-messenger neutrino and ultra-high-energy cosmic-ray data. Aims. The closest composite Seyfert–starburst galaxies (NGC 1068, the Circinus galaxy, and NGC 4945) are observed with the High Energy Stereoscopic System (H.E.S.S.) to provide constraints on cosmic-ray populations in these systems. Methods. Data obtained in H.E.S.S. observations have been analyzed to search for VHE γ -ray counterparts to the GeV γ -ray signals detected with Fermi -LAT and for potential spectral components in the VHE range. Results. No significant signals have been found in these H.E.S.S. data. Upper limits on the VHE γ -ray fluxes were applied to constrain theoretical models involving different spectral components.

Abstract:

The blazar PKS 0903 - 57 has recently been classified as a flat spectrum radio quasar at a redshift of z = 0.2621 . In March and April 2020, Fermi-LAT and AGILE reported tremendous activity in high-energy γ rays with the flux increasing by  ∼ 2 orders of magnitude compared to quiescence. The flare was observed with H.E.S.S. in very-high-energy γ rays for six nights with a total observation time of 13.1 h, resulting in the discovery of PKS 0903 - 57 in this energy band with an average flux of 1.5 × 10 − 10 ph cm − 2 s − 1 above an energy threshold of  ∼ 180 GeV corresponding to 60% of the Crab Nebula flux above the same threshold. The very-high-energy γ-ray flux was strongly variable. X-ray and optical data were collected with Swift and ATOM, and also indicate significant variability. The observed multiwavelength flux and spectral variability during the H.E.S.S. observation window suggest variability time scales on the order of a few hours and reveal complex correlation patterns. The lack of absorption beyond that of the extragalactic background light in the γ-ray domain suggests that the emission region was located outside of the broad-line region. A leptonic one-zone modeling of the six H.E.S.S. observation nights using the dusty torus as seed photons for the inverse-Compton scattering, results in a low magnetization of the emission region. This implies that shock acceleration is likely the main driver during the event.

Abstract:

DIPLODOCUS (Distribution-In-PLateaux methODOlogy for the CompUtation of transport equationS) is a framework being developed for the general transport of particle distribution functions through the seven dimensions of phase space, including forcing terms and interactions between particles. Following Paper I, which details the mathematical background, this second paper provides an overview of the numerical implementation in the form of the code package Diplodocus $.$ jl, written in Julia, including the description of a novel Monte-Carlo sampling technique for the pre-computation of anisotropic collision integrals. In addition to the discussion of numerical implementation, a selection of test cases are presented to examine the package’s capabilities. These test cases focus on micro-scale physical effects: binary collisions, emissive interactions and external forces that are relevant to the modelling of jetted astrophysical sources, such as Active Galactic Nuclei and X-Ray Binaries.