We present \textit{Spitzer} full-orbit thermal phase curves of the hot Jupiter Qatar-1b, a planet with the same equilibrium temperature---and intermediate surface gravity and orbital period---as the well-studied planets HD 209458b and WASP-43b. We measure secondary eclipse of $0.21 \pm 0.02 \%$ at $3.6~\mu$m and $0.30 \pm 0.02 \%$ at $4.5~\mu$m, corresponding to dayside brightness temperatures of $1542^{+32}_{-31}$~K and $1557^{+35}_{-36}$~K, respectively, consistent with a vertically isothermal dayside. The respective nightside brightness temperatures are $1117^{+76}_{-71}$~K and $1167^{+69}_{-74}$~K, in line with a trend that hot Jupiters all have similar nightside temperatures. We infer a Bond albedo of $0.12_{-0.16}^{+0.14}$ and a moderate day-night heat recirculation efficiency, similar to HD 209458b. General circulation models for HD 209458b and WASP-43b predict that their bright-spots should be shifted east of the substellar point by tens of degrees, and these predictions were previously confirmed with \textit{Spitzer} full-orbit phase curve observations. The phase curves of Qatar-1b are likewise expected to exhibit eastward offsets. Instead, the observed phase curves are consistent with no offset: $11^{\circ}\pm 7^{\circ}$ at $3.6~\mu$m and $-4^{\circ}\pm 7^{\circ}$ at $4.5~\mu$m. The discrepancy in circulation patterns between these three otherwise similar planets points to the importance of secondary parameters like rotation rate and surface gravity, and the presence or absence of clouds, in determining atmospheric conditions on hot Jupiters.

Phase curve observations provide an opportunity to study the full energy budgets of exoplanets by quantifying the amount of heat redistributed from their daysides to their nightsides. Theories explaining the properties of phase curves for hot Jupiters have focused on the balance between radiation and dynamics as the primary parameter controlling heat redistribution. However, recent phase curves have shown deviations from the trends that emerge from this theory, which has led to work on additional processes that may affect hot Jupiter energy budgets. One such process, molecular hydrogen dissociation and recombination, can enhance energy redistribution on ultra-hot Jupiters with temperatures above $\sim2000$ K. In order to study the impact of H$_{2}$ dissociation on ultra-hot Jupiters, we present a phase curve of KELT-9b observed with the Spitzer Space Telescope at 4.5 $\mu$m. KELT-9b is the hottest known transiting planet, with a 4.5-$\mu$m dayside brightness temperature of $4566^{+140}_{-136}$ K and a nightside temperature of $2556^{+101}_{-97}$ K. We observe a phase curve amplitude of $0.609 \pm 0.020$ and a hot spot offset of $18.7^{+2.1}_{-2.3}$ degrees. The observed amplitude is too small to be explained by a simple balance between radiation and advection. General circulation models (GCMs) and an energy balance model that include the effects of H$_{2}$ dissociation and recombination provide a better match to the data. The GCMs, however, predict a maximum hot spot offset of $5$ degrees, which disagrees with our observations at $>5\sigma$ confidence. This discrepancy may be due to magnetic effects in the planet's highly ionized atmosphere.

The exoplanet WASP-12b is the prototype for the emerging class of ultra-hot, Jupiter-mass exoplanets. Past models have predicted---and near ultra-violet observations have shown---that this planet is losing mass. We present an analysis of two sets of 3.6 $\mu$m and 4.5 $\mu$m $\textit{Spitzer}$ phase curve observations of the system which show clear evidence of infrared radiation from gas stripped from the planet, and the gas appears to be flowing directly toward or away from the host star. This accretion signature is only seen at 4.5 $\mu$m, not at 3.6 $\mu$m, which is indicative either of CO emission at the longer wavelength or blackbody emission from cool, $\lesssim$ 600 K gas. It is unclear why WASP-12b is the only ultra-hot Jupiter to exhibit this mass loss signature, but perhaps WASP-12b's orbit is decaying as some have claimed, while the orbits of other exoplanets may be more stable; alternatively, the high energy irradiation from WASP-12A may be stronger than the other host stars. We also find evidence for phase offset variability at the level of $6.4\sigma$ ($46.2^{\circ}$) at 3.6 $\mu$m.

A new class of exoplanets is beginning to emerge: planets whose dayside atmospheres more closely resemble stellar atmospheres as most of their molecular constituents dissociate. The effects of the dissociation of these species will be varied and must be carefully accounted for. Here we take the first steps towards understanding the consequences of dissociation and recombination of molecular hydrogen (H$_2$) on atmospheric heat recirculation. Using a simple energy balance model with eastward winds, we demonstrate that H$_2$ dissociation/recombination can significantly increase the day$-$night heat transport on ultra-hot Jupiters (UHJs): gas giant exoplanets where significant H$_2$ dissociation occurs. The atomic hydrogen from the highly irradiated daysides of UHJs will transport some of the energy deposited on the dayside towards the nightside of the planet where the H atoms recombine into H$_2$; this mechanism bears similarities to latent heat. Given a fixed wind speed, this will act to increase the heat recirculation efficiency; alternatively, a measured heat recirculation efficiency will require slower wind speeds after accounting for H$_2$ dissociation/recombination.

We present an optical eclipse observation of the hot Jupiter WASP-12b using the Space Telescope Imaging Spectrograph on board the Hubble Space Telescope. These spectra allow us to place an upper limit of $A_g < 0.064$ (97.5% confidence level) on the planet's white light geometric albedo across 290--570 nm. Using six wavelength bins across the same wavelength range also produces stringent limits on the geometric albedo for all bins. However, our uncertainties in eclipse depth are $\sim$40% greater than the Poisson limit and may be limited by the intrinsic variability of the Sun-like host star --- the solar luminosity is known to vary at the $10^{-4}$ level on a timescale of minutes. We use our eclipse depth limits to test two previously suggested atmospheric models for this planet: Mie scattering from an aluminum-oxide haze or cloud-free Rayleigh scattering. Our stringent nondetection rules out both models and is consistent with thermal emission plus weak Rayleigh scattering from atomic hydrogen and helium. Our results are in stark contrast with those for the much cooler HD 189733b, the only other hot Jupiter with spectrally resolved reflected light observations; those data showed an increase in albedo with decreasing wavelength. The fact that the first two exoplanets with optical albedo spectra exhibit significant differences demonstrates the importance of spectrally resolved reflected light observations and highlights the great diversity among hot Jupiters.

We measure the mass function for a sample of 840 young star clusters with ages between 10-300 Myr observed by the Panchromatic Hubble Andromeda Treasury (PHAT) survey in M31. The data show clear evidence of a high-mass truncation: only 15 clusters more massive than $10^4$ $M_{\odot}$ are observed, compared to $\sim$100 expected for a canonical $M^{-2}$ pure power-law mass function with the same total number of clusters above the catalog completeness limit. Adopting a Schechter function parameterization, we fit a characteristic truncation mass of $M_c = 8.5^{+2.8}_{-1.8} \times 10^3$ $M_{\odot}$. While previous studies have measured cluster mass function truncations, the characteristic truncation mass we measure is the lowest ever reported. Combining this M31 measurement with previous results, we find that the cluster mass function truncation correlates strongly with the characteristic star formation rate surface density of the host galaxy, where $M_c \propto$ $\langle \Sigma_{\mathrm{SFR}} \rangle^{\sim1.1}$. We also find evidence that suggests the observed $M_c$-$\Sigma_{\mathrm{SFR}}$ relation also applies to globular clusters, linking the two populations via a common formation pathway. If so, globular cluster mass functions could be useful tools for constraining the star formation properties of their progenitor host galaxies in the early Universe.

The interaction of ionizing and far-ultraviolet radiation with the interstellar medium is of great importance. It results in the formation of regions in which the gas is ionized, beyond which are photodissociation regions (PDRs) in which the gas transitions to its atomic and molecular form. Several numerical codes have been implemented to study these two main phases of the interstellar medium either dynamically or chemically. In this paper we present TORUS-3DPDR, a new self-consistent code for treating the chemistry of three-dimensional photoionization and photodissociation regions. It is an integrated code coupling the two codes TORUS, a hydrodynamics and Monte Carlo radiation transport code, and 3D-PDR, a photodissociation regions code. The new code uses a Monte Carlo radiative transfer scheme to account for the propagation of the ionizing radiation including the diffusive component as well as a ray-tracing scheme based on the HEALPix package in order to account for the escape probability and column density calculations. Here, we present the numerical techniques we followed and we show the capabilities of the new code in modelling three-dimensional objects including single or multiple sources. We discuss the effects introduced by the diffusive component of the UV field in determining the thermal balance of PDRs as well as the effects introduced by a multiple sources treatment of the radiation field. We find that diffuse radiation can positively contribute to the formation of CO. With this new code, three-dimensional synthetic observations for the major cooling lines are possible, for making feasible a detailed comparison between hydrodynamical simulations and observations.

