Source author record

Manuel Aravena

Manuel Aravena appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

43works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

43 published item(s)

preprint2026arXiv

A census of star-formation and gas mass tracers in two lensed $z \sim 4$ dusty star-forming galaxies

We present new and archival Atacama Large Millimeter/submillimeter Array (ALMA) observations of two strongly lensed dusty star-forming galaxies (DSFGs) selected from the South Pole Telescope survey, SPT0418-47 $(z = 4.225)$ and SPT2147-50 $(z = 3.760)$. We study the [C II], CO(7-6), [C I](2-1), and, in SPT0418-47, $p$-H$_2$O emission, which along with the underlying continuum (rest-frame 160 $μ$m and 380 $μ$m) are routinely used as tracers of gas mass and/or star-formation rate (SFR). We perform a pixel-by-pixel analysis of both sources in the image plane to study the resolved Kennicutt-Schmidt relation, finding generally good agreement between the slopes of the SFR versus gas mass surface density using the different tracers. Using lens modeling methods, we find that the dust emission is more compact than the line emission in both sources, with CO(7-6) and [C I](2-1) similar in extent and [C II] the most extended, reminiscent of recent findings of extended [C II] spatial distributions in galaxies at similar cosmic epochs. We develop the [C I](2-1) / CO(7-6) flux density ratio as an observable proxy for gas depletion timescale ($τ_{\rm dep}$), which can be applied to large samples of DSFGs, in lieu of more detailed inferences of this timescale which require analysis of observations at multiple wavelengths. Furthermore, the extended [C II] emission in both sources, compared to the total continuum and line emission, suggests that [C II], used in recent years as a molecular gas mass and SFR tracer in high-$z$ galaxies, may not always be a suitable tracer of these physical quantities.

preprint2026arXiv

REBELS-MOSFIRE: Weak CIII] Emission is Typical Among Extremely UV-bright, Massive Galaxies at $z\sim7$

We present Keck/MOSFIRE H-band spectroscopic measurements covering the [CIII]1907, CIII]1909 doublet for a sample of 8 z~7 spectroscopically-confirmed star-forming galaxies drawn from the Reionization Era Bright Emission Line Survey (REBELS). This REBELS-MOSFIRE sample is notable for its bright median UV luminosity (Muv=-22.5 AB) and large median stellar mass (log(Mstar/Msun)=9.2). Although three sources show tentative evidence of a CIII] detection, we obtain no confident detections for any of the 8 REBELS-MOSFIRE sources. The median [CIII]1907+CIII]1909 3-sigma upper limit in equivalent width (EW) for the REBELS-MOSFIRE sample is 6.5 AA, and a stack of their H-band MOSFIRE spectra yields a non-detection with an associated 3-sigma upper limit of 2.6 AA. These upper limits fall significantly below the CIII] EW measured in a composite spectrum of representative z~7 star-forming galaxies, as well as those measured for notable early star-forming galaxies such as GN-z11, GHZ2, GS-z12, and RXCJ2248-ID. The lack of strong CIII] emission can be understood within the context of the stellar populations of the REBELS galaxies, as well as the ionization conditions and gas-phase metallicity implied by rest-frame optical spectroscopic properties ([OIII]+Hb EWs, and [OIII]5007/[OII]3727 and [NeIII]3869/[OII]3727 line ratios). The REBELS-MOSFIRE sample represents the higher-mass, higher-metallicity, lower-excitation tail of the z~7 galaxy population, whose ionizing properties must be fully characterized to constrain the role of star-forming galaxies during cosmic reionization.

preprint2024arXiv

A benchmark for extreme conditions of the multiphase interstellar medium in the most luminous hot dust-obscured galaxy at z = 4.6

WISE J224607.6-052634.9 (W2246-0526) is a hot dust-obscured galaxy at $z$ = 4.601, and the most luminous obscured quasar known to date. W2246-0526 harbors a heavily obscured supermassive black hole that is most likely accreting above the Eddington limit. We present observations with the Atacama Large Millimeter/submillimeter Array (ALMA) in seven bands, including band 10, of the brightest far-infrared (FIR) fine-structure emission lines of this galaxy: [OI]$_{63μm}$, [OIII]$_{88μm}$, [NII]$_{122μm}$, [OI]$_{145μm}$, [CII]$_{158μm}$, [NII]$_{205μm}$, [CI]$_{370μm}$, and [CI]$_{609μm}$. A comparison of the data to a large grid of Cloudy radiative transfer models reveals that a high hydrogen density ($n_{H}\sim3\times10^3$ cm$^{-3}$) and extinction ($A_{V}\sim300$ mag), together with extreme ionization ($log(U)=-0.5$) and a high X-ray to UV ratio ($α_{ox}\geq-0.8$) are required to reproduce the observed nuclear line ratios. The values of $α_{ox}$ and $U$ are among the largest found in the literature and imply the existence of an X-ray-dominated region (XDR). In fact, this component explains the a priori very surprising non-detection of the [OIII]$_{88μm}$ emission line, which is actually suppressed, instead of boosted, in XDR environments. Interestingly, the best-fitted model implies higher X-ray emission and lower CO content than what is detected observationally, suggesting the presence of a molecular gas component that should be further obscuring the X-ray emission over larger spatial scales than the central region that is being modeled. These results highlight the need for multiline infrared observations to characterize the multiphase gas in high redshift quasars and, in particular, W2246-0526 serves as an extreme benchmark for comparisons of interstellar medium conditions with other quasar populations at cosmic noon and beyond.

preprint2023arXiv

Brightest Cluster Galaxy Formation in the z=4.3 Protocluster SPT2349-56: Discovery of a Radio-Loud AGN

We have observed the z=4.3 protocluster SPT2349-56 with ATCA with the aim of detecting radio-loud active galactic nuclei (AGN) amongst the ~30 submillimeter galaxies identified in the structure. We detect the central complex of SMGs at 2.2\,GHz with a luminosity of L_2.2=(4.42pm0.56)x10^{25} W/Hz. The ASKAP also detects the source at 888 MHz, constraining the radio spectral index to alpha=-1.6pm0.3, consistent with ATCA non-detections at 5.5 and 9GHz, and implying L_1.4(rest)=(2.4pm0.3)x10^{26}W/Hz. This radio luminosity is about 100 times higher than expected from star formation, assuming the usual FIR-radio correlation, which is a clear indication of an AGN driven by a forming brightest cluster galaxy (BCG). None of the SMGs in SPT2349-56 show signs of AGN in any other diagnostics available to us (notably 12CO out to J=16, OH163um, CII/IR, and optical spectra), highlighting the radio continuum as a powerful probe of obscured AGN in high-z protoclusters. No other significant radio detections are found amongst the cluster members, consistent with the FIR-radio correlation. We compare these results to field samples of radio sources and SMGs, along with the 22 SPT-SMG gravitational lenses also observed in the ATCA program, as well as powerful radio galaxies at high redshifts. Our results allow us to better understand the effects of this gas-rich, overdense environment on early supermassive black hole (SMBH) growth and cluster feedback. We estimate that (3.3pm0.7)x10^{38} W of power are injected into the growing ICM by the radio-loud AGN, whose energy over 100Myr is comparable to the binding energy of the gas mass of the central halo. The AGN power is also comparable to the instantaneous energy injection from supernova feedback from the 23 catalogued SMGs in the core region of 120kpc projected radius. The SPT2349-56 radio-loud AGN may be providing strong feedback on a nascent ICM.

preprint2023arXiv

Extended Lyman-$α$ emission towards the SPT2349-56 protocluster at $z=4.3$

Context. Deep spectroscopic surveys with the Atacama Large Millimeter/submillimeter Array (ALMA) have revealed that some of the brightest infrared sources in the sky correspond to concentrations of dusty star-forming galaxies (DSFG) at high redshift. Among these, the SPT2349-56 protocluster system at z = 4.304 is amongst the most extreme examples due to its high source density and integrated star formation rate. Aims. We conducted a deep Lyman-$α$ line emission survey around SPT2349-56 using the Multi-Unit Spectroscopic Explorer (MUSE) at Very Large Telescope (VLT) in order to characterize this uniquely dense environment. Methods. Taking advantage of the deep three-dimensional nature of this survey, we performed a sensitive search for Lyman-$α$ emitters (LAEs) toward the core and northern extension of the protocluster, which correspond to the brightest infrared regions in this field. Using a smoothed narrowband image extracted from the MUSE datacube around the protocluster redshift, we searched for possible extended structures. Results. We identify only three LAEs at z = 4.3 in this field, in concordance with expectations for blank-fields, and an extended Lyman-$α$ structure spatially associated with core of the protocluster. All the previously-identified DSFGs in this field are undetected in Lyman-$α$ emission, consistent with the conspicuous dust obscuration in these systems. We find an extended Lyman-$α$ structure, about $60 \times 60$ kpc$^{2}$ in size, and located 56 kpc west of the protocluster core. Three DSFGs coincide spatially with the location of this structure. We conclude that either the three co-spatial DSFGs or the protocluster core itself are feeding ionizing photons to the Lyman-$α$ structure.

