No Arabic abstract
We examine the star formation histories (SFHs) of galaxies in smoothed particle hydrodynamics (SPH) simulations, compare them to parametric models that are commonly used in fitting observed galaxy spectral energy distributions, and examine the efficacy of these parametric models as practical tools for recovering the physical parameters of galaxies. The commonly used tau-model, with SFR ~ exp(-t/tau), provides a poor match to the SFH of our SPH galaxies, with a mismatch between early and late star formation that leads to systematic errors in predicting colours and stellar mass-to-light ratios. A one-parameter lin-exp model, with SFR ~ t*exp(-t/tau), is much more successful on average, but it fails to match the late-time behavior of the bluest, most actively star-forming galaxies and the passive, red and dead galaxies. We introduce a 4-parameter model, which transitions from lin-exp to a linear ramp after a transition time, which describes our simulated galaxies very well. We test the ability of these parametrised models to recover (at z=0, 0.5, and 1) the stellar mass-to-light ratios, specific star formation rates, and stellar population ages from the galaxy colours, computed from the full SPH star formation histories using the FSPS code of Conroy et al. (2009). Fits with tau-models systematically overestimate M/L by ~ 0.2 dex, overestimate population ages by ~ 1-2 Gyr, and underestimate sSFR by ~ 0.05 dex. Fits with lin-exp are less biased on average, but the 4-parameter model yields the best results for the full range of galaxies. Marginalizing over the free parameters of the 4-parameter model leads to slightly larger statistical errors than 1-parameter fits but essentially removes all systematic biases, so this is our recommended procedure for fitting real galaxies.
We adapt the L-Galaxies semi-analytic model to follow the star-formation histories (SFH) of galaxies -- by which we mean a record of the formation time and metallicities of the stars that are present in each galaxy at a given time. We use these to construct stellar spectra in post-processing, which offers large efficiency savings and allows user-defined spectral bands and dust models to be applied to data stored in the Millennium data repository. We contrast model SFHs from the Millennium Simulation with observed ones from the VESPA algorithm as applied to the SDSS-7 catalogue. The overall agreement is good, with both simulated and SDSS galaxies showing a steeper SFH with increased stellar mass. The SFHs of blue and red galaxies, however, show poor agreement between data and simulations, which may indicate that the termination of star formation is too abrupt in the models. The mean star-formation rate (SFR) of model galaxies is well-defined and is accurately modelled by a double power law at all redshifts: SFR proportional to $1/(x^{-1.39}+x^{1.33})$, where $x=(t_a-t)/3.0,$Gyr, $t$ is the age of the stars and $t_a$ is the loopback time to the onset of galaxy formation; above a redshift of unity, this is well approximated by a gamma function: SFR proportional to $x^{1.5}e^{-x}$, where $x=(t_a-t)/2.0,$Gyr. Individual galaxies, however, show a wide dispersion about this mean. When split by mass, the SFR peaks earlier for high-mass galaxies than for lower-mass ones, and we interpret this downsizing as a mass-dependence in the evolution of the quenched fraction: the SFHs of star-forming galaxies show only a weak mass dependence.
We introduce ProSpect, a generative galaxy spectral energy distribution (SED) package that encapsulates the best practices for SED methodologies in a number of astrophysical domains. ProSpect comes with two popular families of stellar population libraries (BC03 and EMILES), and a large variety of methods to construct star formation and metallicity histories. It models dust through the use of a Charlot & Fall attenuation model, with re-emission using Dale far-infrared templates. It also has the ability to model AGN through the inclusion of a simple AGN and hot torus model. Finally, it makes use of MAPPINGS-III photoionisation tables to produce line emission features. We test the generative and inversion utility of ProSpect through application to the Shark galaxy formation semi-analytic code, and informed by these results produce fits to the final ultraviolet to far-infrared photometric catalogues produces by the Galaxy and Mass Assembly Survey (GAMA). As part of the testing of ProSpect, we also produce a range of simple photometric stellar mass approximations covering a range of filters for both observed frame and rest frame photometry.
We study the links between star formation history and structure for a large mass-selected galaxy sample at 0.05 < z_phot < 0.30. The galaxies inhabit a very broad range of environments, from cluster cores to the field. Using HST images, we quantify their structure following Hoyos et al. (2012), and divide them into disturbed and undisturbed. We also visually identify mergers. Additionally, we provide a quantitative measure of the degree of disturbance for each galaxy (roughness). The majority of elliptical and lenticular galaxies have relaxed structure, showing no signs of ongoing star formation. Structurally-disturbed galaxies, which tend to avoid the lowest-density regions, have higher star-formation activity and younger stellar populations than undisturbed systems. Cluster spirals with reduced/quenched star formation have somewhat less disturbed morphologies than spirals with normal star-formation activity, suggesting that these passive spirals have started their morphological transformation into S0s. Visually identified mergers and galaxies not identified as mergers but with similar roughness have similar specific star formation rates and stellar ages. The degree of enhanced star formation is thus linked to the degree of structural disturbance, regardless of whether it is caused by major mergers or not. This suggests that merging galaxies are not special in terms of their higher-than-normal star-formation activity. Any physical process that produces roughness, or regions of enhanced luminosity density, will increase the star-formation activity in a galaxy with similar efficiency. An alternative explanation is that star formation episodes increase the galaxies roughness similarly, regardless of whether they are merger-induced or not.
We discuss the spectral energy distributions and physical properties of six galaxies whose photometric redshifts suggest they lie beyond a redshift $zsimeq$9. Each was selected on account of a prominent excess seen in the Spitzer/IRAC 4.5$mu$m band which, for a redshift above $z=9.0$, likely indicates the presence of a rest-frame Balmer break and a stellar component that formed earlier than a redshift $zsimeq10$. In addition to constraining the earlier star formation activity on the basis of fits using stellar population models with BAGPIPES, we have undertaken the necessary, but challenging, follow-up spectroscopy for each candidate using various combinations of Keck/MOSFIRE, VLT/X-shooter, Gemini/FLAMINGOS2 and ALMA. Based on either Lyman-$alpha$ or [OIII] 88 $mu$m emission, we determine a convincing redshift of $z$=8.78 for GN-z-10-3 and a likely redshift of $z$=9.28 for the lensed galaxy MACS0416-JD. For GN-z9-1, we conclude the case remains promising for a source beyond $zsimeq$9. Together with earlier spectroscopic data for MACS1149-JD1, our analysis of this enlarged sample provides further support for a cosmic star formation history extending beyond redshifts $zsimeq$10. We use our best-fit stellar population models to reconstruct the past rest-frame UV luminosities of our sources and discuss the implications for tracing earlier progenitors of such systems with the James Webb Space Telescope.
The majority of galaxies with current star-formation rates (SFRs), SFRo >= 10^-3 Msun/yr, in the Local Cosmological Volume where observations should be reliable, have the property that their observed SFRo is larger than their average star formation rate. This is in tension with the evolution of galaxies described by delayed-tau models, according to which the opposite would be expected. The tension is apparent in that local galaxies imply the star formation timescale tau approx 6.7 Gyr, much longer than the 3.5-4.5 Gyr obtained using an empirically determined main sequence at several redshifts. Using models where the SFR is a power law in time of the form propto (t - t1)^eta for t1 = 1.8 Gyr (with no stars forming prior to t1) implies that eta = 0.18 +- 0.03. This suggested near-constancy of a galaxys SFR over time raises non-trivial problems for the evolution and formation time of galaxies, but is broadly consistent with the observed decreasing main sequence with increasing age of the Universe.