We present a model of the photodissociation regions of NGC 4038, which is part of the Antennae galaxies. We have considered one-dimensional slabs of uniform density all having a maximum $A_V=10\,{\rm mag}$, interacting with plane-parallel radiation. The density range in our simulations spans four orders of magnitude ($100\le n\le 10^6\,{\rm cm}^{-3}$) and the UV field strength spans more than three orders of magnitude ($10\le\chi\le10^{4.5}$ multiples of the ${\rm Draine}$ field), from which we generated a grid of about 1400 simulations. We compare our results with Herschel SPIRE-FTS, CSO and ISO-LWS observations of eight CO transition lines ($J=1-0$ to $8-7$) and the [CI] 609\mu m and [OI] 146\mu m fine structure lines. We find that the molecular and atomic emission lines trace different gas components of NGC 4038, thus single emission models are insufficient to reproduce the observed values. In general, low-$J$ CO transition lines correspond to either low density regions interacting with low UV field strengths, or high density regions interacting with high UV field strengths. Higher $J$ CO transition lines are less dependent on the UV field strength and are fitted by gas with density $n\sim10^{4.5}-10^{5.2}\,{\rm cm}^{-3}$. We find that the observed fine structure line ratio of [CI] 609\mu m/[OI] 146\mu m is reproduced by clouds subject to weaker UV fields compared to the CO lines. We make estimates of the $X_{\rm CO}$ factor which relates the CO emission with the column density of molecular hydrogen, and find that it is less than the canonical Milky Way value.

We present a comprehensive analysis of a broad band spectral line survey of the Orion Kleinmann-Low nebula (Orion KL), one of the most chemically rich regions in the Galaxy, using the HIFI instrument on board the Herschel Space Observatory. This survey spans a frequency range from 480 to 1907 GHz at a resolution of 1.1 MHz. These observations thus encompass the largest spectral coverage ever obtained toward this high-mass star-forming region in the sub-mm with high spectral resolution, and include frequencies $>$ 1 THz where the Earth's atmosphere prevents observations from the ground. In all, we detect emission from 39 molecules (79 isotopologues). Combining this dataset with ground based mm spectroscopy obtained with the IRAM 30 m telescope, we model the molecular emission from the mm to the far-IR using the XCLASS program which assumes local thermodynamic equilibrium (LTE). Several molecules are also modeled with the MADEX non-LTE code. Because of the wide frequency coverage, our models are constrained by transitions over an unprecedented range in excitation energy. A reduced $\chi^{2}$ analysis indicates that models for most species reproduce the observed emission well. In particular, most complex organics are well fit by LTE implying gas densities are high ($>$10$^6$ cm$^{-3}$) and excitation temperatures and column densities are well constrained. Molecular abundances are computed using H$_{2}$ column densities also derived from the HIFI survey. The distribution of rotation temperatures, $T_{\rm rot}$, for molecules detected toward the hot core is significantly wider than the compact ridge, plateau, and extended ridge $T_{\rm rot}$ distributions, indicating the hot core has the most complex thermal structure.

In order to study the temperature distribution of the extended gas within the Orion Kleinmann-Low nebula, we have mapped the emission by methyl cyanide (CH3CN) in its J=6_K-5_K, J=12_K-11_K, J=13_K-12_K, and J=14_K-13_K transitions at an average angular resolution of ~10 arcsec (22 arcsec for the 6_K-5_K lines), as part of a new 2D line survey of this region using the IRAM 30m telescope. These fully sampled maps show extended emission from warm gas to the northeast of IRc2 and the distinct kinematic signatures of the hot core and compact ridge source components. We have constructed population diagrams for the four sets of K-ladder emission lines at each position in the maps and have derived rotational excitation temperatures and total beam-averaged column densities from the fitted slopes. In addition, we have fitted LVG model spectra to the observations to determine best-fit physical parameters at each map position, yielding the distribution of kinetic temperatures across the region. The resulting temperature maps reveal a region of hot (T > 350 K) material surrounding the northeastern edge of the hot core, whereas the column density distribution is more uniform and peaks near the position of IRc2. We attribute this region of hot gas to shock heating caused by the impact of outflowing material from active star formation in the region, as indicated by the presence of broad CH3CN lines. This scenario is consistent with predictions from C-shock chemical models that suggest that gas-phase methyl cyanide survives in the post-shock gas and can be somewhat enhanced due to sputtering of grain mantles in the passing shock front.

We explore the utility of CI as an alternative high-fidelity gas mass tracer for Galactic molecular clouds. We evaluate the X$_{\rm CI}$-factor for the 609 $\mu$m carbon line, the analog of the CO X-factor, which is the ratio of the H$_2$ column density to the integrated $^{12}$CO(1-0) line intensity. We use 3D-PDR to post-process hydrodynamic simulations of turbulent, star-forming clouds. We compare the emission of CI and CO for model clouds irradiated by 1 and 10 times the average background and demonstrate that CI is a comparable or superior tracer of the molecular gas distribution for column densities up to $6 \times 10^{23}$ cm$^{-2}$. Our results hold for both reduced and full chemical networks. For our fiducial Galactic cloud we derive an average $X_{\rm CO}$ of $3.0\times 10^{20}$ cm$^{-2}$K$^{-1}$km$^{-1}$s and $X_{\rm CI}$ of $1.1\times 10^{21}$ cm$^{-2}$K$^{-1}$km$^{-1}$s.

We present a study of cyanoacetylene (HC3N) and cyanodiacetylene (HC5N) in Orion KL, through observations from two line surveys performed with the IRAM 30m telescope and the HIFI instrument on board the Herschel telescope. The frequency ranges covered are 80-280 GHz and 480-1906 GHz. We model the observed lines of HC3N, HC5N, their isotopologues (including DC3N), and vibrational modes, using a non-LTE radiative transfer code. To investigate the chemical origin of HC3N and DC3N in Orion KL, we use a time-dependent chemical model. We detect 40 lines of the ground state of HC3N and 68 lines of its 13C isotopologues. We also detect 297 lines of six vibrational modes of this molecule (nu_7, 2nu_7, 3nu_7, nu_6, nu_5, and nu_6+nu_7) and 35 rotational lines of the ground state of HC5N. We report the first tentative detection of DC3N in a giant molecular cloud with a DC3N/HC3N abundance ratio of 0.015. We provide column densities and isotopic and molecular abundances. We also perform a 2x2" map around Orion IRc2 and we present maps of HC3N lines and maps of lines of the HC3N vibrational modes nu_6 and nu_7. In addition, a comparison of our results for HC3N with those in other clouds allows us to derive correlations between the column density, the FWHM, the mass, and the luminosity of the clouds. The high column densities of HC3N obtained in the hot core, make this molecule an excellent tracer of hot and dense gas. In addition, the large frequency range covered reveals the need to consider a temperature and density gradient in the hot core in order to obtain better line fits. The high D/H ratio (comparable to that obtained in cold clouds) that we derive suggests a deuterium enrichment. Our chemical models indicate that the possible deuterated HC3N present in Orion KL is formed during the gas-phase. This fact provides new hints concerning the processes leading to deuteration.

We present a study of the sulfur-bearing species detected in a line confusion-limited survey towards Orion KL performed with the IRAM 30m telescope in the range 80-281 GHz. The study is part of an analysis of the line survey divided into families of molecules. Our aim is to derive accurate physical conditions and molecular abundances in the different components of Orion KL from observed SO and SO2 lines. First we assumed LTE conditions obtain rotational temperatures. We then used a radiative transfer model, assuming either LVG or LTE excitation to derive column densities of these molecules in the different components of Orion KL. We have detected 68 lines of SO, 34SO, 33SO, and S18O and 653 lines of SO2, 34SO2, 33SO2, SO18O and SO2 v2=1. We provide column densities for all of them and also upper limits for the column densities of S17O, 36SO, 34S18O, SO17O and 34SO2 v2=1 and for several undetected sulfur-bearing species. In addition, we present 2'x2' maps around Orion IRc2 of SO2 transitions with energies from 19 to 131 K and also maps with four transitions of SO, 34SO and 34SO2. We observe an elongation of the gas along the NE-SW direction. An unexpected emission peak appears at 20.5 km/s in most lines of SO and SO2. A study of the spatial distribution of this emission feature shows that it is a new component ~5" in diameter, which lies ~4" west of IRc2. We suggest the emission from this feature is related to shocks associated to the BN object. The highest column densities for SO and SO2 are found in the high-velocity plateau (a region dominated by shocks) and in the hot core. These values are up to three orders of magnitude higher than the results for the ridge components. We also find high column densities for their isotopologues in both components. Therefore, we conclude that SO and SO2 are good tracers, not only of regions affected by shocks, but also of regions with warm dense gas.

