Peer Revenue Fully funded

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 1 September

High Energy Astrophysical Phenomena

Connecting radio pulsars, magnetars, and XDINSs in a unified evolutionary framework using simulation-based inference

Understanding the Galactic population of isolated neutron stars within a unified framework provides key insights into their birth properties, evolutionary pathways, and the connections between different neutron star classes. In this work, we aim to reproduce both the observed radio and quiescent X-ray emission from isolated neutron stars, including radio pulsars, magnetars and X-ray dim isolated neutron stars (XDINSs). We develop a comprehensive population synthesis framework that models neutron star birth properties, their dynamical, rotational, and magneto-thermal evolution, as well as radio and X-ray emission, and the selection effects of corresponding surveys. To simulate realistic X-ray spectra, we account for magnetospheric resonant cyclotron scattering and interstellar absorption. Additionally, we model the observational bias introduced by magnetar outbursts, by linking the outburst rate to magnetic stresses in the stellar crust. We then employ a simulation-based inference method, namely truncated sequential neural posterior estimation, to reconstruct the birth properties, such as the initial magnetic field distribution. We find that the Galactic neutron star population can be described by a two-component log-normal distribution of birth magnetic fields with peaks at $5 \times 10^{12}$ G and $10^{14}$ G, respectively. We further find that the Galactic neutron star birth rate is around 3-4 per century. Our results help establish the contribution of neutron stars as central engines of powerful astrophysical transient phenomena, such as fast radio bursts, super-luminous supernovae and gamma ray bursts.
Read the full abstract

All authors Michele Ronchi, Celsa Pardo-Araujo, Vanessa Graber, Nanda Rea, Clara Dehman, Davide De Grandis, Alessio Marino, Francesco Coti Zelati, Joeri van Leeuwen, José A. Pons

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

Constraining AGN Disk Properties with Gravitational Waves from Inspiraling Stellar-Mass Binary Black Holes in Hierarchical Triple Systems

Space-based gravitational-wave detectors can observe stellar-mass binary black holes (BBHs) long before merger, allowing weak environmental perturbations to accumulate. For binaries embedded in active galactic nucleus (AGN) disks, the local gas density characterizes the environment of the supermassive black hole (SMBH) and compact-object migration. We study whether such signals can constrain this density when a stellar-mass BBH orbits a Kerr SMBH. We evolve the outer orbit with relativistic corrections and gaseous dynamical friction (DF), and construct the detector-frame waveform including BBH inspiral, de Sitter precession, DF phase correction, and moving-source effects. Using Fisher-matrix calculations for sampled systems, we estimate statistical uncertainties and systematic errors. Larger gas densities generally improve the statistical precision of several source and outer-orbit parameters, but also increase systematic errors when DF is omitted. For favorable GW190521-like systems observed by LISA for one year, the disk density can be constrained at the level of $\sigma_{\rho}\sim10^{-12}\text{--}10^{-10}\,{\rm g\,cm^{-3}}$. Such constraints would connect BBH merger environments to the gas structure of galactic nuclei and the conditions that support black hole growth. These results indicate that hierarchical BBH inspirals can probe AGN disk environments, provided that gas effects are modeled consistently.
Read the full abstract

All authors Jie Wu, Mengfei Sun, Jin Li, Zhoujian Cao

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

Polarization Images of Neutron Stars Illuminated by a Thin Accretion Disk: A Comparison with Black Holes

We investigate linear polarization images of static, spherically symmetric neutron stars illuminated by a geometrically and optically thin accretion disk. Neutron star equilibrium configurations are constructed with a polytropic equation of state, and the null geodesic equations and the parallel transport equation for the linear polarization vector are solved in the geometric optics approximation. The numerical results show that the total polarized intensity is generally positively correlated with the total intensity and reaches its maximum near the neutron star surface. As the observer inclination increases, the symmetry of the polarization images is progressively broken. In addition, the magnetic field configuration mainly affects the direction of the polarization vectors, while its influence on the overall polarized intensity distribution is comparatively limited. To further characterize the spatial structure of the polarization direction, we introduce the net electric vector position angle $\chi_{\mathrm{net}}$ and the second azimuthal Fourier mode $\angle\beta_2$. A comparison with polarization images of a Schwarzschild black hole reveals clear differences between the two types of compact objects in the locations of strongly polarized regions and the size of the central region without a polarization signal. These results show that linear polarization images provide information beyond total intensity images for distinguishing neutron stars from black holes.
Read the full abstract

All authors Dan-Dan Feng, Guo-Ping Li, Chen-Yu Yang, Hong Lei

High Energy Astrophysical Phenomena

X-ray Flares in Gamma-Ray Bursts at High and Very-High-Energies

A significant fraction of GRB early afterglows exhibit fast and bright X-ray flares, discovered by Swift X-ray Telescope. Their rapid temporal variability and spectral evolution have suggested an internal shock origin, analogous to prompt MeV emission. However, their physical origin and radiation mechanism remain debated. High-energy (> 100 MeV) and very-high-energy (> 30 GeV) gamma-ray observations provide a powerful probe of flare dissipation region, constraining its size, magnetic field strength, and particle acceleration under the synchrotron self-Compton scenario in optically thin relativistic jets. We present a systematic multi-wavelength study of 66 X-ray flares from 47 GRBs observed by X-ray Telescope over 17 years, all within field of view of Fermi Large Area Telescope. We investigate their GeV counterparts and find that only five flares exhibit significant high-energy emission (> 3 sigma). Broadband spectral modeling indicates that this GeV emission is consistent with standard forward shock afterglow. We further investigate correlations between flare spectral properties and energy fluxes at 1 keV, 10 keV, and 1 GeV. Using a synchrotron self-Compton model, we constrain the physical conditions of the emitting region, including magnetic field strength, bulk Lorentz factor, and emission radius. For the most stringent GeV upper limits, we find a magnetic-to-electron luminosity ratio greater than or equal to 1, implying a highly magnetized emitting region. We predict very-high-energy gamma-ray emission from X-ray flares at early (~500s) and late (~5000s) times and assess their detectability with Cherenkov Telescopes. We find that later X-ray flares provide the most promising targets for follow-up observations owing to improved observational accessibility, sensitivity, and reduced response-time constraints of facilities such as Cherenkov Telescope Array Observatory.
Read the full abstract