preprint2022arXiv

Chaotic and Clumpy Galaxy Formation in an Extremely Massive Reionization-Era Halo

The SPT0311-58 system at z=6.900 is an extremely massive structure within the reionization epoch, and offers a chance to understand the formation of galaxies in an extreme peak in the primordial density field. We present 70mas Atacama Large Millimeter/submillimeter Array observations of the dust continuum and CII 158um emission in the central pair of galaxies and reach physical resolution ~100-350pc, among the most detailed views of any reionization-era system to date. The observations resolve the source into at least a dozen kiloparsec-size clumps. The global kinematics and high turbulent velocity dispersion within the galaxies present a striking contrast to recent claims of dynamically cold thin-disk kinematics in some dusty galaxies just 800Myr later at z~4. We speculate that both gravitational interactions and fragmentation from massive parent disks have likely played a role in the overall dynamics and formation of clumps in the system. Each clump individually is comparable in mass to other 6<z<8 galaxies identified in rest-UV/optical deep field surveys, but with star formation rates elevated by ~3-5x. Internally, the clumps themselves bear close resemblance to greatly scaled-up versions of virialized cloud-scale structures identified in low-redshift galaxies. Our observations are qualitatively similar to the chaotic and clumpy assembly within massive halos seen in simulations of high-redshift galaxies.

preprint2022arXiv

CONCERTO: a breakthrough in wide field-of-view spectroscopy at millimeter wavelengths

CarbON CII line in post-rEionization and ReionizaTiOn (CONCERTO) is a low-resolution spectrometer with an instantaneous field-of-view of 18.6 arcmin, operating in the 130-310 GHz transparent atmospheric window. It is installed on the 12-meter Atacama Pathfinder Experiment (APEX) telescope at 5100 m above sea level. The Fourier transform spectrometer (FTS) contains two focal planes hosting a total of 4304 kinetic inductance detectors. The FTS interferometric pattern is recorded on the fly while continuously scanning the sky. One of the goals of CONCERTO is to characterize the large-scale structure of the Universe by observing the integrated emission from unresolved galaxies. This methodology is an innovative technique and is called line intensity mapping. In this paper, we describe the CONCERTO instrument, the effect of the vibration of the FTS beamsplitter, and the status of the CONCERTO main survey.

preprint2022arXiv

Dual constraints with ALMA: new [O III] 88 ${\rm μ}$m and dust-continuum observations reveal the ISM conditions of luminous LBGs at $z \sim 7$

We present new [O III] 88 ${\rm μ}$m observations of five bright $z \sim 7$ Lyman-break galaxies spectroscopically confirmed by ALMA through the [C II] 158 ${\rm μ}$m line, unlike recent [O III] detections where Lyman-${\rm α}$ was used. This nearly doubles the sample of Epoch of Reionisation galaxies with robust ($5 σ$) detections of [C II] and [O III]. We perform a multi-wavelength comparison with new deep HST images of the rest-frame UV, whose compact morphology aligns well with [O III] tracing ionised gas. By contrast, we find more spatially extended [C II] emission likely produced in neutral gas, as indicated by a [N II] 205 ${\rm μ}$m non-detection in one source. We find a positive correlation between the equivalent width of the optical [O III] and H${\rm β}$ lines and the [O III]/[C II] ratio, as seen in local metal-poor dwarf galaxies. Cloudy models of a nebula of typical density harbouring a young stellar population with a high ionisation parameter appear to adequately reproduce the far-infrared lines. Surprisingly, however, our models fail to reproduce the strength of [O III] 88 ${\rm μ}$m, unless we assume an ${\rm α}$/Fe enhancement and a near-solar nebular oxygen abundance. On spatially resolved scales, we find [O III]/[C II] shows a tentative anti-correlation with infrared excess, $L_{\rm IR}/L_{\rm UV}$, also seen on global scales in the local Universe. Finally, we introduce the far-infrared spectral energy distribution fitting code MERCURIUS to show that dust-continuum measurements of one source appear to favour a low dust temperature coupled with a high dust mass. This implies a high stellar metallicity yield and may point towards the need of dust production or grain-growth mechanisms beyond supernovae.

preprint2022arXiv

Multi-Phase ISM in the z = 5.7 Hyperluminous Starburst SPT0346-52

SPT0346-52 (z=5.7) is the most intensely star-forming galaxy discovered by the South Pole Telescope, with Sigma_SFR ~ 4200 Msol yr^-1 kpc^-2. In this paper, we expand on previous spatially-resolved studies, using ALMA observations of dust continuum, [NII]205 micron, [CII]158 micron, [OI]146 micron, and undetected [NII]122 micron and [OI]63 micron emission to study the multi-phase interstellar medium (ISM) in SPT0346-52. We use pixelated, visibility-based lens modeling to reconstruct the source-plane emission. We also model the source-plane emission using the photoionization code CLOUDY and find a supersolar metallicity system. We calculate T_dust = 48.3 K and lambda_peak = 80 micron, and see line deficits in all five lines. The ionized gas is less dense than comparable galaxies, with n_e < 32 cm^-3, while ~20% of the [CII]158 emission originates from the ionized phase of the ISM. We also calculate the masses of several phases of the ISM. We find that molecular gas dominates the mass of the ISM in SPT0346-52, with the molecular gas mass ~4x higher than the neutral atomic gas mass and ~100x higher than the ionized gas mass.

preprint2022arXiv

Rapid build-up of the stellar content in the protocluster core SPT2349$-$56 at $z\,{=}\,4.3$

The protocluster SPT2349$-$56 at $z\,{=}\,4.3$ contains one of the most actively star-forming cores known, yet constraints on the total stellar mass of this system are highly uncertain. We have therefore carried out deep optical and infrared observations of this system, probing rest-frame ultraviolet to infrared wavelengths. Using the positions of the spectroscopically-confirmed protocluster members, we identify counterparts and perform detailed source deblending, allowing us to fit spectral energy distributions in order to estimate stellar masses. We show that the galaxies in SPT2349$-$56 have stellar masses proportional to their high star-formation rates, consistent with other protocluster galaxies and field submillimetre galaxies (SMGs) around redshift 4. The galaxies in SPT2349$-$56 have on average lower molecular gas-to-stellar mass fractions and depletion timescales than field SMGs, although with considerable scatter. We construct the stellar-mass function for SPT2349$-$56 and compare it to the stellar-mass function of $z\,{=}\,1$ galaxy clusters, finding consistent shapes between the two. We measure rest-frame galaxy ultraviolet half-light radii from our HST-F160W imaging, finding that on average the galaxies in our sample are similar in size to typical star-forming galaxies at these redshifts. However, the brightest HST-detected galaxy in our sample, found near the luminosity-weighted centre of the protocluster core, remains unresolved at this wavelength. Hydrodynamical simulations predict that the core galaxies will quickly merge into a brightest cluster galaxy, thus our observations provide a direct view of the early formation mechanisms of this class of object.

preprint2022arXiv

Revealing the Nature of a Lyman-$α$ Halo in a Strongly Lensed Interacting System at $z=2.92$

Spatially extended halos of H I Ly$α$ emission are now ubiquitously found around high-redshift star-forming galaxies. But our understanding of the nature and powering mechanisms of these halos is still hampered by the complex radiative transfer effects of the Ly$α$ line and limited angular resolution. In this paper, we present resolved Multi Unit Spectroscopic Explorer (MUSE) observations of SGAS J122651.3+215220, a strongly-lensed pair of $L^{*}$ galaxies at $z=2.92$ embedded in a Ly$α$ halo of $L_{Lyα}=(6.2\pm1.3)\times10^{42}$ erg s$^{-1}$. Globally, the system shows a line profile that is markedly asymmetric and redshifted, but its width and peak shift vary significantly across the halo. By fitting the spatially binned Ly$α$ spectra with a collection of radiative transfer galactic wind models, we infer a mean outflow expansion velocity of $\approx 211$ km s$^{-1}$, with higher values preferentially found on both sides of the system's major axis. The velocity of the outflow is validated with the blueshift of low-ionization metal absorption lines in the spectra of the central galaxies. We also identify a faint ($M_{1500} \approx -16.7$) companion detected in both Ly$α$ and the continuum, whose properties are in agreement with a predicted population of satellite galaxies that contribute to the extended Ly$α$ emission. Finally, we briefly discuss the impact of the interaction between the central galaxies on the properties of the halo and the possibility of in situ fluorescent Ly$α$ production.

