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)

Tuesday 8 September

High Energy Astrophysical Phenomena

Dispersion Measure Variability in Fast Radio Bursts from Photoionization

Magnetars became favored engines of fast radio bursts (FRBs) following the discovery of a luminous radio burst coincident with a hard X-ray flare from a Galactic magnetar. Several repeating FRB sources exhibit time-variable rotation measures and compact, spatially coincident persistent radio emission, consistent with energetic-particle nebulae confined by young supernova ejecta. Secular changes in the dispersion measure (DM) of repeating FRBs have also been observed, offering a complementary probe of their local environments; for example, the DM of FRB 121102 rose until 2019 before declining in recent years. Although rising DM evolution has been attributed to shock ionization, shocked ejecta can cool efficiently through metal-line emission and recombine, especially if mixed with cooler gas. Here we argue that DM variations of the observed magnitude and timescale instead arise from changes in the ionization state of supernova ejecta irradiated by X-rays from a time-variable central engine. The dominant rapidly variable contribution comes from the dense, weakly ionized shell swept up by the expanding nebula, whose ionization and recombination times are shorter than those of the more extended ejecta. The model predicts that enhanced FRB activity should be accompanied by rising DM, with a response smoothed over the ionization/recombination time and superposed on a slower secular decline from ejecta expansion. In FRB 121102, the DM maximum occurred close to the burst-rich 2018-2019 activity episodes. If the recently reported renewed activity is sustained, its DM decline should flatten and may reverse into a fresh rise over the coming years.
Read the full abstract

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

Ruling out matter-induced bumps in black-hole effective potentials

Matter fields localized around black holes can produce small secondary bumps in the effective potentials governing black-hole perturbations, potentially leading to quasinormal-mode spectral instabilities and late-time echoes. Considering physically motivated models of anisotropic fluids surrounding a spherically symmetric black hole, we show that the appearance of such a secondary bump is accompanied by either singular matter profiles or divergent characteristic speeds in the fluid. These features indicate that the corresponding matter configurations are unphysical. We further demonstrate that a secondary bump in the effective potential does not necessarily imply the existence of additional light rings in the spacetime.
Read the full abstract

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

Deus Ex Statistica: A Statistical Solution to the binary-binary Outcome of the Chaotic, Non-Hierarchical Four-Body Problem

We present an analytical, statistical solution to the binary-binary (2+2) outcome of the chaotic non-hierarchical four-body problem. The solution is based on the density-of-states formulation pioneered by J. J. Monaghan. The method skips the computationally expensive integration of the equations of motion, and instead samples the outcome from the chaotic phase-space, subject to conservation of energy, momentum, and angular momentum. From the joint distribution, we extract marginal distributions of several key parameters using Monte-Carlo integration, and numerically verify them by comparing to an identical ensemble of scattering experiments produced by the FEWBODY code. From the comparison, we identify a regime not represented by the density-of-states formulation: the hard binary regime, where one binary is much harder then the four-body energy scale, and the system acts as an effective three-body system. We hypothesize that at this regime the system's probability distribution spreads over an effective reduced three-body chaotic phase space. The process of transfer from four-body to three-body phase-space is still not understood, and represents the next natural extension of density-of-states methods.
Read the full abstract

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

Overmassive No More: The Case for Little Red Dots Hosting Black Hole Seeds as Massive as Single Supermassive Stars

Little Red Dots (LRDs) display singular properties unlike any known class of AGN or galaxies, motivating novel mass estimators for their central engines. Inspired by their similarities to stellar phenomena, here we interpret the LRD continuum as being produced by a pseudo-photosphere. We fit tailored stellar atmosphere models to host-subtracted LRD central engines ("black hole stars," BH*s) represented by stacks of $117$ objects. Typical BH* continuum spectra are well fit by models in a narrow range of temperatures ($T_{\rm eff}\approx4200-4800$ K), with bolometric luminosities $\approx10^{43-45}$ erg s$^{-1}$, implying pseudo-photospheric radii $\approx700-2000$ au. Based on these parameters, we explore four different approaches to deriving BH* masses: 1) using the surface gravity from atmosphere models; 2) appealing to the resemblance to super-Eddington phenomena; 3) approximating the escape velocity from the outflowing material; and 4) exploiting the lack of variability to bound the dynamical time. For the typical BH*, all of these methods yield remarkably consistent masses of $\approx10^{4-5}\,M_\odot$, implying a highly super-Eddington luminosity of $L_{\rm{bol}}/L_{\rm{Edd}}\sim5-50$. These mass estimates place BH*s within the scatter of the local scaling relation between black hole mass and host galaxy stellar mass, providing a self-consistent alternative to "overmassive" black holes that lie $2-3$ dex above it. Crucially, our derived masses are consistent with BH*s arising from single supermassive stars (SMSs), whose masses cannot exceed $\approx10^{5-6}\,M_\odot$ due to general relativistic instabilities. Furthermore, for our derived $L_{\rm bol}/L_{\rm Edd}$, the sharp cutoff of the LRD luminosity function matches the maximum theoretical mass of an SMS. With LRDs, we may therefore be directly observing the birth of heavy black hole seeds.
Read the full abstract

All authors Wendy Q. Sun, Rohan P. Naidu, Hanpu Liu, Anna de Graaff, Jenny E. Greene, Jorryt Matthee, Chris Ashall, John Chisholm, Anna-Christina Eilers, Qinyue Fei, Kasper E. Heintz, Daichi Hiramatsu, Vasily Kokorev, Joel Leja, Zhaoran Liu, Priyamvada Natarajan, Pascal A. Oesch, Robert A. Simcoe, Alberto Torralba, Andrea Weibel

High Energy Astrophysical Phenomena

Decoding the Imprints of Energy-Momentum Squared Gravity in Neutron Stars with Machine Learning Analysis

