Exploring star formation in high-z galaxies using atomic and molecular emission lines

The conditions under which stars are formed and the reasons for triggering and quenching of starburst events in high-z galaxies, are still not well understood. Studying the interstellar medium (ISM) and the morphology of high-z galaxies are therefore key points in order to understand galaxy evolution. The cosmic star formation rate density peaks between 1<z<3. This period in the history of the Universe is therefore crucial to investigate in order to know more about the star-formation triggering and quenching mechanisms. Phenomena such as major mergers and galactic nuclear activity are believed to be mechanisms dominating the star formation activity at this period of time. It is therefore necessary to study galaxy populations which show signs of major merger events and active galactic nuclei (AGN). This thesis presents three studies of the ISM in high-z galaxies and their morphologies by: Exploring the physical conditions of the ISM in a sample of dusty star-forming galaxies (DSFGs) using the relative observed line strength of ionised carbon ([CII]) and carbon monoxide (CO). We find that the line ratios can best be described by a medium of [CII] and CO emitting gas with a higher [CII] than CO excitation temperature, high CO optical depth tau(CO)>>1, and low to moderate [CII] optical depth tau(CII)<1. Combining millimetre/sub-millimetre and optical data cubes for the high-z radio galaxy (HzRG) MRC0943-242, has revealed a much more complicated morphology than seen in the individual data sets. The millimetre/sub-millimetre observations data have allowed us to spatially separate of the AGN and starburst dominated components, which ~65 kpc apart. The optical data reveal structures of emitting and absorbing gas at multiple wavelengths. A deep high resolution millimetre/sub-millimetre study of the HzRG MRC1138-262, shows emission from water (\water) and an unusually large amount of neutral atomic carbon ([CI]) relative to highly excited CO compared to lensed DSFGs. The detection of water (H2O) emission, which is not associated with the 246GHz continuum emission, suggests excitation by shocks.The uncommon line ratio between [CI]2-1 and CO(7-6) might be due special conditions of the ISM in MRC1138-262 dominated by cosmic rays or differential lensing in other DSFGs thereby not representing the intrinsic ratio. These three studies of the physical conditions of the ISM and morphology of high-z galaxies at z>2, pave the road for future investigations of the star-forming ISM in high-z galaxies, by illustrating the importance of multi-wavelength, fine structure- and molecular line studies

ALMA finds dew drops in the dusty spider’s web

We present 0.̋5 resolution ALMA detections of the observed 246 GHz continuum, [CI] 3P2→3P1 fine structure line ([CI]2–1), CO(7–6), and H2O lines in the z = 2.161 radio galaxy MRC1138-262, the so-called Spiderweb galaxy. We detect strong [CI]2–1 emission both at the position of the radio core, and in a second component ~4 kpc away from it. The 1100 km s-1 broad [CI]2–1 line in this latter component, combined with its H2 mass of 1.6 × 1010 M⊙, implies that this emission must come from a compact region <60 pc, possibly containing a second active galactic nucleus (AGN). The combined H2 mass derived for both objects, using the [CI]2–1 emission, is 3.3 × 1010 M⊙. The total CO(7–6)/[CI]2–1 line flux ratio of 0.2 suggests a low excitation molecular gas reservoir and/or enhanced atomic carbon in cosmic ray dominated regions. We detect spatially-resolved H2O 211−202 emission – for the first time in a high-z unlensed galaxy – near the outer radio lobe to the east, and near the bend of the radio jet to the west of the radio galaxy. No underlying 246 GHz continuum emission is seen at either position. We suggest that the H2O emission is excited in the cooling region behind slow (10–40 km s-1) shocks in dense molecular gas (103−5 cm-3). The extended water emission is likely evidence of the radio jet’s impact on cooling and forming molecules in the post-shocked gas in the halo and inter-cluster gas, similar to what is seen in low-z clusters and other high-z radio galaxies. These observations imply that the passage of the radio jet in the interstellar and inter-cluster medium not only heats gas to high temperatures, as is commonly assumed or found in simulations, but also induces cooling and dissipation, which can lead to substantial amounts of cold dense molecular gas. The formation of molecules and strong dissipation in the halo gas of MRC1138-262 may explain both the extended diffuse molecular gas and the young stars observed around MRC1138-262

3D tomography of the giant Lyα nebulae of z ≈ 3–5 radio-loud AGN

