We derive molecular-gas-phase $^{12}$C/$^{13}$C isotope ratios for the central few 100 pc of the three nearby starburst galaxies NGC 253, NGC 1068, and NGC 4945 making use of the $\lambda$ $\sim$ 3 mm $^{12}$CN and $^{13}$CN $N$ = 1--0 lines in the ALMA Band 3. The $^{12}$C/$^{13}$C isotopic ratios derived from the ratios of these lines range from 30 to 67 with an average of 41.6 $\pm$ 0.2 in NGC 253, from 24 to 62 with an average of 38.3 $\pm$ 0.4 in NGC 1068, and from 6 to 44 with an average of 16.9 $\pm$ 0.3 in NGC 4945. The highest $^{12}$C/$^{13}$C isotopic ratios are determined in some of the outskirts of the nuclear regions of the three starburst galaxies. The lowest ratios are associated with the northeastern and southwestern molecular peaks of NGC 253, the northeastern and southwestern edge of the mapped region in NGC 1068, and the very center of NGC 4945. In case of NGC 1068, the measured ratios suggest inflow from the outer part of NGC 1068 into the circum-nuclear disk through both the halo and the bar. Low $^{12}$C/$^{13}$C isotopic ratios in the central regions of these starburst galaxies indicate the presence of highly processed material.

When studying the chemistry of PDRs, time dependence becomes important as visual extinction increases, since certain chemical timescales are comparable to the cloud lifetime. Dust temperature is also a key factor, since it significantly influences gas temperature and mobility on dust grains, determining the chemistry occurring on grain surfaces. We present a study of the dust temperature impact and time effects on the chemistry of different PDRs, using an updated version of the Meijerink PDR code and combining it with the time-dependent code Nahoon. We find the largest temperature effects in the inner regions of high $G$$_{\mathrm{0}}$ PDRs, where high dust temperatures favour the formation of simple oxygen-bearing molecules (especially that of O$_2$), while the formation of complex organic molecules is much more efficient at low dust temperatures. We also find that time-dependent effects strongly depend on the PDR type, since long timescales promote the destruction of oxygen-bearing molecules in the inner parts of low $G$$_{\mathrm{0}}$ PDRs, while favouring their formation and that of carbon-bearing molecules in high $G$$_{\mathrm{0}}$ PDRs. From the chemical evolution, we also conclude that, in dense PDRs, CO$_2$ is a late-forming ice compared to water ice, and confirm a layered ice structure on dust grains, with H$_2$O in lower layers than CO$_2$. Regarding steady state, the PDR edge reaches chemical equilibrium at early times ($\lesssim$10$^5$ yr). This time is even shorter ($<$10$^4$ yr) for high $G$$_{\mathrm{0}}$ PDRs. By contrast, inner regions reach equilibrium much later, especially low $G$$_{\mathrm{0}}$ PDRs, where steady state is reached at $\sim$10$^6$-10$^7$ yr.

Spectroscopic studies of ices in nearby star-forming regions indicate that ice mantles form on dust grains in two distinct steps, starting with polar ice formation (H2O rich) and switching to apolar ice (CO rich). We test how well the picture applies to more diffuse and quiescent clouds where the formation of the first layers of ice mantles can be witnessed. Medium-resolution near-infrared spectra are obtained toward background field stars behind the Pipe Nebula. The water ice absorption is positively detected at 3.0 micron in seven lines of sight out of 21 sources for which observed spectra are successfully reduced. The peak optical depth of the water ice is significantly lower than those in Taurus with the same visual extinction. The source with the highest water-ice optical depth shows CO ice absorption at 4.7 micron as well. The fractional abundance of CO ice with respect to water ice is 16+7-6 %, and about half as much as the values typically seen in low-mass star-forming regions. A small fractional abundance of CO ice is consistent with some of the existing simulations. Observations of CO2 ice in the early diffuse phase of a cloud play a decisive role in understanding the switching mechanism between polar and apolar ice formation.

The carriers of the diffuse interstellar bands (DIBs) are largely unidentified molecules ubiquitously present in the interstellar medium (ISM). After decades of study, two strong and possibly three weak near-infrared DIBs have recently been attributed to the C60+ fullerene based on observational and laboratory measurements. There is great promise for the identification of the over 400 other known DIBs, as this result could provide chemical hints towards other possible carriers. In an effort to systematically study the properties of the DIB carriers, we have initiated a new large-scale observational survey: the ESO Diffuse Interstellar Bands Large Exploration Survey (EDIBLES). The main objective is to build on and extend existing DIB surveys to make a major step forward in characterising the physical and chemical conditions for a statistically significant sample of interstellar lines-of-sight, with the goal to reverse-engineer key molecular properties of the DIB carriers. EDIBLES is a filler Large Programme using the Ultraviolet and Visual Echelle Spectrograph at the Very Large Telescope at Paranal, Chile. It is designed to provide an observationally unbiased view of the presence and behaviour of the DIBs towards early-spectral-type stars whose lines-of-sight probe the diffuse-to-translucent ISM. Such a complete dataset will provide a deep census of the atomic and molecular content, physical conditions, chemical abundances and elemental depletion levels for each sightline. Achieving these goals requires a homogeneous set of high-quality data in terms of resolution (R ~ 70000 -- 100000), sensitivity (S/N up to 1000 per resolution element), and spectral coverage (305--1042 nm), as well as a large sample size (100+ sightlines). In this first paper the goals, objectives and methodology of the EDIBLES programme are described and an initial assessment of the data is provided.

In dense and cold molecular clouds dust grains are surrounded by thick icy mantles. It is however not clear if dust growth and coagulation take place before the switch-on of a protostar. This is an important issue, as the presence of large grains may affect the chemical structure of dense cloud cores, including the dynamically important ionization fraction, and the future evolution of solids in protoplanetary disks. To study this further, we focus on L1544, one of the most centrally concentrated pre-stellar cores on the verge of star formation, and with a well-known physical structure. We observed L1544 at 1.2 and 2 mm using NIKA, a new receiver at the IRAM 30 m telescope, and we used data from the Herschel Space Observatory archive. We find no evidence of grain growth towards the center of L1544 at the available angular resolution. Therefore, we conclude that single dish observations do not allow us to investigate grain growth toward the pre-stellar core L1544 and high sensitivity interferometer observations are needed. We predict that dust grains can grow to 200 $\mu$m in size toward the central ~300 au of L1544. This will imply a dust opacity change by a factor of ~2.5 at 1.2 mm, which can be detected using the Atacama Large Millimeter and submillimeter Array (ALMA) at different wavelengths and with an angular resolution of 2".

The temperature of interstellar dust particles is of great importance to astronomers. It plays a crucial role in the thermodynamics of interstellar clouds, because of the gas-dust collisional coupling. It is also a key parameter in astrochemical studies that governs the rate at which molecules form on dust. In 3D (magneto)hydrodynamic simulations often a simple expression for the dust temperature is adopted, because of computational constraints, while astrochemical modelers tend to keep the dust temperature constant over a large range of parameter space. Our aim is to provide an easy-to-use parametric expression for the dust temperature as a function of visual extinction ($A_{\rm V}$) and to shed light on the critical dependencies of the dust temperature on the grain composition. We obtain an expression for the dust temperature by semi-analytically solving the dust thermal balance for different types of grains and compare to a collection of recent observational measurements. We also explore the effect of ices on the dust temperature. Our results show that a mixed carbonaceous-silicate type dust with a high carbon volume fraction matches the observations best. We find that ice formation allows the dust to be warmer by up to 15% at high optical depths ($A_{\rm V}> 20$ mag) in the interstellar medium. Our parametric expression for the dust temperature is presented as $T_{\rm d} = \left[ 11 + 5.7\times \tanh\bigl( 0.61 - \log_{10}(A_{\rm V})\bigr) \right] \, \chi_{\rm uv}^{1/5.9}$, where $\chi_{\rm uv}$ is in units of the Draine (1978) UV field

We present ALMA detections of the [CI] 1-0, CO J=3-2, and CO J=4-3 emission lines, as well as the ALMA band 4 continuum for a compact star-forming galaxy (cSFG) at z=2.225, 3D-HST GS30274. As is typical for cSFGs, this galaxy has a stellar mass of $1.89 \pm 0.47\,\times 10^{11}\,\rm{M}_\odot$, with a star formation rate of $214\pm44\,\rm{M}_\odot\,\rm{yr}^{-1}$ putting it on the star-forming `main-sequence', but with an H-band effective radius of 2.5 kpc, making it much smaller than the bulk of `main-sequence' star-forming galaxies. The intensity ratio of the line detections yield an ISM density (~ 6 $\times 10^{4}\,\rm{cm}^{-3}$) and a UV-radiation field ( ~2 $\times 10^4\,\rm{G}_0$), similar to the values in local starburst and ultra-luminous infrared galaxy environments. A starburst phase is consistent with the short depletion times ($t_{\rm H2, dep} \leq 140$ Myr) we find using three different proxies for the H2 mass ([CI], CO, dust mass). This depletion time is significantly shorter than in more extended SFGs with similar stellar masses and SFRs. Moreover, the gas fraction of 3D-HST GS30274 is smaller than typically found in extended galaxies. We measure the CO and [CI] kinematics and find a FWHM line width of ~$750 \pm 41 $ km s$^{-1}$. The CO and [CI] FWHM are consistent with a previously measured H$\alpha$ FWHM for this source. The line widths are consistent with gravitational motions, suggesting we are seeing a compact molecular gas reservoir. A previous merger event, as suggested by the asymmetric light profile, may be responsible for the compact distribution of gas and has triggered a central starburst event. This event gives rise to the starburst-like ISM properties and short depletion times. The centrally located and efficient star formation is quickly building up a dense core of stars, responsible for the compact distribution of stellar light in 3D-HST GS30274.

We present a multi-wavelength analysis of 52 sub-millimeter galaxies (SMGs), identified using ALMA 870$\mu$m continuum imaging in a pilot program to precisely locate bright SCUBA2-selected sub-mm sources in the UKIDSS Ultra Deep Survey (UDS) field. Using the available deep (especially near-infrared), panoramic imaging of the UDS field at optical-to-radio wavelengths we characterize key properties of the SMG population. The median photometric redshift of the bright ALMA/SCUBA-2 UDS (AS2UDS) SMGs that are detected in a sufficient number of wavebands to derive a robust photometric redshift is $z$=2.65$\pm$0.13. However, similar to previous studies, 27% of the SMGs are too faint at optical-to-near-infrared wavelengths to derive a reliable photometric redshift. Assuming that these SMGs lie at z$\gtrsim$3 raises the median redshift of the full sample to $z$=2.9$\pm$0.2. A subset of 23, unlensed, bright AS2UDS SMGs have sizes measured from resolved imaging of their rest-frame far-infrared emission. We show that the extent and luminosity of the far-infrared emission are consistent with the dust emission arising from regions that are optically thick, on average, at a wavelength of $\lambda_0$$\ge$75$\mu$m (1-$\sigma$ dispersion of 55-90$\mu$m). Using the dust masses derived from our optically-thick spectral energy distribution models we determine that these galaxies have a median hydrogen column density of $N_{H}$=9.8$_{-0.7}^{+1.4}$$\times$10$^{23}$cm$^{-2}$, or a corresponding median $V$-band obscuration of $A_\mathrm{v}$=540$^{+80}_{-40}$mag, averaged along the line of sight to the source of their restframe $\sim$200$\mu$m emission. We discuss the implications of this extreme attenuation by dust for the multiwavelength study of dusty starbursts and reddening-sensitive tracers of star formation.

The physical properties of brown dwarf discs, in terms of their shapes and sizes, are still largely unexplored by observations. To what extent brown dwarf discs are similar to scaled-down T Tauri discs is currently unknown, and this work is a step towards establishing a relationship through the eventual modelling of future observations. We use observations of the brown dwarf disc $\rho$ Oph 102 to infer a fiducial model around which we build a small grid of brown dwarf disc models, in order to model the CO, HCN, and HCO+ line fluxes and the chemistry which drives their abundances. These are the first brown dwarf models to be published which relate detailed, 2D radiation thermochemical disc models to observational data. We predict that moderately extended ALMA antenna configurations will spatially resolve CO line emission around brown dwarf discs, and that HCN and HCO+ will be detectable in integrated flux, following our conclusion that the flux ratios of these molecules to CO emission are comparable to that of T Tauri discs. These molecules have not yet been observed in sub-mm wavelengths in a brown dwarf disc, yet they are crucial tracers of the warm surface-layer gas and of ionization in the outer parts of the disc. We present the prediction that if the physical and chemical processes in brown dwarf discs are similar to those that occur in T Tauri discs -- as our models suggest -- then the same diagnostics that are used for T Tauri discs can be used for brown dwarf discs (such as HCN and HCO+ lines that have not yet been observed in the sub-mm), and that these lines should be observable with ALMA. Through future observations, either confirmation (or refutation) of these ideas about brown dwarf disc chemistry will have strong implications for our understanding of disc chemistry, structure, and subsequent planet formation in brown dwarf discs.

Mid-IR water lines from protoplanetary disks around T Tauri stars have a detection rate of 50\%. Models have identified multiple physical properties of disks such as dust-to-gas mass ratio, dust size power law distribution, disk gas mass, disk inner radius, and disk scale height as potential explanation for the current detection rate. We search for a connection between mid-IR water line fluxes and the strength of the 10~$\mu$m silicate feature. We analyse observed water line fluxes from three blends and compute the 10~$\mu$m silicate feature strength from Spitzer spectra. We use a series of published models, exploring disk dust and gas properties, and the effects of different stars. The models also show that the increasing stellar luminosity enhance simultaneously the strength of this dust feature and the water lines fluxes. No correlation is found between the observed mid-IR water lines and the 10~$\mu$m silicate. Our sample shows the same difference in the peak strength between amorphous and crystalline silicates that was noted in earlier studies, but our models do not support this intrinsic difference in silicate peak strength. Individual properties of our model series are not able to reproduce the most extreme observations, suggesting that more complex dust properties are required. A parametrized settling prescription is able to boost the peak strength by a factor 2 for the standard model. Water line fluxes are unrelated to the composition of the dust. The pronounced regular trends seen in the model results are washed out in the data due to the larger diversity in stellar and disk properties compared to our model series. The disks with with weaker mid-IR water line fluxes are depleted in gas or enhanced in dust in the inner 10~au. In the case of gas depleted disks, settling produces very strong 10~$\mu$m silicate features, with strong peak strength.

We present deep observations at 450 um and 850 um in the Extended Groth Strip field taken with the SCUBA-2 camera mounted on the James Clerk Maxwell Telescope as part of the deep SCUBA-2 Cosmology Legacy Survey (S2CLS), achieving a central instrumental depth of $\sigma_{450}=1.2$ mJy/beam and $\sigma_{850}=0.2$ mJy/beam. We detect 57 sources at 450 um and 90 at 850 um with S/N > 3.5 over ~70 sq. arcmin. From these detections we derive the number counts at flux densities $S_{450}>4.0$ mJy and $S_{850}>0.9$ mJy, which represent the deepest number counts at these wavelengths derived using directly extracted sources from only blank-field observations with a single-dish telescope. Our measurements smoothly connect the gap between previous shallower blank-field single-dish observations and deep interferometric ALMA results. We estimate the contribution of our SCUBA-2 detected galaxies to the cosmic infrared background (CIB), as well as the contribution of 24 um-selected galaxies through a stacking technique, which add a total of $0.26\pm0.03$ and $0.07\pm0.01$ MJy/sr, at 450 um and 850 um, respectively. These surface brightnesses correspond to $60\pm20$ and $50\pm20$ per cent of the total CIB measurements, where the errors are dominated by those of the total CIB. Using the photometric redshifts of the 24 um-selected sample and the redshift distributions of the submillimetre galaxies, we find that the redshift distribution of the recovered CIB is different at each wavelength, with a peak at $z\sim1$ for 450 um and at $z\sim2$ for 850um, consistent with previous observations and theoretical models.

LkCa 15 is an extensively studied star in the Taurus region known for its pre-transitional disk with a large inner cavity in dust continuum and normal gas accretion rate. The most popular hypothesis to explain the LkCa 15 data invokes one or more planets to carve out the inner cavity, while gas continues to flow across the gap from the outer disk onto the central star. We present spatially unresolved HCO+ J=4-3 observations of the LkCa 15 disk from the JCMT and model the data with the ProDiMo code. We find that: (1) HCO+ line-wings are clearly detected, certifying the presence of gas in the cavity within <50 AU of the star. (2) Reproducing the observed line-wing flux requires both a significant suppression of cavity dust (by a factor >10^4 compared to the ISM) and a substantial increase in the gas scale-height within the cavity (H_0/R_0 ~ 0.6). An ISM dust-to-gas ratio (d:g=10^-2) yields too little line-wing flux regardless of the scale-height or cavity gas geometry, while a smaller scale-height also under predicts the flux even with a reduced d:g. (3) The cavity gas mass is consistent with the surface density profile of the outer disk extended inwards to the sublimation radius (corresponding to mass M_d ~ 0.03 M_sun), and masses lower by a factor >10 appear to be ruled out.

(abridged) We report rest-frame submillimeter H2O emission line observations of 11 HyLIRGs/ULIRGs at z~2-4 selected among the brightest lensed galaxies discovered in the Herschel-ATLAS. Using the IRAM NOEMA, we have detected 14 new H2O emission lines. The apparent luminosities of the H2O emission lines are $\mu L_{\rm{H_2O}} \sim 6-21 \times 10^8 L_\odot$, with velocity-integrated line fluxes ranging from 4-15 Jy km s$^{-1}$. We have also observed CO emission lines using EMIR on the IRAM 30m telescope in seven sources. The velocity widths for CO and H2O lines are found to be similar. With almost comparable integrated flux densities to those of the high-J CO line, H2O is found to be among the strongest molecular emitters in high-z Hy/ULIRGs. We also confirm our previously found correlation between luminosity of H2O ($L_{\rm{H_2O}}$) and infrared ($L_{\rm{IR}}$) that $L_{\rm{H_2O}} \sim L_{\rm{IR}}^{1.1-1.2}$, with our new detections. This correlation could be explained by a dominant role of far-infrared (FIR) pumping in the H2O excitation. Modelling reveals the FIR radiation fields have warm dust temperature $T_\rm{warm}$~45-75 K, H2O column density per unit velocity interval $N_{\rm{H_2O}}/\Delta V \gtrsim 0.3 \times 10^{15}$ cm$^{-2}$ km$^{-1}$ s and 100 $\mu$m continuum opacity $\tau_{100} > 1$ (optically thick), indicating that H2O is likely to trace highly obscured warm dense gas. However, further observations of $J\geq4$ H2O lines are needed to better constrain the continuum optical depth and other physical conditions of the molecular gas and dust. We have also detected H2O+ emission in three sources. A tight correlation between $L_{\rm{H_2O}}$ and $L_{\rm{H_2O^+}}$ has been found in galaxies from low to high redshift. The velocity-integrated flux density ratio between H2O+ and H2O suggests that cosmic rays generated by strong star formation are possibly driving the H2O+ formation.

We apply photo-dissociation region (PDR) molecular line emission models, that have varying degrees of enhanced mechanical heating rates, to the gaseous component of simulations of star-forming galaxies taken from the literature. Snapshots of these simulations are used to produce line emission maps for the rotational transitions of the CO molecule and its 13CO isotope up to J = 4-3. We consider two galaxy models: a small disk galaxy of solar metallicity and a lighter dwarf galaxy with 0.2 \zsun metallicity. Elevated excitation temperatures for CO(1 - 0) correlate positively with mechanical feedback, that is enhanced towards the central region of both model galaxies. The emission maps of these model galaxies are used to compute line ratios of CO and 13CO transitions. These line ratios are used as diagnostics where we attempt to match them These line ratios are used as diagnostics where we attempt to match them to mechanically heated single component (i.e. uniform density, Far-UV flux, visual extinction and velocity gradient) equilibrium PDR models. We find that PDRs ignoring mechanical feedback in the heating budget over-estimate the gas density by a factor of 100 and the far-UV flux by factors of ~10 - 1000. In contrast, PDRs that take mechanical feedback into account are able to fit all the line ratios for the central < 2 kpc of the fiducial disk galaxy quite well. The mean mechanical heating rate per H atom that we recover from the line ratio fits of this region varies between $10^{-27}$ -- $10^{-26}$~erg s$^{-1}$.

Far-infrared molecular emission is an important tool used to understand the excitation mechanisms of the gas in the inter-stellar medium of star-forming galaxies. In the present work, we model the emission from rotational transitions with critical densities n >~ 10^4 cm-3. We include 4-3 < J <= 15-14 transitions of CO and 13CO, in addition to J <= 7-6 transitions of HCN, HNC, and HCO+ on galactic scales. We do this by re-sampling high density gas in a hydrodynamic model of a gas-rich disk galaxy, assuming that the density field of the interstellar medium of the model galaxy follows the probability density function (PDF) inferred from the resolved low density scales. We find that in a narrow gas density PDF, with a mean density of ~10 cm-3 and a dispersion \sigma = 2.1 in the log of the density, most of the emission of molecular lines, emanates from the 10-1000 cm-3 part of the PDF. We construct synthetic emission maps for the central 2 kpc of the galaxy and fit the line ratios of CO and 13CO up to J = 15-14, as well as HCN, HNC, and HCO+ up to J = 7-6, using one photo-dissociation region (PDR) model. We attribute the goodness of the one component fits for our model galaxy to the fact that the distribution of the luminosity, as a function of density, is peaked at gas densities between 10 and 1000 cm-3. We explore the impact of different log-normal density PDFs on the distribution of the line-luminosity as a function of density, and we show that it is necessary to have a broad dispersion, corresponding to Mach numbers >~ 30 in order to obtain significant emission from n > 10^4 cm-3 gas. Such Mach numbers are expected in star-forming galaxies, LIRGS, and ULIRGS. By fitting line ratios of HCN(1-0), HNC(1-0), and HCO+(1-0) for a sample of LIRGS and ULIRGS using mechanically heated PDRs, we constrain the Mach number of these galaxies to 29 < M < 77.

We present new Atacama Large Millimeter/Submillimeter Array (ALMA) 850um continuum observations of the original Lyman-alpha Blob (LAB) in the SSA22 field at z=3.1 (SSA22-LAB01). The ALMA map resolves the previously identified submillimeter source into three components with total flux density S_850 = 1.68+/-0.06 mJy, corresponding to a star formation rate of ~150 M_sun/yr. The submillimeter sources are associated with several faint (m~27 mag) rest-frame ultraviolet sources identified in Hubble Space Telescope Imaging Spectrograph (STIS) clear filter imaging (~5850A). One of these companions is spectroscopically confirmed with Keck MOSFIRE to lie within 20 projected kpc and 250 km/s of one of the ALMA components. We postulate that some of these STIS sources represent a population of low-mass star-forming satellites surrounding the central submillimeter sources, potentially contributing to their growth and activity through accretion. Using a high resolution cosmological zoom simulation of a 10^13 M_sun halo at z=3, including stellar, dust and Ly-alpha radiative transfer, we can model the ALMA+STIS observations and demonstrate that Ly-alpha photons escaping from the central submillimeter sources are expected to resonantly scatter in neutral hydrogen, the majority of which is predicted to be associated with halo substructure. We show how this process gives rise to extended Ly-alpha emission with similar surface brightness and morphology to observed giant LABs.

