We report on the neutrino mass measurement result from the first four-week science run of the Karlsruhe Tritium Neutrino experiment KATRIN in spring 2019. Beta-decay electrons from a high-purity gaseous molecular tritium source are energy analyzed by a high-resolution MAC-E filter. A fit of the integrated electron spectrum over a narrow interval around the kinematic endpoint at 18.57 keV gives an effective neutrino mass square value of $(-1.0^{+0.9}_{-1.1})$ eV$^2$. From this we derive an upper limit of 1.1 eV (90$\%$ confidence level) on the absolute mass scale of neutrinos. This value coincides with the KATRIN sensitivity. It improves upon previous mass limits from kinematic measurements by almost a factor of two and provides model-independent input to cosmological studies of structure formation.

The determination of the neutrino mass is one of the major challenges in astroparticle physics today. Direct neutrino mass experiments, based solely on the kinematics of beta-decay, provide a largely model-independent probe to the neutrino mass scale. The Karlsruhe Tritium Neutrino (KATRIN) experiment is designed to directly measure the effective electron antineutrino mass with a sensitivity of 0.2 eV 90% CL. In this work we report on the first operation of KATRIN with tritium which took place in 2018. During this commissioning phase of the tritium circulation system, excellent agreement of the theoretical prediction with the recorded spectra was found and stable conditions over a time period of 13 days could be established. These results are an essential prerequisite for the subsequent neutrino mass measurements with KATRIN in 2019.

Understanding the matter cycle in the interstellar medium of galaxies from the assembly of clouds to star formation and stellar feedback remains an important and exciting field in comtemporary astrophysics. Many open questions regarding cloud and structure formation, the role of turbulence, and the relative importance of the various feedback processes can only be addressed with observations of spectrally resolved lines. We here stress the importance of two specific sets of lines: the finestructure lines of atomic carbon as a tracer of the dark molecular gas and mid-J CO lines as tracers of the warm, active molecular gas in regions of turbulence dissipation and feedback. The observations must cover a wide range of environments (i.e., physical conditions), which will be achieved by large scale surveys of Galactic molecular clouds, the Galactic Center, the Magellanic clouds, and nearby galaxies. To date, such surveys are completely missing and thus constitute an important science opportunity for the next decade and beyond. For the successful interpretation of the observations, it will be essential to combine them with results from (chemical) modelling and simulations of the interstellar medium.

The aim of our study is to investigate the physical properties of the star-forming interstellar medium (ISM) in the Large Magellanic Cloud (LMC) by separating the origin of the emission lines spatially and spectrally. Following Okada et al. (2015, Paper I), we investigate different phases of the ISM traced by carbon-bearing species in four star-forming regions in the LMC, and model the physical properties using the KOSMA-tau PDR model. We mapped 3--13 arcmin$^2$ areas in 30 Dor, N158, N160 and N159 along the molecular ridge of the LMC in [CII]158um with GREAT on board SOFIA, and in CO(2-1) to (6-5), $^{13}$CO(2-1) and (3-2), [CI]3P1-3P0 and 3P2-3P1 with APEX. In all four star-forming regions, the line profiles of CO, $^{13}$CO, and [CI] emission are similar, whereas [CII] typically shows wider line profiles or an additional velocity component. For selected positions in N159 and 30 Dor, we observed the velocity-resolved [OI] 145um and 63um lines for the first time with upGREAT. At some positions, the [OI] line profiles match those of CO, at other positions they are more similar to the [CII] profiles. We interpret the different line profiles of CO, [CII] and [OI] as contributions from spatially separated clouds and/or clouds in different physical phases, which give different line ratios depending on their physical properties. We model the emission from the CO, [CI], [CII], and [OI] lines and the far-infrared continuum emission using the latest KOSMA-tau PDR model, which treats the dust-related physics consistently and computes the dust continuum SED together with the line emission of the chemical species. We find that the line and continuum emissions are not well-reproduced by a single clump ensemble. Toward the CO peak at N159~W, we propose a scenario that the CO, [CII], and [OI] 63um emission are weaker than expected because of mutual shielding among clumps. (abridged for arXiv)

