We present the [CII]158$\mu$m line luminosity functions (LFs) at $z\sim 4-6$ using the ALMA observations of 118 sources, which are selected to have UV luminosity $M_{1500A} < -20.2$ and optical spectroscopic redshifts in COSMOS and ECDF-S. Of the 118 targets, 75 have significant [C II] detections and 43 are upper limits. With this UV-selected sample, we set stringent lower bounds to the volume density of [C II] emitters at $z\sim4-6$. Our derived LFs are consistent with the $z\sim0$ [C II] LF at $\leq 10^9L_\odot$, but show a hint of an excess of [C II] emitters at the luminous end ($>10^9L_\odot$) compared to $z\sim0$. This excess is consistent with a potential population of [C II]-bright UV-faint sources that are largely excluded from the ALPINE target selection. To further constrain such sources, in a companion analysis, Loiacono et al. (2020) estimate the [C II] LF based on serendipitous [C II] sources in the ALPINE maps. Their resulting single data point is consistent with the excess we observe at high-L. Combined, these studies suggest possible redshift evolution in the [C II] LFs at high luminosity between $z\sim 5$ and $0$. This conclusion is further supported by the inferred [C II] LFs from the far-IR and CO LFs at $z\sim4-6$. Integrating our estimate on the best-fit [C II] LF, we set a constraint on molecular gas mass density at $z\sim4-6$, with $\rho_{mol} \sim (2-7)\times10^7M_\odot$ Mpc$^{-3}$. This is broadly consistent with previous studies, although the uncertainties span more than an order of magnitude. Finally, we find that available model predictions tend to underestimate the number densities of [C II] emitters at $z\sim5$.

The ALPINE-ALMA large program targets the [CII] 158 $\mu$m line and the far-infrared continuum in 118 spectroscopically confirmed star-forming galaxies between z=4.4 and z=5.9. It represents the first large [CII] statistical sample built in this redshift range. We present details of the data processing and the construction of the catalogs. We detected 23 of our targets in the continuum. To derive accurate infrared luminosities and obscured star formation rates, we measured the conversion factor from the ALMA 158 $\mu$m rest-frame dust continuum luminosity to the total infrared luminosity (L$_{\rm IR}$) after constraining the dust spectral energy distribution by stacking a photometric sample similar to ALPINE in ancillary single-dish far-infrared data. We found that our continuum detections have a median L$_{\rm IR}$ of 4.4$\times 10^{11}$ L$_\odot$. We also detected 57 additional continuum sources in our ALMA pointings. They are at lower redshift than the ALPINE targets, with a mean photometric redshift of 2.5$\pm$0.2. We measured the 850 $\mu$m number counts between 0.35 and 3.5 mJy, improving the current interferometric constraints in this flux density range. We found a slope break in the number counts around 3 mJy with a shallower slope below this value. More than 40 % of the cosmic infrared background is emitted by sources brighter than 0.35 mJy. Finally, we detected the [CII] line in 75 of our targets. Their median [CII] luminosity is 4.8$\times$10$^8$ L$_\odot$ and their median full width at half maximum is 252 km/s. After measuring the mean obscured SFR in various [CII] luminosity bins by stacking ALPINE continuum data, we find a good agreement between our data and the local and predicted SFR-L$_{\rm [CII]}$ relations of De Looze et al. (2014) and Lagache et al. (2018).

We present the first stellar velocity dispersion measurement of a massive quenching galaxy at z=4.01. The galaxy is first identified as a massive z>~4 galaxy with suppressed star formation from photometric redshifts based on deep multi-band data in the UKIDSS Ultra Deep Survey field. A follow-up spectroscopic observation with MOSFIRE on Keck revealed strong multiple absorption features, which are identified as Balmer absorption lines, giving a secure redshift of z=4.01. Thanks to the high S/N of the spectrum, we are able to estimate the stellar velocity dispersion, sigma=268+/-59 km/s. This velocity dispersion is consistent with that of massive galaxies today, implying no significant evolution in stellar velocity dispersion over the last 12 Gyr. Based on an upper limit on its physical size from deep optical images (r_eff<1.3 kpc), we find that its dynamical mass is consistent with the stellar mass inferred from photometry. Furthermore, the galaxy is located on the mass fundamental plane extrapolated from lower redshift galaxies. Combining all these results, we find that the velocity dispersion does not significantly evolve with redshift, although the size and mass of massive quenched galaxies do. This suggests that the mass in the core of massive galaxies does not evolve significantly, while most of the mass growth occurs in the outskirts of the galaxies, which also increases the size. This picture is consistent with a two-phase formation scenario in which mass and size growth is due to accretion in the outskirts of galaxies via mergers.

The [CII] 158 micron line is one of the strongest IR emission lines, which has been shown to trace the star-formation rate (SFR) of galaxies in the nearby Universe and up to $z \sim 2$. Whether this is also the case at higher redshift and in the early Universe remains debated. The ALPINE survey, which targeted 118 star-forming galaxies at $4.4 < z< 5.9$, provides a new opportunity to examine this question with the first statistical dataset. Using the ALPINE data and earlier measurements from the literature we examine the relation between the [CII] luminosity and the SFR over the entire redshift range from $z \sim 4-8$. ALPINE galaxies, which are both detected in [CII] and dust continuum, show a good agreement with the local L([CII])-SFR relation. Galaxies undetected in the continuum with ALMA are found to be over-luminous in [CII], when the UV SFR is used. After accounting for dust-obscured star formation, by an amount SFR(IR)$\approx$SFR(UV) on average, which results from two different stacking methods and SED fitting, the ALPINE galaxies show an L([CII])-SFR relation comparable to the local one. When [CII] non-detections are taken into account, the slope may be marginally steeper at high-z, although this is still somewhat uncertain. When compared in a homogeneous manner, the $z>6 $ [CII] measurements (detections and upper limits) do not behave very differently from the $z \sim 4-6$ data. We find a weak dependence of L([CII])/SFR on the Lyman-alpha equivalent width. Finally, we find that the ratio L([CII])/LIR $\sim (1-3) \times 10^{-3}$ for the ALPINE sources, comparable to that of "normal" galaxies at lower redshift. Our analysis, which includes the largest sample ($\sim 150$ galaxies) of [CII] measurements at $z>4$ available so far, suggests no or little evolution of the L([CII])-SFR relation over the last 13 Gyr of cosmic time.

We present the physical extent of [CII] 158um line-emitting gas from 46 star-forming galaxies at z=4-6 from the ALMA Large Program to INvestigate CII at Early Times (ALPINE). Using exponential profile fits, we measure the effective radius of the [CII] line (r_e,[CII]) for individual galaxies and compare them with the rest-frame ultra-violet (UV) continuum (r_e,UV) from Hubble Space Telescope images. The effective radius r_e,[CII] exceeds r_e,UV by factors of ~2-3 and the ratio of r_e,[CII]/r_e,UV increases as a function of M_star. We do not find strong evidence that [CII] line, the rest-frame UV, and FIR continuum are always displaced over ~ 1-kpc scale from each other. We identify 30% of isolated ALPINE sources as having an extended [CII] component over 10-kpc scales detected at 4.1$\sigma$-10.9$\sigma$ beyond the size of rest-frame UV and far-infrared (FIR) continuum. One object has tentative rotating features up to ~10-kpc, where the 3D model fit shows the rotating [CII]-gas disk spread over 4 times larger than the rest-frame UV-emitting region. Galaxies with the extended [CII] line structure have high star-formation rate (SFR), stellar mass (M_star), low Lya equivalent-width, and more blue-shifted (red-shifted) rest-frame UV metal absorption (Lya line), as compared to galaxies without such extended [CII] structures. Although we cannot rule out the possibility that a selection bias towards luminous objects may be responsible for such trends, the star-formation driven outflow also explains all these trends. Deeper observations are essential to test whether the extended [CII] line structures are ubiquitous to high-z star-forming galaxies.

Extensive surveys with the \textit{Hubble Space Telescope} (HST) over the past decade, targeting some of the most massive clusters in the sky, have uncovered dozens of galaxy-cluster strong lenses. The massive cluster strong-lens scale is typically $\theta_{E}\sim10\arcsec$ to $\sim30-35\arcsec$, with only a handful of clusters known with Einstein radii $\theta_{E}\sim40\arcsec$ or above (for $z_{source}=2$, nominally). Here we report another very large cluster lens, RXC J0032.1+1808 ($z=0.3956$), the second richest cluster in the redMapper cluster catalog and the 85th most massive cluster in the Planck Sunyaev-Zel'dovich catalog. With our Light-Traces-Mass and fully parametric (dPIEeNFW) approaches, we construct strong lensing models based on 18 multiple images of 5 background galaxies newly identified in the \textit{Hubble} data mainly from the \textit{Reionization Lensing Cluster Survey} (RELICS), in addition to a known sextuply imaged system in this cluster. Furthermore, we compare these models to Lenstool and GLAFIC models that were produced independently as part of the RELICS program. All models reveal a large effective Einstein radius of $\theta_{E}\simeq40\arcsec$ ($z_{source}=2$), owing to the obvious concentration of substructures near the cluster center. Although RXC J0032.1+1808 has a very large critical area and high lensing strength, only three magnified high-redshift candidates are found within the field targeted by RELICS. Nevertheless, we expect many more high-redshift candidates will be seen in wider and deeper observations with \textit{Hubble} or \emph{JWST}. Finally, the comparison between several algorithms demonstrates that the total error budget is largely dominated by systematic uncertainties.

We present dust attenuation properties of spectroscopically confirmed star forming galaxies on the main sequence at redshift ~4.4-5.8. Our analyses are based on the far infrared continuum observations of 118 galaxies at rest-frame $158\,\rm{\mu m}$ obtained with the ALMA large program ALPINE. We study the connection between the UV spectral slope ($\beta$), stellar mass ($M_{\star}$), and infrared excess (IRX$=L_{\rm{IR}}/L_{\rm{UV}}$). Twenty-three galaxies are individually detected in the continuum at $>3.5\,\sigma$ significance. We perform a stacking analysis using both detections and non-detections to study the average dust attenuation properties at z~4.4-5.8. The individual detections and stacks show that the IRX-$\beta$ relation at z~5 is consistent with a steeper dust attenuation curve than typically found at lower redshifts (z<4). The attenuation curve is similar to or even steeper than that of the extinction curve of the Small Magellanic Cloud (SMC). This systematic change of the IRX-$\beta$ relation as a function of redshift suggests an evolution of dust attenuation properties at $z>4$. Similarly, we find that our galaxies have lower IRX values up to 1 dex on average at fixed mass compared to previously studied IRX-$M_{\star}$ relations at $z\lesssim4$, albeit with significant scatter. This implies a lower obscured fraction of star-formation than at lower redshifts. Our results suggest that dust properties of UV-selected star forming galaxies at $z\gtrsim4$ are characterised by (i) a steeper attenuation curve than at $z\lesssim4$, and (ii) a rapidly decreasing dust obscured fraction of star formation as a function of redshift. Nevertheless, even among this UV-selected sample, massive galaxies ($\log M_{\star}/M_\odot > 10$) at z~5-6 already exhibit an obscured fraction of star formation of $\sim45\%$, indicating a rapid build-up of dust during the epoch of reionization.

The molecular gas content of normal galaxies at z>4 is poorly constrained, because the commonly used molecular gas tracers become hard to detect. We use the [CII]158um luminosity, recently proposed as a molecular gas tracer, to estimate the molecular gas content in a large sample of main-sequence star-forming galaxies at z=4.4-5.9, with a median stellar mass of 10^9.7 Msun, drawn from the ALMA Large Program to INvestigate [CII] at Early times (ALPINE) survey. The agreement between molecular gas masses derived from [CII] luminosity, dynamical mass, and rest-frame 850um luminosity, extrapolated from the rest-frame 158um continuum, supports [CII] as a reliable tracer of molecular gas in our sample. We find a continuous decline of the molecular gas depletion timescale from z=0 to z=5.9, which reaches a mean value of (4.6+/-0.8)x10^8 yr at z~5.5, only a factor of 2-3 shorter than in present-day galaxies. This suggests a mild enhancement of star formation efficiency toward high redshifts, unless the molecular gas fraction significantly increases. Our estimates show that the rise in molecular gas fraction as reported previously, flattens off above z~3.7 to achieve a mean value of 63%+/-3 over z=4.4-5.9. This redshift evolution of the gas fraction is in line with the one of the specific star formation rate. We use multi-epoch abundance matching to follow the gas fraction evolution over cosmic time of progenitors of z=0 Milky Way-like galaxies in 10^13 Msun halos and of more massive z=0 galaxies in 10^14 Msun halos. Interestingly, the former progenitors show a monotonic decrease of the gas fraction with cosmic time, while the latter show a constant gas fraction from z=5.9 to z~2 and a decrease at z<2. We discuss three possible effects, namely outflows, halt of gas supplying, and over-efficient star formation, which may jointly contribute to the gas fraction plateau of the latter massive galaxies.

We present the ancillary data and basic physical measurements for the galaxies in the ALMA Large Program to Investigate C+ at Early Times (ALPINE) survey - the first large multi-wavelength survey which aims at characterizing the gas and dust properties of 118 main-sequence galaxies at redshifts 4.4 < z < 5.9 via the measurement of [CII]-emission at 158 micro-meter (64% at >3.5$\sigma$) and the surrounding far-infrared (FIR) continuum in conjunction with a wealth of optical and near-infrared data. We outline in detail the spectroscopic data and selection of the galaxies as well as the ground- and space-based imaging products. In addition, we provide several basic measurements including stellar masses, star formation rates (SFR), rest-frame ultra-violet (UV) luminosities, UV continuum slopes ($\beta$), and absorption line redshifts, as well as H-alpha emission derived from Spitzer colors. We find that the ALPINE sample is representative of the 4 < z < 6 galaxy population selected by photometric methods and only slightly biased towards bluer colors ($\Delta\beta$ ~ 0.2). Using [CII] as tracer of the systemic redshift (confirmed for one galaxy at z=4.5 out of 118 for which we obtained optical [OII]$\lambda$3727A emission), we confirm red shifted Ly-alpha emission and blue shifted absorption lines similar to findings at lower redshifts. By stacking the rest-frame UV spectra in the [CII] rest-frame we find that the absorption lines in galaxies with high specific SFR are more blue shifted, which could be indicative of stronger winds and outflows.

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

Large photometric surveys provide a rich source of observations of quiescent galaxies, including a surprisingly large population at z>1. However, identifying large, but clean, samples of quiescent galaxies has proven difficult because of their near-degeneracy with interlopers such as dusty, star-forming galaxies. We describe a new technique for selecting quiescent galaxies based upon t-distributed stochastic neighbor embedding (t-SNE), an unsupervised machine learning algorithm for dimensionality reduction. This t-SNE selection provides an improvement both over UVJ, removing interlopers which otherwise would pass color selection, and over photometric template fitting, more strongly towards high redshift. Due to the similarity between the colors of high- and low-redshift quiescent galaxies, under our assumptions t-SNE outperforms template fitting in 63% of trials at redshifts where a large training sample already exists. It also may be able to select quiescent galaxies more efficiently at higher redshifts than the training sample.

We present new observations, carried out with IRAM NOEMA, of the atomic neutral carbon transitions [CI](1-0) at 492 GHz and [CI](2-1) at 809 GHz of GN20, a well-studied star-bursting galaxy at $z=4.05$. The high luminosity line ratio [CI](2-1)/[CI](1-0) implies an excitation temperature of $48^{+14}_{-9}$ K, which is significantly higher than the apparent dust temperature of $T_{\rm d}=33\pm2$ K ($\beta=1.9$) derived under the common assumption of an optically thin far-infrared dust emission, but fully consistent with $T_{\rm d}=52\pm5$ K of a general opacity model where the optical depth ($\tau$) reaches unity at a wavelength of $\lambda_0=170\pm23$ $\mu$m. Moreover, the general opacity solution returns a factor of $\sim 2\times$ lower dust mass and, hence, a lower molecular gas mass for a fixed gas-to-dust ratio, than with the optically thin dust model. The derived properties of GN20 thus provide an appealing solution to the puzzling discovery of starbursts appearing colder than main-sequence galaxies above $z>2.5$, in addition to a lower dust-to-stellar mass ratio that approaches the physical value predicted for starburst galaxies.

