Source author record

F. Walter

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

138works
5topics
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

138 published item(s)

preprint2024arXiv

NOEMA reveals the true nature of luminous red JWST z>10 galaxy candidates

The first year of JWST has revealed a surprisingly large number of luminous galaxy candidates beyond $z>10$. While some galaxies have already been spectroscopically confirmed, there is mounting evidence that a subsample of the candidates with particularly red inferred UV colours are, in fact, lower redshift contaminants. These interlopers are often found to be `HST-dark' or `optically faint' galaxies at $z\sim2-6$, a population that is key to improving our understanding of dust-obscured star formation throughout cosmic time. This paper demonstrates the complementarity of ground-based mm-interferometry and JWST infrared imaging to unveil the true nature of red 1.5-2.0 $μ\rm{m}$ dropouts that have been selected as ultra-high-redshift galaxy candidates. We present NOEMA Polyfix follow-up observations of four JWST red 1.5-2.0 $μ\rm{m}$ dropouts selected by Yan et al., 2023 as ultra-high-redshift candidates in the PEARLS-IDF field. The new NOEMA observations constrain the rest-frame far-infrared continuum emission and efficiently discriminate between intermediate- and high-redshift solutions. We report $>10σ$ NOEMA continuum detections of all our target galaxies at observed frequencies of $ν= 236$ and $252\ \rm{GHz}$, with FIR slopes indicating a redshift of $z<5$. We modelled their optical-to-FIR spectral energy distribution (SED) with multiple SED codes, finding that they are not $z>10$ galaxies but dust-obscured, massive star-forming galaxies at $z\sim 2-4$ instead. The contribution to the cosmic star formation rate density (CSFRD) of such sources is not negligible at $z\simeq 3.5$ ($ϕ\gtrsim(1.9-4.4)\times10^{-3}\ \rm{cMpc}^{-3}$; or $>3-6\%$ of the total CSFRD), in line with previous studies of optically faint and sub-millimeter galaxies. This approach opens up a new window onto obscured star formation at intermediate redshifts [abridged].

preprint2022arXiv

A dusty compact object bridging galaxies and quasars at cosmic dawn

Understanding how super-massive black holes form and grow in the early Universe has become a major challenge since the discovery of luminous quasars only 700 million years after the Big Bang. Simulations indicate an evolutionary sequence of dust-reddened quasars emerging from heavily dust-obscured starbursts that then transition to unobscured luminous quasars by expelling gas and dust. Although the last phase has been identified out to a redshift of 7.6, a transitioning quasar has not been found at similar redshifts owing to their faintness at optical and near-infrared wavelengths. Here we report observations of an ultraviolet compact object, GNz7q, associated with a dust-enshrouded starburst at a redshift of z=7.1899+/-0.0005. The host galaxy is more luminous in dust emission than any other known object at this epoch, forming 1,600 solar masses of stars per year within a central radius of 480 parsec. A red point source in the far-ultraviolet is identified in deep, high-resolution imaging and slitless spectroscopy. GNz7q is extremely faint in X-rays, which indicates the emergence of a uniquely ultraviolet compact star-forming region or a Compton-thick super-Eddington black-hole accretion disk at the dusty starburst core. In the latter case, the observed properties are consistent with predictions from cosmological simulations and suggest that GNz7q is an antecedent to unobscured luminous quasars at later epochs.

preprint2022arXiv

Dense Gas History of the Universe: from ASPECS to the ngVLA

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

preprint2022arXiv

The recurrent nova V3890~Sgr: a near-infrared and optical study of the red giant component and its environment

We present an analysis of the red giant component of the recurrent nova V3890 Sgr, using data obtained before and after its 2019 eruption. Its effective temperature is $T_{\rm eff}=3050\pm$200 K for $\log{g}=0.7$, although there are modest changes in $T_{\rm eff}$. There is an overabundance of both carbon ($0.20\pm0.05$~dex) and sodium ($1.0\pm0.3$~dex) relative to their solar values, possibly the result of ejecta from the 1990 nova eruption being entrained into the red giant photosphere. We find $^{12}$C/$^{13}$C $=25\pm2$, a value similar to that found in red giants in other recurrent novae. The interpretation of the quiescent spectrum in the 5--38$\,mu$m region requires the presence of photospheric SiO absorption and cool ($\sim400$~K) dust in the red giant environment. The spectrum in the region of the Na{\sc i} D lines is complex, and includes at least six interstellar components, together with likely evidence for interaction between ejecta from the 2019 eruption and material accumulated in the plane of the binary. Three recurrent novae with giant secondaries have been shown to have environments with different dust content, but photospheres with similar $^{12}$C/$^{13}$C ratios. The SiO fundamental bands most likely have a photospheric origin in the all three stars.

preprint2021arXiv

ExoClock project II: A large-scale integrated study with 180 updated exoplanet ephemerides

The ExoClock project is an inclusive, integrated, and interactive platform that was developed to monitor the ephemerides of the Ariel targets to increase the mission efficiency. The project makes the best use of all available resources, i.e., observations from ground telescopes, mid-time values from the literature and finally, observations from space instruments. Currently, the ExoClock network includes 280 participants with telescopes capable of observing 85\% of the currently known Ariel candidate targets. This work includes the results of $\sim$1600 observations obtained up to the 31st of December 2020 from the ExoClock network. These data in combination with $\sim$2350 mid-time values collected from the literature are used to update the ephemerides of 180 planets. The analysis shows that 40\% of the updated ephemerides will have an impact on future scheduling as either they have a significantly improved precision, or they have revealed biases in the old ephemerides. With the new observations, the observing coverage and rate for half of the planets in the sample has been doubled or more. Finally, from a population perspective, we identify that the differences in the 2028 predictions between the old and the new ephemerides have an STD that is double what is expected from gaussian uncertainties. These findings have implications for planning future observations, where we will need to account for drifts potentially greater than the prediction uncertainties. The updated ephemerides are open and accessible to the wider exoplanet community both from our Open Science Framework (OSF) repository and our website.

preprint2020arXiv

An ALMA survey of the S2CLS UDS field: Optically invisible submillimetre galaxies

We analyse a robust sample of 30 near-infrared-faint (K>25.3, 5 sigma) submillimetre galaxies selected across a 0.96 deg^2 field, to investigate their properties and the cause of their lack of detectable optical/near-infrared emission. Our analysis exploits precise identifications based on ALMA 870um continuum imaging, combined with the very deep near-infrared imaging from the UKIDSS-UDS survey. We estimate that K>25.3 submillimetre galaxies represent 15+/-2 per cent of the total population brighter than S870=3.6mJy, with an expected surface density of ~450/deg^2 above S870>1mJy. As such they pose a source of contamination in surveys for both high-redshift "quiescent" galaxies and very-high-redshift Lyman-break galaxies. We show that these K-faint submillimetre galaxies are simply the tail of the broader submillimetre population, with comparable dust and stellar masses to K<25.3 mag submillimetre galaxies, but lying at significantly higher redshifts (z=3.44+/-0.06 versus z=2.36+/-0.11) and having higher dust attenuation (Av=5.2+/-0.3 versus Av=2.9+/-0.1). We investigate the origin of the strong dust attenuation and find indications that these K-faint galaxies have smaller dust continuum sizes than the K<25.3 galaxies, as measured by ALMA, which suggests their high attenuation is related to their compact sizes. We find a correlation of dust attenuation with star-formation rate surface density (Sigma_SFR), with the K-faint submillimetre galaxies representing the higher-Sigma_SFR and highest-Av galaxies. The concentrated, intense star-formation activity in these systems is likely to be associated with the formation of spheroids in compact galaxies at high redshifts, but as a result of their high obscuration these are completely missed in UV, optical and even near-infrared surveys.

preprint2020arXiv

An Ultra-deep Multi-band VLA Survey of the Faint Radio Sky (COSMOS-XS): Source Catalog and Number Counts

We present ultra-deep, matched-resolution Karl G. Jansky Very Large Array (VLA) observations at 10 and $3$ GHz in the COSMOS field: the COSMOS-XS survey. The final 10 and $3$ GHz images cover $\sim16\rm{arcmin}^{2}$ and $\sim180\rm{arcmin}^{2}$ and reach median rms values of $0.41μ\rm{Jy\,beam}^{-1}$ and $0.53μ\rm{Jy\,beam}^{-1}$, respectively. Both images have an angular resolution of $\sim 2.0''$. To fully account for the spectral shape and resolution variations across the broad bands, we image all data with a multi-scale, multi-frequency synthesis algorithm. We present source catalogs for the 10 and $3$ GHz image with 91 and 1498 sources, respectively, above a peak brightness threshold of $5σ$. We present source counts with completeness corrections included that are computed via Monte Carlo simulations. Our corrected radio counts at $3$ GHz with direct detections down to $\sim2.8μ$Jy are consistent within the uncertainties with other results at 3 and 1.4 GHz, but extend to fainter flux densities than previous direct detections. The ultra-faint $3$ GHz number counts are found to exceed the counts predicted by the semi-empirical radio sky simulations developed in the framework of the SKA Simulated Skies project, consistent with previous P(D) analyses. Our measured source counts suggest a steeper luminosity function evolution for these faint star-forming sources. The semi-empirical Tiered Radio Extragalactic Continuum Simulation (T-RECS) predicts this steeper evolution and is in better agreement with our results. The $10$ GHz radio number counts also agree with the counts predicted by the T-RECS simulation within the expected variations from cosmic variance. In summary, the multi-band, matched-resolution COSMOS-XS survey in the well-studied COSMOS field provides a high-resolution view of the ultra-faint radio sky that can help guide next generation radio facilities.

preprint2019arXiv

Infrared spectroscopy of the recent outburst in V1047 Cen (Nova Centauri 2005)

Fourteen years after its eruption as a classical nova (CN), V1047 Cen (Nova Cen 2005) began an unusual re-brightening in 2019 April. The amplitude of the brightening suggests that this is a dwarf nova (DN) eruption in a CN system. Very few CNe have had DN eruptions within decades of the main CN outburst. The 14 years separating the CN and DN eruptions of V1047 Cen is the shortest of all instances recorded thus far. Explaining this rapid succession of CN and DN outbursts in V1047 Cen may be challenging within the framework of standard theories for DN outbursts. Following a CN eruption, the mass accretion rate is believed to remain high $(\dot{M}\sim10^{-8}$M$_\odot$yr$^{-1})$ for a few centuries, due to the irradiation of the secondary star by the still-hot surface of the white dwarf. Thus a DN eruption is not expected to occur during this high mass accretion phase as DN outbursts, which result from thermal instabilities in the accretion disk, and arise during a regime of low mass accretion rate $(\dot{M}\sim10^{-10}$M$_\odot$yr$^{-1})$. Here we present near-infrared spectroscopy to show that the present outburst is most likely a DN eruption, and discuss the possible reasons for its early occurrence. Even if the present re-brightening is later shown to be due to a cause other than a DN outburst, the present study provides invaluable documentation of this unusual event.

preprint2019arXiv

Modeling Dust and Starlight in Galaxies Observed by Spitzer and Herschel: The KINGFISH Sample

Dust and starlight are modeled for the KINGFISH project galaxies. With data from 3.6 micron to 500 micron, models are strongly constrained. For each pixel in each galaxy we estimate (1) dust surface density; (2) q_PAH, the dust mass fraction in PAHs; (3) distribution of starlight intensities heating the dust; (4) luminosity emitted by the dust; and (5) dust luminosity from regions with high starlight intensity. The models successfully reproduce both global and resolved spectral energy distributions. We provide well-resolved maps for the dust properties. As in previous studies, we find q_PAH to be an increasing function of metallicity, above a threshold Z/Z_sol approx 0.15. Dust masses are obtained by summing the dust mass over the map pixels; these "resolved" dust masses are consistent with the masses inferred from model fits to the global photometry. The global dust-to-gas ratios obtained from this study correlate with galaxy metallicities. Systems with Z/Z_sol > 0.5 have most of their refractory elements locked up in dust, whereas when Z/Z_sol < 0.3 most of these elements tend to remain in the gas phase. Within galaxies, we find that q_PAH is suppressed in regions with unusually warm dust with nu L_nu(70 um) > 0.4L_dust. With knowledge of one long-wavelength flux density ratio (e.g., f_{160}/f_{500}), the minimum starlight intensity heating the dust (U_min) can be estimated to within ~50%. For the adopted dust model, dust masses can be estimated to within ~0.07 dex accuracy using the 500 micron luminosity nu L_nu(500) alone. There are additional systematic errors arising from the choice of dust model, but these are hard to estimate. These calibrated prescriptions may be useful for studies of high-redshift galaxies.

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

Comparing [CII], HI, and CO dynamics of nearby galaxies

The HI and CO components of the interstellar medium (ISM) are usually used to derive the dynamical mass M_dyn of nearby galaxies. Both components become too faint to be used as a tracer in observations of high-redshift galaxies. In those cases, the 158 $μ$m line of atomic carbon [CII] may be the only way to derive M_dyn. As the distribution and kinematics of the ISM tracer affects the determination of M_dyn, it is important to quantify the relative distributions of HI, CO and [CII]. HI and CO are well-characterised observationally, however, for [CII] only very few measurements exist. Here we compare observations of CO, HI, and [CII] emission of a sample of nearby galaxies, drawn from the HERACLES, THINGS and KINGFISH surveys. We find that within R_25, the average [CII] exponential radial profile is slightly shallower than that of the CO, but much steeper than the HI distribution. This is also reflected in the integrated spectrum ("global profile"), where the [CII] spectrum looks more like that of the CO than that of the HI. For one galaxy, a spectrally resolved comparison of integrated spectra was possible; other comparisons were limited by the intrinsic line-widths of the galaxies and the coarse velocity resolution of the [CII] data. Using high-spectral-resolution SOFIA [CII] data of a number of star forming regions in two nearby galaxies, we find that their [CII] linewidths agree better with those of the CO than the HI. As the radial extent of a given ISM tracer is a key input in deriving M_dyn from spatially unresolved data, we conclude that the relevant length-scale to use in determining M_dyn based on [CII] data, is that of the well-characterised CO distribution. This length scale is similar to that of the optical disk.

preprint2016arXiv

He II $λ$4686 emission from the massive binary system in $η$ Car: constraints to the orbital elements and the nature of the periodic minima

η Carinae is an extremely massive binary system in which rapid spectrum variations occur near periastron. Most notably, near periastron the He II $λ4686$ line increases rapidly in strength, drops to a minimum value, then increases briefly before fading away. To understand this behavior, we conducted an intense spectroscopic monitoring of the He II $λ4686$ emission line across the 2014.6 periastron passage using ground- and space-based telescopes. Comparison with previous data confirmed the overall repeatability of EW(He II $λ4686$), the line radial velocities, and the timing of the minimum, though the strongest peak was systematically larger in 2014 than in 2009 by 26%. The EW(He II $λ4686$) variations, combined with other measurements, yield an orbital period $2022.7\pm0.3$ d. The observed variability of the EW(He II $λ4686$) was reproduced by a model in which the line flux primarily arises at the apex of the wind-wind collision and scales inversely with the square of the stellar separation, if we account for the excess emission as the companion star plunges into the hot inner layers of the primary's atmosphere, and including absorption from the disturbed primary wind between the source and the observer. This model constrains the orbital inclination to $135^\circ$-$153^\circ$, and the longitude of periastron to $234^\circ$-$252^\circ$. It also suggests that periastron passage occurred on $T_0 = 2456874.4\pm1.3$ d. Our model also reproduced EW(He II $λ4686$) variations from a polar view of the primary star as determined from the observed He II $λ4686$ emission scattered off the Homunculus nebula.

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

New Constraints on the Molecular Gas in the Prototypical HyLIRGs BRI1202-0725 & BRI1335-0417

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

preprint2016arXiv

The Pan-STARRS1 distant z>5.6 quasar survey: more than 100 quasars within the first Gyr of the universe

Luminous quasars at z>5.6 can be studied in detail with the current generation of telescopes and provide us with unique information on the first gigayear of the universe. Thus far these studies have been statistically limited by the number of quasars known at these redshifts. Such quasars are rare and therefore wide-field surveys are required to identify them and multiwavelength data are needed to separate them efficiently from their main contaminants, the far more numerous cool dwarfs. In this paper, we update and extend the selection for z~6 quasars presented in Banados et al. (2014) using the Pan-STARRS1 (PS1) survey. We present the PS1 distant quasar sample, which currently consists of 124 quasars in the redshift range 5.6<z<6.7 that satisfy our selection criteria. Seventy-seven of these quasars have been discovered with PS1, and 63 of them are newly identified in this paper. We present composite spectra of the PS1 distant quasar sample. This sample spans a factor of ~20 in luminosity and shows a variety of emission line properties. The number of quasars at z>5.6 presented in this work almost double the quasars previously known at these redshifts, marking a transition phase from studies of individual sources to statistical studies of the high-redshift quasar population, which was impossible with earlier, smaller samples.

preprint2015arXiv

B.R.N.O. Contributions #39 Times of minima

This paper presents 1463 times of minima for 455 objects acquired by 46 members and cooperating observers of the Variable Star and Exoplanet Section of the Czech Astronomical Society (B.R.N.O. Observing project). Observations were carried out between October 2013 - September 2014. Some neglected southern eclipsing binaries and newly discovered stars by the observers of project B.R.N.O. are included in the list.

preprint2015arXiv

Bright [CII] 158$μ$m emission in a quasar host galaxy at $z=6.54$

The [CII] 158$μ$m fine-structure line is known to trace regions of active star formation and is the main coolant of the cold, neutral atomic medium. In this \textit{Letter}, we report a strong detection of the [CII] line in the host galaxy of the brightest quasar known at $z>6.5$, the Pan-STARRS1 selected quasar PSO J036.5078+03.0498 (hereafter P036+03), using the IRAM NOEMA millimeter interferometer. Its [CII] and total far-infrared luminosities are $(5.8 \pm 0.7) \times 10^9 \,L_\odot$ and $(7.6\pm1.5) \times 10^{12}\,L_\odot$, respectively. This results in a $L_{[CII]} /L_{TIR}$ ratio of $\sim 0.8\times 10^{-3}$, which is at the high end for those found for active galaxies, though it is lower than the average found in typical main sequence galaxies at $z\sim 0$. We also report a tentative additional line which we identify as a blended emission from the $3_{22} - 3_{13}$ and $5_{23} - 4_{32}$ H$_2$O transitions. If confirmed, this would be the most distant detection of water emission to date. P036+03 rivals the current prototypical luminous J1148+5251 quasar at $z=6.42$, in both rest-frame UV and [CII] luminosities. Given its brightness and because it is visible from both hemispheres (unlike J1148+5251), P036+03 has the potential of becoming an important laboratory for the study of star formation and of the interstellar medium only $\sim 800\,$Myr after the Big Bang.

preprint2015arXiv

Bright [CII] and dust emission in three z>6.6 quasar host galaxies observed by ALMA