All authors Pawan Tiwari, Biswajit Banerjee, Gor Oganesyan, Ansh Chopra, Annarita Ierardi, Marica Branchesi, Stefano Ascenzi, Maria Edvige Ravasio

High Energy Astrophysical Phenomena

Search for Neutrinos from Tidal Disruption Events with IceCube

Tidal disruption events (TDEs) are theorized to produce high-energy neutrinos through photohadronic interactions between accelerated protons and multi-wavelength photons in the accretion disk and outflows. Detecting these neutrinos would provide insight into the dynamics of TDEs. Taking advantage of the recent increase in observed TDEs from wide field-of-view telescopes, we conduct a dedicated search for neutrinos coincident in optical/UV and X-ray wavelengths. We searched for neutrino emission from 89 TDEs selected based on X-ray and optical/UV observations using time-dependent likelihood analysis methods in two parts. First, we searched for emission from individual sources, where we fit the time window of expected neutrino emission. Second, we performed a study of jetted and non-jetted TDE subpopulations using a stacking search with a fixed one year time window. No significant neutrino excess was observed in either search. We set upper limits to the contribution of jetted and non-jetted TDEs detected in optical/UV and X-ray wavelengths to the diffuse astrophysical neutrino flux assuming TDEs are standard candles.
Read the full abstract

All authors R. Abbasi, M. Ackermann, J. Adams, J. A. Aguilar, M. Ahlers, J.M. Alameddine, S. Ali, N. M. Amin, K. Andeen, C. Argüelles, S. Athanasiadou, S. N. Axani, R. Babu, X. Bai, A. Balagopal V., S. W. Barwick, V. Basu, R. Bay, J. J. Beatty, J. Becker Tjus, P. Behrens, J. Beise, C. Bellenghi, S. Benkel, S. BenZvi, D. Berley, E. Bernardini, D. Z. Besson, E. Blaufuss, L. Bloom, S. Blot, F. Bontempo, J. Y. Book Motzkin, C. Boscolo Meneguolo, S. Böser, O. Botner, J. Böttcher, J. Braun, B. Brinson, Z. Brisson-Tsavoussis, L. Brusa, R. T. Burley, D. Butterfield, K. Carloni, J. Carpio, N. Chau, Y. C. Chen, Z. Chen, D. Chirkin, S. Choi, A. Chubarov, B. A. Clark, D. A. Coloma Borja, A. Connolly, J. M. Conrad, D. F. Cowen, C. De Clercq, J. J. DeLaunay, D. Delgado, T. Delmeulle, S. Deng, P. Desiati, K. D. de Vries, G. de Wasseige, T. DeYoung, J. C. Díaz-Vélez, S. DiKerby, T. Ding, M. Dittmer, A. Domi, L. Draper, L. Dueser, D. Durnford, K. Dutta, M. A. DuVernois, T. Ehrhardt, L. Eidenschink, A. Eimer, C. Eldridge, P. Eller, E. Ellinger, D. Elsässer, R. Engel, H. Erpenbeck, W. Esmail, S. Eulig, J. Evans, P. A. Evenson, K. L. Fan, K. Fang, K. Farrag, A. Fattorini, A. R. Fazely, A. Fedynitch, N. Feigl, C. Finley, D. Fox, A. Franckowiak, S. Fukami, P. Fürst, J. Gallagher, E. Ganster, A. Garcia, M. Garcia, E. Genton, L. Gerhardt, A. Ghadimi, C. Glaser, T. Glüsenkamp, J. G. Gonzalez, S. Goswami, A. Granados, D. Grant, S. J. Gray, S. Griffin, S. Griswold, K. M. Groth, D. Guevel, C. Günther, P. Gutjahr, C. Ha, A. Hallgren, L. Halve, F. Halzen, L. Hamacher, M. Handt, K. Hanson, J. Hardin, A. A. Harnisch, P. Hatch, A. Haungs, J. Häußler, K. Helbing, J. Hellrung, B. Henke, L. Hennig, F. Henningsen, L. Heuermann, R. Hewett, N. Heyer, S. Hickford, A. Hidvegi, C. Hill, G. C. Hill, R. Hmaid, K. D. Hoffman, A. Hollnagel, D. Hooper, S. Hori, K. Hoshina, M. Hostert, W. Hou, M. Hrywniak, T. Huber, K. Hultqvist, K. Hymon, A. Ishihara, W. Iwakiri, M. Jacquart, S. Jain, O. Janik, M. Jansson, M. Jin, N. Kamp, D. Kang, W. Kang, A. Kappes, L. Kardum, T. Karg, A. Karle, A. Katil, M. Kauer, J. L. Kelley, M. Khanal, A. Khatee Zathul, A. Kheirandish, T. Kim, H. Kimku, F. Kirchner, J. Kiryluk, C. Klein, S. R. Klein, Y. Kobayashi, S. Koch, A. Kochocki, R. Koirala, H. Kolanoski, T. Kontrimas, L. Köpke, C. Kopper, D. J. Koskinen, P. Koundal, M. Kowalski, T. Kozynets, A. Kravka, N. Krieger, T. Krishnan, K. Kruiswijk, E. Krupczak, E. Kun, N. Kurahashi, C. Lagunas Gualda, L. Lallement Arnaud, M. J. Larson, F. Lauber, J. P. Lazar, K. Leonard DeHolton, A. Leszczyńska, C. Li, J. Liao, C. Lin, Q. R. Liu, Y. T. Liu, M. Liubarska, C. Love, L. Lu, F. Lucarelli, W. Luszczak, Y. Lyu, M. Macdonald, E. Magnus, Y. Makino, E. Manao, S. Mancina, A. Mand, I. C. Mariş, S. Marka, Z. Marka, L. Marten, I. Martinez-Soler, R. Maruyama, J. Mauro, F. Mayhew, F. McNally, K. Meagher, A. Medina, M. Meier, Y. Merckx, L. Merten, J. Mitchell, L. Molchany, S. Mondal, T. Montaruli, R. W. Moore, Y. Morii, A. Mosbrugger, D. Mousadi, E. Moyaux, T. Mukherjee, M. Nakos, U. Naumann, R. Neshat, L. Neste, M. Neumann, M. U. Nisa, K. Noda, A. Noell, A. Novikov, A. Obertacke, V. O'Dell, A. Olivas, R. Orsoe, J. Osborn, E. O'Sullivan, B. Owens, V. Palusova, H. Pandya, A. Parenti, C. Parisel, N. Park, V. Parrish, E. N. Paudel, L. Paul, C. Pérez de los Heros, T. Pernice, T. C. Petersen, J. Peterson, S. Pick, M. Plum, A. Pontén, V. Poojyam, B. Pries, R. Procter-Murphy, G. T. Przybylski, L. Pyras, C. Raab, J. Rack-Helleis, N. Rad, M. Ravn, K. Rawlins, Z. Rechav, A. Rehman, I. Reistroffer, E. Resconi, C. D. Rho, W. Rhode, L. Ricca, B. Riedel, A. Rifaie, E. J. Roberts, S. Rodan, M. Rongen, A. Rosted, C. Rott, T. Ruhe, L. Ruohan, D. Ryckbosch, J. Saffer, D. Salazar-Gallegos, P. Sampathkumar, A. Sandrock, G. Sanger-Johnson, M. Santander, S. Sarkar, M. Scarnera, M. Schaufel, H. Schieler, S. Schindler, L. Schlickmann, B. Schlüter, F. Schlüter, N. Schmeisser, T. Schmidt, F. Schmitt, A. Scholz, F. G. Schröder, S. Schwirn, S. Sclafani, D. Seckel, L. Seen, M. Seikh, S. Seunarine, P. A. Sevle Myhr, R. Shah, S. Shah, A. Sharma, S. Shefali, N. Shimizu, M. Shin, B. Skrzypek, R. Snihur, J. Soedingrekso, D. Soldin, P. Soldin, G. Sommani, D. Song, C. Spannfellner, G. M. Spiczak, C. Spiering, J. Stachurska, M. Stamatikos, T. Stanev, T. Stezelberger, T. Stürwald, T. Stuttard, G. W. Sullivan, I. Taboada, S. Ter-Antonyan, A. Terliuk, A. Thakuri, M. Thiesmeyer, W. G. Thompson, J. Thwaites, W. Tian, S. Tilav, K. Tollefson, J. A. Torres, S. Toscano, D. Tosi, K. Upshaw, A. Vaidyanathan, N. Valtonen-Mattila, J. Valverde, J. Vandenbroucke, T. Van Eeden, N. van Eijndhoven, L. Van Rootselaar, J. van Santen, J. Vara, F. Varsi, M. Velazquez, M. Venugopal, M. Vereecken, S. Vergara Carrasco, S. Verpoest, A. Vijai, J. Villarreal, C. Walck, A. Wang, E. H. S. Warrick, C. Weaver, P. Weigel, A. Weindl, J. Weldert, A. Y. Wen, C. Wendt, J. Werthebach, M. Weyrauch, N. Whitehorn, C. H. Wiebusch, D. R. Williams, L. Witthaus, J. Woodward, G. Wrede, X. W. Xu, J. P. Yanez, Y. Yao, E. Yildizci, S. Yoshida, F. Yu, S. Yu, T. Yuan, S. Yun-Cárcamo, A. Zander Jurowitzki, A. Zegarelli, S. Zhang, Z. Zhang, P. Zhelnin, P. Zilberman, C. Zilleruelo Cañas

