Search

We explore non-redundant aperture masking to derive the size and shape of the ALBA synchrotron light source at optical wavelengths using synchrotron radiation interferometry. We show that non-redundant masks are required due to phase fluctuations arising within the experimental set-up. We also show, using closure phase, that the phase fluctuations are factorizable into element-based errors. We employ multiple masks, including 2, 3, 5, and 6 hole configurations. We develop a process for self-calibration of the element-based amplitudes (square root of flux through the aperture), which corrects for non-uniform illumination over the mask, in order to derive visibility coherences and phases, from which the source size and shape can be derived. We explore the optimal procedures to obtain the most reliable results with the 5-hole mask, based on the temporal scatter in measured coherences and closure phases. We find that the closure phases are very stable, and close to zero (within $2^o$). Through uv-modeling, we consider the noise properties of the experiment and conclude that our visibility measurements per frame are likely accurate to an rms scatter of $\sim 1\%$., Comment: 27 pages, 32 figures, 4 tables. For ALBA Synchrotron Light Source technical memo series

The Global mm-VLBI Array (GMVA) has demonstrated the ability to resolved what may be the general relativistic shadow of the supermassive black hole in M87 at 86 GHz, as well as delineate the inner jet to $\sim 1$~mas distance. We investigate the ability of the planned ngEHT, and the ngEHT + ngVLA at 85 GHz, to image such a nuclear 'ring', and the associated jet, using a constructed model based on the current estimate of the ring size, and a scaled version of the best VLBA image of the M87 jet at 43 GHz. While the resolution does not improve due to the limit set by the diameter of the Earth, the ngEHT alone should provide both a higher fidelity image of the ring on scales $\le 0.1$~mas, and a good image of a more extended jet to $\sim 1$~mas. Adding the ngVLA improves substantially the dynamic range (factor 3.5), as well as adds the ability to image structures on larger scales, in this case out to at least 5~mas, and potentially to much larger scales given the $\sim 10^5$ range in spatial scales covered by the ngVLA itself. Both arrays provide good image fidelity ($\le 0.1$), in the inner $\sim 1$~mas, but the ngEHT-only image does not reproduce the outer jet well, or at all, with fidelity values greater than unity. The combined array reproduces much of the outer jet with good fidelity ($\le 0.3$). Adding the ngVLA also decreases the susceptibility to antenna-based phase errors by a similar factor, and improves the ability for fringe fitting and subsequent phase and amplitude self-calibration. As for scales $< 100~\mu$as, ie. the ring itself, adding the ngVLA makes little change for very bright sources, where uniform weighting can be employed. But for faint sources, adding the ngVLA adds potentially an order-of magnitude sensitivity improvement (Issaoun et al. 2023)., Comment: 14 pages, 5 figures, ngVLA memo 118

We present results of dust continuum and [CII]$\,158\,{\rm \mu m}$ emission line observations of a remarkably UV-luminous ($M_{\rm UV}=-21.6$) galaxy at $z=10.603$: GN-z11. Using the Northern Extended Millimeter Array (NOEMA), observations have been carried out over multiple observing cycles. We achieved a high sensitivity resulting in a $\lambda_{\rm rest}=160\,{\rm \mu m}$ continuum $1\,\sigma$ depth of $13.0\,\rm{\mu Jy/beam}$ and a [CII] emission line $1\,\sigma$ sensitivity of $31\,\rm{mJy/beam\,km/s}$ using $50\,\rm{km/s}$ binning with a $\sim 2\,{\rm arcsec}$ synthesized beam. Neither dust continuum nor [CII]$\,158\,{\rm \mu m}$ line emission are detected at the expected frequency of $\nu_{\rm [CII]} = 163.791\,\rm{GHz}$ and the sky location of GN-z11. The upper limits show that GN-z11 is neither luminous in $L_{\rm IR}$ nor $L_{\rm [CII]}$, with a dust mass $3\,\sigma$ limit of ${\rm log}(M_{\rm dust}/{\rm M_{\odot}}) < 6.5-6.9$ and with a [CII] based molecular gas mass $3\,\sigma$ limit of ${\rm log}(M_{\rm mol,[CII]}/{\rm M_{\odot}}) < 9.3$. Together with radiative transfer calculations, we also investigated the possible cause of the dust poor nature of the GN-z11 showed by the blue color in the UV continuum of GN-z11 ($\beta_{\rm UV}=-2.4$), and found that $\gtrsim3\times$ deeper observations are crucial to study dust production at very high-redshift. Nevertheless, our observations show the crucial role of deep mm/submm observations of very high redshift galaxies to constrain multiple phases in the interstellar medium., Comment: submitted to MNRAS, 7 pages, 4 figures, 1 table

We develop a new process of image plane self-calibration for interferometric imaging data. The process is based on Shape-Orientation-Size (SOS) conservation for the principal triangle in an image generated from the three fringes made from a triad of receiving elements, in situations where interferometric phase errors can be factorized into element-based terms. The basis of the SOS conservation principle is that, for a 3-element array, the only possible image corruption due to an element-based phase screen is a tilt of the aperture plane, leading to a shift in the image plane. Thus, an image made from any 3-element interferometer represents a true image of the source brightness, modulo an unknown translation. Image plane self-calibration entails deriving the unknown translations for each triad image via cross-correlation of the observed triad image with a model image of the source brightness. After correcting for these independent shifts, and summing the aligned triad images, a good image of the source brightness is generated from the full array, recovering source structure at diffraction-limited resolution. The process is iterative, using improved source models based on previous iterations. We demonstrate the technique in the high signal-to-noise context, and include a configuration based on radio astronomical facilities, and simple models of double sources. We show that the process converges for the simple models considered, although convergence is slower than for aperture-plane self-calibration for large-$N$ arrays. As currently implemented, the process is most relevant for arrays with a small number of elements. More generally, the technique provides geometric insight into closure phase and the self-calibration process. The technique is generalizable to non-astronomical interferometric imaging applications across the electromagnetic spectrum., Comment: 17 pages 4 figures. To appear in Journal Optical Society America - A

