Source author record

E. Schinnerer

E. Schinnerer 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

121works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

121 published item(s)

preprint2022arXiv

Linking stellar populations to HII regions across nearby galaxies: I. Constraining pre-supernova feedback from young clusters in NGC1672

One of the fundamental factors regulating the evolution of galaxies is stellar feedback. However, we still do not have strong observational constraints on the relative importance of the different feedback mechanisms (e.g. radiation, ionised gas pressure, stellar winds) in driving HII region evolution and molecular cloud disruption. In this letter, we constrain the relative importance of the various feedback mechanisms from young massive star populations by resolving HII regions across the disk of the nearby star-forming galaxy NGC 1672. We combine measurements of ionised gas nebular lines obtained by PHANGS-MUSE, with high-resolution imaging from the HST in both the narrow-band Hα and broad-band filters. We identify a sample of 40 isolated, compact HII regions in the HST Hα image, for which we measure the sizes that were previously unresolved in seeing-limited ground-based observations. Additionally, we identify the ionisation source(s) for each HII region from catalogues produced as part of the PHANGS-HST survey. We find that the HII regions investigated are mildly dominated by thermal or wind pressure, yet their elevation above the radiation pressure is within the expected uncertainty range. We see that radiation pressure provides a substantially higher contribution to the total pressure than previously found in the literature over similar size scales. In general, we find higher pressures within more compact HII regions, which is driven by the inherent size scaling relations of each pressure term, albeit with significant scatter introduced by the variation in the stellar population properties (e.g. luminosity, mass, age, metallicity). For nearby galaxies, here we provide a promising approach that could yield the statistics required to map out how the importance of different stellar feedback mechanisms evolve over the lifetime of an HII region.

preprint2022arXiv

Probing star formation and ISM properties using galaxy disk inclination III: Evolution in dust opacity and clumpiness between redshift 0.0 < z < 0.7 constrained from UV to NIR

(Abridged) In this paper, we use the Tuffs et al. attenuation - inclination models in ultraviolet (UV), optical, and near-infrared (NIR) bands to investigate the average global dust properties in galaxies as a function of stellar mass $M_{*}$, stellar mass surface density $μ_{*}$, star-formation rate $SFR$, specific star-formation rate $sSFR$, star-formation main-sequence offset $dMS$, and star-formation rate surface density $Σ_{SFR}$ at redshifts $z \sim 0$ and $z \sim 0.7$. We use star-forming galaxies from SDSS ($\sim$ 20000) and GAMA ($\sim$ 2000) to form our low-z sample at $0.04 < z < 0.1$ and star-forming galaxies from COSMOS ($\sim$ 2000) for the sample at $0.6 <z < 0.8$. We find that galaxies at $z \sim 0.7$ have higher optical depth $τ_{B}^{f}$ and clumpiness $F$ than galaxies at $z \sim 0$. The increase in $F$ hints that the stars of $z \sim 0.7$ galaxies are less likely to escape their birth cloud, which might indicate that the birth clouds are larger. We also found that $τ_{B}^{f}$ increases with $M_{*}$ and $μ_{*}$independent of sample and therefore redshift. We found no clear trends in $τ_{B}^{f}$ or $F$ with $SFR$, which could imply that the dust mass distribution is independent of $SFR$. In turn, this would imply that the balance of dust formation and destruction is independent of the $SFR$. Based on an analysis of the inclination-dependence of the Balmer decrement, we find that reproducing the Balmer line emission requires not only a completely optically thick dust component associated with star forming regions, as in the standard Tuffs et al. model, but an extra component of optically thin dust within the birth clouds. This new component implies the existence of dust inside HII regions that attenuates the Balmer emission before it escapes through gaps in the birth cloud and we find it is more important in high-mass galaxies.

preprint2022arXiv

Variations in the $Σ_{\rm SFR} {-} Σ_{\rm mol} {-} Σ_{\rm \star}$ plane across galactic environments in PHANGS galaxies

There exists some consensus that stellar mass surface density ($Σ_{*}$) and molecular gas mass surface density ($Σ_{\rm mol}$) are the main quantities responsible for locally setting the star formation rate. This regulation is inferred from locally resolved scaling relations between these two quantities and the star formation rate surface density ($Σ_{\rm SFR}$). However, the universality of these relations is debated. Here, we probe the interplay between these three quantities across different galactic environments at a spatial resolution of 150 pc. We perform a hierarchical Bayesian linear regression to find the best set of parameters $C_{*}$, $C_{\rm mol}$, and $C_{\rm norm}$ that describe the star-forming plane conformed by these quantities, such that $\log Σ_{\rm SFR} = C_{*} \log Σ_{*} + C_{\rm mol} \log Σ_{\rm mol} + C_{\rm norm}$, and explore variations in the determined parameters across galactic environments, focusing our analysis on the $C_{*}$ and $C_{\rm mol}$ slopes. We find signs of variations in the posterior distributions of $C_{*}$ and $C_{\rm mol}$ across different galactic environments. Bars show the most negative value of $C_{*}$, a sign of longer depletion times, while spiral arms show the highest $C_{*}$ among all environments. We conclude that systematic variations in the interplay of $Σ_{*}$, $Σ_{\rm mol}$ and $Σ_{\rm SFR}$ across galactic environments exist at a spatial resolution of 150 pc, and we interpret these variations as produced by an additional mechanism regulating the formation of stars that is not captured by either $Σ_{*}$ or $Σ_{\rm mol}$. We find that these variations correlate with changes in the star formation efficiency across environments, which could be linked to the dynamical state of the gas that prevents it from collapsing and forming stars, or to changes in the molecular gas fraction.

preprint2020arXiv

ALMA Resolves Giant Molecular Clouds in a Tidal Dwarf Galaxy

Tidal dwarf galaxies (TDGs) are gravitationally bound condensations of gas and stars formed during galaxy interactions. Here we present multi-configuration ALMA observations of J1023+1952, a TDG in the interacting system Arp 94, where we resolve CO(2-1) emission down to giant molecular clouds (GMCs) at 0.64" ~ 45pc resolution. We find a remarkably high fraction of extended molecular emission (~80-90%), which is filtered out by the interferometer and likely traces diffuse gas. We detect 111 GMCs that give a similar mass spectrum as those in the Milky Way and other nearby galaxies (a truncated power law with slope of -1.76+/-0.13). We also study Larson's laws over the available dynamic range of GMC properties (~2 dex in mass and ~1 dex in size): GMCs follow the size-mass relation of the Milky Way, but their velocity dispersion is higher such that the size-linewidth and virial relations appear super-linear, deviating from the canonical values. The global molecular-to-atomic gas ratio is very high (~1) while the CO(2-1)/CO(1-0) ratio is quite low (~0.5), and both quantities vary from north to south. Star formation is predominantly taking place in the south of the TDG, where we observe projected offsets between GMCs and young stellar clusters ranging from ~50pc to ~200pc; the largest offsets correspond to the oldest knots, as seen in other galaxies. In the quiescent north, we find more molecular clouds and a higher molecular-to-atomic gas ratio (~1.5); atomic and diffuse molecular gas also have a higher velocity dispersion there. Overall, the organisation of the molecular ISM in this TDG is quite different from other types of galaxies on large scales, but the properties of GMCs seem fairly similar, pointing to near universality of the star-formation process on small scales.

preprint2020arXiv

Dense Gas in a Giant Molecular Filament

Recent surveys of the Galactic plane in the dust continuum and CO emission lines reveal that large ($\gtrsim 50$~pc) and massive ($\gtrsim 10^5$~$M_\odot$) filaments, know as giant molecular filaments (GMFs), may be linked to galactic dynamics and trace the mid-plane of the gravitational potential in the Milky Way. We have imaged one entire GMF located at $l\sim$52--54$^\circ$ longitude, GMF54 ($\sim$68~pc long), in the empirical dense gas tracers using the HCN(1--0), HNC(1--0), HCO$^+$(1--0) lines, and their $^{13}$C isotopologue transitions, as well as the N$_2$H$^+$(1--0) line. We study the dense gas distribution, the column density probability density functions (N-PDFs) and the line ratios within the GMF. The dense gas molecular transitions follow the extended structure of the filament with area filling factors between 0.06 and 0.28 with respect to $^{13}$CO(1--0). We constructed the N-PDFs of H$_2$ for each of the dense gas tracers based on their column densities and assumed uniform abundance. The N-PDFs of the dense gas tracers appear curved in log-log representation, and the HCO$^+$ N-PDF has the largest log-normal width and flattest power-law slope index. Studying the N-PDFs for sub-regions of GMF54, we found an evolutionary trend in the N-PDFs that high-mass star forming and Photon-Dominate Regions (PDRs) have flatter power-law indices. The integrated intensity ratios of the molecular lines in GMF54 are comparable to those in nearby galaxies. In particular, the N$_2$H$^+$/$^{13}$CO ratio, which traces the dense gas fraction, has similar values in GMF54 and all nearby galaxies except ULIRGs.

preprint2020arXiv

Extended H$α$ over compact far-infrared continuum in dusty submillimeter galaxies -- Insights into dust distributions and star-formation rates at $z\sim2$

