The depth-dependent abundance of both gas-phase and solid-state water within dense, quiescent, molecular clouds is important to both the cloud chemistry and gas cooling. Where water is in the gas phase, it's free to participate in the network of ion-neutral reactions that lead to a host of oxygen-bearing molecules, and its many ortho and para energy levels make it an effective coolant for gas temperatures greater than 20K. Where water is abundant as ice on grain surfaces, and unavailable to cool the gas, significant amounts of oxygen are removed from the gas phase, suppressing the gas-phase chemical reactions that lead to a number of oxygen-bearing species, including O2. Models of FUV-illuminated clouds predict that the gas-phase water abundance peaks in the range Av ~3 and 8mag of the cloud surface, depending on the gas density and FUV field strength. Deeper within such clouds, water is predicted to exist mainly as ice on grain surfaces. More broadly, these models are used to analyze a variety of other regions, including outflow cavities associated with young stellar objects and the surface layers of protoplanetary disks. In this paper, we report the results of observational tests of FUV-illuminated cloud models toward the Orion Molecular Ridge and Cepheus B using data obtained from the Herschel Space Observatory and the Five College Radio Astronomy Observatory. Toward Orion, 2220 spatial positions were observed along the face-on Orion Ridge in the H2O 110-101 557GHz and NH3 J,K=1,0-0,0 572GHz lines. Toward Cepheus B, two strip scans were made in the same lines across the edge-on ionization front. These new observations demonstrate that gas-phase water exists primarily within a few magnitudes of dense cloud surfaces, strengthening the conclusions of an earlier study based on a much smaller data set, and indirectly supports the prediction that water ice is quite abundant in dense clouds.

We present observations of the Trumpler 14/Carina I region carried out using the Stratospheric Terahertz Observatory 2 (STO2). The Trumpler 14/Carina I region is in the west part of the Carina Nebula Complex, which is one of the most extreme star-forming regions in the Milky Way. We observed Trumpler 14/Carina I in the 158 $\mu$m transition of [C\,{\sc ii}] with a spatial resolution of 48$''$ and a velocity resolution of 0.17 km s$^{-1}$. The observations cover a 0.25$^\circ$ by 0.28$^\circ$ area with central position {\it l} = 297.34$^\circ$, {\it b} = -0.60$^\circ$. The kinematics show that bright [C\,{\sc ii}] structures are spatially and spectrally correlated with the surfaces of CO clouds, tracing the photodissociation region and ionization front of each molecular cloud. Along 7 lines of sight that traverse Tr 14 into the dark ridge to the southwest, we find that the [C\,{\sc ii}] luminosity from the HII region is 3.7 times that from the PDR. In same los we find in the PDRs an average ratio of 1:4.1:5.6 for the mass in atomic gas:dark-CO gas: molecular gas traced by CO. Comparing multiple gas tracers including HI 21cm, [C\,{\sc ii}], CO, and radio recombination lines, we find that the HII regions of the Carina Nebula Complex are well-described as HII regions with one-side freely expanding towards us, consistent with the champagne model of ionized gas evolution. The dispersal of the GMC in this region is dominated by EUV photoevaporation; the dispersal timescale is 20-30 Myr.

Recent papers have found that the inferred slope of the high-mass ($>1.5$ M$_\odot$) IMF for field stars in the solar vicinity has a larger value ($\sim 1.7-2.1$) than the slopes ($\sim 1.2-1.7$; Salpeter= 1.35) inferred from numerous studies of young clusters. We attempt to reconcile this apparent contradiction. Stars mostly form in Giant Molecular Clouds, and the more massive stars ($\gtrsim 3$ M$_\odot$) may have insufficient time before their deaths to uniformly populate the solar circle of the Galaxy. We examine the effect of small sample volumes on the {\it apparent} slope, $\Gamma_{\rm app}$, of the high-mass IMF by modeling the present day mass function (PDMF) over the mass range $1.5-6$ M$_\odot$. Depending on the location of the observer along the solar circle and the size of the sample volume, the apparent slope of the IMF can show a wide variance, with typical values steeper than the underlying universal value $\Gamma$. We show, for example, that the PDMFs observed in a small (radius $\sim 200$ pc) volume randomly placed at the solar circle have a $\sim 15-30$\% likelihood of resulting in $\Gamma_{\rm app} \gtrsim \Gamma+ 0.35$ because of inhomogeneities in the surface densities of more massive stars. If we add the a priori knowledge that the Sun currently lies in an interarm region, where the star formation rate is lower than the average at the solar circle, we find an even higher likelihood ($\sim 50-60\%$ ) of $\Gamma_{\rm app} \gtrsim \Gamma+0.35$, corresponding to $\Gamma_{\rm app} \gtrsim 1.7$ when the underlying $\Gamma= 1.35$.

We present high resolution ($0.11"\times0.06"$) 3mm ALMA observations of the highly inclined transition disk around the star T Cha. Our continuum image reveals multiple dust structures: an inner disk, a spatially resolved dust gap, and an outer ring. When fitting sky-brightness models to the real component of the 3mm visibilities, we infer that the inner emission is compact ($\le1$au in radius), the gap width is between 18-28 au, and the emission from the outer ring peaks at $\sim36$ au. We compare our ALMA image with previously published 1.6$\mu$m VLT/SPHERE imagery. This comparison reveals that the location of the outer ring is wavelength dependent. More specifically, the peak emission of the 3mm ring is at a larger radial distance than that of the 1.6$\mu$m ring, suggesting that millimeter-sized grains in the outer disk are located further away from the central star than micron-sized grains. We discuss different scenarios to explain our findings, including dead zones, star-driven photoevaporation, and planet-disk interactions. We find that the most likely origin of the dust gap is from an embedded planet, and estimate --- for a single planet scenario --- that T Cha's gap is carved by a $1.2M_{jup}$ planet.

Using Keck/HIRES spectra {\Delta}v ~ 7 km/s, we analyze forbidden lines of [O I] 6300 {\AA}, [O I] 5577 {\AA} and [S II] 6731 {\AA} from 33 T Tauri stars covering a range of disk evolutionary stages. After removing a high velocity component (HVC) associated with microjets, we study the properties of the low velocity component (LVC). The LVC can be attributed to slow disk winds that could be magnetically (MHD) or thermally (photoevaporative) driven. Both of these winds play an important role in the evolution and dispersal of protoplanetary material. LVC emission is seen in all 30 stars with detected [O I] but only in 2 out of eight with detected [S II] , so our analysis is largely based on the properties of the [O I] LVC. The LVC itself is resolved into broad (BC) and narrow (NC) kinematic components. Both components are found over a wide range of accretion rates and their luminosity is correlated with the accretion luminosity, but the NC is proportionately stronger than the BC in transition disks. The FWHM of both the BC and NC correlates with disk inclination, consistent with Keplerian broadening from radii of 0.05 to 0.5 AU and 0.5 to 5 AU, respectively. The velocity centroids of the BC suggest formation in an MHD disk wind, with the largest blueshifts found in sources with closer to face-on orientations. The velocity centroids of the NC however, show no dependence on disk inclination. The origin of this component is less clear and the evidence for photoevaporation is not conclusive.

We obtain high spatial and spectral resolution images of the CO J=2-1, CN N=2-1 and CS J=5-4 emission with ALMA in Cycle~2. The radial distribution of the turbulent broadening is derived with three approaches: two `direct' and one modelling. The first requires a single transition and derives \Tex{} directly from the line profile, yielding a \vturb{}. The second assumes two different molecules are co-spatial thus their relative linewidths allow for a calculation of \Tkin{} and \vturb{}. Finally we fit a parametric disk model where physical properties of the disk are described by power laws, to compare our `direct' methods with previous values. The two direct methods were limited to the outer $r > 40$~au disk due to beam smear. The direct method found \vturb{} ranging from $\approx$~\vel{130} at 40~au, dropping to $\approx$~\vel{50} in the outer disk, qualitatively recovered with the parametric model fitting. This corresponds to roughly $0.2 - 0.4~c_s$. CN was found to exhibit strong non-LTE effects outside $r \approx 140$~au, so \vturb{} was limited to within this radius. The assumption that CN and CS are co-spatial is consistent with observed linewidths only within $r \lesssim 100$~au, within which \vturb{} was found to drop from \vel{100} ($\approx~0.4~c_s$) to nothing at 100~au. The parametric model yielded a near constant \vel{50} for CS ($0.2 - 0.4~c_s$). We demonstrate that absolute flux calibration is and will be the limiting factor in all studies of turbulence using a single molecule. The magnitude of the dispersion is comparable with or below that predicted by the magneto-rotational instability theory. A more precise comparison would require to reach an absolute calibration precision of order 3\%, or to find a suitable combination of light and heavy molecules which are co-located in the disk.

We present a detailed analysis of narrow of NaI and KI absorption resonance lines toward nearly 40 T Tauri stars in Taurus with the goal of clarifying their origin. The NaI 5889.95 angstrom line is detected toward all but one source, while the weaker KI 7698.96 angstrom line in about two thirds of the sample. The similarity in their peak centroids and the significant positive correlation between their equivalent widths demonstrate that these transitions trace the same atomic gas. The absorption lines are present towards both disk and diskless young stellar objects, which excludes cold gas within the circumstellar disk as the absorbing material. A comparison of NaI and CO detections and peak centroids demonstrates that the atomic and molecular gas are not co-located, the atomic gas is more extended than the molecular gas. The width of the atomic lines corroborates this finding and points to atomic gas about an order of magnitude warmer than the molecular gas. The distribution of NaI radial velocities shows a clear spatial gradient along the length of the Taurus molecular cloud filaments. This suggests that absorption is associated with the Taurus molecular cloud. Assuming the gradient is due to cloud rotation, the rotation of the atomic gas is consistent with differential galactic rotation while the rotation of the molecular gas, although with the same rotation axis, is retrograde. Our analysis shows that narrow NaI and KI absorption resonance lines are useful tracers of the atomic envelope of molecular clouds. In line with recent findings from giant molecular clouds, our results demonstrate that the velocity fields of the atomic and molecular gas are misaligned. The angular momentum of a molecular cloud is not simply inherited from the rotating Galactic disk from which it formed but may be redistributed by cloud-cloud interactions.

The surface density and vertical distribution of stars, stellar remnants, and gas in the solar vicinity form important ingredients for understanding the star formation history of the Galaxy as well as for inferring the local density of dark matter by using stellar kinematics to probe the gravitational potential. In this paper we review the literature for these baryonic components, reanalyze data, and provide tables of the surface densities and exponential scale heights of main sequence stars, giants, brown dwarfs, and stellar remnants. We also review three components of gas (H2, HI, and HII), give their surface densities at the solar circle, and discuss their vertical distribution. We find a local total surface density of M dwarfs of 17.3 pm 2.3 Mo/pc^2. Our result for the total local surface density of visible stars, 27.0 pm 2.7 Mo/pc^2, is close to previous estimates due to a cancellation of opposing effects: more mass in M dwarfs, less mass in the others. The total local surface density in white dwarfs is 4.9 pm 0.6 Mo/pc^2; in brown dwarfs, it is ~1.2 Mo/pc^2. We find that the total local surface density of stars and stellar remnants is 33.4 pm 3 Mo/pc^2, somewhat less than previous estimates. We analyze data on 21 cm emission and absorption and obtain good agreement with recent results on the local amount of neutral atomic hydrogen obtained with the Planck satellite. The local surface density of gas is 13.7 pm 1.6 Mo/pc^2. The total baryonic mass surface density that we derive for the solar neighborhood is 47.1 pm 3.4 Mo/pc^2. Combining these results with others' measurements of the total surface density of matter within 1-1.1 kpc of the plane, we find that the local density of dark matter is 0.013 pm 0.003Mo/pc^3.The local density of all matter is 0.097 pm 0.013 Mo/pc^3. We discuss limitations on the properties of a possible thin disk of dark matter.

We present spectral data cubes of the [CI] 809GHz, 12CO 115GHz, 13CO 110GHz and HI 1.4GHz line emission from an 1 square degree region along the l = 328{\deg} (G328) sightline in the Galactic Plane. Emission arises principally from gas in three spiral arm crossings along the sight line. The distribution of the emission in the CO and [CI] lines is found to be similar, with the [CI] slightly more extended, and both are enveloped in extensive HI. Spectral line ratios per voxel in the data cubes are found to be similar across the entire extent of the Galaxy. However, towards the edges of the molecular clouds the [CI]/13CO and 12CO/13CO line ratios rise by ~50%, and the [CI]/HI ratio falls by ~10$%. We attribute this to these sightlines passing predominantly through the surfaces of photodissociation regions (PDRs), where the carbon is found mainly as C or C+, while the H2 is mostly molecular, and the proportion of atomic gas also increases. We undertake modelling of the PDR emission from low density molecular clouds excited by average interstellar radiation fields and cosmic-ray ionization to quantify this comparison, finding that depletion of sulfur and reduced PAH abundance is needed to match line fluxes and ratios. Roughly one-third of the molecular gas along the sightline is found to be associated with this surface region, where the carbon is largely not to be found in CO. ~10% of the atomic hydrogen along the sightline is cold gas within PDRs.

We have mapped the superwind/halo region of the nearby starburst galaxy M82 in the mid-infrared with $Spitzer-IRS$. The spectral regions covered include the H$_2 S(1)-S(3)$, [NeII], [NeIII] emission lines and PAH features. We estimate the total warm H$_2$ mass and the kinetic energy of the outflowing warm molecular gas to be between $M_{warm}\sim5-17\times10^6$ M$_{\odot}$ and $E_{K}\sim6-20\times10^{53}$ erg. Using the ratios of the 6.2, 7.7 and 11.3 micron PAH features in the IRS spectra, we are able to estimate the average size and ionization state of the small grains in the superwind. There are large variations in the PAH flux ratios throughout the outflow. The 11.3/7.7 and the 6.2/7.7 PAH ratios both vary by more than a factor of five across the wind region. The Northern part of the wind has a significant population of PAH's with smaller 6.2/7.7 ratios than either the starburst disk or the Southern wind, indicating that on average, PAH emitters are larger and more ionized. The warm molecular gas to PAH flux ratios (H$_2/PAH$) are enhanced in the outflow by factors of 10-100 as compared to the starburst disk. This enhancement in the H$_2/PAH$ ratio does not seem to follow the ionization of the atomic gas (as measured with the [NeIII]/[NeII] line flux ratio) in the outflow. This suggests that much of the warm H$_2$ in the outflow is excited by shocks. The observed H$_2$ line intensities can be reproduced with low velocity shocks ($v < 40$ km s$^{-1}$) driven into moderately dense molecular gas ($10^2 <n_H < 10^4$ cm$^{-3}$) entrained in the outflow.

Protoplanetary disks are dispersed by viscous evolution and photoevaporation in a few million years; in the interim small, sub-micron sized dust grains must grow and form planets. The time-varying abundance of small grains in an evolving disk directly affects gas heating by far-ultraviolet photons, while dust evolution affects photoevaporation by changing the disk opacity and resulting penetration of FUV photons in the disk. Photoevaporative flows, in turn, selectively carry small dust grains leaving the larger particles---which decouple from the gas---behind in the disk. We study these effects by investigating the evolution of a disk subject to viscosity, photoevaporation by EUV, FUV and X-rays, dust evolution, and radial drift using a 1-D multi-fluid approach (gas + different dust grain sizes) to solve for the evolving surface density distributions. The 1-D evolution is augmented by 1+1D models constructed at each epoch to obtain the instantaneous disk structure and determine photoevaporation rates. The implementation of a dust coagulation/fragmentation model results in a marginal decrease in disk lifetimes when compared to models with no dust evolution; the disk lifetime is thus found to be relatively insensitive to the evolving dust opacity. We find that photoevaporation can cause significant reductions in the gas/dust mass ratio in the planet-forming regions of the disk as it evolves, and may result in a corresponding increase in heavy element abundances relative to hydrogen. We discuss implications for theories of planetesimal formation and giant planet formation, including the formation of gas-poor giants. After gas disk dispersal, $\sim 3\times 10^{-4}$ \ms\ of mass in solids typically remain, comparable to the solids inventory of our solar system.

The amount of high-energy stellar radiation reaching the surface of protoplanetary disks is essential to determine their chemistry and physical evolution. Here, we use millimetric and centimetric radio data to constrain the EUV luminosity impinging on 14 disks around young (~2-10Myr) sun-like stars. For each object we identify the long-wavelength emission in excess to the dust thermal emission, attribute that to free-free disk emission, and thereby compute an upper limit to the EUV reaching the disk. We find upper limits lower than 10$^{42}$ photons/s for all sources without jets and lower than $5 \times 10^{40}$ photons/s for the three older sources in our sample. These latter values are low for EUV-driven photoevaporation alone to clear out protoplanetary material in the timescale inferred by observations. In addition, our EUV upper limits are too low to reproduce the [NeII] 12.81 micron luminosities from three disks with slow [NeII]-detected winds. This indicates that the [NeII] line in these sources primarily traces a mostly neutral wind where Ne is ionized by 1 keV X-ray photons, implying higher photoevaporative mass loss rates than those predicted by EUV-driven models alone. In summary, our results suggest that high-energy stellar photons other than EUV may dominate the dispersal of protoplanetary disks around sun-like stars.

