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.

We present a physical characterization of MMJ100026.36+021527.9 (a.k.a. ``MAMBO-9''), a dusty star-forming galaxy (DSFG) at $z=5.850\pm0.001$. This is the highest redshift unlensed DSFG (and fourth most distant overall) found to-date, and is the first source identified in a new 2mm blank-field map in the COSMOS field. Though identified in prior samples of DSFGs at 850$\mu$m-1.2mm with unknown redshift, the detection at 2mm prompted further follow-up as it indicated a much higher probability that the source was likely to sit at $z>4$. Deep observations from the Atacama Large Millimeter and submillimeter Array (ALMA) presented here confirm the redshift through the secure detection of $^{12}$CO($J\!=$6$\rightarrow$5) and p-H$_{2}$O(2$_{1,1}\!\rightarrow$2$_{0,2}$). MAMBO-9 is comprised of a pair of galaxies separated by 6kpc with corresponding star-formation rates of 590M$_\odot$yr$^{-1}$ and 220M$_\odot$yr$^{-1}$ total molecular hydrogen gas mass of (1.7$\pm$0.4)$\times10^{11}$M$_\odot$, dust mass of (1.3$\pm$0.3)$\times10^{9}$M$_\odot$ and stellar mass of (3.2$^{+1.0}_{-1.5}$)$\times10^{9}$M$_\odot$. The total halo mass, (3.3$\pm$0.8)$\times10^{12}$M$_\odot$, is predicted to exceed $>10^{15}$M$_\odot$ by $z=0$. The system is undergoing a merger-driven starburst which will increase the stellar mass of the system tenfold in $\tau_{\rm depl}=40-80$Myr, converting its large molecular gas reservoir (gas fraction of 96$^{+1}_{-2}$%) into stars. MAMBO-9 evaded firm spectroscopic identification for a decade, following a pattern that has emerged for some of the highest redshift DSFGs found. And yet, the systematic identification of unlensed DSFGs like MAMBO-9 is key to measuring the global contribution of obscured star-formation to the star-formation rate density at $z>4$, the formation of the first massive galaxies, and the formation of interstellar dust at early times ($<$1Gyr).

The next generation of electromagnetic and gravitational wave observatories will open unprecedented windows to the birth of the first supermassive black holes. This has the potential to reveal their origin and growth in the first billion years, as well as the signatures of their formation history in the local Universe. With this in mind, we outline three key focus areas which will shape research in the next decade and beyond: (1) What were the "seeds" of the first quasars; how did some reach a billion solar masses before z$\sim7$? (2) How does black hole growth change over cosmic time, and how did the early growth of black holes shape their host galaxies? What can we learn from intermediate mass black holes (IMBHs) and dwarf galaxies today? (3) Can we unravel the physics of black hole accretion, understanding both inflows and outflows (jets and winds) in the context of the theory of general relativity? Is it valid to use these insights to scale between stellar and supermassive BHs, i.e., is black hole accretion really scale invariant? In the following, we identify opportunities for the Canadian astronomical community to play a leading role in addressing these issues, in particular by leveraging our strong involvement in the Event Horizon Telescope, the {\it James Webb Space Telescope} (JWST), Euclid, the Maunakea Spectroscopic Explorer (MSE), the Thirty Meter Telescope (TMT), the Square Kilometer Array (SKA), the Cosmological Advanced Survey Telescope for Optical and ultraviolet Research (CASTOR), and more. We also discuss synergies with future space-based gravitational wave (LISA) and X-ray (e.g., Athena, Lynx) observatories, as well as the necessity for collaboration with the stellar and galactic evolution communities to build a complete picture of the birth of supermassive black holes, and their growth and their influence over the history of the Universe.

