We present details of the Automated Radio Telescope Imaging Pipeline (ARTIP) and results of a sensitive blind search for HI and OH absorbers at $z<0.4$ and $z<0.7$, respectively. ARTIP is written in Python 3.6, extensively uses the Common Astronomy Software Application (CASA) tools and tasks, and is designed to enable the geographically-distributed MeerKAT Absorption Line Survey (MALS) team to collaboratively process large volumes of radio interferometric data. We apply it to the first MALS dataset obtained using the 64-dish MeerKAT radio telescope and 32K channel mode of the correlator. With merely 40 minutes on target, we present the most sensitive spectrum of PKS1830-211 ever obtained and characterize the known HI ($z=0.19$) and OH ($z=0.89$) absorbers. We further demonstrate ARTIP's capabilities to handle realistic observing scenarios by applying it to a sample of 72 bright radio sources observed with the upgraded Giant Metrewave Radio Telescope (uGMRT) to blindly search for HI and OH absorbers. We estimate the numbers of HI and OH absorbers per unit redshift to be $n_{21}(z\sim0.18)<$0.14 and $n_{\rm OH}(z\sim0.40)<$0.12, respectively, and constrain the cold gas covering factor of galaxies at large impact parameters (50 kpc $<\rho<$ 150 kpc) to be less than 0.022. Due to the small redshift path, $\Delta z\sim$13 for HI with column density$>5.4\times10^{19}$ cm$^{-2}$, the survey has probed only the outskirts of star-forming galaxies at $\rho>30$ kpc. MALS with the expected $\Delta z\sim10^{3-4}$ will overcome this limitation and provide stringent constraints on the cold gas fraction of galaxies in diverse environments over $0<z<1.5$.

We investigate how damped Lyman-$\alpha$ absorbers (DLAs) at z ~ 2-3, detected in large optical spectroscopic surveys of quasars, trace the population of star-forming galaxies. Building on previous results, we construct a model based on observed and physically motivated scaling relations in order to reproduce the bivariate distributions of metallicity, Z, and HI column density, N(HI). Furthermore, the observed impact parameters for galaxies associated to DLAs are in agreement with the model predictions. The model strongly favours a metallicity gradient, which scales with the luminosity of the host galaxy, with a value of $\gamma$* = -0.019 $\pm$ 0.008 dex kpc$^{-1}$ for L* galaxies that gets steeper for fainter galaxies. We find that DLAs trace galaxies over a wide range of galaxy luminosities, however, the bulk of the DLA cross-section arises in galaxies with L ~ 0.1 L* at z ~ 2.5 broadly consistent with numerical simulations.

We report on the discovery of a peculiar Broad Absorption Line (BAL) quasar identified in our Gaia-assisted survey of red quasars. The systemic redshift of this quasar was difficult to establish due to the absence of conspicuous emission lines. Based on deep and broad BAL troughs (at least SiIV, CIV, and AlIII), a redshift of z=2.41 was established under the assumption that the systemic redshift can be inferred from the red edge of the BAL troughs. However, we observe a weak and spatially-extended emission line at 4450 AA most likely due to Lyman-alpha emission, which implies a systemic redshift of z=2.66 if correctly identified. There is also evidence for the onset of Lyman-alpha forest absorption bluewards of 4450 AA and evidence for H-alpha emission in the K-band consistent with a systemic redshift of z=2.66. If this redshift is correct, the quasar is an extreme example of a detached low-ionisation BAL quasar. The BAL lines must originate from material moving with very large velocities ranging from 22000 to 40000 km/s. To our knowledge, this is the first case of a systemic-redshift measurement based on extended Lyman-$\alpha$ emission for a BAL quasar, a method that should also be useful in cases of sufficiently distant BL Lac quasars without systemic-redshift information.

We present results of the MUSE-ALMA Halos, an ongoing study of the Circum-Galactic Medium (CGM) of low redshift galaxies (z < 1.4), currently comprising 14 strong HI absorbers in five quasar fields. We detect 43 galaxies associated with absorbers down to star formation rate (SFR) limits of 0.01-0.1 solar masses/yr, found within impact parameters (b) of 250 kpc from the quasar sightline. Excluding the targeted absorbers, we report a high detection rate of 89 per cent and find that most absorption systems are associated with pairs or groups of galaxies (three to eleven members). We note that galaxies with the smallest impact parameters are not necessarily the closest to the absorbing gas in velocity space. Using a multi-wavelength dataset (UVES/HIRES, HST, MUSE), we combine metal and HI column densities, allowing for derivation of the lower limits of neutral gas metallicity as well as emission line diagnostics (SFR, metallicities) of the ionised gas in the galaxies. We find that groups of associated galaxies follow the canonical relations of N(HI) -- b and W_r(2796) -- b, defining a region in parameter space below which no absorbers are detected. The metallicity of the ISM of associated galaxies, when measured, is higher than the metallicity limits of the absorber. In summary, our findings suggest that the physical properties of the CGM of complex group environments would benefit from associating the kinematics of individual absorbing components with each galaxy member.

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

We present the results from VLT/X-shooter spectroscopic observations of 11 extremely strong intervening damped Lyman-alpha absorbers (ESDLAs) initially selected as high N(Hi) (i.e.>=5x10^21 cm-2) candidates from the Sloan Digital Sky Survey (SDSS). We confirm the high Hi column densities which we measure to be in the range log N(Hi) = 21.6-22.4. Molecular hydrogen is detected with high column densities (N(H_2)>=10^18 cm-2) in five out of eleven systems, three of which are reported here for the first time. We compare the chemical properties of this sample of ESDLAs, supplemented with literature measurements, to that of DLAs located at the redshift of long-duration gamma-ray bursts (GRB-DLAs). We confirm that the two populations are almost indistinguishable in terms of chemical enrichment, H_2 column density and gas kinematics. All this suggests that ESDLAs and GRB-DLAs probe similar galactic environments. We search for the galaxy counterparts of ESDLAs and find associated emission lines in three out of eleven systems, two of which are reported here for the first time (towards the quasars SDSS J002503.03+114547.80 and SDSS J114347.21+142021.60, respectively). The measured separations between the quasar sightlines and the emission associated with the ESDLA galaxy are all very small (rho < 3 kpc). While the small impact parameters are similar to what is observed for GRB-DLAs, the associated star-formation rates are on average lower than seen for GRB host galaxies. This is explained by long-duration GRBs being associated with the death of massive stars, hence pinpointing regions of active star formation in the GRB host galaxies. Our observations support the suggestion from the literature that ESDLAs could act as blind analogues of GRB-DLAs, probing high column density neutral gas in the heart of high-redshift galaxies, without any prior on the instantaneous star-formation rate.