Using data from ALMA and near-infrared (NIR) integral field spectrographs including both SINFONI and KMOS on the VLT, we investigate the two-dimensional distributions of H$α$ and rest-frame far-infrared (FIR) continuum in six submillimeter galaxies at $z\sim2$. At a similar spatial resolution ($\sim$0.5" FWHM; $\sim$4.5 kpc at $z=2$), we find that the half-light radius of H$α$ is significantly larger than that of the FIR continuum in half of the sample, and on average H$α$ is a median factor of $2.0\pm0.4$ larger. Having explored various ways to correct for the attenuation, we find that the attenuation-corrected H$α$-based SFRs are systematically lower than the IR-based SFRs by at least a median factor of $3\pm1$, which cannot be explained by the difference in half-light radius alone. In addition, we find that in 40% of cases the total $V$-band attenuation ($A_V$) derived from energy balance modeling of the full ultraviolet(UV)-to-FIR spectral energy distributions (SEDs) is significantly higher than that derived from SED modeling using only the UV-to-NIR part of the SEDs, and the discrepancy appears to increase with increasing total infrared luminosity. Finally, considering all our findings along with the studies in the literature, we postulate that the dust distributions in SMGs, and possibly also in less IR luminous $z\sim2$ massive star-forming galaxies, can be decomposed into three main components; the diffuse dust heated by older stellar populations, the more obscured and extended young star-forming HII regions, and the heavily obscured central regions that have a low filling factor but dominate the infrared luminosity in which the majority of attenuation cannot be probed via UV-to-NIR emissions.

preprint2020arXiv

Molecular gas and star formation activity in LIRGs in clusters at intermediate redshifts

We investigate the role of dense Mpc-scale environments in processing molecular gas of cluster galaxies as they fall into the cluster cores. We consider $\sim20$ luminous infrared galaxies (LIRGs) in intermediate-$z$ clusters, from the Hershel Lensing Survey and the Local Cluster Substructure Survey. They include MACS J0717.5+3745 at $z=0.546$ and Abell 697, 963, 1763, and 2219 at $z=0.2-0.3$. We have performed far infrared to ultraviolet spectral energy distribution modeling of the LIRGs, which span cluster-centric distances within $r/r_{200}\simeq0.2-1.6$. We have observed the LIRGs in CO(1$\rightarrow$0) or CO(2$\rightarrow$1) with the Plateau de Bure interferometer and its successor NOEMA, as part of five observational programs carried out between 2012 and 2017. We have compared the molecular gas to stellar mass ratio $M(H_2)/M_\star$, star formation rate (SFR), and depletion time ($τ_{\rm dep}$) of the LIRGs with those of a compilation of cluster and field star forming galaxies. The targeted LIRGs have SFR, $M(H_2)/M_\star$, and $τ_{\rm dep}$ that are consistent with those of both main sequence (MS) field galaxies and star forming galaxies from the comparison sample. However we find that the depletion time, normalized to the MS value, increases with increasing $r/r_{200}$, with a significance of $2.8σ$, which is ultimately due to a deficit of cluster core LIRGs with $τ_{\rm dep}\gtrsimτ_{\rm dep,MS}$. We suggest that a rapid exhaustion of the molecular gas reservoirs occurs in the cluster LIRGs and is effective in suppressing their star formation. This mechanism may explain the exponential decrease of the fraction of cluster LIRGs with cosmic time. The compression of the gas in LIRGs, possibly induced by intra-cluster medium shocks, may be responsible for the short depletion timescales, observed in a large fraction of cluster core LIRGs.

preprint2020arXiv

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

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

preprint2020arXiv

The Star Formation in Radio Survey: 3 - 33 GHz Imaging of Nearby Galaxy Nuclei and Extranuclear Star-forming Regions

We present 3, 15, and 33 GHz imaging towards galaxy nuclei and extranuclear star-forming regions using the Karl G. Jansky Very Large Array as part of the Star Formation in Radio Survey. With $3-33$ GHz radio spectra, we measured the spectral indices and corresponding thermal (free-free) emission fractions for a sample of 335 discrete regions having significant detections in at least two radio bands. After removing 14 likely background galaxies, we find that the median thermal fraction at 33 GHz is $92 \pm 0.8\%$ with a median absolute deviation of $11\%$, when a two-component power-law model is adopted to fit the radio spectrum. Limiting the sample to 238 sources that are confidently identified as star-forming regions, and not affected by potential AGN contamination (i.e., having galactocentric radii $r_{\rm G} \geq 250$ pc), results in a median thermal fraction of $93 \pm 0.8 \%$ with a median absolute deviation of $10\%$. We further measure the thermal fraction at 33 GHz for 163 regions identified at 7" resolution to be $94 \pm 0.8 \%$ with a median absolute deviation of $8\%$. Together, these results confirm that free-free emission dominates the radio spectra of star-forming regions on scales up to $\sim$500 pc in normal star-forming galaxies. We additionally find a factor of $\sim$1.6 increase in the scatter of the measured spectral index and thermal fraction distributions as a function of decreasing galactocentric radius. This trend is likely reflective of the continuous star-formation activity occurring in the galaxy centers, resulting a larger contribution of diffuse nonthermal emission relative to star-forming regions in the disk.

preprint2020arXiv

Ubiquitous velocity fluctuations throughout the molecular interstellar medium

The density structure of the interstellar medium (ISM) determines where stars form and release energy, momentum, and heavy elements, driving galaxy evolution. Density variations are seeded and amplified by gas motion, but the exact nature of this motion is unknown across spatial scale and galactic environment. Although dense star-forming gas likely emerges from a combination of instabilities, convergent flows, and turbulence, establishing the precise origin is challenging because it requires quantifying gas motion over many orders of magnitude in spatial scale. Here we measure the motion of molecular gas in the Milky Way and in nearby galaxy NGC 4321, assembling observations that span an unprecedented spatial dynamic range ($10^{-1}{-}10^3$ pc). We detect ubiquitous velocity fluctuations across all spatial scales and galactic environments. Statistical analysis of these fluctuations indicates how star-forming gas is assembled. We discover oscillatory gas flows with wavelengths ranging from $0.3{-}400$ pc. These flows are coupled to regularly-spaced density enhancements that likely form via gravitational instabilities. We also identify stochastic and scale-free velocity and density fluctuations, consistent with the structure generated in turbulent flows. Our results demonstrate that ISM structure cannot be considered in isolation. Instead, its formation and evolution is controlled by nested, interdependent flows of matter covering many orders of magnitude in spatial scale.

preprint2019arXiv

A$^3$COSMOS: The Dust Attenuation of Star-Forming Galaxies at $z=2.5-4.0$ from the COSMOS-ALMA Archive

We present an analysis of the dust attenuation of star forming galaxies at $z=2.5-4.0$ through the relationship between the UV spectral slope ($β$), stellar mass ($M_{\ast}$) and the infrared excess (IRX$=L_{\rm{IR}}/L_{\rm{UV}}$) based on far-infrared continuum observations from the Atacama Large Millimeter/sub-millimeter Array (ALMA). Our study exploits the full ALMA archive over the COSMOS field processed by the A$^3$COSMOS team, which includes an unprecedented sample of $\sim1500$ galaxies at $z\sim3$ as primary or secondary targets in ALMA band 6 or 7 observations with a median continuum sensitivity of 126 $\rm{μJy/beam}$ (1$σ$). The detection rate is highly mass dependent, decreasing drastically below $\log (M_{\ast}/M_{\odot})=10.5$. The detected galaxies show that the IRX-$β$ relationship of massive ($\log M_{\ast}/M_{\odot} > 10$) main sequence galaxies at $z=2.5-4.0$ is consistent with that of local galaxies, while starbursts are generally offset by $\sim0.5\,{\rm dex}$ to larger IRX values. At the low mass end, we derive upper limits on the infrared luminosities through stacking of the ALMA data. The combined IRX-$M_{\ast}$ relation at $\rm{log\,(M_{\ast}/M_{\odot})>9}$ exhibits a significantly steeper slope than reported in previous studies at similar redshifts, implying little dust obscuration at $\log M_{\ast}/M_{\odot}<10$. However, our results are consistent with early measurements at $z\sim5.5$, indicating a potential redshift evolution between $z\sim2$ and $z\sim6$. Deeper observations targeting low mass galaxies will be required to confirm this finding.

preprint2019arXiv

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

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

preprint2019arXiv

Discovery of four apparently cold dusty galaxies at z=3.62-5.85 in the COSMOS field: direct evidence of CMB impact on high-redshift galaxy observables

We report Atacama Large Millimetre Array (ALMA) observations of four high-redshift dusty star-forming galaxy candidates selected from far-Infrared (FIR)/submm observations in the COSMOS field. We securely detect all galaxies in the continuum and spectroscopically confirm them at z=3.62--5.85 using ALMA 3mm line scans, detecting multiple CO and/or [CI] transitions. This includes the most distant dusty galaxy currently known in the COSMOS field, ID85001929 at z=5.847. These redshifts are lower than we had expected as these galaxies have substantially colder dust temperatures (i.e., their SEDs peak at longer rest frame wavelengths) than most literature sources at z>4. The observed cold dust temperatures are best understood as evidence for optically thick dust continuum in the FIR, rather than the result of low star formation efficiency with rapid metal enrichment. We provide direct evidence that, given their cold spectral energy distributions, CMB plays a significant role biasing their observed Rayleigh-Jeans (RJ) slopes to unlikely steep values and, possibly, reducing their CO fluxes by a factor of two. We recover standard RJ slopes when the CMB contribution is taken into account. High resolution ALMA imaging shows compact morphology and evidence for mergers. This work reveals a population of cold dusty star-forming galaxies that were under-represented in current surveys, and are even colder than typical Main Sequence galaxies at the same redshift. High FIR dust optical depth might be a widespread feature of compact starbursts at any redshift.

preprint2019arXiv

Mapping metallicity variations across nearby galaxy disks

The distribution of metals within a galaxy traces the baryon cycle and the buildup of galactic disks, but the detailed gas phase metallicity distribution remains poorly sampled. We have determined the gas phase oxygen abundances for 7,138 HII regions across the disks of eight nearby galaxies using VLT/MUSE optical integral field spectroscopy as part of the PHANGS-MUSE survey. After removing the first order radial gradients present in each galaxy, we look at the statistics of the metallicity offset (Delta O/H) and explore azimuthal variations. Across each galaxy, we find low (sigma=0.03-0.05 dex) scatter at any given radius, indicative of efficient mixing. We compare physical parameters for those HII regions that are 1 sigma outliers towards both enhanced and reduced abundances. Regions with enhanced abundances have high ionization parameter, higher Halpha luminosity, lower Halpha velocity dispersion, younger star clusters and associated molecular gas clouds show higher molecular gas densities. This indicates recent star formation has locally enriched the material. Regions with reduced abundances show increased Halpha velocity dispersions, suggestive of mixing introducing more pristine material. We observe subtle azimuthal variations in half of the sample, but can not always cleanly associate this with the spiral pattern. Regions with enhanced and reduced abundances are found distributed throughout the disk, and in half of our galaxies we can identify subsections of spiral arms with clearly associated metallicity gradients. This suggests spiral arms play a role in organizing and mixing the ISM.

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

(Sub)millimetre interferometric imaging of a sample of COSMOS/AzTEC submillimetre galaxies III. Environments

We study the environment of 23 submillimetre galaxies (SMGs) drawn from the JCMT/AzTEC 1.1mm S/N-limited sample in the COSMOS field, as well as 4 COSMOS SMGs at z_spec>4.5, and 1 at z_spec=2.49, yielding a sample of 28 SMGs. We search for overdensities using the COSMOS photometric redshifts based on over 30 UV-NIR photometric bands, reaching an accuracy of sigma(Delta z/(1+z))=0.0067 (0.0155) at z<3.5 (>3.5). To identify overdensities we apply the Voronoi tessellation analysis, and estimate the overdensity estimator delta_g as a function of distance from the SMG and/or overdensity center. We test and validate our approach via simulations, X-ray detected groups, and spectroscopic verifications using VUDS and zCOSMOS catalogues showing that even with photometric redshifts in COSMOS we can efficiently retrieve overdensities out to z~5. Our results yield that 11/23 (48%) JCMT/AzTEC 1.1mm SMGs occupy overdense environments. Considering the entire JCMT/AzTEC 1.1mm S/N>4 sample, and accounting for the expected fraction of spurious detections, yields that 35-61% of the SMGs in the S/N-limited sample occupy overdense environments. We perform an X-ray stacking analysis in the 0.5-2keV band using a 32" aperture and our SMG positions, and find statistically significant detections. For our z<2 [z>2] subsample we find an average flux of (4.0+/-0.8)x10^{-16} [(1.3+/-0.5)x10^{-16}] erg/s/cm^2 and a corresponding total mass of M_200 = 2.8x10^{13} [2x10^{13}] MSol. Our results suggest a higher occurrence of SMGs occupying overdense environments at z>3, than at z<3. This may be understood if highly star forming galaxies can only be formed in the highest peaks of the density field tracing the most massive dark matter haloes at early cosmic epochs, while at later times cosmic structure may have matured sufficiently that more modest overdensities correspond to sufficiently massive haloes to form SMGs.

preprint2016arXiv

A novel look at energy equipartition in globular clusters

Two-body interactions play a major role in shaping the structural and dynamical properties of globular clusters (GCs) over their long-term evolution. In particular, GCs evolve toward a state of partial energy equipartition that induces a mass-dependence in their kinematics. By using a set of Monte Carlo cluster simulations evolved in quasi-isolation, we show that the stellar mass dependence of the velocity dispersion $σ(m)$ can be described by an exponential function $σ^2\propto \exp(-m/m_\mathrm{eq})$, with the parameter $m_\mathrm{eq}$ quantifying the degree of partial energy equipartition of the systems. This simple parametrization successfully captures the behaviour of the velocity dispersion at lower as well as higher stellar masses, that is, the regime where the system is expected to approach full equipartition. We find a tight correlation between the degree of equipartition reached by a GC and its dynamical state, indicating that clusters that are more than about 20 core relaxation times old, have reached a maximum degree of equipartition. This equipartition$-$dynamical state relation can be used as a tool to characterize the relaxation condition of a cluster with a kinematic measure of the $m_\mathrm{eq}$ parameter. Vice versa, the mass-dependence of the kinematics can be predicted knowing the relaxation time solely on the basis of photometric measurements. Moreover, any deviations from this tight relation could be used as a probe of a peculiar dynamical history of a cluster. Finally, our novel approach is important for the interpretation of state-of-the-art Hubble Space Telescope proper motion data, for which the mass dependence of kinematics can now be measured, and for the application of modeling techniques which take into consideration multi-mass components and mass segregation.

preprint2016arXiv

ALMA resolves the torus of NGC 1068: continuum and molecular line emission

We have used the Atacama Large Millimeter Array (ALMA) to map the emission of the CO(6-5) molecular line and the 432 μm continuum emission from the 300 pc-sized circumnuclear disk (CND) of the nearby Seyfert 2 galaxy NGC 1068 with a spatial resolution of ~4 pc. These observations spatially resolve the CND and, for the first time, image the dust emission, the molecular gas distribution, and the kinematics from a 7-10 pc-diameter disk that represents the submillimeter counterpart of the putative torus of NGC 1068. We fitted the nuclear spectral energy distribution of the torus using ALMA and near and mid-infrared (NIR/MIR) data with CLUMPY models. The mass and radius of the best-fit solution for the torus are both consistent with the values derived from the ALMA data alone: Mgas_torus=(1+-0.3)x10^5 Msun and Rtorus=3.5+-0.5 pc. The dynamics of the molecular gas in the torus show non-circular motions and enhanced turbulence superposed on the rotating pattern of the disk. The kinematic major axis of the CO torus is tilted relative to its morphological major axis. By contrast with the nearly edge-on orientation of the H2O megamaser disk, we have found evidence suggesting that the molecular torus is less inclined (i=34deg-66deg) at larger radii. The lopsided morphology and complex kinematics of the torus could be the signature of the Papaloizou-Pringle instability, long predicted to likely drive the dynamical evolution of active galactic nuclei (AGN) tori.

preprint2016arXiv

Characterizing spiral arm and interarm star formation

Interarm star formation contributes significantly to a galaxy's star formation budget, and provides an opportunity to study stellar birthplaces unperturbed by spiral arm dynamics. Using optical integral field spectroscopy of the nearby galaxy NGC 628 with VLT/MUSE, we construct Halpha maps including detailed corrections for dust extinction and stellar absorption to identify 391 HII regions at 35pc resolution over 12 kpc^2. Using tracers sensitive to the underlying gravitational potential, we associate HII regions with either arm (271) or interarm (120) environments. Using our full spectral coverage of each region, we find that most HII region physical properties (luminosity, size, metallicity, ionization parameter) are independent of environment. We calculate the fraction of Halpha luminosity due to the diffuse ionized gas (DIG) background contaminating each HII region, and find the DIG surface brightness to be higher within HII regions compared to the surroundings, and slightly higher within arm HII regions. Use of the temperature sensitive [SII]/Halpha line ratio map instead of the Halpha surface brightness to identify HII region boundaries does not change this result. Using the dust attenuation as a tracer of the gas, we find depletion times consistent with previous work (2 x 10^9 yr) with no differences between the arm and interarm, however this is very sensitive to the DIG correction. Unlike molecular clouds, which can be dynamically affected by the galactic environment, we see fairly consistent HII region properties in both arm and interarm environments. This suggests either a difference in arm star formation and feedback, or a decoupling of dense star forming clumps from the more extended surrounding molecular gas.

preprint2016arXiv

Gas Fraction and Depletion Time of Massive Star Forming Galaxies at z~3.2: No Change in Global Star Formation Process out to z>3

The observed evolution of the gas fraction and its associated depletion time in main sequence (MS) galaxies provides insights on how star formation proceeds over cosmic time. We report ALMA detections of the rest-frame $\sim$300$μ$m continuum observed at 240 GHz for 45 massive ($\rm \langle log(M_{\star}(M_{\odot}))\rangle=10.7$), normal star forming ($\rm \langle log(sSFR(yr^{-1}))\rangle=-8.6$), i.e. MS, galaxies at $\rm z\approx3.2$ in the COSMOS field. From an empirical calibration between cold neutral, i.e. molecular and atomic, gas mass $\rm M_{gas}$ and monochromatic (rest-frame) infrared luminosity, the gas mass for this sample is derived. Combined with stellar mass $\rm M_{\star}$ and star formation rate (SFR) estimates (from {\sc MagPhys} fits) we obtain a median gas fraction of $\rm μ_{gas}=M_{gas}/M_{\star}=1.65_{-0.19}^{+0.18}$ and a median gas depletion time $\rm t_{depl.}(Gyr)=M_{gas}/SFR=0.68_{-0.08}^{+0.07}$; correction for the location on the MS will only slightly change the values. The reported uncertainties are the $\rm 1σ$ error on the median. Our results are fully consistent with the expected flattening of the redshift evolution from the 2-SFM (2 star formation mode) framework that empirically prescribes the evolution assuming a universal, log-linear relation between SFR and gas mass coupled to the redshift evolution of the specific star formation rate (sSFR) of main sequence galaxies. While $\rm t_{dep.}$ shows only a mild dependence on location within the MS, a clear trend of increasing $\rm μ_{gas}$ across the MS is observed (as known from previous studies). Further we comment on trends within the MS and (in)consistencies with other studies.

preprint2016arXiv

Gravitational torques imply molecular gas inflow towards the nucleus of M51

The transport of gas towards the centre of galaxies is critical for black hole feeding and, indirectly, it can control active galactic nucleus (AGN) feedback. We have quantified the molecular gas inflow in the central R<1kpc of M51 to be 1 Msun/yr, using a new gravitational torque map and the molecular gas traced by the PdBI Arcsecond Whirlpool Survey (PAWS). The nuclear stellar bar is responsible for this gas inflow. We have also used torque profiles to estimate the location of dynamical resonances, suggesting a corotation for the bar at R~20", and a corotation for the spiral at R~100". We demonstrate how important it is to correct 3.6um images for dust emission in order to compute gravitational torques, and we carefully examine further sources of uncertainty. Our observational measurement of gas inflow can be compared with nuclear molecular outflow rates and provide useful constraints for numerical simulations.

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

Optical Depth Estimates and Effective Critical Densities of Dense Gas Tracers in the Inner Parts of Nearby Galaxy Discs

High critical density molecular lines like HCN(1-0) or HCO+(1-0) represent our best tool to study currently star-forming, dense molecular gas at extragalactic distances. The optical depth of these lines is a key ingredient to estimate the effective density required to excite emission. However, constraints on this quantity are even scarcer in the literature than measurements of the high density tracers themselves. Here, we combine new observations of HCN, HCO+ and HNC(1-0) and their optically thin isotopologues H13CN, H13CO+ and HN13C(1-0) to measure isotopologue line ratios. We use IRAM 30-m observations from the large program EMPIRE and new ALMA observations, which together target 6 nearby star-forming galaxies. Using spectral stacking techniques, we calculate or place strong upper limits on the HCN/H13CN, HCO+/H13CO+ and HNC/HN13C line ratios in the inner parts of these galaxies. Under simple assumptions, we use these to estimate the optical depths of HCN(1-0) and HCO+(1-0) to be τ~2-11 in the active, inner regions of our targets. The critical densities are consequently lowered to values between 5-20$\times 10^5$, 1-3$\times 10^5$ and 9$\times 10^4$ cm-3 for HCN, HCO+ and HNC, respectively. We study the impact of having different beam-filling factors, $η$, on these estimates and find that the effective critical densities decrease by a factor of $\frac{η_{12}}{η_{13}}\,τ_{12}$. A comparison to existing work in NGC 5194 and NGC 253 shows HCN/H13CN and HCO+/H13CO+ ratios in agreement with our measurements within the uncertainties. The same is true for studies in other environments such as the Galactic Centre or nuclear regions of AGN-dominated nearby galaxies.

preprint2016arXiv

The Chandra COSMOS Legacy survey: overview and point source catalog

The COSMOS-Legacy survey is a 4.6 Ms Chandra program that has imaged 2.2 deg$^2$ of the COSMOS field with an effective exposure of $\simeq$160 ks over the central 1.5 deg$^2$ and of $\simeq$80 ks in the remaining area. The survey is the combination of 56 new observations, obtained as an X-ray Visionary Project, with the previous C-COSMOS survey. We describe the reduction and analysis of the new observations and the properties of 2273 point sources detected above a spurious probability of 2$\times 10^{-5}$. We also present the updated properties of the C-COSMOS sources detected in the new data. The whole survey includes 4016 point sources (3814, 2920 and 2440 in the full, soft and hard band). The limiting depths are 2.2 $\times$ 10$^{-16}$, 1.5 $\times$ 10$^{-15}$ and 8.9$\times$ 10$^{-16}$ ${\rm erg~cm}^{-2}~{\rm s}^{-1}$ in the 0.5-2, 2-10 and 0.5-10 keV bands, respectively. The observed fraction of obscured AGN with column density $> 10^{22}$ cm$^{-2}$ from the hardness ratio (HR) is $\sim$50$^{+17}_{-16}$%. Given the large sample, we compute source number counts in the hard and soft bands, significantly reducing the uncertainties of 5-10%. For the first time, we compute number counts for obscured (HR$>$-0.2) and unobscured (HR$<$-0.2) sources and find significant differences between the two populations in the soft band. Due to the un-precedent large exposure, COSMOS-Legacy area is 3 times larger than surveys at similar depth and its depth is 3 times fainter than surveys covering similar area. The area-flux region occupied by COSMOS-Legacy is likely to remain unsurpassed for years to come.

preprint2016arXiv

The effect of unresolved binaries on globular cluster proper-motion dispersion profiles

High-precision kinematic studies of globular clusters require an accurate knowledge of all possible sources of contamination. Amongst other sources, binary stars can introduce systematic biases in the kinematics. Using a set of Monte Carlo cluster simulations with different concentrations and binary fractions, we investigate the effect of unresolved binaries on proper-motion dispersion profiles, treating the simulations like HST proper-motion samples. Since globular clusters evolve towards a state of partial energy equipartition, more massive stars lose energy and decrease their velocity dispersion. As a consequence, on average, binaries have a lower velocity dispersion, since they are more massive kinematic tracers. We show that, in the case of clusters with high binary fraction (initial binary fraction of 50%) and high concentration (i.e., closer to energy equipartition), unresolved binaries introduce a color-dependent bias in the velocity dispersion of main-sequence stars of the order of 0.1-0.3 km s$^{-1}$ (corresponding to $1-6$% of the velocity dispersion), with the reddest stars having a lower velocity dispersion, due to the higher fraction of contaminating binaries. This bias depends on the ability to distinguish binaries from single stars, on the details of the color-magnitude diagram and the photometric errors. We apply our analysis to the HSTPROMO data set of NGC 7078 (M15) and show that no effect ascribable to binaries is observed, consistent with the low binary fraction of the cluster. Our work indicates that binaries do not significantly bias proper-motion velocity-dispersion profiles, but should be taken into account in the error budget of kinematic analyses.

preprint2016arXiv

The EMPIRE Survey: Systematic Variations in the Dense Gas Fraction and Star Formation Efficiency from Full-Disk Mapping of M51

We present the first results from the EMPIRE survey, an IRAM large program that is mapping tracers of high density molecular gas across the disks of nine nearby star-forming galaxies. Here, we present new maps of the 3-mm transitions of HCN, HCO+, and HNC across the whole disk of our pilot target, M51. As expected, dense gas correlates with tracers of recent star formation, filling the "luminosity gap" between Galactic cores and whole galaxies. In detail, we show that both the fraction of gas that is dense, f_dense traced by HCN/CO, and the rate at which dense gas forms stars, SFE_dense traced by IR/HCN, depend on environment in the galaxy. The sense of the dependence is that high surface density, high molecular gas fraction regions of the galaxy show high dense gas fractions and low dense gas star formation efficiencies. This agrees with recent results for individual pointings by Usero et al. 2015 but using unbiased whole-galaxy maps. It also agrees qualitatively with the behavior observed contrasting our own Solar Neighborhood with the central regions of the Milky Way. The sense of the trends can be explained if the dense gas fraction tracks interstellar pressure but star formation occurs only in regions of high density contrast.

preprint2016arXiv

The XXL Survey: IX. Optical overdensity and radio continuum analysis of a supercluster at $z=0.43$

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

preprint2016arXiv

Unveiling The Physics of Star Formation and Feedback in Galaxies

Recent studies show the importance of feedback in the evolution of the star formation rate in the Universe. However, the nature and physics of the feedback are still pressing questions. Radio continuum observations can provide unique dust-unbiased tracers of massive star formation and of the interstellar medium (ISM) and hence are ideal to address the regulation of star formation in galaxies. Our multi-frequency and multi-resolution radio surveys in nearby galaxies enable us to trace various phases of star formation and dissect the thermal and nonthermal ISM in galaxies. Mapping the cosmic ray electron energy index and magnetic field strength, we have found observational evidence that massive star formation significantly affects the energy balance in the ISM through the injection and acceleration of cosmic rays and the amplification of magnetic fields. How the next generation of stars could form in such a magnetized and turbulent ISM will be addressed in our 'EVLA cloud-scale survey of the local group galaxy M33' and in forthcoming surveys with the SKA.

preprint2015arXiv

A direct constraint on the gas content of a massive, passively evolving elliptical galaxy at z = 1.43

Gas and dust in star-forming galaxies at the peak epoch of galaxy assembly are presently the topic of intense study, but little is known about the interstellar medium (ISM) of distant, passively evolving galaxies. We report on a deep 3 mm-band search with IRAM/PdBI for molecular (H$_2$) gas in a massive ($M_{\star}{\sim}6{\times}10^{11}M_{\odot}$) elliptical galaxy at z=1.4277, the first observation of this kind ever attempted. We place a 3$σ$ upper limit of 0.32 Jy km/s on the flux of the CO($J$=2$\rightarrow$1) line or $L'_{\rm CO}$$<$8.8$\times$10$^{9}$ K km/s pc$^2$, assuming a disk-like CO-morphology and a circular velocity scaling with the stellar velocity dispersion as in local early-type galaxies (ETGs). This translates to an H$_2$ mass of $<$3.9$\times$10$^{10}$($α_{\rm CO}$/4.4)$M_{\odot}$ or a gas fraction of $\lesssim$6% assuming a Salpeter initial mass function (IMF) and an ISM dominated by H$_2$, as observed in many local, high-mass ellipticals. This low value approaches that of local ETGs, suggesting that the low star formation activity in massive, high-z passive galaxies reflects a true dearth of gas and a lesser role for inhibitive mechanisms like morphological quenching.

preprint2015arXiv

Heating and cooling of the neutral ISM in the NGC4736 circumnuclear ring

The manner in which gas accretes and orbits within circumnuclear rings has direct implications for the star formation process. In particular, gas may be compressed and shocked at the inflow points, resulting in bursts of star formation at these locations. Afterwards the gas and young stars move together through the ring. In addition, star formation may occur throughout the ring, if and when the gas reaches sufficient density to collapse under gravity. These two scenarios for star formation in rings are often referred to as the `pearls on a string' and `popcorn' paradigms. In this paper, we use new Herschel PACS observations, obtained as part of the KINGFISH Open Time Key Program, along with archival Spitzer and ground-based observations from the SINGS Legacy project, to investigate the heating and cooling of the interstellar medium in the nearby star-forming ring galaxy, NGC4736. By comparing spatially resolved estimates of the stellar FUV flux available for heating, with the gas and dust cooling derived from the FIR continuum and line emission, we show that while star formation is indeed dominant at the inflow points in NGC 4736, additional star formation is needed to balance the gas heating and cooling throughout the ring. This additional component most likely arises from the general increase in gas density in the ring over its lifetime. Our data provide strong evidence, therefore, for a combination of the two paradigms for star formation in the ring in NGC4736.

preprint2015arXiv

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

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

preprint2015arXiv

The Chandra COSMOS Legacy survey: optical/IR identifications

We present the catalog of optical and infrared counterparts of the Chandra COSMOS-Legacy Survey, a 4.6 Ms Chandra program on the 2.2 square degrees of the COSMOS field, combination of 56 new overlapping observations obtained in Cycle 14 with the previous C-COSMOS survey. In this Paper we report the i, K, and 3.6 micron identifications of the 2273 X-ray point sources detected in the new Cycle 14 observations. We use the likelihood ratio technique to derive the association of optical/infrared (IR) counterparts for 97% of the X-ray sources. We also update the information for the 1743 sources detected in C-COSMOS, using new K and 3.6 micron information not available when the C-COSMOS analysis was performed. The final catalog contains 4016 X-ray sources, 97% of which have an optical/IR counterpart and a photometric redshift, while 54% of the sources have a spectroscopic redshift. The full catalog, including spectroscopic and photometric redshifts and optical and X-ray properties described here in detail, is available online. We study several X-ray to optical (X/O) properties: with our large statistics we put better constraints on the X/O flux ratio locus, finding a shift towards faint optical magnitudes in both soft and hard X-ray band. We confirm the existence of a correlation between X/O and the the 2-10 keV luminosity for Type 2 sources. We extend to low luminosities the analysis of the correlation between the fraction of obscured AGN and the hard band luminosity, finding a different behavior between the optically and X-ray classified obscured fraction.

preprint2015arXiv

The Quasar-LBG Two-point Angular Cross-correlation Function at z ~ 4 in the COSMOS Field

In order to investigate the origin of quasars, we estimate the bias factor for low-luminosity quasars at high redshift for the first time. In this study, we use the two-point angular cross-correlation function (CCF) for both low-luminosity quasars at $-24<M_{\rm 1450}<-22$ and Lyman-break galaxies (LBGs). Our sample consists of both 25 low-luminosity quasars (16 objects are spectroscopically confirmed low-luminosity quasars) in the redshift range $3.1<z<4.5$ and 835 color-selected LBGs with $z^{\prime}_{\rm LBG}<25.0$ at $z\sim4$ in the COSMOS field. We have made our analysis for the following two quasar samples; (1) the spectroscopic sample (the 16 quasars confirmed by spectroscopy), and (2) the total sample (the 25 quasars including 9 quasars with photometric redshifts). The bias factor for low-luminosity quasars at $z\sim4$ is derived by utilizing the quasar-LBG CCF and the LBG auto-correlation function. We then obtain the $86\%$ upper limits of the bias factors for low-luminosity quasars, that are 5.63 and 10.50 for the total and the spectroscopic samples, respectively. These bias factors correspond to the typical dark matter halo masses, log $(M_{\rm DM}/(h^{-1}M_{\odot}))=$$12.7$ and $13.5$, respectively. This result is not inconsistent with the predicted bias for quasars which is estimated by the major merger models.

preprint2015arXiv

The Spitzer Survey of Stellar Structure in Galaxies (S$^4$G): Precise Stellar Mass Distributions from Automated Dust Correction at 3.6 microns

The mid-infrared is an optimal window to trace stellar mass in nearby galaxies and the 3.6$μm$ IRAC band has been exploited to this effect, but such mass estimates can be biased by dust emission. We present our pipeline to reveal the old stellar flux at 3.6$μm$ and obtain stellar mass maps for more than 1600 galaxies available from the Spitzer Survey of Stellar Structure in Galaxies (S$^{4}$G). This survey consists of images in two infrared bands (3.6 and 4.5$μm$), and we use the Independent Component Analysis (ICA) method presented in Meidt et al. (2012) to separate the dominant light from old stars and the dust emission that can significantly contribute to the observed 3.6$μm$ flux. We exclude from our ICA analysis galaxies with low signal-to-noise ratio (S/N < 10) and those with original [3.6]-[4.5] colors compatible with an old stellar population, indicative of little dust emission (mostly early Hubble types, which can directly provide good mass maps). For the remaining 1251 galaxies to which ICA was successfully applied, we find that as much as 10-30% of the total light at 3.6$μm$ typically originates from dust, and locally it can reach even higher values. This contamination fraction shows a correlation with specific star formation rates, confirming that the dust emission that we detect is related to star formation. Additionally, we have used our large sample of mass estimates to calibrate a relationship of effective mass-to-light ratio ($M/L$) as a function of observed [3.6]-[4.5] color: $\log(M/L)=-0.339 (\pm 0.057) \times ([3.6]-[4.5]) -0.336 (\pm 0.002)$. Our final pipeline products have been made public through IRSA, providing the astronomical community with an unprecedentedly large set of stellar mass maps ready to use for scientific applications.

preprint2015arXiv

The Subaru COSMOS 20: Subaru Optical Imaging of the HST COSMOS Field with 20 Filters

We present both the observations and the data reduction procedures of the Subaru COSMOS 20 project that is an optical imaging survey of the HST COSMOS field, carried out by using Suprime-Cam on the Subaru Telescope with the following 20 optical filters: 6 broad-band (B, g', V, r', i', and z'), 2 narrow-band (NB711 and NB816), and 12 intermediate-band filters (IA427, IA464, IA484, IA505, IA527, IA574, IA624, IA679, IA709, IA738, IA767, and IA827). A part of this project is described in Taniguchi et al. (2007) and Capak et al. (2007) for the six broad-band and one narrow-band (NB816) filter data. In this paper, we present details of the observations and data reduction for remaining 13 filters (the 12 IA filters and NB711). In particular, we describe the accuracy of both photometry and astrometry in all the filter bands. We also present optical properties of the Suprime-Cam IA filter system in Appendix.

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

ALMA observations of feeding and feedback in nearby Seyfert galaxies: an AGN-driven outflow in NGC 1433

We report ALMA observations of CO(3-2) emission in the Seyfert 2 double-barred galaxy NGC1433, at the unprecedented spatial resolution of 0.5"=24 pc. Our aim is to probe AGN feeding and feedback phenomena through the morphology and dynamics of the gas inside the central kpc. The CO map, which covers the whole nuclear region (nuclear bar and ring), reveals a nuclear gaseous spiral structure, inside the nuclear ring encircling the nuclear stellar bar. This gaseous spiral is well correlated with the dusty spiral seen in Hubble Space Telescope images. The nuclear spiral winds up in a pseudo-ring at 200 pc radius, which might correspond to the inner ILR. Continuum emission is detected at 0.87 mm only at the very centre, and its origin is more likely thermal dust emission than non-thermal emission from the AGN. It might correspond to the molecular torus expected to exist in this Seyfert 2 galaxy. The HCN(4-3) and HCO+(4-3) lines were observed simultaneously, but only upper limits are derived, with a ratio to the CO(3-2) line lower than 1/60 at 3sigma, indicating a relatively low abundance of very dense gas. The kinematics of the gas over the nuclear disk reveal rather regular rotation only slightly perturbed by streaming motions due to the spiral; the primary and secondary bars are too closely aligned with the galaxy major or minor axis to leave a signature in the projected velocities. Near the nucleus, there is an intense high-velocity CO emission feature redshifted to 200 km/s (if located in the plane), with a blue-shifted counterpart, at 2" (100 pc) from the centre. While the CO spectra are quite narrow in the centre, this wide component is interpreted as an outflow involving a molecular mass of 3.6 10^6 Mo and a flow rate 7 Mo/yr. The flow could be in part driven by the central star formation, but is mainly boosted by the AGN through its wind or radio jets.

preprint2014arXiv

ALMA reveals the feeding of the Seyfert 1 nucleus in NGC 1566

We report ALMA observations of CO(3-2) emission in the Seyfert 1 galaxy NGC 1566, at a spatial resolution of 25 pc. Our aim is to investigate the morphology and dynamics of the gas inside the central kpc, and to probe nuclear fueling and feedback phenomena. NGC 1566 has a nuclear bar of 1.7 kpc radius and a conspicuous grand design spiral starting from this radius. The ALMA field of view, of diameter 0.9 kpc, lies well inside the nuclear bar and reveals a molecular trailing spiral structure from 50 to 300~pc in size, which is contributing to fuel the nucleus, according to its negative gravity torques. The spiral starts with a large pitch angle from the center and then winds up in a pseudo-ring at the inner Lindblad resonance (ILR) of the nuclear bar. This is the first time that a trailing spiral structure is clearly seen driving the gas inwards inside the ILR ring of the nuclear bar. This phenomenon shows that the massive central black hole has a significant dynamical influence on the gas, triggering its fueling. The gaseous spiral is well correlated with the dusty spiral seen through extinction in HST images, and also with a spiral feature emitting 0.87mm continuum. This continuum emission must come essentially from cold dust heated by the interstellar radiation field. The HCN(4-3) and HCO+(4-3) lines were simultaneously mapped and detected in the nuclear spiral. The HCO+(4-3) line is 3 times stronger than the HCN(4-3), as expected when star formation excitation dominates over active galactic nucleus (AGN) heating. The CO(3-2)/HCO+(4-3) integrated intensity ratio is \sim 100. The molecular gas is in remarkably regular rotation, with only slight non-circular motions at the periphery of the nuclear spiral arms. These perturbations are quite small, and no outflow nor AGN feedback is detected.

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

Molecular line emission in NGC1068 imaged with ALMA: II. The chemistry of the dense molecular gas

We present a detailed analysis of ALMA Bands 7 and 9 data of CO, HCO+, HCN and CS, augmented with Plateau de Bure Interferometer (PdBI) data of the ~ 200 pc circumnuclear disk (CND) and the ~ 1.3 kpc starburst ring (SB ring) of NGC~1068, a nearby (D = 14 Mpc) Seyfert 2 barred galaxy. We aim at determining the physical characteristics of the dense gas present in the CND and whether the different line intensity ratios we find within the CND as well as between the CND and the SB ring are due to excitation effects (gas density and temperature differences) or to a different chemistry. We estimate the column densities of each species in Local Thermodynamic Equilibrium (LTE). We then compute large one-dimensional non-LTE radiative transfer grids (using RADEX) by using first only the CO transitions, and then all the available molecules in order to constrain the densities, temperatures and column densities within the CND. We finally present a preliminary set of chemical models to determine the origin of the gas. We find that in general the gas in the CND is very dense (> 10^5 cm^-3) and hot (T> 150K), with differences especially in the temperature across the CND. The AGN position has the lowest CO/HCO+, CO/HCN and CO/CS column density ratios. RADEX analyses seem to indicate that there is chemical differentiation across the CND. We also find differences between the chemistry of the SB ring and some regions of the CND; the SB ring is also much colder and less dense than the CND. Chemical modelling does not succeed in reproducing all the molecular ratios with one model per region, suggesting the presence of multi-gas phase components. The LTE, RADEX and chemical analyses all indicate that more than one gas-phase component is necessary to uniquely fit all the available molecular ratios within the CND.

preprint2014arXiv

Molecular line emission in NGC1068 imaged with ALMA. I An AGN-driven outflow in the dense molecular gas

We investigate the fueling and the feedback of star formation and nuclear activity in NGC1068, a nearby (D=14Mpc) Seyfert 2 barred galaxy, by analyzing the distribution and kinematics of the molecular gas in the disk. We have used ALMA to map the emission of a set of dense molecular gas tracers (CO(3-2), CO(6-5), HCN(4-3), HCO+(4-3) and CS(7-6)) and their underlying continuum emission in the central r ~ 2kpc of NGC1068 with spatial resolutions ~ 0.3"-0.5" (~ 20-35pc). Molecular line and dust continuum emissions are detected from a r ~ 200pc off-centered circumnuclear disk (CND), from the 2.6kpc-diameter bar region, and from the r ~ 1.3kpc starburst (SB) ring. Most of the emission in HCO+, HCN and CS stems from the CND. Molecular line ratios show dramatic order-of-magnitude changes inside the CND that are correlated with the UV/X-ray illumination by the AGN, betraying ongoing feedback. The gas kinematics from r ~ 50pc out to r ~ 400pc reveal a massive (M_mol ~ 2.7 (+0.9, -1.2) x 10^7 Msun) outflow in all molecular tracers. The tight correlation between the ionized gas outflow, the radio jet and the occurrence of outward motions in the disk suggests that the outflow is AGN-driven. The outflow rate estimated in the CND, dM/dt ~ 63 (+21, -37) Msun yr^-1, is an order of magnitude higher than the star formation rate at these radii, confirming that the outflow is AGN-driven. The power of the AGN is able to account for the estimated momentum and kinetic luminosity of the outflow. The CND mass load rate of the CND outflow implies a very short gas depletion time scale of <=1 Myr.

preprint2014arXiv

Multi-molecule ALMA observations towards the Seyfert 1 galaxy NGC 1097

The nearby Sy 1 galaxy NGC 1097 represents an ideal laboratory to explore the molecular chemistry in the presence and surroundings of an active galactic nucleus. Exploring the distribution of different molecular species allows us to understand the physical processes affecting the ISM both in the AGN vicinity as well as in the outer star forming molecular ring. We carried out 3 mm ALMA observations of HCN, HCO+, CCH, CS, HNCO, SiO, HC3N, and SO as well as the 13C isotopologues. All species were imaged over the central 2 kpc (~30") of the galaxy at a resolution of ~2.2"x1.5 (150 pc x 100 pc). HCO+ and CS appear to be slightly enhanced in the star forming ring. CCH, showing the largest variations across NGC 1097, is suggested to be a good tracer of both obscured and early stage star formation. HNCO, SiO and HC3N are significantly enhanced in the inner circumnuclear disk surrounding the AGN. Differences in the molecular abundances are observed between the star forming ring and the inner circumnuclear disk. We conclude that the HCN/HCO+ and HCN/CS differences observed between AGN dominated and starburst galaxies are not due to a HCN enhancement due to X-rays, but rather this enhancement is produced by shocked material at distances of 200 pc from the AGN. Additionally we claim the lower HCN/CS to be a combination of a small under-abundance of CS in AGNs together with excitation effects, where a high dense gas component (~10^6 cm^-3) may be more prominent in SB galaxies. However the most promising are the differences found among the dense gas tracers which, at our modest spatial resolution, seem to outline the physical structure of the molecular disk around the AGN. In this picture, HNCO probes the well shielded gas in the disk, surrounding the dense material moderately exposed to X-ray radiation traced by HC3N. Finally SiO might be the innermost molecule in the disk structure.

preprint2014arXiv

Near-Infrared long-slit spectra of Seyfert Galaxies: gas excitation across the central kiloparsec

The excitation of the gas phase of the interstellar medium can be driven by various mechanisms. In galaxies with an active nucleus, such as Seyfert galaxies, both radiative and mechanical energy from the central black hole, or the stars in the disk surrounding it may play a role. We investigate the relative importance and range of influence of the active galactic nucleus for the excitation of ionized and molecular gas in the central kiloparsec of its host galaxy. We present H- and K-band long-slit spectra for a sample of 21 nearby (D <70 Mpc) Seyfert galaxies obtained with the NIRSPEC instrument on the Keck telescope. For each galaxy, we fit the nebular line emission, stellar continua, and warm molecular gas as a function of distance from the nucleus. Our analysis does not reveal a clear difference between the nucleus proper and off-nuclear environment in terms of excitation mechanisms, suggesting that the influence of an AGN reaches far into the disk of the host galaxy. The radial variations in emission line ratios indicate that, while local mechanisms do affect the gas excitation, they are often averaged out when measuring over extended regions. The fully calibrated long-slit spectra, as well as our fitting results, are made available on-line.

preprint2014arXiv

SKA studies of nearby galaxies: star-formation, accretion processes and molecular gas across all environments

The SKA will be a transformational instrument in the study of our local Universe. In particular, by virtue of its high sensitivity (both to point sources and diffuse low surface brightness emission), angular resolution and the frequency ranges covered, the SKA will undertake a very wide range of astrophysical research in the field of nearby galaxies. By surveying vast numbers of nearby galaxies of all types with $μ$Jy sensitivity and sub-arcsecond angular resolutions at radio wavelengths, the SKA will provide the cornerstone of our understanding of star-formation and accretion activity in the local Universe. In this chapter we outline the key continuum and molecular line science areas where the SKA, both during phase-1 and when it becomes the full SKA, will have a significant scientific impact.

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.

preprint2013arXiv

A Detailed Study of the Radio--FIR Correlation in NGC6946 with Herschel-PACS/SPIRE from KINGFISH

We derive the distribution of the synchrotron spectral index across NGC6946 and investigate the correlation between the radio continuum (synchrotron) and far-infrared (FIR) emission using the KINGFISH Herschel PACS and SPIRE data. The radio--FIR correlation is studied as a function of star formation rate, magnetic field strength, radiation field strength, and the total gas surface brightness. The synchrotron emission follows both star-forming regions and the so-called magnetic arms present in the inter-arm regions. The synchrotron spectral index is steepest along the magnetic arms ($α_n \sim 1$), while it is flat in places of giant H{\sc ii} regions and in the center of the galaxy ($α_n \sim 0.6-0.7$). The map of $α_n$ provides an observational evidence for aging and energy loss of cosmic ray electrons propagating in the disk of the galaxy. Variations in the synchrotron--FIR correlation across the galaxy are shown to be a function of both star formation and magnetic fields. We find that the synchrotron emission correlates better with cold rather than with warm dust emission, when the interstellar radiation field is the main heating source of dust. The synchrotron--FIR correlation suggests a coupling between the magnetic field and the gas density. NGC6946 shows a power-law behavior between the total (turbulent) magnetic field strength B and the star formation rate surface density $Σ_{\rm SFR}$ with an index of 0.14\,(0.16)$\pm$0.01. This indicates an efficient production of the turbulent magnetic field with the increasing gas turbulence expected in actively star forming regions. The scale-by-scale analysis of the synchrotron--FIR correlation indicates that the ISM affects the propagation of old/diffused cosmic ray electrons, resulting in a diffusion coefficient of $D_0=4.6\times 10^{28}$\,cm$^2$\,s$^{-1}$ for 2.2\,GeV CREs.

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

Evolution of Galaxies and their Environments at z = 0.1 to 3 in COSMOS

Large-scale structures (LSS) out to z $< 3.0$ are measured in the Cosmic Evolution Survey (COSMOS) using extremely accurate photometric redshifts (photoz). The Ks-band selected sample (from Ultra-Vista) is comprised of 155,954 galaxies. Two techniques -- adaptive smoothing and Voronoi tessellation -- are used to estimate the environmental densities within 127 redshift slices. Approximately 250 statistically significant overdense structures are identified out to z $= 3.0$ with shapes varying from elongated filamentary structures to more circularly symmetric concentrations. We also compare the densities derived for COSMOS with those based on semi-analytic predictions for a $Λ$CDM simulation and find excellent overall agreement between the mean densities as a function of redshift and the range of densities. The galaxy properties (stellar mass, spectral energy distributions (SEDs) and star formation rates (SFRs)) are strongly correlated with environmental density and redshift, particularly at z $< 1.0 - 1.2$. Classifying the spectral type of each galaxy using the rest-frame b-i color (from the photoz SED fitting), we find a strong correlation of early type galaxies (E-Sa) with high density environments, while the degree of environmental segregation varies systematically with redshift out to z $\sim 1.3$. In the highest density regions, 80% of the galaxies are early types at z=0.2 compared to only 20% at z = 1.5. The SFRs and the star formation timescales exhibit clear environmental correlations. At z $> 0.8$, the star formation rate density (SFRD) is uniformly distributed over all environmental density percentiles, while at lower redshifts the dominant contribution is shifted to galaxies in lower density environments.

preprint2013arXiv

Explaining two circumnuclear star forming rings in NGC5248

The distribution of gas in the central kiloparsec of a galaxy has a dynamically rapid evolution. Nonaxisymmetries in the gravitational potential of the galactic disk, such as a large scale stellar bar or spiral, can lead to significant radial motion of gaseous material from larger radii to the central region. The large influx of gas and the subsequent star formation keep the central region constantly changing. However, the ability of gas to reach the nucleus proper to fuel an AGN phase is not guaranteed. Gas inflow can be halted at a circumnuclear star forming ring several hundred parsec away. The nearby galaxy NGC5248 is especially interesting in this sense since it is said to host 2 circumnuclear star forming rings at 100pc and 370pc from its quiescent nucleus. Here we present new subarcsecond PdBI+30m CO(2-1) emission line observations of the central region. For the first time the molecular gas distribution at the smallest stellar ring is resolved into a gas ring, consistent with the presence of a quiescent nucleus. However, the molecular gas shows no ring structure at the larger ring. We combine analyses of the gaseous and stellar content in the central kiloparsec of this galaxy to understand the gas distribution and dynamics of this star forming central region. We discuss the probability of two scenarios leading to the current observations, given our full understanding of this system, and discuss whether there are really two circumnuclear star forming rings in this galaxy.

preprint2013arXiv

How does star formation proceed in the circumnuclear starburst ring of NGC 6951?

Gas inflowing along stellar bars is often stalled at the location of circumnuclear rings, that form an effective reservoir for massive star formation and thus shape the central regions of galaxies. However, how exactly star formation is proceeding within these circumnuclear starburst rings is subject of debate. Two main scenarios for this process have been put forward: In the first the onset of star formation is regulated by the total amount of gas present in the ring with star forming starting once a mass threshold has reached in a `random' position within the ring like `popcorn'. In the second star formation preferentially takes place near the locations where the gas enters the ring. This scenario has been dubbed `pearls-on-a-string'. Here we combine new optical IFU data covering the full stellar bar with existing multi-wavelength data to study in detail the 580 pc radius circumnuclear starburst ring in the nearby spiral galaxy NGC 6951. Using HST archival data together with Sauron and Oasis IFU data, we derive the ages and stellar masses of star clusters as well as the total stellar content of the central region. Adding information on the molecular gas distribution, stellar and gaseous dynamics and extinction, we find that the circumnuclear ring in NGC 6951 is ~1-1.5 Gyr old and has been forming stars for most of that time. We see evidence for preferred sites of star formation within the ring, consistent with the `pearls-on-a-string' scenario, when focusing on the youngest stellar populations. Due to the ring's longevity this signature is washed out when older stellar populations are included in the analysis.

preprint2013arXiv

Mapping dust through emission and absorption in nearby galaxies

Dust has long been identified as a barrier to measuring inherent galaxy properties. However, the link between dust and attenuation is not straightforward and depends on both the amount of dust and its distribution. Herschel imaging of nearby galaxies undertaken as part of the KINGFISH project allows us to map the dust as seen in emission with unprecedented sensitivity and ~1 kpc resolution. We present here new optical integral field unit spectroscopy for eight of these galaxies that provides complementary 100-200 pc scale maps of the dust attenuation through observation of the reddening in both the Balmer decrement and the stellar continuum. The stellar continuum reddening, which is systematically less than that observed in the Balmer decrement, shows no clear correlation with the dust, suggesting that the distribution of stellar reddening acts as a poor tracer of the overall dust content. The brightest HII regions are observed to be preferentially located in dusty regions, and we do find a correlation between the Balmer line reddening and the dust mass surface density for which we provide an empirical relation. Some of the high-inclination systems in our sample exhibit high extinction, but we also find evidence that unresolved variations in the dust distribution on scales smaller than 500 pc may contribute to the scatter in this relation. We caution against the use of integrated A_V measures to infer global dust properties.

preprint2013arXiv

Multi-scale Radio-IR Correlations in M31 and M33: The Role of Magnetic Fields and Star Formation

Interstellar magnetic fields and the propagation of cosmic ray electrons have an important impact on the radio-infrared (IR) correlation in galaxies. This becomes evident when studying different spatial scales within galaxies. We investigate the correlation between the infrared (IR) and free-free/synchrotron radio continuum emission at 20 cm from the two local group galaxies M31 and M33 on spatial scales between 0.4 and 10 kpc. The multi-scale radio-IR correlations have been carried out using a wavelet analysis. The free-free and IR emission are correlated on all scales, but on some scales the synchrotron emission is only marginally correlated with the IR emission. The synchrotron-IR correlation is stronger in M33 than in M31 on small scales (<1 kpc), but it is weaker than in M31 on larger scales. Taking the smallest scale on which the synchrotron-IR correlation exists as the propagation length of cosmic ray electrons, we show that the difference on small scales can be explained by the smaller propagation length in M33 than in M31. On large scales, the difference is due to the thick disk/halo in M33, which is absent in M31. A comparison of our data with data on NGC6946, the LMC and M51 suggests that the propagation length is determined by the ratio of ordered-to-turbulent magnetic field strength, which is consistent with diffusion of CR electrons in the ISM. As the diffusion length of CR electrons influences the radio-IR correlation, this dependence is a direct observational evidence of the importance of magnetic fields for the radio-IR correlation within galaxies. The star formation rate per surface area only indirectly influences the diffusion length as it increases the strength of the turbulent magnetic field.

preprint2013arXiv

Pathway to the Square Kilometre Array - The German White Paper -

The Square Kilometre Array (SKA) is the most ambitious radio telescope ever planned. With a collecting area of about a square kilometre, the SKA will be far superior in sensitivity and observing speed to all current radio facilities. The scientific capability promised by the SKA and its technological challenges provide an ideal base for interdisciplinary research, technology transfer, and collaboration between universities, research centres and industry. The SKA in the radio regime and the European Extreme Large Telescope (E-ELT) in the optical band are on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI) and have been recognised as the essential facilities for European research in astronomy. This "White Paper" outlines the German science and R&D interests in the SKA project and will provide the basis for future funding applications to secure German involvement in the Square Kilometre Array.

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 circumnuclear environment of NGC613: a nuclear starburst caught in the act?

We present near-infrared ($H$- and $K$-band) integral-field observations of the inner $\sim$700pc of the active spiral galaxy NGC613, obtained with SINFONI on the Very Large Telescope. We use emission-line ratios to determine the dominant excitation mechanisms in different regions within our field-of-view, in particular the active nucleus and the star-forming circum-nuclear ring. Diagnostic diagrams involving [FeII] and H$_2$ fluxes indicate that the gas is not only photoionized by the AGN in the nucleus of NGC613, but also shock-heated. On the other hand, the emission line ratios measured in the ``hot spots'' along the ring are fully consistent with them being young star forming regions. We find no sign of radial gas transport from the ring into the core region dominated by the AGN. The ring morphology appears disturbed by a radial outflow of material from the AGN, which is confirmed by the existence of a weak jet in archival radio maps. However, this jet does not seem to have any significant effect on the morphology of the large ($\sim$8$\times$10$^7$ solar masses) reservoir of molecular gas that has accumulated inside the central $\sim$100pc. Such a concentration of molecular gas around an AGN is unusual, and supports a scenario in which star formation is recurrent and episodic in spiral galaxies. In this context, NGC613 appears to be in final stages of the gas accumulation phase, and is likely to undergo a nuclear starburst in the near future.

preprint2013arXiv

The impact of bars on disk breaks as probed by S4G imaging

We have analyzed the radial distribution of old stars in a sample of 218 nearby face-on disks, using deep 3.6um images from the Spitzer Survey of Stellar Structure in Galaxies (S4G). In particular, we have studied the structural properties of those disks with a broken or down-bending profile. We find that, on average, disks with a genuine single exponential profile have a scale-length and a central surface brightness which are intermediate to those of the inner and outer components of a down-bending disk with the same total stellar mass. In the case of barred galaxies, the ratio between the break and the bar radii (Rbr/Rbar) depends strongly on the total stellar mass of the galaxy. For galaxies more massive than 10^10 Msun, the distribution is bimodal, peaking at Rbr/Rbar~2 and ~3.5. The first peak, which is the most populated one, is linked to the Outer Lindblad Resonance of the bar, whereas the second one is consistent with a dynamical coupling between the bar and the spiral pattern. For galaxies below 10^10 Msun, breaks are found up to ~10 Rbar, but we show that they could still be caused by resonances given the rising nature of rotation curves in these low-mass disks. While not ruling out star formation thresholds, our results imply that radial stellar migration induced by non-axysymmetric features can be responsible not only for those breaks at 2 Rbar, but also for many of those found at larger radii.

preprint2013arXiv

The Plateau de Bure + 30m Arcsecond Whirlpool Survey reveals a thick disk of diffuse molecular gas in the M51 galaxy

We present the data of the Plateau de Bure Arcsecond Whirlpool Survey (PAWS), a high spatial and spectral resolution 12CO(1-0) line survey of the inner 10 x 6 kpc of the M51 system, and the first wide-field imaging of molecular gas in a star-forming spiral galaxy with resolution matched to the typical size of Giant Molecular Clouds (40 pc). We describe the observation, reduction, and combination of the Plateau de Bure Interferometer (PdBI) and IRAM-30m "short spacing" data. The final data cube attains 1.1"-resolution over the 270" x 170" field of view, with sensitivity to all spatial scales from the combination of PdBI and IRAM-30m data, and brightness sensitivity of 0.4 K (1 sigma) in each 5 km/s-wide channel map. We find a CO-luminosity of 9.10^8 K km/s pc^2, corresponding to a molecular gas mass of 4.10^9 Msol for a standard CO-to-H2 conversion factor. Unexpectedly, we find that a large fraction, (50 +/- 10)%, of this emission arises mostly from spatial scales larger than 36" or about 1.3 kpc. Through a series of tests, we demonstrate that this extended emission does not result from a processing artifact. We discuss its origin in light of the stellar component, the 12CO/13CO ratio, and the difference between the kinematics and structure of the PdBI-only and hybrid synthesis (PdBI + IRAM-30m) images. The extended emission is consistent with a thick, diffuse disk of molecular gas with a typical scale height of about 200 pc, substructured in unresolved filaments which fills about 0.1% of the volume.

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

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

Constraints on the faint end of the quasar luminosity function at z~5 in the COSMOS field

We present the result of our low-luminosity quasar survey in the redshift range of 4.5 < z < 5.5 in the COSMOS field. Using the COSMOS photometric catalog, we selected 15 quasar candidates with 22 < i' < 24 at z~5, that are ~ 3 mag fainter than the SDSS quasars in the same redshift range. We obtained optical spectra for 14 of the 15 candidates using FOCAS on the Subaru Telescope and did not identify any low-luminosity type-1 quasars at z~5 while a low-luminosity type-2 quasar at z~5.07 was discovered. In order to constrain the faint end of the quasar luminosity function at z~5, we calculated the 1sigma confidence upper limits of the space density of type-1 quasars. As a result, the 1sigma confidence upper limits on the quasar space density are Phi< 1.33*10^{-7} Mpc^{-3} mag^{-1} for -24.52 < M_{1450} < -23.52 and Phi< 2.88*10^{-7} Mpc^{-3} mag^{-1} for -23.52 < M_{1450} < -22.52. The inferred 1sigma confidence upper limits of the space density are then used to provide constrains on the faint-end slope and the break absolute magnitude of the quasar luminosity function at z~5. We find that the quasar space density decreases gradually as a function of redshift at low luminosity (M_{1450} ~ -23), being similar to the trend found for quasars with high luminosity (M_{1450}<-26). This result is consistent with the so-called downsizing evolution of quasars seen at lower redshifts.

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

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

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

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

preprint2012arXiv

Physical properties of dense molecular gas in centres of Seyfert galaxies

We present new ~1" resolution data of the dense molecular gas in the central 50-100 pc of four nearby Seyfert galaxies. PdBI observations of HCN and, in 2 of the 4 sources, simultaneously HCO+ allow us to carefully constrain the dynamical state of the dense gas surrounding the AGN. Analysis of the kinematics shows large line widths of 100-200 km/s FWHM that can only partially arise from beam smearing of the velocity gradient. The observed morphological and kinematic parameters (dimensions, major axis position angle, red and blue channel separation, and integrated line width) are well reproduced by a thick disk, where the emitting dense gas has a large intrinsic dispersion (20-40 km/s), implying that it exists at significant scale heights (25-30% of the disk radius). To put the observed kinematics in the context of the starburst and AGN evolution, we estimate the Toomre Q parameter. We find this is always greater than the critical value, i.e. Q is above the limit such that the gas is stable against rapid star formation. This is supported by the lack of direct evidence, in these 4 Seyfert galaxies, for on-going star formation close around the AGN. Instead, any current star formation tends to be located in a circumnuclear ring. We conclude that the physical conditions are indeed not suited to star formation within the central ~100 pc.

preprint2012arXiv

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

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

preprint2012arXiv

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

Spectral Energy Distributions of Type 1 AGN in the COSMOS Survey I - The XMM-COSMOS Sample

The "Cosmic Evolution Survey" (COSMOS) enables the study of the Spectral Energy Distributions (SEDs) of Active Galactic Nuclei (AGN) because of the deep coverage and rich sampling of frequencies from X-ray to radio. Here we present a SED catalog of 413 X-ray (\xmm) selected type 1 (emission line FWHM$>2000$ km s$^{-1}$) AGN with Magellan, SDSS or VLT spectrum. The SEDs are corrected for the Galactic extinction, for broad emission line contributions, constrained variability, and for host galaxy contribution. We present the mean SED and the dispersion SEDs after the above corrections in the rest frame 1.4 GHz to 40 keV, and show examples of the variety of SEDs encountered. In the near-infrared to optical (rest frame $\sim 8μm$-- 4000Å), the photometry is complete for the whole sample and the mean SED is derived from detections only. Reddening and host galaxy contamination could account for a large fraction of the observed SED variety. The SEDs are all available on-line.

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 Star Formation in Radio Survey: GBT 33 GHz Observations of Nearby Galaxy Nuclei and Extranuclear Star-Forming Regions

We present 33\,GHz photometry of 103 galaxy nuclei and extranuclear star-forming complexes taken with the Green Bank Telescope (GBT) as part of the Star Formation in Radio Survey (SFRS). Among the sources without evidence for an AGN, and also having lower frequency radio data, we find a median thermal fraction at 33GHz of ~76% with a dispersion of ~24%. For all sources resolved on scales <0.5kpc, the thermal fraction is even larger, being >90%. This suggests that the rest-frame 33GHz emission provides a sensitive measure of the ionizing photon rate from young star-forming regions, thus making it a robust star formation rate indicator. Taking the 33GHz star formation rates as a reference, we investigate other empirical calibrations relying on different combinations of warm 24μm dust, total infrared (IR; 8-1000μm), Hα line, and far-UV continuum emission. The recipes derived here generally agree with others found in the literature, albeit with a large dispersion that most likely stems from a combination of effects. Comparing the 33GHz to total IR flux ratios as a function of the radio spectral index, measured between 1.7 and 33GHz, we find that the ratio increases as the radio spectral index flattens which does not appear to be a distance effect. Consequently, the ratio of non-thermal to total IR emission appears relatively constant, suggesting only moderate variations in the cosmic-ray electron injection spectrum and ratio of synchrotron to total cooling processes among star-forming complexes. Assuming that this trend solely arises from an increase in the thermal fraction sets a maximum on the scatter of the non-thermal spectral indices among the star-forming regions of σ_α^{NT} < 0.13.

preprint2012arXiv

X-ray properties of radio-selected star forming galaxies in the Chandra-COSMOS survey

X-ray surveys contain sizable numbers of star forming galaxies, beyond the AGN which usually make the majority of detections. Many methods to separate the two populations are used in the literature, based on X-ray and multiwavelength properties. We aim at a detailed test of the classification schemes and to study the X-ray properties of the resulting samples. We build on a sample of galaxies selected at 1.4 GHz in the VLA-COSMOS survey, classified by Smolcic et al. (2008) according to their optical colours and observed with Chandra. A similarly selected control sample of AGN is also used for comparison. We review some X-ray based classification criteria and check how they affect the sample composition. The efficiency of the classification scheme devised by Smolcic et al. (2008) is such that ~30% of composite/misclassified objects are expected because of the higher X-ray brightness of AGN with respect to galaxies. The latter fraction is actually 50% in the X-ray detected sources, while it is expected to be much lower among X-ray undetected sources. Indeed, the analysis of the stacked spectrum of undetected sources shows, consistently, strongly different properties between the AGN and galaxy samples. X-ray based selection criteria are then used to refine both samples. The radio/X-ray luminosity correlation for star forming galaxies is found to hold with the same X-ray/radio ratio valid for nearby galaxies. Some evolution of the ratio may be possible for sources at high redshift or high luminosity, tough it is likely explained by a bias arising from the radio selection. Finally, we discuss the X-ray number counts of star forming galaxies from the VLA- and C-COSMOS surveys according to different selection criteria, and compare them to the similar determination from the Chandra Deep Fields. The classification scheme proposed here may find application in future works and surveys.

preprint2011arXiv

Calibrating Extinction-Free Star Formation Rate Diagnostics with 33GHz Free-Free Emission in NGC6946

Abridged: Using free-free emission measured in the Ka-band (26-40GHz) for 10 star-forming regions in the nearby galaxy NGC6946, including its starbursting nucleus, we compare a number of SFR diagnostics that are typically considered to be unaffected by interstellar extinction: i.e., non-thermal radio (i.e., 1.4GHz), total infrared (IR; 8-1000um), and warm dust (i.e., 24um) emission, along with the hybrid (obscured + unobscured) indicators of Hα+24um and UV+IR. The 33GHz free-free emission is assumed to provide the most accurate measure of the current SFR. Among the extranuclear star-forming regions, the 24um, Hα+24um and UV+IR SFR calibrations are in good agreement with the 33GHz free-free SFRs. However, each of the SFR calibrations relying on some form of dust emission overestimate the nuclear SFR by a factor of ~2. This is more likely the result of excess dust heating through an accumulation of non-ionizing stars associated with an extended episode of star formation in the nucleus rather than increased competition for ionizing photons by dust. SFR calibrations using the non-thermal radio continuum yield values which only agree with the free-free SFRs for the nucleus, and underestimate the SFRs from the extranuclear star-forming regions by a factor of ~2. This result likely arises from the CR electrons decaying within the starburst region with negligible escape compared to the young extranuclear star-forming regions. Finally, we find that the SFRs estimated using the total 33GHz emission agree well with the free-free SFRs due to the large thermal fractions present at these frequencies even when local diffuse backgrounds are not removed. Thus, rest-frame 33GHz observations may act as a reliable method to measure the SFRs of galaxies at increasingly high redshift without the need of ancillary radio data to account for the non-thermal emission.

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

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

Molecular gas in NUclei of GAlaxies (NUGA) XV. Molecular gas kinematics in the inner 3kpc of NGC6951

Within the NUclei of GAlaxies project we have obtained IRAM PdBI and 30m 12CO(1-0) and 12CO(2-1) observations of the spiral galaxy NGC 6951. Previous work shows that there is indirect evidence of gas inflow from 3 kpc down to small radii: a large-scale stellar bar, a prominent starburst ring (r~580 pc) and a LINER/Seyfert 2 nucleus. In this paper we study the gas kinematics as traced by the CO line emission in detail. We quantify the influence of the large-scale stellar bar by constructing an analytical model of the evolution of gas particles in a barred potential. From this model gravitational torques and mass accumulation rates are computed. We compare our model-based gravitational torque results with previous observationally-based ones. The model also shows that the large-scale stellar bar is indeed the dominant force for driving the gas inward, to the starburst ring. Inside the ring itself a nuclear stellar oval might play an important role. Detailed analysis of the CO gas kinematics there shows that emission arises from two co-spatial, but kinematically distinct components at several locations. The main emission component can always be related to the overall bar-driven gas kinematics. The second component exhibits velocities that are larger than expected for gas on stable orbits, has a molecular gas mass of 1.8x10^6Msun, is very likely connected to the nuclear stellar oval, and is consistent with inflowing motion towards the very center. This may form the last link in the chain of gas inflow towards the active galactic nucleus in NGC 6951.

preprint2011arXiv

On the occupation of X-ray selected galaxy groups by radio AGN since z=1.3

Previous clustering analysis of low-power radio AGN has indicated that they preferentially live in massive groups. The X-ray surveys of the COSMOS field have achieved a sensitivity at which these groups are directly detected out to z=1.3. Making use of Chandra-, XMM- and VLA-COSMOS surveys we identify radio AGN members (10**23.6 < L_1.4GHz/(W/Hz) < 10**25) of galaxy groups (10**13.2 < M_200/M_sun < 10**14.4; 0.1<z<1.3) and study i) the radio AGN -- X-ray group occupation statistics as a function of group mass, and ii) the distribution of radio AGN within the groups. We find that radio AGN are preferentially associated with galaxies close to the center (< 0.2r_200). Compared to our control sample of group members matched in stellar mass and color to the radio AGN host galaxies, we find a significant enhancement of radio AGN activity associated with 10**13.6 < M_200/M_sun < 10**14 halos. We present the first direct measurement of the halo occupation distribution (HOD) for radio AGN, based on the total mass function of galaxy groups hosting radio AGN. Our results suggest a possible deviation from the usually assumed power law HOD model. We also find an overall increase of the fraction of radio AGN in galaxy groups (<1r_200), relative to that in all environments.

preprint2011arXiv

SMA/PdBI multiple line observations of the nearby Seyfert2 galaxy NGC 1068: Shock related gas kinematics and heating in the central 100pc?

We present high angular resolution (0.5-2.0") observations of the mm continuum and the 12CO(J=3-2), 13CO(J=3-2), 13CO(J=2-1), C18O(J=2-1), HCN(J=3-2), HCO+(J=4-3) and HCO+(J=3-2) line emission in the circumnuclear disk (r=100pc) of the proto-typical Seyfert type-2 galaxy NGC1068, carried out with the Submillimeter Array. We further include in our analysis new 13CO(J=1-0) and improved 12CO(J=2-1) observations of NGC1068 at high angular resolution (1.0-2.0") and sensitivity, conducted with the IRAM Plateau de Bure Interferometer. Based on the complex dynamics of the molecular gas emission indicating non-circular motions in the central ~100pc, we propose a scenario in which part of the molecular gas in the circumnuclear disk of NGC1068 is radially blown outwards as a result of shocks. This shock scenario is further supported by quite warm (Tkin>=200K) and dense (nH2=10^4cm^-3) gas constrained from the observed molecular line ratios. The HCN abundance in the circumnuclear disk is found to be [HCN]/[12CO]=10^-3.5. This is slightly higher than the abundances derived for galactic and extragalactic starforming/starbursting regions. This results lends further support to X-ray enhanced HCN formation in the circumnuclear disk of NGC1068, as suggested by earlier studies. The HCO+ abundance ([HCO+]/[12CO]=10^-5) appears to be somewhat lower than that of galactic and extragalactic starforming/starbursting regions. When trying to fit the cm to mm continuum emission by different thermal and non-thermal processes, it appears that electron-scattered synchrotron emission yields the best results while thermal free-free emission seems to over-predict the mm continuum emission.

preprint2011arXiv

Spectroscopy of luminous z>7 galaxy candidates and sources of contamination in z>7 galaxy searches

We present three bright z+ dropout candidates selected from deep Near-Infrared (NIR) imaging of the COSMOS 2 square degree field. All three objects match the 0.8-8um colors of other published z>7 candidates but are three magnitudes brighter, facilitating further study. Deep spectroscopy of two of the candidates covering 0.64-1.02um with Keck-DEIMOS and all three covering 0.94-1.10um and 1.52-1.80um with Keck-NIRSPEC detects weak spectral features tentatively identified as Ly-alpha at z=6.95 and z=7.69 in two of the objects. The third object is placed at z~1.6 based on a 24um and weak optical detection. A comparison with the spectral energy distributions of known z<7 galaxies, including objects with strong spectral lines, large extinction, and large systematic uncertainties in the photometry yields no objects with similar colors. However, the lambda>1um properties of all three objects can be matched to optically detected sources with photometric redshifts at z~1.8, so the non-detection in the i+ and z+ bands are the primary factors which favors a z>7 solution. If any of these objects are at z~7 the bright end of the luminosity function is significantly higher at z>7 than suggested by previous studies, but consistent within the statistical uncertainty and the dark matter halo distribution. If these objects are at low redshift, the Lyman-Break selection must be contaminated by a previously unknown population of low redshift objects with very strong breaks in their broad band spectral energy distributions and blue NIR colors. The implications of this result on luminosity function evolution at high redshift is discussed. We show that the primary limitation of z>7 galaxy searches with broad filters is the depth of the available optical data.

preprint2011arXiv

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

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

preprint2011arXiv

The redshift and nature of AzTEC/COSMOS 1: A starburst galaxy at z=4.6

Based on broad/narrow-band photometry and Keck DEIMOS spectroscopy we report a redshift of z=4.64-0.08+0.06 for AzTEC/COSMOS 1, the brightest sub-mm galaxy in the AzTEC/COSMOS field. In addition to the COSMOS-survey X-ray to radio data, we report observations of the source with Herschel/PACS (100, 160 micron), CSO/SHARC II (350 micron), CARMA and PdBI (3 mm). We do not detect CO(5-4) line emission in the covered redshift ranges, 4.56-4.76 (PdBI/CARMA) and 4.94-5.02 (CARMA). If the line is within this bandwidth, this sets 3sigma upper limits on the gas mass to <~8x10^9 M_Sol and <~5x10^10 M_Sol, respectively (assuming similar conditions as observed in z~2 SMGs). This could be explained by a low CO-excitation in the source. Our analysis of the UV-IR spectral energy distribution of AzTEC 1 shows that it is an extremely young (<~50 Myr), massive (M*~10^11 M_Sol), but compact (<~2 kpc) galaxy forming stars at a rate of ~1300 M_Sol/yr. Our results imply that AzTEC 1 is forming stars in a 'gravitationally bound' regime in which gravity prohibits the formation of a superwind, leading to matter accumulation within the galaxy and further generations of star formation.

preprint2011arXiv

Variation in the dust spectral index across M33

Using the Herschel PACS and SPIRE FIR/submm data, we investigate variations in the dust spectral index $β$ in the nearby spiral galaxy M33 at a linear resolution of 160 pc. We use an iteration method in two different approaches, single and two-component modified black body models. In both approaches, $β$ is higher in the central disk than in the outer disk similar to the dust temperature. There is a positive correlation between $β$ and Halpha as well as with the molecular gas traced by CO(2-1). A Monte-Carlo simulation shows that the physical parameters are better constrained when using the two-component model.

preprint2010arXiv

A Multiwavelength Study of a Sample of 70 micron Selected Galaxies in the COSMOS Field II: The Role of Mergers in Galaxy Evolution

We analyze the morphological properties of a large sample of 1503 70 micron selected galaxies in the COSMOS field spanning the redshift range 0.01<z< 3.5 with a median redshift of 0.5 and an infrared luminosity range of 10^8<L_IR<10^14L_sun with a median luminosity of 10^11.4 L_sun. In general these galaxies are massive, with a stellar mass range of 10^10-10^12 M_sun, and luminous, with -25<M_K<-20. We find a strong correlation between the fraction of major mergers and L_IR, with the fraction at the highest luminosity being up to 50%. We also find that the fraction of spirals drops dramatically with L_IR. Minor mergers likely play a role in boosting the infrared luminosity for sources with low luminosities. The precise fraction of mergers in any given L_IR bin varies by redshift due to sources at z>1 being difficult to classify and subject to the effects of band pass shifting, therefore, these numbers can only be considered lower limits. At z<1, where the morphological classifications are most robust, major mergers clearly dominate the ULIRG population (50-80%) and are important for the LIRG population (25-40%). At z>1 the fraction of major mergers is at least 30-40% for ULIRGs. Although the general morphological trends agree with what has been observed for local (U)LIRGs, the fraction of major mergers is slightly lower than seen locally. This is in part due to the difficulty of identifying merger signatures at high redshift. We argue that given the number of major gas-rich mergers observed and the relatively short timescale that they would be observable in the (U)LIRG phase that it is plausible for the observed red sequence of massive ellipticals (<10^12 M_sun) to have been formed entirely by gas-rich major mergers.

preprint2010arXiv

Cold Molecular Gas in the Inner Two Kiloparsec of NGC4151

We present the first spatially resolved spectroscopic imaging observations of the 12CO(1-0) line emission in the central 2.5 kpc of the Seyfert 1 galaxy NGC4151, obtained with the IRAM Plateau de Bure Interferometer (PdBI). Most of the cold molecular gas is distributed along two curved gas lanes about 1 kpc north and south of the active nucleus, coincident with the circumnuclear dust ring noted by previous authors. These CO arcs lie within the Inner Lindblad Resonance of the large scale oval bar and have kinematics consistent with those derived from neutral hydrogen observations of the disk and bar. Two additional gas clumps are detected that show non-circular motions - one associated with the southern gas lane and one lying ~600 pc north of the nucleus. Closer to the nucleus, no cold molecular gas is detected in the central 300 pc where abundant near-IR H2 line emission arises. This suggests that the H2 line emission is not a good indicator for a cold gas reservoir in NGC4151 and that the H2 is likely photo-excited by the AGN. The upper limit of the CO mass in the central 300 pc is sufficient to support the AGN activity at its current level for 10^7 yrs. The total cold molecular mass detected by PdBI is 4.3 10^7 Msun. Finally, 3 mm continuum emission arising from the location of the AGN is detected with a flux of S~14 mJy and appears to be unresolved at an angular resolution of 2.8" (~180 pc).

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

High-resolution mapping of the physical conditions in two nearby active galaxies based on 12CO(1-0), (2-1) and (3-2) lines

We present a detailed analysis of high resolution observations of the three lowest CO transitions in two nearby active galaxies, NGC4569 and NGC4826. The CO(1-0) and (2-1) lines were observed with the Plateau de Bure Interferometer and the CO(3-2) line with the Submillimeter Array. Combining these data allows us to compare the emission in the three lines and to map the line ratios, R21=I_{CO(2-1)}/I_{CO(1-0)} and R32=I_{CO(3-2)}/I_{CO(1-0)} at a resolution of ~2", i.e., a linear resolution of 160 pc for NGC4569 and 40 pc for NGC4826. In both galaxies the emission in the three lines is similarly distributed spatially and in velocity, and CO is less excited (R32<0.6) than in the Galactic Center or the centers of other active galaxies studied so far. According to a pseudo-LTE model the molecular gas in NGC4569 is cold and mainly optically thick in the CO(1-0) and (2-1) lines; less than 50% of the gas is optically thin in the CO(3-2) line. LVG modeling suggests the presence of an elongated ring of cold and dense gas coinciding with the ILR of the stellar bar. More excited gas is resolved in the circumnuclear disk of NGC4826. According to our pseudo-LTE model this corresponds to warmer gas with a ~50% of the CO(3-2) emission being optically thin. LVG modeling indicates the presence of a semicircular arc of dense and cold gas centered on the dynamical center and ~70 pc in radius. The gas temperature increases and its density decreases toward the center. A near side/far side asymmetry noticeable in the CO, R32 and Pa-alpha maps suggests that opacity effects play a role. Examining published CO maps of nearby active galaxies we find similar asymmetries suggesting that this could be a common phenomenon in active galaxies. These mainly qualitative results open new perspectives for the study of active galaxies with the future Atacama Large Millimeter/submillimeter Array.

preprint2010arXiv

Identifying dynamically young galaxy groups via wide-angle tail galaxies: A case study in the COSMOS field at z=0.53

We present an analysis of a wide-angle tail (WAT) radio galaxy located in a galaxy group in the COSMOS field at a redshift of z=0.53 (hereafter CWAT-02). We find that the host galaxy of CWAT-02 is the brightest galaxy in the group, although it does not coincide with the center of mass of the system. Estimating a) the velocity of CWAT-02, relative to the intra-cluster medium (ICM), and b) the line-of-sight peculiar velocity of CWAT-02's host galaxy, relative to the average velocity of the group, we find that both values are higher than those expected for a dominant galaxy in a relaxed system. This suggests that CWAT-02's host group is dynamically young and likely in the process of an ongoing group merger. Our results are consistent with previous findings showing that the presence of a wide-angle tail galaxy in a galaxy group or cluster can be used as an indicator of dynamically young non-relaxed systems. Taking the unrelaxed state of CWAT-02's host group into account, we discuss the impact of radio-AGN heating from CWAT-02 onto its environment, in the context of the missing baryon problem in galaxy groups. Our analysis strengthens recent results suggesting that radio-AGN heating may be powerful enough to expel baryons from galaxy groups.

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

Molecular gas chemistry in AGN. II. High-resolution imaging of SiO emission in NGC1068: shocks or XDR?

This paper is part of a multi-species survey of line emission from the molecular gas in the circum-nuclear disk (CND) of the Seyfert 2 galaxy NGC1068. Single-dish observations have provided evidence that the abundance of silicon monoxide(SiO) in the CND of NGC1068 is enhanced by 3-4 orders of magnitude with respect to the values typically measured in quiescent molecular gas in the Galaxy. We aim at unveiling the mechanism(s) underlying the SiO enhancement. We have imaged with the IRAM Plateau de Bure interferometer the emission of the SiO(2-1) and CN(2--1) lines in NGC1068 at 150pc and 60pc spatial resolution, respectively. We have also obtained complementary IRAM 30m observations of HNCO and methanol (CH3OH) lines. SiO is detected in a disk of 400pc size around the AGN. SiO abundances in the CND of (1-5)xE-09 are about 1-2 orders of magnitude above those measured in the starburst ring. The overall abundance of CN in the CND is high: (0.2-1)xE-07. The abundances of SiO and CN are enhanced at the extreme velocities of gas associated with non-circular motions close to the AGN (r<70pc). Abundances measured for CN and SiO, and the correlation of CN/CO and SiO/CO ratios with hard X-ray irradiation, suggest that the CND of NGC1068 has become a giant X-ray dominated region (XDR). The extreme properties of molecular gas in the circum-nuclear molecular disk of NGC1068 result from the interplay between different processes directly linked to nuclear activity. Whereas XDR chemistry offers a simple explanation for CN and SiO in NGC1068, the relevance of shocks deserves further scrutiny. The inclusion of dust grain chemistry would help solve the controversy regarding the abundances of other molecular species, like HCN, which are under-predicted by XDR models.

preprint2010arXiv

No Evolution in the IR-Radio Relation for IR-Luminous Galaxies at z<2 in the COSMOS Field

Previous observational studies of the infrared (IR)-radio relation out to high redshift employed any detectable star forming systems at a given redshift within the restricted area of cosmological survey fields. Consequently, the evolution inferred relies on a comparison between the average IR/radio properties of (i) very IR-luminous high-z sources and (ii) more heterogeneous low(er)-z samples that often lack the strongest IR emitters. In this report we consider populations of objects with comparable luminosities over the last 10 Gyr by taking advantage of deep IR (esp. Spitzer 24 micron) and VLA 1.4 GHz observations of the COSMOS field. Consistent with recent model predictions, both Ultra Luminous Infrared Galaxies (ULIRGs) and galaxies on the bright end of the evolving IR luminosity function do not display any change in their average IR/radio ratios out to z~2 when corrected for bias. Uncorrected data suggested ~0.3 dex of positive evolution.

preprint2010arXiv

Probing the faint end of the quasar luminosity function at z ~ 4 in the COSMOS field

We searched for quasars that are ~ 3 mag fainter than the SDSS quasars in the redshift range 3.7 < z < 4.7 in the COSMOS field to constrain the faint end of the quasar luminosity function. Using optical photometric data, we selected 31 quasar candidates with 22 < i' < 24 at z ~ 4. We obtained optical spectra for most of these candidates using FOCAS on the Subaru telescope, and identified 8 low-luminosity quasars at z ~ 4. In order to derive the quasar luminosity function (QLF) based on our spectroscopic follow-up campaign, we estimated the photometric completeness of our quasar survey through detailed Monte Carlo simulations. Our QLF at z ~ 4 has a much shallower faint-end slope beta = -1.67^{+0.11}_{-0.17} than that obtained by other recent surveys in the same redshift. Our result is consistent with the scenario of downsizing evolution of active galactic nuclei inferred by recent optical and X-ray quasar surveys at lower redshifts.

preprint2010arXiv

Radio galaxy feedback in X-ray selected groups from COSMOS: the effect on the ICM

We quantify the importance of the mechanical energy released by radio-galaxies inside galaxy groups. We use scaling relations to estimate the mechanical energy released by 16 radio-AGN located inside X-ray detected galaxy groups in the COSMOS field. By comparing this energy output to the host groups' gravitational binding energy, we find that radio galaxies produce sufficient energy to unbind a significant fraction of the intra-group medium. This unbinding effect is negligible in massive galaxy clusters with deeper potential wells. Our results correctly reproduce the breaking of self-similarity observed in the scaling relation between entropy and temperature for galaxy groups.

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

The local radio-IR correlation in M51

We observed M51 at three frequencies, 1.4GHz (20cm), 4.9GHz (6cm) and 8.4GHz (3.6cm), with the VLA and the Effelsberg 100m telescope to obtain the highest quality radio continuum images of a nearby spiral galaxy. These radio data were combined with deconvolved Spitzer IRAC 8mum and MIPS 24mum images to search for and investigate local changes in the radio-IR correlation. Utilizing wavelet decomposition, we compare the distribution of the radio and IR emission on spatial scales between 200pc and 30kpc. We show that the radio-IR correlation is not uniform across the galactic disk. It presents a complex behavior with local extrema corresponding to various galactic structures, such as complexes of HII regions, spiral arms and interarm filaments, indicating that the contribution of the thermal and non-thermal radio emission is a strong function of environment. In particular, the relation of the 24mum and 20cm emission presents a linear relation within the spiral arms and globally over the galaxy, while it deviates from linearity in the interarm and outer regions as well in the inner region, with two different behaviours: it is sublinear in the interarm and outer region and over-linear in the central 3.5kpc. Our analysis suggests that the changes in the radio/IR correlation reflect variations of ISM properties between spiral arms and interarm region. The good correlation in the spiral arms implies that 24mum and 20cm are tracing recent star formation, while a change in the dust opacity, 'Cirrus' contribution to the IR emission and/or the relation between the magnetic field strength and the gas density can explain the different relations found in the interarm, outer and inner regions.

preprint2010arXiv

The Opacity of Galactic Disks at z~0.7

We compare the surface brightness-inclination relation for a sample of COSMOS pure disk galaxies at z~0.7 with an artificially redshifted sample of SDSS disks well matched to the COSMOS sample in terms of rest-frame photometry and morphology, as well as their selection and analysis. The offset between the average surface brightness of face-on and edge-on disks in the redshifted SDSS sample matches that predicted by measurements of the optical depth of galactic disks in the nearby universe. In contrast, large disks at z~0.7 have a virtually flat surface brightness-inclination relation, suggesting that they are more opaque than their local counterparts. This could be explained by either an increased amount of optically thick material in disks at higher redshift, or a different spatial distribution of the dust.

preprint2010arXiv

The VLA-COSMOS Perspective on the IR-Radio Relation. I. New Constraints on Selection Biases and the Non-Evolution of the IR/Radio Properties of Star Forming and AGN Galaxies at Intermediate and High Redshift

VLA 1.4 GHz (rms noise ~0.012 mJy) and MIPS 24 and 70 micron (rms noise ~0.02 and ~1.7 mJy, respectively) observations covering the 2 square degree COSMOS field are combined with an extensive multi-wavelength data set to study the evolution of the IR-radio relation at intermediate and high redshift. With ~4500 sources -- of which ~30% have spectroscopic redshifts -- the current sample is significantly larger than previous ones used for the same purpose. Both monochromatic IR/radio flux ratios (q24 & q70), as well as the ratio of the total IR and the 1.4 GHz luminosity (qTIR) are used as indicators for the IR/radio properties of star forming galaxies and AGN. Using a sample jointly selected at IR and radio wavelengths in order to reduce selection biases, we provide firm support for previous findings that the IR-radio relation remains unchanged out to at least z~1.4. Moreover, based on data from ~150 objects we also find that the local relation likely still holds at 2.5<z<5. At redshift z<1.4 we observe that radio-quiet AGN populate the locus of the IR-radio relation in similar numbers as star forming sources. In our analysis we employ the methods of survival analysis in order to ensure a statistically sound treatment of flux limits arising from non-detections. We determine the observed shift in average IR/radio properties of IR- and radio- selected populations and show that it can reconcile apparently discrepant measurements presented in the literature. Finally, we also investigate variations of the IR/radio ratio with IR and radio luminosity and find that it hardly varies with IR luminosity but is a decreasing function of radio luminosity.

preprint2010arXiv

The XMM-Newton Wide-Field Survey in the COSMOS field (XMM-COSMOS): demography and multiwavelength properties of obscured and unobscured luminous AGN

We report the final optical identifications of the medium-depth (~60 ksec), contiguous (2 deg^2) XMM-Newton survey of the COSMOS field. XMM-Newton has detected ~800 X-ray sources down to limiting fluxes of ~5x10^{-16}, ~3x10^{-15}, and ~7x10^{-15} erg/cm2/s in the 0.5-2 keV, 2-10 keV and 5-10 keV bands, respectively. The work is complemented by an extensive collection of multi-wavelength data from 24 micron to UV, available from the COSMOS survey, for each of the X-ray sources, including spectroscopic redshifts for ~50% of the sample, and high-quality photometric redshifts for the rest. The XMM and multiwavelength flux limits are well matched: 1760 (98%) of the X-ray sources have optical counterparts, 1711 (~95%) have IRAC counterparts, and 1394 (~78%) have MIPS 24micron detections. Thanks to the redshift completeness (almost 100%) we were able to constrain the high-luminosity tail of the X-ray luminosity function confirming that the peak of the number density of logL_X>44.5 AGN is at z~2. Spectroscopically-identified obscured and unobscured AGN, as well as normal and starforming galaxies, present well-defined optical and infrared properties. We devised a robust method to identify a sample of ~150 high redshift (z>1), obscured AGN candidates for which optical spectroscopy is not available. We were able to determine that the fraction of the obscured AGN population at the highest (L_X>10^{44} erg s^{-1}) X-ray luminosity is ~15-30% when selection effects are taken into account, providing an important observational constraint for X-ray background synthesis. We studied in detail the optical spectrum and the overall spectral energy distribution of a prototypical Type 2 QSO, caught in a stage transitioning from being starburst dominated to AGN dominated, which was possible to isolate only thanks to the combination of X-ray and infrared observations.

preprint2009arXiv

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

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

preprint2009arXiv

A Multiwavelength Study of a Sample of 70 micron Selected Galaxies in the COSMOS Field I: Spectral Energy Distributions and Luminosities

We present a large robust sample of 1503 reliable and unconfused 70microm selected sources from the multiwavelength data set of the Cosmic Evolution Survey (COSMOS). Using the Spitzer IRAC and MIPS photometry, we estimate the total infrared luminosity, L_IR (8--1000 microns), by finding the best fit template from several different template libraries. The long wavelength 70 and 160 micron data allow us to obtain a reliable estimate of L_IR, accurate to within 0.2 and 0.05 dex, respectively. The 70 micron data point enables a significant improvement over the luminosity estimates possible with only a 24 micron detection. The full sample spans a wide range in L_IR, L_IR ~ 10^8-10^14 L_sun, with a median luminosity of 10^11.4 L_sun. We identify a total of 687 luminous, 303 ultraluminous, and 31 hyperluminous infrared galaxies (LIRGs, ULIRGs, and HyLIRGs) over the redshift range 0.01<z<3.5 with a median redshift of 0.5. Presented here are the full spectral energy distributions for each of the sources compiled from the extensive multiwavelength data set from the ultraviolet (UV) to the far-infrared (FIR). Using SED fits we find possible evidence for a subset of cooler ultraluminous objects than observed locally. However, until direct observations at longer wavelengths are obtained, the peak of emission and the dust temperature cannot be well constrained. We use these SEDs, along with the deep radio and X-ray coverage of the field, to identify a large sample of candidate active galactic nuclei (AGN). We find that the fraction of AGN increases strongly with L_IR, as it does in the local universe, and that nearly 70% of ULIRGs and all HyLIRGs likely host a powerful AGN.

preprint2009arXiv

Cosmic evolution of radio selected active galactic nuclei in the COSMOS field

We explore the cosmic evolution of radio AGN with low radio powers (L_1.4GHz < 5\times10^25 W/Hz) out to z=1.3 using to-date the largest sample of ~600 low luminosity radio AGN at intermediate redshift drawn from the VLA-COSMOS survey. We derive the radio luminosity function for these AGN, and its evolution with cosmic time assuming two extreme cases: i) pure luminosity and ii) pure density evolution. The former and latter yield L_*\propto(1+z)^(0.8+/-0.1), and Phi*\propto (1+z)^(1.1+/-0.1), respectively, both implying a fairly modest change in properties of low radio-power AGN since z=1.3. We show that this is in stark contrast with the evolution of powerful (L_1.4GHz > 5\times10^25 W/Hz) radio AGN over the same cosmic time interval, constrained using the 3CRR, 6CE, and 7CRS radio surveys by Willott et al. (2001). We demonstrate that this can be explained through differences in black hole fueling and triggering mechanisms, and a dichotomy in host galaxy properties of weak and powerful AGN. Our findings suggest that high and low radio-power AGN activity is triggered in different stages during the formation of massive red galaxies. We show that weak radio AGN occur in the most massive galaxies already at z~1, and they may significantly contribute to the heating of their surrounding medium and thus inhibit gas accretion onto their host galaxies, as recently suggested for the `radio mode' in cosmological models.

preprint2009arXiv

Molecular Gas in NUclei of GAlaxies (NUGA) XIII. The interacting Seyfert 2/LINER galaxy NGC 5953

We present $^{12}$CO(1--0) and $^{12}$CO(2--1) maps of the interacting Seyfert 2/LINER galaxy NGC 5953 obtained with the IRAM interferometer at resolutions of 2\farcs1 $\times$ 1\farcs4 and 1\farcs1 $\times$ 0\farcs7, respectively. The CO emission is distributed over a disk of diameter $\sim$16\arcsec ($\sim$2.2 kpc), within which are several, randomly distributed peaks. The strongest peak does not coincide with the nucleus, but is instead offset from the center, $\sim2-3^{\prime\prime}$ ($\sim$340 pc) toward the west/southwest. The kinematics of the molecular component are quite regular, as is typical of a rotating disk. We also compared the $^{12}$CO distribution of NGC 5953 with observations at other wavelengths in order to study correlations between different tracers of the interstellar medium. Using NIR images, we computed the gravity torques exerted by the stellar potential on the gas. The torques are predominantly positive in both $^{12}$CO(1--0) and $^{12}$CO(2--1), suggesting that gas is not flowing into the center, and less than 5% of the gas angular momentum is exchanged in each rotation. This comes from the regular and almost axisymmetric total mass and gas distributions in the center of the galaxy. In NGC 5953, the AGN is apparently not being actively fueled in the current epoch.

preprint2009arXiv

Molecular Gas in NUclei of GAlaxies (NUGA): XII. The head-on collision in NGC1961

We present high-resolution CO(1-0) and CO(2-1) maps of the LINER 2 galaxy NGC1961. This galaxy is unusual among late-type (Sc) disk galaxies in having a very large radial extent and inferred dynamical mass. We propose a head-on collision scenario to explain the perturbed morphology of this galaxy- both the off-centered rings and the inflated radius. This scenario is supported by the detection of a steep velocity gradient in the CO(1-0) map at the position of a southwest peak in radio continuum and near-infrared emission. This peak would represent the remnant of the disrupting companion. We use numerical models to demonstrate the plausibility of the scenario. While ram pressure stripping could in principle be important for shocking the atomic gas and produce the striking head-tail morphology, the non detection of this small galaxy group in X-ray emission suggests that any hot intragroup medium has too low a density. The collision with the companion of mass ratio 1:4 produces a lopsided wave, which accounts for the sharp boundary observed in the atomic gas on the southern side. We argue that NGC1961 has a strongly warped disk, which gives the false impression of a nearly face-on system; the main disk is actually more edge-on, and this error in the true inclination has led to the surprisingly high dynamical mass for a morphologically late-type galaxy. In addition, the outward propagating ring artificially enlarges the disk. The collision de-stabilizes the inner disk and can provide gas inflow to the active nucleus.

preprint2009arXiv

Photometric properties of Ly alpha emitters at z=4.86 in the COSMOS 2 square degree field

We present results of a survey for Ly alpha emitters at z=4.86 based on optical narrowband (lambda_c=7126 angstrom, Delta lambda=73 angstrom) and broadband (B, V, r', i', and z') observations of the Cosmic Evolution Survey (COSMOS) field using Suprime-Cam on the Subaru Telescope. We find 79 LAE candidates at z=4.86 over a contiguous survey area of 1.83 deg^2, down to the Ly alpha line flux of 1.47 x 10^-17 ergs s^-1 cm^-2. We obtain the Ly alpha luminosity function with a best-fit Schechter parameters of log L^*=42.9^+0.5_-0.3 ergs s^-1 and phi^* = 1.2^+8.0_-1.1 x 10^-4 Mpc^-3 for alpha=-1.5 (fixed). The two-point correlation function for our LAE sample is xi(r) = (r/4.4^+5.7_-2.9 Mpc)^-1.90+/-0.22. In order to investigate the field-to-field variations of the properties of Ly alpha emitters, we divide the survey area into nine tiles of 0.5^circ x 0.5^circ each. We find that the number density varies with a factor of ~ 2 from field to field with high statistical significance. However, we find no significant field-to-field variance when we divide the field into four tiles with 0.7^circ x 0.7^circ each. We conclude that at least 0.5 deg^2 survey area is required to derive averaged properties of LAEs at z~5, and our survey field is wide enough to overcome the cosmic variance.

preprint2009arXiv

Properties and environment of Radio Emitting Galaxies in the VLA-zCOSMOS survey

We investigate the properties and the environment of radio sources with optical counterpart from the combined VLA-COSMOS and zCOSMOS samples. The advantage of this sample is the availability of optical spectroscopic information, high quality redshifts, and accurate density determination. By comparing the star formation rates estimated from the optical spectral energy distribution with those based on the radio luminosity, we divide the radio sources in three families, passive AGN, non-passive AGN and star forming galaxies. These families occupy specific regions of the 8.0-4.5 $μ$m infrared color--specific star formation plane, from which we extract the corresponding control samples. Only the passive AGN have a significantly different environment distribution from their control sample. The fraction of radio-loud passive AGN increases from ~2% in underdense regions to ~15% for overdensities (1+delta) greater than 10. This trend is also present as a function of richness of the groups hosting the radio sources. Passive AGN in overdensities tend to have higher radio luminosities than those in lower density environments. Since the black hole mass distribution is similar in both environments, we speculate that, for low radio luminosities, the radio emission is controlled (through fuel disponibility or confinement of radio jet by local gas pressure) by the interstellar medium of the host galaxy, while in other cases it is determined by the structure (group or cluster) in which the galaxy resides.

preprint2009arXiv

Spitzer 70 and 160-micron Observations of the COSMOS Field

We present Spitzer 70 and 160 micron observations of the COSMOS Spitzer survey (S-COSMOS). The data processing techniques are discussed for the publicly released products consisting of images and source catalogs. We present accurate 70 and 160 micron source counts of the COSMOS field and find reasonable agreement with measurements in other fields and with model predictions. The previously reported counts for GOODS-North and the extragalactic First Look Survey are updated with the latest calibration, and counts are measured based on the large area SWIRE survey to constrain the bright source counts. We measure an extragalactic confusion noise level of sigma_c = 9.4+/-3.3 mJy (q=5) for the MIPS 160-micron band based on the deep S-COSMOS data and report an updated confusion noise level of sigma_c = 0.35+/-0.15 mJy (q=5) for the MIPS 70-micron band.

preprint2009arXiv

Star formation and dust obscuration at z~2: galaxies at the dawn of downsizing

We present first results of a study aimed to constrain the star formation rate and dust content of galaxies at z~2. We use a sample of BzK-selected star-forming galaxies, drawn from the COSMOS survey, to perform a stacking analysis of their 1.4 GHz radio continuum as a function of different stellar population properties, after removing AGN contaminants from the sample. Dust unbiased star formation rates are derived from radio fluxes assuming the local radio-IR correlation. The main results of this work are: i) specific star formation rates are constant over about 1 dex in stellar mass and up to the highest stellar mass probed; ii) the dust attenuation is a strong function of galaxy stellar mass with more massive galaxies being more obscured than lower mass objects; iii) a single value of the UV extinction applied to all galaxies would lead to grossly underestimate the SFR in massive galaxies; iv) correcting the observed UV luminosities for dust attenuation based on the Calzetti recipe provide results in very good agreement with the radio derived ones; v) the mean specific star formation rate of our sample steadily decreases by a factor of ~4 with decreasing redshift from z=2.3 to 1.4 and a factor of ~40 down the local Universe. These empirical SFRs would cause galaxies to dramatically overgrow in mass if maintained all the way to low redshifts, we suggest that this does not happen because star formation is progressively quenched, likely starting from the most massive galaxies.