We investigate the physical conditions of the gas, atomic and molecular, in the filaments in the context of Photo-Dissociation Regions (PDRs) using the KOSMA-PDR mode of clumpy clouds. We also compare the [CII] vs. [NII] integrated intensity predictions in Abel et al. 2005 for HII regions and adjacent PDRs in the Galactic disk, and check for their applicability under the extreme physical conditions present in the GC. Our preliminary results show that observed integrated intensities are well reproduced by the PDR model. The gas is exposed to a relatively low Far-UV field between 10^2 - 10^3 Draine fields. The total volume hydrogen density is well constrained between 10^4 - 10^5 cm^-3. The hydrogen ionization rate due to cosmic-rays varies between 10^-15 and 4x10^-15 s^-1, with the highest value ~ 10^-14 s^-1 found towards G0.07+0.04. Our results show that the line-of-sight contribution to the total distance of the filaments to the Arches Cluster is not negligible. The spatial distribution of the [CII]/[NII] ratio shows that the integrated intensity ratios are fairly homogeneously distributed for values below 10 in energy units. Calculations including variation on the [C/N] abundance ratio show that tight constraints on this ratio are needed to reproduce the observations.

Context. Models of photon-dominated regions (PDRs) still fail to fully reproduce some of the observed properties, in particular the combination of the intensities of different PDR cooling lines together with the chemical stratification, as observed e.g. for the Orion Bar PDR. Aims. We aim to construct a numerical PDR model, KOSMA-\tau 3D, to simulate full spectral cubes of line emission from arbitrary PDRs in three dimensions (3D). The model is to reproduce the intensity of the main cooling lines from the Orion Bar PDR and the observed layered structure of the different transitions. Methods. We build up a 3D compound, made of voxels ("3D pixels") that contain a discrete mass distribution of spherical "clumpy" structures, approximating the fractal ISM. To analyse each individual clump the new code is combined with the KOSMA-\tau PDR model. Probabilistic algorithms are used to calculate the local FUV flux for each voxel as well as the voxel-averaged line emissivities and optical depths, based on the properties of the individual clumps. Finally, the computation of the radiative transfer through the compound provides full spectral cubes. To test the new model we try to simulate the structure of the Orion Bar PDR and compare the results to observations from HIFI/Herschel and from the Caltech Submillimetre Observatory (CSO). In this context new Herschel data from the HEXOS guaranteed-time key program is presented. Results. Our model is able to reproduce the line integrated intensities within a factor 2.5 and the observed stratification pattern within 0.016 pc for the [Cii] 158 \mu m and different 12/13 CO and HCO+ transitions, based on the representation of the Orion Bar PDR by a clumpy edge-on cavity wall. In the cavity wall, a large fraction of the total mass needs to be contained in clumps. The mass of the interclump medium is constrained by the FUV penetration. Furthermore, ...

