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)

Wednesday 5 August

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

Novel Signatures of Matter-Induced Dark Matter Decay in Large-Volume Neutrino Telescopes

Large-volume neutrino telescopes offer a unique opportunity to search for decaying dark matter through events containing a pair of energetic, highly non-collimated muon tracks emerging from a common vertex. Such events would have negligible Standard Model backgrounds and would constitute a striking signature of new physics. Conventional dark matter annihilation or decay, however, is too strongly constrained to produce an observable rate of such events. We therefore consider scenarios in which an excited dark matter state is extremely long-lived in vacuum but decays much more rapidly in the presence of ordinary matter. We present two realizations of this mechanism. In the first, a long-range scalar field sourced by ordinary matter modifies the dark-sector mass spectrum, kinematically opening the decay $\chi_2 \rightarrow \chi_1 Z'$ near the Earth while leaving it forbidden in vacuum. In the second, the scalar background induces kinetic mixing between a heavy $Z'$ and the photon, greatly enhancing the three-body decay $\chi_2\to\chi_1\mu^+\mu^-$ in matter-rich environments. We calculate the resulting distributions of muon energies and opening angles and show that viable regions of parameter space can yield observable event rates in IceCube, KM3NeT, and other large-volume neutrino telescopes while remaining consistent with existing constraints. We also briefly consider the sensitivity of IceCube to multi-muon events produced by the decays of cosmologically long-lived charged particles with masses $\gtrsim 1$ TeV.
Read the full abstract

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

Emission line formation in scattering dominated media: implications for LRDs

Recent JWST observations of “Little Red Dots” (LRDs) reveal broad and prominent Balmer emission lines. We present a theoretical framework for intrinsic emission line formation and broadening within static, optically thick, scattering-dominated gas envelopes with thermal populations. Using random-walk and diffusion approximations, we derive analytical line profiles for lines forming intrinsically within the scattering medium. We demonstrate that a geometrically thin planar photosphere produces a shallow line profile characterized by a logarithmic plateau and a $v^{-1}$ wing. A radially extended photosphere yields a broken power-law spectrum transitioning from $v^{-\alpha}$ to $v^{-(\alpha+1)}$, with $0 < \alpha < 1$. This is in contrast to a scattering medium external to the line-forming region, which produces an exponential line profile. We show that this broken power-law model can fit the $\mathrm{H}\alpha$ line profiles observed in LRDs. Higher quality spectra may be able to distinguish between intrinsic and extrinsic models for the line broadening in LRDs. In our LTE models, the high contrast between the $\mathrm{H}\alpha$ and continuum flux cannot be explained. Quantitative comparison to LRD spectra requires expanding our models to non-LTE situations.
Read the full abstract

All authors Elisha Modelevsky, Omri Nitzan, Re'em Sari, Eliot Quataert

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

The glow of eternal black holes

We compute the optical appearance of a maximally extended (eternal) Schwarzschild black hole surrounded by a geometrically thin accretion disk. Unlike astrophysical black holes formed by gravitational collapse, the eternal solution contains a white hole (WH) region connected to a past singularity. Following Markov's hypothesis of a limiting density of matter, we replace the past singularity with a spacelike "surface of last scattering" at $r=r_{\rm m}$ and assume that it emits black-body radiation with temperature $T$. Past-directed rays that cross the past horizon terminate on this surface, resulting in a bright spot at the center of the shadow of an eternal black hole. The radial profile of the spot is set by the gravitational frequency shift, with a blueshifted center when the spacelike surface is sufficiently close to the singularity. We use ray tracing to generate an image of the disk and the spot, and we derive a closed-form interferometric signature of the spot. Its visibility is a pure exponential in the baseline length, in contrast to the power-law envelopes of the disk and the photon ring. We constrain the product $r_{\rm m} T$ from the observed 230 GHz fluxes of M87* and Sgr A*. We also consider a Planck-star scenario in which the value of $r_{\rm m}$ is determined on dimensional grounds, and show that it is out of reach of any current or planned facilities. Nevertheless, the darkness of the observed shadows remains a test of whether these black holes are eternal.
Read the full abstract

High Energy Astrophysical Phenomena

Constraining Lorentz invariance violation from the depth of air-shower maximum

Hypothetical Lorentz invariance violation (LIV) remains strongly constrained but has not been excluded as a possible manifestation of new high-energy physics. In subluminal photon-sector LIV, suppressed Bethe$\unicode{x2013}$Heitler pair production modifies the electromagnetic component of extensive air showers. We present a simulation-based toy analysis of this effect using published Pierre Auger fluorescence-detector distributions of the reconstructed depth of air-shower maximum. The modified distributions of the depth of the shower maximum are folded with the Auger detector-response parametrizations and compared with the data using free composition fractions in each reconstructed-energy bin. Within this approach, the comparison gives a sensitivity to LIV compatible with that expected from ultra-high-energy photon-induced showers. We obtain $M_{\rm LIV}>1.5\times 10^{21}\,\text{GeV}$ ($95\%$ confidence level). This result should not be interpreted as a detector-level experimental limit because the detector response and reconstruction are simplified, and the result remains composition dependent. However, the approach presented here opens the way for constraining LIV in hadronic air showers with more detailed analyses.
Read the full abstract

All authors Nickolay S. Martynenko, Grigory I. Rubtsov, Petr S. Satunin, Andrey K. Sharofeev, Sergey V. Troitsky

High Energy Astrophysical Phenomena

Dissecting the Nuclear Structure of NGC 5548 with XRISM. I. Physical Properties of the Highly Ionized Outflows

We present a detailed spectral analysis of an X-Ray Imaging and Spectroscopy Mission (XRISM) observation of the prototypical Seyfert 1 galaxy NGC 5548. XRISM's Resolve microcalorimeter reveals, for the first time, highly ionized outflows in this active galactic nucleus (AGN) through the detection of Fe XXV and Fe XXVI absorption lines in the Fe K band. Modeling the XRISM/Resolve spectrum alongside XMM-Newton Reflection Grating Spectrometer (RGS) data allows us to probe the ionization and kinematic structure of the outflows in this AGN. We identify four distinct ionization components, with ionization parameters log $\xi$ ranging from 0.9 to 3.4. Three of these components are further resolved into two velocity sub-components, demonstrating the multiphase structure of the outflows. The measured outflow velocities span 240 to 2730 km/s. We find a trend of increasing column density with ionization parameter ($\xi$), along with a general pattern of increasing outflow velocity with $\xi$. The XRISM/Resolve spectrum provides a far more detailed absorption measure distribution (AMD) than was previously possible, revealing two distinct slopes above and below $\log\xi \sim 2.6$. A comparison of the Fe XXV absorption line profile with UV absorption lines (C IV and Ly$\alpha$) observed with the Hubble Space Telescope reveals both overlaps and deviations. The XRISM/Resolve results suggest a multiphase, clumpy outflow in NGC 5548, consistent with a "hybrid wind" scenario in which the observed parameter trends arise from multiple origins and driving mechanisms.
Read the full abstract

