JWST has identified a large population of faint, broad-line active galactic nuclei (AGN) in the early universe that are powered by black holes (BHs) that often appear overmassive relative to their host galaxies. In this study, we examine the relationship between BH mass and galaxy stellar mass at $33\sigma$ above the relationship measured for local broad-line AGN. We derive an intrinsic scatter in this relationship of $0.9$ dex, which does not vary over the redshift range of our sample. We also find that the $M_{\rm BH}/M_{\star}$ ratio increases by $2.3$ dex from $z = 3.5$ and $z = 6.5$ with a confidence level of > 3\sigma. We attribute this trend with the increasing fraction of LRDs in our sample at z>4 as their host masses are $\sim1$ dex lower than the non-LRD AGN in our sample. These results support a picture in which the BHs powering JWST's broad-line AGN are genuinely overmassive and become increasingly so with redshift. We discuss the implications of our findings on early BH growth relative to that of their host galaxies and the constraints it places on BH seeding models.

We used the spectro-photometric information of ~219 million stars from Gaia's DR3 to calculate synthetic, narrow-band, metallicity-sensitive CaHK magnitudes that mimic the observations of the Pristine survey, a survey of photometric metallicities of Milky Way (MW) stars that has been mapping >6,500 deg^2 of the northern sky with CFHT since 2015. These synthetic magnitudes were used for an absolute recalibration of the deeper Pristine photometry and, combined with broadband Gaia information, synthetic and Pristine CaHK magnitudes were used to estimate photometric metallicities over the whole sky. The resulting metallicity catalogue is accurate down to [Fe/H]~-3.5 and is particularly suited for the exploration of the metal-poor MW ([Fe/H]<-1.0). We make available here the catalogue of synthetic CaHK_syn magnitudes for all stars with BP/RP information in Gaia DR3, as well as an associated catalogue of more than ~30 million photometric metallicities for high S/N FGK stars. This paper further provides the first public data release of the Pristine catalogue in the form of higher quality recalibrated Pristine CaHK magnitudes and photometric metallicities for all stars in common with the BP/RP spectro-photometric information in Gaia DR3. When available, the much deeper Pristine data greatly enhance the quality of the derived metallicities, in particular at the faint end of the catalogue (G_BP > 16). Combined, both photometric metallicity catalogues include more than two million metal-poor star candidates ([Fe/H]_phot<-1.0) as well as more than 200,000 and ~8,000 very and extremely metal-poor candidates ([Fe/H]_phot<-2.0 and <-3.0, respectively). Finally, we show that these metallicity catalogues can be used efficiently, among other applications to hunt for the most metal-poor stars, and to study how the structure of the MW varies with metallicity.

H2 is the most abundant molecule in the interstellar medium and is a useful tool to study photodissociation regions, where radiative feedback from massive stars on molecular clouds is dominant. The James Webb Space Telescope, with its high spatial resolution, sensitivity, and wavelength coverage provides unique access to the detection of most of H2 lines and the analysis of its spatial morphology. Our goal is to use H2 line emission detected with the JWST in the Horsehead nebula to constrain the physical parameters (e.g., extinction, gas temperature, thermal pressure) throughout the PDR and its geometry. The study of H2 morphology reveals that FUV-pumped lines peak closer to the edge of the PDR than thermalized lines. From H2 lines, we estimate the value of extinction throughout the PDR. We find that AV is increasing from the edge of the PDR to the second and third H2 filaments. Then, we study the H2 excitation in different regions across the PDR. The temperature profile shows that the observed gas temperature is quite constant throughout the PDR, with a slight decline in each of the dissociation fronts. This study also reveals that the OPR is far from equilibrium. We observe a spatial separation of para and ortho rovibrational levels, indicating that efficient ortho-para conversion and preferential ortho self-shielding are driving the spatial variations of the OPR. Finally, we derive a thermal pressure in the first filament around P > 6x10$^6$ K cm$^{-3}$, about ten times higher than that of the ionized gas. We highlight that template stationary 1D PDR models cannot account for the intrinsic 2D structure and the very high temperature observed in the Horsehead nebula. We argue the highly excited, over-pressurized H2 gas at the edge of the PDR interface could originate from the mixing between the cold and hot phase induced by the photo-evaporation of the cloud.

