This thesis investigates the morphology of early-type galaxies in two rich clusters using 2D surface photometry. In particular, the amount of light in the 'disk' component is focussed upon, as the presence of a disk is the main morphological criterion in distinguishing between the traditional 'elliptical' and 'S0' classes. Extensive and photometric E-band CCD observations of continuous areas of the Coma and Abell 1367 clusters were obtained at the 2.5 m Isaac Newton telescope, La Palma during March 1994. A subset of this large data-set has been used in this study, comprising a magnitude-limited (to R = 15.6) sample of 153 galaxies in the two clusters. Surface photometry measurements, including surface brightness profiles and isophotal shapes, have been made for the sample. Atmospheric seeing is a major problem when measuring light profiles at the distance of Coma from ground-based telescopes. Typical seeing at La Palma (FWHM~1.2") is a significant fraction of the effective radius of many Coma/Abell 1367 galaxies (r(_e)~3" for small ellipticals). An iterative algorithm was developed to deconvolve the effects of seeing from surface brightness profiles. The result of the algorithm is to extend the range of useful surface photometry inwards to within 2 times the FWHM. In order to parametrise the surface brightness profiles and discriminate between different profile-types, further software was developed to fit one- and two-component model profiles to the seeing-corrected data. The following parameters were measured and tabulated for each of the 153 galaxies: total magnitude M(_t); half-light parameters r (_1/2) and (μ)(_1/2); SB at half-light radius μ(r(_1/2)); photometric diameter D(_19.23) (equivalent to D(_n)); ellipticity at R = 21.5 isophote ϵ(_21.5); averaged isophote high-order terms (C(_3)), (S(_3)), (C(_4)) and (S(_4)); effective radii and surface brightnesses of 5 single power-law r(^1)(_n) models, r"e and (^)"^ (n = 1,2,3,4,5); best-fitting power-law index n; bulge effective radii and surface brightnesses from the two-component fit and (/^)\; disk effective parameters r'^e and {nY^] and disk-to-bulge luminosity ratio DjB. The measured parameters have been used to investigate various aspects of early-type galaxy morphology. The conclusions are outlined below. Firstly, a two-component r? plus exponential model is a better fit to most galaxies than a single component law fit. Secondly, the traditional division of early-type galaxies into 'elliptical' and 'SO' classes is severely biased by the viewing angle. In fact, it appears that early-type galaxies comprise a population of objects with smoothly varying bulge-to-disk ratio - although a few ellipticals (less than 13%) do not appear to have a exponential component. Finally, there is a general correlation (with much scatter) between the size and the profile shapes of early-type galaxies. The interpretation is that smaller galaxies are more disk-dominated than larger galaxies, which can be linked to the merging process in rich clusters.

Powerful forces are at work in giant Ly α nebulae, a rare and mysterious population in the high redshift universe. Much like the spatially extended emission line halos around high redshift radio galaxies . but without the strong radio emission . Ly α nebulae (or Ly α 'blobs') boast copious Ly α emission (10⁴⁴ erg s⁻¹), large sizes (∼100 kpc), complex gas morphologies, and the company of numerous compact, star-forming galaxies, and may offer a window into dramatic episodes of massive galaxy formation. The small sample sizes and complex inner workings of Ly α nebulae have limited progress on understanding the their space density, environments, and physical conditions. This thesis strives to answer fundamental questions about Ly α nebulae and pave the way for understanding their role in the build up of massive galaxy systems. To address the frequency of collapse of these massive structures, we carried out the largest systematic Ly α nebula survey to date and measured the Ly α nebula space density. As an unbiased test of the environment of Ly α nebulae, we studied the surroundings of a Ly α nebula and confirmed that Ly α nebulae reside preferentially in overdense regions. To disentangle the sources of ionization, we took a census of all the compact ionization sources within a large Ly α nebula using high resolution imaging. Finally, we used photoionization modeling to put constraints on the physical conditions, the metallicity, and the sources of ionization within Ly α nebulae. Future work will be able to build on this thesis by expanding the systematic search for Ly α nebulae to other existing deep broad-band datasets, mapping the three-dimensional overdense structures in which Ly α nebulae live out to ≥ 50 (comoving) Mpc scales, and disentangling multiple sources of ionization within a larger sample of individual systems using deep optical and near-infrared spectroscopy and detailed photoionization modeling.

