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.

This thesis investigates the dynamical evolution of systems orbiting within deeper potentials. Initially we use a simple satellite-halo interaction to study the dynamical processes that act on orbiting systems and we compare these results to analytical theory. Deep images of the Centaurus cluster reveal a spectacular arc of diffuse light that stretches for over 100 kpc and yet is just a few kpc wide. We use numerical simulations to show that this feature can be produced by the tidal debris of a spiral galaxy that has been disrupted by the potential of one of the central cD galaxies of the cluster. The evolution of sub-halos is then studied in a cosmological context using high resolution N-body simulations of galactic mass halos that form in a cold dark matter (CDM) simulation. CDM halos form via a complex series of mergers, accretion events and violent relaxation. Halos are non-spherical, have steep singular density profiles and contain many thousands of surviving dark matter substructure clumps. This will lead to several unique signatures for experiments that aim to detect dark matter either indirectly, through particle annihilation, or directly in a laboratory. For the first time it is possible to construct maps of the gamma-ray sky that result from the annihilation of dark matter particles within simulated dark matter halo distributions.

We present two self-consistent non-parametric models of the local cosmic velocity field based on the density distribution in the PSCz redshift survey of IRAS galaxies. Two independent methods have been applied, both based on the assumptions of gravitational instability and linear biasing. They give remarkably similar results, with no evidence of systematic differences and an r.m.s discrepancy of only ~ 70 kms(^-1) in each Cartesian velocity component. These uncertainties are consistent with a detailed independent error analysis carried out on mock PSCz catalogues constructed from TV-body simulations. The denser sampling provided by the PSCz survey compared to previous IRAS galaxy surveys allows us to reconstruct the velocity field out to larger distances. The most striking feature of the model velocity field is a coherent large-scale streaming motion along a basehne connecting Perseus-Pisces, the Local Supercluster, the Great Attractor, and the Shapley Concentration. We find no evidence for back-infall onto the Great Attractor. Instead, material behind and around the Great Attractor is inferred to be streaming towards the Shapley Concentration, aided by the expansion of two large neighbouring un- derdense regions. The PSCi model velocities compare well with those predicted from the 1.2-Jy redshift survey of IRAS galaxies and, perhaps surprisingly, with those predicted from the distribution of Abell/ACO clusters, out to 140 h(^-1)Mpc. Comparison of the real-space density fields (or, alternatively, the peculiar velocity fields) inferred from the PSCz and cluster catalogues gives a relative (linear) bias parameter between clusters and IRAS galaxies of b(_c) = 4.4 ± 0.6. In addition, we compare the cumulative bulk flows predicted from the PSCz gravity field with those measured from the MarkIII and SFI catalogues of peculiar velocities. A conservative estimate of β = Ω(_0)(^0.6)/b, where b is the bias parameter for IRAS galaxies, gives β= 0.76 ± 0.13 (1-σ), in agreement with other recent determinations. Finally, we perform a detailed comparison of the IRAS PSCz and 1.2-Jy spherical harmonic coefficients of the density and velocity fields in redshift space. Both the monopole terms of the density and velocity fields predicted from the surveys show some inconsistencies. The mismatch in the velocity monopole terms is resolved by masking the 1.2-Jy survey with the PSCz mask and using the galaxies within the PSCz survey for fluxes larger than 1.2 Jy. Davis, Nusser and Willick (1996) have found a discrepancy between the IRAS 1.2-Jy survey gravity field and the MarkIII peculiar velocity field. We conclude that the use of the deeper IRAS PSCz catalogue cannot alone resolve this mismatch.

The origin, nature and distribution of dark matter in the universe form some of the biggest questions in modern astrophysics. Dark matter is distributed on a wide range of scales in the universe. This thesis concentrates on galactic scales, attempting to lower the veil and probe the structure of dark matter in galaxies. (For complete abstract open document)

In this thesis I investigate galaxy evolution by measuring the luminosity functions of galaxies across a wide range in redshift. I measure the abundances of high redshift galaxies in deep HST imaging of the GOODS-North field from the CANDELS survey. I follow this evolution to lower redshift by measuring the luminosity functions of galaxies in ground-based imaging of the XMM-LSS field from the VIDEO survey with optical data from the CFHTLS. First, at high-redshift, I identify 22 candidate z &asymp; 7 and 6 candidate z &asymp; 8 − 9 galaxies. By comparing the number of candidate galaxies with those found in the GOODS-South field, I determine that cosmic variance is not the dominant source of uncertainty on the number counts. I constrain the Schechter parameters for the UV luminosity function at z &asymp; 7 and z &asymp; 8 − 9, finding evidence for evolution in the number density of high redshift galaxies. Next, I present the K_s-band luminosity functions in the 1 degree² and 4.5 degree² overlaps between the VIDEO-XMM field and the CFHTLS-D1 and CFHTLS-W1 fields. I measure the luminosity functions with the 1/V_max method over the range 0.2 &LT; z &LT; 3 in VIDEO-CFHTLS-D1, and over the range 0.2 &LT; z &LT; 1.5 in the shallower VIDEO-CFHTLS-W1 field. I find the evolution of these luminosity functions is best described by luminosity dependent density evolution, where the characteristic magnitude has dimmed at a constant rate since z = 3, while the density has increased since z = 3, first rapidly from z = 3 to z &asymp; 1.5 and then more slowly from z &asymp; 1.5 to z = 0.2. I measure a significant upturn at the faint end of the luminosity function at low redshift. Finally, I compare the VIDEO-CFHTLS-D1 and VIDEO-CFHTLS-W1 luminosity functions with predicted K-band luminosity functions from the Horizon-AGN simulation. I find both an over-prediction in the numbers of faint galaxies and an under-prediction in the numbers of bright galaxies, implying that the feedback from supernovae is insufficient while the feedback from AGN is over-sufficient.

