We present a detailed analysis of a large sample of spectroscopically confirmed ultra-massive quiescent galaxies (${\rm{log}}(M_{\ast}/M_{\odot})\sim11.5$) at $z\gtrsim2$. This sample comprises 15 galaxies selected in the COSMOS and UDS fields by their bright K-band magnitudes and followed up with VLT/X-shooter spectroscopy and HST/WFC3 $H_{F160W}$ imaging. These observations allow us to unambiguously confirm their redshifts ascertain their quiescent nature and stellar ages, and to reliably assess their internal kinematics and effective radii. We find that these galaxies are compact, consistent with the high mass end of the mass-size relation for quiescent galaxies at $z=2$. Moreover, the distribution of the measured stellar velocity dispersions of the sample is consistent with the most massive local early-type galaxies from the MASSIVE Survey showing that evolution in these galaxies, is dominated by changes in size. The HST images reveal, as surprisingly high, that $40\ \%$ of the sample have tidal features suggestive of mergers and companions in close proximity, including three galaxies experiencing ongoing major mergers. The absence of velocity dispersion evolution from $z=2$ to $0$, coupled with a doubling of the stellar mass, with a factor of four size increase and the observed disturbed stellar morphologies support dry minor mergers as the primary drivers of the evolution of the massive quiescent galaxies over the last 10 billion years.

Galaxy mergers are expected to have a significant role in the mass assembly of galaxies in the early Universe, but there are very few observational constraints on the merger history of galaxies at $z>2$. We present the first study of galaxy major mergers (mass ratios $>$ 1:4) in mass-selected samples out to $z\approx6$. Using all five fields of the HST/CANDELS survey and a probabilistic pair count methodology that incorporates the full photometric redshift posteriors and corrections for stellar mass completeness, we measure galaxy pair-counts for projected separations between 5 and 30 kpc in stellar mass selected samples at $9.7 < \log_{10}(\rm{M}_{*}/\rm{M}_{\odot}) < 10.3$ and $\log_{10}(\rm{M}_{*}/\rm{M}_{\odot}) > 10.3$. We find that the major merger pair fraction rises with redshift to $z\approx6$ proportional to $(1+z)^{m}$, with $m = 0.8\pm0.2$ ($m = 1.8\pm0.2$) for $\log_{10}(\rm{M}_{*} / \rm{M}_{\odot}) > 10.3$ ($9.7 < \log_{10}(\rm{M}_{*}/\rm{M}_{\odot}) < 10.3$). Investigating the pair fraction as a function of mass ratio between 1:20 and 1:1, we find no evidence for a strong evolution in the relative numbers of minor to major mergers out to $z<3$. Using evolving merger timescales we find that the merger rate per galaxy ($\mathcal{R}$) rises rapidly from $0.07\pm 0.01$ Gyr$^{-1}$ at $z < 1$ to $7.6\pm 2.7$ Gyr$^{-1}$ at $z = 6$ for galaxies at $\log_{10}(\rm{M}_{*}/\rm{M}_{\odot}) > 10.3$. The corresponding co-moving major merger rate density remains roughly constant during this time, with rates of $\Gamma \approx 10^{-4}$ Gyr$^{-1}$ Mpc$^{-3}$. Based on the observed merger rates per galaxy, we infer specific mass accretion rates from major mergers that are comparable to the specific star-formation rates for the same mass galaxies at $z>3$ - observational evidence that mergers are as important a mechanism for building up mass at high redshift as in-situ star-formation.

We use a large sample of $\sim 350,000$ galaxies constructed by combining the UKIDSS UDS, VIDEO/CFHT-LS, UltraVISTA/COSMOS and GAMA survey regions to probe the major merging histories of massive galaxies ($>10^{10}\ \mathrm{M}_\odot$) at $0.005 < z < 3.5$. We use a method adapted from that presented in Lopez-Sanjuan et al. (2014) using the full photometric redshift probability distributions, to measure pair $\textit{fractions}$ of flux-limited, stellar mass selected galaxy samples using close-pair statistics. The pair fraction is found to weakly evolve as $\propto (1+z)^{0.8}$ with no dependence on stellar mass. We subsequently derive major merger $\textit{rates}$ for galaxies at $> 10^{10}\ \mathrm{M}_\odot$ and at a constant number density of $n > 10^{-4}$ Mpc$^{-3}$, and find rates a factor of 2-3 smaller than previous works, although this depends strongly on the assumed merger timescale and likelihood of a close-pair merging. Galaxies undergo approximately 0.5 major mergers at $z < 3.5$, accruing an additional 1-4 $\times 10^{10}\ \mathrm{M}_\odot$ in the process. Major merger accretion rate densities of $\sim 2 \times 10^{-4}$ $\mathrm{M}_\odot$ yr$^{-1}$ Mpc$^{-3}$ are found for number density selected samples, indicating that direct progenitors of local massive ($>10^{11}\mathrm{M}_\odot$) galaxies have experienced a steady supply of stellar mass via major mergers throughout their evolution. While pair fractions are found to agree with those predicted by the Henriques et al. (2014) semi-analytic model, the Illustris hydrodynamical simulation fails to quantitatively reproduce derived merger rates. Furthermore, we find major mergers become a comparable source of stellar mass growth compared to star-formation at $z < 1$, but is 10-100 times smaller than the SFR density at higher redshifts.