Galaxy formation is one of the central problems of Physical Cosmology. Neutral hydrogen plays an important role, linking the collapse of cooling gas into haloes with the formation of stars. Lyman alpha, hydrogen's strongest spectral line, can directly probe neutral hydrogen in the high redshift Universe. Lyα can be observed in absorption in Damped Lyman Alpha systems (DLAs): high Hi column density regions that dominate the neutral gas content of the Universe between z ~ 0-5. Lyα in emission is an important signature of early, starforming galaxies. Both populations, however, present significant theoretical challenges. As part of my thesis, I have developed a Monte Carlo Lyα radiative transfer code to investigate models of early galaxies. Rauch et al. (2008) performed an ultra-deep spectroscopic survey and discovered a new population of very faint, spatially extended Lyα emitters, which they claimed to be the long-sought host galaxies of DLAs at z ~ 3. I show here that a simple analytical model, which reproduces the incidence rate and kinematics of DLAs in the context of λCDM models for structure formation, also reproduces the size distribution of the faint Lyα emitters for plausible parameters, which supports their identification as DLA host galaxies. The model suggests that galaxies in haloes with vc ~ 100-150 km s-1 account for the majority of DLA host galaxies, and that these galaxies at z ~ 3 are the building blocks of typical present-day galaxies like our Milky Way. I further use my newly developed Lyα code to perform detailed 1D radiative transfer calculations, investigating the spatial and spectral distribution of Lyα emission due to star formation at the centre of DLAs, and its dependence on the spatial and velocity structure of the gas. The modelling reproduces the observed properties of both DLAs and the faint Lyα emitters, including the velocity width and column density distribution of DLAs and the large observed spatial extent of the faint emitters. In the model, haloes hosting DLAs retain up to 20% of the cosmic baryon fraction in the form of neutral hydrogen. The scattering of Lyα photons at the observed radii, which can be as large as 50 kpc, requires the bulk velocity of the gas at the centre of the haloes to be moderate. I furthermore perform 3D Lyα radiative transfer simulations, building on numerical simulations of galaxy formation that include galactic winds and gas infall. The Lyα emission region is shown to be larger and smoother than the cross-section for damped absorption by ~ 50%, with Lyα photons scattered effectively by gas with column densities >~ 1017 cm⁻². The spectra typically show two peaks, with the relative strength of the red (blue) peak being a reflection of the relative contribution of outflow (inflow) in the velocity profile. There is considerable variation in the observed line profile and spectral intensity with viewing angle. These more realistic models support many of the simplifying assumptions of my previous models, and have the potential to probe the important role of galactic winds in protogalaxies. The main conclusion is that the faint population of Lyα emitters are indeed the long sought host population of DLAs. Ultra-faint observations of Lyα emission have exceptional potential to directly probe the spatial distribution and kinematics of neutral hydrogen in early galaxies.

I review theoretical expectations for the probable appearance of galaxies during their formation phase, placing particular emphasis on the uncertainties in these ideas. Recent models suggest that formation may occur relatively recently, but that young galaxies are less spectacular than previously supposed. They may be analogous to recently discovered high red - shift radio galaxies, and indeed they may have been observed directly in faint galaxy counts. I summarise several other lines of evidence which suggest that galaxy formation may have been a recent process. Finally I give preliminary results from a detailed analytic study of the observable properties of young galaxies in a Cold Dark Matter universe. Predictions are given for faint galaxy counts and redshift distributions, and for the galaxy luminosity function.