We report two secure ($z=3.775, 4.012$) and one tentative ($z\approx3.767$) spectroscopic confirmations of massive and quiescent galaxies through $K$-band observations with Keck/MOSFIRE and VLT/X-Shooter. The stellar continuum emission, the absence of strong nebular emission lines and the lack of significant far-infrared detections confirm the passive nature of these objects, disfavoring the alternative solution of low-redshift dusty star-forming interlopers. We derive stellar masses of $\mathrm{log}(M_{\star}/M_\odot)\sim11$ and ongoing star formation rates placing these galaxies $\gtrsim 1-2$ dex below the main sequence at their redshifts. The adopted parametrization of the star formation history suggests that these sources experienced a strong ($\langle \rm SFR \rangle \sim 1200-3500\,M_\odot\,\mathrm{yr}^{-1}$) and short ($\sim 50$ Myr) burst of star formation, peaking $\sim 150-500$ Myr before the time of observation, all properties reminiscent of the characteristics of sub-millimeter galaxies (SMGs) at $z>4$. We investigate this connection by comparing the comoving number densities and the properties of these two populations. We find a fair agreement only with the deepest sub-mm surveys detecting not only the most extreme starbursts, but also more normal galaxies. We support these findings by further exploring the Illustris-TNG cosmological simulation, retrieving populations of both fully quenched massive galaxies at $z\sim3-4$ and SMGs at $z\sim4-5$, with number densities and properties in agreement with the observations at $z\sim3$, but in increasing tension at higher redshift. Nevertheless, as suggested by the observations, not all the progenitors of quiescent galaxies at these redshifts shine as bright SMGs in their past and, similarly, not all bright SMGs quench by $z\sim3$, both fractions depending on the threshold assumed to define the SMGs themselves.

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 study the efficiency of galactic feedback in the early Universe by stacking the [C II] 158 um emission in a large sample of normal star-forming galaxies at 4 < z < 6 from the ALMA Large Program to INvestigate [C II] at Early times (ALPINE) survey. Searching for typical signatures of outflows in the high-velocity tails of the stacked [C II] profile, we observe (i) deviations from a single-component Gaussian model in the combined residuals and (ii) broad emission in the stacked [C II] spectrum, with velocities of |v|<~ 500 km/s. The significance of these features increases when stacking the subset of galaxies with star formation rates (SFRs) higher than the median (SFRmed = 25 Msun/yr), thus confirming their star-formation-driven nature. The estimated mass outflow rates are comparable to the SFRs, yielding mass-loading factors of the order of unity (similarly to local star-forming galaxies), suggesting that star-formation-driven feedback may play a lesser role in quenching galaxies at z > 4. From the stacking analysis of the datacubes, we find that the combined [C II] core emission (|v|< 200 km/s) of the higher-SFR galaxies is extended on physical sizes of ~ 30 kpc (diameter scale), well beyond the analogous [C II] core emission of lower-SFR galaxies and the stacked far-infrared continuum. The detection of such extended metal-enriched gas, likely tracing circumgalactic gas enriched by past outflows, corroborates previous similar studies, confirming that baryon cycle and gas exchanges with the circumgalactic medium are at work in normal star-forming galaxies already at early epochs.

We present ALMA Band 6 (nu=233GHz, lambda=1.3mm) continuum observations towards 68 'normal' star-forming galaxies within two Coma-like progenitor structures at z=2.10 and 2.47, from which ISM masses are derived, providing the largest census of molecular gas mass in overdense environments at these redshifts. Our sample comprises galaxies with a stellar mass range of 1x10^9M_sun - 4x10^11M_sun with a mean M_*~6x10^10M_sun. Combining these measurements with multiwavelength observations and SED modeling, we characterize the gas mass fraction and the star formation efficiency, and infer the impact of the environment on galaxies' evolution. Most of our detected galaxies (~70%) have star formation efficiencies and gas fractions similar to those found for coeval field galaxies and in agreement with the field scaling relations. However, we do find that the proto-clusters contain an increased fraction of massive, gas-poor galaxies, with low gas fractions (f_gas<6-10%) and red rest-frame ultraviolet/optical colors typical of post-starburst and passive galaxies. The relatively high abundance of passive galaxies suggests an accelerated evolution of massive galaxies in proto-cluster environments. The large fraction of quenched galaxies in these overdense structures also implies that environmental quenching takes place during the early phases of cluster assembly, even before virialization. From our data, we derive a quenching efficiency of E_q~0.45 and an upper limit on the quenching timescale of T_q<1Gyr.

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).

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.

Few detailed investigations of stellar populations in passive galaxies beyond z~1 are based on deep spectroscopic observations, due to the difficulty in obtaining such data. We present a study of stellar populations, structure, and mass-to-light ratios of a large sample of bulge-dominated galaxies in the two z=1.27 clusters Lynx E and Lynx W, based on deep ground-based optical spectroscopy combined with imaging from Hubble Space Telescope. We find that Lynx E has a well-defined core of red passive galaxies, while Lynx W lacks such a core. If all the sample galaxies evolve similarly in size from z=1.27 to the present, the data allow only 0.1 dex size-growth at a fixed dynamical mass. However, to link the Lynx central galaxies to brightest cluster galaxies similar to those of low redshift clusters, the Lynx galaxies would have to grow by at least a factor five, possibly through major merging. The mass-to-light ratios and the Balmer absorption lines of the Lynx galaxies are consistent with passive evolution of the stellar populations from z=1.27 to the present and support ages of 1--3 Gyr. The galaxies in the outskirts of the clusters contain younger stellar populations than found in the cluster cores. However, when evolved passively to z~0 both populations are consistent with the observed populations in the Coma cluster galaxies. The bulge-dominated emission line galaxies in the clusters are dominated by stellar populations with sub-solar metallicities. Thus, additional enrichment of these is required to produce Coma-like stellar populations by z~0.

We have obtained three-dimensional maps of the universe in $\sim200\times200\times80$ comoving Mpc$^3$ (cMpc$^3$) volumes each at $z=5.7$ and $6.6$ based on a spectroscopic sample of 179 galaxies that achieves $\gtrsim80$\% completeness down to the Ly$\alpha$ luminosity of $\log(L_{\rm Ly\alpha}/[\mathrm{erg\ s^{-1}}])=43.0$, based on our Keck and Gemini observations and the literature. The maps reveal filamentary large-scale structures and two remarkable overdensities made out of at least 44 and 12 galaxies at $z=5.692$ (z57OD) and $z=6.585$ (z66OD), respectively, making z66OD the most distant overdensity spectroscopically confirmed to date with $>10$ spectroscopically confirmed galaxies. We compare spatial distributions of submillimeter galaxies at $z\simeq 4-6$ with our $z=5.7$ galaxies forming the large-scale structures, and detect a $99.97\%$ signal of cross correlation, indicative of a clear coincidence of dusty star-forming galaxy and dust unobscured galaxy formation at this early epoch. The galaxies in z57OD and z66OD are actively forming stars with star formation rates (SFRs) $\gtrsim5$ times higher than the main sequence, and particularly the SFR density in z57OD is 10 times higher than the cosmic average at the redshift (a.k.a. the Madau-Lilly plot). Comparisons with numerical simulations suggest that z57OD and z66OD are protoclusters that are progenitors of the present-day clusters with halo masses of $\sim10^{14}\ \mathrm{M_\odot}$.

ALMA 870$\mu$m continuum imaging has uncovered a population of blends of multiple dusty star-forming galaxies (DSFGs) in sources originally detected with the Herschel Space Observatory. However, their pairwise separations are much smaller that what is found by ALMA follow-up of other single-dish surveys or expected from theoretical simulations. Using ALMA and VLA, we have targeted three of these systems to confirm whether the multiple 870$\mu$m continuum sources lie at the same redshift, successfully detecting $^{12}$CO($J = 3$-2) and $^{12}$CO($J = 1$-0) lines and being able to confirm that in the three cases all the multiple DSFGs are likely physically associated within the same structure. Therefore, we report the discovery of two new gas-rich dusty protocluster cores (HELAISS02, $z = 2.171 \pm 0.004$; HXMM20, $z = 2.602 \pm 0.002$). The third target is located in the well known COSMOS overdensity at $z = 2.51$ (named CL J1001+0220 in the literature), for which we do not find any new secure CO(1-0) detection, although some of its members show only tentative detections and require further confirmation. From the gas, dust, and stellar properties of the two new protocluster cores, we find very large molecular gas fractions yet low stellar masses, pushing the sources above the main sequence, while not enhancing their star formation efficiency. We suggest that the sources might be newly formed galaxies migrating to the main sequence. The properties of the three systems compared to each other and to field galaxies may suggest a different evolutionary stage between systems.

We combine high-resolution ALMA and HST/CANDELS observations of 20 submillimeter galaxies (SMGs) predominantly from the AS2UDS survey at z~2 with bright rest-frame optical counterparts (Ks < 22.9) to investigate the resolved structural properties of their dust and stellar components. We derive two-dimensional stellar-mass distributions that are inferred from spatial mass-to-light ratio (M/L) corrections based on rest-frame optical colors. Due to the high central column densities of dust in our SMGs, our mass distributions likely represent a lower limit to the true central mass density. The centroid positions between the inferred stellar-mass and the dust distributions agree within 1.1 kpc, indicating an overall good spatial agreement between the two components. The majority of our sources exhibit compact dust configurations relative to the stellar component (with a median ratio of effective radii Re,dust/Re,Mstar = 0.6). This ratio does not change with specific star-formation rate (sSFR) over the factor of 30 spanned by our targets, sampling the locus of "normal" main sequence galaxies up to the starburst regime, log(sSFR/sSFRMS) > 0.5. Our results imply that massive SMGs are experiencing centrally enhanced star formation unlike typical spiral galaxies in the local Universe. The sizes and stellar densities of our SMGs are in agreement with those of the passive population at z=1.5, consistent with these systems being the descendants of z~2 SMGs.

We present the first systematic study of the radio continuum size evolution of star-forming galaxies (SFGs) over the redshift range $0.35<z<2.25$. We use the VLA COSMOS 3GHz map (noise $\rm rms=2.3\,\mu Jy \,beam^{-1}$, $\theta_{\rm beam}=0.75\,\rm arcsec$) to construct a mass-complete sample of 3184 radio-selected SFGs that reside on and above the main-sequence (MS) of SFGs. We find no clear dependence between the radio size and stellar mass, $M_{\star}$, of SFGs with $10.5\lesssim\log(M_\star/\rm M_\odot)\lesssim11.5$. Our analysis suggests that MS galaxies are preferentially extended, while SFGs above the MS are always compact. The median effective radius of SFGs on (above) the MS of $R_{\rm eff}=1.5\pm0.2$ ($1.0\pm0.2$) kpc remains nearly constant with cosmic time; a parametrization of the form $R_{\rm eff}\propto(1+z)^\alpha$ yields a shallow slope of only $\alpha=-0.26\pm0.08\,(0.12\pm0.14)$ for SFGs on (above) the MS. The size of the stellar component of galaxies is larger than the extent of the radio continuum emission by a factor $\sim$2 (1.3) at $z=0.5\,(2)$, indicating star formation is enhanced at small radii. The galactic-averaged star formation rate surface density $(\Sigma_{\rm SFR})$ scales with the distance to the MS, except for a fraction of MS galaxies ($\lesssim10\%$) that harbor starburst-like $\Sigma_{\rm SFR}$. These "hidden" starbursts might have experienced a compaction phase due to disk instability and/or merger-driven burst of star formation, which may or may not significantly offset a galaxy from the MS. We thus propose to jointly use $\Sigma_{\rm SFR}$ and distance to the MS to better identify the galaxy population undergoing a starbursting phase.

Two decades of effort have been poured into both single-dish and interferometric millimeter-wave surveys of the sky to infer the volume density of dusty star-forming galaxies (DSFGs, with SFR>100M$_\odot$ yr$^{-1}$) over cosmic time. Though obscured galaxies dominate cosmic star-formation near its peak at $z\sim2$, the contribution of such heavily obscured galaxies to cosmic star-formation is unknown beyond $z\sim2.5$ in contrast to the well-studied population of Lyman-break galaxies (LBGs) studied through deep, space- and ground-based pencil beam surveys in the near-infrared. Unlocking the volume density of DSFGs beyond $z>3$, particularly within the first 1 Gyr after the Big Bang is critical to resolving key open questions about early Universe galaxy formation: (1) What is the integrated star-formation rate density of the Universe in the first few Gyr and how is it distributed among low-mass galaxies (e.g. Lyman-break galaxies) and high-mass galaxies (e.g. DSFGs and quasar host galaxies)? (2) How and where do the first massive galaxies assemble? (3) What can the most extreme DSFGs teach us about the mechanisms of dust production (e.g. supernovae, AGB stars, grain growth in the ISM) <1 Gyr after the Big Bang? We summarize the types of observations needed in the next decade to address these questions.

Strong gravitational lensing by clusters of galaxies probes the mass distribution at the core of each cluster and magnifies the universe behind it. MACS J0417.5-1154 at z=0.443 is one of the most massive clusters known based on weak lensing, X-ray, and Sunyaev-Zel'dovich analyses. Here we compute a strong lens model of MACS J0417 based on Hubble Space Telescope imaging observations collected, in part, by the Reionization Lensing Cluster Survey (RELICS), and recently reported spectroscopic redshifts from the MUSE instrument on the Very Large Telescope (VLT). We measure an Einstein radius of ThetaE=36'' at z = 9 and a mass projected within 200 kpc of M(200 kpc) = 1.78+0.01-0.03x10**14Msol. Using this model, we measure a ratio between the mass attributed to cluster-member galaxy halos and the main cluster halo of order 1:100. We assess the probability to detect magnified high-redshift galaxies in the field of this cluster, both for comparison with RELICS HST results and as a prediction for the James Webb Space Telescope (JWST) Guaranteed Time Observations upcoming for this cluster. Our lensing analysis indicates that this cluster has similar lensing strength to other clusters in the RELICS program. Our lensing analysis predicts a detection of at least a few z~6-8 galaxies behind this cluster, at odds with a recent analysis that yielded no such candidates in this field. Reliable strong lensing models are crucial for accurately predicting the intrinsic properties of lensed galaxies. As part of the RELICS program, our strong lensing model produced with the Lenstool parametric method is publicly available through the Mikulski Archive for Space Telescopes (MAST).

Large surveys of galaxy clusters with the Hubble and Spitzer Space Telescopes, including CLASH and the Frontier Fields, have demonstrated the power of strong gravitational lensing to efficiently deliver large samples of high-redshift galaxies. We extend this strategy through a wider, shallower survey named RELICS, the Reionization Lensing Cluster Survey. This survey, described here, was designed primarily to deliver the best and brightest high-redshift candidates from the first billion years after the Big Bang. RELICS observed 41 massive galaxy clusters with Hubble and Spitzer at 0.4-1.7um and 3.0-5.0um, respectively. We selected 21 clusters based on Planck PSZ2 mass estimates and the other 20 based on observed or inferred lensing strength. Our 188-orbit Hubble Treasury Program obtained the first high-resolution near-infrared images of these clusters to efficiently search for lensed high-redshift galaxies. We observed 46 WFC3/IR pointings (~200 arcmin^2) with two orbits divided among four filters (F105W, F125W, F140W, and F160W) and ACS imaging as needed to achieve single-orbit depth in each of three filters (F435W, F606W, and F814W). As previously reported by Salmon et al., we discovered 322 z ~ 6 - 10 candidates, including the brightest known at z ~ 6, and the most distant spatially-resolved lensed arc known at z ~ 10. Spitzer IRAC imaging (945 hours awarded, plus 100 archival) has crucially enabled us to distinguish z ~ 10 candidates from z ~ 2 interlopers. For each cluster, two HST observing epochs were staggered by about a month, enabling us to discover 11 supernovae, including 3 lensed supernovae, which we followed up with 20 orbits from our program. We delivered reduced HST images and catalogs of all clusters to the public via MAST and reduced Spitzer images via IRSA. We have also begun delivering lens models of all clusters, to be completed before the JWST GO call for proposals.