Supermassive black holes (BHs) and their host-galaxies are thought to evolve in tandem, with the energy output from the rapidly-accreting BH regulating star formation and the growth of the BH itself. The goal of better understanding this process has led to much work focussing on the properties of quasars at high redshifts, $z\gtrsim 2$, when cosmic star formation and BH accretion both peaked. At these redshifts, however, ground-based statistical studies of the quasar population generally have no access to the rest-frame optical spectral region, which is needed to measure H$\beta$-based BH masses and narrow line region outflow properties. The cornerstone of this thesis has been a new near-infrared spectroscopic catalogue providing rest-frame optical data on 434 luminous quasars at redshifts $1.5 \lesssim z \lesssim 4$. At high redshift, $z \gtrsim 2$, quasar BH masses are derived using the velocity-width of the CIV broad emission-line, based on the assumption that the observed velocity-widths arise from virial-induced motions. However, CIV exhibits significant asymmetric structure which suggests that the associated gas is not tracing virial motions. By combining near-infrared spectroscopic data (covering the hydrogen Balmer lines) with optical spectroscopy from SDSS (covering CIV), we have quantified the bias in CIV BH masses as a function of the CIV blueshift. CIV BH masses are shown to be over-estimated by almost an order of magnitude at the most extreme blueshifts. Using the monotonically increasing relationship between the CIV blueshift and the mass ratio BH(CIV)/BH(H$\alpha$) we derive an empirical correction to all CIV BH-masses. The correction depends only on the CIV line properties and therefore enables the derivation of un-biased virial BH mass estimates for the majority of high-luminosity, high-redshift, spectroscopically confirmed quasars. Quasars driving powerful outflows over galactic scales is a central tenet of galaxy evolution models involving 'quasar feedback' and significant resources have been devoted to searching for observational evidence of this phenomenon. We have used [OIII] emission to probe ionised gas extended over kilo-parsec scales in luminous $z\gtrsim2$ quasars. Broad [OIII] velocity-widths and asymmetric structure indicate that strong outflows are prevalent in this population. We estimate the kinetic power of the outflows to be up to a few percent of the quasar bolometric luminosity, which is similar to the efficiencies required in recent quasar-feedback models. [OIII] emission is very weak in quasars with large CIV blueshifts, suggesting that quasar-driven winds are capable of sweeping away gas extended over kilo-parsec scales in the host galaxies. Using data from a number of recent wide-field photometric surveys, we have built a parametric SED model that is able to reproduce the median optical to infrared colours of tens of thousands of AGN at redshifts $1 < z < 3$. In individual objects, we find significant variation in the near-infrared SED, which is dominated by emission from hot dust. We find that the hot dust abundance is strongly correlated with the strength of outflows in the quasar broad line region, suggesting that the hot dust may be in a wind emerging from the outer edges of the accretion disc.

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

abstract: I combine, compare, and contrast the results from two different numerical techniques (grid vs. particle methods) studying multi-scale processes in galaxy and structure formation. I produce a method for recreating identical initial conditions for one method from those of the other, and explore methodologies necessary for making these two methods as consistent as possible. With this, I first study the impact of streaming velocities of baryons with respect to dark matter, present at the epoch of reionization, on the ability for small halos to accrete gas at high redshift. With the inclusion of this stream velocity, I find the central density profile of halos is reduced, overall gas condensation is delayed, and infer a delay in the inevitable creation of stars. I then combine the two numerical methods to study starburst outflows as they interact with satellite halos. This process leads to shocks catalyzing the formation of molecular coolants that lead to bursts in star formation, a process that is better captured in grid methods. The resultant clumps of stars are removed from their initial dark matter halo, resemble precursors to modern-day globular clusters, and their formation may be observable with upcoming telescopes. Finally, I perform two simulation suites, comparing each numerical method's ability to model the impact of energetic feedback from accreting black holes at the core of giant clusters. With these comparisons I show that black hole feedback can maintain a hot diffuse medium while limiting the amount of gas that can condense into the interstellar medium, reducing the central star formation by up to an order of magnitude. / Dissertation/Thesis / Doctoral Dissertation Astrophysics 2014

