The velocity dispersion of giant molecular clouds. II - Mathematical and numerical refinements

The energy input to giant molecular clouds is recalculated, using the proper linearized equations of motion, including the Coriolis force and allowing for changes in the guiding center. Perturbation theory yields a result in the limit of distant encounters and small initial epicyclic amplitudes. Direct integration of the motion equations allows the strong encounter regime to be studied. The present perturbation theory result differs by a factor of order unity from that of Jog and Ostriker (1988). The result of present numerical integrations for the 2D (planar) velocity dispersion is presented. The accretion rate for a molecular cloud in the Galactic disk is calculated.

The Dissociative Recombination of Protonated Acrylonitrile, CH2CHCNH+, with Implications for Nitrile Chemistry in Dark Molecular Clouds and the Upper Atmosphere of Titan

Measurements on the dissociative recombination (DR) of protonated acrylonitrile, CH2CHCNH+, have been performed at the heavy ion storage ring CRYRING located in the Manne Siegbahn Laboratory in Stockholm, Sweden. It has been found that at~2meV relative kinetic energy about 50% of the DR events involve only ruptures of X–H bonds (where X=C or N)while the rest leads to the production of a pair of fragments each containing two heavy atoms (alongside H and/or H2). The absolute DR cross section has been investigated for relative kinetic energies ranging from ~1 meV to 1 eV. The thermal rate coefficient has been determined to follow the expression k(T) = 1.78 × 10-6 (T/300)-0.80 cm3 s-1 for electron temperatures ranging from ~10 to 1000 K. Gas-phase models of the nitrile chemistry in the dark molecular cloud TMC-1 have been run and results are compared with observations. Also, implications of the present results for the nitrile chemistry of Titan’s upper atmosphere are discussed.

Modelling of deuterium chemistry and its application to molecular clouds

We have developed new models of the chemistry of deuterium for investigating fractionation in interstellar molecular clouds. We have incorporated the latest information on reactions which affect deuteration, extended previous models to include S-D bonds for the first time and included the gasphase chemistry of some doubly-deuterated species. We present models for a wide range of physical parameters, including density, temperature, elemental abundances, and the freeze out of molecules on to dust grains. We discuss the detailed fractionation of particular species and show how fractionation can be used to probe the history of interstellar matter. The freeze out of molecules onto dust leads to significant enhancement in fractionation ratios and, in particular, to large fractionation in doubly-deuterated species.

SO and CS observations of molecular clouds - II. Analysis and modelling of the abundance ratios - probing O-2/CO with SO/CS?

We here analyse the observational SO and CS data presented in Nilsson ct al. (2000). The SO/CS integrated intensity ratio maps are presented for 19 molecular clouds, together with tables of relevant ratios at strategic positions, where we have also observed (SO)-S-34 and/or (CS)-S-34. The SO/CS abundance ratio as calculated from an LTE analysis is highly varying within and between the sources. Our isotopomer observations and Monte Carlo simulations verify that this is not an artifact due to optical depth problems. The variation of the maximum SO/CS abundance ratio between the clouds is 0.2-7. The largest variations within a cloud are found for the most nearby objects, possibly indicating resolution effects. We have also performed time dependent chemical simulations. We compare the simulations with our observed SO/CS abundance ratios and suggest a varying oxygen to carbon initial abundance, differing temporal evolution, density differences and X-ray sources associated with young stellar objects as possible explanations to the variations. In particular, the observed variation of the maximum SO/CS abundance ratio between the clouds can be explained by using initial O/C+ abundance ratios in the range 1.3-2.5. We finally derive a relationship between the SO/CS and O-2/CO abundance ratios, which may be used as a guide to find the most promising interstellar O-2 search targets.

Desorption processes and the deuterium fractionation in molecular clouds

Several processes have been suggested as ways of returning accreted grain mantles to the gas, thus preventing the total removal of molecules from the gas phase in dark quiescent clouds. We attempt to distinguish between them by considering not only the calculated gas-phase abundances, but also the ratio of the abundances of deuterated species to non-deuterated species. We find that the D/H ratio in molecules is relatively model-independent, but that desorption due to the formation of H-2 on grains gives the best overall agreement with the observations.

