We use Bayesian convolutional neural networks and a novel generative model of Galaxy Zoo volunteer responses to infer posteriors for the visual morphology of galaxies. Bayesian CNN can learn from galaxy images with uncertain labels and then, for previously unlabelled galaxies, predict the probability of each possible label. Our posteriors are well-calibrated (e.g. for predicting bars, we achieve coverage errors of 11.8% within a vote fraction deviation of 0.2) and hence are reliable for practical use. Further, using our posteriors, we apply the active learning strategy BALD to request volunteer responses for the subset of galaxies which, if labelled, would be most informative for training our network. We show that training our Bayesian CNNs using active learning requires up to 35-60% fewer labelled galaxies, depending on the morphological feature being classified. By combining human and machine intelligence, Galaxy Zoo will be able to classify surveys of any conceivable scale on a timescale of weeks, providing massive and detailed morphology catalogues to support research into galaxy evolution.

We use deep, $\mu_{r} \lesssim 28\,\text{mag}\,\text{arcsec}^{-2}$, $r$-band imaging from the Dark Energy Camera Legacy Survey (DECaLS) to search for past, or ongoing, merger activity in a sample of $282$ Low Excitation Radio Galaxies (LERGs) at $z<0.07$. Our principle aim is to assess the the role of mergers in the evolution of LERGs. Exploiting the imaging depth, we classify tidal remnants around galaxies as both minor and major morphological disturbances for our LERG sample and $1,622$ control galaxies matched in redshift, stellar mass, and environment. In groups and in the field, the LERG minor merger fraction is consistent with the control population. In galaxy clusters, $8.8 \pm 2.9\,$ % of LERGs show evidence of recent minor mergers in contrast to $23.0\pm 2.0\,$ % of controls. This $\sim 4 \sigma $ deficit of minor mergers in cluster LERGs suggests these events may inhibit this type of nuclear activity for galaxies within the cluster environment. We observe a $> 4\sigma$ excess of major mergers in the LERGs with $M_{*} \lesssim 10^{11}\,\text{M}_{\odot}$, with $10 \pm 1.5\,$ % of these AGN involved in such large-scale interactions compared to $3.2 \pm 0.4\,$% of control galaxies. This excess of major mergers in LERGs decreases with increasing stellar mass, vanishing by $M_{*} > 10^{11.3}\,\text{M}_{\odot}$. These observations show that minor mergers do not fuel LERGs, and are consistent with typical LERGs being powered by accretion of matter from their halo. Where LERGs are associated with major mergers, these objects may evolve into more efficiently accreting active galactic nuclei as the merger progresses and more gas falls on to the central engine.

LSST and Euclid must address the daunting challenge of analyzing the unprecedented volumes of imaging and spectroscopic data that these next-generation instruments will generate. A promising approach to overcoming this challenge involves rapid, automatic image processing using appropriately trained Deep Learning (DL) algorithms. However, reliable application of DL requires large, accurately labeled samples of training data. Galaxy Zoo Express (GZX) is a recent experiment that simulated using Bayesian inference to dynamically aggregate binary responses provided by citizen scientists via the Zooniverse crowd-sourcing platform in real time. The GZX approach enables collaboration between human and machine classifiers and provides rapidly generated, reliably labeled datasets, thereby enabling online training of accurate machine classifiers. We present selected results from GZX and show how the Bayesian aggregation engine it uses can be extended to efficiently provide object-localization and bounding-box annotations of two-dimensional data with quantified reliability. DL algorithms that are trained using these annotations will facilitate numerous panchromatic data modeling tasks including morphological classification and substructure detection in direct imaging, as well as decontamination and emission line identification for slitless spectroscopy. Effectively combining the speed of modern computational analyses with the human capacity to extrapolate from few examples will be critical if the potential of forthcoming large-scale surveys is to be realized.

