We present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other.

We analyze results from the first eighteen months of monthly sub-mm monitoring of eight star-forming regions in the JCMT Transient Survey. In our search for stochastic variability in 1643 bright peaks, only the previously identified source, EC53, shows behavior well above the expected measurement uncertainty. Another four sources, two disks and two protostars, show moderately-enhanced standard deviations in brightness, as expected for stochastic variables. For the two protostars, this apparent variability is the result of single epochs that are much brighter than the mean. In our search for secular brightness variations that are linear in time, we measure the fractional brightness change per year for 150 bright peaks, fifty of which are protostellar. The ensemble distribution of slopes is well fit by a normal distribution with sigma ~ 0.023. Most sources are not rapidly brightening or fading in the sub-mm. Comparison against time-randomized realizations shows that the width of the distribution is dominated by the uncertainty in the individual brightness measurements of the sources. A toy model for secular variability reveals that an underlying Gaussian distribution of linear fractional brightness change sigma = 0.005 would be unobservable in the present sample, whereas an underlying distribution with sigma = 0.02 is ruled out. Five protostellar sources, 10% of the protostellar sample, are found to have robust secular measures deviating from a constant flux. The sensitivity to secular brightness variations will improve significantly with a larger time sample, with a factor of two improvement expected by the conclusion of our 36-month survey.

We present the molecular cloud properties of N55 in the Large Magellanic Cloud using $^{12}$CO(1-0) and $^{13}$CO(1-0) observations obtained with Atacama Large Millimeter Array. We have done a detailed study of molecular gas properties, to understand how the cloud properties of N55 differ from Galactic clouds. Most CO emission appears clumpy in N55, and molecular cores that have YSOs show larger linewidths and masses. The massive clumps are associated with high and intermediate mass YSOs. The clump masses are determined by local thermodynamic equilibrium and virial analysis of the $^{12}$CO and $^{13}$CO emissions. These mass estimates lead to the conclusion that, (a) the clumps are in self-gravitational virial equilibrium, and (b) the $^{12}$CO(1-0)-to-H$_2$ conversion factor, X$_{\rm CO}$, is 6.5$\times$10$^{20}$cm$^{-2}$(K km s$^{-1}$)$^{-1}$. This CO-to-H$_2$ conversion factor for N55 clumps is measured at a spatial scale of $\sim$0.67 pc, which is about two times higher than the X$_{\rm CO}$ value of Orion cloud at a similar spatial scale. The core mass function of N55 clearly show a turnover below 200M$_{\odot}$, separating the low-mass end from the high-mass end. The low-mass end of the $^{12}$CO mass spectrum is fitted with a power law of index 0.5$\pm$0.1, while for $^{13}$CO it is fitted with a power law index 0.6$\pm$0.2. In the high-mass end, the core mass spectrum is fitted with a power index of 2.0$\pm$0.3 for $^{12}$CO, and with 2.5$\pm$0.4 for $^{13}$CO. This power-law behavior of the core mass function in N55 is consistent with many Galactic clouds.

The early evolutionary stage of brown dwarfs are not very well characterized, specially during the embedded phase. To gain insight into the dominant formation mechanism of very low-mass objects and brown dwarfs, we conducted deep observations at 870$\mu$m with the LABOCA bolometer at the APEX telescope. Our goal was to identify young sub-mm sources in the Barnard 30 dark cloud. We complemented these data with multi-wavelength observations from the optical to the far-IR and. As a result, we have identified 34 submm sources and a substantial number of possible and probable Barnard 30 members within each individual APEX/LABOCA beam. They can be classified in three distinct groups. First, 15 out of these 34 have a clear optical or IR counterpart to the submm peak and nine of them are potential proto-BDs candidates. Moreover, a substantial number of them could be multiple systems. A second group of 13 sources comprises candidate members with significant infrared excesses located away from the central submm emission. All of them include brown dwarf candidates, some displaying IR excess, but their association with submm emission is unclear. In addition, we have found six starless cores and, based on the total dust mass estimate, three might be pre-substellar (or pre-BDs) cores. Finally, the complete characterization of our APEX/LABOCA sources, focusing on those detected at 24 and/or 70 $\mu$m, indicates that in our sample of 34 submm sources there are, at least: two WTTs, four CTTs, five young stellar objects (YSOs), eight proto-BD candidates (with another three dubious cases), and one Very Low Luminosity object (VeLLO).

