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

From DAMPE to LHAASO: Rigidity Scales and Composition Changes in Galactic Cosmic Rays

This commentary briefly reviews recent measurements from DAMPE and LHAASO on cosmic rays to explore why the spectra of different species change shape across the TeV-PeV energy range. DAMPE finds that several primary nuclei show spectral softening at a similar magnetic rigidity, while the combined proton and helium data from DAMPE and LHAASO reveal a more complex evolution at higher energies. By comparing these spectral features and composition changes, we discuss what they may tell us about the sources of Galactic cosmic rays and how these particles propagate through the Galaxy.
Read the full abstract

Monday 7 September

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

Survey of Novel Deep Learning Architectures for Denoising Gravitational-wave Signals

Gravitational-wave denoising must handle the full diversity of spinning, precessing binaries, since the recovered waveform underpins parameter estimation, tests of general relativity, and population studies. Matched filtering achieves this at a cost that becomes prohibitive as next-generation detectors push event rates higher; deep learning offers real-time reconstruction, but current methods are developed on narrow parameter spaces, precluding principled comparison and reliable deployment. We present the first controlled comparison of five neural-network architectures for gravitational-wave denoising, trained identically across the full astrophysically-motivated spinning binary-black-hole parameter space. A unifying principle emerges: matching network structure to the spectral anatomy of a coalescence – inspiral, merger, ringdown – outperforms brute-force model scaling. Our Multi-Scale Frequency-Aware architecture embodies this via dedicated parallel branches per frequency regime, achieving the best fidelity while using fewer parameters than larger models. It generalizes from simulated training to real LIGO-Virgo-KAGRA data without retraining, recovering merger morphology across confirmed events spanning three observing runs and both detectors, despite training on a single detector's simulated noise. We construct population-level uncertainty bands from denoising residuals, validate their calibration, and stress-test the framework on extended mass ratios and a spin population resembling hierarchical-merger remnants. Applied to real noise with no known signal, the network suppresses its output almost everywhere, with rare exceptions traced to noise transients rather than a general weakness – indicating the statistic discriminates signal from noise and motivating a future detection study. Released weights give a deployable, reproducible benchmark for future extensions.
Read the full abstract

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

The Association of Solar Radio Bursts with Eruptive and Confined Flares

We classify the metric radio emission associated with eruptive (CME-associated) and confined flares using the sample of events between 2010 and 2016 identified by Kazachenko (2023). We find striking differences in the occurrence of radio bursts between the two classes of flare: for soft X-ray flare sizes above M1.4, confined flares largely lack Type II (1% association rate) and Type IV (6%) emission. Approximately 15% of the sample of $\geq$M1.4 confined flares are associated with impulsive-phase Type III bursts. On the other hand, eruptive flares are associated with Type II, Type III (both impulsive and late phase), and Type IV bursts 45-60% of the time. The different types of radio burst are associated with different drivers (IIIs with electron beams, IIs with shocks, IVs with post-flare loops), so it is striking that the connection of all three burst types to eruptive flares is so pronounced. These results can be interpreted in terms of the reconnection topology of the principle candidate flare types, viz., reconnection between closed field lines for confined flares and X-point reconnection in a CSHKP model for eruptive flares, with interchange reconnection for jet-type flares and Type III bursts.
Read the full abstract

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

SN 2024afyu interpreted as a Pair Instability Supernova

Pair-instability supernovae (PISNe) are the predicted explosions of very massive stars triggered by electron-positron pair production. Numerous transients have been proposed as PISN candidates, yet none has provided unambiguous confirmation of this explosion mechanism. The predicted strengths of nebular emission lines offer a powerful means of testing the PISN scenario. We investigate the nature of SN 2024afyu, a nearby (z = 0.0085), long-lived (trise = 85 +- 11.7 days) SN with peculiar spectral evolution, with the aim of identifying its powering mechanism. We analyse multi-band photometry and optical and near-infrared spectroscopy from shortly after explosion to the nebular phase (around 500 days past peak). Besides early appearance of [Ca II] features, we identify a number of sulfur and silicon emission lines, for which we estimate electron temperatures and elemental masses. SN 2024afyu has an inferred 56Ni mass of around 0.4 to 1.0 solar masses and an inferred sulfur mass of the order of 3 solar masses, substantially larger than expected for conventional core-collapse explosions. SN 2024afyu is photometrically similar (although fainter, M_Peak(r) = -18.9 +- 0.04 mag) but spectroscopically distinct to other proposed PISNe. Yet, existing PISN models broadly reproduce several key characteristics, including the overall spectral appearance and broad photometric evolution. SN 2024afyu is a strong PISN candidate, since alternative scenarios would struggle to explain the combination of broad light curve, large intermediate-mass-element abundance, and general spectroscopic evolution. The discrepancies between the observations and currently available theoretical models highlight the need for new PISN calculations spanning a wider range of progenitor masses, metallicities, mixing prescriptions, and circumstellar environments.
Read the full abstract

