Peer Revenue Publish or profit

Papers

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

Physics Astrophysics

astro-ph Watch Astrophysics (sign in to use)

Saturday 2 May

Astrophysics of Galaxies

Cross-Comparison of Galaxies Detected in the CSST Spectroscopic Survey and the SKA HI Survey

We present a forward-modeling framework to forecast the galaxies detected in the Chinese Space Station Survey Telescope (CSST) spectroscopic survey and the Square Kilometre Array (SKA) HI survey. Starting from the L-Galaxies 2020 semi-analytic model run on the Millennium-II N-body simulation (MS-II), the cold gas in galaxies is partitioned into atomic and molecular components self-consistently within the model. We further model the emission-lines (H $\alpha$, H $\beta$, O III) relevant for the slitless spectrograph of the CSST in a post-processing step. We construct mock lightcones using the Mock Map Facility (MoMaF) approach, simulating the neutral hydrogen (HI) data cubes representing a 2000 hour SKA-Mid spectral line observation from redshifts 0.25–0.5, and employ the Source Finding Application 2(SOFIA-2) source-finding package to generate an HI galaxy catalog. In parallel, we apply the CSST selection function and noise model to obtain a realistic catalog of emission-line galaxies; the emission-line signal is proportional to the star formation rate. These products allow us to cross compare the galaxy samples and assess the synergy between CSST and SKA. We study the correlations of the HI and the emission-line signal with the halo mass, HI mass, and the stellar mass, and the baryonic Tully-Fisher relation (BTFR). We also perform stacking analysis of the HI signal from the CSST-selected sample, which probes the HI content in galaxies with low HI mass. Finally, we derive the optical-HI cross-correlation power spectrum of the galaxies, and measure the bias of these galaxies. These results can provide useful insight on the cold gas and stellar content of the galaxies.
Read the full abstract

All authors Yingfeng Liu, Furen Deng, Wenxiang Pei, Haitao Miao, Qi Xiong, Shuanghao Shu, Xingchen Zhou, Qi Guo, Yan Gong, Yougang Wang, Xuelei Chen

High Energy Astrophysical Phenomena

Type Ib supernovae are bluer than Type Ic supernovae

Type Ib and Ic supernovae (SNe Ib/Ic) are the bright finale of massive stars that have lost their hydrogen envelopes, making them powerful probes of mass stripping in massive star evolution. The advent of modern large photometric and spectroscopic surveys presents a unique opportunity to investigate systematic differences between these two kinds of SNe. In this study, we analyze a large, homogeneous sample of SNe Ib/Ic light curves from the Zwicky Transient Facility. We find a systematic difference in their apparent optical colors at peak: SNe Ib are, on average, bluer than SNe Ic, with two-sample tests confirming that the distributions differ (p<0.05). The difference in their host galaxy reddening, as currently constrained, is too small to fully explain it. The color difference therefore most likely has an intrinsic origin, reflecting progenitors with different degrees of stripping – helium-rich for SNe Ib and helium-poor for SNe Ic. In addition, we find that SNe Ib/Ic with narrow lines (SNe Ibn/Icn) are bluer than those without, which might originate from circumstellar matter interaction, with potential connection to fast blue optical transients. We demonstrate that SN colors offer a promising probe of mass stripping in massive stars, potentially providing a useful tool for analyzing large photometric data and improving predictions for the final outcomes of stripped massive stars.
Read the full abstract

All authors Harim Jin, Selma E. de Mink, Sebastian Holzner, Jakub Klencki, Géza Csörnyei, Sung-Chul Yoon, Iair Arcavi, Wolfgang E. Kerzendorf

Instrumentation and Methods for Astrophysics

Statistically Significant Linear Alignments Among High-Confidence Transient Candidates on POSS-I Photographic Plates

I report the detection of statistically significant linear alignments and anomalous spatial clustering among high-confidence transient candidates in the VASCO catalog of vanishing sources on Palomar Observatory Sky Survey (POSS-I) photographic plates (1949-1957). A machine learning classifier scores 107,875 candidates by their likelihood of being genuine transients. Searching the 36,215 candidates with probability >= 0.50 for collinear groupings narrower than 3 arcsec, I find 7 plates with alignments of 5-8 sources that exceed Monte Carlo expectations (p < 0.03, 10,000 iterations). The aligned sources are point-like, not streaks, which rules out any continuously luminous object crossing the field during the 45-minute exposures. The implied angular rates (1-15 arcsec/s) overlap with the geosynchronous regime but are inconsistent with low or medium Earth orbits, and no artificial satellites existed during the POSS-I era. When I project each alignment onto Earth's surface assuming a high-altitude object, 6 of 7 maintain constant geographic longitude with sub-degree spread (combined p ~ 3e-10). Four of these cluster near -96 deg longitude (central United States); one falls within 0.3 deg of the longitude of the Hanford nuclear production site on a nuclear test window date. Close pairs (< 30 arcsec) occur at 16.2x the random rate, and the nights with alignments are the same nights with excess close pairs (Fisher exact p < 0.0001). Plate artifacts cluster near the ecliptic plane (26%), but high-confidence transients are depleted there (16%; chi-square test p = 3.3e-82), which rules out asteroids, comets, and zodiacal debris as the dominant source. No transient reappears at the same sky position on a different night. All of these transients predate Sputnik 1.
Read the full abstract

Friday 1 May

High Energy Astrophysical Phenomena

XRISM/Resolve observations of Hercules X-1: a pulsating, highly broadened Fe K emission line from the neutron star accretion column

