Source author record

Fabian Walter

Fabian Walter 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

175works
4topics
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

175 published item(s)

preprint2025arXiv

MIDIS: Quantifying the AGN component of X-ray-detected galaxies

We combine the deepest X-ray survey from the Chandra Deep Field-South (CDF-S) `7-Ms' survey with the deepest mid-infrared (5.6$ μm$) image from the JWST/MIRI Deep Imaging Survey (MIDIS) in the Hubble Ultra-Deep Field (HUDF) to study the infrared counterparts and point-source emission of 31 X-ray sources with a median, intrinsic, rest-frame X-ray luminosity of $\log_{10}(L_{\rm Xc}^{\rm 0.5-7keV})$=42.04$\pm$0.22 erg $\rm s^{-1}$. The sample includes 24 AGN with a redshift range, as set by the X-ray detectability, of $z \simeq 0.5-3$. Through a multi-wavelength morphological decomposition, employing three separate classifications (visual, parametric and non-parametric) we separate (where present) the luminosity of the point-like AGN component from the remainder of the host-galaxy emission. The unprecedented mid-infrared sensitivity and imaging resolution of MIRI allows, in many cases, the direct characterisation of point-like (i.e. unresolved) components in the galaxies' emission. We establish a broad agreement between the three morphological classifications. At least 70% of the X-ray sources, including some classified as galaxies, show unresolved emission in the MIRI images, with the unresolved-to-total flux fraction at rest-frame 2$μm$ ranging from $\sim$0.2 to $\sim$0.9. At high X-ray luminosities ($\log_{10}(L_{\rm Xc}$)>43 erg $\rm s^{-1}$) we derive a consistent rest-frame near-infrared 2$ μm$ point-source luminosity to that derived for local AGN, whilst at lower X-ray luminosity we identify an excess in the 2$ μm$ emission compared to pre-JWST studies. We speculate this offset may be driven by a combination of Compton-thick AGN components and nuclear starburst, merger driven activity. Our observations highlight the complex nature of X-ray sources in the distant Universe and demonstrate the power of JWST/MIRI in quantifying their nuclear infrared emission. (Abridged)

preprint2024arXiv

ALMA 300 pc resolution imaging of a z=6.79 quasar: no evidence for supermassive black hole influence on the [C II] kinematics

We present Atacama Large Millimeter/submillimeter Array (ALMA) [C II] 158 $μ\rm{m}$ and dust continuum observations of the $z=6.79$ quasar J0109--3047 at a resolution of $0."045$ ($\sim$300 pc). The dust and [C II] emission are enclosed within a $\sim 500\, \rm{pc}$ radius, with the central beam ($r<144\ \rm{pc}$) accounting for $\sim$25\% (8\%) of the total continuum ([C II]) emission. The far--infrared luminosity density increases radially from $\sim$5 $\times 10^{11} L_\odot\ \rm{kpc}^{-2}$ to a central value of $\sim$70 $\times 10^{11} L_\odot\ \rm{kpc}^{-2}$ (SFRD $\sim$50-700 $M_\odot\ \rm{yr}^{-1}\ \rm{kpc}^{-2}$). The [C II] kinematics are dispersion-dominated with a constant velocity dispersion of $137 \pm 6 \,\rm{km\ s}^{-1}$. The constant dispersion implies that the underlying mass distribution is not centrally peaked, consistent with the expectations of a flat gas mass profile. The lack of an upturn in velocity dispersion within the central beam is inconsistent with a black hole mass greater than $M_{\rm{BH}}<6.5\times 10^{8}\ M_\odot\ (2σ$ level), unless highly fine-tuned changes in the ISM properties conspire to produce a decrease of the gas mass in the central beam comparable to the black hole mass. Our observations therefore imply either that a) the black hole is less massive than previously measured or b) the central peak of the far-infrared and [C II] emission are not tracing the location of the black hole, as suggested by the tentative offset between the near-infrared position of the quasar and the ALMA continuum emission.

preprint2022arXiv

(Nearly) Model-Independent Constraints on the Neutral Hydrogen Fraction in the Intergalactic Medium at $z\sim 5-7$ Using Dark Pixel Fractions in Ly$α$ and Ly$β$ Forests

Cosmic reionization was the last major phase transition of hydrogen from neutral to highly ionized in the intergalactic medium (IGM). Current observations show that the IGM is significantly neutral at $z>7$, and largely ionized by $z\sim5.5$. However, most methods to measure the IGM neutral fraction are highly model-dependent, and are limited to when the volume-averaged neutral fraction of the IGM is either relatively low ($\bar{x}_{\rm HI} \lesssim 10^{-3}$) or close to unity ($\bar{x}_{\rm HI}\sim 1$). In particular, the neutral fraction evolution of the IGM at the critical redshift range of $z=6-7$ is poorly constrained. We present new constraints on $\bar{x}_{\rm HI}$ at $z\sim5.1-6.8$, by analyzing deep optical spectra of $53$ quasars at $5.73<z<7.09$. We derive model-independent upper limits on the neutral hydrogen fraction based on the fraction of "dark" pixels identified in the Lyman $α$ (Ly$α$) and Lyman $β$ (Ly$β$) forests, without any assumptions on the IGM model or the intrinsic shape of the quasar continuum. They are the first model-independent constraints on the IGM neutral hydrogen fraction at $z\sim6.2-6.8$ using quasar absorption measurements. Our results give upper limits of $\bar{x}_{\rm HI}(z=6.3) < 0.79\pm0.04$ (1$σ$), $\bar{x}_{\rm HI} (z=6.5) < 0.87\pm0.03$ (1$σ$), and $\bar{x}_{\rm HI} (z=6.7) < 0.94^{+0.06}_{-0.09}$ (1$σ$). The dark pixel fractions at $z>6.1$ are consistent with the redshift evolution of the neutral fraction of the IGM derived from the Planck 2018.

preprint2022arXiv

A Massive, Dusty HI-Absorption-Selected Galaxy at $z \approx 2.46$ Identified in a CO Emission Survey

We report a NOrthern Extended Millimeter Array (NOEMA) and Atacama Large Millimeter/submillimeter Array (ALMA) search for redshifted CO emission from the galaxies associated with seven high-metallicity ([M/H] $\geq -1.03$) damped Lyman-$α$ absorbers (DLAs) at $z\approx1.64-2.51$. Our observations yielded one new detection of CO(3-2) emission from a galaxy at $z=2.4604$ using NOEMA, associated with the $z=2.4628$ DLA towards QSO B0201+365. Including previous searches, our search results in detection rates of CO emission of $\approx56^{+38}_{-24}$ % and $\approx11^{+26}_{-9}$ %, respectively, in the fields of DLAs with ${\rm [M/H]}>-0.3$ and ${\rm [M/H]}<-0.3$. Further, the HI-selected galaxies associated with five DLAs with [M/H] $>-0.3$ all have high molecular gas masses, $\gtrsim5\times10^{10}\ {\rm M}_\odot$. This indicates that the highest-metallicity DLAs at $z\approx2$ are associated with the most massive galaxies. The newly-identified $z\approx2.4604$ HI-selected galaxy, DLA0201+365g, has an impact parameter of $\approx7$ kpc to the QSO sightline, and an implied molecular gas mass of $(5.04\pm0.78) \times10^{10}\times(α_{\rm CO}/4.36)\times(r_{31}/0.55)\ {\rm M}_\odot$. Archival Hubble Space Telescope Wide Field and Planetary Camera 2 imaging covering the rest-frame near-ultraviolet (NUV) and far-ultraviolet (FUV) emission from this galaxy yield non-detections of rest-frame NUV and FUV emission, and a $5σ$ upper limit of 2.3 M$_\odot$ yr$^{-1}$ on the unobscured star formation rate (SFR). The low NUV-based SFR estimate, despite the very high molecular gas mass, indicates that DLA0201+365g either is a very dusty galaxy, or has a molecular gas depletion time that is around two orders of magnitude larger than that of star-forming galaxies at similar redshifts.

preprint2022arXiv

After The Fall: Resolving the Molecular Gas in Post-Starburst Galaxies

Post-starburst (PSB), or 'E+A', galaxies represent a rapid transitional phase between major, gas-rich mergers and gas-poor, quiescent early-type galaxies. Surprisingly, many PSBs have been shown to host a significant interstellar medium (ISM), despite theoretical predictions that the majority of star-forming gas should be expelled in AGN- or starburst-driven outflows. To-date, the resolved properties of this surviving ISM have remained unknown. We present high resolution ALMA continuum and CO(2$-$1) observations in six gas- and dust-rich PSBs, revealing for the first time the spatial and kinematic structure of their ISM on sub-kpc scales. We find extremely compact molecular reservoirs, with dust and gas surface densities rivaling those found in (ultra-)luminous infrared galaxies. We observe spatial and kinematic disturbances in all sources, with some also displaying disk-like kinematics. Estimates of the internal turbulent pressure in the gas exceed those of normal star-forming disks by at least 2 orders of magnitude, and rival the turbulent gas found in local interacting galaxies, such as the Antennae. Though the source of this high turbulent pressure remains uncertain, we suggest that the high incidence of tidal disruption events (TDEs) in PSBs could play a role. The star formation in these PSBs' turbulent central molecular reservoirs is suppressed, forming stars only 10% as efficiently as starburst galaxies with similar gas surface densities. "The fall" of star formation in these galaxies was not precipitated by complete gas expulsion or redistribution. Rather, this high-resolution view of PSBs' ISM indicates that star formation in their remaining compact gas reservoirs is suppressed by significant turbulent heating.

preprint2022arXiv

ALMA 200 pc imaging of a z~7 quasar reveals a compact, disk-like host galaxy

We present 0".035 resolution (~200 pc) imaging of the 158 um [CII] line and the underlying dust continuum of the z=6.9 quasar J234833.34-305410.0. The 18 h ALMA observations reveal extremely compact emission (diameter ~1 kpc) that is consistent with a simple, almost face-on, rotation-supported disk with a significant velocity dispersion of ~160 km/s. The gas mass in just the central 200 pc is ~4x10^9 M_sun, about a factor two higher than that of the central supermassive black hole. Consequently we do not resolve the black hole's sphere of influence, and find no kinematic signature of the central supermassive black hole. Kinematic modeling of the [CII] line shows that the dynamical mass at large radii is consistent with the gas mass, leaving little room for a significant mass contribution by stars and/or dark matter. The Toomre-Q parameter is less than unity throughout the disk, and thus is conducive to star formation, consistent with the high infrared luminosity of the system. The dust in the central region is optically thick, at a temperature >132 K. Using standard scaling relations of dust heating by star formation, this implies an unprecedented high star formation rate density of >10^4 M_sun / yr / kpc^2. Such a high number can still be explained with the Eddington limit for star formation under certain assumptions, but could also imply that the central supermassive black hole contributes to the heating of the dust in the central 110 pc.

preprint2022arXiv

Chemical Abundance of z~6 Quasar Broad-Line Regions in the XQR-30 Sample

The elemental abundances in the broad-line regions of high-redshift quasars trace the chemical evolution in the nuclear regions of massive galaxies in the early universe. In this work, we study metallicity-sensitive broad emission-line flux ratios in rest-frame UV spectra of 25 high-redshift (5.8 < z < 7.5) quasars observed with the VLT/X-shooter and Gemini/GNIRS instruments, ranging over $\log(M_{\rm{BH}}/M_{\odot})= 8.4-9.8$ in black hole mass and $\log(L_{\rm{bol}}/\rm{erg\, s}^{-1})= 46.7-47.7$ in bolometric luminosity. We fit individual spectra and composites generated by binning across quasar properties: bolometric luminosity, black hole mass, and blueshift of the \civ\, line, finding no redshift evolution in the emission-line ratios by comparing our high-redshift quasars to lower-redshift (2.0 < z < 5.0) results presented in the literature. Using Cloudy-based locally optimally-emitting cloud photoionisation model relations between metallicity and emission-line flux ratios, we find the observable properties of the broad emission lines to be consistent with emission from gas clouds with metallicity that are at least 2-4 times solar. Our high-redshift measurements also confirm that the blueshift of the CIV emission line is correlated with its equivalent width, which influences line ratios normalised against CIV. When accounting for the CIV blueshift, we find that the rest-frame UV emission-line flux ratios do not correlate appreciably with the black hole mass or bolometric luminosity.

preprint2022arXiv

Co-evolution of massive black holes and their host galaxies at high redshift: discrepancies from six cosmological simulations and the key role of JWST

The James Webb Space Telescope will have the power to characterize high-redshift quasars at z>6 with an unprecedented depth and spatial resolution. While the brightest quasars at such redshift (i.e., with bolometric luminosity L_bol> 10^46 erg/s) provide us with key information on the most extreme objects in the Universe, measuring the black hole (BH) mass and Eddington ratios of fainter quasars with L_bol= 10^45-10^46 erg/s opens a path to understand the build-up of more normal BHs at z>6. In this paper, we show that the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA large-scale cosmological simulations do not agree on whether BHs at z>4 are overmassive or undermassive at fixed galaxy stellar mass with respect to the M_BH-M_star scaling relation at z=0 (BH mass offsets). Our conclusions are unchanged when using the local scaling relation produced by each simulation or empirical relations. We find that the BH mass offsets of the simulated faint quasar population at z>4, unlike those of bright quasars, represent the BH mass offsets of the entire BH population, for all the simulations. Thus, a population of faint quasars with L_bol= 10^45-10^46 erg/s observed by JWST can provide key constraints on the assembly of BHs at high redshift. Moreover, this will help constraining the high-redshift regime of cosmological simulations, including BH seeding, early growth, and co-evolution with the host galaxies. Our results also motivate the need for simulations of larger cosmological volumes down to z=6, with the same diversity of sub-grid physics, in order to gain statistics on the most extreme objects at high redshift.

preprint2022arXiv

Constraining galaxy overdensities around three z~6.5 quasars with ALMA and MUSE

We quantify galaxy overdensities around three high-redshift quasars with known [CII] 158um companions: PJ231-20 (z=6.59), PJ308-21 (z=6.24) and J0305-3150 (z=6.61). Recent SCUBA2 imaging revealed the presence of 17 submillimeter galaxies (SMG) with sky separations $0.7'< θ< 2.4'$ from these three quasars. We present ALMA Band 6 follow-up observations of these SCUBA2-selected SMGs to confirm their nature and redshift. We also search for continuum-undetected [CII] emitters in the ALMA pointings and make use of archival MUSE observations to search for Lyman-$α$ Emitters (LAE) associated with the quasars. While most of the SCUBA2-selected sources are detected with ALMA in the continuum, no [CII] line emission could be detected, indicating that they are not at the quasar redshifts. Based on the serendipitous detection of CO 7-6 and [CI] emission lines, we find that four SMGs in the field of PJ231-20 are at z~2.4, which is coincident with the redshift of a Mg II absorber in the quasar rest-frame UV spectrum. We report the discovery of 2 LAEs within <0.6 cMpc of PJ231-20 at the same redshift, indicating a LAE overdensity around this quasar. Taken together, these observations provide new constraints on the large-scale excess of Lyman-$α$- and [CII]-emitting galaxies around z>6 quasars and suggest that only wide-field observations, such as MUSE, ALMA or JWST mosaics, can reveal a comprehensive picture of large-scale structure around quasars in the first billion years of the Universe.

preprint2022arXiv

Exploring the radio spectral energy distribution of the ultraluminous radio-quiet quasar SDSS J0100+2802 at redshift 6.3

We report deep Karl G. Jansky Very Large Array (VLA) observations of the optically ultraluminous and radio-quiet quasar SDSS J010013.02 + 280225.8 (hereafter J0100+2802) at redshift $z=$6.3. We detected the radio continuum emission at 1.5 GHz, 6 GHz, and 10 GHz. This leads to a radio power-law spectral index of $α= -0.52\pm0.18$ ($S \propto ν^α$). The radio source is unresolved in all VLA bands with an upper limit to the size of $0.2^{\prime \prime}$ (i.e., $\sim$ 1.1 kpc) at 10 GHz. We find variability in the flux density (increase by $\sim 33\%$) and the spectral index (steepened) between observations in 2016 and 2017. We also find that the VLA 1.5 GHz flux density observed in the same year is 1.5 times that detected with the Very Long Baseline Array (VLBA) in 2016 at the same frequency. This difference suggests that half of the radio emission from J0100+2802 comes from a compact core within 40 pc, and the rest comes from the surrounding few kpc area which is diffuse and resolved out in the VLBA observations. The diffuse emission is four times brighter than that would be expected if driven by star formation. We conclude that the central active galactic nucleus is the dominant power engine of the radio emission in J0100+2802.

preprint2022arXiv

Hydrogen reionisation ends by $z=5.3$: Lyman-$α$ optical depth measured by the XQR-30 sample

The presence of excess scatter in the Ly-$α$ forest at $z\sim 5.5$, together with the existence of sporadic extended opaque Gunn-Peterson troughs, has started to provide robust evidence for a late end of hydrogen reionisation. However, low data quality and systematic uncertainties complicate the use of Ly-$α$ transmission as a precision probe of reionisation's end stages. In this paper, we assemble a sample of 67 quasar sightlines at $z>5.5$ with high signal-to-noise ratios of $>10$ per $\leq 15$ km s$^{-1}$ spectral pixel, relying largely on the new XQR-30 quasar sample. XQR-30 is a large program on VLT/X-Shooter which obtained deep (SNR $>20$ per pixel) spectra of 30 quasars at $z>5.7$. We carefully account for systematics in continuum reconstruction, instrumentation, and contamination by damped Ly-$α$ systems. We present improved measurements of the mean Ly-$α$ transmission over $4.9<z<6.1$. Using all known systematics in a forward modelling analysis, we find excellent agreement between the observed Ly-$α$ transmission distributions and the homogeneous-UVB simulations Sherwood and Nyx up to $z\leq5.2$ ($<1 σ$), and mild tension ($\sim 2.5 σ$) at $z=5.3$. Homogeneous UVB models are ruled out by excess Ly-$α$ transmission scatter at $z\geq5.4$ with high confidence ($>3.5 σ$). Our results indicate that reionisation-related fluctuations, whether in the UVB, residual neutral hydrogen fraction, and/or IGM temperature, persist in the intergalactic medium until at least $z=5.3$ ($t=1.1$ Gyr after the Big Bang). This is further evidence for a late end to reionisation.

preprint2022arXiv

Kiloparsec-scale imaging of the CO(1-0)-traced cold molecular gas reservoir in a z~3.4 submillimeter galaxy

We present a high-resolution study of the cold molecular gas as traced by CO(1-0) in the unlensed z$\sim$3.4 submillimeter galaxy SMM J13120+4242, using multi-configuration observations with the Karl G. Jansky Very Large Array (JVLA). The gas reservoir, imaged on 0.39" ($\sim$3 kpc) scales, is resolved into two components separated by $\sim$11 kpc with a total extent of 16 $\sim$3 kpc. Despite the large spatial extent of the reservoir, the observations show a CO(1-0) FWHM linewidth of only 267 $\pm$ 64 km s$^{-1}$. We derive a revised line luminosity of L'$_\mathrm{CO(1-0)}$ = (10 $\pm$ 3) $\times$ 10$^{10}$ K km s$^{-1}$ pc$^2$ and a molecular gas mass of M$_\mathrm{gas}$ = (13 $\pm$ 3) $\times$ 10$^{10}$ ($α_\mathrm{CO}$/1) M$_{\odot}$. Despite the presence of a velocity gradient (consistent with previous resolved CO(6-5) imaging), the CO(1-0) imaging shows evidence for significant turbulent motions which are preventing the gas from fully settling into a disk. The system likely represents a merger in an advanced stage. Although the dynamical mass is highly uncertain, we use it to place an upper limit on the CO-to-H$_2$ mass conversion factor $α_\mathrm{CO}$ of 1.4. We revisit the SED fitting, finding that this galaxy lies on the very massive end of the main sequence at z = 3.4. Based on the low gas fraction, short gas depletion time and evidence for a central AGN, we propose that SMM J13120 is in a rapid transitional phase between a merger-driven starburst and an unobscured quasar. The case of SMM J13120 highlights the how mergers may drive important physical changes in galaxies without pushing them off the main sequence.

preprint2022arXiv

LADUMA: Discovery of a luminous OH megamaser at $z > 0.5$

In the local Universe, OH megamasers (OHMs) are detected almost exclusively in infrared-luminous galaxies, with a prevalence that increases with IR luminosity, suggesting that they trace gas-rich galaxy mergers. Given the proximity of the rest frequencies of OH and the hyperfine transition of neutral atomic hydrogen (HI), radio surveys to probe the cosmic evolution of HI in galaxies also offer exciting prospects for exploiting OHMs to probe the cosmic history of gas-rich mergers. Using observations for the Looking At the Distant Universe with the MeerKAT Array (LADUMA) deep HI survey, we report the first untargeted detection of an OHM at $z > 0.5$, LADUMA J033046.20$-$275518.1 (nicknamed "Nkalakatha"). The host system, WISEA J033046.26$-$275518.3, is an infrared-luminous radio galaxy whose optical redshift $z \approx 0.52$ confirms the MeerKAT emission line detection as OH at a redshift $z_{\rm OH} = 0.5225 \pm 0.0001$ rather than HI at lower redshift. The detected spectral line has 18.4$σ$ peak significance, a width of $459 \pm 59\,{\rm km\,s^{-1}}$, and an integrated luminosity of $(6.31 \pm 0.18\,{\rm [statistical]}\,\pm 0.31\,{\rm [systematic]}) \times 10^3\,L_\odot$, placing it among the most luminous OHMs known. The galaxy's far-infrared luminosity $L_{\rm FIR} = (1.576 \pm 0.013) \times 10^{12}\,L_\odot$ marks it as an ultra-luminous infrared galaxy; its ratio of OH and infrared luminosities is similar to those for lower-redshift OHMs. A comparison between optical and OH redshifts offers a slight indication of an OH outflow. This detection represents the first step towards a systematic exploitation of OHMs as a tracer of galaxy growth at high redshifts.

preprint2022arXiv

Long Dark Gaps in the Ly$β$ Forest at $z<6$: Evidence of Ultra Late Reionization from XQR-30 Spectra

We present a new investigation of the intergalactic medium (IGM) near reionization using dark gaps in the Lyman-$β$ (Ly$β$) forest. With its lower optical depth, Ly$β$ offers a potentially more sensitive probe to any remaining neutral gas compared to commonly used Ly$α$ line. We identify dark gaps in the Ly$β$ forest using spectra of 42 QSOs at $z_{\rm em}>5.5$, including new data from the XQR-30 VLT Large Programme. Approximately $40\%$ of these QSO spectra exhibit dark gaps longer than $10h^{-1}{\rm Mpc}$ at $z\simeq5.8$. By comparing the results to predictions from simulations, we find that the data are broadly consistent both with models where fluctuations in the Ly$α$ forest are caused solely by ionizing ultraviolet background (UVB) fluctuations and with models that include large neutral hydrogen patches at $z<6$ due to a late end to reionization. Of particular interest is a very long ($L=28h^{-1}{\rm Mpc}$) and dark ($τ_{\rm eff} \gtrsim 6$) gap persisting down to $z\simeq 5.5$ in the Ly$β$ forest of the $z_{\rm}=5.85$ QSO PSO J025$-$11. This gap may support late reionization models with a volume-weighted average neutral hydrogen fraction of $ \langle x_{\rm HI}\rangle \gtrsim 5\%$ by $z=5.6$. Finally, we infer constraints on $\langle x_{\rm HI}\rangle$ over $5.5 \lesssim z \lesssim 6.0$ based on the observed Ly$β$ dark gap length distribution and a conservative relationship between gap length and neutral fraction derived from simulations. We find $\langle x_{\rm HI}\rangle \leq 0.05$, 0.17, and 0.29 at $z\simeq 5.55$, 5.75, and 5.95, respectively. These constraints are consistent with models where reionization ends significantly later than $z = 6$.

preprint2022arXiv

Microwave Background Temperature at a Redshift of 6.34 from H2O Absorption

Distortions of the observed cosmic microwave background imprinted by the Sunyaev-Zel'dovich effect toward massive galaxy clusters due to inverse Compton scattering of microwave photons by high-energy electrons provide a direct measurement of the microwave background temperature at redshifts from 0 to 1. Some additional background temperature estimates exist at redshifts from 1.8 to 3.3 based on molecular and atomic line excitation temperatures in quasar absorption line systems, but are model dependent. To date, no deviations from the expected (1+z) scaling behavior of the microwave background temperature have been seen, but the measurements have not extended deeply into the matter-dominated era of the universe at redshifts z>3.3. Here we report the detection of sub-millimeter line absorption from the water molecule against the cosmic microwave background at z=6.34 in a massive starburst galaxy, corresponding to a lookback time of 12.8 Gyr. Radiative pumping of the upper level of the ground-state ortho-H2O(110-101) line due to starburst activity in the dusty galaxy HFLS3 results in a cooling to below the redshifted microwave background temperature, after the transition is initially excited by the microwave background. The strength of this effect implies a microwave background temperature of 16.4-30.2 K (1-sigma range) at z=6.34, which is consistent with a background temperature increase with redshift as expected from the standard CDM cosmology.

preprint2022arXiv

Molecular gas in z~6 quasar host galaxies

We investigate the molecular gas content of z~6 quasar host galaxies using the IRAM / Northern Extended Millimeter Array. We target the 3mm dust continuum, and the line emission from CO(6-5), CO(7-6), [CI]2-1 in 10 infra-red-luminous quasars that have been previously studied in their 1mm dust continuum and [CII] line emission. We detect CO(7-6) at various degrees of significance in all the targeted sources, thus doubling the number of such detections in z~6 quasars. The 3mm to 1mm flux density ratios are consistent with a modified black body spectrum with a dust temperature $T_{dust}$~47 K and an optical depth $τ_ν$=0.2 at the [CII] frequency. Our study provides us with four independent ways to estimate the molecular gas mass, $M_{H2}$, in the targeted quasars. This allows us to set constraints on various parameters used in the derivation of molecular gas mass estimates, such as the mass per luminosity ratios $α_{CO}$ and $α_{[CII]}$, the gas-to-dust ratio $δ_{g/d}$, and the carbon abundance [C]/H2. Leveraging either on the dust, CO, [CI], or [CII] emission yields mass estimates of the entire sample in the range $M_{H2}$~$10^{10}$ to $10^{11}$ M$_{\odot}$. We compare the observed luminosities of dust, [CII], [CI], and CO(7-6) with predictions from photo-dissociation and X-ray dominated regions. We find that the former provide better model fits to our data, assuming that the bulk of the emission arises from dense ($n_H>10^4$ cm$^{-3}$) clouds with a column density $N_{H}$~$10^{23}$ cm$^{-2}$, exposed to a radiation field with intensity $G_0$~$10^3$ (in Habing units). Our analysis reiterates the presence of massive reservoirs of molecular gas fueling star formation and nuclear accretion in $z$~6 quasar host galaxies. It also highlights the power of combined 3mm and 1mm observations for quantitative studies of the dense gas content in massive galaxies at cosmic dawn.

preprint2022arXiv

Physical constraints on the extended interstellar medium of the z=6.42 quasar J1148+5251: [CII] 158um, [NII] 205um and [OI] 146um observations

We report new Northern Extended Millimeter Array (NOEMA) observations of the [CII], [NII] and [OI] atomic fine structure lines and dust continuum emission of J1148+5251, a z=6.42 quasar, that probe the physical properties of its interstellar medium (ISM). The radially-averaged [CII] and dust continuum emission have similar extensions (up to $θ= 2.51^{+0.46}_{-0.25}\ \rm{arcsec}$, corresponding to $r= 9.8^{+3.3}_{-2.1}\ \rm{kpc}$ accounting for beam-convolution), confirming that J1148+5251 is the quasar with the largest [CII]-emitting has reservoir known at these epochs.Moreover, if the [CII] emission is examined only along its NE-SW axis, a significant excess ($>5.8σ$) of [CII] emission (with respect to the dust) is detected. The new wide--bandwidth observations enable us to accurately constrain the continuum emission, and do not statistically require the presence of broad [CII] line wings that were reported in previous studies. We also report the first detection of the [OI] and (tentatively) [NII] emission lines in J1148+5251. Using Fine Structure Lines (FSL) ratios of the [CII], [NII], [OI] and previously measured [CI] emission lines, we show that J1148+5251 has similar ISM conditions compared to lower--redshift (ultra)-luminous infrared galaxies. CLOUDY modelling of the FSL ratios exclude X--ray dominated regions (XDR) and favours photodissociation regions (PDR) as the origin of the FSL emission. We find that a high radiation field ($10^{3.5-4.5}\,G_0$), high gas density ($n \simeq 10^{3.5-4.5}\, \rm{cm}^{-3}$) and HI column density of $10^{23} \,\rm{cm^{-2}}$ reproduce the observed FSL ratios well.

preprint2022arXiv

Spatially resolved molecular interstellar medium in a $z=6.6$ quasar host galaxy

We present high spatial resolution (0.4", 2.2kpc) observations of the CO(6-5), CO(7-6) and [CI] lines and dust continuum emission from the interstellar medium in the host galaxy of the quasar J0305$-$3150 at $z=6.6$. These, together with archival [CII] data at comparable spatial resolution, enable studies of the spatial distribution and kinematics between the ISM in different phases. When comparing the radial profiles of CO, [CII] and the dust continuum, we find that the CO and dust continuum exhibit similar spatial distributions, both of which are less extended than the [CII], indicating that the CO and dust continuum are tracing the same gas component, while the [CII] is tracing a more extended one. In addition, we derive the radial profiles of the [CII]/CO, [CII]/far-infrared (FIR), CO/FIR, and dust continuum $S_{98.7 \rm GHz}/S_{258.1 \rm GHz}$ ratios. We find a decreasing $S_{98.7 \rm GHz}/S_{258.1 \rm GHz}$ ratio with radius, possibly indicating a decrease of dust optical depth with increasing radius. We also detect some of the ISM lines and continuum emission in the companion galaxies previously discovered in the field around J0305$-$3150. Through comparing the line-to-line and line-to-FIR ratios, we find no significant differences between the quasar and its companion galaxies.

preprint2022arXiv

The [CII] and FIR properties of z>6 radio-loud quasars

There are only five radio-loud quasars currently known within 1 Gyr from the Big Bang ($z>6$) and the properties of their host galaxies have not been explored in detail. We present a NOrthern Extended Millimeter Array (NOEMA) survey of [CII] (158 $μ$m) and underlying continuum emission of four $z>6$ radio-loud quasars, revealing their diverse properties. J0309+2717 ($z=6.10$) has a bright [CII] line and underlying continuum, implying a starburst with a star-formation rate SFR=$340-1200$ $M_\odot$ yr$^{-1}$. J1429+5447 ($z=6.18$) has a SFR=$520-870$ $M_{\odot}$yr$^{-1}$ and its [CII] profile is consistent with two Gaussians, which could be interpreted as a galaxy merger. J1427+3312 ($z=6.12$) has a moderate SFR=$30-90$ $M_\odot$ yr$^{-1}$. Notably, this is a broad absorption line quasar and we searched for the presence of high-velocity outflows in the host galaxy. Although the NOEMA data reveal a tentative broad component of the [CII] line as wide as $\sim$1400~km~s$^{-1}$, the sensitivity of our current data is not sufficient to confirm it. Finally, P172+18 ($z=6.82$) is undetected in both [CII] and continuum, implying a SFR$<22-40$ $M_{\odot}$yr$^{-1}$. The broad range of SFRs is similar to what is observed in radio-quiet quasars at similar redshifts. If radio jets do not significantly contribute to both [CII] and IR luminosities, this suggest no feedback from the jet on the star formation in the host galaxy.

preprint2022arXiv

The Decoupled Kinematics of high-z QSO Host Galaxies and their Lya halos

We present a comparison of the interstellar medium traced by [CII] (ALMA), and ionised halo gas traced by Lya (MUSE), in and around QSO host galaxies at z~6. To date, 18 QSOs at this redshift have been studied with both MUSE and high-resolution ALMA imaging; of these, 8 objects display a Lya halo. Using datacubes matched in velocity resolution, we compare and contrast the spatial and kinematic information of the Lya halos and the host galaxies' [CII] (and dust-continuum) emission. We find that the Lya halos extend typically 3-30 times beyond the interstellar medium of the host galaxies. The majority of the Lya halos do not show ordered motion in their velocity fields, whereas most of the [CII] velocity fields do. In those cases where a velocity gradient can be measured in Lya, the kinematics do not align with those derived from the [CII] emission. This implies that the Lya emission is not tracing the outskirts of a large rotating disk that is a simple extension of the central galaxy seen in [CII] emission. It rather suggests that the kinematics of the halo gas are decoupled from those of the central galaxy. Given the scattering nature of Lya, these results need to be confirmed with JWST IFU observations that can constrain the halo kinematics further using the non-resonant Ha line.

preprint2022arXiv

The Morpho-Kinematic Architecture of Super Star Clusters in the Center of NGC253

The center of the nearby galaxy NGC\,253 hosts a population of more than a dozen super star clusters (SSCs) which are still in the process of forming. The majority of the star formation of the burst is concentrated in these SSCs, and the starburst is powering a multiphase outflow from the galaxy. In this work, we measure the 350~GHz dust continuum emission towards the center of NGC\,253 at 47~milliarcsecond (0.8 pc) resolution using data from the Atacama Large Millimeter/submillimeter Array (ALMA). We report the detection of 350~GHz (dust) continuum emission in the outflow for the first time, associated with the prominent South-West streamer. In this feature, the dust emission has a width of $\approx$~8~pc, is located at the outer edge of the CO emission, and corresponds to a molecular gas mass of $\sim~(8-17)\times10^6$~M$_\odot$. In the starburst nucleus, we measure the resolved radial profiles, sizes, and molecular gas masses of the SSCs. Compared to previous work at somewhat lower spatial resolution, the SSCs here break apart into smaller substructures with radii $0.4-0.7$~pc. In projection, the SSCs, dust, and dense molecular gas appear to be arranged as a thin, almost linear, structure roughly 155~pc in length. The morphology and kinematics of this structure can be well explained as gas following $x_2$ orbits at the center of a barred potential. We constrain the morpho-kinematic arrangement of the SSCs themselves, finding that an elliptical, angular momentum-conserving ring is a good description of the both morphology and kinematics of the SSCs.

preprint2022arXiv

Using Artificial Intelligence and real galaxy images to constrain parameters in galaxy formation simulations

Cosmological galaxy formation simulations are still limited by their spatial/mass resolution and cannot model from first principles some of the processes, like star formation, that are key in driving galaxy evolution. As a consequence they still rely on a set of 'effective parameters' that try to capture the scales and the physical processes that cannot be directly resolved in the simulation. In this study we show that it is possible to use Machine Learning techniques applied to real and simulated images of galaxies to discriminate between different values of these parameters by making use of the full information content of an astronomical image instead of collapsing it into a limited set of values like size, or stellar/ gas masses. In this work we apply our method to the NIHAO simulations and the THINGS and VLA-ANGST observations of HI maps in nearby galaxies to test the ability of different values of the star formation density threshold $n$ to reproduce observed HI maps. We show that observations indicate the need for a high value of $n \gtrsim 80$ ,cm$^{-3}$ (although the exact numerical value is model-dependent), which has important consequences for the dark matter distribution in galaxies. Our study shows that with innovative methods it is possible to take full advantage of the information content of galaxy images and compare simulations and observations in an interpretable, non-parametric and quantitative manner.

preprint2021arXiv

A Luminous Quasar at Redshift 7.642

Distant quasars are unique tracers to study the formation of the earliest supermassive black holes (SMBHs) and the history of cosmic reionization. Despite extensive efforts, only two quasars have been found at $z\ge7.5$, due to a combination of their low spatial density and the high contamination rate in quasar selection. We report the discovery of a luminous quasar at $z=7.642$, J0313$-$1806, the most distant quasar yet known. This quasar has a bolometric luminosity of $3.6\times10^{13} L_\odot$. Deep spectroscopic observations reveal a SMBH with a mass of $(1.6\pm0.4) \times10^9M_\odot$ in this quasar. The existence of such a massive SMBH just $\sim$670 million years after the Big Bang challenges significantly theoretical models of SMBH growth. In addition, the quasar spectrum exhibits strong broad absorption line (BAL) features in CIV and SiIV, with a maximum velocity close to 20% of the speed of light. The relativistic BAL features, combined with a strongly blueshifted CIV emission line, indicate that there is a strong active galactic nucleus (AGN) driven outflow in this system. ALMA observations detect the dust continuum and [CII] emission from the quasar host galaxy, yielding an accurate redshift of $7.6423 \pm 0.0013$ and suggesting that the quasar is hosted by an intensely star-forming galaxy, with a star formation rate of $\rm\sim 200 ~M_\odot ~yr^{-1}$ and a dust mass of $\sim7\times10^7~M_\odot$. Followup observations of this reionization-era BAL quasar will provide a powerful probe of the effects of AGN feedback on the growth of the earliest massive 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

Deep XMM-Newton Observations of an X-ray Weak, Broad Absorption Line Quasar at $z=6.5$

We report X-ray observations of the most distant known gravitationally lensed quasar, J0439+1634 at $z=6.52$, which is also a broad absorption line (BAL) quasar, using the XMM-Newton Observatory. With a 130 ks exposure, the quasar is significantly detected as a point source at the optical position with a total of 358$^{+19}_{-19}$ net counts using the EPIC instrument. By fitting a power-law plus Galactic absorption model to the observed spectra, we obtain a spectral slope of $Γ=1.45^{+0.10}_{-0.09}$. The derived optical-to-X-ray spectral slope $α_{\rm{ox}}$ is $-2.07^{+0.01}_{-0.01}$, suggesting that the X-ray emission of J0439+1634 is weaker by a factor of 18 than the expectation based on its 2500 Angstrom luminosity and the average $α_{\rm{ox}}$ vs. luminosity relationship. This is the first time that an X-ray weak BAL quasar at $z>6$ has been observed spectroscopically. Its X-ray weakness is consistent with the properties of BAL quasars at lower redshift. By fitting a model including an intrinsic absorption component, we obtain intrinsic column densities of $N_{\rm{H}}=2.8^{+0.7}_{-0.6}\times10^{23}\,\rm{cm}^{-2}$ and $N_{\rm{H}}= 4.3^{+1.8}_{-1.5}\times10^{23}\,\rm{cm}^{-2}$, assuming a fixed $Γ$ of 1.9 and a free $Γ$, respectively. The intrinsic rest-frame 2--10 keV luminosity is derived as $(9.4-15.1)\times10^{43}\,\rm{erg\,s}^{-1}$, after correcting for lensing magnification ($μ=51.3$). The absorbed power-law model fitting indicates that J0439+1634 is the highest redshift obscured quasar with a direct measurement of the absorbing column density. The intrinsic high column density absorption can reduce the X-ray luminosity by a factor of $3-7$, which also indicates that this quasar could be a candidate of intrinsically X-ray weak quasar.

preprint2021arXiv

HI mapping of the Leo Triplet: Morphologies and kinematics of tails and bridges

A fully-sampled and hitherto highest resolution and sensitivity observation of neutral hydrogen (HI) in the Leo Triplet (NGC 3628, M 65/NGC 3623, and M 66/NGC 3627) reveals six HI structures beyond the three galaxies. We present detailed results of the morphologies and kinematics of these structures, which can be used for future simulations. In particular, we detect a two-arm structure in the plume of NGC 3628 for the first time, which can be explained by a tidal interaction model. The optical counterpart of the plume is mainly associated with the southern arm. The connecting part (base) of the plume (directed eastwards) with NGC 3628 is located at the blueshifted (western) side of NGC 3628. Two bases appear to be associated with the two arms of the plume. A clump with reversed velocity gradient (relative to the velocity gradient of M 66) and a newly detected tail, i.e. M 66SE, is found in the southeast of M 66. We suspect that M 66SE represents gas from NGC 3628 which was captured by M 66 in the recent interaction between the two galaxies. Meanwhile gas is falling toward M 66, resulting in features already previously observed in the southeastern part of M 66, e.g. large line widths and double peaks. An upside-down `Y'-shaped HI gas component (M 65S) is detected in the south of M 65 which suggests that M 65 may also have been involved in the interaction. We strongly encourage modern hydrodynamical simulations of this interacting group of galaxies to reveal the origin of the gaseous debris surrounding all three galaxies.