[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.

Aims: We investigate how much of the [CII] emission in the nucleus of the nearby spiral galaxy IC 342 is contributed by PDRs and by the ionized gas. We examine the spatial variations of starburst/PDR activity and study the correlation of the [CII] line with the [NII] 205{\textmu}m emission line coming exclusively from the HII regions. Methods: We present small maps of [CII] and [NII] lines recently observed with the GREAT receiver on board SOFIA. In particular we present a super-resolution method to derive how unresolved, kinematically correlated structures in the beam contribute to the observed line shapes. Results: We find that the emission coming from the ionized gas shows a kinematic component in addition to the general Doppler signature of the molecular gas. We interpret this as the signature of two bi-polar lobes of ionized gas expanding out of the galactic plane. We then show how this requires an adaptation of our understanding of the geometrical structure of the nucleus of IC~342. Examining the starburst activity we find ratios $I([CII])/I(^{12}\mathrm{CO} (1-0))$ between 400 and 1800 in energy units. Applying predictions from numerical models of HII and PDR regions to derive the contribution from the ionized phase to the total [CII] emission we find that 35-90% of the observed [CII] intensity stems from the ionized gas if both phases contribute. Averaged over the central few hundred parsec we find for the [CII] contribution a HII-to-PDR ratio of 70:30. Conclusions: The ionized gas in the center of IC 342 contributes more strongly to the overall [CII] emission than is commonly observed on larger scales and than is predicted. Kinematic analysis shows that the majority of the [CII] emission is related to the strong but embedded star formation in the nuclear molecular ring and only marginally emitted from the expanding bi-polar lobes of ionized gas.

We observationally investigate the relation between the photoelectric heating efficiency in PDRs and the charge of PAHs, which are considered to play a key role in photoelectric heating. Using PACS onboard Herschel, we observed six PDRs spanning a wide range of FUV radiation fields (G_0=100-10^5). To measure the photoelectric heating efficiency, we obtained the intensities of the main cooling lines, i.e., the [OI]63um, 145um, and [CII]158um, as well as the FIR continuum intensity. We used Spitzer/IRS spectroscopic mapping observations to investigate the MIR PAH features in the same regions. We decomposed the MIR PAH emission into that of neutral (PAH^0) and positively ionized (PAH^+) species to derive the fraction of the positively charged PAHs, and compare it to the photoelectric heating efficiency. The heating efficiency traced by ([OI]63um+[OI]145um+[CII]158um) / TIR, ranges between 0.1% and 0.9% in different sources, and the fraction of PAH^+ relative to (PAH^0 + PAH^+) spans from 0(+11)% to 87(+/-10)%. All positions with a high PAH^+ fraction show a low heating efficiency, and all positions with a high heating efficiency have a low PAH^+ fraction, supporting the scenario in which a positive grain charge results in a decreased heating efficiency. Theoretical estimates of the photoelectric heating efficiency show a stronger dependence on the charging parameter gamma=G_0 T^{1/2}/n_e than the observed efficiency reported in this study, and the discrepancy is significant at low gamma. The photoelectric heating efficiency on PAHs, traced by ([OI]63um+[OI]145um+[CII]158um) / (PAH+[OI]63um+[OI]145um+[CII]158um), shows a much better match between the observations and the theoretical estimates. The good agreement of the photoelectric heating efficiency on PAHs with a theoretical model indicates the dominant contribution of PAHs to the photoelectric heating. (abridged for arXiv)

We upgraded the chemical network from the UMIST Database for Astrochemistry 2006 to include isotopes such as ^{13}C and ^{18}O. This includes all corresponding isotopologues, their chemical reactions and the properly scaled reaction rate coefficients. We study the fractionation behavior of astrochemically relevant species over a wide range of model parameters, relevant for modelling of photo-dissociation regions (PDRs). We separately analyze the fractionation of the local abundances, fractionation of the total column densities, and fractionation visible in the emission line ratios. We find that strong C^+ fractionation is possible in cool C^+ gas. Optical thickness as well as excitation effects produce intensity ratios between 40 and 400. The fractionation of CO in PDRs is significantly different from the diffuse interstellar medium. PDR model results never show a fractionation ratio of the CO column density larger than the elemental ratio. Isotope-selective photo-dissociation is always dominated by the isotope-selective chemistry in dense PDR gas. The fractionation of C, CH, CH^+, and HCO^+ is studied in detail, showing that the fractionation of C, CH and CH^+ is dominated by the fractionation of their parental species. The light hydrides chemically derive from C^+, and, consequently, their fractionation state is coupled to that of C^+. The fractionation of C is a mixed case depending on whether formation from CO or HCO^+ dominates. Ratios of the emission lines of [C II], [C I], ^{13}CO, and H^{13}CO^+ provide individual diagnostics to the fractionation status of C^+, C, and CO.

We revised the treatment of interstellar dust in the KOSMA-\tau\ PDR model code to achieve a consistent description of the dust-related physics in the code. The detailed knowledge of the dust properties is then used to compute the dust continuum emission together with the line emission of chemical species. We coupled the KOSMA-\tau\ PDR code with the MCDRT (multi component dust radiative transfer) code to solve the frequency-dependent radiative transfer equations and the thermal balance equation in a dusty clump under the assumption of spherical symmetry, assuming thermal equilibrium in calculating the dust temperatures, neglecting non-equilibrium effects. We updated the calculation of the photoelectric heating and extended the parametrization range for the photoelectric heating toward high densities and UV fields. We revised the computation of the H2 formation on grain surfaces to include the Eley-Rideal effect, thus allowing for high-temperature H2 formation. We demonstrate how the different optical properties, temperatures, and heating and cooling capabilities of the grains influence the physical and chemical structure of a model cloud. The most influential modification is the treatment of H2 formation on grain surfaces that allows for chemisorption. This increases the total H2 formation significantly and the connected H2 formation heating provides a profound heating contribution in the outer layers of the model clumps. The contribution of PAH surfaces to the photoelectric heating and H2 formation provides a boost to the temperature of outer cloud layers, which is clearly traced by high-J CO lines. Increasing the fraction of small grains in the dust size distribution results in hotter gas in the outer cloud layers caused by more efficient heating and cooler cloud centers, which is in turn caused by the more efficient FUV extinction.

Recently, we introduced detailed isotopic chemistry into the KOSMA-tau model for photon-dominated regions (PDRs) to give theoretical predictions for the abundance of the carbon isotopologues as a function of PDR parameters. Combined with radiative transfer computations for specific geometries, we estimated the possible intensity ratio of the [CII]/[13CII] lines. Here, we compare these predictions with new observations. We performed Herschel/HIFI observations of the [CII] 158micron line in a number of PDRs. In all sources we observed at least two hyperfine components of the [13CII] transition allowing to determine the [CII]/[13CII] intensity ratio, after some revision of the intrinsic hyperfine ratios. Comparing the intensity ratios with the results from the updated KOSMA-tau model, we identify cases dominated by chemical fractionation and cases dominated by the optical depth of the main isotopic line. An observable enhancement of the [CII]/[13CII] intensity ratio due to chemical fractionation depends mostly on geometry and velocity structure, and less on the gas density and radiation field. In our observations the [CII]/[13CII] ratio for the integrated line intensity was always dominated by the optical depth of the main isotopic line. However, an enhanced intensity ratio is found for particular velocity components in a few sources: the red-shifted material in the ultracompact HII region Mon R2, the wings of the turbulent profile in the Orion Bar, and possibly a blue wing in NGC7023. The mapping of the [13CII] lines in the Orion Bar allows to derive a C+ column density map confirming the temperature stratification of the C+ layer, in agreement with the chemical stratification of the Bar. The C+ column densities for all sources show that at the position of the [CII] peak emission, a dominant fraction of the gas-phase carbon is in the form of C+.

The modeling and analysis generic interface for external numerical codes (MAGIX) is a model optimizer developed under the framework of the coherent set of astrophysical tools for spectroscopy (CATS) project. The MAGIX package provides a framework of an easy interface between existing codes and an iterating engine that attempts to minimize deviations of the model results from available observational data, constraining the values of the model parameters and providing corresponding error estimates. Many models (and, in principle, not only astrophysical models) can be plugged into MAGIX to explore their parameter space and find the set of parameter values that best fits observational/experimental data. MAGIX complies with the data structures and reduction tools of ALMA (Atacama Large Millimeter Array), but can be used with other astronomical and with non-astronomical data.

Methods: We used the dual-band receiver GREAT on board the SOFIA airborne telescope to perform observations of the [C II] 158 {\mu}m fine-structure line at the postitions of two giant molecular clouds (GMC) in the center of IC 342 (GMCs C and E) and compared the spectra with corresponding ground-based data for low- and mid-J CO and [C I]. We performed model calculations assuming a clumpy photo-dissociation region (PDR) environment using the KOSMA-tau PDR model code to derive physical parameters of the local medium. Results: The [C II] 158 {\mu}m emission resembles the spectral signature of ground-based atomic and molecular lines, which indicates a common origin. The emission from GMC E can be decomposed into a cool, molecular component with weak far-ultraviolet (FUV) fields and low, mean densities of 103 cm^-3 and a strongly excited starburst/PDR region with higher densities of 104 cm^-3 and FUV intensities of 250-300 Draine fields. The emission from GMC C is consistent with gas densities of 5000 cm^-3, FUV intensities of a few Draine fields and total gas masses of 20\times10^6 M$_\odot$. Conclusions: The high spectral resolution of the GREAT receiver allowed us to decompose the [C II] emission of the GMC E into a strongly excited gas component resembling a PDR/starburst environment and a quieter, less excited gas component and to analyze the different components within a single beam individually.

Molecular globules and pillars are spectacular features, found only in the interface region between a molecular cloud and an HII-region. Impacting Far-ultraviolet (FUV) radiation creates photon dominated regions (PDRs) on their surfaces that can be traced by typical cooling lines. With the GREAT receiver onboard SOFIA we mapped and spectrally resolved the [CII] 158 micron atomic fine-structure line and the highly excited 12CO J=11-10 molecular line from three objects in Cygnus X (a pillar, a globule, and a strong IRAS source). We focus here on the globule and compare our data with existing Spitzer data and recent Herschel Open-Time PACS data. Extended [CII] emission and more compact CO-emission was found in the globule. We ascribe this emission mainly to an internal PDR, created by a possibly embedded star-cluster with at least one early B-star. However, external PDR emission caused by the excitation by the Cyg OB2 association cannot be fully excluded. The velocity-resolved [CII] emission traces the emission of PDR surfaces, possible rotation of the globule, and high-velocity outflowing gas. The globule shows a velocity shift of ~2 km/s with respect to the expanding HII-region, which can be understood as the residual turbulence of the molecular cloud from which the globule arose. This scenario is compatible with recent numerical simulations that emphazise the effect of turbulence. It is remarkable that an isolated globule shows these strong dynamical features traced by the [CII]-line, but it demands more observational studies to verify if there is indeed an embedded cluster of B-stars.

We investigate the gas dynamics and the physical properties of photodissociation regions (PDRs) in IC1396A, which is an illuminated bright-rimmed globule with internal structures created by young stellar objects. Our mapping observations of the [CII] emission in IC1396A with GREAT onboard SOFIA revealed the detailed velocity structure of this region. We combined them with observations of the [CI] 3P_1 - 3P_0 and CO(4-3) emissions to study the dynamics of the different tracers and physical properties of the PDRs. The [CII] emission generally matches the IRAC 8 micron, which traces the polycyclic aromatic hydrocarbon (PAH) emissions. The CO(4-3) emission peaks inside the globule, and the [CI] emission is strong in outer regions, following the 8 micron emission to some degree, but its peak is different from that of [CII]. The [CII] emitting gas shows a clear velocity gradient within the globule, which is not significant in the [CI] and CO(4-3) emission. Some clumps that are prominent in [CII] emission appear to be blown away from the rim of the globule. The observed ratios of [CII]/[CI] and [CII]/CO(4-3) are compared to the KOSMA-tau PDR model, which indicates a density of 10^4-10^5 cm-3.

We present an analysis of the systematic CO(2-1) survey at 12" resolution covering most of the local group spiral M 33 which, at a distance of 840 kpc, is close enough that individual giant molecular clouds (GMCs) can be identified. The goal of this work is to study the properties of the GMCs in this subsolar metallicity galaxy. The CPROPS (Cloud Properties) algorithm (Rosolowsky & Leroy 2006) was used to identify 337 GMCs in M 33, the largest sample to date in an external galaxy. The sample is used to study the GMC luminosity function, or mass spectrum under the assumption of a constant N(H2)/ICO ratio. We find that n(L)dL = K*L^(-2.0\pm0.1) for the entire sample. However, when the sample is divided into inner and outer disk samples, the exponent changes from 1.6 \pm 0.2 for the centre 2 kpc to 2.3 \pm 0.2 for galactocentric distances larger than 2 kpc. Based on the emission in the FUV, Halpha, 8mu, and 24mu bands, each cloud was classified in terms of its star forming activity - no star formation, embedded, or exposed star formation (visible in FUV and Halpha). At least one sixth of the clouds had no (massive) star formation, suggesting that the average time required for star formation to start is about one sixth of the total time for which the object is identifiable as a GMC. The clouds without star formation have significantly lower CO luminosities than those with star formation, whether embedded or exposed, presumably related to the lack of heating sources. Taking the cloud sample as a whole, the main non-trivial correlation was the decrease in cloud CO brightness (or luminosity) with galactocentric radius. The complete cloud catalog, including CO and HI spectra and the CO contours on the FUV, Halpha, 8mu, and 24mu images is presented in the appendix.

Aims. Chemical databases such as the UMIST Database for Astrochemistry (UDFA) are indispensable in the numerical modeling of astrochemical networks. Several of the listed reactions in the UDFA have properties that are problematic in numerical computations: Some are parametrized in a way that leads to extremely divergent behavior for low kinetic temperatures. Other reactions possess multiple entries that are each valid in a different temperature regime, but have no smooth transition when switching from one to another. Numerically, this introduces many difficulties.We present corrected parametrizations for these sets of reactions in the UDFA06 database. Methods. From the tabulated parametrization in UDFA, we created artificial data points and used a Levenberg-Marquardt algorithm to find a set of improved fit parameters without divergent behavior for low temperatures. For reactions with multiple entries in the database that each possess a different temperature regime, we present one joint parametrization that is designed to be valid over the whole cumulative temperature range of all individual reactions. Results. We show that it is possible to parametrize numerically problematic reactions from UDFA in a form that avoids low temperature divergence. Additionally, we demonstrate that it is possible to give a collective parametrization for reaction rate coefficients of reactions with multiple entries in UDFA. We present these new fitted values in tabulated form.

Aims: We aim at deriving the excitation conditions of the interstellar gas as well as the local FUV intensities in the molecular cloud surrounding NGC 3603 to get a coherent picture of how the gas is energized by the central stars. Methods: The NANTEN2-4m submillimeter antenna is used to map the [CI] 1-0, 2-1 and CO 4-3, 7-6 lines in a 2' x 2' region around the young OB cluster NGC 3603 YC. These data are combined with C18O 2-1 data, HIRES-processed IRAS 60 and 100 micron maps of the FIR continuum, and Spitzer/IRAC maps. Results: The NANTEN2 observations show the presence of two molecular clumps located south-east and south-west of the cluster and confirm the overall structure already found by previous CS and C18O observations. We find a slight position offset of the peak intensity of CO and [CI], and the atomic carbon appears to be further extended compared to the molecular material. We used the HIRES far-infrared dust data to derive a map of the FUV field heating the dust. We constrain the FUV field to values of \chi = 3 - 6 \times 10^3 in units of the Draine field across the clouds. Approximately 0.2 to 0.3 % of the total FUV energy is re-emitted in the [CII] 158 {\mu}m cooling line observed by ISO. Applying LTE and escape probability calculations, we derive temperatures (TMM1 = 43 K, TMM2 = 47 K), column densities (N(MM1) = 0.9 \times 10^22 cm^-2, N(MM2) = 2.5 \times 10^22 cm^-2) and densities (n(MM1) = 3 \times 10^3 cm^-3, n(MM2) = 10^3 -10^4 cm^-3) for the two observed molecular clumps MM1 and MM2. Conclusions: The cluster is strongly interacting with the ambient molecular cloud, governing its structure and physical conditions. A stability analysis shows the existence of gravitationally collapsing gas clumps which should lead to star formation. Embedded IR sources have already been observed in the outskirts of the molecular cloud and seem to support our conclusions.

We investigate the physics and chemistry of the gas and dust in dense photon-dominated regions (PDRs), along with their dependence on the illuminating UV field. Using Herschel-HIFI observations, we study the gas energetics in NGC 7023 in relation to the morphology of this nebula. NGC 7023 is the prototype of a PDR illuminated by a B2V star and is one of the key targets of Herschel. Our approach consists in determining the energetics of the region by combining the information carried by the mid-IR spectrum (extinction by classical grains, emission from very small dust particles) with that of the main gas coolant lines. In this letter, we discuss more specifically the intensity and line profile of the 158 micron (1901 GHz) [CII] line measured by HIFI and provide information on the emitting gas. We show that both the [CII] emission and the mid-IR emission from polycyclic aromatic hydrocarbons (PAHs) arise from the regions located in the transition zone between atomic and molecular gas. Using the Meudon PDR code and a simple transfer model, we find good agreement between the calculated and observed [CII] intensities. HIFI observations of NGC 7023 provide the opportunity to constrain the energetics at the surface of PDRs. Future work will include analysis of the main coolant line [OI] and use of a new PDR model that includes PAH-related species.

Using Herschel's HIFI instrument we have observed [C II] along a cut through S140 and high-J transitions of CO and HCO+ at two positions on the cut, corresponding to the externally irradiated ionization front and the embedded massive star forming core IRS1. The HIFI data were combined with available ground-based observations and modeled using the KOSMA-tau model for photon dominated regions. Here we derive the physical conditions in S140 and in particular the origin of [C II] emission around IRS1. We identify three distinct regions of [C II] emission from the cut, one close to the embedded source IRS1, one associated with the ionization front and one further into the cloud. The line emission can be understood in terms of a clumpy model of photon-dominated regions. At the position of IRS1, we identify at least two distinct components contributing to the [C II] emission, one of them a small, hot component, which can possibly be identified with the irradiated outflow walls. This is consistent with the fact that the [C II] peak at IRS1 coincides with shocked H2 emission at the edges of the outflow cavity. We note that previously available observations of IRS1 can be well reproduced by a single-component KOSMA-tau model. Thus it is HIFI's unprecedented spatial and spectral resolution, as well as its sensitivity which has allowed us to uncover an additional hot gas component in the S140 region.

We report the first detection of the ground-state rotational transition of the methylidyne cation CH+ towards the massive star-forming region DR21 with the HIFI instrument onboard the Herschel satellite. The line profile exhibits a broad emission line, in addition to two deep and broad absorption features associated with the DR21 molecular ridge and foreground gas. These observations allow us to determine a CH+ J=1-0 line frequency of 835137 +/- 3 MHz, in good agreement with a recent experimental determination. We estimate the CH+ column density to be a few 1e13 cm^-2 in the gas seen in emission, and > 1e14 cm^-2 in the components responsible for the absorption, which is indicative of a high line of sight average abundance [CH+]/[H] > 1.2x10^-8. We show that the CH+ column densities agree well with the predictions of state-of-the-art C-shock models in dense UV-illuminated gas for the emission line, and with those of turbulent dissipation models in diffuse gas for the absorption lines.

We present high resolution HIFI spectroscopy of the nucleus of the archetypical starburst galaxy M82. Six 12CO lines, 2 13CO lines and 4 fine-structure lines are detected. Besides showing the effects of the overall velocity structure of the nuclear region, the line profiles also indicate the presence of multiple components with different optical depths, temperatures and densities in the observing beam. The data have been interpreted using a grid of PDR models. It is found that the majority of the molecular gas is in low density (n=10^3.5 cm^-3) clouds, with column densities of N_H=10^21.5 cm^-2 and a relatively low UV radiation field (GO = 10^2). The remaining gas is predominantly found in clouds with higher densities (n=10^5 cm^-3) and radiation fields (GO = 10^2.75), but somewhat lower column densities (N_H=10^21.2 cm^-2). The highest J CO lines are dominated by a small (1% relative surface filling) component, with an even higher density (n=10^6 cm^-3) and UV field (GO = 10^3.25). These results show the strength of multi-component modeling for the interpretation of the integrated properties of galaxies.

The molecular gas in the DR21 massive star formation region is known to be affected by the strong UV field from the central star cluster and by a fast outflow creating a bright shock. The relative contribution of both heating mechanisms is the matter of a long debate. By better sampling the excitation ladder of various tracers we provide a quantitative distinction between the different heating mechanisms. HIFI observations of mid-J transitions of CO and HCO+ isotopes allow us to bridge the gap in excitation energies between observations from the ground, characterizing the cooler gas, and existing ISO LWS spectra, constraining the properties of the hot gas. Comparing the detailed line profiles allows to identify the physical structure of the different components. In spite of the known shock-excitation of H2 and the clearly visible strong outflow, we find that the emission of all lines up to > 2 THz can be explained by purely radiative heating of the material. However, the new Herschel/HIFI observations reveal two types of excitation conditions. We find hot and dense clumps close to the central cluster, probably dynamically affected by the outflow, and a more widespread distribution of cooler, but nevertheless dense, molecular clumps.

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.

Over the past few years several studies have provided estimates of the SFR (star-formation rate) or the total infrared luminosity from just one infrared band. However these relations are generally derived for entire galaxies, which are known to contain a large scale diffuse emission that is not necessarily related to the latest star-formation episode. We provide new relations to estimate the SFR from resolved star-forming regions at 100 mum and 160 mum. We select individual star-forming regions in the nearby (840 kpc) galaxy M33. We estimate the SFR combining the emission in Halpha and at 24 mum to calibrate the emission at 100 mum and 160 mum as SFR estimators, as mapped with PACS/Herschel. The data are obtained in the framework of the HERM33ES open time key project. There is less emission in the HII regions at 160 mum than at 100 mum. Over a dynamic range of almost 2 dex in Sigma(SFR) we find that the 100 mum emission is a nearly linear estimator of the SFR, whereas that at 160 mum is slightly superlinear. The behaviour of individual star-forming regions is surprisingly similar to that of entire galaxies. At high Sigma(SFR), star formation drives the dust temperature, whereas uncertainties and variations in radiation-transfer and dust-heated processes dominate at low Sigma(SFR). Detailed modelling of both galaxies and individual star forming regions will be needed to interpret similarities and differences between the two and assess the fraction of diffuse emission in galaxies.

In this paper we discuss the physical conditions of clumpy nature in the IC 348 molecular cloud. We combine new observations of fully sampled maps in [C I] at 492 GHz and 12CO 4--3, taken with the KOSMA 3 m telescope at about 1' resolution, with FCRAO data of 12CO 1--0, 13CO 1--0 and far-infrared continuum data observed by HIRES/IRAS. To derive the physical parameters of the region we analyze the three different line ratios. A first rough estimate of abundance is obtained from an LTE analysis. To understand the [C I] and CO emission from the PDRs in IC 348, we use a clumpy PDR model. With an ensemble of identical clumps, we constrain the total mass from the observed absolute intensities. Then we apply a more realistic clump distribution model with a power law index of 1.8 for clump-mass spectrum and a power law index of 2.3 for mass-size relation. We provide detailed fits to observations at seven representative positions in the cloud, revealing clump densities between 4 10$^{4}$ cm$^{-3}$ and 4 10$^{5}$ cm$^{-3}$ and C/CO column density ratios between 0.02 and 0.26. The derived FUV flux from the model fit is consistent with the field calculated from FIR continuum data, varying between 2 and 100 Draine units across the cloud. We find that both an ensemble of identical clumps and an ensemble with a power law clump mass distribution produce line intensities which are in good agreement (within a factor ~ 2) with the observed intensities. The models confirm the anti-correlation between the C/CO abundance ratio and the hydrogen column density found in many regions.

The fractal structure of the interstellar medium suggests that the interaction of UV radiation with the ISM as described in the context of photon-dominated regions (PDR) dominates most of the physical and chemical conditions, and hence the far-infrared and submm emission from the ISM in the Milky Way. We investigate to what extent the Galactic FIR line emission of the important species CO, C, C+, and O, as observed by the Cosmic Background Explorer (COBE) satellite can be modeled in the framework of a clumpy, UV-penetrated cloud scenario. The far-infrared line emission of the Milky Way is modeled as the emission from an ensemble of clumps with a power law clump mass spectrum and mass-size relation with power-law indices consistent with the observed ISM structure. The individual clump line intensities are calculated using the KOSMA-tau PDR-model for spherical clumps. The model parameters for the cylindrically symmetric Galactic distribution of the mass density and volume filling factor are determined by the observed radial distributions. A constant FUV intensity, in which the clumps are embedded, is assumed. We show that this scenario can explain, without any further assumptions and within a factor of about 2, the absolute FIR-line intensities and their distribution with Galactic longitude as observed by COBE.

Star formation at earlier cosmological times takes place in an interstellar medium with low metallicity. The Large Magellanic Cloud (LMC) is ideally suited to study star formation in such an environment. The physical and chemical state of the ISM in a star forming environment can be constrained by observations of submm and FIR spectral lines of the main carbon carrying species, CO, CI and CII, which originate in the surface layers of molecular clouds illuminated by the UV radiation of the newly formed, young stars. We present high-angular resolution sub-millimeter observations in the N159W region in the LMC obtained with the NANTEN2 telescope of the 12CO J = 4-3, J = 7-6, and 13CO J = 4-3 rotational and [CI] 3P1-3P0 and 3P2-3P1 fine-structure transitions. The 13CO J =4-3 and [CI] 3P2-3P1 transitions are detected for the first time in the LMC. We derive the physical and chemical properties of the low-metallicity molecular gas using an escape probability code and a self-consistent solution of the chemistry and thermal balance of the gas in the framework of a clumpy cloud PDR model. The separate excitation analysis of the submm CO lines and the carbon fine structure lines shows that the emitting gas in the N159W region has temperatures of about 80 K and densities of about 10^4 cm^-3. The estimated C to CO abundance ratio close to unity is substantially higher than in dense massive star-forming regions in the Milky Way. The analysis of all observed lines together, including the [CII] line intensity reported in the literature, in the context of a clumpy cloud PDR model constrains the UV intensity to about \chi ~220 and an average density of the clump ensemble of about 10^5 cm^-3, thus confirming the presence of high density material in the LMC N159W region.

We study the effects of a metallicity variation on the thermal balance and [CII] fine-structure line strengths in interstellar photon dominated regions (PDRs). We find that a reduction in the dust-to-gas ratio and the abundance of heavy elements in the gas phase changes the heat balance of the gas in PDRs. The surface temperature of PDRs decreases as the metallicity decreases except for high density ($n>10^6$ cm$^{-3}$) clouds exposed to weak ($\chi< 100$) FUV fields where vibrational H$_2$-deexcitation heating dominates over photoelectric heating of the gas. We incorporate the metallicity dependence in our KOSMA-$\tau$ PDR model to study the metallicity dependence of [CII]/CO line ratios in low metallicity galaxies. We find that the main trend in the variation of the observed CII/CO ratio with metallicity is well reproduced by a single spherical clump, and does not necessarily require an ensemble of clumps as in the semi-analytical model presented by Bolatto et al. (1999).