We present results from spectroscopic observations with X-shooter at the Very Large Telescope of seven H2-bearing DLAs at high redshifts (z$_{\rm abs}\sim 2.5-3$). These DLAs were originally selected from the presence of strong H$_2$ lines directly seen at the DLA redshift in low-resolution, low S/N SDSS spectra. We confirm the detection of molecular hydrogen in all of them. We measure the column densities of HI, H$_2$ in various rotational levels, and metal species, and associated dust extinction. The metallicities, obtained from undepleted species, are in the range logZ=-0.8 to -0.2. We discuss the chemical enrichment in these clouds and compare their properties with that of other molecular-rich systems selected by other means. In particular, we show that three different methods of pre-selection of H$_2$-bearing DLAs in the SDSS have their own biases but complement each other mostly in terms of chemical enrichment. We use the rotational excitation of H$_2$ molecules together with the fine-structure energy levels of neutral carbon to constrain the physical conditions in the gas with the help of numerical modeling as well as analytical expressions for the surface density at which atomic to molecular conversion happens. We find that the H$_2$-bearing medium revealed by the studied DLAs has typical values for the kinetic temperature, hydrogen density, and UV radiation field of, respectively, T$\sim$100K, $n_{\rm H}\sim 100$ cm$^{-3}$, and I$_{\rm UV}$ about twice the intensity of the Draine field. Detailed studies combining different selections should, therefore, bring important clues to understand the HI-H$_2$ transition at high redshift.

We report the detections of molecular hydrogen (H$_2$), vibrationally-excited H$_2$ (H$^*_2$), and neutral atomic carbon (CI), in two new afterglow spectra of GRBs\,181020A ($z=2.938$) and 190114A ($z=3.376$), observed with X-shooter at the Very Large Telescope (VLT). Both host-galaxy absorption systems are characterized by strong damped Lyman-$\alpha$ absorbers (DLAs) and substantial amounts of molecular hydrogen with $\log N$(HI, H$_2$) = $22.20\pm 0.05,~20.40\pm 0.04$ (GRB\,181020A) and $\log N$(HI, H$_2$) = $22.15\pm 0.05,~19.44\pm 0.04$ (GRB\,190114A). The DLA metallicites, depletion levels and dust extinctions are [Zn/H] = $-1.57\pm 0.06$, [Zn/Fe] = $0.67\pm 0.03$, and $A_V = 0.27\pm 0.02$\,mag (GRB\,181020A) and [Zn/H] = $-1.23\pm 0.07$, [Zn/Fe] = $1.06\pm 0.08$, and $A_V = 0.36\pm 0.02$\,mag (GRB\,190114A). We then examine the molecular gas content of all known H$_2$-bearing GRB-DLAs and explore the physical conditions and characteristics of these systems. We confirm that H$_2$ is detected in all CI- and H$^*_2$-bearing GRB absorption systems, but that these rarer features are not necessarily detected in all GRB H$_2$ absorbers. We find that a large molecular fraction of $f_{\rm H_2} \gtrsim 10^{-3}$ is required for CI to be detected. The defining characteristic for H$^*_2$ to be present is less clear, though a large H$_2$ column density is an essential factor. We then derive the H$_2$ excitation temperatures of the molecular gas and find that they are relatively low with $T_{\rm ex} \approx 100 - 300$\,K, however, there could be evidence of warmer components populating the high-$J$ H$_2$ levels in GRBs\,181020A and 190114A. Finally, we demonstrate that the otherwise successful X-shooter GRB afterglow campaign is hampered by a significant dust bias excluding the most dust-obscured H$_2$ absorbers from identification [Abridged].

We present results from a search for strong H2 absorption systems proximate to quasars (zabs~zem) in the Sloan Digital Sky Survey (SDSS) Data Release 14. The search is based on the Lyman-Werner band signature of damped H2 absorption lines without any prior on the associated metal or neutral hydrogen content. This has resulted in the detection of 81 systems with log N(H2)~19-20 located within a few thousand km/s from the quasar. Compared to a control sample of intervening systems, this implies an excess of proximate H2 systems by about a factor of 4 to 5. The incidence of H2 systems increases steeply with decreasing relative velocity, reaching an order of magnitude higher than expected from intervening statistics at Delta_v<1000 km/s. The most striking feature of the proximate systems compared to the intervening ones is the presence of Ly-alpha emission in the core of the associated damped HI absorption line in about half of the sample. This puts constraints on the relative projected sizes of the absorbing clouds to those of the quasar line emitting regions. Using the SDSS spectra, we estimate the HI, metal and dust content of the systems, which are found to have typical metallicities of one tenth Solar, albeit with a large spread among individual systems. We observe trends between the fraction of leaking Ly-alpha emission and the relative absorber-quasar velocity as well as with the excitation of several metal species, similar to what has been seen in metal-selected proximate DLAs. With the help of theoretical HI-H2 transition relations, we show that the presence of H2 helps to break the degeneracy between density and strength of the UV field as main sources of excitation and hence provides unique constraints on the possible origin and location of the absorbing clouds. We suggest that most of these systems originate from galaxies in the quasar group. [truncated]

We present a systematic study of the impact of a dust bias on samples of damped Lyman-$\alpha$ absorbers (DLAs). This bias arises as an effect of the magnitude and colour criteria utilized in the Sloan Digital Sky Survey (SDSS) quasar target selection up until data release 7 (DR7). The bias has previously been quantified assuming only a contribution from the dust obscuration. In this work, we apply the full set of magnitude and colour criteria used up until SDSS DR7 in order to quantify the full impact of dust biasing against dusty and metal-rich DLAs. We apply the quasar target selection algorithm on a modelled population of intrinsic colours, and by exploring the parameter space consisting of redshift, ($z_{\rm QSO}$ and $z_{\rm abs}$), optical extinction, and HI column density, we demonstrate how the selection probability depends on these variables. We quantify the dust bias on the following properties derived for DLAs at z=3: the incidence rate, the mass density of neutral hydrogen and metals, the average metallicity. We find that all quantities are significantly affected. When considering all uncertainties, the mass density of neutral hydrogen is underestimated by 10 to 50%, and the mass density in metals is underestimated by 30 to 200%. Lastly, we find that the bias depends on redshift. At redshift z=2.2, the mass density of neutral hydrogen and metals might be underestimated by up to a factor of 2 and 5, respectively. Characterizing such a bias is crucial in order to accurately interpret and model the properties and metallicity evolution of absorption-selected galaxies.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

Damped Lyman-{\alpha} Absorbers (DLAs) as a class of QSO absorption-line systems are currently our most important source of detailed information on the cosmic chemical evolution of galaxies. However, the degree to which this information is biased by dust remains to be understood. One strategy is to specifically search for QSOs reddened by metal-rich and dusty foreground absorbers. In this Letter we present the discovery of a z=2.60 QSO strongly reddened by dust in an intervening extremely-strong DLA at z=2.226. This QSO was identified through a novel selection combining the astrometric measurements from ESA's Gaia satellite with extent optical and near/mid-infrared photometry. We infer a total neutral atomic-hydrogen column density of log N(HI)=21.95{\pm}0.15 and a lower limit on the gas-phase metallicity of [Zn/H]>-0.96. This DLA is also remarkable in exhibiting shielded neutral gas witnessed in CI and tentative detections of CO molecular bands. The Spectral Energy Distribution (SED) of the QSO is well-accounted for by a normal QSO-SED reddened by dust from a DLA with a 10%-of-Solar metallicity, dust extinction of A_V=0.82{\pm}0.02mag, and LMC-like extinction curve including the characteristic 2175{\AA} extinction feature. Such QSO absorption-line systems have shown to be very rare in previous surveys, which have mostly revealed sight-lines with low extinction. The present case therefore suggests that previous samples have under-represented the fraction of dusty absorbers. Building a complete sample of such systems is needed to assess the significance of this effect.