Cosmological hydrodynamic simulations have shown that galaxies are fed by dense, cold gas streams at high redshift. However, the presence of such gas has never been observationally confirmed. Using the Horizon- MareNostrum simulation, I examined whether cold flows are detectable with low-ionisation metal absorption lines, such as C II 1334. It is concluded that due to their low metallicity and density, it is extremely difficult to prove/disprove the presence of cold flows using the metal absorption lines. Revisiting the acquisition of angular momentum in disc galaxies using high resolution simulations, I found that at the time of accretion, gas and dark matter do carry a similar amount of specific angular momentum which is systematically and significantly higher (at minimum by a factor of 2) than that of the dark matter halo as a whole. Whereas cold streams directly deposit this large amount of angular momentum within a sphere of radius r~0.1 Rvir, dark matter particles easily pass through the central region, depositing their angular momentum over a much more spatially extended region. As a result, in our simulations neither the total specific angular momentum of the baryons nor its radial profile ever follows that of the virialised dark matter halo, contrary to what is typically assumed in the standard theory of disc galaxy formation. In order to better understand the formation of disc galaxies and the missing baryon problem in a LCDM universe, continuous, collective galactic winds are implemented. It is demonstrated that stellar feedback processes are able to suppress star formation by ~30% at z=3, compared to that from the run without feedback sources, but it still produces an unrealistic central peak in the rotation curve. Although inclusion of hypernovae further suppresses star formation, it is unable to quench the formation of low-angular momentum stars enough to remove the peaked rotation curves at high redshift. Finally, feedback from active galactic nuclei turns out to be effective at suppressing star formation in massive galaxies at 1<z<2, reproducing their observed number densities in the redshift range. However, further suppression of residual star formation is required to form quiescent galaxies at z=2.

The majority of galaxy formation studies which explore beyond local redshifts do not typically probe down to the dwarf galaxy stellar mass range of ∼ 10^9 Msun . Thus trends in the observed evolution or characteristics of galaxy formation at a particular epoch are based upon populations of massive galaxies. However the currently favored Λ-Cold Dark Matter (Λ-CDM) theory is based upon hierarchical clustering and merging of lower mass systems, which proceed to make the higher mass, complex morphology of galaxies we observe. Thus it is clear that within the dwarf galaxy mass regime there should be a significant phase of galaxy formation and evolution. This work aims to uncover the influence of local environment on the formation and evolution of dwarf and massive galaxies beyond local redshift, probing down to a mass range lower than that which has been explored by previous studies. A previously successful study titled the Redshift One LDSS-3 Emission line Survey (ROLES), released results for a redshift of z ∼ 1, which compared the [OII] luminosity and galaxy stellar mass functions ([OII] LF and GSMF respectively), star formation rate density (SFRD), and specific star formation rate (sSFR) relations, with a local SDSS dataset. This led to the expansion of the study to lower redshift (this work) which explored low stellar mass galaxies at a redshift of z ∼ 0.75. This follow-up study referred to as ROLES75 (z ∼ 0.75) targeted the same two deep fields explored by the z ∼ 1 study (GOODS-South, MS1054-03 FIRES), which have extensive public photometry. Low mass targets were selected for study by their K-magnitudes (22.5 < KAB < 24) leading to a dwarf mass range of 8.5 < Log(M∗/Msun) < 9.5, and which were most likely to be within our redshift range (0.62 < z < 0.885). Follow-up multi-object spectroscopy targeted the [OII]λ3727A emission line star formation tracer in these targets allowing us to identify and obtain secure spectroscopic redshifts, SED-fit stellar masses and observed [OII] luminosity calibrated star formation rates down to limits of Log(M∗/Msun) ∼ 8.85 and SFR ∼ 0.1 Msun/yr . Science results presented here are similar to those published by the ROLES z ∼ 1 study, however we also studied the influence of the high versus low density environment in which the galaxy populations reside. This study confirmed that while the [OII] luminosity was higher in earlier times, environment does not influence galaxy formation at z ∼ 0.75. The faint-end slope of the [OII] LF, α ∼ 1.25 measured here, is also observed to become increasingly more steep with increasing redshift. The [OII] luminous GSMF is observed to not have significantly evolved since z ∼ 2.75, confirming the result of the previous ROLES work. However the impact of environment on the GSMF is apparent in the high mass end where the imprint of structure from the CDFS field enhances the stellar mass function above the field population. There is also weak evidence of a bi-modal [OII] luminous GSMF indicated by an ‘upturn’ near ∼ 10^9 Msun in the low density field population. The SFRD at z ∼ 0.75 does not confirm the picture presented by the ROLES z ∼ 1 study where a constant scale factor was applicable to the local SDSS SFRD to obtain the z ∼ 1 SFRD. The SFRD in the high mass end at z ∼ 0.75 is lower than would be expected based upon a constant scale factor, while the low stellar mass end exhibits some consistency with this picture. In the high density environment, this dominant SFRD (over the low density field population) is driven by the high density [OII] luminous GSMF in the high stellar mass end, rather than through an enhancement of the SFR. The normalization of the sSFR − M∗ relation at z ∼ 0.75 is found to lie between those corresponding to z ∼ 1 and present day. There is a subtle ‘upturn’ in the sSFR − M∗ relation confirming this observation which was also present in the ROLES z ∼ 1 study but not present in the local SDSS sSFR − M∗ relation. The sSFR of active galaxies does not depend upon the local density in which they are forming, confirming the same conclusion based upon the [OII] LF. However, there is redshift evolution of the sSFR − M∗ relation with respect to local density. The high density sSFR − M∗ relation for star forming galaxies was dominant over its low density counterpart at early times, with the opposite the case at present day. There is suggestion of the crossover or rollover transition occurring at z ∼ 0.75.