We present a catalogue of nearly 3,000 submillimetre sources detected at 850um over ~5 square degrees surveyed as part of the James Clerk Maxwell Telescope (JCMT) SCUBA-2 Cosmology Legacy Survey (S2CLS). This is the largest survey of its kind at 850um, probing a meaningful cosmic volume at the peak of star formation activity and increasing the sample size of submillimetre galaxies selected at 850um by an order of magnitude. We describe the wide 850um survey component of S2CLS, which covers the key extragalactic survey fields: UKIDSS-UDS, COSMOS, Akari-NEP, Extended Groth Strip, Lockman Hole North, SSA22 and GOODS-North. The average 1-sigma depth of S2CLS is 1.2 mJy/beam, approaching the SCUBA-2 850um confusion limit, which we determine to be ~0.8 mJy/beam. We measure the single dish 850um number counts to unprecedented accuracy, reducing the Poisson errors on the differential counts to approximately 4% at S_850~3mJy. With several independent fields, we investigate field-to-field variance, finding that the number counts on 0.5-1 degree scales are generally within 50% of the S2CLS mean for S_850>3mJy, with scatter consistent with the Poisson and estimated cosmic variance uncertainties, although there is a marginal (2-sigma) density enhancement in the GOODS-North field. The observed number counts are in reasonable agreement with recent phenomenological and semi-analytic models. Finally, the large solid angle of S2CLS allows us to measure the bright-end counts: at S_850>10mJy there are approximately ten sources per square degree, and we detect the distinctive up-turn in the number counts indicative of the detection of local sources of 850um emission and strongly lensed high-redshift galaxies. Here we describe the data collection and reduction procedures and present calibrated maps and a catalogue of sources; these are made publicly available.

Now that ALMA is reaching its full capabilities, observations of sub-mm emission line deep fields become feasible. We couple a semi-analytic model of galaxy formation with a radiative transfer code to make predictions for the luminosity function of CO J=1-0 up to CO J=6-5 and [CII] at redshifts z=0-6. We find that: 1) our model correctly reproduces the CO and [CII] emission of low- and high-redshift galaxies and reproduces the available constraints on the CO luminosity function at z<2.75; 2) we find that the CO and [CII] luminosity functions of galaxies increase from z = 6 to z = 4, remain relatively constant till z = 1 and rapidly decrease towards z = 0. The galaxies that are brightest in CO and [CII] are found at z~2; 3) the CO J=3-2 emission line is most favourable to study the CO luminosity and global H2 mass content of galaxies, because of its brightness and observability with currently available sub-mm and radio instruments; 4) the luminosity functions of high-J CO lines show stronger evolution than the luminosity functions of low-J CO lines; 5) our model barely reproduces the available constraints on the CO and [CII] luminosity function of galaxies at z>1.5 and the CO luminosity of individual galaxies at intermediate redshifts. We argue that this is driven by a lack of cold gas in galaxies at intermediate redshifts as predicted by cosmological simulations of galaxy formation.

Ly-alpha blobs (LABs) offer insight into the complex interface between galaxies and their circumgalactic medium. Whilst some LABs have been found to contain luminous star-forming galaxies and active galactic nuclei that could potentially power the Ly-alpha emission, others appear not to be associated with obvious luminous galaxy counterparts. It has been speculated that LABs may be powered by cold gas streaming on to a central galaxy, providing an opportunity to directly observe the `cold accretion' mode of galaxy growth. Star-forming galaxies in LABs could be dust obscured and therefore detectable only at longer wavelengths. We stack deep SCUBA-2 observations of the SSA22 field to determine the average 850um flux density of 34 LABs. We measure S_850 = 0.6 +/- 0.2mJy for all LABs, but stacking the LABs by size indicates that only the largest third (area > 1794 kpc^2) have a mean detection, at 4.5 sigma, with S_850 = 1.4 +/- 0.3mJy. Only two LABs (1 and 18) have individual SCUBA-2 > 3.5 sigma detections at a depth of 1.1mJy/beam. We consider two possible mechanisms for powering the LABs and find that central star formation is likely to dominate the emission of Ly-alpha, with cold accretion playing a secondary role.

The presence of dust can strongly affect the chemical composition of the interstellar medium. We model the chemistry in photodissociation regions (PDRs) using both gas-phase and dust-phase chemical reactions. Our aim is to determine the chemical compositions of the interstellar medium (gas/dust/ice) in regions with distinct (molecular) gas densities that are exposed to radiation fields with different intensities. We have significantly improved the Meijerink PDR code by including 3050 new gas-phase chemical reactions and also by implementing surface chemistry. In particular, we have included 117 chemical reactions occurring on grain surfaces covering different processes, such as adsorption, thermal desorption, chemical desorption, two-body reactions, photo processes, and cosmic-ray processes on dust grains. We obtain abundances for different gas and solid species as a function of visual extinction, depending on the density and radiation field. We also analyse the rates of the formation of CO2 and H2O ices in different environments. In addition, we study how chemistry is affected by the presence/absence of ice mantles (bare dust or icy dust) and the impact of considering different desorption probabilities. The type of substrate (bare dust or icy dust) and the probability of desorption can significantly alter the chemistry occurring on grain surfaces, leading to differences of several orders of magnitude in the abundances of gas-phase species, such as CO, H2CO, and CH3OH. The type of substrate, together with the density and intensity of the radiation field, also determine the threshold extinction to form ices of CO2 and H2O. We also conclude that H2CO and CH3OH are mainly released into the gas phase of low, far-ultraviolet illuminated PDRs through chemical desorption upon two-body surface reactions, rather than through photodesorption.

The understanding of hydrogen attachment to carbonaceous surfaces is essential to a wide variety of research fields and technologies such as hydrogen storage for transportation, precise localization of hydrogen in electronic devices and the formation of cosmic H2. For coronene cations as prototypical Polycyclic Aromatic Hydrocarbon (PAH) molecules, the existence of magic numbers upon hydrogenation was uncovered experimentally. Quantum chemistry calculations show that hydrogenation follows a site-specific sequence leading to the appearance of cations having 5, 11, or 17 hydrogen atoms attached, exactly the magic numbers found in the experiments. For these closed-shell cations, further hydrogenation requires appreciable structural changes associated with a high transition barrier. Controlling specific hydrogenation pathways would provide the possibility to tune the location of hydrogen attachment and the stability of the system. The sequence to hydrogenate PAHs, leading to PAHs with magic numbers of H atoms attached, provides clues to understand that carbon in space is mostly aromatic and partially aliphatic in PAHs. PAH hydrogenation is fundamental to assess the contribution of PAHs to the formation of cosmic H2.

In the earliest phases of star-forming clouds, stable molecular species, such as CO, are important coolants in the gas phase. Depletion of these molecules on dust surfaces affects the thermal balance of molecular clouds and with that their whole evolution. For the first time, we study the effect of grain surface chemistry (GSC) on star formation and its impact on the initial mass function (IMF). We follow a contracting translucent cloud in which we treat the gas-grain chemical interplay in detail, including the process of freeze-out. We perform 3d hydrodynamical simulations under three different conditions, a pure gas-phase model, a freeze-out model, and a complete chemistry model. The models display different thermal evolution during cloud collapse. The equation of state (EOS) of the gas becomes softer with CO freeze-out and the results show that at the onset of star formation, the cloud retains its evolution history such that the number of formed stars differ (by 7%) between the three models. While the stellar mass distribution results in a different IMF when we consider pure freeze-out, with the complete treatment of the GSC, the divergence from a pure gas-phase model is minimal. We find that the impact of freeze-out is balanced by the non-thermal processes; chemical and photodesorption. We also find an average filament width of 0.12 pc ($\pm$0.03 pc), and speculate that this may be a result from the changes in the EOS caused by the gas-dust thermal coupling. We conclude that GSC plays a big role in the chemical composition of molecular clouds and that surface processes are needed to accurately interpret observations, however, that GSC does not have a significant impact as far as star formation and the IMF is concerned.

We propose a set of standard assumptions for the modelling of Class II and III protoplanetary disks, which includes detailed continuum radiative transfer, thermo-chemical modelling of gas and ice, and line radiative transfer from optical to cm wavelengths. We propose new standard dust opacities for disk models, we present a simplified treatment of PAHs sufficient to reproduce the PAH emission features, and we suggest using a simple treatment of dust settling. We roughly adjust parameters to obtain a model that predicts typical Class II T Tauri star continuum and line observations. We systematically study the impact of each model parameter (disk mass, disk extension and shape, dust settling, dust size and opacity, gas/dust ratio, etc.) on all continuum and line observables, in particular on the SED, mm-slope, continuum visibilities, and emission lines including [OI] 63um, high-J CO lines, (sub-)mm CO isotopologue lines, and CO fundamental ro-vibrational lines. We find that evolved dust properties (large grains) often needed to fit the SED, have important consequences for disk chemistry and heating/cooling balance, leading to stronger emission lines in general. Strong dust settling and missing disk flaring have similar effects on continuum observations, but opposite effects on far-IR gas emission lines. PAH molecules can shield the gas from stellar UV radiation because of their strong absorption and negligible scattering opacities. The observable millimetre-slope of the SED can become significantly more gentle in the case of cold disk midplanes, which we find regularly in our T Tauri models. We propose to use line observations of robust chemical tracers of the gas, such as O, CO, and H2, as additional constraints to determine some key properties of the disks, such as disk shape and mass, opacities, and the dust/gas ratio, by simultaneously fitting continuum and line observations.

The mid-IR detection rate of water lines in disks around Herbig stars disks is about 5\%, while it is around 50\% for disks around TTauri stars. The reason for this is still unclear. In this study, we want to find an explanation for the different detection rates between low mass and high mass pre-main-sequence stars (PMSs) in the mid-IR regime. We run disk models with stellar parameters adjusted to spectral types B9 through M2, using the radiation thermo-chemical disk modeling code ProDiMo. We produce convolved spectra at the resolution of Spitzer IRS, JWST MIRI and VLT VISIR spectrographs. We apply random noise derived from typical Spitzer spectra for a direct comparison with observations. The strength of the mid-IR water lines correlates directly with the luminosity of the central star. We explored a small parameter space around a standard disk model, considering dust-to-gas mass ratio, disk gas mass, mixing coefficient for dust settling, flaring index, dust maximum size and size power law distribution index. The models show that it is possible to suppress the water emission, however, current observations are not sensitive enough to detect mid-IR lines in disks for most of the explored parameters. The presence of noise in the spectra, combined with the high continuum flux (noise level is proportional to the continuum flux), is the most likely explanation for the non detections towards Herbig stars. Mid-IR spectra with resolution higher than 20000 are needed to investigate water in protoplanetary disks. Intrinsic differences in disk structure, e.g. inner gaps, gas-to-dust ratio, dust size and distribution, and inner disk scale height, between Herbig and TTauri star disks are able to explain a lower water detection rate in disks around Herbig stars.

We investigate which properties of protoplanetary disks around T Tauri stars affect the physics and chemistry in the regions where mid- and far-IR water lines originate and their respective line fluxes. We search for diagnostics for future observations. With the code ProDiMo, we build a series of models exploring a large parameter space, computing rotational and rovibrational transitions of water in nonlocal thermodynamic equilibrium (non-LTE). We select a sample of transitions in the mid- IR regime and the fundamental ortho and para water transitions in the far-IR. We investigate the chemistry and the local physical conditions in the line emitting regions. We calculate Spitzer spectra for each model and compare far-IR and mid-IR lines. In addition, we use mid-IR colors to tie the water line predictions to the dust continuum. Parameters affecting the water line fluxes in disks by more than a factor of three are : the disk gas mass, the dust-to-gas mass ratio, the dust maximum grain size, ISM(InterStellarMedium) UV radiation field, the mixing parameter of Dubrulle settling, the disk flaring parameter, and the dust size distribution. The first four parameters affect the mid-IR lines much more than the far-IR lines. A key driver behind water spectroscopy is the dust opacity, which sets the location of the water line emitting region. We identify three types of parameters. Parameters, such as dust-to-gas ratio, ISM radiation field, and dust size distribution, affect the mid-IR lines more, while the far-IR transitions are more affected by the flaring index. The gas mass greatly affects lines in both regimes. Higher spectral resolution and line sensitivities, like from the James Webb Space Telescope, are needed to detect a statistically relevant sample of individual water lines to distinguish further between these types of parameters.

We present early results from the JCMT Plane Survey (JPS), which has surveyed the northern inner Galactic plane between longitudes l=7 and l=63 degrees in the 850-{\mu}m continuum with SCUBA-2, as part of the James Clerk Maxwell Telescope Legacy Survey programme. Data from the l=30 degree survey region, which contains the massive star-forming regions W43 and G29.96, are analysed after approximately 40% of the observations had been completed. The pixel-to-pixel noise is found to be 19 mJy/beam, after a smooth over the beam area, and the projected equivalent noise levels in the final survey are expected to be around 10 mJy/beam. An initial extraction of compact sources was performed using the FellWalker method resulting in the detection of 1029 sources above a 5-{\sigma} surface-brightness threshold. The completeness limits in these data are estimated to be around 0.2 Jy/beam (peak flux density) and 0.8 Jy (integrated flux density) and are therefore probably already dominated by source confusion in this relatively crowded section of the survey. The flux densities of extracted compact sources are consistent with those of matching detections in the shallower ATLASGAL survey. We analyse the virial and evolutionary state of the detected clumps in the W43 star-forming complex and find that they appear younger than the Galactic-plane average.

We probe the chemical and energetic conditions in dense gas created by radiative feedback through observations of multiple CO, HCN and HCO$^+$ transitions toward the dense core of M17 SW. We used the dual band receiver GREAT on board the SOFIA airborne telescope to obtain maps of the $J=16-15$, $J=12-11$, and $J=11-10$ transitions of $^{12}$CO. We compare these maps with corresponding APEX and IRAM 30m telescope data for low- and mid-$J$ CO, HCN and HCO$^+$ emission lines, including maps of the HCN $J=8-7$ and HCO$^+$ $J=9-8$ transitions. The excitation conditions of $^{12}$CO, HCO$^+$ and HCN are estimated with a two-phase non-LTE radiative transfer model of the line spectral energy distributions (LSEDs) at four selected positions. The energy balance at these positions is also studied. We obtained extensive LSEDs for the CO, HCN and HCO$^+$ molecules toward M17 SW. The LSED shape, particularly the high-$J$ tail of the CO lines observed with SOFIA/GREAT, is distinctive for the underlying excitation conditions. The critical magnetic field criterion implies that the cold cloudlets at two positions are partially controlled by processes that create and dissipate internal motions. Supersonic but sub-Alfv\'enic velocities in the cold component at most selected positions indicates that internal motions are likely MHD waves. Magnetic pressure dominates thermal pressure in both gas components at all selected positions, assuming random orientation of the magnetic field. The magnetic pressure of a constant magnetic field throughout all the gas phases can support the total internal pressure of the cold components, but it cannot support the internal pressure of the warm components. If the magnetic field scales as $B \propto n^{2/3}$, then the evolution of the cold cloudlets at two selected positions, and the warm cloudlets at all selected positions, will be determined by ambipolar diffusion.

We present high resolution (0."4) IRAM PdBI and ALMA mm and submm observations of the (ultra) luminous infrared galaxies ((U)LIRGs) IRAS17208-0014, Arp220, IC860 and Zw049.057 that reveal intense line emission from vibrationally excited ($\nu_2$=1) J=3-2 and 4-3 HCN. The emission is emerging from buried, compact (r<17-70 pc) nuclei that have very high implied mid-infrared surface brightness $>$$5\times 10^{13}$ L$_{\odot}$ kpc$^{-2}$. These nuclei are likely powered by accreting supermassive black holes (SMBHs) and/or hot (>200 K) extreme starbursts. Vibrational, $\nu_2$=1, lines of HCN are excited by intense 14 micron mid-infrared emission and are excellent probes of the dynamics, masses, and physical conditions of (U)LIRG nuclei when H$_2$ column densities exceed $10^{24}$ cm$^{-2}$. It is clear that these lines open up a new interesting avenue to gain access to the most obscured AGNs and starbursts. Vibrationally excited HCN acts as a proxy for the absorbed mid-infrared emission from the embedded nuclei, which allows for reconstruction of the intrinsic, hotter dust SED. In contrast, we show strong evidence that the ground vibrational state ($\nu$=0), J=3-2 and 4-3 rotational lines of HCN and HCO$^+$ fail to probe the highly enshrouded, compact nuclear regions owing to strong self- and continuum absorption. The HCN and HCO$^+$ line profiles are double-peaked because of the absorption and show evidence of non-circular motions - possibly in the form of in- or outflows. Detections of vibrationally excited HCN in external galaxies are so far limited to ULIRGs and early-type spiral LIRGs, and we discuss possible causes for this. We tentatively suggest that the peak of vibrationally excited HCN emission is connected to a rapid stage of nuclear growth, before the phase of strong feedback.

We obtained an ALMA Cycle 0 spectral scan of the dusty LIRG NGC 4418, spanning a total of 70.7 GHz in bands 3, 6, and 7. We use a combined local thermal equilibrium (LTE) and non-LTE (NLTE) fit of the spectrum in order to identify the molecular species and derive column densities and excitation temperatures. We derive molecular abundances and compare them with other Galactic and extragalactic sources by means of a principal component analysis. We detect 317 emission lines from a total of 45 molecular species, including 15 isotopic substitutions and six vibrationally excited variants. Our LTE/NLTE fit find kinetic temperatures from 20 to 350 K, and densities between 10$^5$ and 10$^7$ cm$^{-3}$. The spectrum is dominated by vibrationally excited HC$_3$N, HCN, and HNC, with vibrational temperatures from 300 to 450 K. We find high abundances of HC$_3$N, SiO, H$_2$S, and c-HCCCH and a low CH$_3$OH abundance. A principal component analysis shows that NGC 4418 and Arp 220 share very similar molecular abundances and excitation, which clearly set them apart from other Galactic and extragalactic environments. The similar molecular abundances observed towards NCG 4418 and Arp 220 are consistent with a hot gas-phase chemistry, with the relative abundances of SiO and CH$_3$OH being regulated by shocks and X-ray driven dissociation. The bright emission from vibrationally excited species confirms the presence of a compact IR source, with an effective diameter $<$5 pc and brightness temperatures $>$350 K. The molecular abundances and the vibrationally excited spectrum are consistent with a young AGN/starburst system. We suggest that NGC 4418 may be a template for a new kind of chemistry and excitation, typical of compact obscured nuclei (CON). Because of the narrow line widths and bright molecular emission, NGC 4418 is the ideal target for further studies of the chemistry in CONs.

Molecular hydrogen is the most abundant molecule in the Universe. It is thought that a large portion of H2 forms by association of hydrogen atoms to polycyclic aromatic hydrocarbons (PAHs). We model the influence of PAHs on total H2 formation rates in photodissociation regions (PDRs) and assess the effect of these formation rates on the total cloud structure. We set up a chemical kinetic model at steady state in a PDR environment and included adiative transfer to calculate the chemistry at different depths in the PDR. This model includes known dust grain chemistry for the formation of H2 and a H2 formation mechanism on PAHs. Since H2 formation on PAHs is impeded by thermal barriers, this pathway is only efficient at higher temperatures (T > 200 K). At these temperatures the conventional route of H2 formation via H atoms physisorbed on dust grains is no longer feasible, so the PAH mechanism enlarges the region where H2 formation is possible. We find that PAHs have a significant influence on the structure of PDRs. The extinction at which the transition from atomic to molecular hydrogen occurs strongly depends on the presence of PAHs, especially for PDRs with a strong external radiation field. A sharp spatial transition between fully dehydrogenated PAHs on the outside of the cloud and normally hydrogenated PAHs on the inside is found. As a proof of concept, we use coronene to show that H2 forms very efficiently on PAHs, and that this process can reproduce the high H2 formation rates derived in several PDRs.

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

We probe the column densities and masses traced by the ionized and neutral atomic carbon with spectrally resolved maps, and compare them to the diffuse and dense molecular gas traced by [C I] and low-$J$ CO lines toward the star-forming region M17SW. We mapped a 4.1pc x 4.7pc region in the [C I] 609 m$\mu$ line using the APEX telescope, as well as the CO isotopologues with the IRAM 30m telescope. We analyze the data based on velocity channel maps that are 1 km/s wide. We correlate their spatial distribution with that of the [C II] map obtained with SOFIA/GREAT. Optically thin approximations were used to estimate the column densities of [C I] and [C II] in each velocity channel. The spatial distribution of the [C I] and all CO isotopologues emission was found to be associated with that of [C II] in about 20%-80% of the mapped region, with the high correlation found in the central (15-23 km/s ) velocity channels. The excitation temperature of [C I] ranges between 40 K and 100 K in the inner molecular region of M17 SW. Column densities in 1 km/s channels between ~10$^{15}$ and ~10$^{17}$ cm$^{-2}$ were found for [C I]. Just ~20% of the velocity range (~40 km/s) that the [C II] line spans is associated with the star-forming material traced by [C I] and CO. The total gas mass estimated from the [C II] emission gives a lower limit of ~4.4x10$^3$ $M_{\odot}$. At least 64% of this mass is not associated with the star-forming material in M17SW. We also found that about 36%, 17%, and 47% of the [C II] emission is associated with the HII, HI, and H_2 regimes, respectively. Comparisons with the H41$\alpha$ line shows an ionization region mixed with the neutral and part of the molecular gas, in agreement with the clumped structure and dynamical processes at play in M17SW. These results are also relevant to extra-galactic studies in which [C II] is often used as a tracer of star-forming material.

Herschel/PACS observations of 29 local (Ultra-)Luminous Infrared Galaxies, including both starburst and AGN-dominated sources as diagnosed in the mid-infrared/optical, show that the equivalent width of the absorbing OH 65 um Pi_{3/2} J=9/2-7/2 line (W_{eq}(OH65)) with lower level energy E_{low}~300 K, is anticorrelated with the [C ii]158 um line to far-infrared luminosity ratio, and correlated with the far-infrared luminosity per unit gas mass and with the 60-to-100 um far-infrared color. While all sources are in the active L_{IR}/M_{H2}>50 Lsun/Msun mode as derived from previous CO line studies, the OH65 absorption shows a bimodal distribution with a discontinuity at L_{FIR}/M_{H2}~100 Lsun/Msun. In the most buried sources, OH65 probes material partially responsible for the silicate 9.7 um absorption. Combined with observations of the OH 71 um Pi_{1/2} J=7/2-5/2 doublet (E_{low}~415 K), radiative transfer models characterized by the equivalent dust temperature, Tdust, and the continuum optical depth at 100 um, tau_{100}, indicate that strong [C ii]158 um deficits are associated with far-IR thick (tau_{100}>~0.7, N_{H}>~10^{24} cm^{-2}), warm (Tdust>~60 K) structures where the OH 65 um absorption is produced, most likely in circumnuclear disks/tori/cocoons. With their high L_{FIR}/M_{H2} ratios and columns, the presence of these structures is expected to give rise to strong [C ii] deficits. W_{eq}(OH65) probes the fraction of infrared luminosity arising from these compact/warm environments, which is >~30-50% in sources with high W_{eq}({OH65}). Sources with high W_{eq}({OH65}) have surface densities of both L_{IR} and M_{H2} higher than inferred from the half-light (CO or UV/optical) radius, tracing coherent structures that represent the most buried/active stage of (circum)nuclear starburst-AGN co-evolution.

We present high-resolution (0.3'') ALMA 870um imaging of 52 sub-millimeter galaxies (SMGs) in the Ultra Deep Survey (UDS) field and investigate the size and morphology of the sub-millimeter (sub-mm) emission on 2-10kpc scales. We derive a median intrinsic angular size of FWHM=0.30$\pm$0.04'' for the 23 SMGs in the sample detected at a signal-to-noise ratio (SNR) >10. Using the photometric redshifts of the SMGs we show that this corresponds to a median physical half-light diameter of 2.4$\pm$0.2kpc. A stacking analysis of the SMGs detected at an SNR <10 shows they have sizes consistent with the 870um-bright SMGs in the sample. We compare our results to the sizes of SMGs derived from other multi-wavelength studies, and show that the rest-frame ~250um sizes of SMGs are consistent with studies of resolved 12CO (J=3-2 to 7-6) emission lines, but that sizes derived from 1.4GHz imaging appear to be approximately two times larger on average, which we attribute to cosmic ray diffusion. The rest-frame optical sizes of SMGs are around four times larger than the sub-millimeter sizes, indicating that the star formation in these galaxies is compact relative to the pre-existing stellar distribution. The size of the starburst region in SMGs is consistent with the majority of the star formation occurring in a central region, a few kpc in extent, with a median star formation rate surface density of 90$\pm$30Msol/yr/kpc$^2$, which may suggest that we are witnessing an intense period of bulge growth in these galaxies.

CO, $^{13}$CO and C$^{18}$O ${\it J}$ = 3--2 observations are presented of the Ophiuchus molecular cloud. The $^{13}$CO and C$^{18}$O emission is dominated by the Oph A clump, and the Oph B1, B2, C, E, F and J regions. The optically thin(ner) C$^{18}$O line is used as a column density tracer, from which the gravitational binding energy is estimated to be $4.5 \times 10^{39}$ J (2282 $M_\odot$ km$^2$ s$^{-2}$). The turbulent kinetic energy is $6.3 \times 10^{38}$ J (320 $M_\odot$ km$^2$ s$^{-2}$), or 7 times less than this, and therefore the Oph cloud as a whole is gravitationally bound. Thirty protostars were searched for high velocity gas, with eight showing outflows, and twenty more having evidence of high velocity gas along their lines-of-sight. The total outflow kinetic energy is $1.3 \times 10^{38}$ J (67 $M_\odot$ km$^2$ s$^{-2}$), corresponding to 21$\%$ of the cloud's turbulent kinetic energy. Although turbulent injection by outflows is significant, but does ${\it not}$ appear to be the dominant source of turbulence in the cloud. 105 dense molecular clumplets were identified, which had radii $\sim$ 0.01--0.05 pc, virial masses $\sim$ 0.1--12 $M_\odot$, luminosities $\sim$ 0.001--0.1 K~km s$^{-1}$ pc$^{-2}$, and excitation temperatures $\sim$ 10--50K. These are consistent with the standard GMC based size-line width relationships, showing that the scaling laws extend down to size scales of hundredths of a parsec, and to sub solar-mass condensations. There is however no compelling evidence that the majority of clumplets are undergoing free-fall collapse, nor that they are pressure confined.

Observations of high redshift quasars at $z>6$ indicate that they harbor supermassive black holes (SMBHs) of a billion solar masses. The direct collapse scenario has emerged as the most plausible way to assemble SMBHs. The nurseries for the direct collapse black holes are massive primordial halos illuminated with an intense UV flux emitted by population II (Pop II) stars. In this study, we compute the critical value of such a flux ($J_{21}^{\rm crit}$) for realistic spectra of Pop II stars through three-dimensional cosmological simulations. We derive the dependence of $J_{21}^{\rm crit}$ on the radiation spectra, on variations from halo to halo, and on the impact of X-ray ionization. Our findings show that the value of $J_{21}^{\rm crit}$ is a few times $\rm 10^4$ and only weakly depends on the adopted radiation spectra in the range between $T_{\rm rad}=2 \times 10^4-10^5$ K. For three simulated halos of a few times $\rm 10^{7}$~M$_{\odot}$, $J_{21}^{\rm crit}$ varies from $\rm 2 \times 10^4 - 5 \times 10^4$. The impact of X-ray ionization is almost negligible and within the expected scatter of $J_{21}^{\rm crit}$ for background fluxes of $J_{\rm X,21} \leq 0.1$. The computed estimates of $J_{21}^{\rm crit}$ have profound implications for the quasar abundance at $z=10$ as it lowers the number density of black holes forming through an isothermal direct collapse by a few orders of magnitude below the observed black holes density. However, the sites with moderate amounts of $\rm H_2$ cooling may still form massive objects sufficient to be compatible with observations.

Context. The seeds of the first supermassive black holes may have resulted from the direct collapse of hot primordial gas in $\gtrsim 10^4$ K haloes, forming a supermassive or quasistar as an intermediate stage. Aims. We explore the formation of a protostar resulting from the collapse of primordial gas in the presence of a strong Lyman-Werner radiation background. Particularly, we investigate the impact of turbulence and rotation on the fragmentation behaviour of the gas cloud. We accomplish this goal by varying the initial turbulent and rotational velocities. Methods. We performed 3D adaptive mesh refinement simulations with a resolution of 64 cells per Jeans length using the ENZO code, simulating the formation of a protostar up to unprecedentedly high central densities of $10^{21}$ cm$^{-3}$, and spatial scales of a few solar radii. To achieve this goal, we employed the KROME package to improve modelling of the chemical and thermal processes. Results. We find that the physical properties of the simulated gas clouds become similar on small scales, irrespective of the initial amount of turbulence and rotation. After the highest level of refinement was reached, the simulations have been evolved for an additional ~5 freefall times. A single bound clump with a radius of $2 \times 10^{-2}$ AU and a mass of ~$7 \times 10^{-2}$ M$_{\odot}$ is formed at the end of each simulation, marking the onset of protostar formation. No strong fragmentation is observed by the end of the simulations, regardless of the initial amount of turbulence or rotation, and high accretion rates of a few solar masses per year are found. Conclusions. Given such high accretion rates, a quasistar of $10^5$ M$_{\odot}$ is expected to form within $10^5$ years.

In the last decade, the growth of supermassive black holes (SMBHs) has been intricately linked to galaxy formation and evolution and is a key ingredient in the assembly of galaxies. To investigate the origin of SMBHs, we perform cosmological simulations that target the direct collapse black hole (DCBH) seed formation scenario in the presence of two different strong Lyman-Werner (LW) background fields. These simulations include the X-ray irradiation from a central massive black hole (MBH), $\rm{H}_2$ self-shielding and stellar feedback from metal-free and metal-enriched stars. We find in both simulations that local X-ray feedback induces metal-free star formation $\sim 0.5$ Myr after the MBH forms. The MBH accretion rate reaches a maximum of $10^{-3}$ $M_{\odot}$ yr$^{-1}$ in both simulations. However, the duty cycle differs which is derived to be $6\%$ and $50\%$ for high and low LW cases, respectively. The MBH in the high LW case grows only $\sim 6\%$ in 100 Myr compared to $16\%$ in the low LW case. We find that the maximum accretion rate is determined by the local gas thermodynamics whereas the duty cycle is determined by the large scale gas dynamics and gas reservoir. We conclude that radiative feedback from the central MBH plays an important role in star formation in the nuclear regions and stifling initial MBH growth, relative to the typical Eddington rate argument, and that initial MBH growth might be affected by the local LW radiation field.

We combine a semi-analytic model of galaxy formation, tracking atomic and molecular phases of cold gas, with a three-dimensional radiative-transfer and line tracing code to study the sub-mm emission from atomic and molecular species (CO, HCN, [CI], [CII], [OI]) in galaxies. We compare the physics that drives the formation of stars at the epoch of peak star formation (SF) in the Universe (z = 2.0) with that in local galaxies. We find that normal star-forming galaxies at high redshift have much higher CO-excitation peaks than their local counterparts and that CO cooling takes place at higher excitation levels. CO line ratios increase with redshift as a function of galaxy star-formation rate, but are well correlated with H2 surface density independent of redshift. We find an increase in the [OI]/[CII] line ratio in typical star-forming galaxies at z = 1.2 and z = 2.0 with respect to counterparts at z = 0. Our model results suggest that typical star-forming galaxies at high redshift consist of much denser and warmer star-forming clouds than their local counterparts. Galaxies belonging to the tail of the SF activity peak at z = 1.2 are already less dense and cooler than counterparts during the actual peak of SF activity (z = 2.0). We use our results to discuss how future ALMA surveys can best confront our predictions and constrain models of galaxy formation.

We report the detection at 850um of the central source in SSA22-LAB1, the archetypal Lyman-alpha Blob (LAB), a 100kpc-scale radio-quiet emission-line nebula at z=3.1. The flux density of the source, $S_{850}=4.6\pm1.1$mJy implies the presence of a galaxy, or group of galaxies, with a total luminosity of $L_{\rm IR}\approx10^{12}L_\odot$. The position of an active source at the center of a ~50kpc-radius ring of linearly polarized Ly-alpha emission detected by Hayes et al. (2011) suggests that the central source is leaking Ly-alpha photons preferentially in the plane of the sky, which undergo scattering in HI clouds at large galactocentric radius. The Ly-alpha morphology around the submillimeter detection is reminiscent of biconical outflow, and the average Ly-alpha line profiles of the two `lobes' are dominated by a red peak, expected for a resonant line emerging from a medium with a bulk velocity gradient that is outflowing relative to the line center. Taken together, these observations provide compelling evidence that the central active galaxy (or galaxies) is responsible for a large fraction of the extended Ly-alpha emission and morphology. Less clear is the history of the cold gas in the circumgalactic medium being traced by Ly-alpha: is it mainly pristine material accreting into the halo that has not yet been processed through an interstellar medium (ISM), now being blown back as it encounters an outflow, or does it mainly comprise gas that has been swept-up within the ISM and expelled from the galaxy?

