Peer Revenue Publish or profit

Papers

Fresh from arXiv, typeset properly. Highlighted words trade on the market.

Physics Astrophysics High Energy Astrophysical Phenomena

astro-ph.HE Watch High Energy Astrophysical Phenomena (sign in to use)

Friday 18 September

High Energy Astrophysical Phenomena

Geometrically Thin Sub-Eddington AGN Accretion Disks show a suppressed Lyman Edge

The observed ultraviolet continua of active galactic nuclei (AGN) generally lack a strong intrinsic HI Lyman edge predicted by classical optically thick accretion disk atmosphere models. We revisit this long-standing problem using our previous sub-Eddington ($L/L_{\rm Edd}\sim 0.03$) geometrically thin 3D accretion disk simulation evolved with gray angle-dependent radiation magnetohydrodynamics (MHD) around a $10^8M_{\odot}$ black hole. This disk is magnetic pressure dominated, and has surface densities more than two orders of magnitude below the corresponding radiation-pressure-supported $\alpha$-disk prediction. By restarting the simulation with multiple frequency groups, we compute the emergent continuum directly from the frequency-dependent radiation fluxes and find no sharp HI Lyman edge at 13.6 eV. This suppressed Lyman edge is mainly caused by the low densities in a magnetic pressure supported disk where the Lyman opacity jump is far weaker than in standard disk models. These results indicate that magnetic pressure support may provide a possible solution to reconciling optically thick AGN disks with the smooth observed continuum emission around the Lyman edge. However, in this simulation, a substantial fraction of the emission near 13.6 eV is expected to originate at smaller radii closer to the central black hole, which is outside our simulation domain. More detailed predictions of the spectrum will require an extension to the inner disk, non-LTE radiation transfer, higher frequency resolution and a survey across a broader range of accretion rates.
Read the full abstract

All authors Ish Kaul, Yan-Fei Jiang, Omer Blaes, Lizhong Zhang

High Energy Astrophysical Phenomena · Cross-listed from astro-ph.GA

The AGN Channel in 3D: Scattering Belts and the Importance of Eccentricity in the Black Hole Population

Active galactic nuclei (AGN) are a promising origin for observed gravitational wave mergers. Current population synthesis models are limited to 1D N-body or Monte Carlo methods which rely on statistical approaches to resolving dynamical scatterings. We present three-dimensional hybrid $N$-body simulations of a population of black holes (BHs) surrounding an AGN using a new code in development, AGNBI, where close interactions are directly simulated for both single and binary objects. Our results show binary formations occur throughout the AGN disk, particularly where the migration time is long. While BHs efficiently migrate and pair up in migration traps, they are frequently disrupted by ionisation or exchanges from binary-single and binary-binary interactions. We report a broadened radial distribution of formations and mergers about the trap radius when accounting for binary-single and binary-binary interactions, which we describe as a "scattering belt”. We assess the influence of the pre-existing BH population. When we model a non-zero eccentricity distribution in the BH population, the binary formations and mergers are suppressed by an order of magnitude compared to a non-eccentric population as fewer BHs embed in the AGN disk. We report an approximate merger rate of $\mathcal{R}_\mathrm{GW}\approx2\text{--}12\,\mathrm{Gpc}^{-3}\,\mathrm{yr}^{-1}$ with maximum BH masses of $10^2\text{--}10^3\,M_\odot$ through hierarchical mergers, where lower initial velocity dispersions produce larger BHs. The mass and mass ratio distributions are comparatively flat, suggesting AGN may be more relevant for high-mass or unequal-mass gravitational wave sources.
Read the full abstract

High Energy Astrophysical Phenomena

Euclid: Quick Data Release (Q1) – Optical hotspots in powerful radio galaxies

The hotspots of powerful radio galaxies are the sites of high-energy particle acceleration, but radio observations give us limited information on the electron energy distributions that they generate, which affects our overall understanding of their physics and energetics. Optical synchrotron counterparts to radio galaxy hotspots, which probe high electron energies, have been known for many years, but have only been studied in a small number of objects selected as being very bright in the radio. By combining Low-Frequency Array (LOFAR) observations with the Euclid Q1 data release in the Euclid Deep Field North, we identify 77 powerful radio galaxies with hotspots and show that a small but significant fraction of them, 7–15%, may have detected infrared and optical counterparts to their hotspots. The hotspot spectral energy distributions in our clearest cases are consistent with particle acceleration up to (but not beyond) 1\,TeV in energy on the assumption of field strengths close to equipartition, whereas non-detected hotspots either do not reach such high electron energies or would have to have very different electron energy spectra or significantly sub-equipartition magnetic field strengths. Two new candidate optical hotspots are clearly detected matching high-resolution LOFAR data; one is remarkable as both the highest-redshift optical hotspot known to date ($z \approx 1.6$) and also the first optical detection of the hotspot of the inner lobe of a restarting source. We discuss approaches to performing a large-scale study of hotspots with the full wide-area Euclid surveys, particularly when combined with the international baselines of LOFAR, enabling for the first time a statistical analysis of high-energy particle acceleration in hotspots.
Read the full abstract

