it is by now well established that baryonic matter represents only about @xmath1 of the total matter density of the universe ( e.g. * ? ? ?
* ) , while the majority is in the form of non - baryonic dark matter ( dm )
. since galaxies form through the accretion of baryonic material onto dynamically dominant dm structures ( halos ) , it would be reasonable to assume that the baryon mass fraction of present day galaxies approximately equals the cosmic value ( @xmath2 ) . despite this expectation , observations point to the fact that galaxies are not able to retain their cosmic `` fair share '' of baryons , and that the resulting baryon deficit depends strongly on the mass of their host halo .
the first line of evidence is provided by observational estimates of the abundance of galaxies as a function of their total stellar mass , a distribution referred to as the galactic stellar mass function ( smf ) .
thanks to the advent of wide area optical surveys with multiband photometric and spectroscopic information , such as the two degree field galaxy redshift survey ( 2dfgrs ) and the sloan digital sky survey ( sdss ) , the smf has been measured over the mass range @xmath3 , using statistical samples of tens of thousands of galaxies and a variety of stellar mass estimation techniques ( @xcite , @xcite , @xcite , @xcite , @xcite , @xcite , @xcite , to name a few ) .
the smf displays an exponential cutoff at masses @xmath4 and an approximate power - law behavior at low masses ( @xmath5 ) , with a `` shallow '' exponent of @xmath6 . on the other hand
, the halo mass function ( hmf ) predicted in the lambda cold dark matter ( @xmath0cdm ) model , follows a much `` steeper '' power - law ( @xmath7 ) over the mass range of interest .
this observation alone excludes the possibility that the stellar mass of a galaxy is simply a fixed fraction of the host halo mass .
one can furthermore statistically derive an average relation between the stellar mass of a galaxy ( @xmath8 ) and the mass of its host halo ( @xmath9 ) , through the technique of abundance matching ( see  [ sec : am ] for details ) .
@xmath8 - @xmath9 relations based on abundance matching ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) have shown that the `` stellar conversion efficiency '' , @xmath10 , never exceeds 25 - 30% . furthermore , @xmath11 peaks for milky way - sized galaxies ( @xmath12 ) , and declines rapidly on either side of the peak ( e.g. figure 2 in @xcite ) .
the second line of evidence comes from direct halo mass measurements , obtained through weak lensing or kinematics studies ( e.g. * ? ? ?
* ; * ? ? ?
for example , @xcite used stacked weak lensing measurements to estimate the average host halo mass of galaxies in different stellar mass bins , and found that @xmath11 never exceeds @xmath13% .
direct halo mass measurements can circumvent a number of assumptions inherent in the application of abundance matching , but such techniques can presently only be applied to a restricted range of stellar mass ( @xmath14 ) , and are affected by their own set of systematics .
stellar mass is not always the dominant baryonic component in a galaxy . in fact , the hi  to  stellar mass ratio ( `` hi fraction '' ; @xmath15 ) tends to increase with decreasing stellar mass , and hi often dominates the baryonic content of low - mass galaxies .
the transition from stellar - mass  dominated to hi
dominated systems takes place at @xmath16 for hi - selected samples ( e.g. * ? ? ?
* see also fig .
[ fig : gsmr ] in this work ) , or at @xmath17 for optically - selected samples ( e.g. * ? ? ?
as a result , it is presently not clear what is the behavior of the `` baryon retention fraction '' @xmath18 in low - mass galaxies , when both stars and cold gas are taken into account . in particular
it is not well understood whether the very low average value of @xmath11 inferred for low - mass halos is a result of poor retention of baryonic material , of the low efficiency of gas - to - stars conversion , or of a combination of both .
for example , @xcite argue that the increasing gas fraction in low - mass galaxies should approximately offset the decreasing stellar - to - halo mass ratio , and result in a roughly constant @xmath19% .
this conclusion was based on an indirect estimate of the cold gas content of galaxies , based on the average @xmath20 relation observed in a set of samples in the literature
. an early work by @xcite , based on the same indirect method , also reached a qualitatively similar conclusion .
@xcite found an approximately constant @xmath21 at the low - mass end using a different indirect method , which involves the comparison of the stellar and hi mass distributions of two different galaxy samples .
these results would imply that low - mass galaxies are relatively efficient at retaining baryonic mass , but very inefficient in converting their gas into stars .
this conclusion , however , would require a `` steep '' hi mass function ( himf ) in the local universe , in contrast to what is measured @xcite .
moreover , the recent work of @xcite , also based on using the average @xmath20 relations for blue and red galaxies separately , found no signs for a flat @xmath21 at low masses .    in this article
we _ directly _ measure the abundance of galaxies as a function of their `` baryonic mass '' ( throughout this article the term baryonic refers to the combined stellar and atomic gas components of galaxies , and baryonic mass is calculated as @xmath22 , where the 1.4 factor accounts for the presence of helium ) .
we use optical data from the seventh data release of the sdss ( sdss dr7 ) to estimate stellar masses , and hi - line flux measurements from the arecibo legacy fast alfa ( alfalfa ) survey to measure atomic gas masses .
the resulting distribution , referred to hereafter as the baryon mass function ( bmf ) of galaxies , can be used in abundance matching to derive a robust @xmath21 - @xmath9 relation . in order to investigate sample selection effects , we employ both an hi - selected and an optically - selected sample drawn from the _ same _ volume to derive the mass distributions for the stellar , atomic hydrogen and baryonic components .
the paper is organized as follows : in section [ sec : data ] , we introduce the datasets used to measure the stellar , hi and baryon mass distributions .
we describe the methodology used to measure atomic hydrogen masses and we estimate stellar masses for our galaxy samples . in section [ sec : bmf ] , we present our measurements of the smf , himf & bmf from both the hi - selected and the optically - selected samples , and compare them against one another as well as against other published results . in section [ sec : biases ] we consider the impact of possible systematics on our measurements , such as stellar mass estimation method , distance uncertainties and the exclusion of some baryonic components ( e.g. molecular gas ) in the calculation of the bmf . in section
[ sec : bar_frac ] , we present the @xmath11 - @xmath9 and @xmath21 - @xmath9 relation in a @xmath0cdm universe . in section [ sec : conclusion ] , we discuss the implications of the result and summarize our main conclusions . throughout this paper
, we use a hubble constant of @xmath23 .
we select galaxies from the current data release of the alfalfa survey , which covers 40% of the planned final survey area ( `` @xmath24.40 '' catalog ; @xcite ) .
we restrict ourselves to two rectangular areas of the `` spring '' coverage of @xmath24.40  ( @xmath25 , @xmath26 & @xmath27 ) , which encompass the virgo cluster as well as the supergalactic plane at low velocities .
we restrict ourselves to galaxies with @xmath28 15000 km s@xmath29  ( @xmath30 ) , in order to avoid the strong radio frequency interference ( rfi ) present at frequencies that correspond to @xmath31 km s@xmath29 .
we discard the nearest extragalactic sources with @xmath32 mpc , because they can carry extreme fractional uncertainties on their distances ( see  [ sec : derived ] for details on the distance assignment method ) .
we furthermore select only hi sources designated as `` code 1 '' in @xmath24.40 , i.e. extragalactic sources detected at high significance ( @xmath33 ) .
in addition , we exclude sources with integrated fluxes below the 50% completeness limit of the alfalfa survey ( see @xcite , section 6 for the derivation of the alfalfa completeness limits ) .
the above requirements are satisfied by 7618 galaxies .
we remove from our sample 204 @xmath24.40  sources which are not crossmatched with an optical source in sdss , as well as 208 additional sources which have been flagged as having problematic sdss photometry ( crossmatch code `` p '' in @xmath24.40 , for details see section 4 in @xcite ) .
this quality cut on the sdss photometry introduces some bias against faint , low surface brightness galaxies of irregular morphology ; such sources are often `` shredded '' ( i.e. assigned multiple photometric objects ) by the sdss magnitude extraction process , and are usually assigned a `` p '' ( `` photometry suspect '' ) crossmatch code in @xmath24.40 .
lastly , 11 additional objects were discarded , in cases where the stellar mass computation method described in  [ sec : derived ] failed to produce physically plausible results .
our final sample thus consists of 7195 extragalactic objects , detected over @xmath342000 deg@xmath35 of high galactic latitude sky and out to @xmath36 mpc .
the upper panel of figure [ fig : coneplot ] displays the spatial distribution of our hi - selected galaxies , and puts in evidence the complex large scale structure in the survey volume .
note that all objects in our hi - selected sample have 21 cm redshifts and line fluxes as well as multi - band optical photometry , and hence estimates of both their stellar and atomic hydrogen masses .
we draw an optically - selected sample from the sdss dr7 @xcite spectroscopic database , in the same volume used to define our hi - selected sample . specifically , we select galaxies that lie within the same sky area ( @xmath25 , @xmath26 & @xmath27 ) and satisfy the same velocity and distance restrictions ( @xmath37 15000 km s@xmath29 , @xmath38 mpc ) .
cmb velocities for our optically - selected galaxies are calculated based on their sdss spectroscopic redshifts ( @xmath39 ) .
we restrict ourselves to objects spectroscopically classified as galaxies in sdss ( ` specclass = 2 ` ) that also have an apparent petrosian magnitude brighter than 17.5 in the @xmath40-band ( @xmath41 ) .
this initial cut results in 22707 galaxies .
due to their large number , it is not practical to inspect all galaxies individually for the quality of their sdss photometry / spectroscopy . as a result
, we expect a fraction of our sources to have issues with their sdss photometry , most often related to `` shredding '' ( i.e. assignment of multiple photometric objects to a single galaxy ) .
this issue affects mostly extended sources with structure in their light distribution , such as low surface brightness ( lsb ) galaxies with irregular morphology . in such cases
, the sdss magnitude will underestimate the true flux of the galaxy , which in turn will result in an underestimate of its stellar mass . a second issue related to shredding ,
is that bright star forming knots in the disks of nearby spiral galaxies can sometimes be cataloged as separate spectroscopic objects , and hence interpreted as low - mass satellites of the main spiral .
we find that applying a color cut on our sample , @xmath42 , removes a fair fraction of these unwanted cases . on the other hand , cuts based on the quality of the sdss spectrum ( such as cuts on ` zconf ` , ` zstatus ` or ` zwarning ` ) are ineffective , since they exclude mostly legitimate faint or lsb dwarf galaxies with noisy spectra .
lastly , we exclude objects for which the stellar mass computation described in  [ sec : derived ] failed to produce physically plausible results .    our final optically - selected sample consists of 22587 galaxies , occupying the same volume as our hi - selected sample .
we crossmatch the optical sample with the full @xmath24.40  catalog ( including code 1 & 2 sources ) , and find 7551 hi source counterparts .
the crossmatch rate is thus approximately @xmath43 , as reported in @xcite .
note that alfalfa non - detected galaxies are not necessarily hi - poor objects ; due to the low emissivity of atomic hydrogen in the 21 cm line , even moderately gas - rich galaxies in the outer portion of the survey volume can be missed by alfalfa .
this point is illustrated by the lower panel of figure [ fig : coneplot ] , which compares the spatial distribution of galaxies in the hi - selected and optically - selected samples .
note that all optically - selected galaxies have multiband optical photometry as well as optical redshifts , and hence an estimate of their stellar mass .
however , only galaxies crossmatched with an @xmath24.40  source have a 21 cm flux measurement , and hence a value for their atomic hydrogen mass
.      -band absolute magnitude ( _ panel a _ ) , @xmath44 color ( _ panel b _ ) , stellar mass ( _ panel c _ ) and hi mass ( _ panel d _ ) , for the optically - selected ( _ red solid line _ ) and hi - selected ( _ blue dashed line _ ) samples . as evident in panel _
b _ , the hi - selected sample is strongly biased against the red galaxy population and as a result it is skewed towards lower luminosity and stellar mass systems .
conversely , the fractional contribution of bright and massive galaxies ( @xmath45 ) is larger for the optically - selected sample ( panels _ a _ and _ c _ ) . ]
we calculate hi masses from the measured 21 cm integrated flux reported in @xmath24.40 . assuming optically thin emission    @xmath46    where @xmath47 is the hi mass in units of the solar mass ( @xmath48 ) , @xmath49 is the integrated flux in jy@xmath50km@xmath51s@xmath29  and @xmath52 is the distance in mpc . distances in this article are calculated according to the method used in @xmath24.40  @xcite : nearby galaxies ( @xmath53 km s@xmath29 ) are assigned distances through the use of a peculiar velocity flow model developed by @xcite , while for more distant galaxies simple hubble distances are used ( @xmath54 , with @xmath55 km@xmath51s@xmath29mpc@xmath29 ) .
moreover , group and cluster membership information is taken into account when available , as well as primary distance measurements published in the literature .
we would like to point out that most of the galaxies in our optically - selected sample are not included in @xmath24.40 , and hence lack the systematic group assignments and primary distance information contained in the catalog .
nevertheless , optically - selected galaxies that lie within the sky area and redshift range of the virgo cluster are placed collectively at the virgo distance ( @xmath56 mpc ) , in order to minimize the effects of peculiar motions on the inferred distances of galaxies in the region .
we would also like to note that the distance assignment method can have a large impact on the determination of mass functions , especially at the low - mass end .
we illustrate this issue in  [ sec : dist_unc ] , where we consider the effect on the himf of using uniformly hubble distances for all galaxies .
we compute stellar masses for our galaxies based on fitting all 5 sdss photometric bands ( @xmath57 ) , with model spectral energy distributions ( seds ) .
the full details of the method can be found in @xcite , but here we summarize the main points : a library of model seds are generated , using the @xcite stellar population synthesis code and assuming a @xcite stellar initial mass function ( imf ) .
models with an extensive range of internal extinction , metallicity and star formation histories are considered .
in particular , star formation history templates include both an exponentially declining component as well as random starburst episodes .
the final physical properties ( e.g. stellar mass , star formation rate , internal extinction etc . ) are computed as the average of all model values , where each model is weighted according to its fit likelihood .
in addition to mean values , `` @xmath58 '' uncertainties of the physical properties can also be derived , as one quarter of the 2.5 - 97.5 percentile range of model values .
the median @xmath58 uncertainty in @xmath59 is 0.086 dex , or about 22% ( excluding uncertainties on the distance ) .
it is important to note that stellar mass estimates of the same galaxy obtained with different methods can have systematic offsets of up to factors of a few . in
[ sec : stellar_est ] we address issues related to stellar mass estimation , and consider alternative methods for calculating stellar masses @xcite .    figure [ fig : histogram ] compares the distributions of @xmath60-band absolute magnitude ( @xmath61 ) and @xmath44 color ( both corrected for milky way extinction ) , stellar mass ( @xmath8 ) and hi mass ( @xmath47 ) for the hi - selected and optically - selected samples .
the most notable difference is in the @xmath44 color distribution , where the hi - selected sample shows a strong bias against the red galaxy population ( see also @xcite ) . as a result
, the optically - selected sample contains a larger proportion of high luminosity and stellar mass systems compared to the hi - selected sample .
by contrast , the @xmath47 distribution of the two samples is very similar , but remember that only those optically - selected galaxies that are detected in alfalfa are included in the histogram .
stellar masses for all galaxies , and hi masses for all alfalfa - detected galaxies are calculated as described in
[ sec : derived ] . for the @xmath3415000 galaxies in the optically - selected sample that lack an alfalfa detection , we assign a lower and an upper limit on their atomic hydrogen content ( @xmath62 ) .
the lower limit is simply @xmath63 , which corresponds to an hi - devoid galaxy .
the upper limit is calculated by assuming that the hi flux of the non - detected galaxy lies just below the alfalfa `` detection limit '' , as defined by the 25% completeness limit of the @xmath24.40  catalog when both code 1 & 2 sources are considered .
more precisely    @xmath64    where @xmath52 is the galaxy distance in mpc determined by its sdss optical redshift , and @xmath65 is the flux level at which the completeness of the @xmath24.40  catalog falls to 25% , in jy km s@xmath29 .
according to eqns .
6 & 7 of @xcite    @xmath66    where @xmath67 is the hi - line profile width in km s@xmath29 , measured at the 50% flux level of the profile peak .
since alfalfa non - detected galaxies lack a measurement of @xmath67 , we assign a value based on the average @xmath8 @xmath68 relation ( i.e. the stellar mass tully - fisher relation ) of @xmath24.40   galaxies .
we then project the @xmath68 value on the line - of - sight according to the sdss @xmath40-band axial ratio , and assuming an intrinsic axial ratio of @xmath69 for all galaxies .
baryonic masses ( i.e. stellar mass + atomic gas mass ) for all galaxies are calculated as @xmath22 , where the 1.4 factor is used to account for the cosmic abundance of helium .
note that alfalfa non - detected galaxies have two assigned values for their hi mass , and consequently two values for their baryonic mass , @xmath70 and @xmath71 .
we calculate cumulative mass functions in logarithmic mass bins for all three components ( i.e. stellar mass , atomic hydrogen mass , baryonic mass ) , separately for the hi - selected and optically - selected samples . since neither sample is volume - limited , mass functions have to take into account the sample selection criteria as well as the large - scale structure in the survey volume .
hi selection is based on a combination of galactic hi integrated flux , @xmath49 , and profile width , @xmath67 ( see  [ sec : hi_selected ] & discussion in section 6 of @xcite ) ; as a result , galaxies of different hi masses and linewidths are detected out to different distances .
similarly , our optically - selected sample is a flux - limited sample , which results in galaxies with different @xmath40-band absolute magnitudes being detected in different volumes . as a result ,
mass distributions are calculated by summing up the number of detections in a given mass bin ( see fig .
[ fig : histogram ] ) , with each detection weighted by an appropriate volume factor .
individual weighting factors are calculated via the `` @xmath72 '' method , as implemented in @xcite .
this is a non - parametric , maximum - likelihood method , which reduces to the standard @xmath73 method @xcite when applied to a spatially homogeneous galactic sample .
the advantage of @xmath72 consists in the fact that it is insensitive to local density fluctuations , and hence mostly immune to structure - induced bias .
we refer the interested reader to the literature for the definition , implementation and details of the method .
the method definition and basic setup is described in @xcite ; details of the implementation of the method for hi data can be found in @xcite ; the specifics of the application of the method on alfalfa data can be found in @xcite , while a shorter qualitative description can be found in @xcite .
lastly , a fraction of galaxies that satisfy all criteria for sdss spectroscopic followup can not be observed for technical reasons ( mostly fiber collisions ) , and are therefore not included in the sdss spectroscopic database .
we therefore correct the normalization of all optically - selected distributions by @xmath74 , using the average spectroscopic completeness value reported in @xcite , @xmath75 .
similarly , a fraction of the alfalfa volume is `` lost '' due to rfi contamination of certain frequency bands in the alfalfa passband .
we correct the normalization of all hi - selected distributions by @xmath76 , where @xmath77 .
we would also like to note that , due to the @xmath78 beam size of the alfa receiver , a number of hi sources are expected to be blended .
we do not attempt to correct for blending but , given that hi - selected galaxies are a weakly clustered population ( e.g. * ? ? ?
* ) , we anticipate the effect on the hi - selected distributions to be small .
( bmf , _ black symbols _ ) , derived from the hi - selected galaxy sample .
error bars represent just the poisson counting error assuming independent errors among different mass bins . ]     ( _ black lines _ ) , derived from the optically - selected galaxy sample .
the atomic hydrogen and baryonic distributions are represented as allowed ranges , based on estimates of the minimum and maximum hi mass for galaxies that are not detected by alfalfa ( \{@xmath62 } , see  [ sec : derived ] for details ) .
error bars again represent just the poisson counting error assuming independent errors among different mass bins . ]
figure [ fig : mfs_hi ] shows the cumulative distribution of stellar mass ( smf ; _ gold symbols _ ) , hi mass ( himf ; _ cyan symbols _ ) and baryonic mass ( bmf ; _ black symbols _ ) , derived from the hi - selected galaxy sample .
the hi - selected bmf follows closely the hi - selected smf at high masses , while at the low - mass end the contribution of the himf becomes dominant ; this is because hi - selected galaxies become more gas - rich as their stellar mass decreases .
figure [ fig : mfs_opt ] shows the same distributions ( smf , _ gold line _ ; himf , _ cyan lines _ ; bmf , _ black lines _ ) derived from the optically - selected sample .
recall that , in the case of the optically - selected sample , a lower and upper limit of the himf and bmf are shown , since sdss galaxies that are undetected by alfalfa are assigned an upper and lower limit on their hi content , and therefore also on their baryonic mass ( see  [ sec : derived ] ) .    ) and low - mass ( @xmath79 ) ends .
this difference is mostly due to the bias of the hi - selected sample against the red galaxy population ( see  [ sec : results ] for a detailed discussion ) . ]     for a detailed discussion . ]    ) .
[ sec : results ] for a detailed discussion . ]    in figure [ fig : opt_vs_hi_gsmf ] we compare the smfs derived from the optically - selected and hi - selected samples .
the two smfs are consistent at intermediate masses , but the optically - selected smf ( gold line ) is systematically higher at the high - mass and low - mass ends . at high masses
the discrepancy is due to the bias of the hi - selected sample against the most massive galaxies , which are usually red passive systems .
the discrepancy at the low - mass end is mostly due to the population of red - sequence dwarf galaxies in the nearby virgo cluster that are undetected by alfalfa ; these are mostly dwarfs with early - type morphologies and very low hi content ( see @xcite , for example ) . on the other hand , the hi - selected and optically - selected himfs ( fig .
[ fig : opt_vs_hi_himf ] ) are in good agreement with one another , with the hi - selected himf having a slightly steeper low - mass end slope than what suggested by the range of the optically - selected distribution .
the two bmfs are mostly in agreement with one another , except at the high - mass end ( a factor of @xmath344 at @xmath80 ) .
this is a direct consequence of the discrepancy between the optically - selected and hi - selected stellar mass functions at high masses .
note that there is little difference between the two bmfs at low masses , which suggests that the low - mass end of the bmf has been measured robustly .      ) .
the red diamonds correspond to the @xcite smf in the local universe ( @xmath81 ) , while the purple triangles correspond to the smf of @xcite , extracted over a larger volume ( @xmath82 ) . both the baldry at el . and
yang et al .
smfs are based on the nyu - vagc galaxy catalog .
[ sec : comparison ] for a discussion of the comparison . ]
figure [ fig : gsmf_vs_other ] compares the optically - selected smf presented in this work ( same as gold line in fig .
[ fig : opt_vs_hi_gsmf ] ) with the local - universe smf of @xcite and the @xcite smf , which are both based on the new york university value - added galaxy catalog ( nyu - vagc ; * ? ? ?
there is excellent agreement between the baldry et al .
smf and our optically - selected smf , especially at intermediate and low stellar masses ( @xmath83 ) .
the deviations at high masses are due to the fact that stellar masses in baldry et al .
are calculated differently than in this work ( see sec . 3 of @xcite for details ) ; note that a systematic difference of just 0.1 dex ( 26% ) in stellar mass would be enough to explain the observed difference in abundance at the high - mass end .
the @xcite smf is systematically higher than our optically - selected smf at high masses , and displays a more pronounced `` plateau '' at intermediate masses .
it is not clear what the cause of the difference at the high - mass end of the distributions is , but it may relate to the fact that the yang et al .
smf is extracted from a significantly larger volume than our measurement ( the maximum redshift is @xmath84 for the yang et al . sample and @xmath85 for the sample used in this work ) .
moreover , yang et al . use the prescription of @xcite to estimate stellar masses , which is based on the galactic @xmath86 color . as discussed in more detail in
[ sec : stellar_est ] , the use of different stellar mass estimators can significantly affect the shape of the measured smf .    ) .
the solid purple line corresponds to the schechter function fit to the @xcite himf , which is based on the full @xmath24.40  catalog of alfalfa sources and without any optical selection cuts .
the green dashed line corresponds to the schechter fit to the @xcite himf , based on 4315 galaxies detected by the hipass survey .
[ sec : comparison ] for a discussion of the comparison . ]
figure [ fig : himf_vs_other ] compares the hi - selected himf presented in this work ( same as the cyan symbols in fig .
[ fig : opt_vs_hi_himf ] ) with the himf of @xcite derived from 10119 galaxies detected by alfalfa over @xmath342600 deg@xmath35 of sky ( purple solid line ) .
there is excellent agreement at intermediate and large hi masses ( @xmath87 ) between the @xcite himf and the himf derived in this work , while at lower masses the @xcite himf is slightly steeper than ours .
this disagreement may be due to the set of additional optical requirements imposed on our hi - selected sample .
as argued in  [ sec : hi_selected ] , these requirements are expected to reduce the number of low - mass systems in the sample and therefore decrease the inferred space density at the low - mass end .
the dashed green line represents the himf of @xcite based on 4315 sources detected by the hi parkes all - sky survey ( hipass ) over the whole southern celestial hemisphere ( @xmath34 29000 deg@xmath35 ) .
there is disagreement between both alfalfa - based himfs and the hipass - based himf of @xcite , at the high - mass end . as argued in @xcite , the higher sensitivity of the alfalfa survey compared to hipass , which enables alfalfa to detect hi - massive systems over a larger volume , should give a statistical advantage to the alfalfa survey in determining the high - mass end of the himf .
however , the difference is too large to be explained by counting statistics or cosmic variance ( e.g. according to the estimates of @xcite or @xcite ) . on the other hand , due to the exponential drop - off of the himf at high masses ,
a flux calibration difference of as low as 0.1 dex could give rise to a similar discrepancy .
color and the @xmath60-band luminosity , according to the widely used @xcite calibration .
the gold dotted line represents the smf computed using stellar masses derived from the galactic @xmath44 color and the @xmath60-band luminosity , according to the the more recent calibration of @xcite .
[ sec : stellar_est ] for a discussion of the comparison . ]    a variety of methods exist to estimate stellar masses from spectra or broadband photometric measurements of galaxies ( * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* to name a few ) .
most methods rely on comparing the actual galactic emission to the light output of a set of galactic stellar population models .
the models that most closely reproduce the observed data are then used to estimate the galactic properties of interest ( e.g. stellar mass , present star formation rate , internal extinction etc . ) ; it is therefore very important to consider model stellar populations which span as large a range of physical parameters as the galaxies in the sample being studied .
inferred galactic properties depend not only on the particular type of data employed by each method ( e.g. spectroscopy vs. broadband photometry or optical vs. near - infrared photometry ) , but also by differences in the way in which the model stellar populations are constructed .
this means that different methods can yield different estimates of a galactic property even when the _ same _ observational measurements are used .
for example , @xcite find that unbiased stellar masses can only be recovered when the true star formation history ( sfh ) of a galaxy is known . in practice
however a restricted set of sfhs is considered ( often in the form of a parametrized function ) , which may introduce systematics for galaxies with sfhs that are not well described by the assumed general form .
additional complications can be introduced by the different treatment of dust reddening among different models .
in general , stellar mass estimates can differ systematically by as much as 0.3 dex , while for individual galaxies the scatter can be as large as 0.6 dex @xcite .
bands used in this work ( see  [ sec : derived ] ) , and those derived from the galactic @xmath86 color and the @xmath60-band luminosity according to the widely used @xcite calibration .
each datapoint corresponds to a galaxy in the optically - selected sample , and the symbol color represents the g - r color of the galaxy .
the two stellar mass estimates agree fairly well for red passive galaxies , while for blue star - forming galaxies bell et al .
masses are systematically larger by up to a factor of @xmath342 ( see  [ sec : stellar_est ] for a detailed discussion ) . ]    here we compare the optically - selected smf presented in this work ( computed from stellar masses derived from sed - fitting , see  [ sec : derived ] ) with the smf obtained using stellar masses derived from a single galactic color , according to the widely - used @xcite calibration as well as the more recent calibration of @xcite .
more specifically , we compute bell et al
. masses by multiplying the @xmath60-band luminosity of each galaxy by a mass - to - light ratio inferred from its @xmath86 color .
we choose this particular combination of bands because it is relatively immune to contamination of galactic spectra by bright nebular emission lines @xcite .
we use a similar procedure to calculate taylor et al .
masses , by using their calibration of @xmath60-band mass - to - light ratio versus @xmath44 color .
@xcite argue that using @xmath44 colors best constrains the galactic stellar mass estimates .
figure [ fig : sed_vs_color ] shows the impact that different methods of estimating stellar mass have on the measurement of the smf .
when bell et al .
stellar masses are used ( _ gold dashed line _ ) , the smf becomes systematically higher at low and intermediate masses ( @xmath88 ) , while it remains mostly unchanged at the high - mass end .
the reason for this pattern becomes evident in fig .
[ fig : bell_vs_sed ] , where we see that bell et al .
masses agree fairly well with the masses derived in this work for red passive galaxies , but are systematically larger ( by up to a factor of @xmath342 ) for blue star - forming galaxies .
@xcite argue that the difference can be primarily attributed to the fact that the stellar population models used for the bell et al .
calibration do not consider `` bursty '' star formation histories which are typical of low - mass galaxies with blue colors .
this leads to systematically older stellar populations for blue galaxies according to the bell et al .
method , which in turn results in systematically higher stellar mass estimates .
note that including models with bursty sfhs in a stellar population library does not by itself guarantee a correct estimate of stellar mass ; overestimating the effect of bursts would result in systematically low stellar masses for blue galaxies .
conversely , when taylor et al .
masses are used ( gold dotted line ) , the smf becomes systematically lower at the high - mass end , while it is mostly unchanged at low and intermediate masses .
again , this is a result of the fact that taylor et al .
masses agree well with the sed - fitting masses used in this work for blue star - forming galaxies , but are systematically lower ( by up to a factor of @xmath341.4 ) for red passive galaxies .       for further discussion ) . ]
stellar , hi and baryonic masses are distance - dependent quantities and hence their statistical distributions are affected by distance uncertainties .
this is particularly true for the low - mass end of the distributions , which is determined by the properties of nearby galaxies ; neglecting the peculiar velocity of some of these objects can cause fractional distance and mass errors of order @xmath89 , especially in a volume with complex large - scale structure such as the one surveyed by alfalfa ( see fig . [
fig : coneplot ] ) .    for this reason ,
nearby galaxies ( @xmath90 km s@xmath29 ) in the @xmath24.40  catalog are assigned distances based on a parametric peculiar velocity flow model @xcite , and only more distant galaxies ( @xmath91 km s@xmath29 ) are assigned simple hubble distances according to their cmb recessional velocity . the @xcite flow model includes a dipole and a quadrupole component ( local group bulk motion & asymmetric expansion ) and two local attractors ( virgo cluster & great attractor ) , and is calibrated against the sfi++ catalog of galaxies with redshift and independent distance distance information ( from tully - fisher ) .
the residuals are then attributed to random thermal motions , which are estimated to have a magnitude of @xmath92 km s@xmath29 . in addition , distances reported in the @xmath24.40  catalog take into account known group and cluster membership as well as primary distance information published in the literature .
this latter information is not available for the majority of the galaxies in our optically - selected sample ( which are not included in @xmath24.40 ) , and we only make an attempt to assign all probable virgo members to the virgo cluster distance ( @xmath93 ) .    here
we re - evaluate the himf for our hi - selected sample using uniformly hubble distances for all galaxies , in order to illustrate the impact of the distance assignment scheme on the derived distributions .
figure [ fig : hubble_vs_flow ] shows that the himf computed using hubble distances ( cyan dashed line ) has a much shallower low - mass end slope compared to the himf presented in this work , which uses flow model distances for most nearby galaxies ( cyan solid line ) .
this result is in agreement with the work of @xcite , who find that neglecting the local peculiar velocity field when estimating distances to nearby galaxies in the alfalfa volume will lead to a systematically shallower low - mass end slope for the himf .      .
cyan diamonds represent 14 galaxies of the heracles survey @xcite detected in co line emission .
the solid and dotted blue lines represent the average and @xmath94 scatter of the distribution found for 125 co detected galaxies in the cold gass survey @xcite .
both samples show that , at least for @xmath95 , molecular gas is almost always a subdominant mass component ( see  [ sec : h2 ] for details ) . ]    throughout this article we use the term `` baryonic mass '' to refer to the sum of the stellar and atomic gas mass ( @xmath22 ) , a convention that is common in the literature .
this definition however excludes a number of baryonic components that are definitely present in galaxies , most notably molecular and ionized -warm or hot- gas .
figure [ fig : h2_mass_frac ] displays the fraction of molecular hydrogen ( h@xmath96 ) mass ( accounting for helium ) to `` baryonic mass '' as defined above , as a function of stellar mass .
the cyan diamonds represent 14 galaxies from the heracles survey @xcite with h@xmath96  masses measured from interferometric co line observations ( using a fixed @xmath97 conversion factor ) . over the probed stellar mass range ( @xmath98 ) , molecular gas rarely contributes more than 10% of the `` baryonic mass '' .
the same conclusion is reached when h@xmath96  mass measurements from the cold gass survey are considered @xcite .
the blue solid and dotted lines show the mean and @xmath94 scatter of the relation between the molecular and `` baryonic '' mass components , based on 125 galaxies detected in co emission with the iram 30 m telescope . in this latter case
we have estimated the atomic hydrogen mass of galaxies indirectly , using the average @xmath99 vs. @xmath8 relation of the cold gass parent sample @xcite .
again , over the stellar mass range probed by the survey ( @xmath100 ) molecular gas is always a subdominant mass component . at lower stellar masses
there is large uncertainty on the fractional contribution of h@xmath96 , as it is not precisely known how well the galactic co emission traces molecular hydrogen mass .
in particular , the @xmath97 conversion factor may vary by about an order of magnitude as we consider less luminous and more metal poor late - type galaxies ( e.g. * ? ? ?
determining the contribution of ionized gas to the total baryonic mass budget of galaxies is much more challenging .
for example @xcite argue that warm ionized hydrogen ( hii ) may amount to about 1/3 of the mass of atomic hydrogen ( hi ) in the disk of the milky way .
if the ratio of ionized - to - neutral hydrogen mass ( @xmath101 ) were fairly constant among galaxies , then the baryonic mass of a galaxy would be given by the expression @xmath102 . if @xmath103 , then the peak value of the @xmath21 - @xmath9 relation ( see fig .
[ fig : am2 ] ) would shift to lower halo mass and the peak value would slightly increase . however , since the precise value and scatter of @xmath101 is not well constrained -and its dependence on galaxy size is not known- we choose not to include the contribution of warm ionized gas in the calculation of @xmath104 .    assessing the contribution of the hot ionized medium ( him ) to the total baryonic mass of a galaxy
is even more challenging .
the coronal him may be the dominant baryonic mass component in galactic halos , especially in massive ellipticals .
determining the mass contribution of the him for less massive galaxies however is observationally challenging . in any case , the tightness of the `` baryonic tully - fisher relation '' when computed just from the stellar and hi mass ( e.g. * ? ? ?
* ; * ? ? ?
* ) implies that the him never dominates the total baryonic mass budget of late - type galaxies , at least within the extent of the galactic hi disk .
let @xmath105 be the cosmic number density of galaxies with baryon mass greater than @xmath104 and let @xmath106 be the cosmic number density of haloes with mass greater than @xmath9 .
the fundamental assumption of the abundance matching method ( @xcite ; @xcite ; also see @xcite for a review ) is that @xmath104 is a monotonically increasing function of @xmath9 . with this assumption
, @xmath107 can be determined by solving the equation    @xmath108    in reality , the baryon content of a halo will depend not only on its mass but also on other parameters , such as its formation history . as a result , a scatter in the distribution of @xmath104 at a given @xmath9 is expected .
neglecting the scatter is , nevertheless , justifiable because the aim of abundance matching is precisely to determine the average value of @xmath104 within a halo of mass @xmath9 .
we evaluate the right hand side of equation  [ eq : am ] , using a halo mass function extracted from one of the cosmological n - body simulations of the horizon project .
the simulation was run with a public version of the gadget code @xcite , and uses @xmath109 particles of mass @xmath110 to simulate the formation and evolution of dm structures in a comoving volume of 100 @xmath111 mpc on a side .
it assumes a cosmology and initial conditions which are consistent with wilkinson microwave anisotropy probe ( wmap ) third year results @xcite , namely @xmath112 , @xmath113 , @xmath114 , and @xmath115 .
the identification of dm haloes and sub - haloes was done with the adaptahop algorithm presented in @xcite .
haloes are identified as groups of particles above a threshold over - density of 80 times the mean density of the universe , which corresponds to a mean overdensity contrast of about 200 .
the identification of subhaloes within haloes is done using the method described in @xcite .
we only keep haloes and sub - haloes with more than 20 particles , i.e. we introduce a minimum halo mass of @xmath116 .
- @xmath9 relation .
the black solid line shows the @xmath8 - @xmath9 relation obtained by matching the stellar mass function of central galaxies to the mass function of haloes excluding subhaloes .
this matching should reproduce the `` true '' relation for central galaxies .
the blue dot - dashed line is the relation obtained by matching the total galaxy stellar mass function , including central and satellite galaxies , to the mass function of haloes excluding subhaloes .
the red dashed line is the relation obtained by matching the total galaxy stellar mass function with the mass function of haloes including subhaloes .
note that , unlike all other figures , this figure uses the stellar mass functions of @xcite , who have separately measured the smf for central and satellite galaxies . ]        an important issue for our analysis is determining whether subhaloes should be included or excluded in the calculation of the hmf used in the abundance matching procedure .
let @xmath117 be the stellar mass of the central galaxy in a halo of mass @xmath9 and let @xmath118 be the stellar mass of the satellite galaxy in a subhalo of mass @xmath9 . if subhaloes are excluded , then we implicitly assume that all galaxies in our samples are central galaxies , that is , @xmath119 . on the other hand ,
if subhaloes are included , we implicitly assume that the @xmath104 - @xmath9 relation is the same for both central and satellite galaxies , that is , @xmath120 .    in order to understand how much this choice could affect our results , we consider the stellar mass functions of @xcite , for which separate distributions for central and satellite galaxies have been presented .
this allows us to calculate the @xmath8 - @xmath9 relation in three ways : firstly , we can match the smf of central galaxies with the halo mass function excluding subhaloes ; this method will give us the correct @xmath8 - @xmath9 for central galaxies , shown by the black solid line in figure  [ censat ] .
then , we consider the total smf of central plus satellite galaxies and match it with the halo mass function excluding subhaloes .
this is equivalent to assuming that all galaxies are central and overestimates the @xmath8 - @xmath9 for central galaxies ( blue dot - dashed line in fig .
[ censat ] ) .
finally , we can match the total smf of central plus satellite galaxies with the total halo mass function including subhaloes .
the result , shown by the red dashed line in figure  [ censat ] , lies below the black solid line because it is effectively a weighted average of the relation for the dominant central galaxy population ( the black solid line ) and the relation for the satellite population , which has lower @xmath8 for a given @xmath9 .
quantitative comparisons of @xmath118 and @xmath117 have received the attention of much recent literature ( see e.g. cattaneo et al .
2012 , submitted , @xcite , @xcite ) .    as we do not distinguish between central and satellite galaxies in the samples used in this work , we shall choose the third method ( the one that corresponds to the red dashed line ) as our best estimator of the @xmath8 - @xmath9 relation and the @xmath104 - @xmath9 relation .
this choice will introduce some systematic bias which is , nonetheless , smaller than the typical uncertainty involved in the determination of the @xmath8 - @xmath9 relation .
we also considered whether abundance matching our hi - selected smf with the hmf excluding subhaloes would give consistent results with our fiducial abundance matching result , obtained by matching our optically - selected smf and the hmf including subhaloes .
physically , this consideration was motivated by the fact that satellite galaxies tend to be redder than central galaxies of the same mass ( e.g. * ? ? ?
* ) , and so hi - selection may be equivalent to the exclusion of satellite galaxies from an observational sample . however , the comparison in fig .
[ censatbis ] of our hi - selected smf with the smf for central galaxies of @xcite shows that this argument may not be valid .
we note that the comparison between the yang et al .
smf for central galaxies and our hi - selected smf is subject to observational systematics , such as the distance assignment scheme or the stellar mass estimator ; for example using hubble distances for the galaxies in our hi - selected sample ( see fig .  [
fig : hubble_vs_flow ] ) would bring the two distributions in fair agreement .
relation ) .
the thick yellow line shows our main result , obtained from abundance matching the stellar mass function of our optically - selected sample with the halo mass function including subhaloes .
the yellow dashed and dash - dotted lines correspond to variations of our main result , obtained by considering the @xcite smf and excluding subhaloes from the hmf respectively , and are shown to illustrate uncertainties .
the magenta dotted , cyan dashed , blue dot - dashed , green solid and red solid lines correspond to the abundance matching results of @xcite , @xcite , @xcite , @xcite and @xcite , respectively .
the big green circles are the results of a stacked weak lensing study of sdss galaxies by @xcite .
all other data points refer to measurements for individual galaxies : the small circles correspond to galaxies with halo mass measurements from weak lensing studies ( @xcite : cyan circles ; @xcite : red circles ; @xcite : blue circles ) .
the star symbols show galaxies in which the halo mass was determined from studies of stellar dynamics ( @xcite : cyan stars ; @xcite : red stars ) .
the triangles show galaxies with halo masses determined from the disc rotation speed ( @xcite : cyan triangles ; @xcite : red triangles ; @xcite : blue triangles ) .
the dotted - dashed horizontal line shows the cosmic baryon fraction @xmath121 . ]    our abundance
matching analysis produces two main results :    1 .
we determine @xmath122 as a function of @xmath9 by comparing the galaxy stellar mass function ( smf ) from our optical sample to the total halo mass function ( including subhaloes ) .
[ item1 ] 2 .
we determine @xmath123 as a function of @xmath9 by comparing the baryonic mass function ( bmf ) from our optically selected sample to the total halo mass function ( including subhaloes ) .
[ item2 ]    the latter relation is the focus of this paper , but we first discuss the former because the results can be compared to an extensive literature of previous studies .
the consistency of our findings with previous work on point [ item1 ] boosts our confidence that our conclusions on point [ item2 ] are robust
.    figure  [ fig : am1 ] shows our result for @xmath122 as a function of @xmath9 , obtained from our optically selected sample ( gold thick solid line ) .
our analysis extends to halo masses as low as @xmath124 , since both our optically - selected and hi - selected samples probe stellar masses down to @xmath125 .
at the high mass end , our relation stops at @xmath126 because our galactic samples are drawn from a relatively small volume , and are not appropriate to measure the abundance of the most massive galaxies and clusters .
we find good agreement with previous estimates from abundance matching obtained in the same manner ( @xcite : cyan dashed line ; @xcite : blue dotted - dashed line ; @xcite : green solid line ; @xcite : red solid line ; @xcite : purple dotted line ) .
the thick yellow dash - dotted and dashed lines are shown to illustrate the systematics introduced by the choice of hmf and smf : the former represents the abundance matching result when our fiducial smf is matched the hmf excluding subhaloes ; the latter is the result of matching the yang et al .
smf with our fiducial hmf , which includes subhaloes .
we also compare our average @xmath8 - @xmath9 relation with values measured for individual galaxies .
the small cyan , red and blue circles correspond to galaxies with measurements of their halo mass @xmath9 from weak lensing @xcite .
the star symbols correspond to galaxies for which @xmath9 was estimated form stellar dynamics @xcite .
the triangles correspond to disc galaxies for which @xmath9 was determined from their rotation speed @xcite .
we remark that , while results for individual galaxies have large scatter , they seem to be systematically lower than any of the relations inferred from abundance matching ( at least over the halo mass range @xmath127 ) .
furthermore , the halo mass for which @xmath128 has a maximum appears to be higher when inferred from measurements of individual galaxies compared to the value derived from abundance matching : in the first case the peak is at @xmath129 , while in the second case it is at @xmath130 .
there is also slight disagreement of all abundance matching results with the results of @xcite ( large green circles ) , who used a stacked weak lensing analysis of over a hundred thousand disk galaxies in sdss separated in three bins of stellar mass .
regardless of the method used , however , there is a clear consensus that @xmath122 is much lower than the universal baryon fraction @xmath121 ( horizontal black dash - dotted line in figure  [ fig : am1 ] ) , at all halo masses .
relation ) .
the gray shaded area shows the results of an abundance matching analysis of our optically - selected sample .
its boundaries correspond to two extreme assumptions for the gas content of galaxies detected optically but not in hi : i ) galaxies that are not detected in hi contain no gas ( lower boundary ) and ii ) galaxies that are not detected in hi contain the largest amount of gas that could have escaped detection from alfalfa ( upper boundary ) .
the thick solid yellow line is the @xmath128-@xmath9 relation ( same as in figure  [ fig : am1 ] ) , and is shown for comparison .
the dotted - dashed line shows the baryon fraction that @xcite predict based on hydrodynamic simulations that include cosmic reionisation . ]
let us now examine the results for the @xmath123 - @xmath9 relation ( fig .
[ fig : am2 ] ) .
the gray shaded area shows the relation derived from our optically - selected sample , matched to a halo mass function that includes both haloes and subhaloes : its upper and lower envelopes correspond to the distribution for @xmath131 and @xmath132 , respectively , as defined in   [ sec : derived ] .
the thick gold solid line represents the relation for the stellar mass ( same as in figure  [ fig : am1 ] ) and has been added for reference .
figure  [ fig : am2 ] puts in evidence the fact that @xmath21 decreases monotonically with decreasing halo mass , despite the fact that atomic gas contributes progressively more to the baryonic mass budget of galaxies .     - @xmath9 relation derived in this work from our optically - selected sample ( same as fig .
[ fig : am2 ] ) . the magenta dotted line and the green solid lines show the @xmath123 - @xmath9 relations that @xcite and @xcite derived from their @xmath122 - @xmath9 relations ( fig .
[ fig : am1 ] ) , using the mean gas - to - stellar mass ratio as a function of stellar mass to account for the gas content of galaxies . the black solid line correspond to the results obtained from our optically - selected sample when we use the same method .
the red solid line represents the result of @xcite , who used separate @xmath133 - @xmath8 relations for blue and red central galaxies .
lastly , the thin black dotted - dashed line is the same as fig . [
fig : am3 ] and shows the baryon fraction that @xcite predict based on hydrodynamic simulations that include cosmic reionisation . ]    in fig .
[ fig : am3 ] we compare our results to those by @xcite and @xcite , who derived their @xmath123 - @xmath9 relations from the equation    @xmath134    in eqn .
[ eqn : mbmh ] the baryonic mass is computed from the stellar mass , using the mean gas - to - stellar mass ratio for galaxies as a function of stellar mass . to enable a cleaner comparison ,
we have made the exercise of re - deriving the @xmath123 - @xmath9 relation from our @xmath122 - @xmath9 relation , using the same procedure followed by @xcite ( as in eq .
[ eqn : mbmh ] ) .
the result is shown by the black solid line in fig .
[ fig : am3 ] .
the main difference between the results obtained by using individual galaxy gas masses ( gray shaded region ) and by adopting a mean gas - to - stellar mass ratio ( black solid line ) seems to be an artificial flattening of the @xmath135 relation at low masses .
this is probably due to the fact that the latter method ignores the large scatter of galactic @xmath136 values from the mean , and leads to the incorrect interpretation that the baryon retention fraction ( @xmath21 ) of low - mass halos asymptotes to some fixed value .
the red solid line in fig .  [ fig : am3 ] corresponds to the result of @xcite , who used a separate mean @xmath136 - @xmath8 relation for blue and red galaxies .
their result shows no signs of flattening at low masses , however the relation only extends down to @xmath137 .
independently of the used method however , all results point to values of @xmath21 that are well below unity , and can not be explained by the effects of cosmic reionization alone ( black dotted - dashed line in fig .
[ fig : am3 ] ; see sec .  [
sec : discussion ] for details ) .    ) for the galaxies in our optically - selected sample that are detected by alfalfa , while the red inverted triangles are maximum hi - fractions ( @xmath138 ) for a representative subset of our optically - selected galaxies that are not detected by alfalfa . ]    the @xmath139 relation that we use to evaluate the right hand side of eqn .
[ eqn : mbmh ] ( black solid line in fig .
[ fig : gsmr ] ) is a power - law fit to the @xmath99 data by @xcite , @xcite , @xcite and @xcite . the best fit relation , plotted as a thin solid black line in fig .
[ fig : gsmr ] , is given by @xmath140 .
the blue contours in fig .
[ fig : gsmr ] represent the distribution of @xmath99 values for the galaxies in our optically - selected sample that are detected by alfalfa .
since the alfalfa survey is a blind hi survey with a fixed integration time per pointing , the alfalfa distribution is expected to lie above the data obtained by pointed observations of optically - selected galaxies ( as in @xcite ) .
the inverted red triangles correspond to the maximum hi - fraction ( @xmath138 ) for a representative subsample of our optically - selected galaxies that lack an alfalfa detection .
note that these upper limits are also systematically higher than the relationship indicated by the power - law fit .
the main result of this paper is the large `` gap '' between the present - day baryon fraction of galaxies in low - mass halos and the cosmic value ( @xmath121 ) , which is present even when the atomic gas content of the galaxies is taken into account .
this result is not contrived , given that atomic gas dominates the baryonic mass budget of galaxies with @xmath141 . moreover ,
the low - mass behavior is in disagreement with previous studies ( e.g. * ? ? ?
* ; * ? ? ?
* ) , who find that the @xmath21 - @xmath9 relation flattens out at low masses , and approaches a roughly constant value of @xmath142 .
these previous results would then require an exceptionally low efficiency of gas - to - stars conversion in low mass systems to explain the observed values of @xmath11 , which decrease monotonically .    in fig .
[ fig : am2 ] we compare our result for the baryon fraction of halos to the predictions at @xmath143 of a cosmological hydrodynamic simulation that includes heating from a photoionizing uv background ( * ? ? ?
* black dotted - dashed curve ) . at the high - mass end ( @xmath144 )
, the low values of @xmath21 can be explained by the fact that infalling gas is shock - heated @xcite , and is kept hot by feedback from active galactic nuclei @xcite .
this picture is supported by considerable observational evidence in the case of x - ray groups and clusters , but the situation within individual galaxies is not so well understood ( see @xcite for a review ) .    at the low mass end
, photoionization heating is expected to become important , since the intergalactic medium is too hot to fall into the shallow potential wells of haloes with @xmath145 , and their baryon fraction is heavily suppressed . however , this process alone can not account for the shape of the @xmath21 - @xmath9 relation at low masses , especially for the onset of a sharp decrease in the @xmath123 ratio at relatively large halo masses ( @xmath146 ) .
additional feedback is therefore needed , usually attributed to the ejection of baryons by stellar - driven winds .
semianalytic models of galaxy formation based on this assumption reproduce a good fit to luminosity functions in the local universe @xcite but the implied outflow rates are enormous . to reproduce the result presented in figure  [ fig : am2 ] ,
the outflow rate in a halo with @xmath147 must be of order a hundred times higher than the star formation rate .
it is difficult at present to reproduce such outflow rates in hydrodynamic simulations .
moreover , observational estimates place the total mass of outflows in normal star - forming galaxies at approximately the same level as the galaxies final stellar mass @xcite . even in the most extreme observational cases ,
the `` mass - loading factor '' ( the ratio of mass loss rate due to outflows over the star formation rate ) is estimated to be @xmath148 @xcite .
therefore , explaining in detail the mechanisms responsible for the very low @xmath21 values found in low - mass galaxies seems to be a fundamental challenge for models of galaxy formation in a @xmath0cdm cosmological context .
recent cosmological hydrodynamic simulations have challenged this statement , by managing to reproduce `` realistic '' galaxies whose properties satisfy a number of observational constraints . for example , the high - resolution eris simulation @xcite has managed to produce a milky - way type object with values of @xmath11 in agreement with those that we see in fig .
[ fig : am2 ] .
more recently , @xcite has managed to reproduce a population of @xmath1491000 simulated galaxies with low stellar - to - halo mass ratios ( @xmath150 at @xmath151 ) , in accordance to observations .
notice , however , that this simulation uses a kinetic rather than thermal wind model , in which the initial wind speed is @xmath152 and the initial mass - loading factor is a factor of four , by construction .
lastly , the work of @xcite has been extended to lower masses by @xcite , who managed to produce a pair of dwarf galaxies ( @xmath153 ) , which obey the observed `` baryonic tully - fisher '' relation and are therefore expected to have the correct baryon - to - halo mass fractions .    while these studies indicate that we may be heading toward a solution of the discrepancy between the observed and the expected baryon content of dark matter halos , many questions ( e.g. the expected outflow rates and re - accretion timescales ) remain unanswered .
therefore , explaining in detail the mechanisms responsible for the very low @xmath21 values found in low - mass galaxies remains an open problem for studies of galaxy formation in a @xmath0cdm universe .
furthermore , our measurement of @xmath154 provides an additional constraint with which cosmological hydrodynamical simulations and semianalytic models will have to confront .
we use optical data from the seventh data release of the sloan digital sky survey ( sdss dr7 ) and 21 cm emission - line data from the arecibo legacy fast alfa ( alfalfa ) survey to measure the `` baryonic mass '' ( defined as @xmath22 ) of galaxies in the local universe , and determine the @xmath143 baryon mass function ( bmf ) .
we use both an hi - selected and an optically - selected sample ( 7195 and 22587 galaxies respectively ) drawn from the same volume , in order to address the effects of sample selection on the derived distributions .
we find that the main difference consists of the optically - selected stellar mass function ( smf ) being systematically larger at high - masses than the hi - selected smf , and find that this difference carries over to the high - mass end of the bmf ( see fig .
[ fig : opt_vs_hi_gsmf ] & [ fig : opt_vs_hi_bmf ] ) .
we combine the obtained mass distributions with the halo mass function in a wmap3 @xmath0cdm cosmology , to obtain average values of @xmath128 and @xmath135 as a function of halo mass ( fig . [
fig : am1 ] & [ fig : am2 ] ) .
our most important result is that low - mass halos seem to have very low galactic baryon fractions compared to the cosmic value ( @xmath2 ) , even when their atomic gas content is taken into account ; for example , the average baryon fraction of halos with @xmath155 is just 2% of the cosmic value ( @xmath156 ) , and displays a monotonically decreasing trend .
this result contrasts with previous indirect measurements of the bmf @xcite , which pointed to an approximately constant value of @xmath157 at the low halo - mass end .
such very low values of @xmath21 are difficult to reconcile with current models of galaxy formation .
photoionization heating in the early universe suppresses the baryonic content of halos only at @xmath158 @xcite , but this mass is more than an order of magnitude smaller than what is required by our result .
therefore , additional feedback mechanisms , such as baryon blowout by supernova explosions , must be present and must be extremely efficient .
it is not yet clear whether hydrodynamic simulations or observational results can accommodate such intense galactic outflows in low mass halos . as a result
, the observed @xmath21 - @xmath9 relation remains difficult to explain , and may represent a challenge to our understanding of galaxy formation and/or the properties of dark matter .
the authors acknowledge the work of the entire alfalfa collaboration team in observing , flagging , and extracting the catalog of galaxies used in this work .
the alfalfa team at cornell is supported by nsf grants ast-0607007 and ast-1107390 to rg and mph and by grants from the brinson foundation .
ep would also like to thank the onassis foundation for partial support .
ac thanks alexie leauthaud for providing many of the data points used to make fig .
[ fig : am1 ] .    funding for the sdss and sdss - ii has been provided by the alfred p. sloan foundation , the participating institutions , the national science foundation , the us department of energy , the nasa , the japanese monbukagakusho , the max planck society and the higher education funding council for england .
the sdss web site is ` http://www.sdss.org ` .
the sdss is managed by the astrophysical research consortium for the participating institutions .