We implement a detailed dust model into the L-Galaxies semi-analytical model which includes: injection of dust by type II and type Ia supernovae (SNe) and AGB stars; grain growth in molecular clouds; and destruction due to supernova-induced shocks, star formation, and reheating. Our grain growth model follows the dust content in molecular clouds and the inter-cloud medium separately, and allows growth only on pre-existing dust grains. At early times, this can make a significant difference to the dust growth rate. Above $z\sim8$, type II SNe are the primary source of dust, whereas below $z\sim8$, grain growth in molecular clouds dominates, with the total dust content being dominated by the latter below $z\sim6$. However, the detailed history of galaxy formation is important for determining the dust content of any individual galaxy. We introduce a fit to the dust-to-metal (DTM) ratio as a function of metallicity and age, which can be used to deduce the DTM ratio of galaxies at any redshift. At $z\lesssim3$, we find a fairly flat mean relation between metallicity and the DTM, and a positive correlation between metallicity and the dust-to-gas (DTG) ratio, in good agreement with the shape and normalisation of the observed relations. We also match the normalisation of the observed stellar mass -- dust mass relation over the redshift range of $0-4$, and to the dust mass function at $z=0$. Our results are important in interpreting observations on the dust content of galaxies across cosmic time, particularly so at high redshift.

Measurements of the UV-continuum slopes provide valuable information on the physical properties of galaxies forming in the early universe, probing the dust reddening, age, metal content, and even the escape fraction. While constraints on these slopes generally become more challenging at higher redshifts as the UV continuum shifts out of the Hubble Space Telescope bands (particularly at z>7), such a characterisation actually becomes abruptly easier for galaxies in the redshift window z=9.5-10.5 due to the Spitzer/IRAC 3.6um-band probing the rest-UV continuum and the long wavelength baseline between this Spitzer band and the Hubble H-band. Higher S/N constraints on the UV slope are possible at z~10 than at z=8. Here we take advantage of this opportunity and five recently discovered bright z=9.5-10.5 galaxies to present the first measurements of the mean slope for a multi-object sample of galaxy candidates at z~10. We find the measured observed slopes of these candidates are $-2.1\pm0.3\pm0.2$ (random and systematic), only slightly bluer than the measured slopes at 3.5<z<7.5 for galaxies of similar luminosities. Small increases in the stellar ages, metallicities, and dust content of the galaxy population from z~10 to z~7 could easily explain the apparent evolution in slopes.

We present high-redshift predictions of the star-formation-rate distribution function (SFR DF), UV luminosity function (UV LF), galactic stellar mass function (GSMF), and specific star-formation rates (sSFRs) of galaxies from the latest version of the Munich semi-analytic model L-Galaxies. We find a good fit to both the shape and normalisation of the SFR DF at $z=4-7$, apart from a slight under-prediction at the low SFR end at $z=4$. Likewise, we find a good fit to the faint number counts for the observed UV LF; at brighter magnitudes our predictions lie below the observations, increasingly so at higher redshifts. At all redshifts and magnitudes, the raw (unattenuated) number counts for the UV LF lie above the observations. Because of the good agreement with the SFR we interpret our under-prediction as an over-estimate of the amount of dust in the model for the brightest galaxies, especially at high-redshift. While the shape of our GSMF matches that of the observations, we lie between (conflicting) observations at $z=4-5$, and under-predict at $z=6-7$. The sSFRs of our model galaxies show the observed trend of increasing normalisation with redshift, but do not reproduce the observed mass dependence. Overall, we conclude that the latest version of L-Galaxies, which is tuned to match observations at $z\leq3$, does a fair job of reproducing the observed properties of galaxies at $z\geq4$. More work needs to be done on understanding observational bias at high-redshift, and upon the dust model, before strong conclusions can be drawn on how to interpret remaining discrepancies between the model and observations.

We used Herschel Hi-GAL survey data to determine whether massive young stellar objects (MYSOs) are resolved at 70$\mu$m and to study their envelope density distribution. Our analysis of three relatively isolated sources in the l=30{\deg} and l=59{\deg} Galactic fields show that the objects are partially resolved at 70$\mu$m. The Herschel Hi-GAL survey data have a high scan velocity which makes unresolved and partially resolved sources appear elongated in the 70$\mu$m images. We analysed the two scan directions separately and examine the intensity profile perpendicular to the scan direction. Spherically symmetric radiative transfer models with a power law density distribution were used to study the circumstellar matter distribution. Single dish sub-mm data were also included to study how different spatial information affects the fitted density distribution. The density distribution which best fits both the 70$\mu$m intensity profile and SED has an average index of ~0.5. This index is shallower than expected and is probably due to the dust emission from bipolar outflow cavity walls not accounted for in the spherical models. We conclude that 2D axisymmetric models and Herschel images at low scan speeds are needed to better constrain the matter distribution around MYSOs.