The study of X-ray pulsar accretion columns helps us characterize accretion physics in this extreme regime of strong gravity and strong magnetic fields. Previous observations of the X-ray pulsar Hercules X-1 revealed a highly broadened Fe K emission line, associated with Doppler motions exceeding 0.1c, suggesting its origin in the accretion column. We obtained a high-spectral resolution view of the Fe K energy band of Hercules X-1 thanks to a 200 ks observation with the XRISM observatory. The XRISM/Resolve microcalorimeter spectra allow us to separate the different spectral components and accurately model them with phenomenological models. We confirm the presence of a broad line near 6.5 keV with a typical $1\sigma$ width of 1 keV. Performing a pulse-phase-resolved analysis, we find that the feature is strongly variable with Her X-1 pulse phase. This is consistent with the proposed origin due to collisional recombination or by reprocessing of the primary X-ray emission in the accretion column, where strong variability with pulse phase is expected due to the rotation of the columns alongside with the neutron star. Additionally, the Fe K line pulsation pattern evolves with the 35-day cycle of Hercules X-1, supporting the scenario that the neutron star and its accretion columns undergo precession, in agreement with recent polarimetric results from the IXPE observatory. We discuss the future applications of modeling of this broad line in X-ray pulsars with physical spectral models. This could be used to detect and track neutron star precession, advancing our understanding of neutron star interiors.
Read the full abstract

All authors Peter Kosec, Laura Brenneman, Erin Kara, Ciro Pinto, Daniele Rogantini, Rudiger Staubert, Dominic Walton, Francesco Barra, Andrew Fabian, Teruaki Enoto, Jon M. Miller, Takuto Narita, Koh Sakamoto, Yutaro Nagai

High Energy Astrophysical Phenomena

Detectability of Polarized Gamma-ray Emission from Blazar Flares with COSI

We investigate the detectability of polarized gamma-ray emission from blazar flares with the Compton Spectrometer and Imager (COSI). Using 17 years of Fermi Large Area Telescope observations, we analyze light curves for 1413 blazars and identify a maximum of 787 sources with flaring episodes through Bayesian block analysis. For each flare, we estimate the minimum detectable polarization MDP99 in the COSI energy band (0.2-5 MeV) using instrument response functions under a range of spectral assumptions and background conditions. Under baseline background levels (1 counts/s), and assuming that blazar flare statistics in the MeV band are comparable to those observed at GeV energies, we find that COSI can realistically detect polarization in up to ~6 flares with MDP99<50% over its two-year prime mission depending on different spectral and flare identification assumptions, with only a few most powerful ones reaching MDP99<20%. These expectations are shown to improve when shorter intervals around bright peaks within long flares are considered. We provide a ranked list of the most promising targets, finding that flat-spectrum radio quasars dominate the population of polarization-detectable events. Through its continuous all-sky monitoring in the largely unexplored MeV band, COSI will open a new observational window on blazar variability and deliver the first direct measurements of MeV polarization, offering unique insights into jet geometry and high-energy emission processes.
Read the full abstract

All authors Garrett A. Latiolais, Jorge Otero-Santos, Michela Negro, Lea Marcotulli, Mohammad Ali Boroumand, Savitri Gallego, Christopher M. Karwin, Israel Martinez-Castellanos, Daniel Kocevski, Marco Ajello, Sara Capecchiacci, Ioannis Liodakis, Srinadh R. Bhavanam, Steven E. Boggs, Dieter H. Hartmann, Carolyn A. Kierans, Tiffany R. Lewis, Alberto Sciaccaluga, John A. Tomsick, Haocheng Zhang, Andreas Zoglauer

Astrophysics of Galaxies

Resolving the Unresolved Galactic Winds in Multi-phase Models. I. Methodology and Application

Galactic winds shape galaxy evolution; however, the outflowing gas is complex: it consists of multiple ionization phases, and its properties vary spatially. Therefore, methods that combine high-fidelity observations with state-of-the-art galactic-wind models are limited. Here we investigate methods for fitting the column density profiles derived from high-quality outflow observations with the multiphase, multiscale wind model from Fielding & Bryan 2022. We identify three key outflow parameters: the initial hot-phase mass-loading factor ($\eta_\text{ M,hot,0}$), the initial cool-phase mass-loading factor ($\eta_\text{ M,cool,0}$), and the initial cool-cloud mass. We obtain good fits for most galaxies, with tight constraints on $\eta_\text{ M,cool,0}$ and moderate constraints on the other two parameters. We find the inferred $\eta_\text{ M,cool,0}$ and $\eta_\text{ M,hot,0}$ are mostly of order unity, with significant scatter. The constraints on $\eta_\text{ M,hot,0}$ suggest that the interaction between the cool and hot phases allows us to indirectly constrain the properties of the hot wind from cool-outflow observations. The model also predicts various radial trends. First, for all galaxies, the cool-phase outflow velocity increases between $1-2$ times of the half-light radius, then reaches a plateau. Second, most galaxies exhibit increasing $\eta_\text{ M,cool}$ and decreasing $\eta_\text{ M,hot}$ with radius, with a few showing the reverse trends. These results are effective, model-conditional constraints, and are consistent with other recent multiphase simulations and observations. This highlights that the velocity-radius mapping encoded in UV absorption profiles enables recovery of outflow spatial structures from spatially integrated spectra. Our method paves the way for future broad parameter studies and guides updates of outflow simulations in future work.
Read the full abstract

All authors Xinfeng Xu, Drummond Fielding, Timothy Heckman, Greg L. Bryan, Alaina Henry, Karla Z. Arellano-Cordova, Cody Carr, John Chisholm, Claude-Andre Faucher-Giguere, Matthew Hayes, Mason Huberty, Michael Jennings, Crystal L. Martin, Claudia Scarlata, Allison L. Strom

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

Three-flavor supernova neutrino simulation using a hybrid quantum-classical algorithm with qutrits

We simulate a self-interacting three-flavor neutrino system within a core-collapse supernova using a hybrid classical-quantum algorithm on a qutrit computer. Based on the Dirac-Frenkel evolution equations, we employ a variation of the quantum-assisted simulator to calculate the system's time evolution operator by performing qutrit Hadamard tests to find expectation values of unitary operators in the Hamiltonian. The time evolution simulation is then done classically. We find that the hybrid algorithm produces results comparable to an exact numerical integration out to times of $t \approx 30 \,\omega_0^{-1}$ with time step $\delta t = 0.005 \,\omega_0^{-1}$, where $\omega_0$ is the energy scale of the single neutrino vacuum oscillations. We discuss the lessons learned in simulating neutrino systems using this hybrid quantum-classical algorithm, along with the advantages it offers over quantum Trotterization.
Read the full abstract

Astrophysics of Galaxies

Validating z > 7.5 Lyman Break Galaxy candidates in the COSMOS field with JWST/PASSAGE

