the initial stellar mass function ( imf ) has been observed directly in dozens of clusters in our galaxy and the large magellanic clouds , and inferred indirectly for our galaxy and other galaxies from metallicities , supernova rates , ionization levels , luminosities , dynamical masses , colors , and other tracers .
there is some convergence in these observations to an imf that is approximately a power law above 1 m@xmath0 , having a slope between @xmath1 and @xmath2 on a @xmath3 histogram ( see reviews in massey 1998 and scalo 1998 ) , with a flattening below this down to @xmath4 m@xmath0 or lower .
the power - law part is in rough agreement with the salpeter ( 1955 ) function , which has a slope of @xmath5 on such a plot , and also with the scalo ( 1986 ) determination , which revised the often - used miller & scalo ( 1979 ) function at intermediate and high mass ( the miller - scalo function is too steep ) .
scalo ( 1998 ) recently emphasized the variations in the derived power law slopes from region to region , which amount to about @xmath6 .
there are no obvious correlations between these variations and physical properties of the star - forming regions , including the cluster density , galactocentric radius , and metallicity ( scalo 1998 ) .
there are also `` gaps '' in the power law parts of some cluster imfs , according to scalo ( e.g. , ngc 103 ; phelps & janes 1993 ) , and a `` ledge '' in the orion imf by hillenbrand ( 1997 ) .
these features led scalo to suggest that physical effects might be involved , making imf variations larger than poisson statistics .
an emphasis on the similarities , rather than the differences , in the imf has been made by massey and collaborators ( see review in massey 1998 ) .
they find an imf power law slope that is fairly constant , around the salpeter value , for clusters with a factor of @xmath7 range in densities ( massey & hunter 1998 ; luhman & rieke 1998 ) and a factor of @xmath8 range in metallicities ( massey , johnson & degioia - eastwood 1995a ; see also the globular cluster study by bellazzini et al . 1995 ) .
this slope was also found by kennicutt , tamblyn & congdon ( 1994 ) in galaxy - wide surveys , and by bresolin & kennicutt ( 1997 ) for individual hii regions in galaxies , using the equivalent widths of hydrogen emission lines , and it was found by greggio et al .
( 1993 ) , marconi et al .
( 1995 ) , holtzman , et al .
( 1997 ) , and grillmair , et al .
( 1998 ) from color magnitude diagrams of stars in the lmc and local dwarf galaxies .
chemical evolution also tends to suggest an invariant imf , at least at intermediate and high mass , in the stellar halo ( nissen et al .
1994 ) and disk ( tsujimoto et al . 1997 ; wyse 1998 ) , in qso damped ly@xmath9 ( lu et al .
1996 ) and ly@xmath9 forest ( wyse 1998 ) lines , the intracluster medium ( renzini et al . 1993 ; wyse 1997 , 1998 ; however , see loewenstein & mushotzky 1996 ) , and elliptical galaxies ( see review in wyse 1998 ) .
the observable in this case is the ratio of the iron to the oxygen abundance , which is always about the same .
since this ratio scales with the ratio of intermediate mass to high mass stars , the imf must be about the same too .
a different result was obtained by heikkil , johansson & olofsson ( 1998 ) , who conjectured that a low c@xmath10o / c@xmath11o ratio in the lmc results from a relatively low proportion of high - mass stars , as observed , for example , by massey et al .
( 1995b ) ; clearly the chemical and fractionation processes involved with this important isotope ratio have to be investigated further .
gibson & matteucci ( 1997 ) also obtained a different result from evolution and metallicity models of elliptical galaxies , including infall and wind ejection of gas .
they suggested the imf had a higher proportion of massive stars than the solar neighborhood , as it appears to in some starburst models ( discussed below ) .
angeletti & giannone ( 1997 ) modeled elliptical galaxies and fit the mass - to - light and [ fe / h ] ratios with an imf that has a mass - dependent high mass slope , different from the salpeter value .
the imf in local ob associations might differ from the salpeter value too .
brown ( 1998 ) reviewed observations , including those with distances from the hipparcos satellite , which suggest a slightly steeper slope , @xmath12 to @xmath13 for local associations , but the stellar types in these studies were obtained from photometry , and massey ( 1998 ) discussed how such data can artificially steepen the imf without spectroscopic observations too .    at the low mass end
, there is no _
systematic _ indication yet of a turnover at masses lower than 0.1 m@xmath0 , but there are some regions where such a turnover has been reported , e.g. , in nearby field stars ( reid & gazis 1997a ) , orion ( hillenbrand 1997 ) , the pleiades ( reid 1998 ) , ngc 6231 ( sung , bessell , & lee 1998 ) , 30 dor ( nota 1998 ) , and the globular cluster ngc 6397 ( king et al .
there are also regions where the imf appears to rise below 0.1 m@xmath0 , again for local stars ( mra et al .
1996 ; gould et al .
1997 ) . scalo ( 1998 ) suggests that the local field results at low mass contain uncertainties in the star formation rate , pre - main - sequence contribution , and mass - luminosity relation .
the older cluster results are uncertain also because of evaporation of low mass members , although reid ( 1998 ) suggests evaporation is not important for the pleiades . in any case , scalo ( 1998 ) and reid ( 1998 ) both provide evidence that the low mass imf , below 1 m@xmath0 , varies slightly in shape from region to region , although no sensible correlations with properties of the environment have yet been noticed .
variations in the low mass imf are also seen for globular clusters , without any obvious dependence on metallicity ( santiago et al .
1996 ; piotto et al .
1997 ; see review in cool 1998 ) , but there may be some dependence on the cluster evaporation rate and susceptibility to orbit shocking ( cool 1998 ) , since low mass stars leave the cluster most easily .
the imf for the young open cluster mentioned above , ngc 6231 ( sung et al . 1998 ) , is somewhat peculiar , because the turnover mass is around 2.5 m@xmath0 , much higher than in ob associations and other clusters . however , ngc 6231 is a dense cluster @xmath8 my old , and a high turnover mass in the central regions might be expected from rapid thermalization leading to mass segregation ( e.g. , scalo 1986 ) .
the low mass imf is also probed by observations of embedded clusters , where most of the stars are pre - main sequence ( zinnecker et al .
1993 ; fletcher & stahler 1994a , b ; lada & lada 1995 ; lada , alves & lada 1996 ; megeath 1996 ; horner , lada , & lada 1997 ) .
pre - main - sequence ( pms ) clusters have not yet evolved kinematically , so differential evaporation and segregation of high and low mass members should not be as much of a problem as it is for older clusters .
pms stars are also much brighter than they will be on the main sequence , so even stars less massive than the hydrogen burning limit can be observed in local clusters ( see review in lada , lada , & muench 1998 ) .
lada et al .
( 1998 ) combined herbig s ( 1998 ) optical data on the cluster ic 348 , which was used to determine the history of star formation in the region , with their own k - band luminosity function to determine the low mass imf .
they found that the imf in this cluster turns over below 0.25 m@xmath0 .
recent searches for brown dwarfs ( zuckerman & becklin 1992 ; reid & gizis 1997b ; reid 1998 ) suggest these objects are so infrequent that the imf can not continue to rise at low mass .
the number of possible brown dwarfs is consistent with an extension of the flattening seen between 1 m@xmath0 and 0.1 m@xmath0 ( reid 1998 ) .      in spite of these approximate similarities in the imf from region to region , there have been two clear and consistent deviations that , if real , would have important consequences for star formation .
one is the apparent shift toward relatively higher masses in regions of extreme activity , such as starburst galaxies ( e.g. rieke et al .
1980 ) , and the other is the apparent shift toward relatively lower masses in regions of extreme _
in_activity , the so - called `` field '' population , away from ob associations , as discussed by massey et al .
( 1995b ) .
indeed , there may be a monotonic transition from low - mass weighted imfs to high - mass weighted imfs with increasing star formation rate , ism pressure , radiation field , etc .
, perhaps because of a change in cloud properties , cloud destruction capabilities of young stars , and other things .
the `` normal '' imf occupies what seems to be a broad range of star formation rates between these two limits , but this broad range is not so broad in terms of _ possible _ ism conditions .
it may only be broad in the sense that it is common in main galaxy disks where some type of feedback between star formation and gas stability is maintained for long periods of time .
reviews of imfs in starburst galaxies are in scalo ( 1990 ) , zinnecker ( 1996 ) , and leitherer ( 1998 ) . the first observation that suggested a `` top - heavy '' imf was that the luminosity of the starburst regions are so high , and the dynamical masses from the rotation curves so low , that there can not be a full complement of low mass stars , which dominate the mass , to go with the high mass stars , which dominate the light ( rieke et al . 1980 , 1993 ; kronberg , biermann , & schwab 1985 ; wright et al .
1988 ; see reviews in telesco 1988 , joseph 1991 , rieke 1991 , thuan 1991 ) .
thus the imf contains a higher proportion of high mass stars than the solar neighborhood .
doane & matthews ( 1993 ) reached the same conclusion for m82 based on many properties of galactic evolution models , including the remaining gas mass , total mass , star formation rate , and supernova rate .
doyon , joseph , & wright ( 1994 ) and smith et al .
( 1995 ) found top - heavy imfs from spectroscopic line ratios for the starburst galaxies ngc 3256 and ugc 8387 , respectively .
smith , herter & haynes ( 1998 ) found the same result from infrared excesses in 20 starburst galaxies .
the distribution of stars was often not resolved in the early starburst studies .
satyapal et al . ( 1995 , 1997 ) derived a smaller extinction - corrected k band luminosity for the inner 30  of the starburst galaxy m82 , and suggested the imf could have the salpeter slope all the way from 0.1 to 100 m@xmath0 with only 30% of the dynamical mass in the star formation event . ho & filippenko ( 1996a , b ) measured the velocity dispersions in two dense clusters in starburst galaxies , clusters that may evolve into globulars like those in the halo of our galaxy ( meurer 1995 ) .
they found that the young globulars have large masses like old globulars , and that the young clusters would evolve into old clusters with a mass - to - light ratio like the old globulars after a hubble time .
because of this , they suggested that the imf might be the same in the young starburst globulars as it is in old globular clusters , and in particular , that it contains low mass stars ( e.g. , de marchi & paresce 1997 ) .
other considerations have led to revisions of the top - heavy models as well .
devereux ( 1989 ) observed 20 nearby starbursts like m82 with infrared photometry , and found no evidence for an unusual imf in any of these ; he pointed out that a revision of extinction corrections for m82 could change the rieke et al .
( 1980 ) result too , as recently confirmed by satyapal et al .
schaerer ( 1996 ) considered detailed evolutionary models of young stellar populations in starburst galaxies and also found that the imf was normal .
calzetti ( 1997 ) modeled multiwavelength spectroscopy and broad - band infrared photometry of 19 starburst galaxies and found that the salpeter imf between 0.1 and 100 m@xmath0 fits the observations .
stasiska & leitherer ( 1996 ) modeled the emission line spectra of giant hii region and starburst galaxies , having a factor - of - ten range in metallicities , and found a salpeter imf up to 100 m@xmath0 ; the lower mass limit could not be obtained , however .
the proposed peculiarity of starburst imfs is therefore not certain , a conclusion reached by scalo ( 1990 ) as well .
maybe only the extreme examples of starbursts require a high proportion of massive stars , while other starbursts , perhaps more `` ordinary , '' do not .
the exact form of the peculiarity in a starburst imf is not known either .
it could be `` truncated '' in the sense that the lower mass limit of the power law , where the power law slope begins to flatten , is higher than usual , perhaps as high as several m@xmath0 instead of several tenths of a m@xmath0 .
such truncation may be observed already for 30 dor ( nota 1998 ) . or
, it could have a power - law slope that is shallower than @xmath14 , or possibly , a higher upper mass limit .
these differences are often hard to distinguish , even though all of them affect the mass - to - light ratio .
fortunately , the exact form of the imf sometimes affects other things as well .
for example , charlot et al .
( 1993 ) suggested that a truncated imf will produce a very red population of red giants , without the corresponding main sequence stars , after the turnoff age reaches the stellar lifetime at the truncation mass .
wang & silk ( 1993 ) suggested that the truncated model gives an oxygen abundance that is too high when the star formation process is over .
presumably these other constraints will help narrow down the possible choices .    the extreme field studied by massey et al .
( 1995b ) is another example where the salpeter imf seems to break down . in regions of the lmc and our galaxy that are more than 30 parsecs from catalogued ob associations , massey et al
. found o - type stars that probably did not drift there during their short lifetimes .
massey et al
. showed with an h - r diagram that there was an approximately steady rate of star formation there for the last 10 my years , and they derived a time - averaged imf .
the result was very steep , with @xmath15 instead of @xmath16 for a salpeter function .
there are several other peculiarities of the imf that may or may not have a specific physical origin .
for example , it has been known for a long time ( larson 1982 ; myers & fuller 1993 ) that high mass clouds produce high mass stars , and low mass clouds produce low mass stars .
two decades ago , this observation and others were related to the concept of a bimodel imf ( eggen 1976 , 1977 ; elmegreen 1978 ; gsten & mezger 1983 ; larson 1977 , 1986 ) , and to the impression that spiral arms preferentially formed massive stars ( mezger & smith 1977 ) , presumably because the molecular clouds in spiral arms are more massive , or somehow different from interarm clouds . today , we know that massive clouds also produce low mass stars in great abundance ( herbig & terndrup 1986 ) , and there are essentially no local regions that produce high mass stars exclusively . thus the `` observation '' that the maximum mass of a star correlates with the mass of the cloud could be only a statistical sampling effect ( larson 1982 ; elmegreen 1983 ; walter & boyd 1991 ; massey & hunter 1998 ) : high mass clouds produce more stars of all types , and so they are more likely to sample at the tail of the imf to produce a few massive stars . in the same sense , high mass clouds should also produce more dwarf stars than low mass clouds .    an imf model based on random sampling from a hierarchical cloud ( elmegreen 1997a ) shows this effect clearly in figure [ fig : mvsm ] .
each curve on the left is from a different cluster , with a number of stars in the clusters equal to 2500 , 5000 , 10000 , 25000 , 50000 , and 100000 .
as the number of model stars increases , the mass of the largest star increases too , as shown on the right of the figure .
the expected slope of the power law relation between cluster mass , @xmath17 , and maximum star mass , @xmath18 , is @xmath19 , i.e. , @xmath20 ; actually the calculation gives a slightly steeper slope because of storage limitations to the number of hierarchical levels allowed in the simulation ( these models use 8 levels with an average of 3 subclumps per clump at each level ) .
this property of increasing maximum star mass with cloud mass should not depend on any specific imf model or on the star formation process .    to be more precise about the stochastic interpretation of this effect
, we can consider an imf @xmath21 , which gives the salpeter function when @xmath22 .
the largest star , with mass @xmath18 , typically satisfies @xmath23 .
this gives @xmath24 .
because the total cluster mass is @xmath25 , the largest star is related to the total cluster mass as @xmath26.\ ] ] for @xmath27 and @xmath22 , this gives @xmath28 , and for @xmath29 @xmath30 , this gives @xmath31 .
thus , a cluster requires several hundred m@xmath0 before a small o - type star is likely to appear , and a cluster with @xmath32 m@xmath0 is likely to contain @xmath8 modest o - type stars or a few very massive ones .
this is the mass range for the largest open clusters in the milky way disk , indicating that when a few o stars appear , a cluster is likely to terminate its formation in an unbound state , presumably as a result of rapid cloud destruction , and make an expanding ob association ( elmegreen 1983 ) .
larson ( 1991 , 1992 ) and khersonsky ( 1997 ) interpret this observation in the opposite way , saying that there is a physical effect which tends to produce higher mass stars in clouds of higher mass , and that the imf partly results from this physical effect .
larson ( 1992 ) attributes the physical effect to a larger accretion length in larger clouds for material going into a star .
if there is such an effect , and the observation is not entirely statistical as discussed above , then the superposition of imfs from clouds of different mass would be different from the imf in any one cloud .
that is , clouds of low mass would produce an imf with a particular slope up to some maximum mass , and clouds of higher mass would produce an imf , perhaps with the same slope , up to higher masses .
the sum of these two imfs would be an imf with greater slope than either separate cloud , because the summed imf contains low mass stars from both clouds , but only high mass stars from the more massive cloud .
larson ( 1991 ) proposed his model by saying , in effect , that only one star forms in each cloud fragment at each level in the hierarchy , and that the star mass is always proportional to the square root of the fragment mass .
then the sum over all fragments is the final imf , and this has a slope of @xmath33 , which is not too bad .
stated this way , there is no problem with the model because each component in the sum is only a single star , not a whole cluster of stars with a separate imf .
however , even in this model , if we add together the stellar populations from many separate clouds with different total masses , then the summed imf will be steeper than the imf in each . at this point , it is important to recall that the imf obtained from large - scale surveys of galaxies and their metallicities is about the same as the imf observed in individual clusters .
thus the summed imf from many regions of star formation has to be about the same as the imf from each separate cluster .
this would seem to suggest that the observed increase in maximum star mass with cloud mass is not the result of a physical process that specifically increases the upper stellar mass in larger clouds .
a second peculiarity about star formation is that the high mass stars tend to form after the low mass stars ( herbig 1962a , b ; iben & talbot 1966 ) .
this has been determined most recently for 30 dor ( massey & hunter 1998 ) .
again one can think of physical reasons for this , such as a gradual warming of the cloud following the formation of low mass stars , and an increase in the jeans mass with temperature ( silk 1977 ; yoshii & saio 1985 ) , but statistical effects can explain it too .    the decreasing nature of the imf at high mass implies that massive stars are likely to form only after a lot of low mass stars have already formed . for a constant star formation rate ,
the average time between the formation of stars in a logarithmic mass range centered on @xmath34 is proportional to the rate at which these stars form , which is inversely proportional to the relative number of the stars , or @xmath35\right)^{-1}\propto m^{x}$ ] . for a constant star formation rate , this is also the average time after star formation begins for the _ first _ appearance of a star with this mass .
thus the most common , low - mass stars form first , followed by the intermediate and then the high mass stars . in terms of the proportion of all stars formed , and for a constant star formation rate in units of mass per year , the maximum stellar mass increases with time as @xmath36 for star formation rate @xmath37 .
thus , a star of mass @xmath38 forms after the proportional time given by @xmath39 for the power law portion of the imf only , with @xmath40 m@xmath0 and @xmath22 , this becomes , for masses in m@xmath0 , @xmath41 for example , in a cluster with a maximum stellar mass of 30 @xmath42 , stars with 10 m@xmath0 form in the last 20% of the time .
note that this increase in time with stellar mass is only for the first appearance of a star with that mass .
the _ average _ time of appearance is independent of stellar mass in the stochastic interpretation .
this suggests a way to check the stochastic model of birth orders .
if the average time of appearance of all stars with mass @xmath34 is the average age of the cloud , and this is true for any @xmath34 , and if the first time of appearance of a star with mass @xmath34 increases with time , as described above , then the birth order is dominated by stochastic effects .
figure [ fig : mt ] shows the relative birth order of stars in the random sampling imf model by elmegreen ( 1998 ) .
each point marks the earliest time of formation of a model star with that mass , among all of the stars that formed in 11 clusters on the left and 6 clusters on the right .
each cluster contains 200 to 500 stars and is independently generated .
the plot shows a distribution of points that traces the inverse of the imf in log intervals .
the fiducial line has a slope of @xmath43 for comparison .
the increase at low mass occurs because the model uses a gaussian probability for failure to form a star at low mass , which is an arbitrary assumption at this time ( it causes the model imf to turn over at low mass ) .
there are no quantitative measures yet of how the earliest appearance of a star with a particular mass increases with time relative to the birth order of all stars , so the predictions of theory can not yet be checked . in a real star - forming region , the star formation rate is probably not constant in time , so the ordinate in figure [ fig : mt ] should not be time , but star number in its proper birth order ; i.e. one should plot , for the first appearance of a star of mass @xmath34 , its mass versus the number of stars that formed before it
.    an interesting implication of the above two discussions is that _ massive stars must not be able to severely limit or halt star formation in their primordial clouds . _
if they did , then , considering their relatively late appearance , each cloud would produce stars with the same local imf up to some maximum stellar mass at which point the cloud is finally destroyed .
this maximum stellar mass at the time of cloud destruction is presumably larger for more massive clouds , because it is harder to destroy a massive cloud than a low mass cloud .
thus the summed imf from many regions of star formation would be steeper than the local imf in each region , for the reasons discussed above , and we could not explain why the integral imf in a galaxy is the same as the cluster imf .
perhaps massive stars only shred the low density parts of these clouds , allowing ( or perhaps stimulating ) , most of the cores that have already formed in these regions to continue their evolution toward stars . in view of this
, it seems that stars of all mass are probably equally likely to form in any cloud ( that contains at least this much mass ) , regardless of what other stars have formed before them .
moreover , the maximum stellar mass is not physically limited , it is only statistically limited in the sense that extremely massive stars are rare .
an exception to this conclusion may arise in the `` field '' population studied by massey et al .
( 1995b ) , in which the summed imf of presumed dispersed clusters appears to be steeper than the imf in each .
maybe massive stars can more easily disrupt their clouds in low pressure `` field '' environments . on the other hand , maybe the massey et al .
observation can be explained in other ways .
a third peculiarity of the imf is that massive stars are often closer to a cluster center than low mass stars .
this effect has been discussed for a long time ( e.g. , sagar et al . 1988 ) , but it was never certain whether the observations reflected a birth position , or just a relaxation effect after the cluster formed .
thermal relaxation makes the high mass stars sink to the center of the cluster , on average , in a thermalization time scale .
however , hillenbrand & hartmann ( 1998 ) and bonnell & davies ( 1998 ) have recently shown that this peculiar mass distribution , which is observed in the trapezium cluster , could not result from thermalization because the orion cluster is too young ( see also jones & walker 1988 ) .
similar mass segregation has been observed in ngc 2157 ( fischer et al .
thus it is an effect that should be explained by an imf theory , and it is likely to result from a physical process during star formation rather than stochastic sampling .    there have been many attempts to explain this relative birth position ( see reviews in larson 1991 , and zinnecker , mccaughrean & wilking 1993 ) , including excess gas clump collisions in the center of a cloud ( larson 1990 ; stahler , palla , & ho 1998 ) , excess gas accretion onto protostars in the cloud core ( larson 1978 ; larson 1982 ; zinnecker 1982 ; bonnell et al .
1997 ) , enhanced gas drag for massive protostars ( larson 1990 , 1991 ) , late - time formation of high mass stars in a collapsing cloud ( murray & lin 1996 ) , and hierarchical cloud structure in which the massive trunk of a hierarchical tree for the gas cloud is closer to the center of the tree than the low - mass branches ( elmegreen 1998 ) . a combination of these effects might be involved .
there are evidently four distinct physical effects that have to be explained by a theory of the imf : + ( 1 ) the power - law slope at intermediate to high stellar mass , with a value close to that found by salpeter ( 1955 ) ; + ( 2 ) the flattening at low mass for clusters and star - forming regions in our galaxy ; + ( 3 ) the increase in the transition mass between power - law and flat regions of the imf in starburst galaxies , + ( 4 ) the preferential birth of high mass stars close to cluster centers .
there are apparently three additional observations that can be explained by stochastic sampling effects : + ( 1 ) the tendency for the largest mass star in a region to increase with the total number of stars ; + ( 2 ) the tendency for the largest mass stars to form last , and + ( 3 ) the seemingly random variations in the intermediate and high mass imf slopes from region to region , along with the appearance of gaps , ledges , and other peculiar departures from a power law .
in addition to these observations , there are tentative indications that other imf features may be present , including : + ( 1 ) a turnover at low mass , and + ( 2 ) a steepening of the power - law part in field regions .
+ future observations should be able to illuminate these uncertain features .
the theory of the imf depends on both the theory of star formation and the theory of cloud structure . because of the complexity of these issues
, the imf may not be understood for a long time .
here we review some recent developments .
previous reviews were in cayrel ( 1990 ) and clarke ( 1998 ) .
one of the key problems in understanding the origin of the imf is to determine the processes that limit the pre - stellar gas accretion onto a star and define the star s final mass .
there have been several ideas on this .
the final mass could result from the star s own ability to limit the accretion of new gas onto its surface , perhaps because of an intense proto - stellar wind ( larson 1982 ; shu , adams , & lizano 1987 ) , which is known to exist at this phase of evolution ( lada 1985 ) .
this possibility has led to several theories of the imf .
one , proposed by nakano , hasegawa , & norman ( 1995 ) , considered three mass - limiting agents : protostellar winds , ionization , and depletion of the gas clump in which the star forms . they concluded that under normal conditions , the protostellar wind would limit the accretion .
then they calculated the fraction of the clump that would get into the star , finding that the final star mass scaled with the clump mass to the @xmath44 power . to get an imf slope of @xmath45 , which they got from the field - star imf in scalo ( 1986 )
, they required a very steep clump mass spectrum , @xmath46 in linear intervals , instead of the usual observation of @xmath47 with @xmath48 ( blitz 1993 ; kramer et al .
they justified this steep clump spectrum by noting that if the linewidths found in a cs survey of orion clumps ( tatematsu et al .
1993 ) were assumed to be virial velocities , and the masses calculated accordingly , then the orion cs clumps would have such a spectrum .    another imf model with wind - limited stellar masses was proposed by adams & fatuzzo ( 1996 ) .
they conjectured that the stellar luminosity scales with the wind mass loss rate , which , at the time when the accretion stops , is proportional to the direct accretion rate onto the star , independent of what goes onto the disk .
this gave them a relation between the stellar properties , i.e. , luminosity and mass , and the gas cloud properties , sound speed and angular rotation rate .
the angular rotation rate entered because they had to determine the fraction of the accreting gas that goes into the star and not the disk . finally , they related the stellar mass to the total luminosity of the star , from the sum of the accretion - driven luminosity , which depends on mass , and the luminosity of a main sequence star at intermediate mass , from the standard mass - luminosity relation .
the result is a relation between stellar mass and the cloud properties , i.e. , sound speed and rotation rate .
taking the mass distribution for clumps with the observed slope @xmath9 discussed above , and the empirical relationship between clump mass and velocity dispersion from larson s ( 1981 ) laws , they got a distribution function for the clump velocity dispersion , which then led to a distribution function for the final star mass if the angular rate of all the clumps is the same .
adams & fatuzzo ( 1996 ) also discussed more general equations giving the stellar mass from cloud properties , considering random variations in these cloud properties , and derived a log - normal final mass distribution , as in larson ( 1973 ) , elmegreen & mathieu ( 1983 ) , zinnecker ( 1984 ) , and elmegreen ( 1985 ) .
silk ( 1995 ) also got a relation between cloud core mass and turbulent linewidth , different from adams & fatuzzo s , considering centrifugally supported cores with luminosities equal to the accretion luminosities .
he considered an equality between the turbulent linewidth and the thermal velocity in the core , which entered into the luminosity through the temperature using the usual radiative transfer equation for a star and the rossland mean opacity .
then , after substituting the empirical correlations between clump rotation rate and size , and between linewidth and size , he got the cloud core mass as a function of the turbulent linewidth .
the distribution function for cloud core mass then followed from a distribution function for linewidth , which came from a theory for the time - dependent deceleration of wind - driven bubbles in a cloud .
the star mass is then taken to be proportional to the cloud core mass .
silk ( 1995 ) made several comments that protostellar outflows limit the mass of a star , but this assumption was not present in any of the theory in his paper , nor in the imf that resulted , particularly considering that a fixed fraction of the cloud core mass was assumed to go into the star .
this was unlike the results of the wind - limited accretion models proposed by nakano , hasegawa , & norman ( 1995 ) and adams & fatuzzo ( 1996 ) .
one of the most popular and persistent methods for determining the mass of a star has been with a combination of cloud fragmentation , accretion , and clump collisions .
larson ( 1978 ) showed with three - dimensional n - body experiments that gas clouds fragment hierarchically by self - gravity , and the fragments accrete material in competition with each other .
the resulting mass distribution was modeled after a fractal , which was a remarkably prescient concept in astronomy considering that mandelbrot ( 1977 ) began to popularize his fractal geometry only a year earlier .
larson showed that the fractal mass function would have a slope of about @xmath49 , in reasonable agreement with observations .
larson also reasoned that the lower mass limit to a star is the thermal jeans mass in the cloud because smaller condensations are not likely to form in the initial collapse .
this experiment was an important departure from standard fragmentation models of the time , since it was previously believed , following hoyle ( 1953 ) and hunter ( 1962 ) , that fragmentation decreased the jeans mass , leading to ever more fragmentation . in larson
s result , the number of clumps formed was about equal to the number of jeans masses in the initial cloud , with no subfragmentation into smaller pieces .
this result was seconded by tohline ( 1980 ) , but investigated again by silk ( 1982 ) , who concluded that initial cloud fragments , particularly elongated fragments , formed by dynamical collapse could in fact fragment again ( see also bonnell & bastien 1993 , burkert et al .
1997 ) .    silk ( 1977 ) considered a different picture of cloud fragmentation , using a thermal jeans mass that increased with time as a result of heating from the stars that already formed .
his model of fragmentation was simpler than larson s , so there was no built - in , power - law mass spectrum from the fragmentation process itself .
the power law in silk s model came from the time increase in the thermal jeans mass and the identification of this mass with the stellar mass at all times .
yoshii & saio ( 1985 ) followed this model by considering the additional influence of coalescence among the opacity - limited fragments .
they showed that the opaque fragments are usually so small that coalescence is unimportant ; thus fragmentation and heating alone determine the imf as a time sequence of increasing thermal jeans masses , as in the silk ( 1977 ) model .
bastien ( 1981 ) obtained a different result by using the jeans length to determine the fragment collision cross section .
the initial jeans length is much larger than the size of the opaque fragments considered by yoshii & saio ( 1985 ) , so bastien ( 1981 ) found that fragment collisions were important .
subsequent work by lejeune & bastien ( 1986 ) , and allen & bastien ( 1995 , 1996 ) considered time - dependent coalescence , and concluded again that both fragmentation and coalescence were important , with coalescence dominating the formation of massive stars .
price & podsiadlowski ( 1995 ) used protostar collisions in a different way .
they proposed that stars grow by accretion at some more - or - less uniform rate , and this accretion stops when two protostars collide , disrupting the gas reservoirs from each .
the final stellar mass was then determined by the product of the accretion rate and the time interval between collisions .
the imf was built up over time as the final stellar mass decreased in the presence of an increasing collision rate that resulted from the continuous formation of more and more protostars .
another coagulation model was proposed by murray & lin ( 1996 ) .
they suggested that fragmentation driven by thermal instabilities leads to clumps that fall in the potential well of a cloud , after which they collide and coalesce with other clumps that have fallen too .
a star forms when the clump mass exceeds the thermal jeans mass , although further coalescence can increase that mass afterwards .
the power - law distribution of masses follows from the coagulation process , in a manner similar to that proposed by silk & takahashi ( 1979 ) .
there is no assumption about wind - limited accretion here ; the star mass is identified with the core mass directly .
an interesting aspect of their model is that they considered positionally correlated velocities , as in a turbulent fluid , but they showed this had no important effect on the model imf .
numerical sph simulations of an evolving imf made of accreting , pseudo - star particle clumps was done by bonnell et al .
the advantages of numerical solutions like this is that they can treat well the competition for gas among all the nucleated centers .
they found that nucleating centers originally close to the center of the whole cloud grow faster and to larger masses because of the larger gas density there .
a different class of theory considers only the statistical aspects of the imf .
the hoyle ( 1953 ) picture of fragmentation led naturally to these ideas , because successive fragmentation produces hierarchical structure and power law mass spectra regardless of many physical details .
larson ( 1973 ) considered a modification to these ideas by proposing that only part of each fragment undergoes further fragmentation , and got from this a log - normal mass distribution instead of a power - law .
when miller & scalo ( 1979 ) proposed that the imf actually had a log - normal form , the statistical implications of this were developed further by elmegreen & mathieu ( 1983 ) , zinnecker ( 1984 ) , and elmegreen ( 1985 ) , who showed that a variety of processes , acting together , would combine to make a log - normal function , regardless of the distribution functions resulting from each separate process .
this point was made again more recently by adams & fatuzzo ( 1996 ) .
the most recent development in this area has been by elmegreen ( 1997a , 1998 ) , who considered the random selection of masses from a hierarchically structured cloud .
these papers depart from the usual scenarios by asserting that hierarchical cloud structure is independent of star formation  that it is set up long before star formation begins in the diffuse cloud stage , and then continuously reestablished during the molecular cloud stage as a result of supersonic turbulence compression . there is no gravitationally - driven fragmentation at all , and no significant clump coalescence .
in fact , the clumps in this model need not even be constant objects , able to move around and coalesce ; they can be ever - changing gas compressions and wavepackets in the chaotic turbulent flow .
this model is motivated by the pervasive appearance of hierarchical structure and spatially - correlated velocities in interstellar clouds .
these features are reminiscent of structures and flows in laboratory turbulence ( e.g. , see reviews in sreenivasan 1991 , and falgarone & phillips 1991 ) .
the basic point of the elmegreen ( 1997a , 1998 ) model is that essentially all local star formation processes that have been considered in detail have a time scale for evolution that varies approximately as the inverse square root of local density , including contraction or collapse from self - gravity , turbulence compression , magnetic diffusion in virialized cores , and coalescence at each level in the hierarchy .
this means that in a model like this , stars will appear here and there , randomly , alone or with neighbors , at a rate that depends almost exclusively on the local square root of density . on a large scale
, this implies that dense clouds in high pressure regions of galaxies will form all of their stars quickly , while low density clouds in low pressure regions will take much more time . on a small scale , within any one cloud , it means that different parts of the hierarchy form stars at different times .
the density always increases at lower levels in a hierarchical structure , and this is where the clumps , contained by other clumps on larger scales , are small and have low mass .
thus the low mass pieces in a cloud are likely to form stars first . as a result ,
there is slightly less mass for other stars that form later in the same or higher levels .
this competition for mass tends to steepen the power in the power - law mass function by several tenths , i.e. , from @xmath49 to @xmath22 , as in the salpeter imf .
perhaps a more important difference in this new model is the way in which cloud structure is measured .
typically , the power law index for molecular and diffuse cloud structure is determined from the mass distribution of separate clumps , resolved in large - scale surveys by various telescopes and then separated into discrete objects either by eye or by computer algorithms that do about the same thing as the eye .
these mass distributions are always much flatter than the imf , having slopes @xmath9 in the range from 1.5 to 1.8 , when the imf has a slope of about @xmath50 . to understand the imf , however , we have to try to view interstellar clouds from the perspective of a forming star
. a star does not care about the clumps that our telescopes resolve and our eyes choose to label as discrete , but only about the general distribution and motion of gas in all forms . in pre - star - forming clouds , where this structure is first established ,
it is largely hierarchical , from scales that are much smaller than mm - wave telescopes can resolve , up to perhaps a galactic scale height .
hierarchical means that most of the small clouds are contained inside larger clouds , and most small stellar groupings are inside larger stellar groupings ( see review in elmegreen & efremov 1998 ) .
for example , efremov ( 1995 ) and colleagues have estimated that 90% of the ob associations in the milky way ( efremov & sitnik 1988 ) , m31 ( efremov , ivanov , & nikolov 1987 ; battinelli 1991 , 1992 ; magnier et al .
1993 ; battinelli , efremov & magnier 1996 ) , m33 ( ivanov 1987 , 1992 ) , and the lmc ( feitzinger & braunsfurth 1984 ) are inside larger star complexes .
scalo ( 1985 ) has reviewed hierarchical structure for interstellar clouds .
it appears as if the hierarchical embedding of interstellar and young stellar structures is nearly all - inclusive .
emission line surveys leading to `` discrete '' clump masses do not consider this aspect of cloud structure .
if two clumps are close together and part of a larger structure , the algorithms call them two separate objects , and do not tabulate the larger `` object '' that contains them . for this reason , all emission line or extinction surveys catalogue clouds that are within a factor of @xmath51 to 10 of the angular resolution , regardless of the cloud distance or the wavelength of the observation
. this factor of @xmath52 for recognized clump size corresponds to a factor of @xmath53 for cloud mass ( which scales as the size squared or cubed in co surveys ) , so the mass functions look reasonably well sampled , but in fact only a small part of the cloud structure is included .
everything smaller than the telescope resolution is not seen , and everything larger than about 3 to 10 resolution elements is subdivided into its component parts and not called a separate object .
this is how every cloud or clump mass function has been evaluated since the beginning of this exercise ( i.e. , since before field & saslaw 1965 ) .
the elmegreen ( 1997a , 1998 ) model takes a different point of view .
it considers the structure on all scales and asks for the probability that any particular mass is chosen from anywhere in the whole hierarchy .
this is presumably what happens as a result of the combination of physical processes that leads to star formation in a real cloud : because of the self - similar nature of turbulent flows , each level in the hierarchy looks the same as any other level ( for masses above the thermal jeans mass ) , so each choice of level for a star - formation event would have the same likelihood as any other choice , modulated only by the density - dependent local evolution rate discussed in the previous paragraph . aside from this variable rate ,
the instantaneous distribution of masses for all structures , and the instantaneous probability of selecting any particular mass , is proportional to @xmath54 for linear intervals in mass ( i.e. , @xmath55 in the notation of the clump spectra given above ) .
such a distribution for hierarchical structure was also recognized by larson ( 1978 ) and fleck ( 1996 ) .
the result gives a steeper instantaneous mass function than emission line or extinction surveys , but there is no conservation of mass as in these standard surveys ( i.e. , the sum of the masses of all possible structures is larger than the mass of the whole cloud , because nearly all of the structures are contained in other structures , and so are multiply counted ) .
this way of viewing cloud structure seems to be closer to what a real star - forming process `` sees '' before it actually begins the sequence of events that makes a star .
the density dependent rate of star formation then steepens the mass function from @xmath54 to @xmath56 .
this result is then identified with the salpeter imf .    for the purposes of understanding the imf
, there are probably better ways to measure cloud structures than with clump - finding algorithms .
the structure on the edges of clouds has been analyzed in terms of fractals , rather than clumps , for many years ( beech 1987 ; bazell & dsert 1988 ; scalo 1990 ; dickman , horvath , & margulis 1990 ; falgarone , phillips , & walker 1991 ; zimmermann & stutzki 1992 , 1993 ; henriksen , 1991 ; hetem & lepine 1993 ; vogelaar & wakker 1994 ; pfenniger & combes 1994 ) .
for the structures inside clouds , stutzki et al .
( 1998 ) showed that fourier transform power spectra of emission line intensity scans across molecular clouds gave power laws , and concluded that the internal structure was scale - invariant , as in a fractal .
they reproduced this structure with a random fractal - generation model .
other models that made cloud fractals were in hetem & lepine ( 1993 ) and elmegreen ( 1997b ) .
the concept that the cloud mass distribution function is the result of fractal structure , presumably generated by turbulence , began with papers by fleck ( 1996 ) , elmegreen & falgarone ( 1996 ) , and stutzki et al .
this seems to be a more reasonable explanation than the older collisional - build - up models of clump structure , particularly since off - center supersonic collisions between clumps should not be sticky ( scalo & pumphrey 1982 ; kimura & tosa 1996 ; fujimoto , & kumai 1997 ) .
the probability of selecting structures for star formation in hierarchical clouds seems to show up more directly in the distribution of masses for open clusters .
when a bound cluster forms , a high fraction of the gas mass has to go into stars ( see verschueren 1990 and references therein ) , so the cluster mass distribution should reflect the mass distribution of cloud structures pretty well . in this case , the result is obviously independent of telescope resolution or cloud - clump recognition bias , because clusters are seen optically , even at great distances .
also , the range of masses for clusters can be reasonably large , exceeding a factor of 100 , so a power law can be measured fairly well if it exists .
cluster selection suffers from various other biases , however , such as extinction , age limitations , and a loss of low - mass clusters with increasing distance . nevertheless , there are some studies that get around these biases , and they confirm the expected result . in two samples of nearby clusters , each with calibrated masses , battinelli et al .
( 1994 ) found power law mass functions with slopes of @xmath57 and @xmath58 .
these clusters are close enough to the sun to be relatively free of distance and extinction effects . also , in the lmc , where @xmath59 clusters are catalogued ( bica et al .
1996 ) , calibrated photometrically for age ( girardi et al 1995 ) , and all at about the same distance , elmegreen & efremov ( 1997 ) found @xmath60 mass functions for discrete age groups ( e.g. , @xmath61 yrs ; dividing the clusters into age groups is important when there are no direct measures of cluster mass , because the conversion from luminosity to mass depends on age ) .    the imf should be steeper than the cluster mass function because of the density - dependent formation rate of individual stars , and the competition for mass in the imf .
these effects are not important for clusters .
open clusters take a high fraction of all the gas mass originally available to them ( or else they would not have ended up bound ) , and the rate at which they do this is not important for their final masses . in the case of stars , new objects forming at one level
steal mass away from the objects that form later at a higher level .
but clusters do not do this : those which form inside other clusters in the general hierarchy can just merge together to make a single cluster of higher mass , appropriate to the higher level in the hierarchy . and if they do not merge , then they stay with their original masses , which are appropriate for the levels in which they form .
other imf models based on fractal or hierarchical cloud structure were developed earlier by henriksen ( 1986 , 1991 ) and larson ( 1992 ) .
henriksen used the size , @xmath62 , distribution function of structures in a fractal of dimension @xmath63 , given as @xmath64 by mandelbrot ( 1983 ) , and assumed that the density @xmath65 varied with size too , not as in a fractal ( which would be @xmath66 ) but as actually observed in self - gravitating clouds , namely @xmath67 .
this gave a mass function for cloud structure that agreed well with observations if @xmath68 .
note that if henriksen used the density dependence for a fractal , in a self - consistent manner with the size distribution , then the resulting mass function would have been independent of the fractal dimension , @xmath60 , as discussed above for hierarchical clouds in general .
larson ( 1992 ) started again with the size distribution for a fractal , and assumed that final stellar mass is directly proportional to the linear size of the cloud structure , suggesting that such structures were filamentary anyway .
then the slope of the imf in linear intervals became @xmath69 for @xmath70 , as implied by the fractal dimension of cloud perimeters , namely @xmath71 , added to 1 to account for the higher dimension of the surface .
this imf , with a slope of @xmath72 , is significantly steeper than the salpeter imf , with a slope of @xmath73 , but larson suggested that maybe stars form in subparts of clouds where @xmath63 is smaller than 2.3 , or perhaps larger structures have larger temperatures , which would break the assumed linear relationship between star mass and scale size .
these models differ from the elmegreen ( 1997a , 1998 ) model , even though both employ `` fractal '' cloud structure , because the latter uses general cloud shapes , not necessarily filaments or sheets , random sampling from all levels in the hierarchy of structures ( giving @xmath60 from sampling alone ) , a density - dependent rate that steepens the imf by preferred sampling at lower mass ( i.e. , steepening the result to @xmath74 ) , and a competition for mass , which steepens the imf further , to @xmath75 , as in the salpeter imf . with these different assumptions ,
the imf in the elmegreen model hardly depends on the fractal dimension at all ; it enters into the exponent of the imf power law as approximately @xmath76 , which is a very small increment for @xmath63 in the likely range of 2 to 3 .      understanding the lower mass limit to star formation is an important part of any imf theory .
various assumptions about this have already been mentioned in the review of theories presented above .
three views on this subject seem to pervade most of them . in the wind - limited accretion models
, the assumption is made that accretion onto a dense core continues until a protostellar wind develops , which presumably follows the onset of deuterium burning ( shu et al .
thus low mass objects do not stay that way very long , they just keep accreting until they become active protostars . in most of the older fragmentation - coalescence models , the minimum mass is the mass of a gravitationally bound fragment that is optically thick . in the various scale - free models ,
based on gravitationally - driven fragmentation or turbulence , the minimum mass is usually the thermal jeans mass , particularly in the larson ( 1992 ) and elmegreen ( 1997a , 1998 ) models .
this assumption is made because , without further fragmentation during the collapse itself , the thermal jeans mass is the smallest structure that can ever become self - gravitating , even in the presence of turbulence and magnetic fields .
the first and the third of these mass limits equals about the observed limit , @xmath4 m@xmath0 . in the first case , it is self - defining : objects accrete until they become star - like , and then they stop accreting , so naturally the lowest possible mass is the mass of a star .
this concept has the powerful advantage over the others that the basic star mass , and perhaps even the imf itself , can be the same everywhere , independent of cloud conditions . in the third case
, the minimum theoretical mass just happens to equal the observed minimum in standard cloud conditions .
an expression for this mass is the bonner - ebert critical mass , @xmath77 here we have normalized the result to standard cloud conditions , taking the pressure inside a cloud equal to the typical self - gravitating core pressure , not the boundary pressure ( which is sometimes lower by a factor of @xmath78 ) .
our use of this form for the thermal jeans mass is motivated by the near constancy of @xmath79 and total @xmath80 in molecular clouds , according to the scaling relations ( larson 1981 ) .
this is much more sensible than writing the thermal jeans mass in terms of @xmath79 and density , for example , when the density varies by several orders of magnitude inside a cloud .
our use of the total pressure for @xmath80 also allows for a decrease of jeans mass in compressive flows , as in hunter & fleck ( 1982 ) .
nevertheless , another theory of the imf in which all of the stellar masses have the local thermal jeans mass has been advanced by padoan , nordlund , & jones ( 1997 ) , with variations in mass entirely the result of variations in the density in a turbulent fluid . in this model ,
the high mass stars form in the low density regions , because the jeans mass is high there for a given temperature .
the second minimum fragment mass discussed above , from the opacity limit of a self - gravitating condensation , is typically much smaller than the minimum star mass , so coalescence has to bring the star mass up to the observed value in these models .
larson has looked for signatures of the thermal jeans mass in several ways , predicting that it should vary with cloud temperature ( larson 1985 , 1992 ) . in larson ( 1992 )
, he suggested that such variations were observed by myers ( 1991 ) in different molecular clouds .
larson ( 1995 ) also suggested that the length scale at a break in the separation distribution for t tauri stars was the jeans length ( see also simon 1997 ; but see bate , clarke , & mccaughrean 1998 and nakajima et al .
1998 for other interpretations . )
elmegreen ( 1997a , 1998 ) considered @xmath81 from a different point of view , stating that most star - forming regions in normal galaxies actually do have nearly constant values of @xmath82 , with small variations in this quantity possibly leading to the observed variations in cluster imfs at low mass ( see section [ sect : intro ] above ) .
this constancy is a result of the numerator in @xmath81 being approximately proportional to the cooling rate per unit mass of molecular gas ( neufeld , lepp , & melnick 1995 ) , and the denominator being approximately proportional to the local surface density of stars and gas in the galaxy , which determine the heating rate per unit mass from cosmic rays and background starlight . as long as cooling roughly balances background heating
, this ratio will be about constant .
other , more local variations in @xmath79 and @xmath80 tend to cancel out when combined as the quantity @xmath82 .
for example , when @xmath80 goes up in a triggered region next to a previous generation of star - formation , @xmath79 generally goes up too because of the enhanced radiation field .
conversely , when @xmath80 is low , as in low - mass clouds like taurus , then @xmath79 is low as well .
elmegreen ( 1997a , 1998 ) also points out that the lower mass limit for stars actually _ does _ change , perhaps by a factor of 10 in extreme starburst regions ( see sect .
[ sect : intro ] above ) where the average @xmath79 can be much larger than in local molecular clouds , e.g. , 100 k instead of 10 k ( cf .
aalto et al .
1995 ) , and the pressure much larger too ( e.g. , @xmath83 ) . but
this variation should be rare in normal galaxies because of feedback processes in the interstellar medium that tend to connect young stellar activity , interstellar pressure , and gas temperature .
if we consider @xmath79 to be a measure of the rate of cooling per unit cloud mass , then @xmath79 in a star - forming region will be high when there is an unusually large embedded or nearby luminosity of stars per unit gas mass .
this occurs when the efficiency for star formation is high in a high mass cluster ( the cluster has to be fairly high mass to contain the luminous stars that dominate the total radiation field ) .
a high efficiency in a low mass cluster would not necessarily increase @xmath81 because low mass clusters are not very luminous per unit mass .
a low efficiency in a high mass cluster , leading to an ob association , might not change @xmath81 much either , because the radiation luminosity per unit gas mass is not particularly large there either .
thus the formation of globular clusters in regions with only moderately large pressures , not the enormous pressures typical of early galactic halos , might be elevate @xmath81 , but this conclusion is very uncertain .
figure [ fig : mj ] shows two numerical experiments that made model imfs from random selections in a hierarchically structured gas cloud , using a rate of selection proportional to the square root of the gas density , as discussed above ( elmegreen 1997a ) .
the models differ only in the value of @xmath81 , which enters into the probability for failure to form a star , taken to be @xmath84 in these cases ( @xmath81 in the figure is written in terms of the computer - program mass , which is scaled to a physical mass after multiplication by @xmath85 .
the imfs are identical except for the larger value of the lower mass limit in the case with high @xmath81 .
the two main contenders for the lower mass limit of star formation , wind - regulation and the thermal jeans mass , have obvious differences in predicted behaviors of the shape of the imf at low mass .
if this shape varies a lot , with some clusters having clearly higher lower mass limits than other clusters , then wind self - regulation might be ruled out . moreover ,
if cluster imfs in regions with large turnover masses also have larger pre - stellar @xmath81 , then the thermal jeans mass would seem to be important .
curiously , both of these models , the wind - regulated imf , and the @xmath81fractal cloud model , have their strong points where the other is weakest .
the wind - regulated imf model seems to have difficulty explaining the nearly universal slope of the imf without some sensitivity to cloud parameters , yet it seems certain that wind self - regulation must play some role in determining the mass that gets onto a star , at least near the lower mass limit for star formation .
similarly , models that rely on @xmath81 have an uncomfortable susceptibility to changes in cloud properties near the lower mass limit , but get the salpeter slope above this limit with essentially no free parameters .
perhaps a combination of the two models would be better .
there are several ways this could occur .
first , the fractal model supposes that a star s mass can be identified with the mass of some structure in a turbulent cloud , but this is obviously too simple  there is a lot that can happen during star formation that will vary the fraction of the cloud mass that gets into a star .
if this fraction is randomly distributed , then the same imf slope and scatter about that slope result , as shown in connection with figure [ fig : scatter ] below , which is discussed in more detail in elmegreen ( 1998 ) .
if it is not randomly distributed , then the power - law slope could change .
but the entire power - law slope is not likely to come from wind - regulation , that would give the slope too great a sensitivity to cloud properties and the randomness of wind - clearing .
this means that , perhaps to within a factor of 2 to 5 , star mass should really be identified with the mass of the clump in which it forms , i.e. , that bigger stars really do form in bigger clumps ( casoli et al .
1986 ; myers , ladd , & fuller 1991 ; myers & fuller 1993 ; pound & blitz 1995 ; motte , andr , & neri 1998 ; not to be confused with the statistical sampling statement discussed in section [ sect : intro ] , that bigger stars form in bigger _ whole _ clouds ) .
this factor of 2 to 5 may contain all the detailed physics of wind clearing .
protostellar winds should also play some role in regulating the smallest stars that form . according to the models by shu and collaborators , as reviewed in shu et al .
( 1987 ) , the onset of the wind is triggered by deuterium burning in a young , pre - stellar object .
if the object is too small , there will be no wind , and presumably accretion will continue until deuterium burning begins .
this concept combines well with the arguments based on a thermal jean mass .
for example , if @xmath81 is much smaller than the deuterium burning limit , then a large number of small brown dwarfs or jupiters should form , and these could either leave the cloud in that state , forever drifting as brown dwarfs , or they could accrete more mass over time and end up finally able to start a wind . in the latter case ,
the effective lower mass limit to the imf would be the deuterium burning mass , not @xmath81 , and the larger masses would follow the fractal mass distribution as before , giving the salpeter function .
on the other hand , if there are a lot of free brown dwarfs , which does not appear to be the case in the solar vicinity ( cf .
[ sect : intro ] ) , and if there are also a lot of regions where @xmath81 is low , which does not appear to be the case locally either , then a model with accretion up to a wind - limited mass would not seem to work .    if @xmath81 is much larger than the deuterium burning limit , then large , self - gravitating clumps will form in the cloud , but the rapid onset of a wind in these clouds might limit the gas mass that actually gets into the star to only a small fraction of what is available . then all of the excess mass between the deuterium limit and m@xmath86 would get dispersed back into the cloud for recycling into other stars , and the lower mass limit would be the deuterium limit . in this second case ,
the mass fraction that goes into each star will be small , but if it is about the same fraction for all clumps , then the higher mass stars will still map out the fractal properties of the cloud and produce the salpeter imf . again
the details of how the imf gets established in a cloud depend strongly on unknown processes during the final collapse and dispersal phase of star formation , but the basic shape of the imf could be relatively independent of these details .
there have been some fundamental changes lately in how we view molecular clouds and star formation , and these changes inevitably affect the various models for the imf .    for example , we are beginning to think that star formation is relatively rapid in cloud cores , occurring in only one or two crossing times regardless of scale ( elmegreen & efremov 1996 ; efremov & elmegreen 1998 ) .
this means that we can not generally wait for magnetic diffusion to occur as a cloud core slowly accretes across field lines
. diffusion typically takes about 10 free fall times for cosmic ray ionization ( shu et al .
1987 ) , longer for ionization by starlight ( myers & khersonsky 1995 ) , and longer still if the gas is clumpy ( elmegreen & combes 1992 ) .
indeed , nakano ( 1998 ) suggested that all star - forming cloud cores are magnetically _ super_critical ; i.e. , they would collapse dynamically across the field lines if turbulent motions were not present .
there are two important implications of this change in thinking , if correct .
first , there would no longer be a strong motivation for imf theories that consider bimodel star formation in the sense that low mass stars come from subcritical cores , and high mass stars come from magnetically supercritical cores ( lizano & shu 1988 ) .
the imf is so remarkably uniform anyway , it does not seem possible to have widely different modes of star formation at high and low mass .
any change from low mass dominance to high mass dominance in the imf could more easily be accommodated by an upward shift in the thermal jeans mass ( cf .
section [ sect : theorylm ] ) .    secondly , this time scale for star formation may be too short to allow multiple interactions between protostellar clumps .
this would then rule out a broad class of coalescence and slow accretion models .
such coalescence is unlikely anyway for clumps that move supersonically relative to each other ; they will fragment or disperse upon collision , rather than stick together .
clump collisions may trigger star formation when they occur ( bhattal et al .
1998 ) , but multiple collisions are probably not the source of the cloud or star mass distributions .    another observation that suggests the same thing is the extremely high stellar densities in young globular clusters , which form today in starburst regions .
as suggested in section [ sect : intro ] , these globulars look like they will evolve into objects similar to our milky way s halo globulars , and the old milky way globulars have an apparently normal imf .
if coalescence is important in local `` normal '' embedded clusters , it would seem to be devastating in globular clusters , where the stellar density can be @xmath87 times higher .
conversely , if these globulars have normal imfs , then protostellar coalescence must be unimportant nearly everywhere .
a similar constraint would come from the imf in small regions , where only a few stars form .
there , protostar interactions are not likely to be important either .
a second major shift in our view of molecular clouds is that most astronomers now think that pre - stellar structure comes from turbulence , not gravitationally - driven fragmentation .
this is because the same character of structure is observed in both self - gravitating and non - self - gravitating clouds , and both types of clouds have correlated motions reminiscent of laboratory turbulence .
a related point is that the cloud cores in which clusters form are not obviously collapsing as a whole : there are no inverse p - cygni profiles indicative of collapse motions for whole cores , as there are in some individual protostars .
the stars are apparently forming inside more - or - less stable cloud cores .
once again , the general collapse models seem untenable
.    the presence of supersonic turbulence also seems to rule out models involving thermal instabilities .
turbulent motions in cluster - forming cores are generally much larger than thermal , and so the dominant forces that structure the gas are the turbulent forces , not the thermal .
this means that turbulence causes the observed structures in pre - stellar clouds , not thermal instabilities .
thermal processes are probably important on the smallest scales , but these are at and below @xmath81 .
turbulence in pre - stellar clouds also appears to generate structure that is much lower in mass than even the smallest stars , perhaps as low as @xmath88 m@xmath0 ( heithausen et al . 1998 ; kramer et al .
this means that star formation occurs in the middle range of all cloud structures , not at the bottom end .
it also means that , while turbulence may generate a cascade of structures down to smaller and smaller scales , the end result of this cascade is not the formation of a star .
turbulence makes cloud structures independent of star formation .
self - gravity , sonic wave generation , magnetic reconnection , and many other aspects of interstellar cloud dynamics are likely to be important too by the time star formation begins .
the remarkable near - uniformity of the imf in regions with a wide range of stellar densities , cluster masses , metallicities , ages , and galactic types strongly suggests there is a _ unifying process in star formation _ that makes the relative probability of forming various masses somewhat independent of physical parameters .
for example , the imf seems independent of how star formation begins , i.e. , whether the region is triggered by some external pressure or forms quiescently ( e.g. see parker et al .
1992 and the contrary opinion by oey & massey 1995 ) , independent of the general cloud shape ( shell , layer , filament , etc . ) , magnetic field strength , ionization level , presence of other stars , and so on .
this is remarkable because these physical variables do control the _ rate _ of star formation in some regions , and they certainly control _ where and when _ star formation begins .
they just do not appear to influence the imf .
the near - uniformity of the _ slope _ of the power law also implies several other things .
( 1 ) there is a _ universal aspect of cloud structure and evolution that produces the power - law portion of the imf independently of the temperature , pressure , and metallicity .
_ this means that the power - law probably does not depend on the protostellar collapse process , e.g. , on the fragmentation or collapse of isothermal clumps , the accretion of gas in thermal equilibrium , or the thermal jeans mass . for all of these processes , thermal temperature and magnetic diffusion are probably important at some stage ( see reviews in nakano 1984 ; mouschovias 1991 ; shu , adams , & lizano 1987 ) .
( 2 ) the imf power - law is also likely to be a _
highly reduced average over many physical processes that either all give about the same result , or are combined so finely in every region that the same proportion of each process is always present . _
the variations in slope of the imf from region to region ( scalo 1998 ) seem at first to suggest something different , that the imf is not in fact uniform but depends sensitively on physical conditions . but are these variations statistically significant ?
rarely do studies of cluster imfs contain more than several hundred stars .
figure [ fig : fluctuations ] shows sample imfs from the model of random selection in a hierarchical cloud ( elmegreen 1998 ) .
the panel in the upper left shows imf slopes for 100 different models , each with a different number of stars and therefore different upper mass limit . in all cases ,
the range of stars chosen for the power - law slope determination is such that there are 200 stars in the fit , starting from the histogram bin that is twenty filled bins ( a factor of 3 in mass ) away from the most massive star , and going to lower mass bins .
this starting point for the fit ensures that there are enough stars to get an accurate slope , and it coincides with an astronomer s decision to avoid fitting the imf to the scarce few highest - mass stars in the cluster .
we pick a different number of total stars for each imf , and use only the highest mass parts of each model , because this is what an astronomer would do as well .
generally an observer catalogues only the brightest stars in a cluster and misses the low mass members . if a cluster has an imf made only from low mass stars , then that cluster is generally so small that there are no high mass stars at all .
evidently the imf slope varies a lot around the salpeter value , shown in the figure by the dashed line at @xmath22 .
this average slope is nearly independent of physical cloud properties in this model ( see the discussion in section [ sect : theory ] ) .
the other three panels in figure [ fig : fluctuations ] show particular imfs and the mass ranges and fits used for the `` x '' marks in the top left panel .
figure [ fig : scatter ] shows the rms deviations in the imf slopes around the mean values for power - law fits that include different numbers of stars in the same 100 imfs ( from elmegreen 1998 ) .
each plotted value , connected by a smooth line , represents the rms deviation for the 100 models , calculated in the mass range from 1 to 10 m@xmath0 ( the characteristic physical mass in the model comes from @xmath81 , which is taken to be 0.35m@xmath0 from equation [ eq : mj ] ) .
the rms deviation around the salpeter slope decreases systematically as more and more stars are included in the fit
. the deviation can be as high as @xmath6 when only 80 stars are used , decreasing to @xmath4 when 1000 stars are used .
this implies that an imf calculated with , say , 150 stars , and having a slope of @xmath89 instead of the salpeter slope of @xmath14 , is in fact statistically consistent with this salpeter slope in a `` universal '' imf .
the fluctuations found by scalo ( 1998 ) in his review correspond to @xmath6 around the value @xmath90 .
these fluctuations could be statistical , considering the number of stars included in typical imf surveys and the likely presence of measurement and star selection errors in the real data .
the dashed line in figure [ fig : scatter ] is for another 100 models , independent of the first , and evaluated with a large ( @xmath91 ) and random fluctuation in the ratio of the star mass to the clump mass .
such variations do not affect the average imf slope or the statistical fluctuations around it .
the imf has been observed directly and indirectly for many years , with many different results , but there seems to be some convergence now in the basic form of the imf , and also some tantalizing indications that this form changes a little when the physical conditions for star formation change a lot .
the main observational results were summarized in section [ sect : sum1 ] .
the theory of the imf seems to be all over the map , even in recent years .
this variation reflects our ignorance in the processes that determine a star s mass and in the processes that cause the complex spatial structures seen in interstellar clouds .
the common aspects in many of these models , including clump and protostellar interactions , wind erosion of protostellar gas concentrations , minimum masses from the lack of self - gravity , and hierarchical cloud structure , are all likely to play some role in generating the imf .
the worry is that the imf is such a highly reduced average over many physical processes that each process does not a have clear signature in the final result .    in view of the observations
, we should perhaps consider `` acceptable '' those models that do not have much sensitivity to either the local aspects of cloud structure , or the large scale aspects of time and place in the universe .
this drives the recent appeal towards models based on common turbulence in one way or another , considering that turbulent structures are likely to be robust .
unfortunately , we do not understand compressible mhd turbulence much either .
there are several key observations that would help .
one is the determination of star formation time scales , always in units of the cloud or clump crossing time .
there are three important timescales : the rise time of the star formation rate in a cloud , the duration of star formation , and the decay time scale .
if the rise time is less than a crossing time , and the duration only one or two crossing times , then imf models based on multiple clump interactions would seem to be ruled out .
there is also a need to know the mass dependence of the mass fraction of clump gas that goes into a star .
this might be determined from the ratio of luminosity to gas mass for type 0 sources , plotted as a function of gas mass for the dense cores where the stars are actually forming .
theory should tell us what gas structures and evolutionary timescales are expected from self - gravitating mhd turbulence in the supersonic - subalfvnic regime of molecular clouds .
if turbulence only makes moving waves and wavepackets , and star formation actually occurs in such regions , then star formation has to be relatively quick in terms of the local crossing time .
it may be that in the absence of strong self - gravity , turbulent clumps are transient and amorphous , but when self - gravity becomes important , the clumps get some integrity and persist for a relatively long time .