All authors P. J. Pessi, S. Barmentloo, Ł. Wyrzykowski, S. Schulze, J. Sollerman, A. Gangopadhyay, P. J. Mikołajczyk, S. Rose, K. Kotysz, C. Fremling, A. Jerkstrand, R. Lunnan, L. Yan, C. Fransson, J. L. Wise, M. Dubey, K. Misra, A. P. Cotter, D. Athanasopoulos, A. Bochenek, M. Bronikowski, U. Burgaz, W. Burzynski, J. M. Carrasco, E. Concepcion, C. Galdies, A. F. Gillan, A. Gkini, E. Gkiokas, E. Głowacki, F-J. Hambsch, J. Japelj, N. Karaman, P. King, S. Kurowski, C. Liu, M. Maskoliunas, Z. McGrath, M. Mihelčič, M. Nikolajuk, W. Ogloza, E. Pakštienė, U. Pylypenko, A. Pucek, N. Rehemtulla, R. M. Rich, A. O. Simon, A. Singh, C. Skoglund, E. Stonkute, Y. Stsefanenka, P. Szober, K. Tsalapatas, B. F. A. van Baal, C. Ventura, K. Vrontaki, A. Wozniak, J. Zdanavicius, S. Zola, E. C. Bellm, J. Castaneda Jaimes, M. J. Graham, G. Helou, F. J. Masci, J. N. Purdum, R. Riddle

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

Nonradial oscillations of realistic anisotropic neutron stars: Polar modes

In this work, we study the polar perturbations of static, spherically symmetric neutron stars with anisotropic pressure in full general relativity, including linear-order perturbations of both the metric and the fluid. We calculate the $f$-mode frequencies and the corresponding damping times using a consistent treatment of the perturbation of the radial vector $k^\alpha$. In particular, its Lagrangian perturbation $\Delta k^\alpha$ is determined by the metric perturbations and the fluid Lagrangian displacement and is constrained to the $(\tilde{u},\tilde{k})$ plane, where $\tilde{u}^\alpha$ is the normalized fluid four-velocity. This constraint introduces an additional dynamical degree of freedom into the perturbation equations. Considering three equations of state and the Horvat and Bowers-Liang prescriptions for pressure anisotropy, we find that the $f$-mode frequency increases with stellar mass, ranging from $1$ to $3$~kHz, while the damping time decreases, ranging from $0.5$ to $1.25$~s. Increasing anisotropy, in the sense of tangential pressure exceeding radial pressure, generally lowers the oscillation frequency, while its effect on the damping time depends on the anisotropy prescription: the damping time decreases with increasing anisotropy for the Horvat model but increases with increasing anisotropy for the Bowers-Liang model. We further find quasi-universal relations between the real and imaginary parts of the $f$-mode frequency, $M\omega_R$ and $M\omega_I$, and the stellar compactness $\mathcal{C}=M/R$, which are largely insensitive to the equation of state. Polynomial fits to these relations achieve an accuracy better than $10\%$, providing a simple phenomenological framework for constraining neutron-star pressure anisotropy through future asteroseismology observations.
Read the full abstract

All authors L. M. Becerra, José F. Rodríguez-Ruiz, E.A. Becerra-Vergara, F. D. Lora-Clavijo

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

A pilot submillimeter search for IceCube neutrino counterparts: JCMT follow-up and dusty-galaxy catalog associations

The submillimeter band offers a complementary route to identifying electromagnetic counterparts of high-energy neutrinos, tracing relativistic transients and persistent dust-obscured galaxies. We present a pilot search combining JCMT/SCUBA-2 follow-up of IceCube events with catalog searches for bright high-redshift dusty systems. Among nine rapidly observed IceCube fields, sources of particular interest were found in three: JCMT0402-0424 ("Shadow Blaster") in IC 210922A, a transient-like source following IC 220115A, and two known blazars in IC 240105A. ALMA resolved the Shadow Blaster into a strongly lensed dusty star-forming galaxy at z=2.988 with a compact gas-rich starburst core. Motivated by this result, we cross-matched IceCube localizations with the Herschel Bright Sources (HerBS) and Planck High-Redshift Source Candidate (PHZ) catalogs. Among 51 historical IceCube events with A_bbox,90 <= pi deg^2, four fields contain at least one HerBS source, compared with 0.723 expected from right-ascension scrambling (fixed-cut p=4.18e-3). Scanning the localization-area threshold gives p_global,A=0.0135. Applying the same IceCube selection to PHZ gives four associated fields versus 3.544 expected (p=0.478). Because the focused HerBS selection followed a broader exploratory analysis with trials-corrected p_global=0.130, the focused result is suggestive rather than definitive. The four well-localized HerBS fields contain six high-redshift DSFGs with apparent L_IR~10^13 L_sun. ALMA imaging and published spectroscopy reveal compact or multiple-component morphologies, extreme molecular-gas kinematics, and possible obscured nuclear activity. These results support a hierarchical strategy: identify luminous submillimeter candidates through wide-field imaging and catalogs, then characterize them at high angular resolution to test for compact gas-rich calorimetric conditions and obscured AGN activity.
Read the full abstract

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

An EOS-Driven Extension of NSCool for Compact Star Cooling with Hadronic and Quark Degrees of Freedom

We present an EOS-driven extension of the NSCool thermal evolution code that enables complete tabulated equations of state to be treated within a single composition-based input structure. The generalized NEW interface includes additional baryonic fractions, a bosonic composition variable, the hadronic volume fraction, and quark fractions, allowing nucleonic, hyperonic, resonant-baryonic and hybrid hadron-quark configurations to be handled within the same workflow. The neutrino sector is extended accordingly, including updated nucleonic modified Urca and bremsstrahlung rates, additional baryonic direct Urca channels, hyperonic processes, baryonic pair breaking and formation and quark direct Urca, modified Urca, bremsstrahlung and PBF contributions. In mixed phases, emissivities and the corresponding core heat capacity contributions are evaluated from phase-local quantities and combined using the hadronic volume fraction. The implementation is validated against the original NSCool calculation for a controlled nucleonic benchmark and demonstrated with representative hadronic and quark-containing EOSs. These calculations are intended as software validation tests of the generalized workflow rather than as observational fits or statistical EOS inference.
Read the full abstract