In this work we present ALMA continuum observations at 880 $\mu$m of 30 sub-mm cores previously identified with APEX/LABOCA at 870$\mu$m in the Barnard 30 cloud. The main goal is to characterize the youngest and lowest mass population in the cloud. As a result, we report the detection of five (out of 30) spatially unresolved sources with ALMA, with estimated masses between 0.9 and 67 M$_{\rm Jup}$. From these five sources, only two show gas emission. The analysis of multi-wavelength photometry from these two objects, namely B30-LB14 and B30-LB19, is consistent with one Class II- and one Class I low-mass stellar object, respectively. The gas emission is consistent with a rotating disk in the case of B30-LB14, and with an oblate rotating envelope with infall signatures in the case of LB19. The remaining three ALMA detections do not have infrared counterparts and can be classified as either deeply embedded objects or as starless cores if B30 members. In the former case, two of them (LB08 and LB31) show internal luminosity upper limits consistent with Very Low Luminosity objects, while we do not have enough information for LB10. In the starless core scenario, and taking into account the estimated masses from ALMA and the APEX/LABOCA cores, we estimate final masses for the central objects in the substellar domain, so they could be classified as pre-BD core candidates.

Investigating variability at the earliest stages of low-mass star formation is fundamental in understanding how a protostar assembles mass. While many simulations of protostellar disks predict non-steady accretion onto protostars, deeper investigation requires robust observational constraints on the frequency and amplitude of variability events characterised across the observable SED. In this study, we develop methods to robustly analyse repeated observations of an area of the sky for submillimetre variability in order to determine constraints on the magnitude and frequency of deeply embedded protostars. We compare \mbox{850 $\mu$m} JCMT Transient Survey data with archival JCMT Gould Belt Survey data to investigate variability over 2-4 year timescales. Out of 175 bright, independent emission sources identified in the overlapping fields, we find 7 variable candidates, 5 of which we classify as \textit{Strong} and the remaining 2 as \textit{Extended} to indicate the latter are associated with larger-scale structure. For the \textit{Strong} variable candidates, we find an average fractional peak brightness change per year of |4.0|\% yr$^{-1}$ with a standard deviation of $2.7\%\mathrm{\:yr}^{-1}$. In total, 7\% of the protostars associated with \mbox{850 $\mu$m} emission in our sample show signs of variability. Four of the five \textit{Strong} sources are associated with a known protostar. The remaining source is a good follow-up target for an object that is anticipated to contain an enshrouded, deeply embedded protostar. In addition, we estimate the \mbox{850 $\mu$m} periodicity of the submillimetre variable source, EC 53, to be \mbox{567 $\pm$ 32 days} based on the archival Gould Belt Survey data.

Most protostars have luminosities that are fainter than expected from steady accretion over the protostellar lifetime. The solution to this problem may lie in episodic mass accretion -- prolonged periods of very low accretion punctuated by short bursts of rapid accretion. However, the timescale and amplitude for variability at the protostellar phase is almost entirely unconstrained. In "A JCMT/SCUBA-2 Transient Survey of Protostars in Nearby Star Forming Regions", we are monitoring monthly with SCUBA-2 the sub-mm emission in eight fields within nearby (<500 pc) star forming regions to measure the accretion variability of protostars. The total survey area of ~1.6 sq.deg. includes ~105 peaks with peaks brighter than 0.5 Jy/beam (43 associated with embedded protostars or disks) and 237 peaks of 0.125-0.5 Jy/beam (50 with embedded protostars or disks). Each field has enough bright peaks for flux calibration relative to other peaks in the same field, which improves upon the nominal flux calibration uncertainties of sub-mm observations to reach a precision of ~2-3% rms, and also provides quantified confidence in any measured variability. The timescales and amplitudes of any sub-mm variation will then be converted into variations in accretion rate and subsequently used to infer the physical causes of the variability. This survey is the first dedicated survey for sub-mm variability and complements other transient surveys at optical and near-IR wavelengths, which are not sensitive to accretion variability of deeply embedded protostars.

The Atacama Large Millimeter/submillimeter Array(ALMA) Band 1 receiver covers the 35-50 GHz frequency band. Development of prototype receivers, including the key components and subsystems has been completed and two sets of prototype receivers were fully tested. We will provide an overview of the ALMA Band 1 science goals, and its requirements and design for use on the ALMA. The receiver development status will also be discussed and the infrastructure, integration, evaluation of fully-assembled band 1 receiver system will be covered. Finally, a discussion of the technical and management challenges encountered will be presented.