Neutron stars (NSs) provide a unique laboratory for testing gravity in the strong-field regime and for searching for deviations from General Relativity (GR). In this work, we investigate the effects of Energy-Momentum Squared Gravity (EMSG) on NS structure and examine whether its signatures can be identified from observable stellar properties using supervised machine learning (ML). We solve the modified Tolman-Oppenheimer-Volkoff equations for approximately $10^{4}$ nuclear equations of state (EOSs) for EMSG coupling parameters $\alpha=\{-5.01,-2.50,0,+2.50,+5.01\}\times10^{-38}\,\mathrm{erg}^{-1}\mathrm{cm}^{3}$, and calculate the gravitational mass $M$, radius $R$, dimensionless tidal deformability $\Lambda$, and fundamental $f$-mode oscillation frequency for each stellar configuration. Imposing observational constraints on $M$, $R$, and $\Lambda$, we split our datasets into train and test sets, and we employ Random Forest (RF), K-Nearest Neighbors (KNN), Support Vector Machine (SVM), Logistic Regression (LR), and Gaussian Naive Bayes (GNB) to classify the representative sectors $\alpha=\{-5.01,0,+5.01\}\times10^{-38}\,\mathrm{erg}^{-1}\mathrm{cm}^{3}$ using $(M, R,\Lambda,f)$. The RF classifier performs best, achieving an accuracy of approximately $99.85\%$ with precision, recall, and F1-scores exceeding $99.8\%$, while KNN also achieves accuracy above $99\%$. The nearly diagonal confusion matrices demonstrate that the observationally viable NS configurations associated with different EMSG sectors remain highly separable in the multidimensional observable space. Our results show that NS observables retain robust signatures of EMSG even after observational filtering, establishing ML-assisted NS observations as a promising complementary approach for probing modified gravity with current and future multi-messenger observations.
Read the full abstract

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

NRHJSur3dq8: a spectral numerical-relativity surrogate family for non-eccentric aligned-spin binary-black-hole waveforms

We present NRHJSur3dq8, a family of numerical-relativity (NR) surrogates for non-eccentric aligned-spin binary-black-hole waveforms inspired by an action-angle parameterization. The model provides both the dynamics: the Hamiltonian values, azimuthal action $J_\phi(\Phi)$, the orbital clock $\tau(\Phi)$ as functions of its phasing, and the radiative pieces: the waveforms, their near-analytic derivatives, and modeling errors associated with each. The surrogated variables are conditioned into pieces that vary slowly on the radiation-reaction timescale, represented on shared $hp$-adaptive Chebyshev elements, and regressed across the intrinsic parameters with a Gaussian-process model. A time-parameterized merger-ringdown set of elements is attached where the adiabatic description breaks down. Leave-one-out validation of the central NRHJSur3dq8_AA flavor against NR reaches a $(2,2)$ unweighted mismatch of $1.3^{+32}_{-1.2}\times10^{-6}$ and an all-mode mismatch of $2.3^{+59}_{-2.2}\times10^{-6}$ over the full inspiral-merger-ringdown span. Sky-averaged and SNR-weighted over $32$ orientations per simulation, the released model's full-IMR in-sample mismatch against the NR simulations used for training is $8.9\times10^{-7}$ to $2.8\times10^{-5}$ across total masses $40$-$120\,M_\odot$. Waveforms are hybridized using post-Newtonian waveforms. Beyond accuracy, NRHJSur3dq8 provides two capabilities uncommon in NR surrogates: nearly analytic derivatives of the strain with respect to the physical parameters, and a calibrated Gaussian-process predictive uncertainty of the model. A full non-hybridized inspiral-merger-ringdown waveform evaluates in $12$ ms on a single CPU call. Batched evaluation amortizes this to $0.39$ ms per waveform at a batch of $4096$, measured at a starting frequency of $20$ Hz.
Read the full abstract

High Energy Astrophysical Phenomena

Hierarchical Inference of the Supermassive Black Hole Binary Merger Rates from Joint Searches using Pulsar Timing Arrays

Gravitational-wave searches do more than identify individual sources - they provide a way to infer the underlying astrophysical populations that produce them. Currently, Pulsar Timing Arrays (PTAs) constrain the supermassive black-hole binary (SMBHB) merger-rate density through both the stochastic gravitational-wave background (SGWB) and searches for individual SMBHB signals. The latter implies converting source upper limits into rate upper limits using detection efficiencies estimated empirically. We instead develop a hierarchical Bayesian framework that places the SMBHB merger-rate density model directly inside the PTA likelihood. The unknown catalog of individual SMBHB merger signals is modeled as a Poisson point process on source-parameter space, so that the number of merger signals is inferred from the data rather than imposed by a detection threshold. We describe two computational routes: a catalog-marginal likelihood based on analyses with fixed numbers of candidate sources, and an explicit transdimensional sampling approach that jointly samples population hyperparameters, latent merger signals, and noise. We validate the method with a toy merger-only population model and then apply it to an astrophysical SMBHB merger-rate density model that jointly predicts the SGWB amplitude and the expected merger-catalog size. In simulations, an individual merger signal adds complementary information to the SGWB and can tighten constraints on the merger-rate population.
Read the full abstract

High Energy Astrophysical Phenomena

Wolf-Rayet stars and black holes: mass distributions and evolutionary connection

The form and relationship between the mass distributions of Wolf-Rayet stars and black holes are analyzed using a direct method for estimating the distribution function from a sample of mass values with individual distributions obtained from current observations. It is shown that, taking into account the mass loss of WR stars by the end of this evolutionary stage, the mass distribution of black holes is close to that of the CO-cores of WR stars before gravitational collapse, $M_\mathrm{BH}\simeq M_\mathrm{CO}$.
Read the full abstract

High Energy Astrophysical Phenomena

Detection of episodic Ultra-fast Outflows in the Low-luminosity quasar Mrk 205 using XMM-Newton data

