Trending in Astrophysics
Dissecting the Hubble tension: Insights from a diverse set of Sound Horizon-free H0 measurements
The Hubble tension is commonly framed as a discrepancy between local, late-time measurements favoring $H_0 \approx 73$ km s$^{-1}$ Mpc$^{-1}$ and early-time, Sound-Horizon-based measurements favoring $H_0 \approx 67$ km s$^{-1}$ Mpc$^{-1}$. We challenge this viewpoint by analyzing 83 Sound-Horizon-independent $H_0$ measurements, categorized into four classes: Distance Ladder measurements using local calibrators; Local One-Step $Λ$CDM measurements assuming the standard expansion history; Pure Local One-Step measurements independent of $H(z)$ shape; and CMB Sound Horizon free measurements using CMB data without the Sound Horizon scale. We find that the 29 Distance Ladder measurements yield $H_0 = 72.74 \pm 0.40$ km s$^{-1}$ Mpc$^{-1}$ ($χ^2_ν\equiv χ^2/d.o.f= 0.74$), while the 54 One-Step measurements collectively yield $H_0 = 68.67 \pm 0.46$ km s$^{-1}$ Mpc$^{-1}$ ($χ^2_ν= 0.85$), a $6.7σ$ tension exceeding the Planck--SH0ES discrepancy. This tension remains significant at $4.5σ$ after accounting for correlations. Among One-Step categories, Local One-Step $Λ$CDM measurements favor the lowest value ($H_0 = 67.18 \pm 0.90$ km s$^{-1}$ Mpc$^{-1}$), Pure Local One-Step yield an intermediate value ($H_0 = 70.38 \pm 1.00$ km s$^{-1}$ Mpc$^{-1}$), and CMB Sound Horizon Free measurements give $H_0 = 68.71 \pm 0.63$ km s$^{-1}$ Mpc$^{-1}$. Thus, that the Hubble tension is better characterized as a discrepancy between the Distance Ladder and all other methodologies, rather than an early-vs-late-time split. We also identify a $2.4σ$ internal tension among One-Step measurements: analyses assuming $Λ$CDM systematically recover lower $H_0$ values by about 3.2 km s$^{-1}$ Mpc$^{-1}$ compared to model-independent methods. This suggests either unrecognized systematics/physics in the Distance Ladder or deviations from $Λ$CDM in the late-time Universe.
2601.00650
Jan 2026Cosmology and Nongalactic Astrophysics
Dynamical Dark Energy Imprints in the Lyman-Alpha Forest
The nature of dark energy (DE) remains elusive, even though it constitutes the dominant energy-density component of the Universe and drives the late-time acceleration of cosmic expansion. By combining measurements of the expansion history from baryon acoustic oscillations, supernova surveys, and cosmic microwave background data, the Dark Energy Spectroscopic Instrument (DESI) Collaboration has inferred that the DE equation of state may evolve over time. The profound implications of a time-variable, ``dynamical" DE (DDE) that departs from a cosmological constant motivate the need for independent observational tests. In this work, we use cosmological hydrodynamical simulations of structure formation to investigate how DDE affects the properties of the Lyman-Alpha ``forest'' of absorption features produced by neutral hydrogen in the cosmic web. We find that DDE models consistent with the DESI constraints induce a spectral tilt in the forest transmitted flux power spectrum, imprinting a scale- and redshift-dependent signature relative to standard Lambda-CDM cosmologies. These models also yield higher intergalactic medium temperatures and reduced Lyman-Alpha opacity compared to Lambda-CDM. We discuss the observational implications of these trends as potential avenues for independent confirmation of DDE.
2601.00767
Jan 2026Cosmology and Nongalactic Astrophysics
Irradiated Atmosphere V: Effects of Vertical-Mixing induced Energy Transport on the Inhomogeneity
Atmospheric variations over time and space boost planetary cooling, as outgoing internal flux responds to stellar radiation and opacity. Vertical mixing regulates this cooling. Our study examines how gravity waves or large-scale induced mixing interact with radiation transfer, affecting temperature inhomogeneity and internal flux. Through the radiative-convective-mixing equilibrium, mixing increases temperature inhomogeneity in the middle and lower atmospheres, redistributing internal flux. Stronger stellar radiation and mixing significantly reduce outgoing flux, slowing cooling. With constant infrared (IR) opacity, lower visible opacity and stronger mixing significantly reduce outgoing flux. Jensen's inequality implies that greater spatial disparities in stellar flux and opacity elevate the ratio of the average internal flux in inhomogeneous columns relative to that in homogeneous columns. This effect, particularly pronounced under high opacity contrasts, amplifies deep-layer temperature inhomogeneity and may enhance cooling. However, with mixing, overall cooling is weaker than without, as both the averaged internal flux of the inhomogeneous columns and that of the homogeneous column decline more sharply for the latter. Thus, while vertical mixing-induced inhomogeneity can enhance cooling, the overall cooling effect remains weaker than in the non-mixing case. Therefore, vertical mixing, by regulating atmospheric structure and flux, is key to understanding planetary cooling.
2601.00606
Jan 2026Earth and Planetary Astrophysics
Callisto's Nonresonant Orbit as an Outcome of Circum-Jovian Disk Substructure
The Galilean moons of Io, Europa, and Ganymede exhibit a 4:2:1 commensurability in their mean motions, a configuration known as the Laplace resonance. The prevailing view for the origin of this three-body resonance involves the convergent migration of the moons, resulting from gas-driven torques in the circum-Jovian disk wherein they accreted. To account for Callisto's exclusion from the resonant chain, a late and/or slow accretion of the fourth and outermost Galilean moon is typically invoked, stalling its migration. Here, we consider an alternative scenario in which Callisto's nonresonant orbit is a consequence of disk substructure. Using a suite of N-body simulations that self-consistently account for satellite-disk interactions, we show that a pressure bump can function as a migration trap, isolating Callisto and alleviating constraints on its timing of accretion. Our simulations position the bump interior to the birthplaces of all four moons. In exploring the impact of bump structure on simulation outcomes, we find that it cannot be too sharp nor flat to yield the observed orbital architecture. In particular, a "Goldilocks" zone is mapped in parameter space, corresponding to a well-defined range in bump aspect ratio. Within this range, Io, Europa, and Ganymede are sequentially trapped at the bump, and ushered across it through resonant lockstep migration with their neighboring, exterior moon. The implications of our work are discussed in the context of uncertainties regarding Callisto's interior structure, arising from the possibility of non-hydrostatic contributions to its shape and gravity field, unresolved by the Galileo spacecraft.
2601.00786
Jan 2026Earth and Planetary Astrophysics
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 free-floating-planet microlensing event caused by a Saturn-mass object
A population of free-floating planets is known from gravitational microlensing surveys. None have a directly measured mass, owing to a degeneracy with the distance, but the population statistics indicate that many are less massive than Jupiter. We report a microlensing event -- KMT-2024-BLG-0792/OGLE-2024-BLG-0516, which was observed from both ground- and space-based telescopes -- that breaks the mass-distance degeneracy. The event was caused by an object with 0.219^{+0.075}_{-0.046} Jupiter masses that is either gravitationally unbound or on a very wide orbit. Through comparison with the statistical properties of other observed microlensing events and predictions from simulations, we infer that this object likely formed in a protoplanetary disk (like a planet), not in isolation (like a brown dwarf), and dynamical processes then ejected it from its birth place, producing a free-floating object.
2601.000571
Dec 2025Earth and Planetary Astrophysics
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
Betelgeuse: Detection of the Expanding Wake of the Companion Star
Recent analyses conclude that Betelgeuse, a red supergiant star (HD 39801), likely has a companion object with a period of about 2000 days orbiting at only 2.3 stellar radii, deep in the chromosphere of the supergiant. A probable detection of such a companion, named Siwarha, has just occurred from speckle imaging. This study finds that Betelgeuse spectra in the optical region and ultraviolet exhibit signatures of variable circumstellar absorption and chromospheric outflows. These variations are consistent with the ~ 2000-day period of the companion object. Circumstellar absorption evident in optical Mn I lines, and mass outflow marked by ultraviolet Fe II, Si I, and Mg I lines increase after the transit of the companion across the disk of Betelgeuse. Following the eclipse of the companion, the absorption and outflow slowly decrease in advance of the next transit. The occurrence and variation of this plasma appear consistent with the presence of a trailing and expanding wake caused by a companion star orbiting within the atmosphere of Betelgeuse.
2601.00470
Jan 2026Solar and Stellar Astrophysics
Hadronic Origin of Sub-PeV Gamma-Ray Emission from LHAASO J0621+3755
Very High Energy (VHE) gamma rays are primarily estimated to arise from high-energy electromagnetic interactions in pulsars and their halo through electron inverse Compton (IC) scattering. Hadronic channels like neutral pion decay have also been proposed to produce TeV-PeV gamma rays from the Pulsar halo. The neutral pions are expected to be generated from cosmic ray (CR) protons interacting with the ambient/cloud. The recent observations of sub-PeV gamma rays from the halo of pulsar PSR J0622+3749 by the Large High Altitude Air Shower Observatory Kilometre-Square Array (LHAASO-KM2A) detector provide a platform to explore different channels of their production. Previous studies support consistency with the leptonic modeling of the halo, which attributes its origin to slow diffusion in the interstellar medium. In this work, we investigated the possibility of proton-proton channel as the origin of these photons. To explain the observed gamma rays with energy $\sim 4$ TeV by the High-Altitude Water Cherenkov (HAWC) telescope till 200 TeV by the LHAASO observatory, one requires the CR proton luminosity to be $η_p\sim 0.1$ of the pulsar PSR J0622+3749 spin-down luminosity. In this case, we have considered the protons propagating in a one-zone superdiffusion environment, specifically $α= 1$ in a cloud of gas density 1 per cm$^{3}$.
2601.00692
Jan 2026High Energy Astrophysical Phenomena
Readdressing the contribution of photonuclear reactions to the muon content of extensive air showers: a heuristic approach
The indirect ground-based observations of cosmic rays through extensive air showers in modern experiments typically involve the use of Monte Carlo simulations to determine the characteristics of the primary particles. These simulations necessitate assumptions about particle interactions at energies that have not yet been experimentally probed, which introduces systematic uncertainties in key observables, particularly the number of muons. Current research on this uncertainty primarily focuses on hadronic interaction models, the dominant source of muon production. This study presents an approach that takes into account another significant mechanism for muon generation: photonuclear reactions. A robust heuristic technique has been developed to estimate the contribution of these interactions to the total number of muons over a wide range of extensive air shower parameters (including primary particle type, energy, and slant atmospheric depth) and photonuclear interaction models, with an absolute percentage error on the order of $10\%$ in the estimated number of muons. Furthermore, several potential applications of the suggested method in relation to modern challenges in extensive air shower physics are discussed.
2512.12481
Dec 2025High Energy Astrophysical Phenomena
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
One-loop power spectrum corrections in interacting dark energy cosmologies
Interacting Dark Energy (IDE) models offer a promising avenue to explore possible exchanges of energy and momentum between dark matter and dark energy, providing a dynamical extension of the standard $Λ$CDM paradigm. Such interactions modify the growth of cosmic structures, imprinting distinctive signatures on the matter power spectrum that can be tested through large-scale structure (LSS) observations. In this work, we compute the one-loop corrections to the matter power spectrum in IDE models. We then reinterpret these results within the standard framework of the Effective Field Theory of Large-Scale Structure (EFTofLSS), which provides a consistent description of mildly non-linear scales and allows for reliable comparisons with observational data. We investigate two commonly studied forms of the coupling function, $Q$, namely $Q = ξ\mathcal{H} ρ_{\rm m}$ and $Q = ξ\mathcal{H} ρ_{\rm DE}$, and introduce a novel interaction term, $Q = Γ\, ρ_{\rm m} \, ρ_{\rm DE} \, θ_{\rm m}$, characterized by the non-linear coupling constant $Γ$, which links the interaction strength to the velocity divergence of dark matter. This coupling function is proposed to isolate the effects solely of the IDE model on mildly non-linear scales. Using Full-Shape (FS) measurements of the galaxy power spectrum from BOSS DR12, we constrain the interaction rate $Γ$, the cosmological parameters, and the bias parameters. We find $Γ= 0.0039 \pm 0.0082$, which is highly consistent with the $Λ$CDM model. This work opens the possibility of testing IDE models at mildly non-linear scales, potentially providing new insights for this class of models beyond the standard $Λ$CDM framework.
2512.11678
Dec 2025Cosmology and Nongalactic Astrophysics
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.031301
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