preprint2021arXiv

ALMA reveals large overdensity and strong clustering of galaxies in quasar environments at $z\sim4$

We present an Atacama Large Millimeter/submillimeter Array (ALMA) survey of CO(4--3) line emitting galaxies in 17 quasar fields at $z\sim4$ aimed to perform the first systematic search of dusty galaxies in high$-z$ quasar environments. Our blind search of galaxies around the quasars results in 5 CO emitters with S/N$\geq5.6$ within a projected radius of $R\lesssim1.5\,h^{-1}\,$cMpc and a velocity range of $\rmΔ\,v=\pm1000\,$km\,s$^{-1}$ around the quasar. In blank fields, we expect to detect only 0.28 CO emitters within the same volume, implying a total overdensity of $17.6^{+11.9}_{-7.6}$ in our fields, and indicating that quasars trace massive structures in the early universe. We quantify this overdensity by measuring the small-scale clustering of CO emitters around quasars, resulting in a cross-correlation length of $r_{\rm 0,QG}=8.37^{+2.42}_{-2.04}\,h^{-1}\,$cMpc, assuming a fixed slope $γ=1.8$. This contradicts the reported mild overdensities (x1.4) of Lyman alpha emitters (LAEs) in the same fields at scales of $R\lesssim7\,h^{-1}\,$cMpc which is well described by a cross-correlation length 3 times lower than that measured for CO emitters. We discuss some possibilities to explain this discrepancy, including low star formation efficiency, and excess of dust in galaxies around quasars. Finally, we constrain, for the first time, the clustering of CO emitters at $z\sim4$, finding an auto-correlation length of $r_{\rm 0,CO}=3.14\pm1.71\,h^{-1}\,$cMpc (with $γ=1.8$). Our work together with the previous study of LAEs around quasars traces simultaneously the clustering properties of both optical and dusty galaxy populations in quasars fields, stressing the importance of multi-wavelength studies, and highlighting important questions about galaxy properties in high$-z$ dense environments.

preprint2021arXiv

Characterization of Two 2mm-detected Optically-Obscured Dusty Star-Forming Galaxies

The 2mm Mapping Obscuration to Reionization with ALMA (MORA) Survey was designed to detect high redshift ($z\gtrsim4$), massive, dusty star-forming galaxies (DSFGs). Here we present two, likely high redshift sources, identified in the survey whose physical characteristics are consistent with a class of optical/near-infrared (OIR) invisible DSFGs found elsewhere in the literature. We first perform a rigorous analysis of all available photometric data to fit spectral energy distributions and estimate redshifts before deriving physical properties based on our findings. Our results suggest the two galaxies, called MORA-5 and MORA-9, represent two extremes of the "OIR-dark" class of DSFGs. MORA-5 ($z_{\rm phot}=4.3^{+1.5}_{-1.3}$) is a significantly more active starburst with a star-formation rate of 830$^{+340}_{-190}$M$_\odot$yr$^{-1}$ compared to MORA-9 ($z_{\rm phot}=4.3^{+1.3}_{-1.0}$) whose star-formation rate is a modest 200$^{+250}_{-60}$M$_\odot$yr$^{-1}$. Based on the stellar masses (M$_{\star}\approx10^{10-11}$M$_\odot$), space density ($n\sim(5\pm2)\times10^{-6}$Mpc$^{-3}$, which incorporates two other spectroscopically confirmed OIR-dark DSFGs in the MORA sample at $z=4.6$ and $z=5.9$), and gas depletion timescales ($<1$Gyr) of these sources, we find evidence supporting the theory that OIR-dark DSFGs are the progenitors of recently discovered $3<z<4$ massive quiescent galaxies.

preprint2021arXiv

COLDz: Probing Cosmic Star Formation With Radio Free-free Emission

Radio free-free emission is considered to be one of the most reliable tracers of star formation in galaxies. However, as it constitutes the faintest part of the radio spectrum -- being roughly an order of magnitude less luminous than radio synchrotron emission at the GHz frequencies typically targeted in radio surveys -- the usage of free-free emission as a star formation rate tracer has mostly remained limited to the local Universe. Here we perform a multi-frequency radio stacking analysis using deep Karl G. Jansky Very Large Array observations at 1.4, 3, 5, 10 and 34 GHz in the COSMOS and GOODS-North fields to probe free-free emission in typical galaxies at the peak of cosmic star formation. We find that $z \sim 0.5 - 3$ star-forming galaxies exhibit radio emission at rest-frame frequencies of $\sim 65 - 90$ GHz that is $\sim 1.5 - 2\times$ fainter than would be expected from a simple combination of free-free and synchrotron emission, as in the prototypical starburst galaxy M82. We interpret this as a deficit in high-frequency synchrotron emission, while the level of free-free emission is as expected from M82. We additionally provide the first constraints on the cosmic star formation history using free-free emission at $0.5 \lesssim z \lesssim 3$, which are in good agreement with more established tracers at high redshift. In the future, deep multi-frequency radio surveys will be crucial in order to accurately determine the shape of the radio spectrum of faint star-forming galaxies, and to further establish radio free-free emission as a tracer of high-redshift star formation.

preprint2021arXiv

Mapping Obscuration to Reionization with ALMA (MORA): 2mm Efficiently Selects the Highest-Redshift Obscured Galaxies

We present the characteristics of 2mm-selected sources from the largest Atacama Large Millimeter and submillimeter Array (ALMA) blank-field contiguous survey conducted to-date, the Mapping Obscuration to Reionization with ALMA (MORA) survey covering 184arcmin$^2$ at 2mm. Twelve of the thirteen detections above 5$σ$ are attributed to emission from galaxies, eleven of which are dominated by cold dust emission. These sources have a median redshift of $\langle z_{\rm 2mm}\rangle=3.6^{+0.4}_{-0.3}$ primarily based on optical/near-infrared (OIR) photometric redshifts with some spectroscopic redshifts, with 77$\pm$11% of sources at $z>3$ and 38$\pm$12% of sources at $z>4$. This implies that 2mm selection is an efficient method for identifying the highest redshift dusty star-forming galaxies (DSFGs). Lower redshift DSFGs ($z<3$) are far more numerous than those at $z>3$ yet likely to drop out at 2mm. MORA shows that DSFGs with star-formation rates in excess of 300M$_\odot$ yr$^{-1}$ and relative rarity of $\sim$10$^{-5}$ Mpc$^{-3}$ contribute $\sim$30% to the integrated star-formation rate density between $3<z<6$. The volume density of 2mm-selected DSFGs is consistent with predictions from some cosmological simulations and is similar to the volume density of their hypothesized descendants: massive, quiescent galaxies at $z>2$. Analysis of MORA sources' spectral energy distributions hint at steeper empirically-measured dust emissivity indices than typical literature studies, with $\langleβ\rangle=2.2^{+0.5}_{-0.4}$. The MORA survey represents an important step in taking census of obscured star-formation in the Universe's first few billion years, but larger area 2mm surveys are needed to more fully characterize this rare population and push to the detection of the Universe's first dusty galaxies.

preprint2020arXiv

A Comparison of the Stellar, CO and Dust-Continuum Emission from Three, Star-Forming HUDF Galaxies at $z\sim 2$

We compare the extent of the dust, molecular gas and stars in three star-forming galaxies, at $z= 1.4, 1.6$ and $2.7$, selected from the Hubble Ultra Deep Field based on their bright CO and dust-continuum emission as well as their large rest-frame optical sizes. The galaxies have high stellar masses, $\mathrm{M}_*>10^{11}\mathrm{M}_\odot$, and reside on, or slightly below, the main sequence of star-forming galaxies at their respective redshifts. We probe the dust and molecular gas using subarcsecond Atacama Large Millimeter/submillimeter Array observations of the 1.3 mm continuum and CO line emission, respectively, and probe the stellar distribution using \emph{Hubble Space Telescope} observations at 1.6 \textmu m. We find that for all three galaxies the CO emission appears $\gtrsim 30\%$ more compact than the stellar emission. For the $z= 1.4$ and $2.7$ galaxies, the dust emission is also more compact, by $\gtrsim 50\%$, than the stellar emission, whereas for the $z=1.6$ galaxy, the dust and stellar emission have similar spatial extents. This similar spatial extent is consistent with observations of local disk galaxies. However, most high redshift observations show more compact dust emission, likely due to the ubiquity of central starbursts at high redshift and the limited sensitivity of many of these observations. Using the CO emission line, we also investigate the kinematics of the cold interstellar medium in the galaxies, and find that all three have kinematics consistent with a rotation-dominated disk.

