M33 is a gas rich spiral galaxy of the Local Group. We investigate the relationship between the two major gas cooling lines and the total infrared (TIR) dust continuum. We mapped the emission of gas and dust in M33 using the far-infrared lines of [CII] and [OI](63um) and the TIR. The line maps were observed with Herschel/PACS. These maps have 50pc resolution and form a ~370pc wide stripe along its major axis covering the sites of bright HII regions, but also more quiescent arm and inter-arm regions from the southern arm at 2kpc galacto-centric distance to the south out to 5.7kpc distance to the north. Full-galaxy maps of the continuum emission at 24um from Spitzer/MIPS, and at 70um, 100um, and 160um from PACS were combined to obtain a map of the TIR. TIR and [CII] intensities are correlated over more than two orders of magnitude. The range of TIR translates to a range of far ultraviolet (FUV) emission of G0,obs~2 to 200 in units of the average Galactic radiation field. The binned [CII]/TIR ratio drops with rising TIR, with large, but decreasing scatter. Fits of modified black bodies (MBBs) to the continuum emission were used to estimate dust mass surface densities and total gas column densities. A correction for possible foreground absorption by cold gas was applied to the [OI] data before comparing it with models of photon dominated regions (PDRs). Most of the ratios of [CII]/[OI] and ([CII]+[OI])/TIR are consistent with two model solutions. The median ratios are consistent with one solution at n~2x10^2 cm-3, G0~60, and and a second low-FUV solution at n~10^4 cm-3, G0~1.5. The bulk of the gas along the lines-of-sight is represented by a low-density, high-FUV phase with low beam filling factors ~1. A fraction of the gas may, however, be represented by the second solution.

The far-infrared (FIR) regime is one of the few wavelength ranges where no astronomical data with sub-arcsecond spatial resolution exist. Neither of the medium-term satellite projects like SPICA, Millimetron nor O.S.T. will resolve this malady. For many research areas, however, information at high spatial and spectral resolution in the FIR, taken from atomic fine-structure lines, from highly excited carbon monoxide (CO), light hydrids, and especially from water lines would open the door for transformative science. A main theme will be to trace the role of water in proto-planetary disks, to observationally advance our understanding of the planet formation process and, intimately related to that, the pathways to habitable planets and the emergence of life. Furthermore, key observations will zoom into the physics and chemistry of the star-formation process in our own Galaxy, as well as in external galaxies. The FIR provides unique tools to investigate in particular the energetics of heating, cooling and shocks. The velocity-resolved data in these tracers will reveal the detailed dynamics engrained in these processes in a spatially resolved fashion, and will deliver the perfect synergy with ground-based molecular line data for the colder dense gas.

The Origins Space Telescope (Origins) traces our cosmic history, from the formation of the first galaxies and the rise of metals to the development of habitable worlds and present-day life. Origins does this through exquisite sensitivity to infrared radiation from ions, atoms, molecules, dust, water vapor and ice, and observations of extra-solar planetary atmospheres, protoplanetary disks, and large-area extragalactic fields. Origins operates in the wavelength range 2.8 to 588 microns and is 1000 times more sensitive than its predecessors due to its large, cold (4.5 K) telescope and advanced instruments. Origins was one of four large missions studied by the community with support from NASA and industry in preparation for the 2020 Decadal Survey in Astrophysics. This is the final study report.

The astrophysical agenda of the 21st century requires a very sharp view of celestial objects. High angular resolution studies are essential for fundamental studies of a broad variety of astrophysical phenomena ranging from relativistic physics of black holes, their gravitational and electromagnetic imprints, violent transient processes, including those producing detectable gravitational waves, birth and evolution of planetary systems. Over the past decades, radio astronomy made huge leap in achieving ground-breaking angular resolution measured in tens of microarcseconds (one tenth of nanoradian and better). Recently a global Event Horizon Telescope (EHT) collaboration obtained first direct images of the shadow of a super-massive black hole in the nucleus of the active galaxy M87. These observations were conducted at 230 GHz. The two first generation Space Very Long Baseline Interferometry (VLBI) missions, VSOP/HALCA led by the Japan Aerospace Exploration Agency (JAXA) and RadioAstron led by the Russia Roscosmos State Corporation and Russia Academy of Sciences, achieved the highest angular resolution at frequencies from 0.3 to 22 GHz in observations conducted in the period 1997 - 2019. The next step in advancing high angular resolution radio astronomy is in combining high frequency (millimeter and sub-millimeter wavelengths) and interferometric baselines exceeding the Earth diameter. The present THEZA White Paper describes a combination which would unify technology developments in giga-/tera-hertz instrumentation and space-borne radio astronomy. The current preprint version of the THEZA White Paper is slightly re-formatted and edited comparing to the official submitted version.

(abridged) The Atacama Large Aperture Submillimeter Telescope (AtLAST) project aims to build a 50-m-class submm telescope with $>1^\circ$ field of view, high in the Atacama Desert, providing fast and detailed mapping of the mm/submm sky. It will thus serve as a strong complement to existing facilities such as ALMA. ALMA's small field of view ($<15^{\prime\prime}$ at 350 GHz) limits its mapping speed for large surveys. Instead, a single dish with a large field of view such as the AtLAST concept can host large multi-element instruments that can more efficiently map large portions of the sky. Small aperture survey instruments (typically much smaller than $<3\times$ the size of an interferometric array element) can mitigate this somewhat but lack the resolution for accurate recovery of source location and have small collecting areas. Furthermore, small aperture survey instruments do not provide sufficient overlap in the spatial scales they sample to provide a complete reconstruction of extended sources (i.e.\ the zero-spacing information is incomplete in $u,v$-space.) The heterodyne instrumentation for the AtLAST telescope that we consider here will take advantage of extensive developments in the past decade improving the performance and pixel count of heterodyne focal plane arrays. Such instrumentation, with higher pixel counts, has alredy begun to take advantage of integration in the focal planes to increase packaging efficiency over simply stacking modular mixer blocks in the focal plane. We extrapolate from the current state-of-the-art to present concept first-generation heterodyne designs for AtLAST.

