Determining cosmological parameters from current observational data requires knowledge of the primordial density perturbations generated during inflation. We begin by examining a model of inflation along a flat direction of the minimal supersymmetric Standard Model (MSSM) and the power spectrum of perturbations it can produce. We consider the fine-tuning issues associated with this model and discuss a modification of the potential to include a hybrid transition that reduces the fine-tuning, without affecting the viability of the model. However, supersymmetric flat directions might play a role in other models of inflation as well. In particular, they may cause a feature in the primordial power spectrum of perturbations, unlike the scale-free spectrum assumed in the standard Lambda Cold Dark Matter (LCDM) cosmological model. We then show that in the presence of such a feature, an alternative cosmological model with a large local void and no dark energy provides a good fit to both Type Ia supernovae and the cosmic microwave background (CMB) data from the WMAP satellite. Constraints from the locally measured Hubble parameter, baryon acoustic oscillations and primordial nucleosynthesis are also satisfied. This degeneracy motivates a search for other independent observational tests of LCDM. The integrated Sachs-Wolfe (ISW) imprint of large-scale structure on the CMB is one such test. The ISW imprint of superstructures of size ~100 Mpc/h at redshift z~0.5 has been detected with >4 sigma significance, however it has been noted that the signal is much larger than expected. We revisit the calculation using linear theory predictions in a LCDM cosmology and find the theoretical prediction is inconsistent by >3 sigma with the observation. If the observed signal is indeed due to the ISW effect then huge, extremely underdense voids are far more common in the observed universe than predicted by LCDM.

In this thesis, I present the results of a new search for, and study of, luminous galaxies in the first billion years of cosmic time. This work is primarily based on a new selection of bright (L≫L*) Lyman-break galaxies (LBGs) at z ≅6 and z≅7 in the UltraVISTA first and second data releases (DR1, DR2) and the UKIDSS (UKIRT Infrared Deep Survey) UDS DR10 (Ultra Deep Survey). The UltraVISTA survey provides deep Y, J,H andKs near-imaging over 1.5 deg² of the Cosmic Evolution Survey (COSMOS) field and the UKIDSS UDS provides J,H and K band data overlapping with the Subaru XMM-Newton Deep Survey (SXDS), with both fields also containing deep optical and mid-infrared imaging essential for the clean detection of z > 5 galaxies. The fields combined provide an unprecedented 1.65 deg² of deep multiwavelength data with which to securely select LBGs using a photometric redshift fitting technique, which can additionally remove probable low-redshift galaxy interlopers and galactic dwarf stars that can contaminate ground-based samples. At z ≅7, the DR1 of the UltraVISTA survey was used to select a sample of ten high-redshift galaxy candidates, which extended to a 5σ limiting magnitude of Y + J ~ 25 (AB magnitude, 2-arcsec diameter circular aperture) over 1 deg². A stack of the four most robust objects from the sample indicated that they were massive (M*≅5 × 109M ʘ), had blue rest-frame UV slopes (β ≅−2.0±0.2) and were highly star-forming (SFR ≅25–50Mʘ yr−1) when compared to previous, fainter, samples of galaxies at z = 7. The number counts of z≅7 galaxies selected within the UltraVISTA DR1 survey was higher than that expected from extrapolations of the rest-frame UltraViolet (UV) luminosity function (LF) from fainter data, a result that was strongly confirmed with an improved search for z ≅7 galaxies using the UltraVISTA DR2 imaging and the UDS field. A total of 34 galaxies at 6.5 < z < 7.5 were found in the combined fields, which included the previously identified robust galaxies from the DR1 imaging. This expanded sample allowed the first determination of the rest-frame UV LF in the range −23.0 < MUV < −21.5 at z ≅7, and the results reveal a power-law decline to bright magnitudes in contrast to the commonly assumed exponentially declining Schechter function extrapolated from fainter data. The excess of galaxies observed at bright magnitudes cannot be accounted for by gravitational lensing or by significant contamination of the sample by Active Galactic Nuclei (AGN) . The observed LF is well described by a double power law, which at the bright end follows the form of the underlying dark matter halo mass function, suggesting that the physical mechanism that inhibits star formation activity in massive galaxies (e.g. AGN feedback or some other form of ‘mass quenching’) has yet become efficient at z ≅7. The deeper imaging data confirm that the z ≅7 LBGs show blue rest-frame UV slopes (median β = −2.0) and are massive (up to M*≅1010M ʘ). Furthermore, an analysis of the ground-based imaging shows that the majority are resolved consistent with larger sizes (r1/2 ≅1–1.5 kpc) than displayed by less massive galaxies. Finally, a new search for z ≅6 galaxies within the UltraVISTA and UDS datasets was undertaken, resulting in a sample of 266 LBGs (−22.7 < MUV < −20.5) galaxies with which to investigate the rest-frame UV LF. The potential contamination by galactic brown dwarfs was investigated quantitatively using a simple model of the Galaxy, showing that the expected contamination rate of the sample was < 3 per cent, and that the stars can be effectively removed by fitting standard stellar spectra to the observed photometry. The galaxy surface density in the UltraVISTA/COSMOS field exceeds that in the UDS/SXDS by a factor of ≅1.4, indicating strong cosmic variance between the two fields. The number counts of galaxies we find are a factor of 2 lower than predicted by the recent LF determination by Bouwens et al., and the derived rest-frame UV LF at z ≅6 revealed that an under dense UDS field can account for some of the observed differences between previous analyses. An evolution in the characteristic magnitude between z ≅5 and z≅7 of ∆M*~ 0.5 was found in contrast to other smaller area surveys, and a double power law was shown to equally well describe the LF at z = 6 as compared to the commonly assumed Schechter function. The bright-end of the LF at z ≅6 tentatively shows a steeper decline than found at z ≅7, which could indicate the onset of mass quenching of the most massive galaxies or the rise of dust obscuration. Comparison with the predictions of the latest theoretical models and simulations of galaxies reveals that most models require substantial (A1500 ~ 1.5–2) average dust extinction at the bright end to reproduce the shape of the galaxy UV LF at z ≅7.