Optically-compact star-forming galaxies (SFGs) have been proposed as immediate progenitors of quiescent galaxies, although their origin and nature are debated. Were they formed in slow secular processes or in rapid merger-driven starbursts? Addressing this question would provide fundamental insight into how quenching occurs. We explore the location of the general population of galaxies with respect to fundamental star-forming and structural relations, identify compact SFGs based on their stellar core densities, and study three diagnostics of the burstiness of star formation: 1) Star formation efficiency, 2) interstellar medium (ISM), and 3) radio emission. The overall distribution of galaxies in the fundamental relations points towards a smooth transition towards quiescence while galaxies grow their stellar cores, although some galaxies suddenly increase their specific star-formation rate when they become compact. From their star formation efficiencies compact and extended SFGs appear similar. In relation to the ISM diagnostic, by studying the CO excitation, the density of the neutral gas, and the strength of the ultraviolet radiation field, compact SFGs resemble galaxies located in the upper envelope of the SFGs main sequence, although yet based on a small sample size. Regarding the radio emission diagnostic we find that galaxies become increasingly compact as the starburst ages, implying that at least some compact SFGs are old starbursts. We suggest that compact SFGs could be starburts winding down and eventually crossing the main sequence towards quiescence.

We use MUSE/VLT to conduct a survey of $z\sim3$ physical quasar pairs at close separation with a fast observation strategy. Our aim is twofold: (i) explore the Ly$\alpha$ glow around the faint-end of the quasar population; (ii) take advantage of the combined illumination of a quasar pair to unveil large-scale intergalactic structures extending between the two quasars. Here, we report the results for a quasar pair ($z=3.020,3.008$; $i=21.84,22.15$), separated by 11.6 arcsec (or 89 projected kpc). MUSE reveals filamentary Ly$\alpha$ structures extending between the two quasars with an average surface brightness of SB$_{\rm Ly\alpha}=1.8\times10^{-18}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$. Photoionization models of the constraints in the Ly$\alpha$, HeII, and CIV line emissions show that the emitting structures are intergalactic bridges with an extent between $\sim89$ and up to $\sim600$ kpc. Our models rule out the possibility that the structure extends for $\sim 2.9$ Mpc, i.e., the separation inferred from the uncertain systemic redshift difference of the quasars if the difference was only due to the Hubble flow. At the current spatial resolution and surface brightness limit, the average projected width of an individual bridge is about 35 kpc. We also detect a strong absorption in HI, NV, and CIV along the background sight-line at higher $z$, which we interpret as due to at least two components of cool, metal enriched, and relatively ionized CGM or IGM surrounding the quasar pair. Two additional HI absorbers are detected along both quasar sight-lines at $\sim -900$ and $-2800$ km s$^{-1}$ from the system, with the latter having associated CIV absorption only along the foreground quasar sight-line. The absence of galaxies in the MUSE field of view at the redshifts of these two absorbers suggests that they trace large-scale structures or expanding shells in front of the quasar pair.

We have used the Multi-Unit Spectroscopic Explorer (MUSE) to study the circumgalactic medium (CGM) of a z = 2.92 radio galaxy, MRC 0943-242 by parametrising its emitting and absorbing gas. In both Ly$\alpha$ $\lambda$1216 and He II $\lambda$1640 lines, we observe emission with velocity shifts of $\Delta v \simeq-1000$ km s$^{-1}$ from the systemic redshift of the galaxy. These blueshifted components represent kinematically perturbed gas that is aligned with the radio axis, which we interpret as jet-driven outflows. Three of the four known Ly$\alpha$ absorbers are detected at the same velocity as C IV $\lambda\lambda1548,1551$ and N V $\lambda\lambda1239,1243$ absorbers, proving that the gas is metal enriched more so than previously thought. At the velocity of a strong Ly$\alpha$ absorber with an HI column of $N_{\rm HI}/{\rm cm}^{-2} = 10^{19.2}$ and velocity shift of $\Delta v \simeq -400$ km s$^{-1},$ we also detect Si II $\lambda$1260 and Si II $\lambda$1527 absorption, which suggests that the absorbing gas is ionisation bounded. With the added sensitivity of this MUSE observation, we are more capable of adding constraints to absorber column densities and consequently determining what powers their ionisation. To do this, we obtain photoionisation grid models in \pkg{cloudy} which show that AGN radiation is capable of ionising the gas and producing the observed column densities in a gas of metallicity of Z/Z$_\odot \simeq$ 0.01 with a nitrogen abundance a factor of 10 greater than that of hydrogen. This metal-enriched absorbing gas, which is also spatially extended over a projected distance of $r \gtrsim 60$ kpc, is likely to have undergone chemical enrichment through stellar winds that have swept up metals from the interstellar-medium and deposited them in the outer regions of the galaxy's halo.