All authors Missagh Mehdipour, Jon M. Miller, Jelle S. Kaastra, Gerard A. Kriss, Keigo Fukumura, Liyi Gu, Doyee Byun, Xin Xiang, Ehud Behar, Laura W. Brenneman, Elisa Costantini, Maryam Dehghanian, Jacobo Ebrero, Massimo Gaspari, Anna Juráňová, Erin Kara, Chen Li, Junjie Mao, Hirofumi Noda, Anna Ogorzalek, Ioanna Psaradaki, Daniele Rogantini, Sascha T. Zeegers, Keqin Zhao

High Energy Astrophysical Phenomena

Bayesian Inference of fine-features of dense matter EOS from future high-precision data of neutron star radii

Future high-precision X-ray and gravitational wave observatories are expected to measure the radii of neutron stars (NSs) with an accuracy better than about 0.1 km. However, it remains unclear what particular aspects of the Equation of State (EOS) and to what precision they will be better constrained. Within a Bayesian framework using a meta-model EOS and mock high-precision NS data, the posterior probability distribution functions (PDFs) of NS matter EOS parameters for both hadronic and quark phases and the transition between them were recently studied. We report here a few highlights of these studies.
Read the full abstract

All authors Bao-An Li, Xavier Grundler, Wen-Jie Xie, Nai-Bo Zhang

High Energy Astrophysical Phenomena

X-ray Reflection as Diagnostic of Supermassive Black Hole Binary Properties

We investigate correlations between prominent features in the relativistic X-ray reflection spectrum emitted by an accreting supermassive black hole (SMBH) binary with the underlying properties of the binary system. Model-independent measurements of the relativistic Fe K$\alpha$ line ($\sim$6.4 keV) and the Compton reflection hump ($\sim$20-30 keV) are shown to be useful in constraining binary parameters. We compute 24,570 X-ray reflection spectra from two mini-disks attached to SMBHs at a fixed separation of 100 $GM/c^2$ on circular orbits, by varying its mass ratio, spin parameters, inclination, orbital phase and total mass accretion rate. We find that the location of the blue peak in the relativistic Fe K$\alpha$ is a relatively robust diagnostic of the binary inclination, which could be obtained from a single-epoch X-ray spectrum. Given a few epochs of spectra, one may be able to determine the orbital phase of the binary and place constraints on its mass ratio by monitoring the Fe K$\alpha$ centroid. Of all parameters, SMBH spin efects are most subtle and prone to measurement degeneracies. Some markers of high SMBH spin may nevertheless surface in the composite spectrum due to increased radiative efficiency. The approach developed here can be used to place preliminary constraints on binary parameters before a full parametrized X-ray spectral fitting method is available. It complements gravitational wave measurements by the Pulsar Timing Arrays (PTAs) and the Laser Interferometer Space Antenna by providing independent constraints on binary parameters that may be prone to degeneracy (inclination), or otherwise inaccessible (spin for PTAs).
Read the full abstract

All authors Julie Malewicz, Tamara Bogdanović, David R. Ballantyne, Laura Brenneman, Thomas Dauser

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

Ambiguity in matter sector for modified gravity involving $\delta^2 \mathcal{L}_{m}/\delta g^{\mu\nu}\delta g^{\alpha\beta}$ and its implications to astrophysics and cosmology

Matter density ($\rho$) and radial pressure ($p$) are often used as the matter Lagrangian density ($\mathcal{L}_{m}$) because both are thermodynamically consistent and produce the same Einstein field equation (EFE) in general relativity (GR). New gravity models with explicit links between matter and geometry instead involve second-order derivatives of $\mathcal{L}_{m}$ relative to the metric tensor. So, picking either $p$ or $-\rho$ for $\mathcal{L}_{m}$ gives different effective EFEs. This confusion appears because one usually treats the four-velocity ($u_\mu$) and the metric tensor ($g_{\mu \nu}$) as independent. Here, we revisit the basics and offer a consistent framework by relaxing that assumption, thereby making the modified gravity theory independent of the choice of $\mathcal{L}_{m}$. Finally, we test this approach on neutron and quark stars (ultraviolet region) and on cosmological situations with radiation-dominated ($p=\rho/3$) equations of state (infrared region), showing how it clarifies the ambiguity in picking $\mathcal{L}_{m}$ for gravity models.
Read the full abstract

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

Non-Thermal Pressure due to Gas Motions in the Intracluster Medium: Confronting XRISM/Resolve with TNG-Cluster Simulations

Intracluster medium (ICM) gas motions probe cluster assembly, feedback, and non-thermal pressure support, but recent XRISM observations reveal velocity dispersions and non-thermal pressure fractions systematically lower than simulations predict, with extreme systems such as Abell 2029 falling below nearly all simulated clusters. Using the TNG-Cluster simulations, we show that the non-thermal pressure fraction depends sensitively on cool-core state and formation history, and provide a two-scale fitting function capturing both the inner cool-core suppression and outer rise of the radial profile. By forward-modeling mock XRISM observations and comparing them with both projected and intrinsic three-dimensional quantities, we find that azimuthal variations and projection effects contribute to the deficit in the observed velocity dispersion and non-thermal pressure fraction. This bias increases with radius and partially offsets the intrinsic outward rise in the true three-dimensional non-thermal pressure fraction. However, these effects cannot explain the extremely low values of the non-thermal pressure fraction observed in Abell 2029, which fall below approximately the 0th-6th percentiles of the simulated cool-core cluster distribution at every measured radius under both the turbulence-only and turbulence-plus-bulk definitions. The remaining tension points to rare dynamical conditions or missing physics affecting the amplitude of gas motions in current ICM models.
Read the full abstract

High Energy Astrophysical Phenomena

Observation of the Moon and Sun shadows with cosmic rays at an average energy of $\text{7}{\times}\text{10}^\text{17}\,$eV

Interactions of cosmic rays with the Moon and the Sun produce deficits in their arrival-direction distributions relative to an isotropic flux. Such shadows have been observed previously at energies between $10^{12}\,$eV and $10^{16}\,$eV. We report the first observation of the Moon and Sun shadows at cosmic-ray energies larger than about $10^{16}\,$eV (average energy of $7\times10^{17}\,$eV), using data collected by the Pierre Auger Observatory. We employ data from three detector arrays covering $3000\,\text{km}^2$, $27\,\text{km}^2$, and $2\,\text{km}^2$, with spacings of $1500\,\text{m}$, $750\,\text{m}$, and $433\,\text{m}$, respectively. The data amount to over 10.6 million events. The Moon and Sun shadows are detected with a combined significance of approximately $3\sigma$. These observations confirm the pointing accuracy of the Surface Detector of the Pierre Auger Observatory using celestial bodies. From the combined Sun and Moon shadows, we infer an overall angular resolution of $(0.59^{+0.15}_{-0.11})^\circ$.
Read the full abstract