We have designed the Extremely Luminous Quasar Survey (ELQS) to provide a highly complete census of unobscured UV-bright quasars during the cosmic noon, $z=2.8-5.0$. Here we report the discovery of 70 new quasars in the ELQS South Galactic Cap (ELQS-S) quasar sample, doubling the number of known extremely luminous quasars in $4,237.3\,\rm{deg}^2$ of the SDSS footprint. These observations conclude the ELQS and we present the properties of the full ELQS quasar catalog, containing 407 quasars over $11,838.5\,\rm{deg}^2$. Our novel ELQS quasar selection strategy resulted in unprecedented completeness at the bright end and allowed us to discover 109 new quasars in total. This marks an increase of $\sim36\%$ (109/298) to the known population at these redshifts and magnitudes, while we further are able to retain a selection efficiency of $\sim80\%$. On the basis of 166 quasars from the full ELQS quasar catalog, who adhere to the uniform criteria of the 2MASS point source catalog, we measure the bright-end quasar luminosity function (QLF) and extend it one magnitude brighter than previous studies. Assuming a single power law with exponential density evolution for the functional form of the QLF, we retrieve the best fit parameters from a maximum likelihood analysis. We find a steep bright-end slope of $\beta\approx-4.1$ and we can constrain the bright-end slope to $\beta\leq-3.4$ with $99\%$ confidence. The density is well modeled by the exponential redshift evolution, resulting in a moderate decrease with redshift ($\gamma\approx-0.4$).

We present a survey for neutral atomic-carbon (CI) along gamma-ray burst (GRB) sightlines, which probes the shielded neutral gas-phase in the interstellar medium (ISM) of GRB host galaxies at high redshift. We compile a sample of 29 medium- to high-resolution GRB optical afterglow spectra spanning a redshift range through most of cosmic time from $1 < z < 6$. We find that seven ($\approx 25\%$) of the GRBs entering our statistical sample have CI detected in absorption. It is evident that there is a strong excess of cold gas in GRB hosts compared to absorbers in quasar sightlines. We investigate the dust properties of the GRB CI absorbers and find that the amount of neutral carbon is positively correlated with the visual extinction, $A_V$, and the strength of the 2175 \AA\ dust extinction feature, $A_{\mathrm{bump}}$. GRBs with CI detected in absorption are all observed above a certain threshold of $\log N$(HI)$/\mathrm{cm}^{-2}$ + [X/H] > 20.7 and a dust-phase iron column density of $\log N$(Fe)$_{\mathrm{dust}}/\mathrm{cm}^{-2}$ > 16.2. In contrast to the SED-derived dust properties, the strength of the CI absorption does not correlate with the depletion-derived dust properties. This indicates that the GRB CI absorbers trace dusty systems where the dust composition is dominated by carbon-rich dust grains. The observed higher metal and dust column densities of the GRB CI absorbers compared to H$_2$- and CI-bearing quasar absorbers is mainly a consequence of how the two absorber populations are selected, but is also required in the presence of intense UV radiation fields in actively star-forming galaxies.