preprint2021arXiv

Molecular Gas Properties and CO-to-H2 Conversion Factors in the Central Kiloparsec of NGC 3351

The CO-to-H$_2$ conversion factor ($α_\rm{CO}$) is critical to studying molecular gas and star formation in galaxies. The value of $α_\rm{CO}$ has been found to vary within and between galaxies, but the specific environmental conditions that cause these variations are not fully understood. Previous observations on $\sim$kpc scales revealed low values of $α_\rm{CO}$ in the centers of some barred spiral galaxies, including NGC 3351. We present new ALMA Band 3, 6, and 7 observations of $^{12}$CO, $^{13}$CO, and C$^{18}$O lines on 100 pc scales in the inner $\sim$2 kpc of NGC 3351. Using multi-line radiative transfer modeling and a Bayesian likelihood analysis, we infer the H$_2$ density, kinetic temperature, CO column density per line width, and CO isotopologue abundances on a pixel-by-pixel basis. Our modeling implies the existence of a dominant gas component with a density of $2{-}3\times10^3$ $\rm{cm^{-3}}$ in the central ${\sim}$1 kpc and a high temperature of 30$-$60 K near the nucleus and near the contact points that connect to the bar-driven inflows. Assuming a CO/H$_2$ abundance of $3\times10^{-4}$, our analysis yields $α_\rm{CO}{\sim}0.5{-}2.0$ $\rm{M_\odot\,(K~km~s^{-1}~pc^2)^{-1}}$ with a decreasing trend with galactocentric radius in the central $\sim$1 kpc. The inflows show a substantially lower $α_\rm{CO} < 0.1$ $\rm{M_\odot\,(K~km~s^{-1}~pc^2)^{-1}}$, likely due to lower optical depths caused by turbulence or shear in the inflows. Over the whole region, this gives an intensity-weighted $α_\rm{CO}$ of ${\sim}1.5$ $\rm{M_\odot\,(K~km~s^{-1}~pc^2)^{-1}}$, which is similar to previous dust modeling based results at kpc scales. This suggests that low $α_\rm{CO}$ on kpc scales in the centers of some barred galaxies may be due to the contribution of low optical depth CO emission in bar-driven inflows.

preprint2021arXiv

Probing Early Super-massive Black Hole Growth and Quasar Evolution with Near-infrared Spectroscopy of 37 Reionization-era Quasars at 6.3 < z <= 7.64

We report the results of near-infrared spectroscopic observations of 37 quasars in the redshift range $6.3< z\le7.64$, including 32 quasars at $z>6.5$, forming the largest quasar near-infrared spectral sample at this redshift. The spectra, taken with Keck, Gemini, VLT, and Magellan, allow investigations of central black hole mass and quasar rest-frame ultraviolet spectral properties. The black hole masses derived from the MgII emission lines are in the range $(0.3-3.6)\times10^{9}\,M_{\odot}$, which requires massive seed black holes with masses $\gtrsim10^{3-4}\,M_{\odot}$, assuming Eddington accretion since $z=30$. The Eddington ratio distribution peaks at $λ_{\rm Edd}\sim0.8$ and has a mean of 1.08, suggesting high accretion rates for these quasars. The CIV - MgII emission line velocity differences in our sample show an increase of CIV blueshift towards higher redshift, but the evolutionary trend observed from this sample is weaker than the previous results from smaller samples at similar redshift. The FeII/MgII flux ratios derived for these quasars up to $z=7.6$, compared with previous measurements at different redshifts, do not show any evidence of strong redshift evolution, suggesting metal-enriched environments in these quasars. Using this quasar sample, we create a quasar composite spectrum for $z>6.5$ quasars and find no significant redshift evolution of quasar broad emission lines and continuum slope, except for a blueshift of the CIV line. Our sample yields a strong broad absorption line quasar fraction of $\sim$24%, higher than the fractions in lower redshift quasar samples, although this could be affected by small sample statistics and selection effects.

preprint2021arXiv

The radio spectral turnover of radio-loud quasars at $z>5$

We present Karl G. Jansky Very Large Array (VLA) S- (2--4 GHz), C- (4--8 GHz), and X-band (8--12 GHz) continuum observations toward seven radio-loud quasars at $z>5$. This sample has previously been found to exhibit spectral peaks at observed-frame frequencies above $\sim$1 GHz. We also present upgraded Giant Metrewave Radio Telescope (uGMRT) band-2 (200 MHz), band-3 (400 MHz), and band-4 (650 MHz) radio continuum observations toward eight radio-loud quasars at $z>5$, selected from our previous GMRT survey, in order to sample their low-frequency synchrotron emission. Combined with archival radio continuum observations, all ten targets show evidence for spectral turnover. The turnover frequencies are $\sim$1--50 GHz in the rest frame, making these targets gigahertz-peaked-spectrum (GPS) or high-frequency-peaker (HFP) candidates. For the nine well-constrained targets with observations on both sides of the spectral turnover, we fit the entire radio spectrum with absorption models associated with synchrotron self-absorption and free-free absorption (FFA). Our results show that FFA in an external inhomogeneous medium can accurately describe the observed spectra for all nine targets, which may indicate an FFA origin for the radio spectral turnover in our sample. As for the complex spectrum of J114657.79+403708.6 at $z=5.00$ with two spectral peaks, it may be caused by multiple components (i.e., core-jet) and FFA by the high-density medium in the nuclear region. However, we cannot rule out the spectral turnover origin of variability. Based on our radio spectral modeling, we calculate the radio loudness $R_{2500\rm\, Å}$ for our sample, which ranges from 12$^{+1}_{-1}$ to 674$^{+61}_{-51}$.

preprint2020arXiv

A Closer Look at Two of the Most Luminous Quasars in the Universe

Ultra-luminous quasars ($M_{1450} \leq -29$) provide us with a rare view into the nature of the most massive and most rapidly accreting supermassive black holes (SMBHs). Following the discovery of two of these extreme sources, J0341${+}$1720 ($M_{1450}=-29.56$, $z=3.71$) and J2125${-}$1719 ($M_{1450}=-29.39$, $z=3.90$), in the Extremely Luminous Quasar Survey (ELQS) and its extension to the Pan-STARRS\,1 footprint (PS-ELQS), we herein present an analysis of their rest-frame UV to optical spectroscopy. Both quasars harbor very massive SMBHs with $M_{\rm{BH}}=6.73_{-0.83}^{+0.75}\times10^{9}\,M_{\odot}$ and $M_{\rm{BH}}=5.45_{-0.55}^{+0.60}\times10^{9}\,M_{\odot}$, respectively, showing evidence of accretion above the Eddington limit ($L_{\rm{bol}}/L_{\rm{Edd}}=2.74_{-0.27}^{+0.39}$ and $L_{\rm{bol}}/L_{\rm{Edd}}=3.01_{-0.30}^{+0.34}$). NOEMA 3 millimeter observations of J0341${+}$1720 reveal a highly star-forming ($\rm{SFR}\approx1500\,M_{\odot}\,\rm{yr}^{-1}$), ultra-luminous infrared galaxy ($L_{\rm{TIR}}\approx1.0\times10^{13}\,L_{\odot}$) host, which, based on an estimate of its dynamical mass, is only ${\sim}30$ times more massive than the SMBH it harbors at its center. As examples of luminous super-Eddington accretion, these two quasars provide support for theories, which explain the existence of billion solar mass SMBHs ${\sim}700$ million years after the Big Bang by moderate super-Eddington growth from standard SMBH seeds.

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

A metal-poor damped Ly-alpha system at redshift 6.4

We identify a strong Ly-alpha damping wing profile in the spectrum of the quasar P183+05 at z=6.4386. Given the detection of several narrow metal absorption lines at z=6.40392, the most likely explanation for the absorption profile is that it is due to a damped Ly-alpha system. However, in order to match the data a contribution of an intergalactic medium 5-38% neutral or additional weaker absorbers near the quasar is also required. The absorption system presented here is the most distant damped Ly-alpha system currently known. We estimate an HI column density ($10^{20.68\pm0.25}\,$cm$^{-2}$), metallicity ([O/H]$=-2.92\pm 0.32$), and relative chemical abundances of a system consistent with a low-mass galaxy during the first Gyr of the universe. This object is among the most metal-poor damped Ly-alpha systems known and, even though it is observed only ~850 Myr after the big bang, its relative abundances do not show signatures of chemical enrichment by Population III stars.

preprint2020arXiv

A search for dust and molecular gas in enormous Ly$α$ nebulae at $z\approx 2$

Enormous Ly$α$ nebulae, extending over 300-500\,kpc around quasars, represent the pinnacle of galaxy and cluster formation. Here we present IRAM Plateau de Bure Interferometer observations of the enormous Ly$α$ nebulae `Slug' ($z$=$2.282$) and `Jackpot' ($z$=$2.041$). Our data reveal bright, synchrotron emission associated with the two radio-loud AGN embedded in the targeted nebulae, as well as molecular gas, as traced via the CO(3-2) line, in three galaxies (two sources in the Slug, and one in the Jackpot). All of the CO emission is associated with galaxies detected in their rest-frame UV stellar emission. The total mass in molecular gas of these three galaxies [$\sim (3-5)\times10^{10}$ M$_\odot$] is comparable with the total ionized gas mass responsible for the diffuse nebular emission. Our observations place limits on the molecular gas emission in the nebulae: The molecular gas surface density is $Σ_{\rm H2}<12-25$ M$_\odot$ pc$^{-2}$ for the Slug nebula and $Σ_{\rm H2}<34-68$ M$_\odot$ pc$^{-2}$ for the Jackpot nebula. These are consistent with the expected molecular gas surface densities, as predicted via photoionization models of the rest-frame UV line emission in the nebulae, and via Ly$α$ absorption in the Jackpot nebula. Compared to other radio--loud quasars at $z>1$, and high-redshift radio--loud galaxies, we do not find any strong trends relating the molecular gas reservoirs, the radio power, and the Ly$α$ luminosities of these systems. The significant step in sensitivity required to achieve a detection of the molecular gas from the nebulae, if present, will require a substantial time investment with JVLA, NOEMA, or ALMA.

preprint2020arXiv

A Significantly Neutral Intergalactic Medium Around the Luminous z=7 Quasar J0252-0503

Luminous $z\ge7$ quasars provide direct probes of the evolution of supermassive black holes (SMBHs) and the intergalactic medium (IGM) during the epoch of reionization (EoR). The Ly$α$ damping wing absorption imprinted by neutral hydrogen in the IGM can be detected in a single EoR quasar spectrum, allowing the measurement of the IGM neutral fraction towards that line of sight. However, damping wing features have only been detected in two $z>7$ quasars in previous studies. In this paper, we present new high quality optical and near-infrared spectroscopy of the $z=7.00$ quasar DES J025216.64--050331.8 obtained with Keck/NIRES and Gemini/GMOS. By using the MgII single-epoch virial method, we find that it hosts a $\rm (1.39\pm0.16) \times10^{9} ~M_\odot$ SMBH accreting at an Eddington ratio of $λ_{\rm Edd}=0.7\pm0.1$, consistent with the values seen in other luminous $z\sim 7$ quasars. Furthermore, the Ly$α$ region of the spectrum exhibits a strong damping wing absorption feature. The lack of associated metal absorption in the quasar spectrum indicates that this absorption is imprinted by a neutral IGM. Using a state-of-the-art model developed by Davies et al., we measure a volume-averaged neutral hydrogen fraction at $z=7$ of $\langle x_{\rm HI} \rangle = 0.70^{+0.20}_{-0.23} (^{+0.28}_{-0.48})$ within 68% (95%) confidence intervals when marginalizing over quasar lifetimes of $10^3\le t_{\rm Q}\le10^8$ yr. This is the highest IGM neutral fraction yet measured using reionization-era quasar spectra.

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

Deceptively cold dust in the massive starburst galaxy GN20 at $z\sim4$

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

preprint2020arXiv

Detecting and Characterizing Young Quasars I: Systemic Redshifts and Proximity Zones Measurements

In a multi-wavelength survey of $13$ quasars at $5.8\lesssim z\lesssim6.5$, that were pre-selected to be potentially young, we find five objects with extremely small proximity zone sizes that may imply UV-luminous quasar lifetimes of $\lesssim 100,000$ years. Proximity zones are regions of enhanced transmitted flux in the vicinity of the quasars that are sensitive to the quasars' lifetimes because the intergalactic gas has a finite response time to their radiation. We combine sub-mm observations from the Atacama Large Millimetre Array (ALMA) and the NOrthern Extended Millimeter Array (NOEMA), as well as deep optical and near-infrared spectra from medium-resolution spectrograph on the Very Large Telescope (VLT) and on the Keck telescopes, in order to identify and characterize these new young quasars, which provide valuable clues about the accretion behavior of supermassive black holes (SMBHs) in the early universe, and pose challenges on current black hole formation models to explain the rapid formation of billion solar mass black holes. We measure the quasars' systemic redshifts, black hole masses, Eddington ratios, emission line luminosities, and star formation rates of their host galaxies. Combined with previous results we estimate the fraction of young objects within the high-redshift quasar population at large to be $5\%\lesssim f_{\rm young}\lesssim 10\%$. One of the young objects, PSO J158-14, shows a very bright dust continuum flux ($F_{\rm cont}=3.46\pm 0.02\,\rm mJy$), indicating a highly star-bursting host galaxy with a star formation rate of approximately $1420\,M_{\odot}\,\rm yr^{-1}$.

preprint2020arXiv

Ionized and atomic interstellar medium in the z = 6.003 quasar SDSS J2310+1855

Observing the interstellar medium (ISM) in $z \gtrsim 6$ quasars host galaxies is essential for understanding the co-evolution between the supermassive black holes and their hosts. To probe the gas physical conditions and search for imprints of Active Galactic Nuclei (AGN) on the ISM, we report ALMA observations of the $\rm [N\ II]_{122 μm}$ and $\rm [O\ I]_{146 μm}$ lines and the underlying continuum from the $z=6.003$ quasar SDSS J231038.88+185519.7. Together with previous $\rm [C\ II]_{158 μm}$ and $\rm [O\ III]_{88 μm}$ observations, we use the ratios of these fine-structure lines to probe the ISM properties. Similar to other high-$z$ systems, this object exhibits a $\rm [C\ II]_{158 μm}$/$\rm [O\ I]_{146 μm}$ ratio comparable to the lowest values found in local (Ultra) luminous infrared galaxies, suggesting a "warmer" and "denser" gas component compared to typical local systems. The $\rm [O\ III]_{88 μm}$/$\rm [O\ I]_{146 μm}$ ratio is lower than that of other local and high-$z$ systems, indicating a smaller ionized gas fraction in this quasar. The $\rm [O\ III]_{88 μm}$/$\rm [N\ II]_{122 μm}$ ratio is comparable to that of local systems, and suggests a metallicity of $Z/Z_{\odot}$=1.5$-$2.1. Based on the $\rm [N\ II]_{122 μm}$ detection, we estimate that $17\%$ of the $\rm [C\ II]_{158 μm}$ emission is associated with ionized gas. The $\rm [N\ II]_{122 μm}$ line shows a "flux deficit" comparable to local systems. The $\rm [O\ I]_{146 μm}$ line, with a $\rm [O\ I]_{146 μm}$/FIR ratio $\ge 2\times$ than expected from the local relation, indicates no $\rm [O\ I]_{\rm 146 μm}$ deficit. The low $\rm [C\ II]_{158 μm}$/$\rm [O\ I]_{146 μm}$ ratio, together with the high $\rm [O\ I]_{146 μm}$/FIR ratio in J2310+1855, reveals that the warm and dense gas is likely a result of AGN heating to the ISM.

preprint2020arXiv

No Redshift Evolution in the Broad Line Region Metallicity up to z=7.54: Deep NIR Spectroscopy of ULAS J1342+0928

We present deep (9 hours) Gemini-N/GNIRS near-infrared spectroscopic observations of ULAS J1342+0928, a luminous quasar at z=7.54. Various broad emission lines were detected, as well as the underlying continuum and iron forests over the rest-frame wavelength 970-2930A. There is a clear trend that higher-ionization emission lines show larger blueshifts with CIV1549 exhibiting 5510^{+240}_{-110} km s-1 blueshift with respect to the systematic redshift from the far-infrared [CII] 158um emission line. Those high ionization lines have wide profiles with FWHM more than 10000 km s-1. A modest blueshift of 340^{+110}_{-80} km s-1 is also seen in MgII, the lowest ionization line identified in the spectrum. The updated MgII-based black hole mass of M_BH=9.1_{-1.3}^{+1.4} x 10^8 M_sun and the Eddington ratio of L_bol/L_Edd=1.1_{-0.2}^{+0.2} confirm that ULAS J1342+0928 is powered by a massive and actively accreting black hole. There is no significant difference in the emission line ratios such as SiIV/CIV and AlIII/CIV when compared to lower-redshift quasars in a similar luminosity range, which suggests early metal pollution of the broad-line-region clouds. This trend also holds for the FeII/MgII line ratio, known as a cosmic clock that traces the iron enrichment in the early universe. Different iron templates and continuum fitting ranges were used to explore how the FeII/MgII measurement changes as a function of spectral modeling. Quasars at even higher redshift or at fainter luminosity range (L_bol<10^46 erg s-1) are needed to probe the sites of early metal enrichment and a corresponding change in the FeII/MgII ratio.

preprint2020arXiv

Plateau de Bure High-z Blue-Sequence Survey 2 (PHIBSS2): Search for Secondary Sources, CO Luminosity Functions in the Field, and the Evolution of Molecular Gas Density through Cosmic Time

We report on the results of a search for serendipitous sources in CO emission in 110 cubes targeting CO(2-1), CO(3-2), and CO(6-5) at z ~ 1-2 from the second Plateau de Bure High-z Blue-Sequence Survey (PHIBSS2). The PHIBSS2 observations were part of a 4-year legacy program at the IRAM Plateau de Bure Interferometer aimed at studying early galaxy evolution from the perspective of molecular gas reservoirs. We present a catalog of 67 candidate secondary sources from this search, with 45 out of the 110 data cubes showing sources in addition to the primary target that appear to be field detections, unrelated to the central sources. This catalog includes the redshifts, line widths, fluxes, as well as an estimation of their reliability based on their false positive probability. We perform a search in the 3D-HST/CANDELS catalogs for the secondary CO detections and tentatively find that ~64% of these have optical counterparts, which we use to constrain their redshifts. Finally, we use our catalog of candidate CO detections to derive the CO(2-1), CO(3-2), CO(4-3), CO(5-4), and CO(6-5) luminosity functions over a range of redshifts, as well as the molecular gas mass density evolution. Despite the different methodology, these results are in very good agreement with previous observational constraints derived from blind searches in deep fields. They provide an example of the type of "deep field" science that can be carried out with targeted observations.

preprint2020arXiv

Pōniuā'ena: A Luminous $z=7.5$ Quasar Hosting a 1.5 Billion Solar Mass Black Hole

We report the discovery of a luminous quasar, J1007+2115 at $z=7.515$ ("Pōniuā'ena"), from our wide-field reionization-era quasar survey. J1007+2115 is the second quasar now known at $z>7.5$, deep into the reionization epoch. The quasar is powered by a $(1.5\pm0.2)\times10^9$ $M_{\odot}$ supermassive black hole (SMBH), based on its broad MgII emission-line profile from Gemini and Keck near-IR spectroscopy. The SMBH in J1007+2115 is twice as massive as that in quasar J1342+0928 at $z=7.54$, the current quasar redshift record holder. The existence of such a massive SMBH just 700 million years after the Big Bang significantly challenges models of the earliest SMBH growth. Model assumptions of Eddington-limited accretion and a radiative efficiency of 0.1 require a seed black hole of $\gtrsim 10^{4}$ $M_{\odot}$ at $z=30$. This requirement suggests either a massive black hole seed as a result of direct collapse or earlier periods of rapid black hole growth with hyper-Eddington accretion and/or a low radiative efficiency. We measure the damping wing signature imprinted by neutral hydrogen absorption in the intergalactic medium (IGM) on J1007+2115's Ly$α$ line profile, and find that it is weaker than that of J1342+0928 and two other $z\gtrsim7$ quasars. We estimate an IGM volume-averaged neutral fraction $\langle x\rm_{HI}\rangle=0.39^{+0.22}_{-0.13}$. This range of values suggests a patchy reionization history toward different IGM sightlines. We detect the 158 $μ$m [C II] emission line in J1007+2115 with ALMA; this line centroid yields a systemic redshift of $z=7.5149\pm0.0004$ and indicates a star formation rate of $\sim210$ $M_{\odot}$ yr$^{-1}$ in its host galaxy.