Ultra-fast outflows (UFOs) from Active Galactic Nuclei (AGN) are winds which reach velocities up to 0.3c, and show mass outflow rates of 0.01-1 $M_{\odot} yr^{-1}$ and kinetic energies of $10^{42-45}$ erg/s. Their extreme properties allow UFOs to carry substantial kinetic power, sufficient to drive AGN feedback by injecting energy and momentum into the interstellar medium, regulating the host galaxy evolution and the supermassive black hole growth. In this work, we present a detailed analysis of XMM-Newton X-ray spectra of the low-luminosity quasar Mrk205. We find evidence of multi-component ultra-fast outflows, traced by two components detected in the high-resolution data, showing velocities in the range of $\sim$0.07c to $\sim$0.1c. In addition, we marginally detected an ultra-fast outflow feature in the Fe K band. According to our analysis, Mrk 205 shows episodes of wind launches, rather than persistent winds. We further discuss whether the energy carried by the UFOs in Mrk 205 is capable of producing significant feedback on its host galaxy.
Read the full abstract

All authors César Ivan Victoria-Ceballos, Anna Lia Longinotti, Yair Krongold, Giorgio Lanzuisi, Miriam Gudiño, Armando Lara-DI

High Energy Astrophysical Phenomena · Cross-listed from hep-ph

A Vector-Like Lepton Interpretation of the High-Energy Nuclear Recoil Candidate in LUX-ZEPLIN

The LUX-ZEPLIN (LZ) experiment has recently reported a nuclear-recoil candidate at $E_R=248\pm23_{\rm stat}\pm23_{\rm sys}\,{\rm keV}$ in an exposure of $2.84\,{\rm tonne\,yr}$. At this unusually high recoil energy, the spin-independent (SI) xenon response is strongly suppressed, motivating dark matter interactions with a harder recoil spectrum. We show that elastic spin-dependent (SD) scattering through the nonrelativistic operator $\mathcal O_4=\mathbf S_{\rm DM}\cdot\mathbf S_N$ provides a simple realization of such a spectrum. We consider singlet-doublet Majorana dark matter, for which a diagonal $Z$ coupling generates the required SD interaction. The accompanying Higgs-mediated SI interaction would generically produce too many low-energy recoils, but can be suppressed along a Higgs blind spot while retaining a nonzero $Z$ coupling. We find a region of the Higgs blind spot that simultaneously predicts an $\mathcal O(1)$ event rate in the LZ high-recoil window and reproduces the observed thermal relic abundance. The same SD interaction leads to solar capture, allowing this region to be tested independently with solar-neutrino searches. Current IceCube limits already probe part of the LZ-compatible thermal parameter space, while a viable region remains. Additional xenon exposure and improved solar-neutrino searches can therefore provide complementary tests of this interpretation.
Read the full abstract

High Energy Astrophysical Phenomena

Resolving structure within the iron line profile of Serpens X-1 with XRISM and NuSTAR

We present the first simultaneous XRISM and NuSTAR observations of the neutron star low mass X-ray binary Serpens X-1. We perform spectral modeling on the joint observations of XRISM/Resolve, XRISM/Xtend, and NuSTAR, testing Comptonization and double thermal continuum model prescriptions. We find that a hybrid double thermal model adequately describes the underlying continuum of the source while the Comptonization models predict largely unphysical parameter values. We perform reflection modeling and confirm with tight constraints an inclination of 5^{\circ} \pm 1^{\circ} and inner disk radius of 6.6 \pm 0.6 Rg for the source. We show that the spectral resolution of XRISM/Resolve allows for the determination of a unique inner disk radius, which has the potential to constrain neutron star spin in other systems. We discuss discrepancies in spectral shape between XRISM/Resolve and NuSTAR above \sim 8 keV. We also present lightcurve and spectral analyses of 8 Type-I X-ray bursts that occurred during the observations.
Read the full abstract

All authors Hayden Hall, Renee M. Ludlam, Edward Cackett, Jon M. Miller, Javier A. García

High Energy Astrophysical Phenomena

Strong Polarization Signatures from Magnetically Stabilized Luminous Thin Accretion Disks

We utilize the Monte Carlo radiation transport capabilities of Cosmos++ to explore the polarization and Faraday rotation of radiation emitted from a set of general relativistic radiation magnetohydrodynamic simulations of magnetically stabilized, thin, black hole accretion disks. The guiding question is whether or not the Faraday rotation depolarizes the radiation to such a degree as to be inconsistent with the relatively high polarization measurements coming from the Imaging X-ray Polarimetry Explorer (IXPE). After first confirming that our code reproduces expected polarization results for electron clouds, planar scattering atmospheres, and AGN disks, we demonstrate the polarization and Faraday rotation capabilities using a Novikov-Thorne disk threaded with a purely toroidal magnetic field. We then analyze temporal snapshots from three different simulations of luminous, thin accretion disks threaded with different magnetic field configurations. We find that the effects of Faraday rotation are generally modest over the energy range of interest, since the strongest fields present in these simulations ($\gtrsim 10^8$ G) are mostly hidden beneath the photosphere. We can easily produce polarization degrees $\ge 4$% for sources seen mostly edge-on. As one example of a relevant target, we compare our results to an IXPE observation of Cyg X-1, finding that we can match the polarization degree quite well, though our polarization angle is rotated approximately $90^\circ$ with respect to the observed data. We speculate on a few possible explanations for this.
Read the full abstract

All authors P. Chris Fragile, Peter Anninos, Madeline Breiling, Derrick Pickrel, Gabrielle Ruminski

High Energy Astrophysical Phenomena · Cross-listed from nucl-th

How strange: Phase diagrams with 3 critical points

We show that the Chiral Mean-Field model (CMF) can produce a phase diagram with three critical points: the nuclear liquid-gas transition, quark deconfinement, and a strangeness driven transition. The strangeness driven transition separates a mainly nucleonic phase from one dominated by hyperons and baryon resonances. We discuss the compositional change in these transitions and their possible signatures in heavy-ion collisions and neutron star mergers.
Read the full abstract

All authors Mateus Reinke Pelicer, Nikolas Cruz Camacho, Rajesh Kumar, Veronica Dexheimer, Jacquelyn Noronha-Hostler

High Energy Astrophysical Phenomena · Cross-listed from hep-ph

Axion Hair and Pulsar Electrodynamics: modelling, discharge dynamics, and particle-in-cell simulations