We present a study of cold gas absorption from a damped Lyman-$\alpha$ absorber (DLA) at redshift $z_{\rm abs}=1.946$ towards two lensed images of the quasar J144254.78+405535.5 at redshift $z_{\rm QSO} = 2.590$. The physical separation of the two lines of sight at the absorber redshift is $d_{\rm abs}=0.7$~kpc based on our lens model. We observe absorption lines from neutral carbon and H$_2$ along both lines of sight indicating that cold gas is present on scales larger than $d_{\rm abs}$. We measure column densities of HI to be $\log N(\rm H\,\i) = 20.27\pm0.02$ and $20.34\pm0.05$ and of H$_2$ to be $\log N(\rm H_2) = 19.7\pm0.1$ and $19.9\pm0.2$. The metallicity inferred from sulphur is consistent with Solar metallicity for both sightlines: $[{\rm S/H}]_A = 0.0\pm0.1$ and $[{\rm S/H}]_B = -0.1\pm0.1$. Based on the excitation of low rotational levels of H$_2$, we constrain the temperature of the cold gas phase to be $T=109\pm20$ and $T=89\pm25$ K for the two lines of sight. From the relative excitation of fine-structure levels of CI, we constrain the hydrogen volumetric densities in the range of $40-110$ cm$^{-3}$. Based on the ratio of observed column density and volumetric density, we infer the average individual `cloud' size along the line of sight to be $l\approx0.1$ pc. Using the transverse line-of-sight separation of 0.7 kpc together with the individual cloud size, we are able to put an upper limit to the volume filling factor of cold gas of $f_{\rm vol} < 0.2$ %. Nonetheless, the projected covering fraction of cold gas must be large (close to unity) over scales of a few kpc in order to explain the presence of cold gas in both lines of sight. Compared to the typical extent of DLAs (~10-30 kpc), this is consistent with the relative incidence rate of CI absorbers and DLAs.

Projected quasar galaxy pairs provide powerful means to study the circumgalactic medium (CGM) that maintains the relics of galactic feedback and the accreted gas from the intergalactic medium. Here, we study the nature of a Lyman Limit system (LLS) with N(HI)=10$^{19.1\pm0.3}$ cm$^{-2}$ and a dust-uncorrected metallicity of [Fe/H]$=-1.1\pm0.3$ at $z=0.78$ towards Q0152$-020$. The MgII absorption profiles are composed of a main saturated and a few weaker optically thin components. Using MUSE observations we detect one galaxy close to the absorption redshift at an impact parameter of 54 kpc. This galaxy exhibits nebular emission lines from which we measure a dust-corrected star formation rate of $10^{+8}_{-5}$ M$_\odot$ yr$^{-1}$ and an emission metallicity of [O/H]$=-0.1\pm0.2$. By combining the absorption line kinematics with the host galaxy morphokinematics we find that while the main absorption component can originate from a galactic wind at $V_{\rm w}=110\pm4$ km s$^{-1}$ the weaker components cannot. We estimate a mass ejection rate of $\dot M\gtrsim0.8$ M$_\odot$ yr$^{-1}$ that translates to a loading factor of $\eta\gtrsim0.1$. Since the local escape velocity of the halo, $V_{\rm esc}\simeq430$ km s$^{-1}$, is a few times larger than $V_{\rm w}$, we expect this gas will remain bound to the host galaxy. These observations provide additional constraints on the physical properties of winds predicted by galaxy formation models. We also present the VLT/X-Shooter data analysis of 4 other absorbing systems at $1.1<z<1.5$ in this sightline with their host galaxies identified in the MUSE data.

We present the detection analysis of a diffuse molecular cloud at z$_{abs}$=2.4636 towards the quasar SDSS J1513+0352(z$_{em}\,\simeq$ 2.68) observed with the X-shooter spectrograph(VLT). We measure very high column densities of atomic and molecular hydrogen, with log N(HI,H$_2$)$\simeq$21.8,21.3. This is the highest H$_2$ column density ever measured in an intervening damped Lyman-alpha system but we do not detect CO, implying log N(CO)/N(H$_2$) < -7.8, which could be due to a low metallicity of the cloud. From the metal absorption lines, we derive the metallicity to be Z $\simeq$ 0.15 Z$_{\odot}$ and determine the amount of dust by measuring the induced extinction of the background quasar light, A$_V$ $\simeq$ 0.4. We also detect Ly-$\alpha$ emission at the same redshift, with a centroid located at a most probable impact parameter of only $\rho\,\simeq$ 1.4 kpc. We argue that the line of sight is therefore likely passing through the ISM of a galaxy as opposed to the CGM. The relation between the surface density of gas and that of star formation seems to follow the global empirical relation derived in the nearby Universe although our constraints on the star formation rate and on the galaxy extent remain too loose to be conclusive. We study the transition from atomic to molecular hydrogen using a theoretical description based on the microphysics of molecular hydrogen. We use the derived chemical properties of the cloud and physical conditions (T$_k\,\simeq$90 K and n$\simeq$250 cm$^{-3}$ derived through the excitation of H$_2$ rotational levels and neutral carbon fine structure transitions to constrain the fundamental parameters that govern this transition. By comparing the theoretical and observed HI column densities, we are able to bring an independent constraint on the incident UV flux, which we find to be in agreement with that estimated from the observed star formation rate.

We present here a survey of high-ionization absorption lines in the afterglow spectra of long-duration gamma-ray bursts (GRBs) obtained with the VLT/X-shooter spectrograph. Our main goal is to investigate the circumburst medium in the natal regions of GRBs. Our primary focus is on the NV 1238,1242 line transitions, but we also discuss other high-ionization lines such as OVI, CIV and SiIV. We find no correlation between the column density of NV and the neutral gas properties such as metallicity, HI column density and dust depletion, however the relative velocity of NV, typically a blueshift with respect to the neutral gas, is found to be correlated with the column density of HI. This may be explained if the NV gas is part of an HII region hosting the GRB, where the region's expansion is confined by dense, neutral gas in the GRB's host galaxy. We find tentative evidence (at 2-sigma significance) that the X-ray derived column density, N_H,X, may be correlated with the column density of NV, which would indicate that both measurements are sensitive to the column density of the gas located in the vicinity of the GRB. We investigate the scenario where NV (and also OVI) is produced by recombination after the corresponding atoms have been stripped entirely of their electrons by the initial prompt emission, in contrast to previous models where highly-ionized gas is produced by photoionization from the GRB afterglow.

We here apply a novel technique selecting quasar candidates purely as sources with zero proper motions in the Gaia data release 2 (DR2). We demonstrate that this approach is highly efficient toward high Galactic latitudes with < 25% contamination from stellar sources. Such a selection technique offers a very pure sample completeness, since all cosmological point sources are selected regardless of their intrinsic spectral properties within the limiting magnitude of Gaia. We carry out a pilot-study by defining a sample compiled by including all Gaia-DR2 sources within one degree of the North Galactic Pole (NGP) selected to have proper motions consistent with zero within 2-sigma uncertainty. By cross-matching the sample to the optical Sloan Digital Sky Survey (SDSS) and the mid-infrared AllWISE photometric catalogues we investigate the colours of each of our sources. Together with already spectroscopically confirmed quasars we are therefore able to determine the efficiency of our selection. The majority of the zero proper motion sources have optical to mid-infrared colours consistent with known quasars. The remaining population may be contaminating stellar sources, but some may also be quasars with colours similar to stars. Spectroscopic follow-up of the zero proper motion sources is needed to unveil such a hitherto hidden quasar population. This approach has the potential to allow substantial progress on many important questions concerning quasars such as determining the fraction of dust-obscured quasars, the fraction of broad absorption line (BAL) quasars, and the metallicity distribution of damped Lyman-$\alpha$ absorbers. The technique could also potentially reveal new types of quasars or even new classes of cosmological point sources.

The cosmic chemical enrichment as measured from damped Ly$\alpha$ absorbers (DLAs) will be underestimated if dusty and metal-rich absorbers have evaded identification. Here we report the discovery and present the spectroscopic observations of a quasar, KV-RQ\,1500-0031, at $z=2.520$ reddened by a likely dusty DLA at $z=2.428$ and a strong MgII absorber at $z=1.603$. This quasar was identified as part of the KiDS-VIKING Red Quasar (KV-RQ) survey, specifically aimed at targeting dusty absorbers which may cause the background quasars to escape the optical selection of e.g. the SDSS quasar survey. For the DLA we find an HI column density of $\log N$(HI) = $21.2\pm 0.1$ and a metallicity of [X/H] = $-0.90\pm 0.20$ derived from an empirical relation based on the equivalent width of SiII$\lambda$1526. We observe a total visual extinction of $A_V=0.16$ mag induced by both absorbers. We compile a sample of 17 additional dusty ($A_V > 0.1$ mag) DLAs toward quasars (QSO-DLAs) from the literature for which we characterize the properties of HI column density, metallicity and dust. From this sample we also estimate a correction factor to the overall DLA metallicity budget. We demonstrate that the dusty QSO-DLAs have high metal column densities ($\log N$(HI) + [X/H]) and are more similar to gamma-ray burst (GRB)-selected DLAs (GRB-DLAs) than regular QSO-DLAs. We evaluate the effect of dust reddening in DLAs as well as illustrate how the induced color excess of the underlying quasars can be significant (up to $\sim 1$ mag in various optical bands), even for low to moderate extinction values ($A_V \lesssim 0.6$ mag). Finally we discuss the direct and indirect implications of a significant dust bias in both QSO- and GRB-DLA samples. [Abridged]

We study a sample of 17 z>1.5 absorbers selected based on the presence of strong CI absorption lines in SDSS spectra and observed with the ESO-VLT spectrograph X-shooter. We derive metallicities, depletion onto dust, and extinction by dust, and analyse the absorption from MgII, MgI, CaII and NaI that are redshifted into the near infrared wavelength range. We show that most of these CI absorbers have high metallicity and dust content. We detect nine CaII absorptions with $W$(CaII$\lambda$3934) >0.23 \AA out of 14 systems where we have appropriate wavelength coverage. The observed equivalent widths are similar to what has been measured in other lower redshift surveys of CaII systems. We detect ten NaI absorptions in the 11 systems where we could observe this absorption. The median equivalent width ($W$(NaI$\lambda$5891) = 0.68 \AA) is larger than what is observed in local clouds with similar HI column densities but also in z<0.7 CaII systems detected in the SDSS. The systematic presence of NaI absorption in these CI systems strongly suggests that the gas is neutral and cold, maybe part of the diffuse molecular gas in the ISM of high-redshift galaxies. Most of the systems (12 out of 17) have $W$(MgII$\lambda$2796) > 2.5 \AA and six of them have log N(HI) < 20.3, with the extreme case of J1341+1852 that has log N(HI) = 18.18. The MgII absorptions are spread over more than $\Delta v$ $\sim$ 400 km s$^{-1}$ for half of the systems; three absorbers have $\Delta v$ > 500 km s$^{-1}$. The kinematics are strongly perturbed for most of these systems, which probably do not arise in quiet disks and must be close to regions with intense star-formation activity and/or are part of interacting objects. All this suggests that a large fraction of the cold gas at high redshift arises in disturbed environments.

I present a Python package developed for fitting Voigt profiles to absorption lines. The software fits multiple components for various atomic lines simultaneously allowing parameters to be tied and fixed. Moreover, the code is able to automatically fit a polynomial continuum model together with the line profiles. Lastly, a physical model can readily be used to constrain thermal and turbulent broadening of absorption lines as well as implementing molecular excitation models. The code can be run with interactive features such as manual continuum placement locally around each line, manual masking of undesired fitting regions, and interactive definition of velocity components for various elements. This greatly improves the ease by which the initial guesses can be estimated. Since the code is written in pure Python, it can easily be scripted and modified to fit the user's needs. The code uses a $\chi^2$ minimization approach to find the best solution. The code and a set of test-data together with the full documentation is available on GitHub.

We present a preparatory spectroscopic survey to identify radio-bright, high-redshift quasars for the MeerKAT Absorption Line Survey (MALS). The candidates have been selected on the basis of a single flux density limit at 1.4 GHz (>200 mJy) together with mid-infrared color criteria from the Wide-field Infrared Survey Explorer (WISE). Through spectroscopic observations using the Nordic Optical Telescope, we identify 72 quasars out of 99 candidates targeted. We measure the spectroscopic redshifts based on characteristic, broad emission lines present in the spectra. Of these 72 quasars, 64 and 48 objects are at sufficiently high redshift (z>0.6 and z>1.4) to be used for the L-band and UHF-band spectroscopic follow-up with the Square Kilometre Array (SKA) precursor in South Africa: the MeerKAT.

Deep galaxy surveys have revealed that the global star formation rate (SFR) density in the Universe peaks at 1 < z < 2 and sharply declines towards z = 0. But a clear picture of the underlying processes, in particular the evolution of cold atomic (~100 K) and molecular gas phases, that drive such a strong evolution is yet to emerge. MALS is designed to use MeerKAT's L- and UHF-band receivers to carry out the most sensitive (N(HI)>10$^{19}$ cm$^{-2}$) dust-unbiased search of intervening HI 21-cm and OH 18-cm absorption lines at 0 < z < 2. This will provide reliable measurements of the evolution of cold atomic and molecular gas cross-sections of galaxies, and unravel the processes driving the steep evolution in the SFR density. The large sample of HI and OH absorbers obtained from the survey will (i) lead to tightest constraints on the fundamental constants of physics, and (ii) be ideally suited to probe the evolution of magnetic fields in disks of galaxies via Zeeman Splitting or Rotation Measure synthesis. The survey will also provide an unbiased census of HI and OH absorbers, i.e. cold gas associated with powerful AGNs (>10$^{24}$ W Hz$^{-1}$) at 0 < z < 2, and will simultaneously deliver a blind HI and OH emission line survey, and radio continuum survey. Here, we describe the MALS survey design, observing plan and the science issues to be addressed under various science themes.

To fully exploit the potential of quasars as probes of cosmic chemical evolution and the internal gas dynamics of galaxies it is important to understand the selection effects behind the quasar samples and in particular if the selection criteria exclude foreground galaxies with certain properties (most importantly a high dust content). Here we present spectroscopic follow-up from the 10.4-m GTC telescope of a dust-reddened quasar, eHAQ0111+0641, from the extended High A_V Quasar (HAQ) survey. We find that the z=3.21 quasar has a foreground Damped Lyman-alpha Absorber (DLA) at z=2.027 along the line of sight. The DLA has very strong metal lines due to a moderately high metallicity (with an inferred lower limit of 25% of the solar metallicity), but a very large gas column density along the line-of-sight in its host galaxy. This discovery is further evidence that there is a dust bias affecting the census of metals, caused by the combined effect of dust obscuration and reddening, in existing samples of z>2 DLAs. The case of eHAQ0111+0641 illustrates that dust bias is not only caused by dust obscuration, but also dust reddening.

Starting from a summary of detection statistics of our recent X-shooter campaign, we review the major surveys, both space and ground based, for emission counterparts of high-redshift damped Ly$\alpha$ absorbers (DLAs) carried out since the first detection 25 years ago. We show that the detection rates of all surveys are precisely reproduced by a simple model in which the metallicity and luminosity of the galaxy associated to the DLA follow a relation of the form, ${\rm M_{UV}} = -5 \times \left(\,[{\rm M/H}] + 0.3\, \right) - 20.8$, and the DLA cross-section follows a relation of the form $\sigma_{DLA} \propto L^{0.8}$. Specifically, our spectroscopic campaign consists of 11 DLAs preselected based on their equivalent width of SiII $\lambda1526$ to have a metallicity higher than [Si/H] > -1. The targets have been observed with the X-shooter spectrograph at the Very Large Telescope to search for emission lines around the quasars. We observe a high detection rate of 64% (7/11), significantly higher than the typical $\sim$10% for random, HI-selected DLA samples. We use the aforementioned model, to simulate the results of our survey together with a range of previous surveys: spectral stacking, direct imaging (using the `double DLA' technique), long-slit spectroscopy, and integral field spectroscopy. Based on our model results, we are able to reconcile all results. Some tension is observed between model and data when looking at predictions of Ly$\alpha$ emission for individual targets. However, the object to object variations are most likely a result of the significant scatter in the underlying scaling relations as well as uncertainties in the amount of dust which affects the emission.