I study mid-infrared and star formation properties of AGN samples using infrared observations, and star-forming galaxies using radio observations in order to investigate the link between star formation, AGN activity and radio luminosity. I present the results of these investigations in this thesis. I carried out an analysis of four complete samples of radio-loud AGN (3CRR, 2Jy, 6CE and 7CE) using near- and mid-IR data taken by the Wide-Field Infrared Survey Explorer (WISE). The combined sample consists of quasars and radio galaxies, and covers a redshift range 0:003 < z < 3:395. The dichotomy in the mid-IR properties of low- and high-excitation radio galaxies (LERGs - HERGs) is analysed using large complete samples. The results show that a division in the accretion modes of powerful LERGs and HERGs clearly stands out in the mid-IR radio plane. Evaluation of the positions of the sample objects in WISE colour-colour diagrams shows that widely used WISE colour cuts are not completely reliable in selecting AGN. I examined the link between AGN activity and star formation by constructing matched samples of local (0 < z < 0:6) radio-loud and radio-quiet AGN in the Herschel-ATLAS fields. AGN accretion and jet powers in these active galaxies are traced by [OIII] emission-line and radio luminosity, respectively. Star formation properties were derived using Herschel 250-_m and stellar mass measurements are taken from the SDSS-MPA-JHU catalogue. The stacking analyses show that star formation rates (SFRs) and specific star formation rate (SSFRs) of both radio-loud and radio-quiet AGN increase with increasing AGN power but that radio-loud AGN tend to have lower SFR. Additionally, radio-quiet AGN are found to have approximately an order of magnitude higher SSFRs than radio-loud AGN for a given level of AGN power. The difference between the star formation properties of radio-loud and -quiet AGN is also seen in samples matched in stellar mass. I also investigated the relationship between SFR and low-frequency radio luminosity observed in star-forming galaxies. I used a sample of star-forming galaxies in the 19 local Universe selected from the SDSS-MPA-JHU catalogue. LOFAR observations of the Herschel-ATLAS North Galactic Pole field (NGP) were carried out as part of the LOFAR surveys Key Science Project at an effective frequency of 150 MHz, which provided low-frequency radio luminosity of sample galaxies. SFRs of galaxies in the sample were derived using MAGPHYS spectral energy distribution (SED) fitting. The results of this study show that the slope of L150/SFR is less than unity and not universal for all star-forming galaxies (SFGs) in the local Universe (0 < z < 0:3). The slope of the L150/SFR relation is also found to be steeper than the L1:4/SFR relation, probably due to the contribution from thermal radio emission at 1.4 GHz. If the L150=SFR relation for strongly star-forming objects is explained naively by electron calorimetry, I conclude that low luminosity sources are not ideal calorimeters and differ from the main locus of SFGs at low redshifts. The different gradients we obtain for the far- IR/radio correlation using samples selected at different frequencies reveal the selection effects on relations derived in this thesis.

Philosophiae Doctor - PhD / In this thesis, we present a multi-wavelength analysis of star-forming galaxies to shed new light on the evolution of the far-IR-radio relations in intermediate redshift (0.3 < z < 0.6) galaxy clusters and galaxy groups. The far-infrared (far-IR) emission from galaxies is dominated by thermal dust emission. The radio emission at 1.4 GHz is predominantly produced by non-thermal synchrotron radiation. The underlying mechanisms, which drive the far-IR-radio correlation, are believed to arise from massive star formation. A number of studies have investigated the relationship as a function of redshift in the field and have found no evolution out to at least z _ 2, however few works have been done in galaxy clusters. In nearby clusters, the median logarithmic ratio of the far-IR to radio luminosity is qFIR = 2.07_0.74, which is lower than the value found in the field, and there is an indication of an enhancement of radio emission relative to the far-IR emission. Understanding the properties of the far-IR-radio correlation in a sample of distant and massive cluster and groups plays an important role in understanding the physical processes in these systems. We have derived total infrared luminosities for a sample of cluster, group, and field galaxies through an empirical relation based on Spitzer MIPS 24 _m photometry. The radio flux densities were measured from deep Very Large Array 1.4 GHz radio continuum observations. We have studied the properties of the far-IR-radio correlation of galaxies at intermediate redshift clusters by comparing the relationship of these galaxies to that of low redshift clusters. We have also examined the properties of the galaxies showing radio excess to determine the extent that galaxy type or environment may explain the radio excess in galaxy clusters. We find that the ratio of far-IR to radio luminosity for galaxies in an intermediate redshift cluster to be qFIR = 1.72_0.63. This value is comparable to that measured in low redshift clusters. A higher fraction of galaxies in clusters show an excess in their radio fluxes when compared to low redshift clusters, and corroborates previous evidence of a cluster enhancement of radio excess sources at this earlier epoch as well. We have also investigated the properties of the far-IR-radio correlation for a sample of galaxy groups in the COSMOS field. We find a lower percentage of radio-excess sources in groups as compared to clusters. This provides preliminary evidence that the number of radioexcess sources may depend on galaxy environment. We also find that a larger fraction of radio-excess sources in clusters are red sequence galaxies.