All authors M. J. Hardcastle, L. Bisigello, M. Bondi, M. Magliocchetti, L. K. Morabito, J. Petley, F. Ricci, H. J. A. Rottgering, K. Rubinur, D. J. B. Smith, D. Stern, B. Altieri, S. Andreon, N. Auricchio, C. Baccigalupi, M. Baldi, A. Balestra, S. Bardelli, P. Battaglia, A. Biviano, E. Branchini, M. Brescia, S. Camera, V. Capobianco, C. Carbone, J. Carretero, S. Casas, M. Castellano, G. Castignani, S. Cavuoti, K. C. Chambers, A. Cimatti, C. Colodro-Conde, G. Congedo, C. J. Conselice, L. Conversi, Y. Copin, A. Costille, F. Courbin, H. M. Courtois, M. Cropper, J.-G. Cuby, A. Da Silva, H. Degaudenzi, G. De Lucia, C. Dolding, H. Dole, F. Dubath, X. Dupac, M. Farina, R. Farinelli, S. Ferriol, P. Fosalba, S. Fotopoulou, M. Frailis, E. Franceschi, M. Fumana, S. Galeotta, K. George, B. Gillis, C. Giocoli, J. Gracia-Carpio, A. Grazian, F. Grupp, S. V. H. Haugan, J. Hoar, W. Holmes, I. M. Hook, F. Hormuth, A. Hornstrup, K. Jahnke, M. Jhabvala, B. Joachimi, S. Kermiche, A. Kiessling, B. Kubik, M. Kümmel, M. Kunz, H. Kurki-Suonio, A. M. C. Le Brun, S. Ligori, P. B. Lilje, V. Lindholm, I. Lloro, G. Mainetti, O. Mansutti, O. Marggraf, M. Martinelli, N. Martinet, F. Marulli, R. J. Massey, E. Medinaceli, S. Mei, M. Meneghetti, E. Merlin, G. Meylan, A. Mora, M. Moresco, L. Moscardini, R. Nakajima, C. Neissner, S.-M. Niemi, J. W. Nightingale, C. Padilla, S. Paltani, F. Pasian, K. Pedersen, W. J. Percival, V. Pettorino, S. Pires, G. Polenta, L. A. Popa, F. Raison, A. Renzi, J. Rhodes, G. Riccio, E. Romelli, M. Roncarelli, B. Rusholme, R. Saglia, Z. Sakr, D. Sapone, B. Sartoris, M. Schirmer, P. Schneider, A. Secroun, G. Seidel, E. Sihvola, P. Simon, C. Sirignano, G. Sirri, J. Skottfelt, L. Stanco, P. Tallada-Crespí, A. N. Taylor, H. I. Teplitz, I. Tereno, N. Tessore, S. Toft, R. Toledo-Moreo, F. Torradeflot, I. Tutusaus, J. Valiviita, T. Vassallo, G. Verdoes Kleijn, A. Veropalumbo, Y. Wang, J. Weller, O. R. Williams, A. Zacchei, G. Zamorani, F. M. Zerbi, E. Zucca, J. García-Bellido, V. Scottez, F. Tarsitano, J. H. Knapen

High Energy Astrophysical Phenomena

Nonlinearly Causal General-Relativistic Two-Fluid Dissipative Magnetohydrodynamics

We present a formulation of general-relativistic (GR) 19-moment dissipative magnetohydrodynamics (MHD), capable of handling all first-order dissipative terms (viscosity, heat conductivity, resistivity and Hall terms), as well as all ideal electron degrees of freedom (number density, momentum and energy) in a GR astrophysical two-fluid plasma. We derive necessary and sufficient nonlinear causality conditions for the constrained first-order 19-moment system, as well as both necessary and sufficient conditions for strong hyperbolicity. To assess the formulation numerically, we develop a high-resolution shock-capturing scheme that solves these equations in a performance-portable fashion. The scheme expresses all evolution equations, including the dissipative and electron sectors, in flux-divergence form, allowing us to model systems with scales separated by orders of magnitude without resolving kinetic scales everywhere on the grid. In addition, we use implicit integration methods to systematically overstep kinetic scales in MHD regions, such as cyclotron and plasma frequencies. To make the scheme as robust as GRMHD codes, we construct necessary and sufficient conditions for a conserved state to be physically admissible, and based on these construct a new physicality-enforcement scheme. As an exact validation comparison, we formulate and derive a full solution to dissipative two-fluid Bondi accretion in general relativity. We then validate the equations against a series of results from kinetic particle-in-cell models of black hole accretion and magnetospheric dynamics, demonstrating that our formulation and scheme can correctly capture major features relevant for feedback on global scales of these solutions, including dimensionless reconnection rates of order $0.1$ and Braginskii-like anisotropic pressures and heat fluxes.
Read the full abstract

High Energy Astrophysical Phenomena · Cross-listed from gr-qc

Structure and astrophysical properties of dark energy admixed neutron stars

The presence of a dark energy (DE) core, made of a Chaplygin dark fluid, is capable of influencing the diverse observable properties of neutron stars (NSs). However, there is also the possibility that NSs could contain a mixture of ordinary nuclear matter and DE. In this perspective, we propose an equation of state (EoS) that includes both normal matter and DE, where a free parameter $f$ quantifies the ratio of DE density with respect to the total energy density. As expected, in the limit $f \rightarrow 0$ we recover the already known results corresponding to pure nuclear matter. Remarkably, increasing this fraction $f$ leads to a softer EoS and thus strongly favors the fulfillment of the causality condition. As a consequence of our mass–radius results, the central compact object in HESS J1731-347 can be satisfactorily described as a DE admixed NS. For sufficiently high stellar masses, our calculations further indicate that both the gravitational redshift and the fundamental nonradial oscillation frequency increase significantly with an increasing DE fraction, whereas the opposite trend is observed for the total gravitational mass and tidal deformability. As is typical for compact stars composed purely of ordinary matter, we find that in our DE admixed stellar models the central density corresponding to the maximum-mass configuration coincides with the point at which the squared radial vibration frequency vanishes. Moreover, higher DE concentrations lead to enhanced stability of these configurations.
Read the full abstract

All authors Juan M. Z. Pretel, Sergio B. Duarte, Mariana Dutra, Odilon Lourenço

High Energy Astrophysical Phenomena · Cross-listed from astro-ph.SR

SN 2026dix: a nearby, transitional Type IIb/Ib supernova consistent with a warm supergiant progenitor

