Peer Revenue Past citations are no guarantee of future returns

Papers

Fresh from arXiv, typeset properly. Highlighted words trade on the market.

Physics Astrophysics

astro-ph Watch Astrophysics (sign in to use)

Thursday 23 April

Cosmology and Nongalactic Astrophysics · Cross-listed from gr-qc

Relativistic frequency shifts in gravitational waves from axion clouds

Superradiant instability of ultralight bosons can produce clouds around rotating black holes, whose continuous gravitational wave (GW) emission is a promising observational target. Precise predictions of the signal frequency and its evolution are essential for detecting such continuous GWs. For axions, self-interactions can populate multiple superradiant modes via nonlinear couplings, and GW emission can occur through various channels. To calculate the frequency shifts of GWs emitted through these channels, we employ relativistic perturbation theory based on a bilinear form. We apply this framework to self-interaction effects for the first time, and also revisit the treatment of the self-gravity contribution. Our results provide a simple and unified framework for calculating frequency shifts, including cases in which multiple modes are excited, and are relevant for next-generation GW observations.
Read the full abstract

High Energy Astrophysical Phenomena

Geometry, Not Calorimetry, Drives the Radio/Infrared/Gamma-Ray Correlation

We investigate whether the observed radio-infrared-$\gamma$-ray correlation in star-forming galaxies is a geometric effect rather than a signature of local cosmic-ray (CR) calorimetry. Using the GALPROP framework, we generate synthetic observations for external viewers from a grid of 3D Milky Way models with varied CR source, gas, interstellar radiation, and magnetic field distributions, all normalised to reproduce local CR data. We find that a tight, quasi-linear correlation arises naturally from line-of-sight integration through the extended, radially-structured disc, even when local calorimetry is absent. The correlation's properties depend strongly on viewing geometry, preserving its form under moderate inclination but breaking down in edge-on views where galactic components are stratified. We conclude that the correlation is primarily an emergent property of geometric projection and viewing angle, not a direct tracer of local calorimetry. While these findings directly apply to systems with similar structure and star formation properties, their extension to the diverse population of galaxies that form the global radio-infrared-$\gamma$-ray correlation requires further investigation. This geometric perspective implies that the scatter in these relations can be a powerful diagnostic of galactic structure and CR escape.
Read the full abstract

Astrophysics of Galaxies

Early metal-enriched baryon cycling before the midpoint of cosmic reionization

Models predict that chemical enrichment and gas redistribution should begin rapidly once star formation starts, but direct constraints at the earliest epochs have been scarce. Here we show that metal-enriched gas in multiple ionic phases was already present around galaxies before the midpoint of cosmic reionization. Using JWST/NIRSpec rest-frame ultraviolet spectroscopy from SPURS, we detect blueshifted metal absorption in three galaxies at $7.2<z<9.3$. The detected transitions span neutral, low-ionization, and high-ionization species, including O I, Si II, C II, Si IV, and C IV, with velocity offsets of $|\Delta v|\sim 50$–$250\,\mathrm{km\,s^{-1}}$ relative to nebular systemic redshifts. The ionic coexistence, overlapping velocity structure, and equivalent-width ratios are consistent with outflowing or otherwise kinematically disturbed galaxy-associated gas, implying rapid metal enrichment. These results show that key conditions for baryon cycling were established in at least a subset of luminous galaxies within the first several hundred million years of cosmic time, well before the completion of reionization.
Read the full abstract

All authors Yongda Zhu, Zhiyuan Ji, George D. Becker, Jiani Ding, Eiichi Egami, Xiaohui Fan, Xiangyu Jin, Weizhe Liu, Jianwei Lyu, Zheng Ma, Suprabhas Narisetty, George H. Rieke, Yunjing Wu, Minghao Yue, Junyu Zhang, Marcia J. Rieke

Earth and Planetary Astrophysics · Cross-listed from physics.space-ph

Properties and Radial Evolution of Solar Wind Turbulence Near Mercury's Orbit

We present a comprehensive statistical study of the radial evolution of solar wind turbulence near Mercury's orbit using long-term magnetic field measurements from the MESSENGER mission. Owing to Mercury's highly elliptical orbit and the spacecraft's repeated, extended residence in the upstream solar wind, the data set provides more than 17,000 hours of observations, enabling robust statistics across well-defined heliocentric distance intervals (0.31-0.47 au). We find that inertial-range spectral slopes remain close to -3/2 throughout Mercury's orbit, showing no significant radial evolution. Combined with low magnetic compressibility, this result indicates a stable, predominantly Alfvenic inertial-range cascade already established here. In contrast, kinetic-range spectral slopes exhibit clear radial evolution, becoming progressively shallower with increasing heliocentric distance, highlighting the greater sensitivity of kinetic-scale turbulence to heliocentric conditions. The ion-scale spectral break frequency decreases with distance in the spacecraft frame, while its normalized form increases relative to the local proton cyclotron frequency, demonstrating that the break is not tied to proton cyclotron frequency but reflects evolving local plasma conditions. Magnetic compressibility shows a similar frequency dependence at all distances, with a subtle radial enhancement of compressive fluctuations at kinetic scales. Autocorrelation analysis reveals strong anisotropy, with the correlation times of field-aligned magnetic fluctuations increasing with heliocentric distance, while those of perpendicular fluctuations remain shorter and nearly invariant. Together, these results demonstrate a clear scale-dependent radial evolution of solar wind turbulence near Mercury-s orbit, providing new constraints on the development of kinetic processes in the inner heliosphere.
Read the full abstract

All authors Xinmin Li, Chuanfei Dong, Lina Z. Hadid, Sae Aizawa, Chi Zhang, Hongyang Zhou, Liang Wang, Jiawei Gao, James A. Slavin

