In recent years, exciting developments have taken place in the identification of the role of cosmic rays in star-forming environments. Observations from radio to infrared wavelengths and theoretical modelling have shown that low-energy cosmic rays (<1 TeV) play a fundamental role in shaping the chemical richness of the interstellar medium, determining the dynamical evolution of molecular clouds. In this review we summarise in a coherent picture the main results obtained by observations and by theoretical models of propagation and generation of cosmic rays, from the smallest scales of protostars and circumstellar discs, to young stellar clusters, up to Galactic and extragalactic scales. We also discuss the new fields that will be explored in the near future thanks to new generation instruments, such as: CTA, for the $\gamma$-ray emission from high-mass protostars; SKA and precursors, for the synchrotron emission at different scales; and ELT/HIRES, JWST, and ARIEL, for the impact of cosmic rays on exoplanetary atmospheres and habitability.

Context. Submillimeter Array (SMA) 870 micron polarization observations of the hot molecular core G31.41+0.31 revealed one of the clearest examples up to date of an hourglass-shaped magnetic field morphology in a high-mass star-forming region. Aims. To better establish the role that the magnetic field plays in the collapse of G31.41+0.31, we carried out Atacama Large Millimeter/submillimeter Array (ALMA) observations of the polarized dust continuum emission at 1.3 mm with an angular resolution four times higher than that of the previous (sub)millimeter observations to achieve an unprecedented image of the magnetic field morphology. Methods. We used ALMA to perform full polarization observations at 233 GHz (Band 6). The resulting synthesized beam is 0.28"x0"20 which, at the distance of the source, corresponds to a spatial resolution of ~875 au. Results. The observations resolve the structure of the magnetic field in G31.41+0.31 and allow us to study the field in detail. The polarized emission in the Main core of G31.41+0.41is successfully fit with a semi-analytical magnetostatic model of a toroid supported by magnetic fields. The best fit model suggests that the magnetic field is well represented by a poloidal field with a possible contribution of a toroidal component of ~10% of the poloidal component, oriented southeast to northwest at ~ -44 deg and with an inclination of ~-45 degr. The magnetic field is oriented perpendicular to the northeast to southwest velocity gradient detected in this core on scales from 1E3-1E4 au. This supports the hypothesis that the velocity gradient is due to rotation and suggests that such a rotation has little effect on the magnetic field. The strength of the magnetic field estimated in the central region of the core with the Davis-Chandrasekhar-Fermi method is ~8-13 mG and implies that the mass-to-flux ratio in this region is slightly supercritical ...

Polarized dust emission outside of disks reveal the magnetic field morphology of molecular clouds. Within disks, however, polarized dust emission can arise from very different mechanisms (e.g., self-scattering), and each of them are useful for constraining physical properties in the disk. For example, these mechanisms allow us to constrain the disk grain size distributions and grain/disk geometries, independent from current methods of measuring these parameters. To accurately model these features and disentangle the various polarization mechanisms, multiwavelength observations at very high resolution and sensitivity are required. With significant upgrades to current interferometric facilities, we can understand how grains evolve in disks during the planet formation process.

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.

From molecular clouds to protoplanetary disks, non-ideal magnetic effects are important in many astrophysical environments. Indeed, in star and disk formation processes, it has become clear that these effects are critical to the evolution of the system. The efficacy of non-ideal effects are, however, determined by the complex interplay between magnetic fields, ionising radiation, cosmic rays, microphysics, and chemistry. In order to understand these key microphysical parameters, we present a one-dimensional non-ideal magnetohydrodynamics code and apply it to a model of a time-dependent, oblique, magnetic shock wave. By varying the microphysical ingredients of the model, we find that cosmic rays and dust play a major role, and that, despite the uncertainties, the inclusion of microphysics is essential to obtain a realistic outcome in magnetic astrophysical simulations.

We have recently obtained polarimetric data at mm wavelengths with ALMA for the young systems DG Tau and CW Tau, for which the rotation properties of jet and disk have been investigated in previous high angular resolution studies. The motivation was to test the models of magneto-centrifugal launch of jets via the determination of the magnetic configuration at the disk surface. The analysis of these data, however, reveals that self-scattering of dust thermal radiation dominates the polarization pattern. It is shown that even if no information on the magnetic field can be derived in this case, the polarization data are a powerful tool for the diagnostics of the properties and the evolution of dust in protoplanetary disks.

Understanding the role of magnetic fields in star-forming regions is of fundamental importance. In the near future, the exceptional sensitivity of SKA will offer a unique opportunity to evaluate the magnetic field strength in molecular clouds and cloud cores through synchrotron emission observations. The most recent Voyager 1 data, together with Galactic synchrotron emission and Alpha Magnetic Spectrometer data, constrain the flux of interstellar cosmic-ray electrons between $\approx3$ MeV and $\approx832$ GeV, in particular in the energy range relevant for synchrotron emission in molecular cloud cores at SKA frequencies. Synchrotron radiation is entirely due to primary cosmic-ray electrons, the relativistic flux of secondary leptons being completely negligible. We explore the capability of SKA in detecting synchrotron emission in two starless molecular cloud cores in the southern hemisphere, B68 and FeSt 1-457, and we find that it will be possible to reach signal-to-noise ratios of the order of $2-23$ at the lowest frequencies observable by SKA ($60-218$ MHz) with one hour of integration.

The interplay between chemistry and thermodynamics determines the final outcome of the process of gravitational collapse and sets the conditions for the formation of the first cosmological objects, including primordial supermassive black holes. In this chapter, we will review the main chemical reactions and the most important heating/cooling processes taking place in a gas of primordial composition, including the effects of local and cosmological radiation backgrounds.

The presence of small amounts of atomic hydrogen, detected as absorption dips in the 21 cm line spectrum, is a well-known characteristic of dark clouds. The abundance of hydrogen atoms measured in the densest regions of molecular clouds can be only explained by the dissociation of H$_2$ due to cosmic rays. We want to assess the role of Galactic cosmic rays in the formation of atomic hydrogen, by using recent developments in the characterisation of the low-energy spectra of cosmic rays and advances in the modelling of their propagation in molecular clouds. We model the attenuation of the interstellar cosmic rays entering a cloud and compute the dissociation rate of molecular hydrogen due to collisions with cosmic-ray protons and electrons as well as fast hydrogen atoms. We compare our results with the available observations. The cosmic-ray dissociation rate is entirely determined by secondary electrons produced in primary ionisation collisions. These secondary particles constitute the only source of atomic hydrogen at column densities above $\sim10^{21}$ cm$^{-2}$. We also find that the dissociation rate decreases with column density, while the ratio between the dissociation and ionisation rates varies between about 0.6 and 0.7. From comparison with observations we conclude that a relatively flat spectrum of interstellar cosmic-ray protons, as the one suggested by the most recent Voyager 1 data, can only provide a lower bound for the observed atomic hydrogen fraction. An enhanced spectrum of low-energy protons is needed to explain most of the observations. Our findings show that a careful description of molecular hydrogen dissociation by cosmic rays can explain the abundance of atomic hydrogen in dark clouds. An accurate characterisation of this process at high densities is crucial for understanding the chemical evolution of star-forming regions.

We present polarimetric data of CW Tau and DG Tau, two well-known Class II disk/jet systems, obtained with the Atacama Large Millimeter/submillimeter Array at 870 $\mu$m and 0."2 average resolution. In CW Tau, the total and polarized emission are both smooth and symmetric, with polarization angles almost parallel to the minor axis of the projected disk. In contrast, DG Tau displays a structured polarized emission, with an elongated brighter region in the disk's near side and a belt-like feature beyond about 0."3 from the source. At the same time the total intensity is spatially smooth, with no features. The polarization pattern, almost parallel to the minor axis in the inner region, becomes azimuthal in the outer belt, possibly because of a drop in optical depth. The polarization fraction has average values of 1.2% in CW Tau and 0.4% in DG Tau. Our results are consistent with polarization from self-scattering of the dust thermal emission. Under this hypothesis, the maximum size of the grains contributing to polarization is in the range 100 - 150 $\mu$m for CW Tau and 50 - 70 $\mu$m for DG Tau. The polarization maps combined with dust opacity estimates indicate that these grains are distributed in a geometrically thin layer in CW Tau, representing a settling in the disk midplane. Meanwhile, such settling is not yet apparent for DG Tau. These results advocate polarization studies as a fundamental complement to total emission observations, in investigations of the structure and the evolution of protoplanetary disks.

We use polarization observations of a circumbinary disk to investigate how the polarization properties change at distinct frequency bands. Our goal is to discern the main mechanism responsible for the polarization through comparison between our observations and model predictions. We used ALMA to perform full polarization observations at 97.5 GHz, 233 GHz and 343.5 GHz. The target is the Class I object BHB07-11, which is the youngest object in the Barnard 59 protocluster. Complementary VLA observations at 34.5 GHz revealed a binary system within the disk. We detect an extended and structured polarization pattern remarkably consistent among all three bands. The distribution of polarized intensity resembles a horseshoe shape with polarization angles following this morphology. From the spectral index between bands 3 and 7, we derive a dust opacity index $\beta \sim 1$ consistent with maximum grain sizes larger than expected to produce self-scattering polarization in each band. The polarization morphology do not match predictions from self-scattering. On the other hand, marginal correspondence is seen between our maps and predictions from radiation field assuming the brightest binary component as main radiation source. Molecular line data from BHB07-11 indicates disk rotation. We produced synthetic polarization maps from a rotating magnetized disk model assuming combined poloidal and toroidal magnetic field components. The magnetic field vectors (i. e., the polarization vectors rotated by 90$\degr$) are better represented by a model with poloidal magnetic field strength about 3 times the toroidal one. The similarity of our polarization patterns among the three bands provides a strong evidence against self-scattering and radiation fields. On the other hand, our data are reasonably well reproduced by a model of disk with toroidal magnetic field components slightly smaller than poloidal ones.

As part of our effort to search for circumstellar disks around high-mass stellar objects, we observed the well-known core G31.41+0.31 with ALMA at 1.4 mm with an angular resolution of~0.22" (~1700 au). The dust continuum emission has been resolved into two cores namely Main and NE. The Main core, which has the stronger emission and is the more chemically rich, has a diameter of ~5300 au, and is associated with two free-free continuum sources. The Main core looks featureless and homogeneous in dust continuum emission and does not present any hint of fragmentation. Each transition of CH3CN and CH3OCHO, both ground and vibrationally excited, as well as those of CH3CN isotopologues, shows a clear velocity gradient along the NE-SW direction, with velocity linearly increasing with distance from the center, consistent with solid-body rotation. However, when comparing the velocity field of transitions with different upper level energies, the rotation velocity increases with increasing energy of the transition, which suggests that the rotation speeds up towards the center. Spectral lines towards the dust continuum peak show an inverse P-Cygni profile that supports the existence of infall in the core. The infall velocity increases with the energy of the transition suggesting that the infall is accelerating towards the center of the core, consistent with gravitational collapse. Despite the monolithic appearance of the Main core, the presence of red-shifted absorption, the existence of two embedded free-free sources at the center, and the rotational spin-up are consistent with an unstable core undergoing fragmentation with infall and differential rotation due to conservation of angular momentum. Therefore, the most likely explanation for the monolithic morphology is that the large opacity of the dust emission prevents the detection of any inhomogeneity in the core.

Galactic cosmic rays are a ubiquitous source of ionisation of the interstellar gas, competing with UV and X-ray photons as well as natural radioactivity in determining the fractional abundance of electrons, ions and charged dust grains in molecular clouds and circumstellar discs. We model the propagation of different components of Galactic cosmic rays versus the column density of the gas. Our study is focussed on the propagation at high densities, above a few g cm$^{-2}$, especially relevant for the inner regions of collapsing clouds and circumstellar discs. The propagation of primary and secondary CR particles (protons and heavier nuclei, electrons, positrons, and photons) is computed in the continuous slowing down approximation, diffusion approximation, or catastrophic approximation, by adopting a matching procedure for the different transport regimes. A choice of the proper regime depends on the nature of the dominant loss process, modelled as continuous or catastrophic. The CR ionisation rate is determined by CR protons and their secondary electrons below $\approx 130$ g cm$^{-2}$ and by electron/positron pairs created by photon decay above $\approx600$ g cm$^{-2}$. We show that a proper description of the particle transport is essential to compute the ionisation rate in the latter case, since the electron/positron differential fluxes depend sensitively on the fluxes of both protons and photons. Our results show that the CR ionisation rate in high-density environments, like, e.g., the inner parts of collapsing molecular clouds or the mid-plane of circumstellar discs, is larger than previously assumed. It does not decline exponentially with increasing column density, but follows a more complex behaviour due to the interplay of different processes governing the generation and propagation of secondary particles.