Strong gravitational lensing (SL) is a powerful means to map the distribution of dark matter. In this work, we perform a SL analysis of the prominent X-ray cluster RXJ0152.7-1357 (z=0.83, also known as CL 0152.7-1357) in \textit{Hubble Space Telescope} images, taken in the framework of the Reionization Lensing Cluster Survey (RELICS). On top of a previously known $z=3.93$ galaxy multiply imaged by RXJ0152.7-1357, for which we identify an additional multiple image, guided by a light-traces-mass approach we identify seven new sets of multiply imaged background sources lensed by this cluster, spanning the redshift range [1.79-3.93]. A total of 25 multiple images are seen over a small area of ~0.4 $arcmin^2$, allowing us to put relatively high-resolution constraints on the inner matter distribution. Although modestly massive, the high degree of substructure together with its very elongated shape make RXJ0152.7-1357 a very efficient lens for its size. This cluster also comprises the third-largest sample of z~6-7 candidates in the RELICS survey. Finally, we present a comparison of our resulting mass distribution and magnification estimates with those from a Lenstool model. These models are made publicly available through the MAST archive.

We present the results of a search for galaxy clusters and groups in the $\sim2$ square degree of the COSMOS field using all available X-ray observations from the XMM-Newton and Chandra observatories. We reach an X-ray flux limit of $3\times10^{-16}\;ergs\;cm^{-2}\;s^{-1}$ in 0.5--2 keV range, and identify 247 X-ray groups with $M_{200c}=8\times10^{12}-3\times10^{14}\;M_{\odot}$ at a redshift range of $0.08\leq z<1.53$, using the multiband photometric redshift and the master spectroscopic redshift catalogues of the COSMOS. The X-ray centres of groups are determined using high-resolution Chandra imaging. We investigate the relations between the offset of the brightest group galaxies (BGGs) from halo X-ray centre and group properties and compare with predictions from semi-analytic models and hydrodynamical simulations. We find that BGG offset decreases with both increasing halo mass and decreasing redshift with no strong dependence on the X-ray flux and SNR. We show that the BGG offset decreases as a function of increasing magnitude gap with no considerable redshift dependent trend. The stellar mass of BGGs in observations extends over a wider dynamic range compared to model predictions. At $z<0.5$, the central dominant BGGs become more massive than those with large offsets by up to 0.3dex, in agreement with model prediction. The observed and predicted lognormal scatter in the stellar mass of both low- and large-offset BGGs at fixed halo mass is $\sim0.3$dex.

[Abridged] We combine new CO(1-0) line observations of 24 intermediate redshift galaxies (0.03 < z < 0.28) along with literature data of galaxies at 0<z<4 to explore scaling relations between the dust and gas content using PAH 6.2 $\mu$m ($L_{6.2}$), CO ($L'_{\rm CO}$), and infrared ($L_{\rm IR}$) luminosities for a wide range of redshifts and physical environments. Our analysis confirms the existence of a universal $L_{6.2}-L_{\rm CO}$ correlation followed by normal star-forming galaxies (SFGs) and starbursts (SBs) at all redshifts. This relation is also followed by local ULIRGs that appear as outliers in the $L_{6.2}-L_{\rm IR}$ and $L_{\rm IR}-L'_{\rm CO}$ relations from the sequence defined by normal SFGs. The emerging tight ($\sigma \approx 0.26$ dex) and linear ($\alpha = 1.03$) relation between $L_{6.2}$ and $L_{\rm CO}$ indicates a $L_{6.2}$ to molecular gas ($M_{\rm H_2}$) conversion factor of $\alpha_{6.2} = M_{\rm H2}/L_{6.2} = (2.7\pm1.3) \times \alpha_{\rm CO}$, where $\alpha_{\rm CO}$ is the $L'_{\rm CO}$ to $M_{\rm H_2}$ conversion factor. We also find that on galaxy integrated scales, PAH emission is better correlated with cold rather than with warm dust emission, suggesting that PAHs are associated with the diffuse cold dust, which is another proxy for $M_{\rm H_2}$. Focusing on normal SFGs among our sample, we employ the dust continuum emission to derive $M_{\rm H_2}$ estimates and find a constant $M_{\rm H_2}/L_{6.2}$ ratio of $\alpha_{6.2} = 12.3 \ M_{\rm H_2}/{\rm L}_{\odot}$ ($\sigma\approx 0.3$ dex). We propose that the presented $L_{6.2}-L'_{\rm CO}$ and $L_{6.2}-M_{\rm H_2}$ relations will serve as useful tools for the determination of the physical properties of high-$z$ SFGs, for which PAH emission will be routinely detected by the James Webb Space Telescope.

We present the rest-frame optical sizes of massive quiescent galaxies (QGs) at $z\sim4$ measured at $K'$-band with the Infrared Camera and Spectrograph (IRCS) and AO188 on the Subaru telescope. Based on a deep multi-wavelength catalog in the Subaru XMM-Newton Deep Survey Field (SXDS), covering a wide wavelength range from the $u$-band to the IRAC $8.0\mu m$ over 0.7 deg$^2$, we evaluate photometric redshift to identify massive ($M_{\star}\sim10^{11}\ M_\odot$) galaxies with suppressed star formation. These galaxies show a prominent 4000$\rm \AA$ break feature at $z\sim4$, suggestive of an evolved stellar population. We then conduct follow-up $K'$-band imaging with adaptive optics for the five brightest galaxies ($K_{AB,total}=22.5\sim23.4$). Compared to lower redshift ones, QGs at $z\sim4$ have smaller physical sizes of effective radii $r_{eff}=0.2$ to $1.8$ kpc. The mean size measured by stacking the four brightest objects is $r_{eff}=0.7\rm\ kpc$. This is the first measurement of the rest-frame optical sizes of QGs at $z\sim4$. We evaluate the robustness of our size measurements using simulations and find that our size estimates are reasonably accurate with an expected systematic bias of $\sim0.2$ kpc. If we account for the stellar mass evolution, massive QGs at $z\sim4$ are likely to evolve into the most massive galaxies today. We find their size evolution with cosmic time in a form of $\log(r_e/{\rm kpc})= -0.44+1.77 \log(t/\rm Gyr)$. Their size growth is proportional to the square of stellar mass, indicating the size-stellar mass growth driven by minor dry mergers.

We present a strong-lensing analysis of four massive galaxy clusters imaged with the Hubble Space Telescope in the Reionization Lensing Cluster Survey. We use a Light-Traces-Mass technique to uncover sets of multiply images and constrain the mass distribution of the clusters. These mass models are the first published for Abell S295 and MACS J0159.8-0849, and are improvements over previous models for Abell 697 and MACS J0025.4-1222. Our analysis for MACS J0025.4-1222 and Abell S295 shows a bimodal mass distribution supporting the merger scenarios proposed for these clusters. The updated model for MACS J0025.4-1222 suggests a substantially smaller critical area than previously estimated. For MACS J0159.8-0849 and Abell 697 we find a single peak and relatively regular morphology, revealing fairly relaxed clusters. Despite being less prominent lenses, three of these clusters seem to have lensing strengths, i.e. cumulative area above certain magnification, similar to the Hubble Frontier Fields clusters (e.g., A($\mu>5$) $\sim 1-3$ arcmin$^2$, A($\mu>10$) $\sim 0.5-1.5$ arcmin$^2$), which in part can be attributed to their merging configurations. We make our lens models publicly available through the Mikulski Archive for Space Telescopes. Finally, using Gemini-N/GMOS spectroscopic observations we detect a single emission line from a high-redshift $J_{125}\simeq25.7$ galaxy candidate lensed by Abell 697. While we cannot rule out a lower-redshift solution, we interpret the line as Ly$\alpha$ at $z=5.800\pm 0.001$, in agreement with its photometric redshift and dropout nature. Within this scenario we measure a Ly$\alpha$ rest-frame equivalent width of $52\pm22$ \AA, and an observed Gaussian width of $117\pm 15$ km/s.

We present morphological properties of dusty star-forming galaxies at z=1-3 determined with high-resolution (FWHM~0"19) Atacama Large Milllimeter/submilimeter Array (ALMA) 1-mm band maps of our ASAGAO survey covering a 26-arcmin^2 area in GOODS-S. In conjunction with the ALMA archival data, the present sample consists of 42 ALMA sources with a wide rest-frame far-infrared (FIR) luminosity L_FIR range of ~10^11-10^13 Lo. To obtain an average rest-frame FIR profile, we perform individual measurements and careful stacking of the ALMA sources using the uv-visibility method that includes positional-uncertainty and smoothing-effect evaluations through Monte-Carlo simulations. We find that the dusty star-forming galaxies have the average FIR-wavelength Sersic index and effective radius of n_FIR=1.2+/-0.2 and R_e,FIR=1.0-1.3 kpc, respectively, additionally with a point source at the center, indicative of the existence of AGN. The average FIR profile agrees with a morphology of an exponential-disk clearly distinguished from a spheroidal profile (Sersic index of 4). We also examine the rest-frame optical Sersic index n_opt and effective radius R_e,opt with the deep Hubble Space Telescope (HST) images. Interestingly, we obtain n_opt=0.9+/-0.3 (~n_FIR) and R_e,opt=3.2+/-0.6 kpc (>R_e,FIR), suggesting that the FIR-emitting disk is embedded within a larger stellar disk. The rest-frame UV and FIR data of HST and ALMA provide us a radial surface density profile of the total star-formation rate (SFR), where the FIR SFR dominates over the UV SFR at the center. Under the simple assumption of a constant SFR, a compact stellar distribution found in z~1-2 compact quiescent galaxies (cQGs) is well reproduced, while a spheroidal stellar morphology of cQGs (n_opt=4) cannot, suggestive of other important mechanisms such as dynamical dissipation.

We investigate Lya, [OIII]5007, Ha, and [CII]158um emission from 1124 galaxies at z=4.9-7.0. Our sample is composed of 1092 Lya emitters (LAEs) at z=4.9, 5.7, 6.6, and 7.0 identified by Subaru/Hyper Suprime-Cam (HSC) narrowband surveys covered by Spitzer large area survey with Subaru/HSC (SPLASH) and 34 galaxies at z=5.148-7.508 with deep ALMA [CII]158um data in the literature. Fluxes of strong rest-frame optical lines of [OIII] and Ha (Hb) are constrained by significant excesses found in the SPLASH 3.6 and 4.5um photometry. At z=4.9, we find that the rest-frame Ha equivalent width and the Lya escape fraction f_Lya positively correlate with the rest-frame Lya equivalent width EW^0_Lya. The f_Lya-EW^0_Lya correlation is similarly found at z~0-2, suggesting no evolution of the correlation over z~0-5. The typical ionizing photon production efficiency of LAEs is logxi_ion/[Hz erg^-1]~25.5 significantly (60-100%) higher than those of LBGs at a given UV magnitude. At z=5.7-7.0, there exists an interesting turn-over trend that the [OIII]/Ha flux ratio increases in EW^0_Lya~0-30 A, and then decreases out to EW^0_Lya~130 A. We also identify an anti-correlation between a [CII] luminosity to star-formation rate ratio (L_[CII]/SFR) and EW^0_Lya at the >99% confidence level. We carefully investigate physical origins of the correlations with stellar-synthesis and photoionization models, and find that a simple anti-correlation between EW_Lya^0 and metallicity explains self-consistently all of the correlations of Lya, Ha, [OIII]/Ha, and [CII] identified in our study, indicating detections of metal-poor (~0.03 Zo) galaxies with EW^0_Lya~200 A.

Strong gravitational lensing by galaxy clusters has become a powerful tool for probing the high-redshift Universe, magnifying distant and faint background galaxies. Reliable strong lensing (SL) models are crucial for determining the intrinsic properties of distant, magnified sources and for constructing their luminosity function. We present here the first SL analysis of MACS J0308.9+2645 and PLCK G171.9-40.7, two massive galaxy clusters imaged with the Hubble Space Telescope in the framework of the Reionization Lensing Cluster Survey (RELICS). We use the Light-Traces-Mass modeling technique to uncover sets of multiply imaged galaxies and constrain the mass distribution of the clusters. Our SL analysis reveals that both clusters have particularly large Einstein radii ($\theta_E>30"$ for a source redshift of $z_s=2$), providing fairly large areas with high magnifications, useful for high-redshift galaxy searches ($\sim2$ arcmin$^{2}$ with $\mu>5$ to $\sim1$ arcmin$^{2}$ with $\mu>10$, similar to a typical \textit{Hubble Frontier Fields} cluster). We also find that MACS J0308.9+2645 hosts a promising, apparently bright (J$\sim23.2-24.6$ AB), multiply imaged high-redshift candidate at $z\sim6.4$. These images are amongst the brightest high-redshift candidates found in RELICS. Our mass models, including magnification maps, are made publicly available for the community through the Mikulski Archive for Space Telescopes.

Dust-enshrouded, starbursting, submillimeter galaxies (SMGs) at $z \geq 3$ have been proposed as progenitors of $z \geq 2$ compact quiescent galaxies (cQGs). To test this connection, we present a detailed spatially resolved study of the stars, dust and stellar mass in a sample of six submillimeter-bright starburst galaxies at $z \sim 4.5$. The stellar UV emission probed by HST is extended, irregular and shows evidence of multiple components. Informed by HST, we deblend Spitzer/IRAC data at rest-frame optical finding that the systems are undergoing minor mergers, with a typical stellar mass ratio of 1:6.5. The FIR dust continuum emission traced by ALMA locates the bulk of star formation in extremely compact regions (median $r_{\rm{e}} = 0.70 \pm 0.29$ kpc) and it is in all cases associated with the most massive component of the mergers (median $\log (M_{*}/M_{\odot}) = 10.49 \pm 0.32$). We compare spatially resolved UV slope ($\beta$) maps with the FIR dust continuum to study the infrared excess ($\rm{IRX} = L_{\rm{IR}}/L_{\rm{UV}}$)-$\beta$ relation. The SMGs display systematically higher $\rm{IRX}$ values than expected from the nominal trend, demonstrating that the FIR and UV emissions are spatially disconnected. Finally, we show that the SMGs fall on the mass-size plane at smaller stellar masses and sizes than cQGs at $z = 2$. Taking into account the expected evolution in stellar mass and size between $z = 4.5$ and $z = 2$ due to the ongoing starburst and mergers with minor companions, this is in agreement with a direct evolutionary connection between the two populations.

We studied the molecular gas properties of AzTEC/C159, a star-forming disk galaxy at $z=4.567$. We secured $^{12}$CO molecular line detections for the $J=2\to1$ and $J=5\to4$ transitions using the Karl G. Jansky VLA and the NOEMA interferometer. The broad (FWHM$\sim750\,{\rm km\,s}^{-1}$) and tentative double-peaked profiles of both $^{12}$CO lines are consistent with an extended molecular gas reservoir, which is distributed in a rotating disk as previously revealed from [CII] 158$\mu$m line observations. Based on the $^{12}$CO(2$\to$1) emission line we derived $L'_{\rm{CO}}=(3.4\pm0.6)\times10^{10}{\rm \,K\,km\,s}^{-1}{\rm \,pc}^{2}$, that yields a molecular gas mass of $M_{\rm H_2 }(\alpha_{\rm CO}/4.3)=(1.5\pm0.3)\times 10^{11}{\rm M}_\odot$ and unveils a gas-rich system with $\mu_{\rm gas}(\alpha_{\rm CO}/4.3)\equiv M_{\rm H_2}/M_\star=3.3\pm0.7$. The extreme star formation efficiency (SFE) of AzTEC/C159, parametrized by the ratio $L_{\rm{IR}}/L'_{\rm{CO}}=(216\pm80)\, {\rm L}_{\odot}{\rm \,(K\,km\,s}^{-1}{\rm \,pc}^{2})^{-1}$, is comparable to merger-driven starbursts such as local ultra-luminous infrared galaxies (ULIRGs) and SMGs. Likewise, the $^{12}$CO(5$\to$4)/CO(2$\to$1) line brightness temperature ratio of $r_{52}= 0.55\pm 0.15$ is consistent with high excitation conditions, similar to that observed in SMGs. We constrained the value for the $L'_{\text{CO}}-{\rm H}_2$ mass conversion factor in AzTEC/C159, i.e. $\alpha_{\text{CO}}=3.9^{+2.7}_{-1.3}{\rm \,M}_{\odot}{\rm \,K}^{-1}{\rm \,km}^{-1}{\rm \,s\,pc}^{-2}$, that is consistent with a self-gravitating molecular gas distribution as observed in local star-forming disk galaxies. Cold gas streams from cosmological filaments might be fueling a gravitationally unstable gas-rich disk in AzTEC/C159, which breaks into giant clumps forming stars as efficiently as in merger-driven systems and generate high gas excitation.

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.

Recent literature suggests that there are two modes through which galaxies grow their stellar mass - a normal mode characterized by quasi-steady star formation, and a highly efficient starburst mode possibly triggered by stochastic events such as galaxy mergers. While these differences are established for extreme cases, the population of galaxies in-between these two regimes is poorly studied and it is not clear where the transition between these two modes of star formation occurs. We utilize ALMA observations of the CO J=3-2 line luminosity in a sample of 20 infrared luminous galaxies that lie in the intermediate range between normal and starburst galaxies at z ~ 0.25-0.6 in the COSMOS field to examine the gas content and star formation efficiency of these galaxies. We compare these quantities to the galaxies' deviation from the well-studied "main sequence" correlation between star formation rate and stellar mass (MS) and find that at log($SFR/SFR_{MS}$) < 0.6, a galaxy's distance to the main sequence is mostly driven by increased gas content, and not a more efficient star formation process.

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 measure the stellar mass function (SMF) of galaxies in the COSMOS field up to $z\sim6$. We select them in the near-IR bands of the COSMOS2015 catalogue, which includes ultra-deep photometry from UltraVISTA-DR2, SPLASH, and Subaru/Hyper-SuprimeCam. At $z>2.5$ we use new precise photometric redshifts with error $\sigma_z=0.03(1+z)$ and an outlier fraction of $12\%$, estimated by means of the unique spectroscopic sample of COSMOS. The increased exposure time in the DR2, along with our panchromatic detection strategy, allow us to improve the stellar mass completeness at high $z$ with respect to previous UltraVISTA catalogues. We also identify passive galaxies through a robust colour-colour selection, extending their SMF estimate up to $z=4$. Our work provides a comprehensive view of galaxy stellar mass assembly between $z=0.1$ and 6, for the first time using consistent estimates across the entire redshift range. We fit these measurements with a Schechter function, correcting for Eddington bias. We compare the SMF fit with the halo mass function predicted from $\Lambda$CDM simulations. We find that at $z>3$ both functions decline with a similar slope in the high-mass end. This feature could be explained assuming that the mechanisms that quench star formation in massive haloes become less effective at high redshift; however further work needs to be done to confirm this scenario. Concerning the SMF low-mass end, it shows a progressive steepening as moving towards higher redshifts, with $\alpha$ decreasing from $-1.47_{-0.02}^{+0.02}$ at $z\simeq0.1$ to $-2.11_{-0.13}^{+0.30}$ at $z\simeq5$. This slope depends on the characterisation of the observational uncertainties, which is crucial to properly remove the Eddington bias. We show that there is currently no consensus on the method to quantify such errors: different error models result in different best-fit Schechter parameters. [Abridged]

We present millimetre dust emission measurements of two Lyman Break Galaxies at z~3 and construct for the first time fully sampled infrared spectral energy distributions (SEDs), from mid-IR to the Rayleigh-Jeans tail, of individually detected, unlensed, UV-selected, main sequence (MS) galaxies at $z=3$. The SED modelling of the two sources confirms previous findings, based on stacked ensembles, of an increasing mean radiation field <U> with redshift, consistent with a rapidly decreasing gas metallicity in z > 2 galaxies. Complementing our study with CO[3-2] emission line observations, we measure the molecular gas mass (M_H2) reservoir of the systems using three independent approaches: 1) CO line observations, 2) the dust to gas mass ratio vs metallicity relation and 3) a single band, dust emission flux on the Rayleigh-Jeans side of the SED. All techniques return consistent M_H2 estimates within a factor of ~2 or less, yielding gas depletion time-scales (tau_dep ~ 0.35 Gyrs) and gas-to-stellar mass ratios (M_H2/M* ~ 0.5-1) for our z~3 massive MS galaxies. The overall properties of our galaxies are consistent with trends and relations established at lower redshifts, extending the apparent uniformity of star-forming galaxies over the last 11.5 billion years.