Pre-stellar cores within molecular clouds provide the very initial conditions in which stars are formed. We use the IRAM 30m telescope and the PdBI to study the chemical and physical properties of the starless core FeSt 1-457 (Core 109), in the Pipe nebula. We fit the hyperfine structure of the N2H+(1-0) IRAM 30m data. This allow us to measure with high precision the velocity field, line widths and opacity and derive the excitation temperature and column density in the core. We use a modified Bonnor-Ebert sphere model adding a temperature gradient towards the center to fit the 1.2 mm continuum emission and visual extinction maps. Using this model, we estimate the abundances of the N2H+ and the rest of molecular lines detected in the 30 GHz wide line survey performed at 3 mm with IRAM 30m using ARTIST software. The core presents a rich chemistry with emission from early (C3H2, HCN, CS) and late-time molecules (e.g., N2H+), with a clear chemical spatial differentiation for nitrogen, oxygen and sulphurated molecules. For most of the molecules detected (HCN, HCO+, CH3OH, CS, SO, 13CO and C18O), abundances are best fitted with three values, presenting a clear decrease of abundance of at least 1 or 2 orders of magnitude towards the center of the core. The Bonnor-Ebert analysis indicates the core is gravitationally unstable and the magnetic field is not strong enough to avoid the collapse. Depletion of molecules onto the dust grains occurs at the interior of the core, where dust grain growth and dust depolarization also occurs. This suggests that these properties may be related. On the other hand, some molecules exhibit asymmetries in their integrated emission maps, which appear to be correlated with a previously reported submillimetre polarization asymmetry. These asymmetries could be due to a stronger interstellar radiation field in the western side of the core.

We present millimetre (SMA) and sub-millimetre (SCUBA-2) continuum observations of the peculiar star KIC 8462852 which displayed several deep and aperiodic dips in brightness during the Kepler mission. Our observations are approximately confusion-limited at 850 $\mu$m and are the deepest millimetre and sub-millimetre photometry of the star that has yet been carried out. No significant emission is detected towards KIC 8462852. We determine upper limits for dust between a few 10$^{-6}$ M$_{\oplus}$ and 10$^{-3}$ M$_{\oplus}$ for regions identified as the most likely to host occluding dust clumps and a total overall dust budget of $<$7.7 M$_{\oplus}$ within a radius of 200 AU. Such low limits for the inner system make the catastrophic planetary disruption hypothesis unlikely. Integrating over the Kepler lightcurve we determine that at least 10$^{-9}$ M$_{\oplus}$ of dust is required to cause the observed Q16 dip. This is consistent with the currently most favoured cometary breakup hypothesis, but nevertheless implies the complete breakup of $\sim$ 30 Comet 1/P Halley type objects. Finally, in the wide SCUBA-2 field-of-view we identify another candidate debris disc system that is potentially the largest yet discovered.

Context. Chamaeleon II molecular cloud is an active star forming region that offers an excellent opportunity for studying the formation of brown dwarfs in the southern hemisphere. Aims. Our aims are to identify a population of pre- and proto- brown dwarfs (5 sigma mass limit threshold of ~0.015 Msun) and provide information on the formation mechanisms of substellar objects. Methods. We performed high sensitivity observations at 870 microns using the LABOCA bolometer at the APEX telescope towards an active star forming region in Chamaeleon II. The data are complemented with an extensive multiwavelength catalogue of sources from the optical to the far-infrared to study the nature of the LABOCA detections. Results. We detect fifteen cores at 870 microns, and eleven of them show masses in the substellar regime. The most intense objects in the surveyed field correspond to the submillimeter counterparts of the well known young stellar objects DK Cha and IRAS 12500-7658. We identify a possible proto-brown dwarf candidate (ChaII-APEX-L) with IRAC emission at 3.6 and 4.5 microns. Conclusions. Our analysis indicates that most of the spatially resolved cores are transient, and that the point-like starless cores in the sub-stellar regime (with masses between 0.016 Msun and 0.066 Msun) could be pre-brown dwarfs cores gravitationally unstable if they have radii smaller than 220 AU to 907 AU (1.2" to 5" at 178 pc) respectively for different masses. ALMA observations will be the key to reveal the energetic state of these pre-brown dwarfs candidates.

Ground-based observations at terahertz (THz) frequencies are a newly explorable area of astronomy for the next ten years. We discuss science cases for a first-generation 10-m class THz telescope, focusing on the Greenland Telescope as an example of such a facility. We propose science cases and provide quantitative estimates for each case. The largest advantage of ground-based THz telescopes is their higher angular resolution (~ 4 arcsec for a 10-m dish), as compared to space or airborne THz telescopes. Thus, high-resolution mapping is an important scientific argument. In particular, we can isolate zones of interest for Galactic and extragalactic star-forming regions. The THz windows are suitable for observations of high-excitation CO lines and [N II] 205 um lines, which are scientifically relevant tracers of star formation and stellar feedback. Those lines are the brightest lines in the THz windows, so that they are suitable for the initiation of ground-based THz observations. THz polarization of star-forming regions can also be explored since it traces the dust population contributing to the THz spectral peak. For survey-type observations, we focus on ``sub-THz'' extragalactic surveys, whose uniqueness is to detect galaxies at redshifts z ~ 1--2, where the dust emission per comoving volume is the largest in the history of the Universe. Finally we explore possibilities of flexible time scheduling, which enables us to monitor active galactic nuclei, and to target gamma-ray burst afterglows. For these objects, THz and submillimeter wavelength ranges have not yet been explored.