All authors The Pierre Auger Collaboration, A. Abdul Halim, P. Abreu, M. Aglietta, M. Ahmed, I. Allekotte, K. Almeida Cheminant, R. Aloisio, J. Alvarez-Muñiz, A. Ambrosone, J. Ammerman Yebra, L. Anchordoqui, B. Andrada, L. Andrade Dourado, L. Apollonio, C. Aramo, J.C. Arteaga Velázquez, P. Assis, G. Avila, E. Avocone, A. Bakalova, Y. Balibrea, A. Baluta, F. Barbato, A. Bartz Mocellin, O. Batalla Cruz, J.P. Behler, C. Berat, M.E. Bertaina, M. Bianciotto, P.L. Biermann, V. Binet, K. Bismark, T. Bister, J. Biteau, J. Blazek, J. Blümer, M. Boháčová, D. Boncioli, C. Bonifazi, N. Borodai, J. Brack, P.G. Brichetto Orquera, S. Buitink, A. Bwembya, T.R. Caba Pineda, K.S. Caballero-Mora, S. Cabana-Freire, L. Caccianiga, J. Caraça-Valente, R. Caruso, A. Castellina, F. Catalani, G. Cataldi, L. Cazon, M. Cerda, B. Čermáková, A. Cermenati, K. Cerny, J.A. Chinellato, J. Chudoba, L. Chytka, R.W. Clay, A.C. Cobos Cerutti, R. Colalillo, R. Conceição, G. Consolati, M. Conte, F. Convenga, D. Correia dos Santos, P.J. Costa, C.E. Covault, M. Cristinziani, C.S. Cruz Sanchez, S. Dasso, K. Daumiller, B.R. Dawson, R.M. de Almeida, E.-T. de Boone, B. de Errico, J. de Jesús, S.J. de Jong, J.R.T. de Mello Neto, I. De Mitri, D. de Oliveira Franco, F. de Palma, V. de Souza, E. De Vito, A. Del Popolo, O. Deligny, N. Denner, K. Denner Syrokvas, L. Deval, A. di Matteo, C. Dobrigkeit, J.C. D'Olivo, L.M. Domingues Mendes, T. Dominguez, Y. Dominguez Ballesteros, Q. Dorosti, R.C. dos Anjos, J. Ebr, F. Ellwanger, R. Engel, M. Erdmann, A. Etchegoyen, C. Evoli, H. Falcke, G. Farrar, A.C. Fauth, T. Fehler, F. Feldbusch, A. Fernandes, M. Fernández Alonso, B. Fick, J.M. Figueira, P. Filip, A. Filipčič, B. Flaggs, A. Franco, M. Freitas, T. Fujii, A. Fuster, C. Galea, B. García, C. Gaudu, P.L. Ghia, U. Giaccari, M. Giammarco, C. Glaser, F. Gobbi, F. Gollan, G. Golup, P.F. Gómez Vitale, J.P. Gongora, N. González, D. Góra, A. Gorgi, M. Gottowik, F. Guarino, G.P. Guedes, Y.C. Guerra, L. Gülzow, S. Hahn, P. Hamal, M.R. Hampel, P. Hansen, A. Haungs, M. Havelka, T. Hebbeker, C. Hojvat, J.R. Hörandel, P. Horvath, M. Hrabovský, T. Huege, A. Insolia, P.G. Isar, M. Ismaiel, P. Janecek, V. Jilek, K.-H. Kampert, B. Keilhauer, V.V. Kizakke Covilakam, H.O. Klages, M. Kleifges, A. Klingel, J. Köhler, F. Krieger, M. Kubatova, N. Kunka, B.L. Lago, N. Langner, N. Leal, M.A. Leigui de Oliveira, Y. Lema-Capeans, A. Letessier-Selvon, I. Lhenry-Yvon, L. Lopes, M. Mallamaci, S. Mancuso, D. Mandat, P. Mantsch, A.G. Mariazzi, C. Marinelli, I.C. Mariş, G. Marsella, D. Martello, S. Martinelli, O. Martínez Bravo, A. Martínez-Mendez, M.A. Martins, H.-J. Mathes, J. Matthews, G. Matthiae, E. Mayotte, S. Mayotte, P.O. Mazur, G. Medina-Tanco, D. Melo, A. Menshikov, C. Merx, S. Michal, M.I. Micheletti, L. Miramonti, M. Mogarkar, S. Mollerach, F. Montanet, L. Morejon, K. Mulrey, R. Mussa, W.M. Namasaka, S. Negi, L. Nellen, K. Nguyen, G. Nicora, M. Niechciol, D. Nitz, D. Nosek, A. Novikov, V. Novotny, L. Nožka, A. Nucita, L.A. Núñez, S.E. Nuza, J. Ochoa, M. Olegario, C. Oliveira, L. Östman, M. Palatka, J. Pallotta, G. Parente, T. Paulsen, M. Pech, J. Pękala, R. Pelayo, C. Pérez Bertolli, L. Perrone, S. Petrera, T. Pierog, M. Pimenta, M. Platino, P. Privitera, C. Priyadarshi, M. Prouza, K. Pytel, S. Querchfeld, J. Rautenberg, D. Ravignani, J.V. Reginatto Akim, M.Z. Rennó, A. Reuzki, J. Ridky, F. Riehn, M. Risse, V. Rizi, B. Rocha Moldes, E. Rodriguez, G. Rodriguez Fernandez, J. Rodriguez Rojo, S. Rossoni, M. Roth, E. Roulet, A.C. Rovero, A. Saftoiu, M. Saharan, F. Salamida, H. Salazar, G. Salina, P. Sampathkumar, N. San Martin, J.D. Sanabria Gomez, F. Sánchez, F.M. Sánchez Rodriguez, E. Santos, F. Sarazin, R. Sarmento, R. Sato, P. Savina, V. Scherini, H. Schieler, M. Schimp, D. Schmidt, O. Scholten, H. Schoorlemmer, P. Schovánek, F.G. Schröder, J. Schulte, T. Schulz, S.J. Sciutto, M. Scornavacche, A. Sedoski, S. Sehgal, S.U. Shivashankara, G. Sigl, K. Simkova, F. Simon, R. Šmída, S. Soares Sippert, P. Sommers, S. Stanič, J. Stasielak, P. Stassi, S. Strähnz, M. Straub, T. Suomijärvi, A.D. Supanitsky, Z. Svozilikova, Z. Szadkowski, F. Tairli, A. Tapia, C. Taricco, C. Timmermans, O. Tkachenko, P. Tobiska, C.J. Todero Peixoto, B. Tomé, A. Travaini, P. Travnicek, C. Trimarelli, M. Tueros, M. Unger, R. Uzeiroska-Geyik, L. Vaclavek, M. Vacula, I. Vaiman, J.F. Valdés Galicia, L. Valore, P. van Dillen, E. Varela, V. Vašíčková, A. Vásquez-Ramírez, D. Veberič, I.D. Vergara Quispe, S. Verpoest, V. Verzi, J. Vicha, S. Vorobiov, J.B. Vuta, A.A. Watson, A. Weindl, M. Weitz, L. Wiencke, H. Wilczyński, B. Wundheiler, B. Yue, A. Yushkov, E. Zas, D. Zavrtanik, M. Zavrtanik