Motivated by recent observational searches of sub-structure in starless molecular cloud cores, we investigate the evolution of density perturbations on scales smaller than the Jeans length embedded in contracting isothermal clouds, adopting the same formalism developed for the expanding Universe and the solar wind. We find that initially small amplitude, Jeans-stable perturbations (propagating as sound waves in the absence of a magnetic field), are amplified adiabatically during the contraction, approximately conserving the wave action density, until they either become nonlinear and steepen into shocks at a time $t_{\rm nl}$, or become gravitationally unstable when the Jeans length decreases below the scale of the perturbations at a time $t_{\rm gr}$. We evaluate analytically the time $t_{\rm nl}$ at which the perturbations enter the non-linear stage using a Burgers' equation approach, and we verify numerically that this time marks the beginning of the phase of rapid dissipation of the kinetic energy of the perturbations. We then show that for typical values of the rms Mach number in molecular cloud cores, $t_{\rm nl}$ is smaller than $t_{\rm gr}$, and therefore density perturbations likely dissipate before becoming gravitational unstable. Solenoidal modes grow at a faster rate than compressible modes, and may eventually promote fragmentation through the formation of vortical structures.

Contributed paper to the Conference "Francesco's Legacy: Star Formation in Space and Time", in memory of Francesco Palla (1954-2016), held in Florence, June 5-9, 2017

The possible occurrence of dead zones in protoplanetary disks subject to the magneto-rotational instability highlights the importance of disk ionization. We present a closed-form theory for the deep-down ionization by X-rays at depths below the disk surface dominated by far-ultraviolet radiation. Simple analytic solutions are given for the major ion classes, electrons, atomic ions, molecular ions and negatively charged grains. In addition to the formation of molecular ions by X-ray ionization of H2 and their destruction by dissociative recombination, several key processes that operate in this region are included, e.g., charge exchange of molecular ions and neutral atoms and destruction of ions by grains. Over much of the inner disk, the vertical decrease in ionization with depth into the disk is described by simple power laws, which can easily be included in more detailed modeling of magnetized disks. The new ionization theory is used to illustrate the non-ideal MHD effects of Ohmic, Hall and Ambipolar diffusion for a magnetic model of a T Tauri star disk using the appropriate Elsasser numbers.

Non-thermal photons deriving from radiative transitions among the internal ladder of atoms and molecules are an important source of photons in addition to thermal and stellar sources in many astrophysical environments. In the present work the calculation of reaction rates for the direct photodissociation of some molecules relevant in early Universe chemistry is presented; in particular, the calculations include non-thermal photons deriving from the recombination of primordial hydrogen and helium atoms for the cases of H2, HD and HeH+. New effects on the fractional abundances of chemical species are investigated and the fits for the HeH+ photodissociation rates by thermal photons are provided.

The radiative cooling of shocked gas with primordial chemical composition is an important process relevant to the formation of the first stars and structures, as well as taking place also in high velocity cloud collisions and supernovae explosions. Among the different processes that need to be considered, the formation kinetics and cooling of molecular hydrogen are of prime interest, since they provide the only way to lower the gas temperature to values well below $\sim$10$^4$~K. In previous works, the internal energy level structure of H$_2$ and its cation has been treated in the approximation of rovibrational ground state at low densities, or trying to describe the dynamics using some arbitrary $v>0$ H$_2$ level that is considered representative of the excited vibrational manifold. In this study, we compute the vibrationally resolved kinetics for the time-dependent chemical and thermal evolution of the post-shock gas in a medium of primordial composition. The calculated non-equilibrium distributions are used to evaluate effects on the cooling function of the gas and on the cooling time. Finally, we discuss the dependence of the results to different initial values of the shock velocity and redshift.

Our knowledge of the shape of radial metallicity gradients in disc galaxies has recently improved. Conversely, the understanding of their time evolution is more complex, since it requires analysis of stellar populations with different ages, or systematic studies of galaxies at different redshifts. In the Local Universe, Hii regions and planetary nebulae (PNe) are important tools to investigate it. We present an in-depth study of all nearby spiral galaxies (M33, M31, NGC300, and M81) with direct-method nebular abundances of both populations. For the first time, we also evaluate the radial migration of PN populations. We analyse Hii region and PN properties to: determine whether oxygen in PNe is a reliable tracer for past interstellar medium (ISM) composition; homogenise the published datasets; estimate the migration of the oldest stellar populations; determine the overall chemical enrichment and slope evolution. We confirm that oxygen in PNe is a reliable tracer for the past ISM metallicity. We find that PN gradients are flatter than or equal to those of Hii regions. When radial motions are negligible, this result provides a direct measurement of the time evolution of the gradient. For galaxies with dominant radial motions, we provide upper limits on the gradient evolution. Finally, the total metal content increases with time in all target galaxies, with early morphological type having a larger increment Delta(O/H) than late-type galaxies. Our findings provide important constraints to discriminate among different galactic evolutionary scenarios, favouring cosmological models with enhanced feedback from supernovae. The advent of extremely large telescopes will allow us to include galaxies in a wider range of morphologies and environments, thus putting firmer constraints to galaxy formation and evolution scenarios.

The local cosmic-ray (CR) spectra are calculated for typical characteristic regions of a cold dense molecular cloud, to investigate two so far neglected mechanisms of dust charging: collection of suprathermal CR electrons and protons by grains, and photoelectric emission from grains due to the UV radiation generated by CRs. The two mechanisms add to the conventional charging by ambient plasma, produced in the cloud by CRs. We show that the CR-induced photoemission can dramatically modify the charge distribution function for submicron grains. We demonstrate the importance of the obtained results for dust coagulation: While the charging by ambient plasma alone leads to a strong Coulomb repulsion between grains and inhibits their further coagulation, the combination with the photoemission provides optimum conditions for the growth of large dust aggregates in a certain region of the cloud, corresponding to the densities $n(\mathrm{H_2})$ between $\sim10^4$ cm$^{-3}$ and $\sim10^6$ cm$^{-3}$. The charging effect of CR is of generic nature, and therefore is expected to operate not only in dense molecular clouds but also in the upper layers and the outer parts of protoplanetary discs.

We investigate the role of the H_2^+ channel on H_2 molecule formation during the collapse of primordial gas clouds immersed in strong radiation fields which are assumed to have the shape of a diluted black-body spectra with temperature T_rad. Since the photodissociation rate of H_2^+ depends on its level population, we take full account of the vibrationally-resolved H_2^+ kinetics. We find that in clouds under soft but intense radiation fields with spectral temperature T_rad < 7000 K, the H_2^+ channel is the dominant H_2 formation process. On the other hand, for harder spectra with T_rad > 7000 K, the H^- channel takes over H_2^+ in the production of molecular hydrogen. We calculate the critical radiation intensity needed for supermassive star formation by direct collapse and examine its dependence on the H_2^+ level population. Under the assumption of local thermodynamic equilibrium (LTE) level population, the critical intensity is underestimated by a factor of a few for soft spectra with T_rad < 7000 K. For harder spectra, the value of the critical intensity is not affected by the level population of H_2^+. This result justifies previous estimates of the critical intensity assuming LTE populations since radiation sources like young and/or metal-poor galaxies are predicted to have rather hard spectra.

As long as magnetic fields remain frozen into the gas, the magnetic braking prevents the formation of protostellar discs. This condition is subordinate to the ionisation fraction characterising the inmost parts of a collapsing cloud. The ionisation level is established by the number and the energy of the cosmic rays able to reach these regions. Adopting the method developed in our previous studies, we computed how cosmic rays are attenuated as a function of column density and magnetic field strength. We applied our formalism to low- and high-mass star formation models obtained by numerical simulations of gravitational collapse that include rotation and turbulence. In general, we found that the decoupling between gas and magnetic fields, condition allowing the collapse to go ahead, occurs only when the cosmic-ray attenuation is taken into account with respect to a calculation in which the cosmic-ray ionisation rate is kept constant. We also found that the extent of the decoupling zone also depends on the dust grain size distribution and is larger if large grains (of radius about 0.1 microns) are formed by compression and coagulation during cloud collapse. The decoupling region disappears for the high-mass case due to magnetic field diffusion that is caused by turbulence and that is not included in the low-mass models. We infer that a simultaneous study of the cosmic-ray propagation during the cloud's collapse may lead to values of the gas resistivity in the innermost few hundred AU around a forming protostar that is higher than generally assumed.

Chemical models used to study the chemical composition of the gas and the ices in the interstellar medium are based on a network of chemical reactions and associated rate coefficients. These reactions and rate coefficients are partially compiled from data in the literature, when available. We present in this paper kida.uva.2014, a new updated version of the kida.uva public gas-phase network first released in 2012. In addition to a description of the many specific updates, we illustrate changes in the predicted abundances of molecules for cold dense cloud conditions as compared with the results of the previous version of our network, kida.uva.2011.

Cosmic rays are an important source of heating in the interstellar medium, in particular in dense molecular cloud cores shielded from the external ultraviolet radiation field. The limits placed on the cosmic-ray ionization rate from measurements of the gas temperature in dense clouds are unaffected by the uncertainties associated to the traditional methods based on the analysis of molecular abundances. However, high-resolution data are required to determine with sufficient accuracy the spatial temperature distribution within prestellar cores. In this contribution we illustrate in detail the case of the well-studied prestellar core L1544, showing that both its thermal structure and chemical composition are consistent with a cosmic ray ionization rate of $\sim 10^{-17}$~s$^{-1}$, significantly smaller than the value measured in the diffuse interstellar medium. We also briefly discuss possible applications of this method to molecular clouds of other galaxies

We discuss the effects of the magnetic field observed in molecular clouds on the process of star formation, concentrating on the phase of gravitational collapse of low-mass dense cores, cradles of sunlike stars. We summarize recent analytic work and numerical simulations showing that a substantial level of magnetic field diffusion at high densities has to occur in order to form rotationally supported disks. Furthermore, newly formed accretion disks are threaded by the magnetic field dragged from the parent core during the gravitational collapse. These disks are expected to rotate with a sub-Keplerian speed because they are partially supported by magnetic tension against the gravity of the central star. We discuss how sub-Keplerian rotation makes it difficult to eject disk winds and accelerates the process of planet migration. Moreover, magnetic fields modify the Toomre criterion for gravitational instability via two opposing effects: magnetic tension and pressure increase the disk local stability, but sub-Keplerian rotation makes the disk more unstable. In general, magnetized disks are more stable than their nonmagnetic counterparts; thus, they can be more massive and less prone to the formation of giant planets by gravitational instability.

The properties of filamentary interstellar clouds observed at sub-millimetre wavelengths, especially by the Herschel Space Observatory, are analysed with polytropic models in cylindrical symmetry. The observed radial density profiles are well reproduced by negative-index cylindrical polytropes with polytropic exponent $1/3\lesssim \gamma_{\rm p} \lesssim 2/3$ (polytropic index $-3\lesssim n \lesssim -3/2$), indicating either external heating or non-thermal pressure components. However, the former possibility requires unrealistically high gas temperatures at the filament's surface and is therefore very unlikely. Non-thermal support, perhaps resulting from a superposition of small-amplitude Alfv\'en waves (corresponding to $\gamma_{\rm p}=1/2$), is a more realistic possibility, at least for the most massive filaments. If the velocity dispersion scales as the square root of the density (or column density) on the filament's axis, as suggested by observations, then polytropic models are characterised by a uniform width. The mass per unit length of pressure-bounded cylindrical polytropes depends on the conditions at the boundary and is not limited as in the isothermal case. However, polytropic filaments can remain stable to radial collapse for values of the axis-to-surface density contrast as large as the values observed only if they are supported by a non-isentropic pressure component.