The C-19 stream is the most metal poor stellar system ever discovered, with a mean metallicity $[Fe/H] = -3.38\pm0.06$. Its low metallicity dispersion ($\sigma_{\rm [Fe/H]}$ &lt; 0.18 at the 95\% confidence level) as well as variations in sodium abundances strongly suggest a globular cluster origin. In this work, we use VLT/UVES spectra of seven C-19 stars to derive more precise velocity measurements for member stars, and to identify two new members with radial velocities and metallicities consistent with the stream's properties. One of these new member stars is located 30 $\deg$ away from the previously identified body of C-19, implying that the stream is significantly more extended than previously known and that more members likely await discovery. In the main part of C-19, we measure a radial velocity dispersion $\sigma_v$ = 6.2$^{+2.0}_{-1.4}$ km s$^{-1}$ from nine members, and a stream width of 0.56$\deg\pm0.08\deg$, equivalent to $\sim$158 pc at a heliocentric distance of 18 kpc. These confirm that C-19 is comparatively hotter, dynamically, than other known globular cluster streams and shares the properties of faint dwarf galaxy streams. On the other hand, the variations in the Na abundances of the three newly observed bright member stars, the variations in Mg and Al for two of them, and the normal Ba abundance of the one star where it can be measured provide further evidence for a globular cluster origin. The tension between the dynamical and chemical properties of C-19 suggests that its progenitor experienced a complex birth environment or disruption history.

The Mid-Infrared Instrument (MIRI) on the James Webb Space Telescope (JWST), has imaging, four coronagraphs and both low and medium resolution spectroscopic modes . Being able to simulate MIRI observations will help commissioning of the instrument, as well as get users familiar with representative data. We designed the MIRI instrument simulator (MIRISim) to mimic the on-orbit performance of the MIRI imager and spectrometers using the Calibration Data Products (CDPs) developed by the MIRI instrument team. The software encorporates accurate representations of the detectors, slicers, distortions, and noise sources along the light path including the telescope's radiative background and cosmic rays. The software also includes a module which enables users to create astronomical scenes to simulate. MIRISim is a publicly available Python package that can be run at the command line, or from within Python. The outputs of MIRISim are detector images in the same uncalibrated data format that will be delivered to MIRI users. These contain the necessary metadata for ingestion by the JWST calibration pipeline.

The 2017 observing campaign of the Event Horizon Telescope (EHT) delivered the first very long baseline interferometry (VLBI) images at the observing frequency of 230 GHz, leading to a number of unique studies on black holes and relativistic jets from active galactic nuclei (AGN). In total, eighteen sources were observed: the main science targets, Sgr A* and M87 along with various calibrators. We investigated the morphology of the sixteen AGN in the EHT 2017 data set, focusing on the properties of the VLBI cores: size, flux density, and brightness temperature. We studied their dependence on the observing frequency in order to compare it with the Blandford-Königl (BK) jet model. We modeled the source structure of seven AGN in the EHT 2017 data set using linearly polarized circular Gaussian components and collected results for the other nine AGN from dedicated EHT publications, complemented by lower frequency data in the 2-86 GHz range. Then, we studied the dependences of the VLBI core flux density, size, and brightness temperature on the frequency measured in the AGN host frame. We compared the observations with the BK jet model and estimated the magnetic field strength dependence on the distance from the central black hole. Our results indicate a deviation from the standard BK model, particularly in the decrease of the brightness temperature with the observing frequency. Either bulk acceleration of the jet material, energy transfer from the magnetic field to the particles, or both are required to explain the observations.

Carbon is an essential element for life but how much can be delivered to young planets is still an open question. The chemical characterization of planet-forming disks is a crucial step in our understanding of the diversity and habitability of exoplanets. Very low-mass stars (&lt;0.2~M_{\odot}) are interesting targets because they host a rich population of terrestrial planets. Here we present the JWST detection of abundant hydrocarbons in the disk of a very low-mass star obtained as part of the MIRI mid-INfrared Disk Survey (MINDS). In addition to very strong and broad emission from C$_2$H$_2$ and its $^{13}$C$^{12}$CH$_2$ isotopologue, C$_4$H$_2$, benzene, and possibly CH$_4$ are identified, but water, PAH and silicate features are weak or absent. The lack of small silicate grains implies that we can look deep down into this disk. These detections testify to an active warm hydrocarbon chemistry with a high C/O ratio in the inner 0.1 au of this disk, perhaps due to destruction of carbonaceous grains. The exceptionally high C$_2$H$_2$/CO$_2$ and C$_2$H$_2$/H$_2$O column density ratios suggest that oxygen is locked up in icy pebbles and planetesimals outside the water iceline. This, in turn, will have significant consequences for the composition of forming exoplanets.