High Energy Astrophysical Phenomena

Nanohertz Gravitational-Wave Constraints on Supermassive Binary Black Holes at Cosmic Dawn

Standard continuous gravitational-wave searches with pulsar timing arrays (PTAs) neglect cosmological redshift, restricting their applicability to the local Universe. We introduce a redshift-aware PTA framework and apply it to the Parkes PTA Data Release 3, deriving the first direct constraints on supermassive binary black holes (SMBBHs) at the cosmic dawn. Evaluating our limits across a broad redshift range, we observationally establish the non-monotonic mass-redshift exclusion boundary driven by the theoretical “redshift bias", demonstrating how PTAs can effectively probe sources at extreme distances. Redshift-aware targeted searches toward high-redshift systems, including the ultraluminous quasar J0100+2802 ($z=6.327$) and the JWST-discovered galaxy JADES-GS-z14-0 ($z=14.32$), rigorously exclude SMBBHs with chirp masses $\mathcal{M}_c \gtrsim 10^{10} M_\odot$ across the nanohertz frequency band. Finally, we demonstrate that high-redshift binaries can be robustly detected and localized, allowing for the accurate measurement of their intrinsic properties. Our results provide the tightest constraints to date on SMBBHs at $z > 6$ and establish a practical framework for probing early-Universe black hole assembly.
Read the full abstract

All authors Yiqin Chen, Shi-Yi Zhao, Xingjiang Zhu, N. D. Ramesh Bhat, Jacob Cardinal Tremblay, Małgorzata Curyło, Valentina Di Marco, George Hobbs, Wenhua Ling, Rami F. Mandow, Richard N. Manchester, Saurav Mishra, Daniel J. Reardon, Sparrow Roch, Christopher J. Russell, Ryan M. Shannon, Jingbo Wang, Shuangqiang Wang, Andrew Zic

Monday 31 August

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

Gravitational waves from core collapse of rotating very-massive stars: 3D numerical relativity computation

We numerically study the collapse of rotating, very-massive stellar cores with masses of $\approx 200$, $300$, $500$, and $1100M_\odot$ into black holes using both axisymmetric and three-dimensional (3D) numerical relativity. Our results indicate that when the dimensionless spin of the resulting black hole exceeds 0.8, a massive disk consistently forms around it. These massive, compact disks, carrying more than about 10% of the black hole's mass, are prone to non-axisymmetric deformations that trigger gravitational-wave bursts with frequencies around 10-50 Hz. Such waves could be detected by the Einstein Telescope and Cosmic Explorer, even from sources a few Gpc away. We also summarize the gravitational-wave signals from axisymmetric collapse, which tend to have lower amplitudes and higher frequencies than those caused by non-axisymmetric instabilities.
Read the full abstract

High Energy Astrophysical Phenomena

GRMHD Simulations of Accreting Proto-Magnetars II. Implications for r-process Nucleosynthesis