preprint2020arXiv

COLDz: A High Space Density of Massive Dusty Starburst Galaxies ~1 Billion Years after the Big Bang

We report the detection of CO($J$=2$\to$1) emission from three massive dusty starburst galaxies at $z$$>$5 through molecular line scans in the NSF's Karl G. Jansky Very Large Array (VLA) CO Luminosity Density at High Redshift (COLDz) survey. Redshifts for two of the sources, HDF 850.1 ($z$=5.183) and AzTEC-3 ($z$=5.298), were previously known. We revise a previous redshift estimate for the third source GN10 ($z$=5.303), which we have independently confirmed through detections of CO $J$=1$\to$0, 5$\to$4, 6$\to$5, and [CII] 158 $μ$m emission with the VLA and the NOrthern Extended Milllimeter Array (NOEMA). We find that two currently independently confirmed CO sources in COLDz are "optically dark", and that three of them are dust-obscured galaxies at $z$$>$5. Given our survey area of $\sim$60 arcmin$^2$, our results appear to imply a $\sim$6-55 times higher space density of such distant dusty systems within the first billion years after the Big Bang than previously thought. At least two of these $z$$>$5 galaxies show star-formation rate surface densities consistent with so-called "maximum" starbursts, but we find significant differences in CO excitation between them. This result may suggest that different fractions of the massive gas reservoirs are located in the dense, star-forming nuclear regions - consistent with the more extended sizes of the [CII] emission compared to the dust continuum and higher [CII]-to-far-infrared luminosity ratios in those galaxies with lower gas excitation. We thus find substantial variations in the conditions for star formation between $z$$>$5 dusty starbursts, which typically have dust temperatures $\sim$57%$\pm$25% warmer than starbursts at $z$=2-3 due to their enhanced star formation activity.

preprint2020arXiv

Megaparsec-scale structure around the proto-cluster core SPT2349$-$56 at $z\,{=}\,4.3$

We present an extensive ALMA spectroscopic follow-up programme of the $z\,{=}\,4.3$ structure SPT2349$-$56, one of the most actively star-forming proto-cluster cores known, to identify additional members using their [C{\sc ii}] 158\,$μ$m and \mbox{CO(4--3)} lines. In addition to robustly detecting the 14 previously published galaxies in this structure, we identify a further 15 associated galaxies at $z\,{=}\,4.3$, resolving 55$\,{\pm}\,$5\,per cent of the 870-$μ$m flux density at 0.5\,arcsec resolution compared to 21\,arcsec single-dish data. These galaxies are distributed into a central core containing 23 galaxies extending out to 300\,kpc in diameter, and a northern extension, offset from the core by 400\,kpc, containing three galaxies. We discovered three additional galaxies in a red {\it Herschel\/}-SPIRE source 1.5\,Mpc from the main structure, suggesting the existence of many other sources at the same redshift as SPT2349$-$56 that are not yet detected in the limited coverage of our data. An analysis of the velocity distribution of the central galaxies indicates that this region may be virialized with a mass of (9$\pm$5)$\,{\times}\,$10$^{12}$\,M$_{\odot}$, while the two offset galaxy groups are about 30 and 60\,per cent less massive and show significant velocity offsets from the central group. We calculate the [C{\sc ii}] and far-infrared number counts, and find evidence for a break in the [C{\sc ii}] luminosity function. We estimate the average SFR density within the region of SPT2349$-$56 containing single-dish emission (a proper diametre of 720\,kpc), assuming spherical symmetry, to be roughly 4$\,{\times}\,10^4$\,M$_{\odot}$\,yr$^{-1}$\,Mpc$^{-3}$; this may be an order of magnitude greater than the most extreme examples seen in simulations.

preprint2020arXiv

The ALMA Spectroscopic Survey in the HUDF: A model to explain observed 1.1 and 0.85 millimeter dust continuum number counts

We present a new semi-empirical model for the dust continuum number counts of galaxies at 1.1 millimeter and 850 \micron. Our approach couples an observationally motivated model for the stellar mass and SFR distribution of galaxies with empirical scaling relations to predict the dust continuum flux density of these galaxies. Without a need to tweak the IMF, the model reproduces the currently available observations of the 1.1 millimeter and 850 \micron number counts, including the observed flattening in the 1.1 millimeter number counts below 0.3 mJy \citep{Gonzalez2019numbercounts} and the number counts in discrete bins of different galaxy properties. Predictions of our work include : (1) the galaxies that dominate the number counts at flux densities below 1 mJy (3 mJy) at 1.1 millimeter (850 $μ$m) have redshifts between $z=1$ and $z=2$, stellar masses of $\sim 5\times10^{10}~\rm{M}_\odot$, and dust masses of $\sim 10^{8}~\rm{M}_\odot$; (2) the flattening in the observed 1.1 millimeter number counts corresponds to the knee of the 1.1 millimeter luminosity function. A similar flattening is predicted for the number counts at 850 $μ$m; (3) the model reproduces the redshift distribution of current 1.1 millimeter detections; (4) to efficiently detect large numbers of galaxies through their dust continuum, future surveys should scan large areas once reaching a 1.1 millimeter flux density of 0.1 mJy rather than integrating to fainter fluxes. Our modeling framework also suggests that the amount of information on galaxy physics that can be extracted from the 1.1 millimeter and 850 $μ$m number counts is almost exhausted.

preprint2020arXiv

The ALMA Spectroscopic Survey in the HUDF: Deep 1.2 mm continuum number counts

We present the results from the 1.2 mm continuum image obtained as part of the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS). The 1.2 mm continuum image has a size of 2.9 (4.2) arcmin$^2$ within a primary beam response of 50% (10%) and a rms value of $9.3\thinspace{\rmμJy\thinspace beam^{-1}}$. We detect 35 sources at high significance (Fidelity $\geq0.5$), 32 of these have well characterized near-infrared HST counterparts. We estimate the 1.2 mm number counts to flux levels of $<30\thinspace{\rmμJy}$ in two different ways: we first use the detected sources to constrain the number counts and find a significant flattening of the counts below $S_ν\sim 0.1$ mJy. In a second approach, we constrain the number counts by using a probability of deflection statistics (P(D)) analysis. For this latter approach, we describe new methods to accurately measure the noise in interferometric imaging (employing jack-knifing in the cube and in the visibility plane). This independent measurement confirms the flattening of the number counts. Our analysis of the differential number counts shows that we are detecting $\sim$93% ($\sim$100% if we include the lower fidelity detections) of the total continuum dust emission associated to galaxies in the HUDF. The ancillary data allows us to study the dependence of the 1.2 mm number counts on redshift ($z=0-4$), galaxy dust mass (${\rm M}_{\rm dust}=10^{7}-10^{9}{\rm M}_{\odot}$), stellar mass (${\rm M}_{*}=10^{9}-10^{12}{\rm M}_{\odot}$), and star-formation rate (${\rm SFR}=1-1000\thinspace{\rm M}_{\odot}\thinspace{\rm yr^{-1}}$). In an accompanying paper we show that the number counts are crucial to constrain galaxy evolution models and the understanding of star-forming galaxies at high redshift.

preprint2020arXiv

The ALMA Spectroscopic Survey in the HUDF: The Cosmic Dust and Gas Mass Densities in Galaxies up to $z\sim3$

Using the deepest 1.2 mm continuum map to date in the Hubble Ultra Deep Field obtained as part of the ALMA Spectroscopic Survey (ASPECS) large program, we measure the cosmic density of dust and implied gas (H$_{2}+$H I) mass in galaxies as a function of look-back time. We do so by stacking the contribution from all $H$-band selected galaxies above a given stellar mass in distinct redshift bins, $ρ_{\rm dust}(M_\ast>M,z)$ and $ρ_{\rm gas}(M_\ast>M,z)$. At all redshifts, $ρ_{\rm dust}(M_\ast>M,z)$ and $ρ_{\rm gas}(M_\ast>M,z)$ grow rapidly as $M$ decreases down to $10^{10}\,M_\odot$, but this growth slows down towards lower stellar masses. This flattening implies that at our stellar mass-completeness limits ($10^8\,M_\odot$ and $10^{8.9}\,M_\odot$ at $z\sim0.4$ and $z\sim3$), both quantities converge towards the total cosmic dust and gas mass densities in galaxies. The cosmic dust and gas mass densities increase at early cosmic time, peak around $z\sim2$, and decrease by a factor $\sim4$ and 7, compared to the density of dust and molecular gas in the local universe, respectively. The contribution of quiescent galaxies -- i.e., with little on-going star-formation -- to the cosmic dust and gas mass densities is minor ($\lesssim10\%$). The redshift evolution of the cosmic gas mass density resembles that of the star-formation rate density, as previously found by CO-based measurements. This confirms that galaxies have relatively constant star-formation efficiencies (within a factor $\sim2$) across cosmic time. Our results also imply that by $z\sim0$, a large fraction ($\sim90\%$) of dust formed in galaxies across cosmic time has been destroyed or ejected to the intergalactic medium.