Gravitational lensing studies the distortions of a distant galaxy’s observed size, shape or flux due to the tidal bending of photons by matter between the source and observer. Such distortions can be used to infer knowledge on the mass distribution of the intervening matter, such as the dark matter halos in which clusters of individual galaxies may reside, or on cosmology through the statistics of the matter density of large scale structure and geometrical factors. In particular, gravitational lensing has the advantage that it is insensitive to the nature of the lensing matter. However, contamination of the signal by correlations between galaxy shape or size and local environment complicate a lensing analysis. Further, measurement of traditional lensing estimators is made more difficult by limitations on observations, in the form of atmospheric distortions or optical limits of the telescope itself. As a result, there has been a large effort within the lensing community to develop methods to either reduce or remove these contaminants, motivated largely by stringent science requirements for current and forthcoming surveys such as CFHTLenS, DES, LSST, HSC, Euclid and others. With the wealth of data from these wide-field surveys, it is more important than ever to understand the full range of independent probes of cosmology at our disposal. In particular, it is desirable to understand how each probe may be used, individually and in conjunction, to maximise the information of a lensing analysis and minimise or mitigate the systematics of each. With this in mind, I investigate the use of galaxy clustering measurements using photometric redshift information, including a contribution from flux magnification, as a probe of cosmology. I present cosmological forecasts when clustering data alone are used, and when clustering is combined with a cosmic shear analysis. I consider two types of clustering analysis: firstly, clustering with only redshift auto-correlations in tomographic redshift bins; secondly, clustering using all available redshift bin correlations. Finally, I consider how inferred cosmological parameters may be biased using each analysis when flux magnification is neglected. Results are presented for a Stage–III ground-based survey, and a Stage–IV space-based survey modelled with photometric redshift errors, and values for the slope of the luminosity function inferred from CFHTLenS catalogues. I find that combining clustering information with shear gives significant improvement on cosmological parameter constraints, with the largest improvement found when all redshift bins are included in the analysis. The addition of galaxy-galaxy lensing gives further improvement, with a full combined analysis improving constraints on dark energy parameters by a factor of > 3. The presence of flux magnification in a clustering analysis does not significantly affect the precision of cosmological constraints when combined with cosmic shear and galaxy-galaxy lensing. However if magnification is neglected, inferred cosmological parameter values are biased, with biases in some cosmological parameters found to be larger than statistical errors. We find that a combination of clustering, cosmic shear and galaxy-galaxy lensing can provide a significant reduction in statistical errors from each analysis individually, however care must be taken to measure and model flux magnification. Finally, I consider how measurements of galaxy size and flux may be used to constrain the dark matter profile of a foreground lens, such as galaxy- or galaxy-cluster-dark matter halos. I present a method of constructing probability distributions for halo profile free parameters using Bayes’ Theorem, provided the intrinsic size-magnitude distribution may be measured from data. I investigate the use of this method on mock clusters, with an aim of investigating the precision and accuracy of returned parameter constraints under certain conditions. As part of this analysis, I quantify the size and significance of inaccuracies in the dark matter reconstruction as a result of limitations in the data from which the sample and size-magnitude distribution is obtained. This method is applied to public data from the Space Telescope A901/902 Galaxy Evolution Survey (STAGES), and results are presented for the four STAGES clusters using measurements of source galaxy size and magnitude, and a combination of both. I find consistent results with existing shear measurements using measurements of galaxy magnitudes, but interesting inconsistent results when galaxy size measurements are used. The simplifying assumptions and limitations of the analysis are discussed, and extensions to the method presented.