The objective of this work is to study how active galactic nuclei (AGN) influence star formation in host galaxies. We present a detailed investigation of the star-formation history and conditions of a $z=2.57$ massive radio galaxy based on VLT/X-SHOOTER and ALMA observations. The deep rest-frame ultraviolet spectrum contains photospheric absorption lines and wind features indicating the presence of OB-type stars. The most significantly detected photospheric features are used to characterize the recent star formation: neither instantaneous nor continuous star-formation history is consistent with the relative strength of the Si II $\lambda$1485 and S V $\lambda$1502 absorption. Rather, at least two bursts of star formation took place in the recent past, at $6^{+1}_{-2}$ Myr and $\gtrsim20$ Myr ago, respectively. We deduce a molecular H$_{2}$ gas mass of $(3.9\pm1.0)\times10^{10}$ M$_{\odot}$ based on ALMA observations of the [C I] $^3$P$_{2}$-$^3$P$_{1}$ emission. The molecular gas mass is only 13 % of its stellar mass. Combined with its high star-formation rate of ($1020^{+190}_{-170}$) Myr, this implies a high star-formation efficiency of $(26\pm8$) Gyr$^{-1}$ and a short depletion time of $(38\pm12)$ Myr. We attribute the efficient star formation to compressive gas motions in order to explain the modest velocity dispersions ($\leqslant$ 55 km s$^{-1}$) of the photospheric lines and of the star-forming gas traced by [C I]. Because of the likely very young age of the radio source, our findings suggest that vigorous star formation consumes much of the gas and works in concert with the AGN to remove any residual molecular gas, and eventually quenching star formation in massive galaxies.

Understanding how and why star formation turns off in massive galaxies is a major challenge for studies of galaxy evolution. Many theoretical explanations have been proposed, but a definitive consensus is yet to be reached.

We present a multiwavelength study of the AGN-host starburst galaxy PKS 0529-549 at z~2.6. We use (1) new ALMA observations of the dust continuum and of the [CI] 370 um line, tracing molecular gas, (2) SINFONI spectroscopy of the [OIII] 5007 Ang line, tracing ionized gas, and (3) ATCA radio continuum images, tracing synchrotron emission. Both [CI] and [OIII] show regular velocity gradients, but their systemic velocities and position angles differ by ~300 km/s and ~30 degrees, respectively. The [CI] is consistent with a rotating disc, aligned with the dust and stellar continuum, while the [OIII] likely traces an outflow, aligned with two AGN-driven radio lobes. We model the [CI] cube using 3D disc models, which give best-fit rotation velocities V~310 km/s and velocity dispersions sigma<30 km/s. Hence, the [CI] disc has V/sigma>10 and is not particularly turbulent, similar to local galaxy discs. The dynamical mass (~10^11 Msun) is comparable to the baryonic mass within the errors. This suggests that baryons dominate the inner galaxy dynamics, similar to massive galaxies at z=0. Remarkably, PKS 0529-549 lies on the local baryonic Tully-Fisher relation, indicating that at least some massive galaxies are already in place and kinematically relaxed at z~2.6. This work highlights the potential of the [CI] line to trace galaxy dynamics at high z, as well as the importance of multiwavelength data to interpret gas kinematics.

