Massive stars are known to have a high multiplicity, with examples of higher order multiples among the nearest and best studied objects. In this paper we study hierarchical multiple systems (an inner binary as a component of a wider binary) of massiv
e stars in a clustered environment, in which a system with a size of 100--1000 au will undergo many close encounters during the short lifetime of a massive star. Using two types of N-body experiment we determine the post-formation collision probabilities of these massive hierarchies. We find that, depending on the specifics of the environment, the hierarchy, and the amount of time that is allowed to pass, tens of percent of hierarchies will experience a collision, typically between the two stars of the inner binary. In addition to collisions, clusters hosting a hierarchical massive system produce high velocity runaways at an enhanced rate. The primordial multiplicity specifics of massive stars appear to play a key role in the generation of these relatively small number events in cluster simulations, complicating their use as diagnostics of a clusters history.
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 $sim$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 $sim$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 times 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 circum-cluster 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 fall-back 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.
With the use of N-body calculations the amount and properties of escaping stars from low-N (N = 100 and 1000) young embedded star clusters prior to gas expulsion are studied over the first 5 Myr of their existence. Besides the number of stars also di
fferent initial radii and binary populations are examined as well as virialised and collapsing clusters. It is found that these clusters can loose substantial amounts (up to 20%) of stars within 5 Myr with considerable velocities up to more than 100 km/s. Even with their mean velocities between 2 and 8 km/s these stars will still be travelling between 2 and 30 pc during the 5 Myr. Therefore can large amounts of distributed stars in star-forming regions not necessarily be counted as evidence for the isolated formation of stars.