We study the properties of magnetised cylindrical polytropes as models for interstellar filamentary clouds, extending the analysis presented in a companion paper (Toci & Galli 2014a). We formulate the general problem of magnetostatic equilibrium in the presence of a helical magnetic field, with the aim of determining the degree of support or compression resulting from the magnetisation of the cloud. We derive scale-free solutions appropriate to describe the properties of the envelopes of filaments at radii larger than the flat-density region. In these solutions, the polytropic exponent determines the radial profiles of the density and the magnetic field. The latter decreases with radius less steeply than the density, and field lines are helices twisted over cylindrical surfaces. A soft equation of state supports magnetic configurations that preferentially compress and confine the filament, whereas in the isothermal limit the field provides support. For each value of the polytropic exponent, the Lorentz force is directed outward or inward depending on whether the pitch angle is below or above some critical value which is a function of the polytropic exponent only.

The formation of protostellar discs is severely hampered by magnetic braking, as long as magnetic fields remain frozen in the gas. The latter condition depends on the levels of ionisation that characterise the innermost regions of a collapsing cloud. The chemistry of dense cloud cores and, in particular, the ionisation fraction is largely controlled by cosmic rays. The aim of this paper is to evaluate whether the attenuation of the flux of cosmic rays expected in the regions around a forming protostar is sufficient to decouple the field from the gas, thereby influencing the formation of centrifugally supported disc. We adopted the method developed in a former study to compute the attenuation of the cosmic-ray flux as a function of the column density and the field strength in clouds threaded by poloidal and toroidal magnetic fields. We applied this formalism to models of low- and high-mass star formation extracted from numerical simulations of gravitational collapse that include rotation and turbulence. For each model we determine the size of the magnetic decoupling zone, where collapse or rotation motion becomes unaffected by the local magnetic field. In general, we find that decoupling only occurs when the attenuation of cosmic rays is taken into account with respect to a calculation in which the cosmic-ray ionisation rate is kept constant. The extent of the decoupling zone also depends on the dust grain size distribution and is larger if large grains (of radius $\sim 10^{-5}$ cm) are formed by compression and coagulation during cloud collapse. We conclude that a realistic treatment of cosmic-ray propagation and attenuation during cloud collapse may lead to a value of the resistivity of the gas in the innermost few hundred AU around a forming protostar that is higher than generally assumed.

The present study aims at characterizing the massive star forming region G35.20N, which is found associated with at least one massive outflow and contains multiple dense cores, one of them recently found associated with a Keplerian rotating disk. We used ALMA to observe the G35.20N region in the continuum and line emission at 350 GHz. The observed frequency range covers tracers of dense gas (e.g. H13CO+, C17O), molecular outflows (e.g. SiO), and hot cores (e.g. CH3CN, CH3OH). The ALMA 870 um continuum emission map reveals an elongated dust structure (0.15 pc long and 0.013 pc wide) perpendicular to the large-scale molecular outflow detected in the region, and fragmented into a number of cores with masses 1-10 Msun and sizes 1600 AU. The cores appear regularly spaced with a separation of 0.023 pc. The emission of dense gas tracers such as H13CO+ or C17O is extended and coincident with the dust elongated structure. The three strongest dust cores show emission of complex organic molecules characteristic of hot cores, with temperatures around 200 K, and relative abundances 0.2-2x10^(-8) for CH3CN and 0.6-5x10^(-6) for CH3OH. The two cores with highest mass (cores A and B) show coherent velocity fields, with gradients almost aligned with the dust elongated structure. Those velocity gradients are consistent with Keplerian disks rotating about central masses of 4-18 Msun. Perpendicular to the velocity gradients we have identified a large-scale precessing jet/outflow associated with core B, and hints of an east-west jet/outflow associated with core A. The elongated dust structure in G35.20N is fragmented into a number of dense cores that may form massive stars. Based on the velocity field of the dense gas, the orientation of the magnetic field, and the regularly spaced fragmentation, we interpret this elongated structure as the densest part of a 1D filament fragmenting and forming massive stars.

Cosmic rays are the main agents in controlling the chemical evolution and setting the ambipolar diffusion time of a molecular cloud. We summarise the processes causing the energy degradation of cosmic rays due to their interaction with molecular hydrogen, focusing on the magnetic effects that influence their propagation. Making use of magnetic field configurations generated by numerical simulations, we show that the increase of the field line density in the collapse region results in a reduction of the cosmic-ray ionisation rate. As a consequence the ionisation fraction decreases, facilitating the decoupling between the gas and the magnetic field.

Cosmic rays (CR) play an important role in dense molecular cores, affecting their thermal and dynamical evolution and initiating the chemistry. Several studies have shown that the formation of protostellar discs in collapsing clouds is severely hampered by the braking torque exerted by the entrained magnetic field on the infalling gas, as long as the field remains frozen to the gas. We examine the possibility that the concentration and twisting of the field lines in the inner region of collapse can produce a significant reduction of the ionisation fraction. To check whether the CR ionisation rate (CRir) can fall below the critical value required to maintain good coupling, we first study the propagation of CRs in a model of a static magnetised cloud varying the relative strength of the toroidal/poloidal components and the mass-to-flux ratio. We then follow the path of CRs using realistic magnetic field configurations generated by numerical simulations of a rotating collapsing core. We find that an increment of the toroidal component of the magnetic field, or, in general, a more twisted configuration of the field lines, results in a decrease in the CR flux. This is mainly due to the magnetic mirroring effect that is stronger where larger variations in the field direction are present. In particular, we find a decrease of the CRir below 10^-18 s-1 in the central 300-400 AU, where density is higher than about 10^9 cm-3. This very low value of the CRir is attained in the cases of intermediate and low magnetisation (mass-to-flux ratio lambda=5 and 17, respectively) and for toroidal fields larger than about 40% of the total field. Magnetic field effects can significantly reduce the ionisation fraction in collapsing clouds. We provide a handy fitting formula to compute approximately the attenuation of the CRir in a molecular cloud as a function of the density and the magnetic configuration.

The cosmological recombination of H and He at z \sim 1000 and the formation of H2 during the dark ages produce a non-thermal photon excess in the Wien tail of the cosmic microwave background (CMB) blackbody spectrum. Here we compute the effect of these photons on the H- photodetachment and H2+ photodissociation processes. We discuss the implications for the chemical evolution of the Universe in the post-recombination epoch, emphasizing how important a detailed account of the full vibrational manifold of H2 and H2+ in the chemical network is. We find that the final abundances of H2, H2+, H3+ and HD are significantly smaller than in previous calculations that neglected the effect of non-thermal photons. The suppression is mainly caused by extra hydrogen recombination photons and could affect the formation rate of first stars. We provide simple analytical approximations for the relevant rate coefficients and briefly discuss the additional effect of dark matter annihilation on the considered reaction rates.

Supermassive stars and quasi-stars (massive stars with a central black hole) are both considered as potential progenitors for the formation of supermassive black holes. They are expected to form from rapidly accreting protostars in massive primordial halos. We explore how long rapidly accreting protostars remain on the Hayashi track, implying large protostellar radii and weak accretion luminosity feedback. We assess the potential role of energy production in the nuclear core, and determine what regulates the evolution of such protostars into quasi-stars or supermassive stars. We follow the contraction of characteristic mass scales in rapidly accreting protostars, and infer the timescales for them to reach nuclear densities. We compare the characteristic timescales for nuclear burning with those for which the extended protostellar envelope can be maintained. We find that the extended envelope can be maintained up to protostellar masses of 3.6x10^8 \dot{m}^3 solar, where \dot{m} denotes the accretion rate in solar masses per year. We expect the nuclear core to exhaust its hydrogen content in 7x10^6 yrs. If accretion rates \dot{m}>>0.14 can still be maintained at this point, a black hole may form within the accreting envelope, leading to a quasi-star. Alternatively, the accreting object will gravitationally contract to become a main-sequence supermassive star. Due to the limited gas reservoir in dark matter halos with 10^7 solar masses, the accretion rate onto the central object may drop at late times, implying the formation of supermassive stars as the typical outcome of direct collapse. However, if high accretion rates are maintained, a quasi-star with an interior black hole may form.

Cosmic-rays constitute the main ionising and heating agent in dense, starless, molecular cloud cores. We reexamine the physical quantities necessary to determine the cosmic-ray ionisation rate (especially the cosmic ray spectrum at E < 1 GeV and the ionisation cross sections), and calculate the ionisation rate as a function of the column density of molecular hydrogen. Available data support the existence of a low-energy component (below about 100 MeV) of cosmic-ray electrons or protons responsible for the ionisation of diffuse and dense clouds. We also compute the attenuation of the cosmic-ray flux rate in a cloud core taking into account magnetic focusing and magnetic mirroring, following the propagation of cosmic rays along flux tubes enclosing different amount of mass and mass-to-flux ratios. We find that mirroring always dominates over focusing, implying a reduction of the cosmic-ray ionisation rate by a factor of 3-4 depending on the position inside the core and the magnetisation of the core.

We report on ALMA observations of continuum and molecular line emission with 0.4" resolution towards the high-mass star forming region G35.20-0.74 N. Two dense cores are detected in typical hot-core tracers, such as CH3CN, which reveal velocity gradients. In one of these cores, the velocity field can be fitted with an almost edge-on Keplerian disk rotating about a central mass of 18 Msun. This finding is consistent with the results of a recent study of the CO first overtone bandhead emission at 2.3mum towards G35.20-0.74 N. The disk radius and mass are >2500 au and 3 Msun. To reconcile the observed bolometric luminosity (3x10^4 Lsun) with the estimated stellar mass of 18 Msun, we propose that the latter is the total mass of a binary system.

Within the precise cosmological framework provided by the Lambda-Cold Dark Matter model and standard Big Bang nucleosynthesis, the chemical evolution of the pregalactic gas can now be followed with accuracy limited only by the uncertainties on the reaction rates. Starting during the recombination era, the formation of the first molecules and molecular ions containing hydrogen, deuterium, helium, and lithium was severely hindered by the low density of the expanding universe, the intensity of the cosmic radiation field, and the absence of solid catalyzers. Molecular hydrogen and deuterated hydrogen, the most abundant species formed in the gas phase prior to structure formation, played a fundamental role in the cooling of the gas clouds that gave birth to the first stellar generation, contributing to determine the scale of fragmentation. Primordial molecules also interacted with the photons of the cosmic background via resonant scattering, absorption and emission. In this review we examine the current status of the chemistry of the early universe and discuss the most relevant reactions for which uncertainties still exist from theory or laboratory experiments. The prospects for detecting spectral distortions or spatial anisotropies due to the first atoms and molecules are also addressed.

Most galaxies follow well-defined scaling relations of metallicity (O/H), star formation rate (SFR), and stellar mass. However, low-metallicity starbursts, rare in the Local Universe but more common at high redshift, deviate significantly from these scaling relations. On the "main sequence" of star formation, these galaxies have high SFR for a given M*; and on the mass-metallicity relation, they have excess M* for their low metallicity. In this paper, we characterize O/H, M*, and SFR for these deviant "low-metallicity starbursts", selected from a sample of ~1100 galaxies, spanning almost two orders of magnitude in metal abundance, a factor of ~10^6 in SFR, and of ~10^5 in stellar mass. Our sample includes quiescent star-forming galaxies and blue compact dwarfs at redshift 0, luminous compact galaxies at redshift 0.3, and Lyman Break galaxies at redshifts 1-3.4. Applying a Principal Component Analysis (PCA) to the galaxies in our sample with M*<10^{10} Msun gives a Fundamental Plane (FP) of scaling relations; SFR and stellar mass define the plane itself, and O/H its thickness. The dispersion for our sample in the edge-on view of the plane is 0.17 dex, independently of redshift and including the metal-poor starbursts. The same FP is followed by 55100 galaxies selected from the Sloan Digital Sky Survey, with a dispersion of 0.06dex. In a companion paper, we develop multi-phase chemical evolution models that successfully predict the observed scaling relations and the FP; the deviations from the main scaling relations are caused by a different (starburst or "active") mode of star formation. These scaling relations do not truly evolve, but rather are defined by the different galaxy populations dominant at different cosmological epochs.

