Most interesting arXiv papers
Monday 19-10-2026
http://arxiv.org/abs/2610.22766 (pdf) (64)
Dust masses and grain properties in z>10 galaxies remain poorly constrained, as JWST UV attenuation and ALMA dust re-emission are rarely modelled consistently. We perform post-processed 3D dust radiative transfer on 208 snapshots of 99 galaxies at z=10-15 from the FirstLight simulations (M_*=10^6.5-10^9.5 M_⊙), applying twelve dust models per snapshot (2496 runs): four dust-to-metal (DTM) prescriptions and three grain size distributions. Metallicity-dependent DTM models reproduce ALMA dust-mass upper limits at z≳10 and detections at z∼7, as dust-to-stellar mass ratios rise sharply near Z∼0.2 Z_⊙, typically at M_*≳10^8.5 M_⊙, where grain growth by metal accretion becomes dominant. A constant DTM of 0.4 instead violates more than half of the upper limits. For dust attenuation, SMC-like dust is favored by the observed M_UV-β_UV relation. The observed UV slopes (β_UV≃-2.3) favor small grains (a_UV/2π≃0.025 μ m) producing steep UV extinction, together with a low abundance of carbonaceous grains or PAHs responsible for the 2175 A bump, consistent with silicate-rich dust. Large-grain stellar dust (a_0=0.1 μ m) instead yields nearly grey extinction and UV slopes close to the intrinsic value (β_UV,0≃-2.6). SMC and MW attenuation curves flatten at A_V≳0.3, showing that relatively flat attenuation does not necessarily imply large grains. For dust re-emission, dust peak temperatures reach 60-80 K, exceeding the 40-50 K commonly assumed in ALMA dust-mass estimates. Predicted IR continuum flux increases with UV luminosity; at current ALMA sensitivities (10-20 μJy), detections generally require M_UV≲-21.5. Roman wide-field surveys may identify such rare UV-bright galaxies.
http://arxiv.org/abs/2610.23016 (pdf) (48)
Assessing the significance of extreme galaxy candidates has typically relied on redshifts and stellar / halo masses inferred through inverse modelling, which carries substantial systematic uncertainty. We combine the halo mass function with a conditional model for galaxy emission and apply the Extreme Value Statistics (EVS) formalism to predict the distribution of the observed flux of the brightest galaxy expected in a given survey area, redshift interval and photometric band. We find the choice of parametric star formation history (SFH) dominates the modelling uncertainty, introducing a spread approaching ∼ 1 dex by z ∼ 18, whilst the stellar population synthesis model, initial mass function and metallicity each contribute ≲ 0.3 dex. The luminous z ≳ 9 candidates in COSMOS-Web cluster about the median prediction for the area they were drawn from, and lie at least an order of magnitude below the 3σ limit, a margin comparable to the SFH systematic. Of the observational samples reaching beyond z ∼ 13, the two spectroscopically confirmed galaxies at z ≃ 14 lie close to the median prediction for their area, whereas photometric candidates at z ≳ 16 reach +4σ (a probability to exceed of ∼ 10^-5), the brightest lying at the limit obtained by converting all halo baryons into stars. Since the prediction bounds the flux far more tightly than the photometry bounds the redshift, such offsets are most naturally resolved by errors in the inferred redshifts. This approach allows extreme candidates to be assessed directly in the observational plane, making explicit the modelling assumptions common to forward and inverse approaches.
http://arxiv.org/abs/2610.23567 (pdf) (40)
Radio galaxies play a significant role in galaxy evolution through feedback mechanisms, making it important to characterize the environments in which they reside. We measure the cross-correlation between the positions of radio sources in the MIGHTEE survey and cosmic microwave background lensing measurements from the Atacama Cosmology Telescope over 20 deg^2 to flux densities S_1.4 GHz>40 μJy, detecting the signal at ∼4σ significance. By combining this with auto-correlation measurements of low- and high-redshift radio samples, we constrain the redshift evolution of the galaxy bias parameter, b(z). Fitting different models for the redshift evolution of the bias, we find that both a model inversely proportional to the linear growth factor, b(z) = (1.63^+0.19_-0.22)/D(z), and a model where bias grows linearly with redshift, b(z) = (0.49^+0.44_-0.32) + (1.62^+0.37_-0.51)z, provide comparably good fits to the data, with the current measurements unable to distinguish between them. Using a radio luminosity function, we estimate the fraction of sources that are active galactic nuclei (AGNs) or star-forming galaxies (SFGs) as a function of redshift, and model the individual contributions of AGNs and SFGs to the total bias. This approach allows us to measure the bias of AGNs and SFGs without individual source classifications. Though no single model is favoured significantly, the best-fitting model gives b_AGN(z) = (2.60^+0.65_-0.91)/D(z) compared to b_SFG(z) = (0.52^+0.47_-0.34)/D(z). Incorporating SFG bias measurements from previous work as a prior tightens these constraints to b_AGN(z) = (1.96^+0.37_-0.40)/D(z) and b_SFG(z) = (0.90^+0.06_-0.06)/D(z).
http://arxiv.org/abs/2610.22860 (pdf) (36)
Recent evidence points to an intermediate-mass black hole (IMBH) in the most massive globular cluster of the Milky Way, ω Centauri. This globular cluster is commonly associated with the Gaia-Sausage/Enceladus (GSE) dwarf galaxy, the last major merger event of our Galaxy, as its nuclear star cluster remnant. Motivated by this, we investigate whether the interaction between a bright GSE-like dwarf galaxy with stellar mass in the range 10^8.5–10^9.7 solar masses and a Milky Way-like host can either suppress or enhance the growth of a central black hole (BH) during the merger process. To this end, we perform high-resolution N-body and smoothed particle hydrodynamics (SPH) simulations with the ChaNGa framework and a modified, angular-momentum-aware Bondi-Hoyle accretion prescription, considering both isolated dwarf galaxies and merger simulations with different stellar feedback efficiencies. The merger simulations show that the interaction does not halt BH accretion; instead, tidal interactions and gas inflows during pericentric passages strongly enhance accretion, leading to BH mass growth by more than two orders of magnitude in both low- and high-feedback models. The enhanced accretion is accompanied by suppressed star formation in the dwarf galaxy. Our results suggest that mergers similar to the GSE event can efficiently fuel IMBH growth in satellite galaxies rather than suppress it. This may be in tension with the estimated low mass of the proposed IMBH in ω Centauri, but it could provide a mechanism for triggering active galactic nucleus activity in dwarf galaxies.
http://arxiv.org/abs/2610.23200 (pdf) (29)
This study combines multi-wavelength continuum data with optical spectral properties for a far-infrared (FIR)-selected sample of Type 1 active galactic nuclei (AGNs) from AKARI data. We explore the relation between AGN physics and concomitant star formation processes using empirical correlations of optical and multiwavelength data. We measured spectral parameters including full widths at half maximum (FWHMs) and velocity shifts of broad and narrow optical emission lines, as well as continuum luminosities from GALEX FUV to AKARI 160 microns. The sample contains 88 AGNs observed at both optical and FIR wavelengths. About 50% are Seyfert 1s with both broad Hβ and Hα lines; almost all remaining objects are S1.9s with only broad Hα. The FWHMs of broad Hβ and Hα are anti-correlated with luminosity in most bands. Shifts of narrow Hβ, Hα, and [O III] 5007 lines and FWHM of [O III] 5007 correlate with optical and IR luminosities. Luminosities at 65 and 90 microns are more strongly connected with optical and other spectral properties than those at 140 and 160 microns. These trends extend the well-known Eigenvector 1 (EV1) correlations into the MIR and FIR domain. The anti-correlation between broad-line FWHMs and NIR and MIR luminosities is probably caused by the inclination of the torus and accretion disc: higher inclinations see less AGN radiation. However, most objects fall into Population A along the quasar EV1 main sequence, where narrower broad Hβ profiles and FIR excess emission are commonly linked to higher Eddington ratios.
http://arxiv.org/abs/2610.23349 (pdf) (28)
A key method of exploring the relationship between active galactic nuclei (AGN) and their host galaxies is to determine the astrophysical properties of galaxies that host AGN, and compare them to those that do not. Star formation rates are believed to be strongly related to nuclear activity, both indirectly via their common fuel supply, and directly via the energy injected by the AGN either suppressing or enhancing star formation. Strong observational evidence of a direct link between nuclear activity and in-situ star formation has proven elusive, likely in part due to the very different timescales of nuclear activity and episodes of star formation. To overcome these challenges, some studies have demonstrated that evidence of the aforementioned link may be found in the distributions of AGN power and star formation rates. While previous attempts to determine these distributions have largely focused on X-ray-selected samples of AGN, in this study we also attempt to determine the star formation rate distributions of infrared-selected AGN. The results vary between our samples, but at low redshift (0.8 < z < 1.2) we find that more powerful X-ray-selected AGN are more likely to occur in galaxies with higher specific star formation rates, consistent with previous studies. At higher redshifts (1.0 < z < 1.8), this trend is reversed with higher power AGN tending to occur in galaxies with lower specific star formation rates. By contrast, infrared-selected galaxies at 1.0 < z < 1.8 show no evidence of a changing star formation rate distribution with AGN power.
http://arxiv.org/abs/2610.22547 (pdf) (26)
We perform a joint Roman–Rubin photometry analysis using the LSST Science Pipelines (LSP) on simulated images, by using the Roman F158 band as the reference band and applying forced photometry on the other bands. We find that the LSP CModel photometry for Roman bands is on average unbiased, but with a bright-end systematic uncertainty floor of 7-8%. The pgauss photometry as part of the Metadetection weak-lensing shape measurement pipeline also exhibits a similar floor despite a substantially different flux-bias behavior. However, this systematic uncertainty floor becomes negligible for the color residuals and both algorithms recover unbiased colors for Roman as well as Rubin. We also show that not accounting for correlated pixel noise in Roman IMCOM coadded images results in underestimation of flux errors by a factor of 4. Based on the measured photometry, we build a machine-learning algorithm that maps the true galaxy properties to empirical flux errors, providing a convenient method to reproduce realistic photometry errors without having to rerun expensive image simulations and photometry pipelines. We apply this empirical error model to a new mock galaxy catalog with improved realism to investigate the impact of different survey depths and configurations on photometric redshift (photo-z) estimation performance. We find that 1) Rubin survey depth significantly impacts the photo-z precision, 2) additional systematic uncertainties on photometry could significantly degrade the photo-z performance especially towards high redshifts, therefore eroding the Roman–Rubin synergy, and 3) Roman photometry becomes important above z 1.5 while adding Rubin photometry to that of Roman improves photo-z performance across all redshifts. Our results provide a foundation for understanding how realistic uncertainties in galaxy photometry may propagate into the Roman–Rubin synergy in photo-z analyses.
http://arxiv.org/abs/2606.02725 (pdf) (25)
After four observational runs, the Ligo-Virgo-Kagra collaboration estimated a local binary black hole (BBH) merger rate density of R_0,LVK≃ 14-26 Gpc^-3 yr^-1 within the 90% credible interval. Some previous studies already pointed out that, when a realistic evolution of the metallicity-dependent cosmic star formation rate density (SFRD) is adopted, theoretical models predict a local BBH merger rate density that exceeds the observed value by at least a factor of ∼ 10. In this paper, we confirm and strengthen this claim by constructing an empirical model for the SFRD and metallicity evolution that includes a correction accounting for iron abundance. The adopted metallicity relation is flexible, enabling us to bracket the wide range of observational uncertainties. We show that, even under the most conservative assumptions regarding both the SFRD and the metallicity relation, the local BBH merger rate density is overestimated by a factor > 10. Attempts to reconcile the predicted and observed merger rates by modifying only the metallicity-dependent SFRD would require unrealistically high metallicities (Z>Z_⊙) even in low-mass galaxies at high redshift. This finding indicates that revisions to the treatment of stellar and binary evolution are necessary to achieve consistency between theoretical predictions and observations. We suggest that even a modest steepening of the delay-time distribution could help alleviate this tension.
http://arxiv.org/abs/2610.23333 (pdf) (24)
Localized fast radio bursts (FRBs) show substantial scatter in dispersion measure (DM) about the mean cosmic DM–redshift relation, also known as the Macquart relation. Both foreground structure and host galaxies contribute to this scatter. In this work, we examine foreground tracers and host galaxy properties using a parent catalog of 136 localized FRBs with measured DMs and spectroscopic host redshifts. Among 77 FRBs within the Gfinder footprint, DM excess relative to a fixed Macquart mean increases with foreground group count excess, with a Spearman coefficient of 0.459. Adding the foreground term to a mean host baseline reduces the leave-one-out root mean square error (RMSE) and mean absolute error (MAE) in extragalactic DM by 13.0
http://arxiv.org/abs/2610.23327 (pdf) (24)
Galaxies resembling the Milky Way in present-day global properties may have followed different evolutionary histories. The Milky Way Near Twins (MWNeT) framework introduced in Paper I addresses this through seven successive stages of present-day compatibility and evolutionary diagnostics. We apply this framework to NGC 3521, NGC 2903, UGC 11487, and CGCG 047-026, spanning different levels of observational completeness and outcomes within the hierarchy. We combine archival measurements with analyses of galaxy structure, spectral-energy-distribution similarity, atomic-gas profiles, and optical nuclear diagnostics. Stages I–IV assess environmental, structural, nuclear, and global compatibility; Stages V–VII examine SEDs, rotation curves, and evolutionary fingerprints. We distinguish supported compatibility, partial agreement, measured incompatibility, and unconstrained diagnostics. We jointly assess these four galaxies within the MWNeT framework for the first time and derive new structural, SED-comparison, and H i profile measurements for UGC 11487. NGC 2903 is a well-constrained benchmark, whereas NGC 3521 combines extensive evolutionary information with uncertain bar identification. UGC 11487 has a securely detected bar, quantitative UV–IR SED similarity, and distinctive transient nuclear tracers, but its nominal stellar mass lies below the adopted interval and its extended rotation curve remains unconstrained. CGCG 047-026 meets the initial environmental, structural, and nuclear-state criteria but remains largely unconstrained at higher stages. The four applications demonstrate the framework&aposs ability to organise complementary evidence and identify observational priorities without equating missing data with physical incompatibility. They establish a basis for future systematic searches and for using candidate near twins as external laboratories of Milky Way evolution.
http://arxiv.org/abs/2610.22483 (pdf) (23)
We present a multi-wavelength analysis of the nearby spiral galaxy NGC 6946, combining optical IFS, infrared, and radio observations. We aim to characterise both the integrated and spatially resolved properties of the stellar populations and ISM, and to interpret them in the context of bar-driven galaxy evolution. We used the SED fitting code SINOPSIS to derive stellar population properties, and different prescriptions to derive ISM properties including atomic and molecular gas masses and the dust mass,temperature, and emissivity index. We find that the central 2.6 kpc hosts 40
http://arxiv.org/abs/2610.22449 (pdf) (23)
Understanding how supermassive black holes (SMBHs) grow and regulate the growth of their host galaxies requires tracing the cold molecular interstellar medium (ISM) throughout the full AGN feeding-feedback cycle. Molecular gas constitutes both the primary fuel for star formation and black-hole accretion and one of the ISM phases most directly affected by AGN-driven radiation, winds, and radio jets. Over the past decade, ALMA has resolved circumnuclear molecular discs, rings, spirals, tori, and jet-driven outflows on scales of tens of parsecs in nearby AGN. These studies have demonstrated that combining spatially resolved gas kinematics with multi-line diagnostics is essential for disentangling the physical processes underlying the observed gas conditions. However, our current picture remains largely based on a limited number of exceptionally nearby systems. Most local AGN have instead been observed with ALMA in a single transition at moderate or low angular resolution, preventing a systematic characterisation of both the dynamics and physical state of the gas. Moreover, existing studies have largely overlooked kinetic-mode AGN hosted by massive early-type galaxies, despite their central importance for understanding maintenance-mode feedback. In this White Paper, we argue that ALMA2040 should transform these detailed case studies into a statistical framework. A ∼10× improvement in line sensitivity relative to the current ALMA, together with fast mosaicking, ultra-wide instantaneous bandwidth, and longer baselines (resolving up to ∼9 mas at 230 GHz), will enable homogeneous, multi-line mapping of the molecular ISM on giant molecular clouds scales across statistically significant samples of both radiative- and kinetic-mode AGN out to z≲0.1.
http://arxiv.org/abs/2610.22598 (pdf) (22)
The interstellar medium (ISM) is magnetized at all scales and phases in our Galaxy and other galaxies across the Universe. Spectro-polarimetric millimetre observations are essential to probe the role of the interstellar magnetic fields in the dynamics of the cold universe, and how this affects the star formation process, and planet formation in nearby star forming regions. Such observations are also essential for understanding magnetic fields in active stars and the shaping and evolution of the circumstellar environments of evolved stars. While ALMA has been revolutionary in advancing our knowledge about magnetic fields in the ISM, the limitations in calibration, especially for circular polarization, and sensitivity for large samples, have hindered robust measurements of field strengths, which are essential for fully assessing their role in ISM dynamics. An ALMA2040 that broadens spectro-polarimetric capabilities to include robust circular calibration and increases sensitivity with more collecting area will enable crucial observations of magnetic fields across the Universe and over a wide range of spatial scales to fully ascertain the influences of magnetic fields from galaxy assembly to forming planets.
http://arxiv.org/abs/2610.22481 (pdf) (22)
Determining whether a satellite galaxy is tidally disrupting from photometric observations alone remains a fundamental challenge in near-field cosmology. Cosmological simulations of Milky Way-mass galaxies predict that satellite populations should be highly disrupted, yet the observed Milky Way satellite population appears largely intact. Understanding this tension requires either uncovering disruption that has gone undetected in Milky Way satellites, or identifying limitations in the simulations themselves. We use the Auriga suite of cosmological simulations to evaluate several commonly used indicators of disruption. We find that, in the simulated satellite population, the shape, size, and proper motion alignment of a satellite&aposs bound component are unreliable indicators of its present-day disruption state despite their widespread use in the literature. In contrast, the extended stellar distribution, specifically the ratio of stellar mass between 3r_1/2 and 10r_1/2 to that within r_1/2, is strongly correlated with disruption. We introduce the extended mass ratio μ_EMR as a metric for comparing extended stellar distributions across satellites of varying masses and exploring their connection to disruption in Milky Way observations. Applying this metric to a sample of Milky Way satellites with measured extended stellar distributions, we find that many are consistent with significant ongoing disruption, even among systems traditionally classified as intact. These results suggest that tidal disruption may be substantially more widespread in the observed Milky Way satellite population than currently recognized, and establish extended stellar distributions as a promising new tool for uncovering this hidden disruption.
http://arxiv.org/abs/2610.23340 (pdf) (20)
A delay time distribution (DTD) quantifies the lag between progenitor formation and the subsequent explosive event. The DTD for a sample of Type Ia supernovae (SNe Ia) from a survey can be modelled through relative rates and the star formation histories (SFH) of all individual galaxies. Modern wide-field surveys require computationally efficient methods to model the DTD without modelling large samples of non-host galaxies. We model the DTD for a volume limited sample from SNe Ia from the Zwicky Transient Facility. The aim is to study the effect of varying SED modelling assumptions, including stellar templates, dust attenuation, SFH models, and photometric coverage. We present a new modelling framework, the Conditional Astral Ancestry Latency Analysis Machine (CAALAM), to infer the DTD of a sample of SNe Ia from the posterior of their host galaxy SFH. A model consisting of a standard stellar template, a non-parametric SFH and a flexible dust attenuation law, is favoured by the data and predicts a non-parametric DTD peaking at 100-330 Myr. When modelled as a power law, the DTD has a characteristic exponent of -1.30 ± 0.04. The exponent is sensitive to dust attenuation, photometry, or SFH. When the DTD is convolved with the SFHs of the host galaxies, we find a bimodal progenitor age distribution, with ∼60
http://arxiv.org/abs/2610.23806 (pdf) (20)
We use the observed abundance and spatial clustering of JWST little red dots (LRDs) over z = 4–8 to confront two formation scenarios: the standard picture in which LRDs form in virialized dark-matter halos, and one in which they are seeded by a dynamically subdominant population of primordial black holes (PBHs). The halo formation scenario accounts for both observables naturally and without fine-tuning. In the PBH scenario, matching the observed LRD abundance requires a negligibly small PBH mass fraction, and a power-spectrum analysis shows that such a Poisson isocurvature component predicts a clustering signal far below what is observed. The same conclusion holds for any population of dynamically subdominant early-forming seeds, unless their initial clustering is finely tuned or their activity is tied to massive host halos. We conclude that the spatial clustering of JWST LRDs strongly favours the standard halo formation picture and disfavours PBHs, or similar early-forming seeds, as a randomly activated LRD population; PBHs remain viable as seeds only if their activity is governed by massive host halos, in which case clustering alone cannot distinguish them from astrophysical seeds.
http://arxiv.org/abs/2610.23339 (pdf) (20)
Short-timescale X-ray variability provides a direct probe of unstable accretion onto compact objects, yet extragalactic X-ray binaries (XRBs) with recurrent flaring remain largely unexplored. Here we present a detailed X-ray study of M81F4 (2CXO J095524.2+690957), a persistently luminous X-ray source towards M81 that exhibits recurrent short X-ray flares. Based on multi-mission X-ray observations spanning two decades, it is found that M81F4 remains persistently bright at L_0.5-8.0 keV∼10^38-10^39 erg s^-1, while producing flares lasting a few hundred seconds and reaching peak luminosities of a few 10^39 erg s^-1, with an observed mean (minimum) recurrence time of ∼0.5 day (∼6 ks). The stacked flare profile shows a rapid rise followed by a slower decay, but no significant hardness evolution. Spectral analysis reveals that both the flare and non-flare states are dominated by a hard power-law continuum, with evidence of highly ionized Fe lines in the non-flare state spectrum. These properties disfavor classical Type I and Type II bursts and have no close analog among known XRBs. We favor a magnetized neutron star accreting from an intermediate mass donor at a persistently high rate with the flares arising from instabilities at the disk-magnetosphere boundary. If confirmed, M81F4 might represent a hidden population of XRBs in which sustained, magnetically-regulated accretion can lead to brief episodes of super-Eddington emission.
http://arxiv.org/abs/2610.22462 (pdf) (19)
How does a star&aposs birth environment set the mass, angular momentum and chemical complexity delivered to its planet-forming disk, and thereby whether it builds planets, brown dwarfs or stellar companions? This requires following material from the cloud to the sub-au to few-au scales where infalling gas, magnetized winds and the disk interact, across all star-forming environments. At frequencies where young disks are optically thin, current facilities, even ALMA after its Wideband Sensitivity Upgrade, fall short by a factor of a few in angular resolution and an order of magnitude in sensitivity. Likewise, resolving how disks around distant high-mass stars fragment into tight multiples requires a linear resolution of a few to a few tens of au, currently out of reach. Yet the disks characterized so far lie almost exclusively in nearby, isolated clouds at near-solar metallicity, although the environment (radiation field, ionization rate, gas and stellar density, pressure and metallicity) varies by orders of magnitude across the Galaxy. Most stars, probably including the Sun, form in clustered regions at kpc distances, and the widest metallicity range lies farther still, in the outer Galaxy, the Central Molecular Zone and the Magellanic Clouds. With three times finer angular resolution and ten times the sensitivity, ALMA2040 will follow the path of material from cloud to disk across the full stellar-mass spectrum and, for the first time, measure how this path changes with environment: resolving the sub-au infall-impact and wind-launching zones nearby, disk fragmentation and the onset of multiplicity in distant high-mass regions, and disk formation at a few tenths of solar metallicity. It will thereby reveal the dust-obscured, embedded origins of planetary systems that optical and infrared facilities, including the ELT, see only once they emerge, and show how the birth environment shapes them.
http://arxiv.org/abs/2610.23636 (pdf) (18)
We present a comprehensive spectroscopic census of the Orion Nebula Cluster (ONC), based primarily on low-resolution optical spectroscopy obtained with MMT/Hectospec. We derive spectral types, youth indicators, and accretion diagnostics for 1513 stars located in the ONC region, enabling the identification and characterization of pre-main-sequence stars. We combine these measurements with a kinematic membership analysis based on Gaia Data Release 3 astrometry to identify and confirm cluster members. Using these combined diagnostics, we identify and characterize 780 bona fide T Tauri stars associated with the ONC. In addition, we construct a census of protoplanetary discs using infrared photometry obtained with the IRAC and MIPS instruments aboard the Spitzer Space Telescope. We further derive stellar masses, ages, and extinctions using the MassAge code, which combines photometric and spectroscopic inputs with stellar evolutionary models, and employs a Monte Carlo approach for uncertainty propagation. The resulting homogeneous catalogue of ONC members provides spectroscopic characterizations, kinematic membership constraints, disc diagnostics, and self-consistent estimates of stellar physical parameters. This dataset offers a refined view of the stellar population of the ONC and provides a robust foundation for investigating its kinematic structure, age distribution, and protoplanetary disc population.
http://arxiv.org/abs/2610.22406 (pdf) (18)
The Wolf-Rayet (WR) phase is a critical stage in the evolution of massive stars. During this phase, the stars experience a large mass loss in the form of fast winds that are enriched with heavy elements. There are only a few WR binaries in the galaxy, with Apep being a truly exceptional system exhibiting a bright colliding-wind region at radio wavelengths and an impressive dust plume in the infrared. Although the angular separation between the two (WN and WC) components has been indirectly inferred to be about 47 mas, it has yet to be confirmed through direct measurements. We aim to image the Apep system at millimetre wavelengths to determine the position of the WR stars and their mass-loss rates. We observed Apep with ALMA in Bands 6 and 7 in a large baseline configuration, achieving an angular resolution of approximately 20 mas. We used a canonical spherical wind model to convert the measured flux densities into wind mass-loss rates. Towards Apep, we detect point-like emission coming from the WN star and elongated emission most likely from the colliding-wind region; the WC star is not detected. We also detect the O star at 0.7” that has been associated with the system. Our analysis favours mass-loss rates for the WN and WC components of Ṁ_WN∼ 2×10^-5 M_⊙ yr^-1 and Ṁ_WC∼10^-5 M_⊙ yr^-1, respectively. The mass-loss rates of the WR stars are significantly lower than previously reported, and we examine several factors that may account for this discrepancy. The diffuse emission from the colliding-wind region cannot be accounted for by the synchrotron flux, and thus additional free-free or dust emission is present.
http://arxiv.org/abs/2610.22707 (pdf) (18)
Unresolved binary stars introduce systematic biases in stellar parameter estimates when analysed assuming a single star. Using ∼917,000 spectra from the Galactic Archaeology with HERMES (GALAH) Data Release 4 (DR4), we develop a binary spectral model that fits each observed spectrum as the superposition of two synthetic spectra generated from the same model architecture used by GALAH DR4. These generative spectral models take as input: effective temperature (T_eff), surface gravity (log g), iron abundance ([Fe/H]), projected rotation velocity (v sin i), microturbulence velocity (v_mic), and radial velocity (v_rad). In our binary model, T_eff and log g are not fitted directly but are instead interpolated from PARSEC isochrones based on stellar mass and age, under the assumption that both components formed contemporaneously with the same compositions. We fit for the primary and secondary masses, a common age and [Fe/H], as well as individual v sin i, v_mic, and v_rad for each component. After iterative optimization and conservative quality cuts to ensure robust detections, we identify ∼12,911 candidate binary systems. We quantify typical biases introduced by unresolved binarity: T_eff shifts of -142 to +200 K, log g shifts of -0.070 to +0.35 dex, and [Fe/H] shifts of -0.055 to +0.13 dex. Despite our conservative approach reducing the number of detections, we recover known trends of increasing binary fraction with higher primary mass and lower metallicity, thus reflecting the broad trends of binary populations. Accounting for unresolved binaries is essential for accurate stellar population modelling and for interpreting spectroscopic surveys–particularly in low-metallicity systems or populations with more massive stars.
http://arxiv.org/abs/2610.23142 (pdf) (17)
Red supergiants (RSGs) provide a near-infrared metallicity diagnostic complementary to optical hot-supergiant spectroscopy, but their low-resolution extragalactic use requires tests of parameter identifiability and modelling robustness. We quantify the information content of extended J-band RSG spectra at R=2000 over 1.155-1.340 um and assess recovery of effective temperature, surface gravity, microturbulence, and metallicity. Using MARCS/TURBOSPECTRUM spectra and an emulator, we perform a six-node Fisher analysis and controlled mock experiments, and analyse 21 SpeX spectra of 20 RSG candidates towards Per OB1. All four parameters produce measurable responses, but only 8-10% of the raw metallicity information and 4-7% of the gravity information are independent of the other parameter directions. Extending beyond 1.22 um provides median fractions of 60% of the raw and 66% of the nuisance-orthogonal metallicity information. Gravity dominates the weakest Fisher mode, while its prior choice shifts association means by at most 0.03 dex. Matched mocks recover metallicity with an rms error of 0.04 dex. Standard continuum treatment biases solar-metallicity mocks by about -0.23 dex, reduced to about -0.01 dex by our adopted treatment. Plausible CNO mismatch adds a zero-point sensitivity of 0.04-0.05 dex, and regularised covariance changes association means by at most 0.03 dex. The preferred analysis gives a mean [Z] = -0.10 dex with a star-to-star dispersion of 0.09 dex. Extended low-resolution J-band spectra therefore yield recoverable metallicity and temperature constraints once continuum placement and weak parameter directions are controlled, although the absolute scale remains conditional on the adopted atmosphere and line-formation framework.
http://arxiv.org/abs/2610.23852 (pdf) (16)
We present early-time optical photometry and spectroscopy for a sample 27 previously unpublished type II supernovae (SNe II) observed by the Precision Observations of Infant Supernova Explosions (POISE) collaboration, combined with a literature sample of 59 events to form a total sample of 86 objects. The resulting light curves comprise near-daily cadence, percent-level precision uBgVri photometry. From spectra acquired within four days of the estimated explosion time, transient narrow emission lines of highly ionised elements, referred to as `flash features&apos, can be identified. Based on the presence or absence of these lines, the SNe were categorised. We estimate key photometric properties, including rise times, peak absolute magnitude, and colour evolution. Our tests indicate that flash SNe are, on average, more luminous at the peak than their non-flash counterparts, which is consistent with previous studies. Interestingly, statistically significant differences were found in the rise times in the B band, with flash SNe exhibiting a tail towards longer rise times. Flash SNe are generally bluer than non-flash events, with distinct colour distributions at 20 days post explosion, likely reflecting intrinsic differences rather than host-galaxy extinction. Colour-curve breaks within 30 days of explosion provide a clear separation between flash and non-flash SNe. The change in slopes across the break, combined with the absolute peak magnitude, yield a strong diagnostic of the flash-spectroscopic properties. This indicates that photometric observations alone may be used to identify SNe experiencing an interaction with dense circumstellar material, even in the absence of early-time spectroscopy.
http://arxiv.org/abs/2607.18315 (pdf) (16)
We report on two deep XRISM observations of the central and northern regions of Hydra-A&aposs X-ray atmosphere, covering the bubbles inflated by jets from the central galaxy&aposs active galactic nucleus. We use spatial-spectral mixing, combining Chandra&aposs high spatial resolution with XRISM&aposs high spectral resolution, to investigate atmospheric kinematics. The velocity dispersion in the northern region, σ_v = 140^+30_-20 km s^-1, is comparable to that in the central region (σ_v = 162 ± 10 km s^-1). We show that motion of the large-scale cocoon shock front could be responsible for the large dispersion toward the north. The dispersion in the northeast quarter of the central pointing, σ_v = 260 ± 50 km s^-1, is among the highest measured. This region contains an X-ray-bright feature previously identified as metal-rich, possibly gas uplifted in the wake of previous-generation cavities. Dispersions in all other quarter regions are low (σ_v ≤ 180 km s^-1) and consistent with previous XRISM results from other objects. The kinetic energy at the center is comparable to the enthalpies of the cavities, while in the north it is roughly an order of magnitude smaller. The jet thus drives gas motion efficiently at smaller scales (r < 95 kpc) and inefficiently at larger scales (95-317 kpc toward the north along the jet). A bulk flow toward our line of sight of -100 ± 30 km s^-1 in the southwest quarter of the central pointing, measured with respect to the central galaxy, is also observed, possibly due to sloshing. No significant bulk motion was detected in any other region.
http://arxiv.org/abs/2606.14849 (pdf) (16)
Motivated by the recent measurements by the Dark Energy Spectroscopic Instrument (DESI) which suggest a late-time matter density approximately 5
http://arxiv.org/abs/2610.23828 (pdf) (15)
Galaxy-galaxy lensing (GGL) is a critical probe of cosmology and galaxy formation. However, predicted GGL amplitudes tend to be ∼ 20-40
http://arxiv.org/abs/2610.22648 (pdf) (14)
The CMB power spectra are a strong probe of the primordial power spectrum (PPS) of the initial inflationary scalar perturbations. There are many proposed models, from the simpler power law form to various well motivated physical models that predict features in the PPS beyond the power law. Given the large parameter space of possible models, it is necessary to have model independent searches for such features. Weak gravitational lensing alters the CMB temperature power spectrum and can bias free form reconstructions of the primordial power spectrum (PPS) if not accounted for precisely, such as when using estimators like the Richardson-Lucy (RL) deconvolution algorithm. We build on our design of the Non linear Iterative Richardson Lucy (NIRL) estimator that updates the lensing contribution using the freeform PPS recovered at the previous iteration, making it truly model independent in the PPS domain while accounting for lensing. Inspired by signal processing and Fourier filtering techniques from radio interferometry, we also design a Fourier space iterative low pass filter (FLIP) to suppress noise overfitting in NIRL reconstructions. We demonstrate the performance of NIRL-FLIP on a range of simulated PPS models with bootstrap estimates of reconstruction uncertainty obtained from simulated realizations, and show that both our estimator and regularization are robust, accurate and precise, and are able to both reduce lensing induced reconstruction bias and suppress noise overfitting by a significant factor in Fourier space. We provide a detailed description of these results and recommend its utility for future model independent searches for hints of novel physics in the form of deviations from a power law PPS. This can further assist in precise and accurate testing of inflation models, in conjunction with joint cosmological parameter estimation.
http://arxiv.org/abs/2610.23564 (pdf) (13)
Motivated by curiosity as to whether the observed Universe at low redshifts is well described by the concordance cosmological model, we investigate the large scale angular clustering of extragalactic neutral hydrogen HI sources in the local Universe. Using the second data release of the FAST All Sky HI Survey (FASHI), we analyze a sample of 142,204 HI sources up to z<0.09, with a median redshift of z̅ = 0.034, covering ∼ 40
http://arxiv.org/abs/2610.23414 (pdf) (13)
We present optical photometric and optical/near-infrared spectroscopic observations of SN 2025wiu, a Type Ib supernova (SN) at z = 0.046927. The r-band light curve exhibits an initial maximum of approximately -18.2 mag, followed by a weaker optical rebrightening reaching approximately -17.8 mag at +38 days in the rest frame relative to the first maximum. However, the reconstructed bolometric light curve does not show a distinct second peak, but instead exhibits a reduced decline rate during the corresponding phase. The spectra consistently show Type Ib features without narrow emission components associated with circumstellar material. The late-time bolometric decline rate of 0.054 +/- 0.008 mag/day is faster than expected for complete trapping of 56Co decay. A radioactive-decay model for the first peak requires nickel mass = 0.47 solar masses, above the range typical of SNe Ib. A combined radioactive-decay and circumstellar-interaction model reproduces the full bolometric light curve with nickel mass 0.38 solar masses, ejecta mass 4.7 solar masses, and an effective circumstellar mass of 0.10 solar masses, with the interaction contributing from the first peak onward. The lack of interaction signatures in the spectra favours an aspherical, disk-like circumstellar structure engulfed by the ejecta, in which the interaction energy is thermalized below the ejecta photosphere; the fitted circumstellar parameters are therefore effective values rather than direct measurements of the circumstellar geometry. SN 2025wiu demonstrates that optical rebrightening does not necessarily imply a genuine bolometric secondary maximum and highlights the importance of multi-band observations for identifying the powering mechanisms of double-peaked stripped-envelope supernovae.
http://arxiv.org/abs/2610.23425 (pdf) (12)
Blazars - active galaxies with the jet pointing at Earth - emit across all electromagnetic wavelengths. As the one-zone model becomes increasingly inadequate to describe blazar emission, more elaborate models are needed. In order to self-consistently describe the entire electromagnetic spectrum, extended jet models are necessary. Notably, kinetic descriptions of extended jets can provide the temporal and spatial evolution of the particle species and the full electromagnetic output. ExHaLe-jet is a lepto-hadronic, extended-jet code allowing for in-depth studies of the particle and radiation evolution on all scales. Here, we present our ongoing study regarding acceleration scenarios of particles within the large-scale jet structure. A single acceleration of particles at the jet base (or elsewhere) is insufficient to produce sufficient photon flux. Thus, re-acceleration through, e.g., consecutive shocks along the jet is needed. We explore two possible scenarios for this.
http://arxiv.org/abs/2610.23218 (pdf) (12)
We investigate why numerical relativity simulations of dust collapse in a matter-dominated era frequently fail before primordial black hole formation due to numerical limitations and inappropriate gauge choices. Using 3-dimensional fully non-linear simulations of a dust fluid, we detail the advantages of advanced numerical methods and the impact of the lapse and shift evolution equations. The challenge is twofold: 1) conventional lapse slicing evolution equations, highly effective in typical black hole simulations, must be adapted in a cosmological context. We demonstrate how different proposed lapse evolution equations perform, enabling fair comparison, and propose our own. 2) The shift evolution equation&aposs role is to spread out the numerical grid observers for a stable simulation. However, the Gamma-driver evolution fails when confronted with shell-focussing of the dust fluid. We show how the shift damping choice depends on the resolution of the shell-focusing region, suggesting that future work should have spatially varying damping. Our improved gauge choices allow simulations to evolve significantly further and potentially reach black hole formation. The best-performing simulations reach close proximity to a black hole; however, further work on the gauge choice is required for this to be stable.
http://arxiv.org/abs/2610.23713 (pdf) (10)
We assess the dark matter interpretation of the 20 GeV gamma-ray halo excess reported by Totani and its proposed connection to the LUX-ZEPLIN (LZ) high-energy nuclear-recoil event, which independently favors dark matter masses above about 200 GeV. Using fifteen years of Fermi-LAT data, we reproduce the excess under the original assumptions and test it against 281 Galactic diffuse emission configurations, three descriptions of the Fermi bubbles, four of Loop I, and uniform and X-ray-brightness-weighted eROSITA bubbles templates. A component correlated with a Navarro-Frenk-White (NFW) annihilation template remains with finer spatial binning and the instrumental response included. Outside the Fermi bubbles footprint, it is detected at nominal significances of about 4-16σ over 5.6-130 GeV, falling to 1-7σ when the enlarged boundaries of a recent substructured bubbles model are also excluded. With the substructured model, which improves the fit substantially, the component remains at 7-18σ in four of six representative foreground configurations. However, its amplitude and spectrum depend strongly on the foreground model, and no tested foreground yields a single NFW normalization across the sky, including with additional spatial freedom in the eROSITA bubbles. Simulations show that foreground mismodeling can produce spurious NFW-correlated emission of the observed amplitude, and that residual-derived bubbles templates can bias the recovered spectrum. The 420 GeV W^+W^- spectrum motivated by the LZ connection describes the data for some foreground models but is disfavored by others. These results disfavor interpreting the full excess as annihilation in a smooth NFW halo and weaken the case for a common origin with the LZ event, but they do not exclude a dark matter contribution to the persistent halo-like residual. [abridged]
http://arxiv.org/abs/2610.23393 (pdf) (10)
We show that horizon-scale magnetic fields around supermassive black holes (SMBHs), only limited by gravity, yield a universal, model-independent ceiling on electromagnetic energy density. For rotating SMBHs, this field enables particle acceleration up to ∼10^27 eV set solely by fundamental constants. We identify a parameter space in accretion rate and magnetic field strength in which quiescent nuclei simultaneously satisfy the transparency conditions for Ultra-High-Energy Cosmic Rays (UHECR) escape and the field requirements for reaching Ultra-High Energies (UHE), naturally complementing photon-rich Active Galactic Nuclei (AGN), which favor neutrino production but suppress UHECR escape. We observe that, when including curvature energy losses, the energy of the primaries scales weakly as E_curv∝M^1/4 and allows to achieve UHE across different compositions and SMBHs masses (10^6–10^10 M_⊙), without requiring extreme spins or super-Eddington accretion. A cross-correlation of the Pierre Auger Observatory (PAO) public data with galaxies hosting reliably measured SMBH masses, separating AGN from quiescent hosts, reveals marginal local excesses (≲2.8σ) consistent with isotropy after correcting for the look-elsewhere effect. The framework provides a direct link between black-hole physics and the maximum energies of cosmic particles, and defines concrete observational tests for next-generation UHECR datasets.
http://arxiv.org/abs/2610.23326 (pdf) (10)
Context. Long-term photometric surveys of the Galactic bulge are vital for characterizing its stellar populations and structure, providing a foundation for next-generation facilities. Aims. We present initial results of the GalActic bulGe with plEiadi (GAGE) project, focusing on variable star discovery in Baade&aposs Window. Methods. Using archival DECam (ugriz) data on the Blanco 4-m Telescope spanning six years, we applied high-precision PSF photometry. We developed a pipeline combining Generalized Lomb-Scargle analysis and polynomial light-curve modelling to derive mean multi-band magnitudes and amplitudes. Results. Quality cuts yield 6,304 high-confidence variable sources. We expand multi-band coverage for 4,568 sources in common with OGLE IV and uncover 1,736 ( 28%) new candidates absent from OGLE IV. Using a Random Forest classifier, we categorize these sources, finding that 1,391 were also missing from both Gaia DR3 and VVV. This population is concentrated at the fainter end (i >= 17.5 mag), enabled by our deeper detection limits. Conclusions. Our pipeline effectively probes deep, highly crowded fields, expanding the variable star inventory of the inner Galaxy as a benchmark for Rubin-LSST.
http://arxiv.org/abs/2610.23258 (pdf) (9)
The CORNISH survey provided the most complete sample of ultra-compact H II (UCHII) regions in the northern Galactic plane to date, but at its 1.5 arcsec resolution the most compact, and plausibly youngest, members of the sample could not be morphologically classified. We present Karl G. Jansky Very Large Array (JVLA) A-configuration observations at 10 GHz of 69 such CORNISH UCHII regions, reaching an angular resolution of 0.2-0.3 arcsec and a typical rms noise of 10 microJy/beam. All but two of the imaged targets are detected; the two non-detections lie in fields that are dynamic-range limited by bright neighbours. Six targets ( 9 per cent) are resolved into pairs of sources, and a further 15 CORNISH UCHII regions detected within the fields reveal complex, multi-component structure. We also report candidate massive young stellar object (MYSO) jet knots and a row of three compact sources coincident with a previously unresolved infrared-dark H II region. The targets have a median deconvolved angular diameter of 0.9 arcsec and physical diameters of 0.004-0.17 pc (median 0.03 pc), with no correlation between angular size and distance, consistent with a population of intrinsically compact, young regions. Spectral indices between 1.4 and 10 GHz place 32 per cent of the targets below the thermal free-free limit, possibly indicating non-thermal emission or variability. These data enable morphological classification of nearly the entire CORNISH UCHII sample, which will be presented in a forthcoming paper.
http://arxiv.org/abs/2610.23825 (pdf) (9)
The evidence, 𝒵 = p(x | ℳ), of a data-generating model ℳ is a powerful tool in Bayesian frameworks: it enables model inference, model misspecification tests, and goodness-of-fit assessment as an alternative to standard discrepancy statistics, e.g. χ^2. While versatile, the evidence is notoriously difficult to compute, even when the likelihood is available. In both likelihood-based and simulation-based inference, estimating 𝒵 provides an immediate solution to high-dimensional goodness-of-fit or model misspecification tests: using 𝒵(x) as the test statistic induces a minimum-volume level-α test. Such tests are essential to assess model fidelity, which is necessary to trust simulation-based parameter inference. Existing learned approaches largely target Bayes factors, which only measure evidence ratios. Instead, we provide accurate, fast, and amortized estimates of the absolute evidence. We present ENFR: Evidence Networks with Flow-matched References. Building on the Evidence Network framework of Jeffrey Wandelt (2024), the proposed method first learns Bayes factors with respect to a flow-matched reference, from which the absolute evidence is recovered. We show improvement in estimation of 𝒵 over the flow-based baseline in all applications considered and the best performance of all density estimation methods considered in our benchmarks. We furthermore introduce and provide the proof for an architecture that jointly learns the individual evidences of a set of models {ℳ_i}. We demonstrate the performance and scaling behavior of these approaches on data from a variety of scientific applications.
http://arxiv.org/abs/2610.23497 (pdf) (9)
DESI DR2 favours dark energy whose equation of state evolves crossing w=-1 in the past. Coupled quintessence can produce such an apparent phantom crossing without a phantom field, since a coupled dark sector analysed under the assumption of separate conservation yields an effective equation of state below -1. We test whether this mechanism can reproduce the DESI result, using two thawing potentials, a cosine and a double-exponential hilltop, confronted with DESI DR2 BAO, Planck PR4 and three supernova compilations. Both models improve on ΛCDM by Δχ^2 = 5–15, recovering 64–92
http://arxiv.org/abs/2610.22790 (pdf) (9)
We report the independent discovery of the cataclysmic variable (CV) SPLUS J102139.31-282247.5 using Southern Photometric Local Universe Survey 12 band photometry, following its initial identification via ASAS-SN variability. Here, we investigate the system using optical photometry alongside optical and X-ray spectroscopy. High-cadence optical light curves from publicly available TESS data and a dedicated campaign at the Pico dos Dias Observatory in Brazil revealed a highly inclined system exhibiting pronounced eclipses with an orbital period of 1.81 hours. We found no evidence of a WD spin period. Despite seasonal gaps, the long-term ASAS-SN light curve revealed several outbursts (15 with Δmag ∼ 1.5 brighter than the mean magnitude) and one superoutburst over the past decade. During Swift X-ray telescope observations, the system exhibited X-ray luminosity of 8.5 × 10^31 erg s^-1 in the 0.3-10 keV band, originating from thermal plasma emission at a temperature consistent with accretion shocks. Optical spectroscopy from the Gemini Multi-Object Spectrographs on the Gemini South telescope confirms the CV nature of the system, showing strong, double-peaked H and HeI emission lines alongside a weak HeII4686 line. Altogether, these features strongly support the classification of SPLUS J102139.31-282247.5 as a quiescent, eclipsing dwarf nova with an unusually high X-ray luminosity.
http://arxiv.org/abs/2610.22488 (pdf) (9)
Giant bubbles and cavities are widespread in star-forming galactic discs, yet the stellar dynamics associated with such structures remain poorly constrained. Here we combine Gaia Data Release 3 with three-dimensional extinction information to identify 3,502 massive-star candidates corresponding to initial masses M_ini≳8 M_⊙ within 5 kpc of the Sun. A search of the stellar-density field alone, independent of interstellar-medium morphology and stellar kinematics, reveals two prominent overdensities, R1 and R2, with characteristic scales of approximately 1.3 kpc. Both are subsequently found to be associated with previously identified giant interstellar cavities in projected position and line-of-sight location. After subtracting locally fitted background velocity fields, stars with full three-dimensional velocities reveal coherent divergent in-plane motions across both structures. R1 extends well beyond a known massive-star association and exhibits strongly anisotropic expansion, with one principal strain rate robustly positive and the other consistent with zero. R2 instead exhibits positive strain rates along both principal directions, indicating broadly two-dimensional internal expansion. These results show that density-selected, cavity-associated stellar structures of comparable scale can host coherent kiloparsec-scale flows with contrasting internal flow geometries and need not share a common stellar kinematic organisation.
http://arxiv.org/abs/2610.22470 (pdf) (9)
Gravitational lensing by rotating black holes is modeled, and the null geodesic equations in the strong field regime are integrated numerically. Observational signatures of the black hole spin parameter are identified by comparing image locations and shapes of extended sources at different a values. For the black hole at the center of our galaxy, angular resolutions of the order of hundreds of micro arc seconds are needed to distinguish between the true source position and gravitationally lensed image position. Angular resolution below one micro arc second is needed to determine the spin parameter of the galactic black hole using gravitational lensing.
http://arxiv.org/abs/2610.23704 (pdf) (8)
We present a comprehensive Sun-to-Earth analysis of the extreme space weather event of 18 -19 January 2026, hereafter the Martin Luther King (MLK) storm, which produced a severe (G4) geomagnetic storm and one of the most intense solar radiation storms of Solar Cycle 25. The geomagnetic disturbance was driven by a fast Earth-directed coronal mass ejection (CME) associated with an X1.9-class flare from NOAA Active Region 14341. We combine multi-viewpoint remote-sensing observations, in situ measurements from Solar Orbiter and Wind, and physics-based magnetohydrodynamic simulations using the CORHEL-CME framework to reconstruct the eruption from the low corona to 1 au and identify the mechanisms governing its geoeffectiveness. Observations and simulations indicate that the eruption originated from a right-handed magnetic flux rope embedded in a complex magnetic environment dominated by a neighboring positive-polarity coronal hole. The simulations reproduce the observed CME morphology, asymmetric expansion, magnetic fields measured at both spacecraft, and Dst index evolution. The main discrepancy is attributed to filament material not included in the model. Favorable magnetic polarity alignment between the CME and coronal-hole field minimized magnetic erosion through interchange reconnection, enabling rapid escape and preservation of the strong southward magnetic field that drove the geomagnetic storm. The prolonged recovery phase resulted from interaction between the trailing CME and a following high-speed stream from the same coronal hole. These results demonstrate that storm intensity was primarily established by the coronal magnetic configuration, whereas storm duration was controlled by the ambient solar wind.
http://arxiv.org/abs/2610.23089 (pdf) (8)
Geometric beaming is an often invoked explanation for the apparent super-Eddington luminosities of ultraluminous X-ray sources. It is a natural consequence of the geometrically thick envelope that is expected to surround the inner accretion disk in the super-Eddington accretion regime; as the radiation can only freely escape through a narrow funnel-shaped cavity, the flux along the funnel axis is enhanced. A fraction of the photons will nevertheless escape the funnel cavity at larger angles, reducing the effectiveness of the beaming. To explore the limits of geometric beaming, we performed Monte Carlo simulations of conservative radiative transfer in an infinitely optically thick funnel while accounting for polarization. We discovered that the beaming becomes inefficient when the funnel opening angle is small (≲15^∘) because the photons tend to remain trapped within the cavity until they reach the upper edge of the funnel where they can escape at larger angles. For most of the photons to escape a narrow funnel through a small solid angle, the funnel has to be unrealistically long. Consequently, geometric beaming could at most boost the apparent luminosity by a factor of a few dozen in our idealized simulations. While larger boost factors on the order of ∼100 are possible for funnels with a length of ≳10^5 Schwarzschild radii, this requires very narrow opening angles of ≲3^∘. Additionally, we found that aberration due to mildly relativistic velocities does not significantly change this result. The maximal polarization degree of light reflected by the funnel remained on the order of a few percent for face-on sources and ≲10
http://arxiv.org/abs/2610.23359 (pdf) (8)
Blazars are typically divided by spectral appearance; may that be strength of optical emission lines or the frequency of the synchrotron peak. However, this may say nothing about the variability properties of the various types. Here, we study the long-term variability properties with nearly 20 years of multiwavelength data of one flat-spectrum radio quasar, PKS 1510-089, and one high-frequency-peaked BL Lac object, PKS 2155-304. Both exhibit strong variability with an energy-dependent variability amplitude. Using the total data sets, we do not find any correlation between the various bands. While the flux distributions of PKS 2155-304 are compatible with log-normal functions, this is not the case for PKS 1510-089. We conclude that the one-zone model is inadequate to describe the long-term evolution of blazars, requiring multi-zone or stratified models. Nevertheless, the underlying jet physics appear similar between the two types and in both sources the variability is, in fact, variable.
http://arxiv.org/abs/2610.22653 (pdf) (8)
Magnetic reconnection around neutron stars (NSs) forms small-scale currents while transferring energy to electron-positron pairs. These dissipating currents cover a wide range of scales, potentially leaving an imprint on time-variable X-ray bursts and pulsed radio emission. We present a reduced, planar model of internal reconnection in a magnetic flux tube with relativistic magnetization, finite resistivity, and strong guide field. Comparing with realistic NS parameters, this model represents the dynamics of a slender magnetic flux tube that is much thinner than the stellar radius. The seed current profile is smooth and does not spontaneously reconnect; instead, a discrete and localized current perturbation is introduced, representing e.g. localized energy release from a NS crust or magnetosphere. Three high-resolution relativistic magnetohydrodynamics simulations are described: one has uniform resistivity and Lundquist number such that reconnection proceeds in the plasmoid-mediated regime; the other two account for the feedback of charge starvation on the resistivity in current layers, where density is depleted; we provide a more general explanation for this depletion. All simulations produce a few dominant magnetic X-points and a secondary Kelvin-Helmholtz instability that excites high-wavenumber currents. We quantify the reconnection rate, the intermittency of dissipation, and the time profile of current spikes. We further present three-dimensional simulations confirming a limited growth of the ideal kink mode in the strong guide field regime. Applications are discussed, including to radio pulse microstructure and current damping below the magnetar surface during an X-ray flare, leading to mass loss and X-ray afterglow.
http://arxiv.org/abs/2610.23410 (pdf) (7)
With the rapid growth of large-scale time-domain surveys, vast amounts of photometric and spectroscopic data are continuously being released. A substantial number of contact binaries have been identified, accompanied by both light curves and radial velocity observations. Traditional modeling of these datasets relies on computationally intensive physical models, which limits their applicability to large-scale analyses. To this end, this paper proposes the use of a fully connected neural network and a SIREN (Sitzmann et al. 2020) model to rapidly generate radial velocity curves. After training, the models can rapidly generate radial velocity curves with high accuracy, significantly reducing computational costs and enabling efficient large-scale studies. This method was further integrated with the light curve computation framework developed by Li et al. (2025), allowing simultaneous fitting of radial velocities and light curves. To assess the reliability of the proposed model, synthetic observational datasets were constructed and used for parameter estimation. The deviations between the derived parameters and the ground-truth values were analyzed quantitatively, enabling a comprehensive evaluation of the model&aposs accuracy.
http://arxiv.org/abs/2610.22662 (pdf) (7)
The Event Horizon Telescope (EHT) has resolved horizon-scale emission in M87* and Sgr A*. Extending these observations to fainter black holes, multiple frequencies, and evolving structures requires greater sensitivity and new observing capabilities. We investigate how ALMA2040 could expand the scientific reach of global millimeter very long baseline interferometry (mm-VLBI), compared with ALMA&aposs Wideband Sensitivity Upgrade (WSU). Phasing a larger ALMA array and enabling simultaneous observations across multiple ALMA bands would improve sensitivity and frequency phase transfer, providing spectral and polarimetric diagnostics of gravity, accretion flows, magnetic fields, and relativistic jets. In the configurations studied, ALMA2040 achieves sub-mJy baseline sensitivity at 230 GHz and improves prospects for faint-source imaging at 345 GHz, opening a path from individual horizon-scale targets toward population studies. Simulated observations of M87* show improved multi-frequency recovery of emission from the compact ring-like region toward the jet base, enabling studies of the accretion-jet connection. For Sgr A*, simulations of interstellar scattering mitigation yield a 27% reduction in the error of the recovered intrinsic source structure relative to the WSU configuration. Beyond black holes, these capabilities benefit geometric-distance measurements, stellar masers, compact transients, and other applications of high-resolution astrophysics. Together with rapid-response capabilities across the VLBI network, ALMA2040 could support a transition from isolated snapshots toward multi-frequency, time-resolved studies of compact sources.
http://arxiv.org/abs/2610.22819 (pdf) (7)
The models of modified gravity that can explain the late-time accelerated expansion of the Universe must pass the local tests of gravity, introducing a screening mechanism so that any additional fifth force in gravity is suppressed at stellar scales. Dark matter halos, in contrast, by their extendedness in scale, may exhibit a transition of screening to small scales, providing a probe of modified gravity and the origin of the late-time acceleration. We perform comprehensive studies of dark matter halos with Vainshtein screening in the Horndeski theory. Given that dark matter is cold and a dark matter halo is spherically symmetric, we analytically clarify the conditions for the halo profile under which the Vainshtein screening works effectively. We find that the Vainshtein screening is less effective even within the Vainshtein radius as the profile approaches a core-like structure around the center of a halo, which motivates us to postulate that some features owing to the fifth force can survive during the dynamical evolution of a halo. We numerically examine our postulate by a N-body simulation for an isolated halo with the assumption of Newtonian two-body interactions with the spherical external fifth force. We make use of three proxies of the mass through hydrostatic dynamics, the Jeans equation, and gravitational lensing, all of which must be identical in the case of Newtonian gravity as a weak-field approximation to General Relativity. We find that these proxies can be biased relative to each other at a (quasi) equilibrium, quantitatively indicating that the fiducial model of gravity is not Newtonian. Our result affirmatively suggests that multiple mass proxies from observations are capable of testing gravity on halo scales.
http://arxiv.org/abs/2607.04901 (pdf) (7)
(abridged) We develop a flexible analytical framework for modeling blastwave evolution in arbitrary environments and use it to investigate how large-scale galactic structure affects the dynamics and morphology of evolved SNRs and SBs. We introduce FRANZ (FRamework for ANalytical one-Zone blastwave dynamics), a modular thin-shell model that follows the local evolution of a shock-surface segment in environments characterized by arbitrary density, velocity and gravitational fields. After validating the model against well-established analytical results, we apply it to study the effects of vertical stratification, galactic shear and dense galactic substructure on blastwave evolution. FRANZ reproduces the classical evolution of blastwaves in uniform media while extending to complex environments. We derive criteria for disk break out in stratified media, characterize the timescales on which differential rotation deforms blastwaves and identify a new mechanism, by which it can suppress the momentum coupling in continuously-driven blastwaves. Interactions with dense filaments modify both shock-surface morphology and dynamics and confound the interpretation of the expansion history of observed remnants, which depends on the density distribution prior to the onset of explosions, which is fundamentally inaccessible from the observed state. In highly structured media with a high volume-filling factor of diffuse gas, ages inferred from the observed state may be systematically overestimated. FRANZ provides a computationally inexpensive and extensible framework for studying blastwaves in realistic galactic environments. It offers a useful complement to numerical simulations for interpreting observations of evolved SNRs and SBs and for developing improved models of stellar feedback in large-scale simulations.