We present a detailed analysis of ALMA Bands 7 and 9 data of CO, HCO+, HCN and CS, augmented with Plateau de Bure Interferometer (PdBI) data of the ~ 200 pc circumnuclear disk (CND) and the ~ 1.3 kpc starburst ring (SB ring) of NGC~1068, a nearby (D = 14 Mpc) Seyfert 2 barred galaxy. We aim at determining the physical characteristics of the dense gas present in the CND and whether the different line intensity ratios we find within the CND as well as between the CND and the SB ring are due to excitation effects (gas density and temperature differences) or to a different chemistry. We estimate the column densities of each species in Local Thermodynamic Equilibrium (LTE). We then compute large one-dimensional non-LTE radiative transfer grids (using RADEX) by using first only the CO transitions, and then all the available molecules in order to constrain the densities, temperatures and column densities within the CND. We finally present a preliminary set of chemical models to determine the origin of the gas. We find that in general the gas in the CND is very dense (> 10^5 cm^-3) and hot (T> 150K), with differences especially in the temperature across the CND. The AGN position has the lowest CO/HCO+, CO/HCN and CO/CS column density ratios. RADEX analyses seem to indicate that there is chemical differentiation across the CND. We also find differences between the chemistry of the SB ring and some regions of the CND; the SB ring is also much colder and less dense than the CND. Chemical modelling does not succeed in reproducing all the molecular ratios with one model per region, suggesting the presence of multi-gas phase components. The LTE, RADEX and chemical analyses all indicate that more than one gas-phase component is necessary to uniquely fit all the available molecular ratios within the CND.