Scaling relations of metallicity (O/H), star formation rate (SFR), and stellar mass give important insight on galaxy evolution. They are obeyed by most galaxies in the Local Universe and also at high redshift. In a companion paper, we compiled a sample of ~1100 galaxies from redshift 0 to ~3, spanning almost two orders of magnitude in metal abundance, a factor of $\sim10^6$ in SFR, and of ~10^5 in stellar mass. We have characterized empirically the star-formation "main sequence" (SFMS) and the mass-metallicity relation (MZR) for this sample, and also identified a class of low-metallicity starbursts, rare locally but more common in the distant universe. These galaxies deviate significantly from the main scaling relations, with high SFR and low metal content for a given M*. In this paper, we model the scaling relations and explain these deviations from them with a set of multi-phase chemical evolution models based on the idea that, independently of redshift, initial physical conditions in a galaxy's evolutionary history can dictate its location in the scaling relations. Our models are able to successfully reproduce the O/H, M*, and SFR scaling relations up to z~3, and also successfully predict the molecular cloud fraction as a function of stellar mass. These results suggest that the scaling relations are defined by different modes of star formation: an "active" starburst mode, more common at high redshift, and a quiescent "passive" mode that is predominant locally and governs the main trends.

Cosmic-ray and X-ray heating are derived from the electron energy loss calculations of Dalgarno, Yan and Liu for hydrogen-helium gas mixtures. These authors treated the heating from elastic scattering and collisional de-excitation of rotationally excited hydrogen molecules. Here we consider the heating that can arise from all ionization and excitation processes, with particular emphasis on the reactions of cosmic-ray and X-ray generated ions with the heavy neutral species, which we refer to as chemical heating. In molecular regions, chemical heating dominates and can account for 50 per cent of the energy expended in the creation of an ion pair. The heating per ion pair ranges in the limit of negligible electron fraction from about 4.3 eV for diffuse atomic gas, to about 13 eV for the moderately dense regions of molecular clouds and to about 18 eV for the very dense regions of protoplanetary disks. An important general conclusion of this study is that cosmic-ray and X-ray heating depends on the physical properties of the medium, i.e., on the molecular and electron fractions, the total density of hydrogen nuclei, and to a lesser extent on the temperature. It is also noted that chemical heating, the dominant process for cosmic-ray and X-ray heating, plays a role in UV irradiated molecular gas.

Context: Supersonic nonthermal motions in molecular clouds are often interpreted as long-lived magnetohydrodynamic (MHD) waves. The propagation and amplitude of these waves is affected by local physical characteristics, most importantly the gas density and the ionization fraction. Aims: We study the propagation, reflection and dissipation of Alfv\'en waves in molecular clouds deriving the behavior of observable quantities such as the amplitudes of velocity fluctuations and the rate of energy dissipation. Methods: We formulated the problem in terms of Els\"asser variables for transverse MHD waves propagating in a one-dimensional inhomogeneous medium, including the dissipation due to collisions between ions and neutrals and to a nonlinear turbulent cascade treated in a phenomenological way. We considered both steady-state and time-dependent situations and solved the equations of the problem numerically with an iterative method and a Lax-Wendroff scheme, respectively. Results: Alfv\'en waves incident on overdense regions with density profiles typical of cloud cores embedded in a diffuse gas suffer enhanced reflection in the regions of the steepest density gradient, and strong dissipation in the core's interior. These effects are especially significant when the wavelength is intermediate between the critical wavelength for propagation and the typical scale of the density gradient. For larger wave amplitudes and/or steeper input spectra, the effects of the perpendicular turbulent cascade result in a stronger energy dissipation in the regions immediately surrounding the dense core. Conclusions: The results may help to interpret the sharp decrease of line width observed in the environments of low-mass cloud cores in several molecular transitions.

We present the results of quantum calculations based on entirely ab initio methods for a variety of molecular processes and chemical reactions involving the LiHe$^+$ ionic polar molecule. With the aid of these calculations we derive accurate reaction rates and fitting expressions valid over a range of gas temperatures representative of the typical conditions of the pregalactic gas. With the help of a full chemical network, we then compute the evolution of the abundance of LiHe$^+$ as function of redshift in the early Universe. Finally, we compare the relative abundance of LiHe$^+$ with that of other polar cations formed in the same redshift interval.

We present a novel chemical database for gas-phase astrochemistry. Named the KInetic Database for Astrochemistry (KIDA), this database consists of gas-phase reactions with rate coefficients and uncertainties that will be vetted to the greatest extent possible. Submissions of measured and calculated rate coefficients are welcome, and will be studied by experts before inclusion into the database. Besides providing kinetic information for the interstellar medium, KIDA is planned to contain such data for planetary atmospheres and for circumstellar envelopes. Each year, a subset of the reactions in the database (kida.uva) will be provided as a network for the simulation of the chemistry of dense interstellar clouds with temperatures between 10 K and 300 K. We also provide a code, named Nahoon, to study the time-dependent gas-phase chemistry of 0D and 1D interstellar sources.

Energy exchange processes play a crucial role in the early Universe, affecting the thermal balance and the dynamical evolution of the primordial gas. In the present work we focus on the consequences of a non-thermal distribution of the level populations of H$_2$: first, we determine the excitation temperatures of vibrational transitions and the non-equilibrium heat transfer; second, we compare the modifications to chemical reaction rate coefficients with respect to the values obtained assuming local thermodynamic equilibrium; third, we compute the spectral distortions to the cosmic background radiation generated by the formation of H$_2$ in vibrationally excited levels. We conclude that non-equilibrium processes cannot be ignored in cosmological simulations of the evolution of baryons, although their observational signatures remain below current limits of detection. New fits to the equilibrium and non-equilibrium heat transfer functions are provided.

Young stars typically have strong magnetic fields, so that the magnetospheres of newly formed close binaries can interact, dissipate energy, and produce synchrotron radiation. The V773 Tau A binary system, a pair of T Tauri stars with a 51 day orbit, displays such a signature, with peak emission taking place near periastron. This paper proposes that the observed emission arises from the change in energy stored in the composite magnetic field of the system. We model the fields using the leading order (dipole) components and show that this picture is consistent with current observations. In this model, the observed radiation accounts for a fraction of the available energy of interaction between the magnetic fields from the two stars. Assuming antisymmetry, we compute the interaction energy $E_{\rm int}$ as a function of the stellar radii, the stellar magnetic field strengths, the binary semi-major axis, and orbital eccentricity, all of which can be measured independently of the synchrotron radiation. The variability in time and energetics of the synchrotron radiation depend on the details of the annihilation of magnetic fields through reconnection events, which generate electric fields that accelerate charged particles, and how those charged particles, especially fast electrons, are removed from the interaction region. However, the major qualitative features are well described by the background changes in the global magnetic configuration driven by the orbital motion. The theory can be tested by observing a collection of pre-main-sequence binary systems.

Observations of dust polarized emission toward star forming regions trace the magnetic field component in the plane of the sky and provide constraints to theoretical models of cloud collapse. We compare high-angular resolution observations of the submillimeter polarized emission of the low-mass protostellar source NGC 1333 IRAS 4A with the predictions of three different models of collapse of magnetized molecular cloud cores. We compute the Stokes parameters for the dust emission for the three models. We then convolve the results with the instrumental response of the Submillimeter Array observation toward IRAS 4A. Finally, we compare the synthetic maps with the data, varying the model parameters and orientation, and we assess the quality of the fit by a \chi^2 analysis. High-angular resolution observations of polarized dust emission can constraint the physical properties of protostars. In the case of IRAS 4A, the best agreements with the data is obtained for models of collapse of clouds with mass-to-flux ratio >2 times the critical value, initial uniform magnetic field of strength ~0.5 mG, and age of the order of a few 10^4 yr since the onset of collapse. Magnetic dissipation, if present, is found to occur below the resolution level of the observations. Including a previously measured temperature profile of IRAS 4A leads to a more realistic morphology and intensity distribution. We also show that ALMA has the capability of distinguishing among the three different models adopted in this work. Our results are consistent with the standard theoretical scenario for the formation of low-mass stars, where clouds initially threaded by large-scale magnetic fields become unstable and collapse, trapping the field in the nascent protostar and the surrounding circumstellar disk. In the collapsing cloud, the dynamics is dominated by gravitational and magnetic forces.

Low-energy cosmic rays are the dominant source of ionization for molecular cloud cores. The ionization fraction, in turn, controls the coupling of the magnetic field to the gas and hence the dynamical evolution of the cores. The purpose of this work is to compute the attenuation of the cosmic-ray flux rate in a cloud core taking into account magnetic focusing, magnetic mirroring, and all relevant energy loss processes. We adopt a standard cloud model characterized by a mass-to-flux ratio supercritical by a factor of about 2 to describe the density and magnetic field distribution of a low-mass starless core, and we follow the propagation of cosmic rays through the core along flux tubes enclosing different amount of mass. We then extend our analysis to cores with different mass-to-flux ratios. We find that mirroring always dominates over focusing, implying a reduction of the cosmic-ray ionization rate by a factor of about 2-3 over most of a solar-mass core with respect to the value in the intercloud medium outside the core. For flux tubes enclosing larger masses the reduction factor is smaller, since the field becomes increasingly uniform at larger radii and lower densities. We also find that the cosmic-ray ionization rate is further reduced in clouds with stronger magnetic field, e.g. by a factor of about 4 for a marginally critical cloud. The magnetic field threading molecular cloud cores affects the penetration of low-energy cosmic rays and reduces the ionization rate by a factor 3-4 depending on the position inside the core and the magnetization of the core.

We formulate a vibrationally resolved kinetics for molecular hydrogen and its cation in the primordial Universe chemistry. Formation, destruction and relaxation processes for each vibrational level are studied and included as chemical pathways of the present model. The fractional abundance of each vibrational level as a function of the redshift is given: a strong deviation from the Boltzmann distribution is found at low $z$. A discussion of the results is provided, also evaluating the effects of relaxation processes on the level populations. Analytical fits for some LTE rate coefficients are given in the Appendix.

The role of HeH+ has been newly assessed with the aid of newly calculated rates which use entirely ab initio methods, thereby allowing us to compute more accurately the relevant abundances within the global chemical network of the early universe. A comparison with the similar role of the ionic molecule LiH+ is also presented. Quantum calculations have been carried out for the gas-phase reaction of HeH+ with H atoms with our new in-house code, based on the negative imaginary potential method. Integral cross sections and reactive rate coefficients obtained under the general conditions of early universe chemistry are presented and discussed. With the new reaction rate, the abundance of HeH+ in the early universe is more than one order of magnitude larger than in previous studies. Our more accurate findings further buttress the possibility to detect cosmological signatures of HeH+.

The relative efficiencies of the chemical pathways that can lead to the destruction of LiH and LiH+ molecules, conjectured to be present in the primordial gas and to control molecular cooling processes in the gravitational collapse of the post-recombination era, are revisited by using accurate quantum calculations for the several reactions involved. The new rates are employed to survey the behavior of the relative abundance of these molecules at redshifts of interest for early universe conditions. We find significant differences with respect to previous calculations, the present ones yielding LIH abundances higher than LiH+ at all redshifts.

This paper considers gravitational perturbations in geometrically thin disks with rotation curves dominated by a central object, but with substantial contributions from magnetic pressure and tension. The treatment is general, but the application is to the circumstellar disks that arise during the gravitational collapse phase of star formation. We find the dispersion relation for spiral density waves in these generalized disks and derive the stability criterion for axisymmetric $(m=0)$ disturbances (the analog of the Toomre parameter $Q_T$) for any radial distribution of the mass-to-flux ratio $\lambda$. The magnetic effects work in two opposing directions: on one hand, magnetic tension and pressure stabilize the disk against gravitational collapse and fragmentation; on the other hand, they also lower the rotation rate making the disk more unstable. For disks around young stars the first effect generally dominates, so that magnetic fields allow disks to be stable for higher surface densities and larger total masses. These results indicate that magnetic fields act to suppress the formation of giant planets through gravitational instability. Finally, even if gravitational instability can form a secondary body, it must lose an enormous amount of magnetic flux in order to become a planet; this latter requirement represents an additional constraint for planet formation via gravitational instability and places a lower limit on the electrical resistivity.