The embedded phase is a crucial period in the development of a young star. Mid-IR observations, now possible with JWST with unprecedented sensitivity, spectral resolution and sharpness are key for probing many physical and chemical processes on sub-arcsecond scales. JOYS addresses a wide variety of questions, from protostellar accretion and the nature of primeval jets, winds and outflows, to the chemistry of gas and ice, and the characteristics of embedded disks. We introduce the program and show representative results. MIRI-MRS data of 17 low-mass and 6 high-mass protostars show a wide variety of features. Atomic line maps differ among refractory (e.g., Fe), semi-refractory (e.g., S) and volatile elements (e.g., Ne), linked to their different levels of depletion and local (shock) conditions. Nested, stratified jet structures consisting of an inner ionized core seen in [Fe II] with an outer H2 layer are commonly seen. Wide-angle winds are found in low-J H2 lines. [S I] follows the jet in the youngest protostars, but is concentrated on source when more evolved. [Ne II] reveals a mix of jet shock and photoionized emission. H I lines measure accretion, but are also associated with jets. Molecular emission (CO2, C2H2, HCN, H2O, ..) is cool compared with disks, and likely associated with hot cores. Deep ice absorption features reveal not just the major ice components but also ions (as part of salts) and complex organic molecules, with comparable abundances from low- to high-mass sources. A second detection of HDO ice in a solar-mass source is presented with HDO/H2O ~ 0.4%, providing a link with disks and comets. A deep search for solid O2 suggests it is not a significant oxygen reservoir. Only few embedded Class I disks show the same forest of water lines as Class II disks do, perhaps due to significant dust extinction of the upper layers [abridged].

Fundamental questions about the physics responsible for fragmenting molecular parsec-scale clumps into cores of ~1000 au are still open, that only a statistically significant investigation with ALMA is able to address: what are the dominant agents that determine the core demographics, mass, and spatial distribution as a function of the physical properties of the hosting clumps, their evolutionary stage and the different Galactic environments in which they reside? To what extent extent is fragmentation driven by clumps dynamics or mass transport in filaments? With ALMAGAL we observed the 1.38 mm continuum and lines toward more than 1000 dense clumps in our Galaxy, with M>500M_sun, surface density > 0.1 g/cm2 and d<7.5 kpc. The ACA and two 12-m array setups were used to deliver a minimum resolution of ~1000 au over the entire sample distance range. The sample covers all evolutionary stages from infrared dark clouds (IRDCs) to HII regions from the tip of the Galactic bar to the outskirts of the Galaxy. The spectral setup includes several molecular lines to trace the multiscale physics and dynamics of gas, notably CH3CN, H2CO, SiO, CH3OH, DCN, HC3N, SO etc. We present an initial overview of the observations and the early science product and results, with a first characterization of the morphological properties of the continuum emission. We use "perimeter-versus-area" and convex hull-versus-area metrics to classify the different morphologies. More extended and morphologically complex shapes are found toward clumps that are relatively more evolved and have higher surface densities.

The detection of low surface brightness galaxies beyond the Local Group poses significant observational challenges, yet these faint systems are fundamental to our understanding of dark matter, hierarchical galaxy formation, and cosmic structure. Their abundance and distribution provide crucial tests for cosmological models, particularly regarding the small-scale predictions of $\Lambda$CDM. We present a systematic detection framework for dwarf galaxy candidates in Ultraviolet Near Infrared Optical Northern Survey (UNIONS) data covering 4,861 deg$^{2}$. Our pipeline preprocesses UNIONS gri-band data through binning, artifact removal, and stellar masking, then employs MTObjects (MTO) for low surface brightness detection. After parameter cuts and cross-matching, we obtain $\sim$360 candidates per deg$^{2}$, totaling $\sim$1.5 million candidates forming our GOBLIN (Galaxies OBserved as Low-luminosity Identified Nebulae) catalog. We fine-tuned the deep learning model Zoobot, pre-trained on Galaxy Zoo labels, for classification. Training data came from visual inspection of literature candidates with probability labels from expert assessments, capturing consensus and uncertainty. Applied to all MTO objects, our method identifies 42,965 dwarf candidates with probability $>$ 0.8, including 23,072 with probability $>$ 0.9. High-probability candidates correlate spatially with massive galaxies (log$(M_{*}/M_{\odot}) \geq$ 10) within 120 Mpc. While some of these objects may have been previously identified in other surveys, we present this extensive catalog of candidates, including their positions, structural parameter estimates, and classification probabilities, as a resource for the community to enable studies of galaxy formation, evolution, and the distribution of dwarf galaxies in different environments.

