astro-ph updates on arXiv.org
astro-ph updates on the arXiv.org e-print archive.
Supermassive Black Hole Winds in X-rays -- SUBWAYS. VI. Interplay between the soft X-ray excess and ultra-fast outflows
oai:arXiv.org:2610.26936v1
arXiv:2610.26936v1 Announce Type: new
Abstract: We present a broadband, multiwavelength study of the interplay between the soft X-ray excess (SE), ultra-fast outflows (UFOs), and warm absorbers (WAs) in 21 active galactic nuclei (AGN; 0.1 <= z <= 0.4) from the SUBWAYS project. We combine simultaneous or quasi-simultaneous UV, optical, and X-ray data from XMM-Newton (EPIC pn and OM), NuSTAR, and the Hubble Space Telescope. We model the SE using a two-coronae framework, in which a warm, optically thick corona (WC) produces the optical/UV and soft X-ray emission, while a hot, optically thin corona produces the hard X-ray continuum. We also test a hybrid model that adds relativistic reflection. The optical/UV coverage constrains the Comptonization seed photon temperature to 0.1-5.5 eV. The two-coronae model consistently describes the broadband emission, with WC electron temperatures kT_e ~ 0.1-0.9 keV and photon indices $\Gamma > 2$, consistent with previous studies. However, 9/21 AGN require the hybrid model, suggesting that the SE has multiple origins whose relative contributions vary. This is especially evident in UFO hosts, where 4/7 observations favor the hybrid model. Compared with sources without UFOs, UFO hosts show lower WC temperatures, harder hot-corona spectra, and positive hardness ratios. These trends should be treated cautiously because of the small sample and their dependence on spectral parameters and instrumental response. We find no significant differences between AGN with and without WAs. Our correlation analysis finds links between the amplification factor and WC temperature, WC optical depth and hardness ratio, and WC slope and optical depth with UV luminosity. We also find a tentative anticorrelation between WC optical depth and UFO outflow velocity.
On the Nature of the BeWD System Nova FQ Circinus
oai:arXiv.org:2610.26957v1
arXiv:2610.26957v1 Announce Type: new
Abstract: The classical nova FQ Cir (Schaefer et al. 2026) has recently been shown to have a Be star donor, and has been identified as a High Mass Cataclysmic Variable (HMCV), with an orbital period of $\sim 2.04d$ (twice the 1.02d periodogram peak). Using the extensive OGLE optical archive, together with existing observations from ASAS-SN we instead find significant modulation between $\sim 50-60$~days and propose that the orbital period is likely closer to $\sim 62d$. We interpret the low amplitude 1.02 d modulation as a non-radial pulsation of the Be star. We also reassess the HMCV classification based on properties of known Be star-white dwarf systems.
Inferring the shapes of low-mass galaxies at $z>6$ from FIREbox$^{HR}$. Where are the disks?
oai:arXiv.org:2610.26969v1
arXiv:2610.26969v1 Announce Type: new
Abstract: Context. The three-dimensional shapes of galaxies provide clues to their formation and evolution. Inferring intrinsic shapes from projected sizes and axis ratios requires assumptions about the population distribution of shape parameters.
Aims. We test commonly adopted assumptions in Bayesian shape inference for high-redshift galaxies.
Methods. We analyze FIREboxHR galaxies with $\log(M_\star/M_\odot)=6.5-9.5$ at $z\sim6-15$. We measure intrinsic shapes from stellar particles using the minimum covariance determinant (MCD) estimator. We also create mock JWST/NIRCam images, measure half-light radii and axis ratios, and infer intrinsic distributions from the joint distribution of projected axis ratio $q$ and size $\log a$.
Results. Low-mass galaxies are predominantly prolate, with a substantial spheroidal component. Intrinsic ellipticity $E$ and triaxiality $T$ are well described by Johnson SB distributions, while intrinsic size $A$ is approximately log-normal. We therefore model $(\mathrm{logit}(E),\mathrm{logit}(T),\log A)$ with a multivariate normal distribution and allow general correlations. This model recovers the particle-measured shape distributions more faithfully than a multivariate normal model in $(E,T,\log A)$. Triaxiality-related parameters remain difficult to constrain from projected data without large samples or informative priors.
Conclusions. Bayesian forward modeling can recover intrinsic shape distributions when its population model adequately describes the underlying galaxies. An unsuitable distribution can bias inferences about morphological evolution at high redshift.
Euclid Bypass end-to-end simulations for spectroscopic galaxy clustering performance assessment
oai:arXiv.org:2610.26990v1
arXiv:2610.26990v1 Announce Type: new
Abstract: We present a simulation framework for Euclid galaxy clustering measurements that reproduces the essential behaviour of the spectroscopic reduction chain, bypassing the computation required by full pixel-level simulations. The framework comprises two complementary tools, FastSpec, which simulates 1-dimensional slitless spectra with a simplified instrument and noise model, and Pypelid, which computes detection probabilities on input catalogues, accounting for the survey dithering sequence and spatially varying noise. We develop calibrated prescriptions for the survey selection function, redshift errors, and redshift interloper populations. Applied to the Euclid Wide Survey, the baseline selection yields a redshift purity of 74% and a sample completeness of 45%, which improves with the use of Euclid photometric colours. The redshift error distribution is found to have a Gaussian core of width $\sigma_z = 9.4 \times 10^{-4}$, with no significant trend over $0.9 < z < 1.8$. Emission lines [OIII]5008 and [SIII]9530 misidentified as H$\alpha$ contribute a few to ten per cent depending on redshift, while catastrophic redshift measurement failures dominate at low signal-to-noise ratio. Simulating the Data Release 1 footprint, we find that complex survey selection effects can be corrected with a forward-modelled random catalogue, recovering the monopole and quadrupole of the galaxy 2-point correlation function. The simplified noise model does not include residual contamination from overlapping spectra, or other artefacts seen in early Euclid data; this will be revisited with DR1 data. The framework provides an efficient route for characterising the Euclid spectroscopic selection function at the precision required for galaxy clustering science, and is the basis for systematic error assessments in Euclid galaxy clustering analyses.
PSR~B0355+54: An Unexpected Compendium of Nearly All Scintillation Arc Morphologies
oai:arXiv.org:2610.27005v1
arXiv:2610.27005v1 Announce Type: new
Abstract: Scintillation arcs are sensitive probes of multipath propagation, offering a window into the distribution of scattering material in the interstellar medium. Contrary to the common observational finding that individual pulsars exhibit only a limited range of scintillation arc morphologies, we report that PSR~B0355+54, observed simultaneously at 2.25/8.60~GHz with the Shanghai Tianma Radio Telescope, displays nearly the full range of known arc patterns, many of which evolve rapidly. Across 54 observations from 2020 to 2026, we detect scintillation arcs in 53/53 usable observations at 2.25~GHz and 33/49 at 8.60~GHz. At 2.25~GHz, the secondary spectra exhibit five broad morphological types: forward arcs, inverted arclets, diffuse arcs, diffuse power, and irregular structures, with prominent displaced arclets appearing in several observations. In contrast, the 8.60-GHz spectra are generally dominated by faint forward arcs. Despite these marked morphological differences, the scaled curvatures and their annual variations agree between the two frequencies, showing that both trace the same dominant scattering screen at a fractional distance of $s=0.069^{+0.015}_{-0.014}$ near the pulsar, further supporting a possible association with the \ion{H}{ii} region Sh~2-205. Simulations show that this diversity can arise from variations in the scattered-image distribution within this screen, while brightness reconstructions at 2.25~GHz resolve localized structure down to $\sim0.02~\mathrm{mas}$. Tracking four arclets over 10~h at 8.60~GHz, including two displaced arclets, reveals a common evolutionary trend consistent with the line of sight sweeping across small-scale scattering structure, helping explain the rapid arc evolution.
Gamma-ray Bursts as Cosmological Probes Calibrated with Strongly Lensed Quasars and Neural Networks
oai:arXiv.org:2610.27009v1
arXiv:2610.27009v1 Announce Type: new
Abstract: Gamma-ray bursts (GRBs) extend distance measurements to redshifts beyond those routinely probed by Type Ia supernovae, but their use in cosmology is limited by calibration circularity and possible redshift-dependent effects. We investigate strongly lensed quasars as external calibrators of the Amati and Combo GRB relations. We reconstruct the luminosity-distance relation using an artificial neural network trained on TDCOSMO strong-lensing distances, with the absolute scale inherited from the lensing-only calibration. To avoid extrapolation, only GRBs within the lens-supported redshift interval are used for calibration. We then infer cosmological parameters through a hierarchical Bayesian analysis of flat $\Lambda$CDM, with $w_0w_a$CDM examined as a consistency test. For flat $\Lambda$CDM, Amati yields $H_0=78.26^{+17.87}_{-16.30}\,\mathrm{km\,s^{-1}\,Mpc^{-1}}$ and $\Omega_m=0.670^{+0.223}_{-0.276}$, while Combo yields $H_0=78.71^{+13.83}_{-11.29}\,\mathrm{km\,s^{-1}\,Mpc^{-1}}$ and $\Omega_m=0.635^{+0.235}_{-0.242}$. Amati is relatively insensitive to the phenomenological redshift term, whereas Combo shows a stronger dependence. Injection--recovery tests indicate that an unmodelled positive Combo trend can bias $H_0$ and $\Omega_m$ upward. Strong-lens calibration therefore offers a route to extending the GRB Hubble diagram, but the inferred $H_0$ is not an independent GRB measurement: it retains the TDCOSMO distance scale. The precision is limited by GRB scatter, redshift-dependent systematics, and lens-modelling uncertainties.
The science case for an in-situ measurement platform in the Venus clouds
oai:arXiv.org:2610.27038v1
arXiv:2610.27038v1 Announce Type: new
Abstract: The next step for Venus exploration after the remote sensing missions of the early 2030s will be in-situ aerial platforms, the Venus equivalent of rovers for Mars. Sophisticated Venus balloon technology is already under development and presents an opportunity for the United Kingdom to act as a bilateral partner in missions exploring key questions about the volcanic and geological history of Venus, the potential habitability of the Venus clouds, and Venus as an analogue for exoplanets. White paper written for the UKSA Space Frontiers 2035 programme.
Graph Neural Networks for Astrophysics and Cosmology
oai:arXiv.org:2610.27051v1
arXiv:2610.27051v1 Announce Type: new
Abstract: Physics seldom makes breakthroughs when we consider objects only in isolation. Instead we model and measure how objects sit among their neighbours, spatial, temporal, or otherwise: a galaxy in its environment, a halo in its history, a photon hit in a detector. Graphs are ideally suited for encoding exactly these relations, and Graph Neural Networks (GNNs) learn the geometric functions that respect them. This review presents the use of GNNs for astrophysics and cosmology, and argues that GNNs often suit physics and physical data better than any alternative: they fit the native geometry of the domain and the signal, enable the incorporation of the symmetries and prior knowledge we already trust and have toiled to build, and stay flexible enough to capture what we do not yet understand. GNNs' physically intuitive workings also render them more interpretable than other deep-learning architectures. We show how we can begin from general principles of the geometry of deep learning and build the graph and the graph neural network from the ground up, and how other architectures, like the convolutional network and the transformer, can be seen as special cases of one unified geometric blueprint. We survey applications across cosmology, galaxy formation, dynamics, and instrumentation; and we chart the frontier. From this survey, we find that only one in six articles outside instrumentation applies GNNs to data rather than just simulations. We also distil a set of best practices, including when not to reach for a graph. Throughout, we will prioritise physical intuition and foundational principles over mathematical machinery whenever possible.
The ALMA REBELS Survey: The [C II] Luminosity Function of UV-Selected Galaxies at $z\sim7$
oai:arXiv.org:2610.27052v1
arXiv:2610.27052v1 Announce Type: new
Abstract: Luminosity functions (LF) are some of the most basic quantities providing insights into galaxy evolution across cosmic time and enabling a direct comparison between observations and simulations. Here, we present the first [CII]158$\mu$m LF at $z\sim7$ from a sizeable UV-selected ALMA sample. Our analysis makes use of 27 galaxies with significant [CII] emission line detections in the epoch of reionization from the REBELS (Reionization Era Bright Emission Line Survey) ALMA Large Program and its pilot programs. We used these line detections and similar from ALPINE to calibrate UV luminosity-$L_{\mathrm{[CII]}}$ and stellar mass-$L_{\mathrm{[CII]}}$ relations at $z\sim7$. Given the challenge in precisely defining a selection volume for the REBELS sample itself, we rely on published UV LF and stellar mass function (SMF) determinations at $z\sim7$ and three complementary approaches to derive the [CII] LF. The UV LF-based method provides a lower limit on the [CII] LF, due to its missing the most dust-obscured systems. Our SMF-based estimate yields a Schechter function with a characteristic luminosity of $L^*_{\mathrm{[CII]}}=10^{9.11^{+0.11}_{-0.27}} L_\odot$ and a steep faint-end slope. From $z\sim7$ to $z\sim5$, our inferred [CII] LF shows an evolution of +0.3 dex in normalisation and +0.1 dex in $L^*_{\mathrm{[CII]}}$. Our results indicate a general underprediction from models of the abundance of [CII] emitters at high-redshift and thus potentially the number densities of massive yet UV-faint systems implied by our SMF-based method. Additionally, they suggest that 63% of the star formation rate density occurs in massive [CII]-bright systems and that $\sim$60-90% of [CII] emitters at $z\sim7$ may be missed by bright UV selections. This highlights the importance of surveys focusing on especially massive, UV-faint sources at $z>6$.
Planets surviving common envelope evolution should be rare, massive, and orbit very close to their host star
oai:arXiv.org:2610.27055v1
arXiv:2610.27055v1 Announce Type: new
Abstract: The discovery of exoplanets around white dwarfs uniquely informs our understanding of planetary survival through stellar evolution. Common Envelope (CE) evolution is one of the proposed mechanisms that could shape the orbital architecture of these planets around stellar remnants. In this work, we predict the occurrence rate and orbital separation/planet mass distributions for planets that could form around white dwarfs through CE evolution. We simulate populations of giant-planet systems, sampling initial conditions from observed mass and semi-major axis distributions for main sequence stars. We evolve the host star from main sequence to white dwarf to assess which planetary configurations could survive, and we predict the resulting post-CE planet population. Our results suggest that such planets should be rare, massive, and have very close orbital separations. We find that most of the surviving planets have masses of $>10\,M_{\rm{Jup}}$ and orbital separations $<1.4$ times the tidal disruption radius, $a_{{disrupt}}$ ($\lesssim0.006$ AU, also called the Roche limit), which is in contrast to high eccentricity migration, which produces planets with orbital separations $>2a_{disrupt}$. The properties of CE survivors from our simulations contrast with the known closely-orbiting planet WD 1856+534 b, suggesting a different mechanism may be responsible for its formation.
How does a galaxy's quenching probability depend on its black hole's mass?
oai:arXiv.org:2610.27061v1
arXiv:2610.27061v1 Announce Type: new
Abstract: We develop a parametric toy model to explore how black hole mass scatter $\sigma_{\rm BH}$ and the dependence of the quenching probability per unit time on black hole mass, parameterized as a power law with index $q$, jointly control the faint-end slope $\alpha$ of the stellar mass function of quenched central galaxies. By varying $\sigma_{\rm BH}$ and $q$, we identify two regimes. At low scatter, $\sigma_{\rm BH} \lesssim 0.4$ dex, $\alpha$ is strongly sensitive to $q$. Gradual quenching, represented by low $q$, produces shallow slopes, while sharp thresholds, with high $q$, produce steep mass functions, making $\alpha$ a powerful diagnostic of AGN feedback prescriptions. At large scatter, $\sigma_{\rm BH} \gtrsim 0.6$ dex, models with vastly different $q$ converge to similar quenched mass functions, because the intrinsic scatter erases the imprint of the feedback implementation. Comparing to cosmological hydrodynamical simulations and semianalytic models, only models with $\sigma_{\rm BH} \gtrsim 0.4$ dex match the observed faint-end slope of $\alpha \approx 0.7$. {\sc SAGE}, {\sc Dark Sage}, and {\sc COLIBRE}, with $\sigma_{\rm BH}$ of 2.0, 0.6, and 0.4 dex, respectively, reproduce this slope, whereas {\sc IllustrisTNG} and {\sc EAGLE}, with $\sigma_{\rm BH} \lesssim 0.2$ dex, produce steeper slopes. Current observational estimates of $\sigma_{\rm BH}$ allow a wide range of $q \gtrsim 0.8$, so the observed slope alone cannot identify the quenching prescription. If $\sigma_{\rm BH}$ is small, simultaneous measurements of $\alpha$ and $\sigma_{\rm BH}$ can distinguish between competing AGN feedback models.
Metal-poor molecular clouds in dwarf galaxies at parsec-scale resolution
oai:arXiv.org:2610.27064v1
arXiv:2610.27064v1 Announce Type: new
Abstract: Understanding how stars form at low metallicity is challenging because carbon monoxide (CO) emission, which usually probes the molecular gas reservoir, is faint, and tends to be concentrated in parsec-scale clumps. To overcome this difficulty, in order to study the molecular clouds in metal-poor environments, we have acquired CO(2-1) ALMA maps at 3-4 pc resolution of four nearby dwarf galaxies with nebular (O/H) metallicities ranging from roughly 1/9 Z_sun to 1/3 Z_sun. Significant CO signal in the cubes has been identified through two masks, a ``strict'', high-confidence, mask and a ``broad'' one for better completeness. The former traces the cores of the clouds, while the latter includes fainter emission in the vicinity of the cores. Through a cloud-finding algorithm applied to the two sets of masked cubes, we have identified more than 200 molecular clouds across the sample, with roughly 50 in NGC3109, at 1/9 Z_sun, the lowest metallicity parsec-scale cloud population yet observed. The best-fit size-linewidth relations for the resolved clouds are similar among all galaxies, with a slope like the Milky Way (MW), but with lower velocity dispersion at a given cloud size. For a fixed CO-to-H_2 MW conversion factor, the cloud virial ratios tend to be large, and show a weak increasing trend with increasing cloud size, and decreasing apparent surface densities and metallicity. We find evidence for faint but significant CO emission outside CO-bright cores, spanning similar velocities to those observed in the cores. We attribute this to effective shielding by gas and dust, despite the low metallicities of our sample, since the CO tends to be cospatial with dust as traced by far-infrared emission.
Small Collisions, Big Impact: The Sensitivity of Atmospheric Temperatures to Collision-Induced Absorption
oai:arXiv.org:2610.27067v1
arXiv:2610.27067v1 Announce Type: new
Abstract: Collision-induced absorption (CIA) data are relevant to determining atmospheric temperatures and, in turn, the habitability of a diverse set of planetary atmospheres. However, these data remain sparse across broad wavelength and temperature ranges. Nevertheless, multiple, sometimes disagreeing, sources of experimental and theoretical data can be available for the same CIA pair. In this paper, we assess the impact of these varying sources on atmospheric temperatures and the limitations that the available data coverage can introduce. In order to inform future measurements and calculations, we test additional cases with extrapolated wavelength and temperature ranges. We first introduce a measure of the importance of CIA opacities depending on pressure and temperature by calculating the Rosseland mean with and without each CIA pair's contribution. Secondly, we employ each CIA source in a radiative transfer model of idealised atmospheres which maximise the effect of the CIA pair of interest. Our findings highlight the uncertainties related to the current data. Temperature differences between CIA data sources range from $\sim$10\,K when sources overlap in wavelength and temperature coverage to >100\,K when spectral gaps or temperature limits differ. Spectral windows or narrow temperature coverage can also hinder the model's convergence of the temperature profile. This work aims to provide guidance on the selection of CIA opacity datasets for atmospheric modelling and to identify the wavelength and temperature regimes in which further laboratory measurements and theoretical calculations are most urgently needed. We publish all our findings in detail in the form of an accessible website.
ASTROLABE: A 3D Particle-Level Surrogate for Cosmology Simulation with Physical Flow Matching
oai:arXiv.org:2610.27068v1
arXiv:2610.27068v1 Announce Type: new
Abstract: The microphysics of dark matter (DM) is imprinted on cosmological structure across an enormous dynamic range. Nonlinear cosmological observables are most reliably modeled with N-body and hydrodynamic simulations. However, systematic N-body scans over DM theory space are prohibitively expensive. Machine learning emulators provide a possible solution, but are currently limited by either low resolution or only the ability to emulate specific observables. We introduce Astrolabe, a flow-matching-based generative emulator designed to be free of both constraints. Astrolabe acts directly on the 6N phase-space coordinates of N individual particles. With no voxel grid in the output, there is no voxel-resolution ceiling, and in principle any observable can be measured from inference. Astrolabe is computationally tractable via a field-particle hybrid transformer in which every particle cross-attends to multiscale field tokens, so compute and memory are linear in N. Astrolabe incorporates physical flow matching by taking the prior distribution for the flow to be approximate COmoving Lagrangian Acceleration (COLA)-based coordinates rather than unphysical noise, with momentum and mass conservation enforced at all points along flow time. In the Lambda-CDM cases we test, Astrolabe produces realistic fields with accurate power spectra, halo mass functions, velocity profiles, and halo density profiles. We also observe out-of-distribution generalization, with accurate emulation of dark-acoustic oscillation initial conditions. Astrolabe is also capable of being fine-tuned to simulate very large N and we demonstrate this by emulating a TNG50-3-dark simulation at a fraction of the cost.This method represents a first step toward a general and powerful cosmology simulation surrogate, which may eventually permit the inclusion of baryonic physics and a systematic exploration of richer dark sector interactions.
Beyond static torques: Self-consistent migration and dynamical trapping of low-mass planets in dusty disks
oai:arXiv.org:2610.27071v1
arXiv:2610.27071v1 Announce Type: new
Abstract: The rapid inward migration of low-mass planets remains a primary theoretical challenge in planet formation. While recent studies show that the solid component of protoplanetary disks can exert strong positive torques to counteract the negative gas torques, traditional static fixed-orbit models fail to capture the fully coupled dynamics of a migrating planet. In this work, we present 2D high-resolution, two-fluid hydrodynamical simulations using a comoving reference frame to track the self-consistent orbital evolution of embedded planets across a comprehensive six-dimensional parameter space. By allowing the planet to dynamically respond to disk-exerted torques, we identify a wide diversity of migration regimes, involving inward/outward migration and oscillatory radial confinement. We also find that massive cores can overcome dust-induced barriers, experiencing dust-mediated runaway outward or inward migration. Conversely, in dust-rich environments continuous gas-dust momentum exchange combined with planet torques dynamically regenerates co-orbital dust asymmetries, creating robust dynamical attractors where super-Earths undergo stable, oscillatory confinement, effectively halting inward drift without requiring pre-existing pressure bumps or disk substructures. This self-induced trapping mechanism could be relevant in the outer disk, beyond the water-ice line in environments with a enhanced local dust-to-gas ratio, where icy pebbles are moderately coupled to the gas. In low-metallicity or highly diffusive disks, planets generally adopt monotonic inward decay, apart from cases with outward migration. Ultimately, the diversity of these self-consistent trajectories demonstrates that static torque approximations are insufficient. Evaluating fully coupled kinematic feedback is essential for accurately predicting planetary survival and the final architecture of planetary systems.
Optical and Ultraviolet Emission from the Gas Photoionized by Gaia BH3: Archival Constraints and Predictions for the 2029 Periastron
oai:arXiv.org:2610.27075v1
arXiv:2610.27075v1 Announce Type: new
Abstract: Gaia BH3, a $32.7\,M_\odot$ black hole (BH) on an 11.6-yr, $e=0.73$ orbit around a metal-poor giant at $\approx590$ pc, accretes from its wind. Neither its mass-loss rate $\dot M_{\rm w}$ nor the fraction of captured gas reaching the BH is known. We ask whether optical and ultraviolet lines from gas photoionized by the hot flow are observable and what they reveal. Unlike X-ray and radio data, they count the ionizing photons the gas absorbs and trace its density and velocity. With Cloudy we compute the H$\alpha$ and He II 4686 emission of the wind and inflow along the orbit for $\dot M_{\rm w}=10^{-11}$-$10^{-8}\,M_\odot\,{\rm yr^{-1}}$, Eddington ratios $f_{\rm Edd}=10^{-8}$-$10^{-5}$ and wind speeds of 30 (slow) and 100 km s$^{-1}$ (fast), and test each model against archival UVES and ESPRESSO spectra. $f_{\rm Edd}$ mainly sets the lines; in dense models they come from gas falling freely onto the BH, spread over thousands of km s$^{-1}$. The archival spectra exclude $\dot M_{\rm w}>2\times10^{-9}\,M_\odot\,{\rm yr^{-1}}$ at $f_{\rm Edd}=2\times10^{-6}$ and $>2.6\times10^{-10}\,M_\odot\,{\rm yr^{-1}}$ at $10^{-5}$ (95%, slow wind), Eddington ratios already excluded by the apastron Chandra limit, $f_{\rm Edd}<7.6\times10^{-7}$. Of two simulated campaigns, only 25 epochs around the 2029 periastron go beyond existing limits. For a slow wind, whose capture rate is constant, a $\ge5\sigma$ H$\alpha$ detection would imply an accretion state 2-3 times above the 2024 limit. For a fast wind, whose captured fraction rises 14-fold toward periastron, the campaign detects $\dot M_{\rm w}\simeq10^{-8}\,M_\odot\,{\rm yr^{-1}}$ for $f_{\rm Edd}\gtrsim1.2\times10^{-6}$, which the apastron limit allows. Below the Chandra limit, a slow wind produces no detectable line.
Balmer breaks in accreting super-massive black hole binaries
oai:arXiv.org:2610.27076v1
arXiv:2610.27076v1 Announce Type: new
Abstract: We show that truncated circumbinary accretion disks around super-massive black hole binaries can produce Balmer breaks. For binaries interacting with surrounding disks in a truncated, non-accreting phase, gas accumulates at the cavity edge, and can create conditions conducive to generating spectral discontinuities from $n=2$ hydrogen bound-free absorption. We develop a simple model for predicting the presence of a Balmer edge and estimating the break strength, and we find that the strongest discontinuities are expected for binaries with periods of months to years and total masses $\gtrsim10^6\,\Msun$. These systems share considerable overlap with those binaries anticipated to appear in upcoming time-domain surveys---like LSST and the Roman HLTDS---as periodically variable quasars. We estimate that for typical feeding rates of a few tenths of their Eddington rate, a few-to-10\% of periodic binaries at $z\lesssim1$ should show prominent Balmer breaks, and that this fraction increases for the more massive black holes that can be discovered at higher redshift. We suggest that black hole binaries that truncate their accretion flows will appear as periodically variable sources, and will additionally be intrinsically X-ray weak, deficient in ionizing emission, and in favorable cases, produce notable Balmer edges in their spectra. This combination provides a complementary observational test for massive black hole binary candidates.
Localizing binary neutron stars with the Einstein Telescope: a full Bayesian comparison of the triangular and 2L configurations
oai:arXiv.org:2610.27081v1
arXiv:2610.27081v1 Announce Type: new
Abstract: The Einstein Telescope (ET), the European next-generation gravitational-wave observatory, is at a decisive stage: the choice of its detector configuration is still open, with profound implications for its scientific output. We provide a quantitative basis for this choice by comparing two candidate configurations at their proposed sites: the triangular design ($\Delta$) at the Euregio Meuse-Rhine site, and the 2L design, with interferometers in Sardinia, Italy, and Lusatia, Germany. We present the first realistic estimates of sky-localization uncertainties for binary neutron star (BNS) mergers, moving beyond Fisher-matrix forecasts to full Bayesian parameter estimation that captures the multimodal structure of the sky posterior. We focus mainly on the high-frequency interferometers alone, providing an expected scenario for the first years of ET operations, as the more challenging low-frequency instruments will come online later. BNS injections span different sky positions and luminosity distances between 40 and 2000 Mpc, corresponding to signal-to-noise ratios (SNRs) from 12 to 720. Our analysis relaxes the long-wavelength approximation and includes Earth's rotation and higher-order modes, all of which improve sky localization. We find that the $\Delta$ configuration localizes low- and moderate-SNR sources in eight degenerate sky modes spread across both hemispheres and over the day and night sides of the Earth, with 90% credible areas reaching $\sim 2 \times 10^4 \, \mathrm{deg}^2$. In contrast, the 2L configuration recovers a single dominant mode that contains at least 50% of the posterior probability and is one to two orders of magnitude smaller in area. This represents a substantial advantage for multimessenger follow-up campaigns.
Dual and lensed AGN at sub-arcsecond separations uncovered by ERIS
oai:arXiv.org:2610.27086v1
arXiv:2610.27086v1 Announce Type: new
Abstract: We present the first results of the ongoing INAF GTO program focused on the search of sub-arcsecond dual AGN, carried out with the near-infrared integral-field unit SPIFFIER of the Enhanced Resolution Imager and Spectrograph (ERIS) at the ESO Very Large Telescope, operating in adaptive-optics mode with a laser guide star. We report integral field spectroscopy with the low-resolution gratings (R ~ 5000 - 5600) of 34 multiple AGN candidates selected from the Gaia Data Release 3 source catalog via the Gaia Multi-Peak (GMP) method, which preferentially identifies luminous unobscured AGN at sub-arcsecond separations. The targets span the redshift range 0.7 < z < 3.4 and exhibit projected angular separations between ~0.2'' and ~0.9'', corresponding to physical separations of ~1.7-7.8 kpc. Of the 34 observed systems,13 (~40 %) are resolved into multiple AGN-like components; 18 (~53 %) are classified as AGN plus foreground star superpositions, one as a single AGN, and two as systems with an AGN primary and an unclassified secondary component. Among the multiple AGN systems, we confirm six dual AGN at projected separations of 0.26'' - 0.80'' (2.2 - 6.7 kpc), in the redshift range z = 0.68 - 1.57, and with velocity offsets in the range {\Delta}v ~ 290 - 950 km/s. We also identify seven gravitationally lensed quasars, only one of which was previously known, including four doubly imaged and three quadruply imaged systems. These results substantially enlarge the census of spectroscopically confirmed compact multiple AGN systems at z > 0.5 and demonstrate the efficiency of ERIS-SPIFFIER for the systematic characterization of sub-arcsecond dual and lensed AGN candidates.
DDELenS: Tomographic Weak Lensing Simulations with Dynamical Dark Energy
oai:arXiv.org:2610.27088v1
arXiv:2610.27088v1 Announce Type: new
Abstract: We present a suite of 128 cosmological simulations in $w_0w_a$CDM cosmologies with massive neutrinos, designed to enable simulation-based inference of dynamical dark energy parameters from weak gravitational lensing and other cosmological observables. Each simulation follows $1024^3$ dark matter particles in a $512\,h^{-1}$Mpc periodic box, with cosmological parameters sampled using a Sobol sequence over the ranges $\Omega_m \in [0.28, 0.35]$, $\sigma_8 \in [0.74, 0.90]$, $w_0 \in [-1.0, -0.5]$, and $w_a \in [-1.5, 0.0]$, motivated by recent results from the Dark Energy Spectroscopic Instrument survey. The nonlinear evolution is computed with $\texttt{PKDGRAV3}$ using the N-body gauge to consistently incorporate relativistic corrections from massive neutrinos, photons, and dark energy. For each cosmology, we generate 10,000 ray-traced convergence maps at six source redshifts $z_s \in \{0.5, 1.0, 1.5, 2.0, 2.5, 1100\}$ using $\texttt{LensTools}$, and additionally provide halo catalogues and projected halo overdensity fields for galaxy--lensing cross-correlation analyses. We generate 100,000 realizations at the fiducial cosmology for the construction of the covariance matrix. This simulation suite provides a foundation for field-level inference and non-Gaussian statistics with next-generation weak lensing surveys such as the Rubin Observatory LSST, Euclid, and the Nancy Grace Roman Space Telescope, as well as CMB lensing experiments including the Simons Observatory.
The SIMPLIFI Survey: Polarization of Dense Filaments with SOFIA III. Relative orientations of antenna temperature gradients with the magnetic field in the DR21 molecular cloud
oai:arXiv.org:2610.27090v1
arXiv:2610.27090v1 Announce Type: new
Abstract: We investigate magnetohydrodynamic (MHD) turbulence in the DR21 molecular cloud using far infrared polarimetry from SOFIA HAWC+ and molecular spectroscopy from the Nobeyama 45~m telescope. The Velocity Channel Gradient (VChG) method is applied to thin velocity slices of the 12CO and 13CO J=1-0 emission to determine the antenna temperature gradients and their orientations in the plane of the sky. The relative orientations of the gradients with respect to the magnetic field position angles enable an evaluation of theoretically predicted alignment of elongated, turbulent eddies for Alfvenic, fast and slow modes and self-gravity. The projected Rayleigh statistic is calculated for both the full cloud and 6 regions of interest to provide a quantitative assessment of these alignments. Our analysis finds mostly parallel alignment of the gradient direction and the local magnetic field orientation in the Main Ridge of DR21 that comprises most of the cloud mass. While parallel alignment can result from fast MHD modes or cloud collapse, previous studies show spectroscopic evidence for the large scale collapse of the DR21 cloud, which supports self-gravity as the primary cause of the observed parallel alignment. One of the six sub-filaments that are connected to the Main Ridge exhibits mostly perpendicular alignment of the position angles of the gradient and magnetic field. These perpendicular relative orientations may be a product of Alfven and slow modes in a magnetized, sub-Alfvenic medium even if the sub-filament is accreting onto the Main Ridge.
Leviathan: a very massive, extended sub-millimetre galaxy spectroscopically confirmed at $\rm z = 4.05$
oai:arXiv.org:2610.27101v1
arXiv:2610.27101v1 Announce Type: new
Abstract: We present Leviathan, a sub-millimetre galaxy (SMG) spectroscopically confirmed at $z=4.051$ with JWST/NIRSpec. Modelling 31 bands of photometry from the rest-frame ultraviolet to the millimetre, we find a stellar mass of $\log(M_\star/M_\odot)=11.34^{+0.15}_{-0.22}$, making Leviathan one of the most massive spectroscopically confirmed galaxies at $z>4$, yet consistent with $\Lambda$CDM for a star-formation efficiency of approximately 20%. Although SMGs are often thought to be building the dense cores of compact quiescent galaxies, Leviathan is extended, with a rest-frame $0.9\,\mu\mathrm{m}$ half-light radius of 4.3 kpc. Spatially resolved SED fitting gives a half-mass radius of $3.1\pm0.4$ kpc, consistent with its size in rest-frame $1.5\,\mu\mathrm{m}$ imaging. This is three times larger than quiescent galaxies of the same mass, but consistent with the size-mass relation of star-forming galaxies at $4
Slowly Rotating Fully Convective Stars at $\lesssim$ 1 Gyr: Evidence from the Hyades and Praesepe
oai:arXiv.org:2610.27112v1
arXiv:2610.27112v1 Announce Type: new
Abstract: We report the discovery of 10 stars in the Praesepe and Hyades clusters with unusually slow rotation that are likely not field contamination. Even though our search was general, all slow rotators are in the fully convective domain. We find that the fraction of slow rotators in the clusters is more than $5\sigma$ above the distribution expected from field contamination. They provide a potential link between the period distribution of the young open clusters and the field. These stars could be a distinct birth population, or they could be an indication that the onset of strong magnetic breaking is earlier than previously realized in this domain. Photometric variability and X-ray luminosity suggest the later but more data is needed to confirm. Searches for similar populations in younger clusters will help determine when slow rotators first emerge and constrain the timescale and physical mechanisms governing angular-momentum loss in fully convective stars. The discovery of these stars also highlights the importance of long-baseline photometric observations for detecting the long rotation periods of low-mass stars.
20,000 Rotation Periods in Nearby Open Clusters and Moving Groups from ZTF and ASAS-SN
oai:arXiv.org:2610.27113v1
arXiv:2610.27113v1 Announce Type: new
Abstract: Open clusters provide uniquely powerful laboratories for studying stellar evolution because their members share a common age and chemical composition. Stellar rotation of low-mass GKM dwarfs in open clusters provides direct empirical constraints on magnetic braking, with implications for stellar interiors, magnetic activity, and star--planet interactions. However, existing cluster rotation samples are dominated by space-based observations, whose limited spatial coverage and short observing baselines can make it difficult to construct an unbiased sample. Long-baseline ground-based time-domain surveys therefore provide an important complement to space-based photometry. I present a catalog of $\sim$20,000 periodic photometric signals, including stellar rotations, eclipsing-binaries, and classical pulsators, measured from the Zwicky Transient Facility (ZTF) and the All-Sky Automated Survey for Supernovae (ASAS-SN) for members of nearby open clusters and moving groups within 1.5 kpc. This catalog provides a new empirical foundation for studying the mass, chemical composition, and age dependence of stellar angular-momentum evolution and for testing models of magnetic braking across a broad range of stellar masses and evolutionary states.
Twenty-five years of XMM-Newton EPIC-MOS2 Unfocused Particle Background: Solar Cycle Modulation, Galactic Cosmic-Ray Correlation, and Implications for NewAthena
oai:arXiv.org:2610.27121v1
arXiv:2610.27121v1 Announce Type: new
Abstract: We present a study of the unfocused particle-induced background measured by the XMM-Newton EPIC-MOS2 detector, extending previous analyses to the end of 2025 and covering approximately twenty-five years of observations. Using out-of-field-of-view (outFOV) data, we investigate the long-term evolution of the particle background over two complete solar cycles and compare it with independent tracers of the Galactic cosmic-rays. We find an extremely tight correlation between the MOS2 outFOV count rate and the AMS proton flux, with an intrinsic scatter of only about 1.1\%, confirming that variations in the unfocused background are primarily driven by changes in the Galactic cosmic-rays flux. The extended temporal baseline clearly reveals the modulation associated with the solar cycle and the characteristic signatures of the 22 year Hale cycle. A Lomb-Scargle analysis identifies a dominant periodicity at 11.28 yr, while structure-function and wavelet analyses reveal a persistent $\sim$27 day modulation associated with solar rotation. We also confirm the presence of a residual background component that does not scale with the Galactic cosmic-ray flux, consistent with previous studies and possibly associated with Compton-scattered high-energy X-ray photons. This work places the particle-induced background measured by the XMM-Newton EPIC-MOS2 detector in the broader context of Galactic cosmic-ray studies, while also providing direct implications for the development of background-monitoring and mitigation strategies for future X-ray missions such as \textit{NewAthena}.
Ammonia-Hydrogen Phase Equilibria in Planetary Envelopes from First-Principles-Based Machine-Learning Molecular Dynamics
oai:arXiv.org:2610.27136v1
arXiv:2610.27136v1 Announce Type: new
Abstract: Ammonia and hydrogen are major constituents of planets within and beyond the solar system, yet how they interact in planetary envelopes is poorly known. We perform large-scale two-phase coexistence simulations of the ammonia-hydrogen system from 500 to 4000 K and 0.1 to 100 GPa, driven by machine-learning interatomic potentials trained on density functional theory calculations. From these we map the compositions of coexisting hydrogen-rich and ammonia-rich phases, and the critical curve that separates conditions of complete miscibility from those of phase separation. The critical temperature increases with pressure, from 500 K at 0.1 GPa to 2000 K near 30 GPa, and we find no phase separation at 3000 and 4000 K up to 100 GPa. The slope of the critical curve decreases substantially between 0.3 and 30 GPa, and the curve nearly flattens above about 30 GPa and 2000 K, where molecules begin to give way to atomic species. The coexistence curves are asymmetric with an ammonia-rich critical composition. For the expected envelope compositions, the phase boundary lies close to, and intersects, published thermal profiles of Uranus and Neptune, implying ammonia may be raining out in these planets today. In contrast, in warm sub-Neptunes such as TOI-270 d, ammonia and hydrogen are likely miscible in all proportions. Together with previous studies, our results suggest that silicates, water and ammonia separate from hydrogen in sequence as planets cool, with implications for their compositional, thermal and structural evolution.
Numerical shock broadening and supernova neutrino conversion: implications for DUNE
oai:arXiv.org:2610.27158v1
arXiv:2610.27158v1 Announce Type: new
Abstract: Hydrodynamical simulations spread the supernova shock front over a few grid zones, a width set by the numerical scheme rather than by the physics of the shock. Where the front crosses the H resonance, that width can make the computed MSW evolution artificially adiabatic. We quantify the effect on six three-dimensional blue-supergiant explosion models from the Garching archive by propagating individual rays twice: through the stored profile, and with the front replaced by a narrow reconstruction, 2.2 km wide, that approximates a sudden crossing while keeping the densities on either side. The stored fronts are 1.6 to 1.7 grid zones wide. Across twelve configurations, each chosen inside the energy band in which that model's front meets the resonance, the reconstruction raises the ray-averaged crossing probability $P_H$ by 0.44 to 0.76 in absolute terms, median 0.56; narrowing the reconstruction thirty-fold changes each difference by at most 0.055. Along all 12878 rays of one configuration, the sudden-limit formula evaluated on each ray's own densities agrees closely with the propagated reconstructed result. The change is largest when the density jump spans the resonance, over a band of energies that can be estimated from the presupernova density once the shock radius and the compression are known, with tails outside it. In a separate calculation with a parametrised shock and a published emission model, time-dependent conversion changes the $\nu_e$ charged-current count in a 40 kt argon detector by -4.4%, with a sign set by the flavour composition of the source. Folding the reconstructed three-dimensional fronts with a consistent emission model remains to be done.
A MUSE Study of Stellar and Ionized-Gas Kinematics in Central Group Galaxies
oai:arXiv.org:2610.27169v1
arXiv:2610.27169v1 Announce Type: new
Abstract: Understanding gas accretion and feedback in low-mass galaxy groups bridges the gap between isolated ellipticals and cluster-dominant galaxies. We present new Multi Unit Spectroscopic Explorer (MUSE) observations of 17 nearby brightest group galaxies (BGGs), revealing diverse origins for their gas reservoirs. Warm ionized gas, traced by Halpha and [N II], is detected in 77% (13/17) of the sample, with morphologies including rotating disks with or without tails and diffuse emission, rings, and filaments. Unlike brightest cluster galaxies, most BGGs host regularly rotating stars, whereas their gas displays diverse kinematics. Nearly half of the rotating gas disks are aligned with the stars, consistent with internal recycling. However, tails, diffuse emission, or rings in three aligned systems may retain signatures of past accretion or interactions. The remaining disks, ring-like structures, and most filamentary systems are misaligned, supporting gas acquisition through mergers, interactions, or halo cooling. Compact disks are consistent with centrally concentrated, low-angular-momentum gas, possibly linked to residual cooling. Most filamentary systems lie near the precipitation threshold (tcool/teddy ~ 1-2) and preferentially occupy the turbulence-dominated region of the k-plot, supporting condensation consistent with chaotic cold accretion. All galaxies with detected X-ray halos, including those without warm gas, exhibit low central entropies (<30 keV cm^2) and short central cooling times (<3 Gyr). These conditions alone cannot predict multiphase gas, suggesting that AGN- or interaction-driven perturbations and turbulence-induced thermal instabilities help regulate condensation. No single mechanism dominates: BGGs are transitional systems in which internal recycling, environmental accretion, and cooling-driven condensation coexist, linking isolated ellipticals to cluster-dominant galaxies.
The GLEAM 4-Jy (G4Jy) Sample: V. Broadband spectral analysis
oai:arXiv.org:2610.27171v1
arXiv:2610.27171v1 Announce Type: new
Abstract: The G4Jy Sample consists of 1,863 of the brightest sources in the southern sky, as selected at 151 MHz using the GaLactic and Extragalactic All-sky MWA (GLEAM) survey. We supplement existing, low-frequency radio data by completing visual inspection of radio images from the Rapid ASKAP Continuum Survey (RACS) at 887.5, 1367.5 and 1655.5 MHz. Doing so on a source-by-source basis allows us to assess the quality of the association of the radio emission, which we indicate via a flag system. Having multiple radio flux-densities allows us to construct broadband radio spectral energy distributions (SEDs), and these extend up to 20 GHz for sources that also have data from the Australia Telescope Compact Array (ATCA). We fit polynomial functions to these SEDs, and determine the distribution of sources that are best-fit by absorption models, these accounting for 11 per cent of the sample. The SEDs are further modelled using the Radio AGN in Semi-analytic Environments (RAiSE) framework, for three assumed density-profiles (set by typical halo masses of 10^12.5, 10^13.5, and 10^14.5 Msol). The derived source ages show the expected distribution in the radio-power--size (P--D) diagram, whilst derived jet powers show non-linear behaviour. The latter is believed to be the result of higher densities, associated with higher magnetic fields, restricting the expansion of the radio jets/lobes.
The Imprint of Core-collapse Asymmetries on Supernova Shock Breakout
oai:arXiv.org:2610.27183v1
arXiv:2610.27183v1 Announce Type: new
Abstract: Shock breakout is the first electromagnetic signal of a supernova, generating a highly luminous and fast-timescale burst of thermal radiation. Motivated by recent theoretical progress in producing successful three-dimensional supernova explosions of massive stars, we undertake a detailed numerical study on the impact that the core-collapse asymmetries might have on the breakout signal. We perform the first fully 3D multi-group radiation hydrodynamics simulations of shock breakout for the case of a red supergiant explosion. We find that large-scale asymmetries formed during the core-collapse process can extend the breakout flash duration to days, with the timescale set by the supernova shock fully wrapping around the progenitor star. Using detailed 3D radiative transfer post-processing, we generate synthetic multi-wavelength light curves and spectra to more closely connect to observations. The light curves show a significant viewing-angle dependence of the overall shape and peak brightness of shock breakout. We find that core-collapse asymmetries have the largest impact on the breakout light curve when observed in shorter UV wavelengths. On the other hand, at optical wavelengths, we find a diminished impact of the asymmetric breakout that more closely resembles a 1D spherical breakout. Our results show how multi-dimensional effects originating from the core-collapse mechanism can play a major role in setting the shock breakout properties, and are particularly relevant to upcoming UV space missions such as ULTRASAT and UVEX.
Evolution of Massive Main-Sequence Stars in Binary Population Synthesis. II. Consequences for Black Hole High-Mass X-ray Binaries and Binary Black Holes
oai:arXiv.org:2610.27186v1
arXiv:2610.27186v1 Announce Type: new
Abstract: Interactions and mass loss during the main-sequence (MS) phase of massive binary stars can significantly alter their subsequent evolution and the properties of the black holes (BHs) they ultimately produce. We use the rapid binary population synthesis code \textsc{COMPAS}, incorporating a new treatment of MS evolution that self-consistently tracks the convective core under arbitrary mass-loss or mass-gain histories. We explore the properties of black hole high-mass X-ray binaries (BH-HMXBs) and binary black holes (BBHs) formed under this new framework. We find that a more detailed treatment of MS evolution is necessary to reproduce the high BH masses observed in BH-HMXBs, which are otherwise absent under prescriptions that do not take the internal structure of MS stars into account. Matching the observed rates and orbital periods favors case A mass transfer that is more non-conservative than commonly assumed. We further find that mass transfer of helium-enriched material during case A mass transfer is common in massive binaries, with potential consequences for the structure of the accreting star. For the BBH population, the new MS framework alters the chirp-mass and mass-ratio distributions of detectable mergers, and case A mass transfer can make a substantial contribution to the BBHs formed through the stable mass transfer channel. These systems show a strong preference for near-equal mass ratios, broadly consistent with recent constraints from gravitational-wave observations.
Analytical Chemistry at the Planetary Frontier
oai:arXiv.org:2610.27197v1
arXiv:2610.27197v1 Announce Type: new
Abstract: Planetary science depends on chemical measurements made under conditions very different from those in a terrestrial laboratory. This review treats the analytical methods of planetary science as one connected chain, from orbital remote sensing through in-situ measurements by landers and rovers and sample collection and curation, to the micro-analysis of returned material on Earth. Orbital methods, including visible and near-infrared reflectance, thermal emission and gamma-ray and neutron spectroscopy, are discussed first. In-situ techniques are then reviewed, including laser-induced breakdown spectroscopy (LIBS), micro-focus X-ray fluorescence (XRF), deep-ultraviolet Raman and fluorescence spectroscopy, evolved gas analysis with gas chromatography-mass spectrometry (GC-MS) and miniaturized mass spectrometry. Contamination control and curation practices are described, followed by laboratory techniques such as secondary ion mass spectrometry (SIMS and NanoSIMS), atom probe tomography (APT), thermal ionization mass spectrometry (TIMS), synchrotron micro-XRF and X-ray absorption near-edge structure (XANES) spectroscopy. Accuracy, calibration and reference materials are discussed as the link between in-situ and laboratory data, and the role of chemometrics and data fusion is examined. Flight and laboratory instruments are compared in terms of spatial resolution, detection limit, sample consumption and throughput. Active neutron-gamma spectroscopy planned for the Dragonfly mission to Titan is also discussed. This comparison shows that in-situ and laboratory techniques are complementary. In-situ instruments provide context and select samples, while laboratory techniques provide isotopic, molecular and speciation data that cannot yet be obtained in situ.
Ionized Gas Kinematics and Star-forming Clumps in Extremely Metal-poor Galaxies
oai:arXiv.org:2610.27201v1
arXiv:2610.27201v1 Announce Type: new
Abstract: Extremely Metal Poor (XMP) galaxies are the least chemically evolved objects in the local Universe and provide important constraints on galaxy formation in low-metallicity environments. This work investigates the ionised gas kinematics and star formation properties of XMP galaxies to identify disturbed or multi-component gas motions and explore their connection with star formation. We analyse long-slit spectroscopic observations of nine XMP galaxies obtained with the Intermediate Dispersion Spectrograph on the 2.5m Isaac Newton Telescope. The study uses the H$\alpha$ emission line to trace ionised gas and emission-line ratios to derive gas-phase metallicities. The kinematic analysis reveals that some galaxies exhibit ordered rotation, with kinematic centres not always coincident with the luminous star-forming heads. Most systems show disturbed, asymmetric, or declining velocity profiles, while localized distortions near bright star-forming heads suggest that massive clumps influence the surrounding gas dynamics. We derive dynamical masses for the dominant clumps and investigate the mass--metallicity and size--metallicity relations. The clumps exhibit low metallicities consistent with the XMP regime and comparable to those observed in low-mass dwarf galaxies. Three galaxies show counter-rotating or kinematically decoupled ionised-gas components, indicating complex gas kinematics and possible external gas acquisition. Our sample is predominantly composed of isolated systems, with metal-poor star formation concentrated in compact head regions. These properties may reflect gas accretion through mergers, interactions, or the acquisition of metal-poor gas from the intergalactic medium and cosmic filaments. Internal processes such as outflows, gas mixing, recycling, and stochastic star formation may also contribute to the observed properties.
Are Dynamical Dark-Energy Cosmologies Old Enough?
oai:arXiv.org:2610.27209v1
arXiv:2610.27209v1 Announce Type: new
Abstract: Recent analysis combining baryon acoustic oscillation measurements with cosmic microwave background and Type Ia supernova data show an apparent preference for a time-dependent dark-energy equation of state within $w_0w_a$ parameterizations. Here we test the robustness of this result against an independent cosmological consistency requirement, namely that the expansion history must provide sufficient time to accommodate the oldest observed stellar populations. Motivated by recent stellar and globular-cluster age determinations, we impose progressively stronger lower bounds, $t_0 \geq 13.6$, $13.75$, $13.8$, and $13.85$ Gyr, on the expansion histories allowed by DESI DR2+CMB+Pantheon+ data. We find a systematic displacement of the posterior in the $(w_0,w_a)$ plane toward the $\Lambda$CDM prediction, $(w_0,w_a)=(-1,0)$, as the minimum allowed age is increased. For $t_0 \geq 13.85$ Gyr, the Mahalanobis separation from $\Lambda$CDM decreases from $2.56\sigma$ to $1.72\sigma$, bringing the $\Lambda$CDM prediction inside the 95\% credible region. Our results show that the current evidence for dynamical dark energy is sensitive to independent cosmic-age information and that stellar chronometry provides a direct robustness test of the DESI-favored expansion histories. Requiring sufficiently old cosmologies can appreciably reduce the statistical preference for departures from the standard $\Lambda$CDM model.
Large-scale Vortices in the Solar Photosphere: Statistical Properties and the Reliability of Local Correlation Tracking
oai:arXiv.org:2610.27222v1
arXiv:2610.27222v1 Announce Type: new
Abstract: Vortices are observed throughout the solar atmosphere, spanning a wide range of spatial scales. Although small-scale vortices have been studied extensively, mesoscale vortices remain less explored, and it is still unclear whether they can be reliably recovered from observations. Mesogranular-scale motions are recovered from the R2D2 simulations by applying Gaussian spatial smoothing to the horizontal velocity field with full width at half maximum (FWHM) values of 2.5-4.0 Mm. To evaluate the accuracy of local correlation tracking (LCT) and its suitability for identifying mesogranular vortices, we apply LCT to synthetic intensity maps at $\tau = 1$ using the same spatial windows and temporal averaging windows of 10-60 minutes. Vortices are identified in both the reference and LCT-derived velocity fields using the $\Gamma$ method. For all tested parameter combinations, LCT underestimates the number of mesogranular-scale vortices. The best overall recovery is obtained for relatively small spatial smoothing scales (FWHM $\leq 3.5$ Mm) and short averaging windows ($\leq 20$ minutes), which provide the most reliable vorticity and lifetime statistics while maintaining strong agreement in vortex radius and area distributions. Our analysis shows that only a modest fraction of LCT-detected vortices correspond to vortices identified in the reference velocity field (3.54%-17.88%), depending on FWHM and temporal averaging. This indicates that vortex-driven dynamics inferred from LCT-detected vortices should be interpreted with caution. Our results suggest that applying the $\Gamma$ method to LCT-derived velocity fields can recover useful statistical properties of mesoscale vortical structures, especially their temporal trends and characteristic sizes.
Searching for Dense Perturbers in Cold Dark Matter Simulations
oai:arXiv.org:2610.27236v1
arXiv:2610.27236v1 Announce Type: new
Abstract: Recently, several low-mass, dense perturbers have been inferred from astrophysical observations of gravitational lensing, stellar streams, and satellite galaxies. If interpreted as dark matter subhalos, their inner densities exceed the predictions of the standard cold dark matter (CDM) model. We systematically investigate this tension using a CDM-only zoom-in simulation of a Milky Way analog. We analyze more than 5000 subhalos with masses down to $10^6\,M_\odot$ and reconstruct their inner density profiles using a tidal truncation model. The simulated subhalos are then compared with the observed compact perturbers in terms of their density profiles, enclosed masses, concentrations, and projected cylindrical masses. None of the simulated CDM subhalos is sufficiently dense to reproduce the inferred mass distributions of the observed perturbers. We further examine the redshift evolution of subhalo concentrations and derive a fitting function to characterize this evolution. Our findings motivate further investigations into exceptionally rare formation histories that could produce such dense perturbers, as well as the possibility of physics beyond the standard CDM model.
AI and The Enclosure of the Scientific Commons
oai:arXiv.org:2610.27251v1
arXiv:2610.27251v1 Announce Type: new
Abstract: The scientific literature can be viewed as an "information commons:" an intangible resource that "no one owns and is shared by everyone." There are well-established economic and political frameworks for understanding such commons, emerging from the history of agricultural common lands found in places like preindustrial England. This article provides a primer on the history and these frameworks, and considers a key difference between the information commons and the fields of yore: while information commons are robust and immensely valuable --- some can be said to rank among humanity's greatest achievements --- the English agricultural commons no longer meaningfully exist. The latter are gone not because of the well-known "tragedy of the commons," but rather a phenomenon known as the enclosure movement (~1600--1800) that transformed the common lands into private property. Many observers, then and now, judge enclosure to have been theft at a societal scale. This article argues that artificial intelligence (AI) provides precisely the tools needed to enable a previously-infeasible enclosure of the information commons. Although AI-enabled enclosure does not operate identically to physical enclosure, it shares the key features of turning the commons into private wealth and undermining the social systems that support them (via "slop spam"). The concepts presented in this article can provide a scholarly basis for interrogating the foreboding feeling that AI may represent not just a smash-and-grab, but an attack upon a whole way of life.
J-PLUS: Identifying Outlier Spectral Energy Distributions in the Fourth Data Release Using Unsupervised Machine Learning
oai:arXiv.org:2610.27252v1
arXiv:2610.27252v1 Announce Type: new
Abstract: The selection of candidates belonging to rare or underrepresented stellar populations, as well as novel objects, is essential for spectroscopic follow-up. In this work, we demonstrate how combining broad- and narrow-band photometric information with machine-learning techniques enables this selection. Using photometric data in twelve filters from the Javalambre Photometric Local Universe Survey (J-PLUS) as input, we applied an anomaly detection technique based on autoencoders to identify spectral energy distributions (SEDs) that deviate from the dominant pattern in the J-PLUS objects sample. Subsequently, we applied the t-distributed stochastic neighbor embedding (t-SNE) algorithm to the J-PLUS colors of the outlier SEDs population. Combining information from the SIMBAD database and other astronomical surveys to classify the outliers identified by the autoencoder, we are able to characterize over-densities in the t-SNE map. Our study specifically focuses on well-defined islands in the t-SNE space, corresponding to specific populations. We identified a large sample of white dwarfs and low-mass stars, including a small subset of polluted white dwarfs, and report potential candidates for spectroscopic follow-up, including 19 carbon-rich stars, 11 cataclysmic variables and 170 white-dwarfs-main-sequence binaries, all previously unreported in the literature.
The Future of Gamma-Ray Bursts as Probes of the Epoch of Reionization in the Era of JWST and the Extremely Large Telescope
oai:arXiv.org:2610.27282v1
arXiv:2610.27282v1 Announce Type: new
Abstract: The Epoch of Reionization (EoR) is a key era of cosmological history that remains poorly understood. An examination of the Ly-alpha damping wing of high-redshift GRBs can provide direct measurements of the fraction of neutral hydrogen (i.e., neutral fraction) of the IGM in the vicinity of the GRB host galaxy, in addition to the neutral column in the host along the sight-line to the GRB. The latter quantity can help infer the fraction of ionizing photons escaping the galaxy (i.e. escape fraction); the average escape fraction being a key input to EoR models. New and upcoming observatories such as the JWST and the Extremely Large Telescope (ELT) can obtain more robust estimates of the hydrogen column density and neutral fraction around high-redshift GRBs. We simulate spectra for the JWST NIRSpec and ELT HARMONI spectrographs using different GRB magnitudes and exposure times to determine their capabilities for measuring column densities and neutral fractions with the Ly-alpha damping wings of high-redshift GRBs. We find that both spectrographs will be highly effective for obtaining column density estimates. Despite the presence of telluric lines, we find that ELT HARMONI will be a powerful instrument for retrieving accurate neutral fraction estimates, but JWST NIRSpec will be vital for obtaining neutral fraction measurements for bright GRBs that fall in redshift ranges where ELT is less effective. Additionally, shorter JWST Target of Opportunity (ToO) follow-up times would greatly increase the accuracy of neutral fraction estimates obtained with NIRSpec spectra.
The presence and order of QCD phase transitions from gravitational-wave observations of neutron star mergers
oai:arXiv.org:2610.27287v1
arXiv:2610.27287v1 Announce Type: new
Abstract: Gravitational-wave observations of neutron star mergers allows us to probe the dense portion of the nuclear equation of state. We develop a formalism for testing the existence and order of phase transitions in the cores of neutron stars using only gravitational-wave observations of their inspiral. We test our formalism on data from the first binary neutron star mergers GW170817 and GW190425, and use synthetic injections assuming the next and uber-next upgrades of the LIGO observatories, A+ and A#, as well as with Cosmic Explorer. We show that populations of mergers in either A+ or A# will be able to measure the existence of a first-order phase transition in observations lasting less than one PhD (less than half a US PhD!). On the other hand, with only the inspiral, we cannot distinguish populations of neutron stars with no phase transition from those with second- or third-order phase transitions or with first-order phase transitions with relatively large transition densities.
Globular Clusters GMRT Pulsar Survey (GCGPS) III: ORION - a multi-beam acceleration search pipeline
oai:arXiv.org:2610.27303v1
arXiv:2610.27303v1 Announce Type: new
Abstract: Globular clusters host rich millisecond pulsar populations, but the known population remains incomplete. We present the observing strategy for Phase II of the Globular Clusters GMRT Pulsar Survey (GCGPS), which uses the SPOTLIGHT system to form 160 synthesised beams simultaneously, substantially expanding the instantaneous field of view compared with the two phased-array beams of Phase I. The configuration is tunable through the beam overlap and the number of contributing antennas, allowing spatial coverage to be balanced against sensitivity. We quantify the sensitivity losses from intra-channel dispersion smearing and finite sampling, which motivate the adopted dispersion-measure limits of 50 and 150 pc/cc for Bands 3 and 4. To process the resulting data, we introduce ORION, a GPU-accelerated pipeline that combines the high-throughput Fourier-domain filtering of PULSCAN with localised parameter refinement using PRESTO. In single-beam tests with simulated binary pulsars, ORION reduces processing time by a factor of approximately 15 to 30 relative to the CPU-based Phase I pipeline, with comparable or higher folded detection significance. On-sky validation with Cycle 49 observations recovered four known globular cluster pulsars with spin periods of approximately 5.4 to 406 ms. Using eight GPUs and 64 CPU cores, ORION processes full 160-beam observations in approximately 0.25 to 0.33 times the observation duration, corresponding to three to four times real-time throughput. Together, the Phase II strategy and ORION provide a scalable framework for expanding globular cluster pulsar searches, with future development targeting highly accelerated and jerked binary systems.
Source velocities of fast radio bursts constrained by their dispersion measure variations
oai:arXiv.org:2610.27311v1
arXiv:2610.27311v1 Announce Type: new
Abstract: Growing evidence suggests that many fast radio bursts (FRBs) originate from extragalactic magnetars. Like all neutron stars, magnetars may be kicked at birth to large velocities, the values of which directly inform the physics governing stellar explosions and compact object formation. A subset of repeating FRBs show significant variability in their dispersion measures (i.e., electron column densities) over time, often exhibiting apparently secular trends. Here we show that all repeating FRBs with significant dispersion measure variations can be generically explained by motion of the source. We use existing measurements to constrain the minimum space velocities of the FRB sources, which are found to be $\gtrsim 20-2800\,{\rm km\,s}^{-1}$. We simultaneously constrain the plasma properties of the circumsource media, which are comparable to those of compact HII regions. For sources with compact radio nebulae, we directly infer the maximum source ages, and find them to be consistent with the ages of known Galactic magnetars. Our findings suggest that significant DM variability is a hallmark of young sources born with natal kicks through prompt formation channels. Applying our methods to a larger sample with continued monitoring of repeating FRB dispersion measures will reveal key physics of compact object formation.
SN 2023iop: An Extreme 1987A-like Supernova and the Death of Very Massive Stars
oai:arXiv.org:2610.27312v1
arXiv:2610.27312v1 Announce Type: new
Abstract: Supernova (SN) 1987A showed that some hydrogen-rich core-collapse supernovae (SNe) rise slowly to peak, defining a distinct subclass. Wide-field surveys have uncovered more extreme members, such as OGLE-2014-SN-73, SN 2020faa, and SN 2023vbw, that challenge neutrino-driven explosion models. We present ultraviolet-to-infrared, radio, X-ray, and spectropolarimetric observations of SN 2023iop at $z = 0.1020$, the most extreme 87A-like event to date. It rose for $\gtrsim88$ rest-frame days to $M_r \approx -20.1$ mag, also placing it among the hydrogen-rich superluminous SNe, and radiated $\gtrsim4\times10^{50}$ erg. Its spectra resemble those of typical Type II SNe and remain photospheric at $\sim330$ days post-maximum. Continuum polarization of 0.8-1.0% indicates an aspherical system. Hydrogen-rich pair-instability supernova (PISN) models of a $190~M_\odot$ star reproduce the light-curve up to 150 days after peak, but not the late-time emission. A late-time spectrum shows a feature near 7300 {\AA}, tentatively identified as [Ni II] $\lambda$7378. The implied $\approx0.02~M_\odot$ of stable $^{58}$Ni is reached only in the most massive PISN models. The low-mass host galaxy ($\sim10^{8.2}~M_\odot$) is consistent with a metal-poor progenitor as required by the PISN channel. Fallback accretion reproduces the light-curve morphology, but not the colors. We consider a pair-instability explosion the most plausible scenario, with the late-time light-curve excess powered by circumstellar interaction, which shows no signatures in the spectra. SN 2023iop shows that the most luminous 87A-like events occupy an underexplored region of massive-star evolution.
Core-Collapse Supernova Rates in the Dark Energy Survey
oai:arXiv.org:2610.27328v1
arXiv:2610.27328v1 Announce Type: new
Abstract: Over its five-year campaign (2013-2017), the Dark Energy Survey (DES) supernova (SN) survey recorded 31,636 light curves. Using a custom-built, two-stage random forest classifier, we identified 2,445 likely core-collapse SNe (CC SNe) within the survey rate analysis footprint. With this sample, we measure the volumetric CC SN rate from redshift z = 0 to 1.3 in 11 redshift bins, representing the most detailed single-survey coverage to date across nine billion years of lookback time. Over this interval, the CC SN rates are well fit by a power law R_CC(z) = R_0 (1 + z)^\alpha, with R_0 = 1.58 \pm 0.17 \times 10^{-4} \text{ yr}^{-1} \text{ Mpc}^{-3} and \alpha = 2.85 \pm 0.25. The rate turnover, expected at cosmic noon, is not yet apparent. Because CC SNe trace recent star formation, we compare our measurements with recent cosmic star-formation history (CSFH) models. Our measurements are ~70 per cent higher than CSFH models using a Salpeter initial mass function (IMF) to predict CC SN rates. Replacing the Salpeter IMF with a Baldry and Glazebrook IMF and a CC SN progenitor mass range of 8-40 M_\odot bridges the gap. Alternatively, CSFHs may miss a hitherto unidentified population of low-luminosity galaxies. We further find that the local ratio of stripped-envelope (SE) SNe to all CC SNe is 31 \pm 5 per cent, remaining constant from z = 0 to 1.3. Combined with an IMF analysis accounting for the "missing progenitor problem" of hydrogen-rich SNe, our results favour a mixture of single massive stars and binaries as SE SN progenitors, with the redshift-invariant ratio favouring binary channels.
Size, Albedo, and Taxonomy of Asteroid (1201) Strenua from Stellar Occultation and ATLAS Photometry
oai:arXiv.org:2610.27453v1
arXiv:2610.27453v1 Announce Type: new
Abstract: Context. A multi-chord stellar occultation by (1201) Strenua observed on 26 April 2026 is combined with ATLAS $c$- and $o$-band photometry.
Aims. We aim to determine the size and geometric albedo of Strenua and reassess its taxonomic classification.
Methods. Four positive chords from a 47-station campaign across Indonesia were fitted within a Bayesian forward-modeling framework to derive an ellipse cross-section and equivalent diameter. ATLAS photometry was modeled using the $H, G_1, G_2$ phase-function system to determine band-specific absolute magnitudes, converted to $H_{\mathrm{V}}$, and combined with the occultation diameter to calculate geometric albedo and update taxonomy.
Results. The chords yield semi-axes $a' = 27.49^{+3.90}_{-1.65}$~km and $b' = 14.33^{+2.37}_{-2.69}$~km (position angle $95.3^{+4.9}_{-2.7}$ deg), giving an area-equivalent diameter $D_{\mathrm{eq}} = 39.86^{+2.06}_{-1.66}$~km. Photometric modeling gives $H_c = 11.50 \pm 0.05$~mag ($G_{1,c} = 0.81 \pm 0.08$, $G_{2,c} = 0.04 \pm 0.02$) and $H_o = 11.29^{+0.02}_{-0.03}$~mag ($G_{1,o} = 0.94^{+0.03}_{-0.05}$, $G_{2,o} = 0.01 \pm 0.01$), yielding $H_{\mathrm{V}} = 11.46 \pm 0.05$~mag. The resulting geometric albedo $p_{\mathrm{V}} = 0.029 \pm 0.003$ places Strenua among the darkest primitive asteroids, securely reclassifying it as P-class and resolving its legacy X-complex ambiguity.
Exomoonlight Sonata: On the Fate of Exomoons Around Migrating White Dwarf Planets
oai:arXiv.org:2610.27484v1
arXiv:2610.27484v1 Announce Type: new
Abstract: The survival of exomoons around close-in white dwarf planets offers a potential constraint on their migration histories and an alternative pathway to post-main-sequence habitability. We investigate whether moons can survive the tidal circularization of planets via high-eccentricity migration, and consider a hierarchical three-body system (white dwarf-planet-moon) after an initial eccentricity excitation of the planetary orbit. Using an analytic equilibrium tidal model, we compare the evolution of the planetary and satellite orbits and examine the conditions required for an evection resonance. We show that the moons of pseudosynchronously rotating planets always migrate inwards since the planetary corotation radius always exceeds the Hill radius, independent of mass and orbital eccentricities. Comparing the migration timescales for the moon and the planet shows that compact moons undergo rapid tidal decay relative to the outer orbit. However, the wider separations favored by slower migration timescales experience pericenter perturbations whose work constitutes a substantial fraction of the moon's orbital energy. Together, our results indicate that moons likely do not survive during a planet's high-eccentricity migration phase in a three-body configuration.
Constraining the FRB Central Engine from Flux Variations of Persistent Radio Sources: FRB 20240114A as the First Attempt
oai:arXiv.org:2610.27489v1
arXiv:2610.27489v1 Announce Type: new
Abstract: Some repeating fast radio bursts (FRBs) are associated with persistent radio sources (PRSs), likely driven by the long-term activity of a central engine. We investigate FRB 20240114A-currently the most active repeating FRB-to explore the link between its burst activity and the flux variability of its associated PRS. We model electrons responsible for FRB emission being continuously injected from a central magnetar into its magnetar wind nebula (MWN), where synchrotron radiation produces the PRS. Comparing our model with the latest PRS spectra and flux variations, we constrain the physical properties of the FRB source. We find that the total energy injected into the MWN during FAST observations from January to August 2024 is $\sim 10^{48} \, \rm{erg}$, consistent with the magnetic energy budget of a magnetar with an internal magnetic field of $\sim 10^{16} \, \rm G$. The bulk Lorentz factor of the FRB-emitting plasma is mildly constrained to $\lesssim 10^3$, and FRB 20240114A is inferred to be powered by a young magnetar aged $\sim 10$-$20 \, \rm {yr}$. These results show that the PRS flux variability effectively probes FRB emission properties. We predict that while PRS luminosity increases rapidly with enhanced FRB activity, it will slowly fade over decades via adiabatic cooling after FRB activity ceases. Energetics suggest that the hyper-active phase observed in FRB 20240114A (burst rate $\sim 700 \, \rm hr^{-1}$) is a rare state with a duty cycle of $\lesssim 10^{-2}$, which is broadly consistent with its detection rate in current FRB-PRS samples.
A $\wedge$-like metallicity profile in the disk of the Sagittarius dwarf galaxy
oai:arXiv.org:2610.27526v1
arXiv:2610.27526v1 Announce Type: new
Abstract: Unveiling the properties of the progenitors that formed the tidal structures of the Milky Way (MW) is essential for exploring the MW's assembly history. Such studies further afford a unique opportunity to dissect the formation and evolution of dwarf galaxies by contrasting surviving and disrupted systems. However, reconstructing the intrinsic structure of progenitors solely from the positions and velocities of their tidal debris remains a formidable challenge. Here, we decode the chemo-kinematic signatures embedded in the stellar streams of the Sagittarius (Sgr) dwarf galaxy through N-body simulations. Our simulation reproduces a vertical metallicity gradient in the Sgr streams. We show that this feature is naturally explained if the stream originated from a disk component with a $\wedge$-like metallicity profile. Remarkably, similar metallicity profiles are identified in MaNGA survey data for a few (3/217) field dwarf galaxies, suggesting this feature is rare but not unique to Sgr. The vertical metallicity gradient of stellar streams preserves extra detailed galactic fossil records, enabling quantitative reconstruction of the peculiar structure of the satellite galaxy. The similar metallicity profiles of Sgr and MW raise the question of a possible special coevolution.
A Radiation-Hardened, Fault-Aware Quantized Neural Network That Self-Heals
oai:arXiv.org:2610.27621v1
arXiv:2610.27621v1 Announce Type: new
Abstract: Deep learning (DL) models are increasingly being deployed into radiation-challenging environments in spacecraft, particle accelerators, and, soon, space-based data centers. Instantiated as a learned orchestration of numbers stored in persistent memory, models are susceptible to continual degradation as high-energy particles corrupt data. At inference time, new data passed through the network fundamentally break the off-line assumptions and guarantees of a static model determined during training. To preserve the quality and efficacy of the models, fault tolerance against bit flips is a clear objective. In this work we developed a runtime method for repairing the quantized weights of neural networks (NNs) after bit flips on commodity, commercial off-the-shelf (COTS) hardware, without storing extra replicas. Our work extends the well-established Fault-Aware Training (FAT) paradigm by adding an alignment penalty to concentrate undamaged activations and weights into a low-rank subspace, and we then reconstruct flagged corrupted coordinates from the uncorrupted survivors. We call this method {\it Projection Healing} and apply it to six regression benchmarks on quantized NNs. Holding a FAT checkpoint fixed and toggling only the runtime heal gives a further ${\sim}23\%$ reduction in the native-unit excess integral. Across those datasets the method extends the operational scrub interval by a median of $2.04\times$ at a $+5\%$ RMSE degradation and $1.96\times$ at $+50\%$, compared with $1.62\times$ and $1.07\times$ for classical Single Error Correction, Double Error Detection (SECDED). This introduces a new path towards designing resilient data-driven models, complementing established hardware redundancy, ECC, and FAT-based techniques.
Probing Nearby Electron Sources and Local Diffusion with the TeV All-Electron Spectrum
oai:arXiv.org:2610.27635v1
arXiv:2610.27635v1 Announce Type: new
Abstract: Radiative losses restrict multi-TeV cosmic-ray electrons to young nearby accelerators, allowing the all-electron spectrum (electrons plus positrons) to retain information about particle release and local transport. We investigate how a nearby source can shape this spectrum using Vela as a representative example. We compare a birth-time burst with continuous release for the same total injected electron energy and examine the dependence on the diffusion coefficient, injection index, and source energetics. The burst spectrum develops a cooling-limited boundary near 28 TeV at the age of Vela, whereas continuous release produces a broader high-energy tail through the contribution of recently injected particles. The flux at Earth is largest over a finite range of diffusion coefficients, reflecting the competition between inefficient arrival and spatial dilution. Harder injection spectra and larger released energies increase the source contribution in the multi-TeV region, and we map the combinations of injection index and released electron energy that can produce an observable spectral feature. For anisotropic diffusion, the source visibility also depends on the angle between the Vela-Earth direction and the local ordered magnetic field. Future direct measurements extending the all-electron spectrum beyond several TeV could therefore distinguish different release histories and constrain both the energetics of nearby sources and local cosmic-ray transport.
Short-timescale X-ray variability of the blazar 1ES 1959+650 observed with AstroSat/SXT
oai:arXiv.org:2610.27708v1
arXiv:2610.27708v1 Announce Type: new
Abstract: We analyse the AstroSat Soft X-ray Telescope (SXT) observation of the high-synchrotron-peaked blazar 1ES 1959+650 taken on 2023 September 2, which covers 362 ks with about 95 ks of good exposure. We build 0.6-7 keV light curves in six sub-bands, binned at 5830 s (the orbital period) and 1000 s. All sub-bands show excess variance, with fractional variability of 3.7 to 6.5 per cent that increases with energy ($\propto E^{0.32\pm0.07}$). The power spectrum of the total band, fitted with a model that includes the window function of the orbital gaps, is a red-noise power law with slope $-1.80$ (95 per cent range $-2.8$ to $-1.4$). A Lomb-Scargle search finds no significant periodicity (global $p>0.57$) and would have detected, with 90 per cent probability, a sinusoid with an rms amplitude above 1.9 to 3.7 per cent at periods of 3 to 50 ks. The soft (0.6-1.4 keV) and hard (1.9-7 keV) light curves are correlated near zero lag, but only marginally (global $p=0.04$ and 0.07 for power-spectrum slopes of $-1.7$ and $-2.0$), and the lag is $-1^{+6}_{-10}$ ks, consistent with zero. No pair of sub-bands that includes 0.6-0.8 keV shows a significant correlation or lag, so an energy-dependent lag of a few ks can be neither confirmed nor excluded. The 0.6-7 keV light curve contains one hump with rise and decay times of about 0.4 to 0.6 d, which is only marginally distinguishable from red noise ($p\approx0.05$ to $0.07$). The count rate distribution is consistent with a lognormal distribution, and the hardness ratio does not depend significantly on flux.
Gaussian process-based Bayesian imaging of bolometric (sub-)millimetre data
oai:arXiv.org:2610.27733v1
arXiv:2610.27733v1 Announce Type: new
Abstract: Widefield single-dish telescopes operating at (sub-)mm wavelengths ($\lambda \sim 4$~mm to 350~$\mu$m) tend to feature much faster mapping speeds than radio interferometers at the same wavelengths, and are more suitable for measuring signals on larger spatial scales. However, the data are dominated by atmospheric fluctuations on large scales and by photon and instrumental noise on smaller scales, making it difficult to recover the true astronomical sky. We build on maria-nifty, a Gaussian process-based Bayesian imaging framework for reconstructing sky maps from single-dish telescope data, implemented on top of Numerical Information Field Theory (NIFTy). It uses modular generative models to efficiently separate the astronomical signal from atmospheric emission while providing pixel-wise uncertainty quantification for the reconstructed sky maps. Moreover, our method cleanly separates individual astrophysical signals, even though they are superposed on the sky. Following its successful validation on synthetic time-ordered data generated using the maria software, we apply it to real data observed with the MUSTANG-2 bolometric array on the 100-meter Green Bank Telescope. For the galaxy cluster RXJ1347.5-1145, we separate the extended cluster emission from the compact source at its centre and obtain a smoother map with fewer artefacts compared to a standard maximum-likelihood reconstruction using the same data. The JAX-based implementation and GPU acceleration provide the basis for applying the method to the vastly larger datasets that the next generation of single-dish telescopes will deliver.
Estimator Precision and Selection for Up-the-Ramp Sampling with a High-Read-Noise InGaAs Camera
oai:arXiv.org:2610.27742v1
arXiv:2610.27742v1 Announce Type: new
Abstract: The noise properties of an InGaAs camera vary with gain, requiring a trade-off between implementation complexity and precision. Using generalized least squares (GLS) with known covariance as the precision benchmark, we examine the applicability and precision penalties of ordinary least squares (OLS) and correlated double sampling (CDS) at different noise operating points. We define $\eta=fT/\sigma_\mathrm{R}^2$, where $f$ is the total pixel count rate, $T$ the fitting span, and $\sigma_\mathrm{R}$ the single-read noise. For constant flux, uniform sampling, and independent white read noise plus Poisson noise, established finite-sampling variance expressions give equal OLS and CDS variances at $\eta=10$ for $N\ge4$ reads; OLS is favored below this boundary and CDS above it. Analytical and Monte Carlo comparisons for three gain modes quantify the precision loss of the better simplified estimator relative to GLS. In high-gain on-sky data, choosing OLS at 20 s and CDS at 90 s increases the background residual standard deviation by 3.55% and 3.01%, respectively, relative to two-iteration feasible GLS (FGLS). Photometry of 21 common sources does not clearly resolve most estimator differences. The observed OLS-CDS crossover occurs earlier than the ideal prediction, indicating that temporal covariance beyond the ideal model also affects the selection boundary. Observational validation is limited to high gain; medium- and low-gain choices rely on theory and simulations using calibrated parameters. OLS or CDS can replace GLS with a small precision penalty at suitable operating points. The ideal criterion provides an initial selection rule, whose practical boundary requires validation against observed temporal covariance.
An improved self-gravity treatment in general relativistic smoothed particle hydrodynamics
oai:arXiv.org:2610.27775v1
arXiv:2610.27775v1 Announce Type: new
Abstract: General relativistic smoothed particle hydrodynamics (GRSPH) on a fixed black hole (BH) spacetime is suitable for fluid bodies that are much less massive than the BH, but the self-gravity of the fluid is added as an extra term in the equation of motion. Standard GRSPH implementations add this as a Newtonian source term to the momentum equation, which is no longer valid when the fluid moves at a sizeable fraction of the speed of light, close to the BH. We develop an improved prescription by deriving the weak gravitational field of each fluid element in its local rest frame and transforming it to the global frame where the momentum equation is solved. The resulting force uses rest-frame distances and softening lengths, accounts for the time dilation of the receiving particle, and its motion relative to the source particle. The derivation applies to any weakly self-gravitating flow on a fixed BH background, provided the body is small compared with its distance from the BH. Implemented in a tree-based GRSPH code, the prescription reproduces the self-gravity of a boosted polytrope with a median error below 1\% up to a Lorentz factor of 2, and keeps a binary of two polytropes stable over ten orbits. As an application, we study tidal disruption events (TDEs) of solar-type stars by supermassive BHs. The treatment of self-gravity shifts the disruption threshold and changes both the remnant core mass and the debris energy distribution. Newtonian self-gravity leaves substantial surviving cores in encounters that our prescription predicts to be nearly complete disruptions.
Cosmic-ray melt windows in interstellar ice make a four-carbon ketose but not the sugar pattern of G+0.693-0.027
oai:arXiv.org:2610.27778v1
arXiv:2610.27778v1 Announce Type: new
Abstract: Erythrulose (ERU), a four-carbon (C4) ketose, is reported towards the Galactic-centre cloud G+0.693-0.027 at a column density comparable to that of glycolaldehyde (GA), while the three-carbon (C3) sugars are not detected. We test whether melt windows opened in interstellar ice by cosmic-ray heavy ions can reproduce this pattern, with each rate fixed to a published value where available. A track model for eleven cosmic-ray nuclei determines the rate and duration of the liquid windows and the energy deposited inside and outside the melts. A radical network including C1 radicals from formaldehyde and methanol follows the inventory through repeated events. Radiolysis between events uses measured half-life doses of organics in cold ice. Hydrated GA opens a channel that preferentially forms the C4 ketose, so the melt-window chemistry can account for the ketose identity of ERU, although cold surface chemistry also yields it. The overall pattern, however, is not reproduced. The C1 radicals add one carbon at a time: dihydroxyacetone is at 15 times its limit when ERU/GA reaches the observed value. The ions deposit 0.41 eV per molecule per event outside the melts, which destroys the sugars between events; the ratio is reached only for half-life doses above 1000 eV, against measured values below 5 eV. The exposure is less certain: 2940 whole-grain melting events, 65 Myr at an ionisation rate $\zeta=10^{-14}$ s$^{-1}$. Within 7 Myr it needs $\zeta\gtrsim9.3\times10^{-14}$ s$^{-1}$, 9.3 times the dense-gas limit in the Central Molecular Zone, or $5.8\times10^{-15}$ s$^{-1}$ if the trapped radicals recombine. With standard assumptions the sugars of G+0.693 are not made in cosmic-ray melt windows. The pattern would need C1 radicals that do not add to the sugar radicals, sugars that survive the solid-phase dose and a local source of MeV ions with $\zeta\gtrsim2\times10^{-13}$ s$^{-1}$.
When the core ate the atmosphere
oai:arXiv.org:2610.27788v1
arXiv:2610.27788v1 Announce Type: new
Abstract: The initial volatile element inventories in a planet's atmosphere, mantle, and core are established by delivery and transport during planetary growth. Conversely, the volatile abundances of Earth's planetary reservoirs record its accretion pathway. We model here volatile transport during planet formation for two different delivery scenarios, via pebble and planetesimal accretion, to predict the volatile inventories of proto-Earth at the end of its main accretion phase. We find that the distributions of volatiles between the atmosphere, mantle, and core exhibit strong sensitivity to volatile delivery timing, metal-silicate partition coefficients, atmosphere-magma solubilities, and planet mass. In particular, the early volatile delivery expected in pebble accretion scenarios substantially enhances core sequestration of hydrogen and carbon relative to continuous delivery during incremental planetesimal accretion. With our fiducial partition coefficients, the pebble accretion model produces a 760 bar carbon-dominated atmosphere with 39 bar of hydrogen and water and a few bars of nitrogen, consistent with volatile concentration estimates for bulk silicate Earth including the atmosphere. The corresponding planetesimal accretion model produces a more massive ($\sim$$1350$ bar) atmosphere with similar composition. The core stores $\sim$$70$% of both the H and the C budget in the early-delivery scenario, but $\lesssim$$50$% in the corresponding continuous delivery scenario. Overall our results demonstrate that substantial amounts of volatiles are sequestered into rocky planet cores during planet formation and that the volatile inventory of bulk silicate Earth was likely established already during the magma ocean phase of planet formation.
Probing the Optical/UV Continuum Variability of Active Galactic Nuclei with the Damped Random Walk Model
oai:arXiv.org:2610.27806v1
arXiv:2610.27806v1 Announce Type: new
Abstract: Optical/UV continuum variability over baselines of up to 20 yr is investigated using the damped random walk (DRW) model for 1214 active galactic nuclei (AGNs) compiled from Y. Shen et al. and C. J. Burke et al. Among these, 115 AGNs with reliable DRW fitting and unbiased DRW parameters (timescale $\tau_{2500}$ and amplitude $\sigma_{2500}$ at the rest-frame 2500~\AA) are used to examine the relations between the variability and the accretion properties. We found a significant correlation of $\tau_{2500}$ with black hole mass ($M_{\rm BH}$) and dimensionless accretion rate ($\dot{\mathscr{M}}$) of the form $\tau_{2500} \propto M_{\rm BH}^{0.50\pm0.07} \dot{\mathscr{M}}^{0.27\pm0.04}$, while the dependence on the continuum luminosity at 5100~\AA\ and redshift becomes insignificant after controlling for $M_{\rm BH}$. $\tau_{2500}$ also correlate with the mass accretion rate ($\dot{M}$) as $\tau_{2500} \propto \dot{M}^{0.41\pm0.04}$, which is consistent with the thermal timescale expected from the standard thin-disk model within $3\sigma$ uncertainties; while the average value is consistent with the expectations of the model, the scatter in the relation can be due to differences in the viscosity parameter ($\alpha$) or in disk temperature profiles. Assuming that \tauuv\ corresponds to the accretion disk thermal timescale, the disk-emitting radius at rest-frame 2500~\AA\ ($R_{2500}$) is found to correlate with $M_{\rm BH}$ as $R_{2500} \propto M_{\rm BH}^{0.53\pm0.05}$. Using the microlensing-based $R_{2500}-M_{\rm BH}$ relation of C. W. Morgan et al., $\alpha$ is statistically estimated to be $0.04\pm0.01$.
What Shapes the Ultraviolet Spectra of Type Ia Supernovae?
oai:arXiv.org:2610.27807v1
arXiv:2610.27807v1 Announce Type: new
Abstract: Ultraviolet (UV) observations of Type Ia supernovae (SNe Ia) exhibit significant diversity, but their physical origin remains poorly understood. We investigate UV spectra of SNe Ia observed with the Swift UltraViolet and Optical Telescope. Our expanded pre-maximum sample probes dependencies on decline rate, ejecta properties, and host environment. We recover the correlation between UV flux and decline rate near maximum in the near-UV (NUV), and find correlations at pre-maximum phases in both the NUV and mid-UV (MUV). This trend is strongly wavelength dependent, with slow decliners showing higher NUV but lower MUV flux than fast decliners. Synthetic TARDIS spectra varying the 56Ni abundance show that enhanced 56Ni suppresses the MUV while increasing the NUV, broadly reproducing the trend. We find no significant correlation between UV flux and ejecta velocity, while Si II 6355 pEW shows trends similar to those with decline rate. We also find no significant correlations with host-galaxy properties, although host environment may contribute to the remaining UV diversity among SNe with similar decline rates. Early-excess SNe show enhanced flux mainly in the NUV rather than a broad UV excess, suggesting that strong circumstellar interaction is not necessarily required to explain the observed excess emission. Peculiar subclasses show distinct UV properties, with brighter and hotter subclasses showing higher UV flux, although the current sample is small and does not cover all SN Ia subtypes. Our results suggest that the UV diversity of SNe Ia is shaped by multiple mechanisms that contribute differently across the UV wavelength range.
The chemical evolution of the LMC from APOGEE abundances: a hybrid two-infall OMEGA+ analysis
oai:arXiv.org:2610.27831v1
arXiv:2610.27831v1 Announce Type: new
Abstract: The Large Magellanic Cloud (LMC) provides resolved stellar abundances and star-formation histories that can be used together to test models of dwarf-galaxy chemical evolution. We investigated if the non-monotonic APOGEE $[\mathrm{Mg}/\mathrm{Fe}]$-$[\mathrm{Fe}/\mathrm{H}]$ sequence can be explained by a hybrid history that combines a second gas-infall episode with a temporary increase in star-formation efficiency (SFE). We selected 4933 LMC red-giant-branch stars from APOGEE-2 DR17. We fitted an adaptively binned median Mg sequence with the two-zone open-box code OMEGA+. The formal fitting objective contained only the Mg sequence. The metallicity distribution function (MDF), one-dimensional abundance distributions, star-formation history (SFH), present-day mass partition, SFR, and supernova (SN) rate were retained as posterior diagnostics. The best-fitting model reproduces the Mg sequence with an RMS residual of 0.01 dex and also gives a close posterior match to O. It contains two delayed-tau infall episodes, with the second beginning at $t_2=10.39$ Gyr and with integrated mass ratio $M_2/M_1=0.159$, followed by a temporary Gaussian SFE enhancement centred at $t=11.91$ Gyr. The second infall is represented by a mixture of enriched and primordial gas. The model reproduces the timing of the late SFH enhancement and yields a present SFR within the adopted comparison range. A hybrid second-infall and SFE-enhancement history is a plausible scenario for reproducing the observed LMC Mg median sequence. Our fit explains the LMC evolution fairly well, but the remaining small MDF, mass, and SN-rate discrepancies suggest that a selection-function-aware treatment of the stellar distributions is required for a stricter comparison.
Single-Parameter Exponential Model Outperforms Standard Dark Matter and Modified Newtonian Frameworks in Predicting Galactic Rotation Curves
oai:arXiv.org:2610.27876v1
arXiv:2610.27876v1 Announce Type: new
Abstract: The functional forms of galactic rotation curves provide a fundamental test for models of gravity and dark matter. In this work, we conduct a rigorous statistical comparison of five distinct models against the complete sample of 175 galaxies from the SPARC database. The models under consideration are: (1) a one-parameter empirical law, physically constrained by each galaxy's measured flat disc rotation velocity ($V_{\text{flat}}$); (2) a classical two-parameter Freeman disk; (3) a three-parameter Navarro-Frenk-White (NFW) dark matter halo plus baryons model; (4) a one-parameter Modified Newtonian Dynamics (MOND) model; and (5) a two-parameter Bosma/Dark Matter Disk model. We employ a Bayesian framework with Markov Chain Monte Carlo (MCMC) for parameter estimation, and use the Akaike Information Criterion (AICc) to assess the goodness-of-fit and statistical parsimony of the best-fit models. Our results reveal a decisive preference for the constrained one-parameter empirical model. It is selected as the best descriptor for a clear majority of the sample (60.6\%) and, critically, also achieves the best median goodness-of-fit, with a reduced chi-squared of $\chi^2_{\nu} = 0.41$. In contrast, the standard NFW framework is preferred in only 12.6\% of cases, while the MOND and Freeman disk models are selected as optimal for just 10.3\% and 9.1\%, respectively. The Bosma model, despite linking dynamics to the gas distribution, is preferred in only 7.4\% of cases. The pronounced statistical and descriptive success of this simple, physically-anchored law over more complex, established theories suggests that its functional form represents a fundamental and highly efficient organizing principle in the dynamics of disk galaxies.
A bright X-ray flaring event from the intermediate supergiant fast X-ray transient IGR J00370+6122 observed with NuSTAR
oai:arXiv.org:2610.27925v1
arXiv:2610.27925v1 Announce Type: new
Abstract: We present an analysis of a long NuSTAR observation of the intermediate supergiant fast X-ray transient (SFXT) IGR J00370+6122, obtained around the peak of its orbital modulation. The observation has a net exposure ~93 ks. During the observation, the source displayed a bright flaring event that was resolved into numerous short flares with durations of 10-1000 s. The 3-70 keV luminosity spanned a dynamic range of ~200, from ~1E34 to ~1.8E36 erg/s; the maximum was reached for only ~40 s in total. The spectrum is well described by an absorbed double-Comptonisation model, with parameters consistent with those found in other X-ray pulsars accreting at similarly low luminosities. This is the first detection of a double-hump continuum in IGR J00370+6122, which extends this spectral state to the class of intermediate SFXTs. We searched for the previously reported ~346 s and ~674 s periodicities and find no significant signal at either period, which makes their interpretation as the neutron star spin period less likely. Using ~21 yr of Swift/BAT monitoring, we refine the orbital period of the system to P_orb=15.6654 +/- 0.0008 d. Finally, we discuss the observed flaring activity in the framework of quasi-spherical settling accretion, which is expected at the observed luminosities.
The environments of type Ibn supernovae: evidence for a strong Wolf-Rayet star component
oai:arXiv.org:2610.27931v1
arXiv:2610.27931v1 Announce Type: new
Abstract: Characterisation of progenitor stars of supernovae (SNe) poses several challenges. Direct observations of progenitors are based on serendipitious archival observations and are only possible for nearby events. Indirect methods include inferring nickel mass from the light curves or oxygen masses from the SN's nebular phase, both of which are tools to constrain the progenitor's initial mass. Alternatively, environments of SNe have been used as a proxy for the progenitor's age and metallicity. The progenitor channel of type Ibn SNe is a hotly debated topic. The mass-loss and the formation of He-rich circumstellar medium requires strong envelope stripping of the progenitor star. For core-collapse supernovae (CCSNe), this requires strong stellar wind mass losses, i.e. a Wolf-Rayet (WR) star, or binary interaction with a companion star. Here, we analyse a sample of host galaxies of type Ibn SNe that have been observed with MUSE. We use H$\alpha$ emission to infer the transient's association to star formation and derive the local metallicity with commonly used strong line diagnostic. We find that the sample of SNe are highly associated with star formation in the host galaxy, typically situated within one of the brightest HII regions in the host galaxy. The NCR distribution of H$\alpha$ pixels in SN Ibn environments strongly resembles that of WR stars, and shows a stronger preference to SF compared to other core-collapse SN types. Metallicity measurements show typically sub-solar metallicity in the SN environments.
Probing Minihalos in the Late Dark Ages through Velocity Acoustic Oscillations in the 21-cm Power Spectrum
oai:arXiv.org:2610.27933v1
arXiv:2610.27933v1 Announce Type: new
Abstract: In the Dark Ages, minihalos above the Jeans mass accrete gas from the intergalactic medium (IGM). The accreted gas would be heated to virial temperature and produces 21 cm signal in emission, while the majority of gas retained in the IGM is cold and has 21 cm signal in absorption. The 21 cm signal of the IGM generally dominates over that of minihalos, making it difficult to extract the contribution of minihalos from the detected 21 cm signal. However, in the late Dark Ages, in regions where the gas has not yet been influenced by the radiation from Pop III stars, the 21 cm signal from the IGM is weak because collisional coupling has become less efficient. More importantly, the relic supersonic dark matter-baryon relative streaming velocities suppress the gas accretion of smaller minihalos and encode their large-scale coherent features into the 21 cm signal, generating the velocity acoustic oscillation (VAO) features on the 21 cm power spectrum. This offers a promising avenue for detecting minihalos. We found that, at $z\sim 15$ ($\sim 20$), the first peak of the net VAO wiggle at $k=0.05$ Mpc$^{-1}$ is $\sim 1.4\times 10^{-3}$ ($\sim 3.1\times 10^{-4}$) mK$^2$. With a total of 24 lunar nights/8064 hr (60 lunar nights/20160 hr) of integration, the signal from $z\sim 15$ ($\sim 20$) is detectable for an array composed of 82944 (278784) dipoles, spread over $\sim 3.5$ ($\sim 6.3$) km. If the small-scale matter power spectrum is enhanced and minihalos are much more numerous than in the standard $\Lambda$CDM model, the signal is detectable for an even smaller scale array.
Constraining Solid HNCO and Its Environment in Open JWST Data with Laboratory Spectra
oai:arXiv.org:2610.27949v1
arXiv:2610.27949v1 Announce Type: new
Abstract: The presence of HNCO in the gaseous phase is well known for the interstellar medium, but the molecule is not detected yet in the solid phase. We report the first laboratory-confirmed identification of solid HNCO at 2260 cm$^{-1}$ in open JWST data using laboratory infrared spectra. For that, we grew binary mixtures of HNCO with main constituents of interstellar ice using ISEAge setup at 10 K: HNCO:H$_2$O = 1:10, HNCO:CO = 1:10, HNCO:CO$_2$ = 1:10, HNCO:CH$_3$OH~=~1:10. We present transmission infrared spectra of pure HNCO and HNCO-bearing ices along with the corresponding band strengths at 10 K. The fitting of open JWST observational data with the laboratory reference spectra in a 10--23 K temperature range show that solid HNCO is predominantly located in the apolar layer rich in CO. The abundance of solid HNCO is approximately estimated as 0.04--0.15\% relative to solid H$_2$O. The positive correlation between HNCO and N$_2$O suggests their co-formation via the chemical route that includes solid N$_2$. The results are also supported by the numerical modeling of HNCO formation in protostellar objects with predicted abundance of 0.06--0.11\%. We calculated the upper limits for CH$_3$CN toward two protostars and evaluated its abundance to solid water using the spectroscopic fitting and modeling.
VVV CL001: A strong candidate for the most ancient star cluster in the Galaxy
oai:arXiv.org:2610.27969v1
arXiv:2610.27969v1 Announce Type: new
Abstract: Finding the oldest globular cluster of the Galaxy is important to constrain a lower limit on the age of the Galaxy and the Universe. VVV CL001 ([l,b]= [5.26,0.77]) is perhaps the most metal-poor ([Fe/H]~-2.2 dex) globular cluster in the Galactic bulge and a promising candidate to be the oldest globular cluster in the entire Milky Way. Due to its recent discovery and high extinction, its physical parameters, especially the distance and age, remain highly uncertain, yet they are essential to define it as an ancient native stellar system of the Milky Way. We aim to refine the distance and age of VVV CL001 using deep and high-quality photometry, to study its role in the formation of the Milky Way, and to determine whether it is the oldest Galactic globular cluster. We used J and KS observations to construct the deepest CMD of VVV CL001 to date, based on FourStar/Baade PSF photometry combined with HST near-IR wide-filter PSF photometry. We selected cluster members through proper motion calculations and the HDBSCAN algorithm, and assigned membership probabilities via Markov-Chain Monte Carlo iterations in a 5-dimensional parameter space. Using r band data from DECam/Blanco, we optically identified horizontal branch members and applied an isochrone-fitting technique to determine the cluster distance, while members around the main-sequence turn-off yielded the age from different isochrone sets. Finally, we refined the cluster orbital parameters. We obtained a distance of d = 7.34 +- 0.46 kpc, the most precise to date. We derived an age of 13.4 +- 0.4 Gyr, indicating that VVV CL001 is one of the oldest globular clusters in the Galaxy. The orbit confirms that the cluster belongs to the Galactic bulge, which, together with its age and metallicity, could link it to the early proto-Milky Way assembly, extending the metal-poor end of the known in situ globular clusters.
Understanding the formation of alpha-Fe bimodality in the Milky Way using constrained cosmological simulations of the Local Group
oai:arXiv.org:2610.27985v1
arXiv:2610.27985v1 Announce Type: new
Abstract: We use a novel suite of constrained cosmological simulations of the Local Group to investigate the formation of alpha--Fe bimodalities in stellar discs of Milky Way analogues. In agreement with past studies using the Auriga model, we show that in our simulations the high-alpha sequence forms in a merger-induced starburst, which further quenches star formation, allowing the new stellar populations to form at low [alpha/Fe]. Producing a bimodality that is prominent at large radii, with a significant [alpha/Fe] gap and overlapping metallicity as in the Milky Way, requires the quenching of star formation across all radii where the bimodality is observed. The quenching can be achieved either with a galaxy-wide starburst, or with a retrograde merger reducing angular momentum in the outer disc and allowing more efficient dilution by the metal-poor gas in the circumgalactic medium (CGM). We show that mergers can also tilt the disc, changing the direction of galactic outflows and enriching a larger volume of the CGM with galactic metals, leading to low [alpha/Fe] scatter in the subsequently formed low-alpha sequence. We argue that the slopes in [alpha/Fe] of the two alpha sequences are crucial in narrowing down the mechanism that leads to their formation. We additionally uncover an indirect role that mergers play in enhancing the accretion of metal-poor gas onto galaxies. The smoothly accreted gas can dominate the gas content in the CGM and aid in the formation of distinct chemical sequences, especially at the low-metallicity end.
Phakoseismology: Asteroseismology of gravitationally microlensed stars
oai:arXiv.org:2610.28005v1
arXiv:2610.28005v1 Announce Type: new
Abstract: Asteroseismology relies on disk-integrated photometry or radial velocities, in which geometric cancellation suppresses modes of high spherical degree ($l\geq4$). Observations are thus limited to deeply penetrating, low-degree modes, leaving the outer envelope inaccessible. We investigate whether a high-magnification microlensing event can temporarily lift this limitation by differentially magnifying the surface of an oscillating source star. We modelled the visibility of solar-like oscillations during microlensing, weighting the surface by limb darkening and the local magnification. We computed multiplet visibilities and the power distribution among azimuthal components, and built a likelihood model for recovering the 3D spin-axis orientation from resolved rotational multiplets. When the lens-source separation becomes comparable to the source's angular radius, the magnification gradient acts as a transient spatial filter that breaks the symmetry of disk integration, raising high-degree mode visibilities by up to three orders of magnitude. The same non-axisymmetric weighting makes the relative power of azimuthal components depend on both the inclination and sky-projected azimuth of the spin axis, imprinting the full spin orientation onto the power spectrum. With $5\%$ relative-power uncertainties, we recover the orientation to $\sigma_i=2.9^\circ$ and $\sigma_\gamma=11.5^\circ$. Microlensing thus does more than amplify distant stars: its magnification pattern enables spatially selective asteroseismology. The effect is strongest for rare close encounters and unlikely to support large samples, and the newly visible high-degree modes will overlap with rotationally split and mixed modes, complicating mode identification. Nevertheless, individual high-quality events may open a window on envelope physics and spin geometries for stars in the Galactic bulge and beyond.
An improvement to the spatial discretisation in the Astrophysical Weighted Particle Hydrodynamics method
oai:arXiv.org:2610.28006v1
arXiv:2610.28006v1 Announce Type: new
Abstract: The discretisation of space in the Astrophysical Weighted Particle Hydrodynamics method is attained through a partition-of-unity-based cubature rule. The node (particle) weights in such rules must in general be obtained using approximations, for which only one has been thus far proposed in the literature. In this work, new approximations for these weights are proposed, their properties are studied, and their impact on the accuracy of the cubature rule and the overall method is assessed. Approximations for node centroids are also proposed and studied, and their possible applications are laid out. Additionally, for schemes in which the kernel-function compact-support size is set to be constant, rather than adaptive, weight and centroid approximations are provided as well. Finally, for the cubature rule arising from the weight approximation previously proposed in the literature, an asymptotic expression for its mean relative error is derived, and its usefulness for understanding the finite-sample properties of the cubature rule is showcased with some simple examples.
Kinematically isolating intracluster light in AS1063: a distinct stellar population tracing the dark matter core
oai:arXiv.org:2610.28013v1
arXiv:2610.28013v1 Announce Type: new
Abstract: We present a technique using integral-field spectroscopy to kinematically isolate intracluster light (ICL) from the brightest cluster galaxy (BCG), enabling insights into their contrasting stellar populations and the inner dark matter halo profile. Applying this technique to MUSE data of Abell S1063, we find that both components are old (~9 Gyr), but the ICL is more metal-poor than the BCG, implying distinct progenitor pathways. Within 100 kpc, the ICL comprises 55 per cent of the BCG+ICL stellar mass - substantially more than inferred from photometric decompositions. The BCG and ICL possess distinct velocity distributions, which remain nearly constant with radius, with an ICL velocity dispersion more than twice that of the BCG, and a line-of-sight velocity offset of 86 +- 18 km/s. Two-dimensional S\'ersic fits yield an ICL centroid offset of 5.2 +- 1.5 kpc from the BCG. These phase-space offsets are indicative of the BCG oscillating in a broad potential minimum. Furthermore, the ICL surface-brightness profile matches that of strong-lensing estimates assuming a cored halo. Simulations indicate that the inner halo profile varies systematically for different dark matter models. Applied to larger samples, our method therefore offers a complementary route to testing the nature of dark matter, while simultaneously revealing the stellar origins of the ICL.
Rapid Variability of Radiative Signatures from Driven Relativistic Reconnection
oai:arXiv.org:2610.28034v1
arXiv:2610.28034v1 Announce Type: new
Abstract: Rapid variability in energetic flares is a distinctive feature of high-energy emission from many magnetically dominated astrophysical sources. We investigate the radiative variability of driven relativistic magnetic reconnection using Particle-in-Cell simulations of strongly magnetized flux tubes pushed to merge at driving speed $v_{\rm push}$. This setup captures the onset of reconnection and provides a controlled model for single reconnection-driven flares near compact objects, where magnetic structures can be twisted, sheared, and forced to merge by turbulent accretion flows or relativistic frame dragging. We find that the merger efficiently accelerates particles into nonthermal distributions, producing synchrotron emission whose variability tracks the reconnection phases. We identify two regimes: a quasi-steady, weakly driven regime for $v_{push}\leq0.1\,c$, and strongly driven regime for $v_{\rm push}>0.1\,c$. The transition occurs when the driving exceeds the typical inflow-to-outflow velocity ratio of $\sim0.1$ in steady collisionless relativistic reconnection, i.e., the average rate at which reconnection can process incoming flux, potentially providing an additional source of variability. The lightcurves rise on $1$--$3$ light-crossing times of the reconnecting region, and the particle spectra develop power laws with index $p\simeq1.8$--$2.4$. These timescales are compatible with observed variability in high-energy sources, including the very-high-energy flares of M87$^*$ and the $\gamma$-ray flares of the Crab Nebula, making this model a viable framework for rapid flaring in active galactic nuclei and pulsar wind nebulae.
On the superbursts as photodisintegrators of ashes from the rp-process -- Numerical simulations of superbursts with MESA II
oai:arXiv.org:2610.28049v1
arXiv:2610.28049v1 Announce Type: new
Abstract: Context: Superbursts are energetic explosions in the crusts of neutron stars in Low-Mass X-ray Binaries (LMXBs). At temperatures above 4 GK, photodisintegration of heavy rp-process ashes may play a major role alongside $^{12}$C fusion. Aims: We investigate how photodisintegration of rp-process ashes affects superburst properties and the initial $^{12}$C abundance required for ignition. Methods: We simulate superbursts with MESA v.24.08.1. Besides the canonical fuel of $20\%$ $^{12}$C and $80\%$ $^{56}$Fe, we test different $^{12}$C abundances mixed with scaled rp-process ashes. Results: Superbursts can ignite with only $10\%$ $^{12}$C when the remaining $90\%$ consists of rp-process ashes. Their luminosity is comparable to that of the canonical fuel, whereas $10\%$ $^{12}$C mixed with $90\%$ $^{56}$Fe yields lower luminosity. $^{56}$Fe remains the most abundant ash nuclide. Although photodisintegration drives nucleosynthesis, neutron, proton and alpha captures release energy at late times. We also find that $A=64$ nuclides or $^{20}$Ne help reduce the $^{12}$C abundance needed for ignition. Conclusions: The composition of accreted material - not just its $^{12}$C content - substantially affects superburst ignition, luminosity and nucleosynthesis.
Constraints on Dark Matter Annihilation from Starlight-Induced Inverse Compton Emission in Omega Centauri
oai:arXiv.org:2610.28068v1
arXiv:2610.28068v1 Announce Type: new
Abstract: The search for the particle identity of dark matter (DM) remains one of the foremost challenges in modern physics. In this work, we derive new constraints on the DM annihilation cross section $\langle\sigma v\rangle$ by considering the inverse Compton (IC) scattering of DM annihilation injected $e^+e^-$ off the starlight radiation field in Omega Centauri. As a likely remnant core of a dwarf galaxy, Omega Centauri provides a dense stellar environment that enhances the predict IC signal. Using multi-wavelength observations from Chandra and Fermi-LAT, we establish upper limits on $\langle\sigma v\rangle$ over a wide DM mass range from MeV to TeV. Notably, in the tens-of-MeV range, our results surpass existing constraints by approximately two orders of magnitude, exceeding $\langle\sigma v\rangle \sim 10^{-31}~{\rm cm^3\,s^{-1}}$. At several GeV to the TeV DM masses range, our constraints also remain stronger than existing limits. Our results demonstrate that starlight induced IC emission is an important component for indirect DM searches in dense stellar systems, establishing Omega Centauri as a powerful probe of the nature of DM.
Pulsar J1909+0749 and Supernova Remnant G42.1-0.2: A case of a pulsar that has exited the shell of its associated supernova remnant
oai:arXiv.org:2610.28073v1
arXiv:2610.28073v1 Announce Type: new
Abstract: (Abridged) The ATNF pulsar catalogue features only 63 pulsars that are possibly associated with a supernova remnant (SNR). Even rarer is the scenario in which the pulsar has escaped the physical extent of its corresponding SNR. We aim to identify the serendipitously discovered shell-type object G42.1-0.2 as a SNR. The known pulsar (PSR) J1909+0749 lies well outside the shell of SNR G42.1-0.2, but it has a cometary tail pointing towards the centre of the SNR; we propose that the SNR and the pulsar originated from the same supernova explosion. By measuring the flux densities of the shell G42.1-0.2 from ${\sim}$200 MHz to ${\sim}$6 GHz, we establish the nonthermal nature of its emission. With the characteristic age of the pulsar, its distance from the Sun estimated from the dispersion measure (DM), and the angular separation between the pulsar and the centre of the shell, we infer the evolutionary state of the SNR and the transverse velocity of the pulsar. Using proper motion measurements, we compare the SNR-pulsar direction with the velocity of the pulsar. The radio spectral index of the shell is $-0.36 \pm 0.08$, confirming that the emission is primarily nonthermal synchrotron. The proper motion measurements are not well constrained, but the amplitude and the direction are consistent with the interpretation that the SNR and the pulsar are associated. If the characteristic age and the DM-based distance represent the true age and the true distance of the pulsar and the SNR, then the standard evolution equations suggest that the SNR must still be in the Sedov-Taylor phase, evolving in a medium with a density of about 1.5 cm$^{-3}$. In addition, as a result of the PSR-SNR association, the transverse velocity of the pulsar PSR J1909+0749 is estimated to be ${\sim}1460~\mathrm{km\,s}^{-1}$, ranking it among the highest transverse velocities observed for pulsars in the Milky Way.
The T-matrix method for modeling the optical properties of grains in exoplanetary atmospheres
oai:arXiv.org:2610.28117v1
arXiv:2610.28117v1 Announce Type: new
Abstract: Observations of exoplanetary atmospheres have revealed a rather complex co-existence between the gaseous and condensed phases, ranging from atmospheres apparently covered in clouds to clear ones, across a myriad of intermediate regimes. The interaction between solids and stellar radiation is of fundamental importance when forward or retrieval models are compared to observations. The simplest approach is to assume spherical particles, whose scattering and absorption properties are readily computed with Mie theory. More advanced computational methods, however, can also treat non-spherical condensates. Among these, the transition matrix approach has proven to provide a rigorous and efficient calculation of the extinction by non-spherical particles over a broad range of sizes and refractive indices. We here focus on the absorption and scattering cross sections of forsterite, spinel and silicon monoxide from the near to the far infrared range for particles of different sizes, with shapes ranging from compact ensembles of spheres to porous aggregates. We find that, depending on the size and the configuration chosen, the differences in extinction can be as high as a factor of 2-3, with structured grains being more efficient absorbers in the near infrared for sizes greater than 1 $\mu$m. In a fiducial transmission spectrum we show that in the $\lambda = 1-5\,\mu$m range the amount of spherical material needed to match the contribution of the structured grains can be a factor 2 higher. Such calculations, available for the community through a Zenodo repository, prove to be useful for interpreting recent and future observations of exoplanetary atmospheres.
Deriving Interplanetary Sector Boundary Crossings using Potential Field Source Surface Modeling
oai:arXiv.org:2610.28123v1
arXiv:2610.28123v1 Announce Type: new
Abstract: The Heliospheric Current Sheet (HCS) and associated interplanetary sector boundaries fundamentally structure the solar wind environment, modulating coronal mass ejection (CME) evolution and solar energetic particle transport. In this study, we evaluate the capability of a static Potential Field Source Surface (PFSS; $R_{ss}=2.5R_\odot$) coronal model to directly map interplanetary sector boundary crossings from the Sun to 1 au without dynamic heliospheric propagation models. We compare ballistically transported HCS crossings derived from GONG-ADAPT and SDO/HMI magnetic maps against sector boundary crossing observations from the SBCross catalog. PFSS modeling reproduces the large-scale interplanetary magnetic field topology well ($J\approx0.50$), achieving hit rates (POD) of 0.64 for GONG-ADAPT and 0.76 for HMI. Carrington synoptic maps outperform flux-transport modeled maps, suggesting that the short-term observational noise from data assimilation degrades topological matching more than map "aging". Predictive skill shows a strong solar cycle dependence, performing significantly better during solar maxima (POD=0.70-0.81) than solar minima (POD=0.56-0.67), attributed to more distinct crossings due to higher HCS inclination angles during maximum. Matched events exhibit travel times close to ballistic estimates of 4.8-4.9 days with an average timing difference centered near zero ($\pm1.1$ days) and 90% of events within $\pm 1.8$ days. The hit rate of HCS structures that interact with CMEs does not degrade (POD = 0.67-0.83), despite propagating $\approx7$ hours faster on average, implying that the CME interaction causes only minor structural changes to the HCS. Larger timing differences (>24 hours) are associated with flatter current sheets. Overall, static PFSS extrapolations offer a computationally inexpensive, effective approach for mapping global sector boundary structure.
Detection and Statistics of large-scale Sector Boundary Crossings in the Solar Wind: The SBCross Catalog
oai:arXiv.org:2610.28139v1
arXiv:2610.28139v1 Announce Type: new
Abstract: Heliospheric sector boundary crossings (SBCs), composed of the Heliospheric Current Sheet (HCS), Heliospheric Plasma Sheet (HPS), and surrounding plasma halo, represent fundamental structures in the solar wind that influence interplanetary dynamics and space weather. We present the \textsc{SBCross} catalog, an automated, publicly available database containing 1,076 SBC events derived from in-situ observations spanning more than 30 years (1995--2026), covering Solar Cycles 23 and 24 alongside the rising and maximum phases of Solar Cycle 25. By combining interplanetary magnetic field (IMF) and 255~eV suprathermal electron pitch angle distributions (PAD) data, our algorithm reliably captures large-scale polarity transitions, with $76\%$ of detected SBCs exhibiting clear HPS signatures ($\beta > 1$, $n_{\mathrm{peak}} > 10\mathrm{\,cm^{-3}}$). Across the catalog, we find an average occurrence rate of $\approx 2.6$ events per Carrington rotation, a median event duration of $\approx 9.0$\,hours, and a median physical thickness of $\approx 7.4 \times 10^6$\,km. Furthermore, approximately $\approx 15\%$ of detected SBCs interact with interplanetary coronal mass ejections (CMEs) on average, with this proportion rising up to $33\%$ during solar maximum. These CME-associated events exhibit systematically longer durations (median of $10.0$\,hours), $54\%$ larger physical thicknesses, and higher average solar wind velocities. Notably, all high-velocity SBC events ($> 600\mathrm{\,km\,s^{-1}}$) are exclusively driven by CME interactions, demonstrating local plasma pile-up and structural entanglement resulting from fast CMEs overtaking the ambient slow background. The \textsc{SBCross} catalog serves as an updated successor to the classic SBC list by L. Svalgaard, delivering a standardized, open-access resource for the heliophysics community to systematically investigate sector boundary transitions.
Hidden engines in tidal disruption events: multimessenger signatures in a lepto-hadronic framework
oai:arXiv.org:2610.28182v1
arXiv:2610.28182v1 Announce Type: new
Abstract: Tidal disruption events (TDEs) provide a unique laboratory for studying transient accretion, outflows, particle acceleration, and the possible connection between stellar disruptions and high-energy neutrinos. Motivated by recent candidate TDE-neutrino associations and by the forthcoming sensitivity of multimessenger facilities, we investigate under which physical conditions TDEs can be detected in MeV, GeV, and TeV gamma rays, and in high-energy neutrinos. We develop a lepto-hadronic emission framework and apply it to three representative sources: AT2019dsg, AT2019ehz, and ASASSN-14li. For each source we explore a set of physically motivated models that vary the proton and electron spectral indices, the minimum Lorentz factor of non-thermal protons, the particle energy partition, the acceleration efficiency, the magnetic compactness, and the presence of a compact corona. We compute the broadband spectral energy distributions, the associated muon-neutrino spectra, and compare them with reference sensitivities for COSI, Fermi/LAT, CTAO, and KM3NeT/ARCA. We find that hard proton spectra and/or large minimum proton Lorentz factors are the most favourable conditions for TeV gamma-ray and neutrino detectability, while GeV-bright models can be predominantly leptonic and neutrino quiet. AT2019dsg and AT2019ehz can reach the CTAO and KM3NeT benchmarks in hadronically favourable configurations, whereas ASASSN-14li is more naturally bright in the MeV-GeV band and constrained by its radio synchrotron emission. Our results highlight the importance of combining MeV, GeV, TeV, radio, and neutrino observations to distinguish between transparent, hidden, and lepton-dominated TDE emission regimes.
Orbital Parameters for 65 Wide Binary Systems: no evidence for gravitational anomaly at low acceleration
oai:arXiv.org:2610.28220v1
arXiv:2610.28220v1 Announce Type: new
Abstract: We use a Monte Carlo analysis to study whether wide binaries can be represented by bound Newtonian orbits, and, if this is the case, whether the distribution of the orbital parameters agrees with expectations. The potential orbits are classified according to the $\chi^2$ probability of representing star coordinates, parallaxes, velocity components and masses. We apply this analysis to 44 systems previously published, to which we add 26 wide binaries, consisting of twin low-mass (M) stars. We find orbital solutions for 65 systems, and derive the distribution of their orbital parameters. For one system composed of M-type stars, no Newtonian bound solution can be found; however, this pair is a member of the young Argus association and therefore not a suitable target. The remaining four systems that cannot be fit were already identified as likely disrupted binaries. The uncertainty in the separation between the components, given by the parallax error, is by far the dominant error source, and the main limitation of the present analysis, affecting especially the determination of the orbits' inclination angle and phase. Newtonian dynamics describe well our sample. When dividing the sample into two subsamples, according to the mutual acceleration, the eccentricity and periastron distributions agree with theoretical expectations; inclination and phase are distorted by the parallax-limited separation, more strongly in the high-acceleration (37 systems) than in the low-acceleration ($A_N < 10^{-9} m s^{-2}$, 28 systems) subsample, which is opposite to what a gravitational signal would produce. In particular the ratio between observed velocity difference and circular velocity does not show any excess for the low acceleration sample. When parametrized as a boost $\gamma$ parameter to the velocity, $0.95 < \gamma <1.05 (1 \sigma)$, fully consistent with Newton also for the low acceleration sample.
The Effect of Field Line Tracing Height on Wang-Sheeley-Arge Solar Wind Predictions at Earth
oai:arXiv.org:2610.28274v1
arXiv:2610.28274v1 Announce Type: new
Abstract: (Abridged abstract): A key component of the widely used Wang-Sheeley-Arge model is the tracing of magnetic field lines in the solar corona, as that largely determines the modeled open/closed field regions and in turn the predicted large-scale solar wind structure. We demonstrate that by slightly raising the height from which magnetic field lines are traced from the photosphere (R$_{su}$), we can improve the prediction of solar wind high-speed enhancements (HSEs). We find a higher R$_{su}$ produces more accurate estimates of the HSE start time, duration, and peak speed, and enables the detection of HSEs that were missed with the standard R$_{su}$. These improvements are primarily associated with changes in the great-circle angular distance, while the flux tube expansion factors remain largely unchanged. A higher R$_{su}$ produces smoother open-field boundaries and reduces small-scale fragmentation within open-field regions. These morphological changes lead to more stable magnetic footpoint mappings and consequently more robust GCD calculations. To assess the statistical relevance of these findings, we extend the analysis to 53 Carrington rotations and perform a multi-metric evaluation of model performance. The results indicate that higher R$_{su}$ improves HSE predictions while reducing false alarms. We thus find that the choice of R$_{su}$ represents an important and previously underexplored parameter that can help enhance our ambient solar wind modeling and forecasting capabilities.
Real-Time 24-hour Ahead Forecasting of Solar Energetic Particle Events with SEPNET-PRISM-R
oai:arXiv.org:2610.28292v1
arXiv:2610.28292v1 Announce Type: new
Abstract: Solar energetic particle (SEP) events, operationally defined as proton flux exceeding 10 pfu in the >10 MeV channel, threaten spacecraft systems, astronaut safety, and high latitude aviation, yet continuous real-time forecasting with adequate lead time remains difficult. We present SEPNET-PRISM-R, a near-real-time system that issues hourly, 24 hour ahead probabilistic SEP forecasts. The model integrates multimodal solar data - including active region patches, flare features, coronal mass ejection (CME) characteristics, and historical proton flux - to concurrently predict SEP occurrences, peak proton flux, soft X-ray flux, and CME velocity, with the primary objective being SEP forecasting. Building on existing literature, SEPNET-PRISM-R enhances operational performance by incorporating a postprocessing module that enables cross-validated decision making and interpretable quantification of feature importance. The training dataset spans from February 1986 to September 2025, with subsequent forecasts generated sequentially through the present via a real-time space weather prediction platform that provides predictive probabilities, feature attributions, and historical analogs, enabling more reliable and informed operational decision-making. The proposed framework substantially suppresses false alarm rates while maintaining robust event detection capabilities, as evidenced by case studies and aggregate statistics.
Coevolution of Dust Growth and Gravitational Instability in Massive Protostellar Disks
oai:arXiv.org:2610.28331v1
arXiv:2610.28331v1 Announce Type: new
Abstract: Massive protostellar disks are expected to sustain high accretion rates through gravitational instability, while dust growth may modify their opacity, cooling efficiency, and stability. Conversely, the time-dependent and nonaxisymmetric structures produced by gravitational instability may regulate the local conditions for dust growth. However, the coupled evolution of dust growth and gravitational instability in massive protostellar disks has not yet been explored self-consistently. We present two-dimensional hydrodynamic simulations that follow gas dynamics, dust growth, and opacity-dependent radiative cooling in a massive protostellar disk. We find that dust grains rapidly grow to millimeter-centimeter sizes shortly after disk formation, substantially reducing the infrared opacity. The enhanced cooling promotes disk fragmentation and produces more frequent accretion-driven luminosity outbursts than in a model without dust growth, with durations and peak luminosities comparable to those observed in massive protostars. Dust growth is broadly regulated by the fragmentation-limited size in dense disk regions, whereas newly formed outer disk regions and low-density regions deviate from this limit because of finite collision timescales and dust transport. These results suggest that dust growth may facilitate episodic accretion and fragmentation, while gravitationally unstable disk structures in turn regulate the dust growth process.
Data-Driven Ages from Asteroseismology for FGK Dwarfs
oai:arXiv.org:2610.28345v1
arXiv:2610.28345v1 Announce Type: new
Abstract: Age is the one of the most difficult stellar properties to constrain. However, upcoming surveys such as PLATO are about to revolutionize our ability to measure FGK dwarf ages with asteroseismology. Here, we use The Cannon to construct a data-driven model for inferring asteroseismic-quality age measurements from large spectroscopic survey observations and demonstrate that forthcoming asteroseismology datasets from future surveys can be leveraged, together with our model, to obtain stellar age samples larger than can be obtained by those surveys alone. We train The Cannon on 323 FGK dwarf stars with asteroseismic age labels from NASA Kepler photometry and high-resolution spectra from the Apache Point Observatory Galactic Evolution Experiment (APOGEE). We apply our model to FGK dwarfs from the Kepler survey that exhibit solar-like oscillations and reproduce their reported asteroseismic-quality ages with an RMS scatter of 1.0 Gyr. We then apply our model to the sample of 3810 Kepler FGK dwarfs with APOGEE spectra to infer data-driven asteroseismic ages. We repeat this for the larger sample of 48,828 FGK dwarfs with APOGEE spectra. Our work demonstrates that data-driven methods can robustly reproduce asteroseismic ages for FGK dwarfs. These methods will be invaluable tools for extending asteroseimology measurements from upcoming surveys (e.g., PLATO) to much larger stellar samples.
Revisiting the Unknown Unknowns in Transient Astrophysics with the SKA Telescopes
oai:arXiv.org:2610.28347v1
arXiv:2610.28347v1 Announce Type: new
Abstract: It was already recognised in the 2004 and 2014 SKA Science Case Chapters: "The Exploration of the Unknown" and "The SKA and the Unknown Unknowns", respectively, that the exploration of the unknown is an extremely important concept in realising the potential of any new instrument. This is especially true in an era of limited budgets but expanding analysis algorithms and capabilities. While this Chapter considers transients and environments, the ethos of these previous Chapters should be considered and embraced broadly and can be applied more broadly across SKAO science. We do not repeat in detail the important background these Chapters provide but encourage the reader to review them before reading this Chapter. In the last decade of radio transient astronomy, we have seen the value of these approaches with the expansion of the FRB population, the discovery of long period pulsars and transients and radio emission from white dwarfs to name just a few. With the SKAO, we should push the explored parameter space further and in new directions in all of the sensitivity, field-of-view, time, spectral, spectro-dynamic, polarisation and imaging domains. We should utilise approaches such as machine learning in allowing us to address the complexity that arises from the expanded parameter space and increased data rates and candidate sources. We will briefly summarise some of the impressive progress that the SKA pathfinders and precursors have made in these areas and what these mean for the SKA telescopes. We will also discuss ways in which it might be possible to push the current capabilities.
Photometric variability of blue supergiants: Constraints on stellar masses and evolutionary stages
oai:arXiv.org:2610.28352v1
arXiv:2610.28352v1 Announce Type: new
Abstract: Blue supergiants (BSGs) are crucial to understanding the evolution of massive stars. Although distinguishing between pre- and post-red supergiant (RSG) phases is challenging, stellar pulsations provide valuable constraints on their masses and evolutionary status. This study analyzes TESS light curves of six BSGs with high L/M to identify periodicities. We performed theoretical modeling of the radial unstable modes to characterize the strange modes. We aim to estimate asteroseismic masses, compare them with spectroscopic and evolutionary mass estimates, and assess the evolutionary stage of the stars. We employed the Lomb-Scargle periodogram and the weighted wavelet Z-transform to address quasiperiodic behavior and time-varying signals. A linear nonadiabatic stability analysis was performed to test for radial perturbations, generating eigenfrequencies that enable stellar mass estimation. In addition, we characterized the stochastic low-frequency variability (SLFV). We derive a first guess at the asteroseismic masses. We find that HD 74371 and HD 79186 are evolving blueward following the RSG phase, while the evolutionary stage of HD 115842 cannot be uniquely determined with the available mass estimates. We also find that HD 52382 and HD 99953 are likely in the post-RSG phase; however, their low masses are difficult to reconcile with the available evolutionary models. The extreme supergiant HD 80077 shows the largest mass discrepancy in the sample, suggesting a more complex evolutionary history. The SLFV properties do not provide a clear diagnostic of the evolutionary stage. These results highlight the potential of asteroseismology as a complementary tool for constraining the evolutionary status of BSGs. However, the mass discrepancies found for some objects indicate that additional observational and theoretical efforts are required.
KMOS Spectroscopic Follow-up of 16 Star Cluster Candidates of the Sagittarius Dwarf Galaxy
oai:arXiv.org:2610.28374v1
arXiv:2610.28374v1 Announce Type: new
Abstract: The Sagittarius (Sgr) dwarf spheroidal galaxy hosts several known globular clusters and numerous additional cluster candidates identified toward its crowded central regions. Establishing the nature of these candidates is difficult because apparent spatial and astrometric structure can arise from the underlying Sgr population. We investigate sixteen Minni stellar-cluster candidates with KMOS spectroscopy to determine whether the observed stars are compatible with Sgr and whether any candidate field shows evidence for a distinct stellar system. We analyse H-band KMOS spectra of 341 targets from the VVVX-GalCen Spectroscopic Survey, deriving radial velocities (RVs) and empirical metallicities. We compare these with an independent Sgr spectroscopic sample and use VVVX photometry and Gaia DR3 astrometry to test for localized stellar excesses relative to local control fields. Of the 341 KMOS targets, 256 (75%) have RVs within the adopted Sgr interval, and their metallicities broadly overlap the Sgr comparison population. No field shows statistically compelling kinematic evidence for a distinct system after accounting for the multiple candidates examined. Minni 329 shows the strongest spatial excess, but this is not significant after multiple-testing correction and is not accompanied by significant central concentration. The KMOS targets predominantly trace the Sgr population, but Sgr compatibility alone does not confirm the individual candidates as distinct stellar clusters. None of the sixteen fields can be confirmed as a separate stellar system with the present data, although this does not establish that the original candidates are spurious. Deeper imaging and larger spectroscopic samples are required to determine their physical nature.
Identifying a jet-shaped point-symmetric morphology in the supernova remnant HB9
oai:arXiv.org:2610.28428v1
arXiv:2610.28428v1 Announce Type: new
Abstract: I identify a point-symmetric morphology in the supernova remnant HB9 (G160.9+2.6), supposedly a core-collapse supernova (CCSN) remnant (CCSNR), and attribute this structure to two or more pairs of jets that exploded the star within the jittering-jets explosion mechanism (JJEM). The pairs of jets that shape the point-symmetric morphology of HB9 are part of a larger number of pairs that exploded the star. The point-symmetric structure is composed of a pair of ears, a pair of ears/nozzles, a pair of dents, and a pair of bright bars. Although these pairs of morphological features are small and the point-symmetric structure is subtle, I consider the identification robust. I use visual inspection to identify the point-symmetric morphology of HB9, a method commonly used to identify jet-shaped, point-symmetric morphologies in other astrophysical objects, such as cooling-flow clusters of galaxies and planetary nebulae. I present arguments for jet-driven shaping rather than CSM shaping. These findings support earlier claims that the JJEM is the primary explosion mechanism of CCSNe.
Tight asteroseismic upper limits on hot subdwarf magnetism
oai:arXiv.org:2610.28443v1
arXiv:2610.28443v1 Announce Type: new
Abstract: The apparent absence of magnetism in subdwarf B stars (sdBs) has emerged as an outstanding observational puzzle, especially in light of their likely origin in binary interactions. Gravity-mode pulsations in some sdBs offer a unique opportunity to extend recent successes in asteroseismic magnetic field detection in red giants. In the presence of strong magnetic fields, gravity modes are expected to be significantly damped, suppressing their observed amplitudes. We place stringent upper limits on magnetic fields in 30 hot subdwarfs by using the observation of low-frequency dipole and quadrupole modes to rule out magnetic suppression. Even under conservative assumptions, our upper limits constrain near-surface magnetic fields to below a kilogauss in roughly half of our sample, with the strongest individual constraint reaching down to $\sim20\,\mathrm{G}$. These upper limits reinforce the apparent dearth of magnetic fields in hot subdwarfs and demonstrate asteroseismology as a powerful new probe of magnetism in these stars.
The "Rubin 2036: Opportunities for Discovery" Workshop: A Community Roadmap for the Future of the Vera C. Rubin Observatory After the LSST
oai:arXiv.org:2610.28482v1
arXiv:2610.28482v1 Announce Type: new
Abstract: The NSF-DOE Vera C. Rubin Observatory will complete its ten-year Legacy Survey of Space and Time (LSST) in the mid-2030s. Its 8.4 m aperture and 9.6 square degree field of view give it an \'etendue unmatched by any facility planned for that era. The Rubin 2036: Opportunities for Discovery workshop (KIPAC/SLAC, 3-5 August 2026) brought together approximately 100 participants to discuss Rubin's future after LSST. This paper summarizes the workshop's structure and discussions across five working groups. The continued-surveying group explored a machine-learning foundation model pretrained on Rubin imaging for detection, deblending, and classification. The cadence group considered allocating post-LSST observing time and strategies for transiting exoplanets, Type Ia supernova cosmology, time-delay cosmography, and second-scale transients. The filter group explored broadband, medium-band, and narrowband filters for photometric redshifts, high-redshift large-scale structure, AGN reverberation mapping, Galactic archaeology, and emission-line imaging. The multi-messenger group explored follow-up of next-generation gravitational-wave and neutrino sources, including standard sirens, kilonovae, and Galactic supernovae. The technology group organized options by cost and lead time, including fast readout, new filters, a faster shutter, and summit edge computing. Cross-cutting topics included AI/ML, Chilean and Latin-American co-leadership, satellite megaconstellations, and weather and climate trends at Cerro Pach\'on. White papers and workshop talks are available through https://zenodo.org/communities/rubin2036workshop/ and https://rubin-observatory.github.io/Rubin-Beyond-LSST/pages/Rubin2036proceedings/.
A Solar Optical Clock
oai:arXiv.org:2610.28485v1
arXiv:2610.28485v1 Announce Type: new
Abstract: Precision spectroscopy underpins quantum sensing, atomic timekeeping, and the search for terrestrial-mass exoplanets. Finding an Earth-Sun analog requires fractional splitting of spectral features that rivals the precision of the best laboratory spectroscopy. This corresponds to resolution of a few parts in $10^5$ of a typical stellar spectral line, or a Doppler velocity precision of less than 10 cm s$^{-1}$. Only a few astrophysical spectrometers reach this regime, and only by dispersing starlight with large gratings and averaging thousands of lines. Here, we show that comparable precision can be achieved on a single solar line by leveraging techniques used in quantum sensors and optical clocks. A laser is frequency-locked to the 1565 nm solar Fe I line by coupling to a single mode of the solar thermal field and counted against an optical frequency comb referenced to the SI second, thus creating a solar optical clock. Operating at the fundamental photon-noise limit of thermal heterodyne spectroscopy, we achieve centroid precision of one part in $10^9$ of the optical frequency, or 30 cm s$^{-1}$ in radial velocity. We recover the Sun's five-minute acoustic oscillations and the expected Rossiter-McLaughlin signal during a solar eclipse. Combining these techniques from quantum sensing with a thirty-meter-class telescope could enable our demonstrated precision on the nearest Sun-like stars in a few hours. Multiplexing such absolute frequency measurements across lines with different magnetic sensitivity and formation depth would reduce averaging time and help disentangle stellar granulation and activity from the Doppler signature of terrestrial exoplanets, providing new approaches to address some of the most demanding measurements in astrophysical spectroscopy.
Resonant three-mode coupling in asteroseismology I: Bifurcations, chaos, and DAV white-dwarf outbursts
oai:arXiv.org:2610.28513v1
arXiv:2610.28513v1 Announce Type: new
Abstract: Resonant three-mode coupling is a fundamental nonlinear process in stellar oscillations and other wave systems. We investigate the global dynamics of a linearly overstable parent mode coupled to two damped daughter modes. The principal control parameters are frequency detuning $\Delta$ and daughter damping rate $\Gamma$, normalized by the parent growth rate. We combine high-resolution maps of the largest Lyapunov exponent in the $(\Gamma,\Delta)$ plane with linear stability analysis, Floquet calculations, and numerical continuation. The resulting parameter space includes stable equilibria, simple and period-doubled limit cycles, chaotic attractors, formally unbounded solutions, and higher-period `shrimp' islands within chaotic regions. Successive period doublings lead to chaos, while boundaries of isolated periodic windows exhibit type-I intermittency.
At large $\Gamma$, the main bifurcation boundaries approach straight lines, indicating that the dynamics is governed primarily by $r=\Delta/(2\Gamma)$. We connect finite-$\Gamma$ bifurcations to an earlier one-dimensional return map. As $\epsilon\rightarrow0$, this map recovers the principal period-doubling sequence, the onset of chaos near $r\simeq0.38$, and the asymptotic positions of prominent periodic windows. We also examine the effects of unequal daughter damping rates.
Applying these results to sporadic outbursts in cool DAV white dwarfs, we find that estimated $g$-mode growth rates, damping rates, and detunings place resonant triplets in both periodic and chaotic regimes. Thus, instability of the three-mode equilibrium need not produce a regular limit cycle: deterministic chaos and intermittency may contribute to irregular outburst amplitudes and recurrence times.
The clumpy fraction of galaxies at $z=1-9$ with JADES: evidence for a mass-dependent transition in clump origin
oai:arXiv.org:2610.28515v1
arXiv:2610.28515v1 Announce Type: new
Abstract: We measure the fraction of clumpy galaxies ($f_{\rm clumpy}$) over the redshift range $z=1-9$ and stellar masses $\log_{10}(M_\star/{\rm M_\odot}) \geq 8$ using JWST/NIRCam imaging from JADES. We define clumps as compact off-centre surface-brightness peaks with characteristic sizes of a few hundred parsecs, and develop a dendrogram-based clump-finding method with astrodendro to identify them directly in rest-frame UV, blue, and optical images. Clump-detection completeness is quantified through injection-recovery simulations, which we use to define ($M_\star$, $z$)-dependent 50 per cent completeness-limited $f_{\rm clumpy}$, alongside a probability-weighted completeness-corrected statistic for comparison. For the full sample, the blue and optical $f_{\rm clumpy}$ reach $\approx55$ per cent at $z\approx1-2$ and decline to $\lesssim5$ per cent by $z\approx8-9$, while the UV fraction remains higher at $\approx12$ per cent. At $z\approx1-2$, $f_{\rm clumpy}$ increases strongly with stellar mass, from $\approx35$ per cent in the lowest-mass bin to $\approx85$ per cent at $\log_{10}(M_\star/{\rm M_\odot})\geq10$. Most clumpy galaxies host only one or two robustly detected off-centre clumps, while higher-multiplicity systems occur preferentially in massive, strongly star-forming galaxies at low to intermediate redshift. Comparison with major merger rates suggests that mergers could account for a substantial fraction of low-mass ($\log_{10}(M_\star/{\rm M_\odot})<8.5$), high-redshift clumpy galaxies, while the observed fractions in more massive systems ($\log_{10}(M_\star/{\rm M_\odot})>9$) require an additional in situ clump-formation channel. Our results are consistent with a transition from merger-dominated clumpy morphologies at early times and low stellar masses to an increasing contribution from in situ clump-formation channels, such as disc instabilities, in more massive systems.
Two Validated Planets and 179 Candidates From a TESS Faint Star Transit Search in the ULTRASAT Fields
oai:arXiv.org:2610.28518v1
arXiv:2610.28518v1 Announce Type: new
Abstract: The TESS mission is designed to survey the entire sky to find planets most amenable to follow-up observations, and transit searches have therefore focused specifically on bright targets. Faint stars ($T_{mag}>13.5$) have not been systematically searched for planets, despite TESS having the sensitivity necessary to do so, therefore representing a promising avenue for increasing TESS's scientific yield. We utilized TESS's sensitivity for faint stars in particular to find more planets for joint optical-UV studies of exoplanet transits. Observing transits in both UV and optical gives insight into exoplanet atmosphere composition and evolution and remains relatively unexplored due to the lack of UV-capable instruments. ULTRASAT is an upcoming UV space telescope that will perform a high cadence survey suitable for observing exoplanet transits. To increase the sample size of planets capable of observations with ULTRASAT, we performed a transit search of faint stars down to $T_{mag}=16$ in the ULTRASAT high cadence survey fields using TESS light curves optimized for faint stars. Our search recovered 262 planet candidates, of which 181 are new to this work. We also statistically validated 2 of those candidates as bona fide planets. Finally, we estimated the ULTRASAT SNR of both new and existing planet candidates, and found that 113 targets will be detectable with ULTRASAT, which would facilitate the largest study of transiting planets in the visible and UV. This work demonstrates the power of TESS to detect exoplanets around faint stars and provides a catalog of well-vetted transiting planet candidates to observe with ULTRASAT.
TOI-3686: A Hot Jupiter and Sub-Neptune In A 2:1 Mean-Motion Resonance
oai:arXiv.org:2610.28533v1
arXiv:2610.28533v1 Announce Type: new
Abstract: Hot Jupiters are typically thought to be "lonely", but a small and growing subset of them have been found to have nearby planets in the same system. We present the discovery of a new member of this class of systems, TOI-3686. TOI-3686 c is a hot Jupiter with $P_\mathrm{orb} = 5.078$ days detected by NASA's TESS mission that shows peak-to-peak transit timing variations (TTVs) of 40 minutes, amongst the largest ever observed for a hot Jupiter. A subsequent search of the TESS data revealed the transits of a second planet, TOI-3686 b, a sub-Neptune with $P_\mathrm{orb} = 2.54$~days and peak-to-peak TTVs of 10 hours, $\sim 16\%$ of its orbital period. We characterized the system through photodynamical modeling of photometry from TESS and ground-based facilities, as well as radial velocity measurements from the NEID and MAROON-X spectrographs. We find that the two planets in the TOI-3686 system are deep in a 2:1 mean-motion resonance, making this the only hot Jupiter system currently known to be in a truly resonant configuration, which may be related to its $\lesssim1$~Gyr system age. Systems like TOI-3686 likely formed through a dynamically quiescent pathway and were captured into resonance after migration through a disk, and may represent the high-mass end of processes that occur in compact systems of smaller planets.
A sub-Earth-sized TESS planet candidate in the nearby M-dwarf binary GJ 774, most likely orbiting GJ 774 A (TIC 351339274)
oai:arXiv.org:2610.28575v1
arXiv:2610.28575v1 Announce Type: new
Abstract: Small planets transiting the nearest M dwarfs are the terrestrial planets most accessible to detailed characterisation. We report a sub-Earth-sized transiting planet candidate in the wide (17.6") M-dwarf binary GJ 774 at 12.86 pc, most likely orbiting the M2.5 primary GJ 774 A (TIC 351339274). The signal was a SPOC multi-sector threshold-crossing event (MES = 7.95) and is not a TESS Object of Interest. In eight TESS sectors (2019-2026) a 1.28-h transit recurs every 4.225138 +/- 0.000010 d with a box-mean depth of 365 +/- 32 ppm in A's light. The 2026 data, which entered neither the SPOC search nor our selection, show it at the ephemeris predicted from 2019-2023 data (304 +/- 43 ppm; none of 10^5 phase-scrambled versions of the 2026 data is as deep). The per-target false-alarm probability is < 0.03 at 95%. Odd/even, secondary-eclipse, archival-plate and Gaia radial-velocity tests exclude unrelated blended eclipsing binaries and eclipsing companions on 4.2-d orbits, and hierarchical eclipsing binaries are improbable. The stars are 0.85 TESS pixel apart, so the TRICERATOPS false-positive probability (0.26) is set by a planet on GJ 774 B. B is disfavoured because its own light does not dim enough: PRF-deblended photometry gives 109 +/- 115 ppm in B's light, against 1089 ppm required if B were the host. Conditional on a point source at A or at B, the joint host model gives P(A)=0.59-0.97 (0.97 for the native target-pixel baseline). Over all sectors the PRF photometry is consistent with A, but in the out-of-sample 2026 data 40% of the deficit appears at B's position, an unexplained 2-3 sigma (PSF-model and depth dependent) tension with host A. On A the planet has R_p=0.70 +/- 0.06 R_Earth, S=15.0 +/- 1.8 S_Earth and T_eq=548 +/- 16 K. It is not validated; resolved ground-based photometry of the pair would test host B directly.
Large fluctuations and cusp-core transformation in dwarf galaxies
oai:arXiv.org:2610.28583v1
arXiv:2610.28583v1 Announce Type: new
Abstract: Dwarf galaxies are abundant and dark matter dominated, rendering them particularly relevant for testing a cosmological model dominated by a dark sector. Whether their density profiles agree with $\Lambda$CDM expectations has long been a point of contention, dubbed the cusp/core problem. While numerical simulations have suggested that the erasure of central dark matter cusps may arise from potential fluctuations in feedback-driven gas, whether the required fluctuation levels are achieved in actual galaxies has remained uncertain. By measuring the power spectra of neutral gas density contrasts in a sample of gas-rich dwarf galaxies, we show that the observed fluctuations match theoretical expectations when their impact is derived from a stochastic dynamical model. Specifically, the relaxation times for core formation at the observationally estimated radii are generally on the order of a Hubble time or less. Core sizes are therefore consistent with conservative estimates assuming their generation through steady fluctuations acting over cosmic history. The relaxation time increases with halo mass, suggesting that core formation via this mechanism becomes less efficient in more massive systems.
Resolving the extended interstellar medium of infrared galaxies II. Fine-structure line emission
oai:arXiv.org:2610.28597v1
arXiv:2610.28597v1 Announce Type: new
Abstract: We investigate how star formation and active galactic nuclei (AGNs) shape the ionized and neutral interstellar medium in 67 nearby luminous and ultraluminous infrared galaxies from the Great Observatories All-sky LIRG Survey (GOALS). Combining mid-infrared (MIR) spectral mapping from the Spitzer Infrared Spectrograph with far-infrared (FIR) spectroscopy from the Herschel Photodetector Array Camera and Spectrometer, we measure fine-structure emission lines in aperture grids extending to projected galactocentric distances of about 14 kpc. We compare the measurements with CLOUDY photoionization models to constrain the ionization parameter, hydrogen density, stellar age, and metallicity. In starburst-dominated galaxies, the combined MIR+FIR fits yield a decreasing ionization parameter with radius, approximately constant hydrogen densities, uniformly young stellar populations, and flat metallicity profiles. MIR-only fits instead yield an increasing ionization parameter with radius, highlighting the importance of FIR diagnostics, particularly the [C II] emission associated with photodissociation regions, in constraining the multiphase gas. In the strongest AGN-dominated system, the nucleus and inner extended regions show elevated [Ne III]/[Ne II] and suppressed [C II]/[O I] ratios, while the outermost regions, at projected distances of about 6 kpc and beyond, show values typical of star-forming regions. This suggests that the AGN may affect gas excitation out to about 5 kpc. Our results demonstrate the importance of combining MIR and FIR diagnostics to recover radial trends that can be hidden in galaxy-integrated measurements. Together with Paper I, they provide a benchmark for spatially resolved studies of gas and dust in starbursts and AGN-hosting galaxies.
Resolving the extended interstellar medium of infrared galaxies I. Polycyclic aromatic hydrocarbon emission
oai:arXiv.org:2610.28599v1
arXiv:2610.28599v1 Announce Type: new
Abstract: We investigate spatial variations in polycyclic aromatic hydrocarbon (PAH) properties and photoelectric heating efficiency in 67 nearby luminous infrared galaxies from the Great Observatories All-sky LIRG Survey (GOALS). Using mid-infrared spectroscopy from Spitzer and far-infrared cooling-line maps from Herschel, we extract spectra in concentric annuli and matched aperture grids out to about 14 kpc. PAH features are measured with the CAFE spectral decomposition tool. We trace photoelectric heating efficiency by combining PAH emission with cooling lines after isolating the neutral-gas contribution. In starburst galaxies, the PAH 11.3/7.7 ratio rises out to about 5 kpc while 6.2/7.7 remains nearly constant, indicating predominantly large, neutral PAHs. Between 5 and 8 kpc, 11.3/7.7 continues to rise while inferred PAH sizes decrease. Beyond 8 kpc, 11.3/7.7 falls and 6.2/7.7 rises, indicating smaller, more ionized PAHs. These trends are consistent with PAH photodestruction near the central starburst. Strong active galactic nucleus (AGN) systems show more neutral PAHs and harder radiation fields in their nuclear and inner regions. Beyond about 5 kpc, PAHs become more ionized and approach the starburst locus as the radiation field weakens. Strong AGNs influence average PAH charge states at 6-14 kpc as effectively as moderate AGNs at 1-6 kpc. Photoelectric heating efficiency increases with radius in both starburst and AGN-hosting galaxies as the radiation field weakens. In starbursts, higher efficiency is associated with more neutral PAHs, while AGN-dominated systems show the opposite relation. These results provide quantitative, spatially resolved constraints on the extent of AGN effects on dust properties in a population of luminous infrared galaxies.
Automated optimization of solar energetic particle acceleration and transport parameters against five-spacecraft observations of the 9 October 2021 event
oai:arXiv.org:2610.28603v1
arXiv:2610.28603v1 Announce Type: new
Abstract: Solar energetic particle (SEP) events threaten space infrastructure, yet physics-based prediction models depend on numerous poorly constrained free parameters. We introduce the Fujitsu Solar Energetic particle Intelligent Tuning and Evaluation Network (SEITEN), which couples the SUSANOO MHD model of CME propagation to SEPIF, a diffusive shock acceleration and focused transport model, through the shock conditions SUSANOO supplies, and automatically optimizes the free parameters of the particle model against the SEP event of 9 October 2021. The mean free path and the injection efficiency vary independently along each magnetic field line and heliolongitudinal bin, giving a 19-dimensional problem solved by quasi-Monte Carlo sampling followed by a Tree-structured Parzen Estimator. The optimized model reproduces proton flux profiles at five spacecraft with a mean absolute error of $\approx 0.45$ in $\log_{10}$ flux, and the inferred mean free paths show a clear east-west structure; these are effective parameters conditional on the assumed transport geometry. Six independent runs show that whether mean free paths or injection efficiencies dominate depends on the MHD shock strength, and that injection efficiencies differ by orders of magnitude between MHD configurations despite nearly identical fit quality. Joint MHD-SEP parameter estimation is therefore a necessary next step towards generalized, operational SEP forecasting.
Community Lessons from builders of infrared multi-object spectrographs
oai:arXiv.org:2610.28605v1
arXiv:2610.28605v1 Announce Type: new
Abstract: We present a summary of lessons learned from our community on the design and construction of infrared multi-object spectrographs. As members of the Magellan Infrared Multi-Object Spectrograph (MIRMOS) development team, we were motivated to learn about the most common troubles faced by our colleagues and incorporate measures in our design and assembly plan to address known issues. The lessons are compiled as a result of interviews we carried out with members of the instrument development teams for several multi-object infrared spectrographs. We found considerable overlap in the lessons learned for building this type of instrument; many of which relate to the cryogenic operation of the instruments. In this paper, we will present the compilation of general and specific lessons which we found relevant for the successful development of MIRMOS or any other future cryogenic instruments.
Slow-Roll Inflation at Arbitrary Order
oai:arXiv.org:2610.28608v1
arXiv:2610.28608v1 Announce Type: new
Abstract: We derive explicit expressions for the primordial scalar and tensor spectra at arbitrary finite order $p$ in the slow-roll expansion, valid both for canonical single-field inflation and $k$-inflation. At any order, the spectra are determined by known combinations of slow-roll parameters weighted by universal coefficients entirely fixed by the free de Sitter solutions, without requiring perturbatively corrected mode functions. Using these spectra, we obtain closed-form consistency relations that connect the tensor tilt $n_t$ and its runnings $\alpha_n^t\equiv d^n n_t/d(\ln k)^n$ to the tensor-to-scalar ratio $r$, the scalar spectral index $n_s$, and its runnings $\alpha_n^s\equiv d^n(n_s-1)/d(\ln k)^n$. The consistency relations are valid for arbitrary $n\geq 0$ and any finite $p\geq 0$, without recursion over lower-order relations. For $k$-inflation, they additionally depend on the scalar sound speed and its scale derivatives. We develop an algorithm that generates both the spectra and consistency relations for any pair $(n,p)$. We reproduce known results through next-to-next-to-next-to-leading order (N3LO) and extend the explicit calculations up to N10LO.
Dynamical stability of Magnetic Supramassive Neutron Stars
oai:arXiv.org:2610.28609v1
arXiv:2610.28609v1 Announce Type: new
Abstract: We investigate the dynamical stability of strongly magnetized supramassive neutron stars close to the mass-shedding limit using fully general-relativistic magnetohydrodynamic simulations in three spatial dimensions. Equilibrium configurations are constructed using a microphysical equation of state and considering poloidal magnetic fields up to maximum strengths ${\sim} 10^{17}\,\mathrm{G}$. The initially unstable magnetic fields reorganize into mixed poloidal--toroidal configurations on Alfv\'en timescales shorter than the stellar rotation period. The subsequent dynamical evolution is governed by timescales set by the Lehnert number. For stars with masses sufficiently lower than the maximum mass, magnetic stresses redistribute angular momentum and drive the development of differential rotation, leading to a new quasi-stationary configuration that remains dynamically stable over several rotation periods. Magnetically-driven winds develop from the outer layers around the equatorial plane and remain baryon-matter dominated with magnetization ${\lesssim}0.2$. Over the simulated timescale, winds carry a total rest mass of order ${\sim}10^{-3}\mathrm{M}_\odot$ at outflow rates of up to a few $M_\odot/\mathrm{s}$. For stars close to (or beyond) the maximum-mass turning point along the stellar equilibrium sequence, strong magnetic fields trigger gravitational collapse by redistributing angular momentum away from the stellar core. The resulting black hole remnants are characterized by angular velocities up to ${\sim}49\%$ larger and angular momentum ${\lesssim}8\%$ smaller than those of the initial equilibrium stars. Our results show that extreme magnetic-field instabilities can significantly alter the structure of supramassive neutron stars, with possible implications for magnetars and long-lived binary neutron-star merger remnants.
Constraints on Millimeter-wave Extragalactic Transient Populations from 4 Years of SPT-3G data
oai:arXiv.org:2610.28633v1
arXiv:2610.28633v1 Announce Type: new
Abstract: Wide-field cosmic microwave background (CMB) instruments have reached sufficient sensitivity and cadence to conduct untargeted millimeter-wave searches for astronomical transients. To date, such searches have primarily discovered stellar flares, while extragalactic detections remain limited to a few flaring blazars and active galactic nuclei. In contrast, recent simulations suggest that current and future CMB experiments should detect a few to dozens of extragalactic transients per year from sources including gamma-ray bursts (GRBs), tidal disruption events, and fast blue optical transients. At timescales under one week, predicted detections are dominated by flares from reverse shocks in long-duration GRBs, with some benchmark models predicting rates of ~10 detections per year for the SPT-3G camera on the South Pole Telescope. In this work, we conduct a transient search optimized for the bright and brief millimeter-wave emission from GRB reverse shocks. We update previous simulations by incorporating a more realistic model of the SPT-3G sensitivity, observing cadence, and detection algorithm, resulting in nearly a factor of two improvement in detection sensitivity compared to prior estimates. We subsequently carry out a search over four seasons of the 1500 deg2 SPT-3G Main field observations from 2019-22, finding no candidate extragalactic events. We use the lack of detections to constrain the parameters of the synthesized GRB populations, setting upper limits on the fraction of GRBs with reverse shock components, f_RS, and the all-sky rate of reverse shocks beamed at Earth, N_all-sky. For a standard 10^51 erg GRB with Lorentz factor Gamma = 100, we constrain f_RS <0.57 and N_all-sky <1226 bursts/year (at 90% confidence). The most readily detectable model we consider is a 10^52 erg "dirty" GRB with Gamma = 50, whose population we constrain to f_RS <0.09 and N_all-sky <120 bursts/year.
Search for Compact Binary Coalescences with Sub-solar Mass Components in Data from the Second Part of LIGO-Virgo-KAGRA's Fourth Observing Run
oai:arXiv.org:2610.20276v1
arXiv:2610.20276v1 Announce Type: cross
Abstract: We report on a gravitational wave search for compact binary coalescences with secondary mass between 0.2 M and 1 M and primary mass between 0.2 M and 10 M, restricting the component mass ratio to lie between 0.1 and 1. The analysis uses data collected by the LIGO and Virgo detectors between 2024 April 10 15:00:00 and 2025 January 28 17:00:00, the second part of the fourth observing run of the LIGO-Virgo-KAGRA Collaboration. Results from two out of three search pipelines show outlier candidates, but not at a sufficient level of significance to provide evidence for detection. We report the combined sensitivity of the current searches over the first and second parts of the observing run to the target sub-solar mass (SSM) black hole (BH) population, as well as to neutron star (NS) binaries with one or more sub-solar component. We place upper limits on the merger rate of SSM BHs, ranging from $50 {\rm Gpc}^{-3} {\rm yr}^{-1}$ to~$5000 {\rm Gpc}^{-3} {\rm yr}^{-1}$ (90% confidence) for different component mass combinations. We use the merger rate limits to constrain two illustrative dark matter (DM) scenarios that can form SSM compact objects: primordial BHs, and dark BHs forming in a dissipative DM model. Our search limits the fraction of DM in primordial BHs, assuming a monochromatic mass distribution, to be $\leq$ (0.018-0.032) at 1 M, and $\leq$ (0.10-0.18) at 0.35 M for early-forming binaries; or to be less than 1 at masses above (0.62-0.74) M for late-forming binaries. For the dissipative model, we constrain the fraction of DM for a range of minimum dark black hole masses below $1\,M$ to be as low as (1.0-1.1) x $10^{-5}$; we also exclude minimum masses above (0.62-0.74) M, if dark BHs constitute all of the DM. For binary NSs that include sub-solar mass components, assuming a simple, fixed population model we find a merger rate upper limit of $249 {\rm Gpc}^{-3} {\rm yr}^{-1}$.
Adaptive Model-Switch Thermodynamic Integration: A Bayes Factor Estimator
oai:arXiv.org:2610.24354v1
arXiv:2610.24354v1 Announce Type: cross
Abstract: Bayes factors are a widely used model comparison criterion, but are often quite challenging and computationally expensive to estimate well. We propose a variation of the model-switch thermodynamic integration method for estimating Bayes factors that decreases the variance of the estimator and essentially eliminates integration error. We choose samples whose density in $\beta$ is proportional to the standard deviation of previously drawn samples. Rather than averaging many samples at each $\beta$, we use a very dense set of $\beta$ values, collecting one sample at each in order of increasing $\beta$ and one sample at each in order of decreasing $\beta$ to minimize "thermic lag". We study our procedure with a toy model (based loosely on some signal recovery problems). In this model we find that adaptive sampling gives us a 2.7-fold decrease in the spread of the estimated Bayes factor. We test our procedure both drawing directly from the posterior distribution and using a Markov Chain Monte Carlo sampler.
Ion Heating in Pair-loaded Relativistic Turbulence
oai:arXiv.org:2610.26731v1
arXiv:2610.26731v1 Announce Type: cross
Abstract: Pair-loaded relativistic turbulence is expected in compact high-energy astrophysical environments and may sustain plasma heating. Using self-consistent three-dimensional radiative particle-in-cell simulations, we identify the mechanism of ion heating mediated by pair-enhanced radiative feedback in relativistic turbulence. We derive the theoretical criterion for the transition to a regime dominated by radiative feedback and show that pair loading promotes the growth and persistence of large, connected optically thick structures. These long-lived structures sustain intense localized ion heating, repeatedly generating prominent hot-ion clumps and substantially enhancing the mean ion temperature. Our results establish a physical link between pair creation and sustained ion heating, providing a theoretical basis for investigating intermittent hot-ion production in black-hole coronae and gamma-ray burst outflows.
Floquet Instability of Photon-Axion Conversion
oai:arXiv.org:2610.26908v1
arXiv:2610.26908v1 Announce Type: cross
Abstract: The discovery of axions would provide evidence for physics beyond the Standard Model. In this work, we investigate axion-photon conversion in the presence of a background magnetic field with simultaneous spatial and temporal periodicity. By performing a Lorentz transformation, the coupled axion-photon system is mapped onto a time-periodic linear system that can be analyzed using Floquet theorem. We find unstable modes associated with resonant crossing between Floquet replicas. The resulting complex Floquet exponents lead to exponentially growing modes, which we identify as a Floquet instability. Then, we determine the parameter regions in which this instability appears in a perturbative analysis. We further show that the Floquet instability enhances the axion-field amplitude generated from an initial photon field. Our results demonstrate that spatiotemporally periodic magnetic backgrounds provide a distinct mechanism for modifying axion-photon conversion through Floquet instability.
Elementary-Symmetric Ricci Gravity: Metric Formulation
oai:arXiv.org:2610.26980v1
arXiv:2610.26980v1 Announce Type: cross
Abstract: We investigate a four-dimensional metric theory of gravity whose Lagrangian is a linear combination of the elementary symmetric polynomials of the Ricci tensor. We derive the field equations and characterize Einstein vacua and a real, diagonalizable class of parallel-Ricci vacua. The latter includes direct products of two constant-curvature surfaces and can exhibit a hierarchy of curvature radii. The Einstein vacua possess multiple branches associated with the distinct roots of a polynomial, whose multiple roots define boundaries in parameter space. We derive the linearized equations on maximally symmetric backgrounds, distinguish the scalar and spin-2 degrees of freedom, and then analyze their classical linear stability. We also obtain the generalized Friedmann constraint and the associated homogeneous evolution equations. An illustrative spatially flat solution evolves from the neighborhood of an unstable de Sitter solution toward a stable de Sitter solution. This behavior may motivate further investigation into early- and late-time accelerated expansion and provide a gravity-driven mechanism for automatic end of inflation.
Annihilation feedback on elastic scattering in self-interacting dark matter
oai:arXiv.org:2610.27056v1
arXiv:2610.27056v1 Announce Type: cross
Abstract: Sizable dark-matter self-interactions, motivated by observations of galactic structure, often arise from the exchange of light force mediators. Light mediators can render scattering nonperturbative and make inelastic processes important. Unitarity, via the optical theorem, implies that elastic and inelastic interactions cannot be treated independently: inelasticity induces absorptive contributions to the two-particle self-energy kernel that determines the elastic amplitude, while the resulting incoming-state dynamics also modifies the inelastic amplitudes. We study the feedback of dark-matter annihilation on elastic scattering in the presence of a long-range force. Although annihilation is a short-distance process, it can be strongly enhanced by such a force. We find that its feedback becomes significant near Sommerfeld resonances, where it suppresses the low-velocity growth of the baseline elastic cross section. Near baseline elastic-scattering zeros, it may instead generate a nonzero elastic floor required by unitarity. A companion paper studies the feedback of radiative bound-state formation on elastic scattering, which can generate sizable self-interactions.
A Natural Dark Matter-Dark Energy Interaction from Axion-Monopole Dynamics
oai:arXiv.org:2610.27059v1
arXiv:2610.27059v1 Announce Type: cross
Abstract: We propose a natural framework for interacting dark matter and dark energy based on axion-monopole dynamics in a dark gauge sector. Dark matter is composed of dark 't Hooft-Polyakov monopoles, while dark energy is an ultralight axion whose coupling to the dark gauge field induces an axion-dependent dyon mass through the Witten effect. This generates a direct dark matter-dark energy interaction. We show that generic models with a significant dark matter-dark energy coupling suffer from a naturalness problem: the Coleman-Weinberg potential induced by the interaction is typically much larger than the dark energy scale. In our model, the corrections arise from a monodromy tower of dyon states, allowing them to be exponentially suppressed. In the minimal setup, the interaction strength is predicted to be $\sim 0.1\%$, below the sensitivity of current observations. We discuss scenarios in which cosmological population of higher-charge dyons can enhance the interaction to the percent level, producing an apparent phantom dark energy equation of state. We further show that the parameter space of the axion-monopole model preferred by fitting the Dark Energy Spectroscopic Instrument (DESI) results can also naturally explain the cosmic birefringence with a small anomaly coefficient.
Deglobalising the fit: how to search for physics beyond vacuum General Relativity in LISA catalogues
oai:arXiv.org:2610.27069v1
arXiv:2610.27069v1 Announce Type: cross
Abstract: The Laser Interferometer Space Antenna will observe up to $\sim 10^{7}$ overlapping signals, which have to be characterised simultaneously with a global fit. This is a challenge for model comparison since rerunning the global fit for every physical model of interest, including deviations from vacuum General Relativity, is computationally demanding. An alternative is a two-step post-processing strategy: a first global fit with vacuum models only, after which a selected group of target sources is reanalysed under a different model with all other sources subtracted. We show that in this strategy a beyond-vacuum effect in a target source can be partly absorbed by the overlapping nuisance sources in the first fit, subtracted with them, and never reach the reanalysis. In the linear-signal approximation, the loss is captured by a single number, the masking ratio $\mathcal{R}\in[0,1]$. A deviation the global fit would detect at $n_\sigma$ is reported by post-processing at $\mathcal{R}\,n_\sigma$, and bounds in null tests are also reported too tight by the same factor. We show how to compute $\mathcal{R}$ from the posterior samples of the vacuum global fit, and give an algorithm that ranks the nuisance sources by how much of the effect each can absorb, so that only those most correlated with it need to be refitted with the target in a down-scaled reanalysis with a few sources. We validate our method with full Markov-chain Monte Carlo inference on a toy model of 21 overlapping chirping sinusoids, where one carries a small non-linear chirp that represents the missed physics, while the other 20 act as nuisance sources. This method can serve as a basis for tests of General Relativity in realistic LISA mock datasets and can be adapted to study the impact of systematics between different gravitational-wave sources.
Cosmic birefringence from a stiff axion field and its cosmological implications
oai:arXiv.org:2610.27082v1
arXiv:2610.27082v1 Announce Type: cross
Abstract: We propose a model in which the observed cosmic microwave background birefringence arises from the interaction of photons with a stiff ($w=1$) pseudo-scalar field that today contributes a subdominant fraction of the total energy density, distinct from both dark matter and dark energy. Big Bang Nucleosynthesis constrains the cosmological density of the stiff field, leading to a lower bound for its coupling with photons, $g_{\phi\gamma} > 8.9 \times 10^{-14}\,\mathrm{GeV}^{-1}$. Photon conversion in cluster-hosted quasars sets an upper bound, $g_{\phi\gamma} < 5.0 \times 10^{-13}\,\mathrm{GeV}^{-1}$. Therefore, there is a narrow window for the discovery of this cosmological component that can be explored in the forthcoming years with astrophysical axion searches. Moreover, we show that this scenario leads to gravitational-wave signals that can be tested with planned experiments.
Modular Bayesian Inference for Partially Misspecified Mixtures with Uncertain Class Labels
oai:arXiv.org:2610.27085v1
arXiv:2610.27085v1 Announce Type: cross
Abstract: Scientific inference often involves observations from multiple latent populations whose class memberships are uncertain and must be inferred from available data. While full Bayesian mixture modelling provides a principled approach to joint uncertainty quantification in this setting, often the model is misspecified or less trustworthy for populations that are not the target of the inference, which can distort posterior inference for the parameters of interest. We propose a modular Bayesian framework that safeguards against such biases. We decouple classification from downstream inference using fixed class probabilities, for example from a classifier trained on a subset of labelled data. We then protect inference on the parameters of interest by preventing the likelihood contributions of less trustworthy populations from influencing them. A key benefit of this principled modularization is that it greatly simplifies inference, reducing a fully joint mixture analysis to sampling label allocations, fitting the protected downstream model for each sample, and averaging the resulting posteriors. We establish calibration conditions on these class probabilities under which the proposed procedure preserves the inferential target of the analysis with known labels, and demonstrate the method in synthetic experiments and in cosmological analysis of photometric supernova samples.
Self-interacting dark matter from radiative bound-state-formation feedback
oai:arXiv.org:2610.27089v1
arXiv:2610.27089v1 Announce Type: cross
Abstract: Observations of galactic structure motivate sizable dark-matter self-interactions at low velocities, naturally realized through light force mediators. Such mediators can also make inelastic processes very efficient. Unitarity then requires their absorptive contributions to the two-particle self-energy kernel to be resummed, which in turn feeds back into both the elastic and inelastic scattering amplitudes. We study the feedback of radiative bound-state formation on elastic scattering and its implications for self-interacting dark matter. We focus on capture processes through emission of a charged scalar; such transitions are particularly efficient because the scattering and bound states are governed by different effective potentials, a feature also shared by radiative capture through emission of a non-Abelian gauge boson. We show that their feedback can generate sizable elastic scattering relevant for halo dynamics, even when the baseline elastic amplitude vanishes. Notably, the generated elastic cross section can grow toward low velocities, resembling long-range scattering. A companion paper examines the feedback from annihilation processes.
Illuminating Light Dark Matter with Fixed-Target Experiments
oai:arXiv.org:2610.27131v1
arXiv:2610.27131v1 Announce Type: cross
Abstract: Astronomical and cosmological observations provide overwhelming evidence for a non-luminous matter component, dark matter (DM), whose particle nature remains unknown. Experimental searches for non-gravitational DM interactions have largely targeted DM with GeV-TeV masses via nuclear recoils, yielding no discovery thus far. DM in the MeV-GeV mass range, however, is theoretically motivated, reproducing the observed relic abundance while remaining invisible to nuclear recoil searches. Such light DM is better probed via electron recoils or at fixed-target experiments. This thesis advances the theoretical and statistical toolkit for light DM searches at fixed-target experiments, spanning model building, statistical inference, global fits, and signal simulation. We survey the landscape of fixed-target searches for light DM, focusing on the Light Dark Matter eXperiment (LDMX). We also simulate DM signatures at the next-generation experiment SHiP to assess its sensitivity. Within minimal DM models, we present a likelihood-based statistical framework for LDMX using the full two-dimensional recoil electron distribution, spanning exclusion, discovery, parameter estimation, and model comparison in frequentist and Bayesian formulations. This establishes the discovery potential of LDMX, and we quantify the direct detection exposure needed to independently confirm a potential signal. Beyond this paradigm, we perform global fits of scalar and fermionic light DM, combining cosmological, astrophysical, and laboratory constraints. We then extend fixed-target searches to spin-1 DM and higher electromagnetic moment dark photons, yielding rich signatures with strong complementarity to other probes. Finally, we go beyond standard nuclear modelling of dark photon production, incorporating elastic nuclear form factors and spectral functions from many-body ab initio methods in place of analytic treatments.
A Multipole Technique for the Simulation and Analysis of Gravitational Waveforms from Extreme Mass Ratio Inspirals
oai:arXiv.org:2610.27153v1
arXiv:2610.27153v1 Announce Type: cross
Abstract: We introduce EMRICode, a Python-based numerical framework for evolving generic extreme mass-ratio inspirals (EMRIs) in Kerr spacetime. The code combines adaptive numerical integration of test-particle geodesics with a prescribed model for the secular evolution of the orbital constants of motion, allowing eccentric and inclined trajectories to be followed over long inspiral timescales. We focus in particular on the treatment of the Carter constant, whose evolution is not directly determined by the leading-order fluxes for energy and axial angular momentum used in the present model. We investigate several prescriptions for determining its evolution and find that a radial-potential-based prescription provides the most stable behavior among those considered, while also exhibiting a systematic residual bias for nearly planar orbits. We compare the resulting inspirals with established radiation-reaction approximations, including the weak-field treatment of Peters and Kerr-based prescriptions from the literature, and use these comparisons to assess the qualitative behavior and limitations of the implementation. The modular structure of EMRICode allows the radiation-reaction prescription to be replaced independently of the underlying geodesic evolution, providing a framework for incorporating more accurate flux models and additional physical effects. We present two simple examples of how interactions with accretion disks and partial tidal disruptions can drive the evolution of a stellar EMRI orbit. EMRICode therefore provides a flexible computational framework for studying the orbital dynamics of EMRIs and for exploring extensions beyond the idealized vacuum, point-particle approximation.
High-cadence GW astrometry in the $\mu$Hz gap
oai:arXiv.org:2610.27221v1
arXiv:2610.27221v1 Announce Type: cross
Abstract: We forecast the sensitivity of a dedicated, fast-spinning astrometric survey to a stochastic gravitational-wave background (SGWB) in the microhertz gap between pulsar timing arrays (PTAs) and LISA. Spinning the satellite faster samples each star more often, raising the Nyquist frequency and opening higher frequencies at no cost in noise: the shorter, noisier exposures are exactly compensated by their number, so the white-noise floor of each position stream is independent of cadence. High frequencies are where astrometry is most competitive, because the deflection measures the strain directly: its strain noise is white, whereas the equivalent noise of pulsar timing grows as $f^2$. With a cross-correlation estimator over $N$ reference stars and the astrometric overlap computed in the distant-source limit, we find that a Gaia-like satellite scanning the ecliptic poles, with $N\sim10^4$ bright stars over $10$~yr, falls short of the background measured by PTAs by a factor $\sim\!60$. A larger, $8$-m-class aperture, or a constellation of satellites covering the whole sky, reaches it, while extending the reach continuously across the $\mu$Hz band.
P-ONE-1: the first detector line of the Pacific Ocean Neutrino Experiment
oai:arXiv.org:2610.27415v1
arXiv:2610.27415v1 Announce Type: cross
Abstract: Neutrino telescopes instrument optically transparent media with three-dimensional arrays of photosensors to detect Cherenkov radiation induced by neutrino interactions within or near the instrumented volume. The small cross-sections of the weak interaction and the low neutrino fluxes above TeV energies require vast detector volumes, which only natural media on Earth can realistically provide. Spurred by the detection of the first neutrino sources by the IceCube Neutrino Observatory, a new generation of neutrino telescopes is emerging to increase the global instrumented volume and detection sensitivity. The Pacific Ocean Neutrino Experiment (P-ONE) aims to install such a telescope in the Northeast Pacific Ocean, building on established neutrino telescope design targets in collaboration with Ocean Networks Canada (ONC). This article presents the design, integration, and testing of the first P-ONE line, P-ONE-1, in preparation for its deployment and integration into ONC's deep-sea observatory NEPTUNE. Alongside the detector design and integration methods, we detail in-air and in-water system-level validation performed at several test sites in Canada, culminating in the successful construction of P-ONE-1 and its preparation for deployment.
Stability Analysis of Axion Dark Energy Model
oai:arXiv.org:2610.27514v1
arXiv:2610.27514v1 Announce Type: cross
Abstract: We investigate the stability of an axion dark energy model characterized by a periodic scalar-field potential motivated by pseudo Nambu Goldstone bosons. While such models are well suited to drive late-time cosmic acceleration and can exhibit rich dynamical features, including oscillatory behavior of the dark-energy equation-of-state parameter, their viability crucially depends on the existence of stable cosmological solutions. By reformulating the background evolution equations as an autonomous dynamical system, we identify the fixed points corresponding to radiation-dominated, matter-dominated, and dark energy dominated epochs. Using both analytical and numerical techniques, we analyze the stability properties of these fixed points within the observationally allowed region of parameter space. Our results show that the radiation- and matter-dominated phases emerge naturally as saddle points, ensuring a standard cosmological evolution, while the late-time dark energy dominated solution acts as a stable attractor. This guarantees the robustness of cosmic acceleration against variations in initial conditions. Overall, our study establishes the axion dark energy model as a dynamically consistent and observationally viable candidate for explaining the late-time accelerated expansion of the Universe.
Untangling Extreme Mass Ratio Inspirals and Massive Black Hole Binaries in LISA Data
oai:arXiv.org:2610.27642v1
arXiv:2610.27642v1 Announce Type: cross
Abstract: The Laser Interferometer Space Antenna (LISA) will observe sources in the millihertz regime, such as Galactic Binaries (GBs), Massive Black Hole Binaries (MBHBs), and Extreme Mass Ratio Inspirals (EMRIs), so its data stream will be a superposition of an unknown number of gravitational-wave sources. We consider the joint inference of various representative combinations of MBHB and EMRI signals in LISA data (with assumed instrumental noise and stochastic GB foreground). Our Bayesian analyses confirm the expectation of weak inter-source coupling between EMRIs, and between EMRIs and MBHBs, despite strong intra-source parameter correlations. Weak EMRI-MBHB coupling enables a hierarchically sequential approach to EMRI and MBHB inference, rather than fully joint inference. We further show that imperfect MBHB subtraction, arising from missing modes or overlapping MBHBs, has only a weak impact on EMRI inference. For multiple MBHBs, we mitigate inference errors through a joint MBHB fit. We demonstrate this strategy on simulated LISA data containing noise, three MBHBs, and three EMRIs, jointly fitting and subtracting the MBHBs before sequentially inferring the EMRIs - ultimately observing little to no bias on all recovered astrophysical parameters. The notions described throughout our work are highly relevant to informing optimal strategies for the global-fit-based pipelines that are currently being developed for LISA.
Axion detection in the non-Markovian noise
oai:arXiv.org:2610.27731v1
arXiv:2610.27731v1 Announce Type: cross
Abstract: The detection of axion dark matter poses a significant challenge due to its exceptionally weak coupling to standard-model particles. Here we develop a noise-tolerant quantum-metrology scheme based on coupled $^{129}$Xe spin ensembles to search for such weakly interacting particles. Within an experimentally accessible parameter regime, our scheme achieves an ideal exclusion limit on the axion coupling constant that is several orders of magnitude more stringent than that of conventional detection methods. We further investigate the performance of our scheme under dissipative noise and find that its robustness is fundamentally governed by the energy-spectrum structure of the composite spin-environment system. Remarkably, when a bound state emerges in the spectrum, the exclusion limit asymptotically recovers its ideal value in the non-Markovian dynamics, whereas under the Markovian approximation the limit deteriorates significantly due to noise. Our results establish a universal mechanism for noise-resilient exploration of beyond-standard-model physics, which can be practically implemented by engineering bound-state formation via tuning the Larmor frequency of the spin ensembles.
Two-Component Dark Matter and Strongly First-Order Electroweak Phase Transition via a Two-Step Mechanism in the Inert Doublet plus Real Singlet Model
oai:arXiv.org:2610.28014v1
arXiv:2610.28014v1 Announce Type: cross
Abstract: We study a two-component dark matter model in which a real singlet scalar and an inert doublet scalar are simultaneously present alongside the Standard Model, each stabilized by an independent discrete symmetry. The real singlet undergoes a strong first-order transition at high temperature, establishing an intermediate phase with a non-zero singlet vacuum expectation value, which in turn generates the barrier required for the Higgs field to complete a strongly first-order electroweak transition at a lower temperature. This two-step mechanism naturally enforces the condition $v_c/T_c \geq 1$, the sphaleron washout-avoidance criterion necessary for electroweak baryogenesis. The gravitational wave spectrum from both steps is computed; the second step produces a decihertz signal, peaking near $0.08$--$0.34~\mathrm{Hz}$ across the viable parameter space and $0.14$--$0.22~\mathrm{Hz}$ at the benchmark points, that reaches the sensitivity of DECIGO for the strongest transitions and of BBO more widely, while the first remains sub-threshold for all planned detectors. At the coupling values required to drive the phase transition, the singlet is generically underabundant as a thermal relic, so the inert doublet carries the dominant fraction of the observed relic density. Confronting the model with current constraints, we find that the LUX-ZEPLIN 2024 result excludes the entire electroweak-phase-transition-viable parameter space of the singlet sector across the full mass range, even after rescaling by the singlet relic fraction, establishing a complete no-go for the singlet as a viable thermal dark matter candidate in this framework.
Direct Urca process facilitated by non-Fermi-liquid quantum criticality
oai:arXiv.org:2610.28023v1
arXiv:2610.28023v1 Announce Type: cross
Abstract: Conventional Fermi-liquid theories of dense nuclear matter predict that the direct Urca process requires a proton fraction above a threshold of $Y_p \approx 0.11$--$0.15$, which is rarely reached in low-mass neutron stars. Here we show that this stringent kinematic constraint is drastically relaxed by the non-Fermi liquid behavior arising from superfluid and superconducting quantum criticality in neutron stars. In the non-Fermi-liquid regimes, nucleons have a smooth momentum occupation without a sharp Fermi surface, and their chemical potentials develop a logarithmic temperature dependence. These intrinsic features fundamentally alter the standard $\beta$-equilibrium relation. Consequently, direct Urca neutrino emission can be activated even at low nuclear densities. This modifies the long-standing Urca-threshold paradigm and indicates that direct Urca process has the opportunity to occur in low-mass neutron stars.
On the Stability of Tolman Oppenheimer Volkoff Fireholes
oai:arXiv.org:2610.28140v1
arXiv:2610.28140v1 Announce Type: cross
Abstract: We investigate the radial stability of singular Tolman--Oppenheimer--Volkoff fireholes with a radiation equation of state. These configurations contain a naked, negative point-mass singularity surrounded by a positive-density fluid, whose mass profile develops a Schwarzschild-like plateau before approaching an extended, self-similar atmosphere. Using the Pr\"ufer phase method to identify the fundamental radial mode, we consider two outer boundary conditions: vanishing Lagrangian pressure perturbation (open firehole) and vanishing radial displacement (boxed firehole). For both, the fundamental squared frequency changes from positive to negative as the outer boundary exceeds a critical radius $R_c$. We find that $R_c$ follows approximately the same scaling with the central mass parameter as the equilibrium firehole radius $R$. However, the stability thresholds differ significantly: $R_c^{\rm O}
Hybrid Stars with Different Amount of Quark Matter
oai:arXiv.org:2610.28408v1
arXiv:2610.28408v1 Announce Type: cross
Abstract: In a previous work, we explored the possibility of neutron stars with large quark cores. Through Bayesian inference, we obtained large sets of equations of state imposing constraints regarding the properties of nuclear matter, mass-radius data from NICER, the GW170817 data from Virgo-LIGO collaboration, perturbative QCD and causality. We also imposed that the maximum mass neutron star (NS$\text{max}$) must have a large quark core. In this work, we obtain a new set of equations of state (SmallQ set) with the same constraints as before. However, we relaxed the last constraint, imposing only that NS$\text{max}$ must have quark core, regardless of its size. We compare the results with the set from our previous paper constrained by NICER and GW170817 observations NJL-GW. For both sets, we use the Relativistic Mean Field model with nonlinear meson interactions for the hadron phase and the Nambu--Jona-Lasinio model with four- and eight-quark interaction terms for the quark phase. The SmallQ set has higher values of phase transition density, with 3 -- 4.5 $\rho_0$, compared to the NJL-GW, with 2 -- 2.5 $\rho_0$. The new set describes NSs with smaller radius and smaller tidal deformability, and is more compatible with NICER and GW170817 data. A difference in the quark matter composition is observed, with a larger quark strange fraction occurring in the SmallQ set. Within this set, $2 M_\odot$ NSs may contain little to no quark matter. It is shown that there is a non negligible fraction of EOS whose mass-radius curve presents a positive slope for 1.8 $M_\odot$ stars, while for nucleonic stars this fraction is negligible and for hyperonic stars the slope is always negative.
Exponential Growth of Turbulent Cross Helicity in Shear-Flow-Driven Turbulence
oai:arXiv.org:2610.28424v1
arXiv:2610.28424v1 Announce Type: cross
Abstract: Cross helicity, the velocity--magnetic-field correlation, is an ideal invariant of incompressible magnetohydrodynamics and constrains energy cascades and magnetic-field generation. Although cross helicity typically decays via visco-resistive effects, we analytically demonstrate that it grows exponentially in sheared inhomogeneous turbulence. This cross-helical shear instability is confirmed using three-dimensional numerical simulations. Growth requires a nonzero gradient of the Reynolds--Maxwell stress to produce turbulent cross helicity; this production overcomes effective losses when shear exceeds a threshold set by the ratio of turbulent magnetic to momentum diffusivities.
MadStree: automated evaluation of tree-level de Sitter correlators
oai:arXiv.org:2610.28526v1
arXiv:2610.28526v1 Announce Type: cross
Abstract: We first derived the formulae of recursive reduction and dlog-form differential equations for arbitrary tree-level cosmological correlator, taking the external relations of massless bulk propagators into consideration. The differential equations could be systematically solved and provide the efficient evaluation of them. We have implemented this methodology as a user-friendly Mathematica package MadStree for computing arbitrary tree-level cosmological correlators.
Lagrangian perturbation theory for modified gravity
oai:arXiv.org:1705.10719v4
arXiv:1705.10719v4 Announce Type: replace
Abstract: We present a formalism to compute Lagrangian displacement fields for a wide range of cosmologies in the context of perturbation theory up to third order. We emphasize the case of theories with scale dependent gravitational strengths, such as chameleons, but our formalism can be accommodated to other modified gravity theories. In the non-linear regime two qualitative features arise. One, as is well known, is that nonlinearities lead to a screening of the force mediated by the scalar field. The second is a consequence of the transformation of the Klein-Gordon equation from Eulerian to Lagrangian coordinates, producing frame-lagging terms that are important especially at large scales, and if not considered, the theory does not reduce to the $\Lambda$CDM model in that limit. We apply our formalism to compute the 1-loop power spectrum and the correlation function in $f(R)$ gravity by using different resummation schemes. We further discuss the IR divergences of these formalisms.
Don't miss the forest for the trees: the Lyman alpha forest power spectrum in effective field theory
oai:arXiv.org:2309.10133v3
arXiv:2309.10133v3 Announce Type: replace
Abstract: We derive an effective field theory (EFT) for cosmological Lyman alpha forest fluctuations valid for the power spectrum at the one-loop order. The ``bottom-up'' EFT expansion at the level of the transmitted flux is identical to the line-of-sight dependent bias model first derived by Desjacques et al. We confirm this result by a ``top-down'' derivation based on the exponential map of the optical depth field. Specifically, we show that the combination of the exponential map and conditions of renormalizability generates the same EFT expansion as the ``bottom-up'' approach. In passing, we point out inconsistencies of the tree-level perturbative expansion of the exponential map without counterterms. To facilitate practical applications, we generalize the FFTLog method for efficient calculations of one-loop integrals from new line-of-sight dependent operators. Finally, we compare the one-loop EFT model against data from the Sherwood hydrodynamic simulation. The theory fits the data with sub percent accuracy up to $k= $ 3 $h$Mpc$^{-1}$ at $z= 2.8$ for both 3D and 1D correlations. Our model can be readily used for cosmological full-shape analyses of Lyman alpha forest data.
Unveiling dark forces with measurements of the Large Scale Structure of the Universe
oai:arXiv.org:2309.11496v3
arXiv:2309.11496v3 Announce Type: replace
Abstract: Cosmology offers opportunities to test Dark Matter independently of its interactions with the Standard Model. We study the imprints of long-range forces acting solely in the dark sector on the distribution of galaxies, the so-called Large Scale Structure (LSS). We derive the strongest constraint on such forces from a combination of Planck and BOSS data. Along the way we consistently develop, for the first time, the Effective Field Theory of LSS in the presence of new dynamics in the dark sector. We forecast that future surveys will improve the current bound by an order of magnitude.
Causality and perturbative-QCD constraints on first-order phase transitions in neutron stars
oai:arXiv.org:2410.21877v3
arXiv:2410.21877v3 Announce Type: replace
Abstract: Causality bounds the compactness and maximum mass of neutron stars. We generalise this bound to equations of state (EOSs) with a first-order phase transition between maximally stiff phases. In units of the transition energy density, the problem depends only on $q=\ptr/\etr$ and $j=\Delta\ve/\etr$. The maximum compactness equals the smooth value $\mathcal{C}=0.354$ for $q\to0$ or $j\to0$, is unaffected for $q\gtrsim0.67$, and is otherwise reduced by typically $\lesssim0.04$ when the maximum-mass star is a hybrid star and by up to $0.13$ when it is the last hadronic star. Matching to a Baym--Pethick--Sutherland crust and a chiral-EFT polytrope at $n_m=1.1\,n_0$ and inverting the bound, no causal EOS supports $2.01\,\Msun$ if the energy-density jump exceeds $\simeq460$--$495~{\rm MeV\,fm^{-3}}$. Additionally, imposing radius constraint of $R\ge11.6$~km it shrinks to $\simeq235$--$395~{\rm MeV\,fm^{-3}}$. These limits hold whatever the matter above the transition, for onsets within $\simeq0.2\,n_0$ of $n_m$. Relaxing maximal stiffness, we build EOSs with a transition at $\mu_{\rm tr}=1000$--$1300$~MeV, matched to perturbative QCD (pQCD) at $2.6$~GeV. pQCD limits the jump; EOSs compatible with PSR J0740+6620, J1614$-$2230, J0437$-$4715, J0614$-$3329, GW170817, and XTE J1814-338 appear only at $\mu_{\rm tr}=1200$ and $1250$~MeV; and inside the causal window, the largest jump agrees with the bound within the grid resolution. Finally, inserting a transition at $2.01\,\Msun$ into hadronic EOSs that satisfy five pulsar observations, the EOS with the lowest break density gives a narrow stable hybrid window ($\Delta M\simeq0.005\,\Msun$), compatible with XTE J1814-338 for jumps of $6$--$8\,n_0$ and with pQCD; three EOSs with higher break densities give none.
Kratos-polrad: Novel GPU system for Monte-Carlo simulations with consistent polarization calculations
oai:arXiv.org:2512.01283v2
arXiv:2512.01283v2 Announce Type: replace
Abstract: Polarized radiation serves as a vital diagnostic tool in astrophysics, providing unique insights into magnetic field geometries, scattering processes, and three-dimensional structures in diverse astrophysical scenarios. To address these applications, we present Kratos_polrad, a novel GPU-accelerated Monte Carlo Radiative Transfer code built upon the heterogeneous computing framework of Kratos, designed for self-consistent and efficient polarization calculations. It utilizes a comprehensive treatment of the Stokes parameters throughout photon propagation, featuring grain-lab frame transformations using quaternion algebra and a consistent analytical integration of polarized extinction and emission within each cell, which are useful in modeling radiative transfer processes with scatterings by aligned dust grains. The code implements two-step polarimetry imaging that decouples Monte Carlo sampling of scattering physics from imaging geometry, enabling efficient synthesis maximizing the utilization of photon packets. Extensive validation against analytical solutions and established codes demonstrates accurate treatment of diverse polarization phenomena, including self-scattering polarization, dichroic extinction in aligned dust grains, and complex polarization patterns in twisted magnetic field configurations. By leveraging massive GPU parallelism, optimized memory access patterns, and analytical approaches for optically thick cells, Kratos_polrad achieves performance improvements of ~10^2 times compared to CPU-based methods, enabling previously prohibitive studies in polarimetric astrophysics.
MAMMOTH-LyC: Investigating the Role of Galaxy Mergers in a Strong Lyman Continuum Leaker at $z=2.39$
oai:arXiv.org:2512.23202v3
arXiv:2512.23202v3 Announce Type: replace
Abstract: The MAMMOTH-LyC survey is a cycle 30 Hubble Space Telescope (HST) medium program obtaining 18-orbit-deep WFC3/UVIS F225W imaging in two massive galaxy protocluster fields at $z\sim2.2$. We introduce this survey by reporting the discovery of J1244-LyC1, a strong Lyman continuum (LyC) leaker at $z = 2.39$, exhibiting clear merger signatures. J1244-LyC1 has a highly significant ($6\sigma$) LyC detection, corresponding to an absolute escape fraction of $f_{\mathrm{esc}} \! =\!36\%\pm6\%$ ($1\sigma$). The LyC emission is spatially resolved into multiple peaks that coincide with the system's disturbed morphology, confirming genuine multi-site LyC leakage. With a stellar mass of $10^{10.2}{M_\odot}$, J1244-LyC1 is both the first confirmed high-redshift LyC-leaking merger and the most massive LyC emitter known to date. We interpret J1244-LyC1 as a merger-driven starburst system in which tidal interactions have disrupted the interstellar medium, creating multiple low-column-density pathways that facilitate LyC escape. This discovery provides the first direct evidence of spatially resolved LyC escape in a merging system, offering new insight into the potential role of major mergers in driving the cosmic reionization.
Formation and Evolution of [Wolf-Rayet] Planetary Nebulae through a Late Thermal Pulse
oai:arXiv.org:2603.26536v2
arXiv:2603.26536v2 Announce Type: replace
Abstract: We present the first radiation-hydrodynamical simulations of the formation of a born-again planetary nebula (PN) triggered by a late thermal pulse (LTP). The 2D radiation-hydrodynamic simulations, performed with the {\sc pluto} code, have been consistently coupled to stellar evolution calculations using the Modules for Experiments in Stellar Astrophysics ({\sc mesa}) code. Very particularly the stellar evolution model uses (i) updated opacity tables for H-deficient, C-rich mixtures during the LTP, and (ii) a mass-loss prescription tailored for H-deficient [Wolf-Rayet]([WR])-type winds during the post-LTP phase. Our stellar model reproduces the nearly complete depletion of H expected after an LTP event, while matching the observed abundances and spectral types of iconic [WR]-type central stars of PNe. The simulations show for the first time that the H-deficient LTP ejecta forms a transient double-shell structure which, after $\sim$1000 yr, becomes fully mixed with the H-rich PN. The ejecta mass ($\sim3.4\times10^{-4}$~M$_\odot$) is too small to leave a lasting imprint on the nebular abundances, predicting H-rich PNe around [WR] central stars. The injection of LTP material into the hot bubble drives turbulence, clump formation, and enhanced mixing, providing an explanation to the larger expansion velocities and larger turbulent nebular structures of PNe with [WR] central stars compared to those with H-rich central stars. These results provide robust support for the born-again scenario as the origin of H-deficient [WR] central stars within H-rich PNe.
An empirical derivation of the average dust attenuation law at 2oai:arXiv.org:2604.09763v2
arXiv:2604.09763v2 Announce Type: replace
Abstract: Dust attenuation affects galaxy SEDs and introduces uncertainties in star formation rates, stellar masses, and metallicities. Direct spectroscopic constraints at early cosmic epochs were limited before JWST. We derive the average dust attenuation law of massive (log10(Mstar/Msun) > 9) star-forming galaxies at 2 < z < 7 using JWST spectroscopy. We combine JWST/NIRSpec spectroscopy from JADES with deep multi-wavelength photometry from the ASTRODEEP-JWST catalogues. For approximately 120 mass-selected galaxies with reliable Balmer decrement (Halpha/Hbeta) measurements, we stack SEDs in bins of Balmer optical depth and derive the selective attenuation curve following Calzetti et al. (1994). JWST/MIRI photometry extends the wavelength coverage into the near-infrared (NIR). The resulting curve spans 0.16-1.14 microns in the rest frame and is well described by a smooth function. We derive R_V approximately 3.98, with Delta R_V = 0.16 indicating its sensitivity to the adopted long-wavelength extrapolation. The average attenuation law agrees remarkably well with the local starburst relation in slope and normalisation, but appears systematically flatter in the ultraviolet than several intermediate-redshift determinations. We find no significant evidence for a 2175 Angstrom UV bump. Our results provide a direct empirical constraint on the average dust attenuation law of massive star-forming galaxies at 2 < z < 7. This analysis exploits the combined spectroscopic and photometric reach of JWST. Despite variations among individual galaxies, the ensemble-average behaviour agrees with the local starburst relation, suggesting that the main mechanisms governing dust attenuation were already operating at early cosmic times. The relatively shallow UV slope may reflect differences in dust-star geometry, radiative-transfer effects, or effective grain-size distributions among galaxy populations.
Silicate cosmic dust grain collisions in the interstellar medium: A molecular dynamics study
oai:arXiv.org:2605.22225v3
arXiv:2605.22225v3 Announce Type: replace
Abstract: (abridged) We aim to predict the most important parameters for grain-grain collision outcomes in models of interstellar grain population evolution: the threshold velocities for shattering and vaporization, and the resulting distributions of grain sizes. We use molecular dynamics simulations which evolve the dynamics of each atom in a dust grain to explore the outcomes of collisions between silicate grains of radii $a \in [5,50]~\AA$ at velocities 0.1-20 km/s. We run simulations of grains with two materials: amorphous SiO$_2$ and an amorphous silicate of composition suggested by Draine \& Hensley (2021). With these simulations, we quantify the collision velocity dependence of shattered and vaporized mass fractions, and the resulting size distributions of shattering products. We find grain shattering thresholds are $\sim$6 km/s for both materials, a factor of $\sim$2 higher than the canonical value for silicates of 2.7 km/s from Tielens et al. (1994) and Jones et al. (1996). The size distributions of shattered products are generally not consistent with the power law distributions predicted by Jones et al. (1996), but the physical reasons for these discrepancies remain uncertain. The Tielens et al. (1994) expressions fail to predict the fraction of shattered or vaporized material observed in our numerical simulations, over-predicting the shattered fraction and under-predicting the vaporized fraction, except for the most massive grains we simulate, where the predicted vaporized fraction agrees reasonably well. The model of Hirashita \& Kobayashi (2013) similarly fails to match the simulations. We provide updated shattering velocity thresholds for nanoscale silicate grains, which suggest that the smallest astrophysical dust grains composed of silicate materials may be more robust to shattering in the interstellar medium than previously assumed.
The Life and Death of Stars That Capture Primordial Black Holes
oai:arXiv.org:2606.02700v2
arXiv:2606.02700v2 Announce Type: replace
Abstract: Primordial black holes (PBHs) in the asteroid mass window ($10^{17}-10^{23}\,{\rm g}$) remain viable dark matter candidates and can be captured by stars. We develop the first global framework for the evolution of stars that capture PBHs, combining analytic calculations, stellar evolution models, 3D general-relativistic magnetohydrodynamic simulations, and Monte Carlo population synthesis. We find that the fate of these systems bifurcates: PBHs that form an accretion disk before consuming the host drive explosive disruption, whereas PBHs captured too late or growing too slowly consume the star quietly. Capture is dominated by three-body interactions with planetary or stellar companions. For a solar-type host with a Jupiter analog, inspiral within a main-sequence lifetime requires $M_{\rm BH}^{\rm crit}\gtrsim 10^{22}\,{\rm g}$, while lighter PBHs generally require tighter companions. Once deposited at the center, the PBH grows through inefficient quasi-spherical Bondi accretion; if it reaches the angular-momentum threshold before consuming the host, the inflow circularizes into a disk. Our Monte Carlo calculations yield sizable quiet-consumption and explosive-disruption populations, with final PBH masses $M_{\rm BH}\sim0.01-1\,M_\odot$ and disk-forming PBH spins $a_\ast\approx0.8$. Disk formation is the point of no return: disk winds and relativistic jets of $\sim10^{45}-10^{50}\,{\rm erg\,s^{-1}}$ disrupt the star within minutes. The resulting transients may include a $\sim$day-long UV/blue signal, radio afterglow, and, if the jet escapes, an X-ray-flash/low-luminosity gamma-ray-burst (XRF/llGRB) signal. For an $O(1)$ PBH dark matter fraction and optimistic capture assumptions, the event rate can reach that of llGRBs. The low-mass, high-spin remnants offer a complementary PBH probe and possible source for subsolar BH mergers.
The Effect of Adiabatic Index on Radius Evolution and the Mass Loss
oai:arXiv.org:2606.17143v2
arXiv:2606.17143v2 Announce Type: replace
Abstract: Models that track the size evolution of exoplanets often assume a prescribed initial thermal state or a single adiabatic index to describe the planetary interior structure, the latter of which is taken to be $\gamma \sim1.4$ which is likely appropriate for evolved planets ($\gtrsim$1 Gyr). Extrapolating this high $\gamma$ to earlier ages (down to $\sim$million years old) is problematic since, according to thermodynamics, the adiabatic index of young planets is $\sim$1.2, which is low enough to drastically change how interior mass is distributed. We quantify the effect of varying the adiabatic index from 1.2 to 1.4 on the expected radius of the exoplanet over time. We find that envelopes of larger adiabatic indices start puffier with all else equal and undergo faster radius contraction with accelerated mass loss. Assumption of high $\gamma$ can therefore overestimate the effect of mass loss in shaping the exoplanetary population, especially when young planets are considered. We highlight the need for a more careful consideration of the initial thermal condition of planets in evolutionary models to properly interpret the radii measurements of exoplanets.
A New Scaling of Neutron Star Tidal Deformability for Directly Probing the Central Equation of State Parameter
oai:arXiv.org:2606.21402v2
arXiv:2606.21402v2 Announce Type: replace
Abstract: The dimensionless tidal deformability, $\Lambda$, of neutron stars (NSs), inferred from gravitational-wave (GW) observations, has thus far been used primarily to constrain the pressure of dense matter near twice nuclear saturation density, leaving the core equation of state (EOS) largely inaccessible to inspiral-phase GW observations. We show that a fundamental central EOS parameter can be probed directly through $\Lambda$ using a perturbative analysis of the dimensionless stellar-structure and tidal-response equations formulated in terms of scaled intrinsic variables, without invoking any specific EOS model. We uncover a remarkable EOS-insensitive scaling relation between $\Lambda$ and the central EOS parameter $\rm{X}\equiv P_{\rm c}/\varepsilon_{\rm c}$, where $P_{\rm c}$ and $\varepsilon_{\rm c}$ denote the central pressure and energy density, respectively. The relation is validated against a broad ensemble of physically viable EOSs. Applying it to tidal deformabilities inferred from events such as GW170817 enables a direct determination of $\rm{X}$. By adopting the previously established effective upper bound $\rm{X}\lesssim0.381$, we infer a tight lower bound, $\Lambda_{\rm{TOV}}\gtrsim9.1^{+0.5}_{-0.5}$, for maximum-mass NSs along stable mass-radius sequences, quantitatively demonstrating that even the most compact stable NSs remain distinctly separated from black holes, for which $\Lambda_{\rm{BH}}=0$. These findings reveal a previously unrecognized connection between inspiral-phase tidal deformability and the central EOS parameter $\rm{X}$, establishing a direct link between GW observables and the microphysics of ultradense matter in the strong-gravity regime. The resulting scaling establishes inspiral-phase tidal deformability as a direct and largely model-insensitive probe of the central EOS parameter $\rm{X}$, and hence of key properties of NS cores.
Search for Diffuse Supernova Neutrino Background in the Full KamLAND Dataset with Neural-Network-Based Event Classification
oai:arXiv.org:2606.29381v2
arXiv:2606.29381v2 Announce Type: replace
Abstract: We report a search for the diffuse supernova neutrino background (DSNB) with the KamLAND detector, targeting electron antineutrinos via inverse beta decay in the neutrino energy range of 8.3 to 30.8 MeV. Using liquid-scintillator exposures of 9.02 kton-year for 8.3 to 9.3 MeV and 9.42 kton-year for 9.3 to 30.8 MeV, we observe seven candidate events after applying a new deep-neural-network-based event classification technique. This result is consistent with the background-only expectation of 16.2 plus or minus 9.4 events, including systematic uncertainties associated with the neural-network selection. A spectral analysis of the energy and radial distributions finds no significant excess attributable to the DSNB. We therefore set 90 percent confidence-level upper limits on the DSNB flux of 38 to 43 per square centimeter per second, depending on the assumed DSNB model. We also derive model-independent 90 percent confidence-level upper limits on the electron-antineutrino flux, obtaining some of the most stringent constraints below 13.3 MeV. Beyond the DSNB search itself, this work demonstrates neural-network-based event classification as a promising approach for suppressing neutron-associated backgrounds in liquid-scintillator neutrino detectors.
Unveiling the Milky Way with a Gaia DR3 census of OB-type stars within 2 kpc. I. Tracing local Galactic structure, massive star-forming regions and core-collapse supernova progenitors
oai:arXiv.org:2607.07068v2
arXiv:2607.07068v2 Announce Type: replace
Abstract: O- and B-type stars are young and hot, thereby serving as vital tracers of the star formation and spiral arm structure of the Milky Way. At the dusk of the \textit{Gaia} DR3 era, a high-confidence and accurate catalogue appears timely. We have characterized 105,971 OB-type stars (T$_{\rm eff} >$ 10,000 K) within 2 kpc from the Sun, using an astro-photometric Bayesian inference tool. Our map unveils a complex view of the young stellar populations across the thin disk, with prominent large-scale features such as the Cepheus Spur, the Giant Oval Cavity, and a segment the Sagittarius-Carina arm all visible. Their inhomogeneous spatial distribution implies that massive star formation has taken place clustered across a few highly concentrated regions. We find a correlation between the overdensities of OB-type stars and young open clusters, and attribute departure from it to extinction effects. We identify 4200 stars as model-dependent core-collapse supernova (ccSN) or direct-collapse black hole (BH) progenitor candidates, including 6 ccSN progenitor candidates within 125 pc -- the maximum radius within which a ccSN explosion could be harmful to Earth's biosphere -- but our analysis suggests none of them exploding within the next 5 Myrs. Finally, we anticipate more evolved OB-type stars to collapse to form BHs within the next 1 Myr than ccSN progenitors will explode, despite the former's much scarcer number. A possible interpretation is a recent star formation burst in the local Milky Way, consistent with overdensities of ccSN progenitor candidates in Cyg OB2 and the Sagittarius-Carina arm.
Infrasound Newtonian Noise Estimation at the Einstein Telescope Candidate Site Sos Enattos
oai:arXiv.org:2607.21154v2
arXiv:2607.21154v2 Announce Type: replace
Abstract: We investigate the seasonal variability of atmospheric infrasound and its contribution to Newtonian noise (NN) at the Sos Enattos site, a leading candidate location for the Einstein Telescope (ET). Infrasound data recorded at three stations, SOE0 (surface), SOE1 ($-84$ m) and SOE3 ($-160$ m), are analysed over multiple seasons to characterise both the temporal variability and the depth dependence of the acoustic field. At the surface, the median amplitude spectral density (ASD) at 1 Hz shows no consistent seasonal modulation, and its variability is driven primarily by wind; the apparent dependence on sea state largely reflects its correlation with wind speed. Underground, where local wind noise is suppressed, a winter maximum of 4--7 dB above the quietest season is observed in both years of data. Using the measured sound-pressure spectra, we estimate the corresponding NN contribution for an underground cavern with the dimensions foreseen for ET, assuming an isotropic field of acoustic plane waves. The surface pressure field, acting through the atmosphere above a test mass at 200 m depth, exceeds the ET-D design sensitivity in median conditions between 1.6 and 7.1 Hz, by up to a factor of 4.7. Inside the cavern, the median NN remains a factor of 1.4 below ET-D with the pressure measured at SOE1, but exceeds it between 2.5 and 8.7 Hz, by up to a factor of 3.5, with the pressure measured at SOE3. The attenuation is not monotonic with depth: the deeper station SOE3 shows higher levels than SOE1 above a few hertz, pointing to local, anthropogenic sources in the underground infrastructure. These results show that underground placement is essential but not sufficient on its own, and that the acoustic emission of the underground infrastructure and future laboratory equipment will require careful evaluation and mitigation to keep atmospheric NN below the ET design sensitivity.
Observational Evidence for Anisotropic Metal Excess around Galaxies
oai:arXiv.org:2607.24929v2
arXiv:2607.24929v2 Announce Type: replace
Abstract: The exchange of matter and energy between galaxies and their surroundings drives the cosmic baryon cycle, yet mapping metal transport remains an observational challenge. While simulations predict that galactic winds escape anisotropically along minor axes, evidence for chemical enrichment in neighboring galaxies is limited. We analyze 1,433 galaxy pairs from the Dark Energy Spectroscopic Instrument survey and detect a gas-phase metallicity excess of 14.6% $\pm$ 3.7% to 24.2% $\pm$ 2.6% in neighbors aligned with the minor axis of massive primary at projected separations of 15--60 kpc, qualitatively consistent with IllustrisTNG simulation. Permutation tests indicate that this excess varies from a marginal detection ($>$92% confidence) at closer separations (15--30 kpc) to a significant signal ($>$98% confidence) at intermediate separations (30--60 kpc). Here, we show that these results provide empirical constraints on galactic feedback models and are consistent with a scenario where anisotropic outflows can enrich nearby galaxies, a mechanism that complements other environmental processes.
VAPOLA - A multiyear, multiband polarization survey of AGNs and Sgr A* at millimeter wavelengths with ALMA II. Spectropolarimetric properties and their evolution from 2017 to 2023
oai:arXiv.org:2608.14900v2
arXiv:2608.14900v2 Announce Type: replace
Abstract: We present a systematic analysis of the spectropolarimetric properties of a sample of 39 active galactic nuclei and Sagittarius A* observed with the Atacama Large Millimeter/submillimeter Array during five VLBI campaigns between 2017 and 2023. We characterize the compact cores in total intensity and polarization, focusing on the behavior or the linear polarization fraction (LP), electric vector position angle (EVPA), and Faraday rotation measure (RM) over time and spectral domains. We investigate both individual objects--such as M87, Sgr A*, 3C273, and 3C279--and ensemble properties of different source classes, including flat-spectrum radio quasars, BLLac objects, and other active galaxies. While total intensity and spectral index are generally stable on weekly timescales, polarization properties often exhibit strong variability, with significant day-to-day changes in LP, EVPA, and RM. Several sources display large EVPA rotations accompanied by variations in LP and RM, in some cases coinciding with flaring activity. We observe that the magnitude of RM increases with observing frequency for all sources for which we have reliable multi-band measurements, consistent with Faraday rotation arising in a magnetized sheath surrounding the relativistic jet, although an origin in the accretion flow--particularly in the case of Sgr A*--cannot be ruled out.
Periodic activity on the (sub)mm surface of the AGB star Mira A
oai:arXiv.org:2608.23724v2
arXiv:2608.23724v2 Announce Type: replace
Abstract: Evolved stars on the asymptotic giant branch (AGB) lose mass through stellar winds launched from their dynamic extended atmospheres. Recent high-angular-resolution observations of the archetypal AGB star Mira A (o Ceti) have revealed compact hotspots on its (sub)millimetre surface, whose origin remains unclear. We aim to characterise the (sub)millimetre photosphere of Mira A and investigate the properties and temporal behaviour of the compact hotspots detected on its surface. We analysed 19 epochs of high-angular-resolution observations obtained with ALMA in Bands 4, 6, 7, and 9. The complex visibilities were fitted with parametric models to derive the flux density, size, morphology, and brightness temperature of the stellar disc and an additional, more compact component across multiple pulsation phases and observing frequencies. Our results show that, near maximum light, the flux density of the stellar disc increases with frequency while its apparent size decreases, in agreement with the expected stratified radio photosphere. In Band 6, both the flux density and the size of the stellar disc appear to vary with pulsation phase, with the largest measured values occurring close to optical maximum light and the lowest observed around minimum light. Compact hotspots are detected in Bands 6 and 7 near maximum light, with their brightness temperatures exceeding those of the stellar disc by several thousand kelvin and reaching a temperature excess of ~8200 K in Band 6. They are highly energetic and are detected over timescales of weeks, while the available observations also suggest potential longer-term changes in both their position and intensity. Hotspot activity appears to occur repeatedly across multiple pulsation cycles and may be linked to strong shocks, magnetic activity, and/or episodic mass ejections, and could therefore play an important role in the mass-loss process in AGB stars.
Wide Field Localization Using Fixed Modulation Collimators
oai:arXiv.org:2609.10384v2
arXiv:2609.10384v2 Announce Type: replace
Abstract: Rapid, precise localization of gamma-ray bursts (GRBs) requires continuous wide-field imaging at hard X-ray and soft gamma-ray energies. We develop a fixed (non-rotating) modulation-collimator architecture to address this challenge. A bigrid collimator, a pair of aligned, periodic slit grids in front of a simple, non-imaging photon-counting detector, transmits an off-axis source through a triangular fringe pattern whose local angular position, the source phase, is recovered from simultaneous count rates in several grid pairs at fixed relative phase offsets, rather than by mechanically rotating a single grid pair as in the classical rotating modulation collimator. Four-phase demodulation of these count rates robustly cancels detector background and grid leakage without requiring rotation. Because the fringe pattern repeats across the sky, a phase measurement corresponds to many candidate source positions; we resolve this ambiguity over an arbitrarily wide field of view using a cascade of progressively coarser-pitched grid stages, generalizing the two-grid vernier-scale technique originally used to localize Sco X-1. This vernier cascade technique reduces the required photon counts and grid-calibration accuracy required to achieve fine angular resolution over a broad field of view. Generalized to two dimensions with orthogonal cascades, the resulting architecture achieves wide-field localization with no moving parts, no time modulation, and only simple photon-counting detectors, a promising candidate for simple, inexpensive, continuous-monitoring instruments supporting time-domain and multi-messenger astrophysics.
Thermal Evolution of Lava Planets Across System Ages: Predictions for Hell of a Survey
oai:arXiv.org:2609.20155v2
arXiv:2609.20155v2 Announce Type: replace
Abstract: Ultra-short-period (USP) rocky exoplanets can have dayside temperatures high enough to maintain permanent magma oceans, sitting at the intersection of interior geophysics and atmospheric chemistry. Coupled feedbacks between the molten surface and outgassed atmosphere can sustain or enhance a volatile envelope, while stellar interactions can erode it. Understanding which outcome prevails, and its observable imprint, requires a multi-target approach across planets at different stages of thermal evolution. We present predictions for the five targets of JWST Cycle 4 program 8864: TOI-1807 b, TOI-2260 b, TOI-431 b, TOI-6255 b, and TOI-2431 b. Using the PROTEUS coupled interior-atmosphere framework, we construct a simulation grid and classify outcomes into six categories, defined by the final interior melt state and by whether the planet retains a detectable atmosphere thick enough to redistribute heat to the nightside. For targets retaining a non-negligible volatile envelope, our models predict higher partial pressures for most species when the surface is molten, except for S$_{2}$, whose enhancement in the solid regime suggests it may serve as a tracer of interior melt state. Despite some targets showing outcomes across multiple scenarios, most tend toward a bare-rock end-member, with global melt fraction $\leq$ 20\% and atmospheric retention sensitive to escape efficiency. Our analysis reveals a minimum escape efficiency threshold below which volatile envelopes survive under energy-limited escape, constraining the conditions required for atmosphere survival on irradiated rocky planets. These predictions will guide interpretation of MIRI-LRS phase curve observations and identify which targets and features best discriminate between competing geophysical states.
Solar spicule population properties inferred from IRIS Mg II k off-limb spectra using likelihood-free Bayesian inference
oai:arXiv.org:2609.37811v2
arXiv:2609.37811v2 Announce Type: replace
Abstract: Solar spicules are ubiquitous and strongly overlap in off-limb observations, which complicates the interpretation of individual events. We aim to infer population-level properties of spicules from IRIS Mg II k off-limb spectra using a forward model that explicitly accounts for line-of-sight superposition and instrumental effects. We constructed a stochastic 3D spicule-population model and synthesized Mg II k spectra with radiative transfer, including instrumental degradation and photon noise. We then constrained the hyperparameters that describe the underlying spicule population with likelihood-free Bayesian inference. This approach combines approximate Bayesian computation and an adaptive population Monte Carlo sampler, comparing observed and simulated datasets through principal-component summary statistics at five heights above the limb. The inferred posterior constrains effective population-level properties of the spicules that dominate the observed Mg II k variability and height dependence. These constraints provide a route to connect off-limb spectroscopy with forthcoming spectropolarimetric diagnostics of the chromosphere--corona interface.
Galaxy Protoclusters as Drivers of Cosmic Reionization: II. Te-Based Metallicities of Lyman-{\alpha} Emitters
oai:arXiv.org:2610.01811v2
arXiv:2610.01811v2 Announce Type: replace
Abstract: Context. Protoclusters at reionization are thought to host a significant fraction of the cosmic star formation rate density, making them crucial for the timing and topology of reionization. The gas-phase metallicity, a key tracer of galaxy evolutionary stage, can now be measured directly in the early Universe with JWST. Aims. We measure the gas-phase metallicity of Lyman-$\alpha$ emitters (LAEs) in the LAGER-z7OD1 protocluster at $z\approx6.93$ via the direct ($T_e$-based) method using the [O III]$\lambda$4363 auroral line, and combine it with stellar masses to constrain the mass-metallicity relation (MZR) in dense environments. Methods. We use JWST/NIRSpec MSA G395H spectroscopy, derive metallicities with self-consistent analytic prescriptions, and propagate the dust-extinction uncertainty into an asymmetric error budget via Monte Carlo simulations based on the observed Balmer decrements. Stellar masses come from BEAGLE SED fitting of JWST/NIRCam and UltraVISTA photometry. Results. We obtain five secure [O III]$\lambda$4363 detections and one upper limit, with 12+log(O/H) $\approx$ 7.2-8.0. The weighted spectral stack yields 12+log(O/H) = $7.76^{+0.06}_{-0.05}$, ~0.26 dex higher than stacks of LAEs at similar redshifts from JWST surveys. However, the MZR fit gives a slope $\gamma$ = 0.232 +/- 0.214 and normalization $Z_0$ = 7.888 +/- 0.199, placing our sources below the established MZR of normal galaxies and recent protocluster baselines at $z\sim5-7$. Conclusions. The higher metallicity relative to general LAEs suggests accelerated assembly driven by the overdensity, while the offset below the MZR points to a complex interplay between mass build-up and gas accretion. Since we target actively star-forming LAEs, we are likely capturing a phase of pristine gas inflow from the intergalactic medium that fuels star formation but dilutes the interstellar media of these protocluster members.
Collisional excitation of SiC$_6$ by helium: cross sections, rate coefficients, and astrophysical implications
oai:arXiv.org:2610.06525v3
arXiv:2610.06525v3 Announce Type: replace
Abstract: Rotational transitions of SiC$_6$ with $j=26$--$40$ have been detected in IRC+10216, but collisional rate coefficients for SiC$_6$ with He or H$_2$ were previously unavailable. We present rate coefficients for SiC$_6$ rotational excitation by He, calculated using a two-dimensional potential energy surface (PES) at the RCCSD(T)-F12 level under the rigid-rotor approximation. Close-coupling (CC) calculations were performed for low rotational levels, while the coupled-states (CS) method was used for higher levels. Convergence tests were conducted, and CC and CS scattering cross sections were compared. The CS rate coefficients cover rotational levels up to $j=40$ and temperatures of 3--60\,K. Although CC and CS rate coefficients differ by up to a factor of $\sim 2$ for some transitions, particularly those affected by weak coupling and resonances, they yield nearly identical excitation and brightness temperatures for the high-$j$ transitions examined, including $19\rightarrow18$. These results support the reliability of the CS method for non-LTE modelling of the high-$j$ transitions considered, for which CC calculations are computationally impractical. We use reduced-mass scaling of the SiC$_6$--He rate coefficients to estimate those for SiC$_6$--H$_2$ and model the observed transitions. A $\chi^2$ analysis of the observed line fluxes, excluding heavily blended features, yields a best-fitting of $n(\mathrm{H_2})=4.33\times10^5$\,cm$^{-3}$ and $N(\mathrm{SiC_6})=8.11\times10^{11}$\,cm$^{-2}$. The latter agrees excellently with the previous observational estimate of $(7.9\pm1.5)\times10^{11}$\,cm$^{-2}$, supporting the consistency of our non-LTE modelling with observational constraints.
A faint dust ring beyond the compact disk of Sz 76
oai:arXiv.org:2610.09881v2
arXiv:2610.09881v2 Announce Type: replace
Abstract: We report new ALMA Band 7 high-resolution ($\sim$0.05$^{''}$) continuum and $^{12}\text{CO}\ (J=3\text{--}2)$ observations of the disk around the young low-mass star Sz 76. Visibility-plane modeling reveals a previously undetected, faint outer dust ring centered at $\sim$45 au surrounding the compact central disk. Including this low-surface-brightness feature increases the characteristic dust radius from $R_{\text{dust,68}} \approx 7$ au to $41$ au. Combined with an outer gas extent of $R_{\text{CO,68}} \approx 77$ au, the resulting gas-to-dust size ratio decreases from an extreme value of $\sim$11 to $\sim$2, placing Sz 76 in line with typical protoplanetary disks and contradicting a scenario of unimpeded inward radial drift. These results demonstrate that spatial filtering or limited surface-brightness sensitivity can severely underestimate disk sizes, highlighting the necessity of multi-scale observations to properly identify dust traps and accurately constrain disk evolution.
Structure and evolution of cosmic voids within DESI DR1
oai:arXiv.org:2610.18348v2
arXiv:2610.18348v2 Announce Type: replace
Abstract: Ongoing stage IV galaxy surveys such as DESI or Euclid offer an unprecedented amount of information to study the large-scale structure of the Universe by mapping tens of millions of galaxies. Within the context of cosmic voids, already released and future catalogues will help better understand the void structure and evolution as well as the cosmological model constraints derived from their statistics. With the recent DESI Data Release 1, our aim is to obtain a new void catalogue that allows us to study the properties of cosmic voids at different redshifts, hence providing insight not only on their structure but also on their evolution.
We select the DESI BGS and LRG samples for this study. Before the void-finding process, the data goes through a preprocessing pipeline tackling redshift space distortions and reconstructing the mildly non-linear velocity field. After that, the new AVISM void finder is applied to the preprocessed data, obtaining two void catalogues in different redshift bins (0.1 < z < 0.3 and 0.4 < z < 0.8), both within a 1400^3 cMpc^3 cube. We successfully recover N_BGS = 1114 and N_LRG = 216 clean voids with sizes up to ~ 65 cMpc. The void size function, 3D shape variables, density profiles, and linear bias are obtained and analysed. Both the size function and the density profiles obtained are in agreement with previous theoretical studies. The expected mass depletion in void evolution is shown by the evolution of density profiles. We obtain a roughly constant linear bias inside r < Re with mean values of 1.31 +/- 0.02 and 2.14 +/- 0.04 for BGS and LRG, respectively. Both void catalogues are publicly released.
Optical constraints on intranight variability in TeV HBLs from a two-year Asiago Schmidt monitoring campaign
oai:arXiv.org:2610.20065v2
arXiv:2610.20065v2 Announce Type: replace
Abstract: We present a two-year homogeneous optical monitoring campaign of three TeV-emitting high-synchrotron-peaked BL Lac objects, 1ES 1959+650, H 1426+428, and PG 1553+113, carried out with the Asiago Schmidt 67/92 cm telescope. TeV blazars are key targets in high-energy astrophysics, relevant for relativistic jet physics, particle acceleration, multiwavelength variability studies, and, in selected cases, multimessenger programmes. In this context, dense optical monitoring provides a low-energy synchrotron tracer of jet activity over temporal baselines that are difficult to cover continuously with high-energy facilities. The dataset consists of more than 17,000 R- and I-band photometric measurements obtained over 957 source-filter observing sessions between 2024 February 1 and 2026 January 31, for a total monitoring time of about 420 h. We applied a uniform differential-photometry and classical F-test framework based on two independent comparison-star pairs. Robust IDV is detected only in H 1426+428, with 15 out of 43 analysed sessions classified as variable in the R band and 5 out of 40 in the I band. No significant IDV is detected in 1ES 1959+650 or PG 1553+113 within the adopted statistical framework. For H 1426+428, the corresponding duration-weighted duty cycles, computed using only variable sessions with a measured local characteristic timescale, are 28.0% in the R band and 7.0% in the I band. The inferred local extrapolated doubling/halving timescales are of the order of 5-15 min in the observer frame. These results indicate that rapid optical synchrotron variability in TeV HBLs is intermittent and state-dependent. They also show that robotic small- and medium-aperture optical facilities can provide valuable long-term support for high-energy monitoring programmes by supplying dense, uniform, and repeatable optical coverage.
Consistently Reproducing Early Galaxy Observations and 21-cm Signal during the Epoch of Reionization
oai:arXiv.org:2610.24064v2
arXiv:2610.24064v2 Announce Type: replace
Abstract: In the standard cosmological model, galaxies forming within dark matter haloes are the primary drivers of the Epoch of Reionization (EoR). The complex physics of galaxy formation introduces scatter in the relation between halo mass and star-formation activity, with implications for the production and escape of ionising photons. In this work, we implement a stochastic source model in our cosmological radiative transfer simulations. We calibrate our model against the high-redshift ultraviolet luminosity function observed by JWST and observational constraints on Lyman-continuum leakage to constrain the ionising escape fraction. We find that stochastic star formation redistributes ionising emissivity toward lower-mass haloes, producing a more fragmented ionisation topology and delaying the growth of the largest ionised regions relative to the deterministic, observationally calibrated model. However, when the two models match total ionising photon production, their global neutral fractions and 21-cm power spectrum remain nearly indistinguishable, despite small differences in the morphology of individual ionised regions. Additionally, we establish a unified framework that lets us apply source models derived from computationally expensive hydrodynamical simulations (SPICE, FIRE-2 and Thesan) and investigate their impact on large-volume reionisation simulations, extending their prescriptions to scales of >100 Mpc. We find that these prescriptions produce distinct ionisation topologies and 21-cm statistical signatures, even when their global reionisation histories are similar. Our results highlight the importance of combining physically motivated galaxy formation prescriptions from hydrodynamical simulations with large-volume radiative transfer simulations to interpret forthcoming 21-cm observations and constrain the sources of cosmic reionisation.
Beyond the Disc--Spheroid Dichotomy at $0.6 < z < 8.0$: Evidence for a Morphological Indeterminacy Regime from JWST
oai:arXiv.org:2610.24078v2
arXiv:2610.24078v2 Announce Type: replace
Abstract: Galaxy morphology is a fundamental tracer of galaxy assembly, but its characterization becomes increasingly uncertain in the early Universe. We investigate morphological classifications across different methodologies over the redshift range $0.6 < z < 8.0$, using JWST observations of 15,567 galaxies with 6,099 synthetic galaxies from IllustrisTNG. To quantify morphology without relying on visual labels or prescribed light profiles, we use an unsupervised-supervised hybrid machine-learning approach based on the non-parametric MEGG system, which combines the second moment of the light distribution, Entropy, Gini, and Gradient Pattern Analysis. We find that parametric and non-parametric classifications remain broadly consistent at $z < 2.5$, but diverge systematically at higher redshift. We identify this epoch as the Morphological Indeterminacy Regime (MIR), in which classical disc and spheroidal categories no longer provide a unique description of the galaxy population. The fraction of mixed-interpretation galaxies increases with redshift, reaching $F_{\rm MIR} \sim 0.55$ at $z > 2.5$. Cross-domain and projection-invariance tests show that the classifications are stable and insensitive to data set and viewing angle. In the interval $2.5 < z < 8.0$, MIR galaxies occupy an intermediate stellar-mass regime (median; $\log(M_\star/M_\odot) \sim 9.7$), between low-mass discs (median; $\log(M_\star/M_\odot) \sim 9.5$) and massive spheroids (median; $\log(M_\star/M_\odot) \sim 10.1$), while irregular systems have stellar-mass distributions similar to those of MIR galaxies. We interpret these results as evidence that morphological indeterminacy is a characteristic feature of galaxy assembly in the early Universe, reflecting a population of dynamically young systems that have not yet settled into the classical Hubble sequence.
The DARWIN Project: Cosmological Simulations of Galaxy Formation down to Dwarf Scales in Diverse Environments
oai:arXiv.org:2610.25786v2
arXiv:2610.25786v2 Announce Type: replace
Abstract: We present the DARWIN project, a suite of cosmological simulations with RAMSES to study dwarf-galaxy formation and evolution across diverse environments. DARWIN-1 is a full-box simulation in a cubic volume of $(65.5\,{\rm cMpc})^3$ with a minimum cell size of 500pc, while DARWIN-2 comprises zoom-in simulations of 24 field regions and a Milky Way analog connected to a filament, reaching 62.5pc resolution. The simulations employ baryonic physics tailored to their aims. For DARWIN-2, we include radiative transfer in three photon energy bands to trace hydrogen and helium ionization. DARWIN-1 reproduces the cosmic star formation rate density and galaxy stellar mass functions reasonably well at $z>1$. Both simulations broadly reproduce observed stellar mass-halo mass relations, chemical abundances, star-formation surface densities, and black hole-galaxy scaling relations. At lower redshift, DARWIN-1 shows an excess in the number density of massive galaxies, although the inferred stellar content of massive halos depends on the definition of galaxy stellar mass. Other systematic differences arise from spatial resolution, feedback efficiency, and environmental selection. To facilitate comparison with observations, we generate mock images using the radiative transfer codes SKIRT and Powderday. We also explore a machine-learning approach to mitigate resolution-driven biases in the dwarf galaxy properties of DARWIN-1 using the higher-resolution DARWIN-2 simulations. The DARWIN simulation suite therefore provides a comprehensive numerical framework for studying the formation, evolution, and internal physics of dwarf galaxies in a cosmological context.
The Central Molecular Zone of our Galaxy, a chemical complexity factory: first interstellar detections of molecules towards G+0.693
oai:arXiv.org:2610.26308v2
arXiv:2610.26308v2 Announce Type: replace
Abstract: The Central Molecular Zone (CMZ) of our Galaxy is a prime laboratory for interstellar chemistry because of its extreme conditions, such as shocks and high cosmic-ray ionisation rates. It is an efficient factory of molecules, with many species first discovered in CMZ sources. This document reviews key results from the study of the CMZ G+0.693-0.027 molecular cloud. Using ultra-sensitive broadband spectral surveys together with laboratory spectroscopy, more than 140 molecules have been identified, including 28 detected for the first time in the interstellar medium (ISM). These discoveries include prebiotic precursors of biomolecules, increasingly complex members of known molecular families, molecules containing periodic third-row elements (S, P, Na, Mg), high-energy stereo- and structural isomers, and new cations. Overall, these findings extend the known limits of interstellar chemical complexity and highlight the CMZ's potential to provide the molecular ingredients for prebiotic chemistry.
New Bounds on Heavy QCD Axions from Big Bang Nucleosynthesis
oai:arXiv.org:2510.23695v3
arXiv:2510.23695v3 Announce Type: replace-cross
Abstract: We study Big Bang Nucleosynthesis (BBN) constraints on heavy QCD axions. BBN offers a powerful probe of new physics that modifies the neutron-to-proton ratio during the process, thanks to the precisely measured primordial Helium-4 abundance. A heavy QCD axion provides an attractive target for this probe, because not only is it a well-motivated hypothetical particle by the strong CP problem, but also it dominantly decays to hadrons if kinematically allowed. A range of its lifetime is thus excluded where the hadronic decays would significantly alter the neutron-to-proton ratio. We compute axion-induced modification of the neutron-to-proton ratio, and obtain robust upper bounds on the axion lifetimes, as low as 0.017 s for the axion mass higher than 300 MeV. Remarkably, this is stronger than projected future CMB bounds via $N_{\rm eff}$. Our bounds are largely insensitive to uncertainties in hadronic cross sections and the axion's branching fractions into various hadrons, as well as to the precise value of the initial axion abundance. We also incorporate, for the first time, several key improvements, such as scattering processes by energetic $K_L$ and secondary hadrons, that can also be important for studying general hadronic injections during BBN, not limited to those from axion decays.
Inflationary Particle Production and the Swampland
oai:arXiv.org:2512.07930v3
arXiv:2512.07930v3 Announce Type: replace-cross
Abstract: We investigate the impact of particle production during inflation in scenarios where an infinite tower of states features a mass scale that decreases exponentially along the inflationary trajectory. Such couplings naturally arise in string effective field theories and are in fact motivated by the Swampland Distance Conjecture (SDC). We show that the corrections to inflationary observables sourced by the tower scale as $(H/\Lambda_{\text{sp}})^{2+p}$, with $H$ being the Hubble scale, $\Lambda_{\text{sp}}$ being the species scale, that is the quantum gravity cut-off, and $p\geq 1$ characterizes the density of states in the tower. As a result, in gravitationally weakly coupled cosmological effective theories, the tower-induced contributions are suppressed relative to the standard single-field predictions, leaving the inflationary phenomenology essentially unchanged. We demonstrate this explicitly across a set of well-motivated inflationary potentials, and we compare the resulting predictions with the most recent observational constraints, including those from the Atacama Cosmology Telescope.
Probing freeze-in dark matter using Bose-Einstein condensate in neutron star
oai:arXiv.org:2605.22245v2
arXiv:2605.22245v2 Announce Type: replace-cross
Abstract: Neutron star (NS) is one of the most promising astrophysical targets to probe non-gravitational interaction of dark matter (DM) with visible matter. Their compactness makes them an ideal object which can capture particle DM efficiently over its lifetime using the DM-nucleon scattering cross-section. If DM particles are bosonic, then the captured DM population may form a Bose-Einstein condensate at the center of the NS, increasing the DM density significantly. In this work, we study the phenomenology of such scenario with enhanced DM annihilation rate due to the increased density in a condensate. The enhanced DM annihilation makes the NS surface `hotter' than in the standard cooling scenario. We show that the annihilation rate is enhanced by a factor of $\mathcal{O}(10^{15}-10^{20})$ if DM forms a condensate without large self-interaction, and DM with freeze-in value annihilation cross-section can heat up the NS to higher temperatures, bringing it within the reach of James Webb Space Telescope. It also allows us to probe DM-nucleon scattering cross section within the neutrino fog regime which will complement the terrestrial direct detection searches. Moreover, the enhanced annihilation from the condensate changes the lower limits on s-wave DM annihilation cross-section for capture-annihilation equilibrium and the formation of a black hole inside the NS. Finally, we show an example of a scalar DM model where such small annihilation and DM-nucleon scattering cross sections can generically arise.
Energy extraction from NED-deformed rotating black holes via the Comisso-Asenjo reconnection process
oai:arXiv.org:2605.26369v2
arXiv:2605.26369v2 Announce Type: replace-cross
Abstract: We study rotating black holes in general relativity coupled to nonlinear electrodynamics (NED), focusing on an axisymmetric solution with deformation parameter $g$. On the spherical seed, weak-field lensing via the Gauss-Bonnet method and the shadow radius yield a spin-insensitive bound by enforcing a conservative $\sim 10\%$ tolerance on the Sgr A* ring size, namely $g/M \lesssim 1.26$. In the eikonal regime we derive analytic quasinormal-mode shifts, even in $g$, and obtain an independent ceiling consistent with the shadow constraint. For the rotating geometry, we provide closed-form ZAMO scalars, chart horizons and the ergoregion, and analyze equatorial geodesics (photon orbits and ISCO). We then formulate the Comisso-Asenjo reconnection channel in the ZAMO frame, identify the negative-energy window, and integrate the extracted power over the allowed radii; from the tolerated fractional departure from the Kerr power we define a spin-dependent extraction bound $g_\delta(a\mid\sigma_0,\xi)$. Taken together, the QNM/shadow ceiling and the extraction bound appreciably narrow the admissible region for $g/M$ in the $(a/M,g/M)$ plane, so even within our deliberately simplified, single-layer equatorial setup, the two complementary probes already provide informative constraints on NED deformations, testable with present data and upcoming horizon-scale and ringdown campaigns.
Detecting Tidal Resonances in Binary Neutron Stars
oai:arXiv.org:2606.06376v2
arXiv:2606.06376v2 Announce Type: replace-cross
Abstract: As a binary neutron star inspirals due to the emission of gravitational waves, the rising tidal frequency resonantly excites vibrational modes. These oscillations are seismological probes of the rich stellar interior, yet it remains to be established whether gravitational-wave interferometers can measure them. Here, we present the first fully Bayesian study of the capability of the Einstein Telescope to detect tidal resonances. We simulate one year of observations and analyse the 200 loudest signals. We find that the Einstein Telescope can identify resonant modes and is sensitive to gravitational-wave phase shifts as small as $\Delta \Phi \approx 0.03$ for favourable events. We further show that neglecting resonances can bias the inferred tidal deformabilities. These results establish tidal resonances as a measurable route for asteroseismology with future detectors.
A dark-matter-sensitive orientation clock near rotating black holes
oai:arXiv.org:2606.06557v3
arXiv:2606.06557v3 Announce Type: replace-cross
Abstract: A coherent triaxial structure has an intrinsic orientation clock because rotation about its intermediate principal axis is unstable. I ask how this clock changes when the structure moves on a circular equatorial orbit near a rotating black hole in an effective dark matter (DM) environment. The calculation uses a comoving orthonormal tetrad, the electric part of the Riemann tensor, Euler--quaternion dynamics, and the conversion from proper time to stationary coordinate time. The clock shift separates exactly into a local proper-time response and a geometric time-conversion response. I compute a direct grid of 600 DM models for Einasto, regularized cored Navarro--Frenk--White, and Hernquist profiles. Across the scanned domain, the total fractional frequency shift lies between about $-3\times10^{-7}$ and $-3\times10^{-4}$. For the fiducial body model, most of the shift comes from the time conversion. Synthetic timing tests show that a shift near $10^{-4}$ is recoverable only when the clock is highly coherent and its Kerr baseline is known independently. This gives a theoretical timing diagnostic and a synthetic observational forecast. I do not fit real data.
Resonant lunar tides of Earth's core and basal magma ocean
oai:arXiv.org:2606.16340v2
arXiv:2606.16340v2 Announce Type: replace-cross
Abstract: Earth's magnetic field is generated by fluid motion in the liquid-metal core and has been active for billions of years. However, prior to the onset of inner-core growth, the sources of mechanical power that drove the geodynamo remain uncertain. During this period, the core may have been overlain by a basal magma ocean (BMO), creating two immiscible fluid layers separated by a density interface, beneath the solid mantle. We develop a theory for lunar tides in this core--BMO system, in which the tidal potential acts on the density contrast between the two layers. The resulting dynamics differ fundamentally from previous models driven by imposed deformation of the core's bounding envelope. In the inviscid limit, the response transitions from an equilibrium tide to an interfacial gravity wave as forcing frequency increases. The two regimes are separated by a resonance that occurs when the forcing frequency matches the natural frequency of the degree-2, order-2 interfacial mode. The inviscid core flow is formally identical to the canonical elliptical vortex, linking the problem to the theory of elliptical instability. Finite viscosity regularises the resonance, introduces phase lag and generates oscillatory boundary layers. Combined with parameterised models of lunar recession and BMO crystallisation, the theory predicts enhanced core-boundary ellipticity, core-flow speed, magnetic Reynolds number and instability metrics, particularly near resonance. These results identify a previously unexplored mechanism for tidally driven flow in differentiated planetary bodies and show that a BMO can enhance tidal coupling to the core, potentially contributing to dynamo action.
Gravitational-wave responses and noise correlations in geometric time-delay interferometry
oai:arXiv.org:2606.28173v2
arXiv:2606.28173v2 Announce Type: replace-cross
Abstract: Space-based gravitational-wave observations will probe the formation and growth of massive black holes and the evolution of compact binaries in the millihertz band. We develop a geometric-link representation of time-delay interferometry using six round-trip and six non-round-trip two-link blocks. Starting from the one-way Doppler response, we derive the elementary block responses and a dictionary that expresses 45 known second-generation combinations in terms of these blocks and residual one-way measurements. Representative Michelson and Relay observables admit compact forms, from which analytical sky-averaged tensor responses are obtained for a fixed, equal-arm constellation. We also derive the noise covariance of the blocks, retaining cross-spectra between blocks that share measurement streams. The resulting noise spectra agree with direct one-way assembly. We apply this formulation to a narrowband compact-binary signal in an equal-arm LISA instrument model. For the Relay example, neglecting the shared-link cross-spectrum changes the low-frequency noise scaling and overestimates the conditional Fisher uncertainties in amplitude or distance by factors of 10.71 at 1~mHz and 3.37 at 3~mHz.
Exact 1D Nonlinear Solutions for Proton-Driven Plasma Wakefields
oai:arXiv.org:2607.06458v3
arXiv:2607.06458v3 Announce Type: replace-cross
Abstract: We extend the classical 1D nonlinear relativistic plasma wakefield paradigm to evaluate the pristine fluid profiles driven by intense relativistic proton trains under the comoving coordinate approximation. Using an intense two-bunch scheme, our adaptive integration model is validated down to machine precision against exact closed-form cold-fluid invariants, confirming that interface boundary corners represent physical step discontinuities in the second potential derivative. Scaling to an $N=100$ micro-bunch train, we resolve macroscopic phase-accumulation and subsequent nonlinear detuning by deriving a closed-form, parameter-free piecewise analytical envelope expression. This solution successfully maps the initial unattenuated quadratic coherent growth alongside the flat asymptotic horizontal saturation horizon governed by nonlinear energy conservation. Furthermore, the framework analytically evaluates the transient envelope overshoot, tracking its decay onto the saturation boundary via an exact phase-slip integration length scale. Providing immediate, ready macro-envelope bounds, this computationally inexpensive tool eliminates empirical adjustments and yields a robust mathematical baseline designed to guide and accelerate complex, multi-dimensional particle-in-cell (PIC) simulations through rapid parameter space exploration before deploying heavy computational assets. It establishes an efficient foundation to optimize asymmetric, custom-shaped driver profiles required to break the fundamental symmetric transformer ratio limit of 2 to make current particle accelerators more efficient.
Entropic Chaplygin-Gas cosmology and late-universe tension diagnostics
oai:arXiv.org:2607.22802v2
arXiv:2607.22802v2 Announce Type: replace-cross
Abstract: We investigate an entropic generalized Chaplygin-gas (gCg) cosmology in which a single dark fluid accounts for both dark matter and dark energy. We constrain the model using DESI baryon acoustic oscillations, redshift-space-distortion measurements, compressed Planck 2018 distance and lensing information, and Type Ia supernova likelihoods with PantheonPlus+SH0ES (PPS) and DES--Dovekie. For Planck+DESI+RSD+PPS, the entropic gCg model yields a Hubble constant of $H_{0}=68.64\pm0.52\,{\rm km\,s^{-1}\,Mpc^{-1}}$ and the amplitude of matter fluctuations parameter of $S_{8}=0.797\pm0.021$, compared with $H_{0}=67.31\pm0.28\,{\rm km\,s^{-1}\,Mpc^{-1}}$ and $S_{8}=0.809\pm0.021$ in $\Lambda$CDM. The simultaneous upward shift in $H_{0}$ and downward shift in $S_{8}$ occurs in the directions required for a partial alleviation of both tensions. We also show that when PPS is replaced by DES--Dovekie, the parameter shifts are smaller, the best-fit improvement is negligible, and the Bayesian evidence favours $\Lambda$CDM. Thus, the entropic gCg model provides a potentially interesting simultaneous shift of the Hubble constant and the matter-clustering amplitude, but its statistical support depends on the adopted supernova likelihood.
Penrose process in radiating Kerr-Vaidya black holes: efficiency, energy transport, and the viability boundary
oai:arXiv.org:2609.19175v2
arXiv:2609.19175v2 Announce Type: replace-cross
Abstract: We study single-particle Penrose extraction from a rotating black hole whose mass decreases in time, described by the ingoing Kerr-Vaidya metric, imposing exact four-momentum conservation at the decay. We prove that the local efficiency is bounded by the instantaneous Kerr value, through a covariant bound that reduces to the Bardeen-Press-Teukolsky limit at the horizon, and we derive an exact transport law, dE/dtau = -dm/dt (r/Sigma) p_r^2 for fixed spin parameter, by which every timelike geodesic gains energy while the hole loses mass. The net energy gain therefore splits into a local Penrose term and two transport terms, which we compute numerically as functions of spin and mass-loss rate for two mass-loss laws. With fixed spin parameter, extraction is possible only below a critical rate alpha_crit(a) = (M-a)/Delta v, set by the approach of the hole to extremality; with fixed a/m no such boundary exists. The transport law extends to the Kerr-Newman-Vaidya family, and the regime where the dynamics matter is shown to lie far outside astrophysical mass-loss rates.
oai:arXiv.org:2604.09763v2
arXiv:2604.09763v2 Announce Type: replace
Abstract: Dust attenuation affects galaxy SEDs and introduces uncertainties in star formation rates, stellar masses, and metallicities. Direct spectroscopic constraints at early cosmic epochs were limited before JWST. We derive the average dust attenuation law of massive (log10(Mstar/Msun) > 9) star-forming galaxies at 2 < z < 7 using JWST spectroscopy. We combine JWST/NIRSpec spectroscopy from JADES with deep multi-wavelength photometry from the ASTRODEEP-JWST catalogues. For approximately 120 mass-selected galaxies with reliable Balmer decrement (Halpha/Hbeta) measurements, we stack SEDs in bins of Balmer optical depth and derive the selective attenuation curve following Calzetti et al. (1994). JWST/MIRI photometry extends the wavelength coverage into the near-infrared (NIR). The resulting curve spans 0.16-1.14 microns in the rest frame and is well described by a smooth function. We derive R_V approximately 3.98, with Delta R_V = 0.16 indicating its sensitivity to the adopted long-wavelength extrapolation. The average attenuation law agrees remarkably well with the local starburst relation in slope and normalisation, but appears systematically flatter in the ultraviolet than several intermediate-redshift determinations. We find no significant evidence for a 2175 Angstrom UV bump. Our results provide a direct empirical constraint on the average dust attenuation law of massive star-forming galaxies at 2 < z < 7. This analysis exploits the combined spectroscopic and photometric reach of JWST. Despite variations among individual galaxies, the ensemble-average behaviour agrees with the local starburst relation, suggesting that the main mechanisms governing dust attenuation were already operating at early cosmic times. The relatively shallow UV slope may reflect differences in dust-star geometry, radiative-transfer effects, or effective grain-size distributions among galaxy populations.
Silicate cosmic dust grain collisions in the interstellar medium: A molecular dynamics study
oai:arXiv.org:2605.22225v3
arXiv:2605.22225v3 Announce Type: replace
Abstract: (abridged) We aim to predict the most important parameters for grain-grain collision outcomes in models of interstellar grain population evolution: the threshold velocities for shattering and vaporization, and the resulting distributions of grain sizes. We use molecular dynamics simulations which evolve the dynamics of each atom in a dust grain to explore the outcomes of collisions between silicate grains of radii $a \in [5,50]~\AA$ at velocities 0.1-20 km/s. We run simulations of grains with two materials: amorphous SiO$_2$ and an amorphous silicate of composition suggested by Draine \& Hensley (2021). With these simulations, we quantify the collision velocity dependence of shattered and vaporized mass fractions, and the resulting size distributions of shattering products. We find grain shattering thresholds are $\sim$6 km/s for both materials, a factor of $\sim$2 higher than the canonical value for silicates of 2.7 km/s from Tielens et al. (1994) and Jones et al. (1996). The size distributions of shattered products are generally not consistent with the power law distributions predicted by Jones et al. (1996), but the physical reasons for these discrepancies remain uncertain. The Tielens et al. (1994) expressions fail to predict the fraction of shattered or vaporized material observed in our numerical simulations, over-predicting the shattered fraction and under-predicting the vaporized fraction, except for the most massive grains we simulate, where the predicted vaporized fraction agrees reasonably well. The model of Hirashita \& Kobayashi (2013) similarly fails to match the simulations. We provide updated shattering velocity thresholds for nanoscale silicate grains, which suggest that the smallest astrophysical dust grains composed of silicate materials may be more robust to shattering in the interstellar medium than previously assumed.
The Life and Death of Stars That Capture Primordial Black Holes
oai:arXiv.org:2606.02700v2
arXiv:2606.02700v2 Announce Type: replace
Abstract: Primordial black holes (PBHs) in the asteroid mass window ($10^{17}-10^{23}\,{\rm g}$) remain viable dark matter candidates and can be captured by stars. We develop the first global framework for the evolution of stars that capture PBHs, combining analytic calculations, stellar evolution models, 3D general-relativistic magnetohydrodynamic simulations, and Monte Carlo population synthesis. We find that the fate of these systems bifurcates: PBHs that form an accretion disk before consuming the host drive explosive disruption, whereas PBHs captured too late or growing too slowly consume the star quietly. Capture is dominated by three-body interactions with planetary or stellar companions. For a solar-type host with a Jupiter analog, inspiral within a main-sequence lifetime requires $M_{\rm BH}^{\rm crit}\gtrsim 10^{22}\,{\rm g}$, while lighter PBHs generally require tighter companions. Once deposited at the center, the PBH grows through inefficient quasi-spherical Bondi accretion; if it reaches the angular-momentum threshold before consuming the host, the inflow circularizes into a disk. Our Monte Carlo calculations yield sizable quiet-consumption and explosive-disruption populations, with final PBH masses $M_{\rm BH}\sim0.01-1\,M_\odot$ and disk-forming PBH spins $a_\ast\approx0.8$. Disk formation is the point of no return: disk winds and relativistic jets of $\sim10^{45}-10^{50}\,{\rm erg\,s^{-1}}$ disrupt the star within minutes. The resulting transients may include a $\sim$day-long UV/blue signal, radio afterglow, and, if the jet escapes, an X-ray-flash/low-luminosity gamma-ray-burst (XRF/llGRB) signal. For an $O(1)$ PBH dark matter fraction and optimistic capture assumptions, the event rate can reach that of llGRBs. The low-mass, high-spin remnants offer a complementary PBH probe and possible source for subsolar BH mergers.
The Effect of Adiabatic Index on Radius Evolution and the Mass Loss
oai:arXiv.org:2606.17143v2
arXiv:2606.17143v2 Announce Type: replace
Abstract: Models that track the size evolution of exoplanets often assume a prescribed initial thermal state or a single adiabatic index to describe the planetary interior structure, the latter of which is taken to be $\gamma \sim1.4$ which is likely appropriate for evolved planets ($\gtrsim$1 Gyr). Extrapolating this high $\gamma$ to earlier ages (down to $\sim$million years old) is problematic since, according to thermodynamics, the adiabatic index of young planets is $\sim$1.2, which is low enough to drastically change how interior mass is distributed. We quantify the effect of varying the adiabatic index from 1.2 to 1.4 on the expected radius of the exoplanet over time. We find that envelopes of larger adiabatic indices start puffier with all else equal and undergo faster radius contraction with accelerated mass loss. Assumption of high $\gamma$ can therefore overestimate the effect of mass loss in shaping the exoplanetary population, especially when young planets are considered. We highlight the need for a more careful consideration of the initial thermal condition of planets in evolutionary models to properly interpret the radii measurements of exoplanets.
A New Scaling of Neutron Star Tidal Deformability for Directly Probing the Central Equation of State Parameter
oai:arXiv.org:2606.21402v2
arXiv:2606.21402v2 Announce Type: replace
Abstract: The dimensionless tidal deformability, $\Lambda$, of neutron stars (NSs), inferred from gravitational-wave (GW) observations, has thus far been used primarily to constrain the pressure of dense matter near twice nuclear saturation density, leaving the core equation of state (EOS) largely inaccessible to inspiral-phase GW observations. We show that a fundamental central EOS parameter can be probed directly through $\Lambda$ using a perturbative analysis of the dimensionless stellar-structure and tidal-response equations formulated in terms of scaled intrinsic variables, without invoking any specific EOS model. We uncover a remarkable EOS-insensitive scaling relation between $\Lambda$ and the central EOS parameter $\rm{X}\equiv P_{\rm c}/\varepsilon_{\rm c}$, where $P_{\rm c}$ and $\varepsilon_{\rm c}$ denote the central pressure and energy density, respectively. The relation is validated against a broad ensemble of physically viable EOSs. Applying it to tidal deformabilities inferred from events such as GW170817 enables a direct determination of $\rm{X}$. By adopting the previously established effective upper bound $\rm{X}\lesssim0.381$, we infer a tight lower bound, $\Lambda_{\rm{TOV}}\gtrsim9.1^{+0.5}_{-0.5}$, for maximum-mass NSs along stable mass-radius sequences, quantitatively demonstrating that even the most compact stable NSs remain distinctly separated from black holes, for which $\Lambda_{\rm{BH}}=0$. These findings reveal a previously unrecognized connection between inspiral-phase tidal deformability and the central EOS parameter $\rm{X}$, establishing a direct link between GW observables and the microphysics of ultradense matter in the strong-gravity regime. The resulting scaling establishes inspiral-phase tidal deformability as a direct and largely model-insensitive probe of the central EOS parameter $\rm{X}$, and hence of key properties of NS cores.
Search for Diffuse Supernova Neutrino Background in the Full KamLAND Dataset with Neural-Network-Based Event Classification
oai:arXiv.org:2606.29381v2
arXiv:2606.29381v2 Announce Type: replace
Abstract: We report a search for the diffuse supernova neutrino background (DSNB) with the KamLAND detector, targeting electron antineutrinos via inverse beta decay in the neutrino energy range of 8.3 to 30.8 MeV. Using liquid-scintillator exposures of 9.02 kton-year for 8.3 to 9.3 MeV and 9.42 kton-year for 9.3 to 30.8 MeV, we observe seven candidate events after applying a new deep-neural-network-based event classification technique. This result is consistent with the background-only expectation of 16.2 plus or minus 9.4 events, including systematic uncertainties associated with the neural-network selection. A spectral analysis of the energy and radial distributions finds no significant excess attributable to the DSNB. We therefore set 90 percent confidence-level upper limits on the DSNB flux of 38 to 43 per square centimeter per second, depending on the assumed DSNB model. We also derive model-independent 90 percent confidence-level upper limits on the electron-antineutrino flux, obtaining some of the most stringent constraints below 13.3 MeV. Beyond the DSNB search itself, this work demonstrates neural-network-based event classification as a promising approach for suppressing neutron-associated backgrounds in liquid-scintillator neutrino detectors.
Unveiling the Milky Way with a Gaia DR3 census of OB-type stars within 2 kpc. I. Tracing local Galactic structure, massive star-forming regions and core-collapse supernova progenitors
oai:arXiv.org:2607.07068v2
arXiv:2607.07068v2 Announce Type: replace
Abstract: O- and B-type stars are young and hot, thereby serving as vital tracers of the star formation and spiral arm structure of the Milky Way. At the dusk of the \textit{Gaia} DR3 era, a high-confidence and accurate catalogue appears timely. We have characterized 105,971 OB-type stars (T$_{\rm eff} >$ 10,000 K) within 2 kpc from the Sun, using an astro-photometric Bayesian inference tool. Our map unveils a complex view of the young stellar populations across the thin disk, with prominent large-scale features such as the Cepheus Spur, the Giant Oval Cavity, and a segment the Sagittarius-Carina arm all visible. Their inhomogeneous spatial distribution implies that massive star formation has taken place clustered across a few highly concentrated regions. We find a correlation between the overdensities of OB-type stars and young open clusters, and attribute departure from it to extinction effects. We identify 4200 stars as model-dependent core-collapse supernova (ccSN) or direct-collapse black hole (BH) progenitor candidates, including 6 ccSN progenitor candidates within 125 pc -- the maximum radius within which a ccSN explosion could be harmful to Earth's biosphere -- but our analysis suggests none of them exploding within the next 5 Myrs. Finally, we anticipate more evolved OB-type stars to collapse to form BHs within the next 1 Myr than ccSN progenitors will explode, despite the former's much scarcer number. A possible interpretation is a recent star formation burst in the local Milky Way, consistent with overdensities of ccSN progenitor candidates in Cyg OB2 and the Sagittarius-Carina arm.
Infrasound Newtonian Noise Estimation at the Einstein Telescope Candidate Site Sos Enattos
oai:arXiv.org:2607.21154v2
arXiv:2607.21154v2 Announce Type: replace
Abstract: We investigate the seasonal variability of atmospheric infrasound and its contribution to Newtonian noise (NN) at the Sos Enattos site, a leading candidate location for the Einstein Telescope (ET). Infrasound data recorded at three stations, SOE0 (surface), SOE1 ($-84$ m) and SOE3 ($-160$ m), are analysed over multiple seasons to characterise both the temporal variability and the depth dependence of the acoustic field. At the surface, the median amplitude spectral density (ASD) at 1 Hz shows no consistent seasonal modulation, and its variability is driven primarily by wind; the apparent dependence on sea state largely reflects its correlation with wind speed. Underground, where local wind noise is suppressed, a winter maximum of 4--7 dB above the quietest season is observed in both years of data. Using the measured sound-pressure spectra, we estimate the corresponding NN contribution for an underground cavern with the dimensions foreseen for ET, assuming an isotropic field of acoustic plane waves. The surface pressure field, acting through the atmosphere above a test mass at 200 m depth, exceeds the ET-D design sensitivity in median conditions between 1.6 and 7.1 Hz, by up to a factor of 4.7. Inside the cavern, the median NN remains a factor of 1.4 below ET-D with the pressure measured at SOE1, but exceeds it between 2.5 and 8.7 Hz, by up to a factor of 3.5, with the pressure measured at SOE3. The attenuation is not monotonic with depth: the deeper station SOE3 shows higher levels than SOE1 above a few hertz, pointing to local, anthropogenic sources in the underground infrastructure. These results show that underground placement is essential but not sufficient on its own, and that the acoustic emission of the underground infrastructure and future laboratory equipment will require careful evaluation and mitigation to keep atmospheric NN below the ET design sensitivity.
Observational Evidence for Anisotropic Metal Excess around Galaxies
oai:arXiv.org:2607.24929v2
arXiv:2607.24929v2 Announce Type: replace
Abstract: The exchange of matter and energy between galaxies and their surroundings drives the cosmic baryon cycle, yet mapping metal transport remains an observational challenge. While simulations predict that galactic winds escape anisotropically along minor axes, evidence for chemical enrichment in neighboring galaxies is limited. We analyze 1,433 galaxy pairs from the Dark Energy Spectroscopic Instrument survey and detect a gas-phase metallicity excess of 14.6% $\pm$ 3.7% to 24.2% $\pm$ 2.6% in neighbors aligned with the minor axis of massive primary at projected separations of 15--60 kpc, qualitatively consistent with IllustrisTNG simulation. Permutation tests indicate that this excess varies from a marginal detection ($>$92% confidence) at closer separations (15--30 kpc) to a significant signal ($>$98% confidence) at intermediate separations (30--60 kpc). Here, we show that these results provide empirical constraints on galactic feedback models and are consistent with a scenario where anisotropic outflows can enrich nearby galaxies, a mechanism that complements other environmental processes.
VAPOLA - A multiyear, multiband polarization survey of AGNs and Sgr A* at millimeter wavelengths with ALMA II. Spectropolarimetric properties and their evolution from 2017 to 2023
oai:arXiv.org:2608.14900v2
arXiv:2608.14900v2 Announce Type: replace
Abstract: We present a systematic analysis of the spectropolarimetric properties of a sample of 39 active galactic nuclei and Sagittarius A* observed with the Atacama Large Millimeter/submillimeter Array during five VLBI campaigns between 2017 and 2023. We characterize the compact cores in total intensity and polarization, focusing on the behavior or the linear polarization fraction (LP), electric vector position angle (EVPA), and Faraday rotation measure (RM) over time and spectral domains. We investigate both individual objects--such as M87, Sgr A*, 3C273, and 3C279--and ensemble properties of different source classes, including flat-spectrum radio quasars, BLLac objects, and other active galaxies. While total intensity and spectral index are generally stable on weekly timescales, polarization properties often exhibit strong variability, with significant day-to-day changes in LP, EVPA, and RM. Several sources display large EVPA rotations accompanied by variations in LP and RM, in some cases coinciding with flaring activity. We observe that the magnitude of RM increases with observing frequency for all sources for which we have reliable multi-band measurements, consistent with Faraday rotation arising in a magnetized sheath surrounding the relativistic jet, although an origin in the accretion flow--particularly in the case of Sgr A*--cannot be ruled out.
Periodic activity on the (sub)mm surface of the AGB star Mira A
oai:arXiv.org:2608.23724v2
arXiv:2608.23724v2 Announce Type: replace
Abstract: Evolved stars on the asymptotic giant branch (AGB) lose mass through stellar winds launched from their dynamic extended atmospheres. Recent high-angular-resolution observations of the archetypal AGB star Mira A (o Ceti) have revealed compact hotspots on its (sub)millimetre surface, whose origin remains unclear. We aim to characterise the (sub)millimetre photosphere of Mira A and investigate the properties and temporal behaviour of the compact hotspots detected on its surface. We analysed 19 epochs of high-angular-resolution observations obtained with ALMA in Bands 4, 6, 7, and 9. The complex visibilities were fitted with parametric models to derive the flux density, size, morphology, and brightness temperature of the stellar disc and an additional, more compact component across multiple pulsation phases and observing frequencies. Our results show that, near maximum light, the flux density of the stellar disc increases with frequency while its apparent size decreases, in agreement with the expected stratified radio photosphere. In Band 6, both the flux density and the size of the stellar disc appear to vary with pulsation phase, with the largest measured values occurring close to optical maximum light and the lowest observed around minimum light. Compact hotspots are detected in Bands 6 and 7 near maximum light, with their brightness temperatures exceeding those of the stellar disc by several thousand kelvin and reaching a temperature excess of ~8200 K in Band 6. They are highly energetic and are detected over timescales of weeks, while the available observations also suggest potential longer-term changes in both their position and intensity. Hotspot activity appears to occur repeatedly across multiple pulsation cycles and may be linked to strong shocks, magnetic activity, and/or episodic mass ejections, and could therefore play an important role in the mass-loss process in AGB stars.
Wide Field Localization Using Fixed Modulation Collimators
oai:arXiv.org:2609.10384v2
arXiv:2609.10384v2 Announce Type: replace
Abstract: Rapid, precise localization of gamma-ray bursts (GRBs) requires continuous wide-field imaging at hard X-ray and soft gamma-ray energies. We develop a fixed (non-rotating) modulation-collimator architecture to address this challenge. A bigrid collimator, a pair of aligned, periodic slit grids in front of a simple, non-imaging photon-counting detector, transmits an off-axis source through a triangular fringe pattern whose local angular position, the source phase, is recovered from simultaneous count rates in several grid pairs at fixed relative phase offsets, rather than by mechanically rotating a single grid pair as in the classical rotating modulation collimator. Four-phase demodulation of these count rates robustly cancels detector background and grid leakage without requiring rotation. Because the fringe pattern repeats across the sky, a phase measurement corresponds to many candidate source positions; we resolve this ambiguity over an arbitrarily wide field of view using a cascade of progressively coarser-pitched grid stages, generalizing the two-grid vernier-scale technique originally used to localize Sco X-1. This vernier cascade technique reduces the required photon counts and grid-calibration accuracy required to achieve fine angular resolution over a broad field of view. Generalized to two dimensions with orthogonal cascades, the resulting architecture achieves wide-field localization with no moving parts, no time modulation, and only simple photon-counting detectors, a promising candidate for simple, inexpensive, continuous-monitoring instruments supporting time-domain and multi-messenger astrophysics.
Thermal Evolution of Lava Planets Across System Ages: Predictions for Hell of a Survey
oai:arXiv.org:2609.20155v2
arXiv:2609.20155v2 Announce Type: replace
Abstract: Ultra-short-period (USP) rocky exoplanets can have dayside temperatures high enough to maintain permanent magma oceans, sitting at the intersection of interior geophysics and atmospheric chemistry. Coupled feedbacks between the molten surface and outgassed atmosphere can sustain or enhance a volatile envelope, while stellar interactions can erode it. Understanding which outcome prevails, and its observable imprint, requires a multi-target approach across planets at different stages of thermal evolution. We present predictions for the five targets of JWST Cycle 4 program 8864: TOI-1807 b, TOI-2260 b, TOI-431 b, TOI-6255 b, and TOI-2431 b. Using the PROTEUS coupled interior-atmosphere framework, we construct a simulation grid and classify outcomes into six categories, defined by the final interior melt state and by whether the planet retains a detectable atmosphere thick enough to redistribute heat to the nightside. For targets retaining a non-negligible volatile envelope, our models predict higher partial pressures for most species when the surface is molten, except for S$_{2}$, whose enhancement in the solid regime suggests it may serve as a tracer of interior melt state. Despite some targets showing outcomes across multiple scenarios, most tend toward a bare-rock end-member, with global melt fraction $\leq$ 20\% and atmospheric retention sensitive to escape efficiency. Our analysis reveals a minimum escape efficiency threshold below which volatile envelopes survive under energy-limited escape, constraining the conditions required for atmosphere survival on irradiated rocky planets. These predictions will guide interpretation of MIRI-LRS phase curve observations and identify which targets and features best discriminate between competing geophysical states.
Solar spicule population properties inferred from IRIS Mg II k off-limb spectra using likelihood-free Bayesian inference
oai:arXiv.org:2609.37811v2
arXiv:2609.37811v2 Announce Type: replace
Abstract: Solar spicules are ubiquitous and strongly overlap in off-limb observations, which complicates the interpretation of individual events. We aim to infer population-level properties of spicules from IRIS Mg II k off-limb spectra using a forward model that explicitly accounts for line-of-sight superposition and instrumental effects. We constructed a stochastic 3D spicule-population model and synthesized Mg II k spectra with radiative transfer, including instrumental degradation and photon noise. We then constrained the hyperparameters that describe the underlying spicule population with likelihood-free Bayesian inference. This approach combines approximate Bayesian computation and an adaptive population Monte Carlo sampler, comparing observed and simulated datasets through principal-component summary statistics at five heights above the limb. The inferred posterior constrains effective population-level properties of the spicules that dominate the observed Mg II k variability and height dependence. These constraints provide a route to connect off-limb spectroscopy with forthcoming spectropolarimetric diagnostics of the chromosphere--corona interface.
Galaxy Protoclusters as Drivers of Cosmic Reionization: II. Te-Based Metallicities of Lyman-{\alpha} Emitters
oai:arXiv.org:2610.01811v2
arXiv:2610.01811v2 Announce Type: replace
Abstract: Context. Protoclusters at reionization are thought to host a significant fraction of the cosmic star formation rate density, making them crucial for the timing and topology of reionization. The gas-phase metallicity, a key tracer of galaxy evolutionary stage, can now be measured directly in the early Universe with JWST. Aims. We measure the gas-phase metallicity of Lyman-$\alpha$ emitters (LAEs) in the LAGER-z7OD1 protocluster at $z\approx6.93$ via the direct ($T_e$-based) method using the [O III]$\lambda$4363 auroral line, and combine it with stellar masses to constrain the mass-metallicity relation (MZR) in dense environments. Methods. We use JWST/NIRSpec MSA G395H spectroscopy, derive metallicities with self-consistent analytic prescriptions, and propagate the dust-extinction uncertainty into an asymmetric error budget via Monte Carlo simulations based on the observed Balmer decrements. Stellar masses come from BEAGLE SED fitting of JWST/NIRCam and UltraVISTA photometry. Results. We obtain five secure [O III]$\lambda$4363 detections and one upper limit, with 12+log(O/H) $\approx$ 7.2-8.0. The weighted spectral stack yields 12+log(O/H) = $7.76^{+0.06}_{-0.05}$, ~0.26 dex higher than stacks of LAEs at similar redshifts from JWST surveys. However, the MZR fit gives a slope $\gamma$ = 0.232 +/- 0.214 and normalization $Z_0$ = 7.888 +/- 0.199, placing our sources below the established MZR of normal galaxies and recent protocluster baselines at $z\sim5-7$. Conclusions. The higher metallicity relative to general LAEs suggests accelerated assembly driven by the overdensity, while the offset below the MZR points to a complex interplay between mass build-up and gas accretion. Since we target actively star-forming LAEs, we are likely capturing a phase of pristine gas inflow from the intergalactic medium that fuels star formation but dilutes the interstellar media of these protocluster members.
Collisional excitation of SiC$_6$ by helium: cross sections, rate coefficients, and astrophysical implications
oai:arXiv.org:2610.06525v3
arXiv:2610.06525v3 Announce Type: replace
Abstract: Rotational transitions of SiC$_6$ with $j=26$--$40$ have been detected in IRC+10216, but collisional rate coefficients for SiC$_6$ with He or H$_2$ were previously unavailable. We present rate coefficients for SiC$_6$ rotational excitation by He, calculated using a two-dimensional potential energy surface (PES) at the RCCSD(T)-F12 level under the rigid-rotor approximation. Close-coupling (CC) calculations were performed for low rotational levels, while the coupled-states (CS) method was used for higher levels. Convergence tests were conducted, and CC and CS scattering cross sections were compared. The CS rate coefficients cover rotational levels up to $j=40$ and temperatures of 3--60\,K. Although CC and CS rate coefficients differ by up to a factor of $\sim 2$ for some transitions, particularly those affected by weak coupling and resonances, they yield nearly identical excitation and brightness temperatures for the high-$j$ transitions examined, including $19\rightarrow18$. These results support the reliability of the CS method for non-LTE modelling of the high-$j$ transitions considered, for which CC calculations are computationally impractical. We use reduced-mass scaling of the SiC$_6$--He rate coefficients to estimate those for SiC$_6$--H$_2$ and model the observed transitions. A $\chi^2$ analysis of the observed line fluxes, excluding heavily blended features, yields a best-fitting of $n(\mathrm{H_2})=4.33\times10^5$\,cm$^{-3}$ and $N(\mathrm{SiC_6})=8.11\times10^{11}$\,cm$^{-2}$. The latter agrees excellently with the previous observational estimate of $(7.9\pm1.5)\times10^{11}$\,cm$^{-2}$, supporting the consistency of our non-LTE modelling with observational constraints.
A faint dust ring beyond the compact disk of Sz 76
oai:arXiv.org:2610.09881v2
arXiv:2610.09881v2 Announce Type: replace
Abstract: We report new ALMA Band 7 high-resolution ($\sim$0.05$^{''}$) continuum and $^{12}\text{CO}\ (J=3\text{--}2)$ observations of the disk around the young low-mass star Sz 76. Visibility-plane modeling reveals a previously undetected, faint outer dust ring centered at $\sim$45 au surrounding the compact central disk. Including this low-surface-brightness feature increases the characteristic dust radius from $R_{\text{dust,68}} \approx 7$ au to $41$ au. Combined with an outer gas extent of $R_{\text{CO,68}} \approx 77$ au, the resulting gas-to-dust size ratio decreases from an extreme value of $\sim$11 to $\sim$2, placing Sz 76 in line with typical protoplanetary disks and contradicting a scenario of unimpeded inward radial drift. These results demonstrate that spatial filtering or limited surface-brightness sensitivity can severely underestimate disk sizes, highlighting the necessity of multi-scale observations to properly identify dust traps and accurately constrain disk evolution.
Structure and evolution of cosmic voids within DESI DR1
oai:arXiv.org:2610.18348v2
arXiv:2610.18348v2 Announce Type: replace
Abstract: Ongoing stage IV galaxy surveys such as DESI or Euclid offer an unprecedented amount of information to study the large-scale structure of the Universe by mapping tens of millions of galaxies. Within the context of cosmic voids, already released and future catalogues will help better understand the void structure and evolution as well as the cosmological model constraints derived from their statistics. With the recent DESI Data Release 1, our aim is to obtain a new void catalogue that allows us to study the properties of cosmic voids at different redshifts, hence providing insight not only on their structure but also on their evolution.
We select the DESI BGS and LRG samples for this study. Before the void-finding process, the data goes through a preprocessing pipeline tackling redshift space distortions and reconstructing the mildly non-linear velocity field. After that, the new AVISM void finder is applied to the preprocessed data, obtaining two void catalogues in different redshift bins (0.1 < z < 0.3 and 0.4 < z < 0.8), both within a 1400^3 cMpc^3 cube. We successfully recover N_BGS = 1114 and N_LRG = 216 clean voids with sizes up to ~ 65 cMpc. The void size function, 3D shape variables, density profiles, and linear bias are obtained and analysed. Both the size function and the density profiles obtained are in agreement with previous theoretical studies. The expected mass depletion in void evolution is shown by the evolution of density profiles. We obtain a roughly constant linear bias inside r < Re with mean values of 1.31 +/- 0.02 and 2.14 +/- 0.04 for BGS and LRG, respectively. Both void catalogues are publicly released.
Optical constraints on intranight variability in TeV HBLs from a two-year Asiago Schmidt monitoring campaign
oai:arXiv.org:2610.20065v2
arXiv:2610.20065v2 Announce Type: replace
Abstract: We present a two-year homogeneous optical monitoring campaign of three TeV-emitting high-synchrotron-peaked BL Lac objects, 1ES 1959+650, H 1426+428, and PG 1553+113, carried out with the Asiago Schmidt 67/92 cm telescope. TeV blazars are key targets in high-energy astrophysics, relevant for relativistic jet physics, particle acceleration, multiwavelength variability studies, and, in selected cases, multimessenger programmes. In this context, dense optical monitoring provides a low-energy synchrotron tracer of jet activity over temporal baselines that are difficult to cover continuously with high-energy facilities. The dataset consists of more than 17,000 R- and I-band photometric measurements obtained over 957 source-filter observing sessions between 2024 February 1 and 2026 January 31, for a total monitoring time of about 420 h. We applied a uniform differential-photometry and classical F-test framework based on two independent comparison-star pairs. Robust IDV is detected only in H 1426+428, with 15 out of 43 analysed sessions classified as variable in the R band and 5 out of 40 in the I band. No significant IDV is detected in 1ES 1959+650 or PG 1553+113 within the adopted statistical framework. For H 1426+428, the corresponding duration-weighted duty cycles, computed using only variable sessions with a measured local characteristic timescale, are 28.0% in the R band and 7.0% in the I band. The inferred local extrapolated doubling/halving timescales are of the order of 5-15 min in the observer frame. These results indicate that rapid optical synchrotron variability in TeV HBLs is intermittent and state-dependent. They also show that robotic small- and medium-aperture optical facilities can provide valuable long-term support for high-energy monitoring programmes by supplying dense, uniform, and repeatable optical coverage.
Consistently Reproducing Early Galaxy Observations and 21-cm Signal during the Epoch of Reionization
oai:arXiv.org:2610.24064v2
arXiv:2610.24064v2 Announce Type: replace
Abstract: In the standard cosmological model, galaxies forming within dark matter haloes are the primary drivers of the Epoch of Reionization (EoR). The complex physics of galaxy formation introduces scatter in the relation between halo mass and star-formation activity, with implications for the production and escape of ionising photons. In this work, we implement a stochastic source model in our cosmological radiative transfer simulations. We calibrate our model against the high-redshift ultraviolet luminosity function observed by JWST and observational constraints on Lyman-continuum leakage to constrain the ionising escape fraction. We find that stochastic star formation redistributes ionising emissivity toward lower-mass haloes, producing a more fragmented ionisation topology and delaying the growth of the largest ionised regions relative to the deterministic, observationally calibrated model. However, when the two models match total ionising photon production, their global neutral fractions and 21-cm power spectrum remain nearly indistinguishable, despite small differences in the morphology of individual ionised regions. Additionally, we establish a unified framework that lets us apply source models derived from computationally expensive hydrodynamical simulations (SPICE, FIRE-2 and Thesan) and investigate their impact on large-volume reionisation simulations, extending their prescriptions to scales of >100 Mpc. We find that these prescriptions produce distinct ionisation topologies and 21-cm statistical signatures, even when their global reionisation histories are similar. Our results highlight the importance of combining physically motivated galaxy formation prescriptions from hydrodynamical simulations with large-volume radiative transfer simulations to interpret forthcoming 21-cm observations and constrain the sources of cosmic reionisation.
Beyond the Disc--Spheroid Dichotomy at $0.6 < z < 8.0$: Evidence for a Morphological Indeterminacy Regime from JWST
oai:arXiv.org:2610.24078v2
arXiv:2610.24078v2 Announce Type: replace
Abstract: Galaxy morphology is a fundamental tracer of galaxy assembly, but its characterization becomes increasingly uncertain in the early Universe. We investigate morphological classifications across different methodologies over the redshift range $0.6 < z < 8.0$, using JWST observations of 15,567 galaxies with 6,099 synthetic galaxies from IllustrisTNG. To quantify morphology without relying on visual labels or prescribed light profiles, we use an unsupervised-supervised hybrid machine-learning approach based on the non-parametric MEGG system, which combines the second moment of the light distribution, Entropy, Gini, and Gradient Pattern Analysis. We find that parametric and non-parametric classifications remain broadly consistent at $z < 2.5$, but diverge systematically at higher redshift. We identify this epoch as the Morphological Indeterminacy Regime (MIR), in which classical disc and spheroidal categories no longer provide a unique description of the galaxy population. The fraction of mixed-interpretation galaxies increases with redshift, reaching $F_{\rm MIR} \sim 0.55$ at $z > 2.5$. Cross-domain and projection-invariance tests show that the classifications are stable and insensitive to data set and viewing angle. In the interval $2.5 < z < 8.0$, MIR galaxies occupy an intermediate stellar-mass regime (median; $\log(M_\star/M_\odot) \sim 9.7$), between low-mass discs (median; $\log(M_\star/M_\odot) \sim 9.5$) and massive spheroids (median; $\log(M_\star/M_\odot) \sim 10.1$), while irregular systems have stellar-mass distributions similar to those of MIR galaxies. We interpret these results as evidence that morphological indeterminacy is a characteristic feature of galaxy assembly in the early Universe, reflecting a population of dynamically young systems that have not yet settled into the classical Hubble sequence.
The DARWIN Project: Cosmological Simulations of Galaxy Formation down to Dwarf Scales in Diverse Environments
oai:arXiv.org:2610.25786v2
arXiv:2610.25786v2 Announce Type: replace
Abstract: We present the DARWIN project, a suite of cosmological simulations with RAMSES to study dwarf-galaxy formation and evolution across diverse environments. DARWIN-1 is a full-box simulation in a cubic volume of $(65.5\,{\rm cMpc})^3$ with a minimum cell size of 500pc, while DARWIN-2 comprises zoom-in simulations of 24 field regions and a Milky Way analog connected to a filament, reaching 62.5pc resolution. The simulations employ baryonic physics tailored to their aims. For DARWIN-2, we include radiative transfer in three photon energy bands to trace hydrogen and helium ionization. DARWIN-1 reproduces the cosmic star formation rate density and galaxy stellar mass functions reasonably well at $z>1$. Both simulations broadly reproduce observed stellar mass-halo mass relations, chemical abundances, star-formation surface densities, and black hole-galaxy scaling relations. At lower redshift, DARWIN-1 shows an excess in the number density of massive galaxies, although the inferred stellar content of massive halos depends on the definition of galaxy stellar mass. Other systematic differences arise from spatial resolution, feedback efficiency, and environmental selection. To facilitate comparison with observations, we generate mock images using the radiative transfer codes SKIRT and Powderday. We also explore a machine-learning approach to mitigate resolution-driven biases in the dwarf galaxy properties of DARWIN-1 using the higher-resolution DARWIN-2 simulations. The DARWIN simulation suite therefore provides a comprehensive numerical framework for studying the formation, evolution, and internal physics of dwarf galaxies in a cosmological context.
The Central Molecular Zone of our Galaxy, a chemical complexity factory: first interstellar detections of molecules towards G+0.693
oai:arXiv.org:2610.26308v2
arXiv:2610.26308v2 Announce Type: replace
Abstract: The Central Molecular Zone (CMZ) of our Galaxy is a prime laboratory for interstellar chemistry because of its extreme conditions, such as shocks and high cosmic-ray ionisation rates. It is an efficient factory of molecules, with many species first discovered in CMZ sources. This document reviews key results from the study of the CMZ G+0.693-0.027 molecular cloud. Using ultra-sensitive broadband spectral surveys together with laboratory spectroscopy, more than 140 molecules have been identified, including 28 detected for the first time in the interstellar medium (ISM). These discoveries include prebiotic precursors of biomolecules, increasingly complex members of known molecular families, molecules containing periodic third-row elements (S, P, Na, Mg), high-energy stereo- and structural isomers, and new cations. Overall, these findings extend the known limits of interstellar chemical complexity and highlight the CMZ's potential to provide the molecular ingredients for prebiotic chemistry.
New Bounds on Heavy QCD Axions from Big Bang Nucleosynthesis
oai:arXiv.org:2510.23695v3
arXiv:2510.23695v3 Announce Type: replace-cross
Abstract: We study Big Bang Nucleosynthesis (BBN) constraints on heavy QCD axions. BBN offers a powerful probe of new physics that modifies the neutron-to-proton ratio during the process, thanks to the precisely measured primordial Helium-4 abundance. A heavy QCD axion provides an attractive target for this probe, because not only is it a well-motivated hypothetical particle by the strong CP problem, but also it dominantly decays to hadrons if kinematically allowed. A range of its lifetime is thus excluded where the hadronic decays would significantly alter the neutron-to-proton ratio. We compute axion-induced modification of the neutron-to-proton ratio, and obtain robust upper bounds on the axion lifetimes, as low as 0.017 s for the axion mass higher than 300 MeV. Remarkably, this is stronger than projected future CMB bounds via $N_{\rm eff}$. Our bounds are largely insensitive to uncertainties in hadronic cross sections and the axion's branching fractions into various hadrons, as well as to the precise value of the initial axion abundance. We also incorporate, for the first time, several key improvements, such as scattering processes by energetic $K_L$ and secondary hadrons, that can also be important for studying general hadronic injections during BBN, not limited to those from axion decays.
Inflationary Particle Production and the Swampland
oai:arXiv.org:2512.07930v3
arXiv:2512.07930v3 Announce Type: replace-cross
Abstract: We investigate the impact of particle production during inflation in scenarios where an infinite tower of states features a mass scale that decreases exponentially along the inflationary trajectory. Such couplings naturally arise in string effective field theories and are in fact motivated by the Swampland Distance Conjecture (SDC). We show that the corrections to inflationary observables sourced by the tower scale as $(H/\Lambda_{\text{sp}})^{2+p}$, with $H$ being the Hubble scale, $\Lambda_{\text{sp}}$ being the species scale, that is the quantum gravity cut-off, and $p\geq 1$ characterizes the density of states in the tower. As a result, in gravitationally weakly coupled cosmological effective theories, the tower-induced contributions are suppressed relative to the standard single-field predictions, leaving the inflationary phenomenology essentially unchanged. We demonstrate this explicitly across a set of well-motivated inflationary potentials, and we compare the resulting predictions with the most recent observational constraints, including those from the Atacama Cosmology Telescope.
Probing freeze-in dark matter using Bose-Einstein condensate in neutron star
oai:arXiv.org:2605.22245v2
arXiv:2605.22245v2 Announce Type: replace-cross
Abstract: Neutron star (NS) is one of the most promising astrophysical targets to probe non-gravitational interaction of dark matter (DM) with visible matter. Their compactness makes them an ideal object which can capture particle DM efficiently over its lifetime using the DM-nucleon scattering cross-section. If DM particles are bosonic, then the captured DM population may form a Bose-Einstein condensate at the center of the NS, increasing the DM density significantly. In this work, we study the phenomenology of such scenario with enhanced DM annihilation rate due to the increased density in a condensate. The enhanced DM annihilation makes the NS surface `hotter' than in the standard cooling scenario. We show that the annihilation rate is enhanced by a factor of $\mathcal{O}(10^{15}-10^{20})$ if DM forms a condensate without large self-interaction, and DM with freeze-in value annihilation cross-section can heat up the NS to higher temperatures, bringing it within the reach of James Webb Space Telescope. It also allows us to probe DM-nucleon scattering cross section within the neutrino fog regime which will complement the terrestrial direct detection searches. Moreover, the enhanced annihilation from the condensate changes the lower limits on s-wave DM annihilation cross-section for capture-annihilation equilibrium and the formation of a black hole inside the NS. Finally, we show an example of a scalar DM model where such small annihilation and DM-nucleon scattering cross sections can generically arise.
Energy extraction from NED-deformed rotating black holes via the Comisso-Asenjo reconnection process
oai:arXiv.org:2605.26369v2
arXiv:2605.26369v2 Announce Type: replace-cross
Abstract: We study rotating black holes in general relativity coupled to nonlinear electrodynamics (NED), focusing on an axisymmetric solution with deformation parameter $g$. On the spherical seed, weak-field lensing via the Gauss-Bonnet method and the shadow radius yield a spin-insensitive bound by enforcing a conservative $\sim 10\%$ tolerance on the Sgr A* ring size, namely $g/M \lesssim 1.26$. In the eikonal regime we derive analytic quasinormal-mode shifts, even in $g$, and obtain an independent ceiling consistent with the shadow constraint. For the rotating geometry, we provide closed-form ZAMO scalars, chart horizons and the ergoregion, and analyze equatorial geodesics (photon orbits and ISCO). We then formulate the Comisso-Asenjo reconnection channel in the ZAMO frame, identify the negative-energy window, and integrate the extracted power over the allowed radii; from the tolerated fractional departure from the Kerr power we define a spin-dependent extraction bound $g_\delta(a\mid\sigma_0,\xi)$. Taken together, the QNM/shadow ceiling and the extraction bound appreciably narrow the admissible region for $g/M$ in the $(a/M,g/M)$ plane, so even within our deliberately simplified, single-layer equatorial setup, the two complementary probes already provide informative constraints on NED deformations, testable with present data and upcoming horizon-scale and ringdown campaigns.
Detecting Tidal Resonances in Binary Neutron Stars
oai:arXiv.org:2606.06376v2
arXiv:2606.06376v2 Announce Type: replace-cross
Abstract: As a binary neutron star inspirals due to the emission of gravitational waves, the rising tidal frequency resonantly excites vibrational modes. These oscillations are seismological probes of the rich stellar interior, yet it remains to be established whether gravitational-wave interferometers can measure them. Here, we present the first fully Bayesian study of the capability of the Einstein Telescope to detect tidal resonances. We simulate one year of observations and analyse the 200 loudest signals. We find that the Einstein Telescope can identify resonant modes and is sensitive to gravitational-wave phase shifts as small as $\Delta \Phi \approx 0.03$ for favourable events. We further show that neglecting resonances can bias the inferred tidal deformabilities. These results establish tidal resonances as a measurable route for asteroseismology with future detectors.
A dark-matter-sensitive orientation clock near rotating black holes
oai:arXiv.org:2606.06557v3
arXiv:2606.06557v3 Announce Type: replace-cross
Abstract: A coherent triaxial structure has an intrinsic orientation clock because rotation about its intermediate principal axis is unstable. I ask how this clock changes when the structure moves on a circular equatorial orbit near a rotating black hole in an effective dark matter (DM) environment. The calculation uses a comoving orthonormal tetrad, the electric part of the Riemann tensor, Euler--quaternion dynamics, and the conversion from proper time to stationary coordinate time. The clock shift separates exactly into a local proper-time response and a geometric time-conversion response. I compute a direct grid of 600 DM models for Einasto, regularized cored Navarro--Frenk--White, and Hernquist profiles. Across the scanned domain, the total fractional frequency shift lies between about $-3\times10^{-7}$ and $-3\times10^{-4}$. For the fiducial body model, most of the shift comes from the time conversion. Synthetic timing tests show that a shift near $10^{-4}$ is recoverable only when the clock is highly coherent and its Kerr baseline is known independently. This gives a theoretical timing diagnostic and a synthetic observational forecast. I do not fit real data.
Resonant lunar tides of Earth's core and basal magma ocean
oai:arXiv.org:2606.16340v2
arXiv:2606.16340v2 Announce Type: replace-cross
Abstract: Earth's magnetic field is generated by fluid motion in the liquid-metal core and has been active for billions of years. However, prior to the onset of inner-core growth, the sources of mechanical power that drove the geodynamo remain uncertain. During this period, the core may have been overlain by a basal magma ocean (BMO), creating two immiscible fluid layers separated by a density interface, beneath the solid mantle. We develop a theory for lunar tides in this core--BMO system, in which the tidal potential acts on the density contrast between the two layers. The resulting dynamics differ fundamentally from previous models driven by imposed deformation of the core's bounding envelope. In the inviscid limit, the response transitions from an equilibrium tide to an interfacial gravity wave as forcing frequency increases. The two regimes are separated by a resonance that occurs when the forcing frequency matches the natural frequency of the degree-2, order-2 interfacial mode. The inviscid core flow is formally identical to the canonical elliptical vortex, linking the problem to the theory of elliptical instability. Finite viscosity regularises the resonance, introduces phase lag and generates oscillatory boundary layers. Combined with parameterised models of lunar recession and BMO crystallisation, the theory predicts enhanced core-boundary ellipticity, core-flow speed, magnetic Reynolds number and instability metrics, particularly near resonance. These results identify a previously unexplored mechanism for tidally driven flow in differentiated planetary bodies and show that a BMO can enhance tidal coupling to the core, potentially contributing to dynamo action.
Gravitational-wave responses and noise correlations in geometric time-delay interferometry
oai:arXiv.org:2606.28173v2
arXiv:2606.28173v2 Announce Type: replace-cross
Abstract: Space-based gravitational-wave observations will probe the formation and growth of massive black holes and the evolution of compact binaries in the millihertz band. We develop a geometric-link representation of time-delay interferometry using six round-trip and six non-round-trip two-link blocks. Starting from the one-way Doppler response, we derive the elementary block responses and a dictionary that expresses 45 known second-generation combinations in terms of these blocks and residual one-way measurements. Representative Michelson and Relay observables admit compact forms, from which analytical sky-averaged tensor responses are obtained for a fixed, equal-arm constellation. We also derive the noise covariance of the blocks, retaining cross-spectra between blocks that share measurement streams. The resulting noise spectra agree with direct one-way assembly. We apply this formulation to a narrowband compact-binary signal in an equal-arm LISA instrument model. For the Relay example, neglecting the shared-link cross-spectrum changes the low-frequency noise scaling and overestimates the conditional Fisher uncertainties in amplitude or distance by factors of 10.71 at 1~mHz and 3.37 at 3~mHz.
Exact 1D Nonlinear Solutions for Proton-Driven Plasma Wakefields
oai:arXiv.org:2607.06458v3
arXiv:2607.06458v3 Announce Type: replace-cross
Abstract: We extend the classical 1D nonlinear relativistic plasma wakefield paradigm to evaluate the pristine fluid profiles driven by intense relativistic proton trains under the comoving coordinate approximation. Using an intense two-bunch scheme, our adaptive integration model is validated down to machine precision against exact closed-form cold-fluid invariants, confirming that interface boundary corners represent physical step discontinuities in the second potential derivative. Scaling to an $N=100$ micro-bunch train, we resolve macroscopic phase-accumulation and subsequent nonlinear detuning by deriving a closed-form, parameter-free piecewise analytical envelope expression. This solution successfully maps the initial unattenuated quadratic coherent growth alongside the flat asymptotic horizontal saturation horizon governed by nonlinear energy conservation. Furthermore, the framework analytically evaluates the transient envelope overshoot, tracking its decay onto the saturation boundary via an exact phase-slip integration length scale. Providing immediate, ready macro-envelope bounds, this computationally inexpensive tool eliminates empirical adjustments and yields a robust mathematical baseline designed to guide and accelerate complex, multi-dimensional particle-in-cell (PIC) simulations through rapid parameter space exploration before deploying heavy computational assets. It establishes an efficient foundation to optimize asymmetric, custom-shaped driver profiles required to break the fundamental symmetric transformer ratio limit of 2 to make current particle accelerators more efficient.
Entropic Chaplygin-Gas cosmology and late-universe tension diagnostics
oai:arXiv.org:2607.22802v2
arXiv:2607.22802v2 Announce Type: replace-cross
Abstract: We investigate an entropic generalized Chaplygin-gas (gCg) cosmology in which a single dark fluid accounts for both dark matter and dark energy. We constrain the model using DESI baryon acoustic oscillations, redshift-space-distortion measurements, compressed Planck 2018 distance and lensing information, and Type Ia supernova likelihoods with PantheonPlus+SH0ES (PPS) and DES--Dovekie. For Planck+DESI+RSD+PPS, the entropic gCg model yields a Hubble constant of $H_{0}=68.64\pm0.52\,{\rm km\,s^{-1}\,Mpc^{-1}}$ and the amplitude of matter fluctuations parameter of $S_{8}=0.797\pm0.021$, compared with $H_{0}=67.31\pm0.28\,{\rm km\,s^{-1}\,Mpc^{-1}}$ and $S_{8}=0.809\pm0.021$ in $\Lambda$CDM. The simultaneous upward shift in $H_{0}$ and downward shift in $S_{8}$ occurs in the directions required for a partial alleviation of both tensions. We also show that when PPS is replaced by DES--Dovekie, the parameter shifts are smaller, the best-fit improvement is negligible, and the Bayesian evidence favours $\Lambda$CDM. Thus, the entropic gCg model provides a potentially interesting simultaneous shift of the Hubble constant and the matter-clustering amplitude, but its statistical support depends on the adopted supernova likelihood.
Penrose process in radiating Kerr-Vaidya black holes: efficiency, energy transport, and the viability boundary
oai:arXiv.org:2609.19175v2
arXiv:2609.19175v2 Announce Type: replace-cross
Abstract: We study single-particle Penrose extraction from a rotating black hole whose mass decreases in time, described by the ingoing Kerr-Vaidya metric, imposing exact four-momentum conservation at the decay. We prove that the local efficiency is bounded by the instantaneous Kerr value, through a covariant bound that reduces to the Bardeen-Press-Teukolsky limit at the horizon, and we derive an exact transport law, dE/dtau = -dm/dt (r/Sigma) p_r^2 for fixed spin parameter, by which every timelike geodesic gains energy while the hole loses mass. The net energy gain therefore splits into a local Penrose term and two transport terms, which we compute numerically as functions of spin and mass-loss rate for two mass-loss laws. With fixed spin parameter, extraction is possible only below a critical rate alpha_crit(a) = (M-a)/Delta v, set by the approach of the hole to extremality; with fixed a/m no such boundary exists. The transport law extends to the Kerr-Newman-Vaidya family, and the regime where the dynamics matter is shown to lie far outside astrophysical mass-loss rates.