We report observations of molecular oxygen (O$_2$) rotational transitions at 487 GHz, 774 GHz, and 1121 GHz toward Orion Peak A. The O2 lines at 487 GHz and 774 GHz are detected at velocities of 10-12 km/s with line widths 3 km/s; however, the transition at 1121 GHz is not detected. The observed line characteristics, combined with the results of earlier observations, suggest that the region responsible for the O$_2$ emission is 9" (6e16 cm) in size, and is located close to the H2 Peak 1position (where vibrationally-excited H$_2$ emission peaks), and not at Peak A, 23" away. The peak O2 column density is 1.1e18/cm2. The line velocity is close to that of 621 GHz water maser emission found in this portion of the Orion Molecular Cloud, and having a shock with velocity vector lying nearly in the plane of the sky is consistent with producing maximum maser gain along the line-of-sight. The enhanced O$_2$ abundance compared to that generally found in dense interstellar clouds can be explained by passage of a low-velocity C-shock through a clump with preshock density 2e4/cm3, if a reasonable flux of UV radiation is present. The postshock O$_2$ can explain the emission from the source if its line of sight dimension is ~10 times larger than its size on the plane of the sky. The special geometry and conditions required may explain why O$_2$ emission has not been detected in the cores of other massive star-forming molecular clouds.

We present spectral line images of [CI] 809 GHz, CO J=1-0 115 GHz and HI 1.4 GHz line emission, and calculate the corresponding C, CO and H column densities, for a sinuous, quiescent Giant Molecular Cloud about 5 kpc distant along the l=328{\deg} sightline (hereafter G328) in our Galaxy. The [CI] data comes from the High Elevation Antarctic Terahertz (HEAT) telescope, a new facility on the summit of the Antarctic plateau where the precipitable water vapor falls to the lowest values found on the surface of the Earth. The CO and HI datasets come from the Mopra and Parkes/ATCA telescopes, respectively. We identify a filamentary molecular cloud, ~75 x 5 pc long with mass ~4 x 10E4 Msun and a narrow velocity emission range of just 4 km/s. The morphology and kinematics of this filament are similar in CO, [CI] and HI, though in the latter appears as self-absorption. We calculate line fluxes and column densities for the three emitting species, which are broadly consistent with a PDR model for a GMC exposed to the average interstellar radiation field. The [C/CO] abundance ratio averaged through the filament is found to be approximately unity. The G328 filament is constrained to be cold (Tdust < 20K) by the lack of far-IR emission, to show no clear signs of star formation, and to only be mildly turbulent from the narrow line width. We suggest that it may represent a GMC shortly after formation, or perhaps still be in the process of formation.

We present the first results from a new carbon monoxide (CO) survey of the southern Galactic plane being conducted with the Mopra radio telescope in Australia. The 12CO, 13CO and C18O J=1-0 lines are being mapped over the l = 305-345 deg, b = +/- 0.5 deg portion of the 4th quadrant of the Galaxy, at 35" spatial and 0.1 km/s spectral resolution. The survey is being undertaken with two principal science objectives: (i) to determine where and how molecular clouds are forming in the Galaxy and (ii) to probe the connection between molecular clouds and the "missing" gas inferred from gamma-ray observations. We describe the motivation for the survey, the instrumentation and observing techniques being applied, and the data reduction and analysis methodology. In this paper we present the data from the first degree surveyed, l = 323-324 deg, b = +/- 0.5 deg. We compare the data to the previous CO survey of this region and present metrics quantifying the performance being achieved; the rms sensitivity per 0.1 km/s velocity channel is ~1.5K for 12CO and ~0.7K for the other lines. We also present some results from the region surveyed, including line fluxes, column densities, molecular masses, 12CO/13CO line ratios and 12CO optical depths. We also examine how these quantities vary as a function of distance from the Sun when averaged over the 1 square degree survey area. Approximately 2 x 10E6 MSun of molecular gas is found along the G323 sightline, with an average H2 number density of nH2 ~ 1 cm-3 within the Solar circle. The CO data cubes will be made publicly available as they are published.

We present a model in which the 22 GHz H$_2$O masers observed in star-forming regions occur behind shocks propagating in dense regions (preshock density $n_0 \sim 10^6 - 10^8$ cm$^{-3}$). We focus on high-velocity ($v_s > 30$ km s$^{-1}$) dissociative J shocks in which the heat of H$_2$ re-formation maintains a large column of $\sim 300-400$ K gas; at these temperatures the chemistry drives a considerable fraction of the oxygen not in CO to form H$_2$O. The H$_2$O column densities, the hydrogen densities, and the warm temperatures produced by these shocks are sufficiently high to enable powerful maser action. The observed brightness temperatures (generally $\sim 10^{11} - 10^{14}$ K) are the result of coherent velocity regions that have dimensions in the shock plane that are 10 to 100 times the shock thickness of $\sim 10^{13}$ cm. The masers are therefore beamed towards the observer, who typically views the shock "edge-on", or perpendicular to the shock velocity; the brightest masers are then observed with the lowest line of sight velocities with respect to the ambient gas. We present numerical and analytic studies of the dependence of the maser inversion, the resultant brightness temperature, the maser spot size and shape, the isotropic luminosity, and the maser region magnetic field on the shock parameters and the coherence path length; the overall result is that in galactic H$_2$O 22 GHz masers these observed parameters can be produced in J shocks with $n_0\sim 10^6 - 10^8$ cm$^{-3}$ and $v_s \sim 30 -200$ km s$^{-1}$. A number of key observables such as maser shape, brightness temperature, and global isotropic luminosity depend only on the particle flux into the shock, $j=n_0v_s$, rather than on $n_0$ and $v_s$ separately.

The formation time, masses, and location of planets are strongly impacted by the physical mechanisms that disperse protoplanetary disks and the timescale over which protoplanetary material is cleared out. Accretion of matter onto the central star, protostellar winds/jets, magnetic disk winds, and photoevaporative winds operate concurrently. Hence, disentangling their relative contribution to disk dispersal requires identifying diagnostics that trace different star-disk environments. Here, we analyze the low velocity component (LVC) of the Oxygen optical forbidden lines, which is found to be blueshifted by a few km/s with respect to the stellar velocity. We find that the [OI] LVC profiles are different from those of [NeII] at 12.81\mu m and CO at 4.7\mu m lines pointing to different origins for these gas lines. We report a correlation between the luminosity of the [OI] LVC and the accretion luminosity Lacc. We do not find any correlation with the X-ray luminosity, while we find that the higher is the stellar FUV luminosity, the higher is the luminosity of the [OI] LVC. In addition, we show that the [OI]6300\AA/5577\AA\ ratio is low (ranging between 1 and 8). These findings favor an origin of the [OI] LVC in a region where OH is photodissociated by stellar FUV photons and argue against thermal emission from an X-ray-heated layer. Detailed modeling of two spectra with the highest S/N and resolution shows that there are two components within the LVC: a broad, centrally peaked component that can be attributed to gas arising in a warm disk surface in Keplerian rotation (with FWHM between ~40 and ~60 km/s), and a narrow component (with FWHM ~10 km/s and small blueshifts of ~2 km/s) that may arise in a cool (<1,000 K) molecular wind.

We present Herschel PACS spectroscopy of the [OI] 63 micron gas-line for three circumstellar disk systems showing signs of significant disk evolution and/or planet formation: HR 8799, HD 377 and RX J1852.3-3700. [OI] is undetected toward HR 8799 and HD 377 with 3 sigma upper limits of 6.8 x 10^-18 W m^-2 and 9.9 x 10^-18 W m^-2 respectively. We find an [OI] detection for RX J1852.3-3700 at 12.3 +- 1.8 x 10^-18 W m^-2. We use thermo-chemical disk models to model the gas emission, using constraints on the [OI] 63 micron, and ancillary data to derive gas mass upper limits and constrain gas-to-dust ratios. For HD 377 and HR 8799, we find 3 sigma upper limits on the gas mass of 0.1-20 Mearth. For RX J1852.3-3700, we find two distinct disk scenarios that could explain the detection of [OI] 63 micron and CO(2-1) upper limits reported from the literature: (i) a large disk with gas co-located with the dust (16-500 AU), resulting in a large tenuous disk with ~16 Mearth of gas, or (ii) an optically thick gas disk, truncated at ~70 AU, with a gas mass of 150 Mearth. We discuss the implications of these results for the formation and evolution of planets in these three systems.

We model the production of OH+, H2O+, and H3O+ in interstellar clouds, using a steady state photodissociation region code that treats the freeze-out of gas species, grain surface chemistry, and desorption of ices from grains. The code includes PAHs, which have important effects on the chemistry. All three ions generally have two peaks in abundance as a function of depth into the cloud, one at A_V<~1 and one at A_V~3-8, the exact values depending on the ratio of incident ultraviolet flux to gas density. For relatively low values of the incident far ultraviolet flux on the cloud ({\chi}<~ 1000; {\chi}= 1= local interstellar value), the columns of OH+ and H2O+ scale roughly as the cosmic ray primary ionization rate {\zeta}(crp) divided by the hydrogen nucleus density n. The H3O+ column is dominated by the second peak, and we show that if PAHs are present, N(H3O+) ~ 4x10^{13} cm^{-2} independent of {\zeta}(crp) or n. If there are no PAHs or very small grains at the second peak, N(H3O+) can attain such columns only if low ionization potential metals are heavily depleted. We also model diffuse and translucent clouds in the interstellar medium, and show how observations of N(OH+)/N(H) and N(OH+)/N(H2O+) can be used to estimate {\zeta}(crp)/n, {\chi}/n and A_V in them. We compare our models to Herschel observations of these two ions, and estimate {\zeta}(crp) ~ 4-6 x 10^-16 (n/100 cm^-3) s^-1 and \chi/n = 0.03 cm^3 for diffuse foreground clouds towards W49N.

Recent infrared observations have demonstrated that photoevaporation driven by high-energy photons from the central star contributes to the dispersal of protoplanetary disks. Here, we show that photoevaporative winds should produce a detectable free-free continuum emission given the range of stellar ionizing photons and X-ray luminosities inferred for young sun-like stars. We point out that VLA observations of the nearby disk around TWHya might have already detected this emission at centimeter wavelengths and calculate the wind electron density and mass flow rate. We also estimate the intensities of H radio recombination lines tracing the wind and discuss which ones could be detected with current instrumentation. The detection and profiles of these recombination lines would unambiguously prove our inference of free-free emission from photoevaporating disks like TWHya. In addition, radio/millimeter data can help constraining wind parameters such as temperature and electron density that are fundamental in measuring mass flow rates.

We report the results of a search for molecular oxygen (O2) toward the Orion Bar, a prominent photodissociation region at the southern edge of the HII region created by the luminous Trapezium stars. We observed the spectral region around the frequency of the O2 N_J = 3_3 - 1_2 transition at 487 GHz and the 5_4 - 3_4 transition at 774 GHz using the Heterodyne Instrument for the Far Infrared on the Herschel Space Observatory. Neither line was detected, but the 3sigma upper limits established here translate to a total line-of-sight O2 column density < 1.5 10^16 cm^-2 for an emitting region whose temperature is between 30K and 250 K, or < 1 10^16 cm^-2 if the O2 emitting region is primarily at a temperature of ~< 100 K. Because the Orion Bar is oriented nearly edge-on relative to our line of sight, the observed column density is enhanced by a factor estimated to be between 4 and 20 relative to the face-on value. Our upper limits imply that the face-on O2 column density is less than 4 10^15 cm^-2, a value that is below, and possibly well below, model predictions for gas with a density of 10^4 - 10^5 cm^-3 exposed to a far ultraviolet flux 10^4 times the local value, conditions inferred from previous observations of the Orion Bar. The discrepancy might be resolved if: (1) the adsorption energy of O atoms to ice is greater than 800 K; (2) the total face-on Av of the Bar is less than required for O2 to reach peak abundance; (3) the O2 emission arises within dense clumps with a small beam filling factor; or, (4) the face-on depth into the Bar where O2 reaches its peak abundance, which is density dependent, corresponds to a sky position different from that sampled by our Herschel beams.

Models of pure gas-phase chemistry in well-shielded regions of molecular clouds predict relatively high levels of molecular oxygen, O2, and water, H2O. Contrary to expectation, the space missions SWAS and Odin found only very small amounts of water vapour and essentially no O2 in the dense star-forming interstellar medium. Only toward rho Oph A did Odin detect a weak line of O2 at 119 GHz in a beam size of 10 arcmin. A larger telescope aperture such as that of the Herschel Space Observatory is required to resolve the O2 emission and to pinpoint its origin. We use the Heterodyne Instrument for the Far Infrared aboard Herschel to obtain high resolution O2 spectra toward selected positions in rho Oph A. These data are analysed using standard techniques for O2 excitation and compared to recent PDR-like chemical cloud models. The 487.2GHz line was clearly detected toward all three observed positions in rho Oph A. In addition, an oversampled map of the 773.8GHz transition revealed the detection of the line in only half of the observed area. Based on their ratios, the temperature of the O2 emitting gas appears to vary quite substantially, with warm gas (> 50 K) adjacent to a much colder region, where temperatures are below 30 K. The exploited models predict O2 column densities to be sensitive to the prevailing dust temperatures, but rather insensitive to the temperatures of the gas. In agreement with these model, the observationally determined O2 column densities seem not to depend strongly on the derived gas temperatures, but fall into the range N(O2) = (3 to >6)e15/cm^2. Beam averaged O2 abundances are about 5e-8 relative to H2. Combining the HIFI data with earlier Odin observations yields a source size at 119 GHz of about 4 - 5 arcmin, encompassing the entire rho Oph A core.

We use the IRS on Spitzer to observe the southern part of the reflection nebula NGC 2023, including the Southern Ridge, which is a photodissociation region par excellence excited by HD 37903. Five pure-rotational H2 emission lines are detected and mapped over and around the Southern Ridge in order to compare with predicted level column densities from theoretical PDR models. We find very good agreement between PDR model predictions and emission line intensities and ratios measured with Spitzer, leading us to conclude that grain photoelectric heating sufficiently warms the gas to produce the observed H2 line emission via collisional excitation. On the Southern Ridge, we infer a hydrogen nucleus density n_H approx 2x10^5 cm^-3 and radiation field strength \chi\ approx 10^4 relative to the local Galactic interstellar radiation field. This high value for \chi\ independently predicts a distance toward HD 37903 of 300 pc, and is consistent with the most recent HIPPARCOS results. Over the map we find that both n_H and \chi\ vary by a factor of sim 3. Such 2-D variations provide clues about the underlying 3-D structure of the Southern Ridge field, which appears to be the tip of a molecular cloud. We also map variations in excitation temperature and the ortho-to-para ratio, the latter attaining values of sim 1.5---2.0 on the Southern Ridge, and find that PDR modeling can readily reproduce observed ortho-to-para ratios that are < 3 for rotational excitation dominated by collisional processes. Last, the stars Sellgren C and G are discovered to be resolved on archival HST images into two point sources each, with separations of lesssim 0.5".

