Source author record

C. Kramer

C. Kramer 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

110works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

110 published item(s)

preprint2026arXiv

Resolved ISM properties and scaling relations in the barred galaxy NGC 3627: constraints from NIKA2 observations

We investigate the interplay between star formation, interstellar medium (ISM) components, and dust properties in NGC 3627 using new NIKA2 1.15 and 2 mm observations from the IMEGIN Large Program. Our goal is to analyze dust and radio emission, decompose contributions in the millimeter-centimeter regime, and explore ISM properties within the galaxy. We perform spectral energy distribution fitting, at both global and spatial scales, using the THEMIS dust model within the HerBIE code, applied to data from 3.4 $μ$m to 6 cm. We decompose emission into dust, free-free, and synchrotron components, and examine correlations with gas surface density and star formation activity. Additionally, we analyze the small dust grain fraction and its variation across the galaxy. We find $\sim$10% radio emission at 2 mm, peaking at 18% in the southern bar-end, which hosts the highest star formation activity. However, an isolated star-forming region beyond this bar-end is the most efficient, as indicated by its elevated dust production efficiency and effective yield, predicted by our simplistic dust evolution model. The 160 $μ$m emission shows the strongest correlation with molecular gas, while 1.15 mm better traces the dust mass surface density. Small grains, which make up $\sim$13% of dust mass (2 $\times$ 10$^{7}$ M$_{\odot}$), are depleted in intense radiation fields, with a notable deficit in the southern tidal tail. ISM properties and chemical evolution indicate that dynamical processes, such as bar-driven gas flows and tidal interactions, are crucial in shaping the galactic structure, influencing star formation efficiency, and dust distribution.

preprint2022arXiv

Massive merging cluster PSZ2G091 as seen by the NIKA2 camera

PSZ2 G091.83+26.11 is a galaxy cluster with M500 = 7.43 x 10^14 Msun at z = 0.822 1. This object exhibits a complex morphology with a clear bimodality observed in X-rays. However, it was detected and analysed in the Planck sample as a single, spherical cluster following a universal profile 2. This model can lead to miscalculations of thermodynamical quantities, like the pressure profile. As future multiwavelength cluster experiments will detect more and more objects at high redshifts, it is crucial to quantify this systematic effect. In this work, we use high-resolution observations of the NIKA2 camera3,4,5,6 to integrate the morphological characteristics of the cluster in our modelling. This is achieved by fitting a two-halo model to the SZ image and then by reconstruction of the resulting projected pressure profile. We then compare these results with the spherical assumption.

preprint2021arXiv

[CII] 158$μ$m emission from Orion A. II. Photodissociation region physics

The [CII] 158$μ$m fine-structure line is the dominant cooling line of moderate-density photodissociation regions (PDRs) illuminated by moderately bright far-ultraviolet (FUV) radiation fields. We aim to understand the origin of [CII] emission and its relation to other tracers of gas and dust in PDRs. One focus is a study of the heating efficiency of interstellar gas as traced by the [CII] line to test models of the photoelectric heating of neutral gas by polycyclic aromatic hydrocarbon (PAH) molecules and very small grains. We make use of a one-square-degree map of velocity-resolved [CII] line emission toward the Orion Nebula complex, and split this out into the individual spatial components, the expanding Veil Shell, the surface of OMC4, and the PDRs associated with the compact HII region of M43 and the reflection nebula NGC 1977. We employed Herschel far-infrared photometric images to determine dust properties. Moreover, we compared with Spitzer mid-infrared photometry to trace hot dust and large molecules, and velocity-resolved IRAM 30m CO(2-1) observations of the molecular gas. The [CII] intensity is tightly correlated with PAH emission in the IRAC 8$μ$m band and far-infrared emission from warm dust. The correlation between [CII] and CO(2-1) is very different in the four subregions and is very sensitive to the detailed geometry. Constant-density PDR models are able to reproduce the observed [CII], CO(2-1), and integrated far-infrared (FIR) intensities. We observe strong variations in the photoelectric heating efficiency in the Veil Shell behind the Orion Bar and these variations are seemingly not related to the spectral properties of the PAHs. The [CII] emission from the Orion Nebula complex stems mainly from moderately illuminated PDR surfaces. Future observations with the James Webb Space Telescope can shine light on the PAH properties that may be linked to these variations.

preprint2021arXiv

Crab nebula at 260 GHz with the NIKA2 polarimeter. Implications for the polarization angle calibration of future CMB experiments

The quest for primordial gravitational waves enclosed in the Cosmic Microwave Background (CMB) polarization B-modes signal motivates the development of a new generation of high sensitive experiments (e.g. CMB-S4, LiteBIRD) that would allow them to detect its imprint.Neverthless, this will be only possible by ensuring a high control of the instrumental systematic effects and an accurate absolute calibration of the polarization angle. The Crab nebula is known to be a polarization calibrator on the sky for CMB experiments, already used for the Planck satellite it exhibits a high polarized signal at microwave wavelengths. In this work we present Crab polarization observations obtained at the central frequency of 260 GHz with the NIKA2 instrument and discuss the accuracy needed on such a measurement to improve the constraints on the absolute angle calibration for CMB experiments.

preprint2021arXiv

Dust Emission in Galaxies at Millimeter Wavelengths: Cooling of star forming regions in NGC6946

Interstellar dust plays an important role in the formation of molecular gas and the heating and cooling of the interstellar medium. The spatial distribution of the mm-wavelength dust emission from galaxies is largely unexplored. The NIKA2 Guaranteed Time Project IMEGIN (Interpreting the Millimeter Emission of Galaxies with IRAM and NIKA2) has recently mapped the mm emission in the grand design spiral galaxy NGC6946. By subtracting the contributions from the free-free, synchrotron, and CO line emission, we map the distribution of the pure dust emission at 1:15mm and 2mm. Separating the arm/interarm regions, we find a dominant 2mm emission from interarms indicating the significant role of the general interstellar radiation field in heating the cold dust. Finally, we present maps of the dust mass, temperature, and emissivity index using the Bayesian MCMC modeling of the spectral energy distribution in NGC6946.

preprint2021arXiv

Exploring the millimetre emission in nearby galaxies: analysis of the edge-on galaxy NGC 891

New observations of the edge-on galaxy NGC 891, at 1.15 and 2 mm obtained with the IRAM 30-m telescope and the NIKA2 camera, within the framework of the IMEGIN (Interpreting the Millimetre Emission of Galaxies with IRAM and NIKA2) Large Program, are presented in this work. By using multiwavelength maps (from the mid-IR to the cm wavelengths) we perform SED fitting in order to extract the physical properties of the galaxy on both global and local ($\sim$kpc) scales. For the interpretation of the observations we make use of a state-of-the-art SED fitting code, HerBIE (HiERarchical Bayesian Inference for dust Emission). The observations indicate a galaxy morphology, at mm wavelengths, similar to that of the cold dust emission traced by sub-mm observations and to that of the molecular gas. The contribution of the radio emission at the NIKA2 bands is very small (negligible at 1.15 mm and $\sim10\%$ at 2 mm) while it dominates the total energy budget at longer wavelengths (beyond 5 mm). On local scales, the distribution of the free-free emission resembles that of the dust thermal emission while the distribution of the synchrotron emission shows a deficiency along the major axis of the disc of the galaxy.

preprint2021arXiv

Galactic Star Formation with NIKA2 (GASTON): Evidence of mass accretion onto dense clumps

High-mass stars ($m_* \gtrsim 8 \, M_\odot$) play a crucial role in the evolution of galaxies, and so it is imperative that we understand how they are formed. We have used the New IRAM KIDs Array 2 (NIKA2) camera on the Institut de Radio Astronomie Millimétrique (IRAM) 30-m telescope to conduct high-sensitivity continuum mapping of $\sim2$ deg$^2$ of the Galactic plane (GP) as part of the Galactic Star Formation with NIKA2 (GASTON) large program. We have identified a total of 1467 clumps within our deep 1.15 mm continuum maps and, by using overlapping continuum, molecular line, and maser parallax data, we have determined their distances and physical properties. By placing them upon an approximate evolutionary sequence based upon 8 $μ$m $\textit{Spitzer}$ imaging, we find evidence that the most massive dense clumps accrete material from their surrounding environment during their early evolution, before dispersing as star formation advances, supporting clump-fed models of high-mass star formation.

preprint2021arXiv

Galactic star formation with NIKA2 (GASTON): Filament convergence and its link to star formation

In the past decade filaments have been recognised as a major structural element of the interstellar medium, the densest of these filaments hosting the formation of most stars. In some star-forming molecular clouds converging networks of filaments, also known as hub filament systems, can be found. These hubs are believed to be preferentially associated to massive star formation. As of today, there are no metrics that allow the systematic quantification of a filament network convergence. Here, we used the IRAM 30m NIKA2 observations of the Galactic plane from the GASTON large programme to systematically identify filaments and produce a filament convergence parameter map. We use such a map to show that: i. hub filaments represent a small fraction of the global filament population; ii. hubs host, in proportion, more massive and more luminous compact sources that non-hubs; iii. hub-hosting clumps are more evolved that non-hubs; iv. no discontinuities are observed in the properties of compact sources as a function of convergence parameter. We propose that the rapid global collapse of clumps is responsible for (re)organising filament networks into hubs and, in parallel, enhancing the mass growth of compact sources.

preprint2021arXiv

GASTON: Galactic Star Formation with NIKA2. Evidence for the mass growth of star-forming clumps

Determining the mechanism by which high-mass stars are formed is essential for our understanding of the energy budget and chemical evolution of galaxies. By using the New IRAM KIDs Array 2 (NIKA2) camera on the Institut de Radio Astronomie Millimétrique (IRAM) 30-m telescope, we have conducted high-sensitivity and large-scale mapping of a fraction of the Galactic plane in order to search for signatures of the transition between the high- and low-mass star-forming modes. Here, we present the first results from the Galactic Star Formation with NIKA2 (GASTON) project, a Large Programme at the IRAM 30-m telescope which is mapping $\approx$2 deg$^2$ of the inner Galactic plane (GP), centred on $\ell$=23.9$^\circ$, $b$=0.05$^\circ$, as well as targets in Taurus and Ophiuchus in 1.15 and 2.00 mm continuum wavebands. In this paper we present the first of the GASTON GP data taken, and present initial science results. We conduct an extraction of structures from the 1.15 mm maps using a dendrogram analysis and, by comparison to the compact source catalogues from Herschel survey data, we identify a population of 321 previously-undetected clumps. Approximately 80 per cent of these new clumps are 70 $μ$m-quiet, and may be considered as starless candidates. We find that this new population of clumps are less massive and cooler, on average, than clumps that have already been identified. Further, by classifying the full sample of clumps based upon their infrared-bright fraction - an indicator of evolutionary stage - we find evidence for clump mass growth, supporting models of clump-fed high-mass star formation.

preprint2021arXiv

Mapping the intracluster medium temperature in the era of NIKA2 and MUSTANG-2

We present preliminary results from an on-going program that aims at mapping the intracluster medium (ICM) temperature of high redshift galaxy clusters from the MaDCoWS sample using a joint analysis of shallow X-ray data obtained by $Chandra$ and high angular resolution Sunyaev-Zel'dovich (SZ) observations realized with the NIKA2 and MUSTANG-2 cameras. We also present preliminary results from an on-going Open Time program within the NIKA2 collaboration that aims at mapping the ICM temperature of a galaxy cluster at $z=0.45$ from the resolved detection of the relativistic corrections to the SZ spectrum. These studies demonstrate how high angular resolution SZ observations will play a major role in the coming decade to push the investigation of ICM dynamics and non-gravitational processes to high redshift before the next generation X-ray observatories come into play.

preprint2021arXiv

Multi-probe analysis of the galaxy cluster CL J1226.9+3332: hydrostatic mass and hydrostatic-to-lensing bias

We present a multi-probe analysis of the well-known galaxy cluster CL J1226.9+3332 as a proof of concept for multi-wavelength studies within the framework of the NIKA2 Sunyaev-Zeldovich Large Program (LPSZ). CL J1226.9+3332 is a massive and high redshift (z = 0.888) cluster that has already been observed at several wavelengths. A joint analysis of the thermal SZ (tSZ) effect at millimeter wavelength with the NIKA2 camera and in X-ray with the XMM-Newton satellite permits the reconstruction of the cluster thermodynamical properties and mass assuming hydrostatic equilibrium. We test the robustness of our mass estimates against different definitions of the data analysis transfer function. Using convergence maps reconstructed from the data of the CLASH program we obtain estimates of the lensing mass, which we compare to the estimated hydrostatic mass. This allows us to measure the hydrostatic-to-lensing mass bias and the associated systematic effects related to the NIKA2 measurement. We obtain M500HSE = (7.65 +- 1.03) 1014 Msun and M500lens = (7.35 +- 0.65) 1014 Msun, which implies a HSE-to-lensing bias consistent with 0 within 20 percent.

preprint2021arXiv

Overdensity of SubMillimiter Galaxies in the GJ526 Field mapped with the NIKA2 Camera

Using the NIKA2 dual band millimeter camera installed on the IRAM30m telescope, we have mapped a relatively large field (~70 arcmin^2) in the direction of the star GJ526 to investigate the nature of the sources found with the MAMBO camera at 1.2 mm ten years earlier. We have found that they must be dust-obscured galaxies (SMGs) in the background beyond the star. The new NIKA2 map at 1.15 mm reveals additional sources and, in fact, an overdensity of SMGs predominantly distributed along a filament-like structure in projection on the sky across the whole observed field. We speculate this might be a cosmic filament at high redshift as revealed in cosmological hydrodynamical simulations. Measurement of spectroscopic redshifts of the SMGs in the candidate filament is required now for a definitive confirmation of the nature of the structure.

preprint2021arXiv

Probing the role of magnetic fields in star-forming filaments: NIKA2-Pol commissioning results toward OMC-1

Dust polarization observations are a powerful, practical tool to probe the geometry (and to some extent, the strength) of magnetic fields in star-forming regions. In particular, Planck polarization data have revealed the importance of magnetic fields on large scales in molecular clouds. However, due to insufficient resolution, Planck observations are unable to constrain the B-field geometry on prestellar and protostellar scales. The high angular resolution of 11.7 arcsec provided by NIKA2-Pol 1.15 mm polarimetric imaging, corresponding to $\sim$ 0.02 pc at the distance of the Orion molecular cloud (OMC), makes it possible to advance our understanding of the B-field morphology in star-forming filaments and dense cores (IRAM 30m large program B-FUN). The commissioning of the NIKA2-Pol instrument has led to several challenging issues, in particular, the instrumental polarization or intensity-to-polarization (leakage) effect. In the present paper, we illustrate how this effect can be corrected for, leading to reliable exploitable data in a structured, extended source such as OMC-1. We present a statistical comparison between NIKA2-Pol and SCUBA2-Pol2 results in the OMC-1 region. We also present tentative evidence of local pinching of the B-field lines near Orion-KL, in the form of a new small-scale hourglass pattern, in addition to the larger-scale hourglass already seen by other instruments such as Pol2.

preprint2021arXiv

