No Arabic abstract
We study the effect of the instantaneous gas expulsion on star clusters wherein the residual gas has a density profile shallower than that of the embedded cluster. This is expected if star formation proceeds with a given SFE per free-fall time in a centrally-concentrated molecular clump. We perform direct N-body simulations whose initial conditions are generated by the program mkhalo falcON adapted for our models. Our model clusters initially have a Plummer profile and are in virial equilibrium with the gravitational potential of the cluster-forming clump. The residual gas contribution is computed based on the model of Parmentier&Pfalzner(2013). Our simulations include mass loss by stellar evolution and the tidal field of the Galaxy. We find that a star cluster with a minimum global SFE of 15% is able to survive instantaneous gas expulsion and to produce a bound cluster. Its violent relaxation lasts no longer than 20 Myr, independently of its global SFE and initial stellar mass. At the end of violent relaxation the bound fractions of surviving clusters with the same global SFEs are similar regardless of their initial stellar mass. Their subsequent lifetime in the gravitational field of the Galaxy depends on their bound stellar masses. We therefore conclude that the critical SFE needed to produce a bound cluster is 15%, which is twice smaller than earlier estimates of 33%. Thus we have improved the survival likelihood of young clusters after instantaneous gas expulsion. Those can now survive instantaneous gas expulsion with global SFEs as low as those observed for embedded clusters of Solar Neighbourhood (15-30%). This is the consequence of the star cluster having a density profile steeper than that of the residual gas. However, in terms of the effective SFE, measured by the virial ratio of the cluster at gas expulsion, our results are in agreement with previous studies.
We have studied the long-term evolution of star clusters of the solar neighborhood, starting from their birth in gaseous clumps until their complete dissolution in the Galactic tidal field. We have combined the local-density-driven cluster formation model of Parmentier & Pfalzner (2013) with direct N-body simulations of clusters following instantaneous gas expulsion. We have studied the relation between cluster dissolution time, $t_{dis}$, and cluster initial mass, $M_{init}$, defined as the cluster {mass at the end of the dynamical response to gas expulsion (i.e. violent relaxation), when the cluster age is 20-30 Myr}. We consider the initial mass to be consistent with other works which neglect violent relaxation. The model clusters formed with a high star formation efficiency (SFE -- i.e. gas mass fraction converted into stars) follow a tight mass-dependent relation, in agreement with previous theoretical studies. However, the low-SFE models present a large scatter in both the initial mass and the dissolution time, and a shallower mass-dependent relation than high-SFE clusters. Both groups differ in their structural properties on the average. Combining two populations of clusters, high- and low-SFE ones, with domination of the latter, yields a cluster dissolution time for the solar neighborhood in agreement with that inferred from observations, without any additional destructive processes such as giant molecular cloud encounters. An apparent mass-independent relation may emerge for our low-SFE clusters when we neglect low-mass clusters (as expected for extra-galactic observations), although more simulations are needed to investigate this aspect.
Base on Gaia Second Data Release and the combination of nonparametric bivariate density estimation with the least square ellipse fitting, we derive the shape parameters of the sample clusters. By analyzing the dislocation of the sample clusters, the dislocation $d$ is related to the X-axis pointing toward the Galactic center, Y-axis pointing in the direction of Galactic rotation, and the Z-axis (log(|H|/pc)) that is positive toward the Galactic north pole. This finding underlines the important role of the dislocation of clusters in tracking the external environment of the Milky Way. The orientation ($q_{pm}$) of the clusters with $e_{pm}$ $geq$ 0.4 presents an aggregate distribution in the range of -45$degr$ to 45$degr$, about 74% of them. This probably suggests that these clusters tend to deform heavily in the direction of the Galactic plane. NGC 752 is in a slight stage of expansion in the two-dimensional space and will deform itself morphology along the direction perpendicular to the original stretching direction in the future if no other events occur. The relative degree of deformation of the sample clusters in the short-axis direction decreases as their ages increase. On average, the severely distorted sample clusters in each group account for about 26% $pm$ 9%. This possibly implies a uniform external environment in the range of $|$H$|$ $leq$ 300 pc if the sample completeness of each group is not taken into account.
The early evolution of a dense young star cluster (YSC) depends on the intricate connection between stellar evolution and dynamical processes. Thus, N-body simulations of YSCs must account for both aspects. We discuss N-body simulations of YSCs with three different metallicities (Z=0.01, 0.1 and 1 Zsun), including metallicity-dependent stellar evolution recipes and metallicity-dependent prescriptions for stellar winds and remnant formation. We show that mass-loss by stellar winds influences the reversal of core collapse. In particular, the post-collapse expansion of the core is faster in metal-rich YSCs than in metal-poor YSCs, because the former lose more mass (through stellar winds) than the latter. As a consequence, the half-mass radius expands more in metal-poor YSCs. We also discuss how these findings depend on the total mass and on the virial radius of the YSC. These results give us a clue to understand the early evolution of YSCs with different metallicity.
We explore how the star formation efficiency in a protocluster clump is regulated by metallicity dependent stellar winds from the newly formed massive OB stars (Mstar >5 Msol). The model describes the co-evolution of the mass function of gravitationally bound cores and of the IMF in a protocluster clump. Dense cores are generated uniformly in time at different locations in the clump, and contract over lifetimes that are a few times their free fall times. The cores collapse to form stars that power strong stellar winds whose cumulative kinetic energy evacuates the gas from the clump and quenches further core and star formation. This sets the final star formation efficiency, SFEf. Models are run with various metallicities in the range Z/Zsol=[0.1,2]. We find that the SFEf decreases strongly with increasing metallicity.The SFEf-metallicity relation is well described by a decaying exponential whose exact parameters depend weakly on the value of the core formation efficiency. We find that there is almost no dependence of the SFEf-metallicity relation on the clump mass. This is due to the fact that an increase (decrease) in the clump mass leads to an increase (decrease) in the feedback from OB stars which is opposed by an increase (decrease) in the gravitational potential of the clump. The clump mass-cluster mass relations we find for all of the different metallicity cases imply a negligible difference between the exponent of the mass function of the protocluster clumps and that of the young clusters mass function. By normalizing the SFEs to their value for the solar metallicity case, we compare our results to SFE-metallicity relations derived on galactic scales and find a good agreement. As a by-product of this study, we also provide ready-to-use prescriptions for the power of stellar winds of main sequence OB stars in the mass range [5,80] Msol in the metallicity range we have considered
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 $muG$ show that only supercritical or marginal critical streams lead to reasonable star forming rates. 3) The GMCs 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)