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

Sunday 6 September

High Energy Astrophysical Phenomena

Detectability of bulk viscosity effects on post-merger gravitational wave signals from binary neutron star mergers

Zhenyu Zhu, Michail Chabanov and Richard O'Shaughnessy

We investigate bulk viscosity (BV) effects on post-merger gravitational wave (PMGW) signals from binary neutron star (BNS) mergers and their detectability using realistic equation of state (EoS) posteriors constrained by recent progress in gravitational-wave (GW) and X-ray observations, as well as neutron skin thickness measurements. A linear fitting formula based on recent BNS simulations with BV is used to estimate the peak frequency shift of the PMGWs caused by BV effects. Subsequently, we evaluate this peak frequency shift for our set of EoS posterior samples. We use the Fisher information matrix and the dataset of simulated observable events from the ET and CE detector network in an optimistic scenario to estimate the measurement accuracy of the peak frequency. We find that BV effects on PMGWs may be detectable only for EoS models with large values of the symmetry energy slope ($L_{\rm sym}$), as favored by neutron-skin thickness measurements. Thus, the detection of BV effects on PMGWs itself can provide abundant information about the symmetry energy and its slope. However, BV effects on PMGWs are generally weak and may be observed only in some extreme and optimistic cases. Therefore, observations of PMGWs and their analyses may not be able to provide very accurate measurements of the symmetry energy through BV effects.
Read the full abstract

High Energy Astrophysical Phenomena

Constraining the Equation of State of Neutron Stars with third-generation Gravitational Wave detectors

Zhenyu Zhu and Richard O'Shaughnessy

We investigated the impact of the number of binary neutron star merger events and neutron star (NS) mass distributions on constraining the equation of state (EoS), tidal deformability and radius of NSs, as well as the nuclear parameters, using binary neutron star inspiral gravitational wave signals with third-generation detectors. We generate simulated gravitational wave signals and compute the Fisher information matrix for each event after the number of events and mass distribution models are given. The covariance of EoS parameters is obtained by holding the non-EoS waveform parameters fixed at their injected values and accumulating the Fisher matrix over all events. Finally, the posterior samples of EoS, tidal deformability, radius and nuclear parameters are generated based on this covariance matrix. We find that larger number of events lead to tighter constraints due to more observed events and data accumulation, as expected given the increased number of detections. Meanwhile, we note that the mass distribution plays a more important role in constraining the EoS. We compare a realistic bimodal Gaussian distribution, a uniform distribution and another uniform including sub-solar mass NSs. The results show that the uniform distribution yields tighter constraints than bimodal models because it includes more low-mass and massive NSs. This also implies that events near $1.4\,\Msun$ provide partly redundant information about the EoS. Additionally, including sub-solar mass NSs can further improve the constraints by significantly reducing the EoS uncertainties at sub-saturation densities and inner-crust, highlighting the importance of sub-saturation density EoS.
Read the full abstract

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

IPHAS J190812.63+045728.1: A deeply eclipsing, X-ray bright cataclysmic variable with frequent outbursts

Raimundo Lopes de Oliveira, Albert Bruch, Koji Mukai et al. (+3 more)

We present a multi-wavelength characterization of IPHAS J190812.63+045728.1. Independently discovered by us after being serendipitously observed by XMM-Newton, followed up with optical photometry with SPARC4/OPD and spectroscopy with GMOS/Gemini South, and exploring archival TESS and ZTF observations, we reveal the system to be an eclipsing accreting white dwarf (WD) system with an orbital period of 5.073 h. The deep optical and flat-bottomed X-ray eclipses completely occult the WD and the accretion disk, providing a model-independent, uncontaminated view of the donor star which we classify as an M2.3 dwarf with a mass of $\sim 0.50\,M_\odot$. The system is a modestly luminous (L$_{x;0.3-10 keV} \sim$ 1.7 $\times$ 10$^{32}$ erg s$^{-1}$) hard thermal (consistent with a shock-heated plasma cooling down from $k$T $\sim$ 26 keV) X-ray source. While these X-ray properties are typical of low luminosity intermediate polars (IPs), long-term archival photometry reveals frequent, large-amplitude outbursts ($\Delta m \sim 2$ mag in the ZTF $zr$ band) that are the typical characteristics of dwarf novae, and would be unusual for IPs. The weakness of the HeII $\lambda$4686 emission line does not allow us to decide between the IP and the dwarf nova scenarios. On the other hand, the optical and X-ray light curves do not exhibit periodic variations caused by the WD rotation. If the system is indeed an IP, the accretion geometry is likely unusual while the dwarf nova-like outbursts indicating disk instabilities would make IPHAS J190812.63+045728.1 a benchmark for testing disk instability models in the presence of magnetically truncated accretion disks. Otherwise, the system is a dwarf nova with X-ray properties at the extreme end of its class.
Read the full abstract

All authors Raimundo Lopes de Oliveira, Albert Bruch, Koji Mukai, A. S. de Araujo, Guilherme P. Guimarães, Ted Leandro de Almeida

Saturday 5 September

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

Connecting the little dots in polarized light

Piero Madau, Roberto Maiolino and Francesco D'Eugenio

Recent observations of the local Little Red Dot (LRD) analog SDSS~J1025+1402 have revealed that the optical continuum is polarized at the ~1.5% level, while broad Halpha has lower polarization (~0.7%) and a different polarization angle. We explain these properties in a scenario in which LRDs are dust-obscured Little Blue Dots (LBDs) viewed at high inclination and powered by super-Eddington accretion. In this framework, electron scattering in the geometrically thick inner accretion flow produces intrinsic continuum polarization, while unpolarized emission from the outer thin disk dilutes the signal, particularly at optical wavelengths. A circumnuclear dust screen contributes dichroic polarization to both the continuum and broad-line emission. Broad Halpha lacks the intrinsic disk component and therefore has lower polarization and a different position angle. We show that this model quantitatively reproduces the observed polarization properties, although degeneracies between accretion rate and inclination remain. These results illustrate the potential of spectropolarimetry as a probe of accretion-flow geometry and orientation in the little-dot population, complementary to spectral-energy-distribution fitting and line diagnostics.
Read the full abstract

High Energy Astrophysical Phenomena