Solar and Stellar Astrophysics

Predicting the thermodynamics in the chromosphere from the translation of SDO data into the IRIS$^{2}$ inversion results using a visual transformer model

We present SDO2IRIS$^2$: a visual transformer model that translates a combination of images of the chromosphere and transition region (TR), observed by AIA, and a line-of-sight magnetogram, provided by HMI, into temperature, line-of-sight velocity (v$_{los}$), velocity of the turbulent motions (v$_{turb}$), and electron density (n$_{e}$) in the chromosphere. Using the thermodynamic variables obtained from the inversion of the chromospheric lines Mg II h&k, observed by IRIS, as the target of the model, and the intensity images in the chromosphere and TR, and the photospheric magnetogram as the input, the predicted T and n$_{e}$ show a strong correlation ($\approx 0.80$) for $\approx$80% of the test inverted data, a moderate-to-strong correlation ($\approx0.63$) for 70% of the v$_{turb}$ of the target test inverted data, while for the $v_{los}$, the correlation is weak. Therefore, the predicted values by SDO2IRIS$^2$ may be used as an estimation of the thermodynamics in the chromosphere, either as a stand-alone result or as complementary information to other chromospheric data observed simultaneously. The execution time employed by SDO2IRIS$^2$ to obtain the thermodynamic values in the chromosphere is of the order of a few minutes, being $\le10$ minutes when using a CPU, and $\le5$ minutes when using a GPU. SDO2IRIS$^2$ opens a new avenue for the use of SDO data thanks to the inversions provided by IRIS observables.
Read the full abstract

All authors Alberto Sainz Dalda, Vishal Upendran, Juno Kim, Kyuhyoun Cho, Paul S. Killam, Viggo Hansteen, Bart De Pontieu

Instrumentation and Methods for Astrophysics

Joint Estimation of Properties of the Lunar Subsurface and Galactic Foregrounds with LuSEE-Night

The Lunar Surface Electromagnetics Experiment (LuSEE-Night) is a joint NASA-DOE-ESA low-frequency radio telescope that will reach the lunar far side in 2027. The unknown dielectric properties of the subsurface at the LuSEE-Night landing site impose the most significant limitation for precision instrument calibration, as reflections from the lunar subsurface can change the primary beam at the 10-20% level. Simulations of these effects have provided insight and concern, showing that the lunar subsurface modeled as a lossy dielectric can absorb a large amount of the power of the sky signal. While this absorption may not strongly impact the signal-to-noise ratio in a sky-noise-dominated regime, it could complicate the beam pattern and make the signal more difficult to model and interpret. We have simulated the far-field properties of the LuSEE-Night beam for varying dielectric profiles of the lunar subsurface. We find that varying the properties of the lunar subsurface has the most significant impact around the antenna resonance, impacting its amplitude, position and width. Conversely, changing the properties of the foreground impacts the data across the band. We use a Bayesian inference pipeline to jointly estimate parameters of a galactic foreground model and dielectric properties of the lunar subsurface around the LuSEE-Night landing site and find that parameters of both the galaxy and subsurface properties can be estimated jointly. While the modeling is somewhat idealized, we believe that the results are largely robust owing to the fact that spectral variations for plausible subsurface and galaxy models have very different spectral signatures.
Read the full abstract

All authors Fatima Yousuf, Zack Li, Stuart D. Bale, David W. Barker, Jack Burns, Christian H. Bye, Hugo Camacho, Cristina-Maria Cordun, Johnny Dorigo Jones, Adam Fahs, Sonia Ghosh, Keith Goetz, Robert Grimm, Sven Herrmann, Joshua J. Hibbard, Oliver Jeong, Marc Klein-Wolt, Léon V.E. Koopmans, Joel Krajewski, Corentin Louis, Milan Maksimović, Ryan McLean, Raul A. Monsalve, Arnur Nigmetov, Paul O'Connor, Aaron Parsons, Michel Piat, Marc Pulupa, Rugved Pund, David Rapetti, Kaja M. Rotermund, Benjamin Saliwanchik, Anže Slosar, Graham Speedie, Nikolai Stefanov, David Sundkvist, Aritoki Suzuki, Harish K. Vedantham, Philippe Zarka

Cosmology and Nongalactic Astrophysics · Cross-listed from hep-ph

Neutron Portal and Dark Matter-Baryon Coincidence: from UV Completion to Phenomenology

We present a dynamical solution to the dark matter-baryon coincidence problem based on the neutron portal operator connecting the visible and dark sector asymmetries. This framework is motivated by the possibility that a strongly supercooled dark confinement phase transition accounts for the nano-Hz stochastic gravitational wave signal observed by pulsar timing arrays, while also generating the dark matter and baryon asymmetry in the Universe. We show that the GeV-scale mass of asymmetric dark matter can be naturally correlated with the (multi-)TeV scale cut-off for the neutron portal through its ultraviolet completion. The dark sector is governed by an approximate fixed point and confines once the heavy portal states are integrated out, dynamically generating a scale of $\mathcal{O} ({\rm GeV})$. We analyze both tree and loop-level ultraviolet completions and demonstrate how the resulting confinement scale is linked to the effective neutron portal scale. We also discuss cosmological constraints and experimental prospects in beam dump searches and colliders for probing the neutron portal.
Read the full abstract

All authors Sudhakantha Girmohanta, Yuichiro Nakai, Yoshihiro Shigekami, Zhihao Zhang

Wednesday 22 April

Astrophysics of Galaxies

Orbital evolution of highly eccentric bodies embedded in a ringed accretion disc

