Researcher profile

P. Gratier

P. Gratier contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
33works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

33 published item(s)

preprint2023arXiv

Constraints on the non-thermal desorption of methanol in the cold core LDN 429-C

Cold cores are an early step of star formation, characterized by densities > 10$^4$ cm$^{-3}$, low temperatures (< 15 K), and very low external UV radiation. We investigate the physico-chemical processes at play to tracing the origin of molecules that are predominantly formed via reactions on dust grain surfaces. We observed the cold core LDN 429-C with the NOEMA interferometer and the IRAM 30m single dish telescope in order to obtain the gas-phase abundances of key species, including CO and CH$_3$OH. Comparing the observed gas phase of methanol to its solid phase previously observed with Spitzer allows us to put quantitative constraints on the efficiency of the non-thermal desorption of this species. With physical parameters determined from available Herschel data, we computed abundance maps of 11 detected molecules with a non-local thermal equilibrium radiative transfer model. These observations allowed us to probe the molecular abundances as a function of density and visual extinction, with the variation in temperature being restrained between 12 and 18 K. We then compared the observed abundances to the predictions of the Nautilus astrochemical model. We find that all molecules have lower abundances at high densities and visual extinctions with respect to lower density regions, except for methanol. Comparing these observations with a grid of chemical models based on the local physical conditions, we were able to reproduce these observations, allowing only the parameter time to vary. Comparing the observed gas-phase abundance of methanol with previous measurements of the methanol ice, we estimate a non-thermal desorption efficiency between 0.002% and 0.09%, increasing with density. The apparent increase in the desorption efficiency cannot be reproduced by our model unless the yield of cosmic-ray sputtering is altered due to the ice composition varying as a function of density.

preprint2022arXiv

From the Circumnuclear Disk in the Galactic Center to thick, obscuring tori of AGNs -- Modelling the molecular emission of a parsec-scale torus as found in NGC1068

The accretion rates needed to fuel the central black hole in a galaxy can be achieved via viscous torques in thick disks and rings, which can be resolved by millimetre interferometry within the inner ~20pc of the active galaxy NGC1068 at comparable scales and sensitivity to single dish observations of the Circumnuclear Disk (CND) in the Galactic Center. To interpret observations of these regions and determine the physical properties of their gas distribution, we present a modelling effort that includes (i) a simple dynamical simulations involving partially inelastic collisions between disk gas clouds, (ii) an analytical model of a turbulent clumpy gas disk calibrated by the dynamical model and observations, (iii) local turbulent and cosmic ray gas heating and cooling via H2O, H2, and CO emission, and (iv) determination of the molecular abundances. We also consider photodissociation regions (PDR) where gas is directly illuminated by the central engine. We compare the resulting model datacubes of the CO, HCN, HCO+, and CS brightness temperatures to available observations. In both cases the kinematics can be explained by one or two clouds colliding with a pre-existing ring, in a prograde sense for the CND and retrograde for NGC1068. And, with only dense disk clouds, the line fluxes can be reproduced to within a factor of about two. To avoid self-absorption of the intercloud medium, turbulent heating at the largest scales, comparable to the disk height, has to be decreased by a factor of 50-200. Our models indicate that turbulent mechanical energy input is the dominant gas heating mechanism within the thick gas disks. In N1068, while the bulk of the AGN X-ray radiation is absorbed in a layer of Compton-thick gas inside the dust sublimation radius, the optical/UV radiation may enhance the molecular line emission from photodissociation regions by ~50% at the inner edge of the gas ring.

preprint2022arXiv

Predicting HCN, HCO + , multi-transition CO, and dust emission of star-forming galaxies: Constraining the properties of resolved gas and dust disks of local spiral galaxies