The First Multi-Messenger Nova: External Shocks, TeV Photons, and Neutrinos in the Next Eruption of T CrB

E. Aydi, P. Craig, K. V. Sokolovsky et al. (+10 more)

T Coronae Borealis is the most compelling anticipated Galactic recurrent nova and one of the best opportunities to test whether nova eruptions can become genuine multi-messenger transients. Recent observational work on classical novae has shown that the GeV $\gamma$-ray luminosity correlates strongly with the differential velocity between interacting outflows, suggesting that shock power depends strongly on velocity contrast. We apply this empirical framework to T CrB and compare it with the 2021 eruption of RS Ophiuchi, the only nova securely detected at TeV energies. We argue that T CrB may lie at the extreme high-$\Delta v$ end of the nova population and that its strongest shock should arise at the external interaction between the nova ejecta and the slow circumbinary medium. This makes the external shock the most likely site of any TeV emission, while the GeV signal may contain contributions from both internal, if present, and external shocks. Because T CrB is substantially closer than RS Oph, it is an especially favorable target for TeV detection. The neutrino case is promising but more uncertain because it depends sensitively on the density and structure of the circumbinary material. We also revisit the recurrence time of T CrB. Using the observed scatter of better-sampled recurrent novae, we infer a representative population-scale recurrence of $T_{\rm rec}\simeq80\pm12$ yr. Taken together, its proximity and potentially extreme shock velocities make T CrB perhaps the best opportunity for a nova to become the first securely established multi-messenger nova.
Read the full abstract

All authors E. Aydi, P. Craig, K. V. Sokolovsky, L. Chomiuk, C. C. Cheung, L. Izzo, J. D. Linford, B. D. Metzger, S. Mohamed, I. Molina, K. Mukai, K. J. Shen, J. L. Sokoloski

High Energy Astrophysical Phenomena

An HLLD Implementation for General Relativistic Magnetohydrodynamics in AthenaK

Jacob Fields, George N. Wong and James M. Stone

We present an implementation of an HLLD approximate Riemann solver for the AthenaK astrophysics code with support for full general relativistic magnetohydrodynamics via a tetrad frame transformation. Our implementation uses an initial guess for the HLLD iterative solve which eliminates the need for an additional conserved-to-primitive inversion, which greatly accelerates performance without affecting accuracy. Additionally, by coupling the method with a first-order flux correction, we are able to use the method reliably even when the magnetization exceeds $10^4$, which we achieve in a SANE accretion disk. Our SANE disk shows that HLLD leads to a more strongly magnetized funnel and more accurate horizon fluxes when compared with HLLE. We further apply the new HLLD implementation to an equal-mass binary neutron star merger. For our long-lived remnant, HLLD enhances the magnetic shear stresses in the outer layers and leads to weaker differential rotation. However, due to weaker gravitational wave emissions in the post-merger phase, the remnant is consistently less compact while producing more dynamical ejecta and a more massive disk. The cost of this new solver is relatively modest thanks to the improved initial guess: our accretion disk tests are only ${\sim}10-25\%$ slower than HLLE, and for our binary neutron star runs with a microphysical equation of state, we find that HLLD is only ${\sim}3\%$ slower than HLLE across all runs.
Read the full abstract

High Energy Astrophysical Phenomena

The accuracy of merger times in hierarchical black hole triples

Giulia Fumagalli and Alejandro Vigna-Gómez

Orbit-averaged equations are widely used to model the gravitational-wave-driven evolution of compact binaries. Their validity relies on the radiation-reaction timescale being much longer than the orbital period, an assumption that can break down for highly eccentric systems. In this work, we investigate the impact of orbit averaging on merger-time estimates by comparing orbit-averaged and non-orbit-averaged prescriptions. We apply our analysis to two different astrophysically motivated populations of black hole binaries in hierarchical triple systems, where von Zeipel-Lidov-Kozai oscillations can drive the inner binary eccentricities close to unity. We find that merger-time estimates obtained with orbit-averaged equations remain remarkably accurate across a wide range of eccentricities, even when the orbital timescale exceeds the radiation-reaction timescale. Significant discrepancies arise only for the most extreme and compact binaries, simultaneously characterized by very high eccentricities and small semi-latus recta. In this regime, the merger time becomes strongly dependent on the initial orbital phase, with different phases yielding predictions that can differ by several orders of magnitude. We also compare our results with an orbit-averaged merger-time prescription specifically developed for hierarchical triple systems and find substantial differences between this prescription and both the generic orbit-averaged and non-orbit-averaged calculations in the high-eccentricity regime. While these results demonstrate the robustness of orbit-averaged merger-time estimates for most systems, they also highlight the potential importance of non-orbit-averaged effects when merger times are used to infer compact-binary delay-time distributions and formation channels.
Read the full abstract

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

Physics-based phenomenological modeling of binary black hole hierarchical formation 1: Synthetic universes from globular cluster simulations for GWTC

R. O'Shaughnessy, R. Mechum, M. Qazalbash et al. (+2 more)

We iteratively model the GWTC-5.0 binary-black-hole census with three components. First, a physically normalized population from Rapster globular-cluster simulations spans cluster mass, metallicity, formation redshift, compactness, and natal black-hole spin. Hierarchical mergers reproduce higher-mass events and their larger effective-spin dispersion, provided black holes are born with zero natal spin. Second, a phenomenological field (isolated-binary) channel supplies the low-mass, preferentially aligned population. Its inclusion flattens the cluster compactness likelihood, permitting ordinary dense globular-cluster birth radii without requiring nuclear-cluster-like conditions. Finally, residual tension at $15$–$30 M_\odot$ motivates an intermediate-mass isotropic component representing field remnants reprocessed in clusters. This is a hypothesis, not a detection: current numerical support precludes a reliable evidence comparison. Physical normalization converts the high-mass rate into $\hat{f}_{\rm GC}\simeq0.39\%$, with local rates $R_{\rm cl}\simeq9.1$ and $R_{\rm field}\simeq16.0$ Gpc$^{-3}$ yr$^{-1}$. The model predicts linked mass-spectrum breaks near $35$ and $70 M_\odot$, a $q\simeq0.5$ feature from first-plus-second-generation pairings, and a symmetric effective-spin distribution that broadens sharply above $45 M_\odot$. These correlated, mass-resolved predictions can be tested as the gravitational-wave census grows.
Read the full abstract