We present ALMA detections of the [CII] 158 micron emission line and the underlying far-infrared continuum of three quasars at 6.6<z<6.9 selected from the VIKING survey. The [CII] line fluxes range between 1.6-3.4 Jy km/s ([CII] luminosities ~(1.9-3.9)x10^9 L_sun). We measure continuum flux densities of 0.56-3.29 mJy around 158 micron (rest-frame), with implied far-infrared luminosities between (0.6-7.5)x10^12 L_sun and dust masses M_d=(0.7-24)x10^8 M_sun. In one quasar we derive a dust temperature of 30^+12_-9 K from the continuum slope, below the canonical value of 47 K. Assuming that the [CII] and continuum emission are powered by star formation, we find star-formation rates from 100-1600 M_sun/yr based on local scaling relations. The L_[CII]/L_FIR ratios in the quasar hosts span a wide range from (0.3-4.6)x10^-3, including one quasar with a ratio that is consistent with local star-forming galaxies. We find that the strength of the L_[CII] and 158 micron continuum emission in z>~6 quasar hosts correlate with the quasar's bolometric luminosity. In one quasar, the [CII] line is significantly redshifted by ~1700 km/s with respect to the MgII broad emission line. Comparing to values in the literature, we find that, on average, the MgII is blueshifted by 480 km/s (with a standard deviation of 630 km/s) with respect to the host galaxy redshift, i.e. one of our quasars is an extreme outlier. Through modeling we can rule out a flat rotation curve for our brightest [CII] emitter. Finally, we find that the ratio of black hole mass to host galaxy (dynamical) mass is higher by a factor 3-4 (with significant scatter) than local relations.

preprint2015arXiv

CO excitation of normal star forming galaxies out to z=1.5 as regulated by the properties of their interstellar medium

We investigate the CO excitation of normal star forming disk galaxies at z=1.5 using IRAM PdBI observations of the CO[2-1], CO[3-2] and CO[5-4] transitions for 4 galaxies, including VLA observations of CO[1-0] for 3 of them, with the aim of constraining the average state of H2 gas. Exploiting prior knowledge of the velocity range, spatial extent and size of the CO emission we measure reliable line fluxes with S/N>4-7 for individual transitions. While the average CO Spectral Line Energy Distribution (SLED) has a sub-thermal excitation similar to the Milky Way (MW) up to CO[3-2], we show that the average CO[5-4] emission is four times stronger than assuming MW excitation. This demonstrates the presence of an additional component of more excited, denser and possibly warmer molecular gas. The ratio of CO[5-4] to lower-J CO emission is lower than in local (U)LIRGs and high-redshift SMGs, correlating closely with the average intensity of the radiation field <U> and with the star formation surface density, but not with the SF efficiency (SFE). The CO[5-4] luminosity correlates linearly with LIR over 4 orders of magnitudes, with z=1.5 BzK galaxies following the same trend as local spirals and (U)LIRGs and high redshift star bursting SMGs. The CO[5-4] luminosity is thus empirically related to the dense gas, and might be a more convenient way to probe it than standard high--density tracers that are much fainter than CO. We see excitation variations among our sample galaxies, that can be linked to their evolutionary state and clumpiness in optical rest frame images. In one galaxy we see spatially resolved excitation variations, where the more highly excited part of the galaxy corresponds to the location of massive SF clumps. This supports to models that suggest that giant clumps are the main source of the high excitation CO emission in high redshift disk-like galaxies.

preprint2015arXiv

Constraining the radio-loud fraction of quasars at z>5.5

Radio-loud Active Galactic Nuclei at z~2-4 are typically located in dense environments and their host galaxies are among the most massive systems at those redshifts, providing key insights for galaxy evolution. Finding radio-loud quasars at the highest accessible redshifts (z~6) is important to study their properties and environments at even earlier cosmic time. They would also serve as background sources for radio surveys intended to study the intergalactic medium beyond the epoch of reionization in HI 21 cm absorption. Currently, only five radio-loud ($R=f_{ν,5{\rm GHz}}/f_{ν,4400Å}>10$) quasars are known at z~6. In this paper we search for 5.5 < z < 7.2 quasars by cross-matching the optical Pan-STARRS1 and radio FIRST surveys. The radio information allows identification of quasars missed by typical color-based selections. While we find no good 6.4 < z <7.2 quasar candidates at the sensitivities of these surveys, we discover two new radio-loud quasars at z~6. Furthermore, we identify two additional z~6 radio-loud quasars which were not previously known to be radio-loud, nearly doubling the current z~6 sample. We show the importance of having infrared photometry for z>5.5 quasars to robustly classify them as radio-quiet or radio-loud. Based on this, we reclassify the quasar J0203+0012 (z=5.72), previously considered radio-loud, to be radio-quiet. Using the available data in the literature, we constrain the radio-loud fraction of quasars at z~6, using the Kaplan--Meier estimator, to be $8.1^{+5.0}_{-3.2}\%$. This result is consistent with there being no evolution of the radio-loud fraction with redshift, in contrast to what has been suggested by some studies at lower redshifts.

preprint2015arXiv

Flash-ionization of pre-existing circumstellar material around Nova Oph 2015

We have obtained daily high resolution Echelle spectroscopy of Nova Oph 2015 during its initial evolution. It reveals the presence of pre-existing circumstellar material around the nova, which is best interpreted as the wind of an evolved companion. On earliest observations, the emission line profiles of Nova Oph 2015 displayed a very narrow emission component (FWHM~60 km/s), recombining over a time scale of 5 days and showing constant low velocity (RV_sun= -4.5 km/s), that we interpret as coming from the wind of the secondary recombining after the ionization from the initial UV-flash of the nova. The underlying broad component underwent a marked reduction in both FWHM and width at zero intensity (the latter declining from 4000 to 2000 km/s in ten days) while increasing by 6 times in flux, that we believe is the result of the high velocity ejecta of the nova being slowed down while trying to expand within the surrounding wind of the companion. Novae with evolved secondaries are very rare in the Galaxy, amounting 3% of the total according to recent estimates. Among them Nova Oph 2015 is perhaps unique in having displayed a long rise to maximum brightness and a slow decline from it, a FeII-type classification (contrary to prevailing He/N) and a probable sub-giant luminosity class for the secondary (instead of the giant (e.g. RS Oph) or supergiant (e.g. V407 Cyg) class for the others).

preprint2015arXiv

HI and CO Velocity Dispersions in Nearby Galaxies

We analyze the velocity dispersions of individual HI and CO profiles in a number of nearby galaxies from the high-resolution HERACLES CO and THINGS HI surveys. Focusing on regions with bright CO emission, we find a CO dispersion value: 7.3 $\pm$ 1.7 km/s. The corresponding HI dispersion is 11.7 $\pm$ 2.3 km/s, yielding a mean HI/CO dispersion ratio of 1.4 $\pm$ 0.2, independent of radius. We find that the CO velocity dispersion increases towards lower peak fluxes. This is consistent with previous work where we showed that when using spectra averaged ("stacked") over large areas, larger values for the CO dispersion are found, and a lower dispersion ratio: 1.0 $\pm$ 0.2. The stacking method is more sensitive to low-level diffuse emission, whereas individual profiles trace narrow-line, GMC-dominated, bright emission. These results provide further evidence that disk galaxies contain not only a thin, low velocity dispersion, high density CO disk that is dominated by GMCs, but also a fainter, higher dispersion, diffuse disk component.

preprint2015arXiv

Spatially Resolved Spitzer-IRS Spectral Maps of the Superwind in M82

We have mapped the superwind/halo region of the nearby starburst galaxy M82 in the mid-infrared with $Spitzer-IRS$. The spectral regions covered include the H$_2 S(1)-S(3)$, [NeII], [NeIII] emission lines and PAH features. We estimate the total warm H$_2$ mass and the kinetic energy of the outflowing warm molecular gas to be between $M_{warm}\sim5-17\times10^6$ M$_{\odot}$ and $E_{K}\sim6-20\times10^{53}$ erg. Using the ratios of the 6.2, 7.7 and 11.3 micron PAH features in the IRS spectra, we are able to estimate the average size and ionization state of the small grains in the superwind. There are large variations in the PAH flux ratios throughout the outflow. The 11.3/7.7 and the 6.2/7.7 PAH ratios both vary by more than a factor of five across the wind region. The Northern part of the wind has a significant population of PAH's with smaller 6.2/7.7 ratios than either the starburst disk or the Southern wind, indicating that on average, PAH emitters are larger and more ionized. The warm molecular gas to PAH flux ratios (H$_2/PAH$) are enhanced in the outflow by factors of 10-100 as compared to the starburst disk. This enhancement in the H$_2/PAH$ ratio does not seem to follow the ionization of the atomic gas (as measured with the [NeIII]/[NeII] line flux ratio) in the outflow. This suggests that much of the warm H$_2$ in the outflow is excited by shocks. The observed H$_2$ line intensities can be reproduced with low velocity shocks ($v < 40$ km s$^{-1}$) driven into moderately dense molecular gas ($10^2 <n_H < 10^4$ cm$^{-3}$) entrained in the outflow.

preprint2015arXiv

The Herschel Comprehensive (U)LIRG Emission Survey (HerCULES): CO Ladders, fine structure lines, and neutral gas cooling

(Ultra) Luminous Infrared Galaxies ((U)LIRGs) are objects characterized by their extreme infrared (8-1000 $μ$m) luminosities ($L_{LIRG}>10^{11} $L$_\odot$ and $L_{ULIRG}>10^{12}$ L$_\odot$). The Herschel Comprehensive ULIRG Emission Survey (HerCULES; PI van der Werf) presents a representative flux-limited sample of 29 (U)LIRGs that spans the full luminosity range of these objects (10$^{11}\leq L_\odot \geq10^{13}$). With the \emph{Herschel Space Observatory}, we observe [CII] 157 $μ$m, [OI] 63 $μ$m, and [OI] 145 $μ$m line emission with PACS, CO J=4-3 through J=13-12, [CI] 370 $μ$m, and [CI] 609 $μ$m with SPIRE, and low-J CO transitions with ground-based telescopes. The CO ladders of the sample are separated into three classes based on their excitation level. In 13 of the galaxies, the [OI] 63 $μ$m emission line is self absorbed. Comparing the CO excitation to the IRAS 60/100 $μ$m ratio and to far infrared luminosity, we find that the CO excitation is more correlated to the far infrared colors. We present cooling budgets for the galaxies and find fine-structure line flux deficits in the [CII], [SiII], [OI], and [CI] lines in the objects with the highest far IR fluxes, but do not observe this for CO $4\leq J_{upp}\leq13$. In order to study the heating of the molecular gas, we present a combination of three diagnostic quantities to help determine the dominant heating source. Using the CO excitation, the CO J=1-0 linewidth, and the AGN contribution, we conclude that galaxies with large CO linewidths always have high-excitation CO ladders, and often low AGN contributions, suggesting that mechanical heating is important.

preprint2015arXiv

The identification of z-dropouts in Pan-STARRS1: three quasars at 6.5<z<6.7

Luminous distant quasars are unique probes of the high redshift intergalactic medium (IGM) and of the growth of massive galaxies and black holes in the early universe. Absorption due to neutral Hydrogen in the IGM makes quasars beyond a redshift of z~6.5 very faint in the optical $z$-band, thus locating quasars at higher redshifts require large surveys that are sensitive above 1 micron. We report the discovery of three new z>6.5 quasars, corresponding to an age of the universe of <850 Myr, selected as z-band dropouts in the Pan-STARRS1 survey. This increases the number of known z>6.5 quasars from 4 to 7. The quasars have redshifts of z=6.50, 6.52, and 6.66, and include the brightest z-dropout quasar reported to date, PSO J036.5078+03.0498 with M_1450=-27.4. We obtained near-infrared spectroscopy for the quasars and from the MgII line we estimate that the central black holes have masses between 5x10^8 and 4x10^9 M_sun, and are accreting close to the Eddington limit (L_Bol/L_Edd=0.13-1.2). We investigate the ionized regions around the quasars and find near zone radii of R_NZ=1.5-5.2 proper Mpc, confirming the trend of decreasing near zone sizes with increasing redshift found for quasars at 5.7<z<6.4. By combining R_NZ of the PS1 quasars with those of 5.7<z<7.1 quasars in the literature, we derive a luminosity corrected redshift evolution of R_NZ,corrected=(7.2+/-0.2)-(6.1+/-0.7)x(z-6) Mpc. However, the large spread in R_NZ in the new quasars implies a wide range in quasar ages and/or a large variation in the neutral Hydrogen fraction along different lines of sight.

preprint2015arXiv

The Impact of Molecular Gas on Mass Models of Nearby Galaxies

We present CO velocity fields and rotation curves for a sample of nearby galaxies, based on data from the HERACLES survey. We combine our data with literature THINGS, SINGS and KINGFISH results to provide a comprehensive sample of mass models of disk galaxies inclusive of molecular gas. We compare the kinematics of the molecular (CO from HERACLES) and atomic (${\rm H{\scriptstyle I}}$ from THINGS) gas distributions to determine the extent to which CO may be used to probe the dynamics in the inner part of galaxies. In general, we find good agreement between the CO and ${\rm H{\scriptstyle I}}$ kinematics with small differences in the inner part of some galaxies. We add the contribution of the molecular gas to the mass models in our galaxies by using two different conversion factors $\mathrm{α_{CO}}$ to convert CO luminosity to molecular gas mass surface density - the constant Milky Way value and the radially varying profiles determined in recent work based on THINGS, HERACLES and KINGFISH data. We study the relative effect that the addition of the molecular gas has upon the halo rotation curves for Navarro-Frenk-White (NFW) and the observationally motivated pseudo-isothermal halos. The contribution of the molecular gas varies for galaxies in our sample - for those galaxies where there is a substantial molecular gas content, using different values of $\mathrm{α_{CO}}$ can result in significant differences to the relative contribution of the molecular gas and, hence, the shape of the dark matter halo rotation curves in the central regions of galaxies.

preprint2015arXiv

THOR - The HI, OH, Recombination Line Survey of the Milky Way - The pilot study: HI observations of the giant molecular cloud W43

To study the atomic, molecular and ionized emission of Giant Molecular Clouds (GMCs), we have initiated a Large Program with the VLA: 'THOR - The HI, OH, Recombination Line survey of the Milky Way'. We map the 21cm HI line, 4 OH lines, 19 H_alpha recombination lines and the continuum from 1 to 2 GHz of a significant fraction of the Milky Way (l=15-67deg, |b|<1deg) at ~20" resolution. In this paper, we focus on the HI emission from the W43 star-formation complex. Classically, the HI 21cm line is treated as optically thin with column densities calculated under this assumption. This might give reasonable results for regions of low-mass star-formation, however, it is not sufficient to describe GMCs. We analyzed strong continuum sources to measure the optical depth, and thus correct the HI 21cm emission for optical depth effects and weak diffuse continuum emission. Hence, we are able to measure the HI mass of W43 more accurately and our analysis reveals a lower limit of M=6.6x10^6 M_sun, which is a factor of 2.4 larger than the mass estimated with the assumption of optically thin emission. The HI column densities are as high as N(HI)~150 M_sun/pc^2 ~ 1.9x10^22 cm^-2, which is an order of magnitude higher than for low mass star formation regions. This result challenges theoretical models that predict a threshold for the HI column density of ~10 M_sun/pc^2, at which the formation of molecular hydrogen should set in. By assuming an elliptical layered structure for W43, we estimate the particle density profiles. While at the cloud edge atomic and molecular hydrogen are well mixed, the center of the cloud is strongly dominated by molecular hydrogen. We do not identify a sharp transition between hydrogen in atomic and molecular form. Our results are an important characterization of the atomic to molecular hydrogen transition in an extreme environment and challenges current theoretical models.

preprint2014arXiv

[CII] 158 $μ$m Emission as a Star Formation Tracer

The [CII] 157.74 $μ$m transition is the dominant coolant of the neutral interstellar gas, and has great potential as a star formation rate (SFR) tracer. Using the Herschel KINGFISH sample of 46 nearby galaxies, we investigate the relation of [CII] surface brightness and luminosity with SFR. We conclude that [CII] can be used for measurements of SFR on both global and kiloparsec scales in normal star-forming galaxies in the absence of strong active galactic nuclei (AGN). The uncertainty of the $Σ_{\rm [CII]}-Σ_{\rm SFR}$ calibration is $\pm$0.21 dex. The main source of scatter in the correlation is associated with regions that exhibit warm IR colors, and we provide an adjustment based on IR color that reduces the scatter. We show that the color-adjusted $Σ_{\rm[CII]}-Σ_{\rm SFR}$ correlation is valid over almost 5 orders of magnitude in $Σ_{\rm SFR}$, holding for both normal star-forming galaxies and non-AGN luminous infrared galaxies. Using [CII] luminosity instead of surface brightness to estimate SFR suffers from worse systematics, frequently underpredicting SFR in luminous infrared galaxies even after IR color adjustment (although this depends on the SFR measure employed). We suspect that surface brightness relations are better behaved than the luminosity relations because the former are more closely related to the local far-UV field strength, most likely the main parameter controlling the efficiency of the conversion of far-UV radiation into gas heating. A simple model based on Starburst99 population-synthesis code to connect SFR to [CII] finds that heating efficiencies are $1\%-3\%$ in normal galaxies.

preprint2014arXiv

A far-IR view of the starburst driven superwind in NGC 2146

NGC 2146, a nearby luminous infrared galaxy (LIRG), presents evidence for outflows along the disk minor axis in all gas phases (ionized, neutral atomic and molecular). We present an analysis of the multi-phase starburst driven superwind in the central 5 kpc as traced in spatially resolved spectral line observations, using far-IR Herschel PACS spectroscopy, to probe the effects on the atomic and ionized gas, and optical integral field spectroscopy to examine the ionized gas through diagnostic line ratios. We observe an increased ~250 km/s velocity dispersion in the [OI] 63 micron, [OIII] 88 micron, [NII] 122 micron and [CII] 158 micron fine-structure lines that is spatially coincident with high excitation gas above and below the disk. We model this with a slow ~200 km/s shock and trace the superwind to the edge of our field of view 2.5 kpc above the disk. We present new SOFIA 37 micron observations to explore the warm dust distribution, and detect no clear dust entrainment in the outflow. The stellar kinematics appear decoupled from the regular disk rotation seen in all gas phases, consistent with a recent merger event disrupting the system. We consider the role of the superwind in the evolution of NGC 2146 and speculate on the evolutionary future of the system. Our observations of NGC 2146 in the far-IR allow an unobscured view of the wind, crucial for tracing the superwind to the launching region at the disk center, and provide a local analog for future ALMA observations of outflows in high redshift systems.

preprint2014arXiv

A low HI column density filament in NGC 2403: signature of interaction or accretion