The ISM is a turbulent, multi-phase, and multi-scale medium following scaling relations. Analytical models of galactic gaseous disks need to take into account the multi-scale and multi-phase nature of the interstellar medium. They can be described as clumpy star-forming accretion disks in vertical hydrostatic equilibrium, with the mid-plane pressure balancing the gravity of the gaseous and stellar disk. ISM turbulence is taken into account by applying Galactic scaling relations to the cold atomic and molecular gas phases. Turbulence is maintained through energy injection by supernovae. With the determination of the gas mass fraction at a given spatial scale, the equilibrium gas temperature between turbulent heating and line cooling, the molecular abundances, and the molecular line emission can be calculated. The resulting model radial profiles of IR, H{\sc i}, CO, HCN, and HCO$^+$ emission are compared to THINGS, HERACLES, EMPIRE, SINGS, and GALEX observations of 17 local spiral galaxies. The Toomre parameter, which measures the stability against star formation (cloud collapse), exceeds unity in the inner disk of a significant number of galaxies. In two galaxies it also exceeds unity in the outer disk. Therefore, in spiral galaxies $Q_{\rm tot}=1$ is not ubiquitous. The model gas velocity dispersion is consistent with the observed H{\sc i} velocity dispersion where available. Within our model HCN and HCO$^+$ is already detectable in relatively low-density gas ($\sim 1000$~cm$^{-3}$). CO and HCN conversion factors and molecular gas depletion time were derived. Both conversion factors are consistent with values found in the literature. Whereas in the massive galaxies the viscous timescale greatly exceeds the star formation timescale, the viscous timescale is smaller than the star formation timescale within $\rm{R}~\sim~2~\rm{R}_{\rm d}$, the disk scale length, in the low-mass galaxies.

preprint2020arXiv

Influence of galactic arm scale dynamics on the molecular composition of the cold and dense ISM III. Elemental depletion and shortcomings of the current physico-chemical models

We present a study of the elemental depletion in the interstellar medium. We combined the results of a Galatic model describing the gas physical conditions during the formation of dense cores with a full-gas-grain chemical model. During the transition between diffuse and dense medium, the reservoirs of elements, initially atomic in the gas, are gradually depleted on dust grains (with a phase of neutralisation for those which are ions). This process becomes efficient when the density is larger than 100~cm$^{-3}$. If the dense material goes back into diffuse conditions, these elements are brought back in the gas-phase because of photo-dissociations of the molecules on the ices followed by thermal desorption from the grains. Nothing remains on the grains for densities below 10~cm$^{-3}$ or in the gas-phase in a molecular form. One exception is chlorine, which is efficiently converted at low density. Our current gas-grain chemical model is not able to reproduce the depletion of atoms observed in the diffuse medium except for Cl which gas abundance follows the observed one in medium with densities smaller than 10~cm$^{-3}$. This is an indication that crucial processes (involving maybe chemisorption and/or ice irradiation profoundly modifying the nature of the ices) are missing.

preprint2016arXiv

A New Reference Chemical Composition for TMC-1

Recent detections of complex organic molecules in dark clouds have rekindled interest in the astrochemical modeling of these environments. Because of its relative closeness and rich molecular complexity, TMC-1 has been extensively observed to study the chemical processes taking place in dark clouds. We use local thermodynamical equilibrium radiative transfer modeling coupled with a Bayesian statistical method which takes into account outliers to analyze the data from the Nobeyama spectral survey of TMC-1 between 8 and 50 GHz. We compute the abundance relative to molecular hydrogen of 57 molecules, including 19 isotopologues in TMC-1 along with their associated uncertainty. The new results are in general agreement with previous abundance determination from Ohishi & Kaifu and the values reported in the review from Agundez & Wakelam. However, in some cases, large opacity and low signal to noise effects allow only upper or lower limits to be derived, respectively.

preprint2016arXiv

Inferring asymmetric limb cloudiness on exoplanets from transit light curves

Clouds have been shown to be present in many exoplanetary atmospheres. Cloud formation modeling predicts considerable inhomogeneities of cloud cover, consistent with optical phase curve observations. However, optical phase curves cannot resolve some existing degeneracies between cloud location and cloud optical properties. We present a conceptually simple technique to detect inhomogeneous cloud cover on exoplanets. Such an inhomogeneous cloud cover produces an asymmetric primary transit of the planet in front of the host star. Asymmetric transits produce characteristic residuals compared to a standard symmetric model. Furthermore, bisector spans can be used to determine asymmetries in the transit light curve. We apply a model of asymmetric transits to the light curves of HAT-P-7b, Kepler-7b and HD209458b and search for possible cloud signatures. The nearly uninterrupted Kepler photometry is particularly well-suited for this method since it allows for a very high time resolution. We do not find any statistically sound cloud signature in the data of the considered planets. For HAT-P-7b, a tentative detection of an asymmetric cloud cover is found, consistent with analysis of the optical phase curve. Based on Bayesian probability arguments, a symmetric model with an offset in the transit ephemeris remains, however, the most viable model. Still, this work demonstrates that for suitable targets, namely low-gravity planets around bright stars, the method can be used to constrain cloud cover characteristics and is thus a helpful additional tool to study exoplanetary atmospheres.

preprint2015arXiv