Stars like our Sun form in self-gravitating dense and cold structures within interstellar clouds, called pre-stellar cores. Although much is known about the physical structure of dense clouds just before and soon after the switch-on of a protostar, the central few thousand astronomical units (au) of pre-stellar cores are unexplored. It is within these central regions that stellar systems assemble and fragmentation may take place, with the consequent formation of binaries and multiple systems. We present ALMA Band 6 observations (ACA and 12m array) of the dust continuum emission of the 8 Msun pre-stellar core L1544, with angular resolution of 2'' x 1.6'' (linear resolution 270 au x 216 au). Within the primary beam, a compact region of 0.1 Msun, which we call a "kernel", has been unveiled. The kernel is elongated, with a central flat zone with radius Rker ~ 10'' (~ 1400 au). The average number density within Rker is ~1 x 10^6 cm^{-3}, with possible local density enhancements. The region within Rker appears to have fragmented, but detailed analysis shows that similar substructure can be reproduced by synthetic interferometric observations of a smooth centrally concentrated dense core with a similar central flat zone. The presence of a smooth kernel within a dense core is in agreement with non-ideal magneto-hydro-dynamical simulations of a contracting cloud core with a peak number density of 1 x 10^7 cm^{-3}. Dense cores with lower central densities are completely filtered out when simulated 12m-array observations are carried out. These observations demonstrate that the kernel of dynamically evolved dense cores can be investigated at high angular resolution with ALMA.

Faint tidal features around galaxies record their merger and interaction histories over cosmic time. Due to their low surface brightnesses and complex morphologies, existing automated methods struggle to detect such features and most work to date has heavily relied on visual inspection. This presents a major obstacle to quantitative study of tidal debris features in large statistical samples, and hence the ability to be able to use these features to advance understanding of the galaxy population as a whole. This paper uses convolutional neural networks (CNNs) with dropout and augmentation to identify galaxies in the CFHTLS-Wide Survey that have faint tidal features. Evaluating the performance of the CNNs against previously-published expert visual classifications, we find that our method achieves high (76%) completeness and low (20%) contamination, and also performs considerably better than other automated methods recently applied in the literature. We argue that CNNs offer a promising approach to effective automatic identification of low surface brightness tidal debris features in and around galaxies. When applied to forthcoming deep wide-field imaging surveys (e.g. LSST, Euclid), CNNs have the potential to provide a several order-of-magnitude increase in the sample size of morphologically-perturbed galaxies and thereby facilitate a much-anticipated revolution in terms of quantitative low surface brightness science.

The processes leading to the birth of high-mass stars are poorly understood. We characterise here a sample of 430 massive clumps from the ATLASGAL survey, which are representative of different evolutionary stages. To establish a census of molecular tracers of their evolution we performed an unbiased spectral line survey covering the 3-mm atmospheric window between 84-117 GHz with the IRAM 30m. A smaller sample of 128 clumps has been observed in the SiO (5-4) transition with the APEX telescope to complement the SiO (2-1) line and probe the excitation conditions of the emitting gas, which is the main focus of the current study. We report a high detection rate of >75% of the SiO (2-1) line and a >90% detection rate from the dedicated follow-ups in the (5-4) transition. The SiO (2-1) line with broad line profiles and high detection rates, is a powerful probe of star formation activity, while the ubiquitous detection of SiO in all evolutionary stages suggests a continuous star formation process in massive clumps. We find a large fraction of infrared-quiet clumps to exhibit SiO emission, the majority of them only showing a low-velocity component (FWHM~5-6 km/s) centred at the rest velocity of the clump. In the current picture, where this is attributed to low-velocity shocks from cloud-cloud collisions, this can be used to pinpoint the youngest, thus, likely prestellar massive structures. Based on the line ratio of the (5-4) to the (2-1) line, our study reveals a trend of changing excitation conditions that lead to brighter emission in the (5-4) line towards more evolved sources. Our analysis delivers a more robust estimate of SiO column density and abundance than previous studies and questions the decrease of jet activity in massive clumps as a function of age.

