Source author record

F. Bertoldi

F. Bertoldi 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

72works
6topics
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

72 published item(s)

preprint2021arXiv

[CII] line intensity mapping the epoch of reionization with the Prime-Cam on FYST

We predict the three-dimensional intensity power spectrum (PS) of the [CII] 158$\,μ$m line throughout the epoch of (and post) reionization at redshifts from $\approx$ 3.5 to 8. We study the detectability of the PS in a line intensity mapping (LIM) survey with the Fred Young Submillimeter Telescope (FYST). We created mock [CII] tomographic scans in redshift bins at $z\approx$ 3.7, 4.3, 5.8, and 7.4 using the Illustris TNG300-1 $Λ$CDM simulation and adopting a relation between the star formation activity and the [CII] luminosity ($L_{[CII]}$) of galaxies. A star formation rate (SFR) was assigned to a dark matter halo in the Illustris simulation in two ways: (i) we adopted the SFR computed in the Illustris simulation and, (ii) we matched the abundance of the halos with the SFR traced by the observed dust-corrected ultraviolet luminosity function of high-redshift galaxies. The $L_{[CII]}$ is related to the SFR from a semi-analytic model of galaxy formation, from a hydrodynamical simulation of a high-redshift galaxy, or from a high-redshift [CII] galaxy survey. The [CII] intensity PS was computed from mock tomographic scans to assess its detectability with the anticipated observational capability of the FYST. The amplitude of the predicted [CII] intensity power spectrum varies by more than a factor of 10, depending on the choice of the halo-to-galaxy SFR and the SFR-to-$L_{[CII]}$ relations. In the planned $4^{\circ} \times 4^{\circ}$ FYST LIM survey, we expect a detection of the [CII] PS up to $z \approx$ 5.8, and potentially even up to $z \approx $ 7.4. The design of the envisioned FYST LIM survey enables a PS measurement not only in small (<10 Mpc) shot noise-dominated scales, but also in large (>50 Mpc) clustering-dominated scales making it the first LIM experiment that will place constraints on the SFR-to-$L_{[CII]}$ and the halo-to-galaxy SFR relations simultaneously.

preprint2021arXiv

Large format imaging spectrograph for the Large Submillimeter Telescope (LST)

We present a conceptual study of a large format imaging spectrograph for the Large Submillimeter Telescope (LST) and the Atacama Large Aperture Submillimeter Telescope (AtLAST). Recent observations of high-redshift galaxies indicate the onset of earliest star formation just a few 100 million years after the Big Bang (i.e., z = 12--15), and LST/AtLAST will provide a unique pathway to uncover spectroscopically-identified first forming galaxies in the pre-reionization era, once it will be equipped with a large format imaging spectrograph. We propose a 3-band (200, 255, and 350 GHz), medium resolution (R = 2,000) imaging spectrograph with 1.5 M detectors in total based on the KATANA concept (Karatsu et al. 2019), which exploits technologies of the integrated superconducting spectrometer (ISS) and a large-format imaging array. A 1-deg2 drilling survey (3,500 hr) will capture a large number of [O III] 88 um (and [C II] 158 um) emitters at z = 8--9, and constrain [O III] luminosity functions at z > 12.

preprint2020arXiv

The redshift and star formation mode of AzTEC2: a pair of massive galaxies at $z=4.63$

We combine observations from the Atacama Large Millimeter/submillimeter Array (ALMA) and the NOrthern Extended Millimeter Array (NOEMA) to assess the redshift and to study the star formation conditions in AzTEC2: one of the brightest sub-millimeter galaxies (SMGs) in the COSMOS field ($S_{\rm 1.1mm}=10.5\pm1.4$mJy). Our high-resolution observations confirm that AzTEC2 splits into two components (namely AzTEC2-A and AzTEC2-B) for which we detect [C$\,$II] and $^{12}$CO(5$\to$4) line emission, implying a redshift of $4.626\pm0.001$ ($4.633\pm0.001$) for AzTEC2-A (AzTEC2-B) and ruling out previous associations with a galaxy at $z\sim1$. We use the $^{12}$CO(5$\to$4) line emission and adopt typical SMG-like gas excitation conditions to estimate the molecular gas mass, which is $M_{\rm gas}(α_{\rm CO}/2.5)=2.1\pm0.4 \times10^{11}{\rm M}_\odot$ for AzTEC2-A, and a factor four lower for AzTEC2-B. With the infrared-derived star formation rate of AzTEC2-A ($1920\pm100 \,M_\odot{\rm \, yr}^{-1}$) and AzTEC2-B ($710\pm 35\,M_\odot{\rm \,yr}^{-1}$), they both will consume their current gas reservoir within $(30-200)\,$Myr. We find evidence of a rotation-dominated [C$\,$II] disk in AzTEC2-A, with a de-projected rotational velocity of $v_{\rm rot}(i=39^\circ)=660\pm130{\rm \,km\,s}^{-1}$, velocity dispersion $\lesssim100{\rm \,km\,s}^{-1}$, and dynamical mass of $M_{\rm dyn}(i=39^\circ)=2.6^{+1.2}_{-0.9}\times10^{11}\,M_\odot$. We propose that an elevated gas accretion rate from the cosmic web might be the main driver of the intense levels of star formation in AzTEC2-A, which might be further enhanced by gravitational torques induced by its minor companion (AzTEC2-B). These results strengthen the picture whereby the population of single-dish selected SMGs is rather heterogeneous, including a population of pairs of massive, highly-active galaxies in a pre-coalescence phase.

preprint2020arXiv

Turbulent Gas in Lensed Planck-selected Starbursts at redshifts 1-3.5

Dusty star-forming galaxies at high redshift (1 < z < 3) represent the most intense star-forming regions in the Universe. Key aspects to these processes are the gas heating and cooling mechanisms. Although it is well known that these galaxies are gas-rich, little is known about the gas excitation conditions. Here we examine these processes in a sample of 24 strongly lensed star-forming galaxies identified by the \textit{Planck} satellite (LPs) at z ~ 1.1 - 3.5. We analyze 162 CO rotational transitions (ranging from Jupper = 1 - 12) and 37 atomic carbon fine-structure lines ([CI]) in order to characterize the physical conditions of the gas in sample of LPs. We simultaneously fit the CO and [CI] lines, and the dust continuum emission, using two different non-LTE, radiative transfer models. The first model represents a two component gas density, while the second assumes a turbulence driven log-normal gas density distribution. These LPs are among the most gas-rich, infrared (IR) luminous galaxies ever observed ($μ_{\rm L}$L$_{\rm IR(8-1000μm) } \sim 10^{13-14.6} $\Lsun; $< μ_{\rm L}$M$_{\rm ISM}> = 2.7 \pm 1.2 \times 10^{12}$ \Msun, with $μ_{\rm L} \sim 10-30$ the average lens magnification factor). Our results suggest that the turbulent ISM present in the LPs can be well-characterized by a high turbulent velocity dispersion ($<ΔV_{\rm turb}> \sim 100 $ \kms) and gas kinetic temperature to dust temperature ratios $<T_{\rm kin}$/$T_{\rm d}> \sim 2.5$, sustained on scales larger than a few kpc. We speculate that the average surface density of the molecular gas mass and IR luminosity $Σ_{\rm M_{\rm ISM}}$ $\sim 10^{3 - 4}$ \Msun pc$^{-2}$ and $Σ_{\rm L_{\rm IR}}$ $\sim 10^{11 - 12}$ \Lsun kpc$^{-2}$, arise from both stellar mechanical feedback and a steady momentum injection from the accretion of intergalactic gas.

preprint2019arXiv

Automated Mining of the ALMA Archive in the COSMOS Field (A3COSMOS): II. Cold Molecular Gas Evolution out to Redshift 6

We present new measurements of the cosmic cold molecular gas evolution out to redshift 6 based on systematic mining of the ALMA public archive in the COSMOS deep field (A3COSMOS). Our A3COSMOS dataset contains ~700 galaxies (0.3 < z < 6) with high-confidence ALMA detections in the (sub-)millimeter continuum and multi-wavelength spectral energy distributions (SEDs). Multiple gas mass calibration methods are compared and biases in band conversions (from observed ALMA wavelength to rest-frame Rayleigh-Jeans(RJ)-tail continuum) have been tested. Combining our A3COSMOS sample with ~1,000 CO-observed galaxies at 0 < z < 4 (75% at z < 0.1), we parameterize galaxies' molecular gas depletion time and molecular gas to stellar mass ratio (gas fraction) each as a function of the stellar mass, offset from the star-forming main sequence (Delta MS) and cosmic age (or redshift). Our proposed functional form provides a statistically better fit to current data (than functional forms in the literature), and implies a "downsizing" effect (i.e., more-massive galaxies evolve earlier than less-massive ones) and "mass-quenching" (gas consumption slows down with cosmic time for massive galaxies but speeds up for low-mass ones). Adopting galaxy stellar mass functions and applying our function for gas mass calculation, we for the first time infer the cosmic cold molecular gas density evolution out to redshift 6 and find agreement with CO blind surveys as well as semi-analytic modeling. These together provide a coherent picture of cold molecular gas, SFR and stellar mass evolution in galaxies across cosmic time.

preprint2016arXiv

ALMA spectroscopic survey in the Hubble Ultra Deep Field: Continuum number counts, resolved 1.2-mm extragalactic background, and properties of the faintest dusty star forming galaxies