We report observations of three rotational transitions of molecular oxygen (O2) in emission from the H2 Peak 1 position of vibrationally excited molecular hydrogen in Orion. We observed the 487 GHz, 774 GHz, and 1121 GHz lines using HIFI on the Herschel Space Observatory, having velocities of 11 km s-1 to 12 km s-1 and widths of 3 km s-1. The beam-averaged column density is N(O2) = 6.5\times1016 cm-2, and assuming that the source has an equal beam filling factor for all transitions (beam widths 44, 28, and 19"), the relative line intensities imply a kinetic temperature between 65 K and 120 K. The fractional abundance of O2 relative to H2 is 0.3 - 7.3\times10-6. The unusual velocity suggests an association with a ~ 5" diameter source, denoted Peak A, the Western Clump, or MF4. The mass of this source is ~ 10 M\odot and the dust temperature is \geq 150 K. Our preferred explanation of the enhanced O2 abundance is that dust grains in this region are sufficiently warm (T \geq 100 K) to desorb water ice and thus keep a significant fraction of elemental oxygen in the gas phase, with a significant fraction as O2. For this small source, the line ratios require a temperature \geq 180 K. The inferred O2 column density \simeq 5\times1018 cm-2 can be produced in Peak A, having N(H2) \simeq 4\times1024 cm-2. An alternative mechanism is a low-velocity (10 to 15 km s-1) C-shock, which can produce N(O2) up to 1017 cm-2.

Photoevaporation driven by the central star is expected to be a ubiquitous and important mechanism to disperse the circumstellar dust and gas from which planets form. Here, we present a detailed study of the circumstellar disk surrounding the nearby star TW Hya and provide observational constraints to its photoevaporative wind. Our new high-resolution (R ~ 30,000) mid-infrared spectroscopy in the [Ne II] 12.81 {\mu}m line confirms that this gas diagnostic traces the unbound wind component within 10AU from the star. From the blueshift and asymmetry in the line profile, we estimate that most (>80%) of the [Ne II] emission arises from disk radii where the midplane is optically thick to the redshifted outflowing gas, meaning beyond the 1 or 4AU dust rim inferred from other observations. We re-analyze high-resolution (R ~ 48, 000) archival optical spectra searching for additional transitions that may trace the photoevaporative flow. Unlike the [Ne II] line, optical forbidden lines from OI, SII, and MgI are centered at the stellar velocity and have symmetric profiles. The only way these lines could trace the photoevaporative flow is if they arise from a disk region physically distinct from that traced by the [Ne II] line, specifically from within the optically thin dust gap. However, the small (~10 km/s) FWHM of these lines suggest that most of the emitting gas traced at optical wavelengths is bound to the system rather than unbound. We discuss the implications of our results for a planet-induced versus a photoevaporation-induced gap.

We compare line emission calculated from theoretical disk models with optical to sub-millimeter wavelength observational data of the gas disk surrounding TW Hya and infer the spatial distribution of mass in the gas disk. The model disk that best matches observations has a gas mass ranging from $\sim10^{-4}-10^{-5}$\ms\ for $0.06{\rm AU} <r<3.5$AU and $\sim 0.06$\ms\ for $ 3.5 {\rm AU} <r<200$AU. We find that the inner dust hole ($r<3.5$AU) in the disk must be depleted of gas by $\sim 1-2$ orders of magnitude compared to the extrapolated surface density distribution of the outer disk. Grain growth alone is therefore not a viable explanation for the dust hole. CO vibrational emission arises within $r\sim 0.5$AU from thermal excitation of gas. [OI] 6300\AA\ and 5577\AA\ forbidden lines and OH mid-infrared emission are mainly due to prompt emission following UV photodissociation of OH and water at $r\lesssim0.1$AU and at $r\sim 4$AU. [NeII] emission is consistent with an origin in X-ray heated neutral gas at $r\lesssim 10$AU, and may not require the presence of a significant EUV ($h\nu>13.6$eV) flux from TW Hya. H$_2$ pure rotational line emission comes primarily from $r\sim 1-30$AU. [OI]63$\mu$m, HCO$^+$ and CO pure rotational lines all arise from the outer disk at $r\sim30-120$AU. We discuss planet formation and photoevaporation as causes for the decrease in surface density of gas and dust inside 4 AU. If a planet is present, our results suggest a planet mass $\sim 4-7$M$_J$ situated at $\sim 3$AU. Using our photoevaporation models and the best surface density profile match to observations, we estimate a current photoevaporative mass loss rate of $4\times10^{-9}$\ms\ yr$^{-1}$ and a remaining disk lifetime of $\sim 5$ million years.

We report the results of a large-area study of water vapor along the Orion Molecular Cloud ridge, the purpose of which was to determine the depth-dependent distribution of gas-phase water in dense molecular clouds. We find that the water vapor measured toward 77 spatial positions along the face-on Orion ridge, excluding positions surrounding the outflow associated with BN/KL and IRc2, display integrated intensities that correlate strongly with known cloud surface tracers such as CN, C2H, 13CO J =5-4, and HCN, and less well with the volume tracer N2H+. Moreover, at total column densities corresponding to Av < 15 mag., the ratio of H2O to C18O integrated intensities shows a clear rise approaching the cloud surface. We show that this behavior cannot be accounted for by either optical depth or excitation effects, but suggests that gas-phase water abundances fall at large Av. These results are important as they affect measures of the true water-vapor abundance in molecular clouds by highlighting the limitations of comparing measured water vapor column densities with such traditional cloud tracers as 13CO or C18O. These results also support cloud models that incorporate freeze-out of molecules as a critical component in determining the depth-dependent abundance of water vapor.

We derive a semi-empirical galactic initial mass function (IMF) from observational constraints. We assume that the star formation rate in a galaxy can be expressed as the product of the IMF, $\psi (m)$, which is a smooth function of mass $m$ (in units of \msun), and a time- and space-dependent total rate of star formation per unit area of galactic disk. The mass dependence of the proposed IMF is determined by five parameters: the low-mass slope $\gamma$, the high-mass slope $-\Gamma$, the characteristic mass $m_{ch}$ (which is close to the mass $m_{\rm peak}$ at which the IMF turns over), and the lower and upper limits on the mass, $m_l$ (taken to be 0.004) and $m_u$ (taken to be 120). The star formation rate in terms of number of stars per unit area of galactic disk per unit logarithmic mass interval, is proportional to $m^{-\Gamma} \left\{1-\exp\left[{-(m/m_{ch})^{\gamma +\Gamma}}\right]\right\}$, where $\cal N_*$ is the number of stars, $m_l<m<m_u$ is the range of stellar masses. The values of $\gamma$ and $\emch$ are derived from two integral constraints: i) the ratio of the number density of stars in the range $m=0.1-0.6$ to that in the range $m=0.6-0.8$ as inferred from the mass distribution of field stars in the local neighborhood, and ii) the ratio of the number of stars in the range $m=0.08 - 1$ to the number of brown dwarfs in the range $m=0.03-0.08$ in young clusters. The IMF satisfying the above constraints is characterized by the parameters $\gamma=0.51$ and $\emch=0.35$ (which corresponds to $m_{\rm peak}=0.27$). This IMF agrees quite well with the Chabrier (2005) IMF for the entire mass range over which we have compared with data, but predicts significantly more stars with masses $< 0.03\, M_\odot$; we also compare with other IMFs in current use.

We report the first detection of chloronium, H$_2$Cl$^+$, in the interstellar medium, using the HIFI instrument aboard the \emph{Herschel} Space Observatory. The $2_{12}-1_{01}$ lines of ortho-H$_2^{35}$Cl$^+$ and ortho-H$_2^{37}$Cl$^+$ are detected in absorption towards NGC~6334I, and the $1_{11}-0_{00}$ transition of para-H$_2^{35}$Cl$^+$ is detected in absorption towards NGC~6334I and Sgr~B2(S). The H$_2$Cl$^+$ column densities are compared to those of the chemically-related species HCl. The derived HCl/H$_2$Cl$^+$ column density ratios, $\sim$1--10, are within the range predicted by models of diffuse and dense Photon Dominated Regions (PDRs). However, the observed H$_2$Cl$^+$ column densities, in excess of $10^{13}$~cm$^{-2}$, are significantly higher than the model predictions. Our observations demonstrate the outstanding spectroscopic capabilities of HIFI for detecting new interstellar molecules and providing key constraints for astrochemical models.

The mass of molecular gas in an interstellar cloud is often measured using line emission from low rotational levels of CO, which are sensitive to the CO mass, and then scaling to the assumed molecular hydrogen H_2 mass. However, a significant H_2 mass may lie outside the CO region, in the outer regions of the molecular cloud where the gas phase carbon resides in C or C+. Here, H_2 self-shields or is shielded by dust from UV photodissociation, where as CO is photodissociated. This H_2 gas is "dark" in molecular transitions because of the absence of CO and other trace molecules, and because H_2 emits so weakly at temperatures 10 K < T < 100 K typical of this molecular component. This component has been indirectly observed through other tracers of mass such as gamma rays produced in cosmic ray collisions with the gas and far-infrared/submillimeter wavelength dust continuum radiation. In this paper we theoretically model this dark mass and find that the fraction of the molecular mass in this dark component is remarkably constant (~ 0.3 for average visual extinction through the cloud with mean A_V ~ 8) and insensitive to the incident ultraviolet radiation field strength, the internal density distribution, and the mass of the molecular cloud as long as mean A_V, or equivalently, the product of the average hydrogen nucleus column and the metallicity through the cloud, is constant. We also find that the dark mass fraction increases with decreasing mean A_V, since relatively more molecular H_2 material lies outside the CO region in this case.

We report sensitive Spitzer IRS spectroscopy in the 10-20 micron region of TW Hya, a nearby T Tauri star. The unusual spectral energy distribution of the source, that of a transition object, indicates that the circumstellar disk in the system has experienced significant evolution, possibly as a result of planet formation. The spectrum we measure is strikingly different from that of other classical T Tauri stars reported in the literature, displaying no strong emission features of water, acetylene, or HCN. The difference indicates that the inner planet formation region (within 5 AU) of the gaseous disk has evolved physically and/or chemically away from the classical T Tauri norm. Nevertheless, TW Hya does show a rich spectrum of emission features of atoms (HI, [NeII], and [NeIII]) and molecules (H2, OH, CO2, HCO+, and possibly CH3), some of which are also detected in classical T Tauri spectra. The properties of the neon emission are consistent with an origin for the emission in a disk irradiated by X-rays (with a possible role for additional irradiation by stellar EUV). The OH emission we detect, which also likely originates in the disk, is hot, arising from energy levels up to 23,000 K above ground, and may be produced by the UV photodissociation of water. The HI emission is surprisingly strong, with relative strengths that are consistent with case B recombination. While the absence of strong molecular emission in the 10-20 micron region may indicate that the inner region of the gaseous disk has been partly cleared by an orbiting giant planet, chemical and/or excitation effects may be responsible instead. We discuss these issues and how our results bear on our understanding of the evolutionary state of the TW Hya disk.

We present the time evolution of viscously accreting circumstellar disks as they are irradiated by ultraviolet and X-ray photons from a low-mass central star. Our model is a hybrid of a 1D time-dependent viscous disk model coupled to a 1+1D disk vertical structure model used for calculating the disk structure and photoevaporation rates. We find that disks of initial mass 0.1M_o around 1M_o stars survive for 4x10^6 years, assuming a viscosity parameter $\alpha=0.01$, a time-dependent FUV luminosity $L_{FUV}~10^{-2}-10^{-3}$ L_o and with X-ray and EUV luminosities $L_X \sim L_{EUV} ~ 10^{-3}$L_o. We find that FUV/X-ray-induced photoevaporation and viscous accretion are both important in depleting disk mass. Photoevaporation rates are most significant at ~ 1-10 AU and at >~ 30 AU. Viscosity spreads the disk which causes mass loss by accretion onto the central star and feeds mass loss by photoevaporation in the outer disk. We find that FUV photons can create gaps in the inner, planet-forming regions of the disk (~ 1-10 AU) at relatively early epochs in disk evolution while disk masses are still substantial. EUV and X-ray photons are also capable of driving gaps, but EUV can only do so at late, low accretion-rate epochs after the disk mass has already declined substantially. Disks around stars with predominantly soft X-ray fields experience enhanced photoevaporative mass loss. We follow disk evolution around stars of different masses, and find that disk survival time is relatively independent of mass for stars with M <~ 3M_o; for M >~ 3M_o the disks are short-lived(~10^5 years).

Extreme ultraviolet (EUV, 13.6 eV $< h\nu \lta 100$ eV) and X-rays in the 0.1-2 keV band can heat the surfaces of disks around young, low mass stars to thousands of degrees and ionize species with ionization potentials greater than 13.6 eV. Shocks generated by protostellar winds can also heat and ionize the same species close to the star/disk system. These processes produce diagnostic lines (e.g., [NeII] 12.8 $\mu$m and [OI] 6300 \AA) that we model as functions of key parameters such as EUV luminosity and spectral shape, X-ray luminosity and spectral shape, and wind mass loss rate and shock speed. Comparing our models with observations, we conclude that either internal shocks in the winds or X-rays incident on the disk surfaces often produce the observed [NeII] line, although there are cases where EUV may dominate. Shocks created by the oblique interaction of winds with disks are unlikely [NeII] sources because these shocks are too weak to ionize Ne. Even if [NeII] is mainly produced by X-rays or internal wind shocks, the neon observations typically place upper limits of $\lta 10^{42}$ s$^{-1}$ on the EUV photon luminosity of these young low mass stars. The observed [OI] 6300 \AA line has both a low velocity component (LVC) and a high velocity component. The latter likely arises in internal wind shocks. For the former we find that X-rays likely produce more [OI] luminosity than either the EUV layer, the transition layer between the EUV and X-ray layer, or the shear layer where the protostellar wind shocks and entrains disk material in a radial flow across the surface of the disk. Our soft X-ray models produce [OI] LVCs with luminosities up to $10^{-4}$ L$_\odot$, but may not be able to explain the most luminous LVCs.

We present a model for the dispersal of protoplanetary disks by winds from either the central star or the inner disk. These winds obliquely strike the flaring disk surface and strip away disk material by entraining it in an outward radial-moving flow at the wind-disk interface which lies several disk scale heights above the mid-plane. The disk dispersal time depends on the entrainment velocity at which disk material flows into this turbulent shear layer interface. If the entrainment efficiency is ~10% of the local sound speed, a likely upper limit, the dispersal time at 1 AU is ~6 Myr for a disk with a surface density of 10^3 g cm^{-2}, a solar mass central star, and a wind with an outflow rate 10^{-8} Msun/yr and terminal velocity 200 km/s. When compared to photoevaporation and viscous evolution, wind stripping can be a dominant mechanism only for the combination of low accretion rates (< 10^{-8} Msun/yr) and wind outflow rates approaching these accretion rates. This case is unusual since generally outflow rates are < 0.1 of of accretion rates.

This white paper is submitted to the Astronomy and Astrophysics 2010 Decadal Survey (Astro2010)1 Committee on the State of the Profession to emphasize the potential of the Stratospheric Observatory for Infrared Astronomy (SOFIA) to contribute to the training of instrumentalists and observers, and to related technology developments. This potential goes beyond the primary mission of SOFIA, which is to carry out unique, high priority astronomical research. SOFIA is a Boeing 747SP aircraft with a 2.5 meter telescope. It will enable astronomical observations anywhere, any time, and at most wavelengths between 0.3 microns and 1.6 mm not accessible from ground-based observatories. These attributes, accruing from the mobility and flight altitude of SOFIA, guarantee a wealth of scientific return. Its instrument teams (nine in the first generation) and guest investigators will do suborbital astronomy in a shirt-sleeve environment. The project will invest $10M per year in science instrument development over a lifetime of 20 years. This, frequent flight opportunities, and operation that enables rapid changes of science instruments and hands-on in-flight access to the instruments, assure a unique and extensive potential - both for training young instrumentalists and for encouraging and deploying nascent technologies. Novel instruments covering optical, infrared, and submillimeter bands can be developed for and tested on SOFIA by their developers (including apprentices) for their own observations and for those of guest observers, to validate technologies and maximize observational effectiveness.

We model the temperature and chemical structure of molecular clouds as a function of depth into the cloud, assuming a cloud of constant density n illuminated by an external FUV (6 eV < E < 13.6 eV) flux G_0 (scaling factor in multiples of the local interstellar field). Extending previous photodissociation region models, we include the freezing of species, simple grain surface chemistry, and desorption (including FUV photodesorption) of ices. We also treat the opaque cloud interior with time-dependent chemistry. Here, under certain conditions, gas phase elemental oxygen freezes out as water ice and the elemental C/O abundance ratio can exceed unity, leading to complex carbon chemistry. Gas phase H2O and O2 peak in abundance at intermediate depth into the cloud, roughly A_V~3-8 from the surface, the depth proportional to ln(G_0/n). Closer to the surface, molecules are photodissociated. Deeper into the cloud, molecules freeze to grain surfaces. At intermediate depths photodissociation rates are attenuated by dust extinction, but photodesorption prevents total freezeout. For G_0 < 500, abundances of H2O and O2 peak at values ~10^(-7), producing columns ~10^(15) per cm^2, independent of G_0 and n. The peak abundances depend primarily on the product of the photodesorption yield of water ice and the grain surface area per H nucleus. At higher values of G_0, thermal desorption of O atoms from grains enhances the gas phase H2O peak abundance and column slightly, whereas the gas phase O2 peak abundance rises to ~10^(-5) and the column to ~2x10^(16) per cm^2. We present simple analytic equations for the abundances as a function of depth which clarify the dependence on parameters. The models are applied to observations of H2O, O2, and water ice in a number of sources, including B68, NGC 2024, and Rho Oph.

High resolution mid-infrared spectra are presented for 155 nuclear and extranuclear regions from the Spitzer Infrared Nearby Galaxies Survey (SINGS). The fluxes for nine atomic forbidden and three molecular hydrogen mid-infrared emission lines are also provided, along with upper limits in key lines for infrared-faint targets. The SINGS sample shows a wide range in the ratio of [SIII]18.71um/[SIII]33.48um, but the average ratio of the ensemble indicates a typical interstellar electron density of 300-400 cm^{-3} on ~23"x15" scales and 500-600 cm^{-3} using ~11"x9" apertures, independent of whether the region probed is a star-forming nuclear, a star-forming extranuclear, or an AGN environment. Evidence is provided that variations in gas-phase metallicity play an important role in driving variations in radiation field hardness, as indicated by [NeIII]15.56um/[NeII]12.81um, for regions powered by star formation. Conversely, the radiation hardness for galaxy nuclei powered by accretion around a massive black hole is independent of metal abundance. Furthermore, for metal-rich environments AGN are distinguishable from star-forming regions by significantly larger [NeIII]15.56um/[NeII]12.81um ratios. Finally, [FeII]25.99um/[NeII]12.81um versus [SiII]34.82um/[SIII]33.48um also provides an empirical method for discerning AGN from normal star-forming sources. However, similar to [NeIII]15.56um/[NeII]12.81um, these mid-infrared line ratios lose their AGN/star-formation diagnostic powers for very low metallicity star-forming systems with hard radiation fields.

