<p>The rough answer is: just like the sun makes it hard or impossible to see planets and stars during the day, it dominates the gravity sky. </p>

<p><a href="https://i.sstatic.net/p4C99.png" rel="noreferrer"><img src="https://i.sstatic.net/p4C99.png" alt="Linear plot of the gravity sky"></a></p>

<p>But there are interesting patterns there if we view the sky using a logarithmic scale.</p>

<p><a href="https://i.sstatic.net/UkCuN.png" rel="noreferrer"><img src="https://i.sstatic.net/UkCuN.png" alt="Log plot of the gravity sky"></a></p>

<p>The above plots are a combination of the gravitational force of the sun, moon and planets, the stars in the Hipparchos catalog, the gravitational field of the Milky Way, and the nearby galaxies. Each component have a somewhat different effect.</p>

<p><strong>The solar system</strong></p>

<p>The force from the sun is shockingly big: 0.0059 N (per kilogram - I will omit that henceforth). At first this might appear impossible since we never notice the force, but remember that the setup for the question is a fixed point in space rather than the free-falling orbit of Earth. Because of the orbit this force exactly cancels.</p>

<p>Note that for spherical objects we only need to calculate the force from the direction of their centre of mass, not for their full extent. Hence regardless of what resolution I use for my map planets and stars will only change one box. </p>

<p>The force from the moon is <span class="math-container">$3.3198\cdot 10^{-5}$</span> N, a factor of 179 less. (Why do we experience tides from the moon more than from the sun? Because the gradient for the gravitational field is steeper for the nearby moon than for the remote sun). The other planets have much less effect, from <span class="math-container">$10^{-6.8}$</span> N for Jupiter down to <span class="math-container">$10^{-9.5}$</span> N for Neptune (and Pluto at <span class="math-container">$10^{-13.4}$</span>).</p>

<p>For positions I used the Matlab Aerospace toolbox planetEphemeris function. Incidentally, I might well have messed up some of the many different coordinate frames here, so don't use these maps to navigate. The important thing for this question is the magnitudes of the forces anyway.</p>

<p><strong>Stars</strong></p>

<p>I added the field of the stars in the Hipparchos catalog, but they do not show up even on my logarithmic map because they are overpowered by the galaxy at large. If we just plot their contributions to different boxes we get the following map:</p>

<p><a href="https://i.sstatic.net/KLlQy.png" rel="noreferrer"><img src="https://i.sstatic.net/KLlQy.png" alt="Log plot of the Hipparchos map"></a></p>

<p>Note that the peak force in each box is on the order of <span class="math-container">$10^{-14}$</span> N: even Pluto matters more than nearby stars. This map also shows an annoying problem: if we just sum the contributions in each 3.6 times 3.6 degree rectangle in the coordinate system the polar rectangles will be narrow and have smaller area; hence there are fewer stars included and less force. Dividing by <span class="math-container">$\cos(\delta)$</span> removes this bias and shows that the stars indeed pull on us roughly equally in all directions. </p>

<p>(I used a rough estimate of their masses: assuming them all to be on the main sequence I calculated their luminosity from their absolute magnitude, and then estimated the mass as <span class="math-container">$M\approx M_\odot (L/L_\odot)^{1/3.5}$</span>.)</p>

<p><strong>The Milky Way</strong></p>

<p><a href="https://i.sstatic.net/JOJXU.png" rel="noreferrer"><img src="https://i.sstatic.net/JOJXU.png" alt="Log picture of the galaxy force"></a></p>