High Energy Astrophysical Phenomena

Abundance Stratification in Type Ia Supernovae – VIII. The 03fg-like SN\,2012dn: a double-degenerate merger interpretation

A detailed study using abundance tomography of SN 2012dn, a peculiar Type Ia supernova, is presented. Despite exhibiting a normal peak luminosity, it retains early-phase super-Chandrasekhar (03fg-like) characteristics: weak Fe III, narrow intermediate-mass element (IME), persistent carbon, and the absence of high-velocity Ca II features. Its nebular spectrum is unusually faint due to grey dimming beginning $\sim$60 days after maximum and shows [O I] emission, rare in SNe Ia. While a Chandrasekhar-mass density profile reproduces the photospheric phase, it fails at late times. The [O I] emission requires additional low-velocity mass, implying a total ejecta mass of $1.66$ $M_\odot$, including $0.33$ $M_\odot$ of oxygen (with $0.1$ $M_\odot$ in the core). The $^{56}$Ni mass is estimated to be $0.45$–$0.49$ $M_\odot$, a range set by uncertainties in the grey-extinction correction applied at the nebular phase, and is insufficient on its own to account for the peak luminosity. Stable iron is confined to intermediate layers and absent from both the core and outer ejecta, the latter indicating sub-solar progenitor metallicity. Silicon and sulfur span the full ejecta, while carbon extends down to $v \sim 6000$ ${\rm km,s^{-1}}$. These properties favour a double-degenerate CO–CO white dwarf merger scenario. The remaining luminosity deficit may be explained by additional energy input from a weak interaction with a low-mass, carbon-rich circumstellar shell.
Read the full abstract

All authors Charles J. Aouad, Paolo A. Mazzali, Stephan Hachinger, Chris Ashall

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

A New Sample of $\sim$ 100 Intermediate-mass Black Holes Reaching $z \approx 1$

We present a systematic search for intermediate-mass black hole (IMBH) active galactic nuclei (AGNs) at $0.5 < z \lesssim 1$ using DESI DR1 spectroscopy.We identify 98 broad-line IMBH AGNs with black hole masses $M_\mathrm{BH}<10^6$ $M_{\odot}$ through quantitative spectral decomposition and broad-H$\beta$ selection. This sample spans $M_\mathrm{BH}=10^{5.5}-10^{6.0}$ $M_{\odot}$ and Eddington ratios from 1.2 to 9.0 extending systematic IMBH AGN searches to intermediate redshift. Compared with a consistently selected $z<0.6$ IMBH sample, the $0.5<z<1$ sources reach higher broad-H$\beta$ luminosities and substantially higher Eddington ratios, indicating more extreme accretion states at earlier cosmic times. They also show broader [O III] profiles and stronger blueshifted wing components, with these kinematic differences persisting after matching in black hole mass and Eddington ratio. These results reveal two distinct signatures of evolution among IMBH AGNs at $z<1$: ability of IMBHs to reach increasingly extreme accretion states toward higher redshift, and systematic changes in their ionized-gas kinematics. The former demonstrates that rapid, including super-Eddington, growth of IMBHs can persist to relatively late cosmic times, while the latter may indicate evolution in the ionized-gas environment and associated outflow activity of actively growing IMBHs. Together, these findings provide new constraints on the evolutionary pathways of IMBHs at $z<1$, and show that seed-mass black holes can continue to undergo rapid growth well after the cosmic dawn.
Read the full abstract

All authors Wen-Juan Liu, Luis C. Ho, Su Yao, Xiao-Bo Dong, Yaqi Zhao

High Energy Astrophysical Phenomena

Two-Phase Structure of Synchrotron-Cooling-Unstable Relativistic Plasma

Using analytic theory, radiative particle-in-cell (PIC) simulations, and fluid simulations, we show that relativistic, synchrotron-cooling, collisionless, high-$\beta$ pair plasmas filament into a two-phase medium. This process occurs through the interplay of the synchrotron cooling instability (SCI) with the synchrotron firehose instability (SFHI). One phase has high plasma~$\beta$ and is infested with small-scale firehose fluctuations, which scatter particles and pin the pressure anisotropy to the firehose-marginal level. The other phase has much lower~$\beta$, causing the suppression of firehose modes and thus allowing large pressure anisotropies. We propose a fluid model for this two-phase plasma, which we use to study the linear and nonlinear evolution of the SCI and SFHI, and to predict the emergence time of the two-phase structure.
Read the full abstract

All authors Agnieszka Wierzchucka, Pablo J. Bilbao, Robert J. Ewart, Dmitri A. Uzdensky, Alexander A. Schekochihin

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

Search for Long-Transient Gravitational Waves from Supernova SN2023ixf using GFH-v2 Pipeline

We present a directed search for long-transient gravitational waves from the possible newborn magnetar remnant of SN2023ixf, a nearby Type II core-collapse supernova in the M101 galaxy. The analysis uses LIGO Hanford and Livingston data from Engineering Run 15, using coincident data lying within the on-source window associated with the supernova. We target signals from a rapidly rotating, non-axisymmetric neutron star whose spin-down is dominated by gravitational-wave emission, producing a power-law decrease in frequency and a corresponding decrease in strain amplitude. The search is performed with the GFH-v2 pipeline, based on the Generalized Frequency Hough transform. No candidate survives the coincidence and follow-up analysis. We therefore set upper limits on the maximum detectable distance as a function of initial frequency and ellipticity. For the highest ellipticity interval, the 90% upper limits reach distances of about 1-2.5~Mpc across most of the analysed band. Although these limits are below the distance to M101, the search provides the first application of GFH-v2 to a nearby core-collapse supernova and characterizes its performance on real detector data.
Read the full abstract