Various processes can induce long-lived overdense rings and arcs in protoplanetary and AGN accretion discs, such as the accumulation of gas at the outer edge of the dead zone, or the infall of material. Using the local approximation of dynamical friction, we investigate the orbital evolution of a low-mass highly-eccentric point-mass accretor (perturber) embedded in an isothermal disc hosting a density ring. We specifically consider the regime in which the eccentricity exceeds four times the disc aspect ratio. For prograde perturbers, orbits that cross the ring progressively circularize while their semi-major axes converge toward the ring radius. As a result, perturbers accumulate, forming a population ring superimposed on the gaseous ring. The ring therefore acts as a migration trap for these eccentric orbits. We also find that prograde orbits tangent to the ring, either at apocentre or pericentre, remain tangential throughout their evolution; perturbers confined to these trajectories experience the highest accretion rates. In contrast, retrograde perturbers always migrate inward. Once the semi-major axis becomes smaller than the ring radius, the eccentricity grows, but not enough for the orbit to intersect the ring again. We also discuss how feedback effects, such as jet launching and thermal torques, could modify the effective forces acting on the perturbers.
Read the full abstract

Astrophysics of Galaxies

Constraining Dark Matter Density Profiles in UFDs with Wide Binaries: Forecast for the Chinese Space Station Survey Telescope

The internal structure of dark matter halos on sub-galactic scales remains a key open question, particularly in the context of the core-cusp problem. Ultra-faint dwarf galaxies (UFDs), owing to their extreme dark matter dominance, provide a promising laboratory to probe these density profiles through stellar tracers. In this work, we assess the capability of the Chinese Space Station Telescope (CSST) to detect and characterize wide binary stars in the nearby UFD Segue 1, using mock observations. We generate mock binary populations based on our existing $N$-body simulations and incorporate realistic CSST observational conditions, including the expected deep-field limiting magnitude ($g \sim 27.5$ mag) and a photometric completeness of approximately $90\%$. The two-point correlation function (2PCF) of stellar pairs is used as a statistical tool to recover the binary fraction under these assumptions. We find that CSST can robustly detect wide binaries at the $3\sigma$ level for binary fractions as low as $f_b \gtrsim 0.01$, provided a stellar sample size of $N_{\mathrm{star}} \gtrsim 2300$. However, distinguishing between cusped and cored dark matter profiles is significantly more demanding, requiring $N_{\mathrm{star}} \gtrsim 6000$ and $f_b \gtrsim 0.1$ within $\sim 40\mathrm{kpc}$.
Read the full abstract

All authors Yixi Tao, Haijun Tian, Bin Yue, Jorge Peñarrubia

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

Ultra High Energy Cosmic Rays from the Local Void

Ultra high energy cosmic rays have been see coming from the direction of the local cosmic void. We use this fact to argue that at least some of these these cosmic rays are relatively light magnetic monopoles and that their relative fraction above 1020 eV can be found from full sky observations.
Read the full abstract

High Energy Astrophysical Phenomena

An Analytic Threshold for LESA-Driven Negative ELN Flux Directions in Core-Collapse Supernovae: Derivation and Population Census

In core-collapse supernovae (CCSNe), deleptonization normally favors $\nu_e$ over $\bar{\nu}_e$ emission. However, lepton-number emission self-sustained asymmetry (LESA) can make the energy-integrated emitted lepton-number flux negative along some directions. We derive a simple diagnostic for this transition and test it in 33 independent 3D CCSN simulations: 25 Princeton/Fornax models ($8.1$–$100\,M_\odot$) and 8 Garching models, including non-, slow-, and fast-rotating $15\,M_\odot$ cases. Of 23 non-black-hole-forming Princeton models, 22 cross the threshold, with median onset $t_c=225\,\mathrm{ms}$, IQR $162$–$264\,\mathrm{ms}$, and cross-model scatter $\mathrm{CV}=18.6\%$. Full-sky flux-sign searches show that the threshold identifies the anti-LESA-pole transition, distinguishing the global LESA-driven crossing from early localized turbulent crossings. The fast-rotating Garching $15\,M_\odot$ model, where rapid rotation suppresses the LESA dipole, is correctly classified as a non-crosser without using any rotation parameter. Both black-hole-forming Princeton models cross near $250\,\mathrm{ms}$ post-bounce and remain above threshold for $1807$ and $2463\,\mathrm{ms}$ before collapse. Thus, in the next nearby CCSN, the emitted $\bar{\nu}_e$ energy flux may exceed the $\nu_e$ flux along some lines of sight. Such directions may also correlate with sustained fast flavor instability, although testing this requires local phase-space distributions or dedicated linear stability analysis. The relevant quantity here is the energy-integrated emitted flux field, i.e. a luminosity difference per steradian, not a neutrino number flux.
Read the full abstract

Cosmology and Nongalactic Astrophysics

Cosmological Gravitational Waves from Ultralight Vector Dark Matter

We compute the abundance of cosmological gravitational waves produced during the evolution of an ultralight vector (spin-1) dark matter field. A homogeneous background vector field breaks spatial isotropy, requiring a Bianchi I geometry and inducing a mixing between the scalar, vector, and tensor perturbation sectors. We derive the perturbation equations in this background and show that, as a consequence of this mixing, scalar perturbations act as a source of tensor modes, generating a stochastic GW background. The production and cosmological evolution of these gravitational waves are implemented in class.VFDM, a modified version of CLASS, from which we obtain the present-day spectrum.
Read the full abstract

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

Multimessenger probes of Axions from Compact Objects

Astrophysics plays a pivotal role in the quest for axions and axion-like particles, offering guidance to experimental efforts and enabling the investigation of axion properties that cannot be probed otherwise. In this context, the extreme conditions in the interiors of compact stellar objects – such as core-collapse supernovae, neutron stars, and binary neutron star mergers – significantly enhance axion production, providing unparalleled sensitivity to extremely feeble couplings to Standard Model particles. In this context, the techniques of multimessenger astrophysics deepens the understanding of powerful transient events, maximizing the capabilities of current instruments to identify possible signatures of axion emission.
Read the full abstract