We present a 52-671um spectral scan toward SgrA* taken with the PACS and SPIRE spectrometers onboard Herschel. The achieved angular resolution allows us to separate, for the first time at far-IR wavelengths, the emission toward the central cavity (gas in the inner central parsec of the galaxy) from that of the surrounding circum-nuclear disk. The spectrum toward SgrA* is dominated by strong [OIII], [OI], [CII], [NIII], [NII], and [CI] fine structure lines (in decreasing order of luminosity) arising in gas irradiated by UV-photons from the central stellar cluster. In addition, rotationally excited lines of 12CO (from J=4-3 to 24-23), 13CO, H2O, OH, H3O+, HCO+ and HCN, as well as ground-state absorption lines of OH+, H2O+, H3O+, CH+, H2O, OH, HF, CH and NH are detected. The excitation of the 12CO ladder is consistent with a hot isothermal component at Tk ~ 10^{3.1} K and n(H2)< 10^4 cm^{-3}. It is also consistent with a distribution of temperature components at higher density with most CO at Tk<300 K. The detected molecular features suggest that, at present, neither very enhanced X-ray, nor cosmic-ray fluxes play a dominant role in the heating of the hot molecular gas. The hot CO component (either the bulk of the CO column or just a small fraction depending on the above scenario) results from a combination of UV- and shock-driven heating. If irradiated dense clumps/clouds do not exist, shocks likely dominate the heating of the hot molecular gas. This is consistent with the high-velocity gas detected toward SgrA*.

We use 3D-PDR, a three-dimensional astrochemistry code for modeling photodissociation regions (PDRs), to post-process hydrodynamic simulations of turbulent, star-forming clouds. We focus on the transition from atomic to molecular gas, with specific attention to the formation and distribution of H, C+, C, H2 and CO. First, we demonstrate that the details of the cloud chemistry and our conclusions are insensitive to the simulation spatial resolution, to the resolution at the cloud edge, and to the ray angular resolution. We then investigate the effect of geometry and simulation parameters on chemical abundances and find weak dependence on cloud morphology as dictated by gravity and turbulent Mach number. For a uniform external radiation field, we find similar distributions to those derived using a one-dimensional PDR code. However, we demonstrate that a three-dimensional treatment is necessary for a spatially varying external field, and we caution against using one-dimensional treatments for non-symmetric problems. We compare our results with the work of Glover et al. (2010), who self-consistently followed the time evolution of molecule formation in hydrodynamic simulations using a reduced chemical network. In general, we find good agreement with this in situ approach for C and CO abundances. However, the temperature and H2 abundances are discrepant in the boundary regions (Av < 5), which is due to the different number of rays used by the two approaches.

Photodissociation regions (PDRs) define the transition zone between an ionized and a dark molecular region. They consist of neutral gas which interacts with far-ultraviolet radiation and are characterized by strong infrared line emission. Various numerical codes treating one-dimensional PDRs have been developed in the past, simulating the complexity of chemical reactions occurring and providing a better understanding of the structure of a PDR. In this paper we present the three-dimensional code, 3D-PDR, which can treat PDRs of arbitrary density distribution. The code solves the chemistry and the thermal balance self-consistently within a given three-dimensional cloud. It calculates the total heating and cooling functions at any point in a given PDR by adopting an escape probability method. It uses a HEALPix-based ray-tracing scheme to evaluate the attenuation of the far-ultraviolet radiation in the PDR and the propagation of the far-infrared/submm line emission out of the PDR. We present benchmarking results and apply 3D-PDR to i) a uniform-density spherical cloud interacting with a plane-parallel external radiation field, ii) a uniform-density spherical cloud interacting with a two-component external radiation field, and iii) a cometary globule interacting with a plane-parallel external radiation field. We find that the code is able to reproduce the benchmarking results of various other one-dimensional numerical codes treating PDRs. We also find that the accurate treatment of the radiation field in the fully three-dimensional treatment of PDRs can in some cases leads to different results when compared to a standard one-dimensional treatment.

Ultra-high spectral resolution observations of time-varying interstellar absorption towards {\kappa} Vel are reported, using the Ultra-High Resolution Facility on the Anglo-Australian Telescope. Detections of interstellar Ca I, Ca II, K I, Na I and CH are obtained, whilst an upper limit on the column density is reported for C_2. The results show continued increases in column densities of K I and Ca I since observations ~ 4 yr earlier, as the transverse motion of the star carried it ~ 10 AU perpendicular to the line of sight. Line profile models are fitted to the spectra and two main narrow components (A & B) are identified for all species except CH. The column density N(K I) is found to have increased by 82 +10-9 % between 1994 and 2006, whilst N(Ca I) is found to have increased by 32 +- 5 % over the shorter period of 2002-2006. The line widths are used to constrain the kinetic temperature to T_k,A < 671 +18-17 K and T_k,B < 114 +15-14 K. Electron densities are determined from the Ca I / Ca II ratio, which in turn place lower limits on the total number density of n_A > 7 * 10^3 cm^-3 and n_B > 2 * 10^4 cm^-3. Calcium depletions are estimated from the Ca I / K I ratio. Comparison with the chemical models of Bell et al. (2005) confirms the high number density, with n = 5 * 10^4 cm^-3 for the best-fitting model. The first measurements of diffuse interstellar bands (DIBs) towards this star are made at two epochs, but only an upper limit of < 40 % is placed on their variation over ~ 9 years. The DIBs are unusually weak for the measured E(B-V) and appear to exhibit similar behaviour to that seen in Orion. The ratio of equivalent widths of the {\lambda}5780 to {\lambda}5797 DIBs is amongst the highest known, which may indicate that the carrier of {\lambda}5797 is more sensitive to UV radiation than to local density.

We have used the Herschel-HIFI instrument to observe interstellar nitrogen hydrides along the sight-lines towards W49N and G10.6-0.4 in order to elucidate the production pathways leading to nitrogen-bearing species in diffuse gas. All detections show absorption by foreground material over a wide range of velocities, as well as absorption associated directly with the hot-core source itself. As in the previously published observations towards G10.6-0.4, the NH, NH2 and NH3 spectra towards W49N show strikingly similar and non-saturated absorption features. We decompose the absorption of the foreground material towards W49N into different velocity components in order to investigate whether the relative abundances vary among the velocity components, and, in addition, we re-analyse the absorption lines towards G10.6-0.4 in the same manner. Abundances, with respect to molecular hydrogen, in each velocity component are estimated using CH. The analysis points to a co-existence of the nitrogen hydrides in diffuse or translucent interstellar gas with a high molecular fraction. Towards both sources, we find that NH is always at least as abundant as both o-NH2 and o-NH3, in sharp contrast to previous results for dark clouds. We find relatively constant N(NH)/N(o-NH3) and N(o-NH2)/N(o-NH3) ratios with mean values of 3.2 and 1.9 towards W49N, and 5.4 and 2.2 towards G10.6-0.4, respectively. The mean abundance of o-NH3 is ~2x10^-9 towards both sources. The nitrogen hydrides also show linear correlations with CN and HNC towards both sources, and looser correlations with CH. The upper limits on the NH+ abundance indicate column densities < 2 - 14 % of N(NH). Surprisingly low values of the ammonia ortho-to-para ratio are found in both sources, ~0.5 - 0.7 +- 0.1. This result cannot be explained by current models as we had expected to find a value of unity or higher.

We combine the galaxy formation model GALFORM with the Photon Dominated Region code UCL_PDR to study the emission from the rotational transitions of 12CO (CO) in galaxies from z=0 to z=6 in the Lambda CDM framework. GALFORM is used to predict the molecular (H2) and atomic hydrogen (HI) gas contents of galaxies using the pressure-based empirical star formation relation of Blitz & Rosolowsky. From the predicted H2 mass and the conditions in the interstellar medium, we estimate the CO emission in the rotational transitions 1-0 to 10-9 by applying the UCL_PDR model to each galaxy. We find that deviations from the Milky-Way CO-H2 conversion factor come mainly from variations in metallicity, and in the average gas and star formation rate surface densities. In the local universe, the model predicts a CO(1-0) luminosity function (LF), CO-to-total infrared (IR) luminosity ratios for multiple CO lines and a CO spectral line energy distribution (SLED) which are in good agreement with observations of luminous and ultra-luminous IR galaxies. At high redshifts, the predicted CO SLED of the brightest IR galaxies reproduces the shape and normalization of the observed CO SLED. The model predicts little evolution in the CO-to-IR luminosity ratio for different CO transitions, in good agreement with observations up to z~5. We use this new hybrid model to explore the potential of using colour selected samples of high-redshift star-forming galaxies to characterise the evolution of the cold gas mass in galaxies through observations with the Atacama Large Millimeter Array.