High Energy Astrophysical Phenomena

Population synthesis and detection prospects for Galactic long-period transients with LISA

The recently-discovered long-period radio transients represent a puzzling new class of astrophysical sources, some of which are thought to be compact binary systems emitting a pulse once per orbit. If this interpretation is correct, the orbital period–and therefore the gravitational-wave frequency–is directly encoded in the radio signal, which typically lies in the millihertz band. In this work, we explore whether these systems can be detected by the Laser Interferometer Space Antenna (LISA). Assuming that this phase-locking between radio pulses and orbital motion applies broadly across the population, we construct synthetic source catalogues informed by both observations and theoretical models of related systems, such as cataclysmic variables. We estimate that between $\sim 0.05\%$ and $\sim 6\%$ of long-period transients will be detectable within four years of observation with LISA, depending on astrophysical assumptions. For detectable systems, we find injected frequencies are recoverable to one part in $\sim 10^{5}$, amplitudes to within a factor $\sim 2$, and the sky positions to within $\sim$30 square degrees. Our results demonstrate that gravitational-wave observations can provide direct evidence for the binary nature of these sources and, importantly, can guide future radio surveys by predicting pulse periods, sky positions, and orbital properties. The mismatch between extracted frequency derivatives and that imposed by gravitational-wave decay is also computed to show that the likelihood of electromotive losses driving orbital evolution can be assessed for each detectable candidate.
Read the full abstract

High Energy Astrophysical Phenomena

Type I X-Ray Burst Models With Rotation

Type I X-ray bursts are powered by unstable thermonuclear burning on the surface of accreting neutron stars in close binary systems. These brief X-ray flashes, with light curves featuring rise times of $1-10$ s, durations of $10-100$ s, and recurrence periods of hours to days, represent the most frequent stellar explosions in our Galaxy. With typical energies of $\sim 10^{39}-10^{40}$ erg, they rank among the most powerful astrophysical transients after supernovae and classical novae. To date, roughly 120 bursting X-ray binaries have been identified in the Milky Way. Several studies have been conducted to characterize the dynamics of these events, with emphasis on reproducing the observed recurrence periods and light curve shapes. In this paper we show, for the first time, that rotation is a key factor shaping the properties of Type I X-ray bursts in rapidly spinning systems. The inclusion of centrifugal forces, together with a suite of rotationally-induced mixing mechanisms, such as meridional circulation and shear-induced turbulent diffusion, reduce surface gravity, shortening the recurrence times and lowering burst energies. Rotation also modifies the extent of the nuclear activity during these events and affects the morphology of their light curves, which are distinctly broader for rapidly rotating neutron stars.
Read the full abstract

High Energy Astrophysical Phenomena

A Real-Time Jet Laboratory in Swift J1727.8-1613

We present the results of our intensive VLBI campaign on the black-hole low-mass X-ray binary (LMXB) Swift J1727.8-1613 during its 2023-2024 outburst. We observed the repeated quenching and re establishment of the highly-extended continuous jet during several transitions between hard-intermediate and soft-intermediate states, and the repeated ejection of transient jets. Using time-dependent visibility model fitting, we tracked the motion of nine discrete jet knots, obtaining some of the most precise measurements of transient jet proper motions and ejection dates in an LMXB. These ejecta were only detectable for a short time with VLBI, and some showed rapid intra-observation flux density variability that was not captured in image reconstructions. For the first time, we use time-dependent visibility modelling to fit a piecewise model for the jet knot flux densities, allowing us to create complex, non-parametric light curves of their intra-observation variability. We observed the launching of multiple ejecta across several state transitions, however, we could not identify a consistent signature of jet ejection in the available X-ray intensity or hardness data. We constrained the intrinsic speeds and bulk Lorentz factors of the jet knots, finding that Swift J1727.8-1613 launched both mildly relativistic and highly relativistic ejecta throughout its outburst. We used their proper motions to constrain a posterior distribution for the maximum inclination angle of the jet axis, which had 50th, 84th, and 99th percentiles of 40{\deg}, 50{\deg}, and 66{\deg}, respectively. These unique observations of the repeated ejection of transient jets by a single LMXB reveal that fixed parameters such as black-hole mass, black-hole spin, and spin-orbit misalignment do not uniquely determine the varying properties of transient jets, particularly their speeds and Lorentz factors.
Read the full abstract

All authors Callan M. Wood, James C. A. Miller-Jones, Arash Bahramian, Steven J. Tingay, Sara E. Motta, Hongmin Cao, Thomas D. Russell, Francesco Carotenuto, Pikky Atri, Diego Altamirano, Alexandra J. Tetarenko, Rob Fender, Elmar Körding, Dipankar Maitra, Sera Markoff, David M. Russell, Gregory R. Sivakoff, Roberto Soria, Valeriu Tudose

Tuesday 4 August

High Energy Astrophysical Phenomena

Multimessenger Constraints on the Central Engine in SN 2023ixf: Neutron Star or Black Hole?

The physical trigger of core-collapse supernova (CC-SNe) is a central theme in modern surveys of the transient Universe. A multi-messenger approach may provide a unique window to CC-SNe produced by a newly born neutron star (NS) or black hole (BH). Here we apply joint electromagnetic-gravitational wave (EM-GW) analysis to the nearby event SN 2023ixf using a model-independent broadband search with a GW-energy detection threshold of a few times $10^{-5} M_\odot c^2$ for chirps and $10^{-2}$ $M_\odot c^2$ for noise-like emission below 1\,kHz. No joint LIGO Hanford-Livingston (H1-L1) signal associated with SN 2023ixf is observed over a 2.5 day window proceeding the first discovery. Together with EM photometric-spectroscopic data, we interpret this null result to indicate the absence of BH formation, in part by mass-scaling of the BH powering GW170817B following GW170817. This multi-messenger outcome supports, hereby, scenarios consistent with NS formation, consistent with a progenitor mass likely below $\sim 20 M_\odot$ estimated by EM interpretations.
Read the full abstract

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