The formation processes and the evolutionary stages of high-mass stars are poorly understood compared to low-mass stars. Large-scale surveys are needed to provide an unbiased census of high column density sites which can potentially host precursors to high-mass stars. Here we use the ATLASGAL survey covering 420 sq. degree of the Galactic plane at 870 $\mu$m; and use the MRE-GLC method to identify the population of embedded sources throughout the inner Galaxy. We identify in total 10861 compact sub-millimeter sources with fluxes above 5 sigma. Completeness tests show that our catalogue is 97% complete above 5 sigma and >99% complete above 7$\sigma$. We correlate this sample with mid-infrared point source catalogues (MSX at 21.3 $\mu$m and WISE at 22 $\mu$m) and determine a lower limit of ~33% that are associated with embedded protostellar objects. We note that the proportion of clumps associated with mid-infrared sources increases with increasing flux density, achieving a rather constant fraction of ~75% of all clumps with fluxes over 5 Jy/beam being associated with star-formation. Examining the source counts as a function of Galactic longitude we are able to identify the most prominent star forming regions in the Galaxy. From the fraction of the likely massive quiescent clumps (~25%) we estimate a formation time-scale of ~7.5+/-2.5 $\times$ 10$^4$yr for the deeply embedded phase before the emergence of luminous YSOs. Such a short duration for the formation of high-mass stars in massive clumps clearly proves that the earliest phases have to be dynamic with supersonic motions.

To understand the origin of water line emission and absorption during high-mass star formation, we decompose high-resolution Herschel-HIFI line spectra toward 19 high-mass star-forming regions into three distinct physical components. Protostellar envelopes are usually seen as narrow absorptions or emissions in the H2O 1113 and 1669 GHz ground-state lines, the H2O 987 GHz excited-state line, and the H2O-18 1102 GHz ground-state line. Broader features due to outflows are usually seen in absorption in the H2O 1113 and 1669 GHz lines, in 987 GHz emission, and not seen in H2O-18, indicating a low column density and a high excitation temperature. The H2O 1113 and 1669 GHz spectra show narrow absorptions by foreground clouds along the line of sight, which have a low column density and a low excitation temperature, although their H2O ortho/para ratios are close to 3. The intensities of the H2O 1113 and 1669 GHz lines do not show significant trends with luminosity, mass, or age. In contrast, the 987 GHz line flux increases with luminosity and the H2O-18 line flux decreases with mass. Furthermore, appearance of the envelope in absorption in the 987 GHz and H2O-18 lines seems to be a sign of an early evolutionary stage. We conclude that the ground state transitions of H2O trace the outer parts of the envelopes, so that the effects of star formation are mostly noticeable in the outflow wings. These lines are heavily affected by absorption, so that line ratios of H2O involving the ground states must be treated with caution. The average H2O abundance in high-mass protostellar envelopes does not change much with time. The 987 GHz line appears to be a good tracer of the mean weighted dust temperature of the source, which may explain why it is readily seen in distant galaxies.

G18.93-0.03 is a prominent dust complex within an 0.8deg long filament, with the molecular clump G18.93/m being IR dark from near IR wavelength up to 160mu. Spitzer composite images show an IR bubble spatially associated with G18.93. We use GRS 13CO and IRAM 30m H13CO+ data to disentangle the spatial structure of the region. From ATLASGAL submm data we calculate the gas mass, while we use the H13CO+ line width to estimate its virial mass. Using HERSCHEL data we produce temperature maps from fitting the SED. With the MAGPIS 20cm and SuperCOSMOS Halpha data we trace the ionized gas, and the VGPS HI survey provides information on the atomic hydrogen gas. We show that the bubble is spatially associated with G18.93, located at a kinematic near distance of 3.6kpc. With 280Msun, the most massive clump within G18.93 is G18.93/m. The virial analysis shows that it may be gravitationally bound and has neither Spitzer young stellar objects nor mid-IR point sources within. Fitting the SED reveals a temperature distribution that decreases towards its center, but heating from the ionizing source puts it above the general ISM temperature. We find that the bubble is filled by HII gas, ionized by an O8.5 star. Between the ionizing source and the IR dark clump G18.93/m we find a layered structure, from ionized to atomic to molecular hydrogen, revealing a PDR. Furthermore, we identify an additional velocity component within the bubble's 8mu emission rim at the edge of the infrared dark cloud and speculate that it might be shock induced by the expanding HII region. While the elevated temperature allows for the build-up of larger fragments, and the shock induced velocity component may lead to additional turbulent support, we do not find conclusive evidence that the massive clump G18.93/m is prone to collapse because of the expanding HII region.

