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) M⊙ 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(-385+770) M⊙ contained within a source diameter of 0.05 pc, 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 Vinf = 0.7( ± 0.2) km s-1, and a total mass infall rate Ṁinf ≃ 2.5(±1.0) × 10-3 M⊙ yr-1 towards the central pc-size region of SDC335. This infall rate brings 750(±300) M⊙ of gas to the centre of the cloud per free-fall time (tff = 3 × 105 yr). This is enough to double the mass already present in the central pc-size region in 3.5-1.0+2.2 × tff. 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.
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 three-dimensional analysis we describe the relative orientation of the magnetic field 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 nT in the highest magnetization case, and (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.
Context.
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.
Aims.
We ...study the impact of grain–grain coagulation on the evolution of magnetohydrodynamic (MHD) resistivities during the collapse of a prestellar core.
Methods.
We use a 1D model to follow the evolution of the dust size distribution, out-of-equilibrium ionisation 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.
Results.
Starting with a Mathis-Rumpl-Nordsieck (MRN) size distribution (Mathis et al. 1977, ApJ, 217, 425), 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 have a large effect on the evolution of the Hall and ambipolar diffusion MHD resistivities, which still drop significantly 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
H
~ 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
μ
m would have been removed in the parent cloud before the collapse.
Conclusions.
When grain drift velocities induced by ambipolar diffusion are included, dust coagulation happening during the collapse of a prestellar core starting from an initial MRN dust size distribution appears to be efficient enough to increase the MHD resistivities to the values necessary to strongly modify the magnetically regulated formation of a planet-forming disk. A consistent treatment of the competition between fragmentation and coagulation is, however, necessary before reaching firm conclusions.
Context: To understand the star formation process, it is important to study the collapse of a prestellar dense core. Aims: We investigate the effect of the magnetic field during the first collapse up ...to the formation of the first core, focusing particularly on the magnetic braking and the launching of outflows. Methods: 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. 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 20 a centrifugally supported disk forms. The strong differential rotation triggers the growth of a slowly expanding magnetic tower. For higher magnetic field strengths corresponding to mu=2, the collapse occurs primarily along the field lines, therefore delivering weaker angular momentum into 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. Conclusions: 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 collapsing envelope is denser and flatter than in the hydrodynamical case and no centrifugally supported disk forms. For values mu < 20, the magnetic field drastically modifies the disk evolution. In a companion paper, the influence of the magnetic field on the dense core fragmentation is studied.
Context.
Understanding the initial properties of star forming material and how they affect the star formation process is a key question. The infalling gas must redistribute most of its initial ...angular momentum inherited from prestellar cores before reaching the central stellar embryo. 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.
Aims.
In the framework of the IRAM CALYPSO survey, we obtained observations of the dense gas kinematics that we used to quantify the amount and distribution of specific angular momentum at all scales in collapsing-rotating Class 0 protostellar envelopes.
Methods.
We used the high dynamic range C
18
O (2−1) and N
2
H
+
(1−0) datasets to produce centroid velocity maps and probe the rotational motions in the sample of 12 envelopes from scales ~50 to ~5000 au.
Results.
We identify differential rotation motions at scales ≲1600 au in 11 out of the 12 protostellar envelopes of our sample by measuring the velocity gradient along the equatorial axis, which we fit with a power-law model v ∝
r
α
. This suggests that coherent motions dominate the kinematics in the inner protostellar envelopes. The radial distributions of specific angular momentum in the CALYPSO sample suggest the following two distinct regimes within protostellar envelopes: the specific angular momentum decreases as
j
∝
r
1.6±0.2
down to ~1600 au and then tends to become relatively constant around ~6 × 10
−4
km s
−1
pc down to ~50 au.
Conclusions.
The values of specific angular momentum measured in the inner Class 0 envelopes suggest that material directly involved in the star formation process (<1600 au) has a specific angular momentum 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 these gradients may not be directly related to rotational motions of the envelopes. We conclude that the specific angular momentum observed at these scales could find its origin in other mechanisms, such as core-forming motions (infall, turbulence), or trace an imprint of the initial conditions for the formation of protostellar cores.
The interstellar medium (ISM) is a complex nonlinear system governed by the interplay between 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 of the emerging structures. The goal of this paper is to build such a description, with the purpose of permitting an efficient comparison that is independent of any specific prior or model. We started from the wavelet scattering transform (WST), a low-variance statistical description of non-Gaussian processes, which was developed in data science and encodes long-range interactions through a hierarchical multiscale approach based on the wavelet transform. We performed a reduction of the WST through a fit of its angular dependencies. This allowed us to gather most of the information it contains into a few components whose physical meanings are identified and describe for instance 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 nonlinear phenomena, in particular interstellar magnetohydrodynamical (MHD) turbulence, independently of any specific priors. 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 applied it successfully to the following three kinds of processes that are a priori very different: fractional Brownian motions, MHD simulations, and Herschel observations of the dust thermal continuum in a molecular cloud. With fewer than 100 RWST coefficients when probing six scales and eight angles on 256 by 256 maps, we were able to perform quantitative comparisons, infer relevant physical properties, and produce realistic synthetic fields.
The angular resolution of a telescope is the primary observational parameter, along with the detector sensitivity in defining the quality of the observed images and of the subsequent scientific ...exploitation of the data. During the last decade in star formation research, many studies have targeted low- and high-mass star formation regions located at different distances, with different telescopes having specific angular resolution capabilities. However, no dedicated studies of the spatial resolution effects on the derived sizes and masses of the sources extracted from the observed images have been published. We present a systematic investigation of the angular resolution effects, with special attention being paid to the derived masses of sources as well as the shape of the resulting source mass functions (SMFs) and to their comparison with the initial stellar mass function. For our study, we chose two star-forming regions observed with
Herschel
, NGC 6334 and Aquila distant of 1750 and 460 pc respectively, and three (magneto)-hydrodynamical simulations, virtually positioned at the same distances as the observed regions. We built surface density maps with different angular resolutions by convolving the surface density images of the five regions to a set of four resolutions differing by a factor of two (9, 18, 36, and 72′′), which allowed us to cover spatial resolutions from 0.6 down to 0.02 pc. Then we detected and measured sources in each of the images at each resolution using getsf and we analysed the derived masses and sizes of the extracted sources. We find that the number of sources does not converge from 0.6 to ≳0.05 pc. It increases by about two when the angular resolution increases with a similar factor, which confirms that these large sources are cluster-forming clumps. Below 0.05 pc, the number of source still increases by about 1.3 when the angular resolution increases by two, suggesting that we are close to, but not yet at, convergence. In this regime of physical scales, we find that the measured sizes and masses of sources linearly depend on the angular resolution with no sign of convergence to a resolution-independent value, implying that these sources cannot be assimilated to isolated prestellar cores. The corresponding SMF peak also shifts with angular resolution, while the slope of the high-mass tail of the SMFs remains almost invariant. We propose that these angular resolution effects could be caused by the underestimated background of the unresolved sources observed against the sloping, hill-like backgrounds of the molecular clouds. If prestellar cores physically distinct from their background exist in cluster-forming molecular clouds, we conclude that their mass must be lower than reported so far in the literature. We discuss various implications for the studies of star formation: the problem of determining the mass reservoirs involved in the star-formation process; the inapplicability of the Gaussian beam deconvolution to infer source sizes; and the impossibility to determine the efficiency of the mass conversion from the cores to the stars. Our approach constitutes a simple convergence test to determine whether an observation is affected by angular resolution.
We present a self-consistent hydrodynamical simulation of a Milky Way-like galaxy at a resolution of 0.05 pc. The model includes star formation and a new implementation of stellar feedback through ...photoionization, radiative pressure and supernovae. The simulation resolves the structure of the interstellar medium at sub-parsec resolution for a few cloud lifetimes and at 0.05 pc for about a cloud-crossing time. The turbulence cascade and gravitation from 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 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 organization of structures from the galactic scale down to a few pc.
Context. 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. Aims. In particular, ...the influence of the magnetic field on the collapse of massive clumps is relatively unexplored, it is therefore of great relevance in the context of the formation of massive stars to investigate its impact. Methods. We perform high resolution, MHD simulations of the collapse of one hundred solar masses, turbulent and magnetized clouds, with 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. Results. 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 1−3 km s-1, although the highest velocities can be as high as 20−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. Conclusions. 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 and therefore could possibly lead to an overestimation of the number of fragments.
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 bidimensional 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 quantitatively reproduce 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.