We analyze spectroscopy from one NIRISS pointing in the JWST-PASSAGE program for seven candidate $z \gtrsim 7.5$ photometrically-selected COSMOS-Web sources. We spectroscopically confirm one out of seven sources as a Lyman break galaxy (LBG) at $z=7.962^{+0.003}_{-0.006}$, with $m_{F150W} = 25.9$ (AB). The remaining sources are too faint in the continuum (i.e., $m_{F150W} \gtrsim 26$ AB) to provide a redshift measurement from the Lyman break, and do not show emission lines in their spectra. Although this study contains only one spectroscopically confirmed source, the confirmation of a luminous $z \sim 8$ galaxy within this $\sim4.8$ arcmin$^2$ field implies a surface density of $\sim 0.21^{+0.59}_{-0.17}$ arcmin$^{-2}$, $\approx 10\times$ higher than inferred from wide-area photometric surveys, suggesting a potential overdensity at $z\sim8$ in the COSMOS field.
Read the full abstract

All authors Zoke W. Sackih, Mason S. Huberty, Claudia Scarlata, Peter J. Watson, Andrew J. Battisti, Farhanul Hasan, Matthew J. Hayes, Matthew A. Malkan, Vihang Mehta, Kalina V. Nedkova, Marc Rafelski, Benedetta Vulcani

Astrophysics of Galaxies

Analysis of the Space Velocities of Several Samples of Open Star Clusters

An analysis of the kinematics of open star clusters (OSCs) using their characteristics from the new Hunt and Reffert catalog was conducted. Based on 4003 OSCs younger than 200 million years, the following values for the angular velocity of the Galaxy's rotation were found: $\Omega_0 = 28.99\pm0.11$ km/s/kpc, $\Omega^{'}_0 = -3.909\pm0.026$ km/s/kpc$^{2}$ and $\Omega^{''}_0 = 0.5662\pm0.018$ km/s/kpc$^{3}$, where $V_0=234.8\pm3.0$ km/s for $R_0=8.1\pm0.1$ kpc. It was found that periodicity in the radial velocities of OSCs is manifested in clusters younger than 600 Myr, while a wave in residual tangential velocities is observed only in the youngest ones, younger than 40 Myr. A spectral Fourier analysis of the radial velocities of three OSC samples with average ages of 18, 72, and 143 Myr was used to obtain the following values of the wavelength $\lambda$ and the velocity perturbation amplitude $f_R$: $\lambda=2.0$ kpc and $f_R=4.3$ km/s, $\lambda=2.2$ kpc and $f_R=8.2$ km/s, $\lambda=2.1$ kpc and $f_R=9.6$ km/s, respectively. A systematic change in the positions of the maxima and minima of the waves in the radial velocities of OSCs was found depending on the age of the sample. From the analysis of these shifts, the value of the absolute value of the difference $|\Delta\Omega|$ between the angular velocity of rotation of the spiral pattern $\Omega_p$ and the rotation velocity of the Galaxy was found, $|\Delta\Omega|=2.0\pm0.5_{stat}\pm2.3_{syst}$ km/s/kpc. Based on this, an estimate of two possible values of the corotation radius was obtained: $8.6\pm0.2$ kpc and $7.6\pm0.2$ kpc, which indicates that the Sun is very close to the corotation.
Read the full abstract

Instrumentation and Methods for Astrophysics

Status of the COSmological Microwave Observations CALibrator

As the sensitivity of CMB telescopes increases, the need for precise calibration becomes critical. Started in 2022, the COSMOCal project aims to place an artificial polarized source in geostationary orbit, which will serve as a reference for CMB telescopes. This source will emit at 90, 150 and 270 GHz and will be linearly polarized with a highly precise orientation smaller than 0.1 deg. This proceeding presents the scientific motivations for the project, the current status of the development of the instrument and the results of a calibration campaign performed in March 2026 at the Institut d'Astrophysique Spatiale.
Read the full abstract

All authors L. Mousset, J. Aumont, S. Berta, L. Bizzarri, F. Boulanger, A. Catalano, F. Cuttaia, A. Denis, A.-L. Fontana, D. González Ovejero, B. Kalemi, S. Leclercq, J.-F. Macías-Pérez, B. Maffei, T. Mcnamara, S. Micheli, M. Migliaccio, L. Montier, P. Morfin, M. Moulin, M. Murgia, I. Myserlis, F. Nati, P. Ortu, M. Pérault, G. Pisano, T. Pisanu, N. Ponthieu, A. Ritacco, S. Savorgnano, L. Terenzi, L. Thibaut, J. Treuttel, C. Vescovi, L. Vacher, M. Zannoni

Cosmology and Nongalactic Astrophysics

Formation and Redshift Evolution of Dark Matter Spikes

Dark matter density spikes forming around adiabatically growing black holes can dramatically enhance indirect and direct detection signals. Canonical predictions, however, assume a zero-mass seed in a purely dark matter environment and do not track the long-term dynamical impact of surrounding stars. We present a semi-analytic framework that first generalizes adiabatic spike formation to include finite seed masses, stellar cusps, and non-circular orbits, and then studies the subsequent cosmic evolution by solving coupled Fokker-Planck equations for the dark matter and stellar phase-space distributions, with a heating rate modulated by the cosmic star formation rate. Starting conservatively from canonical Gondolo-Silk spikes and marginalizing over astrophysical uncertainties, we find that stellar gravitational heating drives the inner slope towards $\gamma_\chi \simeq 1.5$ within a few Gyrs (e.g by $z \lesssim 2$ for spikes formed at $z\simeq 10$), yielding overdensities two to four orders of magnitude below canonical expectations but still well above an NFW-like cusp. We provide redshift-dependent benchmarks for the column density and $J$-factor relevant to scattering, decay and annihilation signatures. Any robust interpretation of indirect dark matter signals from galactic nuclei must account for this evolution.
Read the full abstract

All authors Gonzalo Herrera, Abdelaziz Hussein, Lina Necib, Elliot Y. Davies, Xuejian Shen

Solar and Stellar Astrophysics

Millimeter-wave Detections of Symbiotic Stars in SPT and ACT Data