A survey of ultra-compact high-state AM CVn binaries with ZTF and Gaia: New discoveries and observational constraints on Galactic space density

Ultra-compact AM CVn binaries in the high-state of stable mass transfer ($\dot{M}\gtrsim10^{-9}M_\odot/\rm{yr}$) are expected to be among the loudest persistent sources for upcoming space-based GW observatories. These systems typically have orbital periods $P_{\rm{orb}}\lesssim30$ minutes. We present a systematic search for high-state AM CVn binaries in the Milky Way by targeting a region within the Gaia color-magnitude diagram, combined with ZTF time-domain photometry. We focus on variable targets within that sample with 5-30 minute periods. Our survey discovered three new high-state AM CVn binaries (ZTF J1840-1742, ZTF J2007-0527 and ZTF J2111+3158) with orbital periods of 16.86, 18.63, and 16.71 minutes, and recovered three known systems. We obtain high-speed photometry and phase-resolved spectroscopy to confirm their nature. Their spectra show helium emission lines with no detectable hydrogen. In two targets, the lines are double-peaked, and Doppler tomograms confirm an accretion disk. We estimate mass accretion rates of $\sim10^{-9}M_\odot/\rm{yr}$, and $3\sigma$ X-ray luminosity upper limits of $\sim10^{33}\rm{erg/s}$. Based on this sample, we infer a local space density of high-state AM CVn population in the Milky Way of $\sim(1.0-2.7)\times10^{-8}\rm{pc}^{-3}$ (for disk scale heights $h_z=300-200$ pc), representing 2-5% of the total Galactic AM CVn density. Their birth rate, $\sim(2.4-4.3)\times10^{-4}\rm{yr}^{-1}$ (for $h_z=300-200$ pc), is consistent with the total AM CVn birth rate, implying most systems entering the high-state phase survive and evolve to longer orbital periods. This local space density suggests LISA and TianQin will detect about 36% and 19% of the total Galactic population, respectively, during their nominal 4-yr missions ($S/N\ge5$). The Vera C. Rubin Observatory's LSST will detect about 34-44% of these systems over its 10-yr survey ($m_r\lesssim26.9$ mag).
Read the full abstract

All authors Ilkham Galiullin, Antonio C. Rodriguez, Kareem El-Badry, Valery Suleimanov, Vladislav Dodon, Askar Sibgatullin, Warren R. Brown, Kevin Burdge, Jan van Roestel, Edo Berger, V. Ashley Villar, Ilaria Caiazzo

High Energy Astrophysical Phenomena

Profile Reconstruction from Temporally Stable Emission Components for Timing PSR J1713+0747

The assumption of long-term pulse-profile stability underpins high-precision pulsar timing and forms the basis of pulsar timing array experiments. However, several millisecond pulsars exhibit temporal profile variability that can introduce systematic biases in pulse time of arrival measurements and compromise timing precision. We present a profile-domain analysis of PSR J1713+0747 at low radio frequencies, in the 300-500 MHz band, using upgraded GMRT observations for the Indian Pulsar Timing Array experiment. We model frequency-resolved pulse profiles using a Bayesian Gaussian decomposition framework in which individual Gaussian components are associated with persistent emission regions through informative phase priors that permit modest temporal variations. By tracking the evolution of the decomposed components across observing epochs and frequency sub-bands, we identify central Gaussian components that remain precisely localized despite changes in the integrated pulse morphology. We then reconstruct pulse profiles with realistic noise using these central components and perform timing analysis. Our approach provides a physically motivated framework for mitigating pulse-profile variability and offers a generic methodology for recovering robust timing information from pulsars exhibiting profile evolution.
Read the full abstract

All authors Shaswata Chowdhury, M. A. Krishnakumar, Sharika Dhakappa, Vidit Singh, Debabrata Deb, Jyotijwal Debnath, Kaustubh Rai, Pratik Tarafdar, Abhimanyu Susobhanan, Churchil Dwivedi, Bhal Chandra Joshi, Shantanu Desai, Neelam Dhanda Batra, Jaikhomba Singha, Himanshu Grover, Manjari Bagchi, Mayuresh Surnis, Avinash Kumar Paladi, Aman Srivastava, Arul Pandian B., Suruj Jyoti Das, Jibin Jose, Kuldeep Meena, Sushovan Mondal, K Nobleson, Keitaro Takahashi, Hemanga Tahbildar, Kunjal Vara, Zenia Zuraiq

High Energy Astrophysical Phenomena

From Multimessenger Inference to Simulations: A Ranked Ensemble of Finite-Temperature Equations of State

We construct a set of microphysical, finite-temperature equations of state (EOSs) for numerical simulations of neutron star mergers and core-collapse supernovae that is consistent with modern constraints from nuclear theory and multimessenger astronomy and systematically spans the posterior distribution of allowed EOSs. The EOSs are based on a simplified Skyrme functional whose inputs are nuclear matter saturation properties, extended to supra-saturation densities through the speed of sound at a set of reference densities. The models are assigned continuous likelihood weights from chiral effective field theory and perturbative quantum chromodynamics calculations, the gravitational wave signal GW170817, the complete set of NICER mass-radius measurements, and the Shapiro-delay mass measurement of the radio pulsar J0348. From the resulting catalogue of 5.2 million EOSs, the posterior yields $R_{1.4} = 11.8^{+0.8}_{-0.7}$ km, $\Lambda_{1.4} = 334^{+193}_{-113}$, and $M_{\rm max}^{\rm TOV} = 2.18^{+0.22}_{-0.14}\,M_\odot$ (medians with 90% credible intervals). From this posterior we select a ranked 12-member ensemble, headed by a fiducial, central EOS, whose members are individually plausible while jointly bracketing the posterior spread of neutron-star observables. For all ensemble members we generate general-purpose finite-temperature tables with the SROEOS code, each accompanied by nucleon effective-mass variants that bracket the ensemble's thermal-sector uncertainty, to be made publicly available upon publication.
Read the full abstract

All authors Max Jacobi, Giulia Huez, Sebastiano Bernuzzi, David Radice

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

Column-Density Estimation from the Equivalent Widths of Absorption Doublets