In a companion paper, we demonstrated that static axion field gradients sourced by dense nuclear matter (axion hair) can dominate the near-field electrodynamics of old rotation-powered pulsars, leading to new constraints on light QCD axions and on CP-violating axion-nucleon interactions. This article provides the extended theoretical and numerical framework underlying those results. We begin by providing a detailed description of the sourcing of axion hair from dense nuclear matter, computing self-consistent field profiles for each interaction across the relevant parameter space. We then study the modification to the electrodynamics induced in the polar gap region by these axion gradients; this is done at the analytic level by studying the modification induced by axion field gradients on the effective discharge parameter (computed in the force-free limit of the split monopole magnetic field configuration, and looking at leading deviations from the force-free limit for dipolar field configurations), and numerically by developing dedicated 1D particle-in-cell simulations which capture the leading order dynamical behavior near the star. Our results demonstrate that axion hair serves to either enhance acceleration, or enhance screening, where the relevant effect changes between the northern and southern hemispheres of the star, and between the field lines which support out-flowing and return currents.
Read the full abstract

All authors Samuel J. Witte, Andrea Caputo, Stefan Stelzl, Alexander Chernoglazov, Alexander A. Philippov, Surjeet Rajendran

High Energy Astrophysical Phenomena

Phase-Resolved Ultra-High-Energy Emission Defines an Energetic Boundary for Compact-Object Engines and Dense-Matter Structure

We show that phase-resolved ultra-high-energy emission from LS I \(+61^\circ303\) defines a critical energetic boundary linking the compact-object engine, circumstellar interaction depth, and dense-matter structure. Combining the GeV–TeV–UHE spectral transition with an orbital acceleration–transparency sieve and branch-resolved stellar sequences for 116 certified EoS, we find that the measured \(25\)–\(100\) TeV component remains energetically viable for all 116 canonical-mass EoS and therefore does not discriminate among them. EoS discrimination arises only in the deliberately more demanding pure-hadronic \(>1\) TeV stress-test construction: at 2.0 kpc and \(1.4\,M_\odot\), this boundary gives \(I_{45,\rm crit}=1.5504\), with 82 of 116 EoS satisfying the energetic condition at\(\dot P_{\rm ref}=3\times10^{-15}\,{\rm s\,s^{-1}}\). A separately applied tidal diagnostic, \(\Lambda_{1.4}\leq580\), is satisfied by 43 EoS, with nine satisfying both conditions. These models span several microscopic classes but occupy a narrow structural region, indicating structural overlap rather than selection of a unique EoS family. A separate GeV power-budget test requires \(\dot P_{\rm rot}\gtrsim 5.65\times10^{-14}f_\Omega\,{\rm s\,s^{-1}}\) at 2.0 kpc. For \(f_\Omega\sim1\), the source lies well above the UHE crossing and the EoS discrimination disappears. The UHE observation therefore identifies the engine–environment regime in which differences in neutron-star structure can become observationally discriminating, rather than defining a universal dense-matter selector.
Read the full abstract

High Energy Astrophysical Phenomena

Enabling new sampling strategies for continuous-wave analyses by integrating Bilby into PyFstat

Continuous-wave searches for rotating neutron stars usually cover wide parameter spaces and produce a large number of candidate signals. Stochastic sampling methods therefore play an important role in candidate follow-up and parameter estimation. We present the integration of the Bayesian inference library BILBY into the open-source continuous-wave analysis package PyFstat, enabling access to the broad range of stochastic samplers available through BILBY while preserving the existing PyFstat analysis infrastructure. The implementation is validated through injection studies in Gaussian noise for two representative single-stage follow-up configurations: candidates from a directed search and candidates from an all-sky search. Using the nested sampler DYNESTY, we obtain detection efficiencies consistent with the theoretical sensitivity predictions. In this setup, DYNESTY performs comparably to the existing PTEMCEE-based implementation of PyFstat for the directed case, and also provides an effective single-stage approach for the all-sky candidate follow-up considered here, where we did not identify a PTEMCEE setup with similar recovery performance. Percentile-percentile tests for DYNESTY further show well-calibrated credible intervals, demonstrating reliable parameter estimation. While a full exploration of DYNESTY settings, multi-stage setups, and other samplers is left for future work, these results demonstrate that this flexible and publicly available framework opens useful new possibilities for Bayesian continuous-wave analyses.
Read the full abstract

High Energy Astrophysical Phenomena

Compton-induced HE – VHE $\gamma$-rays from pair cascade emissions in the SEDs of radio galaxies: NGC 1275

In the scheme of active galactic nuclei (AGNi), the blazar subclass is particularly bright in $\gamma$-rays. This is because of their emissions being beamed into our line of sight, as opposed to their misaligned parent population, radio galaxies. This work presents results of a Monte-Carlo code that propagates jet-collimated $\gamma$-rays through the magnetized AGN environment, leading to secondary cascade emissions from relativistic electron-positron pairs. The resulting multi-wavelength emissions are used to fit the broadband spectral energy distribution (SED) of the radio galaxy NGC 1275. Using known physical parameters of NGC 1275 from the literature, we find that cascade emissions from stochastic processes have the potential to reproduce the broadband SEDs of radio galaxies. We also find that the interplay between accretion disk and broad-line region (BLR) photons may provide insights into the production region(s) of $\gamma$-rays in AGNi. While dense accretion disk and BLR radiation fields are vital to the development of cascades, excessive densities lead to the suppression of cascade synchrotron emission and attenuation of \,TeV $\gamma$-rays. Furthermore, the magnetic field is important in shaping the cascade $\gamma$-ray spectrum.
Read the full abstract

High Energy Astrophysical Phenomena

A Two-Parameter Framework for the Diffuse Supernova Neutrino Background

Following important experimental developments, tests of the Diffuse Supernova Neutrino Background (DSNB) are entering a new phase that might lead to a discovery. This motivates the development of simple and broadly applicable tools for comparing theoretical predictions and interpreting experimental results. We show that, over the energy range relevant to current detectors, the DSNB is well described by an exponential spectrum, characterized by the differential flux at 20 MeV, $\phi_{20}$, and a spectral energy scale, $E_0$. We fit the spectra of 24 representative models and demonstrate that the exponential description provides an excellent approximation. We then recast the recent Super-Kamiokande results in the $(\phi_{20},E_0)$ plane, showing that model-specific best fits and upper bounds can be interpreted naturally within this framework. Our results indicate that, above realistic detection thresholds, and until the precision phase is reached, the observable DSNB is effectively a two-parameter phenomenon. We advocate that the $(\phi_{20},E_0)$ parameterization be adopted as a standard, to be used routinely to report future theoretical predictions and experimental results.
Read the full abstract