It is usually assumed that the stellar initial mass function (IMF) takes a universal form and that there exists a direct mapping between this and the distribution of natal core masses (the core mass function, CMF). The IMF and CMF have been best characterized in the Solar neighborhood. Beyond 500~pc from the Sun, in diverse environments where metallicity varies and massive star feedback may dominate, the IMF has been measured only incompletely and imprecisely, while the CMF has hardly been measured at all. In order to establish if the IMF and CMF are indeed universal and related to each other, it is necessary to: 1) perform multi-wavelength large-scale imaging and spectroscopic surveys of different environments across the Galaxy; 2) require an angular resolution of < 0.1'' in the optical/near-IR for stars and < 5'' in the far-IR for cores; 3) achieve far-IR sensitivities to probe 0.1~Msol cores at 2--3 kpc.

Models of the chemical evolution of the interstellar medium, galaxies, and the Universe rely on our understanding of the amounts and chemical composition of the material returned by stars and supernovae. Stellar yields are obtained from stellar-evolution models, which currently lack predictive prescriptions of stellar mass loss, although it significantly affects stellar lifetimes, nucleosynthesis, and chemical ejecta. Galaxy properties are derived from observations of the integrated light of bright member stars. Stars in the late stages of their evolution are among the infrared-brightest objects in galaxies. An unrealistic treatment of the mass-loss process introduces significant uncertainties in galaxy properties derived from their integrated light. We describe current efforts and future needs and opportunities to characterize AGB outflows: driving mechanisms, outflow rates, underlying fundamental physical and chemical processes such as dust grain formation, and dependency of these on metallicity.

Water is fundamental to our understanding of the evolution of planetary systems and the delivery of volatiles to the surfaces of potentially habitable planets. Yet, we currently have essentially no facilities capable of observing this key species comprehensively. With this white paper, we argue that we need a relatively large, cold space-based observatory equipped with a high-resolution spectrometer, in the mid- through far-infrared wavelength range (25-600~$\mu$m) in order to answer basic questions about planet formation, such as where the Earth got its water, how giant planets and planetesimals grow, and whether water is generally available to planets forming in the habitable zone of their host stars.

The Origins Space Telescope, one of four large Mission Concept studies sponsored by NASA for review in the 2020 US Astrophysics Decadal Survey, will open unprecedented discovery space in the infrared, unveiling our cosmic origins. We briefly describe in this article the key science themes and architecture for OST. With a sensitivity gain of up to a factor of 1,000 over any previous or planned mission, OST will open unprecedented discovery space, allow us to peer through an infrared window teeming with possibility. OST will fundamentally change our understanding of our cosmic origins - from the growth of galaxies and black holes, to uncovering the trail of water, to life signs in nearby Earth-size planets, and discoveries never imagined. Built to be highly adaptable, while addressing key science across many areas of astrophysics, OST will usher in a new era of infrared astronomy.

We analyze 3 mm emission of the ultraluminous infrared galaxy Arp 220 for spatially-resolved structure and spectral properties of the merger nuclei. ALMA archival data at ~0.05" resolution are used for extensive visibility fitting and deep imaging of continuum emission. The data are fitted well with two concentric components for each nucleus, such as two Gaussians or one Gaussian plus one exponential disk. The larger components in individual nuclei are similar in shape and extent, ~100-150 pc, to the cm-wave emission due to supernovae. They are therefore identified with the known starburst nuclear disks. The smaller components in both nuclei have about a few 10 pc sizes and peak brightness temperatures (Tb) more than twice higher than in previous single-Gaussian fitting. They correspond to the dust emission that we find centrally concentrated in both nuclei by subtracting the plasma emission measured at 33 GHz. The dust emission in the western nucleus is found to have a peak Tb ~ 530 K and a full width at half maximum of about 20 pc. This component is estimated to have a bolometric luminosity on the order of 10^{12.5} Lsun and a 20 pc-scale luminosity surface density 10^{15.5} Lsun/kpc^2. A luminous AGN is a plausible energy source for these high values while other explanations remain to be explored. Our continuum image also reveals a third structural component of the western nucleus --- a pair of faint spurs perpendicular to the disk major axis. We attribute it to a bipolar outflow from the highly inclined (i ~ 60 deg) western nuclear disk.

[Abridged] Do some environments favor efficient conversion of molecular gas into stars? To answer this, we need to be able to estimate the H2 mass. Traditionally, this is done using CO and a few assumptions but the Herschel observations in the FIR make it possible to estimate the molecular gas mass independently of CO. Previous attempts to derive gas masses from dust emission suffered from biases. Generally, dust surface densities, HI column densities, and CO intensities are used to derive a gas-to-dust ratio (GDR) and the local CO intensity to H2 column density ratio (XCO), sometimes allowing for an additional CO-dark gas component (Kdark). We tested earlier methods, revealing degeneracies among the parameters, and then used a Bayesian formalism to derive the most likely values for each of the parameters mentioned above as a function of position in the nearby low metallicity spiral galaxy M33. The data are from the IRAM 30m CO(2-1) line, high-resolution HI and Herschel dust continuum observations. Solving for GDR, XCO, and Kdark in macro pixels 500 pc in size, we find that (i) allowing for CO-dark gas significantly improves fits; (ii) Kdark decreases with galactocentric distance; (iii) GDR is slightly higher than initially expected and increases with galactocentric distance; (iv) the total amount of dark gas closely follows the radially decreasing CO emission, as might be expected if the dark gas is H2 where CO is photodissociated. The total amount of H2, including dark gas, yields an average XCO of twice the galactic value of 2e20 cm^-2/(K km/s), 55% of this traced directly through CO. The rather constant fraction of dark gas suggests that there is no large population of diffuse H2 clouds (unrelated to GMCs) without CO emission. Unlike in large spirals, we detect no systematic radial trend in XCO, possibly linked to the absence of a radial decrease in CO line ratios.