Newly formed, rapidly spinning, strongly magnetized neutron stars ("millisecond proto-magnetars") can arise in collapsars, neutron star mergers, or white-dwarf accretion-induced collapse, and are often surrounded by compact accretion disks. At accretion rates of ~0.1 Msun/s, these disks can become neutron rich and power outflows capable of rapid neutron-capture (r-process) nucleosynthesis. In Paper I, we presented axisymmetric GRMHD simulations of accretion onto such proto-magnetars and showed how the disk-magnetosphere interaction regulates jet power, variability, and neutron star torques. Here we use the same simulations to study how this interaction regulates the mass, composition, and velocity of the baryon-rich ejecta, comparing proto-magnetar models to otherwise similar black hole accretion. A magnetized neutron star qualitatively changes both the amount and composition of the ejecta. Stronger neutron star magnetic fields suppress accretion and redirect more inflowing material into unbound outflows, even when the magnetosphere remains strongly compressed by the disk. Once the field produces magnetic channeling or centrifugal acceleration, mass loss is enhanced further. Reaction-network calculations show that these magnetically driven outflows can synthesize the full range of r-process nuclei, including the heaviest elements. Moderate neutrino irradiation substantially reduces the third-peak yield, but magnetically accelerated neutron star outflows retain a heavy component more effectively than black hole disk winds; sufficiently strong early-time irradiation suppresses it altogether. Accreting proto-magnetars may therefore be important heavy r-process sources once their neutrino emission has sufficiently declined.
Read the full abstract

All authors Tejas Prasanna, Kyle Parfrey, Brian D. Metzger, Ore Gottlieb, Andrei Beloborodov, Koushik Chatterjee

High Energy Astrophysical Phenomena

Portraits of the young and old in X-ray: New millisecond pulsar detections and updated characterization of young neutron stars

We present analyses of new and archival Chandra and NuSTAR data on two old millisecond pulsars, a young highly-energetic pulsar, and a young neutron star. PSR J1810-0623 is a recently-discovered 4.55 ms radio pulsar in a binary with a possible carbon-oxygen white dwarf companion. We tentatively detect the pulsar for the first time in X-ray using a new short Chandra HRC-I observation. Its faintness suggests a distance more in accord with the further (~1 kpc) of its dispersion measure-inferred distances. PSR J1933-6211 is a 3.54 ms radio pulsar also with a carbon-oxygen white dwarf companion. Our first-time X-ray detection in an archival HRC-S image indicates its X-ray luminosity is similar to that of the 1.9 Msun pulsar PSR J1614-2230. Next, PSR J1813-1749 is a 1-2 kyr pulsar with one of the highest known spin-down luminosities. We measure its spin period of 44.7 ms using 2025 NuSTAR and 2026 Chandra ACIS-S data and update its spin evolution model. We show for the first time the pulse profile of PSR J1813-1749 at 3-79 keV, and we find that the spectra of the pulsar and pulsar wind nebula at these energies did not change significantly between observations in 2018 and 2025. Finally, for the 5-10 kyr neutron star CXOU J182913.1-125113 in the supernova remnant G18.95-1.1, we combine ACIS-I observations from 2009 and 2024 to better measure the neutron star's spectrum.
Read the full abstract

High Energy Astrophysical Phenomena

Evidence for a luminosity-class dichotomy in the supergiant donors of persistent high-mass X-ray binaries and supergiant fast X-ray transients

We have investigated the properties of the supergiant donors in two types of high-mass X-ray binaries (HMXBs, i.e. the so-called supergiant fast X-ray transients (SFXTs) and the classical X-ray sources), by collecting all available information from the literature. A prevalence of O-type versus B-type supergiant donors in SFXTs was previously suggested in the literature, possibly indicative of faster winds in SFXTs than in classical, persistent sources. Our aim is to investigate whether there is evidence for a dichotomy in the stellar luminosity class, since B-type supergiants are also known to be hosted in SFXTs. Our collection of the relevant literature allowed us to build a robust sample of both types of HMXBs with a secure stellar classification. We have found that the SFXTs showing the most extreme, transient X-ray emission host only O-type or B-type stars with a Ib and Iab luminosity class, whereas classical sources have mostly Ia donors. This is robust evidence that might explain the two different behaviours of HMXBs at X-rays since the stellar luminosity class appears to be linked to the wind properties (the supergiant wind terminal velocity, the mass-loss rate, and the wind strength Q parameter), indicative of faster and less dense winds in SFXTs than in classical sources.
Read the full abstract

2609.00175 · 31 Aug · Accepted for publication in Astronomy & Astrophysics (received May 18, 2026; accepted August 11, 2026). The manuscript comprises 11 pages, 5 figures, and 4 tables. The version submitted here also includes all supplementary tables (B.1, C.1, C.2, and D.1), which in the published A&A article are provided exclusively in electronic form via the CDS

High Energy Astrophysical Phenomena

GRMHD Simulations of Accreting Proto-Magnetars I. Implications for Gamma-Ray Burst Jets and Energetic Explosions

Newly formed, rapidly rotating, strongly magnetized neutron stars ("millisecond proto-magnetars") are promising central engines for gamma-ray bursts (GRBs) and luminous supernovae. Although often modeled in isolation, they can be born surrounded by accretion disks in stellar collapse, neutron-star mergers, or accretion-induced collapse. We present axisymmetric GRMHD simulations of hyperaccretion onto such objects, including a physical equation of state and charged-current weak interactions. Holding the weakly magnetized accretion torus fixed, we vary the stellar dipole field strength to span crushed-magnetosphere, magnetically channeled accretion, and centrifugal-propeller regimes, and compare with an otherwise similar accreting black hole. Accretion compresses the stellar magnetosphere and opens additional magnetic flux, producing relativistic jet powers that exceed isolated-dipole spin-down estimates by factors of a few to ~10. Even while the magnetosphere remains compressed against the stellar surface, stronger fields increasingly impede accretion and enhance outflows. Channeled-accretion models show strong jet variability driven by plasmoid eruptions and intermittent magnetospheric accretion, whereas the propeller model produces a steadier, more powerful jet and rapid spin-down. The disk-magnetosphere interaction also regulates how efficiently the neutron star grows and whether it spins up or down; near spin equilibrium, inefficient accretion can delay collapse to a black hole relative to estimates based on the external mass-supply rate. Accreting proto-magnetars can therefore power relativistic jets and baryon-rich outflows with energetics comparable to those inferred for long GRBs and GRB-supernovae. A companion paper explores implications for neutron-rich ejecta and r-process nucleosynthesis.
Read the full abstract