preprint2009arXiv

The Detection of Anomalous Dust Emission in the Nearby Galaxy NGC 6946

We report on the Ka-band (26-40 GHz) emission properties for 10 star-forming regions in the nearby galaxy NGC 6946. From a radio spectral decomposition, we find that the 33 GHz flux densities are typically dominated by thermal (free-free) radiation. However, we also detect excess Ka-band emission for an outer-disk star-forming region relative to what is expected given existing radio, submillimeter, and infrared data. Among the 10 targeted regions, measurable excess emission at 33 GHz is detected for half of them, but in only one region is the excess found to be statistically significant ($\approx7σ$). We interpret this as the first likely detection of so called `anomalous' dust emission outside of the Milky Way. We find that models explaining this feature as the result of dipole emission from rapidly rotating ultrasmall grains are able to reproduce the observations for reasonable interstellar medium conditions. While these results suggest that the use of Ka-band data as a measure of star formation activity in external galaxies may be complicated by the presence of anomalous dust, it is unclear how significant a factor this will be for globally integrated measurements as the excess emission accounts for $\la$10% of the total Ka-band flux density from all 10 regions.

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 X-ray to optical-UV luminosity ratio of X-ray selected Type 1 AGN in XMM-COSMOS

We present a study of the X-ray to optical properties of a sample of 545 X-ray selected Type 1 AGN, from the XMM-COSMOS survey, over a wide range of redshifts ($0.04<\z<4.25$) and X-ray luminosities ($40.6 \leq \Log \Lhard \leq 45.3$). About 60% of them are spectroscopically identified Type 1 AGN, while the others have a reliable photometric redshift and are classified as Type 1 AGN on the basis of their multi-band Spectral Energy Distributions. We discuss the relationship between UV and X-ray luminosity, as parameterized by the $\alphaox$ spectral slope, and its dependence on redshift and luminosity. We compare our findings with previous investigations of optically selected broad-line AGN (mostly from SDSS). A highly significant correlation between $\alphaox$ and $\lo$ is found, in agreement with previous investigations of optically selected samples. We calculate bolometric corrections, $\kbol$, for the whole sample using hard X-ray luminosities ($\Lhard$), and the Eddington ratios for a subsample of 150 objects for which black hole mass estimates are available. We confirm the trend of increasing bolometric correction with increasing Eddington ratio as proposed in previous works. A tight correlation is found between $\alphaox$ and $\kbol$, which can be used to estimate accurate bolometric corrections using only optical and X-ray data. We find a significant correlation between $\alphaox$ and Eddington ratio, in which $\alphaox$ increases for increasing Eddington ratios.

preprint2008arXiv

AGN Host Galaxy Morphologies in COSMOS

We use HST/ACS images and a photometric catalog of the COSMOS field to analyze morphologies of the host galaxies of approximately 400 AGN candidates at redshifts 0.3 < z < 1.0. We compare the AGN hosts with a sample of non-active galaxies drawn from the COSMOS field to match the magnitude and redshift distribution of the AGN hosts. We perform 2-D surface brightness modeling with GALFIT to yield host galaxy and nuclear point source magnitudes. X-ray selected AGN host galaxy morphologies span a substantial range that peaks between those of early-type, bulge-dominated and late-type, disk-dominated systems. We also measure the asymmetry and concentration of the host galaxies. Unaccounted for, the nuclear point source can significantly bias results of these measured structural parameters, so we subtract the best-fit point source component to obtain images of the underlying host galaxies. Our concentration measurements reinforce the findings of our 2-D morphology fits, placing X-ray AGN hosts between early- and late-type inactive galaxies. AGN host asymmetry distributions are consistent with those of control galaxies. Combined with a lack of excess companion galaxies around AGN, the asymmetry distributions indicate that strong interactions are no more prevalent among AGN than normal galaxies. In light of recent work, these results suggest that the host galaxies of AGN at these X-ray luminosities may be in a transition from disk-dominated to bulge-dominated, but that this transition is not typically triggered by major mergers.

preprint2008arXiv

Chasing highly obscured QSOs in the COSMOS field

(abridged) We take advantage of the deep Chandra and Spitzer coverage of a large area (more than 10 times the area covered by the Chandra deep fields, CDFs in the COSMOS field, to extend the search of highly obscured, Compton-thick active nuclei to higher luminosity. These sources have low surface density and large samples can be provided only through large area surveys, like the COSMOS survey. We analyze the X-ray properties of COSMOS MIPS sources with 24$μ$m fluxes higher than 550$μ$Jy. For the MIPS sources not directly detected in the Chandra images we produce stacked images in soft and hard X-rays bands. To estimate the fraction of Compton-thick AGN in the MIPS source population we compare the observed stacked count rates and hardness ratios to those predicted by detailed Monte Carlo simulations including both obscured AGN and star-forming galaxies. The density of lower luminosity Compton-thick AGN (logL(2-10keV)=43.5-44) at z=0.7--1.2 is $(3.7\pm1.1) \times10^{-5}$ Mpc$^{-3}$, corresponding to $\sim67%$ of that of X-ray selected AGN. The comparison between the fraction of infrared selected, Compton thick AGN to the X-ray selected, unobscured and moderately obscured AGN at high and low luminosity suggests that Compton-thick AGN follow a luminosity dependence similar to that discovered for Compton-thin AGN, becoming relatively rarer at high luminosities. We estimate that the fraction of AGN (unobscured, moderately obscured and Compton thick) to the total MIPS source population is $49\pm10%$, a value significantly higher than that previously estimated at similar 24$μ$m fluxes. We discuss how our findings can constrain AGN feedback models.

preprint2007arXiv

COSMOS : Hubble Space Telescope Observations

The Cosmic Evolution Survey (COSMOS) was initiated with an extensive allocation (590 orbits in Cycles 12-13) using the Hubble Space Telescope (HST) for high resolution imaging. Here we review the characteristics of the HST imaging with the Advanced Camera for Surveys (ACS) and parallel observations with NICMOS and WFPC2. A square field (1.8$\sq$°) has been imaged with single-orbit ACS I-F814W exposures with 50% completeness for sources 0.5\arcsec in diameter at I$_{AB} $ = 26.0 mag. The ACS imaging is a key part of the COSMOS survey, providing very high sensitivity and high resolution (0.09\arcsec FWHM, 0.05\arcsec pixels) imaging and detecting 1.2 million objects to a limiting magnitude of 26.5 (AB). These images yield resolved morphologies for several hundred thousand galaxies. The small HST PSF also provides greatly enhanced sensitivity for weak lensing investigations of the dark matter distribution.

preprint2007arXiv

NUclei of GAlaxies: V. Radio emission in 7 NUGA sources

We present high angular resolution radio snap-shot observations of seven nearby low-luminosity active galaxies (LLAGN) from the NUclei of GAlaxies (NUGA) survey. The observations were conducted with MERLIN and EVN/VLBI at 18cm and 6cm. At all observed angular resolutions and frequencies, we find indications for extended emission in about ~40% of the sources, consistent with the decrease of flux with increasing angular resolution. The extended components resemble jet emission in a majority of cases, consistent with the optically thin synchrotron emission implied by their steep spectra. We consider the compact 6cm EVN/VLBI radio emission of our sources in the context of the ``fundamental plane'' that previous LLAGN studies identified within the three-dimensional parameter space of radio luminosity, X-ray luminosity, and black hole mass. We demonstrate, using NGC7217 and NGC1068 as particular examples, that high-resolution, multi-epoch radio observations offer useful information about the origin of offsets from the fundamental plane.

preprint2007arXiv

The First Release COSMOS Optical and Near-IR Data and Catalog

We present imaging data and photometry for the COSMOS survey in 15 photometric bands between 0.3um and 2.4um. These include data taken on the Subaru 8.3m telescope, the KPNO and CTIO 4m telescopes, and the CFHT 3.6m telescope. Special techniques are used to ensure that the relative photometric calibration is better than 1% across the field of view. The absolute photometric accuracy from standard star measurements is found to be 6%. The absolute calibration is corrected using galaxy spectra, providing colors accurate to 2% or better. Stellar and galaxy colors and counts agree well with the expected values. Finally, as the first step in the scientific analysis of these data we construct panchromatic number counts which confirm that both the geometry of the universe and the galaxy population are evolving.

preprint2007arXiv

The Spitzer Legacy Survey of the HST-ACS 2 sq. deg. COSMOS Field: survey strategy and first analysis

The Spitzer-COSMOS survey (S-COSMOS) is a Legacy program (Cycles 2+3) designed to carry out a uniform deep survey of the full 2 sq deg COSMOS field in all seven Spitzer bands (3.6, 4.5, 5.6, 8.0, 24.0, 70.0, 160.0 u). This paper describes the survey parameters, mapping strategy, data reduction procedures, achieved sensitivities to date, and the complete data set for future reference. We show that the observed infrared backgrounds in the S-COSMOS field are within 10% of the predicted background levels. The fluctuations in the background at 24u have been measured and do not show any significant contribution from cirrus, as expected. In addition, we report on the number of asteroid detections in the low galactic latitude COSMOS field. We use the Cycle 2 S-COSMOS data to determine preliminary number counts, and compare our results with those from previous Spitzer Legacy surveys (e.g. SWIRE, GOODS). The results from this "first analysis" confirm that the S-COSMOS survey will have sufficient sensitivity with IRAC to detect ~ L* disks and spheroids out to z ~ 3, and with MIPS to detect ultraluminous starbursts and AGN out to z ~3 at 24u and out to z ~1.5-2 at 70u and 160u.

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.

preprint2004arXiv

Using Bars As Signposts of Galaxy Evolution at High and Low Redshifts

An analysis of the NICMOS Deep Field shows that there is no evidence of a decline in the bar fraction beyond z~0.7, as previously claimed; both bandshifting and spatial resolution must be taken into account when evaluating the evolution of the bar fraction. Two main caveats of this study were a lack of a proper comparison sample at low redshifts and a larger number of galaxies at high redshifts. We address these caveats using two new studies. For a proper local sample, we have analyzed 134 spirals in the near-infrared using 2MASS (main results presented by Menendez-Delmestre in this volume) which serves as an ideal anchor for the low-redshift Universe. In addition to measuring the mean bar properties, we find that bar size is correlated with galaxy size and brightness, but the bar ellipticity is not correlated with these galaxy properties. The bar length is not correlated with the bar ellipticity. For larger high redshift samples we analyze the bar fraction from the 2-square degree COSMOS ACS survey. We find that the bar fraction at z~0.7 is ~50%, consistent with our earlier finding of no decline in bar fraction at high redshifts.

preprint1998arXiv

Molecular Gas and Star Formation in the Host Galaxy of the QSO I Zw 1

We have investigated the ISM of the I Zw 1 QSO host galaxy with Plateau de Bure mm-interferometry and high angular resolution near-infrared imaging spectroscopy. We have detected a circumnuclear gas ring of diameter ~ 1.5'' (1.8 kpc) in its millimetric CO line emission and have mapped the disk and the spiral arms of the host galaxy in the ^{12}CO(1-0) line at 115 GHz as well as in the H (1.65 um) and K (2.2 um) band . Combining our new mm- and NIR-data with available estimates of the radio- and far-infrared contributions to the nuclear emission, we find strong evidence for a nuclear starburst ring. A comparison to other sources with nuclear activity indicates that these rings may be a common phenomenon and contribute a large fraction of the central luminosity.