In this work we have studied the relationship between the kinematics and shapes of Early Type Galaxies (ETGs) in rich clusters. In particular we were interested to extend the kinematic morphology density relation to the richest clusters. We obtained data from FLAMES/GIRAFFE to probe the stellar kinematics of a sample of 30 ETGs in the massive cluster Abell 1689 at z = 0.183, to classify them as Slow Rotators (SRs) or Fast Rotators (Frs). To date, this is the highest redshift cluster studied in this way. We simulated FLAMES/GIRAFFE observations of the local SAURON galaxies to account for the bias introduced compared to the ATLAS3D sample, which we used as a local comparison. We find that the luminosity function of SRs in Abell 1689 is the same as that in ATLAS^3D, down to the faintest objects probed (M_K ≈ -23). The number fraction of SRs over the ETG population in Abell 1689 is f_SR = 0.15 +/- 0.03, consistent with the value found in the Virgo Cluster. However, within the cluster, f_SR rises sharply with the projected number density of galaxies, rising from f_SR = 0.01 in the least dense bin to f_SR = 0.58 in the densest bin. We conclude that the fraction of SRs is not determined by the local number density of galaxies, but rather by the physical location within the cluster. This might be due to dynamical processes which cause SRs (on average more massive) to sink in the gravitational potential of the cluster. Next we explore the distribution of projected ellipticity &epsilon; in galaxies belonging to a sample of clusters from SDSS (z </~ 0.1) and the CLASH survey (z ≈ 0.2). We were interested to establish whether the fraction of galaxies flatter than &epsilon; = 0.4 (a proxy for FRs) varies from cluster to cluster. We find some significant variations. We go on to probe the projected shape as a function of projected cluster-centric radius. In both samples we find that on average galaxies have progressively rounder projected shapes at lower cluster-centric projected distance. In the SDSS sample we show that this trend exists above and beyond the trend for brighter galaxies to be more common near the centre of clusters (bright galaxies are on average rounder). In order to disentangle the trend for SRs (which are rounder) to be more common near the centre of clusters, we isolate a subsample of FRs only, by considering only galaxies with &epsilon; > 0.4. We find that even the intrinsically flat FRs are on average rounder at lower projected cluster-centric distance. We conclude that the observed trend is due either to the dynamic heating of the stellar discs being strongest near the centre of clusters, or due to an anti-correlation of the bulge fractions with the cluster-centric distance.