Observed HI accretion around nearby galaxies can only account for a fraction of the gas supply needed to sustain the currently observed star formation rates. It is possible that additional accretion happens in the form of low column density cold flows, as predicted by numerical simulations of galaxy formation. To contrain the presence and properties of such flows, we present deep HI observations obtained with the NRAO Green Bank Telescope of an area measuring 4 by 4 degrees around NGC 2403. These observations, with a 5 sigma detection limit of 2.4 x 10^18 cm^-2 over a 20 km/s linewidth, reveal the presence of a low-column density, extended cloud outside the main HI disk, about 17' (~16 kpc or ~2R25) to the NW of the center of the galaxy. The total HI mass of the cloud is 6.3 x 10^6 Msun, or 0.15 percent of the total HI mass of NGC 2403. The cloud is associated with an 8-kpc anomalous-velocity HI filament in the inner disk, previously observed in deep VLA observations by Fraternali et al. (2001, 2002). We discuss several scenarios for the origin of the cloud, and conclude that it is either accreting from the intergalactic medium, or is the result of a minor interaction with a neigbouring dwarf galaxy.

preprint2014arXiv

ALLSMOG: an APEX Low-redshift Legacy Survey for MOlecular Gas. I - molecular gas scaling relations, and the effect of the CO/H2 conversion factor

We present ALLSMOG, the APEX Low-redshift Legacy Survey for MOlecular Gas. ALLSMOG is a survey designed to observe the CO(2-1) emission line with the APEX telescope, in a sample of local galaxies (0.01 < z < 0.03), with stellar masses in the range 8.5 < log(M*/Msun) < 10. This paper is a data release and initial analysis of the first two semesters of observations, consisting of 42 galaxies observed in CO(2-1). By combining these new CO(2-1) emission line data with archival HI data and SDSS optical spectroscopy, we compile a sample of low-mass galaxies with well defined molecular gas masses, atomic gas masses, and gas-phase metallicities. We explore scaling relations of gas fraction and gas consumption timescale, and test the extent to which our findings are dependent on a varying CO/H2 conversion factor. We find an increase in the H2/HI mass ratio with stellar mass which closely matches semi-analytic predictions. We find a mean molecular gas fraction for ALLSMOG galaxies of MH2/M* = (0.09 - 0.13), which decreases with stellar mass. We measure a mean molecular gas consumption timescale for ALLSMOG galaxies of 0.4 - 0.7 Gyr. We also confirm the non-universality of the molecular gas consumption timescale, which varies (with stellar mass) from ~100 Myr to ~2 Gyr. Importantly, we find that the trends in the H2/HI mass ratio, gas fraction, and the non-universal molecular gas consumption timescale are all robust to a range of recent metallicity-dependent CO/H2 conversion factors.

preprint2014arXiv

An ALMA Survey of Submillimeter Galaxies in the Extended Chandra Deep Field South: Near-infrared morphologies and stellar sizes

We analyse HST WFC3/$H_{160}$-band observations of a sample of 48 ALMA-detected submillimeter galaxies (SMGs) in the Extended Chandra Deep Field South field, to study their stellar morphologies and sizes. We detect 79$\pm$17% of the SMGs in the $H_{160}$-band imaging with a median sensitivity of 27.8 mag, and most (80%) of the non-detections are SMGs with 870$μ$m fluxes of $S_{870} < $3 mJy. With a surface brightness limit of $μ_H \sim $26 mag arcsec$^{-2}$, we find that 82$\pm$9% of the $H_{160}$-band detected SMGs at $z =$ 1-3 appear to have disturbed morphologies, meaning they are visually classified as either irregulars or interacting systems, or both. By determining a Sérsic fit to the $H_{160}$ surface-brightness profiles we derive a median Sérsic index of $n = $1.2$\pm$0.3 and a median half-light radius of $r_e = $4.4$^{+1.1}_{-0.5}$ kpc for our SMGs at $z = $1-3. We also find significant displacements between the positions of the $H_{160}$-component and 870$μ$m emission in these systems, suggesting that the dusty star-burst regions and less-obscured stellar distribution are not co-located. We find significant differences in the sizes and the Sérsic index between our $z = $2-3 SMGs and $z \sim $2 quiescent galaxies, suggesting a major transformation of the stellar light profile is needed in the quenching processes if SMGs are progenitors of the red-and-dead $z\sim$2 galaxies. Given the short-lived nature of SMGs, we postulate that the majority of the $z = $2-3 SMGs with $S_{870} \gtrsim $2 mJy are early/mid-stage major mergers.

preprint2014arXiv

An ALMA survey of submillimetre galaxies in the Extended Chandra Deep Field South: radio properties and the far-infrared/radio correlation

We present a study of the radio properties of 870$μ$m-selected submillimetre galaxies (SMGs), observed at high resolution with ALMA in the Extended Chandra Deep Field South. From our initial sample of 76 ALMA SMGs, we detect 52 SMGs at $>3σ$ significance in VLA 1400MHz imaging, of which 35 are also detected at $>3σ$ in new 610MHz GMRT imaging. Within this sample of radio-detected SMGs, we measure a median radio spectral index $α_{610}^{1400} = -0.79 \pm 0.06$, (with inter-quartile range $α=[-1.16,-0.56]$) and investigate the far-infrared/radio correlation via the parameter $q_{\rm IR}$, the logarithmic ratio of the rest-frame 8-1000$μ$m flux and monochromatic radio flux. Our median $q_{\rm IR} = 2.56 \pm 0.05$ (inter-quartile range $q_{\rm IR}=[2.42,2.78]$) is higher than that typically seen in single-dish 870$μ$m-selected sources ($q_{\rm IR} \sim 2.4$), which may reflect the fact that our ALMA-based study is not biased to radio-bright counterparts, as previous samples were. Finally, we search for evidence that $q_{\rm IR}$ and $α$ evolve with age in a co-dependent manner, as predicted by starburst models: the data populate the predicted region of parameter space, with the stellar mass tending to increase along tracks of $q_{\rm IR}$ versus $α$ in the direction expected, providing the first observational evidence in support of these models.

preprint2014arXiv

Black hole mass estimates and emission-line properties of a sample of redshift z>6.5 quasars

We present the analysis of optical and near-infrared spectra of the only four $z>6.5$ quasars known to date, discovered in the UKIDSS-LAS and VISTA-VIKING surveys. Our data-set consists of new VLT/X-Shooter and Magellan/FIRE observations. These are the best optical/NIR spectroscopic data that are likely to be obtained for the $z>6.5$ sample using current $6$ - $10$ m facilities. We estimate the black hole mass, the Eddington ratio, and the SiIV/CIV, CIII]/CIV, and FeII/MgII emission-line flux ratios. We perform spectral modeling using a procedure that allows us to derive a probability distribution for the continuum components and to obtain the quasar properties weighted upon the underlying distribution of continuum models. The $z>6.5$ quasars show the same emission properties as their counterparts at lower redshifts. The $z>6.5$ quasars host black holes with masses of $\sim 10^9$ M$_{\odot}$ that are accreting close to the Eddington luminosity ($\langle{\rm log} (L_{\rm Bol}/L_{\rm Edd})\rangle= -0.4\pm0.2$), in agreement with what has been observed for a sample of $4.0<z<6.5$ quasars. By comparing the SiIV/CIV and CIII]/CIV flux ratios with the results obtained from luminosity-matched samples at $z\sim6$ and $2\leq z\leq4.5$, we find no evidence of evolution of the line ratios with cosmic time. We compare the measured FeII/MgII flux ratios with those obtained for a sample of $4.0<z<6.4$ sources. The two samples are analyzed using a consistent procedure. There is no evidence that the FeII/MgII flux ratio evolves between $z=7$ and $z=4$. Under the assumption that the FeII/MgII traces the Fe/Mg abundance ratio, this implies the presence of major episodes of chemical enrichment in the quasar hosts in the first $\sim0.8$ Gyr after the Big Bang.

preprint2014arXiv

CO(1-0) line imaging of massive star-forming disc galaxies at z=1.5-2.2

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

preprint2014arXiv

COLDz: Karl G. Jansky Very Large Array discovery of a gas-rich galaxy in COSMOS

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

preprint2014arXiv

Combined CO & Dust Scaling Relations of Depletion Time and Molecular Gas Fractions with Cosmic Time, Specific Star Formation Rate and Stellar Mass

We combine molecular gas masses inferred from CO emission in 500 star forming galaxies (SFGs) between z=0 and 3, from the IRAM-COLDGASS, PHIBSS1/2 and other surveys, with gas masses derived from Herschel far-IR dust measurements in 512 galaxy stacks over the same stellar mass/redshift range. We constrain the scaling relations of molecular gas depletion time scale (tdepl) and gas to stellar mass ratio (Mmolgas/M*) of SFGs near the star formation main-sequence with redshift, specific star formation rate (sSFR) and stellar mass (M*). The CO- and dust-based scaling relations agree remarkably well. This suggests that the CO-H2 mass conversion factor varies little within 0.6dex of the main sequence (sSFR(ms,z,M*)), and less than 0.3dex throughout this redshift range. This study builds on and strengthens the results of earlier work. We find that tdepl scales as (1+z)^-0.3 *(sSFR/sSFR(ms,z,M*))^-0.5, with little dependence on M*. The resulting steep redshift dependence of Mmolgas/M* ~(1+z)^3 mirrors that of the sSFR and probably reflects the gas supply rate. The decreasing gas fractions at high M* are driven by the flattening of the SFR-M* relation. Throughout the redshift range probed a larger sSFR at constant M* is due to a combination of an increasing gas fraction and a decreasing depletion time scale. As a result galaxy integrated samples of the Mmolgas-SFR rate relation exhibit a super-linear slope, which increases with the range of sSFR. With these new relations it is now possible to determine Mmolgas with an accuracy of 0.1dex in relative terms, and 0.2dex including systematic uncertainties.

preprint2014arXiv

Constraining the Nature of Two Ly-alpha Emitters detected by ALMA at z = 4.7

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

preprint2014arXiv

Discovery of eight z ~ 6 quasars from Pan-STARRS1

High-redshift quasars are currently the only probes of the growth of supermassive black holes and potential tracers of structure evolution at early cosmic time. Here we present our candidate selection criteria from the Panoramic Survey Telescope & Rapid Response System 1 and follow-up strategy to discover quasars in the redshift range 5.7<z<6.2. With this strategy we discovered eight new 5.7<z<6.0 quasars, increasing the number of known quasars at z>5.7 by more than 10%. We additionally recovered 18 previously known quasars. The eight quasars presented here span a large range of luminosities (-27.3 < M_{1450} < -25.4; 19.6 < z_ps1 < 21.2) and are remarkably heterogeneous in their spectral features: half of them show bright emission lines whereas the other half show a weak or no Ly$α$ emission line (25% with rest-frame equivalent width of the Ly$α$ + Nv line lower than 15Å). We find a larger fraction of weak-line emission quasars than in lower redshift studies. This may imply that the weak-line quasar population at the highest redshifts could be more abundant than previously thought. However, larger samples of quasars are needed to increase the statistical significance of this finding.

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

Enabling the next generation of cm-wavelength studies of high-redshift molecular gas with the SKA

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

preprint2014arXiv

High-resolution [C II] imaging of HDF850.1 reveals a merging galaxy at z=5.185

New high-resolution maps with the IRAM Interferometer of the redshifted [C II] 158 micron line and the 0.98mm dust continuum of HDF850.1 at z = 5.185 show the source to have a blueshifted northern component and a redshifted southern component, with a projected separation of 0.3 arcsec, or 2 kpc. We interpret these components as primordial galaxies that are merging to form a larger galaxy. We think it is the resulting merger-driven starburst that makes HDF850.1 an ultraluminous infrared galaxy, with an L(IR) of 1E13 Lsun. The observed line and continuum brightness temperatures and the constant line-to-continuum ratio across the source imply (1) high [C II] line optical depth, (2) a [C II] excitation temperature of the same order as the dust temperature, and (3) dust continuum emission that is nearly optically thick at 158 microns. These conclusions for HDF850.1 probably also apply to other high-redshift submillimeter galaxies and quasar hosts in which the [C II] 158 micron line has been detected, as indicated by their roughly constant [C II]-to-158 micron continuum ratios, in sharp contrast to the large dispersion in their [C II]-to-FIR luminosity ratios. In brightness temperature units, the [C II] line luminosity is about the same as the predicted CO(1-0) luminosity, implying that the [C II] line can also be used to estimate the molecular gas mass, with the same assumptions as for CO.

preprint2014arXiv

Karl G. Jansky Very Large Array observations of cold dust and molecular gas in starbursting quasar host galaxies at z~4.5

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

preprint2014arXiv

Measuring Quasar Variability with Pan-STARRS1 and SDSS

We measure quasar variability using the Panoramic Survey Telescope and Rapid Response System 1 Survey (Pan-STARRS1 or PS1) and the Sloan Digital Sky Survey (SDSS) and establish a method of selecting quasars via their variability in 10,000 square degree surveys. We use 100,000 spectroscopically confirmed quasars that have been well measured in both PS1 and SDSS and take advantage of the decadal time scales that separate SDSS measurements and PS1 measurements. A power law model fits the data well over the entire time range tested, 0.01 to 10 years. Variability in the current PS1-SDSS dataset can efficiently distinguish between quasars and non-varying objects. It improves the purity of a griz quasar color cut from 4.1% to 48% while maintaining 67% completeness. Variability will be very effective at finding quasars in datasets with no u band and in redshift ranges where exclusively photometric selection is not efficient. We show that quasars' rest-frame ensemble variability, measured as a root mean squared in delta magnitudes, is consistent with V(z, L, t) = A0 (1+z)^0.37 (L/L0)^-0.16 (t/1yr)^0.246 , where L0 = 10^46 ergs^-1 and A0 = 0.190, 0.162, 0.147 or 0.141 in the gP1 , rP1 , iP1 or zP1 filter, respectively. We also fit across all four filters and obtain median variability as a function of z, L and lambda as V(z, L, lambda, t) = 0.079(1 + z)^0.15 (L/L0 )^-0.2 (lambda/1000 nm)^-0.44 (t/1yr)^0.246 .

preprint2014arXiv

Spectral Energy Distributions of QSOs at z>5: common AGN-heated dust and occasionally strong star-formation

We present spectral energy distributions (SEDs) of 69 QSOs at z>5, covering a rest frame wavelength range of 0.1mu to ~80mu, and centered on new Spitzer and Herschel observations. The detection rate of the QSOs with Spitzer is very high (97% at lambda_rest ~< 4mu), but drops towards the Herschel bands with 30% detected in PACS (rest frame mid-infrared) and 15% additionally in the SPIRE (rest frame far-infrared; FIR). We perform multi-component SED fits for Herschel-detected objects and confirm that to match the observed SEDs, a clumpy torus model needs to be complemented by a hot (~1300K) component and, in cases with prominent FIR emission, also by a cold (~50K) component. In the FIR detected cases the luminosity of the cold component is on the order of 10^13 L_sun which is likely heated by star formation. From the SED fits we also determine that the AGN dust-to-accretion disk luminosity ratio declines with UV/optical luminosity. Emission from hot (~1300K) dust is common in our sample, showing that nuclear dust is ubiquitous in luminous QSOs out to redshift 6. However, about 15% of the objects appear under-luminous in the near infrared compared to their optical emission and seem to be deficient in (but not devoid of) hot dust. Within our full sample, the QSOs detected with Herschel are found at the high luminosity end in L_UV/opt and L_NIR and show low equivalent widths (EWs) in H_alpha and in Ly_alpha. In the distribution of H_alpha EWs, as determined from the Spitzer photometry, the high-redshift QSOs show little difference to low redshift AGN.

preprint2014arXiv

Star formation relations and CO SLEDs across the J-ladder and redshift