We present the results from our search for HI 21-cm absorption in a sample of 16 strong FeII systems ($W_{\rm r}$(MgII $\lambda2796$) $\ge1.0$ \AA\ and $W_{\rm r}$(FeII $\lambda2600$) or $W_{\rm FeII}$ $\ge1$ \AA) at $0.5<z<1.5$ using the Giant Metrewave Radio Telescope and the Green Bank Telescope. We report six new HI 21-cm absorption detections from our sample, which have increased the known number of detections in strong MgII systems at this redshift range by $\sim50$%. Combining our measurements with those in the literature, we find that the detection rate of HI 21-cm absorption increases with $W_{\rm FeII}$, being four times higher in systems with $W_{\rm FeII}$ $\ge1$ \AA\ compared to systems with $W_{\rm FeII}$ $<1$ \AA. The $N$(HI) associated with the HI 21-cm absorbers would be $\ge 2 \times 10^{20}$ cm$^{-2}$, assuming a spin temperature of $\sim500$ K (based on HI 21-cm absorption measurements of damped Lyman-$\alpha$ systems at this redshift range) and unit covering factor. We find that HI 21-cm absorption arises on an average in systems with stronger metal absorption. We also find that quasars with HI 21-cm absorption detected towards them have systematically higher $E(B-V)$ values than those which do not. Further, by comparing the velocity widths of HI 21-cm absorption lines detected in absorption- and galaxy-selected samples, we find that they show an increasing trend (significant at $3.8\sigma$) with redshift at $z<3.5$, which could imply that the absorption originates from more massive galaxy haloes at high-$z$. Increasing the number of HI 21-cm absorption detections at these redshifts is important to confirm various trends noted here with higher statistical significance.