All authors R. O'Shaughnessy, R. Mechum, M. Qazalbash, Z. Rosenberg, M. Zeeshan

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

Magnetically Confined Mountains on Accreting White Dwarfs

Pedro H. B. Rossetto, Manoel F. Sousa and Diego A. Falceta-Gonçalves

The hydromagnetic structure of magnetically confined mountains on accreting white dwarfs is computed, alongside the effects of accretion on the reduction of the star's magnetic field. The equilibrium structure of the mountain is obtained by numerically solving the Grad-Shafranov equation with a self-consistent scheme that enforces magnetic flux freezing. For characteristic system parameters, it is shown that magnetic field lines are significantly deformed and dragged equatorwards by mountains with masses $\sim10^{-3}M_\odot$ that have a characteristic height of $\sim10^5\,\mathrm{cm}$. A mathematical relation is obtained for the maximum mass that a white dwarf can confine magnetically given the star's magnetic field and temperature. It is found that the accretion buries the star's magnetic field by reducing it by up to $10\%$ of its pre-accretion value. The star's magnetic field remains dominantly dipolar, but with added multipolar components.
Read the full abstract

Friday 4 September

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

Hadron resonance gas with density-dependent interactions for neutron stars and heavy-ion collisions

Volodymyr Vovchenko, Volodymyr Kuznietsov and Tripp Moss

We present a density-dependent generalization of the van der Waals hadron resonance gas model (DD-HRG) for describing both the hot hadronic matter created in heavy-ion collisions and the cold, dense matter inside neutron stars. Non-resonant interactions are incorporated through a generalized excluded-volume prescription with a density-dependent available-volume fraction, supplemented by an arbitrary density-dependent mean field. With isospin-dependent interaction parameters constrained by empirical properties of nuclear matter, the resulting equation of state extends the causality range to include neutron-star interiors and supports two-solar-mass stars. It also improves the description of lattice QCD thermodynamics and conserved-charge susceptibilities at vanishing baryochemical potential, with lattice data favoring reduced repulsion among strange baryons. This DD-HRG framework is available within the latest version of the open-source Thermal-FIST package.
Read the full abstract

High Energy Astrophysical Phenomena

A Possible X-ray and Gamma-ray Quasi-Periodic Oscillation in GRB 241030A

Noel Klingler, Cecilia Chirenti, M. Coleman Miller et al. (+2 more)

Quasiperiodic oscillations (QPOs) in gamma-ray intensity have been reported from a few gamma-ray bursts. These QPOs could be related to fundamental frequencies in the sources, which would lead to new insights about these systems. Here we report on our analysis of GRB~241030A, a several-minute-long burst with $\sim 10$ gamma-ray flares that have approximate periodicities at f ~ 0.04 Hz and f ~ 0.08 Hz. This quasi-periodicity is seen in X-rays with the Swift X-Ray Telescope (XRT) as well as in gamma-rays with the Swift Burst Alert Telescope (BAT) and the Fermi Gamma-ray Burst Monitor (GBM). Compared with a power spectral model that has only red noise, the Bayes factor in favor of a Lorentzian-shaped QPO is more than $10^{200}$, even when the red noise is described by segments with up to four different slopes. However, the interpretation of this power spectral feature as a QPO rather than as a fluctuation in red noise is complicated by the low frequency of the feature as well as the extremely large number of counts, which amplifies differences from smooth power spectral distributions. After discussing the properties of the apparent QPO, we discuss a few processes that could plausibly modulate the flux at this frequency. Among the mechanisms considered, we find the most plausible explanation to be a GRB jet passing through regularly spaced shells of circumburst material, possibly produced by binary interactions near periastron in an eccentric progenitor system.
Read the full abstract

All authors Noel Klingler, Cecilia Chirenti, M. Coleman Miller, Amy Lien, Simone Dichiara

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

Effective Field Theory of Gravity in Relativistic Media

Beka Modrekiladze

We develop an effective field theory of gravity in relativistic media. Integrating out the medium leaves vacuum gravity with updated Feynman rules: a graviton propagator dressed by the stress-energy two-point function of the environment, medium-induced bulk graviton vertices from higher correlators, and generalized worldline couplings encoded in matched Wilson coefficients. The diagram topologies are unchanged from vacuum, so different media (perfect fluids, collisionless matter, coherent scalar fields such as wave dark matter) are not different theories but different correlators inserted into the same diagrams. We derive the in-medium rules for a relativistic fluid and obtain the full 1PN Einstein-Infeld-Hoffmann potential, the 1PN Stokes drag, the gravitational self-energy and a generalized Christodoulou memory whose new tensor structure records the orientation of the medium's anisotropy, turning the permanent strain into an $\textit{astrophysical weathervane}$. The same rules activate phenomena forbidden in vacuum: the sound pole converts the symmetry-protected, non-running black-hole Love number into a resonant, running one, and graviton splitting $h\to hh$ is open into the longitudinal branch of the dressed propagator, and gives the transverse graviton a width, $\Gamma=G_N\omega^3(1-c_s^2)^2/120c_s^3$, a $\textit{gravitational opacity}$ set by the sound speed. The same vertex makes a graviton suffer dynamical friction. We assess observational prospects, from dephasing and tidal resonances within reach of the Einstein Telescope and LISA to proof-of-principle memory and opacity signatures.
Read the full abstract

High Energy Astrophysical Phenomena

Rapid radio variability in the $z\sim7$ blazar VLASS J0410$-$0139: Indication of a hidden population of weak radio jets at Cosmic Dawn

Silvia Belladitta, Eduardo Bañados, Anniek J. Gloudemans et al. (+7 more)