All authors Sandhya Sajith Menon, Lorenzo Pierini, Pia Astone, Cristiano Palomba, Francesco Safai Tehrani, Lorenzo Silvestri, Ornella Juliana Piccinni, Simone Dall'Osso, Stefano Dal Pra, Sabrina D'Antonio, Sergio Frasca, Dafne Guetta, Paola Leaci, Michal Bejger, Alicia Calafat, Evan Goetz, Rafel Jaume Amengual, Ian Jones, David Keitel, Andrzej Krolak, Iuri La Rosa, Andrew Melatos, Joan-Rene Merou, Lorenzo Mirasola, Claudio Salvadore, Alicia M. Sintes, Karl Wette

High Energy Astrophysical Phenomena

Bridging the Population Synthesis of Supermassive Binary Black Holes and the Gravitational Wave Background

Pulsar Timing Arrays (PTAs) constrain population properties of supermassive binary black holes (SMBHBs) through the observation of the gravitational wave background (GWB). Unlike other approaches that interpolate population-synthesis libraries or only consider the mean of the strain spectrum, here we capture its full strain probability density directly from semi-analytic population models. We apply our new method to the semi-analytic SMBHB population model, independently reproducing the parameter estimation for this model performed by the NANOGrav Collaboration with their 15-yr data. We also show the extent to which discrete SMBHB contributions to the GWB resolve degeneracies in the population parameter space. Finally, using the source-count intensity as the intermediate product in our calculation, we map PTA observations, as a proof of principle, to the SMBHB model based on the galaxy merger prescriptions from numerical hydrodynamical simulations "Illustris". We find the effect of delay times $\tau$ between kiloparsec and subparsec SMBHB separations following galaxy mergers, finding $h_{\rm c}$ spanning $(1$-$6)\times10^{-16}$ and $N_{\rm c}$ spanning $(0.3$-$7.1)\times10^{-3}$ for $\tau$ up to 8 Gyr.
Read the full abstract

All authors Kanyuni Iemoto, Boris Goncharov, Gabriela Sato-Polito, Xiaoming Bi

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

Approximating the statistics of a gravitational wave background

The astrophysical origin of the gravitational wave background (GWB) reported by pulsar timing array (PTA) collaborations has yet to be confirmed. A GWB signal made of the sum of individual gravitational wave (GW) from the population of supermassive black hole binaries (SMBHB) would show the imprint of a discrete Poissonian statistics in its spectral properties. In this work, we propose a tool based on the saddlepoint approximation method to estimate the distribution of characteristic strain for any given population model. This tool can be used for Bayesian inference or for a quick visualization of the statistics of the GWB from the output of more realistic semi-analytical models. We introduce the general family of variance mixture Gaussian distributions that models heavy-tailed behavior distributions that is expected for a non-Gaussian GWB signal. We show that setting the correct hierarchical priors to a Gaussian free spectrum PTA likelihood is effectively equivalent to a non-Gaussian likelihood. Using ideal simulations, we compare the performance of the saddlepoint approximation with a log-Normal distribution to infer the parameters of the astrophysical model from the statistics of the GWB and show that correctly modeling higher order moments is essential. Future PTA analyses should include the statistics of the GWB in their pipelines.
Read the full abstract

High Energy Astrophysical Phenomena

Interior Magnetic Fields in Magnetars and Radio Pulsars

Magnetic fields are fundamental for neutron star physics and play a central role in powering the extreme phenomenology of magnetars, including Soft Gamma Repeaters and Anomalous X-ray Pulsars. However, the structure and stability of their internal magnetic fields remain largely unconstrained, as they cannot be directly probed by electromagnetic observations. Using 3D general-relativistic magnetohydrodynamics simulations across a range of rotation rates and magnetic field strengths, we identify two distinct evolutionary regimes leading towards dynamically stable magnetic configurations. In rapidly rotating stars, the Alfvén crossing timescale exceeds the rotation period, allowing differential winding to amplify a strong toroidal magnetic component before the onset of instabilities, leading to long-lived, stable configurations. In contrast, in magnetically dominated stars, instabilities in the poloidal field drive rapid field decay, leaving only a comparatively weak toroidal component. These results imply that the internal magnetic structure of neutron stars depends sensitively on their rotational state: rotation-dominated stars like radio pulsars develop strong toroidal fields, while magnetars are characterized by predominantly poloidal configurations. Our findings therefore show that a neutron star's rotational history shapes its internal magnetic structure, providing a unifying physical picture that connects the observed diversity of neutron star classes to their hidden field configurations.
Read the full abstract

All authors Raj Kishor Joshi, Brynmor Haskell, William Cook, Sebastiano Bernuzzi

High Energy Astrophysical Phenomena

A cosmic-ray database update: CRDB v4.2

The cosmic-ray database, CRDB, has provided data and meta-data to the community since 2013. We present a new release, CRDB v4.2, which includes data and time series published in the last three years (from AMS-02, CALET, DAMPE, LHAASO, etc.), datasets missed in previous releases, and all peer-reviewed ultra-high energy cosmic-ray data (both historical and recent datasets). These data and meta-data are available from the web interface or from the CRDB pip-installable Python library. This release includes results from more than 60 experiments and publications, corresponding to 4534 new data points and 131624 time-series data points.
Read the full abstract

High Energy Astrophysical Phenomena

Mining the time axis with TRON - III. Discovery of a nulling pulsar towards the SMC with MeerKAT