We present FIR-CO luminosity relations ($\log L_{\rm FIR} = α\log L'_{\rm CO} + β$) for the full CO rotational ladder from J=1-0 to J=13-12 for 62 local (z < 0.1) (Ultra) Luminous Infrared Galaxies (LIRGs) using data from Herschel SPIRE-FTS and ground-based telescopes. We extend our sample to high redshifts (z > 1) by including 35 (sub)-millimeter selected dusty star forming galaxies from the literature with robust CO observations. The addition of luminous starbursts at high redshifts enlarge the range of the FIR-CO luminosity relations towards the high-IR-luminosity end while also significantly increasing the small amount of mid-/high-J CO line data available prior to Herschel. This new data-set (both in terms of IR luminosity and J-ladder) reveals linear FIR-CO luminosity relations ($α\sim 1$) for J=1-0 up to J=5-4, with a nearly constant normalisation ($β\sim 2$). This is expected from the (also) linear FIR-(molecular line) relations found for the dense gas tracer lines (HCN and CS), as long as the dense gas mass fraction does not vary strongly within our (merger/starburst)-dominated sample. However from J=6-5 and up to J=13-12 we find an increasingly sub-linear slope and higher normalization constant with increasing J. We argue that these are caused by a warm (~100K) and dense ($>10^4{\rm cm^{-3}}$) gas component whose thermal state is unlikely to be maintained by star formation powered far-UV radiation fields (and thus is no longer directly tied to the star formation rate). We suggest that mechanical heating (e.g., supernova driven turbulence and shocks), and not cosmic rays, is the more likely source of energy for this component. The global CO spectral line energy distributions (SLEDs), which remain highly excited from J=6-5 up to J=13-12, are found to be a generic feature of the (U)LIRGs in our sample, and further support the presence of this gas component.

preprint2014arXiv

The impact of the gas distribution on the determination of dynamical masses of galaxies using unresolved observations

The dynamical mass (M_dyn) is a key property of any galaxy, yet a determination of M_dyn is not straight-forward if spatially resolved measurements are not available. This situation occurs in single-dish HI observations of the local universe, but also frequently in high-redshift observations. M_dyn-measurements in high-redshift galaxies are commonly obtained through observations of the CO line, the most abundant tracer of the molecular medium. Even though the CO linewidth can in most cases be determined with reasonable accuracy, a measurement of the size of the emitting region is typically challenging given current facilities. We show how the integrated spectra (`global profiles') of a variety of galaxy models depend on the spatial distribution of the tracer gas as well as its velocity dispersion. We demonstrate that the choice of tracer emission line significantly affects the shape of the global profiles. In particular, in the case of high (~50 kms-1) velocity dispersions, compact tracers (such as CO) result in Gaussian-like (non-double-horned) profiles, as is indeed frequently seen in high-redshift observations. We determine at which radii the rotation curve reaches the rotation velocity corresponding to the velocity width, and find that for each tracer this happens at a well-defined radius: HI velocity widths typically originate at ~5 optical scale lengths, while CO velocity widths trace the rotation velocity at ~2 scale lengths. We additionally explore other distributions to take into account that CO distributions at high redshift likely differ from those at low redshift. Our models, while not trying to reproduce individual galaxies, define characteristic radii that can be used in conjunction with the measured velocity widths in order to define dynamical masses consistent with the assumed gas distribution.

preprint2014arXiv

The kiloparsec-scale star formation law at redshift 4: wide-spread, highly efficient star formation in the dust-obscured starburst galaxy GN20

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

preprint2014arXiv

The Survey of Lines in M31 (SLIM): Investigating the Origins of [CII] Emission

The [CII] 158 micron line is one of the strongest emission lines observed in star-forming galaxies, and has been empirically measured to correlate with the star formation rate (SFR) globally and on ~kpc scales. However, due to the multi-phase origins of [CII], one might expect this relation to break down at small scales. We investigate the origins of [CII] emission by examining high spatial resolution observations of [CII] in M31, with the Survey of Lines in M31 (SLIM). We present five ~700x700 pc (3"x3") Fields mapping the [CII] emission, Halpha emission, combined with ancillary infrared (IR) data. We spatially separate star-forming regions from diffuse gas and dust emission on ~50 pc scales. We find that the [CII] - SFR correlation holds even at these scales, although the relation typically has a flatter slope than found at larger (~kpc) scales. While the Halpha emission in M31 is concentrated in the SFR regions, we find that a significant amount (~20-90%) of the [CII] emission comes from outside star-forming regions, and that the total IR (TIR) emission has the highest diffuse fraction of all SFR tracers. We find a weak correlation of the [CII]/TIR to dust color in each Field, and find a large scale trend of increasing [CII]/TIR with galactocentric radius. The differences in the relative diffuse fractions of [CII], Halpha and IR tracers are likely caused by a combination of energetic photon leakage from HII regions and heating by the diffuse radiation field arising from older (B-star) stellar populations. However, we find that by averaging our measurements over ~kpc scales, these effects are minimized, and the relation between [CII] and SFR found in other nearby galaxy studies is retrieved.

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 sub-millimeter galaxies in the Extended Chandra Deep Field South: Sub-millimeter properties of color-selected galaxies

We study the sub-mm properties of color-selected galaxies via a stacking analysis applied for the first time to interferometric data at sub-mm wavelengths. We base our study on 344 GHz ALMA continuum observations of ~20''-wide fields centered on 86 sub-mm sources detected in the LABOCA Extended Chandra Deep Field South Sub-mm Survey (LESS). We select various classes of galaxies (K-selected, star-forming sBzK galaxies, extremely red objects and distant red galaxies) according to their optical/NIR fluxes. We find clear, >10-sigma detections in the stacked images of all these galaxy classes. We include in our stacking analysis Herschel/SPIRE data to constrain the dust SED of these galaxies. We find that their dust emission is well described by a modified black body with Tdust~30 K and beta=1.6 and IR luminosities of (5-11)x10^{11} Lsun, or implied star formation rates of 75-140 Msun/yr. We compare our results with those of previous studies based on single-dish observations at 870 micron and find that our flux densities are a factor 2-3 higher than previous estimates. The discrepancy is observed also after removing sources individually detected in ALESS maps. We report a similar discrepancy by repeating our analysis on 1.4\,GHz observations of the whole ECDFS. Hence we find tentative evidence that galaxies that are associated in projected and redshift space with sub-mm bright sources are brighter than the average population. Finally, we put our findings in the context of the cosmic star formation rate density as a function of redshift.

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

An ALMA Survey of Submillimeter Galaxies in the Extended Chandra Deep Field South: The AGN Fraction and X-ray Properties of Submillimeter Galaxies

The large gas and dust reservoirs of submm galaxies (SMGs) could potentially provide ample fuel to trigger an Active Galactic Nucleus (AGN), but previous studies of the AGN fraction in SMGs have been controversial largely due to the inhomogeneity and limited angular resolution of the available submillimeter surveys. Here we set improved constraints on the AGN fraction and X-ray properties of the SMGs with ALMA and Chandra observations in the Extended Chandra Deep Field-South (E-CDF-S). This study is the first among similar works to have unambiguously identified the X-ray counterparts of SMGs; this is accomplished using the fully submm-identified, statistically reliable SMG catalog with 99 SMGs from the ALMA LABOCA E-CDF-S Submillimeter Survey (ALESS). We found 10 X-ray sources associated with SMGs (median redshift z = 2.3), of which 8 were identified as AGNs using several techniques that enable cross-checking. The other 2 X-ray detected SMGs have levels of X-ray emission that can be plausibly explained by their star-formation activity. 6 of the 8 SMG-AGNs are moderately/highly absorbed, with N_H > 10e23 cm-2. An analysis of the AGN fraction, taking into account the spatial variation of X-ray sensitivity, yields an AGN fraction of 17+16-6% for AGNs with rest-frame 0.5-8 keV absorption-corrected luminosity >7.8x10e42 erg s-1; we provide estimated AGN fractions as a function of X-ray flux and luminosity. ALMA's high angular resolution also enables direct X-ray stacking at the precise positions of SMGs for the first time, and we found 4 potential SMG-AGNs in our stacking sample.

preprint2013arXiv

An X-rays Survey of the Young Stellar Population of the Lynds 1641 and Iota Orionis Regions

We present an XMM-Newton survey of the part of Orion A cloud south of the Orion Nebula. This survey includes the Lynds 1641 (L1641) dark cloud, a region of the Orion A cloud with very few massive stars and hence a relatively low ambient UV flux, and the region around the O9 III star Iota Orionis. In addition to proprietary data, we used archival XMM data of the Orion Nebula Cluster (ONC) to extend our analysis to a major fraction of the Orion A cloud. We have detected 1060 X-ray sources in L1641 and Iota Ori region. About 94% of the sources have 2MASS & Spitzer counterparts, 204 and 23 being Class II and Class I or protostars objects, respectively. In addition, we have identified 489 X-ray sources as counterparts to Class III candidates, given they are bright in X-rays and appear as normal photospheres at mid-IR wavelengths. The remaining 205 X-ray sources are likely distant AGNs or other galactic sources not related to Orion A. We find that Class III candidates appear more concentrated in two main clusters in L1641. The first cluster of Class III stars is found toward the northern part of L1641, concentrated around Iota Ori. The stars in this cluster are more evolved than those in the Orion Nebula. We estimate a distance of 300-320 pc for this cluster and thus it is closer than the Orion A cloud. Another cluster rich in Class III stars is located in L1641 South and appears to be a slightly older cluster embedded in the Orion A cloud. Furthermore, other evolved Class III stars are found north of the ONC toward NGC 1977.

preprint2013arXiv

Complete infrared spectral energy distributions of mm detected quasars at z>5

We present Herschel far-infrared (FIR) photometry of eleven quasars at redshift z>5 that have previously been detected at 1.2mm. We perform full spectral energy distribution (SED) fits over the wavelength range lambda_rest ~0.1-400mu for those objects with good Herschel detections. These fits reveal the need for an additional far-infrared (FIR) component besides the emission from a dusty AGN-powered torus. This additional FIR component has temperatures of T_FIR ~ 40-60K with luminosities of L_(8-1000mu) ~ 10^13 L_sun (accounting for 25-60% of the bolometric FIR luminosity). If the FIR dust emission is due to star formation it would suggest star formation rates in excess of 1000 solar masses per year. We show that at long wavelengths (lambda_rest > 50mu) the contribution of the AGN-powered torus emission is negligible. This explains how previous FIR studies of high-redshift quasars that relied on single component fits to (ground-based) observations at lambda_obs > 350mu reached T_FIR and L_FIR values similar to our complete SED fits. Stacking the Herschel data of four individually undetected sources reveals a significant average signal in the PACS bands but not in SPIRE. The average SED of sources with individual Herschel detections shows a striking surplus in near- and mid-infrared emission when compared to common AGN templates. The comparison between two average SEDs (sources with and without individual Herschel detections) matched in the UV/optical indicates that for these objects the strength of the MIR emission may correlate with the strength of the FIR emission.

preprint2013arXiv

High-resolution spectroscopic imaging of CO in a z=4.05 proto-cluster

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

preprint2013arXiv

Physical Conditions of the Gas in an ALMA [CII]-identified Submillimetre Galaxy at z = 4.44

We present CO(2-1) observations of the submillimetre galaxy ALESS65.1 performed with the Australia Telescope Compact Array at 42.3 GHz. A previous ALMA study of submillimetre galaxies in the Extended Chandra Deep Field South detected [CII] 157.74 micron emission from this galaxy at a redshift of z = 4.44. No CO(2-1) emission was detected but we derive a firm upper limit to the cold gas mass in ALESS65.1 of M_gas < 1.7 x 10^10 M_odot. The estimated gas depletion timescale is <50 Myr, which is similar to other high redshift SMGs, and consistent with z > 4 SMGs being the likely progenitors of massive red-and-dead galaxies at z > 2. The ratio of the [CII], CO and far-infrared luminosities implies a strong far-ultraviolet field of G_0 > 10^3, as seen in Galactic star forming regions or local ULIRGs. The observed L_[CII]/L_FIR = 2.3 x 10^{-3} is high compared to local ULIRGs and, combined with L_[CII]/L_CO > 2700, it is consistent with ALESS65.1 either having an extended (several kpc) [CII] emitting region or lower than solar metallicity.

preprint2013arXiv

Shock Excited Molecules in NGC 1266: ULIRG conditions at the center of a Bulge Dominated Galaxy

We investigate the far infrared spectrum of NGC 1266, a S0 galaxy that contains a massive reservoir of highly excited molecular gas. Using the SPIRE-FTS, we detect the $^{12}$CO ladder up to J=(13-12), [C I] and [N II] lines, and also strong water lines more characteristic of UltraLuminous IR Galaxies (ULIRGs). The 12CO line emission is modeled with a combination of a low-velocity C-shock and a PDR. Shocks are required to produce the H2O and most of the high-J 12CO emission. Despite having an infrared luminosity thirty times less than a typical ULIRG, the spectral characteristics and physical conditions of the ISM of NGC 1266 closely resemble those of ULIRGs, which often harbor strong shocks and large-scale outflows.

preprint2013arXiv

The Highest Redshift Quasar at $z=7.085$: A Radio Quiet Source

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

preprint2013arXiv

Towards A Removal of Temperature Dependencies from Abundance Determinations

The metal content of a galaxy is a key property for distinguishing between viable galaxy evolutionary scenarios, and it strongly influences many of the physical processes in the interstellar medium. An absolute and robust determination of extragalactic metallicities is essential in constraining models of chemical enrichment and chemical evolution, however, current gas phase abundance determinations from optical fine-structure lines are uncertain to 0.8 dex as conversion of these optical line fluxes to abundances is strongly dependent on the electron temperature of the ionized gas. In contrast, the far-IR emission lines can be used to derive an O++ abundance that is relatively insensitive to temperature, while the ratio of the optical to far-IR lines provides a consistent temperature to be used in the derivation of an O$^+$ abundance. We present observations of the [O III] 88 \micron\ fine-structure line in NGC 628 that were obtained as part of the KINGFISH program. These data are combined with optical IFU data to derive oxygen abundances for seven HII regions. We find the abundance of these regions to all lie between the high and low values of strong line calibrations and in agreement with estimates that assume temperature fluctuations are present in the HII regions.

preprint2012arXiv

[CII] line emission in massive star-forming galaxies at z=4.7

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of the [CII] 157.7micron fine structure line and thermal dust continuum emission from a pair of gas-rich galaxies at z=4.7, BR1202-0725. This system consists of a luminous quasar host galaxy and a bright submm galaxy (SMG), while a fainter star-forming galaxy is also spatially coincident within a 4" (25 kpc) region. All three galaxies are detected in the submm continuum, indicating FIR luminosities in excess of 10^13 Lsun for the two most luminous objects. The SMG and the quasar host galaxy are both detected in [CII] line emission with luminosities, L([CII]) = (10.0 +/- 1.5)x10^9 Lsun and L([CII]) = (6.5+/-1.0)x10^9 Lsun, respectively. We estimate a luminosity ratio, L([CII])/L(FIR) = (8.3+/-1.2)x10^-4 for the starburst SMG to the North, and L([CII])/L(FIR) = (2.5+/-0.4)x10^-4 for the quasar host galaxy, in agreement with previous high-redshift studies that suggest lower [CII]-to-FIR luminosity ratios in quasars than in starburst galaxies. The third fainter object with a flux density, S(340GHz) = 1.9+/-0.3 mJy, is coincident with a Ly-Alpha emitter and is detected in HST ACS F775W and F814W images but has no clear counterpart in the H-band. Even if this third companion does not lie at a similar redshift to BR1202-0725, the quasar and the SMG represent an overdensity of massive, infrared luminous star-forming galaxies within 1.3 Gyr of the Big Bang.

preprint2012arXiv

A Study of Heating and Cooling of the ISM in NGC 1097 with Herschel-PACS and Spitzer-IRS

NGC 1097 is a nearby Seyfert 1 galaxy with a bright circumnuclear starburst ring, a strong large-scale bar and an active nucleus. We present a detailed study of the spatial variation of the far infrared (FIR) [CII]158um and [OI]63um lines and mid-infrared H2 emission lines as tracers of gas cooling, and of the polycyclic aromatic hydrocarbon (PAH) bands as tracers of the photoelectric heating, using Herschel-PACS, and Spitzer-IRS infrared spectral maps. We focus on the nucleus and the ring, and two star forming regions (Enuc N and Enuc S). We estimated a photoelectric gas heating efficiency ([CII]158um+[OI]63um)/PAH in the ring about 50% lower than in Enuc N and S. The average 11.3/7.7um PAH ratio is also lower in the ring, which may suggest a larger fraction of ionized PAHs, but no clear correlation with [CII]158μm/PAH(5.5 - 14um) is found. PAHs in the ring are responsible for a factor of two more [CII]158um and [OI]63um emission per unit mass than PAHs in the Enuc S. SED modeling indicates that at most 25% of the FIR power in the ring and Enuc S can come from high intensity photodissociation regions (PDRs), in which case G0 ~ 10^2.3 and nH ~ 10^3.5 cm^-3 in the ring. For these values of G0 and nH PDR models cannot reproduce the observed H2 emission. Much of the the H2 emission in the starburst ring could come from warm regions in the diffuse ISM that are heated by turbulent dissipation or shocks.

preprint2012arXiv

Cool Dust in the Outer Ring of NGC 1291

We examine Herschel Space Observatory images of one nearby prototypical outer ring galaxy, NGC 1291, and show that the ring becomes more prominent at wavelengths longer than 160um. The mass of cool dust in the ring dominates the total dust mass of the galaxy, accounting for at least 70% of it. The temperature of the emitting dust in the ring (T=19.5+/-0.3K) is cooler than that of the inner galaxy (T=25.7+/-0.7K). We discuss several explanations for the difference in dust temperature, including age and density differences in the stellar populations of the ring versus the bulge.

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

Detection of atomic carbon [CII] 158 micron and dust emission from a z=7.1 quasar host galaxy

Using the IRAM Plateau de Bure Interferometer, we report the detection of the 158 micron [CII] emission line and underlying dust continuum in the host galaxy of the quasar ULAS J112001.48+064124.3 (hereafter J1120+0641) at z=7.0842+/-0.0004. This is the highest redshift detection of the [CII] line to date, and allows us to put first constraints on the physical properties of the host galaxy. The [CII] line luminosity is (1.2+/-0.2)x10^9 Lsun, which is a factor ~4 lower than observed in a luminous quasar at z=6.42 (SDSS J1148+5251). The underlying far-infrared (FIR) continuum has a flux density of 0.61+/-0.16 mJy, similar to the average flux density of z~6 quasars that were not individually detected in the rest-frame FIR. Assuming the FIR luminosity of L_FIR = 5.8x10^11-1.8x10^12 Lsun is mainly powered by star-formation, we derive a star-formation rate in the range 160-440 Msun/yr and a total dust mass in the host galaxy of 6.7x10^7-5.7x10^8 Msun (both numbers have significant uncertainties given the unknown nature of dust at these redshifts). The [CII] line width of sigma_V=100+/-15 km/s is among the smallest observed when compared to the molecular line widths detected in z~6 quasars. Both the [CII] and dust continuum emission are spatially unresolved at the current angular resolution of 2.0x1.7 arcsec^2 (corresponding to 10x9 kpc^2 at the redshift of J1120+0641).

preprint2012arXiv

Evidence for a clumpy, rotating gas disk in a submillimeter galaxy at z=4

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

preprint2012arXiv

Evidence for CO shock excitation in NGC 6240 from Herschel SPIRE spectroscopy

We present Herschel SPIRE FTS spectroscopy of the nearby luminous infrared galaxy NGC 6240. In total 20 lines are detected, including CO J=4-3 through J=13-12, 6 H2O rotational lines, and [CI] and [NII] fine-structure lines. The CO to continuum luminosity ratio is 10 times higher in NGC 6240 than Mrk 231. Although the CO ladders of NGC 6240 and Mrk 231 are very similar, UV and/or X-ray irradiation are unlikely to be responsible for the excitation of the gas in NGC 6240. We applied both C and J shock models to the H2 v=1-0 S(1) and v=2-1 S(1) lines and the CO rotational ladder. The CO ladder is best reproduced by a model with shock velocity v_s=10 km s^-1 and a pre-shock density n_H=5 * 10^4 cm^-3. We find that the solution best fitting the H2 lines is degenerate: The shock velocities and number densities range between v_s = 17 - 47 km s^-1 and n_H=10^7 - 5 * 10^4 cm^-3, respectively. The H2 lines thus need a much more powerful shock than the CO lines. We deduce that most of the gas is currently moderately stirred up by slow (10 km s^-1) shocks while only a small fraction (< 1 percent) of the ISM is exposed to the high velocity shocks. This implies that the gas is rapidly loosing its highly turbulent motions. We argue that a high CO line-to-continuum ratio is a key diagnostic for the presence of shocks.

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

Evidence of strong quasar feedback in the early Universe

Most theoretical models invoke quasar driven outflows to quench star formation in massive galaxies, this feedback mechanism is required to account for the population of old and passive galaxies observed in the local universe. The discovery of massive, old and passive galaxies at z=2, implies that such quasar feedback onto the host galaxy must have been at work very early on, close to the reionization epoch. We have observed the [CII]158um transition in SDSSJ114816.64+525150.3 that, at z=6.4189, is one of the most distant quasars known. We detect broad wings of the line tracing a quasar-driven massive outflow. This is the most distant massive outflow ever detected and is likely tracing the long sought quasar feedback, already at work in the early Universe. The outflow is marginally resolved on scales of about 16 kpc, implying that the outflow can really affect the whole galaxy, as required by quasar feedback models. The inferred outflow rate, dM/dt > 3500 Msun/yr, is the highest ever found. At this rate the outflow can clean the gas in the host galaxy, and therefore quench star formation, in a few million years.

preprint2012arXiv

Gas fraction and star formation efficiency at z < 1.0

After new observations of 39 galaxies at z = 0.6-1.0 obtained at the IRAM 30m telescope, we present our full CO line survey covering the redshift range 0.2 < z < 1. Our aim is to determine the driving factors accounting for the steep decline in the star formation rate during this epoch. We study both the gas fraction, defined as Mgas/(Mgas+Mstar), and the star formation efficiency (SFE) defined by the ratio between far-infrared luminosity and molecular gas mass (LFIR/M(H2), i.e. a measure for the inverse of the gas depletion time. The sources are selected to be ultra-luminous infrared galaxies (ULIRGs), with LFIR greater than 10^12 Lo and experiencing starbursts. When we adopt a standard ULIRG CO-to-H2 conversion factor, their molecular gas depletion time is less than 100 Myr. Our full survey has now filled the gap of CO observations in the 0.2<z<1 range covering almost half of cosmic history. The detection rate in the 0.6 < z < 1 interval is 38% (15 galaxies out of 39), compared to 60% for the 0.2<z<0.6 interval. The average CO luminosity is L'CO = 1.8 10^10 K km/s pc^2, corresponding to an average H2 mass of 1.45 10^10 Mo. From observation of 7 galaxies in both CO(2-1) and CO(4-3), a high gas excitation has been derived; together with the dust mass estimation, this supports the choice of our low ULIRG conversion factor between CO luminosity and H2, for our sample sources. We find that both the gas fraction and the SFE significantly increase with redshift, by factors of 3 +-1 from z=0 to 1, and therefore both quantities play an important role and complement each other in cosmic star formation evolution.

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

Modeling Dust and Starlight in Galaxies Observed by Spitzer and Herschel: NGC 628 and NGC 6946

We characterize the dust in NGC628 and NGC6946, two nearby spiral galaxies in the KINGFISH sample. With data from 3.6um to 500um, dust models are strongly constrained. Using the Draine & Li (2007) dust model, (amorphous silicate and carbonaceous grains), for each pixel in each galaxy we estimate (1) dust mass surface density, (2) dust mass fraction contributed by polycyclic aromatic hydrocarbons (PAH)s, (3) distribution of starlight intensities heating the dust, (4) total infrared (IR) luminosity emitted by the dust, and (5) IR luminosity originating in regions with high starlight intensity. We obtain maps for the dust properties, which trace the spiral structure of the galaxies. The dust models successfully reproduce the observed global and resolved spectral energy distributions (SEDs). The overall dust/H mass ratio is estimated to be 0.0082+/-0.0017 for NGC628, and 0.0063+/-0.0009 for NGC6946, consistent with what is expected for galaxies of near-solar metallicity. Our derived dust masses are larger (by up to a factor 3) than estimates based on single-temperature modified blackbody fits. We show that the SED fits are significantly improved if the starlight intensity distribution includes a (single intensity) "delta function" component. We find no evidence for significant masses of cold dust T<12K. Discrepancies between PACS and MIPS photometry in both low and high surface brightness areas result in large uncertainties when the modeling is done at PACS resolutions, in which case SPIRE, MIPS70 and MIPS160 data cannot be used. We recommend against attempting to model dust at the angular resolution of PACS.

preprint2012arXiv

On the Variations of Fundamental Constants and AGN feedback in the QSO host galaxy RXJ0911.4+0551 at z=2.79

We report on sensitive observations of the CO(7-6) and CI(2-1) transitions in the z=2.79 QSO host galaxy RXJ0911.4+0551 using the IRAM Plateau de Bure interferometer (PdBI). Our extremely high signal to noise spectra combined with the narrow CO line width of this source (FWHM = 120 km/s) allows us to estimate sensitive limits on the space-time variations of the fundamental constants using two emission lines. Our observations show that the CI and CO line shapes are in good agreement with each other but that the CI line profile is of order 10% narrower, presumably due to the lower opacity in the latter line. Both lines show faint wings with velocities up to +/-250 km/s, indicative of a molecular outflow. As such the data provide direct evidence for negative feedback in the molecular gas phase at high redshift. Our observations allow us to determine the observed frequencies of both transitions with so far unmatched accuracy at high redshift. The redshift difference between the CO and CI lines is sensitive to variations of dF/F with F=alpha^2/mu where alpha is the fine structure constant and mu the proton-to-electron mass ratio. We find dF/F=6.9 +/-3.7 x 10^-6 at a lookback time of 11.3 Gyr, which within the uncertainties, is consistent with no variations of the fundamental constants.

preprint2012arXiv

Resolved [CII] emission in a lensed quasar at z=4.4

We present one of the first resolved maps of the [CII] 158 micron line, a powerful tracer of the star forming inter-stellar medium, at high redshift. We use the new IRAM PdBI receivers at 350 GHz to map this line in BRI 0952-0115, the host galaxy of a lensed quasar at z=4.4 previously found to be very bright in [CII] emission. The [CII] emission is clearly resolved and our data allow us to resolve two [CII] lensed images associated with the optical quasar images. We find that the star formation, as traced by [CII], is distributed over a region of ~ 1 kpc in size near the quasar nucleus, and we infer a star formation surface density >150 Msun/yr/kpc^2, similar to that observed in local ULIRGs. We also reveal another [CII] component, extended over ~ 12 kpc, and located at ~ 10 kpc from the quasar. We suggest that this component is a companion disk galaxy, in the process of merging with the quasar host, whose rotation field is distorted by the interaction with the quasar host, and where star formation, although intense, is more diffuse. These observations suggest that galaxy merging at high-z can enhance star formation at the same time in the form of more compact regions, in the vicinity of the accreting black hole, and in more extended star forming galaxies.

preprint2012arXiv

Resolving the Far-IR Line Deficit: Photoelectric Heating and Far-IR Line Cooling in NGC 1097 and NGC 4559

The physical state of interstellar gas and dust is dependent on the processes which heat and cool this medium. To probe heating and cooling of the ISM over a large range of infrared surface brightness, on sub-kiloparsec scales, we employ line maps of [C \ii] 158 $μ$m, [O \one] 63 $μ$m, and [N \ii] 122 $μ$m in NGC 1097 and NGC 4559, obtained with the PACS spectrometer onboard {\it Herschel}. We matched new observations to existing Spitzer-IRS data that trace the total emission of polycyclic aromatic hydrocarbons (PAHs). We confirm at small scales in these galaxies that the canonical measure of photoelectric heating efficiency, ([C \ii] + [O \one])/TIR, decreases as the far-infrared color, $νf_ν$(70 $μ$m)/$νf_ν$(100 $μ$m), increases. In contrast, the ratio of far-infrared (far-IR) cooling to total PAH emission, ([C \ii] + [O \one])/PAH, is a near constant $\sim$6% over a wide range of far-infrared color, 0.5 \textless\ $νf_ν$(70 $μ$m)/$νf_ν$(100 $μ$m) $\lesssim$ 0.95. In the warmest regions, where $νf_ν$(70 $μ$m)/$νf_ν$(100 $μ$m) $\gtrsim$ 0.95, the ratio ([C \ii] + [O \one])/PAH drops rapidly to 4%. We derived representative values of the local UV radiation density, $G_0$, and the gas density, $n_H$, by comparing our observations to models of photodissociation regions. The ratio $G_0/n_H$, derived from fine-structure lines, is found to correlate with the mean dust-weighted starlight intensity, $<U>$ derived from models of the IR SED. Emission from regions that exhibit a line deficit is characterized by an intense radiation field, indicating that small grains are susceptible to ionization effects. We note that there is a shift in the 7.7 / 11.3 $μ$m PAH ratio in regions that exhibit a deficit in ([C \ii] + [O \one])/PAH, suggesting that small grains are ionized in these environments.

preprint2012arXiv

Spatially resolved kinematics of the multi-phase interstellar medium in the inner disk of M82

We present spatially resolved kinematics of the interstellar NaI D 5890,5896 doublet absorption and 12CO(1-0) emission across the inner ~2x1 kpc of the disk of M82. These data were obtained with the DensePak IFU on the WIYN telescope and the Caltech Owens Valley Radio Observatory (OVRO) millimetre array. By measuring the NaI and CO (and Halpha kinematics from a previous study) at the same spatial resolution, and employing the same line fitting method, we have been able to make meaningful comparisons between the ionized, neutral and molecular gas phases. We detect a component of the NaI line throughout the inner disk with velocities that are forbidden by the known galactic rotation. We interpret this as originating in counter-rotating or perhaps inflowing material. In the southern plume, we find clear evidence of entrained CO gas with kinematics consistent with that of Halpha. On the northern side, the CO kinematics appear to trace more static clouds in the inner halo that could be pre-existing or tidal in origin. We find no evidence that NaI absorption is kinematically associated with the outflow. We conclude that a combination of lack of velocity resolution and confusion of due to the high inclination of the system is acting to prevent detection. Thus, in the search for neutral outflows from galaxies, the signature high velocity components may easily be missed in observations at low spectral resolution and/or sensitivity, and particularly so in highly inclined systems.

preprint2012arXiv

Star Formation Efficiency at Intermediate Redshift

Star formation is evolving very fast in the second half of the Universe, and it is yet unclear whether this is due to evolving gas content, or evolving star formation efficiency (SFE). We have carried out a survey of ultra-luminous galaxies (ULIRG) between z=0.2 and 1, to check the gas fraction in this domain of redshift which is still poorly known. Our survey with the IRAM-30m detected 33 galaxies out of 69, and we derive a significant evolution of both the gas fraction and SFE of ULIRGs over the whole period, and in particular a turning point around z=0.35. The result is sensitive to the CO-to-H2, conversion factor adopted, and both gas fraction and SFE have comparable evolution, when we adopt the low starburst conversion factor of α=0.8 Mo/(K km/s pc^2). Adopting a higher αwill increase the role of the gas fraction. Using α=0.8, the SFE and the gas fraction for z=0.2-1.0 ULIRGs are found to be significantly higher, by a factor 3, than for local ULIRGs, and are comparable to high redshift ones. We compare this evolution to the expected cosmic H2 abundance and the cosmic star formation history.

preprint2012arXiv

The anatomy of an extreme starburst within 1.3Gyr of the Big Bang revealed by ALMA

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

preprint2012arXiv

The Evolving Interstellar Medium of Star Forming Galaxies Since z=2 as Probed by Their Infrared Spectral Energy Distributions

Using data from the mid-infrared to millimeter wavelengths for individual galaxies and for stacked ensembles at 0.5<z<2, we derive robust estimates of dust masses (Mdust) for main sequence (MS) galaxies, which obey a tight correlation between star formation rate (SFR) and stellar mass (M*), and for star-bursting galaxies that fall outside that relation. Exploiting the correlation of gas to dust mass with metallicity (Mgas/Mdust -Z), we use our measurements to constrain the gas content, CO-to-H2 conversion factors (a_co) and star formation efficiencies (SFE) of these distant galaxies. Using large statistical samples, we confirm that a_co and SFE are an order of magnitude higher and lower, respectively, in MS galaxies at high redshifts compared to the values of local galaxies with equivalently high infrared luminosities. For galaxies within the MS, we show that the variations of specific star formation rates (sSFR=SFR/M*) are driven by varying gas fractions. For relatively massive galaxies like those in our samples, we show that the hardness of the radiation field, <U>, which is proportional to the dust mass weighted luminosity (LIR/Mdust), and the primary parameter defining the shape of the SED, is equivalent to SFE/Z. For MS galaxies we measure this quantity, <U>, showing that it does not depend significantly on either the stellar mass or the sSFR. This is explained as a simple consequence of the existing correlations between SFR-M*, M*-Z and Mgas-SFR. Instead, we show that <U> (or LIR/Mdust) does evolve, with MS galaxies having harder radiation fields and thus warmer temperatures as redshift increases from z=0 to 2, a trend which can also be understood based on the redshift evolution of the M*-Z and SFR-M* relations. These results motivate the construction of a universal set of SED templates for MS galaxies which vary as a function of redshift with only one parameter, <U>.

preprint2012arXiv

Variations in the Fundamental constants in the QSO Host J1148+5251 at z = 6.4 and the BR1202-0725 System at z = 4.7

We use sensitive observations of three high redshift sources; [CII] fine structure and CO(2-1) rotational transitions for the z=6.4 Quasar host galaxy (QSO) J1148+5251, and [CII] and CO(5-4) transitions from the QSO BR1202-0725 and its sub-millimeter companion (SMG) galaxy at z=4.7. We use these observations to place constraints on the quantity Dz = z(CO) - z(CII) for each source where z(CO) and z(CII) are the observed redshifts of the CO rotational transition and [CII] fine structure transition respectively, using a combination of approaches; 1) By modelling the emission line profiles using `shapelets' to compare both the emission redshifts and the line profiles themselves, in order to make inferences about the intrinsic velocity differences between the molecular and atomic gas, and 2) By performing a marginalisation over all model parameters in order to calculate a non-parametric estimate of Dz. We derive 99% confidence intervals for the marginalised posterior of Dz of (-1.9 pm 1.3) x10^-3, (-3 pm 8) x10^-4 and (-2 pm 4) x10^-3 for J1148+5251, and the BR1202-0725 QSO and SMG respectively. We show the [CII] and CO(2-1) line profiles for J1148+5251 are consistent with each other within the limits of the data, whilst the [CII] and CO(5-4) line profiles from the BR1202-0725 QSO and SMG respectively have 65 and >99.9% probabilities of being inconsistent, with the CO(5-4) lines ~ 30% wider than the [CII] lines. Therefore whilst the observed values of Dz can correspond to variations in the quantity Delta F/F with cosmic time, where F=alpha^2/mu, with alpha the fine structure constant, and mu the proton-to-electron mass ratio, of both (-3.3 pm 2.3) x10^-4 for a look back time of 12.9 Gyr and of (-5 pm 15) x10^-5 for a look back time of 12.4 Gyr we propose that they are the result of the two species of gas being spatially separated as indicated by the inconsistencies in their line profiles.

preprint2011arXiv

A Constant Molecular Gas Depletion Time in Nearby Disk Galaxies

We combine new sensitive, wide-field CO data from the HERACLES survey with ultraviolet and infrared data from GALEX and Spitzer to compare the surface densities of H2, Sigma_H2, and recent star formation rate, Sigma_SFR, over many thousands of positions in 30 nearby disk galaxies. We more than quadruple the size of the galaxy sample compared to previous work and include targets with a wide range of galaxy properties. Even though the disk galaxies in this study span a wide range of properties, we find a strong and approximately linear correlation between Sigma_SFR and Sigma_H2 at our common resolution of 1kpc. This implies a roughly constant median H2 consumption time, tau_H2 = Sigma_H2 / Sigma_SFR, of ~2.35Gyr (including heavy elements) across our sample. At 1kpc resolution, there is only a weak correlation between Sigma_H2 and tau_H2 over the range Sigma_H2~5-100M_sun/pc^2, which is probed by our data. We compile a broad set of literature measurements that have been obtained using a variety of star formation tracers, sampling schemes and physical scales and show that overall, these data yield almost exactly the same results, although with more scatter. We interpret these results as strong, albeit indirect evidence that star formation proceeds in a uniform way in giant molecular clouds in the disks of spiral galaxies.

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

CO in GN20: The Nature of a z=4 Submillimeter Galaxy

We present a study of the formation of clustered, massive galaxies at large look--back times, via high resolution spectroscopic imaging of CO in the unique GN20 proto-cluster. The data reveal evidence for rich structure and gas dynamics in unprecedented detail, allowing us to image the molecular gas with a resolution of only 1 kpc just 1.5 Gyr after the Big Bang. These state-of-the-art data give new insight into the detailed physical processes involved in early massive galaxy formation, and they provide a first glimpse of the morphological studies that will become feasible on a regular basis with ALMA.

preprint2011arXiv

EVLA observations of a proto-cluster of molecular gas rich galaxies at z = 4.05

We present observations of the molecular gas in the GN20 proto-cluster of galaxies at $z =4.05$ using the Expanded Very Large Array (EVLA). This group of galaxies is the ideal laboratory for studying the formation of massive galaxies via luminous, gas-rich starbursts within 1.6 Gyr of the Big Bang. We detect three galaxies in the proto-cluster in CO 2-1 emission, with gas masses (H$_2$) between $10^{10}$ and $10^{11} \times (α/0.8)$ M$_\odot$. The emission from the brightest source, GN20, is resolved with a size $\sim 2"$, and has a clear north-south velocity gradient, possibly indicating ordered rotation. The gas mass in GN20 is comparable to the stellar mass ($1.3\times 10^{11} \times (α/0.8)$ M$_\odot$ and $2.3\times 10^{11}$ M$_\odot$, respectively), and the sum of gas plus stellar mass is comparable to the dynamical mass of the system ($\sim 3.4\times 10^{11} [sin(i)/sin(45^o)]^{-2}$ M$_\odot$), within a 5kpc radius. There is also evidence for a tidal tail extending another $2"$ north of the galaxy with a narrow velocity dispersion. GN20 may be a massive, gas rich disk that is gravitationally disturbed, but not completely disrupted. There is one Lyman-break galaxy (BD29079) in the GN20 proto-cluster with an optical spectroscopic redshift within our search volume, and we set a 3$σ$ limit to the molecular gas mass of this galaxy of $1.1\times 10^{10} \times (α/0.8)$ M$_\odot$.

preprint2011arXiv

Galaxy Evolution and Star Formation Efficiency at 0.2 < z < 0.6

We present the results of a CO line survey of 30 galaxies at moderate redshift (z \sim 0.2-0.6), with the IRAM 30m telescope, with the goal to follow galaxy evolution and in particular the star formation efficiency (SFE) as defined by the ratio between far-infrared luminosity and molecular gas mass (LFIR/M(H2)). The sources are selected to be ultra-luminous infrared galaxies (ULIRGs), with LFIR larger than 2.8 10^{12} Lsol, experiencing starbursts; adopting a low ULIRG CO-to-H2 conversion factor, their gas consumption time-scale is lower than 10^8 yr. To date only very few CO observations exist in this redshift range that spans nearly 25% of the universe's age. Considerable evolution of the star formation rate is already observed during this period. 18 galaxies out of our sample of 30 are detected (of which 16 are new detections), corresponding to a detection rate of 60%. The average CO luminosity for the 18 galaxies detected is L'CO = 2 10^{10} K km/s pc^2, corresponding to an average H2 mass of 1.6 10^{10} Msol. The FIR luminosity correlates well with the CO luminosity, in agreement with the correlation found for low and high redshift ULIRGs. Although the conversion factor between CO luminosity and H2 mass is uncertain, we find that the maximum amount of gas available for a single galaxy is quickly increasing as a function of redshift. Using the same conversion factor, the SFEs for z\sim 0.2-0.6 ULIRGs are found to be significantly higher, by a factor 3, than for local ULIRGs, and are comparable to high redshift ones. We compare this evolution to the expected cosmic H2 abundance and the cosmic star formation history.

preprint2011arXiv

GOODS-Herschel~: Gas-to-dust mass ratios and CO-to-H_2 conversion factors in normal and starbursting galaxies at high-z

We explore the gas-to-dust mass ratio (G/D) and the CO luminosity-to-Mgas conversion factor (a_co) of two well studied galaxies in the GOODS-N field, that are expected to have different star forming modes, the starburst GN20 at z=4.05 and the normal star-forming galaxy BzK-21000 at z=1.52. Detailed sampling is available for their Rayleigh-Jeans emission via ground based mm interferometry (1.1-6.6mm) along with Herschel, PACS and SPIRE data that probe the peak of their infrared emission. Using the physically motivated Draine & Li (2007) models, as well as a modified black body function, we measure the dust mass (Md) of the sources and find 2.0^{+0.7}_{-0.6} x 10^{9} Msun for GN20 and 8.6^{+0.6}_{-0.9} x 10^{8} Msun for BzK-21000. The addition of mm data reduces the uncertainties of the derived Md by a factor of ~2, allowing the use of the local G\D vs metallicity relation to place constraints on the a_co values of the two sources. For GN20 we derive a conversion factor of a_co < 1.0 Msun pc^{-2}(K km s^{-1})^{-1}, consistent with that of local ULIRGs, while for BzK-21000 we find a considerably higher value, a_co ~4.0 Msun pc^{-2}(K km s^{-1})^{-1}, in agreement with an independent kinematic derivation reported previously. The implied star formation efficiency is ~25 Lsun/Msun for BzK-21000, a factor of ~5-10 lower than that of GN20. The findings for these two sources support the existence of different disk-like and starburst star-formation modes in distant galaxies, although a larger sample is required to draw statistically robust results

preprint2011arXiv

Herschel Far-Infrared and Sub-millimeter Photometry for the KINGFISH Sample of Nearby Galaxies

New far-infrared and sub-millimeter photometry from the Herschel Space Observatory is presented for 61 nearby galaxies from the Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel (KINGFISH) sample. The spatially-integrated fluxes are largely consistent with expectations based on Spitzer far-infrared photometry and extrapolations to longer wavelengths using popular dust emission models. Dwarf irregular galaxies are notable exceptions, as already noted by other authors, as their 500um emission shows evidence for a sub-millimeter excess. In addition, the fraction of dust heating attributed to intense radiation fields associated with photo-dissociation regions is found to be (21+/-4)% larger when Herschel data are included in the analysis. Dust masses obtained from the dust emission models of Draine & Li are found to be on average nearly a factor of two higher than those based on single-temperature modified blackbodies, as single blackbody curves do not capture the full range of dust temperatures inherent to any galaxy. The discrepancy is largest for galaxies exhibiting the coolest far-infrared colors.

preprint2011arXiv

KINGFISH -- Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel: Survey Description and Image Atlas

The KINGFISH project (Key Insights on Nearby Galaxies: a Far-Infrared Survey with Herschel) is an imaging and spectroscopic survey of 61 nearby (d < 30 Mpc) galaxies, chosen to cover a wide range of galaxy properties and local interstellar medium (ISM) environments found in the nearby Universe. Its broad goals are to characterize the ISM of present-day galaxies, the heating and cooling of their gaseous and dust components, and to better understand the physical processes linking star formation and the ISM. KINGFISH is a direct descendant of the Spitzer Infrared Nearby Galaxies Survey (SINGS), which produced complete Spitzer imaging and spectroscopic mapping and a comprehensive set of multi-wavelength ancillary observations for the sample. The Herschel imaging consists of complete maps for the galaxies at 70, 100, 160, 250, 350, and 500 microns. The spectal line imaging of the principal atomic ISM cooling lines ([OI]63um, [OIII]88um, [NII]122,205um, and [CII]158um) covers the subregions in the centers and disks that already have been mapped in the mid-infrared with Spitzer. The KINGFISH and SINGS multi-wavelength datasets combined provide panchromatic mapping of the galaxies sufficient to resolve individual star-forming regions, and tracing the important heating and cooling channels of the ISM, across a wide range of local extragalactic ISM environments. This paper summarizes the scientific strategy for KINGFISH, the properties of the galaxy sample, the observing strategy, and data processing and products. It also presents a combined Spitzer and Herschel image atlas for the KINGFISH galaxies, covering the wavelength range 3.6 -- 500 microns. All imaging and spectroscopy data products will be released to the Herschel user generated product archives.

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

The BLAST 250-micron selected galaxy population in GOODS-South

We identify and investigate the nature of the 20 brightest 250-micron sources detected by the BLAST within the central 150 sq. arcmin of the GOODS-South field. Aided by the available deep VLA radio imaging, reaching S_1.4 = 30 micro-Jy, we have identified radio counterparts for 17/20 of the 250-micron sources. The resulting enhanced positional accuracy of ~1 arcsec has then allowed us to exploit the deep multi-frequency imaging of GOODS-South to establish secure galaxy counterparts for the 17 radio-identified sources, and plausible galaxy candidates for the 3 radio-unidentified sources. Confusion is a serious issue for this deep BLAST 250-micron survey, due to the large size of the beam. Nevertheless, we argue that our chosen counterparts are significant, and often dominant contributors to the measured BLAST flux densities. For all of these 20 galaxies we have been able to determine spectroscopic (8) or photometric (12) redshifts. The result is the first `complete' redshift distribution for a deep 250-micron selected galaxy sample. This reveals that 250-micron surveys reaching detection limits of ~30 mJy contain not only low-redshift spirals/LIRGs, but also the extreme z~2 dust-enshrouded starburst galaxies previously discovered at sub-millimetre wavelengths. Based on their IRAC colours, we find that virtually all of the BLAST galaxy identifications appear better described as analogues of the M82 starburst galaxy, or Sc star-forming discs rather than highly obscured ULIRGs. Inspection of the LABOCA 870-micron imaging of the GOODS-South field yields detections of 7/11 of the z>1 BLAST sources, and reveals 250/870 flux-density ratios consistent with a standard 40K modified black-body fit with a dust emissivity index beta=1.5.

preprint2010arXiv

Black hole accretion and star formation as drivers of gas excitation and chemistry in Mrk231

We present a full high resolution SPIRE FTS spectrum of the nearby ultraluminous infrared galaxy Mrk231. In total 25 lines are detected, including CO J=5-4 through J=13-12, 7 rotational lines of H2O, 3 of OH+ and one line each of H2O+, CH+, and HF. We find that the excitation of the CO rotational levels up to J=8 can be accounted for by UV radiation from star formation. However, the approximately flat luminosity distribution of the CO lines over the rotational ladder above J=8 requires the presence of a separate source of excitation for the highest CO lines. We explore X-ray heating by the accreting supermassive black hole in Mrk231 as a source of excitation for these lines, and find that it can reproduce the observed luminosities. We also consider a model with dense gas in a strong UV radiation field to produce the highest CO lines, but find that this model strongly overpredicts the hot dust mass in Mrk231. Our favoured model consists of a star forming disk of radius 560 pc, containing clumps of dense gas exposed to strong UV radiation, dominating the emission of CO lines up to J=8. X-rays from the accreting supermassive black hole in Mrk231 dominate the excitation and chemistry of the inner disk out to a radius of 160 pc, consistent with the X-ray power of the AGN in Mrk231. The extraordinary luminosity of the OH+ and H2O+ lines reveals the signature of X-ray driven excitation and chemistry in this region.

preprint2010arXiv

Different star formation laws for disks versus starbursts at low and high redshifts

We present evidence that 'bona fide' disks and starburst systems occupy distinct regions in the gas mass versus star formation (SF) rate plane, both for the integrated quantities and for the respective surface densities. This result is based on CO observations of galaxy populations at low and high redshifts, and on the current consensus for the CO luminosity to gas mass conversion factors. The data suggest the existence of two different star formation regimes: a long-lasting mode for disks and a more rapid mode for starbursts, the latter probably occurring during major mergers or in dense nuclear SF regions. Both modes are observable over a large range of SF rates. The detection of CO emission from distant near-IR selected galaxies reveals such bimodal behavior for the first time, as they allow us to probe gas in disk galaxies with much higher SF rates than are seen locally. The different regimes can potentially be interpreted as the effect of a top-heavy IMF in starbursts. However, we favor a different physical origin related to the fraction of molecular gas in dense clouds. The IR luminosity to gas mass ratio (i.e., the SF efficiency) appears to be inversely proportional to the dynamical (rotation) timescale. Only when accounting for the dynamical timescale, a universal SF law is obtained, suggesting a direct link between global galaxy properties and the local SF rate.

preprint2010arXiv

Enhanced dust heating in the bulges of early-type spiral galaxies

Stellar density and bar strength should affect the temperatures of the cool (T ~ 20-30 K) dust component in the inner regions of galaxies, which implies that the ratio of temperatures in the circumnuclear regions to the disk should depend on Hubble type. We investigate the differences between cool dust temperatures in the central 3 kpc and disk of 13 nearby galaxies by fitting models to measurements between 70 and 500 microns. We attempt to quantify temperature trends in nearby disk galaxies, with archival data from Spitzer/MIPS and new observations with Herschel/SPIRE, which were acquired during the first phases of the Herschel observations for the KINGFISH (key insights in nearby galaxies: a far-infrared survey with Herschel) sample. We fit single-temperature modified blackbodies to far-infrared and submillimeter measurements of the central and disk regions of galaxies to determine the temperature of the component(s) emitting at those wavelengths. We present the ratio of central-region-to-disk-temperatures of the cool dust component of 13 nearby galaxies as a function of morphological type. We find a significant temperature gradient in the cool dust component in all galaxies, with a mean center-to-disk temperature ratio of 1.15 +/- 0.03. The cool dust temperatures in the central ~3 kpc of nearby galaxies are 23(+/-3)% hotter for morphological types earlier than Sc, and only 9(+/-3)% hotter for later types. The temperature ratio is also correlated with bar strength, with only strongly barred galaxies having a ratio over 1.2. The strong radiation field in the high stellar density of a galactic bulge tends to heat the cool dust component to higher temperatures, at least in early-type spirals with relatively large bulges, especially when paired with a strong bar.

preprint2010arXiv

Far-Infrared Line Imaging of the Starburst Ring in NGC 1097 with the Herschel/PACS Spectrometer

NGC 1097 is a nearby SBb galaxy with a Seyfert nucleus and a bright starburst ring. We study the physical properties of the interstellar medium (ISM) in the ring using spatially resolved far-infrared spectral maps of the circumnuclear starburst ring of NGC 1097, obtained with the PACS spectrometer on board the Herschel Space Telescope. In particular, we map the important ISM cooling and diagnostic emission lines of [OI] 63 $μ$m, [OIII] 88 $μ$m, [NII] 122 $μ$m, [CII] 158 $μ$m and [NII] 205 $μ$m. We observe that in the [OI] 63 $μ$m, [OIII] 88 $μ$m, and [NII] 122 $μ$m line maps, the emission is enhanced in clumps along the NE part of the ring. We observe evidence of rapid rotation in the circumnuclear ring, with a rotation velocity of ~220$ km s$^{-1}$ (inclination uncorrected) measured in all lines. The [OI] 63 $μ$m/[CII] 158 $μ$m ratio varies smoothly throughout the central region, and is enhanced on the northeastern part of the ring, which may indicate a stronger radiation field. This enhancement coincides with peaks in the [OI] 63 $μ$m and [OIII] 88 $μ$m maps. Variations of the [NII] 122 $μ$m/[NII] 205 $μ$m ratio correspond to a range in the ionized gas density between 150 and 400 cm$^{-3}$.

preprint2010arXiv

Herschel observations of water vapour in Markarian 231

The Ultra Luminous InfraRed Galaxy Mrk 231 reveals up to seven rotational lines of water (H2O) in emission, including a very high-lying (E_{upper}=640 K) line detected at a 4sigma level, within the Herschel/SPIRE wavelength range, whereas PACS observations show one H2O line at 78 microns in absorption, as found for other H2O lines previously detected by ISO. The absorption/emission dichotomy is caused by the pumping of the rotational levels by far-infrared radiation emitted by dust, and subsequent relaxation through lines at longer wavelengths, which allows us to estimate both the column density of H2O and the general characteristics of the underlying far-infrared continuum source. Radiative transfer models including excitation through both absorption of far-infrared radiation emitted by dust and collisions are used to calculate the equilibrium level populations of H2O and the corresponding line fluxes. The highest-lying H2O lines detected in emission, with levels at 300-640 K above the ground state, indicate that the source of far-infrared radiation responsible for the pumping is compact (radius=110-180 pc) and warm (T_{dust}=85-95 K), accounting for at least 45% of the bolometric luminosity. The high column density, N(H2O)~5x10^{17} cm^{-2}, found in this nuclear component, is most probably the consequence of shocks/cosmic rays, an XDR chemistry, and/or an "undepleted chemistry" where grain mantles are evaporated. A more extended region, presumably the inner region of the 1-kpc disk observed in other molecular species, could contribute to the flux observed in low-lying H2O lines through dense hot cores, and/or shocks. The H2O 78 micron line observed with PACS shows hints of a blue-shifted wing seen in absorption, possibly indicating the occurrence of H2O in the prominent outflow detected in OH (Fischer et al., this volume).

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

Imaging the molecular gas in a submm galaxy at z = 4.05: cold mode accretion or a major merger?

We present a high resolution (down to 0.18"), multi-transition imaging study of the molecular gas in the z = 4.05 submillimeter galaxy GN20. GN20 is one of the most luminous starburst galaxy known at z > 4, and is a member of a rich proto-cluster of galaxies at z = 4.05 in GOODS-North. We have observed the CO 1-0 and 2-1 emission with the VLA, the CO 6-5 emission with the PdBI Interferometer, and the 5-4 emission with CARMA. The H_2 mass derived from the CO 1-0 emission is 1.3 \times 10^{11} (α/0.8) Mo. High resolution imaging of CO 2-1 shows emission distributed over a large area, appearing as partial ring, or disk, of ~ 10kpc diameter. The integrated CO excitation is higher than found in the inner disk of the Milky Way, but lower than that seen in high redshift quasar host galaxies and low redshift starburst nuclei. The VLA CO 2-1 image at 0.2" resolution shows resolved, clumpy structure, with a few brighter clumps with intrinsic sizes ~ 2 kpc. The velocity field determined from the CO 6-5 emission is consistent with a rotating disk with a rotation velocity of ~ 570 km s^{-1} (using an inclination angle of 45^o), from which we derive a dynamical mass of 3 \times 10^{11} \msun within about 4 kpc radius. The star formation distribution, as derived from imaging of the radio synchrotron and dust continuum, is on a similar scale as the molecular gas distribution. The molecular gas and star formation are offset by ~ 1" from the HST I-band emission, implying that the regions of most intense star formation are highly dust-obscured on a scale of ~ 10 kpc. The large spatial extent and ordered rotation of this object suggests that this is not a major merger, but rather a clumpy disk accreting gas rapidly in minor mergers or smoothly from the proto-intracluster medium. ABSTRACT TRUNCATED

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

Mapping far-IR emission from the central kiloparsec of NGC 1097

Using photometry of NGC 1097 from the Herschel PACS (Photodetector Array Camera and Spectrometer) instrument, we study the resolved properties of thermal dust continuum emission from a circumnuclear starburst ring with a radius ~ 900 pc. These observations are the first to resolve the structure of a circumnuclear ring at wavelengths that probe the peak (i.e. lambda ~ 100 micron) of the dust spectral energy distribution. The ring dominates the far-infrared (far-IR) emission from the galaxy - the high angular resolution of PACS allows us to isolate the ring's contribution and we find it is responsible for 75, 60 and 55% of the total flux of NGC 1097 at 70, 100 and 160 micron, respectively. We compare the far-IR structure of the ring to what is seen at other wavelengths and identify a sequence of far-IR bright knots that correspond to those seen in radio and mid-IR images. The mid- and far-IR band ratios in the ring vary by less than +/- 20% azimuthally, indicating modest variation in the radiation field heating the dust on ~ 600 pc scales. We explore various explanations for the azimuthal uniformity in the far-IR colors of the ring including a lack of well-defined age gradients in the young stellar cluster population, a dominant contribution to the far-IR emission from dust heated by older (> 10 Myr) stars and/or a quick smoothing of local enhancements in dust temperature due to the short orbital period of the ring. Finally, we improve previous limits on the far-IR flux from the inner ~ 600 pc of NGC 1097 by an order of magnitude, providing a better estimate of the total bolometric emission arising from the active galactic nucleus and its associated central starburst.

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

Scaling Relations between Gas and Star Formation in Nearby Galaxies

High resolution, multi-wavelength maps of a sizeable set of nearby galaxies have made it possible to study how the surface densities of HI, H2 and star formation rate (Sigma_HI, Sigma_H2, Sigma_SFR) relate on scales of a few hundred parsecs. At these scales, individual galaxy disks are comfortably resolved, making it possible to assess gas-SFR relations with respect to environment within galaxies. Sigma_H2, traced by CO intensity, shows a strong correlation with Sigma_SFR and the ratio between these two quantities, the molecular gas depletion time, appears to be constant at about 2Gyr in large spiral galaxies. Within the star-forming disks of galaxies, Sigma_SFR shows almost no correlation with Sigma_HI. In the outer parts of galaxies, however, Sigma_SFR does scale with Sigma_HI, though with large scatter. Combining data from these different environments yields a distribution with multiple regimes in Sigma_gas - Sigma_SFR space. If the underlying assumptions to convert observables to physical quantities are matched, even combined datasets based on different SFR tracers, methodologies and spatial scales occupy a well define locus in Sigma_gas - Sigma_SFR space.

preprint2010arXiv

The Calibration of Monochromatic Far-Infrared Star Formation Rate Indicators

(Abridged) Spitzer data at 24, 70, and 160 micron and ground-based H-alpha images are analyzed for a sample of 189 nearby star-forming and starburst galaxies to investigate whether reliable star formation rate (SFR) indicators can be defined using the monochromatic infrared dust emission centered at 70 and 160 micron. We compare recently published recipes for SFR measures using combinations of the 24 micron and observed H-alpha luminosities with those using 24 micron luminosity alone. From these comparisons, we derive a reference SFR indicator for use in our analysis. Linear correlations between SFR and the 70 and 160 micron luminosity are found for L(70)>=1.4x10^{42} erg/s and L(160)>=2x10^{42} erg/s, corresponding to SFR>=0.1-0.3 M_sun/yr. Below those two luminosity limits, the relation between SFR and 70 micron (160 micron) luminosity is non-linear and SFR calibrations become problematic. The dispersion of the data around the mean trend increases for increasing wavelength, becoming about 25% (factor ~2) larger at 70 (160) micron than at 24 micron. The increasing dispersion is likely an effect of the increasing contribution to the infrared emission of dust heated by stellar populations not associated with the current star formation. The non-linear relation between SFR and the 70 and 160 micron emission at faint galaxy luminosities suggests that the increasing transparency of the interstellar medium, decreasing effective dust temperature, and decreasing filling factor of star forming regions across the galaxy become important factors for decreasing luminosity. The SFR calibrations are provided for galaxies with oxygen abundance 12+Log(O/H)>8.1. At lower metallicity the infrared luminosity no longer reliably traces the SFR because galaxies are less dusty and more transparent.

preprint2010arXiv

The evolution of stellar structures in dwarf galaxies

We present a study of the variation of spatial structure of stellar populations within dwarf galaxies as a function of the population age. We use deep Hubble Space Telescope/Advanced Camera for Surveys imaging of nearby dwarf galaxies in order to resolve individual stars and create composite colour-magnitude diagrams (CMDs) for each galaxy. Using the obtained CMDs, we select Blue Helium Burning stars (BHeBs), which can be unambiguously age-dated by comparing the absolute magnitude of individual stars with stellar isochrones. Additionally, we select a very young (<10 Myr) population of OB stars for a subset of the galaxies based on the tip of the young main-sequence. By selecting stars in different age ranges we can then study how the spatial distribution of these stars evolves with time. We find, in agreement with previous studies, that stars are born within galaxies with a high degree of substructure which is made up of a continuous distribution of clusters, groups and associations from parsec to hundreds of parsec scales. These structures disperse on timescales of tens to hundreds of Myr, which we quantify using the two-point correlation function and the Q-parameter developed by Cartwright & Whitworth (2004). On galactic scales, we can place lower limits on the time it takes to remove the original structure (i.e., structure survives for at least this long), tevo, which varies between ~100~Myr (NGC~2366) and ~350 Myr (DDO~165). This is similar to what we have found previously for the SMC (~80~Myr) and the LMC (~175 Myr). We do not find any strong correlations between tevo and the luminosity of the host galaxy.

preprint2010arXiv

The Fine-Scale Structure of the neutral Interstellar Medium in nearby Galaxies

We present an analysis of the properties of HI holes detected in 20 galaxies that are part of "The HI Nearby Galaxy Survey" (THINGS). We detected more than 1000 holes in total in the sampled galaxies. Where they can be measured, their sizes range from about 100 pc (our resolution limit) to about 2 kpc, their expansion velocities range from 4 to 36 km/s, and their ages are estimated to range between 3 and 150 Myr. The holes are found throughout the disks of the galaxies, out to the edge of the HI; 23% of the holes fall outside R25. We find that shear limits the age of holes in spirals (shear is less important in dwarf galaxies) which explains why HI holes in dwarfs are rounder, on average than in spirals. Shear, which is particularly strong in the inner part of spiral galaxies, also explains why we find that holes outside R25 are larger and older. We derive the scale height of the HI disk as a function of galactocentric radius and find that the disk flares up in all galaxies. We proceed to derive the surface and volume porosity (Q2D and Q3D) and find that this correlates with the type of the host galaxy: later Hubble types tend to be more porous. The size distribution of the holes in our sample follows a power law with a slope of a ~ -2.9. Assuming that the holes are the result of massive star formation, we derive values for the supernova rate (SNR) and star formation rate (SFR) which scales with the SFR derived based on other tracers. If we extrapolate the observed number of holes to include those that fall below our resolution limit, down to holes created by a single supernova, we find that our results are compatible with the hypothesis that HI holes result from star formation.

preprint2010arXiv

The LABOCA survey of the Extended Chandra Deep Field South - radio and mid-infrared counterparts to submillimetre galaxies

We present radio and infrared (3.6-24um) counterparts to submillimetre galaxies (SMGs) detected in the Extended Chandra Deep Field South with the LABOCA 870-um bolometer camera on the 12-m Atacama Pathfinder Experiment. Using the Very Large Array at 1.4GHz and Spitzer we have identified secure counterparts to 79 of the 126 SMGs (SNR>3.7, S_870>4.4mJy) in the field, 62 via their radio and/or 24-um emission, the remainder using a colour-flux cut on IRAC 3.6- and 5.8-um sources chosen to maximise the number of secure, coincident radio and 24-um counterparts. In constructing our radio catalogue, we have corrected for the effects of `flux boosting', then used the corrected flux densities to estimate the redshifts of the SMGs based on the radio/submm spectral indices. The effect of the boosting correction is to increase the median redshift by 0.2 resulting in a value of z=2.2 (+0.7-0.8) (1-sigma errors) for the secure radio counterparts, in agreement with other studies, both spectroscopic and photometric.

preprint2010arXiv

The LABOCA survey of the Extended Chandra Deep Field South: Two modes of star formation in AGN hosts?

We study the co-existence of star formation and AGN activity in X-ray selected AGN by analyzing stacked 870um submm emission from a deep and wide map of the Extended Chandra Deep Field South, obtained with LABOCA at the APEX telescope. The total X-ray sample of 895 sources with median redshift z~1 is detected at a mean submm flux of 0.49+-0.04mJy, corresponding to a typical star formation rate around 30Msun/yr for a T=35K, beta=1.5 greybody far-infrared SED. The good S/N permits stacking analyses for subgroups. We observe a trend of star formation rate increasing with redshift. An increase of star formation rate with AGN luminosity is indicated at the highest L_2-10>~1E44erg/s luminosities only. Increasing trends with X-ray obscuration as expected in some AGN evolutionary scenarios are not observed for the bulk of the X-ray AGN sample but may be present for the highest intrinsic luminosity objects. This suggests a transition between two modes in the coexistence of AGN activity and star formation. For the bulk of the sample, the X-ray luminosity and obscuration of the AGN are not intimately linked to the global star formation rate of their hosts. The hosts are likely massive and forming stars secularly, at rates similar to the pervasive star formation seen in massive galaxies without an AGN at similar redshifts. The change indicated towards more intense star formation, and a more pronounced increase in star formation rates between unobscured and obscured AGN at highest luminosities suggests that luminous AGN follow an evolutionary path on which obscured AGN activity and intense star formation are linked, possibly via merging. Comparison to local hard X-ray selected AGN supports this interpretation. [Abridged]

preprint2010arXiv

Tightly Correlated HI and FUV Emission in the Outskirts of M83

We compare sensitive HI data from The HI Nearby Galaxy Survey (THINGS) and deep far UV (FUV) data from GALEX in the outer disk of M83. The FUV and HI maps show a stunning spatial correlation out to almost 4 optical radii (r25), roughly the extent of our maps. This underscores that HI traces the gas reservoir for outer disk star formation and it implies that massive (at least low level) star formation proceeds almost everywhere HI is observed. Whereas the average FUV intensity decreases steadily with increasing radius before leveling off at ~1.7 r25, the decline in HI surface density is more subtle. Low HI columns (<2 M_solar/pc^2) contribute most of the mass in the outer disk, which is not the case within r25. The time for star formation to consume the available HI, inferred from the ratio of HI to FUV intensity, rises with increasing radius before leveling off at ~100 Gyr, i.e., many Hubble times, near ~1.7 r25. Assuming the relatively short H2 depletion times observed in the inner parts of galaxies hold in outer disks, the conversion of HI into bound, molecular clouds seems to limit star formation in outer galaxy disks. The long consumption times suggest that most of the extended HI observed in M83 will not be consumed by in situ star formation. However, even these low star formation rates are enough to expect moderate chemical enrichment in a closed outer disk.

preprint2010arXiv

Transit timing variation and activity in the WASP-10 planetary system

Transit timing analysis may be an effective method of discovering additional bodies in extrasolar systems which harbour transiting exoplanets. The deviations from the Keplerian motion, caused by mutual gravitational interactions between planets, are expected to generate transit timing variations of transiting exoplanets. In 2009 we collected 9 light curves of 8 transits of the exoplanet WASP-10b. Combining these data with published ones, we found that transit timing cannot be explained by a constant period but by a periodic variation. Simplified three-body models which reproduce the observed variations of timing residuals were identified by numerical simulations. We found that the configuration with an additional planet of mass of $\sim$0.1 $M_{\rm{J}}$ and orbital period of $\sim$5.23 d, located close to the outer 5:3 mean motion resonance, is the most likely scenario. If the second planet is a transiter, the estimated flux drop will be $\sim$0.3 per cent and can be observable with a ground-based telescope. Moreover, we present evidence that the spots on the stellar surface and rotation of the star affect the radial velocity curve giving rise to spurious eccentricity of the orbit of the first planet. We argue that the orbit of WASP-10b is essentially circular. Using the gyrochronology method, the host star was found to be $270 \pm 80$ Myr old. This young age can explain the large radius reported for WASP-10b.

preprint2010arXiv

Unusually Luminous Giant Molecular Clouds in the Outer Disk of M33

We use high spatial resolution (~7pc) CARMA observations to derive detailed properties for 8 giant molecular clouds (GMCs) at a galactocentric radius corresponding to approximately two CO scale lengths, or ~0.5 optical radii (r25), in the Local Group spiral galaxy M33. At this radius, molecular gas fraction, dust-to-gas ratio and metallicity are much lower than in the inner part of M33 or in a typical spiral galaxy. This allows us to probe the impact of environment on GMC properties by comparing our measurements to previous data from the inner disk of M33, the Milky Way and other nearby galaxies. The outer disk clouds roughly fall on the size-linewidth relation defined by extragalactic GMCs, but are slightly displaced from the luminosity-virial mass relation in the sense of having high CO luminosity compared to the inferred virial mass. This implies a different CO-to-H2 conversion factor, which is on average a factor of two lower than the inner disk and the extragalactic average. We attribute this to significantly higher measured brightness temperatures of the outer disk clouds compared to the ancillary sample of GMCs, which is likely an effect of enhanced radiation levels due to massive star formation in the vicinity of our target field. Apart from brightness temperature, the properties we determine for the outer disk GMCs in M33 do not differ significantly from those of our comparison sample. In particular, the combined sample of inner and outer disk M33 clouds covers roughly the same range in size, linewidth, virial mass and CO luminosity than the sample of Milky Way GMCs. When compared to the inner disk clouds in M33, however, we find even the brightest outer disk clouds to be smaller than most of their inner disk counterparts. This may be due to incomplete sampling or a potentially steeper cloud mass function at larger radii.

preprint2010arXiv

Very High Gas Fractions and Extended Gas Reservoirs in z=1.5 Disk Galaxies

We present evidence for very high gas fractions and extended molecular gas reservoirs in normal, near-infrared selected (BzK) galaxies at z~1.5, based on multi-configuration CO[2-1] observations obtained at the IRAM PdBI. Six of the six galaxies observed were securely detected. High resolution observations resolve the CO emission in four of them, implying sizes of order of 6-11 kpc and suggesting the presence of rotation. The UV morphologies are consistent with clumpy, unstable disks, and the UV sizes are consistent with the CO sizes. The star formation efficiencies are homogeneously low and similar to local spirals - the resulting gas depletion times are ~0.5 Gyr, much higher than what is seen in high-z submm galaxies and quasars. The CO luminosities can be predicted to within 0.15 dex from the star formation rates and stellar masses, implying a tight correlation of the gas mass with these quantities. We use dynamical models of clumpy disk galaxies to derive dynamical masses. These models are able to reproduce the peculiar spectral line shapes of the CO emission. After accounting for the stellar and dark matter masses we derive gas masses of 0.4-1.2x10^11 Msun. The conversion factor is very high: alpha_CO=3.6+-0.8, consistent with the Galaxy but four times higher than that of local ultra-luminous IR galaxies. The gas accounts for an impressive 50-65% of the baryons within the galaxies' half light radii. We are witnessing truly gas-dominated galaxies at z~1.5, a finding that explains the high specific SFRs observed for z>1 galaxies. The BzK galaxies can be viewed as scaled-up versions of local disk galaxies, with low efficiency star formation taking place inside extended, low excitation gas disks. They are markedly different than local ULIRGs and high-z submm galaxies, which have more excited and compact gas.

preprint2010arXiv

Weak 13CO in the Cloverleaf Quasar: evidence for a young, early generation starburst

Observations of 12CO at high redshift indicate rapid metal enrichment in the nuclear regions of at least some galaxies in the early universe. However, the enrichment may be limited to nuclei that are synthesized by short-lived massive stars, excluding classical secondary nuclei like 13CO. Testing this idea, we tentatively detect the 13CO J=3-2 line at a level of 0.3 Jy km/s toward the Cloverleaf Quasar at redshift 2.5. This is the first observational evidence for 13CO at high redshift. The 12CO/13CO J=3-2 luminosity ratio is with at least 40 much higher than ratios observed in molecular clouds of the Milky Way and in the ultraluminous galaxy Arp 220, but may be similar to that observed toward NGC 6240. Large Velocity Gradient (LVG) models simulating seven 12CO transitions and the 13CO line yield 12CO/13CO abundance ratios in excess of 100 for the Cloverleaf. It is possible that the measured ratio is affected by a strong submillimeter radiation field, which reduces the contrast between the 13CO line and the background. It is more likely, however, that the ratio is caused by a real deficiency of 13CO. A potential conflict with optical data, indicating high abundances also for secondary nuclei in quasars of high redshift, may be settled if the bulk of the CO emission is originating sufficiently far from the active galactic nucleus.

preprint2009arXiv

A CO emission line from the optical and near-IR undetected submillimeter galaxy GN10

We report the detection of a CO emission line from the submillimiter galaxy (SMG) GN10 in the GOODS-N field. GN10 lacks any counterpart in extremely deep optical and near-IR imaging obtained with the Hubble Space Telescope and ground-based facilities. This is a prototypical case of a source that is extremely obscured by dust, for which it is practically impossible to derive a spectroscopic redshift in the optical/near-IR. Under the hypothesis that GN10 is part of a proto-cluster structure previously identified at z~4.05 in the same field, we searched for CO[4-3] at 91.4 GHz with the IRAM Plateau de Bure Interferometer, and successfully detected a line. We find that the most likely redshift identification is z=4.0424+-0.0013, based on: 1) the very low chance that the CO line is actually serendipitous from a different redshift; 2) a radio-IR photometric redshift analysis; 3) the identical radio-IR SED, within a scaling factor, of two other SMGs at the same redshift. The faintness at optical/near-IR wavelengths requires an attenuation of A_V~5-7.5 mag. This result supports the case that a substantial population of very high-z SMGs exists that had been missed by previous spectroscopic surveys. This is the first time that a CO emission line has been detected for a galaxy that is invisible in the optical and near-IR. Our work demonstrates the power of existing and planned facilities for completing the census of star formation and stellar mass in the distant Universe by measuring redshifts of the most obscured galaxies through millimeter spectroscopy.

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

