Galaxy formation and evolution

Abstract:

Abstract Symbolic regression (SR) has emerged as a powerful method for uncovering interpretable mathematical relationships from data, offering a novel route to both scientific discovery and efficient empirical modelling. This article introduces the Special issue on symbolic regression for the physical sciences, motivated by the Royal Society discussion meeting held in April 2025. The contributions collected here span applications from automated equation discovery and emergent-phenomena modelling to the construction of compact emulators for computationally expensive simulations. The introductory review outlines the conceptual foundations of SR, contrasts it with conventional regression approaches and surveys its main use cases in the physical sciences, including the derivation of effective theories, empirical functional forms and surrogate models. We summarize methodological considerations such as search-space design, operator selection, complexity control, feature selection and integration with modern AI approaches. We also highlight ongoing challenges, including scalability, robustness to noise, overfitting and computational complexity. Finally, we emphasize emerging directions, particularly the incorporation of symmetry constraints, asymptotic behaviour and other theoretical information. Taken together, the papers in this Special issue illustrate the accelerating progress of SR and its growing relevance across the physical sciences. This article is part of the discussion meeting issue ‘Symbolic regression in the physical sciences’.

Abstract:

Understanding how galaxies assemble their mass during the first billion years of cosmic time is a central goal of extragalactic astrophysics, yet joint constraints on their sizes and kinematics remain scarce. We present one of the first statistical studies of the size–mass relation at high redshift with a sample of 213 galaxies at spectroscopic redshifts of from the FRESCO and CONGRESS NIRCam grism surveys. We measure the morphology and kinematics of our sample using the novel forward modelling Bayesian inference tool geko, and complement them with stellar continuum sizes in the rest-frame far ultraviolet (FUV), near ultraviolet (NUV), and optical, obtained from modelling of imaging data from the JADES survey with Pysersic. At , we find that the average H sizes are larger than the stellar continuum (FUV, NUV, and optical), with kpc and kpc for galaxies with . However, we find no significant differences between the stellar continuum sizes at different wavelengths, suggesting that galaxies are not yet steadily growing inside–out at these epochs. Instead, we find that the ratio increases with the distance above the star-forming main sequence (), consistent with an expansion of H sizes during episodes of enhanced star formation caused by an increase in ionizing photons. As galaxies move above the star-forming main sequence, we find an increase of their rotational support , which could be tracing accreting gas illuminated by the emission. Finally, we find that about half of the elongated systems () are not rotationally supported, indicating a potential flattened/prolate galaxy population at high redshift.

Abstract:

We present a measurement of the Lyα intensity mapping power spectrum from the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). We measure the cross-power spectrum of the Lyα intensity and Lyα-emitting galaxies (LAEs) in a redshift range of 1.9 ≤ z ≤ 3.5. We calculate the intensity from HETDEX spectra that do not contain any detected LAEs above a signal-to-noise ratio of 5.5. To produce a power spectrum model and its covariance matrix, we simulate the data using lognormal mocks for the LAE catalog and Lyα intensity in redshift space. The simulations include the HETDEX sensitivity, selection function, and mask. The measurements yield the product of the LAE bias, the intensity bias, the mean intensity of undetected sources, and the ratio of the actual and fiducial redshift-space distortion parameters, bgbI〈I〉F¯RSD/F¯RSDfid= (6.7 ± 3.1), (11.7 ± 1.4), and (8.3 ± 1.5) × 10−22 erg s−1 cm−2 arcsec−2 Å−1 in three redshift bins centered at z¯=2.1 , 2.6, and 3.2, respectively. The results are reasonably consistent with cosmological hydrodynamical simulations that include Lyα radiative transfer. They are, however, significantly smaller than previous results from cross-correlations of quasars with Lyα intensity. These results demonstrate the statistical power of HETDEX for Lyα intensity mapping and pave the way for a more comprehensive analysis. They will also be useful for constraining models of Lyα emission from galaxies used in modern cosmological simulations of galaxy formation and evolution.

Abstract:

