Galaxy formation and evolution

Abstract:

The effect of supernovae (SNe) on star formation in the interstellar medium (ISM) depends sensitively on where SNe explode with respect to ISM clouds. Observationally, SN ISM environments characterized by spatially resolved gas maps can empirically guide the placement of SNe in subgrid models, but unfortunately such measurements remain scarce, as SNe are rare and often distant. Here we demonstrate a new approach—mapping the ISM around evolved massive stars that are soon to explode. These provide a substantially larger sample of “explosion sites” (than just historical SNe) in nearby galaxies that have high-resolution atomic and molecular ISM maps from the Jansky Very Large Array and Atacama Large Millimeter/submillimeter Array. We demonstrate this technique in the well-resolved Local Group spiral M33 by analyzing the 50 pc scale projected ISM densities around red supergiants (RSGs; 8–30 M⊙ stars) Wolf–Rayet stars (W-Rs; >30M⊙ stars), and supernova remnants. We find a mass-dependent correlation between stars and gas clouds, with at least 45% of W-Rs and up to 77% of RSGs having no detectable H2 at their pixel locations. In the sample with H2 detections, we find that more-massive younger progenitors are coincident with denser gas. We show that the density distributions for stars >15 M⊙ are statistically distinct from random alignment of stars and gas in M33. Our work provides the first observationally derived estimate of the fraction of the SN-producing stellar population correlated with ISM density peaks. We demonstrate how this can be compared with galaxy simulations, and advocate similar comparisons to the community for constraining subgrid models.

Abstract:

Radio observations offer a dust-independent probe of star formation and active galactic nucleus (AGN) activity, but sufficiently deep data are required to access the cross-over luminosity between these processes at high redshift (). We present three spectroscopically confirmed high-redshift radio sources (HzRSs) detected at 1.3 GHz at –5.6, with radio luminosities spanning –. These sources were first identified as high-redshift candidates through spectral energy distribution (SED) fitting of archival Hubble, James Webb Space Telescope (JWST) NIRCam + MIRI, and ground-based photometry, and then spectroscopically confirmed via the emission line using wide-field slitless spectroscopy from JWST COSMOS-3D. The star formation rates (SFRs) measured from SED fitting, the flux, and the 1.3 GHz luminosity, span –, demonstrating broad agreement between these SFR tracers. We find that these three sources lie either on or 0.5–1.0 dex above the star-forming main sequence at –6 and have undergone a recent burst of star formation. The sources have extended rest-ultraviolet (UV)/optical morphologies with no evidence for a dominant point source component, indicating that an AGN is unlikely to dominate their rest-UV and optical emission. Two of the sources have complex, multicomponent rest-frame UV/optical morphologies, suggesting that their starbursts may be triggered by merging activity. These HzRSs open up a new window towards probing radio emission powered by star formation alone at , representing a remarkable opportunity to begin tracing star formation, independent of dust, in the early Universe.

Abstract:

Radio-AGNs (active galactic nuclei) are observed to be more strongly clustered than non-active galaxies, though it is unclear whether this is simply due to their preference for massive host galaxies, or if they reside in distinct environments beyond this mass dependence. Using data from three fields covered by the MIGHTEE survey, we measure the angular two-point cross-correlation functions with a large, stellar mass-limited population of near-infrared selected galaxies, overcoming limitations of previous single-deep-field studies. By fitting halo occupation distribution models, we infer the galaxy bias parameters, b, for radio-AGN in three redshift ranges with median redshifts of , , and , finding , , and , respectively. The typical dark matter halo mass decreases with increasing redshift: , , and , which we attribute to the increased abundance of cold gas required to fuel AGN activity at earlier times. The AGN duty cycle is determined to be per cent, and we estimate that the total energy radiated by radio-jets over is per halo, which is sufficient to account for the observed excess heating of gas beyond that of gravitational collapse. Comparing the typical dark matter halo masses to the values obtained for the control sample, we find that the halo masses of radio-AGN are , , and times greater than those of the stellar mass- and redshift-matched galaxies. This difference could arise because AGN feedback suppresses stellar mass growth while leaving halo mass unchanged, or because radio-AGN preferentially reside in earlier forming haloes which are more strongly clustered.

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:

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’.