DIPLODOCUS II: Implementation of transport equations and test cases relevant to micro-scale physics of jetted astrophysical sources

(2026)

Authors:

Christopher N Everett, Marc Klinger-Plaisier, Garret Cotter

Black Holes as Telescopes: Discovering Supermassive Binaries through Quasiperiodic Lensed Starlight

Physical Review Letters American Physical Society (APS) 136:6 (2026) 061403

Authors:

Hanxi Wang, Miguel Zumalacárregui, Bence Kocsis

Abstract:

Supermassive black hole (SMBH) binary systems are an unavoidable outcome of galaxy mergers. Their dynamics encode valuable information about their formation and growth, the composition of their host galactic nuclei, the evolution of galaxies, and the nature of gravity. Many SMBH binaries with separations pc-kpc have been found, but closer (subparsec) binaries remain to be confirmed. Identifying these systems may elucidate how binaries evolve past the “final parsec” until gravitational radiation drives them to coalescence. Methods to discover and characterize SMBH binaries can shed light on these important questions and potentially open new multimessenger channels. Here we show that SMBH binaries in nonactive galactic nuclei can be identified and characterized by the gravitational lensing of individual bright stars, located behind them in the host galaxy. The rotation of “caustics”—regions where sources are hugely magnified due to the SMBH binary’s orbit and inspiral—leads to quasiperiodic lensing of starlight (QPLS). The extreme lensing magnification of individual bright stars produces a significant variation in the host galaxies’ luminosity; their lightcurve traces the orbit of the SMBH binary and its evolution, analogous to the waveforms recorded by gravitational-wave (GW) detectors. QPLS probes the population of sources observable by pulsar timing arrays and space detectors (LISA, TianQin), offering advance warning triggers for merging SMBHs for coincident or follow-up GW detections. SMBH population models predict 1–50 [ 190 5000 ] ( n / pc 3 ) QPLS binaries with period less than 10[40] yr with comparable masses and redshift z < 0.3 , where n is the stellar number density. Additionally, stellar and orbital motion will lead to frequent instances of single or double flares caused by SMBHBs with longer periods. This novel signature can be searched for in a wealth of existing and upcoming time-domain photometric data: identifying quasiperiodic variability in galactic lightcurves will reveal an ensemble of binary systems and illuminate outstanding questions around them.

Detection of an extremely luminous radio counterpart to the Be/X-ray binary A0538−66

Monthly Notices of the Royal Astronomical Society Oxford University Press 548:1 (2026) stag224

Authors:

Justine Crook-Mansour, Rob Fender, Alex Andersson, Hao Qiu, Andrew K Hughes, Jakob van den Eijnden, Fraser J Cowie, Sara Motta, Itumeleng Monageng, Lorenzo Ducci, Sandro Mereghetti, Andries Mathiba, Dougal Dobie, Tara Murphy, David L Kaplan, Francesco Carotenuto, Phil Charles

Abstract:

We present the discovery of radio emission from the Be/X-ray binary A0538−66 with the Australian Square Kilometre Array Pathfinder, and results from a subsequent weekly monitoring campaign with the MeerKAT radio telescope. A0538−66, located in the Large Magellanic Cloud, hosts a neutron star with a short spin period ( ms) in a highly eccentric -d orbit . Its rare episodes of super-Eddington accretion, rapid optical and X-ray flares, and other peculiar properties make it an interesting system among high-mass X-ray binaries. Our MeerKAT data reveal that it is also one of the most radio-luminous neutron star X-ray binaries observed to date, reaching (at 1.28 GHz), with radio emission that appears to be orbitally modulated. We consider several possible mechanisms for the radio emission, and place A0538−66 in context by comparing it to similar systems.

Constraining the nature of the most extreme Galactic particle accelerator

Astronomy & Astrophysics EDP Sciences 706 (2026) a8

Authors:

A Acharyya, F Aharonian, H Ashkar, M Backes, R Batzofin, D Berge, K Bernlöhr, M Böttcher, C Boisson, J Bolmont, F Brun, B Bruno, C Burger-Scheidlin, T Bylund, S Casanova, J Celic, M Cerruti, A Chen, M Chernyakova, JO Chibueze, O Chibueze, B Cornejo, G Cotter, J de Assis Scarpin, M de Bony de Lavergne, M de Naurois, E de Oña Wilhelmi, AG Delgado Giler, J Devin, A Djannati-Ataï, A Dmytriiev, K Egberts, K Egg, J-P Ernenwein, C Escañuela Nieves, P Fauverge, K Feijen, MD Filipovic, G Fontaine, S Funk, S Gabici, YA Gallant, JF Glicenstein, J Glombitza, P Goswami, M-H Grondin, L Heckmann, B Heß, JA Hinton, W Hofmann, TL Holch, M Holler, M Jamrozy, F Jankowsky, A Jardin-Blicq, I Jaroschewski, D Jimeno, I Jung-Richardt, K Katarzyński, D Kerszberg, B Khélifi, N Komin, K Kosack, D Kostunin, RG Lang, S Lazarević, A Lemière, M Lemoine-Goumard, J-P Lenain, P Liniewicz, A Luashvili, J Mackey, D Malyshev, V Marandon, MGF Mayer, A Mehta, AMW Mitchell, R Moderski, L Mohrmann, A Montanari, E Moulin, J Niemiec, L Olivera-Nieto, MO Moghadam, S Panny, RD Parsons, U Pensec, P Pichard, T Preis, G Pühlhofer, M Punch, A Quirrenbach, A Reimer, O Reimer, I Reis, Q Remy, HX Ren, B Reville, F Rieger, G Roellinghoff, G Rowell, B Rudak, K Sabri, S Safi-Harb, V Sahakian, A Santangelo, M Sasaki, F Schüssler, JNS Shapopi, W Si Said, H Sol, Ł Stawarz, S Steinmassl, T Tanaka, AM Taylor, GL Taylor, R Terrier, Y Tian, A Timmermans, M Tsirou, N Tsuji, T Unbehaun, C van Eldik, M Vecchi, C Venter, J Vink, V Voitsekhovskyi, SJ Wagner, A Wierzcholska, M Zacharias, AA Zdziarski, A Zech, W Zhong, S Takekawa

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.

H.E.S.S. detection and multi-wavelength study of the z  ∼  1 blazar PKS 0346−27

Astronomy & Astrophysics EDP Sciences 706 (2026) a246

Authors:

F Aharonian, M Backes, R Batzofin, Y Becherini, D Berge, K Bernlöhr, B Bi, M Böttcher, C Boisson, J Bolmont, F Brun, B Bruno, C Burger-Scheidlin, D Cecchin Momesso, J Celic, M Cerruti, A Chen, M Chernyakova, JO Chibueze, O Chibueze, B Cornejo, G Cotter, J Damascene Mbarubucyeye, ID Davids, J de Assis Scarpin, M de Bony de Lavergne, M de Naurois, E de Oña Wilhelmi, AG Delgado Giler, J Devin, A Djannati-Ataï, A Dmytriiev, K Egberts, K Egg, J-P Ernenwein, C Esca nuela Nieves, K Feijen, MD Filipovic, G Fontaine, S Funk, S Gabici, M Genaro, JF Glicenstein, J Glombitza, P Goswami, L Heckmann, B Hess, JA Hinton, W Hofmann, TL Holch, M Holler, D Horns, M Jamrozy, F Jankowsky, I Jaroschewski, I Jung-Richardt, K Kasprzak, K Katarzyński, D Kerszberg, B Khélifi, N Komin, K Kosack, D Kostunin, RG Lang, S Lazarević, A Lemière, J-P Lenain, P Liniewicz, A Luashvili, J Mackey, D Malyshev, V Marandon, M Mayer, A Mehta, AMW Mitchell, R Moderski, L Mohrmann, A Montanari, E Moulin, J Niemiec, MO Moghadam, S Panny, RD Parsons, U Pensec, P Pichard, T Preis, G Pühlhofer, M Punch, A Quirrenbach, A Reimer, O Reimer, I Reis, B Rudak, K Sabri, V Sahakian, D Jimeno, A Santangelo, M Sasaki, F Schüssler, JNS Shapopi, W Si Said, Ł Stawarz, S Steinmassl, T Takahashi, T Tanaka, AM Taylor, GL Taylor, R Terrier, T Unbehaun, C van Eldik, M Vecchi, C Venter, J Vink, T Wach, SJ Wagner, A Wierzcholska, M Zacharias, A Zech, W Zhong

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.