Understanding the final evolution stages of (very) massive stars and their explosive outcomes, stripped-envelope supernovae (SESNe), represent a long-term challenge for astrophysics. The latest results support a continuous distribution within the traditional SESN subclasses (IIb, Ib, Ic). The nearby ($D \approx 17.5$ Mpc) SN 2026dix seems to be another member of the recently identified group of transitional Type IIb/Ib explosions. A point source is located at the SN position within the uncertainties on multiple pre-explosion HST images. Together with post-explosion photometry and spectroscopy, it provides a good opportunity to study a rare type of SN in detail. We carried out a thorough comparative light-curve (LC) and spectral analysis of SN 2026dix. We also constructed its bolometric LC and modeled it semi-analytically. In addition, we constructed the spectral energy distribution of the presumed progenitor and compared this to model stellar atmospheres and binary stellar evolution tracks. We infer that SN 2026dix arose from an explosion in an interacting binary system, consistent with the properties of a warm ($T_\textrm{eff} \approx 6750-$7750 K; spectral type F2 to A7), luminous ($\log(L_{\rm bol}/L_{\odot})\approx 4.9-5.3$) supergiant primary and a less luminous, less massive hot dwarf companion. The nature of the identified progenitor closely resembles that of some other known cases of SNe IIb and also well aligns with both the results of our semi-analytical LC modeling and the comparison of the LCs and spectra of SN 2026dix with the output of radiative-transfer models of an exploding star with $M_\textrm{ini}$ = 18 $M_{\odot}$.
Read the full abstract

All authors Tamás Szalai, Schuyler D. Van Dyk, Dominik Bánhidi, Andrea P. Nagy, Réka Könyves-Tóth, Alexei V. Filippenko, Thomas G. Brink, WeiKang Zheng, Imre Barna Bíró, István Csányi, Zsófia Bora, Zsóka Horváth, Ágoston Horti-Dávid, András. P. Joó, Dávid Koller, Levente Kriskovics, Klára Lelkes, Katalin Nagy, Krisztián Sárneczky, Ádám Sódor, Róbert Szakáts, András Pál, József Vinkó

High Energy Astrophysical Phenomena

Core-collapse supernovae 2022prr, 2023ucy, and 2024ljc in the luminous infrared galaxy NGC 6745