High Energy Astrophysical Phenomena

A Unified Timescale Relation for Quasi-Periodic Eruptions and Repeated Nuclear Transients

Quasi-periodic eruptions (QPEs) and recurrent nuclear transients (RNTs) exhibit recurrent high-amplitude flares from galactic nuclei, yet their characteristic timescales remain poorly understood. In this work, we compile a sample of these systems and investigate empirical scaling relations between flare timescales and black hole masses. We find that the recurrence timescale exhibits a positive but highly scattered dependence on black hole mass in the combined QPE and RNT sample, approximately following $t_{\rm rec}\propto M_{\rm BH}^{1.19^{+0.52}_{-0.47}}$ with an intrinsic scatter of 0.95 dex. Remarkably, we uncover a tight nearly linear relation between recurrence time and flare timescale for QPEs and RNTs, described by $t_{\rm rec}\propto t_{\rm rise}^{1.00\pm0.07}$ with an intrinsic scatter of 0.29 dex. This relation extends across timescales from hours for QPEs to months and years for nuclear transients. We further find that QPEs and RNTs approximately follow a common empirical relation between recurrence time and flare rise time, although the physical origin of this relation remains uncertain. Our findings reveal a common phenomenological timescale link across RNTs, providing a practical framework for characterizing their temporal behavior.
Read the full abstract

All authors Shifeng Huang, Tinggui Wang, Ning Jiang, Yibo Wang, Zhenfeng Sheng, Tian-Yu Xia, Jiazheng Zhu, Zheyu Lin, Jie Lin, Ji-an Jiang, Kenta Taguchi, Keiichi Maeda

High Energy Astrophysical Phenomena

Stratified Wind Geometries of GRS 1915+105: Density Profiles and Radii

The low-mass X-ray binary GRS 1915+105 exhibits dramatic spectral and flux variability, often accompanied by highly ionized outflows. In this work, we analyze four archival Chandra/HETG observations spanning over a decade, which vary greatly in their ionizing spectrum from soft (low inner disk temperature) to hard (high inner disk temperature). We measure the wind properties, using the insightful tool of the absorption measure distribution (AMD, a distribution of column density with ionization). While the soft epoch reveals a steep, monotonically rising AMD, hard epochs feature a distinct turnover at high ionization parameter ($\log\xi>3.3$). This is the first time we observe a negative AMD slope in an outflow, which we interpret as a localized density inversion, which could point to radiation-pressure compression zones. We utilize the AMD shape and total column density measurements to calculate the density profile and distance of the wind from the central source. For the hard states, we obtain densities of $10^{12}-10^{13}$ cm$^{-3}$ at radii of $10^{10}-10^{11}$ cm. Such measurements are usually not directly possible in X-ray spectroscopy. Here we constrain the inner radius and density of the outflow to a factor of two. The high-$\xi$ values up to $\log\xi\sim6.5$ with no turnover in the soft state imply a smaller radius of a few $10^{8}$cm, which is in tension with the slow velocities observed. We demonstrate that this is a generic tension between compact launch radii and sub-escape velocities of high-column, high-ionization outflows observed in stellar-mass black hole binaries.
Read the full abstract

All authors Noa Keshet, Ehud Behar, Jon M. Miller, Joey Neilsen

High Energy Astrophysical Phenomena

Supermassive black hole binaries in the multi-messenger context of ground and spaceborne VLBI

Evidence of a gravitational wave background suggests the existence of a population of sub-parsec supermassive black hole binaries (SMBHBs), with characteristic angular separations on the order of 1-10 $\mu$as. Spaceborne extensions of Very Long Baseline Interferometry (VLBI) and the next generation ground arrays introduce the possibility of directly imaging these SMBHBs. In this work, a binary spectral energy distribution model is used to predict the detectability of SMBHBs with ground and spaceborne VLBI. We consider the Black Hole Explorer (BHEX), a proposed spaceborne VLBI mission, as our primary case study. We explore the detectable SMBHB parameter space and identify distinguishable binary signatures in the visibility domain. We find that for a flux-density-limited sample, ground array observations are more effective at detecting a wider region of the binary parameter space, with $M_\mathrm{tot} \gtrsim 10^9$ solar mass systems detectable out to redshift $z=0.075$ and beyond, depending on mass ratio. Conversely, inclusion of a spaceborne element such as BHEX, offering finer angular resolution ($\sim6\,\mu$as) and sampling of the (u,v) plane not limited by Earth rotation synthesis, will provide significant benefits in constraining binary properties. An illustrative Fisher analysis shows improvements in characterisation of the separation and position angle of SMBHBs by a factor of $\sim4$. Near-future ground and/or spaceborne VLBI may achieve the first direct observation of a SMBHB, contributing significantly to multi-messenger studies of such systems with pulsar timing arrays and observations across the electromagnetic spectrum.
Read the full abstract

All authors B. Hudson, L. I. Gurvits, D. J. D'Orazio, C. Tiede, T. M. Eubanks, E. Mooij

High Energy Astrophysical Phenomena

Constraints on the Low-frequency Radio Emission of the Galactic FRB Source SGR 1935+2154