We present an analysis of a deep (1$σ$=13 $μ$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σ$ 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.

preprint2016arXiv

ALMA Spectroscopic Survey in the Hubble Ultra Deep Field: implications for spectral line intensity mapping at millimeter wavelengths and CMB spectral distortions

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$μ$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$ $μ$K. In the 242GHz band, the mean brightness is: $T_B = 0.55\pm 0.033$ $μ$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.

preprint2016arXiv

ALMA Spectroscopic Survey in the Hubble Ultra Deep Field: Search for [CII] line and dust emission in 6<z<8 galaxies

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<z<8$ for [CII] line emission and reach a limiting luminosity of L$_{\rm [CII]}\sim$(1.6-2.5)$\times$10$^{8}$ L$_{\odot}$. We identify fourteen [CII] line emitting candidates in this redshift range with significances $>$4.5 $σ$, 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μ$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<z<8$ that is significantly higher than predicted by current models and some extrapolations from galaxies in the local universe.

preprint2016arXiv

Kiloparsec-scale dust disks in high-redshift luminous submillimeter galaxies

We present high-resolution (0.16$"$) 870um Atacama Large Millimeter/submillimeter Array (ALMA) imaging of 16 luminous (L_IR ~ 4 x 10^12 L_sun) submillimeter galaxies (SMGs) from the ALESS survey of the Extended Chandra Deep Field South. This dust imaging traces the dust-obscured star formation in these z~2.5 galaxies on ~1.3 kpc scales. The emission has a median effective radius of $R_e=0.24" \pm 0.02"$, corresponding to a typical physical size of $R_{e}=1.8\pm$0.2 kpc. We derive a median Sérsic index of $n=0.9\pm0.2$, implying that the dust emission is remarkably disk-like at the current resolution and sensitivity. We use different weighting schemes with the visibilities to search for clumps on 0.12$"$ (~1.0 kpc) scales, but we find no significant evidence for clumping in the majority of cases. Indeed, we demonstrate using simulations that the observed morphologies are generally consistent with smooth exponential disks, suggesting that caution should be exercised when identifying candidate clumps in even moderate S/N interferometric data. We compare our maps to comparable-resolution HST H$_{160}$-band images, finding that the stellar morphologies appear significantly more extended and disturbed, and suggesting that major mergers may be responsible for driving the formation of the compact dust disks we observe. The stark contrast between the obscured and unobscured morphologies may also have implications for SED fitting routines that assume the dust is co-located with the optical/near-IR continuum emission. Finally, we discuss the potential of the current bursts of star formation to transform the observed galaxy sizes and light profiles, showing that the $z\sim0$ descendants of these SMGs are expected to have stellar masses, effective radii, and gas surface densities consistent with the most compact massive (M_* ~ 1-2 x 10^11 M_sun) early-type galaxies observed locally.

preprint2016arXiv

Microphysics and dynamics of the Gamma-Ray Burst 121024A

Aims. The aim of the study is to constrain the physics of gamma-ray bursts (GRBs) by analysing the multi-wavelength afterglow data set of GRB 121024A that covers the full range from radio to X-rays. Methods. Using multi-epoch broad-band observations of the GRB 121024A afterglow, we measured the three characteristic break frequencies of the synchrotron spectrum. We used six epochs of combined XRT and GROND data to constrain the temporal slopes, the dust extinction, the X-ray absorption, and the spectral slope with high accuracy. Two more epochs of combined data from XRT, GROND, APEX, CARMA, and EVLA were used to set constraints on the break frequencies and therefore on the micro-physical and dynamical parameters. Results. The XRT and GROND light curves show a simultaneous and achromatic break at around 49 ks. As a result, the crossing of the synchrotron cooling break is no suitable explanation for the break in the light curve. The multi wavelength data allow us to test two plausible scenarios explaining the break: a jet break, and the end of energy injection. The jet-break scenario requires a hard electron spectrum, a very low cooling break frequency, and a non-spreading jet. The energy injection avoids these problems, but requires $ε_e > 1 (k = 2)$, spherical outflow, and $ε_B < 10^{-9}$. Conclusions. In light of the extreme microphysical parameters required by the energy-injection model, we favour a jet-break scenario where $ν_m < ν_{sa}$ to explain the observations. This scenario gives physically meaningful microphysical parameters, and it also naturally explains the reported detection of linear and circular polarisation.

preprint2015arXiv

A blind CO detection of a Distant Red Galaxy in the HS1700+64 proto-cluster

We report the blind detection of 12CO emission from a Distant Red Galaxy, HS1700.DRG55. We have used the IRAM PdBI-WIDEX, with its 3.6GHz of instantaneous dual-polarization bandwidth, to target 12CO(3--2) from galaxies lying in the proto-cluster at z=2.300 in the field HS1700+64. If indeed this line in DRG55 is 12CO(3--2), it's detection at 104.9GHz indicates a z_CO=2.296. None of the other eight known z~2.30 proto-cluster galaxies lying within the primary beam (PB) are detected in 12CO, although the limits are ~2x worse towards the edge of the PB where several lie. The optical/near-IR magnitudes of DRG55 (R_AB>27, K_AB=22.3) mean that optical spectroscopic redshifts are difficult with 10m-class telescopes, but near-IR redshifts would be feasible. The 24um-implied SFR (210 M_odot yr-1), stellar mass (~10^11 M-odot) and 12CO line luminosity (3.6x10^10 K km s-1 pc^2) are comparable to other normal 12CO-detected star forming galaxies in the literature, although the galaxy is some ~2 mag (~6x) fainter in the rest-frame UV than 12CO-detected galaxies at z>2. The detection of DRG55 in 12CO complements three other 12CO detected UV-bright galaxies in this proto-cluster from previous studies, and suggests that many optically faint galaxies in the proto-cluster may host substantial molecular gas reservoirs, and a full blind census of 12CO in this overdense environment is warranted.

preprint2015arXiv

A millimeter-wave redshift search for the unlensed HyLIRG, HS1700.850.1

We report the redshift of an unlensed, highly obscured submillimetre galaxy (SMG), HS1700.850.1, the brightest SMG (S850um =19.1 mJy) detected in the JCMT/SCUBA-2 Baryonic Structure Survey, based on the detection of its CO line emission. Using the IRAM PdBI-WIDEX with 3.6GHz band width, we serendipitously detect an emission line at 150.6 GHz. From a search over 14.5 GHz in the 3-mm and 2-mm atmospheric windows, we confirm the identification of this line as CO(5-4) at z = 2.816, meaning that it does not reside in the z~2.30 proto-cluster in this field. Measurement of the 870um source size (<0.85") from the Sub-Millimeter Array (SMA) confirms a compact emission in a S870um =14.5mJy, LIR~10^13 Lsun component, suggesting an Eddington-limited starburst. We use the double-peaked CO line profile measurements along with the SMA size constraints to study the gas dynamics of a HyLIRG, estimating the gas and dynamical masses of HS1700.850.1. While HS1700.850.1 is one of the most extreme galaxies known in the Universe, we find that it occupies a relative void in the Lyman-Break Galaxy distribution in this field. Comparison with other extreme objects at similar epochs (HyLIRG Quasars), and cosmological simulations, suggests such an anti-bias of bright SMGs could be relatively common, with the brightest SMGs rarely occupying the most overdense regions at z=2-4.

preprint2015arXiv

An Overview of the 2014 ALMA Long Baseline Campaign

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

preprint2015arXiv

Constraining the intracluster pressure profile from the thermal SZ power spectrum

The angular power spectrum of the thermal Sunyaev-Zel'dovich (tSZ) effect is highly sensitive to cosmological parameters such as sigma_8 and Omega_m, but its use as a precision cosmological probe is hindered by the astrophysical uncertainties in modeling the gas pressure profile in galaxy groups and clusters. In this paper we assume that the relevant cosmological parameters are accurately known and explore the ability of current and future tSZ power spectrum measurements to constrain the intracluster gas pressure or the evolution of the gas mass fraction, f_gas. We use the CMB bandpower measurements from the South Pole Telescope and a Bayesian Markov Chain Monte Carlo (MCMC) method to quantify deviations from the standard, universal gas pressure model. We explore analytical model extensions that bring the predictions for the tSZ power into agreement with experimental data. We find that a steeper pressure profile in the cluster outskirts or an evolving f_gas have mild-to-severe conflicts with experimental data or simulations. Varying more than one parameter in the pressure model leads to strong degeneracies that cannot be broken with current observational constraints. We use simulated bandpowers from future tSZ survey experiments, in particular a possible 2000 deg^2 CCAT survey, to show that future observations can provide almost an order of magnitude better precision on the same model parameters. This will allow us to break the current parameter degeneracies and place simultaneous constraints on the gas pressure profile and its redshift evolution, for example.

preprint2015arXiv

Physical properties of z>4 submillimeter galaxies in the COSMOS field

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<z<2 galaxies. We study the IR-radio correlation for our sources and find a deviation from that derived for z<3 ULIRGs (<q_IR>=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.

preprint2014arXiv

Dissecting the origin of the submillimeter emission in nearby galaxies with Herschel and LABOCA

We model the infrared to submillimeter spectral energy distribution of 11 nearby galaxies of the KINGFISH sample using Spitzer and Herschel data and compare model extrapolations at 870um (using different fitting techniques) with LABOCA 870um observations. We investigate how the differences between predictions and observations vary with model assumptions or environment. At global scales, we find that modified blackbody models using realistic cold emissivity indices (beta_c=2 or 1.5) are able to reproduce the 870um observed emission within the uncertainties for most of the sample. Low values (beta_c<1.3) would be required in NGC0337, NGC1512 and NGC7793. At local scales, we observe a systematic 870um excess when using beta_=2.0. The beta_c=1.5 or the Draine and Li (2007) models can reconcile predictions with observations in part of the disks. Some of the remaining excesses occur towards the centres and can be partly or fully accounted for by non-dust contributions such as CO(3-2) or, to a lesser extent, free-free or synchrotron emission. In three non-barred galaxies, the remaining excesses rather occur in the disk outskirts. This could be a sign of a flattening of the submm slope (and decrease of the effective emissivity index) with radius in these objects.

preprint2014arXiv

The far-infrared/radio correlation and radio spectral index of galaxies in the SFR-M* plane up to z 2

[Abridged] We study the evolution of the radio spectral index and far-infrared/radio correlation (FRC) across the star-formation rate-stellar masse (i.e. SFR-M*) plane up to z 2. We start from a M*-selected sample of galaxies with reliable SFR and redshift estimates. We then grid the SFR-M* plane in several redshift ranges and measure the infrared luminosity, radio luminosity, radio spectral index, and ultimately the FRC index (i.e. qFIR) of each SFR-M*-z bin. The infrared luminosities of our SFR-M*-z bins are estimated using their stacked far-infrared flux densities inferred from observations obtained with Herschel. Their radio luminosities and radio spectral indices (i.e. alpha, where Snu nu^-alpha) are estimated using their stacked 1.4GHz and 610MHz flux densities from the VLA and GMRT, respectively. Our far-infrared and radio observations include the most widely studied blank extragalactic fields -GOODS-N/S, ECDFS, and COSMOS- covering a sky area of 2deg^2. Using this methodology, we constrain the radio spectral index and FRC index of star-forming galaxies with M*>10^10Msun and 0<z<2.3. We find that alpha^1.4GHz_610MHz does not evolve significantly with redshift or with the distance of a galaxy with respect to the main sequence (MS) of the SFR-M* plane (i.e. Delta_log(SSFR)_MS=log[SSFR(galaxy)/SSFR_MS(M*,z)]). Instead, star-forming galaxies have a radio spectral index consistent with a canonical value of 0.8, which suggests that their radio spectra are dominated by non-thermal optically thin synchrotron emission. We find that qFIR displays a moderate but statistically significant redshift evolution as qFIR(z)=(2.35+/-0.08)*(1+z)^(-0.12+/-0.04), consistent with some previous literature. Finally, we find no significant correlation between qFIR and Delta_log(SSFR)_MS, though a weak positive trend, as observed in one of our redshift bins, cannot be firmly ruled out using our dataset.

preprint2014arXiv

Varying [CII]/[NII] line ratios in the interacting system BR1202-0725 at z=4.7

We study the properties of the interstellar medium in the interacting system BR1202-0725 at z=4.7 via its [NII] and [CII] fine-structure line emission. This system consists of a QSO, a sub-mm galaxy (SMG), and two Ly-alpha emitters (LAEs). Such a diversity in galaxy properties makes BR1202-0725 a unique laboratory of star formation and galaxy evolution at high redshift. We present ionized nitrogen ([NII] 205 micron) observations of this system, obtained with the IRAM Plateau de Bure Interferometer. We find no [NII] emission at the quasar location, but tentative [NII] line detections associated with the SMG and one of the LAEs. Together with available ionized carbon ([CII] 158 micron) ALMA observations of this system, we find the following: The [CII]/[NII] luminosity ratio is >5.5 for the QSO and the SMG, but it is as low as ~2 in the LAE, suggesting that, in this source, most of the [CII] emission is associated with the ionized medium (HII regions) rather than the neutral one (PDRs). This study demonstrates the importance of combined studies of multiple fine--structure lines in order to pin down the physical properties of the interstellar medium in distant galaxies.

preprint2013arXiv

A molecular line scan in the Hubble Deep Field North

We present a molecular line scan in the Hubble Deep Field North (HDF-N) that covers the entire 3mm window (79-115 GHz) using the IRAM Plateau de Bure Interferometer. Our CO redshift coverage spans z<0.45, 1<z<1.9 and all z>2. We reach a CO detection limit that is deep enough to detect essentially all z>1 CO lines reported in the literature so far. We have developed and applied different line searching algorithms, resulting in the discovery of 17 line candidates. We estimate that the rate of false positive line detections is ~2/17. We identify optical/NIR counterparts from the deep ancillary database of the HDF-N for seven of these candidates and investigate their available SEDs. Two secure CO detections in our scan are identified with star-forming galaxies at z=1.784 and at z=2.047. These galaxies have colors consistent with the `BzK' color selection and they show relatively bright CO emission compared with galaxies of similar dust continuum luminosity. We also detect two spectral lines in the submillimeter galaxy HDF850.1 at z=5.183. We consider an additional 9 line candidates as high quality. Our observations also provide a deep 3mm continuum map (1-sigma noise level = 8.6 $μ$Jy/beam). Via a stacking approach, we find that optical/MIR bright galaxies contribute only to <50% of the SFR density at 1<z<3, unless high dust temperatures are invoked. The present study represents a first, fundamental step towards an unbiased census of molecular gas in `normal' galaxies at high-z, a crucial goal of extragalactic astronomy in the ALMA era.

preprint2013arXiv

A Molecular Line Scan in the Hubble Deep Field North: Constraints on the CO Luminosity Function and the Cosmic H2 Density

We present direct constraints on the CO luminosity function at high redshift and the resulting cosmic evolution of the molecular gas density, $ρ_{\rm H2}$(z), based on a blind molecular line scan in the Hubble Deep Field North (HDF-N) using the IRAM Plateau de Bure Interferometer. Our line scan of the entire 3mm window (79-115 GHz) covers a cosmic volume of ~7000 Mpc$^3$, and redshift ranges z<0.45, 1.01<z<1.89 and z>2. We use the rich multiwavelength and spectroscopic database of the HDF-N to derive some of the best constraints on CO luminosities in high redshift galaxies to date. We combine the blind CO detections in our molecular line scan (presented in a companion paper) with stacked CO limits from galaxies with available spectroscopic redshifts (slit or mask spectroscopy from Keck and grism spectroscopy from HST) to give first blind constraints on high-z CO luminosity functions and the cosmic evolution of the H2 mass density $ρ_{\rm H2}$(z) out to redshifts z~3. A comparison to empirical predictions of $ρ_{\rm H2}$(z) shows that the securely detected sources in our molecular line scan already provide significant contributions to the predicted $ρ_{\rm H2}$(z) in the redshift bins <z>~1.5 and <z>~2.7. Accounting for galaxies with CO luminosities that are not probed by our observations results in cosmic molecular gas densities $ρ_{\rm H2}$(z) that are higher than current predictions. We note however that the current uncertainties (in particular the luminosity limits, number of detections, as well as cosmic volume probed) are significant, a situation that is about to change with the emerging ALMA observatory.

preprint2013arXiv

An ALMA survey of submillimeter galaxies in the Extended Chandra Deep Field South: Source catalog and multiplicity

We present an Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 0 survey of 126 submillimeter sources from the LABOCA ECDFS Submillimeter Survey (LESS). Our 870 micron survey with ALMA (ALESS) has produced maps ~3X deeper and with a beam area ~200X smaller than the original LESS observations, doubling the current number of interferometrically-observed submillimeter sources. The high resolution of these maps allows us to resolve sources that were previously blended and accurately identify the origin of the submillimeter emission. We discuss the creation of the ALESS submillimeter galaxy (SMG) catalog, including the main sample of 99 SMGs and a supplementary sample of 32 SMGs. We find that at least 35% (possibly up to 50%) of the detected LABOCA sources have been resolved into multiple SMGs, and that the average number of SMGs per LESS source increases with LESS flux density. Using the (now precisely known) SMG positions, we empirically test the theoretical expectation for the uncertainty in the single-dish source positions. We also compare our catalog to the previously predicted radio/mid-infrared counterparts, finding that 45% of the ALESS SMGs were missed by this method. Our ~1.6" resolution allows us to measure a size of ~9 kpc X 5 kpc for the rest-frame ~300 um emission region in one resolved SMG, implying a star formation rate surface density of 80 M_sol yr^-1 kpc^-2, and we constrain the emission regions in the remaining SMGs to be <10 kpc. As the first statistically reliable survey of SMGs, this will provide the basis for an unbiased multiwavelength study of SMG properties.

preprint2013arXiv

H2O emission in high-z ultra-luminous infrared galaxies

Using IRAM PdBI we report the detection of H2O in six new lensed ultra-luminous starburst galaxies at high redshift, discovered in the Herschel H-ATLAS survey. The sources are detected either in the 2_{02}-1_{11} or 2_{11}-2_{02} H_2O emission lines with integrated line fluxes ranging from 1.8 to 14 Jy.km/s. The corresponding apparent luminosities are mu x L_H2O ~ 3-12 x 10^8 Lo, where mu is the lensing magnification factor (3 < mu < 12). These results confirm that H2O lines are among the strongest molecular lines in such galaxies, with intensities almost comparable to those of the high-J CO lines, and same profiles and line widths (200-900 km/s) as the latter. With the current sensitivity of PdBI, H2O can therefore easily be detected in high-z lensed galaxies (with F(500um) > 100 mJy) discovered in the Herschel surveys. Correcting the luminosities for lensing amplification, L_H2O is found to have a strong dependence on the IR luminosity, varying as ~L_IR^{1.2}. This relation which needs to be confirmed with better statistics, may indicate a role of radiative (IR) excitation of the H2O lines, and implies that high-z galaxies with L_IR >~ 10^13 Lo tend to be very strong emitters in H2O, that have no equivalent in the local universe.

preprint2013arXiv

HerMES: Candidate High-Redshift Galaxies Discovered with Herschel/SPIRE

We present a method for selecting $z>4$ dusty, star forming galaxies (DSFGs) using Herschel/SPIRE 250/350/500 $μm$ flux densities to search for red sources. We apply this method to 21 deg$^2$ of data from the HerMES survey to produce a catalog of 38 high-$z$ candidates. Follow-up of the first 5 of these sources confirms that this method is efficient at selecting high-$z$ DSFGs, with 4/5 at $z=4.3$ to $6.3$ (and the remaining source at $z=3.4$), and that they are some of the most luminous dusty sources known. Comparison with previous DSFG samples, mostly selected at longer wavelengths (e.g., 850 $μm$) and in single-band surveys, shows that our method is much more efficient at selecting high-$z$ DSFGs, in the sense that a much larger fraction are at $z>3$. Correcting for the selection completeness and purity, we find that the number of bright ($S_{500\,μm} \ge 30$ mJy), red Herschel sources is $3.3 \pm 0.8$ deg$^{-2}$. This is much higher than the number predicted by current models, suggesting that the DSFG population extends to higher redshifts than previously believed. If the shape of the luminosity function for high-$z$ DSFGs is similar to that at $z\sim2$, rest-frame UV based studies may be missing a significant component of the star formation density at $z=4$ to $6$, even after correction for extinction.

preprint2013arXiv

Herschel-ATLAS: a binary HyLIRG pinpointing a cluster of starbursting proto-ellipticals

Panchromatic observations of the best candidate HyLIRG from the widest Herschel extragalactic imaging survey have led to the discovery of at least four intrinsically luminous z=2.41 galaxies across a ~100-kpc region - a cluster of starbursting proto-ellipticals. Via sub-arcsecond interferometric imaging we have measured accurate gas and star-formation surface densities. The two brightest galaxies span ~3 kpc FWHM in submm/radio continuum and CO J=4-3, and double that in CO J=1-0. The broad CO line is due partly to the multitude of constituent galaxies and partly to large rotational velocities in two counter-rotating gas disks -- a scenario predicted to lead to the most intense starbursts, which will therefore come in pairs. The disks have M(dyn) of several x 10^11 Msun, and gas fractions of ~40%. Velocity dispersions are modest so the disks are unstable, potentially on scales commensurate with their radii: these galaxies are undergoing extreme bursts of star formation, not confined to their nuclei, at close to the Eddington limit. Their specific star-formation rates place them ~>5x above the main sequence, which supposedly comprises large gas disks like these. Their high star-formation efficiencies are difficult to reconcile with a simple volumetric star-formation law. N-body and dark matter simulations suggest this system is the progenitor of a B(inary)-type ~10^14.6-Msun cluster.

preprint2013arXiv

The unusual afterglow of the Gamma-Ray Burst 100621A

In order to constrain the broad-band spectral energy distribution of the afterglow of GRB 100621A, dedicated observations were performed in the optical/near-infrared with the 7-channel "Gamma-Ray Burst Optical and Near-infrared Detector" (GROND) at the 2.2m MPG/ESO telescope, in the sub-millimeter band with the large bolometer array LABOCA at APEX, and at radio frequencies with ATCA. Utilizing also Swift X-ray observations, we attempt an interpretation of the observational data within the fireball scenario. The afterglow of GRB 100621A shows a very complex temporal as well as spectral evolution. We identify three different emission components, the most spectacular one causing a sudden intensity jump about one hour after the prompt emission. The spectrum of this component is much steeper than the canonical afterglow. We interpret this component using the prescription of Vlasis et al. (2011) for a two-shell collision after the first shell has been decelerated by the circumburst medium. We use the fireball scenario to derive constraints on the microphysical parameters of the first shell. Long-term energy injection into a narrow jet seems to provide an adequate description. Another noteworthy result is the large ($A_V$ = 3.6 mag) line-of-sight host extinction of the afterglow in an otherwise extremely blue host galaxy.

preprint2012arXiv

A Herschel view of the far-infrared properties of submillimetre galaxies

We study a sample of 61 submillimetre galaxies (SMGs) selected from ground-based surveys, with known spectroscopic redshifts and observed with Herschel as part of the PACS Evolutionary Probe (PEP) and the Herschel Multi-tiered Extragalactic Survey (HerMES) key programmes. We use the broad far-infrared wavelength coverage (100-600um) provided by the combination of PACS and SPIRE observations. Using a power-law temperature distribution model to derive infrared luminosities and dust temperatures, we measure a dust emissivity spectral index for SMGs of beta=2.0+/-0.2. Our results unveil the diversity of the SMG population. Some SMGs exhibit extreme infrared luminosities of ~10^13 Lsun and relatively warm dust components, while others are fainter (~10^12 Lsun) and are biased towards cold dust temperatures. The extreme infrared luminosities of some SMGs (LIR>10^12.7 Lsun, 26/61 systems) imply SFRs of >500Msun yr^-1. Such high SFRs are difficult to reconcile with a secular mode of star formation, and may instead correspond to a merger-driven stage in the evolution of these galaxies. Another observational argument in favour of this scenario is the presence of dust temperatures warmer than that of SMGs of lower luminosities (~40K as opposed to ~25K), consistent with observations of local ULIRGs triggered by major mergers and with results from hydrodynamic simulations of major mergers combined with radiative transfer calculations. Luminous SMGs are also offset from normal star-forming galaxies in the stellar mass-SFR plane, suggesting that they are undergoing starburst events with short duty cycles, compatible with the major merger scenario. On the other hand, a significant fraction of the low infrared luminosity SMGs have cold dust temperatures, are located close to the main sequence of star formation, and thus might be evolving through a secular mode of star formation. [abridged]

preprint2012arXiv

A survey of molecular gas in luminous sub-millimetre galaxies

We present the results from a survey for 12CO emission in 40 luminous sub-millimetre galaxies (SMGs), with 850um fluxes of S850 = 4 - 20 mJy, conducted with the Plateau de Bure Interferometer. We detect 12CO emission in 32 SMGs at z~1.2 - 4.1, including 16 SMGs not previously published. Using multiple 12CO line (J_up =2 - 7) observations, we derive a median spectral line energy distribution for luminous SMGs and use this to estimate a mean gas mass of (5.3 +/- 1.0) \times 10^10 Msun. We report the discovery of a fundamental relationship between 12CO FWHM and 12CO line luminosity in high-redshift starbursts, which we interpret as a natural consequence of the baryon-dominated dynamics within the regions probed by our observations. We use far-infrared luminosities to assess the star-formation efficiency in our SMGs, finding a steepening of the L'CO-LFIR relation as a function of increasing 12CO J_up transition. We derive dynamical masses and molecular gas masses, and use these to determine the redshift evolution of the gas content of SMGs, finding that they do not appear to be significantly more gas rich than less vigorously star-forming galaxies at high redshifts. Finally, we collate X-ray observations, and study the interdependence of gas and dynamical properties of SMGs with their AGN activity and supermassive black hole masses (MBH), finding that SMGs lie significantly below the local M_BH-sigma relation. We conclude that SMGs represent a class of massive, gas-rich ultraluminous galaxies with somewhat heterogeneous properties, ranging from starbursting disc-like systems with L~10^12 L_sun, to the most highly star-forming mergers in the Universe.

preprint2012arXiv

Adaptable Radiative Transfer Innovations for Submillimetre Telescopes (ARTIST) - Dust polarisation module (DustPol)

We present a new publicly available tool (DustPol) aimed to model the polarised thermal dust emission. The module DustPol, which is publicly available, is part of the ARTIST (Adaptable Radiative Transfer Innovations for Submillimetre Telescopes) package, which also offers tools for modelling the polarisation of line emission together with a model library and a Python-based user interface. DustPol can easily manage analytical as well as pre-gridded models to generate synthetic maps of the Stokes I, Q, and U parameters. These maps are stored in FITS format which is straightforwardly read by the data reduction software used, e.g., by the Atacama Large Millimeter Array (ALMA). This turns DustPol into a powerful engine for the prediction of the expected polarisation features of a source observed with ALMA or the Planck satellite as well as for the interpretation of existing submillimetre observations obtained with other telescopes. DustPol allows the parameterisation of the maximum degree of polarisation and we find that, in a prestellar core, if there is depolarisation, this effect should happen at densities of 10^6 cm-3 or larger. We compare a model generated by DustPol with the observational polarisation data of the low-mass Class 0 object NGC 1333 IRAS 4A, finding that the total and the polarised emission are consistent.

preprint2012arXiv

Cool and warm dust emission from M33 (HerM33es)

We study the far-infrared emission from the nearby spiral galaxy M33 in order to investigate the dust physical properties such as the temperature and the luminosity density across the galaxy. Taking advantage of the unique wavelength coverage (100, 160, 250, 350 and 500 micron) of the Herschel Space Observatory and complementing our dataset with Spitzer-IRAC 5.8 and 8 micron and Spitzer-MIPS 24 and 70 micron data, we construct temperature and luminosity density maps by fitting two modified blackbodies of a fixed emissivity index of 1.5. We find that the 'cool' dust grains are heated at temperatures between 11 and 28 K with the lowest temperatures found in the outskirts of the galaxy and the highest ones in the center and in the bright HII regions. The infrared/submillimeter total luminosity (5 - 1000 micron) is estimated to be 1.9x10^9 Lsun. 59% of the total luminosity of the galaxy is produced by the 'cool' dust grains (~15 K) while the rest 41% is produced by 'warm' dust grains (~55 K). The ratio of the cool-to-warm dust luminosity is close to unity (within the computed uncertainties), throughout the galaxy, with the luminosity of the cool dust being slightly enhanced in the center of the galaxy. Decomposing the emission of the dust into two components (one emitted by the diffuse disk of the galaxy and one emitted by the spiral arms) we find that the fraction of the emission in the disk in the mid-infrared (24 micron) is 21%, while it gradually rises up to 57% in the submillimeter (500 micron). We find that the bulk of the luminosity comes from the spiral arm network that produces 70% of the total luminosity of the galaxy with the rest coming from the diffuse dust disk. The 'cool' dust inside the disk is heated at a narrow range of temperatures between 18 and 15 K (going from the center to the outer parts of the galaxy).

preprint2012arXiv

Deep observations of CO line emission from star-forming galaxies in a cluster candidate at z=1.5

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.

preprint2012arXiv

Dust and gas power-spectrum in M33 (HERM33ES)

Power spectra of de-projected images of late-type galaxies in gas and/or dust emission are very useful diagnostics of the dynamics and stability of their interstellar medium. Previous studies have shown that the power spectra can be approximated as two power-laws, a shallow one at large scales (larger than 500 pc) and a steeper one at small scales, with the break between the two corresponding to the line-of-sight thickness of the galaxy disk. We present a thorough analysis of the power spectra of the dust and gas emission at several wavelengths in the nearby galaxy M33. In particular, we use the recently obtained images at five wavelengths by PACS and SPIRE onboard Herschel. The large dynamical range (2-3 dex in scale) of most images allow us to determine clearly the change in slopes from -1.5 to -4, with some variations with wavelength. The break scale is increasing with wavelength, from 100 pc at 24 and 100micron to 350 pc at 500micron, suggesting that the cool dust lies in a thicker disk than the warm dust, may be due to star formation more confined to the plane. The slope at small scale tends to be steeper at longer wavelength, meaning that the warmer dust is more concentrated in clumps. Numerical simulations of an isolated late-type galaxy, rich in gas and with no bulge, like M33, are carried out, in order to better interpret these observed results. Varying the star formation and feedback parameters, it is possible to obtain a range of power-spectra, with two power-law slopes and breaks, which nicely bracket the data. The small-scale power-law is indeed reflecting the 3D behaviour of the gas layer, steepening strongly while the feedback smoothes the structures, by increasing the gas turbulence. M33 appears to correspond to a fiducial model with an SFR of $\sim$ 0.7 Mo/yr, with 10% supernovae energy coupled to the gas kinematics.

preprint2012arXiv

Evidence for Low Extinction in Actively Star Forming Galaxies at z$>$6.5

We present a search for the [CII] 158micron fine structure line (a main cooling line of the interstellar medium) and the underlying far-infrared (FIR) continuum in three high-redshift (6.6<z<8.2) star-forming galaxies using the IRAM Plateau de Bure interferometer. We targeted two Lyman-Alpha-selected galaxies (Lyman-Alpha-Emitters, LAEs) with moderate UV-based star formation rates (SFR~20 M_sun/yr; Himiko at z=6.6 and IOK-1 at z=7.0) and a Gamma Ray Burst (GRB) host galaxy (GRB 090423 at z~8.2). Based on our 3 sigma rest-frame FIR continuum limits, previous (rest-frame) UV continuum measurements and spectral energy distribution (SED) fitting, we rule out SED shapes similar to highly obscured galaxies (e.g. Arp220, M82) and less extreme dust-rich nearby spiral galaxies (e.g. M51) for the LAEs. Conservatively assuming a SED shape typical of local spiral galaxies we derive upper limits for the FIR-based star formation rates (SFRs) of ~70 M_sun/yr, ~50 M_sun/yr and ~40 M_sun/yr for Himiko, IOK-1 and GRB 090423, respectively. For the LAEs these limits are only a factor ~3 higher than the published UV-based SFRs (uncorrected for extinction). This indicates that the dust obscuration in the z>6 LAEs studied here is lower by a factor of a few than what has recently been found in some LAEs at lower redshift (2<z<3.5) with similar UV-based SFRs. A low obscuration in our z>6 LAE sample is consistent with recent rest-frame UV studies of z~7 Lyman-Break-Galaxies (LBGs).

preprint2012arXiv

Ionized nitrogen at high redshift

We present secure [NII[ detections in two mm-bright, strongly lensed objects at high redshift, APM08279+5255 (z=3.911) and MM18423+5938 (z=3.930), using the IRAM Plateau de Bure Interferometer. Due to its ionization energy [NII] is a good tracer of the ionized gas phase in the interstellar medium. The measured fluxes are S([NII])=(4.8+/-0.8) Jy km/s and (7.4+/-0.5) Jy km/s respectively, yielding line luminosities of L([NII]) =(1.8+/-0.3) x 10^9 μ^{-1} Lsun for APM08279+5255 and L([NII]) =(2.8+/-0.2) x 10^9 μ^{-1} Lsun for MM18423+5938. Our high-resolution map of the [NII] and 1 mm continuum emission in MM18423+5938 clearly resolves an Einstein ring in this source, and reveals a velocity gradient in the dynamics of the ionized gas. A comparison of these maps with high-resolution EVLA CO observations enables us to perform the first spatially-resolved study of the dust continuum-to-molecular gas surface brightness (Sigma_{FIR} Sigma_{CO}^N, which can be interpreted as the star formation law) in a high-redshift object. We find a steep relation (N=1.4+/-0.2), consistent with a starbursting environment. We measure a [NII]/FIR luminosity ratio in APM0828+5255 and MM18423+5938 of 9.0 x 10^{-6} and 5.8 x 10^{-6}, respectively. This is in agreement with the decrease of the [NII]/FIR ratio at high FIR luminosities observed in local galaxies.

preprint2012arXiv

Mapping the cold dust temperatures and masses of nearby Kingfish galaxies with Herschel

Taking advantage of the sensitivity and angular resolution of the Herschel Space Observatory at far-infrared and submm wavelengths, we aim to characterize the physical properties of cold dust within nearby galaxies and study the robustness of the parameters we derive using different modified blackbody models. For a pilot subsample of the KINGFISH program, we perform 2 temperature fits of the Spitzer and Herschel photometric data (24 to 500um), with a warm and a cold component, globally and in each resolution element.At global scales, we observe ranges of values for beta_c(0.8 to 2.5) and Tc(19.1 to 25.1K).We compute maps of our parameters with beta fixed or free to test the robustness of the temperature and dust surface density maps we deduce. When the emissivity is fixed, we observe temperature gradients as a function of radius.When the emissivity is fitted as a free parameter, barred galaxies tend to have uniform fitted emissivities.Gathering resolved elements in a Tc-beta_c diagram underlines an anti-correlation between the two parameters.It remains difficult to assess whether the dominant effect is the physics of dust grains, noise, or mixing along the line of sight and in the beam. We finally observe in both cases that the dust column density peaks in central regions of galaxies and bar ends (coinciding with molecular gas density enhancements usually found in these locations).We also quantify how the total dust mass varies with our assumptions about the emissivity index as well as the influence of the wavelength coverage used in the fits. We show that modified blackbody fits using a shallow emissivity (beta_c < 2.0) lead to significantly lower dust masses compared to the beta_c < 2.0 case, with dust masses lower by up to 50% if beta_c=1.5 for instance.The working resolution affects our total dust mass estimates: masses increase from global fits to spatially-resolved fits.

preprint2012arXiv

Measurements of CO redshifts with Z-Spec for lensed submillimeter galaxies discovered in the H-ATLAS survey

We present new observations from Z-Spec, a broadband 185-305 GHz spectrometer, of five sub-millimeter bright lensed sources selected from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) science demonstration phase (SDP) catalog. We construct a redshift finding algorithm using combinations of the signal-to-noise of all the lines falling in the Z-Spec bandpass to determine redshifts with high confidence, even in cases where the signal-to-noise in individual lines is low. We measure the dust continuum in all sources and secure CO redshifts for four out of five (z~1.5-3). In one source, SDP.17, we tentatively identify two independent redshifts and a water line, confirmed at z=2.308. Our sources have properties characteristic of dusty starburst galaxies, with magnification-corrected star formation rates of 10^(2-3) M_sun/yr. Lower limits for the dust masses (~a few 10^8 M_sun) and spatial extents (~1 kpc equivalent radius) are derived from the continuum spectral energy distributions, corresponding to dust temperatures between 54 and 69K. In the LTE approximation, we derive relatively low CO excitation temperatures (< 100 K) and optical depths (tau<1). Performing a non-LTE excitation analysis using RADEX, we find that the CO lines measured by Z-Spec (from J=4->3 to 10->9, depending on the galaxy) localize the best solutions to either a high-temperature / low-density region, or a low-temperature / high-density region near the LTE solution, with the optical depth varying accordingly. Observations of additional CO lines, CO(1-0) in particular, are needed to constrain the non-LTE models.

preprint2012arXiv

Millimeter imaging of submillimeter galaxies in the COSMOS field: Redshift distribution

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 (<z>=3.1+/-0.4) and 870um (<z>=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.

preprint2012arXiv

Quest for COSMOS submillimeter galaxy counterparts using CARMA and VLA: Identifying three high-redshift starburst galaxies

We report on interferometric observations at 1.3 mm at 2"-3" resolution using the Combined Array for Research in Millimeter-wave Astronomy (CARMA). We identify multi-wavelength counterparts of three submillimeter galaxies (SMGs; F(1mm)>5.5 mJy) in the COSMOS field, initially detected with MAMBO and AzTEC bolometers at low, ~10"-30", resolution. All three sources -- AzTEC/C1, Cosbo-3 and Cosbo-8 -- are identified to coincide with positions of 20 cm radio sources. Cosbo-3, however, is not associated with the most likely radio counterpart, closest to the MAMBO source position, but that further away from it. This illustrates the need for intermediate-resolution (~2") mm-observations to identify the correct counterparts of single-dish detected SMGs. All of our three sources become prominent only at NIR wavelengths, and their mm-to-radio flux based redshifts suggest that they lie at redshifts z>~2. As a proof of concept, we show that photometric redshifts can be well determined for SMGs, and we find photometric-redshifts of 5.6+/-1.2, 1.9+0.9(-0.5), and ~4 for AzTEC/C1, Cosbo-3, and Cosbo-8, respectively. Using these we infer that these galaxies have radio-based star formation rates of >~1000 Msol/yr, and IR luminosities of ~10^13 Lsol consistent with properties of high-redshift SMGs. In summary, our sources reflect a variety of SMG properties in terms of redshift and clustering, consistent with the framework that SMGs are progenitors of z~2 and today's passive galaxies.

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.

preprint2011arXiv

A Molecular Einstein Ring Toward the z=3.93 Submillimeter Galaxy MM18423+593

We present high resolution imaging of the low order (J=1 and 2) CO line emission from the z = 3.93 submillimeter galaxy (SMG) MM18423+5938 using the Expanded Very Large Array, and optical and near-IR imaging using the Canada-France-Hawaii Telescope. This SMG with a spectroscopic redshift was thought to be gravitationally lensed given its enormous apparent brightness. We find that the CO emission is consistent with a complete Einstein ring with a major axis diameter of ~ 1.4", indicative of lensing. We have also identified the lens galaxy as a very red elliptical coincident with the geometric center of the ring and estimated its photometric redshift z~1.1. A first estimate of the lens magnification factor is m~12. The luminosity L'_CO(1-0) of the CO(1-0) emission is 2.71+/-0.38 x10^11/m K km s^-1 pc^2, and, adopting the commonly used conversion factor for ULIRGs, the molecular gas mass is M(H_2) = 2.2 x10^11/m M_sol, comparable to unlensed SMGs if corrected by m ~ 12. Our revised estimate of the far-IR luminosity of MM18423+5938 is 2 x10^13/m < L_FIR < 3 x10^14/m L_sol, comparable to that of ULIRGs. Further observations are required to quantify the star formation rate in MM18423+5938 and to constrain the mass model of the lens in more detail.

preprint2011arXiv

Confirmation of the VeLLO L1148-IRS: Star Formation at very low (Column) Density

We report the detection of a compact (of order 5 arcsec; about 1800 AU projected size) CO outflow from L1148-IRS. This confirms that this Spitzer source is physically associated with the nearby (about 325 pc) L1148 dense core. Radiative transfer modeling suggests an internal luminosity of 0.08 to 0.13 L_sun. This validates L1148-IRS as a Very Low Luminosity Object (VeLLO; L < 0.1 L_sun). The L1148 dense core has unusually low densities and column densities for a star-forming core. It is difficult to understand how L1148-IRS might have formed under these conditions. Independent of the exact final mass of this VeLLO (which is likely < 0.24 M_sun), L1148-IRS and similar VeLLOs might hold some clues about the isolated formation of brown dwarfs.

preprint2011arXiv

Molecular Gas in Submillimetre-Faint, Star-Forming Ultraluminous Galaxies at z>1

[abridged] We present interferometric CO observations of twelve z~2 submillimetre-faint, star-forming radio galaxies (SFRGs) which are thought to be ultraluminous infrared galaxies (ULIRGs) possibly dominated by warmer dust (T_dust ~> 40 K) than submillimetre galaxies (SMGs) of similar luminosities. Four other CO-observed SFRGs are included from the literature, and all observations are taken at the Plateau de Bure Interferometer (PdBI) in the compact configuration. Ten of the sixteen SFRGs observed in CO (63%) are detected at >4sigma with a mean inferred molecular gas mass of ~2*10^10 M_sun. SFRGs trend slightly above the local ULIRG L_FIR-L'_CO relation. Since SFRGs are about two times fainter in radio luminosity but exhibit similar CO luminosities to SMGs, this suggests SFRGs are slightly more efficient star formers than SMGs at the same redshifts. SFRGs also have a narrow mean CO line width, 320+-80km/s. SFRGs bridge the gap between properties of very luminous >5*10^12 L_sun SMGs and those of local ULIRGs and are consistent with intermediate stage major mergers. We suspect that more moderate-luminosity SMGs, not yet surveyed in CO, would show similar molecular gas properties to SFRGs. The AGN fraction of SFRGs is consistent with SMGs and is estimated to be 0.3+-0.1, suggesting that SFRGs are observed near the peak phase of star formation activity and not in a later, post-SMG enhanced AGN phase. This CO survey of SFRGs serves as a pilot project for the much more extensive survey of Herschel and SCUBA-2 selected sources which only partially overlap with SMGs. Better constraints on CO properties of a diverse high-z ULIRG population are needed from ALMA to determine the evolutionary origin of extreme starbursts, and what role ULIRGs serve in catalyzing the formation of massive stellar systems in the early Universe.

preprint2011arXiv

Observation of H2O in a strongly lensed Herschel-ATLAS source at z=2.3

The Herschel survey, H-ATLAS, with its large areal coverage, has recently discovered a number of bright, strongly lensed high-z submillimeter galaxies. The strong magnification makes it possible to study molecular species other than CO, which are otherwise difficult to observe in high-z galaxies. Among the lensed galaxies already identified by H-ATLAS, the source J090302.9-014127B (SDP.17b) at z = 2.305 is remarkable due to its excitation conditions and a tentative detection of the H2O 202-111 emission line (Lupu et al. 2010). We report observations of this line in SDP.17b using the IRAM interferometer equipped with its new 277- 371GHz receivers. The H2O line is detected at a redshift of z = 2.3049+/-0.0006, with a flux of 7.8+/-0.5 Jy km s-1 and a FWHM of 250+/-60 km s-1. The new flux is 2.4 times weaker than the previous tentative detection, although both remain marginally consistent within 1.6-sigma. The intrinsic line luminosity and ratio of H2O(202-111)/CO8-7 seem comparable with those of the nearby starburst/enshrouded-AGN Mrk 231, suggesting that SDP.17b could also host a luminous AGN. The detection of a strong H2O 202-111 line in SDP.17b implies an efficient excitation mechanism of the water levels that must occur in very dense and warm interstellar gas.

preprint2011arXiv

Radio studies of galaxy formation: Dense Gas History of the Universe

Line and continuum studies at centimeter through submillimeter wavelengths address probe deep into the earliest, most active and dust obscured phases of galaxy formation, and reveal the molecular and cool atomic gas. We summarize the techniques of radio astronomy to perform these studies, then review the progress on radio studies of galaxy formation. The dominant work over the last decade has focused on massive, luminous starburst galaxies (submm galaxies and AGN host galaxies). The far infrared luminosities are ~ 1e13 Lsun, implying star formation rates, SFR > 1e3 Msun/year. Molecular gas reservoirs are found with masses: M(H_2) > 1e10 (alpha/0.8}) Msun. The CO excitation in these luminous systems is much higher than in low redshift spiral galaxies. Imaging of the gas distribution and dynamics suggests strongly interacting and merging galaxies, indicating gravitationally induced, short duration (~ 1e7 year) starbursts. These systems correspond to a major star formation episode in massive galaxies in proto-clusters at intermediate to high redshift. Recently, radio observations have probed the more typical star forming galaxy population (SFR ~ 100 Msun/year), during the peak epoch of Universal star formation (z ~ 1.5 to 2.5). These observations reveal massive gas reservoirs without hyper-starbursts, and show that active star formation occurs over a wide range in galaxy stellar mass. The conditions in this gas are comparable to those found in the Milky Way disk. A key result is that the peak epoch of star formation in the Universe also corresponds to an epoch when the baryon content of star forming galaxies was dominated by molecular gas, not stars. We consider the possibility of tracing out the dense gas history of the Universe, and perform initial, admittedly gross, calculations. ABRIDGED

preprint2011arXiv

Redshift evolution of the 1.4 GHz volume averaged radio luminosity function in clusters of galaxies

By cross-correlating large samples of galaxy clusters with publicly available radio source catalogs, we construct the volume-averaged radio luminosity function (RLF) in clusters of galaxies, and investigate its dependence on cluster redshift and mass. In addition, we determine the correlation between the cluster mass and the radio luminosity of the brightest source within 50 kpc from the cluster center. We use two cluster samples: the optically selected maxBCG cluster catalog and a composite sample of X-ray selected clusters. The radio data come from the VLA NVSS and FIRST surveys. We use scaling relations to estimate cluster masses and radii to get robust estimates of cluster volumes. We determine the projected radial distribution of sources, for which we find no dependence on luminosity or cluster mass. Background and foreground sources are statistically accounted for, and we account for confusion of radio sources by adaptively degrading the resolution of the radio source surveys. We determine the redshift evolution of the RLF under the assumption that its overall shape does not change with redshift. Our results are consistent with a pure luminosity evolution of the RLF in the range 0.1 < z < 0.3 from the optical cluster sample. The X-ray sample extends to higher redshift and yields results also consistent with a pure luminosity evolution. We find no direct evidence of a dependence of the RLF on cluster mass from the present data, although the data are consistent with the most luminous sources only being found in high-mass systems.

preprint2011arXiv

Spitzer Imaging of Herschel-ATLAS Gravitationally Lensed Submillimeter Sources

We present physical properties of two submillimeter selected gravitationally lensed sources, identified in the Herschel Astrophysical Terahertz Large Area Survey. These submillimeter galaxies (SMGs) have flux densities > 100 mJy at 500 um, but are not visible in existing optical imaging. We fit light profiles to each component of the lensing systems in Spitzer IRAC 3.6 and 4.5 um data and successfully disentangle the foreground lens from the background source in each case, providing important constraints on the spectral energy distributions (SEDs) of the background SMG at rest-frame optical-near-infrared wavelengths. The SED fits show that these two SMGs have high dust obscuration with Av ~4 to 5 and star formation rates of ~100 M_sun/yr. They have low gas fractions and low dynamical masses compared to 850 um selected galaxies.

preprint2010arXiv

100 mum and 160 mum emission as resolved star-formation rate estimators in M33 (HERM33ES)

Over the past few years several studies have provided estimates of the SFR (star-formation rate) or the total infrared luminosity from just one infrared band. However these relations are generally derived for entire galaxies, which are known to contain a large scale diffuse emission that is not necessarily related to the latest star-formation episode. We provide new relations to estimate the SFR from resolved star-forming regions at 100 mum and 160 mum. We select individual star-forming regions in the nearby (840 kpc) galaxy M33. We estimate the SFR combining the emission in Halpha and at 24 mum to calibrate the emission at 100 mum and 160 mum as SFR estimators, as mapped with PACS/Herschel. The data are obtained in the framework of the HERM33ES open time key project. There is less emission in the HII regions at 160 mum than at 100 mum. Over a dynamic range of almost 2 dex in Sigma(SFR) we find that the 100 mum emission is a nearly linear estimator of the SFR, whereas that at 160 mum is slightly superlinear. The behaviour of individual star-forming regions is surprisingly similar to that of entire galaxies. At high Sigma(SFR), star formation drives the dust temperature, whereas uncertainties and variations in radiation-transfer and dust-heated processes dominate at low Sigma(SFR). Detailed modelling of both galaxies and individual star forming regions will be needed to interpret similarities and differences between the two and assess the fraction of diffuse emission in galaxies.

preprint2010arXiv

Cool gas and dust in M33: Results from the Herschel M33 extended survey (HERM33ES)

We present an analysis of the first space-based far-IR-submm observations of M 33, which measure the emission from the cool dust and resolve the giant molecular cloud complexes. With roughly half-solar abundances, M33 is a first step towards young low-metallicity galaxies where the submm may be able to provide an alternative to CO mapping to measure their H$_2$ content. In this Letter, we measure the dust emission cross-section $σ$ using SPIRE and recent CO and \HI\ observations; a variation in $σ$ is present from a near-solar neighborhood cross-section to about half-solar with the maximum being south of the nucleus. Calculating the total H column density from the measured dust temperature and cross-section, and then subtracting the \HI\ column, yields a morphology similar to that observed in CO. The H$_2$/\HI\ mass ratio decreases from about unity to well below 10% and is about 15% averaged over the optical disk. The single most important observation to reduce the potentially large systematic errors is to complete the CO mapping of M 33.

preprint2010arXiv

Discovery of an extremely bright submillimeter galaxy at z=3.93

Serendipitously we have discovered a rare, bright submillimeter galaxy (SMG) with a flux density of 30 +/- 2 mJy at lambda=1.2mm, using MAMBO2 at the IRAM 30-meter millimeter telescope. Although no optical counterpart is known for MM18423+5938, we were able to measure the redshift z=3.92960 +/- 0.00013 from the detection of CO lines using the IRAM Eight MIxer Receiver (EMIR). In addition, by collecting all available photometric data in the far-infrared and radio to constrain its spectral energy distribution, we derive the FIR luminosity 4.8 10^14/m Lsol and mass 6.0 10^9/m Msol for its dust, allowing for a magnification factor m caused by a probable gravitational lens. The corresponding star-formation rate is 8.3 10^4/m Msol/yr. The detection of three lines of the CO rotational ladder, and a significant upper limit for a fourth CO line, allow us to estimate an H2 mass of between 1.9 10^11/m Msol and 1.1 10^12/m Msol. The two lines CI(3p1-3p0) and CI(3p2-3p1) were clearly detected and yield a [CI]/[H2] number abundance between 1.4 10^-5 and 8.0 10^-5. Upper limits are presented for emission lines of HCN, HCO^+, HNC, H_2O and other molecules observed. The moderate excitation of the CO lines is indicative of an extended starburst, and excludes the dominance of an AGN in heating this high-redshift SMG.

preprint2010arXiv

Herschel ATLAS: The cosmic star formation history of quasar host galaxies

We present a derivation of the star formation rate per comoving volume of quasar host galaxies, derived from stacking analyses of far-infrared to mm-wave photometry of quasars with redshifts 0<z<6 and absolute I-band magnitudes -22>I_AB>-32. We use the science demonstration observations of the first ~16 deg^2 from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) in which there are 240 quasars from the Sloan Digital Sky Survey (SDSS) and a further 171 from the 2dF-SDSS LRG and QSO (2SLAQ) survey. We supplement this data with a compilation of data from IRAS, ISO, Spitzer, SCUBA and MAMBO. H-ATLAS alone statistically detects the quasars in its survey area at >5sigma at 250, 350 and 500um. From the compilation as a whole we find striking evidence of downsizing in quasar host galaxy formation: low-luminosity quasars with absolute magnitudes in the range -22>I_AB>-24 have a comoving star formation rate (derived from 100um rest-frame luminosities) peaking between redshifts of 1 and 2, while high-luminosity quasars with I_AB<-26 have a maximum contribution to the star formation density at z~3. The volume-averaged star formation rate of -22>I_AB>-24 quasars evolves as (1+z)^{2.3 +/- 0.7} at z<2, but the evolution at higher luminosities is much faster reaching (1+z)^{10 +/- 1} at -26>I_AB>-28. We tentatively interpret this as a combination of a declining major merger rate with time and gas consumption reducing fuel for both black hole accretion and star formation.

preprint2010arXiv

Herschel-ATLAS: blazars in the SDP field

To investigate the poorly constrained sub-mm counts and spectral properties of blazars we searched for these in the Herschel-ATLAS (H-ATLAS) science demostration phase (SDP) survey catalog. We cross-matched 500$μ$m sources brighter than 50 mJy with the FIRST radio catalogue. We found two blazars, both previously known. Our study is among the first blind blazar searches at sub-mm wavelengths, i.e., in the spectral regime where little is still known about the blazar SEDs, but where the synchrotron peak of the most luminous blazars is expected to occur. Our early results are consistent with educated extrapolations of lower frequency counts and question indications of substantial spectral curvature downwards and of spectral upturns at mm wavelengths. One of the two blazars is identified with a Fermi/LAT $γ$-ray source and a WMAP source. The physical parameters of the two blazars are briefly discussed.These observations demonstrate that the H-ATLAS survey will provide key information about the physics of blazars and their contribution to sub-mm counts.

preprint2010arXiv

Herschel-ATLAS: Statistical Properties of Galactic Cirrus in the GAMA-9 Hour Science Demonstration Phase Field

We study the Spectral Energy Distribution (SED) and the power spectrum of Galactic cirrus emission observed in the 14 deg^2 Science Demonstration Phase field of the Herschel-ATLAS using Herschel and IRAS data from 100 to 500 um. We compare the SPIRE 250, 350 and 500um maps with IRAS 100um emission, binned in 6' pixels. We assume a modified black-body SED with dust emissivity parameter beta (F ~ lambda^(-beta)) and a single dust temperature T_d, and find that the dust temperature and emissivity index varies over the science demonstration field as 10< T_rm < 25 K and 1 < beta< 4. The latter values are somewhat higher than the range of beta often quoted in the literature (1< beta< 2). We estimate the mean values of these parameters to be T_d=19.0 +/- 2.4 K and beta = 1.4 +/- 0.4. In regions of bright cirrus emission, we find that the dust has similar temperatures with T_d = 18.0 +/- 2.5 K, and similar values of beta, ranging from 1.4 +- 0.5 to 1.9+/- 0.5. We show that T_d and beta associated with diffuse cirrus emission are anti-correlated and can be described by the relationship: beta(T_d) = NT_d^alpha with [N=116+/-38, alpha=-1.4+/1 0.1]. The strong correlation found in this analysis is not just limited to high density clumps of cirrus emission as seen in previous studies, but is also seen in diffuse cirrus in low density regions. To provide an independent measure of $T_{\rm d}$ and $β$, we obtain the angular power spectrum of the cirrus emission in the {\it IRAS} and SPIRE maps, which is consistent with a power spectrum of the form P(k)=P_0(k/k_0)^gamma where gamma = â H R2.6+/-m 0.2 for scales of 50-200' in the SPIRE maps. The cirrus rms fluctuation amplitude at angular scales of 100' is consistent with a modified blackbody SED with T_d = 20.1+/- 0.9 K and beta = 1.3+/- 0.2, in agreement with the values obtained above.

preprint2010arXiv

Herschel/PACS far-infrared photometry of two z>4 quasars

We present Herschel far-infrared (FIR) observations of two sub-mm bright quasars at high redshift: SDSS J1148+5251 (z=6.42) and BR 1202-0725 (z=4.69) obtained with the PACS instrument. Both objects are detected in the PACS photometric bands. The Herschel measurements provide additional data points that constrain the FIR spectral energy distributions (SEDs) of both sources, and they emphasise a broad range of dust temperatures in these objects. For lambda_rest ~< 20mu, the two SEDs are very similar to the average SEDs of quasars at low redshift. In the FIR, however, both quasars show excess emission compared to low-z QSO templates, most likely from cold dust powered by vigorous star formation in the QSO host galaxies. For SDSS J1148+5251 we detect another object at 160mu with a distance of ~10 arcseconds from the QSO. Although no physical connection between the quasar and this object can be shown with the available data, it could potentially confuse low-resolution measurements, thus resulting in an overestimate of the FIR luminosity of the z=6.42 quasar.

preprint2010arXiv

High-resolution CO and radio imaging of z~2 ULIRGs: extended CO structures, and implications for the universal star formation law

We present high spatial resolution (0.4", ~3.5 kpc) PdBI interferometric data on three ultra-luminous infrared galaxies (ULIRGs) at z~2: two submillimetre galaxies and one submillimetre faint star forming radio galaxy. The three galaxies have been ro- bustly detected in CO rotational transitions, either 12CO(J=4-3) or 12CO(J=3-2), allowing their sizes and gas masses to be accurately constrained. These are the highest spatial resolution observations observed to date (by a factor of ~2) for intermediate-excitation CO emission in z~2 ULIRGs. The galaxies appear extended over several resolution elements, having a mean radius of 3.7 kpc. High-resolution (0.3") combined MERLIN-VLA observations of their radio continua allow an analysis of the star formation behaviour of these galaxies, on comparable spatial scales to that of the CO observations. This 'matched beam' approach sheds light on the spatial distribution of both molecular gas and star formation, and we can therefore calculate accurate star formation rates and gas surface densities: this allows us to place the three systems in the context of a Kennicutt-Schmidt (KS)-style star formation law. We find a difference in size between the CO and radio emission regions, and as such we suggest that using the spatial extent of the CO emission region to estimate the surface density of star formation may lead to error. This size difference also causes the star formation efficiencies within systems to vary by up to a factor of 5. We also find, with our new accurate sizes, that SMGs lie significantly above the KS relation, indicating that stars are formed more efficiently in these extreme systems than in other high-z star forming galaxies.

preprint2010arXiv

Identification of two bright z>3 submillimeter galaxy candidates in the COSMOS field

We present high-resolution interferometric Submillimeter Array (SMA) imaging at 890 microns (~2" resolution) of two millimeter selected galaxies -- MMJ100015+021549 and MMJ100047+021021 -- discovered with the Max-Planck Millimeter Bolometer (MAMBO) on the IRAM 30 m telescope and also detected with Bolocam on the CSO, in the COSMOS field. The first source is significantly detected at the ~11 sigma level, while the second source is tentatively detected at the ~4 sigma level, leading to a positional accuracy of ~0.2-0.3". MM100015+021549 is identified with a faint radio and K-band source. MMJ100047+021021 shows no radio emission and is tentatively identified with a very faint K-band peak which lies at ~1.2" from a clumpy optical source. The submillimeter-to-radio flux ratio for MM100015+021549 yields a redshift of ~4.8, consistent with the redshift implied by the UV-to-submillimeter photometry, z~3.0-5.0. We find evidence for warm dust in this source with an infrared luminosity in the range ~0.9-2.5x10^{13} L_sun, supporting the increasing evidence for a population luminous submillimeter galaxies at z>3. Finally, the lack of photometric data for MMJ100047+021021 does not allow us to investigate its properties in detail, however its submillimeter-to-radio ratio implies z>3.5.

preprint2010arXiv

Ionization near-zones associated with quasars at z ~ 6

We analyze the size evolution of HII regions around 27 quasars between z=5.7 to 6.4 ('quasar near-zones' or NZ). We include more sources than previous studies, and we use more accurate redshifts for the host galaxies, with 8 CO molecular line redshifts and 9 MgII redshifts. We confirm the trend for an increase in NZ size with decreasing redshift, with the luminosity normalized proper size evolving as: R_{NZ,corrected} = (7.4 \pm 0.3) - (8.0 \pm 1.1) \times (z-6) Mpc. While derivation of the absolute neutral fraction remains difficult with this technique, the evolution of the NZ sizes suggests a decrease in the neutral fraction of intergalactic hydrogen by a factor ~ 9.4 from z=6.4 to 5.7, in its simplest interpretation. Alternatively, recent numerical simulations suggest that this rapid increase in near-zone size from z=6.4 to 5.7 is due to the rapid increase in the background photo-ionization rate at the end of the percolation or overlap phase, when the average mean free path of ionizing photons increases dramatically. In either case, the results are consistent with the idea that z ~ 6 to 7 corresponds to the tail end of cosmic reionization. The scatter in the normalized NZ sizes is larger than expected simply from measurement errors, and likely reflects intrinsic differences in the quasars or their environments. We find that the near-zone sizes increase with quasar UV luminosity, as expected for photo-ionization dominated by quasar radiation.

preprint2010arXiv

Most Submillimetre Galaxies are Major Mergers

We analyse subarcsecond resolution interferometric CO line data for twelve sub-millimetre-luminous (S850um > 5mJy) galaxies with redshifts between 1 and 3, presenting new data for four of them. Morphologically and kinematically most of the twelve systems appear to be major mergers. Five of them are well-resolved binary systems, and seven are compact or poorly resolved. Of the four binary systems for which mass measurements for both separate components can be made, all have mass ratios of 1:3 or closer. Furthermore, comparison of the ratio of compact to binary systems with that observed in local ULIRGs indicates that at least a significant fraction of the compact SMGs must also be late-stage mergers. In addition, the dynamical and gas masses we derive are most consistent with the lower end of the range of stellar masses published for these systems, favouring cosmological models in which SMGs result from mergers. These results all point to the same conclusion, that likely most of the bright SMGs with L_IR > 5x10e12L_sun are major mergers.

preprint2010arXiv

Non-parametric modeling of the intra-cluster gas using APEX-SZ bolometer imaging data

We demonstrate the usability of mm-wavelength imaging data obtained from the APEX-SZ bolometer array to derive the radial temperature profile of the hot intra-cluster gas out to radius r_500 and beyond. The goal is to study the physical properties of the intra-cluster gas by using a non-parametric de-projection method that is, aside from the assumption of spherical symmetry, free from modeling bias. We use publicly available X-ray imaging data from the XMM-Newton observatory and our Sunyaev-Zel'dovich Effect (SZE) imaging data from the APEX-SZ experiment at 150 GHz to de-project the density and temperature profiles for the relaxed cluster Abell 2204. We derive the gas density, temperature and entropy profiles assuming spherical symmetry, and obtain the total mass profile under the assumption of hydrostatic equilibrium. For comparison with X-ray spectroscopic temperature models, a re-analysis of the recent Chandra observation is done with the latest calibration updates. Using the non-parametric modeling we demonstrate a decrease of gas temperature in the cluster outskirts, and also measure the gas entropy profile. These results are obtained for the first time independently of X-ray spectroscopy, using SZE and X-ray imaging data. The contribution of the SZE systematic uncertainties in measuring T_e at large radii is shown to be small compared to the Chandra systematic spectroscopic errors. The upper limit on M_200 derived from the non-parametric method is consistent with the NFW model prediction from weak lensing analysis.

preprint2010arXiv

PACS and SPIRE photometer maps of M33: First results of the Herschel M33 extended survey (HERM33ES)

Within the framework of the HERM33ES key project, we are studying the star forming interstellar medium in the nearby, metal-poor spiral galaxy M33, exploiting the high resolution and sensitivity of Herschel. We use PACS and SPIRE maps at 100, 160, 250, 350, and 500 micron wavelength, to study the variation of the spectral energy distributions (SEDs) with galacto-centric distance. Detailed SED modeling is performed using azimuthally averaged fluxes in elliptical rings of 2 kpc width, out to 8 kpc galacto-centric distance. Simple isothermal and two-component grey body models, with fixed dust emissivity index, are fitted to the SEDs between 24 and 500 micron using also MIPS/Spitzer data, to derive first estimates of the dust physical conditions. The far-infrared and submillimeter maps reveal the branched, knotted spiral structure of M33. An underlying diffuse disk is seen in all SPIRE maps (250-500 micron). Two component fits to the SEDs agree better than isothermal models with the observed, total and radially averaged flux densities. The two component model, with beta fixed at 1.5, best fits the global and the radial SEDs. The cold dust component clearly dominates; the relative mass of the warm component is less than 0.3% for all the fits. The temperature of the warm component is not well constrained and is found to be about 60K plus/minus 10K. The temperature of the cold component drops significantly from about 24K in the inner 2 kpc radius to 13K beyond 6 kpc radial distance, for the best fitting model. The gas-to-dust ratio for beta=1.5, averaged over the galaxy, is higher than the solar value by a factor of 1.5 and is roughly in agreement with the subsolar metallicity of M33.

preprint2010arXiv

Photon Dominated Regions in NGC 3603

Aims: We aim at deriving the excitation conditions of the interstellar gas as well as the local FUV intensities in the molecular cloud surrounding NGC 3603 to get a coherent picture of how the gas is energized by the central stars. Methods: The NANTEN2-4m submillimeter antenna is used to map the [CI] 1-0, 2-1 and CO 4-3, 7-6 lines in a 2' x 2' region around the young OB cluster NGC 3603 YC. These data are combined with C18O 2-1 data, HIRES-processed IRAS 60 and 100 micron maps of the FIR continuum, and Spitzer/IRAC maps. Results: The NANTEN2 observations show the presence of two molecular clumps located south-east and south-west of the cluster and confirm the overall structure already found by previous CS and C18O observations. We find a slight position offset of the peak intensity of CO and [CI], and the atomic carbon appears to be further extended compared to the molecular material. We used the HIRES far-infrared dust data to derive a map of the FUV field heating the dust. We constrain the FUV field to values of χ= 3 - 6 \times 10^3 in units of the Draine field across the clouds. Approximately 0.2 to 0.3 % of the total FUV energy is re-emitted in the [CII] 158 μm cooling line observed by ISO. Applying LTE and escape probability calculations, we derive temperatures (TMM1 = 43 K, TMM2 = 47 K), column densities (N(MM1) = 0.9 \times 10^22 cm^-2, N(MM2) = 2.5 \times 10^22 cm^-2) and densities (n(MM1) = 3 \times 10^3 cm^-3, n(MM2) = 10^3 -10^4 cm^-3) for the two observed molecular clumps MM1 and MM2. Conclusions: The cluster is strongly interacting with the ambient molecular cloud, governing its structure and physical conditions. A stability analysis shows the existence of gravitationally collapsing gas clumps which should lead to star formation. Embedded IR sources have already been observed in the outskirts of the molecular cloud and seem to support our conclusions.

preprint2010arXiv

Radio observations of the cool gas, dust, and star formation in the first galaxies

We summarize cm through submm observations of the host galaxies of z ~ 6 quasars. These observations reveal the cool molecular gas (the fuel for star formation), the warm dust (heated by star formation), the fine structure line emission (tracing the CNM and PDRs), and the synchrotron emission. Our results imply active star formation in ~ 30% of the host galaxies, with star formation rates ~ 10^3 M_sun/year, and molecular gas masses ~ 10^10 M_sun. Imaging of the [CII] emission from the most distant quasar reveals a 'maximal starburst disk' on a scale ~ 1.5 kpc. Gas dynamical studies suggest a departure of these galaxies from the low-z M_{BH} -- M_{bulge} relation, with the black holes being, on average, 15 times more massive than expected. Overall, we are witnessing the co-eval formation of massive galaxies and supermassive black holes within 1 Gyr of the Big Bang.

preprint2010arXiv

Star Forming Dense Cloud Cores in the TeV γ-ray SNR RX J1713.7-3946

RX J1713.7-3946 is one of the TeV γ-ray supernova remnants (SNRs) emitting synchrotron X rays. The SNR is associated with molecular gas located at ~1 kpc. We made new molecular observations toward the dense cloud cores, peaks A, C and D, in the SNR in the 12CO(J=2-1) and 13CO(J=2-1) transitions at angular resolution of 90". The most intense core in 13CO, peak C, was also mapped in the 12CO(J=4-3) transition at angular resolution of 38". Peak C shows strong signs of active star formation including bipolar outflow and a far-infrared protostellar source and has a steep gradient with a r^{-2.2$\pm$0.4} variation in the average density within radius r. Peak C and the other dense cloud cores are rim-brightened in synchrotron X rays, suggesting that the dense cloud cores are embedded within or on the outer boundary of the SNR shell. This confirms the earlier suggestion that the X rays are physically associated with the molecular gas (Fukui et al. 2003). We present a scenario where the densest molecular core, peak C, survived against the blast wave and is now embedded within the SNR. Numerical simulations of the shock-cloud interaction indicate that a dense clump can indeed survive shock erosion, since shock propagation speed is stalled in the dense clump. Additionally, the shock-cloud interaction induces turbulence and magnetic field amplification around the dense clump that may facilitate particle acceleration in the lower-density inter-clump space leading to the enhanced synchrotron X rays around dense cores.

preprint2010arXiv

The Detection of a Population of Submillimeter-Bright, Strongly-Lensed Galaxies

Gravitational lensing is a powerful astrophysical and cosmological probe and is particularly valuable at submillimeter wavelengths for the study of the statistical and individual properties of dusty starforming galaxies. However the identification of gravitational lenses is often time-intensive, involving the sifting of large volumes of imaging or spectroscopic data to find few candidates. We used early data from the Herschel Astrophysical Terahertz Large Area Survey to demonstrate that wide-area submillimeter surveys can simply and easily detect strong gravitational lensing events, with close to 100% efficiency.

preprint2010arXiv

The VLA-COSMOS Survey. IV. Deep Data and Joint Catalog

In the context of the VLA-COSMOS Deep project additional VLA A array observations at 1.4 GHz were obtained for the central degree of the COSMOS field and combined with the existing data from the VLA-COSMOS Large project. A newly constructed Deep mosaic with a resolution of 2.5" was used to search for sources down to 4 sigma with 1 sigma ~ 12 microJy/beam in the central 50'x50'. This new catalog is combined with the catalog from the Large project (obtained at 1.5"x1.4" resolution) to construct a new Joint catalog. All sources listed in the new Joint catalog have peak flux densities of >5 sigma at 1.5" and/or 2.5" resolution to account for the fact that a significant fraction of sources at these low flux levels are expected to be slighty resolved at 1.5" resolution. All properties listed in the Joint catalog such as peak flux density, integrated flux density and source size are determined in the 2.5" resolution Deep image. In addition, the Joint catalog contains 43 newly identified multi-component sources.

preprint2009arXiv

A LABOCA survey of the Extended Chandra Deep Field South - submillimeter properties of near-IR selected galaxies

Using the 330hr ESO-MPG 870-micron survey of the Extended Chandra Deep Field South (ECDF-S) obtained with the Large Apex BOlometer CAmera (LABOCA) on the Atacama Pathfinder EXperiment (APEX), we have carried out a stacking analysis at submillimeter (submm) wavelengths of a sample of 8266 near-infra-red (near-IR) selected (K_vega <= 20) galaxies, including 893 BzK galaxies, 1253 extremely red objects (EROs) and 737 distant red galaxies (DRGs), selected from the Multi-wavelength Survey by Yale-Chile (MUSYC). We measure average 870-micron fluxes of 0.20+-0.01mJy (20.0sigma), 0.45+-0.04mJy (11.3sigma), 0.42+-0.03mJy (14.0sigma), and 0.41+-0.04mJy (10.3sigma) for the K_vega <= 20, BzK, ERO and DRG samples, respectively. For the BzK, ERO and DRG subsamples, which overlap to some degree and are like to be at z ~ 1-2, this implies an average far-IR luminosity of ~2-6x10^{11} Lsolar and star formation rate of ~40-100Msolar. Splitting the BzK galaxies up into star-forming (sBzK) and passive (pBzK) galaxies, the former is significantly detected (0.48+-0.04mJy, 12.0sigma) while the latter is only marginally detected (0.27+-0.10mJy, 2.7sigma), thus confirming that the sBzK/pBzK criteria do isolate obscured, star forming and truly passive galaxies. The K_vega <= 20 galaxies are found to contribute with 6.6+-0.3Jy deg^{-2} (~15%) to the 870-micron extragalactic background light (EBL). sBzK and pBzK galaxies contribute 1.7+-0.2Jy deg^{-2} (~4%) and 0.2+-0.1 Jy deg^{-2} (< 0.5%) to the EBL. [Abridged]

preprint2009arXiv

A Search for Neutral Carbon towards two z=4.05 Submillimetre Galaxies, GN20 and GN20.2

Using the IRAM Plateau-de-Bure Interferometer (PdBI) we have searched for the upper fine structure line of neutral carbon (CI(2->1), nu_rest = 809 GHz) and 12CO(J=7->6) (nu_rest=806GHz) towards the submillimetre galaxies (SMGs) GN20 (SMM J123711.9+622212, z = 4.055) and GN20.2 (SMM J123708.8+622202, z = 4.051). The far-infrared (FIR) continuum is detected at 8sigma significance in GN20, with a flux density of S_1.8mm = 1.9+-0.2mJy, while no continuum is detected in GN20.2. Both sources are statistically undetected in both CI and CO(7-6) lines; we derive line luminosity limits for both CI and CO of L'<~2x10^10K km s^-1 pc^2. Assuming carbon excitation temperatures of T_ex = 30K (the galaxies' measured dust temperatures), we infer CI mass limits of M_CI < 5.4x10^6 Msun (GN20) and M_CI < 6.8x10^6 Msun (GN20.2). The derived CI abundance limits are <1.8x10^-5 for GN20 and <3.8x10^-5 for GN20.2 implying that the systems have Milky Way level neutral carbon enrichment (X[CI]/X[H2]) or lower, similar to high-redshift carbon-detected systems (at 5x10^-5) but about 50 times less than the neutral carbon enrichment of local starburst galaxies. Observations of GN20 and GN20.2 in high-resolution MERLIN+VLA radio maps of GOODS-N are used to further constrain the sizes and locations of active regions. We conclude that the physical gas properties of young rapidly evolving systems like GN20 and GN20.2 are likely significantly different than starburst/ULIRG environments in the local Universe yet similar to z~2 SMGs. Unless gravitationally amplified examples can be found, observations of galaxies like GN20 will require the order of magnitude increase in sensitivity of the Atacama Large Millimetre Array (ALMA) to constrain their CI and high-J CO content, despite the fact that they are the brightest systems at z~4.

preprint2009arXiv

Photometric monitoring of the young star Par 1724 in Orion

We report new photometric observations of the 200000 year old naked weak-line run-away T Tauri star Par 1724, located north of the Trapezium cluster in Orion. We observed in the broad band filters B, V, R, and I using the 90cm Dutch telescope on La Silla, the 80cm Wendelstein telescope, and a 25cm telescope of the University Observatory Jena in Grossschwabhausen near Jena. The photometric data in V and R are consistent with a 5.7 day rotation period due to spots, as observed before between 1960ies and 2000. Also, for the first time, we present evidence for a long-term 9 or 17.5 year cycle in photometric data (V band) of such a young star, a cycle similar to that to of the Sun and other active stars.

preprint2009arXiv

Search for Cold Debris Disks around M-dwarfs. II

Although 70 % of the stars in the Galaxy are M-dwarfs, thermal emission searches for cold debris disks have been conducted mostly for A-type and solar-type stars. We report on new lambda=1.2 mm continuum observations of thirty M-dwarfs, using the MAMBO-2 bolometer array camera at the IRAM 30m telescope. For a statistical analysis, we combine these data with our prior SCUBA and MAMBO-2 observations of 20 other M-dwarfs. Our total sample divides in M-dwarfs in moving groups, with relatively young ages, and in nearby M-dwarfs with unknown ages. Only one cold debris disk (GJ842.2) was detected significantly. We compare the implied disk abundance constraints with those found in two comparable submillimeter surveys of 10 to 190 Myr old A- and FGK-type stars. For the 19 youngest (ages less than 200 Myr) M-dwarfs in our sample, we derive a cold disk fraction of 5.3^{+10.5}_{-5.0} %, compared to 15 +/-11.5 % for FGK-stars and 22^{+33}_{-20} % for A-stars. Hence, for this age group, there is an apparent trend of fewer cold disks for later stellar types, i.e., lower star masses. Although its statistical significance is marginal, this trend is strengthened by the deeper sensitivity of observations in the M-dwarf sample. We derive a cold disk fraction of < 10 % for the older (likely a few Gyr) M-dwarfs in our sample. Finally, although inconclusively related to a debris disk, we present the complex millimeter structure found around the position of the M1.5 dwarf GJ526 in our sample.

preprint2009arXiv

The Herschel ATLAS

The Herschel ATLAS is the largest open-time key project that will be carried out on the Herschel Space Observatory. It will survey 510 square degrees of the extragalactic sky, four times larger than all the other Herschel surveys combined, in five far-infrared and submillimetre bands. We describe the survey, the complementary multi-wavelength datasets that will be combined with the Herschel data, and the six major science programmes we are undertaking. Using new models based on a previous submillimetre survey of galaxies, we present predictions of the properties of the ATLAS sources in other wavebands.

preprint2009arXiv

The LABOCA Survey of the Extended Chandra Deep Field South

We present a sensitive 870 micron survey of the Extended Chandra Deep Field South (ECDFS) using LABOCA on the APEX telescope. The LABOCA ECDFS Submillimetre Survey (LESS) covers the full 30' x 30' field size of the ECDFS and has a uniform noise level of 1.2 mJy/beam. LESS is thus the largest contiguous deep submillimetre survey undertaken to date. The noise properties of our map show clear evidence that we are beginning to be affected by confusion noise. We present a catalog of 126 SMGs detected with a significance level above 3.7 sigma. The ECDFS exhibits a deficit of bright SMGs relative to previously studied blank fields but not of normal star-forming galaxies that dominate the extragalactic background light (EBL). This is in line with the underdensities observed for optically defined high redshift source populations in the ECDFS (BzKs, DRGs,optically bright AGN and massive K-band selected galaxies). The differential source counts in the full field are well described by a power law with a slope of alpha=-3.2, comparable to the results from other fields. We show that the shape of the source counts is not uniform across the field. The integrated 870 micron flux densities of our source-count models account for >65% of the estimated EBL from COBE measurements. We have investigated the clustering of SMGs in the ECDFS by means of a two-point correlation function and find evidence for strong clustering on angular scales <1'. Assuming a power law dependence for the correlation function and a typical redshift distribution for the SMGs we derive a spatial correlation length of r_0=13+/-6 h^-1 Mpc.

preprint2008arXiv

A strong optical flare before the rising afterglow of GRB 080129

We report on GROND observations of a 40 sec duration (rest-frame) optical flare from GRB 080129 at redshift 4.349. The rise- and decay time follow a power law with indices +12 and -8, respectively, inconsistent with a reverse shock and a factor 10$^5$ faster than variability caused by ISM interaction. While optical flares have been seen in the past (e.g. GRB 990123, 041219B, 060111B and 080319B), for the first time, our observations not only resolve the optical flare into sub-components, but also provide a spectral energy distribution from the optical to the near-infrared once every minute. The delay of the flare relative to the GRB, its spectral energy distribution as well as the ratio of pulse widths suggest it to arise from residual collisions in GRB outflows \cite{liw08}.If this interpretation is correct and can be supported by more detailed modelling or observation in further GRBs, the delay measurement provides an independent, determination of the Lorentz factor of the outflow.

preprint2008arXiv

Submillimeter Galaxies at z~2: Evidence for Major Mergers & Constraints on Lifetimes, IMF and CO-H2 Conversion Factor

We report sub-arcsecond resolution IRAM PdBI millimeter CO interferometry of four z~2 submillimeter galaxies (SMGs), and sensitive CO (3-2) flux limits toward three z~2 UV-/optically selected star forming galaxies. The new data reveal for the first time spatially resolved CO gas kinematics in the observed SMGs. Two of the SMGs show double or multiple morphologies, with complex, disturbed gas motions. The other two SMGs exhibit CO velocity gradients of ~500 km/s across 0.2 arcsec (1.6 kpc) diameter regions, suggesting that the star forming gas is in compact, rotating disks. Our data provide compelling evidence that these SMGs represent extreme, short-lived 'maximum' star forming events in highly dissipative mergers of gas rich galaxies. The resulting high mass surface and volume densities of SMGs are similar to those of compact quiescent galaxies in the same redshift range, and much higher than those in local spheroids. From the ratio of the comoving volume densities of SMGs and quiescent galaxies in the same mass and redshift ranges, and from the comparison of gas exhaustion time scales and stellar ages, we estimate that the SMG phase duration is about 100 Myrs. Our analysis of SMGs and optically/UV selected high redshift star forming galaxies supports a 'universal' Chabrier IMF as being valid over the star forming history of these galaxies. We find that the 12CO luminosity to total gas mass conversion factors at z~2-3 are probably similar to those assumed at z~0. The implied gas fractions in our sample galaxies range from 20 to 50%.

preprint2005arXiv

A robust sample of submillimetre galaxies: constraints on the prevalence of dusty, high-redshift starbursts

The modest significance of most sources detected in current submm surveys can potentially compromise some analyses due to the inclusion of spurious sources in catalogues typically selected at >3.0-3.5sigma. Here, we develop and apply a dual-survey extraction technique to SCUBA and MAMBO images of the Lockman Hole. Cut above 5sigma, our catalogue of SMGs is more robust than previous samples, with a reduced likelihood of real, but faint SMGs (beneath and around the confusion limit) entering via superposition with noise. The effective flux limit of the survey is well matched to our deep 1.4-GHz image. The former is sensitive to luminous, dusty galaxies at extreme redshifts whilst the latter probes the z<3 regime. A high fraction of our robust SMGs (~80%) have radio counterparts which, given the ~10% contamination by spurious sources, suggests that very distant SMGs (z>>3) are unlikely to make up more than ~10% of the bright SMG population. This implies that almost all of the bright SMG population is amenable to study via the deepest current radio imaging. We use these radio counterparts to provide an empirical calibration of the positional uncertainty in SMG catalogues. We then go on to outline the acquisition of redshifts for radio-identified SMGs, from sample selection in the submm, to counterpart selection in the radio and optical/IR, to slit placement on spectrograph masks. We determine a median of z=2.05 +/- 0.41 from a sample of six secure redshifts for unambigious radio-identified submm sources and z=2.14 +/- 0.27 when we include submm sources with multiple radio counterparts and/or less reliable redshifts. These figures are consistent with previous estimates, suggesting that our knowledge of the median redshift of bright SMGs population has not been biased by the low significance of the source catalogues employed.