Most interesting arXiv papers
Wednesday 17-12-2025
http://arxiv.org/abs/2512.14486 (pdf) (143)
We present deep ALMA Band 7 observations of the [OIII] 88μm line and underlying dust continuum emission in four UV-bright, gravitationally lensed (magnification μ= 1.4-3.8), JWST-selected galaxies at z = 8.5 - 10.3, with observed magnitudes -22.5 ≲ M_UV≲ -20.5. [OIII] 88μm is confidently detected in UNCOVER-10646 at z=8.5080 ± 0.0011 (15σ) and DHZ1 at z=9.3113 ± 0.0005 (6σ), with both being intrinsically luminous systems [L_[OIII] = (1.1 - 1.6) × 10^9 L_⊙] that follow the local [OIII]-SFR relation. [OIII] 88μm remains undetected in the two z>10 targets, including in the z=10.07 X-ray AGN UHZ1, where we obtain a deep limit of L_[OIII] < 6 × 10^7 L_⊙. Dust emission is not detected in any individual source nor in a stack (<3σ). The high S/N [OIII] 88μm detection in UNCOVER-10646 uniquely reveals an additional broad component (FWHM = 1366_-329^+473 km/s; ΔBIC≈20) indicative of an ionized outflow. We infer a high outflow rate of Ṁ_out = 128_-46^+80 M_⊙yr^-1, corresponding to a mass loading factor η= Ṁ_out/SFR = 2.9_-1.0^+1.8 that is broadly consistent with theoretical predictions and JWST-based studies of ionized outflows at high redshift. While high-resolution ALMA follow-up is required to confirm the outflow, this first systematic study at z>8 highlights the unique diagnostic power of [OIII] 88μm in characterizing galaxies in the early Universe.
http://arxiv.org/abs/2512.13778 (pdf) (85)
Lyman α emission is highly sensitive to dust and neutral hydrogen and is therefore expected to be strongly suppressed in dusty or gas-rich galaxies during the epoch of reionization (EoR). Nevertheless, numerous moderately dusty Lyα emitters (LAEs) are observed at this epoch, suggesting that complex interstellar medium (ISM) geometries and feedback-driven outflows may facilitate Lyα escape. We investigate the dust, gas, and stellar properties of the gravitationally lensed LAE HCM 6A at z=6.5676 to characterize its multiphase ISM and the physical conditions regulating Lyα escape. We combine JWST/NIRISS slitless spectroscopy, HST+JWST/NIRCam imaging, and JWST/NIRSpec slit spectra from the CANUCS program. Using a customized BAGPIPES SED-fitting framework with a flexible attenuation law, we derive stellar, nebular, and dust properties on integrated (∼kpc), slit-level (∼0.1 kpc), and pixel-level (∼25 pc) scales, enabled by strong lensing (μ≈ 8.3-9.1). A Lyα map from SLEUTH traces the spatial distribution of Lyα emission. We measure an unlensed stellar mass of log M_∗ = 8.3-8.4 and an intrinsic UV magnitude of M_UV = -19.8 ± 0.1. The older region (S1) is moderately dusty with consistent stellar and nebular attenuation indicators, implying a uniform ISM geometry, while the youngest region (S3) shows strong discrepancies among dust tracers, indicating a feedback-shaped multiphase ISM. Lyα emission arises primarily from S3, where a dust-cleared central clump enables efficient Lyα escape. We find Calzetti-like attenuation curves with a UV bump that strengthens with stellar age and decreasing A_V. Our observations provide a uniquely detailed, spatially resolved view of a moderately dusty LAE during the EoR, demonstrating how feedback and multiphase ISM structure govern Lyα escape.
http://arxiv.org/abs/2512.14212 (pdf) (65)
Observational campaigns with JWST have revealed a higher-than-expected abundance of UV-bright galaxies at z≳10, with various proposed theoretical explanations. A powerful complementary constraint to break degeneracies between different models is galaxy clustering. In this paper, we combine PANORAMIC pure parallel and legacy imaging along 34 independent sightlines to measure the cosmic variance (σ_CV) in the number count of Lyman break galaxies at z∼10 which is directly related to their clustering strength. We find σ_CV=0.96^+0.20_-0.18, 1.46^+0.54_-0.44, and 1.71^+0.72_-0.59 per NIRCam pointing (∼9.7 arcmin^2, ≲1.5 pMpc at z∼10) for galaxies with M_UV<-19.5, -20, and -20.5. Comparing to galaxies in the UniverseMachine, we find that σ_CV is consistent with our measurements, but that the number densities are a factor ≳5 lower. We implement simple models in the UniverseMachine that represent different physical mechanisms to enhance the number density of UV-bright galaxies. All models decrease σ_CV by placing galaxies at fixed M_UV in lower mass halos, but they do so to varying degrees. Combined constraints on σ_CV and the UVLF tentatively disfavor models that globally increase the star formation efficiency (SFE) or the scatter in the M_UV-M_halo relation, while models that decrease the mass-to-light ratio, or assume a power-law scaling of the SFE with M_halo agree better with the data. We show that with sufficient additional independent sightlines, robust discrimination between models is possible, paving the way for powerful constraints on the physics of early galaxy evolution through NIRCam pure parallel imaging.
http://arxiv.org/abs/2512.13896 (pdf) (42)
Strong gravitational lensing is a powerful probe of cosmology, dark matter (DM), and high-redshift galaxy evolution, but current samples of strongly lensed galaxies (SLGs) remain far too small to exploit its full potential. Herschel&aposs submillimeter (submm) surveys demonstrated that submm selection provides the most efficient and least biased route to identifying high-redshift SLGs, yet produced only a few hundred systems over limited, heterogeneous fields. Achieving the thousands of SLGs required for precision cosmology and detailed studies of distant dusty star-forming galaxies demands a new, wide-area, homogeneous sub-mm survey. A facility like AtLAST, capable of extending Herschel-like methodologies to much larger areas, is uniquely positioned to deliver the order-of-magnitude increase in SLG numbers needed for transformative progress.
http://arxiv.org/abs/2512.13932 (pdf) (41)
Dust and cold gas are not uncommon in nearby early-type galaxies (ETGs), and represent an important aspect of their evolution. However, their origin has been debated for decades. Potential sources include internal processes (e.g., mass loss from evolved stars), external mechanisms (e.g., minor mergers or cooling flows), or a combination of both. Gas-rich minor mergers have long been proposed as an important channel for cold gas fueling in both observations and simulations, but direct evidence of cold gas transportation via gas-rich minor mergers remains elusive, particularly in galaxy groups and clusters where environmental effects are prevalent. In this letter, we present the first unambiguous case of direct cold gas transportation onto a brightest group galaxy (BGG) at z=0.25, driven by an ongoing close-separation gas-rich minor merger with a mass ratio of ∼1:56. High-resolution JWST imaging reveals a heavily obscured, low-mass satellite that is barely visible at restframe optical wavelengths. Tidal stripping from this satellite deposits gas and dust onto the BGG, forming prominent ∼10 kpc dust lanes in situ. Cosmological simulations indicate that such interactions preferentially occur in gas-rich satellites undergoing their first infall in highly eccentric orbits. Our results highlight the pivotal role of gas-rich minor mergers in replenishing cold gas reservoirs and shaping the evolution of central ETGs in galaxy groups.
http://arxiv.org/abs/2512.14066 (pdf) (39)
The James Webb Space Telescope (JWST) has recently identified Abell 2744-QSO1 as a compact, metal-poor, black hole (BH) dominated galaxy at z≃ 7. This system exhibits an extreme black-hole-to-stellar mass ratio and unusually low metallicity, posing significant challenges to BH seeding models. Motivated by these discoveries, we perform high-resolution cosmological simulations with a massive primordial black hole (PBH; M_BH=5×10^7 M_⊙) seed, incorporating for the first time a fully coupled treatment of PBH accretion, BH feedback, and Population III/II star formation and stellar feedback. Although PBHs accelerate structure formation through the seed effect, the associated strong thermal feedback from the accretion delays the onset of star formation to z≲ 10, producing short, bursty episodes throughout the subsequent evolution. PBH-driven outflows expel enriched gas from the nucleus, while sustained inflows from the intergalactic medium continuously replenish pristine material. This feedback-regulated cycle naturally yields low accretion rates (ṁ_BH/ṁ_edd∼ 1-10
http://arxiv.org/abs/2512.13943 (pdf) (38)
The epoch of reionization (EoR), marking the Universe&aposs transition from a neutral to ionized state, represents a pivotal phase for understanding the formation of the first stars and galaxies. Intensity mapping of atomic and molecular lines, such as [CII] and CO J-ladder transitions, across a broad redshift range is a powerful tool for investigating star formation history, metallicity, the distribution of gas and dust, and the physical conditions within galaxies. Additionally, 21 cm line intensity mapping directly probes the neutral hydrogen content in the intergalactic medium, offering a unique window into the timing and morphology of reionization. In this study, we explore the cross-correlation between the 21 cm signal and multi-line intensity mapping (LIM) to forecast their detectability for next-generation experiments. Our analysis emphasizes the complementary potential of these techniques to constrain parameters such as the minimum mass of ionizing sources and the ionization fraction x_e. Cross-correlations with LIM also enable constraints on physical properties like metal enrichment and the relationship between star formation rates and multi-line luminosities. Using mock observations from Square Kilometre Array (SKA)-low 21 cm and Fred Young Submillimeter Telescope (FYST)-like LIM experiments, we find that the [CII]–21 cm cross-correlation can constrain reionization history by measuring x_e across multiple redshift bins with significance levels ranging from 9 to 40σ. We extend our analysis to CO transitions, showing that the CO(1-0)–21 cm cross-correlation provides competitive constraints on reionization parameters. The synergies explored here will enable robust constraints on both reionization and LIM parameters, maximizing the scientific return of current and next-generation intensity mapping experiments.
http://arxiv.org/abs/2512.13881 (pdf) (36)
Submillimetre molecular lines (e.g., CO, HCN, SiO) provide a uniquely powerful view of the physical and chemical processes that govern star formation (SF) and galaxy evolution. Yet, our current picture of the molecular universe beyond the Milky Way remains strikingly incomplete: broad chemical inventories exist for only a handful of galaxies, typically more extreme than the Milky Way, constrained by sensitivity limits and narrow survey strategies. In the 2040s, surveying galaxies with multi-species, multi-transitions observations across diverse galactic environments will be crucial to establish effective chemical diagnostics of the various ISM processes from the early universe to z=0. Extragalactic astrochemistry provides a uniquely sensitive probe of the physical processes shaping galaxies, allowing us to understand, species by species, how gas responds to its local environment and how galaxies grow, transform, and recycle matter over cosmic time.
http://arxiv.org/abs/2512.14281 (pdf) (35)
Galaxy evolution is driven by processes occurring across a wide range of scales, from star formation within giant molecular clouds (parsec scales) to outflows and secular evolution across entire galaxies (kpc scales), and the interplay between galaxies, their dark matter haloes, and large-scale structures (Mpc scales). Connecting the distribution of baryonic matter and energy across these scales will remain one of the key challenges for both theoretical and observational astrophysics in the coming decade. A major development towards meeting this challenge has been the growing ability to obtain highly spatially resolved (parsec-scale) integral-field spectroscopic observations (e.g. with VLT/MUSE), as well as to probe the extremely low-surface brightness outskirts of galaxies at large radii and high vertical scale heights. To combine the two regimes, we need a paradigm shift in the way we do spectroscopy on galaxies, especially considering the ongoing and future photometric surveys. The next decade will also bring a revolution in extensive photometric surveys of large areas of the sky, uncovering low surface brightness features around nearby galaxies. However, to fully understand the processes that connect galaxies to their haloes, shape low surface brightness features, and drive secular evolution, spatially resolved spectroscopy will be essential. Here, we outline the need for wide-field spectroscopic observations of statistically significant samples of nearby galaxies and highlight the key questions that can only be addressed with such data
http://arxiv.org/abs/2512.13822 (pdf) (33)
We investigate how galactic disk structures connect to the detailed properties of their host dark-matter halos using the TNG50 simulation. From the hydrodynamic and matched dark-matter-only runs, we measure a comprehensive list of halo properties describing density structure, angular momentum, shape, assembly history, and environment. Using the morphological decomposition developed in Paper I, we quantify the sizes, scale heights, and mass fractions of the disk components for galaxies at 0 ≤ z ≤ 4. Random Forest (RF) regression shows halo properties alone predict disk size and thickness with high accuracy, while Symbolic Regression (SR) provides compact empirical relations with slightly lower accuracy. Disk height is consistently easier to predict than disk size, and lower-mass halos yield higher accuracy than massive halos. Predictions based on halo properties measured in the full-physics hydro simulations outperform those based on the matched dark-matter-only halos, reflecting the imprint of baryonic restructuring on the inner halo. SHAP analysis reveals that the most informative halo parameters include concentration, Einasto shape, global and inner spin, and recent mass accretion, though their importance varies across disk properties. We show that correlations between disk size and the density-profile shape arise primarily from disk-induced modification of the inner halo, rather than a primordial connection. Finally, we point out that disks become more extended with respect to their host halos at higher redshift in low-mass halos, while massive high-redshift halos show the opposite trend. We provide SR-based prescriptions that accurately map halo properties to disk structures, offering practical tools for galaxy-halo modeling.
http://arxiv.org/abs/2512.14599 (pdf) (32)
We investigate the quenching of galaxies using a mock observational lightcone generated from the Semi-Analytic Model (SAM) L-Galaxies, closely matched to observations from the UKIDSS Ultra Deep Survey (UDS). The sample is used to study merging, rejuvenation, and visibility times for star-forming, quiescent, and post-starburst (PSB) galaxies, to assess the impact on the build-up of the passive galaxy mass functions. We find, for example, that a typical PSB (M_∗∼10^10 M_⊙) at z≈1 has a 15 per cent likelihood of merging and around a 25 per cent likelihood of rejuvenating within 1 Gyr of being identified. Applying these rates and timescales to the observational data, we estimate the fraction of quiescent galaxies that passed through a PSB phase. We find that 18 - 28 per cent of the build-up in the massive end (M_∗>10^10 M_⊙) of the passive mass function at 1<z<2 can be explained by PSBs, with the contribution declining to ∼5 per cent by z ≃ 0.5. Accounting for mergers and rejuvenation reduces the inferred PSB contribution by approximately a factor of two. At lower stellar masses (M_∗ < 10^10 M_⊙), rapid quenching through a PSB phase explains a significantly larger fraction of the growth in the passive mass function. With a visibility time of ∼ 0.75 Gyr, we find that around 60-80 per cent of low-mass passive galaxies underwent a PSB phase. Our findings provide further evidence that low- and high-mass galaxies follow different quenching pathways.
http://arxiv.org/abs/2512.13783 (pdf) (32)
We use the Dark Energy Spectroscopic Instrument Data Release 1 (DESI DR1) to identify stellar members of the Local Group dwarf galaxies. We cross-match DESI targets with candidate members that are based on Gaia proper motions, positions, and photometry. The addition of DESI radial velocities enables secure membership determination in 16 systems, spanning ultra-faint, classical, and irregular dwarf galaxies. Our results confirm that Gaia-based selection algorithms are effective in minimising foreground contamination. With DR1, we find the largest number of member stars for Draco and Sextans. While Draco is discussed in detail in other works, we focus on Sextans and, for the first time with DESI, trace its stellar kinematics to large radii (up to ∼10 half-light radii). Our derived systemic velocity and velocity dispersion are consistent with previous works. We find that the more metal-poor population exhibits a higher velocity dispersion and extends to larger radii, whereas the more metal-rich population is kinematically colder and centrally concentrated. The metallicity gradient is steeper in the inner regions of Sextans (∼ -9× 10^-3 dex arcmin^-1 or ∼-0.29 dex kpc^-1), while almost no gradient in the outskirts, hinting for an ex-situ halo or for an “outside-in” star formation. Although DESI [α/Fe] ratios for Sextans stars with [Fe/H] ≳-2.0 are in line with literature values, those for very metal-poor stars ([Fe/H] ≲-2.0) present a large scatter and strong anti-correlation with metallicity, warranting a caution for using DESI abundances in this regime. We find one ultra metal-poor candidate in Canes Venatici I, which requires higher resolution spectroscopic follow-ups.
http://arxiv.org/abs/2512.14646 (pdf) (31)
Increasing evidence shows that AGN with radio detections have more perturbed ionized gas kinematics and higher outflow detection rates, suggesting a link between radio emission and these processes. In galaxies with weak or ambiguous AGN signatures, some studies attribute the radio emission to star formation, while others propose AGN-driven winds or weak, unresolved jets as the dominant mechanism. To investigate this connection, we take a step back and analyze a sample of star-forming (SF) galaxies with no clear current AGN signatures. Using low-(LOFAR, 144MHz) and high-frequency (FIRST, 1.4GHz) surveys, combined with spatially resolved spectroscopy from the MaNGA survey, we compare SF galaxies with 144 MHz detections that either do or do not have GHz detections. Despite being matched in stellar mass, redshift, and radio (MHz) luminosity, GHz-detected SF galaxies systematically differ from their non-GHz-detected counterparts. The former display enhanced ionized gas-emission line widths, higher central outflow fractions, redder colors, increased central obscuration, and offset emission-line ratios that shift towards (or closer to) the AGN regime (in the [NII] BPT diagram). Furthermore, the non-GHz galaxies are likely undetected due to their extended radio morphologies, while the GHz-detected ones are significantly more radio compact. Most of the properties from the GHz-detected (compared to non-detected) remarkably resemble the behavior found in many studies of radio-detected AGNs. This suggests that the underlying physical mechanisms shaping GHz-detected SF galaxies properties are fundamentally similar. This raises intriguing questions about whether some compact SF galaxies represent a precursor phase of AGN evolution or a form of low-power AGN. The radio compact characteristic sizes of GHz-detected SF galaxies also suggest a connection between AGN and old starburst galaxies.
http://arxiv.org/abs/2512.13801 (pdf) (31)
A critical step to create realistic mock catalogs that support large-scale photometric and spectroscopic sky surveys is the production of cosmological simulations that accurately model the survey observables taking into account the redshift-dependent galaxy formation and evolution processes. Here we develop an efficient framework, HOLCon (Halo Optimised Lightcone Constructor), for post-facto construction of dark matter halo lightcones from simulations and use them to generate a mock galaxy catalogue. HOLCon includes a module to optimise the lightcone&aposs orientation within the simulation box, minimising repeated structures when the survey volume exceeds a single box – a common challenge in modern surveys. A linear interpolation scheme tracks the evolution of halo properties across snapshots. Applied to the publicly available Uchuu simulation, we construct a lightcone of 50 deg^2 and extending up to z = 10, providing representative coverage of deep fields of Stage IV surveys. We validate the lightcone for cosmological applications by comparing the dark matter halo clustering in the lightcone with those from the original simulation snapshots. Subsequently, we make the galaxy-halo connection on the lightcone with a redshift extended version of the SciPIC algorithm producing a comprehensive set of descriptive galaxy attributes. HOLCon leverages Dask, a scalable parallel computing pythonic framework for fast construction of dark matter halo lightcones enabling rapid creation of multiple statistical realizations essential for robust cosmological inference. The produced galaxy mock makes predictions for clustering of Hα emitters, making it a useful cosmology resource.
http://arxiv.org/abs/2512.14372 (pdf) (30)
Present-day galaxy clusters are the largest virialized structures in the Universe, yet their early assembly remains poorly understood. At z>2, clusters in formation span tens of Mpc and host gas-rich, dust-obscured galaxies embedded in extended, low-surface-brightness gaseous environments. Current (sub-)millimeter facilities lack the mapping speed, sensitivity, and contiguous field of view needed to trace the cold gas and dust driving rapid galaxy growth across such scales. A future large single-dish observatory with degree-scale coverage, broad spectral access, and high-multiplex capability would enable comprehensive and uniform mapping of entire protoclusters, revealing where star formation is triggered or quenched, and quantifying the cold gas budget, thus providing information on gas cooling within protocluster environments. In synergy with wide-sky optical/NIR surveys such as Euclid, LSST, and Roman, this facility would provide the missing multi-scale and multiphase submillimeter view needed to uncover how the stellar, gaseous, and dark-matter components assemble in protoclusters, completing our view of early structure formation.
http://arxiv.org/abs/2512.13911 (pdf) (30)
The Large and Small Magellanic Clouds (LMC, SMC) are nearby dwarf galaxies whose proximity uniquely enables molecular cloud-scale resolution observations across the entire Magellanic system, a capability unmatched in any other external galaxy. Their low metallicities resemble conditions near the peak of cosmic star formation, allowing resolved studies of interstellar medium (ISM) phases, molecular cloud lifecycles, and feedback processes that regulate galaxy evolution. Comprehensive, wide-field spectroscopic mapping of CO and its isotopologues in different transitions, complemented by [CI] observations, and combined with already existing HI and HII surveys, will calibrate star-formation laws and gas-phase partition under low-metallicity conditions and furnish benchmarks for interpreting high-redshift galaxies and cosmological simulations. This science requires a large-aperture, wide-field submillimeter single dish with multi-pixel spectroscopic capabilities, operated from a high, dry site such as Chajnantor, to deliver fast, sensitive, high-resolution, degree-scale mapping with total-power fidelity. We present the case for a full molecular atlas of the Magellanic system, the enabling facility requirements, and the transformative impact on galaxy evolution studies
http://arxiv.org/abs/2512.13957 (pdf) (29)
We leverage the overmassive black holes (M_BH/M_∗≈0.1) present in a realization of the BRAHMA cosmological hydrodynamic simulation suite to investigate their role in the emission of the unique “little red dot” (LRD) objects identified by the James Webb Space Telescope (JWST). We find that these black holes can produce LRD-like observables when their emission is modeled with a dense gas cloud shrouding the active galactic nucleus (AGN). Between redshifts 5 and 8, we find the number density of LRDs in this simulation to be 2.04 ± 0.32 × 10^-4 Mpc^-3, which is broadly consistent with current estimates for the total LRD population from JWST. Their emission in the rest-frame visible spectrum is dominated by their AGN, which induces the red color indicative of LRDs via a very strong Balmer break. Additionally, the elevated mass of the black holes reduces the temperature of their accretion discs. This shifts the peak of the AGN emission towards longer wavelengths, and increases their brightness in the rest-frame visible spectrum relative to lower mass black holes accreting at the same rate. These simulated LRDs have very minimal dust attenuation (A_V = 0.21 ± 0.12), limiting the amount of dust re-emission that would occur in the infrared, making them very likely to fall below the observed detection limits from observatories like the Atacama Large Millimeter Array (ALMA). In contrast to the BRAHMA box, the ASTRID simulation produces systematically smaller black holes and predicts LRD number densities that are more than two orders of magnitude lower than current measurements. We therefore conclude that the presence of black holes that are overmassive relative to their host galaxy, and enshrouded in dense gas, is necessary for AGN-dominated LRD models to reproduce both the observed properties and abundances of JWST LRD populations.
http://arxiv.org/abs/2512.14695 (pdf) (26)
A recent study of the compact elliptical galaxy NGC 4486B using JWST-NIRSpec IFU kinematics confirmed a supermassive black hole (SMBH) of mass M_BH=3.6±0.7×10^8 ( 8% of the stellar mass). In addition to its double nucleus, the nuclear kinematics show pronounced asymmetries: a velocity-dispersion peak displaced by 6 pc from the galaxy center and a 16 km/s offset in the mean stellar line-of-sight velocity near the SMBH. We examine the origin of the 12 pc double nucleus and these asymmetries and show that the observations favor an SMBH surrounded by an eccentric nuclear disk (END). END formation models require the SMBH to experience a gravitational wave (GW) recoil following a binary SMBH merger. Our orbit-superposition models contain 50% retrograde stars at the edge of the nuclear region, in striking agreement with END-formation simulations. We infer a pre-merger mass ratio q>0.15 and a recoil kick of 340 km/s. Our N-body simulations show that with such a kick, the SMBH returns to the center within 30 Myr. Its flat central core is also consistent with earlier binary black hole scouring. We test two alternative mechanisms-buoyancy-driven oscillations and a pre-merger SMBH binary-but neither reproduces the observed offsets, favoring the GW-kick scenario. Our direct N-body simulations further show that a prograde SMBH binary in a rotating host can stall in a corotation resonance, delaying coalescence. Thus, although NGC 4486B is an old, relaxed galaxy near the Virgo cluster center, its SMBH appears to have merged only recently, making its nucleus a rare nearby laboratory for studying post-merger SMBH dynamics.
http://arxiv.org/abs/2512.14176 (pdf) (26)
The NIKA2 Cosmological Legacy Survey (N2CLS) is a large programme using the NIKA2 dual-band camera on the IRAM 30 m telescope. Its goal is to improve our understanding of the physics of distant Dusty Star Forming Galaxies (DSFGs) by carrying out deep surveys of two fields, GOODS-North and COSMOS. This work is focussed on GOODS-North, which was observed for 78.2 hours, simultaneously at 1.2 and 2 mm, with a field of view of ∼240 arcmin^2. With such a deep integration, we were able to measure, for the first time, the confusion noise limits at the 30 m telescope using the best sampled ∼ 62 arcmin^2 and masking sources with a flux greater than 0.54 or 0.17 mJy at 1.2 or 2 mm, respectively. We found a confusion noise of 139.1^+ 15.9_- 19.2±11.9 μJy/beam at 1.2 mm and 38.6^+ 9.6_- 13.1±3.7 μJy/beam at 2 mm (the first uncertainty is statistical, the second is the cosmic variance). In this region, this corresponds to half the instrumental noise. To derive these estimates, we devised a novel estimator, referred to as the cross variance, which also enabled us to estimate the correlated confusion noise between the two bands. Thus, we obtained a result of 49.6^+ 15.9_- 24.8± 6.4 μJy/beam. These values are consistent with the state of the art Simulated Infrared Dusty Extragalactic Sky (SIDES) model.
http://arxiv.org/abs/2512.14081 (pdf) (26)
The origin of the multiphase gas within the Fermi/eROSITA bubbles is crucial for understanding Galactic Center (GC) feedback. We use HI4PI data to investigate the kinematics and physical properties of high-velocity clouds (HVCs) toward the GC. Our results reveal that the HVCs exhibit a distinct asymmetric distribution, closely associated with the bar-driven tilted dust lanes and the distorted overshooting streams. We propose that powerful nuclear outflows interact with these gas-rich, off-plane structures, striping and entraining cold gas from the outer Galactic regions (R_GC 0.5–1.7 kpc) rather than solely from the region of the central molecular zone (CMZ; R_GC<0.3 kpc). In this scenario, as the Galactic bar drives gas inflows along the dust lanes, nuclear outflows simultaneously break through the CMZ, sweeping up and ablating cold gas from the boundary layer of these pre-existing structures. This process naturally accounts for the observed high turbulence, complex spectral signatures, and anomalous spatial-kinematic gas patterns, as well as multiwavelength asymmetries of the bubbles. The HVCs are accelerated to about 230–340 km/s over a dynamical time of 3–6 Myr. When the multiphase, inhomogeneous composition of the gas is included, the estimated gas outflow rate reaches 1 Msun/yr. This value is comparable to the bar-driven inflow rate, indicating a tightly coupled gas cycle in the inner Galaxy. Our research highlights the critical role of bar-driven gas dynamics and nuclear feedback in the secular evolution of the Milky Way, offering a valuable paradigm for investigating gas cycles in external galaxies.
http://arxiv.org/abs/2501.18559 (pdf) (25)
We investigate the variability of the UV luminosity function (UVLF) at z > 5 using the SPICE suite of cosmological, radiation-hydrodynamic simulations, which include three distinct supernova (SN) feedback models: bursty-sn, smooth-sn, and hyper-sn. The bursty-sn model, driven by intense and episodic SN explosions, produces the highest fluctuations in the star formation rate (SFR). Conversely, the smooth-sn model, characterized by gentler SN feedback, results in minimal SFR variability. The hyper-sn model, featuring a more realistic prescription that incorporates hypernova (HN) explosions, exhibits intermediate variability, closely aligning with the smooth-sn trend at lower redshifts. These fluctuations in SFR significantly affect the M_UV - M_halo relation, a proxy for UVLF variability. Among the models, bursty-sn produces the highest UVLF variability, with a maximum value of 2.5. In contrast, the smooth-sn and hyper-sn models show substantially lower variability, with maximum values of 1.3 and 1.5, respectively. However, in all cases, UVLF variability strongly correlates with host halo mass, with lower-mass halos showing greater variability due to more effective SN feedback in their shallower gravitational wells. The bursty-sn model, though, results in higher amplitudes. Variability decreases in lower mass haloes with decreasing redshift for all feedback models. This study underscores the critical role of SN feedback in shaping the UVLF, and highlights the mass and redshift dependence of its variability, suggesting that UVLF variability may alleviate the bright galaxy tension observed by JWST at high redshifts.
http://arxiv.org/abs/2512.13823 (pdf) (25)
Intermediate-mass black holes (IMBHs), with masses between roughly 10^2 M_⊙ and 10^5 M_⊙, represent a largely uncharted component of the black-hole (BH) population. They are theoretically predicted to form in several early-Universe pathways, including the remnants of massive Population III stars, the runaway collapse of dense stellar clusters, and the direct collapse of metal-poor gas. Establishing whether IMBHs are present in dwarf galaxy satellites of the Milky Way (MW), and with what occupation fraction - i.e. the fraction of galaxies with a certain stellar mass that host a central BH - provides one of the most incisive tests of BH seed formation models. Despite their importance, present dynamical constraints on IMBHs remain weak. Dynamical IMBH masses or upper limits are available for very few such systems, with secure detections in less than ten cases. A next-generation wide-field spectroscopic facility, capable of combining deep multiplexed stellar spectroscopy with high-resolution integral-field observations of galaxy centers, would open access to IMBH masses in the ≤10 M_⊙ regime. Such an advance would make possible - for the first time - a robust measurement of the IMBH occupation fraction in dwarf galaxies. A key scientific requirement for the coming decades is to establish the observational and instrumental capabilities needed to detect or tightly constrain IMBHs in nearby dwarf galaxies, particularly in the ≃10^3-10^4 M_⊙ mass range, and thereby enable a measurement of their occupation fraction. Such a measurement is fundamental for distinguishing between competing scenarios for the formation of BH seeds in the early Universe.
http://arxiv.org/abs/2512.14254 (pdf) (24)
Polluted white dwarfs provide unique constraints on the elemental compositions of planetary bodies. The tidal disruption of accreting bodies is thought to form circumstellar dusty disks, whose emission spectra could offer additional insights into the mineral phases of the accreted solid material. Silicates are detected in the mid-infrared spectra of several disks, but do not fully account for the near-infrared excess in the disks&apos spectra. Conductive materials, such as metallic iron, are potential sources of near-infrared emissivity. We investigate the role of metallic iron within silicate dust in the observed spectra of the white dwarfs G29-38 and GD56. Using thermal emission spectra calculations, we analyze the abundance of metallic iron in the dust and the disk structure parameters that best fit the observed spectra. We find that metallic-iron-bearing dust enhances the near-infrared opacity, thereby providing a better fit to the G29-38 spectrum for various silicate compositions than metallic-iron-free dust. The best-fit metal-to-silicate mixing ratio is approximately unity, and for Mg-rich pyroxenes, this value is also consistent with G29-38&aposs stellar atmospheric composition within 1-σ observational uncertainties. Based on the spectral fitting and compositional consistency, Fe-rich silicates without metallic iron cannot be ruled out. The observed GD56 spectrum also favors iron-bearing dust. However, the large observational uncertainties of GD56&aposs stellar elemental abundances hinder a precise comparison between the stellar and dust iron abundances. Upcoming high-precision JWST observations will provide a larger sample, enabling statistical analysis of the correlation between the iron abundances in the atmospheres and circumstellar dust of polluted white dwarfs.
http://arxiv.org/abs/2512.13954 (pdf) (24)
We construct samples of optically-selected clusters and photometric luminous red galaxies (LRGs) from the Hyper Suprime-Cam Subaru Strategic Program final year dataset covering ∼ 1200 deg^2. The cluster catalogs extend out to the redshift of 1.38 and contain more than 10000 clusters with richness larger than 15, where the richness is defined to be a membership probability weighted number of galaxies above the stellar masses of approximately 10^10.2M_⊙. The total number of probable red cluster member galaxies in these clusters are more than 6× 10^5. Photometric redshifts of the clusters are shown to be precise with the scatter better than ∼ 0.01 for a wide redshift range. We detect stacked weak lensing signals of clusters out to the redshift of 1, and use them to update constraints on the mass-richness relation. Our catalog of about 6 million photometric LRGs extend out to the redshift of 1.25, and have the scatter of the photometric redshift better than ∼ 0.02 for the redshift range between 0.4 and 1.0.
http://arxiv.org/abs/2512.13782 (pdf) (24)
Radiation-driven winds are ubiquitous in massive stars, but in Very Massive Stars (VMSs), mass loss dominates their evolution, chemical yields, and ultimate fate. Theoretical predictions have often relied on extrapolations of O star prescriptions, likely underestimating true VMS mass-loss rates. In the first of a series of papers on VMS wind properties, we investigate a feature predicted by Monte Carlo (MC) simulations: a mass-loss `kink&apos or upturn where the single-scattering limit is breached and winds transition from optically thin to optically thick. We calculate hydrodynamically consistent non-LTE atmosphere models using the PoWR^HD code, with a grid spanning 40-135M_⊙ and 12-50 kK at fixed log(L_⋆/L_⊙) = 6.0 and solar-like metallicity with Z=0.02. Our models confirm the existence of the kink, where the wind optical depth crosses unity and spectral morphology shifts from O star to WNh types. The predicted location of the kink coincides with the transition stars in the Galactic Arches cluster and reproduces the model-independent transition mass-loss rate of log(Ṁ_trans) ≈ -5.16 from Vink Gräfener (2012). For the first time, we locate the kink at Γ_e≈ 0.43 (M_⋆≈ 60M_⊙) without relying on uncertain stellar masses. Above the kink, mass-loss rates scale much more steeply with decreasing mass (slope 10), in qualitative agreement with MC predictions. We additionally identify two bistability jumps in the mass loss driven by Fe ionisation shifts: the first from FeIV to FeIII near 25 kK and the second from FeIII to FeII near 15 kK. Our models thus provide the first comprehensive confirmation of the VMS mass-loss kink while establishing a mass-loss relation with complex mass and temperature dependencies with consequences for stellar evolution, chemical yields, and the black hole mass spectrum.
http://arxiv.org/abs/2512.14238 (pdf) (23)
Observations in the visible and near infrared are transforming our view of the processes affecting galaxy evolution, much of which is dominated by interactions with the large scale environment. Yet a complete picture is missing, as no corresponding high resolution view of the warm/hot intracluster, circumgalactic, and intergalactic media exists over large areas and a comparably broad range of redshifts. Combined with wide-field optical IFU surveys such as CATARSIS, a large diameter sub-mm telescope with a degree-scale field of view would enable a joint view of galaxy dynamics and gas thermodynamics, transforming our understanding of environmental processes.
http://arxiv.org/abs/2512.13827 (pdf) (21)
Unveiling the fundamental nature of non-baryonic dark matter (DM) has profound implications for our understanding of the Universe and of the physical laws that govern it. Its manifestation as an additional source of matter necessary to explain astrophysical and cosmological observations indicates either a breakdown of General Relativity or that the current Standard Model of Particle Physics is incomplete. In the standard Cold DM (CDM) paradigm, DM consists of collisionless non-relativistic particles with negligible non-gravitational interactions. This simple hypothesis is very successful on large and intermediate scales, but faces challenges on small galactic scales. Local Group (LG) dwarf galaxies can play a fundamental role to elucidate whether these challenges stem from poorly understood fundamental baryonic processes or instead indicate that alternative DM scenarios need to be considered. In particular, a systematic determination of their DM halo properties as a function of stellar mass and star formation histories (SFH) will provide crucial observational benchmarks for models to deal with the trickiest issue that prevents us from advancing in our understanding of DM nature, i.e. the impact of baryonic processes in altering the properties of the inner regions of DM haloes. Such systematic study would require assembling accurate l.o.s. velocities (and metallicities) for several thousands of stars per galaxy, for an heterogeneous sample of target galaxies, spanning order of magnitudes in stellar mass and covering distances from about 100 kpc to more than 1 Mpc. This calls for both multi-objects spectrographs on 12m-class telescopes with fov of a few deg2 and a multiplex power in the several 1000s with the capability of providing dense sampling of the innermost regions, as well as for wide-area multi-objects spectrographs with fov of several arcmin2 on 30-40m class telescopes.
http://arxiv.org/abs/2512.13818 (pdf) (19)
Baryons heat to temperatures above >10^5 K as they accrete onto massive overdensities – galaxies, groups, clusters, and filaments – where they ionize and become optically transparent. Deep mm-wave observations such as those with ALMA have begun to probe a handful (∼4) of massive systems at z∼2-4, while low-resolution mm-wave surveys have detected thousands of objects at arcminute resolution out to z≈2. To truly advance the field of the evolution of large-scale structures, mapping the warm/hot distribution of ionized gas out to the redshift of their formation, the ESO community requires a large-aperture single-dish (sub-)mm telescope. This will need to provide several orders of magnitude higher mapping speeds than currently available while preserving the few arcsecond resolution required for imaging the gas and removing contaminating radio and dusty thermal signals across the full (sub-)mm wavelength range.
http://arxiv.org/abs/2512.13789 (pdf) (19)
We present the chemical abundances of Fe, alpha- and neutron-capture elements in 12 metal-poor Small Magellanic Cloud (SMC) giant stars, observed with the high-resolution spectrographs UVES/VLT and MIKE/Magellan. These stars have [Fe/H] between -2.3 and -1.4 dex, 10 of them with [Fe/H]<-1.8 dex. According to theoretical age-metallicity relations for this galaxy, these stars formed in the first Gyr of life of the SMC and represent the oldest SMC stars known so far. [alpha/Fe] abundance ratios are enhanced but at a lower level than MW metal-poor stars, as expected according to the slow star formation rate of the SMC. The sample exhibits a large star-to-star scatter in all the neutron-capture elements. The two r-process elements measured in this work (Eu and Sm) have abundance ratios from solar up to +1 dex, three of them with [Eu/Fe]>+0.7 dex and labeled as r-II stars. This [r/Fe] distribution indicates that the r-process in the SMC can be extremely efficient but is still largely affected by the stochastic nature of the main sites of production and the inefficient gas mixing in the early SMC evolution. A similar scatter is observable also for the s-process elements (Y, Ba, La, Ce, Nd), with the stars richest in Eu also being rich in these s-elements. Also, all the stars exhibit subsolar [s/Eu] abundance ratios. At the metallicities of these stars, the production of neutron-capture elements is driven by r-process, because the low-mass AGB stars have not yet evolved and left their s-process signature in the interstellar medium. We also present stochastic chemical evolution models tailored for the SMC that confirm this scenario.
http://arxiv.org/abs/2512.14636 (pdf) (18)
We investigate how the dynamical state of galaxy clusters influences their galaxy populations and mass distributions. Using photometrically selected clusters from the DESI Legacy Imaging Survey cross-matched with the UNIONS galaxy shear catalogue, we classify clusters as evolved or evolving based on their rest-frame r-band magnitude gaps and stellar mass ratios between the brightest cluster galaxies (BCGs) and bright satellites. We measure the stellar mass functions, weak-lensing profiles, and radial number density and red-fraction profiles of stacked clusters in both subsamples. Evolved clusters exhibit more concentrated lensing profiles, bimodal stellar mass functions dominated by massive BCGs, and a deficit of intermediate-mass satellites, while evolving clusters show flatter central lensing signals and an excess of massive satellites. Applying the same selection to IllustrisTNG clusters reproduces these trends and links the observed differences to distinct mass accretion histories. These results demonstrate the close link between cluster galaxy populations and the overall dynamical state of their underlying dark matter halo.
http://arxiv.org/abs/2512.14441 (pdf) (17)
Absorption spectroscopy toward high-redshift quasars provides strong constraints on the putative variation of fundamental constants of physics on cosmological time scales. The submillimeter ground-state transitions of methylidyne (CH) and water (H2O), both molecules widespread and coeval in the interstellar medium, provide a sensitive test for variations of alpha, the fine structure constant, and mu, the proton-to-electron mass ratio, taking advantage of the unmatched spectral resolution and frequency reliability of radio techniques. We used ALMA simultaneous observations of the two species to constrain any velocity offset between their absorption profiles toward the radio-bright lensed quasars PKS1830-211 (z_abs=0.88582) and B0218+357 (z_abs=0.68466). Our observational setup minimizes instrumental errors and known sources of systematics, such as time variability of the absorption profile and frequency-dependent morphology of the background quasar. The excellent correlation between CH and H2O opacities, the large number of individual narrow velocity components, and the number of independent spectra obtained due to the intrinsic time variability of the absorption profiles ensure that even the chemical segregation bias is minimized. We obtained bulk velocity shifts delta_v = -0.048 pm 0.028 km/s and -0.13 pm 0.14 km/s (1 sigma confidence level) between CH and H2O in the direction of PKS1830-211(NE) and B0218+357(SW), respectively. These values convert into the 3sigma upper limits |Delta_alpha/alpha| < 0.55 ppm and 1.5 ppm, respectively, taking into account the independent upper limits on |Delta_mu/mu| previously obtained for these systems. These constraints on |Delta_alpha/alpha|, at look-back times of about half the present age of the Universe, are two to four times deeper than previous constraints on any other single high-z system.
http://arxiv.org/abs/2512.13894 (pdf) (17)
Magnification bias offers a uniquely powerful and independent route to cosmological information. As a gravitational-lensing observable, it probes galaxy-matter correlations without relying on galaxy shapes, PSF modelling, or intrinsic-alignment corrections. Its sensitivity spans both geometry and growth: magnification bias simultaneously responds to the matter density, the amplitude of structure, and the redshift evolution of dark energy (DE) below z ≤ 1. Importantly, its parameter degeneracy directions differ from those of shear, Baryon Acoustic Oscillations (BAO), and Cosmic Microwave Background (CMB) data, making it a complementary and consistency-check probe with substantial diagnostic value for the next decade of precision cosmology. However, the current potential of magnification bias is restricted by limited sky coverage, catalogue inhomogeneities, and insufficiently precise redshift or number-count characterisation. A next-generation wide-field submillimetre facility like AtLAST – capable of uniform, deep surveys and spectroscopic mapping – would overcome these limitations and transform magnification bias into a competitive, high-precision cosmological tool.
http://arxiv.org/abs/2512.14664 (pdf) (16)
Milky Way analogues (MWAs) are usually selected from structural and kinematic properties, but robust SED-based similarity criteria are limited by heterogeneous photometry and incomplete wavelength coverage. We present a homogeneous, aperture-photometry SED of the Milky Way near-twin NGC 3521 from the ultraviolet to the radio decameter range. Fluxes are measured within a fixed elliptical isophotal aperture using GALEX, SDSS, WISE, Spitzer/MIPS, Herschel/PACS+SPIRE, and VLA data, and supplemented by meter/decameter constraints. We report new observations obtained in Jan-Feb 2022 with the Ukrainian T-shape radio telescope and derive, for the first time, an upper limit in the 24–32 MHz band. The UV-to-decameter SED (27 points) is modelled with CIGALE, including a dedicated low-frequency radio prescription () that accounts for emission and absorption effects. Using ZTF and NEOWISE data (2014–2025), we detect genuine nuclear variability; optical trends at ∼2^&apos&apos primarily trace the compact nucleus, while NEOWISE variability reflects a mix of nuclear changes and warm-dust emission within the larger aperture. The preferred fit yields M_⋆≃ 6.0×10^10,M_⊙, SFR≃1.65,M_⊙, yr^-1, M_dust≃1.3×10^8,M_⊙, and an effective dust temperature of ∼23 K. The decameter constraint gives S_28, MHz<11.22 Jy, consistent with expectations for a Milky Way-like system placed at 10.7 Mpc. We conclude that an integrated, homogeneous SED, especially below 100 MHz, provides a complementary diagnostic for identifying and validating MWAs and for interpreting how Milky Way properties would appear to an external observer.
http://arxiv.org/abs/2512.14146 (pdf) (16)
We construct mass models of Milky Way dwarf spheroidal galaxies to place constraints on the central black hole (BH) masses they can host. We model the galaxies as a three-component system consisting of the stars, dark matter halo, and a central black hole, using the Osipkov–Merritt–Cuddeford class of anisotropic distribution function. The posterior distribution of black hole mass remains flat toward the low-mass end, indicating that the kinematic data places an upper limit on the black hole mass. Our analysis yields a 95% credible upper limit of log(M_∙/M_⊙) < 6. We combine our results with black hole mass measurements and upper limits from the literature to construct a unified M_∙–σ_relation spanning σ_∼ 10–300,km,s^-1, described by log(M_∙) = 8.32 + 4.08,log(σ_/200,km,s^-1), with an intrinsic scatter of σ_int = 0.55. We compare the inferred limits to models of black hole growth via momentum-driven accretion and stellar capture, which predict black hole masses in the range 10^3–10^4,M_⊙ for the range σ_∼ 6–12,km,s^-1, in close agreement with the M_∙–σ_* relation within the 95% credible upper limits on the black hole masses derived in this work.
http://arxiv.org/abs/2512.13895 (pdf) (14)
Magnification bias, the lensing-induced modification of background source number counts, provides a uniquely powerful probe of the mass density profiles of galaxies and clusters down to kpc scales. Unlike shear-based weak lensing, magnification bias does not rely on galaxy shapes and thus avoids dominant small-scale systematics. Existing studies, however, are limited by sky coverage, positional uncertainty, and insufficiently deep, confusion-limited submillimetre (submm) surveys. A next-generation wide-field, high-throughput submm facility like the proposed 50m-telescope AtLAST is required to unlock this technique&aposs full diagnostic power.
http://arxiv.org/abs/2512.14644 (pdf) (12)
With large numbers of transients discovered by current and future imaging surveys, machine learning is increasingly applied to light curve and host galaxy properties to select events for follow-up. However, finding rare types of transients remains difficult due to extreme class imbalances in training sets, and extracting features from host images is complicated by the presence of bright foreground sources, particularly if the true host is faint or distant. Here we present a data augmentation pipeline for images and light curves that mitigates these issues, and apply this to improve classification of Superluminous Supernovae Type I (SLSNe-I) and Tidal Disruption Events (TDEs) with our existing NEEDLE code. The method uses a Similarity Index to remove image artefacts, and a masking procedure that removes unrelated sources while preserving the transient and its host. This focuses classifier attention on the relevant pixels, and enables arbitrary rotations for class upsampling. We also fit observed multi-band light curves with a two-dimensional Gaussian Process and generate data-driven synthetic samples by resampling and redshifting these models, cross-matching with galaxy images in the same class to produce unique but realistic new examples for training. Models trained with the augmented dataset achieve substantially higher purity: for classifications with a confidence of 0.8 or higher, we achieve 75% (43%) purity and 75% (66%) completeness for SLSNe-I (TDEs).
http://arxiv.org/abs/2512.14203 (pdf) (12)
(Sub-)millimeter spectral lines can be used not only to understand the chemical complexity and enrichment history of an observed portion of our Galaxy, but with spectrally resolved lines, they reveal the physical conditions, dynamics, and even the ionisation state and magnetic field strengths of the gas component of our Galaxy. They are prime tracers of mass assembly and structure formation across scales.
http://arxiv.org/abs/2512.13785 (pdf) (11)
The eponymous BL Lac object BL Lacertae is one of the most well-monitored active galactic nuclei, frequently observed from radio to gamma rays. Its relatively soft γ-ray spectrum peaks near 500 MeV, and since 2020 it has undergone an exceptional series of flaring episodes. The observed emission is well described by synchrotron self-Compton (SSC) models, with negligible contribution from external seed photons. We investigate the physical origin of BL Lacertae&aposs γ-ray temporal and spectral variability using data from the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope, and show that this variability can be explained by a single varying parameter, namely the electrons&apos peak energy, γ_p, under a single-zone SSC scenario with a log-parabolic electron distribution. We use a Markov chain Monte Carlo to estimate the spectral parameters of BL Lacertae over time, selected from an adaptive-binned gamma-ray light curve. We then study the correlation between the inverse Compton peak luminosity, L_IC, and the position of this peak on the SED energy axis, E_p, and compare it with what is expected for a single-zone SSC scenario when only one parameter is free to vary. We find a correlation L_IC = 10^42.33±0.15±0.18_sysE_p^0.98±0.05±0.06_sys, consistent, within the errors, with the linear relation L_IC∝ E_p, expected when γ_p is the only free parameter in the assumed SSC model. This result supports a minimalist SSC scenario in which changes in γ_p dominate the observed temporal and spectral variability of BL Lacertae.
http://arxiv.org/abs/2512.14669 (pdf) (10)
Stellar multiple systems are the norm, not the exception, with many systems undergoing interaction phases during their lifetimes. A subset of these interactions can lead to stellar mergers, where the two components of a stellar binary system come close enough to coalesce into a single star. Accurately modeling stellar mergers requires computationally expensive 3D methods, which are not suited for exploring large parameter spaces as required e.g., by population synthesis studies. In this work, we compare two 1D prescriptions based on the concept of entropy sorting to their 3D counterparts. We employ a basic entropy sorting method (&aposES&apos), which builds the merger remnant by sorting the progenitor stars&apos shells by increasing entropy, and a Python version of the &aposMake Me A Massive Star&apos code (&aposPM&apos), which additionally applies a shock-heating prescription calibrated on SPH simulations of stellar head-on collisions. Comparing to a set of 39 more recent SPH head-on collisions different from the ones used for PM calibration, we find that PM reproduces the outcome of these mergers a lot better than ES in terms of thermal and composition structure post-merger. Both 1D methods produce remnants that are rejuvenated more strongly than expected for massive stars, indicating that increased amounts of hydrogen are being mixed into the core. In an effort to further improve PM, we introduce a scaling factor for the shock-heating. We compare 1D models with both down- and up-scaled heating to a 3D MHD 9 + 8 M_⊙ merger of main-sequence stars. Decreasing the shock-heating improves the agreement in terms of the entropy profile, but has only a minor impact on the subsequent stellar evolution of the remnant. We find that 1D methods are able to approximate 3D stellar merger simulations well, and that shock-heating has to be considered to properly predict the post-merger structures.
http://arxiv.org/abs/2512.14348 (pdf) (10)
Magnetic flux ropes (MFRs) are fundamental magnetic structures in solar eruptions, whose formation is generally attributed to (1) the emergence of subsurface flux tubes or (2) flux cancellation driven by photospheric horizontal flows and magnetic reconnection. Both mechanisms can operate simultaneously during active region evolution, making their relative contributions challenging to quantify. Here, we analyze the formation of a flux rope in a MURaM radiative magnetohydrodynamic (RMHD) simulation, which formed and evolved for approximately three hours before an M-class flare. The formation process is quantified by magnetic helicity flux, which drives the non-potential evolution of magnetic field, with its advection and shear terms on the photosphere corresponding to the emergence and photospheric horizontal flows, respectively. Examining the helicity injected into the flux rope through the photosphere, we find both terms increase significantly as the eruption approaches, with the shear term prevailing overall. Height-dependent analysis of helicity flux, together with magnetic field and velocity distributions, further reveals a gradual transition from the shear to the advection term with an increasing altitude, which is driven by magnetic reconnection above the photosphere. Our results provide quantitative evidence that flux cancellation governs flux rope formation, arising naturally from magnetic field reorganization during active region evolution: as flux emergence transports magnetic flux upward, photospheric shearing motions adjust magnetic field and inject helicity into solar atmosphere, and magnetic reconnection ultimately assembles the main body of flux ropes.
http://arxiv.org/abs/2512.14307 (pdf) (10)
The large-scale limits on the relic signals of gravitational radiation complement the bounds coming from the interferometric detectors (in the audio band) and from the pulsar timing arrays (in the nHz range). Within this inclusive perspective the spectral energy density of the gravitons is sharply suppressed in the aHz region even though the high frequency signal can be comparatively much larger both in the kHz and GHz domains. For there are no direct tests on the expansion rate prior to the formation of the light nuclei, a modified postinflationary timeline affects the total number of e-folds and additionally suppresses the tensor to scalar ratio by making the relic signals effectively invisible in the aHz range. The expansion rate prior to nucleosynthesis is further bounded by the evolution of the hypercharge field and the large-scale magnetism also constrains the decelerated expansion rate. The magnetogenesis requirements are compatible with a potentially detectable spectral energy density of the relic gravitons between the MHz and the THz while the tensor to scalar ratio remains suppressed in the aHz region. A maximum of the spectral energy density of the gravitons in the audio domain leads instead to a larger magnetic field when the scale of the gravitational collapse of the protogalaxy (of the order of the Mpc) gets comparable with the Hubble radius before equality. Along a converse viewpoint the results obtained here imply that a long decelerated stage expanding faster than radiation does not affect the high frequency range but reduces the effective number of e-folds by so enhancing the tensor to scalar ratio, possibly beyond its observational limit.
http://arxiv.org/abs/2512.13814 (pdf) (9)
HD100546 represents a particularly interesting target to study dynamical planet-disk interactions as various features have been observed in both the dust and gas that provide direct and indirect evidence for ongoing planet formation. In this work, we aim to characterize the gas kinematics of five molecular CO emission lines observed with ALMA in HD 100546, to reveal deviations from Keplerian rotation as well as substructures in the peak intensity and line width. We fit the molecular intensity channels with the Discminer package to model the line profiles. Aside from fitting the full cube, we also conduct runs where the blue- and redshifted sides are modeled separately to search for possible asymmetries. Our analysis reveals prominent kinematical spiral features in all five lines on large scales of the disk and we reproduce their morphology with both a linear and logarithmic spiral. In 12CO 2-1, spirals are also seen in the peak intensity residuals, the line width residuals exhibit a prominent ring of enhanced line widths around 125-330 au. The models further show, that the emission from the redshifted side may originate from higher disk layers than that from the blueshifted side. The pitch angles of the spirals are consistent with those driven by an embedded companion inside of 50 au and they suggest a dynamical mechanism rather than gravitational instabilities. We further find indications of a companion around 90-150 au, where tentative dips are present in the radial profiles of the integrated intensity of 13CO and C18O and pressure minima are observed in the azimuthal velocities. For the first time, we also detect downward vertical flows in this region, which coincide with the observed dust gap. The asymmetry in the emission heights may be a result of infall from the disk&aposs environment. Another explanation is provided by a warped inner disk, casting a shadow onto one side of the disk.
http://arxiv.org/abs/2512.13804 (pdf) (9)
The internal structure of dark matter haloes encodes their assembly history and offers critical insight into the nature of dark matter and structure formation. Analytical studies and high-resolution simulations have recently predicted the formation of &aposprompt cusps&apos - steep power-law density profiles that emerge rapidly from the monolithic collapse of primordial peaks. Yet, by z=0 most haloes are well described by Navarro-Frenk-White (NFW) density profiles, raising the question of how these early cusps are transformed in a cosmological setting. We address this problem using 64 zoom-in N-body simulations of eight haloes, each resimulated with eight different free-streaming wavenumbers to control the abundance of small-scale structure while keeping the large-scale environment fixed. To mitigate numerical discreteness effects, we classify artificial fragments and genuine subhaloes with a physically motivated procedure based on matching subhaloes to their progenitor peaks. At the population level, we demonstrate that haloes initially form prompt cusps, and their profiles subsequently transition towards the NFW form. Our study reveals three distinct pathways by which prompt cusps evolve: major mergers, accretion of artificial fragments, and interactions with large-scale filaments, each having a distinct impact on the inner density profile. In particular, we show that the original power-law cusp remains visible in the profile of particles associated with the primordial peak even when the total halo profile is already NFW-like. This work highlights the imprint of early collapse on present-day halo structure and provides new insights into the origin of the universality of the NFW profile.
http://arxiv.org/abs/2512.14593 (pdf) (9)
The ESA&aposs PLAnetary Transits and Oscillations of Stars (PLATO) mission aims to detect planets orbiting around dwarfs and subgiant stars with spectral type F5 or later, including M-dwarfs. The PLATO Input Catalogue (PIC) contains all targets available for observation by the nominal science. The latest version, PIC2.1.0.1, focuses on the Southern PLATO field, named LOPS2, selected as the first long observation field, and includes the P4 sample, one of the four target samples outlined in the Science Requirement Document. P4 includes the M-dwarfs with magnitudes V < 16 located within LOPS2. A characterisation of the M-dwarfs in the PIC is essential for assessing their potentiality to host exoplanets, and eventually estimate the hosted planet(s) properties. The purpose of this paper is to describe how we selected the P4 M-dwarf targets, and obtained their fundamental parameters and properties. In this work, we introduce the P4 sample and detail the methodologies adopted for the measurement of their stellar parameters. Based on a statistical analysis of the P4 sample, we assess both the photometric and volume completeness, and classify the stellar populations according to their Galactic spatial-velocity components. The adopted stellar parameters are validated by comparison with independent methods from the literature used to estimate stellar radii. The P4 sample is compliant with the PLATO science requirements. Being magnitude limited, its volume completeness decreases going towards distances larger than 30 pc, where late-type targets are progressively less covered. The observed large spread in the colour-magnitude diagram is likely due to the combination of several effects such as metallicity, age, binarity and activity. The strategy we adopted for deriving stellar parameters provides results consistent with those obtained in the literature with different and independent methods.
http://arxiv.org/abs/2512.13838 (pdf) (9)
Stars spend most of their lifetime on the main sequence (MS), where hydrogen burning establishes the internal chemical structure that governs the subsequent evolution. In massive stars, mass loss through winds and binary interactions can significantly modify this structure during the MS. We present a new MS evolution framework suitable for rapid binary population synthesis, implemented in the COMPAS code. Building on the semi-analytical model of Shikauchi et al. (2025), our framework captures the evolution of the convective core on the MS under arbitrary mass-loss or mass-gain histories, including a treatment for stellar rejuvenation and MS mergers. This new framework yields more massive helium cores at terminal-age MS, more compact radii in stripped MS stars, and systematically higher black hole masses than commonly used prescriptions. By providing a more realistic treatment of MS evolution, this framework improves the physical consistency of massive stars and binary evolution in rapid population synthesis.
http://arxiv.org/abs/2512.14455 (pdf) (8)
We study primordial perturbations in Higgs-R^2 inflation in the presence of a non-minimal kinetic mixing between the Higgs and the scalaron. By performing a full numerical integration of the multifield background and linear perturbations, we identify distinct dynamical regimes controlled by the non-minimal coupling ξ_h. For ξ_h ∼𝒪(0.1), the turning rate of the inflationary trajectory induces transient multifield effects that imprint localized features in the scalar power spectrum. In the limit ξ_h≪ 1, the adiabatic spectrum becomes featureless, but isocurvature modes fail to fully decay, leaving a residual contribution at the end of inflation. We compute CMB angular power spectra TT, TE, EE and show how these regimes yield distinct observational signatures, delineating the parameter space compatible with current bounds. Our results highlight that suppressing spectral features does not guarantee the elimination of isocurvature perturbations, placing non-trivial constraints on Higgs-R^2 multifield realizations and motivating future probes sensitive to primordial isocurvature.
http://arxiv.org/abs/2512.13802 (pdf) (8)
Understanding how mass and angular momentum flow through accretion discs remains a fundamental unsolved problem in astrophysics. Accreting white dwarfs offer an ideal laboratory for addressing this question: their variability occurs on accessible timescales of seconds to minutes, and their optical spectra contain continuum and emission-line components that trace distinct disc regions. Broad-band timing studies have revealed time-lags similar to those observed in X-ray binaries and active galactic nuclei, suggesting propagating fluctuations and possible coupling to an inner hot flow. However, the blending of line and continuum light in broad filters prevents a physical interpretation of these signals. The 2040s will bring an unprecedented number of disc-accreting systems discovered by Rubin-LSST, space-based gravitational-wave observatories, and third-generation ground and space-based detectors. To extract disc physics from these sources, high-cadence optical spectral-timing, simultaneously resolving continuum and individual lines, is essential. Such measurements would directly map how variability propagates through discs, determine how the outer disc responds to changes in the inner flow, and test whether accretion physics is scale-invariant from white dwarfs to supermassive black holes. This white paper outlines the scientific motivation and observational capabilities required to realise this vision. It highlights the opportunity for ESO to enable a transformative new window on accretion physics in the coming decade.
http://arxiv.org/abs/2512.13791 (pdf) (8)
Appearing as secondary higher-velocity absorption components, high-velocity features (HVFs) have been observed in several absorption lines in many Type Ia supernovae (SNe Ia). The frequency and ubiquity of these components in silicon and calcium features specifically indicates that the mechanism through which they form must be common occurrence among the majority of SNe Ia. Here we present modelling of the HVF evolution in a sample of six well observed SNe Ia with the radiative-transfer code tardis. A base model is derived for each of the SNe to reproduce the photospheric velocity components, followed by a grid of simulations with Gaussian enhancements to the density profile at high velocities. We train a set of neural networks to emulate the impact of these density enhancements upon the simulated silicon line profile. These networks are subsequently used within a Markov-Chain Monte Carlo (MCMC) framework to infer the density enhancement parameters that most closely reproduce the HVF evolution. While we obtain good matches for the silicon profile, we find that a single density enhancement alone cannot simultaneously produce the observed silicon and calcium HVF evolution. Our findings indicate that neither the delayed-detonation mechanism, nor the double-detonation mechanism can produce these HVFs, suggesting that something may be missing from the models.
http://arxiv.org/abs/2512.14625 (pdf) (8)
We present a harmonic-space estimator for the cross-correlation between the kinematic Sunyaev-Zel&aposdovich effect and the reconstructed galaxy momentum field that offers several practical advantages over the traditional stacking approach. The estimator is easy to deploy using relatively modest computational resources and recovers all information available in the galaxy-kSZ cross-correlation. In particular, by using well-understood power spectrum estimation techniques, its statistical uncertainties, including potential correlated uncertainties with other large-scale structure observables, can be easily and accurately estimated. Moreover, standard kSZ stacking measurements can be reconstructed exactly from the estimator at a lower computational cost, employing harmonic-space, catalog-level techniques to recover all small-scale information.
http://arxiv.org/abs/2512.14479 (pdf) (8)
Blazars, a subclass of radio-loud active galactic nuclei with relativistic jets aligned close to our line of sight, emit highly variable non-thermal radiation across the electromagnetic spectrum. The physical origin of their emission and the blazar sequence remain open questions. We present a self-consistent two-zone leptonic model in which relativistic electrons accelerate in a compact region, losing energy via synchrotron and inverse Compton processes, and escape into a larger zone permeated by an external photon field associated with magnetohydrodynamic winds from the accretion disk. By varying only the mass accretion rate onto the central black hole, the model naturally reproduces the blazar sequence, including Compton Dominance, γ-ray spectral indices, and the positions of synchrotron and inverse Compton peaks, while variations in secondary parameters account for the observed spread in the data. Flat Spectrum Radio Quasars exhibit strong external Compton emission from the extended zone, whereas BL Lac objects are dominated by synchrotron and synchrotron self-Compton emission from the compact acceleration region. This framework highlights the key role of accretion rate and spatially structured emission zones in shaping blazar spectra and provides a unified interpretation of their diverse phenomenology.
http://arxiv.org/abs/2512.14198 (pdf) (8)
Context. Numerous classical novae have been observed to emit γ-rays (E > 100 MeV) detected by the Fermi-LAT. The prevailing hypothesis attributes this emission to the interaction of accelerated particles within shocks in the nova ejecta. However, the lack of non-thermal X-ray detection coincident with the γ-rays remains a challenge to this theory. Methods. We performed similar analyses of the Fermi-LAT data for both novae to determine the duration, localization, and spectral properties of the γ-ray emission. These results were compared with optical data from the AAVSO database and X-ray observations from NuSTAR, available for V1723 Sco 2024 only, to infer the nature of the accelerated particles. Finally, we used a physical emission model to extract key parameters related to particle acceleration. Results. V1723 Sco 2024 was found to be a very bright γ-ray source with an emission duration of 15 days allowing us to constrain the spectral index and the total energy of accelerated protons. Despite early NuSTAR observations, no non-thermal X-ray emission was detected simultaneously with the γ-rays. However, unexpected γ-ray and thermal hard X-ray emission were observed more than 40 days after the nova outburst, suggesting that particle acceleration can occur even several weeks post-eruption. V6598 Sgr 2023, on the other hand, was detected by the Fermi-LAT at a significance level of 4σover just two days, one of the shortest γ-ray emission durations ever recorded, coinciding with a rapid decline in optical brightness. Finally, the high ratio of γ-ray to optical luminosities and γ-ray to X-ray luminosities for both novae, as well as the curvature of the γ-ray spectrum of V1723 Sco below 500 MeV, are all more consistent with the hadronic than the leptonic scenario for γ-ray generation in novae.