We present the discovery of a molecular cloud at zabs=2.5255 along the line of sight to the quasar J0000+0048. We perform a detailed analysis of the absorption lines from ionic, neutral atomic and molecular species in different excitation levels, as well as the broad-band dust extinction. We find that the absorber classifies as a Damped Lyman-alpha system (DLA) with logN(HI)(cm^-2)=20.8+/-0.1. The DLA has super-Solar metallicity with a depletion pattern typical of cold gas and an overall molecular fraction ~50%. This is the highest f-value observed to date in a high-z intervening system. Most of the molecular hydrogen arises from a clearly identified narrow (b~0.7 km/s), cold component in which CO molecules are also found, with logN(CO)~15. We study the chemical and physical conditions in the cold gas. We find that the line of sight probes the gas deep after the HI-to-H2 transition in a ~4-5 pc-size cloud with volumic density nH~80 cm^-3 and temperature of only 50 K. Our model suggests that the presence of small dust grains (down to about 0.001 {\mu}m) and high cosmic ray ionisation rate (zeta_H a few times 10^-15 s^-1) are needed to explain the observed atomic and molecular abundances. The presence of small grains is also in agreement with the observed steep extinction curve that also features a 2175 A bump. The properties of this cloud are very similar to what is seen in diffuse molecular regions of the nearby Perseus complex. The high excitation temperature of CO rotational levels towards J0000+0048 betrays however the higher temperature of the cosmic microwave background. Using the derived physical conditions, we correct for a small contribution (0.3 K) of collisional excitation and obtain TCMB(z = 2.53)~9.6 K, in perfect agreement with the predicted adiabatic cooling of the Universe. [abridged]

We present the results of a new spectroscopic survey for dusty intervening absorption systems, particularly damped Ly$\alpha$ absorbers (DLAs), towards reddened quasars. The candidate quasars are selected from mid-infrared photometry from the Wide-field Infrared Survey Explorer combined with optical and near-infrared photometry. Out of 1073 candidates, we secure low-resolution spectra for 108 using the Nordic Optical Telescope on La Palma, Spain. Based on the spectra, we are able to classify 100 of the 108 targets as quasars. A large fraction (50 %) is observed to have broad absorption lines (BALs). Moreover, we find 6 quasars with strange breaks in their spectra, which are not consistent with regular dust reddening. Using template fitting we infer the amount of reddening along each line of sight ranging from A(V)$\approx$0.1 mag to 1.2 mag (assuming an SMC extinction curve). In four cases, the reddening is consistent with dust exhibiting the 2175{\AA} feature caused by an intervening absorber, and for two of these, a MgII absorption system is observed at the best-fit absorption redshift. In the rest of the cases, the reddening is most likely intrinsic to the quasar. We observe no evidence for dusty DLAs in this survey. However, the large fraction of BAL quasars hampers the detection of absorption systems. Out of the 50 non-BAL quasars only 28 have sufficiently high redshift to detect Ly$\alpha$ in absorption.

The sub-population of quasars reddened by intrinsic or intervening clouds of dust are known to be underrepresented in optical quasar surveys. By defining a complete parent sample of the brightest and spatially unresolved quasars in the COSMOS field, we quantify to which extent this sub-population is fundamental to our understanding of the true population of quasars. By using the available multiwavelength data of various surveys in the COSMOS field, we built a parent sample of 33 quasars brighter than $J=20$ mag, identified by reliable X-ray to radio wavelength selection techniques. Spectroscopic follow-up with the NOT/ALFOSC was carried out for four candidate quasars that had not been targeted previously to obtain a 100\% redshift completeness of the sample. The population of high $A_V$ quasars (HAQs), a specific sub-population of quasars selected from optical/near-infrared photometry, is found to contribute $21\%^{+9}_{-5}$ of the parent sample. The full population of bright spatially unresolved quasars represented by our parent sample consists of $39\%^{+9}_{-8}$ reddened quasars defined by having $A_V>0.1$, and $21\%^{+9}_{-5}$ of the sample having $E(B-V)>0.1$ assuming the extinction curve of the Small Magellanic Cloud. We show that the HAQ selection works well for selecting reddened quasars, but some are missed because their optical spectra are too blue to pass the $g-r$ color cut in the HAQ selection. This is either due to a low degree of dust reddening or anomalous spectra. We find that the fraction of quasars with contributing light from the host galaxy is most dominant at $z \lesssim 1$. At higher redshifts the population of spatially unresolved quasars selected by our parent sample is found to be representative of the full population at $J<20$ mag. This work quantifies the bias against reddened quasars in studies that are based solely on optical surveys.

We present the serendipitous discovery of a projected pair of quasi-stellar objects (QSOs) with an angular separation of $\Delta\theta =4.50$ arcsec. The redshifts of the two QSOs are widely different: one, our programme target, is a QSO with a spectrum consistent with being a narrow line Seyfert 1 AGN at $z=2.05$. For this target we detect Lyman-$\alpha$, \ion{C}{4}, and \ion{C}{3]}. The other QSO, which by chance was included on the spectroscopic slit, is a Type 1 QSO at a redshift of $z=1.68$, for which we detect \ion{C}{4}, \ion{C}{3]} and \ion{Mg}{2}. We compare this system to previously detected projected QSO pairs and find that only about a dozen previously known pairs have smaller angular separation.

Understanding how galaxies evolved from the early Universe through cosmic time is a fundamental part of modern astrophysics. In order to study this evolution it is important to sample the galaxies at various times in a consistent way through time. In regular luminosity selected samples, our analyses are biased towards the brightest galaxies at all times (as these are easier to observe and identify). A complementary method relies on the absorption imprint from neutral gas in galaxies, the so-called damped Ly-alpha absorbers (DLAs) seen towards distant bright objects. This thesis seeks to understand how the absorption selected galaxies relate to the emission selected galaxies by identifying the faint glow from the absorbing galaxies at redshift z~2. In the last Chapter, a study of the more evolved, massive galaxies is presented. These galaxies are observed to be a factor of 2 to 6 times smaller than local galaxies of similar masses. A new spectroscopically selected sample is presented and the increased precision of the redshifts allows a more detailed measurement of the scatter in the mass-size relation. The size evolution of massive, quiescent galaxies is modelled by a `dilution' scenario, in which progressively larger galaxies at later times are added to the population of denser galaxies, causing an increase of the mean size of the population. This model describes the evolution of both sizes and number densities very well, however, the scatter in the model increases with time, contrary to the data. It is thus concluded that a combination of `dilution' and individual growth, e.g., through mergers, is needed. For brevity, the individual chapters based on published peer-reviewed articles are omitted and a link to the given article is given instead. Principal supervisor: Johan P. U. Fynbo

