Astrophysics papers — 2026-10-09

The Morning Preprint — this day as one edition: the day's papers, the five things worth knowing, and the episodes behind them.

Today we are diving into the work that matters most concerning our cosmic microwave background polarization data. This involves better separating foreground contamination from the faint cosmological signals to understand early universe physics.

The main effort is creating a new JAX-powered implementation for parametric component separation and clustering optimization for CMB polarization satellites. This method handles spatially varying foreground spectral energy distributions by scanning thousands of pixel subset configurations AKMPC. This approach achieves up to a hundred times speed-up over older optimizers when run on GPUs. It has shown promise in reducing the upper limit on the tensor-to-scalar ratio r by about thirty percent in LiteBIRD-like simulations.

Building upon this acceleration, there is work providing tools for modeling large scale structure, such as Sym-EFT: Accelerating Effective Field Theory of Large Scale Structure with Symbolic Regression Sym-EFT. This tool emulates the one- and two-loop cold dark matter power spectrum using symbolic regression to create versatile functions for fitting counterterms. This emulator suite is fast, achieving evaluations in less than five times ten to the fourth seconds on a single core while maintaining errors better than zero point five percent within its valid range.

On a different scale, research looks at how dark matter models affect the Lyman-alpha forest using Lyman-alpha Forest Signatures of Mixed Fuzzy and Cold Dark Matter Lyman-alpha Forest Signatures of Mixed Fuzzy and Cold Dark Matter. This study compares different dark matter treatments, finding that the flux statistics depend on the velocity prescription rather than just the matter power spectrum. These differences grow toward higher redshift and smaller scales.

This contrasts with work focusing on galaxy clusters, such as Cosmic gas accretion from filaments onto galaxy clusters using the IllustrisTNG simulation Cosmic gas accretion from filaments onto galaxy clusters using the IllustrisTNG simulation. This highlights how filamentary structures drive anisotropic gas inflows into massive systems.

Finally, there are insights into high redshift star formation from JADES: Evolution of nitrogen abundances in star-forming galaxies from z 1.5-7 JADES. These measurements using JWST data reveal significant scatter in nitrogen abundance. Extreme enhancements are tied to the most intense starbursts, which is important context for understanding chemical enrichment in these early systems.

The work concerning gravitational wave emission from nonspherical collapse in an early matter-dominated era is particularly important. It offers a way to probe the pre-Big Bang thermal history by mapping present-day gravitational wave spectra to detector sensitivity bands. This potentially reveals information about reheating temperatures.

This study uses a semirelativistic N-body framework to evolve the collapse of an overdense patch initialized via a Zel'dovich deformation in an Einstein--de Sitter background. It shows that reliable signal prediction demands fully numerical treatment of nonlinear dynamics.

A key finding is that fitting-based procedures and Zel'dovich-based estimates fail to capture post-shell-crossing evolution, leading to incorrect power estimations for gravitational waves. After averaging results from Doroshkevich and BBKS distributions, the two spectra show similar shapes at the peak amplitude. However, the BBKS result is systematically smaller because its dominant contribution comes from peaks of modest height around nu three. This work connects to a broader understanding of cosmological structure formation as one can map these GW spectra across different frequency ranges to constrain early universe physics.

In contrast, research into galaxy clusters focuses on disentangling turbulence from bulk motions in the intracluster medium. This is crucial for interpreting observational data and simulations. A paper using GPU-accelerated filtering pipelines on a major galaxy cluster merger simulation found that turbulent pressure fraction on physical scales of fifty kiloparsecs peaks at five percent during the merger. It then drops to two percent after about one point three gigayears, suggesting turbulence levels remain low unless a cluster has been recently perturbed by a major merger, which aligns with recent XRISM observations.

Another area of investigation involves understanding wave propagation in astrophysical environments. This looks at how solar vortices act as conduits for magnetoacoustic waves. Using information-theoretic metrics like Normalised Mutual Information and Jensen-Shannon Divergence, the authors provided the first observational identification of a solar vortex footpoint by locating it where there is highest shared information with the chromospheric swirl. The Bifrost simulation indicated that vortices primarily support compressive magnetoacoustic waves in the lower atmosphere, and these modes carry a thermal signature relevant for localized heating.

Furthermore, analysis of binary black hole mergers in gravitational wave catalogs reveals distinct mass regions characterized by sharp transitions in mass-ratio and effective spin properties. For example, below fifteen solar masses, the population favors equal-mass binaries with a narrow positive spin distribution. This suggests that different primary-mass ranges trace distinct formation channels, linking the merger statistics to hierarchical assembly processes.

