The degree of coupling between the gas and the magnetic field during the collapse of a core and the subsequent formation of a disk depends on the assumed dust size distribution. We study the impact of grain-grain coagulation on the evolution of magnetohydrodynamic (MHD) resistivities during the collapse of a prestellar core. We use a 1-D model to follow the evolution of the dust size distribution, out-of-equilibrium ionization state and gas chemistry during the collapse of a prestellar core. To compute the grain-grain collisional rate, we consider models for both random and systematic, size-dependent, velocities. We include grain growth through grain-grain coagulation and ice accretion, but ignore grain fragmentation. Starting with a MRN (Mathis et al. 1977) size distribution, we find that coagulation in grain-grain collisions generated by hydrodynamical turbulence is not efficient at removing the smallest grains, and as a consequence does not affect much the evolution of the Hall and ambipolar diffusion MHD resistivities which still severly drop during the collapse like in models without coagulation. The inclusion of systematic velocities, possibly induced by the presence of ambipolar diffusion, increases the coagulation rate between small and large grains, removing small grains earlier in the collapse and therefore limiting the drop in the Hall and ambipolar diffusion resistivities. At intermediate densities ($n_{\rm H} \sim 10^8\,$cm$^{-3}$), the Hall and ambipolar diffusion resistivities are found to be higher by 1 to 2 orders of magnitude in models with coagulation than in models where coagulation is ignored, and also higher than in a toy model without coagulation where all grains smaller than $0.1\,{\mu}$m would have been removed in the parent cloud before the collapse. When grain drift velocities induced by ambipolar diffusion are included, dust coagulation ... (abridged)

The question of what regulates star formation is a long standing issue. To investigate this issue, we run simulations of a kiloparsec cube section of a galaxy with three kinds of stellar feedback: the formation of HII regions, the explosion of supernovae, and the UV heating. We show that stellar feedback is sufficient to reduce the averaged star formation rate (SFR) to the level of the Schmidt- Kennicutt law in Milky-Way like galaxies but not in high-redshift gas rich galaxies suggesting that another type of support should be added. We investigate whether an external driving of the turbulence such as the one created by the large galactic scales could diminish the SFR at the observed level. Assuming that the Toomre parameter is close to 1 as suggested by the observations, we infer a typical turbulent forcing that we argue should be applied parallel to the plane of the galactic disc. When this forcing is applied in our simulations, the SFR within our simulations closely follows the Schmidt- Kennicutt relation. We found that the velocity dispersion is strongly anisotropic with the velocity dispersion alongside the galactic plane being up to 10 times larger than the perpendicular velocity.

In this chapter, we review some historical understanding and recent advances on the Initial Mass Function (IMF) and the Core Mass Function (CMF), both in terms of observations and theories. We focus mostly on star formation in clustered environment since this is suggested by observations to be the dominant mode of star formation. The statistical properties and the fragmentation behaviour of turbulent gas is discussed, and we also discuss the formation of binaries and small multiple systems.

Molecular clouds are a fundamental ingredient of galaxies: they are the channels that transform the diffuse gas into stars. The detailed process of how they do it is not completely understood. We review the current knowledge of molecular clouds and their substructure from scales $\sim~$1~kpc down to the filament and core scale. We first review the mechanisms of cloud formation from the warm diffuse interstellar medium down to the cold and dense molecular clouds, the process of molecule formation and the role of the thermal and gravitational instabilities. We also discuss the main physical mechanisms through which clouds gather their mass, and note that all of them may have a role at various stages of the process. In order to understand the dynamics of clouds we then give a critical review of the widely used virial theorem, and its relation to the measurable properties of molecular clouds. Since these properties are the tools we have for understanding the dynamical state of clouds, we critically analyse them. We finally discuss the ubiquitous filamentary structure of molecular clouds and its connection to prestellar cores and star formation.

[abridged] Understanding how the infalling gas redistribute most of its initial angular momentum inherited from prestellar cores before reaching the stellar embryo is a key question. Disk formation has been naturally considered as a possible solution to this "angular momentum problem". However, how the initial angular momentum of protostellar cores is distributed and evolves during the main accretion phase and the beginning of disk formation has largely remained unconstrained up to now. In the framework of the IRAM CALYPSO survey, we used high dynamic range C$^{18}$O (2-1) and N$_2$H$^+$ (1-0) observations to quantify the distribution of specific angular momentum along the equatorial axis in a sample of 12 Class 0 protostellar envelopes from scales ~50 to 10000 au. The radial distributions of specific angular momentum in the CALYPSO sample suggest two distinct regimes within protostellar envelopes: the specific angular momentum decreases as $j \propto r^{1.6 \pm 0.2}$ down to ~1600 au and then tends to become relatively constant around 6 $\times$ 10$^{-4}$ km s$^{-1}$ pc down to ~50 au. The values of specific angular momentum measured in the inner Class 0 envelopes, namely that of the material directly involved in the star formation process ($<$1600 au), is on the same order of magnitude as what is inferred in small T-Tauri disks. Thus, disk formation appears to be a direct consequence of angular momentum conservation during the collapse. Our analysis reveals a dispersion of the directions of velocity gradients at envelope scales $>$1600 au, suggesting that they may not be related to rotational motions of the envelopes. We conclude that the specific angular momentum observed at these scales could find its origin in core-forming motions (infall, turbulence) or trace an imprint of the initial conditions for the formation of protostellar cores.

Star formation is a complex multi-scale phenomenon that is of significant importance for astrophysics in general. Stars and star formation are key pillars in observational astronomy from local star forming regions in the Milky Way up to high-redshift galaxies. From a theoretical perspective, star formation and feedback processes (radiation, winds, and supernovae) play a pivotal role in advancing our understanding of the physical processes at work, both individually and of their interactions. In this review we will give an overview of the main processes that are important for the understanding of star formation. We start with an observationally motivated view on star formation from a global perspective and outline the general paradigm of the life-cycle of molecular clouds, in which star formation is the key process to close the cycle. After that we focus on the thermal and chemical aspects in star forming regions, discuss turbulence and magnetic fields as well as gravitational forces. Finally, we review the most important stellar feedback mechanisms.

Giant molecular clouds (GMCs) and their stellar offspring are the building blocks of galaxies. The physical characteristics of GMCs and their evolution are tightly connected to galaxy evolution. The macroscopic properties of the interstellar medium propagate into the properties of GMCs condensing out of it, with correlations between e.g. the galactic and GMC scale gas pressures, surface densities and volume densities. That way, the galactic environment sets the initial conditions for star formation within GMCs. After the onset of massive star formation, stellar feedback from e.g. photoionisation, stellar winds, and supernovae eventually contributes to dispersing the parent cloud, depositing energy, momentum and metals into the surrounding medium, thereby changing the properties of galaxies. This cycling of matter between gas and stars, governed by star formation and feedback, is therefore a major driver of galaxy evolution. Much of the recent debate has focused on the durations of the various evolutionary phases that constitute this cycle in galaxies, and what these can teach us about the physical mechanisms driving the cycle. We review results from observational, theoretical, and numerical work to build a dynamical picture of the evolutionary lifecycle of GMC evolution, star formation, and feedback in galaxies.

Recent observations have suggested that circumstellar disks may commonly form around young stellar objects. Although the formation of circumstellar disks can be a natural result of the conservation of angular momentum in the parent cloud, theoretical studies instead show disk formation to be difficult from dense molecular cores magnetized to a realistic level, owing to efficient magnetic braking that transports a large fraction of the angular momentum away from the circumstellar region. We review recent progress in the formation and early evolution of disks around young stellar objects of both low-mass and high-mass, with an emphasis on mechanisms that may bridge the gap between observation and theory, including non-ideal MHD effects and asymmetric perturbations in the collapsing core (e.g., magnetic field misalignment and turbulence). We also address the associated processes of outflow launching and the formation of multiple systems, and discuss possible implications in properties of protoplanetary disks.

Planet-forming disks are fundamental objects thought to be inherited from large scale rotation, through the conservation of angular momentum during the collapse of a prestellar dense core. We investigate the possibility for a protostellar disk to be formed from a motionless dense core which contains non-axisymmetric density fluctuations. The rotation is thus generated locally by the asymmetry of the collapse. We study the evolution of the angular momentum in a non-axisymmetric collapse of a dense core from an analytical point of view. To test the theory, we perform three-dimensional simulations of a collapsing prestellar dense core using adaptative mesh refinement. We start from a non-axisymmetrical situation, considering a dense core with random density perturbations that follow a turbulence spectrum. We analyse the emerging disk comparing the angular momentum it contains with the one expected from our analytic development. We study the velocity gradients at different scales in the simulation as it is done with observations. We show that the angular momentum in the frame of a stellar object which is not located at the center of mass of the core is not conserved, due to inertial forces. Our simulations of such non-axisymmetrical collapse quickly produce accretion disks at the small scales in the core. The analysis of the kinematics at different scales in the simulated core reveals projected velocity gradients of amplitudes similar to the ones observed in protostellar cores, and which directions vary, sometimes even reversing when small and large scales are compared. These complex kinematics patterns appear in recent observations, and could be a discriminating feature with models where rotation is inherited from large scales. Our results from simulations without initial rotation are more consistent with these recent observations than when solid-body rotation is initially [abridged]

The formation of protoplanetary disks is not well understood. To understand how and when these disks are formed, it is crucial to characterize the kinematics of the youngest protostars at a high angular resolution. Here we study a sample of 16 Class 0 protostars to measure their rotation profile at scales from 50 to 500~au and search for Keplerian rotation. We used high-angular-resolution line observations obtained with the Plateau de Bure Interferometer as part of the CALYPSO large program. From $\mathrm{^{13}CO}$ ($J=2-1$), $\mathrm{C^{18}O}$ ($J=2-1$) and SO ($N_{j}=5_{6}-4_{5}$) moment maps, we find that seven sources show rotation about the jet axis at a few hundred au scales: SerpS-MM18, L1448-C, L1448-NB, L1527, NGC1333- IRAS2A, NGC1333-IRAS4B, and SVS13-B. We analyzed the kinematics of these sources in the $uv$ plane to derive the rotation profiles down to 50~au scales. We find evidence for Keplerian rotation in only two sources, L1527 and L1448-C. Overall, this suggests that Keplerian disks larger than 50 au are uncommon around Class 0 protostars. However, in some of the sources, the line emission could be optically thick and dominated by the envelope emission. Due to the optical thickness of these envelopes, some of the disks could have remained undetected in our observations.

Planets form in protoplanetary discs. Their masses, distribution, and orbits sensitively depend on the structure of the protoplanetary discs. However, what sets the initial structure of the discs in terms of mass, radius and accretion rate is still unknown. We perform non-ideal MHD numerical simulations using the adaptive mesh refinement code Ramses, of a collapsing, one solar mass, molecular core to study the disc formation and early, up to 100 kyr, evolution, paying great attention to the impact of numerical resolution and accretion scheme. We found that while the mass of the central object is almost independent of the numerical parameters such as the resolution and the accretion scheme onto the sink particle, the disc mass, and to a lower extent its size, heavily depend on the accretion scheme, which we found, is itself resolution dependent. This implies that the accretion onto the star and through the disc are largely decoupled. For a relatively large domain of initial conditions (except at low magnetisation), we found that the properties of the disc do not change too significantly. In particular both the level of initial rotation and turbulence do not influence the disc properties provide the core is sufficiently magnetized. After a short relaxation phase, the disc settles in a stationary state. It then slowly grows in size but not in mass. The disc itself is weakly magnetized but its immediate surrounding is on the contrary highly magnetized. Our results show that the disc properties directly depend on the inner boundary condition, i.e. the accretion scheme onto the central object, suggesting that the disc mass is eventually controlled by the small scale accretion process, possibly the star-disc interaction. Because of ambipolar diffusion and its significant resistivity, the disc diversity remains limited and except for low magnetisation, their properties are (abridged).

On the basis of its low luminosity, its chemical composition, and the absence of a large-scale outflow, the dense core Cha-MMS1 located in the Chamaeleon I molecular cloud was proposed as a first hydrostatic core (FHSC) candidate a decade ago. Our goal is to test this hypothesis by searching for a slow, compact outflow driven by Cha-MMS1 that would match the predictions of MHD simulations for this short phase of star formation. We use the Atacama Large Millimetre/submillimetre Array (ALMA) to map Cha-MMS1 at high angular resolution in CO 3-2 and 13CO 3-2 as well as in continuum emission. We report the detection of a bipolar outflow emanating from the central core, along a (projected) direction roughly parallel to the filament in which Cha-MMS1 is embedded and perpendicular to the large-scale magnetic field. The morphology of the outflow indicates that its axis lies close to the plane of the sky. We measure velocities corrected for inclination of more than 90km/s which is clearly incompatible with the expected properties of a FHSC outflow. Several properties of the outflow are determined and compared to previous studies of Class 0 and Class I protostars. The outflow of Cha-MMS1 has a much smaller momentum force than the outflows of other Class 0 protostars. In addition, we find a dynamical age of 200-3000yr indicating that Cha-MMS1 might be one of the youngest ever observed Class 0 protostars. While the existence of the outflow suggests the presence of a disk, no disk is detected in continuum emission and we derive an upper limit of 55au to its radius. We conclude that Cha-MMS1 has already gone through the FHSC phase and is a young Class 0 protostar, but it has not brought its outflow to full power yet.

The interstellar medium (ISM) is a complex non-linear system governed by gravity and magneto-hydrodynamics, as well as radiative, thermodynamical, and chemical processes. Our understanding of it mostly progresses through observations and numerical simulations, and a quantitative comparison between these two approaches requires a generic and comprehensive statistical description. The goal of this paper is to build such a description, with the purpose to permit an efficient comparison independent of any specific prior or model. We start from the Wavelet Scattering Transform (WST), a low-variance statistical description of non-Gaussian processes, developed in data science, that encodes long-range interactions through a hierarchical multiscale approach based on the Wavelet transform. We perform a reduction of the WST through a fit of its angular dependencies, allowing to gather most of the information it contains into a few components whose physical meanings are identified, and that describe, e.g., isotropic and anisotropic behaviours. The result of this paper is the Reduced Wavelet Scattering Transform (RWST), a statistical description with a small number of coefficients that characterizes complex structures arising from non-linear phenomena, free from any specific prior. The RWST coefficients encode moments of order up to four, have reduced variances, and quantify the couplings between scales. To show the efficiency and generality of this description, we apply it successfully to three kinds of processes: fractional Brownian motions, MHD simulations, and Herschel observations in a molecular cloud. With fewer than 100 coefficients when probing 6 scales and 8 angles on 256*256 maps, we were able with the RWST to perform quantitative comparisons, to infer relevant physical properties, and to produce realistic synthetic fields.

We present 0.25 arcsec resolution (35 au) ALMA 1.3 mm dust polarization observations for 37 young stellar objects (YSOs) in the Ophiuchus cloud. These data encompass all the embedded protostars in the cloud and several Flat and Class II objects to produce the largest, homogeneous study of dust polarization on disk scales to date. The goal of this study is to study dust polarization down to disk scales. We find that 14/37 (38%) of the YSOs are detected in polarization. Nine of these sources have uniform polarization angles and four sources have azimuthal polarization. The sources with uniform polarization tend to have steeper inclinations (> 60 degree) than those with azimuthal polarization (< 60 degree). The majority (9/14) of the detected sources have polarization morphologies and disk properties consistent with dust self-scattering in optically thick disks. The remaining sources may be instead tracing magnetic fields. Their inferred field directions from rotating the polarization vectors by 90 degree are mainly poloidal or hourglass shaped. We find no evidence of a strong toroidal field component toward any of our disks. For the 23 YSOs that are undetected in polarization, roughly half of them have 3-sigma upper limits of < 2%. These sources also tend to have inclinations < 60 degree and they are generally compact. Since lower inclination sources tend to have azimuthal polarization, these YSOs may be undetected in polarization due to unresolved polarization structure within our beam. We propose that disks with inclinations > 60 degree are the best candidates for future polarization studies of dust self-scattering as these systems will generally show uniform polarization vectors that do not require very high resolution to resolve. We release the continuum and polarization images for all the sources with this publication. Data from the entire survey can be obtained from Dataverse.

The stellar initial mass function (IMF) is playing a critical role in the history of our universe. We propose a theory that is based solely on local processes, namely the dust opacity limit, the tidal forces and the properties of the collapsing gas envelope. The idea is that the final mass of the central object is determined by the location of the nearest fragments, which accrete the gas located further away, preventing it to fall onto the central object. To estimate the relevant statistics in the neighbourhood of an accreting protostar, we perform high resolution numerical simulations. We also use these simulations to further test the idea that fragmentation in the vicinity of an existing protostar is determinant in setting the peak of the stellar mass spectrum. We develop an analytical model, which is based on a statistical counting of the turbulent density fluctuations, generated during the collapse, that are at least equal to the mass of the first hydrostatic core, and sufficiently important to supersede tidal and pressure forces to be self-gravitating. The analytical mass function presents a peak located at roughly 10 times the mass of the first hydrostatic core in good agreement with the numerical simulations. Since the physical processes involved are all local, i.e. occurs at scales of a few 100 AU or below, and do not depend on the gas distribution at large scale and global properties such as the mean Jeans mass, the mass spectrum is expected to be relatively universal.

How and when in the star formation sequence do dust grains start to grow into pebbles is a cornerstone question to both star and planet formation. We compute the polarized radiative transfer from a model solar-type protostellar core, using the POLARIS code, aligning the dust grains with the local magnetic field, following the radiative torques (RATs) theory. We test the dependency of the resulting dust polarized emission with the maximum grain size of the dust size distribution at the envelope scale, from amax = 1 micron to 50 micron. Our work shows that, in the framework of RAT alignment, large dust grains are required to produce polarized dust emission at levels similar to those currently observed in solar-type protostellar envelopes at millimeter wavelengths. Considering the current theoretical dificulties to align a large fraction of small ISM-like grains in the conditions typical of protostellar envelopes, our results suggest that grain growth (typically > 10 micron) might have already significantly progressed at scales 100-1000 au in the youngest objects, observed less than 10^5 years after the onset of collapse. Observations of dust polarized emission might open a new avenue to explore dust pristine properties and describe, for example, the initial conditions for the formation of planetesimals.

High-mass stars and star clusters commonly form within hub-filament systems. Monoceros R2, harbors one of the closest such systems, making it an excellent target for case studies. We investigate the morphology, stability and dynamical properties of the hub-filament system on basis of 13CO and C18O observations obtained with the IRAM-30m telescope and H2 column density maps derived from Herschel dust emission observations. We identified the filamentary network and characterized the individual filaments as either main (converging into the hub) or secondary (converging to a main filament) filaments. The main filaments have line masses of 30-100 Msun/pc and show signs of fragmentation. The secondary filaments have line masses of 12-60 Msun/pc and show fragmentation only sporadically. In the context of Ostriker's hydrostatic filament model, the main filaments are thermally super-critical. If non-thermal motions are included, most of them are trans-critical. Most of the secondary filaments are roughly trans-critical regardless of whether non-thermal motions are included or not. From the main filaments, we estimate a mass accretion rate of 10(-4)-10(-3) Msun/pc into the hub. The secondary filaments accrete into the main filaments with a rate of 0.1-0.4x10(-4) Msun/pc. The main filaments extend into the hub. Their velocity gradients increase towards the hub, suggesting acceleration of the gas. We estimate that with the observed infall velocity, the mass-doubling time of the hub is ~2.5 Myr, ten times larger than the free-fall time, suggesting a dynamically old region. These timescales are comparable with the chemical age of the HII region. Inside the hub, the main filaments show a ring- or a spiral-like morphology that exhibits rotation and infall motions. One possible explanation for the morphology is that gas is falling into the central cluster following a spiral-like pattern.

Molecular clouds are turbulent structures whose star formation efficiency (SFE) is strongly affected by internal stellar feedback processes. In this paper we determine how sensitive the SFE of molecular clouds is to randomised inputs in the star formation feedback loop, and to what extent relationships between emergent cloud properties and the SFE can be recovered. We introduce the yule suite of 26 radiative magnetohydrodynamic (RMHD) simulations of a 10,000 solar mass cloud similar to those in the solar neighbourhood. We use the same initial global properties in every simulation but vary the initial mass function (IMF) sampling and initial cloud velocity structure. The final SFE lies between 6 and 23 percent when either of these parameters are changed. We use Bayesian mixed-effects models to uncover trends in the SFE. The number of photons emitted early in the cluster's life and the length of the cloud provide are the strongest predictors of the SFE. The HII regions evolve following an analytic model of expansion into a roughly isothermal density field. The more efficient feedback is at evaporating the cloud, the less the star cluster is dispersed. We argue that this is because if the gas is evaporated slowly, the stars are dragged outwards towards surviving gas clumps due to the gravitational attraction between the stars and gas. While star formation and feedback efficiencies are dependent on nonlinear processes, statistical models describing cloud-scale processes can be constructed.