We aim to understand the contribution of the ionized, atomic and molecular phases of the ISM to the [CII] emission from clouds near the dynamical center and the BCLMP302 HII region in the north of the nearby galaxy M33 at a spatial resolution of 50pc. We combine high resolution [CII] spectra taken with the HIFI spectrometer onboard the Herschel satellite with [CII] Herschel-PACS maps and ground-based observations of CO(2-1) and HI. All data are at a common spatial resolution of 50pc. Typically, the [CII] lines have widths intermediate between the narrower CO(2-1) and broader HI line profiles. We decomposed the [CII] spectra in terms of contribution from molecular and atomic gas detected in CO(2-1) and HI, respectively. We find that the relative contribution of molecular and atomic gas traced by CO(2-1) and HI varies depends mostly on the local physical conditions and geometry. We estimate that 11-60% and 5-34% of the [CII] intensities in the center and in BCLMP302, respectively, arise at velocities showing no CO(2-1) or HI emission and could arise in CO-dark molecular gas. The deduced strong variation in the [CII] emission not associated with CO and HI cannot be explained in terms of differences in A_v, far-ultraviolet radiation field, and metallicity between the two studied regions. Hence the relative amounts of diffuse (CO-dark) and dense molecular gas possibly vary on spatial scales smaller than 50pc. Based on the emission measure observed at radio wavelengths we estimate the contribution of ionized gas at a few positions to lie between 10-25%. The correlations between the intensities of tracers corresponding to the same velocity range as [CII], differ from the correlation derived from PACS data. The results in this paper emphasize the need for velocity-resolved observations to discern the contribution of different components of the ISM to [CII] emission. (abridged)

(Ultra) Luminous Infrared Galaxies ((U)LIRGs) are objects characterized by their extreme infrared (8-1000 $\mu$m) luminosities ($L_{LIRG}>10^{11} $L$_\odot$ and $L_{ULIRG}>10^{12}$ L$_\odot$). The Herschel Comprehensive ULIRG Emission Survey (HerCULES; PI van der Werf) presents a representative flux-limited sample of 29 (U)LIRGs that spans the full luminosity range of these objects (10$^{11}\leq L_\odot \geq10^{13}$). With the \emph{Herschel Space Observatory}, we observe [CII] 157 $\mu$m, [OI] 63 $\mu$m, and [OI] 145 $\mu$m line emission with PACS, CO J=4-3 through J=13-12, [CI] 370 $\mu$m, and [CI] 609 $\mu$m with SPIRE, and low-J CO transitions with ground-based telescopes. The CO ladders of the sample are separated into three classes based on their excitation level. In 13 of the galaxies, the [OI] 63 $\mu$m emission line is self absorbed. Comparing the CO excitation to the IRAS 60/100 $\mu$m ratio and to far infrared luminosity, we find that the CO excitation is more correlated to the far infrared colors. We present cooling budgets for the galaxies and find fine-structure line flux deficits in the [CII], [SiII], [OI], and [CI] lines in the objects with the highest far IR fluxes, but do not observe this for CO $4\leq J_{upp}\leq13$. In order to study the heating of the molecular gas, we present a combination of three diagnostic quantities to help determine the dominant heating source. Using the CO excitation, the CO J=1-0 linewidth, and the AGN contribution, we conclude that galaxies with large CO linewidths always have high-excitation CO ladders, and often low AGN contributions, suggesting that mechanical heating is important.

We study the wavelength dependence of the dust emission as a function of position and environment across the disk of M33 at a linear resolution of 160 pc using Spitzer and Herschel photometric data. Expressing the emissivity of the dust as a power law, the power-law exponent (beta) is estimated from two independent approaches designed to properly treat the degeneracy between beta and the dust temperature. Both beta and the dust temperature are higher in the inner disk than in the outer disk, contrary to reported beta-T anti-correlations found in other sources. In the cold + warm dust model, the warm component and the ionized gas (Halpha) have a very similar distribution across the galaxy, demonstrating that the model separates the components in an appropriate fashion. The flocculent spiral arms and the dust lanes are evident in the map of the cold component. Both cold and warm dust column densities are high in star forming regions and reach their maxima toward the giant star forming complexes NGC604 and NGC595. beta declines from close to 2 in the center to about 1.3 in the outer disk. beta is positively correlated with star formation and with molecular gas column, as traced by Halpha and CO emission. The lower dust emissivity index in the outer parts of M33 is likely related to the reduced metallicity (different grain composition) and possibly different size distribution. It is not due to the decrease in stellar radiation field or temperature in a simple way because the FIR-bright regions in the outer disk also have a low beta. Like most spirals, M33 has a (decreasing) radial gradient in star formation and molecular-to-atomic gas ratio such that the regions bright in Halpha or CO tend to trace the inner disk, making it difficult to distinguish between their effects on the dust.