All authors Tejas Prasanna, Kyle Parfrey, Brian D. Metzger, Ore Gottlieb, Andrei Beloborodov, Koushik Chatterjee

High Energy Astrophysical Phenomena

Equation of state and neutron star properties with new mass-radius constraints from PSR~J1614–2230, PSR~J2124–3358, and 47~Tuc~X7

We study the impact of new mass-radius information from PSR J1614-2230, PSR J2124-3358, and 47 Tuc X7 in a combined equation of state inference based on chiral effective field theory constraints at nuclear densities and using different high-density extensions, including perturbative QCD constraints. The largest impact stems from the heavy-mass PSR J1614-2230 star, which shifts heavy neutron stars to smaller radii by around 0.4 km. Moreover, the combined astrophysical NICER, LIGO/Virgo, and X-ray information drives the radius posterior to a more data-driven distribution, which is less sensitive to the high-density extension. For the equation of state, the new mass-radius information significantly tightens the pressure and speed-of-sound posterior distributions, especially around three times saturation density. Finally, we make predictions for the poorly constrained masses of PSR J2124-3358 and 47 Tuc X7 based on the combined equation of state analysis and the other astrophysical sources.
Read the full abstract

All authors Melissa Mendes, Isak Svensson, Hannah Gottling, Kai Hebeler, Achim Schwenk, Nathan Rutherford, Lucien Mauviard, Christine Kazantsev, Yves Kini, Denis Gonzalez-Caniulef, Sebastien Guillot, Anna Watts

High Energy Astrophysical Phenomena

A NICER view of PSR J1614$-$2230: a massive and compact millisecond pulsar

Using pulse profile modeling, we obtain the mass-radius measurement of a millisecond pulsar (MSP) with data from the Neutron Star Interior Composition ExploreR, XMM-Newton and the Chandra X-ray Observatory. We report here the radius of PSR J1614$-$2230, the second most massive MSP confirmed by radio timing. All of the data sets are well described by a simple model composed of two circular hot spots. The final result yields an equatorial radius of $R_{\rm eq}=10.06 ^{+1.25}_{-0.87}\,$km, and a gravitational mass of $M=1.937^{+0.012}_{-0.013}\,M_{\odot}$ (equally tailed 68% credible intervals). Although a non-thermal component was previously reported at higher energies, we find no sign of it in either our phase-averaged or phase-resolved spectral analyses. Using new relations linking the compactness to oblateness or surface gravity, and tailored to the spin frequency of PSR J1614$-$2230, we infer a configuration with one hot spot near the pole, and another near the equator. The tight mass posterior is essentially informed by radio timing, while the radius constraint is not as tight due to the low source signal (8.5$\sigma$ X-ray pulse significance). However, over all geometries and atmosphere models tested, the radius posterior tends toward low values ($\lesssim12.08\,$km, 90th percentile in all cases).
Read the full abstract

All authors Lucien Mauviard, Sebastien Guillot, Lami Suleiman, Yves Kini, Denis González-Caniulef, Christine Kazantsev, Pierre Stammler, Devarshi Choudhury, Bas Dorsman, Mariska Hoogkamer, Daniela Huppenkothen, Tuomo Salmi, Anna L. Watts, Jérôme Novak, Natalie A. Webb, Jean-Francois Olive, Matthew Kerr, Lucas Guillemot, Ismaël Cognard, Gilles Theureau

High Energy Astrophysical Phenomena

Redshift Identifiability from Gamma-Ray Burst Prompt Emission

We examine which prompt-emission observables of gamma-ray bursts (GRBs) carry identifiable redshift information in realistic structured-jet scenarios. By examining the photon-level Lorentz and cosmological transformations, we demonstrate any spectral observable reflecting the dominant-patch emission and any temporal observable reflecting the co-moving frame physics in the jet suffers from the same $z$-$\theta_v$ degeneracy. The observable reflecting the absolute flux or count rate will also suffer from this $z$-$\theta_v$ degeneracy in a different form. Due to the large variation of the intrinsic luminosity, combining the observables in the above-mentioned three classes cannot break this degeneracy. However, we demonstrate that two other classes of observables contain the potential to break this degeneracy: (i)~central engine-imprinted time scales, (ii)~emission from out of the dominant patch.
Read the full abstract

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

Magnetic field-induced enhancement and quenching of Urca emission in quark matter

Using first-principles field-theoretic methods, we investigate neutrino emission from strongly magnetized dense quark matter under conditions relevant to compact stars. We account for Landau-level quantization of both electron and quark states and show that it strongly modifies the kinematics of Urca processes. In particular, quark quantization restricts the available phase space in which both quark and electron energies can simultaneously lie near their respective Fermi surfaces. At moderately strong magnetic fields, before pronounced quark quantization sets in, the emission rate tends to increase on average with increasing field strength. In the regime of very strong fields, however, the increasingly restricted phase space first gives rise to Shubnikov–de Haas-type oscillations and then to resonance-like spikes near a discrete sequence of Urca-resonant magnetic field values, separated by regions of strong suppression. Finally, the emission rate becomes nearly completely quenched once $|eB| \gtrsim 6\mu_{e}\mu_{u}$, corresponding to approximately $B\gtrsim 1.5\times 10^{19}~\mbox{G}$ for the representative set of model parameters considered. We also find significant anisotropy in the longitudinal momentum emission near the Urca-resonant magnetic field values.
Read the full abstract

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

Extended radio emission in and around the merging cluster A520 as seen by MeerKAT