We present the results of a joint targeted search of candidate symbiotic stars at millimeter wavelengths using the South Pole Telescope (SPT) and the Atacama Cosmology Telescope (ACT). Candidates are selected from the New Online Database of Symbiotic Variables, restricting to objects that are within either the SPT-3G or ACT~DR6 footprint, covering most of the southern hemisphere and up to a declination of $+20^\circ$. Forced photometry on the 828 candidate symbiotic star locations in SPT and ACT data results in 31 unique objects detected with more than a $3\sigma$ significance using two frequency bands: 18 confirmed and 13 suspected symbiotic stars. We provide the SPT and ACT 95/98, 150, and 220~GHz light curves, along with optical and infrared light curves from 2016–2026, as well as spectral energy distributions, physical parameters from the literature, and brief summaries regarding the nature of each individual object. Using Herschel SPIRE data from 2013, we place upper limits on millimeter flux for CN Cha near the beginning of the optical rise in its 2012/2013 nova, which suggests a strong variability and lag at millimeter wavelengths and results in a rare observance of a Galactic millimeter slow transient. In addition, we provide coadded thumbnails and light curves for the remaining 797 candidate symbiotic stars that did not pass our detection thresholds. Millimeter-wave emission from symbiotic stars is primarily a combination of free-free emission of the ionization region and optically thick blackbody emission of the cooler dust components of the system. When combined with contemporaneous multi-wavelength observations, millimeter-wave observations can be used to test binary models of symbiotic stars and provide insight on the geometry and physical properties of these systems.
Read the full abstract

All authors C. Tandoi, A. Foster, T. J. Maccarone, A. J. Anderson, B. Ansarinejad, M. Archipley, L. Balkenhol, D. R. Barron, K. Benabed, A. N. Bender, B. A. Benson, F. Bianchini, L. E. Bleem, S. Bocquet, F. R. Bouchet, E. Camphuis, M. G. Campitiello, J. E. Carlstrom, J. Carron, C. L. Chang, P. M. Chichura, A. Chokshi, T.-L. Chou, A. Coerver, T. M. Crawford, C. Daley, T. de Haan, K. R. Dibert, M. A. Dobbs, M. Doohan, D. Dutcher, C. Feng, K. R. Ferguson, N. C. Ferree, K. Fichman, S. Galli, A. E. Gambrel, A. K. Gao, F. Ge, F. Guidi, S. Guns, N. W. Halverson, E. Hivon, G. P. Holder, W. L. Holzapfel, J. C. Hood, A. Hryciuk, N. Huang, T. Jhaveri, F. Keruzore, A. R. Khalife, L. Knox, K. Kornoelje, C.-L. Kuo, K. Levy, Y. Li, A. E. Lowitz, C. Lu, G. P. Lynch, A. S. Maniyar, E. S. Martsen, F. Menanteau, M. Millea, J. Montgomery, Y. Nakato, T. Natoli, A. Ouellette, Z. Pan, P. Paschos, K. A. Phadke, A. W. Pollak, K. Prabhu, W. Quan, M. Rahimi, A. Rahlin, C. L. Reichardt, M. Rouble, J. E. Ruhl, A. C. Silva Oliveira, A. Simpson, J. A. Sobrin, A. A. Stark, J. Stephen, C. Trendafilova, J. D. Vieira, A. G. Vieregg, A. Vitrier, Y. Wan, N. Whitehorn, W. L. K. Wu, M. R. Young, J. A. Zebrowski

Instrumentation and Methods for Astrophysics

The Outer Space Treaty Won't Save Us From Ourselves

The rapid growth of human activities in outer space sounds urgent alarms around ethical and philosophical issues, particularly concerning space militarization. The present international legal framework governing activities in space, the Outer Space Treaty (OST), views the peaceful exploration of space for scientific research as co-equal to other 'uses' entitled to "due regard" with respect to "potentially harmful interference" on the part of other space actors. The OST is deeply aspirational but has weak enforcement mechanisms, relying at its core on the goodwill of all involved parties as the fundamental basis for accountability. But that framework now faces unsustainable pressures from both public and private interests, and current agreements like the OST may be unable to exert timely, material protections. Terrestrial frameworks of "ethics of deterrence" versus the "ethics of agreements" are quickly expanding into cosmic environments. We argue for the legal recognition of space as an environment as the basis of any future approach to securing its integrity, and share examples of agreements grounded in peaceful cooperative stewardship of shared environments. These represent potential pathways forward that are ethical and also serve rational self-interest and self-preservation at this crucial juncture for humanity.
Read the full abstract

Earth and Planetary Astrophysics

Chemistry and Isotope Ratios of Substellar Atmospheres in the $\beta$ Pictoris Young Moving Group and Vicinity

Measuring the chemical and isotopic compositions of gas giants and brown dwarfs provides insights into their formation pathways and birth environments. 2MASS J0249-0557 c is an L2-type planetary mass companion ($\sim 12 M_{\mathrm{Jup}}$) orbiting a pair of brown dwarfs in the $\beta$ Pic young moving group and vicinity. Its mass places it at the intersection of planets and brown dwarfs, making it an interesting target for constraining formation pathways at the planet-brown-dwarf boundary. Using high-resolution spectroscopic data of the planet acquired with CRIRES+ mounted on VLT, we conduct atmospheric retrieval with the radiative transfer code petitRADTRANS and the nested sampling tool PyMultiNest. We retrieve a C/O ratio of $0.57\pm0.01$, a metallicity of [M/H] = $0.18\pm0.05$, and a $^{12}$CO/$^{13}$CO ratio of $95^{+23}_{-17}$. We also retrieve atmospheric compositions for two benchmark brown dwarfs in the $\beta$ Pic YMG, 2MASSI J0443+0002 and SIPS J2000-7523, using CRIRES+ data and find consistent compositions. Together with 2MASS J0249-0557 c's wide separation from its host, its compositional consistency with benchmark brown dwarfs supports gravitational collapse in a star-like manner as its most likely formation mechanism. These results deliver a homogeneous comparison of three substellar members in the $\beta$ Pic YMG and vicinity. Their solar-like abundances provide a baseline for exoplanet members in the same moving group, such as $\beta$ Pic b, 51 Eri b, and AF Lep b, whose host stellar compositions are difficult to measure. Future comparisons of atmospheric compositions among this moving group offer the potential to distinguish between formation mechanisms for its planetary members.
Read the full abstract