The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI). The Spanish SKA White Book is the result of the coordinated effort of 120 astronomers from 40 different research centers. The book shows the enormous scientific interest of the Spanish astronomical community in the SKA and warrants an optimum scientific exploitation of the SKA by Spanish researchers, if Spain enters the SKA project.

A two dimensional hydrochemical hybrid code, KM2, is constructed to deal with astrophysical problems that would require coupled hydrodynamical and chemical evolution. The code assumes axisymmetry in cylindrical coordinate system, and consists of two modules: a hydrodynamics module and a chemistry module. The hydrodynamics module solves hydrodynamics using a Godunov-type finite volume scheme and treats included chemical species as passively advected scalars. The chemistry module implicitly solves non-equilibrium chemistry and change of the energy due to thermal processes with transfer of external ultraviolet radiation. Self-shielding effects on photodissociation of CO and H$_2$ are included. In this introductory paper, the adopted numerical method is presented, along with code verifications using the hydrodynamics modules, and a benchmark on the chemistry module with reactions specific to a photon-dominated region (PDR). Finally, as an example of the expected capability, the hydrochemical evolution of a PDR is presented based on the PDR benchmark.

We observed with the JVLA at 3.6 and 1.3 cm a sample of 11 proto-brown dwarf candidates in Taurus in a search for thermal radio jets driven by the most embedded brown dwarfs. We detected for the first time four thermal radio jets in proto-brown dwarf candidates. We compiled data from UKIDSS, 2MASS, Spitzer, WISE and Herschel to build the Spectral Energy Distribution (SED) of the objects in our sample, which are similar to typical Class~I SEDs of Young Stellar Objects (YSOs). The four proto-brown dwarf candidates driving thermal radio jets also roughly follow the well-known trend of centimeter luminosity against bolometric luminosity determined for YSOs, assuming they belong to Taurus, although they present some excess of radio emission compared to the known relation for YSOs. Nonetheless, we are able to reproduce the flux densities of the radio jets modeling the centimeter emission of the thermal radio jets using the same type of models applied to YSOs, but with corresponding smaller stellar wind velocities and mass-loss rates, and exploring different possible geometries of the wind or outflow from the star. Moreover, we also find that the modeled mass outflow rates for the bolometric luminosities of our objects agree reasonably well with the trends found between the mass outflow rates and bolometric luminosities of YSOs, which indicates that, despite the "excess" centimeter emission, the intrinsic properties of proto-brown dwarfs are consistent with a continuation of those of very low mass stars to a lower mass range. Overall, our study favors the formation of brown dwarfs as a scaled-down version of low-mass stars.

We present the detection and analysis of molecular hydrogen emission toward ten interstellar regions in the Large Magellanic Cloud. We examined low-resolution infrared spectral maps of twelve regions obtained with the Spitzer infrared spectrograph (IRS). The pure rotational 0--0 transitions of H$_2$ at 28.2 and 17.1${\,\rm \mu m}$ are detected in the IRS spectra for ten regions. The higher level transitions are mostly upper limit measurements except for three regions, where a 3$\sigma$ detection threshold is achieved for lines at 12.2 and 8.6${\,\rm \mu m}$. The excitation diagrams of the detected H$_2$ transitions are used to determine the warm H$_2$ gas column density and temperature. The single-temperature fits through the lower transition lines give temperatures in the range $86-137\,{\rm K}$. The bulk of the excited H$_2$ gas is found at these temperatures and contributes $\sim$5-17% to the total gas mass. We find a tight correlation of the H$_2$ surface brightness with polycyclic aromatic hydrocarbon and total infrared emission, which is a clear indication of photo-electric heating in photodissociation regions. We find the excitation of H$_2$ by this process is equally efficient in both atomic and molecular dominated regions. We also present the correlation of the warm H$_2$ physical conditions with dust properties. The warm H$_2$ mass fraction and excitation temperature show positive correlations with the average starlight intensity, again supporting H$_2$ excitation in photodissociation regions.