We present a search for z ≃ 7 Lyman-break galaxies using the 1.72 deg 2 near-infrared (NIR) UltraVISTA survey in the COSMOS field, reaching 5 σ depths in Y of 26.2. We incorporated deep Euclid optical and Euclid + Spitzer NIR imaging for a full spectral energy distribution (SED) fitting analysis. We found 289 candidate galaxies at 6.5 ≤ z ≤ 7.5 covering −22.6 ≤ M UV ≤ −20.2, faint enough to overlap with Hubble Space Telescope studies. We conducted a separate selection by including complementary Euclid performance verification imaging (reaching 5 σ depths of 26.3), yielding 140 galaxies in 0.65 deg 2 , with 38 sources unique to this sample. We computed the rest-frame UV luminosity function (UV LF) from our samples, extending below the knee ( M ∗ = 21.14 +0.28 −0.25 ). We find that the shape of the UV LF is consistent with both a Schechter function and a double power law (DPL) at the magnitudes probed by this sample, with a DPL preferred at M UV < −22.5 when bright-end results are included. The UltraVISTA + Euclid sample provides a clean measurement of the LF due to the overlapping NIR filters identifying molecular absorption features in the SEDs of ultra-cool dwarf interlopers, and additional faint galaxies were recovered. A comparison with JWST LFs at z > 7 suggests a gentle evolution in the bright-end slope, although this is limited by a lack of robust bright-end measurements at z > 9. We forecast that in the Euclid Deep Fields, the removal of contaminant ultra-cool dwarfs as point sources will be possible at J E < 24.5. Finally, we present a high-equivalent-width Lyman- α emitter candidate identified by combining HSC, VISTA, and Euclid broadband photometry, highlighting the synergistic power these instruments will have in the Euclid Auxiliary Fields for identifying extreme sources in the epoch of reionisation.

Abstract:

Context. This paper presents new results from the ESO Large Programme Looking into the faintEst WIth MUSE (LEWIS). The LEWIS sample consists of low-surface-brightness galaxies (LSBs) and ultra-diffuse galaxies (UDGs) located inside 0.4 R vir of the Hydra I cluster. Integral-field spectroscopy is acquired for 24 galaxies with the MUSE spectrograph mounted on the Very Large Telescope (VLT). Aims. Our main objectives are to analyse possible correlations between the environment and the integrated stellar population properties of our targets, based on which we infer clues about their formation. Methods. For each galaxy in the sample, we extracted the 1D stacked spectrum in an aperture of one effective radius R e and adopted previously published stellar kinematics to derive the age, metallicity, and [Mg/Fe] through a full spectral fitting technique. Results. We find that the analysed LEWIS sample has a mean metallicity of ⟨[M/H]⟩= − 0.9 ± 0.2 dex and a mean age of 10 ± 2 Gyr, comparable to previous results of UDGs in other clusters. According to their position in the projected phase space, galaxies can be classified into two groups: very early infaller galaxies, which on average have slightly higher metallicities (⟨[M/H]⟩ early = −0.8 ± 0.1 dex), and late infaller galaxies, with slightly lower values (⟨[M/H]⟩ late = −1.0 ± 0.1 dex). According to their properties, late-infallers tend to be rotation-supported systems. Conversely, two types of galaxies are found in the early-infall region. Roughly half have metallicities consistent with the dwarf galaxy mass–metallicity relation. The other half show higher metallicities (with ⟨[M/H]⟩≥ − 1.0 dex) and are located outside the 1 σ scatter of the mass-metallicity relation. The two subgroups of early-infallers also display different timescales for stellar mass assembly. Metal-rich galaxies reached 50% of their stellar mass in less than 1 Gyr and show a prolonged and almost constant star formation over more than 12 Gyr. The other galaxies exhibit a star formation history similar to that found for galaxies in the late-infall region. Both early and late infallers show solar-like α abundances. Conclusions. From the analysis of stellar population properties presented in this work and of stellar kinematics previously obtained from LEWIS, we identified different classes of UDGs within the Hydra I cluster – as shown by metallicities, quenching timescales, and kinematics – which suggest different formation mechanisms. Almost all of the UDGs and LSBs in this cluster are consistent with the puffed-up dwarf formation scenario, having dwarf-like metallicities and being consistent with the mass-metallicity relation for dwarfs. In the innermost regions of the cluster, where more metal-rich UDGs are found, tidal effects or the environment might have influenced their formation and evolution.