Inferring heat recirculation and albedo for exoplanetary atmospheres: Comparing optical phase curves and secondary eclipse data

Basic atmospheric properties such as albedo and heat redistribution between day and nightside have been inferred for a number of planets using observations of secondary eclipses and thermal phase curves. Optical phase curves have not yet been used to constrain these atmospheric properties consistently. We re-model previously published phase curves of CoRoT-1b, TrES-2b and HAT-P-7b and infer albedos and recirculation efficiencies. These are then compared to previous estimates based on secondary eclipse data. We use a physically consistent model to construct optical phase curves. This model takes Lambertian reflection, thermal emission, ellipsoidal variations and Doppler boosting into account. CoRoT-1b shows a non-negligible scattering albedo (0.11<AS<0.3 at 95% confidence) as well as small day-night temperature contrasts, indicative of moderate to high re-distribution of energy between dayside and nightside. These values are contrary to previous secondary eclipse and phase curve analyses. In the case of HAT-P-7b, model results suggest relatively high scattering albedo (AS=0.3). This confirms previous phase curve analysis, however, it is in slight contradiction to values inferred from secondary eclipse data. For TrES-2b, both approaches yield very similar estimates of albedo and heat recirculation. Discrepancies between recirculation and albedo values as inferred from secondary eclipse and optical phase curve analyses might be interpreted as a hint that optical and IR observations probe different atmospheric layers, and hence temperatures.

preprint2015arXiv

Spatially resolved l-c3h+ emission in the horsehead photodissociation region: Further evidence for a top-down hydrocarbon chemistry

