Statistical laws in chaotic dynamics of multiple stars

Statistical analysis of the modeled stellar systems dynamical evolution is performed. The initial global parameters (amount of stars, system size, virial ratio, mass spectrum) are varied. Final state distribution, final binaries and stable triples orbital elements are analysed. It is shown that the probability of the stable triple formation is rather high (about 10-15%). The eccentricity distribution of the final binaries satisfies the f(e) = 2e law. The hierarchy in the stable triple systems is rather strong (the mean ratio of the outer and inner binary sernimajor axes is about 20 : 1). In stable triples the eccentricities of internal binaries are in average greater than the ones of external binaries ((e(in)) over bar approximate to 0.7, (e(ex)) over bar approximate to 0.5). Stable triples with prograde motions are preferable.