Retrieval methods are a powerful analysis technique for modelling exoplanetary atmospheres by estimating the bulk physical and chemical properties that combine in a forward model to best-fit an observed spectrum, and they are increasingly being applied to observations of directly-imaged exoplanets. We have adapted TauREx3, the Bayesian retrieval suite, for the analysis of near-infrared spectrophotometry from directly-imaged gas giant exoplanets and brown dwarfs. We demonstrate TauREx3's applicability to sub-stellar atmospheres by presenting results for brown dwarf benchmark GJ 570D which are consistent with previous retrieval studies, whilst also exhibiting systematic biases associated with the presence of alkali lines. We also present results for the cool exoplanet 51 Eri b, the first application of a free chemistry retrieval analysis to this object, using spectroscopic observations from GPI and SPHERE. While our retrieval analysis is able to explain spectroscopic and photometric observations without employing cloud extinction, we conclude this may be a result of employing a flexible temperature-pressure profile which is able to mimic the presence of clouds. We present Bayesian evidence for an ammonia detection with a 2.7$\sigma$ confidence, the first indication of ammonia in an exoplanetary atmosphere. This is consistent with this molecule being present in brown dwarfs of a similar spectral type. We demonstrate the chemical similarities between 51 Eri b and GJ 570D in relation to their retrieved molecular abundances. Finally, we show that overall retrieval conclusions for 51 Eri b can vary when employing different spectral data and modelling components, such as temperature-pressure and cloud structures.

ERIS (Enhanced Resolution Imager and Spectrograph) is a new adaptive optics instrument installed at the Cassegrain focus of the VLT-UT4 telescope at the Paranal Observatory in Chile. ERIS consists of two near-infrared instruments: SPIFFIER, an integral field unit (IFU) spectrograph covering J to K bands, and NIX, an imager covering J to M bands. ERIS has an adaptive optics system able to work with both LGS and NGS. The Assembly Integration Verification (AIV) phase of ERIS at the Paranal Observatory was carried out starting in December 2021, followed by several commissioning runs in 2022. This contribution will describe the first preliminary results of the on-sky performance of ERIS during its commissioning and the future perspectives based on the preliminary scientific results.

The Medium-Resolution Spectrometer (MRS) provides one of the four operating modes of the Mid-Infrared Instrument (MIRI) on board the James Webb Space Telescope (JWST). The MRS is an integral field spectrometer, measuring the spatial and spectral distributions of light across the 5-28 $\mu m$ wavelength range with a spectral resolving power between 3700-1300. We present the MRS's optical, spectral, and spectro-photometric performance, as achieved in flight, and we report on the effects that limit the instrument's ultimate sensitivity. The MRS flight performance has been quantified using observations of stars, planetary nebulae, and planets in our Solar System. The precision and accuracy of this calibration was checked against celestial calibrators with well-known flux levels and spectral features. We find that the MRS geometric calibration has a distortion solution accuracy relative to the commanded position of 8 mas at 5 $\mu m$ and 23 mas at 28 $\mu m$. The wavelength calibration is accurate to within 9 km/sec at 5 $\mu m$ and 27 km/sec at 28 $\mu m$. The uncertainty in the absolute spectro-photometric calibration accuracy was estimated at 5.6 +- 0.7 %. The MIRI calibration pipeline is able to suppress the amplitude of spectral fringes to below 1.5 % for both extended and point sources across the entire wavelength range. The MRS point spread function (PSF) is 60 % broader than the diffraction limit along its long axis at 5 $\mu m$ and is 15 % broader at 28 $\mu m$. The MRS flight performance is found to be better than prelaunch expectations. The MRS is one of the most subscribed observing modes of JWST and is yielding many high-profile publications. It is currently humanity's most powerful instrument for measuring the mid-infrared spectra of celestial sources and is expected to continue as such for many years to come.