We present results from a deep MeerKAT 1.28 GHz observation of Abell 520, a galaxy cluster in the middle of a merger and a treasure trove of merger-generated phenomena. MeerKAT's combination of high angular resolution and sensitivity to extended radio emission lets us uncover unprecedented detail of the cluster's giant halo and puzzling features beyond the halo. We find that the main radio halo is sharply bounded along the direction of the NE-SW merger axis by the two X-ray shock fronts. In the perpendicular direction, the halo wings smoothly span over 2 Mpc, with the spectrum steepening in the outskirts. We detect faint radio emission upstream of the main (SW) bow shock; the emissivity jump across this shock can be used to constrain the physical mechanisms at work there. We discover a tail of radio emission extending past the NE shock to the gas-poor, high-entropy “dark subcluster” in the NE outskirts. The tail and subcluster have a higher radio/X-ray brightness ratio than the main halo's, suggesting efficient reacceleration of cosmic rays. We also discover two peripheral radio relics beyond the cluster virial radius, and a faint, $\sim2$ Mpc long, broad radio bridge connecting the dark subcluster to one of the relics along the direction perpendicular to the merger axis. The spectral slopes of the dark subcluster and the bridge are consistent with that of the main halo. This bridge and both relics may be caused by another giant merger-generated shock front at the virial radius.
Read the full abstract

All authors S. Giacintucci, W. D. Cotton, M. Markevitch, T. Venturi, T. E. Clarke, H. Bourdin

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

Exploring stellar evolution with Luminous Red Novae and Interacting Gap Transients

Interacting gap transients encompass several families of astrophysical phenomena where eruptive mass loss and interaction play an important role in shaping the observational properties and energetics of the events. Composed of Luminous Red Novae (LRNe), Intermediate Luminosity Red Transients (ILRTs), Luminous Blue Variable (LBV) giant eruptions, and pre-supernova precursors, interacting gap transients offer crucial insights into unstable mass transfer, extreme stellar mass loss, and binary evolution. In these proceedings, I review the observational signatures, progenitors, physical origin, and classification challenges of these heterogeneous phenomena. Finally, I highlight how ongoing and future wide-field time-domain missions will revolutionize this field by compiling complete volume-limited samples, tracking faint pre-outburst precursors, and establishing statistical links between gap transient families and their progenitor populations.
Read the full abstract

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

Dynamics and Frequency Conversion of Accreting Axion Clouds

Axion fields can form exponentially growing gravitational clouds around compact objects through self-interaction-driven relaxation of ambient axion waves. As the field amplitude approaches the axion decay constant, nonlinear effects become important. We identify two distinct regimes of late-time evolution, determined by the gravitational fine-structure constant and the cloud growth rate: a Bosenova regime, characterized by collapse accompanied by explosive axion bursts, and a saturation regime, in which self-interaction-induced axion emission balances accretion. In the latter regime, the emitted axion radiation exhibits stable discrete spectral lines at odd multiples of the bound-state energy, directly probing the global structure of the axion potential beyond its quadratic minimum. We show that single-cosine potentials and QCD axion-like potentials predict distinct emission spectra, enabling probes of the underlying axion self-interaction structure and its ultraviolet completion through terrestrial detection of relativistic axion fluxes from compact objects.
Read the full abstract

All authors Ximeng Li, Zhen Zhong, Yifan Chen, Yue Zhao, Hidetoshi Omiya, Vitor Cardoso

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

Solar Spin-Dependent Dark Matter-Neutron Cross Section Constraints: The Lost Case

The Sun has long served as a natural dark matter detector, capturing halo particles that scatter with solar nuclei and electrons and subsequently produce indirect signals from dark matter annihilation that can be probed by either neutrino or $\gamma$-ray instruments. Previous efforts have computed cross-section constraints for dark matter scattering with electrons, as well as spin-independent and spin-dependent dark matter-proton scattering. However, the spin-dependent dark matter-neutron scattering scenario has been neglected, despite its importance in direct detection. For the first time, we compute the capture and evaporation rates for spin-dependent dark matter-neutron scattering in the Sun, including odd-neutron isotope targets in the current standard solar model. Comparing the predicted annihilation signals against current neutrino and $\gamma$-ray observations, and projecting the reach of upcoming detectors, we find that solar constraints exceed direct detection limits for several annihilation channels. For models in which the dark matter annihilates into long-lived mediators, these constraints can extend below the neutrino fog.
Read the full abstract

High Energy Astrophysical Phenomena

Causal self-consistency of the Blandford–McKee self-similar solution

The Blandford–McKee (BM) solution describes the ultra-relativistic self-similar flow behind a strong spherical blast wave propagating into an external density profile $\propto r^{-k}$, where $k<4$ and $r$ is the distance from the center, but it applies only to the hot shell adjacent to the shock, sufficiently deep behind it, the fluid leaves the BM regime. Since the similarity profiles are determined solely by the shock conditions, with no conditions imposed at the inner end, the validity of the solution near the shock depends on whether the flow beyond the BM regime can imprint on this hot shell. For shallow density profiles, $k<k_g\simeq 2.062$, the characteristic structure of the BM solution already prevents forward-going acoustic information from the non-BM interior from reaching the shock during the ultra-relativistic stage. For steeper profiles, $k_g<k<4$, the hot equation of state fails while the flow is still relativistic and shock-connected. There, we derive a new self-similar solution for the cooling relativistic flow, which has its own similarity scale, equations, and characteristic structure. It overlaps with the hot BM solution toward the shock, and reaches a sonic point beyond this overlap. This critical point separates the shock-connected BM-plus-cooling composite from the deeper downstream region, so acoustic signals generated beyond it cannot propagate toward the shock. The mechanism is explicit at $k=7/2$, for which the cooling solution is obtained analytically. This causal structure is therefore what enables the BM solution, which remains self-consistent near the shock even though it does not globally describe the full downstream flow.
Read the full abstract

All authors Gilad Sadeh, Tamar Faran, Doron Kushnir, Eli Waxman

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

Weighing Little Red Dots with Transient Events

