Newtonian explanation of galaxy rotation curves based on distribution of baryonic matter Konstantin Pavlovich, Alex Pavlovich & Alan Sipols June 9, 2014
360 W Illinois St., 8F, Chicago IL 60654, USA email: [email protected]
ABSTRACT Circular velocities of stars and gas in galaxies generally do not decline in accordance with widely expected Keplerian fall off in velocities further from the galactic nucleus. Two main groups of theories were proposed to explain the supposed discrepancy—first, the most commonly accepted, is the suggestion of the existence of large non baryonic dark matter halo, and, second are theories advocating some modification to the law of gravity. So far however, there is no empirical evidence for either dark matter or modified gravity. Here we show that a broad range of galaxy rotation curves can be explained solely in accordance with Newton’s law of gravity by modeling the distribution of baryonic matter in a galaxy. We demonstrate that the expectation of Keplerian fall off is incorrect, and that a large number of likely galaxy mass distribution profiles should in fact produce flat or accelerating rotation curves similar to those observed in reality. We further support our theoretical findings with the model fit of 47 real galaxies’ rotation curves, representing a broad range of galactic types and sizes, and achieving correlation of expected and observed velocities of over 0.995 for all cases. Our results make theories of exotic dark matter or modified gravity unnecessary for the explanation of galaxy rotation curves.
INTRODUCTION Since Fritz Zwicky’s observations 1 of galaxy rotation curves in 1930s astronomers were puzzled by the fact that circular velocities of stars and gas further from the nucleus of the galaxy generally do not decline in accordance with widely expected Keplerian fall off, similar to the decline in the Solar system 2. Multiple studies 3,4,5,6,7 confirmed that galaxy rotation curves are mostly flat, with some galaxies showing modestly declining 8 and some accelerating circular velocities further away from a nucleus. There are two main groups of theories proposed to explain these phenomena. First, the most commonly accepted, is the suggestion of existence of large non baryonic dark matter galactic halo 9,10,11,12 . Second, there are theories advocating some modification to the law of gravity, with Modified Newtonian Dynamic 13,14,15 and Scalar– Tensor–Vector Gravity 16,17 , being top contenders. In spite of intensive search, all attempts to detect non baryonic dark matter failed. There is also no empirical support for modifications to the law of gravity.
Most of the research in the past approached the problem of galaxy rotation curves assuming spherically symmetric gravitational field due to a point source, using the gravity law in the form F = G M m / R2 inside the galaxy disk with the entire galactic mass as a parameter. However, gravity law in this canonical form is only appropriate for two point masses or spherical objects, and only outside of the sphere’s radius. We believe it is incorrect to apply it in this form to the entire galaxy mass for calculations within the galaxy disk, where gravitational forces have an integral effect from multiple directions. In our opinion, this oversimplification created expectation of Keplerian fall off, which, in turn, led to the rise of theories of non baryonic dark matter or modified gravity to explain the supposed discrepancy.
Kepler’s laws are applicable for a system with a massive star in the center orbited by a few planets with immaterial mass. In such system, one can ignore mutual gravitational pull of planets, although even in
Page 1 of 12 this simple case Newton’s adjustments to Kepler’s 3rd law are required for precise calculation. Once we increase the number of objects in a system to billions, especially if the aggregate mass of these objects exceeds the central nucleus mass, Kepler’s laws become inapplicable. Galaxy is this type of system—the mass of a central black hole is several orders of magnitude lower than the mass of a galaxy. Such system should be viewed as complex and uneven swarm of objects, which behavior can only be analyzed via integration of mutual gravitational forces between all those objects.
INTEGRATION OF GRAVITATIONAL FORCES IN A GALAXY
To calculate a cumulative gravitational force exerted upon any selected point in a galaxy, we have to integrate gravitational forces between selected point and all other points within a galaxy. In a two dimensional 1 model of a galaxy the gravitational attraction force between random point A on x axis with mass MA at a distance r from the center (Figure 1) and random point B(x,y) with MB is determined by the gravity law as:
2 FA = G MB MA / (AB) (1)
Figure 1 | Theoretical illustration of gravity forces in a galaxy
The vector of this force is at an angle α to the X axis (Figure 1). The horizontal component FAcos( α) generates a centripetal gravity force and vertical component FAsin( α) generates a tangential gravity force. In general, the tangential component is significantly weaker than the centripetal as the gravitational pull of the mass above and below the X axis largely cancel each other out. In the case of axisymmetric systems the tangential component F Asin( α) is completely compensated by an equal force FAsin( α) from a mirror point B’. Integrated over galactic disk area D with radius R, the centripetal and tangential components of the gravity force are calculated as:
M x, y cos α Centripetal force F = G M dx dy 2.1 AB
M x, y sin α Tangential force F = G M dx dy 2.2 AB 2 2 2 Where АВ = (r – xb) + y b ; Cos(α) = (r – xb)/(AB); Sin(α) = y b/(AB);
1 The two dimensional simplification of the model is warranted given that distances in the galactic plane are significantly greater than the distances perpendicular to the plane. Thicker profile in a bulge can be accounted for in a two dimensional model by increased mass density in the center.
Page 2 of 12
2 Given that centripetal force F Ac=МАVА /r we find orbital velocity V A as:
M x, y r – x V = G r dx dy 3 AB Except for the uniform mass distribution case, integral (3) is unsolvable with direct methods for even simple functions of mass distribution M=f(x,y). However, by breaking the galaxy disk into a sufficiently large number of small areas as a matrix, we can approximate this integral with a sum:
, – = 4 AB This approximation of integral (3) with a sum (4) is a crucial step. Some earlier researchers 18,19 acknowledged the integral (3), yet largely resorted to models with spherically symmetric gravitational field due to a point source, not reflecting galaxy disk structure and ignoring the impact of cumulative gravitational pull from low mass density galactic periphery. As we demonstrate below accounting for galaxy disk structure and peripheral mass outside of the measurement orbit makes a profound impact on galaxy rotation curves.
GALAXY MASS DENSITY MODEL
In order to solve equation (4) for a large number of galaxies of various types and sizes, we created a galaxy Mass Density Model (MDM). The model is based on N by N matrix in which each cell represents a region of space scaled automatically to selected galaxy size . Setting matrix resolution at 41 by 41 enables approximation of the integral (3) with error tolerance of no more than 0.35%, which for practical purposes is sufficient to model key features and behavior of galaxy rotation curves. In principle, however, it is easy to make this model as granular as desired. Each cell in the matrix is ascribed a nominal mass density number from zero to 100, where zero is vacuum and 100 the densest area (Figure 2). The mass of each cell is then calculated as a total galaxy mass multiplied by respective cell nominal density number and divided by the sum of nominal density numbers of all cells. Calculating gravity forces between selected cell on the axis of measurement and all other cells we solve equation (4). By repeating this procedure for other cells along the axis of measurement we determine rotation curve for selected mass distribution and axis.
Figure 2 | Galaxy mass density model setup
Page 3 of 12
For an axisymmetric mass density distribution the whole matrix can be setup with a one dimensional array from nucleus to the edge of the galaxy (Figure 2a). Importantly, axisymmetric model setup allows us to run MDM in both directions—from mass distribution to rotation curve, and from rotation curve to implied axisymmetric mass distribution. As we demonstrate below, such simple axisymmetric setup is largely sufficient to explain general behavior and shapes of galaxy rotation curves even for spiral and other non axisymmetric galaxies. This is due to the fact that in most cases centripetal forces are significantly stronger than tangential forces generated by non axisymmetric structures such as spiral arms. When necessary however, the two dimensional functionality of our matrix enables us to model any non axisymmetric distributions including all types of irregular features such as spirals and bars (Figure 2b) and account for tangential gravitational forces. It also allows us to utilize galaxy photographic data to infer approximate mass distribution and estimate circular velocities for various axes for galaxies with non axisymmetric mass distributions.
To translate galaxy photograph into an approximate mass density distribution matrix, it is possible to measure average relative luminance for each rectangular area of the source image corresponding to a target matrix cell, and correct resulting luminance matrix values for sky background light and the mass density to relative luminance ratio (MD/RL). The MD/RL adjustment function can be generalized from correlations between the distribution profiles of relative luminance and mass density, evaluated on a sample of (a) galaxy images with known view angle and (b) mass density distribution estimates generated by MDM from rotation velocity inputs. For conversion of edge on to face on 3D projection the matrix can be further transformed taking into account galaxy spatial orientation parameters. Lastly, automatic heuristic correction algorithm and/or precise manual override can be applied to those individual cells where value is likely distorted by objects present on the image but known or suspected to be foreign to the galaxy.
SHAPES OF GALAXY ROTATION CURVES IN RELATION TO MASS DISTRIBUTION
To explore theoretical behavior of galaxy rotation curves depending on a range of axisymmetric mass distributions we used a model of a hypothetical galaxy with a mass 10 11 MÍ and a radius 20 kpc. We setup several axisymmetric mass density distribution profiles as a function of a distance from the nucleus 1/(x+1) k, where x is the distance (1, 2, …, n kpc) and k is a real number (Figure 3b). As we can see (Figure 3a), the rotation curve flattens rapidly with even very low mass density at a galactic periphery. At k=1.5, which implies 1% periphery to center density ratio, the rotation curve is practically flat. For galaxies with relatively denser periphery at k<1.4 (implying edge to center density ratio above 1.4%) we should expect accelerating circular velocities further from the nucleus. Of course, galaxy mass distribution does not need to follow a specific function. Thus, for illustrative purposes we also constructed a hypothetical “Ideal Flat Curve” with the circular velocity set at a constant 178 km/sec at every point of the curve (Figure 3; red dotted line) and using MDM calculated implied mass distribution required to produce this rotation curve. This hypothetical Ideal Flat Curve implies even lower 0.2% periphery to center density ratio. It requires faster initial density drop from the nucleus, followed by slower density decline in the galactic disk. In general, Keplerian decline in velocities can only be achieved with almost instant density drop (Figure 3; green dashed line) which implies practically all galaxy mass in the center. Density distributions around the Ideal Flat Curve and up to k=1.5 produce approximately flat rotation curves. Density profiles above this area generate accelerating rotation curves.
Page 4 of 12
Figure 3 Expected rotation curves of a hypothetical galaxy depending on mass density profile
3a Galaxy Rotation Curves 3b Mass Density Profiles 250 100 k =1 k = 1,2 90 k = 1,4 k = 1,5 200 80
70 k = 1,6 k = 1,7
150 60 k = 1,8 k = 2
50 Hypothetical Keplerian 100 40 "Ideal Flat Curve"
Mass density(% of max) Mass 30 Circular velocity (km/sec) velocity Circular
50 20
10
0 0 0 2 4 6 8 10 12 14 16 18 20 0 2 4 6 8 10 12 14 16 18 20 Distance from galactic nucleus (kpc) Distance from galactic nucleus (kpc)
Conceptually, there are two reasons why very low peripheral density is sufficient to produce flat and even accelerating rotation curves. First, galactic periphery covers much larger area than the bulge and even at a low density cumulative mass of the disk become comparable with that of the bulge. Second, a low mass close to the point of measurement can exert more significant gravitational force locally than much heavier but distant galactic nucleus. The conclusion is inescapable—a system with mass distribution resembling that of any typical galaxy should in fact produce a rotation curve similar to those we observe in reality, without exotic dark matter, but strictly in accordance with known law of gravity.
MODEL FIT TO ROTATION CURVES OF REAL GALAXIES
To test our Mass Density Model on real galaxies we selected a sample of 47 galaxies covering twelve galaxy types (Sa, SAb, SBc, etc.), broad spectrum of shapes of rotation curves and velocities, and a wide range of galactic diameters from 2kpc to 80kpc. Using MDM, we calculated galactic mass and its distribution solving for the best fit to respective galaxy rotation curve and diameter. For 46 of 47 galaxies MDM automatically achieved extremely close fit with correlation over 0.995 between observed velocities and those generated by the model (Figure 4). Such high correlation means that not only MDM fit general shapes and velocities of rotation curves, but also replicated all nuanced features of rotation curves such as local maxima and minima in velocities. Consistent with our theoretical observations, these results were achieved with mass distribution profiles typically showing rapid drop in density from a nucleus (median density of 5.9% at a distance 15% of R max from the nucleus), followed by a more gradual density decline in a galactic disk reaching median mass density of 0.1% at the optical edge of a galaxy.
The only instance where MDM automatic fit initially failed to achieve correlation over 0.995 was the case of galaxy NGC 5194. Specifically, the tail end of the rotation curve did not fit (Figure 4, dotted line), suggesting an unlikely high mass density at a galactic periphery. Further investigation 20,21 revealed that NGC 5194 has a nearby satellite galaxy NGC 5195 stealing matter from NGC 5194 periphery (Figure 2c). Once the adjustment for the satellite galaxy mass and distance was made, the rotation curve matched perfectly with correlation over 0.9999 (Figure 4). We believe this case illustrates MDM predictive and explanatory power, as it forces us to consider close galactic neighborhood that impacts the behavior of a particular galaxy.
Page 5 of 12
Figure 4 | Model fit to rotation curves of real galaxies Name of each galaxy and estimated mass is shown on each chart. Blue circles represent observed rotation curve (data from http://www.ioa.s.u tokyo.ac.jp/~sofue/); Red line denotes MDM generated expected rotation curve; and Green dashed line shows implied axisymmetric mass distribution. Left Y axis shows speed in km/sec; Right Y axis is mass density in % of maximum; and X axis is distance from the nucleus in kpc. Includes 12 of total 47 galaxies; the remaining charts can be found in the extended figures.
400 100 400 100 400 100 350 350 350 80 80 80 300 300 300
250 60 250 60 250 60 200 200 200 40 40 40 150 NGC 2841 (Sb) 150 NGC 1097 (SBb) 150 NGC 2590 (SBc) 100 Est. Mass = 433 x 10 9 MÍ 100 Est. Mass = 384 x 10 9 MÍ 100 Est. Mass = 180 x 10 9 M 20 20 20 50 50 50 0 0 0 0 0 0 0.0 6.0 12.0 18.0 24.0 30.0 36.0 0.0 5.0 10.0 15.0 20.0 25.0 30.0 0.0 4.0 8.0 12.0 16.0
300 100 300 Dotted line is best fit 100 300 100 without adjustment 250 80 250 for the satellite 80 250 80 200 200 NGC 5195 200 60 60 60 150 M31 "Andromeda" (SAb) 150 150 Est. Mass = 280 x 10 9 M 40 NGC 5194 (Sc) 40 NGC 2708 (SBc) 40 100 100 9 100 Est. Mass = 97 x 10 M Est. Mass = 38 x 10 9 M 50 20 50 20 50 20
0 0 0 0 0 0 0.0 5.0 10.0 15.0 20.0 25.0 30.0 0.0 2.0 4.0 6.0 8.0 10.0 12.0 14.0 0.0 1.0 2.0 3.0 4.0 5.0
250 100 250 100 250 100
200 80 200 80 200 80
150 60 150 60 150 60 Milky Way (SBc) NGC 4527 (SBc) 100 Est. Mass = 149 x 10 9 M 40 100 NGC 1808 (SAB) 40 100 Est. Mass = 73 x 10 9 M 40 Est. Mass = 60 x 10 9 M 50 20 50 20 50 20
0 0 0 0 0 0 0.0 5.0 10.0 15.0 20.0 0.0 4.0 8.0 12.0 16.0 0.0 2.5 5.0 7.5 10.0 12.5
100 NGC 4448 (Sba) 300 100 160 NGC 3495 (Scd) 250 100 Est. Mass = 12 x 10 9 M 140 Est. Mass = 13 x 10 9 M 250 80 80 200 80 120 200 60 100 60 150 60 150 80 Observed circular velocity NGC 6951 (SBbc) 40 60 40 100 40 100 Est. Mass = 93 x 10 9 M Expected circular velocity 40 Mass density distribution 50 20 20 50 20 20 0 0 0 0 0 0 0.0 2.0 4.0 6.0 8.0 10.0 0.0 1.0 2.0 3.0 4.0 5.0 0.0 0.5 1.0 1.5 2.0
It is important to make a few observations about galactic masses. For the majority, even well researched galaxies, there is no consensus on galaxy mass between various sources. For example, the range of estimates for our close neighbor Andromeda galaxy M31 includes anything from 20*10 10 MÍ(16) to 1,230*10 10 MÍ(22). Given that galaxy mass is not directly measured, but estimated from observed velocities on the basis of various galactic models, it is not surprising to see such broad range of estimates. As we explained above, our mass estimates and their respective distributions are computed automatically from our MDM model solving for the best fit to rotation curve. In the vast majority of cases, our mass estimates are below those commonly accepted. This leads to a conclusion that observed rotation curves can be readily explained without extra matter (dark or not) traditionally ascribed to galaxies.
Page 6 of 12
TANGENTIAL GRAVITY FORCES IN SPIRAL GALAXIES The additional advantage of MDM over earlier research methods is that MDM enables the analysis of tangential gravity forces and their impact on galaxy dynamics. In most cases tangential gravity force will be significantly weaker than centripetal force (Figure 5). This is due to the fact that sides of the galactic disk pull in opposite directions and largely cancel itself even with non symmetric mass distribution. However, there are instances, when significant mass such as a stellar cluster is located close to the point of measurement, as is the case of galaxy M101 at 160 degrees axis (Figure 5b). In such situations, particularly on galactic periphery where centripetal force is weaker, the tangential force generated by nearby mass may exceed the centripetal force, significantly impacting the vector of the gravity force (Figure 5c) and consequently the shape of the rotation curve. This modeling of gravity forces at various axes allows us to make a testable prediction that rotation curves should differ in a spiral galaxy depending on the axis of measurement.
Figure 5 | Tangential gravity forces in a spiral galaxy. 5a Model of Pinwheel galaxy M101: Axis at 0 degrees 5b Model of Pinwheel galaxy M101: Axis at 160 degrees
00000000000000012223520000000000000000000 0000000001 7 7 8 6 4 410161824151210 9 9 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
00000000000000666865422210000000000000000 0000000002 31423 8 913 91025181515101414 8 8 7 0 0 0 0 0 0 0 0 0 0 0 0 0
000000000000272330131175432210000000000000000 000000 0 0 2 3 3 821 81013 91328281713131412111110 4 1 0 0 0 0 0 0 0 0 0 0 0
000000000001311212013965422210000000000000000 000000 0 0 4 6 51524151313142033351912111410 9 6 4 7 5 4 0 0 0 0 0 0 0 0 0 0
00000000005582326186874221200000000000000000 0000000 0 3 6 7 81815252726302116151511 9 2 7 4 5 2 3 4 2 1 1 0 0 0 0 0 0 0
00000000011110142921167983220001000000000000000 0000 0 0 0 0 6 8 710231821283439161515 8 4 2 3 71212 510 4 4 3 1 2 0 0 0 0 0 0
00000000021218323417116672211122100000000000000 0000 0 0 0 1 7 9 9 26 20 16 24 45 33 30 14 11 11 6 5 6 3 10 19 10 24 18 11 10 10 5 2 0 0 0 0 0 0
000000000514273835171010775555654333200000000000 000 0 0 0 0 2 7 13 11 19 25 19 21 45 27 14 13 12 11 6 8 11 12 21 34 16 20 24 20 17 13 5 3 0 0 0 0 0 0
0 0 0 0 0 0 0 0 3 7153131281413141410 910121917151513 8 6 4 3 2 1 2 0 0 0 0 0 0 0 0 0 0 0 0 0 4 81314142125302331211516131211162033442726191919212111 9 4 0 0 0 0 0
0 0 0 0 0 1 0 1 5 9 14 19 26 21 19 20 23 17 17 16 17 20 27 24 26 34 29 14 9 7 7 5 3 0 1 0 0 0 0 0 0 0 0 0 0 0 0 4 9141526312617202323172116171420223841 24 21 14 13 14 14 24 31 25 21 2 2 0 0 0
0 0 0 0 0 0 0 2 18 20 17 15 19 17 22 22 23 28 36 23 25 33 39 40 34 32 30 27 19 1411 8 6 3 0 0 0 0 0 0 0 0 0 0 0 1 1 511141618272419202121252431372424 33 33 25 21 17 14 12 12 7 24 23 23 19 8 2 0 0 0
0 0 0 0 0 0 0 2 9 17 18 13 14 17 21 27 28 29 31 28 26 31 37 39 39 31 31 27 25 19 191511 5 2 0 0 0 0 0 0 0 0 0 0 0 312142120232325241823283833474439 47 37 33 41 18 18 16 12 14 12 11 17 19 23 31 13 0 0 0
0 0 0 0 0 0 0 4 10 26 22 12 14 22 34 37 33 33 27 28 30 29 33 33 31 32 40 34 27 21171512 9 5 2 2 0 0 0 0 0 0 0 1 2 3 9152025282623172634575750524446 41 37 49 51 31 26 22 13 7 5 11 11 8 24 27 19 5 1 0
0 0 0 0 0 0 2 4 14 28 24 19 22 35 36 39 37 39 35 36 35 29 28 27 24 27 32 33 37 3322191417 8 6 0 0 0 0 0 0 0 0 3 3 8101419272825241726416554524546 43 43 34 39 45 45 39 30 19 20 16 5 11 11 22 22 19 4 4 1
0 0 0 0 0 0 2 9 16 40 27 23 30 36 34 38 35 39 41 46 46 47 40 32 29 23 23 24 25 30252218 161013 3 0 0 0 0 0 0 0 5 5 813203030 22252517214275 585744 39 39 35 33 43 35 41 40 40 23 21 20 12 8 10 14 15 21 9 7 1
0 0 0 0 0 0 3 7 17 42 38 35 39 46 38 37 39 50 47 61 59 66 52 43 37 28 26 19 20 25303424 191412 5 0 1 0 2 0 2 0 2 8 818172323 20212122346554 565246 42 39 40 29 29 24 34 31 40 33 27 35 8 11 11 9 7 21 10 8 1
0 0 0 0 0 0 4 6 18 31 36 37 50 47 36 38 45 54 52 59 63 56 50 53 45 31 24 18 16 19 21 26 26 31 17 16 7 2 0 0 3 0 0 0 4 8 14 14 17 16 20 18 22 27 22 35 69 62 57 58 65 55 60 61 40 42 39 20 21 26 33 22 35 29 15 11 8 7 8 6 6 4
0 0 0 0 0 0 2 7 14 30 40 48 56 51 39 43 49 59 66 71 72 65 62 63 57 39 25 21 13 14 16 21 25 22 19 15 8 2 0 0 0 3 0 0 2 8 15 18 21 17 16 14 23 23 30 29 59 64 61 69 75 69 69 61 60 51 30 28 26 26 33 29 24 23 13 10 8 7 7 11 8 0
0 0 0 0 0 0 1 7 13 33 39 50 53 49 51 51 53 63 73 72 71 69 64 56 55 47 36 26 16 13 14 16 19 21 18 15 8 4 0 0 2 0 0 0 2 7 15 19 22 19 16 14 18 25 36 53 76 63 63 80 75 78 79 65 64 63 42 35 36 33 42 27 17 10 8 14 8 7 16 18 15 1
0 0 0 0 0 0 2 12 13 28 38 42 52 59 51 56 57 75 80 84 87 82 73 67 66 63 44 36 25 18 14 16 17 19 17 14 10 8 2 0 0 0 0 0 5 12 17 31 31 25 21 14 13 16 36 44 63 63 75 80 91 95 86 69 69 66 63 48 36 53 60 36 15 10 8 8 4 6 37 54 20 9
0 0 0 0 0 0 2 4 11 27 42 36 41 50 53 59 58 76 83 90100 91 80 75 63 76 53 39 30 23 21 18 20 16 17 18 8 5 5 0 0 0 0 0 3 12 14 19 26 26 21 14 13 13 26 36 47 66 67 73 91100 90 78 69 65 58 44 36 42 42 41 20 12 7 7 4 2 10 31 21 8
0 0 0 0 0 0 1 9 10 22 34 49 39 45 55 56 65 69 78 86 95 91 84 80 63 64 59 29 22 27 22 20 23 23 20 13 8 3 3 1 0 0 0 0 0 13 10 16 24 34 30 19 16 18 25 39 57 56 64 64 82 87 84 83 76 58 56 55 45 37 38 27 24 18 5 5 3 2 2 4 21 12
0 0 0 0 0 0 0 3 14 19 27 38 37 36 44 58 65 66 69 73 79 78 75 69 61 62 69 35 22 21 21 22 30 30 19 14 10 3 2 1 0 0 0 2 2 8 17 14 18 22 25 25 20 19 24 31 45 63 62 65 71 72 73 75 57 57 59 53 45 39 49 27 19 14 4 2 0 0 0 0 9 9
5 0 0 0 0 0 2 4 18 24 28 41 42 39 36 48 63 63 64 65 69 69 65 58 57 54 65 34 19 17252528 272015 9 3 0 0 0 0 0 0 0 5 912192233 25242326283753 505672 71 73 63 53 56 51 59 41 36 36 34 22 10 3 0 0 0 0 0 6 4
4 6 0 0 0 1 3 5 7 12 20 36 60 53 36 35 42 51 60 61 60 55 53 46 52 56 75 42 21 17 24 26 28 25 21 14 12 4 1 0 0 0 0 0 0 2 5 11 15 17 21 37 33 23 23 29 43 52 66 56 63 59 66 59 49 51 51 49 52 42 38 27 11 4 1 0 0 0 0 0 0 5
6 9 12 4 2 2 4 7 8 10 15 27 42 33 26 28 30 39 42 40 40 39 42 44 57 58 65 41 26 17 23 24 23 20 16 14 11 5 0 0 0 0 0 0 1 3 6 15 15 19 21 27 34 32 27 24 32 40 47 59 61 47 54 45 38 43 39 51 53 50 39 28 13 12 2 0 0 0 0 0 0 0
2 12 21 31 10 6 4 8 8 10 17 25 29 33 26 21 20 24 29 29 35 39 39 45 52 54 57 34 3018242527 181614 9 4 0 0 0 0 0 0 1 0 3 8111425 27313232312728 2946 46 41 50 39 37 36 47 50 48 40 30 33 13 7 2 0 0 0 0 0 0 0
0 8 20 54 37 12 7 8 14 13 23 24 22 33 40 31 34 35 43 33 41 43 46 54 52 50 38 282320192431261513 8 3 0 0 0 0 0 0 0 2 1 5 7 719273031393933293036 3635393538 38463937363014 7 1 0 0 0 0 0 0 0 0
0 1 9 15 18 16 7 8 10 15 29 35 27 33 38 40 41 45 45 39 34 41 46 44 47 33 25 21 2120172628211413 7 2 0 0 0 0 0 0 0 0 0 3 4 4 914293234403731262827 33373934 34363935383118 6 2 0 0 0 0 0 0 0 0
0 0 1 8 11 7 7 8 11 11 8 18 20 21 23 30 39 45 51 49 37 47 39 37 31 24 21 17 23 282330251911 9 7 1 0 0 0 0 0 0 0 0 0 0 2 3 8131526242439332628313333 373635222327401717 7 4 0 0 0 0 0 0 0 0
0 0 0 3 6 8 8 7 9 10 8 12 16 20 19 22 26 31 41 33 28 33 24 24 17 16 21 15 21 31 21192026 9 8 6 0 0 0 0 0 0 0 0 0 0 0 0 0 3 315151727201723363629282734 221419242816 9 2 0 0 0 0 0 0 0 0 0
0 0 0 4 8 10 21 15 14 11 11 5 5 7 13 16 16 18 21 25 33 22 20 14 12 13 16 14 27 452416212310 7 6 3 0 0 0 0 0 0 0 0 0 0 0 0 0 1 4 5 619121016171723222217 141412222614 4 2 0 0 0 0 0 0 0 0 0
0 0 0 1 7 9 19 22 22 18 8 11 7 12 14 12 14 17 24 41 38 20 16 11 11 12 13 14 30 3328211818 8 7 6 4 0 0 0 0 0 0 0 0 0 0 0 0 0 1 2 1 1 5 5 5 9101423201917191518179400000000000
0 0 0 0 1 4 9 19 27 24 19 17 11 7 12 13 14 21 27 44 33 18 11 8 6 11 11 16 39 34 2725152415 5 3 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 2 2 7 71014131421261917209200000000000
0 0 0 0 0 1 5 13 31 23 23 23 24 14 14 16 19 19 26 34 21 12 6 5 8 15 15 21 30 26 1813151721 8 3 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2 1 2 2 3 7 61010172831311914 9 18200000000000
0 0 0 0 0 0 0 2 8 19 21 25 31 21 21 19 24 20 16 19 10 3 2 4 8 15 16 27 33 20 14 13 108231471000 000000000012222238961117353827157310100 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 2 2 4 91113172018241012 7 3 2 9111219352813 913 9 6 8 7 41000 0000000000122245786161621343218145000000000000 0 0
0 0 0 0 0 0 0 0 0 0 0 3 5 810111010 5 5 4 71014111317282510 910 4 211 3 0 0000 00000000001223456913182623141011220000000000000 0
00000000000025344324691214131215182418168202200000 0 000000000002455711132021118551000000000000000
00000000000021124577101191414101215141422200000000 0 000000000000532866232373150000000000000000
000000000000000000147889988118000000000000 000000000 00000002266020000000000000000000 Legend: Axis of measument Local mass generating strong tangential gravity force
5c Centripetal and tangential gravity forces in M101 at 0 and 160 degrees axes
3.0E-10 Tangential gravity force Centripetal gravity 2.5E-10 exceeds centripetal gravity force; 0 degrees axis force at the 160 degrees axis 2.0E-10 Tangential gravity force; 0 degrees axis 1.5E-10 Centripetal gravity force; 160 degrees axis 1.0E-10 Tangential gravity 5.0E-11 force; 160 degrees axis
0.0E+00 0 2 4 6 8 10 12 14 16 -5.0E-11 Gravity force (in N) exerted on force(inon mass exertedkgN) a 1 Gravity Distance from galactic nucleus (kpc) -1.0E-10
Page 7 of 12
CONCLUSION In this article we demonstrated that a broad spectrum of galaxy rotation curves can be explained without either dark matter or modified gravity, but by accounting for the distribution of baryonic matter and modeling its behavior strictly in accordance with the known law of gravity. Our Mass Density Model provides a theoretical framework for understanding of galaxy rotation and practical tool for modeling galaxy dynamics in relationship to its mass distribution. We also pointed two related mistakes of earlier research which led to the “invention” of dark matter and modified gravity theories—first, the expectation of Keplerian fall off in circular velocities in galaxies similar to that in the Solar system, and, second, the assumption of the entire galaxy mass as a spherically symmetric gravitational field due to a point source for calculations inside the galaxy disk. Our results show that both are fundamentally incorrect.
ACKNOWLEGEMENTS We want to sincerely thank Professor Yoshiaki Sofue of the University of Tokyo for making detailed rotation curve data on many galaxies publicly available.
REFERENCES 1. Zwicky, F. On the Masses of Nebulae and of Clusters of Nebulae. Astrophysical Journal, 86, 217 (1937). 2. Persic, M., Salucci, P., Stel, F. The universal rotation curve of spiral galaxies —I. The dark matter connection. Monthly Notices of the Royal Astronomical Society, 281, 27–47 (1995). 3. Rubin V.C., Burstein D., Ford W.K., Thonnard N. Rotation velocities of 16 Sa Galaxies and a comparison of Sa, Sb and Sc rotation properties. The Astrophysical Journal , 289, 81–104 (1984). 4. Rubin, V.C., Ford, W.K., Strom, K. M., Strom, S. E., Romanishin, W. Extended rotation curves of high luminosity spiral galaxies. Astrophysical Journal, 224, 782–795 (1978). 5. Sofue, Y., et al. Central Rotation Curves of Spiral Galaxies. University of Tokyo Astrophysics, 99, 03 (1999). 6. Rubin, V.C., Ford, W.K., Thonnard, N. Rotational properties of 21 SC galaxies with a large range of luminosities and radii, from NGC 4605 /R = 4kpc/ to UGC 2885 /R = 122 kpc/. Astrophysical Journal, 238, 471–487 (1980). 7. Binney, J. & Tremaine, S. Galactic Dynamics, 2nd Edition, Princeton University Press (2008). 8. Casertano, S. & van Gorkom, J. H. Declining rotation curves —The end of a conspiracy? Astronomical Journal, 101, 1231–1241 (1991). 9. Freeman, K.C. On the Disks of Spiral and S0 Galaxies. Astrophysical Journal, 160, 811 (1970). 10. Governato, F. et al. Bulgeless dwarf galaxies and dark matter cores from supernova driven outflows. Nature, 463, 203–206 (2010) 11. Van Albada T.S., Bahcall, J.H., Begeman, K. & Sanscisi R. Distribution of dark matter in the spiral galaxy NGC 3198. The Astrophysical Journal, 295, 305–313 (1985). 12. Memola, E., Salucci, P. & Babic, A. Dark Matter halos around isolated ellipticals. Astronomy & Astrophysics, 534, A50 (2011). 13. Milgrom M. A modification of the Newtonian dynamics as a possible alternative to the hidden mass hypothesis. Astrophysical Journal, 270, 365–370 ( 1983). 14. Begeman K.G., Broeils A.H. & Sanders R.H. Extended rotation curves of spiral galaxies: dark haloes and modified dynamics. Monthly Notices of the Royal Astronomical Society, 249, 523–537 (1991). 15. Milgrom, M. & Sanders R.H. Modified Newtonian Dynamics Rotation Curves of Very Low Mass Spiral Galaxies. The Astrophysical Journal, 658, L17–L20 (2007).
Page 8 of 12
16. Brownstein J.R. & Moffat, J.W. Galaxy rotation curves without nonbaryonic dark matter. The Astrophysical Journal, 636, 721–741, (2006). 17. Sanders R.H. A tensor vector scalar framework for modified dynamics and cosmic dark matter. Monthly Notices of the Royal Astronomical Society, 363, 459 (2005). 18. Sofue Y. The Mass Distribution and Rotation Curve in the Galaxy. Planets, Stars and Stellar Systems, Springer, Berlin, Vol. 5, ed. G. Gilmore., Chap 19 (2013). 19. Sofue Y. Rotation Curve and Mass Distribution in the Galactic Center—From Black Hole to Entire Galaxy. Astronomical Society of Japan, 65, 118 (2013). 20. Durrell P.R., Mihos C.J., Feldmeier J.J. & Ciardullo R. Kinematics of planetary nebulae in M51’s tidal debris. Astrophysical Journal, 582, 170 183, (2003). 21. Joon Hyeop Lee, et al. Hubble Space Telescope Pixel Analysis of the Interacting S0 Galaxy NGC 5195 (M51B), Astrophysical Journal, 754, 80, (2011). 22. Evans N.W., Wilkinson, M.I. The Mass of the Andromeda Galaxy. Monthly Notices of the Royal Astronomical Society, 316, 929 (2000).
EXTENDED FIGURES See pages 10–12 below for extended figures.
Page 9 of 12
Extended Figure 1 | Model fit to rotation curves of real galaxies Name of each galaxy and estimated mass is shown on each chart. Blue circles represent observed rotation curve (data from http://www.ioa.s.u tokyo.ac.jp/~sofue/); Red line denotes MDM generated expected rotation curve; and Green dashed line shows implied axisymmetric mass distribution. Left Y axis shows speed in km/sec; Right Y axis is mass density in % of maximum; and X axis is distance from the nucleus in kpc.
160 100 300 100 250 100 140 80 250 80 200 80 120 200 100 60 60 150 60 80 NGC 2403 (Sc) 150 NGC 253 (Sc) 9 9 60 Est. Mass = 59 x 10 MÍ 40 NGC 3521 (SBbc) 40 100 Est. Mass = 69 x 10 MÍ 40 100 9 40 Est. Mass = 129 x 10 MÍ 20 50 20 50 20 20 0 0 0 0 0 0 0.0 5.0 10.0 15.0 20.0 0.0 4.0 8.0 12.0 16.0 20.0 24.0 0.0 2.0 4.0 6.0 8.0 10.0
120 100 200 100 250 100
100 80 80 200 80 150 80 60 60 150 60 60 NGC 598 (Sc) 100 9 NGC 891 (Sb) Est. Mass = 12 x 10 MÍ 40 NGC 660 (SBa) 40 100 40 9 40 Est. Mass = 81 x 10 9 MÍ Est. Mass = 126 x 10 MÍ 50 20 20 20 50 20
0 0 0 0 0 0 0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 0.0 5.0 10.0 15.0 20.0 0.0 5.0 10.0 15.0 20.0 25.0
300 100 160 100 300 100 UGC 03691 (Sc) 140 9 250 80 Est. Mass = 17 x 10 MÍ 80 250 80 120 200 200 60 100 60 60 150 NGC 1068 (Sb) 80 150 NGC 1365 (SBb) 9 9 Est. Mass = 72 x 10 MÍ 40 60 40 Est. Mass = 250 x 10 MÍ 40 100 100 40 20 50 20 50 20 20 0 0 0 0 0 0 0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 0.0 2.0 4.0 6.0 0.0 10.0 20.0 30.0
300 100 250 100 250 100
250 80 200 80 200 80 200 60 150 60 150 60 150 NGC 1417 (SBb) IC 342 (SAB) Est. Mass = 100 x 10 9 MÍ 40 100 Est. Mass = 122 x 10 9 MÍ 40 100 UGC 2855 (SBc) 40 100 Est. Mass = 27 x 10 9 MÍ 50 20 50 20 50 20
0 0 0 0 0 0 0.0 2.0 4.0 6.0 8.0 10.0 0.0 5.0 10.0 15.0 20.0 0.0 1.0 2.0 3.0 4.0 5.0
Page 10 of 12
Extended Figure 2 | Model fit to rotation curves of real galaxies Name of each galaxy and estimated mass is shown on each chart. Blue circles represent observed rotation curve (data from http://www.ioa.s.u tokyo.ac.jp/~sofue/); Red line denotes MDM generated expected rotation curve; and Green dashed line shows implied axisymmetric mass distribution. Left Y axis shows speed in km/sec; Right Y axis is mass density in % of maximum; and X axis is distance from the nucleus in kpc.
250 100 350 100 300 100 NGC 3031 (Sb) 300 Est. Mass = 108 x 10 9 MÍ 200 80 80 250 80 250 200 150 60 200 60 60 150 NGC 2903 (Sc) 150 100 Est. Mass = 130 x 10 9 MÍ 40 40 NGC 3079 (SBc) 40 100 9 100 Est. Mass = 132 x 10 MÍ 50 20 20 20 50 50
0 0 0 0 0 0 0.0 4.0 8.0 12.0 16.0 20.0 24.0 0.0 3.0 6.0 9.0 12.0 15.0 18.0 21.0 0.0 3.0 6.0 9.0 12.0 15.0 18.0 21.0
200 100 300 100 300 100
80 250 80 250 80 150 200 200 60 60 60 100 NGC 3198 (SBc) 150 NGC 4565 (Sb) 150 NGC 7331 (Sa) Est. Mass = 120 x 10 9 MÍ 40 Est. Mass = 290 x 10 9 MÍ 40 Est. Mass = 270 x 10 9 MÍ 40 100 100 50 20 20 50 50 20
0 0 0 0 0 0 0.0 5.0 10.0 15.0 20.0 25.0 30.0 0.0 5.0 10.0 15.0 20.0 25.0 30.0 0.0 5.0 10.0 15.0 20.0 25.0 30.0
250 100 250 100 180 100 160 200 80 200 80 140 80 120 150 60 150 60 60 NGC 3672 (Sc) 100 NGC 3628 (Sb) Est. Mass = 85 x 10 9 MÍ 80 NGC 4062 (SBc) 100 Est. Mass = 95 x 10 9 MÍ 40 100 40 Est. Mass = 14 x 10 9 MÍ 40 60 50 20 50 20 40 20 20 0 0 0 0 0 0 0.0 2.0 4.0 6.0 8.0 10.0 12.0 14.0 0.0 3.0 6.0 9.0 12.0 0.0 1.0 2.0 3.0 4.0
300 100 180 100 300 100 160 250 250 80 140 80 80 200 120 200 60 60 60 100 150 NGC 4303 (SBbc) 150 NGC 4321 (SABbc) 80 40 Est. Mass = 55 x 10 9 MÍ 40 Est. Mass = 252 x 10 9 MÍ 40 100 NGC 4258 (SBbc) 60 100 Est. Mass = 177 x 10 9 MÍ 40 50 20 20 50 20 20 0 0 0 0 0 0 0.0 10.0 20.0 30.0 0.0 3.0 6.0 9.0 12.0 0.0 5.0 10.0 15.0 20.0 25.0
Page 11 of 12
Extended Figure 3 | Model fit to rotation curves of real galaxies Name of each galaxy and estimated mass is shown on each chart. Blue circles represent observed rotation curve (data from http://www.ioa.s.u tokyo.ac.jp/~sofue/); Red line denotes MDM generated expected rotation curve; and Green dashed line shows implied axisymmetric mass distribution. Left Y axis shows speed in km/sec; Right Y axis is mass density in % of maximum; and X axis is distance from the nucleus in kpc.
300 100 200 100 250 100
250 80 80 200 80 150 200 60 60 150 60 150 100 NGC 4631 (SBcd) 40 Est. Mass = 62 x 10 9 MÍ 40 100 40 100 NGC 4569 (SABab) NGC 4736 (Sab) 9 9 Est. Mass = 78 x 10 MÍ 50 Est. Mass = 35 x 10 MÍ 50 20 20 50 20
0 0 0 0 0 0 0.0 2.0 4.0 6.0 8.0 0.0 5.0 10.0 15.0 0.0 2.5 5.0 7.5 10.0
250 100 300 100 250 100
200 80 250 80 200 80 200 150 60 60 150 60 150 NGC 5033 (Sc) NGC 5055 (SBc) 100 NGC 4945 (SBcd) 40 Est. Mass = 357 x 10 9 MÍ 40 100 40 100 Est. Mass = 210 x 10 9 MÍ Est. Mass = 94 x 10 9 MÍ 50 20 50 20 50 20
0 0 0 0 0 0 0.0 5.0 10.0 15.0 20.0 0.0 10.0 20.0 30.0 0.0 10.0 20.0 30.0 40.0
250 100 250 100 250 100
200 80 200 80 200 80
150 60 150 60 150 60 NGC 6946 (Sc) NGC 5457 (Sc) 100 40 100 9 40 100 Est. Mass = 92 x 10 9 M 40 NGC 5236 (SBc) Est. Mass = 116 x 10 MÍ Í Est. Mass = 160 x 10 9 MÍ 50 20 50 20 50 20
0 0 0 0 0 0 0.0 10.0 20.0 30.0 40.0 0.0 5.0 10.0 15.0 0.0 4.0 8.0 12.0
300 100 300 100
250 80 250 80 200 200 60 60 NGC 6674 (SBb) 150 NGC 5907 (Sc) 150 Est. Mass = 11 x 10 9 MÍ Est. Mass = 260 x 10 9 MÍ 40 40 100 100
50 20 50 20
0 0 0 0 0.0 10.0 20.0 30.0 0.0 0.2 0.4 0.6 0.8 1.0
Page 12 of 12