PROTONATED HCN IN MOLECULAR CLOUDS

Observations of protonated HCN (HCNH+) in a selection of galactic molecular clouds are reported. This species plays a key role in understanding the chemistry of the important high density tracer HCN. HCNH+ has been detected in the nearby cold dust cloud TMC-1 with an ratio relative to HCN of [HCNH+]/[HCN] between 0.015 and 0.26 (preferred value 0.03) and tentatively in DR21(OH) with a ratio of approximately 0.01. This is about 100 times higher than the ratio of protonated carbon monoxide to CO [HCO+]/[CO], but comparable to the [HCS+]/[CS] ratio. Possible explanations of these high abundance ratios are discussed in the light of model calculations.

GEMINI near-infrared spectroscopic observations of young massive stars embedded in molecular clouds

K-band spectra of young stellar candidates in four Southern hemisphere clusters have been obtained with the Gemini Near-Infrared Spectrograph in Gemini South. The clusters are associated with IRAS sources that have colours characteristic of ultracompact H II regions. Spectral types were obtained by comparison of the observed spectra with those of a near-infrared (NIR) library; the results include the spectral classification of nine massive stars and seven objects confirmed as background late-type stars. Two of the studied sources have K-band spectra compatible with those characteristic of very hot stars, as inferred from the presence of C IV, N III and N V emission lines at 2.078, 2.116 and 2.100 mu m, respectively. One of them, I16177_IRS1, has a K-band spectrum similar to that of Cyg OB2 7, an O3If* supergiant star. The nebular K-band spectrum of the associated Ultra-Compact (UC) H II region shows the s-process [Kr III] and [Se IV] high excitation emission lines, previously identified only in planetary nebula. One young stellar object was found in each cluster, associated with either the main IRAS source or a nearby resolved Midecourse Space eXperiment (MSX) component, confirming the results obtained from previous NIR photometric surveys. The distances to the stars were derived from their spectral types and previously determined JHK magnitudes; they agree well with the values obtained from the kinematic method, except in the case of IRAS 15408-5356, for which the spectroscopic distance is about a factor of 2 smaller than the kinematic value.

SPH simulations of clumps formation by dissipative collisions of molecular clouds - II. Magnetic case

We performed computer simulations of interstellar cloud-cloud collisions using the three-dimensional smoothed particle magnetohydrodynamics method. In order to study the role of the magnetic field on the process of collision-triggered fragmentation, we focused our attention on head-on supersonic collisions between two identical spherical molecular-clouds. Two extreme configurations of the magnetic field were adopted: parallel and perpendicular to the initial clouds motion. The initial magnetic field strength was approximately 12.0 muG. In the parallel case, much more of the collision debris were retained in the shocking region than in the non-magnetic case where gas escaped freely throughout the symmetry plane. Differently from the non-magnetic case, eddy-like vortices were formed. The regions of highest vorticity and the the regions of highest density are offset. We found clumps formation only in the parallel case, however, they were larger, hotter and less dense than in the analogous non-magnetic case. In the perpendicular case, the compressed field works as a magnetic wall, preventing a stronger compression of the colliding clouds. This last effect inhibits direct contact of the two clouds. In both cases, we found that the field lines show a chaotic aspect in large scales. Also, the field magnitude is considerably amplified in the shock layer. However, the field distribution is almost coherent in the higher density regions.

SPH simulations of clumps formation by dissipative collision of molecular clouds - I. Non magnetic case