We present a search for radio pulses from the Galactic magnetar SGR 1935+2154, a well-known source of fast radio bursts (FRBs), at $\sim$110 MHz using the Large Phased Array (LPA) of the Pushchino Radio Astronomy Observatory. Data from two active periods in 2020 (March – May and September – November, with $\sim 3.5$ minutes of daily coverage) were analyzed with new methods tailored to both FRB-like single pulses and pulsar-like periodic signals. No significant FRB-like pulses were found. Using Monte Carlo simulations, $3\sigma$ upper limits were derived for the burst rate: for a log-normal energy distribution the limit is $\sim$${10}^{1.5}~{\rm{d}}^{-1}$ for a mean of average monochromatic isotropic luminosity $L_{\nu{\rm ,mean}}\sim1.3\times{10}^{29}~{\rm{erg~s^{-1}~ {Hz}^{-1}}}$ and a natural log-space scatter of $\sigma\sim0.85$; while for a power-law distribution it is $\sim$${10}^{1.8}~{\rm{d}}^{-1}$ for an index $\beta\lesssim3.0$ and a minimum average monochromatic isotropic luminosity $L_{\nu{\rm{,min}}}\lesssim0.7\times{10}^{25}~{\rm{erg~s^{-1}~{Hz}^{-1}}}$. When folded at the known 3.24781628 s period of SGR 1935+2154, a weak pulse was noted (S/N $<$ 3.16), but the significance is insufficient for a secure detection of the pulsar-like emission signal. A conservative upper limit on the average monochromatic isotropic luminosity of any possible periodic emission is $2.08\times{10}^{19}~{\rm{erg~s^{-1}~{Hz}^{-1}}}$. Our results offer meaningful low-frequency upper limits on the burst rate of SGR 1935+2154, and hint for very faint pulsar-like radiation at meter wavelengths.
Read the full abstract

All authors Chen-Ran Hu, Jinhuang Cao, S.A. Tyul'bashev, Pei Wang, Yong-Feng Huang, E.A. Brylyakova, G.E. Tyul'basheva, Jin-Jun Geng, Orkash Amat, Ze-Cheng Zou, Chen Du, Nurimangul Nurmamat, Lang Cui, Fan Xu, Xiao-Fei Dong, Chen Deng

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

A Missing Tool for Calculating Auto/Cross-correlation Function under Nonuniform Sampling Observations

Nonuniform sampling presents a long-standing challenge in astrophysical time-domain analysis, invalidating the standard autocorrelation and cross-correlation functions and forcing researchers to adopt ad-hoc methods like interpolation or binning, which introduce unquantified biases and lack rigorous error estimation. Here we introduce a new method for calculating the nonuniform autocorrelation function (NUACF) and nonuniform cross-correlation function (NUCCF) for irregularly sampled time series. Instead of relying on interpolation, it naturally evaluates the correlation function by incorporating time-interval weights and misalignment penalties. Monte Carlo simulations provide confidence bands for significance assessment and a complete error budget for the time delays that accounts for both flux uncertainties and sampling irregularity (essential but generally absent from existing methods). Through extensive simulations, we demonstrate our method outperforms traditional methods across various conditions, from strictly periodic to complex repeating variability patterns (e.g., intermittent but aperiodic). Its effectiveness is demonstrated via various real astrophysical data sets, revealing repetitive variability in stellar light curves, measuring time delays for multi-band disc reverberation in the AGN Fairall 9, and providing model-independent validation of time delays for the gravitationally lensed quasar HE 0435-1223. The method provides a rigorous and general solution to the ubiquitous problem of nonuniform sampling, positioning it as a useful tool for large-scale time-domain survey data analysis. The framework is also directly applicable to emerging time-domain phenomena such as fast radio bursts (FRBs), enabling, e.g., the study of correlations between persistent radio source luminosity and repeating FRB activity, or among the multi-parameter variability curves of FRB emission itself.
Read the full abstract

All authors Chen-Ran Hu, Yong-Feng Huang, Jin-Jun Geng, Orkash Amat, Ze-Cheng Zou, Chen Deng, Fan Xu, Xiao-Fei Dong, Chen Du, Nurimangul Nurmamat, Pei Wang, Lang Cui, Cheng-Ming Li

High Energy Astrophysical Phenomena

Nuclear equation-of-state effects on the two-dimensional post-outburst thermal evolution of magnetized neutron-star crusts

We present a controlled two-dimensional study of nuclear-equation-of-state (EOS) effects on the post-outburst thermal relaxation of magnetized neutron-star crusts. Six EOS models are evolved at a fixed gravitational mass of $1.4 M_{\odot}$ with EOS-specific TOV backgrounds, crust compositions, and transport inputs under identical magnetic-field and heating prescriptions. Our analysis combines self-consistent multi-EOS evolution with a BSk-family factorization that separates structural, microphysical, and interaction contributions, together with an accepted-step energy ledger and a numerical-sensitivity budget. The models show EOS-dependent changes in both the early response and the later redistribution of heat between the surface, crust, and inner boundary; light-curve crossings near $10^2$ days demonstrate that the EOS effect is not a simple luminosity rescaling. By 1000 days the cumulative surface-photon energy fraction differs by less than one percentage point, whereas the internal energy partition differs much more strongly. The updated sensitivity tests show that the late-time luminosity differences and the largest peak contrasts exceed the corresponding BSk24 numerical-sensitivity scales, while the smallest early peak contrast remains less securely resolved. The calculations are intended as a reproducible EOS-sensitivity benchmark rather than an observational fit.
Read the full abstract

All authors Wei-Feng Zhang, Zhi-Fu Gao, Wen-Qi Ma, Na Wang, Yan-ling Liu

High Energy Astrophysical Phenomena

Sharpening the Mass Measurement of PSR J1048+2339 via High-Cadence Photometry

Neutron stars are among the densest objects in the universe, and the mass measurements are crucial to the understanding of their equation of state. Pulsars in tight binaries can be weighed dynamically, but accurate mass estimates are often hampered by uncertainties in orbital inclinations. In this paper, we report on multi-epoch optical monitoring of the pulsar binary, PSR J1048+2339, using the 0.5-m RIFT telescope at Lulin Observatory. Our observations reveal frequent optical flares with eclipses caused by the companion. Using the eclipse duration with archival photometric data and orbital ephemeris in the literature, we constrain the inclination to 78.0+/-1.0 degrees and the pulsar mass to 1.79 (+0.09, -0.08) solar masses. This measurement is among the most precise for any redback system to date, demonstrating a novel way to improve the pulsar mass measures with small aperture optical telescopes. Furthermore, the optical eclipse duration is significantly shorter than that measured in M/GeV gamma-rays, possibly suggesting that a fraction of the eclipsed gamma-ray emission originates from the apex of the intrabinary shock. Based on these constraints, we estimate the magnetic field strength of the companion to be 30-170 G. This is consistent with the high field strength implied by the X-ray orbital modulation observed.
Read the full abstract