The work on SAETASS is particularly important because it offers a novel, open-source numerical tool for solving the time-dependent transport equation in spherical symmetry. This solver is crucial for modeling how non-thermal particles move through astrophysical plasmas like supernova remnants. It uses a conservative finite-volume framework and an operator splitting architecture to handle advection, losses, and diffusion simultaneously.

The application of this solver to cosmic ray proton transport in a real astrophysical scenario shows that it can recover established steady-state limits while also revealing pre-equilibrium temporal dynamics across different diffusion regimes like Kolmogorov, Kraichnan, and Bohm. This capability is complemented by the work on AMS-02 proton and antiproton fluxes during solar minimum. This develops a unified charge-dependent solar modulation model by solving the three-dimensional Parker transport equation.

This modulation model fits observed proton and antiproton fluxes using a local interstellar spectrum constrained by Voyager data. It employs neural network surrogate models to speed up parameter scans, providing a single framework for understanding how solar activity affects charged particle propagation across different charge states.

Another piece of research addresses the hadronic origin of sub-PeV gamma-ray emission from pulsars. This explores an alternative explanation based on proton-proton interactions within a one-zone superdiffusive environment characterized by a diffusion index of one point zero five. This scenario requires a cosmic ray proton luminosity about zero point one four times the spin-down luminosity of PSR J0622+3749 to model the observed gamma-ray spectrum.

Meanwhile, the study on wide binary stars tests Modified Newtonian Dynamics by measuring differential radial velocities with high precision. These measurements achieve improvements over Gaia DR3 in resolving velocity differences. This technique allows researchers to constrain the global MOND acceleration scale by combining these kinematic measurements with astrometric data.

The analysis of Laniakea also shows that modeling the bias induced on cosmological observables by this supercluster structure can shift inferred values for the Hubble constant, potentially worsening existing tensions in cosmology. This work uses peculiar velocity data from Cosmicflows-4 to model the triaxial expansion of Laniakea and its impact on luminosity distance corrections for SN Ia and SBF data.

The most significant finding today relates to how we can probe fundamental physics through cosmic structure. This involves detecting the Integrated Sachs Wolfe effect or ISWRS using cross-correlations between cosmological probes. This work explores linking current and future Cosmic Microwave Background experiments like Simons Observatory and CMB-S4 with large scale structure surveys such as Euclid and LSST to forecast results.

The analysis models how the ISWRS effect correlates with gravitational potential tracers like galaxy clustering, cosmic shear, and CMB lensing potential to test the neutrino mass sum M nu in the ν Λ CDM model. This suggests that future experiments might be able to differentiate between different cosmological models by observing this subtle effect. While CMB-HD's reconstructed lensing potential is expected to yield strong detections even with conservative noise assumptions, other probes like galaxy clustering can also provide significant signals, although cosmic shear results are limited by the need for accurate non-linear modeling.

Building on this cosmological work, there is research into how different dark matter models affect structure formation. A study on halo mass functions in mixed cold and fuzzy dark matter cosmologies investigates the impact of an ultralight axion-like component with a mass between ten and twenty-four point five electron volts constituting up to three tenths of the total dark matter. This work shows that this fuzzy dark matter traces the large scale distribution while suppressing small scale structure, which leads to a reduction in low mass haloes and modifies the halo mass function depending on redshift and the fraction of fuzzy dark matter.

This framework is then used to create a phenomenological model that maps cold dark matter halo mass functions to their mixed cold and fuzzy counterparts. This reproduces simulated results within about zero point one to zero point two dex across the explored parameter space. This method offers a computationally efficient way to predict structure formation in these mixed models without needing dedicated simulations for every parameter combination.

In terms of instrumental modeling and inference, there is work identifying an intrinsic degeneracy in simplified semi-analytic two spot models for thermal X-ray pulse-profile modeling. This matters because it suggests that even a simplified model can produce two peaks on the likelihood surface when more complex factors like Doppler effects are included. This analytic degeneracy might offer intuition for multi modal posterior distributions seen in other studies.

This finding is less significant than the cosmological or dark matter structure work, but it points to limitations in parameter inference methods. The paper suggests that while the exact algebraic degeneracy exists without Doppler effects, including Doppler effects, statistic errors, and instrument response can still result in two peaks on the likelihood surface. This analytic degeneracy may provide intuition for multi modal posterior distributions discovered in previous works.