ALMA observations of the long wavelength dust continuum are used to estimate the interstellar medium (ISM) masses in a sample of 708 galaxies at z = 0.3 to 4.5 in the COSMOS field. The galaxy sample has known far-infrared luminosities and, hence, star formation rates (SFRs), and stellar masses (M$_{\rm *}$) from the optical-infrared spectrum fitting. The galaxies sample SFRs from the main sequence (MS) to 50 times above the MS. The derived ISM masses are used to determine the dependence of gas mass on redshift, M$_{\rm *}$, and specific SFR (sSFR) relative to the MS. The ISM masses increase approximately 0.63 power of the rate of increase in SFRs with redshift and the 0.32 power of the sSFR/sSFR$_MS$. The SF efficiencies also increase as the 0.36 power of the SFR redshift evolutionary and the 0.7 power of the elevation above the MS; thus the increased activities at early epochs are driven by both increased ISM masses and SF efficiency. Using the derived ISM mass function we estimate the accretion rates of gas required to maintain continuity of the MS evolution ($>100$ \msun yr$^{-1}$ at z $>$ 2.5). Simple power-law dependences are similarly derived for the gas accretion rates. We argue that the overall evolution of galaxies is driven by the rates of gas accretion. The cosmic evolution of total ISM mass is estimated and linked to the evolution of SF and AGN activity at early epochs.

We present an overview of the "KIFF" project, which provides ultra-deep Ks-band imaging of all six of the Hubble Frontier Fields clusters Abell 2744, MACS-0416, Abell S1063, Abell 370, MACS-0717 and MACS-1149. All of these fields have recently been observed with large allocations of Directors' Discretionary Time with the HST and Spitzer telescopes covering 0.4 < lambda < 1.6 microns and 3.6--4.5 microns, respectively. VLT/HAWK-I integrations of the first four fields reach 5-sigma limiting depths of Ks~26.0 (AB, point sources) and have excellent image quality (FWHM ~ 0."4). Shorter Keck/MOSFIRE integrations of the MACS-0717 (MACS-1149) field better observable in the north reach limiting depths Ks=25.5 (25.1) with seeing FWHM ~0."4 (0."5). In all cases the Ks-band mosaics cover the primary cluster and parallel HST/ACS+WFC3 fields. The total area of the Ks-band coverage is 490 arcmin^2. The Ks-band at 2.2 microns crucially fills the gap between the reddest HST filter (1.6 micron ~ H-band) and the IRAC 3.6 micron passband. While reaching the full depths of the space-based imaging is not currently feasible from the ground, the deep Ks-band images provide important constraints on both the redshifts and the stellar population properties of galaxies extending well below the characteristic stellar mass across most of the age of the universe, down to, and including, the redshifts of the targeted galaxy clusters (z < 0.5).

ALMA Cycle 2 observations of the long wavelength dust emission in 145 star-forming galaxies are used to probe the evolution of star-forming ISM. We also develop the physical basis and empirical calibration (with 72 low-z and z ~ 2 galaxies) for using the dust continuum as a quantitative probe of interstellar medium (ISM) masses. The galaxies with highest star formation rates (SFRs) at <z> = 2.2 and 4.4 have gas masses up to 100 times that of the Milky Way and gas mass fractions reaching 50 to 80%, i.e. gas masses 1 - 4 times their stellar masses. We find a single high-z star formation law: SFR = 35 M_ mol^0.89 x (1+z)_{z=2}^0.95 x (sSFR)_{MS}^0.23 \msun yr^-1 -- an approximately linear dependence on the ISM mass and an increased star formation efficiency per unit gas mass at higher redshift. Galaxies above the Main Sequence (MS) have larger gas masses but are converting their ISM into stars on a timescale only slightly shorter than those on the MS -- thus these 'starbursts' are largely the result of having greatly increased gas masses rather than and increased efficiency for converting gas to stars. At z $> 1$, the entire population of star-forming galaxies has $\sim$ 2 - 5 times shorter gas depletion times than low-z galaxies. These shorter depletion times indicate a different mode of star formation in the early universe -- most likely dynamically driven by compressive, high-dispersion gas motions -- a natural consequence of the high gas accretion rates.