We report on the discovery of a nulling pulsar, PSR J0052$-$7551, in archival MeerKAT L-band observations of the Small Magellanic Cloud (SMC). The discovery was made using the image-plane transient detection pipeline, TRON, which creates a residual image for each time sample and uses a heuristic pruning strategy to identify promising candidates in the resulting temporal image cube. The TRON pipeline triggered on the approximately minute timescale variability associated with the nulling behaviour of a canonical pulsar. Follow-up MeerKAT UHF-band beamformer observations revealed a pulsar with a spin period of 0.55 s and a measured period derivative of 6.1(6)$\times$10$^{-17}$s$\cdot$s$^{-1}$. The source is seen in projection toward the periphery of the SMC, and has a dispersion measure value of 29.24(2) pc$\cdot$cm$^{-3}$, which is on the order of the expected Milky Way contribution along the line of sight. We use a Gaussian Mixture Model to fit the nulling behaviour which yielded a nulling fraction of $\sim 32$%. We find that the null lengths are consistent with an exponential distribution indicating a stochastic process. Additionally, power spectra of the UHF visibility data reveal an unexplained quasi-periodic feature at ~5$-$7 s that cannot be attributed to the source itself, RFI, or standard instrumental effects. Following on from previous image-plane detections of pulsars by the TRON approach, which relied on either scintillation or eclipsing behaviour, this is the first detection triggered by nulling. This discovery confirms TRON's sensitivity to a wide range of astrophysical phenomena of medium timescale duration.
Read the full abstract

All authors T. Blecher, M. Geyer, O. M. Smirnov, A. Coelen, I. Heywood, I. Pelisoli, A. K. Hughes, J. S. Bright, V. Prayag, T. Myburgh, C. Tasse, B. Hugo, J. S. Kenyon, H. L. Bester, S. J. Perkins, J. Dawson, V. G. Samboco, A. J. T. Ramaila, B. Ngcebetsha, I. Sihlangu, C. Choza, V. S. Dhillon, N. Castro Segura, J. McCormac

High Energy Astrophysical Phenomena

Oblique Collision of a Relativistic Cold Shell with an Ideal Reflecting Wall

Relativistic flows are common in astrophysics and often form shocks when different parts of the flow collide at relativistic relative velocities. Such collisions are often oblique, forming two shocks whose shocked fluids are separated by a contact discontinuity, which is treated here as an ideal reflecting “wall” where the flow on either side is modeled separately. The latter is modeled in the lab frame $S$ as a uniform cold planar shell propagating into vacuum at velocity $v_1=\beta_1c$ normal to its vacuum interface, colliding with the wall at an incidence angle $\alpha_1$. The collision point $P$ moves along the wall at a velocity $v_p=v_{1}/\sin\alpha_1$, and a boost along the wall at $v_p$ leads to a steady-state frame $S'$ where this problem is highly simplified. However, a “super-luminal” regime exists where $v_p>c\Leftrightarrow\tan\alpha_1<\Gamma_{1}\beta_{1}=(1-\beta_{1}^2)^{-1/2}\beta_1$ and no steady-state frame $S'$ exists. It corresponds to only very small $\alpha_1$ in the Newtonian regime, but nearly all $\alpha_1$ in the relativistic regime. We solve this problem \textit{fully analytically} using integral conservation laws, in the attachmrnt region where point $P$ is attached to the wall. This region of parameter space is bound at high $\alpha_1$ by the detachment line, which coincides with the sonic line for a cold initial shell. A weak-shock solution exist in all this region, while a strong-shock solution exists only in the sub-luminal attachment region – between the luminal line and the detachment/sonic line where the two solutions coincide and beyond which point $P$ detaches from the wall and shocked fluid spills into the vacuum.
Read the full abstract

All authors Jonathan Granot, Prasanta Bera, Michael Rabinovich, Paz Beniamini

High Energy Astrophysical Phenomena

Timing Gravity with Pulsars in the Strong Field

We propose a novel approach for the timing of pulsars orbiting a supermassive black hole, which implements the fully relativistic calculations of the photon travel time into a robust timing model. We generate realistic mock catalogues of pulsar times-of-arrival for several putative pulsars on tight orbits around the Galactic Center supermassive black hole, Sagittarius A* (Sgr A*). Then, we perform a proof-of-concept sensitivity analysis to forecast the accuracy that future observational facilities, like the Squared Kilometer Array, will achieve in the characterization of the parameters of our timing model. Our analysis shows how the observation of pulsars at the Galactic Center will open an incredibly promising avenue for the characterization of the physical properties of Sgr A*, which can improve by at least three orders of magnitude the current constraints on the black hole's mass achieved with the S-stars and event-horizon scale observations.
Read the full abstract

High Energy Astrophysical Phenomena

Correlation between X-ray and gamma data of Swift measurements

Several studies over the last two decades have used canonical correlation analysis (CCA) to study the relationships between main {\gamma}-ray (e.g. fluence, peak flux, and duration) and main X-ray (flux, decay and spectral index, and hydrogen column density) data from gamma-ray bursts (GRBs). In this paper, we revisit this approach using a much larger dataset to identify potential new insights into these relationships. We used CCA to investigate the interrelationship of the aforementioned gamma-ray and X-ray parameters. Using the derived canonical variables, we calculated their correlations (canonical loadings) with the original data. Consistently with previous research, the analysis revealed that gamma-ray fluence and X-ray flux have the strongest correlation, while the X-ray decay index and spectral index have a lower contribution. Interestingly, our analysis of a much larger dataset reveals that the HI column density makes a significant contribution to the overall correlation. This finding, in the context of the collapsar model for long GRBs, could be interpreted as an indication that the progenitor star ejected an HI envelope during the GRB.
Read the full abstract

