Garret Cotter

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.

Abstract:

Context. Microquasars have emerged as promising candidates to explain the cosmic-ray flux at petaelectronvolt energies. LHAASO observations revealed V4641 Sgr as the most extreme example so far. Its gamma-ray spectrum extends up to 800 TeV, which requires particles with multi-PeV energy. The TeV emission is highly extended, which challenges expectations given the reported low-inclination angle of the V4641 Sgr jets. Aims. We spatially and spectrally resolved the gamma-ray emission from V4641 Sgr and investigated the particle acceleration in the system. Methods. Using ≈100 h of H.E.S.S. data, we performed a spectro-morphological study of the gamma-ray emission around V4641 Sgr. We employed HI and dedicated CO observations of the region to infer the target material for cosmic-ray interactions. Results. We detected multi-TeV emission around V4641 Sgr with a high significance. The emission region is elongated, and its major and minor axes are 0.34° ±0.01 syst ±0.04 stat and 0.06° ±0.01 syst ±0.01 stat , respectively. We found a power-law spectrum with an index ≈1.8, and together with results from other gamma-ray instruments, this reveals a spectral energy distribution (SED) that peaks at energies of ≈100 TeV for the first time. We found indications (3 σ ) of a two-component morphology, with indistinguishable spectral properties. The position of V4641 Sgr is inconsistent with the best-fit position of the single-component model and with the dip between the two components. We found no significant evidence of an energy-dependent morphology. No dense gas was found at any distance towards V4641 Sgr, which places an upper limit of n gas ≲ 0.2 cm −3 within the gamma-ray emission region. Conclusions. The peak of the SED at ≈100 TeV identifies V4641 Sgr as a candidate cosmic-ray accelerator beyond the so-called knee. The absence of dense target gas places stringent energetic constraints on hadronic interpretations, however. The H.E.S.S. measurement requires an unusually hard (≈1.5) spectral index for the protons. A leptonic scenario faces fewer obstacles if the particle transport is fast enough to avoid losses and to reproduce the observed energy-independent morphology. The absence of bright X-ray emission across the gamma-ray emission region requires a magnetic field strength ≲3 μG, however. Our findings favour a leptonic origin of the gamma-ray emission. This conclusion does not exclude hadron acceleration in the V4641 Sgr system.

Abstract:

Context. PKS 0346-27 is a low synchrotron peaked blazar at redshift 0.991. The very high energy (VHE; E > 100 GeV) spectra of blazars are always affected by γγ absorption by the extragalactic background light (EBL), and subsequently no blazars have been detected in VHE γ -rays at redshifts exceeding 1. Aims. This is the goal of a target-of-opportunity (ToO) programme by H.E.S.S.: to observe flaring high-redshift ( z ≳ 1) blazars. Importantly, extending the redshift range of VHE-detected blazars to z ≳ 1 will yield insights into the cosmological evolution of both the VHE blazar population and the EBL. Methods. We report H.E.S.S. ToO and multi-wavelength observations of the blazar PKS 0346−27. We analysed and modelled the H.E.S.S. data together with simultaneous data from Fermi -LAT, Swift (XRT and UVOT), using single-zone leptonic and hadronic models. Results. PKS 0346-27 was detected by H.E.S.S. at a significance of 6.3 σ during one night on 3 November 2021, while for other nights before and after this day, upper limits on the VHE flux have been determined. No evidence for intra-night γ -ray variability has been found. A flare in high-energy ( E > 100 MeV) γ -rays detected by Fermi -LAT preceded the H.E.S.S. detection by 2 days. A fit with a single-zone emission model to the contemporaneous spectral energy distribution during the detection night was possible with a proton-synchrotron-dominated hadronic model, requiring a proton-kinetic-energy-dominated jet power temporarily exceeding the source’s Eddington limit, although alternative (e.g. multi-zone) models cannot be ruled out. A one-zone leptonic model is, in principle, also able to fit the flare-state spectral energy distribution. However, it requires implausible parameter choices, in particular, extreme Doppler and bulk Lorentz factors of ≳80.

Abstract:

DIPLODOCUS (Distribution-In-PLateaux methODOlogy for the CompUtation of transport equationS) is a novel framework being developed for the mesoscopic modelling of astrophysical systems via the transport of particle distribution functions through the seven dimensions of phase space, including continuous forces and discrete interactions between particles. This first paper in a series provides an overview of the analytical framework behind the model, consisting of an integral formulation of the relativistic transport equations (Boltzmann equations) and a discretisation procedure for the particle distribution function (Distribution-In-Plateaux). The latter allows for the evaluation of anisotropic interactions, and generates a conservative numerical scheme for a distribution function’s transport through phase space.