<p>I modelled the Milky Way using the approximate values and formulas in Carroll &amp; Ostlie (p. 884), essentially a neutral gas disk, a thin disk, a thick disk, a central bulge, a stellar halo plus a dark matter halo. The disks for example have density scaling as <span class="math-container">$$\rho(r,z)=\rho_0 e^{-r/h_r}e^{-|z|/h_z}$$</span> while the dark matter halo has the <a href="https://en.wikipedia.org/wiki/Navarro%E2%80%93Frenk%E2%80%93White_profile" rel="noreferrer">Navarro-Frenck-White profile</a> <span class="math-container">$$\rho(r)=\frac{\rho_0}{\frac{r}{R_s}\left(1+\frac{r}{R_s}\right)^2}.$$</span> I integrate these functions along a ray extending from the part of the sky I am interested in. Note that the integral of the force will be <span class="math-container">$$F=\int_0^\infty \frac{G\rho(\mathbf{x}_0 + \mathbf{d} r)}{r^2} r^2 dr = G \int_0^\infty \rho(\mathbf{x}_0 + \mathbf{d} r) dr$$</span> since it is really an integral of an ever widening pyramid of area <span class="math-container">$r^2$</span>: the total force is only dependent on the total amount of mass there. Since these density fields change slowly I ignore lateral variations. Again there should be a correction for declination, but it does not make a great difference in the picture.</p>

<p>The end result is a fair amount of force, especially along the galactic plane. This causes "galactic tides" in stellar orbits and can sometimes rival the force from remote planets in the solar system. Again, we do not notice most of the gravitational field because the sun is free-falling in the galaxy.</p>

<p>Sag A* does not show up in this image because it is not very massive compared to the full galaxy. The force is just <span class="math-container">$9.3678\cdot 10^{-15}$</span> N.</p>

<p><strong>Nearby galaxies</strong></p>

<p><a href="https://i.sstatic.net/wsFKN.png" rel="noreferrer"><img src="https://i.sstatic.net/wsFKN.png" alt="Log near galaxy force"></a></p>

<p>I used the <a href="https://heasarc.gsfc.nasa.gov/W3Browse/galaxy-catalog/neargalcat.html" rel="noreferrer">NEARGALCAT catalogue</a> to estimate the effect of nearby galaxies. This is in many ways the most dodgy estimate of all. While we know positions well we have weak data on actual masses. I used the estimated masses whenever they existed, and a linear fit between angular extent and mass derived from those to impute masses for galaxies lacking values. I am a bit uncertain if the mass estimates are just photometric or include dark matter; if not including DM, they should be an order of magnitude higher.</p>

<p>Typical forces are on the order of <span class="math-container">$10^{-19}$</span> N, although heavier galaxies might pull us with <span class="math-container">$10^{-15}$</span> N. They are hence on par with the nearby stars: totally overpowered by the galactic gravitational field and the solar system. </p>

<p>Technically, since galaxies are extended non-spherical objects, I should integrate over areas here. But there are limits to my energy. </p>

<p>In short, the gravity sky has nontrivial patterns if one has gravity-eyes with a logarithmic response.</p>