We report the discovery of eMACSJ1341-QG-1, a quiescent galaxy at $z=1.594$ located behind the massive galaxy cluster eMACSJ1341.9$-$2442 ($z=0.835$). The system was identified as a gravitationally lensed triple image in Hubble Space Telescope images obtained as part of a snapshot survey of the most X-ray luminous galaxy clusters at $z>0.5$ and spectroscopically confirmed in ground-based follow-up observations with the ESO/X-Shooter spectrograph. From the constraints provided by the triple image, we derive a first, crude model of the mass distribution of the cluster lens, which predicts a gravitational amplification of a factor of $\sim$30 for the primary image and a factor of $\sim$6 for the remaining two images of the source, making eMACSJ1341-QG-1 by far the most strongly amplified quiescent galaxy discovered to date. Our discovery underlines the power of SNAPshot observations of massive, X-ray selected galaxy clusters for lensing-assisted studies of faint background populations.

At redshift z = 2, when the Universe was just three billion years old, half of the most massive galaxies were extremely compact and had already exhausted their fuel for star formation(1-4). It is believed that they were formed in intense nuclear starbursts and that they ultimately grew into the most massive local elliptical galaxies seen today, through mergers with minor companions(5,6), but validating this picture requires higher-resolution observations of their centres than is currently possible. Magnification from gravitational lensing offers an opportunity to resolve the inner regions of galaxies(7). Here we report an analysis of the stellar populations and kinematics of a lensed z = 2.1478 compact galaxy, which surprisingly turns out to be a fast-spinning, rotationally supported disk galaxy. Its stars must have formed in a disk, rather than in a merger-driven nuclear starburst(8). The galaxy was probably fed by streams of cold gas, which were able to penetrate the hot halo gas until they were cut off by shock heating from the dark matter halo(9). This result confirms previous indirect indications(10-13) that the first galaxies to cease star formation must have gone through major changes not just in their structure, but also in their kinematics, to evolve into present-day elliptical galaxies.

We present ALMA detections of the [CI] 1-0, CO J=3-2, and CO J=4-3 emission lines, as well as the ALMA band 4 continuum for a compact star-forming galaxy (cSFG) at z=2.225, 3D-HST GS30274. As is typical for cSFGs, this galaxy has a stellar mass of $1.89 \pm 0.47\,\times 10^{11}\,\rm{M}_\odot$, with a star formation rate of $214\pm44\,\rm{M}_\odot\,\rm{yr}^{-1}$ putting it on the star-forming `main-sequence', but with an H-band effective radius of 2.5 kpc, making it much smaller than the bulk of `main-sequence' star-forming galaxies. The intensity ratio of the line detections yield an ISM density (~ 6 $\times 10^{4}\,\rm{cm}^{-3}$) and a UV-radiation field ( ~2 $\times 10^4\,\rm{G}_0$), similar to the values in local starburst and ultra-luminous infrared galaxy environments. A starburst phase is consistent with the short depletion times ($t_{\rm H2, dep} \leq 140$ Myr) we find using three different proxies for the H2 mass ([CI], CO, dust mass). This depletion time is significantly shorter than in more extended SFGs with similar stellar masses and SFRs. Moreover, the gas fraction of 3D-HST GS30274 is smaller than typically found in extended galaxies. We measure the CO and [CI] kinematics and find a FWHM line width of ~$750 \pm 41 $ km s$^{-1}$. The CO and [CI] FWHM are consistent with a previously measured H$\alpha$ FWHM for this source. The line widths are consistent with gravitational motions, suggesting we are seeing a compact molecular gas reservoir. A previous merger event, as suggested by the asymmetric light profile, may be responsible for the compact distribution of gas and has triggered a central starburst event. This event gives rise to the starburst-like ISM properties and short depletion times. The centrally located and efficient star formation is quickly building up a dense core of stars, responsible for the compact distribution of stellar light in 3D-HST GS30274.