All authors Istvan I. Racz, Lajos G. Balazs, Istvan Horvath, Sandor Pinter

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

Exact Cancellation of Horizon Modes in Kerr Ringdown

Ringdown gravitational waves probe black hole spacetimes via quasinormal modes (QNMs). Recent studies have reported oscillatory features in QNM-filtered waveforms, and have led to conflicting interpretations as to whether they represent a near-horizon "direct wave". Reexamining a particle plunging into a Kerr black hole, we show that the saddle-point approximation underlying the proposed interpretation fails near the horizon. Our exact analysis reveals that every source-induced horizon-mode pole is canceled by an infinite tower of Matsubara zeros, independently of orbital and black-hole parameters. The first four stages of this sequential screening are also confirmed by time-domain numerical integration. Thus, within the near-horizon perturbative framework, source-induced horizon-mode signals do not survive at late times.
Read the full abstract

All authors Kei-ichiro Kubota, Yusuke Manita, Hayato Motohashi, Daiki Watarai

High Energy Astrophysical Phenomena

LOFAR follow-up of sources in the First LHAASO Catalogue

The First LHAASO Catalogue lists 90 sources of very-high-energy gamma-ray emission, many without an established counterpart at other frequencies, reopening the question of which Galactic accelerators produce the highest-energy cosmic rays. We search for low-frequency radio counterparts to the 55 1LHAASO sources that fall within the Galactic fields of the LOFAR Two-metre Sky Survey Data Release Three (LoTSS-DR3), using 144~MHz continuum maps at $6''$ and $20''$ resolution, complemented with archival data from other radio surveys where available. We present an overview of the radio emission towards all 55 sources, and identify the structures potentially associated with the gamma-ray emission. We report 24 new supernova remnant (SNR) candidates, nine flat-spectrum radio shells with faint or no infrared counterpart, and six sources with no clear Galactic radio emission present in the LOFAR maps. A bootstrap analysis shows that 23 of the 50 sources in the Galactic plane footprint (46%, for $|b| \leq 5$) overlap a catalogued SNR, whereas random placement yields only $10.0\pm2.8$ matches even if the underlying Galactic SNR population is ten times larger than currently catalogued; in the outer Galaxy, where chance alignment is rarest, 4 of 12 sources coincide with an SNR against a null expectation of $\lesssim2$. We estimate that at least a quarter, and up to half, of the 1LHAASO Galactic-plane sources genuinely reside in SNR environments, consistent with a population dominated by pulsar wind nebulae, SNR–PWN composites, and cosmic-ray-illuminated clouds. Individual results include a new SNR candidate inside the TeV shell HESS~J1912+101, a radio candidate for the SNR hypothesised to explain the GeV emission of 1LHAASO~J1945+2424, and the first proposed counterparts of 1LHAASO~JJ0056+6343u and 1LHAASO~J2200+5643u.
Read the full abstract

All authors Maria Arias, Timothy Shimwell, Martin Hardcastle, Roland Timmerman

High Energy Astrophysical Phenomena

NEXUS: Transient Searches and First Results from Year One Observations

We describe ongoing efforts of high-redshift transient searches using multi-epoch NIRCam imaging (F200W+F444W) and NIRSpec MSA/PRISM spectroscopy from the NEXUS JWST multi-cycle Treasury program targeting the north ecliptic pole region. The transient search area covers $\sim 61.5\,{\rm arcmin^2}$ between the reference epoch and each subsequent NEXUS-Deep epoch at a cadence of $\sim2$~months. In the first year of observations from NEXUS, we detect 68 robust transients, with the host photometric redshift distribution declining rapidly at $z>2$ but extending to $z_{\rm phot}\approx 6$. In addition, we obtained secure spectroscopic redshifts for 37 transients ($\sim54\%$) from NIRSpec/PRISM and NIRCam/WFSS, with seventeen at $1 < z < 2$, eight at $2 < z < 3$, two at $3 < z < 4$, and one tentative host association at z = 6.151. Overall, the NEXUS program recovers observed supernova (SN) rates broadly consistent with other SN search programs with JWST. While NEXUS achieves the highest transient detection efficiency, 6.6 SNe per imaging hr ($2.6\times$ COSMOS-SN and $26\times$ JADES-SN), the limited filter coverage (F200W+F444W only) limits robust identification and classification of high-$z$ SNe for deep spectroscopic follow-up. We describe the details of the data reduction and transient detection pipeline, and report the Year-1 transient sample along with their light curves, available MSA spectroscopy, host association, and the raw detection rate. We also describe and release a new PSF photometry package that properly accounts for correlated pixel noise from combining drizzled images.
Read the full abstract

All authors Ming-Yang Zhuang, Lei Hu, Justin D. R. Pierel, Zhiwei Pan, Yue Shen, Adam J. Burgasser, Xiaohui Fan, Luis C. Ho, Junyao Li, Alice E. Shapley, Zachary Stone, Feige Wang, Lifan Wang

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

A New Window for Testing the PPN in the Strong-Field Regime: Stellar Rotation of S301