We present for the first time a detailed study of the properties of molecular gas in metal-rich environments such as early-type galaxies (ETGs). We have explored Photon-Dominated Region (PDR) chemistry for a wide range of physical conditions likely to be appropriate for these sources. We derive fractional abundances of the 20 most chemically reactive species as a function of the metallicity, as a function of the optical depth and for various volume number gas densities, Far-Ultra Violet (FUV) radiation fields and cosmic ray ionisation rates. We also investigate the response of the chemistry to the changes in $\alpha-$element enhancement as seen in ETGs. We find that the fractional abundances of CS, H$_{2}$S, H$_{2}$CS, H$_{2}$O, H$_{3}$O$^{+}$, HCO$^{+}$ and H$_{2}$CN seem invariant to an increase of metallicity whereas C$^{+}$, CO, C$_{2}$H, CN, HCN, HNC and OCS appear to be the species most sensitive to this change. The most sensitive species to the change in the fractional abundance of $\alpha-$elements are C$^{+}$, C, CN, HCN, HNC, SO, SO$_{2}$, H$_{2}$O and CS. Finally, we provide line brightness ratios for the most abundant species, especially in the range observable with ALMA. Discussion of favorable line ratios to use for the estimation of super-solar metallicities and $\alpha$-elements are also provided.

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.

Herschel/HIFI observations have revealed the presence of widespread absorption by hydrogen fluoride (HF) J = 1-0 rotational transition, toward a number of Galactic sources. We present observations of HF J = 1-0 toward the high-mass star-forming region Sagittarius B2(M). The spectrum obtained shows a complex pattern of absorption, with numerous features covering a wide range of local standard of rest velocities (-130 to 100 km -1). An analysis of this absorption yields HF abundances relative to H2 of ~1.3 {\times}10-8, in most velocity intervals. This result is in good agreement with estimates from chemical models, which predict that HF should be the main reservoir of gas-phase fluorine under a wide variety of interstellar conditions. Interestingly, we also find velocity intervals in which the HF spectrum shows strong absorption features that are not present, or are very weak, in spectra of other molecules, such as 13CO (1-0) and CS (2-1). HF absorption reveals components of diffuse clouds with small extinction that can be studied for the first time. Another interesting observation is that water is significantly more abundant than hydrogen fluoride over a wide range of velocities toward Sagittarius B2(M), in contrast to the remarkably constant H2O/HF abundance ratio with average value close to unity measured toward other Galactic sources.

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.

As part of a far-infrared (FIR) spectral scan with Herschel/PACS, we present the first detection of the hydroxyl radical (OH) towards the Orion Bar photodissociation region (PDR). Five OH rotational Lambda-doublets involving energy levels out to E_u/k~511 K have been detected (at ~65, ~79, ~84, ~119 and ~163um). The total intensity of the OH lines is I(OH)~5x10^-4 erg s^-1 cm^-2 sr^-1. The observed emission of rotationally excited OH lines is extended and correlates well with the high-J CO and CH^+ J=3-2 line emission (but apparently not with water vapour), pointing towards a common origin. Nonlocal, non-LTE radiative transfer models including excitation by the ambient FIR radiation field suggest that OH arises in a small filling factor component of warm (Tk~160-220 K) and dense (n_H~10^{6-7} cm^-3) gas with source-averaged OH column densities of ~10^15 cm^-2. High density and temperature photochemical models predict such enhanced OH columns at low depths (A_V<1) and small spatial scales (~10^15 cm), where OH formation is driven by gas-phase endothermic reactions of atomic oxygen with molecular hydrogen. We interpret the extended OH emission as coming from unresolved structures exposed to far-ultraviolet (FUV) radiation near the Bar edge (photoevaporating clumps or filaments) and not from the lower density "interclump" medium. Photodissociation leads to OH/H2O abundance ratios (>1) much higher than those expected in equally warm regions without enhanced FUV radiation fields.

The influence of turbulent mixing on the chemistry of the interstellar medium has so far received little attention. Previous studies of this effect have suggested that it might play an important role in mixing the various phases of the interstellar medium. In this paper we examine the potential effects of turbulent diffusion on the deuterium chemistry within molecular clouds. We find that such mixing acts to reduce the efficiency of deuteration in these clouds by increasing the ionization fraction and reducing freeze-out of heavy molecules. This leads to lower abundances for many deuterated species. We also examine the influence of turbulent mixing on the transition from atomic hydrogen to H2 and from atomic deuterium to HD near the cloud edge. We find that including turbulent diffusion in our models serves to push these transitions deeper into the cloud and helps maintain a higher atomic fraction throughout the cloud envelope. Based on these findings, we propose a new process to account for the significant scatter in the observed atomic D/H ratio for galactic sightlines extending beyond the Local Bubble. Although several mechanisms have been put forward to explain this scatter, they are unable to fully account for the range in D/H values. We suggest a scenario in which turbulent mixing of atomic and molecular gas at the edges of molecular clouds causes the observed atomic D/H ratio to vary by a factor of ~2.

We have examined methanol emission from Orion KL with of the {\em Herschel}/HIFI instrument, and detected two methanol bands centered at 524 GHz and 1061 GHz. The 524 GHz methanol band (observed in HIFI band 1a) is dominated by the isolated $\Delta$J$=$0, K$=-4\rightarrow$-3, v$_t$$=$0 Q branch, and includes 25 E-type and 2 A-type transitions. The 1061 GHz methanol band (observed in HIFI band 4b) is dominated by the $\Delta$J$=$0, K$=7\rightarrow$6, v$_t$$=$0 Q branch transitions which are mostly blended. We have used the isolated E-type v$_t$$=$0 methanol transitions to explore the physical conditions in the molecular gas. With HIFI's high velocity resolution, the methanol emission contributed by different spatial components along the line of sight toward Orion KL (hot core, low velocity flow, and compact ridge) can be distinguished and studied separately. The isolated transitions detected in these bands cover a broad energy range (upper state energy ranging from 80 K to 900 K), which provides a unique probe of the thermal structure in each spatial component. The observations further show that the compact ridge is externally heated. These observations demonstrate the power of methanol lines as probes of the physical conditions in warm regions in close proximity to young stars.

We report observations of the reactive molecular ions OH$^+$, H$_2$O$^+$, and H$_3$O$^+$ towards Orion KL with Herschel/HIFI. All three $N=1-0$ fine-structure transitions of OH$^+$ at 909, 971, and 1033GHz and both fine-structure components of the doublet {\it ortho}-H$_2$O$^+$ $1_{11}-0_{00}$ transition at 1115 and 1139GHz were detected; an upper limit was obtained for H$_3$O$^+$. OH$^+$ and H$_2$O$^+$ are observed purely in absorption, showing a narrow component at the source velocity of 9 kms$^{-1}$, and a broad blueshifted absorption similar to that reported recently for HF and {\it para}-H$_{2}^{18}$O, and attributed to the low velocity outflow of Orion KL. We estimate column densities of OH$^+$ and H$_2$O$^+$ for the 9 km s$^{-1}$ component of $9 \pm 3 \times 10^{12}$cm$^{-2}$ and $7 \pm 2 \times 10^{12}$cm$^{-2}$, and those in the outflow of $1.9 \pm 0.7 \times 10^{13}$cm$^{-2}$ and $1.0 \pm 0.3 \times 10^{13}$cm$^{-2}$. Upper limits of $2.4\times 10^{12}$cm$^{-2}$ and $8.7\times 10^{12}$cm$^{-2}$ were derived for the column densities of {\it ortho} and {\it para}-H$_3$O$^+$ from transitions near 985 and 1657GHz. The column densities of the three ions are up to an order of magnitude lower than those obtained from recent observations of W31C and W49N. The comparatively low column densities may be explained by a higher gas density despite the assumption of a very high ionization rate.

In contrast with numerous studies on the physical and chemical structure of low- and high-mass protostars, much less is known about their intermediate-mass counterparts, a class of objects that could help to elucidate the mechanisms of star formation on both ends of the mass range. We present the first results from a rich HIFI spectral dataset on an intermediate-mass protostar, OMC2-FIR4, obtained in the CHESS (Chemical HErschel SurveyS of star forming regions) key programme. The more than 100 methanol lines detected between 554 and 961 GHz cover a range in upper level energy of 40 to 540 K. Our physical interpretation focusses on the hot core, but likely the cold envelope and shocked regions also play a role in reality, because an analysis of the line profiles suggests the presence of multiple emission components. An upper limit of 10^-6 is placed on the methanol abundance in the hot core, using a population diagram, large-scale source model and other considerations. This value is consistent with abundances previously seen in low-mass hot cores. Furthermore, the highest energy lines at the highest frequencies display asymmetric profiles, which may arise from infall around the hot core.

We report the detection of more than 48 velocity-resolved ground rotational state transitions of H2(16)O, H2(18)O, and H2(17)O - most for the first time - in both emission and absorption toward Orion KL using Herschel/HIFI. We show that a simple fit, constrained to match the known emission and absorption components along the line of sight, is in excellent agreement with the spectral profiles of all the water lines. Using the measured H2(18)O line fluxes, which are less affected by line opacity than their H2(16)O counterparts, and an escape probability method, the column densities of H2(18)O associated with each emission component are derived. We infer total water abundances of 7.4E-5, 1.0E-5, and 1.6E-5 for the plateau, hot core, and extended warm gas, respectively. In the case of the plateau, this value is consistent with previous measures of the Orion-KL water abundance as well as those of other molecular outflows. In the case of the hot core and extended warm gas, these values are somewhat higher than water abundances derived for other quiescent clouds, suggesting that these regions are likely experiencing enhanced water-ice sublimation from (and reduced freeze-out onto) grain surfaces due to the warmer dust in these sources.