PSZ2G091:A massive double cluster at z=0.822 observed by the NIKA2 camera

PSZ2 G091.83+26.11 is a massive galaxy cluster with M500 = 7.43 x 10^14 Msun at z = 0.822. This object exhibits a complex morphology with a clear bimodality observed in X-rays. However, it was detected and analysed in the Planck sample as a single, spherical cluster following a universal profile [1]. This model can lead to miscalculations of thermodynamical quantities, like the pressure profile. As future multiwavelength cluster experiments will detect more and more objects at higher redshifts (where we expect the fraction of merging objects to be higher), it is crucial to quantify this systematic effect. In this work, we use high-resolution observations of PSZ2 G091.83+26.11 by the NIKA2 camera to integrate the morphological characteristics of the cluster in our modelling. This is achieved by fitting a two-halo model to the SZ image and then by reconstruction of the resulting projected pressure profile. We then compare these results with the spherical assumption.

preprint2021arXiv

Searching for high-z DSFGs with NIKA2 and NOEMA

As the possible progenitors of passive galaxies at z=2-3, dusty star-forming galaxies (DSFGs) at z>4 provide a unique perspective to study the formation, assembly, and early quenching of massive galaxies in the early Universe. The extreme obscuration in optical-IR makes (sub)mm spectral scans the most universal and unbiased way to confirm/exclude the high-z nature of candidate dusty star-forming galaxies. We present here the status of the NIKA2 Cosmological Legacy Survey (N2CLS), which is the deepest wide-area single-dish survey in the millimeter searching for high-z DSFGs. We also introduce a joint-analysis method to efficiently search for the spectroscopic redshift of high-z DSFGs with noisy spectra and photometric data and present its success in identifying the redshift of DSFGs found in NIKA2 science verification data.

preprint2021arXiv

The LPSZ-CLASH galaxy cluster sample: combining lensing and hydrostatic mass estimates

Starting from the clusters included in the NIKA sample and in the NIKA2 Sunyaev-Zel'dovich Large Program (LPSZ) we have selected a sample of six common objects with the Cluster Lensing And Supernova survey with Hubble (CLASH) lensing data. For the LPSZ clusters we have at our disposal both high-angular resolution observations of the thermal SZ with NIKA and NIKA2 and X-ray observations with XMM-Newton from which hydrostatic mass estimates can be derived. In addition, the CLASH dataset includes lensing convergence maps that can be converted into lensing estimates of the total mass of the cluster. One-dimensional mass profiles are used to derive integrated mass estimates accounting for systematic effects (data processing, modeling, etc.). Two-dimensional analysis of the maps can reveal substructures in the cluster and, therefore, inform us about the dynamical state of each system. Moreover, we are able to study the hydrostatic mass to lensing mass bias, across different morphology and a range of redshift clusters to give more insight on the hydrostatic mass bias. The analysis presented in this proceeding follows the study discussed in Ferragamo et al. 2021.

preprint2021arXiv

The NIKA2 Sunyaev-Zeldovich Large Program

The NIKA2 Guaranteed-Time SZ Large Program (LPSZ) is dedicated to the high-angular resolution SZ mapping of a representative sample of 45 SZ-selected galaxy clusters drawn from the catalogues of the Planck satellite, or of the Atacama Cosmology Telescope. The LPSZ sample spans a mass range from $3$ to $11 \times 10^{14} M_{\odot}$ and a redshift range from $0.5$ to $0.9$, extending to higher redshift and lower mass the previous samples dedicated to the cluster mass calibration and universal properties estimation. The main goals of the LPSZ are the measurement of the average radial profile of the ICM pressure up to $R_{500}$ by combining NIKA2 with Planck or ACT data, and the estimation of the scaling law between the SZ observable and the mass using NIKA2, XMM-Newton and Planck/ACT data. Furthermore, combining LPSZ data with existing or forthcoming public data in lensing, optical/NIR or radio domains, we will build a consistent picture of the cluster physics and further gain knowledge on the mass estimate as a function of the cluster morphology and dynamical state.

preprint2020arXiv

Calibration and Performance of the NIKA2 camera at the IRAM 30-meter Telescope

NIKA2 is a dual-band millimetric continuum camera of 2900 Kinetic Inductance Detectors (KID), operating at $150$ and $260\,\rm{GHz}$, installed at the IRAM 30-meter telescope. We present the performance assessment of NIKA2 after one year of observation using a dedicated point-source calibration method, referred to as the \emph{baseline} method. Using a large data set acquired between January 2017 and February 2018 that span the whole range of observing elevations and atmospheric conditions encountered at the IRAM 30-m telescope, we test the stability of the performance parameters. We report an instantaneous field of view (FOV) of 6.5' in diameter, filled with an average fraction of $84\%$ and $90\%$ of valid detectors at $150$ and $260\,\rm{GHz}$, respectively. The beam pattern is characterized by a FWHM of $17.6'' \pm 0.1''$ and $11.1''\pm 0.2''$, and a beam efficiency of $77\% \pm 2\%$ and $55\% \pm 3\%$ at $150$ and $260\,\rm{GHz}$, respectively. The rms calibration uncertainties are about $3\%$ at $150\,\rm{GHz}$ and $6\%$ at $260\,\rm{GHz}$. The absolute calibration uncertainties are of $5\%$ and the systematic calibration uncertainties evaluated at the IRAM 30-m reference Winter observing conditions are below $1\%$ in both channels. The noise equivalent flux density (NEFD) at $150$ and $260\,\rm{GHz}$ are of $9 \pm 1\, \rm{mJy}\cdot s^{1/2}$ and $30 \pm 3\, \rm{mJy}\cdot s^{1/2}$. This state-of-the-art performance confers NIKA2 with mapping speeds of $1388 \pm 174$ and $111 \pm 11 \,\rm{arcmin}^2\cdot \rm{mJy}^{-2}\cdot \rm{h}^{-1}$ at $150$ and $260\,\rm{GHz}$. With these unique capabilities of fast dual-band mapping at high (better that 18'') angular resolution, NIKA2 is providing an unprecedented view of the millimetre Universe.

preprint2020arXiv

Cluster cosmology with the NIKA2 SZ Large Program

The main limiting factor of cosmological analyses based on thermal Sunyaev-Zel'dovich (SZ) cluster statistics comes from the bias and systematic uncertainties that affect the estimates of the mass of galaxy clusters. High-angular resolution SZ observations at high redshift are needed to study a potential redshift or morphology dependence of both the mean pressure profile and of the mass-observable scaling relation used in SZ cosmological analyses. The NIKA2 camera is a new generation continuum instrument installed at the IRAM 30-m telescope. With a large field of view, a high angular resolution and a high-sensitivity, the NIKA2 camera has unique SZ mapping capabilities. In this paper, we present the NIKA2 SZ large program, aiming at observing a large sample of clusters at redshifts between 0.5 and 0.9, and the characterization of the first cluster oberved with NIKA2.

preprint2020arXiv

Exploiting NIKA2/XMM-Newton imaging synergy for intermediate mass, high-$z$ galaxy clusters within the NIKA2 SZ Large Program

High-resolution mapping of the intra-cluster medium (ICM) up to high redshift and down to low masses is crucial to derive accurate mass estimates of the galaxy cluster and to understand the systematic effects affecting cosmological studies based on galaxy clusters. We present a spatially-resolved Sunyaev-Zel'dovich (SZ)/X-ray analysis of ACT-CL J0215.4+0030, a high redshift ($z=0.865$) galaxy cluster of intermediate mass ($M_{500}\simeq3.5\times10^{14}\;\mathrm{M_\odot}$) observed as part of the ongoing NIKA2 SZ Large Program, a follow up of a representative sample of objects at $0.5 \leqslant z \leqslant 0.9$. In addition to the faintness and small angular size induced by its mass and redshift, the cluster is contaminated by point sources that significantly affect the SZ signal. Therefore, this is an interesting case study for the most challenging sources of the NIKA2 cluster sample. We present the NIKA2 observations of this cluster and the resulting data. We reconstruct the ICM pressure profile by performing a joint analysis of the SZ signal and of the point sources in the NIKA2 150 GHz map. We obtain high-quality estimates of the ICM thermodynamical properties with NIKA2. We compare the pressure profile extracted from the NIKA2 map to the pressure profile obtained from X-ray data only by deprojecting XMM-Newton observations of the cluster. We combine the NIKA2 pressure profile with the X-ray deprojected density to extract detailed information on the ICM. The radial distribution of its thermodynamic properties indicate that the cluster has a disturbed core. The hydrostatic mass of the cluster is to be compatible with estimations from SZ and X-rays scaling relations. We conclude that the NIKA2 SZ large program can deliver quality information on the thermodynamics of the ICM even for one of its faintest clusters, after a careful treatment of the contamination by point sources.

preprint2020arXiv

Gas phase Elemental abundances in Molecular cloudS (GEMS). II. On the quest for the sulphur reservoir in molecular clouds: the $H_{2}S$ case

Sulphur is one of the most abundant elements in the Universe. Surprisingly, sulphuretted molecules are not as abundant as expected in the interstellar medium, and the identity of the main sulphur reservoir is still an open question. Our goal is to investigate the H$_{2}$S chemistry in dark clouds, as this stable molecule is a potential sulphur reservoir. Using millimeter observations of CS, SO, H$_{2}$S, and their isotopologues, we determine the physical conditions and H$_{2}$S abundances along the cores TMC 1-C, TMC 1-CP, and Barnard 1b. The gas-grain model Nautilus is then used to model the sulphur chemistry and explore the impact of photo-desorption and chemical desorption on the H$_2$S abundance. Our model shows that chemical desorption is the main source of gas-phase H$_2$S in dark cores. The measured H$_{2}$S abundance can only be fitted if we assume that the chemical desorption rate decreases by more than a factor of 10 when $n_{\rm H}>2\times10^{4}$. This change in the desorption rate is consistent with the formation of thick H$_2$O and CO ice mantles on grain surfaces. The observed SO and H$_2$S abundances are in good agreement with our predictions adopting an undepleted value of the sulphur abundance. However, the CS abundance is overestimated by a factor of $5-10$. Along the three cores, atomic S is predicted to be the main sulphur reservoir. We conclude that the gaseous H$_2$S abundance is well reproduced, assuming undepleted sulphur abundance and chemical desorption as the main source of H$_2$S. The behavior of the observed H$_{2}$S abundance suggests a changing desorption efficiency, which would probe the snowline in these cores. Our model, however, overestimates the observed gas-phase CS abundance. Given the uncertainty in the sulphur chemistry, our data are consistent with a cosmic elemental S abundance with an uncertainty of a factor of 10.

preprint2020arXiv

Gas-phase Elemental abundances in Molecular cloudS (GEMS) III. Unlocking the CS chemistry: the CS+O reaction

CS is among the most abundant gas-phase S-bearing molecules in cold dark molecular clouds. It is easily observable with several transitions in the millimeter wavelength range, and has been widely used as a tracer of the gas density in the interstellar medium in our Galaxy and external galaxies. Chemical models fail to account for the observed CS abundances when assuming the cosmic value for the elemental abundance of sulfur. The CS+O -> CO + S reaction has been proposed as a relevant CS destruction mechanism at low temperatures, and could explain the discrepancy between models and observations. Its reaction rate has been experimentally measured at temperatures of 150-400 K, but the extrapolation to lower temperatures is doubtful. Here we calculate the CS+O reaction rate at temperatures <150 K which are prevailing in the interstellar medium. We performed ab initio calculations to obtain the three lowest PES of the CS+O system. These PESs are used to study the reaction dynamics, using several methods to eventually calculate the CS+O thermal reaction rates. We compare the results of our theoretical calculations for 150-400 K with those obtained in the laboratory. Our detailed theoretical study on the CS+O reaction, which is in agreement with the experimental data obtained at 150-400 K, demonstrates the reliability of our approach. After a careful analysis at lower temperatures, we find that the rate constant at 10 K is negligible, which is consistent with the extrapolation of experimental data using the Arrhenius expression. We use the updated chemical network to model the sulfur chemistry in TMC1 based on molecular abundances determined from GEMS project observations. In our model, we take into account the expected decrease of the cosmic ray ionization rate along the cloud. The abundance of CS is still overestimated when assuming the cosmic value for the sulfur abundance.

preprint2020arXiv

Large-scale Map of Millimeter-wavelength Hydrogen Radio Recombination Lines around a Young Massive Star Cluster

We report the first map of large-scale (10 pc in length) emission of millimeter-wavelength hydrogen recombination lines (mm-RRLs) toward the giant H II region around the W43-Main young massive star cluster (YMC). Our mm-RRL data come from the IRAM 30 m telescope and are analyzed together with radio continuum and cm-RRL data from the Karl G. Jansky Very Large Array and HCO$^{+}$ 1-0 line emission data from the IRAM 30 m. The mm-RRLs reveal an expanding wind-blown ionized gas shell with an electron density ~70-1500 cm$^{-3}$ driven by the WR/OB cluster, which produces a total Ly$α$ photon flux of 1.5 x 10$^{50}$ s$^{-1}$. This shell is interacting with the dense neutral molecular gas in the W43-Main dense cloud. Combining the high spectral and angular resolution mm-RRL and cm-RRL cubes, we derive the two-dimensional relative distributions of dynamical and pressure broadening of the ionized gas emission and find that the RRL line shapes are dominated by pressure broadening (4-55 km s$^{-1}$) near the YMC and by dynamical broadening (8-36 km s$^{-1}$) near the shell's edge. Ionized gas clumps hosting ultra-compact H II regions found at the edge of the shell suggest that large-scale ionized gas motion triggers the formation of new star generation near the periphery of the shell.

preprint2020arXiv

LEGO II: A 3 mm molecular line study covering 100 pc of one of the most actively star-forming portions within the Milky Way Disc

The current generation of (sub)mm-telescopes has allowed molecular line emission to become a major tool for studying the physical, kinematic, and chemical properties of extragalactic systems, yet exploiting these observations requires a detailed understanding of where emission lines originate within the Milky Way. In this paper, we present 60$^{\prime\prime}$ ($\sim$3pc) resolution observations of many 3mm-band molecular lines across a large map of the W49 massive star-forming region ($\sim$100$\times$100pc at 11kpc), which were taken as part of the "LEGO" IRAM-30m large project. We find that the spatial extent or brightness of the molecular line transitions are not well correlated with their critical densities, highlighting abundance and optical depth must be considered when estimating line emission characteristics. We explore how the total emission and emission efficiency (i.e. line brightness per H$_{2}$ column density) of the line emission vary as a function of molecular hydrogen column density and dust temperature. We find that there is not a single region of this parameter space responsible for the brightest and most efficiently emitting gas for all species. For example, we find that the HCN transition shows high emission efficiency at high column density ($10^{22}$cm$^{-2}$) and moderate temperatures (35K), whilst e.g. N$_2$H$^+$ emits most efficiently towards lower temperatures ($10^{22}$cm$^{-2}$; <20K). We determine $X_{\mathrm{CO} (1-0)} \sim 0.3 \times 10^{20} \mathrm{cm^{-2}(Kkms^{-1})^{-1}}$, and $α_{\mathrm{HCN} (1-0)} \sim 30\mathrm{M_\odot(Kkms^{-1}pc^2)^{-1}}$, which both differ significantly from the commonly adopted values. In all, these results suggest caution should be taken when interpreting molecular line emission.