We analyze the near-infrared to UV data of 16 quasars with redshifts ranging from 0.71 $<$ $z$ $<$ 2.13 to investigate dust extinction properties. The sample presented in this work is obtained from the High $A_V$ Quasar (HAQ) survey. The quasar candidates were selected from the Sloan Digital Sky Survey (SDSS) and the UKIRT Infrared Deep Sky Survey (UKIDSS), and follow-up spectroscopy was carried out at the Nordic Optical Telescope (NOT) and the New Technology Telescope (NTT). To study dust extinction curves intrinsic to the quasars, from the HAQ survey we selected 16 cases where the Small Magellanic Cloud (SMC) law could not provide a good solution to the spectral energy distributions (SEDs). We derived the extinction curves using Fitzpatrick & Massa 1986 (FM) law by comparing the observed SEDs to the combined quasar template from Vanden Berk et al. 2001 and Glikman et al. 2006. The derived extinction, $A_V$, ranges from 0.2-1.0 mag. All the individual extinction curves of our quasars are steeper ($R_V=2.2$-2.7) than that of the SMC, with a weighted mean value of $R_V=2.4$. We derive an `average quasar extinction curve' for our sample by fitting SEDs simultaneously by using the weighted mean values of the FM law parameters and a varying $R_V$. The entire sample is well fit with a single best-fit value of $R_V=2.2\pm0.2$. The `average quasar extinction curve' deviates from the steepest Milky Way and SMC extinction curves at a confidence level $\gtrsim95\%$. Such steep extinction curves suggest a significant population of silicates to produce small dust grains. Moreover, another possibility could be that the large dust grains may have been destroyed by the activity of the nearby active galactic nuclei (AGN), resulting in steep extinction curves.

Quasi-stellar object (QSO) spectral templates are important both to QSO physics and for investigations that use QSOs as probes of intervening gas and dust. However, combinations of various QSO samples obtained at different times and with different instruments so as to expand a composite and to cover a wider rest frame wavelength region may create systematic effects, and the contribution from QSO hosts may contaminate the composite. We have constructed a composite spectrum from luminous blue QSOs at 1 < z < 2.1 selected from the Sloan Digital Sky Survey (SDSS). The observations with X-shooter simultaneously cover ultraviolet (UV) to near- infrared (NIR) light, which ensures that the composite spectrum covers the full rest-frame range from Ly$\beta$ to 11350 $\AA$ without any significant host contamination. Assuming a power-law continuum for the composite we find a spectral slope of $\alpha_\lambda$ = 1.70+/-0.01, which is steeper than previously found in the literature. We attribute the differences to our broader spectral wavelength coverage, which allows us to effectively avoid fitting any regions that are affected either by strong QSO emissions lines (e.g., Balmer lines and complex [Fe II] blends) or by intrinsic host galaxy emission. Finally, we demonstrate the application of the QSO composite spectrum for evaluating the reddening in other QSOs.

We present a detailed analysis of a red quasar at z=2.32 with an intervening damped Lyman-alpha absorber (DLA) at z=2.13. Using high quality data from the X-shooter spectrograph at ESO Very Large Telescope we find that the absorber has a metallicity consistent with Solar. We observe strong C I and H$_2$ absorption indicating a cold, dense absorbing medium. Partial coverage effects are observed in the C I lines, from which we infer a covering fraction of $27 \pm 6$ % and a physical diameter of the cloud of 0.1 pc. From the covering fraction and size, we estimate the size of the background quasar's broad line region. We search for emission from the DLA counterpart in optical and near-infrared imaging. No emission is observed in the optical data. However, we see tentative evidence for a counterpart in the H and K' band images. The DLA shows high depletion (as probed by [Fe/Zn]=-1.22) indicating that significant amounts of dust must be present in the DLA. By fitting the spectrum with various dust reddened quasar templates we find a best-fitting amount of dust in the DLA of $A(V)_{\rm DLA}=0.28 \pm 0.01|_{\rm stat} \pm 0.07|_{\rm sys}$. We conclude that dust in the DLA is causing the colours of this intrinsically very luminous background quasar to appear much redder than average quasars, thereby not fulfilling the criteria for quasar identification in the Sloan Digital Sky Survey. Such chemically enriched and dusty absorbers are thus underrepresented in current samples of DLAs.

The reionisation of the Universe is thought to have ended around z~6, as inferred from spectroscopy of distant bright background sources, such as quasars (QSO) and gamma-ray burst (GRB) afterglows. Furthermore, spectroscopy of a GRB afterglow provides insight in its host galaxy, which is often too dim and distant to study otherwise. We present the high S/N VLT/X-shooter spectrum of GRB130606A at z=5.913. We aim to measure the degree of ionisation of the IGM between 5.02<z<5.84 and to study the chemical abundance pattern and dust content of its host galaxy. We measured the flux decrement due to absorption at Ly$\alpha$, $\beta$ and $\gamma$ wavelength regions. The hydrogen and metal absorption lines formed in the host galaxy were fitted with Voigt profiles to obtain column densities. Our measurements of the Ly$\alpha$-forest optical depth are consistent with previous measurements of QSOs, but have a much smaller uncertainty. The analysis of the red damping wing yields a neutral fraction $x_{HI}<0.05$ (3$\sigma$). We obtain column density measurements of several elements. The ionisation corrections due to the GRB is estimated to be negligible (<0.03 dex), but larger corrections may apply due to the pre-existing radiation field (up to 0.4 dex based on sub-DLA studies). Our measurements confirm that the Universe is already predominantly ionised over the redshift range probed in this work, but was slightly more neutral at z>5.6. GRBs are useful probes of the ionisation state of the IGM in the early Universe, but because of internal scatter we need a larger statistical sample to draw robust conclusions. The high [Si/Fe] in the host can be due to dust depletion, alpha-element enhancement, or a combination of both. The very high value of [Al/Fe]=2.40+/-0.78 might connected to the stellar population history. We estimate the host metallicity to be -1.7<[M/H]<-0.9 (2%-13% of solar). (trunc.)

Quasi-stellar objects (QSOs) whose spectral energy distributions (SEDs) are reddened by dust either in their host galaxies or in intervening absorber galaxies are to a large degree missed by optical color selection criteria like the one used by the Sloan Digital Sky Survey (SDSS). To overcome this bias against red QSOs, we employ a combined optical and near-infrared color selection. In this paper, we present a spectroscopic follow-up campaign of a sample of red candidate QSOs which were selected from the SDSS and the UKIRT Infrared Deep Sky Survey (UKIDSS). The spectroscopic data and SDSS/UKIDSS photometry are supplemented by mid-infrared photometry from the Wide-field Infrared Survey Explorer. In our sample of 159 candidates, 154 (97%) are confirmed to be QSOs. We use a statistical algorithm to identify sightlines with plausible intervening absorption systems and identify nine such cases assuming dust in the absorber similar to Large Magellanic Cloud sightlines. We find absorption systems toward 30 QSOs, 2 of which are consistent with the best-fit absorber redshift from the statistical modeling. Furthermore, we observe a broad range in SED properties of the QSOs as probed by the rest-frame 2 {\mu}m flux. We find QSOs with a strong excess as well as QSOs with a large deficit at rest-frame 2 {\mu}m relative to a QSO template. Potential solutions to these discrepancies are discussed. Overall, our study demonstrates the high efficiency of the optical/near-infrared selection of red QSOs.

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.

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.