Observations of HDO are an important complement for studies of water, because they give strong constraints on the formation processes -- grain surfaces versus energetic process in the gas phase, e.g. in shocks. The HIFI observations of multiple transitions of HDO in Sgr~B2(M) presented here allow the determination of the HDO abundance throughout the envelope, which has not been possible before with ground-based observations only. The abundance structure has been modeled with the spherical Monte Carlo radiative transfer code RATRAN, which also takes radiative pumping by continuum emission from dust into account. The modeling reveals that the abundance of HDO rises steeply with temperature from a low abundance ($2.5\times 10^{-11}$) in the outer envelope at temperatures below 100~K through a medium abundance ($1.5\times 10^{-9}$) in the inner envelope/outer core, at temperatures between 100 and 200~K, and finally a high abundance ($3.5\times 10^{-9}$) at temperatures above 200~K in the hot core.

To investigate the accretion and feedback processes in massive star formation, we analyze the shapes of emission lines from hot molecular cores, whose asymmetries trace infall and expansion motions. The high-mass star forming region SgrB2(M) was observed with Herschel/HIFI (HEXOS key project) in various lines of HCN and its isotopologues, complemented by APEX data. The observations are compared to spherically symmetric, centrally heated models with density power-law gradient and different velocity fields (infall or infall+expansion), using the radiative transfer code RATRAN. The HCN line profiles are asymmetric, with the emission peak shifting from blue to red with increasing J and decreasing line opacity (HCN to H$^{13}$CN). This is most evident in the HCN 12--11 line at 1062 GHz. These line shapes are reproduced by a model whose velocity field changes from infall in the outer part to expansion in the inner part. The qualitative reproduction of the HCN lines suggests that infall dominates in the colder, outer regions, but expansion dominates in the warmer, inner regions. We are thus witnessing the onset of feedback in massive star formation, starting to reverse the infall and finally disrupting the whole molecular cloud. To obtain our result, the THz lines uniquely covered by HIFI were critically important.

In contrast to extensively studied dense star-forming cores, little is known about diffuse gas surrounding star-forming regions. We study molecular gas in the high-mass star-forming region NGC6334I, which contains diffuse, quiescent components that are inconspicuous in widely used molecular tracers such as CO. We present Herschel/HIFI observations of CH toward NGC6334I observed as part of the CHESS key program. HIFI resolves the hyperfine components of its J=3/2-1/2 transition, observed in both emission and absorption. The CH emission appears close to the systemic velocity of NGC6334I, while its measured linewidth of 3 km/s is smaller than previously observed in dense gas tracers such as NH3 and SiO. The CH abundance in the hot core is 7 10^-11, two to three orders of magnitude lower than in diffuse clouds. While other studies find distinct outflows in, e.g., CO and H2O toward NGC6334I, we do not detect outflow signatures in CH. To explain the absorption signatures, at least two absorbing components are needed at -3.0 and +6.5 km/s with N(CH)=7 10^13 and 3 10^13 cm^-2. Two additional absorbing clouds are found at +8.0 and 0.0 km/s, both with N(CH)=2 10^13 cm^-2. Turbulent linewidths for the four absorption components vary between 1.5 and 5.0 km/s in FWHM. We constrain physical properties of our CH clouds by matching our CH absorbers with other absorption signatures. In the hot core, molecules such as H2O and CO trace gas that is heated and dynamically influenced by outflow activity, whereas CH traces more quiescent material. The four CH absorbers have column densities and turbulent properties consistent with diffuse clouds: two are located near NGC6334, and two are unrelated foreground clouds. Local density and dynamical effects influence the chemical composition of physical components of NGC6334, causing some components to be seen in CH but not in other tracers, and vice versa.

In the past decade, much progress has been made in characterising the processes leading to the enhanced deuterium fractionation observed in the ISM and in particular in the cold, dense parts of star forming regions such as protostellar envelopes. Very high molecular D/H ratios have been found for saturated molecules and ions. However, little is known about the deuterium fractionation in radicals, even though simple radicals often represent an intermediate stage in the formation of more complex, saturated molecules. The imidogen radical NH is such an intermediate species for the ammonia synthesis in the gas phase. Herschel/HIFI represents a unique opportunity to study the deuteration and formation mechanisms of such species, which are not observable from the ground. We searched here for the deuterated radical ND in order to determine the deuterium fractionation of imidogen and constrain the deuteration mechanism of this species. We observed the solar-mass Class 0 protostar IRAS16293-2422 with the heterodyne instrument HIFI as part of the Herschel key programme CHESS (Chemical HErschel Surveys of Star forming regions). The deuterated form of the imidogen radical ND was detected and securely identified with 2 hyperfine component groups of its fundamental transition in absorption against the continuum background emitted from the nascent protostar. The 3 groups of hyperfine components of its hydrogenated counterpart NH were also detected in absorption. We derive a very high deuterium fractionation with an [ND]/[NH] ratio of between 30 and 100%. The deuterium fractionation of imidogen is of the same order of magnitude as that in other molecules, which suggests that an efficient deuterium fractionation mechanism is at play. We discuss two possible formation pathways for ND, by means of either the reaction of N+ with HD, or deuteron/proton exchange with NH.

We present the first high spectral resolution observations of Orion KL in the frequency ranges 1573.4 - 1702.8 GHz (band 6b) and 1788.4 - 1906.8 GHz (band 7b) obtained using the HIFI instrument on board the Herschel Space Observatory. We characterize the main emission lines found in the spectrum, which primarily arise from a range of components associated with Orion KL including the hot core, but also see widespread emission from components associated with molecular outflows traced by H2O, SO2, and OH. We find that the density of observed emission lines is significantly diminished in these bands compared to lower frequency Herschel/HIFI bands.

Despite the low elemental deuterium abundance in the Galaxy, enhanced molecular D/H ratios have been found in the environments of low-mass star forming regions, and in particular the Class 0 protostar IRAS 16293-2422. The CHESS (Chemical HErschel Surveys of Star forming regions) Key Program aims at studying the molecular complexity of the interstellar medium. The high sensitivity and spectral resolution of the HIFI instrument provide a unique opportunity to observe the fundamental 1,1,1 - 0,0,0 transition of the ortho-D2O molecule, inaccessible from the ground, and to determine the ortho-to-para D2O ratio. We have detected the fundamental transition of the ortho-D2O molecule at 607.35 GHz towards IRAS 16293-2422. The line is seen in absorption with a line opacity of 0.62 +/- 0.11 (1 sigma). From the previous ground-based observations of the fundamental 1,1,0 - 1,0,1 transition of para-D2O seen in absorption at 316.80 GHz we estimate a line opacity of 0.26 +/- 0.05 (1 sigma). We show that the observed absorption is caused by the cold gas in the envelope of the protostar. Using these new observations, we estimate for the first time the ortho to para D2O ratio to be lower than 2.6 at a 3 sigma level of uncertainty, to be compared with the thermal equilibrium value of 2:1.

We present observations of twelve rotational transitions of H2O-16, H2O-18, and H2O-17 toward the massive star-forming region NGC 6334 I, carried out with Herschel/HIFI as part of the guaranteed time key program Chemical HErschel Surveys of Star forming regions (CHESS). We analyze these observations to obtain insights into physical processes in this region. We identify three main gas components (hot core, cold foreground, and outflow) in NGC 6334 I and derive the physical conditions in these components. The hot core, identified by the emission in highly excited lines, shows a high excitation temperature of 200 K, whereas water in the foreground component is predominantly in the ortho- and para- ground states. The abundance of water varies between 4 10^-5 (outflow) and 10^-8 (cold foreground gas). This variation is most likely due to the freeze-out of water molecules onto dust grains. The H2O-18/H2O-17 abundance ratio is 3.2, which is consistent with the O-18/O-17 ratio determined from CO isotopologues. The ortho/para ratio in water appears to be relatively low 1.6(1) in the cold, quiescent gas, but close to the equilibrium value of three in the warmer outflow material (2.5(0.8)).

The HIFI instrument on board the Herschel Space Observatory has been used to observe interstellar nitrogen hydrides along the sight-line towards G10.6-0.4 in order to improve our understanding of the interstellar chemistry of nitrogen. We report observations of absorption in NH N=1-0, J=2-1 and ortho-NH2 1_1,1-0_0,0. We also observed ortho-NH3 1_0-0_0, and 2_0-1_0, para-NH3 2_1-1_1, and searched unsuccessfully for NH+. All detections show emission and absorption associated directly with the hot-core source itself as well as absorption by foreground material over a wide range of velocities. All spectra show similar, non-saturated, absorption features, which we attribute to diffuse molecular gas. Total column densities over the velocity range 11-54 km/s are estimated. The similar profiles suggest fairly uniform abundances relative to hydrogen, approximately 6*10^-9, 3*10^-9, and 3*10^-9 for NH, NH2, and NH3, respectively. These abundances are discussed with reference to models of gas-phase and surface chemistry.