SPICA, the cryogenic infrared space telescope recently pre-selected for a `Phase A' concept study as one of the three remaining candidates for ESA's fifth medium class (M5) mission, is foreseen to include a far-infrared polarimetric imager (SPICA-POL, now called B-BOP), which would offer a unique opportunity to resolve major issues in our understanding of the nearby, cold magnetized Universe. This paper presents an overview of the main science drivers for B-BOP, including high dynamic range polarimetric imaging of the cold interstellar medium (ISM) in both our Milky Way and nearby galaxies. Thanks to a cooled telescope, B-BOP will deliver wide-field 100-350 micron images of linearly polarized dust emission in Stokes Q and U with a resolution, signal-to-noise ratio, and both intensity and spatial dynamic ranges comparable to those achieved by Herschel images of the cold ISM in total intensity (Stokes I). The B-BOP 200 micron images will also have a factor ~30 higher resolution than Planck polarization data. This will make B-BOP a unique tool for characterizing the statistical properties of the magnetized interstellar medium and probing the role of magnetic fields in the formation and evolution of the interstellar web of dusty molecular filaments giving birth to most stars in our Galaxy. B-BOP will also be a powerful instrument for studying the magnetism of nearby galaxies and testing galactic dynamo models, constraining the physics of dust grain alignment, informing the problem of the interaction of cosmic rays with molecular clouds, tracing magnetic fields in the inner layers of protoplanetary disks, and monitoring accretion bursts in embedded protostars.

The connection between the pre-stellar core mass function (CMF) and the stellar initial mass function (IMF) lies at the heart of all star formation theories. In this paper, we study the earliest phases of star formation with a series of high-resolution numerical simulations that include the formation of sinks. In particular, we focus on the transition from cores to sinks within a massive molecular filament. We compare the CMF and IMF between magnetized and unmagnetized simulations, and between different resolutions. We find that selecting cores based on their kinematic virial parameter excludes collapsing objects because they host large velocity dispersions. Selecting only the thermally unstable magnetized cores, we observe that their mass-to-flux ratio spans almost two orders of magnitude for a given mass. We also see that, when magnetic fields are included, the CMF peaks at higher core mass values with respect to pure hydrodynamical simulations. Nonetheless, all models produce sink mass functions with a high-mass slope consistent with Salpeter. Finally, we examine the effects of resolution and find that, in isothermal simulations, even models with very high dynamical range fail to converge in the mass function. Our main conclusion is that, although the resulting CMFs and IMFs have similar slopes in all simulations, the cores have slightly different sizes and kinematical properties when a magnetic field is included. However, a core selection based on the mass-to-flux ratio alone is not enough to alter the shape of the CMF, if we do not take thermal stability into account. Finally, we conclude that extreme care should be given to resolution issues when studying sink formation with an isothermal equation of state.

We review the role that magnetic field may have on the formation and evolution of molecular clouds. After a brief presentation and main assumptions leading to ideal MHD equations, their most important correction, namely the ion-neutral drift is described. The nature of the multi-phase interstellar medium (ISM) and the thermal processes that allows this gas to become denser are presented. Then we discuss our current knowledge of compressible magnetized turbulence, thought to play a fundamental role in the ISM. We also describe what is known regarding the correlation between the magnetic and the density fields. Then the influence that magnetic field may have on the interstellar filaments and the molecular clouds is discussed, notably the role it may have on the prestellar dense cores as well as regarding the formation of stellar clusters. Finally we briefly review its possible effects on the formation of molecular clouds themselves. We argue that given the magnetic intensities that have been measured, it is likely that magnetic field is i) responsible of reducing the star formation rate in dense molecular cloud gas by a factor of a few, ii) strongly shaping the interstellar gas by generating a lot of filaments and reducing the numbers of clumps, cores and stars, although its exact influence remains to be better understood. % by a factor on the order of at least 2. Moreover at small scales, magnetic braking is likely a dominant process that strongly modifies the outcome of the star formation process. Finally, we stress that by inducing the formation of more massive stars, magnetic field could possibly enhance the impact of stellar feedback.

We introduce the histogram of oriented gradients (HOG), a tool developed for machine vision that we propose as a new metric for the systematic characterization of observations of atomic and molecular gas and the study of molecular cloud formation models. In essence, the HOG technique takes as input extended spectral-line observations from two tracers and provides an estimate of their spatial correlation across velocity channels. We characterize HOG using synthetic observations of HI and $^{13}$CO(J=1-0) emission from numerical simulations of MHD turbulence leading to the formation of molecular gas after the collision of two atomic clouds. We find a significant spatial correlation between the two tracers in velocity channels where $v_{HI}\approx v_{^{13}CO}$, independent of the orientation of the collision with respect to the line of sight. We use HOG to investigate the spatial correlation of the HI, from the THOR survey, and the $^{13}$CO(J=1-0) emission, from the GRS, toward the portion of the Galactic plane 33.75$\lt l\lt$35.25$^{o}$ and $|b|\lt$1.25$^{o}$. We find a significant spatial correlation between the tracers in extended portions of the studied region. Although some of the regions with high spatial correlation are associated with HI self-absorption features, suggesting that it is produced by the cold atomic gas, the correlation is not exclusive to this kind of region. The HOG results also indicate significant differences between individual regions: some show spatial correlation in channels around $v_{HI}\approx v_{^{13}CO}$ while others present this correlation in velocity channels separated by a few km/s. We associate these velocity offsets to the effect of feedback and to the presence of physical conditions that are not included in the atomic-cloud-collision simulations, such as more general magnetic field configurations, shear, and global gas infall.

The stellar mass spectrum is an important property of the stellar cluster and a fundamental quantity to understand our Universe. The fragmentation of diffuse molecular cloud into stars is subject to physical processes such as gravity, turbulence, thermal pressure, and magnetic field. The final mass of a star is believed to be a combined outcome of a virially unstable reservoir and subsequent accretion. We aim to clarify the roles of different supporting energies, notably the thermal pressure and the magnetic field, in determining the stellar mass. Following previous studies by Lee & Hennebelle (2018a,b), we perform a series of numerical experiments of stellar cluster formation inside an isolated molecular clump. By changing the effective equation of state (EOS) of the diffuse gas (that is to say gas whose density is below the critical density at which dust becomes opaque to its radiation) and the strength of the magnetic field, we investigate whether any characteristic mass is introduced into the fragmentation processes. The EOS of the diffuse gas, including the bulk temperature and the polytropic index, does not affect significantly the shape of the stellar mass spectrum. The presence of magnetic field slightly modifies the shape of the mass spectrum only when extreme values are applied. This study confirms that the peak of the IMF is primarily determined by the adiabatic high-density end of the EOS that mimics the radiation inside the high-density gas. Furthermore, the shape of the mass spectrum is mostly sensitive to the density PDF, and the magnetic field has likely only a secondary role. In particular, we stress that the Jeans mass at the mean cloud density and at the critical density are not responsible of setting the peak.

Understanding the formation mechanisms of protoplanetary disks and multiple systems, and their pristine properties, is a key question for modern astrophysics. The properties of the youngest disks, embedded in rotating infalling protostellar envelopes, have largely remained unconstrained up to now. In the framework of the IRAM-PdBI CALYPSO survey, we have obtained sub-arcsecond observations of the dust continuum emission at 231 GHz and 94 GHz, for a sample of 16 solar-type Class 0 protostars. In an attempt to identify disk-like structures embedded at small scales in the protostellar envelopes, we model the dust continuum emission visibility profiles using both Plummer-like envelope models and envelope models including additional Gaussian disk-like components. Our analysis shows that in the CALYPSO sample, 11 of the 16 Class 0 protostars are better reproduced by models including a disk-like dust continuum component contributing to the flux at small scales, but less than 25% of these candidate protostellar disks are resolved at radii > 60 au. Including all available literature constraints on Class 0 disks at subarcsecond scales, we show that our results are representative: most (> 72% in a sample of 26 protostars) Class 0 protostellar disks are small and emerge only at radii < 60 au. Our multiplicity fraction at scales 100-5000 au is in global agreement with the multiplicity properties of Class I protostars at similar scales. We confront our observational constraints on the disk size distribution in Class 0 protostars to the typical disk properties from protostellar formation models. Because they reduce the centrifugal radius, and produce a disk size distribution peaking at radii <100 au during the main accretion phase, magnetized models of rotating protostellar collapse are favored by our observations.

The dynamics of massive clumps, the environment where massive stars originate, is still unclear. Many theories predict that these regions are in a state of near-virial equilibrium, or near energy equi-partition, while others predict that clumps are in a sub-virial state. Observationally, the majority of the massive clumps are in a sub-virial state with a clear anti-correlation between the virial parameter $\alpha_{vir}$ and the mass of the clumps $M_{c}$, which suggests that the more massive objects are also the more gravitationally bound. Although this trend is observed at all scales, from massive clouds down to star-forming cores, theories do not predict it. In this work we show how, starting from virialized clumps, an observational bias is introduced in the specific case where the kinetic and the gravitational energies are estimated in different volumes within clumps and how it can contribute to the spurious $\alpha_{vir}-M_{c}$ anti-correlation in these data. As a result, the observed effective virial parameter $\tilde{\alpha}_{eff}<\alpha_{vir}$, and in some circumstances it might not be representative of the virial state of the observed clumps.

A numerical shearing box is used to perform three-dimensional simulations of a 1 kpc stratified cubic box of turbulent and self-gravitating interstellar medium (in a rotating frame) with supernovae and HII feedback. We vary the value of the velocity gradient induced by the shear and the initial value of the galactic magnetic field. Finally the different star formation rates and the properties of the structures associated with this set of simulations are computed. We first confirm that the feedback has a strong limiting effect on star formation. The galactic shear has also a great influence: the higher the shear, the lower the SFR. Taking the value of the velocity gradient in the solar neighbourhood, the SFR is too high compared to the observed Kennicutt law, by a factor approximately three to six. This discrepancy can be solved by arguing that the relevant value of the shear is not the one in the solar neighbourhood, and that in reality the star formation efficiency within clusters is not 100%. Taking into account the fact that star-forming clouds generally lie in spiral arms where the shear can be substantially higher (as probed by galaxy-scale simulations), the SFR is now close to the observed one. Different numerical recipes have been tested for the sink particles, giving a numerical incertitude of a factor of about two on the SFR. Finally we have also estimated the velocity dispersions in our dense clouds and found that they lie below the observed Larson law by a factor of about two. Conclusions. In our simulations, magnetic field, shear, HII regions, and supernovae all contribute significantly to reduce the SFR. In this numerical setup with feedback from supernovae and HII regions and a relevant value of galactic shear, the SFRs are compatible with those observed, with a numerical incertitude factor of about two.

The multifractal geometry remains an under-exploited approach to describe and quantify the large-scale structure of interstellar clouds. In this paper, the typical tools of multifractal analysis are applied to Herschel far-infrared (70-500 $\mu$m) dust continuum maps, which represent an ideal case of study. This dust component is a relatively optically thin tracer at these wavelengths and the size in pixel of the maps is well suitable for this statistical analysis. We investigate the so-called multifractal spectrum and generalised fractal dimension of six Hi-GAL regions in the third Galactic quadrant. We measure systematic variations of the spectrum at increasing wavelength, which generally correspond to an increasing complexity of the image, and we observe peculiar behaviours of the investigated fields, strictly related to the presence of high-emission regions, which in turn are connected to star formation activity. The same analysis is applied to synthetic column density maps, generated from numerical turbulent molecular cloud models and from fractal Brownian motion (fBm), allowing for the confrontation of the observations with models with well controlled physical parameters. This comparison shows that, in terms of multifractal descriptors of the structure, fBm images exhibit a very different, quite unrealistic behaviour when compared with Hi-GAL observations, whereas the numerical simulations appear in some cases (depending on the specific model) more similar to the observations. Finally, the link between mono-fractal parameters (commonly found in the literature) and multifractal indicators is investigated: the former appear to be only partially connected with the latter, with the multifractal analysis offering a more extended and complete characterization of the cloud structure.

The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration.

The role of the magnetic field during protostellar collapse is poorly constrained from an observational point of view, although it could be significant if we believe state-of-the-art models of protostellar formation. We present polarimetric observations of the 233 GHz thermal dust continuum emission obtained with ALMA in the B335 Class 0 protostar. Linearly polarized dust emission arising from the circumstellar material in the envelope of B335 is detected at all scales probed by our observations, from radii of 50 to 1000 au. The magnetic field structure producing the dust polarization has a very ordered topology in the inner envelope, with a transition from a large-scale poloidal magnetic field, in the outflow direction, to strongly pinched in the equatorial direction. This is probably due to magnetic field lines being dragged along the dominating infall direction since B335 does not exhibit prominent rotation. Our data and their qualitative comparison to a family of magnetized protostellar collapse models show that, during the magnetized collapse in B335, the magnetic field is maintaining a high level of organization from scales 1000 au to 50 au: this suggests the field is dynamically relevant and capable of influencing the typical outcome of protostellar collapse, such as regulating the disk size in B335.

Molecular clouds are essentially made up of atomic and molecular hydrogen, which in spite of being the simplest molecule in the ISM plays a key role in the chemical evolution of molecular clouds. Since its formation time is very long, the H2 molecules can be transported by the turbulent motions within the cloud toward low density and warm regions, where its enhanced abundance can boost the abundances of molecules with high endothermicities. We present high resolution simulations where we include the evolution of the molecular gas under the effect of the dynamics, and we analyze its impact on the abundance of CH+.

We investigate the dependence of the peak of the IMF on the physics of the so-called first Larson core, which corresponds to the point where the dust becomes opaque to its own radiation. We perform numerical simulations of collapsing clouds of $1000 M_\odot$ for various gas equation of state (eos), paying great attention to the numerical resolution and convergence. The initial conditions of these numerical experiments are varied in the companion paper. We also develop analytical models that we confront to our numerical results. If an isothermal eos is used, we show that the peak of the IMF shifts to lower masses with improved numerical resolution. When an adiabatic eos is employed, numerical convergence is obtained. The peak position varies with the eos and we find that the peak position is about ten times the mass of the first Larson core. By analyzing the stability of non-linear density fluctuations in the vicinity of a point mass and then summing over a reasonable density distribution, we find that tidal forces exert a strong stabilizing effect and likely lead to a preferential mass several times larger than that of the first Larson core. We propose that in a sufficiently massive and cold cloud, the peak of the IMF is determined by the thermodynamics of the high density adiabatic gas as well as the stabilizing influence of tidal forces. The resulting characteristic mass is about ten times the mass of the first Larson core, which altogether leads to a few tenths of solar masses. Since these processes are not related to the large scale physical conditions and to the environment, our results suggest a possible explanation for the apparent universality of the peak of the IMF.

We conduct numerical experiments in which we systematically vary the initial density over four orders of magnitude and the turbulent velocity over a factor ten. In a companion paper, we investigate the dependence of this distribution on the gas thermodynamics. We performed a series of hydrodynamical numerical simulations using adaptive mesh refinement, with special attention to numerical convergence. We also adapted an existing analytical model to the case of collapsing clouds by employing a density probability distribution function (PDF) $\propto \rho^{-1.5}$ instead of a lognormal distribution. Simulations and analytical model both show two support regimes, dominated by either thermal energy or turbulence. For the first regime, we infer that $dN/d \log M \propto M^0$, while for the second, we obtain $dN/d \log M \propto M^{-3/4}$. This is valid up to about ten times the mass of the first Larson core, as explained in the companion paper, leading to a peak of the mass spectrum at $\sim 0.2 M_\odot$. From this point, the mass spectrum decreases with decreasing mass. Although the mass spectra we obtain for the most compact clouds qualitatively resemble the observed initial mass function, the distribution exponent is shallower than the expected Salpeter exponent of -1.35. Nonetheless, we observe a possible transition toward a slightly steeper value that is broadly compatible with the Salpeter exponent for masses above a few solar masses. This change in behavior is associated with the change in density PDF, which switches from a power-law to a lognormal distribution. Our results suggest that while gravitationally induced fragmentation could play an important role for low masses, it is likely the turbulently induced fragmentation that leads to the Salpeter exponent.

Abridged. Understanding the detailed structure of the interstellar gas is essential for our knowledge of the star formation process. The small-scale structure of the interstellar medium (ISM) is a direct consequence of the galactic scales and making the link between the two is essential. We perform adaptive mesh simulations that aim to bridge the gap between the intermediate galactic scales and the self-gravitating prestellar cores. For this purpose we use stratified supernova regulated ISM magneto-hydrodynamical (MHD) simulations at the kpc scale to set up the initial conditions. We then zoom, performing a series of concentric uniform refinement and then refining on the Jeans length for the last levels. This allows us to reach a spatial resolution of a few $10^{-3}$ pc. The cores are identified using a clump finder and various criteria based on virial analysis. Their most relevant properties are computed and, due to the large number of objects formed in the simulations, reliable statistics are obtained. The cores properties show encouraging agreements with observations. The mass spectrum presents a clear powerlaw at high masses with an exponent close to $\simeq -1.3$ and a peak at about 1-2 $M_\odot$. The velocity dispersion and the angular momentum distributions are respectively a few times the local sound speed and a few $10^{-2}$ pc km s$^{-1}$. We also find that the distribution of thermally supercritical cores present a range of magnetic mass-to-flux over critical mass-to-flux ratio which typically ranges between $\simeq$0.3 and 3.

Observations suggest that star formation in filamentary molecular clouds occurs in a two-step process, with the formation of filaments preceding that of prestellar cores and stars. Here, we apply the gravo-turbulent fragmentation theory of Hennebelle & Chabrier 08, 09, 13 to a filamentary environment, taking into account magnetic support. We discuss the induced geometrical effect on the cores, with a transition from 3D geometry at small scales to 1D at large ones. The model predicts the fragmentation behavior of a filament for a given mass per unit length (MpL) and level of magnetization. This CMF for individual filaments is then convolved with the distribution of filaments to obtain the final system CMF. The model yields two major results: (i) the filamentary geometry naturally induces a hierarchical fragmentation process, first into groups of cores, separated by a length equal to a few filament Jeans lengths, i.e. a few times the filament width. These groups then fragment into individual cores. (ii) Non-magnetized filaments with high MpL are found to fragment excessively, at odd with observations. This is resolved by taking into account the magnetic field treated simply as additional pressure support). The present theory suggests two complementary modes of star formation: while small (spherical or filamentary) structures will collapse directly into prestellar cores, according to the standard Hennebelle-Chabrier theory, the large (filamentary) ones, the dominant population according to observations, will follow the afore-described two-step process.

Recent observations suggest that intensive molecular cloud collision can trigger massive star/cluster formation. The most important physical process caused by the collision is a shock compression. In this paper, the influence of a shock wave on the evolution of a molecular cloud is studied numerically by using isothermal magnetohydrodynamics (MHD) simulations with the effect of self-gravity. Adaptive-mesh-refinement and sink particle techniques are used to follow long-time evolution of the shocked cloud. We find that the shock compression of turbulent inhomogeneous molecular cloud creates massive filaments, which lie perpendicularly to the background magnetic field as we have pointed out in a previous paper. The massive filament shows global collapse along the filament, which feeds a sink particle located at the collapse center. We observe high accretion rate dot{M}_acc > 10^{-4} M_sun/yr that is high enough to allow the formation of even O-type stars. The most massive sink particle achieves M>50 M_sun in a few times 10^5 yr after the onset of the filament collapse.

Aims. We wish to clarify the origin of the multiple jet features emanating from the binary protostar SVS 13A (= VLA4A/VLA4B). Methods. We used the Plateau de Bure Interferometer to map at 0.3-0.8" (~70-190 au) dust emission at 1.4 mm, CO(2-1), SiO(5-4), SO(65-54). Revised proper motions for VLA4A/4B and jet wiggling models are computed to clarify their respective contribution. Results. VLA4A shows compact dust emission suggestive of a disk < 50 au, and is the hot corino source, while CO/SiO/SO counterparts to the small-scale H2 jet originate from VLA4B and reveal the jet variable velocity structure. This jet exhibits ~ 3" wiggling consistent with orbital motion around a yet undetected ~ 20-30 au companion to VLA4B, or jet precession. Jet wiggling combined with velocity variability can explain the large apparent angular momentum in CO bullets. We also uncover a synchronicity between CO jet bullets and knots in the HH7-11 chain demonstrating that they trace two distinct jets. Their ~ 300 yr twin outburst period may be triggered by close perihelion approach of VLA4A in an eccentric orbit around VLA4B. A third jet is tentatively seen at PA ~ 0 degrees. Conclusions. SVS13 A harbors at least 2 and possibly 3 distinct jet sources. The CO and HH7-11 jets are launched from quasi-coplanar disks, separated by 20-70 au. Their synchronous major events every 300 yr favor external triggering by close binary interactions, a scenario also invoked for FU Or outbursts.

We investigate the origin of observed local star formation relations using radiative magnetohydrodynamic simulations with self-consistent star formation and ionising radiation. We compare these clouds to the density distributions of local star-forming clouds and find that the most diffuse simulated clouds match the observed clouds relatively well. We then compute both observationally-motivated and theoretically-motivated star formation efficiencies (SFEs) for these simulated clouds. By including ionising radiation, we can reproduce the observed SFEs in the clouds most similar to nearby Milky Way clouds. For denser clouds, the SFE can approach unity. These observed SFEs are typically 3 to 10 times larger than the "total" SFEs, i.e. the fraction of the initial cloud mass converted to stars. Converting observed to total SFEs is non-trivial. We suggest some techniques for doing so, though estimate up to a factor of ten error in the conversion.

We investigate the conditions of ideal magnetohydrodynamic (MHD) turbulence responsible for the relative orientation between density structures, characterized by their gradient, $\vec{\nabla}\rho$, and the magnetic field, $\vec{B}$, in molecular clouds (MCs). For that purpose, we construct an expression for the time evolution of the angle, $\phi$, between $\vec{\nabla}\rho$ and $\vec{B}$ based on the transport equations of MHD turbulence. Using this expression, we find that the configuration where $\vec{\nabla}\rho$ and $\vec{B}$ are mostly parallel, $\cos\phi=1$, and where $\vec{\nabla}\rho$ and $\vec{B}$ are mostly perpendicular, $\cos\phi=0$, constitute attractors, that is, the system tends to evolve towards either of these configurations and they are more represented than others. This fact would explain the predominant alignment or anti-alignment between column density, $N_H$, structures and the projected magnetic field orientation, $\hat{B}_\perp$, reported in observations. Additionally, we find that departures from the $\cos\phi=0$ configurations are related to convergent flows, quantified by the divergence of the velocity field, $\vec{\nabla}\cdot\vec{v}$, in the presence of a relatively strong magnetic field. This would explain the observed change in relative orientation between $N_H$-structures and $\hat{B}_\perp$ towards MCs, from mostly parallel at low $N_H$ to mostly perpendicular at the highest $N_H$, as the result of the gravitational collapse and/or convergence of flows. Finally, we show that the density threshold that marks the observed change in relative orientation towards MCs, from $N_H$ and $\hat{B}_\perp$ being mostly parallel at low $N_H$ to mostly perpendicular at the highest $N_H$, is related to the magnetic field strength and constitutes a crucial piece of information for determining the role of the magnetic field in the dynamics of MCs.

Galaxy evolution and star formation are two multi-scale problems tightly linked to each other. To understand the interstellar cycle, which triggers galaxy evolution, it is necessary to describe simultaneously the large-scale evolution widely induced by the feedback processes and the details of the gas dynamics that controls the star formation process through gravitational collapse. We perform a set of three-dimensional high-resolution numerical simulations of a turbulent, self-gravitating and magnetized interstellar medium within a $1\ \mathrm{kpc}$ stratified box with supernova feedback correlated with star-forming regions. In particular, we focus on the role played by the magnetic field and the feedback on the galactic vertical structure, the star formation rate (SFR) and the flow dynamics. For this purpose we vary their respective intensities. We extract properties of the dense clouds arising from the turbulent motions and compute power spectra of various quantities. Using a distribution of supernovae sufficiently correlated with the dense gas, we find that supernova explosions can reproduce the observed SFR, particularly if the magnetic field is on the order of a few $\mu G$. The vertical structure, which results from a dynamical and an energy equilibrium is well reproduced by a simple analytical model, which allows us to estimate the coupling between the gas and the supernovae. We found the coupling to be rather low and on the order of 1.5$\%$. Strong magnetic fields may help to increase this coupling by a factor of about 2-3. To characterize the flow we compute the power spectra of various quantities in 3D but also in 2D in order to account for the stratification of the galactic disc.

Stellar clusters form by gravitational collapse of turbulent molecular clouds, with up to several thousand stars per cluster. They are thought to be the birthplace of most stars and therefore play an important role in our understanding of star formation, a fundamental problem in astrophysics. The initial conditions of the molecular cloud establish its dynamical history until the stellar cluster is born. However, the evolution of the cloud's angular momentum during cluster formation is not well understood. Current observations have suggested that turbulence scrambles the angular momentum of the cluster-forming cloud, preventing spin alignment amongst stars within a cluster. Here we use asteroseismology to measure the inclination angles of spin axes in 48 stars from the two old open clusters NGC~6791 and NGC~6819. The stars within each cluster show strong alignment. Three-dimensional hydrodynamical simulations of proto-cluster formation show that at least 50 % of the initial proto-cluster kinetic energy has to be rotational in order to obtain strong stellar-spin alignment within a cluster. Our result indicates that the global angular momentum of the cluster-forming clouds was efficiently transferred to each star and that its imprint has survived after several gigayears since the clusters formed.

We statistically evaluate the relative orientation between gas column density structures, inferred from Herschel submillimetre observations, and the magnetic field projected on the plane of sky, inferred from polarized thermal emission of Galactic dust observed by BLASTPol at 250, 350, and 500 micron, towards the Vela C molecular complex. First, we find very good agreement between the polarization orientations in the three wavelength-bands, suggesting that, at the considered common angular resolution of 3.0 arcminutes that corresponds to a physical scale of approximately 0.61 pc, the inferred magnetic field orientation is not significantly affected by temperature or dust grain alignment effects. Second, we find that the relative orientation between gas column density structures and the magnetic field changes progressively with increasing gas column density, from mostly parallel or having no preferred orientation at low column densities to mostly perpendicular at the highest column densities. This observation is in agreement with previous studies by the Planck collaboration towards more nearby molecular clouds. Finally, we find a correspondence between the trends in relative orientation and the shape of the column density probability distribution functions. In the sub-regions of Vela C dominated by one clear filamentary structure, or "ridges", we find a sharp transition from preferentially parallel or having no preferred relative orientation at low column densities to preferentially perpendicular at highest column densities. In the sub-regions of Vela C dominated by several filamentary structures with multiple orientations, or "nests", such a transition is also present, but it is clearly less sharp than in the ridge-like sub-regions. Both of these results suggest that the magnetic field is dynamically important for the formation of density structures in this region.

Large-scale shocks formed by clustered feedback of young OB stars are considered an important source of mechanical energy for the ISM and a trigger of molecular cloud formation. Their interaction sites are locations where kinetic energy and magnetic fields are redistributed between ISM phases. In this work we study the effect of the magnetic field on the expansion and fragmentation of supershells and look for the signatures of supershell collisions on dense structures and on the kinetic and magnetic energy distribution of the ISM. We performed a series of high-resolution, three-dimensional simulations of colliding supershells. The shocks are created by time-dependent feedback and evolve in a diffuse turbulent environment that is either unmagnetized or has different initial magnetic field configurations. In the hydrodynamical situation, the expansion law of the superbubbles is consistent with the radius-time relation that is theoretically predicted for wind-blown bubbles. The supershells fragment over their entire surface into small dense clumps that carry more than half of the total kinetic energy in the volume. However, this is not the case when a magnetic field is introduced, either in the direction of the collision or perpendicular to the collision. In magnetized situations, the shell surfaces are more stable to dynamical instabilities. When the magnetic field opposes the collision, the expansion law of the supershells also becomes significantly flatter than in the hydrodynamical case. Although a two-phase medium arises in all cases, in the MHD simulations the cold phase is limited to lower densities.

This paper assesses the roles of the presence of warm H2, and the increased formation rate due to the ion-neutral drift. We performed ideal MHD simulations that include the heating and cooling of the multiphase ISM, and where we treat dynamically the formation of H2. In a post-processing step we compute the abundances of species at chemical equilibrium. We show that CH+ is efficiently formed at the edge of clumps, in regions where the H2 fraction is low, but nevertheless higher than its equilibrium value, and where the gas temperature is high. We show that warm and out of equilibrium H2 increases the integrated column densities of CH+ by one order of magnitude, up to values still 3-10 times lower than those observed in the diffuse ISM. We balance the Lorentz force with the ion-neutral drag to estimate the ion-drift velocities (vd). We find that the vd distribution peaks around 0.04 km s-1, and that high vd are too rare to have a significant statistical impact on the abundances of CH+. Compared to previous works, our multiphase simulations reduce the spread in vd, and our self-consistent treatment of the ionisation leads to much reduced vd. Nevertheless, our resolution study shows that this velocity distribution is not converged: the ion-neutral drift has a higher impact on CH+ at higher resolution. On the other hand, our ideal MHD simulations do not include ambipolar diffusion, which would yield lower drift velocities. Within these limitations, we conclude that warm H2 is a key ingredient in the efficient formation of CH+ and that the ambipolar diffusion has very little influence on the abundance of CH+, mainly due to the small drift velocities obtained. However, we point out that small-scale processes and other non-thermal processes not included in our MHD simulation may be of crucial importance, and higher resolution studies with better controlled dissipation processes are needed.

Understanding how accretion proceeds in proto-planetary discs and more generally their dynamics is a crucial issue for explaining the conditions in which planets form. The role that accretion of gas from the surrounding molecular cloud onto the disc may have on its structure needs to be quantified. We perform tri-dimensional simulations using the Cartesian AMR code RAMSES of an accretion disc subject to infalling material. For the aspect ratio of $H/R \simeq 0.15$ and disk mass $M_d \simeq 10^{-2}$ M$_\odot$ used in our study, we find that for typical accretion rates on the order of a few 10$^{-7}$ M$_\odot$ yr$^{-1}$, values of the $\alpha$ parameter as high as a few 10$^{-3}$ are inferred. The mass that is accreted in the inner part of the disc is typically at least $50\%$ of the total mass that has been accreted onto the disc. Our results suggest that external accretion of gas at moderate values, onto circumstellar discs may trigger prominent spiral arms, reminiscent of recent observations made with various instruments, and lead to significant transport through the disc. If confirmed from observational studies, such accretion may therefore influence disc evolution.

Fast jets are thought to be a crucial ingredient of star formation because they might extract angular momentum from the disk and thus allow mass accretion onto the star. However, it is unclear whether jets are ubiquitous, and likewise, their contribution to mass and angular momentum extraction during protostar formation remains an open question. Our aim is to investigate the ejection process in the low-mass Class 0 protostar L1157. This source is associated with a spectacular bipolar outflow, and the recent detection of high-velocity SiO suggests the occurrence of a jet. Observations of CO 2-1 and SiO 5-4 at 0.8" resolution were obtained with the IRAM Plateau de Bure Interferometer as part of the CALYPSO large program. The jet and outflow structure were fit with a precession model. We derived the column density of CO and SiO, as well as the jet mass-loss rate and mechanical luminosity. High-velocity CO and SiO emission resolve for the first time the first 200 au of the outflow-driving molecular jet. The jet is strongly asymmetric, with the blue lobe 0.65 times slower than the red lobe. This suggests that the large-scale asymmetry of the outflow is directly linked to the jet velocity and that the asymmetry in the launching mechanism has been at work for the past 1800 yr. Velocity asymmetries are common in T Tauri stars, which suggests that the jet formation mechanism from Class 0 to Class II stages might be similar. Our model simultaneously fits the inner jet and the clumpy 0.2 pc scale outflow by assuming that the jet precesses counter-clockwise on a cone inclined by 73 degree to the line of sight with an opening angle of 8 degree on a period of 1640 yr. The estimated jet mass flux and mechanical luminosity are 7.7e-7 Msun/yr, and 0.9 Lsun, indicating that the jet could extract at least 25% of the gravitational energy released by the forming star.

The formation of protoplanetary discs during the collapse of molecular dense cores is significantly influenced by angular momentum transport, notably by the magnetic torque. In turn, the evolution of the magnetic field is determined by dynamical processes and non-ideal MHD effects such as ambipolar diffusion. Considering simple relations between various timescales characteristic of the magnetized collapse, we derive an expression for the early disc radius, $ r \simeq 18 \, {\rm AU} \, \left({\eta_{\rm AD} / 0.1 \, {\rm s}} \right)^{2/9} \left({B_z / 0.1\, {\rm G}} \right) ^{-4/9} \left({M / 0.1 \msol} \right) ^{1/3},$ where $M$ is the total disc plus protostar mass, $\eta_\mathrm{AD}$ is the ambipolar diffusion coefficient and $B_z$ is the magnetic field in the inner part of the core. This is about significantly smaller than the discs that would form if angular momentum was conserved. The analytical predictions are confronted against a large sample of 3D, non-ideal MHD collapse calculations covering variations of a factor 100 in core mass, a factor 10 in the level of turbulence, a factor 5 in rotation, and magnetic mass-to-flux over critical mass-to-flux ratios 2 and 5. The disc radius estimates are found to agree with the numerical simulations within less than a factor 2. A striking prediction of our analysis is the weak dependence of circumstellar disc radii upon the various relevant quantities, suggesting weak variations among class-0 disc sizes. In some cases, we note the onset of large spiral arms beyond this radius.

Molecular cloud structure is regulated by stellar feedback in various forms. Two of the most important feedback processes are UV photoionisation and supernovae from massive stars. However, the precise response of the cloud to these processes, and the interaction between them, remains an open question. In particular, we wish to know under which conditions the cloud can be dispersed by feedback, which in turn can give us hints as to how feedback regulates the star formation inside the cloud. We perform a suite of radiative magnetohydrodynamic simulations of a 10^5 solar mass cloud with embedded sources of ionising radiation and supernovae, including multiple supernovae and a hypernova model. A UV source corresponding to 10% of the mass of the cloud is required to disperse the cloud, suggesting that the star formation efficiency should be on the order of 10%. A single supernova is unable to significantly affect the evolution of the cloud. However, energetic hypernovae and multiple supernovae are able to add significant quantities of momentum to the cloud, approximately 10^{43} g cm/s of momentum per 10^{51} ergs of supernova energy. This is on the lower range of estimates in other works, since dense gas clumps that remain embedded inside the HII region cause rapid cooling in the supernova blast. We argue that supernovae alone are unable to regulate star formation in molecular clouds, and that strong pre-supernova feedback is required to allow supernova blastwaves to propagate efficiently into the interstellar medium

We develop a detailed chemical network relevant to the conditions characteristic of prestellar core collapse. We solve the system of time-dependent differential equations to calculate the equilibrium abundances of molecules and dust grains, with a size distribution given by size-bins for these latter. These abundances are used to compute the different non-ideal magneto-hydrodynamics resistivities (ambipolar, Ohmic and Hall), needed to carry out simulations of protostellar collapse. For the first time in this context, we take into account the evaporation of the grains, the thermal ionisation of Potassium, Sodium and Hydrogen at high temperature, and the thermionic emission of grains in the chemical network, and we explore the impact of various cosmic ray ionisation rates. All these processes significantly affect the non-ideal magneto-hydrodynamics resistivities, which will modify the dynamics of the collapse. Ambipolar diffusion and Hall effect dominate at low densities, up to n_H = 10^12 cm^-3, after which Ohmic diffusion takes over. We find that the time-scale needed to reach chemical equilibrium is always shorter than the typical dynamical (free fall) one. This allows us to build a large, multi-dimensional multi-species equilibrium abundance table over a large temperature, density and ionisation rate ranges. This table, which we make accessible to the community, is used during first and second prestellar core collapse calculations to compute the non-ideal magneto-hydrodynamics resistivities, yielding a consistent dynamical-chemical description of this process.

We present a new set of analytic models for the expansion of HII regions powered by UV photoionisation from massive stars and compare them to a new suite of radiative magnetohydrodynamic simulations of turbulent, self-gravitating molecular clouds. To perform these simulations we use the Eulerian adaptive mesh magnetohydrodynamics code RAMSES-RT, including radiative transfer of UV photons. Our analytic models successfully predict the global behaviour of the HII region provided the density and velocity structure of the cloud is known. We give estimates for the HII region behaviour based on a power law fit to the density field assuming that the system is virialised. We give a radius at which the ionisation front should stop expanding ("stall"). If this radius is smaller than the distance to the edge of the cloud, the HII region will be trapped by the cloud. This effect is more severe in collapsing clouds than in virialised clouds, since the density in the former increases dramatically over time, with much larger photon emission rates needed for the HII region to escape a collapsing cloud. We also measure the response of Jeans unstable gas to the HII regions to predict the impact of UV radiation on star formation in the cloud. We find that the mass in unstable gas can be explained by a model in which the clouds are evaporated by UV photons, suggesting that the net feedback on star formation should be negative

Within four nearby (d < 160 pc) molecular clouds, we statistically evaluate the structure of the interstellar magnetic field, projected on the plane of the sky and integrated along the line of sight, as inferred from the polarized thermal emission of Galactic dust observed by Planck at 353 GHz and from the optical and NIR polarization of background starlight. We compare the dispersion of the field orientation directly in vicinities with an area equivalent to that subtended by the Planck effective beam at 353 GHz (10') and using the second-order structure functions of the field orientation angles. We find that the average dispersion of the starlight-inferred field orientations within 10'-diameter vicinities is less than 20 deg, and that at these scales the mean field orientation is on average within 5 deg of that inferred from the submillimetre polarization observations in the considered regions. We also find that the dispersion of starlight polarization orientations and the polarization fractions within these vicinities are well reproduced by a Gaussian model of the turbulent structure of the magnetic field, in agreement with the findings reported by the Planck collaboration at scales greater than 10' and for comparable column densities. At scales greater than 10', we find differences of up to 14.7 deg between the second-order structure functions obtained from starlight and submillimetre polarization observations in the same positions in the plane of the sky, but comparison with a Gaussian model of the turbulent structure of the magnetic field indicates that these differences are small and are consistent with the difference in angular resolution between both techniques.

Stars are often observed to form in clusters. It is therefore important to understand how such a region of concentrated mass is assembled out of the diffuse medium and its properties eventually prescribe the important physical mechanisms and determine the characteristics of the stellar cluster. We study the formation of a gaseous proto-cluster inside a molecular cloud by performing high resolution MHD simulations and associate its internal properties to those of the parent cloud by varying the level of the initial turbulence of the cloud, with a view to better characterize the subsequent stellar cluster formation. The gaseous proto-cluster is formed out of global collapse of a molecular cloud, and has non-negligible rotation due to angular momentum conservation during the collapse of the object. Most of the star formation occurs in this region which occupies only a small volume fraction of the whole cloud. We identify such regions in simulations and compare the gas and sink particles to observations. The gaseous proto-cluster inferred from simulation results present a mass-size relation that is compatible with observations. We stress that the stellar cluster radius, although clearly correlated with the gas cluster radius, depends sensitively on its definition. Energy analysis is performed to confirm that the gaseous proto-cluster is a product of gravo-turbulent reprocessing and that the support of turbulent and rotational energy against self-gravity yields a state of global virial equilibrium although collapse is occurring at smaller scale and the cluster is forming stars actively. This object then serves as the antecedent of the stellar cluster, to which the energy properties are passed on.

Most stars are born in the gaseous proto-cluster environment. The knowledge of this intermediate stage gives more accurate constraints on star formation characteristics. We demonstrate that a virialized globally supported structure, in which star formation happens, is formed out of a collapsing molecular cloud, and derive a mapping from the parent cloud parameters to the proto-cluster to predict its properties, with a view to confront analytical calculations with observations and simulations. The virial theorem is decomposed into two dimensions to account for the rotation and the flattened geometry. Equilibrium is found by balancing rotation, turbulence and self-gravity, while turbulence is maintained by accretion driving and dissipates in one crossing time. The angular momentum and the accretion rate of the proto-cluster are estimated from the parent cloud properties. The two-dimensional virial model predicts the size and velocity dispersion given the mass of the proto-cluster and that of the parent cloud. The gaseous proto-clusters lie on a sequence of equilibrium with the trend $R \sim M^{0.5}$, with limited variations depending on the evolutionary stage, the parent cloud, and the parameters not well known like turbulence driving efficiency by accretion and the turbulence anisotropy. The model reproduces successfully observations and simulation results. The properties of proto-clusters follow universal relations and they can be derived from that of the parent cloud. Using simple estimates to infer the peak position of the core mass function (CMF) we find a weak dependence on the cluster mass suggesting that the physical conditions inside proto-clusters may contribute to set a CMF, and by extension an IMF, that looks independent of the environment.

Galactic cosmic rays (CR) are particles presumably accelerated in supernova remnant shocks that propagate in the interstellar medium up to the densest parts of molecular clouds, losing energy and their ionisation efficiency because of the presence of magnetic fields and collisions with molecular hydrogen. Recent observations hint at high levels of ionisation and at the presence of synchrotron emission in protostellar systems, which leads to an apparent contradiction. We want to explain the origin of these CRs accelerated within young protostars as suggested by observations. Our modelling consists of a set of conditions that has to be satisfied in order to have an efficient CR acceleration through diffusive shock acceleration. We analyse three main acceleration sites, then we follow the propagation of these particles through the protostellar system up to the hot spot region. We find that jet shocks can be strong accelerators of CR protons, which can be boosted up to relativistic energies. Other promising acceleration sites are protostellar surfaces, where shocks caused by impacting material during the collapse phase are strong enough to accelerate CR protons. In contrast, accretion flow shocks are too weak to efficiently accelerate CRs. Though CR electrons are weakly accelerated, they can gain a strong boost to relativistic energies through re-acceleration in successive shocks. We suggest a mechanism able to accelerate both CR protons and electrons through the diffusive shock acceleration mechanism, which can be used to explain the high ionisation rate and the synchrotron emission observed towards protostellar sources. The existence of an internal source of energetic particles can have a strong and unforeseen impact on the ionisation of the protostellar disc, on the star and planet formation processes, and on the formation of pre-biotic molecules.

The filamentary structure of the molecular interstellar medium and the potential link of this morphology to star formation have been brought into focus recently by high resolution observational surveys. An especially puzzling matter is that local interstellar filaments appear to have the same thickness, independent of their column density. This requires a theoretical understanding of their formation process and the physics that governs their evolution. In this work we explore a scenario in which filaments are dissipative structures of the large-scale interstellar turbulence cascade and ion-neutral friction (also called ambipolar diffusion) is affecting their sizes by preventing small-scale compressions. We employ high-resolution, 3D MHD simulations, performed with the grid code RAMSES, to investigate non-ideal MHD turbulence as a filament formation mechanism. We focus the analysis on the mass and thickness distributions of the resulting filamentary structures. Simulations of both driven and decaying MHD turbulence show that the morphologies of the density and the magnetic field are different when ambipolar diffusion is included in the models. In particular, the densest structures are broader and more massive as an effect of ion-neutral friction and the power spectra of both the velocity and the density steepen at a smaller wavenumber. The comparison between ideal and non-ideal MHD simulations shows that ambipolar diffusion causes a shift of the filament thickness distribution towards higher values. However, none of the distributions exhibit the pronounced peak found in the observed local filaments. Limitations in dynamical range and the absence of self-gravity in these numerical experiments do not allow us to conclude at this time whether this is due to the different filament selection or due to the physics inherent of the filament formation.

Accretion discs are ubiquitous in the universe and it is a crucial issue to understand how angular momentum and mass are being radially transported in these objects. Here, we study the role played by non-linear spiral patterns within hydrodynamical and non self-gravitating accretion disc assuming that external disturbances such as infall onto the disc may trigger them. To do so, we computed self-similar solutions that describe discs in which a spiral wave propagates. Such solutions present both shocks and critical sonic points that we carefully analyze. For all allowed temperatures and for several spiral shocks, we calculated the wave structure. In particular we inferred the angle of the spiral patern, the stress it exerts on the disc as well as the associated flux of mass and angular momentum as a function of temperature. We quantified the rate of angular momentum transport by means of the dimensionless $\alpha$ parameter. For the thickest disc we considered (corresponding to $h/r$ values of about 1/3), we found values of $\alpha$ as high as $0.1$, and scaling with the temperature $T$ such that $\alpha \propto T^{3/2} \propto (h/r)^3$. The spiral angle scales with the temperature as $\arctan(r/h)$. The existence of these solutions suggests that perturbations occurring at disc outer boundaries, such as for example perturbations due to infall motions, can propagate deep inside the disc and therefore should not be ignored, even when considering small radii.

H2 is the simplest and the most abundant molecule in the ISM, and its formation precedes the formation of other molecules. Understanding the dynamical influence of the environment and the interplay between the thermal processes related to the formation and destruction of H2 and the structure of the cloud is mandatory to understand correctly the observations of H2. We perform high resolution MHD colliding flow simulations with the AMR code RAMSES in which the physics of H2 has been included. We compare the simulation results with various observations including the column densities of excited rotational levels. Due to a combination of thermal pressure, ram pressure and gravity, the clouds produced at the converging point of HI streams are highly inhomogeneous. H2 molecules quickly form in relatively dense clumps and spread into the diffuse interclump gas. This in particular leads to the existence of significant abundances of H2 in the diffuse and warm gas that lies in between clumps. Simulations and observations show similar trends, specially for the HI-to-H2 transition. The abundances of excited rotational levels, calculated at equilibrium in the simulations are very similar to the observed abundances inferred from FUSE results. This is a direct consequence of the presence of the H2 enriched diffuse and warm gas. Our simulations show that H2 rapidly forms in the dense clumps and, due to the complex structure of molecular clouds, quickly spreads at lower densities. Consequently a significant fraction of warm H2 exists in the low density gas. This warm H2 leads to column densities of excited rotational levels close to the observed ones likely revealing the complex intermix between the warm and the cold gas in molecular clouds. This suggests that the 2-phase structure of molecular clouds is an essential ingredient to fully understand molecular hydrogen in these objects.

We performed high-resolution, 3D MHD simulations and we compared to observations of translucent molecular clouds. We show that the observed populations of rotational levels of H2 can arise as a consequence of the multi-phase structure of the ISM.

Owing to the paucity of sub-arcsecond (sub)mm observations required to probe the innermost regions of newly forming protostars, several fundamental questions are still being debated, such as the existence and coevality of close multiple systems. We study the physical and chemical properties of the jets and protostellar sources in the NGC1333-IRAS4A proto-binary system using continuum emission and molecular tracers of shocked gas. We observed NGC1333-IRAS4A in the SiO(6-5), SO(6_5-5_4), and CO(2-1) lines and the continuum emission at 1.3, 1.4, and 3 mm using the IRAM Plateau de Bure Interferometer in the framework of the CALYPSO large program. We clearly disentangle for the first time the outflow emission from the two sources A1 and A2. The two protostellar jets have very different properties: the A1 jet is faster, has a short dynamical timescale (<10^3 yr), and is associated with H2 shocked emission, whereas the A2 jet, which dominates the large-scale emission, is associated with diffuse emission, bends, and emits at slower velocities. The observed bending of the A2 jet is consistent with the change of propagation direction observed at large scale and suggests jet precession on very short timescales (~200-600 yr). In addition, a chemically rich spectrum with emission from several COMs (e.g. HCOOH, CH3OCHO, CH3OCH3) is only detected towards A2. Finally, very high-velocity shocked emission (~50 km s^-1) is observed along the A1 jet. An LTE analysis shows that SiO, SO, and H2CO abundances in the gas phase are enhanced up to (3-4)x10^{-7}, (1.4-1.7)x10^{-6}, and (3-7.9)x10^{-7}, respectively. The intrinsic different properties of the jets and driving sources in NGC1333-IRAS4A suggest different evolutionary stages for the two protostars, with A1 being younger than A2, in a very early stage of star formation previous to the hot-corino phase.

The surroundings of HII regions can have a profound influence on their development, morphology, and evolution. This paper explores the effect of the environment on H II regions in the MonR2 molecular cloud. We aim to investigate the density structure of envelopes surrounding HII regions and to determine their collapse and ionisation expansion ages. The Mon R2 molecular cloud is an ideal target since it hosts an H II region association. Column density and temperature images derived from Herschel data were used together to model the structure of HII bubbles and their surrounding envelopes. The resulting observational constraints were used to follow the development of the Mon R2 ionised regions with analytical calculations and numerical simulations. The four hot bubbles associated with H II regions are surrounded by dense, cold, and neutral gas envelopes. The radial density profiles are reminiscent of those of low-mass protostellar envelopes. The inner parts of envelopes of all four HII regions could be free-falling because they display shallow density profiles. As for their outer parts, the two compact HII regions show a density profile, which is typical of the equilibrium structure of an isothermal sphere. In contrast, the central UCHii region shows a steeper outer profile, that could be interpreted as material being forced to collapse. The size of the heated bubbles, the spectral type of the irradiating stars, and the mean initial neutral gas density are used to estimate the ionisation expansion time, texp, 0.1Myr,for the dense UCHII and compact HII regions and 0.35 Myr for the extended HII region. The envelope transition radii between the shallow and steeper density profiles are used to estimate the time elapsed since the formation of the first proto stellar embryo, Tinf : 1Myr, for the ultra-compact, 1.5 / 3Myr for the compact, and greater than 6Myr for the extended HII regions.

We have implemented non-ideal Magneto-Hydrodynamics (MHD) effects in the Adaptive Mesh Refinement (AMR) code RAMSES, namely ambipolar diffusion and Ohmic dissipation, as additional source terms in the ideal MHD equations. We describe in details how we have discretized these terms using the adaptive Cartesian mesh, and how the time step is diminished with respect to the ideal case, in order to perform a stable time integration. We have performed a large suite of test runs, featuring the Barenblatt diffusion test, the Ohmic diffusion test, the C-shock test and the Alfven wave test. For the latter, we have performed a careful truncation error analysis to estimate the magnitude of the numerical diffusion induced by our Godunov scheme, allowing us to estimate the spatial resolution that is required to address non-ideal MHD effects reliably. We show that our scheme is second-order accurate, and is therefore ideally suited to study non-ideal MHD effects in the context of star formation and molecular cloud dynamics.

In this paper, we provide a more accurate description of the evolution of the magnetic flux redistribution during prestellar core collapse by including resistive terms in the magnetohydrodynamics (MHD) equations. We focus more particularly on the impact of ambipolar diffusion. We use the adaptive mesh refinement code RAMSES to carry out such calculations. The resistivities required to calculate the ambipolar diffusion terms were computed using a reduced chemical network of charged, neutral and grain species. The inclusion of ambipolar diffusion leads to the formation of a magnetic diffusion barrier in the vicinity of the core, preventing accumulation of magnetic flux in and around the core and amplification of the field above 0.1G. The mass and radius of the first Larson core remain similar between ideal and non-ideal MHD models. This diffusion plateau has crucial consequences on magnetic braking processes, allowing the formation of disk structures. Magnetically supported outflows launched in ideal MHD models are weakened when using non-ideal MHD. Contrary to ideal MHD misalignment between the initial rotation axis and the magnetic field direction does not significantly affect the results for a given mu, showing that the physical dissipation truly dominate over numerical diffusion. We demonstrate severe limits of the ideal MHD formalism, which yield unphysical behaviours in the long-term evolution of the system. This includes counter rotation inside the outflow, interchange instabilities, and flux redistribution triggered by numerical diffusion, none observed in non-ideal MHD. Disks with Keplerian velocity profiles form in all our non-ideal MHD simulations, with final mass and size which depend on the initial magnetisation. This ranges from a few 0.01 solar masses and 20-30 au for the most magnetised case (mu=2) to 0.2 solar masses and 40-80 au for a lower magnetisation (mu=5).

The main signature of the interaction between cosmic rays and molecular clouds is the high ionisation degree. This decreases towards the densest parts of a cloud, where star formation is expected, because of energy losses and magnetic effects. However recent observations hint to high levels of ionisation in protostellar systems, therefore leading to an apparent contradiction that could be explained by the presence of energetic particles accelerated within young protostars. Our modelling consists of a set of conditions that has to be satisfied in order to have an efficient particle acceleration through the diffusive shock acceleration mechanism. We find that jet shocks can be strong accelerators of protons which can be boosted up to relativistic energies. Another possibly efficient acceleration site is located at protostellar surfaces, where shocks caused by impacting material during the collapse phase are strong enough to accelerate protons. Our results demonstrate the possibility of accelerating particles during the early phase of a proto-Solar-like system and can be used as an argument to support available observations. The existence of an internal source of energetic particles can have a strong and unforeseen impact on the star and planet formation process as well as on the formation of pre-biotic molecules.

We numerically investigate the dynamics of a 2D non-magnetised protoplanetary disc surrounded by an inflow coming from an external envelope. We find that the accretion shock between the disc and the inflow is unstable, leading to the generation of large-amplitude spiral density waves. These spiral waves propagate over long distances, down to radii at least ten times smaller than the accretion shock radius. We measure spiral-driven outward angular momentum transport with 1e-4 < alpha < 1e-2 for an inflow accretion rate Mout>1e-8 Msun/yr. We conclude that the interaction of the disc with its envelope leads to long-lived spiral density waves and radial angular momentum transport with rates that cannot be neglected in young non-magnetised protostellar discs.

Within ten nearby (d < 450 pc) Gould Belt molecular clouds we evaluate statistically the relative orientation between the magnetic field projected on the plane of sky, inferred from the polarized thermal emission of Galactic dust observed by Planck at 353 GHz, and the gas column density structures, quantified by the gradient of the column density, $N_H$. The selected regions, covering several degrees in size, are analyzed at an effective angular resolution of 10' FWHM, thus sampling physical scales from 0.4 to 40 pc in the nearest cloud. The column densities in the selected regions range from $N_H \approx 10^{21}$ to $10^{23}$ cm$^{-2}$, and hence they correspond to the bulk of the molecular clouds. The relative orientation is evaluated pixel by pixel and analyzed in bins of column density using the novel statistical tool called "Histogram of Relative Orientations". Throughout this study, we assume that the polarized emission observed by Planck at 353 GHz is representative of the projected morphology of the magnetic field in each region, i.e., we assume a constant dust grain alignment efficiency, independent of the local environment. Within most clouds we find that the relative orientation changes progressively with increasing $N_H$, from preferentially parallel or having no preferred orientation to preferentially perpendicular. In simulations of magnetohydrodynamic turbulence in molecular clouds this trend in relative orientation is a signature of Alfv\'enic or sub-Alfv\'enic turbulence, implying that the magnetic field is significant for the gas dynamics at the scales probed by Planck. We compare the deduced magnetic field strength with estimates we obtain from other methods and discuss the implications of the Planck observations for the general picture of molecular cloud formation and evolution.

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.

We present methods and results from "21-cm Spectral Line Observations of Neutral Gas with the EVLA" (21-SPONGE), a large survey for Galactic neutral hydrogen (HI) absorption with the Karl G. Jansky Very Large Array (VLA). With the upgraded capabilities of the VLA, we reach median root-mean-square (RMS) noise in optical depth of $\sigma_{\tau}=9\times 10^{-4}$ per $0.42\rm\,km\,s^{-1}$ channel for the 31 sources presented here. Upon completion, 21-SPONGE will be the largest HI absorption survey with this high sensitivity. We discuss the observations and data reduction strategies, as well as line fitting techniques. We prove that the VLA bandpass is stable enough to detect broad, shallow lines associated with warm HI, and show that bandpass observations can be combined in time to reduce spectral noise. In combination with matching HI emission profiles from the Arecibo Observatory ($\sim3.5'$ angular resolution), we estimate excitation (or spin) temperatures ($\rm T_s$) and column densities for Gaussian components fitted to sightlines along which we detect HI absorption (30/31). We measure temperatures up to $\rm T_s\sim1500\rm\,K$ for individual lines, showing that we can probe the thermally unstable interstellar medium (ISM) directly. However, we detect fewer of these thermally unstable components than expected from previous observational studies. We probe a wide range in column density between $\sim10^{16}$ and $>10^{21}\rm\,cm^{-2}$ for individual HI clouds. In addition, we reproduce the trend between cold gas fraction and average $\rm T_s$ found by synthetic observations of a hydrodynamic ISM simulation by Kim et al. (2014). Finally, we investigate methods for estimating HI $\rm T_s$ and discuss their biases.

Understanding the formation of binary and multiple stellar systems largely comes down to studying the circumstances for the fragmentation of a condensing core during the first stages of the collapse. However, the probability of fragmentation and the number of fragments seem to be determined to a large degree by the initial conditions. In this work we study the fate of the linear perturbations of a homogeneous gas sphere both analytically and numerically. In particular, we investigate the stability of the well-known homologous solution that describes the collapse of a uniform spherical cloud. The difficulty of the mathematical singularity in the perturbation equations is surpassed here by explicitly introducing a weak shock next to the sonic point. In parallel, we perform adaptive mesh refinement (AMR) numerical simulations of the linear stages of the collapse and compared the growth rates obtained by each method. With this combination of analytical and numerical tools, we explore the behavior of both spherically symmetric and non-axisymmetric perturbations. The numerical experiments provide the linear growth rates as a function of the core's initial virial parameter and as a function of the azimuthal wave number of the perturbation. The overlapping regime of the numerical experiments and the analytical predictions is the situation of a cold and large cloud, and in this regime the analytically calculated growth rates agree very well with the ones obtained from the simulations. The use of a weak shock as part of the perturbation allows us to find a physically acceptable solution to the equations for a continuous range of growth rates. The numerical simulations agree very well with the analytical prediction for the most unstable cores, while they impose a limit of a virial parameter of 0.1 for core fragmentation in the absence of rotation.

Context. Molecular clouds are known to be turbulent and strongly affected by stellar feedback. Moreover, stellar feedback is believed to drive turbulence at large scales in galaxies. Aims. We study the role played by supernovae in molecular clouds and the influence of the magnetic field on this process. Methods. We performed three-dimensional numerical simulations of supernova explosions, in and near turbulent self-gravitating molecular clouds. In order to study the influence of the magnetic field, we performed both hydrodynamical and magnetohydrodynamical simulations. We also ran a series of simple uniform density medium simulations and developed a simple analytical model. Results. We find that the total amount of momentum that is delivered during supernova explosions typically varies by a factor of about 2, even when the gas density changes by 3 orders of magnitude. However, the amount of momentum delivered to the dense gas varies by almost a factor of 10 if the supernova explodes within or outside the molecular cloud. The magnetic field has little influence on the total amount of momentum injected by the supernova explosions but increases the momentum injected into the dense gas. Conclusions. Supernovae that explode inside molecular clouds remove a significant fraction of the cloud mass. Supernovae that explode outside have a limited influence on the cloud. It is thus essential to know sufficiently well the correlation between supernovae and the surrounding dense material in order to know whether supernovae can regulate star formation effectively.

Context. Ultraviolet radiation plays a crucial role in molecular clouds. Radiation and matter are tightly coupled and their interplay influences the physical and chemical properties of gas. In particular, modeling the radiation propagation requires calculating column densities, which can be numerically expensive in high-resolution multidimensional simulations. Aims. Developing fast methods for estimating column densities is mandatory if we are interested in the dynamical influence of the radiative transfer. In particular, we focus on the effect of the UV screening on the dynamics and on the statistical properties of molecular clouds. Methods. We have developed a tree-based method for a fast estimate of column densities, implemented in the adaptive mesh refinement code RAMSES. We performed numerical simulations using this method in order to analyze the influence of the screening on the clump formation. Results. We find that the accuracy for the extinction of the tree-based method is better than 10%, while the relative error for the column density can be much more. We describe the implementation of a method based on precalculating the geometrical terms that noticeably reduces the calculation time. To study the influence of the screening on the statistical properties of molecular clouds we present the probability distribution function (PDF) of gas and the associated temperature per density bin and the mass spectra for different density thresholds. Conclusions. The tree-based method is fast and accurate enough to be used during numerical simulations since no communication is needed between CPUs when using a fully threaded tree. It is then suitable to parallel computing. We show that the screening for far UV radiation mainly affects the dense gas, thereby favoring low temperatures and affecting the fragmentation.

We examine variations of the stellar initial mass function (IMF) in extreme environments within the formalism derived by Hennebelle \& Chabrier. We focus on conditions encountered in progenitors of massive early type galaxies and starburst regions. We show that, when applying the concept of turbulent Jeans mass as the characteristic mass for fragmentation in a turbulent medium, instead of the standard thermal Jeans mass for purely gravitational fragmentation, the peak of the IMF in such environments is shifted towards smaller masses, leading to a bottom-heavy IMF, as suggested by various observations. In very dense and turbulent environments, we predict that the high-mass tail of the IMF can become even steeper than the standard Salpeter IMF, with a limit for the power law exponent $\alpha\simeq -2.7$, in agreement with recent observational determinations. This steepening is a direct consequence of the high densities and Mach values in such regions but also of the time dependence of the fragmentation process, as incorporated in the Hennebelle-Chabrier theory. We provide analytical parametrizations of these IMFs in such environments, to be used in galaxy evolution calculations. We also calculate the star formation rates and the mass-to-light ratios expected under such extreme conditions and show that they agree well with the values inferred in starburst environments and massive high-redshift galaxies. This reinforces the paradigm of star formation as being a universal process, i.e. the direct outcome of gravitationally unstable fluctuations in a density field initially generated by large scale shock-dominated turbulence. This globally enables us to infer the variations of the stellar IMF and related properties for atypical galactic conditions.

We present a new algorithm, named Autonomous Gaussian Decomposition (AGD), for automatically decomposing spectra into Gaussian components. AGD uses derivative spectroscopy and machine learning to provide optimized guesses for the number of Gaussian components in the data, and also their locations, widths, and amplitudes. We test AGD and find that it produces results comparable to human-derived solutions on 21cm absorption spectra from the 21cm SPectral line Observations of Neutral Gas with the EVLA (21-SPONGE) survey. We use AGD with Monte Carlo methods to derive the HI line completeness as a function of peak optical depth and velocity width for the 21-SPONGE data, and also show that the results of AGD are stable against varying observational noise intensity. The autonomy and computational efficiency of the method over traditional manual Gaussian fits allow for truly unbiased comparisons between observations and simulations, and for the ability to scale up and interpret the very large data volumes from the upcoming Square Kilometer Array and pathfinder telescopes.

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.

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.

Context: Star formation efficiency (SFE) theories are currently based on statistical distributions of turbulent cloud structures and a simple model of star formation from cores. They remain poorly tested, especially at the highest densities. Aims: We investigate the effects of gas density on the SFE through measurements of the core formation efficiency (CFE). With a total mass of $\sim2\times10^4$ M$_\odot$, the W43-MM1 ridge is one of the most convincing candidate precursor of starburst clusters and thus one of the best place to investigate star formation. Methods: We used high-angular resolution maps obtained at 3 mm and 1 mm within W43-MM1 with the IRAM Plateau de Bure Interferometer to reveal a cluster of 11 massive dense cores (MDCs), and, one of the most massive protostellar cores known. An Herschel column density image provided the mass distribution of the cloud gas. We then measured the 'instantaneous' CFE and estimated the SFE and the star formation rate (SFR) within subregions of the W43-MM1 ridge. Results: The high SFE found in the ridge ($\sim$6% enclosed in $\sim$8 pc$^3$) confirms its ability to form a starburst cluster. There is however a clear lack of dense cores in the northern part of the ridge, which may be currently assembling. The CFE and the SFE are observed to increase with volume gas density while the SFR steeply decreases with the virial parameter, $\alpha_{vir}$. Statistical models of the SFR may well describe the outskirts of the W43-MM1 ridge but struggle to reproduce its inner part, which corresponds to measurements at low $\alpha_{vir}$. It may be that ridges do not follow the log-normal density distribution, Larson relations, and stationary conditions forced in the statistical SFR models.

What exactly controls star formation in the Galaxy remains controversial. In particular, the role of feedback and magnetic field are still partially understood. We investigate the role played by supernovae feedback and magnetic field onto the star formation and the structure of the Galactic disk. We perform numerical simulations of the turbulent, magnetized, self-gravitating, multi-phase, supernovae regulated ISM within a 1 kpc stratified box. We implemented various schemes for the supernovae. This goes from a random distribution at a fixed rate to distributions for which the supernovae are spatially and temporally correlated to the formation of stars. To study the influence of magnetic field on star formation, we perform both hydrodynamical and magneto-hydrodynamical simulations. We find that supernovae feedback has a drastic influence on the galactic evolution. The supernovae distribution is playing a very significant role. When the supernovae are not correlated with star formation events, they do not modify significantly the very high star formation rate obtained without feedback. When the supernovae follow the accretion, the star formation rate can be reduced by a factor up to 30. Magnetic field is also playing a significant role. It reduces the star formation rate by a factor up to 2-3 and reduces the number of collapse sites by a factor of about 2. The exact correlation between the supernovae and the dense gas appears to have significant consequences on the galactic disk evolution and the star formation. This implies that small scale studies are necessary to understand and quantify the feedback efficiency. Magnetic field does influence the star formation at galactic scales by reducing the star formation rate and the number of star formation sites.

We review current theories for the origin of the Stellar Initial Mass Function (IMF) with particular focus on the extent to which the IMF can be considered universal across various environments. To place the issue in an observational context, we summarize the techniques used to determine the IMF for different stellar populations, the uncertainties affecting the results, and the evidence for systematic departures from universality under extreme circumstances. We next consider theories for the formation of prestellar cores by turbulent fragmentation and the possible impact of various thermal, hydrodynamic and magneto-hydrodynamic instabilities. We address the conversion of prestellar cores into stars and evaluate the roles played by different processes: competitive accretion, dynamical fragmentation, ejection and starvation, filament fragmentation and filamentary accretion flows, disk formation and fragmentation, critical scales imposed by thermodynamics, and magnetic braking. We present explanations for the characteristic shapes of the Present-Day Prestellar Core Mass Function and the IMF and consider what significance can be attached to their apparent similarity. Substantial computational advances have occurred in recent years, and we review the numerical simulations that have been performed to predict the IMF directly and discuss the influence of dynamics, time-dependent phenomena, and initial conditions.

In order to shed light on the main physical processes controlling fragmentation of massive dense cores, we present a uniform study of the density structure of 19 massive dense cores, selected to be at similar evolutionary stages, for which their relative fragmentation level was assessed in a previous work. We inferred the density structure of the 19 cores through a simultaneous fit of the radial intensity profiles at 450 and 850 micron (or 1.2 mm in two cases) and the Spectral Energy Distribution, assuming spherical symmetry and that the density and temperature of the cores decrease with radius following power-laws. We find a weak (inverse) trend of fragmentation level and density power-law index, with steeper density profiles tending to show lower fragmentation, and vice versa. In addition, we find a trend of fragmentation increasing with density within a given radius, which arises from a combination of flat density profile and high central density and is consistent with Jeans fragmentation. We considered the effects of rotational-to-gravitational energy ratio, non-thermal velocity dispersion, and turbulence mode on the density structure of the cores, and found that compressive turbulence seems to yield higher central densities. Finally, a possible explanation for the origin of cores with concentrated density profiles, which are the cores showing no fragmentation, could be related with a strong magnetic field, consistent with the outcome of radiation magnetohydrodynamic simulations.

We use the Karl G. Jansky Very Large Array (VLA) to conduct a high-sensitivity survey of neutral hydrogen (HI) absorption in the Milky Way. In combination with corresponding HI emission spectra obtained mostly with the Arecibo Observatory, we detect a widespread warm neutral medium (WNM) component with excitation temperature <Ts>= 7200 (+1800,-1200) K (68% confidence). This temperature lies above theoretical predictions based on collisional excitation alone, implying that Ly-{\alpha} scattering, the most probable additional source of excitation, is more important in the interstellar medium (ISM) than previously assumed. Our results demonstrate that HI absorption can be used to constrain the Ly-{\alpha} radiation field, a critical quantity for studying the energy balance in the ISM and intergalactic medium yet notoriously difficult to model because of its complicated radiative transfer, in and around galaxies nearby and at high redshift.

SPLASH (the Southern Parkes Large-Area Survey in Hydroxyl) is a sensitive, unbiased and fully-sampled survey of the Southern Galactic Plane and Galactic Centre in all four ground-state transitions of the hydroxyl (OH) radical. The survey provides a deep census of 1612-, 1665-, 1667- and 1720-MHz OH absorption and emission from the Galactic ISM, and is also an unbiased search for maser sources in these transitions. We present here first results from the SPLASH pilot region, which covers Galactic longitudes 334 to 344 degrees and latitudes of -2 to +2 degrees. Diffuse OH is widely detected in all four transitions, with optical depths that are always small (averaged over the Parkes beam), and with departures from LTE common even in the 1665- and 1667-MHz main lines. To a 3$\sigma$ sensitivity of 30 mK, we find no evidence of OH envelopes extending beyond the CO-bright regions of molecular cloud complexes, and conclude that the similarity of the OH excitation temperature and the level of the continuum background is at least partly responsible for this. We detect masers and maser candidates in all four transitions, approximately 50 per cent of which are new detections. This implies that SPLASH will produce a substantial increase in the known population of ground-state OH masers in the Southern Galactic Plane.

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.

Recent Herschel observations have confirmed that filaments are ubiquitous in molecular clouds and suggest that irrespectively of the column density, there is a characteristic width of about 0.1 pc whose physical origin remains unclear. We develop an analytical model that can be applied to self-gravitating accreting filaments. It is based on one hand on the virial equilibrium of the central part of the filament and on the other hand on energy balance between the turbulence driven by accretion onto the filament and dissipation. We consider two dissipation mechanisms the turbulent cascade and the ion-neutral friction. Our model predicts that the width of the filament inner part is almost independent of the column density and leads to values comparable to what is inferred observationally if dissipation is due to ion-neutral friction. On the contrary turbulent dissipation leads to a structure that is bigger and depends significantly on the column density. Our model provides a reasonable physical explanation which could explain the observed filament width when they are self-gravitating. It predicts the correct order or magnitude though hampered by some uncertainties.

SPIRou is a near-infrared (nIR) spectropolarimeter / velocimeter proposed as a new-generation instrument for CFHT. SPIRou aims in particular at becoming world-leader on two forefront science topics, (i) the quest for habitable Earth-like planets around very- low-mass stars, and (ii) the study of low-mass star and planet formation in the presence of magnetic fields. In addition to these two main goals, SPIRou will be able to tackle many key programs, from weather patterns on brown dwarf to solar-system planet atmospheres, to dynamo processes in fully-convective bodies and planet habitability. The science programs that SPIRou proposes to tackle are forefront (identified as first priorities by most research agencies worldwide), ambitious (competitive and complementary with science programs carried out on much larger facilities, such as ALMA and JWST) and timely (ideally phased with complementary space missions like TESS and CHEOPS). SPIRou is designed to carry out its science mission with maximum efficiency and optimum precision. More specifically, SPIRou will be able to cover a very wide single-shot nIR spectral domain (0.98-2.35 \mu m) at a resolving power of 73.5K, providing unpolarized and polarized spectra of low-mass stars with a ~15% average throughput and a radial velocity (RV) precision of 1 m/s.

We present a self-consistent hydrodynamical simulation of a Milky Way-like galaxy, at the resolution of 0.05 pc. The model includes star formation and a new implementation of stellar feedback through photo-ionization, radiative pressure and supernovae. The simulation resolves the structure of the interstellar medium at subparsec resolution for a few cloud lifetimes, and at 0.05 pc for about a cloud crossing time. Turbulence cascade and gravitation from the kpc scales are de facto included in smaller structures like molecular clouds. We show that the formation of a bar influences the dynamics of the central ~100 pc by creating resonances. At larger radii, the spiral arms host the formation of regularly spaced clouds: beads on a string and spurs. These instabilities pump turbulent energy into the gas, generally in the supersonic regime. Because of asymmetric drift, the supernovae explode outside of their gaseous nursery, which diminishes the effect of feedback on the structure of clouds. The evolution of clouds is thus mostly due to fragmentation and gas consumption, regulated mainly by supersonic turbulence. The transition from turbulence supported to self-gravitating gas is detected in the gas density probability distribution function at ~2000 cm^-3. The power spectrum density suggests that gravitation governs the hierarchical organisation of structures from the galactic scale down to a few parsecs.

We describe a morphological imprint of magnetization found when considering the relative orientation of the magnetic field direction with respect to the density structures in simulated turbulent molecular clouds. This imprint was found using the Histogram of Relative Orientations (HRO): a new technique that utilizes the gradient to characterize the directionality of density and column density structures on multiple scales. We present results of the HRO analysis in three models of molecular clouds in which the initial magnetic field strength is varied, but an identical initial turbulent velocity field is introduced, which subsequently decays. The HRO analysis was applied to the simulated data cubes and mock-observations of the simulations produced by integrating the data cube along particular lines of sight. In the 3D analysis we describe the relative orientation of the magnetic field $\mathbf{B}$ with respect to the density structures, showing that: 1.The magnetic field shows a preferential orientation parallel to most of the density structures in the three simulated cubes. 2.The relative orientation changes from parallel to perpendicular in regions with density over a critical density $n_{T}$ in the highest magnetization case. 3.The change of relative orientation is largest for the highest magnetization and decreases in lower magnetization cases. This change in the relative orientation is also present in the projected maps. In conjunction with simulations HROs can be used to establish a link between the observed morphology in polarization maps and the physics included in simulations of molecular clouds.

The relative importance of primordial molecular cloud fragmentation versus large-scale accretion still remains to be assessed in the context of massive core/star formation. Studying the kinematics of the dense gas surrounding massive-star progenitors can tell us the extent to which large-scale flow of material impacts the growth in mass of star-forming cores. Here we present a comprehensive dataset of the 5500(+/-800) Msun infrared dark cloud SDC335.579-0.272 (hereafter SDC335) which exhibits a network of cold, dense, parsec-long filaments. Atacama Large Millimeter Array (ALMA) Cycle 0 observations reveal two massive star-forming cores, MM1 and MM2, sitting at the centre of SDC335 where the filaments intersect. With a gas mass of 545(+770,-385) Msun contained within a source diameter of 0.05pc, MM1 is one of the most massive, compact protostellar cores ever observed in the Galaxy. As a whole, SDC335 could potentially form an OB cluster similar to the Trapezium cluster in Orion. ALMA and Mopra single-dish observations of the SDC335 dense gas furthermore reveal that the kinematics of this hub-filament system are consistent with a global collapse of the cloud. These molecular-line data point towards an infall velocity V_{inf} =0.7(+/-0.2) km/s, and a total mass infall rate \dot{M}_{inf} = 2.5(+/-1.0) x 10^{-3} Msun/yr towards the central pc-size region of SDC335. This infall rate brings 750(+/-300) Msun of gas to the centre of the cloud per free-fall time (t_{ff}=3x10^5 yr). This is enough to double the mass already present in the central pc-size region in 3.5(+2.2,-1.0) x t_{ff}. These values suggest that the global collapse of SDC335 over the past million year resulted in the formation of an early O-type star progenitor at the centre of the cloud's gravitational potential well.

{Filaments are ubiquitous in the interstellar medium as recently emphasized by Herschel, yet their physical origin remains elusive} {It is therefore important to understand the physics of molecular clouds to investigate how filaments form and what is the role played by various processes such as turbulence and magnetic field.} {We use ideal MHD simulations to study the formation of clumps in various conditions including different magnetization and Mach numbers as well as two completely different setup. We then perform several analysis to compute the shape of the clumps and their link to velocities and forces using various approaches.} {We find that on average, clumps in MHD simulations are more filamentary that clumps in hydrodynamical simulations. Detailed analyses reveal that the filaments are in general preferentially aligned with the strain which means that these structures simply result from the strech induced by turbulence. Moreover filaments tend to be confined by the Lorentz force which therefore lead them to survive longer in magnetized flows. We show that they have in all simulations a typical thickness equal to a few grid cells suggesting that they are primarily associated to the energy dissipation within the flow. We estimate the order of magnitude of the dissipation length associated to the ion-neutral friction and conclude that in well UV shielded regions it is of the order of 0.1 pc and therefore could possibly set the typical size of non self-gravitating filaments.} {Filaments are ubiquitous because they are the results of the very generic turbulent strain and because magnetic field help to keep them coherent. We suggest that energy dissipation is playing a determinant role in their formation.}

The present paper extends our previous theory of the stellar initial mass function (IMF) by including the time-dependence, and by including the impact of magnetic field. The predicted mass spectra are similar to the time independent ones with slightly shallower slopes at large masses and peak locations shifted toward smaller masses by a factor of a few. Assuming that star-forming clumps follow Larson type relations, we obtain core mass functions in good agreement with the observationally derived IMF, in particular when taking into account the thermodynamics of the gas. The time-dependent theory directly yields an analytical expression for the star formation rate (SFR) at cloud scales. The SFR values agree well with the observational determinations of various Galactic molecular clouds. Furthermore, we show that the SFR does not simply depend linearly on density, as sometimes claimed in the literature, but depends also strongly on the clump mass/size, which yields the observed scatter. We stress, however, that {\it any} SFR theory depends, explicitly or implicitly, on very uncertain assumptions like clump boundaries or the mass of the most massive stars that can form in a given clump, making the final determinations uncertain by a factor of a few. Finally, we derive a fully time-dependent model for the IMF by considering a clump, or a distribution of clumps accreting at a constant rate and thus whose physical properties evolve with time. In spite of its simplicity, this model reproduces reasonably well various features observed in numerical simulations of converging flows. Based on this general theory, we present a paradigm for star formation and the IMF.

[Abridged] Theoretical and numerical studies of star formation have shown that magnetic field (B) has a strong influence on both disk formation and fragmentation; even a relatively low B can prevent these processes. However, very few studies investigated the combined effects of B and turbulence. We study the effects of turbulence in magnetized core collapse, focusing on the magnetic diffusion, the orientation of the angular momentum (J) of the protostellar core, and on its consequences on disk formation, fragmentation and outflows. We perform 3D, AMR, MHD simulations of magnetically supercritical collapsing dense cores of 5 Msun using the MHD code RAMSES. A turbulent velocity field is imposed as initial conditions, characterised by a Kolmogorov power spectrum. Different levels of turbulence and magnetization are investigated, as well as 3 realisations for the turbulent velocity field. Magnetic diffusion, orientation of the rotation axis with respect to B, transport of J, disk formation, fragmentation and outflows formation are studied. The turbulent velocity field imposed as initial conditions contains a non-zero J, responsible for a misalignment of the rotation axis. Turbulence is also responsible for an effective turbulent diffusivity in the vicinity of the core. Both effects are responsible for a significant decrease of the magnetic braking, and facilitate the formation of early massive disks for not too high magnetization. Fragmentation can occur even with mu ~ 5 at late time in contrast with 1 Msun cores for which fragmentation is prevented for such values of mu. Slow asymmetric outflows are launched. They carry a mass which is comparable to the mass within the core. Because of misalignment and turbulent diffusion, massive disk formation is possible though their mass and size are still reduced compared to the hydrodynamical case. We find that for mu >= 5, fragmentation can happen.

This paper is a numerical study of the condensation of the warm neutral medium (WNM) into cold neutral medium (CNM) structures under the effect of turbulence and thermal instability. Using low resolution simulations we explored the impact of the WNM initial density and properties of the turbulence (stirring in Fourier with a varying mix of solenoidal and compressive modes) on the cold gas formation. Two sets of initial conditions which match the observations were selected to produce high resolution simulations (1024^3) allowing to study in details the properties of the produced dense structures. For typical values of the density, pressure and velocity dispersion of the WNM in the solar neighborhood, the turbulent motions of the HI can not provoque the phase transition from WNM to CNM, whatever their amplitude and their distribution in solenoidal and compressive modes. On the other hand we show that a quasi-isothermal increase in WNM density of a factor of 2 to 4 is enough to induce the phase transition, leading to the transition of about 40 percent of the gas to the cold phase within 1 Myr. Given the observed properties of the HI in the local ISM, the WNM and individual CNM structures in the local ISM are sub or transsonic and their dynamics are tightly interwoven. The velocity field bears the evidence of subsonic turbulence with a 2D power spectrum following the Kolmogorov law as P(k) \sim k^{-8/3} while the density is highly contrasted with a singificantly shallower power spectrum, reminiscent of what is observed in the cold ISM. Supra-thermal line width observed for CNM might be the result of relative velocity between cold structures. Finally, the cold structures denser than 5 cm^{-3} reproduce well the laws M \sim L^{2.25-2.28} and sigma(v) \sim 0.5-0.8L^{1/3} generally observed in molecular clouds.

In order to study the fragmentation of massive dense cores, which constitute the cluster cradles, we observed with the PdBI in the most extended configuration the continuum at 1.3 mm and the CO(2-1) emission of four massive cores. We detect dust condensations down to ~0.3 Msun and separate millimeter sources down to 0.4" or ~1000 AU, comparable to the sensitivities and separations reached in optical/infrared studies of clusters. The CO(2-1) high angular resolution images reveal high-velocity knots usually aligned with previously known outflow directions. This, in combination with additional cores from the literature observed at similar mass sensitivity and spatial resolution, allowed us to build a sample of 18 protoclusters with luminosities spanning 3 orders of magnitude. Among the 18 regions, ~30% show no signs of fragmentation, while 50% split up into ~4 millimeter sources. We compiled a list of properties for the 18 massive dense cores, such as bolometric luminosity, total mass, and mean density, and found no correlation of any of these parameters with the fragmentation level. In order to investigate the combined effects of magnetic field, radiative feedback and turbulence in the fragmentation process, we compared our observations to radiation magneto-hydrodynamic simulations, and found that the low-fragmented regions are well reproduced in the magnetized core case, while the highly-fragmented regions are consistent with cores where turbulence dominates over the magnetic field. Overall, our study suggests that the fragmentation in massive dense cores could be determined by the initial magnetic field/turbulence balance in each particular core.

Stars form within molecular clouds but our understanding of this fundamental process remains hampered by the complexity of the physics that drives their evolution. We review our observational and theoretical knowledge of molecular clouds trying to confront the two approaches wherever possible. After a broad presentation of the cold interstellar medium and molecular clouds, we emphasize the dynamical processes with special focus to turbulence and its impact on cloud evolution. We then review our knowledge of the velocity, density and magnetic fields. We end by openings towards new chemistry models and the links between molecular cloud structure and star--formation rates.

We briefly review the recent numerical works that have been performed to understand the formation of massive stars. After a brief description of the classical works, we review more specifically $i)$ the problem of building stars more massive than 20 $M_\odot$ and $ii)$ how to prevent the massive cores to fragment in many objects. Multi-D simulations succeed in circumventing the radiative pressure leading to the formation of massive stars although some questions are still debated regarding how is accretion exactly proceeding. While the core fragmentation is slightly reduced by the radiative feedback and the magnetic field when they are treated separately, it is almost entirely suppressed when both of them are included. This is because, magnetic field by removing angular momentum focusses the flow in a compact region. This makes the radiative feedback very efficient leading to a significant increase of the temperature.

It is believed that the majority of stars form in clusters. Therefore it is likely that the gas physical conditions that prevail in forming clusters, largely determine the properties of stars that form and in particular the initial mass function. We develop an analytical model to account for the formation of low mass clusters and the formation of stars within clusters. The formation of clusters is determined by an accretion rate, the virial equilibrium and energy as well as thermal balance. For this latter both molecular and dust cooling are considered using published rates. The star distribution is computed within the cluster using the physical conditions inferred from this model and the Hennebelle & Chabrier theory. Our model reproduces well the mass-size relation of low mass clusters (up to few $\simeq 10^3$ M$_\odot$ of stars corresponding to about 5 times more gas) and an initial mass function which is $i)$ very close to the Chabrier's IMF, $ii)$ weakly dependent on the mass of the clusters, $iii)$ relatively robust to (i.e. not too steeply dependent on) variations of physical quantities as accretion rate, radiation and cosmic rays abundances. The weak dependence of the mass distribution of stars with the cluster mass results from the compensation between varying clusters densities, velocity dispersions and temperatures all inferred from first physical principles. This constitutes a possible explanation for the apparent universality of the IMF within the Galaxy though variations with the local conditions could certainly be observed.

Theoretical studies of collapsing clouds have found that even a relatively weak magnetic field (B) may prevent the formation of disks and their fragmentation. However, most previous studies have been limited to cases where B and the rotation axis of the cloud are aligned. We study the transport of angular momentum, and its effects on disk formation, for non-aligned initial configurations and a range magnetic intensities. We perform 3D AMR MHD simulations of magnetically supercritical collapsing dense cores using the code Ramses. We compute the contributions of the processes transporting angular momentum (J), in the envelope and the region of the disk. We clearly define what could be defined as centrifugally supported disks and study their properties. At variance with earlier analyses, we show that the transport of J acts less efficiently in collapsing cores with non-aligned rotation axis and B. Analytically, this result can be understood by taking into account the bending of field lines occurring during the gravitational collapse. For the transport of J, we conclude that magnetic braking in the mean direction of B tends to dominate over both the gravitational and outflow transport of J. We find that massive disks, containing at least 10% of the initial core mass, can form during the earliest stages of star formation even for mass-to-flux ratios as small as 3 to 5 times the critical value. At higher field intensities, the early formation of massive disks is prevented. Given the ubiquity of Class I disks, and because the early formation of massive disks can take place at moderate magnetic intensities, we speculate that for stronger fields, disks will form later, when most of the envelope will have been accreted. In addition, we speculate that some observed early massive disks may actually be outflow cavities, mistaken for disks by projection effects. (Abridged version of the abstract.)

Our main purpose is to estimate the effect of assuming uniform density on the line-of-sight in PDR chemistry models, compared to a more realistic distribution for which total gas densities may well vary by several orders of magnitude. A secondary goal of this paper is to estimate the amount of molecular hydrogen which is not properly traced by the CO (J = 1 -> 0) line, the so-called "dark molecular gas". We use results from a magnetohydrodynamical (MHD) simulation as a model for the density structures found in a turbulent diffuse ISM with no star-formation activity. The Meudon PDR code is then applied to a number of lines of sight through this model, to derive their chemical structures. It is found that, compared to the uniform density assumption, maximal chemical abundances for H2, CO, CH and CN are increased by a factor 2 to 4 when taking into account density fluctuations on the line of sight. The correlations between column densities of CO, CH and CN with respect to those of H2 are also found to be in better overall agreement with observations. For instance, at N(H2) > 2.10^{20} cm-2, while observations suggest that d[log N(CO)]=d[log N(H2)] = 3.07 +/- 0.73, we find d[log N(CO)]=d[log N(H2)] =14 when assuming uniform density, and d[log N(CO)]=d[log N(H2)] = 5.2 when including density fluctuations.

We present a new publicly available tool (DustPol) aimed to model the polarised thermal dust emission. The module DustPol, which is publicly available, is part of the ARTIST (Adaptable Radiative Transfer Innovations for Submillimetre Telescopes) package, which also offers tools for modelling the polarisation of line emission together with a model library and a Python-based user interface. DustPol can easily manage analytical as well as pre-gridded models to generate synthetic maps of the Stokes I, Q, and U parameters. These maps are stored in FITS format which is straightforwardly read by the data reduction software used, e.g., by the Atacama Large Millimeter Array (ALMA). This turns DustPol into a powerful engine for the prediction of the expected polarisation features of a source observed with ALMA or the Planck satellite as well as for the interpretation of existing submillimetre observations obtained with other telescopes. DustPol allows the parameterisation of the maximum degree of polarisation and we find that, in a prestellar core, if there is depolarisation, this effect should happen at densities of 10^6 cm-3 or larger. We compare a model generated by DustPol with the observational polarisation data of the low-mass Class 0 object NGC 1333 IRAS 4A, finding that the total and the polarised emission are consistent.

These lectures attempt to expose the most important ideas, which have been proposed to explain the formation of stars with particular emphasis on the formation of brown dwarfs and low-mass stars. We first describe the important physical processes which trigger the collapse of a self-gravitating piece of fluid and regulate the star formation rate in molecular clouds. Then we review the various theories which have been proposed along the years to explain the origin of the stellar initial mass function paying particular attention to four models, namely the competitive accretion and the theories based respectively on stopped accretion, MHD shocks and turbulent dispersion. As it is yet unsettled whether the brown dwarfs form as low-mass stars, we present the theory of brown dwarfs based on disk fragmentation stressing all the uncertainties due to the radiative feedback and magnetic field. Finally, we describe the results of large scale simulations performed to explain the collapse and fragmentation of molecular clouds.

We present an analytical determination of the star formation rate (SFR) in molecular clouds, based on a time-dependent extension of our analytical theory of the stellar initial mass function (IMF). The theory yields SFR's in good agreement with observations, suggesting that turbulence {\it is} the dominant, initial process responsible for star formation. In contrast to previous SFR theories, the present one does not invoke an ad-hoc density threshold for star formation; instead, the SFR {\it continuously} increases with gas density, naturally yielding two different characteristic regimes, thus two different slopes in the SFR vs gas density relationship, in agreement with observational determinations. Besides the complete SFR derivation, we also provide a simplified expression, which reproduces reasonably well the complete calculations and can easily be used for quick determinations of SFR's in cloud environments. A key property at the heart of both our complete and simplified theory is that the SFR involves a {\it density-dependent dynamical time}, characteristic of each collapsing (prestellar) overdense region in the cloud, instead of one single mean or critical freefall timescale. Unfortunately, the SFR also depends on some ill determined parameters, such as the core-to-star mass conversion efficiency and the crossing timescale. Although we provide estimates for these parameters, their uncertainty hampers a precise quantitative determination of the SFR, within less than a factor of a few.

We report the results of radiation-magneto-hydrodynamics calculations in the context of high mass star formation, using for the first time a self-consistent model for photon emission (i.e. via thermal emission and in radiative shocks) and with the high resolution necessary to resolve properly magnetic braking effects and radiative shocks on scales <100 AU. We investigate the combined effects of magnetic field, turbulence, and radiative transfer on the early phases of the collapse and the fragmentation of massive dense cores. We identify a new mechanism that inhibits initial fragmentation of massive dense cores, where magnetic field and radiative transfer interplay. We show that this interplay becomes stronger as the magnetic field strength increases. Magnetic braking is transporting angular momentum outwards and is lowering the rotational support and is thus increasing the infall velocity. This enhances the radiative feedback owing to the accretion shock on the first core. We speculate that highly magnetized massive dense cores are good candidates for isolated massive star formation, while moderately magnetized massive dense cores are more appropriate to form OB associations or small star clusters.

We present a simple dimensional argument to illustrate the impact of nonthermal support from turbulent velocity dispersion on the shape of the prestellar core mass function (CMF), precursor of the stellar initial mass function (IMF). The argument demonstrates the need to invoke such support to recover the Salpeter slope in the high-mass part of the CMF/IMF, whereas pure thermal support leads to a much steeper slope. This simple dimensional argument clearly highlights the results obtained in the complete Hennebelle-Chabrier theory of the IMF.

Radiative transfer has a strong impact on the collapse and the fragmentation of prestellar dense cores. We present the radiation-hydrodynamics solver we designed for the RAMSES code. The method is designed for astrophysical purposes, and in particular for protostellar collapse. We present the solver, using the co-moving frame to evaluate the radiative quantities. We use the popular flux limited diffusion approximation, under the grey approximation (one group of photon). The solver is based on the second-order Godunov scheme of RAMSES for its hyperbolic part, and on an implicit scheme for the radiation diffusion and the coupling between radiation and matter. We report in details our methodology to integrate the RHD solver into RAMSES. We test successfully the method against several conventional tests. For validation in 3D, we perform calculations of the collapse of an isolated 1 M_sun prestellar dense core, without rotation. We compare successfully the results with previous studies using different models for radiation and hydrodynamics. We have developed a full radiation hydrodynamics solver in the RAMSES code, that handles adaptive mesh refinement grids. The method is a combination of an explicit scheme and an implicit scheme, accurate to the second-order in space. Our method is well suited for star formation purposes. Results of multidimensional dense core collapse calculations with rotation are presented in a companion paper.

Stars and more particularly massive stars, have a drastic impact on galaxy evolution. Yet the conditions in which they form and collapse are still not fully understood. In particular, the influence of the magnetic field on the collapse of massive clumps is relatively unexplored, it is thus of great relevance in the context of the formation of massive stars to investigate its impact. We perform high resolution, MHD simulations of the collapse of hundred solar masses, turbulent and magnetized clouds, using the adaptive mesh refinement code RAMSES. We compute various quantities such as mass distribution, magnetic field and angular momentum within the collapsing core and study the episodic outflows and the fragmentation that occurs during the collapse. The magnetic field has a drastic impact on the cloud evolution. We find that magnetic braking is able to substantially reduce the angular momentum in the inner part of the collapsing cloud. Fast and episodic outflows are being launched with typical velocities of the order of 3-5 km s$^{-1}$ although the highest velocities can be as high as 30-40 km s$^{-1}$. The fragmentation in several objects, is reduced in substantially magnetized clouds with respect to hydrodynamical ones by a factor of the order of 1.5-2. We conclude that magnetic fields have a significant impact on the evolution of massive clumps. In combination with radiation, magnetic fields largely determine the outcome of massive core collapse. We stress that numerical convergence of MHD collapse is a challenging issue. In particular, numerical diffusion appears to be important at high density therefore possibly leading to an over-estimation of the number of fragments.

We investigate the formation and evolution of giant molecular clouds (GMCs) by the collision of convergent warm neutral medium (WNM) streams in the interstellar medium, in the presence of magnetic fields and ambipolar diffusion (AD), focusing on the evolution of the star formation rate (SFR) and efficiency (SFE), as well as of the mass-to-magnetic-flux ratio (M2FR) in the forming clouds. We find that: 1) Clouds formed by supercritical inflow streams proceed directly to collapse, while clouds formed by subcritical streams first contract and then re-expand, oscillating on the scale of tens of Myr. 2) Our suite of simulations with initial magnetic field strength of 2, 3, and 4 $\mu\G$ show that only supercritical or marginal critical streams lead to reasonable star forming rates. 3) The GMC's M2FR is a generally increasing function of time, whose growth rate depends on the details of how mass is added to the GMC from the WNM. 4) The M2FR is a highly fluctuating function of position in the clouds. 5) In our simulations, the SFE approaches stationarity, because mass is added to the GMC at a similar rate at which it converts mass to stars. In such an approximately stationary regime, the SFE provides a proxy of the supercritical mass fraction in the cloud. 6) We observe the occurrence of buoyancy of the low-M2FR regions within the gravitationally-contracting GMCs, so that the latter naturally segregate into a high-density, high-M2FR "core" and a low-density, low-M2FR "envelope", without the intervention of AD. (Abridged)

We review the various theories which have been proposed along the years to explain the origin of the stellar initial mass function. We pay particular attention to four models, namely the competitive accretion and the theories based respectively on stopped accretion, MHD shocks and turbulent dispersion. In each case, we derive the main assumptions and calculations that support each theory and stress their respective successes and failures or difficulties.

Radiative feedback and magnetic field are understood to have a strong impact on the protostellar collapse. We present high resolution numerical calculations of the collapse of a 1 solar mass dense core in solid body rotation, including both radiative transfer and magnetic field. Using typical parameters for low-mass cores, we study thoroughly the effect of radiative transfer and magnetic field on the first core formation and fragmentation. We show that including the two aforementioned physical processes does not correspond to the simple picture of adding them separately. The interplay between the two is extremely strong, via the magnetic braking and the radiation from the accretion shock.

We present a simple statistical analysis of recent numerical simulations exploring the correlation between the core mass function obtained from the fragmentation of a molecular cloud and the stellar mass function which forms from these collapsing cores. Our analysis shows that the distributions of bound cores and sink particles obtained in the simulations are consistent with the sinks being formed predominantly from their parent core mass reservoir, with a statistical dispersion of the order of one third of the core mass. Such a characteristic dispersion suggests that the stellar initial mass function is relatively tightly correlated to the parent core mass function, leading to two similar distributions, as observed. This in turn argues in favor of the IMF being essentially determined at the early stages of core formation and being only weakly affected by the various environmental factors beyond the initial core mass reservoir, at least in the mass range explored in the present study. Accordingly, the final IMF of a star forming region should be determined reasonably accurately, statistically speaking, from the initial core mass function, provided some uniform efficiency factor. The calculations also show that these statistical fluctuations, due e.g. to variations among the core properties, broaden the low-mass tail of the IMF compared with the parent CMF, providing an explanation for the fact that this latter appears to underestimate the number of "pre brown dwarf" cores compared with the observationally-derived brown dwarf IMF.

We present 21-cm absorption measurements towards 12 radio continuum sources with previously identified thermally-unstable warm neutral medium (WNM). These observations were obtained with the Expanded Very Large Array (EVLA) and were complemented with the HI emission spectra obtained with the Arecibo Observatory. Out of 12 sources, HI absorption was detected along 5 lines of sight (seven new absorption features in total), resulting in a detection rate of ~42%. While our observations are sensitive to the WNM with a spin temperature T_s<3000 K, we detected only two wide absorption lines with T_s=400-900 K. These temperatures lie above the range allowed for the cold neutral medium (CNM) by the thermal equilbrium models and signify the thermally unstable WNM. Several absorption features have an optical depth of only a few x10^{-3}. While this is close or lower than what is theoretically expected for the CNM, we show that these weak lines are important for constraining the fraction of the thermally unstable WNM. Our observations demonstrate that, for the first time, high bandpass stability can be achieved with the VLA, allowing detection of absorption lines with a peak optical depth of ~10^{-3}.

Outflows and jets are intimately related to the formation of stars, and play an important role in redistributing mass, energy and angular momentum within the dense core and parent cloud. The interplay between magnetic field and rotation is responsible for launching these outflows, whose formation has been generally carried out for idealized systems where the angle $\alpha$ between the rotation axis and large-scale magnetic field is zero. Here we explore, through three-dimensional ideal magneto-hydrodynamic simulations, the effects of a non-zero $\alpha$ on the formation of outflows during the collapse of dense pre-stellar cores. We find that mass ejection is less efficient for increasing angle $\alpha$, and that outflows are essentially suppressed for $\alpha\sim90^{\circ}$. An important consequence is a corresponding increase of the mass accreted onto the adiabatic (first) core. In addition, mean flow velocities tend to increase with $\alpha$, and misaligned configurations produce clumpy, heterogeneous outflows that undergo precession, and are more prone to instabilities.

In contrast with numerous studies on the physical and chemical structure of low- and high-mass protostars, much less is known about their intermediate-mass counterparts, a class of objects that could help to elucidate the mechanisms of star formation on both ends of the mass range. We present the first results from a rich HIFI spectral dataset on an intermediate-mass protostar, OMC2-FIR4, obtained in the CHESS (Chemical HErschel SurveyS of star forming regions) key programme. The more than 100 methanol lines detected between 554 and 961 GHz cover a range in upper level energy of 40 to 540 K. Our physical interpretation focusses on the hot core, but likely the cold envelope and shocked regions also play a role in reality, because an analysis of the line profiles suggests the presence of multiple emission components. An upper limit of 10^-6 is placed on the methanol abundance in the hot core, using a population diagram, large-scale source model and other considerations. This value is consistent with abundances previously seen in low-mass hot cores. Furthermore, the highest energy lines at the highest frequencies display asymmetric profiles, which may arise from infall around the hot core.

In the past decade, much progress has been made in characterising the processes leading to the enhanced deuterium fractionation observed in the ISM and in particular in the cold, dense parts of star forming regions such as protostellar envelopes. Very high molecular D/H ratios have been found for saturated molecules and ions. However, little is known about the deuterium fractionation in radicals, even though simple radicals often represent an intermediate stage in the formation of more complex, saturated molecules. The imidogen radical NH is such an intermediate species for the ammonia synthesis in the gas phase. Herschel/HIFI represents a unique opportunity to study the deuteration and formation mechanisms of such species, which are not observable from the ground. We searched here for the deuterated radical ND in order to determine the deuterium fractionation of imidogen and constrain the deuteration mechanism of this species. We observed the solar-mass Class 0 protostar IRAS16293-2422 with the heterodyne instrument HIFI as part of the Herschel key programme CHESS (Chemical HErschel Surveys of Star forming regions). The deuterated form of the imidogen radical ND was detected and securely identified with 2 hyperfine component groups of its fundamental transition in absorption against the continuum background emitted from the nascent protostar. The 3 groups of hyperfine components of its hydrogenated counterpart NH were also detected in absorption. We derive a very high deuterium fractionation with an [ND]/[NH] ratio of between 30 and 100%. The deuterium fractionation of imidogen is of the same order of magnitude as that in other molecules, which suggests that an efficient deuterium fractionation mechanism is at play. We discuss two possible formation pathways for ND, by means of either the reaction of N+ with HD, or deuteron/proton exchange with NH.

In contrast to extensively studied dense star-forming cores, little is known about diffuse gas surrounding star-forming regions. We study molecular gas in the high-mass star-forming region NGC6334I, which contains diffuse, quiescent components that are inconspicuous in widely used molecular tracers such as CO. We present Herschel/HIFI observations of CH toward NGC6334I observed as part of the CHESS key program. HIFI resolves the hyperfine components of its J=3/2-1/2 transition, observed in both emission and absorption. The CH emission appears close to the systemic velocity of NGC6334I, while its measured linewidth of 3 km/s is smaller than previously observed in dense gas tracers such as NH3 and SiO. The CH abundance in the hot core is 7 10^-11, two to three orders of magnitude lower than in diffuse clouds. While other studies find distinct outflows in, e.g., CO and H2O toward NGC6334I, we do not detect outflow signatures in CH. To explain the absorption signatures, at least two absorbing components are needed at -3.0 and +6.5 km/s with N(CH)=7 10^13 and 3 10^13 cm^-2. Two additional absorbing clouds are found at +8.0 and 0.0 km/s, both with N(CH)=2 10^13 cm^-2. Turbulent linewidths for the four absorption components vary between 1.5 and 5.0 km/s in FWHM. We constrain physical properties of our CH clouds by matching our CH absorbers with other absorption signatures. In the hot core, molecules such as H2O and CO trace gas that is heated and dynamically influenced by outflow activity, whereas CH traces more quiescent material. The four CH absorbers have column densities and turbulent properties consistent with diffuse clouds: two are located near NGC6334, and two are unrelated foreground clouds. Local density and dynamical effects influence the chemical composition of physical components of NGC6334, causing some components to be seen in CH but not in other tracers, and vice versa.

Despite the low elemental deuterium abundance in the Galaxy, enhanced molecular D/H ratios have been found in the environments of low-mass star forming regions, and in particular the Class 0 protostar IRAS 16293-2422. The CHESS (Chemical HErschel Surveys of Star forming regions) Key Program aims at studying the molecular complexity of the interstellar medium. The high sensitivity and spectral resolution of the HIFI instrument provide a unique opportunity to observe the fundamental 1,1,1 - 0,0,0 transition of the ortho-D2O molecule, inaccessible from the ground, and to determine the ortho-to-para D2O ratio. We have detected the fundamental transition of the ortho-D2O molecule at 607.35 GHz towards IRAS 16293-2422. The line is seen in absorption with a line opacity of 0.62 +/- 0.11 (1 sigma). From the previous ground-based observations of the fundamental 1,1,0 - 1,0,1 transition of para-D2O seen in absorption at 316.80 GHz we estimate a line opacity of 0.26 +/- 0.05 (1 sigma). We show that the observed absorption is caused by the cold gas in the envelope of the protostar. Using these new observations, we estimate for the first time the ortho to para D2O ratio to be lower than 2.6 at a 3 sigma level of uncertainty, to be compared with the thermal equilibrium value of 2:1.

We present observations of twelve rotational transitions of H2O-16, H2O-18, and H2O-17 toward the massive star-forming region NGC 6334 I, carried out with Herschel/HIFI as part of the guaranteed time key program Chemical HErschel Surveys of Star forming regions (CHESS). We analyze these observations to obtain insights into physical processes in this region. We identify three main gas components (hot core, cold foreground, and outflow) in NGC 6334 I and derive the physical conditions in these components. The hot core, identified by the emission in highly excited lines, shows a high excitation temperature of 200 K, whereas water in the foreground component is predominantly in the ortho- and para- ground states. The abundance of water varies between 4 10^-5 (outflow) and 10^-8 (cold foreground gas). This variation is most likely due to the freeze-out of water molecules onto dust grains. The H2O-18/H2O-17 abundance ratio is 3.2, which is consistent with the O-18/O-17 ratio determined from CO isotopologues. The ortho/para ratio in water appears to be relatively low 1.6(1) in the cold, quiescent gas, but close to the equilibrium value of three in the warmer outflow material (2.5(0.8)).

The HIFI instrument on board the Herschel Space Observatory has been used to observe interstellar nitrogen hydrides along the sight-line towards G10.6-0.4 in order to improve our understanding of the interstellar chemistry of nitrogen. We report observations of absorption in NH N=1-0, J=2-1 and ortho-NH2 1_1,1-0_0,0. We also observed ortho-NH3 1_0-0_0, and 2_0-1_0, para-NH3 2_1-1_1, and searched unsuccessfully for NH+. All detections show emission and absorption associated directly with the hot-core source itself as well as absorption by foreground material over a wide range of velocities. All spectra show similar, non-saturated, absorption features, which we attribute to diffuse molecular gas. Total column densities over the velocity range 11-54 km/s are estimated. The similar profiles suggest fairly uniform abundances relative to hydrogen, approximately 6*10^-9, 3*10^-9, and 3*10^-9 for NH, NH2, and NH3, respectively. These abundances are discussed with reference to models of gas-phase and surface chemistry.

We discuss the detection of absorption by interstellar hydrogen fluoride (HF) along the sight line to the submillimeter continuum sources W49N and W51. We have used Herschel's HIFI instrument in dual beam switch mode to observe the 1232.4762 GHz J = 1 - 0 HF transition in the upper sideband of the band 5a receiver. We detected foreground absorption by HF toward both sources over a wide range of velocities. Optically thin absorption components were detected on both sight lines, allowing us to measure - as opposed to obtain a lower limit on - the column density of HF for the first time. As in previous observations of HF toward the source G10.6-0.4, the derived HF column density is typically comparable to that of water vapor, even though the elemental abundance of oxygen is greater than that of fluorine by four orders of magnitude. We used the rather uncertain N(CH)-N(H2) relationship derived previously toward diffuse molecular clouds to infer the molecular hydrogen column density in the clouds exhibiting HF absorption. Within the uncertainties, we find that the abundance of HF with respect to H2 is consistent with the theoretical prediction that HF is the main reservoir of gas-phase fluorine for these clouds. Thus, hydrogen fluoride has the potential to become an excellent tracer of molecular hydrogen, and provides a sensitive probe of clouds of small H2 column density. Indeed, the observations of hydrogen fluoride reported here reveal the presence of a low column density diffuse molecular cloud along the W51 sight line, at an LSR velocity of ~ 24kms-1, that had not been identified in molecular absorption line studies prior to the launch of Herschel.

We report the first detection of chloronium, H$_2$Cl$^+$, in the interstellar medium, using the HIFI instrument aboard the \emph{Herschel} Space Observatory. The $2_{12}-1_{01}$ lines of ortho-H$_2^{35}$Cl$^+$ and ortho-H$_2^{37}$Cl$^+$ are detected in absorption towards NGC~6334I, and the $1_{11}-0_{00}$ transition of para-H$_2^{35}$Cl$^+$ is detected in absorption towards NGC~6334I and Sgr~B2(S). The H$_2$Cl$^+$ column densities are compared to those of the chemically-related species HCl. The derived HCl/H$_2$Cl$^+$ column density ratios, $\sim$1--10, are within the range predicted by models of diffuse and dense Photon Dominated Regions (PDRs). However, the observed H$_2$Cl$^+$ column densities, in excess of $10^{13}$~cm$^{-2}$, are significantly higher than the model predictions. Our observations demonstrate the outstanding spectroscopic capabilities of HIFI for detecting new interstellar molecules and providing key constraints for astrochemical models.

We present Herschel/HIFI observations of the fundamental rotational transitions of ortho- and para-H$_2^{16}$O and H$_2^{18}$O in absorption towards Sagittarius~B2(M) and W31C. The ortho/para ratio in water in the foreground clouds on the line of sight towards these bright continuum sources is generally consistent with the statistical high-temperature ratio of 3, within the observational uncertainties. However, somewhat unexpectedly, we derive a low ortho/para ratio of $2.35 \pm 0.35$, corresponding to a spin temperature of $\sim$27~K, towards Sagittarius~B2(M) at velocities of the expanding molecular ring. Water molecules in this region appear to have formed with, or relaxed to, an ortho/para ratio close to the value corresponding to the local temperature of the gas and dust.

The formation of massive stars is a highly complex process in which it is not clear whether the star-forming gas is in global gravitational collapse or in an equilibrium state, supported by turbulence. By studying one of the most massive and dense star-forming regions in the Galaxy at a distance of less than 3 kpc, the filament containing the well-known sources DR21 and DR21(OH), we expect to find observational signatures that allow to discriminate between the two views. We use molecular line data from our 13CO 1-0, CS 2-1, and N2H+ 1-0 survey of the Cygnus X region obtained with the FCRAO and high-angular resolution observations of CO, CS, HCO+, N2H+, and H2CO, obtained with the IRAM 30m telescope. We observe a complex velocity field and velocity dispersion in the DR21 filament in which regions of highest column-density, i.e. dense cores, have a lower velocity dispersion than the surrounding gas and velocity gradients that are not (only) due to rotation. Infall signatures in optically thick line profiles of HCO+ and 12CO are observed along and across the whole DR21 filament. From modelling the observed spectra, we obtain a typical infall speed of 0.6 km/s and mass accretion rates of the order of a few 10^-3 Msun/yr for the two main clumps constituting the filament. These massive (4900 and 3300 Msun) clumps are both gravitationally contracting. All observed kinematic features in the DR21 filament can be explained if it is formed by the convergence of flows at large scales and is now in a state of global gravitational collapse. Whether this convergence of flows originated from self-gravity at larger scales or from other processes can not be settled with the present study. The observed velocity field and velocity dispersion are consistent with results from (magneto)-hydrodynamic simulations where the cores lie at the stagnation points of convergent turbulent flows.

We report the detection of absorption by interstellar hydroxyl cations and water cations, along the sight-line to the bright continuum source W49N. We have used Herschel's HIFI instrument, in dual beam switch mode, to observe the 972 GHz N = 1 - 0 transition of OH+ and the 1115 GHz 1(11) - 0(00) transition of ortho-H2O+. The resultant spectra show absorption by ortho-H2O+, and strong absorption by OH+, in foreground material at velocities in the range 0 to 70 km/s with respect to the local standard of rest. The inferred OH+/H2O+ abundance ratio ranges from ~ 3 to ~ 15, implying that the observed OH+ arises in clouds of small molecular fraction, in the 2 - 8% range. This conclusion is confirmed by the distribution of OH+ and H2O+ in Doppler velocity space, which is similar to that of atomic hydrogen, as observed by means of 21 cm absorption measurements, and dissimilar from that typical of other molecular tracers. The observed OH+/H abundance ratio of a few E-8 suggests a cosmic ray ionization rate for atomic hydrogen of (0.6 - 2.4) E-16 s-1, in good agreement with estimates inferred previously for diffuse clouds in the Galactic disk from observations of interstellar H3+ and other species.

We study the rotational properties of magnetized and self-gravitating molecular cloud cores formed in 2 very high resolution 3D molecular cloud simulations.The simulations have been performed using the code RAMSES at an effective resolution of 4096^3.One simulation represents a mildly magnetically-supercritical cloud and the other a strongly magnetically-supercritical cloud.A noticeable difference between the 2 simulations is the core formation efficiency (CFE) of the high density cores.In the strongly supercritical simulations the CFE is ~33 % per free-fall time of the cloud tff,cl, whereas in the mildly supercritical simulations this value goes down to ~6%/tff,cl. A comparison of the intrinsic specific angular momentum j3D distributions of the cores with the distribitions of j2D derived using synthetic 2D velocity maps of the cores,shows that the synthetic observations tend to overestimate the true value of j by a factor of ~10.The origin of this discrepancy lies in the fact that contrary to the intrinsic determination which sums up the individual gas parcels contributions to j, the determination of j using the observational procedure which is based on a measurement on the global velocity gradient under the hypothesis of uniform rotation smoothes out the complex fluctuations present in the 3D velocity field. Our results provide a natural explanation for the discrepancy by a factor ~10 observed between the intrinsic 3D distributions of j and the corresponding distributions derived in real observations.We suggest that measurements of j which are based on the measurement of the observed global velocity gradients may need to be reduced by a factor of ~10 in order to derive a more accurate estimate of j in the cores.

The outflow driven by the low-mass class 0 protostar L1157 is the prototype of the so-called chemically active outflows. The bright bowshock B1 in the southern outflow lobe is a privileged testbed of magneto-hydrodynamical (MHD) shock models, for which dynamical and chemical processes are strongly interdependent. We present the first results of the unbiased spectral survey of the L1157-B1 bowshock, obtained in the framework of the key program "Chemical Herschel Surveys of Star Forming Regions" (CHESS). The main aim is to trace the warm and chemically enriched gas and to infer the excitation conditions in the shock region. The CO 5-4 and H2O lines have been detected at high-spectral resolution in the unbiased spectral survey of the HIFI-Band 1b spectral window (555-636 GHz), presented by Codella et al. in this volume. Complementary ground-based observations in the submm window help establish the origin of the emission detected in the main-beam of HIFI, and the physical conditions in the shock.}{Both lines exhibit broad wings, which extend to velocities much higher than reported up to now. We find that the molecular emission arises from two regions with distinct physical conditions: an extended, warm (100K), dense (3e5 cm-3) component at low-velocity, which dominates the water line flux in Band~1; a secondary component in a small region of B1 (a few arcsec) associated with high-velocity, hot (> 400 K) gas of moderate density ((1.0-3.0)e4 cm-3), which appears to dominate the flux of the water line at 179mu observed with PACS. The water abundance is enhanced by two orders of magnitude between the low- and the high-velocity component, from 8e-7 up to 8e-5. The properties of the high-velocity component agree well with the predictions of steady-state C-shock models.

We report the detection of absorption lines by the reactive ions OH+, H2O+ and H3O+ along the line of sight to the submillimeter continuum source G10.6$-$0.4 (W31C). We used the Herschel HIFI instrument in dual beam switch mode to observe the ground state rotational transitions of OH+ at 971 GHz, H2O+ at 1115 and 607 GHz, and H3O+ at 984 GHz. The resultant spectra show deep absorption over a broad velocity range that originates in the interstellar matter along the line of sight to G10.6$-$0.4 as well as in the molecular gas directly associated with that source. The OH+ spectrum reaches saturation over most velocities corresponding to the foreground gas, while the opacity of the H2O+ lines remains lower than 1 in the same velocity range, and the H3O+ line shows only weak absorption. For LSR velocities between 7 and 50 kms$^{-1}$ we estimate total column densities of $N$(OH+) $> 2.5 \times 10^{14}$ cm$^{-2}$, $N$(H2O+) $\sim 6 \times 10^{13}$ cm$^{-2}$ and $N$(H3O+) $\sim 4.0 \times 10^{13}$ cm$^{-2}$. These detections confirm the role of O$^+$ and OH$^+$ in initiating the oxygen chemistry in diffuse molecular gas and strengthen our understanding of the gas phase production of water. The high ratio of the OH+ by the H2O+ column density implies that these species predominantly trace low-density gas with a small fraction of hydrogen in molecular form.

Complex turbulent motions are ubiquitously observed in many astrophysical systems. Their origin, however, is still poorly understood. When cosmic structures form, they grow in mass via accretion from the surrounding environment. We propose that this accretion is able to drive internal turbulent motions in a wide range of astrophysical objects and study this process in the case of galaxies, molecular clouds and protoplanetary disks. We use a combination of numerical simulations and analytical arguments to predict the level of turbulence as a function of the accretion rate, the dissipation scale, and the density contrast, and compare with observational data. We find that in Milky Way type galaxies the observed level of turbulence in the interstellar medium can be explained by accretion, provided that the galaxies gain mass at a rate comparable to the rate at which they form stars. This process is particularly relevant in the extended outer disks beyond the star-forming radius. We also calculate the rate at which molecular clouds grow in mass when they build up from the atomic component of the galactic gas and find that their internal turbulence is likely to be driven by accretion as well. It is the very process of cloud formation that excites turbulent motions on small scales by establishing the turbulent cascade. In the case of T Tauri disks, we show that accretion can drive subsonic turbulence at the observed level if the rate at which gas falls onto the disk is comparable to the rate at which disk material accretes onto the central star. This also explains the observed relation of accretion rate and stellar mass, dM/dt ~ M^1.8. The efficiency required to convert infall motion into turbulence is of the order of a few percent in all three cases. We conclude that accretion-driven turbulence is a universal concept with far-reaching implications for a wide range of astrophysical objects.

We report the detection of strong absorption by interstellar hydrogen fluoride along the sight-line to the submillimeter continuum source G10.6-0.4 (W31C). We have used Herschel's HIFI instrument, in dual beam switch mode, to observe the 1232.4763 GHz J=1-0 HF transition in the upper sideband of the Band 5a receiver. The resultant spectrum shows weak HF emission from G10.6-0.4 at LSR velocities in the range -10 to -3 km/s, accompanied by strong absorption by foreground material at LSR velocities in the range 15 to 50 km/s. The spectrum is similar to that of the 1113.3430 GHz 1(11)-0(00) transition of para-water, although at some frequencies the HF (hydrogen fluoride) optical depth clearly exceeds that of para-H2O. The optically-thick HF absorption that we have observed places a conservative lower limit of 1.6E+14 cm-2 on the HF column density along the sight-line to G10.6-0.4. Our lower limit on the HF abundance, 6E-9 relative to hydrogen nuclei, implies that hydrogen fluoride accounts for between ~ 30 and 100% of the fluorine nuclei in the gas phase along this sight-line. This observation corroborates theoretical predictions that - because the unique thermochemistry of fluorine permits the exothermic reaction of F atoms with molecular hydrogen - HF will be the dominant reservoir of interstellar fluorine under a wide range of conditions.

The formation process of binary stars and multiple systems is poorly understood. Here, we seek to determine the typical outcome of protostellar collapse and to constrain models of binary formation by core fragmentation during collapse, using high-resolution millimeter continuum imaging of very young (Class 0) protostars observed at the beginning of the main accretion phase. We carried out a pilot high-resolution study of 5 Class 0 objects, using the most extended (A) configuration of the IRAM PdBI at 1.3 mm, which allow us to probe the multiplicity of Class 0 protostars down to separations a ~50 AU and circumstellar mass ratios q ~0.07. We show that our PdBI observations revealed only wide (>1500 AU) protobinary systems and/or outflow-generated features. When combined with previous millimeter interferometric observations of Class 0 protostars, our pilot PdBI study tentatively suggests that the binary fraction in the ~ 75-1000 AU range increases from the Class 0 to the Class I stage. It also seems to argue against purely hydrodynamic models of binary star formation. We briefly discuss possible alternative scenarios to reconcile the low multiplicity rate of Class 0 protostars on small scales with the higher binary fraction observed at later evolutionary stages.

It is established that both radiative transfer and magnetic field have a strong impact on the collapse and the fragmentation of prestellar dense cores, but no consistent calculation exists yet at such scales. We present original AMR calculations including magnetic field (in the ideal MHD limit) and radiative transfer, within the Flux Limited Diffusion approximation, of the collapse of a 1 solar mass dense core. We compare the results with calculations performed with a barotropic EOS. We show that radiative transfer has an important impact on the collapse and the fragmentation, through the cooling or heating of the gas, and is complementary of the magnetic field. A larger field yields a stronger magnetic braking, increasing the accretion rate and thus the effect of the radiative feedback. Even for a strongly magnetized core, where the dynamics of the collapse is dominated by the magnetic field, radiative transfer is crucial to determine the temperature and optical depth distributions, two potentially accessible observational diagnostics. A barotropic EOS cannot account for realistic fragmentation. The diffusivity of the numerical scheme, however, is found to strongly affect the output of the collapse, leading eventually to spurious fragmentation. Both radiative transfer and magnetic field must be included in numerical calculations of star formation to obtain realistic collapse configurations and observable signatures. Nevertheless, the numerical resolution and the robustness of the solver are of prime importance to obtain reliable results. When using an accurate solver, the fragmentation is found to always remain inhibited by the magnetic field, at least in the ideal MHD limit, even when radiative transfer is included.

{It is well established that the atomic interstellar hydrogen is filling the galaxies and constitutes the building blocks of molecular clouds.} {To understand the formation and the evolution of molecular clouds, it is necessary to investigate the dynamics of turbulent and thermally bistable as well as barotropic flows.} {We perform high resolution 3-dimensional hydrodynamical simulations of 2-phase, isothermal and polytropic flows.} {We compare the density PDF and Mach number density relation in the various simulations and conclude that 2-phase flows behave rather differently than polytropic flows. We also extract the clumps and study their statistical properties such as mass spectrum, mass-size relation and internal velocity dispersion. In each case, it is found that the behaviour is well represented by a simple powerlaw. While the various exponents inferred are very similar for the 2-phase, isothermal and polytropic flows, we nevertheless find significant differences, such as for example the internal velocity dispersion which is smaller for 2-phase flows.} {The structures statistics are very similar to what has been inferred from observations, in particular the mass spectrum, the mass-size relation and the velocity dispersion-size relation are all powerlaws whose indices well agree with the observed values. Our results suggest that in spite of various statistics being similar for 2-phase and polytropic flows, they nevertheless present significant differences, stressing the necessity to consider the proper thermal structure of the interstellar atomic hydrogen for computing its dynamics as well as the formation of molecular clouds.}

We investigate the case of CII 158 micron observations for SPICA/SAFARI using a three-dimensional magnetohydrodynamical (MHD) simulation of the diffuse interstellar medium (ISM) and the Meudon PDR code. The MHD simulation consists of two converging flows of warm gas (10,000 K) within a cubic box 50 pc in length. The interplay of thermal instability, magnetic field and self-gravity leads to the formation of cold, dense clumps within a warm, turbulent interclump medium. We sample several clumps along a line of sight through the simulated cube and use them as input density profiles in the Meudon PDR code. This allows us to derive intensity predictions for the CII 158 micron line and provide time estimates for the mapping of a given sky area.

In the context of star and planet formation, understanding the formation of disks is of fundamental importance. Previous studies found that the magnetic field has a very strong impact on the collapse of a prestellar cloud, particularly in possibly suppressing the formation of a disk even for relatively modest values of the magnetic intensity. Since observations infer that cores have a substantial level of magnetization, this raises the question of how disks form. However, most studies have been restricted to the case in which the initial angle, $\alpha$, between the magnetic field and the rotation axis equals 0$^\circ$. We explore and analyse the influence of non aligned configurations on disk formation. We perform 3D ideal MHD, AMR numerical simulations for various values of $\mu$, the ratio of the mass-to-flux to the critical mass-to-flux, and various values of $\alpha$. We find that disks form more easily as $\alpha$ increases from 0 to 90$^\circ$. We propose that as the magnetized pseudo-disks become thicker with increasing $\alpha$, the magnetic braking efficiency is lowered. We also find that even small values of $\alpha$ ($\simeq$ 10-20$^\circ$) show significant differences with the alligned case. Within the framework of ideal MHD and for our choice of initial conditions, centrifugally supported disks cannot form for values of $\mu$ smaller than $\simeq$3, if the magnetic field and the rotation axis are perpendicular, and smaller than about $\simeq$5-10 when they are perfectly aligned.

Recently, Hennebelle and Chabrier (2008) derived an analytical theory for the mass spectrum of non self-gravitating clumps associated with overdensities in molecular clouds and for the initial mass function of gravitationally bound prestellar cores, as produced by the turbulent collapse of the cloud. In this companion paper, we examine the effects of the non-isothermality of the flow, of the turbulence forcing and of local fluctuation of the velocity dispersion, on the mass function. In particular, we investigate the influence of a polytropic equation of state and of the effective adiabatic exponent $\gamma$ and find that it has a drastic influence on the low mass part of the IMF. We also consider a barotropic equation of state (i.e. a piecewise polytropic eos) that mimics the thermal behaviour of the molecular gas and compare the prediction of our theory with the results of numerical simulations and with the observationally-derived IMF, for cloud parameters which satisfy Larson's type relations. We find that for clouds whose density is, at all scales, almost an order of magnitude larger than the density inferred for the CO clumps in the Galaxy, a good agreement is obtained between the theory and the observed IMF, suggesting that star formation preferentially occurs in high density environments. We derive an analytical expression for the IMF which generalizes the expression previously obtained for the isothermal case. This easy-to-implement analytical IMF should serve as a template to compare observational or numerical results with the theory.

Abridged: We study the properties of clumps formed in three-dimensional weakly magnetized magneto-hydrodynamic simulations of converging flows in the thermally bistable, warm neutral medium (WNM). We find that: (1) Similarly to the situation in the classical two-phase medium, cold, dense clumps form through dynamically-triggered thermal instability in the compressed layer between the convergent flows, and are often characterised by a sharp density jump at their boundaries though not always. (2) However, the clumps are bounded by phase-transition fronts rather than by contact discontinuities, and thus they grow in size and mass mainly by accretion of WNM material through their boundaries. (3) The clump boundaries generally consist of thin layers of thermally unstable gas, but these layers are often widened by the turbulence, and penetrate deep into the clumps. (4) The clumps are approximately in both ram and thermal pressure balance with their surroundings, a condition which causes their internal Mach numbers to be comparable to the bulk Mach number of the colliding WNM flows. (5) The clumps typically have mean temperatures 20 < T < 50 K, corresponding to the wide range of densities they contain (20 < n < 5000 pcc) under a nearly-isothermal equation of state. (6) The turbulent ram pressure fluctuations of the WNM induce density fluctuations that then serve as seeds for local gravitational collapse within the clumps. (7) The velocity and magnetic fields tend to be aligned with each other within the clumps, although both are significantly fluctuating, suggesting that the velocity tends to stretch and align the magnetic field with it. (8) The typical mean field strength in the clumps is a few times larger than that in the WNM. (9) The magnetic field strength has a mean value of B ~ 6 mu G ...

The radioactive isotope $^{60}$Fe ($T_{1/2} = 1.5 $ Myr) was present in the early solar system. It is unlikely that it was injected directly into the nascent solar system by a single, nearby supernova. It is proposed instead that it was inherited during the molecular cloud stage from several supernovae belonging to previous episodes of star formation. The expected abundance of $^{60}$Fe in star forming regions is estimated taking into account the stochasticity of the star-forming process, and it is showed that many molecular clouds are expected to contain $^{60}$Fe (and possibly $^{26}$Al [$T_{1/2} = 0.74 $ Myr]) at a level compatible with that of the nascent solar system. Therefore, no special explanation is needed to account for our solar system's formation.

This White Paper, submitted to the National Academy of Sciences' Astro2010 Decadal Review Committee, focuses on 2 central themes in the study of young brown dwarfs -- their formation mechanism and disk characteristics -- which are of direct relevance to fundamental questions of stellar and planetary origins and properties.

We review the theories and observations of high-mass star formation emphasizing the differences with those of low-mass star formation. We hereafter describe the progress expected to be achieved with Herschel, thanks notably to Key Programmes dedicated to the earliest phases of high-mass star formation.

{It has recently been proposed that giant molecular complexes form at the sites where streams of diffuse warm atomic gas collide at transonic velocities.} {We study the global statistics of molecular clouds formed by large scale colliding flows of warm neutral atomic interstellar gas under ideal MHD conditions. The flows deliver material as well as kinetic energy and trigger thermal instability leading eventually to gravitational collapse.} {We perform adaptive mesh refinement MHD simulations which, for the first time in this context, treat self-consistently cooling and self-gravity.} {The clouds formed in the simulations develop a highly inhomogeneous density and temperature structure, with cold dense filaments and clumps condensing from converging flows of warm atomic gas. In the clouds, the column density probability density distribution (PDF) peaks at $\sim 2 \times 10^{21} \psc$ and decays rapidly at higher values; the magnetic intensity correlates weakly with density from $n \sim 0.1$ to $10^4 \pcc$, and then varies roughly as $n^{1/2}$ for higher densities.} {The global statistical properties of such molecular clouds are reasonably consistent with observational determinations. Our numerical simulations suggest that molecular clouds formed by the moderately supersonic collision of warm atomic gas streams.}

We derive an analytical theory of the prestellar core initial mass function based on an extension of the Press-Schechter statistical formalism. With the same formalism, we also obtain the mass spectrum for the non self-gravitating clumps produced in supersonic flows. The mass spectrum of the self-gravitating cores reproduces very well the observed initial mass function and identifies the different mechanisms responsible for its behaviour. The theory predicts that the shape of the IMF results from two competing contributions, namely a power-law at large scales and an exponential cut-off (lognormal form) centered around the characteristic mass for gravitational collapse. The cut-off exists already in the case of pure thermal collapse, provided that the underlying density field has a lognormal distribution. Whereas pure thermal collapse produces a power-law tail steeper than the Salpeter value, dN/dlog M\propto M^{-x}, with x=1.35, this latter is recovered exactly for the (3D) value of the spectral index of the velocity power spectrum, n\simeq 3.8, found in observations and in numerical simulations of isothermal supersonic turbulence. Indeed, the theory predicts that x=(n+1)/(2n-4) for self-gravitating structures and x=2-n'/3 for non self-gravitating structures, where n' is the power spectrum index of log(rho). We show that, whereas supersonic turbulence promotes the formation of both massive stars and brown dwarfs, it has an overall negative impact on star formation, decreasing the star formation efficiency. This theory provides a novel theoretical foundation to understand the origin of the IMF and to infer its behaviour in different environments. It also provides a complementary approach and useful guidance to numerical simulations exploring star formation, while making testable predictions.

The development of parallel supercomputers allows today the detailed study of the collapse and the fragmentation of prestellar cores with increasingly accurate numerical simulations. Thanks to the advances in sub-millimeter observations, a wide range of observed initial conditions enable us to study the different modes of low-mass star formation. The challenge for the simulations is to reproduce the observational results. Two main numerical methods, namely AMR and SPH, are widely used to simulate the collapse and the fragmentation of prestellar cores. We compare thoroughly these two methods within their standard framework. We use the AMR code RAMSES and the SPH code DRAGON. Our physical model is as simple as possible and consists of an isothermal sphere rotating around the z-axis. We first study the conservation of angular momentum as a function of the resolution. Then, we explore a wide range of simulation parameters to study the fragmentation of prestellar cores. There seems to be a convergence between the two methods, provided resolution in each case is sufficient. Resolution criteria adapted to our physical cases, in terms of resolution per Jeans mass, for an accurate description of the formation of protostellar cores are deduced from the present study. This convergence is encouraging for future work in simulations of low-mass star formation, providing the aforementioned criteria are fulfilled. Higher resolution figures can be downloaded at http://www-dapnia.cea.fr/Projets/COAST/paper_amrvssph.pdf

(Abridged) The formation of molecular clouds (MCs) from the diffuse interstellar gas appears to be a necessary step for star formation, as young stars invariably occur within them. However, the mechanisms controlling the formation of MCs remain controversial. In this contribution, we focus on their formation in compressive flows driven by interstellar turbulence and large-scale gravitational instability. Turbulent compression driven by supernovae appears insufficient to explain the bulk of cloud and star formation. Rather, gravity must be important at all scales, driving the compressive flows that form both clouds and cores. Cooling and thermal instability allow the formation of dense gas out of moderate, transonic compressions in the warm diffuse gas, and drive turbulence into the dense clouds. MCs may be produced by an overshoot beyond the thermal-pressure equilibrium between the cold and warm phases of atomic gas, caused by some combination of the ram pressure of compression and the self-gravity of the compressed gas. In this case, properties of the clouds such as their mass, mass-to-magnetic flux ratio, and total kinetic and gravitational energies are in general time-variable quantities. MCs may never enter a quasi-equilibrium or virial equilibrium state but rather continuously collapse to stars.

Abridged. A large fraction of stars are found in binary systems. It is therefore important for our understanding of the star formation process, to investigate the fragmentation of dense molecular cores. We study the influence of the magnetic field, ideally coupled to the gas, on the fragmentation in multiple systems of collapsing cores. We present high resolution numerical simulations performed with the RAMSES MHD code starting with a uniform sphere in solid body rotation and a uniform magnetic field parallel to the rotation axis. We pay particular attention to the strength of the magnetic field and interpret the results using the analysis presented in a companion paper. The results depend much on the amplitude, $A$, of the perturbations seeded initially. For a low amplitude, $A=0.1$, we find that for values of the mass-to-flux over critical mass-to-flux ratio, $\mu$, as high as $\mu = 20$, the centrifugally supported disk which fragments in the hydrodynamical case, is stabilized and remains axisymmetric. Detailed investigations reveals that this is due to the rapid growth of the toroidal magnetic field induced by the differential motions within the disk. For values of $\mu$ smaller $\simeq 5$, corresponding to larger magnetic intensities, there is no centrifugally supported disk because of magnetic braking. When the amplitude of the perturbation is equal to $A=0.5$, each initial peak develops independently and the core fragments for a large range of $\mu$. Only for values of $\mu$ close to 1 is the magnetic field able to prevent the fragmentation. Since a large fraction of stars are binaries, the results of low magnetic intensities preventing the fragmentation in case of weak perturbations, is problematic. We discuss three possible mechanisms...

Abridged. It is important for the star formation process to understand the collapse of a prestellar dense core. We investigate the effect of the magnetic field during the first collapse up to the formation of the firstcore, focusing particularly on the magnetic braking and the launching of outflows. We perform 3D AMR high resolution numerical simulations of a magnetically supercritical collapsing dense core using the RAMSES MHD code and develop semi-analytical models that we compare with the numerical results. We study in detail the various profiles within the envelope of the collapsing core for various magnetic field strengths. Even modest values of magnetic field strength modify the collapse significantly. This is largely due to the amplification of the radial and toroidal components of the magnetic field by the differential motions within the collapsing core. For a weak magnetic intensity corresponding to an initial mass-to-flux over critical mass-to-flux ratio, $\mu$ equals to 20, a centrifugally supported disk forms. The strong differential rotation triggers the growth of a slowly expanding magnetic tower. For a higher magnetic field strengths corresponding to $\mu=2$, the collapse occurs primarily along the field lines, therefore delivering weaker angular momentum in the inner part whereas at the same time, strong magnetic braking occurs. As a consequence no centrifugally supported disk forms. An outflow is launched from the central thermally supported core. Detailed comparisons with existing analytical predictions indicate that it is magneto-centrifugally driven. For cores having a mass-to-flux over critical mass-to-flux radio $\mu < 5$, the magnetic field appears to have a significant impact.....

The development of parallel supercomputers allows today the detailed study of the collapse and the fragmentation of prestellar cores with increasingly accurate numerical simulations. Thanks to the advances in sub-millimeter observations, a wide range of observed initial conditions enable us to study the different modes of low-mass star formation. The challenge for the simulations is to reproduce the observational results. Two main numerical methods, namely AMR and SPH, are widely used to simulate the collapse and the fragmentation of prestellar cores. We compare here thoroughly these two methods with numerical resolution requirements deduced from previous studies. Our physical model is as simple as possible, and consists of an isothermal sphere rotating around the z-axis. We first study the conservation of angular momentum as a function of the resolution. Then, we explore a wide range of simulation parameters to study the fragmentation of prestellar cores. There seems to be a convergence between the two methods, provided resolution in each case is sufficient. Resolution criteria adapted to our physical cases, in terms of resolution per Jeans mass, for an accurate description of the formation of protostellar cores are deduced from the present study. This convergence is encouraging for future work in simulations of low-mass star formation, providing the aforementioned criteria are fulfilled.

Aims. We probe the effect of turbulent diffusion on the chemistry at the interface between a cold neutral medium (CNM) cloudlet and the warm neutral medium (WNM). Methods. We perform moving grid, multifluid, 1D, hydrodynamical simulations with chemistry including thermal and chemical diffusion. The diffusion coefficients are enhanced to account for turbulent diffusion. We post-process the steady-states of our simulations with a crude model of radiative transfer to compute line profiles. Results. Turbulent diffusion spreads out the transition region between the CNM and the WNM. We find that the CNM slightly expands and heats up: its CH and H$_2$ content decreases due to the lower density. The change of physical conditions and diffusive transport increase the H$^+$ content in the CNM which results in increased OH and H$_2$O. Diffusion transports some CO out of the CNM. It also brings H$_2$ into contact with the warm gas with enhanced production of CH$^+$, H$_3^+$, OH and H$_2$O at the interface. O lines are sensitive to the spread of the thermal profile in the intermediate region between the CNM and the WNM. Enhanced molecular content at the interface of the cloud broadens the molecular line profiles and helps exciting transitions of intermediate energy. The relative molecular yield are found higher for bigger clouds. Conclusions. Turbulent diffusion can be the source of additional molecular production and should be included in chemical models of the interstellar medium (ISM). It also is a good candidate for the interpretation of observational problems such as warm H$_2$, CH$^+$ formation and presence of H$_3^+$.

The interstellar atomic hydrogen is known to be a 2-phase medium in which turbul ence plays an important r\^ole. Here we present high resolution numerical simulations describing the gas from tens of parsec down to hundreds of AU. This high resolut ion allows to probe numerically, the small scale structures which naturally arises from the turbulence and the 2-phase physics.

Detailed modelling of individual protostellar condensations, is important to test the various theories. Here we present comparisons between strongly induced collapse models with one young class-0, IRAS4A in the Perseus cloud and one prestellar cloud observed in the Coalsack molecular cloud.

The neutral atomic interstellar medium (HI) is well known to be strongly magnetized. Since HI is a 2-phase medium the questions arise of what is the effect of the magnetic field on a 2-phase medium, how magnetic field affects thermal instability ? Here we summarize analytical and numerical results which have been obtained previously to investigate this question.

{It is necessary to understand the dynamics of the atomic gas to use complex modeling and to carry out detailed comparisons between theoretical models and observations.}{In a companion paper, we present high resolution bidimensional numerical simulations of the interstellar atomic hydrogen. Here, we further characterize these simulations and we compare our results with various observations.}{We give statistics of the column density and velocity along the line of sight and show that they compare favorably with observations of high-latitude lines of sight. We compute synthetic HI spectra and qualitatively discuss the information that could be inferred if these spectra were observed. Finally, we extract CNM clouds and study their physical properties finding strong similarities with real clouds. In particular, we find that the clouds follow Larson-type relations, i.e. $M \propto L^\gamma$, where $\gamma \simeq 1.7$ (we speculate that in 3D, $\gamma \simeq 2.5$) and $\sqrt{<\delta v^2>} \propto L^{0.4}$. We also find that the distribution, ${\cal N}(N)$, of the column density, $N$, of the CNM structures formed in the simulation follows ${\cal N}(N) \propto N^{-1.2}$ which is marginally compatible with the observational result obtained by Heiles & Troland (2005). From the mass-size relation and the mass spectrum, we derive an exponent for the column density distribution close to the value obtained in the numerical simulation.} {We conclude that the simulations reproduce various observational features reasonably well....}

{We study in some details the statistical properties of the turbulent 2-phase interstellar atomic gas.{We present high resolution bidimensional numerical simulations of the interstellar atomic hydrogen which describe it over 3 to 4 orders of magnitude in spatial scales.}{The simulations produce naturally small scale structures having either large or small column density. It is tempting to propose that the former are connected to the tiny small scale structures observed in the ISM. We compute the mass spectrum of CNM structures and find that ${\cal N}(M) dM \propto M ^{-1.7} dM$, which is remarkably similar to the mass spectrum inferred for the CO clumps. We propose a theoretical explanation based on a formalism inspired from the Press & Schecter (1974) approach and used the fact that the turbulence within WNM is subsonic. This theory predicts ${\cal N}(M) \propto M ^{-5/3}$ in 2D and ${\cal N}(M) \propto M ^{-16/9}$ in 3D. We compute the velocity and the density power-spectra and conclude that, although the latter is rather flat, as observed in supersonic isothermal simulations, the former follows the Kolmogorov prediction and is dominated by its solenoidal component. This is due to the bistable nature of the flow which produces large density fluctuations even when the rms Mach number (of WNM) is not large. We also find that, whereas the energy at large scales is mainly in the WNM, at smaller scales, it is dominated by the kinetic energy of the CNM fragments.}

The earliest phases of massive star formation in clusters are still poorly understood. Here, we test the hypothesis for high-mass star formation proposed in our earlier paper (Peretto et al. 2006). In order to confirm the physical validity of this hypothesis, we carried out IRAM Plateau de Bure interferometer observations of NGC 2264-C and performed SPH numerical simulations of the collapse of a Jeans-unstable, prolate dense clump. Our Plateau de Bure observations reveal the presence of a new compact source (C-MM13) located only \~ 10000 AU away, but separated by ~ 1.1 km/s in (projected) velocity, from the most massive Class 0 object (C-MM3) lying at the very center of NGC 2264-C. Detailed comparison with our numerical SPH simulations supports the view that NGC 2264-C is an elongated cluster-forming clump in the process of collapsing and fragmenting along its long axis, leading to a strong dynamical interaction and possible protostar merger in the central region of the clump. The present study also sets several quantitative constraints on the initial conditions of large-scale collapse in NGC 2264-C. Our hydrodynamic simulations indicate that the observed velocity pattern characterizes an early phase of protocluster collapse which survives for an only short period of time (i.e., < 10^5 yr). To provide a good match to the observations the simulations require an initial ratio of turbulent to gravitational energy of only ~ 5 %, which strongly suggests that the NGC 2264-C clump is structured primarily by gravity rather than turbulence. The required "cold'' initial conditions may result from rapid compression by an external trigger.

In this paper, we present a new method to perform numerical simulations of astrophysical MHD flows using the Adaptive Mesh Refinement framework and Constrained Transport. The algorithm is based on a previous work in which the MUSCL--Hancock scheme was used to evolve the induction equation. In this paper, we detail the extension of this scheme to the full MHD equations and discuss its properties. Through a series of test problems, we illustrate the performances of this new code using two different MHD Riemann solvers (Lax-Friedrich and Roe) and the need of the Adaptive Mesh Refinement capabilities in some cases. Finally, we show its versatility by applying it to two completely different astrophysical situations well studied in the past years: the growth of the magnetorotational instability in the shearing box and the collapse of magnetized cloud cores. We have implemented this new Godunov scheme to solve the ideal MHD equations in the AMR code RAMSES. It results in a powerful tool that can be applied to a great variety of astrophysical problems, ranging from galaxies formation in the early universe to high resolution studies of molecular cloud collapse in our galaxy.

Recent high resolution numerical simulations have suggested that the interstellar atomic hydrogen clouds have a complex two-phase structure. Since molecular clouds form through the contraction of HI gas, the question arises as to whether this structure is maintained in the molecular phase or not. Here we investigate whether the warm neutral atomic hydrogen (WNM) can exist in molecular clouds. We calculate how far a piece of WNM which is not heated by the UV photons could penetrate into the cloud, and find that in the absence of any heating it is unlikely that large fraction of WNM survives inside high pressure molecular clouds. We then consider two possible heating mechanisms, namely dissipation of turbulent energy and dissipation of MHD waves propagating in the WNM inside the cloud. We find that the second one is sufficient to allow the existence of WNM inside a molecular cloud of size $\simeq$ 1 pc having pressure equal to $\simeq 10 \times P_{\rm ISM}$. This result suggests the possibility that channels of magnetised WMN may provide efficient energy injection for sustaining internal turbulence which otherwise decays in a crossing time.

The onset of gravitational collapse in cluster-forming clouds is still poorly known. Our goal is to use the Class 0 protostar IRAS 4A, which is undergoing collapse in the active molecular cloud NGC 1333, to set constraints on this process. In particular we want to measure the mass infall rate and investigate whether the collapse could have been triggered by a strong external perturbation. We analyze existing continuum observations to derive the density structure of the envelope, and use our new molecular line observations done with the IRAM 30m telescope to probe its velocity structure. We perform a detailed comparison of this set of data with a numerical model of collapse triggered by a fast external compression. Both the density and velocity structures of the envelope can be well fitted by this model of collapse induced by a fast external compression for a time elapsed since point mass formation of 1-2 x 10**4 yr. We deduce a large mass infall rate of 0.7-2 x 10**-4 Msun/yr. The momentum required for the perturbation to produce this large mass infall rate is of the same order as the momenta measured for the NGC 1333 numerous outflows. Our analysis shows also that the turbulence is highly non uniform in the envelope, dropping from supersonic to subsonic values toward the center. The inner subsonic turbulence is most likely a relic of the conditions prevailing in the dense core before the onset of collapse. The vigorous collapse undergone by IRAS 4A was triggered by a fast external compression, probably related to the expansion of a nearby cavity, which could have triggered the collapse of the nearby Class 0 protostar IRAS 4B simultaneously. This cavity could have been generated by an outflow but we have not found a good protostellar candidate yet.

We review recent theoretical progress aimed at understanding the formation and the early stages of evolution of giant planets, low-mass stars and brown dwarfs. Calculations coupling giant planet formation, within a modern version of the core accretion model, and subsequent evolution yield consistent determinations of the planet structure and evolution. Because of the uncertainties in the initial conditions, however, it is not possible to say whether young planets are faint or bright compared with low-mass young brown dwarfs. We review the effects of irradiation and evaporation on the evolution of short period planets and argue that substantial mass loss may have occurred for these objects. Concerning star formation, geometrical effects in protostar core collapse are examined by comparing 1D and 3D calculations. Spherical collapse is shown to overestimate the core inner density and temperature and thus to yield incorrect initial conditions for PMS or young brown dwarf evolution. Accretion is also shown to occur over a very limited fraction of the protostar surface. Accretion affects the evolution of young brown dwarfs and yields more compact structures for a given mass and age, thus fainter luminosities. This can lead to severe misinterpretations of the mass and/or age of young accreting objects from their location in the HR diagram. We argue that newborn stars and brown dwarfs should appear rapidly over an extended area in the HR diagram, depending on their accretion history, rather than on a well defined birth line. Finally, we suggest that the distinction between planets and brown dwarfs be based on an observational diagnostic, reflecting the different formation mechanisms between these two distinct populations, rather than on an arbitrary, confusing definition.

Lada et al. recently presented a detailed near-infrared extinction map of Globule G2 in the Coalsack molecular cloud complex, showing that this starless core has a well-defined central extinction minimum. We propose a model for G2 in which a rapid increase in external pressure is driving an approximately symmetric compression wave into the core. The rapid increase in external pressure could arise because the core has recently been assimilated by the Coalsack cloud complex, or because the Coalsack has recently been created by two large-scale converging flows. The resulting compression wave has not yet converged on the centre of the core, so there is a central rarefaction. The compression wave has increased the density in the swept-up gas by about a factor of ten, and accelerated it inwards to speeds of order $0.4 {\rm km} {\rm s}^{-1}$. It is shown that even small levels of initial turbulence destroy the ring seen in projection almost completely. In the scenario of strong external compression that we are proposing this implies that the initial turbulent energy in this globule is such that $E_{{\rm turb}} / E_{{\rm grav}} \le 2 %$. Protostar formation should occur in about $40,000 {\rm years}$.

The effect of Alfv\'en waves on the thermal instability of the Interstellar Medium (ISM) is investigated both analytically and numerically. A stability analysis of a finite amplitude circularly polarized Alfv\'en wave propagating parallel to an ambient magnetic field in a thermally unstable gas at thermal equilibrium is performed, leading to a dispersion relation which depends on 3 parameters, namely the square ratio of the sonic and Alfv\'en velocities ($\beta$), the wave amplitude and the ratio between the wave temporal period and the cooling time. Depending on the values of these 3 parameters, the Alfv\'en waves can stabilize the large-scale perturbations, destabilize those whose wavelength is a few times the Alfv\'en wavelength $\lambda_{AW}$ or leave the growth rate of the short scales unchanged. In order to investigate the non-linear regime, two different numerical experiments are performed in a slab geometry. (TBC)

We present a numerical and analytical study of the thermal fragmentation of a turbulent flow of interstellar hydrogen. We first present the different dynamical processes and the large range of spatial (and temporal) scales that need to be adequately represented in numerical simulations. Next, we present bidimensionnal simulations of turbulent converging flows which induce the dynamical condensation of the warm neutral phase into the cold phase. We then analyse the cold structures and the fraction of unstable gas in each simulation, paying particular attention to the influence of the degree of turbulence. When the flow is very turbulent a large fraction of the gas remains in the thermally unstable domain. This unstable gas forms a filamentary network. We show that the fraction of thermally unstable gas is strongly correlated with the level of turbulence of the flow. We then develop a semi-analytical model to explain the origin of this unstable gas. This simple model is able to reproduce quantitatively the fraction of unstable gas observed in the simulations and its correlation with turbulence. Finally, we stress the fact that even when the flow is very turbulent and in spite of the fact that a large fraction of the gas is maintained dynamically in the thermally unstable domain, the classical picture of a 2-phase medium with stiff thermal fronts and local pressure equilibrium turns out to be still relevant in the vicinity of the cold structures.

Improving our understanding of the earliest stages of star formation is crucial to gain insight into the origin of stellar masses, multiple systems, and protoplanetary disks. We discuss recent advances made in this area thanks to detailed mapping observations at infrared and (sub)millimeter wavelengths. Although ambipolar diffusion appears to be too slow to play a direct role in the formation of dense cores, there is nevertheless good evidence that the gravitational collapse of isolated protostellar cores is strongly magnetically controlled. We also argue that the beginning of protostellar collapse is much more violent in cluster-forming clouds than in regions of distributed star formation.

We investigate numerically and semi-analytically the collapse of low-mass, rotating prestellar cores. Initially, the cores are in approximate equilibrium with low rotation (the initial ratio of thermal to gravitational energy is $\alpha_0 \simeq 0.5$, and the initial ratio of rotational to gravitational energy is $\beta_0 = 0.02 {\rm to} 0.05$). They are then subjected to a steady increase in external pressure. Fragmentation is promoted -- in the sense that more protostars are formed -- both by more rapid compression, and by higher rotation (larger $\beta_0$). In general, the large-scale collapse is non-homologous, and follows the pattern described in Paper I for non-rotating clouds, viz. a compression wave is driven into the cloud, thereby increasing the density and the inflow velocity. The effects of rotation become important at the centre, where the material with low angular momentum forms a central primary protostar (CPP), whilst the material with higher angular momentum forms an accretion disc around the CPP. More rapid compression drives a stronger compression wave and delivers material more rapidly into the outer parts of the disc.

We propose an extension of the semi-analytical solutions derived by Lin et al. (1965) describing the two-dimensional homologous collapse of a self-gravitating rotating cloud having uniform density and spheroidal shape, which includes magnetic field (with important restrictions) and thermal pressure. The evolution of the cloud is reduced to three time dependent ordinary equations allowing to conduct a quick and preliminary investigation of the cloud dynamics during the precollapse phase, for a wide range of parameters. We apply our model to the collapse of a rotating and magnetized oblate and prolate isothermal core. Hydrodynamical numerical simulations are performed and comparison with the semi-analytical solutions is discussed. Under the assumption that all cores are similar, an apparent cloud axis ratio distribution is calculated from the sequence of successive evolutionary states for each of a large set of initial conditions. The comparison with the observational distribution of the starless dense cores belonging to the catalog of Jijina et al. (1999) shows a good agreement for the rotating and initially prolate cores (aspect ratio $\simeq 0.5$) permeated by an helical magnetic field ($\simeq 17-20 \mu$G for a density of $\simeq 10^4$ cm$^{-3}$).

The ISOGAL project is an infrared survey of specific regions sampling the Galactic Plane selected to provide information on Galactic structure,stellar populations,stellar mass-loss and the recent star formation history of the inner disk and Bulge of the Galaxy. ISOGAL combines 7 and 15 micron ISOCAM observations - with a resolution of 6'' at worst - with DENIS IJKs data to determine the nature of the sources and theinterstellar extinction. We have observed about 16 square degrees with a sensitivity approaching 10-20mJy, detecting ~10^5 sources,mostly AGB stars,red giants and young stars. The main features of the ISOGAL survey and the observations are summarized in this paper,together with a brief discussion of data processing and quality. The primary ISOGAL products are described briefly (a full description is given in Schuller et al. 2003, astro-ph/0304309): viz. the images and theISOGAL-DENIS five-wavelength point source catalogue. The main scientific results already derived or in progress are summarized. These include astrometrically calibrated 7 and 15um images,determining structures of resolved sources; identification and properties of interstellar dark clouds; quantification of the infrared extinction law and source dereddening; analysis of red giant and (especially) AGB stellar populations in the central Bulge,determining luminosity,presence of circumstellar dust and mass--loss rate,and source classification,supplemented in some cases by ISO/CVF spectroscopy; detection of young stellar objects of diverse types,especially in the inner Bulge with information about the present and recent star formation rate; identification of foreground sources with mid-IR excess. These results are the subject of about 25 refereed papers published or in preparation.

Improving our understanding of the initial conditions and earliest stages of protostellar collapse is crucial to gain insight into the origin of stellar masses, multiple systems, and protoplanetary disks. Observationally, there are two complementary approaches to this problem: (1) studying the structure and kinematics of prestellar cores observed prior to protostar formation, and (2) studying the structure of young (e.g. Class 0) accreting protostars observed soon after point mass formation. We discuss recent advances made in this area thanks to (sub)millimeter mapping observations with large single-dish telescopes and interferometers. In particular, we argue that the beginning of protostellar collapse is much more violent in cluster-forming clouds than in regions of distributed star formation. Major breakthroughs are expected in this field from future large submillimeter instruments such as Herschel and ALMA.

Ordinary differential equations describing the self-similar collapse of a rotating, magnetized, self-gravitating and isothermal filament are derived. Explicit homologous solutions are studied with special emphasis on the bifurcation that occurs at the magnetosonic critical point. It is shown that there is a critical value for the toroidal magnetic field slope at the origin above which no bifurcation occurs, the solution remains homologous, and below which the density and the poloidal magnetic field tend to zero at large radius (envelope) whereas the toroidal magnetic field and azimuthal velocity relax towards a constant value. A series of spatial profiles of density, velocity and magnetic field, potentially useful for comparison with numerical or observational studies, is obtained numerically and discussed.

We present numerical simulations of the evolution of low-mass, isothermal, molecular cores which are subjected to an increase in external pressure $P\xt$. If $P\xt$ increases very slowly, the core approaches instability quite quasistatically. However, for larger (but still quite modest) $dP\xt/dt$ a compression wave is driven into the core, thereby triggering collapse from the outside in. If collapse of a core is induced by increasing $P\xt$, this has a number of interesting consequences. (i) The density profile is approximately flat in the centre during the prestellar phase (i.e. before the compression wave converges on the centre creating the central protostar). (ii) During the prestellar phase there are (subsonic) inward velocities in the outer layers of the core, whilst the inner parts are still approximately at rest. (iii) There is an initial short phase of rapid accretion (notionally the Class 0 phase), followed by a longer phase of slower accretion (the Class I phase). All these features accord well with observation, but are at variance with the predictions of the standard theory of star formation based on the inside-out collapse of a singular isothermal sphere. We note that the setting up of a coherent inward velocity field appears to be a generic feature of compression waves; and we speculate that interactions and interference between such velocity fields may play a crucial r\^ole in initiating the fragmentation of cores and the genesis of multiple star systems.

This paper presents follow-up observations of the mid-Infrared dark clouds selected from the ISOGAL inner Galaxy sample. On-the-fly maps of 13CO, C18O and the 1.2 mm continuum emission have been obtained a the IRAM 30-m telescope. A preliminary analysis of the physico-chemical conditions at works in the clouds is given and molecular depletion onto cold grains is discussed.

We find exact and explicit solutions of the axisymmetric MHD equations of a self-gravitating polytropic gas. These solutions are able to describe a flat (uniform density) subsonic internal core, contracting homologously, of a collapsing cloud with aspherical pressure, aspherical gravitationnal potential, magnetic field and rotation. Two of the solutions describing the collapse of a rotating and magnetized cloud present outflow at the pole or at the equator.