High Energy Astrophysical Phenomena

Archival Multiband Gravitational-Wave Signals from Massive Black Hole Binary Mergers

While massive black hole binaries (MBHBs) merge at gravitational-wave frequencies above the pulsar timing array (PTA) sensitivity band, we show that they leave orphaned low-frequency contributions in the PTA pulsar term. Due to the light-propagation time between each pulsar in the array and Earth, the pulsar term acts as a time-delayed probe of a chirping merger with a specific frequency response determined by the direction of origin and intrinsic properties of the MBHB. We provide a detailed consideration of how such a multiband signal would manifest in a full PTA, demonstrate an approach to stack these orphaned pulsar terms across the array, and discuss prospects for an archival, multiband search in conjunction with MBHB mergers observed in astrometric data or spaceborne interferometers like LISA.
Read the full abstract

All authors Alexander W. Criswell, Stephen R. Taylor, Kris Pardo, Alberto Sesana, David Izquierdo, Silvia Bonoli, Daniele Spinoso

High Energy Astrophysical Phenomena

Gravity Echoes from Supermassive Black Hole Binaries

Pulsar timing arrays record gravitational waves from supermassive black hole binaries at two spacetime points: an Earth term, measured when the wave passes the Earth, and a pulsar term, measured when the wave passed each pulsar at an earlier epoch. We show that a future $\mu$Hz-band detection of a nearby massive binary by a mission such as $\mu$Ares would turn PTA pulsar terms into targeted probes of binary evolution. In analogy with supernova light echoes, each pulsar term acts as a gravity echo: a dated snapshot of the binary at an earlier stage of its inspiral. Together, the $\mu$Hz Earth-term measurement and the nHz pulsar-term echoes provide a temporal baseline that neither detector could access alone. For a fiducial equal-mass binary with total mass $10^9\,M_\odot$ at 80~Mpc, we find a combined pulsar timing array echo signal-to-noise ratio of 33, with up to 24 pulsars individually resolving the signal among pulsars with 50-year baselines. The angular dependence of the single-pulsar echo sensitivity alone enables independent sky localization of the source to $\sim$10–100~deg$^2$, and the resolved pulsar-term frequencies directly measure the binary inspiral rate hundreds to thousands of years ago. With sufficient pulsar distance precision, a small set of anchor pulsars could additionally phase-connect the array and trace the post-Newtonian evolution coherently over kpc baselines. The source population required for gravity echoes is drawn from the same massive-end census responsible for the observed nanoHertz stochastic background.
Read the full abstract

Astrophysics of Galaxies

DeepDive: Simultaneous Formation of Massive Quiescent Galaxies in High-Redshift Galaxy Overdensities

We report on the spectroscopic confirmation of overdense regions of massive quiescent galaxies (QGs) in the early Universe with JWST/NIRSpec. Based on data from the DeepDive NIRSpec program and archival data from the Dawn JWST Archive, we confirm three QGs in the vicinity of Jekyll & Hyde, a pair of massive QG and a dusty star-forming galaxy, at $z=3.71$ and two QGs around SXDS-27434 at $z=4.01$. According to the analysis of galaxy number density with photometric redshifts, Jekyll & Hyde (SXDS-27434) are in an overdense region, where the number density of galaxies is three (four) times higher than the average in the COSMOS (SXDS) field. SED fitting suggests that most of the QGs follow similar star formation histories and have consistent formation and quenching epochs. The same trend is observed in other proto-clusters hosting QGs that were already identified by ground-based telescopes, indicating that the large-scale environment plays an important role in the formation of QGs. In addition, JWST spectra reveal a broad H$\alpha$ emission line from SXDS-27434 and faint emission lines from other three QGs, which are identified as AGN-driven based on their emission line ratios. The overdensity is also reproduced by the Illustris TNG300 simulation at $z=3.71$, in which the member QGs also have similar quenching epochs. These results are consistent with a scenario in which the large-scale structure enhances merger activity and/or gas accretion and triggers AGN feedback, thereby driving simultaneous formation and quenching of QGs in the overdensity.
Read the full abstract

All authors Takumi Kakimoto, Masayuki Tanaka, Kei Ito, Francesco Valentino, Makoto Ando, Gabriel Brammer, Massissilia L. Hamadouche, Vasily Kokorev, Jacqueline Antwi-Danso, William M. Baker, Daniel Ceverino, Andreas L. Faisst, Marion Farcy, Michaela Hirschmann, Christian Kragh Jespersen, Mariko Kubo, Allison W. S. Man, Masato Onodera, Rhythm Shimakawa, Sune Toft, John R. Weaver, Po-Feng Wu, Pengpei Zhu

Astrophysics of Galaxies · Cross-listed from physics.chem-ph

Chaos Gated Tunneling Drives Molecular Reactivity in Astrophysical Environments

Accurate modeling of ion-molecule reaction networks is essential for understanding the chemical evolution of planetary ionospheres, particularly for giant planets where proton-transfer chains drive atmospheric composition. However, predicting reaction rates in these ultracold environments remains a challenge due to the non-trivial interplay between vibrational dynamics and quantum tunneling. In this work we present a chaos-diagnostic framework that integrates multireference electronic structure theory, Adiabatic Gauge Potentials (AGP), and Random Matrix Theory (RMT) to characterize the microscopic dynamics of proton transport. Using the formation of H+3 and the proton-bound cluster H+5 as representative model systems relevant to Jovian atmospheres, we demonstrate that the transition state acts as a dynamical bottleneck where quantum chaos is notably suppressed, effectively enhancing tunneling probabilities. We introduce a fragility index based on the AGP slope to quantify how specific vibrational modes reintroduce chaos and suppress reactivity. This diagnostic approach offers a generalizable, data-driven metric for identifying vibrationally gated pathways in complex astrochemical networks, providing a theoretical basis for refining kinetic models of planetary and interstellar plasmas
Read the full abstract