This thesis explores the expansive world of General Relativity, and its role to play in modern cosmology and quantum field theory. We begin with a pedagogical approach to relativity, in particular, highlighting upon the ambiguity that arises with the conventions used in different textbooks. A brief introduction to tensor calculus has also been provided in the appendix. The preliminary chapters are also complimented with examples of numerical relativity via simulation. We then move on to discuss examples of non-linear systems, and their exact solutions. Such systems will be analogous to those we shall encounter later, upon considering scalar field theories as a means of modelling dark energy. We shall introduce the axion as our highly motivated dark matter candidate, since this will ultimately determine the behaviour of the scalar field. Coupled to a scaling factor across the spatial domain, it is found that this scalar field will ultimately determine the evolution of our universe. The key result of this thesis has been the possibility to screen both the cosmological constant, and flatness of the universe, to within observable parameters. These results will be explicitly derived from first principles. Also included is a tentative approach to holographic theory, in which strongly correlated systems may be modelled within the asymptotic domain of Anti-de Sitter (AdS) space. Ultimately, our aspirations are to bridge the gap with condensed matter theory, in particular with the publications included within the latter appendices. These publications discuss graphene as a revolutionary new material, for inclusion in both transistor-based and optoelectronic devices.

We address the issue of kinematic heating in disc galaxies by analysing a suite of cosmological Milky Way-type disc simulations run with different particle-and grid-based hydrodynamical codes and different resolution, and compare them with observations of the Milky Way. By studying the kinematics of disc stars in these simulations, we seek to determine whether or not the existence of a fragile thin disc is possible within a cosmological framework, where multiple mergers and interactions are the essence of galaxy formation. We study the velocity dispersion-age relation for disc stars at $z=0$ and find that four of the simulations, the stellar disc appears to undergo continual/secular heating. Two other simulations suggest a “saturation” in the heating profile for young stars in the disc. None of the simulations have thin discs as old as that of the Milky Way. We also analyse the kinematics of disc stars at the time of their birth, and find that in some simulations old stars are born cold within the disc and are subsequently heated, while other simulations possess old stellar populations, which are born relatively hot. The models which are in better agreement with observations of the Milky Way's stellar disc undergo significantly lower minor-merger/assembly activity after the last major merger. By running a set of isolated Milky Way-type simulations with different resolution and different density thresholds for star formation we conclude that, on top of the effects of mergers, there exists a ``floor'' in the dispersion that is related to the underlying treatment of the heating and cooling of the interstellar medium, and the low density threshold which such codes use for star formation. A persistent issue in simulations of disc galaxies is the formation of large spheroidal components, and disc galaxies with larger bulge to disc ratios than is observed. This problem is alleviated by supernova feedback. We found that by increasing the feedback in the simulations, we decrease the amount of stars that are accreted onto the main galaxy. The star formation is quenched more efficiently in low mass satellites when stronger feedback is implemented as well as in the main halo. These effects result in a disc galaxy, which has formed less stars overall, but more importantly, contains less accreted stars. As the strong stellar feedback quenches the star formation in the small building blocks, the metallicity of the accreted stars is lower than in the case where less feedback was used. In the context of hierarchical formation, mass assembly is expected to be scale free. Yet the properties of galaxies depend strongly on their mass. We examine how baryonic physics has different effects at different mass scales by analysing three cosmological simulations using the same initial conditions that are scaled to three different masses. Despite their identical dark matter merger history, we show that the simulated galaxies have significantly different stellar accretion histories. As we go down in mass, the lowest mass progenitors are unable to form stars, resulting in a low mass galaxy with less accreted stars. The overall chemical properties are also distinct at the different mass scales, as one might expect from the mass-metallicity relation of observed galaxies. We examine gradients of chemical abundances with radius and with height above the disc, and look for properties that are retained at different mass scales and properties which change, often dramatically. We analyse the kinematic and chemical properties of their accreted and in-situ populations. Again, trends can be found that persist at all mass scales, providing signatures of hierarchical structure formation. We find that accreted populations in the high mass simulation did not resemble any of the populations in the lower mass galaxies, showing that the chemical properties of proto-galaxies, which merge at high redshift to form massive galaxies, differ from the properties of low mass galaxies that survive at z=0. We probe further the signatures of hierarchical structure formation at smaller scales, in dwarf galaxies. We analysed the morphologies, kinematics and chemical properties of two simulated dwarf galaxies with different merger histories. We again analyse the accreted and in-situ populations. Observations of dwarf galaxies have found that they are comprised of multiple components. Our simulated dwarfs indicate that such populations may indeed be a manifestation of the hierarchical formation process in action in these lower mass galaxies. In one simulated dwarf, the in-situ stellar component forms a thin disc and a thick disc. We show that the thick disc in this simulation forms from in-situ stars that are born kinematically hot in the disc from early gas-rich mergers. The thin disc is formed quiescently from the later infall of gas. The accreted stars in the simulation were found to form an extended stellar halo. Chemical signatures of the three populations are also explored. The second dwarf we analysed has different galactic components, a result found to be due to the different merger history of this galaxy. The last major merger in this simulation occurs early on in the formation process between two proto-galaxies of similar mass. The result is a dwarf galaxy comprised of a disc formed of in-situ stars and a flattened rotating stellar halo formed of accreted stars. The angular momentum of the accreted and old in-insitu stars is obtained from the last major merger. We discuss the resemblance of this flattened rotating stellar halo to fast rotating flattened elliptical galaxies, and propose that such structures may explain some of the observed extra-galactic thick discs. These studies show that galactic properties emerge through the complex inter-play between hierarchical structure formation, star formation, and feedback from supernovae. Different modelling of these processes will alter the simulated galaxy's properties, and detailed comparisons with observations can then be made to determine the dominant processes responsible for different galactic properties. We remain optimistic that further improvement in modelling will allow deeper insights into the processes of galaxy formation and evolution.