Na I D absorption is widely used to estimate the dust extinction toward Galactic and extragalactic sources, but the inferred column densities are unreliable once the lines saturate. We present an analytical framework for estimating column densities from the equivalent widths ($EW$s) of absorption-line doublets. The method requires only that the two doublet members be measured separately – the velocity structure need not be resolved – so a resolving power $\lambda/\Delta\lambda\gtrsim10^{3}$ suffices for Na I D, two orders of magnitude below velocity-resolved methods. For the Gaussian curve of growth, the classical inversion reduces to a universal saturation correction depending only on the doublet ratio, $R=EW_2/EW_1$: at fixed $R$, the inferred column density is linear in the measured $EW$. The formulation yields the exact Gaussian inversion, $N_{\tilde\kappa}$, which provides the Doppler parameter; a closed-form second-order approximation; and a strict lower bound, $N_{\rm min}$, valid for an arbitrary velocity structure. We recommend reporting $N_{\rm min}$: an optimal, assumption-free lower bound, within a few per cent of $N_{\tilde\kappa}$ for $R\gtrsim1.4$ and at most a factor $3.5$ below it at $R=1.1$. The ratio $N_{\tilde\kappa}/N_{\rm 2nd}$ provides a per-object saturation diagnostic. We validate the framework against published profile-fitting Na I D column densities toward Type Ia supernovae, using only the integrated $EW$s. An explicit closed-form formula converts the two $EW$s into an extinction estimate through the observed $\log N(\mathrm{Na\,I})$–$A_V$ relation. A public implementation with full uncertainty propagation accompanies the paper.
Read the full abstract

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

The Most Dangerous Seed: Nonlinear Optimal Perturbations in Rayleigh-Taylor Instability

Long-term instabilities in astrophysical fluids are inherently nonlinear, where even small-amplitude perturbations can trigger dramatic instability. However, owing to the complex interactions among non-normal modes, the perturbation structures responsible for the greatest growth remain poorly understood. In this paper, we employ the nonlinear optimization method known as the conditional nonlinear optimal perturbation (CNOP) to identify the most dangerous initial velocity perturbation, i.e., the perturbation that maximizes the kinetic energy growth in the two-dimensional compressible Rayleigh-Taylor instability in astrophysical hydrodynamics. Compared with random perturbations, the optimal perturbation forms a coherent wave-packet structure localized around the density interface. We investigate its dependence on spatial resolution and optimization time horizon through two sets of numerical experiments. For a fixed optimization time, increasing the spatial resolution produces a more sharply localized wave packet, whereas for a fixed spatial resolution, increasing the optimization time causes the wave packet to become progressively more dispersed. Furthermore, we analyze the optimal perturbations in Fourier space using the fast Fourier transform (FFT), which provides a clearer characterization of the spectral distribution. Higher-resolution simulations concentrate most of the perturbation energy into only a few dominant modes, while longer optimization times distribute the energy over a broader range of modes. These results indicate that short optimization time horizons involve relatively weak modal interactions and remain closer to the linear regime, whereas longer optimization times enhance nonlinear modal interactions, broaden the spectral distribution, and reduce the predictive capability of linear stability theory.
Read the full abstract

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

The most extreme high-z radio quasars in the eRASS:1 X-ray survey

We present the selection and identification of high-z radio quasars from a combination of X-ray (eROSITA All-Sky Survey), optical/NIR (DECam Local Volume Exploration), and radio (Rapid ASKAP Continuum Survey) surveys. After building a sample of 68 new high-z quasar candidates, we performed follow-up spectroscopic observations on 46 sources, confirming 39 to have redshifts z>3.5. For the subset of radio quasars at z>4, the 14 new sources identified here represent an increase of $\sim$65% in the known population of objects above the same flux limits in the same area. Using X-ray-to-optical relative intensity, we identified blazars in our sample, which we then used to compare the total number of z>4, X-ray selected blazars to the predictions of a fractional IC/CMB model previously proposed in the literature. Overall we find a relatively good agreement, even considering the potential biases in the redshift distribution from the incompleteness of the spectroscopic follow-up and the blazar classification criteria. The sources presented in this work represent a sample well suited for investigating the evolution of the most extreme systems in the early Universe, including their relativistic jets and accretion properties
Read the full abstract

All authors L. Ighina, L. Palmieri, L. N. Martínez-Ramírez, A. Caccianiga, T. Connor, A. Moretti, B. Arsioli, Y. Beletsky, E. Marini, A. Rossi

High Energy Astrophysical Phenomena

First fast radio burst search campaign at the Argentine Institute of Radio Astronomy

Fast radio bursts (FRBs) are intense millisecond-duration radio transients of extragalactic origin whose physical nature remains under active investigation, and which also serve as probes of the intergalactic medium. We report on the first FRB search campaign carried out at the Argentine Institute of Radio Astronomy (IAR) between December 2024 and March 2026, targeting nearby galaxy superclusters in the southern sky. We observed fields in the Ophiuchus, Shapley, and Sculptor/Phoenix supercluster regions with one of the two 30~m antennas of the IAR, using a ROACH-based backend with 400~MHz of bandwidth centred at 1400~MHz and a time resolution of 41–82~$\mu$s, for a total net observing time of 212~h. The data were searched for dispersed single pulses with PRESTO in the dispersion measure range $100 \leq \mathrm{DM} \leq 500$, and candidates were classified with the FETCH machine learning classifier. The pipeline was validated on archival Parkes data containing known FRBs and on synthetic bursts injected into IAR observations. One FRB candidate, FRB~20251018, was identified on 18 October 2025 in an observation pointed towards the galaxy cluster A2870, in the Phoenix supercluster, with a dispersion measure of $243$, a signal-to-noise ratio of 8.2, and a FETCH probability of $p=0.99$. To the best of our knowledge, this would be the first FRB detected from South America. A set of more marginal candidates is also presented. These results demonstrate the capability of the IAR antennas to detect FRBs and support the continuation and extension of the monitoring campaign, including coincident dual-antenna observations and cross-matches with gravitational-wave events and electromagnetic transients.
Read the full abstract

All authors Carlos O. Lousto, Harsh Prajapati, Ezequiel Zubieta, Susana B. Araujo Furlan, Guillermo Gancio, Joaquín Pelle, Gustavo E. Romero, Federico García, Santiago del Palacio, Nathan Cahill

High Energy Astrophysical Phenomena

Probing Axion Like Particle-Proton Interactions through the Interplay of Dark Matter and Cosmic Rays in Large-Scale Structure