All authors Yurou Liu, Yapeng Zhang, Jerry W. Xuan, Dimitri Mawet, Ignas Snellen, Rico Landman, Tomas Stolker, Sam de Regt, Aurora Kesseli, Malena Rice

Earth and Planetary Astrophysics

A Coordinate System for Dynamical Instabilities in Hierarchical Systems in REBOUND

We implement coordinates suitable for studying wide binary systems in TRACE, a hybrid integrator in the widely used open-source N-body integration package REBOUND. This is a regime in which traditional hybrid integrators perform poorly. The coordinate system supports close encounters between any pair of bodies in the system. We describe the implementation of this coordinate system and benchmark its performance against other integrators in the REBOUND ecosystem. In tests of planet-planet scattering, stellar flybys, and ZLK oscillations. TRACE in wide binary coordinates is qualitatively correct when other hybrid methods fail, and in many cases returns statistically similar results to the high-precision IAS15 integrator with up to 9x speedups. We also provide some guidelines for when use of these coordinates are appropriate.
Read the full abstract

High Energy Astrophysical Phenomena

Neural Posterior Estimation for UHECR source inference from 3D propagation simulations

The identification of ultra-high energy cosmic ray sources is one of the open challenges of high-energy astrophysics. As charged particles travel through the Universe, they are deflected by extragalactic magnetic fields and lose energy through interactions with background radiation, making source inference highly non-trivial. Existing approaches either rely on simplified propagation models or on computationally prohibitive Monte Carlo methods. Here we present a simulation-based inference framework trained on three-dimensional CRPropa~3 propagation simulations that produces calibrated posterior distributions over source energy, distance, direction, and primary composition for individual UHECR events. The model combines a Deep Set encoder, handling the variable number of detected secondary particles, with a normalizing flow, and is trained on approximately 5 million simulated events covering a broad range of extragalactic magnetic field configurations. Validated on held-out simulations, all source parameters are recovered without systematic bias, with directional parameters best constrained and source distance most uncertain, consistent with the underlying propagation physics. Primary composition classification achieves $\geq$~98.2% accuracy across all mass groups. This framework provides a scalable and physically interpretable interface between detailed propagation simulations and Bayesian source inference relevant for current UHECR data.
Read the full abstract

High Energy Astrophysical Phenomena

Nonuniform Particle Injection into Black Hole Jets by Radiative Magnetic Reconnection

Active galactic nuclei often exhibit highly collimated relativistic plasma outflows launched from the vicinity of their central black holes. One of the key theoretical challenges in understanding black hole jet formation is the origin of the plasma that feeds the jet, which remains poorly understood, particularly in explaining the observed jet emission. In this study, we focus on electron positron pair production generated by high energy photons from non axisymmetric magnetic reconnection near the black hole, as suggested by recent three dimensional general relativistic magnetohydrodynamics simulations. By employing general relativistic ray tracing, we calculate the spatial distribution of the pair production rate in the jet, taking into account photon propagation and collision angles in curved spacetime. We find that our scenario can naturally supply a sufficient amount of plasma to explain the observed radio emission from the M87 jet, even when photon anisotropy is considered. Furthermore, we show that a spinning black hole plays a crucial role in shaping the spartial dsitribution of the pairs, which in turn affects jet acceleration and very high energy emission from the jet base.
Read the full abstract

Astrophysics of Galaxies

Prospects for Observing Galaxy Spectral Energy Distribution from the Radio to the far-Infrared in the Era of Next-Generation Radio Telescopes

The superb sensitivity and angular resolution of the next-generation radio telescopes with combined frequency coverage of approximately over three orders of magnitude (100 MHz–100 GHz) will sample the radio and far-infrared (FIR) spectral energy distribution (SED) of galaxies and revolutionize the galaxy formation study at the epoch of re-ionization and beyond. We present a prospect of observing the radio–FIR continuum SEDs of galaxies in the redshift of up to $z\approx 20$ based on an ensemble of the simulated `energy balanced' panchromatic SED (from UV to FIR) extended to the radio. For `realistic' populations of UV star-forming galaxies and dusty star-forming galaxies, we simulate their SEDs by accounting for the CMB effect and the radio–IR correlation. The flux density evolution of the UV-bright star-forming galaxies and the dusty star-forming galaxies at the selected observing frequencies covered by the current (ALMA) and next generation (SKA and ngVLA) radio-millimeter telescopes, suggest that massive galaxies (M$_* \gtrsim 10^{10}$M$_{\odot}$) are detectable at any redshift ($0<z<20$) in high frequency ($\nu>90$GHz). In particular, when operating, the ngVLA high-frequency ($\approx 100$ GHz) band is capable of detecting galaxies with M$_* \gtrsim 10^{9}$M$_{\odot}$ almost independently from redshift and the SKA low-frequency observing window ($\lesssim1$ GHz) has sufficient sensitivity to detect M$_* \gtrsim 10^{10}$M$_{\odot}$ dusty star-forming galaxies up to the epoch of reionization ($z=5\sim7$). We also show that the brightness of anomalous microwave emission (AME) in the galaxy SED is insignificant if the galaxies are beyond the local Universe (e.g., $z\gtrsim 0.1$).
Read the full abstract

All authors Ilsang Yoon, Jonathan Letai, Hansung B. Gim, Eric F. Jiménez-Andrade, Intae Jung, Caitlin Casey, Eric J. Murphy, Min S. Yun

Astrophysics of Galaxies

The complex history of NGC 1427A revealed by its star clusters and star formation history