First redshift determination of an optically/UV faint submillimeter galaxy using CO emission lines

We report the redshift of a distant, highly obscured submm galaxy (SMG), based entirely on the detection of its CO line emission. We have used the newly commissioned Eight-MIxer Receiver (EMIR) at the IRAM 30m telescope, with its 8 GHz of instantaneous dual-polarization bandwidth, to search the 3-mm atmospheric window for CO emission from SMMJ14009+0252, a bright SMG detected in the SCUBA Lens Survey. A detection of the CO(3--2) line in the 3-mm window was confirmed via observations of CO(5--4) in the 2-mm window. Both lines constrain the redshift of SMMJ14009+0252 to z=2.9344, with high precision (dz=2 10^{-4}). Such observations will become routine in determining redshifts in the era of the Atacama Large Millimeter/submillimeter Array (ALMA).

preprint2009arXiv

HERACLES: The HERA CO-Line Extragalactic Survey

We present the HERA CO-Line Extragalactic Survey (HERACLES), an atlas of CO emission from 18 nearby galaxies that are also part of The HI Nearby Galaxy Survey (THINGS) and the Spitzer Infrared Nearby Galaxies Survey (SINGS). We used the HERA multi-pixel receiver on the IRAM 30-m telescope to map the CO J=2-1 line over the full optical disk (defined by the isophotal radius r_25) of each target, at 13" angular resolution and 2.6 km/s velocity resolution. Here we describe the observations and reduction of the data and show channel maps, azimuthally averaged profiles, integrated intensity maps, and peak intensity maps. The implied H2 masses range from 7 \times 10^6 to 6 \times 10^9 M_sun, with four low metallicity dwarf irregular galaxies yielding only upper limits. In the cases where CO is detected, the integrated H2-to-HI ratios range from 0.02 - 1.13 and H2-to-stellar mass ratios from 0.01 to 0.25. Exponential scale lengths of the CO emission for our targets are in the range 0.8 - 3.2 kpc, or 0.2 \pm 0.05 r_25. The intensity-weighted mean velocity of CO matches that of HI very well, with a 1σscatter of only 6 km/s. The CO J=2-1 / J=1-0 line ratio varies over a range similar to that found in the Milky Way and other nearby galaxies, \sim 0.6 - 1.0, with higher values found in the centers of galaxies. The typical line ratio, \sim 0.8, could be produced by optically thick gas with an excitation temperature of \sim 10 K.

