Trending in Astrophysics of Galaxies
Little Red Dots as Globular Clusters in Formation
Little Red Dots (LRDs), among the most enigmatic high-redshift discoveries by JWST, are commonly believed to be powered by accreting supermassive black holes. Here, we explore the possibility that these sources are globular clusters in formation, with rest-frame UV arising from a very young stellar population and rest-frame optical from a short-lived supermassive ($>10^4$ M$_\odot$) star. The spectral profiles of LRDs are broadly consistent with this scenario, though the observed temperatures and bolometric luminosities favor emission reprocessed by optically thick, continuum-driven winds not fully captured by current models. The LRD $z\sim5-7$ UV luminosity function naturally evolves, under standard evolutionary and mass-loss prescriptions, into a present-day mass function with a turnover at $\log_{10}(M_\ast$/$M_\odot)=5.3$ and an exponential cutoff at high masses, consistent with local globular-cluster populations. We estimate the total present-day number density of LRDs formed across all redshifts to be $\approx0.3$ Mpc$^{-3}$, similar to local globular clusters. The observed LRD redshift range matches the age distribution of metal-poor globular clusters, without current LRD counterparts to the metal-rich population. If LRDs are globular clusters in formation, we predict chemical abundance patterns characteristic of multiple stellar populations, including enhanced He and N, and potential Na-O and Al-Mg anti-correlations. These results offer a local perspective to explore this surprisingly abundant population of distant sources, and a potential new window into extreme stellar astrophysics in the early Universe.
2602.15935
Feb 2026Astrophysics of Galaxies
The large cores of dark matter and globular clusters in AS1063. Possible evidence of self-interacting dark matter. Or not
Deep JWST images of AS1063 reveals tens of thousands of globular clusters in the galaxy cluster AS1063. When compared with the lensing model based on the same JWST data, the distribution of globular clusters traces closely the distribution of lensing mass (mostly composed of dark matter). Interestingly, both the distributions of dark matter and globular clusters have large central cores. However the size of the core in the distribution of globular clusters is about half the size the core of the dark matter distribution. We argue that the standard cold dark matter and fuzzy dark matter models struggle to explain these large cores. Meanwhile, the self interacting dark matter with a velocity dependent cross section, combined with core stalling, offers a natural explanation to the existence of these cores if $σ_{\rm SI}\approx 0.3$ cm$^2$ g$^{-1}$ for galaxy cluster halos. But we also discuss how the lack of hydrodynamical N-body simulations capable of resolving globular clusters in galaxy cluster scale halos, hinders the possibility of ruling out the standard non-collisional dark matter scenario. Future high-resolution hydrodynamical simulations of galaxy clusters, with several trillion particles, and containing over a hundred thousand globular clusters, can provide the insight needed to transform the epistemic nature of dark matter into an ontological one
2602.15940
Feb 2026Astrophysics of Galaxies
Physics-Informed Neural Networks for Modeling Galactic Gravitational Potentials
We introduce a physics-informed neural framework for modeling static and time-dependent galactic gravitational potentials. The method combines data-driven learning with embedded physical constraints to capture complex, small-scale features while preserving global physical consistency. We quantify predictive uncertainty through a Bayesian framework, and model time evolution using a neural ODE approach. Applied to mock systems of varying complexity, the model achieves reconstruction errors at the sub-percent level ($0.14\%$ mean acceleration error) and improves dynamical consistency compared to analytic baselines. This method complements existing analytic methods, enabling physics-informed baseline potentials to be combined with neural residual fields to achieve both interpretable and accurate potential models.
2602.01806
Feb 2026Astrophysics of Galaxies
Variations in the Milky Way's Stellar Mass Function at [Fe/H] < -1
We present the first determination of the Galactic stellar mass function (MF) for low-mass stars (0.2-0.5 M_sun) at metallicities [Fe/H] < -1. A sample of ~53,000 stars was selected as metal-poor on the basis of both their halo-like orbits and their spectroscopic [Fe/H] from Gaia DR3 BP/RP (XP) spectra. These metallicity estimates for low-mass stars were enabled by calibrating Gaia XP spectra with stellar parameters from SDSS-V. For -1.5 < [Fe/H] < -1, we find that the MF below 0.5 M_sun exhibits a "bottom-heavy" power-law slope of alpha ~ -1.6. We tentatively find that at even lower metallicities, the MF becomes very bottom-light, with a near-flat power-law slope of alpha ~ 0 that implies a severe deficit of low-mass stars. This metallicity-dependent variation is insensitive to the adopted stellar evolution model. These results show that the Galactic low-mass MF is not universal, with variations in the metal-poor regime. A further calibration of XP metallicities in the regime of M < 0.5 M_sun and [Fe/H] < -1.5 will be essential to verify these tentative low-metallicity trends.
2601.19522
Jan 2026Astrophysics of Galaxies
Unveiling the 3D structure of the central molecular zone from stellar kinematics and photometry: The 50 and 20 km/s clouds
The central molecular zone (CMZ), surrounding the Galactic centre, is the largest reservoir of dense molecular gas in the Galaxy. Despite its relative proximity, the 3D structure of the CMZ remains poorly constrained, primarily due to projection effects. We aim to constrain the line-of-sight location of two molecular clouds in the CMZ -- the 50 and 20 km/s clouds -- and to investigate their possible physical connection using stellar kinematics and photometry. This study serves as a pilot for future applications across the full CMZ. We estimated the line-of-sight position of the clouds by analysing stellar kinematics, stellar densities, and stellar populations towards the cloud regions and a control field. We find an absence of westward moving stars in the cloud regions, which indicates that they lie on the near side of the CMZ. This interpretation is supported by the stellar density distributions. The similar behaviour observed in the two clouds, as well as in the region between them (the ridge), suggests that they are located at comparable distances and are physically linked. We also identified an intermediate-age stellar population (2-7 Gyr) in both regions, consistent with that observed on the near side of the CMZ. We estimated the line-of-sight distances at which the clouds and the ridge become kinematically detectable (i.e. where the proper motion component parallel to the Galactic plane differs from that of the control field at the 3 sigma level) by converting their measured proper motions parallel to the Galactic plane using a theoretical model of the stellar distribution. We find that the 50 and 20 km/s clouds are located at $43\pm8$ pc and $56\pm11$ pc from Sgr A*, respectively, and that the ridge lies at $56\pm11$ pc; this supports the idea that the clouds are physically connected through the ridge.
2601.05252
Jan 2026Astrophysics of Galaxies
The impact of cosmic voids on AGN activity
From the Eagle project, we study the properties of galaxies hosting AGN in cosmic voids and their surrounding structures, filaments and walls, at $z=0$, comparing them to non-AGN galaxies in similar environments. We found that the AGN fraction decreases as a function of void-centric distance, with void galaxies displaying the highest AGN fraction (12\%), and galaxies in denser environments, showing the lowest AGN fraction (6.7\%), consistent with observations. The AGN fraction is particularly high in most massive void galaxies when controlling for stellar mass. When comparing AGN host galaxies to inactive ones, we find that AGN galaxies tend to have slightly more massive SMBHs, higher specific star formation rates, and reside in higher-mass haloes at a given stellar mass than non-AGN galaxies. At $\rm M_{*} > \rm 10^{10.2} \rm M_{\odot}$, AGN hosts in voids tend to have slightly more massive SMBHs than those in denser environments. Otherwise, the AGN population does not show a clear trend in relation to the global environment. In contrast, non-AGN void galaxies host more massive SMBHs, slightly higher sSFRs, and are located in more massive haloes than those in denser environments. Analysing the recent merger histories of both AGN and non-AGN populations, we find that a larger fraction of massive AGN galaxies have undergone major mergers compared to non-AGN galaxies, regardless of environment. Notably, AGN galaxies in voids show a higher frequency of recent mergers, especially major mergers, than their counterparts in other environments, especially at high stellar mass. Our results suggest that the evolution of SMBHs in voids is closely related to that of their host galaxies and their surrounding environment, while the most recent AGN activity is more strongly linked to recent interactions.
2601.00594
Jan 2026Astrophysics of Galaxies
A New Empirical Fit to Galaxy Rotation Curves
We present a new empirical model for galaxy rotation curves that introduces a velocity correction term ω, derived from observed stellar motion and anchored to Keplerian baselines. Unlike parametric halo models or modified gravity theories, this approach does not alter Newtonian dynamics or invoke dark matter distributions. Instead, it identifies a repeatable kinematic offset that aligns with observed rotation profiles across a wide range of galaxies. Using SPARC data [1], we demonstrate that this model consistently achieves high fidelity fits, often outperforming MOND and CDM halo models in RMSE and R-squared metrics without parametric tuning. The method is reproducible, minimally dependent on mass modeling, and offers a streamlined alternative for characterizing galactic dynamics. While the velocity correction ω lacks a definitive physical interpretation, its empirical success invites further exploration. We position this model as a local kinematic tool rather than a cosmological framework, and we welcome dialogue on its implications for galactic structure and gravitational theory. Appendix B presents RMSE and R2 comparisons showing that this method consistently outperforms MOND and CDM halo models across a representative galaxy sample.
2601.00522
Jan 2026Astrophysics of Galaxies
A Spatially Resolved Evolutionary Sequence of Multi-wavelength AGN Host Galaxies
We study the spatially resolved star formation, gas ionisation, and outflow properties of 1813 active galactic nuclei (AGNs) from the MaNGA survey, which we classify into infrared (IR), broad-line (BL), narrow-line (NL), and radio (RD) AGNs based on their mid-infrared colours, optical spectra, and/or radio photometry. We also provide estimations of AGN power at different wavelengths. AGN incidence is found to increase with stellar mass following a power-law, with the high-mass end dominated by RDAGNs and the low-mass end dominated by NLAGNs. Compared to their mass-matched non-AGN counterparts, we find that IRAGNs, BLAGNs, and NLAGNs on average show enhanced specific star formation rates, younger stellar populations, and harder ionisation towards the centre. RDAGNs, in contrast, show radial profiles similar to quiescent galaxies. [OIII] outflows are more common and stronger in BL/IRAGNs, while RDAGNs on average show no outflow features. The outflow incidence increases with [OIII] luminosity, and the features in BL/IRAGNs on average extend to ~2 kpc from the nuclei. We further discuss a possible evolutionary sequence of AGNs and their host galaxies, where AGNs with strong emission lines or dust tori are present in star-forming galaxies. Later, young compact radio jets emerge, the host galaxies gradually quench, and the AGN hosts eventually evolve into globally quiescent systems with larger radio jets that prevent further gas cooling.
2512.11694
Dec 2025Astrophysics of Galaxies
GLIMPSE-D: An Exotic Balmer-Jump Object at z=6.20? Revisiting Photometric Selection and the Cosmic Abundance of Pop III Galaxies
We present deep JWST/NIRSpec G395M spectroscopy of GLIMPSE-16043, a promising $z\sim6$ Pop III candidate originally identified through NIRCam photometry as having weak [OIII]$λ\lambda4959,5007$ emission. Our follow-up reveals clear [OIII] emission, ruling out a genuine zero-metallicity nature. However, the combination of the measured line fluxes and photometry indicates that its spectral energy distribution requires an extraordinarily strong Balmer jump ($-1.66 \pm 0.47$ mag) and H$α$ equivalent width ($3750\pm1800$ Å), features that cannot be reproduced by current stellar+nebular or pure nebular photoionization models. The only models approaching the observations to almost within $1σ$ involve a hot ($T_{\rm eff}\!\simeq\!10^{4.7}$ K) single blackbody embedded in a low-$T_{\rm e}$ nebular environment, suggestive of scenarios such as a tidal-disruption event or a microquasar with strong disk winds. This cautions that photometric Pop~III selections are vulnerable to contamination when the rest-frame optical continuum is undetected. Motivated by this, we refine the photometric Pop III selection criteria to exclude the locus of extreme Balmer-jump objects. The revised criteria also recover the recently reported spectroscopic candidate AMORE6, demonstrating that the updated selection preserves sensitivity to genuine Pop III-like sources while removing key contaminants. Applying the refined criteria across legacy survey fields and five newly released CANUCS lensing cluster fields, we revisit the Pop III UV luminosity function and estimate the Pop III cosmic star-formation rate density to be $\approx[10^{-6}$--$10^{-4}]$~$M_{\odot}$~yr$^{-1}$~cMpc$^{-3}$ at $z\simeq6$--7, falling in the range of current theoretical predictions.
2512.11790
Dec 2025Astrophysics of Galaxies
The Effect of a Non-universal Extinction Curve on the Wesenheit Function and Cepheid Distances
The Wesenheit function is widely used to reduce the effects of interstellar reddening in distance measurements. Its construction, however, relies on the assumption of a universal extinction curve and on fixed values of the total-to-selective extinction ratio, Rv. Recent studies have shown that Rv varies significantly across the Milky Way and between different galaxies, raising concerns about systematic biases in Wesenheit magnitudes and period-Wesenheit relations. In this work, we discuss the impact of non-universal extinction on Wesenheit indices by combining the Rv-dependent extinction curve with a grid of stellar atmosphere models. We compute the integrated extinction in optical and near-infrared passbands, derive Rv-dependent R coefficients for multiple Wesenheit indices, and examine how changes in Rv propagate into Wesenheit magnitudes and Cepheid distances in our Galaxy. We find that the R coefficients in the Wesenheit functions vary strongly with Rv. For classical Cepheids in the Milky Way disk, variations of Rv within the typical observed range (2.6-3.6) can lead to substantial differences in the Wesenheit function, reaching +-0.7 mag from the mean for the Gaia-based Wesenheit index W_G and resulting in distance errors of almost 40%. Near-infrared Wesenheit indices are much less sensitive to Rv changes. Our results clearly show that accounting for variable Rv is essential when applying period-Wesenheit relations, particularly in the optical regime, or that near or mid infrared based distances should be used. While we present this effect for classical Cepheids, it applies to all pulsating stars for which period-Wesenheit relations are used to infer distances.
2512.11758
Dec 2025Astrophysics of Galaxies
Rapid sinking and efficient mergers of supermassive black holes in compact high-redshift galaxies
We present a cosmological zoom-in simulation targeting the high redshift compact progenitor phase of massive galaxies, with the most massive galaxy reaching a stellar mass of $M_{\star}=8.5\times 10^{10} \ M_{\odot}$ at $z=5$. The dynamics of supermassive black holes (SMBHs) is modelled from seeding down to their coalescence at sub-parsec scales due to gravitational wave (GW) emission by utilising a new version of the KETJU code, which combines regularised integration of sufficiently massive SMBHs with a dynamical friction subgrid model for lower-mass SMBHs. All nine massive galaxies included in this study go through a gas-dominated phase of early compaction in the redshift range of $z\sim 7-9$, starting at stellar masses of $M_\star\gtrsim 10^8\ \mathrm{M}_\odot$ and ending at a few times $M_{\star}\sim 10^9\ \mathrm{M}_\odot$. The sizes, masses and broad band fluxes of these compact systems are in general agreement with the population of systems observed with JWST known as `Little Red Dots'. In the compact phase, the stellar and SMBH masses grow rapidly, leading to a sharp decline in the central gas fractions. The outer regions, however, remain relatively gas-rich, leading to subsequent off-centre star formation and size growth. Due to the very high central stellar densities ($ρ_{\star}\gtrsim 10^{13}\,\mathrm{M_\odot/kpc^3}$), the SMBHs merge rapidly, typically just $\sim 4-35\ \mathrm{Myr}$ after the SMBH binaries have become bound. Combining KETJU with the phenomenological PhenomD model resolves the complete evolution of the GW emission from SMBH binaries through the Pulsar Timing Array frequency waveband up to the final few orbits that produce GWs observable with the future LISA mission.
2512.11665
Dec 2025Astrophysics of Galaxies
A Critical Evaluation of the Physical Nature of the Little Red Dots
Little Red Dots (LRDs) are a newly identified class of active galactic nuclei (AGNs) uncovered by JWST deep surveys. Their enigmatic properties challenge the canonical AGN paradigm and have stimulated ideas on early massive black hole (BH) formation. In this review, we summarize how early BHs shape the characteristic features of LRDs, how their nuclear environments differ from those of normal AGNs, and how future observations can distinguish between competing scenarios. Our main conclusions are as follows: (1) LRDs show broad-line emission consistent with mass accretion onto BHs with $M_{\rm BH}\simeq 10^{6-7}~M_\odot$, suggesting that AGN activity is a plausible origin of their dominant red optical emission. (2) Stellar components can reproduce the continuum energetics through dusty star formation. However, the required stellar mass would be too large to remain consistent with other LRD properties. Therefore, a purely stellar origin is unlikely to be the dominant power source, although star formation may still contribute to the UV emission. (3) The coexistence of broad-line emission with Balmer absorption and break features on LRD spectra, neither of which can be explained by stellar populations, suggests that nuclear BHs are enshrouded by dense gas with a high covering fraction. (4) Gas-enshrouded AGNs can produce red optical spectra without requiring dust reddening through a combination of gas attenuation and thermal self-emission with an effective temperature of $T_{\rm eff}\simeq 5000~{\rm K}$, which also accounts for the flat infrared continuum. (5) From the spectral features and redshift evolution, LRDs are likely a transient phase in early BH growth, possibly the first accretion episodes of newborn BHs. (6) Testing models for LRD spectra and origins through time variability, ionizing sources, post-LRD objects, and low-redshift analogs is particularly promising.
2512.031303
Dec 2025Astrophysics of Galaxies
Tracing Quenching in Nearby Galaxies Through Inner Surface Mass Density and Cold Gas Content
The inner stellar mass surface density within 1 kpc, Sigma1, has emerged as a suitable proxy for bulge growth and galaxy quenching. However, the dependence of cold gas content on Sigma1 has not been thoroughly explored. In this paper, we examine the relationship between Sigma1, as well as the mass-relative parameter delSigma1, and the atomic (fHI) and molecular (fH2) cold gas fractions in massive, nearby galaxies. We utilize a sample of 341 galaxies with HI data and 201 galaxies with H2 data from the xGASS and xCOLDGASS surveys, spanning 0.02 <= z <= 0.05 and a stellar mass range of 10^10 <= M_*/M_odot <= 10^11.5. While we observe that a decline in both fHI and fH2 is associated with increasing Sigma1, we find that fH2 shows a sharper decline above a threshold value of delSigma1 = 0. In addition, the fraction of galaxies with AGN activity (Seyferts and LINERs) increases with delSigma1, with the greatest increase occurring between 0 <= delSigma1 <= 0.2 dex. We propose an evolutionary track in the plane of fH2-delSigma1, whereby molecular gas depletion at fixed mass coincides with a rise in AGN activity. Our results suggest that central bulge growth is more tightly coupled to the depletion of molecular gas rather than atomic gas, with AGN feedback possibly contributing to this process. Our work highlights the utility of Sigma1 and delSigma1 as tracers of quenching in massive galaxies.
2511.18227
Nov 2025Astrophysics of Galaxies
Euclid Quick Data Release (Q1): Hunting for luminous z > 6 galaxies in the Euclid Deep Fields -- forecasts and first bright detections
The evolution of the rest-frame ultraviolet luminosity function (UV LF) is a powerful probe of early star formation and stellar mass build-up. At z > 6, its bright end (MUV < -21) remains poorly constrained due to the small volumes of existing near-infrared (NIR) space-based surveys. The Euclid Deep Fields (EDFs) will cover 53 deg^2 with NIR imaging down to 26.5 AB, increasing area by a factor of 100 over previous space-based surveys. They thus offer an unprecedented opportunity to select bright z > 6 Lyman break galaxies (LBGs) and constrain the UV LF's bright end. With NIR coverage extending to 2um, Euclid can detect galaxies out to z = 13. We present forecasts for the number densities of z > 6 galaxies expected in the final EDF dataset. Using synthetic photometry from spectral energy distribution (SED) templates of z = 5--15 galaxies, z = 1--4 interlopers, and Milky Way MLT dwarfs, we explore optimal selection methods for high-z LBGs. A combination of S/N cuts with SED fitting (from optical to MIR) yields the highest-fidelity sample, recovering >76% of input z > 6 LBGs while keeping low-z contamination <10%. This excludes instrumental artefacts, which will affect early Euclid releases. Auxiliary data are critical: optical imaging from the Hyper Suprime-Cam and Vera C. Rubin Observatory distinguishes genuine Lyman breaks, while Spitzer/IRAC data help recover z > 10 sources. Based on empirical double power-law LF models, we expect >100,000 LBGs at z = 6-12 and >100 at z > 12 in the final Euclid release. In contrast, steeper Schechter models predict no z > 12 detections. We also present two ultra-luminous (MUV < -23.5) candidates from the EDF-N Q1 dataset. If their redshifts are confirmed, their magnitudes support a DPL LF model at z > 9, highlighting Euclid's power to constrain the UV LF's bright end and identify the most luminous early galaxies for follow-up.
2511.02926
Nov 2025Astrophysics of Galaxies