The core-collapse supernovae (CCSNe) 2022prr, 2023ucy, and 2024ljc exploded within the central regions of the nearby luminous infrared galaxy (LIRG) NGC6745. We present our follow-up data with the aim to characterise the spectrophotometric evolution of these CCSNe, determine their host-galaxy extinction, and to define their physical nature. The presented spectrophotometric data set covers the optical and near-infrared region for our three supernovae (SNe). The data are analysed and compared to other hydrogen-rich CCSNe, and a host-galaxy extinction value Av is derived based on broadband photometry for each of the targets. A spectral energy distribution (SED) for the host galaxy is constructed, and a combination of different galactic components are fitted to the SED to estimate the intrinsic CCSN rate and the age of the starburst of NGC6745. 2022prr is a SN2009ip-like event, whose light curve evolution consists of two consecutive luminous events A and B with peak magnitudes of Mr = -14.8 and -18.3 mag, respectively. The high degree of V and R-band polarisation (P=2.8%) near maximum light suggests interaction with highly aspherical circumstellar material. The spectral evolution of the SN is dominated by narrow Balmer emission lines, which have a complex evolution over time. 2022prr joins the growing sample of SN2009ip-like events that show early `flash ionisation' features of CIII, NIII, and HeII. 2023ucy is found to be a normal Type IIP SN. However, 2024ljc is a rapidly evolving Type IIb SN with a post-maximum r-band decline of roughly 3.2 mag by +41 d whereas normal IIb SNe decline typically by 1.5+/-0.3 mag during the same timescale. While the events exploded in a LIRG, the SNe are obscured by only moderate V-band host extinctions ranging from 0 to 1 mag. Based on the SED modelling we estimate an age of the host starburst to be 40Myr, and the intrinsic CCSN rate to be 0.4+/-0.1 SN/year.
Read the full abstract

All authors K. K. Matilainen, E. Kankare, T. Nagao, T. Reynolds, A. Efstathiou, S. Mattila, C. R. Angus, A. Reguitti, Y.-Z. Cai, M. Fraser, L. Galbany, M. González-Bañuelos, C. P. Gutiérrez, T. Kangas, T. L. Killestein, T. Kravtsov, P. Lundqvist, S. Moran, A. Pastorello, A. Popowicz, I. Salmaso, M. Stritzinger

High Energy Astrophysical Phenomena

Beyond Disk Truncation: X-ray Reverberation Signatures of an Outflowing Corona

The corona in black hole X-ray binaries (BHXRBs) is likely dynamic during outbursts. A mildly relativistic outflowing corona has been proposed to explain hard spectra, weak disk reflection, and the higher than expected polarization degree observed in the hard state of X-ray binaries. In this work, we investigate the spectral and timing effects of coronal outflow using Monte Carlo radiative-transfer simulations. The model consists of a geometrically thin, optically thick truncated disk and an ellipsoidal corona with a prescribed bulk outflow velocity. We calculate the Comptonized continuum, the disk reflection component, and the lag-frequency spectra for different outflow velocities and disk truncation radii. We find that increasing the outflow velocity reduces the reflection fraction through relativistic beaming, because fewer Comptonized photons irradiate the disk. For $\beta \lesssim 0.5$, the high-frequency soft lag is only weakly affected by the outflow velocity, whereas increasing the disk truncation radius shifts the zero-crossing frequency($\nu_0$) to substantially lower frequencies. Thus, timing properties can help distinguish the coronal-outflow scenario from the disk-truncation scenario. We further apply the model to MAXI J1820+070 and find that its unusual $R-\Gamma$ anti-correlation during the plateau phase is qualitatively consistent with a contracting corona accompanied by increasing bulk velocity. Disk recession alone would predict the opposite evolution of $\nu_0$ under the assumptions of our model. Future polarization calculations will provide an additional test of extended outflowing-corona models.
Read the full abstract

All authors Ke Qin, Bei You, Adam Ingram, Barbara De Marco

High Energy Astrophysical Phenomena

X-ray polarization of the Compton-thin Seyfert galaxy NGC 5506 with IXPE

We report on the observation of the Compton-thin Seyfert galaxy NGC~5506 performed with the Imaging X-ray Polarimetry Explorer (IXPE). The observation, with a net observing time of 975~ks, was partly simultaneous with NuSTAR (52~ks) and XMM-Newton (37~ks). Using a simple baseline model with IXPE data only, we find an upper limit on the 2–8~keV polarization fraction at the $99\%$ confidence level: $\Pi<3.1\%$. We then performed a spectropolarimetric analysis combining data from the three X-ray observatories. While pure spectral modeling confirms the presence of both relativistic and cold distant reflection, in addition to the primary continuum power law, their contributions lead to somewhat different constraints on the coronal $\Pi$, depending on the specific assumptions. We also performed Monte Carlo simulations with the radiative transfer code MONK, exploring various coronal geometries. Comparing these results with the IXPE data, and assuming a radially extended corona, as suggested by previous observations of similar sources, the corresponding polarization upper limits require inclination angles below 40{\deg}.
Read the full abstract

All authors Daniele Tagliacozzo, Alessandro Leonardo Lai, Andrea Gnarini, Giorgio Matt, Stefano Bianchi, Francesco Ursini, Andrea Marinucci, Alessia Tortosa, Sudip Chakraborty, Federico Ferretti, Vittoria Elvezia Gianolli, Elias Kammoun, Frédéric Marin, Herman Marshall, Riccardo Middei, Romana Mikusincova, Dawoon Kim

High Energy Astrophysical Phenomena

COBALT-BLUE: a fast and accurate model for gravitational self-lensing binary systems

We present COBALT-BLUE, a fast public model for the optical light curve of a luminous star in a non-accreting binary with a compact object. The model includes gravitational self-lensing, Doppler beaming, ellipsoidal variations, eccentric orbital geometry, and wavelength-dependent limb darkening. In sufficiently edge-on systems, it predicts sharp periodic self-lensing flares together with phase-locked orbital modulation outside the flare. We improve upon earlier models by including non-linear limb darkening laws based on stellar atmosphere simulations to calculate the self-lensing flare profile. This enables light curves to be self-consistently calculated in multiple wavebands from one set of physical parameters, thus offering the prospect of constraining binary parameters such as compact object mass from photometry alone if multiple self-lensing flares are observed. We also include the option to utilise simpler limb darkening prescriptions to save computational expense. As a proof of principle, we fit our model to ZTF data of 14 candidate self-lensing sources. The model describes the flare morphology well for most candidates, but we find that the best fitting parameters are not physically plausible. We conclude that the observed flares are most likely stellar flares and not self-lensing events. In future, we plan to use the model to systematically search optical survey data for signatures of binary motion.
Read the full abstract

All authors Danny J. Dixon, Adam Ingram, Matthew J. Middleton, Adam McMaster, Grzegorz Wiktorowicz, Hugh Dickinson, Allison Crossland, Eric C. Bellm

High Energy Astrophysical Phenomena

Wave-particle acceleration during supernova explosions

We investigate particle acceleration during the early stages of a supernova explosion via a resonant interaction between electromagnetic waves and quantum waves associated with particles. For typical progenitor masses of $(8; 12; 20)\times M_s$, we find that this mechanism can accelerate protons to energies of up to $3\times 10^{14}$ eV. Interestingly, despite the continuously declining luminosity, the particle energy initially increases, reaches a maximum, and subsequently decreases, reflecting the characteristic behavior of the resonant acceleration mechanism.
Read the full abstract

High Energy Astrophysical Phenomena

A method for enhanced gamma-proton discrimination with imaging atmospheric Cherenkov telescopes

Purpose: Very-high-energy gamma-ray astronomy is an important window to study the extreme astrophysical processes in the universe, and how to effectively distinguish between gamma photons and background signals (mainly protons) is a key technical challenge to realize very-high-energy gamma-ray detection. Spaceborne calorimeters record secondary particle spatial distributions, achieving a background rejection power of 10000-100000 near 1 TeV. Ground-based imaging Cherenkov telescopes measure Cherenkov photon angular distributions, yet they only reach a rejection factor of about 10 in this energy range, leaving substantial room for optimization. Methods: We use CORSIKA and sim_telarray to simulate H.E.S.S. detector responses. Photon emission positions are reconstructed from Cherenkov shower images. Inspired by analysis techniques for space-borne calorimeter, we propose new observables derived from the transverse distribution of Cherenkov photon emission positions to separate gamma-ray signals from proton backgrounds. We employ a likelihood-ratio method to combine single-telescope variables across multiple telescopes and compare its performance with conventional Hillas-based analysis. Results: The results show that observables built from the transverse distribution of Cherenkov photon emission positions are less sensitive to Cherenkov image size and outperform the Hillas width parameter for single-telescope gamma-proton discrimination. Furthermore, the likelihood ratio combining width yields substantially better performance than the conventional mean reduced scaled width: the improvement reaches roughly one order of magnitude near 1 TeV, and this enhancement remains above 50% for energies above 10 TeV.
Read the full abstract

High Energy Astrophysical Phenomena · Cross-listed from astro-ph.IM

Combining astrometry with pulsar timing: the first joint analysis of very low frequency gravitational waves

The pHz to sub-nHz GW regime remains largely unexplored but is crucial for mapping the early inspiral stage of SMBHBs and probing early-Universe physics. Astrometry and pulsar timing offer orthogonal and deeply complementary secular observables to investigate this frequency band. We aim to present the first joint data analysis combining real astrometric proper motions with binary pulsar timing, to search for and constrain continuous gravitational waves (CWs) sourced by SMBHBs in the ultra-low-frequency regime ($10^{-12} \le f_{GW} \le10^{-9} Hz$). We employ a Bayesian model selection and upper limit estimation framework to combine apparent proper motion displacements of $\sim 1.5 \times 10^6$ quasars from the Gaia CRF3 catalog with the line-of-sight orbital period derivatives ($\dot{P}_b$) of 11 high-precision binary pulsars. To prevent spurious detections, we heavily model instrumental and astrophysical systematics: we propagate the Galactic potential uncertainty for pulsars via Monte Carlo simulations and perform a Vector Spherical Harmonics (VSH) decomposition up to the octupole order (l=3) for quasars. We find no statistically significant evidence for a CWs signal in the joint analysis ($\ln B_{joint} = -0.42 \pm 0.03$). In the absence of detection, we set the tightest constraints to date on CW strain in the pHz band, yielding a 95% upper limit of $h_0\le6.4x10^{-11}$ at a reference frequency of $f_{ref} = 4 \times 10^{-10}$ Hz. The combined dataset achieves full sky coverage and improves single-dataset upper limits by 20%-30%. Combining orthogonal observables successfully breaks spatial degeneracies intrinsic to isolated searches. Furthermore, forecasts from inj.-rec. indicate that with the extended temporal baseline and reduced uncertainties of the upcoming Gaia DR4, this joint framework is poised to break the $h_0 < 10^{-11}$ upper limit barrier for sub-nHz CWs.
Read the full abstract

All authors Lorenzo Filipello, William Beordo, Volodymyr Akhmetov, Mariateresa Crosta, Mario Gilberto Lattanzi

High Energy Astrophysical Phenomena

Gamma-ray Quasi-periodic Oscillations in the Changing Look Blazar B2 1420+326 (OQ 334)

In the Fermi-LAT gamma-ray light curve from January 1 to August 25, 2026, of the changing look blazar B2 1420+326, we report here the detection of transient quasi-periodic oscillations (QPOs) of period ranging around 10-13 d. The source entered the flare state around May 7, 2026, and lasted till around mid-August, 2026 during which it depicted strong QPO signatures. We utilized standard techniques, including weighted wavelet Z-transform analysis, generalized Lomb-Scargle periodogram, and phase dispersion minimization, to search for and confirm these findings. Prior to the flare, the source displayed a mild damped oscillatory feature of similar periodicity with a damping timescale of around 20 d. Based on the period timescale and nature of the QPOs detected, we discuss possible physical mechanisms for the observed salient features.
Read the full abstract

High Energy Astrophysical Phenomena

A Gamma-ray Emitting Head-Tail Radio Galaxy Swimming Through the Hot Medium of the Merging Cluster Abell 3627

Head-tail radio sources are jetted active galactic nuclei (AGN) moving through the dense intracluster medium of galaxy clusters. They are extremely rare in the gamma-ray sky, likely because they are usually observed at large viewing angles. Therefore, every new gamma-ray detection of head-tail radio galaxies allows us to probe the origin of the high-energy emission in this unique class of AGN. Here we report, for the first time, the association of ESO~0137-G007 (z=0.016), a head-tail radio galaxy located at the outskirts of the merging cluster Abell 3627, with the gamma-ray source FL16Y~J1615.5$-$6034. The motion of the galaxy through the dense cluster environment and subsequent ram pressure and/or intracluster medium turbulence have led to the formation of a $>$500 kpc-long jet exhibiting a twisted, wiggly radio morphology. Its optical spectrum, obtained with the 4.1 m SOAR telescope, is devoid of emission lines, indicating that the underlying accretion activity is radiatively inefficient. From the measured line-of-sight stellar velocity dispersion, we estimate the mass of the central supermassive black hole to be $(3.36\pm1.32)\times10^9$ \Msun. We briefly discuss several radiative models to explain the observed broadband spectral energy distribution of ESO~0137$-$G007. We conclude that deeper multi-wavelength observations of this unique gamma-ray emitting head-tail radio galaxy will enable us to better understand the interaction of the jet with the hot cluster environment, thus setting the stage for the research of this class of AGN with the upcoming Square Kilometer Array.
Read the full abstract

All authors Vaidehi S. Paliya, M. Böttcher, Swayamtrupta Panda, Kiran Wani, D. J. Saikia, Mahadev Pandge, C. S. Stalin

High Energy Astrophysical Phenomena

Why most neutron star low-mass X-ray binaries accrete transiently: an evolutionary study of transient and persistent phases

A neutron star (NS) low-mass X-ray binary (LMXB), in which an NS accretes matter from a low-mass donor star, is an ideal source for probing some fundamental aspects of physics and astronomy, such as strong gravity, superdense matter, and the accretion-ejection processes. However, to reliably achieve these goals, one must adequately understand NS LMXBs, including why some accrete persistently and others transiently. Focused models, such as those based on a thermal-viscous instability in the accretion disk, are considered to explain transient accretion. However, broader perspectives, including which LMXB parameter values and phases cause transients and why there are more transients than persistents, remain poorly understood. Here, our computation of the long-term evolution of NS LMXBs addresses these questions, providing insight into LMXB parameters and phases, naturally producing more transients than persistents, and being partially consistent with the known properties of observed sources. For example, we typically find a greater fraction of persistent phase at lower orbital periods from the LMXB evolution computation, which is somewhat consistent with observations. However, a lack of full consistency calls for improving the aforementioned focused models, and our computations provide a new way to discriminate among these models.
Read the full abstract

High Energy Astrophysical Phenomena

Neutrino Emission from Proton Photon Jet Interactions in Microquasars

Microquasars are candidate sites for high energy particle production within our galaxy. Because their distances and compact emission regions limit direct observational constraints, numerical simulations are useful for connecting jet dynamics with possible multi messenger signatures. This work models neutrino production from relativistic magnetohydrodynamic microquasar jets, focusing on the proton photon channel associated with interactions between accelerated protons and ambient photon fields. Proton proton interactions are acknowledged as a possible process in dense environments, but they are not modeled in the present calculation. A ray tracing procedure is applied to the hydrodynamic simulation output to produce synthetic neutrino images from the perspective of a stationary observer. Synthetic spectra and intensity maps are presented and compared with representative sensitivities of current and future neutrino detectors.
Read the full abstract

High Energy Astrophysical Phenomena

CRAFT HTR2: Polarimetry of 64 non-repeating fast radio bursts from the updated CRAFT catalogue

We present high-time resolution spectro-polarimetric data for 34 new fast radio bursts (FRBs) discovered by the Commensal Real-time Fast Transients (CRAFT) survey on the Australian Square Kilometer Array Pathfinder (ASKAP) during the period May 2024 to June 2026. Most of these were detected by the higher-sensitivity CRAFT COherent (CRACO) detection system that was commissioned on the telescope during this period. This new sample doubles the size of the CRAFT HTR catalogue and probes a fainter population of FRBs thanks to the improved sensitivity of CRACO. We compare the distribution of extragalactic rotation measure (RM) and polarisation fraction to the CHIME and DSA catalogues. While no significant differences were seen between CRAFT and DSA, the extragalactic RM distribution seen in CHIME FRBs (which are detected at lower frequency) was substantially lower. Surprisingly, we find no significant differences in the linear polarisation fraction distribution between the three FRB catalogues, suggesting an indifference to the different telescope observing frequencies. We show tentative evidence for wider and fainter bursts possessing lower polarisation fractions; this is predominantly driven by the growing sample of unpolarised bursts that are, in almost all cases, wider ($\gg$10 ms) and fainter ($\ll$10$^{34}$ ergs s$^{-1}$ Hz$^{-1}$) than the median ASKAP detection.
Read the full abstract

All authors T. Dial, A. T. Deller, Marcin Glowacki, A. Bera, R. M. Shannon, Clancy W. James, Joscha N. Jahns-Schindler, Ziteng Wang, Akhil Jaini

Thursday 17 September

High Energy Astrophysical Phenomena · Cross-listed from astro-ph.SR

COPAS: Compact Objects and the Physics of Accretion Survey I. Accreting Binaries with P < 83 Minutes

We present the first results from the Compact Objects and the Physics of Accretion Survey (COPAS), a high-cadence photometric survey focused on the detection of accreting white dwarf binaries mainly through their outbursts using TESS. From an initial sample of 2,054 cataclysmic variable (CV) candidates observed during TESS Cycle 6 and 7, we identify 11 systems with periods below 83 minutes. This includes 7 AM CVns and 4 CVs. Our Gemini spectroscopic data confirm ZTF18aaxuusk and ASASSN-21eo as AM CVns, while Gaia21akb is confirmed as a CV. We provide the first period measurements for 6 out of the 11 systems. Complementing our TESS data with ground-based observations, we measured superoutburst recurrence times of approximately one year for the AM CVn systems ASASSN-21in and ZTF18aaxuusk, and 1.8-2.0 years for the short-period CVs ASASSN-18rd and ASASSN-15ev. Using TESS, we also estimate the mass-ratio for ASASSN-19ct (q = 0.05 +/- 0.02), which is consistent with a highly evolved low-mass donor star. Though the number of ultracompact systems detected among the total sample is low (0.54%), and likely the result of observational biases and intrinsic population rareness, the reclassification of 5 systems as AM CVns highlights the importance of continuous monitoring to identify more systems through their outbursts and short periods.
Read the full abstract

All authors Wendy Mendoza, Liliana Rivera Sandoval, Jan Kára, Ryan J. Oelkers, Yuri Cavecchi, Thomas J. Maccarone, Manuel Pichardo Marcano, Meryem K. Dağ, Simone Scaringi, Kieran O'Brien, Martina Veresvarska, Christian Knigge

High Energy Astrophysical Phenomena

Effects of Physical Cooling on the Structure of Circumbinary Disks

Radiation pressure and radiative cooling often play an important role in deciding the structure and dynamics of astrophysical objects. In this paper, we present a radiative MHD simulation of a circumbinary disk (CBD) around an equal-mass SMBBH. It employs an approximate thermodynamic model based on local thermodynamic equilibrium and diffusion cooling which captures how the internal energy evolves in an optically thick region. The disk reaches a quasi-steady state within $r=5a$ whose mass distribution differs substantially from the one obtained from simulations without physical thermodynamics: the local maximum in the azimuthally-averaged surface density and the lump are completely erased; the time-averaged inner edge becomes more eccentric; and the vertical density distribution depends strongly on the radius. These changes are attributable to the enhancement of the magnetic fields caused by the cooling-driven vertical compression of the disk. Our result emphasizes that the basic properties of the inner CBD are sensitive to both radiation and magnetic physics.
Read the full abstract

High Energy Astrophysical Phenomena · Cross-listed from gr-qc

Non-Conservative Mechanics in the Compact-Body Mass-Shell Model and its Waveform Structuring

A novel reinterpretation on the effective one-body (EOB) mapping of compact-body binary (CBB) coalescence was previously considered. Instead of the conventional “marble in a funnel" description of EOB, CBB coalescence is modeled as the concurrent rotation and contraction of a hollow shell with a constant reduced mass measure. Merger is achieved, therefore, when the shell radius reaches the total mass radius. While past works analyzed the energetics sector of the model and the radiated gravitational waves (GWs), this work considers the geometric waveform modeling aspect. More specifically, we mature the mass-shell waveform analysis initially discussed with toy-model conditions by imposing necessary non-conservative mechanics, leading to the analytical expression of usable mass-shell modeled GW-forms. This is to set the foundation for eventual computational endeavors with the analytical waveforms presented here, enabling a comparison against current and state-of-the-art models.
Read the full abstract

High Energy Astrophysical Phenomena

Energy Partition in Relativistic Electron-Positron-Ion Reconnection

Relativistic magnetic reconnection, a prime mechanism for particle acceleration, has been extensively studied in electron–ion and pair plasmas, but how the dissipated energy is shared among species in electron–positron–ion plasmas, and in particular how ions are energized within a pair-dominated layer, remains poorly understood. We develop an analytic theory that couples single-particle orbit dynamics to the collective plasma response, and we validate it with particle-in-cell simulations spanning pair-dominated to electron–ion compositions. Particles are directly energized in electric-dominated regions until the reconnected magnetic field deflects them into the outflow; the characteristic escape length is in turn set by the relativistic skin depth of the reconnection-heated mixture, including the ion contribution. Writing $\sigma_e^*$ for the magnetic energy available per particle in units of $m_ec^2$ and $\mu=m_i/m_e$, the mean lepton energy scales as $\sigma_e^*$, whereas the mean ion energy scales as $\mu^{1/3}(\sigma_e^*)^{2/3}$ when ions escape non-relativistically from the electric-dominated regions and approaches the lepton energy when their escape is ultrarelativistic. The predicted energy fractions agree with simulations across composition, mass ratio, and magnetization. The ion share of the dissipated energy can fall far below the commonly assumed ion–lepton equipartition, reducing the energy available for hadronic and neutrino emission in reconnection-powered astrophysical sources.
Read the full abstract

High Energy Astrophysical Phenomena

The Auger Radio Infill SKALA Extension (ARISE): First Air-Shower Measurements

The Auger Radio Infill SKALA Extension (ARISE) was installed at the Pierre Auger Observatory in 2025 to detect the radio emission in the band of $50-350\,$MHz of near-vertical cosmic-ray air showers of energies around and above $100\,$PeV. The array comprises 18 SKALA~v2 antennas deployed around one of the surface-detector stations of the $433\,$m array (SD-433), which provides the trigger for ARISE. We present first measurements with ARISE and show the detection of radio pulses associated with extensive air showers. Radio signals are found in coincidence with SD-433 events, and a comparison of the arrival directions reconstructed from the radio and surface detector signals confirms the measurement of air showers with ARISE. These observations provide a first assessment of the detector setup and demonstrate the potential of ARISE for radio detection of near-vertical air showers.
Read the full abstract

High Energy Astrophysical Phenomena

A mechanism for fast radio bursts and their narrow spectra

We find an emission mechanism for fast radio bursts (FRBs) capable of producing broad and narrow spectra. It operates in magnetar explosions, which launch an electromagnetic pulse from the inner magnetosphere at radii $r_0\sim 10^{7}$ cm. Plasma inside the pulse accelerates outward with Lorentz factor $\Gamma(r)\approx r/r_0$. Before exiting the magnetosphere at radius $r\sim 10^{10}$ cm, the pulse passes through ambient kHz waves that populate the magnetospheres of active magnetars. We show that the waves experience stimulated scattering inside the pulse, with the boost factor $\Gamma^2$ converting them into radial GHz radiation. This process exponentially amplifies tiny GHz seeds to extreme luminosities consistent with observed FRBs. The produced GHz burst rides inside the explosion pulse, which helps it escape to distant observers. The FRB duration and energy are controlled by the pulse duration and energy. In a saturated emission regime, the predicted burst contains an ultrastrong microsecond component.
Read the full abstract

High Energy Astrophysical Phenomena

Self-Similar Mass Spectra of Hierarchical Black Hole Mergers

Black holes that merge repeatedly carry a record of that history in their mass distribution. We treat hierarchical merging as a coagulation problem, in which a single kernel encodes how the merger rate depends on the masses involved and on cosmic time, and we derive the late-time spectra that such populations approach. The framework yields closed scaling laws, an exactly solvable benchmark, and a classification of merger environments by a single exponent measuring how strongly mergers favor or suppress massive participants. Applying it to primordial black holes requires care in translating published merger rates into kernel form; done correctly, three of the four standard binary-formation channels map onto superlinear kernels and are therefore candidates for runaway growth, in which the heaviest objects would dominate and no steady mass-conserving spectrum exists; establishing physical gelation, however, requires the full population-dependent kernel and its finite coagulation history, neither of which we settle here. Only the early three-body channel maps into the nongelling regime and admits slow self-similar growth. Numerical solutions across thirteen suppressive, nongelling kernels show that the shape of the spectrum is not fixed by the scaling exponent alone, and that accounting for the energy radiated at each merger measurably changes the high-mass cutoff and drains the population's total black-hole mass.
Read the full abstract

All authors Liam Blum, Stefano Profumo, Ryan Schantz, Quoc Ha Tran

High Energy Astrophysical Phenomena

Systematic Uncertainties and Their Impact on Gamma-Ray Searches for Dark Matter in Dwarf Galaxies

Gamma-ray observations of Milky Way dwarf galaxies provide some of the strongest constraints on dark matter annihilation and are approaching the sensitivity required to test dark matter interpretations of the Galactic Center Gamma-Ray Excess. The results, however, depend critically on the inferred dark matter distributions in these systems, conventionally described by their $J$-factors and spatial profiles. We review several sources of systematic uncertainty in dwarf galaxy $J$-factor determinations, including stellar membership, unresolved binaries, departures from equilibrium and spherical symmetry, and assumptions regarding the stellar and dark matter density profiles. Most $J$-factor determinations reported in the literature neglect several of these uncertainties, and the resulting measurements are often more dispersed than their quoted uncertainties would suggest. This indicates uncertainties of order $\Delta\log_{10}J\sim0.5$ or larger may be warranted for many of these systems. Using an ensemble of simulated Fermi-LAT observations of a stacked sample of 27 dwarf galaxies, we quantify how these unmodeled uncertainties affect the resulting constraints on the dark matter annihilation cross section. We find that adding an uncertainty of $\Delta\log_{10}J\sim 0.5$ or 1.0 in quadrature weakens the median upper limits of a stacked analysis by factors of 1.4 or 3.2, respectively. Furthermore, underestimating these uncertainties can cause such analyses to exclude the true annihilation cross section more frequently than the nominal confidence level would imply. Such undercoverage is especially pronounced in stacked analyses and persists even when the adopted and true $J$-factor uncertainties agree. Mismodeling spatially extended halos as point sources, as is sometimes done, further strengthens the resulting limits and increases the probability that true signals will be incorrectly excluded.
Read the full abstract

High Energy Astrophysical Phenomena

Little Red Dots As Super-Eddington Fountain Flows

Little red dots (LRDs) may be powered by supermassive black holes (SMBHs) accreting above the Eddington limit. Spectroscopy of LRDs often shows absorption troughs blueshifted by $\sim100-300\,{\rm km\,s^{-1}}$, implying a slow wind. This is puzzling: super-Eddington disks drive much faster winds, $\gtrsim10^3-10^4\,{\rm km\,s^{-1}}$. We argue that LRDs are super-Eddington SMBHs viewed off-axis and engulfed in a slow wind that covers most sight-lines; the fast wind escapes near the poles. Trapped light puffs the inner disk into a quasi-spherical envelope that launches the slow wind. The wind may be “photon-tired”, meaning the light only barely unbinds it. Such marginally unbound winds lead to fountain flows, where some gas escapes and the rest falls back. We model the envelope with idealized, spherically symmetric “marginally unbound” ($v\sim v_{\rm esc}$) and “photon-tired” winds. We feed these profiles into the radiative transfer code Sirocco to study how the wind reprocesses the light from an accreting $10^6\,M_\odot$ SMBH. We describe most of the wind as a “Balmer cocoon” – a Compton-thick region in which depletion of Balmer continuum photons ($h\nu >3.4\,{\rm eV}$) rather than Lyman continuum photons ($h\nu > 13.6\,{\rm eV}$) keeps the gas ionized. The spectra span the range of LRD-like sources: “little blue dots” at lower outflow rates ($\sim2.5\,M_\odot\,{\rm yr}^{-1}$); V-shaped LRDs with a Balmer break ($\sim5-10\,M_\odot\,{\rm yr}^{-1}$); and red LRDs with full breaks ($\sim15\,M_\odot\,{\rm yr}^{-1}$). Our Balmer line profiles show P~Cygni features atop broad, exponential wings, as is observed. The break is possible at lower wind densities than in LTE models because Lyman~$\alpha$ trapping sustains our $n=2$ hydrogen population and electron scattering enhances the optical depth. Our arguments are also applicable to winds from supermassive stars or quasi-stars.
Read the full abstract

All authors Nicholas Kaaz, Eliot Quataert, Taya Govreen-Segal, Hanpu Liu

High Energy Astrophysical Phenomena

Systematic Effects of Hydrogen and Helium Atmosphere Mismatch on Radius Inference in PSR~J0740+6620-like Synthetic NICER Data

Constraints on neutron star radii provide insight into the properties of the cold, dense matter in their interiors. Previous studies using synthetic Neutron star Interior Composition Explorer (NICER) pulse waveform data have demonstrated that radius inferences derived therefrom are robust against several classes of modeling systematics. Here we explore the consequences of assuming the wrong atmospheric composition, using synthetic data based on the $\sim 2.1~M_\odot$ pulsar PSR~J0740$+$6620. We find that the assumption of a hydrogen atmosphere when the synthetic data assumed a helium atmosphere, or vice versa, produces little bias in the inferred radius at the spot-to-background ratio of the actual PSR~J0740$+$6620 data. However, when we increase the spot-to-background ratio by a factor of $\sim20$ while keeping the total number of counts fixed at the observed $\sim 5.5\times 10^5$, we find that composition mismatch can produce significantly biased radius estimates while remaining hidden inside a fit with statistically acceptable residuals. Even in these cases, the Bayesian evidence consistently identifies the correct atmospheric model. Our findings reinforce the importance of using the Bayesian evidence for model comparison and goodness-of-fit tests.
Read the full abstract

All authors Isiah M. Holt, M. Coleman Miller, Alexander J. Dittmann, Frederick K. Lamb

High Energy Astrophysical Phenomena

A systematic search for long GRBs without supernovae: global properties and rate of events

Some long duration gamma-ray bursts (LGRBs) are not associated with bright supernovae (SNe), indicating that their explosion mechanism differs from the traditional collapsar channel. SN-less LGRBs are rare, thus it remains unclear whether they are sporadic outliers or a substantial component of the GRB population. In this work, we use two decades of Swift observations to conduct a systematic search for nearby ($z\lesssim$0.35) LGRBs without a confirmed SN association. Our analysis reveals a heterogeneous population linked, in part, to star-forming or dusty environments and, in part, to quiescent galaxies. Through the analysis of their high-energy properties, afterglows, and environments, we identify a sample of LGRBs with robust SN limits and distance scales, more than doubling the number of candidate SN-less events. From our fiducial sample, we derive a volumetric rate, not corrected for beaming, of $R = (0.5 \pm 0.2) {\rm~Gpc^{-3}\ yr^{-1}}$, comparable to the apparent rate of short duration GRBs (SGRBs) and standard LGRBs with SNe with luminosities $L \gtrsim\,10^{50}\,{\rm erg\,s^{-1}}$. For plausible beaming factors, the inferred rate of SN-less LGRBs is consistent with several compact binary merger channels, although the allowed parameter space is already narrow if the entire population is attributed to a single progenitor channel. Owing to the small number of candidate events and the uncertainties in their distance scale, our sample leaves the low-luminosity tail unconstrained. A larger sample will be needed to determine whether they are the faint tail of known progenitor channels or evidence of a distinct route to relativistic explosions.
Read the full abstract

All authors Yu-Han Yang, Eleonora Troja, Nicolas Monsalves, Massine El Kabir, Narjes Shahamat Dehsorkh, Hira Waseem

High Energy Astrophysical Phenomena

Low Frequency Radio Imaging Study of PSR J1818-1607

We report on the low-frequency radio observations of magnetar PSR J1818$-$1607 carried out with the upgraded Giant Metrewave Radio Telescope at band 3 (300$-$500 MHz) and band 4 (550-750 MHz). We have identified the continuum source associated with the magnetar and report variations in flux density and its spectral index. The flux density timeseries for the magnetar reveals the presence of two potential radio flaring episodes with varying spectral properties. We discuss the implications of the spectral index changes on the potential emission mechanisms for radio-loud magnetars. We also report non-detections of the continuum source as well as no direct indication for an associated diffuse emission from multiple archival radio imaging surveys.
Read the full abstract

All authors A. Yadav, M. P. Surnis, B. C. Joshi, M. Bagchi