Star-forming low-mass galaxies in the dense environments of galaxy clusters provide opportunities to study how environmental effects such as ram-pressure stripping, tidal interactions, or galaxy mergers shape a galaxy's star formation history. We combined integral-field spectroscopic observations with the Multi Unit Spectroscopic Explorer (MUSE) and available multi-band imaging of the star-forming galaxy NGC 1427A, located near the centre of the Fornax galaxy cluster, at a distance of 20 Mpc. Our aim was to trace the evolutionary history of NGC 1427A using the star formation history reconstructed from the integrated spectra and employing star clusters as surviving tracers of past star formation episodes. We fitted the spectral energy distribution of 222 star cluster candidates using archival $u,g,r$, and $i$ photometry to derive the ages and masses. For 58 clusters, we additionally incorporated their MUSE spectra in the fits and found an encouraging agreement between the photometric and spectroscopic results. The comparison of the age distribution of star clusters with star formation histories from a full spectrum fitting of the MUSE data found a reasonable agreement, with evidence for multiple episodes of star formation throughout the history of NGC 1427A. In particular, we found a population of young clusters ($\sim$ 10 Myr) that is located along the star formation edge and within the northern object, and a population of intermediate-age clusters ($\sim$ 100 - 300 Myr) with corresponding peaks in the star formation history of NGC 1427A. We interpret these populations in the context of the orbital evolution of NGC 1427A in the Fornax cluster and conclude that this galaxy has experienced not only ram-pressure stripping, but also tidal interactions or even a minor galaxy merger. The northern object is likely a regular component of the galaxy.
Read the full abstract

All authors Katja Fahrion, Michael Hilker, Avinash Chaturvedi, Juan P. Carvajal, Thomas H. Puzia

Instrumentation and Methods for Astrophysics

Recommendations for the Astronomy Graduate Admissions Process

As the AAS Working Group on Graduate Admissions (WGGA) we are sharing brief recommendations for improving and standardizing key elements of the graduate admissions process in astronomy. Most astronomy graduate programs have large areas of overlap in their admissions processes; however, the existing small variations in requirements and mismatches in communication and transparency make admissions more challenging for students and programs alike. To improve this situation, and building on the work presented in the AAS Graduate Admissions Task Force (GATF) report we recommend a few simple and straightforward changes for application content, communication, and timelines. These include an application format that consists of 1) two 500-word recommendation letters, 2) one 1500-word application essay, 3) an applicant CV, and 4) unofficial transcripts; and an admissions timeline that includes effective and transparent communication from programs and encouraging an April 1st "down-select date" for applicants.
Read the full abstract

All authors Emily M. Levesque, Courtney D. Dressing, Rachel Ivie, Grace Krahm, Meredith A. MacGregor, Daniel R. Piacitelli, Tom Rice

Cosmology and Nongalactic Astrophysics

Implications of SARAS3 data for Coulomb-like interacting dark matter

The 21-cm signal from cosmic dawn is a potentially sensitive probe of interactions between dark matter (DM) and baryons. We investigate the implications of the SARAS3 non-detection in the 55.5-84.4 MHz band for Coulomb-like interacting DM (IDM). In contrast to earlier constraint analyses that focused primarily on baryon cooling, we model the interaction self-consistently by including both excess cooling of the gas and the suppression of structure formation, which delays the onset of star formation and hence suppresses the Ly$\alpha$, X-ray, and ionizing backgrounds at early times. We perform a joint Bayesian fit of a global 21-cm signal model and a flexible foreground model to the SARAS3 antenna temperature, and find that the signal parameters remain weakly constrained after marginalizing over the foregrounds. The null result is nonetheless informative: the data disfavour deep absorption features within the observed band, with the strongest bound at $z = 23.6$ ($\nu \approx 57.7$ MHz), where $T_{21} \gtrsim -277.6$ mK at $3\sigma$. Comparing the IDM and standard cold dark matter scenarios, we find no statistically significant preference for IDM (Bayes factor $B \approx 1.7$). While we do not constrain the strength of baryon-DM interactions, the SARAS3 non-detection places a meaningful upper bound on the amplitude of the global 21-cm signal in this class of models.
Read the full abstract

All authors Shikhar Mittal, Prakhar Bansal, Harry Bevins, Saurabh Singh

Astrophysics of Galaxies

NEFERTITI: Linking early galaxy formation to the assembly of the Milky Way

We use a new implementation of the NEFERTITI galaxy formation model, coupled to $\sim 30$ high-resolution Caterpillar dark-matter simulations of Milky Way (MW) analogues, to connect early galaxy formation with the MW's assembly down to $z=0$. Our locally-constrained model resolves minihaloes hosting the first PopIII stars and self-consistently tracks inhomogeneous ionization and chemical enrichment. PopIII star formation begins at $z\simeq27$, peaks at $z\simeq10-15$, and persists down to $z\lesssim5$, producing PopIII systems with $M_*\sim10-5\times10^5\:{\rm M_\odot}$. The present-day descendants of PopIII stars span ${\rm [Fe/H]<-9}$ to ${\rm [Fe/H]\approx-1}$, with the most metal-poor stars typically enriched by a few (1-4) low-energy supernova progenitors. Pair-instability supernova descendants more commonly form in massive haloes ($M_{\rm vir}>10^8\:{\rm M_\odot}$), often externally enriched, reflecting the strong feedback and delayed recovery following energetic explosions. These early systems serve as building blocks for the present-day Galaxy's metal-poor component: although 90$\%$ of the total stellar mass formed in situ, the accreted component dominates at $[{\rm Fe/H}]<-1$ and accounts for nearly all stars with $[{\rm Fe/H}]<-3$. This accreted population is largely built by a few ($\sim5$) massive ($M_*>10^8\:{\rm M_\odot}$) destroyed dwarfs, but lower-mass systems become increasingly important at low metallicities, with ultra-faint and classical dSph analogues contributing $\sim25\%$ at $[{\rm Fe/H}]<-3$. Our model simultaneously reproduces the properties of metal-poor MW stars and the JWST "Hebe" galaxy at $z\sim11$, supporting its identification as a pure PopIII system. Ultimately, NEFERTITI is a key tool to interpret upcoming local and high-$z$ observations linking the near- and far-field cosmology.
Read the full abstract

All authors Ioanna Koutsouridou, Stefania Salvadori, Ása Skúladóttir, Viola Gelli, Elka Rusta, Lapo Querci, David S. Aguado, Alice Mori

Astrophysics of Galaxies

Tidal pre-conditioning and ram-pressure stripping in NGC 1427A. Deep VLT/MUSE spectroscopy and FUV-to-radio observations trace a Fornax Cluster dwarf in transformation