In the framework of the study of the Galactic metallicity gradient and its time evolution, we present new high-resolution spectroscopic observations obtained with FLAMES and the fiber link to UVES at VLT of three open clusters (OCs) located within $\sim$7~kpc from the Galactic Center (GC): NGC~6192, NGC~6404, NGC~6583. We also present new orbit determination for all OCs with Galactocentric distances (R$_{\rm{GC}}) \leq$8~kpc and metallicity from high-resolution spectroscopy. We aim to investigate the slope of the inner disk metallicity gradient as traced by OCs and at discussing its implication on the chemical evolution of our Galaxy. We have derived memberships of a group of evolved stars for each clusters, obtaining a sample of 4, 4, and 2 member stars in NGC~6192, NGC~6404, and NGC~6583, respectively. Using standard LTE analysis we derived stellar parameters and abundance ratios for the iron-peak elements Fe, Ni, Cr, and for the $\alpha$-elements Al, Mg, Si, Ti, Ca. We calculated the orbits of the OCs currently located within 8~kpc from the GC, and discuss their implication on the present-time radial location. {The average metallicities of the three clusters are all oversolar: [Fe/H]= $+0.12\pm0.04$ (NGC~6192), $+0.11\pm0.04$ (NGC 6404), $+0.37\pm0.03$ (NGC 6583). They are in qualitative agreement with their Galactocentric distances, being all internal OCs, and thus expected to be metal richer than the solar neighborhood. The abundance ratios of the other elements over iron [X/Fe] are consistent with solar values. The clusters we have analysed, together with other OC and Cepheid data, confirm a steep gradient in the inner disk, a signature of an evolutionary rate different than in the outer disk.

We analyze the chemical abundances of planetary nebulae and HII regions in the M81 disk for insight on galactic evolution, and compare it with that of other galaxies, including the Milky Way. We acquired Hectospec/MMT spectra of 39 PNe and 20 HII regions, with 33 spectra viable for temperature and abundance analysis. Our PN observations represent the first PN spectra in M81 ever published, while several HII region spectra have been published before, although without a direct electron temperature determination. We determine elemental abundances of helium, nitrogen, oxygen, neon, sulfur, and argon in PNe and HII regions, and determine their averages and radial gradients. The average O/H ratio of PNe compared to that of the HII regions indicates a general oxygen enrichment in M81 in the last ~10 Gyr. The PN metallicity gradient in the disk of M81 is -0.055+-0.02 dex/kpc. Neon and sulfur in PNe have a radial distribution similar to that of oxygen, with similar gradient slopes. If we combine our HII sample with the one in the literature we find a possible mild evolution of the gradient slope, with results consistent with gradient steepening with time. Additional spectroscopy is needed to confirm this trend. There are no Type I PNe in our M81 sample, consistently with the observation of only the brightest bins of the PNLF, the galaxy metallicity, and the evolution of post-AGB shells. Both the young and the old populations of M81 disclose shallow but detectable negative radial metallicity gradient, which could be slightly steeper for the young population, thus not excluding a mild gradients steepening with the time since galaxy formation. During its evolution M81 has been producing oxygen; its total oxygen enrichment exceeds that of other nearby galaxies.

Nearby galaxies are ideal places to study in detail metallicity gradients and their time evolution. We consider chemical abundances of a new sample of \hii\ regions complemented with previous literature data-sets. We compare \hii\ region and PN abundances obtained with a common set of observations taken at MMT. With an updated theoretical model, we follow the time evolution of the baryonic components and chemical abundances in the disk of M33, assuming that the galaxy is accreting gas from an external reservoir. Supported by a uniform sample of nebular spectroscopic observations, we conclude that: {\em i}) the metallicity distribution in M33 is very complex, showing a central depression in metallicity probably due to observational bias; {\em ii}) the metallicity gradient in the disk of M33 has a slope of -0.037$\pm$ 0.009 dex kpc$^{-1}$ in the whole radial range up to $\sim$8 kpc, and -0.044$\pm$ 0.009 dex kpc$^{-1}$ excluding the central kpc; {\em iii}) there is a small evolution of the slope with time from the epoch of PN progenitor formation to the present-time.}

The single stable isotope of beryllium is a pure product of cosmic-ray spallation in the ISM. Assuming that the cosmic-rays are globally transported across the Galaxy, the beryllium production should be a widespread process and its abundance should be roughly homogeneous in the early-Galaxy at a given time. Thus, it could be useful as a tracer of time. In an investigation of the use of Be as a cosmochronometer and of its evolution in the Galaxy, we found evidence that in a log(Be/H) vs. [alpha/Fe] diagram the halo stars separate into two components. One is consistent with predictions of evolutionary models while the other is chemically indistinguishable from the thick-disk stars. This is interpreted as a difference in the star formation history of the two components and suggests that the local halo is not a single uniform population where a clear age-metallicity relation can be defined. We also found evidence that the star formation rate was lower in the outer regions of the thick disk, pointing towards an inside-out formation.

Beryllium stellar abundances were suggested to be a good tracer of time in the early Galaxy. In an investigation of its use as a cosmochronometer, using a large sample of local halo and thick-disk dwarfs, evidence was found that in a log(Be/H) vs. [alpha/Fe] diagram the halo stars separate into two components. One is consistent with predictions of evolutionary models while the other is chemically indistinguishable from the thick-disk stars. This is interpreted as a difference in the star formation history of the two components and suggests that the local halo is not a single uniform population where a clear age-metallicity relation can be defined.

We study the abundance of CCH in prestellar cores both because of its role in the chemistry and because it is a potential probe of the magnetic field. We also consider the non-LTE behaviour of the N=1-0 and N=2-1 transitions of CCH and improve current estimates of the spectroscopic constants of CCH. We used the IRAM 30m radiotelescope to map the N=1-0 and N=2-1 transitions of CCH towards the prestellar cores L1498 and CB246. Towards CB246, we also mapped the 1.3 mm dust emission, the J=1-0 transition of N2H+ and the J=2-1 transition of C18O. We used a Monte Carlo radiative transfer program to analyse the CCH observations of L1498. We derived the distribution of CCH column densities and compared with the H2 column densities inferred from dust emission. We find that while non-LTE intensity ratios of different components of the N=1-0 and N=2-1 lines are present, they are of minor importance and do not impede CCH column density determinations based upon LTE analysis. Moreover, the comparison of our Monte-Carlo calculations with observations suggest that the non-LTE deviations can be qualitatively understood. For L1498, our observations in conjunction with the Monte Carlo code imply a CCH depletion hole of radius 9 x 10^{16} cm similar to that found for other C-containing species. We briefly discuss the significance of the observed CCH abundance distribution. Finally, we used our observations to provide improved estimates for the rest frequencies of all six components of the CCH(1-0) line and seven components of CCH(2-1). Based on these results, we compute improved spectroscopic constants for CCH. We also give a brief discussion of the prospects for measuring magnetic field strengths using CCH.

We explore the effects of magnetic energy dissipation on the formation of the first stars. For this purpose, we follow the evolution of primordial chemistry in the presence of magnetic fields in the post-recombination universe until the formation of the first virialized halos. From the point of virialization, we follow the protostellar collapse up to densities of $\sim10^{12}$ cm$^{-3}$ in a one-zone model. In the intergalactic medium (IGM), comoving field strengths of $\gtrsim0.1$ nG lead to Jeans masses of $10^8 M_\odot$ or more and thus delay gravitational collapse in the first halos until they are sufficiently massive. During protostellar collapse, we find that the temperature minimum at densities of $\sim10^3$ cm$^{-3}$ does not change significantly, such that the characteristic mass scale for fragmentation is not affected. However, we find a significant temperature increase at higher densities for comoving field strengths of $\gtrsim0.1$ nG. This may delay gravitational collapse, in particular at densities of $\sim10^9$ cm$^{-3}$, where the proton abundance drops rapidly and the main contribution to the ambipolar diffusion resistivity is due to collisions with Li$^+$. After the formation of the protostar, the increased gas temperatures may enhance the protostellar accretion rate. Our model confirms that initial weak magnetic fields may be amplified considerably during gravitational collapse and become dynamically relevant. For instance, a comoving field strength above $10^{-5}$ nG will be amplified above the critical value for the onset of jets which can magnetize the IGM.

Low-energy cosmic rays are a fundamental source of ionization for molecular clouds, influencing their chemical, thermal and dynamical evolution. The purpose of this work is to explore the possibility that a low-energy component of cosmic-rays, not directly measurable from the Earth, can account for the discrepancy between the ionization rate measured in diffuse and dense interstellar clouds. We collect the most recent experimental and theoretical data on the cross sections for the production of H2+ and He+ by electron and proton impact, and we discuss the available constraints on the cosmic-ray fluxes in the local interstellar medium. Starting from different extrapolations at low energies of the demodulated cosmic-ray proton and electron spectra, we compute the propagated spectra in molecular clouds in the continuous slowing-down approximation taking into account all the relevant energy loss processes. The theoretical value of the cosmic-ray ionization rate as a function of the column density of traversed matter is in agreement with the observational data only if either the flux of cosmic-ray electrons or of protons increases at low energies. The most successful models are characterized by a significant (or even dominant) contribution of the electron component to the ionization rate, in agreement with previous suggestions. However, the large spread of cosmic-ray ionization rates inferred from chemical models of molecular cloud cores remains to be explained. Available data combined with simple propagation models support the existence of a low-energy component (below about 100 MeV) of cosmic-ray electrons or protons responsible for the ionization of molecular cloud cores and dense protostellar envelopes.

[abridged] Beryllium is a pure product of cosmic ray spallation. This implies a relatively simple evolution in time of the beryllium abundance and suggests its use as a time-like observable. We study the evolution of Be in the early Galaxy and its dependence on kinematic and orbital parameters. We investigate the formation of the halo and the thick disk of the Galaxy and the use of Be as a cosmochronometer. Beryllium abundances are determined from high resolution, high signal to noise UVES spectra with spectrum synthesis in the largest sample of halo and thick disk stars analyzed to date. We present our observational results in various diagrams. 1) In a log(Be/H) vs [Fe/H] diagram we find a marginal statistical detection of a real scatter, above what expected from measurement errors, with a larger scatter among halo stars. The detection of the scatter is further supported by the existence of pairs of stars with identical atmospheric parameters and different Be abundances. 2) In an log(Be/H) vs [alpha/Fe] diagram, the halo stars separate into two components; one is consistent with predictions of evolutionary models, while the other has too high alpha and Be abundances and is chemically indistinguishable from thick disk stars. This suggests that the halo is not a single uniform population where a clear age-metallicity relation can be defined. 3) In diagrams of Rmin vs [alpha/Fe] and log(Be/H) the thick disk stars show a possible decrease of [alpha/Fe] with Rmin, whereas no dependence of Be with Rmin is seen. This anticorrelation suggests that the star formation rate was lower in the outer regions of the thick disc, pointing towards an inside-out formation. The lack of correlation for Be indicates that it is insensitive to the local conditions of star formation.

We report the discovery of a trend of increasing barium abundance with decreasing age for a large sample of Galactic open clusters. The observed pattern of [Ba/Fe] vs. age can be reproduced with a Galactic chemical evolution model only assuming a higher Ba yield from the $s$-process in low-mass stars than the average one suggested by parametrized models of neutron-capture nucleosynthesis. We show that this is possible in a scenario where the efficiency of the extra-mixing processes producing the neutron source $^{13}$C is anti-correlated with the initial mass, with a larger efficiency for lower masses. This is similar to the known trend of extended mixing episodes acting in H-rich layers and might suggest a common physical mechanism.