Recent JWST observations have revealed a large population of compact red sources at $z \gtrsim 4$, known as Little Red Dots (LRDs), many of which show signatures of accreting massive black holes (BHs). The physical nature of these sources and their connection to host galaxies are under debate. We propose an independent avenue for constraining their nature through transient phenomena, such as tidal disruption events (TDEs) and quasi-periodic eruptions (QPEs), arising from interactions between a star and the gas envelope surrounding the BH. These event rates depend sensitively on BH mass and provide a way to "weigh" LRDs. We calculate the expected TDE and QPE rates in LRDs under three distinct scenarios: (1) LRDs are truly overmassive BHs, (2) LRDs have BH masses following the classical local scaling relations (and the reported BH masses in observations are overestimated), and (3) the currently observed LRDs are only the tip of the iceberg of a larger population of low-mass BHs. We find that the predicted TDE and QPE rates differ dramatically across scenarios, especially in the presence of steep stellar cusps. The expected TDE rates per degree-square, assuming a Hernquist stellar distribution with a Bahcall-Wolf cusp embedded, are $2.78 \times 10^{-3}$, $1.96 \times 10^{-3}$, and $3.37 \times 10^{-2} \, {\rm yr}^{-1} \, {\rm deg}^{-2}$ for the three scenarios, respectively, while the QPE rates are $1.64 \times 10^{-2}$, $4.72 \times 10^{-2}$, and $4.96 \times 10^{-1} \, {\rm yr}^{-1} \, {\rm deg}^{-2}$. Upcoming wide-field surveys with Euclid, Roman, and LSST may be capable of detecting these high-redshift transient events and obtaining light curves, which encode additional information about the BH mass and the gas structure of LRDs. Stellar transient events will provide valuable insight into the early assembly of massive BHs.
Read the full abstract

All authors Vinh Tran, Xuejian Shen, Oliver Zier, Anna de Graaff, Rohan P. Naidu, Mark Vogelsberger

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

Ad Astra White Paper: A Pitch for the Next 25 Years of NASA's Physics of the Cosmos Program

Astrophysical observations of our universe have been key to our understanding of how the universe works. Shortly after the turn of the millennium, the National Research Council delivered Connecting Quarks with the Cosmos: Eleven Science Questions for the New Century. In the subsequent quarter-century, we have made substantial progress in answering each question. These advancements have, in part, arisen because of the success of major US facilities across several domains of physics, guided by long-term planning documents which still largely focus on these questions. This report seeks to provide a status update on each question, and to outline what space-based facilities are crucial for future progress, intended to guide NASA's preparatory work for the Astro2030 Decadal.
Read the full abstract

All authors Eric Burns, Ivan Agullo, Igor Andreoni, Catherine M. Deibel, Christopher L. Fryer, Natasha Latouf, M. Coleman Miller, Jillian C. Rastinejad, Breann N. Sitarski, Zorawar Wadiasingh

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

The Analysis, not the Aperture: End-to-End Transformer Reconstruction for Imaging Atmospheric Cherenkov Telescopes

Imaging Atmospheric Cherenkov Telescopes (IACTs) detect very-high-energy gamma rays by imaging the nanosecond Cherenkov flash of the air shower they initiate in the Earth's atmosphere. For four decades the first steps of IACT event reconstruction have been essentially unchanged, relying on a heavy parameterisation and dimensionality reduction of the recorded images. This is reasonable when the image is bright, but discards important information when only a few tens of Cherenkov photons are recorded, which is a primary reason why small telescopes perform poorly at sub-TeV energies. We show that this limitation is a property of the analysis rather than of the hardware. We simulate a deliberately simple and idealised compact telescope and treat each event as a short movie that is passed directly to a video vision transformer with a factorised spatio-temporal encoder. A single composite network with a gradient-normalised multi-task loss performs gamma/hadron classification, energy regression and arrival-direction regression at once. This is the first application of a video vision transformer to IACT data. We compare it against an optimised standard analysis on the same dataset. The transformer lowers the energy threshold by a factor of three, from 0.22 to 0.07 TeV, and reconstructs arrival directions down to 0.05 TeV. At 0.2 TeV it increases the effective collection area by a factor of three, and raises the gamma/hadron separation power from an area under the receiver operating characteristic curve of 0.80 to 0.91. At 0.05 TeV, where the standard analysis retains almost nothing, that area grows by nearly two orders of magnitude. These results show promising new opportunities for compact and affordable telescopes operating at sub-TeV energies, paving the way for a broader exploration of time-domain astrophysics.
Read the full abstract

All authors Elli Jobst, Lea Heckmann, Lukas Heinrich, David Paneque

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

Bekenstein–Hod Bound: A $3.3\sigma$ Confirmation from GW250114

Black holes link gravity, thermodynamics and information via limits on the relaxation rate of a perturbed system. The Bekenstein-Hod bound imposes a minimum relaxation time at fixed temperature, but its observational test demands both black hole thermodynamic characterization and decay time measurement. Here we test this bound with GW250114, the loudest gravitational-wave signal yet observed from a binary black-hole merger. We infer the remnant temperature from pre-merger data truncated at least $10\,M$ before the peak and its longest-lived decay time from post-merger data, thereby avoiding direct reuse of the same strain samples. The ringdown frequencies and damping times are allowed to vary independently rather than being fixed to the Kerr spectrum. For the primary $t_{<}=-10\,M$ analysis, the bound is verified at $3.3-3.6\sigma$ across the focal ringdown start times, representing a substantial improvement over the $91\%$ confidence level set by GW150914. The conclusion remains robust under varied pre-merger cutoffs and explicit inclusion of the short-lived first overtone in waveform modelling. This separated-data measurement converts an information-theoretic relaxation bound into a precision test of a single astrophysical black hole.
Read the full abstract

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

The Auger Radio Infill SKALA Extension (ARISE): Science Case and Instrumentation (ARENA 2026)

The Auger Radio Infill SKALA Extension (ARISE) at the Pierre Auger Observatory in Argentina was deployed in 2025 and measures cosmic-ray air showers in the energy region of the Galactic-to-extragalactic transition. ARISE is comprised of 18 SKALA-2 antennas featuring two polarization channels each, deployed within $100\,$m of a surface detector station in the enhancement area of the Pierre Auger Observatory. This area of the surface array features a denser spacing of $433\,$m between surface stations, each equipped with underground muon detectors. One of these surface detector stations provides a trigger for simultaneous readout of all ARISE antenna channels. The wide frequency range of ARISE, from $50$ to $350\,$MHz, includes the sub-band of optimum signal-to-noise ratio for air-shower radio emission against the Galactic radio background. In combination with the dense antenna spacing, this enables a relatively low detection threshold, and ARISE aims at demonstrating full detection efficiency for near-vertical air showers above $100\,$ PeV. As an advantage over the current radio detectors at Auger, which are more efficient for inclined air showers, this would enable low systematic uncertainties for physics analysis combining ARISE radio measurements with coincident measurements of the underground muon detectors in the same area. In this presentation, we will provide an overview over the ARISE instrumentation operating at the Pierre Auger Observatory and will outline the science goals.
Read the full abstract

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