Doppler boosting allows blazars to be detected out to high-$z$, making them promising probes of the intergalactic medium through the 21 cm forest. We report 0.144$-$11 GHz observations of the most distant known blazar, VLASS J041009.05$-$013919.88 at z$\sim$7, obtained with the upgraded Giant Metrewave Radio Telescope (uGMRT), the LOw Frequency ARray (LOFAR) and the Very Large Array (VLA). The first uGMRT epoch (300$-$820 MHz, April 2023) revealed an inverted radio spectrum which, combined with earlier (2021$-$2022) VLA data (1.5$-$11 GHz), unveiled a double-peaked spectrum potentially indicative of multi-epoch jet activity. A second uGMRT campaign (August 2023), simultaneous with new VLA observations (1.5$-$11 GHz), instead revealed a flat low-frequency and a peaked high-frequency spectrum, ruling out this interpretation. While limited by the two uGMRT epochs, variability analysis favors intrinsic jet processes, indicating a highly relativistic, closely aligned jet ($\theta<3$ deg, $\delta>19.3$, $\Gamma>9.7$). The equipartition magnetic field ($> 1$ mG) exceeds the equivalent Cosmic Microwave Background field at $z\sim7$ (0.2 mG), indicating synchrotron losses dominate. The inferred Doppler boosting ($\delta>19.3$) implies that J0410$-$0139 is intrinsically radio-weak. As a blazar, it traces a much larger parent population of radio quasars at $z\sim7$, detectable only with deep ($\leq \mu$Jy) next-generation radio observations. LOFAR 144 MHz observations (April and July 2024) yielded $\sim$2 mJy, well below the $\sim$8 mJy predicted from uGMRT epochs, confirming strong variability at rest-frame $\sim$1 GHz. Low-frequency monitoring will be crucial for identifying high radio intensity states suitable for future 21 cm forest studies.
Read the full abstract

All authors Silvia Belladitta, Eduardo Bañados, Anniek J. Gloudemans, Timothy W. Shimwell, Emmanuel Momjian, Huib Intema, Chiara Mazzucchelli, Christian Fendt, Bhargav Vaidya, Fabian Walter

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

The youngest white dwarf companion to a millisecond pulsar: Insights from NGC 362D

Greta Ettorre, Emanuele Dalessandro, Mario Cadelano et al. (+7 more)

We report on the identification of the optical counterpart to the recently discovered millisecond pulsar (MSP) NGC362D in the Galactic globular cluster NGC362 based on deep, multi-band, and multi-epoch Hubble Space Telescope observations. Our analysis robustly shows that this object is a very low-mass ($\sim0.18 M_{\odot}$) He white dwarf (WD) still in the pre-cooling phase. Interestingly, a detailed comparison with updated binary evolution models indicates that this object completed the mass-transfer phase only recently ($\sim0.6$ Gyr ago), thus strongly suggesting that COM-NGC362D is the youngest WD companion to a MSP identified to date. Remarkably, in this respect, the photometric properties of NGC362D show, for the first time in this class of objects, significant wavelength-dependent variations that are consistent with the presence of residual circumstellar material. This system therefore provides a valuable test case to directly probe the immediate aftermath of the MSP recycling process and to constrain the early evolutionary stages of proto-WD companions and their radio and optical properties. Our results suggest that residual material can produce radio and optical signatures that mimic those of systems with non-degenerate companions, potentially leading to the misclassification of young MSPs.
Read the full abstract

All authors Greta Ettorre, Emanuele Dalessandro, Mario Cadelano, Alessandro Ridolfi, Cristina Pallanca, Paulo C. C. Freire, Vivek Venkatraman Krishnan, Maurizio Salaris, Franca D'Antona, Rouhin Nag

High Energy Astrophysical Phenomena

Electromagnetic alignment and jet precession around supermassive black holes: Quasi-periodic oscillations in tidal disruption events

Pau Amaro-Seoane, Leif Lui, Alejandro Torres-Orjuela et al. (+1 more)

We evaluate quasi-periodic oscillations and jet formation in tidal disruption events using the covariant formulation of electromagnetic angular-momentum transfer. General-relativistic frame-dragging tears apart misaligned transient accretion flows, forming an isolated inner mini-disk. The accumulation of magnetic flux on the event horizon powers a relativistic jet via the Blandford-Znajek mechanism. Because the magnetic field anchors to the precessing mini-disk, the jet axis rotates, generating geometric modulations in the observed X-ray and radio fluxes. To ensure physical consistency with the force-free magnetosphere required to launch a Blandford-Znajek jet, we model the electromagnetic back-reaction using a split-monopole magnetic field topology. By performing a small-spin expansion of the Noether current density over the event horizon, we derive a closed-form analytical reaction torque exerted by the electromagnetic field on the accretion plasma. We evaluate the resulting kinematics to show that the electromagnetic back-reaction induces a retrograde precession of the mini-disk, coupling with the prograde Lense-Thirring precession to dictate the global oscillation frequency. We formulate explicit predictions for observable transient signals and predict a monotonic attenuation of the peak-to-trough flux ratio as the mini-disk aligns, as well as a specific frequency drift signature characterized by an initial lengthening followed by an asymptotic shortening of the time interval between consecutive flares. We establish an analytical mechanism where magnetic flux depletion stalls alignment, predicting a constant residual modulation amplitude at late times. We formulate a methodology to extract the black hole spin and the magnetic flux density directly from the temporal derivatives of this predicted frequency drift, operating independently of spectral continuum fitting.
Read the full abstract

All authors Pau Amaro-Seoane, Leif Lui, Alejandro Torres-Orjuela, Xian Chen

High Energy Astrophysical Phenomena

Exploring the connection between Fast Radio Bursts and binary neutron star mergers

Barbara Patricelli, Maria Grazia Bernardini, Cecilia Sgalletta et al. (+1 more)

Fast Radio Bursts (FRBs) are highly energetic radio sources whose duration is of the order of milliseconds. The physical origin of these sources is still unknown. Many models suggest magnetars as possible progenitors of FRBs, and this is supported by the association between FRBs and the Galactic magnetar SGR 1935+2154; other proposed progenitors include binary neutron star (BNS) mergers, that are themselves linked to magnetar formation. In this work we investigate the possible connection between FRBs and BNS mergers, including magnetars that might be produced in such events, by comparing the detection rates inferred from synthetic BNS and associated FRB populations with the rates observed by CHIME. We produce a synthetic catalog of BNS mergers by combining recent theoretically predicted BNS merger rate as a function of redshift and the neutron star mass distribution inferred from measurements of Galactic BNSs. Using this catalog we predict the number of BNS systems ending as magnetars (stable or supramassive neutron star) or black holes (formed promptly or after the collapse of a hypermassive neutron star) for different equations of state. We then simulate for each BNS (and therefore for each magnetar remnant) an associated FRB and estimate how many of them can be potentially detected by CHIME. We find that the rate of BNS mergers and the rate of magnetars produced after BNS represents a non-negligible fraction of the FRBs detected by CHIME, both repeating and non-repeating. Although additional formation channels need to be considered to account for the entire population of FRBs, the existence of a fraction of FRBs that might genuinely be associated to BNS mergers has profound implications in the context of multi- messenger astronomy, supporting the systematic searches of coincident detections of FRBs and GWs from a BNS merger with current and future facilities.
Read the full abstract