preprint2020arXiv

The properties of the interstellar medium of galaxies across time as traced by the neutral atomic carbon [CI]

We report ALMA observations of the neutral atomic carbon transitions [CI] and multiple CO lines in a sample of $\sim30$ main sequence galaxies at $z\sim1$, including novel information on [CI](2-1) and CO(7-6) for 7 of such normal objects. We complement our observations with a collection of $>200$ galaxies with coverage of similar transitions, spanning the $z=0-4$ redshift interval and a variety of ambient conditions from local to high-redshift starbursts. We find systematic variations in the [CI]/IR and [CI]/high-$J$ ($J=7$) CO luminosity ratios among the various samples. We interpret these differences as increased dense molecular gas fractions and star formation efficiencies in the strongest high-redshift starbursts with respect to normal main sequence galaxies. We further report constant $L'_{\rm [CI]2-1}$/$L'_{\rm [CI]1-0}$ ratios across the galaxy populations and redshifts, suggesting that gas temperatures $T_{\rm exc}$ traced by [CI] do not strongly vary. We find only a mild correlation with $T_{\rm dust}$ and that, generally, $T_{\rm exc} \lesssim T_{\rm dust}$. We fit the line ratios with classical PDR models, retrieving consistently larger densities and intensities of the UV radiation fields in submm galaxies than in main sequence and local objects. However, these simple models fall short in representing the complexity of a multiphase interstellar medium and should be treated with caution. Finally, we compare our observations with the Santa Cruz semi-analytical model of galaxy evolution, recently extended to simulate submm emission. While we confirm the success in reproducing the CO lines, we find systematically larger [CI] luminosities at fixed IR luminosity than predicted theoretically. This highlights the necessity of improving our understanding of the mechanisms regulating the [CI] emission on galactic scales. We release our data compilation to the community.

preprint2020arXiv

VLA-ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (VLASPECS): Total Cold Gas Masses and CO Line Ratios for z=2-3 "Main Sequence" Galaxies

Using the NSF's Karl G. Jansky Very Large Array (VLA), we report six detections of CO(J=1-0) emission and one upper limit in z=2-3 galaxies originally detected in higher-J CO emission in the Atacama Large submillimeter/Millimeter Array (ALMA) Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS). From the CO(J=1-0) line strengths, we measure total cold molecular gas masses of M_gas = 2.4-11.6 x 10^10 (alpha_CO/3.6) Msun. We also measure a median CO(J=3-2) to CO(J=1-0) line brightness temperature ratio of r_31 = 0.84 +/- 0.26, and a CO(J=7-6) to CO(J=1-0) ratio range of r_71 <0.05 to 0.17. These results suggest that CO(J=3-2) selected galaxies may have a higher CO line excitation on average than CO(J=1-0) selected galaxies, based on the limited, currently available samples from the ASPECS and VLA CO Luminosity Density at High Redshift (COLDz) surveys. This implies that previous estimates of the cosmic density of cold gas in galaxies based on CO(J=3-2) measurements should be revised down by a factor of ~=2 on average based on assumptions regarding CO excitation alone. This correction further improves the agreement between the best currently existing constraints on the cold gas density evolution across cosmic history from line scan surveys, and the implied characteristic gas depletion times.

preprint2019arXiv

Physical Characterization of an Unlensed Dusty Star-Forming Galaxy at $z=5.85$

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$μ$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 $τ_{\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).

preprint2019arXiv

The ALMA Spectroscopic Survey in the HUDF: CO emission lines and 3 mm continuum sources

The ALMA SPECtroscopic Survey in the {\it Hubble} Ultra Deep Field is an ALMA large program that obtained a frequency scan in the 3\,mm band to detect emission lines from the molecular gas in distant galaxies. We here present our search strategy for emission lines and continuum sources in the HUDF. We compare several line search algorithms used in the literature, and critically account for the line-widths of the emission line candidates when assessing significance. We identify sixteen emission lines at high fidelity in our search. Comparing these sources to multi-wavelength data we find that all sources have optical/infrared counterparts. Our search also recovers candidates that have lower significance that can be used statistically to derive, e.g. the CO luminosity function. We apply the same detection algorithm to obtain a sample of six 3 mm continuum sources. All of these are also detected in the 1.2 mm continuum with optical/near-infrared counterparts. We use the continuum sources to compute 3 mm number counts in the sub-mJy regime, and find them to be higher by an order of magnitude than expected for synchrotron-dominated sources. However, the number counts are consistent with those derived at shorter wavelengths (0.85--1.3\,mm) once extrapolating to 3\,mm with a dust emissivity index of $β=1.5$, dust temperature of 35\,K and an average redshift of $z=2.5$. These results represent the best constraints to date on the faint end of the 3 mm number counts.

preprint2019arXiv

The ALMA Spectroscopic Survey in the HUDF: Constraining cumulative CO emission at $1 \lesssim z \lesssim 4$ with power spectrum analysis of ASPECS LP data from 84 to 115 GHz

We present a power spectrum analysis of the ALMA Spectroscopic Survey Large Program (ASPECS LP) data from 84 to 115 GHz. These data predominantly probe small-scale fluctuations ($k=10$-$100$ h Mpc$^{-1}$) in the aggregate CO emission in galaxies at $1 \lesssim z \lesssim 4$. We place an integral constraint on CO luminosity functions (LFs) in this redshift range via a direct measurement of their second moments in the three-dimensional (3D) auto-power spectrum, finding a total CO shot noise power $P_{\textrm{CO,CO}}(k_{\textrm{CO(2-1)}}) \leq 1.9\times10^2$ $μ$K$^2$ (Mpc h$^{-1}$)$^3$. This upper limit ($3σ$) is consistent with the observed ASPECS CO LFs in Decarli et al. 2019, but rules out a large space in the range of $P_{\textrm{CO,CO}}(k_{\textrm{CO(2-1)}})$ inferred from these LFs, which we attribute primarily to large uncertainties in the normalization $Φ_*$ and knee $L_*$ of the Schechter-form CO LFs at $z > 2$. Also, through power spectrum analyses of ASPECS LP data with 415 positions from galaxies with available optical spectroscopic redshifts, we find that contributions to the observed mean CO intensity and shot noise power of MUSE galaxies are largely accounted for by ASPECS blind detections, though there are $\sim20$% contributions to the CO(2-1) mean intensity due to sources previously undetected in the blind line search. Finally, we sum the fluxes from individual blind CO detections to yield a lower limit on the mean CO surface brightness at 99 GHz of $\langle T_{\textrm{CO}} \rangle = 0.55\pm0.02$ $μ$K, which we estimate represents $68$-$80$% of the total CO surface brightness at this frequency.

preprint2019arXiv

The interstellar medium content of galaxies in the ALMA era

The advent of the Atacama Large Millimeter/submillimeter Array (ALMA) has enabled a new era for studies of the formation and assembly of distant galaxies. Cosmological deep field surveys with ALMA and other interferometers have flourished in the last few years covering wide ranges of galaxy properties and redshift, and allowing us to gain critical insights into the physical mechanisms behind the galaxy growth. Here, we present a brief review of recent studies that aim to characterize the interstellar medium properties of galaxies at high redshift ($z>1$), focusing on blank-field ALMA surveys of dust continuum and molecular line emission. In particular, we show recent results from the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS) large program.

preprint2016arXiv

(Sub)millimetre interferometric imaging of a sample of COSMOS/AzTEC submillimetre galaxies IV. Physical properties derived from spectral energy distributions