Constraining equation of state of neutron star using neutron star-black hole mergers

In this work, we study tidal effects that can be observed in gravitational wave transient from neutron star-black hole (NSBH) mergers. It has been well known that neutron star's tidal deformation has imprint of equation of state and hence composition of neutron star matter. Compared to binary neutron star mergers, where we observe a combined effect of two neutron stars by measuring the parameter known as lambda tilde. In neutron star black hole mergers, we measure the tidal deformability of the neutron star component clearly. In this work, we study the feasibility of measuring the bare tidal deformability of a neutron star with planned gravitational detectors in the near and far future. We use the parameter estimation computation framework, Bilby is used to perform Bayesian inference for hundreds of NSBH mergers in current and future GW detector networks and obtain the posterior probability distribution functions of binary parameters. We find that we need at least 20 events to clearly distinguish the equation state of neutron stars in future detectors, including the Einstein Telescope and Cosmic Explorer.
Read the full abstract

High Energy Astrophysical Phenomena

The First Array-Wide Diffuse Flux Search for UHE Neutrinos with the Askaryan Radio Array

The Askaryan Radio Array (ARA) is an ultrahigh energy (UHE) neutrino detector at the South Pole that searches for impulsive broadband radio signals from neutrino-induced particle showers in glacial ice. ARA consists of five autonomous stations with receiving antennas deployed up to 200 m deep and has accumulated the largest livetime of any in-ice radio array. We present the first array-wide diffuse UHE neutrino search using data collected from January 2013 to December 2023, incorporating improved detector characterization and simulation, including data-driven noise and electronics models, revised antenna responses, and updated neutrino interaction and lepton propagation modeling, within a unified framework for event processing, background rejection, and cut optimization across independently operating stations. This search is expected either to identify the first UHE neutrino candidates observed by an in-ice radio detector or to set the most stringent UHE neutrino diffuse flux limits above a few EeV, while informing analysis strategies for under-construction and future radio arrays, such as RNO-G and IceCube-Gen2 Radio.
Read the full abstract

High Energy Astrophysical Phenomena

The DSA/Chronoscope fast radio burst survey: forecasts and science overview

Fast radio bursts (FRBs) are bright extragalactic transients with several mysteries surrounding their origins. Large FRB samples enable accurate measurements of the cosmic matter distribution, in particular on scales $\lesssim 10$ Mpc. These measurements will impact cosmological inference and our understanding of astrophysical feedback, from the circumgalactic medium to cluster scales. Here we forecast the expected yields, redshifts, and host galaxies of FRBs as observed by the Deep Synoptic Array (DSA), and describe the key science cases enabled by the large FRB sample. The DSA will be an interferometer consisting of 1650$\times$6.15 m antennas, operating between 0.7–2 GHz, to be located in Nevada, USA. The Chronoscope backend on the DSA, hereafter DSA/Chronoscope, is designed to search for FRBs across the field of view in real time, enabling the storage of full-polarization voltage data. Extrapolating from existing FRB surveys, we expect roughly $10^4$ FRB detections per year in each of three search sub-bands. Combining across sub-bands, the survey could produce $\sim$ 10$^5$ FRBs over the nominal 5-year DSA survey, assuming Euclidean source counts and a baseline compute backend that can search $6\times10^6$ beams at 1 ms sampling. The well-characterized DSA synthesized beam and deep simultaneous reference images will enable localization precisions of $\lesssim$ 250 milliarcseconds. Key science cases include the use of FRB propagation effects in probing cosmic baryons, and studies of the FRB phenomenon using FRB host galaxies and their local environment, as well as multiwavelength counterparts.
Read the full abstract

All authors Liam Connor, Kaitlyn Shin, Vikram Ravi, Stella Koch Ocker, Casey J. Law, Kritti Sharma, Samuel McCarty, Gregg Hallinan, Shami Chatterjee, James M. Cordes, Dean Howarth, Fabian Walter, Elisabeth Krause, Vishnu Balakrishnan, Alexa C. Gordon, Calvin Leung, Shion Andrew

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

Radio measurements of air showers with the IceCube-Gen2 surface prototype station at the Pierre Auger Observatory

The design of the IceCube-Gen2 observatory, proposed as a next-generation extension of IceCube, includes a surface array consisting of scintillators and radio antennas. In addition to several such detectors already deployed at IceCube's surface array at the South Pole, a complete prototype station including three SKALA antennas has been operating at the Pierre Auger Observatory for several years. This setup has been used to successfully observe radio signals from air showers, demonstrated through coincident observations with the Auger Surface Detector. In this contribution, we present an updated analysis of these radio signals, and compare them to CoREAS simulations using the reconstruction from the Auger Surface Detector as input.
Read the full abstract

High Energy Astrophysical Phenomena

STOKES tables and KY codes

The family of KY codes is an XSPEC-compatible package of fully relativistic models for spectro-polarimetric modeling of X-ray binary systems and active galactic nuclei. Within Kerr space-time, they compute accretion disc line profiles, reflection spectra from discs illuminated by Comptonized coronal radiation, and thermal disc emission including disc-corona interaction and disc self-irradiation. Many of these models include polarimetric computations, and more advanced timing variants compute X-ray, UV and optical reverberation. The primary source of Comptonized X-rays is assumed to be in the lamp-post on-axis or slab above the disc geometry. The KY codes operate at various levels of complexity; some models use analytical approximations for local re-processing, while others rely on detailed numerical tables. For the most accurate treatment of rest-frame re-processing in partially ionized accretion disc atmospheres, the KY codes employ pre-computed spectro-polarimetric tables generated with the STOKES code. These transmission and reflection tables are also designed for standalone use in XSPEC and as inputs to higher-level models beyond KY. The chapter first introduces the STOKES tables in various variants and shows their non-relativistic usage in XSPEC models for (i) an extended source illuminating a distant geometrically thin disc, and (ii) a compact source illuminating a vertically extended equatorial obscurer (e.g., a torus or thick wind). The relativistic KY ray-tracing is then described, followed by an introduction to each KY model. Examples of fits to the IXPE data of X-ray binaries and active galactic nuclei using these models are shown. Together, these models form a comprehensive and publicly available framework for spectro-polarimetric modeling and data fitting of accretion. We provide links to their online repositories.
Read the full abstract