We present one of the first resolved maps of the [CII] 158 micron line, a powerful tracer of the star forming inter-stellar medium, at high redshift. We use the new IRAM PdBI receivers at 350 GHz to map this line in BRI 0952-0115, the host galaxy of a lensed quasar at z=4.4 previously found to be very bright in [CII] emission. The [CII] emission is clearly resolved and our data allow us to resolve two [CII] lensed images associated with the optical quasar images. We find that the star formation, as traced by [CII], is distributed over a region of ~ 1 kpc in size near the quasar nucleus, and we infer a star formation surface density >150 Msun/yr/kpc^2, similar to that observed in local ULIRGs. We also reveal another [CII] component, extended over ~ 12 kpc, and located at ~ 10 kpc from the quasar. We suggest that this component is a companion disk galaxy, in the process of merging with the quasar host, whose rotation field is distorted by the interaction with the quasar host, and where star formation, although intense, is more diffuse. These observations suggest that galaxy merging at high-z can enhance star formation at the same time in the form of more compact regions, in the vicinity of the accreting black hole, and in more extended star forming galaxies.

We investigated the chemistry of nitrogen--containing species, principally isotopomers of CN, HCN, and HNC, in a sample of pre-protostellar cores. We used the IRAM 30 m telescope to measure the emission in rotational and hyperfine transitions of CN, HCN, 13CN, H13CN, HN13C, and HC15N, in L 1544, L 183, Oph D, L 1517B, L 310. The observations were made along axial cuts through the dust emission peak, at a number of regularly--spaced offset positions. The observations were reduced and analyzed to obtain the column densities, using the measurements of the less abundant isotopic variants in order to minimize the consequences of finite optical depths in the lines. The observations were compared with the predictions of a free--fall gravitational collapse model, which incorporates a non-equilibrium treatment of the relevant chemistry. We found that CN, HCN, and HNC remain present in the gas phase at densities well above that at which CO depletes on to grains. The CN:HCN and the HNC:HCN abundance ratios are larger than unity in all the objects of our sample. Furthermore, there is no observational evidence for large variations of these ratios with increasing offset from the dust emission peak and hence with density. Whilst the differential freeze--out of CN and CO can be understood in terms of the current chemistry, the behaviour of the CN:HCN ratio is more difficult to explain. Models suggest that most nitrogen is not in the gas phase but may be locked in ices. Unambiguous conclusions require measurements of the rate coefficients of the key neutral--neutral reactions at low temperatures.

We report the detection of the [CII]157.74um fine-structure line in the lensed galaxy BRI 0952-0115 at z=4.43, using the APEX telescope. This is the first detection of the [CII] line in a source with L_FIR < 10^13 L_sun at high redshift. The line is very strong compared to previous [CII] detections at high-z (a factor of 5-8 higher in flux), partly due to the lensing amplification. The L_[CII]/L_FIR ratio is 10^-2.9, which is higher than observed in local galaxies with similar infrared luminosities. Together with previous observations of [CII] at high redshift, our result suggests that the [CII] emission in high redshift galaxies is enhanced relative to local galaxies of the same infrared luminosity. This finding may result from selection effects of the few current observations of [CII] at high redshift, and in particular the fact that non detections may have not been published (although the few published upper limits are still consistent with the [CII] enhancement scenario). If the trend is confirmed with larger samples, it would indicate that high-z galaxies are characterized by different physical conditions with respect to their local counterparts. Regardless of the physical origin of the trend, this effect would increase the potential of the [CII]158um line to search and characterize high-z sources.