Environment is known to have a significant impact on the evolution of galaxies. This is most evident in the local Universe, where the oldest and most massive galaxies are found at the of massive galaxy clusters. Current theory predicts that galaxies will form earlier and evolve more rapidly in the densest regions of the Universe. What is not clear is how rapidly the of environment start to have an impact on galaxies, at what stage can we detect physical differences between galaxies in dense regions and those in the field? By the time galaxies are assembled in virialised clusters the effects are clear, but at higher redshift (z ≳ 2), in the unvirialised progenitors of clusters (protoclusters) the effects are harder to detect. In this thesis I study the impact of environment in a z =3.1 protocluster in the SSA22 field. I consider the fraction of mergers in the protocluster, comparing it to the fraction of mergers in field at a similar redshift. My classification is based on the morphology of Lyman break (LBGs), using HST ACS/F814W imaging, which probes the rest frame UV. I find a marginal enhancement of the merger fraction, 48±10 per cent for LBGs in the protocluster compared 30±6 per cent in the field, suggesting that galaxy-galaxy mergers are one of the key driving accelerated star formation and AGN growth in protocluster environments. Having considered the fraction of mergers in the protocluster I then turn my attention to the physical properties of LBGs. I use multiwavelength data and spectral energy distribution fitting to determine the mass of LBGs in the protocluster and in the field. I find no statistical evidence for an enhancement of mass in the protocluster, suggesting that the protocluster environment has not impacted the average mass of LBGs at this redshift. It is possible that the protocluster LBG population may become more massive than LBGs in the field at lower redshift, or the galaxies may cease to be detectable by the Lyman break method before a mass difference between the protocluster galaxies and field is observable. Finally I consider the Lyman-α blobs (LABs) within the protocluster. These are large (~10- 100kpc) scale regions of diffuse Lyman-α emission, thought to be associated with overdense regions. 35 LABs have been detected in the SSA22 protocluster, indicating the presence of large clouds of gas in the circumgalactic medium. A debate has arisen regarding the powering mechanism of the LABs, particularly between star forming processes (e.g. Lyman-α escaping from a star forming galaxy or photoionizing radiation escaping from a star forming galaxy or active galactic nuclei) and a cold accretion model. The latter involves gas gravitationally cooling as it falls into the centre of a dark matter halo to feed a central galaxy. Some of this energy heats the cold gas, which then emits Lyman-α as it cools. The cold gas accretion theory gained popularity because some LABs appear not to contain a luminous galaxy or AGN which could explain the observed emission. One suggestion is that the central galaxy could be hidden by dust and that this could explain the lack of a detection in UV or optical. I therefore use SCUBA2 850μm imaging to search for submm sources (dusty star forming galaxies) in the LABs. I detect submm sources in only two of the LABs at 3.5δ, however, stacking all the LABs gives an average flux density of S850 = 0.6±0.2mJy. This suggests that on average the LABs do contain a submm source which could be a dusty galaxy. However, stacking the LABs by size indicates that only the largest third (area ≳ 1800kpc²) have a mean detection, at 4.5δ, with S850 = 1.4±0.3mJy, suggesting that different mechanisms may dominate the larger and smaller LAB populations. I explore two possible mechanisms for powering the LABs, cold accretion and central star forming galaxies. I find that central star formation is more likely to be the dominant source of emission, with cold accretion playing a secondary role.