We present FIR-CO luminosity relations ($\log L_{\rm FIR} = \alpha \log L'_{\rm CO} + \beta$) for the full CO rotational ladder from J=1-0 to J=13-12 for 62 local (z < 0.1) (Ultra) Luminous Infrared Galaxies (LIRGs) using data from Herschel SPIRE-FTS and ground-based telescopes. We extend our sample to high redshifts (z > 1) by including 35 (sub)-millimeter selected dusty star forming galaxies from the literature with robust CO observations. The addition of luminous starbursts at high redshifts enlarge the range of the FIR-CO luminosity relations towards the high-IR-luminosity end while also significantly increasing the small amount of mid-/high-J CO line data available prior to Herschel. This new data-set (both in terms of IR luminosity and J-ladder) reveals linear FIR-CO luminosity relations ($\alpha \sim 1$) for J=1-0 up to J=5-4, with a nearly constant normalisation ($\beta \sim 2$). This is expected from the (also) linear FIR-(molecular line) relations found for the dense gas tracer lines (HCN and CS), as long as the dense gas mass fraction does not vary strongly within our (merger/starburst)-dominated sample. However from J=6-5 and up to J=13-12 we find an increasingly sub-linear slope and higher normalization constant with increasing J. We argue that these are caused by a warm (~100K) and dense ($>10^4{\rm cm^{-3}}$) gas component whose thermal state is unlikely to be maintained by star formation powered far-UV radiation fields (and thus is no longer directly tied to the star formation rate). We suggest that mechanical heating (e.g., supernova driven turbulence and shocks), and not cosmic rays, is the more likely source of energy for this component. The global CO spectral line energy distributions (SLEDs), which remain highly excited from J=6-5 up to J=13-12, are found to be a generic feature of the (U)LIRGs in our sample, and further support the presence of this gas component.

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants

The ability of metal free gas to cool by molecular hydrogen in primordial halos is strongly associated with the strength of ultraviolet (UV) flux produced by the stellar populations in the first galaxies. Depending on the stellar spectrum, these UV photons can either dissociate $\rm H_{2}$ molecules directly or indirectly by photo-detachment of $\rm H^{-}$ as the latter provides the main pathway for $\rm H_{2}$ formation in the early universe. In this study, we aim to determine the critical strength of the UV flux above which the formation of molecular hydrogen remains suppressed for a sample of five distinct halos at $z>10$ by employing a higher order chemical solver and a Jeans resolution of 32 cells. We presume that such flux is emitted by PopII stars implying atmospheric temperatures of $\rm 10^{4}$~K. We performed three-dimensional cosmological simulations and varied the strength of the UV flux below the Lyman limit in units of $\rm J_{21}$. Our findings show that the value of $\rm J_{21}^{crit}$ varies from halo to halo and is sensitive to the local thermal conditions of the gas. For the simulated halos it varies from 400-700 with the exception of one halo where $\rm J_{21}^{crit} \geq 1500$. This has important implications for the formation of direct collapse black holes and their estimated population at z > 6. It reduces the number density of direct collapse black holes by almost three orders of magnitude compared to the previous estimates.

We investigate the fueling and the feedback of star formation and nuclear activity in NGC1068, a nearby (D=14Mpc) Seyfert 2 barred galaxy, by analyzing the distribution and kinematics of the molecular gas in the disk. We have used ALMA to map the emission of a set of dense molecular gas tracers (CO(3-2), CO(6-5), HCN(4-3), HCO+(4-3) and CS(7-6)) and their underlying continuum emission in the central r ~ 2kpc of NGC1068 with spatial resolutions ~ 0.3"-0.5" (~ 20-35pc). Molecular line and dust continuum emissions are detected from a r ~ 200pc off-centered circumnuclear disk (CND), from the 2.6kpc-diameter bar region, and from the r ~ 1.3kpc starburst (SB) ring. Most of the emission in HCO+, HCN and CS stems from the CND. Molecular line ratios show dramatic order-of-magnitude changes inside the CND that are correlated with the UV/X-ray illumination by the AGN, betraying ongoing feedback. The gas kinematics from r ~ 50pc out to r ~ 400pc reveal a massive (M_mol ~ 2.7 (+0.9, -1.2) x 10^7 Msun) outflow in all molecular tracers. The tight correlation between the ionized gas outflow, the radio jet and the occurrence of outward motions in the disk suggests that the outflow is AGN-driven. The outflow rate estimated in the CND, dM/dt ~ 63 (+21, -37) Msun yr^-1, is an order of magnitude higher than the star formation rate at these radii, confirming that the outflow is AGN-driven. The power of the AGN is able to account for the estimated momentum and kinetic luminosity of the outflow. The CND mass load rate of the CND outflow implies a very short gas depletion time scale of <=1 Myr.

Radiative feedback from populations II stars played a vital role in early structure formation. Particularly, photons below the Lyman limit can escape the star forming regions and produce a background ultraviolet (UV) flux which consequently may influence the pristine halos far away from the radiation sources. These photons can quench the formation of molecular hydrogen by photo-detachment of $\rm H^{-}$. In this study, we explore the impact of such UV radiation on fragmentation in massive primordial halos of a few times $\rm 10^{7}$~M${_\odot}$. To accomplish this goal, we perform high resolution cosmological simulations for two distinct halos and vary the strength of the impinging background UV field in units of $\rm J_{21}$. We further make use of sink particles to follow the evolution for 10,000 years after reaching the maximum refinement level. No vigorous fragmentation is observed in UV illuminated halos while the accretion rate changes according to the thermal properties. Our findings show that a few 100-10, 000 solar mass protostars are formed when halos are irradiated by $\rm J_{21}=10-500$ at $\rm z>10$ and suggest a strong relation between the strength of UV flux and mass of a protostar. This mode of star formation is quite different from minihalos, as higher accretion rates of about $\rm 0.01-0.1$ M$_{\odot}$/yr are observed by the end of our simulations. The resulting massive stars are the potential cradles for the formation of intermediate mass black holes at earlier cosmic times and contribute to the formation of a global X-ray background.

We report Plateau de Bure Interferometer (PdBI) 1.1 mm continuum imaging towards two extremely red H-[4.5]>4 (AB) galaxies at z>3, which we have previously discovered making use of Spitzer SEDS and Hubble Space Telescope CANDELS ultra-deep images of the Ultra Deep Survey field. One of our objects is detected on the PdBI map with a 4.3 sigma significance, corresponding to Snu(1.1mm)=(0.78 +/- 0.18) mJy. By combining this detection with the Spitzer 8 and 24 micron photometry for this source, and SCUBA2 flux density upper limits, we infer that this galaxy is a composite active galactic nucleus/star-forming system. The infrared (IR)-derived star formation rate is SFR~(200 +/- 100) Msun/yr, which implies that this galaxy is a higher-redshift analogue of the ordinary ultra-luminous infrared galaxies (ULIRGs) more commonly found at z~2-3. In the field of the other target, we find a tentative 3.1 sigma detection on the PdBI 1.1 mm map, but 3.7 arcsec away of our target position, so it likely corresponds to a different object. In spite of the lower significance, the PdBI detection is supported by a close SCUBA2 3.3 sigma detection. No counterpart is found on either the deep SEDS or CANDELS maps, so, if real, the PdBI source could be similar in nature to the sub-millimetre source GN10. We conclude that the analysis of ultra-deep near- and mid-IR images offers an efficient, alternative route to discover new sites of powerful star formation activity at high redshifts.

We report on two regularly rotating galaxies at redshift z=2, using high resolution spectra of the bright [CII] 158 micron emission line from the HIFI instrument on the Herschel Space Observatory. Both SDSS090122.37+181432.3 ("S0901") and SDSS J120602.09+514229.5 ("the Clone") are strongly lensed and show the double-horned line profile that is typical of rotating gas disks. Using a parametric disk model to fit the emission line profiles, we find that S0901 has a rotation speed v sin(i) = 120 +/- 7 km/s and gas velocity dispersion sigma < 23 km/s. The best fitting model for the Clone is a rotationally supported disk having v sin(i) = 79 +/- 11 km/s and sigma < 4km/s. However the Clone is also consistent with a family of dispersion-dominated models having sigma = 92 +/- 20 km/s. Our results showcase the potential of the [CII] line as a kinematic probe of high redshift galaxy dynamics: [CII] is bright; accessible to heterodyne receivers with exquisite velocity resolution; and traces dense star-forming interstellar gas. Future [CII] line observations with ALMA would offer the further advantage of spatial resolution, allowing a clearer separation between rotation and velocity dispersion.

CO observations in active galactic nuclei and star-bursts reveal high kinetic temperatures. Those environments are thought to be very turbulent due to dynamic phenomena such as outflows and high supernova rates. We investigate the effect of mechanical heating (MH) on atomic fine-structure and molecular lines, and their ratios. We use those ratios as a diagnostic to constrain the amount of MH in an object and also study its significance on estimating the H2 mass. Equilibrium PDRs models were used to compute the thermal and chemical balance for the clouds. The equilibria were solved for numerically using the optimized version of the Leiden PDR-XDR code. Large velocity gradient calculations were done as post-processing on the output of the PDR models using RADEX. High-J CO line ratios are very sensitive to MH. Emission becomes at least one order of magnitude brighter in clouds with n~10^5~cm^-3 and a star formation rate of 1 Solar Mass per year (corresponding to a MH rate of 2 * 10^-19 erg cm^-3 s^-1). Emission of low-J CO lines is not as sensitive to MH, but they do become brighter in response to MH. Generally, for all of the lines we considered, MH increases excitation temperatures and decreases the optical depth at the line centre. Hence line ratios are also affected, strongly in some cases. Ratios involving HCN are a good diagnostic for MH, such as HCN(1-0)/CO(1-0) and HCN(1-0)/HCO^+(1-0). Both ratios increase by a factor 3 or more for a MH equivalent to > 5 percent of the surface heating, as opposed to pure PDRs. The first major conclusion is that low-J to high-J intensity ratios will yield a good estimate of the MH rate (as opposed to only low-J ratios). The second one is that the MH rate should be taken into account when determining A_V or equivalently N_H, and consequently the cloud mass. Ignoring MH will also lead to large errors in density and radiation field estimates.

The structure of protoplanetary disks is thought to be linked to the temperature and chemistry of their dust and gas. Whether the disk is flat or flaring depends on the amount of radiation that it absorbs at a given radius, and on the efficiency with which this is converted into thermal energy. The understanding of these heating and cooling processes is crucial to provide a reliable disk structure for the interpretation of dust continuum emission and gas line fluxes. Especially in the upper layers of the disk, where gas and dust are thermally decoupled, the infrared line emission is strictly related to the gas heating/cooling processes. We aim to study the thermal properties of the disk in the oxygen line emission region, and to investigate the relative importance of X-ray (1-120 Angstrom) and far-UV radiation (FUV, 912-2070 Angstrom) for the heating balance there. We use [OI] 63 micron line fluxes observed in a sample of protoplanetary disks of the Taurus/Auriga star forming region and compare it to the model predictions presented in our previous work. The data were obtained with the PACS instrument on board the Herschel Space Observatory as part of the Herschel Open Time Key Program GASPS (GAS in Protoplanetary diskS), published in Howard et al. (2013). Our theoretical grid of disk models can reproduce the [OI] absolute fluxes and predict a correlation between [OI] and the sum Lx+Lfuv. The data show no correlation between the [OI] line flux and the X-ray luminosity, the FUV luminosity or their sum. The data show that the FUV or X-ray radiation has no notable impact on the region where the [OI] line is formed. This is in contrast with what is predicted from our models. Possible explanations are that the disks in Taurus are less flaring than the hydrostatic models predict, and/or that other disk structure aspects that were left unchanged in our models are important. ..abridged..

We present spectrally resolved Herschel/HIFI observations of the young multiple system T Tau in atomic and molecular lines. While CO, H2O, [C II], and SO lines trace the envelope and the outflowing gas up to velocities of 33 km/s with respect to systemic, the CN 5-4 hyperfine structure lines at 566.7, 566.9 GHz show a narrow double-peaked profile centered at systemic velocity, consistent with an origin in the outer region of the compact disk of T Tau N. Disk modeling of the T Tau N disk with the thermo-chemical code ProDiMo produces CN line fluxes and profiles consistent with the observed ones and constrain the size of the gaseous disk (R_out = 110 (+10, -20) AU) and its inclination (i = 25 \pm 5 degree). The model indicates that the CN lines originate in a disk upper layer at 40-110 AU from the star, which is irradiated by the stellar UV field and heated up to temperatures of 50-700 K. With respect to previously observed CN 2-1 millimeter lines, the CN 5-4 lines appear to be less affected by envelope emission, due to their larger critical density and excitation temperature. Hence, high-J CN lines are a unique confusion-free tracer of embedded disks, such as the disk of T Tau N.

Atoms and molecules, and in particular CO, are important coolants during the evolution of interstellar star-forming gas clouds. The presence of dust grains, which allow many chemical reactions to occur on their surfaces, strongly impacts the chemical composition of a cloud. At low temperatures, dust grains can lock-up species from the gas phase which freeze out and form ices. In this sense, dust can deplete important coolants. Our aim is to understand the effects of freeze-out on the thermal balance and the evolution of a gravitationally bound molecular cloud. For this purpose, we perform 3D hydrodynamical simulations with the adaptive mesh code FLASH. We simulate a gravitationally unstable cloud under two different conditions, with and without grain surface chemistry. We let the cloud evolve until one free-fall time is reached and track the thermal evolution and the abundances of species during this time. We see that at a number density of 10$^4$ cm$^{-3}$ most of the CO molecules are frozen on dust grains in the run with grain surface chemistry, thereby depriving the most important coolant. As a consequence, we find that the temperature of the gas rises up to $\sim$25 K. The temperature drops once again due to gas-grain collisional cooling when the density reaches a few$\times$10$^4$ cm$^{-3}$. We conclude that grain surface chemistry not only affects the chemical abundances in the gas phase, but also leaves a distinct imprint in the thermal evolution that impacts the fragmentation of a star-forming cloud. As a final step, we present the equation of state of a collapsing molecular cloud that has grain surface chemistry included.

The effects of Wheeler's quantum foam on black hole growth are explored from an astrophysical perspective. Quantum fluctuations in the form of mini (10^-5 g) black holes can couple to macroscopic black holes and allow the latter to grow exponentially in mass on a time scale of ~10^9 years. Consequently, supermassive black holes can acquire a lot of their mass through these quantum contributions over the life time of the universe. This alleviates the need for very efficient forms of baryonic matter accretion more recent than a redshift z~6. Sgr A* in the Milky Way center is a candidate to verify this quantum space-time effect, with a predicted mass growth rate of 4x10^-3 Mo yr^-1. A few comments on the possibility and consequences of dark matter as quantum grown black holes are made, with a big crunch fate of the universe.

At the distance of 99-116 pc, HD141569A is one of the nearest HerbigAe stars that is surrounded by a tenuous disk, probably in transition between a massive primordial disk and a debris disk. We observed the fine-structure lines of OI at 63 and 145 micron and the CII line at 157 micron with the PACS instrument onboard the Herschel Space Telescope as part of the open-time large programme GASPS. We complemented the atomic line observations with archival Spitzer spectroscopic and photometric continuum data, a ground-based VLT-VISIR image at 8.6 micron, and 12CO fundamental ro-vibrational and pure rotational J=3-2 observations. We simultaneously modeled the continuum emission and the line fluxes with the Monte Carlo radiative transfer code MCFOST and the thermo-chemical code ProDiMo to derive the disk gas- and dust properties assuming no dust settling. The models suggest that the oxygen lines are emitted from the inner disk around HD141569A, whereas the [CII] line emission is more extended. The CO submillimeter flux is emitted mostly by the outer disk. Simultaneous modeling of the photometric and line data using a realistic disk structure suggests a dust mass derived from grains with a radius smaller than 1 mm of 2.1E-7 MSun and from grains with a radius of up to 1 cm of 4.9E-6 MSun. We constrained the polycyclic aromatic hydrocarbons (PAH) mass to be between 2E-11 and 1..4E-10 MSun assuming circumcircumcoronene (C150H30) as the representative PAH. The associated PAH abundance relative to hydrogen is lower than those found in the interstellar medium (3E-7) by two to three orders of magnitude. The disk around HD141569A is less massive in gas (2.5 to 4.9E-4 MSun or 67 to 164 MEarth) and has a flat opening angle (<10%). [abridged]

We investigate the multi-wavelength properties of a sample of 450-\mu m selected sources from the SCUBA-2 Cosmology Legacy Survey (S2CLS). A total of 69 sources were identified above 4\sigma\ in deep SCUBA-2 450-\mu m observations overlapping the UDS and COSMOS fields and covering 210 sq. arcmin to a typical depth of \sigma 450=1.5 mJy. Reliable cross identification are found for 58 sources (84 per cent) in Spitzer and Hubble Space Telescope WFC3/IR data. The photometric redshift distribution (dN/dz) of 450\mu m-selected sources is presented, showing a broad peak in the redshift range 1<z<3, and a median of z=1.4. Combining the SCUBA-2 photometry with Herschel SPIRE data from HerMES, the submm spectral energy distribution (SED) is examined via the use of modified blackbody fits, yielding aggregate values for the IR luminosity, dust temperature and emissivity of <LIR>=10^12 +/- 0.8 L_sol, <T_D>=42 +/- 11 K and <\beta_D>=1.6 +/- 0.5, respectively. The relationship between these SED parameters and the physical properties of galaxies is investigated, revealing correlations between T_D and LIR and between \beta_D and both stellar mass and effective radius. The connection between star formation rate and stellar mass is explored, with 24 per cent of 450 \mu m sources found to be ``star-bursts'', i.e. displaying anomalously high specific SFRs. However, both the number density and observed properties of these ``star-burst'' galaxies are found consistent with the population of normal star-forming galaxies.

Context. This paper discusses the sensitivity of water lines to chemical processes and radiative transfer for the protoplanetary disk around TW Hya. The study focuses on the Herschel spectral range in the context of new line detections with the PACS instrument from the Gas in Protoplanetary Systems project (GASPS). Aims. The paper presents an overview of the chemistry in the main water reservoirs in the disk around TW Hya. It discusses the limitations in the interpretation of observed water line fluxes. Methods. ... (abbreviated) Results. We report new line detections of p-H2O (3_22-2_11) at 89.99 micron and CO J=18-17 at 144.78 micron for the disk around TW Hya. Disk modeling shows that the far-IR fine structure lines ([OI], [CII]) and molecular submm lines are very robust to uncertainties in the chemistry, while the water line fluxes can change by factors of a few. The water lines are optically thick, sub-thermally excited and can couple to the background continuum radiation field. The low-excitation water lines are also sensitive to uncertainties in the collision rates, e.g. with neutral hydrogen. The gas temperature plays an important role for the [OI] fine structure line fluxes, the water line fluxes originating from the inner disk as well as the high excitation CO, CH+ and OH lines. Conclusions. Due to their sensitivity on chemical input data and radiative transfer, water lines have to be used cautiously for understanding details of the disk structure. Water lines covering a wide range of excitation energies provide access to the various gas phase water reservoirs (inside and outside the snow line) in protoplanetary disks and thus provide important information on where gas-phase water is potentially located. Experimental and/or theoretical collision rates for H2O with atomic hydrogen are needed to diminish uncertainties from water line radiative transfer.

The growth of luminous structures and the building blocks of life in the Universe began as primordial gas was processed in stars and mixed at galactic scales. The mechanisms responsible for this development are not well understood and have changed over the intervening 13 billion years. To follow the evolution of matter over cosmic time, it is necessary to study the strongest (resonance) transitions of the most abundant species in the Universe. Most of them are in the ultraviolet (UV; 950A-3000A) spectral range that is unobservable from the ground. A versatile space observatory with UV sensitivity a factor of 50-100 greater than existing facilities will revolutionize our understanding of the Universe. Habitable planets grow in protostellar discs under ultraviolet irradiation, a by-product of the star-disk interaction that drives the physical and chemical evolution of discs and young planetary systems. The electronic transitions of the most abundant molecules are pumped by the UV field, providing unique diagnostics of the planet-forming environment that cannot be accessed from the ground. Earth's atmosphere is in constant interaction with the interplanetary medium and the solar UV radiation field. A 50-100 times improvement in sensitivity would enable the observation of the key atmospheric ingredients of Earth-like exoplanets (carbon, oxygen, ozone), provide crucial input for models of biologically active worlds outside the solar system, and provide the phenomenological baseline to understand the Earth atmosphere in context. In this white paper, we outline the key science that such a facility would make possible and outline the instrumentation to be implemented.

We present an extension of the code ProDiMo that allows for a modeling of processes pertinent to active galactic nuclei and to an ambient chemistry that is time dependent. We present a proof-of-concept and focus on a few astrophysically relevant species, e.g., H+, H2+ and H3+; C+ and N+; C and O; CO and H2O; OH+, H2O+ and H3O+; HCN and HCO+. We find that the freeze-out of water is strongly suppressed and that this affects the bulk of the oxygen and carbon chemistry occurring in AGN. The commonly used AGN tracer HCN/HCO+ is strongly time-dependent, with ratios that vary over orders of magnitude for times longer than 10^4 years. Through ALMA observations this ratio can be used to probe how the narrow-line region evolves under large fluctuations in the SMBH accretion rate. Strong evolutionary trends, on time scales of 10^4-10^8 years, are also found in species such as H3O+, CO, and H2O. These reflect, respectively, time dependent effects in the ionization balance, the transient nature of the production of molecular gas, and the freeze-out/sublimation of water.

Context. The seeds of the supermassive black holes (SMBHs) with masses of ~10^9 M_Sun observed already at z ~ 6 may have formed through the direct collapse of primordial gas in T_vir >~ 10^4 K halos, whereby the gas must stay hot (~10^4 K) in order to avoid fragmentation. Aims. The interplay between magnetic fields, turbulence, and a UV radiation background during the gravitational collapse of primordial gas in a halo is explored; in particular, the possibilities for avoiding fragmentation are examined. Methods. Using an analytical one-zone model, the evolution of a cloud of primordial gas is followed from its initial cosmic expansion through turnaround, virialization, and collapse up to a density of 10^7 cm-3. Results. It was found that in halos with no significant turbulence, the critical UV background intensity (J_21^crit) for keeping the gas hot is lower by a factor ~10 for an initial comoving magnetic field B_0 ~ 2 nG than for the zero-field case, and even lower for stronger fields. In turbulent halos, J_21^crit is found to be a factor ~10 lower than for the zero-field-zero-turbulence case, and the stronger the turbulence (more massive halo and/or stronger turbulent heating) the lower J_21^crit. Conclusions. The reduction in J_21^crit is particularly important, since it exponentially increases the number of halos exposed to a supercritical radiation background.

It is argued that a change in particle rest mass must involve a multiply connected space-time topology. The LHC can probe these topological effects through the particle leakage (~18%) that it experiences in particle mass changing interactions. With a probability of ~82%, this 4-space leakage causes a downward shift in the observable Higgs mass, from its physical value of 131.6 GeV to 125.2 GeV.

Observational studies show that there is a strong link between the formation and evolution of galaxies and the growth of supermassive black holes (SMBH) at their centers. However, the underlying physics of this observed relation is poorly understood. In order to study the effects of X-ray radiation on the surroundings of the black hole, we implement X-ray Dominated Region (XDR) physics into Enzo and use the radiation transport module Moray to calculate the radiative transfer for a polychromatic spectrum. In this work, we investigate the effects of X-ray irradiation, produced by a central massive black hole (MBH) with a mass of M = 5x10^4 M_(solar), on ambient gas with solar and zero metallicity. We find that in the solar metallicity case, due to high opacity of the metals, the energy deposition rate in the central region (< 20 pc) is high and hence the temperatures in the center are on the order of 10^(5-7) K. Moreover, due to the cooling ability and high intrinsic opacity of solar metallicity gas, column densities of 10^(24) cm^(-2) are reached at a radius of 20 pc from the MBH. These column densities are about 3 orders of magnitudes higher than in the zero metallicity case. Furthermore, in the zero metallicity case an X-ray induced H II region is formed already after 5.8 Myr. This causes a significant outflow of gas (~8x10^6 M_(solar) from the central region, with the gas reaching outflow velocities up to ~100 km s^(-1). At later times, ~23 Myr after we insert the MBH, we find that the solar metallicity case also develops an X-ray induced H II region, but delayed by ~17 Myr.

Water is key in the evolution of protoplanetary disks and the formation of comets and icy/water planets. While high excitation water lines originating in the hot inner disk have been detected in several T Tauri stars (TTSs), water vapor from the outer disk, where most of water ice reservoir is stored, was only reported in the closeby TTS TW Hya. We present spectrally resolved Herschel/HIFI observations of the young TTS DG Tau in the ortho- and para- water ground-state transitions at 557, 1113 GHz. The lines show a narrow double-peaked profile, consistent with an origin in the outer disk, and are ~19-26 times brighter than in TW Hya. In contrast, CO and [C II] lines are dominated by emission from the envelope/outflow, which makes H2O lines a unique tracer of the disk of DG Tau. Disk modeling with the thermo-chemical code ProDiMo indicates that the strong UV field, due to the young age and strong accretion of DG Tau, irradiates a disk upper layer at 10-90 AU from the star, heating it up to temperatures of 600 K and producing the observed bright water lines. The models suggest a disk mass of 0.015-0.1 Msun, consistent with the estimated minimum mass of the solar nebula before planet formation, and a water reservoir of ~1e2-1e3 Earth oceans in vapour, and ~100 times larger in the form of ice. Hence, this detection supports the scenario of ocean delivery on terrestrial planets by impact of icy bodies forming in the outer disk.

A topological extension of general relativity is presented. The superposition principle of quantum mechanics, as formulated by the Feynman path integral, is taken as a starting point. It is argued that the trajectories that enter this path integral are distinct, despite any quantum uncertainty in geometry, and thus that space-time topology is multiply connected. Specifically, space-time at the Planck scale consists of a lattice of three-tori that facilitates many distinct paths for particles to travel along. To add gravity, mini black holes are attached to this lattice. These mini black holes represent Wheeler's quantum foam and result from the fact that GR is not conformally invariant. The stable creation of such mini black holes is found to be caused by the existence of macroscopic (so long-lived) black holes. This connection, by which macroscopic black holes induce mini black holes, is a topological expression of Mach's principle. The proposed topological extension of GR can be tested because, if correct, the dark energy density of the universe should be linearly proportional to the total number of macroscopic black holes in the universe at any time. This prediction, although strange, agrees with current astrophysical observations.

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

The first deep blank-field 450um map (1-sigma~1.3mJy) from the SCUBA-2 Cosmology Legacy Survey (S2CLS), conducted with the James Clerk Maxwell Telescope (JCMT) is presented. Our map covers 140 arcmin^2 of the COSMOS field, in the footprint of the HST CANDELS area. Using 60 submillimetre galaxies (SMGs) detected at >3.75-sigma, we evaluate the number counts of 450um-selected galaxies with flux densities S_450>5mJy. The 8-arcsec JCMT beam and high sensitivity of SCUBA-2 now make it possible to directly resolve a larger fraction of the cosmic infrared background (CIB, peaking at ~200um) into the individual galaxies responsible for its emission than has previously been possible at this wavelength. At S_450>5mJy we resolve (7.4[+/-]0.7)x10^-2 MJy/sr of the CIB at 450um (equivalent to 16[+/-]7% of the absolute brightness measured by COBE at this wavelength) into point sources. A further ~40% of the CIB can be recovered through a statistical stack of 24um emitters in this field, indicating that the majority (~60%) of the CIB at 450um is emitted by galaxies with S_450>2mJy. The average redshift of 450um emitters identified with an optical/near-infrared counterpart is estimated to be <z>=1.3, implying that the galaxies in the sample are in the ultraluminous class (L_IR~1.1x10^12 L_sun). If the galaxies contributing to the statistical stack lie at similar redshifts, then the majority of the CIB at 450um is emitted by galaxies in the LIRG class with L_IR>3.6x10^11 L_sun.

Molecular hydrogen is the most abundant molecule in the universe. A large fraction of H2 forms by association of hydrogen atoms adsorbed on polycyclic aromatic hydrocarbons (PAHs), where formation rates depend crucially on the H sticking probability. We have experimentally studied PAH hydrogenation by exposing coronene cations, confined in a radiofrequency ion trap, to gas phase atomic hydrogen. A systematic increase of the number of H atoms adsorbed on the coronene with the time of exposure is observed. Odd coronene hydrogenation states dominate the mass spectrum up to 11 H atoms attached. This indicates the presence of a barrier preventing H attachment to these molecular systems. For the second and fourth hydrogenation, barrier heights of 72 +- 6 meV and 40 +- 10 meV, respectively are found which is in good agreement with theoretical predictions for the hydrogenation of neutral PAHs. Our experiments however prove that the barrier does not vanish for higher hydrogenation states. These results imply that PAH cations, as their neutral counterparts, exist in highly hydrogenated forms in the interstellar medium. Due to this catalytic activity, PAH cations and neutrals seem to contribute similarly to the formation of H2.

Most of the mass in protoplanetary disks is in the form of gas. The study of the gas and its diagnostics is of fundamental importance in order to achieve a detailed description of the thermal and chemical structure of the disk. The radiation from the central star (from optical to X-ray wavelengths) and viscous accretion are the main source of energy and dominates the disk physics and chemistry in its early stages. This is the environment in which the first phases of planet formation will proceed. We investigate how stellar and disk parameters impact the fine-structure cooling lines [NeII], [ArII], [OI], [CII] and H2O rotational lines in the disk. These lines are potentially powerful diagnostics of the disk structure and their modelling permits a thorough interpretation of the observations carried out with instrumental facilities such as Spitzer and Herschel. Following Aresu et al. (2011), we computed a grid of 240 disk models, in which the X-ray luminosity, UV-excess luminosity, minimum dust grain size, dust size distribution power law and surface density distribution power law, are systematically varied. We solve self-consistently for the disk vertical hydrostatic structure in every model and apply detailed line radiative transfer to calculate line fluxes and profiles for a series of well known mid- and far-infrared cooling lines. The [OI] 63 micron line flux increases with increasing FUV luminosity when Lx < 1e30 erg/s, and with increasing X-ray luminosity when LX > 1e30 erg/s. [CII] 157 micron is mainly driven by FUV luminosity via C+ production, X-rays affect the line flux to a lesser extent. [NeII] 12.8 micron correlates with X-rays; the line profile emitted from the disk atmosphere shows a double-peaked component, caused by emission in the static disk atmosphere, next to a high velocity double-peaked component, caused by emission in the very inner rim. (abridged)

Context. Planets are thought to eventually form from the mostly gaseous (~99% of the mass) disks around young stars. The density structure and chemical composition of protoplanetary disks are affected by the incident radiation field at optical, FUV, and X-ray wavelengths, as well as by the dust properties. Aims. The effect of FUV and X-rays on the disk structure and the gas chemical composition are investigated. This work forms the basis of a second paper, which discusses the impact on diagnostic lines of, e.g., C+, O, H2O, and Ne+ observed with facilities such as Spitzer and Herschel. Methods. A grid of 240 models is computed in which the X-ray and FUV luminosity, minimum grain size, dust size distribution, and surface density distribution are varied in a systematic way. The hydrostatic structure and the thermo-chemical structure are calculated using ProDiMo. Results. The abundance structure of neutral oxygen is stable to changes in the X-ray and FUV luminosity, and the emission lines will thus be useful tracers of the disk mass and temperature. The C+ abundance distribution is sensitive to both X-rays and FUV. The radial column density profile shows two peaks, one at the inner rim and a second one at a radius r=5-10 AU. Ne+ and other heavy elements have a very strong response to X-rays, and the column density in the inner disk increases by two orders of magnitude from the lowest (LX = 1e29 erg/s) to the highest considered X-ray flux (LX = 1e32 erg/s). FUV confines the Ne+ ionized region to areas closer to the star at low X-ray luminosities (LX = 1e29 erg/s). H2O abundances are enhanced by X-rays due to higher temperatures in the inner disk and higher ionization fractions in the outer disk. The line fluxes and profiles are affected by the effects on these species, thus providing diagnostic value in the study of FUV and X-ray irradiated disks around T Tauri stars. (abridged)

We present Herschel-PACS spectroscopy of the [OI]63um far-infrared cooling line from a sample of six unlensed and spectroscopically-confirmed 870um-selected submillimetre (submm) galaxies (SMGs) at 1.1<z<1.6 from the LABOCA Extended Chandra Deep Field South (ECDFS) Submm Survey (LESS). This is the first survey of [OI]63um, one of the main photodissociation region (PDR) cooling lines, in SMGs. New high-resolution ALMA interferometric 870um continuum imaging confirms that these six Herschel-targeted SMG counterparts are bona fide sources of submm emission. We detect [OI]63um in two SMGs with a SNR >3, tentatively detect [OI]63um in one SMG, and constrain the line flux for the non-detections. We also exploit the combination of submm continuum photometry from 250-870um and our new PACS continuum measurements to constrain the far-infrared (FIR) luminosity, L_FIR, in these SMGs to < 30%. We find that SMGs do not show a deficit in their [OI]63um-to-far-infrared continuum luminosity ratios (with ratios ranging from ~0.5-1.5%), similar to what was seen previously for the [CII]158um-to-FIR ratios in SMGs. These observed ratios are about an order of magnitude higher than what is seen typically for local ultra luminous infrared galaxies (ULIRGs), which adds to the growing body of evidence that SMGs are not simply `scaled up' versions of local ULIRGs. Rather, the PDR line-to-L_FIR ratios suggest that the star formation modes of SMGs are likely more akin to that of local normal (lower-luminosity) star-forming galaxies, with the bulk of the star formation occurring in extended regions, galaxy-scale (~kpc) in size. These observations represent the first step towards a census of the major PDR cooling lines in typical SMGs that will be attainable with ALMA, enabling detailed modelling to probe the global properties of the star formation and the evolutionary status of SMGs.

Star formation in the centers of galaxies is thought to yield massive stars with a possibly top-heavy stellar mass distribution. It is likely that magnetic fields play a crucial role in the distribution of stellar masses inside star-forming molecular clouds. In this context, we explore the effects of magnetic fields, with a typical field strength of 38 microG, such as in RCW 38, and a field strength of 135 microG, similar to NGC 2024 and the infrared dark cloud G28.34+0.06, on the initial mass function (IMF) near (< 10 pc) a 10^7 solar mass black hole. Using these conditions, we perform a series of numerical simulations with the hydrodynamical code FLASH to elucidate the impact of magnetic fields on the IMF and the star-formation efficiency (SFE) emerging from an 800 solar mass cloud. We find that the collapse of a gravitationally unstable molecular cloud is slowed down with increasing magnetic field strength and that stars form along the field lines. The total number of stars formed during the simulations increases by a factor of 1.5-2 with magnetic fields. The main component of the IMF has a lognormal shape, with its peak shifted to sub-solar (< 0.3 M_sun) masses in the presence of magnetic fields, due to a decrease in the accretion rates from the gas reservoir. In addition, we see a top-heavy, nearly flat IMF above ~2 solar masses, from regions that were supported by magnetic pressure until high masses are reached. We also consider the effects of X-ray irradiation if the central black hole is active. X-ray feedback inhibits the formation of sub-solar masses and decreases the SFEs even further. Thus, the second contribution is no longer visible. We conclude that magnetic fields potentially change the SFE and the IMF both in active and inactive galaxies, and need to be taken into account in such calculations.

In 1997, an extension of general relativity was proposed that predicts the dark energy density \Lambda\ to vary linearly with the total number of macroscopic black holes in the universe. We explore this prediction and find that \Lambda\ must be roughly constant after the bulk of the stellar mass black holes are in place, so for a redshift z smaller than unity. Conversely, the highest black hole formation rate corresponds to the peak in the cosmic star formation history, earlier than z=1. This yields a fast declining \Lambda, by a factor of about 5, from redshift 1 to 3. At even earlier times, before many stars were formed, the value of \Lambda\ should be much smaller than its current value. These predicted effects are all consistent with current data, and near future observations can definitively confirm or disproof the link between the dark energy density and the total number of black holes in the universe.

The star formation rates in starburst galaxies are orders of magnitude higher than in local star-forming regions, and the origin of this difference is not well understood. We use sub-mm spectral line maps to characterize the physical conditions of the molecular gas in the luminous Galactic star-forming region W49A and compare them with the conditions in starburst galaxies. We probe the temperature and density structure of W49A using H_2CO and HCN line ratios over a 2'x2' (6.6x6.6 pc) field with an angular resolution of 15" (~0.8 pc) provided by the JCMT Spectral Legacy Survey. We analyze the rotation diagrams of lines with multiple transitions with corrections for optical depth and beam dilution, and estimate excitation temperatures and column densities. Comparing the observed line intensity ratios with non-LTE radiative transfer models, our results reveal an extended region (about 1'x1', equivalent to ~3x3 pc at the distance of W49A) of warm (> 100 K) and dense (>10^5 cm^-3) molecular gas, with a mass of 2x10^4 - 2x10^5 M_Sun (by applying abundances derived for other regions of massive star-formation). These temperatures and densities in W49A are comparable to those found in clouds near the center of the Milky Way and in starburst galaxies. The highly excited gas is likely to be heated via shocks from the stellar winds of embedded, O-type stars or alternatively due to UV irradiation, or possibly a combination of these two processes. Cosmic rays, X-ray irradiation and gas-grain collisional heating are less likely to be the source of the heating in the case of W49A.

Numerical simulations suggest that the first galaxies are formed in protogalactic halos with virial temperatures >= 10^4 K. It is likely that such halos are polluted with trace amounts of metals produced by the first generation of stars. The presence of dust can significantly change the chemistry and dynamics of early galaxies. In this article, we aim to assess the role of dust on the thermal and dynamical evolution of the first galaxies in the presence of a background UV flux, and its implications for the observability of Lyman alpha emitters and sub-mm sources. We have performed high resolution cosmological simulations using the adaptive mesh refinement code FLASH to accomplish this goal. We have developed a chemical network appropriate for these conditions and coupled it with the FLASH code. The main ingredients of our chemical model include the formation of molecules, a multi-level treatment of atomic hydrogen, line trapping of Lyman alpha photons and, photoionization and photodissociation processes in a UV background. We found that the formation of molecules (H_{2} and HD) is significantly enhanced in the presence of dust grains. The presence of a background UV flux strongly influences the formation of molecules by photodissociating them. We explore the evolution after a major merger, leading to the formation of a binary disk. These disks have gas masses of ~10^{7} M_sun at a redshift of 5.4. Each disk lies in a separate subhalo as a result of the merger event. The disks are supported by turbulent pressure due to the highly supersonic turbulence present in the halo. The presence of dust does not significantly reduce the Lyman alpha emission. The emission of Lyman alpha is extended and originates from the envelope of the halo due to line trapping effects. We also find that dust masses of 10^8 M_sun are required to observe the dust continuum emission from z 5 galaxies with ALMA.

Context. Multitransition CO observations of galaxy centers have revealed that significant fractions of the dense circumnuclear gas have high kinetic temperatures, which are hard to explain by pure photon excitation, but may be caused by dissipation of turbulent energy. Aims. We aim to determine to what extent mechanical heating should be taken into account while modelling PDRs. To this end, the effect of dissipated turbulence on the thermal and chemical properties of PDRs is explored. Methods. Clouds are modelled as 1D semi-infinite slabs whose thermal and chemical equilibrium is solved for using the Leiden PDR-XDR code. Results. In a steady-state treatment, mechanical heating seems to play an important role in determining the kinetic temperature of the gas in molecular clouds. Particularly in high-energy environments such as starburst galaxies and galaxy centers, model gas temperatures are underestimated by at least a factor of two if mechanical heating is ignored. The models also show that CO, HCN and H2 O column densities increase as a function of mechanical heating. The HNC/HCN integrated column density ratio shows a decrease by a factor of at least two in high density regions with n \sim 105 cm-3, whereas that of HCN/HCO+ shows a strong dependence on mechanical heating for this same density range, with boosts of up to three orders of magnitude. Conclusions. The effects of mechanical heating cannot be ignored in studies of the molecular gas excitation whenever the ratio of the star formation rate to the gas density is close to, or exceeds, 7 \times 10-6 M yr-1 cm4.5 . If mechanical heating is not included, predicted column densities are underestimated, sometimes even by a few orders of magnitude. As a lower bound to its importance, we determined that it has non-negligible effects already when mechanical heating is as little as 1% of the UV heating in a PDR.

We present detailed model fits to observations of the disc around the Herbig Ae star HD 163296. This well-studied object has an age of ~ 4 Myr, with evidence of a circumstellar disc extending out to ~ 540AU. We use the radiation thermo-chemical disc code ProDiMo to model the gas and dust in the circumstellar disc of HD 163296, and attempt to determine the disc properties by fitting to observational line and continuum data. These include new Herschel/PACS observations obtained as part of the open-time key program GASPS (Gas in Protoplanetary Systems), consisting of a detection of the [OI]63mic line and upper limits for several other far infrared lines. We complement this with continuum data and ground-based observations of the 12CO 3-2, 2-1 and 13CO J=1-0 line transitions, as well as the H2 S(1) transition. We explore the effects of stellar ultraviolet variability and dust settling on the line emission, and on the derived disc properties. Our fitting efforts lead to derived gas/dust ratios in the range 9-100, depending on the assumptions made. We note that the line fluxes are sensitive in general to the degree of dust settling in the disc, with an increase in line flux for settled models. This is most pronounced in lines which are formed in the warm gas in the inner disc, but the low excitation molecular lines are also affected. This has serious implications for attempts to derive the disc gas mass from line observations. We derive fractional PAH abundances between 0.007 and 0.04 relative to ISM levels. Using a stellar and UV excess input spectrum based on a detailed analysis of observations, we find that the all observations are consistent with the previously assumed disc geometry.

The ULIRG Mrk 231 exhibits very strong water rotational lines between \lambda = 200-670\mu m, comparable to the strength of the CO rotational lines. High redshift quasars also show similar CO and H2O line properties, while starburst galaxies, such as M82, lack these very strong H2O lines in the same wavelength range, but do show strong CO lines. We explore the possibility of enhancing the gas phase H2O abundance in X-ray exposed environments, using bare interstellar carbonaceous dust grains as a catalyst. Cloud-cloud collisions cause C and J shocks, and strip the grains of their ice layers. The internal UV field created by X-rays from the accreting black hole does not allow to reform the ice. We determine formation rates of both OH and H2O on dust grains, having temperature T_dust=10-60 K, using both Monte Carlo as well as rate equation method simulations. The acquired formation rates are added to our X-ray chemistry code, that allows us to calculate the thermal and chemical structure of the interstellar medium near an active galactic nucleus. We derive analytic expressions for the formation of OH and H2O on bare dust grains as a catalyst. Oxygen atoms arriving on the dust are released into the gas phase under the form of OH and H2O. The efficiencies of this conversion due to the chemistry occurring on dust are of order 30 percent for oxygen converted into OH and 60 percent for oxygen converted into H_2O between T_dust=15-40 K. At higher temperatures, the efficiencies rapidly decline. When the gas is mostly atomic, molecule formation on dust is dominant over the gas-phase route, which is then quenched by the low H2 abundance. Here, it is possible to enhance the warm (T> 200 K) water abundance by an order of magnitude in X-ray exposed environments. This helps to explain the observed bright water lines in nearby and high-redshift ULIRGs and Quasars.

The conditions that affect the formation of stars in radiatively and mechanically active environments are quite different than the conditions that apply to our local interstellar neighborhood. In such galactic environments, a variety of feedback processes can play a significant role in shaping the initial mass function (IMF). Here, we present a numerical study on the effects of an accreting black hole and the influence of nearby massive stars to a collapsing, 800 M_sun, molecular cloud at 10 pc distance from the black hole. We parametrize and study radiative feedback effects of hard X-rays emanating from the black hole broad line region, increased cosmic ray rates due to supernovae in starbursts, and strong UV radiation produced by nearby massive stars. We also investigate the importance of shear from the supermassive, 10^6-10^8 M_sun, black hole as the star-forming cloud orbits around it. We find that thermal pressure from X-rays compresses the cloud, which induces a high star formation rate early on, but reduces the overall star formation efficiency to about 7% due to gas depletion by evaporation. We see that the turn-over mass of the IMF increases up to a factor of 2.3, M_turn = 1-1.5 M_sun, for the model with the highest X-ray flux (160 erg s^-1 cm^-2), while the high-mass slope of the IMF becomes Gamma > -1. This results in more high mass stars and a non-Salpeter IMF. Cosmic rays penetrate deeply into the cloud and increase the gas temperature (50-200 K), which leads to a reduced formation efficiency of low mass stars. High cosmic ray rates increase the average mass of stars, thereby shifting the turn-over mass to higher values, i.e., up to several solar masses. Due to this process, the onset of star formation is also delayed. We conclude that the IMF inside active galaxies is different than the one obtained from local environments.

Observations of star forming environments revealed that the abundances of some deuterated interstellar molecules are markedly larger than the cosmic D/H ratio of 10-5. Possible reasons for this pointed to grain surface chemistry. How- ever, organic molecules and water, which are both ice constituents, do not enjoy the same deuteration. For example, deuterated formaldehyde is very abundant in comets and star forming regions, while deuterated water rarely is. In this article, we explain this selective deuteration by following the formation of ices (using the rate equation method) in translucent clouds, as well as their evolu- tion as the cloud collapses to form a star. Ices start with the deposition of gas phase CO and O onto dust grains. While reaction of oxygen with atoms (H or D) or molecules (H2) yields H2O (HDO), CO only reacts with atoms (H and D) to form H2CO (HDCO, D2CO). As a result, the deuteration of formaldehyde is sensitive to the gas D/H ratio as the cloud undergoes gravitational collapse, while the deuteration of water strongly depends on the dust temperature at the time of ice formation. These results reproduce well the deuterium fractionation of formaldehyde observed in comets and star forming regions and can explain the wide spread of deuterium fractionation of water observed in these environments.

We present the detection of four rotational emission lines of water vapor, from energy levels Eu/k= 101 - 454 K, in the gravitationally lensed z=3.9 QSO host galaxy APM08279+5255. While the lowest H2O lines are collisionally excited in clumps of warm, dense gas (density of hydrogen nuclei n_H=(3.1 +/- 1.2) x 10^6 cm^-3, gas temperature T_g ~ 105 +/- 21 K), we find that the excitation of the higher lines is dominated by the intense local infrared radiation field. Since only collisionally excited emission contributes to gas cooling, we conclude that H2O is not a significant coolant of the warm molecular gas. Our excitation model requires the radiatively excited gas to be located in an extended region of high 100 micron opacity (tau_100 = 0.9 +/- 0.2). Locally, such extended infrared-opaque regions are found only in the nuclei of ultraluminous infrared galaxies. We propose a model where the infrared-opaque circumnuclear cloud, which is penetrated by the X-ray radiation field of the QSO nucleus, contains clumps of massive star formation where the H2O emission originates. The radiation pressure from the intense local infrared radiation field exceeds the thermal gas pressure by about an order of magnitude, suggesting close to Eddington-limited star formation in these clumps.

We study the effects of local environmental conditions affecting the diffuse interstellar band (DIB) carriers within the Upper Scorpius subgroup of the Sco OB2 association. We obtained optical spectra with FEROS, and measured the equivalent widths of five DIBs (at 5780, 5797, 6196, 6379, and 6613 Angstrom) as well as those of absorption lines of di-atomic molecules (CH, CH+, CN) and atoms (KI, CaI) towards 89 targets in the direction of Upper Scorpius. We construct a simple radiative transfer and chemical network model of the diffuse interstellar medium (ISM) sheet in front of Upp Sco to infer the effective radiation field. Both the interstellar radiation field strength, I_UV, and the molecular hydrogen fraction, f_H2, have been derived for 55 sightlines probing the Upp Sco ISM. We discuss the relations between DIB strengths, CH and CH+ line strengths, E(B-V), I_UV, and f_H2. The ratio between the 5780 and 5797 DIBs reveals a (spatial) dependence on the local environment in terms of cloud density and exposure to the interstellar radiation field, reflecting the molecular nature of these DIB carriers.

The initial mass function (IMF) of the first (Population III) stars and Population II (Pop II) stars is poorly known due to a lack of observations of the period between recombination and reionization. In simulations of the formation of the first stars, it has been shown that, due to the limited ability of metal-free primordial gas to cool, the IMF of the first stars is a few orders of magnitude more massive than the current IMF. The transition from a high-mass IMF of the first stars to a lower-mass current IMF is thus important to understand. To study the underlying physics of this transition, we performed several simulations using the cosmological hydrodynamic adaptive mesh refinement code Enzo for metallicities of 10^{-4}, 10^{-3}, 10^{-2}, and 10^{-1} Z_{\odot}. In our simulations we include a star formation prescription that is derived from a metallicity dependent multi-phase ISM structure, an external UV radiation field, and a mechanical feedback algorithm. We also implement cosmic ray heating, photoelectric heating and gas-dust heating/cooling, and follow the metal enrichment of the ISM. It is found that the interplay between metallicity and UV radiation leads to the co-existence of Pop III and Pop II star formation in non-zero metallicity (Z/Z_{\odot} \geq10^{-2}) gas. A cold (T<100 K) and dense (\rho>10^{-22} g cm^{-3}) gas phase is fragile to ambient UV radiation. In a metal-poor (Z/Z_{\odot} \leq10^{-3}) gas, the cold and dense gas phase does not form in the presence of a radiation field of F_{0}\sim10^{-5}-10^{-4} erg cm^{-2} s^{-1}. Therefore, metallicity by itself is not a good indicator of the Pop III-Pop II transition. Metal-rich (Z/Z_{\odot}\geq10^{-2}) gas dynamically evolves two to three orders of magnitude faster than metal poor gas (Z/Z_{\odot}\leq10^{-3}). The simulations including SNe show that pre-enrichment of the halo does not affect the mixing of metals.

Aims. Interstellar dust grains, because of their catalytic properties, are crucial to the formation of H2, the most abundant molecule in the Universe. The formation of molecular hydrogen strongly depends on the ability of H atoms to stick on dust grains. In this study we determine the sticking coefficient of H atoms chemisorbed on graphitic surfaces, and estimate its impact on the forma- tion of H2. Methods. The sticking probability of H atoms chemisorbed onto graphitic surfaces is obtained using a mixed classical-quantum dynamics method. In this, the H atom is treated quantum- mechanically and the vibrational modes of the surface are treated classically. The implications of sticking for the formation of H2 are addressed by using Kinetic Monte Carlo simulations that follow how atoms stick, move and associate with each other on dust surfaces of different temper- ature. Results. In our model, molecular hydrogen forms very efficiently for dust temperatures lower than 15 K through the involvement of physisorbed H atoms. At dust temperatures higher than 15 K and gas temperatures lower than 2000 K, H2 formation differs strongly if the H atoms com- ing from the gas phase have to cross a square barrier (usually considered in previous studies) or a barrier obtained by DFT calculations to become chemisorbed. The product of sticking times efficiency can be increased by many orders of magnitude when realistic barriers are considered. If graphite phonons are taken into account in the dynamics calculations, then H atoms stick better on the surface at high energies, but the overall H2 formation efficiency is only slightly affected. Our results suggest that H2 formation can proceed efficiently in photon dominated regions, X-ray dominated regions, hot cores and in the early Universe when the first dust is available.

The nearly edge-on galaxy NGC4945 is one of the closest galaxies where an AGN and starburst coexist, and is one of the brightest sources at 100 keV. Near and mid-infrared spectroscopy have shown very strong obscuration of its central region, rivaled only in strength by some of the most deeply obscured ULIRGs. We aim to determine the spatial distribution of ISM features in the central 426x426 pc^2 of NGC4945. We map the central region of NGC4945 in three Spitzer-IRS modules (SH, SL and LL). We produce maps of the flux distribution of the starburst tracers [Ne II], [Ne III], [S III] and [S IV] at 12.81, 15.56, 18.71 and 10.51 mum, respectively, and a map of the AGN narrow-line region tracer [Ne V] at 14.32 mum. We also mapped the S(1), S(2) and S(3) pure rotational lines of H2, which trace the distribution of warm molecular hydrogen. We obtained an extinction map (A_V) based on the apparent strength of the 9.7 mum silicate absorption feature. Our A_V map traces the contours of the starburst ring but the highest extinction (A_V(9.85 mum)~60) is found at the nucleus. Within the uncertainty of the astrometry all emission lines are found to peak on the nucleus, except for the warm molecular hydrogen emission which shows a maximum 60-100 pc NW of the nucleus. We favour a scenario in which the lower H2 0-0 S(1) and S(2) rotational lines originate mainly from an unobscured extra-nuclear component associated with the super-wind cone observed in the HST NICMOS map of the H2 1-0 S(1) vibrational line. For the [Ne V] emission we infer an attenuation of a factor 12-160 (A_V=55-112) based on a comparison of the ratio of our [Ne V] flux and the absorption-corrected 14-195 keV Swift-BAT flux to the average [Ne V]/BAT ratio for Seyfert 1 nuclei. The high attenuation indicates that [Ne V] and [O IV] cannot be used as extinction-free tracers of AGN power in galaxies with deeply buried nuclei.

The Lyman alpha line is a robust tracer of high redshift galaxies. We present estimates of Lyman alpha emission from a protogalactic halo illuminated by UV background radiation fields with various intensities. For this purpose, we performed cosmological hydrodynamics simulations with the adaptive mesh refinement code FLASH, including a detailed network for primordial chemistry,comprising the formation of primordial molecules, a multi-level model for the hydrogen atom as well as the photo-ionization and photo-dissociation processes in a UV background. We find that the presence of a background radiation field J_21 excites the emission of Lyman alpha photons, increasing the Lyman alpha luminosity up to two orders of magnitude. For a halo of \sim 10^10 M_sun, we find that a maximum flux of 5 \times 10^-15 erg cm^-2 s^-1 is obtained for J21 \times f_esc = 0.1, where f_esc is the escape fraction of the ionizing radiation. Depending on the environmental conditions, the flux may vary by three orders of magnitude. For J_21 \times f_esc > 0.1 the Lyman alpha luminosity decreases as the atomic hydrogen abundance becomes rather small. The fluxes derived here can be probed using Subaru and the upcoming James Webb Space Telescope. The emission of Lyman alpha photons is extended and comes from the envelope of the halo rather than its core. In the center of the halo, line trapping becomes effective above columns of 10^22 cm^-2 and suppresses the emission of Lyman alpha. In addition, cooling by primordial molecules may decrease the gas temperature in the central region, which further reduces Lyman alpha emission. In the central core, H_2 is photo-dissociated for a background flux of J_21 \geq 1000. For weaker radiation fields, i.e. J_21 < 0.1, H_2 and HD cooling are particularly strong in the center of the halo, leading to gas temperatures as low as \sim 100 K.

Recent FERMI observations provide a lower limit of 10^{-15} G for the magnetic field strength in the intergalactic medium (IGM). This is consistent with theoretical expectations based on the Biermann battery effect, which predicts such IGM fields already at redshifts z~10. During gravitational collapse, such magnetic fields can be amplified by compression and by turbulence, giving rise to the small-scale dynamo. On scales below the Jeans length, the eddy turnover timescale is much shorter than the free-fall timescale, so that saturation can be reached during collapse. This scenario has been tested and confirmed with magneto-hydrodynamical simulations following the collapse of a turbulent, weakly magnetized cloud. Based on a spectral analysis, we confirm that turbulence is injected on the Jeans scale. For the power spectrum of the magnetic field, we obtain the Kazantsev slope which is characteristic for the small-scale dynamo. A calculation of the critical length scales for ambipolar diffusion and Ohmic dissipation shows that these scales are always small enough to allow significant amplification of the magnetic field by small-scale eddies. We discuss potential implications for the protostellar accretion disk, with particular focus on the magneto-rotational instability, which may change the morphology of the disk and reduce the accretion rate by a factor of a few.

Many efforts have been made to model the mass distribution and dynamical evolution of the circumnuclear gas in active galactic nuclei (AGNs). However, chemical evolution is not included in detail in three-dimensional (3-D) hydrodynamic simulations. The X-ray radiation from the AGN can drive the gas chemistry and affect the thermodynamics, as well as the excitation of the interstellar medium (ISM). Therefore, we estimate the effects (on chemical abundances and excitation) of X-ray irradiation by the AGN, for atomic and molecular gas in a 3-D hydrodynamic model of an AGN torus. We obtain the abundances of various species from an X-ray chemical model. A 3-D radiative transfer code estimates the level populations, which result in line intensity maps. Predictions for the CO J=1-0 to J=9-8 lines indicate that mid-J CO lines are excellent probes of density and dynamics in the central (<60 pc) region of the AGN, in contrast to the low-J CO lines. Analysis of the X_CO/\alpha conversion factors shows that only the higher-J CO lines can be used for gas mass determination in AGN tori. The [C II] 158 um emission traces mostly the hot (T_k>1000m K) central (<60 pc) region of the AGN torus. The [C II] 158 um line will be useful for ALMA observations of high redshift (z>1) AGNs. The spatial scales (>0.25 pc) probed with our simulations match the size of the structures that ALMA will resolve in nearby (<45 Mpc at 0.01") galaxies.

The formation of the first galaxies is accompanied by large accretion flows and virialization shocks, during which the gas is shock-heated to temperatures of $\sim10^4$ K, leading to potentially strong fluxes in the Lyman alpha line. Indeed, a number of Lyman alpha blobs has been detected at high redshift. In this letter, we explore the origin of such Lyman alpha emission using cosmological hydrodynamical simulations that include a detailed model of atomic hydrogen as a multi-level atom and the effects of line trapping with the adaptive mesh refinement code FLASH. We see that baryons fall into the center of a halo through cold streams of gas, giving rise to a Lyman alpha luminosity of at least $\rm 10^{44} erg s^{-1}$ at $\rm z=4.7$, similar to observed Lyman alpha blobs. We find that a Lyman alpha flux of $\rm 5.0\times 10^{-17} erg cm^{-2} s^{-1}$ emerges from the envelope of the halo rather than its center, where the photons are efficiently trapped. Such emission can be probed in detail with the upcoming James Webb Space Telescope (JWST) and will constitute an important probe of gas infall and accretion.

The understanding of the formation process of massive stars (>8 Msun) is limited, due to theoretical complications and observational challenges. We investigate the physical structure of the large-scale (~10^4-10^5 AU) molecular envelope of the high-mass protostar AFGL2591 using spectral imaging in the 330-373 GHz regime from the JCMT Spectral Legacy Survey. Out of ~160 spectral features, this paper uses the 35 that are spatially resolved. The observed spatial distributions of a selection of six species are compared with radiative transfer models based on a static spherically symmetric structure, a dynamic spherical structure, and a static flattened structure. The maps of CO and its isotopic variations exhibit elongated geometries on scales of ~100", and smaller scale substructure is found in maps of N2H+, o-H2CO, CS, SO2, CCH, and methanol lines. A velocity gradient is apparent in maps of all molecular lines presented here, except SO, SO2, and H2CO. We find two emission peaks in warm (Eup~200K) methanol separated by 12", indicative of a secondary heating source in the envelope. The spherical models are able to explain the distribution of emission for the optically thin H13CO+ and C34S, but not for the optically thick HCN, HCO+, and CS, nor for the optically thin C17O. The introduction of velocity structure mitigates the optical depth effects, but does not fully explain the observations, especially in the spectral dimension. A static flattened envelope viewed at a small inclination angle does slightly better. We conclude that a geometry of the envelope other than an isotropic static sphere is needed to circumvent line optical depth effects. We propose that this could be achieved in envelope models with an outflow cavity and/or inhomogeneous structure at scales smaller than ~10^4 AU. The picture of inhomogeneity is supported by observed substructure in at least six species.

Context: T Tauri stars have X-ray luminosities ranging from L_X = 10^28-10^32 erg/s. These luminosities are similar to UV luminosities (L_UV 10^30-10^31 erg/s) and therefore X-rays are expected to affect the physics and chemistry of the upper layers of their surrounding protoplanetary disks. Aim: The effects and importance of X-rays on the chemical and hydrostatic structure of protoplanetary disks are investigated, species tracing X-ray irradiation (for L_X >= 10^29 erg/s) are identified and predictions for [OI], [CII] and [NII] fine structure line fluxes are provided. Methods: We have implemented X-ray physics and chemistry into the chemo-physical disk code ProDiMo. We include Coulomb heating and H2 ionization as heating processes and primary and secondary ionization due to X-rays in the chemistry. Results: X-rays heat up the gas causing it to expand in the optically thin surface layers. Neutral molecular species are not much affected in their abundance and spatial distribution, but charged species such as N+, OH+, H2O+ and H3O+ show enhanced abundances in the disk surface. Conclusions: Coulomb heating by X-rays changes the vertical structure of the disk, yielding temperatures of ~ 8000 K out to distances of 50 AU. The chemical structure is altered by the high electron abundance in the gas in the disk surface, causing an efficient ion-molecule chemistry. The products of this, OH+, H2O+ and H3O+, are of great interest for observations of low-mass young stellar objects with the Herschel Space Observatory. [OI] (at 63 and 145 mic) and [CII] (at 158 mic) fine structure emission are only affected for L_X > 10^30 erg/s.

Context: Molecular data of extreme environments, such as Arp 220, but also NGC 253, show evidence for extremely high cosmic ray (CR) rates (10^3-10^4 * Milky Way) and mechanical heating from supernova driven turbulence. Aims: The consequences of high CR rates and mechanical heating on the chemistry in clouds are explored. Methods: PDR model predictions are made for low, n=10^3, and high, n=10^5.5 cm^-3, density clouds using well-tested chemistry and radiation transfer codes. Column densities of relevant species are discussed, and special attention is given to water related species. Fluxes are shown for fine-structure lines of O, C+, C, and N+, and molecular lines of CO, HCN, HNC, and HCO+. A comparison is made to an X-ray dominated region model. Results: Fine-structure lines of [CII], [CI], and [OI] are remarkably similar for different mechanical heating and CR rates, when already exposed to large amounts of UV. HCN and H2O abundances are boosted for very high mechanical heating rates, while ionized species are relatively unaffected. OH+ and H2O+ are enhanced for very high CR rates zeta > 5 * 10^-14 s^-1. A combination of OH+, OH, H2O+, H2O, and H3O+ trace the CR rates, and are able to distinguish between enhanced cosmic rays and X-rays.

Numerous cosmological simulations have been performed to study the formation of the first objects. We present the results of high resolution 3-D cosmological simulations of primordial objects formation using the adaptive mesh refinement code FLASH by including in an approximate manner the radiative transfer effects of Lyman alpha photons. We compare the results of a Lyman alpha trapping case inside gas clouds with atomic and molecular hydrogen cooling cases.The principal objective of this research is to follow the collapse of a zero metallicity halo with an effective equation of state (that accounts for the trapping) and to explore the fate of a halo in each of the three cases, specifically, the impact of thermodynamics on fragmentation of halos.Our results show that in the case of Lyman alpha trapping, fragmentation is halted and a massive object is formed at the center of a halo. The temperature of the gas remains well above $10^{4}$ K and the halo is not able to fragment to stellar masses. In the atomic cooling case, gas collapses into one or two massive clumps in contrast to the Lyman alpha trapping case. For the molecular hydrogen cooling case, gas cools efficiently and fragments.The formation of massive primordial objects is thus strongly dependent on the thermodynamics of the gas. A salient feature of our results is that for the formation of massive objects, e.g. intermediate mass black holes, feedback effects are not required to suppress $H_{2}$ cooling, as molecular hydrogen is collisionally dissociated at temperatures higher than $10^{4}$ K as a consequence of Lyman alpha trapping.

Star formation is regulated through a variety of feedback processes. In this study, we treat feedback by X-rays and discuss its implications. Our aim is to investigate whether star formation is significantly affected when a star forming cloud resides in the vicinity of a strong X-ray source. We perform an Eulerian grid simulation with embedded Lagrangian sink particles of a collapsing molecular cloud near a massive, 10^7 M_o black hole. The chemical and thermal changes caused by radiation are incorporated into the FLASH code. When there is strong X-ray feedback the star forming cloud fragments into larger clumps whereby fewer but more massive protostellar cores are formed. Competitive accretion has a strong impact on the mass function and a near-flat, non-Salpeter IMF results.

We explore the amplification of magnetic seed fields during the formation of the first stars and galaxies. During gravitational collapse, turbulence is created from accretion shocks, which may act to amplify weak magnetic fields in the protostellar cloud. Numerical simulations showed that such turbulence is sub-sonic in the first star-forming minihalos, and highly supersonic in the first galaxies with virial temperatures larger than 10^4 K. We investigate the magnetic field amplification during the collapse both for Kolmogorov and Burgers-type turbulence with a semi-analytic model that incorporates the effects of gravitational compression and small-scale dynamo amplification. We find that the magnetic field may be substantially amplified before the formation of a disk. On scales of 1/10 of the Jeans length, saturation occurs after ~10^8 yr. Although the saturation behaviour of the small-scale dynamo is still somewhat uncertain, we expect a saturation field strength of the order ~10^{-7} n^{0.5} G in the first star-forming halos, with n the number density in cgs units. In the first galaxies with higher turbulent velocities, the magnetic field strength may be increased by an order of magnitude, and saturation may occur after 10^6 to 10^7 yr. In the Kolmogorov case, the magnetic field strength on the integral scale (i.e. the scale with most magnetic power) is higher due to the characteristic power-law indices, but the difference is less than a factor of 2 in the saturated phase. Our results thus indicate that the precise scaling of the turbulent velocity with length scale is of minor importance. They further imply that magnetic fields will be significantly enhanced before the formation of a protostellar disk, where they may change the fragmentation properties of the gas and the accretion rate.

Context. Water together with O2 are important gas phase ingredients to cool dense gas in order to form stars. On dust grains, H2 O is an important constituent of the icy mantle in which a complex chemistry is taking place, as revealed by hot core observations. The formation of water can occur on dust grain surfaces, and can impact gas phase composition. Aims. The formation of molecules such as OH, H2 O, HO2, H2 O2, as well as their deuterated forms and O2 and O3 is studied in order to assess how the chemistry varies in different astrophysical environments, and how the gas phase is affected by grain surface chemistry. Methods. We use Monte Carlo simulations to follow the formation of molecules on bare grains as well as the fraction of molecules released into the gas phase. We consider a surface reaction network, based on gas phase reactions, as well as UV photo-dissociation of the chemical species. Results. We show that grain surface chemistry has a strong impact on gas phase chemistry, and that this chemistry is very different for different dust grain temperatures. Low temperatures favor hydrogenation, while higher temperatures favor oxygenation. Also, UV photons dissociate the molecules on the surface, that can reform subsequently. The formation-destruction cycle increases the amount of species released into the gas phase. We also determine the time scales to form ices in diffuse and dense clouds, and show that ices are formed only in shielded environments, as supported by observations.

Star formation is regulated through a variety of feedback processes. In this study, we treat feedback by metal injection and a UV background as well as by X-ray irradiation. Our aim is to investigate whether star formation is significantly affected when the ISM of a proto-galaxxy enjoys different metallicities and when a star forming cloud resides in the vicinity of a strong X-ray source. We perform cosmological Enzo simulations with a detailed treatment of non-zero metallicity chemistry and thermal balance. We also perform FLASH simulations with embedded Lagrangian sink particles of a collapsing molecular cloud near a massive, 10^{7} M\odot, black hole that produces X-ray radiation. We find that a multi-phase ISM forms for metallicites as small as 10^{-4} Solar at z = 6, with higher (10^{-2}Z\odot) metallicities supporting a cold (< 100 K) and dense (> 10^{3} cm^{-3}) phase at higher (z = 20) redshift. A star formation recipe based on the presence of a cold dense phase leads to a self-regulating mode in the presence of supernova and radiation feedback. We also find that when there is strong X-ray feedback a collapsing cloud fragments into larger clumps whereby fewer but more massive protostellar cores are formed. This is a consequence of the higher Jeans mass in the warm (50 K, due to ionization heating) molecular gas. Accretion processes dominate the mass function and a near-flat, non-Salpeter IMF results.

(abridged) When preplanetary bodies reach proportions of ~1 km or larger in size, their accretion rate is enhanced due to gravitational focusing (GF). We have developed a new numerical model to calculate the collisional evolution of the gravitationally-enhanced growth stage. We validate our approach against existing N-body and statistical codes. Using the numerical model, we explore the characteristics of the runaway growth and the oligarchic growth accretion phases starting from an initial population of single planetesimal radius R_0. In models where the initial random velocity dispersion (as derived from their eccentricity) starts out below the escape speed of the planetesimal bodies, the system experiences runaway growth. We find that during the runaway growth phase the size distribution remains continuous but evolves into a power-law at the high mass end, consistent with previous studies. Furthermore, we find that the largest body accretes from all mass bins; a simple two component approximation is inapplicable during this stage. However, with growth the runaway body stirs up the random motions of the planetesimal population from which it is accreting. Ultimately, this feedback stops the fast growth and the system passes into oligarchy, where competitor bodies from neighboring zones catch up in terms of mass. Compared to previous estimates, we find that the system leaves the runaway growth phase at a somewhat larger radius. Furthermore, we assess the relevance of small, single-size fragments on the growth process. In classical models, where the initial velocity dispersion of bodies is small, these do not play a critical role during the runaway growth; however, in models that are characterized by large initial relative velocities due to external stirring of their random motions, a situation can emerge where fragments dominate the accretion.

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

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

We discuss the formation of supermassive black holes in the early universe, and how to probe their subsequent evolution with the upcoming mm/sub-mm telescope ALMA. We first focus on the chemical and radiative conditions for black hole formation, in particular considering radiation trapping and molecular dissociation effects. We then turn our attention towards the magnetic properties in the halos where the first black holes form, and show that the presence of turbulence may lead to a magnetic dynamo, which could support the black hole formation process by providing an efficient means of transporting the angular momentum. We finally focus on observable properties of high-redshift black holes with respect to ALMA, and discuss how to distinguish between chemistry driven by the starburst and chemistry driven by X-rays from the black hole.

Molecules that trace the high-density regions of the interstellar medium may be used to evaluate the changing physical and chemical environment during the ongoing nuclear activity in (Ultra-)Luminous Infrared Galaxies. The changing ratios of the HCN(1-0), HNC(1-0), HCO+(1-0), CN(1-0) and CN(2-1), and CS(3-2) transitions were compared with the HCN(1-0)/CO(1-0) ratio, which is proposed to represent the progression time scale of the starburst. These diagnostic diagrams were interpreted using the results of theoretical modeling using a large physical and chemical network to describe the state of the nuclear ISM in the evolving galaxies. Systematic changes are seen in the line ratios as the sources evolve from early stage for the nuclear starburst (ULIRGs) to later stages. These changes result from changing environmental conditions and particularly from the lowering of the average density of the medium. A temperature rise due to mechanical heating of the medium by feedback explains the lowering of the ratios at later evolutionary stages. Infrared pumping may affect the CN and HNC line ratios during early evolutionary stages. Molecular transitions display a behavior that relates to changes of the environment during an evolving nuclear starburst. Molecular properties may be used to designate the evolutionary stage of the nuclear starburst. The HCN(1-0)/CO(1-0) and HCO+(1-0)/HCN(1-0) ratios serve as indicators of the time evolution of the outburst.

Accretion among macroscopic bodies of ~km size or larger is enhanced significantly due to gravitational focusing. Two regimes can be distinguished. Initially, the system experiences runaway growth, in which the gravitational focusing factors increase, and bodies at the high-mass tail of the distribution grow fastest. However, at some point the runaway body dynamically heats its environment, gravitational focusing factors decrease, and runaway growth passes into oligarchic growth. Based on the results of recent simulations, we reconsider the runaway growth-oligarchy transition. In contrast to oligarchy, we find that runaway growth cannot be approximated with a two component model (of small and large bodies) and that the criterion of Ida & Makino (1993), which is frequently adopted as the start of oligarchy, is not a sufficient condition to signify the transition. Instead, we propose a new criterion based on timescale arguments. We then find a larger value for the runaway growth-oligarchy transition: from several hundreds of km in the inner disk regions up to ~10^3 km for the outer disk. These findings are consistent with the view that runaway growth has been responsible for the size distribution of the present day Kuiper belt objects. Our finding furthermore outlines the proper initial conditions at the start of the oligarchy stage.

Photon-dominated regions (PDRs) are expected to show a layered structure in molecular abundances and emerging line emission, which is sensitive to the physical structure of the region as well as the UV radiation illuminating it. We aim to study this layering in the Orion Bar, a prototypical nearby PDR with a favorable edge-on geometry. We present new maps of 2 by 2 arcminute fields at 14-23 arcsecond resolution toward the Orion Bar in the SO 8_8-9_9, H2CO 5_(1,5)-4_(1,4), 13CO 3-2, C2H 4_(9/2)-3_(7/2) and 4_(7/2)-3_(5/2), C18O 2-1 and HCN 3-2 transitions. The data reveal a clear chemical stratification pattern. The C2H emission peaks close to the ionization front, followed by H2CO and SO, while C18O, HCN and 13CO peak deeper into the cloud. A simple PDR model reproduces the observed stratification, although the SO emission is predicted to peak much deeper into the cloud than observed while H2CO is predicted to peak closer to the ionization front than observed. In addition, the predicted SO abundance is higher than observed while the H2CO abundance is lower than observed. The discrepancies between the models and observations indicate that more sophisticated models, including production of H2CO through grain surface chemistry, are needed to quantitatively match the observations of this region.

In massive primordial galaxies, the gas may directly collapse and form a single central massive object if cooling is suppressed. Line cooling by molecular hydrogen can be suppressed in the presence of a strong soft-ultraviolet radiation field, but the role played by other cooling mechanisms is less clear. In optically thin gas, Lyman-Alpha cooling can be very effective, maintaining the gas temperature below 10^4 K over many orders of magnitude in density. However, the large neutral hydrogen column densities present in primordial galaxies render them highly optically thick to Lyman-Alpha photons. In this letter, we examine in detail the effects of the trapping of these Lyman-Alpha photons on the thermal and chemical evolution of the gas. We show that despite the high optical depth in the Lyman series lines, cooling is not strongly suppressed, and proceeds via other atomic hydrogen transitions, in particular the 2s-1s and the 3-2 transitions. At densities larger than 10^9 cm^{-3}, collisional dissociation of molecular hydrogen becomes the dominant cooling process and decreases the gas temperature to about 5000 K. The gas temperature evolves with density as $T \propto \rho^{\gamma_{\rm eff} - 1}$, with $\gamma_{\rm eff} = 0.97-0.98$. The evolution is thus very close to isothermal, and so fragmentation is possible, but unlikely to occur during the initial collapse. However, after the formation of a massive central object, we expect that later-infalling, higher angular momentum material will form an accretion disk that may be unstable to fragmentation, which may give rise to star formation with a top-heavy IMF.

(abridged) We explore how ALMA observations can probe high-redshift galaxies in unprecedented detail. We discuss the main observables that are excited by the large-scale starburst, and formulate expectations for the chemistry and the fluxes in the center of active galaxies, where chemistry may be driven by the absorption of X-ray photons. We show that such X-ray dominated regions (XDRs) should be large enough to be resolved with ALMA, and predict the expected amount of emission in CO and various fine-structure lines. We discuss how such XDRs can be distinguished from a strong starburst on the same spatial scales based on the CO line SED. Our models are compared to known sources like NGC 1068 and APM 08279. We also analyze the properties of the z=6.42 quasar SDSS J114816.64+525150.3, and find that the observed emission in CO, [CII] and [CI] requires a dense warm and a low-density cold gas component. We estimate the expected number of sources at redshifts higher than 6, finding that one could expect one black hole with $10^6$ solar masses per arcmin$^2$.

Properties of candidate stars, forming out of molecular clouds, depend on the ambient conditions of the parent cloud. We present a series of 2D and 3D simulations of fragmentation of molecular clouds in starburst regions as well as clouds under conditions in dwarf galaxies, leading to the formation of protostellar cores. We explore in particular the metallicity dependence of molecular cloud fragmentation and the possible variations in the dense core mass function, as the expression of a multi-phase ISM, due to dynamic and thermodynamic effects in starburst and metal-poor environments. To study the level of fragmentation during the collapse, the adaptive mesh refinement code FLASH is used. With this code, including self-gravity, thermal balance, turbulence and shocks, collapse is simulated with four different metallicity dependent cooling functions. Turbulent and rotational energies are considered as well. During the simulations, number densities of 10^8-10^9 cm^-3 are reached. The influences of dust and cosmic ray heating are investigated and a comparison to isothermal cases is made. The presented results indicate that fragmentation increases with metallicity, while cosmic ray and gas-grain collisional heating counteract this. We also find that modest rotation and turbulence can impact the cloud evolution as far as fragmentation is concerned. In this light, we conclude that radiative feedback, in starburst regions, will inhibit fragmentation, while low-metallicity dwarfs also should enjoy a star formation mode in which fragmentation is suppressed.

Over 300 extrasolar planets (exoplanets) have been detected orbiting nearby stars. We now hope to conduct a census of all planets around nearby stars and to characterize their atmospheres and surfaces with spectroscopy. Rocky planets within their star's habitable zones have the highest priority, as these have the potential to harbor life. Our science goal is to find and characterize all nearby exoplanets; this requires that we measure the mass, orbit, and spectroscopic signature of each one at visible and infrared wavelengths. The techniques for doing this are at hand today. Within the decade we could answer long-standing questions about the evolution and nature of other planetary systems, and we could search for clues as to whether life exists elsewhere in our galactic neighborhood.

We discuss the use of SPICA to study the cosmic history of star formation and accretion by supermassive black holes. The cooling lines, in particular the high-J rotational lines of CO, provide a clear-cut and unique diagnostic for separating the contributions of star formation and AGN accretion to the total infrared luminosity of active, gas-rich galaxies. We briefly review existing efforts for studying high-J CO emission from galaxies at low and high redshift. We finally comment on the detectability of cooling radiation from primordial (very low metallicity) galaxies containing an accreting supermassive black hole with SPICA/SAFARI.

Since the main cooling lines of the gas phase are important tracers of the interstellar medium in Galactic and extragalactic sources, proper and detailed understanding of their emission, and the ambient conditions of the emitting gas, is necessary for a robust interpretation of the observations. With high resolution (7"-9") maps (~3x3 pc^2) of mid-J molecular lines we aim to probe the physical conditions and spatial distribution of the warm (50 to few hundred K) and dense gas (n(H_2)>10^5 cm^-3) across the interface region of M17 SW nebula. We have used the dual color multiple pixel receiver CHAMP+ on APEX telescope to obtain a 5'.3x4'.7 map of the J=6-5 and J=7-6 transitions of 12CO, the 13CO J=6-5 line, and the {^3P_2}-{^3P_1} 370 um fine-structure transition of [C I] in M17 SW. LTE and non-LTE radiative transfer models are used to constrain the ambient conditions. The warm gas extends up to a distance of ~2.2 pc from the M17 SW ridge. The 13CO J=6-5 and [C I] 370 um lines have a narrower spatial extent of about 1.3 pc along a strip line at P.A=63 deg. The structure and distribution of the [C I] {^3P_2}-{^3P_1} 370 um map indicate that its emission arises from the interclump medium with densities of the order of 10^3 cm^-3. The warmest gas is located along the ridge of the cloud, close to the ionization front. An LTE approximation indicates that the excitation temperature of the embedded clumps goes up to ~120 K. The non-LTE model suggests that the kinetic temperature at four selected positions cannot exceed 230 K in clumps of density n(H_2)~5x10^5 cm^-3, and that the warm T_k>100 K and dense (n(H_2)>10^4 cm^-3) gas traced by the mid-J 12CO lines represent just about 2% of the bulk of the molecular gas. The clump volume filling factor ranges between 0.04 and 0.11 at these positions.

We have used the Odin satellite to obtain strip maps of the ground-state rotational transitions of ortho-water and ortho-ammonia, as well as CO(5-4) and 13CO(5-4) across the PDR, and H218O in the central position. A physi-chemical inhomogeneous PDR model was used to compute the temperature and abundance distributions for water, ammonia and CO. A multi-zone escape probability method then calculated the level populations and intensity distributions. These results are compared to a homogeneous model computed with an enhanced version of the RADEX code. H2O, NH3 and 13CO show emission from an extended PDR with a narrow line width of ~3 kms. Like CO, the water line profile is dominated by outflow emission, however, mainly in the red wing. The PDR model suggests that the water emission mainly arises from the surfaces of optically thick, high density clumps with n(H2)>10^6 cm^-3 and a clump water abundance, with respect to H2, of 5x10^-8. The mean water abundance in the PDR is 5x10^-9, and between ~2x10^-8 -- 2x10^-7 in the outflow derived from a simple two-level approximation. Ammonia is also observed in the extended clumpy PDR, likely from the same high density and warm clumps as water. The average ammonia abundance is about the same as for water: 4x10^-9 and 8x10^-9 given by the PDR model and RADEX, respectively. The similarity of water and ammonia PDR emission is also seen in the almost identical line profiles observed close to the bright rim. Around the central position, ammonia also shows some outflow emission although weaker than water in the red wing. Predictions of the H2O(110-101) and (111-000) antenna temperatures across the PDR are estimated with our PDR model for the forthcoming observations with the Herschel Space Observatory.

We present new high resolution numerical simulations of the ISM in a central R ~32 parsecs region around a supermassive black hole at a galactic center. Three-dimensional hydrodynamic modeling of the ISM (Wada & Norman 2002) with the nuclear starburst now includes tracking of the formation of molecular hydrogen out of the neutral hydrogen phase. In a quasi equilibrium state, mass fraction of H2 is about 0.4 (total H2 mass is ~1.5 10^6 Msun) of the total gas mass for the uniform far UV (FUV) with G_0 = 10. The gas forms an inhomogeneous disk, whose scale-height becomes larger in the outer region. H2 forms a thin nuclear disk in the inner ~ 5 pc, which is surrounded by molecular clouds swelled up toward h < 10 pc. The velocity field of the disk is highly turbulent in the torus region, whose velocity dispersion is ~ 20 km/s on average. Average supernova rate (SNR) of ~ 5 10^-5/yr is large enough to energize these structures. Gas column densities toward the nucleus larger than 10^22 cm^-2 are observed if the viewing angle is smaller than \theta_v ~ 50 deg from the edge-on. However, the column densities are distributed over almost two orders of magnitude around the average for any given viewing angle due to the clumpy nature of the torus. For a stronger FUV (G_0 =100), the total H2 mass in an equillibrium is only slightly smaller (~ 0.35). Finally the morphology and kinematics of the circumnuclear molecular gas disks emerging from our models is similar to that revealed by recent near infrared observations using VLTI/Keck.

With the incorporation of high-J molecular lines, we aim to constrain the physical conditions of the dense gas in the central region of the Seyfert 2 galaxy NGC 1068 and to determine signatures of the AGN or the starburst contribution. We used the James Clerk Maxwell Telescope to observe the J=4-3 transition of HCN, HNC, and HCO+, as well as the CN N_J=2_{5/2}-1_{3/2} and N_J=3_{5/2}-2_{5/2}, in NGC 1068. We estimate the excitation conditions of HCN, HNC, and CN, based on the line intensity ratios and radiative transfer models. We find that the bulk emission of HCN, HNC, CN, and the high-J HCO+ emerge from dense gas n(H_2)>=10^5 cm^-3). However, the low-J HCO+ lines (dominating the HCO+ column density) trace less dense (n(H_2)<10^5 cm^-3) and colder (T_K<=20 K) gas, whereas the high-J HCO+ emerges from warmer (>30 K) gas than the other molecules. The HCO+ J=4-3 line intensity, compared with the lower transition lines and with the HCN J=4-3 line, support the influence of a local XDR environment. The estimated N(CN)/N(HCN)~1-4 column density ratios are indicative of an XDR/AGN environment with a possible contribution of grain-surface chemistry induced by X-rays or shocks.

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

Context: The HD and H2 molecules play important roles in the cooling of primordial and very metal-poor gas at high redshift. Aims: Grain surface and gas phase formation of HD and H2 is investigated to assess the importance of trace amounts of dust, 10^{-5}-10^{-3} Zo, in the production of HD and H2. Methods: We consider carbonaceous and silicate grains and include both physisorption and chemisorption, tunneling, and realistic grain surface barriers. We find, for a collapsing gas cloud environment with coupled chemical and thermal balance, that dust abundances as small as 10^{-5} solar lead to a strong boost in the H2 formation rate due to surface reactions. As a result of this enhancement in H2, HD is formed more efficiently in the gas phase through the D+ + H2 reaction. Direct formation of HD on dust grains cannot compete well with this gas phase process for dust temperatures below 150 K. We also derive up-to-date analytic fitting formulae for the grain surface formation of H2 and HD, including the different binding energies of H and D. Results: Grain surface reactions are crucial to the availability of H2 and HD in very metal-poor environments. Above metallicities of 10^{-5} solar, the grain surface route dominates the formation of H2, which in turn, drives the formation of HD in the gas phase. At dust temperatures above 150 K, laboratory experiments and theoretical modelling suggest that H2 formation on grains is suppressed while HD formation on grains is not.

We study the prospects of finding the first quasars in the universe with ALMA and JWST. For this purpose, we derive a model for the high-redshift black hole population based on observed relations between the black hole mass and the host galaxy. We re-address previous constraints from the X-ray background with particular focus on black hole luminosities below the Eddington limit as observed in many local AGN. For such luminosities, up to 20% of high-redshift black holes can be active quasars. We then discuss the observables of high-redshift black holes for ALMA and JWST by adopting NGC 1068 as a reference system. We calculate the expected flux of different fine-structure lines for a similar system at higher redshift, and provide further predictions for high-J CO lines. We discuss the expected fluxes from stellar light, the AGN continuum and the Lyman $\alpha$ line for JWST. Line fluxes observed with ALMA can be used to derive detailed properties of high-redshift sources. We suggest two observational strategies to find potential AGN at high redshift and estimate the expected number of sources, which is between 1-10 for ALMA with a field of view of $\sim(1')^2$ searching for line emission and 100-1000 for JWST with a field of view of $(2.16')^2$ searching for continuum radiation. We find that both telescopes can probe high-redshift quasars down to redshift 10 and beyond, and therefore truely detect the first quasars in the universe.

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

Circumstellar disks provide the material reservoir for the growth of young stars and for planet formation. We combine a high-level radiative transfer program with a thermal-chemical model of a typical T Tauri star disk to investigate the diagnostic potential of the far-infrared lines of water for probing disk structure. We discuss the observability of pure rotational H2O lines with the Herschel Space Observatory, specifically the residual gas where water is mainly frozen out. We find that measuring both the line profile of the ground 110-101 ortho-H2O transition and the ratio of this line to the 312-303 and 221-212 line can provide information on the gas phase water between 5-100 AU, but not on the snow line which is expected to occur at smaller radii.

Aims: Molecular emission lines originating in the nuclei of luminous infra-red galaxies are used to determine the physical properties of the nuclear ISM in these systems. Methods: A large observational database of molecular emission lines is compared with model predictions that include heating by UV and X-ray radiation, mechanical heating, and the effects of cosmic rays. Results: The observed line ratios and model predictions imply a separation of the observedsystems into three groups: XDRs, UV-dominated high-density (n>=10^5 cm-3) PDRs, and lower-density (n=10^4.5 cm-3) PDRs that are dominated by mechanical feedback. Conclusions: The division of the two types of PDRs follows naturally from the evolution of the star formation cycle of these sources, which evolves from deeply embedded young stars, resulting in high-density (n>=10^5 cm-3) PDRs, to a stage where the gas density has decreased (n=10^4.5 cm-3) and mechanical feedback from supernova shocks dominates the heating budget.

Infrared Dark Clouds (IRDCs) represent the earliest observed stages of clustered star formation, characterized by large column densities of cold and dense molecular material observed in silhouette against a bright background of mid-IR emission. Up to now, IRDCs were predominantly known toward the inner Galaxy where background infrared emission levels are high. We present Spitzer observations with the Infrared Camera Array toward object G111.80+0.58 (G111) in the outer Galactic Plane, located at a distance of ~3 kpc from us and ~10 kpc from the Galactic center. Earlier results show that G111 is a massive, cold molecular clump very similar to IRDCs. The mid-IR Spitzer observations unambiguously detect object G111 in absorption. We have identified for the first time an IRDC in the outer Galaxy, which confirms the suggestion that cluster-forming clumps are present throughout the Galactic Plane. However, against a low mid-IR back ground such as the outer Galaxy it takes some effort to find them.

An overview is given of the chemical processes that occur in primordial systems under the influence of radiation, metal abundances and dust surface reactions. It is found that radiative feedback effects differ for UV and X-ray photons at any metallicity, with molecules surviving quite well under irradiation by X-rays. Starburst and AGN will therefore enjoy quite different cooling abilities for their dense molecular gas. The presence of a cool molecular phase is strongly dependent on metallicity. Strong irradiation by cosmic rays (>100x the Milky Way value) forces a large fraction of the CO gas into neutral carbon. Dust is important for H2 and HD formation, already at metallicities of 10^-4-10^-3 solar, for electron abundances below 10^-3.

We present a method which extends Monte Carlo studies to situations that require a large dynamic range in particle number. The underlying idea is that, in order to calculate the collisional evolution of a system, some particle interactions are more important than others and require more resolution, while the behavior of the less important, usually of smaller mass, particles can be considered collectively. In this approximation groups of identical particles, sharing the same mass and structural parameters, operate as one unit. The amount of grouping is determined by the zoom factor -- a free parameter that determines on which particles the computational effort is focused. Two methods for choosing the zoom factors are discussed: the `equal mass method,' in which the groups trace the mass density of the distribution, and the `distribution method,' which additionally follows fluctuations in the distribution. Both methods achieve excellent correspondence with analytic solutions to the Smoluchowski coagulation equation. The grouping method is furthermore applied to simulations involving runaway kernels, where the particle interaction rate is a strong function of particle mass, and to situations that include catastrophic fragmentation. For the runaway simulations previous predictions for the decrease of the runaway timescale with the initial number of particles ${\cal N}$ are reconfirmed, extending ${\cal N}$ to $10^{160}$. Astrophysical applications include modeling of dust coagulation, planetesimal accretion, and the dynamical evolution of stars in large globular clusters. The proposed method is a powerful tool to compute the evolution of any system where the particles interact through discrete events, with the particle properties characterized by structural parameters.

Many present-day galaxies are known to harbor supermassive, >10^6 Mo, black holes. These central black holes must have grown through accretion from less massive seeds in the early universe. The molecules CO and H2 can be used to trace this young population of accreting massive black holes through the X-ray irradiation of ambient gas. The X-rays drive a low-metallicity ion-molecule chemistry that leads to the formation and excitation of CO and H2 in 100<T<1,000 K gas. H2 traces very low metallicity gas, ~10^-3 solar or less, while some pollution by metals, ~10^-2 solar or more, must have taken place to form CO. Strong CO J>15 and H2 S(0) & S(1) emission is found that allows one to constrain ambient conditions. Comparable line strengths cannot be produced by FUV or cosmic ray irradiation. Weak, but perhaps detectable, H3+ (2,2) - (1,1) emission is found and discussed. The models predict that black hole masses larger than 10^5 Mo can be detected with ALMA, over a redshift range of 5-20, provided that the black holes radiate close to Eddington.

The nuclei of ULIRGs harbor massive young stars, an accreting central black hole, or both. Results are presented for molecular gas that is exposed to X-rays (1-100 keV, XDRs) and far-ultraviolet radiation (6-13.6 eV, PDRs). Attention is paid to species like HCO+, HCN, HNC, OH, H2O and CO. Line ratios of HCN/HCO+ and HNC/HCN discriminate between PDRs and XDRs. Very high J (>10) CO lines, observable with HIFI/Herschel, discriminate very well between XDRs and PDRs. In XDRs, it is easy to produce large abundances of warm (T>100 K) H2O and OH. In PDRs, only OH is produced similarly well.

The observed molecular properties of a sample of FIR-luminous and OH megamaser (OH-MM) galaxies have been investigated. The ratio of high and low-density tracer lines is found to be determined by the progression of the star formation in the system. The HCO+/HCN and HCO+/HNC line ratios are good proxies for the density of the gas, and PDR and XDR sources can be distinguished using the HNC/HCN line ratio. The properties of the OH-MM sources in the sample can be explained by PDR chemistry in gas with densities higher than 10^5.5 cm^-3, confirming the classical OH-MM model of IR pumped amplification with (variable) low gains.

In the inner parts of the Galaxy the Infrared Dark Clouds (IRDCs) are presently believed to be the progenitors of massive stars and star clusters. Many of them are predominantly devoid of active star formation and for now they represent the earliest observed stages of massive star formation. Their Outer Galaxy counterparts, if present, are not easily identified because of a low or absent mid-IR background. We characterize the ambient conditions in the Outer Galaxy IRDC candidate G111.80+0.58, a relatively quiescent molecular core complex in the vicinity of NGC7538. We conduct molecular line observations on a number of dense cores and analyze the data in terms of excitation temperature, column and volume density, mass and stability. The temperatures (15-20K) are higher than expected from only cosmic ray heating, but comparable to those found in massive cores. Star forming activity could be present in some cores, as indicated by the presence of warm gas and YSO candidates. The observed super-thermal line-widths are typical for star forming regions. The velocity dispersion is consistent with a turbulent energy cascade over the observed size scales. We do not find a correlation between the gas temperature and the line-width. The LTE masses we find are much larger than the thermal Jeans mass and fragmentation is expected. In that case the observed lines represent the combined emission of multiple unresolved components. We conclude that G111.80+0.58 is a molecular core complex with bulk properties very similar to IRDCs in an early, but not pristine, star forming state. The individual cores are close to virial equilibrium and some contain sufficient material to form massive stars and star clusters. The ambient conditions suggest that turbulence is involved in supporting the cores against gravitational collapse.

This paper describes a James Clerk Maxwell Telescope (JCMT) legacy survey that has been awarded roughly 500 hrs of observing time to be carried out from 2007 to 2009. In this survey we will map with SCUBA-2 (Submillimetre Common User Bolometer Array 2) almost all of the well-known low-mass and intermediate-mass star-forming regions within 0.5 kpc that are accessible from the JCMT. Most of these locations are associated with the Gould Belt. From these observations we will produce a flux-limited snapshot of star formation near the Sun, providing a legacy of images, as well as point-source and extended-source catalogues, over almost 700 square degrees of sky. The resulting images will yield the first catalogue of prestellar and protostellar sources selected by submillimetre continuum emission, and should increase the number of known sources by more than an order of magnitude. We will also obtain CO maps with the array receiver HARP (Heterodyne Array Receiver Programme), in three CO isotopologues, of a large typical sample of prestellar and protostellar sources. We will then map the brightest hundred sources with the SCUBA-2 polarimeter (POL-2), producing the first statistically significant set of polarization maps in the submillimetre. The images and source catalogues will be a powerful reference set for astronomers, providing a detailed legacy archive for future telescopes, including ALMA, Herschel and JWST.

Based on its strong CO+ emission it is argued that the M 82 star-burst galaxy is exposed to a combination of FUV and X-ray radiation. The latter is likely to be the result of the star-burst superwind, which leads to diffuse thermal emission at ~0.7 keV, and a compact hard, 2-10 keV, source (but not an AGN). Although a photon-dominated region (FUV) component is clearly present in the nucleus of M 82, and capable of forming CO+, only X-ray irradiated gas of density 10^3-10^5 cm^-3 can reproduce the large, ~(1-4)x10^13 cm^-2, columns of CO+ that are observed toward the proto-typical star-burst M 82. The total X-ray luminosity produced by M 82 is weak, ~10^41 erg s^-1, but this is sufficient to drive the formation of CO+.

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.

Context: Existing SWAS observations and future HIFI/Herschel data require a clear sense of the information content of water emission and absorption lines. Aims: Investigate wether the ground-state transition of ortho-H2O, 1_(10)-1_(01), at 557GHz can be used to measure the column density throughout an interstellar cloud. Methods: We make use of a multi-zone escape probability code suitable for the treatment of molecular line emission. Results: For low abundances, i.e., X(H2O)=<10^(-9), the intensity of the 1_(10)-1_(01) transition scales with the total column of H2. However, this relationship breaks down with increasing abundance, i.e., optical depth, due to line trapping and -- for T_(dust)>=25K, X(H2O)=<10^(-8) and n~10^4 cm^(-3) -- absorption of the dust continuum. Conclusions: An observed decline in intensity per column, expected if H2O is a surface tracer, does not necessarily mean that the water is absent in the gas phase at large columns, but can be caused by line trapping and subsequent collisional de-excitation. To determine the amount of water vapour in the interstellar medium, multiple line measurements of optically thin transitions are needed to disentangle radiative transfer and local excitation effects.

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

Stars form in the densest, coldest, most quiescent regions of molecular clouds. Molecules provide the only probes which can reveal the dynamics, physics, chemistry and evolution of these regions, but our understanding of the molecular inventory of sources and how this is related to their physical state and evolution is rudimentary and incomplete. The Spectral Legacy Survey (SLS) is one of seven surveys recently approved by the JCMT Board. Starting in 2007, the SLS will produce a spectral imaging survey of the content and distribution of all the molecules detected in the 345 GHz atmospheric window (between 332 GHz and 373 GHz) towards a sample of 5 sources. Our intended targets are: a low mass core (NGC1333 IRAS4), 3 high mass cores spanning a range of star forming environments and evolutionary states (W49, AFGL2591, and IRAS20126), and a PDR (the Orion Bar). The SLS will use the unique spectral imaging capabilities of HARP-B/ACSIS to study the molecular inventory and the physical structure of these objects, which span different evolutionary stages and physical environments, to probe their evolution during the star formation process. As its name suggests, the SLS will provide a lasting data legacy from the JCMT that is intended to benefit the entire astronomical community. As such, the entire data set (including calibrated spectral datacubes, maps of molecular emission, line identifications, and calculations of the gas temperature and column density) will be publicly available.

Understanding the origin of stellar masses is a key problem in astrophysics. In the solar neighborhood, the mass distribution of stars follows a seemingly universal pattern. In the centre of the Milky Way, however, there are indications for strong deviations and the same may be true for the nuclei of distant starburst galaxies. Here we present the first numerical hydrodynamical calculations of stars formed in a molecular region with chemical and thermodynamic properties similar to those of warm and dusty circum-nuclear starburst regions. The resulting IMF is top-heavy with a peak at ~ 15 Msun, a sharp turn-down below 7 Msun and a power-law decline at high masses. We find a natural explanation for our results in terms of the temperature dependence of the Jeans mass, with collapse occuring at a temperature of ~ 100 K and an H2 density of a few times 10^5 cm^-3, and discuss possible implications for galaxy formation and evolution.

\abridged The nuclei of active galaxies harbor massive young stars, an accreting central black hole, or both. In order to determine the physical conditions that pertain to molecular gas close to the sources of radiation, numerical models are constructed. These models determine the thermal and chemical balance of molecular gas that is exposed to X-rays and far-ultraviolet radiation, as a function of depth. We present a grid of XDR and PDR models that span ranges in density, irradiation and column density. We find that the fine-structure line ratios of e.g. [SiII] 35mum/[CII] 158 mum are larger in XDRs than in PDRs, for a given density, column and irradiation strength. We find that the line ratios HCN/HCO+ and HNC/HCN, as well as the column density ratio CN/HCN, discriminate between PDRs and XDRs. The HCN/HCO+ 1-0 ratio is <1 (>1) for XDRs (PDRs) if the density exceeds 10^5 cm^-3 and if the column density is >10^23 cm^-2. For columns <10^ 22.5 cm^-2 the XDR HCN/HCO+ 1-0 ratio becomes larger than one, although the individual HCN 1-0 and HCO+ 1-0 line intensities are weaker. For modest densities, n=10^4-10^5 cm^-3, and strong radiation fields (>100 erg s^-1 cm^-2), HCN/HCO+ ratios can become larger in XDRs than PDRs as well. The HCN/CO 1-0 ratio is typically smaller in XDRs, and the HCN emission in XDRs is boosted with respect to CO only for columns >10^{23} cm^{-2} and densities <10^4 cm^-3. CO is typically warmer in XDRs than in PDRs, for the same total energy input. This leads to higher CO J=N+1-N/CO 1-0, N>=1, line ratios in XDRs. Lines with N>=10, like CO(16-15) and CO(10-9) observable with HIFI/Herschel, discriminate very well between XDRs and PDRs. Column density ratios indicate that CH, CH+, NO, HOC+ and HCO are good PDR/XDR discriminators.

Context: Sticking of colliding dust particles through van der Waals forces is the first stage in the grain growth process in protoplanetary disks, eventually leading to the formation of comets, asteroids and planets. A key aspect of the collisional evolution is the coupling between dust and gas motions, which depends on the internal structure (porosity) of aggregates. Aims: To quantify the importance of the internal structure on the collisional evolution of particles, and to create a new coagulation model to investigate the difference between porous and compact coagulation in the context of a turbulent protoplanetary disk. Methods: We have developed simple prescriptions for the collisional evolution of porosity of grain-aggregates in grain-grain collisions. Three regimes can then be distinguished: `hit-and-stick' at low velocities, with an increase in porosity; compaction at intermediate velocities, with a decrease of porosity; and fragmentation at high velocities. (..) Results: (..) We can discern three different stages in the particle growth process (..) We find that when compared to standard, compact models of coagulation, porous growth delays the onset of settling, because the surface area-to-mass ratio is higher, a consequence of the build-up of porosity during the initial stages. As a result, particles grow orders of magnitudes larger in mass before they rain-out to the mid-plane. Depending on the turbulent viscosity and on the position in the nebula, aggregates can grow to (porous) sizes of ~ 10 cm in a few thousand years. We also find that collisional energies are higher than in the limited PCA/CCA fractal models, thereby allowing aggregates to restructure. It is concluded that the microphysics of collisions plays a key role in the growth process.

The ISM of active galaxy centers is exposed to a combination of cosmic ray, FUV and X-ray radiation. We apply PDR models to this ISM with both `normal' and highly elevated (5\times 10^{-15}s^-1) cosmic-ray rates and compare the results to those obtained for XDRs. Our existing PDR-XDR code is used to construct models over a 10^3-10^5 cm^-3 density range and for 0.16-160 erg s^-1 cm^-2 impingent fluxes. We obtain larger high J (J>10) CO ratios in PDRs when we use the highly elevated cosmic ray rate, but these are always exceeded by the corresponding XDR ratios. The [CI] 609 mum/13CO(2-1) line ratio is boosted by a factor of a few in PDRs with n~10^3 cm^-3 exposed to a high cosmic ray rate. At higher densities ratios become identical irrespective of cosmic ray flux, while XDRs always show elevated [CI] emission per CO column. The HCN/CO and HCN/HCO+ line ratios, combined with high J CO emission lines, are good diagnostics to distinguish between PDRs under either low or high cosmic ray irradiation conditions, and XDRs. Hence, the HIFI instrument on Herschel, which can detect these CO lines, will be crucial in the study of active galaxies.

The polytropic equation of state of an atomic hydrogen gas is examined for primordial halos with baryonic masses of M_h~10^7-10^9 Mo. For roughly isothermal collapse around 10^4 K, we find that line trapping of Lyman alpha (HI and HeII) photons causes the polytropic exponent to stiffen to values significantly above unity. Under the assumptions of zero H2 abundance and very modest pollution by metals (<10^-4 Solar), fragmentation is likely to be inhibited for such an equation of state. We argue on purely thermodynamic grounds that a single black hole of ~0.02-0.003M_h can form at the center of a halo for z=10-20 when the free-fall time is less than the time needed for a resonantly scattered Lyman alpha photon to escape from the halo. The absence of H2 follows naturally from the high, 10^4 K, temperatures that are attained when Lyman alpha photons are trapped in the dense and massive halos that we consider. An H2 dissociating UV background is needed if positive feedback effects on H2 formation from X-rays occur. The black hole to baryon mass fraction is suggestively close to what is required for these intermediate mass black holes, of mass M_BH~10^4-10^6 Mo, to act as seeds for forming the supermassive black holes of mass ~0.001M_spheroid found in galaxies today.

An evolutionary model of star formation in ULIRGs was created using existing star formation and radiative transfer codes (STARBURST99, RADMC and RADICAL) as building blocks. The results of the simulations are compared to data from two IRAS catalogs. From the simulations it is found that the dust surrounding the starburst region is made up from two components. There is a low optical depth (tau=0.1, which corresponds to 0.1% of the total dust mass), hot (T~400K) non-grey component close to the starburst (scale size 10pc) and a large scale, colder grey component (100pc, 75K) with a much larger column (tau=10). The simulations also show that starburst galaxies can be powered by massive star formation. The parameters for this star forming region are difficult to determine, since the IR continuum luminosity is only sensitive to the total UV input. Therefore, there is a degeneracy between the total starburst mass and the initial mass function (IMF) slope. A less massive star formation with a shallower IMF will produce the same amount of OB stars and therefore the same amount of irradiating UV flux. Assuming the stars are formed according to a Salpeter IMF (Psi(M) ~ M^-2.35), the star formation region should produce 10^9 Msun of stars (either in one instantaneous burst, or in a continuous process) in order to produce enough IR radiation.

We model the spectra and surface brightness distributions for the Lyman alpha (Lya) radiation expected from protogalaxies that are caught in the early stages of their assembly. We use the results of a companion paper to characterize the radiation emerging from spherically collapsing gas clouds. We then modify the intrinsic spectra to incorporate the effect of subsequent resonant scattering in the intergalactic medium (IGM). Using these models, we interpret a number of recent observations of extended Lya blobs (LABs) at high redshift. We suggest, based on the angular size, energetics, as well as the relatively shallow surface brightness profiles, and double-peaked spectra, that several of these LABs may be associated with collapsing protogalaxies. We suggest two follow-up observations to diagnose the presence of gas infall. High S/N spectra of LABs should reveal a preferential flattening of the surface brightness profile at the red side of the line. Complementary imaging of the blobs at redshifted Balmer alpha wavelengths should reveal the intrinsic Lya emissivity and allow its separation from radiative transfer effects. We show that Lya scattering by infalling gas can reproduce the observed spectrum of Steidel et al's LAB2 as accurately as a recently proposed outflow model. Finally, we find similar evidence for infall in the spectra of point-like Lyman alpha emitters. The presence of scattering by the infalling gas implies that the intrinsic Lya luminosities, and derived quantities, such as the star-formation rate, in these objects may have been underestimated by about an order of magnitude.

We present Monte Carlo calculations of Lyman alpha (Lya) radiative transfer through collapsing gas clouds, representing proto-galaxies that are caught in the process of their assembly. Such galaxies produce Lya flux over an extended solid angle from a spatially extended Lya emissivity and/or from scattering effects. We study the effect of the gas distribution and kinematics, and of the Lya emissivity profile, on the emergent spectrum and surface brightness distribution. The emergent Lya spectrum is typically double-peaked and asymmetric. In practice, the blue peak is significantly enhanced and the red peak, in most cases, will be undetectable. The resulting effective blueshift, combined with scattering in the intergalactic medium, will render extended Lya emission from collapsing protogalaxies difficult to detect beyond redshift z=4. The surface brightness distribution is typically flat, and a strong wavelength dependence of its slope (with preferential flattening at the red side of the line) would be a robust indication that Lya photons are being generated (rather than just scattered) in a spatially extended region around the galaxy. We also find that for self-ionized clouds whose effective Lya optical depth is less than 10^3, infall and outflow models can produce nearly identical spectra and surface brightness distributions, and are difficult to distinguish from one another. The presence of deuterium with a cosmic abundance may produce a narrow but detectable dip in the spectra of systems with moderate hydrogen column densities, in the range 10^18-10^20 cm^-2. Finally, we present a new analytic solution for the emerging Lya spectrum in the limiting case of a static uniform sphere, extending previous solutions for static plane-parallel slabs.

A radiative transfer method for the treatment of molecular lines is presented. We apply this method to previous SWAS and ISO observations of water vapor in the source S140 in order to make models to plan for, and to interpret, HIFI data. Level populations are calculated with the use of a three-dimensional (multi-zone) escape probability method and with a long characteristic code that uses Monte Carlo techniques with fixed directions. Homogeneous and inhomogeneous models are used to compute the differences between water line profiles across the S140 region. We find that when an outflow or infall velocity field with a gradient of a few kms^{-1} is adopted, line profiles with a FWHM of 6 kms^{-1} are found, in agreement with observations. Inhomogeneous models are favoured to produce a single-peaked line profile. When zooming in on smaller regions within the PDR, the shapes of the line profiles start to differ due to the different temperature and density distributions there. The embedded sources are traced by high excitation lines of, e.g., 3_{21}-2_{21}, 3_{03}-2_{12}, 2_{12}-1_{01} and 2_{20}-1_{11}. The computed intensities are roughly consistent with existing ISO observations. Water emission in a PDR source like S140 requires a combination of a pure PDR and an embedded source in order to match the observations. Because of its good angular resolution, HIFI will be able to distinguish between a dense star forming region or a more diffuse gas component. It is therefore important for future observing programs to consider both in their predictions of the emission characteristics of water in these environments.

Stars form by gravoturbulent fragmentation of interstellar gas clouds. The supersonic turbulence ubiquitously observed in Galactic molecular gas generates strong density fluctuations with gravity taking over in the densest and most massive regions. Collapse sets in to build up stars and star clusters. Turbulence plays a dual role. On global scales it provides support, while at the same time it can promote local collapse. Stellar birth is thus intimately linked to the dynamic behavior of parental gas clouds, which governs when and where protostellar cores form, and how they contract and grow in mass via accretion from the surrounding cloud material to build up stars. The equation of state plays a pivotal role in the fragmentation process. Under typical cloud conditions, massive stars form as part of dense clusters following the "normal" mass function observed, e.g. in the solar neighborhood. However, for gas with an effective polytropic index greater than unity star formation becomes biased towards isolated massive stars. This is relevant for understanding the properties of zero-metallicity stars (Population III) or stars that form under extreme environmental conditions like in the Galactic center or in luminous starbursts.

We consider the excitation of water in the Photon Dominated Region (PDR). With the use of a three-dimensional escape probability method we compute the level populations of ortho- and para-H_2O up to 350 K (i.e., 8 levels), as well as line intensities for various transitions. Homogeneous and inhomogeneous models are presented with densities of 10^4-10^5 cm^{-3} and the differences between the resulting intensities are displayed. Density, temperature, and abundance distributions inside the cloud are computed with the use of a self-consistent physi-chemical (in)homogeneous model in order to reproduce the line intensities observed with SWAS, and to make predictions for various lines that HIFI will probe in the future. Line intensities vary from 10^{-13} erg cm^{-2} s^{-1} sr^{-1} to a few times 10^{-6} erg cm^{-2} s^{-1} sr^{-1}. We can reproduce the intensity for the 1_{10}-1_{01} line observed by the SWAS satellite. It is found that the 2_{12}-1_{01} line is the strongest, whereas the 3_{12}-2_{21} line is the weakest, in all the models. It is found that the 1_{10}-1_{01} line probes the total column, while higher excitation lines probe the higher density gas (e.g., clumps).

NGC 4258 is well known for its 22 GHz water maser features that appear to reside in a thin Keplerian disk orbiting at R_0=0.13 pc. The influence of gravitational redshift on the level of maser saturation and the line strength and width of systemic water maser emission is studied. The effective gravitational velocity shift along the gain path through the disk amounts to about twice the thermal velocity width of H2O molecules at 500 K for a masing ring of thickness 0.2R_0. This GR induced frequency shift has a sign opposite to the shift caused by the projected velocity difference vector of Keplerian motion on the blue side of the systemic maser emission and has the same sign on the red systemic side. Consequently, the negative optical depth that a line photon sees as it traverses the masing disk along a path of minimal Keplerian velocity gradient is larger on the blue side. Note that one is not concerned here with an overall shift in the maser spectrum caused by the fact that the masing ring as a whole is located in the black hole's potential well, i.e., it is the change in the gravitational potential along the gain path that matters. It is found, for a Schwarzschild metric, that saturation is enhanced (suppressed) for systemic blue-shifted (red-shifted) maser features by gravitational redshift of line photons. Approaching peak flux systemic features are, on average, brighter, due to saturation, and wider (~2 km/s), due to rebroadening, than their receding counterparts (~1 km/s). The effects are consistent with observations that span several years, i.e., over periods longer than needed for disk material to move through the systemic emission region. Further confirmation of this effect can follow by combining all available systemic 22 GHz data over the past two decades (abridged).

The polytropic equation of state (EOS) of primordial gas clouds with modest enrichment is computed, motivated by the recent observations of very Fe-deficient stars, [Fe/H]~10^{-3.5}-10^{-5}, such as HE 0107-5240 and CS 29498-043. These stars are overabundant, relative to Fe, in C and O. We assume that the observed abundances of species like C, O, Si and Fe are representative of the gas from which the currently observed metal-deficient stars formed. Under this assumption, we find that this primordial metal abundance pattern has profound consequences for the thermal balance and chemical composition of the gas, and hence for the EOS of the parental cloud. The polytopic EOS is soft for low, [O/H]<10^{-3}, oxygen abundances, but stiffens to a polytropic index gamma large than unity for [O/H]>10^{-2} due to the large opacity in the CO and H2O cooling lines. It is further found that a regulating role is played by the presence and temperature of the dust, even when the overall carbon abundance is only [C/H]~10^{-2}. When the dust is warmer than the gas, a region with gamma~1.2 results around a density of ~10^4 cm^{-3}. When the dust is colder than the gas, a region with gamma~0.8 is found for a density of ~10^6 cm^{-3}. Implications for the primordial initial mass function (IMF) as well as the IMF in starburst galaxies, where the metallicity is super-solar, are explored and related to processes that influence the temperature of the ambient dust.

We present a far-ultraviolet (PDR) and an X-ray dominated region (XDR) code. We include and discuss thermal and chemical processes that pertain to irradiated gas. An elaborate chemical network is used and a careful treatment of PAHs and H2 formation, destruction and excitation is included. For both codes we calculate four depth-dependent models for different densities and radiation fields, relevant to conditions in starburst galaxies and active galactic nuclei. A detailed comparison between PDR and XDR physics is made for total gas column densities between ~10^20 and ~10^25 cm^-2. We show cumulative line intensities for a number of fine-structure lines (e.g., [CII], [OI], [CI], [SiII], [FeII]), as well as cumulative column densities and column density ratios for a number of species (e.g., CO/H2, CO/C, HCO+/HCN, HNC/HCN). The comparison between the results for the PDRs and XDRs shows that column density ratios are almost constant up to N_H=10^22 cm^-2 for XDRs, unlike those in PDRs. For example, CO/C in PDRs changes over four orders of magnitude from the edge to N_H=10^22 cm^-2. The CO/C and CO/H2 ratios are lower in XDRs at low column densities and rise at N_H > 10^23 cm^-2. At most column densities N_H > 10^21.5 cm^-2, HNC/HCN ratios are lower in XDRs too, but they show a more moderate increase at higher N_H.

We study the formation of molecular hydrogen on dust grain surfaces and apply our results to the high redshift universe. We find that a range of physical parameters, in particular dust temperature and gas temperature, but not so much dust surface composition, influence the formation rate of H$_2$. The H$_2$ formation rate is found to be suppressed above gas kinetic temperatures of a few hundred K and for dust temperatures above 200-300 K and below 10 K. We highlight the differences between our treatment of the H$_2$ formation process and other descriptions in the literature. We also study the relative importance of H$_2$ formation on dust grains with respect to molecular hydrogen formation in the gas phase, through the H$^-$ route. The ratio of formation rates of these two routes depends to a large part on the dust abundance, on the electron abundance, and also on the relative strength of the FUV (extra-)galactic radiation field. We find that for a cosmological evolution of the star formation rate and dust density consistent with the Madau plot, a positive feedback effect on the abundance of H$_2$ due to the presence of dust grains can occur at redshifts $z\ge 3$. This effect occurs for a dust-to-gas mass ratio as small as 10$^{-3}$ of the galactic value.

A prediction of standard inflationary cosmology is that the elemental composition of the medium out of which the earliest stars and galaxies condensed consisted primarily of hydrogen and helium 4He with small admixtures of deuterium, lithium 7Li, and 3He. The most red-shifted quasars, galaxies, and Ly-alpha absorbers currently observed, however, all exhibit at least some admixture of heavier elements, as do the most ancient stars in the Galaxy. Here we examine ways in which the abundance of these same elements, if present before the epoch of population III formation, might be observationally established or ruled out.

We studied the HST-WFPC2 observations of the inner kpc of the interacting galaxy M51 in six bands from 2550 to 8140 Angstrom. The images show an oval shaped area ("bulge") of 11x16 arcsec or 450x650 pc around the nucleus, dominated by a smooth population of old stars with overimposed dust lanes. We found 30 bright point-like sources in the bulge of M51; many of these are bright in the UV. They are located in elongated "strings" which follow the general pattern of the dust lanes. The spectral energy distributions of the pointsources are compared with those of clusters or single stars. We present arguments that the pointsources are isolated massive stars or small groups of a few massive stars. The energy distributions are best fitted with those of single stars with 4000<Teff<50000 K, 4.2<log L/Lsun<7.2, and 12<Mstar<200 Msun. We have derived upper limits to the total mass of low mass stars that could be "hiding" within the point sources. For the bluest sources the upper limit is only a few hundred Msun. We conclude that massive stars are formed outside clusters (or in very low mass clusters) in the bulge of M51. The estimated star formation rate in the bulge is 1 to 2x10^(-3) Msun/yr. This star formation rate could be sustained for about 200 to 400 Myrs, suggesting that the ongoing massive star formation in the bulge of M51 is fed/triggered by the interaction with its companion about 400 Mrs ago. Theoretical cloud-models show that CO is destroyed in low extinction clouds near the bright nucleus of M51, so that H_2, [OI] and [CII] are the dominants coolants. This condition favors the formation of massive stars, in agreement with the observations. This mode of massive-star formation in the bulge of M51 may resemble the star formation in the early Universe, when the CO and dust contents were low due to the low metallicity.

Future observations with ALMA will be able to detect star-forming primordial gal axies, and perhaps even their dwarf spheroidal satellites, in CO emission lines at redshifts approaching 10. These observations will compliment other tools designed to study the dawn of galaxy formation, such as NGST and FIRST.

Recent observations with the Submillimeter Wave Astronomy Satellite indicate abundances of gaseous H2O and O2 in dense molecular clouds which are significantly lower than found in standard homogeneous chemistry models. We present here results for the thermal and chemical balance of inhomogeneous molecular clouds exposed to ultraviolet radiation in which the abundances of H2O and O2 are computed for various density distributions, radiation field strengths and geometries. It is found that an inhomogeneous density distribution lowers the column densities of H2O and O2 compared to the homogeneous case by more than an order of magnitude at the same A_V. O2 is particularly sensitive to the penetrating ultraviolet radiation, more so than H2O. The S140 and rho Oph clouds are studied as relevant test cases of star-forming and quiescent regions. The SWAS results of S140 can be accommodated naturally in a clumpy model with mean density of 2x10^3 cm-3 and enhancement I_UV=140 compared with the average interstellar radiation field, in agreement with observations of [CI] and 13CO of this cloud. Additional radiative transfer computations suggest that this diffuse H2O component is warm, ~60-90 K, and can account for the bulk of the 1_10-1_01 line emission observed by SWAS. The rho Oph model yields consistent O2 abundances but too much H2O, even for [C]/[O]=0.94, if I_UV<10 respectively <40 for a mean density of 10^3 respectively 10^4 cm-3. It is concluded that enhanced photodissociation in clumpy regions can explain the low H2O and O2 abundances and emissivities found in the large SWAS beam for extended molecular clouds, but that additional freeze-out of oxygen onto grains is needed in dense cold cores.

A thought experiment is proposed to unify quantum mechanics and general relativity. The central paradigm is that space-time {\it topology} is ultimately responsible for the Heisenberg uncertaintly principle. It is found that Plankian space-time exhibits a complicated, but also definite, multiply connected character. In this framework, an analysis of the interactions in Nature is presented. I. The Universal ground state of the constructed theory derives from the properties of the topological manifold $Q=2T^3\oplus 3S^1\times S^2$, which has 23 intrinsic degrees of freedom, discrete $Z_3$ and $Z_2\times Z_3$ internal groups, an SU(5) gauge group, and leads to a U(1) symmetry on a lattice. The structure of $Q$ provides a unique equation motion for the mass-energy and particle rest mass wave functions. In its excited state the Universe is characterized by a lattice of three-tori, $L(T^3)$. The topological identifications present in this structure, a direct reflection of the Heisenberg uncertainty principle, provide the boundary conditions for solutions to the equation of motion, and suggest an interpretation for the conceptually difficult concept of quantum mechanical entanglement. II. In the second half of the paper the (observable) properties of $Q$ and $L(T^3)$ are investigated. One reproduces the standard model, and the theory naturally contains a Higgs field with possible inflation. The electron and its neutrino are identified as particle ground states and their masses, together with those of all other known particles, are predicted. A mass of $m_{\rm H}=131.6$ GeV is found for the Higgs boson. [Abridged]

Based on our high resolution, two-dimensional hydrodynamical simulations, we propose that large cavities may be formed by the nonlinear development of the combined thermal and gravitational instabilities, without need for stellar energy injection in a galaxy modeling the Large Magellanic Clouds (LMC). Our numerical model of the star formation allows us to follow the evolution of the blastwaves due to supernovae in the inhomogeous, multi-phase, and turbulent-like media self-consistently. Formation of kpc-scale inhomogeneity, such as cavities, observed HI map of the LMC, is suppressed by frequent supernovae (average supernova rate for the whole disk is ~0.001/yr). However the supernova explosions are necessary for the hot component (T_g > 10^{6-7} K). Position-velocity maps show that kpc-scale shells/arcs formed through the nonlinear evolution in a model without stellar energy feedback has similar kinematics to explosional phenomena, such as supernovae. We also find that dense clumps and filamentary structure are formed due to a natural consequence of the non-linear evolution of the multi-phage ISM. Although the ISM in a small scale looks turbulent-like and transient, the global structure of the ISM is quasi-stable. We compare the observations of HI and molecular gas of the LMC with the numerically obtained HI and CO map. The morphology and statistical properties of the numerical HI and CO maps are discussed.

The baryons inside high-redshift halos with virial temperatures above T=10^4 K cool radiatively as they condense inside dark matter potential wells. We show that the release of the gravitational binding energy, over the halo assembly time-scale, results in a significant and detectable Ly alpha flux. At the limiting line flux of approximately 10^-19 erg/s/cm^2/asec^2 of the Next Generation Space Telescope, several sufficiently massive halos, with velocity dispersions above sigma=120 km/s, would be visible per 4'x4' field. The halos would have characteristic angular sizes of approximately 10 asec, would be detectable in a broad-band survey out to redshifts z=6-8, and would provide a direct probe of galaxies in the process of forming. They may be accompanied by HeII Ly alpha emission at the 10 percent level, but remain undetectable at other wavelengths. Our predictions are in good agreement with the recent finding of two Ly alpha ``blobs'' at z=3.1 by Steidel et al. (1999).

Models are presented for the polytropic equation of state of self-gravitating, quiescent interstellar gas clouds. A detailed analysis, including chemistry, thermal balance, and radiative transfer, is performed for the physical state of the gas as a function of density, metallicity, velocity field, and background radiation field. It is found that the stiffness of the equation of state strongly depends on all these physical parameters, and the adiabatic index varies between 0.2-1.4. The implications for star formation, in particular at high redshift and in starburst galaxies, and the initial stellar mass function are discussed.

We present models for dusty high-redshift Lyman alpha emitting galaxies by combining the Press-Schechter formalism with a treatment of inhomogeneous dust distribution inside galaxies. These models reproduce the surface density of emitters inferred from recent observations, and also agree with previous non-detections. Although a detailed determination of the individual model parameters is precluded by uncertainties, we find that (i) the dust content of primordial galaxies builds up in no more than 5x10^8 years, (ii) the galactic HII regions are inhomogeneous with a cloud covering factor of order unity, and (iii) the overall star formation efficiency is at least 5 percent. Future observations should be able to detect Lyman alpha galaxies upto redshifts of about z=8. If the universe is reionized at z(reion)<8, the corresponding decline in the number of Lyman alpha emitters at z>z(reion) could prove to be a useful probe of the reionization epoch.

ISO [C II] 158 micron, [O I] 63 micron, and H_2 9 and 17 micron observations are presented of the reflection nebula Ced 201, which is a photon-dominated region illuminated by a B9.5 star with a color temperature of 10,000 K (a cool PDR). In combination with ground based [C I] 609 micron, CO, 13CO, CS and HCO+ data, the carbon budget and physical structure of the reflection nebula are constrained. The obtained data set is the first one to contain all important cooling lines of a cool PDR, and allows a comparison to be made with classical PDRs. To this effect one- and three-dimensional PDR models are presented which incorporate the physical characteristics of the source, and are aimed at understanding the dominant heating processes of the cloud. The contribution of very small grains to the photo-electric heating rate is estimated from these models and used to constrain the total abundance of PAHs and small grains. Observations of the pure rotational H_2 lines with ISO, in particular the S(3) line, indicate the presence of a small amount of very warm, approximately 330 K, molecular gas. This gas cannot be accommodated by the presented models.

We present models for the interstellar medium in disk galaxies. In particular, we investigate whether the ISM in low surface brightness galaxies can support a significant fraction of molecular gas given their low metallicity and surface density. It is found that the abundance and line brightness of CO in LSB galaxies is small and typically below current observational limits. Still, depending on physical details of the ISM, the fraction of gas in the form of molecular hydrogen can be significant in the inner few kpc of a low surface brightness galaxy. This molecular gas would be at temperatures of $\sim$ 30 -- 50 K, rather higher than in high surface brightness galaxies. These results may help explain the star forming properties and inferred evolutionary history of LSB galaxies.

Models are presented for the Lyman alpha emission of dusty high-redshift galaxies by combining the Press-Schechter formalism with a treatment of the inhomogeneous dust distribution inside galaxies. It is found that the amount of Lyman alpha radiation escaping from the galaxies strongly depends on the time over which the dust is produced through stellar activity, and on the ambient inhomogeneity of the HII regions that surround the ionizing OB stars. Good agreement is found with recent observations, as well as previous non-detections. Our models indicate that the dust content builds up in no more than approximately 5x10^8 yr, the galactic HII regions are inhomogeneous with a cloud covering factor of order unity, and the overall star formation efficiency is at least about 5%. It is predicted that future observations can detect these Lyman alpha galaxies upto redshifts of about 8.

We report an upper limit of 9 x 10^{12} cm-2 on the column density of water in the translucent cloud along the line of sight toward HD 154368. This result is based upon a search for the C-X band of water near 1240 \AA carried out using the Goddard High Resolution Spectrograph of the Hubble Space Telescope. Our observational limit on the water abundance together with detailed chemical models of translucent clouds and previous measurements of OH along the line of sight constrain the branching ratio in the dissociative recombination of H_3O+ to form water. We find at the $3\sigma$ level that no more than 30% of dissociative recombinations of H_3O+ can lead to H_2O. The observed spectrum also yielded high-resolution observations of the Mg II doublet at 1239.9 \AA and 1240.4 \AA, allowing the velocity structure of the dominant ionization state of magnesium to be studied along the line of sight. The Mg II spectrum is consistent with GHRS observations at lower spectral resolution that were obtained previously but allow an additional velocity component to be identified.

Effective star formation rates in tabular form are computed which yield a prescription for the star formation activity in model galaxies as a function of ambient density, metallicity, and stellar feedback. The effects of supernova explosions on the thermal balance of the Interstellar Medium (ISM) and the presence of a multi-phase ISM are explicitly included. The resulting grid of models can be implemented easily in N-body codes for the computation of star formation processes in merging galaxies and cosmological simulations.

Models are presented for CO rotational line emission by high redshift starburst galaxies. The influence of the cosmic microwave background on the thermal balance and the level populations of atomic and molecular species is explicitly included. Predictions are made for the observability of starburst galaxies through line and continuum emission between z=5 and z=30. It is found that the Millimeter Array could detect a starburst galaxy with ~10^5 Orion regions, corresponding to a star formation rate of about 30 Mo yr^{-1}, equally well at z=5 or z=30 due to the increasing cosmic microwave background temperature with redshift. Line emission is a potentially more powerful probe than dust continuum emission of very high redshift objects.

We have constructed three-dimensional models for the equilibrium abundance of molecular hydrogen within diffuse interstellar clouds of arbitrary geometry that are illuminated by ultraviolet radiation. The position-dependent photo- dissociation rate of H$_2$ within such clouds was computed using a 26-ray approximation to model the attenuation of the incident ultraviolet radiation field by dust and by H$_2$ line absorption. We have applied our modeling technique to the isolated diffuse cloud G236+39, assuming that the cloud has a constant density and that the thickness of the cloud along the line of sight is at every point proportional to the 100 um continuum intensity measured by IRAS. We find that our model can successfully account for observed variations in the ratio of 100 umu continuum intensity to HI column density, with larger values of that ratio occurring along lines of sight in which the molecular hydrogen fraction is expected to be largest. Using a standard chi^2 analysis to assess the goodness of fit of our models, we find (at the 60sigma level) that a three-dimensional model is more successful in matching the observational data than a one-dimensional model in which the geometrical extent of the cloud along the line of sight is assumed to be very much smaller than its extent in the plane-of-the-sky. If D is the distance to G236+39, and given standard assumptions about the rate of grain-catalysed H_2 formation, we find that the cloud has an extent along the line of sight that is 0.9+-0.1 times its mean extent projected onto the plane of the sky; a gas density of (53+-8)(100 pc/D) cm^-3; and is illuminated by a radiation field of (1.1+-0.2) (100 pc/D) times the mean interstellar radiation field estimated by Draine (1978). The derived 100 um emissivity per nucleon is (1.13+-0.06)x10^-20 MJy sr^-1 cm^2.

We present models of the cosmological star formation and metal production history of (proto-)galaxies with varying axis ratios. More massive and/or roughly spherical systems reach the threshold-metallicity for a transition to a multi-phase interstellar medium earlier than less massive, more flattened systems. Therefore, more flattened, lower-mass systems start to form stars actively at smaller redshifts. A natural explanation is found in the overall robustness of the interstellar medium against complete expulsion (blow-away) at high total masses, and in the prevention of metal enrichment in the outer regions due to axial outflow along the symmetry axis of a non-spherical proto-galaxy (blow-out). We suggest that the observed predominance of spheroidal systems observed at high redshift, e.g. in the Hubble Deep Field, is due to this effect: At z>2, roundish (proto-)galaxies with total (dark+baryonic) masses of ~10^11 M_o and/or the inner spheroidal cores of similarly massive flattened systems sustain a multi-phase interstellar medium, and therefore a high star-formation rate, whose magnitude depends on the fraction of baryonic matter in the systems. Conversely, the peak at z~1-2 in the observed cosmological metal production rate coincides with the epochs of star formation of lower mass spheroidals, as well as of massive proto-galactic disks.

A set of algebraic equations for the topological properties of space-time is derived, and used to extend general relativity into the Planck domain. A unique basis set of three-dimensional prime manifolds is constructed which consists of $S^3$, $S^1\times S^2$, and $T^3$. The action of a loop algebra on these prime manifolds yields topological invariants which constrain the dynamics of the four-dimensional space-time manifold. An extended formulation of Mach's principle and Einstein's equivalence of inertial and gravitational mass is proposed which leads to the correct classical limit of the theory. It is found that the vacuum possesses four topological degrees of freedom corresponding to a lattice of three-tori. This structure for the quantum foam naturally leads to gauge groups O(n) and SU(n) for the fields, a boundary condition for the universe, and an initial state characterized by local thermal equilibrium. The current observational estimate of the cosmological constant is reproduced without fine-tuning and found to be proportional to the number of macroscopic black holes. The black hole entropy follows immediately from the theory and the quantum corrections to its Schwarzschild horizon are computed.

We study the formation of molecular hydrogen, after the epoch of re-ionization, in the context of canonical galaxy formation theory due to hierarchical clustering. There is an initial epoch of $H_2$ production in the gas phase through the $H^-$ route which ends at a redshift of order unity. Star formation in the protogalactic disks can become self-regulated. The process responsible for the feedback is heating of the gas by the internal stellar radiation field which can dominate the background radiation field at various epochs. It is possible to define a maximum star formation rate during this epoch. Plausible estimates give a rate of 0.2-2 Mo yr-1 for condensations corresponding to 1 sigma and 2 sigma initial density fluctuations. Therefore, the production of metals and dust proceeds slowly in this phase. This moderate epoch is terminated by a phase transition to a cold dense and warm neutral/ionized medium once the metals and dust have increased to a level Z=0.03-0.1 Z_o. Then: (1) atoms and molecules such as C, O and CO become abundant and cool the gas to below $300 K$ ; (2) the dust abundance has become sufficiently high to allow shielding of the molecular gas and; (3) molecular hydrogen formation can occur rapidly on grain surfaces. This phase transition occurs at a redshift of approximately 1.5, with a fiducial range of 1.2<z<2. The combination of feedback and a phase transition provides a natural resolution of the G-dwarf problem.