We present the various science cases for building Band 1 receivers as part of ALMA's ongoing Development Program. We describe the new frequency range for Band 1 of 35-52 GHz, a range chosen to maximize the receiver suite's scientific impact. We first describe two key science drivers: 1) the evolution of grains in protoplanetary disks and debris disks, and 2) molecular gas in galaxies during the era of re-ionization. Studies of these topics with Band 1 receivers will significantly expand ALMA's Level 1 Science Goals. In addition, we describe a host of other exciting continuum and line science cases that require ALMA's high sensitivity and angular resolution. For example, ALMA Band 1 continuum data will probe the Sunyaev-Zel'dovich Effect in galaxy clusters, Very Small Grains and spinning dust, ionized jets from young stars, spatial and flaring studies of Sgr A*, the acceleration sites of solar flares, pulsar wind nebulae, radio supernovae, and X-ray binaries. Furthermore, ALMA Band 1 line data will probe chemical differentiation in cloud cores, complex carbon chain molecules, extragalactic radio recombination lines, masers, magnetic fields through Zeeman effect measurements, molecular outflows from young stars, the co-evolution of star formation and active galactic nuclei, and the molecular content of galaxies at z ~ 3. ALMA provides similar to better sensitivities than the JVLA over 35-50 GHz, with differences increasing with frequency. ALMA's smaller antennas and shorter baselines, greater number of baselines, and single-dish capabilities, however, give it a significant edge for observing extended emission, making wide-field maps (mosaics), or attaining high image fidelity, as required by the described science cases.

The transition between the proto-star, Class I, and the pre-main sequence star, Class II, phases is still one of the most uncertain, and important, stages in the knowledge of the process of formation of an individual star, because it is the stage that determines the final mass of the star. We observed the YSO "Holoea", associated with IRAS 05327+3404, which was classified as an object in transition between the Class I and Class II phases with several unusual properties, and appears to be surrounded by large amounts of circumstellar material. We used the SMA and BIMA telescopes at millimeter and sub-millimeter wavelengths to observe the dust continuum emission and the CO (1-0) and (2-1), HCO+ (1-0) and (3-2), and HCN (1-0) transitions in the region around IRAS 05327+3404. We detected two continuum emission peaks at 1.1-mm: SMM 1, the sub-mm counterpart of IRAS 05327+3404, and SMM 2, ~6 arcsec to the West. The emissions of the three molecules show marked differences. The CO emission near the systemic velocity is filtered out by the telescopes, and CO mostly traces the high-velocity gas. The HCO+ and HCN emissions are more centrally concentrated around the central parts of the region, and show several intensity peaks coincident with the sub-mm continuum peaks. We identify two main molecular outflows: a bipolar outflow in an E-W direction that would be powered by SMM 1 and another one in a NE direction, which we associate with SMM 2. We propose that the SMM sources are probably Class I objects, with SMM 1 in an earlier evolutionary stage.

Physical conditions in dense and cold regions of interstellar clouds favour the formation of ice mantles on the surfaces of interstellar grains. It is predicted that most of the gaseous species heavier than H2 or He will adsorb onto the grains and will disappear from the gas-phase, changing its chemistry, within 10^9/n_H years. Nonetheless, many molecules in molecular clouds are not completely depleted in timescales of 10^5 yr. Several speculative mechanisms have been proposed to explain why molecules stay in the gas phase, but up to now none are fully convincing. At the same time, these mechanisms are not mutually exclusive and we can still explore the effects of other possible processes. We speculate on the consequences of H2 coating of grains on the evaporation rates of adsorbed species. More experiments and simulations are needed to calculate the evaporation rate Eevap(X-H2).

The Pipe nebula is a massive, nearby, filamentary dark molecular cloud with a low star-formation efficiency threaded by a uniform magnetic field perpendicular to its main axis. It harbors more than a hundred, mostly quiescent, very chemically young starless cores. The cloud is, therefore, a good laboratory to study the earliest stages of the star-formation process. We aim to investigate the primordial conditions and the relation among physical, chemical, and magnetic properties in the evolution of low-mass starless cores. We used the IRAM 30-m telescope to map the 1.2 mm dust continuum emission of five new starless cores, which are in good agreement with previous visual extinction maps. For the sample of nine cores, which includes the four cores studied in a previous work, we derived a Av to NH2 factor of (1.27$\pm$0.12)$\times10^{-21}$ mag cm$^{2}$ and a background visual extinction of ~6.7 mag possibly arising from the cloud material. We derived an average core diameter of ~0.08 pc, density of ~10$^5$ cm$^{-3}$, and mass of ~1.7 Msun. Several trends seem to exist related to increasing core density: (i) diameter seems to shrink, (ii) mass seems to increase, and (iii) chemistry tends to be richer. No correlation is found between the direction of the surrounding diffuse medium magnetic field and the projected orientation of the cores, suggesting that large scale magnetic fields seem to play a secondary role in shaping the cores. The full abstract is available in the pdf.