We present Spitzer photometric (IRAC and MIPS) and spectroscopic (IRS low resolution) observations for 314 stars in the Formation and Evolution of Planetary Systems (FEPS) Legacy program. These data are used to investigate the properties and evolution of circumstellar dust around solar-type stars spanning ages from approximately 3 Myr to 3 Gyr. We identify 46 sources that exhibit excess infrared emission above the stellar photosphere at 24um, and 21 sources with excesses at 70um. Five sources with an infrared excess have characteristics of optically thick primordial disks, while the remaining sources have properties akin to debris systems. The fraction of systems exhibiting a 24um excess greater than 10.2% above the photosphere is 15% for ages < 300 Myr and declines to 2.7% for older ages. The upper envelope to the 70um fractional luminosity appears to decline over a similar age range. The characteristic temperature of the debris inferred from the IRS spectra range between 60 and 180 K, with evidence for the presence of cooler dust to account for the strength of the 70um excess emission. No strong correlation is found between dust temperature and stellar age. Comparison of the observational data with disk models containing a power-law distribution of silicate grains suggest that the typical inner disk radius is > 10 AU. Although the interpretation is not unique, the lack of excess emission shortwards of 16um and the relatively flat distribution of the 24um excess for ages <300~Myr is consistent with steady-state collisional models.

We calculate the rate of photoevaporation of a circumstellar disk by energetic radiation (FUV, 6eV $<h\nu<$13.6eV; EUV, 13.6eV $<h\nu<$0.1keV; and Xrays, $h\nu>0.1$keV) from its central star. We focus on the effects of FUV and X-ray photons since EUV photoevaporation has been treated previously, and consider central star masses in the range $0.3-7 {\rm M}_{\odot}$. Contrary to the EUV photoevaporation scenario, which creates a gap at about $r_g\sim 7\ (M_*/1{\rm M}_{\odot})$ AU and then erodes the outer disk from inside out, we find that FUV photoevaporation predominantly removes less bound gas from the outer disk. Heating by FUV photons can cause significant erosion of the outer disk where most of the mass is typically located. X-rays indirectly increase the mass loss rates (by a factor $\sim 2$) by ionizing the gas, thereby reducing the positive charge on grains and PAHs and enhancing FUV-induced grain photoelectric heating. FUV and X-ray photons may create a gap in the disk at $\sim 10$ AU under favourable circumstances. Photoevaporation timescales for M$_* \sim 1{\rm M}_{\odot}$ stars are estimated to be $\sim 10^6$ years, after the onset of disk irradiation by FUV and X-rays. Disk lifetimes do not vary much for stellar masses in the range $0.3-3$M$_{\odot}$. More massive stars ($\gtrsim 7 {\rm M}_{\odot}$) lose their disks rapidly (in $\sim 10^5$ years) due to their high EUV and FUV fields. Disk lifetimes are shorter for shallow surface density distributions and when the dust opacity in the disk is reduced by processes such as grain growth or settling. The latter suggests that the photoevaporation process may accelerate as the dust disk evolves.

We report the results of a search for pure rotational molecular hydrogen emission from the circumstellar environments of young stellar objects with disks using the Texas Echelon Cross Echelle Spectrograph (TEXES) on the NASA Infrared Telescope Facility and the Gemini North Observatory. We searched for mid-infrared H2 emission in the S(1), S(2), and S(4) transitions. Keck/NIRSPEC observations of the H2 S(9) transition were included for some sources as an additional constraint on the gas temperature. We detected H2 emission from 6 of 29 sources observed: AB Aur, DoAr 21, Elias 29, GSS 30 IRS 1, GV Tau N, and HL Tau. Four of the six targets with detected emission are class I sources that show evidence for surrounding material in an envelope in addition to a circumstellar disk. In these cases, we show that accretion shock heating is a plausible excitation mechanism. The detected emission lines are narrow (~10 km/s), centered at the stellar velocity, and spatially unresolved at scales of 0.4 arcsec, which is consistent with origin from a disk at radii 10-50 AU from the star. In cases where we detect multiple emission lines, we derive temperatures > 500 K from ~1 M_earth of gas. Our upper limits for the non-detections place upper limits on the amount of H2 gas with T > 500 K of less than a few Earth masses. Such warm gas temperatures are significantly higher than the equilibrium dust temperatures at these radii, suggesting that the gas is decoupled from the dust in the regions we are studying and that processes such as UV, X-ray, and accretion heating may be important.

We investigate the relations among the stellar continuum-subtracted 8 micron polycyclic aromatic hydrocarbon (PAH 8 micron) emission, 24 micron hot dust emission, and 160 micron cold dust emission in fifteen nearby face-on spiral galaxies in the Spitzer Infrared Nearby Galaxies Survey sample. The relation between PAH 8 and 24 micron emission measured in ~2 kpc regions is found to exhibit a significant amount of scatter, and strong spatial variations are observed in the (PAH 8 micron)/24 micron surface brightness ratio. In particular, the (PAH 8 micron)/24 micron surface brightness ratio is observed to be high in the diffuse interstellar medium and low in bright star-forming regions and other locations with high 24 micron surface brightness. PAH 8 micron emission is found to be well-correlated with 160 micron emission on spatial scales of ~2 kpc, and the (PAH 8 micron)/160 micron surface brightness ratio is generally observed to increase as the 160 micron surface brightness increases. These results suggest that the PAHs are associated with the diffuse, cold dust that produces most of the 160 micron emission in these galaxies, and the variations in the (PAH 8 micron)/160 micron ratio may generally be indicative of either the intensity or the spectrum of the interstellar radiation field that is heating both the PAHs and the diffuse interstellar dust.

We present self-consistent models of gas in optically-thick dusty disks and calculate its thermal, density and chemical structure. The models focus on an accurate treatment of the upper layers where line emission originates, and at radii $\gtrsim 0.7$ AU. We present results of disks around $\sim 1{\rm M}_{\odot}$ stars where we have varied dust properties, X-ray luminosities and UV luminosities. We separately treat gas and dust thermal balance, and calculate line luminosities at infrared and sub-millimeter wavelengths from all transitions originating in the predominantly neutral gas that lies below the ionized surface of the disk. We find that the [ArII] 7$\mu$m, [NeII] 12.8$\mu$m, [FeI] 24$\mu$m, [SI] 25$\mu$m, [FeII] 26$\mu$m, [SiII] 35 $\mu$m, [OI] 63$\mu$m and pure rotational lines of H$_2$, H$_2$O and CO can be quite strong and are good indicators of the presence and distribution of gas in disks. We apply our models to the disk around the nearby young star, TW Hya, and find good agreement between our model calculations and observations. We also predict strong emission lines from the TW Hya disk that are likely to be detected by future facilities. A comparison of CO observations with our models suggests that the gas disk around TW Hya may be truncated to $\sim 120 $ AU, compared to its dust disk of 174 AU. We speculate that photoevaporation due to the strong stellar FUV field from TW Hya is responsible for the gas disk truncation.

We use observations of the CI, CII, HI, and H_2 column densities along lines of sight in the Galactic plane to determine the formation rate of H_2 on grains and to determine chemical reaction rates with Polycyclic Aromatic Hydrocarbons. Photodissociation region models are used to find the best fit parameters to the observed columns. We find the H_2 formation rate on grains has a low rate (R ~ 1 x 10^(-17) cm^(3) s^(-1)) along lines of sight with low column density (A_V < 0.25) and low molecular fraction (f_(H_2) < 10^(-4)). At higher column densities (0.25 < A_V <2.13), we find a rate of R ~ 3.5x10^(-17) cm^(3) s^(-1). The lower rate at low column densities could be the result of grain processing by interstellar shocks which may deplete the grain surface area or process the sites of H +H formation, thereby inhibiting H_2 production. Alternatively, the formation rate may be normal, and the low molecular fraction may be the result of lines of sight which graze larger clouds. Such lines of sight would have a reduced H_2 self-shielding compared to the line-of-sight column. We find the reaction C^+ +PAH^- --> C + PAH^0 is best fit with a rate 2.4 x 10^(-7) \Phi_PAH T_2^(-0.5) cm^(3) s^(-1) with T_2= T/100 K and the reaction C^+ + PAH^0 --> C + PAH^+ is best fit with a rate 8.8x 10^(-9)\Phi_PAH cm^(3) s^(-1). In high column density gas we find \Phi_PAH ~ 0.4. In low column density gas, \Phi_PAH is less well constrained with \Phi_PAH ~ 0.2 - 0.4.

This paper is one in a series presenting results obtained within the Formation and Evolution of Planetary Systems (FEPS) Legacy Science Program on the Spitzer Space Telescope. Here we present a study of dust processing and growth in seven protoplanetary disks. Our spectra indicate that the circumstellar silicate dust grains have grown to sizes at least 10 times larger than observed in the interstellar medium, and show evidence for a non-negligible (~5 % in mass fractions) contribution from crystalline species. These results are similar to those of other studies of protoplanetary disks. In addition, we find a correlation between the strength of the amorphous silicate feature and the shape of the spectral energy distribution. This latter result is consistent with the growth and subsequent gravitational settling of dust grains towards the disk mid-plane. Further, we find a change in the relative abundance of the different crystalline species: more enstatite relative to forsterite is observed in the inner warm dust population at ~1 AU, while forsterite dominates in the colder outer regions at ~5 to 15 AU. This change in the relative abundances argues for a localized crystallization process rather than a radial mixing scenario where crystalline silicates are being transported outwards from a single formation region in the hot inner parts of the disk. Last, we report the detection of emission from polycyclic aromatic hydrocarbon molecules in five out of seven sources. We find a tentative PAH band at 8.2 micron, previously undetected in the spectra of disks around low-mass pre-main-sequence stars.

(abbreviated) We report detection with the Spitzer Space Telescope of cool dust surrounding solar type stars. The observations were performed as part of the Legacy Science Program, ``Formation and Evolution of Planetary Systems'' (FEPS). From the overall FEPS sample (Meyer et al. 2006) of 328 stars having ages ~0.003-3 Gyr we have selected sources with 70 um flux densities indicating excess in their spectral energy distributions above expected photospheric emission........ .....The rising spectral energy distributions towards - and perhaps beyond - 70 um imply dust temperatures T_dust <45-85 K for debris in equilibrium with the stellar radiation field. We infer bulk properties such as characteristic temperature, location, fractional luminosity, and mass of the dust from fitted single temperature blackbody models. For >1/3 of the debris sources we find that multiple temperature components are suggested, implying a spatial distribution of dust extending over many tens of AU. Because the disks are dominated by collisional processes, the parent body (planetesimal) belts may be extended as well. Preliminary assessment of the statistics of cold debris around sun-like stars shows that ~10% of FEPS targets with masses between 0.6 and 1.8 Msun and ages between 30 Myr and 3 Gyr exhibit 70 um emission in excess of the expected photospheric flux density. We find that fractional excess amplitudes appear higher for younger stars and that there may be a trend in 70 um excess frequency with stellar mass.

We report observations from the Spitzer Space Telescope (SST) regarding the frequency of 24 micron excess emission toward sun-like stars. Our unbiased sample is comprised of 309 stars with masses 0.7-2.2 Msun and ages from <3 Myr to >3 Gyr that lack excess emission at wavelengths <=8 microns. We identify 30 stars that exhibit clear evidence of excess emission from the observed 24/8 micron flux ratio. The implied 24 micron excesses of these candidate debris disk systems range from 13 % (the minimum detectable) to more than 100 % compared to the expected photospheric emission. The frequency of systems with evidence for dust debris emitting at 24 micron ranges from 8.5-19 % at ages <300 Myr to < 4 % for older stars. The results suggest that many, perhaps most, sun-like stars might form terrestrial planets.

We present SWAS detections of the ground-state 1(10)-1(01) transition of o-H2O at 557 GHz in 18 molecular outflows. These results are combined with ground-based observations of the J=1-0 transitions of 12CO and 13CO obtained at the FCRAO and, for a subset of the outflows, data from ISO. Assuming the SWAS water line emission originates from the same gas traced by CO emission, we find that the outflowing gas in most outflows has an o-H2O abundance relative to H2 of between 10(-7) and 10(-6). Analysis of the water abundance as a function of outflow velocity reveals a strong dependence. The water abundance increases with velocity, and at the highest outflow velocities some outflows have relative o-H2O abundances of order 10(-4). However the mass of gas with such elevated water abundances represents less that 1% of the total outflow gas mass. The ISO LWS observations of high-J rotational lines of CO and the 179.5 micron transition of o-H2O provide evidence for a warmer outflow component than required to produce either the SWAS or FCRAO lines. The mass associated with the ISO emission is similar to that responsible for the highest velocity water emission detected by SWAS. The bulk of the outflowing gas has an abundance of o-H2O well below what would be expected if the gas has passed through a C-shock with shock velocities greater than 10 km/s. Gas-phase water can be depleted in the post-shock gas due to freeze-out onto grain mantles, however the rate of freeze-out is too slow to explain our results. Therefore we believe that only a small fraction of the outflowing molecular gas has passed through shocks strong enough to fully convert the gas-phase oxygen to water. This result has implications for the acceleration mechanism of the molecular gas in these outflows.

HST NICMOS narrowband images of the shocked molecular hydrogen emission in OMC-1 are analyzed to reveal new information on the BN/KL outflow. The outstanding morphological feature of this region is the array of molecular hydrogen ``fingers'' emanating from the general vicinity of IRc2 and the presence of several Herbig-Haro objects. The NICMOS images appear to resolve individual shock fronts. This work is a more quantitative and detailed analysis of our data from a previous paper (Schultz etal. 1999). Line strengths for the H_2 1--0 S(4) plus 2--1 S(6) lines at 1.89 micron are estimated from measurements with the Paschen_alpha continuum filter F190N at 1.90 micron, and continuum measurements at 1.66 and 2.15 micron. We compare the observed H_2 line strengths and ratios of the 1.89 micron and 2.12 micron 1--0 S(1) lines with models for molecular cloud shock waves. Most of the data cannot be fit by J-shocks, but are well matched by C-shocks with shock velocities in the range of 20--45 km/s and preshock densities of 10^{4} - 10^{6} cm^{-3}, similiar to values obtained in larger beam studies which averaged over many shocks. There is also some evidence that shocks with higher densities have lower velocities.

We present Spitzer IRS spectra (R ~600, 10 - 38 micron) of 38 positions in the Galactic Center (GC), all at the same Galactic longitude and spanning plus/minus 0.3 degrees in latitude. Our positions include the Arches Cluster, the Arched Filaments, regions near the Quintuplet Cluster, the ``Bubble'' lying along the same line-of-sight as the molecular cloud G0.11-0.11, and the diffuse interstellar gas along the line-of-sight at higher Galactic latitudes. From measurements of the [O IV], [Ne II], [Ne III], [Si II], [S III], [S IV], [Fe II], [Fe III], and H_2 S(0), S(1), and S(2) lines we determine the gas excitation and ionic abundance ratios. The Ne/H and S/H abundance ratios are ~ 1.6 times that of the Orion Nebula. The main source of excitation is photoionization, with the Arches Cluster ionizing the Arched Filaments and the Quintuplet Cluster ionizing the gas nearby and at lower Galactic latitudes including the far side of the Bubble. In addition, strong shocks ionize gas to O^{+3} and destroy dust grains, releasing iron into the gas phase (Fe/H ~ 1.3 times 10^{-6} in the Arched Filaments and Fe/H ~ 8.8 times 10^{-6} in the Bubble). The shock effects are particularly noticeable in the center of the Bubble, but O$^{+3}$ is present in all positions. We suggest that the shocks are due to the winds from the Quintuplet Cluster Wolf-Rayet stars. On the other hand, the H_2 line ratios can be explained with multi-component models of warm molecular gas in photodissociation regions without the need for H_2 production in shocks.

We have studied the relationship between the star formation rate (SFR) surface density and gas surface density in the spiral galaxy M51a (NGC 5194), using multi-wavelength data obtained as part of the Spitzer Infrared Nearby Galaxies Survey (SINGS). We introduce a new SFR index based on a linear combination of H-alpha emission-line and 24 micron continuum luminosities, that provides reliable extinction-corrected ionizing fluxes and SFR densities over a wide range of dust attenuations. The combination of these extinction-corrected SFR densities with aperture synthesis HI and CO maps has allowed us to probe the form of the spatially-resolved star formation law on scales of 0.5 to 2 kpc. We find that the resolved SFR vs gas surface density relation is well represented by a Schmidt power law, which is similar in form and dispersion to the disk-averaged Schmidt law. We observe a comparably strong correlation of the SFR surface density with the molecular gas surface density, but no significant correlation with the surface density of atomic gas. The best-fitting slope of the Schmidt law varies from N = 1.37 to 1.56, with zeropoint and slope that change systematically with the spatial sampling scale. We tentatively attribute these variations to the effects of areal sampling and averaging of a nonlinear intrinsic star formation law. Our data can also be fitted by an alternative parametrization of the SFR surface density in terms of the ratio of gas surface density to local dynamical time, but with a considerable dispersion.