The recently discovered S301 star, with an orbital period of 8.7 years and eccentricity $e = 0.982$, reaches a pericenter of just $140\,R_{\mathrm{s}}$ — approximately nine times closer than S2's $\sim 1200R_{\mathrm{s}}$ — offering an unprecedented laboratory for testing strong-field gravity. We demonstrate that S301's stellar spin precession provides a novel probe of the PPN parameter $\gamma$, scaling as $(2\gamma_{\mathrm{PPN}} + 1)/3$ relative to general relativity (GR). For maximum projected velocity shifts of $46.1~\mathrm{km\,s^{-1}}$, next-generation spectrographs could constrain $\gamma_{\mathrm{PPN}}$ to $\sim 10\%$ precision. When combined with S301's orbital precession — a substantial $\sim 1.95^\circ$ pericenter advance per orbit (GR prediction), approximately an order of magnitude larger than S2's $\sim 0.2^\circ$, constraining $(2\beta_{\mathrm{PPN}} + 2\gamma_{\mathrm{PPN}} - 1)/3$ — the spin measurement breaks the $\beta$-$\gamma$ degeneracy. A MCMC analysis of projected S2 and S301 observations yields $\sigma_\beta \approx 2.3 \times 10^{-4}$, with the $\gamma$ precision limited by the spin precession measurement to $\sigma_\gamma \approx 0.1$. While the spin precession provides an independent consistency check, its true significance lies in probing gravity at $\phi/c^2 \sim 10^{-4}$ and $v/c \sim 0.08$ — four orders of magnitude stronger than Solar System tests and an order of magnitude stronger than current S-star constraints. This complementary approach opens a new window for testing the PPN framework in the strong-field regime and demonstrates the power of combining multiple stellar probes at the Galactic Center.
Read the full abstract

Sunday 6 September

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

Multimessenger Spectroscopy

A multimessenger method is developed for scheduling optical spectroscopy of ultracompact double white dwarf binaries from Laser Interferometer Space Antenna (LISA) data. The orbital parameters inferred from the gravitational wave (GW) signal are used to predict the argument of latitude $u(t)$, which sets the phase dependence of the radial velocity. With $u(t)$ known in advance, spectra can be taken at selected orbital phases. The GW and spectroscopic descriptions are linked by $u(t)=\Phi_{\rm GW}(t)/2+\delta_{e,\omega}(t)+\mathcal{O}(e^2)$, where $\Phi_{\rm GW}$ is the GW phase, $e$ the eccentricity, and $\delta_{e,\omega}$ the equation of the centre to $\mathcal{O}(e)$, which vanishes for a circular orbit. Radial velocities are regressed on $\cos u$; for double-lined binaries, the two slopes determine the semi-amplitudes $K_1$ and $K_2$. With the inclination and orbital period measured by LISA, both component masses follow from the semi-amplitudes through Kepler's third law, and the resulting chirp mass with the GW amplitude gives the source distance. We derive waveforms to first order in the Laplace–Lagrange parameters; the circular quadrupole waveform is recovered as $e\to0$. Using simulated LISA data for one circular and one eccentric system, we forecast how long $u(t)$ remains accurate enough for phase-resolved spectroscopy. Masses derived from the velocity amplitudes are free of the biases that tides, mass transfer, or acceleration along the line of sight introduce when $\dot f_{\rm GW}$ alone is used. For a four-year LISA observation, scheduling at quadrature remains viable for roughly a decade after the observation ends. Four years after the observation midpoint, the ephemeris predicts the orbital phase to $0.3$~s for the eccentric system and $0.9$~s for the circular one, compared with $1$–$30$~s for optical ephemerides propagated over the same baseline.
Read the full abstract

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

Solar Capture Tests of Inelastic Dark Matter after the LZ High-Recoil Event

The LUX-ZEPLIN (LZ) Collaboration has reported a $248$ keV nuclear-recoil candidate for which endothermic dark matter (DM) gives some of the largest local significances. We study Solar-capture constraints on three interpretations: a thermal Higgsino, a thermal pseudo-Dirac fermion with off-diagonal vector interactions, and neutron-philic endothermic spin-dependent scattering through $\mathcal{O}_4$. The canonical full-density thermal Higgsino is excluded: its LZ-preferred splitting near $0.38\,\mathrm{MeV}$ lies well below the splitting required to suppress Solar capture sufficiently to satisfy the IceCube upper limits on DM annihilation in the Sun, $\delta\gtrsim0.51$-$0.56\,\mathrm{MeV}$. For the pseudo-Dirac benchmark parameters that fit the LZ event, we find $C_\odot=1.53\times10^{20}\,\mathrm{s}^{-1}$, while two-state kinetics limits the fixed-orbit annihilation rate to $\Gamma_A\lesssim4.4\times10^{14}\,\mathrm{s}^{-1}$, about $1.7\times10^5$ below the IceCube upper limit for the $b\bar b$ channel. A semi-analytic treatment indicates that re-excitation cycles further cool the captured population, although a full phase-space calculation is required for its final distribution. For neutron-philic $\mathcal{O}_4$ scattering at $m_\chi=1$ TeV and $\delta=300$ keV, finite-temperature capture gives $C_\odot=3.98\times10^{17}\,\mathrm{s}^{-1}$, with nuclear-structure uncertainties giving an envelope $(1.77$-$9.72)\times10^{17}\,\mathrm{s}^{-1}$. Under the equilibrium assumption, the upper edge remains about a factor $154$ below the IceCube upper limit for the $b\bar b$ channel. We use $b\bar b$ only as a soft-hadronic proxy because an exact constraint requires the model-specific annihilation spectrum. Solar capture therefore excludes the thermal-Higgsino interpretation but not endothermic explanations generically.
Read the full abstract

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

Physics-based phenomenological modeling of binary black hole hierarchical formation 2: Autodifferentiable functional inference of hierarchical compact-binary populations