All authors Kwan-Lok Li, Ka-Yui Au, Lupin C. C. Lin, C. Y. Hui, Albert K. H. Kong, Jumpei Takata

High Energy Astrophysical Phenomena

Near-Horizon Tidal Disruption Events

Tidal disruption events (TDEs) offer a unique dynamical probe of the spacetime geometry around supermassive black holes (SMBHs). While conventional TDEs occur around SMBHs of $M_\bullet \sim 10^6 - 10^7 M_\odot$, where stars are disrupted far from the event horizon, here we identify a special class of TDEs around rapidly spinning SMBHs with masses $M_\bullet \gtrsim 10^8 M_\odot$, where the tidal-disruption radius approaches the gravitational radius. We term these events "near-horizon TDEs" and, by calculating geodesics in Kerr spacetime, investigate how the proximity of the horizon modifies the debris evolution and subsequent fallback. We find that for stars encountering the SMBH on parabolic orbits, up to $\sim (90 \%-95 \%)$ of the stellar debris either plunges directly into the SMBH or escapes the system. Bound orbits, by contrast, retain a substantially larger fraction of the debris. Using the resulting debris distribution, we calculate the fallback rates and find that bound orbits produce intrinsically higher peak fallback rates and shorter decay timescales, by factors of $\sim 10^3$ relative to parabolic orbits. However, when the stellar orbital angular momentum is particularly low, the peak fallback rate can be substantially suppressed by debris lost to plunge orbits. This combination of rapid fallback and a mass deficit naturally explains overluminous TDEs such as ASASSN-15lh, which standard TDE models have struggled to reproduce. Our work establishes near-horizon TDEs as a new probe of strong-field gravity and a promising tool for identifying massive, rapidly spinning SMBHs.
Read the full abstract

High Energy Astrophysical Phenomena

SN 2025fhm: A central-engine powered Ic-BL supernova associated with X-ray transient EP250304a

We present X-ray, optical, and radio follow-up observations of EP250304a, an extragalactic fast X-ray transient (EFXT) discovered by the Einstein Probe. Its X-ray light curve exhibits two broad pulses with comparable peak fluxes within the first $\sim$1~ks, a feature rarely seen among low-luminosity gamma-ray bursts or EFXTs. Optical follow-up observations were carried out using the Korea Microlensing Telescope Network, the Thai Robotic Telescope, the Las Cumbres Observatory 1~m global network, the Gemini Multi-Object Spectrograph on Gemini south telescope, and the Global Supernova Network. The fast-cooling phase (within 3 days) of optical data can be well fitted by a shocked cocoon model. However, during the supernova phase (SN 2025fhm, from 3 to 88 days), the late-time light curve cannot be explained solely by radioactive $^{56}$Ni decay, as demonstrated by a grid of simulations using the one-dimensional Lagrangian radiation hydrodynamics code SNEC, which reveals a significant energy excess at late epochs. To account for this excess, a central engine like a rapidly spinning, highly magnetized neutron star is needed to provide additional energy injection. This model yields a best-fit spin period of $\sim$12.60~ms and magnetic field strength of $\sim 3.52\times10^{15} \rm G$, and it successfully explains both the late-time bolometric light curve and the early X-ray pulse structures. Our results indicate that EP250304a/SN 2025fhm is likely powered by a central magnetar rather than by radioactive decay alone, offering new insights into the energy budget and physical origin of EFXTs and their associated supernovae.
Read the full abstract

All authors Cui-Ying Song, Nan Jiang, Xiaofeng Wang, Yi-Han Iris Yin, Lingzhi Wang, Wenxiong Li, Shengyu Yan, Dae-Sik Moon, Tao An, Aleksandar Cikota, Samaporn Tinyanont, Liang-Duan Liu, Cui-Yuan Dai, Christopher D. Matzner, Bin-Bin Zhang, Lixin Yu, Qinyu Wu, Hong Soo Park, Sang Chul Kim, Youngdae Lee, Yu-Hao Zhang, Haowei Peng, Franz E. Bauer, Joseph R. Farah, Moira Andrews, Kathryn Wynn, Yuan Qi Ni, D. Andrew Howell, Curtis McCully, Ning-Chen Sun, Danfeng Xiang, Yuan Liu, Wenxin Wang, Yijia Zhang, Wei Chen

High Energy Astrophysical Phenomena

Radio Emission with Dust: A Bow-Shock Interpretation for the TDE Candidate AT 2019avd

Forward shocks produced by interactions between a single ejected blob and the circumnuclear medium (CNM) around a supermassive black hole are widely invoked to explain radio emission from tidal disruption events (TDEs). However, recent observational evidence such as rapidly rising radio emission, double-peaked broadband spectral energy distributions, and declining shock energies at late times poses challenges to this picture. Alternative scenarios have therefore been proposed, including bow shocks arising from outflows colliding with dense clouds, the production of new ejecta at late times, and variations in microphysical parameters. To investigate the roles of these factors, we analyze the long-lived radio flare of the TDE candidate AT 2019avd, for which dust has been identified through mid-infrared observations, using a radio dataset spanning more than five years. We find that a bow-shock scenario can account for the majority of the observed radio emission, while a forward shock, or an additional bow shock involving different clouds, may begin to dominate at late times. However, a forward-shock scenario remains viable when allowing for evolving microphysical parameters and more complex CNM structures. We discuss the remaining tensions in both the outflow-CNM and outflow-cloud interaction scenarios.
Read the full abstract

All authors Shuyuan Wei, Yanan Wang, Guobin Mou, Linhui Wu, Santiago del Palacio, Ranieri D. Baldi, Xiaolong Yang, Yang Huang

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

Low-Redshift Interlopers in the SDSS DR16 $z > 5$ Quasar Catalogue