With the goal of investigating the degree to which the mid-infrared emission traces the star formation rate (SFR), we analyze Spitzer 8 um and 24 um data of star-forming regions in a sample of 33 nearby galaxies with available HST/NICMOS images in the Paschen-alpha (1.8756 um) emission line. The galaxies are drawn from the Spitzer Infrared Nearby Galaxies Survey (SINGS) sample, and cover a range of morphologies and a factor ~10 in oxygen abundance. Published data on local low-metallicity starburst galaxies and Luminous Infrared Galaxies are also included in the analysis. Both the stellar-continuum-subtracted 8 um emission and the 24 um emission correlate with the extinction-corrected Pa-alpha line emission, although neither relationship is linear. Simple models of stellar populations and dust extinction and emission are able to reproduce the observed non-linear trend of the 24 um emission versus number of ionizing photons, including the modest deficiency of 24 um emission in the low metallicity regions, which results from a combination of decreasing dust opacity and dust temperature at low luminosities. Conversely, the trend of the 8 um emission as a function of the number of ionizing photons is not well reproduced by the same models. The 8 um emission is contributed, in larger measure than the 24 um emission, by dust heated by non-ionizing stellar populations, in agreement with previous findings. Two SFR calibrations, one using the 24 um emission and the other using a combination of the 24 um and H-alpha luminosities (Kennicutt et al. 2007), are presented. No calibration is presented for the 8 um emission, because of its significant dependence on both metallicity and environment. The calibrations presented here should be directly applicable to systems dominated by on-going star formation.

(simplified) Results on the properties of warm H2 in 57 normal galaxies are derived from H2 rotational transitions, obtained as part of SINGS. This study extends previous extragalactic surveys of H2, the most abundant constituent of the molecular ISM, to more common systems (L_FIR = e7 to 6e10 L_sun) of all morphological and nuclear types. The S(1) transition is securely detected in the nuclear regions of 86% of SINGS galaxies with stellar masses above 10^9.5 M_sun. The derived column densities of warm H2 (T > ~100 K), even though averaged over kiloparsec-scale areas, are commensurate with those of resolved PDRs; the median of the sample is 3e20 cm-2. They amount to between 1% and >30% of the total H2. The power emitted in the sum of the S(0) to S(2) transitions is on average 30% of the [SiII] line power, and ~4e-4 of the total infrared power (TIR) within the same area for star-forming galaxies, which is consistent with excitation in PDRs. The fact that H2 emission scales tightly with PAH emission, even though the average radiation field intensity varies by a factor ten, can also be understood if both tracers originate predominantly in PDRs, either dense or diffuse. A large fraction of the 25 LINER/Sy targets, however, strongly depart from the rest of the sample, in having warmer H2 in the excited states, and an excess of H2 emission with respect to PAHs, TIR and [SiII]. We propose a threshold in H2 to PAH power ratios, allowing the identification of low-luminosity AGNs by an excess H2 excitation. A dominant contribution from shock heating is favored in these objects. Finally, we detect, in nearly half the star-forming targets, non-equilibrium ortho to para ratios, consistent with FUV pumping combined with incomplete ortho-para thermalization by collisions, or possibly non-equilibrium PDR fronts advancing into cold gas.

Using the new capabilities of the Spitzer Space Telescope and extensive multiwavelength data from the Spitzer Infrared Nearby Galaxies Survey (SINGS), it is now possible to study the infrared properties of star formation in nearby galaxies down to scales equivalent to large HII regions. We are therefore able to determine what fraction of large, infrared-selected star-forming regions in normal galaxies are highly obscured and address how much of the star formation we miss by relying solely on the optical portion of the spectrum. Employing a new empirical method for deriving attenuations of infrared-selected star-forming regions we investigate the statistics of obscured star formation on 500pc scales in a sample of 38 nearby galaxies. We find that the median attenuation is 1.4 magnitudes in H-alpha and that there is no evidence for a substantial sub-population of uniformly highly-obscured star-forming regions. The regions in the highly-obscured tail of the attenuation distribution (A_H-alpha > 3) make up only ~4% of the sample of nearly 1800 regions, though very embedded infrared sources on the much smaller scales and lower luminosities of compact and ultracompact HII regions are almost certainly present in greater numbers. The highly-obscured cases in our sample are generally the bright, central regions of galaxies with high overall attenuation but are not otherwise remarkable. We also find that a majority of the galaxies show decreasing radial trends in H-alpha attenuation. The small fraction of highly-obscured regions seen in this sample of normal, star-forming galaxies suggests that on 500pc scales the timescale for significant dispersal or break up of nearby, optically-thick dust clouds is short relative to the lifetime of a typical star-forming region.

HD 38529 is a post-main sequence G8III/IV star (3.5 Gyr old) with a planetary system consisting of at least two planets having Msin(i) of 0.8 MJup and 12.2 MJup, semimajor axes of 0.13 AU and 3.74 AU, and eccentricities of 0.25 and 0.35, respectively. Spitzer observations show that HD 38529 has an excess emission above the stellar photosphere, with a signal-to-noise ratio (S/N) at 70 micron of 4.7, a small excess at 33 micron (S/N=2.6) and no excess <30 micron. We discuss the distribution of the potential dust-producing planetesimals from the study of the dynamical perturbations of the two known planets, considering in particular the effect of secular resonances. We identify three dynamically stable niches at 0.4-0.8 AU, 20-50 AU and beyond 60 AU. We model the spectral energy distribution of HD 38529 to find out which of these niches show signs of harboring dust-producing plantesimals. The secular analysis, together with the SED modeling resuls, suggest that the planetesimals responsible for most of the dust emission are likely located within 20-50 AU, a configuration that resembles that of the Jovian planets + Kuiper Belt in our Solar System. Finally, we place upper limits (8E-6 lunar masses of 10 micron particles) to the amount of dust that could be located in the dynamically stable region that exists between the two planets (0.25--0.75 AU).

We report the results of spectroscopic mapping observations carried out toward small (1 x 1 arcmin) regions within the supernova remnants W44, W28, IC443, and 3C391 using the Infrared Spectrograph of the Spitzer Space Telescope. These observations, covering the 5.2 - 37 micron spectral region, have led to the detection of a total of 15 fine structure transitions of Ne+, Ne++, Si+, P+, S, S++, Cl+, Fe+, and Fe++; the S(0) - S(7) pure rotational lines of molecular hydrogen; and the R(3) and R(4) transitions of hydrogen deuteride. In addition to these 25 spectral lines, the 6.2, 7.7, 8.6, 11.3 and 12.6 micron PAH emission bands were also observed. Most of the detected line transitions have proven strong enough to map in several sources, providing a comprehensive picture of the relative distribution of the various line emissions observable in the Spitzer/IRS bandpass. A principal component analysis of the spectral line maps reveals that the observed emission lines fall into five distinct groups, each of which may exhibit a distinct spatial distribution: (1) lines of S and H2 (J > 2); (2) the H2 S(0) line; (3) lines of ions with appearance potentials less than 13.6 eV; (4) lines of ions with appearance potentials greater than 13.6 eV, not including S++; (5) lines of S++. Lines of group (1) likely originate in molecular material subject to a slow, nondissociative shock that is driven by the overpressure within the supernova remnant, and lines in groups (3) - (5) are associated primarily with dissociative shock fronts with a range of (larger) shock velocities. The H2 S(0) line shows a low-density diffuse emission component, and - in some sources - a shock-excited component.

We report the detection of [Ne II] emission at 12.81 micron in four out of the six optically thick dust disks observed as part of the FEPS Spitzer Legacy program. In addition, we detect a H I(7-6) emission line at 12.37 micron from the source RXJ1852.3-3700. Detections of [Ne II] lines are favored by low mid-infrared excess emission. Both stellar X-rays and extreme UV (EUV) photons can sufficiently ionize the disk surface to reproduce the observed line fluxes, suggesting that emission from Ne+ originates in the hot disk atmosphere. On the other hand, the H I(7-6) line is not associated with the gas in the disk surface and magnetospheric accretion flows can account only for at most ~30% of the observed flux. We conclude that accretion shock regions and/or the stellar corona could contribute to most of the H I(7-6)emission. Finally, we discuss the observations necessary to identify whether stellar X-rays or EUV photons are the dominant ionization mechanism for Ne atoms. Because the observed [Ne II] emission probes very small amounts of gas in the disk surface (~10^{-6} Jupiter masses) we suggest using this gas line to determine the presence or absence of gas in more evolved circumstellar disks.

Physical dust models are presented for 65 galaxies in the SINGS survey that are strongly detected in the four IRAC bands and three MIPS bands. For each galaxy we estimate (1) the total dust mass, (2) the fraction of the dust mass contributed by PAHs, and (3) the intensity of the starlight heating the dust grains. We find that spiral galaxies have dust properties resembling the dust in the local region of the Milky Way, with similar dust-to-gas ratio, and similar PAH abundance. The observed SEDs, including galaxies with SCUBA photometry, can be reproduced by dust models that do not require "cold" (T<10K) dust. The dust-to-gas ratio is observed to be dependent on metallicity. In the interstellar media of galaxies with A_O=12+log(O/H)>8.1, grains contain a substantial fraction of interstellar Mg, Si and Fe. Galaxies with A_O<8.1 and extended HI envelopes in some cases appear to have global dust-to-gas ratios that are low for their measured oxygen abundance, but in the regions where infrared emission is detected, the dust-to-gas ratio generally appears to be consistent with a substantial fraction of interstellar Mg, Si, and Fe being contained in dust. The PAH index q_PAH -- the fraction of the dust mass in the form of PAHs -- correlates with metallicity. The nine galaxies in our sample with A_O<8.1 have a median q_PAH=1.0%, whereas galaxies with A_O>8.1 have a median q_PAH=3.55%. The derived dust masses favor a value X_CO approx 4e20 cm^{-2}(K kms)^{-1} for the CO to H_2 conversion factor. Except for some starbursting systems (Mrk33, Tolo89, NGC3049), dust in the diffuse ISM dominates the IR power.

[abridged] We present observations of the dust and atomic gas phase in seven dwarf irregular galaxies of the M81 group from the SINGS and THINGS surveys. The Spitzer observations provide a first glimpse of the nature of the non-atomic ISM in these metal-poor (Z~0.1 Z_sun), quiescent (SFR~0.001-0.1 M_sun/yr) dwarf galaxies. Dust emission is detected in five out of the seven targets. Most detected dust emission is restricted to HI column densities >1x10^21 cm^-2. Spitzer spectroscopy of two regions in the brightest galaxies (IC 2574 and Holmberg II) show distinctly different spectral shapes. The spectrum of IC 2574 shows aromatic features that are less luminous (relative to the FIR luminosity) compared to an average SINGS spiral galaxy by a factor of \~7 . The aromatic features in Holmberg~II (which has only a slightly lower gas-phase metallicity) are fainter than in IC 2574 by an order of magnitude. This result emphazises that the strength of the aromatic features is not a simple linear function of metallicity. We estimate dust masses of ~10^4-10^6 M_sun for the M81 dwarf galaxies, resulting in an average dust--to--gas ratio (M_dust/M_HI) of ~3x10^-4 (1.5x10^-3 if only the HI that is associated with dust emission is considered); this is an order of magnitude lower than the typical value derived for the SINGS spirals. The dwarf galaxies are underluminous per unit star formation rate at 70um as compared to the more massive galaxies in SINGS by a factor of ~2. However, the average 70um/160um ratio in the sample dwarf galaxies is higher than what is found in the other galaxies of the SINGS sample. This can be explained by a combination of a lower dust content in conjunction with a higher dust temperature in the dwarfs.

We provide an overview of the Spitzer Legacy Program ``Formation and Evolution of Planetary Systems" (FEPS) which was proposed in 2000, begun in 2001, and executed aboard the Spitzer Space Telescope between 2003 and 2006. This program exploits the sensitivity of Spitzer to carry out mid-infrared spectrophotometric observations of solar-type stars. With a sample of ~ 328 stars ranging in age from ~ 3 Myr to ~ 3 Gyr, we trace the evolution of circumstellar gas and dust from primordial planet-building stages in young circumstellar disks through to older collisionally generated debris disks. When completed, our program will help define the time scales over which terrestrial and gas giant planets are built, constrain the frequency of planetesimal collisions as a function of time, and establish the diversity of mature planetary architectures. In addition to the observational program, we have coordinated a concomitant theoretical effort aimed at understanding the dynamics of circumstellar dust with and without the effects of embedded planets, dust spectral energy distributions, and atomic and molecular gas line emission. Together with the observations, these efforts will provide astronomical context for understanding whether our Solar System - and its habitable planet - is a common or a rare circumstance. Additional information about the FEPS project can be found on the team website: feps.as.arizona.edu

Using data from the Spitzer Space Telescope Legacy Science Program ``Formation and Evolution of Planetary Systems'' (FEPS), we have searched for debris disks around 9 FGK stars (2-10 Gyr), known from radial velocity (RV) studies to have one or more massive planets. Only one of the sources, HD 38529, has excess emission above the stellar photosphere; at 70 micron the signal-to-noise ratio in the excess is 4.7 while at wavelengths < 30 micron there is no evidence of excess. The remaining sources show no excesses at any Spitzer wavelengths. Applying survival tests to the FEPS sample and the results for the FGK survey published in Bryden et al. (2006), we do not find a significant correlation between the frequency and properties of debris disks and the presence of close-in planets. We discuss possible reasons for the lack of a correlation.

We present a sample of low-resolution 5-38um Spitzer IRS spectra of the inner few square kiloparsecs of 59 nearby galaxies spanning a large range of star formation properties. A robust method for decomposing mid-infrared galaxy spectra is described, and used to explore the behavior of PAH emission and the prevalence of silicate dust extinction. Evidence for silicate extinction is found in ~1/8 of the sample, at strengths which indicate most normal galaxies undergo A_V < ~3 magnitudes averaged over their centers. The contribution of PAH emission to the total infrared power is found to peak near 10% and extend up to ~20%, and is suppressed at metallicities Z < ~Z_sun/4, as well as in low-luminosity AGN environments. Strong inter-band PAH feature strength variations (2-5x) are observed, with the presence of a weak AGN and, to a lesser degree, increasing metallicity shifting power to the longer wavelength bands. A peculiar PAH emission spectrum with markedly diminished 5-8um features arises among the sample solely in systems with relatively hard radiation fields harboring low-luminosity AGN. The AGN may modify the emitting grain distribution and provide the direct excitation source of the unusual PAH emission, which cautions against using absolute PAH strength to estimate star formation rates in systems harboring active nuclei. Alternatively, the low star formation intensity often associated with weak AGN may affect the spectrum. The effect of variations in the mid-infrared spectrum on broadband infrared surveys is modeled, and points to more than a factor of two uncertainty in results which assume a fixed PAH emission spectrum, for redshifts z=0-2.5.

The ultraviolet-to-radio continuum spectral energy distributions are presented for all 75 galaxies in the Spitzer Infrared Nearby Galaxies Survey (SINGS). A principal component analysis of the sample shows that most of the sample's spectral variations stem from two underlying components, one representative of a galaxy with a low infrared-to-ultraviolet ratio and one representative of a galaxy with a high infrared-to-ultraviolet ratio. The influence of several parameters on the infrared-to-ultraviolet ratio is studied (e.g., optical morphology, disk inclination, far-infrared color, ultraviolet spectral slope, and star formation history). Consistent with our understanding of normal star-forming galaxies, the SINGS sample of galaxies in comparison to more actively star-forming galaxies exhibits a larger dispersion in the infrared-to-ultraviolet versus ultraviolet spectral slope correlation. Early type galaxies, exhibiting low star formation rates and high optical surface brightnesses, have the most discrepant infrared-to-ultraviolet correlation. These results suggest that the star formation history may be the dominant regulator of the broadband spectral variations between galaxies. Finally, a new discovery shows that the 24 micron morphology can be a useful tool for parametrizing the global dust temperature and ultraviolet extinction in nearby galaxies. The dust emission in dwarf/irregular galaxies is clumpy and warm accompanied by low ultraviolet extinction, while in spiral galaxies there is typically a much larger diffuse component of cooler dust and average ultraviolet extinction. For galaxies with nuclear 24 micron emission, the dust temperature and ultraviolet extinction are relatively high compared to disk galaxies.

We present Spitzer imaging of the metal-deficient (Z ~30% Z_sun) Local Group dwarf galaxy NGC 6822. On spatial scales of ~130 pc, we study the nature of IR, H alpha, HI, and radio continuum emission. Nebular emission strength correlates with IR surface brightness; however, roughly half of the IR emission is associated with diffuse regions not luminous at H alpha (as found in previous studies). The global ratio of dust to HI gas in the ISM, while uncertain at the factor of ~2 level, is ~25 times lower than the global values derived for spiral galaxies using similar modeling techniques; localized ratios of dust to HI gas are about a factor of five higher than the global value in NGC 6822. There are strong variations (factors of ~10) in the relative ratios of H alpha and IR flux throughout the central disk; the low dust content of NGC 6822 is likely responsible for the different H alpha/IR ratios compared to those found in more metal-rich environments. The H alpha and IR emission is associated with high-column density (> ~1E21 cm^-2) neutral gas. Increases in IR surface brightness appear to be affected by both increased radiation field strength and increased local gas density. Individual regions and the galaxy as a whole fall within the observed scatter of recent high-resolution studies of the radio-far IR correlation in nearby spiral galaxies; this is likely the result of depleted radio and far-IR emission strengths in the ISM of this dwarf galaxy.