Small hydrocarbons, such as C2H, C3H and C3H2 are more abundant in photo-dissociation regions (PDRs) than expected based on gas-phase chemical models. To explore the hydrocarbon chemistry further, we observed a key intermediate species, the hydrocarbon ion l-C3H+, in the Horsehead PDR with the Plateau de Bure Interferometer at high-angular resolution (6&#39;&#39;). We compare with previous observations of C2H and c-C3H2 at similar angular resolution and new gas-phase chemical model predictions to constrain the dominant formation mechanisms of small hydrocarbons in low-UV flux PDRs. We find that, at the peak of the HCO emission (PDR position), the measured l-C3H+, C2H and c-C3H2 abundances are consistent with current gas-phase model predictions. However, in the first PDR layers, at the 7.7 mum PAH band emission peak, which are more exposed to the radiation field and where the density is lower, the C2H and c-C3H2 abundances are underestimated by an order of magnitude. At this position, the l-C3H+ abundance is also underpredicted by the model but only by a factor of a few. In addition, contrary to the model predictions, l-C3H+ peaks further out in the PDR than the other hydrocarbons, C2H and c-C3H2. This cannot be explained by an excitation effect. Current gas-phase photochemical models thus cannot explain the observed abundances of hydrocarbons, in particular in the first PDR layers. Our observations are consistent with a top-down hydrocarbon chemistry, in which large polyatomic molecules or small carbonaceous grains are photo-destroyed into smaller hydrocarbon molecules/precursors.

preprint2015arXiv

The 2014 KIDA network for interstellar chemistry

Chemical models used to study the chemical composition of the gas and the ices in the interstellar medium are based on a network of chemical reactions and associated rate coefficients. These reactions and rate coefficients are partially compiled from data in the literature, when available. We present in this paper kida.uva.2014, a new updated version of the kida.uva public gas-phase network first released in 2012. In addition to a description of the many specific updates, we illustrate changes in the predicted abundances of molecules for cold dense cloud conditions as compared with the results of the previous version of our network, kida.uva.2011.

preprint2015arXiv

Widespread Galactic CF+ absorption: detection toward W49 with the Plateau de Bure Interferometer

To study the usefulness of \CFP\ as a tracer of the regions where C\p\ and \HH\ coexist in the interstellar medium. We used the Plateau de Bure Interferometer to synthesize \CFP\ J=1-0 absorption at 102.6 GHz toward the core of the distant HII region W49N at l = 43.2\degr, b=0.0\degr, and we modeled the fluorine chemistry in diffuse/translucent molecular gas. We detected \CFP\ absorption over a broad range of velocity showing that \CFP\ is widespread in the \HH-bearing Galactic disk gas. Originally detected in dense gas in the Orion Bar and Horsehead PDR, \CFP\ was subsequently detected in absorption from diffuse and translucent clouds seen toward \bll\ and 3C111. Here we showed that \CFP\ is distributed throughout the diffuse and translucent molecular disk gas with N(\CFP)/N(\HH) $= 1.5-2.0\times10^{-10}$, increasing to N(\CFP)/N(\HH) $= 3.5\times10^{-10}$ in one cloud at 39 \kms\ having higher N(\HH) $\approx 3\times10^{21}\pcc$. Models of the fluorine chemistry reproduce the observed column densities and relative abundance of HF, from which \CFP\ forms, but generally overpredict the the column density of \CFP\ by factors of 1.4-4. We show that a free space photodissociation rate $Γ\ga 10^{-9}\ps$, comparable to that of CH, might account for much of the discrepancy but a recent calculation finds a value about ten times smaller. In the heavily blended and kinematically complex spectra seen toward W49, \CFP\ absorption primarily traces the peaks of the \HH\ distribution.

preprint2014arXiv

Dense molecular globulettes and the dust arc towards the runaway O star AE Aur (HD 34078)

Some runaway stars are known to display IR arc-like structures around them, resulting from their interaction with surrounding interstellar material. The properties of these features as well as the processes involved in their formation are still poorly understood. We aim at understanding the physical mechanisms that shapes the dust arc observed near the runaway O star AEAur (HD34078). We obtained and analyzed a high spatial resolution map of the CO(1-0) emission that is centered on HD34078, and that combines data from both the IRAM interferometer and 30m single-dish antenna. The line of sight towards HD34078 intersects the outer part of one of the detected globulettes, which accounts for both the properties of diffuse UV light observed in the field and the numerous molecular absorption lines detected in HD34078&#39;s spectra, including those from highly excited H2 . Their modeled distance from the star is compatible with the fact that they lie on the 3D paraboloid which fits the arc detected in the 24 μm Spitzer image. Four other compact CO globulettes are detected in the mapped area. These globulettes have a high density and linewidth, and are strongly pressure-confined or transient. The good spatial correlation between the CO globulettes and the IR arc suggests that they result from the interaction of the radiation and wind emitted by HD 34078 with the ambient gas. However, the details of this interaction remain unclear. A wind mass loss rate significantly larger than the value inferred from UV lines is favored by the large IR arc size, but does not easily explain the low velocity of the CO globulettes. The effect of radiation pressure on dust grains also meets several issues in explaining the observations. Further observational and theoretical work is needed to fully elucidate the processes shaping the gas and dust in bow shocks around runaway O stars. (Abridged)

preprint2014arXiv

Modeling Complex Organic Molecules in dense regions: Eley-Rideal and complex induced reaction

Recent observations have revealed the existence of Complex Organic Molecules (COMs) in cold dense cores and prestellar cores. The presence of these molecules in such cold conditions is not well understood and remains a matter of debate since the previously proposed &#34;warm- up&#34; scenario cannot explain these observations. In this article, we study the effect of Eley- Rideal and complex induced reaction mechanisms of gas-phase carbon atoms with the main ice components of dust grains on the formation of COMs in cold and dense regions. Based on recent experiments we use a low value for the chemical desorption efficiency (which was previously invoked to explain the observed COM abundances). We show that our introduced mechanisms are efficient enough to produce a large amount of complex organic molecules in the gas-phase at temperatures as low as 10K.

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&#34; 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.

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

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

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 IRAM-30m line survey of the Horsehead PDR: III. High abundance of complex (iso-)nitrile molecules in UV-illuminated gas

Complex (iso-)nitrile molecules, such as CH3CN and HC3N, are relatively easily detected in our Galaxy and in other galaxies. We constrain their chemistry through observations of two positions in the Horsehead edge: the photo-dissociation region (PDR) and the dense, cold, and UV-shielded core just behind it. We systematically searched for lines of CH3CN, HC3N, C3N, and some of their isomers in our sensitive unbiased line survey at 3, 2, and 1mm. We derived column densities and abundances through Bayesian analysis using a large velocity gradient radiative transfer model. We report the first clear detection of CH3NC at millimeter wavelength. We detected 17 lines of CH3CN at the PDR and 6 at the dense core position, and we resolved its hyperfine structure for 3 lines. We detected 4 lines of HC3N, and C3N is clearly detected at the PDR position. We computed new electron collisional rate coefficients for CH3CN, and we found that including electron excitation reduces the derived column density by 40% at the PDR position. While CH3CN is 30 times more abundant in the PDR than in the dense core, HC3N has similar abundance at both positions. The isomeric ratio CH3NC/CH3CN is 0.15+-0.02. In the case of CH3CN, pure gas phase chemistry cannot reproduce the amount of CH3CN observed in the UV-illuminated gas. We propose that CH3CN gas phase abundance is enhanced when ice mantles of grains are destroyed through photo-desorption or thermal-evaporation in PDRs, and through sputtering in shocks. (abridged)

preprint2013arXiv

The IRAM-30m line survey of the Horsehead PDR: IV. Comparative chemistry of H2CO and CH3OH

Aims. We investigate the dominant formation mechanism of H2CO and CH3OH in the Horsehead PDR and its associated dense core. Methods. We performed deep integrations of several H2CO and CH3OH lines at two positions in the Horsehead, namely the PDR and dense core, with the IRAM-30m telescope. In addition, we observed one H2CO higher frequency line with the CSO telescope at both positions. We determine the H2CO and CH3OH column densities and abundances from the single-dish observations complemented with IRAM-PdBI high-angular resolution maps (6&#34;) of both species. We compare the observed abundances with PDR models including either pure gas-phase chemistry or both gas-phase and grain surface chemistry. Results. We derive CH3OH abundances relative to total number of hydrogen atoms of ~1.2e-10 and ~2.3e-10 in the PDR and dense core positions, respectively. These abundances are similar to the inferred H2CO abundance in both positions (~2e-10). We find an abundance ratio H2CO/CH3OH of ~2 in the PDR and ~1 in the dense core. Pure gas-phase models cannot reproduce the observed abundances of either H2CO or CH3OH at the PDR position. Both species are therefore formed on the surface of dust grains and are subsequently photodesorbed into the gas-phase at this position. At the dense core, on the other hand, photodesorption of ices is needed to explain the observed abundance of CH3OH, while a pure gas-phase model can reproduce the observed H2CO abundance. The high-resolution observations show that CH3OH is depleted onto grains at the dense core. CH3OH is thus present in an envelope around this position, while H2CO is present in both the envelope and the dense core itself. Conclusions. Photodesorption is an efficient mechanism to release complex molecules in low FUV-illuminated PDRs, where thermal desorption of ice mantles is ineffective.

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

A detailed view of a Molecular Cloud in the far outer disk of M33

The amount of H$_2$ present in spiral galaxies remains uncertain, particularly in the dim outer regions and in low-metallicity environments. We present high-resolution CO(1--0) observations with the Plateau de Bure interferometer of the most distant molecular cloud in the local group galaxy M 33. The cloud is a single entity rather than a set of smaller clouds within the broad beam of the original single-dish observations. The interferometer and single-dish fluxes are very similar and the line widths are indistinguishable, despite the difference in beamsize. At a spatial resolution of 10 pc, beyond the optical radius of the M 33, the CO brightness temperature reaches 2.4 Kelvins. A virial mass estimate for the cloud yields a mass of $4.3 \times 10^4$ \msun and a ratio $\ratio \simeq 3.5 \times 10^{20} \Xunit$. While no velocity gradient is seen where the emission is strong, the velocity is redshifted to the extreme SW and blue-shifted to the far NE. If the orientation of the cloud is along the plane of the disk (i.e. not perpendicular), then these velocities correspond to slow infall or accretion. The rather modest infall rate would be about $2 \times 10^{-4}$\moyr.

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 &#39;cool&#39; 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 &#39;cool&#39; dust grains (~15 K) while the rest 41% is produced by &#39;warm&#39; 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 &#39;cool&#39; 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

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&#34;$ 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&#34;$ 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 hyperfine structure in the rotational spectrum of CF+

Context. CF+ has recently been detected in the Horsehead and Orion Bar photo-dissociation regions. The J=1-0 line in the Horsehead is double-peaked in contrast to other millimeter lines. The origin of this double-peak profile may be kinematic or spectroscopic. Aims. We investigate the effect of hyperfine interactions due to the fluorine nucleus in CF+ on the rotational transitions. Methods. We compute the fluorine spin rotation constant of CF+ using high-level quantum chemical methods and determine the relative positions and intensities of each hyperfine component. This information is used to fit the theoretical hyperfine components to the observed CF+ line profiles, thereby employing the hyperfine fitting method in GILDAS. Results. The fluorine spin rotation constant of CF+ is 229.2 kHz. This way, the double-peaked CF+ line profiles are well fitted by the hyperfine components predicted by the calculations. The unusually large hyperfine splitting of the CF+ line therefore explains the shape of the lines detected in the Horsehead nebula, without invoking intricate kinematics in the UV-illuminated gas.

preprint2012arXiv

The IRAM-30m line survey of the Horsehead PDR: I. CF+ as a tracer of C+ and a measure of the Fluorine abundance

C+ is a key species in the interstellar medium but its 158 μm fine structure line cannot be observed from ground-based telescopes. Current models of fluorine chemistry predict that CF+ is the second most important fluorine reservoir, in regions where C+ is abundant. We detected the J = 1-0 and J = 2-1 rotational lines of CF+ with high signal-to-noise ratio towards the PDR and dense core positions in the Horsehead. Using a rotational diagram analysis, we derive a column density of N(CF+) = (1.5 - 2.0) \times 10^12 cm^-2. Because of the simple fluorine chemistry, the CF+ column density is proportional to the fluorine abundance. We thus infer the fluorine gas-phase abundance to be F/H = (0.6 - 1.5) \times 10^-8. Photochemical models indicate that CF+ is found in the layers where C+ is abundant. The emission arises in the UV illuminated skin of the nebula, tracing the outermost cloud layers. Indeed, CF+ and C+ are the only species observed to date in the Horsehead with a double peaked line profile caused by kinematics. We therefore propose that CF+, which is detectable from the ground, can be used as a proxy of the C+ layers.

preprint2012arXiv

The IRAM-30m line survey of the Horsehead PDR: II. First detection of the l-C3H+ hydrocarbon cation

We present the first detection of the l-C3H+ hydrocarbon in the interstellar medium. The Horsehead WHISPER project, a millimeter unbiased line survey at two positions, namely the photo-dissociation region (PDR) and the nearby shielded core, revealed a consistent set of eight unidentified lines toward the PDR position. Six of them are detected with a signal-to-noise ratio from 6 to 19, while the two last ones are tentatively detected. Mostly noise appears at the same frequency toward the dense core, located less than 40&#34; away. We simultaneously fit 1) the rotational and centrifugal distortion constants of a linear rotor, and 2) the Gaussian line shapes located at the eight predicted frequencies. The observed lines can be accurately fitted with a linear rotor model, implying a 1Sigma ground electronic state. The deduced rotational constant value is Be= 11244.9512 +/- 0.0015 MHz, close to that of l-C3H. We thus associate the lines to the l-C3H+ hydrocarbon cation, which enables us to constrain the chemistry of small hydrocarbons. A rotational diagram is then used to infer the excitation temperature and the column density. We finally compare the abundance to the results of the Meudon PDR photochemical model.

preprint2011arXiv

Giant Molecular Clouds in the Local Group Galaxy M33

We present an analysis of the systematic CO(2-1) survey at 12&#34; 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

Star formation in M33 (HerM33es)

Within the key project &#34;Herschel M33 extended survey&#34; (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&#39;x2&#39; 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&#39;x2&#39; 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.

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

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&#34;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&#34; to 25&#34;. 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

The Molecular Interstellar Medium of the Local Group Dwarf NGC6822

Do molecular clouds collapse to form stars at the same rate in all environments? In large spiral galaxies, the rate of transformation of H2 into stars (hereafter SFE) varies little. However, the SFE in distant objects (z~1) is much higher than in the large spiral disks that dominate the local universe. Some small local group galaxies share at least some of the characteristics of intermediate-redshift objects, such as size or color. Recent work has suggested that the Star Formation Efficiency (SFE, defined as the SFRate per unit H2) in local Dwarf galaxies may be as high as in the distant objects. A fundamental difficulty in these studies is the independent measure of the H2 mass in metal-deficient environments. At 490 kpc, NGC6822 is an excellent choice for this study; it has been mapped in the CO(2-1) line using the multibeam receiver HERA on the 30 meter IRAM telescope, yielding the largest sample of giant molecular clouds (GMCs) in this galaxy. Despite the much lower metallicity, we find no clear difference in the properties of the GMCs in NGC 6822 and those in the Milky Way except lower CO luminosities for a given mass. Several independent methods indicate that the total H2 mass in NGC 6822 is about 5 x 10^6 Msun in the area we mapped and less than 10^7 Msun in the whole galaxy. This corresponds to a NH2/ICO ~ 4 x 10^{21} cm^-2 /(Kkm/s) over large scales, such as would be observed in distant objects, and half that in individual GMCs. No evidence was found for H2 without CO emission. Our simulations of the radiative transfer in clouds are entirely compatible with these NH2/ICO values. The SFE implied is a factor 5 - 10 higher than what is observed in large local universe spirals.