High-redshift quasars ($z>5$) are crucial for studying early black hole growth and reionisation. SDSS Data Release 16 catalogues 655 such sources with multiple redshift estimators. We systematically examined this subsample and found 480 sources (73.2%) where pipeline redshift (Z\_DR16Q) disagrees with the derived estimate (Z\_FIT) significantly. Visual inspection confirms that 125 of these (19.1% of the total) are actually low-redshift active galactic nuclei. In 114 cases, broad Mg\,ii with C\,iii and C\,iv was misidentified as Ly$\alpha$; in 11 cases, H$\alpha$ was misidentified as Ly$\alpha$. This error is immediately revealed by the absence of Gunn–Peterson trough blueward of Ly$\alpha$. Crucially, Z\_FIT recovers correct redshift in the affected sources where emission lines are prominent, indicating the algorithm is sound; the issue lies in community preference for Z\_DR16Q or Z\_SPEC. We urge caution against uncritical use of pipeline redshifts for $z>5$ DR16Q sources and recommend Z\_FIT as the preferred estimator.
Read the full abstract

High Energy Astrophysical Phenomena

Revisiting EBL Constraints from Gamma-Ray Observations: A Critical Assessment and Methodological Improvements

The extragalactic background light (EBL), ranging from the infrared to the ultraviolet bands, is the second most intense photon field in the universe, surpassed only by the cosmic microwave background (CMB). It is primarily generated by starlight in galaxies, either directly or through absorption by dust and re-emission at longer wavelengths. Very high energy (VHE, $E$ > $100\,$GeV) photons can be absorbed via $\gamma \gamma$ interactions with the EBL during their propagation across cosmological distances, providing an indirect method to probe the EBL by studying its impact on the energy spectra of distant VHE sources. This paper examines the robustness of EBL constraints derived from gamma-ray data, critically assessing the assumptions made in previous studies about the intrinsic source spectra, the uncertainties of the observations and the validity of the analysis tools. We find that earlier studies likely underestimated the uncertainties in the EBL intensity constraints, underscoring the need to account for systematic uncertainties comprehensively. By employing a Monte Carlo (MC) simulation and a plausible model for unknown systematic errors, we compute more realistic uncertainties. Additionally, we discuss possible alternatives to set EBL constraints, relaxing the assumptions on the intrinsic spectra of gamma-ray sources, with the goal of achieving more robust constraints.
Read the full abstract

All authors K. Abe, S. Abe, J. Abhir, A. Abhishek, V. A. Acciari, A. Aguasca-Cabot, I. Agudo, I. Albanese, T. Aniello, S. Ansoldi, A. Arbet Engels, C. Arcaro, T. T. H. Arnesen, A. Babić, C. Bakshi, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Baxter, J. Becerra González, W. Bednarek, E. Bernardini, J. Bernete, A. Berti, J. Besenrieder, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, Ž. Bošnjak, E. Bronzini, I. Burelli, A. Campoy-Ordaz, A. Carosi, R. Carosi, M. Carretero-Castrillo, A. J. Castro-Tirado, D. Cerasole, G. Ceribella, A. Cerviño, Y. Chai, G. Chon, A. Cifuentes Santos, J. L. Contreras, J. Cortina, S. Covino, G. D'Amico, P. Da Vela, F. Dazzi, A. De Angelis, B. De Lotto, M. Delfino, J. Delgado, C. Delgado Mendez, F. Di Pierro, R. Di Tria, L. Di Venere, A. Dinesh, D. Dominis Prester, A. Donini, D. Dorner, M. Doro, L. Eisenberger, D. Elsaesser, L. Foffano, L. Font, F. Frías García-Lago, S. Fröse, Y. Fukazawa, S. García Soto, M. Gaug, J. G. Giesbrecht Paiva, N. Giglietto, F. Giordano, P. Gliwny, T. Gradetzke, R. Grau, J. G. Green, P. Günther, D. Hadasch, A. Hahn, G. Harutyunyan, T. Hassan, J. Herrera Llorente, D. Hrupec, D. Israyelyan, J. Jahanvi, I. Jiménez Martínez, J. Jiménez Quiles, S. Kankkunen, T. Kayanoki, G. W. Kluge, J. Konrad, P. M. Kouch, H. Kubo, J. Kushida, M. Láinez, A. Lamastra, E. Lindfors, S. Lombardi, F. Longo, R. López-Coto, M. López-Moya, A. López-Oramas, S. Loporchio, L. Lulić, E. Lyard, P. Majumdar, M. Makariev, M. Mallamaci, G. Maneva, M. Manganaro, S. Mangano, K. Mannheim, M. Mariotti, M. Martínez, P. Maruševec, D. Mazin, S. Menchiari, J. Méndez Gallego, S. Menon, D. Miceli, J. M. Miranda, R. Mirzoyan, M. Molero González, E. Molina, H. A. Mondal, A. Moralejo, C. Nanci, A. Negro, V. Neustroev, M. Nievas Rosillo, C. Nigro, L. Nikolić, K. Nilsson, S. Nozaki, A. Okumura, J. Otero-Santos, S. Paiano, D. Paneque, R. Paoletti, J. M. Paredes, M. Peresano, M. Persic, M. Pihet, F. Podobnik, P. G. Prada Moroni, E. Prandini, W. Rhode, M. Ribó, J. Rico, A. Roy, N. Sahakyan, F. G. Saturni, F. Schiavone, K. Schmitz, F. Schmuckermaier, A. Sciaccaluga, G. Silvestri, A. Simongini, J. Sitarek, V. Sliusar, D. Sobczynska, A. Stamerra, J. Strišković, D. Strom, M. Strzys, Y. Suda, M. Takahashi, R. Takeishi, J. Tartera Barberà, P. Temnikov, T. Terzić, A. Tutone, S. Ubach, M. Vazquez Acosta, S. Ventura, G. Verna, I. Viale, A. Vigliano, C. F. Vigorito, E. Visentin, V. Vitale, M. Vorbrugg, I. Vovk, R. Walter, C. Walther, F. Wersig