CC BY-SA 4.0
2019-12-05T22:32:14.053Z
+1. I've often daydreamed about what an 'incident graviton flux' map would look like. This answers that brilliantly.
Ingolifs
CC BY-SA 4.0
2019-12-05T23:42:26.060Z
This is a great piece of work, thankyou.  From NEARGALCAT I see that the catchment dataset contains 869 nearby galaxies which are either < 11 Mpc from us (or with radial velocities <600 kps rel to centre of the Local Group of ~51 galaxies). It is interesting to see how much stronger the Milky Way is than these nearby galaxies.  I would love to see the contribution of the other estimated 2 trillion galaxies in the observable universe. And LCDM posits that 85% of the universe's matter is Dark (not in the form of stars, planets or visible gas):- [Universe](https://en.wikipedia.org/wiki/Universe).
steveOw
CC BY-SA 4.0
2019-12-06T00:19:27.737Z
@steveOw - Thanks! Yes, it would be lovely to have a bigger extragalactic map. The SDSS catalogue is tempting, but only covers parts of the sky and gets incomplete as one goes further out. Overall, in an infinite constant density universe the total force will diverge (as already noted by Newton) since $F=\int_0^\infty G\rho r^2/r^2 dr$ diverges but in the real universe with finite age and GR the picture is a bit better. Integrating that to a distance corresponding to the age of the universe gives $1.1\cdot 10^{-11}$ N/kg (but that limit is handwavy due to nonconstant past expansion).
Anders Sandberg
CC BY-SA 4.0
2019-12-06T01:58:11.903Z
That $1.1*10^{-11}$N/Kg is interestingly large - stronger than nearby stars and local galaxies and only an order of magnitude less than the outer Milky Way.  I came up with a meagre $3*10^{-27}$N/Kg (per 3.6 degree * 3.6 degree pixel) assuming $10^{14}$ solar masses spread evenly through a spherical, static, visible universe of radius $1.4*10^{10}$pc ($4.2*10^{23}$Km). But I could well be in error and will re-check tomorrow!
steveOw
CC BY-SA 4.0
2019-12-06T18:49:00.870Z
I revised my rough estimate of average background Force per target mass per pixel to: $0.86*10^{-14}N/Kg$ which is much closer to your value. (The details are in a new appendix to my original question).
steveOw
CC BY-SA 4.0
2019-12-07T18:56:54.293Z
It doesn't matter how strong the contribution of anything is if it is more-or-less uniformly distributed on the sky. To make a "sky map of gravity", then for each direction you need to make a sum over the whole sky of the scalar product of the thing you have calculated with the radial vector you are interested in. That would be a sky map of gravity.
ProfRob
CC BY-SA 4.0
2019-12-07T20:41:22.313Z
I should add though - you deserve the +100.
ProfRob
CC BY-SA 4.0
2019-12-10T03:29:19.107Z
@RobJeffries I have written some stuff in response to your comments in the [chatroom](https://chat.stackexchange.com/rooms/101931/discussion-between-steveow-and-rob-jeffries).
steveOw
CC BY-SA 4.0
2019-12-05T19:44:39.217Z
2019-12-05T19:44:39.217Z
Anders Sandberg
<p><strong>PREAMBLE</strong></p>

<p>Sky Maps exist which show the distribution and intensity of electromagnetic radiation received on Earth at various wavelengths.  An obvious form is a <a href="https://en.wikipedia.org/wiki/Star_chart" rel="noreferrer">map of the stars</a> (Star Chart). </p>

<p>Maps have also been made for the distribution/intensity of <a href="https://en.wikipedia.org/wiki/Cosmic_microwave_background" rel="noreferrer">Cosmic Microwave Background Radiation</a> and <a href="https://www.chandra.harvard.edu/incl/Skymap/chandra_skymap.html" rel="noreferrer">"to visualize the X-ray Universe"</a>.</p>

<p>I find it interesting to consider the appearance of a Sky Map which showed the variation in intensity of the Gravitational Attraction Force (as currently experienced on Earth) attributable to masses lying in different directions (looking radially outwards from the Earth).</p>

<p>Essentially (if we start with a Newtonian Model) this would be a map which indicated the magnitude of the sum of "mass/distance_squared" for all the masses falling inside each particular pixel of the map.  This then raises the questions of cosmological nature e.g. "how far do gravitational influences extend away from the Earth", along the lines of <a href="https://en.wikipedia.org/wiki/Olbers%27_paradox" rel="noreferrer">Olbers Paradox</a> and "how was matter distributed through the Universe through time" - see <a href="https://en.wikipedia.org/wiki/Cosmological_principle" rel="noreferrer">The Cosmological Principle</a>.</p>

<p><strong>QUESTION</strong></p>

<p>To bring this into the form of a question that can be answered succinctly, my question is: "What would a SkyMap of Radial Gravitational Force, looking outwards from Earth, look like (excluding that due to masses in the Solar System and Milky Way galaxy)."</p>

<p><strong>APPENDIX</strong></p>

<p><strong>Calculation of Average Background Gravitational Force from Visible Universe</strong></p>

<p>(Assuming, for simplicity, a static, spherical, universe of uniformly-distributed ordinary matter with Newtonian gravity and a test mass at the centre of the universe).</p>

<p><span class="math-container">$\rho$</span> Average density of visible matter  <span class="math-container">$=3*10^{-28} Kg/m^{3}$</span></p>

<p><span class="math-container">$R$</span> Radius of visible universe   <span class="math-container">$=1.7*10^{26} m$</span></p>

<p>Note:- values of <span class="math-container">$\rho$</span> and <span class="math-container">$R$</span> are obtained from <a href="https://people.cs.umass.edu/~immerman/stanford/universe.html" rel="noreferrer">here</a>.</p>

<p><span class="math-container">$G$</span> Universal Gravitational Constant  <span class="math-container">$=6.673*10^{-11} Nm^{2}Kg^{-2}$</span></p>

<p><span class="math-container">$f$</span> Fraction of sky represented by <span class="math-container">$3.6 * 3.6$</span> degree rectangular pixel <span class="math-container">$=1/(50*100)$</span> <span class="math-container">$= 0.0002$</span>.</p>

<p>The formula for Force per unit mass <span class="math-container">$F_m$</span> is obtained by integrating (over the interval <span class="math-container">$0,R$</span>) the incremental force per unit mass <span class="math-container">$F_i/m = G.M_i/R_i^2)$</span>.</p>