We have observed the nucleus of the nearby luminous infrared galaxy NGC 4418 with subarcsec resolution at 860 and 450 micron for the first time to characterize its hidden power source. A ~20 pc (0.1") hot dusty core was found inside a 100 pc scale concentration of molecular gas at the galactic center. The 860 micron continuum core has a deconvolved (peak) brightness temperature of 120-210 K. The CO(3-2) peak brightness temperature there is as high as 90 K at 50 pc resolution. The core has a bolometric luminosity of about 10^{11} Lsun, which accounts for most of the galaxy luminosity. It is Compton thick (N_H >~ 10^{25} cm^{-2}) and has a high luminosity-to-mass ratio ~500 Lsun/Msun as well as a high luminosity surface density 10^{8.5+-0.5} Lsun pc^{-2}. These parameters are consistent with an AGN to be the main luminosity source (with an Eddington ratio about 0.3) while they can be also due to a young starburst near its maximum L/M. We also found an optical color (reddening) feature that we attribute to an outflow cone emanating from the nucleus. The hidden hot nucleus thus shows evidence of both an inflow, previously seen with absorption lines, and the new outflow reported here in a different direction. The nucleus must be rapidly evolving with these gas flows.

We present Herschel SPIRE FTS spectroscopy of the nearby luminous infrared galaxy NGC 6240. In total 20 lines are detected, including CO J=4-3 through J=13-12, 6 H2O rotational lines, and [CI] and [NII] fine-structure lines. The CO to continuum luminosity ratio is 10 times higher in NGC 6240 than Mrk 231. Although the CO ladders of NGC 6240 and Mrk 231 are very similar, UV and/or X-ray irradiation are unlikely to be responsible for the excitation of the gas in NGC 6240. We applied both C and J shock models to the H2 v=1-0 S(1) and v=2-1 S(1) lines and the CO rotational ladder. The CO ladder is best reproduced by a model with shock velocity v_s=10 km s^-1 and a pre-shock density n_H=5 * 10^4 cm^-3. We find that the solution best fitting the H2 lines is degenerate: The shock velocities and number densities range between v_s = 17 - 47 km s^-1 and n_H=10^7 - 5 * 10^4 cm^-3, respectively. The H2 lines thus need a much more powerful shock than the CO lines. We deduce that most of the gas is currently moderately stirred up by slow (10 km s^-1) shocks while only a small fraction (< 1 percent) of the ISM is exposed to the high velocity shocks. This implies that the gas is rapidly loosing its highly turbulent motions. We argue that a high CO line-to-continuum ratio is a key diagnostic for the presence of shocks.

We report 350 and 230 GHz observations of molecular gas and dust in the starburst nucleus of NGC 253 at 20-40 pc (1"-2") resolution. The data contain CO(3-2), HCN(4-3), CO(2-1), 13CO(2-1), C18O(2-1), and continuum at 0.87 mm and 1.3 mm toward the central kiloparsec. The CO(2-1) size of the galaxy's central molecular zone (CMZ) is measured to be about 300 pc x 100 pc at the half maximum of intensity. Five clumps of dense and warm gas stand out in the CMZ at arcsecond resolution, and they are associated with compact radio sources due to recent massive star formation. They contribute one third of the CO emission in the central 300 pc and have 12CO peak brightness temperatures around 50 K, molecular gas column densities on the order of 10^{4} Msun pc^{-2}, gas masses on the order of 10^{7} Msun in the size scale of 20 pc, volume-averaged gas densities of n(H2) ~ 4000 cm^{-3}, and high HCN-to-CO ratios suggestive of higher fractions of dense gas than in the surrounding environment. It is suggested that these are natal molecular cloud complexes of massive star formation. The CMZ of NGC 253 is also compared with that of our Galaxy in CO(2-1) at the same 20 pc resolution. Their overall gas distributions are strikingly similar. The five molecular cloud complexes appear to be akin to such molecular complexes as Sgr A, Sgr B2, Sgr C, and the l=1.3deg cloud in the Galactic center. On the other hand, the starburst CMZ in NGC 253 has higher temperatures and higher surface (and presumably volume) densities than its non-starburst cousin.

Infrared pumping and its effect on the excitation of HCN molecules can be important when using rotational lines of HCN to probe dense molecular gas in galaxy nuclei. We report the first extragalactic detection of (sub)millimeter rotational lines of vibrationally excited HCN, in the dust-enshrouded nucleus of the luminous infrared galaxy NGC 4418. We estimate the excitation temperature of T_vib ~ 230 K between the vibrational ground and excited (v_2=1) states. This excitation is most likely due to infrared radiation. At this high vibrational temperature the path through the v_2=1 state must have a strong impact on the rotational excitation in the vibrational ground level, although it may not be dominant for all rotational levels. Our observations also revealed nearly confusion limited lines of CO, HCN, HCO+, H13CN, HC15N, CS, N2H+, and HC3N at lambda ~ 1 mm. Their relative intensities may also be affected by the infrared pumping.

We present high resolution large scale observations of the molecular and atomic gas in the Local Group Galaxy M33. The observations were carried out using the HERA at the 30m IRAM telescope in the CO(2-1) line achieving a resolution of 12"x2.6 km/s, enabling individual GMCs to be resolved. The observed region mainly along the major axis out to a radius of 8.5 kpc, and covers the strip observed with HIFI/PACS Spectrometers as part of the HERM33ES Herschel key program. The achieved sensitivity in main beam temperature is 20-50 mK at 2.6 km/s velocity resolution. The CO(2-1) luminosity of the observed region is 1.7\pm0.1x10^7 Kkm/s pc^2, corresponding to H2 masses of 1.9x10^8 Msun (including He), calculated with a NH2/ICO twice the Galactic value due to the half-solar metallicity of M33. HI 21 cm VLA archive observations were reduced and the mosaic was imaged and cleaned using the multi-scale task in CASA, yielding a series of datacubes with resolutions ranging from 5" to 25". The HI mass within a radius of 8.5 kpc is estimated to be 1.4x10^9 Msun. The azimuthally averaged CO surface brightness decreases exponentially with a scale length of 1.9\pm0.1 kpc whereas the atomic gas surface density is constant at Sigma_HI=6\pm2 Msun/pc^2 deprojected to face-on. The central kiloparsec H_2 surface density is Sigma_H2=8.5\pm0.2 Msun/pc^2. The star formation rate per unit molecular gas (SF Efficiency, the rate of transformation of molecular gas into stars), as traced by the ratio of CO to Halpha and FIR brightness, is constant with radius. The SFE appears 2-4 times greater than of large spiral galaxies. A morphological comparison of molecular and atomic gas with tracers of star formation shows good agreement between these maps both in terms of peaks and holes. A few exceptions are noted. Several spectra, including those of a molecular cloud situated more than 8 kpc from the galaxy center, are presented.

The Ultra Luminous InfraRed Galaxy Mrk 231 reveals up to seven rotational lines of water (H2O) in emission, including a very high-lying (E_{upper}=640 K) line detected at a 4sigma level, within the Herschel/SPIRE wavelength range, whereas PACS observations show one H2O line at 78 microns in absorption, as found for other H2O lines previously detected by ISO. The absorption/emission dichotomy is caused by the pumping of the rotational levels by far-infrared radiation emitted by dust, and subsequent relaxation through lines at longer wavelengths, which allows us to estimate both the column density of H2O and the general characteristics of the underlying far-infrared continuum source. Radiative transfer models including excitation through both absorption of far-infrared radiation emitted by dust and collisions are used to calculate the equilibrium level populations of H2O and the corresponding line fluxes. The highest-lying H2O lines detected in emission, with levels at 300-640 K above the ground state, indicate that the source of far-infrared radiation responsible for the pumping is compact (radius=110-180 pc) and warm (T_{dust}=85-95 K), accounting for at least 45% of the bolometric luminosity. The high column density, N(H2O)~5x10^{17} cm^{-2}, found in this nuclear component, is most probably the consequence of shocks/cosmic rays, an XDR chemistry, and/or an "undepleted chemistry" where grain mantles are evaporated. A more extended region, presumably the inner region of the 1-kpc disk observed in other molecular species, could contribute to the flux observed in low-lying H2O lines through dense hot cores, and/or shocks. The H2O 78 micron line observed with PACS shows hints of a blue-shifted wing seen in absorption, possibly indicating the occurrence of H2O in the prominent outflow detected in OH (Fischer et al., this volume).

We present an analysis of the first space-based far-IR-submm observations of M 33, which measure the emission from the cool dust and resolve the giant molecular cloud complexes. With roughly half-solar abundances, M33 is a first step towards young low-metallicity galaxies where the submm may be able to provide an alternative to CO mapping to measure their H$_2$ content. In this Letter, we measure the dust emission cross-section $\sigma$ using SPIRE and recent CO and \HI\ observations; a variation in $\sigma$ is present from a near-solar neighborhood cross-section to about half-solar with the maximum being south of the nucleus. Calculating the total H column density from the measured dust temperature and cross-section, and then subtracting the \HI\ column, yields a morphology similar to that observed in CO. The H$_2$/\HI\ mass ratio decreases from about unity to well below 10% and is about 15% averaged over the optical disk. The single most important observation to reduce the potentially large systematic errors is to complete the CO mapping of M 33.

We present a full high resolution SPIRE FTS spectrum of the nearby ultraluminous infrared galaxy Mrk231. In total 25 lines are detected, including CO J=5-4 through J=13-12, 7 rotational lines of H2O, 3 of OH+ and one line each of H2O+, CH+, and HF. We find that the excitation of the CO rotational levels up to J=8 can be accounted for by UV radiation from star formation. However, the approximately flat luminosity distribution of the CO lines over the rotational ladder above J=8 requires the presence of a separate source of excitation for the highest CO lines. We explore X-ray heating by the accreting supermassive black hole in Mrk231 as a source of excitation for these lines, and find that it can reproduce the observed luminosities. We also consider a model with dense gas in a strong UV radiation field to produce the highest CO lines, but find that this model strongly overpredicts the hot dust mass in Mrk231. Our favoured model consists of a star forming disk of radius 560 pc, containing clumps of dense gas exposed to strong UV radiation, dominating the emission of CO lines up to J=8. X-rays from the accreting supermassive black hole in Mrk231 dominate the excitation and chemistry of the inner disk out to a radius of 160 pc, consistent with the X-ray power of the AGN in Mrk231. The extraordinary luminosity of the OH+ and H2O+ lines reveals the signature of X-ray driven excitation and chemistry in this region.

Within the framework of the HERM33ES key project, we are studying the star forming interstellar medium in the nearby, metal-poor spiral galaxy M33, exploiting the high resolution and sensitivity of Herschel. We use PACS and SPIRE maps at 100, 160, 250, 350, and 500 micron wavelength, to study the variation of the spectral energy distributions (SEDs) with galacto-centric distance. Detailed SED modeling is performed using azimuthally averaged fluxes in elliptical rings of 2 kpc width, out to 8 kpc galacto-centric distance. Simple isothermal and two-component grey body models, with fixed dust emissivity index, are fitted to the SEDs between 24 and 500 micron using also MIPS/Spitzer data, to derive first estimates of the dust physical conditions. The far-infrared and submillimeter maps reveal the branched, knotted spiral structure of M33. An underlying diffuse disk is seen in all SPIRE maps (250-500 micron). Two component fits to the SEDs agree better than isothermal models with the observed, total and radially averaged flux densities. The two component model, with beta fixed at 1.5, best fits the global and the radial SEDs. The cold dust component clearly dominates; the relative mass of the warm component is less than 0.3% for all the fits. The temperature of the warm component is not well constrained and is found to be about 60K plus/minus 10K. The temperature of the cold component drops significantly from about 24K in the inner 2 kpc radius to 13K beyond 6 kpc radial distance, for the best fitting model. The gas-to-dust ratio for beta=1.5, averaged over the galaxy, is higher than the solar value by a factor of 1.5 and is roughly in agreement with the subsolar metallicity of M33.

We report ~100 pc (0.3") resolution observations of (sub)millimeter HCO+ and CO lines in the ultraluminous infrared galaxy Arp 220. The lines peak at two merger nuclei, with HCO+ being more spatially concentrated than CO. Asymmetric line profiles with blueshifted absorption and redshifted emission are discovered in HCO+(3-2) and (4-3) toward the two nuclei and in CO(3-2) toward one nucleus. We suggest that these P-Cygni profiles are due to ~100 km/s outward motion of molecular gas from the nuclei. This gas is most likely outflowing from the inner regions of the two nuclear disks rotating around individual nuclei, clearing the shroud around the luminosity sources there.

NGC 2024, a sites of massive star formation, have complex internal structures caused by cal heating by young stars, outflows, and stellar winds. These complex cloud structures lead to intricate emission line shapes. The goal of this paper is to show that the complex line shapes of 12 CO lines in NGC 2024 can be explained consistently with a model, whose temperature and velocity structure are based on the well-established scenario of a PDR and the Blister model. We present velocity-resolved spectra of seven CO lines ranging from J=3 to J=13, and we combined these data with CO high-frequency data from the ISO satellite. We find that the bulk of the molecular cloud associated with NGC 2024 consists of warm (75 K) and dense (9e5 cm-3) gas. An additional hot (~ 300 K) component, located at the interface of the HII region and the molecular cloud, is needed to explain the emission of the high-J CO lines. Deep absorption notches indicate that very cold material (20 K) exists in front of the warm material, too. A temperature and column density structure consistent with those predicted by PDR models, combined with the velocity structure of a Blister model, appropriately describes the observed emission line profiles of this massive star forming region. This case study of NGC 2024 shows that, with physical insights into these complex regions and careful modeling, multi-line observations of CO can be used to derive detailed physical conditions in massive star forming regions.

We have used the Submillimeter Array (SMA) to make the first interferometric observations (beam size ~1") of the 12CO J=6-5 line and 435 micron (690 GHz) continuum emission toward the central region of the nearby ULIRG Arp 220. These observations resolve the eastern and western nuclei from each other, in both the molecular line and dust continuum emission. At 435 micron, the peak intensity of the western nucleus is stronger than the eastern nucleus, and the difference in peak intensities is less than at longer wavelengths. Fitting a simple model to the dust emission observed between 1.3 mm and 435 micron suggests that dust emissivity power law index in the western nucleus is near unity and steeper in the eastern nucleus, about 2, and that the dust emission is optically thick at the shorter wavelength. Comparison with single dish measurements indicate that the interferometer observations are missing ~60% of the dust emission, most likely from a spatially extended component to which these observations are not sensitive. The 12CO J=6-5 line observations clearly resolve kinematically the two nuclei. The distribution and kinematics of the 12CO J=6-5 line appear to be very similar to lower J CO lies observed at similar resolution. Analysis of multiple 12CO line intensities indicates that the molecular gas in both nuclei have similar excitation conditions, although the western nucleus is warmer and denser. The excitation conditions are similar to those found in other extreme environments, including M82, Mrk 231, and BR 1202-0725. Simultaneous lower resolution observations of the 12CO, 13CO, and C18O J=2-1 lines show that the 13CO and C18O lines have similar intensities, which suggests that both of these lines are optically thick, or possibly that extreme high mass star formation has produced in an overabundance of C18O.

We present high resolution (0."4) observations of HNC J=3-2 with the SubMillimeter Array (SMA). We find luminous HNC 3-2 line emission in the western part of Arp220, centered on the western nucleus, while the eastern side of the merger shows relatively faint emission. A bright (36 K), narrow (60 km/s) emission feature emerges from the western nucleus, superposed on a broader spectral component. A possible explanation is weak maser emission through line-of-sight amplification of the background continuum source. There is also a more extended HNC 3-2 emission feature north and south of the nucleus. This feature resembles the bipolar OH maser morphology around the western nucleus. Substantial HNC abundances are required to explain the bright line emission from this warm environment. We discuss this briefly in the context of an X-ray chemistry and radiative excitation. We conclude that the luminous and possibly amplified HNC emission of the western nucleus of the Arp220 merger reflects the unusual, and perhaps transient, environment of the starburst/AGN activity there. The faint HNC line emission towards Arp220-east reveals a real difference in physical conditions between the two merger nuclei.

Deuterated ions are abundant in cold (T=10 K), dense (n=10^5 cm^-3) regions, in which CO is frozen out onto dust grains. In such environments, the deuterium fractionation of such ions can exceed the elemental abundance ratio of D/H by a factor of 10^4. In this paper we use the deuterium fractionation to investigate the evolutionary state of Class 0 protostars. In a sample of 20 protostellar objects, we found a clear correlation between the N2D+/N2H+ ratio and evolutionary tracers. As expected, the coolest, i.e. the youngest, objects show the largest deuterium fractionation. Furthermore, we find that sources with a high N2D+/N2H+ ratio show clear indication for infall.

We used the Submillimeter Array (SMA) to image 860 micron continuum and CO(3-2) line emission in the ultraluminous merging galaxy Arp 220, achieving a resolution of 0.23" (80 pc) for the continuum and 0.33" (120 pc) for the line. The CO emission peaks around the two merger nuclei with a velocity signature of gas rotation around each nucleus, and is also detected in a kpc-size disk encompassing the binary nucleus. The dust continuum, in contrast, is mostly from the two nuclei. The beam-averaged brightness temperature of both line and continuum emission exceeds 50 K at and around the nuclei, revealing the presence of warm molecular gas and dust. The dust emission morphologically agrees with the distribution of radio supernova features in the east nucleus, as expected when a starburst heats the nucleus. In the brighter west nucleus, however, the submillimeter dust emission is more compact than the supernova distribution. The 860 micron core, after deconvolution, has a size of 50-80 pc, consistent with recent 1.3 mm observations, and a peak brightness temperature of (0.9-1.6)x10^2 K. Its bolometric luminosity is at least 2x10^{11} Lsun and could be ~10^{12} Lsun depending on source structure and 860 micron opacity, which we estimate to be of the order of tau_{860} ~ 1 (i.e., N_{H_2} ~ 10^{25} cm^{-2}). The starbursting west nuclear disk must have in its center a dust enshrouded AGN or a very young starburst equivalent to hundreds of super star clusters. Further spatial mapping of bolometric luminosity through submillimeter imaging is a promising way to identify the heavily obscured heating sources in Arp 220 and other luminous infrared galaxies.

We find that the HNC J=3-2 emission is brighter than the HCN 3-2 emission by factors of 1.5 to 2.3 in the ultraluminous mergers Arp220 and Mrk231 and the luminous IR galaxy NGC4418. We furthermore report the detection of HNC J=4-3 in Mrk231. Overluminous HNC emission is unexpected in warm molecular gas in ultraluminous galaxies since I(HNC)>I(HCN) is usually taken as a signature of cold (10 - 20 K) dark clouds. Since the molecular gas of the studied galaxies is warm (T_k > 40 K) we present two alternative explanations to the overluminous HNC: a) HNC excitation is affected by pumping of the rotational levels through the mid-infrared continuum and b) XDRs (X-ray Dominated Regions) influence the abundances of HNC. HNC may become pumped at 21.5 micron brightness temperatures of 50 K, suggesting that HNC-pumping could be common in warm, ultraluminous galaxies with compact IR-nuclei.On the other hand, all three galaxies are either suspected of having buried AGN - or the presence of AGN is clear (Mrk231) - indicating that X-rays may affect the ISM chemistry. We conclude that both the pumping and XDR alternatives imply molecular cloud ensembles distinctly different from those of typical starforming regions in the Galaxy, or the ISM of less extreme starburst galaxies. The HNC molecule shows the potential of becoming an additional important tracer of extreme nuclear environments.

The THz atmospheric windows centered at roughly 1.3 and 1.5~THz, contain numerous spectral lines of astronomical importance, including three high-J CO lines, the N+ line at 205 microns, and the ground transition of para-H2D+. The CO lines are tracers of hot (several 100K), dense gas; N+ is a cooling line of diffuse, ionized gas; the H2D+ line is a non-depleting tracer of cold (~20K), dense gas. As the THz lines benefit the study of diverse phenomena (from high-mass star-forming regions to the WIM to cold prestellar cores), we have built the CO N+ Deuterium Observations Receiver (CONDOR) to further explore the THz windows by ground-based observations. CONDOR was designed to be used at the Atacama Pathfinder EXperiment (APEX) and Stratospheric Observatory For Infrared Astronomy (SOFIA). CONDOR was installed at the APEX telescope and test observations were made to characterize the instrument. The combination of CONDOR on APEX successfully detected THz radiation from astronomical sources. CONDOR operated with typical Trec=1600K and spectral Allan variance times of 30s. CONDOR's first light observations of CO 13-12 emission from the hot core Orion FIR4 (= OMC1 South) revealed a narrow line with T(MB) = 210K and delta(V)=5.4km/s. A search for N+ emission from the ionization front of the Orion Bar resulted in a non-detection. The successful deployment of CONDOR at APEX demonstrates the potential for making observations at THz frequencies from ground-based facilities.

The presence of H2D+ in dense cloud cores underlies ion-molecule reactions that strongly enhance the deuterium fractionation of many molecular species. We determine the H2D+ abundance in one starless core, Barnard 68, that has a particularly well established physical, chemical, and dynamical structure. We observed the ortho-H2D+ ground-state line 1_10-1_11, the N2H+ J=4-3 line, and the H13CO+ 4-3 line with the APEX telescope. We report the probable detection of the o-H2D+ line at an intensity Tmb=0.22 +- 0.08 K and exclusively thermal line width, and find only upper limits to the N2H+ 4-3 and H13CO+ 4-3 intensities. Within the uncertainties in the chemical reaction rates and the collisional excitation rates, chemical model calculations and excitation simulations reproduce the observed intensities and that of o-H2D+ in particular.

Prestellar cores are unique laboratories for studies of the chemical and physical conditions preceding star formation. We observed the prestellar core L1544 in the fundamental transition of ortho-H2D+ (1_1,0-1_1,1) at different positions over 100", and found a strong correlation between its abundance and the CO depletion factor. We also present a tentative detection of the fundamental transition of para-D2H+ (1_1,0-1_0,1) at the dust emission peak. Maps in N2H+, N2D+, HCO+ and DCO+ are used, and interpreted with the aid of a spherically symmetric chemical model that predicts the column densities and abundances of these species as a function of radius. The correlation between the observed deuterium fractionation of H3+, N2H+ and HCO+ and the observed integrated CO depletion factor across the core can be reproduced by this chemical model. In addition a simpler model is used to study the H2D+ ortho-to-para ratio. We conclude that, in order to reproduce the observed ortho-H2D+ observations, the grain radius should be larger than 0.3 microns.

The central 2x1 kpc of the starburst galaxy NGC 253 has been imaged using the Submillimeter Array at a 60 pc resolution in the J=2-1 transitions of 12CO, 13CO, and C18O as well as in the 1.3 mm continuum. Molecular gas and dust are mainly in the circumnuclear disk of ~500 pc radius, with warm (~40 K) and high area-filling factor gas in its central part. Two gas shells or cavities have been discovered in the circumnuclear disk. They have ~100 pc diameters and have large velocity widths of 80-100 km/s, suggestive of expansion at ~50 km/s. Modeled as an expanding bubble, each shell has an age of ~0.5 Myr and needed kinetic energy of ~1E46 J as well as mean mechanical luminosity of ~1E33 W for its formation. The large energy allows each to be called a superbubble. A ~10^6 Msun super star cluster can provide the luminosity, and could be a building block of the nuclear starburst in NGC 253. Alternatively, a hypernova can also be the main source of energy for each superbubble, not only because it can provide the mechanical energy and luminosity but also because the estimated rate of superbubble formation and that of hypernova explosions are comparable. Our observations indicate that the circumnuclear molecular disk harboring the starburst is highly disturbed on 100 pc or smaller scales, presumably by individual young clusters and stellar explosions, in addition to globally disturbed in the form of the well-known superwind.

We present the first interferometric CO(J=3-2) observations (beam size of 3.9"x1.6" or 160pc x 65pc) with the Submillimeter Array (SMA) toward the center of the Seyfert 2 galaxy M51. The image shows a strong concentration at the nucleus and weak emission from the spiral arm to the northwest. The integrated intensity of the central component in CO(J=3-2) is almost twice as high as that in CO(J=1-0), indicating that the molecular gas within an ~80 pc radius of the nucleus is warm (>~100 K) and dense (~10^4 cm^-3). Similar intensity ratios are seen in shocked regions in our Galaxy, suggesting that these gas properties may be related to AGN or starburst activity. The central component shows a linear velocity gradient (~1.4 km/s/pc) perpendicular to the radio continuum jet, similar to that seen in previous observations and interpreted as a circumnuclear molecular disk/torus around the Seyfert 2 nucleus. In addition, we identify a linear velocity gradient (~0.7 km/s/pc) along the jet. Judging from the energetics, the velocity gradient can be explained by supernova explosions or energy and momentum transfer from the jet to the molecular gas via interaction, which is consistent with the high intensity ratio.

The Submillimeter Array (SMA) has been under construction at the 4100 m summit of Mauna Kea, Hawaii. The array is going to allow imaging of lines and continuum at sub-arcsecond resolution in submillimeter wavelengths. The status of the array and the results from recent commissioning observations of nearby galaxies are reported.

The center of M83, a barred starburst galaxy with a double nucleus, has been observed in the CO(J=2-1) and CO(J=3-2) lines with the Submillimeter Array. The molecular gas shows a distribution and kinematics typical for barred galaxies at kpc radii, but reveals unusual kinematics around the double nucleus in the central 300 pc. Our CO velocity data show that the visible nucleus in M83 is at least 3" (65 pc) away from the galaxy's dynamical center, which most likely coincides with the center of symmetry previously determined in K band and is suggested to host another nucleus. We discovered high-velocity molecular gas associated with the visible off-center nucleus, and also found a steep velocity gradient across it. We attribute these features to a gas disk rotating around the off-center nucleus, which may be the remnant of a small galaxy accreted by M83. The dynamical mass of this component is estimated to be 3x10^8 Msun within a radius of 40 pc. The dynamical perturbation from the off-center nucleus may have played a key role in shaping the lopsided starburst.

We report the first submillimeter interferometric observations of an ultraluminous infrared galaxy. We observed Arp220 in the CO J=3-2 line and 342GHz continuum with the single baseline CSO-JCMT interferometer consisting of the Caltech Submillimeter Observatory (CSO) and the James Clerk Maxwell Telescope (JCMT). Models were fit to the measured visibilities to constrain the structure of the source. The morphologies of the CO J=3-2 line and 342GHz continuum emission are similar to those seen in published maps at 230 and 110GHz. We clearly detect a binary source separated by about 1 arcsec in the east-west direction in the 342GHz continuum. The CO J=3-2 visibility amplitudes, however, indicate a more complicated structure, with evidence for a compact binary at some velocities and rather more extended structure at others. Less than 30% of the total CO J=3-2 emission is detected by the interferometer, which implies the presence of significant quantities of extended gas. We also obtained single-dish CO J=2-1, CO J=3-2 and HCN J=4-3 spectra. The HCN J=4-3 spectrum, unlike the CO spectra, is dominated by a single redshifted peak. The HCN J=4-3/CO J=3-2, HCN J=4-3/HCN J=1-0 and CO J=3-2/2-1 line ratios are larger in the redshifted (eastern) source, which suggests that the two sources may have different physical conditions. This result might be explained by the presence of an intense starburst that has begun to deplete or disperse the densest gas in the western source, while the eastern source harbors undispersed high density gas.