Computer experiments of interstellar cloud collisions were performed with a new smoothed-particle-hydrodynamics (SPH) code. The SPH quantities were calculated by using spatially adaptive smoothing lengths and the SPH fluid equations of motion were solved by means of a hierarchical multiple time-scale leapfrog. Such a combination of methods allows the code to deal with a large range of hydrodynamic quantities. A careful treatment of gas cooling by H, H(2), CO and H II, as well as a heating mechanism by cosmic rays and by H(2) production on grains surface, were also included in the code. The gas model reproduces approximately the typical environment of dark molecular clouds. The experiments were performed by impinging two dynamically identical spherical clouds onto each other with a relative velocity of 10 km s(-1) but with a different impact parameter for each case. Each object has an initial density profile obeying an r(-1)-law with a cutoff radius of 10 pc and with an initial temperature of 20 K. As a main result, cloud-cloud collision triggers fragmentation but in expense of a large amount of energy dissipated, which occurred in the head-on case only. Off-center collision did not allow remnants to fragment along the considered time (similar to 6 Myr). However, it dissipated a considerable amount of orbital energy. Structures as small as 0.1 pc, with densities of similar to 10(4) cm(-3), were observed in the more energetic collision.

Molecular line and continuum studies of the early stages of star formation

New stars form in dense interstellar clouds of gas and dust called molecular clouds. The actual sites where the process of star formation takes place are the dense clumps and cores deeply embedded in molecular clouds. The details of the star formation process are complex and not completely understood. Thus, determining the physical and chemical properties of molecular cloud cores is necessary for a better understanding of how stars are formed. Some of the main features of the origin of low-mass stars, like the Sun, are already relatively well-known, though many details of the process are still under debate. The mechanism through which high-mass stars form, on the other hand, is poorly understood. Although it is likely that the formation of high-mass stars shares many properties similar to those of low-mass stars, the very first steps of the evolutionary sequence are unclear. Observational studies of star formation are carried out particularly at infrared, submillimetre, millimetre, and radio wavelengths. Much of our knowledge about the early stages of star formation in our Milky Way galaxy is obtained through molecular spectral line and dust continuum observations. The continuum emission of cold dust is one of the best tracers of the column density of molecular hydrogen, the main constituent of molecular clouds. Consequently, dust continuum observations provide a powerful tool to map large portions across molecular clouds, and to identify the dense star-forming sites within them. Molecular line observations, on the other hand, provide information on the gas kinematics and temperature. Together, these two observational tools provide an efficient way to study the dense interstellar gas and the associated dust that form new stars. The properties of highly obscured young stars can be further examined through radio continuum observations at centimetre wavelengths. For example, radio continuum emission carries useful information on conditions in the protostar+disk interaction region where protostellar jets are launched. In this PhD thesis, we study the physical and chemical properties of dense clumps and cores in both low- and high-mass star-forming regions. The sources are mainly studied in a statistical sense, but also in more detail. In this way, we are able to examine the general characteristics of the early stages of star formation, cloud properties on large scales (such as fragmentation), and some of the initial conditions of the collapse process that leads to the formation of a star. The studies presented in this thesis are mainly based on molecular line and dust continuum observations. These are combined with archival observations at infrared wavelengths in order to study the protostellar content of the cloud cores. In addition, centimetre radio continuum emission from young stellar objects (YSOs; i.e., protostars and pre-main sequence stars) is studied in this thesis to determine their evolutionary stages. The main results of this thesis are as follows: i) filamentary and sheet-like molecular cloud structures, such as infrared dark clouds (IRDCs), are likely to be caused by supersonic turbulence but their fragmentation at the scale of cores could be due to gravo-thermal instability; ii) the core evolution in the Orion B9 star-forming region appears to be dynamic and the role played by slow ambipolar diffusion in the formation and collapse of the cores may not be significant; iii) the study of the R CrA star-forming region suggests that the centimetre radio emission properties of a YSO are likely to change with its evolutionary stage; iv) the IRDC G304.74+01.32 contains candidate high-mass starless cores which may represent the very first steps of high-mass star and star cluster formation; v) SiO outflow signatures are seen in several high-mass star-forming regions which suggest that high-mass stars form in a similar way as their low-mass counterparts, i.e., via disk accretion. The results presented in this thesis provide constraints on the initial conditions and early stages of both low- and high-mass star formation. In particular, this thesis presents several observational results on the early stages of clustered star formation, which is the dominant mode of star formation in our Galaxy.