All authors Saptarshi G. Dastider, K. Prashant, P. Shruti, C. Sudheesh, Jobin Cyriac

Cosmology and Nongalactic Astrophysics

Relativistic effects in k-essence

Relativistic effects are sensitive to subtle changes in dark energy. These effects grow on very large scales and at high redshifts, which will be the reach of upcoming surveys. We investigate these effects in both the linear and the angular galaxy power spectra in a late-time universe dominated by cold dark matter and k-essence, focusing on three core models (dilaton, tachyon, and DBI scalar fields) and contrasting their predictions with those of the concordance model. By enforcing identical present-day cosmological parameters, we isolate the imprints of k-essence dynamics and perturbations on very large scales. We found that relativistic corrections dominate on very large scales and grow with redshift, but are largely insensitive to k-essence microphysics in Fourier space, leading to strong degeneracies among the models. However, in the angular power spectrum, where line-of-sight integrals are naturally included, relativistic effects are significantly amplified, yielding better sensitivity to clustering k-essence. In particular, the tachyon exhibits clear deviations across multipoles and redshifts, with distinct imprints in the Doppler and the combined (velocity and gravitational) potentials contributions. Furthermore, our results show that neglecting relativistic corrections can lead to systematic misestimation of deviations of k-essence from the cosmological constant. Our results show the relativistic angular galaxy power spectrum as a more consistent and robust probe of ultra-large-scale physics. These findings underscore the need for full relativistic modelling in next-generation surveys that are targeting horizon-scale modes, where the imprint of non-standard dark energy is most pronounced.
Read the full abstract

All authors Didam Duniya, Isaac Opio, Bishop Mongwane, Hassan Abdalla

High Energy Astrophysical Phenomena

Probing past mergers of supermassive black holes with pulsar timing arrays: The role of pulsar terms

By monitoring the times of arrival of radio pulses from millisecond pulsars, Pulsar Timing Arrays (PTAs) serve as unique gravitational wave (GW) laboratories in the nanohertz band. To date, the primary astrophysical sources of GWs targeted in this frequency range have been inspiraling supermassive black hole binaries (SMBHBs) on circular and eccentric orbits. In this work, we demonstrate that, thanks to the so-called pulsar term in the timing residual waveform of GW signals, PTAs can probe individual SMBHBs that merged before timing observations began. We refer to the latter as zombie binaries. Using SMBHB population models consistent with current PTA constraints, we find that while the probability of detecting such systems in existing PTA datasets remains low, the Square Kilometer Array observatory is expected to achieve sufficient sensitivity to have a few zombie binaries with optimal matched filter signal-to-noise ratios exceeding 3 in its data. Although their confident identification might be challenging, this new class of PTA sources opens a novel window for studying the most massive SMBHBs in our local universe.
Read the full abstract

Astrophysics of Galaxies

Chaotic migration of LISA Extreme Mass Ratio Inspirals in a turbulent accretion disk: effect on waveform de-phasing

Gravitational wave (GW) detector LISA will observe near-coalescence extreme mass ratio inspirals (EMRIs), which typically form in galactic central accretion disks. Torques from the disk can alter the GW-driven inspiral trajectory of an embedded EMRI from the vacuum expectation, leading to potentially observable GW dephasing ($\Delta\psi_{\rm gas}$). So far, all studies compute $\Delta\psi_{\rm gas}$ for a thin, laminar disk, with negligible flow turbulence, whereby the disk exerts the well-understood linear torque ($T_{\rm lin}$). However, these disks must be turbulent due to magneto-rotational instability in the inner regions. Hence, we present a proof-of-concept general prescription for the turbulent torque ($T_{\rm turb}$) acting on an EMRI by modeling it as a Gaussian distribution around $T_{\rm lin}$, inspired by recent global simulations that study such torques. We compute $\Delta\psi_{\rm gas}$ for the “golden” circular EMRI with total source mass $M=10^6~{\rm M}_\odot$ and mass ratio $q=5\times10^{-5}$ in its final four-year evolution at redshift $z=0.276$ and signal-to-noise ratio (SNR) $=50$ by varying turbulence amplitude $C$ ($=1$ in the aforementioned study), maximum correlation timescale ($N_{\rm max}$), Eddington ratio ${\rm f}_{\rm Edd}$, disk aspect ratio $h_0$, and turbo-viscous coefficient $\alpha$ in a reasonable parameters space. For $N_{\rm max}=100$ orbits, we find that for $C\gtrsim{10}$, ${\rm f}_{\rm Edd}\gtrsim0.3$, $h_0\gtrsim0.03$, and $\alpha\gtrsim0.1$, dephasings due to $T_{\rm lin}$ are unobservable but could become detectable ($\Delta\psi_{\rm gas}>8/$SNR) if EMRIs experience turbulent torques. Hence, this work motivates running MHD simulations of accretion disks with embedded early-inspiral LISA EMRIs over long timescales to understand the imprint of the turbulent environment on their orbital parameters and gravitational waveforms.
Read the full abstract

All authors Mudit Garg, Lucio Mayer, Yinhao Wu, Yacine Ali-Haïmoud, Douglas N.C. Lin

Astrophysics of Galaxies

Light, heavy, primordial: exploring the diversity of black hole seeding and growth mechanisms in the JWST era

