Most interesting arXiv papers
Tuesday 17-02-2026
http://arxiv.org/abs/2412.15027 (pdf) (106)
We present JWST/NIRSpec IFU observations of the z=7.152 galaxy system B14-65666, as part of the GA-NIFS survey. Line and continuum emission in this massive system (log10(M*/Msol)=9.8+/-0.2) is resolved into two strong cores surrounded by diffuse emission, as seen in recent JWST/NIRCam imaging. Our dataset contains detections of [OII]3726,3729, [NeIII]3869,3968, Balmer lines, [OIII]4959,5007, HeI5875, and weak [OIII]4363. Each spectrum is fit with a model that consistently incorporates interstellar medium conditions (i.e., electron temperature, Te, electron density, ne, and colour excess, E(B-V)). The resulting line fluxes are used to constrain the gas-phase metallicity (Zg 0.2-0.3 solar) and HBeta-based SFR for each region. Common line ratio diagrams (O32-R23, R3-R2, Ne3O2-R23) reveal that each line-emitting region lies at the intersection of low- and high-redshift galaxies, suggesting low ionisation and higher metallicity compared to the predominantly lower-mass galaxies studied with the JWST/NIRSpec IFU so far at z>5.5. Spaxel-by-spaxel fits reveal evidence for both narrow (FWHM<400 km s-1) and broad (FWHM>500 km s-1) line emission, the latter of which likely represents tidal interaction or outflows. Comparison to ALMA [CII]158um and [OIII]88um data shows a similar velocity structure, and we explore optical-far-infrared diagnostics. The two core galaxies both lie on the mass-metallicity relation at z>4, but show contrasting properties (e.g., M*, Zg), suggesting distinct evolutionary pathways. Combining the NIRSpec IFU and ALMA datasets, our analysis opens new windows into the merging system B14-65666.
http://arxiv.org/abs/2602.12325 (pdf) (95)
We present a panchromatic optical-mm characterization of GNz7q, a recently identified X-ray weak, rapidly growing red quasar embedded within a dusty starburst galaxy at z=7.1899, using the full suite of JWST/NIRCam, NIRSpec, MIRI, and archival NOEMA observations. Our deep NIRSpec/G395M spectroscopy reveals unambiguous broad Balmer emission (FWHM =2221±20kms^-1), confirming a super-Eddington accreting black hole (λ_Edd=2.7±0.4) with a mass of log(M_BH/M_⊙)=7.55±0.34, using accretion-rate corrected BH mass estimators. After subtracting the point source, we robustly detect stellar emission from the host galaxy across multiple NIRCam and MIRI filters. Out joint morphological-spectral analysis yields a stellar mass of log (M_*/M_⊙)=10.5±0.4 and an intense star formation rate of SFR=330±97 M_⊙ yr^-1, confirming the host as a massive, dusty starburst galaxy. We find that GNz7q lies on the local M_BH-M_* relation (M_BH/M_*≃ 0.001) and is well positioned to evolve into the locus of massive SDSS quasars with log (M_BH/M_⊙)≈ 9 and M_*≈ 10^11 M_⊙ at z∼ 6, owing to its remarkably rapid growth in both the black hole and its host galaxy. This stands in stark contrast to many recently reported JWST AGN populations at similar redshifts, including the little red dots (LRDs), whose weak or undetected star formation makes it difficult for them to grow into the massive galaxies hosting SDSS-like quasars. These results suggest that GNz7q marks as a rare, pivotal phase of early BH-galaxy co-eolution, plausibly providing a crucial direct pathway to the supermassive black hole systems within the first billion years of the Universe.
http://arxiv.org/abs/2602.12548 (pdf) (63)
We present an analysis of medium/high-resolution JWST/NIRSpec spectra for 11 LRDs at z ∼ 4. By decomposing the broad and narrow components of the Balmer emission lines, we investigate the connection between line emission and UV/optical continua for the LRD population. We find that the broad Hα luminosity strongly correlates with the optical continuum (but not with the UV), indicating a common AGN origin for both. In contrast, the [O III] line strength is correlated with the UV continuum rather than the optical. Using the width and luminosity of the broad Hα line, we estimate central black hole masses of 10^6-10^8 M_⊙ accreting at high Eddington ratios, consistent with an early (λ_Edd∼ 0.6), rapid-growth phase of AGN evolution. Assuming a constant mass accretion rate in the framework of slim-disk models, we infer growth timescales of ∼ 10^5-10^7 yr, and suggest LRDs may evolve into narrow-line Seyfert 1 galaxies. Upper limits from our spectra indicate that LRDs exhibit intrinsically weak optical Fe II emission compared to typical AGN. To simultaneously account for the inferred broad-line region size and observed luminosity, we propose a "Clumpy Envelope" model in which the optical emission arises from an extended, clumpy gas with a characteristic radius of tens of light-days. The diversity in observed optical continuum shapes can be explained by radial temperature gradients and self-absorption effects within this structure. Our results demonstrate the power of JWST high-resolution spectroscopy in probing the central engines and physical nature of the LRD population.
http://arxiv.org/abs/2602.12885 (pdf) (62)
Wolf-Rayet (WR) stars have recently attracted attention as possible drivers of early chemical enrichment, including the production of fluorine, whose nucleosynthetic origin remains debated. To test the contribution of massive stars to fluorine production in the early Universe, we conducted Atacama Large Millimeter/submillimeter Array Band 5 spectroscopy of the HF(1-0) absorption line toward a dusty star-forming galaxy at z=6.024. This galaxy has a known gas-phase metallicity and is too young for low-mass AGB stars to have contributed significantly, providing a clean environment to isolate massive-star yields. We do not detect significant HF absorption (∼2σ) and derive a conservative 5σ upper limit of N_HF/N_H_2 < 2.2×10^-9. This limit is about an order of magnitude below typical local measurements, indicating inefficient fluorine enrichment ∼0.9 Gyr after the Big Bang. Comparison with chemical evolution models shows that our constraint is consistent with scenarios without WR yields at this epoch. Expanding the sample of HF absorption measurements in high-redshift galaxies with well-characterized metallicities will be crucial for tracing the onset of WR enrichment and fluorine production across cosmic time.
http://arxiv.org/abs/2602.13389 (pdf) (58)
The Lyman α (Lyα) line from high-redshift galaxies is a powerful probe of the Epoch of Reionization (EoR). Neutral hydrogen in the intergalactic medium (IGM) can significantly attenuate the emergent Lyα line, even in the damping wing of the cross-section. However, interpreting this damping wing imprint relies on our prior knowledge of the spectrum that escapes from the galaxy and its environs into the IGM. This emergent spectrum is highly sensitive to the composition and geometry of the interstellar and circumgalactic media, and so exhibits a large galaxy to galaxy scatter. Characterizing this scatter is further complicated by non-trivial selection effects introduced by observational surveys. Here we build a flexible, empirical model for the emergent Lyα spectra. Our model characterizes the emergent Lyα luminosity, the velocity offset of the Lyα line with respect to the systemic redshift, and the Hα luminosity, with multivariate probability distributions conditioned on the UV magnitude. We constrain these distributions using z∼5-6 galaxy observations with VLT MUSE and JWST NIRCam, forward-modeling observational selection functions together with galaxy parameters. Our model results in Lyα equivalent width distributions that are a better match to (independent) Subaru observations than previous empirical models. The extended distributions of Lyα equivalent widths and velocity offsets we obtain could facilitate Lyα transmission during the early stages of the EoR. We also illustrate how our model can be used to identify GN-z11-like outliers, potentially originating from merging systems. We publish fitting functions and make our model publicly available.
http://arxiv.org/abs/2602.13578 (pdf) (51)
We present a systematic study of the internal mass structure of early-type galaxies (ETGs) based on 106 galaxy-scale strong gravitational lenses with background quasars, all having spectroscopic redshifts. From this parent sample, we select 24 systems with high-quality ancillary data for joint analysis of strong lensing geometry and stellar kinematics. A key contribution is the derivation of new single-aperture stellar velocity dispersions for 11 lens galaxies via an iterative spectroscopic fitting procedure that mitigates quasar contamination, providing previously unavailable data. We model the total mass-density profile as a power law, ρ∝ r^-γ, and parameterise its logarithmic slope as γ= γ_0 + γ_z · z_l + γ_s ·log, where z_l is the lens redshift and the surface mass density. Within a flat ΛCDM framework and using DESI BAO measurements as a prior, we constrain the parameters via Monte Carlo nested sampling to γ_0 = 1.62^+0.11_-0.12, γ_z = -0.35^+0.08_-0.09, and γ_s = 0.37^+0.08_-0.07 (68
http://arxiv.org/abs/2602.13007 (pdf) (50)
The unprecedented rest-frame UV and optical coverage provided by JWST enables simultaneous constraints on the electron density (n_e) and temperature (T_e) of ionized gas in galaxies at z>5. We present a self-consistent direct method based on multiple OIII]1661,66) and [OIII] (λ4363, and λ5007) transitions to characterize the physical conditions of the high-ionization zone. This new approach is insensitive to a wide range of n_e due to the high critical densities of the OIII] and [OIII] transitions. Applying this technique to six galaxies at z=5-9, we find electron densities up to n_e∼ 3× 10^5 cm^-3 and temperatures of T_e ∼ 20,000 K in systems at z>6. Accounting for these self-consistent densities changes the derived T_e and modifies the inferred metallicities by up to 0.29 dex relative to previous estimates. We discuss the reported N/O overabundances in the high-z galaxies from our sample, which arise entirely from the high N^3+/H^+ values inferred from NIV] lines. We point out that a T_e-stratification, in which the N^3+ zone has a slightly higher T_e than T_e([OIII]), could substantially reduce the inferred N/O. Quantitatively, if T_e(N^3+) were 10
http://arxiv.org/abs/2602.14074 (pdf) (47)
We present a parametric strong-lensing model for the galaxy cluster MACS J1931.8-2635 (z_l = 0.35), accompanying the detection of the spectroscopically confirmed SN Eos at z = 5.13 (Coulter et al. 2026). We identify 10 new multiple-image systems in recent VENUS JWST/NIRCam imaging, so that the model is constrained with a total of 19 robust multiple-image systems – nine of which also have a spectroscopic redshift. For the point-like source corresponding to SN Eos, our model predicts a total of five images, with the observed radial image pair having a similar magnification of μ≃ 25 - 30 and a small time delay of < 5 days, in agreement with their simultaneous observation. According to the model, the other three predicted images arrived earlier, with time delays of 3.7 ± 0.7, 3.5 ± 0.7 and 54.0 ± 10.8 years prior to the two observed images, and with magnifications of 12.9 ± 2.6, 13.0 ± 2.9 and 2.2 ± 0.4, respectively. The absence of detections at the predicted positions, where the host galaxy&aposs images are also visible, confirms the transient nature of the source. SN Eos and its host galaxy are studied in separate articles, and we here focus on the lens model. The final model reaches a very good r.m.s. distance between model and observations of 0.44”. We present the lens-modeling results, including newly identified systems such as a triply imaged, grand-design spiral galaxy candidate at z ≃ 3.65_-0.09^+0.04, and briefly discuss the potential of using high-redshift lensed SNe for cosmography.
http://arxiv.org/abs/2602.13114 (pdf) (39)
Post-starburst (PSB) galaxies, having recently experienced a starburst followed by rapid quenching, are excellent laboratories to probe physical mechanisms that drive starbursts and shutting down of star formation. Integral-field spectroscopy reveals the galaxies&apos spatially-resolved properties, where observed directional patterns can be linked to the galaxies&apos past evolution. We measure the resolved star-formation histories (SFHs), stellar metallicity evolution and dust properties of three local PSBs from the MaNGA survey, down to 0.5" resolution (∼0.3kpc) using a hierarchical Bayesian model. Local parameters were constrained simultaneously with parameters describing spatial trends. We found that all three galaxies first experienced an outer, weaker and slower quenching starburst, followed by a central, stronger and faster quenching starburst that peaked ∼ 1Gyr after the first. The central starbursts induced a significantly stronger rise in stellar metallicity compared to the outer starbursts. These results are consistent with the effects of a recent gas-rich (wet) merger, where the first pericentre passage triggered starbursts in the outer regions, while the later coalescence triggers a stronger centralised starburst. We find non-axisymmetric features in the maps of burst mass fraction and dust attenuation in all galaxies, which could be caused by tidal effects during the recent merger. Comparisons with literature binary merger simulations suggests that the galaxies&apos rapid quenching was driven by gas consumption and the stabilisation against gas gravitational collapse by a growing spheroid, while AGN feedback was not necessarily a primary cause.
http://arxiv.org/abs/2602.14496 (pdf) (38)
The James Webb Space Telescope (JWST) has enabled the discovery of hundreds of supermassive black holes (SMBHs) at redshifts z≳ 4-7. A non-negligible fraction of these SMBHs are hosted in galaxies with BH-to-galaxy mass ratios (M_BH/M_⋆) being excessively larger than that for local SMBHs by ∼ 1-2 dex. The origin of these “overmassive” BHs remains elusive, demanding either a heavy seed formation scenario or rapid growth of seed BHs. Their deviation from local scaling relations challenges our understanding of how SMBHs and their host galaxies coevolve across cosmic time. In this paper, we apply phenomenological modelings for BHs and galaxies to dark matter halo merger histories from N-body simulations to investigate the subsequent evolution of JWST-discovered “overmassive” SMBHs. We find that early evolution of “overmassive” SMBHs is dominated by stunted accretion leading to gradual decreases in M_BH/M_⋆ ratios. In contrast, less massive SMBHs experience super-Eddington accretion during their early evolution, resulting in a slow increase of mass ratios toward M_BH/M_⋆∼ 0.01. Convergence occurs at M_BH∼ 10^8 M_⊙ with M_BH/M_⋆∼ 0.01. At lower redshift, nearly all SMBHs evolve onto local relations, as expected given that our models adopt empirical relations derived from low-redshift observations. This suggests that the global feedback mechanisms regulating the coevolution of M_BH/M_⋆ ratios are implicitly encoded in local relations in terms of star-formation rate distribution, black hole accretion rate distribution and their active (quiescent) fractions.
http://arxiv.org/abs/2602.14197 (pdf) (33)
We report the Fe K relativistic reverberation feature for the first time in the Narrow-line Seyfert 1 galaxy Ton S180. Using a long observation from XMM-Newton we find that the Fe K emission lag peaks at 117±49 s in the lag energy spectrum computed for frequencies (0.3-8.5) × 10^-4 Hz. The lag amplitude drops to 22.85±14.20 s as the frequency increases to (8.5-30) × 10^-4 Hz. The time-averaged spectrum of the source shows a relatively narrow Fe K line at ∼6.4 keV, resulting in black hole spin to be low (a= 0.43_-0.14^+0.10) found from the reflection modelling. We perform general relativistic transfer function modelling of the lag energy spectra individually. This provides an independent timing-based measure of the spin at a=0.30_-0.17^+0.34, and black hole mass M_BH = 0.29_-0.16^+0.01×10^8M_⊙, comparable to the previous measurement, and height of the corona h = 2.59_-0.33^+5.17r_g. Further, we observe that the Fe K lag and the black hole mass fit well in the linear lag-mass relation shown by other Seyfert 1 galaxies.
http://arxiv.org/abs/2512.15652 (pdf) (32)
Over the last 15 years, the Atacama Large Millimeter/submillimeter Array (ALMA) has revolutionized astrophysics by providing unprecedented resolution and sensitivity in observing the cold universe, including the formation of stars, planets, and galaxies. With groundbreaking discoveries ranging from the first detailed images of protoplanetary disks to the kinematics of galaxies in the Epoch of Reionization, ALMA has showcased the vast discovery potential of the (sub-)mm wavelength regime. However, in another 15 years from now–in the 2040s–the science landscape will have changed dramatically as new major observational facilities will have started their operations or have come towards advanced maturity in their scientific outcome (e.g., JWST, Rubin Observatory, ELT, Euclid, Gaia, Plato, Ariel, Roman Space Telescope, SPHEREx, LiteBIRD, LISA, SKA and others). At the same time, ALMA&aposs current Wideband Sensitivity Upgrade will have been in place for 10 years, and ALMA itself will have been operational for 30 years. To fully exploit this era, the community needs a next-generation facility operating at (sub-)mm wavelengths with capabilities far beyond those possible within ALMA&aposs current infrastructure. To this end, ALMA2040 is a community-driven initiative to define the key scientific questions of the 2040s and translate them into a technical vision for a next-generation transformational (sub-)millimeter facility. Our goal with this document is to summarize the current status of the effort, synthesize outcomes from the 2025 workshops, outline next steps toward a reference design concept, and invite broad participation from the global mm/sub-mm community to help shape this future facility. In the following we provide details on the process and scope. We invite everyone who wishes to join the effort and/or contribute to the dedicated White Papers planned for 2026.
http://arxiv.org/abs/2602.14628 (pdf) (32)
To explore the environmental mechanisms causing quenching in nearby star-forming galaxies, we study the variation with local and large-scale environments of a star formation concentration index, C-index ≡log(r_50, Hα/r_50, cont), that traces the spatially-resolved distribution of Hα emission. Our analysis combines (i) GAMA spectroscopic redshift survey data to optically select galaxy groups and reconstruct the cosmic web, (ii) eROSITA data to identify X-ray-emitting groups, and (iii) SAMI Galaxy Survey data to characterise spatially-resolved star formation. We find that galaxies in X-ray+optical groups exhibit the lowest median C-index and the highest fraction of centrally-concentrated star-forming galaxies relative to optical groups and the field (independently of group or stellar mass). Star-forming galaxies in more X-ray luminous groups at fixed dynamical mass show more concentrated star formation. At large scales, nodes show the lowest median C-index and the highest fraction of centrally-concentrated star-forming galaxies relative to filaments and voids, which have similar C-index distributions. C-index correlates most strongly with the distance to the closest node, leaving no significant role for other local or large-scale environment metrics. Finally, regular star-forming galaxies tend to have spins aligned parallel to filaments, consistent with smooth gas accretion, while centrally-concentrated galaxies tend have spins aligned perpendicular to filaments, likely driven by mergers and associated with bulge growth. These results suggest that multi-scale environmental processes, i.e. locally and at large-scale, act to concentrate star formation toward galaxy centres, via gas-related mechanisms in nodes and ram-pressure stripping in X-ray+optical groups.
http://arxiv.org/abs/2602.13396 (pdf) (32)
At low redshifts, hydroxyl megamasers (OHMs) have been shown to trace galaxy mergers, obscured starbursts, high molecular gas densities, and candidate dual supermassive black hole systems. Given this astrophysical utility, exploring these sources at larger cosmological look-back times is therefore of key interest. While previous OHM surveys have been limited to redshifts of z ≲ 0.25, the ability to expand the OHM frontier is significantly enhanced with new high-sensitivity radio facilities such as MeerKAT. In this Letter, we report the discovery of an OHM in the gravitational lens system HATLAS J142935.3-002836 at z = 1.027, the most distant OHM source yet detected. The spectrum has blended 1667 and 1665 MHz emission and exhibits a highly complex profile, with spectral components ranging in widths of <8 km s^-1 to ∼300 km s^-1. The integrated (magnification uncorrected) luminosity of log(L_OH / L_⊙) = 5.51 ± 0.67 makes this the most apparently luminous OHM known to date. In the same wide-band dataset, we have also detected a previously unknown H I absorption line. The signal-to-noise ratio of over 150 with just a 4.7 h observation highlights the potential that MeerKAT and the future Square Kilometre Array mid-frequency array offer to explore the high-redshift OHM universe.
http://arxiv.org/abs/2602.12460 (pdf) (31)
We present a new Pop III + Pop II subgrid framework implemented in the moving-mesh code arepo, designed to study the impact of Pop III feedback on star formation in the early universe. The framework combines primordial non-equilibrium chemistry, metal-line cooling, IMF-sampled stellar evolution with SN feedback, and approximate Lyman-Werner (LW) and ionizing radiation transport. We run a suite of 1c Mpc/h box simulations with different initial conditions and resolutions from z=127 to z=10. The highest gas mass and spatial resolution in the fiducial simulation reach ∼10 M_⊙ and ∼4 pc, respectively. The model successfully reproduces the UV-inferred Pop II star formation rate density (SFRD) from recent JWST observations across all initial conditions, with only minor variation driven by local halo interactions and LW irradiation. We find that the volume filling factor of metal-enriched gas converges to ∼1
http://arxiv.org/abs/2602.12329 (pdf) (31)
Combining deep Hubble Space Telescope (HST) images and extensive data from the Multi-Unit Spectroscopic Explorer, we present new mass models of the cluster MACS J1149.5+2223, strongly lensing the supernova (SN) Refsdal, fully exploiting the source surface-brightness distribution of the SN host for the first time. In detail, we incorporated 77,000 HST pixels, in addition to the known 106 point-like multiple images, in our modeling. We considered four different models to explore the effect of the relative weighting of the point-like multiple image positions and flux distribution of the SN host on the model optimization. When the SN host&aposs extended image is included, we find that the statistical uncertainties of all 34 free model parameters are reduced by factors ranging from one to two orders of magnitude compared to the statistical uncertainty of the point-like only model, irrespective of the adopted different image weights. We quantified the remarkably increased level of precision with which the cluster&aposs total mass and the predicted time delays of the SN Refsdal multiple image positions can be reconstructed. We also show the delensed image of the SN host, a spiral galaxy at zSN = 1.49, in multiple HST bands. In all those applications, we obtain a significant reduction of the statistical uncertainty, which is now below the level of even the small systematic uncertainty on the mass model that could be assessed by the different approaches. These results demonstrate that with extended image models of lensing clusters it is possible to measure the cluster&aposs total mass distribution, the values of the cosmological parameters, and the physical properties of high-redshift sources with an unparalleled precision, making the typically not-quantified systematic uncertainties now crucial.
http://arxiv.org/abs/2602.14882 (pdf) (29)
We present imaging and spectroscopic analyses of Chandra and XMM-Newton observations of ACT-CL J0123.5-0428, one of the most massive, highest redshift galaxy clusters detected within the survey fields of the Atacama Cosmology Telescope. The Chandra data are sufficient to characterize the morphology of this cluster and constrain the geometrically deprojected temperature in 2 spatial bins out to r_2500, revealing a dynamically relaxed system whose temperature drops to kT = 1.8±0.6 keV in the inner ∼40kpc. Within this same inner radius, the surface brightness and density of the ICM is sharply peaked, and the cooling time falls to t_cool=280^+150_-120 Myr. A novel forward-modeling analysis of the XMM data extends imaging and spectroscopic measurements of this system out to r_500, constraining the redshift to z=1.50±0.03, with a mean temperature of kT = 7.3±1.1 keV and an emission-weighted mean metallicity of Z/Z_⊙ = 0.43^+0.46_-0.25. We also utilize the limited optical/IR photometric coverage of the cluster to characterize the properties of the brightest cluster galaxy (BCG), which is coincident with the X-ray peak. Despite the high redshift and strong cool core, the BCG exhibits no signs of recent or ongoing star formation, suggesting AGN feedback has been acting persistently to stem star formation since z∼ 2.5. These measurements identify ACT-CL J0123.5-0428 as the highest redshift, dynamically relaxed, cool core galaxy cluster discovered to date, making it a premier target for future astrophysical and cosmological studies.
http://arxiv.org/abs/2602.13902 (pdf) (28)
Modern studies in astrophysics and cosmology increasingly rely on simulations and cross-survey analyses, yet differences in data generation, instrumentation, calibration, and unmodeled physics introduce distribution mismatches between datasets (“domain shift”). In machine-learning pipelines, this occurs when the joint distribution of inputs and labels differs between the training (source) and application (target) domains, causing source-trained models to underperform on the target. Transfer learning and domain adaptation provide principled ways to mitigate this effect. We study a concrete simulation-to-observation case: semi-supervised domain adaptation (SSDA) to transfer a four-class spectral classifier – high-redshift quasars, low-redshift quasars, galaxies, and stars – from J-PAS mock catalogs based on DESI spectra to real J-PAS observations. Our pipeline pretrains on abundant labeled DESI→J-PAS mocks and adapts to the target domain using a small labeled J-PAS subset. We benchmark SSDA against two baselines: a J-PAS–only supervised model trained with the same target-label budget, and a mocks-only model evaluated on held-out J-PAS data. On this held-out J-PAS data, SSDA achieves a macro-F1 score (balancing precision and recall) of 0.82 and an overall true positive rate of 0.89, compared to 0.79/0.85 for the J-PAS–only baseline and 0.73/0.87 for the mocks-only model. The gains are driven primarily by improved quasar classification, especially in the high-redshift subclass (F1=0.66 vs. 0.55/0.37), yielding better-calibrated candidate lists for spectroscopic targeting (e.g., WEAVE-QSO) and AGN searches. This study shows how modest target supervision enables robust, data-efficient simulation-to-observation transfer when simulations are plentiful but target labels are scarce.
http://arxiv.org/abs/2602.13392 (pdf) (28)
Galaxy evolution is shaped by star formation and stellar feedback at scales unresolved by current high-resolution cosmological simulations. Precise subgrid models are thus necessary, and different approaches have been developed. However, they are degenerate and often primarily calibrated to reproduce stellar masses from observations. To explore these degeneracies, we perform three cosmological zoom-in radiation-hydrodynamics simulations of the same galaxy within a 5×10^11 M_⊙ dark matter halo at z∼1, each with a different subgrid model: mechanical feedback, a combination of mechanical feedback and thermal feedback, and delayed cooling. We calibrate the simulations to match in stellar mass, isolating the effect of the models on the circumgalactic medium (CGM). Our findings demonstrate that despite producing galaxies with comparable stellar masses, the three models lead to distinct feedback modes, resulting in notable variations in the CGM properties. The delayed cooling run is dominated by ejective feedback and exhibits high burstiness, whereas mechanical and the hybrid models primarily feature preventive feedback, respectively acting at the galaxy and halo scales. Delayed cooling reduces the baryon mass to half the universal baryon fraction while mechanical feedback retains most baryons, with the hybrid model standing in between. Delayed cooling also ejects significantly more metals into the CGM than both other models. While for delayed cooling and mechanical feedback metals are almost evenly distributed in the CGM, they are concentrated around satellites in the hybrid model. These discrepancies emphasize the need to design an appropriate subgrid model to understand how stellar feedback regulates galaxy growth.
http://arxiv.org/abs/2602.12839 (pdf) (27)
The main objective is to calibrate the OPT-IR luminosity relation for quasars, focusing on accurate estimations of dusty torus and accretion disk luminosities. We analyzed contaminations related to host galaxies, particularly from polar dust, the interstellar medium, and stellar emission that affect the optical and infrared. We used a sample of nearly 400 quasars with photometrical observations and spectroscopical redshift divided into four redshift bins (0.7-2.4). Full spectral energy distribution (SED) fitting was performed with the CIGALE code, and results were compared with simplified photometric luminosity estimates. The impact of non-active galactic nucleus components and the role of polar dust in the fitting process were assessed. We show that for sources with a disk luminosity above 10^45 [erg/s], the photometric estimates are consistent with SED-based values. While polar dust contributes marginally to luminosity, its presence significantly alters SED fitting, particularly the torus opening angle and cold dust properties. In the optical domain, stellar emission is the dominant contamination. In the infrared, disk emission and cold dust play major roles. We propose two empirical calibrations for the OPT-IR relation. We conclude that the optical band is dominated by the accretion disk component above 10^45 or 10^46 [erg/s] depending on redshift, while IR luminosity is dominated by the dusty torus emission above 1.6 × 10^45 or 2 × 10^46 [erg/s] depending on the redshift. In this high-luminosity regime, simplified photometric methods yield reliable disk and torus luminosity estimates. The aim of the analysis we present is to test the parameter space in order to build a well behaving OPT-IR nonlinear luminosity relation for quasars that could serve as a cosmological probe.
http://arxiv.org/abs/2602.12332 (pdf) (27)
Based on deep high resolution JWST images of AS1063, and after a careful masking of artifacts, extended features in the cluster, and background galaxies (including known lensed ones), we have identified tens of thousands of unresolved point sources in the central region of the galaxy cluster. We extended the identification of these point sources up to 1.18 Mpc from the center of the cluster using data in the second module. Most of these sources are expected to be globular clusters orbiting in the deep potential well of the cluster, but also the surviving compact cores of satellite galaxies. We study the distribution of the globular clusters and compared it with the distribution of mass from a lens model derived from the same JWST data. We find a very tight correlation between the two distributions, but also some differences, including a more concentrated distribution for the globular clusters than for dark matter. We explored the possibility of using the distribution of globular clusters as a proxy for the lensing mass. We find that a simple smoothing kernel can transform the discrete distribution of point sources into a continuous two-dimensional distribution that matches well the lensing convergence. This suggests that globular clusters can be used as tracers of the dark matter distribution in other massive clusters where gravitational lensing constraints are scarce but globular clusters can be detected more easily, for instance in low redshift galaxy clusters.
http://arxiv.org/abs/2602.12324 (pdf) (26)
We investigate the physical origin of star formation suppression in gas-rich early-type galaxies using five high-resolution hydrodynamical idealized galaxy simulations, performed with the moving-mesh code AREPO. These simulations include one Milky Way-like galaxy and four early-type galaxies, of which one early-type galaxy is found to have significantly less star formation despite a substantial molecular gas reservoir. We apply a modified virial theorem to the overdensities in each galaxy to quantify the forces regulating their stability and thus star formation. We find evidence that, in the suppressed galaxy, strong Coriolis forces driven by elevated galactic shear may inhibit gravitational collapse. This is caused by the galaxy&aposs high central compactness, providing a physical mechanism for the suppression of star formation that does not require the removal of molecular gas. In contrast, less compact early-type galaxies host more gravity-dominated clouds and therefore exhibit higher star formation rates. However, we find that this gravitational stability occurs without significantly increasing the classical Toomre-Q parameter, and therefore a new criterion for suppressed star formation may be needed. We also discuss the impact of our choice of overdensity scale and connections to observations of molecular clouds.
http://arxiv.org/abs/2602.14873 (pdf) (25)
We investigate the bar fraction in galaxy pairs from the SDSS to assess how galaxy interactions affect bar structures. Compared to isolated galaxies, close pairs exhibit a significantly reduced bar fraction at projected separations within 25 kpc. This reduction is driven almost entirely by systems showing clear merger or disturbance signatures, indicating that tidal interactions suppress bars. The decline is dominated by a decrease in weak bars, while the fraction of strong bars remains largely unchanged. Bar suppression is primarily associated with major mergers and is strongest in massive host galaxies. A weaker but statistically significant suppression is detected in minor mergers only for massive galaxies with small bulges. In contrast, no significant dependence of bar suppression on the relative orientation between pair members is found. These findings provide observational evidence that tidal perturbations in major mergers play a key role in regulating bar evolution.
http://arxiv.org/abs/2602.13445 (pdf) (24)
Measurements from the Absolute Radiometer for Cosmology, Astrophysics, and Diffuse Emission 2 (ARCADE-2) reveal a strong radio background in the GHz frequency range. Since the cosmological 21-cm signal is measured relative to the background radiation temperature, the presence of a radio excess can significantly alter its characteristics. Previous studies have explored the impact of an inhomogeneous radio background on the global 21-cm signal and 21-cm power spectrum. This non-uniform radio background is also expected to introduce substantial non-Gaussianity. In this work, using the bispectrum, we analyze the non-Gaussianity in the 21-cm signal in the presence of an excess galactic radio background and investigate how line-of-sight radio fluctuations from early galaxies influence its nature. We find that even a moderate enhancement in radio efficiency in early galaxies significantly affects the small-scale 21-cm bispectrum. Furthermore, the delayed heating transition caused by a galactic radio background shifts the sign change in the squeezed-limit bispectrum to lower redshifts (z∼11), providing a potential observational signature for distinguishing different radio background models. These results demonstrate that the 21-cm bispectrum, particularly in the squeezed limit, is highly sensitive to radio background fluctuations, making it a powerful tool for probing high-redshift radio-loud sources and the physics of the early cosmic epoch.
http://arxiv.org/abs/2602.12479 (pdf) (24)
The physical nature and origins of the newly discovered class of Odd Radio Circles (ORCs) remain unclear. We investigate a model whereby ORCs are synchrotron-emitting vortex rings formed by the Richtmyer-Meshkov instability (RMI) when a shock interacts with a low-density fossil radio lobe in the intergalactic medium using 3D magnetohydrodynamic simulations. These rings initially exhibit oscillatory behavior that damps over time. We implement a new method to model Inverse-Compton cooling and synchrotron cooling at high frequencies in a scale-free manner, enabling us to test a wide range of model parameters against the observational constraints. We find that shock strengths of Mach 2-4 are consistent with the data, as expected in accretion, merger-driven, or active galactic nuclei-driven shocks. We find that the initial size of the bubbles required to explain the rings ranges from 140 to 250 kpc, with initial energy in the bubble of order 10^57-10^59 erg, consistent with fossil lobes inflated by moderately powerful radio galaxies. Derived ambient pressures and densities place ORCs in low density environments, such as the outskirts of galaxy groups with ages of order 70-200 Myr. Our synthetic radio maps match the polarization properties of ORC1 and predict a dependency of the tangential magnetic field angle on the aspect ratio of ORCs. A key distinguishing trait of the RMI-driven vortex ring model is that it does not require the ORC to be centered on its host galaxy and is therefore redshift agnostic.
http://arxiv.org/abs/2602.12647 (pdf) (23)
Radio relics in merging galaxy clusters are widely interpreted as synchrotron emission from relativistic electrons accelerated at large-scale shocks. However, the efficiency of diffusive shock acceleration (DSA) is expected to be reduced in the low-Mach-number, weakly turbulent environments characteristic of cluster merger shocks, and recent results suggest that DSA itself may not constitute a viable physical mechanism. In this work, we investigate ballistic surfing acceleration (BSA) as an electrodynamically grounded mechanism for electron energization that does not rely on prescribed diffusion coefficients. We formulate BSA under typical cluster shock conditions and derive the balance between coherent acceleration by the shock convection electric field and radiative losses due to synchrotron and inverse-Compton cooling. This balance determines both the maximum electron energy and the resulting steady-state spectrum. By forward-modeling the associated synchrotron emission and comparing it with integrated radio observations of the Sausage and Toothbrush relics, we find that the observed spectral curvature and high-frequency steepening can be reproduced when only a very small fraction (∼ 10^-9 - 10^-8) of the available BSA acceleration capacity contributes to systematic electron energization. Despite this extremely small efficiency, it is sufficient to accelerate electrons to Lorentz factors γ∼ 10^4 - 10^5 under cluster conditions. These results suggest that radio relics provide a promising astrophysical laboratory for probing coherent acceleration, and that the BSA framework may account for the production of relativistic electrons in cluster shocks.
http://arxiv.org/abs/2602.13384 (pdf) (23)
Dark sectors consisting of atomic constituents (electrons, protons, and photons) offer a well-motivated extension to the Standard Model while providing multiple avenues for phenomenological study. In this work, we explore the impact of conversions between the dark and Standard Model photons in the primordial CMB spectral distortion epoch (10^3 ≲ z ≲ 10^6). These conversions are resonantly enhanced when the induced thermal masses of both photonic species are equal, thus leading to the possibility that sizeable distortions can be produced. To this end, we solve the Boltzmann equation at early times to determine the (irreducible) freeze-in or freeze-out abundance of dark photons. This procedure also allows us to update the limits on generic milli-charged dark sectors using the ACT DR6 bound on the number of effective radiative degrees of freedom (N_eff). By then modeling the evolution of the thermal masses in both sectors, we compute the primordial CMB distortion using the Landau-Zener formalism. We find that when the dark electron and proton are roughly similar in mass (the positronium limit), current spectral distortion data from the COBE/FIRAS instrument is able to rule out novel regions of parameter space. We also forecast bounds from the proposed FOSSIL satellite, finding that spectral distortions can also be used to probe the ultra-low dark electric charge regions of parameter space, which are difficult to investigate by other means.
http://arxiv.org/abs/2602.12343 (pdf) (23)
We present a simple analytic approximation for the covariance between pre-reconstruction galaxy power spectrum measurements and post-reconstruction two-point correlation functions. This cross-covariance is essential for joint analyses that combine full-shape clustering information with baryon acoustic oscillation (BAO) measurements, as commonly performed in modern spectroscopic surveys. Our model builds on the disconnected contribution to the covariance and accounts for the damping of correlations due to the BAO reconstruction process. We validate our analytic prescription against numerical simulations from the Dark Energy Spectroscopic Instrument (DESI), testing both idealized cubic geometries and realistic survey configurations including complex footprints and fiber assignment effects. Despite neglecting survey window functions in the analytic calculation, we find excellent agreement with simulation-based covariances and demonstrate that cosmological parameter constraints are virtually unchanged when using our approximation. Our results show that the pre-post cross-covariance is sufficiently small that even approximate treatments are adequate for cosmological inference, opening a pathway toward fully analytic covariance matrices for next-generation galaxy surveys.
http://arxiv.org/abs/2602.14304 (pdf) (22)
Hubble Space Telescope (HST) observations of 92 galaxies that have a strong showing of I-band Asymptotic Giant Branch (IAGB) stars in their color-magnitude diagrams (CMDs), are used to measure the relative offset between the mean apparent I-band magnitudes of the IAGB population and the corresponding apparent I-band magnitudes of the TRGB as measured in the same frames (and CMDs) of those individual galaxies. This first exploratory, large-sample comparison is independent of any extinction (foreground or internal) that may be shared by these two populations. The marginalized luminosity functions used to determine the modal value of the IAGB population are well fit by a single, symmetric Gaussian. The difference in the two apparent magnitudes (in the sense IAGB minus TRGB) is -0.589 mag, with a combined standard deviation of +/- 0.119 mag. Adopting M_I = -4.05 mag for the TRGB stars, the modal absolute magnitude of the IAGB is then calculated to be M_I(IAGB) = -4.64 +/- 0.12 mag. The ensemble dispersion quoted above gives a standard error on the mean of +/- 0.012 mag (based on the full sample of 92 galaxies). Independently, the three geometry-based zero points for I-band AGB stars are found (in Paper I) to be M_I = -4.49 +/- 0.003 mag in the LMC (4204 stars), M_I = -4.67 +/- 0.008 mag, for the SMC (916 stars) and M_I = -4.78 +/- 0.030 mag for NGC4258 (62 stars), leading to a global zero-point (weighted) average of <M_I> = -4.64 +/- 0.15 mag (stat). The scatter found in the anchors is comparable to the scatter in the field sample discussed here, but the calibration sample is small. The application of this method to galaxies well outside of the Local Group, shows that these standard candles can readily be found and measured out to at least 9 Mpc, using already available archival data
http://arxiv.org/abs/2602.12415 (pdf) (22)
The inner Galaxy is the most complex region of the Milky Way, comprising the bulge, inner thin and thick discs, and inner halo; the formation of the bar transferred gas and stars from the disc inward. Accretion of dwarf galaxies also occurred over the Galaxy&aposs lifetime, merging with the original bulge. In this work, we constrain the metal-rich stars of the earliest spheroidal bulge. To study the oldest bulge stars, distributed in a spheroid, we applied kinematical and dynamical criteria in the metal-rich range [Fe/H] > -0.8. This complements our previous analysis of a symmetric sample with [Fe/H] < -0.8. We derived individual abundances through spectral synthesis for C, N, O, Al, P, S, K, Mn, and Ce using stellar parameters from APOGEE DR17, and compared the results with literature data and chemical-evolution models. The alpha elements Mg, Si, and Ca, and iron-peak elements V, Cr, Co, and Ni follow the expected trends relative to the models. Mn shows secondary behaviour. S and K display significant star-to-star scatter but remain broadly compatible with predictions. Phosphorus and cerium show an excess around [Fe/H] ∼ -0.7, more pronounced than in the metal-poor sample, suggesting a distinctive signature of the earliest bulge population. Diagrams of [Mg/Mn] versus [Al/Fe] and [Ni/Fe] versus [(C+N)/O] indicate an in situ origin for most stars. At super-solar metallicities, a subset shows enhanced K and Mn (possibly S) with low [Ce/Fe], hinting at enrichment linked to the nuclear disc and bar, and tracing a chemically distinct population shaped by the innermost Galaxy.
http://arxiv.org/abs/2602.13394 (pdf) (22)
Galaxy evolution is driven by star formation and stellar feedback on scales unresolved by current high-resolution cosmological simulations, requiring robust subgrid models. However, these models remain degenerate, often calibrated primarily to match observed stellar masses. To explore these degeneracies, we conduct three state-of-the-art cosmological zoom-in simulations of the same galaxy, each incorporating different subgrid models: mechanical feedback, a combination of mechanical and thermal feedback, and delayed cooling. We compare their circumgalactic media (CGM) through quasar absorption sightlines of HI, MgII, CIV, and OVI. Our findings demonstrate that despite producing galaxies with the same stellar masses, the models lead to distinct feedback modes and CGM properties. Column densities and covering fractions serve as effective diagnostics of subgrid models, with all four ions providing strong constraints as they trace diverse gas phases, exhibit complementary spatial distributions, and originate from different mechanisms. Although all simulations bracket observed column density distributions, direct comparisons are limited by scarce detections and significant scatter in absorption strengths. Covering fractions of weak absorbers provides the most robust constraints. All models fail to reproduce HI and MgII covering fractions, and delayed cooling overproduces OVI covering fractions, while the other models underproduce them. The simulation including mechanical feedback reproduces the observed CIV covering fractions well, whereas the other models show slight offsets. We argue that this discrepancy is likely driven by unresolved thermal structures for HI and MgII, and insufficient metals for CIV and OVI, arising from missing physics such as AGNs or cosmic rays.
http://arxiv.org/abs/2602.14747 (pdf) (21)
We investigate whether a non-minimal coupling between dark matter and gravity can influence the internal dynamics of ultra-diffuse galaxies. Within this framework, the gravitational potential is modified by an additional term that captures the interaction between spacetime curvature and the dark matter with a coupling constant determined by a length scale L. Using spherical Jeans modeling, we analyze the kinematic data of three ultra-diffuse galaxies namely: NGC 1052-DF2, NGC 1052-DF4, and Dragonfly 44, which span the observational extremes from dark matter deficient to dark matter dominated systems. For each galaxy we explore several dark matter halo profiles, two orbital anisotropy models, and both with and without Stellar-to-Halo Mass Relation scenarios, and we perform a Bayesian parameter inference. Our results show that across all the considered configurations, the constrained astrophysical parameters are consistent with standard ones from General Relativity. The posterior distributions of L show no preference for non-zero values and result only in upper limits, suggesting that any non-minimal coupling contribution must be small and perturbative on this scales. Future high precision velocity measurements will be essential to determine whether non-minimal coupling effects can become observationally distinguishable in low-acceleration systems.
http://arxiv.org/abs/2602.14087 (pdf) (21)
We combine narrow/medium-band filter photometry from the Southern Photometric Local Universe Survey (S-PLUS) DR4 with ultra broad-band filter photometry from Gaia EDR3 to derive fundamental stellar parameters (T_eff, log g, [Fe/H], ages) and elemental-abundance ratios ([C/Fe] and [α/Fe]) for 5.4 million stars in the Galaxy (4.9 million dwarfs and 0.5 million giants), as well as for over 0.7 million red giant stars in the Large and Small Magellanic Clouds (LMC and SMC). The precisions of the abundance estimates range from 0.05-0.10 dex for metallicity in the relatively metal-rich range ([Fe/H] > -1.0) to 0.10-0.30 dex in the metal-poor regime ([Fe/H] <-1.0), 0.10-0.20 dex for [C/Fe], and 0.05 dex for [α/Fe]. The stellar parameters for LMC and SMC member stars are somewhat less precise than those from the S-PLUS main survey, primarily because of the effect of high reddening. The use of both metallicity- and carbon-sensitive filters provides unbiased measurements of both [Fe/H] and [C/Fe], of particular importance for very low-metallicity ([Fe/H] < -2.0) stars, where carbon enhancement can lead to systematically high estimates of [Fe/H] when only a single metallicity-sensitive filter is employed. Furthermore, multiple narrow-band filters enable metallicity estimates down to [Fe/H] ∼ -4.0 with an accuracy of around 0.3 dex, exceeding the precision typically achieved by low/medium-resolution spectroscopy. This extensive photometric dataset, combined with the other three datasets in this series, will serve as a valuable legacy resource for Milky Way and Magellanic Clouds studies.
http://arxiv.org/abs/2602.13678 (pdf) (21)
We present photometric and spectroscopic observations of luminous red nova (LRN) AT 2025abao, the fourth discovered in M 31. The LRN, associated to the AGB star WNTR23bzdiq, was discovered during the fast rise following the minimum phase. It reached the peak at g=15.1 mag (M_g=-9.5±0.1 mag), and then it settled onto a long-duration plateau in the red bands, lasting 70 days, while it was slowly linearly declining in the blue bands. The object showed similarities at peak with the canonical LRNe V838 Monocerotis, V1309 Scorpii, and with the faint and fast-evolving AT 2019zhd, the third LRN in M31, though the later evolution is different. Spectroscopically, AT 2025abao evolved as a canonical LRN: the early spectra present a blue continuum with narrow Balmer lines in emission; at peak, the spectral continuum has cooled to a yellow colour, with a photospheric temperature of 6000 K. Balmer lines have weakened while absorption lines from metals (Fe I, Fe II, Sc II, Ba II, Ti II) have developed, and in particular broad (FWHM∼700 km/s) from the UV Ca II H K lines. Medium- and high-resolution spectra reveal a counter-P Cygni absorption profile in Hα. Finally, late time spectra show an orange continuum (T∼4000-5000 K), a return in strength of the Balmer lines and the formation of molecular absorption bands. AT 2025abao is the rare case of a LRN with detailed archival information regarding the progenitor system. For the first time, we obtained the spectral energy distribution in the infrared of the precursor of a LRN, which is consistent with that of an M giant/AGB. We propose that the dichotomy of light curve behaviour in LRNe (two peaks vs. plateau) can be explained by the extent and H-richness of the common envelope.
http://arxiv.org/abs/2602.12970 (pdf) (21)
The recent discovery of three off-nuclear tidal disruption events (EP240222a, AT2024tvd, and AT2025abcr) - following the first such source, 3XMM J2150-05 - reveals a small but robust population of off-nuclear, or `wandering&apos, black holes (WBHs) with masses M_∙ > 10^4 M_⊙. Two demographic trends are already apparent: (i) all events occur in massive, early-type parent galaxies with stellar masses 10.8 ≲log_10(M_⋆/M_⊙) ≲ 11.1; and (ii) events at larger halo-centric radii (R_TDE/R_200) are associated with dwarf satellites (M_⋆∼ 10^7 M_⊙), while those closer to halo centers lack detected stellar counterparts. Using results from the cosmological simulation, we show that both trends naturally arise from hierarchical galaxy formation. By combining the simulation with empirical constraints on the local galaxy population, we compute the volumetric density of WBHs, ϕ_WBH(M_⋆), finding that it peaks at log_10(M_⋆/M_⊙)=11.10^+0.05_-0.10 and that more than half of all WBHs in the local Universe reside in galaxies with 10.7 ≲log_10(M_⋆/M_⊙) ≲ 11.2, explaining (i) and predicting its persistence as the sample grows. We further show that ii), i.e., the observed link between detection of stellar counterparts and R_TDE/R_200, is also expected from tidal stripping. These results demonstrate that off-nuclear TDEs are powered by the population of WBHs long predicted by cosmological simulations.
http://arxiv.org/abs/2602.13918 (pdf) (20)
We present the absolute flux and wavelength calibration of the James Webb Space Telescope (JWST) Near-Infrared Camera (NIRCam) Wide Field Slitless Spectroscopy (WFSS) mode. Each of NIRCam&aposs two modules (A and B) provides independent long wavelength (LW) grism spectroscopy over the 2.4-5.0 micron range, with orthogonally oriented R and C grisms. Using commissioning and calibration data from programs 01076, 01536, 01537, 01538, 01479, 01480, 04449, 04498, 06606, and 06628, we have measured the field-dependent geometry and wavelength dispersion of both first and second order spectra across the full detector area. The trace geometry was modeled using two-dimensional third-order polynomials that reproduce the observed spectral positions with an RMS accuracy better than 0.1 pixel. Wavelength calibration, derived from observations of the planetary nebula SMP LMC 58, achieves a precision of 1.6-1.9A for the +1 orders and 0.5A for the +2 orders. Absolute flux calibration, established from observations of the G-type star standard P330E, provides a consistent sensitivity function across all grisms and modules with an absolute flux accuracy of 3
http://arxiv.org/abs/2602.12358 (pdf) (20)
[Abridged] We investigate how the identification of Hickson-like CGs depends on the observer&aposs reference frame, quantifying how frequently the same system would be recognised from different vantage points. Using a mock lightcone built from the Millennium I Simulation plus a semi-analytic model of galaxy formation, we identified 7709 CGs when applying the standard Hickson-like criteria. For each CG, we placed 1000 random observers on a surrounding sphere and reapplied the velocity and compactness requirements to test recoverability. We also examined the variation of population and local isolation. The velocity concordance criterion shows modest sensitivity to the observer&aposs location: 10% of CGs fail for some observers, typically groups with members with high peculiar velocities (>1000 km/s). The compactness requirement is far more fragile, as 44% of CGs are missed by most observers, and these systems are very elongated or are chance alignments in real space. Tightening selection limits reduces this dependence. Lowering the surface brightness threshold to μ≤ 23 mag/arcsec^2 reduces the compactness dependence to 16%, while reducing the velocity limit to ΔV≤ 250 km/s lowers velocity-driven failures to less than 4%. Applying both cuts simultaneously yields up to 84% observer-independent groups, although with a substantially smaller sample. Population and isolation are affected by bright interlopers seen from different directions. While such interlopers are common, they have only a minor effect on the compactness and velocity concordance criteria; however, the local isolation is commonly broken. Observer frame effects, dominated by the compactness criterion, can significantly bias Hickson-like CG samples. However, adjusting surface brightness and velocity difference thresholds allows users to balance the physical reliability according to their specific scientific goals.
http://arxiv.org/abs/2412.13595 (pdf) (19)
The origin of cool-core (CC) and non-cool-core (NCC) dichotomy of galaxy clusters remains uncertain. Previous simulations have found that cluster mergers are effective in destroying CCs but fail to prevent overcooling in cluster cores when radiative cooling is included. Feedback from active galactic nuclei (AGN) is a promising mechanism for balancing cooling in CCs; however, the role of AGN feedback in CC/NCC transitions remains elusive. In this work, we perform three-dimensional binary cluster merger simulations incorporating AGN feedback and radiative cooling, aiming to investigate the heating effects from mergers and AGN feedback on CC destruction. We vary the mass ratio and impact parameter to examine the entropy evolution of different merger scenarios. We find that AGN feedback is essential in regulating the merging clusters, and that CC destruction depends on the merger parameters. Our results suggest three scenarios regarding CC/NCC transitions: (1) CCs are preserved in minor mergers or mergers that do not trigger sufficient heating, in which cases AGN feedback is crucial for preventing the cooling catastrophe; (2) CCs are transformed into NCCs by major mergers during the first core passage, and AGN feedback is subdominant; (3) in major mergers with a large impact parameter, mergers and AGN feedback operate in concert to destroy the CCs.
http://arxiv.org/abs/2602.14895 (pdf) (19)
TiO plays an important role in determining the atmospheric structure of M-type stars and in shaping the visual part of their spectra. We compute synthetic spectra for late-type M giants and identify systematic discrepancies in the wavelength range 5810-5850 A. To investigate the origin of these discrepancies, we analyse experimental TiO absorption cross-sections. We report the detection of a molecular band of the singlet system ^1Π-a^1Δ of the TiO, with an R-head located at 5814.8 A, overlapping the 1-3 band of the C^3Δ-X^3Δ system. The lower state of the band is identified as the a^1Δ v"=0 state, while the upper state is most likely the ^1Π in its ground vibrational state. The empirical band intensity is derived by comparing the relative strengths of neighbouring bands from the C^3Δ-X^3Δ and B^3Π-X^3Δ systems in the experimental cross-section. The band intensity is further validated by synthetic spectrum calculations for the late-type giant 30 Her (M6 III) and comparison with its observed spectrum from the MELCHIORS library. The newly identified band is sufficiently strong to affect the flux distribution in the spectra of cool stars.
http://arxiv.org/abs/2602.14148 (pdf) (19)
The center of the Galaxy is a prominent source in X-rays and gamma-rays. The study of its high-energy (HE) emission is crucial in understanding the physical phenomena taking place in this dense and extreme environment, where the closest supermassive black hole (SMBH) to us, Sgr A*, is lurking nearly invisible, today, in most of the energy spectrum. These phenomena are probably common to other galactic nuclei and may explain the feedback processes between nuclear regions and galaxies, so important for the overall evolution of the Universe. The Galactic center HE emission is very complex and consists of both thermal and non thermal radiation produced by compact and extended sources, surrounded by more diffuse components. All these objects and media are interacting with each other in the narrow and dense Central Molecular Zone (CMZ). Some of them also show relevant extensions towards the Galactic poles, indicating energetic outflows that seem to link the center to the recently observed large Galactic polar structures. In spite of the fundamental advances obtained in the last twenty five years with the most sensitive X-ray and gamma-ray observatories, several questions remain open to investigations. We review here the main observational results and the open issues on the high-energy diagnostics of the Galactic nuclear activity, focusing on processes that take place in the CMZ, and in particular discussing the role of the present and past SMBH activities in powering this region and possibly the whole Galaxy.
http://arxiv.org/abs/2602.14925 (pdf) (18)
The [C/N]-age relation has become a powerful tool for reconstructing the formation history of the Milky Way (MW), providing the largest age sample for field giant stars. However, at metallicities below [Fe/H] < -0.5, stellar surfaces are altered by a poorly understood process known as extra mixing, which modifies [C/N] in a mass- and metallicity-dependent manner. This effect complicates the application of the traditional [C/N]-age relation in metal-poor regimes. Within the MW, constraining the mass dependence of extra mixing is particularly challenging because stars at [Fe/H] < -0.5 are predominantly old and therefore low-mass, leading to strong degeneracies between mass and metallicity. In this work, we explore the potential of the Magellanic Clouds (MCs) to disentangle these effects and constrain extra mixing as a function of age and metallicity. By comparing empirical corrections calibrated in the MW with predictions from thermohaline mixing models, we isolate the mass dependence of extra mixing in the MCs down to [Fe/H] ∼-0.7. We find that the empirical calibration performs well for lower-mass stars (< 1.25 M_⊙), while theoretical models successfully reproduce the observed mass dependence down to ∼ 1.25 M_⊙. We further present the first observational evidence that extra mixing becomes ineffective above ∼ 1.8 M_⊙ at [Fe/H] ∼ -0.7. Our results demonstrate the feasibility of deriving [C/N]-based ages for individual stars in external galaxies. Future observations targeting higher-log g or fainter stars in the MCs will provide stronger constraints on extra-mixing processes and enable the calibration of [C/N]-age relation that can be applied to low-metallicity individual stars in the MW or external galaxies.
http://arxiv.org/abs/2602.14630 (pdf) (18)
Cosmic voids contain higher-order cosmological information and are of interest for astroparticle physics. Finding genuine matter underdensities in sparse galaxy surveys is, however, an underconstrained problem. Traditional void finding algorithms produce deterministic void catalogs, neglecting the probabilistic nature of the problem. We present a method to sample from the stochastic mapping from galaxy catalogs to arbitrary void definitions. Our algorithm uses a deep graph neural network to evolve "test particles" according to a flow-matching objective. We demonstrate the method in a simplified example setting but outline steps to generalize it towards practically usable void finders. Trained on a deterministic teacher, the model performs well but has considerable stochasticity which we interpret as regularization. Cosmological information in the predicted void catalogs outperforms the teacher. On the one hand, our method can cheaply emulate existing void finders with apparently useful regularization. More importantly, it also allows us to find the Bayes-optimal mapping between observed galaxies and any void definition. This includes definitions operating at the level of simulated matter density and velocity fields.
http://arxiv.org/abs/2602.13382 (pdf) (18)
I consider recent observations of ultramassive black holes. These appear to confirm theoretical predictions that the relation between central black hole mass M and spheroid velocity dispersion σ has the same form M ∝ σ^4 in spiral galaxies, elliptical galaxies, and cluster ellipticals, but has differing normalizations. These arise from the need for longer black hole accretion episodes to expel the gas otherwise potentially able to feed the holes in the latter two types of host. In a sample drawn from a mixture of galaxy host types the fitted power of σ will slightly exceed the theoretically-derived value of 4 because of the differing normalizations. The observed hole masses do not currently reach the theoretical maximum values possible for disc accretion, set by the equality of the ISCO and self-gravity radii, probably because the host galaxies have insufficient gas.
http://arxiv.org/abs/2602.12929 (pdf) (18)
The Magellanic Clouds (MCs) offer a unique laboratory for studying galaxy interaction and the evolution of dwarf galaxies. By investigating when and how stars formed, the star formation history (SFH) is a powerful tool to provide constraints for dynamical modeling of the system&aposs past interactions and understand the processes of stripping and triggered star formation in tidally influenced environments. We aim to reconstruct the SFH of the Magellanic Bridge, the gaseous and stellar stream connecting the two Clouds. We used data from the deep optical STEP survey, which covers 54 deg ^2 across the Small Magellanic Cloud (SMC) and the Bridge, reaching stars below the oldest main sequence turnoff at the distance of the MCs. We applied the synthetic color-magnitude diagram (CMD) technique to 14 deg^2 of STEP data. We constructed two libraries of synthetic stellar populations based on the PARSEC-COLIBRI and BaSTI stellar evolutionary models, with metallicities in the range -2.0≤[Fe/H]≤0 across the whole Hubble time. We find a clear peak of recent star formation ∼100 Myr ago in the Magellanic Bridge, which becomes increasingly pronounced toward the SMC. The low metallicity of this population suggests that it formed from gas stripped from the SMC during its most recent close encounter with the LMC. In the eastern part of the Bridge (LMC side), the star formation peaks at earlier times, around 10 Gyr and 2 Gyr ago. We estimate a total stellar mass in the Bridge of (5.1 ± 0.2) × 10^5 M_⊙ and a present-day stellar metallicity of [Fe/H]∼-0.6 dex, close to SMC value.
http://arxiv.org/abs/2602.13754 (pdf) (17)
The presence of dust in supernova remnants (SNRs) is confirmed by extensive infrared data from observatories such as Spitzer, Herschel, and JWST, alongside theoretical models of dust formation. This study explores the existence of dust and ices, particularly water ice via 62 μm in SNRs such as the Crab Nebula and N49, using observational data and preliminary modeling with Cloudy. Observations suggest that water ice may be present in IC 443 and possibly other remnants, though the 63 μm band could also indicate [OI] emission. Theoretical models indicate that water ice could survive under certain conditions in SNRs, with densities and temperatures analyzed. Further observations and refined simulations are needed to confirm these findings.
http://arxiv.org/abs/2602.13470 (pdf) (17)
In this eighth paper of the SOFIA-FORCAST series on Milky Way GHII regions, we present an analysis of the massive star-forming complex W43 Main. We compared our 11 - 37 micron maps with multi-wavelength observations from the near-infrared to radio, and investigated the physical nature of compact sources and dust substructures. We applied SED fitting to constrain properties of the compact infrared objects, and examined the evolutionary states of the extended subregions. We identified 20 compact infrared objects, 16 (80%) of which we classify as massive young stellar objects (MYSOs) or candidate MYSOs. W43 Main resides at the junction of the Scutum spiral arm and the Galactic Bar, a location where enhanced turbulence is anticipated and has been proposed as a potential influence on star-formation activity. Nevertheless, our analysis shows that its Lyman continuum photon production rate, the mass of its most massive MYSO, and its MYSO density are all consistent with the survey-wide median values. We therefore conclude that, despite W43 Main&aposs unique Galactic environment, its present star formation activity appears broadly consistent with that of an average Galactic GHII region.
http://arxiv.org/abs/2602.12677 (pdf) (17)
Spectral-line results from a new cryogenic phased array feed (cryoPAF) on the Murriyang telescope at Parkes are presented. This array offers a significant improvement in field of view, aperture efficiency, bandwidth, chromaticity and survey speed compared with conventional horn-fed receivers. We demonstrate this with measurements of sky calibrators and observations of 21-cm neutral hydrogen (HI) in the LMC and the nearby galaxy NGC 6744. Within 0.3 deg of the optical axis, the ratio of system temperature to dish aperture efficiency is 25 K and the ratio with beam efficiency is 21 K (at 1.4 GHz). For the previously measured T_sys = 17 K, respective efficiency values 0.7 and 0.8 are derived. Our HI observational results are in good agreement with previous results, although detailed comparison with multibeam observations of the LMC suggests that the earlier observations may have missed an extended component of low-column-density gas (8× 10^18 cm^-2). We use the cryoPAF zoom-band and wideband data to make a preliminary investigation of whether the large number of simultaneous beams (72) permits the use of novel data reduction methods to reduce the effects of foreground/background continuum contamination and RFI. We also investigate if these methods can better protect against signal loss for the detection of faint, extended cosmological signals such as HI intensity maps. Using robust higher-order singular value decomposition (SVD) techniques, we find encouraging results for the detection of both compact and extended sources, including challenging conditions with high RFI occupancy and significant sky continuum structure. Examples are shown that demonstrate that 3D SVD techniques offer a significant improvement in noise reduction and signal capture compared with more traditional layered 2D techniques.
http://arxiv.org/abs/2512.24085 (pdf) (15)
This paper presents simulations of the synchrotron self-Compton (SSC) spectrum within the Internal-Collision-induced Magnetic Reconnection and Turbulence (ICMART) model. We investigate how key parameters like the magnetization σ_0 shape the broadband spectral energy distribution by regulating the electron distribution and magnetic field strength. The overall spectrum typically comprises two components: synchrotron radiation peaking at E_p with a low-energy spectral index α between -1 and -1.5, and an SSC component peaking at E_ssc. At high energies, Klein-Nishina suppression causes an exponential cutoff. The flux ratio Y between these components is critical: when Y is small, the SSC peak can be suppressed. Spectral features of the synchrotron component reveal the underlying physical conditions: harder spectra with α∼-1 indicate a large Y parameter and strong KN suppression. We find a positive correlation between Y and σ_0, contrasting with internal shock model predictions. Applied to GRB 221009A, our model suggests σ_0≤20 can reproduce the MeV-TeV observations. This study underscores the value of combined MeV-TeV observations in probing GRB emission mechanisms.
http://arxiv.org/abs/2602.14604 (pdf) (14)
Context. Solid methane (CH_4) is an important molecule in interstellar and planetary environments, serving as a precursor to complex organic compounds and a potential biosignature in exoplanetary studies. Despite its significance, laboratory data on low-temperature phase of methane below 10 K remain limited. Aims. We aim to obtain spectra of methane in binary mixtures at 10 K and compare it to the spectra obtained at 6.7 K. These temperatures correspond to phases II and II* of pure methane and are representative of dark molecular clouds and protostars at early stages. We also aim to test the obtained data applicability to JWST data interpretation. Methods. Laboratory reference spectra were obtained on the ISEAge setup via FTIR spectroscopy in transmission mode. A weighted χ^2 minimization is used for the fitting. Results. We present infrared spectra with corresponding band strengths of pure methane and binary mixtures with methane: CH_4:H_2O,CH_4:CO_2, CH_4:CH_3OH, CH_4:NH_3 at 6.7 K and 10 K showing a 20 Conclusions. Our results reveal distinct spectral features for methane in non-H_2O environments, enabling more accurate interpretation of JWST observations. The dataset of spectra, publicly available on Zenodo, can be used for fitting JWST data.
http://arxiv.org/abs/2602.14192 (pdf) (14)
We present a stellar stream found in images of the nearby, nearly face-on, late-type galaxy, NGC 3938 obtained for the LBT Imaging of Galactic Halos and Tidal Structures (LIGHTS) survey that is thin, has very low mean surface brightness (_g⟩≈ 28.7 mag arcsec^-2 and _r⟩≈ 28.1 mag arcsec^-2), appears to lie nearly on the plane of the sky, and wraps more than half way around a host galaxy that is otherwise apparently isolated. We estimate that the progenitor had a stellar mass of ∼ 3.7× 10^7 M_⊙. Despite an intriguing apparent offset between the centroid of the host galaxy and the apparent center of the stream orbit, we find that we can reproduce the morphology, including this apparent off-centering, with simple models and standard assumptions about the host (thin disk centered within a canonical spherical dark matter halo) and the progenitor satellite orbit. We identify a number of detailed features of the stream, such as changes in curvature and density, that will require more complex models to reproduce. Even this rather simple system provides a rich set of constraints with which to explore the accretion history and gravitational potential of an otherwise unremarkable late-type galaxy. Given the depth of the LIGHTS images, this system is an example of the types of stellar stream that could be found in a majority of nearby giant galaxies with the 10-year stack of Rubin/LSST data.
http://arxiv.org/abs/2602.13385 (pdf) (14)
Globular cluster stellar streams probe galaxy-formation processes and can potentially reveal the distribution of dark matter in galaxies. In many theoretical studies, streams are modeled with particle-spray or direct N-body codes. But particle-spray methods abstract away the internal dynamics of the progenitor by making strong assumptions about the escape physics, while direct N-body is prohibitively expensive for realistic (N>10^5) systems. In this paper, we present the stream-modeling capabilities of KRIOS, a new basis-expansion N-body code for collisional stellar dynamics, that bridges this runtime vs. accuracy gap. We show that KRIOS reproduces NBODY6++GPU cluster models, and their associated streams, more accurately than particle spray in a fraction of the NBODY6++GPU wall-clock time. We then compare KRIOS to various particle-spray methods on 10 orbits similar to known Milky Way streams. The morphology and kinematics of these streams most disagree when the progenitor is tightly bound to the host, as these systems are often subject to stronger tidal forces. Finally, we discuss which elements of the progenitor physics are most important for modeling stellar streams and how these might be incorporated into particle-spray methods.
http://arxiv.org/abs/2602.12792 (pdf) (14)
We present PrestaLine: Gorynych, a comprehensive numerical tool designed to model the dynamical, chemical, and spectral evolution of collapsing molecular cores from the prestellar phase to protostellar accretion. The code integrates three key components: (1) Kamelung, a 1D hydrodynamics module simulating gravitational collapse up to first hydrostatic core (FHSC) formation followed by an accretion of envelope onto a young star; (2) Presta, a chemical evolution module post-processing density and temperature profiles to compute time-dependent molecular abundances; and (3) Uran(IA), a radiative transfer module generating synthetic molecular line spectra. We apply Gorynych to compare low-mass (5 Msun) and high-mass (50 Msun) cores, finding that their dynamical evolution is remarkably similar, with differences primarily in spatial scaling. Chemical evolution reveals sharp abundance changes during the FHSC transition, particularly for CO, H2O, and HCO+, though pre-collapse chemical initialization has minimal impact on most species. Spectral maps of 13CO(2-1), HCO+(3-2), and H2O(110–101) lines show distinct kinematic signatures of infall and depletion, with high-mass cores exhibiting spatially extended emission. Our results highlight Gorynych&aposs utility to couple theoretical collapse models and observations, providing a framework to diagnose core evolution and initial conditions from molecular line data.
http://arxiv.org/abs/2602.14101 (pdf) (13)
Magnetic properties of dust grains play a fundamental role in their alignment with ambient magnetic fields and magnetic dipole emission. In the radiative torque (RAT) paradigm, superparamagnetic inclusions (SPIs) embedded within dust grains are expected to significantly enhance magnetic susceptibility and alignment efficiency. Previous studies have generally assumed SPIs of a single characteristic size. In this work, we develop an effective superparamagnetism model that explicitly accounts for a power-law size distribution of SPIs. We show that the effective zero-frequency susceptibility can be described by the superparamagnetic susceptibility of uniform-sized inclusions evaluated at the critical blocking size, reduced by a factor F_eff∼ 0.1. It exhibits a slight increase with dust temperature T_d, in contrast to the rapid decrease for the case of single-size SPIs. For rotating grains at angular frequency ω, we identify a characteristic resonance size of SPIs that dominates the magnetic response, N_res = (T_d/T_act) ln (ν_0/ω) with T_act activation temperature and ν_0 the characteristic attempt frequency of SPIs. The frequency-dependent effective susceptibility is well described by the maximum susceptibility χ_eff^ max(ω) at N_res, reduced by a factor G_eff∼ 0.1. Unlike models assuming uniform-sized inclusions, we find that the effective susceptibility exhibits a nearly flat spectrum for frequency below ν_0, arising from the progressive activation of larger inclusions at lower frequencies. This effective superparamagnetism model based on the SPI size distrbution has important implications for magnetic grain alignment, dust polarization, and magnetic dipole emission across diverse environments.
http://arxiv.org/abs/2602.13397 (pdf) (13)
We examine the multi-wavelength period-luminosity-color relations for Cepheid variables in the Large and Small magellanic Clouds and the Milky Way. From first-principles stellar physics, the luminosity of a Cepheid is determined by its radius and surface temperature, yielding a fundamental PLC relation whose observational proxies are pulsation period and intrinsic color. Using Cepheids in the Magellanic Clouds, we show that the PLC relation recovers the known geometries and line-of-sight tilts of their disks, confirming its ability to detect true distance-modulus variations that are achromatic and consistent across all filters. Surprisingly, for Milky Way Cepheids with individually determined reddenings and HST and Gaia parallaxes, the residuals from multi-wavelength PL fits are also found to be achromatic, identical in sign and amplitude across all passbands, in this case indicating that parallax errors are the dominant source of scatter. Applying bandpass-averaged corrections to individual Cepheids recovers the theoretically expected wavelength-dependent narrowing of the instability strip, and results in revised parallaxes with a median improvement in precision of roughly a factor of two. In addition, they show no statistically significant correlation with metallicity over the range -0.2 < Fe/H < 0.05 dex. The final extinction- and reddening-corrected PLC relation yields an rms scatter of 0.04 mag, corresponding to 2 percent precision in distance per star. Use of a physically grounded PLC will provide a more robust foundation for the Cepheid-based extragalactic distance scale and the determination of the Hubble constant.
http://arxiv.org/abs/2602.12814 (pdf) (13)
We report the multi-temperature structure of the intracluster medium (ICM) in the Centaurus cluster core observed with XRISM/Resolve. Thanks to its high energy resolution, Resolve enables us to measure fine structures of highly ionized emission lines from Si to Fe and to directly determine the excitation temperature and the ionization temperature from the emission line ratio diagnostics. The observed spectrum in the Centaurus core is well-represented by a double-temperature thermal plasma at collisional ionization equilibrium state rather than an isothermal one. The line ratio diagnostics also support this biphasic temperature structure. Particularly, the observed line ratios show a trend of increasing ionization temperature with atomic mass, while the ionization and excitation temperatures of Fe show nearly the same temperature. The resultant line ratios, which are well-represented by the two temperatures ICM, 1.6 and 3 keV, are also fairly consistent with the expected numbers when assuming the radial single-temperature ICM was projected in the cluster core along the line of sight. Due to the limited low-energy sensitivity of the Resolve with the gate valve closed, we investigated the effect of the cool component using the XMM-Newton/RGS spectrum, but it ultimately did not affect our results. The observed flux ratio between the Fe XXV He alpha resonance and forbidden lines shows an about 20% reduction, suggesting the presence of resonant scattering.
http://arxiv.org/abs/2602.12331 (pdf) (13)
Polycyclic aromatic hydrocarbons (PAHs) are responsible for a variety of near- and mid-infrared spectral features in Galactic and extragalactic sources. A feature at 1.05 μm arising from electronic transitions in PAH cations is predicted by laboratory experiments but has never been observationally confirmed. We conduct a dedicated search for this feature in absorption on a highly-extinguished sight line toward BD+40 4223, a blue supergiant in Cyg OB2, using the TripleSpec spectrograph at Palomar Observatory. We place a 5σ upper limit on the feature strength of Δτ_1.05/A_V < 5.6 ×10^-3, ruling out theoretical estimates with > 10σ significance. We constrain the effective temperature of BD+40 4223 to be log_10(T_eff)=4.41±0.03 and infer that it is veiled by 6.39±0.05 magnitudes of visual extinction, consistent with but more constraining than previous determinations. As dust on the sight line toward BD+40 4223 appears typical of the diffuse interstellar medium, this non-detection challenges existing models of PAH material properties and/or charge distribution.
http://arxiv.org/abs/2602.14825 (pdf) (13)
The angular, harmonic cross-correlation between gravitational wave (GW) events and galaxy catalogues contains rich information on the large-scale structure and the origin of compact binary mergers. In this work, we study how uncertainties in the binary formation channel affect the predicted cross-correlation signal for both current-generation and next-generation networks of detectors. We generate five mock GW catalogues for which we vary the progenitor-to-remnant mass-transfer function and the time-delay probability distribution between progenitor and remnant. We then cross-correlate these catalogues with galaxy samples modelled on the 2MASS Photometric Redshift catalogue (2MPZ) and the Gaia-unWISE quasar catalogue (Quaia). We find that the mass-transfer function has negligible effect on the cross-correlation signal, with differences remaining within redshift uncertainties. In contrast, the time-delay distribution dramatically affects the redshift distribution of the GW events and, with it, the cross-correlation signal, particularly for shallow galaxy catalogues. In particular, current-generation facilities can achieve significant detections only for the longest time delays when cross-correlated with 2MPZ, whilst all cross-correlations with the deeper Quaia catalogue are marginally detectable or consistent with zero. Our exploratory results thus demonstrate that forecasts on cosmological or astrophysical parameters derived from GW-galaxy cross-correlations are, as expected, strongly sensitive to the assumed binary formation history.
http://arxiv.org/abs/2602.13153 (pdf) (13)
Using long-term observations from the Transiting Exoplanet Survey Satellite (TESS) along with spectroscopic observations from the 3.6 m Devasthal Optical Telescope (DOT), we present a comprehensive time-series and spectral analysis of the intermediate polar V515 And. Our analysis reveals that V515 And resides within the period gap, with the detection of its orbital period of 2.73116 h. Additionally, we confirm the earlier findings of the spin and beat periods to be 465.4721 s and 488.6067 s, respectively. The time-resolved timing analysis reveals that V515 And undergoes changes in its accretion geometry, not only between different TESS sectors but also within individual sector observations. The system exhibits a transition in the dominant accretion mode, switching between disc-fed and stream-fed accretion. In the TESS light curve, we identify two successive outburst-like episodes, each persisting for roughly a day and reaching peak luminosities of 2.7×10^33 and 1.9×10^33 erg s^-1. Our analysis suggests that these bursts belong to the recently proposed class of micronovae. The optical spectrum of V515 And is characterised by strong Balmer and He II emission lines and shows an inverse Balmer decrement indicating the magnetic nature of the source.
http://arxiv.org/abs/2602.12531 (pdf) (12)
Traditional interstellar polycyclic aromatic hydrocarbon (PAH) diagnostics rely on empirical band ratios, which often suffer from information loss and sample-selection bias. We introduce a machine learning framework that bypasses these limitations by treating the complete 2.75-20 micron emission spectrum as a high-dimensional fingerprint. Using a Random Forest classifier trained on a dataset of 23,653 spectra, we achieve a robust classification F1-score of about 0.96 across 12 size and charge categories. Our model maintains high performance on synthetic mixtures of unseen molecules. Feature importance analysis reveals that PAH size diagnostics are not universal but highly charge-dependent; while neutral size is traced mainly by C-H stretching modes, sizing ionized species also relies on the morphology of 6-8 micron C-C complexes, with the 12.5 micron feature emerging as a robust cross-charge tracer. This approach provides a robust, data-driven pathway for decoding the physical conditions of the interstellar medium.
http://arxiv.org/abs/2602.14978 (pdf) (12)
More than 150 measurements of the HeII 4686A emission line in spectra of SS433 were obtained during 388 nights in 2020-2025 with the Transient Double-beam Spectrograph on the 2.5 m telescope of Caucasian Mountain Observatory of Sternberg Astronomical Institute. We found that the HeII emission line formation region is not eclipsed and is significantly larger than both the donor star and the photosphere of the supercritical accretion disk. The HeII radial velocity curve was found to be independent of the precessional phase and inconsistent with the photometric curve. These findings suggest that the HeII line does not reflect the orbital motion of the compact object. Therefore, spectroscopic estimates of the masses of the components in SS433 based on the HeII emission line can be unrealistic.
http://arxiv.org/abs/2602.14235 (pdf) (12)
After distinguishing the role of classical and quantum inhomogeneities in cosmological backgrounds, we constrain the initial states of the relic gravitons as soon as the different wavelengths of the spectrum cross the comoving Hubble radius, without any reference to earlier timescales. According to this pragmatic perspective the quantum states with finite energy density at the crossing time consistently affect the two-point functions and the related power spectra. An initial state different from the vacuum turns out to be marginally permitted in the low frequency range (associated with the largest observable wavelengths that crossed the comoving Hubble radius) while the intermediate and high frequency domains of the spectrum are populated by the gravitons produced from the vacuum. The non classical correlations are expected to dominate between the kHz and the THz since in this region gravitons are produced quantum mechanically with a negligible contribution from the initial state.
http://arxiv.org/abs/2602.14963 (pdf) (11)
Spectral distortion anisotropies of the cosmic microwave background (CMB) provide a new probe of the early Universe that can be accessed using traditional CMB imaging techniques. It is possible to compute the creation and evolution of anisotropic signals for various scenarios using the frequency hierarchy method recently developed for CosmoTherm. However, the current treatment is not perfect and some approximations had to be made. Here, we carefully construct a modified form for the evolution equations that has the full equilibrium solutions built into the formulation. We improve the formalism to account for i) additional stimulated scattering effects, ii) kinematic corrections to the thermalization terms, iii) corrections to the standard perturbation variables and iv) direct photon sources. These effect could not be captured with the original formulation of the frequency hierarchy method but are indeed important for cleanly separating real distortions from temperature signals. However, we show that previous results are not altered significantly when compared to the improved formulation presented here. As a new worked example, which could indeed not be treated before, we also illustrate how possible changes in the temperature-redshift relation would create spectral distortion anisotropies in the pre-recombination era. The theoretical methods presented here are also an important step towards being able to consistently predict the CMB spectral distortion anisotropies in photon-dark photon and photon-axion conversion scenarios.
http://arxiv.org/abs/2602.14891 (pdf) (11)
The sensitivities of STIS first-order, medium resolution modes were redetermined from on-orbit observations and CALSPEC models (version 11) of the primary white-dwarf spectrophotometric standard stars G191-B2B, GD 71, and GD 153. The sensitivity of an additional configuration was updated by comparing observations of the secondary standard BD +75 325 with the STIS low-resolution spectrum that has been calibrated consistently with the version 11 models. The procedures used to derive the sensitivities and verify the PHOTTAB reference files prior to their activation in CRDS (on May 1, 2025) are described. Results are presented in graphical form in an extensive appendix. Issues and uncertainties are discussed briefly, along with recommendations for future work.
http://arxiv.org/abs/2602.14640 (pdf) (11)
HESS J0852-463 is a TeV γ-ray source located in the Galactic plane. The region consists of a supernova remnant (SNR, RX J0852.0-4622) with a shell-like morphology, commonly referred to as Vela Junior, and a pulsar denoted PSR J0855-4644. Pulsars are among the most efficient leptonic accelerators in our Galaxy, making this region particularly interesting to study. We utilise the most recent data taken by the High Energy Stereoscopic System (H.E.S.S.), to investigate any γ-ray emission associated with the pulsar in this region, PSR J0855-4644. We applied a full forward folding method on the H.E.S.S. data. Utilising 3D modelling techniques, we evaluated the TeV γ-ray emission towards the various components of this complex system. The distinct energy-dependent morphology observed in our data motivates further investigation of this source. We resolved the emission in the Vela Junior region into various components, several of which correspond to the SNR itself. In particular, we find a new extended component which is coincident with the position of PSR J0855-4644. The spectrum follows a power-law distribution with a best-fit index of ΓE = 1.81 ±0.07stat which differs from the properties of the surrounding γ-ray emission of the Vela Junior SNR. A one-zone leptonic joint fit between the X-rays (from XMM-Newton) and γ-rays (from H.E.S.S.) leads to a lower limit on the magnetic field of 1.6μG and a spectral index of α = 1.88 ±0.01, in line with expectations of pulsar wind nebulae (PWNe). In this paper, we report the first detection of the PWN of PSR J0855-4644 at TeV energies with the H.E.S.S. experiment at a significance of 12.2σ. This is attributed to the advanced techniques of the 3D analysis. Based on the pulsar&aposs characteristics, its PWN is consistent with the known TeV PWNe population in the Galaxy.
http://arxiv.org/abs/2602.14242 (pdf) (11)
In the I-band color-magnitude diagrams (CMD) of resolved nearby galaxies, the reddest asymptotic giant branch (AGB) stars form a previously unremarked-upon, but nevertheless distinct and easily-identified population of high-luminosity stars. Hereafter we refer to this population as being comprised of I-Band AGB (IAGB) stars. Identifying these stars in the Large Magellanic Cloud (LMC), the Small Magellanic Cloud (SMC) and in NGC4258 (for all three of which there are published geometric distances) we find that the marginalized luminosity functions are each well approximated by single-peaked Gaussians, having one-sigma dispersions of +/- 0.22 mag, +/- 0.25 mag and +/- 0.24 mag, respectively. The zero points for the modal I-band absolute magnitudes of IAGB stars are found to be M_I = -4.49 +/- 0.003 mag (stat) in the LMC (4204 stars), M_I = -4.67 +/- 0.008 mag (stat), for the SMC sample (916 stars), and M_I = -4.78 +/- 0.030 mag (stat) for NGC4258 (62 stars). A global average over these three independent calibrations of the IAGB zero point (weighted inversely by squares of their systematic errors) gives <M_I> = -4.65 +/- 0.119 mag (stat) +/- 0.025 (sys). In Paper II we will show the results of applying the IAGB Method to 92 galaxies additional galaxies resolved by HST, reaching out to distances just short of 10 Mpc.
http://arxiv.org/abs/2602.13580 (pdf) (11)
Cold atomic and molecular gas are commonly observed in the winds of both external galaxies and the Milky Way, yet the survival and origin of these cool phases within hot galactic winds is poorly understood. To help gain insight into these problems, we carry out time-dependent chemical modelling of cool clouds in the Milky Way&aposs nuclear wind, which possess unusual molecularto-atomic hydrogen ratios that are inconsistent with both disc values and predictions from chemical equilibrium models. We confirm that CO and Hi emission comparable to that in the observed nuclear wind clouds cannot be produced by gas in chemical equilibrium, but that such conditions can be produced in a molecule-dominated cloud that has had its atomic envelope rapidly removed and has not yet reached a new chemical equilibrium. Clouds in this state harbour large reservoirs of molecular gas and consequently have anomalously large CO-to-H2 conversion factors, suggesting that the masses of the observed clouds may be significantly larger than suggested by earlier analyses assuming disc-like conversions. These findings provide a new framework for interpreting cold gas in galactic winds, providing strong evidence that cold outflows can originate from the galactic disc molecular clouds that survive acceleration into the wind but lose their diffuse atomic envelopes in the process, and suggesting that the Milky Way&aposs nuclear outflow may be more heavily mass-loaded than previously thought.
http://arxiv.org/abs/2602.12355 (pdf) (11)
Context. Tidal disruption events (TDEs) are rare transients that provide important insights into the physics of galactic nuclei. A recently identified feature in their optical light curves is the presence of early bump-like structures (precursors) that appear before the onset of the main flare or during its rise. Aims. We aim to build and study the first sample of precursor TDEs in order to improve our understanding of these features, which could be key to revealing the origin of the optical emission in TDEs. Methods. We compiled all known precursor TDEs from the literature, searched for additional candidates, and analyzed them as a sample. Results. We find that precursor TDEs predominantly fall within the repeating TDE, fast TDE, and TDE in active galactic nucleus (AGN) subclasses. We reveal a positive correlation between the occurrence time of the precursors relative to the main peak and the central black hole mass. Conclusions. We suggest that the precursors appear due to interactions between the incoming stellar debris and the disk or leftover material from an earlier disruption (repeating and fast TDEs) or a stable pre-existing disk (TDEs in AGNs). Precursors are therefore potentially key signatures of repeating partial TDEs in previously quiescent galaxies.
http://arxiv.org/abs/2602.12822 (pdf) (10)
The `Inverse Distance Ladder&apos uses relative-distance measurements with type Ia supernovae (SNe Ia), anchored to an absolute distance scale from Baryon Acoustic Oscillations (BAO) and the cosmic microwave background (CMB), to provide an alternative measurement technique for the local expansion rate H_0. Using SNe Ia from the Dark Energy Survey and BAO measurements from the Dark Energy Spectroscopic Instrument, we show that the inverse distance ladder is unable to explain the Hubble Tension, anchored as it is to the CMB and its value of H_0 = 67.4 ± 0.5 kms^-1 Mpc^-1. To do so, we first show that the suite of systematics considered in cosmology analyses with SNe Ia only move the inferred H_0 by <0.1kms^-1 Mpc^-1, and second, we investigate the scale of redshift-dependent magnitude changes necessary to change the inferred inverse distance ladder H_0 from 67 to the local network of distance measurements value of 73, and the impact that this would have on other cosmological inferences with SNe Ia. We find that a change of dμ/dz = 0.2 mag would be necessary to infer an H_0 in concordance with local distance measurements, and that this dμ/dz value would result in a Flat ΛCDM Ω_M = 0.23, 10σ discrepant with other cosmological probes, indicating that the precision of SNe Ia measurements preclude the necessary redshift evolution for an H_0 of 74 kms^-1 Mpc^-1. Therefore, we conclude that current SN Ia cosmology leaves little freedom for the inverse distance ladder to yield H_0 values significantly different from 67 kms^-1 Mpc^-1.
http://arxiv.org/abs/2602.14956 (pdf) (10)
The results of interferometric observations of the star-forming region S255IR in the frequency range 210–250 GHz are presented. The observations were carried out with the antenna array SMA (Hawaii, USA). Fifty-three molecules were detected, including complex organic molecules (COMs) such as CH_3CHO, CH_3CN, CH_3CH_2CN, and many others. Typical rotational temperatures in the hot core SMA1 fall in the range 100–200 K. Optical depths in the lines of methanol and some other molecules in the cores SMA1 and SMA2 were estimated. In SMA1, the optical depth of one of the strongest methanol lines, 5_-1-4_-1E, proved to be 23.8 ± 1.5. Based on this value, one can assume that the lines of other oxygen-containing COMs, such as CH_3OCHO, CH_3OCH_3, CH_3CH_2OH, which are typically much less abundant in hot cores than methanol, are optically thin in SMA1. Most of the detected molecules can be roughly divided into two groups. The molecules of the first group emit exclusively toward the hot core SMA1, while some or all lines of the molecules of the second group, in addition to SMA1, can be seen toward a ring-like structure to the west of SMA1. This structure is most likely associated with the walls of a cavity formed by high-velocity outflows driven by young stellar objects (YSOs) in molecular cores SMA1, SMA2, and possibly SMA3. The gas temperature and density in the cavity walls were estimated using methanol lines. The temperature was found to be about 50–60 K, and the density about 10^7-10^8 cm^-3. The column density of methanol near the brightness peaks in the lines of this molecule is about 5× 10^15 cm^-2. The column densities of other COMs in the ring-like structure will be determined in future studies with increased sensitivity achieved by spectral line stacking.
http://arxiv.org/abs/2602.14547 (pdf) (10)
Foreground masking and incomplete sky coverage complicate CMB polarization analyses by inducing mode coupling and imperfect E/B separation, with particularly strong impact on searches for primordial B-modes. We present SkyReconNet-P, a convolutional neural network for inpainting CMB polarization maps that extends the SkyReconNet framework to jointly reconstruct the polarization (Q,U) maps from partial-sky observations. The method combines regional processing with a hybrid design, utilizing standard convolution and dilated convolution to do a multiscale feature integration. We evaluate performance at both the map and power spectrum level using two masking scenarios: a generated random mask and the Planck 2018 common polarization inpainting mask. For both masking scenarios, SkyReconNet-P reproduces the large-scale morphology of the target maps. In power-spectrum space, we find that the reconstructed E-mode spectrum closely tracks the target at low multipoles, while small biases emerge at higher ℓ. For B-mode, the raw reconstructed spectra exhibit a larger multipole-dependent bias, which we mitigate using a simulation-based linear calibration. We show that the calibrated B-mode spectrum preserve more information by comparing it with spectrum estimation using pseudo-C_ℓ. Finally, we demonstrate cosmological parameter inference from calibrated reconstructed spectra by fitting (r, A_lens) with a Gaussian bandpower likelihood, recovering posteriors consistent with injected parameters across three test ensembles down to r ∼ 10^-3. These results support inpainting as a complementary route to cut-sky approaches when downstream pipelines can greatly benefit from statistically well-characterized, gap-filled polarization maps.
http://arxiv.org/abs/2602.14196 (pdf) (10)
We develop a self-consistent nonlinear extension of diffusive shock acceleration that incorporates cosmic ray (CR) backreaction on the shock precursor together with a physically motivated upstream-escape mechanism that produces an exponential high energy cutoff. The CR pressure gradient decelerates the upstream flow facing the shock wave, generating an extended precursor in which higher rigidity particles sample a larger cumulative velocity gradient and thereby acquire a progressively harder spectrum. Finite-size/escape effects are modeled by a momentum-dependent loss term, which naturally terminates acceleration and steepens the spectrum near the cutoff. The precursor compression ratio is not imposed as a closure condition but is determined dynamically by requiring consistency between the injection rate inferred from thermal leakage at the subshock and the injection strength demanded by the nonlinear shock modification, with CR-driven wave heating providing stabilizing negative feedback. Applying the model to young supernova-remnant-like parameters and standard one-zone Galactic diffusion, we reproduce the main features of the latest DAMPE proton spectrum: gradual hardening from hundreds of GeV to multi-TeV energies and a subsequent exponential cutoff at tens of TeV. The resulting spectral evolution follows directly from the competition between precursor-mediated nonlinear feedback and upstream escape.
http://arxiv.org/abs/2602.12948 (pdf) (10)
Aim. We present and interpret the data of the nearby hydrogen-deficient but helium-rich superluminous supernova SN 2021bnw which reached a magnitude of -20.7 at maximum luminosity in g band. Methods. We discuss the light curves and spectra of SN 2021bnw based on its spectro-photometric follow up exploiting different observational facilities. We reproduce the NIR spectrum of SN 2021bnw with TARDIS to inspect the chemical composition at late photospheric phases and identify helium features. We also use a STELLA model coupling hydrodynamics and radiation transport to constrain the physical parameters of the explosion assmunig a 56Ni+CSM scenario. Results. We suggest that SN 2021bnw was mainly powered by the interaction of the ejecta with a previously lost He-rich circumstellar material, coupled with a central power source. Conclusions. This work expands the data sample of He-rich superluminous supernovae rich (SLSNe Ib) and, assuming a single progenitor scenario, can constrain the masses and the physics of their progenitors.
http://arxiv.org/abs/2602.12560 (pdf) (10)
Polycyclic aromatic hydrocarbons (PAHs) are recognized as the primary contributors to the aromatic infrared bands (AIBs) widely observed in space. However, analyzing these AIBs remains challenging because of the immense structural diversity within the PAH family, which makes the computation of reliable reference spectra difficult. To address this, we developed an efficient graph neural network (GNN) framework that can predict PAH absorption spectra up to 10,000 times faster than traditional quantum chemical methods. We evaluated four representative GNN architectures, including graph convolutional network (GCN), graph attention network (GAT), message passing neural network (MPNN), and attentive fingerprint (AFP). The AFP model is found to deliver the best overall performance and is further trained using five different spectral distance metrics as loss functions, among which the Jensen-Shannon divergence yields the most accurate and stable results. The model performs best for PAHs containing 20-40 carbon atoms, while accuracy decreases for larger molecules, reflecting the limited availability of training data. Overall, this framework offers a fast method to generate approximate reference spectra for small- to medium-sized PAHs, supporting future AIB analysis.
http://arxiv.org/abs/2602.14993 (pdf) (9)
ISRO&aposs terrain characterization and hazard mapping from Chandrayaan-2 Orbiter High Resolution Camera (OHRC) stereo imagery were central to the safe landing of Chandrayaan-3 - the first successful landing in the lunar south polar region. However, these elevation products were generated with a proprietary pipeline and have not been publicly released. We present a 0.30 m/pixel digital elevation model (DEM) of the Chandrayaan-3 Vikram landing site using a fully open workflow based on ISIS, the Ames Stereo Pipeline, and ALE, achieving sub-meter resolution comparable to mission-reported products. The reconstruction covers 2.18 x 2.24 km with 91.2% valid pixel coverage, 8.1 cm median triangulation error, and 40-50 cm relative vertical precision. The Vikram lander and Pragyan rover are individually resolved. Geodetic alignment to an LROC NAC stereo DEM achieves approximately 30 m horizontal accuracy; pixel-wise validation at 3 m resolution confirms negligible vertical bias (median dz = +0.28 m) and robust dispersion (NMAD = 2.88 m). Stable OHRC stereo convergence requires Community Sensor Model (CSM) camera models; the legacy ISIS camera model failed across two independent sites. At 0.30 m, these DEMs complement LROC NAC DTMs (approximately 1 m), resolving sub-meter hazards below the NAC detection threshold. Applied to the extensive OHRC south polar archive, this methodology provides independent capability for hazard mapping and landing site analysis for upcoming missions including Chandrayaan-4, LUPEX, and Artemis.
http://arxiv.org/abs/2512.23040 (pdf) (9)
We generalize the short-time Fourier transform (STFT) formalism for radial velocity extraction to cases where the underlying spectral components are unknown. The method factorizes a spectroscopic time series into principal spectra and time-dependent kernels, enabling simultaneous recovery of both. In Fourier space, each inverse-wavelength slice is decomposed by singular value decomposition, and radial velocity shifts are inferred from phase differences between epochs. In the high-SNR regime, this provides a linearized and statistically tractable estimate of differential velocities. The method is validated on synthetic and SOAP simulations and applied to EXPRES observations of HD 34411 and τ Ceti, recovering coherent signals and reaching the instrumental precision limit of 30 cm/s. Apart from p-mode modulation, the residuals show no significant long-term correlations and allow the detection of signals with semi-amplitudes down to 50 cm/s with ≲10 cm/s uncertainty. The framework thus enables extreme-precision radial velocity measurements in the presence of spectral variability, representing a step toward detecting and characterizing Earth-like planets around solar-type stars.
http://arxiv.org/abs/2602.14006 (pdf) (9)
We expand on our previous study of the impact of atmospheric seeing on polarization cross-talk, and show how the formalism that was developed in that work can be applied to treat the case of spatial modulators of polarization. Beside formally demonstrating how the problem of cross-talk is fully eliminated in such devices, we also gain insight on the meaning of polarimetric noise of temporal modulation schemes in the limit of very high modulation frequency. We also describe the problem of spectrograph instabilities, and how the spectral gradients that are naturally associated with a line spectrum feed into the problem of polarimetric errors induced by mechanical vibrations, thermal drifts, and pointing jitter. Finally, we show how this formalism can be used to estimate the contribution of polarization cross-talk to the errors on the elements of the 4×4 Stokes response matrix, for the purpose of producing realistic error budgets for polarimetric instrumentation.
http://arxiv.org/abs/2602.12955 (pdf) (9)
Modern time-domain surveys produce alert streams at a scale that makes exhaustive manual inspection infeasible, requiring automated methods to identify unusual transients for follow-up. In this work, we present an unsupervised anomaly detection pipeline applied to the ZTF alert stream using the Lasair broker. We define normal objects as SN Ia, SN II, and SN Ib/c. Anomalous objects include (i) more exotic transients (AGN, TDEs, SLSNe, CVs, and nuclear transients) and (ii) supernova-labeled objects, either spectroscopically or by Lasair, with anomalous properties, such as incorrect or absent host associations, or non-supernova-like light curves. Our pipeline consists of three independently trained simple autoencoders operating on distinct alert stream data products: object features, triplet image cutouts, and light curves. Each model is trained on predominantly normal transients, and performance is assessed using the recall of exotic objects and the purity of all anomalous objects across both a spectroscopically classified held-out test set and the live alert stream. In the test set, performance is evaluated at a fixed rank corresponding to the top ten scoring candidates, while in the alert stream it is evaluated using an anomaly threshold defined from test set behavior. Across both settings, the algorithms consistently recover exotic transients and anomalous supernovae among their top-ranked candidates. Over 25 days of live alert stream application, we identify 87 unusual supernova candidates for follow-up. The overlap between anomalies flagged by different autoencoders in the test set is non-existent, and in the alert stream is small, with maximum overlap between any two algorithms being 11 objects. The framework is data-efficient, requiring only a few thousand training examples, making it well suited for early and ongoing application to the Rubin Observatory alert stream.
http://arxiv.org/abs/2602.15021 (pdf) (8)
Cross-survey generalization is a critical challenge in stellar spectral analysis, particularly in cases such as transferring from low- to moderate-resolution surveys. We investigate this problem using pre-trained models, focusing on simple neural networks such as multilayer perceptrons (MLPs), with a case study transferring from LAMOST low-resolution spectra (LRS) to DESI medium-resolution spectra (MRS). Specifically, we pre-train MLPs on either LRS or their embeddings and fine-tune them for application to DESI stellar spectra. We compare MLPs trained directly on spectra with those trained on embeddings derived from transformer-based models (self-supervised foundation models pre-trained for multiple downstream tasks). We also evaluate different fine-tuning strategies, including residual-head adapters, LoRA, and full fine-tuning. We find that MLPs pre-trained on LAMOST LRS achieve strong performance, even without fine-tuning, and that modest fine-tuning with DESI spectra further improves the results. For iron abundance, embeddings from a transformer-based model yield advantages in the metal-rich ([Fe/H] > -1.0) regime, but underperform in the metal-poor regime compared to MLPs trained directly on LRS. We also show that the optimal fine-tuning strategy depends on the specific stellar parameter under consideration. These results highlight that simple pre-trained MLPs can provide competitive cross-survey generalization, while the role of spectral foundation models for cross-survey stellar parameter estimation requires further exploration.
http://arxiv.org/abs/2602.13875 (pdf) (8)
Massive stars dominate the radiative and mechanical feedback of young stellar populations, yet their intense ultraviolet fields and strong winds are typically presumed to preclude Earth-like habitability. We quantify this expectation by mapping time dependent habitable zones (HZs) for solar-metallicity stars with initial masses of 0.8-120 M_⊙. From rotating and non-rotating GENEC tracks we derive bolometric “climate” HZ boundaries and enforce XUV energy-limited escape and wind ram-pressure retention constraints for a dipole-magnetized Earth analogue. The operational inner edge is set by the most restrictive limit, and we measure the annulus lifetime, the longest continuous residence at fixed orbit, and the maximum number of dynamically packed terrestrial planets it can host. We find a sharp main-sequence ceiling: while a 9 M_⊙ star sustains an operational HZ for ∼ 30 Myr at ∼ 70-130 AU, the main-sequence annulus becomes brief and extremely narrow by 12 M_⊙ and disappears by 15 M_⊙. Post main-sequence evolution can reopen HZs up to ∼ 25-30 M_⊙, but only for ∼ 0.03-1.5 Myr at hundreds to ∼ 10^3 AU, disappearing by ∼ 40 M_⊙. Rotation modestly increases habitable lifetimes near the upper main sequence without altering the high mass ceiling. Initial Mass Function (IMF) weighting shows that massive stars contribute only ∼ 10^-4 of the habitable planet-time budget. Even so, they still add of order a few 10^5 operationally habitable Earth analogues to the Milky Way at any instant. This implies that massive star systems are unlikely to dominate the Galaxy wide habitability budget, but they may still provide a set of short-lived, observationally distinct targets for biosignature searches.
http://arxiv.org/abs/2602.14804 (pdf) (8)
Astrochemical observations have revealed a surprisingly high level of chemical complexity, including long carbon chains, in the earliest stages of Sun-like star formation. The origin of these species and whether they undergo further growth, possibly contributing to the molecular complexity of planetary systems, remain open questions. We present recent observations performed using the 100-m Green Bank Telescope of the prestellar core L1544 and the protostellar system IRAS 16293-2422. In L1544, we detected several complex carbon-bearing species, including C_2S, C_3S, C_3N, c-C_3H, C_4H, and C_6H, complementing previously reported emission of cyanopolyynes. In IRAS 16293-2422, we detected c-C_3H and, for the first time, HC_7N. Thanks to the high spectral resolution, we refine the rest frequencies of several c-C_3H and C_6H transitions. We perform radiative transfer analysis, highlighting a chemical difference between the two sources: IRAS 16293-2422 shows column densities 10-100 times lower than L1544. We perform astrochemical modeling, employing an up-to-date chemical network with revised reaction rates. The models reproduce the general trends, with cyanopolyyne and polyynyl radical abundances decreasing as molecular size increases, but they underestimate the abundances of cyanopolyynes longer than HC_5N by up to two orders of magnitude. Current models, which include the dominant neutral-neutral formation routes, cannot account for this discrepancy, suggesting that the chemical network is incomplete. We propose that additional ion-molecule reactions are crucial for the formation of these species. Developing a more comprehensive chemical network for long carbon chains is essential for accurately interpreting present and future observations.
http://arxiv.org/abs/2602.14533 (pdf) (8)
A model-agnostic description of the baryon acoustic oscillation (BAO) feature in redshift space requires a number of ingredients. Physically, one must describe the impact of cosmological bulk flows which progressively and anisotropically smear out the feature over time. One must also model the effects of the scale dependence of tracer bias and the mode coupling between short and long scales. All of these can be incorporated using the Zel&aposdovich approximation alone, without reference to any particular cosmological model. On the technical front, one needs a robust, complete and cosmology-independent basis to describe the shape of the real space BAO feature in linear theory, which can then be propagated to the nonlinearly evolved, measured feature in redshift space. In this work, we describe how these ingredients – which we have systematically constructed in recent work – come together in an accurate framework capable of describing the BAO-scale pairwise measurements of state-of-the-art galaxy surveys. Using mock observations and N-body simulations, we show that our template-free framework can produce unbiased and precise cosmological constraints for samples with realistic levels of nonlinearity. This work represents one of the final steps in constructing a data-ready analysis framework for model-agnostic cosmological inference from the BAO feature.
http://arxiv.org/abs/2602.14458 (pdf) (8)
We present a systematic investigation of the radiative properties and optical appearance of rotating SV black holes surrounded by a thin accretion disks, and mainly analyze the influence of the regularization parameter g on related observables. The results show that although the kinetic quantities and the location of the innermost stable circular orbit (ISCO) depend on the regularization parameter g, the radiative efficiency of the rotating SV black hole is the same as its Kerr counterpart. Within the Novikov-Thorne thin-disk model, the radiative flux, effective temperature, and spectral luminosity are studied, and by adopting observational parameters relevant to SgrA* and M87*, concrete examples of the rotating SV black holes are calculated and compared with that of the Kerr black holes. The results show that the parameter g suppresses the maximum values of these quantities. In addition, using a backward ray-tracing technique, we numerically simulate the optical appearance of rotating SV black holes and analyze the corresponding intensity images, redshift and observed flux distributions. Our results show that these quantities are affected by g. In particular, as g increases, the observed intensity is significantly suppressed and the photon ring region has remarkable increase in its width. Our findings suggest that accretion-disk-related observables may provide important avenues to distinguish rotating SV black holes and Kerr black holes, and offer theoretical guidance for future high-resolution observations.
http://arxiv.org/abs/2602.14000 (pdf) (8)
Spectral and spatial fringes in polarized light are produced by the interference of transmitted and reflected waves at the interface between materials with different indexes of refraction. These instrumental artifacts can affect the accuracy of optical designs conceived for high-sensitivity spectroscopy and polarimetry. We review the fundamental mechanisms that are responsible for these artifacts and the possible design pathways that allow us to mitigate them. In order to do so, we also present an approximate treatment of the problem of the transmission and reflection of light through (possibly absorptive) birefringent layers, relying on a few fundamental results that can be found in the already extensive literature on the subject. Unfortunately, many of these results remain the domain of a niche of investigators working in the field of thin films and optical coatings, and are often overlooked even by experienced designers of spectro-polarimetric instrumentation. The treatment presented in this work is limited to isotropic materials and uniaxial crystals, which are the most common types of optics employed in polarimetric instrumentation, and it fundamentally relies on the approximation of small birefringence for its implementation. An extensive set of modeling examples is provided to highlight the salient characteristics of polarization fringes, as well as to assess how approximate treatments such as this compare to exact but more computational expensive formulations of the problem such as Berreman&aposs calculus.
http://arxiv.org/abs/2602.13777 (pdf) (8)
Water ice has been detected in several supernova remnants despite the strong heating and radiation in these environments. This challenges standard expectations for dust survival. Using two dimensional magnetohydrodynamic simulations, we study how a supernova shock interacts with a dense interstellar cloud. The simulations show that the shock naturally compresses the cloud into dense structures similar to those inferred in well known remnants. Although temperatures remain high in the adiabatic phase, simple considerations indicate that cooling would act quickly once included. These results suggest that shock cloud interactions create the physical conditions needed for water ice to form. Future work including radiative cooling and grain surface chemistry will allow direct modelling of ice growth in these compressed regions.
http://arxiv.org/abs/2602.13644 (pdf) (8)
Understanding turbulence in interplanetary coronal mass ejections is fundamental to space plasma research and critical for assessing the impact of space weather on geospace. Turbulence governs energy cascade, plasma heating, magnetic reconnection, and solar wind magnetosphere coupling, thereby influencing both ICME evolution and geoeffectiveness. While previous event-based and statistical studies have examined ICME turbulence and its radial evolution in great detail, no significant measurements of ICME magnetic turbulence at a specific vantage point have been made using multiple observatories separated azimuthally. Here, we present the first multipoint analysis of MHD turbulence across ICME plasma regions, using four spacecraft at the Sun-Earth L1 point, separated by 80 RE along the dawn-dusk direction. Previous studies reveal that ICME shocks, sheaths, and magnetic clouds are highly non-uniform, with strong azimuthal variability. Using high-resolution magnetic field observations from ISRO&aposs Aditya-L1, NASA&aposs Wind and ACE, and NOAA&aposs DSCOVR, we analyze turbulence associated with the 10th October 2024 solar storm, which triggered the second-strongest geomagnetic storm of solar cycle 25. Our results reveal significant variability and differing turbulence maturity across small separations, supported by analysis of field-aligned and perpendicular magnetic field cascades, indicating strong anisotropies. Sheath turbulence is substantially modified by shock induced energy injection. Evidence of compressible turbulence and plasma energization at the flux rope interaction region indicates that internal processes, such as magnetic reconnection, strongly influence ICME plasma evolution, highlighting pronounced spatial variability in turbulence and plasma states observed by multiple L1 monitors near Earth and underscoring their potential role in space weather impacts.
http://arxiv.org/abs/2602.14781 (pdf) (7)
Following the detections of the first extragalactic “Luminous Red Nova” (LRN) M31 RV in 1989, and its first Galactic counterpart V4332 Sgr in 1994, there have been many discoveries of similar, or closely related, objects. They are important because they bridge the luminosity gap between the brightest novae and supernovae, a largely unexplored parameter space. The cause of eruptions in LRNe is still unclear, a stellar merger being the most favored mechanism. However, barring V1309 Sco, there has been no direct evidence for a merger in the other objects. Here we present strong evidence that V4332 Sgr was a merger event. High resolution infrared observations of the CO fundamental band show an unusually small ^12C/^13C ratio of 3.5±1. This indicates that a violent event had occurred, whose effects penetrated deep enough to allow CNO cycle processed ^13C in the inner H burning shell to be brought to the surface. We rule out planetary ingestion, and propose that the eruption was due to a merger between V4332 Sgr and a companion star. It is shown that V4332 Sgr was likely surrounded by an edge-on disk before its eruption. If this disk was a flattened common envelope containing V4332 Sgr and a companion star, then a merger scenario would not be inconsistent. Furthermore, V4332 Sgr had multiple outbursts, previously unreported but an important piece of information, since multiple outbursts are a trait shared by many LRNe.
http://arxiv.org/abs/2602.13291 (pdf) (7)
Artificial Intelligence (AI) has transformed robotics, healthcare, industry, and scientific discovery, yet a major frontier may lie beyond Earth. Space exploration and settlement offer vast environments and resources, but impose constraints unmatched on Earth: delayed/intermittent communications, extreme resource scarcity, heterogeneous expertise, and strict safety, accountability, and command authority. The key challenge is auditable coordination among specialised humans, robots, and digital services in a safety-critical system-of-systems. We introduce Agent Mars, an open, end-to-end multi-agent simulation framework for Mars base operations. Agent Mars formalises a realistic organisation with a 93-agent roster across seven layers of command and execution (human roles and physical assets), enabling base-scale studies beyond toy settings. It implements hierarchical and cross-layer coordination that preserves chain-of-command while allowing vetted cross-layer exchanges with audit trails; supports dynamic role handover with automatic failover under outages; and enables phase-dependent leadership for routine operations, emergencies, and science campaigns. Agent Mars further models mission-critical mechanisms-scenario-aware short/long-horizon memory, configurable propose-vote consensus, and translator-mediated heterogeneous protocols-to capture how teams align under stress. To quantify behaviour, we propose the Agent Mars Performance Index (AMPI), an interpretable composite score with diagnostic sub-metrics. Across 13 reproducible Mars-relevant operational scripts, Agent Mars reveals coordination trade-offs and identifies regimes where curated cross-layer collaboration and functional leadership reduce overhead without sacrificing reliability. Agent Mars provides a benchmarkable, auditable foundation for Space AI.
http://arxiv.org/abs/2602.13074 (pdf) (7)
In this paper, we revisit the linear perturbation (including the comoving curvature perturbation and field perturbation) dynamics in the SR-USR-SR inflation with instantaneous transitions. Using the junction method and asymptotic expansions of Hankel functions, we derive accurate asymptotic expressions for the time evolution of mode functions and the resulting power spectrum, based on three systematic rules for identifying the dominant terms across transitions. Our results reveal that a finite dip of the final power spectrum arises from the cancellation between two growing modes within the linear perturbation theory, rather than between constant and growing terms as previously suggested. We also provide analytical descriptions of the amplitude enhancement and oscillatory features in the linear power spectrum, in agreement with numerical computations. These simple, tractable formulas not only facilitate theoretical calculations but also yield testable predictions for future CMB observations.
http://arxiv.org/abs/2602.13051 (pdf) (7)
Observational appearance of debris disks is largely controlled by collisional grinding of their dust grains. However, the mechanical strength of dust at sizes in the micrometer to millimeter range is poorly known. Recent studies suggested that dust particles in the Solar system might have a higher critical fragmentation energy Q_D^* value than previously anticipated. Another recent study considered the Fomalhaut debris disk and found lower Q_D^* values to provide better fits to the data. In order to constrain the mechanical strength of dust, we investigate collisional evolution of debris disks with Q_D^* prescriptions differing by ∼ 3 orders of magnitude. We find that, above a certain threshold Q_D^* value, the disk&aposs collisional evolution is dominated by rebounding – rather than disruptive or cratering – collisions. Rebounding (a.k.a. bouncing) collisions are those in which both impactors survive, being only slightly eroded and producing fragments that only carry a minor fraction of their mass. We show that disks dominated by rebounding collisions would have brightness profiles increasing outward outside the parent belt. Since such profiles are not observed, this places an upper limit on how hard the debris dust is allowed to be in order not to violate the observations. We derive an approximate analytic expression for this limit: Q_D^* ≈ (1/8) v_K^2(r) for grains close to the radiation pressure blowout size, where v_K in the Keplerian circular speed at a distance r from the star. This implies Q_D^* ≲ 10^9...10 erg g^-1 for micrometer-sized grains in typical debris disks. Even though rebounding collisions are not expected to affect debris disk evolution significantly, we emphasize that these collisions are actually much more frequent than disruptive and cratering ones in all debris disks.
http://arxiv.org/abs/2602.12930 (pdf) (7)
We present time-series medium-band (R 20-40) observations of the third interstellar object 3I/ATLAS (C/2025 N1) obtained with the 7-Dimensional Telescope (7DT), enabling spatially resolved monitoring of its gas and dust activity from 2025 July to September. The m400-band image (lambda_c = 400 nm, Delta lambda approx 25 nm) reveals the emergence of pronounced and spatially extended CN emission at heliocentric distances r_h < 3 au. This onset is consistently identified across multiple diagnostics, including a break in the light-curve evolution, excess reflectance, inward expansion of annular excess beyond 10,000-20,000 km, growth of the coma half-light radius from 11,000 to 19,000 km, and a rapid rise in the CN production rate Q_CN relative to Af rho. We further separate the CN-emitting and dust-scattered components through two-dimensional surface-brightness fitting into inner (dust) and outer (gas) components. The outer component preserves a nearly constant profile shape, varying only in normalization, implying relatively fast expansion of CN-bearing molecules. Together, these results reveal a transition in the optical from dust-dominated scattering at large heliocentric distances to volatile-driven, gas-dominated activity as 3I/ATLAS enters the inner Solar System. The timing and characteristics of the CN activation resemble the volatile enhancement observed in 2I/Borisov, suggesting that both known active interstellar objects exhibit comparable activation behavior at heliocentric distances of 2-3 au.
http://arxiv.org/abs/2602.13441 (pdf) (7)
Re-derivation of the sensitivities of all of the Space Telescope Imaging Spectrograph (STIS) observing modes were required after major updates were introduced to the model atmospheres of the three primary standard stars. The new predicted continuum fluxes were up to 2-3% different from the models used to originally calibrate STIS. This work focuses on the re-derivation of spectral sensitivities for the five STIS low-resolution (L-mode) gratings: G140L, G230L, G230LB, G430L, and G750L, which span wavelengths from the far-ultraviolet through the near infrared. Updated photometric throughput tables were delivered to the Calibration Reference Data System (CRDS) on April 7, 2022 and April 14, 2023, which triggered a recalibration of all historical STIS datasets taken with these modes. The sensitivities derived from each of the standard stars typically agree with one another to better than 1%, though discrepancies as large as 1.5% are found in spectral regions most impacted by hydrogen absorption.
http://arxiv.org/abs/2512.17885 (pdf) (7)
Magnetic fields are critical at many scales to galactic dynamics and structure, including multiphase pressure balance, dust processing, and star formation. Dynamo action determines their dynamical structure and strength. Simulations of combined large- and small-scale dynamos have successfully developed mean fields with strength and topology consistent with observations but with turbulent fields much weaker than observed, while simulations of small-scale dynamos with parameters relevant to the interstellar medium yield turbulent fields an order of magnitude below the values observed or expected theoretically. We use the Pencil Code accelerated on GPUs with Astaroth to perform high-resolution simulations of a supernova-driven galactic dynamo including heating and cooling in a periodic domain. Our models show that the strength of the turbulent field produced by the small-scale dynamo approaches an asymptote at only modest magnetic Prandtl numbers. This allows us to use these models to suggest the essential characteristics of this constituent of the magnetic field for inclusion in global galactic models. The asymptotic limit occurs already at magnetic Prandtl number of only a few hundred, many orders of magnitude below physical values in the the interstellar medium and consistent with previous findings for isothermal compressible flows.
http://arxiv.org/abs/2602.13198 (pdf) (7)
We introduce , a CUDA/C-based Monte Carlo radiative transfer code accelerated using graphics processing units (GPUs). derives from the CPU-based code and offloads the most computationally expensive stages of the calculation – superphoton generation, sampling, tracking, and scattering – to the GPU. Whereas handles photons sequentially, processes large numbers of superphotons concurrently, leveraging the single-instruction, multiple-thread (SIMT) execution model of modern GPUs. Benchmarks demonstrate a speedup of about 12× relative to the original CPU implementation on a single GPU, with runtime limited primarily by register pressure rather than compute or memory bandwidth saturation. We validate the implementation through analytic tests for a optically thin synchrotron sphere, as well as comparisons with for scattering synchrotron sphere and GRMHD simulation data. Relative errors remain below a percent level and convergence is consistent with the expected N_s^-1/2 Monte Carlo scaling. By significantly reducing computational costs, GPUmonty enables the extensive parameter space surveys and faster spectra modeling required to interpret horizon-scale observations of supermassive black holes. is publicly available under the GNU General Public License.
http://arxiv.org/abs/2602.13137 (pdf) (7)
Double white dwarfs (DWDs) are by far the most common compact binaries in the Milky Way, are important low-frequency gravitational-wave sources, and in some cases merge to become Type Ia supernovae. So far, no DWD has been identified solely through relativistic Doppler beaming, even though the beaming amplitude directly relates to the radial velocity semi-amplitude. In this work, we initiate a comprehensive binary population synthesis using SeBa and incorporate the resulting binaries into a tripartite Galaxy model. Our proof-of-concept simulations demonstrate that the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) can reliably recover relatively bright (r ≲20mag) unequal-mass binaries in compact orbits with P ≈ 10-600 minutes with moderate to high inclinations. We find that LSST can detect at least 287 short-period DWDs, of which 47 are LISA-detectable gravitational wave sources. LSST lightcurves allow us to readily determine the period and fully characterize the orbit, in contrast with the challenges of orbit determination for DWDs in spectroscopic searches. The formation of unequal mass, short-period DWDs strongly depends on the assumptions regarding the mass-transfer phases during binary population synthesis, and the total number and characteristics of Doppler-beamed DWD systems observed in LSST will provide new tests of models of stellar binary evolution. Here, we lay the foundation for the comprehensive integration of synthetic Galactic binary population into realistic LSST survey simulations, thereby enabling quantitative forecasts of the number and characteristics of any binary sub-population during the LSST era.
http://arxiv.org/abs/2602.12307 (pdf) (7)
We present GAMERA-OP (Orthogonal-Plus), a three-dimensional finite-volume magnetohydrodynamics (MHD) solver for orthogonal curvilinear geometries. The solver advances magnetic fields using constrained transport to preserve ∇·𝐁=0 to machine precision and employs geometry-consistent high-order reconstruction with an enhanced Partial Donor Cell method (e-PDM) that accounts for geometry curvature. Flexible numerics include various numerical fluxes and time integrators. In axial symmetric coordinates, angular momentum are preserved to round-off, and a ring-averaging treatment near the axis relaxes CFL constraints while maintaining divergence-free magnetic fields. Optional capabilities include the semi-relativistic (Boris) correction, background-field splitting, and an anisotropic MHD formulation. Rewritten in C, the code adopts a modular design that simplifies case setup and facilitates the addition of physics modules and coupling to other first-principles codes. Standard benchmarks across multiple geometries verify the code&aposs high accuracy, low numerical diffusion, and robust handling of coordinate singularities and rotating flows. GAMERA-OP provides a practical, high-order framework for space and astrophysical plasma applications where orthogonal curvilinear coordinates and exact angular-momentum conservation are advantageous.
http://arxiv.org/abs/2602.14515 (pdf) (6)
This paper proposes a novel quantum effect wherein particle spectra show extra broadening due to local short-time acceleration (as well as in a deSitter spacetime background). This effect provides a simple interpretation for the acceleration interpolation relation required to modify the kinematics of a test particle in the Modified Inertial interpretation of Modified Newtonian Dynamics (MOND). This effect can be regarded as a generalization of the thermal blackbody spectrum generated by the Unruh effect (which arises from long-time uniform acceleration in a flat background) to the scenario of local short-time non-uniform acceleration. This effect offers a unified framework for understanding the accelerated expansion of the universe and the anomalies in galactic rotation curves or radial acceleration. It is worth emphasizing that this modified kinematic interpretation of MOND necessitates a quantum equivalence principle as its physical foundation, that is, extending the equivalence at the level of expectation values (first-order moments) in the classical equivalence principle to the equivalence at the level of second-order moment quantum fluctuations. Therefore, the role of this effect within the frameworks of quantum reference frames and quantum gravity is also discussed to some extent.
http://arxiv.org/abs/2602.13463 (pdf) (6)
Protoplanetary disc mass is one of the most fundamental properties of a planet-forming system, as it sets the total mass budget available for planet formation. However, obtaining disc mass measurements remain challenging, since it is not possible to directly detect H_2, and CO abundance ratios are poorly constrained. Dynamical measurements of the disc mass are now possible, but they are not suited to all discs since the measurements typically require well-behaved emission surfaces. A long-standing method is to obtain continuum flux measurements from the dust emission, and convert to a total disc mass by assumption of the dust-to-gas mass ratio, ε. This quantity is poorly constrained in protoplanetary discs. % We investigate the impact of ε on the morphology of planet-containing hydrodynamical simulations of dusty protoplanetary accretion discs, and suggest that if a planet mass estimate can be obtained, then disc morphology could be used to constrain ε in observed systems relative to each other, improving the total disc mass estimates of protoplanetary discs.
http://arxiv.org/abs/2602.12376 (pdf) (6)
X-ray pulse-profile modeling of millisecond pulsars offers a direct route to measuring neutron star masses and radii, thereby constraining the dense-matter equation of state. However, standard analyses typically rely on ad hoc hotspot parameterizations rather than self-consistent physical models. While connecting surface heating directly to the magnetospheric geometry provides a more natural physical pathway, computing global magnetospheric solutions is too computationally expensive to perform on-the-fly during parameter inference. In this work, we bridge this gap by deriving fully analytic, first-principles expressions for surface return currents in mixed dipole–quadrupole magnetospheres. Working within force-free electrodynamics, we generalize the field-aligned current invariant Λ, the crucial scalar that maps the far-zone magnetic structure to the near-zone heating rate, from the standard dipole approximation to arbitrary quadrupolar configurations. We demonstrate that even when the quadrupole component is sub-dominant in the far zone (the mixing regime), using a dipole-based heating prescription fails to capture the significant enhancement or suppression of the return-current density on the polar cap. Our consistent quadrupole-aware framework reveals that these multipolar currents redistribute the surface heating, leading to systematic discrepancies in predicted pulse profiles that are amplified by atmosphere beaming and can reach ∼ 30
http://arxiv.org/abs/2602.14838 (pdf) (6)
It has long been suspected that black hole accretion-outflow coupling is invariant from the stellar to supermassive scales. Stellar mass black hole accretion flows are known to launch jets and outflows as they transition through critical accretion rate thresholds, with values well constrained observationally owing to their short evolutionary timescales. In contrast, accretion flows in typical supermassive black hole (SMBH) systems (those in active galactic nuclei) evolve over thousands of years, making the critical transitions at which jets are launched impossible to constrain in individual systems. Tidal disruption events (TDEs) provide the unique opportunity to witness the birth and evolution of an accretion flow onto a SMBH which evolves on timescales of years. Here we show that TDEs launch outflows during a super-Eddington accretion phase and a second, physically distinct outflow, at a critical accretion rate of L_crit≈0.02 L_Edd, the same as the critical accretion rate for state transitions observed in accreting stellar mass black holes. This work naturally explains the mechanism, observed properties, and detection rate for prompt and delayed outflows observed in TDEs, which until now have been open problems. More broadly, we demonstrate that SMBHs exhibit the same accretion-outflow coupling as stellar mass black holes and that the critical low accretion rate threshold for jet formation in black holes is scale invariant.
http://arxiv.org/abs/2602.14343 (pdf) (6)
In this study, we aim to constrain the atmospheric composition of GJ 1214 b using all available transits observed with the upgraded CRIRES+ spectrograph at the VLT by searching for the signatures of water vapour, methane, and carbon dioxide. We analysed eight CRIRES+ transit datasets covering the K band (1.90-2.45 microns) at a resolving power of R 100,000. We used the SysRem algorithm to correct for telluric and stellar contributions and employed the cross-correlation technique with templates from petitRADTRANS to search for H2O, CH4, and CO2. Injection-recovery tests across a grid of metallicities (Z) and cloud-deck pressures (pc) were performed to quantify detection limits. We also generated predictions for ANDES observations using end-to-end simulated datasets with EXoPLORE. We detect no significant H2O, CH4, or CO2 signatures. Injection-recovery tests show that such non-detections exclude atmospheres with low-altitude clouds and moderate or low metallicities. CH4 yields the tightest empirical limits, with CO2 unexpectedly ruling out intermediate metallicities ( 100xsolar) with clouds deeper due to its rapidly rising opacity in compressed, high-Z atmospheres. Our constraints are in line with either a high-Z or a high-altitude aerosol layer, in agreement with recent JWST inferences. The combined analysis of eight CRIRES+ datasets provides the most stringent high-resolution constraints on the atmospheric properties of GJ 1214 b to date. Simulations of a single transit observed with ANDES on the ELT predict modest improvements for H2O, a substantially expanded detectable region for CH4, and the strongest gains for CO2, making the latter a particularly effective tracer for characterising high-metallicity atmospheres in sub-Neptunes.
http://arxiv.org/abs/2602.14294 (pdf) (6)
The advent of new and more powerful infrared spectrographs has significantly motivated the advancement of the study of atomic and molecular line lists and stellar atmosphere models. While optical abundance determinations rely on extensively validated line lists and modeling frameworks, infrared measurements still face larger uncertainties, largely driven by the choice of atmospheric models and the quality of the available atomic data. In this work, we aim to deliver a homogeneous and reproducible set of atomic absorption lines in the Y, J, and H bands (9800 - 18000 (Angstrom)), based exclusively on laboratory atomic data. We analyse CRIRES spectra of six Gaia FGK Benchmark Stars spanning a wide range in effective temperature, surface gravity, and chemical composition. Synthetic spectra are computed using the benchmark stellar parameters, and each transition is evaluated independently in every star through a quantitative sequence that examines line depth, saturation, blending (purity), and the agreement between observed and synthetic line profiles. We identify a set of robust atomic transitions in these bands that remain consistent across the full range of stellar parameters represented in our sample. Lines of alpha-elements such as Mg I, Si I, and Ca I, together with several Fe I transitions, satisfy all robustness criteria. Among the neutron-capture species explored, only Sr II provides lines that consistently meet our requirements. Beyond the specific list of accepted transitions, this study demonstrates that a fully quantitative, multi-criteria framework provides a transparent and reproducible foundation for near-infrared line validation as laboratory data, stellar atmosphere models, and instrumentation continue to improve.
http://arxiv.org/abs/2602.13992 (pdf) (6)
An analysis of the Millimetre Astronomy Legacy Team 90 GHz (MALT90) survey has produced a sample of 27 candidate dense molecular clumps with large collapse motions, as revealed by large “blue” asymmetrical line profiles of the optically thick line. %with respect to the optically thin line. New, more sensitive molecular line observations of this sample, conducted with the Mopra 22-m telescope, confirm the blue asymmetries in the line profiles, with large, positive values of the asymmetry parameter A (A̅_HCO^+ = 0.69±0.01), and positive, but smaller asymmetries in the and lines: (A̅_HCN = 0.35±0.01 and A̅_HNC = 0.28±0.01), as expected for a less optically thick tracer in collapsing clumps. The small, positive mean asymmetry parameters for and , A̅_C_2H = 0.15±0.02 and A̅_H^13CO^+ = 0.18±0.03, likely indicate slightly optically thick emission for at least some clumps. The hyperfine ratios for are in their optically thin, LTE, values, but for they are not; the F=1 → 1 hyperfine line shows abnormally weak intensities. A simple two-component model shows that self-absorption of the background F = 1 → 1 hyperfine line by the main F = 2 → 1 hyperfine line of a cold, foreground, redshifted cloud can reproduce the observed hyperfine intensities and match the and line profiles. All of these results are consistent with self-absorption of the optically thick lines on the red side of the profile, as expected for collapsing clumps. A simple two-cloud model suggests that this sample represents dense clumps with extreme collapse velocities, V_inf∼ 2.4 .
http://arxiv.org/abs/2602.13677 (pdf) (6)
The systematic variation of solar active region (AR) properties with their magnetic flux has been the subject of numerous studies but the proposed scaling laws still vary rather widely. A correct representation of these laws and the deviations from them is important for modelling the source term in surface flux transport and dynamo models of space climate variation, and it may also help constrain the subsurface origin of active regions. Here we determine active region scaling laws based on the recently constructed ARISE active region data base listing bipolar ARs for cycle 23, 24 and 25. For the area A, pole separation d and tilt angle γ we find scalings against magnetic flux Φ and heliographic latitude λ. Residuals from these relations are also modelled. These scaling relations are recommended for use in space climate research for the modelling of future data or missing past data, as well as for the identification of candidate rogue ARs. We confirm that the tilt angle distribution of non-Hale ARs shows a significant excess at low tilts (anti-Hale ARs). In contrast to earlier studies we show that neither the anti-Hale ARs nor non-Hales in general follow Joy&aposs law: instead, their tilt angle distribution is best represented by vanishing mean tilt. These results are most easily reconciled with a scenario where the AR flux loops originate in the deep convective zone or below, gaining tilt during their rise under the action of the Coriolis force. A small fraction of the loops is subjected to extreme, intermittent torques resulting in either very large tilts or anti-Hale orientation. Anti-Hale ARs are suggested to be fully curled `XO-loops&apos, and their excess is caused by a simple mechanical effect, as the contact of their legs increases resistance against further deformation by the torque.
http://arxiv.org/abs/2602.12890 (pdf) (6)
Kink oscillation frequency is a key parameter for coronal seismology. It is still unclear how gravitational stratification affects the kink frequency in curved coronal loops. This work aims to investigate the effect of gravitational stratification on the frequency of kink oscillations in curved coronal loops and discuss their seismological potential. We conduct numerical computations within the ideal MHD framework to study different kink polarizations and harmonics in a curved, gravitationally stratified coronal loop. The oscillation frequencies derived from the Lagrangian displacement are compared with the WKB approximation. For the vertically polarized fundamental mode, the oscillation frequency deviates from the WKB approximation by about 18
http://arxiv.org/abs/2602.12777 (pdf) (6)
Asymptotic Giant Branch (AGB) stars significantly contribute to the chemical composition of the universe. In their outflows, complex chemistry takes place, which critically depends on the local temperature. Therefore, if we want to accurately model the AGB environment, we need accurate cooling rates. The CO molecule is abundant in AGB outflows, and has a dipole moment, which enables it to cool through emission from its rotational transitions. We therefore expect it to significantly contribute to cooling in this environment, even at low temperatures (10 K ⩽ T⩽ 3000 K). Currently, CO cooling rates are available for ISM-like conditions, which encompasses a different parameter regime, with generally lower densities and velocity gradients, compared to AGB winds. Therefore, these ISM cooling rates might not be applicable to the AGB regime. In this paper, we compute CO cooling rates for hydrodynamics simulations of AGB outflows. To evaluate the net cooling rate, we calculate the energy level distribution of CO self-consistently, using the non Local Thermodynamical Equilibrium (NLTE) line radiative transfer code Magritte. We verify whether already existing CO cooling rate prescriptions for the interstellar medium (ISM) are applicable for this regime. We noticed minor differences between these prescriptions and our calculated cooling rates in general. However, when used far outside their originally intended parameter regimes, significant differences occur. Therefore, we propose a new CO cooling rate prescription for the AGB environment and we study how the computed cooling rate varies depending on input parameters.