Lyα emission nebulae are ubiquitous around high-redshift galaxies and are tracers of the gaseous environment on scales out to ≳100 pkpc (proper kiloparsec). High-redshift radio galaxies (HzRGs, type-2 radio-loud quasars) host large-scale nebulae observed in the ionised gas differ from those seen in other types of high-redshift quasars. In this work, we exploit MUSE observations of Lyα nebulae around eight HzRGs (2.92 < z < 4.51). All of the HzRGs have large-scale Lyα emission nebulae with seven of them extended over 100 pkpc at the observed surface brightness limit (∼2 − 20 × 10−19 erg s−1 cm−2 arcsec−2). Because the emission line profiles are significantly affected by neutral hydrogen absorbers across the entire nebulae extent, we performed an absorption correction to infer maps of the intrinsic Lyα surface brightness, central velocity, and velocity width, all at the last scattering surface of the observed Lyα photons. We find the following: (i) that the intrinsic surface brightness radial profiles of our sample can be described by an inner exponential profile and a power law in the low luminosity extended part; (ii) our HzRGs have a higher surface brightness and more asymmetric nebulae than both radio-loud and radio-quiet type-1 quasars; (iii) intrinsic nebula kinematics of four HzRGs show evidence of jet-driven outflows but we find no general trends for the whole sample; (iv) a relation between the maximum spatial extent of the Lyα nebula and the projected distance between the active galactic nuclei (AGN) and the centroids of the Lyα nebula; and (v) an alignment between radio jet position angles and the Lyα nebula morphology. All of these findings support a scenario in which the orientation of the AGN has an impact on the observed nebular morphologies and resonant scattering may affect the shape of the surface brightness profiles, nebular kinematics, and relations between the observed Lyα morphologies. Furthermore, we find evidence showing that the outskirts of the ionised gas nebulae may be ‘contaminated’ by Lyα photons from nearby emission halos and that the radio jet affects the morphology and kinematics of the nebulae. Overall, this work provides results that allow us to compare Lyα nebulae around various classes of quasars at and beyond cosmic noon (z ∼ 3)

Cosmic Vine: A z=3.44 Large-Scale Structure Hosting Massive Quiescent Galaxies

We report the discovery of a large-scale structure at z=3.44 revealed by JWST data in the EGS field. This structure, dubbed "Cosmic Vine", consists of 20 galaxies with spectroscopic redshifts at $3.43<z<3.45$ and six galaxy overdensities with consistent photometric redshifts, making up a vine-like structure extending over a ~4x0.2 pMpc^2 area. The two most massive galaxies (M*~10^10.9 Msun) of the Cosmic Vine are found to be quiescent with bulge-dominated morphologies ($B/T>70\%$). Comparisons with simulations suggest that the Cosmic Vine would form a cluster with halo mass >10^14 Msun at z=0, and the two massive galaxies are likely forming the brightest cluster galaxies (BCGs). The results unambiguously reveal that massive quiescent galaxies can form in growing large-scale structures at z>3, thus disfavoring the environmental quenching mechanisms that require a virialized cluster core. Instead, as suggested by the interacting and bulge-dominated morphologies, the two galaxies are likely quenched by merger-triggered starburst or AGN feedback before falling into a cluster core. Moreover, we found that the observed specific star formation rates of massive quiescent galaxies in z>3 dense environments are two orders of magnitude lower than that of the BCGs in the TNG300 simulation. This discrepancy potentially poses a challenge to the models of massive cluster galaxy formation. Future studies comparing a large sample with dedicated cluster simulations are required to solve the problem.Comment: Submitted to A&

The mysterious morphology of MRC0943-242 as revealed by ALMA and MUSE

© 2016 ESO. We present a pilot study of the z = 2.923 radio galaxy MRC0943-242, where we combine information from ALMA and MUSE data cubes for the first time. Even with modest integration times, we disentangle the AGN and starburst dominated components. These data reveal a highly complex morphology as the AGN, starburst, and molecular gas components show up as widely separated sources in dust continuum, optical continuum, and CO line emission observations. CO(1-0) and CO(8-7) line emission suggest that there is a molecular gas reservoir offset from both the dust and the optical continuum that is located ~90 kpc from the AGN. The UV line emission has a complex structure in emission and absorption. The line emission is mostly due to a large scale ionisation cone energised by the AGN, and a Lya emitting bridge of gas between the radio galaxy and a heavily star-forming set of components. Strangely, the ionisation cone has no Lya emission. We find this is due to an optically thick layer of neutral gas with unity covering fraction spread out over a region of at least ~100 kpc from the AGN. Other less thick absorption components are associated with Lya emitting gas within a few tens of kpc from the radio galaxy and are connected by a bridge of emission. We speculate that this linear structure of dust, Lya and CO emission, and the redshifted absorption seen in the circum nuclear region may represent an accretion flow feeding gas into this massive AGN host galaxy

Sub-millimetre galaxies with Webb:Near-infrared counterparts and multi-wavelength morphology