High Rydberg state carbon recombination lines toward Cassiopeia A: 332 MHz VLA observations and comparison with H I and molecular lines

We present observations of the C270alpha carbon recombination line, 21 cm neutral hydrogen line and (CO)-C-12 (J = 1 --> 0) molecular line toward Cas A. A comparison of the distribution of recombination line optical depths over the face of Cas A with that of H I optical depths and molecular line emission favors the association of C270alpha regions with H I rather than molecular clouds. The association makes it possible to self-consistently determine several physical parameters of the clouds by combining the recombination line and 21 cm H I measurements.

Dissociative Recombination of Protonated Formic Acid: Implications for Molecular Cloud and Cometary Chemistry

At the heavy ion storage ring CRYRING in Stockholm, Sweden, we have investigated the dissociative recombination of DCOOD2+ at low relative kinetic energies, from ~1 meV to 1 eV. The thermal rate coefficient has been found to follow the expression k(T) = 8.43 × 10-7 (T/300)^-0.78 cm3 s-1 for electron temperatures, T, ranging from ~10 to ~1000 K. The branching fractions of the reaction have been studied at ~2 meV relative kinetic energy. It has been found that ~87% of the reactions involve breaking a bond between heavy atoms. In only 13% of the reactions do the heavy atoms remain in the same product fragment. This puts limits on the gas-phase production of formic acid, observed in both molecular clouds and cometary comae. Using the experimental results in chemical models of the dark cloud, TMC-1, and using the latest release of the UMIST Database for Astrochemistry improves the agreement with observations for the abundance of formic acid. Our results also strengthen the assumption that formic acid is a component of cometary ices.

DIFFUSION OF MAGNETIC FIELD AND REMOVAL OF MAGNETIC FLUX FROM CLOUDS VIA TURBULENT RECONNECTION

The diffusion of astrophysical magnetic fields in conducting fluids in the presence of turbulence depends on whether magnetic fields can change their topology via reconnection in highly conducting media. Recent progress in understanding fast magnetic reconnection in the presence of turbulence reassures that the magnetic field behavior in computer simulations and turbulent astrophysical environments is similar, as far as magnetic reconnection is concerned. This makes it meaningful to perform MHD simulations of turbulent flows in order to understand the diffusion of magnetic field in astrophysical environments. Our studies of magnetic field diffusion in turbulent medium reveal interesting new phenomena. First of all, our three-dimensional MHD simulations initiated with anti-correlating magnetic field and gaseous density exhibit at later times a de-correlation of the magnetic field and density, which corresponds well to the observations of the interstellar media. While earlier studies stressed the role of either ambipolar diffusion or time-dependent turbulent fluctuations for de-correlating magnetic field and density, we get the effect of permanent de-correlation with one fluid code, i.e., without invoking ambipolar diffusion. In addition, in the presence of gravity and turbulence, our three-dimensional simulations show the decrease of the magnetic flux-to-mass ratio as the gaseous density at the center of the gravitational potential increases. We observe this effect both in the situations when we start with equilibrium distributions of gas and magnetic field and when we follow the evolution of collapsing dynamically unstable configurations. Thus, the process of turbulent magnetic field removal should be applicable both to quasi-static subcritical molecular clouds and cores and violently collapsing supercritical entities. The increase of the gravitational potential as well as the magnetization of the gas increases the segregation of the mass and magnetic flux in the saturated final state of the simulations, supporting the notion that the reconnection-enabled diffusivity relaxes the magnetic field + gas system in the gravitational field to its minimal energy state. This effect is expected to play an important role in star formation, from its initial stages of concentrating interstellar gas to the final stages of the accretion to the forming protostar. In addition, we benchmark our codes by studying the heat transfer in magnetized compressible fluids and confirm the high rates of turbulent advection of heat obtained in an earlier study.