Finally, there is work concerning nuclear physics parameter inference using semi agnostic priors for radio pulsar timing and X-ray pulse profile modeling. This research assesses how detections of mass, radius, and tidal deformability from neutron stars constrain nuclear empirical parameters when using semi agnostic equation of state priors. The study shows that not all nuclear empirical parameters significantly correlate with the pressure in different density regimes, which challenges their independent inference.

The work concerning the peculiar motions of galaxies is particularly important because it offers a way to trace how structures grow and where matter is distributed throughout the universe. This helps us probe dark energy and gravity on large scales. The paper AKMPC presents a Bayesian method to estimate the radial component of peculiar velocities for galaxies hosting Type Ia supernovae, doing this without assuming local linearity or a fixed cosmology. This means it avoids bias even when peculiar velocities are large and is consistent with different cosmological models.

This approach builds upon earlier work that looks at how cosmic ray acceleration affects the early evolution of bow shocks around massive runaway stars. AKMPC uses 3D ideal cosmic ray magnetohydrodynamic simulations to see how particle acceleration from these shocks changes over time. It finds that variations in cosmic ray diffusion rates strongly dictate the morphology and non-thermal emission of the bow shock.

Another piece of research explores the impact of finite-temperature quantum gravity effects on cosmological parameters using artificial neural networks. AKMPC modifies systems to incorporate temperature-dependent quantum corrections into the Hubble parameter. It reveals that these new parameters improve fits to 2018 Planck data, suggesting potential refinements for cosmological models despite the persistent Hubble tension.

Moving toward more specific astrophysical phenomena, the study on tidal disruption of stellar binaries provides a robust channel for generating diverse transient events. AKMPC uses simulations to show that tidal separation of stellar binaries by a supermassive black hole can naturally produce eccentric Tidal Disruption Events, classifying them into elliptical or hyperbolic types based on the outcome.

We also have work investigating potential dark matter signals near supermassive black holes, specifically looking for heavy WIMPs in M31. AKMPC analyzes detection prospects using very high energy gamma-ray bands, estimating that the Cherenkov Telescope Array could probe a significant portion of TeV-scale WIMP parameter space in M31 if spike density is optimistic.

Finally, there is research into cometary activity in the outer Solar system. AKMPC models comet activity at large heliocentric distances by proposing that the recombination of radicals causes outbursts. This could explain observed dust contributions to the cosmic extragalactic background light and influence the orbits of small-radius comets in the Oort cloud.

The most significant finding relates to how we can search for neutrino and Hot Dark Matter wakes using 3D maps of a tracer of HDM. This effect is important because it provides a smoking gun for free-streaming, which cannot be faked by changes to the background expansion history, offering a complementary way to find massive neutrinos or hot dark matter components beyond traditional observables.

The work by Nascimento and Loverde on Dark winds on the horizon: Prospects for detecting neutrino and hot dark matter wakes in large-scale structure shows that while these wakes are unlikely to be seen with the most natural tracer, 2D weak lensing surveys can achieve high significance with idealistic 3D maps of an HDM tracer if the effective free-streaming length is sufficiently small. This length needs to be around zero point zero three Mpc inverse to reach a signal-to-noise ratio of one for a dark matter species making up one percent of the total dark matter. This contrasts with the idea that these wakes are generally undetectable using standard tracers.

This potential detection hinges on the nature of free-streaming, which is what AKMPC explores in detail regarding cosmological signatures. Furthermore, this effect serves as an alternative search avenue for hot subcomponents of dark matter beyond just massive neutrinos, complementing traditional observables.

Today's papers

The papers

Important terms

Parametric Component Separation
A new JAX-powered method for separating foreground contamination from faint cosmic microwave background polarization signals by scanning thousands of pixel configurations.
Sym-EFT
A tool using symbolic regression to emulate one- and two-loop cold dark matter power spectra, enabling fast fitting of counterterms for large scale structure models.
Lyman-alpha Forest Signatures
Research comparing different dark matter treatments in the Lyman-alpha forest, showing that flux statistics depend on velocity prescription rather than just the matter power spectrum.
Integrated Sachs Wolfe Effect (ISWRS)
A technique using cross-correlations between CMB experiments and large scale structure surveys to probe fundamental physics, potentially constraining neutrino mass sum.
Dark Matter Halo Mass Functions
Modeling how mixed cold and fuzzy dark matter affects structure formation, showing that fuzzy components suppress small scale structures while tracing large scales.