preprint2020arXiv

Probing the Nature of High Redshift Weak Emission Line Quasars: A Young Quasar with a Starburst Host Galaxy

We present the discovery of PSO J083.8371+11.8482, a weak emission line quasar with extreme star formation rate at $z=6.3401$. This quasar was selected from Pan-STARRS1, UHS, and unWISE photometric data. Gemini/GNIRS spectroscopy follow-up indicates a MgII-based black hole mass of $M_\mathrm{BH}=\left(2.0^{+0.7}_{-0.4}\right)\times10^9~M_\odot$ and an Eddington ratio of $L_\mathrm{bol}/L_\mathrm{Edd}=0.5^{+0.1}_{-0.2}$, in line with actively accreting supermassive black hole (SMBH) at $z\gtrsim6$. HST imaging sets strong constraint on lens-boosting, showing no relevant effect on the apparent emission. The quasar is also observed as a pure point-source with no additional emission component. The broad line region (BLR) emission is intrinsically weak and not likely caused by an intervening absorber. We found rest-frame equivalent widths of EW(Ly$α$+NV) $=5.7\pm0.7$ Angstrom, EW(CIV) $\leq5.8$ Angstrom (3-sigma upper limit), and EW(MgII) $=8.7\pm0.7$ Angstrom. A small proximity zone size ($R_\mathrm{p}=1.2\pm0.4$ pMpc) indicates a lifetime of only $t_\mathrm{Q}=10^{3.4\pm0.7}$ years from the last quasar phase ignition. ALMA shows extended [CII] emission with a mild velocity gradient. The inferred far-infrared luminosity ($L_\mathrm{FIR}=(1.2\pm0.1)\times10^{13}\,L_\odot$) is one of the highest among all known quasar hosts at $z\gtrsim6$. Dust and [CII] emissions put a constraint on the star formation rate of SFR $=900-4900~M_\odot\,\mathrm{yr^{-1}}$, similar to that of hyper-luminous infrared galaxy. Considering the observed quasar lifetime and BLR formation timescale, the weak-line profile in the quasar spectrum is most likely caused by a BLR which is not yet fully formed rather than continuum boosting by gravitational lensing or a soft continuum due to super-Eddington accretion.

preprint2020arXiv

Revealing the Accretion Physics of Supermassive Black Holes at Redshift z~7 with Chandra and Infrared Observations

X-ray emission from quasars has been detected up to redshift $z=7.5$, although only limited to a few objects at $z>6.5$. In this work, we present new Chandra observations of five $z>6.5$ quasars. By combining with archival Chandra observations of six additional $z>6.5$ quasars, we perform a systematic analysis on the X-ray properties of these earliest accreting supermassive black holes (SMBHs). We measure the black hole masses, bolometric luminosities ($L_{\rm bol}$), Eddington ratios ($λ_{\rm Edd}$), emission line properties, and infrared luminosities ($L_{\rm IR}$) of these quasars using infrared and sub-millimeter observations. Correlation analysis indicates that the X-ray bolometric correction (the factor that converts from X-ray luminosity to bolometric luminosity) decreases with increasing $L_{\rm bol}$, and that the UV/optical-to-X-ray ratio, $α_{\rm ox}$, strongly correlates with $L_{\rm 2500}$, and moderately correlates with $λ_{\rm Edd}$ and blueshift of CIV emission lines. These correlations are consistent with those found in lower-$z$ quasars, indicating quasar accretion physics does not evolve with redshift. We also find that $L_{\rm IR}$ does not correlate with $L_{\rm 2-10 keV}$ in these luminous distant quasars, suggesting that the ratio of the SMBH growth rate and their host galaxy growth rate in these early luminous quasars are different from those of local galaxies. A joint spectral analysis of the X-ray detected $z>6.5$ quasars yields an average X-ray photon index of $ Γ=2.32^{+0.31}_{-0.30}$, steeper than that of low-$z$ quasars. By comparing it with the $Γ-λ_{\rm Edd}$ relation, we conclude that the steepening of $Γ$ for quasars at $z>6.5$ is mainly driven by their higher Eddington ratios.

preprint2020arXiv

Strong Mg II and Fe II Absorbers at 2.2 < z < 6.0

We present a study of strong intervening absorption systems in the near-IR spectra of 31 luminous quasars at $z>5.7$. The quasar spectra were obtained with {\it Gemini} GNIRS that provide continuous wavelength coverage from $\sim$0.9 to $\sim$2.5 $μ$m. We detect 32 strong Mg II doublet absorbers with rest-frame equivalent width $W_r$ ($\lambda2796$) $>1.0$ Åat $2.2 < z < 6.0$. Each Mg II absorber is confirmed by at least two associated Fe II absorption lines in the rest-frame wavelength range of $\sim 1600-2600$ Å. We find that the comoving line density ($dN/dX$) of the strong Fe II-bearing Mg II absorbers decreases towards higher redshift at $z>3$, consistent with previous studies. Compared with strong Mg II absorbers detected in damped Ly$α$ systems at 2 $<z<$ 4, our absorbers are potentially less saturated and show much larger rest-frame velocity widths. This suggests that the gas traced by our absorbers are potentially affected by galactic superwinds. We analyze the {\it Hubble Space Telescope} near-IR images of the quasars and identify possible associated galaxies for our strong absorbers. There are a maximum of two galaxy candidates found within 5" radius of each absorber. The median F105W-band magnitude of these galaxy candidates is 24.8 mag, which is fainter than the $L^*$ galaxy luminosity at $z\sim$ 4. By using our observed $dN/dX$ of strong Mg II absorbers and galaxy candidates median luminosity, we suggest that at high redshift, strong Mg II absorbers tend to have a more disturbed environment but smaller halo size than that at $z <$ 1.

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 Molecular ISM in the Super Star Clusters of the Starburst NGC 253

We present submillimeter spectra of the (proto-)super star cluster (SSC) candidates in the starbursting center of the nearby galaxy NGC 253 identified by Leroy et al. (2018). The 2.5pc resolution of our ALMA cycle 3 observations approach the size of the SSCs and allows the study of physical and chemical properties of the molecular gas in these sources. In the 14 SSC sources and in the frequency ranges 342.0-345.8 GHz and 353.9-357.7 GHz we detect 55 lines belonging to 19 different chemical species. The SSCs differ significantly in chemical complexity, with the richest clusters showing 19 species and the least complex showing 4 species. We detect HCN isotopologues and isomers (H$^{13}$CN, HC$^{15}$N, H$^{15}$NC), abundant HC$_3$N, SO and S$^{18}$O, SO$_2$, and H$_2$CS. The gas ratios CO/HCN, CO/HCO$^+$ are low, ~1-10, implying high dense gas fractions in the SSCs. Line ratio analyses suggests chemistry consistent with photon-dominated regions and mechanical heating. None of the SSCs near the galaxy center show line ratios that imply an X-ray dominated region, suggesting that heating by any (still unknown) AGN does not play a major role. The gas temperatures are high in most sources, with an average rotational temperature of ~130 K in SO$_2$. The widespread existence of vibrationally excited HCN and HC$_3$N transitions implies strong IR radiation fields, potentially trapped by a greenhouse effect due to high continuum opacities.

preprint2020arXiv

The RESOLVE project: a multi-physics experiment with a temporary dense seismic array on the Argentière Glacier, French Alps

Recent work in the field of cryo-seismology demonstrates that high frequency (>1 Hz) waves provide key constraints on a wide range of glacier processes such as basal friction, surface crevassing or subglacial water flow. Establishing quantitative links between the seismic signal and the processes of interest however requires detailed characterization of the wavefield, which at the high frequencies of interest necessitates the deployment of large and particularly dense seismic arrays. Although dense seismic array monitoring has recently become routine in geophysics, its application to glaciated environments has yet remained limited. Here we present a dense seismic array experiment made of 98 3-component seismic stations continuously recording during 35 days in early spring on the Argentière Glacier, French Alps. The seismic dataset is supplemented by a wide range of complementary observations obtained from ground penetrating radar, drone imagery, GNSS positioning and in-situ instrumentation of basal glacier sliding velocities and subglacial water discharge. Through applying multiple processing techniques including event detection from template matching and systematic beamforming, we demonstrate that the present dataset provides enhanced spatial resolution on basal stick slip and englacial fracturing sources as well as novel constraints on heterogeneous nature of the noise field generated by subglacial water flow and on the link between crevasse properties and englacial seismic velocities. We finally outline in which ways further work using this dataset could help tackle key remaining questions in the field.

preprint2020arXiv

The turbulent gas structure in the centers of NGC253 and the Milky Way

We compare molecular gas properties in the starbursting center of NGC253 and the Milky Way Galactic Center (GC) on scales of ~1-100 pc using dendograms and resolution-, area- and noise-matched datasets in CO (1-0) and CO (3-2). We find that the size-line width relations in NGC253 and the GC have similar slope, but NGC253 has larger line widths by factors of ~2-3. The $σ^2/R$ dependency on column density shows that, in the GC, on scales of 10-100 pc the kinematics of gas over $N>3\times10^{21}$ cm$^{-2}$ are compatible with gravitationally bound structures. In NGC253 this is only the case for column densities $N>3\times10^{22}$ cm$^{-2}$. The increased line widths in NGC253 originate in the lower column density gas. This high-velocity dispersion, not gravitationally self-bound gas is likely in transient structures created by the combination of high average densities and feedback in the starburst. The high densities turns the gas molecular throughout the volume of the starburst, and the injection of energy and momentum by feedback significantly increases the velocity dispersion at a given spatial scale over what is observed in the GC.

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.

preprint2020arXiv

X-ray Observations of a [C II]-bright, z=6.59 Quasar/Companion System

We present deep Chandra observations of PSO J231.6576$-$20.8335, a quasar at redshift z=6.59 with a nearby (${\sim}8$ proper kpc) companion galaxy. ALMA observed both the quasar and companion to be bright in [C II], and the system has significant extended Ly$α$ emission around the quasar, suggesting that a galaxy merger is ongoing. Unlike previous studies of two similar systems, and despite observing the system with Chandra for 140 ks, we do not detect the companion in X-rays. The quasar itself is detected, but only $13.3^{+4.8}_{-3.7}$ net counts are observed. From a basic spectral analysis, the X-ray spectrum of the quasar is soft (hardness ratio of $\mathcal{HR} = -0.60_{-0.27}^{+0.17}$, power-law index of $Γ=2.6^{+1.0}_{-0.9}$), which results in a rest-frame X-ray luminosity comparable to other bright quasars ($L_{2-10} = 1.09^{+2.20}_{-0.70}\times 10^{45}\ \textrm{erg}\ \textrm{s}^{-1}$) despite the faint observed X-ray flux. We highlight two possible interpretations of this result: the quasar has a steep value of $Γ$ -- potentially related to observed ongoing Eddington accretion -- thereby pushing much of the emission out of our observed band, or the quasar has a more normal spectrum ($Γ{\sim}2$) but is therefore less X-ray luminous ($L_{2-10} \sim 0.6 \times 10^{45}\ \textrm{ erg}\ \textrm{ s}^{-1}$).

preprint2019arXiv

Probing the full CO spectral line energy distribution (SLED) in the nuclear region of a quasar-starburst system at $z=6.003$

We report Atacama Large Millimeter/submillimeter Array (ALMA) observations of CO $(8-7)$, $(9-8)$, $\rm H_{2}O (2_{0,2}-1_{1,1})$ and $\rm OH^{+} (1_{1}-0_{1})$ and NOrthern Extended Millimeter Array (NOEMA) observations of CO $(5-4)$, $(6-5)$, $(12-11)$ and $(13-12)$ towards the $z = 6.003$ quasar SDSS J231038.88+185519.7, aiming to probe the physical conditions of the molecular gas content of this source. We present the best sampled CO spectral line energy distribution (SLED) at $z = 6.003$, and analyzed it with the radiative transfer code MOLPOP-CEP. Fitting the CO SLED to a one-component model indicates a kinetic temperature $T_{\rm kin} = 228 \ \rm K$, molecular gas density $log (n(\rm H_{2})/\rm cm^{-3}$ )=4.75, and CO column density $log(N(\rm CO)/\rm cm^{-2}) =17.5$, although a two-component model better fits the data. In either case, the CO SLED is dominated by a "warm" and "dense" component. Compared to samples of local (Ultra) Luminous Infrared Galaxies ((U)LIRGs), starburst galaxies and high redshift Submillimeter Galaxies (SMGs), J2310+1855 exhibits higher CO excitation at ($J \geq 8$), like other high redshift quasars. The high CO excitation, together with the enhanced $L_{\rm H_{2}O}/ L_{IR} $, $L_{\rm H_{2}O}/ L_{CO} $ and $L_{OH^{+}}/L_{\rm H_{2}O} $ ratios, suggests that besides the UV radiation from young massive stars, other mechanisms such as shocks, cosmic rays and X-rays might also be responsible for the heating and ionization of the molecular gas. In the nuclear region probed by the molecular emissions lines, any of these mechanisms might be present due to the powerful quasar and the starburst activity.

preprint2019arXiv

Resolving the interstellar medium in the nuclear region of two z=5.78 quasar host galaxies with ALMA

We present ALMA observations of the [CII] 158 $μ$m fine structure line and dust continuum emission from two quasars, SDSS J104433.04-012502.2 and SDSS J012958.51-003539.7, at z=5.78. The ALMA observations at 0.2'' resolution map the dust and gas on kpc scales. The spatially resolved emission shows a similar trend of decreasing [CII]-FIR ratios with increasing FIR surface brightnesses as was found in the infrared luminous galaxies with intense star formation. We confirm the velocity gradients of [CII] emission found previously in SDSS J0129-0035. No clear evidence of order motion is detected in SDSS J1044-0125. The velocity maps and PV-diagrams also suggest turbulent gas clumps in both objects. We tentatively detect a [CII] peak offset 4.9 kpc to the East of SDSS J1044-0125. This may be associated with an infalling companion, or node of gas outflow. All these results suggest significant dynamical evolution of the ISM in the nuclear region of these young quasar-starburst systems. We fit the velocity map of the [CII] emission from SDSS J0129-0035 with a rotating disk model. The result suggests a face-on system with an inclination angle of 16 +/- 20 degree and constrains the lower limit of the host galaxy dynamical mass to be 2.6x10^10 Msun within the [CII] emitting region. It is likely that SDSS J0129-0035, as well as other young quasars with super massive black hole masses on order of 10^7 to 10^8 Msun, falls close to the black hole and host galaxy mass relation defined by local galaxies.

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 REQUIEM Survey I: A Search for Extended Ly-Alpha Nebular Emission Around 31 z>5.7 Quasars

The discovery of quasars few hundred megayears after the Big Bang represents a major challenge to our understanding of black holes and galaxy formation and evolution. Their luminosity is produced by extreme gas accretion onto black holes, which already reached masses of 10$^9$ M$_\odot$ by z ~ 6. Simultaneously, their host galaxies form hundreds of stars per year, using up gas in the process. To understand which environments are able to sustain the rapid formation of these extreme sources we started a VLT/MUSE effort aimed at characterizing the surroundings of a sample of 5.7 < z < 6.6 quasars dubbed: the Reionization Epoch QUasar InvEstigation with MUSE (REQUIEM) survey. We here present results of our searches for extended Ly-Alpha halos around the first 31 targets observed as part of this program. Reaching 5-sigma surface brightness limits of 0.1-1.1 x 10$^{-17}$ erg/s/cm$^2$/arcsec$^2$ over a 1 arcsec$^2$ aperture, we were able to unveil the presence of 12 Ly-Alpha nebulae, 8 of which are newly discovered. The detected nebulae show a variety of emission properties and morphologies with luminosities ranging from 8 x 10$^{42}$ to 2 x 10$^{44}$ erg/s, FWHMs between 300 and 1700 km/s, sizes < 30 pkpc, and redshifts consistent with those of the quasar host galaxies. As the first statistical and homogeneous investigation of the circum-galactic medium of massive galaxies at the end of the reionization epoch, the REQUIEM survey enables the study of the evolution of the cool gas surrounding quasars in the first 3 Gyr of the Universe. A comparison with the extended Ly-Alpha emission observed around bright (M$_{1450}$ < -25 mag) quasars at intermediate redshift indicates little variations on the properties of the cool gas from z ~ 6 to z ~ 3 followed by a decline in the average surface brightness down to z ~ 2.

preprint2016arXiv

A Total Molecular Gas Mass Census in z~2-3 Star-forming Galaxies: Low-J CO Excitation Probes of Galaxies' Evolutionary States

We present CO(1-0) observations obtained at the Karl G. Jansky Very Large Array (VLA) for 14 z~2 galaxies with existing CO(3-2) measurements, including 11 galaxies which contain active galactic nuclei (AGN) and three submillimeter galaxies (SMGs). We combine this sample with an additional 15 z~2 galaxies from the literature that have both CO(1-0) and CO(3-2) measurements in order to evaluate differences in CO excitation between SMGs and AGN host galaxies, measure the effects of CO excitation on the derived molecular gas properties of these populations, and to look for correlations between the molecular gas excitation and other physical parameters. With our expanded sample of CO(3-2)/CO(1-0) line ratio measurements, we do not find a statistically significant difference in the mean line ratio between SMGs and AGN host galaxies as found in the literature, instead finding r_3,1=1.03+/-0.50 for AGN host galaxies and r_3,1=0.78+/-0.27 for SMGs (or r_3,1=0.90+/-0.40 for both populations combined). We also do not measure a statistically significant difference between the distributions of the line ratios for these populations at the p=0.05 level, although this result is less robust. We find no excitation dependence on the index or offset of the integrated Schmidt-Kennicutt relation for the two CO lines, and obtain indices consistent with N=1 for the various sub-populations. However, including low-z "normal" galaxies increases our best-fit Schmidt-Kennicutt index to N~1.2. While we do not reproduce correlations between the CO line width and luminosity, we do reproduce correlations between CO excitation and star formation efficiency.

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

preprint2016arXiv

High Resolution Observations of Molecular Lines in Arp 220: Kinematics, Morphology, and Limits on the Applicability of the Ammonia Thermometer