We explore the nature of variations in dust emission within an individual galaxy using 3.6 - 160 micron Spitzer Space Telescope observations and 450 and 850 micron James Clerk Maxwell Telescope observations of the edge-on Sd spiral galaxy NGC 4631 with the goals of understanding the relation between polycyclic aromatic hydrocarbons (PAH) and dust emission, studying the variations in the colors of the dust emission, and searching for possible excess submillimeter emission compared to what is anticipated based on dust models applied to the mid- and far-infrared data. PAH emission at 8 micron is found to correlate best with hot dust emission at 24 micron on kiloparsec scales, although the relation breaks down on scales equal to hundreds of parsecs, possibly because of differences in the mean free paths between the photons that excite the PAHs and heat the dust and possibly because the PAHs are destroyed by the hard radiation fields in the centers of some star formation regions. The ratio of 8 micron PAH emission to 160 micron cool dust emission appears to be a function of radius. The 70/160 and 160/450 micron flux density ratios are remarkably constant even though the surface brightness varies by a factor of 25 in each wave band, which suggests that the emission is from dust heated by a nearly-uniform radiation field. Globally, we find an excess of 850 micron emission relative to what would be predicted by dust models. The 850 micron excess is highest in regions with low 160 micron surface brightness, although the strength and statistical significance of this result depends on the model fit to the data. We rule out variable emissivity functions or ~4 K dust as the possible origins of this 850 micron emission, but we do discuss the other possible mechanisms that could produce the emission.

We have carried out a sensitive search for gas emission lines at infrared and millimeter wavelengths for a sample of 15 young sun-like stars selected from our dust disk survey with the Spitzer Space Telescope. We have used mid-infrared lines to trace the warm (300-100 K) gas in the inner disk and millimeter transitions of 12CO to probe the cold (~20 K) outer disk. We report no gas line detections from our sample. Line flux upper limits are first converted to warm and cold gas mass limits using simple approximations allowing a direct comparison with values from the literature. We also present results from more sophisticated models following Gorti and Hollenbach (2004) which confirm and extend our simple analysis. These models show that the SI line at 25.23 micron can set constraining limits on the gas surface density at the disk inner radius and traces disk regions up to a few AU. We find that none of the 15 systems have more than 0.04 MJ of gas within a few AU from the disk inner radius for disk radii from 1 AU up to ~40 AU. These gas mass upper limits even in the 8 systems younger than ~30 Myr suggest that most of the gas is dispersed early. The gas mass upper limits in the 10-40 AU region, that is mainly traced by our CO data, are <2 Mearth. If these systems are analogs of the Solar System, either they have already formed Uranus- and Neptune-like planets or they will not form them beyond 100 Myr. Finally, the gas surface density upper limits at 1 AU are smaller than 0.01% of the minimum mass solar nebula for most of the sources. If terrestrial planets form frequently and their orbits are circularized by gas, then circularization occurs early.

We report the detection of interstellar hydrogen deuteride (HD) toward the supernova remnant IC443, and the tentative detection of HD toward the Herbig Haro objects HH54 and HH7 and the star forming region GGD37 (Cepheus A West). Our detections are based upon spectral line mapping observations of the R(3) and R(4) rotational lines of HD, at rest wavelengths of 28.502 and 23.034 micron respectively, obtained using the Infrared Spectrograph onboard the Spitzer Space Telescope. The HD R(4)/R(3) line intensity ratio promises to be a valuable probe of the gas pressure in regions where it can be observed. The derived HD/H2 abundance ratios are 1.19(+0.35/-0.24)E-5, 1.80(+0.54/-0.32)E-5, and 1.41(+0.46/-0.33)E-5 respectively (68.3% confidence limits, based upon statistical errors alone) for IC443 (clump C), HH54, and HH7. If HD is the only significant reservoir of gas-phase deuterium in these sources, the inferred HD/H2 ratios are all consistent with a gas-phase elemental abundance [n(D)/n(H)](gas) ~ 7.5E-6, a factor 2 - 3 below the values obtained previously from observations of atomic deuterium in the local bubble and the Galactic halo. However, similarly low gas-phase deuterium abundances have been inferred previously for molecular gas clouds in the Orion region, and in atomic clouds along sight-lines within the Galactic disk to stars more distant than 500 pc from the Sun.

The recent star formation (SF) in the early-type spiral galaxy M81 is characterized using imaging observations from the far-ultraviolet (UV) to the far-infrared (IR). We compare these data with models of the stellar, gas, and dust emission for sub-galactic regions. Our results suggest the existence of a diffuse dust emission not directly linked to the recent SF. We find a radial decrease of the dust temperature and dust mass density, and in the attenuation of the stellar light. The IR emission in M81 can be modeled with three components: 1) cold dust with a temperature <T_c>=18+-2 K, concentrated near the HII regions but also presenting a diffuse distribution; 2) warm dust with T_w=53+-7 K, directly linked with the HII regions; and 3) aromatic molecules, with diffuse morphology peaking around the HII regions. We derive several relationships to obtain total IR luminosities from IR monochromatic fluxes, and we compare five different star formation rate (SFR) estimators for HII regions in M81 and M51: the UV, Halpha, and three estimators based on Spitzer data. We find that the Halpha luminosity absorbed by dust correlates tightly with the 24 microns emission. The correlation with the total IR luminosity is not as good. Important variations from galaxy to galaxy are found when estimating the total SFR with the 24 microns or the total IR emission alone. The most reliable estimations of the total SFRs are obtained by combining the Halpha emission (or the UV) and an IR luminosity (especially the 24 microns emission), which probe the unobscured and obscured SF, respectively. For the entire M81 galaxy, about 50% of the total SF is obscured by dust. The percentage of obscured SF ranges from 60% in the inner regions of the galaxy to 30% in the outer zones.

We present Spitzer observations of the dwarf starburst galaxy NGC 1705 obtained as part of SINGS. The galaxy morphology is very different shortward and longward of ~5 microns: short-wavelength imaging shows an underlying red stellar population, with the central super star cluster (SSC) dominating the luminosity; longer-wavelength data reveals warm dust emission arising from two off-nuclear regions offset by ~250 pc from the SSC. These regions show little extinction at optical wavelengths. The galaxy has a relatively low global dust mass (~2E5 solar masses, implying a global dust-to-gas mass ratio ~2--4 times lower than the Milky Way average). The off-nuclear dust emission appears to be powered by photons from the same stellar population responsible for the excitation of the observed H Alpha emission; these photons are unassociated with the SSC (though a contribution from embedded sources to the IR luminosity of the off-nuclear regions cannot be ruled out). Low-resolution IRS spectroscopy shows moderate-strength PAH emission in the 11.3 micron band in the eastern peak; no PAH emission is detected in the SSC or the western dust emission complex. There is significant diffuse 8 micron emission after scaling and subtracting shorter wavelength data; the spatially variable PAH emission strengths revealed by the IRS data suggest caution in the interpretation of diffuse 8 micron emission as arising from PAH carriers alone. The metallicity of NGC 1705 falls at the transition level of 35% solar found by Engelbracht and collaborators; the fact that a system at this metallicity shows spatially variable PAH emission demonstrates the complexity of interpreting diffuse 8 micron emission. A radio continuum non-detection, NGC 1705 deviates significantly from the canonical far-IR vs. radio correlation. (Abridged)

We analyze extensive data on NGC1377 from the Spitzer Infrared Nearby Galaxies Survey (SINGS). Within the category of nascent starbursts, that we previously selected by their infrared to radio continuum ratios in large excess of the average and their hot dust, NGC1377 has the largest infrared excess yet measured. Optical imaging reveals a morphological distortion suggestive of a recent accretion event. Infrared spectroscopy reveals a compact and opaque source dominated by a hot, self-absorbed continuum (tau ~ 20 in the 10 micron silicate band). We provide physical evidence against non-stellar activity being the heating source. HII regions are detected through the single [NeII] line, probing <1% of the ionizing radiation. Not only is the optical depth very high, but >85% of ionizing photons are suppressed by dust. The only other detected emission features are molecular hydrogen lines, arguably excited mainly by shocks, besides photodissociation regions, and weak aromatic bands. The new observations support our interpretation in terms of an extremely young starburst (<1 Myr). More generally, galaxies deficient in radio synchrotron are likely observed within a few Myr of the onset of a starburst and after a long quiescence, prior to the replenishment of the ISM with cosmic rays. The similar infrared-radio properties of NGC1377 and some infrared-luminous galaxies suggest that NGC1377 constitutes an archetype to better understand starburst evolution. Although rare locally because observed in a brief evolutionary stage, nascent starbursts may represent a non-negligible fraction of merger-induced starbursts dominating deep infrared counts. Since they differ dramatically from usual starburst templates, they bear important consequences for the interpretation of deep surveys.

Mid-infrared diagnostics are presented for a large portion of the Spitzer Infrared Nearby Galaxies Survey (SINGS) sample plus archival data from the Infrared Space Observatory and the Spitzer Space Telescope. Our main result is that these mid-infrared diagnostics effectively constrain a target's dominant power source. The combination of a high ionization line index and PAH strength serves as an efficient discriminant between AGN and star-forming nuclei, confirming progress made with ISO spectroscopy on starbursting and ultraluminous infrared galaxies. The sensitivity of Spitzer allows us to probe fainter nuclei and star-forming regions within galaxy disks. We find that both star-forming nuclei and extranuclear regions stand apart from nuclei that are powered by Seyfert or LINER activity. In fact, we identify areas within four diagnostic diagrams containing >90% Seyfert/LINER nuclei or >90% HII regions/HII nuclei. We also find that, compared to starbursting nuclei, extranuclear regions typically separate even further from AGN, especially for low-metallicity extranuclear environments. In addition, instead of the traditional mid-infrared approach to differentiating between AGN and star-forming sources that utilizes relatively weak high-ionization lines, we show that strong low-ionization cooling lines of X-ray dominated regions like [SiII] 34.82 micron can alternatively be used as excellent discrimants. Finally, the typical target in this sample shows relatively modest interstellar electron density (~400 cm^{-3}) and obscuration (A_V ~ 1.0 mag for a foreground screen), consistent with a lack of dense clumps of highly obscured gas and dust residing in the emitting regions. [Abridged]

We present new images (groundbased optical and mid-infrared (MIR) from the Spitzer Space Telescope) and spectra (from Spitzer) of the archetypal starburst galaxy M 82. The Spitzer data show that the MIR emission extends at least 6 kpc along the minor axis of the galaxy. We use the optical and infrared data to demonstrate that the extended emission is dominated by emission from dust. The colors of the MIR emission and the spectra indicate that there is a strong component of aromatic feature emission (the MIR features commonly attributed to polycyclic aromatic hydrocarbons). The dust continuum and aromatic feature emission are both strong in the well-known superwind region of this galaxy; clearly the carrier of the aromatic features can survive in close proximity to the wind, far from the plane of the galaxy. We also see significant emission by dust well outside the superwind region, providing the clearest picture to date of the dust distribution in the halo of this galaxy.

We present Spitzer 3.6-160 micron images, Spitzer mid-infrared spectra, and JCMT SCUBA 850 micron images of the Sombrero Galaxy (NGC 4594), an Sa galaxy with a 10^9 M_solar low luminosity active galactic nucleus (AGN). The brightest infrared sources in the galaxy are the nucleus and the dust ring. The spectral energy distribution of the AGN demonstrates that, while the environment around the AGN is a prominent source of mid-infrared emission, it is a relatively weak source of far-infrared emission, as had been inferred for AGN in previous research. The weak nuclear 160 micron emission and the negligible polycyclic aromatic hydrocarbon emission from the nucleus also implies that the nucleus is a site of only weak star formation activity and the nucleus contains relatively little cool interstellar gas needed to fuel such activity. We propose that this galaxy may be representative of a subset of low ionization nuclear emission region galaxies that are in a quiescent AGN phase because of the lack of gas needed to fuel circumnuclear star formation and Seyfert-like AGN activity. Surprisingly, the AGN is the predominant source of 850 micron emission. We examine the possible emission mechanisms that could give rise to the 850 micron emission and find that neither thermal dust emission, CO line emission, bremsstrahlung emission, nor the synchrotron emission observed at radio wavelengths can adequately explain the measured 850 micron flux density by themselves. The remaining possibilities for the source of the 850 micron emission include a combination of known emission mechanisms, synchrotron emission that is self-absorbed at wavelengths longer than 850 microns, or unidentified spectral lines in the 850 micron band.

We review advances in the modeling of protoplanetary disks. This review will focus on the regions of the disk beyond the dust sublimation radius, i.e. beyond 0.1 - 1 AU, depending on the stellar luminosity. We will be mostly concerned with models that aim to fit spectra of the dust continuum or gas lines, and derive physical parameters from these fits. For optically thick disks, these parameters include the accretion rate through the disk onto the star, the geometry of the disk, the dust properties, the surface chemistry and the thermal balance of the gas. For the latter we are mostly concerned with the upper layers of the disk, where the gas and dust temperature decouple and a photoevaporative flow may originate. We also briefly discuss optically thin disks, focusing mainly on the gas, not the dust. The evolution of these disks is dominated by accretion, viscous spreading, photoevaporation, and dust settling and coagulation. The density and temperature structure arising from the surface layer models provide input to models of photoevaporation, which occurs largely in the outer disk. We discuss the consequences of photoevaporation on disk evolution and planet formation.

(Abridged) We present an initial look at the far infrared-radio correlation within the star-forming disks of four nearby, nearly face-on galaxies (NGC 2403, NGC 3031, NGC 5194, and NGC 6946). Using Spitzer MIPS imaging and WSRT radio continuum data, observed as part of the Spitzer Infrared Nearby Galaxies Survey (SINGS), we are able to probe variations in the logarithmic 24mu/22cm (q_24) and 70mu/22cm (q_70) surface brightness ratios across each disk at sub-kpc scales. We find general trends of decreasing q_24 and q_70 with declining surface brightness and with increasing radius. The residual dispersion around the trend of q_24 and q_70 versus surface brightness is smaller than the residual dispersion around the trend of q_24 and q_70 versus radius, on average by ~0.1 dex, indicating that the distribution of star formation sites is more important in determining the infrared/radio disk appearance than the exponential profiles of disks. We have also performed preliminary phenomenological modeling of cosmic ray electron (CRe^-) diffusion using an image-smearing technique, and find that smoothing the infrared maps improves their correlation with the radio maps. Exponential kernels tend to work better than Gaussian kernels which suggests that additional processes besides simple random-walk diffusion in three dimensions must affect the evolution of CRe^-s. The best fit smoothing kernels for the two less active star-forming galaxies (NGC 2403 and NGC 3031) have much larger scale-lengths than those of the more active star-forming galaxies (NGC 5194 and NGC 6946). This difference may be due to the relative deficit of recent CRe^- injection into the interstellar medium (ISM) for the galaxies having largely quiescent disks.

Recent observations of the ground state transition of HDO at 464 GHz towards the protoplanetary disk of DM Tau have detected the presence of water vapor in the regions just above the outer disk midplane (Ceccarelli et al 2005). In the absence of non-thermal desorption processes, water should be almost entirely frozen onto the grain mantles and HDO undetectable. In this Letter we present a chemical model that explores the possibility that the icy mantles are photo-desorbed by FUV (6eV < h nu < 13.6eV) photons. We show that the average Interstellar FUV field is enough to create a layer of water vapor above the disk midplane over the entire disk. Assuming a photo-desorption yield of 10^{-3}, the water abundance in this layer is predicted to be ~ 3 x 10^{-7} and the average H2O column density is ~ 1.6x 10^{15} cm^{-2}. The predictions are very weakly dependent on the details of the model, like the incident FUV radiation field, and the gas density in the disk. Based on this model, we predict a gaseous HDO/H2O ratio in DM Tau of ~1%. In addition, we predict the ground state transition of water at 557 GHz to be undetectable with ODIN and/or HSO-HIFI.

We present spatially resolved Spitzer imaging of the supergiant shell region of the M81 group dwarf galaxy IC 2574 obtained as part of the Spitzer Infrared Nearby Galaxies Survey. This region harbors one of the best nearby examples of a kinematically distinct HI shell, with an associated remnant stellar cluster; the shell is initiating sequential star formation as it interacts with the surrounding interstellar medium. This region dominates the infrared luminosity of IC 2574 and is spatially resolved in all Spitzer imaging bands. We study the differences in dust temperature as a function of local environment and compare local star formation rates as inferred from H Alpha and total infrared luminosities. We find that the strong H Alpha sources are associated with regions of warm dust; however, the most luminous infrared and H Alpha sources are not necessarilyco-spatial. The coolest dust is found in the regions farthest from the rim of the shell; these regions show the best agreement between star formation rates derived from H Alpha and from total infrared luminosities (although discrepancies at the factor of 3-4 level still exist). There is considerable variation in the radio-far infrared correlation in different regions surrounding the shell. The low dust content of the region may influence the scatter seen in these relations; these data demonstrate that the expanding shell is dramatically affecting its surroundings by triggering star formation and altering the dust temperature.