The early environmental transformation of low-mass cluster galaxies from gas-rich to gas-poor remains poorly constrained, partly because clear, phase-resolved observations are rare. NGC 1427A, a disturbed star-forming dwarf in the Fornax cluster, offers a favorable case for studying this process. We aim to build a spatially resolved, multi-phase picture of NGC 1427A to constrain the roles of ram-pressure stripping and tidal perturbations. We combine a deep, spatially contiguous VLT/MUSE mosaic with ancillary data from the FUV to the radio. Full-spectrum fitting of the MUSE cube yields maps of stellar kinematics, ages, metallicities, and continuum attenuation, while emission-line modeling provides ionized-gas kinematics, Balmer-decrement reddening, and star-formation-rate surface densities. Ancillary multi-wavelength data trace neutral and molecular gas, dust, and recent star formation, placing the MUSE-based results in a broader multi-phase context. We find a pronounced decoupling between stars and gas: the H I and ionized gas rotate about an axis tilted with respect to the stellar field and are globally blueshifted. Stellar and nebular attenuation, infrared dust tracers, and H I morphology indicate stripping with a strong line-of-sight component that has reached the ISM. At the same time, the asymmetric gas and dust distribution, together with structured and time-dependent star formation, points to an additional gravitational perturbation, with a recent mild fly-by by a nearby dwarf being the favored interpretation. We propose that dwarf-dwarf tidal effects have torqued and pre-conditioned the gas, while the Fornax intracluster medium is driving ram-pressure stripping that now reaches the ISM and coincides with a declining global star formation rate. This places NGC 1427A at the onset of environmental quenching, making it a useful benchmark for early cluster dwarf transformation.
Read the full abstract

All authors J. P. Carvajal, T. H. Puzia, M. Blaña, P. K. Nayak, K. Fahrion, M. Hilker, E. Emsellem, E. D. Skillman, M. D. Mora, B. Miller, Y. Ordenes-Briceño, P. Eigenthaler, R. Rahatgaonkar, J. Chanamé

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

Series solutions to the TOV equations

We present general series solutions to the Tolman-Oppenheimer-Volkoff equations for compact stellar objects. We develop an algorithm to compute the coefficients of the power series in terms of the equation of state and its derivatives with respect to the thermodynamic variables. Using these results, we establish general properties of analytic solutions and their relation to the regularity of the equation of state. Applying the theory of Padé approximants, we derive series representations for meromorphic functions whose domains of convergence may include isolated poles. These analytic solutions are then used to obtain closed-form expressions to approximate the radius and mass of stellar objects. We apply the formalism to specific models, namely fluids with affine equations of state and polytropic fluids, and compare the results with those obtained from numerical integration. Lastly, we extend the formalism to piecewise equations of state, deriving series solutions that can be matched across transition hypersurfaces.
Read the full abstract

High Energy Astrophysical Phenomena

On the isotropy of viscosity in accretion discs

Accretion discs are fundamental to many astrophysical systems, providing the conversion of gravitational potential energy into radiation that we can observe. In many systems there is evidence that discs are warped; from spatially-resolved observations of protoplanetary discs, to the features of lightcurves and line profiles from discs around supermassive black holes in galaxy centres. The dynamics of warped discs is largely controlled by the physical nature of the internal disc viscosity. While typically disc viscosity is hydromagnetic in origin, simulations of magnetized discs cannot match observed rates of angular momentum transport in planar discs and thus cannot be used to determine the ratio of the torques responsible for driving accretion to those responsible for evolving the disc warp. The analytic work of Ogilvie is the most comprehensive model for warped disc evolution, but makes assumptions that need to be tested. In particular, it assumes that the disc viscosity is Navier-Stokes, and therefore small-scale and isotropic. Here we attempt to test this model using the long periods of X-ray binaries that are due to precession of the disc. These systems have well-constrained estimates of the component of viscosity responsible for driving accretion, and by looking at systems with and without evidence for disc misalignment and precession we can constrain the component of viscosity responsible for flattening the disc. We conclude that the observational constraints suggest that the Ogilvie model provides an adequate description of the disc evolution, but that there are indications that the internal disc viscosity might be marginally non-isotropic.
Read the full abstract

Astrophysics of Galaxies

Merge and Strip II: Imprint of galaxy formation physics and viscosity on baryon-dominated dwarf galaxies

Motivated by the discovery of peculiar dwarf galaxies inside galaxy clusters such as blue candidates (BCs), dark galaxies and ultra-diffuse galaxies (UDGs), we present hydrodynamic simulations of galaxy mergers in cluster environments. We vary the viscosity and stellar feedback prescriptions, realistically modelling possible conditions for hydrodynamic drag and fluid instabilities, as well as internal destabilization through stellar feedback-driven heating and gas loss. We find that long-lived tidal dwarf galaxies (TDGs) can form throughout all viscosity values applicable to galaxy clusters if stellar feedback is moderate. Our results expand on studies of cloud crushing simulations, investigating the entrainment problem in intracluster medium ambience. The smallest clouds have gas masses on the order of $M_\text{gas} \sim 10^7 \text{ M}_\odot$ and reach relatively low final drift velocities of $\sim 100 \text{ km/s}$. The lowest possible Reynolds number acting on this class of clouds is $Re \sim 1$ for full Spitzer viscosity. Almost all TDGs display elevated star formation rates of $0.01-0.1 \text{ M}_\odot / \text{yr}$, which are stable across several Gyr. Based on their matching properties, we support that BCs observed in the Virgo cluster are likely stripped TDGs. Similar features are also found in comparison with dark galaxies and baryon-dominated UDGs, implying that a subsample of these objects are also long-lived TDGs. This work provides robust evidence that stripping from galaxy mergers is a viable channel for the formation of stable cold gas clouds and dark matter-deficient galaxies observed in galaxy clusters.
Read the full abstract

All authors Anna Ivleva, Klaus Dolag, Rhea-Silvia Remus, Duncan A. Forbes, Tirso Marin-Gilabert

Astrophysics of Galaxies

Stellar mass and morphology segregation in pairs and multiplets in the cosmic web