We present a series of numerical galaxy formation studies which apply new numerical methods to produce increasingly realistic galaxy formation models. We first investigate the metallicity evolution of a large set of idealized hydrodynamical galaxy merger simulations of colliding galaxies. We find that inflows of metal--poor interstellar gas triggered by galaxy tidal interactions can account for the systematically lower central oxygen abundances observed in local interacting galaxies. We show the central metallicity evolution during merger events is determined by a competition between the inflow of low--metallicity gas and enrichment from star formation. We find a time-averaged depression in the galactic nuclear metallicity of ~0.07 dex for gas--poor disk--disk interactions, which explains the observed close pair mass-metallicity and separation-metallicity relationships. / Astronomy

This thesis presents multi-frequency data, galaxy identifications, estimated redshifts, and derived physical properties for the sub-millimetre source sample produced by the SCUBA HAlf Degree Extragalactic Survey (SHADES). SHADES is the largest, complete, sub-millimetre survey undertaken to date, and the aim of this work is to exploit this survey to study the evolution of sub-mm selected galaxies at high redshift, explore their possible connectionwith localmassive ellipticals, and to test current models of galaxy formation. The SHADES sample was selected from 850 micron images made with the submillimetre camera SCUBA at the James Clerk Maxwell Telescope. These submillimetre maps cover a total area of 720 arcmin2 split between two well-studied extra-galactic survey fields, the Lockman Hole East and the Subaru/XMMNewton Deep Field (SXDF). The resulting sample of 120 sub-millimetre sources is the focus of this thesis. Here the wealth of information provided by deep radio, optical, near-infrared and mid-infrared imaging of the two SHADES fields is exploited to complete the identification of the SHADES sample, and then to derive a robust redshift estimate for every sub-millimetre source. Where possible this is achieved from the optical+ infrared photometry using a new two-component redshift estimation code developed specifically to deal with starbursting galaxies with potentially highly stochastic star-formation histories. The effectiveness of this code is demonstrated via comparison with the small subset of SHADES source which possess robust spectroscopic redshifts. For those galaxies which are too faint for effective redshift constraints to be provided by the existing optical+infrared photometry, the information on the long-wavelength spectral energy distribution provided by the radio+submm photometry is utilised to provide cruder constraints or limits on redshift. The result is the first complete and unbiased estimate of the redshift distribution of the bright extragalactic sub-millimetre galaxy population. It is found that the brightest sub-mm sources are confined to the redshift range 2 < z < 4, while more moderate luminosity sources span the full range of redshift out to z ∼ 5. The fits to themulti-frequency photometry provided by the redshift estimation technique are also used to derive estimates of the stellar mass, and star-formation history of each SHADES galaxy. The average derived stellar mass is ∼ 3 × 1011 M⊙ and it is found that the violent starburst powering the sub-millimetre emission typically contributes less than 10% of the stellar mass of the galaxy which has been assembled prior to the “current” starburst event. The distributions of redshift, stellar mass, and star-burst ages are compared with the predictions of a range of galaxy models, including the suite of models originally used to motivate the SHADES survey in van Kampen et al. (2005), and themost recent incarnation of the Durhamsemi-analytic galaxy formationmodels described by Swinbank et al. (2008). It is found that the redshift distribution and sub-mmflux versus redshift for bright sub-mmgalaxies can be reproduced best by one of the van Kampen models, which is based on semi-analytic modelling with a Chabrier IMF. We can rule out the non-semi-analytic prediction models and the Durham semi-analytic model with a top-heavy IMF. However the stellar masses are systematically underpredicted by all of the models. Either the stellar masses derived from the SHADES data have been systematically over-estimated, or the models need to be modified (perhaps by the inclusion of AGN feedback) to allow larger galaxy masses to assembled prior to z ∼ 2. Finally, it is demonstrated that themass in place prior to the observed starburst cannot have been produced by an analogous super-burst at higher redshift, but rather requires to have been assembledmore gradually over a timescale of ∼ 1−2 Gyr. It is thus concluded thatmassive galaxies undergo theirmost violent phase of star formation at redshifts 2 < z < 4, but that the enormous starbursts which lead to detection in current sub-millimetre surveys can only take place in the potential well provided by an already massive galaxy. This supports a scenario in which bright sub-millimetre galaxies are indeed the progenitors of the massive elliptical galaxies found in the local Universe.

