reliable determination of total mass of spiral galaxies is a crucial step in the search for non - baryonic dark matter . in particular , laboratory experiments aimed at detection of dark matter particles ,
could succeed only if there is sufficient amount of dark matter inside the disk of our galaxy in the sun s vicinity .
present models of mass distribution in spiral galaxies assume that in addition to stellar disk , bulge , luminous halo , there is massive halo of non - baryonic dark matter , most often assumed spherical .
it dominates the total galaxy mass and also governs the dynamics of outer parts of the stellar disk .
this last feature makes it possible to test in some cases , if spherically symmetric rather than flattened dark matter distribution is actually present .
recently , we have analyzed mass distribution in several spiral galaxies @xcite .
we find that for some galaxies it can not be spherical at larger radii . since we expect in this case a flattened mass distribution to better approximate the gravitational potential at larger radii than the spherical one , we apply the global disk model .
we then find that mass distribution of luminous matter accounts for rotation curves of the examined galaxies .
there is also increasing number of studies whose authors , using more involved galaxy models than so far , conclude , contrary to earlier findings , that luminous matter accounts for rotation in the internal galactic regions @xcite .
there is also possibility that rotation of galaxies in the outermost regions can be driven by magnetic fields , not by dark matter @xcite .
these results need to be seen properly in the context of non - baryonic dark matter searches .
currently , there is strong observational evidence of non - baryonic dark matter in clusters of galaxies , in particular , from gravitational lensing by clusters and from colliding clusters , as e.g. the bullet cluster @xcite .
also , the cosmological model provides support to the non - baryonic dark matter hypothesis , albeit at a more theoretical level . on the other hand ,
direct searches give null results , except for dama / libra @xcite , however , these results are disputable and have not been confirmed by other searchers .
the consistent picture might emerge if one assumed that clustering scale of non - baryonic dark matter is that of galaxy clusters . at smaller scales ,
non - baryonic dark matter could be dilute .
in particular , its contribution to spiral galaxy masses could be much less than anticipated in spherical halo models .
it can not also be excluded that spiral galaxies may differ from each other in abundance of dark matter
some of them could be reach in dark matter , the other might be devoid of it .    before presenting results for some spiral galaxies obtained with our model based on a global disk approximation ,
we briefly signalize some important issues concerning rotation curve modeling .    the methods for ascertaining mass distribution in spiral galaxies commonly use very simplified , one would even say naive , one- or many - component galaxy models , that assume symmetries of mass distribution ( axial , spherical ) , properties of rotation fields ( concentric circular ) , and composition characteristics ( constant @xmath0 ratios , _ etc _ ) , which at most approximately agree with observations or sometimes even contradict them .
our model of spiral galaxies is also simplified , however , it possess some features which , according to us , gives it the advantage over other models
e.g. it does not assume in advance any m / l ratio or its constancy , but is obtained at output as a resulting m / l ratio profile , which seems more realistic .    in the past , the prevailing opinion was that one needs dark matter inside galaxies to account for their rotation and to stabilize their allegedly cold disks , whereas there is increasing number of papers modeling the same galaxies without dark matter .
this is achieved by considering sufficiently complex and thus also more `` flexible '' models , for example , those assuming wide , not thin , stellar disks , ect .
there is also too much freedom in choosing galaxy models , and there is no general agreement which models are closer to reality .
but , why to assume circular orbits , axial or spherical symmetry , constant mass - to - light ratios when this is not the case for real galaxies ? as is pointed out by @xcite , most of spirals are not differentially rotating axisymmetric disks with only tiny disturbances due to spiral arms
. there are large scale deviations from the symmetry , often with very large amplitudes or even with warped hi layers @xcite .
the predictions of the simple galaxy models might further be modified if one preferred to reconstruct whole doppler images rather than rotation curves , the more the latter are only some interpretations of doppler images obtained not without additional assumptions .
going even further , what would happen if , instead of using the simplified models , one used a sort of templates for inferring matter distribution in real galaxies based on their doppler images and the observed qualitative characteristics like the number of arms , the presence of a bar _ etc _ ?
such templates , relating mass distribution and its qualitative shape features to simulated doppler images could be catalogued by elaborating a large number of galaxy - like virialized stationary systems from n - body simulations .
would nt this approach be the right step toward realistic modeling of mass distribution in galaxies ? to what extent could this approach change our understanding of spiral galaxies and current opinions about their composition ?
the question how much dark matter is present in galaxy s interiors remains still open @xcite , @xcite @xcite point out that the better are the data used for constraining the properties of galaxy halos and the better is the quality of simulations , there is more difficulties with dark matter models on galaxy scales .
also @xcite says that the case for cdm , having < < a number of stubborn problems > > , is weak judged on the data from galactic scales only .
this is especially evident from the still growing number of examples of morphologically various galaxies whose rotation curves can be accounted for without dark matter hypothesis by using non - canonical , more realistic approach to mass modeling .
for example , even though still assuming the simplifying unrealistic axi - symmetry and that the radial mass profile should follow that of luminosity with constant @xmath0 ratio , @xcite were able , under the strict maximum ( thick ) disk hypothesis , that is , not including a dark halo , to reproduce  with i - band @xmath0 ratio of @xmath1 and consistently with stellar population synthesis models  the overall structure of the optical rotation curves for most of galaxies in their sample of 74 spiral galaxies with various surface brightness profiles and rotation curve morphologies .
they link the poorest fits with non - axisymmetric features , like bars or strong spiral arms , that influence the brightness profiles , the determination of inclination , major - axis position angle , _
etc_. this situation , however , might be cured by considering more realistic models .    that the interiors of spiral galaxies can not be dominated by dark matter is also suggested by dynamical arguments . as noted by @xcite ,
if the average projected surface mass density of dark matter halo was greater than that of the optical disk , the common instabilities would be suppressed @xcite together with lopsided modes @xcite appearing in many disk galaxies .
in addition , as they point out , there are models that reproduce features of rotation curves within the optical disks on different scales and that the models rarely require dark matter halos @xcite .    these interesting observations of absence of dark matter inside galaxies are best expressed by @xcite .
they speculate , having succeeded in fitting some galaxies without dark matter , whether occam s razor argument would not be applicable .
why to introduce something that is not needed ?
the authors also form less authoritative statements about dark matter presence in galaxies than others .
what they merely say is that dark matter presence is only statistically significant .
the presence of a dark matter component that dominates the outer mass distribution in spiral galaxies , seems evident from their optical rotation curves which do not fall off sufficiently fast close to their optical edge @xcite and from their extended hi rotation curves @xcite .
however , in the context of what has been said above , one can also conjecture that the halo would disappear if the galaxy models were more realistic  allowing for changing mass - to - light ratios , assuming non - axisymmetric mass features ( such as spiral arms and bars ) , if the assumed extinction models was altered_etc_. such models should be also taking into account that the velocity field in galaxies is much more complicated than circular and that the published rotation curves are obtained by reducing a two - dimensional doppler image of a galaxy into a single one - dimensional curve ( such an image by itself carries already a reduced , partial information about rotation ) , and that this reduction is performed not without additional assumptions concerning the interpretation , averaging methods , _ etc _ , of the doppler image ( for example , an averaging method which produces rotation curve by calculating pairs @xmath2 $ ] scanning a doppler image , where the averaging at a given @xmath3 is carried out over a window of size @xmath4 and weighted by , say , a function of local intensity of the image within the window , is presumably good for circular orbits of test bodies in a spherically symmetric gravitational potential , but not necessarily for a self - gravitating system which may be flattened ) .    as concerns dark matter in outer galaxy regions , there is also possibility ,
never taken into account in modeling of flattened galaxies , that rotation of slightly ionized gas at large radii may be driven also by magnetic fields whose strength may be comparable with that of gravitation of whole galaxy at these radii @xcite , helping to lower ( or even eliminate ) the amount of dark matter needed in the outskirts of galaxies where rotation curves are usually flat @xcite .
it is true that realistic modeling of galaxy rotation is very challenging and might be very difficult to achieve , but , at the same time it is not clear whether the naive models used so far , lead to trustworthy results  has the opposite ever been rigorously proved ?      in modeling of galactic optical discs , it is assumed they are highly flattened , axisymmetric and cold .
there is a common opinion linked with stability argument of @xcite , that for such disks to be stable within the optical radius , one necessarily needs the galaxies to be immersed within massive dark matter haloes .
this conclusion comes about as follows . as observed in numerical simulations , when random stellar motions are small compared to the streaming rotational motion in flattened galaxies , without such a halo a bar would form in one rotation period , and only afterward axial symmetry would be regained . during this stage ,
approaching stationary state , the disk would heat up and could not remain cold ( this would not agree with the assumption that disks are cold ) .
therefore , cold disks must be immersed in ( usually spherical ) halos of sufficiently large mass .
there are various criteria known for deciding whether a virialized rotating flattened system could be stable .
@xcite conclude that a spherical halo with halo - to - disk mass ratio of @xmath5 to @xmath6 and initial value of total rotational kinetic energy to total gravitational energy ratio of @xmath7 renders the disk stable ( of course , the halo mass beyond the disk s outer radius has no effect and should not be taken into account ) .
their conclusions are supported by many computer simulations and also by linear stability analysis carried out by toomre for mestel disk @xcite .
note , however , that for nonlinear differential equations ( and hydrodynamical description of galaxies leads to such equations ) it is not a priori obvious that linear approximation appropriately describes solutions , the more that the gravitational potential also evolves along with matter , thus also the instability modes and the associated spectrum of eigenvalues .
@xcite criticized the criterion , as well as various other criteria , to be of little use .
nevertheless , he points out , they all lead to similar demands on mass distribution .
as shown in @xcite , based on @xcite models , one comes to the conclusion that , for stability , total spherical mass interior to the half - mass radius of the disc should be more than twice the disc mass in the same region .
@xcite and @xcite find similarly , using their numerical models , that for stability of differentially rotating disks the ratio of spheroidal mass to that of the disk interior to the outer radius of the disk , should be approximately @xmath8 , while for rigid disks the mass could be smaller @xcite .
this spheroidal - to - disk mass criterion is much more handy to use than the ostriker - peebles  like criteria in modeling galaxies , since mass functions of different galaxy components are much simpler to find .
it should be noted , however , that the assumption of cold disks is based on the belief that the local observation of small dispersion of velocities in solar neighborhood in our own galaxy can be extrapolated onto whole galaxy and even onto all other flattened galaxies . even if the assumption was true , there is an argument due to @xcite that the halo is not so efficient a stabilizer as a bulge alone could be .
a massive halo would change rotation field of the bar in a way not measured in real galaxies @xcite . as for modest halo
, it would be difficult to decide about its presence since its effect would be small @xcite .
it is also known that random motions of stars of the disc near the center have stabilizing property as well as a bulge or a thick disk .
neglecting this effect could result in overestimating the bulge+halo mass required for global stability .
@xcite surmises that most galaxies might have only moderate fractions of mass in the bulge+halo . as concerns the internal galaxy regions comprising the optical disk
, the amount of dark matter is very poorly known @xcite . on the other hand , contrary to the common opinion that only spheroidal halos can help to stabilize cold galactic disks ,
bulges @xcite or thick disks @xcite are now considered to be more efficient stabilizers @xcite .
what s more , the presence of dark matter in large amounts , dominating the disk mass , would simply suppress the instabilities required to form the observed structures such as bars and spiral arms @xcite .    because spherical halo mass outside the optical radius has no effect on disk stability , the facts mentioned above , reverse the disk stability argument by @xcite , frequently used as a support for the need of a massive dark matter halo within the optical radius @xcite . here
, it is worth to add that @xcite construct a model disk galaxy with an almost flat rotation curve which is stable without presence of dark matter .      in modeling of mass distribution in galaxies ,
provided one has indeed a model that properly describes a given galaxy , one in principle should check whether the found mass distribution is stable . if disks were indeed cold , then for stability the spherical mass component should be sufficiently large .
it seems however , that for estimating of the mass distribution in galaxies , the stability issue is not so important for astronomers , the more that the galaxy models used in this respect are very simple , one could say even naive .
majority of them assume idealized substructures such as disks and bulges with simple parametric profiles , and assume symmetries that are not observed ( spherical and axial ) .
for example , the spherical mass to disk mass ratio in model used by @xcite(things ) is much less than 2 , _ e.g.,_for galaxies ngc 4736 , ngc 3031 .
despite the fact that a dark matter halo was included , the ratio would suggest that the found mass distribution is instable as it does not meet the stability criteria mentioned previously , that is , the disk - to - halo mass is too large .
there is a big but , however .
all mass models for flattened galaxies are merely some idealizations that do not take into account different processes that make real galaxies stable .
the role of the models is rather to roughly estimate the gravitational potential from rotation curves assuming circular orbits and axial symmetry , and so on , and then to derive the corresponding mass distribution which one believes to some degree approximate the actual amount of matter in real , and otherwise stable galaxies .
for example , maximal disk models give `` effective '' surface density which should be treated rather as a substitution for different galaxy components both those flattened and spheroidal , described altogether by a single function of radius ; only after a suitable de - projection one would infer about various galaxy components , thus it would be at least exaggeration to apply to such models the stability criteria which says that axisymmetric cool disks are instable .    of course , in real galaxies , masses are neither distributed in infinitely thin disks nor in ideally spherically symmetric structures , orbits are not circular , rotation field has a very complicated structure , and mass distribution does not have axial symmetry
axial symmetry is clearly broken by bars and spiral structure
. thus it would be at least inappropriate to qualify or disqualify different models on the basis of stability arguments even if they do not contain a halo of dark matter , since the models , by construction , do not describe real galaxies and ignore many subtleties .
the models are only slightly more accurate than would be a dimensional analysis , which states that the amount of electromagnetic energy radiated by a galaxy should be given by @xmath9 , where @xmath10 is some characteristic velocity derived from doppler image , @xmath11 is its radial size estimated based on the spatial range of the luminous stuff , @xmath12 is some mass - to - light ratio derived on population analysis ground , and @xmath13 is some unknown dimensionless factor related to the particular geometry of mass distribution and of rotation field .
all parameters apart from @xmath13 can be estimated from measurements .
the models only better approximate the parameter @xmath13 by giving its concrete model , more or less accurate .
for example , the estimation should be better when for an elliptical galaxy one uses a spherical rather than a disk - like model , or , for a spiral galaxy , a disk model rather than spherical .
the approximation should improve when one distinguishes spherical and disk - like subcomponents , further , one can add a width to the disk , and so on .
but it would be hard to treat these models seriously when dynamics and stability of a galaxy is concerned ( especially when maximal disk model is used which treats on equal footing flat and spheroidal internal mass components ) .    without realistic modeling such difficult problems as structural stability , can not be answered . by realistic modeling
we understand finding , for a given galaxy , such a spatial rotation field and mass distribution for an n - body system  a computer model of this galaxy  which would be able to reconstruct at least the observed doppler image of the galaxy .
reconstruction of a published rotation curve only , is not sufficient , because different groups obtain for the same galaxy different rotation curves depending on the details of their procedures applied for processing doppler images .
such images are also often nonsymmetric , unlike simple galaxy models , and lot of information about galaxy kinematics they convey is missing by the fact that the images are by construction a sort of projection of a three dimensional rotational field onto the observation plane .      @xmath0 ratio is important for studying the formation and evolution of galaxies , thus already from this standpoint it is rather strange to assume it is constant from the start , while it would be more desirable to determine it as a function of radius .
unfortunately , the theoretical @xmath0 for an individual galaxy is not well known even when assumed constant .
its derived value depends strongly on the assumed mass function for the stellar population , on the star formation history , on metallicity and depends also on extinction .
it is also affected by the assumed internal extinction law , by the errors in luminosity measurements , and by the distance dimming @xcite .
this ratio is generally assumed constant separately for the bulge and for the disk .
accordingly , one tacitly assumes that the contribution of bulge and that of disk to the rotation curve should be controlled only by the light distribution via a single number found by best fits as if brightness profiles were to trace the overall features of rotation curves .
however , the assumption of constant @xmath0 ratio may not be correct , since it requires that extinction and population gradients across the luminous disk should be small @xcite .
constant @xmath0 ratio seems also improbable since it would imply that composition of a galaxy is homogenous .
overestimating the importance of constant @xmath0 may change considerably and in uncontrolled way the contribution of bulge and that of disk to the overall rotation , the more that optical rotation curves are not featureless , while luminosity is usually fitted by simple laws ( exponential , _ etc _ ) .
in particular , this may lead to errors in reconstruction of the relation between mass density and rotation of the disk component , since rotation of a disk is a nonlocal functional of the density profile , sensitive to the detailed local structure of the profile .
this in turn may lead to wrong estimates of mass profiles of other mass components . to give only an example
, many spirals simply do not have exponential disks or any other given by a simple analytical law , and assuming otherwise , results in large discrepancies in disk scale lengths published by different authors for the same galaxies @xcite .
closely related to the non - uniqueness of the estimates of mass profiles is the so called disk - halo degeneracy  the uncertainty in the relative contributions of different galaxy components to the overall rotation ; for unambiguous disk - halo - bulge decomposition one requires additional constraints on mass distribution , in particular , such as the mass - to - light ratios or the assumed mass profiles and shapes of dark haloes . in finding such constraints
may also be helpful measurements of gravitational lensing in spiral galaxies @xcite .
usual rotation - curve fitting methods have too many free parameters and are not unique , require additional assumptions , sometimes quite arbitrary . for example , @xcite present disk - halo decompositions with several different assumptions for stellar @xmath0 ratio for the same galaxy , or assume that only the minimum disk is present and that the observed rotation curve is due entirely to dark matter , or assume that stellar @xmath0 ratio is zero ( ! ) and take into account contribution of the atomic gas ( h i and he ) .
in contrast , our approach does not have the arbitrariness problem , since we predict , not assume mass to light ratios .    as pointed out by @xcite , there is only limited evidence for constant @xmath0 .
what s more , the results of @xcite for m33 suggest that this assumption is not valid at all , @xmath0 ratio of the galaxy s disk increases approximately 5 times over 6 scale lengths of the inner disk .
their findings also suggest that there is very little dark matter in m33 s central regions .
as we have seen above , the case for dark matter inside galaxies is very weak and the amount of ascertained mass of dark matter is strongly model - dependent . it should be stressed that the models are constructed in the same framework of newtonian gravitation . there is increasing number of papers where more and more complex models are considered compared with previous simple models of galaxies .
they account for the observed rotation of internal regions of flattened galaxies without dark matter , and some of them , based on maximal disk hypothesis , account for global rotation also without dark matter .
the use of models with extended disk finds additional support from the fact that dark haloes may be very flattened rather than spherical . for flattened systems
, however , gravitation gets more complicated and new effects difficult to tackle with , such as influence of external masses on internal orbits , must be taken into account and which would be absent if external masses were distributed spherically symmetric .
also stellar mass - to - light ratio need not be constant .
assuming otherwise , galaxy mass models become very stiff leading , for example , to disk - halo conspiracy .
as we saw , this also leads to discrepancies between predictions .    as concerns the presence of dark matter in spiral galaxies , we therefore prefer a more cautious approach than just to assume the presence of the ( unobserved ) dark matter halo from the beginning .
we propose to admit dark matter , only if models that composed of baryonic matter distributed in a flattened disk ( of stars and gas ) and of more spherical bulge and stellar halo , fail to properly account for the dynamics of the disk measured through rotation curves .
such an approach , we believe , could significantly reduce the undesirable model  dependence of the amount of cdm in spiral galaxies @xcite .
most models of galaxy formation predict that @xmath0 should be a declining function of distance , e.g. the inner regions of a galaxy form first and contain older , higher @xmath0 ratio star populations @xcite . as follows from the previous section it is not correct to assume that the ratio is constant throughout galaxies .
this is the reason why we prefer to determine mass distribution in disk based on measurements other than brightness profiles first , and only then determine mass - to - light ratios as a result of comparison of the obtained mass profile with the measured luminosity profile .
in several galaxies we have studied so far using our approach we observe that this ratio indeed decreases while density approaches that of hydrogen and helium known from measurements .
this also suggest that dark matter is not needed for some galaxies even in the outermost regions .
conceptually simple , direct method for determining the gross mass distribution in galaxies , is rotation curve inversion .
maximal disk model is an example .
provided rotation curve is known globally , the corresponding mass distribution can be found via single integral functional .
the method is independent of any conventional assumptions concerning decomposition of a galaxy into various idealized subsystems with their unknown ( usually correlated with each other ) parameters determinable from best fits .
unfortunately , direct rotation curve inversion is not unambiguous if mass distribution does not have special symmetries , e.g. spherical symmetry .
various ambiguities pertinent to gravitation of flattened systems , easily hidden by the use of a particular model for the dark mass @xcite , may influence the actual importance of different galaxy components .
this fact seems not much appreciated . in addition , there is also ambiguity in mass - to - light ratios within optical disks , since extended rotation curves alone give no indication of where the luminous disk could end @xcite .
consequently , as noted in @xcite , there are mass models which fit rotation curves within the optical disk both with and without significant dark component @xcite
. what s more , rotation of a galaxy considered dark matter dominated in one model can be devoid of dark matter in another model @xcite .    for a flattened mass distribution ,
local density is a nonlocal functional of the rotation field . and vice versa , velocity on a given orbit is affected both by gravitation of masses external and internal to that orbit .
for that example , it would be impossible to reconstruct dynamical mass function from rotation of luminous matter moving in a potential well of unseen external very flattened halo of dark matter .
such dark haloes considerably flattened toward the stellar plane , resembling a disk more closely than a sphere @xcite , with axis ratios between 0.1 and 0.3 where suggested by @xcite and @xcite .
@xcite analyzed errors resulting from extrapolation of the rotation curve beyond the last measured point and gave a criterion for their estimation in global disk model .
these errors are the fundamental obstacle in using disk models .
the reconstructed mass density is viable only out to a distance of @xmath14 of the radius of the last measurement point .
the less serious is sensitiveness of this method to noise criticized in @xcite .
however , as @xcite illustrates , this leads only to @xmath15 uncertainty .
therefore one needs a method to minimize such errors .
once one accepts that very flattened haloes of unseen matter can exist then one may lose any predictability .
for example , since dark haloes are not seen , one may assume anything about them , say , that outside the last measurement point there is present a flattened dark matter ring of a given arbitrarily large mass and given shape , encircling the galaxy .
to balance its influence on the luminous matter , one needs correspondingly large amount of mass inside galaxy .
this in turn would made impossible even to determine the mass - to - light ratio of the luminous matter since the ratio would depend on the unknown mass of the flattened dark matter halo .
there is a simple criterion for deciding , whether the mass distribution at larger radii in some galaxy is flattened rather than spherical . in the presence of spherical dark matter halo , the non - spherical component of mass distribution should be negligible .
then , only the radial component of the gravitational force is important . in this case keplerian mass function @xmath16 defined for
rotation curve @xmath17 must be a nondecreasing function of radius @xmath18}\geqslant0.\ ] ] in deriving this inequality , we assumed circular orbits of matter in the vicinity of the galactic plane , which is customary in modeling of rotation curves .
for spherical systems , the keplerian mass function is identical with the true mass function .
thus , if sphericity condition ( [ eq : sphcond ] ) is not satisfied at larger radii , the true mass distribution can not be spherical , and spherical cdm halo is thereby excluded ( as far as the usual assumptions in galaxy modeling are employed ) , or other forms of dark matter are required forming flattened structures .    we can use the global disk model especially for galaxies with rotation curves breaking the sphericity condition at larger radii . note , that the distribution of internal masses is not much important for the overall gravitation at larger radii ( contribution from higher gravitational multipoles of internal masses decreases quickly with radius ) , therefore , the disk model can be used also for the inner parts , in spite of the fact that the center may be rather spheroidal since , after de - projection of a part of the disk density onto a spheroidal component , this does not change the central mass significantly  this is an approximation , a standard procedure in the maximal disk model .
since rotation curves of disks satisfying the sphericity condition are also possible @xcite , one can try to apply the disk model also to spiral galaxies of which rotation curves do not violate this condition .
the disk model assumes that whole matter in a galaxy is distributed in the vicinity of the galactic equatorial plane .
the model also assumes that orbits of stars are circular and that dispersion of velocities is negligible .
the other simplifying assumption is that whole system is axisymmetric .
the relation between surface mass density @xmath19 and velocity is then given by @xmath20}{\mathrm{d}\chi}\right)\nonumber , \end{aligned}\ ] ] where @xmath21 and @xmath22 are elliptic functions of the first and second kind @xcite .
this integral is understood in the principal value sense ( @xmath23 ) , since both summands are divergent , thus the integral must be skilfully calculated .
in practice , however , it suffices to assume that matter is located close to the galactic plane only in the galaxy outskirts where the external orbits are stabilized by gravitational potential of the galactic interior ( this contribution to gravitational potential is almost spherical at large radii ) , while in the vicinity of galactic interior the density can be treated as an effective density understood rather as a column density of masses projected onto the galactic plane  a substitution for otherwise spherical and stable bulges _
etc_in internal regions . interpreted in this way
, this approximated cold disc model of an otherwise stable spiral galaxy also avoids , or at least weakens , the stability counterarguments often raised against global disk models ( see also the discussion in section [ sec : stabarg ] ) .
it is seen form equation ( [ eq:6 ] ) that for a flattened mass distribution of which the disk model is a limiting approximation , rotational velocity of matter on a given orbit is dependent both of masses distributed interior and exterior to that orbit .
the inverse relation is a functional of whole rotation curve , and it reads @xcite @xmath24.\end{aligned}\ ] ] in the literature of this subject , for some mysterious reason , there is always used an ( equivalent ) expression containing derivatives of velocities which are hardly measurable with satisfactory accuracy , e.g. @xcite .
the crucial point now is that a rotation curve is never known globally , thus the mass distribution even in the region where rotation is known , can not be determined from rotation only .
this is the manifestation of `` non - locality '' of gravity of flattened systems  one can not determine their mass functions from the observed rotation and , therefore , additional constraints on this mass distribution is required .
this is particularly important in the regions where measurements of rotation end , while the analogous error is negligible in the internal galaxy regions .
this property is fatal for reconstruction of mass distribution inside a galaxy when outside the luminous part of this galaxy there is present an external ring - like halo composed of dark matter ; if it were spherical there would be no problem at all .    having realized the unpleasant properties of flattened galaxies we devised in @xcite a method which may help to minimize this uncertainty or remove it .
one does not have to assume anything beyond the last measured point , but simply use the available data of mass distribution as a lower bound for the amount of matter in regions where rotation for some reasons could not be determined , but the amount of gas is still measurable . our goal is to examine how good is our model in reproducing the mass distribution in a sample of spiral galaxies .
the thin disk model s role is to find such a global @xmath19 , that the corresponding velocity of rotation @xmath25 defined by integral ( [ eq:6 ] ) , would best overlap with the observed rotation curve @xmath17 .
this task , however , is not unique , since for every @xmath26 , there is an infinite set of functions @xmath19 that give rise to the same @xmath25 for @xmath27 , even though the @xmath19 s may significantly differ from one another for @xmath28 and have different total masses . to minimize this non - uniqueness
, we look for mass distribution in flattened galaxies by iterations .
a particular realization of such iteration proposed in @xcite assumes that mass distribution of gas in a given galaxy is known for radii greater than the range of the rotation curve .
we bring to the reader s attention the fact , that mass distribution found in @xcite for galaxy ngc4736 using this method , is consistent with all measurements , and this was possible without dark matter .
this result is interesting , since another model of this galaxy @xcite predicts almost 70% abundance of dark matter and does not well explain its rotation at large radii .
it is therefore natural to examine whether this galaxy is exceptional or , maybe , if there are other spiral galaxies with lower abundance of dark matter .
below , we present results for several galaxies obtained with the use of global disk model .
we applied the iteration method analogous to that introduced by @xcite .
we are aware this is not a representative sample of galaxies .
however , we think , that even when rotation of only a single galaxy could be explained without dark matter , this would be interesting as such and worthy to be presented .    initially , we planned to examine only galaxies with rotation curves breaking the condition of spherical mass distribution , because such galaxies are presumably flattened and , therefore , natural candidates to which global disk model is applicable .
however , in order to see whether the model produces viable results for other galaxies , we included also galaxies ngc 6946 and ugc 6446 for which we also found good and complete data .    as concerns the error analysis , our method by construction exactly reproduces rotation curves , thus also fits rotation measurements within error bars .
this is the reason why we show rotation curves without error bars .
results obtained in the global disk model for galaxy ngc 4736 at @xmath29 . *
a ) * rotation curve : things measurements _ [ solid circles ] _ , the model _ [ solid line ] _ , and comparison with high resolution rotation curve taken from @xcite .
* b ) * model global surface mass density _ [ solid line ] _ , surface mass density of hi+he _ [ open circles ] _ ; and surface brightness : b - band _ [ solid squares ] _ , v - band _ [ solid triangles ] _ , i - band _ _ [ solid circles]__@xcite ; * c ) * mass - to - light ratio profile without inclusion of hi+he . ]    here , we once again consider galaxy ngc 4736 , which earlier was studied in @xcite .
this time , however , we use another , newer rotation curve published by @xcite . in figure
[ fig:4736 ] , this new curve is compared with that obtained by @xcite we used last time . both the curves are shown assuming distance @xmath30 different from the previous @xmath31 .
although both curves do not overlap we come to comparable conclusions concerning the amount of non - baryonic dark matter .
ngc 7436 still remains a galaxy in which the halo is not necessarily required .
surface density obtained in the global disk model smoothly approaches that of gas in the outskirts of the galaxy , we obtain also low @xmath0 ratios ( @xmath32 for i - band , @xmath33 for v - band , and @xmath34 for b - band ) . in addition , the local @xmath0 ratio decreases with radius for the largest radii .
therefore , despite the new rotation curve does not break sphericity condition so strongly unlike that by @xcite , we still consider this galaxy as containing no dark matter or that abundance of dark matter is at most very low .
_ upper : _ comparison of the ngc 7793 rotation curves obtained in things @xcite , @xcite and @xcite . symbols stand for measurements and lines for model fits .
_ middle : _ surface mass densities obtained in global disk model compared with surface mass densities of neutral hydrogen and helium and with surface brightness in the b - band .
_ bottom : _ mass - to - light ratios derived based on mass densities obtained for different rotation curves . ]
we take rotation curves published by a ) @xcite , b ) @xcite and c ) three rotation curves published by @xcite .
the latter three curves were merged such to obtain a single rotation curve . in effect
we have three distinct rotation curves for the same galaxy .
as seen in figure [ fig:7793 ] the curves are similar which shows rotation measurements for this galaxy are trustworthy .
the corresponding masses are also close to each other , @xmath35 .
surface mass density from the global disk model smoothly converges to that of gas at the outskirts of the galaxy for all three cases .
we obtain also low values for @xmath0 ratios ( @xmath36 in b - band ) .
we stress that the attempt of @xcite to account for the rotation with spherical cdm halo , gives very poor results . by giving up the halo , we obtain surface mass density in agreement with the amount of gas , and low @xmath0 ratios . therefore , we think it is reasonable to assume that the galaxy does not have any dark matter halo .       _
upper : _ comparison of the ngc 1365 rotation curves obtained in @xcite and @xcite .
symbols stand for measurements and lines for model fits . _
middle : _ surface mass densities obtained in global disk model compared with a 2nd order exponential decay fit and with surface mass densities of neutral hydrogen and helium and surface brightness in the b - band . _
bottom : _ mass - to - light ratios derived based on mass densities obtained different rotation curves . ]    now we consider two rotation curves @xcite and @xcite the latter suitably smoothed out in a way to conform with data points and error bars .
we have noticed that surface mass density corresponding to the rotation curve has approximately exponential falloff ( see figure [ fig:1365 ] ) .
therefore , for a while one may assume that surface mass density is a superposition of two exponential profiles , in order to check whether such mass distribution accounts for the observed rotation .
we find , this is indeed the case , however the mass density becomes larger than that of gas ( although quickly falling off ) . on the other hand ,
the surface density obtained with our method converges smoothly to that of gas already in a place where galaxy is still bright in b - band .
this can be interpreted as that apart from hi and he there must be present other luminous stuff . in general , this is very interesting , that for this galaxy already the doubly - exponential disk alone suffices to account for the rotation curve of @xcite .
irrespectively of which rotation curve is considered , the obtained mass - to - light ratio is low ( @xmath37 to @xmath38 in b - band depending on which rotation curve and method is used ) .
this ratio falls - off with increasing radius in every case .
as so , galaxy ngc 1365 can be modeled without non - baryonic dark matter .       _
upper : _ comparison of the ngc 6946 rotation curves obtained in @xcite , @xcite and things @xcite .
symbols stand for measurements and lines for model fits . _
middle : _ surface mass densities obtained in global disk model compared with surface mass densities of neutral hydrogen and helium and surface brightness in the b - band . _
bottom : _ mass - to - light ratios derived based on mass densities obtained different rotation curves . ]    we take three rotation curves published by @xcite , by @xcite and @xcite ( the latter were smoothed out ) .
as is seen , the published rotation data differ significantly from each other , particularly rotation curve of @xcite exceeds much from the other .
by using global disk model we obtained surface densities which smoothly converges to that of gas .
mass - to - light ratio is low based on rotation by sanders ( @xmath39 in b - band ) , but it becomes twice greater for rotation curve of things ( @xmath40 ) @xmath0 ratio for things grows slightly with radius , for the other two curves it is approximately constant . summing up , based on available rotation data it is difficult to uniquely decide about presence of non - baryonic dark matter in this galaxy .
the result depends of which rotation curve is used .
using our modeling method we can explain all these rotation curves .
it is worth of pointing out that the obtained surface densities smoothly converge to that of gas .
however , based on things rotation curve , we obtain @xmath0 ratio much greater than for ngc 4736 , ngc 1365 , ngc 7793 , similarly , for the other two rotation curves .       _ upper : _ comparison of the ngc 6446 rotation curve obtained in @xcite .
symbols stand for measurements and lines for model fits . _
middle : _ surface mass density obtained in global disk model compared with that of obtained in spherical mass model , with that of neutral hydrogen and helium and with surface brightness in the b - band .
_ bottom : _ mass - to - light ratio derived based on mass density obtained different rotation curves . ]
the @xmath0 profile of galaxy ugc 6446 grows with radius attaining values much higher than for the other examined galaxies , even after subtraction of hi and he .
such a singular behavior of the @xmath0 profile is expected at the edge of luminous disk , where the luminosity per unit of surface tends to zero , @xmath41 , while surface mass density remains finite .
then , the local @xmath0 ratio becomes singular , @xmath42 .
after subtraction of mass contribution of gas , the increasing @xmath0 ratio indicates , that the stars at the outskirts of luminous disk , contain possible admixture of died stars , `` black '' dwarfs , neutron stars , or stars of too low a mass for nuclear burning ( so called brown dwarfs ) .
it is important to note , that the mass density of such a died star component needs only to be residual , such that total mass density is infinitesimally bigger than the mass density of shining stars , @xmath43 for the @xmath0 ratio significantly increases at the edge : @xmath44 when both @xmath45 and @xmath41 .
this very small mass density , @xmath46 , of died stars does not show up anywhere except at the very edge where it is strongly enhanced by small @xmath47 .
contribution of died stars to the mass of galaxy is thus negligible .
the following thought experiment provides a strong argument against the existence of a spherical halo around galaxy ugc 6446 .
suppose for a while that this galaxy is dominated by such a halo .
then , the mass distribution could be considered spherical with good accuracy , thus the corresponding mass function could be calculated from @xmath48 ( and this can be done , since the rotation curve of this galaxy satisfies the sphericity condition ) .
then , the corresponding column mass density @xmath49 can be found by projecting the volume density onto the galactic plane by using the formula @xmath50 } } , & \nonumber\end{aligned}\ ] ] where @xmath26 is the distance of the outermost point on the rotation curve .
note , that our formula is very advantageous from the practical point of view , contrary to the ( equivalent ) handbook formula , as it does not require derivatives of velocities , which are measured with very poor accuracy .
we compare the column density @xmath49 with total surface mass density obtained from the disk model and with its component surface mass density of gas ( hi+he ) .
the result is shown in figure [ fig:6446 ] .
we see that close to the `` edge '' , @xmath49 is lower than the amount of gas !
we stress that the two surface mass densities seen in figure [ fig:6446 ] correspond to two extreme idealizations  in the first the galaxy is approximated by a thin disc , in the other , by a spherically symmetric mass distribution .
the realistic mass distribution must be somewhere in between , but , as is seen , closer to that in disk model . in view of these observations ,
we classify ugc 6446 as a galaxy without massive spheroidal dark halo .
we considered 5 galaxies , among which one was reexamined with a new rotation curve .
our aim was to answer the question if the galaxies could be satisfactorily modeled without non - baryonic dark matter . in this
regard we check whether the following conditions are met : the found mass distribution should account for the measured rotation curve , it should approximately converge toward that observed for gas ( hi + he ) in the outer parts of galaxies ( as this is the lower bound for the amount of luminous matter ) , the mass - to - light ratio should be low and its profile should be a declining function of radius close to the galaxy `` edge '' .
we find that these conditions are satisfied by our results for galaxies ngc 4736 , ngc 7793 and ngc 1365 for all rotation curves .
it is reasonable to state that these are galaxies with flattened , disk - like mass distributions , for which the disk model approximation works quite well .
the case of galaxy ngc 6946 is not so clear .
there are significant discrepancies between published rotation curves , and consequently , in the obtained @xmath0 ratios .
galaxy ugc 6446 satisfies some of our criteria without non - baryonic dark matter , however its @xmath0 ratio is highest and , in addition , it grows dramatically toward galaxy `` edge '' . in another test we assumed that the galaxy is dominated by a spherical dark matter halo and , using a spherical model , obtained the resulting bulk mass density .
this density , accounting for the galaxy rotation , is in this case lower than that of gas , which can be considered as a strong argument for absence of spherical non - baryonic dark halo .
this argument does not exclude the presence of a non - spherical halo .
however , the mass needed in the galaxy outskirts to account for the rotation in the disk model is practically that of gas
. it would be therefore worthwhile to consider other than cdm causes of the @xmath0 ratio growth .
we have also signalized the controversial issue , to what extent the very simplified galaxy models like the ones commonly used by astronomers or like that used in the current paper , can be treated as models of realistic galaxies . in particular , whether stability arguments can be applied to favourize or reject a particular galaxy model , since the role of these models , which serve only as an idealized approximation incorporating many assumptions , is to estimate different global characteristics of galaxies like total mass , the amount of gas , global mass - to - light ratio , dark - to - luminous mass ratio _ etc _ , not to reconstruct exact spatial mass distribution , the more that many effects like the influence of global magnetic fields in a galaxy on the rotation of gas in the galaxy outskirts , are normally not taken into account , and these fields may be important for accounting for flat rotation curves @xcite .    global disk model is a natural approximation for flattened galaxies such as spiral ones , as far as it is not assumed in advance , sometimes by force , that each spiral galaxy must necessarily be immersed in an invisible halo of dark matter .
moreover , we are able to find such mass distributions which perfectly account for rotation curves , while in models assuming massive halo these rotation curves often can not be satisfactorily reconstructed , eg . in the case of ngc 7793 in @xcite , or
the dark halo found is too weak to account for stability criteria ( ngc 4736 in @xcite ) ( whereas , historically , the role of such haloes was not only to explain flat rotation , but also to stabilize disks ) .
although our model is very simple , similarly as other galaxy models customarily used by astronomers , it is better by that it does not assume constant mass - to - light ratios but it predicts the ratio as a resultant mass - to - light profile .
our results pertain only to the particular galaxies examined by us so far . from the fact that there is no need of non - baryonic dark matter in galaxies ngc 4736 or ngc 7793 , one can not conclude that one does not need it in other spiral galaxies . to claim otherwise ,
further studies are needed . nonetheless , we find interesting the existence of spiral galaxies without non - baryonic dark matter or galaxies with small abundance of such matter
.    @r@|@c@|@c@|@c name & ngc 7793 & ngc 1365 & ngc 6946 + incl .
angle [ deg ] & @xmath51 @xmath52 & @xmath53 @xmath54 & @xmath55 @xmath54 + distance [ mpc ] & @xmath56 @xmath52 & @xmath57 @xmath58 & @xmath59 @xmath60 + morphological type & sad @xmath52 & sbb @xmath54 & sabc @xmath54 + @xmath61 @xmath62}$ ] & @xmath63 @xmath64 & @xmath65 @xmath66 & @xmath67 @xmath68 + @xmath69 @xmath70}$ ] & @xmath71 @xmath52 & @xmath72 @xmath66 & @xmath73 @xmath60 + @xmath74 @xmath70}$ ] & - & - & @xmath75 @xmath76 +    [ cols= " > " , ]      +
we acknowledge the usage of the hyperleda database ( http://leda.univ-lyon1.fr )      athanas & sellwood ( 1987 ) iaus , 117 , 300a athanassoula , e. , bosma , a. , & papaioannou , s. ( 1987 ) a&a , 179 , 23 bahcall , j. n. , & casertano ( 1985 ) apj , 293 , l7 battaner , e. ; garrido , j. l. ; membrado , m. ; florido , e ( 1992 ) nature , 360 , 6405 , p. 652
bernabei r. , belli p. , montecchia f. at al .
, 2007 , nucl .
b , proc . suppl . , 166 , 87 binney , j. , & tremaine , s. ( 1987 ) , galactic dynamics ( princeton univ . press )
bosma , a. ( 1978 ) , ph.d .
thesis , univ .
groningen bosma , a. ( 1999 ) , 1999aspc .. 182 .. 339b boomsma r. , t. a. oosterloo , f. fraternali , j. m. van der hulst and r. sancisi ( 2008 ) a&a , 490 , 2 , 555 - 570 bratek  . , jaocha j. , kutschera m. , 2008 , mnras , 391 , 1373 buchhorn , m. ( 1992 ) , ph.d
. thesis , australian natl .
carignan c. , 1985 , apjs , 58 , 107 carignan c. , puche d. , 1990 , aj , 100 , 394 ciardullo , durrell , at al ( 2004 ) apj , 614 , 167c clowe d. , brada m. , gonzalez a.h . , markevitch m. , randall s.w . , jones c. , zaritsky d. , 2006 , apj , 648 , l109 de blok , mcgaugh , rubin ( 2001 ) aj , 122 , 2396d de blok , w.  j.  g. and walter , f. and brinks , e. and trachternach , c. and oh , s .- h . and kennicutt , r.  c. , 2008 , aj , 136 , 2648 - 2719 dicaire , i. ; carignan , c. ; amram , p. ; marcelin , m. ; hlavacek - larrondo , j. ; de denus - baillargeon , m .-
m . ; daigle , o. ; hernandez , o. ( 2008 ) aj , 135 , 6 , 2038 - 2047 .
evans ( 2001 ) 2001idm .. conf ... 85e fall , s. m. ; efstathiou , g. ( 1980 ) mnras.193 .. 189f hohl f. , ( 1976 ) aj ..... 81 ... 30h herrmann and ciardullo ( 2005 ) aipc , 804 , 341h hoekstra h. , van albada t.s . , sancisi r. , 2001 , mnras , 323 , 453 hyperleda database ( http://leda.univ-lyon1.fr ) jaocha j. , bratek  . , kutschera m. , kolonko m. , 2007 , acta phys . pol .
b , 38 , 3859 jaocha j. , bratek  . , kutschera m. , 2008 , apj , 679 , 373 jorsater s. , van moorsel g.a . , 1995 , aj , 110 , 2037 kalnajs , a. j. ( 1983 ) in iau symp . 100 , internal kinematics and dynamics of disk galaxies , ed .
e. athanassoula ( dordrecht : reidel ) , 87 kalnajs ( 1987 ) iaus , 117 , 289k kent , s. m. 1986 , aj , 91 , 1301 kent s.m . , 1987 , aj , 93 , 816 kutschera , m. ; jaocha , j. ( 2004 ) acta physica polonica b , 35 , 10 , p.2493 lake , g. , & feinswog , l. ( 1989 ) aj , 98 , 166 maller , ariyeh h. ; simard , luc ; guhathakurta , puragra ; hjorth , jens ; jaunsen , andreas o. ; flores , ricardo a. ; primack , joel r. ( 2000 ) apj ... 533 .. 194 m mulder , p. s. , and van driel , w. ( 1993 ) a&a , 272 , 63 munoz - tunon , c. , prieto , m. , beckman , j. , and cepa , j. , 1989 , ap&ss 156 , 3001 olling , r. p. ( 1996 ) , aj , 112 , 481 ostriker , peebles ( 1973 ) apj , 186 , 467o palunas & williams ( 2000 ) aj , 120 , 2884p rix , h .- w . , &
zaritsky d. ( 1995 ) apj , 447 , 82 rubin , v. c. , thonnard , n. , & ford , w. k. j. ( 1978 ) apj , 225 , l107 sackett , p. d. , rix , h .- w . ,
jarvis , b. j. & freeman , k. c. ( 1994 ) , apj , 436 , 629 sackett ( 1996 ) iaus , 173 , 165s sackett ( 1997 ) apj , 483 , 103s snchez - salcedo , f. j. ; reyes - ruiz , m. 2004 apj ... 607 .. 247s sanders r. h. ( 1996 ) apj , 473 , 117 - 129 sellwood ( 1980 ) a&a .. 89 .. 296s sellwood ( 1983 ) iaus .. 100 .. 197s sellwood , j. ( 1998 ) , in asp conf .
182 , galaxy dynamics , ed .
d. r. merritt , m. valluri , & j. a. sellwood ( san francisco : asp ) , 339 sellwood , kosowsky ( 2000 ) ? ?
arxiv : astro - ph/0009074v1 , gas & galaxy evolution asp conference series , vol . *
* volume * * , 2000 j. e. hibbard , m. p. rupen and j. h. van gorkom , eds .
sellwood , j. a. , & evans , n. w. ( 2001 ) , apj   data sofue y. , 1997 , pasj , 49 , 17 sug - whan k. , muk - suk c. , 1984 , jkas , 17 , 23 toomre , a. ( 1964 ) apj , 139 , 1217 t tully r.b . , verheijen m.a.w .
, pierce m.j . , huang jia - sheng , wainscoat r.j . , 1996 ,
aj , 112 , 2471 van albada , g. d. , & sancisi , r. ( 1986 ) phil .
london a , 320 , 447 verheijen m.a.w . , sancisi r. , 2001 , a&a , 370 , 765 williams m. j. , m. bureau & michele cappellari ( 2009 ) mnras.tmp.1591w , young j.s . , scoville n.z .
, 1991 , ara&a , 29 , 581 znmar snchez r. , sellwood j. a. , weiner b. j. , williams t. b ( 2008 ) apj , 674 , 2 , 797 - 813