The expansion of the universe is currently accelerating, as first inferred by Efstathiou et al. (1990), Ostriker & Steinhardt (1995) and directly determined by Riess et al. (1998) and Perlmutter et al. (1999). Current constraints are consistent with a time independent equation-of-state of w = -1, which is to be expected when a constant vacuum energy density dominates. But the Quantum Field Theory prediction for the magnitude of this vacuum energy is very much larger than that inferred (Weinberg, 1989; Koksma & Prokopec, 2011). It is entirely possible that the cause of the expansion has an alternative explanation, with both the inclusion of a quantum scalar field and modified gravity theories able to reproduce an expansion history close to, but potentially deviating from, that of a cosmological constant and cold dark matter. In this work I investigate the consistency of the VIMOS Public Extragalactic Redshift Survey (VIPERS) v7 census of the galaxy distribution at z = 0:8 with the expansion history and linear growth rate predicted by General Relativity (GR) when a Planck Collaboration et al. (2016) fiducial cosmology is assumed. To do so, I measure the optimally weighted redshift-space power spectrum (Feldman et al., 1994), which is anisotropic due to the coherent infall of galaxies towards overdensities and outflow from voids (Kaiser, 1987). The magnitude of this anisotropy can distinguish between modified theories of gravity as the convergence (divergence) rate of the velocity field depends on the effective strength of gravity on cosmological scales (Guzzo et al., 2008). This motivates measuring the linear growth rate rather than the background expansion, which is indistinguishable for a number of modified gravity theories. In Chapter 6 I place constraints of fσ8(0:76) = 0:44 ± 0:04; fσ8(1:05) = 0:28 ± 0:08; with the completed VIPERS v7 survey; the combination remains consistent with General Relativity at 95% confidence. The dependence of the errors on the assumed priors will be investigated in future work. Further anisotropy is introduced by the Alcock-Paczynski effect - a distortion of the observed power spectrum due to the assumption of a fiducial cosmology differing from the true one. These two sources of anisotropy may be separated based on their distinct scale and angular dependence with sufficiently precise measurements. Doing so degrades the constraints: fσ8(0:76) = 0:31 ± 0:10; fσ8(1:05) = -0:04 ± 0:26; but allows for the background expansion (FAP ≡ (1 + z)DAH=c) to be simultaneously constrained. Galaxy redshift surveys may then directly compare both the background expansion and linear growth rate to the GR predictions I find the VIPERS v7 joint-posterior on (fσ8; FAP ) shows no compelling deviation from the GR expectation although the sizeable errors reduce the significance of this conclusion. In Chapter 4 I describe and outline corrections for the VIPERS spectroscopic selection, which enable these constraints to be made. The VIPERS selection strategy is (projected) density dependent and may potentially bias measures of galaxy clustering. Throughout this work I present numerous tests of possible systematic biases, which are performed with the aid of realistic VIPERS mock catalogues. These also allow for accurate statistical error estimates to be made { by incorporating the sample variance due to both the finite volume and finite number density. Chapter 5 details the development and testing of a new, rapid approach for the forward modelling of the power spectrum multipole moments obtained from a survey with an involved angular mask. An investigation of the necessary corrections for the VIPERS PDR-1 angular mask is recorded. This includes an original derivation for the integral constraint correction for a smoothed, joint-field estimate of ¯n(z) and a description of how the mask should be accounted for in light of the Alcock- Paczynski effect. Chapter 7 investigates the inclusion of a simple local overdensity transform: 'clipping' prior to the redshift-space distortions (RSD) analysis. This tackles the root cause of non-linearity and potentially extends the validity of perturbation theory. Moreover, this marked clustering statistic potentially amplifies signatures of modified gravity and, as a density-weighted two-point statistic, includes information not available to the power spectrum. I show that a linear real-space power spectrum with a Kaiser factor and a Lorentzian damping yields a significant bias without clipping, but that this may be removed with a sufficiently strict transform; similar behaviour is observed for the VIPERS v7 dataset. Estimates of fσ8 for different thresholds are highly correlated due to the overlapping volume, but the bias for insufficient clipping can be calibrated and the correlation obtained using mock catalogues. A maximum likelihood value for the combined constraint of a number of thresholds is shown to achieve a ' 16% decrease in statistical error relative to the most precise single-threshold estimate. The results are encouraging to date but represent a work in progress; the final analysis will be submitted to Astronomy & Astrophysics as Wilson et al. (2016). In addition to this, an original extension of the prediction for a clipped Gaussian field to a clipped lognormal field is presented. The results of tests of this model with a real-space cube populated according to the halo occupation distribution model are also provided.

This thesis presents the discovery of the first luminous z & 6.5 quasars in the VISTA kilo-degree Infrared Galaxy Survey (VIKING). After some basic quality control, quasar selection is investigated via use of initial data supplementedwith detailed modelling of the photometric and spatial distributions of stars of spectral type M, L and T, which are known to be the cause of significant contamination in quasar colour selection spaces. Optimised selection constraints are placed on detection significance and morphology and the performance of a traditional colour selection technique is compared to a Bayesian model comparison technique. The latter is found to offer a ∼10 per cent gain in completeness over traditional colour selection. Quasar candidates are ranked via Bayesian model comparison and a subset of the highest ranked objects are put forward for follow-up imaging. In June 2011, 44 high-z quasar candidates underwent deep optical i- and z- band imaging on the ESO NTT. Just 6 of these candidates were found to have optical colours consistent z & 6.5 quasars. Spectroscopic follow-up of these objects is ongoing, but thus far three new quasars have been discovered at redshifts of z=6.5, 6.7, 6.9. This discovery rate is consistentwith zero evolution in the rate of decline in quasar space density from z & 6.4. This differs fromthe latest results from UKIDSS. Further results expected from these and other surveys will begin to constrain the true nature of quasar space density evolution in the near future. The discovery of three z ≥ 6.5 quasars in VIKING is a significant highlight in the first year of VISTA science operations. These quasars will remain important probes of the high-z universe throughout the next decade.

In this Thesis, I explore the formation and evolution of stellar spiral structures embedded in realistic dark matter haloes with very high resolution simulations. I first study the impact of the shape of the dark matter haloes. I find that non-adiabatic changes to the dark matter halo shape, commonly found in cosmological simulations due to the assembly history of haloes, can trigger strong two-armed grand-design spiral structures extending from the inner disc to the outer region. The nature of the spiral structures is found to be consistent with kinematic density waves based on the study of their power spectra. Such grand-design spiral structures may help the formation of transient multi-armed spiral structures if the self-gravity in disc is strong enough. Evolution of spiral structures is similar when the disc and the halo are misaligned, although warps develop additionally. I further find a strong correlation between the torque strength from the halo and the strength of the corresponding spiral structures. In the second part of my Thesis I then study the influence of subhaloes by including them from realistic cosmological simulations. I identify five different massive subhaloes that hit the central region of the disc, two out of which hit the disc twice. Aside from disc heating, three distinct generations of spiral structures are found in the stellar disc, which can be related to different subhaloes. For each generation, counter-rotating single-armed spiral structures develop first. They wind up very quickly before two-armed spiral structures become prominent. These spiral structures are again identified as kinematic density waves. We find that rather than interacting with the disc through resonances, subhaloes preferentially trigger spiral structures impulsively, due to their relatively short impact time with the disc. The strength of spiral structures can be related to the integrated strength of the torque generated by subhaloes. The correlation between the torque strength exerted by a triaxial dark matter halo and by subhaloes and the spiral strength may provide constraints on the distribution of dark matter.