We present initial results from the Herschel GT key program: Herschel observations of EXtra-Ordinary Sources (HEXOS) and outline the promise and potential of spectral surveys with Herschel/HIFI. The HIFI instrument offers unprecedented sensitivity, as well as continuous spectral coverage across the gaps imposed by the atmosphere, opening up a largely unexplored wavelength regime to high-resolution spectroscopy. We show the spectrum of Orion KL between 480 and 560 GHz and from 1.06 to 1.115 THz. From these data, we confirm that HIFI separately measures the dust continuum and spectrally resolves emission lines in Orion KL. Based on this capability we demonstrate that the line contribution to the broad-band continuum in this molecule-rich source is ~20-40% below 1 THz and declines to a few percent at higher frequencies. We also tentatively identify multiple transitions of HD18O in the spectra. The first detection of this rare isotopologue in the interstellar medium suggests that HDO emission is optically thick in the Orion hot core with HDO/H2O ~ 0.02. We discuss the implications of this detection for the water D/H ratio in hot cores.

We performed a sensitive search for the ground-state emission lines of ortho- and para-water vapor in the DM Tau protoplanetary disk using the Herschel/HIFI instrument. No strong lines are detected down to 3sigma levels in 0.5 km/s channels of 4.2 mK for the 1_{10}--1_{01} line and 12.6 mK for the 1_{11}--0_{00} line. We report a very tentative detection, however, of the 1_{10}--1_{01} line in the Wide Band Spectrometer, with a strength of T_{mb}=2.7 mK, a width of 5.6 km/s and an integrated intensity of 16.0 mK km/s. The latter constitutes a 6sigma detection. Regardless of the reality of this tentative detection, model calculations indicate that our sensitive limits on the line strengths preclude efficient desorption of water in the UV illuminated regions of the disk. We hypothesize that more than 95-99% of the water ice is locked up in coagulated grains that have settled to the midplane.

We report a detection of the fundamental rotational transition of hydrogen fluoride in absorption towards Orion KL using Herschel/HIFI. After the removal of contaminating features associated with common molecules ("weeds"), the HF spectrum shows a P-Cygni profile, with weak redshifted emission and strong blue-shifted absorption, associated with the low-velocity molecular outflow. We derive an estimate of 2.9 x 10^13 cm^-2 for the HF column density responsible for the broad absorption component. Using our best estimate of the H2 column density within the low-velocity molecular outflow, we obtain a lower limit of ~1.6 x 10^-10 for the HF abundance relative to hydrogen nuclei, corresponding to 0.6% of the solar abundance of fluorine. This value is close to that inferred from previous ISO observations of HF J=2--1 absorption towards Sgr B2, but is in sharp contrast to the lower limit of 6 x 10^-9 derived by Neufeld et al. (2010) for cold, foreground clouds on the line of sight towards G10.6-0.4.

We discuss the detection of absorption by interstellar hydrogen fluoride (HF) along the sight line to the submillimeter continuum sources W49N and W51. We have used Herschel's HIFI instrument in dual beam switch mode to observe the 1232.4762 GHz J = 1 - 0 HF transition in the upper sideband of the band 5a receiver. We detected foreground absorption by HF toward both sources over a wide range of velocities. Optically thin absorption components were detected on both sight lines, allowing us to measure - as opposed to obtain a lower limit on - the column density of HF for the first time. As in previous observations of HF toward the source G10.6-0.4, the derived HF column density is typically comparable to that of water vapor, even though the elemental abundance of oxygen is greater than that of fluorine by four orders of magnitude. We used the rather uncertain N(CH)-N(H2) relationship derived previously toward diffuse molecular clouds to infer the molecular hydrogen column density in the clouds exhibiting HF absorption. Within the uncertainties, we find that the abundance of HF with respect to H2 is consistent with the theoretical prediction that HF is the main reservoir of gas-phase fluorine for these clouds. Thus, hydrogen fluoride has the potential to become an excellent tracer of molecular hydrogen, and provides a sensitive probe of clouds of small H2 column density. Indeed, the observations of hydrogen fluoride reported here reveal the presence of a low column density diffuse molecular cloud along the W51 sight line, at an LSR velocity of ~ 24kms-1, that had not been identified in molecular absorption line studies prior to the launch of Herschel.

We have observed CH absorption lines ($J=3/2, N=1 \leftarrow J=1/2, N=1$) against the continuum source Sgr~B2(M) using the \textit{Herschel}/HIFI instrument. With the high spectral resolution and wide velocity coverage provided by HIFI, 31 CH absorption features with different radial velocities and line widths are detected and identified. The narrower line width and lower column density clouds show `spiral arm' cloud characteristics, while the absorption component with the broadest line width and highest column density corresponds to the gas from the Sgr~B2 envelope. The observations show that each `spiral arm' harbors multiple velocity components, indicating that the clouds are not uniform and that they have internal structure. This line-of-sight through almost the entire Galaxy offers unique possibilities to study the basic chemistry of simple molecules in diffuse clouds, as a variety of different cloud classes are sampled simultaneously. We find that the linear relationship between CH and H$_2$ column densities found at lower $A_V$ by UV observations does not continue into the range of higher visual extinction. There, the curve flattens, which probably means that CH is depleted in the denser cores of these clouds.

Optical emission is detected from filaments around the central galaxies of clusters of galaxies. These filaments have lengths of tens of kiloparsecs. The emission is possibly due to heating caused by the dissipation of mechanical energy and by cosmic ray induced ionisation. CO millimeter and submillimeter line emissions as well as H$_{2}$ infrared emission originating in such filaments surrounding NGC~1275, the central galaxy of the Perseus cluster, have been detected. Our aim is to identify those molecular species, other than CO, that may emit detectable millimeter and submillimeter line features arising in these filaments, and to determine which of those species will produce emissions that might serve as diagnostics of the dissipation and cosmic ray induced ionisation. The time-dependent UCL photon-dominated region modelling code was used in the construction of steady-state models of molecular filamentary emission regions at appropriate pressures, for a range of dissipation and cosmic ray induced ionisation rates and incident radiation fields.HCO$^+$ and C$_2$H emissions will potentially provide information about the cosmic ray induced ionisation rates in the filaments. HCN and, in particular, CN are species with millimeter and submillimeter lines that remain abundant in the warmest regions containing molecules. Detections of the galaxy cluster filaments in HCO$^{+}$, C$_{2}$H, and CN emissions and further detections of them in HCN emissions would provide significant constraints on the dissipation and cosmic ray induced ionisation rates.

We present Herschel/HIFI observations of the fundamental rotational transitions of ortho- and para-H$_2^{16}$O and H$_2^{18}$O in absorption towards Sagittarius~B2(M) and W31C. The ortho/para ratio in water in the foreground clouds on the line of sight towards these bright continuum sources is generally consistent with the statistical high-temperature ratio of 3, within the observational uncertainties. However, somewhat unexpectedly, we derive a low ortho/para ratio of $2.35 \pm 0.35$, corresponding to a spin temperature of $\sim$27~K, towards Sagittarius~B2(M) at velocities of the expanding molecular ring. Water molecules in this region appear to have formed with, or relaxed to, an ortho/para ratio close to the value corresponding to the local temperature of the gas and dust.

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.

We report the detection of absorption by interstellar hydroxyl cations and water cations, along the sight-line to the bright continuum source W49N. We have used Herschel's HIFI instrument, in dual beam switch mode, to observe the 972 GHz N = 1 - 0 transition of OH+ and the 1115 GHz 1(11) - 0(00) transition of ortho-H2O+. The resultant spectra show absorption by ortho-H2O+, and strong absorption by OH+, in foreground material at velocities in the range 0 to 70 km/s with respect to the local standard of rest. The inferred OH+/H2O+ abundance ratio ranges from ~ 3 to ~ 15, implying that the observed OH+ arises in clouds of small molecular fraction, in the 2 - 8% range. This conclusion is confirmed by the distribution of OH+ and H2O+ in Doppler velocity space, which is similar to that of atomic hydrogen, as observed by means of 21 cm absorption measurements, and dissimilar from that typical of other molecular tracers. The observed OH+/H abundance ratio of a few E-8 suggests a cosmic ray ionization rate for atomic hydrogen of (0.6 - 2.4) E-16 s-1, in good agreement with estimates inferred previously for diffuse clouds in the Galactic disk from observations of interstellar H3+ and other species.

