Most interesting arXiv papers
Friday 11-09-2026
http://arxiv.org/abs/2606.13359 (pdf) (52)
We analyse the cluster content of the ∼ 4.5 deg^2 XMM-LSS and CDFS VIDEO fields which are expected to be partially covered by the upcoming MOONRISE survey. Using AMICO and WaZP photometric redshift-based cluster finders, we construct a sample of 519 cluster candidates detected by both finders in the redshift range z = 0.1-3, including 74 detections at z > 1.5. For all detections, we identify the Brightest Central Galaxy (BCG) and compute a list of probabilistic cluster memberships. Our photometric redshift measurements of the clusters agree well with spectroscopic redshifts from the literature, when available. From ancillary spectroscopic data, we assign z_spec measurements to 116 cluster candidates based on their spectroscopic members and to 204 based on their likely BCGs. We also show that candidates containing Radio-Loud members are efficiently recovered using the prior-based cluster finder PPM. We perform a preliminary analysis of the galaxy content of these candidates, focusing on the Red-Sequence components of their apparent Colour-Magnitude Diagram. By comparing with models of galaxy evolution, we show that this population is consistent with a model of passive evolution with a formation at high redshift, and is already in place at z = 1.5-2.0. Finally, our cluster sample is used to evaluate how these clusters would be detected and characterised, according to various MOONRISE strategies. We show that cluster spectroscopic confirmation and characterisation could be efficiently achieved up to z∼1.7 even with the shallowest survey strategy. This open unprecedented insight into the physical properties of high-redshift galaxy clusters and into galaxy formation in dense environments.
http://arxiv.org/abs/2605.31052 (pdf) (52)
Early (z ≳ 2) Massive (M_⋆≳ 10^10 M_⊙), Quenched Galaxies (MQGs) challenge current galaxy formation models. In this series, we study these systems using the new COLIBRE cosmological hydrodynamical simulations. Following the broad agreement between their predictions and observations found in the first paper, this second paper explores the processes driving galaxies to become massive and quenched in COLIBRE, identifying Active Galactic Nucleus (AGN) feedback as the primary quenching mechanism in both the thermal (L200m6 simulation) and hybrid (thermal+jet, L200m7h simulation) AGN feedback models implemented. However, the two models behave differently: while the thermal model efficiently quenches massive galaxies at z>3, the hybrid model is less effective because black holes (BHs) grow more slowly in the early Universe, and the jet component, which dominates the feedback energy, acts on longer timescales to impact galaxies. Both models predict quasar-like MQGs (AGN with L_bol≳10^45 erg s^-1), with the most luminous systems associated with more recently quenched galaxies. Compared to star-forming galaxies of similar mass, MQGs host more massive BHs and exhibit higher star formation efficiencies. These differences arise primarily from their environments before quenching, particularly at local (0.3 cMpc) to intermediate scales (1.0 cMpc), where overdense regions are associated with enhanced gas inflows, higher BH accretion and, hence, feedback power. We find that about 54
http://arxiv.org/abs/2609.10803 (pdf) (45)
JWST has uncovered a class of objects called LRDs, whose nature is still widely debated. In this work, we present a comprehensive spectroscopic analysis of nine LRDs in the Extended Groth Strip (EGS) field studied as part of THRILS and C3PO, both JWST Cycle 3 programs. These targets, photometrically selected based on their compact red appearance, are observed with deep spectroscopic exposures (≥8 hours), enabling robust detections of broad Balmer lines, He I emission, and other spectral features characteristic of AGN activity. Using the [SII] λλ6716,6731 doublet, we find electron densities (n_e) between 2.33 < log (n_e) < 2.97 cm^-3, comparable to those in narrow line regions (NLR) of local AGN and high-z galaxies. The spectroscopic depth further enables detailed characterization of broad Balmer line profiles. We fit both Gaussian and convolved exponential models to each source and find that five LRDs are statistically better described by the latter model. We measure optical depths τ_sc = 0.56-0.91, scattering fractions f_SC = 0.40-0.82 which correspond to column densities log(N_e) ∼ 23.93-24.14 cm^-2, and covering fractions c_f = 0.43-0.59. These results indicate a clumpier broad line region (BLR) geometry that deviates from conventional LRD models, which predict covering fractions close to unity. Furthermore, these Compton-thick gas columns may explain the X-ray weakness of LRDs. We also find that THRILS LRDs are narrow-line dominated compared to literature AGN-dominated LRDs, and show that exponential profile fitting corrects for systematic overestimation of black hole masses from Gaussian-based measurements.
http://arxiv.org/abs/2609.11827 (pdf) (40)
We present the first spectral analysis from the JWST Cycle 3 GO Program #6362, the Dyn-Atmo Survey of dynamical-atmospheric benchmarks, characterizing the atmospheres of directly-imaged companions with dynamical masses and measured orbits. This program has obtained moderate-resolution spectra (R∼2700, 3-5μm) using the NIRSpec/G395H IFU. We observe the T4-dwarf companion, HD 13724 B, companion to a G-type star at a projected angular separation of 0”2. Updated orbital monitoring constrains the mass to 38.3±0.6 M_Jup, which we use as a Gaussian prior when inferring the atmospheric state. We fit the continuum-subtracted spectrum using the BT-SETTL, Sonora Diamondback, Sonora Flame Skimmer, and NEWERA-PHOENIX self-consistent model grids, together with petitRADTRANS atmospheric retrievals. We adopt bulk parameter values of T_eff=1150±195 K, log g=5.1±0.3 dex, [M/H]=0.4±0.4, log_10(K_zz)=6.33_-2.22^+2.59, and C/O=0.50±0.13 using a bayesian average across model grids. We find degeneracies in the models between surface gravity, metallicity, and K_zz for Sonora Diamondback and Flame Skimmer. The Flame Skimmer and NEWERA-PHOENIX vertical mixing strengths are incompatible, with Flame Skimmer finding a systematically lower value. The retrievals favor an atmosphere in chemical equilibrium, with enhanced metallicity and a solar C/O ratio, but with no detection of disequilibrium chemistry, and no detection of trace gases other than CH_4, CO_2, ^12CO and H_2O. Ultimately, we identify critical challenges in forward modeling continuum-subtracted G395H spectra in the low signal-to-noise regime. We find that including independent dynamical mass measurements better constrains atmospheric parameters, but degeneracies between self-consistent model parameters limit accurate atmospheric parameter inferences.
http://arxiv.org/abs/2609.11094 (pdf) (38)
We present the first ultra-deep rest-UV spectroscopy of four UV-bright Little Red Dots (LRDs), obtained from the SPURS Cycle 4 Large Program. The spectra reveal broad CIV (FWHM ≈2700-2800 km s^-1) in two LRDs, alongside narrow-line densities elevated above star-forming galaxies (n_e∼10^4-10^5 cm^-3, reaching 10^6 cm^-3 in the most extreme source) and nitrogen-enhancements in all four LRDs. We detect broad HeII emission (FWHM ≈930 km s^-1) in one LRD, and two others with fast P-Cygni absorption (≳2200 km s^-1). Strong interstellar absorption lines and Lyα damping wings reveal the UV continuum is deeply embedded in neutral gas (N_HI≳10^22 cm^-2) in all four LRDs. Detections of fluorescent FeII and OI emission and fine-structure absorption indicate this gas lies close to the UV-emitting region. In archival z>4 samples, we find nitrogen and strong CIII] emission are significantly more common in LRDs than in the galaxy population. The transmission of broad CIV, tracing the broad-line region or cocoon, depends on rest-optical color within our sample, consistent with an orientation-dependent picture in which bluer, less obscured sightlines offer a more direct, polar view of the central engine and its outflows. We find several potential signatures of very massive stars, whose winds may contribute to nitrogen enhancement. We investigate other abundance patterns expected from supermassive stars but our results are inconclusive. Our results place the UV-emitting region within ≲8 pc of the LRD nucleus, consistent with an actively assembling nuclear star cluster. Dynamical interactions in this extremely dense environment, including tidal disruption of stars, may explain the high incidence of nitrogen enhancements in LRDs.
http://arxiv.org/abs/2606.12254 (pdf) (37)
Studying young stellar populations within galaxies can help refine our understanding of recent star formation in galaxies and their evolution. With this motivation, we present a catalog of 25,000 recently formed (within 400 Myr) star-forming clumps (SFCs) in 17 morphologically diverse nearby galaxies, including 8 massive, classic spirals, 6 intermediate-mass, flocculent spirals, and 3 dwarf irregulars. We used far- and near-UV observations from the UltraViolet Imaging Telescope (UVIT), whose 1.5" angular resolution and 28&apos field-of-view allow us to probe SFCs at a mean physical scale of 54 parsec, within the full extent of our galaxies. We adopted a homogeneous SFC detection criterion, corrected for spatially varying dust attenuation (using 6" resolution A_V maps, made by combining FUV with archival infrared observations), and estimated the SFC ages by comparing the observed UV color-magnitude diagrams with Starburst99 simple stellar population models. Using our SFC catalog, we studied the age demographic of the recently formed stellar populations across different galaxy morphologies and observed age trends consistent with several well-known phenomena, such as the inside-out formation of disc galaxies, local gravitational instabilities leading to flocculent spiral arms, and the stochastic nature of star formation in dwarf galaxies. Leveraging full galaxy coverage and FUV data, our catalog complements existing optically-identified star cluster catalogs in the literature towards improving our understanding of star formation across a wide range of galaxy morphologies, masses, and environments. We make the SFC catalog and A_V maps of our 17 galaxies publicly available with this paper.
http://arxiv.org/abs/2609.11879 (pdf) (36)
We present FM-JADES-v1, a self-supervised foundation model for James Webb Space Telescope (JWST) deep-field science, trained with 482,444 objects from the JWST Advanced Deep Extragalactic Survey (JADES) Data Release 5 using multi-band imaging and the photometric catalog. The shared embedding space is trained without class labels. We demonstrate that FM-JADES-v1 can serve as a powerful tool for object discovery and improving property measurements using two experiments, blind active discovery and few-band photometric redshift. For blind object discovery, FM-JADES-v1 identifies rare object populations such as high-redshift galaxies and Little Red Dots (LRDs) without any prior population labels or population-specific selection criteria. These rare populations emerge as isolated islands in the embedding space, which can be identified without prior astrophysical knowledge. For few-band photometric redshift, FM-JADES-v1&aposs learned embeddings achieve σ_NMAD=0.157 in a strictly controlled three-band (F115W/F200W/F356W) photo-z benchmark, compared to σ_NMAD= 0.44 for template fitting. These results demonstrate the potential of self-supervised multi-modal representations as scalable discovery spaces for large astronomical surveys. Applied to ongoing and future wide-field surveys from JWST, Roman, Euclid, and Rubin/LSST, this framework could enable systematic searches for rare populations, as well as enabling multiple downstream tasks such as improving astrophysical property measurements.
http://arxiv.org/abs/2609.11828 (pdf) (34)
We present new JWST/NIRSpec IFU high contrast spectroscopic observations of the brown dwarf companion GJ 758 B. Extensive radial velocity monitoring of the primary has enabled high-precision measurements of the companion mass through characterization of their joint orbital dynamics. The combination of prior space- and ground-based photometry and spectroscopic observations in concert with a spectral model-independent mass makes this target an ideal benchmark for calibrating evolutionary models with known mass-age-luminosity degeneracies. A template-matching analysis on the NIRSpec/IFU data (F290LP,G395H; λ∈ [2.87-5.27] μm) detects the companion at the expected orbital position with incredible significance (265.2σ). Comprehensive spectral analysis including published datasets over a variety of different spectral model families reveal that an inference using a uniform radius prior is plagued by inaccurate posteriors. The atmospheric model degeneracy between metallicity, gravity, and clouds pull the inferred mass/gravity/radius to an unphysical portion of parameter space for all model families tested. A constrained inference with a gaussian prior on the mass can mitigate this issue but models without clouds still struggle to accurately reproduce the K-band photometry and result in mass inferences which are inconsistent with the dynamics. We explore modifying the spectral models with custom post-processed clouds and molecular abundance perturbations of ammonia and hydrogen sulfide. We find the best-fitting model includes clouds with no patchiness (hole fraction f_hole = 0), a small fractional depletion in ammonia and moderate enhancement in hydrogen sulfide.
http://arxiv.org/abs/2609.11654 (pdf) (33)
We analyse the neutral ISM of galaxies at redshift z=0 in the COLIBRE hydrodynamical simulations and investigate the conditions under which these two gas phases coexist in a narrow range of thermal pressures. COLIBRE galaxies are selected based on the metallicity of their ISM. The median and mass-weighted distributions of thermal pressures of the multiphase ISM are analysed and compared to thermal equilibrium models and observations. The ISM in galaxies with gas metallicities similar to solar values exhibits a clear multiphase structure with a warm and cold phase coexisting in a certain range of thermal pressures. The pressures at which the ISM is multiphase depend on the gas metallicity. For COLIBRE galaxies with lower metallicities (Z_ISM≲0.1 Z_⊙), this multiphase structure largely disappears, partly due to resolution. The thermal pressures weighted by the HI mass of the neutral phases in COLIBRE galaxies are lower than, but still comparable to, some theoretical works and observational estimates. The thermal pressures show a strong dependence on the weighting scheme. If weighted by the star formation rate or CI mass, the thermal pressures of the cold phase match those derived from observations. The resulting pressures depend on a combination of the assumed radiation field strength, dust abundance, limited resolution and the weighting scheme. The strong dependence of the thermal pressures on the weighting scheme indicates observational tracers used to estimate the thermal pressure are biased towards high-pressure regions.
http://arxiv.org/abs/2609.11191 (pdf) (32)
We investigate the co-evolution of supermassive black holes (SMBHs) and central galaxies in massive clusters using the GIZMO-SIMBA and TNG-Cluster zoom-in simulations at z=0-5. We find that the distinct subgrid physics of these two models suggest fundamentally different evolutionary pathways. On the one hand, GIZMO-SIMBA, employs torque-limited accretion and predicts a supply-driven scenario where the SMBHs rapidly assemble synchronized with dark matter halo (M_200c) growth (i.e. the halo mass-BH mass relation is set by z=3.0 and similar to the present-day relationship). On the other hand, TNG-Cluster, exhibits a feedback-regulated growth phase delayed by an early thermal suppression. While both models successfully reproduce some local black hole-galaxy scaling relations, they imply significantly different evolution for these relations. Analysis of the BH mass-gas mass ratio relations suggests that TNG-Cluster&aposs isotropic kinetic winds efficiently deplete cold gas, resulting in a "hard quench" of star formation. In the black hole accretion rate (BHAR)-star formation rate (SFR) relation we find that both simulations successfully reproduce the decoupling of BHAR and star formation observed in recent massive cluster ellipticals. The divergent evolutionary trends emphasize the importance of the multiphase intracluster medium; while these subgrid models do not have the necessary resolution and employ distinct formalisms, the sustained BHAR in quenched systems resemble outcomes broadly consistent with modern multiphase feeding paradigms such as chaotic cold accretion in turbulent cluster cores. Furthermore, we demonstrate that for both models, black hole mass is a primary regulator of atomic and molecular gas depletion in galaxy clusters.
http://arxiv.org/abs/2609.11497 (pdf) (32)
Rest-frame ultraviolet emission lines are diagnostics of hard ionizing spectra and highly ionized gas in galaxies resembling the sources of cosmic reionization. Nebular CIV emission traces both intrinsic line production and resonant scattering through its velocity-resolved profile. Previous studies have treated CIV empirically, using integrated fluxes, equivalent widths, and line ratios rather than modeling the radiative transfer (RT) shaping the emergent profiles. We present the first systematic RT modeling of resonantly scattered CIV emission profiles in 18 local star-forming galaxies observed with HST/COS. Using the clumpy RT framework PEACOCK, idealized experiments show that profile morphology is governed by column density, gas kinematics, intrinsic equivalent width, and aperture-dependent recovery of scattered emission. We fit continuum-normalized spectra corrected for stellar contributions with a neural-network-accelerated Bayesian pipeline, reproducing P-Cygni-like, double-peaked, double-peaked with two absorption troughs, and blue-bump plus red-peak morphologies. Emission infilling connects these profile classes and complicates kinematic inferences, while aperture losses decouple observed net equivalent widths from intrinsic line production. Combining RT-inferred intrinsic equivalent widths with rest-UV diagnostics, we find that most galaxies are consistent with stellar photoionization, although some may require harder spectra or more extreme ionization conditions. Comparison with CLASSY reveals stronger associations of stellar mass and star formation rate with gas kinematics than with CIV column density. Our results establish CIV as both a tracer of hard ionizing radiation and a resonant-line probe of gas structure and kinematics, with velocity-resolved profiles providing information beyond integrated line strengths in local reionization analogs.
http://arxiv.org/abs/2609.11328 (pdf) (30)
We constrain the particle mass of the fuzzy dark matter model using globular cluster candidates in dwarf galaxies, from the Euclid Early Release Observations, within the Fornax and Perseus galaxy clusters. We model the orbital evolution of globular clusters by combining fuzzy dark matter density profiles with the corresponding dynamical friction formalism, and integrate their orbits using galpy. In the fuzzy dark matter framework, both the halo density profile and the dynamical friction force are modified by the wave-like nature of the ultra-light dark matter particle, whose mass is parameterized as m_22=m_χ/10^-22,eV and explored over the range 0.1 ≤ m_22≤ 100. This framework naturally alleviates the long-standing globular cluster timing problem originally identified in the Fornax dSph. We find that globular clusters no longer experience orbital decay toward the centers of their host galaxies for m_22≤1.46^+0.51_-0.32 in the Fornax cluster and m_22≤1.09^+0.14_-0.21 in the Perseus cluster. Future Euclid data releases, providing substantially larger samples of globular cluster candidates, will enable significantly tighter constraints on the dark matter particle mass within the fuzzy dark matter scenario.
http://arxiv.org/abs/2609.11104 (pdf) (30)
We present the largest sample to date of nitrogen-loud (N-loud) quasars with strong broad N IV] λ1486 and/or N III] λ1750 emission lines over the redshift range 1.6 < z < 4.3, selected from the Dark Energy Spectroscopic Instrument (DESI) Data Release 1. The final sample contains 1,993 N-loud quasars, corresponding to about 1.2% of the parent quasar sample. The L_1450 distribution of the N-loud quasars is broadly similar to that of the DESI parent sample, but their redshift distribution is distinct, with a stronger concentration around z ∼ 2.5–3. Their composite spectrum displays a broadly similar UV continuum shape to that of the parent quasars, while showing significantly enhanced broad nitrogen emission features, including N V, N IV], and N III]. Other metal emission features also show a moderate enhancement. Relative to a control sample matched in redshift and UV continuum luminosity, the N-loud quasars show systematically narrower broad C IV and Mg II emission lines, lower single-epoch virial black hole masses, and higher Eddington ratios, suggesting that N-loud quasars may preferentially appear during a relatively rapid black hole accretion phase. The radio-loud fraction is 10.1%, with the highest fraction among objects exhibiting both N III] and N IV] emission. The catalog provides a statistical baseline for future studies of nitrogen enhancement and its physical origin.
http://arxiv.org/abs/2609.11510 (pdf) (29)
Extensive Euclid satellite imaging at low surface brightness has revealed that most nearby galaxies are diffuse-looking spheroids, where the stellar radius increases monotonically over three decades in luminosity. We argue this Diffuse Galaxy sequence results from internal stellar diffusion by Wave Dark Matter (ψDM), as wave energy is transferred to star orbits over time. In particular, the soliton random motion scatters central stars onto radial orbits that become enhanced with each passage through the centre, slowly "puffing up" the stellar profile. Heating is greater within massive galaxies as ψDM fluctuations are stronger and more frequent, reproducing the Diffuse Galaxy sequence and also accounts for the rising velocity dispersion along the sequence, from Ultra-Faint to Dwarf Spheroidal and Ultra Diffuse galaxies, favouring a light boson, m_ψ=2.88^+0.14_-0.13×10^-22eV. Winding back this diffusion, we predict the stellar content of Diffuse Galaxies, including globular clusters, formed near the centre, as anticipated by ψDM simulations, where gas cools efficiently within the dense soliton. This predicted ψDM evolution from compact beginnings towards diffuse-looking spheroidal galaxies today can now be fully charted from JWST to Euclid.
http://arxiv.org/abs/2609.11465 (pdf) (27)
Deeply bound stellar substructures of about 1 kpc at the center of bulges are found in a subset of galaxies from the cosmological chemo-hydrodynamical zoom-in CIELO simulation suite. They were identified as stellar overdensities at the lowest binding energies in each galaxy&aposs circularity-energy (ε, E) plane. We refer to these overdensities as stellar knots. We aim to characterize their properties in the context of bulge assembly using CIELO simulated galaxies spanning a wide range of stellar masses (10^8.0-10^10.7 M_⊙), with diverse formation histories. We inspect 54 galactic bulges. Within them, we isolate stellar knot candidates in (ε, E) space, and identify 28 robust knots satisfying successive selection criteria: kinematic, concentration, and morphology. We characterize their chemical enrichment, formation timescales, formation sites, progenitor gas origin, and spatial distributions. Across all galaxy masses, knots are systematically alpha-element enhanced, having assembled the bulk of their stellar mass at earlier epochs and on shorter timescales than other bulge populations, with a median of ∼2Gyr versus ∼5Gyr. Regarding their origin, knots are predominantly in-situ (5
http://arxiv.org/abs/2609.10868 (pdf) (27)
Down-the-barrel absorption line spectra are a powerful probe of the properties of multiphase galactic outflows. In this work, we demonstrate a method to generate realistic synthetic spectra using high resolution (Δx=5-20 pc), 20-kiloparsec-scale outflow simulations from the next generation of the CGOLS project. Aiming to mimic modern observational surveys of nearby star-forming galaxies, we generate a suite of UV absorption lines, from which we measure common observables such as equivalent widths and velocity statistics, and compare them to observed values. The velocity statistics from the simulated data agree well with high-resolution, far-UV HST treasury data, but the equivalent widths (EWs) display more variability. Our fiducial EWs are systematically lower than the empirical predictions by ∼50-97
http://arxiv.org/abs/2609.10672 (pdf) (27)
Within the Hot Big Bang picture, as usually interpreted in the standard LCDM framework and its hierarchical scenario for elliptical galaxy assembly, the cosmic microwave background (CMB) is associated with photon decoupling near z = 1100, and the Planck angular power spectra are taken to constrain primordial 10^-5-level fluctuations that seeded later structure formation. We show that, in addition to the LCDM model, two cosmological models based on Milgromian gravitation are also consistent with the Planck CMB power spectra, such that the spectra do not uniquely favor dark matter. We then argue that the published CMB spectra are not direct observables, but reconstructed quantities obtained after foreground subtraction and corrections for weak lensing, photon-electron scattering, and gravitational energy redshifts along the line of sight. Because these steps depend on the assumed growth history of structure, the resulting spectra depend on the cosmological model adopted in the data reduction. We also note that the reported correlation between the hemispherical CMB asymmetry and the excess on-sky distribution of elliptical galaxies may point to matter inhomogeneities spanning large portions of the observable universe. Motivated by observational tensions in the hierarchical assembly of elliptical galaxies, we review the consequences of elliptical-galaxy formation on a downsizing timescale, which implies a significant extragalactic foreground released at 15< z < 20. Under conservative assumptions, this foreground contributes at least 1.4 percent of the observed CMB energy density and has not been included in CMB analyses. It follows that the reported 10^-5-level fluctuations is introduced by overly aggressive foreground cleaning.
http://arxiv.org/abs/2606.26098 (pdf) (27)
Compact and red sources in the high redshift (z∼5) Universe, known as "Little Red Dots" (LRDs), are among JWST&aposs most intriguing discoveries. These sources have broad Balmer emission lines, weak X-ray emission, and unique spectral energy distributions (SEDs) poorly fit by either stellar or AGN templates. Local LRD-like objects allow for detailed studies of the underlying physical processes with archival multi-wavelength datasets unavailable in the high-z Universe. We show that the SDSS ugriz filters at z≈0.4, 0.8 overlap well with the JWST filters used to select LRDs at z∼5. We use SDSS quasars to define a sample of ∼1300 Local Red Dots (LoRDs) which share the same photometric colors of LRDs. A subset of the LoRD sample selected to have V-shaped continua (N=244) show prominent higher-order Balmer absorption features and [NeV]λ3426 emission, both of which would likely be missed in JWST/PRISM observations given the low spectral resolution. A composite SED of the LoRDs differs from a typical quasar SED in the rest-frame UV/optical, but the two agree with each other in the NIR. The LoRD SED matches well with a stack of LRDs and can be modeled either as a reddened AGN combined with a host galaxy, or as a reddened AGN combined with a host galaxy and a cool blackbody. Interestingly, the LoRDs are X-ray detected at a rate comparable to typical quasars. However, the probability that LoRDs and typical quasars would go undetected, if subject to the LRD X-ray upper limits, is >50
http://arxiv.org/abs/2609.10700 (pdf) (26)
We investigate the internal structure of two galaxies from the LEWIS sample: the ultra-diffuse galaxy UDG-1 and the extended dwarf LSB-6. Both show coherent stellar rotation combined with a non-negligible fraction of random motion, with no signs of ongoing disturbance. This is the first attempt to constrain the dynamics and dark matter physics of rotation-supported UDGs through integral-field stellar kinematics. We model the galaxies as stellar spheroids with typical dwarf-like thickness embedded in spherical dark matter halos, deriving two-dimensional velocity fields that we compare with the observations in a Bayesian framework. Besides cold dark matter, we test fuzzy, self-interacting, and non-minimally coupled dark matter. The data prefer a cuspy halo in UDG-1 and a cuspless one in LSB-6. All the alternative models remain viable, although current data cannot conclusively discriminate among them. For fuzzy dark matter we obtain mutually consistent boson masses, m_α= 4.7^+1.8_-1.1× 10^-22 eV for UDG-1 and m_α= 6.3^+0.4_-0.4× 10^-22 eV for LSB-6. For self-interacting dark matter, LSB-6 yields a robust constraint on the velocity-weighted cross section, ⟨ σv⟩/m = 14.3^+2.0_-1.9 cm^2 km g^-1 s^-1. For non-minimally coupled dark matter we derive upper bounds on the coupling length, implying marginal deviations from ΛCDM. Regardless of the model, both galaxies exhibit a dark matter content comparable with halos of typical dwarf galaxies with similar stellar masses. Together with stellar population analyses and globular cluster information, these results support a scenario in which UDG-1 and LSB-6 originate from puffed-up dwarfs. This work is the first in a series exploiting kinematic data of extremely low-surface-brightness galaxies to test fundamental physics, from dark matter to modified theories of gravity.
http://arxiv.org/abs/2609.10680 (pdf) (26)
Mock catalogues for the Roman Galaxy Redshift Survey benefit from galaxy populations whose physical properties and emission-line observables are predicted within a single galaxy-formation framework. We present an emulator-assisted calibration of the semi-analytic model Galacticus for this purpose. Because direct posterior exploration with Galacticus is computationally prohibitive, we train Gaussian-process emulators to predict the observables entering the calibration likelihood. We calibrate the model simultaneously to stellar mass functions, the star-formation-rate function, size–mass relations, the black-hole mass–velocity-dispersion relation, the mass–metallicity relation, and H-alpha luminosity functions. Combining these observables breaks parameter degeneracies present in individual fits and can reveal tensions where different observables favour different regions of parameter space. The resulting model provides a broadly successful simultaneous description of the calibration data, including the abundance and redshift evolution of H-alpha emitters relevant to the Roman Galaxy Redshift Survey. A direct Galacticus calculation at the maximum a posteriori point gives similarly good agreement with the calibration data, providing a calibrated galaxy-formation model suitable for constructing physically consistent mock catalogues for Roman.
http://arxiv.org/abs/2609.10661 (pdf) (25)
We extend the Dark Matter-Baryon Separability Condition to show that the same framework developed for dark matter deficient galaxies naturally admits a conjugate branch describing preferential baryonic depletion. Using the recently discovered almost-dark galaxy TTT J1237327+143535 as a worked example, we derive a family of dynamical consistency relations parameterized by the unknown progenitor ratio μ_i, including thresholds for σ_los, M_dyn, the enclosed baryonic fraction and the dynamical mass to light ratio. We further connect the separability framework to the expected globular cluster population, providing an independent consistency test of the inferred halo mass. We also provide some bounds on the baryon ratio using the tidal properties of Virgo Cluster. Future measurements of stellar kinematics, gas content and globular clusters can therefore determine whether this galaxy occupies the positive separability branch and also test whether dark matter deficient and baryon depleted systems can be described within a common framework.
http://arxiv.org/abs/2609.10709 (pdf) (23)
Baryonic feedback is a dominant source of systematic uncertainty for upcoming weak-lensing surveys, but current tools for modeling its effect rely on spherical approximations and density profiles calibrated almost entirely on two-point statistics. We introduce BIND (Baryonic INpainting with Deep learning), a conditional flow-matching model that learns a field-level mapping from dark-matter-only halos to their hydrodynamical counterparts. BIND is trained on halos from the 1024 paired hydrodynamical and dark-matter-only simulations of the CAMELS 50 h^-1 Mpc SB35 suite and samples dark matter, gas, and stellar mass fields over redshift across the full 35-dimensional ΛCDM and IllustrisTNG galaxy formation parameter space. BIND recovers dark matter, gas, and stellar masses at the percent level, reproduces azimuthally averaged profiles to ≲10
http://arxiv.org/abs/2605.04781 (pdf) (23)
We search for a population-level signature of gravitational-wave recoiling supermassive black holes: a positive correlation between dust obscuration and the magnitude of the line-of-sight velocity offset of broad emission lines relative to the host. Using the SDSS DR16 quasar catalogue, we estimate the velocity offset, Δv, as the difference between the broad Hβ redshift and a noise-weighted redshift from narrow lines ([O III] 5007, [O II] 3728, and Ca II 3934). We adopt the redshift-relative colour excess Δ(g-i) as a proxy for dust column density. Analysing ∼10^5 quasars that meet basic spectral quality requirements, we find a modest but highly significant positive correlation between |Δv| and Δ(g-i) (Spearman r≃0.12 and Pearson r≃0.13, with p≪10^-10 in both cases). The fraction of highly obscured quasars increases with |Δv|, indicating that the correlation is driven by a dust-reddened subpopulation. The result is robust to the choice of minimum |Δv| threshold and to the line redshift estimator (peak vs. centroid). As expected, the correlation is largely absent when velocity offsets are computed between narrow emission lines. We find systematic differences between redshifted and blueshifted subsamples, which may point to residual velocity biases or additional physical effects (e.g. winds, inflows, orientation-dependent obscuration, or asymmetric broad-line regions). Recoiling massive black holes provide a natural explanation for the observed correlation, but alternative scenarios should be explored. If confirmed, this would enable population-level constraints on massive black hole merger rates, recoil dynamics, and active galactic nuclei disc properties.
http://arxiv.org/abs/2609.10665 (pdf) (19)
We present an emulator of the halo mass function (HMF) based on Gaussian-process regression, that exploits the evolution mapping framework. This framework splits cosmological parameters into shape parameters (_s), which set the shape of the linear power spectrum, and evolution parameters (_e), which only rescale its amplitude, with their net effect entirely captured by the clustering amplitude, σ_12. The emulation acts in two steps: first, an emulator predicts the overall shape of the HMF given _s and σ_12; then, a second emulator applies a small correction accounting for the residual dependence on the integrated growth history of density fluctuations. The emulator is trained on a suite of 100 simulations of flat ΛCDM cosmologies but, owing to evolution mapping, can be applied to predict the abundance of haloes in cosmologies with any extended space of _e, such as curvature, quintessence and dynamical dark energy. We test the emulator against independent ΛCDM runs and on an additional validation suite of dynamical dark energy simulations, achieving per cent-level accuracy over a wide range of masses and redshifts, with the performance only slightly degrading at early times, where the training data points are fewer. Compared to other state-of-the-art HMF emulators, its accuracy is higher in nearly all tested cases. Our results demonstrate the broad applicability of the evolution mapping framework to build very accurate emulators of statistics of the density field. The emulator is publicly available as part of the Aletheia emulator suite.
http://arxiv.org/abs/2609.11609 (pdf) (19)
Beginning in March 2026, the Event Horizon Telescope (EHT) conducted a two-month observing campaign aimed at obtaining the first time-resolved sequence of images of the black hole M87* at the center of the Virgo A galaxy. The primary scientific objectives of this campaign are to constrain the black hole spin and to investigate the temporal variability of the magnetic field and plasma properties within the horizon-scale jet-launching region. In this work, we present theoretical predictions for the spectral index maps expected from these observations. Our predictions are based on Magnetically Arrested Disk simulations performed with the EBHLIGHT code and post-processed using the IPOLE general relativistic ray-tracing code. We investigate the temporal evolution of the spectral index for a range of electron distribution models and five different black hole spin values. The analysis includes both frequency bands currently present within the EHT and neighboring frequencies proposed for future extensions of the EHT array. In addition to extending previous studies to a broader range of observing frequencies, our work systematically investigates the time variability of the predicted spectral index maps. In purely thermal electron models, the temporal variability of the spectral index maps closely traces variations in the magnetic field strength and electron temperature in the immediate vicinity of the event horizon. By contrast, models incorporating a non-thermal electron population exhibit substantially weaker spectral index variability, owing to the contribution of a power-law component with a fixed spectral slope. Our numerical results are consistent with theoretical expectations. Finally, a comparison with the recently measured integrated spectral index of the M87* enables us to place preliminary constraints on the model parameters.
http://arxiv.org/abs/2609.10710 (pdf) (19)
Recent multiwavelength observations of galaxy group and cluster gas suggest stronger baryonic feedback than our best-calibrated hydrodynamical simulations produce, while modeling this feedback remains a primary source of uncertainty in Stage-IV weak-lensing (WL) analyses. We present a suite of ray-traced maps generated with BIND (Baryonic INpainting with Deep learning), a conditional flow-matching model that paints baryonic mass and gas thermodynamics onto the halos of dark-matter-only simulations, and which was developed in a companion paper. Applied to IllustrisTNG300-Dark and ray-traced, we generate convergence, optical depth, and Compton-y maps at five source redshifts, each with 1000 pseudo-independent realizations. We build lightcones across a 256-node Sobol sequence spanning the thirty-dimensional IllustrisTNG galaxy formation prior, with individual parameter variations and at the fiducial model. In validation, the maps match those built from the IllustrisTNG300 halos to within LSST-Y10-like precision for a range of WL statistics. Across the prior, the response to feedback exceeds Stage-IV statistical precision by more than an order of magnitude on small scales, and different statistics respond to different model sectors: galactic winds control the WL power spectrum and gas auto- and cross-spectra, while the stellar initial mass function slope and AGN parameters shape the morphological statistics (PDF, peaks, minima, and Minkowski functionals). Finally, we find that the response of statistics can be compressed into seven halo properties which linearly predict a range of WL and SZ statistics. We publicly release the maps, statistics, and model tables.
http://arxiv.org/abs/2609.11479 (pdf) (16)
4U 0114+65 is a high-mass X-ray binary composed of the B1 Ia supergiant V* V662 Cas and one of the slowest known accreting neutron stars, with a spin period of ∼9.4 ks. In 2025, its long X-ray pulsations became undetectable in Swift/BAT monitoring, motivating a Director&aposs Discretionary Time observation with XMM-Newton. We compare this observation with a 2015 XMM-Newton observation, when the source was brighter and clearly pulsed, and analyze the long-term spin evolution using Swift/BAT data. We performed average and pulse-phase-resolved spectroscopy using the same model as in previous work. The 2025 observation still reveals weak pulsations, with a period of about 9.3 ks, despite their non-detection in Swift/BAT. The overall spectral shape remains similar in both epochs, but the luminosity decreased by about one order of magnitude, mainly due to strong suppression of the bulk-motion Comptonization component. Although the absorbing column is higher in 2025, the inferred wind properties remain broadly compatible with those from 2015, suggesting that no major global change in the donor wind is required. Instead, the results point to a substantial reduction in accretion efficiency close to the neutron-star magnetosphere. We propose that 4U 0114+65 may be evolving toward partial centrifugal inhibition in the strong-coupling regime, where the toroidal magnetic-field component is comparable to the poloidal one. If confirmed, this would represent the first observational evidence of this state. Accretion would become progressively less efficient and more intermittent without reaching a fully developed propeller regime. The apparent disappearance of the pulse in long-term hard X-ray monitoring would then result from the lower luminosity and reduced absolute pulsed flux, rather than from the loss of the underlying spin modulation.
http://arxiv.org/abs/2609.11487 (pdf) (15)
The C IV λλ1548,1551 resonance doublet is a key tracer of warm gas (T∼10^5 K) within and around galaxies. Recent observations have detected this line in both absorption and emission, revealing asymmetric profiles in galaxies and spatially extended haloes around active galactic nuclei (AGNs). Resonance scattering can strongly modify the emergent spectra and spatial distributions, complicating their interpretation. Using 3D Monte Carlo radiative transfer simulations, we study C IV resonance scattering over a broad range of column densities, intrinsic emission-line widths, and outflow velocities. We find that multiple scattering broadens the line profile and, in outflowing media, modifies the doublet ratio, R_CIV, defined as the flux ratio of the K and H components at 1548 and 1551 A, respectively. When the outflow velocity approaches or exceeds the doublet separation (≃500 km s^-1), K-line photons are redistributed around the H component, driving R_CIV below its intrinsic value and, in optically thick fast outflows, even below unity. We also combine photoionization models with resonance scattering to investigate extended C IV haloes around AGNs and compare them with He II λ1640 emission. Simple photoionization models do not produce C IV emission more extended than He II, whereas resonance scattering redistributes locally produced and central-source C IV photons to larger radii. These results demonstrate that the C IV doublet ratio and spatial distribution provide complementary diagnostics of warm gas.
http://arxiv.org/abs/2609.11313 (pdf) (15)
One-halo galactic conformity is the tendency for satellites in halos with quenched centrals to have lower star-formation activity than those in halos with star-forming centrals at fixed halo mass. It is an important probe of the galaxy–halo connection and halo-wide quenching processes that may couple central and satellite evolution. However, its existence remains controversial, because conformity must be measured at fixed halo mass, while halo masses are difficult to estimate accurately. In this Letter, we measure one-halo conformity in SDSS using five stellar-mass-complete samples and three halo-mass estimates: an ML estimate whose star-forming and quenched stellar mass–halo mass relations (SHMRs) agree with independent weak-lensing constraints, and two conventional abundance-matching (AM) estimates. We quantify conformity as the difference in median log( sSFR) between satellites of star-forming and quenched centrals, using both satellite-level and halo-level statistics. The two AM estimates produce strong positive conformity signals, consistent with previous AM-based measurements, but these signals are not reproduced with the ML halo masses. For the halo-level statistic, the representative AM-based signals are +0.38±0.04 dex and +0.23±0.04 dex for the luminosity-ranking and mass-ranking AM halo masses, detected relative to no conformity at about 10σ and 6σ, respectively. In contrast, the ML result is consistent with no conformity, +0.00±0.03 dex; the satellite-level statistic gives a similar result. Thus, with halo-mass estimates consistent with weak-lensing constraints, we find no detectable one-halo conformity signal in the present SDSS sample, suggesting that the strong AM-based signal is largely driven by halo-mass estimation biases.
http://arxiv.org/abs/2609.10662 (pdf) (14)
The recent 248-keV nuclear recoil candidate reported by LUX-ZEPLIN has renewed interest in endothermic inelastic dark matter, motivating us to examine its capture in primordial dark stars. Relative to the conventional elastic-capture picture in dark star evolution, endothermic capture adds two qualitative features. First, it opens only after the growing star crosses a compactness threshold, expressed as a kinematic radius R_kin_χM_⋆/δ set by the stellar mass, the dark-matter–nucleus reduced mass and the mass splitting. The accreting growth carries the star through the threshold, and the capture rate turns on quadratically above it. Second, although newly captured particles generically begin on nonthermal bound orbits, endothermic kinematics can keep this normally transient population spatially extended, turning it into a persistent reservoir rather than an intermediate step toward a thermal core. Following complete chains of state-changing collisions, we find that the reservoir compacts sharply and then stalls, because a ground state particle below a compactness-dependent orbital energy has no allowed up-scatter anywhere in the star. A percent-level radius contraction reopens the relaxation. These kinematic results do not depend on whether the excited state decays promptly or is long-lived. As a result, inelastic capture does not replenish a thermal annihilation core. The captured population forms an evolving orbital distribution that sets up the co-evolutionary dynamics between the star and the dark matter in the core and the reservoir, which we develop in a companion paper.
http://arxiv.org/abs/2609.11593 (pdf) (13)
The MICADO instrument for the ESO ELT will provide Adaptive Optics (AO) diffraction-limited NIR imaging and spectroscopy. A critical component of MICADO is the Point Spread Function Reconstruction (PSF-R) Software, which is designed to perform a blind reconstruction of the PSF for all observations using AO telemetry data only, both with the Single Conjugate (SCAO) and Multi-Conjugate Adaptive Optics (MCAO) modes of MICADO. The SCAO PSF-R Pipeline is currently being completed, and validation with ERIS at VLT and SOUL+LUCI at LBT has demonstrated high accuracy, of percent-level on-axis and 5-20% off-axis. Initial scientific simulations confirm that the PSF-R software meets the necessary precision for high-quality analysis, securing its role for MICADO&aposs first light and future ELT operations.
http://arxiv.org/abs/2609.10985 (pdf) (13)
Primordial black holes (PBHs) are a well-motivated dark matter candidate, and their cosmic abundance is constrained by a variety of observational probes. PBHs in the mass range 10^15 g - 10^17 g are evaporating today via Hawking radiation, a process that can heat interstellar dust and modify its thermal emission. Recent studies have used this effect to place constraints on the abundance of non-spinning PBHs. We extend this approach by investigating the influence of PBH spin on dust-heating constraints. Furthermore, we account for secondary photons that originate not only from the decay of gauge bosons but also from the decay of hadrons produced via the fragmentation of primary quarks and gluons emitted through Hawking radiation. By comparing the dust heating rate induced by spinning PBHs with the maximum cooling rate of dust, considering both silicate and graphite grains, we derive new upper limits on the fraction of dark matter in the form of PBHs, f_PBH. Our results show that the constraints depend on both PBH mass and spin. Smaller PBHs with higher spin yield stronger limits. For example, in the cases we investigated, the strongest constraint is f_PBH∼ 1.5 × 10^-4 for M_PBH = 10^15 g and spin parameter a_* = 0.9999. Although these limits are less stringent than existing constraints in the same mass range, they provide a distinct and complementary approach to constraining the abundance of PBHs.
http://arxiv.org/abs/2609.11181 (pdf) (12)
The aperture mass filter is a widely used tool in weak lensing studies. While previous surveys have typically relied on the flat-sky approximation, forthcoming Stage IV surveys such as Euclid will cover sufficiently large areas to require the construction of curved-sky aperture mass maps. In this paper, we extend the shear-based aperture mass formalism to the sphere and introduce a software package for a fast and precise computation of spherical aperture mass maps from the shear field, defined directly on the celestial sphere. Unlike existing approaches, our method does not rely on planar projections and does not require the reconstruction of convergence maps. The software provides four different implementations, all of which are independent of the choice of aperture mass filter function. Starting from the brute-force approach, whose computational cost renders it impractical for future surveys covering more than 10,000 deg2, we investigate alternative galaxy-grouping strategies to improve computational efficiency. We demonstrate speed-up factors of up to 40 at a resolution of NSIDE = 4096 while achieving significantly higher precision than existing state-of-the-art methods. In particular, we show that the spherical shear-based implementations using a pixel-based approach outperforms the existing methods in terms of both precision and computational time. Finally, we show that, for weak-lensing cluster detection, the cluster sample produced using the spherical shear-based aperture mass method consistently outperforms those obtained using projected shear-based approaches in the purity-completeness plane, while simplifying the resulting selection function. Our spherical mass aperture software, called SMART, is made available on GitHub.
http://arxiv.org/abs/2609.11764 (pdf) (12)
We derive the equations of relativistic magnetohydrodynamics from Lorentz-violating Myers-Pospelov electrodynamics. Starting from the modified fermionic and electromagnetic equations of motion, we employ the covariant Wigner formalism to establish the semiclassical spectral constraints and distribution function, which, combined with Noether&aposs theorem, yield the energy-momentum tensor, particle current, and corresponding conservation laws in the presence of Lorentz violation. The resulting hydrodynamic description contains Lorentz-violating corrections to the thermodynamic quantities and to the fluid currents. In the purely timelike sector, these corrections can be absorbed into effective thermodynamic quantities, preserving the isotropic structure of the fluid equations. In contrast, spacelike backgrounds introduce explicit contributions along the preferred direction, which cannot be absorbed into the standard hydrodynamic variables and give rise to anisotropic forces and torques. These effects lead to departures from standard relativistic magnetohydrodynamics and provide a direct macroscopic manifestation of the Lorentz-violating background. Our results provide a first-principles derivation of Lorentz-violating relativistic magnetohydrodynamics and establish a framework for investigating the macroscopic consequences of Lorentz violation in magnetized relativistic fluids.
http://arxiv.org/abs/2609.11555 (pdf) (12)
The scatter in elemental abundances among stars of similar age and metallicity reflects chemical inhomogeneity in their birth environments, making abundance scatter a key observable for chemical-tagging studies of Galactic formation and evolution. Using a large sample of subgiant stars with precise ages and elemental abundances derived from LAMOST low-resolution spectra, we investigate the intrinsic chemical abundance scatter of the low-α thin disk near the solar neighborhood (7 < R < 10 kpc). We model the abundance ratio [X/Fe], for each of the 19 elements of concern, as a function of age, [Fe/H], and [Mg/Fe], and deduce the intrinsic dispersions with a forward modelling technique. Our results confirm previous findings that the intrinsic scatters are small, typically ≲0.05 dex, for light elements (C, Al), α-elements (O, Mg, Si, Ca, Ti), and iron-peak elements (Mn, Ni). A dedicated analysis of M67 yields similarly small scatter values for these elements, implying limited discriminatory power from light-element abundances alone. In contrast, neutron-capture elements exhibit substantially larger scatters, typically ≳0.1 dex, which are significantly larger than those of M67 member stars (∼0.07 dex). In particular, our analysis suggests that the abundance variations of individual neutron-capture elements cannot be explained by a single tracer such as [Ba/Fe]. These findings clarify the utility of neutron-capture elements for chemical tagging and highlight the potential of low-resolution spectroscopy in such studies.
http://arxiv.org/abs/2609.10997 (pdf) (12)
SOLAR-C is an international solar-observing satellite mission led by Japan Aerospace Exploration Agency (JAXA). It aims to elucidate mass and energy transport in solar atmospheres through extreme ultraviolet (EUV) spectroscopy. The mission carries the EUV High-throughput Spectroscopic Telescope (EUVST) and the Solar Spectral Irradiance Monitor (SoSpIM), enabling comprehensive observations across a wide temperature range (10^4 K to 10^7 K) with minimum temperature gaps and high spatial and temporal resolution. To achieve its science objectives regarding atmospheric heating and solar flare eruptions, SOLAR-C implements a flexible and responsive operational procedure and a data processing system, while building on the heritage of the Hinode and IRIS satellites. The Chief Observer creates observation timelines that include core observation plans and approved proposed plans while taking into account the solar activity levels. A SpaceWire communication architecture is employed onboard the spacecraft, and a method is implemented in which the mission instrument temporarily acts as the network master during data transfer to support the high data rate requirements. Telemetry is downlinked at ground stations worldwide and gathered at the Institute of Space and Astronautical Science, JAXA. The EUVST data are calibrated at the SOLAR-C Science Center at Nagoya University, while the SoSpIM data are calibrated at the Processing and Archiving Facility before being integrated into the science data products. The data will be made publicly available immediately. The integrated operational scheme for this mission is expected to advance our understanding of solar atmospheric heating and flare processes.
http://arxiv.org/abs/2609.11798 (pdf) (11)
Binary stars are common in the Galaxy, and understanding how stellar binarity influences the formation and evolution of planetary systems is an active area of research. In this study, we use metal-enriched white dwarfs in wide binaries as tracers of long-lived planetary systems. With Data Release 1 from the Dark Energy Spectroscopic Instrument (DESI), we find that the fraction of cool metal-enriched white dwarfs in wide binaries is 9.8 ± 2.1
http://arxiv.org/abs/2609.10622 (pdf) (11)
In the second half of 2017, as Solar Cycle 24 approached its minimum, a pronounced and rigidity-dependent depression in the galactic cosmic ray (GCR) proton flux was observed using AMS-02. This event was driven by a sequence of strong coronal mass ejections (CMEs). In this work, we model the modulation of GCR protons during this period using a steady-state Parker transport framework with parameters constrained by AMS-02 observations. Treating the event as a purely long-term modulation interval yields excellent agreement with the observed proton intensities, but implies an inferred rigidity dependence of particle mean free paths that is inconsistent with basic turbulence theory, corresponding to an effective inversion of the scattering regime. We argue that this behaviour arises from the ability of highly flexible transport parameterizations to reproduce integrated modulation signatures through correlated adjustments of diffusion coefficients during transient events. To address this inconsistency, we relax the commonly adopted assumption of a time-independent ratio of perpendicular to parallel diffusion coefficients and introduce a physically motivated reduction of this ratio during the CME-dominated interval. This modification preserves the quality of the spectral fits whilst restoring the expected rigidity ordering of mean free paths. Our results indicate that allowing the perpendicular-to-parallel diffusion ratio to vary in time provides an effective means of constraining transport solutions during the 2017 depression and may help to resolve similar inconsistencies reported in other GCR modulation studies.
http://arxiv.org/abs/2609.11594 (pdf) (10)
MICADO, one of the ELT first-light instruments, will rely on point-spread-function reconstruction (PSF-R) to meet scientific requirements, especially for SCAO observations of targets that lack bright reference stars nearby. PSF-R will be performed in post-processing using only AO telemetry, recorded simultaneously with the science data. For SCAO observations, both on- and off-axis PSFs will be delivered. We present the current status and results from an ESO TTR campaign begun in 2023, targeting a sparse three-star asterism under varying conditions with the ERIS instrument at VLT. Reconstructed PSFs reach 10% accuracy in Strehl ratio and FWHM at separations up to half the size of the isoplanatic patch. The current results suggest the telemetry-driven approach can reproduce not only the core PSF structure but also the extended wings, which are essential for reliable photometry and morphological analysis of bright extended sources. The obtained results also corroborate the key role of wind as a limiting factor of the method, reducing the amount of recovered information on specific directions.
http://arxiv.org/abs/2609.10782 (pdf) (10)
The 3D distribution of dust and gas in the Milky Way&aposs Central Molecular Zone (CMZ) is key to understanding gas inflows toward the Galactic Centre (GC), the process of star formation in this extreme environment, and the propagation of energetic cosmic rays originating from Sgr A*. However, while recent efforts have combined datasets in a Bayesian framework to estimate the near/far positions of individual molecular clouds in the CMZ, conflicts between different methodologies still remain and we are still lacking a comprehensive, model-independent map of all of the gas and dust in the CMZ, which is critical to address key science questions. Here we develop a new methodology to infer the 3D dust distribution of the CMZ. The key idea of the method is to use stellar proper motions to get probabilistic information about the unknown stellar distances through a model of the distribution of star positions and velocities of the nuclear stellar disc (NSD), co-spatial to the CMZ. Taking stellar proper motions and extinctions as input, the latter adopted as a proxy of the dust column density, the method returns the 3D dust distribution. It is non parametric, makes no a-priori assumption on the dust distribution, and is fundamentally distinct and largely independent of all existing methods. We show that the method can robustly and effectively reconstruct the mock 3D CMZ structure by testing it on a range of mock dust distributions, both analytically generated and taken from hydrodynamical simulations. Finally, we discuss the prospects for applying the method to real data.
http://arxiv.org/abs/2609.10674 (pdf) (10)
An eV-scale QCD axion comprising the observed dark matter (DM) abundance is expected to generate a photon line at infrared energies that would be observable or near-observable in data collected by the James Webb Space Telescope (JWST), as might more general axion-like particles (ALPs) over a broader range of masses and couplings. This has motivated a number of efforts to either forecast JWST sensitivities to a QCD axion or realize them through analyses of publicly available datasets. At present, no consensus exists; leading analyses disagree by as much as an order of magnitude in terms of axion-coupling sensitivity, implying an orders-of-magnitude discrepancy in raw flux density sensitivity, and consistency between the analyses and prior forecasts is unclear. We address these outstanding discrepancies with a bespoke data reduction and flexible nonparametric inference procedure that lead to well-controlled and robust limits on the decay of eV-scale axion DM, consistent with previously forecasted sensitivities. We further demonstrate that the strongest previously claimed sensitivities exceed those attainable by any analysis of the datasets from which they were derived. We exclude QCD axion DM for masses between 500 meV and 2.5 eV using NIRSpec data, while setting limits on ALP DM complementary to other astrophysical constraints at masses between 100 meV and 500 meV. However, we find the sensitivities to be systematically limited, and therefore unlikely to be improved upon by ongoing data collection or re-analysis unless instrumental modeling and data reduction pipelines improve considerably.
http://arxiv.org/abs/2609.11787 (pdf) (9)
We present a broadband spectral study of the Sco-like Z source GX 349+2 using AstroSat LAXPC and SXT observations obtained during 2019 and 2024. The X-ray light curves exhibit large variability, and the hardness-intensity diagram traces the characteristic Z-shaped track. Using hardness ratio and intensity, the Z-track is segmented into six regions comprising normal branch, soft apex, and an extended flaring branch. We perform simultaneous broadband spectral fitting in the energy range 0.7-25.0 keV using two model configurations, tbabs*(thcomp*bbodyrad+diskbb) and tbabs*(bbodyrad+thcomp*diskbb), to study the evolution of emission from the neutron-star boundary layer, accretion disk, and Comptonizing corona along the Z-track. The comparative analysis reveals that, irrespective of the adopted Comptonization geometry, the inferred inner disk radius remains large ( 100 km) and exhibits no systematic inward motion as the source luminosity increases. This suggests that the accretion disk remains truncated and the disk luminosity tends to saturate, consistent with a radiation pressure influenced inner disk. The relatively small variation in the inferred mass accretion rate across the branches indicates that changes in Ṁ alone cannot account for the observed Z-track evolution. Instead, the boundary-layer component becomes progressively dominant toward the flaring branch, with both its temperature and flux increasing with total luminosity, while the relative contribution from the disk decreases. These results indicate a redistribution of accretion power toward the boundary layer and the associated Comptonizing region, suggesting that the primary energy dissipation in GX 349+2 occurs within the boundary layer region of the neutron star.
http://arxiv.org/abs/2609.11031 (pdf) (9)
We present optical photometric and spectroscopic observations of the Type II supernova SN2025abyc. During the optically thick phase between approximately 10 and 70 d after explosion, its light curves show strongly wavelength-dependent decline rates of approximately 2.7, 2.1, 0.9, and 0.8 mag/100d in the g, c, r, and o bands, respectively. At approximately 70 d, the light curves begin to depart from their nearly linear plateau evolution and gradually transition toward the radioactive tail. A Fermi-Dirac fit to the well-sampled ATLAS o-band light curve yields a transition midpoint of t_PT 100.5d. The interval between the end of the linear plateau and this transition midpoint is approximately 30 d, indicating a prolonged plateau-to-tail transition. This timescale is comparable to those measured for SN2013by, SN2013ej, and SN2014G. Spectroscopically, at +13 d post-explosion, the Halpha profile appears weak and broad, whereas Hbeta and Hgamma display clear P-Cygni profiles. This morphology can be explained by the normal early spectroscopic evolution of SNe II, although partial filling of the Halpha absorption trough by emission associated with circumstellar interaction cannot be excluded. SN2025abyc otherwise follows the general photospheric velocity evolution of SNe II, while remaining toward the high-velocity side of the comparison distribution in Halpha, Hbeta, and FeII. Exploratory light-curve modelling suggests a synthesized Ni mass of approximately 0.03-0.04 solar mass. We suggest that the extended circumstellar environment, Ni distribution, and hydrogen-envelope structure could all play a role in shaping the observed light-curve evolution, particularly the prolonged plateau-to-tail transition.
http://arxiv.org/abs/2609.10897 (pdf) (9)
The Vera C. Rubin Observatory&aposs Legacy Survey of Space and Time (LSST) is expected to significantly advance the study of Milky Way stellar streams. In particular, the deep, precise photometry from LSST should greatly increase the statistical sensitivity to density fluctuations in stellar streams, which can be used to probe the small-scale distribution of dark matter. However, current forecasts generally neglect the impact of observational systematics that will be imprinted on stream density measurements. In this study, we develop a realistic forward-modeling framework to inject stellar streams into LSST-like observations including photometric uncertainties, survey depth variations, background contamination, and imperfect star-galaxy classification. We develop a likelihood-ratio analysis to assess the detectability of gaps in stellar streams in the presence of these observational systematics. In the presence of realistic survey systematics, we find that after four years of operations, LSST will be sensitive to density reductions of ∼50
http://arxiv.org/abs/2609.10794 (pdf) (9)
We carry out a search for the relative baryon-cold dark matter (CDM) density and velocity modes in the DESI full-shape (FS) data, and investigate whether DESI&aposs preference for evolving dark energy can be affected by these relative perturbations, which are a known contaminant of the baryon acoustic oscillations. For the first time, we use a holistic effective field theory (EFT) treatment of the power spectrum and bispectrum both for the standard terms (without relative perturbations), modeled in EFT to one-loop order, and for the relative perturbations. For the latter we implement an efficient computation of the one-loop power spectrum that allows us to scan over different cosmologies. We obtain robust constraints on the relative density and velocity bias parameters b_δ_bc, and b_v^2_bc for all the DESI DR1 tracers, finding e.g. b_δ_bc=5 ± 6, b_v^2_bc=-0.031± 0.023 (at 68% CL) for the LRG2 sample in a DESI FS-only analysis within ΛCDM, supplemented by CMB priors on the baryon density and spectral tilt. Assuming instead a w_0w_aCDM cosmological model, and combining DESI FS with the CMB primary anisotropies, CMB lensing, and Pantheon+ supernovae, our LRG2 constraints read b_δ_bc=2 ± 5, b_v^2_bc=-0.032± 0.022. Adding the relative velocity divergence biases we find b_δ_bc=-17_-17^+19, b_θ_bc=-760_-680^+700, b_v^2_bc=-0.032± 0.024 for the same sample. The constraints on the dark energy equation of state parameters are virtually unchanged when marginalizing over the baryon-CDM modes. We conclude that DESI&aposs preference for evolving dark energy is robust with respect to the baryon-CDM relative perturbations.
http://arxiv.org/abs/2609.10705 (pdf) (9)
ESA&aposs astrometric space mission Gaia is uniquely positioned for microlensing studies thanks to its all-sky coverage and simultaneous astrometric and photometric observations. In Gaia Data Release 4 (DR4), astrometric time series will be published for the first time, providing a rich dataset for isolated stellar remnant searches. In anticipation of DR4, we prepared tools for fitting photometric and 1D astrometric measurements. We tested the recovery of true event parameters and population distributions with Gaia-like observations. We created a mock survey of microlensing events mirroring Gaia DR4&aposs design. The events were generated using dedicated Galactic simulations to provide realistic expectations on stellar and dark lens yields. For selected events, we simulated observations following Gaia&aposs scanning law and precision, using the astromet package. Such simulations are needed to interpret forthcoming data; they also enable testing inference tools with known ground truth. We used nested sampling to infer full posterior distributions and compare them with the true lens parameters. Our simulation predicts 322 ± 70 microlensing events with astrometric signals above Gaia precision thresholds; 209 ± 57 of them have detectable photometric signals. Among those, 82 ± 9 and 49 ± 7, respectively, have remnant lenses. To assess Gaia&aposs capabilities in constraining remnant mass distributions, we modelled 49 events from this sample. We recovered the true Einstein times and radii for most events. Microlensing parallax measurements are consistent with true values but uncertain, limiting mass determination. We recommend focusing on events with bright sources or anomalous fits to optimize dark lens searches. We make available the mock Gaia time-series dataset, as well as a toolkit repository for working with the data format, fitting models, and visualizing results.
http://arxiv.org/abs/2609.09419 (pdf) (9)
We study linear gravitational perturbations and the corresponding quasinormal-mode spectra of electrically and magnetically charged black holes in the presence of higher-order corrections to the Einstein-Hilbert action. In particular, we consider a four-derivative interaction that couples the Riemann tensor to the electromagnetic field strength in the unique (Horndeski) combination that preserves second-order equations of motion. We provide a comprehensive analysis of the stability conditions. We numerically confirm the absence of exponentially growing modes throughout the physically admissible region of parameter space where the theory remains free of pathologies. We then investigate in detail the properties of the quasinormal-mode spectra and identify several distinctive features induced by the higher-order interaction.
http://arxiv.org/abs/2609.11511 (pdf) (8)
The recent observation of a 220 PeV neutrino event (KM3-230213A) by the KM3NeT observatory presents a formidable challenge to standard astrophysical source models. We investigate the hypothesis that this ultra-high-energy signature originates from the terminal evaporation burst of a Primordial Black Hole (PBH). Since PBH evolution spans cosmic history, static vacuum approximations fail to capture early-universe dynamics. Embedding the PBH in a cosmological background via the McVittie spacetime, we demonstrate that early-universe cosmological accretion and expansion-suppressed Hawking emission shift the required initial mass window for a terminal burst occurring today. We show that although the early universe environment dictates the black hole&aposs overall lifespan, its final explosion today (z ≈ 0) is governed by standard Schwarzschild thermodynamics. This mechanism naturally produces the intense 220 PeV local flux while suppressing early emissions, thereby satisfying diffuse isotropic background limits. Consequently, this dynamical framework alters the mapping between current ultra-high-energy neutrino observables and the primordial curvature perturbations that seeded them.
http://arxiv.org/abs/2512.15870 (pdf) (8)
DESI DR2 data have been widely interpreted as evidence for late-time evolving dark energy (DE) with an apparent phantom crossing. Here we investigate an alternative explanation, based on early-Universe physics. If dark acoustic oscillations (DAO) are close in scale to baryon acoustic oscillations (BAO), they can bias the extraction of the BAO scale from the peak in the galaxy correlation function. This leads to an apparent shift in the inferred distance if the superposition of BAO and DAO features is misinterpreted as being due to BAO only. Taking this shift into account, we find that a DAO with percent-level amplitude can reconcile DESI DR2 with Planck 2018 as well as Pantheon+ supernovae data, with fit improvement at a similar level compared to evolving DE. Notably, a DAO feature with the required properties has been predicted in a previously proposed scenario that resolves the Hubble tension via a pre-recombination decoupling of dark matter and dark radiation (DRMD). The presence of a DAO feature close to the BAO peak can be scrutinized with future full-shape galaxy clustering data from DESI and Euclid.
http://arxiv.org/abs/2609.11306 (pdf) (7)
Initial spatial clustering of primordial black holes (PBHs) induced by local-type non-Gaussianity (LNG) can substantially modify cosmological constraints on PBH abundance. Existing formulations of initial PBH clustering are primarily restricted to either small non-Gaussianity or specific regimes, and a unified theoretical framework has been lacking. Here, we present a general analytical formulation for the PBH two-point correlation function applicable to arbitrary LNGs. Under the assumptions that PBHs form only at the large peaks of perturbations, and that long-wavelength modes weakly modulate the local variance of short-wavelength (PBH-scale) perturbations, we derive an analytic expression for the PBH bias parameter, directly connecting initial clustering to the primordial trispectrum in the collapsed limit. We demonstrate the versatility of our formula by computing the bias parameters in the ultra-slow-roll inflation, curvaton, and modulated reheating scenarios. We also formally generalize the framework to broad power spectra to account for correlations across different PBH mass scales. Because our formulation does not rely on weak or perturbative non-Gaussianity assumptions, our result provides a universal theoretical basis for evaluating initial clustering impact on PBH observables.
http://arxiv.org/abs/2609.11925 (pdf) (7)
WASP-12b&aposs orbit is decaying, for unknown reasons. The planet&aposs period is shrinking more rapidly than can be attributed to equilibrium tides or dynamical tides in a main-sequence star. Planetary obliquity tides could be sufficiently dissipative to drive WASP-12b&aposs inspiral, but would also damp the planet&aposs obliquity, halting the decay. Millholland Laughlin proposed that a nearby, low-mass planet (∼ 10 M_⊕) is maintaining a large obliquity for WASP-12b, sustaining the dissipation. We re-evaluated this hypothesis, finding that the companion must be more massive than originally proposed (≳ 65 M_⊕) to absorb WASP-12b&aposs orbital angular momentum. Radial velocity data allowed us to rule out a companion of this type. Any companions within 3 AU have K ≲ 14 m/s at 95% confidence.
http://arxiv.org/abs/2609.11908 (pdf) (7)
We formulate the theoretical aspects of direct waves (DWs) in the case of extreme-mass merger. A DW is a source-driven waveform characterized by a complex frequency ω_G, which reflects the orbital motion of the particle in the vicinity of the black hole, including a part of the orbit inside the ergoregion: its real part is governed by frame dragging and its imaginary part by the redshift of the source. Using the Green&aposs function technique, we derive the source-driven frequency ω_G, describe its screening by the potential barrier, and discuss its relation to dynamically excited quasinormal modes (QNMs). We also discuss the late-time decay of DW and predict that it decays with the third-order horizon mode. We then introduce a pole-splitting method, which divides the whole waveform into a QNM-pole sector and a non-QNM sector. Unlike QNM filtering, which multiplies the waveform spectrum by a filter function and thereby deforms it through a frequency-dependent time shift (i.e., group delay), our pole-splitting method merely divides the transfer function into pole and non-pole parts, separating the full waveform. Simulating a quasi-circular plunge into a Kerr black hole with medium and rapid spins, we find that the frequency and decay rate of the non-pole sector in the dominant mode, ℓ = m = 2, evolve consistently with ω_G-or with its screened counterpart ω_screen-establishing the DW as a probe of the ergoregion and of the redshift effect around a black hole.
http://arxiv.org/abs/2609.11907 (pdf) (7)
What if dark energy is recurrent throughout cosmic history? In this picture, episodes of scalar-field dark energy become less surprising and more natural features of cosmic evolution. We construct a homogeneous multi-axion cosmology with a transient contribution during Big Bang nucleosynthesis, two early dark energy components before recombination, and a thawing field that supplies the present dark-energy density. As Hubble friction weakens, the fields begin to roll at successive epochs set by their potential-curvature scales. Their initial displacements affect the rolling delays and peak energy fractions. The transient fields have third-power cosine potentials, whose sextic minima allow faster-than-radiation dilution during rapid, small-amplitude oscillations. All four fields are evolved in a common Friedmann background, with their first roll and subsequent dynamics resolved numerically. With reference matter and radiation densities taken from Planck 2018, the benchmark has a nucleosynthesis-era peak fraction of approximately 0.99
http://arxiv.org/abs/2609.11453 (pdf) (7)
Lunar pits, some of which are interpreted as collapse features into underlying lava tubes, expose otherwise inaccessible stratigraphy and may provide entry points to subsurface voids that preserve records of lunar volcanism and offer potential sites for future human exploration. We synthesize the current state of knowledge on lunar pits and lava tubes, covering their morphological characteristics, classification, proposed formation mechanisms, mechanical stability, and detection from orbit. We then review the open science questions that pit and pit-wall investigation is uniquely placed to address, spanning the volcanic stratigraphy of the lunar maria, the structure and lateral variability of the regolith, and the dimensions and accessibility of subsurface conduits. To evaluate how these questions can be tackled in situ, we assess the feasibility and expected performance of geophysical and remote-sensing investigations for subsurface voids and surface exposures, mainly focusing on gravity measurements, ground-penetrating radar, high-resolution imaging, and spectroscopy. Building on this, we present LunarLeaper, a small legged robot mission concept combining a gravimeter, ground-penetrating radar, high-resolution imager, spectrometer, and leg-based geomechanical experiments to deliver the first in situ investigation of a mare pit. The concept targets the Marius Hills Pit and its associated rille, with a mobility architecture optimized for the rugged terrain encountered at pit edges and funnel slopes.
http://arxiv.org/abs/2609.10955 (pdf) (7)
The populations of stellar-mass black holes (BHs) in globular clusters (GCs) play a key role in their dynamical evolution, however the mechanisms surrounding their formation and retention are uncertain. In this work, we extend the analysis of Paper I by fitting coupled rapid cluster evolution and multimass equilibrium models to a large sample of Milky Way GCs, under a variety of prescriptions for stellar evolution, BH formation and supernovae (SN) natal kicks. We explore the impacts of adopting SSE or PARSEC (through SEVN) prescriptions for BH initial-final mass relations, the rapid or delayed SN fallback mechanisms, and an ad hoc grid of kick strengths ejecting between 40 and 80 per cent of all BHs formed. All models reproduce the same present-day conditions despite starting from notably different initial BH populations, due to the correlation found between the initial cluster densities and initial BH mass fractions. A linear relationship is found between the (log) initial half-mass density and the initial BH mass fraction, with the SEVN models resulting in median densities (ρ_h,0∼ 10^7.2±1.1 M_⊙ pc^-3) nearly an order of magnitude higher than those of SSE (ρ_h,0∼ 10^6.4±0.9 M_⊙ pc^-3). We also find that both the bottom-light initial mass functions and the present-day BH mass fractions previously inferred are relatively robust against the stellar evolution models and natal kick prescriptions assumed. Finally, we discuss the implications of these results on the expected numbers and properties of dynamical binary-BH mergers, and the growth of intermediate-mass BHs.