Axion-like particles (ALPs) provide a well-motivated dark matter candidate whose interactions with Standard Model particles may generate observable cosmological signatures. We investigate the diffuse $\gamma$-ray emission produced by the scattering process between ALP dark matter and cosmic-ray (CR) protons inside dark matter halos, and explore its detectability through the cross-correlation between the unresolved $\gamma$-ray background (UGRB) measured by Fermi Large Area Telescope (LAT) and the Two Micron All Sky Survey Redshift Survey (2MRS) galaxy catalog in the local Universe ($z\leq0.1$). To model the cosmological $\mathrm{GeV}$-$\mathrm{PeV}$ CR proton population, we develop a phenomenological bottom-up prescription that links CR injection to the star formation rate through the supernova rate, while the steady-state CR density is determined by diffusive escape. Our forecasts indicate that UGRB-galaxy cross-correlation measurements can probe previously unexplored regions of the ALP parameter space over a broad range of ALP masses. More generally, this work presents a general framework for cosmological observables generated by interactions between dark matter and an additional astrophysical field, opening a new avenue for indirect dark matter searches using large-scale structure observations.
Read the full abstract

High Energy Astrophysical Phenomena

Iron K$\alpha$ signatures from accretion disks around fermionic dark matter cores

The fluorescent iron line and its broadening due to relativistic effects are excellent probes to study the inner part of an accretion disk and the space-time geometry near the compact object. We investigate the iron K$\alpha$ line profile within the extended RAR model, which describes a fermionic dark matter distribution on galaxy scales. The most general solutions are characterized by a compact and highly degenerate core able to mimic the central black hole, transitioning into an extended halo composed of the same particles. We aim to contrast the resulting line morphologies in this scenario with those predicted by the standard Kerr black hole paradigm. Special attention will be given to MCG-06-30-15 galaxy. We compute the line profile using the numerical ray-tracing code, Skylight. We consider two distinct configurations for the emissivity of the cold accretion disk: an irradiation profile based on the lamp-post corona prescription, and a phenomenological power-law profile. The resulting profiles exhibit a diverse phenomenology. In particular, the most compact fermion cores produce a line broadening comparable to that observed in rapidly rotating black holes. The presence of emitting matter at radii smaller than a gravitational radius yields distinctive spectral features that are entirely absent in the black hole scenario. For MCG-06-30-15 galaxy, we find a good agreement with the observed broad features of the iron line profile, provided the compactness of the fermion core is close to critical. These results reinforce the need for independent black hole spin measurements. Combined with such constraints, iron-line spectroscopy may provide a powerful observational tool to distinguish black holes from alternative compact solutions, in particular compact fermionic dark-matter cores.
Read the full abstract

All authors V. Crespi, F. L. Vieyro, C. R. Argüelles, J. A. Rueda

High Energy Astrophysical Phenomena

Thermal evolution of neo-neutron stars. II. Temperature-dependent crusts

The neo-neutron star phase is an intermediate stage in the evolution from a hot, lepton-rich compact object, which is formed in the aftermath of a core-collapse supernova or an accretion-induced collapse of a white dwarf, to a cold deleptonized neutron star (NS). Alternatively, this phase can occur after a binary NS merger if the final compact object does not collapse immediately into a black hole. Radial temperature profiles provided by numerical simulations of proto-NS evolution suggest that the composition and mechanical structure of the star's core at the beginning of the neo-NS phase are approximately identical to those of a cold NS. In contrast, the composition and mechanical structure of the outer layers (i.e., the crust) are expected to change as the star cools. The thermal evolution of neo-NSs with temperature-dependent (inner) envelopes was considered by Beznogov et al. (2020). Here, we investigated what happens if the equation of state (EoS) of the entire crust, including its composition, depends on temperature. We analyzed the thermal, mechanical, and chemical evolution of neo-NSs by further developing and employing neo-NSCool, our NS thermal evolution code. First, we show that thermally dripped neutrons slightly slow down the cooling during the early thermal relaxation stage. Then, we demonstrate that the EoSs featuring exotic light species in the deepest layers of the crust result in significantly slower cooling in the photon-cooling era compared to the EoSs that favor massive nuclei. The crust's composition also impacts the crystallization process as well as the contraction of the crust during cooling. The EoSs with light nuclei cause the innermost layers to remain liquid for a longer time than the EoSs with heavy nuclei.
Read the full abstract

All authors Sagnik Chatterjee, Mikhail V. Beznogov, Dany Page, Adriana R. Raduta

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

Gravitational Enstrophy: Local Geometric Origin and Energy-Cascade Constraints

Two-dimensional incompressible fluids conserve energy and enstrophy in the inviscid limit. Fjørtoft's argument uses this pair to constrain nonlinear energy transfer. We examine the corresponding question in General Relativity using the magnetic Weyl functional $\mathcal Z=\int B_{ab}B^{ab}\sqrt{\gamma}\,d^3x$. Its transverse radiativ e part carries an additional frequency-squared weight relative to gravitational-wave energy: $\delta\mathcal Z_k\simeq\omega_k^2W_k$ in the radiation zone, with $W=8\pi G E_{\rm GW}$. We derive its balance law from the Bianchi identities, obtaining the algebraic cancellations, bulk sources and boundary fluxes. Conservation of positive e nergy and curvature norm through nonlinear transfer gives a Fjørtoft constraint when their weights are ordered. Resonant weak-wave interactions generally change the c urvature norm; $2\leftrightarrow2$ dynamics instead conserves energy and wave action. We examine the competition between transfer and linear losses in Kerr and AdS backg rounds, where long-lived or confined modes provide time for nonlinear evolution. In the hydrodynamic regime of AdS$_4$, a fluid-weighted magnetic Weyl functional inherit s fluid enstrophy conservation in the inviscid limit. We compute its proportionality coefficient to boundary enstrophy in closed form at leading gradient order, giving a bulk diagnostic for holographic turbulence.
Read the full abstract

High Energy Astrophysical Phenomena

Revisiting the XMM-Newton Observations of the Galactic Microquasar SS 433: Implications for the Origin of the Ultrahigh-Energy Emission Detected by LHAASO

Recently, the Large High Altitude Air Shower Observatory (LHAASO) detected ultrahigh-energy (UHE; photon energy E>100TeV) $\gamma$-ray emission toward SS 433, the microquasar embedded in the W50 nebula, making it a promising Galactic PeVatron candidate. We reanalyze the archival XMM-Newton observations covering the bipolar jets and the thermal X-ray shell north of SS 433, and derive spatially resolved profiles of the nonthermal X-ray intensity and photon index along both jets. The jet emission softens with distance from the source, implying a correspondingly evolving electron population. In particular, a hard electron component appears close to the jet bases, which can account for the UHE emission from SS 433 via inverse Compton radiation if the magnetic field remains approximately uniform along the jets. The result, however, is highly sensitive to the magnetic field profile. For flux-conserving configurations in which the field decreases as the jet expands, the stronger field required in the inner regions may reduce the number of X-ray-emitting electrons and suppress their inverse Compton emission. Furthermore, electron transport calculations show that injection only at the jet bases cannot reproduce the observed intensity and spectral evolution, particularly the downstream re-brightening features, indicating additional particle injection and/or re-acceleration within the jets.
Read the full abstract

