Most interesting arXiv papers
Friday 16-01-2026
http://arxiv.org/abs/2601.09948 (pdf) (66)
We present new measurements of galaxy luminosity functions (LFs) from JWST/NIRCam imaging over the redshift range z=4.5-9.7, using photometric catalogs from JADES and public extragalactic fields. Our analysis includes rest-frame UV and B-band LFs, as well as apparent LFs in F090W, F115W, F200W, F356W, and F444W. We present the first constraints on the rest-frame B-band LF at z 7-8 and extend existing measurements at z 5 to M(B) = -18 mag. The B-band LFs evolve more strongly with redshift than UV LFs, though both decline more gradually than predicted by simulations at z>5. No single existing simulation reproduces all observed trends, with discrepancies likely driven by assumptions about binary evolution and stellar population synthesis models. The apparent LFs in F356W and F444W show hints of a bright-end excess at all redshifts, extending to fainter magnitudes at higher redshift. While extreme emission line galaxies may partially account for it, the excess may also indicate a population of moderately red, optically bright sources - potentially dusty star-forming galaxies or obscured AGNs. Finally, we find that rest-frame B-band luminosity correlates more tightly with stellar mass than UV, making it a powerful tracer of mass assembly and reinforcing the diagnostic value of rest-frame optical LFs in uncovering the physical processes that drive early galaxy formation.
http://arxiv.org/abs/2601.10372 (pdf) (54)
We investigate the correlation between ionized gas outflows, jets, and star formation in a sample of 42 local type-I active galactic nuclei (AGNs) exhibiting significant [O III] outflows. This study uses both new submillimeter (sub-mm) observations and archival data from the James Clerk Maxwell Telescope. Our analysis, which includes a correction for jet emission in the sub-mm bands, fitting spectral energy distribution, and analyzing spectra, enables us to derive star-formation rates (SFRs) through various methods. By comparing radio power and SFRs, we select a sub-sample of jetted AGNs of which radio emission is mostly from the jets. We find that jetted AGNs predominantly lie above the main sequence of star-forming galaxies, suggesting a correlation between jet activity and star formation. By comparing dust extinction, we demonstrate that jetted AGNs do not have more dust which is the fuel of both star formation and AGN activity. Therefore, this correlation is more likely to arise from AGN feedback. We also find that the Eddington ratio does not impact the specific SFRs (sSFRs) of our sample. Additionally, for jetted AGNs, stronger radio emission corresponds to higher sSFRs, suggesting that jet emission may promote star formation, i.e., positive feedback. Our results not only shed light on the feedback mechanisms of AGNs but also underscore the complex interplay between black hole activity and star formation in galaxy evolution.
http://arxiv.org/abs/2601.09780 (pdf) (54)
We present [O_III] 88μm observations towards four Herschel-selected dusty star-forming galaxies (DSFGs; log_10 μL_IR/L_⊙ = 13.5 - 14 at z = 2.9 - 4) using the Atacama Compact Array (ACA) in Bands 9 and 10. We detect [O_III] emission in all four targets at >3σ, finding line luminosity ratios (L_[O_III] / L_IR = 10^-4.2 to 10^-3) similar to local spiral galaxies, and an order of magnitude lower when compared with local dwarf galaxies as well as high-redshift Lyman-break galaxies. Using the short-wavelength capabilities of the ACA, these observations bridge the populations of galaxies with [O_III] emission at low redshift from space missions and at high redshift from ground-based studies. The difference in [O_III] emission between these DSFGs and other high-redshift galaxies reflects their more evolved stellar populations (> 10 Myr), larger dust reservoirs (M_dust ∼ 10^9 - 11 M_⊙), metal-rich interstellar medium (Z ∼ 0.5 - 2 Z_⊙), and likely weaker ionization radiation fields. Ancillary [C_II] emission on two targets provide L_[ O_III] / L_[ C_II] ratios at 0.3 - 0.9, suggesting that ionized gas represents a smaller fraction of the total gas reservoir in DSFGs, consistent with theoretical models of DSFGs as transitional systems between gas-rich, turbulent disks and more evolved, gas-poor galaxies. Expanding samples of DSFGs with [O_III] emission will be key to place this heterogeneous, poorly-understood galactic phase in its astrophysical context.
http://arxiv.org/abs/2601.09789 (pdf) (51)
Luminous Red Galaxies (LRGs) are among the most massive galaxies at any epoch, and lack ongoing star formation. As systems hosting most of the baryonic mass in the local Universe, they preserve imprints of the quenching mechanisms in the early Universe. We exploited the large BOSS and DESI spectroscopic datasets to perform the first homogeneous and continuous mapping of the evolution of stellar population properties of a complete sample of the most massive LRGs (log (M_*/M_⊙)> 11.5) at 0.15 < z < 0.8. By consistently fitting the same spectral indices at all redshifts, we measured trends of [Fe/H], [alpha/Fe], and light-weighted age as a function of redshift. These galaxies exhibit a passive light-weighted age evolution and flat [Fe/H] and [alpha/Fe] trends towards lower redshift, indicating genuinely passive evolution. These trends are robust against the choice of stellar population models and analysis assumptions, and they support the predictions from IllustrisTNG, which predict negligible chemical evolution for the most massive quenched systems at z < 0.8. Our results suggest that, despite nearly 5 Gyr of cosmic time and a 3-4x increase in number density, the stellar population properties of massive quiescent galaxies remain essentially unchanged since z 0.8. This shows a negligible progenitor bias below z 0.8, and a genuinely passive evolution. Newly added systems after z ∼ 0.8 were already largely quenched and chemically mature, while subsequent evolution was dominated by dry mergers without altering the bulk of the stellar populations.
http://arxiv.org/abs/2601.10709 (pdf) (48)
The ongoing Euclid mission aims to measure spectroscopic redshifts for approximately two million galaxies using the H α line emission detected in near-infrared slitless spectroscopic data from the Euclid Deep Fields (EDFs). These measurements will reach a flux limit of 5× 10^-17 erg cm^-2 s^-1 in the redshift range 0.4<z<1.8, opening the door to numerous investigations involving galaxy evolution, extending well beyond the mission&aposs core objectives. The achieved H α luminosity depth will lead to a sufficiently high sampling, enabling the reconstruction of the large-scale galaxy environment. We assess the quality of the reconstruction of the galaxy cosmic web environment with the expected spectroscopic dataset in EDFs. The analysis is carried out on the Flagship and GAEA galaxy mock catalogues. The quality of the reconstruction is first evaluated using geometrical and topological statistics measured on the cosmic web, namely the length of filaments, the area of walls, the volume of voids, and its connectivity and multiplicity. We then quantify how accurately gradients in galaxy properties with distance from filaments can be recovered. As expected, the small-scale redshift-space distortions, have a strong impact on filament lengths and connectivity, but can be mitigated by compressing galaxy groups before skeleton extraction. The cosmic web reconstruction is biased when relying solely on H α emitters. This limitation can be mitigated by applying stellar mass weighting during the reconstruction. However, this approach introduces non-trivial biases that need to be accounted for when comparing to theoretical predictions. Redshift uncertainties pose the greatest challenge in recovering the expected dependence of galaxy properties, though the well-established stellar mass transverse gradients towards filaments can still be observed.
http://arxiv.org/abs/2601.09778 (pdf) (33)
JWST&aposs "Little Red Dots" (LRDs) are increasingly interpreted as active galactic nuclei (AGN) obscured by dense thermalized gas rather than dust as evidenced by their X-ray weakness, blackbody-like continua, and Balmer line profiles. A key question is how LRDs connect to standard UV-luminous AGN and whether transitional phases exist and if they are observable. We present the "X-Ray Dot" (XRD), a compact source at z=3.28 observed by the NIRSpec WIDE GTO survey. The XRD exhibits LRD hallmarks: a blackbody-like (T_eff≃ 6400K) red continuum, a faint but blue rest-UV excess, falling mid-IR emission, and broad Balmer lines (FWHM ∼ 2700-3200 km s^-1). Unlike LRDs, however, it is remarkably X-ray luminous (L_2-10keV = 10^44.18ergs^-1) and has a continuum inflection that is bluewards of the Balmer limit. We find that the red rest-optical and blue mid-IR continuum cannot be reproduced by standard dust-attenuated AGN models without invoking extremely steep extinction curves, nor can the weak mid-IR emission be reconciled with well-established X-ray–torus scaling relations. We therefore consider an alternative scenario: the XRD may be an LRD in transition, where the gas envelope dominates the optical continuum but optically thin sightlines allow X-rays to escape. The XRD may thus provide a physical link between LRDs and standard AGN, offering direct evidence that LRDs are powered by supermassive black holes and providing insight into their accretion properties.
http://arxiv.org/abs/2601.10097 (pdf) (32)
We report the discovery of two swallowtail image formations at z=2.91 and z=6.70 behind the galaxy cluster MS 0451.6-0305 in JWST-NIRCam imaging. We find that in both of the above lensed systems, the complex image morphology cannot be reproduced by simple fold/cusp caustics, and detailed lens modeling reveals higher-order swallowtail caustic configurations. In the z=2.91 lens system, a small part of the source galaxy (which itself is part of a galaxy group) containing atleast two compact knots sits inside the swallowtail caustic, producing a quadruply imaged arc. At two of the image positions of these knots, we infer point source magnifications of ≳ 300, implying lensing-corrected effective radii of ≲ 0.8-1.5 pc. The z=6.70 system exhibits even more complex image formation. We therefore only use the most confidently identified counter-images of knots in this system as constraints in our lens modeling. The resulting model predicts magnifications ∼20-200 and lensing-corrected effective radii of ≲ 0.8-18.5 pc for various knots. Together, these two systems represent the first example of observations of multiple swallowtail image formations in a single galaxy cluster and demonstrate the ability of swallowtail caustics to magnify individual substructures at sub-parsec scales, from intermediate redshifts to the first billion years of the Universe.
http://arxiv.org/abs/2601.10573 (pdf) (31)
We present a statistical study on the origins of the UV continuum and narrow/broad emission lines in little red dots (LRDs), a newly discovered class of active galactic nuclei (AGNs). Leveraging all archived JWST/NIRSpec data, we build a sample of 28 spectroscopically-confirmed LRDs at 5<z_spec<7.2, by requiring broad Hα emission, blue UV colors, V-shaped continua, and compact morphologies. We define a control sample of 9 blue, compact, broad-line AGNs without red optical continua (hereafter little blue dots; LBDs), and examine correlations between rest UV and the narrow/broad Hα luminosities in these populations. In LRDs, both narrow and broad Hα components are tightly correlated with the UV continuum, and the luminosity ratios are consistent with those in young starburst galaxies. In contrast, the UV to broad Hα ratios in LBDs closely match local unobscured AGNs and are statistically different from LRDs. The Lyα occurrence rates and strengths do not differ between LRDs and LBDs and are comparable to normal star-forming galaxies. These results are consistent with a scenario where the central BH in LRDs is enshrouded by a dense opaque gas envelope – in this model, the UV continuum as well as narrow and even broad Hα emissions are not powered by AGNs but predominantly by young massive stars surrounding the envelope, while the envelope radiates as a ∼ 5000 K blackbody. As the envelope dissipates, direct AGN emission can emerge, potentially transforming LRDs into LBDs and marking the end of a short-lived phase of rapid black hole growth.
http://arxiv.org/abs/2601.10361 (pdf) (29)
The bending of radio galaxies in galaxy clusters is expected to be caused by interactions with the local environment. The physical processes responsible for jet bending, and their influence on the polarization properties of radio galaxies, remain poorly understood, leading to the question of whether jet properties in bent radio galaxies differ from those in linear radio galaxies. Using a sample of 24 polarized bent radio galaxies, observed with the Karl G. Jansky Very Large Array at 1–2 GHz, we test for correlation of bending angle with polarization parameters measuring Faraday rotation, intrinsic fractional polarization, and Faraday rotation dispersion, used here as a measure of turbulence along the line of sight. We find no statistically significant correlations. At the spatial resolution of our dataset (3–46 kpc, median 18.4 kpc), our results indicate that we are primarily probing larger-scale intracluster medium effects not related to bending angle. The absence of a statistically significant correlation suggests that bent radio galaxies are reliable probes of intracluster magnetic fields, because their intrinsic properties do not appear to introduce systematic biases into measured polarization parameters. We do detect a preference for source magnetic field vectors to align with the direction of jet bending. Finally, we estimate that the POSSUM and SKA surveys will contain ≳300 and ≳1000 polarized radio galaxies, respectively, providing large future samples with a range of bending angles and similar redshift distribution and number of beams per source as in our sample, enabling our results to be tested with greater statistical power.
http://arxiv.org/abs/2601.09810 (pdf) (28)
Luminous infrared galaxies are key sites of obscured stellar mass assembly at z > 0.5. Their star formation rates (SFRs) are often estimated using the luminosities of the 6.2 micron and 11.2 micron polycyclic aromatic hydrocarbon (PAH) features, or those of the [Ne II] and [Ne III] fine-structure lines, as they are minimally affected by obscuration. It is uncertain whether the calibration of these features as SFR tracers depends on the starburst bolometric luminosity or the level of Active Galactic Nucleus (AGN) activity. We here investigate the relationship between the luminosities of PAH and Neon lines with star formation rate for highly luminous objects using radiative transfer modeling and archival observations of 42 local Ultraluminous (>= 10^12 L_sun) Infrared Galaxies (ULIRGs). We find that PAH and [Ne II] features arise mainly in star-forming regions, with small contributions from the AGN or host, but that the [Ne III] line has a mixed contribution from both star formation and AGN activity. We present relations between L_PAH and L_NeII, and both starburst luminosity and SFR. We find relations for lower luminosity (L_IR = 10^10-10^12 L_sun) systems underestimate the SFRs in local ULIRGs by up to 1 dex. The 6.2 micron and 11.2 micron PAH features, and the [Ne II] line, are thus good tracers of SFR in ULIRGs. We do not find that a more luminous AGN affects the relationship between SFR and PAH or Neon luminosity, but that it can make PAH emission harder to discern. Our results and derived relations are relevant to studies of star-forming and composite galaxies at z < 3 with the James Webb Space Telescope.
http://arxiv.org/abs/2601.10640 (pdf) (26)
We present a four-component high-resolution model of thermal dust emission for microwave and sub-mm frequencies derived from Planck HFI, WHAM and Gaia. The resulting high-resolution model derived here forms the basis for the thermal dust model employed in the Cosmoglobe DR2 reanalysis of COBE-DIRBE. The four dust components are called “cold dust”, “hot dust”, “nearby dust”, and “Ha correlated dust”, respectively, and trace different physical environments. The spatial distributions of the nearby dust and Ha dust components are defined by the Edenhofer et al. Gaia 3D extinction model and the WHAM survey, respectively, while the hot and cold dust components are fit freely pixel-by-pixel to the Planck HFI data. We use a global parameter grid search coupled to an amplitude map Gibbs sampler to fit this model to Planck HFI data. In agreement with the companion low-resolution analysis, we find that the hot dust component is strongly correlated with the FIRAS Cii map, while the cold dust component is strongly correlated with the HI4PI Hi map. Despite its fewer degrees of freedom per pixel compared to the Planck 2015 legacy dust model, we find that this new model performs competitively in terms of overall residuals, capturing over 98% of the full-sky dust variance for all channels. When fitting a spatially varying 3-parameter MBB model to the new dust model with isotropic SEDs, we find very similar spatial distributions to those of the official Planck analysis, and this new model thus represents an economical decomposition of previously published spatially varying spectral parameter maps. We conclude that this new model represents both a statistically more efficient summary of thermal dust in the microwave and far-infrared regimes and a physically more realistic decomposition of the sky compared to the traditional 3-parameter MBB model. (abridged)
http://arxiv.org/abs/2601.10304 (pdf) (25)
The 21-cm global signal, a treasure trove of information about the nature of the first luminous sources of the Universe, has traditionally been modelled assuming that these early sources were predominantly star-forming galaxies. However, recent observations by the James Webb Space Telescope (JWST) have revealed several AGNs as early as z 10 - 10.4 . In light of this, it is important to investigate the contribution of such AGNs to the 21-cm signal. Assuming that these AGNs are seeded by Primordial Black Holes (PBHs) and employing an analytical PBH model, consistent with existing cosmological and astrophysical constraints, we show that these exotic objects can have a significant impact on the redshift evolution of the global signal.
http://arxiv.org/abs/2601.10635 (pdf) (24)
We present spatially resolved spectroscopic observations of the planetary nebula NGC 6720, the Ring Nebula, taken during the science verification phase of WEAVE, a new instrument mounted on the William Herschel Telescope on La Palma. We use the instrument&aposs Large Integral Field Unit (LIFU) to obtain spectra of the Ring Nebula, covering its entire optically bright inner regions as well as parts of its much fainter outer molecular halo. We report the discovery of emission from [Fe v] and [Fe vi] confined to a narrow “bar” extending across the central regions of the nebula. No lines of other elements share this morphology or, at the spectral resolving power used (R ∼ 2500), the same radial velocity. The extent to which iron in this bar is depleted is presently unclear; comparison with JWST-detected dust continuum emission suggests that some dust grain destruction may be occurring in the region, but there is currently no observational evidence for the > 50 km s^-1 shock waves or T > 10^6 K X-ray emitting gas needed to enable this. Where the bar is located along the line of sight through the nebula, and how it was created, are new puzzles to be solved for this iconic planetary nebula.
http://arxiv.org/abs/2601.10514 (pdf) (22)
A substantial fraction of metal-poor stars in the local Milky Way halo exhibit large overabundances of carbon. These stars, dubbed Carbon-Enhanced Metal-Poor (CEMP) stars, provide crucial constraints on the nature of the early universe including the earliest nucleosynthetic events. Whether these stars exist at similar rates in nearby galaxies is a major open question with implications for the environmental dependence of early chemical evolution. Here, we present the discovery of the first five CEMP stars in the Milky Way&aposs largest dwarf companion, the LMC, using SDSS-V spectra from the BOSS instrument. We measure metallicities ranging from [Fe/H] = -2.1 to -3.2 and evolutionary state corrected carbon enhancements of [C/Fe] = +1.2 to +2.4, placing these stars among the most metal-poor and carbon-rich ever identified in the LMC. This discovery demonstrates that CEMP stars do exist in the LMC despite previous null detections, and establishes the foundation for measuring the CEMP occurrence rate in this system. Such measurements will provide critical tests of whether environmental differences affect the formation channels and frequencies of these ancient, carbon-rich stars.
http://arxiv.org/abs/2601.10642 (pdf) (21)
The study is devoted to comparing the parameters of the interstellar medium of HII regions in the Kongiu and Groves catalogs for the galaxy NGC 628. The article analyzes the characteristics of star-forming regions, including a comparison of radiation fluxes in the ranges of 7.7 μm and 21 μm and in the Hα, Hβ, OIII and CO lines, calculating the kinematic parameters (FWHM) for the lines, and analyzing the spatial distribution of regions for both catalogs. The results of the study showed that the regions from the Groves catalog demonstrate higher line widths compared to the Kongiu catalog. Signs of possible misidentified classification of some regions from the Groves catalog were revealed: there is a possibility that some of them are not HII regions, but shock ionization regions.
http://arxiv.org/abs/2601.10579 (pdf) (21)
We present new Hubble Space Telescope (HST) UV spectra of the K2 V star HD 166620, the first star clearly recognized to be in a "magnetic grand minimum" state analogous to the Sun&aposs "Maunder Minimum" in the late 1600&aposs. The stellar H I Lyman-alpha surface fluxes are extremely low, about a factor of two below fluxes observed during solar minimum, and also significantly lower than those of Tau Ceti (G8 V) and HD 191408 (K2.5 V), two stars more similar to HD 166620 in spectral type and age ( 10 Gyr) than the Sun. The Tau Ceti data that are compared with HD 166620 include both old archival data and a new HST observation as well. The Lyman alpha data are used to confirm a nondetection of astrospheric Lyman-alpha absorption for this star, suggesting a very weak wind with Mdot<0.1 Mdot_sun. The very compact astrosphere inferred for Tau Ceti indicates that the star&aposs debris disk is at least partly exposed to the ISM, and we discuss possible consequences.
http://arxiv.org/abs/2601.10240 (pdf) (18)
To date, Galactic Astronomy has largely concerned itself with astrophysical processes, and with the locations, space motions and compositions of objects. Consider, for example, the elucidation of the components of the Galaxy over the past decades, its mapping as enabled by Gaia and its predecessors, the photometric and spectroscopic characterization of innumerable astrophysical objects in various wavelength ranges, both from the ground and from space, and the expanding discovery and characterization of exoplanets; all focused on the current, static Galaxy. This White Paper proposes a dedicated program to derive stellar ages from time-series spectroscopy to hasten the transformation of this static conception into a dynamical one with age-labeled objects and events.
http://arxiv.org/abs/2601.10595 (pdf) (15)
The scarcity of high-fidelity extinction measurements remains a bottleneck in deriving accurate stellar properties from Gaia parallaxes. In this work, we aim to derive precision extinction estimates for APOGEE DR19 stars, establishing a new benchmark for Galactic stellar population studies. We first determine reddening by comparing observed colorsr, etrieved from photometric surveys or standardized synthetic magnitudes from Gaia BP/RP spectra, to intrinsic colors predicted via an XGBoost model. The model is trained on minimally reddened stars to infer intrinsic colors and their associated uncertainties, using APOGEE stellar parameters (Teff, logg, [Fe/H], and [alpha/Fe]). The derived reddening values are then converted into extinctions using an anchor ratio of A_BP / A_RP = 1.694 +/- 0.004, derived from red-clump-like stars. Here, we provide extinction measurements in 39 filters across 10 photometric systems and introduce a new empirical extinction curve optimized for broadband passbands. Our extinction estimates (Av) outperform existing results (Bayestar19, StarHorse, SEDEX), achieving a typical precision of 0.03 mag in Av. Notably, we identify systematic deviations of up to 30% between monochromatic and passband-integrated extinction ratios at wavelengths greater than 700 nm. This result highlights the necessity of adopting passband-specific coefficients when correcting extinction to derive stellar parameters. As the foundation for a forthcoming series of papers, these benchmark measurements will be used to (1) revise asteroseismic scaling relations, (2) calibrate differential reddening in open clusters, and (3) reconcile heterogeneous dust maps into a unified, all-sky extinction scheme.
http://arxiv.org/abs/2601.09796 (pdf) (14)
Milky Way halo substructures identified in dynamical space are known to suffer from contamination from the Milky Way in-situ stars, which makes their accreted origins uncertain. We present detailed chemical abundances of 35 stars belonging to two sets of dynamically tagged groups, Rg8 and Rg9, to investigate their accreted nature. Both groups are composed of stars with low orbital energy and very retrograde orbits. We find that Rg8 and Rg9 are chemically indistinguishable across all elements, from C to Eu, strongly indicating that they belong to the same structure. The iron-abundance distribution of this low-E retrograde group has a prominent peak at [Fe/H] ≈-2.1, revealing that its main population is very metal-poor, and a secondary peak at [Fe/H] ≈-1.5, very likely due to contamination from Milky Way in-situ stars. These groups also heavily overlap with the Thamnos substructure in dynamical space, and we thus use them to investigate the chemical properties of Thamnos. The dominant, low-metallicity population provides strong evidence for the ex-situ origin of Thamnos, as well as its very metal-poor nature. We do not see any evidence of an α knee in our sample, which is consistent with previous studies. Comparison with the Cetus-Palca stream in the chemical space shows similar abundance distributions, and thus it suggests that the Thamnos progenitor dwarf galaxy had a truncated star formation history due to its early merger with the Milky Way.
http://arxiv.org/abs/2601.10404 (pdf) (13)
We analyzed optical/near-IR Large Binocular Telescope spectra of 16 sources alerted by Gaia between 2021 and 2024 due to significant photometric variability. Half of the spectra were taken during quiescence and the rest during a burst or at intermediate brightness. Our analysis of their ten-year light curves and photometric/spectroscopic features provide evidence that all 16 sources are accreting Young Stellar Objects (YSOs). One object, Gaia23bab, is a known EXor source. Other light curves either have peaks over a stable baseline, or significant variability throughout the entire observation period, suggesting multiple contributing processes. All spectra exhibit emission lines from accretion columns, and over half of them show atomic forbidden lines as signatures of outflowing gas. We determined stellar parameters, accretion luminosity (Lacc) and mass accretion rate (Macc) at different brightness phases. Only two sources showed variability primarily due to extinction. During quiescence, our sources exhibit Lacc and Macc values typical of T Tauri and Herbig Ae/Be (HAEBE) sources, supporting the hypothesis that any YSO may undergo episodic accretion. In bursts, the Lacc and Macc of sources with photometric variations exceeding 2 mag follow a shallower relation with stellar luminosity and mass, typical of known EXor sources. This group includes one Class I, one flat-spectrum, and two Class II sources. Notably, the other Class I source, Gaia24beh, shows an Lacc value about ten times higher than typical EXor bursts of the same mass. In the other cases, Lacc and Macc align with variability seen in T Tauri and HAEBE sources.
http://arxiv.org/abs/2601.09807 (pdf) (13)
One of the earliest mechanisms proposed for plasma heating was magnetic pumping (MP). However, its significance for astrophysical phenomena, including particle acceleration, has yet to be appreciated. MP-energized particles tap energy from magnetic-field oscillations. A particle&aposs momentum component perpendicular to the local B-field increases during field growth by virtue of the adiabatic invariant p_⊥^2/B=const. The gained p_⊥ is then partially scattered elastically into the parallel momentum, p_∥, with p^2=p_∥^2+p_⊥^2=const, thereby retaining some fraction of the gained energy before the field decreases to its minimum. This scattering breaks the reversibility of energy exchange between particles and oscillating magnetic fields, thereby increasing the particle energy after each MP cycle. Field oscillations are often assumed to be sinusoidal, and the resulting MP is treated perturbatively. These simplifications restrict astrophysical applications, leaving objects with strong magnetic perturbations outside the scope of adequate treatment. We develop a nonperturbative approach to MP that is suitable for a broad spectrum of magnetic turbulence. The treatment comprises two steps. The first step is common: converting a kinetic equation into an infinite hierarchy of moments of the particle distribution function. The second step is new in MP treatments: we find an exact closure at an arbitrary level of the moment system. The heating is treated exactly at the second-moment closure. Particle acceleration generally requires a higher-level closure to determine the power-law index and the maximum energy of accelerated particles. We propose a method for extracting these crucial acceleration data from the second moment for a broadband random field.
http://arxiv.org/abs/2401.02476 (pdf) (13)
Milky Way globular clusters (GCs) display chemical enrichment in a phenomenon called multiple stellar populations (MSPs). While the enrichment mechanism is not fully understood, there is a correlation between a cluster&aposs mass and the fraction of enriched stars found therein. However, present-day GC masses are often smaller than their masses at the time of formation due to dynamical mass loss. In this work, we explore the relationship between mass and MSPs using the stellar stream 300S. We present the chemical abundances of eight red giant branch member stars in 300S with high-resolution spectroscopy from Magellan/MIKE. We identify one enriched star characteristic of MSPs and no detectable metallicity dispersion, confirming that the progenitor of 300S was a globular cluster. The fraction of enriched stars (12.5
http://arxiv.org/abs/2601.10314 (pdf) (13)
We present the implementation of a Self-Calibration of Intrinsic Alignments of galaxies as an extension to the Vera C. Rubin Observatory&aposs Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC)&aposs weak lensing 3x2pt pipeline (TXPipe). As a demonstration, we have run this pipeline on the Dark Energy Survey (DES) year one data set. We find indications of a non-zero intrinsic alignment signal in the higher redshift bins, while in the lower bins our results look more uncertain. We believe this is caused by known issues with the individual galaxies photo-z estimation. This effect is particularly harmful for the self-calibration method, since it has high requirements for reliable estimation of the photo-zs, and the need for individual galaxy point estimates and tomographic binning to match. We show how different methods of recreating the redshift probability distribution can affect the detection of intrinsic alignment.
http://arxiv.org/abs/2601.10598 (pdf) (12)
The Galaxy&aposs most common known planetary systems have several Earth-to-Neptune-size planets in compact orbits. At small orbital separations, larger planets are less common than their smaller counterparts by an order of magnitude. The young star V1298 Tau hosts one such compact planetary system, albeit with four planets that are uncommonly large (5 to 10 Earth radii). The planets form a chain of near-resonances that result in transit-timing variations of several hours. Here we present a multi-year campaign to characterize this system with transit-timing variations, a method insensitive to the intense magnetic activity of the star. Through targeted observations, we first resolved the previously unknown orbital period of the outermost planet. The full 9-year baseline from these and archival data then enabled robust determination of the masses and orbital parameters for all four planets. We find the planets have low, sub-Neptune masses and nearly circular orbits, implying a dynamically tranquil history. Their low masses and large radii indicate that the inner planets underwent a period of rapid cooling immediately after dispersal of the protoplanetary disk. Still, they are much less dense than mature planets of comparable size. We predict the planets will contract to 1.5-4.0 Earth radii and join the population of super-Earths and sub-Neptunes that nature produces in abundance.
http://arxiv.org/abs/2601.10637 (pdf) (11)
In recent years, a two-dimensional (2D) coronal seismology technique applied to spectral-imaging data from the Coronal Multi-channel Polarimeter (CoMP) and UCoMP has enabled routine measurement of the global coronal magnetic field. The technique combines coronal transverse wave phase speed from Doppler measurements with electron densities from the Fe xiii 10798/10747 Å intensity ratio to infer the magnetic field strength, while the wave propagation directions from Doppler measurements trace the magnetic field direction. To validate the accuracy and robustness of this method, we use forward modeling of a MURaM simulation that produces open and closed magnetic structures with excited waves. From the synthetic Doppler velocity, Fe xiii infrared line intensities, and linear polarization signals, we apply the 2D coronal seismology technique to estimate the magnetic field strength and direction. A comparison with the simulation ground truth shows close agreement, indicating that the technique can recover the line-of-sight emissivity-weighted magnetic field direction and strength with high accuracy. We also perform a parameter-space analysis to quantify sensitivities of the method to parameter choice. These findings provide practical guidance for CoMP/UCoMP-like analysis and demonstrate that 2D coronal seismology can deliver reliable, LOS emissivity-weighted measurements of the coronal magnetic field from coronal wave observations.
http://arxiv.org/abs/2601.10631 (pdf) (10)
Atacama Large Millimeter/submillimeter Array observations suggest that the disc of HD 163296 is being actively shaped by embedded, yet unseen protoplanets, as indicated by numerous gas and dust substructures consistent with planet-disc interaction models. We report the first detection of simple organic molecules, HCN and C2H, tracing a candidate circumplanetary disc (CPD) in the HD 163296 system, located at an orbital radius of R=88±7 au and azimuth ϕ=46±3^∘ (or R=0.75”, PA=350^∘ in projected sky coordinates), and originating near the midplane of the circumstellar disc. The signature is localised but spectrally resolved, and it overlaps with a previously reported planet candidate, P94, identified through kinematic perturbations traced by CO lines. We propose a scenario in which the observed chemical anomalies arise from increased heating driven by the forming planet and ongoing accretion through its CPD, facilitating the thermal desorption of species that would otherwise remain frozen out in the disc midplane, and potentially triggering the activation barriers of chemical reactions that lead to enhanced molecular production. Based on a first-order dynamical analysis of the HCN spectrum from the CPD–isolated with a 7σ significance–we infer an upper limit on the planet mass of 1.8 M_Jup, consistent with predictions from CO kinematics and constraints from direct imaging studies. By comparing the CPD sizes derived from our models with theoretical expectations where the CPD radius corresponds to roughly one-third of the planet&aposs Hill radius, we favor CPD gas temperatures T > 150 K, planet masses M_p < 1.0 M_Jup, and CPD radii R_CPD < 2 au.
http://arxiv.org/abs/2601.10121 (pdf) (10)
A central goal in supernova (SN) research is to identify and characterize their progenitors. However, this is very difficult due to the limited archival images with sufficient depth and spatial resolution required for direct progenitor detection and due to the circumstellar dust which often biases the estimate of their intrinsic parameters. This field will be revolutionized by Euclid and the upcoming Chinese Space Station Survey Telescope (CSST), which conduct deep, wide-field, high-resolution and multi-band imaging surveys. We analyze their detection capability by comparing the model magnitudes of red supergiant (RSG) progenitors with the detection limits under different conditions, and we estimate the annual detection rates with Monte-Carlo simulations. We explore how to recover the intrinsic properties of SN progenitors with the help of radiation transfer calculations in circumstellar dust. We find the optical and near-infrared filters of the Euclid and CSST are highly effective for detecting RSG progenitors. We predict that archival images from the completed 2 surveys will enable ≲13 (or 24) progenitor detections per year within the mass range of 8–16 (or 8–25)M_⊙, an order of magnitude higher than the current detection rate of ∼1 detection per year. In the presence of circumstellar dust, the emerging spectral energy distribution (SED) of the progenitor is mainly affected by the optical depth and is almost independent of dust temperature in the Euclid and CSST filters. Our mock tests demonstrate that one can derive the progenitor mass and dust optical depth simultaneously by fitting the observed SED over the 11 filters of the 2 surveys while fixing the dust temperature to a typical value. Euclid and CSST will significantly enlarge the sample of direct progenitor detections with accurate mass measurements, which is crucial to resolve the long-standing RSG problem.
http://arxiv.org/abs/2601.09846 (pdf) (10)
Recently Quenched Ellipticals (RQEs) provide a valuable test case for disentangling intrinsic and environmental quenching, particularly because they are commonly classified as isolated central galaxies in low-mass halos. However, central/satellite assignments and isolation labels can vary across group catalogs, and such misclassifications can strongly bias physical interpretations. We present a uniform, physically motivated reassessment of the environments of 155 RQEs previously identified as centrals in an SDSS-based group catalog. We construct value-added neighbor catalogs via KD-tree searches and apply consistent thresholds in projected separation and line-of-sight velocity to (i) verify centrality, (ii) quantify isolation using a mass-ratio-based companion criterion, and (iii) identify potential pseudo-centrals via proximity to massive clusters. We find that 132/155 (85.2
http://arxiv.org/abs/2601.09768 (pdf) (10)
In data-driven scientific discovery, a challenge lies in classifying well-characterized phenomena while identifying novel anomalies. Current semi-supervised clustering algorithms do not always fully address this duality, often assuming that supervisory signals are globally representative. Consequently, methods often enforce rigid constraints that suppress unanticipated patterns or require a pre-specified number of clusters, rendering them ineffective for genuine novelty detection. To bridge this gap, we introduce CLiMB (CLustering in Multiphase Boundaries), a domain-informed framework decoupling the exploitation of prior knowledge from the exploration of unknown structures. Using a sequential two-phase approach, CLiMB first anchors known clusters using constrained partitioning, and subsequently applies density-based clustering to residual data to reveal arbitrary topologies. We demonstrate this framework on RR Lyrae stars data from the Gaia Data Release 3. CLiMB attains an Adjusted Rand Index of 0.829 with 90% seed coverage in recovering known Milky Way substructures, drastically outperforming heuristic and constraint-based baselines, which stagnate below 0.20. Furthermore, sensitivity analysis confirms CLiMB&aposs superior data efficiency, showing monotonic improvement as knowledge increases. Finally, the framework successfully isolates three dynamical features (Shiva, Shakti, and the Galactic Disk) in the unlabelled field, validating its potential for scientific discovery.
http://arxiv.org/abs/2601.10550 (pdf) (9)
The remnants of binary neutron star mergers are expected to be massive, rapidly rotating stars whose oscillations produce gravitational waves in the kilohertz band. The degree of differential rotation and the rotation profiles strongly influence their structure, stability and oscillation spectrum, and must therefore be taken into account when modeling their dynamics. We extend the pseudospectral code ROXAS (Relativistic Oscillations of non-aXisymmetric neutron stArS) to enable the dynamical evolution of oscillating, differentially rotating neutron stars. Using the updated code, we aim to study the star&aposs oscillation frequencies. We extend the previous formalism, based on primitive variables and the conformal flatness approximation, to differential rotation. Within this framework, we run a series of axisymmetric and non-axisymmetric simulations of perturbed, differentially rotating neutron stars with different rotation rates, and extract their oscillation frequencies. Axisymmetric modes, as well as those under the Cowling approximation, show excellent agreement with published results. We show that the secondary fundamental mode in the Cowling approximation is an artifact that does not appear in dynamical spacetimes. In addition, we provide, for the first time, frequency values for non-axisymmetric modes in differentially rotating configurations evolved in conformal flatness. This extension broadens the range of physical scenarios that can be studied with ROXAS, and represents a step toward more realistic modeling of post-merger remnants and their gravitational-wave emission.
http://arxiv.org/abs/2601.10612 (pdf) (9)
Observations of active galactic nuclei have revealed a correlation between millimeter luminosity, X-ray luminosity, and mass, suggesting the emission in each of these bands is powered by a common source. Starting with a set of five general relativistic magnetohydrodynamic simulations with dynamically important magnetic fields, we perform ray-tracing calculations to produce spectra including synchrotron emission, bremsstrahlung emission, and Compton scattering. Our models with similar Eddington ratios to the objects for which the relationship was inferred naturally reproduce observations without tuning. Our lower Eddington ratio models depart from this relationship, likely attributable to an observational bias against extremely low accretion rates. We find that inverse Compton scattering dominates the production of X-rays over bremsstrahlung radiation in almost all models, and in all models consistent with the observed correlation. We find only a modest spin dependence in this relationship. This study demonstrates that a compact, hot accretion flow with dynamically important magnetic fields can naturally explain observed millimeter and X-ray properties in low-luminosity active galactic nuclei. Future work should explore the impacts of non-thermal electron populations, weaker magnetic fields, and radiative cooling.
http://arxiv.org/abs/2601.10283 (pdf) (8)
Helioseismic and recent asteroseismic observations reveal fine structure in the power spectrum with alternating peaks and troughs above the acoustic cut-off frequency. This structure is interpreted as the interference patterns of high-frequency acoustic waves excited in the solar interior and propagating into the atmosphere, known as pseudomodes. Pseudomodes exhibit clear solar-cycle variability, with frequency shifts that occur predominantly in anti-phase with the activity cycle, although the underlying mechanism remains uncertain. This work investigates how the subsurface excitation source location and the photospheric acoustic cut-off frequency influence the formation, frequency distribution, and solar-cycle variability of pseudomodes. We employ an analytical Klein-Gordon subsurface cavity model, which is shown to act as an effective Fabry-Pérot interferometer for high-frequency waves that experience constructive and destructive interference between the source location and the lower turning point. We derive an effective dispersion relation isolating the effects of the source location and photospheric cut-off on the pseudomode frequency. The model reproduces the observed peak-trough pseudomode spectrum for reasonable parameter values constrained by Bayesian MCMC best-fitting to GONG observations. We also find that solar-cycle-associated 11-year modulations of the source location result in anti-phase pseudomode frequency shifts, whereas similar cyclic variations in the cut-off frequency produce harmonic-dependent behaviour, yielding both in-phase and anti-phase shifts. As the acoustic cut-off and mode excitation relate to stratification and flows in the solar interior, the results highlight pseudomodes as a powerful diagnostic tool for changes in subsurface solar and stellar structure through the solar cycle.
http://arxiv.org/abs/2601.10152 (pdf) (7)
Building on the transformative success of optical redshift surveys, the emerging technique of neutral hydrogen (HI) intensity mapping (IM) offers a novel probe of large-scale structure (LSS) growth and the late-time accelerated expansion of the universe. We present cosmological forecasts for the Baryon Acoustic Oscillations from Integrated Neutral Gas Observations (BINGO), a pioneering HI IM experiment, quantifying its potential to constrain the Planck-calibrated ΛCDM cosmology and extensions to the w_0w_aCDM dark energy model. For BINGO&aposs Phase 1 configuration, we simulate the HI IM signal using a lognormal model and incorporate three dominant systematics: foreground residuals, thermal noise, and beam resolution effects. Using Bayesian inference, we derive joint constraints on six cosmological parameters (Ω_b h^2, Ω_c h^2, 100θ_s, n_s, ln 10^10 A_s, and τ_r) alongside 60 HI parameters (b_HI^i, Ω_HI^i b_HI^i) across 30 frequency channels. Our results demonstrate that combining BINGO with the Planck 2018 CMB dataset tightens the confidence regions of cosmological parameters to ∼40 These results demonstrate BINGO&aposs potential to extract significant cosmological information from the HI distribution and provide constraints competitive with current and future cosmological surveys.
http://arxiv.org/abs/2601.10038 (pdf) (6)
Artificial intelligence is rapidly transforming astronomical research, yet the scientific community has largely treated this transformation as an engineering challenge rather than an epistemological one. This perspective article argues that philosophy of science offers essential tools for navigating AI&aposs integration into astronomy–conceptual clarity about what "understanding" means, critical examination of assumptions about data and discovery, and frameworks for evaluating AI&aposs roles across different research contexts. Drawing on an interdisciplinary workshop convening astronomers, philosophers, and computer scientists, we identify several tensions. First, the narrative that AI will "derive fundamental physics" from data misconstrues contemporary astronomy as equation-derivation rather than the observation-driven enterprise it is. Second, scientific understanding involves more than prediction–it requires narrative construction, contextual judgment, and communicative achievement that current AI architectures struggle to provide. Third, because narrative and judgment matter, human peer review remains essential–yet AI-generated content flooding the literature threatens our capacity to identify genuine insight. Fourth, while AI excels at well-defined problem-solving, the ill-defined problem-finding that drives breakthroughs appears to require capacities beyond pattern recognition. Fifth, as AI accelerates what is feasible, pursuitworthiness criteria risk shifting toward what AI makes easy rather than what is genuinely important. We propose "pragmatic understanding" as a framework for integration–recognizing AI as a tool that extends human cognition while requiring new norms for validation and epistemic evaluation. Engaging with these questions now may help the community shape the transformation rather than merely react to it.
http://arxiv.org/abs/2601.09835 (pdf) (6)
Secular perturbations from binary stars and distant massive planets can drive cold planets onto nearly parabolic orbits with pericenter passages extremely close to their host stars. Meanwhile, short-period super-Earths are frequently observed around nearby stars. Gravitational scattering between these two distinct populations can lead to substantial orbital energy exchange, liberating some intruders from the gravitational confinement of their host systems. This process offers a robust formation channel for a subset of the abundant freely floating planet population. It may also significantly perturb the original orbits of close-in planets, induce collisional trajectories between close-in planets and their host stars, and disrupt the dynamical evolution of cold planets toward close stellar encounters.
http://arxiv.org/abs/2512.20171 (pdf) (6)
Cosmic voids in the large-scale structure are among the useful probes for testing gravity on cosmological scales. In this paper, we investigate the evolution of voids in the Horndeski theory using the effective field theory (EFT) of dark energy. Modeling the void formation with the dynamics of spherical mass shells, we study how modifications of gravity encoded into the EFT of dark energy change the linearly extrapolated critical density contrast that is relevant for the criterion for void formation, with particular focus on the time-dependent parameter characterizing the effect of kinetic braiding. It is found that the change in the critical density contrast is one order of magnitude smaller than the dimensionless EFT parameter because of a slight imbalance between two compensating effects. We then compute the void abundance using the Sheth–van de Weygaert void size function and demonstrate that it exhibits scale-dependent modifications. It is shown that the modifications to the void size function on small scales are almost entirely determined by the modified linear matter power spectrum, while the modifications on large scales are dominated by the contributions from the linear matter spectrum and the critical density contrast.