We present rest-frame optical spectra from the FMOS-COSMOS survey of twelve $z \sim 1.6$ \textit{Herschel} starburst galaxies, with Star Formation Rate (SFR) elevated by $\times$8, on average, above the star-forming Main Sequence (MS). Comparing the H$\alpha$ to IR luminosity ratio and the Balmer Decrement we find that the optically-thin regions of the sources contain on average only $\sim 10$ percent of the total SFR whereas $\sim90$ percent comes from an extremely obscured component which is revealed only by far-IR observations and is optically-thick even in H$\alpha$. We measure the [NII]$_{6583}$/H$\alpha$ ratio, suggesting that the less obscured regions have a metal content similar to that of the MS population at the same stellar masses and redshifts. However, our objects appear to be metal-rich outliers from the metallicity-SFR anticorrelation observed at fixed stellar mass for the MS population. The [SII]$_{6732}$/[SII]$_{6717}$ ratio from the average spectrum indicates an electron density $n_{\rm e} \sim 1,100\ \mathrm{cm}^{-3}$, larger than what estimated for MS galaxies but only at the 1.5$\sigma$ level. Our results provide supporting evidence that high-$z$ MS outliers are the analogous of local ULIRGs, and are consistent with a major merger origin for the starburst event.

We present and analyse an extensive dataset of the superluminous supernova (SLSN) LSQ14mo (z = 0.256), consisting of a multi-colour lightcurve from -30 d to +70 d in the rest-frame and a series of 6 spectra from PESSTO covering -7 d to +50 d. This is among the densest spectroscopic coverage, and best-constrained rising lightcurve, for a fast-declining hydrogen-poor SLSN. The bolometric lightcurve can be reproduced with a millisecond magnetar model with ~ 4 M_sol ejecta mass, and the temperature and velocity evolution is also suggestive of a magnetar as the power source. Spectral modelling indicates that the SN ejected ~ 6 M_sol of CO-rich material with a kinetic energy of ~ 7 x 10^51 erg, and suggests a partially thermalised additional source of luminosity between -2 d and +22 d. This may be due to interaction with a shell of material originating from pre-explosion mass loss. We further present a detailed analysis of the host galaxy system of LSQ14mo. PESSTO and GROND imaging show three spatially resolved bright regions, and we used the VLT and FORS2 to obtain a deep (five-hour exposure) spectra of the SN position and the three star-forming regions, which are at a similar redshift. The FORS spectrum at +300 days shows no trace of SN emission lines and we place limits on the strength of [O I] from comparisons with other Ic SNe. The deep spectra provides a unique chance to investigate spatial variations in the host star-formation activity and metallicity. The specific star-formation rate is similar in all three components, as is the presence of a young stellar population. However, the position of LSQ14mo exhibits a lower metallicity, with 12 + log(O/H) = 8.2 in both the R23 and N2 scales (corresponding to ~ 0.3 Z_sol). We propose that the three bright regions in the host system are interacting, which thus triggers star-formation and forms young stellar populations.

While performing a survey to detect recoiling supermassive black holes, we have identified an unusual source having a projected offset of 800 pc from a nearby dwarf galaxy. The object, SDSS J113323.97+550415.8, exhibits broad emission lines and strong variability. While originally classified as a supernova (SN) because of its nondetection in 2005, we detect it in recent and past observations over 63 yr and find over a magnitude of rebrightening in the last 2 years. Using high-resolution adaptive optics observations, we constrain the source emission region to be <12 pc and find a disturbed host-galaxy morphology indicative of recent merger activity. Observations taken over more than a decade show narrow [O III] lines, constant ultraviolet emission, broad Balmer lines, a constant putative black hole mass over a decade of observations despite changes in the continuum, and optical emission-line diagnostics consistent with an active galactic nucleus (AGN). However, the optical spectra exhibit blueshifted absorption, and eventually narrow Fe II and [Ca II] emission, each of which is rarely found in AGN spectra. While this peculiar source displays many of the observational properties expected of a potential black hole recoil candidate, some of the properties could also be explained by a luminous blue variable star (LBV) erupting for decades since 1950, followed by a Type IIn SN in 2001. Interpreted as an LBV followed by a SN analogous to SN 2009ip, the multi-decade LBV eruptions would be the longest ever observed, and the broad Halpha emission would be the most luminous ever observed at late times (>10 yr), larger than that of unusually luminous supernovae such as SN 1988Z, suggesting one of the most extreme episodes of pre-SN mass loss ever discovered.