The gravitational-wave (GW) census contains mass and spin structure consistent with contributions from black holes assembled through repeated mergers in dense environments. Connecting that structure to formation physics requires models that are both physically interpretable and tractable within population inference. We construct an autodifferentiable, physics-based phenomenological model in which each dense environment is represented by a coagulation response and a population of such environments produces an observable merger-rate density. Embedded in the gwkokab Poisson-likelihood framework, this model enables joint inference of natal-population and interaction parameters from the GW census. Applied to GWTC-5.0, the framework shows why simple pairwise coagulation models struggle to reproduce the observed high-mass, comparable-mass population and tests alternative interaction structures against the data, while retaining an explicitly modeled natal component.
Read the full abstract

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

A Common $\sim$55 GeV Inelastic Dark Matter Origin of the LZ, Fermi-LAT, and AMS-02 Excesses and Paleo-Detector Prospects

We investigate the 248~keV nuclear recoil excess reported by LUX-ZEPLIN (LZ) in a leptophobic inelastic dark matter framework, in which the mass splitting naturally suppresses low-energy recoils. We find two viable branches, including the endothermic solution with a high dark matter mass $m_\chi \gtrsim400$~GeV and the exothermic one with a mass of tens GeV. Remarkably, the exothermic dark matter scenario with $m_\chi\sim55$~GeV can simultaneously account for the LZ excess, the Galactic Center GeV gamma-ray excess, and the AMS-02 antiproton excess. The favored annihilation cross-section is slightly below the one required for the observed relic abundance, resulting in an overabundant relic in the minimal setup. An additional neutrino-portal mediator can provide the required depletion while leaving the gamma-ray and antiproton signals essentially unchanged. We further show that paleo-detectors provide a powerful probe of the inelastic interpretation, requiring exposures of $\sim10$ and $\sim100$~mg Gyr for the endothermic and exothermic scenarios, respectively.
Read the full abstract

All authors Meiwen Yang, Quan-feng Wu, Kai-Kai Duan, Zi-Qing Xia, Yue-Lin Sming Tsai, Yi-Zhong Fan

High Energy Astrophysical Phenomena

Possible Sources of Iron Nuclei in Ultra-High-Energy Cosmic Rays

The origin of a cosmic ray particle with an ultra-high energy of 2.44 x 10^20 eV recorded by the Telescope Array is discussed by assuming that it is a heavy nucleus, namely an iron or silicon nucleus. Iron is chosen as a heaviest element in cosmic rays, while silicon is chosen to find out how much the theoretical results obtained vary with the mass of the heavy nucleus. Assuming extragalactic origin of the particle with such energy, the constraints on the distances to possible sources are analysed using a routine calculation of the cosmic ray propagation in extragalactic space. These possible sources are discussed.
Read the full abstract

High Energy Astrophysical Phenomena

CoportSL: A Contribution-constrained Hybrid Slow-light Framework for Time-dependent Polarized GRMHD Imaging

Fast-light approximations neglect fluid evolution along rays, whereas slow-light modeling is indispensable for recovering the true magnetohydrodynamic state. However, full slow-light radiative transfer for extended general relativistic magnetohydrodynamic (GRMHD) sources requires simultaneous access to many fluid snapshots and is memory-intensive. We introduce CoportSL, the first contribution-constrained hybrid slow-light framework for time-dependent full-Stokes imaging. It uses emission, absorption, and Faraday contributions to identify where fluid evolution must be retained, applies fast light elsewhere, and loads only snapshots spanning the relevant delays. Tests with M87*-like magnetically arrested disk GRMHD data show that the contribution-based region and delay-based snapshot restrictions each keep normalized full-image Stokes differences below $4\times10^{-3}$ relative to the corresponding complete calculation. At this accuracy, CoportSL requires 75.3% and 44.7% fewer snapshot layers for near-horizon and jet images, respectively; its per-frame slow-light transfer time remains comparable to fast light. For the two configurations, source-code estimates place the capacities of the principal data structures at 255-657 GiB for a fixed public ipole version and 20.2-37.3GiB for CoportSL, bringing both configurations within workstation-scale memory. Fast–slow comparisons further show close agreement in near-horizon variability, whereas jet variability follows similar overall trends but differs in local peaks and amplitudes; in both cases, fast light misses substantial full-Stokes spatial structure. As the next-generation Event Horizon Telescope (ngEHT) advances toward dynamical imaging and spatially resolved polarimetry, CoportSL provides a computationally practical way to model full-Stokes finite-light-travel-time signatures in extended black hole systems.
Read the full abstract

All authors Fan Zhou, Jiewei Huang, Yuehang Li, Minyong Guo, Bin Chen

High Energy Astrophysical Phenomena

Chandra X-ray imaging and IC/CMB model for the inner jet of PKS 0637-752

We present the X-ray surface brightness maps of the jet from the quasar PKS 0637-752, using two epochs of archival Chandra observations (1999 and 2017). We present, for the first time, a model of the faint inner jet extending from 3.4 to 7 arcsec from the core, interpreting its X-ray emission as inverse Compton scattering of cosmic microwave background photons by the synchrotron-emitting relativistic electrons. The inner jet contrasts with the bright outer part of the jet, dominated by several discrete knots for which upper limits to the Fermi gamma-ray flux rule out the inverse Compton mechanism. It is the first detailed inverse Compton model published for this inner jet. We examine three physically motivated scenarios, distinguished by the high-energy cutoff of the electron distribution, that can account for the observed radio, millimeter/submillimeter, and X-ray emission while remaining consistent with optical and gamma-ray upper limits. We also report for the first time X-rays from the radio lobe at the end of the receding jet.
Read the full abstract

All authors Jaya Maithil, Daniel A. Schwartz, Aneta Siemiginowska, Diana M. Worrall, Preeti Kharb