Determining the structure of and the velocity field in prestellar cores is essential to understanding protostellar evolution.} {We have observed the dense prestellar cores L 1544 and L 183 in the $N = 1 \to 0$ rotational transition of CN and \thcn in order to test whether CN is depleted in the high--density nuclei of these cores.} {We have used the IRAM 30 m telescope to observe along the major and minor axes of these cores. We compare these observations with the 1 mm dust emission, which serves as a proxy for the hydrogen column density.}{We find that while CN\jone is optically thick, the distribution of \thcn\jone intensity follows the dust emission well, implying that the CN abundance does not vary greatly with density. We derive an abundance ratio of $\rm [CN]/[\hh]=\dix{-9}$ in L 183 and 1-3\tdix{-9} in L 1544, which, in the case of L 183, is similar to previous estimates obtained by sampling lower--density regions of the core.}{We conclude that CN is not depleted towards the high--density peaks of these cores and thus behaves like the N-containing molecules \nnhp and \nhhh. CN is, to our knowledge, the first C--containing molecule to exhibit this characteristic.

Context: The evolution of a young protocluster depends on the relative spatial distributions and dynamics of both stars and gas. Aims: We study the distribution and properties of the gas and stars surrounding the luminous (10^4 L_sun) protocluster IRAS 18511+0146. Methods: IRAS 18511+0146 and the cluster associated with it has been investigated using the sub-millimetre (JCMT-SCUBA), infrared (Spitzer-MIPSGAL, Spitzer-GLIMPSE, Palomar) and radio (VLA) continuum data. Cluster simulations have been carried out in order to understand the properties of clusters as well as to compare with the observations. Results: The central most obscured part of the protocluster coincident with the compact sub-millimetre source found with SCUBA is responsible for at least 2/3 of the total luminosity. A number of cluster members have been identified which are bright in mid infrared and show rising (near to mid infrared) spectral energy distributions suggesting that these are very young stellar sources. In the mid infrared 8.0 micron image, a number of filamentary structures and clumps are detected in the vicinity of IRAS 18511+0146. Conclusions: Based on the luminosity and cluster size as well as on the evolutionary stages of the cluster members, IRAS 18511+0146 is likely to be protocluster with the most massive object being a precursor to a Herbig type star.

Abridged: The thermal structure of a starless core is crucial for our understanding of the physics in these objects and hence for our understanding of star formation. Theory predicts a gas temperature drop in the inner 5000 AU of these objects, but there has been little observational proof of this. We performed VLA observations of the NH3 (1,1) and (2,2) transitions towards the pre-stellar core L 1544 in order to measure the temperature gradient between the high density core nucleus and the surrounding core envelope. Our VLA observation for the first time provide measurements of gas temperature in a core with a resolution smaller than 1000 AU. We have also obtained high resolution Plateau de Bure observations of the 110 GHz 111-101 para-NH2D line in order to further constrain the physical parameters of the high density nucleus. We have estimated the temperature gradient using a model of the source to fit our data in the u,v plane. We find that indeed the temperature decreases toward the core nucleus from 12 K down to 5.5 K resulting in an increase of a factor of 50% in the estimated density of the core from the dust continuum if compared with the estimates done with constant temperature of 8.75 K. We also found a remarkably high abundance of deuterated ammonia with respect to the ammonia abundance (50%+-20%), which proves the persistence of nitrogen bearing molecules at very high densities (2e6 cm-3) and shows that high-resolution observations yield higher deuteration values than single-dish observations. Our analysis of the NH3 and NH2D kinematic fields shows a decrease of specific angular momentum from the large scales to the small scales.

Understanding the depletion of heavy elements is a fundamental step towards determining the structure of pre-protostellar cores just prior to collapse. We study the dependence of the NO abundance on position in the pre-protostellar cores L1544 and L183. We observed the 150 GHz and 250~GHz transitions of NO and the 93 GHz transitions of \NTHP \ towards L1544 and L183 using the IRAM 30 m telescope. We compare the variation of the NO column density with position in these objects with the H column density derived from dust emission measurements. We find that NO behaves differently from \NTHP \ and appears to be partially depleted in the high density core of L1544. Other oxygen-containing compounds are also likely to be partially depleted in dense-core nuclei. The principal conclusions are that: the prestellar core L1544 is likely to be 'carbon-rich'; the nitrogen chemistry did not reach equilibrium prior to gravitational collapse, and nitrogen is initially (at densities of the order of $10^4$~cm$^{-3}$) mainly in atomic form; the grain sticking probabilities of atomic C, N and, probably, O are significantly smaller than unity.

Single-dish spectra and interferometric maps of (sub)mm lines of H2O-18 and HDO are used to study the chemistry of water in eight regions of high-mass star formation. The spectra indicate HDO excitation temperatures of ~110 K and column densities in an 11'' beam of ~2x10^14 cm^-2 for HDO and ~2x10^17 cm^-2 for H2O, with the N(HDO)/N(H2O) ratio increasing with decreasing temperature. Simultaneous observations of CH3OH and SO2 indicate that 20-50% of the single-dish line flux arises in the molecular outflows of these objects. The outflow contribution to the H2O-18 and HDO emission is estimated to be 10-20%. Radiative transfer models indicate that the water abundance is low (~10^-6) outside a critical radius corresponding to a temperature in the protostellar envelope of ~100 K, and `jumps' to H2O/H2 ~ 10^-4 inside this radius. This value corresponds to the observed abundance of solid water and together with the derived HDO/H2O abundance ratios of ~1/1000 suggests that the origin of the observed water is evaporation of grain mantles. This idea is confirmed in the case of AFGL 2591 by interferometer observations of HDO, H2O-18 and SO2 lines, which reveal compact (~800 AU) emission with a systematic velocity gradient. This size is similar to that of the 1.3 mm continuum towards AFGL 2591, from which we estimate a mass of ~0.8 M0, or ~5% of the mass of the central star. We speculate that we may be observing a circumstellar disk in an almost face-on orientation.

We reconsider H2 formation on grain surfaces. We develop a rate equation model which takes into account the presence of both physisorbed and chemisorbed sites on the surface, including quantum mechanical tunnelling and thermal diffusion. In this study, we took into consideration the uncertainties on the characteristics of graphitic surfaces. We calculate the H2 formation efficiency with the Langmuir Hinshelwood and Eley Rideal mechanisms, and discuss the importance of these mechanisms for a wide range of grain and gas temperatures. We also develop a Monte Carlo simulation to calculate the H2 formation efficiency and compare the results to our rate equation models. Our results are the following: (1) Depending on the barrier against chemisorption, we predict the efficiency of H2 formation for a wide range of grain and gas temperatures. (2) The Eley-Rideal mechanism has an impact on the H2 formation efficiency at high grain and gas temperatures. (3) The fact that we consider chemisorption in our model makes the rate equation and Monte Carlo approaches equivalent.

We compute polarization maps for molecular cloud cores modeled as magnetized singular isothermal toroids, under the assumption that the emitting dust grains are aspherical and aligned with the large-scale magnetic field. We show that, depending on the inclination of the toroid with the line-of-sight, the bending of the magnetic field lines resulting from the need to counteract the inward pull of gravity naturally produces a depolarization effect toward the centre of the map. We compute the decrease of polarization degree with increasing intensity for different viewing angles and frequencies, and we show that an outward increasing temperature gradient, as expected in starless cores heated by the external radiation field, enhances the decrease of polarization. We compare our results with recent observations, and we conclude that this geometrical effect, together with other mechanisms of depolarization, may significantly contribute to the decrease of polarization degrees with intensity observed in the majority of molecular cloud cores. Finally, we consider the dependence of the polarization degree on the dust temperature gradient predicted for externally heated clouds, and we briefly comment on the limits of the Chandrasekhar-Fermi formula to estimate the magnetic field strength in molecular cloud cores.

We present a survey of the polarization properties of mm-wavelength methanol masers, comprising both classes, and transitions from 84.5 to 157.0 GHz. Linear polarization is found in more than half of the observed objects, and circular polarization is tentatively detected in two sources. Class I and Class II CH_3OH masers show similar polarization properties.The largest linear polarization is found in the 133 GHz Class I maser towards L379 (39.5 %), and in the 157 GHz Class II maser towards G9.62+0.19 (36.7 %). The spectral profiles of the polarization angle of Class I masers are mostly flat, except for two sources showing a linear slope. Since the mm-line methanol masers are expected to be weakly (or not) saturated, we suggest that the stronger fractional polarizations found by us are enhanced by anisotropic pumping and radiative losses. In NGC 7538, we find, for both maser classes, a good agreement between our polarization angles, and those measured for the submillimeter dust continuum. This can be taken as evidence for magnetic alignment of dust grains. It is also possible that an unsaturated maser with equally populated magnetic substates simply amplifies polarized continuum seed radiation. For Class II masers, the polarization properties of the various velocity components towards a given source with detectable polarization are quite homogeneous. A possible explanation is discussed. Since methanol is non-paramagnetic, the circular polarization of the unsaturated maser emission can only be due to variations of the angle between the magnetic field and the line of sight along the maser propagation path.

We present the results of a three-dimensional Monte Carlo radiative transfer code for starless molecular cloud cores heated by an external isotropic or non-isotropic interstellar radiation field. The code computes the dust temperature distribution inside model clouds with specified but arbitrary density profiles. In particular we examine in detail spherical (Bonnor-Ebert) clouds, axisymmetric and non-axisymmetric toroids, and clouds heated by an external stellar source in addition to the general interstellar field. For these configurations, the code also computes maps of the emergent intensity at different wavelengths and arbitrary viewing angle, that can be compared directly with continuum maps of prestellar cores. In the approximation where the dust temperature is independent of interactions with the gas and where the gas is heated both by collisions with dust grains and ionization by cosmic rays, the temperature distribution of the gas is also calculated. For cloud models with parameters typical of dense cores, the results show that the dust temperature decreases monotonically from a maximum value near the cloud's edge (14-15 K) to a minimum value at the cloud's center (6-7 K). Conversely, the gas temperature varies in a similar range, but, due to efficient dust-gas coupling in the inner regions and inefficient cosmic-ray heating in the outer regions, the gradient is non-monotonic and the gas temperature reaches a maximum value at intermediate radii. The emission computed for these models (at 350 micron and 1.3 mm) shows that deviations from spherical symmetry in the density and/or temperature distributions are generally reduced in the simulated intensity maps (even without beam convolution), especially at the longer wavelengths.

This paper explores the sensitivity of line profiles of H2D+, HCO+ and N2H+, observed towards the center of L1544, to various kinematic and chemical parameters. The total width of the H2D+ line can be matched by a static model and by models invoking ambipolar diffusion and gravitational collapse. The derived turbulent line width is b=0.15 km/s for the static case and <~ 0.05 km/s for the collapse case. However, line profiles of HC18O+ and N2H+ rule out the static solution. The double-peaked H2D+ line shape requires either infall speeds in the center that are much higher than predicted by ambipolar diffusion models, or a shell-type distribution of H2D+, as is the case for HCO+ and N2H+. At an offset of ~20 arcsec from the dust peak, the H2D+ abundance drops by a factor of ~5.

We have observed the dark cloud Barnard 68 with SCUBA at 850 um and with SIMBA at 1.2 mm. The submillimetre and millimetre dust emission correlate well with the extinction map of Alves, Lada and Lada (2001).The A_V/850um correlation is clearly not linear and suggests lower temperatures for the dust in the inner core of the cloud. Assuming a model for the temperature gradient, we derive the cloud-averaged dust emissivities (normalised to the V-Band extinction efficiency) at 850 um and 1.2 mm. We find k_850um/k_V = 4.0 +/- 1.0 x 10^-5 and k_1.2mm/k_V = 9.0 +/- 3.0 x 10^-6. These values are compared with other determinations in this wavelength regime and with expectations for models of diffuse dust and grain growth in dense clouds.

We have considered the thermal equilibrium in pre-protostellar cores in the approximation where the dust temperature is independent of interactions with the gas and where the gas is heated both by collisions with dust grains and ionization by cosmic rays. We have then used these results to study the stability of cores in the limit where thermal pressure dominates over magnetic field and turbulence. We find that for cores with characteristics similar to those observed, the gas and dust temperatures are coupled in the core interior. As a consequence, the gas temperature like the dust temperature decreases towards the center of these objects. The density structure computed taking into account such deviations from isothermality are not greatly different from that expected for an isothermal Bonnor-Ebert sphere. It is impossible in the framework of these models to have a stable equilibrium core with mass above about 5 solar masses and column density compatible with observed values. We conclude from this that observed high mass cores are either supported by magnetic field or turbulence or are already in a state of collapse. Lower mass cores on the other hand have stable states and we conclude that the much studied object B68 may be in a state of stable equilibrium if the internal gas temperature is computed in self-consistent fashion. Finally we note that in molecular clouds such as Ophiuchus and Orion with high radiation fields and pressures, gas and dust temperatures are expected to be well coupled and hence one expects temperatures to be relatively high as compared to low pressure clouds like Taurus.