preprint2020arXiv

Molecular globules in the Veil bubble of Orion. IRAM 30m 12CO, 13CO, and C18O 2-1 expanded maps of Orion A

Strong winds and ultraviolet (UV) radiation from O-type stars disrupt and ionize their molecular core birthplaces, sweeping up material into parsec-size shells. Owing to dissociation by starlight, the thinnest shells are expected to host low molecular abundances and therefore little star formation. Here, we expand previous maps taken with the IRAM 30m telescope and present square-degree 12CO and 13CO (J=2-1) maps of the wind-driven "Veil bubble'' that surrounds the Trapezium cluster and its natal Orion molecular core (OMC). Although widespread and extended CO emission is largely absent from the Veil, we show that several CO "globules'' exist and are embedded in the [CII]158um-bright shell that confines the bubble. This includes the first detection of quiescent CO at negative LSR velocities in Orion. Given the harsh UV irradiation conditions in this translucent material, the detection of CO globules is surprising. These globules are small (R=7,100 AU), not massive (M=0.3M_Sun), and are moderately dense: n_ H=4x10^4 cm^-3 (median values). They are confined by the external pressure of the shell, P_ext/k~10^7 cm^-3 K, and are likely magnetically supported. They are either transient objects formed by instabilities or have detached from pre-existing molecular structures, sculpted by the passing shock associated with the expanding shell and by UV radiation from the Trapezium. Some represent the first stages in the formation of small pillars, others of isolated small globules. Although their masses do not suggest they will form stars, one globule matches the position of a known YSO. The lack of extended CO in the "Veil shell'' demonstrates that feedback from massive stars expels, agitates, and reprocesses most of the disrupted molecular cloud gas, thereby limiting the star-formation rate in the region. The presence of globules is a result of this feedback.

preprint2019arXiv

A low-mass galaxy cluster as a test-case study for the NIKA2 SZ Large Program

High-resolution mapping of the hot gas in galaxy clusters is a key tool for cluster-based cosmological analyses. Taking advantage of the NIKA2 millimeter camera operated at the IRAM 30-m telescope, the NIKA2 SZ Large Program seeks to get a high-resolution follow-up of 45 galaxy clusters covering a wide mass range at high redshift in order to re-calibrate some of the tools needed for the cosmological exploitation of SZ surveys. We present the second cluster analysis of this program, targeting one of the faintest sources of the sample in order to tackle the difficulties in data reduction for such faint, low-SNR clusters. In this study, the main challenge is the precise estimation of the contamination by sub-millimetric point sources, which greatly affects the tSZ map of the cluster. We account for this contamination by performing a joint fit of the SZ signal and of the flux density of the compact sources. A prior knowledge of these fluxes is given by the adjustment of the SED of each source using data from both NIKA2 and the \textit{Herschel} satellite. The first results are very promising and demonstrate the possibility to estimate thermodynamic properties with NIKA2, even in a compact cluster heavily contaminated by point sources.

preprint2019arXiv

Debris disks around stars in the NIKA2 era

The new NIKA2 camera at the IRAM 30m radiotelescope was used to observe three known debris disks in order to constrain the SED of their dust emission in the millimeter wavelength domain. We have found that the spectral index between the two NIKA2 bands (1mm and 2mm) is consistent with the Rayleigh-Jeans regime (lambda^{-2}), unlike the steeper spectra (lambda^{-3}) measured in the submillimeter-wavelength domain for two of the three disks $-$ around the stars Vega and HD107146. We provide a succesful proof of concept to model this spectral inversion in using two populations of dust grains, those smaller and those larger than a grain radius a0 of 0.5mm. This is obtained in breaking the slope of the size distribution and the functional form of the absorption coefficient of the standard model at a0. The third disk - around the star HR8799 - does not exhibit this spectral inversion but is also the youngest.

preprint2019arXiv

GASTON: Galactic Star Formation with NIKA2: A new population of cold massive sources discovered

Understanding where and when the mass of stars is determined is one of the fundamental, mostly unsolved, questions in astronomy. Here, we present the first results of GASTON, the Galactic Star Formation with NIKA2 large programme on the IRAM 30m telescope, that aims to identify new populations of low-brightness sources to tackle the question of stellar mass determination across all masses. In this paper, we focus on the high-mass star formation part of the project, for which we map a $\sim2$ deg$^2$ region of the Galactic plane around $l=24^\circ$ in both 1.2 mm and 2.0 mm continuum. Half-way through the project, we reach a sensitivity of 3.7 mJy/beam at 1.2mm. Even though larger than our target sensitivity of 2 mJy, the current sensitivity already allows the identification of a new population of cold, compact sources that remained undetected in any (sub-)mm Galactic plane survey so far. In fact, about 25% of the $\sim 1600$ compact sources identified in the 1.2 mm GASTON image are new detections. We present a preliminary analysis of the physical properties of the GASTON sources as a function of their evolutionary stage, arguing for a potential evolution of the mass distribution of these sources with time.

preprint2019arXiv

Mapping the gas thermodynamic properties of the massive cluster merger MOO J1142$+$1527 at z = 1.2

We present the results of the analysis of the very massive cluster MOO J1142$+$1527 at a redshift $z = 1.2$ based on high angular resolution NIKA2 Sunyaev-Zel'dovich (SZ) and $Chandra$ X-ray data. This multi-wavelength analysis enables us to estimate the shape of the temperature profile with unprecedented precision at this redshift and to obtain a map of the gas entropy distribution averaged along the line of sight. The comparison between the cluster morphological properties observed in the NIKA2 and $Chandra$ maps together with the analysis of the entropy map allows us to conclude that MOO J1142$+$1527 is an on-going merger hosting a cool-core at the position of the X-ray peak. This work demonstrates how the addition of spatially-resolved SZ observations to low signal-to-noise X-ray data can bring valuable insights on the intracluster medium thermodynamic properties at $z>1$.

preprint2019arXiv

NIKA: a mm camera for Sunyaev-Zel'dovich science in clusters of galaxies

Clusters of galaxies, the largest bound objects in the Universe, constitute a cosmological probe of choice, which is sensitive to both dark matter and dark energy. Within this framework, the Sunyaev-Zel'dovich (SZ) effect has opened a new window for the detection of clusters of galaxies and for the characterization of their physical properties such as mass, pressure and temperature. NIKA, a KID-based dual band camera installed at the IRAM 30-m telescope, was particularly well adapted in terms of frequency, angular resolution, field-of-view and sensitivity, for the mapping of the thermal and kinetic SZ effect in high-redshift clusters. In this paper, we present the NIKA cluster sample and a review of the main results obtained via the measurement of the SZ effect on those clusters: reconstruction of the cluster radial pressure profile, mass, temperature and velocity.

preprint2019arXiv

NIKA2 mapping and cross-instrument SED extraction of extended sources with Scanamorphos

The steps taken to tailor to NIKA2 observations the Scanamorphos algorithm (initially developed to subtract low-frequency noise from Herschel on-the-fly observations) are described, focussing on the consequences of the different instrument architecture and observation strategy. The method, making the most extensive use of the redundancy built in the multi-scan coverage with large arrays of a given region of the sky, is applicable to extended sources, while the pipeline is so far optimized for compact sources. An example of application is given. A related tool to build consistent broadband SEDs from 60 microns to 2 mm, combining Herschel and NIKA2 data, has also been developed. Its main task is to process the data least affected by low-frequency noise and coverage limitations (i.e. the Herschel data) through the same transfer function as the NIKA2 data, simulating the same scan geometry and applying the same noise and atmospheric signal as extracted from the 1 mm and 2 mm data.

preprint2019arXiv

Observing with NIKA2Pol from the IRAM 30m telescope. Early results on the commissioning phase

The NIKA2 polarization channel at 260 GHz (1.15 mm) has been proposed primarily to observe galactic star-forming regions and probe the critical scales between 0.01-0.05 pc at which magnetic field lines may channel the matter of interstellar filaments into growing dense cores. The NIKA2 polarimeter consists of a room temperature continuously rotating multi-mesh HWP and a cold polarizer that separates the two orthogonal polarizations onto two 260 GHz KIDs arrays. We describe in this paper the preliminary results obtained during the most recent commissioning campaign performed in December 2018. We concentrate here on the analysis of the extended sources, while the observation of compact sources is presented in a companion paper [12]. We present preliminary NIKA2 polarization maps of the Crab nebula. We find that the integrated polarization intensity flux measured by NIKA2 is consistent with expectations.In terms of polarization angle, we are still limited by systematic uncertainties that will be further investigated in the forthcoming commissioning campaigns.

preprint2019arXiv

The NIKA polarimeter on science targets. Crab nebula observations at 150 GHz and dual-band polarization images of Orion Molecular Cloud OMC-1

We present here the polarization system of the NIKA camera and give a summary of the main results obtained and performed studies on Orion and the Crab nebula. The polarization system was equipped with a room temperature continuously rotating multi-mesh half wave plate and a grid polarizer facing the NIKA cryostat window. NIKA even though less sensitive than NIKA2 had polarization capability in both 1 and 2 millimiter bands. NIKA polarization observations demonstrated the ability of such a technology in detecting the polarization of different targets, compact and extended sources like the Crab nebula and Orion Molecular Cloud region OMC-1. These measurements together with the developed techniques to deal with systematics, opened the way to the current observations of NIKA2 in polarization that will provide important advances in the studies of galactic and extra-galactic emission and magnetic fields.

preprint2019arXiv

Unveiling the merger dynamics of the most massive MaDCoWS cluster at $z = 1.2$ from a multi-wavelength mapping of its intracluster medium properties

The characterization of the Intra-Cluster Medium (ICM) properties of high-redshift galaxy clusters is fundamental to our understanding of large-scale structure formation processes. We present the results of a multi-wavelength analysis of the very massive cluster MOO J1142$+$1527 at a redshift $z = 1.2$ discovered as part of the Massive and Distant Clusters of WISE Survey (MaDCoWS). This analysis is based on high angular resolution $Chandra$ X-ray and NIKA2 Sunyaev-Zel'dovich (SZ) data. Although the X-ray data have only about 1700 counts, we are able to determine the ICM thermodynamic radial profiles, namely temperature, entropy, and hydrostatic mass. These have been obtained with unprecedented precision at this redshift and up to $0.7R_{500}$, thanks to the combination of high-resolution X-ray and SZ data. The comparison between the galaxy distribution mapped in infrared by $Spitzer$ and the morphological properties of the ICM derived from the combined analysis of the $Chandra$ and NIKA2 data leads us to the conclusion that the cluster is an on-going merger. We measure the hydrostatic mass profile of the cluster in four angular sectors centered on the large-scale X-ray centroid. This allows us to estimate a systematic uncertainty on the cluster total mass that characterizes both the impact of the observed deviations from spherical symmetry and of the core dynamics on the mass profile. We further combine the X-ray and SZ data at the pixel level to obtain maps of the temperature and entropy distributions averaged along the line of sight. We find a relatively low entropy core at the position of the X-ray peak and high temperature regions located on its south and west sides. The increase in ICM temperature at the location of the SZ peak is expected given the merger dynamics. (abridged)

preprint2018arXiv

NIKA 150 GHz polarization observations of the Crab nebula and its spectral energy distribution

The Crab nebula is a supernova remnant exhibiting a highly polarized synchrotron radiation at radio and millimeter wavelengths. It is the brightest source in the microwave sky with an extension of 7 by 5 arcminutes and commonly used as a standard candle for any experiment which aims at measuring the polarization of the sky. Though its spectral energy distribution has been well characterized in total intensity, polarization data are still lacking at millimetre wavelengths. We report in this paper high resolution (18 arcsec FWHM) observations of the Crab nebula in total intensity and linear polarization at 150 GHz with the NIKA camera. NIKA, operated at the IRAM 30 m telescope from 2012 to 2015, is a camera made of Lumped Element Kinetic Inductance Detectors (LEKIDs) observing the sky at 150 and 260 GHz. From these observations we are able to reconstruct the spatial distribution of the polarization degree and angle of the Crab nebula, which is found to be compatible with previous observations at lower and higher frequencies. Averaging across the source and using other existing data sets we find that the Crab nebula polarization angle is consistent with being constant over a wide range of frequencies with a value of -87.7$^\circ$ +- 0.3 in Galactic coordinates. We also present the first estimation of the Crab nebula spectral energy distribution polarized flux in a wide frequency range: 30-353 GHz. Assuming a single power law emission model we find that the polarization spectral index $β_{pol}$ = - 0.347 +- 0.026 is compatible with the intensity spectral index $β$ = - 0.323 +- 0.001.

preprint2017arXiv

Search for High-energy Neutrinos from Binary Neutron Star Merger GW170817 with ANTARES, IceCube, and the Pierre Auger Observatory

The Advanced LIGO and Advanced Virgo observatories recently discovered gravitational waves from a binary neutron star inspiral. A short gamma-ray burst (GRB) that followed the merger of this binary was also recorded by the Fermi Gamma-ray Burst Monitor (Fermi-GBM), and the Anticoincidence Shield for the Spectrometer for the International Gamma-Ray Astrophysics Laboratory (INTEGRAL), indicating particle acceleration by the source. The precise location of the event was determined by optical detections of emission following the merger. We searched for high-energy neutrinos from the merger in the GeV--EeV energy range using the ANTARES, IceCube, and Pierre Auger Observatories. No neutrinos directionally coincident with the source were detected within $\pm500$ s around the merger time. Additionally, no MeV neutrino burst signal was detected coincident with the merger. We further carried out an extended search in the direction of the source for high-energy neutrinos within the 14-day period following the merger, but found no evidence of emission. We used these results to probe dissipation mechanisms in relativistic outflows driven by the binary neutron star merger. The non-detection is consistent with model predictions of short GRBs observed at a large off-axis angle.

preprint2016arXiv

Dense gas in low-metallicity galaxies

Stars form out of the densest parts of molecular clouds. Far-IR emission can be used to estimate the Star Formation Rate (SFR) and high dipole moment molecules, typically HCN, trace the dense gas. A strong correlation exists between HCN and Far-IR emission, with the ratio being nearly constant, over a large range of physical scales. A few recent observations have found HCN to be weak with respect to the Far-IR and CO in subsolar metallicity (low-Z) objects. We present observations of the Local Group galaxies M33, IC10, and NGC6822 with the IRAM 30meter and NRO 45m telescopes, greatly improving the sample of low-Z galaxies observed. HCN, HCO$^+$, CS, C$_2$H, and HNC have been detected. Compared to solar metallicity galaxies, the Nitrogen-bearing species are weak (HCN, HNC) or not detected (CN, HNCO, N$_2$H$^+$) relative to Far-IR or CO emission. HCO$^+$ and C$_2$H emission is normal with respect to CO and Far-IR. While $^{13}$CO is the usual factor 10 weaker than $^{12}$CO, C$^{18}$O emission was not detected down to very low levels. Including earlier data, we find that the HCN/HCO$^+$ ratio varies with metallicity (O/H) and attribute this to the sharply decreasing Nitrogen abundance. The dense gas fraction, traced by the HCN/CO and HCO$^+$/CO ratios, follows the SFR but in the low-Z objects the HCO$^+$ is much easier to measure. Combined with larger and smaller scale measurements, the HCO$^+$ line appears to be an excellent tracer of dense gas and varies linearly with the SFR for both low and high metallicities.

