1 What problem does the eccentricity vector solve?

For ideal two-body motion,

(1)

with

(2)

for a satellite whose mass is negligible relative to the central body.

Angular momentum immediately gives one conserved vector,

(3)

where h is the specific angular momentum. Because the gravitational acceleration is parallel to r,

(4)

Thus h is constant.

This tells us something extremely important: h fixes the orbital plane and the direction of positive orbital rotation. But it does not tell us where periapsis lies inside that plane. Rotating an ellipse about its focus within the same plane leaves h unchanged.

We therefore seek a second conserved vector with two properties:

1.

it should lie in the orbital plane rather than normal to it;

2.

it should identify the preferred direction of the conic, which for an ellipse is the periapsis direction.

The eccentricity vector is precisely that object.

The distinction in Figure 1 is important. The magnitude

(5)

is dimensionless. It is not the distance from the focus to periapsis. For an ellipse that physical distance is

(6)

2 The “discovery” idea: build another conserved vector

The usual formula looks mysterious if we start at the end,

(7)

A more useful way to understand it is to ask how someone could have been led to this combination.

The angular-momentum vector h is perpendicular to the orbital plane. The velocity v lies in the orbital plane. Therefore

(8)

is itself in the orbital plane. That makes it a natural candidate for constructing an in-plane conserved vector.

The question becomes:

What is the time derivative of v × h under an inverse-square central force?

If that derivative can be written as the derivative of some other simple vector, the difference between the two quantities will be constant.

This is exactly what happens.

3 A needed kinematic identity: radial distance and radial direction

Before differentiating v × h, we need two facts about the position vector. Write

(9)

Differentiating gives

(10)

Because r has unit magnitude,

(11)

Differentiating this identity gives

(12)

so

(13)

Thus the derivative of a unit vector is perpendicular to the unit vector itself: changes direction, not magnitude.

Dot Eq. (10) with r:

r ⋅ v	= ṙ (r ⋅r) + r(r ⋅)	(14)
	= ṙ.	(15)

Therefore

(16)

This says that the scalar radial rate is simply the component of velocity along the radial unit vector.

Solving Eq. (10) for the unit-vector derivative,

(17)

Only the transverse part of velocity rotates r.

4 The key cross-product identity

Now evaluate

(18)

Using h = r × v,

(19)

Apply the vector triple-product identity

(20)

Then

r × h	= r(r ⋅ v) − v(r ⋅ r)	(21)
	= (rr)ṙ − rv	(22)
	= −r .	(23)

Using Eq. (17),

(24)

so

(25)

This identity is the algebraic hinge of the eccentricity-vector derivation.

5 Deriving the conserved vector from Newton’s law

Differentiate v × h:

(v × h)

= × h + v ×.

(26)

Because = 0,

(27)

Newton’s inverse-square acceleration is

(28)

Therefore

a × h

= −(r × h).

(29)

Substitute Eq. (25):

a × h	= −	(30)
	= μ.	(31)

Thus

(32)

Now the construction is no longer mysterious. Move the right-hand derivative to the left:

(33)

or

(34)

Therefore the quantity in parentheses is a conserved vector.

Divide by the constant μ and define

(35)

This is the eccentricity vector.

The dimensional check is also revealing. Since

(36)

and

(37)

the ratio (v × h)∕μ is dimensionless, just like r. Therefore e is dimensionless.

6 Why does this vector point toward periapsis?

Conservation alone does not yet prove that e points toward periapsis. That fact comes from projecting Eq. (35) onto the radial direction.

Let α be the instantaneous angle between e and r. Then

(38)

Dot Eq. (35) with r:

(39)

Use the scalar triple-product identity,

(40)

But

(41)

so

(42)

Therefore

(43)

Equation (39) becomes

(44)

Rearrange:

(45)

and solve for r:

(46)

Define the semilatus rectum

(47)

Then

(48)

This is the polar equation of a conic with the central body at the focus.

For e > 0, r is smallest when the denominator is largest, which occurs at

(49)

Thus the direction α = 0, namely the direction of e, is the direction of closest approach. That is periapsis.

Therefore

(50)

Once that fact is established, the angle α is conventionally called the true anomaly ν, so Eq. (48) becomes

(51)

This is one of the most useful results in celestial mechanics: the same vector derived from Newton’s differential equation simultaneously encodes the orbit’s shape and its orientation within the orbital plane.

7 Why is its magnitude the conic eccentricity?

The symbol e in Eq. (51) is exactly the geometric eccentricity parameter of the conic. Since that equation was obtained using

