A statistical model for the initial stellar mass function

The nonlinear stages of the fragmentation of a collapsing molecular cloud are modeled by coagulation theory. Several distinct physical processes are discussed, including protostellar fragment coalescence, gas accretion, and binary formation. This work confirms and extends the earlier result of Nakano (1966) that an approximately self-similar limiting form of the mass spectrum develops after several mean initial collision times. An approximate solution to the velocity-averaged coagulation equation is given for an arbitrary power-law dependence of the coalescence rate on mass, with dimensional dependence proportional to m to the power lambda; i.e., the asymptotic mass spectrum varies as m to the -3 lambda/2 power at small masses and cuts off exponentially at large masses, the characteristic mass depending on the number of collision times elapsed. Simple physical arguments suggest that lambda may increase with increasing mass, but is restricted to the range from 2/3 to 4/3. A large fraction of collisions could result in binary formation.