All authors Barbara Patricelli, Maria Grazia Bernardini, Cecilia Sgalletta, Michela Mapelli

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

Can isolated binaries form unequal-mass binary black-hole mergers with a high-spin primary black hole?

Xiao-Tian Xu, Dong Lai and Bin Liu

GWTC-5 has revealed a subpopulation of merging binary black holes (BBHs) with a high-spin primary black hole (BH) and possibly unequal BH masses. GW241110 additionally exhibits a large spin-orbit misalignment, which suggests a hierarchical-merger origin. However, other formation scenarios are possible or even likely, especially for events without constraints on spin-orbit misalignment. As an alternative to hierarchical mergers, we investigate whether binary evolution can produce unequal-mass BBH mergers with a high-spin primary BH. Rather than performing comprehensive population-synthesis calculations, we examine the evolutionary pathways of forming merging BBHs and assess their uncertainties. We identify two possible pathways for producing unequal-mass BBHs with a high-spin primary. In initially wide binaries, mass-ratio reversal can make the tidally spun-up second-born BH both more massive and more rapidly rotating than the first-born BH; alternatively, in an initially close, unequal-mass binary, the primary star may evolve chemically homogeneously, while the secondary star evolves normally, producing a high-spin first-born BH that is more massive than its companion. Generally, large spin-orbit misalignment can be produced by large natal kicks or tertiary-induced nodal precession and/or Zeipel-Lidov-Kozai oscillations. We conclude that hierarchical mergers are not uniquely required to produce unequal-mass BBHs with a high-spin primary BH, although each isolated-binary pathway faces important theoretical constraints. Future detections of more merger events with primary BH spins around 0.7 would discriminate between binary evolution and hierarchical merger origin.
Read the full abstract

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

Probing Neutron Star Equation of State Universality with Gravitational Waves

Praveer Tiwari, Ajit Kumar Mehta and K. G. Arun