Due to significant galaxy contamination and impurity in stellar mass selected samples (up to 95% from z=0-3), we examine the star formation history, quenching time-scales, and structural evolution of galaxies using a constant number density selection with data from the UKIDSS Ultra-Deep Survey field. Using this methodology we investigate the evolution of galaxies at a variety of number densities from $z=0-3$. We find that samples chosen at number densities ranging from $3\times10^{-4}$ to 10$^{-5}$ galaxies Mpc$^{-3}$ (corresponding to $z\sim0.5$ stellar masses of M$_{*}= 10^{10.95-11.6}$ M$_{0}$) have a star forming blue fraction of $\sim50$\% at $z\sim2.5$, which evolves to a nearly $100$\% quenched red and dead population by $z\sim 1$. We also see evidence for number density downsizing, such that the galaxies selected at the lowest densities (highest masses) become a homogeneous red population before those at higher number densities. Examining the evolution of the colours for these systems furthermore shows that the formation redshift of galaxies selected at these number densities is $z_{\rm form}>3$. The structural evolution through size and Sersic index fits reveal that while there remains evolution in terms of galaxies becoming larger and more concentrated in stellar mass at lower redshifts, the magnitude of the change is significantly smaller than for a mass selected sample. We also find that changes in size and structure continues at $z < 1$, and is coupled strongly to passivity evolution. We conclude that galaxy structure is driving the quenching of galaxies, such that galaxies become concentrated before they become passive.

Connecting galaxies with their descendants (or progenitors) at different redshifts can yield strong constraints on galaxy evolution. Observational studies have historically selected samples of galaxies using a physical quantity, such as stellar mass, either above a constant limit or at a constant cumulative number density. Investigation into the efficacy of these selection methods has not been fully explored. Using a set of four semi-analytical models based on the output of the Millennium Simulation, we find that selecting galaxies at a constant number density (in the range $-4.3 < \log\ n\ [\mathrm{Mpc}^{-3}\ h^{3}] < -3.0$) is superior to a constant stellar mass selected sample, although it still has significant limitations. Recovery of the average stellar mass, stellar mass density and average star-formation rate is highly dependent on the choice of number density but can all be recovered to within $<50\%$ at the commonly employed choice of $\log\ n\ [\mathrm{Mpc}^{-3}\ h^{3}] = -4.0$, corresponding to $\log M_\odot / h \sim 11.2$ at $z=0$, but this increases at lower mass limits. We show that there is a large scatter between the location of a given galaxy in a rank ordering based on stellar mass between different redshifts. We find that the inferred velocity dispersion may be a better tracer of galaxy properties, although further investigation is warranted into simulating this property. Finally, we find that over large redshift ranges selection at a constant number density is more effective in tracing the progenitors of modern galaxies than vice-versa.

We present a study on the stellar mass growth of the progenitors of local massive galaxies with a variety of number density selections with $n\le1\times10^{-4}\,\rm{Mpc^{-3}}$ (corresponding to $M_*=10^{11.24}\rm{M_{\odot}}$ at z=0.3) in the redshift range $0.3<z<3.0$. We select the progenitors of massive galaxies using a constant number density selection, and one which is adjusted to account for major mergers. We find that the progenitors of massive galaxies grow by a factor of four in total stellar mass over this redshift range. On average the stellar mass added via the processes of star formation, major and minor mergers account for $24\pm8\%$, $17\pm15\%$ and $34\pm14\%$, respectively, of the total galaxy stellar mass at $z=0.3$. Therefore $51\pm20\%$ of the total stellar mass in massive galaxies at $z=0.3$ is created externally to their z=3 progenitors. We explore the implication of these results on the cold gas accretion rate and size evolution of the progenitors of most massive galaxies over the same redshift range. We find an average gas accretion rate of $\sim66\pm32\,\rm{M_{\odot}yr^{-1}}$ over the redshift range of $1.5<z<3.0$. We find that the size evolution of a galaxy sample selected this way is on average lower than the findings of other investigations.