preprint2016arXiv

Dust Properties in HII Regions in M33

The conversion of the IR emission into star formation rate can be strongly dependent on the physical properties of the dust, which are affected by the environmental conditions where the dust is embedded. We study here the dust properties of a set of HII regions in the Local Group Galaxy M33 presenting different spatial configurations between the stars, gas and dust to understand the dust evolution under different environments. We model the SED of each region using the DustEM tool and obtain the mass relative to hydrogen for Very Small Grains (YVSG), Polycyclic Aromatic Hydrocarbons (YPAH) and Big Grains (YBG). The relative mass of the VSGs (YVSG/YTOT) is a factor of 1.7 higher for HII regions classified as filled and mixed than for regions presenting a shell structure. The enhancement of VSGs within NGC 604 and NGC 595 is correlated to expansive gas structures with velocities greater than 50 km/s. The gas-to-dust ratio derived for the HII regions in our sample exhibits two regimes related to the HI-H2 transition of the ISM. Regions corresponding to the HI diffuse regime present a gas-to-dust ratio compatible with the expected value if we assume that the gas-to-dust ratio scales linearly with metallicity, while regions corresponding to a H2 molecular phase present a flatter dust-gas surface density distribution. The fraction of VSGs can be affected by the conditions of the interstellar environment: strong shocks of 50-90 km/s existing in the interior of the most luminous HII regions can lead to fragmentation of BGs into smaller ones, while the more evolved shell and clear shell objects provide a more quiescent environment where reformation of dust BG grains might occur. The gas-to-dust variations found in this analysis might imply that grain coagulation and/or gas-phase metals incorporation to the dust mass is occurring in the interior of the HII regions in M33.

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

High angular resolution Sunyaev-Zel'dovich observations of MACS J1423.8+2404 with NIKA: Multiwavelength analysis

The prototype of the NIKA2 camera, NIKA, is an instrument operating at the IRAM 30-m telescope, which can observe simultaneously at 150 and 260GHz. One of the main goals of NIKA2 is to measure the pressure distribution in galaxy clusters at high resolution using the thermal SZ (tSZ) effect. Such observations have already proved to be an excellent probe of cluster pressure distributions even at high redshifts. However, an important fraction of clusters host submm and/or radio point sources, which can significantly affect the reconstructed signal. Here we report on <20" resolution observations at 150 and 260GHz of the cluster MACSJ1424, which hosts both radio and submm point sources. We examine the morphology of the tSZ signal and compare it to other datasets. The NIKA data are combined with Herschel satellite data to study the SED of the submm point source contaminants. We then perform a joint reconstruction of the intracluster medium (ICM) electronic pressure and density by combining NIKA, Planck, XMM-Newton, and Chandra data, focusing on the impact of the radio and submm sources on the reconstructed pressure profile. We find that large-scale pressure distribution is unaffected by the point sources because of the resolved nature of the NIKA observations. The reconstructed pressure in the inner region is slightly higher when the contribution of point sources are removed. We show that it is not possible to set strong constraints on the central pressure distribution without accurately removing these contaminants. The comparison with X-ray only data shows good agreement for the pressure, temperature, and entropy profiles, which all indicate that MACSJ1424 is a dynamically relaxed cool core system. The present observations illustrate the possibility of measuring these quantities with a relatively small integration time, even at high redshift and without X-ray spectroscopy.

preprint2016arXiv

High angular resolution SZ observations with NIKA and NIKA2

NIKA2 (New IRAM KID Arrays) is a dual band (150 and 260 GHz) imaging camera based on Kinetic Inductance Detectors (KIDs) and designed to work at the IRAM 30 m telescope (Pico Veleta, Spain). Built on the experience of the NIKA prototype, NIKA2 has been installed at the 30 m focal plane in October 2015 and the commissioning phase is now ongoing. Through the thermal Sunyaev-Zeldovich (tSZ) effect, NIKA2 will image the ionized gas residing in clusters of galaxies with a resolution of 12 and 18 arcsec FWHM (at 150 and 260 GHz, respectively). We report on the recent tSZ measurements with the NIKA camera and discuss the future objectives for the NIKA2 SZ large Program, 300h of observation dedicated to SZ science. With this program we intend to perform a high angular resolution follow-up of a cosmologically-representative sample of clusters belonging to SZ catalogues, with redshift greater than 0.5. The main output of the program will be the study of the redshift evolution of the cluster pressure profile as well as that of the scaling laws relating the cluster global properties.

preprint2016arXiv

High-resolution tSZ cartography of clusters of galaxies with NIKA at the IRAM 30-m telescope

The thermal Sunyaev-Zeldovich effect (tSZ) is a powerful probe to study clusters of galaxies and is complementary with respect to X-ray, lensing or optical observations. Previous arcmin resolution tSZ observations ({\it e.g.} SPT, ACT and Planck) only enabled detailed studies of the intra-cluster medium morphology for low redshift clusters ($z < 0.2$). Thus, the development of precision cosmology with clusters requires high angular resolution observations to extend the understanding of galaxy cluster towards high redshift. NIKA2 is a wide-field (6.5 arcmin field of view) dual-band camera, operated at $100 \ {\rm mK}$ and containing $\sim 3300$ KID (Kinetic Inductance Detectors), designed to observe the millimeter sky at 150 and 260 GHz, with an angular resolution of 18 and 12 arcsec respectively. The NIKA2 camera has been installed on the IRAM 30-m telescope (Pico Veleta, Spain) in September 2015. The NIKA2 tSZ observation program will allow us to observe a large sample of clusters (50) at redshift ranging between 0.5 and 1. As a pathfinder for NIKA2, several clusters of galaxies have been observed at the IRAM 30-m telescope with the NIKA prototype to cover the various configurations and observation conditions expected for NIKA2.

preprint2016arXiv

Millimeter and Submillimeter Excess Emission in M33 revealed by Planck and LABOCA

Previous studies have shown the existence of an excess of emission at submillimeter (submm) and millimeter (mm) wavelengths in the spectral energy distribution (SED) of many low-metallicity galaxies. The goal of the present study is to model separately the emission from the star forming (SF) component and the emission from the diffuse interstellar medium (ISM) in the nearby spiral galaxy M33. We decomposed the observed SED of M33 into its SF and diffuse components. Mid-infrared (MIR) and far-infrared (FIR) fluxes were extracted from Spitzer and Herschel data. At submm and mm wavelengths, we used ground-based observations from APEX to measure the emission from the SF component and data from the Planck space telescope to estimate the diffuse emission. Both components were separately fitted using radiation transfer models based on standard dust properties and a realistic geometry. Both modeled SEDs were combined to build the global SED of M33. In addition, the radiation field necessary to power the dust emission in our modeling was compared with observations from GALEX, Sloan, and Spitzer. Our modeling is able to reproduce the observations at MIR and FIR wavelengths, but we found a strong excess of emission at submm and mm wavelengths, where the model expectations severely underestimate the LABOCA and Planck fluxes. We also found that the ultraviolet (UV) radiation escaping the galaxy is 70% higher than the model predictions. We determined a gas-to-dust mass ratio Gdust~100, significantly lower than the value expected from the sub-solar metallicity of M33. We discussed different hypotheses to explain the discrepancies found in our study (i.e., excess of emission at submm and mm wavelengths, deficit of UV attenuation, and abnormally low value of Gdust), concluding that different dust properties in M33 is the most plausible explanation.

preprint2016arXiv

Molecular and atomic gas along and across the main sequence of star-forming galaxies

We use spectra from the ALFALFA, GASS and COLD GASS surveys to quantify variations in the mean atomic and molecular gas mass fractions throughout the SFR-M* plane and along the main sequence (MS) of star-forming galaxies. Although galaxies well below the MS tend to be undetected in the Arecibo and IRAM observations, reliable mean atomic and molecular gas fractions can be obtained through a spectral stacking technique. We find that the position of galaxies in the SFR-M* plane can be explained mostly by their global cold gas reservoirs as observed in the HI line, with in addition systematic variations in the molecular-to-atomic ratio and star formation efficiency. When looking at galaxies within +/-0.4 dex of the MS, we find that as stellar mass increases, both atomic and molecular gas mass fractions decrease, stellar bulges become more prominent, and the mean stellar ages increase. Both star formation efficiency and molecular-to-atomic ratios vary little for massive main sequence galaxies, indicating that the flattening of the MS is due to the global decrease of the cold gas reservoirs of galaxies rather than to bottlenecks in the process of converting cold atomic gas to stars.

preprint2016arXiv

NIKA 2: next-generation continuum/polarized camera at the IRAM 30 m telescope and its prototype

NIKA 2 (New Instrument of Kids Array) is a next generation continuum and polarized instrument successfully installed in October 2015 at the IRAM 30 m telescope on Pico-Veleta (Granada, Spain). NIKA 2 is a high resolution dual-band camera, operating with frequency multiplexed LEKIDs (Lumped Element Kinetic Inductance Detectors) cooled at 100 mK. Dual color images are obtained thanks to the simultaneous readout of a 1020 pixels array at 2 mm and 1140 x 2 pixels arrays at 1.15 mm with a final resolution of 18 and 12 arcsec respectively, and 6.5 arcmin of Field of View (FoV). The two arrays at 1.15 mm allow us to measure the linear polarization of the incoming light. This will place NIKA 2 as an instrument of choice to study the role of magnetic fields in the star formation process. The NIKA experiment, a prototype for NIKA 2 with a reduced number of detectors (about 400 LEKIDs) and FoV (1.8 arcmin), has been successfully operated at the IRAM 30 telescope in several open observational campaigns. The performance of the NIKA 2 polarization setup has been successfully validated with the NIKA prototype.

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 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 first CO+ image: Probing the HI/H2 layer around the ultracompact HII region Mon R2

The CO+ reactive ion is thought to be a tracer of the boundary between a HII region and the hot molecular gas. In this study, we present the spatial distribution of the CO+ rotational emission toward the Mon R2 star-forming region. The CO+ emission presents a clumpy ring-like morphology, arising from a narrow dense layer around the HII region. We compare the CO+ distribution with other species present in photon-dominated regions (PDR), such as [CII] 158 mm, H2 S(3) rotational line at 9.3 mm, polycyclic aromatic hydrocarbons (PAHs) and HCO+. We find that the CO+ emission is spatially coincident with the PAHs and [CII] emission. This confirms that the CO+ emission arises from a narrow dense layer of the HI/H2 interface. We have determined the CO+ fractional abundance, relative to C+ toward three positions. The abundances range from 0.1 to 1.9x10^(-10) and are in good agreement with previous chemical model, which predicts that the production of CO+ in PDRs only occurs in dense regions with high UV fields. The CO+ linewidth is larger than those found in molecular gas tracers, and their central velocity are blue-shifted with respect to the molecular gas velocity. We interpret this as a hint that the CO+ is probing photo-evaporating clump surfaces.

preprint2016arXiv

The NIKA2 instrument, a dual-band kilopixel KID array for millimetric astronomy

NIKA2 (New IRAM KID Array 2) is a camera dedicated to millimeter wave astronomy based upon kilopixel arrays of Kinetic Inductance Detectors (KID). The pathfinder instrument, NIKA, has already shown state-of-the-art detector performance. NIKA2 builds upon this experience but goes one step further, increasing the total pixel count by a factor $\sim$10 while maintaining the same per pixel performance. For the next decade, this camera will be the resident photometric instrument of the Institut de Radio Astronomie Millimetrique (IRAM) 30m telescope in Sierra Nevada (Spain). In this paper we give an overview of the main components of NIKA2, and describe the achieved detector performance. The camera has been permanently installed at the IRAM 30m telescope in October 2015. It will be made accessible to the scientific community at the end of 2016, after a one-year commissioning period. When this happens, NIKA2 will become a fundamental tool for astronomers worldwide.

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

First polarised light with the NIKA camera

NIKA is a dual-band camera operating with 315 frequency multiplexed LEKIDs cooled at 100 mK. NIKA is designed to observe the sky in intensity and polarisation at 150 and 260 GHz from the IRAM 30-m telescope. It is a test-bench for the final NIKA2 camera. The incoming linear polarisation is modulated at four times the mechanical rotation frequency by a warm rotating multi-layer Half Wave Plate. Then, the signal is analysed by a wire grid and finally absorbed by the LEKIDs. The small time constant (< 1ms ) of the LEKID detectors combined with the modulation of the HWP enables the quasi-simultaneous measurement of the three Stokes parameters I, Q, U, representing linear polarisation. In this paper we present results of recent observational campaigns demonstrating the good performance of NIKA in detecting polarisation at mm wavelength.

preprint2015arXiv

Measuring star formation with resolved observations: the test case of M33

Context. Measuring star formation at a local scale is important to constrain star formation laws. Yet, it is not clear whether and how the measure of star formation is affected by the spatial scale at which a galaxy is observed. Aims. We want to understand the impact of the resolution on the determination of the spatially resolved star formation rate (SFR) and other directly associated physical parameters such as the attenuation. Methods. We have carried out a multi-scale, pixel-by-pixel study of the nearby galaxy M33. Assembling FUV, Halpha, 8, 24, 70, and 100 micron maps, we have systematically compared the emission in individual bands with various SFR estimators from a resolution of 33 pc to 2084 pc. Results. We have found that there are strong, scale-dependent, discrepancies up to a factor 3 between monochromatic SFR estimators and Halpha+24 micron. The scaling factors between individual IR bands and the SFR show a strong dependence on the spatial scale and on the intensity of star formation. Finally, strong variations of the differential reddening between the nebular emission and the stellar continuum are seen, depending on the specific SFR (sSFR) and on the resolution. At the finest spatial scales, there is little differential reddening at high sSFR. The differential reddening increases with decreasing sSFR. At the coarsest spatial scales the differential reddening is compatible with the canonical value found for starburst galaxies. Conclusions. Our results confirm that monochromatic estimators of the SFR are unreliable at scales smaller than 1 kpc. Furthermore, the extension of local calibrations to high redshift galaxies presents non-trivial challenges as the properties of these systems may be poorly known.

preprint2015arXiv

Pressure distribution of the high-redshift cluster of galaxies CL J1226.9+3332 with NIKA

