No Arabic abstract
The fraction of star formation that results in bound star clusters is influenced by the density spectrum in which stars are formed and by the response of the stellar structure to gas expulsion. We analyse hydrodynamical simulations of turbulent fragmentation in star-forming regions to assess the dynamical properties of the resulting population of stars and (sub)clusters. Stellar subclusters are identified using a minimum spanning tree algorithm. When considering only the gravitational potential of the stars and ignoring the gas, we find that the identified subclusters are close to virial equilibrium (the typical virial ratio Q_vir~0.59, where virial equilibrium would be Q_vir~0.5). This virial state is a consequence of the low gas fractions within the subclusters, caused by the accretion of gas onto the stars and the accretion-induced shrinkage of the subclusters. Because the subclusters are gas-poor, up to a length scale of 0.1-0.2 pc at the end of the simulation, they are only weakly affected by gas expulsion. The fraction of subclusters that reaches the high density required to evolve to a gas-poor state increases with the density of the star-forming region. We extend this argument to star cluster scales, and suggest that the absence of gas indicates that the early disruption of star clusters due to gas expulsion (infant mortality) plays a smaller role than anticipated, and is potentially restricted to star-forming regions with low ambient gas densities. We propose that in dense star-forming regions, the tidal shocking of young star clusters by the surrounding gas clouds could be responsible for the early disruption. This `cruel cradle effect would work in addition to disruption by gas expulsion. We suggest possible methods to quantify the relative contributions of both mechanisms.
We model the dynamical evolution of star forming regions with a wide range of initial properties. We follow the evolution of the regions substructure using the Q-parameter, we search for dynamical mass segregation using the Lambda_MSR technique, and we also quantify the evolution of local density around stars as a function of mass using the Sigma_LDR method. The amount of dynamical mass segregation measured by Lambda_MSR is generally only significant for subvirial and virialised, substructured regions - which usually evolve to form bound clusters. The Sigma_LDR method shows that massive stars attain higher local densities than the median value in all regions, even those that are supervirial and evolve to form (unbound) associations. We also introduce the Q-Sigma_LDR plot, which describes the evolution of spatial structure as a function of mass-weighted local density in a star forming region. Initially dense (>1000 stars pc^{-2}), bound regions always have Q >1, Sigma_LDR > 2 after 5Myr, whereas dense unbound regions always have Q < 1, Sigma_LDR > 2 after 5Myr. Less dense regions (<100 stars pc^{-2}) do not usually exhibit Sigma_LDR > 2 values, and if relatively high local density around massive stars arises purely from dynamics, then the Q-Sigma_LDR plot can be used to estimate the initial density of a star forming region.
We present analyses of the spatial distributions of stars in the young (1 - 3 Myr) star-forming regions IC348 and NGC1333 in the Perseus Giant Molecular Cloud. We quantify the spatial structure using the $mathcal{Q}$-parameter and find that both IC348 and NGC1333 are smooth and centrally concentrated with $mathcal{Q}$-parameters of 0.98 and 0.89 respectively. Neither region exhibits mass segregation ($Lambda_{rm MSR} = 1.1^{+0.2}_{-0.3}$ for IC348 and $Lambda_{rm MSR} = 1.2^{+0.4}_{-0.3}$ for NGC1333, where $Lambda_{rm MSR} sim 1$ corresponds to no mass segregation), nor do the most massive stars reside in areas of enhanced stellar surface density compared to the average surface density, according to the $Sigma_{rm LDR}$ method. We then constrain the dynamical histories and hence initial conditions of both regions by comparing the observed values to $N$-body simulations at appropriate ages. Stars in both regions likely formed with sub-virial velocities which contributed to merging of substructure and the formation of smooth clusters. The initial stellar densities were no higher than $rho sim 100 - 500$M$_odot$pc$^{-3}$ for IC348 and $rho sim 500 - 2000$M$_odot$pc$^{-3}$ for NGC1333. These initial densities, in particular that of NGC1333, are high enough to facilitate dynamical interactions which would likely affect $sim$10 per cent of protoplanetary discs and binary stars.
The Q-parameter is used extensively to quantify the spatial distributions of stars and gas in star-forming regions as well as older clusters and associations. It quantifies the amount of structure using the ratio of the average length of a minimum spanning tree, mbar, to the average length within the complete graph, sbar. The interpretation of the Q-parameter often relies on comparing observed values of Q, mbar and sbar to idealised synthetic geometries, where there is little or no match between the observed star-forming regions and the synthetic regions. We measure Q, mbar, and sbar over 10 Myr in N-body simulations which are compared to IC 348, NGC 1333, and the ONC. For each star-forming region we set up simulations that approximate their initial conditions for a combination of different virial rations and fractal dimensions. We find that dynamical evolution of idealised fractal geometries can account for the observed Q, mbar, and sbar values in nearby star-forming regions. In general, an initially fractal star-forming region will tend to evolve to become more smooth and centrally concentrated. However, we show that initial conditions, as well as where the edge of the region is defined, can cause significant differences in the path that a star-forming region takes across the mbar-sbar plot as it evolves. We caution that the observed Q-parameter should not be directly compared to idealised geometries. Instead, it should be used to determine the degree to which a star-forming region is either spatially substructured or smooth and centrally concentrated.
Maser emission plays an important role as a tool in star formation studies. It is widely used for deriving kinematics, as well as the physical conditions of different structures, hidden in the dense environment very close to the young stars, for example associated with the onset of jets and outflows. We will summarize the recent observational and theoretical progress on this topic since the last maser symposium: the IAU Symposium 242 in Alice Springs.
Recent observations and hydrodynamical simulations of star formation inside a giant molecular cloud have revealed that, within a star forming region, stars do not form evenly distributed throughout this region, but rather in small sub-clumps. It is generally believed that these sub-clumps merge and form a young star cluster. The time-scale of this merging process is crucial for the evolution and the possible survival of the final star cluster. The key issue is whether this merging process happens faster than the time needed to remove the residual gas of the cloud. A merging time-scale shorter than the gas-removal time would enhance the survival chances of the resulting star cluster. In this paper we show by means of numerical simulations that the time-scale of the merging is indeed very fast. Depending on the details of the initial sub-clump distribution, the merging may occur before the gas is expelled from the newly-formed cluster either via supernovae or the winds from massive stars. Our simulations further show that the resulting merger-objects have a higher effective star formation efficiency than the overall star forming region and confirm the results that mass-segregated sub-clumps form mass-segregated merger-objects.