We present the La Silla Schmidt Southern Survey (LS4), a new wide-field, time-domain survey to be conducted with the 1 m ESO Schmidt telescope. The 268 megapixel LS4 camera mosaics 32 2k$\times$4k fully depleted CCDs, providing a $\sim$20 deg$^2$ field of view with $1''$ pixel$^{-1}$ resolution. The LS4 camera will have excellent performance at longer wavelengths: in a standard 45 s exposure the expected 5$\sigma$ limiting magnitudes in $g$, $i$, $z$ are $\sim$21.5, $\sim$20.9, and $\sim$20.3 mag (AB), respectively. The telescope design requires a novel filter holder that fixes different bandpasses over each quadrant of the detector. Two quadrants will have $i$ band, while the other two will be $g$ and $z$ band and color information will be obtained by dithering targets across the different quadrants. The majority (90%) of the observing time will be used to conduct a public survey that monitors the extragalactic sky at both moderate (3 d) and high (1 d) cadence, as well as focused observations within the Galactic bulge and plane. Alerts from the public survey will be broadcast to the community via established alert brokers. LS4 will run concurrently with the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST). The combination of LS4+LSST will enable detailed holistic monitoring of many nearby transients: high-cadence LS4 observations will resolve the initial rise and peak of the light curve while less-frequent but deeper observations by LSST will characterize the years before and after explosion. Here, we summarize the primary science objectives of LS4 including microlensing events in the Galaxy, extragalactic transients, the search for electromagnetic counterparts to multi-messenger events, and cosmology.
We demonstrate that primordial black holes (PBHs) lighter than $10^9 \, \text{g}$, which evaporated before the Big Bang nucleosynthesis, can induce significant isocurvature perturbations due to their biased clustering amplitude and the branching ratio of the Hawking radiation differing from the abundance ratio. By leveraging the upper bound on the isocurvature perturbations from the cosmic microwave background anisotropies reported by the Planck collaboration, we derive a new upper bound on the abundance of these light PBHs.
We present baryon acoustic oscillation (BAO) measurements from more than 14 million galaxies and quasars drawn from the Dark Energy Spectroscopic Instrument (DESI) Data Release 2 (DR2), based on three years of operation. For cosmology inference, these galaxy measurements are combined with DESI Lyman-$\alpha$ forest BAO results presented in a companion paper. The DR2 BAO results are consistent with DESI DR1 and SDSS, and their distance-redshift relationship matches those from recent compilations of supernovae (SNe) over the same redshift range. The results are well described by a flat $\Lambda$CDM model, but the parameters preferred by BAO are in mild, $2.3\sigma$ tension with those determined from the cosmic microwave background (CMB), although the DESI results are consistent with the acoustic angular scale $\theta_*$ that is well-measured by Planck. This tension is alleviated by dark energy with a time-evolving equation of state parametrized by $w_0$ and $w_a$, which provides a better fit to the data, with a favored solution in the quadrant with $w_0>-1$ and $w_a<0$. This solution is preferred over $\Lambda$CDM at $3.1\sigma$ for the combination of DESI BAO and CMB data. When also including SNe, the preference for a dynamical dark energy model over $\Lambda$CDM ranges from $2.8-4.2\sigma$ depending on which SNe sample is used. We present evidence from other data combinations which also favor the same behavior at high significance. From the combination of DESI and CMB we derive 95% upper limits on the sum of neutrino masses, finding $\sum m_\nu<0.064$ eV assuming $\Lambda$CDM and $\sum m_\nu<0.16$ eV in the $w_0w_a$ model. Unless there is an unknown systematic error associated with one or more datasets, it is clear that $\Lambda$CDM is being challenged by the combination of DESI BAO with other measurements and that dynamical dark energy offers a possible solution.
We present the Baryon Acoustic Oscillation (BAO) measurements with the Lyman-alpha (LyA) forest from the second data release (DR2) of the Dark Energy Spectroscopic Instrument (DESI) survey. Our BAO measurements include both the auto-correlation of the LyA forest absorption observed in the spectra of high-redshift quasars and the cross-correlation of the absorption with the quasar positions. The total sample size is approximately a factor of two larger than the DR1 dataset, with forest measurements in over 820,000 quasar spectra and the positions of over 1.2 million quasars. We describe several significant improvements to our analysis in this paper, and two supporting papers describe improvements to the synthetic datasets that we use for validation and how we identify damped LyA absorbers. Our main result is that we have measured the BAO scale with a statistical precision of 1.1% along and 1.3% transverse to the line of sight, for a combined precision of 0.65% on the isotropic BAO scale at $z_{eff} = 2.33$. This excellent precision, combined with recent theoretical studies of the BAO shift due to nonlinear growth, motivated us to include a systematic error term in LyA BAO analysis for the first time. We measure the ratios $D_H(z_{eff})/r_d = 8.632 \pm 0.098 \pm 0.026$ and $D_M(z_{eff})/r_d = 38.99 \pm 0.52 \pm 0.12$, where $D_H = c/H(z)$ is the Hubble distance, $D_M$ is the transverse comoving distance, $r_d$ is the sound horizon at the drag epoch, and we quote both the statistical and the theoretical systematic uncertainty. The companion paper presents the BAO measurements at lower redshifts from the same dataset and the cosmological interpretation.
this https URL 21 pages, 13 figures, 3 tables
We present the Damped Ly$\alpha$ Toolkit for automated detection and characterization of Damped Ly$\alpha$ absorbers (DLA) in quasar spectra. Our method uses quasar spectral templates with and without absorption from intervening DLAs to reconstruct observed quasar forest regions. The best-fitting model determines whether a DLA is present while estimating the redshift and HI column density. With an optimized quality cut on detection significance ($\Delta \chi_{r}^2>0.03$), the technique achieves an estimated 72% purity and 71% completeness when evaluated on simulated spectra with S/N$>2$ that are free of broad absorption lines (BAL). We provide a catalog containing candidate DLAs from the DLA Toolkit detected in DESI DR1 quasar spectra, of which 21,719 were found in S/N$>2$ spectra with predicted $\log_{10} (N_\texttt{HI}) > 20.3$ and detection significance $\Delta \chi_{r}^2 >0.03$. We compare the Damped Ly$\alpha$ Toolkit to two alternative DLA finders based on a convolutional neural network (CNN) and Gaussian process (GP) models. We present a strategy for combining these three techniques to produce a high-fidelity DLA catalog from DESI DR2 for the Ly$\alpha$ forest baryon acoustic oscillation measurement. The combined catalog contains 41,152 candidate DLAs with $\log_{10} (N_\texttt{HI}) > 20.3$ from quasar spectra with S/N$>2$. We estimate this sample to be approximately 76% pure and 71% complete when BAL quasars are excluded.
The Dark Energy Spectroscopic Instrument (DESI) data release 2 (DR2) galaxy and quasar clustering data represents a significant expansion of data from DR1, providing improved statistical precision in BAO constraints across multiple tracers, including bright galaxies (BGS), luminous red galaxies (LRGs), emission line galaxies (ELGs), and quasars (QSOs). In this paper, we validate the BAO analysis of DR2. We present the results of robustness tests on the blinded DR2 data and, after unblinding, consistency checks on the unblinded DR2 data. All results are compared to those obtained from a suite of mock catalogs that replicate the selection and clustering properties of the DR2 sample. We confirm the consistency of DR2 BAO measurements with DR1 while achieving a reduction in statistical uncertainties due to the increased survey volume and completeness. We assess the impact of analysis choices, including different data vectors (correlation function vs. power spectrum), modeling approaches and systematics treatments, and an assumption of the Gaussian likelihood, finding that our BAO constraints are stable across these variations and assumptions with a few minor refinements to the baseline setup of the DR1 BAO analysis. We summarize a series of pre-unblinding tests that confirmed the readiness of our analysis pipeline, the final systematic errors, and the DR2 BAO analysis baseline. The successful completion of these tests led to the unblinding of the DR2 BAO measurements, ultimately leading to the DESI DR2 cosmological analysis, with their implications for the expansion history of the Universe and the nature of dark energy presented in the DESI key paper.
We conduct an extended analysis of dark energy constraints, in support of the findings of the DESI DR2 cosmology key paper, including DESI data, Planck CMB observations, and three different supernova compilations. Using a broad range of parametric and non-parametric methods, we explore the dark energy phenomenology and find consistent trends across all approaches, in good agreement with the $w_0w_a$CDM key paper results. Even with the additional flexibility introduced by non-parametric approaches, such as binning and Gaussian Processes, we find that extending $\Lambda$CDM to include a two-parameter $w(z)$ is sufficient to capture the trends present in the data. Finally, we examine three dark energy classes with distinct dynamics, including quintessence scenarios satisfying $w \geq -1$, to explore what underlying physics can explain such deviations. The current data indicate a clear preference for models that feature a phantom crossing; although alternatives lacking this feature are disfavored, they cannot yet be ruled out. Our analysis confirms that the evidence for dynamical dark energy, particularly at low redshift ($z \lesssim 0.3$), is robust and stable under different modeling choices.
We present the results of the single-component Sérsic profile fitting for the magnitude-limited sample of \IE$<23$ galaxies within the 63.1 deg$^2$ area of the Euclid Quick Data Release (Q1). The associated morphological catalogue includes two sets of structural parameters fitted using \texttt{SourceXtractor++}: one for VIS \IE images and one for a combination of three NISP images in \YE, \JE and \HE bands. We compare the resulting Sérsic parameters to other morphological measurements provided in the Q1 data release, and to the equivalent parameters based on higher-resolution \HST imaging. These comparisons confirm the consistency and the reliability of the fits to Q1 data. Our analysis of colour gradients shows that NISP profiles have systematically smaller effective radii ($R_{\rm e}$) and larger Sérsic indices ($n$) than in VIS. In addition, we highlight trends in NISP-to-VIS parameter ratios with both magnitude and $n_{\rm VIS}$. From the 2D bimodality of the $(u-r)$ colour-$\log(n)$ plane, we define a $(u-r)_{\rm lim}(n)$ that separates early- and late-type galaxies (ETGs and LTGs). We use the two subpopulations to examine the variations of $n$ across well-known scaling relations at $z<1$. ETGs display a steeper size--stellar mass relation than LTGs, indicating a difference in the main drivers of their mass assembly. Similarly, LTGs and ETGs occupy different parts of the stellar mass--star-formation rate plane, with ETGs at higher masses than LTGs, and further down below the Main Sequence of star-forming galaxies. This clear separation highlights the link known between the shutdown of star formation and morphological transformations in the Euclid imaging data set. In conclusion, our analysis demonstrates both the robustness of the Sérsic fits available in the Q1 morphological catalogue and the wealth of information they provide for studies of galaxy evolution with Euclid.
Stellar bars are key structures in disc galaxies, driving angular momentum redistribution and influencing processes such as bulge growth and star formation. Quantifying the bar fraction as a function of redshift and stellar mass is therefore important for constraining the physical processes that drive disc formation and evolution across the history of the Universe. Leveraging the unprecedented resolution and survey area of the Euclid Q1 data release combined with the Zoobot deep-learning model trained on citizen-science labels, we identify 7711 barred galaxies with $M_* \gtrsim 10^{10}M_\odot$ in a magnitude-selected sample $I_E < 20.5$ spanning $63.1 deg^2$. We measure a mean bar fraction of $0.2-0.4$, consistent with prior studies. At fixed redshift, massive galaxies exhibit higher bar fractions, while lower-mass systems show a steeper decline with redshift, suggesting earlier disc assembly in massive galaxies. Comparisons with cosmological simulations (e.g., TNG50, Auriga) reveal a broadly consistent bar fraction, but highlight overpredictions for high-mass systems, pointing to potential over-efficiency in central stellar mass build-up in simulations. These findings demonstrate Euclid's transformative potential for galaxy morphology studies and underscore the importance of refining theoretical models to better reproduce observed trends. Future work will explore finer mass bins, environmental correlations, and additional morphological indicators.
Modern astronomical surveys, such as the Euclid mission, produce high-dimensional, multi-modal data sets that include imaging and spectroscopic information for millions of galaxies. These data serve as an ideal benchmark for large, pre-trained multi-modal models, which can leverage vast amounts of unlabelled data. In this work, we present the first exploration of Euclid data with AstroPT, an autoregressive multi-modal foundation model trained on approximately 300 000 optical and infrared Euclid images and spectral energy distributions (SEDs) from the first Euclid Quick Data Release. We compare self-supervised pre-training with baseline fully supervised training across several tasks: galaxy morphology classification; redshift estimation; similarity searches; and outlier detection. Our results show that: (a) AstroPT embeddings are highly informative, correlating with morphology and effectively isolating outliers; (b) including infrared data helps to isolate stars, but degrades the identification of edge-on galaxies, which are better captured by optical images; (c) simple fine-tuning of these embeddings for photometric redshift and stellar mass estimation outperforms a fully supervised approach, even when using only 1% of the training labels; and (d) incorporating SED data into AstroPT via a straightforward multi-modal token-chaining method improves photo-z predictions, and allow us to identify potentially more interesting anomalies (such as ringed or interacting galaxies) compared to a model pre-trained solely on imaging data.
The extent to which the environment affects galaxy evolution has been under scrutiny by researchers for decades. With the first data from Euclid, we can begin to study a wide range of environments and their effects as a function of redshift, using 63 sq deg of space-based data. In this paper, we present results from the Euclid Q1 Release, where we measure the passive-density and morphology-density relations at z=0.25--1. We determine if a galaxy is passive using the specific star-formation rate, and we classify the morphologies of galaxies using the Sérsic index n and the u-r colours. We measure the local environmental density of each galaxy using the Nth-nearest neighbour method. We find that at fixed stellar mass, the quenched fraction increases with increasing density up to z=0.75. This result shows the separability of the effects from the stellar mass and the environment, at least at z<0.75. At z>0.75, we observe weak environmental effects, with most high mass galaxies being quenched independently of environment. Up to z=0.75, the ETG fraction increases with density at fixed stellar mass, meaning the environment also transforms the morphology of the galaxy independently of stellar mass, at low mass. For high mass galaxies, almost all galaxies are early-types, with low impact from the environment. At z>0.75, the morphology depends mostly on stellar mass, with only low-mass galaxies being affected by the environment. Given that the morphology classifications use u-r colours, these are correlated to the star-formation rate, and as such our morphology results should be taken with caution; future morphology classifications should verify these results. To summarise, we identify the passive-density and morphology-density relations at z<0.75, but at z>0.75 the relations are less strong. At z>0.75, the uncertainties are large, and future Euclid data releases are key to confirm these trends.
The star-forming main sequence (SFMS) is a tight relation observed between stellar masses and star-formation rates (SFR) in a population of galaxies. The relation holds for different redshift, morphological, and environmental domains, and is a key to understanding the underlying relations between a galaxy budget of cold gas and its stellar content. \Euclid Quick Data Release 1 (Q1) gives the opportunity to investigate this fundamental relation in galaxy formation and evolution. We complement the \Euclid release with public IRAC observations of the \Euclid Deep Fields (EDFs), improving the quality of recovered photometric redshifts, stellar masses, and star-formation rates, as shown both from simulations and comparison with available spectroscopic redshifts. From Q1 data alone, we recover more than $\sim 30\,\mathrm{k}$ galaxies with $\log_{10}(M_\ast/M_\odot) > 11$, giving a precise constraint of the SFMS at the high-mass end. We investigated SFMS, in a redshift interval between $0.2$ and $3.0$, comparing our results with the existing literature and fitting them with a parameterisation taking into account the presence of a bending of the relation at the high-mass end, depending on the bending mass $M_0$. We find good agreement with previous results in terms of $M_0$ values. We also investigate the distribution of physical (e.g., dust absorption $A_V$ and formation age) and morphological properties (e.g., Sérsic index and radius) in the SFR--stellar mass plane, and their relation with the SFMS. These results highlight the potential of \Euclid in studying the fundamental scaling relations that regulate galaxy formation and evolution in anticipation of the forthcoming Data Release 1.
Investigating the drivers of the quenching of star formation in galaxies is key to understanding their evolution. The Euclid mission will provide rich spatial and spectral data from optical to infrared wavelengths for millions of galaxies, enabling precise measurements of their star formation histories. Using the first Euclid Quick Data Release (Q1), we developed a probabilistic classification framework, that combines the average specific star-formation rate ($\rm sSFR_\tau$) inferred over two timescales ($\tau={10^8,10^9}$ yr), to categorize galaxies as `Ageing' (secularly evolving), `Quenched' (recently halted star formation), or `Retired' (dominated by old stars). We validated this methodology using synthetic observations from the IllustrisTNG simulation. Two classification methods were employed: a probabilistic approach, integrating posterior distributions, and a model-driven method optimizing sample purity and completeness using IllustrisTNG. At $z<0.1$ and $M_\ast \gtrsim 3\times10^{8}\, M_\odot$, we obtain Euclid class fractions of 68-72%, 8-17%, and 14-19% for Ageing, Quenched, and Retired populations, respectively, consistent with previous studies. The evolution with redshift shows increasing/decreasing fraction of Ageing/Retired galaxies. The fraction of quenched systems shows a weaker dependence on stellar mass and redshift, varying between 5% and 15%. We analysed the mass-size-metallicity relation for each population. Ageing galaxies generally exhibit disc morphologies and low metallicities. Retired galaxies show compact structures and enhanced chemical enrichment, while Quenched galaxies form an intermediate population, more compact and chemically evolved than Ageing systems. This work demonstrates Euclid's great potential for elucidating the physical nature of the quenching mechanisms that govern galaxy evolution.
To better understand the role of active galactic nuclei (AGN) in galaxy evolution, it is crucial to achieve a complete and pure AGN census. X-ray surveys are key to this, but identifying their counterparts (CTPs) at other wavelengths remains challenging due to their larger positional uncertainties and limited availability of deeper, uniform ancillary data. Euclid is revolutionising this effort, offering extensive coverage of nearly the entire extragalactic sky, particularly in the near-infrared bands, where AGN are more easily detected. With the first Euclid Quick Data Release (Q1), we identifyed, classifyed, and determined the redshifts of Euclid CTPs to known point-like sources from major X-ray surveys, including XMM-Newton, Chandra, and eROSITA. Using Bayesian statistics, combined with machine learning (ML), we identify the CTPs to 11 286 X-ray sources from the three X-ray telescopes. For the large majority of 10 194 sources, the associations are unique, with the remaining $\sim$ 10% of multi-CTP cases equally split between XMM-Newton and eROSITA. ML is then used to distinguish between Galactic (8%) and extragalactic (92%) sources. We computed photo-zs using deep learning for the 8617 sources detected in the 10th data release of the DESI Legacy Survey, reaching an accuracy and a fraction of outliers of about 5%. Based on their X-ray luminosities, over 99% of CTPs identified as extragalactic are classified as AGN, most of which appear unobscured given their hardness ratios. With this paper, we release our catalogue, which includes identifiers, basic X-ray properties, the details of the associations, and additional features such as Galactic/extragalactic classifications and photometric/spectroscopic redshifts. We also provide probabilities for sub-selecting the sample based on purity and completeness, allowing users to tailor the sample according to their specific needs.
Galaxy major mergers are a key pathway to trigger AGN. We present the first detection of major mergers in the Euclid Deep Fields and analyse their connection with AGN. We constructed a stellar-mass-complete ($M_*>10^{9.8}\,M_{\odot}$) sample of galaxies from the first quick data release (Q1), in the redshift range z=0.5-2. We selected AGN using X-ray data, optical spectroscopy, mid-infrared colours, and processing \IE observations with an image decomposition algorithm. We used CNNs trained on cosmological simulations to classify galaxies as mergers and non-mergers. We found a larger fraction of AGN in mergers compared to the non-merger controls for all AGN selections, with AGN excess factors ranging from 2 to 6. Likewise, a generally larger merger fraction ($f_{merg}$) is seen in active galaxies than in the non-active controls. We analysed $f_{merg}$ as a function of the AGN bolometric luminosity ($L_{bol}$) and the contribution of the point-source to the total galaxy light in the \IE-band ($f_{PSF}$) as a proxy for the relative AGN contribution fraction. We uncovered a rising $f_{merg}$, with increasing $f_{PSF}$ up to $f_{PSF}=0.55$, after which we observed a decreasing trend. We then derived the point-source luminosity ($L_{PSF}$) and showed that $f_{merg}$ monotonically increases as a function of $L_{PSF}$ at z<0.9, with $f_{merg}>$50% for $L_{PSF}>2\,10^{43}$ erg/s. At z>0.9, $f_{merg}$ rises as a function of $L_{PSF}$, though mergers do not dominate until $L_{PSF}=10^{45}$ erg/s. For X-ray and spectroscopic AGN, we computed $L_{bol}$, which has a positive correlation with $f_{merg}$ for X-ray AGN, while shows a less pronounced trend for spectroscopic AGN due to the smaller sample size. At $L_{bol}>10^{45}$ erg/s, AGN mostly reside in mergers. We concluded that mergers are strongly linked to the most powerful, dust-obscured AGN, associated with rapid supermassive black hole growth.
Active galactic nuclei (AGN) play a key role in galaxy evolution but are challenging to identify due to their varied observational signatures. Furthermore, understanding their impact requires quantifying their strength relative to their host galaxies. We developed a deep learning (DL) model for identifying AGN in imaging data by deriving the contribution of the central point source. Trained on Euclidised mock galaxy images with injected AGN levels, in the form of varying contributions of the point-spread function (PSF), our model can precisely and accurately recover the injected AGN contribution fraction $f_{\rm PSF}$, with a mean difference between the predicted and true $f_{\rm PSF}$ of $-0.0078$ and an overall root mean square error (RMSE) of 0.051. This method moves beyond binary AGN classification, enabling precise AGN contribution measurements. Applying our model to a stellar-mass-limited sample ($M_{\ast} \ge 10^{9.8} M_{\odot}$, $0.5 \le z \le 2.0$) from the first \Euclid quick data release (Q1), we identify $48,840 \pm 78$ AGN over 63.1 deg$^2$ ($7.8\pm0.1$%) using a threshold of $f_{\rm PSF} > 0.2$. We compare our DL-selected AGN with those identified in X-ray, mid-infrared (MIR), and optical spectroscopy and investigate their overlapping fractions depending on different thresholds on the PSF contribution. We find that the overlap increases with increasing X-ray or bolometric AGN luminosity. The AGN luminosity in the $I_{\rm E}$ filter correlates with host galaxy stellar mass, suggesting faster supermassive black hole (SMBH) growth in more massive galaxies. Moreover, the mean relative contribution of the AGN is higher in quiescent galaxies than in star-forming ones. Starburst galaxies and the most massive galaxies (across the star-formation main sequence) tend to host the most luminous AGN, indicating concomitant assembly of the SMBH and the host galaxy.
Red quasars constitute an important but elusive phase in the evolution of supermassive black holes, where dust obscuration can significantly alter their observed properties. They have broad emission lines, like other quasars, but their optical continuum emission is significantly reddened, which is why they were traditionally identified based on near- and mid-infrared selection criteria. This work showcases the capability of the \Euclid space telescope to find a large sample of red quasars, using \Euclid near infrared (NIR) photometry. We first conduct a forecast analysis, comparing a synthetic catalogue of red QSOs with COSMOS2020. Using template fitting, we reconstruct \Euclid-like photometry for the COSMOS sources and identify a sample of candidates in a multidimensional colour-colour space achieving $98\%$ completeness for mock red QSOs with $30\%$ contaminants. To refine our selection function, we implement a probabilistic Random Forest classifier, and use UMAP visualisation to disentangle non-linear features in colour-space, reaching $98\%$ completeness and $88\%$ purity. A preliminary analysis of the candidates in the \Euclid Deep Field Fornax (EDF-F) shows that, compared to VISTA+DECAm-based colour selection criteria, \Euclid's superior depth, resolution and optical-to-NIR coverage improves the identification of the reddest, most obscured sources. Notably, the \Euclid exquisite resolution in the $I_E$ filter unveils the presence of a candidate dual quasar system, highlighting the potential for this mission to contribute to future studies on the population of dual AGN. The resulting catalogue of candidates, including more the 150 000 sources, provides a first census of red quasars in \Euclid Q1 and sets the groundwork for future studies in the Euclid Wide Survey (EWS), including spectral follow-up analyses and host morphology characterisation.
We present a catalogue of candidate active galactic nuclei (AGN) in the $Euclid$ Quick Release (Q1) fields. For each $Euclid$ source we collect multi-wavelength photometry and spectroscopy information from Galaxy Evolution Explorer (GALEX), $Gaia$, Dark Energy Survey (DES), Wise-field Infrared Survey Explorer (WISE), $Spitzer$, Dark Energy Survey (DESI), and Sloan Digital Sky Survey (SDSS), including spectroscopic redshift from public compilations. We investigate the AGN contents of the Q1 fields by applying selection criteria using $Euclid$ colours and WISE-AllWISE cuts finding respectively 292,222 and 65,131 candidates. We also create a high-purity QSO catalogue based on $Gaia$ DR3 information containing 1971 candidates. Furthermore, we utilise the collected spectroscopic information from DESI to perform broad-line and narrow-line AGN selections, leading to a total of 4392 AGN candidates in the Q1 field. We investigate and refine the Q1 probabilistic random forest QSO population, selecting a total of 180,666 candidates. Additionally, we perform SED fitting on a subset of sources with available $z_{\text{spec}}$, and by utilizing the derived AGN fraction, we identify a total of 7766 AGN candidates. We discuss purity and completeness of the selections and define two new colour selection criteria ($JH$_$I_{\text{E}}Y$ and $I_{\text{E}}H$_$gz$) to improve on purity, finding 313,714 and 267,513 candidates respectively in the Q1 data. We find a total of 229,779 AGN candidates equivalent to an AGN surface density of 3641 deg$^{-2}$ for $18<I_{\text{E}}\leq 24.5$, and a subsample of 30,422 candidates corresponding to an AGN surface density of 482 deg$^{-2}$ when limiting the depth to $18<I_{\text{E}}\leq 22$. The surface density of AGN recovered from this work is in line with predictions based on the AGN X-ray luminosity functions.
Our understanding of cosmic star-formation at $z>3$ used to largely rely on rest-frame UV observations. However, these observations overlook dusty and massive sources, resulting in an incomplete census of early star-forming galaxies. Recently, infrared data from Spitzer and the James Webb Space Telescope (JWST) have revealed a hidden population at $z\sim$3-6 with extreme red colours. Taking advantage of the overlap between imaging in the Euclid Deep Fields (EDFs), covering $\sim$ 60 deg$^2$, and ancillary Spitzer observations, we identified 27000 extremely red objects with $H_E-$IRAC2>2.25 (dubbed HIEROs) down to a $10\sigma$ completeness magnitude limit of IRAC2 $=$ 22.5 AB. After a visual inspection to discard artefacts and objects with troubling photometry, we ended up with a final sample of 3900 candidates. We retrieved the physical parameter estimates for these objects from the SED-fitting tool CIGALE. Our results confirm that HIERO galaxies may populate the high-mass end of the stellar mass function at $z>3$, with some reaching extreme stellar masses ($M_*>10^{11}M_\odot$) and exhibiting high dust attenuation ($A_V>3$). However, we consider stellar mass estimates unreliable for $z>3.5$, favouring a lower-z solution. The challenges faced by SED-fitting tools in characterising these objects highlight the need for further studies, incorporating shorter-wavelength and spectroscopic data. Euclid spectra will help resolve degeneracies and better constrain the physical properties of the brightest galaxies. Given the extreme nature of this population, characterising these sources is crucial for understanding galaxy evolution. This work demonstrates Euclid's potential to provide statistical samples of rare, massive, dust-obscured galaxies at $z>3$, which will be prime targets for JWST, ALMA, and ELT.
Recent James Webb Space Telescope (JWST) observations have revealed a population of sources with a compact morphology and a `v-shaped' continuum, namely blue at rest-frame $\lambda<4000$A and red at longer wavelengths. The nature of these sources, called `little red dots' (LRDs), is still debated, since it is unclear if they host active galactic nuclei (AGN) and their number seems to drastically drop at z<4. We utilise the 63 $deg^2$ covered by the quick Euclid Quick Data Release (Q1) to extend the search for LRDs to brighter magnitudes and to lower z than what has been possible with JWST to have a broader view of the evolution of this peculiar galaxy population. The selection is done by fitting the available photometric data (Euclid, Spitzer/IRAC, and ground-based griz data) with two power laws, to retrieve the rest-frame optical and UV slopes consistently over a large redshift range (i.e, z<7.6). We exclude extended objects and possible line emitters, and perform a visual inspection to remove imaging artefacts. The final selection includes 3341 LRD candidates from z=0.33 to z=3.6, with 29 detected in IRAC. Their rest-frame UV luminosity function, in contrast with previous JWST studies, shows that the number density of LRD candidates increases from high-z down to z=1.5-2.5 and decreases at even lower z. Less evolution is apparent focusing on the subsample of more robust LRD candidates having IRAC detections, which is affected by low statistics and limited by the IRAC resolution. The comparison with previous quasar UV luminosity functions shows that LRDs are not the dominant AGN population at z<4. Follow-up studies of these LRD candidates are key to confirm their nature, probe their physical properties and check for their compatibility with JWST sources, since the different spatial resolution and wavelength coverage of Euclid and JWST could select different samples of compact sources.
The Euclid mission aims to survey around 14000 deg^{2} of extragalactic sky, providing around 10^{5} gravitational lens images. Modelling of gravitational lenses is fundamental to estimate the total mass of the lens galaxy, along with its dark matter content. Traditional modelling of gravitational lenses is computationally intensive and requires manual input. In this paper, we use a Bayesian neural network, LEns MOdelling with Neural networks (LEMON), for modelling Euclid gravitational lenses with a singular isothermal ellipsoid mass profile. Our method estimates key lens mass profile parameters, such as the Einstein radius, while also predicting the light parameters of foreground galaxies and their uncertainties. We validate LEMON's performance on both mock Euclid data sets, real Euclidised lenses observed with Hubble Space Telescope (hereafter HST), and real Euclid lenses found in the Perseus ERO field, demonstrating the ability of LEMON to predict parameters of both simulated and real lenses. Results show promising accuracy and reliability in predicting the Einstein radius, axis ratio, position angle, effective radius, Sérsic index, and lens magnitude for simulated lens galaxies. The application to real data, including the latest Quick Release 1 strong lens candidates, provides encouraging results, particularly for the Einstein radius. We also verified that LEMON has the potential to accelerate traditional modelling methods, by giving to the classical optimiser the LEMON predictions as starting points, resulting in a speed-up of up to 26 times the original time needed to model a sample of gravitational lenses, a result that would be impossible with randomly initialised guesses. This work represents a significant step towards efficient, automated gravitational lens modelling, which is crucial for handling the large data volumes expected from Euclid.
The matter distribution around galaxy clusters is distributed over several filaments, reflecting their positions as nodes in the large-scale cosmic web. The number of filaments connected to a cluster, namely its connectivity, is expected to affect the physical properties of clusters. Using the first Euclid galaxy catalogue from the Euclid Quick Release 1 (Q1), we investigate the connectivity of galaxy clusters and how it correlates with their physical and galaxy member properties. Around 220 clusters located within the three fields of Q1 (covering $\sim 63 \ \text{deg}^2$), are analysed in the redshift range $0.2 < z < 0.7$. Due to the photometric redshift uncertainty, we reconstruct the cosmic web skeleton, and measure cluster connectivity, in 2-D projected slices with a thickness of 170 comoving $h^{-1}.\text{Mpc}$ and centred on each cluster redshift, by using two different filament finder algorithms on the most massive galaxies ($M_*\ > 10^{10.3} \ M_\odot$). In agreement with previous measurements, we recover the mass-connectivity relation independently of the filament detection algorithm, showing that the most massive clusters are, on average, connected to a larger number of cosmic filaments, consistent with hierarchical structure formation models. Furthermore, we explore possible correlations between connectivities and two cluster properties: the fraction of early-type galaxies and the Sérsic index of galaxy members. Our result suggests that the clusters populated by early-type galaxies exhibit higher connectivity compared to clusters dominated by late-type galaxies. These preliminary investigations highlight our ability to quantify the impact of the cosmic web connectivity on cluster properties with Euclid.
We report on serendipitous \Euclid observations of previously known transients, using the \Euclid Q1 data release. By cross-matching with the Transient Name Server (TNS) we identify 164 transients that coincide with the data release. Although the \Euclid Q1 release only includes single-epoch data, we are able to make \Euclid photometric measurements at the location of 161 of these transients. \Euclid obtained deep photometric measurements or upper limits of these transients in the \IE, \YE, \JE, and \HE bands at various phases of the transient light-curves, including before, during, and after the observations of ground-based transient surveys. Approximately 70\% of known transients reported in the six months before the \Euclid observation date and with discovery magnitude brighter than 24 were detected in \Euclid $\IE$ images. Our observations include one of the earliest near-infrared detections of a Type~Ia supernova (SN~2024pvw) 15~days prior to its peak brightness, and the late-phase (435.9~days post peak) observations of the enigmatic core-collapse SN~2023aew. \Euclid deep photometry provides valuable information on the nature of these transients such as their progenitor systems and power sources, with late time observations being a uniquely powerful contribution. In addition, \Euclid is able to detect the host galaxies of some transients that were previously classed as hostless. The Q1 data demonstrate the power of the \Euclid data even with only single-epoch observations available, as will be the case for much larger areas of sky in the Euclid Wide Survey.
The Euclid Q1 fields were selected for calibration purposes in cosmology and are therefore relatively devoid of nearby galaxies. However, this is precisely what makes them interesting fields in which to search for dwarf galaxies in local density environments. We take advantage of the unprecedented depth, spatial resolution, and field of view of the Euclid Quick Release (Q1) to build a census of dwarf galaxies in these regions. We have identified dwarfs in a representative sample of 25 contiguous tiles in the Euclid Deep Field North (EDF-N), covering an area of 14.25 sq. deg. The dwarf candidates were identified using a semi-automatic detection method, based on properties measured by the Euclid pipeline and listed in the MER catalogue. A selection cut in surface brightness and magnitude was used to produce an initial dwarf candidate catalogue, followed by a cut in morphology and colour. This catalogue was visually classified to produce a final sample of dwarf candidates, including their morphology, number of nuclei, globular cluster (GC) richness, and presence of a blue compact centre. We identified 2674 dwarf candidates, corresponding to 188 dwarfs per sq. deg. The visual classification of the dwarfs reveals a slightly uneven morphological mix of 58% ellipticals and 42% irregulars, with very few potentially GC-rich (1.0%) and nucleated (4.0%) candidates but a noticeable fraction (6.9%) of dwarfs with blue compact centres. The distance distribution of 388 (15%) of the dwarfs with spectroscopic redshifts peaks at about 400 Mpc. Their stellar mass distribution confirms that our selection effectively identifies dwarfs while minimising contamination. The most prominent dwarf overdensities are dominated by dEs, while dIs are more evenly distributed. This work highlights Euclid's remarkable ability to detect and characterise dwarf galaxies across diverse masses, distances, and environments.
this https URL paper submitted to the special A&A issue
this https URL . Paper submitted as part of the A&A Special Issue `Euclid Quick Data Release (Q1)'. 16 pages, 15 figures, plus appendices
this https URL . Paper submitted as part of the A&A Special Issue `Euclid Quick Data Release (Q1)'. 20 pages, 11 figures, plus appendices