High-angular resolution observations of dense molecular cores show that these cores can be clumpier at smaller scales, and that some of these clumps can also be unbound or transient. The use of chemical models of the evolution of the molecular gas provides a way to probe the physical properties of the clouds. We study the properties of the clump and inter-clump medium in the starless CS core in LDN 673 by carrying out a molecular line survey with the IRAM 30-m telescope toward two clumps and two inter-clump positions. We also observed the 1.2-mm continuum with the MAMBO-II bolometer at IRAM. The dust continuum map shows four condensations, three of them centrally peaked, coinciding with previously identified sub-millimetre sources. We confirm that the denser clump of the region, $n\sim3.6 \times10^5$\cmt, is also the more chemically evolved, and it could still undergo further fragmentation. The inter-clump medium positions are denser than previously expected, likely $n\sim1\times10^3$--1$\times10^4$\cmt\ due to contamination, and are chemically young, similar to the gas in the lower density clump position. We argue that the density contrast between these positions and their general young chemical age would support the existence of transient clumps in the lower density material of the core. We were also able to find reasonable fits of the observationally derived chemical abundances to models of the chemistry of transient clumps.

Context: Observations of small carbon-bearing molecules such as CCH, C4H, c-C3H2, and HCO in the Horsehead Nebula have shown these species to have higher abundances towards the edge of the source than towards the center. Aims: Given the determination of a wide range of values for zeta (s-1), the total ionization rate of hydrogen atoms, and the proposal of a column-dependent zeta(N_H), where N_H is the total column of hydrogen nuclei, we desire to determine if the effects of zeta(N_H) in a single object with spatial variation can be observable. We chose the Horsehead Nebula because of its geometry and high density. Method: We model the Horsehead Nebula as a near edge-on photon dominated region (PDR), using several choices for zeta, both constant and as a function of column. The column-dependent zeta functions are determined by a Monte Carlo model of cosmic ray penetration, using a steep power-law spectrum and accounting for ionization and magnetic field effects. We consider a case with low-metal elemental abundances as well as a sulfur-rich case. Results: We show that use of a column-dependent zeta(N_H) of 5(-15) s-1 at the surface and 7.5(-16) s-1 at Av = 10 on balance improves agreement between measured and theoretical molecular abundances, compared with constant values of zeta.

Localised regions of enhanced emission from HCO+, NH3 and other species near Herbig-Haro objects (HHOs) have been interpreted as arising in a photochemistry stimulated by the HHO radiation on high density quiescent clumps in molecular clouds. Static models of this process have been successful in accounting for the variety of molecular species arising ahead of the jet; however recent observations show that the enhanced molecular emission is widespread along the jet as well as ahead. Hence, a realistic model must take into account the movement of the radiation field past the clump. It was previously unclear as to whether the short interaction time between the clump and the HHO in a moving source model would allow molecules such as HCO+ to reach high enough levels, and to survive for long enough to be observed. In this work we model a moving radiation source that approaches and passes a clump. The chemical picture is qualitatively unchanged by the addition of the moving source, strengthening the idea that enhancements are due to evaporation of molecules from dust grains. In addition, in the case of several molecules, the enhanced emission regions are longer-lived. Some photochemically-induced species, including methanol, are expected to maintain high abundances for ~10,000 years.

The Pipe Nebula is a massive, nearby dark molecular cloud with a low star-formation efficiency which makes it a good laboratory to study the very early stages of the star formation process. The Pipe Nebula is largely filamentary, and appears to be threaded by a uniform magnetic field at scales of few parsecs, perpendicular to its main axis. The field is only locally perturbed in a few regions, such as the only active cluster forming core B59. The aim of this study is to investigate primordial conditions in low-mass pre-stellar cores and how they relate to the local magnetic field in the cloud. We used the IRAM 30-m telescope to carry out a continuum and molecular survey at 3 and 1 mm of early- and late-time molecules toward four selected starless cores inside the Pipe Nebula. We found that the dust continuum emission maps trace better the densest regions than previous 2MASS extinction maps, while 2MASS extinction maps trace better the diffuse gas. The properties of the cores derived from dust emission show average radii of ~0.09 pc, densities of ~1.3x10^5 cm^-3, and core masses of ~2.5 M_sun. Our results confirm that the Pipe Nebula starless cores studied are in a very early evolutionary stage, and present a very young chemistry with different properties that allow us to propose an evolutionary sequence. All of the cores present early-time molecular emission, with CS detections toward all the sample. Two of them, Cores 40 and 109, present strong late-time molecular emission. There seems to be a correlation between the chemical evolutionary stage of the cores and the local magnetic properties that suggests that the evolution of the cores is ruled by a local competition between the magnetic energy and other mechanisms, such as turbulence.