In this work, we investigate whether the location of galaxies within the large-scale structures (LSS) of the Universe affects their stellar mass ($M_\star$) and morphology. To this end, we attempt to disentangle the effects of local and large-scale environments in their distributions. We classify 25309 galaxies in the redshift range ${0.02 < z \leq 0.04}$ with $\log M_\star/\rm{M}_\odot \geq 9.5$ in terms of the main LSS (voids, clusters, and not clusters nor voids, referred to as NCNV) and local environment (singlets and multiplets; galaxies with and without companions). We present the stellar mass and morphology distributions in these environments, and for a subsample of galaxy pairs. Even in voids, we find that $\sim22\%$ of galaxies have companions. Stellar mass distributions show that galaxies are less massive in voids, regardless of their local environment. Satellites in voids are, too, less massive relative to their centrals than in NCNV pairs. In terms of morphology, the denser the LSS, the greater is the proportion of early-type galaxies, even among singlets. In voids and NCNV, late-type multiplets tend to be later-type spirals than singlets. In pairs, centrals tend to be more early-type than satellites. The sample, curated to avoid morphology incompleteness, yields slightly higher fractions of early-type galaxies and multiplets than previous studies. We conclude that the local environment alone is insufficient to explain the distribution of stellar mass and the morphology of galaxies in the local Universe. The observed mass distributions support a scenario in which galaxy assembly depends critically on the host halos, and the properties of these halos are related to their large-scale environment. This would explain the finding of lower-mass galaxies in voids than in denser environments, and provide a basis for considering a common evolutionary origin for multiplets.
Read the full abstract

All authors G. Torres-Ríos, S. Verley, I. Pérez, M. Argudo-Fernández, B. Bidaran, S. Duarte Puertas, Y. K. González-Koda

Cosmology and Nongalactic Astrophysics

Coverage is not enough: Frequentist tests of simulation-based inference for primordial non-Gaussianity

(Abridged) Simulation-based inference (SBI) has emerged as a powerful framework for extracting cosmological information from complex, non-linear data where analytical likelihoods are unavailable. Its reliability is commonly assessed using coverage-based diagnostics under the prior predictive distribution, which probe calibration only in an averaged sense and do not constrain posterior behavior at fixed parameter value, the regime relevant for practical inference. We investigate these limitations in the context of primordial non-Gaussianity, parameterized by $f_\mathrm{NL}$, using simulations of the dark matter halo field. We compare SBI based on contrastive neural ratio estimation (CNRE) with likelihood-based inference (LBI) using the power spectrum, bispectrum, and wavelet scattering transform (WST) coefficients across 1000 realizations. SBI and LBI agree well on posterior means and skewness, while the variance agrees on average but shows weaker realization-by-realization consistency. Larger differences arise in the kurtosis, indicating discrepancies in the posterior tails. These effects are already present for the power spectrum - where the Gaussian likelihood assumed in LBI is best justified - and are most pronounced for the combined power spectrum and bispectrum, where SBI posteriors are often underconfident and can yield weaker constraints than either statistic individually, despite passing coverage tests. WST coefficients further tighten constraints on $f_\mathrm{NL}$, even when restricted to large scales. Our results highlight both the potential of higher-order statistics and the need for validation strategies that probe the posterior shape beyond standard coverage diagnostics.
Read the full abstract

Astrophysics of Galaxies

Discovery of Quasar Variability and Early Accretion Disk Signatures at Cosmic Dawn

In the nearby universe, quasars are well known to exhibit variability in their brightness over time, offering a powerful tool to probe the physics of accretion onto the SMBH and directly measure the mass of the SMBH. However, detecting variability in early quasars remains challenging. Here, we report the detection of multi-wavelength infrared and X-ray variability in a quasar observed just 850 million years after the Big Bang. The infrared variability spans five filters, tracing rest-frame ultraviolet and optical emission from the accretion disk, while the X-ray variability probes the corona. The variable spectrum reveals that the accretion disk has a geometrically thin, optically thick structure. This provides observational constraints on the accretion disk structure at early times, when quasars are accreting at high Eddington ratios and reside in extreme environments. Our findings demonstrate the feasibility of characterizing accretion physics using variability in the early universe, laying the groundwork for studies exploiting upcoming facilities such as the Rubin Observatory and Roman Space Telescope. These facilities will discover large samples of variable high-redshift quasars, enabling population-level variability studies of accretion physics and black hole masses, filling key missing ingredients in understanding early SMBH growth.
Read the full abstract

All authors Gene C. K. Leung, Anna-Christina Eilers, Christos Panagiotou, Julien Wolf, Kishalay De, Luke Weisenbach, Minghao Yue, Xiaohui Fan, Yuzo Ishikawa, Erin Kara, Mirko Krumpe, Andrea Merloni, Robert A. Simcoe, Feige Wang, Jinyi Yang

Astrophysics of Galaxies

Self-acceleration of Hardening Binaries

A Keplerian binary immersed in a bath of lighter particles hardens by ejecting them through gravitational slingshots. This process drives, for example, the evolution of supermassive black hole binaries following galaxy mergers, and has long been described with just two parameters: the hardening rate and the eccentricity growth rate. Here we show that the secular dynamics is substantially richer. Combining symmetry arguments with extensive three-body scattering experiments, we demonstrate that the medium exerts a net force on the binary's center of mass (CoM), induces apsidal precession, and rotates the orbital plane when the CoM velocity has an out-of-plane component. Remarkably, these deterministic effects persist even in a perfectly uniform and isotropic medium, as the binary's own asymmetry provides the propulsion. The interplay of self-acceleration, precession, and dynamical friction drives the CoM along an outward spiral. For supermassive black hole binaries, this displacement dominates over Brownian motion and approaches the radius of influence, suggesting they may be significantly offset from their host galaxies' centers. The displacement also enlarges the stellar loss cone, with direct implications for the final-parsec problem. We further show that the previously reported circularization of small-mass-ratio binaries is a numerical artifact of truncating long-lived encounters: all binaries undergo eccentricity growth. Our results enrich the standard picture of binary hardening and have implications in a variety of astrophysical contexts, including gravitational-wave source populations.
Read the full abstract