Using NIRSpec on JWST, we studied a sample of 15 intermediate-mass (1.8-4.1 Msun) young stellar objects (YSOs) previously identified with MIRI photometry in the low-metallicity NGC 346 star-forming cluster in the Small Magellanic Cloud (SMC). All objects, observed in the 1.7-5.3 micron range, show strong hydrogen recombination lines in the Paschen, Brackett, Pfund, and Humphreys series, confirming their very young ages. The spectra of 11 YSOs show prominent absorption bands from the three most important ice species (H2O, CO2, CO), marking the first detection of these ices in intermediate-mass YSOs beyond our Galaxy. In three YSOs, water ice appears to be in crystalline form. In some objects, we also detect 13CO2 and OCS ices -- never before observed beyond the Milky Way (MW) -- and methanol ice in at least one star. We compared the column densities of H2O, CO2, and CO ices with those measured in more and less massive protostars in the MW and Large Magellanic Cloud, finding that in NGC 346 ice column densities reach values nearly an order of magnitude lower than in more massive objects (~1x10^{17} cm-2 for water and ~1x10^{16} cm-2 for CO2 and CO). However, the relative proportions of the ice species abundances do not differ from those in massive MW YSOs. This suggests that metallicity may not significantly affect ice chemistry in protoplanetary discs and that, shielded by the protostellar envelope or deep in the midplane, circumstellar material is likely impervious to the radiation environment.

MOONS is the new Multi-Object Optical and Near-infrared Spectrograph currently under construction for the Very Large Telescope (VLT) at ESO. This remarkable instrument combines, for the first time, the collecting power of an 8-m telescope, 1000 fibres with individual robotic positioners, and both low- and high-resolution simultaneous spectral coverage across the 0.64-1.8 micron wavelength range. This facility will provide the astronomical community with a powerful, world-leading instrument able to serve a wide range of Galactic, extragalactic and cosmological studies. Construction is now proceeding full steam ahead and this overview article presents some of the science goals and the technical description of the MOONS instrument. More detailed information on the MOONS surveys is provided in the other dedicated articles in this Messenger issue.

We examine the distribution of radio emission from ~42,000 quasars from the Sloan Digital Sky Survey, as measured in the LOFAR Two-Metre Sky Survey (LoTSS). We present a model of the radio luminosity distribution of the quasars that assumes that every quasar displays a superposition of two sources of radio emission: active galactic nuclei (jets) and star-formation. Our two-component model provides an excellent match to the observed radio flux density distributions across a wide range of redshifts and quasar optical luminosities; this suggests that the jet-launching mechanism operates in all quasars but with different powering efficiency. The wide distribution of jet powers allows for a smooth transition between the 'radio-quiet' and 'radio-loud' quasar regimes, without need for any explicit bimodality. The best-fit model parameters indicate that the star-formation rate of quasar host galaxies correlates strongly with quasar luminosity and also increases with redshift at least out to z~2. For a model where star-formation rate scales as $SFR \propto L_{bol}^\alpha (1+z)^\beta$, we find $\alpha = 0.47 \pm 0.01$ and $\beta = 1.61 \pm 0.05$, in agreement with far-infrared studies. Quasars contribute ~0.15 per cent of the cosmic star-formation rate density at z=0.5, rising to 0.4 per cent by z=2. The typical radio jet power is seen to increase with both increasing optical luminosity and black hole mass independently, but does not vary with redshift, suggesting intrinsic properties govern the production of the radio jets. We discuss the implications of these results for the triggering of quasar activity and the launching of jets.

In the era of JWST, with its unprecedented sensitivity and spectral resolution, infrared spectral surveys have revealed a rich inventory of ices, including complex organic molecules (COMs), in young stellar objects (YSOs). However, robust methods to decompose and quantify these absorption features particularly across broad spectral ranges, are still under investigation. We present INDRA (Ice-fitting with NNLS-based Decomposition and Retrieval Algorithm), a fully Python-based tool that performs continuum and silicate removal, global ice fitting using Weighted Non-Negative Least Squares (NNLS), and estimates column densities and statistical significance. We apply INDRA to NGC 1333 IRAS 2A, a target from the JWST Observations of Young protoStars (JOYS+) program previously studied using local fitting. We derive optical depths via polynomial continuum subtraction and remove silicate absorption using a synthetic model, isolating ice features for global MIRI fitting. Our results are consistent with previous local fits, confirming simple species and COMs, and expand the inventory by identifying additional absorption features from CO2 and NH4+. We also propose the presence of organic refractories contributing up to 9.6% in the spectral region of 5-8 microns among the various ice components, whose inclusion significantly improves the global spectral fitting. These broad absorption features, extending across 5.5-11 microns, are likely produced by large, complex molecules containing carbonyl (C=O), hydroxyl (O-H), amine (N-H), and C-H bending modes. Our expanded inventory, now incorporating these organic residues, offers new insights into the chemical evolution of ices in star-forming regions and highlights the importance of global spectral fitting in constraining ice compositions.