We utilised the unprecedented depth and resolution of recent early-release science (ERS) JWST observations to define the near-infrared counterparts of sub-millimetre galaxies (SMGs). We identified 45 SCUBA-2 SMG positions within the Cosmic Evolution Early Release Science (CEERS) survey JWST/NIRCam fields. Through an analysis of multi-wavelength p-values, NIRCam colours and predicted SCUBA-2 fluxes, we define 43 JWST/NIRCam counterparts to the SCUBA-2 SMGs, finding a 63 per cent agreement with those identified in prior HST studies. Using EaZy-py, we fitted the available HST and JWST observations to quantify the photometric redshifts of the NIRCam-SMGs, establishing a broad range of redshift from z 0.2-5.4 with a median of z 2.29, in agreement with other studies of SMGs. We identified significant variations in the morphology of the NIRCam-SMGs from isolated discs and spheroidal galaxies to irregular interacting systems. We analysed their rest-frame optical and near-infrared morphological properties (e.g. effective radius (Re), Sérsic index (n), concentration (C), asymmetry (A), clumpiness (S), as well as the Gini and M20 parameters), finding, on average, late-type disc-like morphologies with large scatter into the intermediate and merger regions of the non-parametric parameter space. For the non-merging galaxies, we find a median rest-frame optical size and Sérsic index (and 1δ scatter) of Re = 3.10 ± 1.67 kpc and n = 0.96 ± 0.66. Whilst in the rest-frame near-infrared, we establish more compact, higher Sérsic index morphologies (Re = 1.64 ± 0.97, n = 1.85 ± 0.63). We further establish that both the rest-frame optical and near-infrared effective radii correlate negatively (at a 2δ level) with redshift, whilst the Sérsic index remains constant with cosmic time. Our results are consistent with the picture of inside-out galaxy evolution, with more centrally concentrated older stellar populations, and more extended, younger star-forming regions whose stellar emission is heavily attenuated in the central regions.</p

Discovery of methanimine (CH2NH) megamasers toward compact obscured galaxy nuclei

We present the first search for the 5.29 GHz methanimine (CH2NH) 110 - 111 transition toward a sample of galaxy nuclei. We target seven galaxies that host compact obscured nuclei (CONs) with the Karl G. Jansky Very Large Array. These galaxies are characterized by Compton-thick cores. CH2NH emission is detected toward six CONs. The brightness temperatures measured toward Arp 220 indicate maser emission. Isotropic luminosities of the CH2NH transition, from all sources where it is detected, exceed 1 Lpdbl and thus may be considered megamasers. We also detect formaldehyde (H2CO) emission toward three CONs. The isotropic CH2NH luminosities are weakly correlated with the infrared luminosity of the host galaxy and strongly correlated with OH megamaser luminosities from the same galaxies. Non-local thermodynamic equilibrium radiative transfer models suggest that the maser is pumped by the intense millimeter-to-submillimeter Our study suggests that CH2NH megamasers are linked to the nuclear processes within 100 pc of the Compton-thick nucleus within CONs

Massive galaxies on the road to quenching: ALMA observations of powerful high redshift radio galaxies

We present 0.″3 (band 6) and 1.″5 (band 3) ALMA observations of the (sub)millimeter dust continuum emission for 25 radio galaxies at 1 <  z <  5.2. Our survey reaches a rms flux density of ∼50 μJy in band 6 (200–250 GHz) and ∼20 μJy in band 3 (100–130 GHz). This is an order of magnitude deeper than single-dish 850 μm observations, and reaches fluxes where synchrotron and thermal dust emission are expected to be of the same order of magnitude. Combining our sensitive ALMA observations with low-resolution radio data from ATCA, higher resolution VLA data, and infrared photometry from Herschel and Spitzer, we have disentangled the synchrotron and thermal dust emission. We determine the star-formation rates and AGN infrared luminosities using our newly developed Multi-resolution and multi-object/origin spectral energy distribution fitting code (MR-MOOSE). We find that synchrotron emission contributes substantially at λ ∼ 1 mm. Through our sensitive flux limits and accounting for a contribution from synchrotron emission in the mm, we revise downward the median star-formation rate by a factor of seven compared to previous estimates based solely on Herschel and Spitzer data. The hosts of these radio-loud AGN appear predominantly below the main sequence of star-forming galaxies, indicating that the star formation in many of the host galaxies has been quenched. Future growth of the host galaxies without substantial black hole mass growth will be needed to bring these objects on the local relation between the supermassive black holes and their host galaxies. Given the mismatch in the timescales of any star formation that took place in the host galaxies and lifetime of the AGN, we hypothesize that a key role is played by star formation in depleting the gas before the action of the powerful radio jets quickly drives out the remaining gas. This positive feedback loop of efficient star formation rapidly consuming the gas coupled to the action of the radio jets in removing the residual gas is how massive galaxies are rapidly quenched