The Spitzer Infrared Nearby Galaxies Survey (SINGS) is carrying out a comprehensive multi-wavelength survey on a sample of 75 nearby galaxies. The 1-850um spectral energy distributions are presented using broadband imaging data from Spitzer, 2MASS, ISO, IRAS, and SCUBA. The infrared colors derived from the globally-integrated Spitzer data are generally consistent with the previous generation of models that were developed based on global data for normal star-forming galaxies, though significant deviations are observed. Spitzer's excellent sensitivity and resolution also allow a detailed investigation of the infrared spectral energy distributions for various locations within the three large, nearby galaxies NGC3031 (M81), NGC5194 (M51), and NGC7331. Strong correlations exist between the local star formation rate and the infrared colors f_nu(70um)/f_nu(160um) and f_nu(24um)/f_nu(160um), suggesting that the 24 and 70um emission are useful tracers of the local star formation activity level. Preliminary evidence indicates that variations in the 24um emission, and not variations in the emission from polycyclic aromatic hydrocarbons at 8um, drive the variations in the f_nu(8.0um)/f_nu(24um) colors within NGC3031, NGC5194, and NGC7331. If the galaxy-to-galaxy variations in spectral energy distributions seen in our sample are representative of the range present at high redshift then extrapolations of total infrared luminosities and star formation rates from the observed 24um flux will be uncertain at the factor-of-five level (total range). The corresponding uncertainties using the redshifted 8.0um flux (e.g. observed 24um flux for a z=2 source) are factors of 10-20. Considerable caution should be used when interpreting such extrapolated infrared luminosities.

(Abridged) Far ultraviolet to far infrared images of the nearby galaxy NGC5194, from Spitzer, GALEX, Hubble Space Telescope and ground--based data, are used to investigate local and global star formation, and the impact of dust extinction in HII-emitting knots. In the IR/UV-UV color plane, the NGC5194 HII knots show the same trend observed for normal star-forming galaxies, having a much larger dispersion than starburst galaxies. We identify the dispersion as due to the UV emission predominantly tracing the evolved, non-ionizing stellar population, up to ages 50-100 Myr. While in starbursts the UV light traces the current SFR, in NGC5194 it traces a combination of current and recent-past SFR. Unlike the UV emission, the monochromatic 24 micron luminosity is an accurate local SFR tracer for the HII knots in NGC5194; this suggests that the 24 micron emission carriers are mainly heated by the young, ionizing stars. However, preliminary results show that the ratio of the 24 micron emission to the SFR varies by a factor of a few from galaxy to galaxy. While also correlated with star formation, the 8 micron emission is not directly proportional to the number of ionizing photons. This confirms earlier suggestions that the carriers of the 8 micron emission are heated by more than one mechanism.

We report infrared spectroscopic observations of HD 105, a nearby ($\sim 40$ pc) and relatively young ($\sim 30$ Myr) G0 star with excess infrared continuum emission, which has been modeled as arising from an optically thin circumstellar dust disk with an inner hole of size $\gtrsim 13$ AU. We have used the high spectral resolution mode of the Infrared Spectrometer (IRS) on the Spitzer Space Telescope to search for gas emission lines from the disk. The observations reported here provide upper limits to the fluxes of H$_2$ S(0) 28$\mu$m, H$_2$ S(1) 17$\mu$m, H$_2$ S(2) 12 $\mu$m, [FeII] 26$\mu$m, [SiII] 35$\mu$m, and [SI] 25$\mu$m infrared emission lines. The H$_2$ line upper limits directly place constraints on the mass of warm molecular gas in the disk: $M({\rm H_2})< 4.6$, 3.8$\times 10^{-2}$, and $3.0\times 10^{-3}$ M$_J$ at $T= 50$, 100, and 200 K, respectively. We also compare the line flux upper limits to predictions from detailed thermal/chemical models of various gas distributions in the disk. These comparisons indicate that if the gas distribution has an inner hole with radius $r_{i,gas}$, the surface density at that inner radius is limited to values ranging from $\lesssim 3$ gm cm$^{-2}$ at $r_{i,gas}=0.5$ AU to 0.1 gm cm$^{-2}$ at $r_{i,gas}= 5-20$ AU. These values are considerably below the value for a minimum mass solar nebula, and suggest that less than 1 M$_J$ of gas (at any temperature) exists in the 1-40 AU planet-forming region. Therefore, it is unlikely that there is sufficient gas for gas giant planet formation to occur in HD 105 at this time.

The nearby spiral galaxy NGC 7331 was spectrally mapped from 5-38um using all modules of Spitzer's IRS spectrograph. A strong new dust emission feature, presumed due to PAHs, was discovered at 17.1um. The feature's intensity is nearly half that of the ubiquitous 11.3um band. The 7-14um spectral maps revealed significant variation in the 7.7 and 11.3um PAH features between the stellar ring and nucleus. Weak [OIV] 25.9um line emission was found to be centrally concentrated in the nucleus, with an observed strength over 10% of the combined neon line flux, indicating an AGN or unusually active massive star photo-ionization. Two [SIII] lines fix the characteristic electron density in the HII regions at n_e < ~200 cm^-3. Three detected H_2 rotational lines, tracing warm molecular gas, together with the observed IR continuum, are difficult to match with standard PDR models. Either additional PDR heating or shocks are required to simultaneously match lines and continuum.

We present 3-160 micron photometry obtained with the IRAC and MIPS instruments for the first five targets from the Spitzer Legacy Science Program "Formation and Evolution of Planetary Systems" and 4-35 micron spectro-photometry obtained with the IRS for two sources. We discuss in detail our observations of the debris disks surrounding HD 105 (G0V, 30 +- 10 Myr) and HD 150706 (G3V, ~ 700 +- 300 Myr). For HD 105, possible interpretations include large bodies clearing the dust inside of 45 AU or a reservoir of gas capable of sculpting the dust distribution. The disk surrounding HD 150706 also exhibits evidence of a large inner hole in its dust distribution. Of the four survey targets without previously detected IR excess, spanning ages 30 Myr to 3 Gyr, the new detection of excess in just one system of intermediate age suggests a variety of initial conditions or divergent evolutionary paths for debris disk systems orbiting solar-type stars.

We model gas and dust emission from regions 0.3-20 AU from a central low mass star in intermediate-aged(~10Myr) disks with optically thin dust. The models treat thermal balance and chemistry self-consistently and calculate the vertical density temperature structure of the gas in a disk. The gas and dust temperatures are calculated separately. We focus on IR gas emission lines and find that the [SI]25.23 micron line is the strongest emission line for a wide range of disk and stellar parameters, followed by emission from [SiII]34.8 um, [FeII]26 um and [OI]63um. [FeI]24um and rotational lines of OH and H2O are strong when gas masses are high (>0.1 M_J). Emission from the rotational lines of H2 is difficult to detect, unless disk gas masses are substantial (>0.1 M_J).

When stars form in small groups (N = 100 - 500 members), their circumstellar disks are exposed to little EUV radiation but a great deal of FUV radiation from massive stars in the group. This paper calculates mass loss rates for circumstellar disks exposed to external FUV radiation. Previous work treated large disks and/or intense radiation fields in which the disk radius exceeds the critical radius (supercritical disks) where the sound speed in the FUV heated layer exceeds the escape speed. This paper shows that significant mass loss still takes place for subcritical systems. Some of the gas extends beyond the disk edge (above the disk surface) to larger distances where the temperature is higher, the escape speed is lower, and an outflow develops. The evaporation rate is a sensitive function of the stellar mass and disk radius, which determine the escape speed, and the external FUV flux, which determines the temperature structure of the flow. Disks around red dwarfs are readily evaporated and shrink to disk radii of 15 AU on short time scales (10 Myr) when exposed to moderate FUV fields with $G_0$ = 3000. Although disks around solar type stars are more durable, these disks shrink to 15 AU in 10 Myr for intense FUV radiation fields with $G_0$ = 30,000; such fields exist in the central 0.7 pc of a cluster with N = 4000 stars. If our solar system formed in the presence of such strong FUV radiation fields, this mechanism could explain why Neptune and Uranus in our solar system are gas poor, whereas Jupiter and Saturn are gas rich. This mechanism for photoevaporation can also limit the production of Kuiper belt objects and can suppress giant planet formation in sufficiently large clusters, such as the Hyades, especially for disks associated with low mass stars.

We have observed 6 clouds along the line-of-sight toward W49A using the Submillimeter Wave Astronomy Satellite (SWAS) and several ground-based observatories. The ortho-H2O 1-0 and OH (1665 and 1667 MHz) transitions are observed in absorption, whereas the low-J CO, 13CO, and C18O lines, as well as the [CI] 1-0 transition, are seen in emission. By using both the o-H218O and o-H2O absorption lines, we are able to constrain the column-averaged o-H_2O abundances in each line-of-sight cloud to within about an order of magnitude. Assuming the standard N(H2)/N(CO) ratio of 10^4, we find N(o-H2O)/N(H2) = 8.1 x 10^-8 - 4 x 10^-7 for three clouds with optically thin water lines. In three additional clouds, the H$_2$O lines are saturated so we have used observations of the H218O ground-state transition to find upper limits to the water abundance of 8.2x 10^-8 - 1.5x10^-6. We measure the OH abundance from the average of the 1665 and 1667 MHz observations and find N(OH)/N(H2) = 2.3x10^-7 - 1.1x10^-6. The o-H2O and OH abundances are similar to those determined for line-of-sight water absorption features towards W51 and Sgr B2 but are higher than those seen from water emission lines in molecular clouds. However, the clouds towards W49 have lower ratios of OH relative to H2O column densities than are predicted by simple models which assume that dissociative recombination is the primary formation pathway for OH and H2O. Building on the work of Neufeld et al. (2002), we present photo-chemistry models including additional chemical effects, which can also explain the observed OH and H2O column densities as well as the observed H2O/CO abundance ratios.

The Palomar Integral Field Spectrograph was used to probe a variety of environments in nine nearby galaxies that span a range of morphological types, luminosities, metallicities, and infrared-to-blue ratios. For the first time, near-infrared spectroscopy was obtained for nuclear or bright HII regions in star-forming galaxies over two spatial dimensions (5.7"x10.0") in the [FeII] (1.257um), [FeII] (1.644um), Paschen beta (1.282um), H2 (2.122um), and Brackett gamma (2.166um) transition lines. These data yield constraints on various characteristics of the star-forming episodes in these regions, including their strength, maturity, spatial variability, and extinction. The HII regions stand out from the nuclei. Unlike observations of nuclear regions, HII region near-infrared observations do not show a spatial coincidence of the line and continuum emission; the continuum and line maps of HII regions usually show distinct and sometimes spatially-separated morphologies. Gauging from Paschen beta and Brackett gamma equivalent widths and luminosities, the HII regions have younger episodes of star formation than the nuclei and more intense radiation fields. Near-infrared line ratio diagnostics suggest that HII regions have "purer" starbursting properties. The correlation between ionizing photon density and mid-infrared color is consistent with the star formation activity level being higher for HII regions than for nuclei. And though the interpretation is complicated, on a purely empirical basis the HII regions show lower Fe1+ abundances than nuclei by an order of magnitude.

The SIRTF Nearby Galaxy Survey is a comprehensive infrared imaging and spectroscopic survey of 75 nearby galaxies. Its primary goal is to characterize the infrared emission of galaxies and their principal infrared-emitting components, across a broad range of galaxy properties and star formation environments. SINGS will provide new insights into the physical processes connecting star formation to the interstellar medium properties of galaxies, and provide a vital foundation for understanding infrared observations of the distant universe and ultraluminous and active galaxies. The galaxy sample and observing strategy have been designed to maximize the scientific and archival value of the data set for the SIRTF user community at large. The SIRTF images and spectra will be supplemented by a comprehensive multi-wavelength library of ancillary and complementary observations, including radio continuum, HI, CO, submillimeter, BVRIJHK, H-alpha, Paschen-alpha, ultraviolet, and X-ray data. This paper describes the main astrophysical issues to be addressed by SINGS, the galaxy sample and the observing strategy, and the SIRTF and other ancillary data products.

Much of the interstellar medium in disk galaxies is in the form of neutral atomic hydrogen, H I. This gas can be in thermal equilibrium at relatively low temperatures, T < 300 K (the cold neutral medium, or CNM) or at temperatures somewhat less than 10^4 K (the warm neutral medium, or WNM). These two phases can coexist over a narrow range of pressures, P_min < P < P_max. We determine P_min and P_max in the plane of the Galaxy as a function of Galactocentric radius R using recent determinations of the gas heating rate and the gas phase abundances of interstellar gas. We provide an analytic approximation for P_min as a function of metallicity, far-ultraviolet radiation field, and the ionization rate of atomic hydrogen. Over most of the disk of the Galaxy, the H I must be in two phases: the weight of the H I in the gravitational potential of the Galaxy is large enough to generate thermal pressures exceeding P_min, so that turbulent pressure fluctuations can produce cold gas that is thermally stable; and the mean density of the H I is too low for the gas to be all CNM. Our models predict the presence of CNM gas to R = 16-18 kpc, somewhat farther than previous estimates. We also examine the potential impact of turbulent heating on our results and provide expressions for the heating rate as a function of Galactic radius.

Observations of W51 with the Submillimeter Wave Astronomy Satellite (SWAS) have yielded the first detection of water vapor in a diffuse molecular cloud. The water vapor lies in a foreground cloud that gives rise to an absorption feature at an LSR velocity of 6 km/s. The inferred H2O column density is 2.5E+13 cm-2. Observations with the Arecibo radio telescope of hydroxyl molecules at ten positions in W51 imply an OH column density of 8E+13 cm-2 in the same diffuse cloud. The observed H2O/OH ratio of ~ 0.3 is significantly larger than an upper limit derived previously from ultraviolet observations of the similar diffuse molecular cloud lying in front of HD 154368. The observed variation in H2O/OH likely points to the presence in one or both of these clouds of a warm (T > 400) gas component in which neutral-neutral reactions are important sources of OH and/or H2O.

(Abridged) We present the analysis of the main FIR fine structure lines emission in NGC1313 and NGC6946. We calculate that a component probably associated with the diffuse disks contributes <~40% in N6946 and ~30 % in N1313 to the total [CII] emission. The main PDR physical parameters responsible for the neutral atomic gas emission in N1313 and N6946 do not significantly differ from what Malhotra etal (2001) found by modelling the integrated emission of a sample of 60 normal galaxies,although there are evidences for a beam averaged contribution of a less active component inside NGC6946 higher than its contribution in the integrated emission of normal galaxies. CO and [CII] in N6946 are well correlated with a mean [CII]/CO ratio similar to that of the normal galaxies sample. In N1313 the [CII]/CO seems to systematically increase from the North to the south, along the S-shaped spiral arm, indicating much more inhomogeneous conditions than in N6946. HI and [CII] in N6946 are completely de-correlated, probably because they arise from different gas components. In N1313 we successfully detect two distinct gas components: a cirrus-like component where HI and [CII] are weakly correlated as observed in our Galaxy, and a component associated with dense PDRs completely de-correlated from HI as observed in N6946.Finally, we find that the HI residing in dense PDRs and presumably recently photo-dissociated, constitutes a few % of the total HI. In turn, this dense gas component produces most of the [CII] emission emitted by the atomic neutral medium.

We present the results of an investigation of the effects of Far Ultraviolet (FUV) radiation from hot early type OB stars on clumps in star-forming molecular clouds. Clumps in Photodissociation regions (PDRs) undergo external heating which, if rapid, creates strong photoevaporative mass flows off the clump surfaces, and drives shocks into the clumps, compressing them to high densities. The clumps lose mass on relatively short timescales. The evolution of an individual clump is found to be sensitive to its initial colunm density, the temperature of the heated surface and the ratio of the ``turn-on time'' $t_{FUV}$ of the heating flux on a clump to its initial sound crossing-time $t_{c}$. In this paper, we use spherical 1-D numerical hydrodynamic models as well as approximate analytical models to study the evolution of turbulence-generated and pressure-confined clumps in PDRs. Turbulent clumps evolve so that their column densities are equal to a critical value determined by the local FUV field, and typically have short photoevaporation timescales, $\sim 10^{4-5}$ years for a 1 M$_{\odot}$ clump in a typical star-forming region. Clumps that are confined by an interclump medium may either get completely photoevaporated, or may preserve a shielded core with a warm, dissociated, protective shell that absorbs the incident FUV flux. We compare our results with observations of some well-studied PDRs: the Orion Bar, M17SW, NGC 2023 and the Rosette Nebula. The data are consistent with both interpretations of clump origin, with a slight indication for favouring the turbulent model for clumps over pressure-confined clumps.

