Astronomy&Astrophysicsmanuscriptno.10905-xar (cid:13)c ESO2009 January26,2009 Imaging galactic diffuse clouds: CO emission, reddening and turbulent flow in the gas around ζ Oph H.S.Liszt1,J.Pety2,3,andK.Tachihara4 1 NationalRadioAstronomyObservatory,520EdgemontRoad,Charlottesville,VA,USA22903-2475 2 InstitutdeRadioastronomieMillime´trique,300RuedelaPiscine,F-38406SaintMartind’He`res,France 3 Obs.deParis,61av.del’Observatoire,75014,Paris,France 4 NationalAstronomicalObservatoryofJapan,2-21-1,Osawa,Mitaka,Tokyo181-8588,Japan receivedJanuary26,2009 9 ABSTRACT 0 0 Context.Mostdiffusecloudsareonlyknownaskinematicfeaturesinabsorptionspectra,butthosewithappreciableH contentmay 2 2 bevisibleintheemissionofsuchsmallmoleculesasCH,OH,andCO. Aims.Weinterpretingreaterdetailtheextensiveobservationsof12COemissionfromdiffusegasseenaroundthearchetypicallineof n sighttoζOph. a Methods.The12COemissionisimagedinpositionandposition-velocityspace,analyzedstatistically,andthencomparedwithmaps J 8 oftotalreddeningE∞B−VandwithmodelsoftheC+-COtransitioninH2-bearingdiffuseclouds. Results. Aroundζ Oph,12COemissionappearsintwodistinctintervalsofreddeningcenterednearE∞ ≈0.4and0.65mag,of B−V which<∼0.2magisbackgroundmaterial.Withineitherinterval,theintegrated12COintensityvariesupto6-12Kkms−1,comparedto ] A 1.5Kkms−1towardζOph.Nearly80%oftheindividualprofileshavevelocitydispersionsσv<0.6kms−1,whicharesubsonicatthe kinetictemperaturederivedfromH towardζOph,55K.Partlyasaresult,12COemissionexposestheinternal,turbulent,supersonic G 2 (1-3kms−1)gasflowswithespecialclarityinthecoresofstronglines.Theflowsaremanifestedasresolvedvelocitygradientsin . narrow,subsonically-broadenedlinecores. h p Conclusions.ThescatterbetweenN(CO)andEB−Vinglobal,COabsorptionlinesurveystowardbrightstarsispresentinthegasseen aroundζOph,reflectingtheextremesensitivityofN(12CO)toambientconditions.Thetwo-componentnatureoftheopticalabsorption - o towardζ Ophiscoincidentalandthestarisoccultedbyasinglebodyofgaswithacomplexinternalstructure,notbytwodistinct r clouds.Theverybright12COlinesindiffusegasariseatN(H2)≈ 1021 cm−2 inregionsofmodestdensityn(H)≈ 200−500cm−3 t andsomewhatmorecompleteC+-COconversion.Giventhevarietyofstructureintheforegroundgas,itisapparentthatonlylarge s a surveysofabsorptionsightlinescanhopetocapturetheintrinsicbehaviorofdiffusegas. [ Keywords. interstellarmedium–molecules 1 v 1 1. Introduction L204seenseveraldegreestothegalacticSouth.L204isclearly 1 outlined against the Hα emission from the ionized gas in the 1 Thelineofsighttothenearby(140-160pc)runawayO9.5Vstar star’sHIIregion(Gaustadetal.,2001).Inthiswork,the12CO 1 HD149757,ζ Oph,hasservedasthearchetypefordetailedob- datacube from Tachihara et al. (2000) is employed to study the 1. servationalstudiesoftheinternalcompositionofdiffuse(AV<∼1 diffusegasatA ≈1magseennearerthestar.Wescrutinizethe 0 mag) clouds (Herbig, 1968; Morton, 1975), for optical/uv de- V entireCOimageoftheabsorption-linehostwhoseoverallprop- 9 tectionofnewmoleculesindiffusegas(Maieretal.,2001)and ertieshavesooftenbeeninferredfromonemicroscopicabsorp- 0 fortheoreticalmodelsofmoleculargasindiffuseclouds(Black tionsightlinetowardthestar,andweinquiretowhatextentthat : v & Dalgarno, 1977; Van Dishoeck & Black, 1986, 1988; Kopp lineofsightfaithfullyrepresentsthehostgas.Moreover,large- i etal.,1996).TheH2-bearingportionsofthegasoccultingζOph scalemapsofreddeningandextinctionhavebecomeavailableat X aredenseenoughtohostappreciablecolumndensitiesofcarbon comparableresolution(thoughonlyalongtheentirelineofsight, r monoxide,N(12CO)≈2.4×1015 cm−2(Morton,1975;Wannier seeSchlegeletal.(1998)andDobashietal.(2005)),andweem- a et al., 1982; Lambert et al., 1994; Sonnentrucker et al., 2007), ploythesetocontrolagainstpossibleconfusionbetweendiffuse and these are readily detectable in mm-wave emission toward and darker sightlines, a source of concern given the strong CO thestar(Knapp&Jura,1976;Liszt,1979;Langeretal.,1987). lineswesee. They were very partially mapped in CO emission (Kopp et al., The plan of this work is as follows. Section 2 summarizes 1996;Liszt,1997),aswellasCHandOH(Crutcher,1979;Liszt, whatisknownobservationallyofthelineofsighttowardthestar 1997). and describes the pre-existing H I, CO, A∞ and E∞ datasets COJ=1-0emissionaroundζOphwasimagedinmuchmore whicharediscussedhere.Section3discussVestheapBp−eVaranceof complete fashion by Tachihara et al. (2000), who focused their theskyaroundζ Ophintermsofthestatisticsof12COemission discussion on the properties of the nearby dark cloud complex andreddening.Section4discussesCOprofilesandlinewidthsin termsoftheturbulentflowswhichareprominentintheemission Sendoffprintrequeststo:H.S.Liszt profiles. Section 5 discusses physical conditions in the CO and Correspondenceto:[email protected] H -bearing host gas, especially the regions of extremely bright 2 2 Liszt,PetyandTachihara:CO,reddeningandturbulencearoundζOph Fig.1. Limiting reddening (Schlegel et al., 1998) and 12CO J=1-0 emission (Tachihara et al., 2000) in the vicinity of ζ Oph. The 12COemissioncontoursareshownatlevels1,2,4,6,8,10,12,16,20,24Kkms−1andhavebeencalculatedoverdifferentvelocity rangesintheupperrightandlowerleftportionsofthemap,asindicatedatthemapcorners.ThedarkcloudcomplextotheSouth ofζ OphisknownasL204andthegasattheupperright,aroundζ Oph,asL121. (11-12K)COemission.Sect.6isasummaryandSect.7(avail- 2. Observations able online) discusses the relationship between the extinction and reddening measurements over the ζ Oph field and presents 2.1. Carbonmonoxide someadditionalviewsofthe12COobservations. ThedatacubeofTachiharaetal.(2000)comprisesnearly11,000 spectrafromtheNANTENtelescopewithabeamwidthHPBW =2.7(cid:48) ona4(cid:48) gridingalacticcoordinates.Thespectrahave0.1 kms−1 resolution and the single-channel rms at this resolution, 0.5 K, is relatively high compared to that in the small numbers Liszt,PetyandTachihara:CO,reddeningandturbulencearoundζOph 3 islikelyhostedinatomicgasoverlongpaths,asopposedtothe more-localizeddiffuseanddarkcloudsofinteresthere. A comparison between the limiting reddening of Schlegel et al. (1998) and maps of the limiting extinction A∞ computed V fromstarcountsbyDobashietal.(2005)canbefoundinSect. 7(onlineonly). 2.4. Somegeneralconditionsinthegasalongthelineof sighttoζ Oph Infrontofthestar,E =0.32mag,N(HI)=5.2×1020cm−2, B−V N(H )=4.5×1020cm−2(Savageetal.,1977),N(12CO)≈2.4× 2 1015 cm−2 (Morton, 1975; Wannier et al., 1982; Sonnentrucker etal.,2007),N(C0)=N(12CO)(Morton,1975)andN(C+)>∼3× 1017 cm−2 (Cardellietal.,1993).Themeankinetictemperature ofthemoleculargasinferredfromH absorptionis54K(Savage 2 etal.,1977). Toward the star the limiting extinction from the work of Schlegeletal.(1998)isE∞ =0.55mag,sothatthebackground B−V reddeningisapproximately0.23mag. The distance to the occulting material is generally taken to be very close to that of the star, just outside the nearer edge of Fig.2. Integrated H I line brightness from the LDSS survey the star’s H II region (Wood et al., 2005). At a distance of 140 of Hartmann & Burton (1997) at 35(cid:48) resolution, converted to pc,1(cid:48)correspondsto0.041pcand1oto2.44pc. equivalent reddening E = N(H I)/5.8 × 1021 cm−2, N(H I) (cid:82) B−V =1.823×1018cm−2 T (HIdv)/(Kkms−1). B 3. Theskyaroundζ Ophviewedinreddeningand COemission ofdemonstrationspectratypicallyshowninearlierwork(Liszt, Figure 1 is a composite image of the limiting reddening E∞ 1997).StatisicallysignificantdetectionsoftheCOrequireW B−V CO andintegrated12COJ=1-0emission,anupdatedversionofFig. >∼1kms−1,whereWCOistheintegratedintensity 1 in Tachihara et al. (2000). The gray-scale underlay is the Toward the star, we show the profile of Liszt (1997) from limiting reddening (Schlegel et al., 1998) normalized to white the then-NRAO Kitt Peak 12m telescope at 1(cid:48) (HPBW) spatial at the minimum value seen over the region, E∞ = 0.23 mag. resolutionand0.12kms−1 spectralresolution(seeSect.3).To B−V Superposedonthereddeningmap,theintegrated12COintensity ensurethatthisprofileiscompatiblewiththosefromNANTEN W hasbeencalculatedseparatelyoverthediffuse/translucent werecentlyusedtheARO12mKittPeaktelescopetore-observe CO northwest region, referred to as L121, and the translucent/dark several positions having comparatively strong emission in the southeast region, L204. For the diffuse gas of L121 at upper NANTENdatacube.The12mspectraagreewiththeNANTEN right,theredandbluecontourscorrespondtothevelocityranges datatobetterthan5%,aremarkablecoincidenceconsideringthe differenceinhardwareandspatialresolution. aboveandbelowv=-0.25kms−1asindicatedintheupperright corner: this division corresponds to the natural separation be- tweenthetwocomponentsofCO,CHandOHemissionfound 2.2. HI towardandaroundζ Oph,asshownintheprofilesandposition- velocitymapsofLiszt(1997).ForthedarkergasofL204tothe To help distinguish between foreground and background mate- southeastW wascalculatedoverthreeintervalsbutnearlyall CO rial, or atomic and molecular gases, we employed the H I pro- of the emission from L204 occurs in the velocity interval 2-6 files from the Leiden-Dwingeloo all-sky H I survey (Hartmann kms−1representedbythewhitecontours(seeFig.5).Thenoise &Burton,1997).Thesedatahave35(cid:48) resolutionona0.5o grid level is somewhat larger for these contours, because W was CO ingalacticcoordinates. calculatedoverasomewhatbroaderinterval. ThestrongCOemissionassociatedwithL204oftenfollows theridgelinesoftheextinctionwithsomethingofasetback(red- 2.3. Reddeningandextinction dening without apparent CO) in the direction of the star. This The reddening toward ζ Oph is known to be E = 0.32 mag is consistent with the edge-on geometry for L204 described by B−V (Morton,1975),butmapsoftheforegroundextinctionoveronly Tachihara et al. (2000). Nearly all of the CO emission in the thefirst140pcareunavailable.Instead,weemploythelimiting L204 dark cloud complex is found in the interval v = 2-6 lsr reddening(fromheretoinfinity)fromtheworkofSchlegeletal. kms−1 and so does not overlap that of the more diffuse gas to (1998),denotedbyE∞B−V,withaspatialresolutionof6.1(cid:48),pub- the north, seen at vlsr <∼ 2 kms−1. The extent to which the two lished on a 2.5(cid:48) grid; the stated global rms error of this dataset regionsareatrestinthedirectionsjoiningthem,andmightpar- is 16% (a percentage at each pixel). Their values are based on takeoftheopticalpumpingexcitationmechanismfordiffusegas a determination of the dust column density estimated from the describedbyWannieretal.(1997)isunknown.Tachiharaetal. IRAS 100 micron flux adopting the temperature variation de- (2000)assumedthatthetwocloudswereco-movingattheedge rivedfromCOBE/DIRBE240microndata.Theminimumlim- oftheHIIregionaroundζOphinordertodiscusstheenergetics itingreddeningintheregion,approximately0.23mag(Fig.1), ofthegas. 4 Liszt,PetyandTachihara:CO,reddeningandturbulencearoundζOph which include all of the unrelated background and foreground gas.Immediatelyaroundζ OphE∞ =0.55mag,comparedto B−V E = 0.32 mag in front of the star, so some 0.55-0.32 mag B−V = 0.23 mag of reddening occurs behind the star. This also cor- responds well to the absolute minimum in Fig. 1, E∞ = 0.23 B−V magfoundsome1.5oNorthofthestar. To estimate the amount of extraneous foreground material, wenotethatthisisexpectedtobeinatomicformand,inabsorp- tionagainstthestar,N(HI)= 5.25×1020 (Savageetal.,1977) correspondingtoEB−V<∼0.1mag.ThisleavesEB−V>∼0.22mag of foreground material associated with the pure H component 2 and gives the general idea that perhaps as much as 0.32 mag should be subtracted from the map of E∞ to infer the local B−V reddeningintrinsictotheH -bearinggas.However,thisisprob- 2 ablyanoverestimateoftherequiredcorrectionbecausesomeof the foreground atomic gas is in the vicinity of L121 and L204, providingshielding. Figure 2 is a map of the integrated H I brightness, scaled to the mean relationship between hydrogen column density and extinction, that is, the map shows 1.823 × (cid:82) 1018 cm−2 T (HI)dv/5.8×1021 cm−2mag−1 where the veloc- B ityintegralisinunitsofKkms−1.Severalofthelocalreddening peakswhichlackCOemissioninFig.1evenatE∞ =0.65mag B−V arepresentaspeaksinFig.2,forinstanceat(l,b)=(7.5o,21.5o) Fig.3.IntegratedintensityoftheCOlineW vs.limitingred- deningE∞ .Top:forthediffuseregionL12C1Oatupperrightin and(4o,20o).Someofthisgasmustbeindigenoustotheregion B−V of interest, but the extinction associated with the H I peaks is Fig. 1; bottom: for the darker L204 gas at lower left in Fig. 1. only ≈ 0.06 mag, judging from the peak levels in Fig. 2 (0.28 Ineachcasethedataareshowntwice,atrightasindividualdata mag)comparedtothenearbybackgroundlevel(0.22mag). points, at left binned and contoured at logarithmic (factor two) TheminimainFig.1andFig.2andthedifferencebetween intervals: the contour representation of the L121 data is super- the foreground and limiting reddening toward the star consis- posed at bottom. Also, in each case the data are divided into tentlyimplyabackgroundreddeningcontributionof≈ 0.2mag regions of lower and higher extinction with the separation oc- curingatthemeanE∞ overdatapointswhichlackstatistically overtheregionofinterest. B−V significantdetectionsofCOemission(0.46magand0.67magat topandbottom,respectively).Thelocusofthesightlinetoward 3.2. QuantitativerelationshipbetweenCOemissionand ζ OphintheB-regionoftheupperpanelisindicatedbya(red) reddening cross.Mean12COemissionprofilesareshowninFig.5. To quantify the relationship between reddening and CO emis- sion,andtocompareandcontrastthediffuseanddarksightlines TotheNorthofL204,oneithersideofthestaristheζ Oph we divided the extent of Fig. 1 along a Northeast-Southwest diffusecloud,describedastheL121complexbyTachiharaetal. diagonal in the trough of reddening between L121 and L204, (2000).Mostoftheemissionfromthatgasoccursat-2.4kms−1 alongalinerunningfrom(l,b)=(3o,19.333o)to(10o,25o).Fig. <vlsr<2.2kms−1.Thered-shiftedkinematiccomponentat-0.2 3showstheintegratedCOintensitiesWCO andlimitingredden- kms−1 <v <2.2kms−1whichisseeninCOemissionandin ing E∞ for all points in both regions: the profile integral was lsr B−V many species in absorption toward the star, is very much con- takenovertherange-2.4kms−1 <v <2.2kms−1 forthedif- lsr finedtothenorthernandnorthwesternedgesofthebroaderdis- fuse L121 gas shown in the top panel, and 2 kms−1< v < 8 lsr tribution of blue-shifted gas. Its separate identity as a second kms−1 for the darker L204 region in the lower portion of the cloud is somewhat marginal in the emission maps and hardly map of Fig. 1 and the lower panel in Fig. 3. The rms noise in supported by a more detailed examination of the turbulent gas integratedintensityis0.5-0.6Kkms−1.Profileintegralsabove1 kinematics(alsoseeFig.7). Kkms−1generallyrepresentrealdetections. The diffuse gas in the L121 complex around ζ Oph is sep- OncewerealizedthattheCOemissioninthediffuseregion arated from the L204 dark cloud by an extended trough in wasbimodal,asillustratedinthetoppanelsofFig.3,wefurther the reddening whose overall mean value is E∞ ≈ 0.43 mag. sub-dividedthediffusegasintoAandBportionscorresponding B−V Immediately below the star is a pronounced minimum whose to the two branches of the emission distribution. The A and B meanisE∞ ≈0.34magandwhoseabsoluteminimumisE∞ portionswereseparatedatE∞ =0.455mag,whichisthemean =0.29maBg−.VTheoverallimpressionisofacylindricalshellBg−eV- E∞ for those sightlines alBo−nVg which W < 1 K kms−1and ometry and perhaps a separate, more circular shell of radius ≈ wBh−icVh therefore lack statistically significaCnOt detections of CO 1o around the star. However, a yet-larger map of the extinction emission.TheA-branchpixelshavestrongCOemissionatE∞ showsthatL204ispartofamuchlargerridge. substantially below the mean of those sightlines lacking BC−OV emission at all. For the diffuse gas in L121 there is actually a substantial spatial segregation of the A and B portion pixels, 3.1. Separatingforegroundandbackgroundreddening with unweighted mean < (l,b) > = (6.0o±0.7o,23.8o±0.7o) and Asnotedabove,wewishtounderstandtherelationshipsbetween <(l,b)>=(4.0o±0.9o,22.3o±1.0o)fortheAandBportions,re- W and reddening, but we only have maps of the reddening spectively. These centroids are on opposite sides of ζ Oph and CO Liszt,PetyandTachihara:CO,reddeningandturbulencearoundζOph 5 3.3. Statisticsofthelinebrightness Figure4presentsthedistributionofintegratedandpeakbright- ness. The differences between the A (lower E∞ ) and B por- B−V tions of the diffuse gas are somewhat more pronounced in the distributionofthelineprofileintegralatleftandsomewhatless so in the peak, so that the line widths differ more than the line heights.Althoughthemedianbrightnessisatmostonly3Kor3 Kkms−1inthediffusegasandthetailofthedistributionseems very poorly populated above, say, 8 K peak brightness: paths which traverse L121 even in the short dimension have a very substantialchanceofcontainingatleastonesuchbrightline,as discussedinSect.5(seeFig.7). It is not the intention here to discuss the dark gas, but it shouldbenotedthatthedifferenceinmeanbrightnessbetween thediffuseanddarkregionsaremodestandcorrespondapprox- imately to the differences in E∞ , thereby preserving the pos- B−V sibilityofacommonCO-H conversionfactor;thesameisalso 2 true for the A and B portions of L121, see Fig. 5. This occurs despitethefactthatmostofthefreegas-phasecarbonisinC+in L121 (99% toward ζ Oph) and in CO in the dark gas of L204 (whereN(CO)≈3×1017cm−2)implyingadifferenceinCOcol- umndensityandW /N(CO)byafactoroforder50.InSection CO 5 we discuss the very different proportionality W ∝ N(CO) CO Fig.4.Cumulativeprobabilitydistributionsoftheintegratedin- whichisobservedwithinthediffuseregimealone(Liszt,2007a). tensityW (left)andpeaklinetemperature(right)overthedif- CO fuse,northeasternL121(attop)anddarksoutheastL204sight- 3.4. AnincidentalboundonN(H )overthediffusegas linesshowninFig.1.OnlythosesightlinesatwhichW ≥1K 2 CO kms−1havebeenconsidered. Becausethelineofsighttoζ Ophoccursatsuchahighvalueof E∞ relativetotherestofL121weinferthatN(H )isneververy B−V 2 muchlargerinL121thantowardthestar.Forinstance,ifwetake E∞ = 0.65 mag characteristic of the B-portion and subtract a separated by more than 1σ+1σ in each coordinate. The red- baBc−kVgroundcontribution0.23magequaltothattowardthestar, shifted gas appears mostly in the higher-extinction B-portion theremaininggascolumnwithE∞ =0.65-0.23=0.42mag B−V while the blue-shifted gas appears more nearly in both the A corresponds to N(H) ≈ 2.4 × 1021 H-nuclei cm−2 and N(H ) 2 andB-portions. < 1.2×1021 cm−2.Conversely,becausethelineofsighttoward In either emission branch, WCO varies widely over a rela- the star has such a high value E∞ = 0.55 mag relative to the tivelynarrowrangeinE∞ ,mimicingtheextremesensitivities rest of the region, emission fromB−Vthe A-branch at E∞ = 0.4 B−V B−V ofN(CO)whichareseeninglobalabsorptionlinesurveysofsin- mag,probablyarisesinregionswhosereddeningandN(H )are 2 glesightlinesoverwidely-separatedregions(seeLiszt(2007a)). actuallybelowthoseseentowardthestar. Thewidthsofthetwobranchesarecomparabletothestatedrms noiseoftheE∞ data(16%)andtheyareseparatedby0.2mag, B−V which is of the same order as the background contribution: it 4. Lineprofiles,linewidthsandturbulentflowsin ispossiblethattheywouldmorenearlycoincideinEB−V ifthe diffusegas background contribution were strongly variable. However, the B-branchathigherE∞ ismoreheavilypopulatedatW >5 4.1. COlineprofilesandprofilewidths B−V CO Kkms−1,suggestingthatitactuallyissomewhatmorestrongly For a pure-H gas, the Doppler temperature equivalent to a 2 shielded,fosteringahigherCOabundanceandbrightness:inthis given linewidth is T = 43.6 K FWHM2 or T = 242 K dopp dopp regime,WCO∝N(CO)(ibidandseeSect.6)andtheconversion σ 2 where σ is the velocity disperson in kms−1. Typical ki- offreecarbonfromC+toCOoccursoveraverynarrowinterval nevtic temperavtures in CO-bearing diffuse gas are above 30 K inE and/orN(H ). B−V 2 (Sonnentrucker et al., 2007; Burgh et al., 2007; Liszt, 2007a) The distribution of brightness in the darker L204 region is andthemeankinetictemperatureseeninsurveysofH absorp- 2 not similarly bimodal at all E∞ and there is less rationale for tionis70-80K(Rachfordetal.,2002;Savageetal.,1977).CO B−V asimpledivisionintosub-portions.Howeverwhenthisisdone, profiles with FWHM <∼ 1 kms−1 are subsonic at diffuse cloud atE∞ =0.67magcorrespondingagaintothemeanoverpixels temperatures and the purely-thermal velocity dispersions of B−V lackingstatisticallysignificantCOemission,theresultantlower- COmolecules,0.1kms−1atT =50K,donotcontributeim- K E∞B−V portioncorresponds(inE∞B−V andWCO)totheentiretyof portantlytotheobservedlinewidths. thediffuseregion.Consequentlythelower-E∞ partofL204is Unweighted mean profiles formed over the sub-portions of B−V labelled L204 AB and the other L204 C. The behavior of WCO thediffuseL121anddarkL204regionsareshowninFig.5.The with E∞ in the dark gas is complex but clearly bimodal for meanprofileshavelinewidthswhicharesupersonic,FWHMof B−V WCO >∼12Kkms−1;thestrongestemissionisbynomeanslim- typically 2-3 kms−1 but the individual sightlines in L121 typi- ited to the darkest regions. Statistics of the brightness distribu- callyhavesubsoniclinewidths;theFWHMofthetwokinematic tion over L121 and L204 are shown in Fig. 4 and discussed in componentsseentowardthestarat1(cid:48) resolution(showninFig. thefollowingsections. 5asthedarkdashedline)are0.6kms−1and1.1kms−1,equiva- 6 Liszt,PetyandTachihara:CO,reddeningandturbulencearoundζOph Fig.5. Mean spectra over the diffuse and dark sub-portions de- Fig.6. Statistics of line profile velocity dispersions in the dif- finedinFig.3:theprofilesinthediffusegas(L121)arethoseat fuse gas. Large outer panel: cumulative probability distribution lower velocity. The integrated brightnesses are W = 2.9 and histograms of the measured velocity dispersion σ shown sep- CO v 4.0Kkms−1 forL121and2.3and7.4Kkms−1 forL204.Also arately for the A and B portions of the diffuse gas (see Sect. 4 shown(dashedblackline)isaspectrumat1(cid:48) resolutiontoward of the text). The one-dimensional thermal velocity dispersions ζ Ophfromthe12mtelescope. correspondingtokinetictemperaturesof20,40and80Kinagas ofpureH areshownatupperleft.Insetatlowerright:variation 2 ofσ withW . v CO lenttoDopplertemperaturesofonly16Kand53KinapureH 2 gas,comparedtoameasuredtemperatureinH of54Kasnoted 2 inSect.2.4. gradientsinthehostmedium,i.e.thecharacteroftheturbulence. ThedistributionofvelocitydispersionsfoundovertheL121 Althoughveryelaborateanalysesoflineprofilecentroids(Pety region is shown in Fig. 6. To produce this figure we used the & Falgarone, 2003) and fluctuations in the wings of optically followingwindowingtechniquetomeasurethewidthsofspec- thickprofiles(Falgarone&Phillips,1996)havebeenusedtoin- tra which might contain more than one kinematic component: ferthepropertiesofturbulenceindenserneutralmedia,turbulent select the overall velocity interval of the diffuse gas; find the flowsinthenearbydiffusegasofL121areimmediatelyvisible peakchannel;selectthosecontiguouschannelsaroundthispeak inthelinecoresand(usually)spatiallyresolvedintotheshiftsof with temperatures above a noise threshold of 0.25 K; calculate individual,subsonicprofiles. thebrightness-weightedvelocitydispersionoverthosechannels; ThisisillustratedforζOphinFig.7,whereweshowacom- maskoffthatportionoftheprofile;repeatuntilnochannelabove prehensiveseriesoflongitude-velocitydiagramsspacedatone- 1 K remains unmasked. The dispersions measured en masse in pixel(4(cid:48))intervalsingalacticlatitude,aslabelledintheindivid- thiswayaresubjecttooverestimatingthewidthincasesofun- ualpanels.Thepeaklinebrightnesses(K)seenintheindividual recognizdblending,buttheyagreetowithin10%foraseriesof panelsofFig.7arelabelledonthebarsshowingthecolorscal- test profiles which were fit with gaussian components (for in- ingineachpanel.Althoughlineswithpeakbrightnessesabove stance,thoseshowninFig.8). 8KseemrelativelyrareinFig.4,theyarecommonenoughthat The general properties of the gas in L121 have often been peak brightnesses 7.5 K and higher appear in 40%, 18 of the inferredfromprofilesseenalongthesingle,microscopicabsorp- 44individualpanelsofFig.7.AsshowninFig.3,themostin- tion line of sight toward the star and a sensitive CO profile to- tenselinesinthediffusegasarebynomeanslimitedtothemore ward ζ Oph at 1(cid:48) spatial resolution is also included in Fig. 5. heavily-reddenedsightlines. How representative is it of the larger-scale distributions of the Figure7showsthatadescriptionofthegasintermsoftwo host gas? The stronger CO component at vlsr = -0.7 kms−1 to- identifiableforegroundcloudsattheabsorptionlineorCOemis- ward ζ Oph is one of the narrowest lines known in a diffuse sionlinevelocitiestowardζ Ophisnotappropriate.Thelineof cloud, with FWHM = 0.60 kms−1. Broad consideration of this sight to the star could just as well have occured behind any of questionisgiveninthediscussionofFig.7whichdisplaysmost theprofilesexhibitedinFig.7,leadingtoawidevarietyofpos- oftheprofilesseenintheL121region,immediatelyfollowing. sible interpretations of the intervening medium. Discussing the mapinFig.1,wesawthatthered-shiftedcomponentgenerally appears as something of a fringe at the northern edges of more 4.2. Velocitygradientsandturbulentflows broadly-distributed, negative-velocity gas. Viewed in position- The connections between the profiles seen at individual pixels velocityspaceinFig.7,thered-shiftedgasisoftenseenasapro- andthemeanprofilesshowninFig.5aretheflowsandvelocity nouncedkinematicexcursionorwing,forinstance,atb=22.8o Liszt,PetyandTachihara:CO,reddeningandturbulencearoundζOph 7 or 23.4o. Further to the North, at b > 23.2o, the kinematic pat- includinggrain-assistedneutralizationofatomicions(including ternundulatesacrossthepositive-velocityportion.Atb=23.4o mostimportantlytheprotons,ibid),toallowequilibriumofthe theresolvedvelocitygradientspanstheentirerangeofvelocities self-shieldingof H formedon grainsandthe formationof CO 2 inthediffusegas:viewedfromadifferentperspective,thissame bythermalelectron-recombinationofafixedquantityofHCO+ regionmighthavebeenseenasasingle,broadline. X(HCO+)=N(HCO+)/N(H )=2×10−9(theactualsecularevo- 2 Figure 8 is an expanded view of the panel in Fig. 7 at b= lutionistracedbyLiszt(2007b));and,tocalculatetherotational 22.6o, where we also show several included line profiles, their excitationofCOassumingmicroturbulentradiativetransferwith gaussiandecompositionandtheresultantFWHM.Exceptatthe a linewidth determined by the local sound speed. Note that the centerofthediagramwherethereisapartiallyspatially-resolved interstellar radiation field in the models has not been increased velocity gradient, the FWHM are small enough to be subsonic abovethemeaninterstellarvaluetoaccountforthepresenceof for an H gas at typical diffuse cloud temperatures above 30 thestarandthatalloftheaspectsshownhereforthemodelre- 2 K.Thustheobservationsshowindetailhowsupersonicprofiles sultsarepresentintheabsorptionlinedatasummarizedbyLiszt mightarisefromthecoincidentalsuperpositionoraddition(for (2007a). instance through beam-smearing) of velocity-shifted quiescent Model results are shown for just two fairly high densities regions. n(H)=256cm−3and512cm−3.ToformthegraphsinFig.9,re- AsinFig.8therearemanyhighly-localized,relativelybroad sultswerederivedbyintegratingalongsightlinesrangingacross linesin4(cid:48)pixels,oftenjoinedtoabruptbutspatially-resolvedve- thefacesofmodelspheresfromcentertoedge.Thegraphssum- locitygradientsandreversalsinnarrow-linedgas.Thissuggests marizeresultsgleanedfrommodelswhosecentralcolumnden- that broader lines are themselves composed of unresolved ve- sitiesN(H)variedwidely,sothatthesamevalueofN(H )might 2 locity gradients and it seems possible that any profile in Fig. 7 occuratdifferentpositionsacrossthefacesofmodelswithdif- withawidthsubstantiallyabovesonicisanunresolvedgradient. fering n(H) and N(H), and therefore have different N(CO) and For instance, compare the velocity span at b=22.87o, which is W .ForN(H )=1×1021 cm−2 themodelshavetypicalsizes CO 2 spatiallyresolved,withthatattheedgeoftheemissionregionat of1.2-2.5pcbuttheCOabundanceisconcentratedintosmaller b=23.33o.Recallingtheextremelynarrowblue-shiftedCOcom- centralportionsofthehostbodiesowingtothechemistryofCO ponenttowardζ OphandthelowimpliedDopplertemperature, and that of H , and the CO emission is more concentrated still 2 wewonderwhetherandatwhatscaleprofileshavinglinewidths owingtogeometryandradiativepumping. thataresubsonicatthehighertemperaturesexpectedfordiffuse As noted above N(H ) = 0.5 − 1 × 1021 cm−2 in the CO- 2 gas,say60K,mightshowvelocityorspatialsubstructure.The emittingregionsaroundthestar.Asshowninpanelcatthelower geometryoftheturbulentflowsproducingthepatternsinFig.7 left,thisispreciselytheregimewherecarbonisabouttorecom- willbeconsideredinasubsequentpaper. binefullytoCOatthequoteddensities:theincreaseofN(CO) with N(H ) is very rapid. Both N(CO) and W vary rapidly 2 CO and have large scatter when plotted against N(H ). Substantial 5. Physicalconditionsinthediffusegas 2 CO column densities can accumulate in gas which is still rela- There are some 1000 sightlines with WCO >∼ 1 K kms−1 in tively warm, 30-50 K, giving rise to 12 K lines as observed in the diffuse, Northwest portion of the map in Fig. 1. At a dis- thebrightestprofiles. tance of 140 pc, the H -mass associated with these sightlines, The most important consideration is as shown at the upper 2 parametrized by their average H2 column density <N(H2)>, is rightinpanelb,WCO ∝N(CO)forWCO <∼10Kkms−1,explic- M=430 M <N(H )/1×1021 cm−2 >wheretypicalvaluesof itlyindependentofdensityandimplicitlyindependentofN(H ) 0 2 2 N(H )alongtheL121sightlinesare0.5−1.2×1021 cm−2 (see and other cloud properties. This is a very general consequence 2 Sect.2.4and3.4).Thismaybecomparedwiththevalue520M of very sub-thermal excitation, as first shown by Goldreich & 0 given by Tachihara et al. (2000) based on a CO-H conversion Kwan(1974)anddoesnotrequirelowopticaldepth.Asshown 2 factorN(H )/W =1.56×1020cm−2/(Kkms−1). in Fig. 4, some 80-90% of the diffuse cloud spectra have W 2 CO CO Nonetheless, deriving the physical parameters of host dif- < 5-6 K kms−1 and virtually all are below 10-12 K kms−1, fuse gas from CO profiles is challenging. Strong fractionation just in the regime characteristic of sightlines studied in uv and of the carbon isotopes (Liszt & Lucas, 1998; Sonnentrucker mm-wave CO absorption work generally, see Fig. 6 of Liszt et al., 2007; Burgh et al., 2007; Liszt, 2007a) causes the (2007a).Therefore,the12CObrightnessmapinFig.1isamap N(12CO)/N(13CO)tovaryintherange15-150,thusmakingit of N(12CO) in the diffuse gas and the same would be true for impossibletoderivethelineexcitationtemperaturesandoptical N(13CO) and the brightness of the 13CO line. This is the one depths,orthekinetictemperaturesandcolumndensities,under unambiguousresultofmappingCOemissioninanydiffusegas. the usual assumption (valid in dark gas) that the relative abun- TheextremesensitivityofW toN(H )thereforearisesbe- CO 2 dances of 12CO and 13CO only reflect the general interstellar causeW ∝N(CO),sothataplotofW vs.N(H )isequiv- CO CO 2 carbonisotoperatio.However,itisstraightforwardtoshowthat alenttoplottingN(CO)againstN(H ).Thenetresultisthatal- 2 thegeneralpropertiesoftheCOobservations,evenincludingthe though commonly-used values of the CO-H conversion factor 2 ratherunexpectedlybrightlines,fiteasilyintotheframeworkof apply to some gas (as indicated in panel a of Fig. 9) the CO- diffusegasattypicaltemperatures30-60Kandmodestdensity. H conversionfactorvarieswidelyindiffusegasandtheactual 2 The real underlying mysteries are the working of the poorly- N(CO)/N(H ) ratio is small but very uncertain. A map of CO 2 understoodpolyatomicchemistrywhichformstheCOandother emissionfromdiffusegasisanimageofthechemistry,notthe speciesatsuchmodestdensities(Lisztetal.,2006),andtheori- massdistribution. ginoftheturbulentflowswhichmaypowerthechemistry. Last,notethattheJ=1-0COlinebrightnessisinsensitiveto Figure 9 reports some results from the models which were densityatfixedN(CO),indicatingthatothertracersarerequired usedtointerpretCOabsorptionlinedatabyLiszt(2007a).The to measure the local density when mm-wave emission profiles underlying physics are: to heat a uniform-density gas sphere are analyzed. The J=2-1/J=1-0 line brightness ratios at lower by the photoelectric effect on small grains, as in the work of right in Fig. 9 are not very sensitve to density, which explains Wolfire et al. (1995, 2003), to calculate the ionization balance why line brightness ratios 0.7 - 0.75 are indeed so commonly 8 Liszt,PetyandTachihara:CO,reddeningandturbulencearoundζOph Fig.7.Longitude-velocitydiagramstraversingtheζ OphdiffusecloudatmanylatitudeswithintheNorthwestportionoftheregion showninFig.1.Theverticalseparationbetweenpanelsis4(cid:48) (1pixel)andthegalacticlatitudewithineachpanelislabelled.The colorscaleineachpanelareseparatelynormalizedwithinarange(Kelvins)shownontheindividualcolorbars. observed in diffuse and translucent gas (Falgarone et al., 1998; 6. Summary Pety et al., 2008). The J=3-2/J=1-0 brightness ratio is a better Thelineofsighttothenearby(140-160pc)runaway09.5Vstar indicator of density in the CO lines, but care must be taken to ζ Ophhasformanyyearsbeenusedasanarchetypeforstudy- matchthespatialresolutionofthetwolines. ingthepropertiesofdiffusecloudsinopticalanduvabsorption. Because the material has an appreciable molecular content, the host diffuse clouds can actually be imaged on the sky in space and radial velocity. Because the gas is well extended and com- parativelyclose,itprovidesanunusualopportunityforstudyof Liszt,PetyandTachihara:CO,reddeningandturbulencearoundζOph 9 Fig.9. Models of 12CO formation and excitation. Results are shown for models having n(H) = n(H I) + 2n(H ) = 256 cm−3 2 and 512 cm−3. Upper left: Variation of integrated 12CO bright- Fig.8.Thelongitude-velocitydiagramatb=22.6o fromFig.7, ness WCO with N(H2); the shaded line labelled “2” is a CO- H conversion factor 2 × 1020 H cm−2 (K kms−1)−1. Upper expanded and decomposed into its constituent profiles at three 2 2 right: variation of W with N(12CO). Bottom left: variation positions (shown inset). Overlaid on each inset spectrum is a CO of N(12CO) with N(H ). Lower right: comparison of 12CO line one-componentgaussianfitwhoseFWHMisshownattheupper 2 brightnessesinthelowesttworotationaltransitions. left. diffuse gas and its interaction with its surroundings, including flows in this gas are in general directly visible as the spatially thestar. andkinematicallyresolvedvelocitygradientsinsimple,narrow, We began by comparing the 12CO J=1-0 emission line dat- brightlineswhosewidths(σ <0.6kms−1)aresubsonicatdif- v acubeofTachiharaetal.(2000)(HPBW2.7(cid:48)observationsona4(cid:48) fuse cloud temperatures T = 30 - 60 K (Fig 4 and 8). Other, K =0.1pcpixelgrid)withamapofthelimitingreddeningE∞ locally-broaderlineprofileswilllikelyberesolvedintosuchve- B−V from the work of Schlegel et al. (1998) having similar resolu- locityshiftsofnarrowlinecoreswithhigher(than4(cid:48))resolution tion6.1(cid:48) (Fig.1).ThereddeningintheL121complexnearand although this remains to be tested. Conversely, it also remains aroundζ OphrangesfromE∞ =0.23-0.75magandtheinte- tobeseenwhetherCOlineprofileswhicharesubsonicbutstill B−V grated CO emission up to W = 12K kms−1, with 12 K peak supra-thermalat2.7(cid:48)resolutionwillshowspatialorvelocitysub- CO temperatures, which are very bright lines indeed. Comparison structurewhenmappedathigherresolution. ofreddeningofthestar(0.32mag)andthroughtheMilkyWay Webrieflydiscussedsomemodellingresultsoftheformation (E∞ = 0.55 mag) in the same direction, and comparison of H and excitation of CO in diffuse media (Fig. 9). The brightness B−V I seen toward the star in uv absorption and around the star in of the strongest CO lines can be understood by noting that the 21cmemission(Fig.2),suggeststhat≈0.2magofE∞ should H columndensitiesinthegasaroundζ Opharenearthepoint B−V 2 beascribedtounrelatedbackgroundmaterial. (N(H ) ≈ 1021 cm−2) where carbon fully recombines to CO at 2 CO emission from diffuse L121 gas seen around ζ Oph is evenmodestdensitiesn(H)=200−500cm−3,sothatsubstantial bimodal in E∞ , clustering around E∞ = 0.4 mag and 0.65 columnsofCOmayformingaswhichisattypicaldiffusecloud B−V B−V mag and varying widely (1 K kms−1 < W < 6-12 K kms−1) temperatures (above 30 K). In turn, such densities will excite CO withE∞ inoneoftwonarrowranges,seeFig.3:thesamelarge CO to the required degree even though they are far too low to B−V scatterinCOcolumndensitywithreddeningandN(H)whichis thermalizethelowerrotationallevelpopulations. seen globally in galactic absorption line surveys also occurs in In the range WCO <∼ 10 K kms−1 it is the case that WCO ∝ the single region studied here. The lower-reddening branch of N(CO), as a consequence of very sub-thermal excitation. The the emission is spatially segregated to the galactic northeast of CO-H conversion factor therefore varies widely in diffuse gas 2 the star and has somewhat smaller mean integrated brightness (becauseN(CO)variesrapidlywithN(H )andwithgreatscat- 2 W (Fig. 4) and velocity dispersion σ (0.35 kms−1 vs. 0.42 ter), but it takes on values N(H )/W = 2×1020H cm−2/(K CO v 2 CO 2 kms−1;Fig.6).However,peaktemperatures8-12Karepresent kms−1)inlimitedcircumstances(Fig.9panelaatupperleft). inbothbranches. If W ∝ N(CO), a sky map of W like that in Fig. 1 is CO CO Thetwomoststrikingobservationalresultsofthisstudyare amapoftheinterstellardiffusecloudchemistry.Thisshouldbe thestronglineswhichemanatefromthediffusegas,uptoW contrastedwiththemoreusualassumptionofaconstantCOcon- CO =12Kkms−1,andthevelocitystructurepresentinthestrongly version factor N(H )/W , in which case a CO map traces the 2 CO emitting 12CO line cores there (Figs. 7 and 8). The turbulent contoursofthemass(i.e.thebulkmolecularmaterial)largelyin- 10 Liszt,PetyandTachihara:CO,reddeningandturbulencearoundζOph dependentofphysicalconditions.Eventhedarkergasseeninthe —.1997,A&A,322,962 L204 complex South of ζ Oph is not immune to the influences —.2007a,A&A,476,291 ofchemistry,whichareclearlyvisibleinthedisplacementsbe- —.2007b,A&A,461,205 tweenregionsofhigherW andE∞ . Liszt,H.S.&Lucas,R.1998,A&A,339,561 CO B−V Maier,J.P.,Lakin,N.M.,Walker,G.A.H.,&Bohlender,D.A.2001,ApJ,553, The wealth of structure seen in the foreground CO bright- 267 ness map has important implications for absorption line study Morton,D.C.1975,ApJ,197,85 ofdiffuseclouds.Sincetheviewinggeometryisanaccidentof Pety,J.&Falgarone,E.2003,A&A,412,417 Pety,J.,Lucas,R.,&Liszt,H.2008,A&A,inpress the relative locations of the Sun and ζ Oph, our line of sight Rachford,B.L.,Snow,T.P.,Tumlinson,J.,Shull,J.M.,Blair,W.P.,Ferlet, to the star could equally well have intercepted any of the wide R.,Friedman,S.D.,Gry,C.,Jenkins,E.B.,Morton,D.C.,Savage,B.D., varietyofprofilesshowninFig.7.Thisbeliesourabilitytoin- Sonnentrucker,P.,Vidal-Madjar,A.,Welty,D.E.,&York,D.G.2002,ApJ, ferthegeneralpropertiesofthehostgasfromstudiesalongany 577,221 single sightline, no matter how comprehensive; studies of indi- Savage,B.D.,Drake,J.F.,Budich,W.,&Bohlin,R.C.1977,ApJ,216,291 Schlegel,D.J.,Finkbeiner,D.P.,&Davis,M.1998,ApJ,500,525 vidual absorption sightlines must be viewed demographically, Sonnentrucker,P.,Welty,D.E.,Thorburn,J.A.,&York,D.G.2007,Astrophys. asdatapointswithinalargefamilyofpossibleoutcomes,evenin J.,Suppl.Ser.,168,58 nominallysimilarconditions.Evenbeyondthis,therearesome Tachihara,K.,Abe,R.,Onishi,T.,Mizuno,A.,&Fukui,Y.2000,Publ.Astron. obviousfundamentallimitstotheuseofabsorptionlinestode- Soc.Jpn.,52,1147 VanDishoeck,E.F.&Black,J.H.1986,Astrophys.J.,Suppl.Ser.,62,109 rivethepropertiesoftheinterveninggas.Inthepresentcaseonly —.1988,ApJ,334,771 amapcouldcorrectthefalseimpressionthatthestarisocculted Wannier,P.,Penprase,B.E.,&Andersson,B.-G.1997,ApJ,487,L165 bytwoseparateforegroundclouds.Likewise,theturbulentflows Wannier,P.G.,Penzias,A.A.,&Jenkins,E.B.1982,ApJ,254,100 in the foreground gas appear clearly in maps of the gas but not Wolfire,M.G.,Hollenbach,D.,McKee,C.F.,Tielens,A.G.G.M.,&Bakes,E. atmostindividualpositions,therebyconveyingthefalseimpres- L.O.1995,ApJ,443,152 Wolfire,M.G.,McKee,C.F.,Hollenbach,D.,&Tielens,A.G.G.M.2003,ApJ, sionofanoverly-quiescentmedium. 587,278 We intend to map with higher angular resolution some re- Wood,K.,Haffner,L.M.,Reynolds,R.J.,Mathis,J.S.,&Madsen,G.2005, gionsoftheL121gaswhoselineprofilesarethermalatthe2.7(cid:48) ApJ,633,295 NANTENresolution,toseewhatkindofsubstructuremightbe presentwhensomeformsofline-broadeningareabsent.Wewill alsomapsomeL121gaswhoseprofilesarebroaderandwhose velocitygradientsarenotfullyresolvedinFig.7at4(cid:48) (0.16pc) beam-spacing. This is the second paper in a loose series (see Pety et al. (2008))whichwillalsoreportobservationsofsimilarkinemat- ics on smaller angular scales and at higher angular resolution 6(cid:48)(cid:48)-22(cid:48)(cid:48)inotherdiffusecloudsofunknowndistancewhosepres- encewasfirstmanifestedinourmm-waveabsorptionstudiesof polyatomicmolecules.Inasubsequentpaperwewilldiscussthe geometricalandphysicalinterpretationoftheinternalstructures responsiblefortheflowsseenthereandinFig.7. Acknowledgements. The National Radio Astronomy Observatory is operated by Associated Universites, Inc. under a cooperative agreement with the US NationalScienceFoundation.IRAMisoperatedbyCNRS(France),theMPG (Germany)andtheIGN(Spain).ThisresearchmadeuseoftheSimbadastro- nomicaldatabaseandtheNASAADSastrophysicaldatabasesystem.Thiswork profitedfromdiscussionsofCOexcitationwithMichelGuelin. References Black,J.H.&Dalgarno,A.1977,Astrophys.J.,Suppl.Ser.,34,405 Burgh,E.B.,France,K.,&McCandliss,S.R.2007,ApJ,658,446 Cardelli,J.A.,Mathis,J.S.,Ebbets,D.C.,&Savage,B.D.1993,ApJ,402,L17 Crutcher,R.M.1979,ApJ,231,L151 Dobashi,K.,Uehara,H.,Kandori,R.,Sakurai,T.,Kaiden,M.,Umemoto,T.,& Sato,F.2005,Publ.Astron.Soc.Jpn.,57,S1 Falgarone,E.,Panis,J.-F.,Heithausen,A.,Perault,M.,Stutzki,J.,Puget,J.-L., &Bensch,F.1998,A&A,331,669 Falgarone,E.&Phillips,T.G.1996,ApJ,472,191 Gaustad,J.E.,McCullough,P.R.,Rosing,W.,&VanBuren,D.2001,Publ. Astron.Soc.Pac.,113,1326 Goldreich,P.&Kwan,J.1974,ApJ,189,441 Hartmann, D. & Burton, W. B. 1997, Atlas of galactic neutral hydrogen (Cambridge;NewYork:CambridgeUniversityPress) Herbig,G.H.1968,ZeitschriftfurAstrophysik,68,243 Knapp,G.R.&Jura,M.1976,ApJ,209,782 Kopp,M.,Gerin,M.,Roueff,E.,&LeBourlot,J.1996,A&A,305,558 Lambert,D.L.,Sheffer,Y.,Gilliland,R.L.,&Federman,S.R.1994,ApJ,420, 756 Langer,W.D.,Glassgold,A.E.,&Wilson,R.W.1987,ApJ,322,450 Liszt,H.,Lucas,R.,&Pety,J.2006,A&A,448,253 Liszt,H.S.1979,ApJ,233,L147