The James Webb Space Telescope (JWST) has revealed a puzzling population of massive black holes in the first billion years, many of which are over-massive compared to their hosts (obese black holes), and reside in metal-poor hosts, posing a challenge for theoretical models at these early epochs. In this work, we compare the observational properties of astrophysically-seeded black holes using the DELPHI semi-analytic model and cosmologically-seeded primordial black holes (PBHs) using the PHANES analytic model. We explore the growth of light ($\sim 100 M_\odot$) and heavy ($\sim 10^{3-5}M_\odot$) seeds through mergers and accretion (both Eddington-limited and at super-Eddington rates) in the astrophysical scenario; PBHs (seeded between $10^{0.5-6}M_\odot$) only grow through accretion at sub-Eddington rates. Comparing to observables at $z \sim 5-10$, the only model that can be ruled out is the one where we allow Eddington-limited accretion onto light seeds. The observed high values of the black hole mass-stellar mass relation ($0.3-1$) can be reproduced by both PBHs and heavy seeds accreting at super-Eddington rates. However, only the PBH and Eddington-limited heavy seeding models can simultaneously reproduce the observed black hole masses (${\rm M_{bh}}$), stellar masses ($M_*$), and extremely low host metallicities ($Z \leq 0.01 Z_\odot$) inferred at $z \sim 7-10$. Crucially, we find PBHs show decrease in the black hole mass-stellar mass ratio with increasing halo mass at all redshifts, contrary to any astrophysical black hole model. Selecting systems at $z \sim 7$ with ${\rm M_{bh}}/M_* > 0.1$ and bolometric luminosities $\sim 10^{44-46} {\rm erg~s^{-1}}$ that show a negative black hole to stellar mass ratio and reside in $10^{9-11}M_\odot$ halos offer a promising clustering-based discriminant of PBH seeding models.
Read the full abstract

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

Turbulence Mode Decomposition and Anisotropy in Magnetically Dominated Collisionless Plasmas

We use the 3D fully kinetic simulation to study different turbulence modes and turbulence anisotropy of relativistic turbulence in magnetically dominated collisionless plasmas. We extend the method developed by Cho & Lazarian (2002) for decomposing non-relativistic magnetohydrodynamic (MHD) turbulence into Alfvén, fast, and slow modes to the regime of collisionless plasmas. We find that Alfvén and slow modes are anisotropic, following the Goldreich & Sridhar (1995) scaling, while fast modes are isotropic. We observe a larger kinetic energy fraction of fast modes compared to that in the non-relativistic MHD turbulence, suggesting a stronger coupling of Alfvén and fast modes in relativistic magnetized turbulence in collisionless plasmas. We further examine the dynamic alignment and find a weaker scale dependence of the alignment angle than previously proposed. The dominant thermal fluctuations in the kinetic range can cause flattening of the turbulent velocity structure function and weakening of the turbulence anisotropy and dynamic alignment near the kinetic scales.
Read the full abstract

All authors Samuel T. Sebastian, Siyao Xu, Yue Hu, Luca Comisso, Saikat Das, Joonas Nättilä

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

Radial adiabatic perturbations of stellar compact objects

We present a covariant and gauge-invariant formulation of the theory of radial adiabatic linear perturbations of self-gravitating, non-dissipative imperfect fluids within the theory of general relativity. By codifying the thermodynamical properties of the source into an equation of state and an ansatz on anisotropic pressure that involves both matter and kinematic variables, we obtain a set of equations that is directly applicable to a wide variety of thermodynamic theories for matter fields. As examples, we evaluate and compare the predictions of the Eckart theory, the Bemfica-Disconzi-Noronha-Kovtun theory, and the Truncated Israel-Stewart theory on the properties and evolution of radial adiabatic perturbations of stellar compact objects modeled by classical equilibrium solutions. Introducing a new solution of the Einstein field equations, and imposing causality, we propose an upper bound for the maximum compactness of dynamically stable stars with non-trivial radial and tangential pressures.
Read the full abstract

Astrophysics of Galaxies

Characterizing the GD-1 Stream with DESI DR2 Data: Thin Stream and Hot Cocoon

GD-1 is among the longest, coldest stellar streams in the Milky Way, making it an ideal target for probing dark matter substructure through dynamical heating. We present a catalog of 608 spectroscopically confirmed GD-1 members from the first three years of Dark Energy Spectroscopic Instrument (DESI) observations. This constitutes the largest homogeneous spectroscopic sample of GD-1, doubling the number of members previously available only through heterogeneous compilations combining multiple surveys with different systematics. Using these data, we derive updated stream tracks in sky position, proper motion, and radial velocity that extend over $100^\circ$ of the stream. We apply a Gaussian mixture model to decompose the stream into a dynamically cold thin component ($\sigma_V = 2.49\pm 0.28$ km s$^{-1}$, width $= 0.23\pm0.01^\circ$) and a kinematically hot cocoon ($\sigma_V = 6.13\pm0.75$ km s$^{-1}$, width $= 2.18\pm0.17^\circ$). The cocoon contains $\sim30\%$ of members and its velocity dispersion is consistent with $\sim11$ Gyr of heating by cold dark matter subhalos. We also detect a large proper motion dispersion ($41.36\pm4.98$ km s$^{-1}$) along the stream direction in the cocoon component. This feature indicates a significant line-of-sight distance spread in the cocoon, and its origin will be further explored in a forthcoming paper. These measurements demonstrate the power of DESI spectroscopy for characterizing the multi-component phase-space structure of stellar streams and constraining small-scale dark matter substructure.
Read the full abstract