We present a detailed study of the emission from a z = 2.35 galaxy that causes damped Lyman-alpha absorption in the spectrum of the background QSO, SDSS J 2222-0946. We present the results of extensive analyses of the stellar continuum covering the rest frame optical-UV regime based on broad-band HST imaging, and of spectroscopy from VLT/X-Shooter of the strong emission lines: Ly-alpha, [OII], [OIII], [NII], H-alpha and H-beta. We compare the metallicity from the absorption lines in the QSO spectrum with the oxygen abundance inferred from the strong-line methods (R23 and N2). The two emission-line methods yield consistent results: [O/H] = -0.30+/-0.13. Based on the absorption lines in the QSO spectrum a metallicity of -0.49+/-0.05 is inferred at an impact parameter of 6.3 kpc from the centre of the galaxy with a column density of hydrogen of log(N_HI)=20.65+/-0.05. The star formation rates of the galaxy from the UV continuum and H-alpha line can be reconciled assuming an amount of reddening of E(B-V) = 0.06+/-0.01, giving an inferred SFR of 13+/-1 M_sun / yr (assuming a Chabrier IMF). From the HST imaging, the galaxy associated with the absorption is found to be a compact (re=1.12 kpc) object with a disc-like, elongated (axis ratio 0.17) structure indicating that the galaxy is seen close to edge on. Moreover, the absorbing gas is located almost perpendicularly above the disc of the galaxy suggesting that the gas causing the absorption is not co-rotating with the disc. We investigate the stellar and dynamical masses from SED-fitting and emission-line widths, respectively, and find consistent results of 2x10^9 M_sun. We suggest that the galaxy is a young proto-disc with evidence for a galactic outflow of enriched gas. This galaxy hints at how star-forming galaxies may be linked to the elusive population of damped Lyman-alpha absorbers.

The quasar Q0918+1636 (z=3.07) has two intervening high-metallicity Damped Lyman-alpha Absorbers (DLAs) along the line of sight, at redshifts of z=2.412 and 2.583. The z=2.583 DLA is located at a large impact parameter of 16.2 kpc, and despite this large impact parameter it has a very high metallicity (consistent with solar), a substantial fraction of H_2 molecules, and it is dusty as inferred from the reddened spectrum of the background QSO. The z=2.412 DLA has a metallicity of [M/H]=-0.6 (based on ZnII and SiII). In this paper we present new observations of this interesting sightline. HST/WFC3 imaging was obtained in the F606W, F105W and F160W bands. This is complemented by ground-based imaging in the u-, g-bands as well as K_s observations in the near-infrared (NIR). In addition, we present further spectroscopy with the ESO/VLT X-Shooter spectrograph. Based on these observations we obtain the following results: By fitting stellar population synthesis models to the photometric SED we constrain the physical properties of the z=2.583 DLA galaxy, and we infer its morphology by fitting a Sersic model to its surface brightness profile. We find it to be a relatively massive (M_star 10^10 M_sun), strongly star-forming (SFR~30 M_sun / yr, dusty (E_(B-V)=0.4) galaxy with a disk-like morphology. We detect most of the strong emission lines from the z=2.583 DLA [OIII],3727, [OIII],4960, [OIII],5007, Hbeta, and Halpha, albeit at low signal-to-noise (SN) ratio except for the [OIII],5007 line. We also detect [OIII],5007 emission from the galaxy counterpart of the z=2.412 DLA at a small impact parameter (<2 kpc). Overall our findings are consistent with the emerging picture that high-metallicity DLAs are associated with relatively (compared to typical DLAs) luminous and massive galaxy counterparts.

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 report the discovery of a unique gravitational lens system, SDSSJ2222+2745, producing five spectroscopically confirmed images of a z_s=2.82 quasar lensed by a foreground galaxy cluster at z_l=0.49. We also present photometric and spectroscopic evidence for a sixth lensed image of the same quasar. The maximum separation between the quasar images is 15.1". Both the large image separations and the high image multiplicity of the lensed quasar are in themselves exceptionally rare, and observing the combination of these two factors is an exceptionally unlikely occurrence in present datasets. This is only the third known case of a quasar lensed by a cluster, and the only one with six images. The lens system was discovered in the course of the Sloan Giant Arcs Survey, in which we identify candidate lenses in the Sloan Digital Sky Survey and target these for follow up and verification with the 2.56m Nordic Optical Telescope. Multi-band photometry obtained over multiple epochs from September 2011 to September 2012 reveal significant variability at the ~10-30% level in some of the quasar images, indicating that measurements of the relative time delay between quasar images will be feasible. In this lens system we also identify a bright (g = 21.5) giant arc corresponding to a strongly lensed background galaxy at z_s=2.30. We fit parametric models of the lens system, constrained by the redshift and positions of the quasar images and the redshift and position of the giant arc. The predicted time delays between different pairs of quasar images range from ~100 days to ~6 years.

We present the results of a search for red QSOs using a selection based on optical imaging from SDSS and near-infrared imaging from UKIDSS. For a sample of 58 candidates 46 (79%) are confirmed to be QSOs. The QSOs are predominantly dust-reddened except a handul at redshifts z>3.5. The dust is most likely located in the QSO host galaxies. 4 (7%) of the candidates turned out to be late-type stars, and another 4 (7%) are compact galaxies. The remaining 4 objects we could not identify. In terms of their optical spectra the QSOs are similar to the QSOs selected in the FIRST-2MASS red Quasar survey except they are on average fainter, more distant and only two are detected in the FIRST survey. We estimate the amount of extinction using the SDSS QSO template reddened by SMC-like dust. It is possible to get a good match to the observed (restframe ultraviolet) spectra, but for nearly all the reddened QSOs it is not possible to match the near-IR photometry from UKIDSS. The likely reasons are that the SDSS QSO template is too red at optical wavelengths due to contaminating host galaxy light and that the assumed SMC extinction curve is too shallow. Our survey has demonstrated that selection of QSOs based on near-IR photometry is an efficent way to select QSOs, including reddened QSOs, with only small contamination from late-type stars and compact galaxies. This will be useful with ongoing and future wide-field near-IR surveys such as the VISTA and EUCLID surveys. [Abridged]

Recently, the number of detected galaxy counterparts of z > 2 Damped Lyman-alpha Absorbers in QSO spectra has increased substantially so that we today have a sample of 10 detections. M{\o}ller et al. in 2004 made the prediction, based on a hint of a luminosity-metallicity relation for DLAs, that HI size should increase with increasing metallicity. In this paper we investigate the distribution of impact parameter and metallicity that would result from the correlation between galaxy size and metallicity. We compare our observations with simulated data sets given the relation of size and metallicity. The observed sample presented here supports the metallicity-size prediction: The present sample of DLA galaxies is consistent with the model distribution. Our data also show a strong relation between impact parameter and column density of HI. We furthermore compare the observations with several numerical simulations and demonstrate that the observations support a scenario where the relation between size and metallicity is driven by feedback mechanisms controlling the star-formation efficiency and outflow of enriched gas.

We report on our serendipitous discovery that the objects Q 01323-4037 and Q 0132-4037, listed in the V\'eron-Cetty & V\'eron catalog (2006) as two different quasars, are actually a quasar and a star. We briefly discuss the origin of the misidentification, and provide a refined measurement of the quasar redshift.

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