So far, we've examined the methods by which several properties of stars can be measured: distance, luminosity, temperature, size. Another fundamental property of a star, or any celestial object, is its mass. How can one measure the mass of a star? It turns out that one must find a star which is in orbit around another star(s) and use gravity as a tool to turn orbital motion into mass.
But it's a complicated subject. Let's start off with a simple example of orbital motion -- planets in our own solar system -- and work our way up to the more complex cases of distant stars.
In 1609, Kepler published a book, Nova Astronomica, in which he made a revolutionary claim: planets don't move in CIRCLES, as everyone had previously thought (including Copernicus); instead,
Kepler's First Law: the planets move in ellipses, with the Sun at one focus.
Let's pause to consider the properties of an ellipse. There are several ways to define this shape:
Below is a sample ellipse drawn on a background grid. Mark on your piece of paper the following quantities; make all measurements in units of the grid spacing.
Kepler's Second Law: A line connecting a planet to the Sun sweeps out equal areas in equal times.
Use this rule on your paper:
Note that Kepler's Second Law is nice, but it deals only with RELATIVE speeds. We can't use it directly to figure out exactly how far a planet will move over a period of, say, three weeks. Is there any way to figure out exactly where a planet will be at some particular time?
Yes -- two ways, in fact:
Thanks to fast, cheap computers, we members of the modern world can simply plug some initial conditions into a computer program and integrate the motion of the bodies numerically to follow their motion as a function of time. One of the problems on this week's homework asks you to do just that.
Back in the old days, however, computing wasn't quite so cheap and easy. Kepler didn't have computers or calculators. In fact, Kepler didn't even have LOGARITHMS; Napier published his first tables in 1632, after Kepler had derived his laws of planetary motion. Kepler did manage to devise a method of computing the position of a planet at any time. It involved a bit of calculating, but nowhere near as much as the modern brute-force approach. Since it's an interesting little mathematical puzzle, let's look at it in some detail.
The problem is this: we know the orbital parameters of a planet's motion around the Sun: period P, semimajor axis a, eccentricity e. We also know the time T when the planet reaches its perihelion passage. Where will the planet be in its orbit at some later time t?
If the orbit is circular, then this is easy: the fraction of a complete orbit is equal to the fraction of a complete period which has elapsed since the last perihelion passage. For example, if (t - T) is exactly one-quarter of the period P, then the planet will have made exactly one-quarter of a full circle around the Sun.
But, as Kepler's Second Law states, planets in elliptical orbits do NOT move with a constant speed, nor with a constant angular speed. In real life, some planetary orbits are significantly non-circular, so the circular approximation won't work. What can we do?
We can measure the position of a planet in its elliptical orbit with the angle between its radius vector and the perihelion position. This angle is called the true anomaly, and is conventionally written as the letter v.
Yes, I've moved the principal focus closer to the center of the circle than it should be, for clarity.
What is v as a function of time t?
Kepler found an answer to this question, but it required a bit of a roundabout journey. He discovered that if he defined a "auxiliary" quantity, he could solve for that quantity as a function of time; and then he could convert from the "auxiliary" quantity to the desired true anomaly v.
Kepler's first step was to draw a circle around the ellipse, and project the position of the planet on its elliptical orbit upwards to meet the circle.
The angle E measured from perihelion position, to center of circle, to projected position of planet, is called the eccentric anomaly. If one can find E, one can go back to the desired variables r and v like so:
Kepler's next step was to find a mathematical relationship between this eccentric anomaly E and time. He computed n, the average angular speed of the planet, also called the mean motion,
You should express the mean motion, and all angular quantities involved in Kepler's equation, in radians.
Finally, he was able to write the following, known as Kepler's equation:
Looks good, right? You have to go through several steps:
Actually, it is. Step 3, solving Kepler's equation for E given some time difference (t - T), turns out to be a doozy. here is no simple closed-form solution, in general, though there are series solutions which converge quickly for small values of the eccentricity. Astronomer Friedrich Bessel devised his eponymous "Bessel Functions" in 1817 as a means to solve Kepler's equation. Letting M = (t - T) n (this is known as the "mean anomaly") one can write
where Jk is the k'th Bessel function of the first kind.
Fortunately, even though there isn't a closed solution, Kepler's equation yields pretty quickly to any number of numerical attacks, especially when the eccentricity e is small. You can program a computer to use Newton's method or the the bisection technique. You can even fall back on trial and error. Give it a try.
Exercise: Your piece of paper shows an elliptical orbit. Suppose that the orbit is P = 500 days in a counterclockwise direction. Where will the planet be at (t - T) = 400 days after perihelion passage? Calculate the true anomaly angle v and use it to mark the position of the planet along the orbit.
Start by finding the mean motion n and the mean anomaly M = n(t - T). Use a starting guess that the eccentric anomaly E is equal to the mean anomaly. Plug that E into the left-hand side of Kepler's equation and see what you get. If you get exactly M, you're done! Otherwise, modify your value of E and try again. See if you can find E to four significant figures.
Copyright © Michael Richmond. This work is licensed under a Creative Commons License.