preprint2009arXiv

Low, Milky-Way like, Molecular Gas Excitation of Massive Disk Galaxies at z~1.5

We present evidence for Milky-Way-like, low-excitation molecular gas reservoirs in near-IR selected massive galaxies at z~1.5, based on IRAM Plateau de Bure Interferometer CO[3-2] and NRAO Very Large Array CO[1-0] line observations for two galaxies that had been previously detected in CO[2-1] emission. The CO[3-2] flux of BzK-21000 at z=1.522 is comparable within the errors to its CO[2-1] flux, implying that the CO[3-2] transition is significantly sub-thermally excited. The combined CO[1-0] observations of the two sources result in a detection at the 3 sigma level that is consistent with a higher CO[1-0] luminosity than that of CO[2-1]. Contrary to what is observed in submillimeter galaxies and QSOs, in which the CO transitions are thermally excited up to J>=3, these galaxies have low-excitation molecular gas, similar to that in the Milky Way and local spirals. This is the first time that such conditions have been observed at high redshift. A Large Velocity Gradient analysis suggests that molecular clouds with density and kinetic temperature comparable to local spirals can reproduce our observations. The similarity in the CO excitation properties suggests that a high, Milky-Way-like, CO to H_2 conversion factor could be appropriate for these systems. If such low-excitation properties are representative of ordinary galaxies at high redshift, centimeter telescopes such as the Expanded Very Large Array and the longest wavelength Atacama Large Millimeter Array bands will be the best tools for studying the molecular gas content in these systems through the observations of CO emission lines.