We present stringent constraints on the average mid-, far-infrared and radio emissions of $\sim$14200 quiescent galaxies (QGs), identified out to $z=3$ in the COSMOS field via their rest-frame NUV$-$r and r$-$J colors, and with stellar masses $M_{\star}=10^{9.8-12.2} \,M_{\odot} $. Stacking in deep Spitzer (MIPS $24\,\mu$m), Herschel (PACS and SPIRE), and VLA (1.4 GHz) maps reveals extremely low dust-obscured star formation rates for QGs (SFR $<0.1-3\,M_{\odot}$yr$^{-1}$ at $z \leqslant 2$ and $<6-18\,M_{\odot}$yr$^{-1}$ at $z > 2$), consistent with the low unobscured SFRs ($<0.01-1.2\,M_{\odot}$yr$^{-1}$) inferred from modeling their ultraviolet-to-near-infrared photometry. The average SFRs of QGs are $>10\times$ below those of star-forming galaxies (SFGs) within the $M_{\star}$- and $z$-ranges considered. The stacked 1.4 GHz signals (S/N $> 5$) are, if attributed solely to star formation, in excess of the total (obscured plus unobscured) SFR limits, suggestive of a widespread presence of low-luminosity active galactic nuclei (AGN) among QGs. Our results reaffirm the existence of a significant population QGs out to $z = 3$, thus corroborating the need for powerful quenching mechanism(s) to terminate star formation in galaxies at earlier epochs.

We present deep, near-infrared HST/WFC3 grism spectroscopy and imaging for a sample of 14 galaxies at z~2 selected from a mass-complete photometric catalog in the COSMOS field. By combining the grism observations with photometry in 30 bands, we derive accurate constraints on their redshifts, stellar masses, ages, dust extinction and formation redshifts. We show that the slope and scatter of the z~2 mass-size relation of quiescent galaxies is consistent with the local relation, and confirm previous findings that the sizes for a given mass are smaller by a factor of two to three. Finally, we show that the observed evolution of the mass-size relation of quiescent galaxies between z=2 and 0 can be explained by quenching of increasingly larger star-forming galaxies, at a rate dictated by the increase in the number density of quiescent galaxies with decreasing redshift. However, we find that the scatter in the mass-size relation should increase in the quenching-driven scenario in contrast to what is seen in the data. This suggests that merging is not needed to explain the evolution of the median mass-size relation of massive galaxies, but may still be required to tighten its scatter, and explain the size growth of individual z=2 galaxies quiescent galaxies.

We measure the merger fraction of massive galaxies using the UltraVISTA/COSMOS $Ks$-band selected catalog, complemented with the deeper, higher resolution 3DHST+CANDELS catalog selected in the HST/WFC3 $H$-band, presenting the largest mass-complete photometric merger sample up to $z\sim3$. We find that selecting mergers using the $H_{160}$-band flux ratio leads to an increasing merger fraction with redshift, while selecting mergers using the stellar mass ratio causes a diminishing redshift dependence. Defining major and minor mergers as having stellar mass ratios of 1:1 - 4:1 and 4:1 - 10:1 respectively, the results imply $\sim$1 major and $\lesssim$1 minor merger for an average massive (log$(M_{\star}/M_{\odot}) \geqslant 10.8$) galaxy during $z=0.1-2.5$. There may be an additional $\sim 0.5(0.3)$ major (minor) merger if we use the $H$-band flux ratio selection. The observed amount of major merging alone is sufficient to explain the observed number density evolution for the very massive (log$(M_{\star}/M_{\odot}) \geqslant 11.1$) galaxies. We argue that these very massive galaxies can put on a maximum of $6\%$ of stellar mass in addition to major and minor merging, so that their number density evolution remains consistent with observations. The observed number of major and minor mergers can increase the size of a massive quiescent galaxy by a factor of two at most. This amount of merging is enough to bring the compact quiescent galaxies formed at $z>2$ to lie at $1\sigma$ below the mean of the stellar mass-size relation as measured in some works (e.g. Newman et al. 2012), but additional mechanisms are needed to fully explain the evolution, and to be consistent with works suggesting stronger evolution (e.g. van der Wel et al. 2014).

