-,4 (NASA-CR-195940) AIRBORNE N94-34705 ASTRONOMY WITH A 150 MICROMETER - 500 MICROMETER HETEROOYNE SPECTROMETER Fina| Report, I Oct. Unclas 1983 - 30 Sep. 1991 (California Univ.) 129 p G3/89 0010009 b • i j / J 'Y _ NASA-CR-195960 e / ./ t/ /a p FINAL REPORT NASA GRANT NAG2-254 AIRBORNE ASTRONOMY WITH A 150gm - 500gm HETERODYNE SPECTROMETER Period: October 1, 1983 through September 30, 1991 A. L. Betz, Principal Investigator Space Sciences Laboratory, University of California, Berkeley ABSTRACT This report summarizes work done under NASA Grant NAG2-254 awarded to the University of California. The project goal was to build a far-infrared heterodyne spec- trometer for NASA's Kuiper Airborne Observatory, and to use this instrument to observe atomic and molecular spectral lines from the interstellar medium. This goal was success- fully achieved; the spectrometer is now in routine use aboard the KAO. Detections of particular note have been the 370gm line of neutral atomic carbon, the 158 l.tm transition of ionized carbon, many of the high-J rotational lines of 12CO and 13CO between J=9-8 and J=22-21, the 119 gm ground-state rotational line of OH, and the 219 gm ground-state rotational line of H2D +. All of these lines were observed at spectral resolutions exceeding 1 part in 106, thereby allowing accurate line shapes and Doppler velocities to be meas- ured. I. Introduction The following synopsis relies on the detailed presentations available from cited pub- lications included here as Appendices. A list of corresponding publications is included on page 13. References to individual entries are noted in the text by a number in brackets (e.g., [1]). The text below begins with a technical review of the spectrometer and then details the specific astronomical observations accomplished with this instrument. lI. Technical Overview Our FIR heterodyne spectrometer has the 4 basic subsystems characteristic of most heterodyne instruments: a local oscillator (LO), a mixer, an intermediate-frequency (IF) amplifier, and a multi-channel radio-frequency signal analyzer (either a filter bank or an acousto-optic spectrometer: AOS). The LO is an optically-pumped laser [1], and the mixer is a cooled Schottky-diode in a special mount [6]. These components, together with the IF amplifier, are the critical "front-end" components of the spectrometer, and they are the only ones actually mounted on the telescope. The "back-end" components such as the filter bank are mounted in the experimenters' rack. They are non-critical in that they affect the resolution but not the sensitivity of the spectrometer. For this reason, a heterodyne spectrometer has a very 18 Igg4 _,_ practical advantage over an incoherent instrument for high resolution spectroscopy. A heterodyne spectrometer does its spectral analysis after "detection" whereas a spectrom- eter using incoherent detection does its resolving before detection. In this latter case, the inevitable losses and extraneous emissions of the spectrometer reduce the signal-to-noise ratio from that obtainable with a perfect instrument. On the other hand, the coherent (heterodyne) spectrometer does have an additional but fundamental source of noise from quantum fluctuations. Realistically however, the quantum-noise level of a heterodyne receiver even at _,= 100ktm is only TQ= 140 K, or about 10-21WHz -1/2, and thus is negligible compared to the non-fundamental noise of current receivers. Perhaps the most obvious advantage of heterodyne spectroscopy is its capability for ultra-high resolution. A velocity resolution <1 km/s can easily be achieved with a velocity-scale accuracy lim- ited only by the knowledge of the laser and source line-frequencies. The front-end components of our airborne spectrometer weigh 200 lbs. and have a moment of 800 ft-lbs, measured from the center of the air bearing on the KAO. The size of the laser-frame dictates the overall size of the front-end hardware. The frame is 40" L x 13.5" W x 8" H and consists of four 1-inch O.D. invar rods that separate 2 insulating granite end-plates. A 40" x 13.5" aluminum plate is attached to the top of the laser frame, and on it are positioned all the optical components of the spectrometer. The entire assembly fits halfway into the instrument cavity at the Naysmith focus position. The other half of the spectrometer protrudes beyond the position of the standard instrument flange, but is braced to the instrument cavity. The standard instrument-mounting flange is removed, of course. The layout of the spectrometer is illustrated in Figure 1 of Appendix A (also refer- ence [1]), which also contains other details of the instrument's configuration and opera- tion. III. Scientific Program (1) FY85 Results At the time of the instrument's observational debut in 1985, there were no high resolution spectrometers aboard the KAO that operated at wavelengths between 200 lxm and 500 ptm. Consequently, this unexplored spectral region was ripe for the detection of many new and important lines. Ground-based observations are for the most part impossi- ble because of obscuration by water vapor in the Earth's atmosphere. The new instrument was successfully inaugurated during two flights in July, 1985, when the 373 lxm fine structure line of atomic carbon was detected in four sources [2]. The results, presented in detail in Appendix B, allow us to derive useful estimates for the excitation and abun- dance of neutral carbon in dense molecular clouds. Since these first observations, the sen- sitivity and wavelength coverage of the instrument have been continually improved. (2) FY86 Results Neutral Atomic Carbon: The 373 _tm Fine Structure Line In FY86, our second year on the KAO, we continued our survey of the 810 GHz (373 I.tm) 3p2-3p 1fine structure line of neutral atomic carbon in dense molecular clouds. Our J = 2-1 results from a year earlier indicated that neutral carbon is very abundant, as originally shown by the Caltech group from observations of the 610 _tm C I J=l-0 line. -2- u" Our analysis, however, indicated that while carbon was abundant, it was probably not optically thick in the fine structure lines, since excitation temperatures between 30 and 60 K were indicated by the observed line intensity ratios. Observations of the J = 2-1 line in additional sources previously detected in the C I J=l-0 line, especially Orion, were needed to confirm the general validity of our conclusions. Three flights were scheduled for the end of the 1986 observing season for this purpose. (Earlier winter flights were not available.) The flights in September, 1986, went very well, and new data on the J=2-1 CI line at selected locations in M42, NGC 2024, S140, DR21, and W3 were obtained to complement our earlier results [4]. Figure 1 in Appendix D shows new J=2-1 spectra in 4 of the above mentioned sources. The 1986 results show a factor of 3 improvement in sensitivity over our 1985 detections. In addition a limited amount of mapping of J=2-1 C I distribution was done for each of the new sources. Our 1986 results on CI generally confirm our earlier analysis but with improved limits, and the updated conclusions were reported at the January, 1987, AAS meeting in Pasadena. The observations and analysis of the CI data from 1985 and 1986 formed a major part of the Ph.D. thesis dissertation of Jonas Zmuidzinas [17]. (3) FY87 Results Singly Ionized Carbon: the 158 Ixm Fine Structure Line In FY87 our work concentrated on observations of the 158 p.m fine structure line of C II at the edges of dense molecular clouds. Our purpose is to study the excitation, abun- dance, and dynamics of gas in photodissociation regions which lie between neutral molecular clouds and the H II regions surrounding nearby O- and B-type stars. In early February 1987, we observed the CII fine structure line in M42, NGC 2024, and W3. These observations were the first to resolve the narrow C II line completely, and to open up the short wavelength (_,<200 Ixm) part of the FIR _pectrum to heterodyne spectros- copy. Our results on M42, in particular, at 0.8 km s- resolution showed that the C II emission has multiple velocity components at some positions, and noticeable changes in line shape and velocity between the BN and Trapezium regions [3]. Several C II com- ponents correlate with C+ recombination line features observed at microwave frequen- cies, while others appear at known velocities of molecular material such as CO. The results, partially presented in Appendices C, I, and K, allow us to derive useful estimates for the excitation, abundance, and velocity distribution of ionized carbon that can be readily compared with similarly accurate microwave data. Following these first hetero- dyne observations, we have gone on to detect C II emission from more than 100 positions in over 40 sources observed. A thorough analysis requires a detailed comparison with carbon recombination line data (when available) and isotopic CO data from microwave observations. Such a com- parison can readily be made because our C II line profiles are measured at the same high spectral resolution (0.8 km s-1) as the microwave data. In general we find that the C II line profiles are very similar to those of J=l-0 and J=2-1 12CO. These CO lines are usu- ally optically thick in the clouds observed and thus only provide information on the near- side cloud boundary. The C II lines are thought to arise from photodissociation regions at these same cloud boundaries. The intensity of C II appears to correlate with the UV flux "- . impinging on the peripheral molecular material. The C II emission is seen to peak on the _-.._loud boundary close to the adjacent H II region whenever the geometry permits a view -3- perpendicularto thedirectionof UV illumination(e.g.,NGC7538,CepB, S140,M17) [11]. Also for agivensource,thevelocitiesofthestrongerCII linecomponents(when thereareseveral)agreewith thoseof thecarbonrecombinationlines,whichin general havevelocitiessimilartothemolecularmaterial(CO)ratherthanthehydrogenrecombi- nation lines. This fact is additionalevidencefor locating the CII in the molecular cloud/HII interfaceregion,asexpected. An interestingphenomenonobservedin a numberof sourcesis the presenceof apparent"self-reversal"intheCII lineemission.Asshownin AppendixI,thespectraof W49 andW51appeartoshowtwovelocitycomponentsw, henin factsomeoftheemis- sionnearline centerhasbeenabsorbedby foregroundC÷in thegroundstate.Similar "self-reversed"profilesareseenin thelow-J COlinesfor thesesources.Othersources showingCII self-reversalsareW3,NGC2024,andtoalesserextentMon R2andM17. We cannotgivea completeexplanationhere,butthegeneralconclusionis thationized carbon in a foregroundcloud or cloud periphery exists in low density regions (nil2< 103cm-3)wherethecollisionalexcitationrateisinsufficienttomaintainthepopu- lation of theupperenergylevelin LTE. Consequentlyw, henthevelocitydispersionof theforegroundmaterialis lessthanthatof thebackgroundemittinggas,weseeabsorp- tion. We knowthatthereis coolerforegroundmaterialin the sourcesmentionedfrom microwavedataonCOandH2CO. Theabsorptioncentersof self-reversedlinesarefrequentlyopticallythick.Thisfact is self-evidentfromtheabsorptiondepths.ForW51,theminimumopticaldepthof CII in absorptionis - 2.0, and and the maximum excitation temperature of the absorbing foreground gas is 29 K. The gas kinetic temperature is higher, and our preliminary cal- culations indicate that in the absorbing cloud Tkin> 40 K, if nFi2< 500 cm-. In W51 the pathlength through the absorbing foreground cloud is -2 parsecs, given the C II column density of 1018 cm -2 we calculate from the integrated strength of the absorption. In such diffuse regions, C U line radiation is thought to be the dominant cooling mechanism; however, when the line is optically thick, then the cloud temperature must rise and CO emission becomes more significant. The fact that all C II sources do not show self-reversals indicates that morphological effects must be taken into consideration. For example, if there is a sufficient velocity gra- dient (1 linewidth) between the background CII emitter and the foreground gas, then no absorption can occur. Note that the sources for which self-absorption is seen generally have relatively wide emission lines > 10 km s-1, which makes a larger "quasi-continuum" available for subsequent absorption. In front- and side-UV-illuminated sources we don't see much evidence for self-absorption; whereas in the interior- or back-illuminated sources we see absorption just like that seen in J=l-0 12CO. The profiles of the C II lines generally follow those of CO; however, from the 12CII emission profiles alone we cannot conclude with certainty whether C II is itself optically thick like CO. Other fine structure lines, such as those from neutral oxygen, need to be observed at similarly high spectral resolution before we can make definitive statements on the relative importance of various cooling processes in photodissociation regions. -4- (4)FY88Results Shock-excited CO: the 153 _tm J=l 7-16 Line In January 1988, we observed the 153 l.tm (1956.018 GHz) J= 17-16 transition of 12CO at 5 positions in the IRc2 region of the Orion molecular cloud. The results, shown as Figure 2 of Appendix E, are the first spectra of a far-infrared CO line taken with sufficient resolution to not only show the lineshape, but also to provide an accurately calibrated Doppler velocity scale. The far-infrared CO lines can now be compared directly with millimeter and submillimeter CO data taken at similarly high resolution of <1 km s-1. A full account of the work has now been published [5] and is reproduced here as Appendix E. We'll just summarize the two principal conclusions here. At 0.8 km s-l resolution we can distin__uish 2 velocity components in the spectra: a broad component with FWHM - 30 km s- which is very similar to the plateau com- ponent seen in millimeter-wave lines, plus a second narrower component with FWHM < 10 km s-1 that seems dynamically and spatially characteristic of the quiescent molecular cloud. The narrow component is strongest at the Trapezium and suggests that UV- heating may provide some of the excitation for high-J CO emission. The existence of such hot quiescent gas is quite interesting, and heretofore unknown. We also see weak narrow h•ne emls•ston• (2 km sml (FWHM)) of this kind in our CO J= 17-16 spectra of AFGL 490. Appendix J is another reference to this work that we presented at the Inter- national Conference on Millimeter and Submillimeter-wave Astronomy, held in Kona, HI [10]. Figure 2 of this latter paper shows a comparison of the peak brightness tempera- tures of CO emission in the plateau component measured with heterodyne instruments all the way up to J = 17-16. A single-component fit to the data indicates that essentially all of the CO plateau emission can be characterized by an excitation temperature of 180 + 50 K over this range of J-values. The question remains, however, whether even higher-J CO emission seen in Orion has a plateau component which can be similarly represented, and whether it is dynamically related to the low-J emission. In addition, we need to under- stand the significance of the narrow component CO emission we see in Orion. During the CO observations we also made a brief search for shock-excited emission from H20. The 954-- 863 line of H20 at 1969.223 GHz was searched for in Orion without success. This transition, predicted to be a very weak maser in the calculations of Neufeld and Melnick in 1987, was nevertheless close enough to the CO J=17-16 line that we could do it at the same time. The failure to detect it places no strong constraints on either the abundance or excitation of H20. (5) FY89 Results (A) Shock-Excited CO: the 119 Ftm J = 22-21 Line In January 1989 we observed the J=22-21 line of 12CO at 2528 GHz (l18.8[tm) with our heterodyne spectrometer. The spectra obtained on the IRc2 region of Orion are the highest frequency heterodyne data obtained so far on the KAO. The high 0.6 km s-l resolution of the spectra allow us to resolve the line clearly. We distinguish 2 velocity components in the high-J CO emission: a broad component with a 35 km s-1 line width (FWHM), and a narrow component with an 8 km s-1 width (FWHM). The 2 components can be unambiguously associated with those seen at millimeter wavelengths - the wide one with the outflow gas from IRc2 and the narrow with the -5- quiescentmolecularcloud. Theseresultsarethe first to demonstratethat the high-J infraredemissionof COcomesfromthesamekinematiccomponentsofgasasthelow-J emission,whichisthoughttoariseprincipallyfrommuchcoolermaterial. By comparingtheseobservationswith dataon lower-Jtransitions,including the J=17-16line which we observedat similarly high resolutionin KAO observationsof January1988,weseethattheJ=22-21emissioncomesfromtheshock-heatedinterface betweenthe outflowing andquiescentgascomponents.While this effect has been predictedby magnetohydrodynamimc odelsof interstellarshocks,the newdataprovide the onlydirectevidencefor thepredictedheatingof quiescentgasto temperaturesnear 600K just beforeit is sweptinto thestreamof outflowinggas. A reprintof theAp. J. Letterspublication[7]describingthisworkisincludedasAppendixG. (B) OH Fundamental Line at 119 _tm in the Orion-lRc2 Region On our second flight in January 1989 we observed the 2514 GHz (119.2 gm) line of OH at 0.6 km s-1 resolution in the same region of Orion. From observations at the IRc2, H2 peak-1 and H2 peak-2 positions, we find the emission to be quite compact, with a size <25". The relatively small size means that the OH emission, although associated with the outflow from IRc2, is not primarily produced in the gas immediately behind the shock-front interface, which has a - 60"diameter. Most theoretical models of the shock region associate both the far-infrared OH and high-J CO emission with the shock-front. However, we see now that only the CO emission is sufficiently extended for this associa- tion to be valid. Another interesting aspect of the OH emission at the IRc2 position is that apparently the entire blue-shifted half of the emission line profile for this ground-state OH transition is absorbed by less excited foreground gas. An additional sharp absorber is present at 5.5 km s-1 VLSR, the characteristic velocity of the "hot core" source. This is the first time the "hot core" has been seen in absorption, and naturally implies that a significant part of the "hot core" lies in the foreground of a more compact OH emission region, and in fact in front of the dust responsible for the FIR continuum. Our analysis places most of the OH emitting region within 5"of IRc2, and most of the absorbing gas more than 5 times further away. The "self-absorption" of half the l l91xm line intensity explains an apparent discrepancy noticed between the relative intensities of the fundamental and various excited-state lines: namely, that the OH 163 _tm/119 l.tm integrated intensity ratio was more than a factor of 2 larger than allowed even by gas at infinite temperature. Now we see that is not the case, and something near LTE likely prevails. Here again high resolu- tion makes a big difference to the interpretation of data; new phenomena are visible when you have the ability to see them. A full account of the observations and analysis may be found in reference [8], reproduced here as Appendix H. (C) 12CO and laCO J=9-8 Observations In August 1989 we observed the 12CO and 13CO J=9-8 lines (at 1037 and 991 GHz, respectively) in a variety of galactic molecular clouds that have been extensively studied in lower-J transitions. The purpose of these 300 l.tm CO observations was to look for high excitation, probably low optical depth, gas that would not be evident in the low-J data. Any such gas would be a likely indicator of UV-excitation from embedded young -6- stars.Whatwefoundwassomethingquitedifferent- averagedoverour 1.4"beammuch ..... 12 13 of the emission in the J=9-8 trans_tlons of CO and CO appears to originate in the dense regions of the molecular cloud with the same moderate excitation temperatures (Tex - 50 K) seen in the low-J (J < 3) CO lines. Our 13CO observations are the first detections of this line in any source, and the first detection of a resolved 13CO line above J=2-1. With J=9-8 data from both the 12CO and 13CO isotopic species taken under identical conditions, we are able to deduce the optical depth of the 12CO J = 9-8 line independently of the excitation conditions. We find the latter line to be significantly optically thick in most of the sources observed. This result is a somewhat surprising given the relatively low line intensities that were observed: Ta* - 10-20 K. The gas must be either cold, significantly clumped, or more likely both. A non-LTE analysis (including line opacity effects) of the 12CO and 13CO J=9-8 data in the context of previous observations of lower-J lines suggests that much of the high-J CO emission arises from the bulk of the molecular cloud material. However, most of the emitting region is clumped, with beam filling factors <1 for the high density gas producing the 13CO emission. Although some additional emission from warmer gas near 100 K is also evident in the 12CO data, it is optically thin and very beam diluted. It is interesting to note that in some sources the J=9-8 lines of 12CO and 13CO are probably produced under quite different physical conditions, as a consequence of the high densi- ties or radiation levels required to maintain level populations at high-J. In the J=9-8 12CO line we do not see the intense extended emission reported by other observers for the J=7-6 12CO line. More details on the J=9-8 results are to be found in reference [13] included here as Appendix M. (D) C H Data Analysis In January 1988 we completed our initial survey of 158 lxm C II line emission in 22 galactic molecular clouds. C II was detected from 40 positions in 17 sources, and the line profiles observed at 0.8 km s-1 resolution were compared with available CO J = 2-1 and 1-0 data. The results on the dynamics of photodissociated gas have been divided into three categories, which are being analyzed separately. The first, on C II emission from side-illuminated molecular clouds, has been completed [11], and is included here as Appendix K. The reprint of the second manuscript [9], on sources which show "self- reversed" C II line profiles, is included here as Appendix I. (6) FY90 Results (A) CH in the Large Magellanic Cloud In May, 1990, we participated in the FY90 Southern Skies Expedition of the KAO. The primary goal of our successful 2-flight series was to observe the 158 lxm C II fine structure line near a number of identifiable H 17regions in the LMC, and to correlate any observed emission with low-J CO emission available from ground-based observations. Because we are interested in dynamical as well as intensity correlations, the high spectral resolution of the infrared heterodyne receiver, comparable to that of a millimeter wave receiver, is essential. The LMC is known to have a low metallicity relative to the Galaxy. Therefore, the UV flux from newly formed stars can penetrate far deeper into molecular clouds before being attenuated by dust. The previously observed low intensity of millimeter-wave CO -7- radiationfrom theLMC hasbeenattributedbothtotherelativelylow abundanceof C andO, andalsoto theenhancedphotodestructionof unshieldedCOby UV radiation. Conversely,theCII radiationcouldbeexpectedtoberelativelymoreintensebecauseof thegreaterpenetrationdepthof ionizingradiation,andhence,thecolumndensityofCII shouldbehigheraroundahotstarintheLMC thaninourgalaxy. OurobservationsshowedthattheCII emittingregionsarewidespreadin theLMC, aswouldbeexpectedfromtherelativepaucityofdust,andthattheintensityofradiation in a givenregionis comparableto (or slightly lessthan)thatwhichwouldbeobserved from aGalacticgiantmolecularcloudatthesamedistance[14]. Thiscomparisonsug- geststhatC+is eithercolderor lessabundantin theLMC thanin our owngalaxy,on average;orthatCII ismarginallyopticallythickor partiallyself-absorbedi,n contrastto currentpredictionsof photodissociationregions.The spectraillustratedin Figure2 of Appendix N showexamplesof our Cil resultsfrom the 30 Doradusregion, oneof numerouspositionsweobservedin theLMC. HereCII emissionis strongattheposi- tionsof the two far-infraredcontinuumpeaksthatstraddlethecentralexcitationsource R136. COemissionfromtheseregionsof UV-heateddustisrelativelyweak. However, no CII emissionwasdetectedatthelocationswhich showthestrongestCOmillimeter emissionin severalLMC molecularcloudregions,in contrasttothe usualcasein our ownGalaxy. The anticorrelationprovidesdirectevidencefor thehypothesisthatUV- photodestructionof CO is themajorreasonwhy COlines areweakin the LMC. The anticorrelationdoesnot supportacurrentmodelof CO radiationfrom galaxieswhich arguesthatmostCOemission,evenfromlow-Jlines,arisesfromthewarmgasaround photodissociationregions.In thecaseoftheLMC, thelow dustabundanceappearstobe themostsignificantdeterminantof bothCOandC13radiation.Thisconclusioncouldbe testedwith futureobservationsof CII andCOin the SMC,a companiongalaxywith evenlowermetallicity. (B) CH in G333.6-0.2 The other major C II source observed during the Southern Hemisphere deployment was the immense RCW 106 complex, which includes the giant H II region G333.6-0.2, which is the most intense C I3 and H2 source in the sky. A total of 11 positions in RCW 106 were observed to provide an overall assessment of the correlation of C II and molec- ular gas dynamics and excitation. The relatively large line width of the emission makes self-absorption from less excited foreground material readily apparent. Similar self- absorption is seen in just about all the wide-line (>15 km s-1 FWHM) sources we have observed so far in the galaxy. The ability of the spectrometer to resolve line profiles allows us to detect absorption components of significant column density (and optical depth) that would otherwise be missed. There are interesting changes in line shapes with position, presumably because of "self-absorption" effects at the 1" offset position. The complex spectra observed at various positions will take some time to analyze fully, but a manuscript detailing the initial results is currently in preparation. It also includes data taken in subsequent years on the J= 9-8, 12-11, and 14-13 transitions of CO. (C) Guest Investigator Flights During FY90 two flights were scheduled for guest investigators. Because of the air- craft delays of March 1990 both flights were carried-over to FY91. Problems with the aircraft cavity-wall in Oct., 1990, and with the OSM (chopper) mechanism in Dec., 1990 -8-