We present simulations of the capabilities of the ngVLA to image at $\sim 0.75$ kpc resolution ($0.085"$), molecular line emission from star forming disk galaxies at high redshift. The results are compared to the current capabilities of ALMA. ALMA can detect the integrated emission, and determine the velocity gradient and size across the brighter emission regions of the galaxy. The ngVLA is a factor $\sim 6$ more sensitive at the adopted spatial and velocity resolution. This sensitivity is needed to recover the detailed column density distribution, velocity field, and velocity dispersion at full resolution. The ngVLA will enable detailed analysis of spectral line profiles at $0.75$~kpc resolution, even in relatively faint regions. The ngVLA will trace the rotation curves to large radii, and recover sub-structure in the disks, such as clumps, spiral arms, bars, and rings. Detection of these features is crucial in order to assess how cold gas precipitates the formation of stars at high redshift., Comment: ngVLA memo series 103 https://ngvla.nrao.edu/page/memos

We present the discovery of a 4 kpc molecular gas lane in the Cygnus A host galaxy, using ALMA CO 2-1 observations. The gas lane is oriented roughly perpendicular to the projected radio jet axis. The CO emission generally follows the clumpy dust lanes seen in HST I-band images. The total molecular gas mass is $30\times 10^8$ M$_\odot$ for Milky Way type clouds, and $3.6 \times 10^8$ M$_\odot$ for starburst conditions. There is a velocity change from the northern to southern CO peaks of about $\pm 175$~km~s$^{-1}$, and an apparently smooth velocity gradient between the peaks, although the emission in the central region is weak. In the inner $\sim 0.5"$ projected distance from the radio core, comparison of the CO velocities to those observed for H$_2$ 2.1218 $\mu$m emission shows higher velocities for the vibrationally excited warm molecular gas than the cooler CO 2-1 line emitting gas at similar projected radii. A possible explanation for these different projected velocities at a given radius is that the cooler CO gas is distributed in a clumpy ring at radius $\sim 1.5"$ to $2"$, while the warm H$_2$ 2.12$\mu$m emitting gas is interior to this ring. Of course, the current data cannot rule-out a clumpy, amorphous molecular gas distribution linearly distributed perpendicular to the radio jet axis. We consider surface brightness properties on scales down to $\sim 265$~pc, and discuss the Cygnus A results in the context of other radio galaxies with CO emission., Comment: 5 Figures, 23 pages, 2 tables, accepted for publication in the ApJ

We review the evolution of the cosmic average molecular gas density to large look-back times, using observations of rotational transitions of CO. Molecular gas is the fuel for star formation in galaxies. Deep searches for CO emission from distant galaxies have delineated the density of molecular gas back to $z \sim 5$, or within 1~Gyr of the Big Bang. The results show a rise and fall in the gas density that parallels, and likely drives, the rise and fall of the cosmic star formation rate density. We present the potential for the next generation Very Large Array to image the distribution and dynamics of the molecular gas in early galaxies, and to make a precise measurement of the dense gas history of the Universe., Comment: 5 pages. 4 figures. to appear in Radio Science

We investigate the imaging performance of an interferometric array in the case of wide field, high resolution, narrow band, snapshot imaging. We find that, when uv-cell sizes are sufficiently small (ie. image sizes are sufficiently large), each instantaneous visibility record is gridded into its own uv-cell. This holds even for dense arrays, like the core of the next generation VLA. In this particular, application, Uniform weighting of the gridded visibilities approaches Natural weighting, with its often deleterious consequences on the resulting synthesized beam. For a core-dominated array, we show that the resulting image noise is highly correlated on scales comparable to the spatial frequencies of the core baselines. In general, this study accentuates the fact that, for imaging applications that require high resolution (Plains array and greater), many of the core antennas can be employed as a separate subarray for low resolution science, without sacrificing the quality of the high resolution science., Comment: 18 pages; Next Generation VLA Memo No. 78; https://ngvla.nrao.edu/page/memos#gen-memo

We present the first direct imaging of what may be the thick torus in the active galactic nucleus (AGN) of the archetype powerful radio galaxy Cygnus A, using the Jansky Very Large Array (VLA) at 18 GHz to 48 GHz, with a resolution down to 45 mas. Such a torus has long been a key component of AGN models, but direct imaging on the relevant physical scales in sources of extreme (quasar-like) luminosities, remains scarce. An elongated structure, perpendicular to the radio jets and centered on the core, is well resolved, with a full length of $0.48"$ (528 pc), and a full width of $0.26"$ (286 pc). The radio emission spectrum is consistent with optically thin free-free emission. We present a toy model of a flaring torus, with a half-opening angle for the poloidal region of $62^o$. The radio jets are oriented along the poles. The observations require a clumpy gas distribution, with the free-free emission dominated by clumps with densities $\ge 4000$ cm$^{-3}$, Comment: Accepted to Astrophysical Journal Letters, 10 pages, 7 figures

We consider the capabilities of ALMA and the ngVLA to detect and image the[CII] 158\,$\mu$m line from galaxies into the cosmic `dark ages' ($z \sim 10$ to 20). The [CII] line may prove to be a powerful tool in determining spectroscopic redshifts, and galaxy dynamics, for the first galaxies. In 40\,hr, ALMA has the sensitivity to detect the integrated [CII] line emission from a moderate metallicity, active star-forming galaxy [$Z_A = 0.2\,Z_{\odot}$; star formation rate (SFR)= 5\,$M_\odot$\,yr$^{-1}$], at $z = 10$ at a significance of 6$\sigma$. The ngVLA will detect the integrated [CII] line emission from a Milky-Way like star formation rate galaxy ($Z_{A} = 0.2\,Z_{\odot}$, SFR = 1\,$M_\odot$\,yr$^{-1}$), at $z = 15$ at a significance of 6$\sigma$. Imaging simulations show that the ngVLA can determine rotation dynamics for active star-forming galaxies at $z \sim 15$, if they exist. The [CII] detection rate in blind surveys will be slow (of order unity per 40\,hr pointing., Comment: 8 pages, 3 figures; to be published in the ASP Monograph Series, "Science with a Next-Generation VLA", ed. E. J. Murphy (ASP, San Francisco, CA). arXiv admin note: substantial text overlap with arXiv:1709.02001

We perform simulations of the capabilities of the next generation Very Large Array to image stellar radio photospheres. For very large (in angle) stars, such as red supergiants within a few hundred parsecs, good imaging fidelity results can be obtained on radio photospheric structures at 38 GHz employing standard techniques, such as disk model fitting and subtraction, with hundreds of resolution elements over the star, even with just the ngVLA-classic baselines to 1000 km. Using the ngVLA Rev B plus long baseline configuration (with baselines out to 9000 km, August 2018), we find for main sequence stars within $\sim$ 10 pc, the photospheres can be easily resolved at 85 GHz, with accurate measures of the mean brightness and size, and possibly imaging large surface structures, as might occur on e.g., active M dwarf stars. For more distant main sequence stars, we find that measurements of sizes and brightnesses can be made using disk model fitting to the u,v-data down to stellar diameters $\sim$ 0.4 mas in a few hours. This size would include M0 V stars to a distance of 15 pc, A0 V stars to 60 pc, and Red Giants to 2.4 kpc. Based on the Hipparcos catalog, we estimate that there are at least 10,000 stars that will be resolved by the ngVLA. While the vast majority of these (95\%) are giants or supergiants, there are still over 500 main sequence stars that can be resolved, with $\sim$ 50 to 150 in each spectral type (besides O stars). Note that these are lower limits, since radio photospheres can be larger than optical, and the Hipparcos catalog might not be complete. Our initial look into the Gaia catalog suggests these numbers might be pessimistic by a factor few., Comment: 12 pages 3 figures 1 table; accepted for ASP 'Science with the Next Generation Very Large Array'

We perform simulations of the capabilities of the next generation Very Large Array in the context of imaging low order CO emission from typical high redshift star forming galaxies at ~ 1 kpc resolution. We adopt as a spatial and dynamical template the CO 1-0 emission from M 51, scaled accordingly for redshift, transition, and total gas mass. The molecular gas masses investigated are factors of 1.4, 3.5, and 12.5 larger that of M 51, at z = 0.5, 2, and 4.2, respectively. The z = 2 galaxy gas mass is comparable to the lowest mass galaxies currently being discovered in the deepest ALMA and NOEMA cosmological CO line surveys, corresponding to galaxies with star formation rates ~ 10 to 100 $M_\odot$ yr$^{-1}$. The ngVLA will perform quality imaging at 1kpc resolution of the gas distribution and dynamics over this disk. We recover the overall rotation curve, galaxy orientation properties, and molecular ISM internal velocity dispersion. The model at z = 4.2 corresponds to a massive star forming main sequence disk (SFR ~ 130 $M_\odot$ yr$^{-1}$). The ngVLA can obtain 1kpc resolution images of such a system in a reasonable integration time, and recover the basic galaxy orientation parameters, and, asymptotically, the maximum rotation velocity. We compare the ngVLA results with capabilities of ALMA and the Jansky VLA. ALMA and the VLA can detect the integrated low order CO emission from these galaxies, but lack the sensitivity to perform the high resolution imaging to recover the dynamics at 1kpc scales. To do so would require of order 1000 hrs per galaxy with these current facilities. We investigate a 'minimal' ngVLA configuration, removing the longest baselines and much of the very compact core, and find good imaging can still be performed at 1 kpc resolution., Comment: 9 pages 9 figures 2 tables to appear in Science with the next generation Very Large Array ASP

New massively redundant low frequency arrays allow for a novel investigation of closure relations in interferometry. We employ commissioning data from the Hydrogen Epoch of Reionization Array to investigate closure quantities in this densely packed grid array of 14m antennas operating at 100 MHz to 200 MHz. We investigate techniques that utilize closure phase spectra for redundant triads to estimate departures from redundancy for redundant baseline visibilities. We find a median absolute deviation from redundancy in closure phase across the observed frequency range of about 4.5deg. This value translates into a non-redundancy per visibility phase of about 2.6deg, using prototype electronics. The median absolute deviations from redundancy decrease with longer baselines. We show that closure phase spectra can be used to identify ill-behaved antennas in the array, independent of calibration. We investigate the temporal behavior of closure spectra. The Allan variance increases after a one minute stride time, due to passage of the sky through the primary beam of the transit telescope. However, the closure spectra repeat to well within the noise per measurement at corresponding local sidereal times (LST) from day to day. In future papers in this series we will develop the technique of using closure phase spectra in the search for the HI 21cm signal from cosmic reionization., Comment: 32 pages. 11 figures. Accepted to Radio Science

We use new and archival Chandra observations of Cygnus A, totalling $\sim$1.9 Ms, to investigate the distribution and temperature structure of gas lying within the projected extent of the cocoon shock and exhibiting a rib-like structure. We confirm that the X-rays are dominated by thermal emission with an average temperature of around 4 keV, and have discovered an asymmetry in the temperature gradient, with the southwestern part of the gas cooler than the rest by up to 2 keV. Pressure estimates suggest that the gas is a coherent structure of single origin located inside the cocoon, with a mass of roughly $2\times10^{10} M_{\odot}$. We conclude that the gas is debris resulting from disintegration of the cool core of the Cygnus A cluster after the passage of the jet during the early stages of the current epoch of activity. The 4 keV gas now lies on the central inside surface of the hotter cocoon rim. The temperature gradient could result from an offset between the centre of the cluster core and the Cygnus A host galaxy at the switch-on of current radio activity., Comment: 15 pages, 12 figures, accepted by MNRAS

Ultra-steep spectrum (USS) radio sources are good tracers of powerful radio galaxies at $z > 2$. Identification of even a single bright radio galaxy at $z > 6$ can be used to detect redshifted 21cm absorption due to neutral hydrogen in the intervening IGM. Here we describe a new sample of high-redshift radio galaxy (HzRG) candidates constructed from the TGSS ADR1 survey at 150 MHz. We employ USS selection ($\alpha \le -1.3$) in $\sim10000$ square degrees, in combination with strict size selection and non-detections in all-sky optical and infrared surveys. We apply flux density cuts that probe a unique parameter space in flux density ($50 < S_{\textrm{150}} < 200$ mJy) to build a sample of 32 HzRG candidates. Follow-up Karl G. Jansky Very Large Array (VLA) observations at 1.4 GHz with an average beam size of $1.3$ arcseconds ($"$) revealed $\sim 48\%$ of sources to have a single radio component. P-band (370 MHz) imaging of 17 of these sources revealed a flattening radio SED for ten sources at low frequencies, which is expected from compact HzRGs. Two of our sources lie in fields where deeper multi-wavelength photometry and ancillary radio data are available and for one of these we find a best-fit photo-z of $4.8 \pm 2.0$. The other source has $z_{\textrm{phot}}=1.4 \pm 0.1$ and a small angular size ($3.7"$), which could be associated with an obscured star forming galaxy or with a "dead" elliptical. One USS radio source not part of the HzRG sample but observed with the VLA nonetheless is revealed to be a candidate giant radio galaxy with a host galaxy photo-z of $1.8\pm0.5$, indicating a size of 875 kpc., Comment: 17 pages, 13 figures, accepted for publication in MNRAS

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

Until recently, determining the rotational properties of galaxies in the early universe (z>4, Universe age <1.5Gyr) was impractical, with the exception of a few strongly lensed systems. Combining the high resolution and sensitivity of ALMA at (sub-) millimeter wavelengths with the typically high strength of the [CII] 158um emission line from galaxies and long-developed dynamical modeling tools raises the possibility of characterizing the gas dynamics in both extreme starburst galaxies and normal star forming disk galaxies at z~4-7. Using a procedure centered around GIPSY's ROTCUR task, we have fit tilted ring models to some of the best available ALMA [CII] data of a small set of galaxies: the MS galaxies HZ9 & HZ10, the Damped Lyman-alpha Absorber (DLA) host galaxy ALMA J0817+1351, the submm galaxies AzTEC/C159 and COSMOS J1000+0234, and the quasar host galaxy ULAS J1319+0950. This procedure directly derives rotation curves and dynamical masses as functions of radius for each object. In one case, we present evidence for a dark matter halo of O(10^11) solar masses. We present an analysis of the possible velocity dispersions of AzTEC/C159 and ULAS J1319+0950 based on matching simulated observations to the integrated [CII] line profiles. Finally, we test the effects of observation resolution and sensitivity on our results. While the conclusions remain limited at the resolution and signal-to-noise ratios of these observations, the results demonstrate the viability of the modeling tools at high redshift, and the exciting potential for detailed dynamical analysis of the earliest galaxies, as ALMA achieves full observational capabilities., Comment: 20 pages, 6 figures, Accepted for publication in ApJ

We consider the capabilities of current and future large facilities operating at 2\,mm to 3\,mm wavelength to detect and image the [CII] 158\,$\mu$m line from galaxies into the cosmic "dark ages" ($z \sim 10$ to 20). The [CII] line may prove to be a powerful tool in determining spectroscopic redshifts, and galaxy dynamics, for the first galaxies. We emphasize that the nature, and even existence, of such extreme redshift galaxies, remains at the frontier of open questions in galaxy formation. In 40\,hr, ALMA has the sensitivity to detect the integrated [CII] line emission from a moderate metallicity, active star-forming galaxy [$Z_A = 0.2\,Z_{\odot}$; star formation rate (SFR) = 5\,$M_\odot$\,yr$^{-1}$], at $z = 10$ at a significance of 6$\sigma$. The next-generation Very Large Array (ngVLA) will detect the integrated [CII] line emission from a Milky-Way like star formation rate galaxy ($Z_{A} = 0.2\,Z_{\odot}$, SFR = 1\,$M_\odot$\,yr$^{-1}$), at $z = 15$ at a significance of 6$\sigma$. Imaging simulations show that the ngVLA can determine rotation dynamics for active star-forming galaxies at $z \sim 15$, if they exist. Based on our very limited knowledge of the extreme redshift Universe, we calculate the count rate in blind, volumetric surveys for [CII] emission at $z \sim 10$ to 20. The detection rates in blind surveys will be slow (of order unity per 40\,hr pointing). However, the observations are well suited to commensal searches. We compare [CII] with the [OIII] 88$\mu$m line, and other ancillary information in high $z$ galaxies that would aid these studies., Comment: 11pages, 8 figures, Accepted for the Astrophysical Journal

We present ALMA observations of the dust continuum and [C II] 158um line emission from the z=6.0695 Lyman Break Galaxy WMH5. These observations at 0.3" spatial resolution show a compact (~3kpc) main galaxy in dust and [C II] emission, with a 'tail' of emission extending to the east by about 5kpc (in projection). The [C II] tail is comprised predominantly of two distinct sub-components in velocity, separated from the core by ~100 and 250km/s, with narrow intrinsic widths of about 80km/s, which we call 'sub-galaxies'. The sub-galaxies themselves are extended east-west by about 3kpc in individual channel images. The [C II] tail joins smoothly into the main galaxy velocity field. The [C II] line to continuum ratios are comparable for the main and sub-galaxy positions, within a factor 2. In addition, these ratios are comparable to z~5.5 LBGs. We conjecture that the WMH5 system represents the early formation of a galaxy through the accretion of smaller satellite galaxies, embedded in a smoother gas distribution, along a possibly filamentary structure. The results are consistent with current cosmological simulations of early galaxy formation, and support the idea of very early enrichment with dust and heavy elements of the accreting material., Comment: 10 pages, 9 figures, Accepted for publication in ApJ

Large-scale vortices in protoplanetary disks are thought to form and survive for long periods of time. Hence, they can significantly change the global disk evolution and particularly the distribution of the solid particles embedded in the gas, possibly explaining asymmetries and dust concentrations recently observed at sub-millimeter and millimeter wavelengths. We investigate the spatial distribution of dust grains using a simple model of protoplanetary disk hosted by a giant gaseous vortex. We explore the dependence of the results on grain size and deduce possible consequences and predictions for observations of the dust thermal emission at sub-millimeter and millimeter wavelengths. Global 2D simulations with a bi-fluid code are used to follow the evolution of a single population of solid particles aerodynamically coupled to the gas. Possible observational signatures of the dust thermal emission are obtained using simulators of ALMA and ngVLA observations. We find that a giant vortex not only captures dust grains with Stokes number St < 1 but can also affect the distribution of larger grains (with St '~' 1) carving a gap associated to a ring composed of incompletely trapped particles. The results are presented for different particle size and associated to their possible signatures in disk observations. Gap clearing in the dust spatial distribution could be due to the interaction with a giant gaseous vortex and their associated spiral waves, without the gravitational assistance of a planet. Hence, strong dust concentrations at short sub-mm wavelengths associated with a gap and an irregular ring at longer mm and cm wavelengths could indicate the presence of an unseen gaseous vortex., Comment: 11 pages, 11 figures, accepted for publication in A&A

We make use of the deep Karl G. Jansky Very Large Array (VLA) COSMOS radio observations at 3 GHz to infer radio luminosity functions of star-forming galaxies up to redshifts of z~5 based on approximately 6000 detections with reliable optical counterparts. This is currently the largest radio-selected sample available out to z~5 across an area of 2 square degrees with a sensitivity of rms=2.3 ujy/beam. By fixing the faint and bright end shape of the radio luminosity function to the local values, we find a strong redshift trend that can be fitted with a pure luminosity evolution L~(1+z)^{(3.16 +- 0.2)-(0.32 +- 0.07) z}. We estimate star formation rates (SFRs) from our radio luminosities using an infrared (IR)-radio correlation that is redshift dependent. By integrating the parametric fits of the evolved luminosity function we calculate the cosmic SFR density (SFRD) history since z~5. Our data suggest that the SFRD history peaks between 21000 Msol/yr) contribute an additional <2% in the entire observed redshift range. We find evidence of a potential underestimation of SFRD based on ultraviolet (UV) rest-frame observations of Lyman break galaxies (LBGs) at high redshifts (z>4) on the order of 15-20%, owing to appreciable star formation in highly dust-obscured galaxies, which might remain undetected in such UV observations., Comment: 18 pages, 8 figures, 2 tables, accepted for publication in A&A

We present the VLA-COSMOS 3 GHz Large Project based on 384 hours of observations with the Karl G. Jansky Very Large Array (VLA) at 3 GHz (10 cm) toward the two square degree Cosmic Evolution Survey (COSMOS) field. The final mosaic reaches a median rms of 2.3 uJy/beam over the two square degrees at an angular resolution of 0.75". To fully account for the spectral shape and resolution variations across the broad (2 GHz) band, we image all data with a multiscale, multifrequency synthesis algorithm. We present a catalog of 10,830 radio sources down to 5 sigma, out of which 67 are combined from multiple components. Comparing the positions of our 3 GHz sources with those from the Very Long Baseline Array (VLBA)-COSMOS survey, we estimate that the astrometry is accurate to 0.01" at the bright end (signal-to-noise ratio, S/N_3GHz > 20). Survival analysis on our data combined with the VLA-COSMOS 1.4~GHz Joint Project catalog yields an expected median radio spectral index of alpha=-0.7. We compute completeness corrections via Monte Carlo simulations to derive the corrected 3 GHz source counts. Our counts are in agreement with previously derived 3 GHz counts based on single-pointing (0.087 square degrees) VLA data. In summary, the VLA-COSMOS 3 GHz Large Project simultaneously provides the largest and deepest radio continuum survey at high (0.75") angular resolution to date, bridging the gap between last-generation and next-generation surveys., Comment: 19 pages, 18 figures, 5 tables, accepted for publication in A&A

The largest galaxies in the Universe reside in galaxy clusters. Using sensitive observations of carbon-monoxide, we show that the Spiderweb Galaxy -a massive galaxy in a distant protocluster- is forming from a large reservoir of molecular gas. Most of this molecular gas lies between the protocluster galaxies and has low velocity dispersion, indicating that it is part of an enriched inter-galactic medium. This may constitute the reservoir of gas that fuels the widespread star formation seen in earlier ultraviolet observations of the Spiderweb Galaxy. Our results support the notion that giant galaxies in clusters formed from extended regions of recycled gas at high redshift., Comment: Published in Science on 2 Dec 2016 (accepted 21 Oct 2016), 20 pages, 5 figures, 3 tables. Contains minor copy-editing differences with published version

We present ALMA observations of cold dust and molecular gas in four high-luminosity, heavily reddened (A$_{\rm{V}} \sim 2.5-6$ mag) Type 1 quasars at $z\sim2.5$ with virial M$_{\rm{BH}} \sim 10^{10}$M$_\odot$, to test whether dusty, massive quasars represent the evolutionary link between submillimetre bright galaxies (SMGs) and unobscured quasars. All four quasars are detected in both the dust continuum and in the $^{12}$CO(3-2) line. The mean dust mass is 6$\times$10$^{8}$M$_\odot$ assuming a typical high redshift quasar spectral energy distribution (T=41K, $\beta$=1.95 or T=47K, $\beta$=1.6). The implied star formation rates are very high - $\gtrsim$1000 M$_\odot$ yr$^{-1}$ in all cases. Gas masses estimated from the CO line luminosities cover $\sim$1-5$\times10^{10}$($\alpha_{\rm{CO}} / 0.8$)M$_\odot$ and the gas depletion timescales are very short - $\sim5-20$Myr. A range of gas-to-dust ratios is observed in the sample. We resolve the molecular gas in one quasar - ULASJ2315$+$0143 ($z=2.561$) - which shows a strong velocity gradient over $\sim$20 kpc. The velocity field is consistent with a rotationally supported gas disk but other scenarios, e.g. mergers, cannot be ruled out at the current resolution of these data. In another quasar - ULASJ1234+0907 ($z=2.503$) - we detected molecular line emission from two millimetre bright galaxies within 200 kpc of the quasar, suggesting that this quasar resides in a significant over-density. The high detection rate of both cold dust and molecular gas in these sources, suggests that reddened quasars could correspond to an early phase in massive galaxy formation associated with large gas reservoirs and significant star formation., Comment: MNRAS Accepted; Very minor changes to match accepted version

We present Karl G Jansky Very Large Array (VLA) observations of CO(2-1) line emission and rest-frame 250GHz continuum emission of the Hyper-Luminous IR Galaxies (HyLIRGs) BRI1202-0725 (z=4.69) and BRI1335-0417 (z=4.41), with an angular resolution as high as 0.15". Our low order CO observations delineate the cool molecular gas, the fuel for star formation in the systems, in unprecedented detail. For BRI1202-0725, line emission is seen from both extreme starburst galaxies: the quasar host and the optically obscured submm galaxy (SMG), in addition to one of the Lyman-alpha emitting galaxies in the group. Line emission from the SMG shows an east-west extension of about 0.6". For Lyalpha-2, the CO emission is detected at the same velocity as [CII] and [NII], indicating a total gas mass ~4.0*10^10 solar masses. The CO emission from BRI1335-0417 peaks at the nominal quasar position, with a prominent northern extension (~1", a possible tidal feature). The gas depletion timescales are ~10^7 years for the three HyLIRGs, consistent with extreme starbursts, while that of Lyalpha-2 may be consistent with main sequence galaxies. We interpret these sources as major star formation episodes in the formation of massive galaxies and supermassive black holes (SMBHs) via gas rich mergers in the early Universe., Comment: 28 pages, 10 figures, Accepted for publication in ApJ

We present a search for [CII] line and dust continuum emission from optical dropout galaxies at $z>6$ using ASPECS, our ALMA Spectroscopic Survey in the Hubble Ultra-Deep Field (UDF). Our observations, which cover the frequency range $212-272$ GHz, encompass approximately the range $6$4.5 $\sigma$, two of which correspond to blind detections with no optical counterparts. At this significance level, our statistical analysis shows that about 60\% of our candidates are expected to be spurious. For one of our blindly selected [CII] line candidates, we tentatively detect the CO(6-5) line in our parallel 3-mm line scan. None of the line candidates are individually detected in the 1.2 mm continuum. A stack of all [CII] candidates results in a tentative detection with $S_{1.2mm}=14\pm5\mu$Jy. This implies a dust-obscured star formation rate (SFR) of $(3\pm1)$ M$_\odot$ yr$^{-1}$. We find that the two highest--SFR objects have candidate [CII] lines with luminosities that are consistent with the low-redshift $L_{\rm [CII]}$ vs. SFR relation. The other candidates have significantly higher [CII] luminosities than expected from their UV--based SFR. At the current sensitivity it is unclear whether the majority of these sources are intrinsically bright [CII] emitters, or spurious sources. If only one of our line candidates was real (a scenario greatly favored by our statistical analysis), we find a source density for [CII] emitters at $6

We present direct estimates of the mean sky brightness temperature in observing bands around 99GHz and 242GHz due to line emission from distant galaxies. These values are calculated from the summed line emission observed in a blind, deep survey for specrtal line emission from high redshift galaxies using ALMA (the 'ASPECS' survey). In the 99 GHz band, the mean brightness will be dominated by rotational transitions of CO from intermediate and high redshift galaxies. In the 242GHz band, the emission could be a combination of higher order CO lines, and possibly [CII] 158$\mu$m line emission from very high redshift galaxies ($z \sim 6$ to 7). The mean line surface brightness is a quantity that is relevant to measurements of spectral distortions of the cosmic microwave background, and as a potential tool for studying large-scale structures in the early Universe using intensity mapping. While the cosmic volume and the number of detections are admittedly small, this pilot survey provides a direct measure of the mean line surface brightness, independent of conversion factors, excitation, or other galaxy formation model assumptions. The mean surface brightness in the 99GHZ band is: $T_B = 0.94\pm 0.09$ $\mu$K. In the 242GHz band, the mean brightness is: $T_B = 0.55\pm 0.033$ $\mu$K. These should be interpreted as lower limits on the average sky signal, since we only include lines detected individually in the blind survey, while in a low resolution intensity mapping experiment, there will also be the summed contribution from lower luminosity galaxies that cannot be detected individually in the current blind survey., Comment: To appear the Astrophysical Journal, incorporating the first and second round of referee comments. Revised version 2 fixed typo in brightness units, equation, and text, and revised acknowledgements

We present an analysis of a deep (1$\sigma$=13 $\mu$Jy) cosmological 1.2-mm continuum map based on ASPECS, the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field. In the 1 arcmin$^2$ covered by ASPECS we detect nine sources at $>3.5\sigma$ significance at 1.2-mm. Our ALMA--selected sample has a median redshift of $z=1.6\pm0.4$, with only one galaxy detected at z$>$2 within the survey area. This value is significantly lower than that found in millimeter samples selected at a higher flux density cut-off and similar frequencies. Most galaxies have specific star formation rates similar to that of main sequence galaxies at the same epoch, and we find median values of stellar mass and star formation rates of $4.0\times10^{10}\ M_\odot$ and $\sim40~M_\odot$ yr$^{-1}$, respectively. Using the dust emission as a tracer for the ISM mass, we derive depletion times that are typically longer than 300 Myr, and we find molecular gas fractions ranging from $\sim$0.1 to 1.0. As noted by previous studies, these values are lower than using CO--based ISM estimates by a factor $\sim$2. The 1\,mm number counts (corrected for fidelity and completeness) are in agreement with previous studies that were typically restricted to brighter sources. With our individual detections only, we recover $55\pm4\%$ of the extragalactic background light (EBL) at 1.2 mm measured by the Planck satellite, and we recover $80\pm7\%$ of this EBL if we include the bright end of the number counts and additional detections from stacking. The stacked contribution is dominated by galaxies at $z\sim1-2$, with stellar masses of (1-3)$\times$10$^{10}$ M$_\odot$. For the first time, we are able to characterize the population of galaxies that dominate the EBL at 1.2 mm., Comment: Accepted for publication on the Astrophysical Journal

We present observations with the Karl G. Jansky Very Large Array (VLA) at 3 GHz toward a sub-field of the XXL-North 25 deg$^{2}$ field targeting the first supercluster discovered in the XXL Survey. The structure has been found at a spectroscopic redshift of 0.43 and extending over $0.35\times0.1$ deg$^{2}$ on the sky. We present the 3 GHz VLA radio continuum observations, the radio mosaic and radio source catalogue, and, perform an analysis of the supercluster in the optical and radio regimes using photometric redshifts from the CFHTLS survey and our new VLA data. Our 3 GHz radio mosaic has a resolution of $3.2\times1.9$ arcsec$^2$, and covers an area of $41\times41$ arcmin$^2$ with rms noise level below $\sim20 \mu$Jy beam$^{-1}$ ($\sim11\mu$Jy beam$^{-1}$ in central $15\times15$ arcmin$^2$). We extract a catalogue of 155 radio sources $\geq6\sigma$, 8 of which are multicomponent sources, and 123 ($79\%$) can be associated with optical sources in the CFHTLS W1 catalogue. Applying Voronoi tessellation analysis (VTA) in the area around the X-ray identified supercluster we identify a total of 17 overdensities at $z=0.35-0.50$, 7 of which are associated with clusters detected in the $XMM-Newton$ XXL data. We find a mean photometric redshift of 0.43 for our overdensities, consistent with the spectroscopic redshifts of 7 X-ray detected clusters. The full VTA-identified structure extends over $\sim0.6\times0.2$ deg$^2$ on the sky, ($\sim12\times4$ Mpc$^{2}$ at $z=0.43$). We associate 8 radio sources with potential group/cluster member galaxies none of which are large radio galaxies. The spatial distribution of the red and blue potential group member galaxies suggests that the clusters are dynamically young and not virialised yet, as expected for hierarchical structure growth in a $\Lambda$CDM universe. Further spectroscopic data are required to analyze the dynamical state of the groups.

We investigate the use of closure phase as a method to detect the HI 21cm signal from the neutral IGM during cosmic reionzation. Closure quantities have the unique advantage of being independent of antenna-based calibration terms. We employ realistic, large area sky models from Sims et al. (2016). These include an estimate of the HI 21cm signal generated using 21cm FAST, plus continuum models of both the diffuse Galactic synchrotron emission and the extragalactic point sources. We employ the CASA simulator and adopt the Dillon-Parsons HERA configuration to generate a uv measurement set. We then use AIPS to calculate the closure phases as a function of frequency ('closure spectra'), and python scripts for subsequent analysis. We find that the closure spectra for the HI signal show dramatic structure in frequency, and based on thermal noise alone, the redundant HERA-331 array should detect these fluctuations easily. Comparatively, the frequency structure in the continuum closure spectra is much smoother than that seen in the HI closure spectra. Unfortunately, when the line and continuum signals are combined, the continuum dominates the visibilities at the level of 10^3 to 10^4, and the line signal is lost. We have investigated fitting and removing smooth curves in frequency to the line plus continuum closure spectra, and find that the continuum itself shows enough structure in frequency in the closure spectra to preclude separation of the continuum and line based on such a process. We have also considered the subtraction of the continuum from the visibilities using a sky model, prior to calculation of the closure spectra. TRUNCATED., Comment: 14 pages. HERA Memo Number 13 located at http://reionization.org/science/memos/ reload -- changed one word

We summarize the design, capabilities, and some of the priority science goals of a next generation Very Large Array (ngVLA). The ngVLA is an interferometric array with 10x larger effective collecting area and 10x higher spatial resolution than the current VLA and the Atacama Large Millimeter Array (ALMA), optimized for operation in the wavelength range 0.3cm to 3cm. The ngVLA opens a new window on the Universe through ultra-sensitive imaging of thermal line and continuum emission down to milliarcecond resolution, as well as unprecedented broad band continuum polarimetric imaging of non-thermal processes. The continuum resolution will reach 9mas at 1cm, with a brightness temperature sensitivity of 6K in 1 hour. For spectral lines, the array at 1" resolution will reach 0.3K surface brightness sensitivity at 1cm and 10 km/s spectral resolution in 1 hour. These capabilities are the only means with which to answer a broad range of critical scientific questions in modern astronomy, including direct imaging of planet formation in the terrestrial-zone, studies of dust-obscured star formation and the cosmic baryon cycle down to pc-scales out to the Virgo cluster, making a cosmic census of the molecular gas which fuels star formation back to first light and cosmic reionization, and novel techniques for exploring temporal phenomena from milliseconds to years. The ngVLA is optimized for observations at wavelengths between the superb performance of ALMA at submm wavelengths, and the future SKA1 at few centimeter and longer wavelengths. This memo introduces the project. The science capabilities are outlined in a parallel series of white papers. We emphasize that this initial set of science goals are simply a starting point for the project. We invite comment on these programs, as well as new ideas, through our public forum link on the ngVLA web page https://science.nrao.edu/futures/ngvla, Comment: 6 pages, 5 figures. See memo series for science working group reports as well as public input forum on ngVLA project page: https://science.nrao.edu/futures/ngvla

In the low-redshift Universe, the most powerful radio sources are often associated with gas-rich galaxy mergers or interactions. We here present evidence for an advanced, gas-rich (`wet') merger associated with a powerful radio galaxy at a redshift of z~2. This radio galaxy, MRC 0152-209, is the most infrared-luminous high-redshift radio galaxy known in the southern hemisphere. Using the Australia Telescope Compact Array, we obtained high-resolution CO(1-0) data of cold molecular gas, which we complement with HST/WFPC2 imaging and WHT long-slit spectroscopy. We find that, while roughly M(H2) ~ 2 x 10$^{10}$ M$_{\odot}$ of molecular gas coincides with the central host galaxy, another M(H2) ~ 3 x 10$^{10}$ M$_{\odot}$ is spread across a total extent of ~60 kpc. Most of this widespread CO(1-0) appears to follow prominent tidal features visible in the rest-frame near-UV HST/WFPC2 imaging. Ly$\alpha$ emission shows an excess over HeII, but a deficiency over L(IR), which is likely the result of photo-ionisation by enhanced but very obscured star formation that was triggered by the merger. In terms of feedback, the radio source is aligned with widespread CO(1-0) emission, which suggests that there is a physical link between the propagating radio jets and the presence of cold molecular gas on scales of the galaxy's halo. Its optical appearance, combined with the transformational stage at which we witness the evolution of MRC 0152-209, leads us to adopt the name `Dragonfly Galaxy'., Comment: Accepted for publication in MNRAS (12 pages, 5 figures)

The Square Kilometre Array will be a revolutionary instrument for the study of gas in the distant Universe. SKA1 will have sufficient sensitivity to detect and image atomic 21 cm HI in individual galaxies at significant cosmological distances, complementing ongoing ALMA imaging of redshifted high-J CO line emission and far-infrared interstellar medium lines such as [CII] 157.7 um. At frequencies below ~50 GHz, observations of redshifted emission from low-J transitions of CO, HCN, HCO+, HNC, H2O and CS provide insight into the kinematics and mass budget of the cold, dense star-forming gas in galaxies. In advance of ALMA band 1 deployment (35 to 52 GHz), the most sensitive facility for high-redshift studies of molecular gas operating below 50~GHz is the Karl G. Jansky Very Large Array (VLA). Here, we present an overview of the role that the SKA could play in molecular emission line studies during SKA1 and SKA2, with an emphasis on studies of the dense gas tracers directly probing regions of active star-formation., Comment: 10 pages, to be published in Proceedings of Science: "Advancing Astrophysics with the Square Kilometre Array", PoS(AASKA14)161

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

We present high-resolution observations of the 880 $\mu$m (rest-frame FIR) continuum emission in the z$=$4.05 submillimeter galaxy GN20 from the IRAM Plateau de Bure Interferometer (PdBI). These data resolve the obscured star formation in this unlensed galaxy on scales of 0.3$^{\prime\prime}$$\times$0.2$^{\prime\prime}$ ($\sim$2.1$\times$1.3 kpc). The observations reveal a bright (16$\pm$1 mJy) dusty starburst centered on the cold molecular gas reservoir and showing a bar-like extension along the major axis. The striking anti-correlation with the HST/WFC3 imaging suggests that the copious dust surrounding the starburst heavily obscures the rest-frame UV/optical emission. A comparison with 1.2 mm PdBI continuum data reveals no evidence for variations in the dust properties across the source within the uncertainties, consistent with extended star formation, and the peak star formation rate surface density (119$\pm$8 M$_{\odot}$ yr$^{-1}$ kpc$^{-2}$) implies that the star formation in GN20 remains sub-Eddington on scales down to 3 kpc$^2$. We find that the star formation efficiency is highest in the central regions of GN20, leading to a resolved star formation law with a power law slope of $\Sigma_{\rm SFR}$ $\sim$ $\Sigma_{\rm H_2}^{\rm 2.1\pm1.0}$, and that GN20 lies above the sequence of normal star-forming disks, implying that the dispersion in the star formation law is not due solely to morphology or choice of conversion factor. These data extend previous evidence for a fixed star formation efficiency per free-fall time to include the star-forming medium on $\sim$kpc-scales in a galaxy 12 Gyr ago., Comment: 6 pages, 5 figures, accepted to ApJL

The broad spectral bandwidth at mm and cm-wavelengths provided by the recent upgrades to the Karl G. Jansky Very Large Array (VLA) has made it possible to conduct unbiased searches for molecular CO line emission at redshifts, z > 1.31. We present the discovery of a gas-rich, star-forming galaxy at z = 2.48, through the detection of CO(1-0) line emission in the COLDz survey, through a sensitive, Ka-band (31 to 39 GHz) VLA survey of a 6.5 square arcminute region of the COSMOS field. We argue that the broad line (FWHM ~570 +/- 80 km/s) is most likely to be CO(1-0) at z=2.48, as the integrated emission is spatially coincident with an infrared-detected galaxy with a photometric redshift estimate of z = 3.2 +/- 0.4. The CO(1-0) line luminosity is L'_CO = (2.2 +/- 0.3) x 10^{10} K km/s pc^2, suggesting a cold molecular gas mass of M_gas ~ (2 - 8)x10^{10}M_solar depending on the assumed value of the molecular gas mass to CO luminosity ratio alpha_CO. The estimated infrared luminosity from the (rest-frame) far-infrared spectral energy distribution (SED) is L_IR = 2.5x10^{12} L_solar and the star-formation rate is ~250 M_solar/yr, with the SED shape indicating substantial dust obscuration of the stellar light. The infrared to CO line luminosity ratio is ~114+/-19 L_solar/(K km/s pc^2), similar to galaxies with similar SFRs selected at UV/optical to radio wavelengths. This discovery confirms the potential for molecular emission line surveys as a route to study populations of gas-rich galaxies in the future.

We present detections of the CO(J= 1-0) emission line in a sample of four massive star-forming galaxies at z~1.5-2.2 obtained with the Karl G. Jansky Very Large Array (VLA). Combining these observations with previous CO(2-1) and CO(3-2) detections of these galaxies, we study the excitation properties of the molecular gas in our sample sources. We find an average line brightness temperature ratios of R_{21}=0.70+\-0.16 and R_{31}=0.50+\-0.29, based on measurements for three and two galaxies, respectively. These results provide additional support to previous indications of sub-thermal gas excitation for the CO(3-2) line with a typically assumed line ratio R_{31}~0.5. For one of our targets, BzK-21000, we present spatially resolved CO line maps. At the resolution of 0.18 arcsec (1.5 kpc), most of the emission is resolved out except for some clumpy structure. From this, we attempt to identify molecular gas clumps in the data cube, finding 4 possible candidates. We estimate that <40 % of the molecular gas is confined to giant clumps (~1.5 kpc in size), and thus most of the gas could be distributed in small fainter clouds or in fairly diffuse extended regions of lower brightness temperatures than our sensitivity limit., Comment: Accepted 26 April 2014; Title modified to match published version

We report optical VLT FORS2 spectroscopy of the two Ly-alpha emitters (LAEs) companions to the quasi-stellar object (QSO) - sub-millimetre galaxy (SMG) system BRI1202-0725 at z = 4.7, which have recently been detected in the [CII]158um line by the Atacama Large Millimetre/Sub-millimetre Array (ALMA). We detect Ly-alpha emission from both sources and so confirm that these Ly-alpha emitter candidates are physically associated with the BRI1202- 0725 system. We also report the lack of detection of any high ionisation emission lines (N V, Si IV, C IV and He II) and find that these systems are likely not photoionised by the quasar, leaving in situ star formation as the main powering source of these LAEs. We also find that both LAEs have Ly-alpha emission much broader (1300 km/s) than the [CII] emission and broader than most LAEs. In addition, both LAEs have roughly symmetric Ly-alpha profiles implying that both systems are within the HII sphere produced by the quasar. This is the first time that the proximity zone of a quasar is probed by exploiting nearby Ly-alpha emitters. We discuss the observational properties of these galaxies in the context of recent galaxy formation models., Comment: 6 pages, 5 figures, accepted for publication in MNRAS

We present Karl G. Jansky Very Large Array (VLA) observations of 44 GHz continuum and CO J=2-1 line emission in BR1202-0725 at z=4.7 (a starburst galaxy and quasar pair) and BRI1335-0417 at z=4.4 (also hosting a quasar). With the full 8 GHz bandwidth capabilities of the upgraded VLA, we study the (rest-frame) 250 GHz thermal dust continuum emission for the first time along with the cold molecular gas traced by the Low-J CO line emission. The measured CO J=2-1 line luminosities of BR1202-0725 are L'(CO) = (8.7+/-0.8)x10^10 K km/s pc^2 and L'(CO) = (6.0+/-0.5)x10^10 K km/s pc^2 for the submm galaxy (SMG) and quasar, which are equal to previous measurements of the CO J=5-4 line luminosities implying thermalized line emission and we estimate a combined cold molecular gas mass of ~9x10^10 Msun. In BRI1335-0417 we measure L'(CO) = (7.3+/-0.6)x10^10 K km/s pc^2. We detect continuum emission in the SMG BR1202-0725 North (S(44GHz) = 51+/-6 microJy), while the quasar is detected with S(44GHz) = 24+/-6 microJy and in BRI1335-0417 we measure S(44GHz) = 40+/-7 microJy. Combining our continuum observations with previous data at (rest-frame) far-infrared and cm-wavelengths, we fit three component models in order to estimate the star-formation rates. This spectral energy distribution fitting suggests that the dominant contribution to the observed 44~GHz continuum is thermal dust emission, while either thermal free-free or synchrotron emission contributes less than 30%., Comment: 16 pages , 4 figures. Accepted for publication in ApJ

We present a CO(1-0) survey for cold molecular gas in a representative sample of 13 high-z radio galaxies (HzRGs) at 1.4

We present the methodology for ``blind'' millimetre-wave surveys for redshifted molecular absorption in the CO/HCO$^+$ rotational lines. The frequency range $30-50$ GHz appears optimal for such surveys, providing sensitivity to absorbers at $z \gtrsim 0.85$. It is critical that the survey is ``blind'', i.e. based on a radio-selected sample, including sources without known redshifts. We also report results from the first large survey of this kind, using the Q-band receiver on the Green Bank Telescope (GBT) to search for molecular absorption towards 36 sources, 3 without known redshifts, over the frequency range $39.6 - 49.5$ GHz. The GBT survey has a total redshift path of $\Delta z \approx 24$, mostly at $0.81 < z < 1.91$, and a sensitivity sufficient to detect equivalent ${\rm H_2}$ column densities $\gtrsim 3 \times 10^{21}$ cm$^{-2}$ in absorption at $5\sigma$ significance (using the CO-to-${\rm H_2}$ and HCO$^+$-to-${\rm H_2}$ conversion factors of the Milky Way). The survey yielded no confirmed detections of molecular absorption, yielding the $2\sigma$ upper limit $n(z=1.2) < 0.15$ on the redshift number density of molecular gas at column densities $N({\rm H_2}) \gtrsim 3 \times 10^{21}$ cm$^{-2}$., Comment: 6 pages, 2 figures. Accepted for publication in The Astrophysical Journal

We present 1-2 GHz Very Large Array A-configuration continuum observations on the highest redshift quasar known to date, the $z=7.085$ quasar ULAS J112001.48+064124.3. The results show no radio continuum emission at the optical position of the quasar or its vicinity at a level of $\geq 3\sigma$ or $23.1 \mu$Jy beam$^{-1}$. This $3\sigma$ limit corresponds to a rest frame 1.4 GHz luminosity density limit of $L_{\nu,1.4\,GHz} < 1.76 \times 10^{24}$ W Hz$^{-1}$ for a spectral index of $\alpha=0$, and $L_{\nu,1.4\,GHz} < 1.42 \times 10^{25}$ W Hz$^{-1}$ for a spectral index of $\alpha=-1$. The rest-frame 1.4 GHz luminosity limits are $L_{rad} < 6.43 \times 10^6 L_{\odot}$ and $L_{\rm rad} < 5.20 \times 10^7 L_{\odot}$ for $\alpha=0$ and $\alpha=-1$, respectively. The derived limits for the ratio of the rest frame 1.4 GHz luminosity density to the $B$-band optical luminosity density are $R\rlap{}_{1.4}^{*} < 0.53$ and $< 4.30$ for the above noted spectral indices, respectively. Given our upper limits on the radio continuum emission and the radio-to-optical luminosity ratio, we conclude that this quasar is radio-quiet and located at the low end of the radio quiet distribution of high redshift ($z \gtrsim 6$) quasars., Comment: 11 pages, 3 figures, accepted for publication in AJ

We present a study of the formation of clustered, massive galaxies at large look-back times via spectroscopic imaging of CO in the unique GN20 proto-cluster at z = 4.05. Existing observations show that this is a dense concentration of gas-rich, very active star forming galaxies, including multiple bright submillimeter galaxies (SMGs). Using deep, high-resolution VLA CO(2-1) observations, we image the molecular gas with a resolution of ~1 kpc just 1.6 Gyr after the Big Bang. The SMGs GN20.2a and GN20.2b have deconvolved sizes of ~5 kpc X 3 kpc and ~8 kpc X 5 kpc (Gaussian FWHM) in CO(2-1), respectively, and we measure gas surface densities up to ~12,700/1,700X(sin i) (\alpha_CO/0.8) M_sun/pc^2 for GN20.2a/GN20.2b in the highest-resolution maps. Dynamical mass estimates allow us to constrain the CO-to-H_2 conversion factor to \alpha_CO = 1.7+/-0.8 M_sun (K km s^{-1} pc^2)^-1 for GN20.2a and \alpha_CO = 1.1+/-^{1.5}_{1.1} M_sun (K km/s pc^2)^-1 for GN20.2b. We measure significant offsets (0.5"-1") between the CO and optical emission, indicating either dust obscuration on scales of tens of kpc or that the emission originates from distinct galaxies. CO spectral line energy distributions imply physical conditions comparable to other SMGs and reveal further evidence that GN20.2a and GN20.2b are in different merging stages. We carry out a targeted search for CO emission from the 14 known B-band Lyman break galaxies (LBGs) in the field, tentatively detecting CO in a previously-undetected LBG and placing 3sigma upper limits on the CO luminosities of those that may lie within our bandpass. A blind search for emission line sources down to a 5sigma limiting CO luminosity of L'_CO(2-1) = 8 X 10^9 K km/s pc^2 and covering \Delta z = 0.0273 (~20 comoving Mpc) produces no other strong contenders associated with the proto-cluster., Comment: Accepted to ApJ

The high-redshift radio galaxy MRC 1138-262 (`Spiderweb Galaxy'; z = 2.16), is one of the most massive systems in the early Universe and surrounded by a dense `web' of proto-cluster galaxies. Using the Australia Telescope Compact Array, we detected CO(1-0) emission from cold molecular gas -- the raw ingredient for star formation -- across the Spiderweb Galaxy. We infer a molecular gas mass of M(H2) = 6x10^10 M(sun) (for M(H2)/L'(CO)=0.8). While the bulk of the molecular gas coincides with the central radio galaxy, there are indications that a substantial fraction of this gas is associated with satellite galaxies or spread across the inter-galactic medium on scales of tens of kpc. In addition, we tentatively detect CO(1-0) in the star-forming proto-cluster galaxy HAE 229, 250 kpc to the west. Our observations are consistent with the fact that the Spiderweb Galaxy is building up its stellar mass through a massive burst of widespread star formation. At maximum star formation efficiency, the molecular gas will be able to sustain the current star formation rate (SFR ~ 1400 M(sun)/yr, as traced by Seymour et al.) for about 40 Myr. This is similar to the estimated typical lifetime of a major starburst event in infra-red luminous merger systems., Comment: 7 pages, 3 figures - accepted for publication in MNRAS

We present further analysis of the [CII] 158$\mu$m fine structure line and thermal dust continuum emission from the archetype extreme starburst/AGN group of galaxies in the early Universe, BRI 1202-0725 at $z=4.7$, using the Atacama Large Millimeter Array. The group is long noted for having a closely separated (26kpc in projection) FIR-hyperluminous quasar host galaxy and an optically obscured submm galaxy (SMG). A short ALMA test observation reveals a rich laboratory for the study of the myriad processes involved in clustered massive galaxy formation in the early Universe. Strong [CII] emission from the SMG and the quasar have been reported earlier by Wagg et al. (2012) based on these observations. In this letter, we examine in more detail the imaging results from the ALMA observations, including velocity channel images, position-velocity plots, and line moment images. We present detections of [CII] emission from two Ly$\alpha$-selected galaxies in the group, demonstrating the relative ease with which ALMA can detect the [CII] emission from lower star formation rate galaxies at high redshift. Imaging of the [CII] emission shows a clear velocity gradient across the SMG, possibly indicating rotation or a more complex dynamical system on a scale $\sim 10$kpc. There is evidence in the quasar spectrum and images for a possible outflow toward the southwest, as well as more extended emission (a 'bridge'), between the quasar and the SMG, although the latter could simply be emission from Ly$\alpha$-1 blending with that of the quasar at the limited spatial resolution of the current observations. These results provide an unprecedented view of a major merger of gas rich galaxies driving extreme starbursts and AGN accretion during the formation of massive galaxies and supermassive black holes within 1.3 Gyr of the Big Bang., Comment: for the astrophysical journal. 5 figs

We present Karl G. Jansky Very Large Array (VLA) observations of the CO(2-1) emission in the z=4.05 submillimeter galaxy (SMG) GN20. These high-resolution data allow us to image the molecular gas at 1.3 kpc resolution just 1.6 Gyr after the Big Bang. The data reveal a clumpy, extended gas reservoir, 14 +/- 4 kpc in diameter, in unprecedented detail. A dynamical analysis shows that the data are consistent with a rotating disk of total dynamical mass 5.4 +/- 2.4 X 10^11 M_sun. We use this dynamical mass estimate to constrain the CO-to-H_2 mass conversion factor (alpha_CO), finding alpha_CO=1.1 +/- 0.6 M_sun (K km s^-1 pc^2)^-1. We identify five distinct molecular gas clumps in the disk of GN20 with masses a few percent of the total gas mass, brightness temperatures of 16-31K, and surface densities of >3,200-4,500 X (alpha_CO/0.8) M_sun pc^-2. Virial mass estimates indicate they could be self-gravitating, and we constrain their CO-to-H_2 mass conversion factor to be <0.2-0.7 M_sun (K km s^-1 pc^2)^-1. A multiwavelength comparison demonstrates that the molecular gas is concentrated in a region of the galaxy that is heavily obscured in the rest-frame UV/optical. We investigate the spatially-resolved gas excitation and find that the CO(6-5)/CO(2-1) ratio is constant with radius, consistent with star formation occuring over a large portion of the disk. We discuss the implications of our results in the context of different fueling scenarios for SMGs., Comment: 15 pages, 9 figures, accepted for publication in ApJ

We report results from a deep Jansky Very Large Array (JVLA) search for CO 1-0 line emission from galaxies in a candidate galaxy cluster at z~1.55 in the COSMOS field. We target 4 galaxies with optical spectroscopic redshifts in the range z=1.47-1.59. Two of these 4 galaxies, ID51613 and ID51813, are nominally detected in CO line emission at the 3-4 sigma level. We find CO luminosities of 2.4x10^10 K km/s pc^2 and 1.3x10^10 K km/s pc^2, respectively. Taking advantage from the clustering and 2-GHz bandwidth of the JVLA, we perform a search for emission lines in the proximity of optical sources within the field of view of our observations. We limit our search to galaxies with K<23.5 (AB) and z_phot=1.2-1.8. We find 2 bright optical galaxies to be associated with significant emission line peaks (>4 sigma) in the data cube, which we identify with the CO line emission. To test the reliability of the line peaks found, we performed a parallel search for line peaks using a Bayesian inference method. Monte Carlo simulations show that such associations are statistically significant, with probabilities of chance association of 3.5% and 10.7% for ID 51207 and ID 51380, respectively. Modeling of their optical/IR SEDs indicates that the CO detected galaxies and candidates have stellar masses and SFRs in the range (0.3-1.1)x10^11 M_sun and 60-160 M_sun/yr, with SFEs comparable to that found in other star-forming galaxies at similar redshifts. By comparing the space density of CO emitters derived from our observations with the space density derived from previous CO detections at z~1.5, and with semi-analytic predictions for the CO luminosity function, we suggest that the latter tend to underestimate the number of CO galaxies detected at high-redshift. Finally, we argue about the benefits of future blind CO searches in clustered fields with upcoming submm/radio facilities., Comment: Accepted for publication in MNRAS. Abstract has been slightly shortened compared to original pdf version

We present new IRAM PdBI 1.3mm continuum observations at ~1.5" resolution of 28 SMGs previously discovered with the 870um bolometer LABOCA at APEX within the central 0.7deg2 of the COSMOS field. 19 out of the 28 LABOCA sources were detected with the PdBI at a >~3sigma level of ~1.4mJy/b. A combined analysis of this new sample with existing interferometrically identified SMGs in the COSMOS field yields the following results: 1) >~15%, and possibly up to ~40% of single-dish detected SMGs consist of multiple sources, 2) statistical identifications of multi-wavelength counterparts to the single-dish SMGs yield that only ~50% of these single-dish SMGs have real radio or IR counterparts, 3) ~18% of interferometric SMGs have only radio or even no multi-wavelength counterpart at all, and 4) ~50-70% of z>~3 SMGs have no radio counterparts down to an rms of 7-12uJy at 1.4GHz. Using the exact interferometric positions to identify proper multi-wavelength counterparts allows us to determine accurate photometric redshifts for these sources. The redshift distributions of the combined and the individual 1.1mm and 870um selected samples have a higher mean and broader width than the redshift distributions derived in previous studies. Our sample supports the previous tentative trend that on average brighter and/or mm-selected SMGs are located at higher redshifts. There is a tentative offset between the mean redshift for the 1.1mm (=3.1+/-0.4) and 870um (=2.6+/-0.4) selected samples, with the 1.1mm sources lying on average at higher redshifts. Based on our nearly complete sample of AzTEC 1.1mm SMGs within a uniform 0.15deg2 area we infer a higher surface density of z>~4 SMGs than predicted by current cosmological models. In summary, our findings imply that (sub-)millimeter interferometric identifications are crucial to build statistically complete and unbiased samples of SMGs., Comment: 35 pages, 18 figures, 10 tables; accepted for publication in A&A