High Energy Astrophysical Phenomena

General-relativistic structure of two-component quantum dark fermion stars

We develop a general-relativistic framework for two-component quantum dark fermion stars: equilibrium configurations of two degenerate fermion species governed by gravity, a Yukawa-mediated dark fifth force and a globally relevant Bohm quantum-pressure correction. The treatment retains the full covariant form of the nonlinear Klein–Gordon equation in the Schwarzschild interior, with closure relations valid at arbitrary compactness, from the ultralight baseline up to densities of order $0.16\,\mathrm{fm}^{-3}$ at which relativistic scalar densities and self-consistent effective fermion masses become unavoidable. Two Lagrangian parameters, the dark-fermion mass and the dimensionless ratio between the Yukawa channel and gravity, together fix the equilibrium structure; a joint measurement of mass and radius therefore constrains the dark-sector microphysics directly from gravitational-wave observables. For particle masses in the band $10^{-11}$–$10^{-10}$ eV, the configurations exhibit radii of $3$-$24\,\mathrm{km}$ and compactness in the range $0.14$-$0.34$, behaving as dual mimickers: at moderate compactness they overlap with neutron stars in mass, radius, and inspiral frequency; at the most compact end ($\kappa\sim 0.34$, $R_T/R_S\simeq 1.5$) they cross the photon sphere and could masquerade as low-mass black holes in the mass-gap region. Tidal-deformability measurements discriminate against both populations: the dimensionless tidal deformability is measurably smaller than the neutron-star value yet remains non-zero, unlike that of a genuine black hole. These objects populate the sensitivity bands of LISA, the Einstein Telescope, and Cosmic Explorer, producing astrometric microlensing signatures individually resolvable by Gaia. The framework thus provides a reproducible and falsifiable template for constraining dark-sector properties through multi-messenger observations.
Read the full abstract

High Energy Astrophysical Phenomena

A JWST redshift for the host galaxy of EP250207b of z=3.2: a collapsar origin is viable

We present James Webb Space Telescope (JWST) and Hubble Space Telescope (HST) observations of the field of the fast X-ray transient (FXT) detected by Einstein Probe, EP250207b, to resolve any ambiguity about the host galaxy and redshift of the FXT. EP250207b was originally associated with a nearby galaxy at z=0.082, based on its low chance alignment probability, and a binary neutron star merger origin was proposed. However, we report the detection of a background galaxy at z=3.2 at the location of EP250207b. Assuming this galaxy is the actual host galaxy, the rest-frame energetics and timescales of the event change. Furthermore, the available data are not able to rule out the presence of a supernova associated with EP250207b if at this redshift. We model the X-ray, optical, near-infrared and radio light curves using a tophat jet model implemented in Redback and find that they are consistent with an on-axis gamma ray burst afterglow. The energetics and host galaxy properties do not allow us to distinguish between a collapsar and a merger driven event.
Read the full abstract

All authors Agnes P. C. van Hoof, Peter G. Jonker, Andrew J. Levan, Nial R. Tanvir, Franz E. Bauer, Joe Bright, Francesco Carotenuto, Ting-Wan Chen, Ashley Chrimes, Gregory Corcoran, Laura Cotter, Joyce N. D. van Dalen, Rob A.J. Eyles-Ferris, Morgan Fraser, Daniele B. Malesani, Daniel Mata Sánchez, Antonio Martin-Carrillo, Paul O'Brien, Francesca Onori, Jonathan Quirola-Vásquez, Maria E. Ravasio, Andrea Rossi, Javi Sánchez-Sierras, Nikhil Sarin, Steve Schulze, Hui Sun, Manuel A.P. Torres

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

From Størmer to Schwarzschild: Analytical Dynamics of Charged Particles in a Dipole Magnetosphere

We investigate the conservative dynamics of relativistic charged particles in a Schwarzschild spacetime threaded by an externally supported dipolar magnetic field. The gravitational sector is treated exactly in isotropic coordinates, while the electromagnetic field is described by the leading asymptotic term of the relativistic exterior dipole solution. This formulation provides an analytically tractable relativistic extension of the classical Størmer problem and allows the circular-orbit and stability structure to be studied explicitly. For the outward-Lorentz branch, characterized by a positive signed magnetic coupling for positive azimuthal motion, we find a non-trivial competition between radial and vertical stability. Weak magnetic coupling shifts the radially marginal orbit inward, whereas at stronger coupling vertical instability becomes the limiting mechanism and the inner boundary of stable circular motion moves outward. The associated vertical epicyclic mode progressively softens and vanishes at the vertical stability boundary. The marginally bound sequence exhibits a distinct strong-coupling behaviour, illustrating that energetic and stability boundaries need not evolve together. The radial and vertical epicyclic frequencies also generate characteristic commensurabilities that may provide natural sites for nonlinear orbital coupling. Because the magnetic field is treated as a test field, these charged-particle effects leave the vacuum Schwarzschild null-geodesic structure unchanged. Magnetic signatures in compact-object observations must therefore arise through the dynamics and radiation of charged matter, plasma propagation effects, or physics beyond the test-field approximation. The framework developed here provides an analytical basis for studying nonlinear dynamics, radiation reaction, and more realistic compact-object magnetospheres.
Read the full abstract

High Energy Astrophysical Phenomena · Cross-listed from physics.plasm-ph

Radiation reaction in the classical relativistic Størmer problem

We extend the classical relativistic Størmer problem by incorporating radiation reaction through the Landau–Lifshitz formulation, thereby providing a self-consistent description of dissipative charged-particle motion in a static dipole magnetic field. We derive the complete dimensionless equations and distinguish them from a reduced drag-only model that preserves the exact energy-loss law while omitting directional effects associated with magnetic field gradients. For planar motion, the reduced system yields exact instantaneous evolution laws for the particle energy and canonical angular momentum, together with averaged transport equations for regular bound librations. When the motion is constrained to the instantaneous circular branch, the secular evolution can be integrated in closed form and approaches a simple large-radius power law. This analytical solution provides a useful benchmark, although the circular branch is radially unstable and therefore does not describe generic nearby trajectories. Numerical integrations further illustrate the nonuniform dissipative deformation of planar rosette-like orbits. In three dimensions, the complete Landau–Lifshitz force produces local exponential damping of small vertical perturbations, with the leading contribution arising from the field-gradient term absent from the reduced model. These exact, averaged, conditional, and local results establish a controlled analytical framework for studying radiation-driven phase-space transport in strongly inhomogeneous magnetic fields and provide a foundation for future global simulations, kinetic descriptions, and calculations of the associated electromagnetic emission.
Read the full abstract