Open clusters offer a unique possibility to study the time evolution of the radial metallicity gradients of several elements in our Galaxy, because they span large intervals in age and Galactocentric distance, and both quantities can be more accurately derived than for field stars. We re-address the issue of the Galactic metallicity gradient and its time evolution by comparing the empirical gradients traced by a sample of 45 open clusters with a chemical evolution model of the Galaxy. At variance with previous similar studies, we have collected from the literature only abundances derived from high--resolution spectra. The clusters have distances $7 < RGC<22$ kpc and ages from $\sim 30$ Myr to 11 Gyr. We also consider the $\alpha$-elements Si, Ca, Ti, and the iron-peak elements Cr and Ni. The data for iron-peak and $\alpha$-elements indicate a steep metallicity gradient for R_GC<12$ kpc and a plateau at larger radii. The time evolution of the metallicity distribution is characterized by a uniform increase of the metallicity at all radii, preserving the shape of the gradient, with marginal evidence for a flattening of the gradient with time in the radial range 7-12 kpc. Our model is able to reproduce the main features of the metallicity gradient and its evolution with an infall law exponentially decreasing with radius and with a collapse time scale of the order of 8 Gyr at the solar radius. This results in a rapid collapse in the inner regions, i.e. $R_{\rm GC}\lesssim 12$ kpc (that we associate with an early phase of disk formation from the collapse of the halo) and in a slow inflow of material per unit area in the outer regions at a constant rate with time.

Magnetic fields are believed to play a crucial role in the process of star formation. We compare high-angular resolution observations of the submillimeter polarized emission of NGC 1333 IRAS 4A, tracing the magnetic field around a low-mass protostar, with models of the collapse of magnetized molecular cloud cores. Assuming a uniform dust alignment efficiency, we computed the Stokes parameters and synthetic polarization maps from the model density and magnetic field distribution by integrations along the line-of-sight and convolution with the interferometric response. The synthetic maps are in good agreement with the data. The best-fitting models were obtained for a protostellar mass of 0.8 solar masses, of age 9e4 yr, formed in a cloud with an initial mass-to-flux ratio ~2 times the critical value. The magnetic field morphology in NGC 1333 IRAS 4A is consistent with the standard theoretical scenario for the formation of solar-type stars, where well-ordered, large-scale, rather than turbulent, magnetic fields control the evolution and collapse of the molecular cloud cores from which stars form.

Previous works have demonstrated that the generation of secondary CMB anisotropies due to the molecular optical depth is likely too small to be observed. In this paper, we examine additional ways in which primordial chemistry and the dark ages might influence the CMB. We present a detailed and updated chemical network and give an overview of the interactions of molecules with the CMB. We consider the optical depth due to line absorption, photoionization, photodissociation and free-free processes, and estimate the resulting changes in the CMB temperature and its power spectrum. The most promising results are obtained for the negative hydrogen ion \HM and the \HeHII molecule. The free-free process of \HM yields a relative change in the CMB temperature of up to $2\times10^{-11}$, and leads to a frequency-dependent change in the power spectrum of the order $10^{-7}$ at 30 GHz. With a change of the order $10^{-10}$ in the power spectrum, our result for the bound-free process of \HM is significantly below a previous suggestion. \HeHII efficiently scatters CMB photons and smears out primordial fluctuations, leading to a change in the power spectrum of the order $10^{-8}$. Improvements in the accuracy of future CMB experiments may thus help to constrain and finally detect these interesting signals from the dark ages of the universe.

Strong magnetization makes the disks surrounding young stellar objects rotate at rates that are too sub-Keplerian to enable the thermal launching of disk winds from their surfaces unless the rate of gas diffusion across field lines is dynamically fast. This under-appreciated implication of disk magnetization poses a considerable challenge for disk-wind theory.

Interstellar magnetic fields influence all stages of the process of star formation, from the collapse of molecular cloud cores to the formation of protostellar jets. This requires us to have a full understanding of the physical properties of magnetized plasmas of different degrees of ionization for a wide range of densities and temperatures. We derive general equations governing the magneto-hydrodynamic evolution of a three-fluid medium of arbitrary ionization, also including the possibility of charged dust grains as the main charge carriers. In a companion paper (Pinto & Galli 2007), we complement this analysis computing accurate expressions of the collisional coupling coefficients. Over spatial and temporal scales larger than the so-called large-scale plasma limit and the collision-dominated plasma limit, and for non-relativistic fluid speeds, we obtain an advection-diffusion for the magnetic field. We derive the general expressions for the resistivities, the diffusion time scales and the heating rates in a three-fluid medium and we use them to estimate the evolution of the magnetic field in molecular clouds and protostellar jets. Collisions between charged particles significantly increase the value of the Ohmic resistivity during the process of cloud collapse, affecting in particular the decoupling of matter and magnetic field and enhancing the rate of energy dissipation. The Hall resistivity can take larger values than previously found when the negative charge is mostly carried by dust grains. In weakly-or mildy-ionized protostellar jets, ambipolar diffusion is found to occur on a time scale comparable to the dynamical time scale, limiting the validity of steady-state and nondissipative models to study the jet's structure.

The charged component of the insterstellar medium consists of atomic and molecular ions, electrons, and charged dust grains, coupled to the local Galactic magnetic field. Collisions between neutral particles (mostly atomic or molecular hydrogen) and charged species, and between the charged species themselves, affect the magnetohydrodynamical behaviour of the medium and the dissipation of electric currents. The friction force due to elastic collisions between particles of different species in the multi-component interstellar plasma is a nonlinear function of the temperature of each species and the Mach number of the relative drift velocity. The aim of this paper is to provide an accurate and, as far as possible, complete set of momentum transfer rate coefficients for magnetohydrodynamical studies of the interstellar medium. Momentum transfer rates are derived from available experimental data and theoretical calculations of cross sections within the classic approach developed by Boltzmann and Langevin for a wide range of values of the temperature and the drift velocity. Accurate numerical values for momentum transfer rates are obtained and fitted to simple analytical formulae expressing the dependence of the results on the gas temperature and the relative drift velocity. The often used polarization approximation is in satisfactory agreement with our results only for collisions between H2and molecular ions (HCO+, H3+). For other kinds of collisions, the polarization approximation fails by large factors, and must be replaced by more accurate expressions.

As part of a large survey of halo and thick disc stars, we found one halo star, HD 106038, exceptionally overabundant in beryllium. In spite of its low metallicity, [Fe/H] = -1.26, the star has log(Be/H) = -10.60, which is similar to the solar meteoritic abundance, log(Be/H) = -10.58. This abundance is more than ten times higher the abundance of stars with similar metallicity and cannot be explained by models of chemical evolution of the Galaxy that include the standard theory of cosmic-ray spallation. No other halo star exhibiting such a beryllium overabundance is known. In addition, overabundances of Li, Si, Ni, Y, and Ba are also observed. We suggest that all these chemical peculiarities, but the Ba abundance, can be simultaneously explained if the star was formed in the vicinity of a hypernova.

We consider the formation of OB associations from two perspectives: (a) the fractional gas consumption in star formation,epsilon, per dynamical time scale t_dyn in a galaxy, and (b) the origin of the so-called Kennicutt-Schmidt law that the rate of star formation per unit area is proportional to a power, alpha, of the surface density in HI and H_2 gas when certain thresholds are crossed. The empirical findings that epsilon is approximately 0.01 and alpha is approximately 1.4 or 1.5 have simple explanations if the rate of star formation is magnetically regulated. An empirical test of the ideas resides in an analysis of why giant OB associations are ``strung out like pearls along the arms" of spiral galaxies.

We consider the accretion process in a disk with magnetic fields that are dragged in from the interstellar medium by gravitational collapse. Two diffusive processes are at work in the system: (1) "viscous" torques exerted by turbulent and magnetic stresses, and (2) "resistive" redistribution of mass with respect to the magnetic flux arising from the imperfect conduction of current. In steady state, self-consistency between the two rates of drift requires that a relationship exists between the coefficients of turbulent viscosity and turbulent resistivity. Ignoring any interactions with a stellar magnetosphere, we solve the steady-state equations for a magnetized disk under the gravitational attraction of a mass point and threaded by an amount of magnetic flux consistent with calculations of magnetized gravitational collapse in star formation. Our model mean-field equations have an exact analytical solution that corresponds to magnetically diluted Keplerian rotation about the central mass point. The solution yields the strength of the magnetic field and the surface density as functions of radial position in the disk and their connection with the departure from pure Keplerian rotation in representative cases. We compare the predictions of the theory with the available observations concerning T Tauri stars, FU Orionis stars, and low- and high-mass protostars. Finally, we speculate on the physical causes for high and low states of the accretion disks that surround young stellar objects. One of the more important results of this study is the physical derivation of analytic expressions for the turbulent viscosity and turbulent resistivity.

The evolution of radial gradients of metallicity in disk galaxies and its relation with the disk formation are not well understood. Theoretical models of galactic chemical evolution make contrasting predictions about the time evolution of metallicity gradients. To test chemical evolution models and trace the star formation and accretion history of low luminosity disk galaxies we focus on the Local Group galaxy M33. We analyze O/H and S/H abundances in planetary nebulae, H{\sc ii} regions, and young stars, together with known [Fe/H] abundances in the old stellar population of M33. With a theoretical model, we follow the time evolution of gas (diffuse and condensed in clouds), stars, and chemical abundances in the disk of M33, assuming that the galaxy is accreting gas from an external reservoir. Our model is able to reproduce the available observational constraints on the distribution of gas and stars in M33 and to predict the time evolution of several chemical abundances. In particular, we find that a model characterized by a continuous infall of gas on the disk, at a rate of $\dot M_{\rm inf}\approx 1$ $M_\odot$ yr$^{-1}$, almost constant with time, can also account for the relatively high rate of star formation and for the shallow chemical gradients. Supported by a large sample of high resolution observations for this nearby galaxy, we conclude that the metallicity in the disk of M33 has increased with time at all radii, with a continuous flattening of the gradient over the last $\sim 8$ Gyr.

Aims: To measure the beryllium abundance in two TO stars of the Globular Cluster NGC 6752, one oxygen rich and sodium poor, the other presumably oxygen poor and sodium rich. Be abundances in these stars are used to put on firmer grounds the hypothesis of Be as cosmochronometer and to investigate the formation of Globular Clusters. Method:We present near UV spectra with resolution R$\sim 45000$ obtained with the UVES spectrograph on the 8.2m VLT Kueyen telescope, analysed with spectrum synthesis based on plane parallel LTE model atmospheres. Results:Be is detected in the O rich star with log(Be/H)=-12.04 $\pm$0.15, while Be is not detected in the other star for which we obtain the upper limit log(Be/H)$<$-12.2. A large difference in nitrogen abundance (1.6 dex) is found between the two stars. Conclusions:The Be measurement is compatible with what found in field stars with the same [Fe/H] and [O/H]. The 'Be age' of the cluster is found to be 13.3 Gyrs, in excellent agreement with the results from main sequence fitting and stellar evolution. The presence of Be confirms the results previously obtained for the cluster NGC 6397 and supports the hypothesis that Be can be used as a clock for the early formation of the Galaxy. Since only an upper limit is found for the star with low oxygen abundance, we cannot decide between competing scenarios of Globular Cluster formation, but we can exclude that 'polluted' stars are substantially younger than 'unpolluted' ones. We stress that the Be test might be the only measurement capable of distinguishing between these scenarios.

We formulate the problem of magnetic field dissipation during the accretion phase of low-mass star formation, and we carry out the first step of an iterative solution procedure by assuming that the gas is in free-fall along radial field lines. This so-called ``kinematic approximation'' ignores the back reaction of the Lorentz force on the accretion flow. In quasi steady-state, and assuming the resistivity coefficient to be spatially uniform, the problem is analytically soluble in terms of Legendre's polynomials and confluent hypergeometric functions. The dissipation of the magnetic field occurs inside a region of radius inversely proportional to the mass of the central star (the ``Ohm radius''), where the magnetic field becomes asymptotically straight and uniform. In our solution, the magnetic flux problem of star formation is avoided because the magnetic flux dragged in the accreting protostar is always zero. Our results imply that the effective resistivity of the infalling gas must be higher by several orders of magnitude than the microscopic electric resistivity, to avoid conflict with measurements of paleomagnetism in meteorites and with the observed luminosity of regions of low-mass star formation.

We study the self-similar collapse of an isothermal magnetized rotating cloud in the ideal magnetohydrodynamic (MHD) regime. In the limit of small distance from the accreting protostar we find an analytic solution that corresponds to free-fall onto a central mass point. The density distribution is not spherically symmetric but depends on the mass loading of magnetic field lines, which can be obtained by matching our inner solution to an outer collapse solution previously computed by Allen, Shu & Li. The concentration of magnetic field trapped by the central mass point under field-freezing, independent on the details of the starting state, creates a split monopole configuration where the magnetic field strength increases as the inverse square of the distance from the center. Under such conditions, the inflow eventually becomes subalfvenic and the outward transfer of angular momentum by magnetic braking very efficient, thus preventing the formation of a centrifugally supported disk. Instead, the azimuthal velocity of the infalling gas decreases to zero at the center, and the gas spirals into the star. Therefore, the dissipation of dynamically important levels of magnetic field is a fundamental requisite for the formation of protoplanetary disks around young stars.