The majority of galaxy formation studies which explore beyond local redshifts do not typically probe down to the dwarf galaxy stellar mass range of ∼ 10^9 Msun . Thus trends in the observed evolution or characteristics of galaxy formation at a particular epoch are based upon populations of massive galaxies. However the currently favored Λ-Cold Dark Matter (Λ-CDM) theory is based upon hierarchical clustering and merging of lower mass systems, which proceed to make the higher mass, complex morphology of galaxies we observe. Thus it is clear that within the dwarf galaxy mass regime there should be a significant phase of galaxy formation and evolution. This work aims to uncover the influence of local environment on the formation and evolution of dwarf and massive galaxies beyond local redshift, probing down to a mass range lower than that which has been explored by previous studies. A previously successful study titled the Redshift One LDSS-3 Emission line Survey (ROLES), released results for a redshift of z ∼ 1, which compared the [OII] luminosity and galaxy stellar mass functions ([OII] LF and GSMF respectively), star formation rate density (SFRD), and specific star formation rate (sSFR) relations, with a local SDSS dataset. This led to the expansion of the study to lower redshift (this work) which explored low stellar mass galaxies at a redshift of z ∼ 0.75. This follow-up study referred to as ROLES75 (z ∼ 0.75) targeted the same two deep fields explored by the z ∼ 1 study (GOODS-South, MS1054-03 FIRES), which have extensive public photometry. Low mass targets were selected for study by their K-magnitudes (22.5 < KAB < 24) leading to a dwarf mass range of 8.5 < Log(M∗/Msun) < 9.5, and which were most likely to be within our redshift range (0.62 < z < 0.885). Follow-up multi-object spectroscopy targeted the [OII]λ3727A emission line star formation tracer in these targets allowing us to identify and obtain secure spectroscopic redshifts, SED-fit stellar masses and observed [OII] luminosity calibrated star formation rates down to limits of Log(M∗/Msun) ∼ 8.85 and SFR ∼ 0.1 Msun/yr . Science results presented here are similar to those published by the ROLES z ∼ 1 study, however we also studied the influence of the high versus low density environment in which the galaxy populations reside. This study confirmed that while the [OII] luminosity was higher in earlier times, environment does not influence galaxy formation at z ∼ 0.75. The faint-end slope of the [OII] LF, α ∼ 1.25 measured here, is also observed to become increasingly more steep with increasing redshift. The [OII] luminous GSMF is observed to not have significantly evolved since z ∼ 2.75, confirming the result of the previous ROLES work. However the impact of environment on the GSMF is apparent in the high mass end where the imprint of structure from the CDFS field enhances the stellar mass function above the field population. There is also weak evidence of a bi-modal [OII] luminous GSMF indicated by an ‘upturn’ near ∼ 10^9 Msun in the low density field population. The SFRD at z ∼ 0.75 does not confirm the picture presented by the ROLES z ∼ 1 study where a constant scale factor was applicable to the local SDSS SFRD to obtain the z ∼ 1 SFRD. The SFRD in the high mass end at z ∼ 0.75 is lower than would be expected based upon a constant scale factor, while the low stellar mass end exhibits some consistency with this picture. In the high density environment, this dominant SFRD (over the low density field population) is driven by the high density [OII] luminous GSMF in the high stellar mass end, rather than through an enhancement of the SFR. The normalization of the sSFR − M∗ relation at z ∼ 0.75 is found to lie between those corresponding to z ∼ 1 and present day. There is a subtle ‘upturn’ in the sSFR − M∗ relation confirming this observation which was also present in the ROLES z ∼ 1 study but not present in the local SDSS sSFR − M∗ relation. The sSFR of active galaxies does not depend upon the local density in which they are forming, confirming the same conclusion based upon the [OII] LF. However, there is redshift evolution of the sSFR − M∗ relation with respect to local density. The high density sSFR − M∗ relation for star forming galaxies was dominant over its low density counterpart at early times, with the opposite the case at present day. There is suggestion of the crossover or rollover transition occurring at z ∼ 0.75.