Far Ultraviolet radiation has been recognized as the main source of heating of the neutral interstellar gas, and, as a consequence, it determines whether the thermal balance of the neutral gas results in cold ($T\sim 50 - 100 K$) clouds (CNM), warm ($T \sim 10^4 K$) clouds (WNM), or a combination of the two. High FUV fields convert the neutral gas to WNM, while low fields result in CNM. The sources of FUV radiation are the short-lived massive stars that generally originate in associations that form in Giant Molecular Clouds present in the galactic disk. Using McKee & Williams' (1997) distribution of birthrates for OB associations in the Galaxy, we determine the expected behavior of the time-dependent FUV field for random positions in the local ISM. The FUV field is calculated in two bands (912-1100 \AA and 912-2070 \AA) and at the wavelength 1400 \AA. Our median value for the [912-2070 \AA] band is $G_0=1.6$ times Habing's (1968) value for the radiation field at the solar circle in this band, and quite close to Draine's (1976) value, $G_0=1.7$. Due to attenuation by dust, only associations within about 500 pc contribute significantly to the energy density at a given point. The FUV field is asymmetric at a given point, and the asymmetry grows for higher fields. The FUV field fluctuates with a variety of amplitudes, the larger ones being less frequent. The mean field is about twice the median field because of these fluctuations, or spikes, in the radiation field. These spikes, which last $\sim 30$ Myr, are caused by the infrequent birth of nearby associations. We also model shorter duration spikes caused by runaway OB stars. The presence of a fluctuating heating rate created by the fluctuating FUV field converts CNM to WNM and vice versa.

We will utilize the sensitivity of SIRTF through the Legacy Science Program to carry out spectrophotometric observations of solar-type stars aimed at (1) defining the timescales over which terrestrial and gas giant planets are built, from measurements diagnostic of dust/gas masses and radial distributions; and (2) establishing the diversity of planetary architectures and the frequency of planetesimal collisions as a function of time through observations of circumstellar debris disks. Together, these observations will provide an astronomical context for understanding whether our solar system - and its habitable planet - is a common or a rare circumstance. Achieving our science goals requires measuring precise spectral energy distributions for a statistically robust sample capable of revealing evolutionary trends and the diversity of system outcomes. Our targets have been selected from two carefully assembled databases of solar-like stars: (1) a sample located within 50 pc of the Sun spanning an age range from 100-3000 Myr for which a rich set of ancillary measurements (e.g. metallicity, stellar activity, kinematics) are available; and (2) a selection located between 15 and 180 pc and spanning ages from 3 to 100 Myr. For stars at these distances SIRTF is capable of detecting stellar photospheres with SNR >30 at lambda < 24 microns for our entire sample, as well as achieving SNR >5 at the photospheric limit for over 50% of our sample at lambda=70 microns. Thus we will provide a complete census of stars with excess emission down to the level produced by the dust in our present-day solar system. More information concerning our program can be found at: http://gould.as.arizona.edu/feps

Since 1995, astronomers have discovered planets with masses comparable to that of Jupiter (318 times Earth's mass) in orbit around approximately 60 stars. Although unseen directly, the presence of these planets is inferred by the small reflex motions that they gravitationally induce on the star they orbit; these result in small periodic wavelength shifts in the stellar spectrum. Since this method favors the detection of massive objects orbiting in close proximity to the star, the question of whether these systems also contain analogs of the smaller constituents of our Solar System has remained unanswered. Using an alternative approach, we report here observations of an aging carbon-star, IRC+10216, that reveal the presence of circumstellar water vapor, a molecule not expected in measurable abundances around such a star and thus a distinctive signature of an orbiting cometary system. The only plausible explanation for this water vapor is that the recent evolution of IRC+10216 - which is accompanied by a prodigious increase in its luminosity - is now causing the vaporization of a collection of orbiting icy bodies, a process first considered in a previous theoretical study.

The most important cooling lines of the neutral interstellar medium (ISM) lie in the far-infrared (FIR). We present measurements by the Infrared Space Observatory Long Wavelength Spectrometer of seven lines from neutral and ionized ISM of 60 normal, star-forming galaxies. The galaxy sample spans a range in properties such as morphology, FIR colors (indicating dust temperature), and FIR/Blue ratios (indicating star-formation activity and optical depth). In two-thirds of the galaxies in this sample, the [CII] line is proportional to FIR dust continuum. The other one-third show a smooth decline in [CII]/FIR with increasing F60/F100 and FIR/B, spanning a range of a factor of more than 50. Two galaxies, at the warm and active extreme of the range have [CII]/FIR < 2 \times 10^{-4} (3-sigma upper limit). This is due to increased positive grain charge in the warmer and more active galaxies, which leads to less efficient heating by photoelectrons from dust grains. The ratio of the two principal photodissociation region (PDR) cooling lines [CII]/[OI] shows a tight correlation with F60/F100, indicating that both gas and dust temperatures increase together. We derive a theoretical scaling between [NII] and [CII] from ionized gas and use it to separate [CII] emission from neutral PDRs and ionized gas. Comparison of PDR models of Kaufman et al. (1999) with observed ratios of (a) [OI]/[CII] and ([CII]+[OI])/FIR and (b) [OI]/FIR and F60/F100 yields far-UV flux G0 and gas density n. The derived G0 scales as n to the power 1.4. We interpret this correlation as arising from Stromgren sphere scalings if much of the line and continuum luminosity arises near star-forming regions. The differences in G0 and n may be due to differences in the physical properties of the star-forming clouds.(Short abstract)

We present mid-infrared imaging and far-infrared (FIR) spectroscopy of 5 IBm galaxies observed by ISO as part of our larger study of the inter- stellar medium of galaxies. Most of the irregulars in our sample are very actively forming stars. The mid-infrared imaging was in a band centered at 6.75 microns that is dominated by polycyclic aromatic hydrocarbons (PAHs) and in a band centered at 15 microns that is dominated by small dust grains. The spectroscopy of 3 of the galaxies includes [CII]158 microns and [OI]63 microns, important coolants of photodissociation regions (PDRs), and [OIII]88 microns and [NII]122 microns, which come from ionized gas. [OI]145 microns and [OIII]52 microns were measured in one galaxy as well. These data are combined with PDR and HII region models to deduce properties of the interstellar medium of these galaxies.

We report a strong correlation between the [CII] 158 micron cooling line and the mid-infrared flux in the 5-10 micron range in a wide variety of star-forming galaxies. The mid-infrared flux is dominated by Aromatic Feature Emission (AFE), which is thought to arise from large polycyclic aromatic hydrocarbon molecules or `PAHs' and generally associated with the smallest interstellar grains. The [CII] line is the dominant gas coolant in most regions of atomic interstellar gas, and therefore reflects the heating input to the gas. The ratio of these two quantities, [CII]/AFE, remains nearly constant with the ratio of the IRAS 60 micron band flux to the 100 micron band flux, R(60/100). This is in contrast to the drop in the [CII]/FIR ratio with increasing R(60/100), which signal higher dust temperatures and more intense radiation fields. We interpret the stable [CII]/AFE ratio as evidence that gas heating is dominated by the PAHs or small grains which are also AFE carriers over a wide range of conditions. The trend of decreasing [CII]/FIR with increasing radiation field suggests a decrease in the importance of PAHs or small grains relative to large grains both in gas heating and in dust cooling. We summarize the observed trends and suggest two plausible scenarios.

We mapped a 13 by 22 pc region in emission from 492 GHz [C I] and, for the first time, 551 GHz ^{13}CO(5-4) in the giant molecular cloud M17SW, using the Submillimeter Wave Astronomy Satellite. The morphologies of the [C I] and ^{13}CO emission are strikingly similar. The extent and intensity of the [C I] and ^{13}CO(5-4) emission is explained as arising from photodissociation regions on the surfaces of embedded molecular clumps. Modeling of the ^{13}CO(5-4) emission in comparison to ^{13}CO(1-0) indicates a temperature gradient across the cloud, peaking to at least 63 K near the M17 ionization front and decreasing to at least 20 K at the western edge of the cloud. We see no correlation between gas density and column density. The beam-averaged column density of C I in the core is 1x10^{18} cm^-2, and the mean column density ratio N(C I)/N(CO) is about 0.4. The variations of N(C I)/N(CO) with position in M17SW indicate a similar clump size distribution throughout the cloud.

Four IRAS-detected early type galaxies were observed with ISO. With the exception of the 15 micron image of NGC1052, the mid-IR emission from NGC1052, NGC1155, NGC5866 and NGC6958 at 4.5, 7 and 15 microns show extended emission. Mid-IR emission from NGC1052, NGC1155, and NGC6958 follows a de Vaucouleurs profile. The ratio of 15/7 micron flux decreases with radius in these galaxies, approaching the values empirically observed for purely stellar systems. In NGC5866, the 7 and 15 micron emission is concentrated in the edge-on dust lane. All the galaxies are detected in the [CII] line, and the S0s NGC1155 and NGC5866 are detected in the [OI] line as well. The ISO-LWS observations of the [CII] line are more sensitive measures of cool, neutral ISM than HI and CO by about a factor of 10-100. Three of four early type galaxies, namely NGC1052, NGC6958 and NGC5866, have low ratio FIR/Blue and show a lower [CII]/FIR, which is due to a softer radiation field from old stellar populations. The low [CII]/CO ratio in NGC5866 ([CII]/CO(1-0) < 570) confirms this scenario. We estimate the UV radiation expected from the old stellar populations in these galaxies and compare it to that needed to heat the gas to account for the cooling observed [CII] and [OI] lines. In three out of four galaxies, NGC1052, NGC5866 and NGC6958, the predicted UV radiation falls short by a factor of 2-3. In view of the observed intrinsic scatter in the "UV-upturn" in elliptical galaxies and its great sensitivity to age and metallicity effects, this is not significant. However, the much larger difference (about a factor of 20) between the UV radiation from old stars and that needed to produce the FIR lines for NGC 1155 is strong evidence for the presence of young stars, in NGC1155.

We present new calculations on the contribution from cooling hot gas to the photoionization of warm ionized gas in the Galaxy. We show that hot gas in cooling supernova remnants (SNRs) is an important source of photoionization, particularly for gas in the halo. We find that in many regions at high latitude this source is adequate to account for the observed ionization so there is no need to find ways to transport stellar photons from the disk. The flux from cooling SNRs sets a floor on the ionization along any line of sight. Our model flux is also shown to be consistent with the diffuse soft X-ray background and with soft X-ray observations of external galaxies. We consider the ionization of the clouds observed towards the halo star HD 93521, for which there are no O stars close to the line of sight. We show that the observed ionization can be explained successfully by our model EUV/soft X-ray flux from cooling hot gas. In particular, we can match the H alpha intensity, the S++/S+ ratio, and the C+* column. From observations of the ratios of columns of C+* and either S+ or H0, we are able to estimate the thermal pressure in the clouds. The slow clouds require high (~10^4 cm^-3 K) thermal pressures to match the N(C+*)/N(S+) ratio. Additional heating sources are required for the slow clouds to maintain their ~7000 K temperatures at these pressures, as found by Reynolds, Hausen & Tufte (1999).

Photodissociation Region (PDR) models are computed over a wide range of physical conditions, from those appropriate to giant molecular clouds illuminated by the interstellar radiation field to the conditions experienced by circumstellar disks very close to hot massive stars. These models use the most up-to-date values of atomic and molecular data, the most current chemical rate coefficients, and the newest grain photoelectric heating rates which include treatments of small grains and large molecules. In addition, we examine the effects of metallicity and cloud extinction on the predicted line intensities. Results are presented for PDR models with densities over the range n=10^1-10^7 cm^-3 and for incident far-ultraviolet radiation fields over the range G_0=10^-0.5-10^6.5, for metallicities Z=1 and 0.1 times the local Galactic value, and for a range of PDR cloud sizes. We present line strength and/or line ratio plots for a variety of useful PDR diagnostics: [C II] 158 micron, [O I] 63 and 145 micron, [C I] 370 and 609 micron, CO J=1-0, J=2-1, J=3-2, J=6-5 and J=15-14, as well as the strength of the far-infrared continuum. These plots will be useful for the interpretation of Galactic and extragalactic far infrared and submillimeter spectra observable with ISO, SOFIA, SWAS, FIRST and other orbital and suborbital platforms. As examples, we apply our results to ISO and ground based observations of M82, NGC 278, and the Large Magellenic Cloud.

This paper presents new near-infrared observations of the planetary nebula NGC 2346. The data include a broad K-band image, an image in the H2 vibrationally excited 1-0S(1) line and K band slit spectra at three positions in the nebula. In the H2 1-0S(1) line, the nebula is characterized by a central, bright torus, surrounded by weaker emission with a typical butterfly shape, as seen in Halpha and CO lines. The K band spectra show 11 H2 lines with excitation energies from 6150 to 12552 K. The H2 data have been compared to the predictions of models which follow the evolution with time of the H2 emission in PNe of different core mass and shell properties (Natta & Hollenbach 1998). These models compute the emission originating in the photodissociation region (PDR) created at the inner edge of the neutral shell by the UV radiation of the central core, as well as the emission in the shock associated to the expansion of the shell inside the precursor red-giant wind. In NGC 2346, a PDR origin of the H2 emission in a low-density molecular shell (n<~10^4 cm^-3) is indicated. At these low densities, time-dependent H2 chemistry and X-ray heating of the neutral gas enhance the predicted PDR H2 line intensity by large factors.

We report measurements of the [CII] fine structure line at 157.714 micron in 30 normal star-forming galaxies with the Long Wavelength Spectrometer (LWS) on the Infrared Space Observatory (ISO). The ratio of the line to total far-infrared luminosity, [CII]/FIR, measures the ratio of the cooling of gas to that of dust; and thus the efficiency of the grain photoelectric heating process. This ratio varies by a factor of 40 in the current sample. About two-thirds of the galaxies have [CII]/FIR = 2-7 x 10^{-3}. The other one-third show trends of decreasing [CII]/FIR with increasing dust temperature, indicated by the ratio of infrared emission at 60 and 100 microns; and with increasing star-formation activity, measured by the ratio of far-infrared and blue band luminosity. We find three FIR bright galaxies with undetected [CII] line at 3-sigma upper limits of [CII]/FIR < 0.5-2 x 10^{-4}. The trend in the [CII]/FIR ratio with the temperature of dust and with star-formation activity may be due to decreased efficiency of photoelectric heating of gas at high UV radiation intensity as dust grains become positively charged, decreasing the yield and the energy of the photoelectrons. The three galaxies with no observed PDR lines have among the highest ratios of FIR to Blue luminosity and the ratio of 60 and 100 micron emission. Their lack of [CII] lines may be due to a continuing trend of decreasing [CII]/FIR with increasing star-formation activity and dust temperature seen in one-third of the sample with warm IRAS colors. In that case the upper limits on [CII]/FIR imply a ratio of UV flux to gas density G_0/n > 10 cm$^3. The low [CII]/FIR could also be due to either weak [CII] because of self-absorption or strong FIR continuum from regions weak in [CII], such as dense HII regions or plasma ionized by hard radiation of AGNs.

In this paper we discuss the propagation of dust through the interstellar medium (ISM), and describe the destructive effects of stellar winds, jets, and supernova shock waves on interstellar dust. We review the probability that grains formed in stellar outflows or supernovae survive processing in and propagation through the ISM, and incorporate themselves relatively unprocessed into meteoritic bodies in the solar system. We show that very large (radii >= 5 micron) and very small grains (radii <= 100 Angstrom) with sizes similar to the pre-solar SiC and diamond grains extracted from meteorites, can survive the passage through 100\kms shock waves relatively unscathed. High velocity (>= 250 km/s) shocks destroy dust efficiently. However, a small (~10%) fraction of the stardust never encountered such fast shocks before incorporation into the solar system. All grains should therefore retain traces of their passage through interstellar shocks during their propagation through the ISM. The grain surfaces should show evidence of processing due to sputtering and pitting due to small grain cratering collisions on the micron-sized grains. This conclusion seems to be in conflict with the evidence from the large grains recovered from meteorites which seem to show little interstellar processing.

We analyze the new mid-infrared maps of NGC 6946 for variations in the color ratio of the 7-to-15 micron emission. Our preliminary findings are that this mid-infrared color is remarkably constant between arms and inter-arm regions, and as a function of radius in the disk, excluding the nuclear region. As surface brightness ranges by more than an order of magnitude and the radius runs from about 0.5 to 3kpc, the color ratio remains constant to about +/-20%. Our interpretation is that (1) hard UV radiation from OB stars does not dominate the heating of the grains radiating in the mid-infrared; and (2) that surface brightness variations are driven primarily by surface-filling fraction in the disk, and by radiation intensity increases in starburst environments, such as the nucleus of NGC 6946.

The nearby spiral galaxy NGC 6946 was observed with ISO-CAM in the mid-infrared, achieving 7 arcsec resolution and sub-MJy/steradian sensitivity. Images taken with CAM filters LW2 (7 microns) and LW3 (15 microns) are analysed to determine the morphology of this galaxy and understand better the emission mechanisms. The mid-infrared emission follows an exponential disk with a scale length 75 arcsec. This is 60% of the scale length in the optical R-band and radio continuum. The nuclear starburst region is too bright for reliable measurement in these images. Its surface brightness exceeds the inner disk by at least a factor of 12. The arms and interarm regions are clearly detected, with each of these components contributing about equally to the disk emission. The arm-interarm contrast is 2-4 in the mid-IR, close to that measured in the visible R band light and lower than the contrast in $H\alpha$, suggesting that non-ionizing radiation contributes significantly to dust heating.