WARNING: MATHS CONTENT!!!

One of the things I had to do as part of the rework of my novel The Tau Ceti Diversion, is to try and work out the surface gravity of my fictional planets. From the Kepler data, there are two exoplanets located in the Tau Ceti system that are likely to be in the system’s habitable zone, or where there is the possibility of liquid water on the surface, and perhaps life as we know it.

To play around with my estimates of gravity, I used ratioed rearrangements of Newton’s law of gravity (law of universal gravitation) and a simple formula relating the density of a spherical planet to its mass and radius (these are at the bottom of the post in the ADDENDUM).

Here’s Newtons famous law:)

The two planets thought to be in Tau Ceti’s habitable zone are denoted Tau Ceti e and Tau Ceti f. What is known about these two planets is their likely orbit, eccentricity, and their mass. All of these properties have been derived by calculation, based on observed data, so are all known to within appropriate error bounds, but I’m leaving the error off my scribblings so things don’t get too messy.

Tau Ceti e is thought to be around 4.3 Earth Masses, or Me (i.e. 4.3 times as heavy as Earth), while Tau Ceti f, the planet that orbits a bit further out, is thought to be around 6.67 Me. For the astronomically minded, these two planets orbit at around 0.55 and 1.35 AU from Tau Ceti respectively.

So, here’s where I cheated a bit, like any good engineer. I started with the answer I wanted and calculated backwards to see if the answer I wanted led to reasonable base assumptions. This is not as cheeky as it sounds, because when you have an insoluble problem (i.e. not enough data is known for an explicit result), an iterative approach is often used.

For my story to work, I needed a surface gravity on my planet of no more than 1.2g – that’s twenty percent higher than Earth’s. But how could I get a gravity that low on a planet that was over 4 times the mass of Earth? The answer is that surface gravity is a function of mass and radius, or going a step further along the calculation path, mass and density.

I used a ratioed form of Newton’s law that allowed me to relate the ratio of two planets gravitational forces to the ratios of their masses and radii. I already knew the ratio of the gravities ( assumed at gTCe/gE= 1.2) and the ratio of the masses (MTCe/ME = 4.3), so could calculate the ratio of radii (rE/rTCe) at 1.89. Using another formula that related the ratio of the two planet’s densities to their ratioed mass and radii, I could then calculate their ratioed densities (dens TCe/ densE) at 0.63. So at the end of all that, to have a surface gravity of 1.2 g, Tau Ceti e would have to have a density of 63% of Earth’s. Is that reasonable?

The density of Earth is 5.514 g/cm3, not too much different from the density of a rocky planet like Mercury (5.427 g/cm3), but a lot higher than other solar system planets like Jupiter and Uranus (1.326 g/cm3 and 1.27 g/cm3 respectively), comprised of lighter materials. A surface gravity of 1.2g on Tau Ceti e would put its density at around 3.5 g/cm3, less dense than our own rocky planets, but certainly in a feasible range.

So what sort of densities would you expect for the Tau Ceti system? One clue is the metallicity of the system, which is a measure of the ratio of iron to hydrogen in the star’s makeup. In the case of Tau Ceti, this is estimated to be around one third of our own sun. This indicates the star is likely to be older than the Sun, made up of stellar remnants left over from less evolved stars that have not had time to form as much of the heavier elements in their internal fusion factories.

So Tau Ceti is made up of lighter elements. Based on this, it was reasonable to assume that the planets in the Tau Ceti system would also be made up of proportionally lighter elements, and quite possibly in the range I had estimated. Tau Ceti e and Tau Ceti f are also large planets – much larger than our own Earth – so having a density in between Earth and our own gas giants also made sense to me.

Using the same planetary density I had calculated for Tau Ceti e, for the larger Tau Ceti f, gave me a surface density of around 1.4g for the bigger planet – just a little too high for feasible human colonisation – and that fit nicely with my story as well.

It was a lot of fun playing with these calculations, and thankfully the known science fit with my story, at least with some comfortable wiggle room!

The Tau Ceti Diversion is due to be launched on September 1^{st} 2016! Read more about what happens in the story here!

ADDENDUM

For those interested in the maths. . .

Density formula: densp= Mp / (4/3*pi()*rp^3)

Where:

densp= Density of Planet (kg/m3)

Mp = mass of planet (kg)

rp = radius of planet (m)

In ratio form: densp1/densp2= Mp1/Mp2 *(rp2/rp1)^3

Ratio of Newtons law relating gravity, mass and radius of two planets:

gp1/gp2= Mp1/Mp2 *(rp2/rp1)^2