Current neutron star equation-of-state inferences assume that a common equation of state describes all neutron stars. However, {poorly understood dense-matter physics may give rise to distinct long-lived stellar subpopulations following different mass–tidal-deformability relations, invalidating this assumption. In this Letter, we demonstrate that next-generation gravitational-wave observatories enable neutron star equation-of-state universality to be transformed from an assumption into a testable hypothesis. By combining hierarchical Bayesian inference with a piecewise-polytropic parametrization of the neutron star equation of state, we establish a framework for testing equation-of-state universality with next-generation gravitational-wave observatories. We find that about 30 (50) loud binary neutron star mergers with network SNRs above 100 are sufficient to provide compelling evidence for departures from equation-of-state universality, even if only 20% (10%) of the population exhibits one of the representative non-universality parametrizations considered here.
Read the full abstract

High Energy Astrophysical Phenomena

Nonstationary Stochastic Timing Signatures in the Prompt Gamma-Ray Light Curve of GRB 170817A

Lin Xie and Dahai Yan

We investigate time-dependent stochastic structure in the weak prompt gamma-ray emission of GRB 170817A using change-point and deep-kernel Gaussian-process (GP) models. The analysis is based on the 10-300 keV Fermi/GBM light curve with 0.10 s time resolution. Two change-point configurations identify similar covariance transitions at 0.269 and 0.237 s after the gravitational-wave merger, with 10-90% transition widths of 0.441 and 0.393 s, respectively. At the representative gate-defined boundary of t_tr = 0.27 s, all four fixed-split assignments yield positive evidence gains over the full-exposure stationary Matern-3/2 reference. The largest gain is Delta ln Z_split = 5.72 +/- 0.16 for the Matern-3/2 to Matern-3/2 assignment, while the gate-matched damped random walk (DRW) to Matern-3/2 assignment gives Delta ln Z_split = 5.50 +/- 0.16. The fixed-split comparison therefore supports segment-specific covariance evolution without requiring a change of covariance family. The deep-kernel models recover localized time-deformation features peaking at 1.65 and 1.75 s after the merger for DRW and Matern-3/2 base kernels, respectively. Their offsets from the adopted gamma-ray onset at 1.74 s are -0.09 and +0.01 s, both within the 0.10 s sampling resolution. The full-band timing features remain stable when the bin width is changed from 0.10 to 0.12 s, whereas energy subdivision produces larger shifts in gate locations than in warp peaks. Residual diagnostics show that the models reproduce the dominant temporal structure, although localized residual dependence remains. We interpret the gate-defined covariance transition and the localized time-deformation feature as model-dependent phenomenological timing diagnostics. Further simulation calibration and count-level modelling are needed to assess the statistical robustness and physical origin of the recovered nonstationary structure.
Read the full abstract

High Energy Astrophysical Phenomena

The October 2022 flare in OJ 287 and the mass of its primary black hole

Mauri J. Valtonen, Staszek Zola, Manpreet Singh et al. (+21 more)

The bright blazar OJ~287 has demonstrated a sequence of flares, which are well explained by a quasi-Keplerian orbit model. The flares are associated with the impact of the secondary on the accretion disk of the primary. The orbit must precess in order to produce the correct sequence of flares, and from the precession rate we calculate the mass of the primary. This precession rate gives the mass of the primary $M_{BH} = (18.35\pm0.05) \times 10^9 M_{\odot}$. Two kinds of flares have been identified: direct flares from the impacts, and tidal flares arising from an increased accretion flow into the jet. The precession rate and the primary black hole mass may be independently determined from both sets of flares; the tidal flare of October 2022 was recommended for an intense campaign for this reason. This paper describes these observations over a wide spectral range. We show that the October 2022 flare fits the expectations for a tidal flare and thus supports the earlier determination of the mass of the binary black hole system in OJ 287. The mass of the primary may also be deduced from secondary indicators such as the correlation with the hydrogen line strength and the black hole mass. These studies require that the mass is above $M_{BH} \sim 10^{10} M_{\odot}$, but do not specify the value more exactly.
Read the full abstract

All authors Mauri J. Valtonen, Staszek Zola, Manpreet Singh, Andrei V. Berdyugin, Kari Nilsson, A. Gopakumar, Alok C. Gupta, Tapio Pursimo, Alexandr E. Volvach, Marek Drozdz, Waldemar Ogloza, Rene Hudec, Martin Jelínek, Jan Štrobl, Michal Zejmo, Stefano Ciprini, Daniel E. Reichart, Vladimir V. Kouprianov, Alberto Sadun, Markus Mugrauer, Katsura Matsumoto, Ryo Imazawa, Makoto Uemura, Lankeswar Dey

High Energy Astrophysical Phenomena

The angular structure of the GW170817 jet from prompt emission alone

R. Moradi and R. Ruffini

We determine the angular structure of the GW170817 jet by the prompt emission alone, without afterglow fitting, circumburst density or microphysical parameters. We assume that GRB 090510 and GW170817 have outflows of the same kind, observed respectively on-axis and at the interferometric viewing angle of $20^\circ$. We support this assumption with independent gravitational-wave data showing compatible binary masses and radiated energies. We show that for an observer whose beaming cone is filled with outflow, the point-source Doppler scalings do not apply: $E_{\rm iso}=4\pi\epsilon(\theta_v)$, with $\epsilon$ the energy radiated per unit solid angle along the line of sight, while the peak energy follows $E_{\rm p,i}\propto\Gamma(\theta_v)$. We obtain $n={\rm d}\ln\epsilon/{\rm d}\ln\Gamma = 3.76\pm0.29$ from the ratio of the two bursts, with no free parameter and no assumed angle. This excludes four structures in common use at $4.7\sigma$ to $13\sigma$; three remain above $4\sigma$ across the full reported range of the peak energy of the GW170817 jet. Two prompt spectra fix no angular scale; supplying it with the core Lorentz factor of GRB 090510 and the viewing angle gives $\epsilon\propto\theta^{-7.4}$ outside a core of $2^\circ$–$5^\circ$, in agreement with the width inferred from $367$ short bursts, the outflow remaining relativistic at $\Gamma=33$ on the line of sight. The exponent exceeds what the Lorentz boost of a uniform comoving flow can produce, so the structure is intrinsic to the outflow and not a consequence of the boost. The same structure fixes the emission radius, which contributes $0.41$~s of the $1.74$~s delay between the gravitational-wave signal and the gamma-rays, the remainder being the launch and breakout of the jet, with no free parameters. We conclude that the faintness lies in the structure of the GW170817 jet, not in the de-beaming of a bright core.
Read the full abstract

High Energy Astrophysical Phenomena

Through a glass, darkly: a combined framework for estimating fast radio burst host galaxy and population properties in an era of uncertain host identification

C.W. James, B.C. Andersen, L. Marnoch et al. (+5 more)

The identification of fast radio burst (FRB) host galaxies, and subsequently their redshifts ($z$), has allowed the FRB dispersion measure (DM) to be used to probe the cosmological distribution of ionised gas, and study the properties of the FRB population itself. However, current methods cannot account for FRBs with uncertain host galaxy associations, leading to underutilisation of data, and potential biases towards nearby, bright hosts. In this work, we develop a methodology which can. We do so by combining three ingredients - the zDM code, for modelling the Macquart relation; PATH, for statistical host galaxy identification; and a set of models able to describe the intrinsic FRB host galaxy distribution - into a single formalism. We prove the fidelity of our formalism by using a synthetic set of FRB observations, with simulated hosts sampled from galaxy catalogues, and show that it reproduces intrinsic host galaxy properties even for FRB localisation uncertainties of 30", where a traditional PATH analysis produces no confidant host associations. When applied to a sample of FRBs localised by the Australian Square Kilometre Array Pathfinder, we confirm previous results showing that FRBs prefer host galaxies fainter than that predicted by star-formation, consistent with an exponential surface density scaling as $0.41^{+0.13}_{-0.09}$ times the half-light radius. We encourage the application of this formalism to data-sets from other FRB-hunting instruments.
Read the full abstract

All authors C.W. James, B.C. Andersen, L. Marnoch, J.L. Hoffmann, N. Loudas, J.X. Prochaska, S.D. Ryder, M. Woodland

High Energy Astrophysical Phenomena

Evidence of Orbital Evolution of Gigahertz-Peaked Spectra in Binary Millisecond Pulsar

Rahul Sharan and Bhaswati Bhattacharyya

We present the first systematic investigation of the continuous orbital evolution of the Gigahertz-Peaked Spectum turnover frequency ($\nu_{peak}$) in the binary millisecond pulsar, PSR J2144$-$5237. Using approximately 50 hours of densely sampled observations spanning the $\sim10$-day orbit with the Parkes Ultra-Wideband Low frequency receiver, we obtained continuous frequency coverage from 704 MHz to 4 GHz, enabling robust measurements of the spectral evolution throughout the binary orbit. We detect a clear and systematic modulation of the turnover frequency with orbital phase for both the $1^{st}$ pulse and $2^{nd}$ pulse components, with both emission components exhibiting remarkably similar evolutionary trends despite originating from distinct emission regions within the pulsar magnetosphere. The correlated evolution of the two pulse components strongly suggests that the observed spectral variability is governed by propagation through the intra-binary environment. The large orbital modulation of turnover frequency can be explained by cyclotron resonance absorption in magnetized plasma associated with the companion, while alternative absorption mechanisms are unable to reproduce the observed behaviour. Modelling the orbital evolution of the turnover frequency within the cyclotron absorption framework, can constrain the electron density distribution in the pulsar wind, estimate the magnetic field strength of the companion, and place limits on the geometry of the absorbing plasma and the binary system. This study establish orbital evolution of the Gigahertz Peaked Spectum as a powerful new diagnostic of the magnetized intra-binary medium and provide a hitherto unexplored probe of pulsar wind-companion interactions in eclipsing millisecond pulsar binaries.
Read the full abstract

High Energy Astrophysical Phenomena

A Differentiable Neural Surrogate for Photon Propagation in Neutrino Telescopes

Felix J. Yu, Berthy T. Feng, Nicholas Kamp et al. (+1 more)

Large-volume neutrino telescopes infer neutrino properties from Cherenkov light, but simulating the transport of billions of photons through highly scattering ice or water is computationally costly. We introduce candela, a differentiable SIREN neural field that learns the photon Green's function of the IceCube Neutrino Observatory, a cubic-kilometer detector embedded in Antarctic glacial ice. Given a point-like energy deposit and sensor, it predicts the expected photon yield and full arrival-time distribution at the sensor. Complete events are simulated by decomposing charged-particle energy deposits into point-like sources and superposing their predicted sensor responses. Trained on Monte-Carlo simulations, candela generates events $50$–$100\times$ faster than existing methods, with cost scaling only weakly with neutrino energy. It keeps median yields within $2\%$ of the MC expectation and timing distributions at the MC statistical floor across six photon-count decades. The model also provides end-to-end gradients with respect to event parameters and opens a path toward optimizing scattering-medium properties, which often dominate systematic uncertainties in neutrino telescopes.
Read the full abstract

All authors Felix J. Yu, Berthy T. Feng, Nicholas Kamp, Carlos A. Argüelles

Thursday 3 September

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

Constraints on Scalar–Tensor–Vector Gravity Theory Parameters Inferred from Quasiperiodic Oscillations

Marco Muccino, Kuantay Boshkayev, Binay Prakash Akhouri et al. (+2 more)

We study circular geodesics of neutral test particles in the static charged solution of Scalar–Tensor–Vector Gravity (STVG), and we use the twin kilohertz quasiperiodic oscillations (QPOs) of twelve accreting compact objects to bound its parameters. Starting from the effective potential we obtain closed forms for the specific energy, the specific angular momentum and the Keplerian angular velocity, and from these the radial and vertical epicyclic frequencies. The horizon, photon sphere, shadow radius and innermost stable circular orbit are given in closed form as well, and each one reduces to the RN and Schwarzschild value in the appropriate limit. Null geodesics are integrated numerically and show how the enhanced coupling widens the capture cross section while leaving the logarithmic divergence of the deflection angle at the photon sphere intact. Within the RP model we then run Metropolis–Hastings MCMC simulations on the QPO pairs of eight neutron stars and four microquasars, comparing the Schwarzschild spacetime, the RN spacetime, STVG with vanishing charge and the full STVG solution. The best fits are obtained mainly for the charged solutions in the neutron star sample, while all four models describe the black hole sample equally well. We show that the orbital dynamics depends on the mass $M$, the coupling $\alpha$ and the charge $Q$ only through the two combinations $\Meff=(1+\alpha)M$ and $\Qeff^{2}=(1+\alpha)(\alpha M^{2}+Q^{2})$, which accounts for the multimodal posteriors found in the neutron star sample and sets a model-independent limit on what QPO timing alone can measure. The astrophysical implications of these bounds, in particular for inferred masses above the Tolman–Oppenheimer–Volkoff limit, are discussed in detail.
Read the full abstract

All authors Marco Muccino, Kuantay Boshkayev, Binay Prakash Akhouri, Izzet Sakall{\i}, Yassine Sekhmani

High Energy Astrophysical Phenomena

Preferential accretion onto eccentric and unequal binary black holes

Stanislav DeLaurentiis, Zoltán Haiman and Magdalena Siwek

Supermassive binary black holes (SMBBHs) are expected to be surrounded by circumbinary disks (CBDs) which affect the binary through gravitational forces and accretion. It has been reported that the binary can experience “preferential accretion” where one black hole (BH) out-accretes the other for hundreds of orbits, but this asymmetry has yet to be fully described or understood. In this work, we utilize a suite of 80 SMBBH hydrodynamical simulations with varying mass ratios ($q_b$) and eccentricity ($e_b$) in order to robustly delineate the behavior of preferential accretion, determine its relationship to the structure of the CBD, and study its observational consequences. We characterize the accretion-rate ratio $\lambda(t) \equiv \dot{M}_2(t)/\dot{M}_1(t)$ and the mass-ratio rate of change $\dot{q}_b \equiv d/dt(M_2/M_1)$ across the suite. We confirm that the secondary tends to out-accrete the primary ($\lambda \geq 1$), and find this preference to be strongest for low-$e_b$, low-$q_b$ binaries and increasingly time-variable toward high $e_b$. We also find that (i) the time-variability of $\lambda$ tracks the precession of the CBD, (ii) there can be sub- and super-Eddington accretion in a single binary, and (iii) the gas-driven approach toward equal mass becomes particularly slow for highly eccentric, high $q_b$ binaries, suggesting that some binaries may not reach $q_b=1$ within the $30\,\mathrm{Myr}$ lifetime of a quasar and therefore allowing LISA to constrain the accretion history of SMBBHs. Our findings also suggest that periodically flickering jets are a potential observable signature of many binaries.
Read the full abstract

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

Waveform models for the gravitational-wave memory effect: III. Phenomenological frequency-domain model for nonspinning binaries

Arwa Elhashash and David A. Nichols

We present a phenomenological frequency-domain model for the gravitational-wave (GW) memory signal from nonspinning binary-black-hole mergers on quasicircular orbits. We develop separate amplitude and phase models for the dominant $(l,m)=(2,0)$ spherical-harmonic mode of the GW memory signal. The amplitude and phase models are built from superpositions of elementary and transcendental functions, which can be evaluated efficiently. Both portions of the model are calibrated against a numerical-relativity surrogate model over mass ratios from one to eight. Their accuracy is assessed by computing their mismatch with the numerical-relativity models using the Advanced LIGO sensitivity curve from the fourth observing run. The mismatches are of the order $10^{-4}\unicode{x2013}10^{-3}$ over the parameter space of total mass covered by LIGO. The resulting frequency-domain model is a more computationally efficient waveform model for the GW memory signal than a related time-domain model earlier produced by the authors. An open-access implementation of both the time-domain and frequency-domain models is provided in the Python package GWMemoryModel, which can be used for relevant analyses of nonspinning binary black holes.
Read the full abstract

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

AGN-DB: A Unified Multi-Wavelength Database of Active Galactic Nuclei

Alessandro Peca, Nico Cappelluti, C. Megan Urry et al. (+29 more)

We present the Active Galactic Nuclei Database (AGN-DB), a comprehensive, multi-wavelength catalog compiled from more than 100 publicly available AGN catalogs and samples released by the end of 2025, spanning radio to $\gamma$-ray wavelengths. The database contains approximately 8.1 million unique sources, approximately 7.8 million of which remain after flagging stellar contaminants, and approximately 6.8 million of these are classified as AGN. Source cross-matching across catalogs is performed using Lyra, a Bayesian likelihood-ratio framework that jointly considers positional uncertainties, source densities, and photometric information to compute posterior match probabilities. The resulting catalog provides astrometric coordinates, redshifts, photometry, and classifications for each unique source. All multi-catalog provenance is preserved. For every property, we store the full array of values and originating catalog identifiers, enabling multi-epoch and multi-survey analyses. In this paper, we describe the AGN-DB pipeline, including the cross-matching methodology, and present the statistical properties of the v1.0 catalog. AGN-DB is designed to enable population studies, spectral energy distribution modeling, AGN classification, and variability analyses at an unprecedented scale. Its pipeline is designed to facilitate the integration of new catalogs, allowing AGN-DB to be updated regularly, with releases planned at least annually.
Read the full abstract

All authors Alessandro Peca, Nico Cappelluti, C. Megan Urry, Zhongtian Hu, Jerry R. Bonnell, Jack McKeown, Giulia Cerini, Xulei Sun, Fabio Pacucci, Tracey Jane Turner, Peter G. Boorman, Aritra Ghosh, Connor Auge, Rebeka L. Böttger, Adi Foord, Massimiliano Galeazzi, Jeyhan S. Kartaltepe, Iver Warburton Kilmarrin, Allison Kirkpatrick, Michael J. Koss, Stephanie LaMassa, Md Mahmudunnobe, Stefano Marchesi, Lea Marcotulli, Isaac Moskowitz, Priyamvada Natarajan, Mitsunori Ogihara, Meredith Powell, David Sanders, Dominic Sicilian, Chuan Tian, Ezequiel Treister

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

Super-Eddington Accretion and Early-Stage Feedback in Ton S180

Pierpaolo Condò, Giacomo Venturi, Eleonora Parlanti et al. (+15 more)

Narrow-line Seyfert 1 (NLSy1) galaxies are key laboratories for studying rapid supermassive black hole (SMBH) growth and active galactic nucleus (AGN) feedback at high accretion rates. We investigate the nearby NLSy1 Ton S180 with VLT-MUSE optical integral field spectroscopy to connect its nuclear accretion properties with the spatially resolved ionized gas and host-galaxy kinematics. We modeled the unresolved nuclear spectrum and applied a custom point spread function subtraction to recover the host-galaxy emission on kiloparsec scales. The nuclear spectrum requires a complex permitted-line decomposition and a blueshifted [O III] outflow component. Single-epoch estimators and the stellar velocity dispersion imply black hole masses in the range $\log(M_{\rm BH}/M_\odot) = 6.5 - 7.7$. Combined with the observed luminosity, this implies a dimensionless mass accretion rate of $\dot{M} / \dot{M}_{\rm Edd} = 4.1 - 980$, confirming the extreme accretion regime. The host galaxy shows a circumnuclear ring, an inner elongated structure consistent with a bar, and rotation-dominated gas and stellar kinematics. Simple inflow models do not significantly better reproduce the observed velocity field. We detect a resolved ionized outflow extending about 2 kpc west of the nucleus, with mildly blueshifted velocities (with a maximum of $v_{\rm max} \sim 340$ km s$^{-1}$). Its mass outflow rate is only $\sim 4.2 \times 10^{-4}\,M_\odot\,\mathrm{yr}^{-1}$, whereas the unresolved nuclear outflow reaches $v_{\rm max} \sim 1140$ km s$^{-1}$ and $\dot{M}_{\rm out} > 1.5\,M_\odot\,\mathrm{yr}^{-1}$. This contrast may reflect either weak ionized coupling from nuclear to galactic scales or different episodes of AGN activity over time. These results show that Ton S180 is undergoing rapid SMBH growth, while the observed ionized outflow remains confined to the inner few kiloparsecs and is weak on host-galaxy scales.
Read the full abstract

All authors Pierpaolo Condò, Giacomo Venturi, Eleonora Parlanti, Marco Berton, Enrico Congiu, Alessia Tortosa, Stefano Carniani, Francesco Tombesi, Tommaso Zana, Claudio Ricci, Ezequiel Treister, Luis C. Ho, Enrico Piconcelli, Miguel Coloma Puga, Emilia Järvelä, Luca Crepaldi, Amelia Vietri, Benedetta Dalla Barba

High Energy Astrophysical Phenomena

A Kinematic Measurement of the Effective Viscosity of the Intracluster Medium: The Cold Front of A3667 Seen by XRISM

Yutaka Fujita and Yuki Omiya

Constraints on the effective viscosity of the intracluster medium (ICM) come almost exclusively from imaging: the power spectrum of X-ray surface brightness fluctuations, and the morphology of Kelvin-Helmholtz rolls at cold fronts. Both indicate a viscosity below the isotropic Spitzer value. XRISM/Resolve opens a complementary, purely kinematic route, in which the viscosity is inferred from the thickness of the viscous shear layer that develops at a sliding contact discontinuity. We apply it to the prototypical cold front of A3667. The interface is the classical Rayleigh problem, in which momentum diffuses to a depth d in a time T, giving a kinematic viscosity \nu ~ d^2/(4T). The XRISM velocity map shows a line-of-sight velocity jump of 535 km/s across the front, completed within the 100 kpc-wide region immediately inside it, the only region with an anomalously broad line. Hence d ~< 100 kpc. Adopting interaction times set by the sizes of the interacting regions divided by this velocity, we obtain dynamic viscosities of 5.2 x 10^3 and 3.3 x 10^3 g/cm/s for the cool and the hot side of the front, or 6.4 and 0.72 times the Spitzer value. Because d is limited by the angular resolution of XRISM, and the viscosity scales as its square, these are upper limits. The kinematic constraint is consistent with, and independent of, the imaging results.
Read the full abstract