We intended to study the incidence and characteristics of water masers in the envelopes of stars in the post-AGB and PN evolutionary stages. We have used the 64-m antenna in Parkes (Australia) to search for water maser emission at 22 GHz, towards a sample of 74 sources with IRAS colours characteristic of post-AGB stars and PNe, at declination $< -32 \deg$. In our sample, 39% of the sources are PNe or PNe candidates, and 50% are post-AGB stars or post-AGB candidates. We have detected four new water masers, all of them in optically obscured sources: three in PNe candidates (IRAS 12405-6219, IRAS 15103-5754, and IRAS 16333-4807); and one in a post-AGB candidate (IRAS 13500-6106). The PN candidate IRAS 15103-5754 has water fountain characteristics, and it could be the first PN of this class found. We confirm the tendency suggested in Paper I that the presence of water masers in the post-AGB phase is favoured in obscured sources with massive envelopes. We propose an evolutionary scenario for water masers in the post-AGB and PNe stages, in which ``water fountain'' masers could develop during post-AGB and early PN stages. Later PNe would show lower velocity maser emission, both along jets and close to the central objects, with only the central masers remaining in more evolved PNe.

All the effort that the astrophysical community has put into the development of the Virtual Observatory (VO) has surpassed the non-return point: the VO is a reality today, and an initiative that will self-sustain, and to which all archival projects must adhere. We have started the design of the scientific archive for the DSS-63 70-m antenna at NASA's DSN station in Robledo de Chavela (Madrid). Here we show how we can use all VO proposed data models to build a VO-compliant single-dish, multiple-feed, radio astronomical archive data model (RADAMS) suitable for the archival needs of the antenna. We also propose an exhaustive list of Universal Content Descriptors (UCDs) and FITS keywords for all relevant metadata. We will further refine this data model with the experience that we will gain from that implementation.

We present a study of the chemistry of a dense photon-dominated region (PDR) using a time-dependent chemical model. Our major interest is to study the spatial distribution of complex molecules such as hydrocarbons and cyanopolyynes in the cool dense material bordering regions where star formation has taken place. Our standard model uses a homogeneous cloud of density 2x10e4 cm-3 and temperature T=40 K, which is irradiated by a far-ultraviolet radiation field of intermediate intensity, given by X=100. We find that over a range of times unsaturated hydrocarbons (e.g., C2H, C4H, C3H2) have relatively high fractional abundances in the more external layers of the PDR, whereas their abundances in the innermost layers are several orders of magnitudes lower. On the other hand, molecules that are typical of late-time chemistry are usually more abundant in the inner parts of the PDR. We also present results for models with different density, temperature, intensity of the radiation field and initial fractional abundances. Our results are compared with both high- and moderate-angular resolution observations of the Horsehead nebula. Our standard model is partially successful in reproducing the observations. Additional models run with different physical parameters are able to reproduce the abundance of many of the observed molecules, but we do not find a single model that fits all the observations at the same time. We discuss the suitability of a time-dependent model of a dense PDR such as ours as an estimator of the age of a PDR, provided that enough observational data exist.

We present Very Large Array (VLA) observations of H2O and OH masers, as well as radio continuum emission at 1.3 and 18 cm toward three sources previously cataloged as planetary nebulae (PNe) and in which single-dish detections of H2O masers have been reported: IRAS 17443-2949, IRAS 17580-3111, and IRAS 18061-2505. Our goal was to unambiguously confirm their nature as water-maser-emitting PNe, a class of objects of which only two bona-fide members were previously known. We detected and mapped H2O maser emission toward all three sources, while OH maser emission is detected in IRAS 17443-2949 and IRAS 17580-3111 as well as in other two objects within the observed fields: IRAS 17442-2942 (unknown nature) and IRAS 17579-3121 (also cataloged as a possible PN). We found radio continuum emission associated only with IRAS 18061-2505. Our results confirm IRAS 18061-2505 as the third known case of a PN associated with H2O maser emission. The three known water-maser-emitting PNe have clear bipolar morphologies, which suggests that water maser emission in these objects is related to non-spherical mass-loss episodes. We speculate that these bipolar PNe would have ``water-fountain'' Asymptotic Giant Branch (AGB) and post-AGB stars as their precursors. A note of caution is given for other objects that have been classified as OHPNe (objects with both OH maser and radio continuum emission, that could be extremely young PNe) based on single-dish observations, since interferometric data of both OH masers and continuum are necessary for a proper identification as members of this class.