Ultra diffuse galaxies, characterized by their low surface brightness and large physical size, constitute a subclass of dwarf galaxies that challenge our current understanding of galaxy formation and evolution. In this paper, we probe the properties of 74 UDGs, identified in the MATLAS survey, based on a comprehensive study of their globular cluster (GC) populations. We obtained high resolution HST imaging of these galaxies using the ACS F606W and F814W filters, allowing us to select GCs based on color and concentration index. After background subtraction and completeness correction, we calculate an overall total of 387 GCs. The number of GCs per galaxy ranges from 0 to 38, with the majority (64%) having low counts (0-2 GCs). On average, the more massive UDGs host a larger number of GCs. We find that our UDGs have specific frequencies (S_N) ranging from 0 to 91, with a small population (9%) with S_N > 30. The median S_N of our sample is similar to the one for the Perseus cluster UDGs, despite the fact that our UDGs are found in lower density environments. The S_N measurements for individual galaxies can extend beyond those found in Perseus, but remain below the values found for UDGs in the Virgo and Coma cluster. Based on a trending analysis of the S_N values with the host galaxy properties, we find trends with host galaxy size, roundness, color, and local density. For the UDGs with sufficiently high statistics, we study 2D density maps of the GC distributions, which show a variety of appearances: symmetric, asymmetric, off-center, and elongated. The UDGs with disturbed density maps also show disturbed stellar light morphologies. We further quantify the distribution by modeling it with a Sersic profile, finding R_{e,GC}/R_{e,gal} ~ 1.0, which indicates that the GCs follow the stellar light of the host galaxy.

Our current understanding is that intermediate- to high-mass stars form in a way similar to low-mass stars, that is, through disk accretion. However, the physical conditions that play a role in disk formation, evolution, and the possibility of (sub)stellar companion formation, are significantly different. We search for the mm counterparts of four intermediate- to high-mass (4-10 Solar mass) young stellar objects (YSOs) in the giant Hii region M17 at a distance of 1.7 kpc. These objects expose their photospheric spectrum such that their location on the pre-main-sequence (PMS) is well established. They have a circumstellar disk that is likely remnant of the formation process. With ALMA we have detected, for the first time, these four YSOs in M17, in Band 6 and 7, as well as four other serendipitous objects. Besides the flux measurements, the source size and spectral index provide important constraints on the physical mechanism(s) producing the observed emission. We apply different models to estimate the dust and gas mass contained in the disks. All our detections are spatially unresolved, constraining the source size to <120 au, and have a spectral index in the range 0.5-2.7. The derived (upper limits on the) disk dust masses are on the order of a few Earth masses and estimations of the upper limits on the gas mass vary between $10^{-5}$ and $10^{-3}$ Solar mass. In two objects (B331 and B268) free-free emission indicates the presence of ionized material around the star. The four serendipitous detections are likely (low-mass) YSOs. We compare the derived disk masses of our M17 targets to those obtained for YSOs in low-mass star-forming regions (SFRs) and Herbig stars, as a function of stellar mass, age, luminosity, and outer disk radius. The M17 sample, though small, is both the most massive and the youngest sample, yet has the lowest mean disk mass.

Source classifications, stellar masses and star formation rates are presented for 80,000 radio sources from the first data release of the Low Frequency Array Two-metre Sky Survey (LoTSS) Deep Fields, which represents the widest deep radio survey ever undertaken. Using deep multi-wavelength data spanning from the ultraviolet to the far-infrared, spectral energy distribution (SED) fitting is carried out for all of the LoTSS-Deep host galaxies using four different SED codes, two of which include modelling of the contributions from an active galactic nucleus (AGN). Comparing the results of the four codes, galaxies that host a radiative AGN are identified, and an optimised consensus estimate of the stellar mass and star-formation rate for each galaxy is derived. Those galaxies with an excess of radio emission over that expected from star formation are then identified, and the LoTSS-Deep sources are divided into four classes: star-forming galaxies, radio-quiet AGN, and radio-loud high-excitation and low-excitation AGN. Ninety-five per cent of the sources can be reliably classified, of which more than two-thirds are star-forming galaxies, ranging from normal galaxies in the nearby Universe to highly-starbursting systems at z>4. Star-forming galaxies become the dominant population below 150-MHz flux densities of about 1 mJy, accounting for 90 per cent of sources at a 150-MHz flux density of 100 microJy. Radio-quiet AGN comprise around 10 per cent of the overall population. Results are compared against the predictions of the SKADS and T-RECS radio sky simulations, and improvements to the simulations are suggested.