preprint2009arXiv

NIR spectroscopy of SDSS J0303-0019: a low luminosity, high Eddington ratio quasar at z~6

We present sensitive near--infrared VLT ISAAC spectroscopic observations of the z=6.08 quasar SDSS J030331.40-001912.9. This QSO is more than a magnitude fainter than other QSOs at z~6 for which NIR spectroscopy has been obtained to date and is therefore presumably more representative of the QSO population at the end of Cosmic Reionization. Combining rest--frame UV continuum luminosity with the width measurements of the Mg II and C IV lines, we derive a black hole mass of 2(+1.0/-0.5) x 10^8 solar masses, the lowest mass observed for z~6 QSOs to date, and derive an Eddington ratio of 1.6(+0.4/-0.6), amongst the highest value derived for QSOs at any redshift. The Spitzer 24 micron non--detection of this QSO does not leave space for a significant hot dust component in its optical/near--infrared SED, in common with one other faint QSO at z=6, but in contrast to more than twenty more z=6 QSOs and all known lower redshift QSOs with sufficiently deep multi-wavelength photometry. We conclude that we have found evidence for differences in the intrinsic properties of at least one z~6 QSO as compared to the lower--redshift population.

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.

preprint2009arXiv

The Spitzer Local Volume Legacy: Survey Description and Infrared Photometry

The survey description and the near-, mid-, and far-infrared flux properties are presented for the 258 galaxies in the Local Volume Legacy (LVL). LVL is a Spitzer Space Telescope legacy program that surveys the local universe out to 11 Mpc, built upon a foundation of ultraviolet, H-alpha, and HST imaging from 11HUGS (11 Mpc H-alpha and Ultraviolet Galaxy Survey) and ANGST (ACS Nearby Galaxy Survey Treasury). LVL covers an unbiased, representative, and statistically robust sample of nearby star-forming galaxies, exploiting the highest extragalactic spatial resolution achievable with Spitzer. As a result of its approximately volume-limited nature, LVL augments previous Spitzer observations of present-day galaxies with improved sampling of the low-luminosity galaxy population. The collection of LVL galaxies shows a large spread in mid-infrared colors, likely due to the conspicuous deficiency of 8um PAH emission from low-metallicity, low-luminosity galaxies. Conversely, the far-infrared emission tightly tracks the total infrared emission, with a dispersion in their flux ratio of only 0.1 dex. In terms of the relation between infrared-to-ultraviolet ratio and ultraviolet spectral slope, the LVL sample shows redder colors and/or lower infrared-to-ultraviolet ratios than starburst galaxies, suggesting that reprocessing by dust is less important in the lower mass systems that dominate the LVL sample. Comparisons with theoretical models suggest that the amplitude of deviations from the relation found for starburst galaxies correlates with the age of the stellar populations that dominate the ultraviolet/optical luminosities.