Disks are a natural outcome of the star formation process in which they play a crucial role. Luminous, massive stars of spectral type earlier than B4 are likely to be those that benefit most from the existence of accretion disks, which may significantly reduce the effect of radiation pressure on the accreting material. The scope of the present contribution is to review the current knowledge about disks in young high-mass (proto)stars and discuss their implications. The issues of disk stability and lifetime are also discussed. We conclude that for protostars of less than ~20 solar masses, disks with mass comparable to that of the central star are common. Above this limit the situation is unclear and there are no good examples of proto O4-O8 stars surrounded by accretion disks: in these objects only huge, massive, toroidal, non-equilibrium rotating structures are seen. It is clear on the other hand that the observed disks in stars of 10-20 solar masses are likely to be unstable and with short lifetimes.

We investigate the evolution of the star formation rate in the early Galaxy using beryllium and oxygen abundances in metal poor stars. Specifically, we show that stars belonging to two previously identified kinematical classes (the so-called ``accretion'' and ``dissipative'' populations) are neatly separated in the [O/Fe] vs. Log(Be/H) diagram. The dissipative population follows the predictions of our model of Galactic evolution for the thick disk component, suggesting that the formation of this stellar population occurred on a timescale significantly longer (by a factor ~5-10) than the accretion component. The latter shows a large scatter in the [O/Fe] vs. Log(Be/H) diagram, probably resulting from the inhomogeneous enrichment in oxygen and iron of the protogalactic gas. Despite the limitation of the sample, the data suggest that the combined use of products of spallation reactions (like beryllium) and elemental ratios of stellar nucleosynthesis products (like [O/Fe]) can constrain theoretical models for the formation and early evolution of our Galaxy.

(Abridged) We develop an analytical spectral method to solve the equations of equilibrium for a self-gravitating, magnetized fluid body, under the only hypotheses that (a) the equation of state is isothermal, (b) the configuration is scale-free, and (c) the body is electrically neutral. All physical variables are represented as series of scalar and vector spherical harmonics of degree l and order m, and the equilibrium equations are reduced to a set of coupled quadratic algebraic equations for the expansion coefficients of the density and the magnetic vector potential. The method is general, and allows to recover previously known hydrostatic and magnetostatic solutions possessing axial symmetry. A linear perturbation analysis of the equations in spectral form show that these basic axisymmetric states, considered as a continuos sequence with the relative amount of magnetic support as control parameter, have in general no neighboring nonaxisymmetric equilibria. This result lends credence to a conjecture originally made by H. Grad and extends early results obtained by E. Parker to the case of self-gravitating magnetized bodies. The only allowed bifurcations of this sequence of axisymmetric equilibria are represented by distortions with dipole-like angular dependence (l=1) that can be continued into the nonlinear regime. These new configurations are either (i) azimuthally asymmetric (m=+-1) or (ii) azimuthally symmetric but without reflection symmetry with respect to the equatorial plane (m=0). To the extent that interstellar clouds can be represented as isolated magnetostatic equilibria, the results of this study suggest that the observed triaxial shapes of molecular cloud cores can be interpreted in terms of weakly damped Alfven oscillations about an equilibrium state.

This paper reports simulations allowing for stochastic accretion and mass loss within closed and open systems modeled using a previously developed multi-population, multi-zone (halo, thick disk, thin disk) treatment. The star formation rate is computed as a function of time directly from the model equations and all chemical evolution is followed without instantaneous recycling. Several types of simulations are presented here: (1) a closed system with bursty mass loss from the halo to the thick disk, and from the thick to the thin disk, in separate events to the thin disk; (2) open systems with random environmental (extragalactic) accretion, e.g. by infall of high velocity clouds directly to the thin disk; (3) schematic open system single and multiple collision events and intracluster stripping. For the open models, the mass of the Galaxy has been explicitly tracked with time. We present the evolution of the star formation rate, metallicity histories, and concentrate on the light elements. We find a wide range of possible outcomes, including an explanation for variations in the Galactic D/H ratio, and highlight the problems for uniquely reconstructing star forming histories from contemporary abundance measurements.

We recall the emergence of the ``3He problem'', its currently accepted solution, and we summarize the presently available constraints on models of stellar nucleosynthesis and studies of Galactic chemical evolution from measurements of the abundance of 3He in the Galaxy.

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 the first detection of beryllium in two turnoff stars of the old, metal-poor globular cluster NGC 6397. The beryllium lines are clearly detected and we determine a mean beryllium abundance of log(Be/H)=-12.35 +/- 0.2. The beryllium abundance is very similar to that of field stars of similar Fe content. We interpret the beryllium abundance observed as the result of primary spallation of cosmic rays acting on a Galactic scale, showing that beryllium can be used as a powerful cosmochronometer for the first stellar generations. With this method, we estimate that the cluster formed 0.2-0.3 Gyr after the onset of star formation in the Galaxy, in excellent agreement with the age derived from main sequence fitting. From the same spectra we also find low O (noticeably different for the two stars) and high N abundances, suggesting that the original gas was enriched in CNO processed material. Our beryllium results, together with the N, O, and Li abundances, provide insights on the formation of this globular cluster, showing that any CNO processing of the gas must have occurred in the protocluster cloud before the formation of the stars we observe now. We encounter, however, difficulties in giving a fully consistent picture of the cluster formation, able to explain the complex overall abundance pattern.

We develop a formalism wherein the solution of the equations of equilibrium for a self-gravitating, magnetized interstellar cloud may be accomplished analytically, under the only hypotheses that (a) the cloud is scale-free, and (b) is electrically neutral. All variables are represented as series of scalar and vector spherical harmonics, and the equilibrium equations are reduced to a set of coupled algebraic equations for the expansion coefficients of the density and the magnetic vector potential. Previously known axisymmetric solutions are recovered and new non-axisymmetric solutions are found and discussed as models for pre-protostellar molecular cloud cores.

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.

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.

(Abridged) We consider the expansion of an initially self-gravitating, static, singular, isothermal cloud core. For t>0, the gas is ionized and heated to a higher uniform temperature by the formation of a luminous, but massless, star in its center. The approximation that the mass and gravity of the central star is negligible for the subsequent motion of the HII region holds for distances much greater than about 100 AU and for the massive cloud cores that give rise to high-mass stars. If the initial ionization and heating is approximated to occur instantaneously at t=0, then the subsequent flow (for r >> 100 AU) caused by the resulting imbalance between self-gravity and thermal pressure is self-similar. Because of the steep density profile, pressure gradients produce a shock front that travels into the cloud, accelerating the gas to supersonic velocities in what has been called the ``champagne phase.'' We then study the self-similar solutions of the expansion of HII regions embedded in molecular clouds characterized by more general power-law density distributions with exponent 3/2 < n < 3. In these cases, the shock velocity is an increasing function of the exponent n, and diverges as n tends to 3. We show that this happens because the model includes an origin, where the pressure driving the shock diverges because the enclosed heated mass is infinite. Our results imply that the continued photoevaporation of massive reservoirs of neutral gas (e.g., surrounding disks and/or globules) nearby to the embedded ionizing source is required in order to maintain over a significant timescale the emission measure observed in champagne flows.

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.

We review and update some aspects of deuterium chemistry in the post-recombination Universe with particular emphasis on the formation and destruction of HD. We examine in detail the available theoretical and experimental data for the leading reactions of deuterium chemistry and we highlight the areas where improvements in the determination of rate coefficients are necessary to reduce the remaining uncertainties. We discuss the cooling properties of HD and the modifications to the standard cooling function introduced by the presence of the cosmological radiation field. Finally, we consider the effects of deuterium chemistry on the dynamical collapse of primordial clouds in a simple ``top-hat'' scenario, and we speculate on the minimum mass a cloud must have in order to be able to cool in a Hubble time.

We present high resolution HST ultraviolet spectroscopy of NGC3242, the only planetary nebula with a measured abundance of 3He. The Space Telescope Imaging Spectrograph (STIS) has been used to observe the CIII multiplet near 1908 A. The presence of 12C and 13C lines at these wavelengths allows a direct estimate of the carbon isotopic ratio in the ionized gas of the nebula. We have detected the 12C doublet and obtained an upper limit on the 13C 1909.6 A line, resulting in a carbon isotopic ratio 12C/13C > 38, in agreement with standard stellar models. The lack of the 13C line and the presence of the 3He line in the spectrum indicate that the progenitor star did not undergo a phase of deep mixing during the last stages of its evolution. The significance of this result for studies of stellar nucleosynthesis and Galactic chemical evolution is discussed.

In this paper we study the production of Li, Be and B nuclei by Galactic cosmic ray spallation processes. We include three kinds of processes: (i) spallation by light cosmic rays impinging on interstellar CNO nuclei (direct processes); (ii) spallation by CNO cosmic ray nuclei impinging on interstellar p and 4He (inverse processes); and (iii) alpha-alpha fusion reactions. The latter dominate the production of 6Li and 7Li. We calculate production rates for a closed-box Galactic model, verifying the quadratic dependence of the Be and B abundances for low values of Z. These are quite general results and are known to disagree with observations. We then show that the multi-zone multi-population model we used previously for other aspects of Galactic evolution produces quite good agreement with the linear trend observed at low metallicities without fine tuning. We argue that reported discrepancies between theory and observations do not represent a nucleosynthetic problem, but instead are the consequences of inaccurate treatments of Galactic evolution.

We have computed the dust temperature distribution to be expected in a pre-protostellar core in the phase prior to the onset of gravitational instability. We have done this in the approximation that the heating of the dust grains is solely due to the attenuated external radiation field and that the core is optically thin to its own radiation. This permits us to consider non spherically symmetric geometries. We predict the intensity distributions of our model cores at millimeter and sub-millimeter wavelengths and compare with observations of the well studied object L1544. We have also developed an analytical approximation for the temperature at the center of spherically symmetric cores and we compare this with the numerical calculations. Our results show (in agreement with Evans et al. 2001) that the temperatures in the nuclei of cores of high visual extinction (> 30 magnitudes) are reduced to values of below ~8 K or roughly half of the surface temperature. This has the consequence that maps at wavelengths shortward of 1.3 mm see predominantly the low density exterior of pre-protostellar cores. It is extremely difficult to deduce the true density distribution from such maps alone. We have computed the intensity distribution expected on the basis of the models of Ciolek & Basu (2000) and compared with the observations of L1544. The agreement is good with a preference for higher inclinations (37 degrees instead of 16) than that adopted by Ciolek & Basu (2000). We find that a simple extension of the analytic approximation allows a reasonably accurate calculation of the dust temperature as a function of radius in cores with density distributions approximating those expected for Bonnor-Ebert spheres and suggest that this may be a useful tool for future calculations of the gas temperature in such cores.

In this paper we study the evolution of 7Li in the Galaxy considering the contributions of various stellar sources: type II supernovae, novae, red giant stars, and asymptotic giant branch (AGB) stars. We present new results for the production of 7Li in AGB stars via the hot bottom burning process, based on stellar evolutionary models by Frost (1997). In the light of recent observations of dense circumstellar shells around evolved stars in the Galaxy and in the Magellanic Clouds, we also consider the impact of a very high mass-loss rate episode (superwind) before the evolution off the AGB phase on the 7Li enrichment in the interstellar medium. We compare the Galactic evolution of 7Li obtained with these new 7Li yields (complemented with a critical re-analysis of the role of supernovae, novae and giant stars) with a selected compilation of spectroscopic observations including halo and disk field stars as well as young stellar clusters. We conclude that even allowing for the large uncertainties in the theoretical calculation of mass-loss rates at the end of the AGB phase, the superwind phase has a significant effect on the 7Li enrichment of the Galaxy.