We characterise the physical nature of a 1.1 mm-selected, flux-limited, and interferometrically followed up sample of SMGs in COSMOS. We used the MAGPHYS code to fit the multiwavelength (UV-radio) SEDs of 16 of the target SMGs. We also constructed the pure radio SEDs of our SMGs using three different radio bands (325 MHz, 1.4 GHz, and 3 GHz). Moreover, since two SMGs in our sample, AzTEC1 and AzTEC3, benefit from previous CO line observations, we studied their properties in more detail. We found that 63% of our target SMGs lie above the galaxy main-sequence by more than a factor of 3, and hence are starbursts. The 3 GHz radio sizes we have previously measured for the target SMGs were compared with the present stellar mass estimates, and we found that the z>3 SMGs are fairly consistent with the mass-size relationship of z~2 compact, quiescent galaxies (cQGs). The median IR-radio correlation parameter is found to be q=2.27, which is lower than measured locally (median q=2.64). AzTEC1 is found to have a sub-Eddington SFR surface density (by a factor of 2.6), while AzTEC3 appears to be an Eddington-limited starburster. The gas reservoir in these two high-z SMGs would be exhausted in only ~86 and 19 Myr at the current SFR, respectively. A comparison of the MAGPHYS-based properties of our SMGs with those of equally bright ALESS SMGs suggests that the two populations share fairly similar physical characteristics, including the q parameter. A hint of negative correlation is found between the 3 GHz size and the level of starburstiness, and hence cosmic-ray electrons in more compact starbursts might be more susceptible to free-free absorption. Some of the derived low and high q values (compared to the local median) could be the result of a specific merger/post-starburst phase of galaxy evolution. Overall, our results support the scenario where z>3 SMGs evolve into today's giant ellipticals.

preprint2016arXiv

ALMA Imaging and Gravitational Lens Models of South Pole Telescope-Selected Dusty, Star-Forming Galaxies at High Redshifts

The South Pole Telescope has discovered one hundred gravitationally lensed, high-redshift, dusty, star-forming galaxies (DSFGs). We present 0.5" resolution 870um Atacama Large Millimeter/submillimeter Array imaging of a sample of 47 DSFGs spanning z=1.9-5.7, and construct gravitational lens models of these sources. Our visibility-based lens modeling incorporates several sources of residual interferometric calibration uncertainty, allowing us to properly account for noise in the observations. At least 70% of the sources are strongly lensed by foreground galaxies (mu_870um > 2), with a median magnification mu_870um = 6.3, extending to mu_870um > 30. We compare the intrinsic size distribution of the strongly lensed sources to a similar number of unlensed DSFGs and find no significant differences in spite of a bias between the magnification and intrinsic source size. This may indicate that the true size distribution of DSFGs is relatively narrow. We use the source sizes to constrain the wavelength at which the dust optical depth is unity and find this wavelength to be correlated with the dust temperature. This correlation leads to discrepancies in dust mass estimates of a factor of 2 compared to estimates using a single value for this wavelength. We investigate the relationship between the [CII] line and the far-infrared luminosity and find that the same correlation between the [CII]L_FIR ratio and Sigma_FIR found for low-redshift star-forming galaxies applies to high-redshift galaxies and extends at least two orders of magnitude higher in Sigma_FIR. This lends further credence to the claim that the compactness of the IR-emitting region is the controlling parameter in establishing the "[CII] deficit."

preprint2016arXiv

ALMA spectroscopic survey in the Hubble Ultra Deep Field: CO luminosity functions and the evolution of the cosmic density of molecular gas