Three billion years after the big bang (at redshift z=2), half of the most massive galaxies were already old, quiescent systems with little to no residual star formation and extremely compact with stellar mass densities at least an order of magnitude larger than in low redshift ellipticals, their descendants. Little is known about how they formed, but their evolved, dense stellar populations suggest formation within intense, compact starbursts 1-2 Gyr earlier (at 3<z<6). Simulations show that gas-rich major mergers can give rise to such starbursts which produce dense remnants. Sub-millimeter selected galaxies (SMGs) are prime examples of intense, gas-rich, starbursts. With a new, representative spectroscopic sample of compact quiescent galaxies at z=2 and a statistically well-understood sample of SMGs, we show that z=3-6 SMGs are consistent with being the progenitors of z=2 quiescent galaxies, matching their formation redshifts and their distributions of sizes, stellar masses and internal velocities. Assuming an evolutionary connection, their space densities also match if the mean duty cycle of SMG starbursts is 42 (+40/-29) Myr (consistent with independent estimates), which indicates that the bulk of stars in these massive galaxies were formed in a major, early surge of star-formation. These results suggests a coherent picture of the formation history of the most massive galaxies in the universe, from their initial burst of violent star-formation through their appearance as high stellar-density galaxy cores and to their ultimate fate as giant ellipticals.

Galaxy interactions/mergers have been shown to dominate the population of IR luminous galaxies (log(LIR)>11.6Lsun) in the local Universe (z<0.25). Recent studies based on the relation between galaxies' star formation rates and stellar mass (the SFR-M relation or the galaxy main sequence (MS)) have suggested that galaxy interaction/mergers may only become significant when galaxies fall well above the galaxy MS. Since the typical SFR at given M increases with redshift, the existence of galaxy MS implies that massive, IR-luminous galaxies at high-z may not necessarily be driven by galaxy interactions. We examine the role of galaxy interactions in the SFR-M relation by carrying out a morphological analysis of 2084 Herschel-selected galaxies at 0.2 < z < 1.5 in the COSMOS field. Herschel-PACS and -SPIRE observations covering the full 2-deg^2 COSMOS field provide one of the largest far-IR selected samples of high-redshift galaxies with well-determined redshifts to date, with sufficient sensitivity at z ~ 1, to sample objects lying on and above the galaxy MS. Using a detailed visual classification scheme, we show that the fraction of "disk galaxies" decreases and the fraction of "irregular" galaxies increases systematically with increasing LIR out to z ~ 1.5 and z ~ 1.0, respectively. At log(LIR) > 11.5 Lsun, >50% of the objects show evident features of strongly interacting/merger systems, where this percentage is similar to the studies of local IR-luminous galaxies. The fraction of interacting/merger systems also systematically increases with the deviation from the SFR-M relation, supporting the view that galaxies fall above the MS are more dominated by mergers than the MS galaxies. Meanwhile, we find that ~18% of massive IR-luminous MS galaxies are classified as interacting systems, where this population may not evolve through the evolutionary track predicted by a simple gas exhaustion model.

Quiescent massive galaxies at z~2 are thought to be the progenitors of present-day massive ellipticals. Observations revealed them to be extraordinarily compact. The determination of stellar ages, star formation rates and dust properties via spectroscopic measurements has up to now only been feasible for the most luminous and massive specimens (~3x M*). Here we present a spectroscopic study of two near-infrared selected galaxies which are close to the characteristic stellar mass M* (~0.9x M* and ~1.3x M*) and whose observed brightness has been boosted by the gravitational lensing effect. We measure the redshifts of the two galaxies to be z=1.71\pm0.02 and z=2.15\pm0.01. By fitting stellar population synthesis models to their spectro-photometric SEDs we determine their ages to be 2.4^{+0.8}_{-0.6} Gyr and 1.7\pm0.3 Gyr, respectively, which implies that the two galaxies have higher mass-to-light ratios than most quiescent z~2 galaxies in other studies. We find no direct evidence for active star-formation or AGN activity in either of the two galaxies, based on the non-detection of emission lines. Based on the derived redshifts and stellar ages we estimate the formation redshifts to be z=4.3^{+3.4}_{-1.2} and z=4.3^{+1.0}_{-0.6}, respectively. We use the increased spatial resolution due to the gravitational lensing to derive constraints on the morphology. Fitting Sersic profiles to the de-lensed images of the two galaxies confirms their compactness, with one of them being spheroid-like, and the other providing the first confirmation of a passive lenticular galaxy at a spectroscopically derived redshift z~2.

