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.Comment: 13 pages, 10 figures; Accepted for publication in MNRA

The explosive BN/KL outflow emerging from OMC1 behind the Orion Nebula may have been powered by the dynamical decay of a non-hierarchical multiple system ∼500 years ago that ejected the massive stars I, BN, and source n, with velocities of about 10 to 30 km s −1 . New proper motion measurements of H 2 features show that within the errors of measurement, the outflow originated from the site of stellar ejection. Combined with published data, these measurements indicate an outflow age of ∼500 years, similar to the time since stellar ejection. The total kinetic energy of the ejected stars and the outflow is about 2 to 6 × 10 47 ergs. It is proposed that the gravitational potential energy released by the formation of a short-period binary, most likely source I, resulted in stellar ejection and powered the outflow. A scenario is presented for the formation of a compact, non-hierarchical multiple star system, its decay into an ejected binary and two high-velocity stars, and launch of the outflow. Three mechanisms may have contributed to the explosion in the gas: (i) Unbinding of the circumcluster envelope following stellar ejection, (ii) disruption of circumstellar disks and high-speed expulsion of the resulting debris during the final stellar encounter, and (iii) the release of stored magnetic energy. Plausible proto-stellar disk end envelope properties can produce the observed outflow mass, velocity, and kinetic energy distributions. The ejected stars may have acquired new disks by fallback or Bondi-Hoyle accretion with axes roughly orthogonal to their velocities. The expulsion of gas and stars from OMC1 may have been driven by stellar interactions.

We investigate the evolution, following gas dispersal, of a star cluster produced from a hydrodynamical calculation of the collapse and fragmentation of a turbulent molecular cloud. We find that when the gas, initially comprising ≈60 per cent of the mass, is removed, the system settles into a bound cluster containing ≈30–40 per cent of the stellar mass surrounding by an expanding halo of ejected stars. The bound cluster expands from an initial radius of <0.05 to 1–2 pc over ≈4–10 Myr, depending on how quickly the gas is removed, implying that stellar clusters may begin with far higher stellar densities than usually assumed. With rapid gas dispersal, the most massive stars are found to be mass segregated for the first ∼1 Myr of evolution, but classical mass segregation only develops for cases with long gas removal time‐scales. Eventually, many of the most massive stars are expelled from the bound cluster. Despite the high initial stellar density and the extensive dynamical evolution of the system, we find that the stellar multiplicity is almost constant during the 10 Myr of evolution. This is because the primordial multiple systems are formed in a clustered environment and, thus, by their nature are already resistant to further evolution. The majority of multiple system evolution is confined to the decay of high‐order systems (particularly quadruple systems) and the formation of a significant population of very wide (104–105 au) multiple systems in the expanding halo. This formation mechanism for wide binaries potentially solves the problem of how most stars apparently form in clusters and yet a substantial population of wide binaries exist in the field. We also find that many of these wide binaries and the binaries produced by the decay of high‐order multiple systems have unequal mass components, potentially solving the problem that hydrodynamical simulations of star formation are found to underproduce unequal mass solar‐type binaries.