The thermal Sunyaev-Zel'dovich (tSZ) effect is expected to provide a low scatter mass proxy for galaxy clusters since it is directly proportional to the cluster thermal energy. The tSZ observations have proven to be a powerful tool for detecting and studying them, but high angular resolution observations are now needed to push their investigation to a higher redshift. In this paper, we report high angular (< 20 arcsec) resolution tSZ observations of the high-redshift cluster CL J1226.9+3332 (z=0.89). It was imaged at 150 and 260 GHz using the NIKA camera at the IRAM 30-meter telescope. The 150 GHz map shows that CL J1226.9+3332 is morphologically relaxed on large scales with evidence of a disturbed core, while the 260 GHz channel is used mostly to identify point source contamination. NIKA data are combined with those of Planck and X-ray from Chandra to infer the cluster's radial pressure, density, temperature, and entropy distributions. The total mass profile of the cluster is derived, and we find $M_{500} = 5.96^{+1.02}_{-0.79} $ x $10^{14} M_{\odot}$ within the radius $R_{500} = 930^{+50}_{-43}$ kpc, at a 68% confidence level. ($R_{500}$ is the radius within which the average density is 500 times the critical density at the cluster's redshift.) NIKA is the prototype camera of NIKA2, a KIDs (kinetic inductance detectors) based instrument to be installed at the end of 2015. This work is, therefore, part of a pilot study aiming at optimizing tSZ NIKA2 large programs.

preprint2015arXiv

Simultaneous multifrequency radio observations of the Galactic Centre magnetar SGR J1745-2900

We report on simultaneous observations of the magnetar SGR J1745-2900 at frequencies $ν= 2.54$ to $225\,\rm{GHz}$ using the Nancay 94-m equivalent, Effelsberg 100-m, and IRAM 30-m radio telescopes. We detect SGR J1745-2900 up to 225 GHz, the highest radio frequency detection of pulsed emission from a neutron star to date. Strong single pulses are also observed from 4.85 up to 154 GHz. At the millimetre band we see significant flux density and spectral index variabilities on time scales of tens of minutes, plus variability between days at all frequencies. Additionally, SGR J1745-2900 was observed at a different epoch at frequencies 296 to 472 GHz using the APEX 12-m radio telescope, with no detections. Over the period MJD 56859.83-56862.93 the fitted spectrum yields a spectral index of $\left<α\right> = -0.4 \pm 0.1$ for a reference flux density $\left< S_{154} \right> = 1.1 \pm 0.2\rm{\,mJy}$ (with $S_ν \propto ν^α)$, a flat spectrum alike those of the other radio-loud magnetars. These results show that strongly magnetized neutron stars can be effective radio emitters at frequencies notably higher to what was previously known and that pulsar searches in the Galactic Centre are possible in the millimetre band.

preprint2015arXiv

Temperatures of dust and gas in S~140

In dense parts of interstellar clouds (> 10^5 cm^-3), dust & gas are expected to be in thermal equilibrium, being coupled via collisions. However, previous studies have shown that the temperatures of the dust & gas may remain decoupled even at higher densities. We study in detail the temperatures of dust & gas in the photon-dominated region S 140, especially around the deeply embedded infrared sources IRS 1-3 and at the ionization front. We derive the dust temperature and column density by combining Herschel PACS continuum observations with SOFIA observations at 37 $μ$m and SCUBA at 450 $μ$m. We model these observations using greybody fits and the DUSTY radiative transfer code. For the gas part we use RADEX to model the CO 1-0, CO 2-1, 13CO 1-0 and C18O 1-0 emission lines mapped with the IRAM-30m over a 4' field. Around IRS 1-3, we use HIFI observations of single-points and cuts in CO 9-8, 13CO 10-9 and C18O 9-8 to constrain the amount of warm gas, using the best fitting dust model derived with DUSTY as input to the non-local radiative transfer model RATRAN. We find that the gas temperature around the infrared sources varies between 35 and 55K and that the gas is systematically warmer than the dust by ~5-15K despite the high gas density. In addition we observe an increase of the gas temperature from 30-35K in the surrounding up to 40-45K towards the ionization front, most likely due to the UV radiation from the external star. Furthermore, detailed models of the temperature structure close to IRS 1 show that the gas is warmer and/or denser than what we model. Finally, modelling of the dust emission from the sub-mm peak SMM 1 constrains its luminosity to a few ~10^2 Lo. We conclude that the gas heating in the S 140 region is very efficient even at high densities, most likely due to the deep UV penetration from the embedded sources in a clumpy medium and/or oblique shocks.

preprint2015arXiv

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

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

preprint2015arXiv

Velocity resolved [CII] spectroscopy of the center and the BCLMP302 region of M33 (HerM33es)

We aim to understand the contribution of the ionized, atomic and molecular phases of the ISM to the [CII] emission from clouds near the dynamical center and the BCLMP302 HII region in the north of the nearby galaxy M33 at a spatial resolution of 50pc. We combine high resolution [CII] spectra taken with the HIFI spectrometer onboard the Herschel satellite with [CII] Herschel-PACS maps and ground-based observations of CO(2-1) and HI. All data are at a common spatial resolution of 50pc. Typically, the [CII] lines have widths intermediate between the narrower CO(2-1) and broader HI line profiles. We decomposed the [CII] spectra in terms of contribution from molecular and atomic gas detected in CO(2-1) and HI, respectively. We find that the relative contribution of molecular and atomic gas traced by CO(2-1) and HI varies depends mostly on the local physical conditions and geometry. We estimate that 11-60% and 5-34% of the [CII] intensities in the center and in BCLMP302, respectively, arise at velocities showing no CO(2-1) or HI emission and could arise in CO-dark molecular gas. The deduced strong variation in the [CII] emission not associated with CO and HI cannot be explained in terms of differences in A_v, far-ultraviolet radiation field, and metallicity between the two studied regions. Hence the relative amounts of diffuse (CO-dark) and dense molecular gas possibly vary on spatial scales smaller than 50pc. Based on the emission measure observed at radio wavelengths we estimate the contribution of ionized gas at a few positions to lie between 10-25%. The correlations between the intensities of tracers corresponding to the same velocity range as [CII], differ from the correlation derived from PACS data. The results in this paper emphasize the need for velocity-resolved observations to discern the contribution of different components of the ISM to [CII] emission. (abridged)

preprint2014arXiv

Chemistry in isolation: High CCH/HCO+ line ratio in the AMIGA galaxy CIG 638

Multi-molecule observations towards an increasing variety of galaxies have been showing that the relative molecular abundances are affected by the type of activity. However, these studies are biased towards bright active galaxies, which are typically in interaction. We study the molecular composition of one of the most isolated galaxies in the local Universe where the physical and chemical properties of their molecular clouds have been determined by intrinsic mechanisms. We present 3 mm broad band observations of the galaxy CIG 638, extracted from the AMIGA sample of isolated galaxies. The emission of the J=1-0 transitions of CCH, HCN, HCO+, and HNC are detected. Integrated intensity ratios between these line are compared with similar observations from the literature towards active galaxies including starburst galaxies (SB), active galactic nuclei (AGN), luminous infrared galaxies (LIRG), and GMCs in M33. A significantly high ratio of CCH with respect to HCN, HCO+, and HNC is found towards CIG 638 when compared with all other galaxies where these species have been detected. This points to either an overabundance of CCH or to a relative lack of dense molecular gas as supported by the low HCN/CO ratio, or both. The data suggest that the CIG 638 is naturally a less perturbed galaxy where a lower fraction of dense molecular gas, as well as a more even distribution could explain the measured ratios. In this scenario the dense gas tracers would be naturally dimmer, while the UV enhanced CCH, would be overproduced in a less shielded medium.

preprint2014arXiv

Deuteration around the ultracompact HII region Mon R2

The massive star-forming region Mon R2 hosts the closest ultra-compact HII region that can be spatially resolved with current single-dish telescopes. We used the IRAM-30m telescope to carry out an unbiased spectral survey toward two important positions (namely IF and MP2), in order to studying the chemistry of deuterated molecules toward Mon R2. We found a rich chemistry of deuterated species at both positions, with detections of C2D, DCN, DNC, DCO+, D2CO, HDCO, NH2D, and N2D+ and their corresponding hydrogenated species and isotopologs. Our high spectral resolution observations allowed us to resolve three velocity components: the component at 10 km/s is detected at both positions and seems associated with the layer most exposed to the UV radiation from IRS 1; the component at 12 km/s is found toward the IF position and seems related to the molecular gas; finally, a component at 8.5 km/s is only detected toward the MP2 position, most likely related to a low-UV irradiated PDR. We derived the column density of all the species, and determined the deuterium fractions (Dfrac). The values of Dfrac are around 0.01 for all the observed species, except for HCO+ and N2H+ which have values 10 times lower. The values found in Mon R2 are well explained with pseudo-time-dependent gas-phase model in which deuteration occurs mainly via ion-molecule reactions with H2D+, CH2D+ and C2HD+. Finally, the [H13CN]/[HN13C] ratio is very high (~11) for the 10 km/s component, which also agree with our model predictions for an age of ~0.01-0.1 Myr. The deuterium chemistry is a good tool for studying star-forming regions. The low-mass star-forming regions seem well characterized with Dfrac(N2H+) or Dfrac(HCO+), but it is required a complete chemical modeling to date massive star-forming regions, because the higher gas temperature together with the rapid evolution of massive protostars.

preprint2014arXiv

First observation of the thermal Sunyaev-Zel'dovich effect with kinetic inductance detectors

Clusters of galaxies provide valuable information on the evolution of the Universe and large scale structures. Recent cluster observations via the thermal Sunyaev-Zel'dovich (tSZ) effect have proven to be a powerful tool to detect and study them. In this context, high resolution tSZ observations (~ tens of arcsec) are of particular interest to probe intermediate and high redshift clusters. Observations of the tSZ effect will be carried out with the millimeter dual-band NIKA2 camera, based on Kinetic Inductance Detectors (KIDs) to be installed at the IRAM 30-meter telescope in 2015. To demonstrate the potential of such an instrument, we present tSZ observations with the NIKA camera prototype, consisting of two arrays of 132 and 224 detectors that observe at 140 and 240 GHz with a 18.5 and 12.5 arcsec angular resolution, respectively. The cluster RX J1347.5-1145 was observed simultaneously at 140 and 240 GHz. We used a spectral decorrelation technique to remove the atmospheric noise and obtain a map of the cluster at 140 GHz. The efficiency of this procedure has been characterized through realistic simulations of the observations. The observed 140 GHz map presents a decrement at the cluster position consistent with the tSZ nature of the signal. We used this map to study the pressure distribution of the cluster by fitting a gNFW model to the data. Subtracting this model from the map, we confirm that RX J1347.5-1145 is an ongoing merger, which confirms and complements previous tSZ and X-ray observations. For the first time, we demonstrate the tSZ capability of KID based instruments. The NIKA2 camera with ~ 5000 detectors and a 6.5 arcmin field of view will be well-suited for in-depth studies of the intra cluster medium in intermediate to high redshifts, which enables the characterization of recently detected clusters by the Planck satellite.

preprint2014arXiv

High resolution SZ observations at the IRAM 30-m telescope with NIKA

High resolution observations of the thermal Sunyaev-Zel'dovich (tSZ) effect are necessary to allow the use of clusters of galaxies as a probe for large scale structures at high redshifts. With its high resolution and dual-band capability at millimeter wavelengths, the NIKA camera can play a significant role in this context. NIKA is based on newly developed Kinetic Inductance Detectors (KIDs) and operates at the IRAM 30m telescope, Pico Veleta, Spain. In this paper, we give the status of the NIKA camera, focussing on the KID technology. We then present observations of three galaxy clusters: RX J1347.5-1145 as a demonstrator of the NIKA capabilities and the recent observations of CL J1226.9+3332 (z = 0.89) and MACS J0717.5+3745 (z = 0.55). We also discuss prospects for the final NIKA2 camera, which will have a 6.5 arcminute field of view with about 5000 detectors in two bands at 150 and 260 GHz.

preprint2014arXiv

Latest NIKA results and the NIKA-2 project

NIKA (New IRAM KID Arrays) is a dual-band imaging instrument installed at the IRAM (Institut de RadioAstronomie Millimetrique) 30-meter telescope at Pico Veleta (Spain). Two distinct Kinetic Inductance Detectors (KID) focal planes allow the camera to simultaneously image a field-of-view of about 2 arc-min in the bands 125 to 175 GHz (150 GHz) and 200 to 280 GHz (240 GHz). The sensitivity and stability achieved during the last commissioning Run in June 2013 allows opening the instrument to general observers. We report here the latest results, in particular in terms of sensitivity, now comparable to the state-of-the-art Transition Edge Sensors (TES) bolometers, relative and absolute photometry. We describe briefly the next generation NIKA-2 instrument, selected by IRAM to occupy, from 2015, the continuum imager/polarimeter slot at the 30-m telescope.

preprint2014arXiv

Performance and calibration of the NIKA camera at the IRAM 30 m telescope

The New IRAM KID Array (NIKA) instrument is a dual-band imaging camera operating with Kinetic Inductance Detectors (KID) cooled at 100 mK. NIKA is designed to observe the sky at wavelengths of 1.25 and 2.14 mm from the IRAM 30 m telescope at Pico Veleta with an estimated resolution of 13\,arcsec and 18 arcsec, respectively. This work presents the performance of the NIKA camera prior to its opening to the astrophysical community as an IRAM common-user facility in early 2014. NIKA is a test bench for the final NIKA2 instrument to be installed at the end of 2015. The last NIKA observation campaigns on November 2012 and June 2013 have been used to evaluate this performance and to improve the control of systematic effects. We discuss here the dynamical tuning of the readout electronics to optimize the KID working point with respect to background changes and the new technique of atmospheric absorption correction. These modifications significantly improve the overall linearity, sensitivity, and absolute calibration performance of NIKA. This is proved on observations of point-like sources for which we obtain a best sensitivity (averaged over all valid detectors) of 40 and 14 mJy.s$^{1/2}$ for optimal weather conditions for the 1.25 and 2.14 mm arrays, respectively. NIKA observations of well known extended sources (DR21 complex and the Horsehead nebula) are presented. This performance makes the NIKA camera a competitive astrophysical instrument.

preprint2014arXiv

The Heating of Mid-Infrared Dust in the Nearby Galaxy M33: A Testbed for Tracing Galaxy Evolution

Because the 8 μm polycyclic aromatic hydrocarbon (PAH) emission has been found to correlate with other well-known star formation tracers, it has widely been used as a star formation rate (SFR) tracer. There are, however, studies that challenge the accuracy and reliability of the 8 μm emission as a SFR tracer. Our study, part of the Herschel M33 Extended Survey (HERM33ES) open time key program, aims at addressing this issue by analyzing the infrared emission from the nearby spiral galaxy M33 at the high spatial scale of 75 pc. Combining data from the Herschel Space Observatory and the Spitzer Space Telescope we find that the 8 μm emission is better correlated with the 250 μm emission, which traces cold interstellar gas, than with the 24 μm emission. The L(8)/L(24) ratio is highly depressed in 24 μm luminous regions, which correlate with known HII regions. We also compare our results with the dust emission models by Draine & Li (2007). We confirm that the depression of 8 μm PAH emission near star-forming regions is higher than what is predicted by models; this is possibly an effect of increased stellar radiation from young stars destroying the dust grains responsible for the 8 μm emission as already suggested by other authors. We find that the majority of the 8 μm emission is fully consistent with heating by the diffuse interstellar medium, similar to what recently determined for the dust emission in M31 by Draine at al. (2013). We also find that the fraction of 8 μm emission associated with the diffuse interstellar radiation field ranges between 60% and 80% and is 40% larger than the diffuse fraction at 24 μm.

preprint2014arXiv

The IRAM M33 CO(2-1) Survey - A complete census of the molecular gas out to 7 kpc

In order to study the ISM and the interplay between the atomic and molecular components in a low-metallicity environment, we present a complete high angular and spectral resolution map and data cube of the 12CO(2-1) emission from the Local Group galaxy M33. Its metallicity is roughly half-solar, such that we can compare its ISM with that of the Milky Way with the main changes being the metallicity and the gas mass fraction. The data have a 12" angular resolution (50pc) with a spectral resolution of 2.6 km/s and a mean noise level of 20 mK per channel in antenna temperature. A radial cut along the major axis was also observed in the 12CO(1-0) line. The CO data cube and integrated intensity map are optimal when using HI data to define the baseline window and the velocities over which the CO emission is integrated. Great care was taken when building these maps, testing different windowing and baseline options and investigating the effect of error beam pickup. The total CO(2-1) luminosity is 2.8e7 K km/s pc2, following the spiral arms in the inner disk. There is no clear variation in the CO(2-1/1-0) intensity ratio with radius and the average value is roughly 0.8. The total molecular gas mass is estimated, using a N(H2)/Ico(1-0)=4e20 cm-2/(K km/s) conversion factor, to be 3.1e8 Msol. The CO spectra in the cube were shifted to zero velocity by subtracting the velocity of the HI peak from the CO spectra. Hence, the velocity dispersion between the atomic and molecular components is extremely low, independently justifying the use of the HI line in building our maps. Stacking the spectra in concentric rings shows that the CO linewidth and possibly the CO-HI velocity dispersion decrease in the outer disk. Using the CO(2-1) emission to trace the molecular gas, the probability distribution function of the H2 column density shows an excess at high column density above a log normal distribution.

preprint2014arXiv

The molecular circumnuclear disk (CND) in Centaurus A: A multi-transition CO and [CI] survey with Herschel, APEX, JCMT, and SEST

We present new CO and C^o line measurements of the compact circumnuclear disk in the center of NGC 128 (Centaurus~A) obtained with the Herschel Space Observatory, as well as SEST, JCMT, and APEX. The Cen A center CO ladder is quite different from those of either star-burst galaxies or AGNs. In addition, the relative intensity of the central Cen A [CI] emission lines is much greate than that in any other galaxy. The CO surface brightness of the compact circumnuclear disk (CND) is significantly higher than that of the much more extended thin disk (ETD) in the same line of sight. Our LVG and PDR/XDR models suggest that much of the CND gas is relatively cool (25 - 80 K) and not very dense (~ 300 cm^{-3}) if the heating is by UV photons, although there is some gas in both the CND and the ETD with a much higher density of ~30 000 cm^{-3}. Finally, there is also high-excitation, high-density phase in the CND (but not in the ETD), either in the form of an extreme PDR but more likely in the form of an XDR. The total gas mass of the Cen A CND is 8.4 x 10^{7} M(sun), uncertain by a factor of two. The CO-H2 conversion factor is 4 x 10^{20} K km/s, also within a factor of two.

preprint2014arXiv

The NIKA 2012-2014 observing campaigns: control of systematic effects and results

The New IRAM KID Array (NIKA) is a dual-band camera operating with frequency multiplexed arrays of Lumped Element Kinetic Inductance Detectors (LEKIDs) cooled to 100 mK. NIKA is designed to observe the intensity and polarisation of the sky at 1.25 and 2.14 mm from the IRAM 30 m telescope. We present the improvements on the control of systematic effects and astrophysical results made during the last observation campaigns between 2012 and 2014.

preprint2013arXiv

Dense Gas in M33 (HerM33es)

We aim to better understand the emission of molecular tracers of the diffuse and dense gas in giant molecular clouds and the influence that metallicity, optical extinction, density, far-UV field, and star formation rate have on these tracers. Using the IRAM 30m telescope, we detected HCN, HCO+, 12CO, and 13CO in six GMCs along the major axis of M33 at a resolution of ~ 114pc and out to a radial distance of 3.4kpc. Optical, far-infrared, and submillimeter data from Herschel and other observatories complement these observations. To interpret the observed molecular line emission, we created two grids of models of photon-dominated regions, one for solar and one for M33-type subsolar metallicity. The observed HCO+/HCN line ratios range between 1.1 and 2.5. Similarly high ratios have been observed in the Large Magellanic Cloud. The HCN/CO ratio varies between 0.4% and 2.9% in the disk of M33. The 12CO/13CO line ratio varies between 9 and 15 similar to variations found in the diffuse gas and the centers of GMCs of the Milky Way. Stacking of all spectra allowed HNC and C2H to be detected. The resulting HCO+/HNC and HCN/HNC ratios of ~ 8 and 6, respectively, lie at the high end of ratios observed in a large set of (ultra-)luminous infrared galaxies. HCN abundances are lower in the subsolar metallicity PDR models, while HCO+ abundances are enhanced. For HCN this effect is more pronounced at low optical extinctions. The observed HCO+/HCN and HCN/CO line ratios are naturally explained by subsolar PDR models of low optical extinctions between 4 and 10 mag and of moderate densities of n = 3x10^3 - 3x10^4 cm^-3, while the FUV field strength only has a small effect on the modeled line ratios. The line ratios are almost equally well reproduced by the solar-metallicity models, indicating that variations in metallicity only play a minor role in influencing these line ratios.

preprint2013arXiv

Detection of the tSZ effect with the NIKA camera

We present the first detection of the thermal Sunyaev-Zel'dovich (tSZ) effect from a cluster of galaxies performed with a KIDs (Kinetic Inductance Detectors) based instrument. The tSZ effect is a distortion of the black body CMB (Cosmic Microwave Background) spectrum produced by the inverse Compton interaction of CMB photons with the hot electrons of the ionized intra-cluster medium. The massive, intermediate redshift cluster RX J1347.5-1145 has been observed using NIKA (New IRAM KIDs arrays), a dual-band (140 and 240 GHz) mm-wave imaging camera, which exploits two arrays of hundreds of KIDs: the resonant frequencies of the superconducting resonators are shifted by mm-wave photons absorption. This tSZ cluster observation demonstrates the potential of the next generation NIKA2 instrument, being developed for the 30m telescope of IRAM, at Pico Veleta (Spain). NIKA2 will have 1000 detectors at 140GHz and 2x2000 detectors at 240GHz, providing in that band also a measurement of the linear polarization. NIKA2 will be commissioned in 2015.

preprint2013arXiv

Gas and dust cooling along the major axis of M33 (HerM33es): ISO/LWS CII observations

We aim to better understand the heating of the gas by observing the prominent gas cooling line [CII] at 158um in the low-metallicity environment of the Local Group spiral galaxy M33 at scales of 280pc. In particular, we aim at describing the variation of the photoelectric heating efficiency with galactic environment. In this unbiased study, we used ISO/LWS [CII] observations along the major axis of M33, in combination with Herschel PACS and SPIRE continuum maps, IRAM 30m CO 2-1 and VLA HI data to study the variation of velocity integrated intensities. The ratio of [CII] emission over the far-infrared continuum is used as a proxy for the heating efficiency, and models of photon-dominated regions are used to study the local physical densities, FUV radiation fields, and average column densities of the molecular clouds. The heating efficiency stays constant at 0.8% in the inner 4.5kpc radius of the galaxy where it starts to increase to reach values of ~3% in the outskirts at about 6kpc radial distance. The rise of efficiency is explained in the framework of PDR models by lowered volume densities and FUV fields, for optical extinctions of only a few magnitudes at constant metallicity. In view of the significant fraction of HI emission stemming from PDRs, and for typical pressures found in the Galactic cold neutral medium (CNM) traced by HI emission, the CNM contributes ~15% to the observed [CII] emission in the inner 2kpc radius of M33. The CNM contribution remains largely undetermined in the south, while positions between 2 and 7.3kpc radial distance in the north of M33 show a contribution of ~40%+-20%.

preprint2013arXiv

Kinematics of the ionized-to-neutral interfaces in Monoceros R2

Context. Monoceros R2 (Mon R2), at a distance of 830 pc, is the only ultra-compact H ii region (UC H ii) where its associated photon-dominated region (PDR) can be resolved with the Herschel Space Observatory. Aims. Our aim is to investigate observationally the kinematical patterns in the interface regions (i.e., the transition from atomic to molecular gas) associated with Mon R2. Methods. We used the HIFI instrument onboard Herschel to observe the line profiles of the reactive ions CH+, OH+ and H2O+ toward different positions in Mon R2. We derive the column density of these molecules and compare them with gas-phase chemistry models. Results. The reactive ion CH+ is detected both in emission (at central and red-shifted velocities) and in absorption (at blue-shifted velocities). OH+ is detected in absorption at both blue- and red-shifted velocities, with similar column densities. H2O+ is not detected at any of the positions, down to a rms of 40 mK toward the molecular peak. At this position, we find that the OH+ absorption originates in a mainly atomic medium, and therefore is associated with the most exposed layers of the PDR. These results are consistent with the predictions from photo-chemical models. The line profiles are consistent with the atomic gas being entrained in the ionized gas flow along the walls of the cavity of the H ii region. Based on this evidence, we are able to propose a new geometrical model for this region. Conclusions. The kinematical patterns of the OH+ and CH+ absorption indicate the existence of a layer of mainly atomic gas for which we have derived, for the first time, some physical parameters and its dynamics.

preprint2013arXiv

Large scale IRAM 30m CO-observations in the giant molecular cloud complex W43

We aim to give a full description of the distribution and location of dense molecular clouds in the giant molecular cloud complex W43. It has previously been identified as one of the most massive star-forming regions in our Galaxy. To trace the moderately dense molecular clouds in the W43 region, we initiated an IRAM 30m large program, named W43-HERO, covering a large dynamic range of scales (from 0.3 to 140 pc). We obtained on-the-fly-maps in 13CO (2-1) and C18O (2-1) with a high spectral resolution of 0.1 km/s and a spatial resolution of 12". These maps cover an area of ~1.5 square degrees and include the two main clouds of W43, as well as the lower density gas surrounding them. A comparison with Galactic models and previous distance calculations confirms the location of W43 near the tangential point of the Scutum arm at a distance from the Sun of approximately 6 kpc. The resulting intensity cubes of the observed region are separated into sub-cubes, centered on single clouds which are then analyzed in detail. The optical depth, excitation temperature, and H2 column density maps are derived out of the 13CO and C18O data. These results are then compared with those derived from Herschel dust maps. The mass of a typical cloud is several 10^4 solar masses while the total mass in the dense molecular gas (>100 cm^-3) in W43 is found to be about 1.9e6 solar masses. Probability distribution functions obtained from column density maps derived from molecular line data and Herschel imaging show a log-normal distribution for low column densities and a power-law tail for high densities. A flatter slope for the molecular line data PDF may imply that those selectively show the gravitationally collapsing gas.

preprint2013arXiv

Low-velocity shocks traced by extended SiO emission along the W43 ridges: witnessing the formation of young massive clusters

The formation of high-mass stars is tightly linked to that of their parental clouds. We here focus on the high-density parts of W43, a molecular cloud undergoing an efficient event of formation. The cloud structure is studied with a column density image derived from Herschel continuum maps obtained at 70, 160, 250, 350, and 500 micron. We identify two high-column density filamentary clouds, quoted as the W43-MM1 and W43-MM2 ridges, which both account for 1.5x10^4 Msun gas mass above 10^23 cm-2 and within areas of 5 and 14pc^2, respectively. We used the N_2H^+ 1--0 line to confirm that the W43-MM1 and W43-MM2 ridges are structures coherent in velocity and gravitationally bound, despite their large velocity dispersion and ~5 kms line widths. The most intriguing result of the W43 large program is the bright wide-spread SiO 2--1 emission: 1--11 K kms$ stretching an area of ~28 pc^2. Concentrated toward the W43-MM1 and W43-MM2 ridges and their immediate surroundings, it leads to a total luminosity of L_SiO 2-1 ~4 10^4 K kms kpc^2pc^2. We measured a steep relation between the luminosity and velocity extent of the SiO~2--1 lines and propose to use it to distinguish the low-velocity shocks observed here from the more classical high-velocity ones associated with outflows of high-mass young stellar objects. We used state-of-the-art shock models to demonstrate that low-velocity (<10 kms^-1) shocks with a small amount (10%) of Si atoms initially in gas phase or in grain mantles can explain the observed SiO column density in W43. The spatial and velocity overlaps between the ridges high-density gas (n_H2>10^4-10^5 cm^-3) and the shocked SiO gas suggests that ridges could be forming via colliding flows driven by gravity and accompanied by low-velocity shocks. This mechanism may be the initial conditions for the formation of young massive clusters in these ridges.

preprint2013arXiv

Spatial distribution of small hydrocarbons in the neighborhood of the Ultra Compact HII region Monoceros R2

We study the chemistry of small hydrocarbons in the photon-dominated regions (PDRs) associated with the ultra-compact HII region Mon R2. Our goal is to determine the variations of the abundance of small hydrocarbons in a high-UV irradiated PDR and investigate their chemistry. We present an observational study of CH, CCH and c-C$_3$H$_2$ in Mon R2 combining data obtained with the IRAM 30m telescope and Herschel. We determine the column densities of these species, and compare their spatial distributions with that of polycyclic aromatic hydrocarbon (PAH). We compare the observational results with different chemical models to explore the relative importance of gas-phase, grain-surface and time-dependent chemistry in these environments. The emission of the small hydrocarbons show different patterns. The CCH emission is extended while CH and c-C$_3$H$_2$ are concentrated towards the more illuminated layers of the PDR. The ratio of the column densities of c-C$_3$H$_2$ and CCH shows spatial variations up to a factor of a few, increasing from $N(c-C$_3$H$_2$)/N(CCH)\approx0.004$ in the envelope to a maximum of $\sim0.015-0.029$ towards the 8$μ$m emission peak. Comparing these results with other galactic PDRs, we find that the abundance of CCH is quite constant over a wide range of G$_0$, whereas the abundance of c-C$_3$H$_2$ is higher in low-UV PDRs. In Mon R2, the gas-phase steady-state chemistry can account relatively well for the abundances of CH and CCH in the most exposed layers of the PDR, but falls short by a factor of 10 to reproduce c-C$_3$H$_2$. In the molecular envelope, time-dependent effects and grain surface chemistry play a dominant role in determining the hydrocarbons abundances. Our study shows that CCH and c-C$_3$H$_2$ present a complex chemistry in which UV photons, grain-surface chemistry and time dependent effects contribute to determine their abundances.

preprint2013arXiv

Spectral Energy Distributions of HII regions in M33 (HerM33es)

Within the framework of the Herschel M 33 extended survey HerM33es we study the Spectral Energy Distribution (SED) of a set of HII regions in M 33 as a function of the morphology. We present a catalogue of 119 HII regions morphologically classified: 9 filled, 47 mixed, 36 shell, and 27 clear shell HII regions. For each object we extract the photometry at twelve available wavelength bands (from FUV-1516A to IR-250mi) and obtain the SED. We also obtain emission line profiles across the regions to study the location of the stellar, ionised gas, and dust components. We find trends for the SEDs related to the morphology, showing that the star and gas-dust configuration affects the ratios of the emission in different bands. The mixed and filled regions show higher emission at 24mi than the shells and clear shells, which could be due to the proximity of the dust to the stellar clusters in the case of mixed and filled regions. The FIR peak for shells and clear shells seems to be located towards longer wavelengths, indicating that the dust is colder for this type of objects.The logarithmic 100/70mi ratio for filled and mixed regions remains constant over one order of magnitude in Halpha and FUV surface brightness, while the shells and clear shells exhibit a wider range of values of almost two orders of magnitude. We derive dust masses and temperatures fitting the individual SEDs with dust models proposed in the literature. The derived dust mass range is between 10^2-10^4 Msun and the cold dust temperature spans T(cold)~12-27 K. The spherical geometrical model proposed for the Halpha clear shells is confirmed by the emission profile obtained from the observations and is used to infer the electron density within the envelope: the typical electron density is 0.7+-0.3 cm^-3, while filled regions can reach values two to five times higher.

preprint2013arXiv

The NIKA instrument: results and perspectives towards a permanent KID based camera for the Pico Veleta observatory

The New IRAM KIDs Array (NIKA) is a pathfinder instrument devoted to millimetric astronomy. In 2009 it was the first multiplexed KID camera on the sky; currently it is installed at the focal plane of the IRAM 30-meters telescope at Pico Veleta (Spain). We present preliminary data from the last observational run and the ongoing developments devoted to the next NIKA-2 kilopixels camera, to be commissioned in 2015. We also report on the latest laboratory measurements, and recent improvements in detector cosmetics and read-out electronics. Furthermore, we describe a new acquisition strategy allowing us to improve the photometric accuracy, and the related automatic tuning procedure.

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

Variation in the dust emissivity index across M33 with Herschel and Spitzer (HerM33es)

We study the wavelength dependence of the dust emission as a function of position and environment across the disk of M33 at a linear resolution of 160 pc using Spitzer and Herschel photometric data. Expressing the emissivity of the dust as a power law, the power-law exponent (beta) is estimated from two independent approaches designed to properly treat the degeneracy between beta and the dust temperature. Both beta and the dust temperature are higher in the inner disk than in the outer disk, contrary to reported beta-T anti-correlations found in other sources. In the cold + warm dust model, the warm component and the ionized gas (Halpha) have a very similar distribution across the galaxy, demonstrating that the model separates the components in an appropriate fashion. The flocculent spiral arms and the dust lanes are evident in the map of the cold component. Both cold and warm dust column densities are high in star forming regions and reach their maxima toward the giant star forming complexes NGC604 and NGC595. beta declines from close to 2 in the center to about 1.3 in the outer disk. beta is positively correlated with star formation and with molecular gas column, as traced by Halpha and CO emission. The lower dust emissivity index in the outer parts of M33 is likely related to the reduced metallicity (different grain composition) and possibly different size distribution. It is not due to the decrease in stellar radiation field or temperature in a simple way because the FIR-bright regions in the outer disk also have a low beta. Like most spirals, M33 has a (decreasing) radial gradient in star formation and molecular-to-atomic gas ratio such that the regions bright in Halpha or CO tend to trace the inner disk, making it difficult to distinguish between their effects on the dust.

preprint2012arXiv

Cool and warm dust emission from M33 (HerM33es)

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

preprint2012arXiv

Dust and gas power-spectrum in M33 (HERM33ES)

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

preprint2012arXiv

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

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

preprint2012arXiv

Herschel / HIFI observations of CO, H2O and NH3 in Mon R2

Context. Mon R2 is the only ultracompact HII region (UCHII) where the associated photon-dominated region (PDR) can be resolved with Herschel. Due to its brightness and proximity, it is the best source to investigate the chemistry and physics of highly UV-irradiated PDRs. Aims. Our goal is to estimate the abundance of H2O and NH3 in this region and investigate their origin. Methods. We present new observations obtained with HIFI and the IRAM-30m telescope. Using a large velocity gradient approach, we model the line intensities and derive an average abundance of H2O and NH3 across the region. Finally, we model the line profiles with a non-local radiative transfer model and compare these results with the abundance predicted by the Meudon PDR code. Results. The variations of the line profiles and intensities indicate complex geometrical and kinematical patterns. The H2O lines present a strong absorption at the ambient velocity and emission in high velocity wings towards the HII region. The spatial distribution of the o-H2^18O line shows that the its emission arises in the PDR surrounding the HII region. By modeling the o-H2^18O emission we derive a mean abundance of o-H2O of ~10^-8 relative to H2. The ortho-H2O abundance is however larger, ~1x10^-7, in the high velocity wings. Possible explanations for this larger abundance include an expanding hot PDR and/or an outflow. Ammonia seems to be present only in the envelope with an average abundance of ~2x10^-9 relative to H2. Conclusions. The Meudon PDR code can account for the measured water abundance in the high velocity gas as long as we assume that it originates from a <1 mag hot expanding layer of the PDR, i.e. that the outflow has only a minor contribution to this emission. To explain the abundances in the rest of the cloud the molecular freeze out and grain surface chemistry would need to be included.

preprint2012arXiv

Spectral line survey of the ultracompact HII region Mon R2

Ultracompact (UC) HII regions constitute one of the earliest phases in the formation of a massive star and are characterized by extreme physical conditions (Go>10^5 Habing field and n>10^6 cm^-3). The UC HII Mon R2 is the closest one and therefore an excellent target to study the chemistry in these complex regions. We carried out a 3mm and 1mm spectral survey using the IRAM 30-m telescope towards three positions that represent different physical environments in Mon R2: (i) the ionization front (IF) at (0",0"); two peaks in the molecular cloud (ii) MP1 at the offset (+15",-15") and (iii) MP2 at the farther offset (0",40"). In addition, we carried out extensive modeling to explain the chemical differences between the three observed regions. We detected more than thirty different species. We detected SO+ and C4H suggesting that UV radiation plays an important role in the molecular chemistry of this region. We detected the typical PDR molecules CN, HCN, HCO, C2H, and c-C3H2. While the IF and the MP1 have a chemistry similar to that found in high UV field and dense PDRs like the Orion Bar, the MP2 is more similar to lower UV/density PDRs like the Horsehead nebula. We also detected complex molecules that are not usually found in PDRs (CH3CN, H2CO, HC3N, CH3OH and CH3C2H). Sulfur compounds CS, HCS+, C2S, H2CS, SO and SO2 and the deuterated species DCN and C2D were also identified. [DCN]/[HCN]=0.03 and [C2D]/[C2H]=0.05, are among the highest in warm regions. Our results show that the high UV/dense PDRs present a different chemistry from that of the low UV case. Abundance ratios like [CO+]/[HCO+] or [HCO]/[HCO+] are good diagnostics to differentiate between them. In Mon R2 we have the two classes of PDRs, a high UV PDR towards the IF and the adjacent molecular bar and a low-UV PDR which extends towards the north-west following the border of the cloud.

preprint2012arXiv

Spectrally resolved CII emission in M~33 ({\tt HerM33es}): Physical conditions and kinematics around BCLMP 691

This work presents high spectral resolution observations of the \CII\ line at 158 \micron, one of the major cooling lines of the interstellar medium, taken with the HIFI heterodyne spectrometer on the Herschel satellite. In BCLMP 691, an \HII\ region far north (3.3 kpc) in the disk of M 33, the \CII\ and CO line profiles show similar velocities within $0.5 \kms$, while the \HI\ line velocities are systematically shifted towards lower rotation velocities by $\sim 5\kms$. Observed at the same $12"$ angular resolution, the \CII\ lines are broader than those of CO by about 50% but narrower than the \HI\ lines. The \CII\ line to far-infrared continuum ratio suggests a photoelectric heating efficiency of 1.1%. The data, together with published models indicate a UV field $G_0 \sim 100$ in units of the solar neighborhood value, a gas density $n_H \sim 1000 \cc$, and a gas temperature $T\sim 200$ K. Adopting these values, we estimate the C$^+$ column density to be $N_{C^+} \approx 1.3 \times 10^{17} \cmt$. The \CII\ emission comes predominantly from the warm neutral region between the \HII\ region and the cool molecular cloud behind it. From published abundances, the inferred C$^+$ column corresponds to a hydrogen column density of $N_H \sim 2 \times 10^{21} \cmt$. The CO observations suggest that $N_H = 2 N_{H_2} \sim 3.2 \times 10^{21} \cmt$ and 21cm measurements, also at $12"$ resolution, yield $N_\HI \approx 1.2 \times 10^{21} \cmt$ within the \CII\ velocity range. Thus, some H$_2$ not detected in CO must be present, in agreement with earlier findings based on the SPIRE 250 -- 500 $μ$m emission.

preprint2012arXiv

The impact of interactions, bars, bulges, and AGN on star formation efficiency in local massive galaxies

Using observations from the GASS and COLD GASS surveys and complementary data from SDSS and GALEX, we investigate the nature of variations in gas depletion time observed across the local massive galaxy population. The large and unbiased COLD GASS sample allows us to assess the relative importance of galaxy interactions, bar instabilities, morphologies and the presence of AGN in regulating star formation efficiency. Both the H2 mass fraction and depletion time vary as a function of the distance of a galaxy from the main sequence in the SFR-M* plane. The longest gas depletion times are found in below-main sequence bulge-dominated galaxies that are either gas-poor, or else on average less efficient than disk-dominated galaxy at converting into stars any cold gas they may have. We find no link between AGN and these long depletion times. The galaxies undergoing mergers or showing signs of morphological disruptions have the shortest molecular gas depletion times, while those hosting strong stellar bars have only marginally higher global star formation efficiencies as compared to matched control samples. Our interpretation is that depletion time variations are caused by changes in the ratio between the gas mass traced by the CO(1-0) observations, and the gas mass in high density star-forming cores, with interactions, mergers and bar instabilities able to locally increase pressure and raise the ratio of efficiently star-forming gas to CO-detected gas. Building a sample representative of the local massive galaxy population, we derive a global Kennicutt-Schmidt relation of slope 1.18+/-0.24, and observe structure within the scatter around this relation, with galaxies having low (high) stellar mass surface densities lying systematically above (below) the mean relation, suggesting that gas surface density is not the only parameter driving the global star formation ability of a galaxy.

preprint2012arXiv

The structure of hot gas in Cepheus B

By observing radiation-affected gas in the Cepheus B molecular cloud we probe whether the sequential star formation in this source is triggered by the radiation from newly formed stars. We used the dual band receiver GREAT onboard SOFIA to map [C II] and CO 13--12 and 11--10 in Cep B and compared the spatial distribution and the spectral profiles with complementary ground-based data of low-$J$ transitions of CO isotopes, atomic carbon, and the radio continuum. The interaction of the radiation from the neighboring OB association creates a large photon-dominated region (PDR) at the surface of the molecular cloud traced through the photoevaporation of C^+. Bright internal PDRs of hot gas are created around the embedded young stars, where we detect evidence of the compression of material and local velocity changes; however, on the global scale we find no indications that the dense molecular material is dynamically affected.

preprint2011arXiv

A Constant Molecular Gas Depletion Time in Nearby Disk Galaxies

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

preprint2011arXiv

Dust heating sources in galaxies: the case of M33 (HERM33ES)

Dust emission is one of the main windows to the physics of galaxies and to star formation as the radiation from young, hot stars is absorbed by the dust and reemitted at longer wavelengths. The recently launched Herschel satellite now provides a view of dust emission in the far-infrared at an unequaled resolution and quality up to 500 μm. In the context of the Herschel HERM33ES open time key project, we are studying the moderately inclined Scd local group galaxy M33 which is located only 840 kpc away. In this article, using Spitzer and Herschel data ranging from 3.6 μm to 500 μm, along with HI, Hα maps, and GALEX ultraviolet data we have studied the emission of the dust at the high spatial resolution of 150 pc. Combining Spitzer and Herschel bands, we have provided new, inclination corrected, resolved estimators of the total infrared brightness and of the star formation rate from any combination of these bands. The study of the colors of the warm and cold dust populations shows that the temperature of the former is, at high brightness, dictated by young massive stars but, at lower brightness, heating is taken over by the evolved populations. Conversely, the temperature of the cold dust is tightly driven by the evolved stellar populations.

preprint2011arXiv

Giant Molecular Clouds in the Local Group Galaxy M33

We present an analysis of the systematic CO(2-1) survey at 12" resolution covering most of the local group spiral M 33 which, at a distance of 840 kpc, is close enough that individual giant molecular clouds (GMCs) can be identified. The goal of this work is to study the properties of the GMCs in this subsolar metallicity galaxy. The CPROPS (Cloud Properties) algorithm (Rosolowsky & Leroy 2006) was used to identify 337 GMCs in M 33, the largest sample to date in an external galaxy. The sample is used to study the GMC luminosity function, or mass spectrum under the assumption of a constant N(H2)/ICO ratio. We find that n(L)dL = K*L^(-2.0\pm0.1) for the entire sample. However, when the sample is divided into inner and outer disk samples, the exponent changes from 1.6 \pm 0.2 for the centre 2 kpc to 2.3 \pm 0.2 for galactocentric distances larger than 2 kpc. Based on the emission in the FUV, Halpha, 8mu, and 24mu bands, each cloud was classified in terms of its star forming activity - no star formation, embedded, or exposed star formation (visible in FUV and Halpha). At least one sixth of the clouds had no (massive) star formation, suggesting that the average time required for star formation to start is about one sixth of the total time for which the object is identifiable as a GMC. The clouds without star formation have significantly lower CO luminosities than those with star formation, whether embedded or exposed, presumably related to the lack of heating sources. Taking the cloud sample as a whole, the main non-trivial correlation was the decrease in cloud CO brightness (or luminosity) with galactocentric radius. The complete cloud catalog, including CO and HI spectra and the CO contours on the FUV, Halpha, 8mu, and 24mu images is presented in the appendix.

preprint2011arXiv

Spectral Energy Distributions of a set of HII regions in M33 (HerM33es)

Within the framework of the HerM33es Key Project for Herschel and in combination with multi-wavelength data, we study the Spectral Energy Distribution (SED) of a set of HII regions in the Local Group Galaxy M33. Using the Halpha emission, we perform a classification of a selected HII region sample in terms of morphology, separating the objects in filled, mixed, shell and clear shell objects. We obtain the SED for each HII region as well as a representative SED for each class of objects. We also study the emission distribution of each band within the regions. We find different trends in the SEDs for each morphological type that are related to properties of the dust and their associated stellar cluster. The emission distribution of each band within the region is different for each morphological type of object.

preprint2011arXiv

Star formation in M33 (HerM33es)

Within the key project "Herschel M33 extended survey" (HerM33es), we are studying the physical and chemical processes driving star formation and galactic evolution in the nearby galaxy M33, combining the study of local conditions affecting individual star formation with properties only becoming apparent on global scales. Here, we present recent results obtained by the HerM33es team. Combining Spitzer and Herschel data ranging from 3.6um to 500um, along with HI, Halpha, and GALEX UV data, we have studied the dust at high spatial resolutions of 150pc, providing estimators of the total infrared (TIR) brightness and of the star formation rate. While the temperature of the warm dust at high brightness is driven by young massive stars, evolved stellar populations appear to drive the temperature of the cold dust. Plane-parallel models of photon dominated regions (PDRs) fail to reproduce fully the [CII], [OI], and CO maps obtained in a first spectroscopic study of one 2'x2' subregion of M33, located on the inner, northern spiral arm and encompassing the HII region BCLMP302.

preprint2011arXiv

The Herschel M33 extended survey (HerM33es): PACS spectroscopy of the star forming region BCLMP 302

Context: The emission line of [CII] at 158 micron is one of the strongest cooling lines of the interstellar medium (ISM) in galaxies. Aims: Disentangling the relative contributions of the different ISM phases to [CII] emission, is a major topic of the HerM33es program, a Herschel key project to study the ISM in the nearby spiral galaxy M33. Methods: Using PACS, we have mapped the emission of [CII] 158 micron, [OI] 63 micron, and other FIR lines in a 2'x2' region of the northern spiral arm of M33, centered on the HII region BCLMP302. At the peak of H-alpha emission, we have observed in addition a velocity resolved [CII] spectrum using HIFI. We use scatterplots to compare these data with PACS 160 micron continuum maps, and with maps of CO and HI data, at a common resolution of 12 arcsec or 50 pc. Maps of H-alpha and 24 micron emission observed with Spitzer are used to estimate the SFR. We have created maps of the [CII] and [OI] 63 micron emission and detected [NII] 122 micron and NIII 57 micron at individual positions. Results: The [CII] line observed with HIFI is significantly broader than that of CO, and slightly blue-shifted. In addition, there is little spatial correlation between [CII] observed with PACS and CO over the mapped region. There is even less spatial correlation between [CII] and the atomic gas traced by HI. Detailed comparison of the observed intensities towards the HII region with models of photo ionization and photon dominated regions, confirms that a significant fraction, 20--30%, of the observed [CII] emission stems from the ionized gas and not from the molecular cloud. The gas heating efficiency, using the ratio between [CII] and the TIR as a proxy, varies between 0.07 and 1.5%, with the largest variations found outside the HII region.

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

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

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

preprint2010arXiv

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

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

preprint2010arXiv

CHESS, Chemical Herschel surveys of star forming regions:Peering into the protostellar shock L1157-B1

The outflow driven by the low-mass class 0 protostar L1157 is the prototype of the so-called chemically active outflows. The bright bowshock B1 in the southern outflow lobe is a privileged testbed of magneto-hydrodynamical (MHD) shock models, for which dynamical and chemical processes are strongly interdependent. We present the first results of the unbiased spectral survey of the L1157-B1 bowshock, obtained in the framework of the key program "Chemical Herschel Surveys of Star Forming Regions" (CHESS). The main aim is to trace the warm and chemically enriched gas and to infer the excitation conditions in the shock region. The CO 5-4 and H2O lines have been detected at high-spectral resolution in the unbiased spectral survey of the HIFI-Band 1b spectral window (555-636 GHz), presented by Codella et al. in this volume. Complementary ground-based observations in the submm window help establish the origin of the emission detected in the main-beam of HIFI, and the physical conditions in the shock.}{Both lines exhibit broad wings, which extend to velocities much higher than reported up to now. We find that the molecular emission arises from two regions with distinct physical conditions: an extended, warm (100K), dense (3e5 cm-3) component at low-velocity, which dominates the water line flux in Band~1; a secondary component in a small region of B1 (a few arcsec) associated with high-velocity, hot (> 400 K) gas of moderate density ((1.0-3.0)e4 cm-3), which appears to dominate the flux of the water line at 179mu observed with PACS. The water abundance is enhanced by two orders of magnitude between the low- and the high-velocity component, from 8e-7 up to 8e-5. The properties of the high-velocity component agree well with the predictions of steady-state C-shock models.

preprint2010arXiv

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

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

preprint2010arXiv

Detection of interstellar oxidaniumyl: abundant H2O+ towards the star-forming regions DR21, Sgr B2, and NGC6334

We identify a prominent absorption feature at 1115 GHz, detected in first HIFI spectra towards high-mass star-forming regions, and interpret its astrophysical origin. The characteristic hyperfine pattern of the H2O+ ground-state rotational transition, and the lack of other known low-energy transitions in this frequency range, identifies the feature as H2O+ absorption against the dust continuum background and allows us to derive the velocity profile of the absorbing gas. By comparing this velocity profile with velocity profiles of other tracers in the DR21 star-forming region, we constrain the frequency of the transition and the conditions for its formation. In DR21, the velocity distribution of H2O+ matches that of the [CII] line at 158μ\m and of OH cm-wave absorption, both stemming from the hot and dense clump surfaces facing the HII-region and dynamically affected by the blister outflow. Diffuse foreground gas dominates the absorption towards Sgr B2. The integrated intensity of the absorption line allows us to derive lower limits to the H2O+ column density of 7.2e12 cm^-2 in NGC 6334, 2.3e13 cm^-2 in DR21, and 1.1e15 cm^-2 in Sgr B2.

preprint2010arXiv

Excitation of the molecular gas in the nuclear region of M82

We present high resolution HIFI spectroscopy of the nucleus of the archetypical starburst galaxy M82. Six 12CO lines, 2 13CO lines and 4 fine-structure lines are detected. Besides showing the effects of the overall velocity structure of the nuclear region, the line profiles also indicate the presence of multiple components with different optical depths, temperatures and densities in the observing beam. The data have been interpreted using a grid of PDR models. It is found that the majority of the molecular gas is in low density (n=10^3.5 cm^-3) clouds, with column densities of N_H=10^21.5 cm^-2 and a relatively low UV radiation field (GO = 10^2). The remaining gas is predominantly found in clouds with higher densities (n=10^5 cm^-3) and radiation fields (GO = 10^2.75), but somewhat lower column densities (N_H=10^21.2 cm^-2). The highest J CO lines are dominated by a small (1% relative surface filling) component, with an even higher density (n=10^6 cm^-3) and UV field (GO = 10^3.25). These results show the strength of multi-component modeling for the interpretation of the integrated properties of galaxies.

preprint2010arXiv

First results on Martian carbon monoxide from Herschel/HIFI observations

We report on the initial analysis of Herschel/HIFI carbon monoxide (CO) observations of the Martian atmosphere performed between 11 and 16 April 2010. We selected the (7-6) rotational transitions of the isotopes ^{13}CO at 771 GHz and C^{18}O at 768 GHz in order to retrieve the mean vertical profile of temperature and the mean volume mixing ratio of carbon monoxide. The derived temperature profile agrees within less than 5 K with general circulation model (GCM) predictions up to an altitude of 45 km, however, show about 12-15 K lower values at 60 km. The CO mixing ratio was determined as 980 \pm 150 ppm, in agreement with the 900 ppm derived from Herschel/SPIRE observations in November 2009.

preprint2010arXiv

Gas morphology and energetics at the surface of PDRs: new insights with Herschel observations of NGC 7023

We investigate the physics and chemistry of the gas and dust in dense photon-dominated regions (PDRs), along with their dependence on the illuminating UV field. Using Herschel-HIFI observations, we study the gas energetics in NGC 7023 in relation to the morphology of this nebula. NGC 7023 is the prototype of a PDR illuminated by a B2V star and is one of the key targets of Herschel. Our approach consists in determining the energetics of the region by combining the information carried by the mid-IR spectrum (extinction by classical grains, emission from very small dust particles) with that of the main gas coolant lines. In this letter, we discuss more specifically the intensity and line profile of the 158 micron (1901 GHz) [CII] line measured by HIFI and provide information on the emitting gas. We show that both the [CII] emission and the mid-IR emission from polycyclic aromatic hydrocarbons (PAHs) arise from the regions located in the transition zone between atomic and molecular gas. Using the Meudon PDR code and a simple transfer model, we find good agreement between the calculated and observed [CII] intensities. HIFI observations of NGC 7023 provide the opportunity to constrain the energetics at the surface of PDRs. Future work will include analysis of the main coolant line [OI] and use of a new PDR model that includes PAH-related species.

preprint2010arXiv

Herschel observations of water vapour in Markarian 231

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

preprint2010arXiv

HIFI observations of warm gas in DR21: Shock versus radiative heating

The molecular gas in the DR21 massive star formation region is known to be affected by the strong UV field from the central star cluster and by a fast outflow creating a bright shock. The relative contribution of both heating mechanisms is the matter of a long debate. By better sampling the excitation ladder of various tracers we provide a quantitative distinction between the different heating mechanisms. HIFI observations of mid-J transitions of CO and HCO+ isotopes allow us to bridge the gap in excitation energies between observations from the ground, characterizing the cooler gas, and existing ISO LWS spectra, constraining the properties of the hot gas. Comparing the detailed line profiles allows to identify the physical structure of the different components. In spite of the known shock-excitation of H2 and the clearly visible strong outflow, we find that the emission of all lines up to > 2 THz can be explained by purely radiative heating of the material. However, the new Herschel/HIFI observations reveal two types of excitation conditions. We find hot and dense clumps close to the central cluster, probably dynamically affected by the outflow, and a more widespread distribution of cooler, but nevertheless dense, molecular clumps.

preprint2010arXiv

Molecular and Atomic Gas in the Local Group Galaxy M33

We present high resolution large scale observations of the molecular and atomic gas in the Local Group Galaxy M33. The observations were carried out using the HERA at the 30m IRAM telescope in the CO(2-1) line achieving a resolution of 12"x2.6 km/s, enabling individual GMCs to be resolved. The observed region mainly along the major axis out to a radius of 8.5 kpc, and covers the strip observed with HIFI/PACS Spectrometers as part of the HERM33ES Herschel key program. The achieved sensitivity in main beam temperature is 20-50 mK at 2.6 km/s velocity resolution. The CO(2-1) luminosity of the observed region is 1.7\pm0.1x10^7 Kkm/s pc^2, corresponding to H2 masses of 1.9x10^8 Msun (including He), calculated with a NH2/ICO twice the Galactic value due to the half-solar metallicity of M33. HI 21 cm VLA archive observations were reduced and the mosaic was imaged and cleaned using the multi-scale task in CASA, yielding a series of datacubes with resolutions ranging from 5" to 25". The HI mass within a radius of 8.5 kpc is estimated to be 1.4x10^9 Msun. The azimuthally averaged CO surface brightness decreases exponentially with a scale length of 1.9\pm0.1 kpc whereas the atomic gas surface density is constant at Sigma_HI=6\pm2 Msun/pc^2 deprojected to face-on. The central kiloparsec H_2 surface density is Sigma_H2=8.5\pm0.2 Msun/pc^2. The star formation rate per unit molecular gas (SF Efficiency, the rate of transformation of molecular gas into stars), as traced by the ratio of CO to Halpha and FIR brightness, is constant with radius. The SFE appears 2-4 times greater than of large spiral galaxies. A morphological comparison of molecular and atomic gas with tracers of star formation shows good agreement between these maps both in terms of peaks and holes. A few exceptions are noted. Several spectra, including those of a molecular cloud situated more than 8 kpc from the galaxy center, are presented.

preprint2010arXiv

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

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

preprint2010arXiv

Photon Dominated Regions in NGC 3603

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

preprint2010arXiv

Properties of compact 250 μm emission and HII regions in M33 (HERM33ES)

Within the framework of the HERM33ES Key Project, using the high resolution and sensitivity of the Herschel photometric data, we study the compact emission in the Local Group spiral galaxy M33 to investigate the nature of the compact SPIRE emission sources. We extracted a catalogue of sources at 250um in order to investigate the nature of this compact emission. Taking advantage of the unprecedented Herschel resolution at these wavelengths, we also focus on a more precise study of some striking Halpha shells in the northern part of the galaxy. We present a catalogue of 159 compact emission sources in M33 identified by SExtractor in the 250um SPIRE band that is the one that provides the best spatial resolution. We also measured fluxes at 24um and Halpha for those 159 extracted sources. The morphological study of the shells also benefits from a multiwavelength approach including Halpha, far-UV from GALEX, and infrared from both Spitzer IRAC 8um and MIPS 24um in order to make comparisons. For the 159 compact sources selected at 250um, we find a very strong Pearson correlation coefficient with the MIPS 24um emission (r24 = 0.94) and a rather strong correlation with the Halpha emission, although with more scatter (rHa = 0.83). The morphological study of the Halpha shells shows a displacement between far-ultraviolet, Halpha, and the SPIRE bands. The cool dust emission from SPIRE clearly delineates the Halpha shell structures. The very strong link between the 250um compact emission and the 24um and Halpha emissions, by recovering the star formation rate from standard recipes for HII regions, allows us to provide star formation rate calibrations based on the 250um compact emission alone. The different locations of the Halpha and far-ultraviolet emissions with respect to the SPIRE cool dust emission leads to a dynamical age of a few Myr for the Halpha shells and the associated cool dust.

preprint2010arXiv

The origin of the [C II] emission in the S140 PDRs - new insights from HIFI

Using Herschel's HIFI instrument we have observed [C II] along a cut through S140 and high-J transitions of CO and HCO+ at two positions on the cut, corresponding to the externally irradiated ionization front and the embedded massive star forming core IRS1. The HIFI data were combined with available ground-based observations and modeled using the KOSMA-tau model for photon dominated regions. Here we derive the physical conditions in S140 and in particular the origin of [C II] emission around IRS1. We identify three distinct regions of [C II] emission from the cut, one close to the embedded source IRS1, one associated with the ionization front and one further into the cloud. The line emission can be understood in terms of a clumpy model of photon-dominated regions. At the position of IRS1, we identify at least two distinct components contributing to the [C II] emission, one of them a small, hot component, which can possibly be identified with the irradiated outflow walls. This is consistent with the fact that the [C II] peak at IRS1 coincides with shocked H2 emission at the edges of the outflow cavity. We note that previously available observations of IRS1 can be well reproduced by a single-component KOSMA-tau model. Thus it is HIFI's unprecedented spatial and spectral resolution, as well as its sensitivity which has allowed us to uncover an additional hot gas component in the S140 region.

preprint2009arXiv

HERACLES: The HERA CO-Line Extragalactic Survey

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

preprint2007arXiv

Gravitational stability in the disk of M51

Star formation laws, like i.e. the Schmidt law relating star formation rate and total gas density, have been studied in several spiral galaxies but the underlying physics are not yet well understood. M51, as a nearby face-on, grand design spiral galaxy studied in many line transitions, is an ideal target to study the connection between physical conditions of the gas and star formation activity. In this contribution we combine molecular, atomic, total gas and stellar surface densities and study the gravitational stability of the gas (Schuster et al.2007, Hitschfeld et al. in prep.). From our IRAM-30m 12 CO2-1 map and complementary HI-, Radio Continuum- and ACS-HST B-band-data we derive maps of the total gas density and the stellar surface density to study the gravitational stability of the gas via the Toomre Q parameter.