H2O+ has been observed in its ortho- and para- states toward the massive star forming core Sgr B2(M), located close to the Galactic center. The observations show absorption in all spiral arm clouds between the Sun and Sgr B2. The average o/p ratio of H2O+ in most velocity intervals is 4.8, which corresponds to a nuclear spin temperature of 21 K. The relationship of this spin temperature to the formation temperature and current physical temperature of the gas hosting H2O+ is discussed, but no firm conclusion is reached. In the velocity interval 0 to 60 km/s, an ortho/para ratio of below unity is found, but if this is due to an artifact of contamination by other species or real is not clear.

We present spectrally resolved observations of triatomic carbon (C3) in several ro-vibrational transitions between the vibrational ground state and the low-energy nu2 bending mode at frequencies between 1654-1897 GHz along the sight-lines to the submillimeter continuum sources W31C and W49N, using Herschel's HIFI instrument. We detect C3 in absorption arising from the warm envelope surrounding the hot core, as indicated by the velocity peak position and shape of the line profile. The sensitivity does not allow to detect C3 absorption due to diffuse foreground clouds. From the column densities of the rotational levels in the vibrational ground state probed by the absorption we derive a rotation temperature (T_rot) of ~50--70 K, which is a good measure of the kinetic temperature of the absorbing gas, as radiative transitions within the vibrational ground state are forbidden. It is also in good agreement with the dust temperatures for W31C and W49N. Applying the partition function correction based on the derived T_rot, we get column densities N(C3) ~7-9x10^{14} cm^{-2} and abundance x(C3)~10^{-8} with respect to H2. For W31C, using a radiative transfer model including far-infrared pumping by the dust continuum and a temperature gradient within the source along the line of sight we find that a model with x(C3)=10^{-8}, T_kin=30-50 K, N(C3)=1.5 10^{15} cm^{-2} fits the observations reasonably well and provides parameters in very good agreement with the simple excitation analysis.