In this paper we use ASPECS, the ALMA Spectroscopic Survey in the {\em Hubble} Ultra Deep Field (UDF) in band 3 and band 6, to place blind constraints on the CO luminosity function and the evolution of the cosmic molecular gas density as a function of redshift up to $z\sim 4.5$. This study is based on galaxies that have been solely selected through their CO emission and not through any other property. In all of the redshift bins the ASPECS measurements reach the predicted `knee' of the CO luminosity function (around $5\times10^{9}$ K km/s pc$^2$). We find clear evidence of an evolution in the CO luminosity function with respect to $z\sim 0$, with more CO luminous galaxies present at $z\sim 2$. The observed galaxies at $z\sim 2$ also appear more gas-rich than predicted by recent semi-analytical models. The comoving cosmic molecular gas density within galaxies as a function of redshift shows a factor 3-10 drop from $z \sim 2$ to $z \sim 0$ (with significant error bars), and possibly a decline at $z>3$. This trend is similar to the observed evolution of the cosmic star formation rate density. The latter therefore appears to be at least partly driven by the increased availability of molecular gas reservoirs at the peak of cosmic star formation ($z\sim2$).

preprint2016arXiv

ALMA Spectroscopic Survey in the Hubble Ultra Deep Field: Molecular gas reservoirs in high-redshift galaxies

We study the molecular gas properties of high-$z$ galaxies observed in the ALMA Spectroscopic Survey (ASPECS) that targets a $\sim1$ arcmin$^2$ region in the Hubble Ultra Deep Field (UDF), a blind survey of CO emission (tracing molecular gas) in the 3mm and 1mm bands. Of a total of 1302 galaxies in the field, 56 have spectroscopic redshifts and correspondingly well-defined physical properties. Among these, 11 have infrared luminosities $L_{\rm{}IR}>10^{11}$ L$_\odot$, i.e. a detection in CO emission was expected. Out these, 7 are detected at various significance in CO, and 4 are undetected in CO emission. In the CO-detected sources, we find CO excitation conditions that are lower than typically found in starburst/SMG/QSO environments. We use the CO luminosities (including limits for non-detections) to derive molecular gas masses. We discuss our findings in context of previous molecular gas observations at high redshift (star-formation law, gas depletion times, gas fractions): The CO-detected galaxies in the UDF tend to reside on the low-$L_{\rm{}IR}$ envelope of the scatter in the $L_{\rm{}IR}-L'_{\rm{}CO}$ relation, but exceptions exist. For the CO-detected sources, we find an average depletion time of $\sim$ 1 Gyr, with significant scatter. The average molecular-to-stellar mass ratio ($M_{\rm{}H2}$/$M_*$) is consistent with earlier measurements of main sequence galaxies at these redshifts, and again shows large variations among sources. In some cases, we also measure dust continuum emission. On average, the dust-based estimates of the molecular gas are a factor $\sim$2-5$\times$ smaller than those based on CO. Accounting for detections as well as non-detections, we find large diversity in the molecular gas properties of the high-redshift galaxies covered by ASPECS.

preprint2016arXiv

ALMA Spectroscopic Survey in the Hubble Ultra Deep Field: Survey Description

We present the rationale for and the observational description of ASPECS: The ALMA SPECtroscopic Survey in the Hubble Ultra-Deep Field (UDF), the cosmological deep field that has the deepest multi-wavelength data available. Our overarching goal is to obtain an unbiased census of molecular gas and dust continuum emission in high-redshift (z$>$0.5) galaxies. The $\sim$1$'$ region covered within the UDF was chosen to overlap with the deepest available imaging from HST. Our ALMA observations consist of full frequency scans in band 3 (84-115 GHz) and band 6 (212-272 GHz) at approximately uniform line sensitivity ($L'_{\rm CO}\sim$2$\times$10$^{9}$ K km/s pc$^2$), and continuum noise levels of 3.8 $μ$Jy beam$^{-1}$ and 12.7 $μ$Jy beam$^{-1}$, respectively. The molecular surveys cover the different rotational transitions of the CO molecule, leading to essentially full redshift coverage. The [CII] emission line is also covered at redshifts $6.0<z<8.0$. We present a customized algorithm to identify line candidates in the molecular line scans, and quantify our ability to recover artificial sources from our data. Based on whether multiple CO lines are detected, and whether optical spectroscopic redshifts as well as optical counterparts exist, we constrain the most likely line identification. We report 10 (11) CO line candidates in the 3mm (1mm) band, and our statistical analysis shows that $<$4 of these (in each band) are likely spurious. Less than 1/3 of the total CO flux in the low-J CO line candidates are from sources that are not associated with an optical/NIR counterpart. We also present continuum maps of both the band 3 and band 6 observations. The data presented here form the basis of a number of dedicated studies that are presented in subsequent papers.

preprint2016arXiv

ALMA Spectroscopic Survey in the Hubble Ultra Deep Field: The Infrared Excess of UV-selected z=2-10 galaxies as a function of UV-continuum Slope and Stellar Mass

We make use of deep 1.2mm-continuum observations (12.7microJy/beam RMS) of a 1 arcmin^2 region in the Hubble Ultra Deep Field to probe dust-enshrouded star formation from 330 Lyman-break galaxies spanning the redshift range z=2-10 (to ~2-3 Msol/yr at 1sigma over the entire range). Given the depth and area of ASPECS, we would expect to tentatively detect 35 galaxies extrapolating the Meurer z~0 IRX-beta relation to z>~2 (assuming T_d~35 K). However, only 6 tentative detections are found at z>~2 in ASPECS, with just three at >3sigma. Subdividing z=2-10 galaxies according to stellar mass, UV luminosity, and UV-continuum slope and stacking the results, we only find a significant detection in the most massive (>10^9.75 Msol) subsample, with an infrared excess (IRX=L_{IR}/L_{UV}) consistent with previous z~2 results. However, the infrared excess we measure from our large selection of sub-L* (<10^9.75 Msol) galaxies is 0.11(-0.42)(+0.32) and 0.14(-0.14)(+0.15) at z=2-3 and z=4-10, respectively, lying below even an SMC IRX-beta relation (95% confidence). These results demonstrate the relevance of stellar mass for predicting the IR luminosity of z>~2 galaxies. We furthermore find that the evolution of the IRX-stellar mass relationship depends on the evolution of the dust temperature. If the dust temperature increases monotonically with redshift (as (1+z)^0.32) such that T_d~44-50 K at z>=4, current results are suggestive of little evolution in this relationship to z~6. We use these results to revisit recent estimates of the z>~3 SFR density. One less obvious implication is in interpreting the high Halpha EWs seen in z~5 galaxies: our results imply that star-forming galaxies produce Lyman-continuum photons at twice the efficiency (per unit UV luminosity) as implied in conventional models. Star-forming galaxies can then reionize the Universe, even if the escape fraction is <10%.

preprint2015arXiv

(Sub)millimetre interferometric imaging of a sample of COSMOS/AzTEC submillimetre galaxies - I. Multiwavelength identifications and redshift distribution

We used the PdBI to map a sample of 15 SMGs in the COSMOS field at the wavelength of 1.3 mm. The target SMGs were originally discovered in the JCMT/AzTEC 1.1 mm continuum survey at S/N=4-4.5. This paper presents, for the first time, interferometric millimetre-wavelength observations of these sources. The angular resolution of our observations, 1.8", allowed us to accurately determine the positions of the target SMGs. Using a detection threshold of S/N>4.5 regardless of multiwavelength counterpart association, and 4<S/N<=4.5 if a multiwavelength counterpart within 1.5" is also present, the total number of detections in our survey is 22. Three of our detected SMGs (AzTEC21, 27, and 28; which corresponds to 20%) are marginally resolved at our angular resolution, and these sources are found to have elongated or clumpy morphologies and/or multiple components. Using optical to NIR photometric redshifts, available spectroscopic redshifts, and redshifts estimated from the radio-to-submm spectral index we infer a median redshift of $\tilde{z}=3.20\pm0.25$ for our sample. To study the overall multiplicity and redshift distribution of flux-limited samples of SMGs we combined these sources with the 15 brightest AzTEC SMGs detected at 1.1 mm, AzTEC1-15, and studied previously. This constitutes a complete, flux- and S/N-limited 1.1-mm selected sample. We find that the median redshift for the 15 brightest AzTEC SMGs ($\tilde{z}=3.05\pm0.44$) is consistent with that for AzTEC16-30. This conforms to recent observational findings that SMGs do not exhibit any significant trend between the redshift and (sub)mm flux density. For the combined AzTEC1-30 sample we derive a median redshift of $\tilde{z}=3.17\pm0.27$, consistent with previous results based on mm-selected samples. We further infer that within the combined AzTEC1-30 sample $\sim25\pm9\%$ of sources separate into multiple components.

preprint2015arXiv

(Sub)millimetre interferometric imaging of a sample of COSMOS/AzTEC submillimetre galaxies - II. The spatial extent of the radio-emitting regions

Radio emission at cm wavelengths from highly star-forming galaxies, such as SMGs, is dominated by synchrotron radiation arising from supernova activity. Using deep, high-resolution ($1σ=2.3$ $μ$Jy beam$^{-1}$; $0.75^{"}$) cm radio-continuum observations taken by the VLA-COSMOS 3 GHz Large Project, we studied the radio-emitting sizes of a flux-limited sample of SMGs in the COSMOS field. Of the 39 SMGs studied here, 3 GHz emission was detected towards 18 of them ($\sim46\pm11\%$) with S/N ratios in the range of ${\rm S/N=4.2-37.4}$. Using 2D elliptical Gaussian fits, we derived a median deconvolved major axis FWHM size of $0.54^{"}\pm 0.11^{"}$ for our 18 SMGs detected at 3 GHz. For the 15 SMGs with known redshift we derived a median linear major axis FWHM of $4.2\pm0.9$ kpc. No clear correlation was found between the radio-emitting size and the 3 GHz or submm flux density, or the redshift of the SMG. However, there is a hint of larger radio sizes at $z\sim2.5-5$ compared to lower redshifts. The sizes we derived are consistent with previous SMG sizes measured at 1.4 GHz and in mid-$J$ CO emission, but significantly larger than those seen in the (sub)mm continuum emission. One possible scenario is that SMGs have i) an extended gas component with a low dust temperature, and which can be traced by low- to mid-$J$ CO line emission and radio continuum emission, and ii) a warmer, compact starburst region giving rise to the high-$J$ line emission of CO, which could dominate the dust continuum size measurements. Because of the rapid cooling of CR electrons in dense starburst galaxies ($\sim10^4-10^5$ yr), the more extended synchrotron radio-emitting size being a result of CR diffusion seems unlikely. Instead, if SMGs are driven by galaxy mergers the radio synchrotron emission might arise from more extended magnetised ISM around the starburst region.

preprint2015arXiv

Globular Cluster Populations: Results Including S$^4$G Late-Type Galaxies

Using 3.6 and 4.5$μ$m images of 73 late-type, edge-on galaxies from the S$^4$G survey, we compare the richness of the globular cluster populations of these galaxies to those of early type galaxies that we measured previously. In general, the galaxies presented here fill in the distribution for galaxies with lower stellar mass, M$_*$, specifically $\log({\rm M}_*/{\rm M}_\odot) < 10$, overlap the results for early-type galaxies of similar masses, and, by doing so, strengthen the case for a dependence of the number of globular clusters per $10^9\ {\rm M}_\odot$ of galaxy stellar mass, T$_{\rm N}$, on M$_*$. For $8.5 < \log ({\rm M}_*/{\rm M}_\odot) < 10.5$ we find the relationship can be satisfactorily described as T$_{\rm N} = ({\rm M}_*/10^{6.7})^{-0.56}$ when M$_*$ is expressed in solar masses. The functional form of the relationship is only weakly constrained and extrapolation outside this range is not advised. Our late-type galaxies, in contrast to our early-types, do not show the tendency for low mass galaxies to split into two T$_{\rm N}$ families. Using these results and a galaxy stellar mass function from the literature, we calculate that in a volume limited, local Universe sample, clusters are most likely to be found around fairly massive galaxies (M$_* \sim 10^{10.8}$ M$_\odot$) and present a fitting function for the volume number density of clusters as a function of parent galaxy stellar mass. We find no correlation between T$_{\rm N}$ and large-scale environment, but do find a tendency for galaxies of fixed M$_*$ to have larger T$_{\rm N}$ if they have converted a larger proportion of their baryons into stars.

preprint2015arXiv

The nature of the [CII] emission in dusty star-forming galaxies from the SPT survey

We present [CII] observations of 20 strongly lensed dusty star forming galaxies at 2.1 < z < 5.7 using APEX and Herschel. The sources were selected on their 1.4 mm flux (S_1.4mm > 20 mJy) from the South Pole Telescope survey, with far-infrared (FIR) luminosities determined from extensive photometric data. The [CII] line is robustly detected in 17 sources, all but one being spectrally resolved. Eleven out of 20 sources observed in [CII] also have low-J CO detections from ATCA. A comparison with mid- and high-J CO lines from ALMA reveals consistent [CII] and CO velocity profiles, suggesting that there is little differential lensing between these species. The [CII], low-J CO and FIR data allow us to constrain the properties of the interstellar medium. We find [CII] to CO(1-0) luminosity ratios in the SPT sample of 5200 +- 1800, with significantly less scatter than in other samples. This line ratio can be best 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. The geometric structure of photodissociation regions allows for such conditions.

preprint2014arXiv

Globular Cluster Populations: First Results from S$^4$G Early-Type Galaxies

Using 3.6$μ$m images of 97 early-type galaxies, we develop and verify methodology to measure globular cluster populations from the S$^4$G survey images. We find that 1) the ratio, T$_{\rm N}$, of the number of clusters, N$_{\rm CL}$, to parent galaxy stellar mass, M$_*$, rises weakly with M$_*$ for early-type galaxies with M$_* > 10^{10}$ M$_\odot$ when we calculate galaxy masses using a universal stellar initial mass function (IMF), but that the dependence of T$_{\rm N}$ on M$_*$ is removed entirely once we correct for the recently uncovered systematic variation of IMF with M$_*$, and 2) for M$_* < 10^{10}$ M$_\odot$ there is no trend between N$_{\rm CL}$ and M$_*$, the scatter in T$_{\rm N}$ is significantly larger (approaching 2 orders of magnitude), and there is evidence to support a previous, independent suggestion of two families of galaxies. The behavior of N$_{\rm CL}$ in the lower mass systems is more difficult to measure because these systems are inherently cluster poor, but our results may add to previous evidence that large variations in cluster formation and destruction efficiencies are to be found among low mass galaxies. The average fraction of stellar mass in clusters is $\sim$ 0.0014 for M$_* > 10^{10}$ M$_\odot$ and can be as large as $\sim 0.02$ for less massive galaxies. These are the first results from the S$^4$G sample of galaxies, and will be enhanced by the sample of early-type galaxies now being added to S$^4$G and complemented by the study of later type galaxies within S$^4$G.

preprint2011arXiv

Unveiling the mask on the ULIRG-to-QSO transition object [H89]1821+643 at z=0.3: a gas-poor/gas-rich galaxy merger and the implications for CO-based dynamical mass estimates

We report the detection of the CO J=1-0 emission line in [H89]1821+643, one of the most optically luminous QSOs in the local Universe and a template ULIRG-to-QSO transition object, located in a rich cool-core cluster at z=0.297. The CO emission is likely to be extended, highly asymmetric with respect to the center of the host elliptical where the QSO resides, and corresponds to a molecular gas mass of ~ 8.0 x 10^{9} M_sun. The dynamical mass enclosed by the CO line emission region could amount up to ~ 1.7x10^{12} M_sun (80% of the total mass of the elliptical host). The bulk of the CO emission is located at ~9 kpc south-east from the nuclei position, and close to a faint optical structure, suggesting that the CO emission could either represent a gas-rich companion galaxy merging with the elliptical host or a tail-like structure reminiscent of a previous interaction. We argue that the first scenario is more likely given the large masses implied for the CO source, which would imply a highly asymmetrical elliptical host. The close alignment between the CO emission major axis and the radio-plume suggests a possible role of the latter excitation of the ambient gas reservoir. The stacking technique was used to search for CO emission and 3mm continuum emission from galaxies in the surrounding cluster, however no detection was found either toward individual galaxies or the stacked ensemble of galaxies, with a 3 sigma limit of < 1.1x10^9 M_sun for the molecular gas.

preprint2010arXiv

Cold molecular gas in massive disk galaxies at z=1.5

We report the detection of the CO J=1-0 emission line in three near-infrared selected star-forming galaxies at z~1.5 with the Very Large Array (VLA) and the Green Bank telescope (GBT). These observations directly trace the bulk of molecular gas in these galaxies. We find H_2 gas masses of 8.3 \pm 1.9 x 10^{10} M_sun, 5.6 \pm 1.4 x 10^{10} M_sun and 1.23 \pm 0.34 x 10^{11} M_sun for BzK-4171, BzK-21000 and BzK-16000, respectively, assuming a conversion alpha_CO=3.6 M_sun (K km s^{-1} pc^{2})^{-1}. We combined our observations with previous CO 2-1 detections of these galaxies to study the properties of their molecular gas. We find brightness temperature ratios between the CO 2-1 and CO 1-0 emission lines of 0.80_{-0.22}^{+0.35}, 1.22_{-0.36}^{+0.61} and 0.41_{-0.13}^{+0.23} for BzK-4171, BzK-21000 and BzK-16000, respectively. At the depth of our observations it is not possible to discern between thermodynamic equilibrium or sub-thermal excitation of the molecular gas at J=2. However, the low temperature ratio found for BzK-16000 suggests sub-thermal excitation of CO already at J=2. For BzK-21000, a Large Velocity Gradient model of its CO emission confirms previous results of the low-excitation of the molecular gas at J=3. From a stacked map of the CO 1-0 images, we measure a CO 2-1 to CO 1-0 brightness temperature ratio of 0.92_{-0.19}^{+0.28}. This suggests that, on average, the gas in these galaxies is thermalized up to J=2, has star-formation efficiencies of ~100 L_sun (K km s^{-1} pc^2)^{-1} and gas consumption timescales of ~0.4 Gyr, unlike SMGs and QSOs at high redshifts.

preprint2010arXiv

The Spitzer Survey of Stellar Structure in Galaxies (S^4G)

The Spitzer Survey of Stellar Structure in Galaxies S^4G is an Exploration Science Legacy Program approved for the Spitzer post-cryogenic mission. It is a volume-, magnitude-, and size-limited (d < 40 Mpc, |b| > 30 degrees, m_(Bcorr) < 15.5, D25>1') survey of 2,331 galaxies using IRAC at 3.6 and 4.5 microns. Each galaxy is observed for 240 s and mapped to > 1.5 x D25. The final mosaicked images have a typical 1 sigma rms noise level of 0.0072 and 0.0093 MJy / sr at 3.6 and 4.5 microns, respectively. Our azimuthally-averaged surface brightness profile typically traces isophotes at mu_3.6 (AB) (1 sigma) ~ 27 mag arcsec^-2, equivalent to a stellar mass surface density of ~ 1 Msun pc^-2. S^4G thus provides an unprecedented data set for the study of the distribution of mass and stellar structures in the local Universe. This paper introduces the survey, the data analysis pipeline and measurements for a first set of galaxies, observed in both the cryogenic and warm mission phase of Spitzer. For every galaxy we tabulate the galaxy diameter, position angle, axial ratio, inclination at mu_3.6 (AB) = 25.5 and 26.5 mag arcsec^-2 (equivalent to ~ mu_B (AB) =27.2 and 28.2 mag arcsec^-2, respectively). These measurements will form the initial S^4G catalog of galaxy properties. We also measure the total magnitude and the azimuthally-averaged radial profiles of ellipticity, position angle, surface brightness and color. Finally, we deconstruct each galaxy using GALFIT into its main constituent stellar components: the bulge/spheroid, disk, bar, and nuclear point source, where necessary. Together these data products will provide a comprehensive and definitive catalog of stellar structures, mass and properties of galaxies in the nearby Universe.

preprint2009arXiv

Environment of MAMBO galaxies in the COSMOS field

Submillimeter galaxies (SMG) represent a dust-obscured high-redshift population undergoing massive star formation activity. Their properties and space density have suggested that they may evolve into spheroidal galaxies residing in galaxy clusters. In this paper, we report the discovery of compact (~10"-20") galaxy overdensities centered at the position of three SMGs detected with the Max-Planck Millimeter Bolometer camera (MAMBO) in the COSMOS field. These associations are statistically significant. The photometric redshifts of galaxies in these structures are consistent with their associated SMGs; all of them are between z=1.4-2.5, implying projected physical sizes of ~170 kpc for the overdensities. Our results suggest that about 30% of the radio-identified bright SMGs in that redshift range form in galaxy density peaks in the crucial epoch when most stars formed.

preprint2009arXiv

Imaging Atomic and Highly Excited Molecular Gas in a z=6.42 Quasar Host Galaxy: Copious Fuel for an Eddington-Limited Starburst at the End of Cosmic Reionization

We have imaged CO(J=7-6) and CI(3P2-3P1) emission in the host galaxy of the z=6.42 quasar SDSS J114816.64+525150.3 (hereafter: J1148+5251) through observations with the Plateau de Bure Interferometer. The region showing CO(J=7-6) emission is spatially resolved, and its size of 5 kpc is in good agreement with earlier CO(J=3-2) observations. In combination with a revised model of the collisional line excitation in this source, this indicates that the highly excited molecular gas traced by the CO J=7-6 line is subthermally excited (showing only 58+/-8% of the CO J=3-2 luminosity), but not more centrally concentrated. We also detect CI(3P2-3P1) emission in the host galaxy of J1148+5251, but the line is too faint to enable a reliable size measurement. From the CI(3P2-3P1) line flux, we derive a total atomic carbon mass of M_CI=1.1x10^7 M_sun, which corresponds to ~5x10^-4 times the total molecular gas mass. We also searched for H2O(J_KaKc=2_12-1_01) emission, and obtained a sensitive line luminosity limit of L'_H2O<4.4x10^9 K kms pc^2, i.e., <15% of the CO(J=3-2) luminosity. The warm, highly excited molecular gas, atomic gas and dust in this quasar host at the end of cosmic reionization maintain an intense starburst that reaches surface densities as high as predicted by (dust opacity) Eddington limited star formation over kiloparsec scales.