All authors Emma Jarvis, Ting S. Li, Sergey E. Koposov, Raymond G. Carlberg, Monica Valluri, Nasser Mohammed, J. Aguilar, S. Ahlen, Carlos Allende Prieto, Leandro Beraldo e Silva, D. Bianchi, D. Brooks, Amanda Byström, T. Claybaugh, A. P. Cooper, A. Cuceu, A. de la Macorra, Arjun Dey, Biprateep Dey, P. Doel, J. E. Forero-Romero, E. Gaztañaga, Oleg Y. Gnedin, Satya Gontcho A Gontcho, G. Gutierrez, K. Honscheid, R. Joyce, R. Kehoe, T. Kisner, Namitha Kizhuprakkat, A. Kremin, Mika Lambert, M. Landriau, L. Le Guillou, Gustavo E. Medina, A. Meisner, R. Miquel, S. Nadathur, Joan Najita, N. Palanque-Delabrouille, W. J. Percival, F. Prada, I. Pérez-Ràfols, Tian Qiu, Alexander H. Riley, Constance M. Rockosi, G. Rossi, E. Sanchez, Nathan Sandford, E. F. Schlafly, D. Schlegel, J. Silber, D. Sprayberry, G. Tarlé, B. A. Weaver, R. Zhou, H. Zou

Astrophysics of Galaxies

Black hole mass, host galaxy mass, and dark matter halos: Testing the environmental connection

We investigate the connection between supermassive black holes (SMBHs), their host galaxies, and large-scale dark-matter halos using broad-line X-ray AGN from the XMM–XXL and Stripe\,82X surveys, together with galaxies from VIPERS and SDSS/Stripe\,82. Building on the homogeneous host-galaxy catalogue presented in Paper~I, we test whether AGN with a given black-hole mass, $M_{\rm BH}$, inhabit different large-scale environments from non-AGN galaxies with similar host properties. We first examine the empirical $M_{\rm BH}$–$M_{\star}$ relation of the AGN sample. We find a shallow trend with substantial scatter, likely driven by flux-limited selection effects and uncertainties in virial black-hole mass estimates. The ratio $M_{\rm BH}/M_{\star}$ decreases with increasing stellar mass, and AGN lying above and below the empirical relation show different median host properties, consistent with non-synchronous SMBH and stellar growth. We then divide the AGN into two black-hole mass bins, $8.0 \le \log(M_{\rm BH}/M_\odot) < 8.5$ and $8.5 \le \log(M_{\rm BH}/M_\odot) < 9.0$, and construct galaxy control samples matched in $M_{\star}$, SFR, and sSFR using a multivariate nearest-neighbour method. From AGN–galaxy cross-correlation functions, we infer the characteristic halo masses of AGN and matched galaxies. In the lower-$M_{\rm BH}$ bin, AGN occupy halos statistically indistinguishable from those of their controls. In the higher-$M_{\rm BH}$ bin, we find a mild indication that AGN may reside in somewhat more massive halos, with a difference of about 0.4 dex, although still consistent within the uncertainties. If confirmed with larger samples, this would suggest that halo-scale processes become important mainly at the highest $M_{\rm BH}$.
Read the full abstract

All authors G. Mountrichas, F. Shankar, F. J. Carrera, A. Georgakakis

Astrophysics of Galaxies

Revisiting radio synchrotron diagnostics in star-forming galaxies

Radio continuum observations are widely used to study cosmic ray (CR) electron populations and transport processes in star-forming galaxies, but their interpretation relies on several simplifying assumptions. Here, we revisit three common assumptions: that some vertical radio profiles can be explained by CR advection alone, that radio spectra directly trace the galaxy-wide CR electron spectrum, and that bremsstrahlung and Coulomb losses are negligible for radio-emitting electrons. We model radio emission using time-dependent CR electron evolution in a magnetohydrodynamical simulation of an isolated Milky Way-mass galaxy. CR electron spectra are evolved self-consistently along Lagrangian tracer particles with the CREST framework, including injection at supernova remnants, advection with the gas, and spatially and temporally varying radiative losses. We compare these results to commonly adopted steady-state models. We find that advection-only transport in self-consistently driven galactic winds fails to reproduce the extended vertical radio intensity profiles observed in edge-on galaxies, despite reproducing the observed steepening of spectral indices with height. This is because slowly accelerating winds keep electrons in strong cooling environments for too long. Matching observed radio haloes with advection alone requires unrealistically high midplane wind velocities, implying that additional transport or re-acceleration processes are required. Although galaxy-integrated CR electron spectra at radio-emitting energies are similar across models, the resulting synchrotron spectra differ systematically because radio emission is biased toward young electrons in dense, strongly magnetised regions. Finally, we show that bremsstrahlung and Coulomb losses significantly shape radio spectra even when their loss rate is subdominant and therefore cannot be neglected.
Read the full abstract

All authors Maria Werhahn, Christoph Pfrommer, Philipp Girichidis, Joseph Whittingham, Léna Jlassi, Rüdiger Pakmor, Rebekka Bieri, Rainer Weinberger, Volker Springel, Freeke van de Voort

Cosmology and Nongalactic Astrophysics

Interpretable Analytic Formulae for GWTC-4 Binary Black Hole Population Properties via Symbolic Regression