(52)

the magnitude of the eccentricity vector is the scalar conic eccentricity.

Its value classifies the Kepler orbit:

(53)

A circle is the special elliptical case e = 0.

This circular limit explains an important singularity in classical orbital elements. If

(54)

then the argument of periapsis ω cannot be defined uniquely because there is no distinguished periapsis direction.

8 Connection with orbital energy

The eccentricity magnitude can also be related directly to the specific mechanical energy

(55)

square Eq. (35):

e2	= 2	(56)
	= + 1 −.	(57)

Because v lies in the orbital plane while h is normal to it,

(58)

so

(59)

We already found

(60)

Hence

(61)

Factor h²∕μ²:

(62)

Using Eq. (55),

(63)

so

(64)

This formula connects the two principal invariants of Kepler motion: angular momentum controls the orbital plane and transverse scale, while energy controls whether the conic is bound or unbound. Their combination determines eccentricity.

For an ellipse,

(65)

Substituting into Eq. (64),

(66)

Therefore

(67)

and with Eq. (47),

(68)

Thus the semilatus rectum is where Newtonian dynamics, conic geometry, and classical orbital elements meet.

9 An equivalent Cartesian formula

Equation (35) is compact, but software often uses an equivalent expression containing only r, v, and scalar dot products. Start with

(69)

Apply the triple-product identity:

(70)

Hence

(71)

Since

(72)

Eq. (35) becomes

(73)

Both forms are mathematically identical:

(74)

The cross-product form is often better for geometric understanding; the expanded form makes the dependence on radial velocity r ⋅ v = rṙ explicit.

10 A physical reading of the two terms

The eccentricity vector can be read as a competition between

(75)

The radial unit vector always has magnitude one and continually rotates as the body moves around the focus. The first term also changes as the velocity direction changes. For an inverse-square force, Eq. (32) says these two changes are exactly synchronized:

(76)

Their difference is therefore frozen in inertial space.

That is the deepest intuition behind Eq. (35): it is not a lucky algebraic combination. It is a cancellation created by the special 1∕r² central force.

If the force law is perturbed, for example by Earth’s oblateness, third-body gravity, or atmospheric drag, the osculating eccentricity vector generally changes slowly with time. Its changing direction describes apsidal precession, and its changing magnitude describes changes in orbital shape.

11 The Laplace–Runge–Lenz vector

A closely related conserved quantity is

(77)

This is commonly called the Laplace–Runge–Lenz vector, or simply the Runge–Lenz vector in many mechanics texts. Comparing Eqs. (35) and (77),

(78)

The two vectors therefore have exactly the same direction. The normalized form e is particularly convenient in astrodynamics because its magnitude is immediately the dimensionless orbital eccentricity.

12 Worked GPS-like numerical example

Consider the Earth-centered inertial state

(79)

(80)

with

(81)

First compute the specific angular momentum,

(82)

which gives approximately

(83)

with

(84)

The position magnitude is

(85)

so

(86)

Now evaluate

(87)

The result is

(88)

Its magnitude is

(89)

The orbit is therefore a low-eccentricity ellipse.

The specific mechanical energy is

(90)

which corresponds to

(91)

The semilatus rectum obtained dynamically is

(92)

while the orbital-element expression gives

(93)

The agreement provides a useful implementation check.

The angle between e and r is

(94)

which recovers the true anomaly used to generate the state.

Finally, the radial speed is

(95)

The positive sign indicates that the spacecraft is moving outward, away from periapsis toward apoapsis.

13 Implementation recipe from a Cartesian state

Given r, v, and μ, the eccentricity-vector calculation is short:

r	= ∥r∥,	(96)
h	= r × v,	(97)
e	= −,	(98)
e	= ∥e∥.	(99)

For e > 0, the unit periapsis direction is

(100)

This is exactly the P axis of the perifocal PQW frame.

In numerical software, one should treat very small e carefully. As e → 0, the magnitude of e becomes small and its normalized direction becomes sensitive to numerical noise because a circular orbit has no physically unique periapsis.

14 What to remember

The eccentricity vector can be understood through a compact chain of ideas:

(101)

Its final form is

(102)

Its magnitude is

(103)

its direction is toward periapsis, and its projection onto r immediately generates the orbit equation

(104)

For an ellipse,

(105)

and energy supplies the equivalent relation

(106)

The vector therefore compresses an unusually large amount of orbital information into one conserved quantity: the orbit’s eccentricity, periapsis direction, conic equation, and connection between angular momentum and energy.

References

References