In this work I present an improved model for following the formation of Pop III stars, their effects on early galaxy evolution, and how we might search for them. I make use of a new subgrid model of turbulent mixing to accurately follow the time scales required to mix supernova (SN) ejecta -- enriched with heavy elements -- into the pristine gas. I implement this model within a large-scale cosmological simulation and follow the fraction of gas with metallicity below a critical value marking the boundary between Pop III and metal enriched Population II (Pop II) star formation. I demonstrate that accounting for subgrid mixing results in a Pop III stars formation rate that is 2-3 times higher than standard models with the same physical resolution.
I also implement and track a new "Primordial metals" (PM) scalar that tracks the metals generated by Pop III SNe. These metals are taken up by second generation stars and likely result in a subclass of carbon-enhanced, metal-poor (CEMP) stars. By tracking both regular metals and PM, I can model, in post-processing, the elemental abundances of simulation stars. I find good agreement between observations of CEMP-no Milky Way halo stars and second generation stars within the simulation when assuming the first stars had a typical mass of 60 M☉, providing clues as to the Pop III initial mass function.
star formation (Dahlem et al. 2006) by comparing maps of 120-240 MHz synchrotron emission and hydrogen alpha (Hα) emission of the tidally-interacting, edge-on, barred spiral galaxy UGC 9665. Synchrotron emission traces magnetic field strength to a rough first order, while Hα emission traces recent massive star formation. UGC 9665 was selected because it was included in the LOw Frequency ARray (LOFAR) TwoMetre Sky Survey (LoTSS; Shimwell et al. (2017)) as well as the Calar Alto Legacy Integral Field Area Survey (CALIFA; Sanchez et al. (2012)). I generated vertical intensity profiles at several distances along the disk from the galactic center for synchrotron emission and Hα in order to measure how the intensity of each falls off with distance from the midplane. In addition to correlating the vertical profiles to see if there is a relationship between star formation and magnetic field strength, I fit the LOFAR vertical profiles to characteristic Gaussian and exponential functions given by Dumke et al. (1995). Fitting these equations have been shown to be good indicators of the main mode of cosmic ray transport, whether it is advection (exponential fit) or diffusion (Gaussian fit) (Heesen et al. 2016). Cosmic rays originate from supernova,
and core collapse supernovae occur in star forming regions, which also produce
advective winds, so I test the correlation between star-forming regions and advective regions as predicted by the Heesen et al. (2016) method. Similar studies should be conducted on different galaxies in the future in order to further test these hypotheses and how well LOFAR and CALIFA complement each other, which will be made possible by the full release of the LOFAR Two-Metre Sky Survey (LoTSS) (Shimwell et al. 2017).
The reionization of the Universe is thought to have completed by redshift z~5.5. To probe this era, galaxy observations in the Subaru Deep Field (SDF) have identified more than 100 galaxies at z~6, many spectroscopically confirmed through follow-up observations. Using available optical/IR data, we model with CIGALE the spectral energy distributions (SEDs) of 43 SDF galaxies, including newly acquired data from the UKIRT WFCAM K-band for seven previously studied objects. In particular, modeling deep IR photometry is sensitive to the galaxy's Lyman continuum (LyC) escape fraction (fesc). We find the median implied fesc value as ~0.4+/-0.1 (mean error). Significant uncertainties in data and fitting result in a large range of fesc for individual objects, but analysis suggests that fesc is likely high enough for galaxies to finish reionization by z~6. More importantly, we find trends between the CIGALE UV slope b, fesc, and dust extinction E(B-V): for a given E(B-V), b appear steeper by ~0.4 than at z=0. Lower fesc values appear to be associated with bluer b and lower E(B-V), but only weakly. This suggests that LyC could have escaped through holes with sufficiently wide opening angles surrounding the ISM from outflows of supernovae and/or weak AGN to escape, but resulting in a large range of implied fesc values depending on the orientation of each galaxy. The current HST, Spitzer and ground-based photometric and model errors for the 43 galaxies are large, so IR spectroscopic observations with the James Webb Space Telescope are needed to better constrain this possibility.