Instrumentation and Methods for Astrophysics
Agentic research is oxymoronic
The use of agentic large language models obviates human interpretation of scientific results, and will lead to substantial distrust in the literature.
Papers
Fresh from arXiv, typeset properly. Highlighted words trade on the market.
Monday 31 August
Instrumentation and Methods for Astrophysics
Ad Astra White Paper: A Pitch for the Next 25 Years of NASA's Physics of the Cosmos Program
Astrophysical observations of our universe have been key to our understanding of how the universe works. Shortly after the turn of the millennium, the National Research Council delivered Connecting Quarks with the Cosmos: Eleven Science Questions for the New Century. In the subsequent quarter-century, we have made substantial progress in answering each question. These advancements have, in part, arisen because of the success of major US facilities across several domains of physics, guided by long-term planning documents which still largely focus on these questions. This report seeks to provide a status update on each question, and to outline what space-based facilities are crucial for future progress, intended to guide NASA's preparatory work for the Astro2030 Decadal.
Read the full abstract
Instrumentation and Methods for Astrophysics
Slysh haloes: the waste heat of cold computing as a submillimetre technosignature
Searches for Dysonian waste heat have operated almost exclusively in the mid-infrared and are therefore sensitive primarily to technology radiating at 100-600K. We argue that mature, computation-dominated civilisations may instead dissipate much of their energy at far lower temperatures. The Landauer cost of irreversible computation scales linearly with temperature, while ambient temperatures at large circumstellar radii approach the 2.7K cosmic microwave background floor. Cold computation is therefore thermodynamically attractive and, because the required radiating area scales as T**-4, potentially conspicuous. These considerations predict a new object class, which we term the SLYSH HALO after the first advocate of cold Dysonian searches. A Slysh halo is a physically motivated partial Dyson swarm producing grey, line-free thermal emission from the cold outer regions of planetary systems. We show that such structures are energetically and materially plausible, and that M dwarfs provide especially favourable search targets. Archival far-infrared and submillimetre surveys of nearby stars (DEBRIS, DUNES and SONS) can in principle be reinterpreted to constrain cold circumstellar dissipation at approximately the 10**20 Watt level, several orders of magnitude below the waste-heat luminosities targeted by previous infrared searches. Additional opportunities are provided by archival observations from Planck and the JCMT, together with the reprocessing of interferometric data from facilities such as ALMA and NOEMA. We assemble six observational discriminants that separate engineered radiators from natural cold sources and outline a three-tier search programme. Even a null result would provide the first temperature-complete assessment of Dysonian technosignatures.
Read the full abstract
Instrumentation and Methods for Astrophysics
The Analysis, not the Aperture: End-to-End Transformer Reconstruction for Imaging Atmospheric Cherenkov Telescopes
Imaging Atmospheric Cherenkov Telescopes (IACTs) detect very-high-energy gamma rays by imaging the nanosecond Cherenkov flash of the air shower they initiate in the Earth's atmosphere. For four decades the first steps of IACT event reconstruction have been essentially unchanged, relying on a heavy parameterisation and dimensionality reduction of the recorded images. This is reasonable when the image is bright, but discards important information when only a few tens of Cherenkov photons are recorded, which is a primary reason why small telescopes perform poorly at sub-TeV energies. We show that this limitation is a property of the analysis rather than of the hardware. We simulate a deliberately simple and idealised compact telescope and treat each event as a short movie that is passed directly to a video vision transformer with a factorised spatio-temporal encoder. A single composite network with a gradient-normalised multi-task loss performs gamma/hadron classification, energy regression and arrival-direction regression at once. This is the first application of a video vision transformer to IACT data. We compare it against an optimised standard analysis on the same dataset. The transformer lowers the energy threshold by a factor of three, from 0.22 to 0.07 TeV, and reconstructs arrival directions down to 0.05 TeV. At 0.2 TeV it increases the effective collection area by a factor of three, and raises the gamma/hadron separation power from an area under the receiver operating characteristic curve of 0.80 to 0.91. At 0.05 TeV, where the standard analysis retains almost nothing, that area grows by nearly two orders of magnitude. These results show promising new opportunities for compact and affordable telescopes operating at sub-TeV energies, paving the way for a broader exploration of time-domain astrophysics.
Read the full abstract
Astrophysics of Galaxies
The impact of hot mode accretion and angular momentum conservation on metallicity gradients of galactic discs
The hot circumgalactic medium (CGM) of a galaxy inevitably rotates more slowly than the cold gas in the disc it surrounds, due to a higher pressure support against gravity. If it accretes vertically onto the disc, angular momentum conservation leads to radial flows, advection of metals inwards and the steepening of a metallicity gradient. The observed gradient hence carries information on the kinematics of the accreting gas. Previous models built on this premise parametrised the angular momentum mismatch as a simple function of radius, constant in time, rather than derived from the properties of the CGM itself. Here, we present a semi-analytic chemical evolution model of a galaxy with a cosmologically motivated evolution, a self-consistent and time-evolving potential and an isothermal rotating CGM. The model includes inside-out formation and the radial flows induced both by the deepening of the potential well (acting on gas and stars) and by accretion, following their combined effect on the gas-phase metallicity gradient over cosmic time. We use a Bayesian framework to fit simultaneously the structural properties and metallicity gradient of a galaxy, inferring both the required inside-out growth and the CGM temperature needed to drive the required radial flows and abundance gradient. Fitting our model to the Milky Way, we find a mild preference for a subvirial corona ($T_\mathrm{CGM}/T_{200}=0.82^{+0.44}_{-0.38}$), in mild tension with observations, and we discuss possible resolutions. The predicted CGM rotation velocity lies between ~100-160 km/s at the present time and is substantially lower at earlier epochs. We further recover an evolution of the gradient over the past ~8Gyr consistent with observations. Our approach is not restricted to the Milky Way and can be directly applied to external galaxies, offering a new handle to constrain hot mode accretion onto galactic discs.
Read the full abstract
Cosmology and Nongalactic Astrophysics
SPT-3G D1: Quadratic-Estimator CMB Lensing Reconstruction and Cosmology
We present a map of the cosmic microwave background (CMB) lensing potential reconstructed from observations taken during the 2019 and 2020 seasons with the third-generation camera on the South Pole Telescope (SPT), covering the $1500\,{\rm deg}^{2}$ SPT-3G Main field, referred to as the SPT-3G D1 dataset. From the multi-frequency temperature and polarization data, we reconstruct the CMB lensing field using a quadratic estimator that jointly accounts for the $T$, $E$, and $B$ fields and their covariance. The resulting lensing map is dominated by polarization information for $L \lesssim 600$ and provides the highest signal-to-noise measurement per mode reported to date. With nuisance parameters fixed to their best-fit values, we measure a lensing amplitude consistent with unity at $2\%$ precision relative to the $\Lambda$CDM model that best fits the combined Planck, ACT DR6, and SPT-3G D1 $TT/TE/EE$ likelihoods (${\rm CMB}_{\rm SPA}$). We further measure the structure-growth parameter $\sigma_{8}\Omega_{\rm m}^{0.25}$ to be $0.6046\pm0.0096$ from the SPT-3G D1 lensing spectrum alone and $0.6020\pm0.0084$ when combined with ACT DR6 and Planck PR4 CMB lensing. By further combining this with ${\rm CMB}_{\rm SPA}$ and the latest DESI DR2 BAO data, we obtain $\sum m_{\nu} < 0.072\,\mathrm{eV}$ (95% C.L.) when allowing the neutrino mass to vary within $\Lambda$CDM. Compared with previous work, the better agreement of our measurement with DESI DR2 BAO yields both this relaxed upper bound and reduced ($\mathord{\sim}2\sigma$) preferences for nonzero spatial curvature and for deviations of $(w_0,w_a)$ from the $\Lambda$CDM expectation. When we combine CMB lensing with the DES Y3 3$\times$2pt analysis, we obtain $S_{8}=0.811\pm0.011$, corresponding to a $1.4\%$ constraint on the late-time clustering amplitude. This precision is competitive with that obtained from the primary CMB within $\Lambda$CDM.
Read the full abstract
Solar and Stellar Astrophysics
Accretion Rates and Thermohaline Convection in Polluted White Dwarfs
Polluted white dwarfs provide a unique window into the composition of exoplanetary material, but interpreting their surface abundances requires a quantitative understanding of the mixing processes in their envelopes. Accretion of metal-rich debris onto hydrogen-rich white dwarfs drives thermohaline (fingering) convection in the underlying radiative zone. Previous studies have modeled this as enhanced diffusive transport, often inferring mass accretion rates, $\dot{M}_{\rm acc}$, that exceed constraints from X-ray observations. Here we develop a time-dependent model that treats thermohaline mixing as a propagating front whose evolution is coupled self-consistently to the surface abundance of heavy elements. We perform hydrodynamical simulations which demonstrate that thermohaline mixing drives the composition gradient towards marginal stability. Using this condition, we derive analytic solutions in shallow and deep stratification regimes, finding that the mixed layer depth, $h$, and surface abundance, $X_{\rm surf}$, both scale as $t^{1/2}$ at early times, transitioning to $h \propto t^{4/27}$ and $X_{\rm surf} \propto t^{10/27}$ at later times. The mixed layer reaches depths of only $\sim 100$ km over $10^5$ yr and, without gravitational settling, never reaches a steady state. Including settling leads to an equilibrium surface abundance scaling as $X_{\rm surf, eq} \propto \dot{M}_{\rm acc}^{30/43}$. Applied to G 29–38, the model matches the observed heavy-element abundance for $\dot{M}_{\rm acc} \sim 2\times 10^9$–$4\times 10^{9}~\mathrm{g~s^{-1}}$, in better agreement with X-ray constraints than prior estimates incorporating thermohaline mixing. We also present a general scaling relation for inferring accretion rates directly from observed surface abundances and stellar properties, giving a practical prescription for other polluted white dwarfs.
Read the full abstract
Instrumentation and Methods for Astrophysics
Feasibility of Capillary-Driven Orbital Liquid Mirror Telescopes
Space-based liquid mirror telescopes (LMTs) have been recently proposed to overcome the scaling limits of traditional rigid mirrors. In these new space architectures, capillary forces would shape a reflective liquid into a spherical optical surface in microgravity. However, the equilibrium, stability, and dynamic behavior of such liquid interfaces under orbital perturbations remain largely uncharacterized. This study investigates the feasibility of space LMTs by combining axisymmetric capillary with thin-film models. Solar radiation pressure and self-gravitation are found to be the dominant perturbations in Sun–Earth L2 orbits, causing equilibrium deformations that exceed the $\lambda/20$ optical threshold at apertures above 0.86 m for visible ($\lambda = 380$ nm) and 1.58 m for near-infrared observatories ($\lambda = 2.5$ $\mu$m). Micrometeoroid impacts self-heal within seconds to minutes, settling after filling requires years to decades at large apertures, and propulsive maneuvers breach optical limits within seconds to hours. These results highlight the need for interface control in large-scale systems such as NASA's 50 m Fluidic Telescope (FLUTE). Thermocapillary actuation supports apertures of 2.3/8.0 m (VIS/NIR) at 100 $\mu$K and 10.6/37.7 m at 10 $\mu$K. Simultaneously, parasitic Marangoni flows impose sub-mK uniformity for NIR and 1–100 $\mu$K for visible operation on any capillary LMT. Departing the capillary regime by imposing a surface-normal settling acceleration, ranging from 39 mm s$^{-2}$ (NIR,5 m) to 17 m s$^{-2}$ (VIS,50 m), is the most promising route to tens of meters in aperture. This approach enables a passive implementation effectively suppressing perturbation induced interface distortions and thermal Marangoni flows. Keywords: Liquid mirror telescope, Optics, Space
Read the full abstract
High Energy Astrophysical Phenomena · Cross-listed from gr-qc
Bekenstein–Hod Bound: A $3.3\sigma$ Confirmation from GW250114
Black holes link gravity, thermodynamics and information via limits on the relaxation rate of a perturbed system. The Bekenstein-Hod bound imposes a minimum relaxation time at fixed temperature, but its observational test demands both black hole thermodynamic characterization and decay time measurement. Here we test this bound with GW250114, the loudest gravitational-wave signal yet observed from a binary black-hole merger. We infer the remnant temperature from pre-merger data truncated at least $10\,M$ before the peak and its longest-lived decay time from post-merger data, thereby avoiding direct reuse of the same strain samples. The ringdown frequencies and damping times are allowed to vary independently rather than being fixed to the Kerr spectrum. For the primary $t_{<}=-10\,M$ analysis, the bound is verified at $3.3-3.6\sigma$ across the focal ringdown start times, representing a substantial improvement over the $91\%$ confidence level set by GW150914. The conclusion remains robust under varied pre-merger cutoffs and explicit inclusion of the short-lived first overtone in waveform modelling. This separated-data measurement converts an information-theoretic relaxation bound into a precision test of a single astrophysical black hole.
Read the full abstract
Instrumentation and Methods for Astrophysics
The Auger Radio Infill SKALA Extension (ARISE): Science Case and Instrumentation (ARENA 2026)
The Auger Radio Infill SKALA Extension (ARISE) at the Pierre Auger Observatory in Argentina was deployed in 2025 and measures cosmic-ray air showers in the energy region of the Galactic-to-extragalactic transition. ARISE is comprised of 18 SKALA-2 antennas featuring two polarization channels each, deployed within $100\,$m of a surface detector station in the enhancement area of the Pierre Auger Observatory. This area of the surface array features a denser spacing of $433\,$m between surface stations, each equipped with underground muon detectors. One of these surface detector stations provides a trigger for simultaneous readout of all ARISE antenna channels. The wide frequency range of ARISE, from $50$ to $350\,$MHz, includes the sub-band of optimum signal-to-noise ratio for air-shower radio emission against the Galactic radio background. In combination with the dense antenna spacing, this enables a relatively low detection threshold, and ARISE aims at demonstrating full detection efficiency for near-vertical air showers above $100\,$ PeV. As an advantage over the current radio detectors at Auger, which are more efficient for inclined air showers, this would enable low systematic uncertainties for physics analysis combining ARISE radio measurements with coincident measurements of the underground muon detectors in the same area. In this presentation, we will provide an overview over the ARISE instrumentation operating at the Pierre Auger Observatory and will outline the science goals.
Read the full abstract
Astrophysics of Galaxies
ALMA CO(2-1) Gas Dynamics in NGC 315: A Multi-Method Benchmark for Supermassive Black Hole Mass Measurement
We present ALMA Cycle~7 \cotwo\ observations of the circumnuclear disk in NGC~315 at an angular resolution of $0\farcs230\times0\farcs175$, improving on past measurements and resolving the sphere of influence (SOI) of the supermassive black hole (SMBH), whose mass has previously been estimated of $M_{\rm BH}= \left(2.08^{+0.33}_{-0.15}\right) \times 10^9$~M$_\odot$ The high spatial resolution and sensitivity enable robust full-cube forward modeling of the molecular gas kinematics and a direct comparison of multiple independent gas-based dynamical modeling techniques. We apply standard Bayesian codes using both MCMC and nested sampling approaches, as well as a frequentist code to the same dataset, exploring systematic uncertainties associated with the stellar mass distribution, gas surface-brightness parameterization, and disk geometry. All methods yield consistent black hole masses, indicating that the inferred $M_{\rm BH}$ is not strongly method-dependent. Combining the ensemble of independent molecular-gas-based models, we derive an ensemble median black hole mass of $M_{\rm BH}/10^9\,\mathrm{M_\odot} = 2.02^{+0.04}_{-0.05}$(stat)$^{+0.05}_{-0.04}$(sys), where the comparable contributions to the full error budget arise from modeling systematics rather than formal fitting uncertainties. Our $M_{\rm BH}$ is consistent with the empirical $M_{\rm BH}$–$\sigma_\star$ and $M_{\rm BH}$–$L_{\rm bulge}$ scaling relations, and lies 32% below an independent stellar-dynamical measurement, a discrepancy we discuss in the context of systematic differences between gas- and stellar-based methods. NGC~315 serves as a benchmark for quantifying molecular gas-dynamical $M_{\rm BH}$ systematic uncertainties and for future cross-comparisons of gaseous and stellar dynamical approaches.
Read the full abstract
High Energy Astrophysical Phenomena · Cross-listed from gr-qc
Constraining equation of state of neutron star using neutron star-black hole mergers
In this work, we study tidal effects that can be observed in gravitational wave transient from neutron star-black hole (NSBH) mergers. It has been well known that neutron star's tidal deformation has imprint of equation of state and hence composition of neutron star matter. Compared to binary neutron star mergers, where we observe a combined effect of two neutron stars by measuring the parameter known as lambda tilde. In neutron star black hole mergers, we measure the tidal deformability of the neutron star component clearly. In this work, we study the feasibility of measuring the bare tidal deformability of a neutron star with planned gravitational detectors in the near and far future. We use the parameter estimation computation framework, Bilby is used to perform Bayesian inference for hundreds of NSBH mergers in current and future GW detector networks and obtain the posterior probability distribution functions of binary parameters. We find that we need at least 20 events to clearly distinguish the equation state of neutron stars in future detectors, including the Einstein Telescope and Cosmic Explorer.
Read the full abstract
High Energy Astrophysical Phenomena
The First Array-Wide Diffuse Flux Search for UHE Neutrinos with the Askaryan Radio Array
The Askaryan Radio Array (ARA) is an ultrahigh energy (UHE) neutrino detector at the South Pole that searches for impulsive broadband radio signals from neutrino-induced particle showers in glacial ice. ARA consists of five autonomous stations with receiving antennas deployed up to 200 m deep and has accumulated the largest livetime of any in-ice radio array. We present the first array-wide diffuse UHE neutrino search using data collected from January 2013 to December 2023, incorporating improved detector characterization and simulation, including data-driven noise and electronics models, revised antenna responses, and updated neutrino interaction and lepton propagation modeling, within a unified framework for event processing, background rejection, and cut optimization across independently operating stations. This search is expected either to identify the first UHE neutrino candidates observed by an in-ice radio detector or to set the most stringent UHE neutrino diffuse flux limits above a few EeV, while informing analysis strategies for under-construction and future radio arrays, such as RNO-G and IceCube-Gen2 Radio.
Read the full abstract
Instrumentation and Methods for Astrophysics
METIS high-contrast imaging simulations: From instrument modelling to science readiness
The Mid-infrared Extremely Large Telescope (ELT) Imager and Spectrograph (METIS) instrument, expected to see first light in early 2030, aims to detect and characterise exoplanets and circumstellar disks through high- contrast imaging (HCI) and spectroscopy. The High-contrast End-to-End Performance Simulator (HEEPS), initially developed to support the design of the METIS HCI modes, has evolved into a crucial tool for the METIS science team to prepare and optimize observations. HEEPS is an open-source Python-based software with a modular architecture, integrating the wavefront Fresnel propagation package PROPER, and HCI image processing with the Vortex Image Processing (VIP) package. Though designed for METIS, its modularity has been applied to other HCI instruments as well. This work presents recent updates to HEEPS, including modelling of the final METIS pupil and Lyot stops, revised quasi-static non-common path aberrations (NCPA) and Talbot effect simulations informed by as-built optical surface errors, and updated METIS Single Conjugated Adaptive Optics (SCAO) simulations. We also discuss advancements in NCPA control strategies focusing on framerate, latency and sensing performance optimization, particularly for mitigating water vapor seeing effects using the asymmetric Lyot wavefront sensor (ALF) algorithm. With these refinements, we present a comprehensive grid of HCI performance simulations for METIS, covering a range of magnitudes in the L, M, and N-bands, and several HCI observing modes. These simulations produce updated 5-sigma sensitivity contrast curves and mock HCI observations, providing key insights on HCI performance for instrument optimization and science observation planning. Our results underscore the key role of end-to-end simulations in bridging instrumental design and scientific readiness in the ELT era.
Read the full abstract
Cosmology and Nongalactic Astrophysics
Resonant Dynamics of Gravitational Wave and Chiral Alfvén wave in the Early Universe
In this work, we investigate the resonant interaction between stochastic gravitational waves (GWs) and chiral Alfvén waves in a magnetized chiral plasma in the early Universe. Starting from covariant chiral magnetohydrodynamics coupled to linearized gravity, we derive a closed system of equations for the coupled chiral-Alfvén velocity and magnetic-field perturbations. Both analytics and numerics show a parametric resonance at the sum frequency of the two CVE-split branches, with an instability band that widens as the GW strain grows. This CVE-to-Alfvén ratio is not a free parameter. In the early Universe, it is fixed by the Standard Model's relativistic degrees of freedom. Including the plasma's backreaction on the GW makes the energy exchange nonlinear, breaking the usual single-frequency Manley-Rowe picture. For a self-consistent choice of parameters, where backreaction is only a small correction, the unsuppressed resonance still drives the system to a finite-time blow-up within a few Hubble times. We then work out the resonant frequencies and growth rates involved, mapping out what current and future gravitational-wave detectors could observe.
Read the full abstract
Earth and Planetary Astrophysics
Öpik-type collision frequency for Kozai-driven projectiles: Target bodies on eccentric and inclined orbits
For high-inclination projectiles undergoing Lidov–Kozai secular evolution, coupled variations in eccentricity and inclination complicate estimates of long-term collision frequencies with a target body. Previous collision-frequency studies for Kozai-driven projectiles have considered target-body orbits that are either circular or confined to the reference plane. The present study considers the general case in which the target body's orbit has both non-zero eccentricity and non-zero inclination. For each complete Lidov–Kozai cycle, the initial relative nodal longitude and the target body's initial argument of periapsis serve as two independent orientation variables. The single-cycle collision frequency is evaluated for each pair of initial values. Under the incommensurability condition considered here, the long-term mean is then obtained by explicitly averaging these values over both orientation variables. For the special cases in which the target body's orbit is circular or lies in the reference plane, the method is shown analytically to reduce to the corresponding formulations of previous studies. Numerical comparisons confirm these reductions. For the general eccentric and inclined case, the resulting collision frequency predicts a projectile survival curve that agrees well with the fraction of projectiles remaining in independent direct dynamical integrations. Finally, it is shown that, for a fixed dynamical setup and under the incommensurability condition, the long-term mean collision frequency associated with a given closed Hamiltonian level curve is independent of both the projectile's initial secular state along that level curve and its initial longitude of ascending node.
Read the full abstract
Solar and Stellar Astrophysics
From spectroscopic abundances to evolutionary [$\alpha$/Fe] in stellar models
Large-scale, high-resolution spectroscopic surveys map detailed chemical abundance patterns for millions of stars. However, libraries of stellar evolution models and isochrones still rely on a simplified description of chemical abundances, generally assuming either solar-scaled compositions or a constant enrichment for all $\alpha$-elements. Moreover, the definition of $\alpha$-enrichment is survey-dependent and is not necessarily the most relevant quantity for stellar structure and evolution. In the era of precision stellar astrophysics, stellar models need to better account for the wealth of chemical information provided by spectroscopic surveys. We develop a perturbative method to compute stellar evolution tracks for stars with arbitrary patterns of $\alpha$-elements. We quantify the response of stellar evolutionary tracks to variations of individual $\alpha$-elements through polynomial expansions, which can be linearly combined to produce synthetic tracks for a given abundance pattern. We also introduce the evolutionary $\alpha$-enhancement, [$\alpha$/Fe]$_{\rm ev}$, a constant enhancement for all $\alpha$-elements that reproduces the impact of an observed abundance pattern on stellar models. We validate both methods using stars with detailed APOGEE abundances and quantify the differences between spectroscopic [$\alpha$/Fe] and [$\alpha$/Fe]$_{\rm ev}$ for different surveys. We identify a minimal set of chemical elements required for accurate stellar models and present ALCHEMY, a simple and fast tool to compute [$\alpha$/Fe]$_{\rm ev}$ for general $\alpha$-element abundance patterns.
Read the full abstract
High Energy Astrophysical Phenomena
The DSA/Chronoscope fast radio burst survey: forecasts and science overview
Fast radio bursts (FRBs) are bright extragalactic transients with several mysteries surrounding their origins. Large FRB samples enable accurate measurements of the cosmic matter distribution, in particular on scales $\lesssim 10$ Mpc. These measurements will impact cosmological inference and our understanding of astrophysical feedback, from the circumgalactic medium to cluster scales. Here we forecast the expected yields, redshifts, and host galaxies of FRBs as observed by the Deep Synoptic Array (DSA), and describe the key science cases enabled by the large FRB sample. The DSA will be an interferometer consisting of 1650$\times$6.15 m antennas, operating between 0.7–2 GHz, to be located in Nevada, USA. The Chronoscope backend on the DSA, hereafter DSA/Chronoscope, is designed to search for FRBs across the field of view in real time, enabling the storage of full-polarization voltage data. Extrapolating from existing FRB surveys, we expect roughly $10^4$ FRB detections per year in each of three search sub-bands. Combining across sub-bands, the survey could produce $\sim$ 10$^5$ FRBs over the nominal 5-year DSA survey, assuming Euclidean source counts and a baseline compute backend that can search $6\times10^6$ beams at 1 ms sampling. The well-characterized DSA synthesized beam and deep simultaneous reference images will enable localization precisions of $\lesssim$ 250 milliarcseconds. Key science cases include the use of FRB propagation effects in probing cosmic baryons, and studies of the FRB phenomenon using FRB host galaxies and their local environment, as well as multiwavelength counterparts.
Read the full abstract
Cosmology and Nongalactic Astrophysics
Measuring peculiar velocity and tomographic redshift dipole with DESI DR1 catalogs
The so-called “cosmic dipole tension” challenges the Cosmological Principle by positing a discrepancy between the Solar System's peculiar velocity inferred from the Cosmic Microwave Background (CMB) dipole and that derived from large-scale structure number-count dipoles. Here we provide a high-precision determination of the kinematic dipole using the redshift-dipole method applied to the first data release (DR1) of the Dark Energy Spectroscopic Instrument (DESI). By exploiting the Doppler-induced modulation of observed redshifts, this estimator is intrinsically less sensitive to imaging systematics and selection-function uncertainties that can bias traditional number-count measurements. We conduct a tomographic analysis of four tracer populations, Bright Galaxy Sample, Luminous Red Galaxies, Emission Line Galaxies, and quasars, spanning $0.1<z<2.1$. Survey geometry and statistical uncertainties are quantified using 1,000 EZmock realizations. We find that the high-redshift QSO sample implies a peculiar velocity of $v = 357.95_{-48.47}^{+55.05}\,\mathrm{km\,s^{-1}}$, in excellent agreement with the CMB-inferred value of $369.82 \pm 0.11\,\mathrm{km\,s^{-1}}$. By contrast, a complementary number-count analysis yields a significantly enhanced dipole amplitude, which we attribute to leakage of large-scale power and to incompleteness within the DESI DR1 footprint. These results indicate that the redshift dipole provides a cleaner and more reliable probe of the kinematic rest frame, offering strong support for the standard kinematic interpretation at high redshift and helping to resolve the apparent dipole anomaly.
Read the full abstract
Instrumentation and Methods for Astrophysics
Beyond $X_\mathrm{max}$ : Reconstructing Air Shower Profiles with Information Field Theory with SKA-Low
While radio measurements of extensive air showers have shown to achieve a high precision of $X_\mathrm{max}$ sensitivity, it has been shown that parameters beyond $X_\mathrm{max}$ can also be reconstructed. These shape parameters contain additional sensitivity to the hadronic physics in the shower as well as its mass composition. In this work, we showcase a reconstruction framework to recover the full longitudinal profile from realistic radio measurements. The framework is based on Information Field Theory that infers the full profile with a forward-based model, which uses a Gaisser-Hillas profile with weakly informative shower priors, SMIET with a template library to synthesise pulses at any event geometry, and a realistic antenna response and noise level emulating that of SKA-Low. We verify the self-consistency of our framework with $\sim 900$ events generated with SMIET with antennas placed on the $\vec{v} \times (\vec{v} \times \vec{B})$ axis. The framework recovers the full profile within uncertainty and capture correlations between shower parameters. We yield an $X_\mathrm{max}$ resolution of $< 9$ g cm$^{-2}$ as well as resolutions of the width and asymmetry with minimal bias. The profile is also recovered with a bias of $< 4$% at all atmospheric depths $< 1200$ g cm$^{-2}$. We aim to apply this framework with pulses simulated from CoREAS with measured noise, ultimately extending the framework to realistic antenna layouts such as from LOFAR or SKA-Low.
Read the full abstract
Solar and Stellar Astrophysics
Extremes of solar spectral irradiance in the SORCE/XPS record
Extreme and rare changes in space mission solar irradiance records are scientifically relevant but difficult to quantify because these records are finite, instrument dependent, and affected by observational gaps and time varying measurement quality. We evaluated extreme daily logarithmic changes in the band integrated 0.1-7.0 nm irradiance measured by photodiode 7 of the Solar Radiation and Climate Experiment/X-Ray Photometer System (SORCE/XPS) from 2005 to 2019. After constructing a regular daily series by linear interpolation, we analyzed daily logarithmic changes in irradiance and fitted stationary Generalized Extreme Value models to 60 day block maxima and transformed 15 day block minima. Block lengths were selected using Ljung-Box diagnostics and sample autocorrelation functions. Maximum likelihood estimation was used as the primary inferential method, with probability weighted moments as a sensitivity check. The maximum likelihood GEV shape estimates were 0.1661 for maxima and 0.2928 for transformed minima, with closely aligned estimates under the two fitting methods. These positive point estimates are compatible with Fréchet type tails under the selected block constructions. Annualized return levels provide interpretable summaries of extreme relative increases and reductions, but estimates for long return periods remain strongly dependent on extrapolation beyond the 15 year record. Reported measurement precision and absolute uncertainty were used to qualify the interpretation of the fitted tails and were not propagated through the likelihood. By combining EVT based tail modeling with explicit consideration of measurement precision, absolute uncertainty, interpolation, and mission data gaps, the analysis provides an uncertainty aware astrostatistical baseline for extreme value inference from processed solar mission records.
Read the full abstract
Cosmology and Nongalactic Astrophysics
Cosmic variance and ergodicity in finite systems with correlations
We consider the difference between ensemble and volume average in cosmology. It is known that for sufficiently weak long-range correlations the root mean square of the difference, which we call ergodicity bias, decays like $R^{-3/2}$ in the limit of large volume $R^3$. We calculate the condition this imposes on the power spectrum of a Gaussian random field. We quantify the bias for finite $R$, and show that the $R\to\infty$ limit is of little relevance for cosmological observations when the measured scales and correlations extend to the size of the observable universe. We consider curvature, density, and velocity perturbations. On large scales the bias is important in all three cases. For the density perturbations, which are observationally the most relevant, the relative bias first exceeds 100% at the separation $r=177$ Mpc, and is larger than 100% for all $r>560$ Mpc. It should be taken into account when comparing ensemble and volume averages for large-scale structure. The bias is also large for the cosmic microwave background temperature perturbations on large angular scales, but this is not relevant for observations, as their analysis does not involve volume averaging.
Read the full abstract
Solar and Stellar Astrophysics
Study of Butterfly-Shaped Coronal Hole Evolution Across the Solar Disk
We studied the evolution of a butterfly-shaped coronal hole (CH) observed from 2025 September 8 to 14 using SDO/AIA 193~Å\ images and SDO/HMI line-of-sight magnetograms. The CH boundary was identified using a fixed 100~DN intensity threshold, and its area, mean EUV intensity, and total unsigned magnetic flux were tracked across the solar disk. As the CH approached the central meridian, its area increased from approximately $7.70\times10^{16}$ to $1.63\times10^{17}$~m$^{2}$ and its unsigned magnetic flux from $8.68\times10^{21}$ to $1.91\times10^{22}$~Mx, while the mean 193~Å\ intensity decreased from about 61 to 45~DN. The area and unsigned flux showed similar temporal trends, whereas the intensity varied oppositely. Near central-meridian passage, these quantities remained relatively stable for several days. A high-speed solar-wind stream was detected at 1~AU about 2–4 days later, reaching 600–750~km~s$^{-1}$, consistent with the CH being its likely solar source. These results link the morphological, radiative, and magnetic evolution of the CH with the subsequent high-speed solar-wind stream.
Read the full abstract
Astrophysics of Galaxies
Stellar streams around dwarf galaxies are observationally rare in the local Universe
The frequency and properties of stellar streams around dwarf galaxies remain observationally under-explored. Building a statistically significant sample is challenging because dwarf galaxies are smaller, fainter, and far more numerous than massive galaxies. We present the first results from a systematic survey for stellar streams around dwarf galaxies. We first introduce a classification metric for satellite accretion features (streams, shells, and asymmetric stellar haloes) using DECam imaging from the DES and DECaLS surveys. Applying this metric to the DES footprint, we find a low observed frequency of accretion features (5.1$\%$) and no bona fide stellar streams (with only one identified in DECaLS). We consider whether the apparent excess of shells relative to streams can be explained by observational biases. Our results indicate that stellar streams are observationally rare in the local Universe, motivating further theoretical work to determine whether this reflects genuine differences in merger rates across galaxy mass regimes.
Read the full abstract
Instrumentation and Methods for Astrophysics
Radio measurements of air showers with the IceCube-Gen2 surface prototype station at the Pierre Auger Observatory
The design of the IceCube-Gen2 observatory, proposed as a next-generation extension of IceCube, includes a surface array consisting of scintillators and radio antennas. In addition to several such detectors already deployed at IceCube's surface array at the South Pole, a complete prototype station including three SKALA antennas has been operating at the Pierre Auger Observatory for several years. This setup has been used to successfully observe radio signals from air showers, demonstrated through coincident observations with the Auger Surface Detector. In this contribution, we present an updated analysis of these radio signals, and compare them to CoREAS simulations using the reconstruction from the Auger Surface Detector as input.
Read the full abstract
Cosmology and Nongalactic Astrophysics
A consistent, physical, and analytical model for CMB observables of reionisation I. Derivation, constraining power, and detectability of the auto-spectra
The Epoch of Reionisation imprints its history and morphology on the Cosmic Microwave Background temperature and polarisation anisotropies through two effects: The kinetic Sunyaev Zel'dovich (kSZ) effect and Thomson scattering. To perform joint analyses of these observables and put tight constraints on the reionisation history and morphology, a common reionisation model is necessary. We present a quick, physically-accurate analytical approach to derive consistently the angular power spectra of the three main resulting imprints that are the spatial fluctuations of the Thomson optical depth, the patchy kSZ effect, and the scattering and screening B-modes. The approach differs from existing analytical models, as it is calibrated on high-resolution hydrodynamical simulations and depends on physical parameters of reionisation. Our predicted power spectra are consistent with the ones derived from cosmological simulations, whilst including more physical information – such as the typical size of ionised bubbles and the duration of reionisation. We provide an updated detectability assessment and show that both $C_\ell^{\tau\tau}$ and $C_\ell^{BB}$ could be measured by experiments similar to CMB Stage 4 and CMB-HD, respectively. We show the complementarity of the three observables in constraining reionisation history and morphology, and illustrate the potential of their joint analysis to get a global picture of reionisation. This analytical model will be a powerful tool in the analysis of upcoming CMB data, either alone or in combination with independent datasets such as measurements of the high-redshift 21cm signal. The code to derive these spectra is publicly available online as the preion Python package.
Read the full abstract
Astrophysics of Galaxies
Magnetohydrodynamics of charged interstellar dust. Multifluid models and study of the linear modes
Interstellar grains play key roles in star and planet formation, including the coupling of the gas to the magnetic field during the protostellar collapse. These roles depend on the local grain size distribution, which requires a multifluid treatment of charged dust. We aim to understand the fundamental physics of the dynamics of a dust distribution in interaction with the gas and the magnetic field. In particular, the purpose is to characterize the (de)coupling conditions of these different components. We provide a multifluid model of charged dust which accounts for the inertia of the grains. A chemical network is used to simulate the charge equilibrium in collapsing protostellar cores. We compute the Alfven modes and the magnetosonic modes to understand the coupling regimes between the gas, the dust fluids, and the magnetic field. We also analyze and compare to the predictions of existing models, that are the neutral dust multifluid and the standard non-ideal magnetohydrodynamics. The charged multifluid model agrees with non-ideal magnetohydrodynamics on the larger scales of a collapsing dense core and the forming disk, while we successfully extend to new regimes where the inertia of dust grains matters. We found that high charge-to-mass dust grains carry the propagation of magnetohydrodynamical waves in protostellar envelopes. We provide analytical expressions of the speed of these waves depending on the dust distribution. The magnetocompressive perturbations lead to local dust-to-gas ratio variations at au scales. A theoretical understanding of the dynamics of a charged dust distribution is provided in the linear regime. The closed set of magnetohydrodynamics equations can be implemented in numerical codes to explore nonlinear effects during the protostellar collapse such as turbulence, angular momentum transport and magnetic dust clumping.
Read the full abstract
Solar and Stellar Astrophysics
The umbrella effect: Magnetic canopy topology modulates spatially averaged chromospheric three-minute power
Three-minute ($\sim$5 mHz) oscillations are a prominent signature of wave propagation in the solar chromosphere, but their measured power can vary substantially between full-disc observations and small, high-resolution fields of view. We investigate whether chromospheric magnetic canopy topology, particularly its lateral continuity and areal filling, modulates field-of-view-averaged three-minute power. We analyse H$\alpha$ and Ca II 8542 Å spectral imaging from the Swedish 1-m Solar Telescope, comparing canopy-poor quiet-Sun scenes with canopy-dominated active-region environments, complemented by radiative-MHD Bifrost simulations of canopy-poor and canopy-rich atmospheres. We compute field-of-view-averaged power spectra from fixed-wavelength intensities and bisector-derived velocities in the observations, and from synthetic intensities, velocities, and temperature proxies in the simulations. We also quantify the relation between integrated 3-5 mHz excess power and canopy filling factor, using HMI-based magnetic-field extrapolations for the observations. In both observations and simulations, canopy-poor scenes show clear 3-5 mHz power enhancements, whereas canopy-dominated scenes show strongly reduced measured power. The trend is present in H$\alpha$ and Ca II 8542 Å intensities and is supported by complementary velocity and temperature proxies. Integrated 3-5 mHz excess power decreases with increasing canopy filling factor. We conclude that chromospheric field-of-view-averaged three-minute power is a topology-weighted observable: dense, laterally continuous magnetic canopies can strongly reduce the measured large-scale three-minute signature without requiring the absence of local wave activity. This provides a practical criterion for interpreting chromospheric oscillation power as a proxy for atmospheric wave energy in the Sun and in magnetised cool-star atmospheres.
Read the full abstract
Astrophysics of Galaxies
Constraining Baryonic Feedback at $z\sim 1$ with the H$\alpha$ Luminosity Function
Baryonic feedback remains one of the primary uncertainties in galaxy formation models because of its impact on the interstellar medium (ISM) and star formation. We investigate its effect on the H$\alpha$ luminosity function (LF) using the Illustris and IllustrisTNG simulations, incorporating direction-dependent dust attenuation and a physically motivated treatment of the propagation of Lyman-series photons. We model three sources of H$\alpha$ emission: star-forming H II regions, collisional excitation and recombination in the diffuse ISM, and AGN-induced photo-excitation of diffuse gas. The models are calibrated against observed H$\alpha$ LFs at $z\sim 1-1.5$ using a single free parameter relating the star formation rate to the H II-region H$\alpha$ luminosity. We find that emission from diffuse gas contributes $\sim 1\%-10\%$ of the total H$\alpha$ luminosity in both simulation sets. AGN photo-excitation produces only a modest enhancement of the bright end of the LF, although it dominates the H$\alpha$ emission of Illustris galaxies with $L_{\text{H}\alpha}\gtrsim 3\times10^{42}\ \rm{erg \ s^{-1}}$. The low-feedback IllustrisTNG simulations provide good agreement with the observed LFs for all emission models considered, whereas the high-feedback Illustris simulation fails to simultaneously reproduce the faint and bright ends of the LF at $z\sim 1$. This tension between H$\alpha$ LFs and stacked kSZ measurements highlights the combined power of these probes for constraining the interplay between baryonic feedback, gas properties, and the physics governing emission-line galaxies.
Read the full abstract
High Energy Astrophysical Phenomena
STOKES tables and KY codes
The family of KY codes is an XSPEC-compatible package of fully relativistic models for spectro-polarimetric modeling of X-ray binary systems and active galactic nuclei. Within Kerr space-time, they compute accretion disc line profiles, reflection spectra from discs illuminated by Comptonized coronal radiation, and thermal disc emission including disc-corona interaction and disc self-irradiation. Many of these models include polarimetric computations, and more advanced timing variants compute X-ray, UV and optical reverberation. The primary source of Comptonized X-rays is assumed to be in the lamp-post on-axis or slab above the disc geometry. The KY codes operate at various levels of complexity; some models use analytical approximations for local re-processing, while others rely on detailed numerical tables. For the most accurate treatment of rest-frame re-processing in partially ionized accretion disc atmospheres, the KY codes employ pre-computed spectro-polarimetric tables generated with the STOKES code. These transmission and reflection tables are also designed for standalone use in XSPEC and as inputs to higher-level models beyond KY. The chapter first introduces the STOKES tables in various variants and shows their non-relativistic usage in XSPEC models for (i) an extended source illuminating a distant geometrically thin disc, and (ii) a compact source illuminating a vertically extended equatorial obscurer (e.g., a torus or thick wind). The relativistic KY ray-tracing is then described, followed by an introduction to each KY model. Examples of fits to the IXPE data of X-ray binaries and active galactic nuclei using these models are shown. Together, these models form a comprehensive and publicly available framework for spectro-polarimetric modeling and data fitting of accretion. We provide links to their online repositories.