[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.

CO isotopologue transitions are routinely observed in molecular clouds to probe the column density of the gas, the elemental ratios of carbon and oxygen, and to trace the kinematics of the environment. We aim at estimating the abundances, excitation temperatures, velocity field and velocity dispersions of the three main CO isotopologues towards a subset of the Orion B molecular cloud. We use the Cramer Rao Bound (CRB) technique to analyze and estimate the precision of the physical parameters in the framework of local-thermodynamic-equilibrium excitation and radiative transfer with an additive white Gaussian noise. We propose a maximum likelihood estimator to infer the physical conditions from the 1-0 and 2-1 transitions of CO isotopologues. Simulations show that this estimator is unbiased and efficient for a common range of excitation temperatures and column densities (Tex > 6 K, N > 1e14 - 1e15 cm-2). Contrary to the general assumptions, the different CO isotopologues have distinct excitation temperatures, and the line intensity ratios between different isotopologues do not accurately reflect the column density ratios. We find mean fractional abundances that are consistent with previous determinations towards other molecular clouds. However, significant local deviations are inferred, not only in regions exposed to UV radiation field but also in shielded regions. These deviations result from the competition between selective photodissociation, chemical fractionation, and depletion on grain surfaces. We observe that the velocity dispersion of the C18O emission is 10% smaller than that of 13CO. The substantial gain resulting from the simultaneous analysis of two different rotational transitions of the same species is rigorously quantified. The CRB technique is a promising avenue for analyzing the estimation of physical parameters from the fit of spectral lines.

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]

Complex organic molecules (COMs) have been detected in a few Class 0 protostars but their origin is not well understood. Going beyond studies of individual objects, we want to investigate the origin of COMs in young protostars on a statistical basis. We use the CALYPSO survey performed with the IRAM PdBI to search for COMs at high angular resolution in a sample of 26 solar-type protostars, including 22 Class 0 and four Class I objects. Methanol is detected in 12 sources and tentatively in one source, which represents half of the sample. Eight sources (30%) have detections of at least three COMs. We find a strong chemical differentiation in multiple systems with five systems having one component with at least three COMs detected but the other component devoid of COM emission. The internal luminosity is found to be the source parameter impacting the most the COM chemical composition of the sources, while there is no obvious correlation between the detection of COM emission and that of a disk-like structure. A canonical hot-corino origin may explain the COM emission in four sources, an accretion-shock origin in two or possibly three sources, and an outflow origin in three sources. The CALYPSO sources with COM detections can be classified into three groups on the basis of the abundances of oxygen-bearing molecules, cyanides, and CHO-bearing molecules. These chemical groups correlate neither with the COM origin scenarii, nor with the evolutionary status of the sources if we take the ratio of envelope mass to internal luminosity as an evolutionary tracer. We find strong correlations between molecules that are a priori not related chemically (for instance methanol and methyl cyanide), implying that the existence of a correlation does not imply a chemical link. [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.

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.

Although, from a theoretical point of view, magnetic fields are believed to have a significant role during the early stages of star formation, especially during the main accretion phase, the magnetic field strength and topology is poorly constrained in the youngest accreting Class 0 protostars that lead to the formation of solar-type stars. We carried out observations of the polarized dust continuum emission with the SMA interferometer at 0.87mm, in order to probe the structure of the magnetic field in a sample of 12 low-mass Class 0 envelopes, including both single protostars and multiple systems, in nearby clouds. Our SMA observations probe the envelope emission at scales of 600-5000 au with a spatial resolution ranging from 600 to 1500 au depending on the source distance. We report the detection of linearly polarized dust continuum emission in all of our targets, with average polarization fractions ranging from 2% to 10% in these protostellar envelopes. The polarization fraction decreases with the continuum flux density, which translates into a decrease with the H2 column density within an individual envelope. Our analysis show that the envelope-scale magnetic field is preferentially observed either aligned or perpendicular to the outflow direction. Interestingly, our results suggest for the first time a relation between the orientation of the magnetic field and the rotational energy of envelopes, with a larger occurrence of misalignment in sources where strong rotational motions are detected at hundreds to thousands of au scales. We also show that the best agreement between the magnetic field and outflow orientation is found in sources showing no small-scale multiplicity and no large disks at about 100 au scales.

Understanding the processes that determine the stellar Initial Mass Function (IMF) is a critical unsolved problem, with profound implications for many areas of astrophysics. In molecular clouds, stars are formed in cores, gas condensations which are sufficiently dense that gravitational collapse converts a large fraction of their mass into a star or small clutch of stars. In nearby star-formation regions, the core mass function (CMF) is strikingly similar to the IMF, suggesting that the shape of the IMF may simply be inherited from the CMF. Here we present 1.3 mm observations, obtained with ALMA, the world's largest interferometer, of the active star-formation region W43-MM1, which may be more representative of the Galactic-disk regions where most stars form. The unprecedented resolution of these observations reveals, for the first time, a statistically robust CMF at high masses, with a slope that is markedly shallower than the IMF. This seriously challenges our understanding of the origin of the IMF.

The role played by gravity in the transfer of interstellar matter from molecular cloud scales to protostellar scales is still highly debated. Only detailed studies on the kinematics of large samples of star-forming clouds will settle the issue. We present new IRAM 30m observations of a sample of 27 infrared dark clouds covering a large range of sizes, masses, and aspect ratios. Preliminary results suggest that gravity is regulating the dynamical evolution of these clouds on a couple of parsec scales.