Towards Multiple-Star Population Synthesis

The multiplicities of stars, and some other properties, were collected recently by Eggleton & Tokovinin, for the set of 4559 stars with Hipparcos magnitude brighter than 6.0 (4558 excluding the Sun). In this paper I give a numerical recipe for constructing, by a Monte Carlo technique, a theoretical ensemble of multiple stars that resembles the observed sample. Only multiplicities up to 8 are allowed; the observed set contains only multiplicities up to 7. In addition, recipes are suggested for dealing with the selection effects and observational uncertainties that attend the determination of multiplicity. These recipes imply, for example, that to achieve the observed average multiplicity of 1.53, it would be necessary to suppose that the real population has an average multiplicity slightly over 2.0. This numerical model may be useful for (a) comparison with the results of star and star cluster formation theory, (b) population synthesis that does not ignore multiplicity above 2, and (c) initial conditions for dynamical cluster simulations.