there is a great variety of astronomical objects in the universe  by increasing average size or mass : asteroids , planets , brown dwarfs , stars and stellar remnants , binary and multiple stars , star clusters , molecular clouds , galaxies , groups and clusters of galaxies .
each of these classes of objects follows a certain ( but generally not well known ) distribution function in size , luminosity , or mass .
the most fundamental characteristic property of all these distribution functions is that smaller or less massive things of a given kind are more abundant than larger or more massive ones .
presumably the same would hold not only within a certain class of objects , but also if one were to glue together the individual distribution functions to build a `` universal '' distribution function , as there would be , _ per unit volume of space _ , more asteroids than planets , more planets than stars , more stars than galaxies , etc .    in this paper
we are doing just that : putting together all known distribution functions for different object classes to see whether there is such a thing like a `` universal distribution function '' .
 `` universality '' has a twofold meaning here : we can ask ( 1 ) whether the all - object distribution function is sufficiently _ continuous _ to render the notion of universality useful in the first place , and ( 2 ) if yes , whether the _ slope _ of the individual distribution functions and that of the grand total is , to a certain degree , universally the same . for the present study we choose _ mass _ as the independent parameter , as mass
is probably the most fundamental  albeit difficult to determine  property of an astronomical object ( luminosity is inappropriate because there are dark objects ) .
we show that by proper normalization we can indeed establish a continuous mass function of the universe from @xmath4 ( 100 meter - size asteroids ) to at least @xmath5 ( the most massive clusters of galaxies ) , and that there is surprising similarity in the slope of the mass function over a large range of mass .
the question of a universal distribution function of astronomical objects has apparently been of interest only to fritz zwicky . in his 1942 physical review paper
entitled `` on the large scale distribution of matter in the universe '' , zwicky argued on purely theoretical grounds that the luminosity function ( lf ) of galaxies had to be essentially exponential , i.e.  ever fainter galaxies being increasingly more numerous , just as observed today ( see below )  but in contradiction with the data available then , which showed a more gaussian - type lf for galaxies .
zwicky s vision was that galaxies constitute an ensemble of `` particles '' in a stationary state of statistical equilibrium .
frequent close encounters would dissolve larger aggregates of galaxies on the one hand , and accumulate new ones on the other hand , resulting in a boltzman - type energy distribution for galaxies and groups and clusters of galaxies _ as a whole_. today , we know that the temperature and density of the `` galaxy gas '' is much too small , and thus the time - scale of galaxy interactions much too long , for a statistical equilibrium to be established or maintained .
but zwicky was perfectly right not only in his prediction of an exponential lf ( which he himself set out to prove , see zwicky 1957 , p.  220ff ) , but also in his view that galaxian systems , from dwarf galaxies up to the richest clusters of galaxies , follow one , continuous distribution function ( which was observationally established for the first time by bahcall 1979 , see also below ) .
this unity of objects on the scale of galaxies and systems of galaxies is now understood in terms of a general , cosmological `` bottom - up scenario '' centred around the hierarchical , gravitational clustering of dark halos in an expanding universe having formed from primordial density fluctuations , as , e.g. , described by the `` press - schechter formalism '' ( press & schechter 1974 ) .
while zwicky considered stars in his 1942 paper as well , he was well aware that the mass distribution of stars must be shaped by totally different physical processes .
put in simple terms , the relation between stars and ( small ) galaxies is plausibly one of gas fragmentation ( from large to small ) , rather than clustering ( from small to large ) .
so why should there be a universal mass function , if there can hardly be universal ( global ) physics behind the frequency distribution of such things like asteroids , stars , and galaxy clusters ?
this must be the reason why the question of a universal mass distribution ( beyond the scale of galaxies and systems of galaxies ) has never been taken up since zwicky ; it simply does not seem to make much sense from a physical point of view .
conversely , the  surprising  degree of universality in the mass distribution of astronomical objects found in the present paper does not entitle us to conclude that there indeed is overarching physics at work , though this can not be excluded either . while we will , in the discussion section , briefly indulge in speculations on what it all could mean , we emphasize that our primary goal is not to find a sort of cosmological principle that governs the mass distribution function of the universe , but to provide a valuable _ piece of cosmography_. a good knowledge of the frequency distribution of things in the universe is an end in itself .
there are conceptual difficulties with the definition of an astronomical object .
what exactly is an astronomical object , i.e.  what kind of objects should we include in such a study ?
the most simple requirement would be that such an object is gravitationally bound . however , sub - kilometer asteroids ( let alone dust particles , if we were to push the distribution function to the extremes ) are hold together by non - gravitational forces . also , on large scales , where gravity does dominate , we wish , for the sake of continuity , to include molecular clouds which are in general unbound but are _ identifiable entities _ none the less .
each class of objects has its own specific technical difficulties in the establishment of a distribution function .
there is a rich history of attempts , e.g. , to determine the stellar lf ( and thus the initial mass function , imf ) , or the lf and mass - to - light ( m / l ) ratio of galaxies .
the present paper is not meant to be a meta - study where all these attempts for the different object classes are critically assessed .
rather , we will take `` reasonable '' distribution functions from the recent literature , with no consideration of the uncertainties , and put them together to construct a _ first - guess universal mass function_. these distribution functions will partially overlap in mass and partially be separated by gaps in mass regions where we lack objects or data .
another principal difficulty of such a study is the normalization .
everything has to be normalized to the same volume of space , but our objects are sampled in vastly different volumes : while clusters of galaxies are seen over a significant part of the observable universe , and the galaxian lf , though sampled in a fairly local volume , can still be regarded as universal , the stellar lf is simply that of the solar neighbourhood , and the mass function of planets and small bodies is entirely drawn from our solar system ( never mind the exponentially growing number of  nearby  exoplanets ) .
so the smaller the scale the larger and more uncertain the extrapolation to a universal mean .
this kind of uncertainty will not soon go away .
rather , we have to rely , as all of cosmology does , on the cosmological principle , i.e. that we are living in a typical , average environment  on all scales .
all our mass functions will in the end be normalized to a volume of 1 mpc@xmath6 and given in units of solar masses ( [ @xmath7 = { \rm mpc}^{-3 } { \rm m_{\sun}}^{-1}$ ] ) .
the paper is organized as follows . in section 2
we present the individual mass functions for the different object classes , starting with galaxies and going up first to supergalactic , then going down to subgalactic structures .
these functions are put together to a universal mass function in section 3 , where the mass distribution of dark matter halos from theory is added for comparison .
finally , in section 4 we discuss possible meanings of the universality found and give our conclusions .
a hubble constant of 70 km s@xmath8 is used throughout .
while the galaxian lf is fairly well known , except at the very faintest end , the mass function ( hereafter mf ) of galaxies is very hard to determine , as galaxies are dominated by non - baryonic dark matter ( dm ) of a principally unknown nature .
the galaxian mf is reconstructed from individual or statistical mass - to - light - ratio determinations ( or rather estimates ! ) that have to be convolved with the luminosity distribution which is usually taken to be a schechter ( 1976 ) function of the form + @xmath9 vale & ostriker ( 2004 ) give the following mass - to - light function for galaxies : + @xmath10^{1/k}}\ ] ] with the best - fitting parameters @xmath11 , @xmath12 , @xmath13 = 4 , @xmath14 = 0.57 , @xmath15 = 3.72 , and @xmath16 = 0.23 . at the low - mass end
this gives @xmath17 , at the high - mass end @xmath18 .
a total galaxy mf is then obtained through @xmath19 by inserting ( 2 ) and the schechter function ( 1 ) with parameters taken from the 2df galaxy lf ( norberg et al .
2002 ) , as used also by vale & ostriker ( 2004 ) : @xmath20 , @xmath21 ( in the @xmath22 band ) , and @xmath23 .
however , we use this vale & ostriker galaxy mf only for galaxy masses smaller than @xmath24 , as the high - mass end ( with its shallow slope , see above ) would produce too many cd - type galaxies ( as vale and ostriker note , the cd halos that went into their mass - to - light function are rather hosting the whole cluster than the cd alone ) .
a more appropriate mf at the high - mass end may be provided by van den bosch et al .
( 2005 ) who work with the statistical distribution of satellite galaxies .
their mass - to - light function is given by @xmath25}\ ] ] with the best - fitting parameters @xmath26 = 115 , log @xmath27 = 10.76 , log @xmath28 = 12.15 , @xmath29 = 2.92 and @xmath30 = 0.29 ( for a @xmath31cdm clustering normalization of @xmath32 = 0.9 ) .
this is somewhat steeper at the high - mass end ( @xmath33 versus @xmath34 ) , giving a total mf that is falling down more rapidly now with increasing mass , while it is roughly compatible with vale & ostriker ( 2004 ) at the low - mass end .
we therefore use the resulting van den bosch et al .
mf for galaxies with masses larger than @xmath24 , keeping however the vale & ostriker mf for smaller masses , with continuity at @xmath35 .
this _ total _
( dm - dominated ) galaxy mf , as a combination of the two functions for larger and smaller masses as just described , is shown in fig .  1 .
the universal normalization ( per mpc@xmath6 ) is of course provided by the adopted normalized lf .    for comparison
we also show the baryonic mf of galaxies in fig .  1 ( stars , stellar remnants , atomic and molecular gas ) , taken from read & trentham ( 2005 ) .
it is based on a host of nir broad - band as well as radio observations and is given in the form of a schechter mass function @xmath36 with parameters @xmath37 , @xmath38 , and @xmath23 .
the baryonic galaxy mass is indeed dominated by stars ; the schechter @xmath39 is the same ( see read & trentham 2005 for a detailed account of the contribution by the gas ) .
we show the galaxy mfs down to a mass of @xmath40 ; this is of course a bold extrapolation ; the faint end of the lf , let alone the mf of galaxies , is not well constrained .
the mass of groups and clusters of galaxies is dominated by dm as well and has to be estimated either from the kinematics of member galaxies or ( for clusters only ) from the x - ray intracluster gas properties or the gravitational lensing of background galaxies .
there are several studies of the group and cluster mf in the literature ( e.g.  bahcall & cen 1993 , biviano et al .
1993 , girardi et al .
1998 ) . for our purposes
we take the mf of girardi et al .
( 1998 ) which is based on optical virial mass estimates of a complete sample of 152 nearby clusters .
it was derived only for masses larger than a few @xmath41 , but we extrapolate it down to @xmath42 , because it is compatible with the mf of bahcall & cen ( 1993 ) which does cover also the low - mass end by including small groups . the adopted mf for groups and clusters , shown in fig .
2 , is a schechter - type mf of the form ( 5 ) , however fixing @xmath43 ( flat end ) , with @xmath44 , and @xmath45 .
next we consider subgalactic structures : star clusters and molecular clouds which can be as massive as extreme dwarf galaxies but in contrast to these do not contain dark matter .
the observed lf of globular clusters is typically bell - shaped ( gaussian ) .
however , the decline towards fainter luminosities ( masses ) is probably an effect of dynamical evolution ( dissolution of smaller clusters in the galactic tidal field ) and may not reflect the initial distribution [ for a different view , however , see parmentier & gilmore 2007 ] .
as suggested by surdin ( 1979 ) , racine ( 1980 ) , and harris & pudritz ( 1994 ) , the initial mass distribution of globular clusters can be modelled by a power law of the form @xmath46 or in our mf notation : @xmath47 such a distribution is indeed observed for young clusters in merging galaxies ( e.g. whitmore & schweizer 1995 ) , and the high - mass tails of the old globular cluster systems of nearby galaxies are still following this form very closely .
for globulars with masses larger than @xmath48 , harris & pudritz ( 1994 ) find a mean best - fitting @xmath49 , which we adopt for our globular cluster mf in the range @xmath50 .
for lower masses we consider only open clusters ( see below ) .
how can this mf now be normalized to a universal density ?
this is done by the `` specific frequency '' @xmath51 which is the number of globular clusters per `` piece '' of mother galaxy of absolute @xmath52 magnitude @xmath53 .
@xmath54 is the estimated total number of globular clusters .
@xmath55 depends on the morphological type of the mother galaxy , but a reasonable mean value is @xmath56 ( see harris 2001 ) , which we adopt . with a universal @xmath52 luminosity density ,
@xmath57 , we have for the mean number of globular clusters per mpc@xmath6 : @xmath58}{8.55 \times 10 ^ 7}\,\,.\ ] ] norberg et al .
( 2002 ) give a @xmath22 band @xmath59 . as the mean @xmath60 colour of galaxies
will not be far from the solar value of 0.65 , we assume the same @xmath57 value in @xmath52 .
the resulting normalization constant from ( 9 ) is then @xmath61 .
as shown , e.g. , by elmegreen & efremov ( 1997 ) from their study of hundreds of open clusters in the lmc , the lf of open clusters systematically depends on age . using the same power law for the lf , @xmath62 @xmath63 is smaller ( i.e.  the lf more flattened ) the younger the cluster .
very young clusters ( age less than @xmath64 years ) have essentially a flat distribution ( @xmath65 1 ) ; somewhat older clusters up to an age of @xmath66 years show @xmath65 1.5 , and for clusters older than @xmath66 years we finally have @xmath65 2 .
the flat lfs found with young clusters is obviously due to the presence ( and dominance , in terms of light ) of short - lived ob stars .
clearly , then , for the _ mass _ distribution ( 7 ) of open clusters , @xmath63 = 2 is a reasonable choice .
a normalization ( fixing constant @xmath67 in ( 7 ) for open clusters ) is achieved by comparing the populations of open and globular clusters in the milky way galaxy . for open clusters we rely on vandenbergh & lafontaine ( 1984 ) who , aside from finding @xmath65 2 as well , give a mean surface density for open clusters in the solar neighbourhood of @xmath68
[ see piskunov et al .
( 2006 ) and references therein for more recent open cluster population studies ] .
this is for clusters of absolute @xmath52 magnitude between @xmath3 and @xmath69 , roughly corresponding to a mass range of @xmath70 to @xmath71 solar masses .
if we assume a more or less constant cluster density in the galactic disk with 10 kpc radius , giving a total disk surface of @xmath72 , we expect a total number of galactic open clusters of ca .
30 @xmath73 300 = 9000 .
as there are about 150 globular clusters in the galaxy , we have @xmath74 with @xmath75 from above , this gives @xmath76 .
although molecular clouds are not gravitationally bound and are not even in a state of long - term equilibrium , being formed and dispersed on a time - scale shorter than the dynamical time - scale of the galaxy , they are objects in the sense that one can measure their sizes and masses and thus determine a distribution function . molecular clouds exhibit a fractal structure , i.e.  smaller ( and denser ) clouds are embedded in larger ones .
the cumulative size distribution for substructures inside a fractal is given by @xmath77 or in differential form @xmath78 where @xmath79 is the fractal dimension ( mandelbrot 1982 ; see elmegreen & falgarone 1996 , whose formalism we follow here ) . on the other hand , there is a size - mass relation which can be assumed to be a power law again : @xmath80 the size distribution can now be transformed into a mass distribution : @xmath81 observations show that @xmath82 ( see elmegreen & falgarone 1996 ) .
this is consistent with another definition of the fractal dimension : @xmath83 ( mandelbrot 1982 ) , suggesting that the size and mass distribution of molecular clouds is , in fact , the result of fractal gas structure .
the mf for molecular clouds can therefore be written again as @xmath84    to fix the normalization constant @xmath85 we assume that the mean density of molecular hydrogen in the universe is entirely made up of clouds with masses between @xmath86 and @xmath87 ( see elmegreen & falgarone 1996 for the mass range ) : @xmath88 with @xmath89 and @xmath90 from read & trentham ( 2005 ) , we get @xmath91 .
the mfs we have derived so far are plotted together in fig .
we note , as has been noted before , that the mfs of molecular clouds and star clusters ( and , in fact , of stars as well , see below ) have very similar power laws
. such would be a natural outcome , if stars and star clusters form(ed ) by the fragmentation of turbulent cloud cores ( e.g.  elmegreen & efremov 1997 ) .      in the spectrum of cosmic objects we would next have to deal with binary and ( the much less frequent ) multiple stars before coming to single stars .
however , we are still lacking reliable binary statistics .
a very rough estimate of the mf of binary stars is the following .
consider the extreme case that every star is in a system of two stars of equal mass .
then the mass function of binaries per unit volume would obviously be lowered in amplitude by at least a factor of two as compared to the mf of single stars , while otherwise being of similar shape and covering essentially the same mass range .
this shows that in our first - guess universal mass function , binary and multiple star systems can be neglected .
in the following we will determine the mf of a universally averaged ,
_ evolved _
, present - day population of single stars _ and stellar remnants_. such a distribution function has never been discussed in the literature , precisely because it is interesting only in the present ( new ) context .
the mass function of single stars at their birth , i.e. the initial mass function ( imf ) , is an imortant datum for all of astronomy and is of course the backbone of our treatment .
the original form given by salpeter ( 1955 ) : @xmath92 is still in use for stars of higher mass . for low - mass stars various modifications of the salpeter function
have been suggested ( e.g. , kroupa 2001 ) , but there is agreement that well below one solar mass the power - law index is around @xmath93 , rather than @xmath94 . for our purposes we adopt an index of @xmath93 for masses smaller than 0.5 m@xmath95 ( down to @xmath96 = 0.1 m@xmath95 ) , and the salpeter index for stellar masses above that limit ( up to @xmath96 = 120 m@xmath95 ) .    a second ,
even less well - known ingredient for our calculation is the relation between the initial mass of a star and the final mass of its remnant ( white dwarf , neutron star , black hole ) .
stars with initial masses smaller than about 8 m@xmath95 end up as white dwarfs ; stars with initial masses above this critical limit end up as neutron stars or black holes . for the initial - mass range 1 to 7 m@xmath95 and solar metallicity ( which we assume here for simplicity )
we adopt the final white dwarf masses given by weidemann ( 2000 , table 2.1 ) ; they range from 0.55 m@xmath95 for a 1 m@xmath95 star initially , to 1.02 m@xmath95 for a 7 m@xmath95 star initially . for subsolar masses , where little is known about the initial - final mass ratio
, we assume a simple linear relation from ( @xmath97 ) = ( 0.1 , 0.1 ) to ( 1 , 0.55 ) .
the exact relation is nearly irrelevant for our purposes because these stars do not evolve in a hubble time , anyway .    for high - mass stars ( @xmath98 m@xmath95 )
the situation is very complex , as the evolution and therefore also the mass loss of these stars strongly depends on metallicity and other parameters , such as rotation etc .
( woosley , heger & weaver 2002 ; heger et al .
we have constructed an initial - final mass function for these stars , drawing on the information contained in the papers just mentioned , by extending the low - mass relation continuously with a linear relation from ( @xmath97 ) = ( 7 , 1.02 ) up to ( 15 , 1.4 ) , and then again linearly from ( 15 , 1.4 ) all the way up to ( 120 , 25 ) .
this is a very rough approximation ( in lack of better knowledge ) and we assume solar metallicity ( metal - poor stars would suffer less mass loss , resulting in higher - mass remnants ) .
the adopted initial - final mass function is shown in fig .
two more ingredients for our calculation of the mf of stars and stellar remnants for an evolved stellar system are needed : the lifetime of a star as a function of its initial mass , and the star formation rate as a function of time ( i.e.  the star formation history ) . the stellar lifetime ( which is essentially the time spent on the main sequence ) is also a function of metallicity ; for simplicity we again assume solar metallicity .
a table with the stellar lifetimes for solar - metallicity stars with initial mass between 0.8 and 120 m@xmath95 is given , e.g. , in the textbook of sparke & gallagher ( 2000 ) ; for stars with mass below 0.8 m@xmath95 the lifetime is larger than a hubble time .
we have fitted a polynomial of second degree to the tabulated pairs ( log stellar mass , log lifetime ) to get the estimated lifetimes for stars of any ( non - tabulated ) mass .
three types of star formation histories have been used as input : ( 1 ) an initial burst of star formation , i.e.  all stars are born at the same instant , approximately applying to elliptical and s0 galaxies , ( 2 ) a constant star formation rate , simulating the case of spiral and irregular galaxies , and ( 3 ) a `` cosmic '' star formation history .
the latter is inferred from high - redshift observations and is usually given as star formation rate per unit volume as a function of redshift .
we have adopted the log ( sfr)z relation of madau , pozzetti & dickinson ( 1998 , fig .  3 ) and turned it into a log ( sfr)t relation with the help of a `` cosmology calculator '' using @xmath99 km s@xmath100 mpc@xmath100 , @xmath101 and @xmath102 .
the result is shown in fig .
the present - day mf of stars and stellar remnants for an assumed star formation history is now calculated in the following way .
the whole mass range from 0.1 to 120 m@xmath95 is divided into bins of width 0.1 m@xmath95 , giving 1200 mass points .
the evolutionary timespan considered is 12 gyr , i.e.  we assume that the first stars were built around 12 gyr ago , in accord with the typical age of milky way globular clusters .
dividing this timespan into intervals of 0.1 myr gives 12000 timesteps .
consider now that we start with a star formation event .
the original mass function @xmath103 is then simply the adopted imf . after every time step
we compare the time elapsed since the birth event with the lifetime of a star of given mass .
stars with lifetimes smaller than the population age are instantly lost , i.e. taken away from the mass function at their initial mass and added at a different place of the mass function at their much smaller remnant ( final ) mass . in this way stars of
subsequently smaller masses are lost from the population , leading to the well - known `` burning down '' of the imf from high to low masses .
non - coeval stellar populations ( constant sfr or cosmic sfh ) can easily be calculated as superposition of many little bursts distributed over the whole timespan ( one burst for a given population fraction at every timestep ) .
the resulting present - day mfs of stars and stellar remnants for the three different star formation histories considered are shown , along with the imf , in figs .  6 and 7 .
the mfs are normalized to the same integrated , _ initial _ stellar mass as contained in the imf ; the normalization of the imf itself is arbitrary .
the mf would in principle be unchanged , i.e.  remain identical to the imf , for masses somewhat smaller than 1 m@xmath95 .
however , we note a little spike on top of the imf around @xmath104 , or @xmath105 0.6 m@xmath95 .
this is due to the `` graveyard '' of white dwarfs of initial mass 1 m@xmath95 or smaller . the depletion at masses larger than about 4 m@xmath95 ( @xmath106 ) is independent of the star formation history
these are the black hole remnants of rapidly evolving high - mass stars ( initial mass larger than ca .
30 m@xmath95 , see fig .  4 ) .
in the intermediate mass range of 1 to 4 m@xmath95 we see the expected differences : the largest depletion of the mf at supersolar masses is found for the `` initial burst '' case ( i.e.  e and s0 galaxy stellar populations ) , the smallest depletion for constant star formation rate ( spiral and irregular galaxy stellar populations ) , and the cosmic sfh case is lying somewhere between , as the two principal galaxy classes roughly contribute the same amount of stellar mass to the universal mean ( see read & trentham 2005 ) .
the ( artificial ) cut - off at @xmath107 corresponds to the maximum final black - hole mass of ca .
25 m@xmath95 originating from the initial - mass cut - off at 120 m@xmath95 . as the high - mass end of the mf
is entirely shaped by the very poorly constrained final - initial mass function ( fig .
4 ) , this part of the distribution function should not be taken at face value .
overall , from the viewpoint of a universal mass function of things , the present - day mf of stars and stellar remnants is not dramatically different from the stellar imf ; for subsolar masses it is essentially the same , for supersolar masses it is down in amplitude by roughly one order of magnitude .
brown dwarfs seem to cover a very narrow mass range , from the lower mass limit for hydrogen burning at @xmath96 = 0.08 m@xmath95 downwards to about 0.01 m@xmath95 , where the realm of giant planets is entered ( see below ) .
the mf of brown dwarfs is much flatter than that of stars , even of low - mass stars .
luhmann et al .
( 2000 ) and bejar et al .
( 2001 ) , studying brown dwarfs in young open clusters , give @xmath108 for the usual powerlaw form @xmath109 .
we adopt @xmath110 for objects with masses between 0.1 and 0.01 m@xmath95 and attach this distribution continuously to the stellar one at 0.1 m@xmath95 .
the normalization of the stellar mf is straightforward .
read & trentham ( 2005 ) give a stellar contribution to the mean mass density of the universe of @xmath111 .
hence we have @xmath112 where @xmath113 is the mf of stars and stellar remnants ( from 0.1 to 25 m@xmath95 ) as determined above for a cosmic star formation history ( fig .
7 ) , and @xmath114 is the running ( mass - dependent ) power - law index of the same distribution .
brown dwarfs are not included ; their contribution to the total mass is negligible .
in fact , for the integrated mf we could also have neglected the stellar remnants . with @xmath89
this gives a normalization constant of @xmath115 .
the normalized mf of stars and stellar remnants and brown dwarfs appears in figs .  9 to 11 .
inspite of the rapidly growing number of newly detected extrasolar planets , our knowledge of the mass function of planetary and small bodies has still to rely entirely on the census of our solar system .      but even for the solar system we are far from having a complete census .
the recent detection of several pluto - sized objects in the kuiper belt means that completeness of the solar system member census can be claimed only to about a mass of @xmath116 , or a diamter of 1000 km , which are roughly pluto s values .
taking data from the latest edition of _
allan s astrophysical quantities _ ( cox 2000 ) for planets , moons , and bright asteroids , and having turned the diameters of bright asteroids and of some `` unweighted '' moons into masses by assuming a reasonable mean mass density of 2 g @xmath117 , we get a mf of solar system bodys as shown in fig .
8 ( indicated as crosses , giving numbers of objects per unit solar mass ) .
this distribution can be fitted by a straight line representing a rather shallow power - law index of @xmath118 .
as mentioned , below @xmath119 ( pluto s mass ) this mf is bound to be incomplete .
we construct a more reliable low - mass extension from faint asteroids , in the following way .
ivezic et al . (
2001 ) give a diameter distribution function for 13000 asteroids found in the _ sloan digital sky survey_. the diameters had to be inferred from magnitudes by assuming an albedo ; ivezic et al .
( 2001 ) assumed an albedo of 0.14 for reddish ( interpreted as mainly silicate , rocky ) and an albedo of 0.04 for blueish ( mainly carbonaceous ) asteroids .
the resulting diameter ( @xmath79 ) distribution ( see ivezic et al . ,
25 ) can be fitted by a power law of the form @xmath120 where parameter @xmath121 , however , is different for different size regimes : for asteroids larger than 5 km the distribution is very steep with @xmath121 = 4 , for smaller asteroids it is a bit more gentle with @xmath121 = 2.3 .
this seems to be in good accord with other studies of the size distribution of asteroids ( see , e.g. , tedesco , cellino & zappala 2005 ; bagatin 2006 ) .
the size distribution is easily turned into a mf by again assuming a mean mass density @xmath122 = 2 g @xmath117 for the asteroids . with @xmath123
we have @xmath124 the new power - law index @xmath125 for the mf is now @xmath3 for asteroids
larger than 5 km (= more massive than 6.6 @xmath126 ) , and @xmath127 for smaller ( less massive ) bodies .
the mass range covered extends from @xmath128 ( corresponding to asteroids of @xmath129 0.3 km ) up to @xmath130 ( @xmath129 100 km ) .
this is the mf adopted for our study , shown in fig .
lacking proper normalization , it is simply attached to the previously determined mf for planets , moons , and bright asteroids by extrapolating the steep , high - mass power law ( @xmath131 ) towards still higher masses , glueing the two distributions together at @xmath132 , the assumed limit of completeness for the planetary mf . by this
we should more or less include the catalogued , brighter ( larger , more massive ) asteroids missed by the _
sdss_-based work of izevic et al .
( 2001 ; see tedesco et al .
kuiper - belt ( trans - neptunian ) objects seem to follow a distribution similar to the asteroids ( e.g. , bagatin 2005 , fig .
2 ) . by their exclusion , for lack of adequate data , we are likely underestimating the mf amplitude for small bodies , but probably not more than by a factor of 10 , which is still not very relevant for our universal mf .
the exploration of extrasolar planets has only just begun .
what we have learned so far is that the maximum planetary mass seems to be around 10 jupiter masses , or 1/100 of a solar mass  which happens to coincide roughly with the _ minimum _ mass of brown dwarfs ( see above ) .
most exoplanets detected so far have  by necessity , given the detection techniques
masses in the narrow range of 1 to 10 jupiters ( in fact , these are minimum masses , due to the unknown inclination of the planetary orbits ) .
so we can use the ( minimum ) mass distribution of the presently known exoplantes ( around 200 , as of late 2006 ) to extend the solar system mf towards 10 jupiter masses . for this purpose we have collected data published in the web - based _ interactive extrasolar planets catalog _ of jean schneider and glued the narrow mass distribution function directly to the high - mass end of the solar system mf ( i.e. to jupiter at 1/1000 of a solar mass )
; this is the dotted curve in fig .
9 indicated by `` exoplanets '' .
its exact form is irrelevant here , because it will change faster than we could publish it !
the universal abundance of planterary systems is of course completely unknown . for solar - type ( g - type )
stars the frequency of planetary systems is found to be in the 10 % range ( e.g. , quirrenbach 2005 , p.39 ) .
but the detection of extrasolar planets is still biased towards high masses and small periods ; so that freqency could be significantly higher . on the other hand , the majority of stars are m dwarfs , for which planetary companions have yet to be discovered .
the ubiquity of protoplanetary disks around low - mass stars seems to suggest the ubiquity of planetary systems ; but we do not yet know whether this is generally true . nor do we have a clue whether the abundance of planetary and small bodies scales with the mass of the star , or whether the shape of the planetary mf is always the same , etc .
so , for want of better knowledge we make a very simple ( solar - centric ) assumption : we assume that , on average in the universe , there is one solar system for every 10 solar masses ( reproducing , of course , the 10 % frequency for g - type stars ) . since , according to read & trentham ( 2005 ) , the mean stellar mass density of the universe is @xmath133 we simply have to shift the solar system mf per unit solar mass shown in fig .
8 upwards ( enhance the frequency ) by a factor of 3.81 @xmath134 to get a `` universal '' mf of substellar bodies per solar mass and per cubic megaparsec .
this is the distribution shown in fig .
it is apparent that the mfs of giant planets and brown dwarfs do not smoothly interconnect , but there is also no reason why they should .
in fact , there seems to be a genuine frequency gap between the two types of objects ( e.g. , quirrenbach 2005 , p.40 ) . in any case
, it should be borne in mind that this planet - star normalization is probably the most unreliable normalization of our whole study .
we are now in a position that we can put together the mf for objects on superstellar scales ( clusters and groups of galaxies , galaxies , star clusters , and molecular clouds , fig .
3 ) with the mf for objects on the stellar and substellar scales ( stars and stellar remnants , brown dwarfs , planets , and asteroids , fig .  9 ) , to get an overall , `` universal '' mf that stretches from the most massive clusters of galaxies  the largest bound aggregates of matter in the universe at @xmath1  all the way down to sub - kilometer - size bodies of the solar system at @xmath0 (= 20 million tons ) , thus covering 36 ( ! ) orders of magnitude in mass .
this is shown in fig .
10 which constitutes the principal outcome of the present paper .
it is gratifying that the two halves almost perfectly connect to each other around one solar mass .
remember that this normalization was achieved on the basis of the mean universal mass density carried by the stars , embodied in the galaxies on the one hand , and comprised by the stars themselves on the other . before discussing the meaning of this and of the whole mf , we will in the following
add one more piece of information : the mf of dark matter halos as derived from theory and numerical simulations .
a highly advanced theory of structure formation should one day be able to reproduce the observed universal mf .
of course , we are very far from having reached this point . however , on larger scales the basic agent of structure formation is the dissipationless gravitational clustering of matter , a process which is fairly well understood and is simulated , with ever increasing numerical resolution , since more than 30 years . in the context of a cold dark matter ( cdm ) scenario , starting with a harrison - zeldovich spectrum of primordial density fluctuations , there is increasing fluctuation power with decreasing mass at recombination , getting asymptotically flat towards smaller masses ( e.g. , longair 1998 , p.  316 ) .
this means that cdm structure is forming first on small scales and is building up in a hierarchical way to ever larger scales .
it is straightforward to show ( e.g. , longair 1998 , p.373 , see also below ) , by using the formalism of press & schechter ( 1974 ) , that the flat low - mass end of the fluctuation spectrum at recombination results in a power - law mf with slope @xmath3 at small masses , i.e. @xmath135 at the high - mass end , the mf rapidly turns down to zero , giving the functional form a characteristic knee reminiscent of the schechter luminosity function whose form was , in fact , motivated by the press - schechter result ( schechter 1976 ) .
the general form and low - mass slope of the dm mass function has been confirmed and strengthened by more recent semi - analytic work ( sheth & thormen ( 1999 , see also vale & ostriker 2004 ) and by numerous numerical simulations of the hierarchical collapse and clustering of dm halos ( e.g. , reed et al .
2003 ) .
recently , diemand , moore & stadel ( 2005a ) performed very - high - resolution numerical cdm simulations and found dm substructure right down to the theoretical limit at @xmath136 , corresponding to the mass of planet earth ( ! ) .
the mf of these `` microhalos '' turns out to have the same slope of @xmath3 as the mf of more massive halos , at least down to @xmath137 where the function starts to bend down , depending on the nature of cdm ( neutralinos or axions ) .
moreover , in the mass range @xmath138 the normalization ( absolute number of halos per solar mass and cubic mpc ) of diemand et al.s simulation is only a factor of two different from the one given by reed et al .
( 2003 ) when the latter mf is extrapolated down by 10 to 15 orders of magnitude , with a forced slope of @xmath3 , to reach the mass range of diemand et al .
( 2005a , see their fig .  3 ) .    for our study
, in order to have a theoretical counterpart of the mf of galaxies and clusters of galaxies , we adopt a mf for dm halos of the form given in ( 24 ) , i.e.  a slope of @xmath3 for the entire mass range of @xmath139 with the normalization of the reed et al .
extrapolation given in fig .  3 of diemand et al .
( 2005a ) .
as mentioned , below @xmath140 the mf is not known but will probably fall down very rapidly . above @xmath141
the mf bends down in the usual `` schechter - press '' manner ( see e.g. , reed et al .  2003 , vale & ostriker 2004 ) ; for this part of the mf we can directly rely on the cluster data .
the resulting , very simple theoretical mf for dm halos is shown in fig .
11 , along with all other mfs shown already in fig .  10 .
obviously , there is good match on the scale of galaxies and clusters of galaxies .
this had also to be expected , given the fact that , with respect to mass , dm is the major constituent of these structures . in adition ,
the normalization of the simulated mf is of course not independent but is drawn from observed quantities , such as the mean mass density of the universe and the density fluctuation amplitude encapsulated in @xmath32 .
for galaxies , however , we note a difference : the observed slope is flatter than the simulated slope .
this means that the simulations predict more substructure on galactic scales than is actually observed , i.e.  most dm minihalos seem to lack a baryonic ( dwarf galaxy ) counterpart .
this discrepancy has manifested itself also as `` satellite catastrophe '' ( e.g. , moore at al .
the relation of the mini / microhalos on still smaller scales to astronomical objects is addressed below .
here is first a recapitulation of our main findings : + ( 1 ) the mass functions of individual classes of astronomical objects can be concatenated to build a nearly _ continuous _ , `` universal '' mass function that stretches from small asteroids all the way up to the richest clusters of galaxies , covering 36 orders of magnitude in mass .
+ ( 2 ) most individual mfs follow , in part or entirely , a power law with index @xmath142 .
a notable exception is the mf of planetary and small bodies which is significantly flatter . as a consequence the combined mf from stars to clusters of galaxies shows a `` universal '' slope of @xmath142 as well .
+ ( 3 ) cdm halos from numerical simulations follow a nearly identical mf , with @xmath143 , holding even at the scale of stars and planets . + how
should this be interpreted ? what can be expected as a natural outcome of structure formation ? what can not , and is in need of a special explanation ?
certainly the mass function(s ) must somehow reflect the physical processes involved with the formation of these objects .
three different regimes of astronomical objects , known to have radically different formation histories , can be distinguished :    \(1 ) _ galaxies and supergalactic structures _
( from dwarf galaxies to rich clusters of galaxies ) were formed by gravitational clustering in a _ bottom - up _ manner .
their basic constituent is thought to be cdm , being clustered in underlying dark halos .
these halos ( according to the most recent simulations ) have a rich substructure down to the scale of stars and planets .
the standard formation scenario envisions the growth of primordial density fluctuations in cdm and the subsequent , hierarchical gravitational instability and collapse of cdm halos from the smallest ( planetary ) to the largest ( galaxy cluster ) scale .
the @xmath144 mass distribution , which is roughly observed on the galactic and supergalactic scales and is implied by numerical simulations to hold for all scales , is in principle well understood as a consequence of the power spectrum of primordial density fluctuations , @xmath145 , and the `` transfer function '' describing the change of spectrum due to the damping of fluctuations in the pre - recombination epoch , @xmath146 in its most simple form .
the resulting power spectrum at recombination is @xmath147 , and the low - mass end ( @xmath148 ) of the mass function of dm halos at the present epoch becomes ( as inferred , e.g. , from efstathiou 1990 , and longair 1998 ) : @xmath149 with a harrison - zeldovich spectrum ( @xmath150 ) and for adiabatic cdm fluctuations ( @xmath151 in the limit for large @xmath16 ) , this turns into the familiar @xmath144 law . clearly , the primordial spectral index is not as important as the transfer function .
so _ the @xmath144 behaviour can essentially be viewed as a property of cdm_. hot dark matter ( hdm ) , e.g. , would behave very differently  all fluctuations on large scales would get completely damped out and structure formation would ( have to ) turn into a top - down scenario .
\(2 ) _ stars and star clusters _ are believed to have formed by the gravoturbulent fragmentation of molecular clouds in a _ top - down _ manner .
it has been noted before that the mass functions of stars , star clusters , and molecular cloud cores are quite similar . for stars we have the essentially still valid salpeter imf with power - law slope @xmath94 ,
star clusters have mass functions with slope @xmath152 ( see sect .
2.3 ) , and molecular gas clouds show a fractal structure with fractal dimension @xmath153 , corresponding again to a mf of slope @xmath3 ( testi & sargent 1998 , elmegreen & falgarone 1996 ) .
in fact , a fractal structure with dimension between 1.5 and 2.5 is generally found for the interstellar medium , galactic star fields , and even dwarf irregular galaxies ( elmegreen & elmegreen 2001 , parodi & binggeli 2003 ) .
this mf universality has been interpreted as direct evidence for a universal formation mechanism for stars and star clusters in turbulent gas clouds ( elmegreen & efremov 1997 ) . the origin of the scale - invariant clumpy structure of the gas itself is not known ; it could result from gravitational fragmentation , collisional agglomeration , or turbulence ( see references in elmegreen & efremov 1997 ) . nor
is it clear what physical parameters determine the slope of the mf , why it is @xmath3 ; the physics involved is exceedingly complex ( see , e.g. , jappsen et al .
2005 ) .
in the realm of stars and star clusters we should in principle also consider the question of lifetimes , in addition to formation processes ( this point was brought to our attention by the referee of this paper ) .
the lifetimes of stars were explicitely dealt with in our calculation of the mf of an evolved stellar population .
low - mass stars can get older than a hubble time , anyway .
the same is true for high - mass globular clusters .
however , _
open clusters _ are known to have only a lifetime of a few 100 myr ( e.g. , piskunov et al .
so one could indeed wonder why they match the extrapolated mf of stars so well ( see fig .
apparently , their smaller lifetimes ( as compared to the bulk of stars ) are to a certain degree compensated by a higher formation efficiency .
this is in accord with the basic fact that star formation efficiency is generally higher in high - density regions , which of course happen to be the loci of bound cluster formation ( see elmegreen 2006 ) .    in any case
, it is a remarkable coincidence that , on the scale of stars and star clusters , the bottom - up cdm formation scenario should produce the same slope , and even similar amplitude , as the top - down star formation scenario .
the underlying physical processes , even if not exactly known , are expected to be totally different .  or _ is it conceivable that the dm substructure in galaxies has some influence on the fragmentation of the gas , thereby imprinting on the cloud structure the same `` universal '' mass distribution ? _ this is difficult to imagine but should not be excluded without further investigation .
it was questioned whether microhalos could survive the tidal perturbations by other dm halos and , within the galaxy , by stars and molecular clouds . but
recent work ( e.g. , goerdt et al . 2007 ) reaffirms that due to their high density ( cuspy structure ) most microhalos should have survived .
if so , they might act on the gas in a manner that `` conserves '' the mf .
although different processes to explain the stellar mf are usually invoked ( see above ) , minihalos ( not microhalos ) have already been causally connected with old stars and globular clusters , if only in terms of their distribution and kinematics , and not yet their mf ( diemand , madau & moore 2005b ) .
\(3 ) _ planets , moons , and small bodies _ were formed by the coagulation of dust particles and the subsequent accretional growth of small planetesimals ( giant planets additionally by the accretion of gas ) in a _ bottom - up _ manner . however , there are also _ top - down _ processes on these scales . for instance , asteroids formed ( and as a population continuously evolve ) by the collisional fragmentation of large planetesimals .
the formation processes for these objects are surely no less complex than for stars , and we can not even try to explain the mf on the scale of planets and small bodies in detail .
we also recall that the normalization of the planetary mf ( their abundance _ per cubic mpc _ )
remains very uncertain , simply because we know them only very nearby .
however , one property of the mf on these small scales is noteworthy : the general slope appears to be significantly flatter than the `` universal '' @xmath3 .
again we can only speculate on the reason for this difference . in the mass range from stars to the largest structures in the universe , gravitation is certainly the basic agent of structure formation , be it by bottom - up clustering or top - down fragmentation . on smaller scales
, however , gravitation becomes less dominant and electromagnetic forces come into play ; very small asteroids ( rocks ) , e.g. , are no longer held together by gravity but by van der waals ( hence electromagnetic ) forces .
so one could speculate that the break in the slope of the mf around the mass of planets has something to do with the _ transition of the main force responsible for the constitution of objects from gravity to electromagnetic forces_.    in conclusion , the title question of this paper can certainly be answered by yes .
there is a universal mass function in the sense that it is possible to put together a _ continuous _
mass function `` of the universe '' : from asteroids and planets , over stars and stellar remnants , star clusters and gas clouds , galaxies , all the way up to the richest custers of galaxies .
there is also a universal mass function in the sense that the slope of most individual mass functions and that of the grand total is very similar : @xmath154 .
we have no explanation for this nearly universal slope for very different scales and processes .
it is not even clear whether an explanation is needed , i.e.  whether it is more than a chance coincidence
.    it would certainly be interesting to extend the mf of astronomical objects from asteroids further down in mass to meter - size and submeter - size bodies ( which would have to include living things such as the reader ! ) , to dust , molecules , atoms , nuclei , and finally elementary particles . such an  even more speculative
extension was not aimed at here .
we thank gustav tammann , helmut jerjen , pavel kroupa , and ralf klessen for inspiration and fruitful discussions , jay gallagher for a crucial advice , and the anonymous referee for interesting comments .
financial support by the swiss national science foundation is gratefully acknowledged .
bahcall n. , 1979 , apj , 232 , 689 bahcall n. , cen r. , 1993 , apj , 407 , l49 bagatin a.c .
2005 , in lazzaro d. et al .
, eds , iau symp . 229 , asteroids , comets , meteors .
bejar v.j . ,
martin e.l .
, zapatero osorio m.r .
et al . , 2001 , apj , 556 , 830 biviano a. , girardi m. , giuricin g. , mardirossian f. , mezzetti m. , 1993 , apj , 411 , l13 cox a.n . ,
ed . , 2000 ,
allen s astrophysical quantities , fourth edition .
springer , new york diemand j. , madau p. , moore b. 2005b , mnras , 364 , 367 diemand j. , moore b. , stadel j. , 2005a , nat , 433 , 389 efstathiou g. , 1990 , in peacock j.a .
, eds , physics of the early universe .
sussp publications , taylor & francis , p. 394
elmegreen b.g .
2006 , in de koter et al .
, eds , mass loss from stars and the evolution of stellar clusters , asp conf . ser .
, in press ( astro - ph/0610679 ) elmegreen b.g .
, efremov y.n . , 1997 , apj , 480 , 235 elmegreen b.g . ,
elmegreen d.m . , 2001 ,
aj , 121 , 1507 elmegreen b.g .
, falgarone e. , 1996 , apj , 471 , 816 girardi m. , borgani s. , giuricin g. , mardirossian f. , mezzetti m. , 1998 , apj , 506 , 45 goerdt t. , genedin o. , moore b. , diemand j. , stadel j. 2007 , mnras , 375 , 191 harris w.e . , 2001 , in labhardt l. , binggeli b. , eds , lecture notes of saas - fee advanced course 1998 , star clusters .
springer , berlin , p. 305
harris w.e .
, pudritz r.e .
, 1994 , apj , 429 , 177 heger a. , fryer c.l . ,
woosley s.e . , langer n. , hartmann d.h . , 2003 ,
apj , 591 , 288 ivezic l. , tabachnik s. , rafikov r. et al . , 2001 , apj , 122 , 2749 jappsen a .-
, klessen r.s . ,
larson r.b . , li y. , mac low m .- m . , 2005 , a&a , 435 , 611 kroupa p. , 2001 ,
mnras , 322 , 231 longair m.s . , 1998 ,
galaxy formation .
springer , berlin luhmann k.l .
, rieke g.h .
, young e.t .
et al . , 2000 ,
apj , 540 , 1016 madau p. , pozzetti l. , dickinson m. , 1998 , apj , 498 , 106 mandelbrot b. , 1982 , the fractal geometry of nature .
freeman , san francisco moore b. , ghigna s. , governato f. et al . , 1999 , apj , 524 , l19 norberg p. et al ( the 2dfgrs team ) , 2002 , mnras , 336 , 907 parmentier , g. , gilmore , g. 2007 , mnras , in press ( astro - ph/0702258 ) parodi b.r . , binggeli b. , 2003 , a&a , 398 , 501 piskunov a.e
. , kharchenko n.v . , rser s. , schilbach e. , scholz r .- d . 2006 , a&a , 445 , 545 press w.h . , schechter p. , 1974
, apj , 187 , 425 quirrenbach a. , 2005 , in queloz d. et al , eds , lecture notes of saas - fee advanced course 31 , extrasolar planets .
springer , berlin racine r. , 1980 , in hesser j.e . , ed .
, star clusters .
reidel , dordrecht , p. 369
, trentham n. , 2005 , phil .
rsa , 363 , no .
1837 , 2693 reed d. , gardner j. , quinn t. et al . , 2003 , mnras , 346 , 565 salpeter e.e . , 1955 , apj , 121 , 161 schechter p. , 1976 ,
apj , 203 , 297 sheth r.k . , tormen g. , 1999 , mnras , 308 , 119 sparke l.s .
, gallagher j.s . , 2000 , galaxies in the universe .
cambridge university press , p. 10
surdin v.g . , 1979 , sva , 23 , 643 tedesco e.f . , cellino a. , zappala v. , 2005 , aj , 129 , 2869 testi l. , sargent a.i . , 1998 , apj , 508 , l91 vale a. , ostriker j.p . , 2004 , mnras , 353 , 189 van den bergh s. , lafontaine a. , 1984 , aj , 89 , 1822 van den bosch f.c .
, yang x. , mo h.j . ,
norberg p. , 2005 , mnras , 356 , 1233 whitmore b.c . , schweizer f. , 1995 , aj 109 , 960 woosley s.e . , heger a. , weaver t.a .
, 2002 , rev .
phys . , 74 , 1015 zwicky f. , 1942 , phys .
review , 61 , 489 zwicky f. , 1957 , morphological astronomy , springer , berlin