Aims: We investigated the possible relationship between the evolutionary stage of post-AGB stars and planetary nebulae (PNe) and the presence of water masers in their envelopes. Methods: We have used NASA's 70-m antenna in Robledo de Chavela (Spain) to search for the water maser transition at 22235.08 MHz, towards a sample of 105 sources with IRAS colour characteristic of post-AGB stars and PNe at declination >-32 deg. 83% of the sources in the sample are post-AGB stars, 15% PNe or PN candidates, while only 2% seem to be HII regions. Results: We have detected five water masers, of which four are reported for the first time: two in PNe (IRAS 17443-2949 and IRAS 18061-2505), a ``water fountain'' in a post-AGB star (IRAS 16552-3050), and one in a source previously catalogued as a PN, but whose classification is uncertain (IRAS 17580-3111). Conclusions: The unexpected detections of water masers in two objects among the small subset of PNe led us to suggest that the PNe harbouring water masers are a special type of massive, rapidly evolving PNe.

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.

Continuous-time, random-walk Monte Carlo simulations of H2 formation on grains have been performed for surfaces that are stochastically heated by photons. We have assumed diffuse cloud conditions and used a variety of grains of varying roughness and size based on olivine. The simulations were performed at different optical depths. We confirmed that small grains (r <= 0.02 micron) have low modal temperatures with strong fluctuations, which have a large effect on the efficiency of the formation of molecular hydrogen. The grain size distribution highly favours small grains and therefore H2 formation on these particles makes a large contribution to the overall formation rate for all but the roughest surfaces. We find that at A_V=0 only the roughest surfaces can produce the required amount of molecular hydrogen, but by A_V=1, smoother surfaces are possible alternatives. Use of a larger value for the evaporation energy of atomic hydrogen, but one still consistent with experiment, allows smoother surfaces to produce more H2.

We present FCRAO maps as well as combined BIMA and FCRAO maps of the high density molecular emission towards the CS core in the L673 region. With the FCRAO telescope, we mapped the emission in the CS(2-1), C34S(2-1), HCO+(1-0), and H13CO+(1-0) lines. The high density molecular emission, which arises from a filamentary structure oriented in the NW-SE direction, shows clear morphological differences for each molecule. We find that HCO+ has an extremely high optical depth, and that the H13CO+ emission is well correlated with submm sources. The BIMA and FCRAO combined maps recover emission from a lot of other structure which was previously undetected or only marginally detected, and show an overall aspect of a filamentary structure connecting several intense clumps. We found a total 15 clumps in our combined data cube, all of them resolved by our angular resolution, with diameters in the 0.03-0.09 pc range. We find a clear segregation between the northern and southern region of the map: the northern section shows the less chemically evolved gas and less massive but more numerous clumps, while the southern region is dominated by the largest and most massive clump, and contains the more evolved gas, as traced by emission of late-time molecules. We find that the derived clump masses are below the virial mass, and that the clumps masses become closer to the virial mass when they get bigger and more massive. This supports the idea that these clumps must be transient, and only the more massive ones have a chance to last long enough to form stars. The clumps we detect are probably in an earlier evolutionary stage than the ``starless cores'' reported recently in the literature. Only the most massive one has properties similar to a ``starless core''.

A multitransitional study with the BIMA interferometric array was carried out toward the starless core found in the L673 region, in order to study the small-size structure of the cores detected with previous single--dish observations, which provides us with a test of the predictions of the chemical model of Taylor et al. (1996; 1998). We detected emission in the CS (2-1), N2H+ (1-0), and HCO+ (1-0) lines. Several clumps of size ~0.08 pc were found for each line distributed all over the region where previous single-dish emission was found (Morata et al. 1997). Each molecular transition traces differently the clump distribution, although in some cases the detected clumps are coincident. The distribution of the N2H+ emission and the single-dish NH3 emission are coincident and compatible with an origin in the same gas. The large fraction of missing flux measured for the CS (2-1) transition can be explained if the cloud is formed by a clumpy and heterogeneous medium. Four positions were selected to derive the abundance ratios [N2H+/CS] and [HCO+/CS] from the molecular column density determinations, and to compare them with the values predicted by the chemical model. The model was able to explain the interferometric observations, and, in particular, the chemical differentiation of the detected clumps and the coincidence of the NH3 and N2H+ emissions. The lack of HCO+ towards the two selected positions that trace the more evolved clumps cannot be accounted for by the model, but it is possibly due to strong self-absorption. We propose a classification of the studied clumps according to the stage of chemical evolution indicated by the molecular abundances.