We review the difficulties of the classical fission and fragmentation hypotheses for the formation of binary and multiple stars. A crucial missing ingredient in previous theoretical studies is the inclusion of dynamically important levels of magnetic fields. As a minimal model for a candidate presursor to the formation of binary and multiple stars, we therefore formulate and solve the problem of the equilibria of isopedically magnetized, singular isothermal disks, without the assumption of axial symmetry. Considerable analytical progress can be made if we restrict our attention to models that are scale-free, i.e., that have surface densities that vary inversely with distance from the rotation axis of the system. In agreement with earlier analysis by Syer and Tremaine, we find that lopsided (M=1) configurations exist at any dimensionless rotation rate, including zero. Multiple-lobed (M = 2, 3, 4, ...) configurations bifurcate from an underlying axisymmetric sequence at progressively higher dimensionless rates of rotation, but such nonaxisymmetric sequences always terminate in shockwaves before they have a chance to fission into M=2, 3, 4, ... separate bodies. On the basis of our experience in this paper, and the preceding Paper I, we advance the hypothesis that binary and multiple star-formation from smooth (i.e., not highly turbulent) starting states that are supercritical but in unstable mechanical balance requires the rapid (i.e., dynamical) loss of magnetic flux at some stage of the ensuing gravitational collapse.

A crucial missing ingredient in previous theoretical studies of fragmentation is the inclusion of dynamically important levels of magnetic fields. As a minimal model for a candidate presursor to the formation of binary and multiple stars, we therefore consider the equilibrium configuration of isopedically magnetized, scale-free, singular isothermal disks, without the assumption of axial symmetry. We find that lopsided (M=1) configurations exist at any dimensionless rotation rate, including zero. Multiple-lobed (M = 2, 3, 4, ...) configurations bifurcate from an underlying axisymmetric sequence at progressively higher dimensionless rates of rotation, but such nonaxisymmetric sequences always terminate in shockwaves before they have a chance to fission into separate bodies. We advance the hypothesis that binary and multiple star-formation from smooth (i.e., not highly turbulent) starting states that are supercritical but in unstable mechanical balance requires the rapid (i.e., dynamical) loss of magnetic flux at some stage of the ensuing gravitational collapse.

We present a stochastic model based on the Monte Carlo technique to study the inhomogeneous chemical evolution of the Galactic halo. In particular, we consider the local enrichment and subsequent mixing of the interstellar gas resulting from bursts of star formation, following explicitly the fragmentation and coalescence of interstellar gas clouds. The model takes into account the mixing of halo gas through cloud-cloud collisions and the delayed mixing of supernova ejecta into the interstellar medium. The consequences of the infall of halo gas onto the Galactic disk are also discussed. The mass spectrum of clouds, the age-metallicity relation, and the G-dwarf distribution are investigated at different times. We analyze in detail the predictions of our model for the abundance of elements like Eu, Ba, and Sr in the Galactic halo, following their production by $r$-process nucleosynthesis for different assumptions on the supernova mass range. Finally, we compare our results with spectroscopic data for the chemical composition of metal-poor halo stars.

We present the results of a study aimed at determining the 12C/13C ratio in two samples of planetary nebulae (PNe) by means of millimeter wave observations of 12CO and 13CO. The first group includes six PNe which have been observed in the 3He+ hyperfine transition by Balser et al. (1997); the other group consists of 22 nebulae with rich molecular envelopes. We have determined the carbon isotopic ratio in 14 objects, 9 of which are new detections. The results indicate a range of values of 12C/13C between 9 and 23. We estimate the mass of the progenitors of the PNe of our sample and combine this information with the derived 12C/13C isotopic ratios to test the predictions of stellar nucleosynthesis models. We find that the majority of PNe have isotopic ratios below the values expected from current standard asymptotic giant branch models in the mass range of interest. We suggest that the progenitors of the PNe must have undergone a non-standard mixing process during their red giant phase and/or asymptotic giant phase, resulting in a significant enhancement of the 13C abundance in the surface layers. Our study confirms a similar behaviour inferred from spectroscopic observations of field population II stars and globular cluster giants, and extends it to the final stages of stellar evolution. Finally, we discuss the implications of our results on models of Galactic chemical evolution of 3He and 12C/13C.

We describe the basic features of a Monte Carlo model specifically designed to follow the inhomogenous chemical evolution of the Galactic halo, taking into account the effects of local enrichment and mixing of the halo gas, and with particular emphasis on elements like Eu produced by r-process nucleosynthesis. We compare our results with spectroscopic data for the chemical composition of metal-poor halo stars and globular clusters like M13, M5, M92, M4, and we infer some constraints on the star formation history of the halo and the rate and mass spectrum of supernovae during the first epoch of Galaxy evolution.

We investigate the equilibrium properties of self-gravitating magnetized clouds with polytropic equations of state with negative index n. In particular, we consider scale-free isopedic configurations that have constant dimensionless spherical mass-to-flux ratio lambda_r and that may constitute ``pivotal'' states for subsequent dynamical collapse to form groups or clusters of stars. For given Gamma = 1 + 1/n, equilibria with smaller values of lambda_r are more flattened, ranging from spherical configurations with lambda_r = infinity to completely flattened states for lambda_r=1. For a given amount of support provided by the magnetic field as measured by the dimensionless parameter H_0, equilibria with smaller values of Gamma are more flattened. However, logatropic (defined by Gamma = 0) disks do not exist. The only possible scale-free isopedic equilibria with logatropic equation of state are spherical uniformly magnetized clouds.

We follow the chemical evolution of the Galaxy for elements from Ba to Eu, using an evolutionary model suitable to reproduce a large set of Galactic (local and non local) and extragalactic constraints. Input stellar yields for neutron-rich nuclei have been separated into their s-process and r-process components. The production of s-process elements in thermally pulsing asymptotic giant branch stars of low mass proceeds from the combined operation of two neutron sources: the dominant reaction 13C(alpha,n)16O, which releases neutrons in radiative conditions during the interpulse phase, and the reaction 22Ne(alpha,n)25Mg, marginally activated during thermal instabilities. The resulting s-process distribution is strongly dependent on the stellar metallicity. For the standard model discussed in this paper, it shows a sharp production of the Ba-peak elements around Z = Z_sun/4. Concerning the r-process yields, we assume that the production of r-nuclei is a primary process occurring in stars near the lowest mass limit for Type II supernova progenitors. The r-contribution to each nucleus is computed as the difference between its solar abundance and its s-contribution given by the Galactic chemical evolution model at the epoch of the solar system formation. We compare our results with spectroscopic abundances of elements from Ba to Eu at various metallicities (mainly from F and G stars) showing that the observed trends can be understood in the light of the present knowledge of neutron capture nucleosynthesis. Finally, we discuss a number of emerging features that deserve further scrutiny.

The chemical evolution of the Galaxy is followed for the elements affected by neutron capture, in particular for those in the atomic number range 56 to 63 (Ba, La, Ce, Pr, Nd, Sm and Eu). Neutrons by the major 13^C source, released in radiative conditions in the interpulse periods of TP-AGB stars, give rise to an efficient s-processing, making low-mass AGB the major contributors to the chemical evolution of heavy elements. The s-process scenario, characterized is analyzed using AGB stellar evolutionary calculations with the FRANEC Code (FRASCATI Raphson-Newton Evolutionary Code), and are applied over a wide range of stellar masses and metallicities. The presence of r-process elements in low metallicity stars is indicative of a prompt enrichment of the Galaxy by early generation of stars, and low mass SNII appear to be good candidates for primary production of r-nuclei. The chemical evolution model used here is organized over three-zone, halo, thick and thin disk.

The process of molecule formation in the primordial gas is considered in the framework of Friedmann cosmological models from redshift z=1e4 to z=0. First, a comprehensive analysis of 87 gas phase reaction rates (both collisional and radiative) relevant in the physical environment of the expanding universe is presented and critically discussed. On this basis, calculations are carried out of the abundance of 21 molecular species as function of redshift evaluating consistently the molecular heating and cooling due to H2, HD and LiH molecules. One of the major improvements of this work is the use of a better treatment of H recombination that leads to a reduction of a factor 2-3 in the abundance of electrons and H+ at freeze-out, with respect to previous studies. We find that in the standard model, the residual fractional ionization at z=1 is 3.02e-4, and the main molecular species fractional abundances H2/H=1.1e-6, HD/H2=1.1e-3, HeH+/H=6.2e-13, LiH+/H=9.4e-18, and LiH/LiH+=7.6e-3. Finally, we provide accurate cooling functions of H2, HD and LiH in a wide range of density and temperature that can be conveniently used in a variety of cosmological applications.

We present the results of a study aimed at determining the 12C/13C ratio in two samples of planetary nebuale (PNe) by means of mm-wave observations of 12CO and 13CO. The first group includes six PNe which have been observed in the 3He+ hyperfine transition; the other group consists of 23 nebulae with rich molecular envelopes. We have determined the isotopic ratio in 14 objects and the results indicate a range of values between 9 and 23. In particular, three PNe have ratios well below the value predicted by standard evolutionary models (about 20), indicating that some extra-mixing process has occurred in these stars. We briefly discuss the implications of our results for standard and nonstandard stellar nucleosynthesis.

We show that gaseous \HI disks of primordial composition irradiated by an external radiation field can develop a multiphase medium with temperatures between 10^2 and 10^4 K due to the formation of molecular hydrogen. For a given \HI column density there is a critical value of the radiation field below which only the cold \HI phase can exist. Due to a time decreasing quasar background, the gas starts cooling slowly after recombination until the lowest stable temperature in the warm phase is reached at a critical redshift $z=z_{cr}$. Below this redshift the formation of molecular hydrogen promotes a rapid transition towards the cold \HI phase. We find that disks of protogalaxies with $10^{20}\simlt N_{HI}\simlt 10^{21} cm^{-2}$ are gravitationally stable at $T\sim 10^4$ K and can start their star formation history only at $z \simlt z_{cr}\sim 2$, after the gas in the central portion of the disk has cooled to temperatures $T\simlt 300$ K. Such a delayed starbust phase in galaxies of low gas surface density and low dynamical mass can disrupt the disks and cause them to fade away. These objects could contribute significantly to the faint blue galaxy population.

We examine the processes of formation and destruction of LiH molecules in the primordial gas at temperatures T \leq 5000 K in the framework of standard Friedmann cosmological models and we compute the optical depth of the Universe due to Thomson scattering of cosmic background photons on LiH. With the help of a simple model of evolution of a spherical density perturbation we follow the linear growth and the subsequent collapse of protogalactic clouds of various masses, evaluating consistently their physical characteristics and chemical composition. We determine the expected level of anisotropy of the cosmic background radiation due to the presence of LiH molecules in high-redshift protogalactic clouds. We conclude that LiH spectral features generated in the primordial gas are hardly detectable with current techniques and instrumentation.

We consider the problem of the apparent lack of old T Tauri stars in low-mass star forming regions in the framework of the standard model of low-mass star formation. We argue that the similarity between molecular cloud lifetime and ambipolar diffusion timescale implies that star formation does not take place instantaneously, nor at a constant rate. We conclude that the probability of finding a large population of old stars in a star forming region is intrinsically very small and that the post-T Tauri problem is by and large not existent.

The discrepancy between the observed abundances of 3He in the ISM and those predicted by stellar and galactic chemical evolution remains largely unexplained. In this paper, we attempt to shed some light on this unsolved problem by presenting a quantitative comparison of the 3He abundances recently measured in six planetary nebulae (PNe) with the corresponding predictions of stellar evolution theory. The determination of the mass of the PNe progenitors allows us to dismiss, to a good degree of confidence, the hypothesis that the abundance of 3He in the envelope of all low-mass stars is strongly reduced with respect to the standard theoretical values by some mixing mechanism acting in the latest phases of stellar evolution. The abundance versus mass correlation, allowance made for the limitation of the sample, is in fact found to be fully consistent with the classical prediction of stellar evolution. We examine the implications of this result on the galactic evolution of 3He with the help of a series of models with standard and non-standard nucleosynthesis prescriptions. The results are found to be consistent with the observed galactic abundances only if the vast majority of low-mass stars follows non-standard prescriptions. This implies that either the sample of PNe nebulae under exam is highly biased, or the solution to the 3He problem lies elsewhere.