We present the COSMOS2015 catalog which contains precise photometric redshifts and stellar masses for more than half a million objects over the 2deg$^{2}$ COSMOS field. Including new $YJHK_{\rm s}$ images from the UltraVISTA-DR2 survey, $Y$-band from Subaru/Hyper-Suprime-Cam and infrared data from the Spitzer Large Area Survey with the Hyper-Suprime-Cam Spitzer legacy program, this near-infrared-selected catalog is highly optimized for the study of galaxy evolution and environments in the early Universe. To maximise catalog completeness for bluer objects and at higher redshifts, objects have been detected on a $\chi^{2}$ sum of the $YJHK_{\rm s}$ and $z^{++}$ images. The catalog contains $\sim 6\times 10^5$ objects in the 1.5 deg$^{2}$ UltraVISTA-DR2 region, and $\sim 1.5\times 10^5$ objects are detected in the "ultra-deep stripes" (0.62 deg$^{2}$) at $K_{\rm s}\leq 24.7$ (3$\sigma$, 3", AB magnitude). Through a comparison with the zCOSMOS-bright spectroscopic redshifts, we measure a photometric redshift precision of $\sigma_{\Delta z/(1+z_s)}$ = 0.007 and a catastrophic failure fraction of $\eta=0.5$%. At $3<z<6$, using the unique database of spectroscopic redshifts in COSMOS, we find $\sigma_{\Delta z/(1+z_s)}$ = 0.021 and $\eta=13.2\% $. The deepest regions reach a 90\% completeness limit of 10$^{10}M_\odot$ to $z=4$. Detailed comparisons of the color distributions, number counts, and clustering show excellent agreement with the literature in the same mass ranges. COSMOS2015 represents a unique, publicly available, valuable resource with which to investigate the evolution of galaxies within their environment back to the earliest stages of the history of the Universe. The COSMOS2015 catalog is distributed via anonymous ftp (ftp://ftp.iap.fr/pub/from_users/hjmcc/COSMOS2015/) and through the usual astronomical archive systems (CDS, ESO Phase 3, IRSA).

We study the environment of 23 submillimetre galaxies (SMGs) drawn from the JCMT/AzTEC 1.1mm S/N-limited sample in the COSMOS field, as well as 4 COSMOS SMGs at z_spec>4.5, and 1 at z_spec=2.49, yielding a sample of 28 SMGs. We search for overdensities using the COSMOS photometric redshifts based on over 30 UV-NIR photometric bands, reaching an accuracy of sigma(Delta z/(1+z))=0.0067 (0.0155) at z<3.5 (>3.5). To identify overdensities we apply the Voronoi tessellation analysis, and estimate the overdensity estimator delta_g as a function of distance from the SMG and/or overdensity center. We test and validate our approach via simulations, X-ray detected groups, and spectroscopic verifications using VUDS and zCOSMOS catalogues showing that even with photometric redshifts in COSMOS we can efficiently retrieve overdensities out to z~5. Our results yield that 11/23 (48%) JCMT/AzTEC 1.1mm SMGs occupy overdense environments. Considering the entire JCMT/AzTEC 1.1mm S/N>4 sample, and accounting for the expected fraction of spurious detections, yields that 35-61% of the SMGs in the S/N-limited sample occupy overdense environments. We perform an X-ray stacking analysis in the 0.5-2keV band using a 32" aperture and our SMG positions, and find statistically significant detections. For our z<2 [z>2] subsample we find an average flux of (4.0+/-0.8)x10^{-16} [(1.3+/-0.5)x10^{-16}] erg/s/cm^2 and a corresponding total mass of M_200 = 2.8x10^{13} [2x10^{13}] MSol. Our results suggest a higher occurrence of SMGs occupying overdense environments at z>3, than at z<3. This may be understood if highly star forming galaxies can only be formed in the highest peaks of the density field tracing the most massive dark matter haloes at early cosmic epochs, while at later times cosmic structure may have matured sufficiently that more modest overdensities correspond to sufficiently massive haloes to form SMGs.

ALMA Cycle 2 observations of the long wavelength dust emission in 180 star-forming (SF) galaxies are used to investigate the evolution of ISM masses at z = 1 to 6.4. The ISM masses exhibit strong increases from z = 0 to $\rm <z>$ = 1.15 and further to $\rm <z>$ = 2.2 and 4.8, particularly amongst galaxies above the SF galaxy main sequence (MS). The galaxies with highest SFRs at $\rm <z>$ = 2.2 and 4.8 have gas masses 100 times that of the Milky Way and gas mass fractions reaching 50 to 80\%, i.e. gas masses 1 - 4$\times$ their stellar masses. For the full sample of galaxies, we find a single, very simple SF law: $\rm SFR \propto M_{\rm ISM}^{0.9}$, i.e. a `linear' dependence on the ISM mass -- on and above the MS. Thus, the galaxies above the MS are converting their larger ISM masses into stars on a timescale similar to those on the MS. At z $> 1$, the entire population of star-forming galaxies has $\sim$5 - 10$\times$ shorter gas depletion times ($\sim0.2$ Gyr) than galaxies at low redshift. These {\bf shorter depletion times are due to a different, dominant mode of SF in the early universe} -- dynamically driven by compressive, high dispersion gas motions and/or galaxy interactions. The dispersive gas motions are a natural consequence of the extraordinarily high gas accretion rates which must occur to maintain the prodigious SF.

We present SIGAME (SImulator of GAlaxy Millimetre/submillimetre Emission), a new numerical code designed to simulate the 12CO rotational line emission spectrum of galaxies. Using sub-grid physics recipes to post-process the outputs of smoothed particle hydrodynamics (SPH) simulations, a molecular gas phase is condensed out of the hot and partly ionized SPH gas. The gas is subjected to far-UV radiation fields and cosmic ray ionization rates which are set to scale with the local star formation rate volume density. Level populations and radiative transport of the CO lines are solved with the 3-D radiative transfer code LIME. We have applied SIGAME to cosmological SPH simulations of three disc galaxies at z=2 with stellar masses in the range ~(0.5-2)x10^11 Msun and star formation rates ~40-140 Msun/yr. Global CO luminosities and line ratios are in agreement with observations of disc galaxies at z~2 up to and including J=3-2 but falling short of the few existing J=5-4 observations. The central 5 kpc regions of our galaxies have CO 3-2/1-0 and 7-6/1-0 brightness temperature ratios of ~0.55-0.65 and ~0.02-0.08, respectively, while further out in the disc the ratios drop to more quiescent values of ~0.5 and <0.01. Global CO-to-H2 conversion (alpha_CO) factors are ~=1.5 Msun*pc^2/(K km s/1), i.e. ~2-3 times below typically adopted values for disc galaxies, and alpha_CO increases with radius, in agreement with observations of nearby galaxies. Adopting a top-heavy Giant Molecular Cloud (GMC) mass spectrum does not significantly change the results. Steepening the GMC density profile leads to higher global line ratios for J_up>=3 and CO-to-H2 conversion factors [~=3.6 Msun*pc^2/(K km/s)].

The COSMOS-Legacy survey is a 4.6 Ms Chandra program that has imaged 2.2 deg$^2$ of the COSMOS field with an effective exposure of $\simeq$160 ks over the central 1.5 deg$^2$ and of $\simeq$80 ks in the remaining area. The survey is the combination of 56 new observations, obtained as an X-ray Visionary Project, with the previous C-COSMOS survey. We describe the reduction and analysis of the new observations and the properties of 2273 point sources detected above a spurious probability of 2$\times 10^{-5}$. We also present the updated properties of the C-COSMOS sources detected in the new data. The whole survey includes 4016 point sources (3814, 2920 and 2440 in the full, soft and hard band). The limiting depths are 2.2 $\times$ 10$^{-16}$, 1.5 $\times$ 10$^{-15}$ and 8.9$\times$ 10$^{-16}$ ${\rm erg~cm}^{-2}~{\rm s}^{-1}$ in the 0.5-2, 2-10 and 0.5-10 keV bands, respectively. The observed fraction of obscured AGN with column density $> 10^{22}$ cm$^{-2}$ from the hardness ratio (HR) is $\sim$50$^{+17}_{-16}$%. Given the large sample, we compute source number counts in the hard and soft bands, significantly reducing the uncertainties of 5-10%. For the first time, we compute number counts for obscured (HR$>$-0.2) and unobscured (HR$<$-0.2) sources and find significant differences between the two populations in the soft band. Due to the un-precedent large exposure, COSMOS-Legacy area is 3 times larger than surveys at similar depth and its depth is 3 times fainter than surveys covering similar area. The area-flux region occupied by COSMOS-Legacy is likely to remain unsurpassed for years to come.

We present S\'IGAME simulations of the [CII]157.7$\mu$ fine structure line emission from cosmological smoothed particle hydrodynamics (SPH) simulations of seven main sequence galaxies at z=2. Using sub-grid physics prescriptions the gas in our simulations is modeled as a multi-phased interstellar medium (ISM) comprised of molecular gas residing in giant molecular clouds, an atomic gas phase associated with photo-dissociation regions (PDRs) at the cloud surfaces, and a diffuse, ionized gas phase. Adopting logotropic cloud density profiles and accounting for heating by the local FUV radiation field and cosmic rays by scaling both with local star formation rate (SFR) volume density, we calculate the [CII] emission using a photon escape probability formalism. The [CII] emission peaks in the central $\lesssim$1 kpc of our galaxies as do the SFR radial profiles, with most [CII] ($\gtrsim$70%) originating in the molecular gas phase, whereas further out ($\gtrsim$2 kpc), the atomic/PDR gas dominates ($\gtrsim$90%) the [CII] emission, no longer tracing on-going star formation. Throughout, the ionized gas contribution is negligible ($\lesssim$3%). The [CII] luminosity vs. SFR ([CII]-SFR) relationship, integrated as well as spatially resolved (on scales of 1 kpc), delineated by our simulated galaxies is in good agreement with the corresponding relations observed locally and at high redshifts. In our simulations, the molecular gas dominates the [CII] budget at SFR$\gtrsim$20 M$_{\odot}$yr$^{-1}$ (SFR density $\gtrsim$0.5 M$_{\odot}$yr$^{-1}$kpc$^{-2}$), while atomic/PDR gas takes over at lower SFRs, suggesting a picture in which [CII] predominantly traces the molecular gas in high-density/pressure regions where star formation is on-going, and otherwise reveals the atomic/PDR gas phase.

Two millimeter observations of the MACS J1149.6+2223 cluster have detected a source that was consistent with the location of the lensed MACS1149-JD galaxy at z=9.6. A positive identification would have rendered this galaxy as the youngest dust forming galaxy in the universe. Follow up observation with the AzTEC 1.1 mm camera and the IRAM NOrthern Extended Millimeter Array (NOEMA) at 1.3 mm have not confirmed this association. In this paper we show that the NOEMA observations associate the 2 mm source with [PCB2012] 2882 ([PCB2012] 2882 is the NED-searchable name for this source.), source number 2882 in the Hubble Space Telescope ( HST) Cluster Lensing and Supernova (CLASH) catalog of MACS J1149.6+2223. This source, hereafter referred to as CLASH 2882, is a gravitationally lensed spiral galaxy at z=0.99. We combine the GISMO 2 mm and NOEMA 1.3 mm fluxes with other (rest frame) UV to far-IR observations to construct the full spectral energy distribution (SED) of this galaxy, and derive its star formation history, and stellar and interstellar dust content. The current star formation rate of the galaxy is 54/mu Msun yr-1, and its dust mass is about 5 10^7/mu Msun, where mu is the lensing magnification factor for this source, which has a mean value of 2.7. The inferred dust mass is higher than the maximum dust mass that can be produced by core collapse supernovae (CCSN) and evolved AGB stars. As with many other star forming galaxies, most of the dust mass in CLASH 2882 must have been accreted in the dense phases of the ISM.

We report a discovery of 6 massive galaxies with both extremely large Lya equivalent width and evolved stellar population at z ~ 3. These MAssive Extremely STrong Lya emitting Objects (MAESTLOs) have been discovered in our large-volume systematic survey for strong Lya emitters (LAEs) with twelve optical intermediate-band data taken with Subaru/Suprime-Cam in the COSMOS field. Based on the SED fitting analysis for these LAEs, it is found that these MAESTLOs have (1) large rest-frame equivalent width of EW_0(Lya) ~ 100--300 A, (2) M_star ~ 10^10.5--10^11.1 M_sun, and (3) relatively low specific star formation rates of SFR/M_star ~ 0.03--1 Gyr^-1. Three of the 6 MAESTLOs have extended Ly$\alpha$ emission with a radius of several kpc although they show very compact morphology in the HST/ACS images, which correspond to the rest-frame UV continuum. Since the MAESTLOs do not show any evidence for AGNs, the observed extended Lya emission is likely to be caused by star formation process including the superwind activity. We suggest that this new class of LAEs, MAESTLOs, provides a missing link from star-forming to passively evolving galaxies at the peak era of the cosmic star-formation history.

Context. A possible correlation between CO luminosity (L_CO ) and its line width (FWHM) has been suggested and denied in the literature. Such claims were often based on a small, or heterogeneous sample of galaxies, and thus inconclusive. Aims. We aim to prove or dis-prove the L_CO -FWHM correlation. Methods. We compile a large sample of submm galaxies at z>2 from the literature, and investigate the L_CO-FWHM relation. Results. After carefully evaluating the selection effects and uncertainties such as inclination and magnification via gravitational lensing, we show that there exist a weak but significant correlation between L_CO and FWHM. We also discuss a feasibility to measure the cosmological distance using the correlation.

We present the first spectroscopic confirmation of an ultra-massive galaxy at redshift z>3 using data from Keck-NIRSPEC, VLT-Xshooter, and GTC-Osiris. We detect strong [OIII] and Ly$\alpha$ emission, and weak [OII], CIV, and HeII, placing C1-23152 at a spectroscopic redshift of $z_{spec}$=3.351. The modeling of the emission-line corrected spectral energy distribution results in a best-fit stellar mass of $M_{*}=3.1^{+0.6}_{-0.7}\times10^{11} M_{\odot}$, a star-formation rate of <7 $M_{\odot} yr^{-1}$, and negligible dust extinction. The stars appear to have formed in a short intense burst ~300-500 Myr prior to the observation epoch, setting the formation redshift of this galaxy at z~4.1. From the analysis of the line ratios and widths, and the observed flux at 24$\mu$m, we confirm the presence of a luminous hidden active galactic nucleus (AGN), with bolometric luminosity of ~$10^{46}erg s^{-1}$. Potential contamination to the observed SED from the AGN continuum is constrained, placing a lower limit on the stellar mass of $2\times10^{11} M_{\odot}$. HST/WFC3 $H_{160}$ and ACS $I_{814}$ images are modeled, resulting in an effective radius of $r_{e}$~1 kpc in the $H_{160}$ band and a Sersic index n~4.4. This object may be a prototype of the progenitors of local most massive elliptical galaxies in the first 2 Gyr of cosmic history, having formed most of its stars at z>4 in a highly dissipative, intense, and short burst of star formation. C1-23152 is completing its transition to a post-starburst phase while hosting a powerful AGN, potentially responsible for the quenching of the star formation activity.

We study the physical properties of a sample of 6 SMGs in the COSMOS field, spectroscopically confirmed to lie at z>4. We use new GMRT 325 MHz and 3 GHz JVLA data to probe the rest-frame 1.4 GHz emission at z=4, and to estimate the sizes of the star-forming (SF) regions of these sources, resp. Combining our size estimates with those available in the literature for AzTEC1 and AzTEC3 we infer a median radio-emitting size for our z>4 SMGs of (0.63"+/-0.12")x(0.35"+/-0.05") or 4.1x2.3 kpc^2 (major times minor axis; assuming z=4.5) or lower if we take the two marginally resolved SMGs as unresolved. This is consistent with the sizes of SF regions in lower-redshift SMGs, and local normal galaxies, yet higher than the sizes of SF regions of local ULIRGs. Our SMG sample consists of a fair mix of compact and more clumpy systems with multiple, perhaps merging, components. With an average formation time of ~280 Myr, derived through modeling of the UV-IR SEDs, the studied SMGs are young systems. The average stellar mass, dust temperature, and IR luminosity we derive are M*~1.4x10^11 M_sun, T_dust~43 K, and L_IR~1.3x10^13L_sun, resp. The average L_IR is up to an order of magnitude higher than for SMGs at lower redshifts. Our SMGs follow the correlation between dust temperature and IR luminosity as derived for Herschel-selected 0.1<z<2 galaxies. We study the IR-radio correlation for our sources and find a deviation from that derived for z<3 ULIRGs (<q_IR>=1.95+/-0.26 for our sample, compared to q~2.6 for IR luminous galaxies at z<2). In summary, we find that the physical properties derived for our z>4 SMGs put them at the high end of the L_IR-T_dust distribution of SMGs, and that our SMGs form a morphologically heterogeneous sample. Thus, further in-depth analyses of large, statistical samples of high-redshift SMGs are needed to fully understand their role in galaxy formation and evolution.

The relationship between galaxy star formation rates (SFR) and stellar masses ($M_\ast$) is re-examined using a mass-selected sample of $\sim$62,000 star-forming galaxies at $z \le 1.3$ in the COSMOS 2-deg$^2$ field. Using new far-infrared photometry from $Herschel$-PACS and SPIRE and $Spitzer$-MIPS 24 $\mu$m, along with derived infrared luminosities from the NRK method based on galaxies' locations in the restframe color-color diagram $(NUV - r)$ vs. $(r - K)$, we are able to more accurately determine total SFRs for our complete sample. At all redshifts, the relationship between median $SFR$ and $M_\ast$ follows a power-law at low stellar masses, and flattens to nearly constant SFR at high stellar masses. We describe a new parameterization that provides the best fit to the main sequence and characterizes the low mass power-law slope, turnover mass, and overall scaling. The turnover in the main sequence occurs at a characteristic mass of about $M_{0} \sim 10^{10} M_{\odot}$ at all redshifts. The low mass power-law slope ranges from 0.9-1.3 and the overall scaling rises in SFR as a function of $(1+z)^{4.12 \pm 0.10}$. A broken power-law fit below and above the turnover mass gives relationships of $SFR \propto M_{*}^{0.88 \pm 0.06}$ below the turnover mass and $SFR \propto M_{*}^{0.27 \pm 0.04}$ above the turnover mass. Galaxies more massive than $M_\ast \gtrsim 10^{10}\ M_{\rm \odot}$ have on average, a much lower specific star formation rate (sSFR) than would be expected by simply extrapolating the traditional linear fit to the main sequence found for less massive galaxies.

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).

RXJ0848.6+4453 (Lynx W) at redshift 1.27 is part of the Lynx Supercluster of galaxies. Our analysis of stellar populations and star formation history in the cluster covers 24 members and is based on deep optical spectroscopy from Gemini North and imaging data from HST. Focusing on the 13 bulge-dominated galaxies for which we can determine central velocity dispersions, we find that these show a smaller evolution of sizes and velocity dispersions than reported for field galaxies and galaxies in poorer clusters. The galaxies in RXJ0848.6+4453 populate the Fundamental Plane similar to that found for lower redshift clusters with a zero point offset corresponding to an epoch of last star formation at z_form= 1.95+-0.2. The spectra of the galaxies in RXJ0848.6+4453 are dominated by young stellar populations at all galaxy masses and in many cases show emission indicating low level on-going star formation. The average age of the young stellar populations (estimated from H-zeta) is consistent with a major star formation episode 1-2 Gyr prior, which in turn agrees with z_form=1.95. Galaxies dominated by young stellar populations are distributed throughout the cluster. We speculate that low level star formation has not yet been fully quenched in the center of this cluster may be because the cluster is significantly poorer than other clusters previously studied at similar redshifts, which appear to have very little on-going star formation in their centers.

The role of major mergers in galaxy and black hole formation is not well constrained. To help address this, we develop an automated method to identify late-stage galaxy mergers before coalescence of the galactic cores. The resulting sample of mergers is distinct from those obtained using pair-finding and morphological indicators. Our method relies on median-filtering of high-resolution images in order to distinguish two concentrated galaxy nuclei at small separations. Using mock images, we derive statistical contamination and incompleteness corrections for the fraction of late-stage mergers. We apply our new method to a magnitude-limited (I < 23) sample of 44,164 galaxies from the COSMOS HST/ACS catalog. Using a mass-complete sample with $\log M_*/M_\odot > 10.6$ and $0.25 < z \leq 1.00$, we find ~5% of systems are late-stage mergers with separations between 2.2 and 8 kpc. Correcting for incompleteness and contamination, the fractional merger rate increases strongly with redshift as $(1+z)^{3.8\pm0.9}$, in agreement with earlier studies and with dark matter halo merger rates. Separating the sample into star-forming and quiescent galaxies shows that the merger rate for star-forming galaxies increases strongly redshift, $(1+z)^{4.5\pm1.3}$, while the merger rate for quiescent galaxies is consistent with no evolution, $(1+z)^{1.1\pm1.2}$. Limiting our sample to galaxies with spectroscopic redshifts from zCOSMOS, we find that the star formation rates and X-ray selected AGN activity in likely late-stage mergers are enhanced by factors of ~2 relative to a control sample. Combining our sample with more widely separated pairs, we find that $8\pm5\%$ of star formation and $20\pm8\%$ of AGN activity is triggered by close encounters (<143 kpc) or mergers, once more suggesting that major mergers are not the only channels for star formation and black hole growth. (abridged)

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.

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 improved constraints on the formation time scale of massive cluster galaxies based on rest-frame optical spectra of galaxies in a forming cluster located at z=2.16. The spectra are obtained with MOIRCS on the Subaru telescope with an integration time of ~7 hours. We achieve accurate redshift measurements by fitting SEDs using the spectra and broad-band photometry simultaneously, allowing us to identify probable cluster members. Clusters at low redshifts are dominated by quiescent galaxies, but we find that quiescent galaxies and star forming galaxies co-exist in this z=2 system. Interestingly, the quiescent galaxies form a weak red sequence in the process of forming. By stacking the spectra of star forming galaxies, we observe strong emission lines such as [OII] and [OIII] and we obtain a tentative hint of AGN activities in these galaxies. On the other hand, the stacked spectrum of the quiescent galaxies reveals a clear 4000A break with a possible CaII H+K absorption feature and strong emission lines such as [OII] are absent in the spectrum, confirming the quiescent nature of these galaxies. We then perform detailed spectral analyses of the stacked spectrum, which suggest that these massive quiescent galaxies formed at redshifts between 3 and 4 on a time scale of <~0.5 Gyr. This short formation time scale is not reproduced in recent numerical simulations. We discuss possible mechanisms for how these galaxies form 10^11 Msun stellar mass on a short time scale and become red and quiescent by z=2.

We estimate the galaxy stellar mass function and stellar mass density for star-forming and quiescent galaxies with 0.2<z<4. We construct a deep K<24 sample of 220000 galaxies selected using the UltraVISTA DR1 data release. Our analysis is based on precise 30-band photometric redshifts. By comparing these photometric redshifts with 10800 spectroscopic redshifts from the zCOSMOS bright and faint surveys, we find a precision of sigma(dz/(1+z))=0.008 at i<22.5 and sigma(dz/(1+zs))=0.03 at 1.5<z<4. We derive the stellar mass function and correct for the Eddington bias. We find a mass-dependent evolution of the global and star-forming populations. This mass-dependent evolution is a direct consequence of the star formation being quenched in galaxies more massive than M>10^10.7Msun. For the mass function of the quiescent galaxies, we do not find any significant evolution of the high-mass end at z<1; however we observe a clear flattening of the faint-end slope. From z~3 to z~1, the density of quiescent galaxies increases over the entire mass range. Their comoving stellar mass density increases by 1.6 dex between z~3 and z~1 and by less than 0.2dex at z<1. We infer the star formation history from the mass density evolution and we find an excellent agreement with instantaneous star formation rate measurements at z<1.5, while we find differences of 0.2dex at z>1.5 consistent with the expected uncertainties. We also develop a new method to infer the specific star formation rate from the mass function of star-forming galaxies. We find that the specific star formation rate of 10^10Msun galaxies increases continuously in the redshift range 1<z<4. Finally, we compare our results with a semi-analytical model and find that these models overestimate the density of low mass quiescent galaxies by an order of magnitude, while the density of low-mass star-forming galaxies is successfully reproduced.

We quantify the presence of Active Galactic nuclei (AGN) in a mass-complete (M_* >5e10 M_sun) sample of 123 star-forming and quiescent galaxies at 1.5 < z < 2.5, using X-ray data from the 4 Ms Chandra Deep Field-South (CDF-S) survey. 41+/-7% of the galaxies are detected directly in X-rays, 22+/-5% with rest-frame 0.5-8 keV luminosities consistent with hosting luminous AGN (L_0.5-8keV > 3e42 ergs/s). The latter fraction is similar for star-forming and quiescent galaxies, and does not depend on galaxy stellar mass, suggesting that perhaps luminous AGN are triggered by external effects such as mergers. We detect significant mean X-ray signals in stacked images for both the individually non-detected star-forming and quiescent galaxies, with spectra consistent with star formation only and/or a low luminosity AGN in both cases. Comparing star formation rates inferred from the 2-10 keV luminosities to those from rest-frame IR+UV emission, we find evidence for an X-ray excess indicative of low-luminosity AGN. Among the quiescent galaxies, the excess suggests that as many as 70-100% of these contain low- or high-luminosity AGN, while the corresponding fraction is lower among star-forming galaxies (43-65%). The ubiquitous presence of AGN in massive, quiescent z ~ 2 galaxies that we find provides observational support for the importance of AGN in impeding star formation during galaxy evolution.

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).

We present a high spatial-resolution HST/NICMOS imaging survey in the field of a known protocluster surrounding the powerful radio galaxy MRC1138-262 at z=2.16. Previously, we have shown that this field exhibits a substantial surface overdensity of red J-H galaxies. Here we focus on the stellar masses and galaxy effective radii in an effort to compare and contrast the properties of likely protocluster galaxies with their field counterparts and to look for correlations between galaxy structure and (projected) distance relative to the radio galaxy. We find a hint that quiescent, cluster galaxies are on average less dense than quiescent field galaxies of similar stellar mass and redshift. In fact, we find only two (of nine) quiescent protocluster galaxies are of simliar density to the majority of the massive, quiescent compact galaxies (SEEDs) found in several field surveys. Furthermore, there is some indication that the structural Sersic n parameter is higher (n ~ 3-4) on average for cluster galaxies compared to the field SEEDs (n ~ 1-2) This result may imply that the accelerated galaxy evolution expected (and observed) in overdense regions also extends to structural evolution presuming that massive galaxies began as dense (low n) SEEDs and have already evolved to be more in line with local galaxies of the same stellar mass.

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.

[Abridged]. Here, we report on the discovery of the galaxy counterpart of the z_abs=2.58 DLA on the line-of-sight to the z=3.07 quasar SDSS J091826.16+163609.0. The galaxy counterpart of the DLA is detected in the OIII 5007 and OII 3726,3729 emission lines redshifted into the NIR at an impact parameter of 16 kpc. Ly-alpha emission is not detected. The upper limit implies that Ly-alpha emission from this galaxy is suppressed by more than an order of magnitude. The DLA is amongst the most metal-rich DLAs studied so far at comparable redshifts. We find evidence for substantial depletion of refractory elements onto dust grains. Fitting the main metal line component of the DLA, which is located at z_abs=2.5832 and accounts for at least 85% of the total column density of low-ionisation species, we measure metal abundances from ZnII, SII, SiII, CrII, MnII, FeII and NiII of -0.12, -0.26, -0.46, -0.88, -0.92, -1.03 and -0.78, respectively. In addition, we detect absorption in the Lyman and Werner bands of hydrogen, which represents the first detection of H_2 molecules with X-shooter. The background quasar Q0918+1636 is amongst the reddest QSOs at redshifts 3.02<z<3.12 from the SDSS catalogue. Its UV to NIR spectrum is well fitted by a composite QSO spectrum reddened by SMC/LMC-like extinction curves at z_abs=2.58 with a significant amount of extinction given by A_V = 0.2 mag. This supports previous claims that there may be more metal-rich DLAs missing from current samples due to dust reddening of the background QSOs. The fact that there is evidence for dust both in the central emitting regions of the galaxy (as evidenced by the lack of Ly-alpha emission) and at an impact parameter of 16 kpc (as probed by the DLA) suggests that dust is widespread in this system.

We have initiated a survey using the newly commissioned X-shooter spectrograph to target candidate relatively metal-rich damped Lyman-alpha absorbers (DLAs). The spectral coverage of X-shooter allows us to search for not only Lyman-alpha emission, but also rest-frame optical emission lines. We have chosen DLAs where the strongest rest-frame optical lines ([OII], [OIII], Hbeta and Halpha) fall in the NIR atmospheric transmission bands. In this first paper resulting from the survey, we report on the discovery of the galaxy counterpart of the z_abs = 2.354 DLA towards the z=2.926 quasar Q2222$-0946. This DLA is amongst the most metal-rich z>2 DLAs studied so far at comparable redshifts and there is evidence for substantial depletion of refractory elements onto dust grains. We measure metallicities from ZnII, SiII, NiII, MnII and FeII of -0.46+/-0.07, -0.51+/-0.06, -0.85+/-0.06, -1.23+/-0.06, and -0.99+/-0.06, respectively. The galaxy is detected in the Lyman-alpha, [OIII] lambda4959,5007 Halpha emission lines at an impact parameter of about 0.8 arcsec (6 kpc at z_abs = 2.354). We infer a star-formation rate of 10 M_sun yr^-1, which is a lower limit due to the possibility of slit-loss. Compared to the recently determined Halpha luminosity function for z=2.2 galaxies the DLA-galaxy counterpart has a luminosity of L~0.1L^*_Halpha. The emission-line ratios are 4.0 (Lyalpha/Halpha) and 1.2 ([OIII]/Halpha). The Lyalpha line shows clear evidence for resonant scattering effects, namely an asymmetric, redshifted (relative to the systemic redshift) component and a much weaker blueshifted component. The fact that the blueshifted component is relatively weak indicates the presence of a galactic wind. The properties of the galaxy counterpart of this DLA is consistent with the prediction that metal-rich DLAs are associated with the most luminous of the DLA-galaxy counterparts.

The presence of massive, compact, quiescent galaxies at z>2 presents a major challenge for theoretical models of galaxy formation and evolution. Using one of the deepest large public near-IR surveys to date, we investigate in detail the correlations between star formation and galaxy structural parameters (size, stellar mass, and surface density) from z=2 to the present. At all redshifts, massive quiescent galaxies (i.e. those with little or no star formation) occupy the extreme high end of the surface density distribution and follow a tight mass-size correlation, while star-forming galaxies show a broad range of both densities and sizes. Conversely, galaxies with the highest surface densities comprise a nearly-homogeneous population with little or no ongoing star formation, while less dense galaxies exhibit high star-formation rates and varying levels of dust obscuration. Both the sizes and surface densities of quiescent galaxies evolve strongly from z=2-0; we parameterize this evolution for both populations with simple power law functions and present best-fit parameters for comparison to future theoretical models. Higher-mass quiescent galaxies undergo faster structural evolution, consistent with previous results. Interestingly, star-forming galaxies' sizes and densities evolve at rates similar to those of quiescent galaxies. It is therefore possible that the same physical processes drive the structural evolution of both populations, suggesting that "dry mergers" may not be the sole culprit in this size evolution.

One of the main science goal of the future European Extremely Large Telescope will be to understand the mass assembly process in galaxies as a function of cosmic time. To this aim, a multi-object, AO-assisted integral field spectrograph will be required to map the physical and chemical properties of very distant galaxies. In this paper, we examine the ability of such an instrument to obtain spatially resolved spectroscopy of a large sample of massive (0.1<Mstellar<5e11Mo) galaxies at 2<z<6, selected from future large area optical-near IR surveys. We produced a set of about one thousand numerical simulations of 3D observations using reasonable assumptions about the site, telescope, and instrument, and about the physics of distant galaxies. These data-cubes were analysed as real data to produce realistic kinematic measurements of very distant galaxies. We then studied how sensible the scientific goals are to the observational (i.e., site-, telescope-, and instrument-related) and physical (i.e., galaxy-related) parameters. We specifically investigated the impact of AO performance on the science goal. We did not identify any breaking points with respect to the parameters (e.g., the telescope diameter), with the exception of the telescope thermal background, which strongly limits the performance in the highest (z>5) redshift bin. We find that a survey of Ngal galaxies that fulfil the range of science goals can be achieved with a ~90 nights program on the E-ELT, provided a multiplex capability M Ngal/8.

In order to investigate the structure and dynamics of the recently discovered massive (M_* > 10^11 M_sun) compact z~2 galaxies, cosmological hydrodynamical/N-body simulations of a proto-cluster region have been undertaken. At z=2, the highest resolution simulation contains ~5800 resolved galaxies, of which 509, 27 and 5 have M_* > 10^10 M_sun, > 10^11 M_sun and > 4x10^11 M_sun, respectively. Effective radii and characteristic stellar densities have been determined for all galaxies. At z=2, for the definitely well resolved mass range of M_* > 10^11 Msun, the mass-size relation is consistent with observational findings for the most compact z~2 galaxies. The very high velocity dispersion recently measured for a compact z~2 galaxy (~510 km/s; van Dokkum et al 2009) can be matched at about the 1-sigma level, although a somewhat larger mass than the estimated M_* ~ 2 x 10^11 M_sun is indicated. For the above mass range, the galaxies have an average axial ratio <b/a> = 0.64 +/- 0.02 with a dispersion of 0.1, an average rotation to 1D velocity dispersion ratio <v/sigma> = 0.46 +/- 0.06 with a dispersion of 0.3, and a maximum value of v/sigma ~ 1.1. Rotation and velocity anisotropy both contribute in flattening the compact galaxies. Some of the observed compact galaxies appear flatter than any of the simulated galaxies. Finally, it is found that the massive compact galaxies are strongly baryon dominated in their inner parts, with typical dark matter mass fractions of order only 20% inside of r=2R_eff.

We study the relation between size and star formation activity in a complete sample of 225 massive (M > 5 x 10^10 Msun) galaxies at 1.5<z<2.5, selected from the FIREWORKS UV-IR catalog of the CDFS. Based on stellar population synthesis model fits to the observed restframe UV-NIR SEDs, and independent MIPS 24 micron observations, 65% of galaxies are actively forming stars, while 35% are quiescent. Using sizes derived from 2D surface brightness profile fits to high resolution (FWHM_{PSF}~0.45 arcsec) groundbased ISAAC data, we confirm and improve the significance of the relation between star formation activity and compactness found in previous studies, using a large, complete mass-limited sample. At z~2, massive quiescent galaxies are significantly smaller than massive star forming galaxies, and a median factor of 0.34+/-0.02 smaller than galaxies of similar mass in the local universe. 13% of the quiescent galaxies are unresolved in the ISAAC data, corresponding to sizes <1 kpc, more than 5 times smaller than galaxies of similar mass locally. The quiescent galaxies span a Kormendy relation which, compared to the relation for local early types, is shifted to smaller sizes and brighter surface brightnesses and is incompatible with passive evolution. The progenitors of the quiescent galaxies, were likely dominated by highly concentrated, intense nuclear star bursts at z~3-4, in contrast to star forming galaxies at z~2 which are extended and dominated by distributed star formation.

We investigate the star formation history of the universe using FIREWORKS, a multiwavelength survey of the CDFS. We study the evolution of the specific star formation rate (sSFR) with redshift in different mass bins from z = 0 to z ~ 3. We find that the sSFR increases with redshift for all masses. The logarithmic increase of the sSFR with redshift is nearly independent of mass, but this cannot yet be verified at the lowest-mass bins at z > 0.8, due to incompleteness. We convert the sSFRs to a dimensionless growth rate to facilitate a comparison with a semi-analytic galaxy formation model that was implemented on the Millennium Simulation. The model predicts that the growth rates and sSFRs increase similarly with redshift for all masses, consistent with the observations. However, we find that for all masses, the inferred observed growth rates increase more rapidly with redshift than the model predictions. We discuss several possible causes for this discrepancy, ranging from field-to-field variance, conversions to SFR, and shape of the IMF. We find that none of these can solve the discrepancy completely. We conclude that the models need to be adapted to produce the steep increase in growth rate between redshift z=0 and z=1.

We present a comparison between the observed color distribution, number and mass density of massive galaxies at 1.5 < z < 3 and a model by Hopkins et al. that relates the quasar and galaxy population on the basis of gas-rich mergers. In order to test the hypothesis that quiescent red galaxies are formed after a gas-rich merger involving quasar activity, we confront photometry of massive (M > 4x10^10 Msun) galaxies extracted from the FIRES, GOODS-South, and MUSYC surveys, together spanning an area of 496 arcmin^2, with synthetic photometry from hydrodynamical merger simulations. As in the Hopkins et al. (2006b) model, we use the observed quasar luminosity function to estimate the merger rate. We find that the synthetic U-V and V-J colors of galaxies that had a quasar phase in their past match the colors of observed galaxies that are best characterized by a quiescent stellar population. At z ~ 2.6, the observed number and mass density of quiescent red galaxies with M > 4x10^10 Msun is consistent with the model in which every quiescent massive galaxy underwent a quasar phase in the past. At z ~ 1.9, 2.8 times less quiescent galaxies are observed than predicted by the model as descendants of higher redshift quasars. The merger model also predicts a large number of galaxies undergoing merger-driven star formation. We find that the predicted number and mass density accounts for 30-50% of the observed massive star-forming galaxies. However, their colors do not match those of observed star-forming galaxies. In particular, the colors of dusty red galaxies are not reproduced by the simulations. Several possible origins of this discrepancy are discussed. The observational constraints on the validity of the model are currently limited by cosmic variance and uncertainties in stellar population synthesis and radiative transfer.

Models of disk galaxy formation commonly predict the existence of an extended reservoir of accreted hot gas surrounding massive spirals at low redshift. As a test of these models, we use X-ray and H-alpha data of the two massive, quiescent edge-on spirals NGC 5746 and NGC 5170 to investigate the amount and origin of any hot gas in their halos. Contrary to our earlier claim, the Chandra analysis of NGC 5746, employing more recent calibration data, does not reveal any significant evidence for diffuse X-ray emission outside the optical disk, with a 3-sigma upper limit to the halo X-ray luminosity of 4e39 erg/s. An identical study of the less massive NGC 5170 also fails to detect any extraplanar X-ray emission. By extracting hot halo properties of disk galaxies formed in cosmological hydrodynamical simulations, we compare these results to expectations for cosmological accretion of hot gas by spirals. For Milky Way-sized galaxies, these high-resolution simulations predict hot halo X-ray luminosities which are lower by a factor of ~2 compared to our earlier results reported by Toft et al. (2002). We find the new simulation predictions to be consistent with our observational constraints for both NGC 5746 and NGC 5170, while also confirming that the hot gas detected so far around more actively star-forming spirals is in general probably associated with stellar activity in the disk. Observational results on quiescent disk galaxies at the high-mass end are nevertheless providing powerful constraints on theoretical predictions, and hence on the assumed input physics in numerical studies of disk galaxy formation and evolution.

We present an analysis of galaxies in the CDF-South. We find a tight relation to z=3 between color and size at a given mass, with red galaxies being small, and blue galaxies being large. We show that the relation is driven by stellar surface density or inferred velocity dispersion: galaxies with high surface density are red and have low specific star formation rates, and galaxies with low surface density are blue and have high specific star formation rates. Surface density and inferred velocity dispersion are better correlated with specific star formation rate and color than stellar mass. Hence stellar mass by itself is not a good predictor of the star formation history of galaxies. In general, galaxies at a given surface density have higher specific star formation rates at higher redshift. Specifically, galaxies with a surface density of 1-3 10^9 Msun/kpc^2 are "red and dead" at low redshift, approximately 50% are forming stars at z=1, and almost all are forming stars by z=2. This provides direct additional evidence for the late evolution of galaxies onto the red sequence. The sizes of galaxies at a given mass evolve like 1/(1+z)^(0.59 +- 0.10). Hence galaxies undergo significant upsizing in their history. The size evolution is fastest for the highest mass galaxies, and quiescent galaxies. The persistence of the structural relations from z=0 to z=2.5, and the upsizing of galaxies imply that a relation analogous to the Hubble sequence exists out to z=2.5, and possibly beyond. The star forming galaxies at z >= 1.5 are quite different from star forming galaxies at z=0, as they have likely very high gas fractions, and star formation time scales comparable to the orbital time.

Using deep near-infrared spectroscopy Kriek et al. (2006) found that ~45% of massive galaxies at z~2.3 have evolved stellar populations and little or no ongoing star formation. Here we determine the sizes of these quiescent galaxies using deep, high-resolution images obtained with HST/NIC2 and laser guide star-assisted Keck/AO. Considering that their median stellar mass is 1.7x10^11 Solar masses the galaxies are remarkably small, with a median effective radius of 0.9 kpc. Galaxies of similar mass in the nearby Universe have sizes of ~5 kpc and average stellar densities which are two orders of magnitude lower than the z~2.3 galaxies. These results extend earlier work at z~1.5 and confirm previous studies at z>2 which lacked spectroscopic redshifts and imaging of sufficient resolution to resolve the galaxies. Our findings demonstrate that fully assembled early-type galaxies make up at most ~10% of the population of K-selected quiescent galaxies at z~2.3, effectively ruling out simple monolithic models for their formation. The galaxies must evolve significantly after z~2.3, through dry mergers or other processes, consistent with predictions from hierarchical models.

We present HST NICMOS+ACS and Spitzer IRAC+MIPS observations of 41 galaxies at 2<z<3.5 in the FIRES MS1054 field with red and blue rest-frame optical colors. About half of the galaxies are very compact (effective radii r_e < 1 kpc) at rest-frame optical wavelengths, the others are extended (1< r_e < 10 kpc). For reference, 1 kpc corresponds to 0.12 arcsec at z=2.5 in the adopted cosmology. We separate actively star forming galaxies from quiescent galaxies by modeling their rest-frame UV-NIR SEDs. The star forming galaxies span the full range of sizes, while the quiescent galaxies all have r_e<2kpc. In the redshift range where MIPS 24 micron imaging is a sensitive probe of re-radiated dust emission (z<2.5), the 24 micron fluxes confirm that the light of the small quiescent galaxies is dominated by old stars, rather than dust-enshrouded star formation or AGN activity. The inferred surface mass densities and velocity dispersions for the quiescent galaxies are very high compared to those in local galaxies. The galaxies follow a Kormendy relation (between surface brightness and size) with approximately the same slope as locally, but shifted to brighter surface brightnesses, consistent with a mean stellar formation redshift of z_f~5. This paper demonstrates a direct relation between star formation activity and size at z~2.5, and the existence of a significant population of massive, extremely dense, old stellar systems without readily identifiable counterparts in the local universe.

We present the rest-frame optical (B, V, and R band) luminosity functions (LFs) of galaxies at 2<z<3.5, measured from a K-selected sample constructed from the deep NIR MUSYC, the ultradeep FIRES, and the GOODS-CDFS. This sample is unique for its combination of area and range of luminosities. The faint-end slopes of the LFs at z>2 are consistent with those at z~0. The characteristic magnitudes are significantly brighter than the local values (e.g., ~1.2 mag in the R band), while the measured values for Phi_star are typically ~5 times smaller. The B-band luminosity density at z~2.3 is similar to the local value, and in the R band it is ~2 times smaller than the local value. We present the LF of Distant Red Galaxies (DRGs), which we compare to that of non-DRGs. While DRGs and non-DRGs are characterized by similar LFs at the bright end, the faint-end slope of the non-DRG LF is much steeper than that of DRGs. The contribution of DRGs to the global densities down to the faintest probed luminosities is 14%-25% in number and 22%-33% in luminosity. From the derived rest-frame U-V colors and stellar population synthesis models, we estimate the mass-to-light ratios (M/L) of the different subsamples. The M/L ratios of DRGs are ~5 times higher (in the R and V bands) than those of non-DRGs. The global stellar mass density at 2<z<3.5 appears to be dominated by DRGs, whose contribution is of order ~60%-80% of the global value. Qualitatively similar results are obtained when the population is split by rest-frame U-V color instead of observed J-K color.

We present deep, high angular-resolution HST/NICMOS imaging in the Hubble Deep Field South (HDF-S), focusing on a subset of 14 Distant Red Galaxies (DRGs) at z ~ 2.5 galaxies that have been pre-selected to have J-K > 2.3. We find a clear trend between the rest-frame optical sizes of these sources and their luminosity-weighted stellar ages as inferred from their broad-band spectral energy distributions (SEDs). Galaxies whose SEDs are consistent with being dusty and actively star forming generally show extended morphologies in the NICMOS images (r_e >~ 2 kpc), while the 5 sources which are not vigorously forming stars are extremely compact (r_e <~ 1 kpc). This trend suggests a direct link between the mean ages of the stars and the size and density of the galaxies and supports the conjecture that early events quench star-formation and leave compact remnants. Furthermore, the compact galaxies have stellar surface mass densities which exceed those of local galaxies by more than an order of magnitude. The existence of such massive dense galaxies presents a problem for models of early-type galaxy formation and evolution. Larger samples of DRGs and higher spatial resolution imaging will allow us to determine the universality of the results presented here for a small sample.

Models of disk galaxy formation commonly predict the existence of an extended reservoir of hot gas surrounding massive spirals at low redshift. As a test of these models, we have obtained X-ray and optical data of the two massive edge-on spirals NGC 5746 and NGC 5170, in order to investigate the amount and origin of hot gas in their disks and halos. Chandra observations of NGC 5746 reveal evidence for diffuse X-ray emission with a total luminosity of ~7 x 10^39 erg/s surrounding this galaxy out to at least ~20 kpc from the disk, whereas an identical study of the less massive NGC 5170 fails to detect any extraplanar X-ray emission. Unlike the case for other disk galaxies with detected X-ray halos, the halo emission around NGC 5746 is not accompanied by extraplanar H-alpha or radio emission, and there is no evidence for significant nuclear or starburst activity in the disk. In contrast to these other cases, the emission around NGC 5746 therefore appears to arise from the cooling of externally accreted material rather than from disk outflows. To verify this idea, we present results of cosmological simulations of galaxy formation and evolution, showing our observations to be in good agreement with expectations for cosmological accretion, while also confirming that the X-ray halos of other spirals do not fit well into an accretion scenario. We find that the estimated cooling rate of hot halo gas around NGC 5746 would provide sufficient material for star formation in the disk to proceed at its present rate. This lends support to the idea that a supply of hot ambient gas is potentially available as fuel for star formation in massive, nearby spirals, and suggests that accretion of hot gas could be important for maintaining the stellar disks of such galaxies.

We present first results of a spectroscopic survey targeting K-selected galaxies at z=2.0-2.7 using the GNIRS instrument on Gemini-South. We obtained near-infrared spectra with a wavelength coverage of 1.0-2.5 micron for 26 K-bright galaxies (K<19.7) selected from the MUSYC survey using photometric redshifts. We successfully derived spectroscopic redshifts for all 26 galaxies using rest-frame optical emission lines or the redshifted Balmer/4000 Angstrom break. Twenty galaxies have spectroscopic redshifts in the range 2.0<z<2.7, for which bright emission lines like Halpha and [OIII] fall in atmospheric windows. Surprisingly, we detected no emission lines for nine of these 20 galaxies. The median 2 sigma upper limit on the rest-frame equivalent width of Halpha for these nine galaxies is ~10 Angstrom. The stellar continuum emission of these same nine galaxies is best fitted by evolved stellar population models. The best-fit star formation rate (SFR) is zero for five out of nine galaxies, and consistent with zero within 1 sigma for the remaining four. Thus, both the Halpha measurements and the independent stellar continuum modeling imply that 45% of our K-selected galaxies are not forming stars intensely. This high fraction of galaxies without detected line emission and low SFRs may imply that the suppression of star formation in massive galaxies occurs at higher redshift than is predicted by current CDM galaxy formation models. However, obscured star formation may have been missed, and deep mid-infrared imaging is needed to clarify this situation.

We present a survey of the extinction properties of ten lensing galaxies, in the redshift range z = 0.04 - 1.01, using multiply lensed quasars imaged with the ESO VLT in the optical and near infrared. The multiple images act as 'standard light sources' shining through different parts of the lensing galaxy, allowing for extinction studies by comparison of pairs of images. We explore the effects of systematics in the extinction curve analysis, including extinction along both lines of sight and microlensing, using theoretical analysis and simulations. In the sample, we see variation in both the amount and type of extinction. Of the ten systems, seven are consistent with extinction along at least one line of sight. The mean differential extinction for the most extinguished image pair for each lens is A(V) = 0.56 +- 0.04, using Galactic extinction law parametrization. The corresponding mean R_V = 2.8 +- 0.4 is consistent with that of the Milky Way at R_V = 3.1, where R_V = A(V)/E(B-V). We do not see any strong evidence for evolution of extinction properties with redshift. Of the ten systems, B1152+199 shows the strongest extinction signal of A(V) = 2.43 +- 0.09 and is consistent with a Galactic extinction law with R_V = 2.1 +- 0.1. Given the similar redshift distribution of SN Ia hosts and lensing galaxies, a large space based study of multiply imaged quasars would be a useful complement to future dark energy SN Ia surveys, providing independent constraints on the statistical extinction properties of galaxies up to z~1.

We present the evolution of the luminosity-size and stellar mass-size relations of luminous (L_V>3.4x10^10h_70^-2L_sun) and of massive (M_*>3x10^10h_70^-2M_sun) galaxies in the last ~11 Gyr. We use very deep near-infrared images of the Hubble Deep Field-South and the MS1054-03 field in the J_s, H and K_s bands from FIRES to retrieve the sizes in the optical rest-frame for galaxies with z>1. We combine our results with those from GEMS at 0.2<z<1 and SDSS at z~0.1 to achieve a comprehensive picture of the optical rest-frame size evolution from z=0 to z=3. Galaxies are differentiated according to their light concentration using the Sersic index n. For less concentrated objects, the galaxies at a given luminosity were typically ~3+-0.5 (+-2 sigma) times smaller at z~2.5 than those we see today. The stellar mass-size relation has evolved less: the mean size at a given stellar mass was \~2+-0.5 times smaller at z~2.5, evolving proportional to (1+z)^{-0.40+-0.06}. Simple scaling relations between dark matter halos and baryons in a hierarchical cosmogony predict a stronger (although consistent within the error bars) than observed evolution of the stellar mass-size relation. The observed luminosity-size evolution out to z~2.5 matches well recent infall model predictions for Milky-Way type objects. For low-n galaxies, the evolution of the stellar mass-size relation would follow naturally if the individual galaxies grow inside-out. For highly concentrated objects, the situation is as follows: at a given luminosity, these galaxies were ~2.7+-1.1 times smaller at z~2.5 (or put differently, were typically ~2.2+-0.7 mag brighter at a given size than they are today), and at a given stellar mass the size has evolved proportional to (1+z)^{-0.45+-0.10}.

We use near-infrared (NIR) spectroscopy obtained with GNIRS on Gemini, NIRSPEC on KECK, and ISAAC on the VLT to study the rest-frame optical continua of three `Distant Red Galaxies' (having Js - Ks > 2.3) at z>2. All three galaxy spectra show the Balmer/4000 Angstrom break in the rest-frame optical. The spectra allow us to determine spectroscopic redshifts from the continuum with an estimated accuracy dz/(1+z) ~ 0.001-0.04. These redshifts agree well with the emission line redshifts for the 2 galaxies with Halpha emission. This technique is particularly important for galaxies that are faint in the rest-frame UV, as they are underrepresented in high redshift samples selected in optical surveys and are too faint for optical spectroscopy. Furthermore, we use the break, continuum shape, and equivalent width of Halpha together with evolutionary synthesis models to constrain the age, star formation timescale, dust content, stellar mass and star formation rate of the galaxies. Inclusion of the NIR spectra in the stellar population fits greatly reduces the range of possible solutions for stellar population properties. We find that the stellar populations differ greatly among the three galaxies, ranging from a young dusty starburst with a small break and strong emission lines to an evolved galaxy with a strong break and no detected line emission. The dusty starburst galaxy has an age of 0.3 Gyr and a stellar mass of 1*10^11 Msun. The spectra of the two most evolved galaxies imply ages of 1.3-1.4 Gyr and stellar masses of 4*10^11 Msun. The large range of properties seen in these galaxies strengthens our previous much more uncertain results from broadband photometry. Larger samples are required to determine the relative frequency of dusty starbursts and (nearly) passively evolving galaxies at z~2.5.

Hot gaseous haloes surrounding galaxies and extending well beyond the distribution of stars are a ubiquitous prediction of galaxy formation scenarios. The haloes are believed to consist of gravitationally trapped gas with a temperature of millions of Kelvin. The existence of such hot haloes around massive elliptical galaxies has been established through their X-ray emission. While gas out-flowing from starburst spiral galaxies has been detected, searches for hot haloes around normal, quiescent spiral galaxies have so far failed, casting doubts on the fundamental physics in galaxy formation models. Here we present the first detection of a hot, large-scale gaseous halo surrounding a normal, quiescent spiral galaxy, NGC 5746, alleviating a long-standing problem for galaxy formation models. In contrast to starburst galaxies, where the X-ray halo can be powered by the supernova energy, there is no such power source in NGC 5746. The only compelling explanation is that we are here witnessing a galaxy forming from gradually in-flowing hot and dilute halo gas.

We determined the K band luminosity function (LF) of cluster galaxies at redshift z~1.2, using near-infrared images of three X-ray luminous clusters at z=1.11,1.24,1.27. The composite LF was derived down to M*+4, by means of statistical background subtraction, and is well described by a Schechter function with K*=20.5 AB mag and alpha=-1. From the K band composite LF we derived the stellar mass function of cluster galaxies. Using available X-ray mass profiles we determined the M/L ratios of these three clusters, which tend to be lower than those measured in the local universe. With these data, no significant difference can be seen between the shapes of the cluster galaxies LF and the LF of field galaxies at similar redshift. We also found no significant evolution out to z ~1.2 in the bright (<M*+4) part of the cluster galaxies LF probed in this study, apart from a brightening of ~1.3 mag of the characteristic magnitude of the high redshift LF. We confirm, and extend to higher redshift, the result from previous work that the redshift evolution of the characteristic magnitude M* is consistent with passive evolution of a stellar population formed at z>2. The results obtained in this work support and extend previous findings that most of the stars in bright galaxies were formed at high redshift, and that K-bright (M>10^11 Msun) galaxies were already in place at z ~ 1.2, at least in the central regions of X-ray luminous clusters. Together with recent results on the field galaxies stellar mass function, this implies that most of the stellar mass is already assembled in massive galaxies by z ~ 1, both in low and high density environments.

Using the deep multi-wavelength MUSYC, GOODS, and FIRES surveys we construct a stellar mass-limited sample of galaxies at 2<z<3. The sample comprises 294 galaxies with M>10^11 Solar masses distributed over four independent fields with a total area of almost 400 sq arcmin. The mean number density of massive galaxies in this redshift range is (2.2+-0.6) x 10^-4 Mpc^-3. We present median values and 25th and 75th percentiles for the distributions of observed R mags, observed J-K colors, and rest-frame UV continuum slopes, M/L(V) ratios, and U-V colors. The galaxies show a large range in all these properties. The ``median galaxy'' is faint in the observer's optical (R=25.9), red in the observed near-IR (J-K=2.48), has a rest-frame UV spectrum which is relatively flat (beta=-0.4), and rest-frame optical colors resembling those of nearby spiral galaxies (U-V=0.62). We determine which galaxies would be selected as Lyman break galaxies (LBGs) or Distant Red Galaxies (DRGs, having J-K>2.3) in this mass-limited sample. By number DRGs make up 69% of the sample and LBGs 20%, with a small amount of overlap. By mass DRGs make up 77% and LBGs 17%. Neither technique provides a representative sample of massive galaxies at 2<z<3 as they only sample the extremes of the population. As we show here, multi-wavelength surveys with high quality photometry are essential for an unbiased census of massive galaxies in the early Universe. The main uncertainty in this analysis is our reliance on photometric redshifts; confirmation of the results presented here requires extensive near-infrared spectroscopy of optically-faint samples.

We present UBVRIz' optical images of the 0.32 square degree Extended Hubble Deep Field South. This is one of four fields comprising the MUSYC survey, which is optimized for the study of galaxies at z=3, AGN demographics, and Galactic structure. We calculate corrected aperture photometry and its uncertainties and find through tests that these provide a significant improvement upon standard techniques. Our photometric catalog of 62968 objects is complete to a total magnitude of R_AB=25. We select z=3 Lyman break galaxy (LBG) candidates from their UVR colors and find a sky surface density of 1.4/sq. arcmin and an angular correlation function w(theta) = 2.3+-1.0 theta^{-0.8}. (Abridged)

This paper describes the observations and data reduction techniques for the version 2.0 images and catalog of the NICMOS Ultra Deep Field Treasury program. All sources discussed in this paper are based on detections in the combined NICMOS F110W and F160W bands only. The NICMOS images are drizzled to 0.09 arc second pixels and aligned to the ACS UDF F850LP image which was rebinned to the same pixel scale. These form the NICMOS version 2.0 UDF images. The catalog sources are chosen with a conservative detection limit to avoid the inclusion of numerous spurious sources. The catalog contains 1293 objects in the 144 x 144 arc sececonds NICMOS subfield of the UDF. The 5 sigma signal to noise level is an average 0.6 arc second diameter aperture AB magnitude of ~27.7 at 1.1 and 1.6 microns. The catalog sources, listed in order of right ascension, satisfy a minimum signal to noise criterion of 1.4 sigma in at least 7 contiguous pixels of the combined F110W and F160W image

We take advantage of the Hubble Ultra Deep Field (UDF) data to study the restframe optical and ultra violet (UV) morphologies of the novel population of Distant Red Galaxies (DRGs). Six galaxies with J-Ks > 2.3 are found to Ks=21.5, five of which have photometric redshifts z_phot > 2, corresponding to a surface density of 0.9/arcmin^2. The surface brightness distributions of the z_phot > 2 galaxies are better represented by exponential disks than R^{1/4}-laws. Two of the z_phot > 2 galaxies are extended, while three have compact morphologies. The restframe optical morphology of the z_phot > 2 galaxies is quite different from the restframe UV morphology: all the galaxies have red central components which dominate in the NICMOS H_{160}-band images, and distinct off-center blue features which show up in (and often dominate) the ACS images. The mean measured effective radius of the z_phot > 2 galaxies is <r_e> =1.9+/-1.4 kpc, similar (within the errors) to the mean size of LBGs at similar redshifts. All the DRGs are resolved in the ACS images, while four are resolved in the NICMOS images. Two of the z_phot > 2 galaxies are bright X-ray sources and hence host AGN. The diverse restframe optical and UV morphological properties of DRGs derived here suggest that they have complex stellar populations, consisting of both evolved populations that dominate the mass and the restframe optical light, and younger populations, which show up as patches of star formation in the restframe UV light; in many ways resembling the properties of normal local galaxies. This interpretation is supported by fits to the broadband SEDs, which for all five z_phot > 2 are best represented by models with extended star formation histories and substantial amounts of dust.

We derive the Ks-band luminosity function (LF) of the z=1.237 massive X-ray luminous cluster of galaxies RDCS J1252.9-2927. Photometric redshifts, derived from deep multi-wavelength BVRIzJKs data, and calibrated using a large subset of galaxies with spectroscopic redshifts, are used to separate the cluster galaxy population from the foreground and background field galaxy population. This allows for a simultaneous determination of the characteristic magnitude Ks* and faint end slope alpha of the LF without having to make an uncertain statistical background subtraction. The shape of the bright end of the derived LF is similar to that measured at similar restframe wavelengths (in the z-band) in local clusters, but the characteristic magnitude is brighter by Delta Mz*=-1.40(+0.49/-0.58) magnitudes, suggesting that the massive cluster ellipticals that dominate the bright end of the LF were already in place at z=1.237. The derived alpha is shallower (at the 2.2 sigma level) than the value measured in the z-band in clusters in the local universe, suggesting that RDCS J1252.9-2927 contains a relatively smaller fraction of low mass galaxies than clusters in the local universe. The results presented in this paper are a challenge for semi analytical hierarchical models of galaxy formation which predict the characteristic magnitude to grow fainter and the faint end slope to steepen with redshift as the massive galaxies break up into their progenitors. The observed evolution is consistent with a scenario in which clusters are composed of a population of massive galaxies which formed at high redshift (z>>1) and subsequently evolved passively, and a population of lower mass galaxies which are gradually accreted from the field, primarily at lower redshift (z<1). (Abridged)

X-ray luminosities and surface brightness profiles of the hot gas haloes of simulated disc galaxies at redshifts z=0-2 are presented. The galaxies are extracted from fully cosmological simulations and correspond in mass to the Milky Way. We find that the bolometric X-ray luminosities of the haloes decrease by a factor 4-10 from z~1 to z~0, reflecting the decrease in the rate at which hot halo gas cools out on to the disc. At all redshifts, most of the emission is found to originate within 10-15 kpc of the disc. When combined with models in which the evolution of disc X-ray luminosity is dominated by X-ray binaries, the predicted halo luminosities at z~1 show good agreement with constraints from spiral galaxies in Chandra Deep Field data. There is an indication that haloes with a metal abundance of 0.3 solar overpredict observed X-ray luminosities at z~1, suggesting that halo metallicities are lower than this value. Prospects for direct detection of the haloes of Milky Way-sized galaxies with current and future X-ray instrumentation are discussed. It is found that XEUS should be able to single out the halo emission of highly inclined Milky Way-sized disc galaxies out to approximately z=0.3. For such galaxies in this redshift interval, we estimate a lower limit to the surface density of detectable haloes on the sky of ~10 deg^(-2). More generally, owing to their luminosity evolution, the optimum redshifts at which to observe such haloes could be 0.5<z<1, depending on their assembly history.

We present deep SofI and ISAAC near-infrared imaging data of the X-ray luminous galaxy cluster RDCS J1252.9-2927. The ISAAC data were taken at the ESO Very Large Telescope under very good seeing conditions and reach limiting Vega magnitudes of 25.6 and 24.1 in the J- and Ks-bands respectively. The image quality is 0.45 arc seconds in both passbands. We use these data to construct a colour-magnitude (C-M) diagram of galaxies that are within 20 arc seconds of the cluster center and brighter than Ks=24, which is five magnitudes fainter than the apparent magnitude of a L* galaxy in this cluster. The C-M relation is clearly identified as an over-density of galaxies with colours near J-Ks=1.85. The slope of the relation is -0.05 and the intrinsic scatter is 0.06 magnitudes. Both the slope and the scatter are consistent with the values measured for clusters at lower redshifts. These quantities have not evolved from z=0 to z=1.24. However, significant evolution in the mean J-Ks colour is detected. On average, the galaxies in RDCS J1252.9-2927 are 0.25 magnitudes bluer than early-type galaxies in the Coma cluster. Using instantaneous single-burst solar-metallicity models, the average age of galaxies in the center of RDCS J1252.9-2927 is 2.7 Gyrs.

We present optical and near-infrared observations of the dim afterglow of GRB 020124, obtained between 2 and 68 hours after the gamma-ray burst. The burst occurred in a very faint (R > 29.5) Damped Ly-alpha Absorber (DLA) at a redshift of z = 3.198 +- 0.004. The derived column density of neutral hydrogen is log(N_H) = 21.7 +- 0.2 and the rest-frame reddening is constrained to be E(B-V) < 0.065, i.e., A_V < 0.20 for standard extinction laws with R_V ~ 3. The resulting dust-to-gas ratio is less than 11 % of that found in the Milky Way, but consistent with the SMC and high-redshift QSO DLAs, indicating a low metallicity and/or a low dust-to-metals ratio in the burst environment. A grey extinction law (large R_V), produced through preferential destruction of small dust grains by the GRB, could increase the derived A_V and dust-to-gas ratio. The dimness of the afterglow is however fully accounted for by the high redshift: If GRB 020124 had been at z = 1 it would have been approximately 1.8 mag brighter--in the range of typical bright afterglows.

We report the discovery of an overdensity of faint red galaxies in the vicinity of the z=1.786 radio galaxy 3C 294. The overdensity, discovered in a 84 min Ks-band ISAAC/VLT image is significant at the 2.4 sigma level (compared to the local field density), and overlaps with the extended X-ray emission around 3C 294 detected with the Chandra X-ray Observatory. The near-infrared colours of the galaxies making up the overdensity show a large scatter and the galaxies do not follow a red sequence in the colour magnitude diagram. If the galaxies are in a cluster at z=1.786 they must be dominated by young stellar populations with different star-formation histories.

A photometric redshift analysis of galaxies in the field of the wide-separation gravitational lens MG2016+112 reveals a population of 69 galaxies with photometric redshifts consistent with being in a cluster at the redshift of the giant elliptical lensing galaxy, z=1.00. The Ks-band luminosity function of the cluster galaxies is well represented by the Schechter function with a characteristic magnitude Ks*=18.90 (-0.57/+0.45) and faint-end slope alpha=-0.60 (-0.33/+0.39), consistent with what is expected for a passively evolving population of galaxies formed at high redshift, zf>2. From the total Ks-band flux of the cluster galaxies and a dynamical estimate of the total mass of the cluster, the restframe Ks-band mass-to-light ratio of the cluster is derived to be M/L=27(-17/+64)(h_50)(M/L)_{sun}, in agreement with the upper limit derived earlier from Chandra X-ray observations and the value derived locally in the Coma cluster. The cluster galaxies span a red sequence with a considerable scatter in the colour-magnitude diagrams, suggesting that they contain young stellar populations in addition to the old populations of main-sequence stars that dominate the Ks-band luminosity function. This is in agreement with spectroscopic observations which show that 5 out of the 6 galaxies in the field confirmed to be at the redshift of the lensing galaxy have emission lines. The projected spatial distribution of the cluster galaxies is filamentary-like rather than centrally concentrated around the lensing galaxy, and show no apparent luminosity segregation. A handful of the cluster galaxies show evidence of merging/interaction. The results presented in this paper suggest that a young cluster of galaxies is assembling around MG2016+112.

The X-ray properties of the haloes of disk galaxies formed in fully cosmological, hydro/gravity simulations are discussed. The results are found to be consistent with observational X-ray detections and upper limits. Disk galaxy haloes are predicted to be about an order of magnitude brighter in soft X-rays at z~1 than at z=0. The Ly-alpha and H-alpha surface brightness of an edge-on, Milky Way like model galaxy has been determined. The emission is found to be quite extended, with a scale height of about 600 pc, neglecting extinction corrections.

Bolometric and 0.2-2 keV X-ray luminosities of the hot gas haloes of simulated disc galaxies have been calculated at redshift z=0. The TreeSPH simulations are fully cosmological and the sample of 44 disc galaxies span a range in characteristic circular speeds of V_c = 130-325 km/s. The galaxies have been obtained in simulations with a considerable range of physical parameters, varying the baryonic fraction, the gas metallicity, the meta-galactic UV field, the cosmology, the dark matter type, and also the numerical resolution. The models are found to be in agreement with the (few) relevant X-ray observations available at present. The amount of hot gas in the haloes is also consistent with constraints from pulsar dispersion measures in the Milky Way. Forthcoming XMM and Chandra observations should enable much more stringent tests and provide constraints on the physical parameters. We find that simple cooling flow models over-predict X-ray luminosities by up to two orders of magnitude for high (but still realistic) cooling efficiencies relative to the models presented here. Our results display a clear trend that increasing cooling efficiency leads to decreasing X-ray luminosities at z=0. The reason is found to be that increased cooling efficiency leads to a decreased fraction of hot gas relative to total baryonic mass inside of the virial radius at present. At gas metal abundances of a third solar this hot gas fraction becomes as low as just a few percent. We also find that most of the X-ray emission comes from the inner parts (inner about 20 kpc) of the hot galactic haloes. Finally, we find for realistic choices of the physical parameters that disc galaxy haloes possibly were more than one order of magnitude brighter in soft X-ray emission at z=1, than at present.

Observations of the extraordinarily bright optical afterglow (OA) of GRB 991208 started 2.1 d after the event. The flux decay constant of the OA in the R-band is -2.30 +/- 0.07 up to 5 d, which is very likely due to the jet effect, and after that it is followed by a much steeper decay with constant -3.2 +/- 0.2, the fastest one ever seen in a GRB OA. A negative detection in several all-sky films taken simultaneously to the event implies either a previous additional break prior to 2 d after the occurrence of the GRB (as expected from the jet effect). The existence of a second break might indicate a steepening in the electron spectrum or the superposition of two events. Once the afterglow emission vanished, contribution of a bright underlying SN is found, but the light curve is not sufficiently well sampled to rule out a dust echo explanation. Our determination of z = 0.706 indicates that GRB 991208 is at 3.7 Gpc, implying an isotropic energy release of 1.15 x 10E53 erg which may be relaxed by beaming by a factor > 100. Precise astrometry indicates that the GRB coincides within 0.2" with the host galaxy, thus given support to a massive star origin. The absolute magnitude is M_B = -18.2, well below the knee of the galaxy luminosity function and we derive a star-forming rate of 11.5 +/- 7.1 Mo/yr. The quasi-simultaneous broad-band photometric spectral energy distribution of the afterglow is determined 3.5 day after the burst (Dec 12.0) implying a cooling frequency below the optical band, i.e. supporting a jet model with p = -2.30 as the index of the power-law electron distribution.

We present UBVRIz' photometry of the gravitational lens candidate CLASS B1152+119 obtained with the Nordic Optical Telescope. The two QSO components are resolved in the B, V, R, I and z' bands confirming the lensing nature of the system. The z=0.44 lens galaxy is clearly detected in B, R, I and z' and its position is found to be almost coincident with the faint QSO image which is heavily extincted (relative to the brighter QSO image) by dust in the lens galaxy. The extinction curve of the lens galaxy derived from the relative photometry is well fitted by a Galactic extinction law with 1.3 < R_V < 2.0 and E(B-V) ~ 1. From a simple model of the system we predict a time delay of ~ 60 days.