abstract: Numerical simulations are very helpful in understanding the physics of the formation of structure and galaxies. However, it is sometimes difficult to interpret model data with respect to observations, partly due to the difficulties and background noise inherent to observation. The goal, here, is to attempt to bridge this gap between simulation and observation by rendering the model output in image format which is then processed by tools commonly used in observational astronomy. Images are synthesized in various filters by folding the output of cosmological simulations of gasdynamics with star-formation and dark matter with the Bruzual- Charlot stellar population synthesis models. A variation of the Virgo-Gadget numerical simulation code is used with the hybrid gas and stellar formation models of Springel and Hernquist (2003). Outputs taken at various redshifts are stacked to create a synthetic view of the simulated star clusters. Source Extractor (SExtractor) is used to find groupings of stellar populations which are considered as galaxies or galaxy building blocks and photometry used to estimate the rest frame luminosities and distribution functions. With further refinements, this is expected to provide support for missions such as JWST, as well as to probe what additional physics are needed to model the data. The results show good agreement in many respects with observed properties of the galaxy luminosity function (LF) over a wide range of high redshifts. In particular, the slope (alpha) when fitted to the standard Schechter function shows excellent agreement both in value and evolution with redshift, when compared with observation. Discrepancies of other properties with observation are seen to be a result of limitations of the simulation and additional feedback mechanisms which are needed. / Dissertation/Thesis / Ph.D. Physics 2012