We observe Arp 220, the nearest Ultra-Luminous Infrared Galaxy (ULIRG), over 4 GHz in the K and Ka bands. We provide constraints for the kinematics,morphology, and identify molecular species on scales resolving both nuclei (0.6" or 230 pc). We detect multiple molecular species, including hydroxyl in both cores. We tentatively detect H2O at 21.84 GHz in both nuclei, indicating the likely presence of maser emission. The observed frequency range also contains metastable ammonia transitions from (J,K) = (1,1) to (5,5), as well as the (9,9) inversion line, which, together are a well-known thermometer of dense molecular gas. Furthermore, the non-metastable (4,2) and (10,9) and possibly the (3,1) lines are also detected. We apply a standard temperature analysis to Arp 220. However, the analysis is complicated in that standard LTE assumptions do not hold. There are indications that a substantial fraction of ammonia could be in the non-metastable transitions as opposed to only the metastable ones. Thus, the non-metastable transitions could be essential to constraining the temperature. We compare all of these data to ALMA observations of this source, confirming the outflow previously observed by other tracers in both nuclei.

preprint2016arXiv

The Spatially Resolved [CII] Cooling Line Deficit in Galaxies

We present [CII] 158um measurements from over 15,000 resolved regions within 54 nearby galaxies of the KINGFISH program to investigate the so-called [CII] "line cooling deficit" long known to occur in galaxies with different luminosities. The [CII]/TIR ratio ranges from above 1% to below 0.1% in the sample, with a mean value of 0.48+-0.21%. We find that the surface density of 24um emission dominates this trend, with [CII]/TIR dropping as nuInu{24um} increases. Deviations from this overall decline are correlated with changes in the gas phase metal abundance, with higher metallicity associated with deeper deficits at a fixed surface brightness. We supplement the local sample with resolved [CII] measurements from nearby luminous infrared galaxies and high redshift sources from z=1.8-6.4, and find that star formation rate density drives a continuous trend of deepening [CII] deficit across six orders of magnitude in SFRD. The tightness of this correlation suggests that an approximate star formation rate density can be estimated directly from global measurements of [CII]/TIR, and a relation is provided to do so. Several low-luminosity AGN hosts in the sample show additional and significant central suppression of [CII]/TIR, but these deficit enhancements occur not in those AGN with the highest X-ray luminosities, but instead those with the highest central starlight intensities. Taken together, these results demonstrate that the [CII] cooling line deficit in galaxies likely arises from local physical phenomena in interstellar gas.

preprint2015arXiv

ALMA Multi-line Imaging of the Nearby Starburst Galaxy NGC 253

We present spatially resolved ($\sim$50 pc) imaging of molecular gas species in the central kiloparsec of the nearby starburst galaxy NGC 253, based on observations taken with the Atacama Large Millimeter/submillimeter Array (ALMA). A total of 50 molecular lines are detected over a 13 GHz bandwidth imaged in the 3 mm band. Unambiguous identifications are assigned for 27 lines. Based on the measured high CO/C$^{17}$O isotopic line ratio ($\gtrsim$350), we show that $^{12}$CO(1-0) has moderate optical depths. A comparison of the HCN and HCO$^{+}$ with their $^{13}$C-substituted isotopologues shows that the HCN(1-0) and HCO$^{+}$(1-0) lines have optical depths at least comparable to CO(1-0). H$^{13}$CN/H$^{13}$CO$^{+}$ (and H$^{13}$CN/HN$^{13}$C) line ratios provide tighter constraints on dense gas properties in this starburst. SiO has elevated abundances across the nucleus. HNCO has the most distinctive morphology of all the bright lines, with its global luminosity dominated by the outer parts of the central region. The dramatic variation seen in the HNCO/SiO line ratio suggests that some of the chemical signatures of shocked gas are being erased in the presence of dominating central radiation fields (traced by C$_{2}$H and CN). High density molecular gas tracers (including HCN, HCO$^+$, and CN) are detected at the base of the molecular outflow. We also detect hydrogen $β$ recombination lines that, like their $α$ counterparts, show compact, centrally peaked morphologies, distinct from the molecular gas tracers. A number of sulfur based species are mapped (CS, SO, NS, C$_{2}$S, H$_{2}$CS and CH$_{3}$SH) and have morphologies similar to SiO.

preprint2015arXiv

An ALMA survey of Sub-millimeter Galaxies in the Extended Chandra Deep Field South: Physical properties derived from ultraviolet-to-radio modelling

[abridged] The ALESS survey has followed-up a sample of 122 sub-millimeter sources in the Extended Chandra Deep Field South at 870um with ALMA, allowing to pinpoint the positions of sub-millimeter galaxies (SMGs) to 0.3'' and to find their precise counterparts at different wavelengths. This enabled the first compilation of the multi-wavelength spectral energy distributions (SEDs) of a statistically reliable survey of SMGs. In this paper, we present a new calibration of the MAGPHYS modelling code that is optimized to fit these UV-to-radio SEDs of z>1 star-forming galaxies using an energy balance technique to connect the emission from stellar populations, dust attenuation and dust emission in a physically consistent way. We derive statistically and physically robust estimates of the photometric redshifts and physical parameters for the ALESS SMGs. We find that they have a median stellar mass $M_\ast=(8.9\pm0.1)\times10^{10} M_\odot$, SFR$=280\pm70 M_\odot$/yr, overall V-band dust attenuation $A_V=1.9\pm0.2$ mag, dust mass $M_\rm{dust}=(5.6\pm1.0)\times10^8 M_\odot$, and average dust temperature Tdust~40 K. The average intrinsic SED of the ALESS SMGs resembles that of local ULIRGs in the IR range, but the stellar emission of our average SMG is brighter and bluer, indicating lower dust attenuation, possibly because they are more extended. We explore how the average SEDs vary with different parameters, and we provide a new set of SMG templates. To put the ALESS SMGs into context, we compare their stellar masses and SFRs with those of less actively star-forming galaxies at the same redshifts. At z~2, about half of the SMGs lie above the star-forming main sequence, while half are at the high-mass end of the sequence. At higher redshifts (z~3.5), the SMGs tend to have higher SFR and Mstar, but the fraction of SMGs that lie significantly above the main sequence decreases to less than a third.

preprint2015arXiv

Discovery of Large Molecular Gas Reservoirs in Post-Starburst Galaxies

Post-starburst (or "E+A") galaxies are characterized by low H$α$ emission and strong Balmer absorption, suggesting a recent starburst, but little current star formation. Although many of these galaxies show evidence of recent mergers, the mechanism for ending the starburst is not yet understood. To study the fate of the molecular gas, we search for CO (1-0) and (2-1) emission with the IRAM 30m and SMT 10m telescopes in 32 nearby ($0.01<z<0.12$) post-starburst galaxies drawn from the Sloan Digital Sky Survey. We detect CO in 17 (53%). Using CO as a tracer for molecular hydrogen, and a Galactic conversion factor, we obtain molecular gas masses of $M(H_2)=10^{8.6}$-$10^{9.8} M_\odot$ and molecular gas mass to stellar mass fractions of $\sim10^{-2}$-$10^{-0.5}$, comparable to those of star-forming galaxies. The large amounts of molecular gas rule out complete gas consumption, expulsion, or starvation as the primary mechanism that ends the starburst in these galaxies. The upper limits on $M(H_2)$ for the 15 undetected galaxies range from $10^{7.7} M_\odot$ to $10^{9.7} M_\odot$, with the median more consistent with early-type galaxies than with star-forming galaxies. Upper limits on the post-starburst star formation rates (SFRs) are lower by $\sim10\times$ than for star-forming galaxies with the same $M(H_2)$. We also compare the molecular gas surface densities ($Σ_{\rm H_2}$) to upper limits on the SFR surface densities ($Σ_{\rm SFR}$), finding a significant offset, with lower $Σ_{\rm SFR}$ for a given $Σ_{\rm H_2}$ than is typical for star-forming galaxies. This offset from the Kennicutt-Schmidt relation suggests that post-starbursts have lower star formation efficiency, a low CO-to-H$_2$ conversion factor characteristic of ULIRGs, and/or a bottom-heavy initial mass function, although uncertainties in the rate and distribution of current star formation remain.

preprint2015arXiv

Faint CO Line Wings in Four Star-Forming (Ultra)luminous Infrared Galaxies

We report the results of a search for large velocity width, low-intensity line wings - a commonly used signature of molecular outflows - in four low redshift (ultra)luminous infrared galaxies (U/LIRGs) that appear to be dominated by star formation. The targets were drawn from a sample of fourteen such galaxies presented in Chung et al. (2011), who showed the stacked CO spectrum of the sample to exhibit 1000 km/s-wide line wings. We obtained sensitive, wide bandwidth imaging of our targets using the IRAM Plateau de Bure Interferometer. We detect each target at very high significance but do not find the claimed line wings in these four targets. Instead, we constrain the flux in the line wings to be only a few percent. Casting our results as mass outflow rates following Cicone et al. (2014) we show them to be consistent with a picture in which very high mass loading factors preferentially occur in systems with high AGN contributions to their bolometric luminosity. We identify one of our targets, IRAS05083 (VII Zw 31), as a candidate molecular outflow.

preprint2015arXiv

High-resolution mass models of dwarf galaxies from LITTLE THINGS

We present high-resolution rotation curves and mass models of 26 dwarf galaxies from LITTLE THINGS. LITTLE THINGS is a high-resolution Very Large Array HI survey for nearby dwarf galaxies in the local volume within 11 Mpc. The rotation curves of the sample galaxies derived in a homogeneous and consistent manner are combined with Spitzer archival 3.6 micron and ancillary optical U, B, and V images to construct mass models of the galaxies. We decompose the rotation curves in terms of the dynamical contributions by baryons and dark matter halos, and compare the latter with those of dwarf galaxies from THINGS as well as Lambda CDM SPH simulations in which the effect of baryonic feedback processes is included. Being generally consistent with THINGS and simulated dwarf galaxies, most of the LITTLE THINGS sample galaxies show a linear increase of the rotation curve in their inner regions, which gives shallower logarithmic inner slopes alpha of their dark matter density profiles. The mean value of the slopes of the 26 LITTLE THINGS dwarf galaxies is alpha =-0.32 +/- 0.24 which is in accordance with the previous results found for low surface brightness galaxies (alpha = -0.2 +/- 0.2) as well as the seven THINGS dwarf galaxies (alpha =-0.29 +/- 0.07). However, this significantly deviates from the cusp-like dark matter distribution predicted by dark-matter-only Lambda CDM simulations. Instead our results are more in line with the shallower slopes found in the Lambda CDM SPH simulations of dwarf galaxies in which the effect of baryonic feedback processes is included. In addition, we discuss the central dark matter distribution of DDO 210 whose stellar mass is relatively low in our sample to examine the scenario of inefficient supernova feedback in low mass dwarf galaxies predicted from recent Lambda SPH simulations of dwarf galaxies where central cusps still remain.

preprint2015arXiv

Imaging the cold molecular gas in SDSS J1148 + 5251 at z = 6.4

We present Karl G. Jansky Very Large Array (VLA) observations of the CO ($J = 2 \rightarrow 1$) line emission towards the $z = 6.419$ quasar SDSS J$114816.64+525150.3$ (J$1148+5251$). The molecular gas is found to be marginally resolved with a major axis of $0.9"$ (consistent with previous size measurements of the CO ($J = 7 \rightarrow 6$) emission). We observe tentative evidence for extended line emission towards the south west on a scale of ~$1.4"$, but this is only detected at $3.3σ$ significance and should be confirmed. The position of the molecular emission region is in excellent agreement with previous detections of low frequency radio continuum emission as well as [C ii] line and thermal dust continuum emission. These CO ($J = 2 \rightarrow 1$) observations provide an anchor for the low excitation part of the molecular line SED. We find no evidence for extended low excitation component, neither in the spectral line energy distribution nor the image. We fit a single kinetic gas temperature model of 50 K. We revisit the gas and dynamical masses in light of this new detection of a low order transition of CO, and confirm previous findings that there is no extended reservoir of cold molecular gas in J$1148+5251$, and that the source departs substantially from the low $z$ relationship between black hole mass and bulge mass. Hence, the characteristics of J$1148+5251$ at $z = 6.419$ are very similar to $z$~$2$ quasars, in the lack of a diffuse cold gas reservoir and kpc-size compactness of the star forming region.

preprint2015arXiv

Next Generation Very Large Array Memo No. 8 Science Working Group 3: Galaxy Assembly through Cosmic Time

The Next-Generation Very Large Array (ngVLA) will be critical for understanding how galaxies are built and evolve at the earliest epochs. The sensitivity and frequency coverage will allow for the detection of cold gas and dust in `normal' distant galaxies, including the low-J transitions of molecular gas tracers such as CO, HNC, and HCO+; synchrotron and free-free continuum emission; and even the exciting possibility of thermal dust emission at the highest (z~7) redshifts. In particular, by enabling the total molecular gas reservoirs to be traced to unprecedented sensitivities across a huge range of epochs simultaneously -- something no other radio or submillimeter facility will be capable of -- the detection of the crucial low-J transitions of CO in a diverse body of galaxies will be the cornerstone of ngVLA's contribution to high-redshift galaxy evolution science. The ultra-wide bandwidths will allow a complete sampling of radio SEDs, as well as the detection of emission lines necessary for spectroscopic confirmation of elusive dusty starbursts. The ngVLA will also deliver unique contributions to our understanding of cosmic magnetism and to science accessible through microwave polarimetry. Finally, the superb angular resolution will move the field beyond detection experiments and allow detailed studies of the morphology and dynamics of these systems, including dynamical modeling of disks/mergers, determining the properties of outflows, measuring black hole masses from gas disks, and resolving multiple AGN nuclei. We explore the contribution of a ngVLA to these areas and more, as well as synergies with current and upcoming facilities including ALMA, SKA, large single-dish submillimeter observatories, GMT/TMT, and JWST.

preprint2015arXiv

Spatially extended and high-velocity dispersion molecular component in spiral galaxies: single-dish vs. interferometric observations

Recent studies of the molecular medium in nearby galaxies have provided mounting evidence that the molecular gas can exist in two phases: one that is clumpy and organized as molecular clouds and another one that is more diffuse. This last component has a higher velocity dispersion than the clumpy one. In order to investigate these two molecular components further, we compare the fluxes and line widths of CO in NGC 4736 and NGC 5055, two nearby spiral galaxies for which high-quality interferometric as well as single-dish data sets are available. Our analysis leads to two main results: 1) Employing three different methods, we determine the flux recovery of the interferometer as compared to the single-dish to be within a range of 35-74% for NGC4736 and 81-92% for NGC5055, and 2) when focusing on high (SNR>5) lines of sight, the single-dish line widths are larger by ~(40+-20)% than the ones derived from interferometric data; which is in agreement with stacking all lines of sight. These results point to a molecular gas component that is distributed over spatial scales larger than 30"(~1kpc), and is therefore filtered out by the interferometer. The available observations do not allow us to distinguish between a truly diffuse gas morphology and a uniform distribution of small clouds that are separated by less than the synthesized beam size (~3" or ~100pc), as they would both be invisible for the interferometer. This high velocity dispersion component has a dispersion similar to what is found in the atomic medium, as traced through observations of the HI line.

preprint2015arXiv

Submillimeter Observations of CLASH 2882 and the Evolution of Dust in this Galaxy

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

preprint2015arXiv

The Multi-Phase Cold Fountain in M82 Revealed by a Wide, Sensitive Map of the Molecular ISM

We present a wide area (~ 8 x 8 kpc), sensitive map of CO (2-1) emission around the nearby starburst galaxy M82. Molecular gas extends far beyond the stellar disk, including emission associated with the well-known outflow as far as 3 kpc from M82's midplane. Kinematic signatures of the outflow are visible in both the CO and HI emission: both tracers show a minor axis velocity gradient and together they show double peaked profiles, consistent with a hot outflow bounded by a cone made of a mix of atomic and molecular gas. Combining our CO and HI data with observations of the dust continuum, we study the changing properties of the cold outflow as it leaves the disk. While H_2 dominates the ISM near the disk, the dominant phase of the cool medium changes as it leaves the galaxy and becomes mostly atomic after about a kpc. Several arguments suggest that regardless of phase, the mass in the cold outflow does not make it far from the disk; the mass flux through surfaces above the disk appears to decline with a projected scale length of ~ 1-2 kpc. The cool material must also end up distributed over a much wider angle than the hot outflow based on the nearly circular isophotes of dust and CO at low intensity and the declining rotation velocities as a function of height from the plane. The minor axis of M82 appears so striking at many wavelengths because the interface between the hot wind cavity and the cool gas produces Halpha, hot dust, PAH emission, and scattered UV light. We also show the level at which a face-on version of M82 would be detectable as an outflow based on unresolved spectroscopy. Finally, we consider multiple constraints on the CO-to-H$_2$ conversion factor, which must change across the galaxy but appears to be only a factor of ~ 2 lower than the Galactic value in the outflow.

preprint2015arXiv

The radial variation of HI velocity dispersions in dwarfs and spirals

Gas velocity dispersions provide important diagnostics of the forces counteracting gravity to prevent collapse of the gas. We use the 21 cm line of neutral atomic hydrogen (HI) to study HI velocity dispersion and HI phases as a function of galaxy morphology in 22 galaxies from The HI Nearby Galaxy Survey (THINGS). We stack individual HI velocity profiles and decompose them into broad and narrow Gaussian components. We study the HI velocity dispersion and the HI surface density, as a function of radius. For spirals, the velocity dispersions of the narrow and broad components decline with radius and their radial profiles are well described by an exponential function. For dwarfs, however, the profiles are much flatter. The single Gaussian dispersion profiles are, in general, flatter than those of the narrow and broad components. In most cases, the dispersion profiles in the outer disks do not drop as fast as the star formation profiles, derived in the literature. This indicates the importance of other energy sources in driving HI velocity dispersion in the outer disks. The radial surface density profiles of spirals and dwarfs are similar. The surface density profiles of the narrow component decline more steeply than those of the broad component, but not as steep as what was found previously for the molecular component. As a consequence, the surface density ratio between the narrow and broad components, an estimate of the mass ratio between cold HI and warm HI, tends to decrease with radius. On average, this ratio is lower in dwarfs than in spirals. This lack of a narrow, cold HI component in dwarfs may explain their low star formation activity.

preprint2015arXiv

The Time Domain Spectroscopic Survey: Variable Object Selection and Anticipated Results

We present the selection algorithm and anticipated results for the Time Domain Spectroscopic Survey (TDSS). TDSS is an SDSS-IV eBOSS subproject that will provide initial identification spectra of approximately 220,000 luminosity-variable objects (variable stars and AGN) across 7,500 square degrees selected from a combination of SDSS and multi-epoch Pan-STARRS1 photometry. TDSS will be the largest spectroscopic survey to explicitly target variable objects, avoiding pre-selection on the basis of colors or detailed modeling of specific variability characteristics. Kernel Density Estimate (KDE) analysis of our target population performed on SDSS Stripe 82 data suggests our target sample will be 95% pure (meaning 95% of objects we select have genuine luminosity variability of a few magnitudes or more). Our final spectroscopic sample will contain roughly 135,000 quasars and 85,000 stellar variables, approximately 4,000 of which will be RR Lyrae stars which may be used as outer Milky Way probes. The variability-selected quasar population has a smoother redshift distribution than a color-selected sample, and variability measurements similar to those we develop here may be used to make more uniform quasar samples in large surveys. The stellar variable targets are distributed fairly uniformly across color space, indicating that TDSS will obtain spectra for a wide variety of stellar variables including pulsating variables, stars with significant chromospheric activity, cataclysmic variables and eclipsing binaries. TDSS will serve as a pathfinder mission to identify and characterize the multitude of variable objects that will be detected photometrically in even larger variability surveys such as LSST.

preprint2015arXiv

Variations in the Star Formation Efficiency of the Dense Molecular Gas across the Disks of Star-Forming Galaxies

We present a new survey of HCN(1-0) emission, a tracer of dense molecular gas, focused on the little-explored regime of normal star-forming galaxy disks. Combining HCN, CO, and infrared (IR) emission, we investigate the role of dense gas in Star Formation (SF), finding systematic variations in both the apparent dense gas fraction and the apparent SF efficiency (SFE) of dense gas. The latter may be unexpected, given the popularity of gas density threshold models to explain SF scaling relations. We used the IRAM 30-m telescope to observe HCN(1-0) across 29 nearby disk galaxies whose CO(2-1) emission has previously been mapped by the HERACLES survey. Because our observations span a range of galactocentric radii, we are able to investigate the properties of the dense gas as a function of local conditions. We focus on how the IR/CO, HCN/CO, and IR/HCN ratios (observational cognates of the SFE, dense gas fraction, and dense gas SFE) depend on the stellar surface density and the molecular/atomic ratio. The HCN/CO ratio correlates tightly with these two parameters across a range of 2.1 dex and increases in the high surface density parts of galaxies. Simultaneously, the IR/HCN ratio decreases systematically with these same parameters and is ~6-8 times lower near galaxy centers than in the outer regions. For fixed line-mass conversion factors, these results are incompatible with a simple model in which SF depends only on the gas mass above some density threshold. Only a specific set of environment-dependent conversion factors can render our observations compatible with such a model. Whole cloud models, such as the theory of turbulence regulated SF, do a better job of matching our data. We explore one such model in which variations in the Mach number and in the mean density would respectively drive the trends within galaxy disks and the differences between disk and merging galaxies (abridged).

preprint2014arXiv

A Bayesian blind survey for cold molecular gas in the Universe

A new Bayesian method for performing an image domain search for line-emitting galaxies is presented. The method uses both spatial and spectral information to robustly determine the source properties, employing either simple Gaussian, or other physically motivated models whilst using the evidence to determine the probability that the source is real. In this paper, we describe the method, and its application to both a simulated data set, and a blind survey for cold molecular gas using observations of the Hubble Deep Field North taken with the Plateau de Bure Interferometer. We make a total of 6 robust detections in the survey, 5 of which have counterparts in other observing bands. We identify the most secure detections found in a previous investigation, while finding one new probable line source with an optical ID not seen in the previous analysis. This study acts as a pilot application of Bayesian statistics to future searches to be carried out both for low-$J$ CO transitions of high redshift galaxies using the JVLA, and at millimeter wavelengths with ALMA, enabling the inference of robust scientific conclusions about the history of the molecular gas properties of star-forming galaxies in the Universe through cosmic time.

preprint2014arXiv

ALMA Imaging of HCN, CS and dust in Arp 220 and NGC 6240

We report ALMA Band 7 (350 GHz) imaging at 0.4 - 0.6arcsec resolution and Band 9 (696 GHz) at ~0.25arcsec resolution of the luminous IR galaxies Arp 220 and NGC 6240. The long wavelength dust continuum is used to estimate ISM masses for Arp 220 East, West and NGC 6240 of 1.9, 4.2 and 1.6x10^9 msun within radii of 69, 65 and 190 pc. The HCN emission was modeled to derive the emissivity distribution as a function of radius and the kinematics of each nuclear disk, yielding dynamical masses consistent with the masses and sizes derived from the dust emission. In Arp 220, the major dust and gas concentrations are at radii less than 50 pc in both counter-rotating nuclear disks. The thickness of the disks in Arp 220estimated from the velocity dispersion and rotation velocities are 10-20 pc and the mean gas densities are n_H2 ~10^5 cm^-3 at R < 50 pc. We develop an analytic treatment for the molecular excitation (including photon trapping), yielding volume densities for both the HCN and CS emission with n_H2 ~2x10^5 cm^-3. The agreement of the mean density from the total mass and size with that required for excitation suggests that the volume is essentially filled with dense gas, i.e. it is not cloudy or like swiss cheese.

preprint2014arXiv

ALMA Observation of 158 micron [CII] Line and Dust Continuum of a z=7 Normally Star-forming Galaxy in the Epoch of Reionization

We present ALMA observations of the [CII] line and far-infrared (FIR) continuum of a normally star-forming galaxy in the reionization epoch, the z=6.96 Ly-alpha emitter (LAE) IOK-1. Probing to sensitivities of sigma_line = 240 micro-Jy/beam (40 km/s channel) and sigma_cont = 21 micro-Jy/beam, we found the galaxy undetected in both [CII] and continuum. Comparison of UV - FIR spectral energy distribution (SED) of IOK-1, including our ALMA limit, with those of several types of local galaxies (including the effects of the cosmic microwave background, CMB, on the FIR continuum) suggests that IOK-1 is similar to local dwarf/irregular galaxies in SED shape rather than highly dusty/obscured galaxies. Moreover, our 3 sigma FIR continuum limit, corrected for CMB effects, implies intrinsic dust mass M_dust < 6.4 x 10^7 M_sun, FIR luminosity L_FIR < 3.7 x 10^{10} L_sun (42.5 - 122.5 micron), total IR luminosity L_IR < 5.7 x 10^{10} L_sun (8 - 1000 micron) and dust-obscured star formation rate (SFR) < 10 M_sun/yr, if we assume that IOK-1 has a dust temperature and emissivity index typical of local dwarf galaxies. This SFR is 2.4 times lower than one estimated from the UV continuum, suggesting that < 29% of the star formation is obscured by dust. Meanwhile, our 3 sigma [CII] flux limit translates into [CII] luminosity, L_[CII] < 3.4 x 10^7 L_sun. Locations of IOK-1 and previously observed LAEs on the L_[CII] vs. SFR and L_[CII]/L_FIR vs. L_FIR diagrams imply that LAEs in the reionization epoch have significantly lower gas and dust enrichment than AGN-powered systems and starbursts at similar/lower redshifts, as well as local star-forming galaxies.

preprint2014arXiv

ALMA Reveals the Molecular Medium Fueling the Nearest Nuclear Starburst

We use ALMA to derive the mass, length, and time scales associated with the nuclear starburst in NGC 253. This region forms ~2 M_sun/yr of stars and resembles other starbursts in scaling relations, with star formation consuming the gas reservoir 10 times faster than in galaxy disks. We present observations of CO, the high effective density transitions HCN(1-0), HCO+(1-0), CS(2-1), and their isotopologues. We identify ten clouds that appear as peaks in line emission and enhancements in the HCN-to-CO ratio. These clouds are massive (~10^7 M_sun) structures with sizes (~30 pc) similar to GMCs in other systems. Compared to disk galaxy GMCs, they show high line widths (~20-40 km/s) given their size, with implied Mach numbers ~90. The clouds also show high surface (~6,000 M_sun/pc^2) and volume densities (n_H2~2,000 cm^-3). Given these, self-gravity can explain the line widths. This short free fall time (~0.7 Myr) helps explain the more efficient star formation in NGC 253. We also consider the starburst region as a whole. The geometry is confused by the high inclination, but simple models support a non-axisymmetric, bar-like geometry with a compact, clumpy region of high gas density embedded in an extended CO distribution. Even for the whole region, the surface density still exceeds that of a disk galaxy GMC. The orbital time (~10 Myr), disk free fall time (<~ 3 Myr), and disk crossing time (<~ 3 Myr) are each much shorter than in a normal spiral galaxy disk. Some but not all aspects of the structure correspond to predictions from assuming vertical dynamical equilibrium or a marginally stable rotating disk. Finally, the CO-to-H2 conversion factor implied by our cloud calculations is approximately Galactic, contrasting with results showing a low value for the whole starburst region. The contrast provides resolved support for the idea of mixed molecular ISM phases in starburst galaxies.

preprint2014arXiv

An ALMA Survey of Submillimeter Galaxies in the Extended Chandra Deep Field South: The Redshift Distribution and Evolution of Submillimeter Galaxies

We present the first photometric redshift distribution for a large unbiased sample of 870um selected submillimeter galaxies (SMGs) with robust identifications based on observations with the Atacama Large Millimeter Array (ALMA). In our analysis we consider 96 SMGs in the Extended Chandra Deep Field South, 77 of which have 4-19 band, optical-near-infrared, photometry. We model the Spectral Energy Distributions (SEDs) for these 77 SMGs, deriving a median photometric redshift of z=2.3+/-0.1. The remaining 19 SMGs have insufficient optical or near-infrared photometry to derive photometric redshifts, but a stacking analysis of IRAC and Herschel observations confirms they are not spurious. Assuming these sources have an absolute H-band magnitude distribution comparable to that of a complete sample of z~1-2 SMGs, we demonstrate that the undetected SMGs lie at higher redshifts, raising the median redshift for SMGs to z=2.5+/-0.2. More critically we show that the proportion of galaxies undergoing an SMG phase at z>3 is 35+/-5% of the total population. We derive a median stellar mass for SMGs of Mstar=(8+/-1)x10^10Mo, but caution that there are significant systematic uncertainties in our stellar mass estimate, up to x5 for individual sources. We compare our sample of SMGs to a volume-limited, morphologically classified sample of ellipticals in the local Universe. Assuming the star formation activity in SMGs has a timescale of ~100Myr we show that their descendants at z~0 would have a space density and M_H distribution which are in good agreement with those of local ellipticals. In addition the inferred mass-weighted ages of the local ellipticals broadly agree with the look-back times of the SMG events. Taken together, these results are consistent with a simple model that identifies SMGs as events that form most of the stars seen in the majority of luminous elliptical galaxies at the present day.

preprint2014arXiv

An ALMA Survey of Submillimetre Galaxies in the Extended Chandra Deep Field South: The Far-Infrared Properties of SMGs

We exploit ALMA 870um (345GHz) observations of submillimetre sources in the Extended Chandra Deep Field South to investigate the far-infrared properties of high-redshift submillimetre galaxies (SMGs). Using the precisely located 870um ALMA positions of 99 SMGs, together with 24um and radio imaging of this field, we deblend the Herschel/SPIRE imaging of this region to extract their far-infrared fluxes and colours. The median photometric redshifts for ALMA LESS (ALESS) SMGs which are detected in at least two SPIRE bands increases with wavelength of the peak in their SEDs, with z=2.3+/-0.2, 2.5+/-0.3 and 3.5+/-0.5 for the 250, 350 and 500-um peakers respectively. We find that 34 ALESS SMGs do not have a >3-sigma counterpart at 250, 350 or 500-um. These galaxies have a median photometric redshift of z=3.3+/-0.5, which is higher than the full ALESS SMG sample; z=2.5+/-0.2. Using the photometric redshifts together with the 250-870um photometry, we estimate the far-infrared luminosities and characteristic dust temperature of each SMG. The median infrared luminosity of the S_870um>2mJy SMGs is L_IR=(3.0+/-0.3)x10^{12}Lo(SFR=300+/-30Mo/yr). At a fixed luminosity, the characteristic dust temperature of these high-redshift SMGs is 2-3K lower than comparably luminous galaxies at z=0, reflecting the more extended star formation occurring in these systems. By extrapolating the 870um number counts to S_ 870um=1mJy, we show that the contribution of S_870um>1mJy SMGs to the cosmic star formation budget is 20% of the total over the redshift range z~1-4. We derive a median dust mass for these SMGs of M_d=(3.6+/-0.3)x10^8Mo and by adopting an appropriate gas-to-dust ratio, we estimate an average molecular mass of M_H2=(4.2+/-0.4)x10^{10}Mo. Finally, we use our estimates of the H2 masses to show that SMGs with S_870um>1mJy contain ~10% of the z~2 volume-averaged H2 mass density at this epoch.

preprint2014arXiv

Andromeda's Dust

Spitzer Space Telescope and Herschel Space Observatory imaging of M31 is used, with a physical dust model, to construct maps of dust surface density, dust-to-gas ratio, starlight heating intensity, and PAH abundance, out to R=25kpc. The global dust mass is M_d=5.4x10^7Msol, the global dust/H mass ratio is M_d/M_H=0.0081, and the global PAH abundance is <q_PAH>=0.039. The dust surface density has an inner ring at R=5.6kpc, a maximum at R=11.2kpc, and an outer ring at R=15.1kpc. The dust/gas ratio varies from M_d/M_H=0.026 at the center to ~0.0027 at R=25kpc. From the dust/gas ratio, we estimate the ISM metallicity to vary by a factor ~10, from Z/Zsol=3 at R=0 to ~0.3 at R=25kpc. The dust heating rate parameter <U> peaks at the center, with <U> approx 35, declining to <U> approx 0.25 at R=20kpc. Within the central kpc, the starlight heating intensity inferred from the dust modeling is close to what is estimated from the stars in the bulge. The PAH abundance reaches a peak q_PAH=0.045 at R=11.2kpc. When allowance is made for the different spectrum of the bulge stars, q_PAH for the dust in the central kpc is similar to the overall value of q_PAH in the disk. The silicate-graphite-PAH dust model used here is generally able to reproduce the observed dust spectral energy distribution across M31, but overpredicts 500um emission at R=2-6kpc, suggesting that at R=2-6kpc, the dust opacity varies more steeply with frequency (with beta approx 2.3 between 200 and 600um) than in the model

preprint2014arXiv

Dust and gas in luminous proto-cluster galaxies at z=4.05: the case for different cosmic dust evolution in normal and starburst galaxies

We measure the dust and gas content of the three sub-millimeter galaxies (SMGs) in the GN20 proto-cluster at $z=4.05$ using new IRAM Plateau de Bure interferometer (PdBI) CO(4-3) and 1.2-3.3 mm continuum observations. All these three SMGs are heavily dust obscured, with UV-based star formation rate (SFR) estimates significantly smaller than the ones derived from the bolometric infrared (IR), consistent with the spatial offsets revealed by HST and CO imaging. Based also on evaluations of the specific SFR, CO-to-H$_2$ conversion factor and gas depletion timescale, we classify all the three galaxies as starbursts (SBs), although with a lower confidence for GN20.2b that might be a later stage merging event. We place our measurements in the context of the evolutionary properties of main sequence (MS) and SB galaxies. ULIRGs have 3-5 times larger $L'_{\rm CO}/M_{\rm dust}$ and $M_{\rm dust}/M_\star$ ratios than $z=0$ MS galaxies, but by $z\sim2$ the difference appears to be blurred, probably due to differential metallicity evolution. SB galaxies appear to slowly evolve in their $L'_{\rm CO}/M_{\rm dust}$ and $M_{\rm dust}/M_\star$ ratios all the way to $z>6$ (consistent with rapid enrichment of SB events), while MS galaxies rapidly increase in $M_{\rm dust}/M_\star$ from $z=0$ to 2 (due to gas fraction increase, compensated by a decrease of metallicities). While no IR/submm continuum detection is available for indisputably normal massive galaxies at $z>2.5$, we show that if metallicity indeed decrease rapidly for these systems at $z>3$ as claimed in the literature, we should expect a strong decrease of their $M_{\rm dust}/M_\star$, consistent with recent PdBI and ALMA upper limits. We conclude that the $M_{\rm dust}/M_\star$ ratio could be a powerful tool for distinguishing starbursts from normal galaxies at $z>4$.

preprint2014arXiv

Dust continuum emission as a tracer of gas mass in galaxies

We use a sample of 36 galaxies from the KINGFISH (Herschel IR), HERACLES (IRAM CO), and THINGS (VLA HI) surveys to study empirical relations between Herschel infrared (IR) luminosities and the total mass of the interstellar gas (H2+HI). Such a comparison provides a simple empirical relationship without introducing the uncertainty of dust model fitting. We find tight correlations, and provide fits to these relations, between Herschel luminosities and the total gas mass integrated over entire galaxies, with the tightest, almost linear, correlation found for the longest wavelength data (SPIRE500). However, we find that accounting for the gas-phase metallicity (affecting the dust-to-gas ratio) is crucial when applying these relations to low-mass, and presumably high-redshift, galaxies. The molecular (H2) gas mass is found to be better correlated with the peak of the IR emission (e.g. PACS160), driven mostly by the correlation of stellar mass and mean dust temperature. When examining these relations as a function of galactocentric radius we find the same correlations, albeit with a larger scatter, up to a radius of 0.7 r_25 (within which most of the galaxy's baryonic mass resides). However, beyond this radius the same correlations no longer hold, with the gas mass (predominantly HI) increasing relative to the infrared emission. The tight relations found for the bulk of the galaxy's baryonic content suggest that the total gas masses of disk-like (non-merging) galaxies can be inferred from far-infrared continuum measurements in situations where only the latter are available, e.g. in ALMA continuum observations of high-redshift galaxies.

preprint2014arXiv

High-Resolution Radio Continuum Measurements of the Nuclear Disks of Arp 220

We present new Karl G. Jansky Very Large Array radio continuum images of the nuclei of Arp 220, the nearest ultra-luminous infrared galaxy. These images have both the angular resolution to study detailed morphologies of the two nuclei that power the system and sensitivity to a wide range of spatial scales. At 33 GHz, and with a resolution of 0".081 x 0".063 (29.9 x 23.3 pc), we resolve the emission surrounding both nuclei and conclude that is mostly synchrotron in nature. The spatial distributions of radio emission in both nuclei are well described by exponential profiles. These have deconvolved half-light radii of 51 and 35 pc for the eastern and western nuclei, and they match the number density profile of radio supernovae observed with very long baseline interferometry. This similarity might be due to the fast cooling of cosmic rays electrons caused by the presence of a strong (~ mG) magnetic field in this system. We estimate high luminosity surface densities of $\mathrm{Σ_{IR} \sim 4.2^{+1.6}_{-0.7} \times 10^{13}}$ (east) and $\mathrm{\sim 9.7^{+3.7}_{-2.4} \times 10^{13}~(west)~L_{\odot}~kpc^{-2}}$, and star formation rate surface densities of $\mathrm{Σ_{SFR} \sim 10^{3.7\pm0.1}}$ (east) and $\mathrm{\sim 10^{4.1\pm0.1}~(west)~M_{\odot}~yr^{-1}~kpc^{-2}}$. These values, especially for the western nucleus are, to our knowledge, the highest luminosity and star formation rate surface densities measured for any star-forming system. Despite these high values, the nuclei lie below the dusty Eddington limit in which radiation pressure is balanced only by self-gravity. The small measured sizes also imply that the nuclei of Arp 220 are only transparent in the frequency range ~ 5 to 350 GHz. Our results offer no clear evidence that an active galactic nucleus dominates the emission from either nucleus at 33 GHz.

preprint2014arXiv

Pinpointing the Molecular Gas within a Lyman Alpha Blob at z ~ 2.7

We present IRAM PdBI observations of the CO(3-2) and CO(5-4) line transitions from a Ly-alpha blob at z~2.7 in order to investigate the gas kinematics, determine the location of the dominant energy source, and study the physical conditions of the molecular gas. CO line and dust continuum emission are detected at the location of a strong MIPS source that is offset by ~1.5" from the Ly-alpha peak. Neither of these emission components is resolved with the 1.7" beam, showing that the gas and dust are confined to within ~7kpc from this galaxy. No millimeter source is found at the location of the Ly-alpha peak, ruling out a central compact source of star formation as the power source for the Ly-alpha emission. Combined with a spatially-resolved spectrum of Ly-alpha and HeII, we constrain the kinematics of the extended gas using the CO emission as a tracer of the systemic redshift. Near the MIPS source, the Ly-alpha profile is symmetric and its line center agrees with that of CO line, implying that there are no significant bulk flows and that the photo-ionization from the MIPS source might be the dominant source of the Ly-alpha emission. In the region near the Ly-alpha peak, the gas is slowly receding (~100km/s) with respect to the MIPS source, thus making the hyper-/superwind hypothesis unlikely. We find a sub-thermal line ratio between two CO transitions, I_CO(5-4)/I_CO(3-2)=0.97+/-0.21. This line ratio is lower than the average values found in high-z SMGs and QSOs, but consistent with the value found in the Galactic center, suggesting that there is a large reservoir of low-density molecular gas that is spread over the MIPS source and its vicinity.

preprint2014arXiv

Polycyclic Aromatic Hydrocarbon and Mid-Infrared Continuum Emission in a z>4 Submillimeter Galaxy

We report the detection of 6.2um polycyclic aromatic hydrocarbon (PAH) and rest-frame 4-7um continuum emission in the z=4.055 submillimeter galaxy GN20, using the Infrared Spectrograph (IRS) on board the Spitzer Space Telescope. This represents the first detection of PAH emission at z>4. The strength of the PAH emission feature is consistent with a very high star formation rate of ~1600Msun/yr. We find that this intense starburst powers at least ~1/3 of the faint underlying 6um continuum emission, with an additional, significant (and perhaps dominant) contribution due to a power-law-like hot dust source, which we interpret to likely be a faint, dust-obscured active galactic nucleus (AGN). The inferred 6um AGN continuum luminosity is consistent with a sensitive upper limit on the hard X-ray emission as measured by the Chandra X-Ray Observatory if the previously undetected AGN is Compton-thick. This is in agreement with the finding at optical/infrared wavelengths that the galaxy and its nucleus are heavily dust-obscured. Despite the strong power-law component enhancing the mid-infrared continuum emission, the intense starburst associated with the photon-dominated regions that give rise to the PAH emission appears to dominate the total energy output in the infrared. GN20 is one of the most luminous starburst galaxies known at any redshift, embedded in a rich protocluster of star-forming galaxies. This investigation provides an improved understanding of the energy sources that power such exceptional systems, which represent the extreme end of massive galaxy formation at early cosmic times.

preprint2014arXiv

Search for [CII] emission in z=6.5-11 star-forming galaxies

We present the search for the [CII] emission line in three $z>6.5$ Lyman-alpha emitters (LAEs) and one J-Dropout galaxy using the Combined Array for Research in Millimeter-wave Astronomy (CARMA) and the Plateau de Bure Interferometer (PdBI). We observed three bright $z\sim6.5-7$ LAEs discovered in the SUBARU deep field (SDF) and the Multiple Imaged lensed $z\sim 11$ galaxy candidate found behind the galaxy cluster MACSJ0647.7+7015. For the LAEs IOK-1 ($z=6.965$), SDF J132415.7+273058 ($z=6.541$) and SDF J132408.3+271543 ($z=6.554$) we find upper limits for the [CII] line luminosity of $<2.05$, $<4.52$ and $<10.56\times10^{8}{\rm L}_{\odot}$ respectively. We find upper limits to the FIR luminosity of the galaxies using a spectral energy distribution template of the local galaxy NGC 6946 and taking into account the effects of the Cosmic Microwave Background on the mm observations. For IOK-1, SDF J132415.7+273058 and SDF J132408.3+271543 we find upper limits for the FIR luminosity of $<2.33$, $3.79$ and $7.72\times10^{11}{\rm L}_{\odot}$ respectively. For the lensed galaxy MACS0647-JD, one of the highest redshift galaxy candidate to date with $z_{\rm ph}=10.7^{+0.6}_{-0.4}$ we put an upper limit in the [CII] emission of $<1.36\times10^{8}\times(μ/15)^{-1}{\rm L}_{\odot}$ and an upper limit in the FIR luminosity of $<6.1\times10^{10}\times(μ/15)^{-1}{\rm L}_{\odot}$ (where $μ$ is the magnification factor). We explore the different conditions relevant for the search for [CII] emission in high redshift galaxies as well as the difficulties for future observations with ALMA and CCAT.

preprint2014arXiv

The GISMO 2-millimeter Deep Field in GOODS-N

We present deep continuum observations using the GISMO camera at a wavelength of 2 mm centered on the Hubble Deep Field (HDF) in the GOODS-N field. These are the first deep field observations ever obtained at this wavelength. The 1 sigma sensitivity in the innermost approx. 4 arcminutes of the 7 utes map is approx. 135 uJy/beam, a factor of three higher in flux/beam sensitivity than the deepest available SCUBA 850 um observations, and almost a factor of four higher in flux / beam sensitivity than the combined MAMBO/AzTEC 1.2 mm observations of this region. Our source extraction algorithm identifies 12 sources directly, and another 3 through correlation with known sources at 1.2 mm and 850 um. Five of the directly detected GISMO sources have counterparts in the MAMBO/AzTEC catalog, and four of those also have SCUBA counterparts. HDF850.1, one of the first blank-field detected submillimeter galaxies, is now detected at 2 mm. The median redshift of all sources with counterparts of known redshifts is med(z) = 2.91 +/- 0.94. Statistically, the detections are most likely real for 5 of the seven 2 mm sources without shorter wavelength counterparts, while the probability for none of them being real is negligible.

preprint2013arXiv

A deep search for molecular gas in two massive Lyman break galaxies at z=3 and 4: vanishing CO-emission due to low metallicity?

We present deep IRAM Plateau de Bure Interferometer (PdBI) observations, searching for CO-emission toward two massive, non-lensed Lyman break galaxies (LBGs) at z=3.216 and 4.058. With one low significance CO detection (3.5 sigma) and one sensitive upper limit, we find that the CO lines are >~ 3-4 times weaker than expected based on the relation between IR and CO luminosities followed by similarly, massive galaxies at z=0-2.5. This is consistent with a scenario in which these galaxies have low metallicity, causing an increased CO-to-H_2 conversion factor, i.e., weaker CO-emission for a given molecular (H_2) mass. The required metallicities at z>3 are lower than predicted by the fundamental metallicity relation (FMR) at these redshifts, consistent with independent evidence. Unless our galaxies are atypical in this respect, detecting molecular gas in normal galaxies at z>3 may thus remain challenging even with ALMA.

preprint2013arXiv

A high-dispersion molecular gas component in nearby galaxies

We present a comprehensive study of the velocity dispersion of the atomic (HI) and molecular (H2) gas components in the disks (R < R25) of a sample of 12 nearby spiral galaxies with moderate inclinations. Our analysis is based on sensitive high resolution data from the THINGS (atomic gas) and HERACLES (molecular gas) surveys. To obtain reliable measurements of the velocity dispersion, we stack regions several kilo-parsecs in size, after accounting for intrinsic velocity shifts due to galactic rotation and large-scale motions. We stack using various parameters: the galacto-centric distance, star formation rate surface density, HI surface density, H2 surface density, and total gas surface density. We fit single Gaussian components to the stacked spectra and measure median velocity dispersions for HI of 11.9 +/- 3.1 km/s and for H2 of 12.0 +/- 3.9 km/s. The CO velocity dispersions are thus, surprisingly, very similar to the corresponding ones of HI, with an average ratio of sigma(HI)/sigma(CO) = 1.0 +/- 0.2 irrespective of the stacking parameter. The measured CO velocity dispersions are significantly higher (factor 2) than the traditional picture of a cold molecular gas disk associated with star formation. The high dispersion implies an additional thick molecular gas disk (possibly as thick as the HI disk). Our finding is in agreement with recent sensitive measurements in individual edge-on and face-on galaxies and points towards the general existence of a thick disk of molecular gas, in addition to the well-known thin disk in nearby spiral galaxies.

preprint2013arXiv

An ALMA survey of submillimetre galaxies in the Extended Chandra Deep Field South: High resolution 870um source counts

We report the first counts of faint submillimetre galaxies (SMG) in the 870-um band derived from arcsecond resolution observations with the Atacama Large Millimeter Array (ALMA). We have used ALMA to map a sample of 122 870-um-selected submillimetre sources drawn from the (0.5x0.5)deg^2 LABOCA Extended Chandra Deep Field South Submillimetre Survey (LESS). These ALMA maps have an average depth of sigma(870um)~0.4mJy, some ~3x deeper than the original LABOCA survey and critically the angular resolution is more than an order of magnitude higher, FWHM of ~1.5" compared to ~19" for the LABOCA discovery map. This combination of sensitivity and resolution allows us to precisely pin-point the SMGs contributing to the submillimetre sources from the LABOCA map, free from the effects of confusion. We show that our ALMA-derived SMG counts broadly agree with the submillimetre source counts from previous, lower-resolution single-dish surveys, demonstrating that the bulk of the submillimetre sources are not caused by blending of unresolved SMGs. The difficulty which well-constrained theoretical models have in reproducing the high-surface densities of SMGs, thus remains. However, our observations do show that all of the very brightest sources in the LESS sample, S(870um)>12mJy, comprise emission from multiple, fainter SMGs, each with 870-um fluxes of <9mJy. This implies a natural limit to the star-formation rate in SMGs of <10^3 M_Sun/yr, which in turn suggests that the space densities of z>1 galaxies with gas masses in excess of ~5x10^10 M_Sun is <10^-5 Mpc^-3. We also discuss the influence of this blending on the identification and characterisation of the SMG counterparts to these bright submillimetre sources and suggest that it may be responsible for previous claims that they lie at higher redshifts than fainter SMGs.

preprint2013arXiv

CARMA Survey Toward Infrared-bright Nearby Galaxies (STING). III. The Dependence of Atomic and Molecular Gas Surface Densities on Galaxy Properties

We investigate the correlation between CO and HI emission in 18 nearby galaxies from the CARMA Survey Toward IR-Bright Nearby Galaxies (STING) at sub-kpc and kpc scales. Our sample, spanning a wide range in stellar mass and metallicity, reveals evidence for a metallicity dependence of the HI column density measured in regions exhibiting CO emission. Such a dependence is predicted by the equilibrium model of McKee & Krumholz, which balances H_2 formation and dissociation. The observed HI column density is often smaller than predicted by the model, an effect we attribute to unresolved clumping, although values close to the model prediction are also seen. We do not observe HI column densities much larger than predicted, as might be expected were there a diffuse HI component that did not contribute to H_2 shielding. We also find that the H_2 column density inferred from CO correlates strongly with the stellar surface density, suggesting that the local supply of molecular gas is tightly regulated by the stellar disk.

preprint2013arXiv

Clumping and the Interpretation of kpc-Scale Maps of the Interstellar Medium: Smooth HI and Clumpy, Variable H2 Surface Density

Many recent models consider the structure of individual interstellar medium (ISM) clouds as a way to explain observations of large parts of galaxies. To compare such models to observations, one must understand how to translate between surface densities observed averaging over large (~kpc) scales and surface densities on the scale of individual clouds (~pc scale), which are treated by models. We define a "clumping factor" that captures this translation as the ratio of the mass-weighted surface density, which is often the quantity of physical interest, to the area-weighted surface density, which is observed. We use high spatial resolution (sub-kpc) maps of CO and HI emission from nearby galaxies to measure the clumping factor of both atomic and molecular gas. The molecular and atomic ISM exhibit dramatically different degrees of clumping. As a result, the ratio H2/HI measured at ~kpc resolution cannot be trivially interpreted as a cloud-scale ratio of surface densities. HI emission appears very smooth, with a clumping factor of only ~1.3. Based on the scarce and heterogeneous high resolution data available, CO emission is far more clumped with a widely variable clumping factor, median ~7 for our heterogeneous data. Our measurements do not provide evidence for a universal mass-weighted surface density of molecular gas, but also cannot conclusively rule out such a scenario. We suggest that a more sophisticated treatment of molecular ISM structure, one informed by high spatial resolution CO maps, is needed to link cloud-scale models to kpc-scale observations of galaxies.

preprint2013arXiv

Cool Gas in High Redshift Galaxies

Over the last decade, observations of the cool interstellar medium in distant galaxies via molecular and atomic fine structure line emission has gone from a curious look into a few extreme, rare objects, to a mainstream tool to study galaxy formation, out to the highest redshifts. Molecular gas has now been observed in close to 200 galaxies at z>1, including numerous AGN host-galaxies out to z~7, highly starforming sub-millimeter galaxies (median redshift z~2.5), and increasing samples of 'main-sequence' star forming galaxies at z~1.5-2.5. Studies have moved well beyond simple detections, to dynamical imaging at kpc-scale resolution, and multi-line, multi-species studies that determine the physical conditions in the interstellar medium. Observations of the cool gas are the required complement to studies of the stellar density and star formation history of the Universe, as they reveal the phase of the interstellar medium that immediately precedes star formation. Current observations suggest that the order of magnitude increase in the cosmic star formation rate density from z~0 to 2 is commensurate with a similar increase in the gas to stellar mass ratio in star forming disk galaxies. Progress has been made on determining the CO luminosity to H_2 mass conversion factor at high-z, and the dichotomy between high versus low depletion time values for main sequence versus starburst galaxies, respectively, with a likely dependence on metallicity and other local physical conditions. Studies of atomic fine structure line emission are rapidly progressing, with some tens of galaxies detected in the exceptionally bright [CII] 158 micron line to date. This line is proving to be a unique tracer of galaxy dynamics in the early Universe and has the potential to be the most direct means of obtaining spectroscopic redshifts for the first galaxies during cosmic reionization.

preprint2013arXiv

Dust May Be More Rare Than Expected in Metal Poor Galaxies

'Normal' galaxies observed at z>6, when the Universe was <1 billion years old, thus far show no evidence of the cold dust that accompanies star formation in the local Universe, where the dust-to-gas mass ratio is 1%. A prototypical example is 'Himiko' (z=6.6), which a mere 840 Myr after the Big Bang is forming stars at a rate of 30-100 Msun/yr, yielding a mass assembly time M^{star}/SFR 150x10^6 yr. Himiko is estimated to have a low fraction (2-3% of the Solar value) of elements heavier than helium (metallicity), and although its gas mass cannot be asserted at this time its dust-to-stellar mass ratio is constrained to be <0.05%. The local galaxy I Zw 18, with a metallicity 4% solar and forming stars less rapidly than Himiko but still vigorously for its mass (M^{star}/SFR 1.6x10^9 yr), is also very dust deficient and perhaps one of the best analogues of primitive galaxies accessible to detailed study. Here we report observations of dust emission from I Zw 18 from which we determine its dust mass to be 450-1800 Msun, yielding a dust-to-stellar mass ratio \approx 10^{-6}-10^{-5} and a dust-to-gas mass ratio 3.2-13x10^{-6}. If I Zw 18 is a reasonable analog of Himiko, then Himiko's dust mass is \approx 50,000 Msun, a factor of 100 below the current upper limit. These numbers are considerably uncertain, but if most high-z galaxies are more like Himiko than like the quasar host SDSS J114816.64+525150.3, then the prospects for detecting the gas and dust in them are much poorer than hitherto anticipated.

preprint2013arXiv

Empirical predictions for (sub-)millimeter line and continuum deep fields

[abridged] Modern (sub-)millimeter/radio interferometers will enable us to measure the dust and molecular gas emission from galaxies that have luminosities lower than the Milky Way, out to high redshifts and with unprecedented spatial resolution and sensitivity. This will provide new constraints on the star formation properties and gas reservoir in galaxies throughout cosmic times through dedicated deep field campaigns targeting the CO/[CII] lines and dust continuum emission. In this paper, we present empirical predictions for such (sub-)millimeter line and continuum deep fields. We base these predictions on the deepest available optical/near-infrared ACS and NICMOS data on the Hubble Ultra Deep Field. Using a physically-motivated spectral energy distribution model, we fit the observed optical/near-infrared emission of 13,099 galaxies with redshifts up to z=5, and obtain median likelihood estimates of their stellar mass, star formation rate, dust attenuation and dust luminosity. We derive statistical constraints on the dust emission in the infrared and (sub-)millimeter which are consistent with the observed optical/near-infrared emission in terms of energy balance. This allows us to estimate, for each galaxy, the (sub-)millimeter continuum flux densities in several ALMA, PdBI/NOEMA and JVLA bands. Using empirical relations between the observed CO/[CII] line luminosities and the infrared luminosity, we infer the flux of the CO(1-0) and [CII] lines from the estimated infrared luminosity of each galaxy in our sample. We then predict the fluxes of higher CO transition lines CO(2-1) to CO(7-6) bracketing two extreme gas excitation scenarios. We use our predictions to discuss possible deep field strategies with ALMA. The predictions presented in this study will serve as a direct benchmark for future deep field campaigns in the (sub-)millimeter regime.

preprint2013arXiv

Investigating the presence of 500 um submillimeter excess emission in local star forming galaxies

Submillimeter excess emission has been reported at 500 microns in a handful of local galaxies, and previous studies suggest that it could be correlated with metal abundance. We investigate the presence of an excess submillimeter emission at 500 microns for a sample of 20 galaxies from the Key Insights on Nearby Galaxies: a Far Infrared Survey with Herschel (KINGFISH) that span a range of morphologies and metallicities (12+log(O/H)=7.8-8.7). We probe the far-infrared (IR) emission using images from the Spitzer Space Telescope and Herschel Space Observatory in the wavelength range 24-500 microns. We model the far-IR peak of the dust emission with a two-temperature modified blackbody and measure excess of the 500 micron photometry relative to that predicted by our model. We compare the submillimeter excess, where present, with global galaxy metallicity and, where available, resolved metallicity measurements. We do not find any correlation between the 500 micron excess and metallicity. A few individual sources do show excess (10-20%) at 500 microns; conversely, for other sources, the model overpredicts the measured 500 micron flux density by as much as 20%, creating a 500 micron "deficit". None of our sources has an excess larger than the calculated 1-sigma uncertainty, leading us to conclude that there is no substantial excess at submillimeter wavelengths at or shorter than 500 microns in our sample. Our results differ from previous studies detecting 500 micron excess in KINGFISH galaxies largely due to new, improved photometry used in this study.

preprint2013arXiv

Molecular Gas and Star Formation in Nearby Disk Galaxies

We compare molecular gas traced by 12CO(2-1) maps from the HERACLES survey, with tracers of the recent star formation rate (SFR) across 30 nearby disk galaxies. We demonstrate a first-order linear correspondence between Sig_mol and Sig_SFR but also find important second-order systematic variations in the apparent molecular gas depletion time, t_dep^mol = Sig_mol / Sig_SFR. At our 1 kpc common resolution, CO correlates closely with many tracers of the recent SFR. Weighting each line of sight equally and using a fixed, Milky Way alpha_CO, our data yield a molecular gas depletion time, t_dep^mol=Sig_mol/Sig_SFR ~ 2.2 Gyr with 0.3 dex scatter, in good agreement with literature data. We apply a forward-modeling approach to constrain the power-law index, N, that relates the SFR surface density and the molecular gas surface density and find N=1+/-0.15 for our full data set with some variation from galaxy to galaxy. However, we caution that a power law treatment oversimplifies the topic given that we observe correlations between t_dep^mol and other local and global quantities. The strongest of these are a decreased t_dep^mol in low-mass, low-metallicity galaxies and a correlation of the kpc-scale t_dep^mol with dust-to-gas ratio, D/G. These correlations can be explained by a CO-to-H2 conversion factor that depends on D/G in the theoretically expected way. This is not a unique interpretation, but external evidence of conversion factor variations makes it a conservative one. After applying a D/G-dependent alpha_CO, some weak correlations between t_dep^mol and local conditions persist. In particular, we observe lower t_dep^mol and enhanced CO excitation associated with some nuclear gas concentrations. These appear to reflect real enhancements in the SFR/H2 and t_dep appears multivalued at fixed Sig_mol, supporting the the idea of "disk" and "starburst" modes driven by environmental factors.

preprint2013arXiv

Nuclear Water Maser Emission in Centaurus A

We report the detection of a 22GHz water maser line in the nearest (D~3.8Mpc) radio galaxy Centaurus A using the Australia Telescope Compact Array (ATCA). The line is centered at a velocity of ~960kms-1, which is redshifted by about 400kms-1 from the systemic velocity. Such an offset, as well as the width of ~120kms-1, could be consistent with either a nuclear maser arising from an accretion disk of the central supermassive black hole, or for a jet maser that is emitted from the material that is shocked near the base of the jet in Centaurus\,A. The best spatial resolution of our ATCA data constrains the origin of the maser feature within <3pc from the supermassive black hole. The maser exhibits a luminosity of ~1Lo, which classifies it as a kilomaser, and appears to be variable on timescales of months. A kilomaser can also be emitted by shocked gas in star forming regions. Given the small projected distance from the core, the large offset from systemic velocity, as well as the smoothness of the line feature, we conclude that a jet maser line emitted by shocked gas around the base of the AGN is the most likely explanation. For this scenario we can infer a minimum density of the radio jet of ~>10cm-3, which indicates substantial mass entrainment of surrounding gas into the propagating jet material.

preprint2013arXiv

On the effect of the Cosmic Microwave Background in high-redshift (sub-)millimeter observations

Modern (sub-)millimeter interferometers enable the measurement of the cool gas and dust emission of high-redshift galaxies (z>5). However, at these redshifts the cosmic microwave background (CMB) temperature is higher, approaching, and even exceeding, the temperature of cold dust and molecular gas observed in the local Universe. In this paper, we discuss the impact of the warmer CMB on (sub-)millimeter observations of high-redshift galaxies. The CMB affects the observed (sub-)millimeter dust continuum and the line emission (e.g. carbon monoxide, CO) in two ways: (i) it provides an additional source of (both dust and gas) heating; and (ii) it is a non-negligible background against which the line and continuum emission are measured. We show that these two competing processes affect the way we interpret the dust and gas properties of high-redshift galaxies using spectral energy distribution models. We quantify these effects and provide correction factors to compute what fraction of the intrinsic dust (and line) emission can be detected against the CMB as a function of frequency, redshift and temperature. We discuss implications on the derived properties of high-redshift galaxies from (sub-)millimeter data. Specifically, the inferred dust and molecular gas masses can be severely underestimated for cold systems if the impact of the CMB is not properly taken into account.

preprint2013arXiv

Star Formation and Gas Kinematics of Quasar Host Galaxies at z~6: New insights from ALMA

We present ALMA observations of the [C II] 158 micron fine structure line and dust continuum emission from the host galaxies of five redshift 6 quasars. We also report complementary observations of 250 GHz dust continuum and CO (6-5) line emission from the z=6.00 quasar SDSS J231038.88+185519.7. The ALMA observations were carried out in the extended array at 0.7" resolution. We have detected the line and dust continuum in all five objects. The derived [C II] line luminosities are 1.6x10^{9} to 8.8x10^{9} Lsun and the [C II]-to-FIR luminosity ratios are 3.0-5.6x10^{-4}, which is comparable to the values found in other high-redshift quasar-starburst systems and local ultra-luminous infrared galaxies. The sources are marginally resolved and the intrinsic source sizes (major axis FWHM) are constrained to be 0.3" to 0.6" (i.e., 1.7 to 3.5 kpc) for the [C II] line emission and 0.2" to 0.4" (i.e., 1.2 to 2.3 kpc) for the continuum. These measurements indicate that there is vigorous star formation over the central few kpc in the quasar host galaxies. The ALMA observations also constrain the dynamical properties of the atomic gas in the starburst nuclei. The intensity-weighted velocity maps of three sources show clear velocity gradients. Such velocity gradients are consistent with a rotating, gravitationally bound gas component, although they are not uniquely interpreted as such. Under the simplifying assumption of rotation, the implied dynamical masses within the [C II]-emitting regions are of order 10^{10} to 10^{11} Msun. Given these estimates, the mass ratios between the SMBHs and the spheroidal bulge are an order of magnitude higher than the mean value found in local spheroidal galaxies, which is in agreement with results from previous CO observations of high redshift quasars.

preprint2013arXiv

The galaxy environment of a QSO at z ~ 5.7

High-redshift quasars are believed to reside in massive halos in the early universe and should therefore be located in fields with overdensities of galaxies, which are thought to evolve into galaxy clusters seen in the local universe. However, despite many efforts, the relationship between galaxy overdensities and z~6 quasars is ambiguous. This can possibly be attributed to the difficulty of finding galaxies with accurate redshifts in the vicinity of z~6 quasars. So far, overdensity searches around z~6 quasars have been based on studies of Lyman break galaxies (LBGs), which probe a redshift range of Delta z ~ 1. This range is large enough to select galaxies that may not be physically related to the quasar. We use deep narrow- and broadband imaging to study the environment of the z=5.72 quasar ULAS J0203+0012. The redshift range probed by our narrow-band selection of Lyman alpha emitters (LAEs) is Delta z ~ 0.1, which is significantly narrower than the LBG searches. This is the first time that LAEs were searched for near a z~6 quasar, in an effort to provide clues about the environments of quasars at the end of the epoch of reionization. We find no enhancement of LAEs in the surroundings of ULAS J0203+0012 in comparison with blank fields. We explore different explanations and interpretations for this non-detection of a galaxy overdensity, including that (1) the strong ionization from the quasar may prevent galaxy formation in its immediate vicinity and (2) high-redshift quasars may not reside in the center of the most massive dark matter halos.

preprint2013arXiv

The Green Bank Telescope Maps the Dense, Star-Forming Gas in the Nearby Starburst Galaxy M82

Observations of the Milky Way and nearby galaxies show that dense molecular gas correlates with recent star formation, suggesting that the formation of this gas phase may help regulate star formation. A key test of this idea requires wide-area, high-resolution maps of dense molecular gas in galaxies to explore how local physical conditions drive dense gas formation, but these observations have been limited because of the faintness of dense gas tracers like HCN and HCO+. Here we demonstrate the power of the Robert C. Byrd Green Bank Telescope -- the largest single-dish millimeter radio telescope -- for mapping dense gas in galaxies by presenting the most sensitive maps yet of HCN and HCO+ in the starburst galaxy M82. The HCN and HCO+ in the disk of this galaxy correlates with both recent star formation and more diffuse molecular gas and shows kinematics consistent with a rotating torus. The HCO+ emission extending to the north and south of the disk is coincident with the outflow previously identified in CO and traces the eastern edge of the hot outflowing gas. The central starburst region has a higher ratio of star formation to dense gas than the outer regions, pointing to the starburst as a key driver of this relationship. These results establish that the GBT can efficiently map the dense molecular gas at 90 GHz in nearby galaxies, a capability that will increase further with the 16 element feed array under construction.

preprint2013arXiv

The LBT Bootes Field Survey: I. The Rest-frame UV and Luminosity Functions and Clustering of Bright Lyman Break Galaxies at z~3

We present a deep LBT/LBC U-band imaging survey (9 deg2) covering the NOAO Bootes field. A total of 14,485 Lyman Break Galaxies (LBGs) at z~3 are selected, which are used to measure the rest-frame UV luminosity function (LF). The large sample size and survey area reduce the LF uncertainties due to Poisson statistics and cosmic variance by >3 compared to previous studies. At the bright end, the LF shows excess power compared to the best-fit Schechter function, which can be attributed to the contribution of $z\sim3$ quasars. We compute the rest-frame near-infrared LF and stellar mass function (SMF) of z~3 LBGs based on the R-band and IRAC [4.5 micro m]-band flux relation. We investigate the evolution of the UV LFs and SMFs between z~7 and z~3, which supports a rising star formation history in the LBGs. We study the spatial correlation function of two bright LBG samples and estimate their average host halo mass. We find a tight relation between the host halo mass and the galaxy star formation rate (SFR),which follows the trend predicted by the baryonic accretion rate onto the halo, suggesting that the star formation in LBGs is fueled by baryonic accretion through the cosmic web. By comparing the SFRs with the total baryonic accretion rates, we find that cosmic star formation efficiency is about 5%-20% and it does not evolve significantly with redshift, halo mass, or galaxy luminosity.

preprint2013arXiv

The Starburst-Driven Molecular Wind in NGC 253 and the Suppression of Star Formation

The under-abundance of very massive galaxies in the universe is frequently attributed to the effect of galactic winds. Although ionized galactic winds are readily observable most of the expelled mass is likely in cooler atomic and molecular phases. Expanding molecular shells observed in starburst systems such as NGC 253 and M 82 may facilitate the entrainment of molecular gas in the wind. While shell properties are well constrained, determining the amount of outflowing gas emerging from such shells and the connection between this gas and the ionized wind requires spatial resolution <100 pc coupled with sensitivity to a wide range of spatial scales, hitherto not available. Here we report observations of NGC 253, a nearby starburst galaxy (D~3.4 Mpc) known to possess a wind, which trace the cool molecular wind at 50 pc resolution. At this resolution the extraplanar molecular gas closely tracks the Hα filaments, and it appears connected to molecular expanding shells located in the starburst region. These observations allow us to directly measure the molecular outflow rate to be > 3 Msun/yr and likely ~9 Msun/yr. This implies a ratio of mass-outflow rate to star formation rate of at least η~1-3, establishing the importance of the starburst-driven wind in limiting the star formation activity and the final stellar content.

preprint2012arXiv

A UV Ultra-luminous Lyman Break Galaxy at Z=2.78 in NDWFS Bootes Field

We present one of the most ultraviolet (UV) luminous Lyman Break Galaxies (LBGs) (J1432+3358) at z=2.78, discovered in the NOAO Deep Wide-Field Survey (NDWFS) Bootes field. The R-band magnitude of J1432+3358 is 22.29 AB, more than two magnitudes brighter than typical L* LBGs at this redshift. The deep z-band image reveals two components of J1432+3358 separated by 1.0" with flux ratio of 3:1. The high signal-to-noise ratio (S/N) rest-frame UV spectrum shows Lya emission line and interstellar medium absorption lines. The absence of NV and CIV emission lines, the non-detection in X-ray and radio wavelengths and mid-infrared (MIR) colors indicate no or weak active galactic nuclei (AGN) (<10%) in this galaxy. The galaxy shows broader line profile with the full width half maximum (FWHM) of about 1000 km/s and larger outflow velocity (~500 km/s) than those of typical z~3 LBGs. The physical properties are derived by fitting the spectral energy distribution (SED) with stellar synthesis models. The dust extinction, E(B-V)=0.12, is similar to that in normal LBGs. The star formation rates (SFRs) derived from the SED fitting and the dust-corrected UV flux are consistent with each other, ~300 Msun/yr, and the stellar mass is 1.3e11 Msun. The SFR and stellar mass in J1432+3358 are about an order of magnitude higher than those in normal LBGs. The SED-fitting results support that J1432+3358 has a continuous star formation history with the star formation episode of 630 Myr. The morphology of J1432+3358 and its physical properties suggest that J1432+3358 is in an early phase of 3:1 merger process. The unique properties and the low space number density (~1e-7 Mpc^{-3})are consistent with the interpretation that such galaxies are either found in a short unobscured phase of the star formation or that small fraction of intensive star-forming galaxies are unobscured.

preprint2012arXiv

An ALMA survey of Sub-millimetre Galaxies in the Extended Chandra Deep Field South: Detection of [C II] at z=4.4

We present ALMA 870-um (345GHz) observations of two sub-millimetre galaxies (SMGs) drawn from an ALMA study of the 126 sub-millimeter sources from the LABOCA Extended Chandra Deep Field South Survey (LESS). The ALMA data identify the counterparts to these previously unidentified sub-millimeter sources and serendipitously detect bright emission lines in their spectra which we show are most likely to be [C II]157.74um emission yielding redshifts of z=4.42 and z=4.44. This blind detection rate within the 7.5-GHz bandpass of ALMA is consistent with the previously derived photometric redshift distribution of SMGs and suggests a modest, but not dominant (<25%), tail of 870-um selected SMGs at z>4. We find that the ratio of L_CII/L_FIR in these SMGs is much higher than seen for similarly far-infrared-luminous galaxies at z~0, which is attributed to the more extended gas reservoirs in these high-redshift ULIRGs. Indeed, in one system we show that the [C II] emission shows hints of extended emission on >3kpc scales. Finally, we use the volume probed by our ALMA survey to show that the bright end of the [C II] luminosity function evolves strongly between z=0 and z~4.4, reflecting the increased ISM cooling in galaxies as a result of their higher star-formation rates. These observations demonstrate that even with short integrations, ALMA is able to detect the dominant fine structure cooling lines from high-redshift ULIRGs, measure their energetics and trace their evolution with redshift.

preprint2012arXiv

Dust-to-Gas Ratio in the Extremely Metal Poor Galaxy I ZW 18

The blue compact dwarf galaxy I Zw 18 is one of the most metal poor systems known in the local Universe (12 + log(O/H) $=$ 7.17). In this work we study I Zw 18 using data from {\it Spitzer}, {\it Herschel Space Telescope} and IRAM Plateau de Bure Interferometer. Our data set includes the most sensitive maps of I Zw 18, to date, in both, the far infrared and the CO $J=1\rightarrow0$ transition. We use dust emission models to derive a dust mass upper limit of only M$_{dust}\leq1.1\times10^4$ M$_{\odot}$ ($3σ$ limit). This upper limit is driven by the non-detection at 160 $μ$m, and it is a factor of 4-10 times smaller than previous estimates (depending upon the model used). We also estimate an upper limit to the total dust-to-gas mass ratio of M$_{Dust}$/M$_{gas}\leq5.0\times10^{-5}$. If a linear correlation between the dust-to-gas mass ratio and metallicity (measure as O/H) were to hold, we would expect a ratio of 3.9$\times10^{-4}$. We also show that the infrared SED is similar to that of starbursting systems.

preprint2012arXiv

Estimating the Star Formation Rate at 1 kpc Scales in Nearby Galaxies

Using combinations of Hα, ultraviolet (UV), and infrared (IR) emission, we estimate the star formation rate (SFR) surface density, Σ_SFR, at 1 kpc resolution for 30 disk galaxies that are targets of the IRAM HERACLES CO survey. We present a new physically-motivated IR spectral energy distribution-based approach to account for possible contributions to 24\mum emission not associated with recent star formation. Considering a variety of "reference" SFRs from the literature, we revisit the calibration of the 24\mum term in hybrid (UV+IR or Hα+IR) tracers. We show that the overall calibration of this term remains uncertain at the factor of two level because of the lack of wide-field, robust reference SFR estimates. Within this uncertainty, published calibrations represent a reasonable starting point for 1 kpc-wide areas of star-forming disk galaxies but we re-derive and refine the calibration of the IR term in these tracers to match our resolution and approach to 24\mum emission. We compare a large suite of Σ_SFR estimates and find that above Σ_SFR \sim 10^-3 M_\odot yr^-1 kpc^-2 the systematic differences among tracers are less than a factor of two across two orders of magnitude dynamic range. We caution that methodology and data both become serious issues below this level. We note from simple model considerations that focusing on a part of a galaxy dominated by a single stellar population the intrinsic uncertainty in Hαand FUV-based SFRs are \sim 0.3 and \sim 0.5 dex.

preprint2012arXiv

Herschel-PACS observations of [OI]63um towards submillimetre galaxies at z~1

We present Herschel-PACS spectroscopy of the [OI]63um far-infrared cooling line from a sample of six unlensed and spectroscopically-confirmed 870um-selected submillimetre (submm) galaxies (SMGs) at 1.1<z<1.6 from the LABOCA Extended Chandra Deep Field South (ECDFS) Submm Survey (LESS). This is the first survey of [OI]63um, one of the main photodissociation region (PDR) cooling lines, in SMGs. New high-resolution ALMA interferometric 870um continuum imaging confirms that these six Herschel-targeted SMG counterparts are bona fide sources of submm emission. We detect [OI]63um in two SMGs with a SNR >3, tentatively detect [OI]63um in one SMG, and constrain the line flux for the non-detections. We also exploit the combination of submm continuum photometry from 250-870um and our new PACS continuum measurements to constrain the far-infrared (FIR) luminosity, L_FIR, in these SMGs to < 30%. We find that SMGs do not show a deficit in their [OI]63um-to-far-infrared continuum luminosity ratios (with ratios ranging from ~0.5-1.5%), similar to what was seen previously for the [CII]158um-to-FIR ratios in SMGs. These observed ratios are about an order of magnitude higher than what is seen typically for local ultra luminous infrared galaxies (ULIRGs), which adds to the growing body of evidence that SMGs are not simply `scaled up' versions of local ULIRGs. Rather, the PDR line-to-L_FIR ratios suggest that the star formation modes of SMGs are likely more akin to that of local normal (lower-luminosity) star-forming galaxies, with the bulk of the star formation occurring in extended regions, galaxy-scale (~kpc) in size. These observations represent the first step towards a census of the major PDR cooling lines in typical SMGs that will be attainable with ALMA, enabling detailed modelling to probe the global properties of the star formation and the evolutionary status of SMGs.

preprint2012arXiv

HST narrow-band search for extended Ly-alpha emission around two z>6 quasars

We search for extended Ly-alpha emission around two z>6 quasars, SDSS J1030+0524 (z=6.309) and SDSS J1148+5251 (z=6.419) using WFC3 narrow-band filters on board the Hubble Space Telescope. For each quasar, we collected two deep, narrow-band images, one sampling the Ly-alpha line+continuum at the quasar redshifts and one of the continuum emission redwards of the line. After carefully modeling the Point Spread Function, we find no evidence for extended Ly-alpha emission. These observations set 2-sigma limits of L(Ly-alpha, extended) < 3.2 x 10^{44} erg/s for J1030+0524 and L(Ly-alpha, extended) < 2.5 x 10^{44} erg/s for J1148+5251. Given the star formation rates typically inferred from (rest-frame) far-infrared measurements of z~6 quasars, these limits are well below the intrinsic bright Ly-alpha emission expected from the recombination of gas photoionized by the quasars or by the star formation in the host galaxies, and point towards significant Ly-alpha suppression or dust attenuation. However, small extinction values have been observed along the line of sight to the nuclei, thus reddening has to be coupled with other mechanisms for Ly-alpha suppression (e.g., resonance scattering). No Ly-alpha emitting companions are found, down to a 5-sigma sensitivity of ~ 1 x 10^{-17} erg/s/cm^2/arcsec^2 (surface brightness) and ~ 5 x 10^{-17} erg/s/cm^2 (assuming point sources).

preprint2012arXiv

Little Things

We present LITTLE THINGS (Local Irregulars That Trace Luminosity Extremes, The HI Nearby Galaxy Survey) that is aimed at determining what drives star formation in dwarf galaxies. This is a multi-wavelength survey of 37 Dwarf Irregular and 4 Blue Compact Dwarf galaxies that is centered around HI-line data obtained with the National Radio Astronomy Observatory (NRAO) Very Large Array (VLA). The HI-line data are characterized by high sensitivity (less than 1.1 mJy/beam per channel), high spectral resolution (less than or equal to 2.6 km/s), and high angular resolution (~6 arcseconds. The LITTLE THINGS sample contains dwarf galaxies that are relatively nearby (less than or equal to 10.3 Mpc; 6 arcseconds is less than or equal to 300 pc), that were known to contain atomic hydrogen, the fuel for star formation, and that cover a large range in dwarf galactic properties. We describe our VLA data acquisition, calibration, and mapping procedures, as well as HI map characteristics, and show channel maps, moment maps, velocity-flux profiles, and surface gas density profiles. In addition to the HI data we have GALEX UV and ground-based UBV and Halpha images for most of the galaxies, and JHK images for some. Spitzer mid-IR images are available for many of the galaxies as well. These data sets are available on-line.

preprint2012arXiv

Low CO Luminosities in Dwarf Galaxies

[Abridged] We present maps of CO 2-1 emission covering the entire star-forming disks of 16 nearby dwarf galaxies observed by the IRAM HERACLES survey. The data have 13 arcsec angular resolution, ~250 pc at our average distance of 4 Mpc, and sample the galaxies by 10-1000 resolution elements. We apply stacking techniques to perform the first sensitive search for CO emission in dwarfs outside the Local Group ranging from single lines-of-sight, stacked over IR-bright regions of embedded star formation, and stacked over the entire galaxy. We detect 5 dwarfs in CO with total luminosities of L_CO = 3-28 1e6 Kkmspc2. The other 11 dwarfs remain undetected in CO even in the stacked data and have L_CO < 0.4-8 1e6 Kkmspc2. We combine our sample of dwarfs with a large literature sample of spirals to study scaling relations of L_CO with M_B and metallicity. We find that dwarfs with metallicities of Z ~ 1/2-1/10 Z_sun have L_CO about 1e2-1e4x smaller than spirals and that their L_CO per unit L_B is 10-100x smaller. A comparison with tracers of star formation (FUV and 24 micron) shows that L_CO per unit SFR is 10-100x smaller in dwarfs. One possible interpretation is that dwarfs form stars much more efficiently, however we argue that the low L_CO/SFR ratio is due to significant changes of the CO-to-H2 conversion factor, alpha_CO, in low metallicity environments. Assuming a constant H2 depletion time of 1.8 Gyr (as found for nearby spirals) implies alpha_CO values for dwarfs with Z ~ 1/2-1/10 Z_sun that are more than 10x higher than those found in solar metallicity spirals. This significant increase of alpha_CO at low metallicity is consistent with previous studies, in particular those which model dust emission to constrain H2 masses. Even though it is difficult to parameterize the metallicity dependence of alpha_CO, our results suggest that CO is increasingly difficult to detect at lower metallicities.

preprint2012arXiv

On the Origin of the Supergiant HI Shell and Putative Companion in NGC 6822

We present new Hubble Space Telescope Advanced Camera for Surveys imaging of six positions spanning 5.8 kpc of the HI major axis of the Local Group dIrr NGC 6822, including both the putative companion galaxy and the large HI hole. The resulting deep color magnitude diagrams show that NGC 6822 has formed >50% of its stars in the last ~5 Gyr. The star formation histories of all six positions are similar over the most recent 500 Myr, including low-level star formation throughout this interval and a weak increase in star formation rate during the most recent 50 Myr. Stellar feedback can create the giant HI hole, assuming that the lifetime of the structure is longer than 500 Myr; such long-lived structures have now been observed in multiple systems and may be the norm in galaxies with solid-body rotation. The old stellar populations (red giants and red clump stars) of the putative companion are consistent with those of the extended halo of NGC 6822; this argues against the interpretation of this structure as a bona fide interacting companion galaxy and against its being linked to the formation of the HI hole via an interaction. Since there is no evidence in the stellar population of a companion galaxy, the most likely explanation of the extended HI structure in NGC 6822 is a warped disk inclined to the line of sight.

preprint2012arXiv

Star Formation in Quasar Host Galaxies at Redshift 6: Millimeter Surveys and New Insights from ALMA

We have been carrying out a systematic survey of the star formation and ISM properties in the host galaxies of z~6 quasars. Our 250 GHz observations, together with available data from the literature, yield a sample of 14 z~6 quasars that are bright in millimeter dust continuum emission with estimated FIR luminosities of a few 10^12 to 10^13 Lsun. Most of these millimeter-detected z~6 quasars have also been detected in molecular CO line emission, indicating molecular gas masses on order of 10^10 Msun. We have searched for [C II] 158 micron fine structure line emission toward four of the millimeter bright z~6 quasars with ALMA and all of them have been detected. All these results suggest massive star formation at rates of about 600 to 2000 Msun/yr over the central few kpc region of these quasar host galaxies.

preprint2012arXiv

The heating of dust by old stellar populations in the Bulge of M31

We use new Herschel multi-band imaging of the Andromeda galaxy to analyze how dust heating occurs in the central regions of galaxy spheroids that are essentially devoid of young stars. We construct a dust temperature map of M31 through fitting modified blackbody SEDs to the Herschel data, and find that the temperature within 2 kpc rises strongly from the mean value in the disk of 17 pm 1K to \sim35K at the centre. UV to near-IR imaging of the central few kpc shows directly the absence of young stellar populations, delineates the radial profile of the stellar density, and demonstrates that even the near-UV dust extinction is optically thin in M31's bulge. This allows the direct calculation of the stellar radiation heating in the bulge, U\ast(r), as a function of radius. The increasing temperature profile in the centre matches that expected from the stellar heating, i.e. that the dust heating and cooling rates track each other over nearly two orders of magnitude in U\ast. The modelled dust heating is in excess of the observed dust temperatures, suggesting that it is more than sufficient to explain the observed IR emission. Together with the wavelength dependent absorption cross section of the dust, this demonstrates directly that it is the optical, not UV, radiation that sets the heating rate. This analysis shows that neither young stellar populations nor stellar near-UV radiation are necessary to heat dust to warm temperatures in galaxy spheroids. Rather, it is the high densities of Gyr-old stellar populations that provide a sufficiently strong diffuse radiation field to heat the dust. To the extent which these results pertain to the tenuous dust found in the centres of early-type galaxies remains yet to be explored.

preprint2012arXiv

The Intense Starburst HDF850.1 in a Galaxy Overdensity at z=5.2 in the Hubble Deep Field

The Hubble Deep Field (HDF) is a region in the sky that provides one of the deepest multi-wavelength views of the distant universe and has led to the detection of thousands of galaxies seen throughout cosmic time. An early map of the HDF at a wavelength of 850 microns that is sensitive to dust emission powered by star formation revealed the brightest source in the field, dubbed HDF850.1. For more than a decade, this source remained elusive and, despite significant efforts, no counterpart at shorter wavelengths, and thus no redshift, size or mass, could be identified. Here we report, using a millimeter wave molecular line scan, an unambiguous redshift determination for HDF850.1 of z=5.183. This places HDF850.1 in a galaxy overdensity at z~5.2 in the HDF, corresponding to a cosmic age of only 1.1 Gyr after the Big Bang. This redshift is significantly higher than earlier estimates and higher than most of the >100 sub-millimeter bright galaxies identified to date. The source has a star formation rate of 850 M_sun/yr and is spatially resolved on scales of 5 kpc, with an implied dynamical mass of ~1.3x10^11 M_sun, a significant fraction of which is present in the form of molecular gas. Despite our accurate redshift and position, a counterpart arising from starlight remains elusive.

preprint2012arXiv

The Shapes of the HI Velocity Profiles of the THINGS Galaxies

We analyze the shapes of the HI velocity profiles of The HI Nearby Galaxy Survey (THINGS) to study the phase structure of the neutral interstellar medium (ISM) and its relation to global galaxy properties. We use a method analogous to the stacking method sometimes used in high redshift HI observations to construct high signal-to-noise (S/N) profiles. We call these high S/N profiles super profiles. We analyze and discuss possible systematics that may change the observed shapes of the super profiles. After quantifying these effects and selecting a sub-sample of unaffected galaxies, we find that the super profiles are best described by a narrow and a broad Gaussian component, which are evidence of the presence of the Cold Neutral Medium (CNM) and the Warm Neutral Medium (WNM). The velocity dispersion of the narrow component range from ~3.4 to ~8.6 km/s with an average of 6.5+/-1.5 km/s, whereas that of the broad component range from ~10.1 to ~24.3 km/s with an average of 16.8+/-4.3 km/s. We find that the super profile parameters correlate with star formation indicators such as metallicity, FUV-NUV colors and H_alpha luminosities. The flux ratio between the narrow and broad components tends to be highest for high metallicity, high star formation rate (SFR) galaxies. We show that the narrow component identified in the super profiles is associated with the presence of star formation, and possibly with molecular hydrogen.

preprint2012arXiv

Tracing Cold HI Gas in Nearby, Low-Mass Galaxies

We analyze line-of-sight atomic hydrogen (HI) line profiles of 31 nearby, low-mass galaxies selected from the Very Large Array - ACS Nearby Galaxy Survey Treasury (VLA-ANGST) and The HI Nearby Galaxy Survey (THINGS) to trace regions containing cold (T $\lesssim$ 1400 K) HI from observations with a uniform linear scale of 200 pc/beam. Our galaxy sample spans four orders of magnitude in total HI mass and nine magnitudes in M_B. We fit single and multiple component functions to each spectrum to isolate the cold, neutral medium given by a low dispersion (<6 km/s) component of the spectrum. Most HI spectra are adequately fit by a single Gaussian with a dispersion of 8-12 km/s. Cold HI is found in 23 of 27 (~85%) galaxies after a reduction of the sample size due to quality control cuts. The cold HI contributes ~20% of the total line-of-sight flux when found with warm HI. Spectra best fit by a single Gaussian, but dominated by cold HI emission (i.e., have velocity dispersions <6 km/s) are found primarily beyond the optical radius of the host galaxy. The cold HI is typically found in localized regions and is generally not coincident with the very highest surface density peaks of the global HI distribution (which are usually areas of recent star formation). We find a lower limit for the mass fraction of cold-to-total HI gas of only a few percent in each galaxy.

preprint2012arXiv

VLA-ANGST: A high-resolution HI Survey of Nearby Dwarf Galaxies

We present the "Very Large Array survey of Advanced Camera for Surveys Nearby Galaxy Survey Treasury galaxies (VLA-ANGST)." VLA-ANGST is a National Radio Astronomy Observatory Large Program consisting of high spectral (0.6-2.6 km/s) and spatial (~6") resolution observations of neutral, atomic hydrogen (HI) emission toward 35 nearby dwarf galaxies from the ANGST survey. ANGST is a systematic HST survey to establish a legacy of uniform multi-color photometry of resolved stars for a volume-limited sample of nearby galaxies (D\lesssim4 Mpc). VLA-ANGST provides VLA HI observations of the sub-sample of ANGST galaxies with recent star formation that are observable from the northern hemisphere and that were not observed in the "The HI Nearby Galaxy Survey" (THINGS). The overarching scientific goal of VLA-ANGST is to investigate fundamental characteristics of the neutral interstellar medium (ISM) of dwarf galaxies. Here we describe the VLA observations, the data reduction, and the final VLA-ANGST data products. We present an atlas of the integrated HI maps, the intensity-weighted velocity fields, the second moment maps as a measure for the velocity dispersion of the HI, individual channel maps, and integrated HI spectra for each VLA-ANGST galaxy. We closely follow the observational setup and data reduction of THINGS to achieve comparable sensitivity and angular resolution. A major difference, however, is the high velocity resolution of the VLA-ANGST observations (0.65 and 1.3km/s for the majority of the galaxies). The VLA-ANGST data products are made publicly available at: https://science.nrao.edu/science/surveys/vla-angst. With available star formation histories from resolved stellar populations and lower resolution ancillary observations from the FIR to the UV, VLA-ANGST will enable detailed studies of the relationship between the ISM and star formation in dwarf galaxies on a ~100 pc scale.

preprint2011arXiv

A Molecular Star Formation Law in the Atomic Gas Dominated Regime in Nearby Galaxies

We use the IRAM HERACLES survey to study CO emission from 33 nearby spiral galaxies down to very low intensities. Using atomic hydrogen (HI) data, mostly from THINGS, we predict the local mean CO velocity from the mean HI velocity. By renormalizing the CO velocity axis so that zero corresponds to the local mean HI velocity we are able to stack spectra coherently over large regions as function of radius. This enables us to measure CO intensities with high significance as low as Ico = 0.3 K km/s (H2_SD = 1 Msun/pc2), an improvement of about one order of magnitude over previous studies. We detect CO out to radii Rgal = R25 and find the CO radial profile to follow a uniform exponential decline with scale length of 0.2 R25. Comparing our sensitive CO profiles to matched profiles of HI, Halpha, FUV, and IR emission at 24um and 70um, we observe a tight, roughly linear relation between CO and IR intensity that does not show any notable break between regions that are dominated by molecular (H2) gas (H2_SD > HI_SD) and those dominated by atomic gas (H2_SD < HI_SD). We use combinations of FUV+24um and Halpha+24um to estimate the recent star formation rate (SFR) surface density, SFR_SD, and find approximately linear relations between SFR_SD and H2_SD. We interpret this as evidence for stars forming in molecular gas with little dependence on the local total gas surface density. While galaxies display small internal variations in the SFR-to-H2 ratio, we do observe systematic galaxy-to-galaxy variations. These galaxy-to-galaxy variations dominate the scatter in relations between CO and SFR tracers measured at large scales. The variations have the sense that less massive galaxies exhibit larger ratios of SFR-to-CO than massive galaxies. Unlike the SFR-to-CO ratio, the balance between HI and H2 depends strongly on the total gas surface density and radius. It must also depend on additional parameters.

preprint2011arXiv

A Survey of Atomic Carbon at High Redshift

We present a survey of atomic carbon (CI) emission in high-redshift (z>2) submillimeter galaxies (SMGs) and quasar host galaxies (QSOs). Sensitive observations of the CI(3P_1->3P_0) and CI(3P_2->3P_1) lines have been obtained at the IRAM Plateau de Bure interferometer and the IRAM 30m telescope. A total of 16 CI lines have been targeted in 10 sources, leading to a total of 10 detected lines --- this doubles the number of CI observations at high redshift to date. We include previously published CI observations (an additional 5 detected sources) in our analysis. Our main finding is that the CI properties of the studied high-redshift systems do not differ significantly from what is found in low-redshift systems, including the Milky Way. The CI(3P_2->3P_1)/CI(3P_1->3P_0) and the CI(3P_1->3P_0)/12CO(3-2) line luminosity (L') ratios change little in our sample, with respective ratios of 0.55+/-0.15 and 0.32+/-0.13. The CI lines are not an important contributor to cooling of the molecular gas (average L_CI/L_FIR ~ (7.7+/-4.6) x 10^-6). We derive a mean carbon excitation temperature of 29.1+/-6.3 K, broadly consistent with dust temperatures derived for high-redshift starforming systems, but lower than gas temperatures typically derived for starbursts in the local universe. The carbon abundance of X_CI/X_H2~8.4+/-3.5 x 10^-5 is of the same order as found in the Milky Way and nearby galaxies. This implies that the high-z galaxies studied here are significantly enriched in carbon on galactic scales, even though the look-back times are considerable (the average redshift of the sample sources corresponds to an age of the universe of ~2 Gyr).

preprint2011arXiv

CARMA Survey Toward Infrared-bright Nearby Galaxies (STING) II: Molecular Gas Star Formation Law and Depletion Time Across the Blue Sequence

We present an analysis of the relationship between molecular gas and current star formation rate surface density at sub-kpc and kpc scales in a sample of 14 nearby star-forming galaxies. Measuring the relationship in the bright, high molecular gas surface density ($\Shtwo\gtrsim$20 \msunpc) regions of the disks to minimize the contribution from diffuse extended emission, we find an approximately linear relation between molecular gas and star formation rate surface density, $\nmol\sim0.96\pm0.16$, with a molecular gas depletion time $\tdep\sim2.30\pm1.32$ Gyr. We show that, in the molecular regions of our galaxies there are no clear correlations between \tdep\ and the free-fall and effective Jeans dynamical times throughout the sample. We do not find strong trends in the power-law index of the spatially resolved molecular gas star formation law or the molecular gas depletion time across the range of galactic stellar masses sampled (\mstar $\sim$$10^{9.7}-10^{11.5}$ \msun). There is a trend, however, in global measurements that is particularly marked for low mass galaxies. We suggest this trend is probably due to the low surface brightness CO, and it is likely associated with changes in CO-to-H2 conversion factor.

preprint2011arXiv

CARMA Survey Toward Infrared-bright Nearby Galaxies (STING): Molecular Gas Star Formation Law in NGC4254

This study explores the effects of different assumptions and systematics on the determination of the local, spatially resolved star formation law. Using four star formation rate (SFR) tracers (Hαwith azimuthally averaged extinction correction, mid-infrared 24 micron, combined Hαand mid-infrared 24 micron, and combined far-ultraviolet and mid-infrared 24 micron), several fitting procedures, and different sampling strategies we probe the relation between SFR and molecular gas at various spatial resolutions and surface densities within the central 6.5 kpc in the disk of NGC4254. We find that in the high surface brightness regions of NGC4254 the form of the molecular gas star formation law is robustly determined and approximately linear and independent of the assumed fraction of diffuse emission and the SFR tracer employed. When the low surface brightness regions are included, the slope of the star formation law depends primarily on the assumed fraction of diffuse emission. In such case, results range from linear when the fraction of diffuse emission in the SFR tracer is ~30% or less (or when diffuse emission is removed in both the star formation and the molecular gas tracer), to super-linear when the diffuse fraction is ~50% and above. We find that the tightness of the correlation between gas and star formation varies with the choice of star formation tracer. The 24 micron SFR tracer by itself shows the tightest correlation with the molecular gas surface density, whereas the Hαcorrected for extinction using an azimuthally-averaged correction shows the highest dispersion. We find that for R<0.5R_25 the local star formation efficiency is constant and similar to that observed in other large spirals, with a molecular gas depletion time ~2 Gyr.

preprint2011arXiv

CO (2-1) Line Emission in Redshift 6 Quasar Host Galaxies

We report new observations of CO (2-1) line emission toward five z~6 quasars using the Ka-band receiver system on the Expanded Very Large Array (EVLA). Strong detections were obtained in two of them, SDSS J092721.82+200123.7 and CFHQS J142952.17+544717.6, and a marginal detection was obtained in another source, SDSS J084035.09+562419.9. Upper limits of the CO (2-1) line emission have been obtained for the other two objects. The CO (2-1) line detection in J0927+2001, together with previous measurements of the CO (6-5) and (5-4) lines, reveals important constraints on the CO excitation in the central ~10 kpc region of the quasar host galaxy. The CO (2-1) line emission from J1429+5447 is resolved into two distinct peaks separated by 1.2" (~6.9 kpc), indicating a possible gas-rich, major merging system, and the optical quasar position is consistent with the west peak. This result is in good agreement with the picture in which intense host galaxy star formation is coeval with rapid supermassive black hole accretion in the most distant universe. The two EVLA detections are ideal targets for further high-resolution imaging (e.g., with ALMA or EVLA observations) to study the gas distribution, dynamics, and SMBH-bulge mass relation in these earliest quasar-host galaxy systems.

preprint2011arXiv

CO(1-0) in z>2 Quasar Host Galaxies: No Evidence for Extended Molecular Gas Reservoirs

We report the detection of CO(1-0) emission in the strongly lensed high-redshift quasars IRAS F10214+4724 (z=2.286), the Cloverleaf (z=2.558), RX J0911+0551 (z=2.796), SMM J04135+10277 (z=2.846), and MG 0751+2716 (z=3.200), using the Expanded Very Large Array and the Green Bank Telescope. We report lensing-corrected CO(1-0) line luminosities of L'(CO) = 0.34-18.4 x 10^10 K km/s pc^2 and total molecular gas masses of M(H2) = 0.27-14.7 x 10^10 Msun for the sources in our sample. Based on CO line ratios relative to previously reported observations in J>=3 rotational transitions and line excitation modeling, we find that the CO(1-0) line strengths in our targets are consistent with single, highly-excited gas components with constant brightness temperature up to mid-J levels. We thus do not find any evidence for luminous extended, low excitation, low surface brightness molecular gas components. These properties are comparable to those found in z>4 quasars with existing CO(1-0) observations. These findings stand in contrast to recent CO(1-0) observations of z~2-4 submillimeter galaxies (SMGs), which have lower CO excitation and show evidence for multiple excitation components, including some low-excitation gas. These findings are consistent with the picture that gas-rich quasars and SMGs represent different stages in the early evolution of massive galaxies.

preprint2011arXiv

Complex Radio Spectral Energy Distributions in Luminous and Ultraluminous Infrared Galaxies

We use the Expanded Very Large Array to image radio continuum emission from local luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs) in 1 GHz windows centered at 4.7, 7.2, 29, and 36 GHz. This allows us to probe the integrated radio spectral energy distribution (SED) of the most energetic galaxies in the local universe. The 4-8 GHz flux densities agree well with previous measurements. They yield spectral indices α\approx -0.67 (where F_ν\propto ν^α) with \pm 0.15 (1σ) scatter, typical of nonthermal (synchrotron) emission from star-forming galaxies. The contrast of our 4-8 GHz data with literature 1.5 and 8.4 GHz flux densities gives further evidence for curvature of the radio SED of U/LIRGs. The SED appears flatter near \sim 1 GHz than near \sim 6 GHz, suggesting significant optical depth effects at the lower frequencies. The high frequency (28-37 GHz) flux densities are low compared to extrapolations from the 4-8 GHz data. We confirm and extend to higher frequency a previously observed deficit of high frequency radio emission for luminous starburst galaxies.

preprint2011arXiv

Constraining Dust and Molecular Gas Properties in Lyman Alpha Blobs at z~3

In order to constrain the bolometric luminosities, dust properties and molecular gas content of giant Lyman alpha nebulae, the so-called Lyman alpha blobs, we have carried out a study of dust continuum and CO line emission in two well-studied representatives of this population at z ~ 3: a Lya blob discovered by its strong Spitzer MIPS 24um detection (LABd05; Dey et al. 2005) and the Steidel blob 1 (SSA22-LAB01; Steidel et al. 2000). We find that the spectral energy distribution of LABd05 is well described by an AGN-starburst composite template with L(FIR) = (4.0 +/- 0.5) x 10^12 Lsun, comparable to high-z sub-millimeter galaxies and ultraluminous infrared galaxies. New APEX/LABOCA 870um measurements rule out the reported SCUBA detection of the SSA22-LAB01 (S[850um] = 16.8 mJy) at the > 4sigma level. Consistent with this, ultra-deep Plateau de Bure Interferometer (PdBI) observations with ~2arcsec spatial resolution also fail to detect any 1.2mm continuum source down to ~0.45mJy per beam (3sigma). Combined with the existing (sub)mm observations in the literature, we conclude that the FIR luminosity of SSA22-LAB01 remains uncertain. No CO line is detected in either case down to integrated flux limits of (Snu dV) < 0.25--1.0 Jy km/s, indicating a modest molecular gas reservoir, M(H_2) < 1--3 x 10^10 Msun. The non-detections exclude, with high significance (12 sigma), the previous tentative detection of a CO(4-3) line in the SSA22-LAB01. The increased sensitivity afforded by ALMA will be critical in studying molecular gas and dust in these interesting systems.

preprint2011arXiv

Dark and luminous matter in THINGS dwarf galaxies

We present mass models for the dark matter component of seven dwarf galaxies taken from "The HI Nearby Galaxy Survey" (THINGS) and compare these with those from numerical Lambda Cold Dark Matter (LCDM) simulations. The THINGS high-resolution data significantly reduce observational uncertainties and thus allow us to derive accurate dark matter distributions in these systems. We here use the bulk velocity fields when deriving the rotation curves of the galaxies. Compared to other types of velocity fields, the bulk velocity field minimizes the effect of small-scale random motions more effectively and traces the underlying kinematics of a galaxy more properly. The "Spitzer Infrared Nearby Galaxies Survey" (SINGS) 3.6 micron and ancillary optical data are used for separating the baryons from their total matter content in the galaxies. The sample dwarf galaxies are found to be dark matter dominated over most radii. We find discrepancies between the derived dark matter distributions of the galaxies and those of LCDM simulations, even after corrections for non-circular motions have been applied. The observed solid body-like rotation curves of the galaxies rise too slowly to reflect the cusp-like dark matter distribution in CDM halos. Instead, they are better described by core-like models such as pseudo-isothermal halo models dominated by a central constant-density core. The mean value of the logarithmic inner slopes of the mass density profiles is alpha = -0.29 +- 0.07. They are significantly different from the steep slope of ~ -1.0 inferred from previous dark-matter-only simulations, and are more consistent with shallower slopes found in recent LCDM simulations of dwarf galaxies in which the effects of baryonic feedback processes are included.

preprint2011arXiv

Evidence for non-evolving FeII/MgII ratios in rapidly accreting z~6 QSOs

[abridged] Quasars (QSOs) at the highest known redshift (z~6) are unique probes of the early growth of supermassive black holes (BHs). Until now, only the most luminous QSOs have been studied, often one object at a time. Here we present the most extensive consistent analysis to date of z>4 QSOs with observed NIR spectra, combining three new z~6 objects from our ongoing VLT-ISAAC program with nineteen 4<z<6.5 sources from the literature. The new sources extend the existing SDSS sample towards the faint end of the QSO luminosity function. Using a maximum likelihood fitting routine optimized for our spectral decomposition, we estimate the black hole mass (MBH), the Eddington ratio (defined as Lbol/LEdd) and the FeII/MgII line ratio, a proxy for the chemical abundance, to characterize both the central object and the broad line region gas. The QSOs in our sample host BHs with masses of ~10^9 M\odot that are accreting close to the Eddington luminosity, consistent with earlier results. We find that the distribution of observed Eddington ratios is significantly different than that of a luminosity-matched comparison sample of SDSS QSOs at lower redshift (0.35<z<2.25): the average <log(Lbol/LEdd)>=-0.37 (Lbol/LEdd~0.43) with a scatter of 0.20 dex for the z>4 sample and the <log(Lbol/LEdd)>=-0.80 (Lbol/LEdd~0.16) with a scatter of 0.24 dex for the 0.35<z<2.25 sample. This implies that, at a given luminosity, the MBH at high-z is typically lower than the average MBH of the lower-redshift population, i.e. the z>4 sources are accreting significantly faster than the lower-redshift ones. We show that the derived FeII/MgII ratios depend sensitively on the performed analysis: our self-consistent, homogeneous analysis significantly reduces the FeII/MgII scatter found in previous studies. The measured FeII/MgII line ratios show no sign of evolution with cosmic time in the redshift range 4<z<6.5 [...]

preprint2011arXiv

Extended Cold Molecular Gas Reservoirs in z~3.4 Submillimeter Galaxies

We report the detection of spatially resolved CO(1-0) emission in the z~3.4 submillimeter galaxies (SMGs) SMM J09431+4700 and SMM J13120+4242, using the Expanded Very Large Array (EVLA). SMM J09431+4700 is resolved into the two previously reported millimeter sources H6 and H7, separated by ~30kpc in projection. We derive CO(1-0) line luminosities of L'(CO 1-0) = (2.49+/-0.86) and (5.82+/-1.22) x 10^10 K km/s pc^2 for H6 and H7, and L'(CO 1-0) = (23.4+/-4.1) x 10^10 K km/s pc^2 for SMM J13120+4242. These are ~1.5-4.5x higher than what is expected from simple excitation modeling of higher-J CO lines, suggesting the presence of copious amounts of low-excitation gas. This is supported by the finding that the CO(1-0) line in SMM J13120+4242, the system with lowest CO excitation, appears to have a broader profile and more extended spatial structure than seen in higher-J CO lines (which is less prominently seen in SMM J09431+4700). Based on L'(CO 1-0) and excitation modeling, we find M_gas = 2.0-4.3 and 4.7-12.7 x 10^10 Msun for H6 and H7, and M_gas = 18.7-69.4 x 10^10 Msun for SMM J13120+4242. The observed CO(1-0) properties are consistent with the picture that SMM J09431+4700 represents an early-stage, gas-rich major merger, and that SMM J13120+4242 represents such a system in an advanced stage. This study thus highlights the importance of spatially and dynamically resolved CO(1-0) observations of SMGs to further understand the gas physics that drive star formation in these distant galaxies, which becomes possible only now that the EVLA rises to its full capabilities.

preprint2011arXiv

Far-infrared and Molecular CO Emission From the Host Galaxies of Faint Quasars at z~6

We present new millimeter and radio observations of nine z~6 quasars discovered in deep optical and near-infrared surveys. We observed the 250 GHz continuum in eight of the nine objects and detected three of them. New 1.4 GHz radio continuum data have been obtained for four sources, and one has been detected. We searched for molecular CO (6-5) line emission in the three 250 GHz detections and detected two of them. We study the FIR and radio emission and quasar-host galaxy evolution with a sample of 18 z~6 quasars that are faint at UV/optical wavelengths (rest-frame 1450A magnitudes of m_1450\ge20.2). The average FIR-to-AGN UV luminosity ratio of this faint quasar sample is about two times higher than that of the bright quasars at z~6 (m_1450<20.2). A fit to the average FIR and AGN bolometric luminosities of both the UV/optically faint and bright z~6 quasars, and the average luminosities of samples of submillimeter /millimeter-observed quasars at z~2 to 5, yields a relationship of L_{FIR} {L_{bol}}^{0.62}. Five of the 18 faint z~6 quasars have been detected at 250 GHz. These 250 GHz detections, as well as most of the millimeter-detected optically bright z~6 quasars, follow a shallower trend of L_{FIR} {L_{bol}}^{0.45} defined by the starburst-AGN systems in local and high-z universe. The millimeter continuum detections in the five objects and molecular CO detections in three of them reveal a few x10^8 M_sun of FIR-emitting warm dust and 10^10 M_sun of molecular gas in the quasar host galaxies. All these results argue for massive star formation in the quasar host galaxies, with estimated star formation rates of a few hundred M_sun yr^{-1}. Additionally, the higher FIR-to-AGN luminosity ratio found in these 250 GHz-detected faint quasars also suggests a higher ratio between star formation rate and supermassive black hole accretion rate than the UV/optically most luminous quasars at z~6.

preprint2011arXiv

Imaging the Molecular Gas Properties of a Major Merger Driving the Evolution of a z=2.5 Submillimeter Galaxy

We report the detection of spatially extended CO 1-0 and 5-4 emission in the z=2.49 submillimeter galaxy (SMG) J123707+6214, using the Expanded Very Large Array and the Plateau de Bure Interferometer. The large molecular gas reservoir is spatially resolved into two CO(1-0) components (north-east and south-west; previously identified in CO 3-2 emission) with gas masses of 4.3 and 3.5 x 10^10 (alpha_CO/0.8) Msun. We thus find that the optically invisible north-east component slightly dominates the gas mass in this system. The total molecular gas mass derived from the CO(1-0) observations is ~2.5 times larger than estimated from CO(3-2). The two components are at approximately the same redshift, but separated by ~20 kpc in projection. The morphology is consistent with that of an early-stage merger. The total amount of molecular gas is sufficient to maintain the intense 500 Msun/yr starburst in this system for at least ~160 Myr. We derive line brightness temperature ratios of r_31=0.39+/-0.09 and 0.37+/-0.10, and r_51=0.26+/-0.07 and 0.25+/-0.08 in the two components, respectively, suggesting that the J>=3 lines are substantially subthermally excited. This also suggests comparable conditions for star formation in both components. Given the similar gas masses of both components, this is consistent with the comparable starburst strengths observed in the radio continuum emission. Our findings are consistent with other recent studies that find evidence for lower CO excitation in SMGs than in high-z quasar host galaxies with comparable gas masses. This may provide supporting evidence that both populations correspond to different evolutionary stages in the formation of massive galaxies.

preprint2011arXiv

Interaction between dark matter sub-halos & galactic gaseous disk

We investigate the idea that the interaction of Dark Matter (DM) sub-halos with the gaseous disk of galaxies can be the origin for the observed holes and shells found in their neutral hydrogen (HI) distributions. We use high resolution hydrodynamic simulations to show that pure DM sub-halos impacting a galactic disk are not able to produce holes; on the contrary, they result in high density regions in the disk. However, sub-halos containing a small amount of gas (a few percent of the total DM mass of the sub-halo) are able to displace the gas in the disk and form holes and shells. The size and lifetime of these holes depend on the sub-halo gas mass, density and impact velocity. A DM sub-halo of mass $10^8$ $M_{\odot}$ and a gas mass fraction of $\sim 3$%, is able to create a kpc scale hole, with a lifetime similar to those observed in nearby galaxies. We also register an increase in the star formation rate at the rim of the hole, again in agreement with observations. Even though the properties off these simulated structures resemble those found in observations we find that the number of predicted holes (based on mass and orbital distributions of DM halos derived from cosmological N-body simulations) falls short compared to the observations. Only a handful of holes are produced per Giga year. This leads us to conclude that DM halo impact is not the major channel through which these holes are formed.

preprint2011arXiv

Molecular Gas at High Redshift

In order to understand galaxy evolution through cosmic times it is critical to derive the properties of the molecular gas content of galaxies, i.e. the material out of which stars ultimately form. The last decade has seen rapid progress in this area, with the detection of massive molecular gas reservoirs at high redshifts in submillimeter-selected galaxies and quasars. In the latter case, molecular gas reservoirs have been quantified out to redshifts z>6, i.e. towards the end of cosmic reionization when the universe was less than one Gyr old. The recent discovery of molecular gas in more normal galaxies have extended these studies from the most extreme objects in the universe (SFR~1000M_sun/yr; quasars and submillimeter galaxies) to more 'normal' starforming systems at redshifts 1.5-2.5 (with SFR~100M_sun/yr). However, detecting the molecular gas reservoirs of high-redshift galaxies that only have moderate star formation rates (~<10 M_sun/yr, similar to the faint galaxies seen in the Hubble Ultra Deep Field) will likely have to await the completion of ALMA.

preprint2011arXiv

Observational Evidence Against Long-Lived Spiral Arms in Galaxies

We test whether the spiral patterns apparent in many large disk galaxies should be thought of as dynamical features that are stationary in a co-rotating frame for > t_{dyn}, as implied by the density wave approach for explaining spiral arms. If such spiral arms have enhanced star formation (SF), observational tracers for different stages of the SF sequence should show a spatial ordering, from up-stream to downstream in the corotating frame: dense HI, CO, tracing molecular hydrogen gas, 24 micron emission tracing enshrouded SF and UV emission tracing unobscured young stars. We argue that such a spatial ordering should be reflected in the angular cross-correlation (CC, in polar coordinates) using all azimuthal positions among pairs of these tracers; the peak of the CC should be offset from zero, in different directions inside and outside the corotation radius. Recent spiral SF simulations by Dobbs & Pringle, show explicitly that for the case of a stationary spiral arm potential such angular offsets between gas and young stars of differing ages should be observable as cross-correlation offsets. We calculate the angular cross-correlations for different observational SF sequence tracers in 12 nearby spiral galaxies, drawing on a data set with high quality maps of the neutral gas HI, THINGS), molecular gas (CO, HERACLES) along with 24 micron emission (Spitzer, SINGS); we include FUV images (GALEX) and 3.6 $μ$m emission (Spitzer, IRAC) for some galaxies, tracing aging stars and longer timescales. In none of the resulting tracer cross-correlations for this sample do we find systematic angular offsets, which would be expected for a stationary dynamical spiral pattern of well-defined pattern speed. This result indicates that spiral density waves in their simplest form are not an important aspect of explaining spirals in large disk galaxies.

preprint2011arXiv

The central slope of dark matter cores in dwarf galaxies: Simulations vs. THINGS

We make a direct comparison of the derived dark matter (DM) distributions between hydrodynamical simulations of dwarf galaxies assuming a LCDM cosmology and the observed dwarf galaxies sample from the THINGS survey in terms of (1) the rotation curve shape and (2) the logarithmic inner density slope alpha of mass density profiles. The simulations, which include the effect of baryonic feedback processes, such as gas cooling, star formation, cosmic UV background heating and most importantly physically motivated gas outflows driven by supernovae (SNe), form bulgeless galaxies with DM cores. We show that the stellar and baryonic mass is similar to that inferred from photometric and kinematic methods for galaxies of similar circular velocity. Analyzing the simulations in exactly the same way as the observational sample allows us to address directly the so-called "cusp/core" problem in the LCDM model. We show that the rotation curves of the simulated dwarf galaxies rise less steeply than CDM rotation curves and are consistent with those of the THINGS dwarf galaxies. The mean value of the logarithmic inner density slopes alpha of the simulated galaxies' dark matter density profiles is ~ -0.4 +- 0.1, which shows good agreement with α= -0.29 +- 0.07 of the THINGS dwarf galaxies. The effect of non-circular motions is not significant enough to affect the results. This confirms that the baryonic feedback processes included in the simulations are efficiently able to make the initial cusps with α~ -1.0 to -1.5 predicted by dark-matter-only simulations shallower, and induce DM halos with a central mass distribution similar to that observed in nearby dwarf galaxies.

preprint2011arXiv

The Displaced Dusty ISM of NGC 3077: Tidal Stripping in the M 81 Triplet

We present the detection of extended (~30 kpc^2) dust emission in the tidal \hi\ arm near NGC 3077 (member of the M\,81 triplet) using SPIRE on board Herschel. Dust emission in the tidal arm is typically detected where the \hi\ column densities are >10^21 cm^-2. The SPIRE band ratios show that the dust in the tidal arm is significantly colder (~13 K) than in NGC 3077 itself (~31 K), consistent with the lower radiation field in the tidal arm. The total dust mass in the tidal arm is ~1.8 x 10^6 M_sun (assuming beta=2), i.e. substantially larger than the dust mass associated with NGC 3077 (~2 x 10^5 M_sun). Where dust is detected, the dust-to-gas ratio is 6+/-3 x 10^-3, consistent within the uncertainties with what is found in NGC 3077 and nearby spiral galaxies with Galactic metallicities. The faint HII regions in the tidal arm can not be responsible for the detected enriched material and are not the main source of the dust heating in the tidal arm. We conclude that the interstellar medium (atomic HI, molecules and dust) in this tidal feature was pre-enriched and stripped off NGC 3077 during its recent interaction (~3 x 10^8 yr ago) with M 82 and M 81. This implies that interaction can efficiently remove heavy elements and enriched material (dust, molecular gas) from galaxies. As interactions were more frequent at large lookback times, it is conceivable that they could substantially contribute (along with galactic outflows) to the enrichment of the intergalactic medium.

preprint2011arXiv

The Emission by Dust and Stars of Nearby Galaxies in the Herschel KINGFISH Survey

Using new far-infrared imaging from the Herschel Space Observatory with ancillary data from ultraviolet to submillimeter wavelengths, we estimate the total emission from dust and stars of 62 nearby galaxies in the KINGFISH survey in a way that is as empirical and model-independent as possible. We collect and exploit these data in order to measure from the spectral energy distributions (SEDs) precisely how much stellar radiation is intercepted and re-radiated by dust, and how this quantity varies with galaxy properties. By including SPIRE data, we are more sensitive to emission from cold dust grains than previous analyses at shorter wavelengths, allowing for more accurate estimates of dust temperatures and masses. The dust/stellar flux ratio, which we measure by integrating the SEDs, has a range of nearly three decades. The inclusion of SPIRE data shows that estimates based on data not reaching these far-IR wavelengths are biased low. We find that the dust/stellar flux ratio varies with morphology and total IR luminosity. We also find that dust/stellar flux ratios are related to gas-phase metallicity, while the dust/stellar mass ratios are less so. The substantial scatter between dust/stellar flux and dust/stellar mass indicates that the former is a poor proxy of the latter. Comparing the dust/stellar flux ratios and dust temperatures, we show that early-types tend to have slightly warmer temperatures than spiral galaxies, which may be due to more intense interstellar radiation fields, or to different dust grain compositions. Finally, we show that early-types and early-type spirals have a strong correlation between the dust/stellar flux ratio and specific star formation rate, which suggests that the relatively bright far-IR emission of some of these galaxies is due to ongoing star formation and the radiation field from older stars.

preprint2011arXiv

The First High Redshift Quasar from Pan-STARRS

We present the discovery of the first high redshift (z > 5.7) quasar from the Panoramic Survey Telescope and Rapid Response System 1 (Pan-STARRS1 or PS1). This quasar was initially detected as an i dropoutout in PS1, confirmed photometrically with the SAO Widefield InfraRed Camera (SWIRC) at Arizona's Multiple Mirror Telescope (MMT) and the Gamma-Ray Burst Optical/Near-Infrared Detector (GROND) at the MPG 2.2 m telescope in La Silla. The quasar was verified spectroscopically with the the MMT Spectrograph, Red Channel and the Cassegrain Twin Spectrograph (TWIN) at the Calar Alto 3.5 m telescope. It has a redshift of 5.73, an AB z magnitude of 19.4, a luminosity of 3.8 x 10^47 erg/s and a black hole mass of 6.9 x 10^9 solar masses. It is a Broad Absorption Line quasar with a prominent Ly-beta peak and a very blue continuum spectrum. This quasar is the first result from the PS1 high redshift quasar search that is projected to discover more than a hundred i dropout quasars, and could potentially find more than 10 z dropout (z > 6.8) quasars.

preprint2011arXiv

The Formation of Kiloparsec-Scale HI Holes in Dwarf Galaxies

The origin of kpc-scale holes in the atomic hydrogen (H I) distributions of some nearby dwarf irregular galaxies presents an intriguing problem. Star formation histories (SFHs) derived from resolved stars give us the unique opportunity to study past star forming events that may have helped shape the currently visible H I distribution. Our sample of five nearby dwarf irregular galaxies spans over an order of magnitude in both total H I mass and absolute B-band magnitude and is at the low mass end of previously studied systems. We use Very Large Array H I line data to estimate the energy required to create the centrally dominant hole in each galaxy. We compare this energy estimate to the past energy released by the underlying stellar populations computed from SFHs derived from data taken with the Hubble Space Telescope. The inferred integrated stellar energy released within the characteristic ages exceeds our energy estimates for creating the holes in all cases, assuming expected efficiencies. Therefore, it appears that stellar feedback provides sufficient energy to produce the observed holes. However, we find no obvious signature of single star forming events responsible for the observed structures when comparing the global SFHs of each galaxy in our sample to each other or to those of dwarf irregular galaxies reported in the literature. We also fail to find evidence of a central star cluster in FUV or Halpha imaging. We conclude that large H I holes are likely formed from multiple generations of star formation and only under suitable interstellar medium conditions.

preprint2011arXiv

The M81 Group Dwarf Irregular Galaxy DDO 165. I. High Velocity Neutral Gas in a Post-Starburst System

We present new multi-configuration VLA HI spectral line observations of the M81 group dIrr post-starburst galaxy DDO 165. The HI morphology is complex, with multiple column density peaks surrounding a large region of very low HI surface density that is offset from the center of the stellar distribution. The bulk of the neutral gas is associated with the southern section of the galaxy; a secondary peak in the north contains ~15% of the total HI mass. These components appear to be kinematically distinct, suggesting that either tidal processes or large-scale blowout have recently shaped the ISM of DDO 165. Using spatially-resolved position-velocity maps, we find multiple localized high-velocity gas features. Cross-correlating with radius-velocity analyses, we identify eight shell/hole structures in the ISM with a range of sizes (~400-900 pc) and expansion velocities (~7-11 km/s). These structures are compared with narrow- and broad-band imaging from KPNO and HST. Using the latter data, recent works have shown that DDO 165's previous "burst" phase was extended temporally (>1 Gyr). We thus interpret the high-velocity gas features, HI holes, and kinematically distinct components of the galaxy in the context of the immediate effects of "feedback" from recent star formation. In addition to creating HI holes and shells, extended star formation events are capable of creating localized high velocity motion of the surrounding interstellar material. A companion paper connects the energetics from the HI and HST data.

preprint2011arXiv

The M81 Group Dwarf Irregular Galaxy DDO 165. II. Connecting Recent Star Formation with ISM Structures and Kinematics

We compare the stellar populations and complex neutral gas dynamics of the M81 group dIrr galaxy DDO 165 using data from the HST and the VLA. Paper I identified two kinematically distinct HI components, multiple localized high velocity gas features, and eight HI holes and shells (the largest of which spans ~2.2x1.1 kpc). Using the spatial and temporal information from the stellar populations in DDO 165, we compare the patterns of star formation over the past 500 Myr with the HI dynamics. We extract localized star formation histories within 6 of the 8 HI holes identified in Paper I, as well as 23 other regions that sample a range of stellar densities and neutral gas properties. From population synthesis modeling, we derive the energy outputs (from stellar winds and supernovae) of the stellar populations within these regions over the last 100 Myr, and compare with refined estimates of the energies required to create the HI holes. In all cases, we find that "feedback" is energetically capable of creating the observed structures in the ISM. Numerous regions with significant energy inputs from feedback lack coherent HI structures but show prominent localized high velocity gas features; this feedback signature is a natural product of temporally and spatially distributed star formation. In DDO 165, the extended period of heightened star formation activity (lasting more than 1 Gyr) is energetically capable of creating the observed holes and high velocity gas features in the neutral ISM.

preprint2011arXiv

Water vapor emission reveals a highly obscured, star forming nuclear region in the QSO host galaxy APM08279+5255 at z=3.9

We present the detection of four rotational emission lines of water vapor, from energy levels Eu/k= 101 - 454 K, in the gravitationally lensed z=3.9 QSO host galaxy APM08279+5255. While the lowest H2O lines are collisionally excited in clumps of warm, dense gas (density of hydrogen nuclei n_H=(3.1 +/- 1.2) x 10^6 cm^-3, gas temperature T_g ~ 105 +/- 21 K), we find that the excitation of the higher lines is dominated by the intense local infrared radiation field. Since only collisionally excited emission contributes to gas cooling, we conclude that H2O is not a significant coolant of the warm molecular gas. Our excitation model requires the radiatively excited gas to be located in an extended region of high 100 micron opacity (tau_100 = 0.9 +/- 0.2). Locally, such extended infrared-opaque regions are found only in the nuclei of ultraluminous infrared galaxies. We propose a model where the infrared-opaque circumnuclear cloud, which is penetrated by the X-ray radiation field of the QSO nucleus, contains clumps of massive star formation where the H2O emission originates. The radiation pressure from the intense local infrared radiation field exceeds the thermal gas pressure by about an order of magnitude, suggesting close to Eddington-limited star formation in these clumps.

preprint2010arXiv

[CII] line emission in BRI1335-0417 at z=4.4

Using the 12m APEX telescope, we have detected redshifted emission from the 157.74micron [CII] line in the z=4.4074 quasar BRI1335-0417. The linewidth and redshift are in good agreement with previous observations of high-J CO line emission. We measure a [CII] line luminosity, L_[CII] = (16.4 +/- 2.6)x10^9 Lsun, making BRI~1335-0417 the most luminous, unlensed [CII] line emitter known at high-redshift. The [CII]-to-FIR luminosity ratio of (5.3+/-0.8)x10^-4 is ~3x higher than expected for an average object with a FIR luminosity L_FIR = 3.1x10^13 Lsun, if this ratio were to follow the trend observed in other FIR-bright galaxies that have been detected in [CII] line emission. These new data suggest that the scatter in the [CII]-to-FIR luminosity ratio could be larger than previously expected for high luminosity objects. BR1335-0417 has a similar FIR luminosity and [CII]/CO luminosity compared to local ULIRGS and appears to be a gas-rich merger forming stars at a rate of a few thousand solar masses per year.

preprint2010arXiv

A joint analysis of BLAST 250--500um and LABOCA 870um observations in the Extended Chandra Deep Field South

We present a joint analysis of the overlapping BLAST 250, 350, 500um, and LABOCA 870um observations (from the LESS survey) of the ECDF-S. Out to z~3, the BLAST filters sample near the peak wavelength of far-infrared (FIR) emission from galaxies (rest-frame wavelengths ~60--200um), primarily produced by dust heated through absorption in star-forming clouds. However, identifying counterparts to individual BLAST peaks is very challenging, given the large beams (FWHM 36--60"). In contrast, the 870um observations have a significantly smaller 19" FWHM beam, and are sensitive to higher redshifts (z~1--5, and potentially beyond) due to the more favourable negative K-correction. We use the LESS data, as well as deep Spitzer and VLA imaging, to identify 118 individual sources that produce significant emission in the BLAST bands. We characterize the temperatures and FIR luminosities for a subset of 69 sources which have well-measured submm SEDs and redshift measurements out to z~3. For flux-limited sub-samples in each BLAST band, and a dust emissivity index β=2.0, we find a median temperature T=30K (all bands) as well as median redshifts: z=1.1 (interquartile range 0.2--1.9) for S250 > 40mJy; z=1.3 (interquartile range 0.6--2.1) for S350 > 30mJy; and z=1.6 (interquartile range 1.3--2.3) for S500 > 20mJy. Taking into account selection effects for our survey (a bias toward detecting lower-temperature galaxies), we find no evidence for evolution in the local FIR-temperature correlation out to z~2.5. Comparing with star-forming galaxy SED templates, about 8% of our sample appears to exhibit significant excesses in the radio and/or mid-IR, consistent with those sources harbouring an AGN. We describe the following techniques in two appendices: our `matched filter' for identifying sources in the presence of point-source confusion; and our approach for identifying counterparts using likelihood ratios.

preprint2010arXiv

Arm & Interarm Star Formation in Spiral Galaxies

We investigate the relationship between spiral arms and star formation in the grand-design spirals NGC 5194 and NGC 628 and in the flocculent spiral NGC 6946. Filtered maps of near-IR (3.6 micron) emission allow us to identify "arm regions" that should correspond to regions of stellar mass density enhancements. The two grand-design spirals show a clear two-armed structure, while NGC 6946 is more complex. We examine these arm and interarm regions, looking at maps that trace recent star formation - far-ultraviolet (GALEX NGS) and 24 micron emission (Spitzer, SINGS) - and cold gas - CO (Heracles) and HI (Things). We find the star formation tracers and CO more concentrated in the spiral arms than the stellar 3.6 micron flux. If we define the spiral arms as the 25% highest pixels in the filtered 3.6 micron images, we find that the majority (60%) of star formation tracers occurs in the interarm regions; this result persists qualitatively even when considering the potential impact of finite data resolution and diffuse interarm 24 micron emission. Even with a generous definition of the arms (45% highest pixels), interarm regions still contribute at least 30% to the integrated star formation rate tracers. We look for evidence that spiral arms trigger star or cloud formation using the ratios of star formation rate (SFR, traced by a combination of FUV and 24 micron emission) to H_2 (traced by CO) and H_2 to HI. Any enhancement of SFR / M(H_2) in the arm region is very small (less than 10%) and the grand design spirals show no enhancement compared to the flocculent target. Arm regions do show a weak enhancement in H_2/HI compared to the interarm regions, but at a fixed gas surface density there is little clear enhancement in the H_2/HI ratio in the arm regions. Thus, it seems that spiral arms may only act to concentrate the gas to higher densities in the arms.

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

Dense Molecular Gas Excitation at High Redshift: Detection of HCO+(J=4-3) Emission in the Cloverleaf Quasar

We report the detection of HCO+(J=4-3) emission in the Cloverleaf Quasar at z=2.56, using the IRAM Plateau de Bure Interferometer. HCO+ emission is a star formation indicator similar to HCN, tracing dense molecular hydrogen gas (n(H2) ~= 10^5 cm^-3) within star-forming molecular clouds. We derive a lensing-corrected HCO+(J=4-3) line luminosity of L'(HCO+(4-3)) = (1.6+/-0.3) x 10^9 (mu_L/11)^-1 K km/s pc^2, which corresponds to only 48% of the HCO+(J=1=0) luminosity, and <~4% of the CO(J=3-2) luminosity. The HCO+ excitation thus is clearly subthermal in the J=4-3 transition. Modeling of the HCO+ line radiative transfer suggests that the HCO+ emission emerges from a region with physical properties comparable to that exhibiting the CO line emission, but 2x higher gas density. This suggests that both HCO+ and CO lines trace the warm, dense molecular gas where star formation actively takes place. The HCO+ lines have only ~2/3 the width of the CO lines, which may suggest that the densest gas is more spatially concentrated. In contrast to the z=3.91 quasar APM08279+5255, the dense gas excitation in the Cloverleaf is consistent with being purely collisional, rather than being enhanced by radiative processes. Thus, the physical properties of the dense gas component in the Cloverleaf are consistent with those in the nuclei of nearby starburst galaxies. This suggests that the conditions in the dense, star-forming gas in active galactic nucleus-starburst systems at early cosmic times like the Cloverleaf are primarily affected by the starburst itself, rather than the central active black hole.

preprint2010arXiv

Dense Molecular Gas Excitation in Nuclear Starbursts at High Redshift: HCN, HNC, and HCO+(J=6-5) Emission in the z=3.91 Quasar Host of APM08279+5255

We report the detection of surprisingly strong HCN, HNC, and HCO+(J=6-5) emission in the host galaxy of the z=3.91 quasar APM08279+5255 through observations with CARMA. HCN, HNC, and HCO+ are typically used as star formation indicators, tracing dense molecular hydrogen gas [n(H2) > 10^5,cm^-3] within star-forming molecular clouds. However, the strength of their respective line emission in the J=6-5 transitions in APM08279+5255 is extremely high, suggesting that they are excited by another mechanism besides collisions in the dense molecular gas phase alone. We derive J=6-5 line luminosities of L'(HCN)=(4.9+/-0.6), L'(HNC)=(2.4+/-0.7), and L'(HCO+)=(3.0+/-0.6)x10^10 (mu_L)^-1 K km/s pc^2 (where mu_L is the lensing magnification factor), corresponding to L' ratios of ~0.23-0.46 relative to CO(J=1-0). Such high line ratios would be unusual even in the respective ground-state (J=1-0) transitions, and indicate exceptional, collisionally and radiatively driven excitation conditions in the dense, star-forming molecular gas in APM08279+5255. Through an expansion of our previous modeling of the HCN line excitation in this source, we show that the high rotational line fluxes are caused by substantial infrared pumping at moderate opacities in a ~220K warm gas and dust component. This implies that standard M_dense/L' conversion factors would substantially overpredict the dense molecular gas mass M_dense. We also find a HCN J=6-5/5-4 L' ratio greater than 1 (1.36+/-0.31) - however, our models show that the excitation is likely not `super-thermal', but that the high line ratio is due to a rising optical depth between both transitions. These findings are consistent with the picture that the bulk of the gas and dust in this source is situated in a compact, nuclear starburst, where both the highly active galactic nucleus and star formation contribute to the heating.

preprint2010arXiv

Detection of molecular gas in a distant submillimetre galaxy at z=4.76 with ATCA

We have detected the CO(2-1) transition from the submillimetre galaxy (SMG) LESSJ033229.4-275619 at z=4.755 using the new Compact Array Broadband Backend system on the Australian Telescope Compact Array. These data have identified a massive gas reservoir available for star formation for the first time in an SMG at z~5. We use the luminosity and velocity width (FWHM of 160 km/s) of the CO(2--1) line emission to constrain the gas and dynamical mass of Mgas~1.6x10^10 Msun and Mdyn(<2kpc)~5x10^10 (0.25/sin^2(i)) Msun, respectively, similar to that observed for SMGs at lower redshifts of z~2-4, although we note that our observed CO FWHM is a factor of ~3 narrower than typically seen in SMGs. Together with the stellar mass we estimate a total baryonic mass of Mbary~1x10^11 Msun, consistent with the dynamical mass for this young galaxy within the uncertainties. Dynamical and baryonic mass limits of high-redshift galaxies are useful tests of galaxy formation models: using the known z~4-5 SMGs as examples of massive baryonic systems, we find that their space density is consistent with that predicted by current galaxy formation models. In addition, these observations have helped to confirm that z~4-5 SMGs possess the baryonic masses and gas consumption timescales necessary to be the progenitors of the luminous old red galaxies seen at z~3. Our results provide a preview of the science that ALMA will enable on the formation and evolution of the earliest massive galaxies in the Universe.

preprint2010arXiv

Dust-free quasars in the early Universe

The most distant quasars known, at redshifts z=6, generally have properties indistinguishable from those of lower-redshift quasars in the rest-frame ultraviolet/optical and X-ray bands. This puzzling result suggests that these distant quasars are evolved objects even though the Universe was only seven per cent of its current age at these redshifts. Recently one z=6 quasar was shown not to have any detectable emission from hot dust, but it was unclear whether that indicated different hot-dust properties at high redshift or if it is simply an outlier. Here we report the discovery of a second quasar without hot-dust emission in a sample of 21 z=6 quasars. Such apparently hot-dust-free quasars have no counterparts at low redshift. Moreover, we demonstrate that the hot-dust abundance in the 21 quasars builds up in tandem with the growth of the central black hole, whereas at low redshift it is almost independent of the black hole mass. Thus z=6 quasars are indeed at an early evolutionary stage, with rapid mass accretion and dust formation. The two hot-dust-free quasars are likely to be first-generation quasars born in dust-free environments and are too young to have formed a detectable amount of hot dust around them.

preprint2010arXiv

Hierarchical Stellar Structures in the Local Group Dwarf Galaxy NGC 6822

We present a comprehensive study of the star cluster population and the hierarchical structure in the clustering of blue stars with ages <~ 500 Myr in the Local Group dwarf irregular galaxy NGC 6822. Our observational material comprises the most complete optical stellar catalog of the galaxy from imaging with the Suprime-Cam at the 8.2-m SUBARU Telescope. We identify 47 distinct star clusters with the application of the nearest-neighbor density method to this catalog for a detection threshold of 3sigma above the average stellar density. The size distribution of the detected clusters can be very well approximated by a Gaussian with a peak at ~ 68 pc. Their cluster mass function is fitted very well by a power-law with index alpha ~ 1.5 +/- 0.7, consistent with other Local Group galaxies and the cluster initial mass function. The application of the nearest-neighbor density method for various density thresholds, other than 3sigma, enabled the identification of stellar concentrations in various length-scales. The stellar density maps constructed with this technique provide a direct proof of hierarchically structured stellar concentrations in NGC 6822. We illustrate this hierarchy by the so-called "dendrogram" of the detected stellar structures, which demonstrates that most of the detected structures split up into several substructures over at least three levels. We quantify the hierarchy of these structures with the use of the minimum spanning tree method. The morphological hierarchy in stellar clustering, which we observe in NGC 6822 resembles that of the turbulent interstellar matter, suggesting that turbulence on pc- and kpc-scales has been probably the major agent that regulated clustered star formation in NGC 6822.

preprint2010arXiv

Molecular Gas in Redshift 6 Quasar Host Galaxies

We report our new observations of redshifted carbon monoxide emission from six z~6 quasars, using the PdBI. CO (6-5) or (5-4) line emission was detected in all six sources. Together with two other previous CO detections, these observations provide unique constraints on the molecular gas emission properties in these quasar systems close to the end of the cosmic reionization. Complementary results are also presented for low-J CO lines observed at the GBT and the VLA, and dust continuum from five of these sources with the SHARC-II bolometer camera at the CSO. We then present a study of the molecular gas properties in our combined sample of eight CO-detected quasars at z~6. The detections of high-order CO line emission in these objects indicates the presence of highly excited molecular gas, with estimated masses on the order of 10^10 M_sun within the quasar host galaxies. No significant difference is found in the gas mass and CO line width distributions between our z~6 quasars and samples of CO-detected $1.4\leq z\leq5$ quasars and submillimeter galaxies. Most of the CO-detected quasars at z~6 follow the far infrared-CO luminosity relationship defined by actively star-forming galaxies at low and high redshifts. This suggests that ongoing star formation in their hosts contributes significantly to the dust heating at FIR wavelengths. The result is consistent with the picture of galaxy formation co-eval with supermassive black hole (SMBH) accretion in the earliest quasar-host systems. We investigate the black hole--bulge relationships of our quasar sample, using the CO dynamics as a tracer for the dynamical mass of the quasar host. The results place important constraints on the formation and evolution of the most massive SMBH-spheroidal host systems at the highest redshift.

preprint2010arXiv

The Scale Dependence of the Molecular Gas Depletion Time in M33

We study the Local Group spiral galaxy M33 to investigate how the observed scaling between the (kpc-averaged) surface density of molecular gas (Σ_H2) and recent star formation rate (Σ_SFR) relates to individual star-forming regions. To do this, we measure the ratio of CO emission to extinction-corrected Halpha emission in apertures of varying sizes centered both on peaks of CO and Halpha emission. We parameterize this ratio as a molecular gas (H_2) depletion time (τ_dep). On large (kpc) scales, our results are consistent with a molecular star formation law (Sigma_SFR \sim Sigma_H2^b) with b \sim 1.1 - 1.5 and a median τ_dep \sim 1 Gyr, with no dependence on type of region targeted. Below these scales, τ_dep is a strong function of adopted angular scale and the type of region that is targeted. Small (\lesssim 300pc) apertures centered on CO peaks have very long τ_dep (i.e., high CO-to-Halpha flux ratio) and small apertures targeted toward Halpha peaks have very short τ_dep. This implies that the star formation law observed on kpc scales breaks down once one reaches aperture sizes of \lesssim 300pc. For our smallest apertures (75pc), the difference in τ_dep between the two types of regions is more than one order of magnitude. This scale behavior emerges from averaging over star-forming regions with a wide range of CO-to-Halpha ratios with the natural consequence that the breakdown in the star formation law is a function of the surface density of the regions studied. We consider the evolution of individual regions the most likely driver for region-to-region differences in τ_dep (and thus the CO-to-Halpha ratio).

preprint2010arXiv

Total Molecular Gas Masses of z~3 Lyman-Break Galaxies: CO(1-0) Emission in MS1512-cB58 and the Cosmic Eye

We report the detection of CO(1-0) emission toward the lensed L*_UV Lyman-break galaxies (LBGs) MS1512-cB58 (z=2.73) and the Cosmic Eye (z=3.07), using the Expanded Very Large Array. The strength of the CO line emission reveals molecular gas reservoirs with masses of (4.6+/-1.1) x 10^8 (mu_L/32)^-1 (alpha_CO/0.8) Msun and (9.3+/-1.6) x 10^8 (mu_L/28)^-1 (alpha_CO/0.8) Msun, respectively. These observations suggest by ~30%-40% larger gas reservoirs than estimated previously based on CO(3-2) observations due to subthermal excitation of the J=3 line. These observations also suggest gas mass fractions of 0.46+/-0.17 and 0.16+/-0.06. The CO(1-0) emission in the Cosmic Eye is slightly resolved on scales of 4.5"+/-1.5", consistent with previous studies of nebular emission lines. This suggests that the molecular gas is associated with the most intensely star-forming regions seen in the ultraviolet (UV). We do not resolve the CO(1-0) emission in cB58 at ~2" resolution, but find that the CO(1-0) emission is also consistent with the position of the UV-brightest emission peak. The gas masses, gas fractions, moderate CO line excitation, and star formation efficiencies in these galaxies are consistent with what is found in nearby luminous infrared galaxies. These observations thus currently represent the best constraints on the molecular gas content of `ordinary' (i.e., ~L*_UV) z~3 star-forming galaxies. Despite comparable star formation rates, the gas properties of these young LBGs seem to be different from the recently identified optical/infrared-selected high-z massive, gas-rich star-forming galaxies, which are more gas-rich and massive, but have lower star formation efficiencies, and presumably trace a different galaxy population.

preprint2009arXiv

BLAST: the far-infrared/radio correlation in distant galaxies

We investigate the correlation between far-infrared (FIR) and radio luminosities in distant galaxies, a lynchpin of modern astronomy. We use data from the Balloon-borne Large Aperture Submillimetre Telescope (BLAST), Spitzer, the Large Apex BOlometer CamerA (LABOCA), the Very Large Array (VLA) and the Giant Metre-wave Radio Telescope (GMRT) in the Extended Chandra Deep Field South (ECDFS). For a catalogue of BLAST 250-micron-selected galaxies, we re-measure the 70--870-micron flux densities at the positions of their most likely 24-micron counterparts, which have a median [interquartile] redshift of 0.74 [0.25, 1.57]. From these, we determine the monochromatic flux density ratio, q_250 = log_10 (S_250micron / S_1400MHz), and the bolometric equivalent, q_IR. At z ~= 0.6, where our 250-micron filter probes rest-frame 160-micron emission, we find no evolution relative to q_160 for local galaxies. We also stack the FIR and submm images at the positions of 24-micron- and radio-selected galaxies. The difference between q_IR seen for 250-micron- and radio-selected galaxies suggests star formation provides most of the IR luminosity in ~< 100-uJy radio galaxies, but rather less for those in the mJy regime. For the 24-micron sample, the radio spectral index is constant across 0 < z < 3, but q_IR exhibits tentative evidence of a steady decline such that q_IR is proportional to (1+z)^(-0.15 +/- 0.03) - significant evolution, spanning the epoch of galaxy formation, with major implications for techniques that rely on the FIR/radio correlation. We compare with model predictions and speculate that we may be seeing the increase in radio activity that gives rise to the radio background.

preprint2009arXiv

Comparison of H-alpha and UV Star Formation Rates in the Local Volume: Systematic Discrepancies for Dwarf Galaxies

(abridged) Using a complete sample of ~300 star-forming galaxies within 11 Mpc, we evaluate the consistency between star formation rates (SFRs) inferred from the far ultraviolet (FUV) non-ionizing continuum and H-alpha nebular emission, assuming standard conversion recipes in which the SFR scales linearly with luminosity at a given wavelength. Our analysis probes SFRs over 5 orders of magnitude, down to ultra-low activities on the order of ~0.0001 M_sun/yr. The data are drawn from the 11 Mpc H-alpha and Ultraviolet Galaxy Survey (11HUGS), which has obtained H-alpha fluxes from ground-based narrowband imaging, and UV fluxes from imaging with GALEX. For normal spiral galaxies (SFR~1 M_sun/yr), our results are consistent with previous work which has shown that FUV SFRs tend to be lower than H-alpha SFRs before accounting for internal dust attenuation, but that there is relative consistency between the two tracers after proper corrections are applied. However, a puzzle is encountered at the faint end of the luminosity function. As lower luminosity dwarf galaxies, roughly less active than the Small Magellanic Cloud, are examined, H-alpha tends to increasingly under-predict the SFR relative to the FUV. Although past studies have suggested similar trends, this is the first time this effect is probed with a statistical sample for galaxies with SFR~<0.1 M_sun/yr. A range of standard explanations does not appear to be able to account for the magnitude of the systematic. Some recent work has argued for an IMF which is deficient in high mass stars in dwarf and low surface brightness galaxies, and we also consider this scenario.

preprint2009arXiv

Does Stellar Feedback Create HI Holes? An HST/VLA Study of Holmberg II

We use deep HST/ACS F555W and F814W photometry of resolved stars in the M81 Group dwarf irregular galaxy Ho II to study the hypothesis that the holes identified in the neutral ISM (HI) are created by stellar feedback. From the deep photometry, we construct color-magnitude diagrams (CMDs) and measure the star formation histories (SFHs) for stars contained in HI holes from two independent holes catalogs, as well as select control fields, i.e., similar sized regions that span a range of HI column densities. Converting the recent SFHs into stellar feedback energies, we find that enough energy has been generated to have created all holes. However, the required energy is not always produced over a time scale that is less than the estimated kinematic age of the hole. The combination of the CMDs, recent SFHs, and locations of young stars shows that the stellar populations inside HI holes are not coherent, single-aged, stellar clusters, as previously suggested, but rather multi-age populations distributed across each hole. From a comparison of the modeled and observed integrated magnitudes, and the locations and energetics of stars inside of HI holes, we propose a potential new model: a viable mechanism for creating the observed HI holes in Ho II is stellar feedback from multiple generations of SF spread out over tens or hundreds of Myr, and thus, the concept of an age for an HI hole is intrinsically ambiguous. We further find that \halpha and 24 micron emission, tracers of the most recent star formation, do not correlate well with the positions of the HI holes. However, UV emission, which traces star formation over roughly the last 100 Myr, shows a much better correlation with the locations of the HI holes.

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.

preprint2008arXiv

A Molecular Einstein Ring at z=4.12: Imaging the Dynamics of a Quasar Host Galaxy Through a Cosmic Lens

We present high-resolution (0.3") Very Large Array (VLA) imaging of the molecular gas in the host galaxy of the high redshift quasar PSS J2322+1944 (z=4.12). These observations confirm that the molecular gas (CO) in the host galaxy of this quasar is lensed into a full Einstein ring, and reveal the internal dynamics of the molecular gas in this system. The ring has a diameter of ~1.5", and thus is sampled over ~20 resolution elements by our observations. Through a model-based lens inversion, we recover the velocity gradient of the molecular reservoir in the quasar host galaxy of PSS J2322+1944. The Einstein ring lens configuration enables us to zoom in on the emission and to resolve scales down to ~1 kpc. From the model-reconstructed source, we find that the molecular gas is distributed on a scale of 5 kpc, and has a total mass of M(H2)=1.7 x 10^10 M_sun. A basic estimate of the dynamical mass gives M_dyn = 4.4 x 10^10 (sin i)^-2 M_sun, that is, only ~2.5 times the molecular gas mass, and ~30 times the black hole mass (assuming that the dynamical structure is highly inclined). The lens configuration also allows us to tie the optical emission to the molecular gas emission, which suggests that the active galactic nucleus (AGN) does reside within, but not close to the center of the molecular reservoir. Together with the (at least partially) disturbed structure of the CO, this suggests that the system is interacting. Such an interaction, possibly caused by a major `wet' merger, may be responsible for both feeding the quasar and fueling the massive starburst of 680 M_sun/yr in this system, in agreement with recently suggested scenarios of quasar activity and galaxy assembly in the early universe.

preprint2008arXiv

A Search for Molecular Gas in the Nucleus of M87 and Implications for the Fueling of Supermassive Black Holes

Supermassive black holes in giant elliptical galaxies are remarkably faint given their expected accretion rates. This motivates models of radiatively inefficient accretion, due to either ion-electron thermal decoupling, generation of outflows that inhibit accretion, or settling of gas to a gravitationally unstable disk that forms stars in preference to feeding the black hole. The latter model predicts the presence of cold molecular gas in a thin disk around the black hole. Here we report Submillimeter Array observations of the nucleus of the giant elliptical galaxy M87 that probe 230 GHz continuum and CO(J=2--1) line emission. Continuum emission is detected from the nucleus and several knots in the jet, including one that has been undergoing flaring behavior. We estimate a conservative upper limit on the mass of molecular gas within ~100pc and +-400km/s line of sight velocity of the central black hole of ~8x10^6Msun, which includes an allowance for possible systematic errors associated with subtraction of the continuum. Ignoring such errors, we have a 3 sigma sensitivity to about 3x10^6Msun. In fact, the continuum-subtracted spectrum shows weak emission features extending up to 4 sigma above the RMS dispersion of the line-free channels. These may be artifacts of the continuum subtraction process. Alternatively, if they are interpreted as CO emission, then the implied molecular gas mass is ~5x10^6Msun spread out over a velocity range of 700km/s. These constraints on molecular gas mass are close to the predictions of the model of self-gravitating, star-forming accretion disks fed by Bondi accretion (Tan & Blackman 2005).

preprint2008arXiv

Formation of a Quasar Host Galaxy through a Wet Merger 1.4 Billion Years after the Big Bang

We present high-resolution Very Large Array imaging of the molecular gas in the host galaxy of the high redshift quasar BRI 1335-0417 (z=4.41). Our CO(2-1) observations have a linear resolution of 0.15" (1.0 kpc) and resolve the molecular gas emission both spatially and in velocity. The molecular gas in BRI 1335-0417 is extended on scales of 5 kpc, and shows a complex structure. At least three distinct components encompassing about two thirds of the total molecular mass of 9.2 x 10^10 M_sun are identified in velocity space, which are embedded in a structure that harbors about one third of the total molecular mass in the system. The brightest CO(2-1) line emission region has a peak brightness temperature of 61+/-9 K within 1 kpc diameter, which is comparable to the kinetic gas temperature as predicted from the CO line excitation. This is also comparable to the gas temperatures found in the central regions of nearby ultra-luminous infrared galaxies, which are however much more compact than 1 kpc. The spatial and velocity structure of the molecular reservoir in BRI 1335-0417 is inconsistent with a simple gravitationally bound disk, but resembles a merging system. Our observations are consistent with a major, gas-rich (`wet') merger that both feeds an accreting supermassive black hole (causing the bright quasar activity), and fuels a massive starburst that builds up the stellar bulge in this galaxy. Our study of this z>4 quasar host galaxy may thus be the most direct observational evidence that `wet' mergers at high redshift are related to AGN activity.

preprint2008arXiv

The Star Formation Efficiency in Nearby Galaxies: Measuring Where Gas Forms Stars Effectively

We measure the star formation efficiency (SFE), the star formation rate per unit gas, in 23 nearby galaxies and compare it to expectations from proposed star formation laws and thresholds. We use HI maps from THINGS and derive H2 maps from HERACLES and BIMA SONG CO. We estimate the star formation rate by combining GALEX FUV maps and SINGS 24 micron maps, infer stellar surface density profiles from SINGS 3.6 micron data, and use kinematics from THINGS. We measure the SFE as a function of: the free-fall and orbital timescales; midplane gas pressure; stability of the gas disk to collapse (including the effects of stars); the ability of perturbations to grow despite shear; and the ability of a cold phase to form. In spirals, the SFE of H2 alone is nearly constant at 5.25 +/- 2.5 x 10^(-10) yr^(-1) (equivalent to an H2 depletion time of 1.9x10^9 yr) as a function of all of these variables at our 800 pc resolution. Where the ISM is mostly HI, on the other hand, the SFE decreases with increasing radius in both spiral and dwarf galaxies, a decline reasonably described by an exponential with scale length 0.2-0.25 r_25. We interpret this decline as a strong dependence of GMC formation on environment. The ratio of H2 to HI appears to be a smooth function of radius, stellar surface density, and pressure spanning from the H2-dominated to HI-dominated ISM. The radial decline in SFE is too steep to be reproduced only by increases in the free-fall time or orbital time. Thresholds for large-scale instability suggest that our disks are stable or marginally stable and do not show a clear link to the declining SFE. We suggest that ISM physics below the scales that we observe - phase balance in the HI, H2 formation and destruction, and stellar feedback - governs the formation of GMCs from HI.

preprint2008arXiv

The Star Formation Law in Nearby Galaxies on Sub-Kpc Scales

(Abridged) We present a comprehensive analysis of the relationship between star formation rate surface density (SFR SD) and gas surface density (gas SD) at sub-kpc resolution in a sample of 18 nearby galaxies. We use high resolution HI data from THINGS, CO data from HERACLES and BIMA SONG, 24 micron data from the Spitzer Space Telescope, and UV data from GALEX. We target 7 spiral galaxies and 11 late-type/dwarf galaxies and investigate how the star formation law differs between the H2-dominated centers of spiral galaxies, their HI-dominated outskirts and the HI-rich late-type/dwarf galaxies. We find that a Schmidt-type power law with index N=1.0+-0.2 relates the SFR SD and the H2 SD across our sample of spiral galaxies, i.e., that H2 forms stars at a constant efficiency in spirals. The average molecular gas depletion time is ~2*10^9 yrs. We interpret the linear relation and constant depletion time as evidence that stars are forming in GMCs with approximately uniform properties and that the H2 SD may be more a measure of the filling fraction of giant molecular clouds than changing conditions in the molecular gas. The relationship between total gas SD and SFR SD varies dramatically among and within spiral galaxies. Most galaxies show little or no correlation between the HI SD and the SFR SD. As a result, the star formation efficiency (SFE = SFR SD / gas SD) varies strongly across our sample and within individual galaxies. We show that in spirals the SFE is a clear function of radius, while the dwarf galaxies in our sample display SFEs similar to those found in the outer optical disks of the spirals. Another general feature of our sample is a sharp saturation of the HI SD at ~9 M_sol/pc^2 in both the spiral and dwarf galaxies.

preprint2007arXiv

Star Formation in NGC 5194 (M51a). II. The Spatially-Resolved Star Formation Law

We have studied the relationship between the star formation rate (SFR) surface density and gas surface density in the spiral galaxy M51a (NGC 5194), using multi-wavelength data obtained as part of the Spitzer Infrared Nearby Galaxies Survey (SINGS). We introduce a new SFR index based on a linear combination of H-alpha emission-line and 24 micron continuum luminosities, that provides reliable extinction-corrected ionizing fluxes and SFR densities over a wide range of dust attenuations. The combination of these extinction-corrected SFR densities with aperture synthesis HI and CO maps has allowed us to probe the form of the spatially-resolved star formation law on scales of 0.5 to 2 kpc. We find that the resolved SFR vs gas surface density relation is well represented by a Schmidt power law, which is similar in form and dispersion to the disk-averaged Schmidt law. We observe a comparably strong correlation of the SFR surface density with the molecular gas surface density, but no significant correlation with the surface density of atomic gas. The best-fitting slope of the Schmidt law varies from N = 1.37 to 1.56, with zeropoint and slope that change systematically with the spatial sampling scale. We tentatively attribute these variations to the effects of areal sampling and averaging of a nonlinear intrinsic star formation law. Our data can also be fitted by an alternative parametrization of the SFR surface density in terms of the ratio of gas surface density to local dynamical time, but with a considerable dispersion.

preprint2006arXiv

A Search for H2O in the Strongly Lensed QSO MG 0751+2716 at z=3.2

We present a search for 183 GHz H_2O(3_13-2_20) emission in the infrared-luminous quasar MG 0751+2716 with the NRAO Very Large Array (VLA). At z=3.200+/-0.001, this water emission feature is redshifted to 43.6 GHz. As opposed to the faint rotational transitions of HCN (the standard high-density tracer at high-z), H_2O(3_13-2_20) is observed with high maser amplification factors in Galactic star-forming regions. It therefore holds the potential to trace high-density star-forming regions in the distant universe. If indeed all star-forming regions in massively star-forming galaxies at z>3 have similar physical properties as e.g. the Orion or W49N molecular cloud cores, the flux ratio between the maser-amplified H_2O(3_13-2_20) and the thermally excited CO(1-0) transitions may be as high as factor of 20 (but has to be corrected by their relative filling factor). MG 0751+2716 is a strong CO(4-3) emitter, and therefore one of the most suitable targets to search for H_2O(3_13-2_20) at cosmological redshifts. Our search resulted in an upper limit in line luminosity of L'(H_2O) < 0.6 x 10^9 K km/s pc^2. Assuming a brightness temperature of T_b(H_2O) ~= 500 K for the maser emission and CO properties from the literature, this translates to a H_2O(3_13-2_20)/CO(4-3) area filling factor of less than 1%. However, this limit is not valid if the H_2O(3_13-2_20) maser emission is quenched, i.e. if the line is only thermally excited. We conclude that, if our results were to hold for other high-z sources, H_2O does not appear to be a more luminous alternative to HCN to detect high-density gas in star-forming environments at high redshift.

preprint2006arXiv

First Detection of HCO+ Emission at High Redshift

We report the detection of HCO+(1-0) emission towards the Cloverleaf quasar (z=2.56) through observations with the Very Large Array. This is the first detection of ionized molecular gas emission at high redshift (z>2). HCO+ emission is a star formation indicator similar to HCN, tracing dense molecular hydrogen gas (n(H_2) ~= 10^5 cm^{-3}) within star-forming molecular clouds. We derive a lensing-corrected HCO+ line luminosity of L'(HCO+) = 3.5 x 10^9 K km/s pc^2. Combining our new results with CO and HCN measurements from the literature, we find a HCO+/CO luminosity ratio of 0.08 and a HCO+/HCN luminosity ratio of 0.8. These ratios fall within the scatter of the same relationships found for low-z star-forming galaxies. However, a HCO+/HCN luminosity ratio close to unity would not be expected for the Cloverleaf if the recently suggested relation between this ratio and the far-infrared luminosity were to hold. We conclude that a ratio between HCO+ and HCN luminosity close to 1 is likely due to the fact that the emission from both lines is optically thick and thermalized and emerges from dense regions of similar volumes. The CO, HCN and HCO+ luminosities suggest that the Cloverleaf is a composite AGN--starburst system, in agreement with the previous finding that about 20% of the total infrared luminosity in this system results from dust heated by star formation rather than heating by the AGN. We conclude that HCO+ is potentially a good tracer for dense molecular gas at high redshift.

preprint2006arXiv

Warm Dust and Spatially Variable PAH Emission in the Dwarf Starburst Galaxy NGC 1705

We present Spitzer observations of the dwarf starburst galaxy NGC 1705 obtained as part of SINGS. The galaxy morphology is very different shortward and longward of ~5 microns: short-wavelength imaging shows an underlying red stellar population, with the central super star cluster (SSC) dominating the luminosity; longer-wavelength data reveals warm dust emission arising from two off-nuclear regions offset by ~250 pc from the SSC. These regions show little extinction at optical wavelengths. The galaxy has a relatively low global dust mass (~2E5 solar masses, implying a global dust-to-gas mass ratio ~2--4 times lower than the Milky Way average). The off-nuclear dust emission appears to be powered by photons from the same stellar population responsible for the excitation of the observed H Alpha emission; these photons are unassociated with the SSC (though a contribution from embedded sources to the IR luminosity of the off-nuclear regions cannot be ruled out). Low-resolution IRS spectroscopy shows moderate-strength PAH emission in the 11.3 micron band in the eastern peak; no PAH emission is detected in the SSC or the western dust emission complex. There is significant diffuse 8 micron emission after scaling and subtracting shorter wavelength data; the spatially variable PAH emission strengths revealed by the IRS data suggest caution in the interpretation of diffuse 8 micron emission as arising from PAH carriers alone. The metallicity of NGC 1705 falls at the transition level of 35% solar found by Engelbracht and collaborators; the fact that a system at this metallicity shows spatially variable PAH emission demonstrates the complexity of interpreting diffuse 8 micron emission. A radio continuum non-detection, NGC 1705 deviates significantly from the canonical far-IR vs. radio correlation. (Abridged)

preprint2005arXiv

Spitzer Observations Of The Supergiant Shell Region In IC 2574

We present spatially resolved Spitzer imaging of the supergiant shell region of the M81 group dwarf galaxy IC 2574 obtained as part of the Spitzer Infrared Nearby Galaxies Survey. This region harbors one of the best nearby examples of a kinematically distinct HI shell, with an associated remnant stellar cluster; the shell is initiating sequential star formation as it interacts with the surrounding interstellar medium. This region dominates the infrared luminosity of IC 2574 and is spatially resolved in all Spitzer imaging bands. We study the differences in dust temperature as a function of local environment and compare local star formation rates as inferred from H Alpha and total infrared luminosities. We find that the strong H Alpha sources are associated with regions of warm dust; however, the most luminous infrared and H Alpha sources are not necessarilyco-spatial. The coolest dust is found in the regions farthest from the rim of the shell; these regions show the best agreement between star formation rates derived from H Alpha and from total infrared luminosities (although discrepancies at the factor of 3-4 level still exist). There is considerable variation in the radio-far infrared correlation in different regions surrounding the shell. The low dust content of the region may influence the scatter seen in these relations; these data demonstrate that the expanding shell is dramatically affecting its surroundings by triggering star formation and altering the dust temperature.

preprint2003arXiv

Molecular Gas in the Host Galaxy of a Quasar at Redshift z=6.42

Observations of the molecular gas phase in quasar host galaxies provide fundamental constraints on galaxy evolution at the highest redshifts. Molecular gas is the material out of which stars form; it can be traced by spectral line emission of carbon--monoxide (CO). To date, CO emission has been detected in more than a dozen quasar host galaxies with redshifts (z) larger 2, the record holder being at z=4.69. At these distances the CO lines are shifted to longer wavelengths, enabling their observation with sensitive radio and millimetre interferometers. Here we present the discovery of CO emission toward the quasar SDSS J114816.64+525150.3 (hereafter J1148+5251) at a redshift of z=6.42, when the universe was only 1/16 of its present age. This is the first detection of molecular gas at the end of cosmic reionization. The presence of large amounts of molecular gas (M(H_2)=2.2e10 M_sun) in an object at this time demonstrates that heavy element enriched molecular gas can be generated rapidly in the earliest galaxies.

preprint1999arXiv

The Violent Interstellar Medium of Nearby Dwarf Galaxies

High resolution HI observations of nearby dwarf galaxies (most of which are situated in the M 81 group at a distance of about 3.2 Mpc) reveal that their neutral interstellar medium (ISM) is dominated by hole-like features most of which are expanding. A comparison of the physical properties of these holes with the ones found in more massive spiral galaxies (such as M 31 and M 33) shows that they tend to reach much larger sizes in dwarf galaxies. This can be understood in terms of the galaxy's gravitational potential. The origin of these features is still a matter of debate. In general, young star forming regions (OB-associations) are held responsible for their formation. This picture, however, is not without its critics and other mechanism such as the infall of high velocity clouds, turbulent motions or even gamma ray bursters have been recently proposed. Here I will present one example of a supergiant shell in IC 2574 which corroborates the picture that OB associations are indeed creating these structures. This particular supergiant shell is currently the most promising case to study the effects of the combined effects of stellar winds and supernova-explosions which shape the neutral interstellar medium of (dwarf) galaxies.