We present dynamical and structural scaling relations of quiescent galaxies at z=2, including the dynamical mass-size relation and the first constraints on the fundamental plane (FP). The backbone of the analysis is a new, very deep VLT/X-shooter spectrum of a massive, compact, quiescent galaxy at z=2.0389. We detect the continuum between 3700-22000A and several strong absorption features (Balmer series, Ca H+K, G-band), from which we derive a stellar velocity dispersion of 318 +/- 53 km/s. We perform detailed modeling of the continuum emission and line indices and derive strong simultaneous constraints on the age, metallicity, and stellar mass. The galaxy is a dusty (A_V=0.77 (+0.36,-0.32)) solar metallicity (log(Z/Zsun) = 0.02 (+0.20,-0.41)) post starburst galaxy, with a mean luminosity weighted log(age/yr) of 8.9 +/- 0.1. The galaxy formed the majority of its stars at z>3 and currently has little or no ongoing star formation. We compile a sample of three other z~2 quiescent galaxies with measured velocity dispersions, two of which are also post starburst like. Their dynamical mass-size relation is offset significantly less than the stellar mass-size relation from the local early type relations, which we attribute to a lower central dark matter fraction. Recent cosmological merger simulations qualitatively agree with the data, but can not fully account for the evolution in the dark matter fraction. The z~2 FP requires additional evolution beyond passive stellar aging, to be in agreement with the local FP. The structural evolution predicted by the cosmological simulations is insufficient, suggesting that additional, possibly non-homologous structural evolution is needed.

In this work I investigate the redshift evolution of pair fraction of a sample of 196 massive galaxies from z = 0 to 3, selected from the COSMOS field. We find that on average a massive galaxy undergoes ~ 1.1 \pm 0.5 major merger since z = 3. I will review the current limitations of using the pair fraction as a probe for quantifying the impact of mergers on galaxy evolution. This work is based on the paper Man et al. (2011).

Using a mass-selected ($M_{\star} \ge 10^{11} M_{\odot}$) sample of 198 galaxies at 0 < z < 3.0 with HST/NICMOS $H_{160}$-band images from the COSMOS survey, we find evidence for the evolution of the pair fraction above z ~ 2, an epoch in which massive galaxies are believed to undergo significant structural and mass evolution. We observe that the pair fraction of massive galaxies is 0.15 \pm 0.08 at 1.7 < z < 3.0, where galaxy pairs are defined as massive galaxies having a companion of flux ratio from 1:1 to 1:4 within a projected separation of 30 kpc. This is slightly lower, but still consistent with the pair fraction measured previously in other studies, and the merger fraction predicted in halo-occupation modelling. The redshift evolution of the pair fraction is described by a power law F(z) = (0.07 \pm 0.04) * (1+z) ^ (0.6 \pm 0.5). The merger rate is consistent with no redshift evolution, however it is difficult to constrain due to the limited sample size and the high uncertainties in the merging timescale. Based on the merger rate calculation, we estimate that a massive galaxy undergoes on average 1.1 \pm 0.5 major merger from z = 3 to 0. The observed merger fraction is sufficient to explain the number density evolution of massive galaxies, but insufficient to explain the size evolution. This is a hint that mechanism(s) other than major merging may be required to increase the sizes of the massive, compact quiescent galaxies from z ~ 2 to 0.