A&A 498, 771–783 (2009)DOI: 10.1051/0004-6361/200811496c○ ESO 2009Astronomy&AstrophysicsThe ionization fraction gradient across the Horsehead edge:an archetype for mo<strong>le</strong>cular clouds ⋆J. R. Goicoechea 1 ,J.Pety 2,3 , M. Gerin 3 , P. Hily-Blant 4 , and J. Le Bourlot 5<strong>tel</strong>-<strong>00726959</strong>, version 1 - 31 Aug 20121 Laboratorio de Astrofísica Mo<strong>le</strong>cular. Centro de Astrobiología. CSIC-INTA. Carretera de Ajalvir, Km 4. Torrejón de Ardoz,28850 Madrid, Spaine-mail: goicoechea@damir.iem.csic.es2 IRAM, 300 rue de la Piscine, 38406 Grenob<strong>le</strong> Cedex, Francee-mail: pety@iram.fr3 LERMA - LRA, UMR 8112, CNRS, Observatoire de Paris and Éco<strong>le</strong> Norma<strong>le</strong> Supérieure, 24 rue Lhomond, 75231 Paris, Francee-mail: maryvonne.gerin@lra.ens.fr4 Laboratoire d’Astrophysique, Observatoire de Grenob<strong>le</strong>, BP 53, 38041 Grenob<strong>le</strong> Cedex 09, Francee-mail: pierre.hilyblant@obs.ujf-grenob<strong>le</strong>.fr5 LUTH, UMR 8102 CNRS, Université Paris 7 and Observatoire de Paris, Place J. Janssen, 92195 Meudon, Francee-mail: Jacques.Lebourlot@obspm.frReceived 10 December 2008 / Accepted 11 February 2009ABSTRACTContext. The ionization fraction (i.e., the e<strong>le</strong>ctron abundance) plays a key ro<strong>le</strong> in the chemistry and dynamics of mo<strong>le</strong>cular clouds.Aims. We study the H 13 CO + ,DCO + and HOC + line emission towards the Horsehead, from the shielded core to the UV irradiatedcloud edge, i.e., the photodissociation region (PDR), as a template to investigate the ionization fraction gradient in mo<strong>le</strong>cular clouds.Methods. We analyze an IRAM Plateau de Bure Interferometer map of the H 13 CO + J = 1–0 line at a 6.8 ′′ × 4.7 ′′ resolution, comp<strong>le</strong>mentedwith IRAM-30 m H 13 CO + and DCO + higher-J line maps and new HOC + and CO + observations. We compare self-consistentlythe observed spatial distribution and line intensities with detai<strong>le</strong>d depth-dependent predictions of a PDR model coup<strong>le</strong>d with a nonlocalradiative transfer calculation. The chemical network includes deuterated species, 13 C fractionation reactions and HCO + /HOC +isomerization reactions. The ro<strong>le</strong> of neutral and charged PAHs in the cloud chemistry and ionization balance is investigated.Results. The detection of the HOC + reactive ion towards the Horsehead PDR proves the high ionization fraction of the outer UVirradiated regions, where we derive a low [HCO + ]/[HOC + ] ≃ 75–200 abundance ratio. In the absence of PAHs, we reproduce theobservations with gas-phase metal abundances, [Fe+Mg+...], lower than 4 × 10 −9 (with respect to H), and a cosmic-ray ionization rateof ζ = (5± 3)× 10 −17 s −1 . The inclusion of PAHs modifies the ionization fraction gradient and increases the required metal abundance.Conclusions. The ionization fraction in the Horsehead edge follows a steep gradient, with a sca<strong>le</strong> <strong>le</strong>ngth of ∼0.05 pc (or ∼25 ′′ ), from[e − ] ≃ 10 −4 (or n e ∼ 1–5 cm −3 )inthePDRtoafewtimes∼10 −9 in the core. PAH − anions play a ro<strong>le</strong> in the charge balance of thecold and neutral gas if substantial amounts of free PAHs are present ([PAH] > 10 −8 ).Key words. astrochemistry – ISM: clouds – radiative transfer – radio lines: ISM – ISM: mo<strong>le</strong>cu<strong>le</strong>s – ISM: abundances1. IntroductionThe e<strong>le</strong>ctron abundance ([e − ] = n e /n H ) plays a fundamentalro<strong>le</strong> in the chemistry and dynamics of inters<strong>tel</strong>lar gas.The degree of ionization determines the preponderance of ionneutralreactions, i.e., the main formation route for most chemicalspecies in mo<strong>le</strong>cular clouds (Herbst & K<strong>le</strong>mperer 1973;Oppenheimer & Dalgarno 1974). In addition, the ionizationfraction constrains the coupling of matter and magnetic fields,which drives the dissipation of turbu<strong>le</strong>nce and the transfer of angularmomentum, thus having crucial implications in protos<strong>tel</strong>larcollapse and accretion disks (e.g., Balbus & Haw<strong>le</strong>y 1991).High-angular resolution observations of inters<strong>tel</strong>lar cloudsreveal steep density, temperature and turbu<strong>le</strong>nce gradients aswell as sharp chemical variations. Accordingly, the e<strong>le</strong>ctron⋆ Based on observations obtained with the IRAM Plateau de Bure interferometerand 30 m <strong>tel</strong>escope. IRAM is supported by INSU/CNRS(France), MPG (Germany), and IGN (Spain).abundance should vary within a cloud depending on the relativeionizing sources and prevailing chemistry.Rotational line emission of mo<strong>le</strong>cular ions such as DCO +and HCO + have been traditionally used to estimate theionization fraction in mo<strong>le</strong>cular clouds because (i) theyare abundant and easily observab<strong>le</strong>; (ii) dissociative recombinationis their main destruction route, and thus theirabundances are roughly inversely proportional to the e<strong>le</strong>ctronabundance (e.g., Guélin et al. 1982; Wootten et al. 1982;de Boisanger et al. 1996; Williams et al. 1998; Caselli et al.1998; Maret & Bergin 2007; Hezareh et al. 2008). On the otherhand, the presence of reactive ions (species such as HOC + orCO + that react rapidly with H 2 ) is predicted to be a sensitiveindicator of high ionization fraction regions, e.g., the UV irradiatedcloud surfaces (e.g., Smith et al. 2002; Fuente et al. 2003).In order to constrain the ionization fraction gradient frommodels, the cloud chemistry and physics cannot be simplifiedmuch because the charge balance depends on parameters suchArtic<strong>le</strong> published by EDP Sciences
772 J. R. Goicoechea et al.: The ionization fraction gradient across the Horsehead edge: an archetype for mo<strong>le</strong>cular cloudsTab<strong>le</strong> 1. Observation parameters of the PdBI maps shown in Fig. 1.Mo<strong>le</strong>cu<strong>le</strong> Transition Frequency Instrument Config. Beam PA Vel. Resol. Int. Time a T sys Noise b,† Obs. dateGHz arcsec ◦ km s −1 hours K KH 13 CO + 1–0 86.754288 PdBI C & D 6.8 × 4.7 13 0.2 6.5 150 0.10 2006-07HCO 1 0,1 3/2, 2–0 0,0 1/2, 1 86.670760 PdBI C & D 6.7 × 4.4 16 0.2 6.5 150 0.09 2006-07a We observed a 7-field mosaic centered on the IR peak at α 2000 = 05 h 40 m 54.27 s , δ 2000 = −02 ◦ 28 ′ 00 ′′ (Abergel et al. 2003) with the followingoffsets: (−5.5 ′′ , −22.0 ′′ ), (5.5 ′′ , −22.0 ′′ ), (11.0 ′′ ,0.0 ′′ ), (0.0 ′′ ,0.0 ′′ ), (−11.0 ′′ ,0.0 ′′ ), (−5.5 ′′ , 22.0 ′′ )and(5.5 ′′ , 22.0 ′′ ). The total field-of-view is80.1 ′′ × 102.1 ′′ and the half power primary beam is 58.1 ′′ . The mosaic was Nyquist samp<strong>le</strong>d in declination at 3.4 mm and largely oversamp<strong>le</strong>d inright ascension. This maximizes the field of view along the PDR edge whi<strong>le</strong> the oversampling in the perpendicular direction eases the deconvolution.On-source time was computed as if the source was always observed with 6 antennae; b the noise values refer to the mosaic phase center(mosaic noise is inhomogeneous due to primary beam correction; it steeply increases at the mosaic edges).Tab<strong>le</strong> 2. Observation parameters of the IRAM-30 m observations.<strong>tel</strong>-<strong>00726959</strong>, version 1 - 31 Aug 2012Mo<strong>le</strong>cu<strong>le</strong> Transition Frequency Instrument F eff B eff Resol. Resol. Int. Time Noise † Observing Obs. dateGHz arcsec km s −1 hours K ModeHCO + J = 1–0 89.188523 30 m/A100 0.95 0.78 27.6 ′′ 0.20 4.7 0.02 ON-OFF 2008HOC + J = 1–0 89.487414 30 m/A100 0.95 0.78 27.5 ′′ 0.20 4.7 0.02 ON-OFF 2008CO + 2, 5/2–1, 3/2 236.062578 30 m/A230 0.91 0.52 10.4 ′′ 0.20 4.7 0.05 ON-OFF 2008H 13 CO + J = 1–0 86.754288 30 m/AB100 0.95 0.78 28.4 ′′ 0.20 2.6 0.10 OTF map 2006-07H 13 CO + J = 3–2 260.255339 30 m/HERA 0.90 0.46 13.5 ′′ 0.20 5.9 0.06 OTF map 2006DCO + J = 2–1 144.077289 30 m/CD150 0.93 0.69 18.0 ′′ 0.08 5.9 0.18 OTF map 2006DCO + J = 3–2 216.112582 30 m/HERA 0.90 0.52 11.4 ′′ 0.11 1.5 0.10 OTF map 2006† The noise (in T mb sca<strong>le</strong>) refers to the channel spacing obtained by averaging adjacent channels to the velocity resolution given in the tab<strong>le</strong>s.as the penetration of UV radiation, the cosmic-ray ionizationrate (ζ) and the abundance of key species (e.g., metals and PAH).Compared to other works, in this paper we determine theionization fraction gradient by direct comparison of H 13 CO +and DCO + high-angular resolution maps and HOC + pointed observations,with detai<strong>le</strong>d depth-dependent chemical and radiativetransfer models covering a broad range of cloud physicalconditions. Indeed, the observed field-of-view contains the famousHorsehead PDR (the UV illuminated edge of the cloud)and a dense and cold core discovered by us from its intenseDCO + line emission (Pety et al. 2007). Due to its simp<strong>le</strong> geometryand moderate distance (d ≃ 400 pc), the Horsehead PDRand associated core are good templates to study the steep gradientsexpected in mo<strong>le</strong>cular clouds (e.g., Pety et al. 2005, 2007;Goicoechea et al. 2006; Gerin et al. 2009).The paper is organized as follows. The observations are presentedin Sect. 2 and the models used to interpret them are describedinSect.3.The chemistry of H 13 CO + ,DCO + and HOC +(our observational probes of the ionization fraction) is analyzedin Sect. 4. In Sect. 5 we investigate the ro<strong>le</strong> of metals, PAHsand ζ on the e<strong>le</strong>ctron abundance determination. The main resultsand constrains are presented in Sect. 6 and discussed in Sect. 7.2. Observations2.1. Observations and data reductionTab<strong>le</strong>s 1 and 2 summarize the observation parameters of the dataobtained with the PdBI and the IRAM–30 m <strong>tel</strong>escope that weshall study in this work. The H 13 CO + J = 1–0 line emission mapwas first presented in Gerin et al. (2009). Frequency-switched,on-the-fly maps (OTF) obtained at the IRAM-30 m were used toproduce the short-spacings needed to comp<strong>le</strong>ment a 7-field mosaicacquired with the 6 PdBI antennae in the CD configuration(baseline <strong>le</strong>ngths from 24 to 176 m). Correlator backends wereused (VESPA for IRAM-30 m observations). The high angularresolution PdBI H 13 CO + J = 1–0 map comp<strong>le</strong>ments our previousH 13 CO + J = 3–2 and DCO + J = 2–1 and 3–2 maps takenwith the IRAM-30 m <strong>tel</strong>escope and first presented in Pety et al.(2007).In this work we present new IRAM–30 m deeper integrationsin the HOC + ,H 13 CO + and HCO + J = 1–0 lines, and an upperlimit for the CO + emission towards the PDR (defined here asthe HCO emission peak; Gerin et al. 2009). The position switchingobserving mode was used. The on-off cyc<strong>le</strong> duration was1mnandtheoff-position offsets were (Δα, Δδ) = (−100 ′′ , 0 ′′ ),i.e., the H ii region ionized by σOri and free of mo<strong>le</strong>cular gasemission. Position accuracy is estimated to be ∼3 ′′ for the 30 mdata and better than 0.5 ′′ for the PdBI data. The data processingwas done with the GILDAS 1 softwares (e.g., Pety 2005b).TheIRAM-30mdatawerefirstcalibratedtotheTA ∗ sca<strong>le</strong> usingthe chopper wheel method (Penzias & Burrus 1973), and finallyconverted to main beam temperatures T mb using the forwardand main beam efficiencies F eff and B eff displayed in Tab<strong>le</strong> 2(e.g., Greve et al. 1998). The amplitude accuracy for heterodyneobservations with the IRAM–30 m <strong>tel</strong>escope is ∼10%. PdBI dataand short-spacing data were merged before imaging and deconvolutionof the mosaic, using standard techniques of GILDAS andused in our previous works (see e.g., Pety et al. 2005).2.2. DCO + and H 13 CO + spatial distribution, HOC + detectionFigure 1 shows H 13 CO + J = 1–0, HCO 1 0,1 -0 0,0 (PdBI) andDCO + J = 2–1, 3–2 integrated line intensity maps (IRAM-30 m;Pety et al. 2007), as well as the aromatic infrared band emission(AIB, observed with ISOCAM, Abergel et al. 2003) thattraces the UV illuminated edge of the cloud, i.e., the PDR. TheDCO + emission is concentrated in a narrow, arclike structure ofdense and cold gas behind the PDR (Pety et al. 2007). Hence,it shows a very different spatial distribution than the emissionof “PDR tracers” such as C 2 H, C 4 H, c-C 3 H 2 (Pety et al. 2005),HCO radicals (Gerin et al. 2009), vibrationally excited H 2(Habart et al. 2005) or the AIB emission (Compiègne et al.2008). The H 13 CO + J = 1–0 emission follows the DCO +1 See http://www.iram.fr/IRAMFR/GILDAS
- Page 1 and 2:
UNIVERSITÉ PIERRE ET MARIE CURIEHA
- Page 3 and 4:
tel-00726959, version 1 - 31 Aug 20
- Page 5 and 6:
Table des matières1 Rapport de sou
- Page 7 and 8:
Rapport après soutenanceHabilitati
- Page 9 and 10:
Chapitre 2Curriculum vitaetel-00726
- Page 11 and 12:
2.7 ANIMATION ET DIFFUSION DE LA CU
- Page 13 and 14:
2.8 PARCOURS 131992-1993 ÉCOLE NOR
- Page 15 and 16:
Chapitre 3Copyright: IRAM/PdBIIntro
- Page 17 and 18:
tel-00726959, version 1 - 31 Aug 20
- Page 19 and 20:
4.2 ETUDES DIRECTES EN ÉMISSION 19
- Page 21 and 22:
4.4 LA LUMINOSITY CO PAR MOLÉCULE
- Page 23 and 24:
356 E. Falgarone et al.: Extreme ve
- Page 25 and 26:
358 E. Falgarone et al.: Extreme ve
- Page 27 and 28:
360 E. Falgarone et al.: Extreme ve
- Page 29 and 30:
362 E. Falgarone et al.: Extreme ve
- Page 31 and 32:
364 E. Falgarone et al.: Extreme ve
- Page 33 and 34:
366 E. Falgarone et al.: Extreme ve
- Page 35 and 36:
368 E. Falgarone et al.: Extreme ve
- Page 37 and 38:
A&A 518, A45 (2010)1001010⌠⌡ τ
- Page 39 and 40:
A&A 518, A45 (2010)If X HCO + is as
- Page 41 and 42:
A&A 518, A45 (2010)tel-00726959, ve
- Page 43 and 44:
A&A 518, A45 (2010)tel-00726959, ve
- Page 45 and 46:
A&A 518, A45 (2010)Table E.2. Data
- Page 47 and 48:
A&A 541, A58 (2012)tel-00726959, ve
- Page 49 and 50:
A&A 541, A58 (2012)tel-00726959, ve
- Page 51 and 52:
A&A 541, A58 (2012)0.08 0.10.20.12
- Page 53 and 54:
A&A 541, A58 (2012)b=240.11 0.1 0.1
- Page 55 and 56:
0.080.09 0.1A&A 541, A58 (2012)0.14
- Page 57 and 58:
A&A 541, A58 (2012)43˚00'0.3 0.44
- Page 59 and 60:
A&A 541, A58 (2012)68˚00'1.1 0.91.
- Page 61 and 62:
A&A 541, A58 (2012)α(J2000)-19H57m
- Page 63 and 64: A&A 541, A58 (2012)tel-00726959, ve
- Page 65 and 66: A&A 541, A58 (2012)Jenkins, E. B.,
- Page 67 and 68: A&A 541, A58 (2012)Galactic Latitud
- Page 69 and 70: tel-00726959, version 1 - 31 Aug 20
- Page 71 and 72: 5.2 UNE PHYSIQUE BIEN CONTRAINTE ET
- Page 73 and 74: 5.3 PERSPECTIVES : DES RELEVÉS DE
- Page 75 and 76: A&A 435, 885-899 (2005)DOI: 10.1051
- Page 77 and 78: tel-00726959, version 1 - 31 Aug 20
- Page 79 and 80: J. Pety et al.: Are PAHs precursors
- Page 81 and 82: J. Pety et al.: Are PAHs precursors
- Page 83 and 84: J. Pety et al.: Are PAHs precursors
- Page 85 and 86: J. Pety et al.: Are PAHs precursors
- Page 87 and 88: J. Pety et al.: Are PAHs precursors
- Page 89 and 90: J. Pety et al.: Are PAHs precursors
- Page 91 and 92: 566 J. R. Goicoechea et al.: Low su
- Page 93 and 94: 568 J. R. Goicoechea et al.: Low su
- Page 95 and 96: 570 J. R. Goicoechea et al.: Low su
- Page 97 and 98: 572 J. R. Goicoechea et al.: Low su
- Page 99 and 100: 574 J. R. Goicoechea et al.: Low su
- Page 101 and 102: 576 J. R. Goicoechea et al.: Low su
- Page 103 and 104: 578 J. R. Goicoechea et al.: Low su
- Page 105 and 106: 580 J. R. Goicoechea et al.: Low su
- Page 107 and 108: tel-00726959, version 1 - 31 Aug 20
- Page 109 and 110: J. R. Goicoechea et al.: Low sulfur
- Page 111 and 112: L42J. Pety et al.: Deuterium fracti
- Page 113: L44J. Pety et al.: Deuterium fracti
- Page 117 and 118: 774 J. R. Goicoechea et al.: The io
- Page 119 and 120: 776 J. R. Goicoechea et al.: The io
- Page 121 and 122: 778 J. R. Goicoechea et al.: The io
- Page 123 and 124: 780 J. R. Goicoechea et al.: The io
- Page 125 and 126: 782 J. R. Goicoechea et al.: The io
- Page 127 and 128: A&A 494, 977-985 (2009)DOI: 10.1051
- Page 129 and 130: M. Gerin et al.: HCO mapping of the
- Page 131 and 132: M. Gerin et al.: HCO mapping of the
- Page 133 and 134: M. Gerin et al.: HCO mapping of the
- Page 135 and 136: M. Gerin et al.: HCO mapping of the
- Page 137 and 138: Table 1. Observation parameters for
- Page 139 and 140: A&A 534, A49 (2011)tel-00726959, ve
- Page 141 and 142: A&A 534, A49 (2011)tel-00726959, ve
- Page 143 and 144: A&A 534, A49 (2011)tel-00726959, ve
- Page 145 and 146: tel-00726959, version 1 - 31 Aug 20
- Page 147 and 148: tel-00726959, version 1 - 31 Aug 20
- Page 149 and 150: tel-00726959, version 1 - 31 Aug 20
- Page 151 and 152: tel-00726959, version 1 - 31 Aug 20
- Page 153 and 154: tel-00726959, version 1 - 31 Aug 20
- Page 155 and 156: tel-00726959, version 1 - 31 Aug 20
- Page 157 and 158: ¡¢£¤¥¢¦¢ §¨ ©
- Page 159 and 160: tel-00726959, version 1 - 31 Aug 20
- Page 161 and 162: 7.2 MODE INTERFÉROMÉTRIQUE 161tel
- Page 163 and 164: 7.3 LE FUTUR DE LA RADIO-ASTRONOMIE
- Page 165 and 166:
CLASS evolution: I. Improved OTF su
- Page 167 and 168:
CLASS evolution: I. Improved OTF su
- Page 169 and 170:
A&A 526, A47 (2011)DOI: 10.1051/000
- Page 171 and 172:
S. Maret et al.: Weeds: a CLASS ext
- Page 173 and 174:
S. Maret et al.: Weeds: a CLASS ext
- Page 175 and 176:
IRAM-30m EMIR time/sensitivity esti
- Page 177 and 178:
IRAM-30m EMIR time/sensitivity esti
- Page 179 and 180:
IRAM-30m HERA time/sensitivity esti
- Page 181 and 182:
¬þ«ðZñÛùþ˜ï³þ«ðWMõÑ
- Page 183 and 184:
¨-A¨-, #f¡§¨\¤qŠ’‹¼€(
- Page 185 and 186:
¤H F!";
- Page 187 and 188:
#O.u\J0$!&%©ÉQ©6¦©¨-.u& (HÌ>
- Page 189 and 190:
¨*.G©6E¦cǦ©¨&ẗØn¨"!-‡
- Page 191 and 192:
¢¡¤£¦¥¨§ P© ¡ £¦¥¨§
- Page 193 and 194:
'¥0vB"(-7£O.'¥!5#&(-7)7)(:9.*d#J
- Page 195 and 196:
¥¤¡ b@‚¥_ƒO@02,.(-EFJbe O@#V
- Page 197 and 198:
02J$#J+'10-ï 3 02EFJ+')EKJ5L5†R(
- Page 199 and 200:
A&A 517, A12 (2010)tel-00726959, ve
- Page 201 and 202:
A&A 517, A12 (2010)tel-00726959, ve
- Page 203 and 204:
A&A 517, A12 (2010)tel-00726959, ve
- Page 205 and 206:
A&A 517, A12 (2010)tel-00726959, ve
- Page 207 and 208:
A&A 517, A12 (2010)tel-00726959, ve
- Page 209 and 210:
A&A 517, A12 (2010)tel-00726959, ve
- Page 211 and 212:
A&A 517, A12 (2010)tel-00726959, ve
- Page 213 and 214:
A&A 517, A12 (2010)tel-00726959, ve
- Page 215 and 216:
A&A 517, A12 (2010)tel-00726959, ve
- Page 217 and 218:
Table C.1. Definition of the symbol
- Page 219 and 220:
IRAM Memo 2011-2WIFISYN:The GILDAS
- Page 221 and 222:
WIFISYN3. practiceWIFISYN3. practic
- Page 223 and 224:
WIFISYNA. implementation planWIFISY
- Page 225 and 226:
tel-00726959, version 1 - 31 Aug 20
- Page 227 and 228:
8.3 PERSPECTIVES 227tel-00726959, v
- Page 229 and 230:
3 REQUIREMENTS 44 CHANGES FOR END-U
- Page 231 and 232:
5 CHANGES FOR PROGRAMMERS 125 CHANG
- Page 233 and 234:
A EXHAUSTIVE DESCRIPTION OF THE CHA
- Page 235 and 236:
Chapitre 9Copyright: Stéphane Guis
- Page 237 and 238:
9.3 ACTIVITÉS 2008-2011 237tel-007
- Page 239 and 240:
Contribution de l'Action Spécique
- Page 241 and 242:
A. des multi-pixels à bure : une s
- Page 243 and 244:
5. besoins en services annexes, bé
- Page 245 and 246:
Articles publiés dans des revues
- Page 247 and 248:
ARTICLES PUBLIÉS DANS DES REVUES
- Page 249 and 250:
Mémos IRAM et ALMAtel-00726959, ve
- Page 251 and 252:
Actes de colloques nationaux et int
- Page 253 and 254:
ACTES DE COLLOQUES NATIONAUX ET INT
- Page 255 and 256:
tel-00726959, version 1 - 31 Aug 20