The outflow driven by the low-mass class 0 protostar L1157 is the prototype of the so-called chemically active outflows. The bright bowshock B1 in the southern outflow lobe is a privileged testbed of magneto-hydrodynamical (MHD) shock models, for which dynamical and chemical processes are strongly interdependent. We present the first results of the unbiased spectral survey of the L1157-B1 bowshock, obtained in the framework of the key program "Chemical Herschel Surveys of Star Forming Regions" (CHESS). The main aim is to trace the warm and chemically enriched gas and to infer the excitation conditions in the shock region. The CO 5-4 and H2O lines have been detected at high-spectral resolution in the unbiased spectral survey of the HIFI-Band 1b spectral window (555-636 GHz), presented by Codella et al. in this volume. Complementary ground-based observations in the submm window help establish the origin of the emission detected in the main-beam of HIFI, and the physical conditions in the shock.}{Both lines exhibit broad wings, which extend to velocities much higher than reported up to now. We find that the molecular emission arises from two regions with distinct physical conditions: an extended, warm (100K), dense (3e5 cm-3) component at low-velocity, which dominates the water line flux in Band~1; a secondary component in a small region of B1 (a few arcsec) associated with high-velocity, hot (> 400 K) gas of moderate density ((1.0-3.0)e4 cm-3), which appears to dominate the flux of the water line at 179mu observed with PACS. The water abundance is enhanced by two orders of magnitude between the low- and the high-velocity component, from 8e-7 up to 8e-5. The properties of the high-velocity component agree well with the predictions of steady-state C-shock models.

We present SHARC-2 350 micron data on 20 luminous z ~ 2 starbursts with S(1.2mm) > 2 mJy from the Spitzer-selected samples of Lonsdale et al. and Fiolet et al. All the sources were detected, with S(350um) > 25 mJy for 18 of them. With the data, we determine precise dust temperatures and luminosities for these galaxies using both single-temperature fits and models with power-law mass--temperature distributions. We derive appropriate formulae to use when optical depths are non-negligible. Our models provide an excellent fit to the 6um--2mm measurements of local starbursts. We find characteristic single-component temperatures T1 ~ 35.5+-2.2 K and integrated infrared (IR) luminosities around 10^(12.9+-0.1) Lsun for the SWIRE-selected sources. Molecular gas masses are estimated at 4 x 10^(10) Msun, assuming kappa(850um)=0.15 m^2/kg and a submillimeter-selected galaxy (SMG)-like gas-to-dust mass ratio. The best-fit models imply >~2 kpc emission scales. We also note a tight correlation between rest-frame 1.4 GHz radio and IR luminosities confirming star formation as the predominant power source. The far-IR properties of our sample are indistinguishable from the purely submillimeter-selected populations from current surveys. We therefore conclude that our original selection criteria, based on mid-IR colors and 24 um flux densities, provides an effective means for the study of SMGs at z ~ 1.5--2.5.

We report the detection of absorption lines by the reactive ions OH+, H2O+ and H3O+ along the line of sight to the submillimeter continuum source G10.6$-$0.4 (W31C). We used the Herschel HIFI instrument in dual beam switch mode to observe the ground state rotational transitions of OH+ at 971 GHz, H2O+ at 1115 and 607 GHz, and H3O+ at 984 GHz. The resultant spectra show deep absorption over a broad velocity range that originates in the interstellar matter along the line of sight to G10.6$-$0.4 as well as in the molecular gas directly associated with that source. The OH+ spectrum reaches saturation over most velocities corresponding to the foreground gas, while the opacity of the H2O+ lines remains lower than 1 in the same velocity range, and the H3O+ line shows only weak absorption. For LSR velocities between 7 and 50 kms$^{-1}$ we estimate total column densities of $N$(OH+) $> 2.5 \times 10^{14}$ cm$^{-2}$, $N$(H2O+) $\sim 6 \times 10^{13}$ cm$^{-2}$ and $N$(H3O+) $\sim 4.0 \times 10^{13}$ cm$^{-2}$. These detections confirm the role of O$^+$ and OH$^+$ in initiating the oxygen chemistry in diffuse molecular gas and strengthen our understanding of the gas phase production of water. The high ratio of the OH+ by the H2O+ column density implies that these species predominantly trace low-density gas with a small fraction of hydrogen in molecular form.

We report the detection of strong absorption by interstellar hydrogen fluoride along the sight-line to the submillimeter continuum source G10.6-0.4 (W31C). We have used Herschel's HIFI instrument, in dual beam switch mode, to observe the 1232.4763 GHz J=1-0 HF transition in the upper sideband of the Band 5a receiver. The resultant spectrum shows weak HF emission from G10.6-0.4 at LSR velocities in the range -10 to -3 km/s, accompanied by strong absorption by foreground material at LSR velocities in the range 15 to 50 km/s. The spectrum is similar to that of the 1113.3430 GHz 1(11)-0(00) transition of para-water, although at some frequencies the HF (hydrogen fluoride) optical depth clearly exceeds that of para-H2O. The optically-thick HF absorption that we have observed places a conservative lower limit of 1.6E+14 cm-2 on the HF column density along the sight-line to G10.6-0.4. Our lower limit on the HF abundance, 6E-9 relative to hydrogen nuclei, implies that hydrogen fluoride accounts for between ~ 30 and 100% of the fluorine nuclei in the gas phase along this sight-line. This observation corroborates theoretical predictions that - because the unique thermochemistry of fluorine permits the exothermic reaction of F atoms with molecular hydrogen - HF will be the dominant reservoir of interstellar fluorine under a wide range of conditions.

We identify a prominent absorption feature at 1115 GHz, detected in first HIFI spectra towards high-mass star-forming regions, and interpret its astrophysical origin. The characteristic hyperfine pattern of the H2O+ ground-state rotational transition, and the lack of other known low-energy transitions in this frequency range, identifies the feature as H2O+ absorption against the dust continuum background and allows us to derive the velocity profile of the absorbing gas. By comparing this velocity profile with velocity profiles of other tracers in the DR21 star-forming region, we constrain the frequency of the transition and the conditions for its formation. In DR21, the velocity distribution of H2O+ matches that of the [CII] line at 158\mu\m and of OH cm-wave absorption, both stemming from the hot and dense clump surfaces facing the HII-region and dynamically affected by the blister outflow. Diffuse foreground gas dominates the absorption towards Sgr B2. The integrated intensity of the absorption line allows us to derive lower limits to the H2O+ column density of 7.2e12 cm^-2 in NGC 6334, 2.3e13 cm^-2 in DR21, and 1.1e15 cm^-2 in Sgr B2.

Far-infrared [C ii] 158 micron, [O i] 145 micron and [O i] 63 micron fine structure emission line fluxes were measured from archival Infrared Space Observatory Long Wavelength Spectrometer spectra of 46 extragalactic sources, with 28 sources providing detections in all three lines. For 12 of the sources, the contribution to the [C ii] 158 micron line flux from H ii regions could be estimated from their detected [N ii] 122 micron line fluxes. The measured [C ii]/[O i] and [O i] 63/145 micron line flux ratios were compared with those from a grid of PDR models previously computed using the UCL PDR code. Persistent offsets between the observed and modelled line ratios could be partly attributed to the effects of [O i] 63 micron self-absorption. Using the SMMOL code, we calculated model [O i] line profiles and found that the strength of the [O i] 63 micron line was reduced by 20-80%, depending on the PDR parameters. We conclude that high PDR densities and radiation field strengths, coupled with the effects of [O i] 63 micron self-absorption, are likely to provide the best match to the observed line flux ratios.

Photon-dominated regions (PDRs) are powerful molecular line emitters in external galaxies. They are expected in galaxies with high rates of massive star formation due to either starburst (SB) events or starburst coupled with active galactic nuclei (AGN) events. We have explored the PDR chemistry for a range of physical conditions representing a variety of galaxy types. Our main result is a demonstration of the sensitivity of the chemistry to changes in the physical conditions. We adopt crude estimates of relevant physical parameters for several galaxy types and use our models to predict suitable molecular tracers of those conditions. The set of recommended molecular tracers differs from that which we recommended for use in galaxies with embedded massive stars. Thus, molecular observations can in principle be used to distinguish between excitation by starburst and by SB+AGN in distant galaxies. Our recommendations are intended to be useful in preparing Herschel and ALMA proposals to identify sources of excitation in galaxies.

Residual star formation at late times in early-type galaxies and their progenitors must be suppressed in order to explain the population of red, passively evolving systems we see today. Likewise, residual or newly accreted reservoirs of molecular gas that are fuelling star formation must be destroyed. This suppression of star formation in early-type galaxies is now commonly attributed to AGN feedback wherein the reservoir of gas is heated and expelled during a phase of accretion onto the central supermassive black hole. However, direct observational evidence for a link between the destruction of this molecular gas and an AGN phase has been missing so far. We present new mm-wavelength observations from the IRAM 30m telescope of a sample of low redshift SDSS early-type galaxies currently undergoing this process of quenching of late-time star formation. Our observations show that the disappearance of the molecular gas coincides within less than 100 Myr with the onset of accretion onto the black hole and is too rapid to be due to star formation alone. Since our sample galaxies are not associated to powerful quasar activity or radio jets, we conclude that low-luminosity AGN episodes are sufficient to suppress residual star formation in early-type galaxies. This `suppression mode' of AGN feedback is very different from the `truncation mode' linked to powerful quasar activity during early phases of galaxy formation.

Massive stars undergo a violent death when the supply of nuclear fuel in their cores is exhausted, resulting in a catastrophic "core-collapse" supernova. Such events are usually only detected at least a few days after the star has exploded. Observations of the supernova SNLS-04D2dc with the Galaxy Evolution Explorer space telescope reveal a radiative precursor from the supernova shock before the shock reached the surface of the star and show the initial expansion of the star at the beginning of the explosion. Theoretical models of the ultraviolet light curve confirm that the progenitor was a red supergiant, as expected for this type of supernova. These observations provide a way to probe the physics of core-collapse supernovae and the internal structures of their progenitor stars

We have carried out a systematic search for the molecular ion CO+ in a sample of 8 protoplanetary and planetary nebulae in order to determine the origin of the unexpectedly strong HCO+ emission previously detected in these sources. An understanding of the HCO+ chemistry may provide direct clues to the physical and chemical evolution of planetary nebulae. We find that the integrated intensity of the CO+ line may be correlated with that of HCO+, suggesting that the reaction of CO+ with molecular hydrogen may be an important formation route for HCO+ in these planetary nebulae.

We present theoretically-established values of the CO-to-H2 and C-to-H2 conversion factors that may be used to estimate the gas masses of external galaxies. We consider four distinct galaxy types, represented by M51, NGC 6946, M82 and SMC N27. The physical parameters that best represent the conditions within the molecular clouds in each of the galaxy types are estimated using a chi^2 analysis of several observed atomic fine structure and CO rotational lines. This analysis is explored over a wide range of density, radiation field, extinction, and other relevant parameters. Using these estimated physical conditions in methods that we have previously established, CO-to-H2 conversion factors are then computed for CO transitions up to J=9-8. For the conventional CO(1-0) transition, the computed conversion factor varies significantly below and above the canonical value for the Milky Way in the four galaxy types considered. Since atomic carbon emission is now frequently used as a probe of external galaxies, we also present, for the first time, the C-to-H2 conversion factor for this emission in the four galaxy types considered.

We present a comparison between independent computer codes, modeling the physics and chemistry of interstellar photon dominated regions (PDRs). Our goal was to understand the mutual differences in the PDR codes and their effects on the physical and chemical structure of the model clouds, and to converge the output of different codes to a common solution. A number of benchmark models have been created, covering low and high gas densities and far ultraviolet intensities. The benchmark models were computed in two ways: one set assuming constant temperatures, thus testing the consistency of the chemical network and photo-processes, and a second set determining the temperature selfconsistently. We investigated the impact of PDR geometry and agreed on the comparison of results from spherical and plane-parallel PDR models. We identified a number of key processes governing the chemical network which have been treated differently in the various codes, and defined a proper common treatment. We established a comprehensive set of reference models for ongoing and future PDR model bench-marking and were able to increase the agreement in model predictions for all benchmark models significantly.

Clouds containing molecular dark matter in quantities relevant for star formation may exist in minihaloes of the type of cold dark matter included in many cosmological simulations or in the regions of some galaxies extending far beyond their currently known boundaries. We have systematically explored parameter space to identify conditions under which plane-parallel clouds contain sufficient column densities of molecular dark matter that they could be reservoirs for future star formation. Such clouds would be undetected or at least appear by current observational criteria to be uninteresting from the perspective of star formation. We use a time-dependent PDR code to produce theoretical models of the chemistry and emission arising in clouds for our chosen region of parameter space. We then select a subset of model clouds with levels of emission that are low enough to be undetectable or at least overlooked by current surveys. The existence of significant column densities of cold molecular dark matter requires that the background radiation field be several or more orders of magnitude weaker than that in the solar neighbourhood. Lower turbulent velocities and cosmic ray induced ionization rates than typically associated with molecular material within a kpc of the Sun are also required for the molecular matter to be dark. We find that there is a large region within the parameter space that results in clouds that might contain a significant mass of molecular gas whilst remaining effectively undetectable or at least not particularly noticeable in surveys. We note briefly conditions under which molecular dark matter may contain a dynamically interesting mass.

A reliable estimate of the molecular gas content in galaxies plays a crucial role in determining their dynamical and star-forming properties. However, H2, the dominant molecular species, is difficult to observe directly, particularly in the regions where most molecular gas is thought to reside. Its mass is therefore commonly inferred by assuming a direct proportionality with the integrated intensity of the CO(J=1-0) emission line, using a CO-to-H2 conversion factor, X. Although a canonical value for X is used extensively in such estimates, there is increasing evidence, both theoretical and observational, that the conversion factor may vary by over an order of magnitude under conditions different to those of the local neighbourhood. In an effort to understand the influence of changing environmental conditions on the conversion factor, we derive theoretical estimates of X for a wide range of physical parameters using a photon-dominated region (PDR) time-dependent chemical model, benchmarking key results against those of an independent PDR code to ensure reliability. Based on these results, the sensitivity of the X factor to change in each physical parameter is interpreted in terms of the chemistry and physical processes within the cloud. In addition to confirming previous observationally derived trends, we find that the time-dependence of the chemistry, often neglected in such models, has a considerable influence on the value of the conversion factor.

The effect of a newly born star cluster inside a giant molecular cloud (GMC) is to produce a hot bubble and a thin, dense shell of interstellar gas and dust swept up by the H II expansion, strong stellar winds, and repeated supernova explosions. Lying at the inner side of the shell is the photodissociation region (PDR), the origin of much of the far-infrared/sub-millimeter/millimeter (FIR/sub-mm/mm) radiation from the interstellar medium (ISM). We present a model for the expanding shell at different stages of its expansion which predict mm/sub-mm and far-IR emission line intensities from a series of key molecular and atomic constituents in the shell. The kinematic properties of the swept-up shell predicted by our model are in very good agreement with the measurements of the supershell detected in the nearby starburst galaxy M 82. We compare the modeling results with the ratio-ratio plots of the FIR/sub-mm/mm line emission in the central 1.0 kpc region to investigate the mechanism of star forming activity in M 82. Our model has yielded appropriate gas densities, temperatures, and structure scales compared to those measured in M 82, and the total H2 content is compatible with the observations. This implies that the neutral ISM of the central star-forming region is a product of fragments of the evolving shells.

Transient microstructure in the diffuse interstellar medium has been observed towards galactic and extragalactic sources for decades, usually in lines of atoms and ions, and, more recently, in molecular lines. Evidently, there is a molecular component to the transient microstructure. In this paper we explore the chemistry that may arise in such microstructure. We use a PDR code to model the conditions of relatively high density, low temperature, very low visual extinction and very short elapsed time that are appropriate for these objects. We find that there is a well defined region of parameter space where detectable abundances of molecular species might be found. The best matching models are those where the interstellar microstructure is young (<100 yr), small (~100 AU), and dense (>10^4 cm^-3).