Recent LIGO-Virgo-KAGRA (LVK) analyses have revealed complex structure in the binary black hole (BBH) population, including distinct features in the primary mass spectrum and nontrivial spin-mass correlations. However, the phenomenological models used to capture these features often lack analytic transparency, making it difficult to isolate robust physical laws from modeling artifacts. To address this, symbolic regression is applied to the posterior inference products of the GWTC-4 catalog, discovering compact, closed-form analytic expressions for four key population relationships: (i) the merger-rate evolution with redshift; (ii) the mass-ratio dependence of the effective-spin distribution; (iii) the redshift evolution of the effective-spin distribution; and (iv) the conditional mass-ratio distributions associated with the 10 solar mass and 35 solar mass primary mass peaks. This framework successfully compresses both rigid and highly flexible models into differentiable phenomenological laws, dynamically recovering a consistent low-redshift merger-rate slope without assuming an a priori power-law form. The exact analytic derivatives provided by symbolic regression show that the mass ratio–effective spin and redshift–effective spin correlations are robustly driven by broadening of the posterior widths rather than shifts in the mean. Furthermore, qualitatively distinct functional forms for the mass-ratio distributions conditioned on the 10 solar mass and 35 solar mass primary mass peaks are identified. These closed-form expressions enable exact analytic gradient diagnostics and compact surrogate summaries, particularly for flexible numerical posteriors that are not otherwise available in low-dimensional analytic form. They also facilitate rapid downstream calculations for rate forecasting, formation channel comparison, and stochastic background estimation.
Read the full abstract

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

Cracking Gravitational Wave Multiple Ringdown Modes in Space

Ringdown signals from perturbed black holes (BHs) offer a clean window into BH spacetime, strong-field gravity, and fundamental physics. Presently the quasi-normal modes of stellar-mass BH ringdowns have been successfully extracted in the ground-based gravitational wave (GW) observations. Looking ahead, the future space-borne observatories will listen to the ringdowns from massive BH binary coalescences more loudly and resolve multiple modes to unprecedented precision, which calls for efficient approaches to mitigate the sharply increasing computational burden. We develop a practical ringdown analysis pipeline for space-borne detectors by implementing FIREFLY, a novel acceleration algorithm validated in ground-based detectors, and for the first time demonstrate its compatibility and effectiveness with the time-delay interferometry (TDI) observables. With high fidelity, we achieve an over two-order-of-magnitude speedup for a simulated ringdown signal including six modes, providing a viable and scalable route for multi-mode ringdown analysis in the space context. This new approach has sound statistical interpretation and is extensible to other GW sources in band.
Read the full abstract

All authors Ziming Wang, Han Wang, Yuxin Yang, Yiming Dong, Hai-Tian Wang, Yi-Ming Hu, Lijing Shao

Instrumentation and Methods for Astrophysics

Planetary Exploration 3.0: A Roadmap for Software-Defined, Radically Adaptive Space Systems

The surface and subsurface of worlds beyond Mars remain largely unexplored. Yet these worlds hold keys to fundamental questions in planetary science - from potentially habitable subsurface oceans on icy moons to ancient records preserved in Kuiper Belt objects. NASA's success in Mars exploration was achieved through incrementalism: 22 progressively sophisticated missions over decades. This paradigm, which we call Planetary Exploration 2.0 (PE 2.0), is untenable for the outer Solar System, where cruise times of a decade or more make iterative missions infeasible. We propose Planetary Exploration 3.0 (PE 3.0): a paradigm in which unvisited worlds are explored by a single or a few missions with radically adaptive space systems. A PE 3.0 mission conducts both initial exploratory science and follow-on hypothesis-driven science based on its own in situ data returns, evolving spacecraft capabilities to work resiliently in previously unseen environments. The key enabler of PE 3.0 is software-defined space systems (SDSSs) - systems that can adapt their functions at all levels through software updates. This paper presents findings from a Keck Institute for Space Studies (KISS) workshop on PE 3.0, covering: (1) PE 3.0 systems engineering including science definition, architecture, design methods, and verification & validation; (2) software-defined space system technologies including reconfigurable hardware, multi-functionality, and modularity; (3) onboard intelligence including autonomous science, navigation, controls, and embodied AI; and (4) three PE 3.0 mission concepts: a Neptune/Triton smart flyby, an ocean world explorer, and an Oort cloud reconnaissance mission.
Read the full abstract

All authors Masahiro Ono, Daniel Selva, Morgan L. Cable, Marie Ethvignot, Margaret Hansen, Andreas M. Hein, Elena-Sorina Lupu, Zachary Manchester, David Murrow, Chad Pozarycki, Pascal Spino, Amanda Stockton, Mathieu Choukroun, Soon-Jo Chung, John Day, Alexander Demagall, Anthony Freeman, Chloe Gentgen, Michel D. Ingham, Charity M. Phillips-Lander, Richard Rieber, Alejandro Salado, Maria Sakovsky, Lori R. Shiraishi, Yisong Yue, Kris Zacny

Astrophysics of Galaxies

Self-regulated galaxy evolution within a self-consistently varying galaxy-wide IMF

Semi-analytical evolution models of galaxies are a useful and computationally inexpensive tool for fast assessment of individual properties and their evolution. In this work, specifically the influence of a metallicity and star-formation rate (SFR) dependent galaxy-wide stellar initial mass function (IGIMF) on the self-regulation of star-formation in a galaxy is of interest. All models – both non-varying gwIMFs and the IGIMF – reproduce reasonable gas fractions, gas depletion timescales and the main sequence of star-forming galaxies. However, only the IGIMF model accurately predicts the mass-metallicity relation and provides a more comprehensive description of quenched elliptical galaxies. For massive ellipticals all models suggest the need for an additional gas heating source to reach a quenched state. Using a different stellar yield table in the IGIMF model does not significantly affect the results. In all models, the galaxies evolve self-regulated, determined by the accretion rate. The self-regulated constancy of the SFR reflects the constant SFRs of nearby star-forming galaxies. The specific gas-accretion rate of all galaxies appears to be comparable to the Hubble constant. The inclusion of outflows improves the results for the canonical gwIMF model, but not significantly, while for the IGIMF model it has no significant impact.
Read the full abstract

All authors Lukas Hof, Pavel Kroupa, Gerhard Hensler, Jan Pflamm-Altenburg