preprint2008arXiv

THINGS: The HI Nearby Galaxy Survey

We present "The HI Nearby Galaxy Survey (THINGS)", a high spectral (<=5.2 km/s) and spatial (~6") resolution survey of HI emission in 34 nearby galaxies obtained using the NRAO Very Large Array (VLA). The overarching scientific goal of THINGS is to investigate fundamental characteristics of the interstellar medium (ISM) related to galaxy morphology, star formation and mass distribution across the Hubble sequence. Unique characteristics of the THINGS database are the homogeneous sensitivity as well as spatial and velocity resolution of the HI data which is at the limit of what can be achieved with the VLA for a significant number of galaxies. A sample of 34 objects at distances 2<D<15Mpc (resulting in linear resolutions of ~100 to 500pc) are targeted in THINGS, covering a wide range of star formation rates (10^-3 to 6 M_sun/yr), total HI masses M_HI (0.01 to 14x10^9 M_sun), absolute luminosities M_B (-11.5 to -21.7 mag) and metallicities (7.5 to 9.2 in units of 12+log[O/H]). We describe the setup of the VLA observations, the data reduction procedures and the creation of the final THINGS data products. We present an atlas of the integrated HI maps, the velocity fields, the second moment (velocity dispersion) maps and individual channel maps of each THINGS galaxy. The THINGS data products are made publicly available through a dedicated webpage. Accompanying THINGS papers address issues such as the small-scale structure of the ISM, the (dark) matter distribution in THINGS galaxies, and the processes leading to star formation.

preprint2007arXiv

Dust and Atomic Gas in Dwarf Irregular Galaxies of the M81 Group: The SINGS and THINGS view

[abridged] We present observations of the dust and atomic gas phase in seven dwarf irregular galaxies of the M81 group from the SINGS and THINGS surveys. The Spitzer observations provide a first glimpse of the nature of the non-atomic ISM in these metal-poor (Z~0.1 Z_sun), quiescent (SFR~0.001-0.1 M_sun/yr) dwarf galaxies. Dust emission is detected in five out of the seven targets. Most detected dust emission is restricted to HI column densities >1x10^21 cm^-2. Spitzer spectroscopy of two regions in the brightest galaxies (IC 2574 and Holmberg II) show distinctly different spectral shapes. The spectrum of IC 2574 shows aromatic features that are less luminous (relative to the FIR luminosity) compared to an average SINGS spiral galaxy by a factor of \~7 . The aromatic features in Holmberg~II (which has only a slightly lower gas-phase metallicity) are fainter than in IC 2574 by an order of magnitude. This result emphazises that the strength of the aromatic features is not a simple linear function of metallicity. We estimate dust masses of ~10^4-10^6 M_sun for the M81 dwarf galaxies, resulting in an average dust--to--gas ratio (M_dust/M_HI) of ~3x10^-4 (1.5x10^-3 if only the HI that is associated with dust emission is considered); this is an order of magnitude lower than the typical value derived for the SINGS spirals. The dwarf galaxies are underluminous per unit star formation rate at 70um as compared to the more massive galaxies in SINGS by a factor of ~2. However, the average 70um/160um ratio in the sample dwarf galaxies is higher than what is found in the other galaxies of the SINGS sample. This can be explained by a combination of a lower dust content in conjunction with a higher dust temperature in the dwarfs.

preprint2007arXiv

The Calibration of Mid-Infrared Star Formation Rate Indicators

With the goal of investigating the degree to which the mid-infrared emission traces the star formation rate (SFR), we analyze Spitzer 8 um and 24 um data of star-forming regions in a sample of 33 nearby galaxies with available HST/NICMOS images in the Paschen-alpha (1.8756 um) emission line. The galaxies are drawn from the Spitzer Infrared Nearby Galaxies Survey (SINGS) sample, and cover a range of morphologies and a factor ~10 in oxygen abundance. Published data on local low-metallicity starburst galaxies and Luminous Infrared Galaxies are also included in the analysis. Both the stellar-continuum-subtracted 8 um emission and the 24 um emission correlate with the extinction-corrected Pa-alpha line emission, although neither relationship is linear. Simple models of stellar populations and dust extinction and emission are able to reproduce the observed non-linear trend of the 24 um emission versus number of ionizing photons, including the modest deficiency of 24 um emission in the low metallicity regions, which results from a combination of decreasing dust opacity and dust temperature at low luminosities. Conversely, the trend of the 8 um emission as a function of the number of ionizing photons is not well reproduced by the same models. The 8 um emission is contributed, in larger measure than the 24 um emission, by dust heated by non-ionizing stellar populations, in agreement with previous findings. Two SFR calibrations, one using the 24 um emission and the other using a combination of the 24 um and H-alpha luminosities (Kennicutt et al. 2007), are presented. No calibration is presented for the 8 um emission, because of its significant dependence on both metallicity and environment. The calibrations presented here should be directly applicable to systems dominated by on-going star formation.

preprint2007arXiv

Variations in 24 micron morphologies among galaxies in the Spitzer Infrared Nearby Galaxies Survey: New insights into the Hubble sequence

To study the distribution of star formation and dust emission within nearby galaxies, we measured five morphological parameters in the 3.6 and 24 micron wave bands for 65 galaxies in the Spitzer Infrared Nearby Galaxies Survey (SINGS) and 8 galaxies that were serendipitously observed by SINGS. The morphological parameters demonstrate strong variations along the Hubble sequence, including statistically significant differences between S0/a-Sab and Sc-Sd galaxies. Early-type galaxies are generally found to be compact, centralized, symmetric sources in the 24 micron band, while late-type galaxies are generally found to be extended, asymmetric sources. These results suggest that the processes that increase the real or apparent sizes of galaxies' bulges also lead to more centralized 24 micron dust emission. Several phenomena, such as strong nuclear star formation, Seyfert activity, or outer ring structures, may cause galaxies to deviate from the general morphological trends observed at 24 microns. We also note that the 24 micron morphologies of Sdm-Im galaxies are quite varied, with some objects appearing very compact and symmetric while others appear diffuse and asymmetric. These variations reflect the wide variation in star formation in irregular galaxies as observed at other wavelengths. The variations in the 24 micron morphological parameters across the Hubble sequence mirror many of the morphological trends seen in other tracers of the ISM and in stellar emission. However, the 24 micron morphological parameters for the galaxies in this sample do not match the morphological parameters measured in the stellar wave bands. This implies that the distribution of dust emission is related to but not equivalent to the distribution of stellar emission.

preprint2007arXiv

Warm molecular hydrogen in the Spitzer SINGS galaxy sample

(simplified) Results on the properties of warm H2 in 57 normal galaxies are derived from H2 rotational transitions, obtained as part of SINGS. This study extends previous extragalactic surveys of H2, the most abundant constituent of the molecular ISM, to more common systems (L_FIR = e7 to 6e10 L_sun) of all morphological and nuclear types. The S(1) transition is securely detected in the nuclear regions of 86% of SINGS galaxies with stellar masses above 10^9.5 M_sun. The derived column densities of warm H2 (T > ~100 K), even though averaged over kiloparsec-scale areas, are commensurate with those of resolved PDRs; the median of the sample is 3e20 cm-2. They amount to between 1% and >30% of the total H2. The power emitted in the sum of the S(0) to S(2) transitions is on average 30% of the [SiII] line power, and ~4e-4 of the total infrared power (TIR) within the same area for star-forming galaxies, which is consistent with excitation in PDRs. The fact that H2 emission scales tightly with PAH emission, even though the average radiation field intensity varies by a factor ten, can also be understood if both tracers originate predominantly in PDRs, either dense or diffuse. A large fraction of the 25 LINER/Sy targets, however, strongly depart from the rest of the sample, in having warmer H2 in the excited states, and an excess of H2 emission with respect to PAHs, TIR and [SiII]. We propose a threshold in H2 to PAH power ratios, allowing the identification of low-luminosity AGNs by an excess H2 excitation. A dominant contribution from shock heating is favored in these objects. Finally, we detect, in nearly half the star-forming targets, non-equilibrium ortho to para ratios, consistent with FUV pumping combined with incomplete ortho-para thermalization by collisions, or possibly non-equilibrium PDR fronts advancing into cold gas.

preprint2006arXiv

Spitzer and JCMT Observations of the Active Galactic Nucleus in the Sombrero Galaxy (NGC 4594)

We present Spitzer 3.6-160 micron images, Spitzer mid-infrared spectra, and JCMT SCUBA 850 micron images of the Sombrero Galaxy (NGC 4594), an Sa galaxy with a 10^9 M_solar low luminosity active galactic nucleus (AGN). The brightest infrared sources in the galaxy are the nucleus and the dust ring. The spectral energy distribution of the AGN demonstrates that, while the environment around the AGN is a prominent source of mid-infrared emission, it is a relatively weak source of far-infrared emission, as had been inferred for AGN in previous research. The weak nuclear 160 micron emission and the negligible polycyclic aromatic hydrocarbon emission from the nucleus also implies that the nucleus is a site of only weak star formation activity and the nucleus contains relatively little cool interstellar gas needed to fuel such activity. We propose that this galaxy may be representative of a subset of low ionization nuclear emission region galaxies that are in a quiescent AGN phase because of the lack of gas needed to fuel circumnuclear star formation and Seyfert-like AGN activity. Surprisingly, the AGN is the predominant source of 850 micron emission. We examine the possible emission mechanisms that could give rise to the 850 micron emission and find that neither thermal dust emission, CO line emission, bremsstrahlung emission, nor the synchrotron emission observed at radio wavelengths can adequately explain the measured 850 micron flux density by themselves. The remaining possibilities for the source of the 850 micron emission include a combination of known emission mechanisms, synchrotron emission that is self-absorbed at wavelengths longer than 850 microns, or unidentified spectral lines in the 850 micron band.

preprint2006arXiv

The Effect of Star Formation on the Far-Infrared--Radio Correlation within Galaxies

Using data obtained for twelve galaxies as part of the {\it Spitzer} Infrared Nearby Galaxies Survey (SINGS) and the Westerbork Synthesis Radio Telescope (WSRT)-SINGS radio continuum survey, we study how star formation activity affects the far-infrared (FIR)--radio correlation {\it within} galaxies by testing a phenomenological model, which describes the radio image as a smeared version of the FIR image. The physical basis of this description is that cosmic-ray (CR) electrons will diffuse measurably farther than the mean free path of dust-heating photons before decaying by synchrotron radiation. This description works well in general. Galaxies with higher infrared surface brightnesses have best-fit smoothing scale-lengths of a few hundred parsecs, substantially shorter than those for lower surface brightness galaxies. We interpret this result to suggest that galaxies with higher disk averaged star formation rates have had a recent episode of enhanced star formation and are characterized by a higher fraction of young CR electrons that have traveled only a few hundred parsecs from their acceleration sites in supernova remnants compared to galaxies with lower star formation activity.

preprint2006arXiv

The Mid-Infrared Spectrum of Star-Forming Galaxies: Global Properties of PAH Emission

We present a sample of low-resolution 5-38um Spitzer IRS spectra of the inner few square kiloparsecs of 59 nearby galaxies spanning a large range of star formation properties. A robust method for decomposing mid-infrared galaxy spectra is described, and used to explore the behavior of PAH emission and the prevalence of silicate dust extinction. Evidence for silicate extinction is found in ~1/8 of the sample, at strengths which indicate most normal galaxies undergo A_V < ~3 magnitudes averaged over their centers. The contribution of PAH emission to the total infrared power is found to peak near 10% and extend up to ~20%, and is suppressed at metallicities Z < ~Z_sun/4, as well as in low-luminosity AGN environments. Strong inter-band PAH feature strength variations (2-5x) are observed, with the presence of a weak AGN and, to a lesser degree, increasing metallicity shifting power to the longer wavelength bands. A peculiar PAH emission spectrum with markedly diminished 5-8um features arises among the sample solely in systems with relatively hard radiation fields harboring low-luminosity AGN. The AGN may modify the emitting grain distribution and provide the direct excitation source of the unusual PAH emission, which cautions against using absolute PAH strength to estimate star formation rates in systems harboring active nuclei. Alternatively, the low star formation intensity often associated with weak AGN may affect the spectrum. The effect of variations in the mid-infrared spectrum on broadband infrared surveys is modeled, and points to more than a factor of two uncertainty in results which assume a fixed PAH emission spectrum, for redshifts z=0-2.5.

preprint2005arXiv

VLA Imaging of the Intriguing HI Cloud HIJASS J1021+6842 in the M81 Group

We present VLA HI 21cm observations of HIJASS J1021+6842 which has been discovered in the direction of the M81 group. Our synthesis imaging reveals that the HI is distributed over a larger angular extent and velocity range than the single dish discovery observations. Assuming that HIJASS J1021+6842 is at the distance of the M81 group, we detect 1.5 x 10^8 M_sun of HI distributed over as much as 30 kpc, i.e., substantially larger than the biggest dwarf galaxies in the same group. At the depth of our imaging, the HI appears to be confined to at least 7 clouds. Peak HI column densities are ~1.8 x 10^20 atoms cm^-2 which is well below the canonical star formation threshold of ~10^21 atoms cm^-2 and therefore consistent with the fact that no optical counterpart has as yet been identified. A gradient in velocity is observed across the extent of the detected HI; assuming that the object is gravitationally bound we derive a dynamical mass of 7 x 10^9 M_sun and a dark-to-luminous mass ratio of >10. Alternatively, a tidal origin may also result in the observed velocity gradient which would lead to a considerably lower dynamical mass. Given the above properties and the absence of evidence of a stellar population, HIJASS J1021+6842 is unique amongst the other systems in the M81 group.

preprint2004arXiv

First Results from THINGS: The HI Nearby Galaxy Survey

We describe The HI Nearby Galaxy Survey (THINGS), the largest program ever undertaken at the VLA to perform 21-cm HI observations of the highest quality (~7'', <= 5 km/s resolution) of nearby galaxies. The goal of THINGS is to investigate key characteristics related to galaxy morphology, star formation and mass distribution across the Hubble sequence. A sample of 34 objects with distances between 3 and 10 Mpc will be observed, covering a wide range of evolutionary stages and properties. Data from THINGS will complement SINGS, the Spitzer Infrared Nearby Galaxy Survey. For the THINGS sample, high-quality observations at comparable resolution will thus be available from the X-ray regime through to the radio part of the spectrum. THINGS data can be used to investigate issues such as the small-scale structure of the ISM, its three-dimensional structure, the (dark) matter distribution and processes leading to star formation. To demonstrate the quality of the THINGS data products, we present some prelimary HI maps here of four galaxies from the THINGS sample.

preprint2004arXiv

Mid-Infrared IRS Spectroscopy of NGC 7331: A First Look at the SINGS Legacy

The nearby spiral galaxy NGC 7331 was spectrally mapped from 5-38um using all modules of Spitzer's IRS spectrograph. A strong new dust emission feature, presumed due to PAHs, was discovered at 17.1um. The feature's intensity is nearly half that of the ubiquitous 11.3um band. The 7-14um spectral maps revealed significant variation in the 7.7 and 11.3um PAH features between the stellar ring and nucleus. Weak [OIV] 25.9um line emission was found to be centrally concentrated in the nucleus, with an observed strength over 10% of the combined neon line flux, indicating an AGN or unusually active massive star photo-ionization. Two [SIII] lines fix the characteristic electron density in the HII regions at n_e < ~200 cm^-3. Three detected H_2 rotational lines, tracing warm molecular gas, together with the observed IR continuum, are difficult to match with standard PDR models. Either additional PDR heating or shocks are required to simultaneously match lines and continuum.

preprint2004arXiv

Spatially Resolved Ultraviolet, H-alpha, Infrared, and Radio Star Formation in M81

We present Multiband Imaging Photometer for Spitzer (MIPS) observations of M81 at 24, 70, and 160 micron. The grand design nature of M81 is clearly seen, showing two well resolved spiral arms containing numerous bright star forming regions. The MIPS images reveal a significant amount of cold dust associated with the spiral arms. We investigate the variation of the ultraviolet (UV), H-alpha, and infrared (IR) luminosities and star formation rate (SFR) indicators across the face of M81 using the MIPS images and archival UV and H-alpha images. For regions in M81, we find that UV and H-alpha SFRs (uncorrected for dust attenuation) are always lower than the IR SFR. The cause of this behavior is dust attenuation and/or using SFR calibrations appropriate for entire galaxies, not regions in galaxies. The characteristics of the dust attenuation for the regions indicate the dust grains and/or geometry are different from those in starburst galaxies. The behavior of the infrared-radio correlation in M81 is seen to vary from the global average, with variations correlated with the morphology of M81.

preprint2003arXiv

SINGS: The SIRTF Nearby Galaxies Survey

The SIRTF Nearby Galaxy Survey is a comprehensive infrared imaging and spectroscopic survey of 75 nearby galaxies. Its primary goal is to characterize the infrared emission of galaxies and their principal infrared-emitting components, across a broad range of galaxy properties and star formation environments. SINGS will provide new insights into the physical processes connecting star formation to the interstellar medium properties of galaxies, and provide a vital foundation for understanding infrared observations of the distant universe and ultraluminous and active galaxies. The galaxy sample and observing strategy have been designed to maximize the scientific and archival value of the data set for the SIRTF user community at large. The SIRTF images and spectra will be supplemented by a comprehensive multi-wavelength library of ancillary and complementary observations, including radio continuum, HI, CO, submillimeter, BVRIJHK, H-alpha, Paschen-alpha, ultraviolet, and X-ray data. This paper describes the main astrophysical issues to be addressed by SINGS, the galaxy sample and the observing strategy, and the SIRTF and other ancillary data products.

preprint2002arXiv

The Broadband Afterglow of GRB980329

We present radio observations of the afterglow of the bright gamma-ray burst GRB980329 made between one month and several years after the burst, a re-analysis of previously published submillimeter data, and late-time optical and near-infrared (NIR) observations of the host galaxy. From the absence of a spectral break in the optical/NIR colors of the host galaxy, we exclude the earlier suggestion that GRB980329 lies at a redshift of z >~5. We combine our data with the numerous multi-wavelength observations of the early afterglow, fit a comprehensive afterglow model to the entire broadband dataset, and derive fundamental physical parameters of the blast-wave and its host environment. Models for which the ejecta expand isotropically require both a high circumburst density and extreme radiative losses from the shock. No low density model (n << 10 cm^{-3}) fits the data. A burst with a total energy of ~ 10^{51} erg, with the ejecta narrowly collimated to an opening angle of a few degrees, driven into a surrounding medium with density ~ 20 cm^{-3}, provides a satisfactory fit to the lightcurves over a range of redshifts.