<p>Here <span class="math-container">$M_i$</span> is the incremental mass contained within the incremental volume <span class="math-container">$V_i$</span> of the spherical shell of radius <span class="math-container">$R_i=R_{i-1}+\text{d}R$</span> and given by <span class="math-container">$M_i = \rho.V_i \approx \rho.4\pi.R_i^2.\text{d}R$</span>.</p>

<p>Thus the total force/unit mass is given by:-</p>

<p><span class="math-container">$$ F_m = G\int_0^R \frac{\rho.4\pi.R_i^2}{R_i^2}\text{d}R  = 4\pi G \rho R $$</span>.</p>

<p>And so the force per unit target mass per pixel is given by:- </p>

<p><span class="math-container">$$F_{m,f} =4\pi G\rho R f  = 0.86*10^{-14} NKg^{-1}.$$</span></p>

CC BY-SA 4.0
2019-12-04T04:13:38.083Z
`+1` for a cool question! Hopefully it averages to zero beyond some radius in spacetime where there's not enough information to model it anyway.
uhoh
CC BY-SA 4.0
2019-12-04T17:19:24.450Z
@barrycarter  In a static Newtonian gravitational universe with homogenous distribution of matter Olbers' theorem indicates that as we move in equal steps outwards radially from Earth the additional mass in an enclosing "catchment volume sphere" of radius $R$ increases $\propto R^3$ while the extra grav force at Earth decreases  $\propto 1/R^2$.  So the force at Earth should increase $\propto R$ as we increase $R$.  I'm interested to see what modern astronomy or cosmology models might have to say about the situation. (I consider the grav forces from the Milky Way and Solar System as known).
steveOw
CC BY-SA 4.0
2019-12-04T21:29:04.680Z
@barrycarter -  Including distant galaxies certainly, but I'm not restricting it to them alone, there could be other significant types/patterns of contributing matter, e.g. from early pre-galactic universe, also Dark Matter.
steveOw
CC BY-SA 4.0
2019-12-05T01:23:12.100Z
[Cosmic Void “Pushes” Milky Way](https://www.skyandtelescope.com/astronomy-news/cosmic-void-pushes-milky-way-3001201723/). This is a cross section, not a sky map. The sky maps are of radiation that radiates from distant sources. Gravitational fields don't work that way.
Keith McClary
CC BY-SA 4.0
2019-12-05T02:04:03.420Z
@Keith McClary - Very interesting article & video thanks.  That model seems consistent with Newtonian gravity to me:- e.g. a region of below-average mass behind me will lead to increase in my acceleration towards a region of above average mass in front of me. Given a map of (static) masses in a 3D volume it is straightforward to compute a Sky Map of the variation in grav attractive force  in different directions acting on a small target mass (e.g. at the centre of the volume). Whereas the resultant acceleration of the target depends on the balance of forces, obtainable by vector summation.
steveOw
CC BY-SA 4.0
2019-12-07T18:41:37.203Z
This makes no sense. Gravity is a vector field. The gravitaitional force directed along a particular line in space is determined by ALL the gravitational sources in the (observable) universe, not just those along that line.
ProfRob
CC BY-SA 4.0
2019-12-07T18:47:23.010Z
It's a sky map of $\sum m_i/r_i^2$, that's all, and certainly not of force.
ProfRob
CC BY-SA 4.0
2019-12-07T18:50:15.370Z
@Rob Jeffries In a Newtonian Model gravitational force operates in a straight line from particle to particle.  The net acceleration of a test particle is the resultant of all the contributing force vectors.
steveOw
CC BY-SA 4.0
2019-12-07T18:52:48.810Z
Yes. Of ALL of them, not just those along a particular line. I agree it is an interesting thing to map, but it isn't a "sky map of gravity".
ProfRob
CC BY-SA 4.0
2019-12-07T18:59:21.703Z
However, it can be done. Just take Anders' maps and sum the scalar product of the thing he has calculated with the radial direction you are interested in and do the sum over the whole sky and for each radial direction. That would be a "sky map of gravity".
ProfRob
CC BY-SA 4.0
2019-12-07T19:00:20.493Z
@RobJeffries It is a SkyMap of "Gravitational Attraction Forces" acting on a test particle. I used "Sky Map of Gravity" in the title for brevity. "Gravity" is an ambiguous entity.  Gravitational Force and Gravitational Acceleration less so.
steveOw
CC BY-SA 4.0
2019-12-07T19:07:42.323Z
@Rob Jeffries I think perhaps that you are presuming a Field model whereas I am assuming a particle-to-particle model as in Newton's Law.
steveOw
CC BY-SA 4.0
2019-12-07T19:10:33.233Z
@steveOw gravity is a vector in Newtonian gravity too. It is also a field. I have no idea what you mean by "Newton's Law". But if you mean the universal gravitation law, that is a vector equation and gravitational fields obey the principle of superposition. There is no reason why you should integrate the force for all sources along a line of sight, but then not include the (resolved) force component along that line for all other sources.
ProfRob
CC BY-SA 4.0
2019-12-07T21:33:24.927Z
Let us [continue this discussion in chat](https://chat.stackexchange.com/rooms/101931/discussion-between-steveow-and-rob-jeffries).
steveOw
<gravity><cosmology>
A Sky Map of Gravity - what would it look like?
steveOw
# A Sky Map of Gravity - what would it look like?

(Asked by: username_0 on 2019-12-04)

**PREAMBLE**

Sky Maps exist which show the distribution and intensity of electromagnetic radiation received on Earth at various wavelengths. An obvious form is a [map of the stars](https://en.wikipedia.org/wiki/Star_chart) (Star Chart).

Maps have also been made for the distribution/intensity of [Cosmic Microwave Background Radiation](https://en.wikipedia.org/wiki/Cosmic_microwave_background) and ["to visualize the X-ray Universe"](https://www.chandra.harvard.edu/incl/Skymap/chandra_skymap.html).

I find it interesting to consider the appearance of a Sky Map which showed the variation in intensity of the Gravitational Attraction Force (as currently experienced on Earth) attributable to masses lying in different directions (looking radially outwards from the Earth).

Essentially (if we start with a Newtonian Model) this would be a map which indicated the magnitude of the sum of "mass/distance\_squared" for all the masses falling inside each particular pixel of the map. This then raises the questions of cosmological nature e.g. "how far do gravitational influences extend away from the Earth", along the lines of [Olbers Paradox](https://en.wikipedia.org/wiki/Olbers%27_paradox) and "how was matter distributed through the Universe through time" - see [The Cosmological Principle](https://en.wikipedia.org/wiki/Cosmological_principle).

**QUESTION**

To bring this into the form of a question that can be answered succinctly, my question is: "What would a SkyMap of Radial Gravitational Force, looking outwards from Earth, look like (excluding that due to masses in the Solar System and Milky Way galaxy)."

**APPENDIX**

**Calculation of Average Background Gravitational Force from Visible Universe**

(Assuming, for simplicity, a static, spherical, universe of uniformly-distributed ordinary matter with Newtonian gravity and a test mass at the centre of the universe).

$\rho$ Average density of visible matter $=3\*10^{-28} Kg/m^{3}$

$R$ Radius of visible universe $=1.7\*10^{26} m$

Note:- values of $\rho$ and $R$ are obtained from [here](https://people.cs.umass.edu/~immerman/stanford/universe.html).

$G$ Universal Gravitational Constant $=6.673\*10^{-11} Nm^{2}Kg^{-2}$

$f$ Fraction of sky represented by $3.6 \* 3.6$ degree rectangular pixel $=1/(50\*100)$ $= 0.0002$.

The formula for Force per unit mass $F\_m$ is obtained by integrating (over the interval $0,R$) the incremental force per unit mass $F\_i/m = G.M\_i/R\_i^2)$.

Here $M\_i$ is the incremental mass contained within the incremental volume $V\_i$ of the spherical shell of radius $R\_i=R\_{i-1}+\text{d}R$ and given by $M\_i = \rho.V\_i \approx \rho.4\pi.R\_i^2.\text{d}R$.

Thus the total force/unit mass is given by:-

$$ F\_m = G\int\_0^R \frac{\rho.4\pi.R\_i^2}{R\_i^2}\text{d}R = 4\pi G \rho R $$.

And so the force per unit target mass per pixel is given by:-

$$F\_{m,f} =4\pi G\rho R f = 0.86\*10^{-14} NKg^{-1}.$$

**Tags:** gravity, cosmology
> **Comments:**
> **username_1**: `+1` for a cool question! Hopefully it averages to zero beyond some radius in spacetime where there's not enough information to model it anyway.  
> **username_0**: @barrycarter  In a static Newtonian gravitational universe with homogenous distribution of matter Olbers' theorem indicates that as we move in equal steps outwards radially from Earth the additional mass in an enclosing "catchment volume sphere" of radius $R$ increases $\propto R^3$ while the extra grav force at Earth decreases  $\propto 1/R^2$.  So the force at Earth should increase $\propto R$ as we increase $R$.  I'm interested to see what modern astronomy or cosmology models might have to say about the situation. (I consider the grav forces from the Milky Way and Solar System as known).  
> **username_0**: @barrycarter -  Including distant galaxies certainly, but I'm not restricting it to them alone, there could be other significant types/patterns of contributing matter, e.g. from early pre-galactic universe, also Dark Matter.  
> **username_2**: [Cosmic Void “Pushes” Milky Way](https://www.skyandtelescope.com/astronomy-news/cosmic-void-pushes-milky-way-3001201723/). This is a cross section, not a sky map. The sky maps are of radiation that radiates from distant sources. Gravitational fields don't work that way.  
> **username_0**: @Keith McClary - Very interesting article & video thanks.  That model seems consistent with Newtonian gravity to me:- e.g. a region of below-average mass behind me will lead to increase in my acceleration towards a region of above average mass in front of me. Given a map of (static) masses in a 3D volume it is straightforward to compute a Sky Map of the variation in grav attractive force  in different directions acting on a small target mass (e.g. at the centre of the volume). Whereas the resultant acceleration of the target depends on the balance of forces, obtainable by vector summation.  
> **username_3**: This makes no sense. Gravity is a vector field. The gravitaitional force directed along a particular line in space is determined by ALL the gravitational sources in the (observable) universe, not just those along that line.  
> **username_3**: It's a sky map of $\sum m_i/r_i^2$, that's all, and certainly not of force.  
> **username_0**: @Rob Jeffries In a Newtonian Model gravitational force operates in a straight line from particle to particle.  The net acceleration of a test particle is the resultant of all the contributing force vectors.  
> **username_3**: Yes. Of ALL of them, not just those along a particular line. I agree it is an interesting thing to map, but it isn't a "sky map of gravity".  
> **username_3**: However, it can be done. Just take Anders' maps and sum the scalar product of the thing he has calculated with the radial direction you are interested in and do the sum over the whole sky and for each radial direction. That would be a "sky map of gravity".  
> **username_0**: @RobJeffries It is a SkyMap of "Gravitational Attraction Forces" acting on a test particle. I used "Sky Map of Gravity" in the title for brevity. "Gravity" is an ambiguous entity.  Gravitational Force and Gravitational Acceleration less so.  
> **username_0**: @Rob Jeffries I think perhaps that you are presuming a Field model whereas I am assuming a particle-to-particle model as in Newton's Law.  
> **username_3**: @steveOw gravity is a vector in Newtonian gravity too. It is also a field. I have no idea what you mean by "Newton's Law". But if you mean the universal gravitation law, that is a vector equation and gravitational fields obey the principle of superposition. There is no reason why you should integrate the force for all sources along a line of sight, but then not include the (resolved) force component along that line for all other sources.  
> **username_0**: Let us [continue this discussion in chat](https://chat.stackexchange.com/rooms/101931/discussion-between-steveow-and-rob-jeffries).  

---

## 1 answers

### Answer by username_4 on 2019-12-05. Score: 16 | Accepted answer

The rough answer is: just like the sun makes it hard or impossible to see planets and stars during the day, it dominates the gravity sky.

[![Linear plot of the gravity sky](https://i.sstatic.net/p4C99.png)](https://i.sstatic.net/p4C99.png)

But there are interesting patterns there if we view the sky using a logarithmic scale.

[![Log plot of the gravity sky](https://i.sstatic.net/UkCuN.png)](https://i.sstatic.net/UkCuN.png)

The above plots are a combination of the gravitational force of the sun, moon and planets, the stars in the Hipparchos catalog, the gravitational field of the Milky Way, and the nearby galaxies. Each component have a somewhat different effect.

**The solar system**

The force from the sun is shockingly big: 0.0059 N (per kilogram - I will omit that henceforth). At first this might appear impossible since we never notice the force, but remember that the setup for the question is a fixed point in space rather than the free-falling orbit of Earth. Because of the orbit this force exactly cancels.

Note that for spherical objects we only need to calculate the force from the direction of their centre of mass, not for their full extent. Hence regardless of what resolution I use for my map planets and stars will only change one box.

The force from the moon is $3.3198\cdot 10^{-5}$ N, a factor of 179 less. (Why do we experience tides from the moon more than from the sun? Because the gradient for the gravitational field is steeper for the nearby moon than for the remote sun). The other planets have much less effect, from $10^{-6.8}$ N for Jupiter down to $10^{-9.5}$ N for Neptune (and Pluto at $10^{-13.4}$).

For positions I used the Matlab Aerospace toolbox planetEphemeris function. Incidentally, I might well have messed up some of the many different coordinate frames here, so don't use these maps to navigate. The important thing for this question is the magnitudes of the forces anyway.

**Stars**

I added the field of the stars in the Hipparchos catalog, but they do not show up even on my logarithmic map because they are overpowered by the galaxy at large. If we just plot their contributions to different boxes we get the following map:

[![Log plot of the Hipparchos map](https://i.sstatic.net/KLlQy.png)](https://i.sstatic.net/KLlQy.png)

Note that the peak force in each box is on the order of $10^{-14}$ N: even Pluto matters more than nearby stars. This map also shows an annoying problem: if we just sum the contributions in each 3.6 times 3.6 degree rectangle in the coordinate system the polar rectangles will be narrow and have smaller area; hence there are fewer stars included and less force. Dividing by $\cos(\delta)$ removes this bias and shows that the stars indeed pull on us roughly equally in all directions.

(I used a rough estimate of their masses: assuming them all to be on the main sequence I calculated their luminosity from their absolute magnitude, and then estimated the mass as $M\approx M\_\odot (L/L\_\odot)^{1/3.5}$.)

**The Milky Way**

[![Log picture of the galaxy force](https://i.sstatic.net/JOJXU.png)](https://i.sstatic.net/JOJXU.png)

I modelled the Milky Way using the approximate values and formulas in Carroll & Ostlie (p. 884), essentially a neutral gas disk, a thin disk, a thick disk, a central bulge, a stellar halo plus a dark matter halo. The disks for example have density scaling as $$\rho(r,z)=\rho\_0 e^{-r/h\_r}e^{-|z|/h\_z}$$ while the dark matter halo has the [Navarro-Frenck-White profile](https://en.wikipedia.org/wiki/Navarro%E2%80%93Frenk%E2%80%93White_profile) $$\rho(r)=\frac{\rho\_0}{\frac{r}{R\_s}\left(1+\frac{r}{R\_s}\right)^2}.$$ I integrate these functions along a ray extending from the part of the sky I am interested in. Note that the integral of the force will be $$F=\int\_0^\infty \frac{G\rho(\mathbf{x}\_0 + \mathbf{d} r)}{r^2} r^2 dr = G \int\_0^\infty \rho(\mathbf{x}\_0 + \mathbf{d} r) dr$$ since it is really an integral of an ever widening pyramid of area $r^2$: the total force is only dependent on the total amount of mass there. Since these density fields change slowly I ignore lateral variations. Again there should be a correction for declination, but it does not make a great difference in the picture.

The end result is a fair amount of force, especially along the galactic plane. This causes "galactic tides" in stellar orbits and can sometimes rival the force from remote planets in the solar system. Again, we do not notice most of the gravitational field because the sun is free-falling in the galaxy.

Sag A\* does not show up in this image because it is not very massive compared to the full galaxy. The force is just $9.3678\cdot 10^{-15}$ N.

**Nearby galaxies**

[![Log near galaxy force](https://i.sstatic.net/wsFKN.png)](https://i.sstatic.net/wsFKN.png)

I used the [NEARGALCAT catalogue](https://heasarc.gsfc.nasa.gov/W3Browse/galaxy-catalog/neargalcat.html) to estimate the effect of nearby galaxies. This is in many ways the most dodgy estimate of all. While we know positions well we have weak data on actual masses. I used the estimated masses whenever they existed, and a linear fit between angular extent and mass derived from those to impute masses for galaxies lacking values. I am a bit uncertain if the mass estimates are just photometric or include dark matter; if not including DM, they should be an order of magnitude higher.

Typical forces are on the order of $10^{-19}$ N, although heavier galaxies might pull us with $10^{-15}$ N. They are hence on par with the nearby stars: totally overpowered by the galactic gravitational field and the solar system.

Technically, since galaxies are extended non-spherical objects, I should integrate over areas here. But there are limits to my energy.

In short, the gravity sky has nontrivial patterns if one has gravity-eyes with a logarithmic response.

> **Comments:**
> **username_5**: +1. I've often daydreamed about what an 'incident graviton flux' map would look like. This answers that brilliantly.  
> **username_0**: This is a great piece of work, thankyou.  From NEARGALCAT I see that the catchment dataset contains 869 nearby galaxies which are either < 11 Mpc from us (or with radial velocities <600 kps rel to centre of the Local Group of ~51 galaxies). It is interesting to see how much stronger the Milky Way is than these nearby galaxies.  I would love to see the contribution of the other estimated 2 trillion galaxies in the observable universe. And LCDM posits that 85% of the universe's matter is Dark (not in the form of stars, planets or visible gas):- [Universe](https://en.wikipedia.org/wiki/Universe).  
> **username_4**: @steveOw - Thanks! Yes, it would be lovely to have a bigger extragalactic map. The SDSS catalogue is tempting, but only covers parts of the sky and gets incomplete as one goes further out. Overall, in an infinite constant density universe the total force will diverge (as already noted by Newton) since $F=\int_0^\infty G\rho r^2/r^2 dr$ diverges but in the real universe with finite age and GR the picture is a bit better. Integrating that to a distance corresponding to the age of the universe gives $1.1\cdot 10^{-11}$ N/kg (but that limit is handwavy due to nonconstant past expansion).  
> **username_0**: That $1.1*10^{-11}$N/Kg is interestingly large - stronger than nearby stars and local galaxies and only an order of magnitude less than the outer Milky Way.  I came up with a meagre $3*10^{-27}$N/Kg (per 3.6 degree * 3.6 degree pixel) assuming $10^{14}$ solar masses spread evenly through a spherical, static, visible universe of radius $1.4*10^{10}$pc ($4.2*10^{23}$Km). But I could well be in error and will re-check tomorrow!  
> **username_0**: I revised my rough estimate of average background Force per target mass per pixel to: $0.86*10^{-14}N/Kg$ which is much closer to your value. (The details are in a new appendix to my original question).  
> **username_3**: It doesn't matter how strong the contribution of anything is if it is more-or-less uniformly distributed on the sky. To make a "sky map of gravity", then for each direction you need to make a sum over the whole sky of the scalar product of the thing you have calculated with the radial vector you are interested in. That would be a sky map of gravity.  
> **username_3**: I should add though - you deserve the +100.  
> **username_0**: @RobJeffries I have written some stuff in response to your comments in the [chatroom](https://chat.stackexchange.com/rooms/101931/discussion-between-steveow-and-rob-jeffries).  

---

