1 Definition

In ideal Newtonian two-body motion, the eccentricity vector is the conserved, dimensionless vector

(1)

where

(2)

Here r is the position of the orbiting body relative to the central body, v = dr∕dt is its relative velocity, h is the specific angular-momentum vector, and μ = G(M + m) ≈ GM is the gravitational parameter.

The eccentricity vector has two defining geometric properties:

(3)

and, for e > 0,

(4)

Thus one vector contains both the shape and an important part of the orientation of a Keplerian orbit. Its magnitude e is the ordinary scalar eccentricity of the conic, while its direction fixes the periapsis line within the orbital plane.

The vector is also the dimensionless form of the Laplace–Runge–Lenz vector. If

(5)

then

(6)

A useful distinction is

(7)

where r_p is the physical focus-to-periapsis distance and c is the center-to-focus distance of an ellipse. These three quantities should not be confused: only e is the magnitude of the eccentricity vector.

For a circular orbit,

(8)

A circle has no unique periapsis, so the zero vector correctly carries no preferred in-plane direction.

2 Why introduce the eccentricity vector?

The governing equation for ideal two-body motion is

(9)

One conserved vector follows immediately:

(10)

Because the gravitational acceleration is parallel to r,

(11)

Thus h is constant.

The angular-momentum vector tells us the orbital plane and its positive normal direction, but it does not identify where periapsis lies inside that plane. Rotating an ellipse about its focus without changing its plane leaves h unchanged. A second conserved vector is therefore useful: one that lies in the orbital plane and singles out the periapsis direction. Equation (1) is that vector.

The remaining question is not what the eccentricity vector is, but why Newton’s inverse-square law produces exactly this combination. The following derivation answers that question.

3 Derivation: searching for an in-plane conserved vector

The formula in Eq. (1) looks arbitrary if it is simply quoted. A more instructive route is to search for a second conserved quantity.

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

(12)

is itself in the orbital plane. This makes it a natural candidate from which to construct a vector that could identify a preferred in-plane direction.

The key question is:

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

If that derivative matches the derivative of another simple in-plane vector, their difference will be conserved. That is exactly what happens.

4 A needed kinematic identity: radial distance and radial direction

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

(13)

Differentiating gives

(14)

Because r has unit magnitude,

(15)

Differentiating this identity gives

(16)

so

(17)

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

Dot Eq. (14) with r:

r ⋅ v	= ṙ (r ⋅r) + r(r ⋅)	(18)
	= ṙ.	(19)

Therefore

(20)

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

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

(21)

Only the transverse part of velocity rotates r.

5 The key cross-product identity

Now evaluate

(22)

Using h = r × v,

(23)

Apply the vector triple-product identity

(24)

Then

r × h	= r(r ⋅ v) − v(r ⋅ r)	(25)
	= (rr)ṙ − rv	(26)
	= −r .	(27)

Using Eq. (21),

(28)

so

(29)

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

6 Deriving the conserved vector from Newton’s law

Differentiate v × h:

(v × h)

= × h + v ×.

(30)

Because = 0,

(31)

Newton’s inverse-square acceleration is

(32)

Therefore

a × h

= −(r × h).

(33)

Substitute Eq. (29):

a × h	= −	(34)
	= μ.	(35)

Thus

(36)

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

(37)

or

(38)

Therefore the quantity in parentheses is a conserved vector.

Divide by the constant μ. The conserved dimensionless vector is therefore

(39)

This is the eccentricity vector.

The dimensional check is also revealing. Since

(40)

and

(41)

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

7 Why does this vector point toward periapsis?

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

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

(42)

Dot Eq. (39) with r:

(43)

Use the scalar triple-product identity,

(44)

But

(45)

so

(46)

Therefore

(47)

Equation (43) becomes

(48)

Rearrange:

(49)

and solve for r:

(50)

The standard polar equation of a conic with a focus at the origin is written

(51)

where p is the semilatus rectum. Comparing this geometric conic form with Eq. (50) identifies

(52)

Thus p = h²∕μ is not an arbitrary orbital-mechanics definition; it is the dynamical value of the conic’s geometrically defined semilatus rectum. Therefore

(53)

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

(54)

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

Therefore

(55)

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

(56)

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.

8 Why is its magnitude the conic eccentricity?

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

(57)

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

Its value classifies the Kepler orbit:

(58)

A circle is the special elliptical case e = 0.

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

(59)

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

9 Connection with orbital energy

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

(60)

square Eq. (39):

e2	= 2	(61)
	= + 1 −.	(62)

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

(63)

so

(64)

We already found

(65)

Hence

(66)

Factor h²∕μ²:

(67)

Using Eq. (60),

(68)

so

(69)

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,

(70)

Substituting into Eq. (69),

(71)

Therefore

(72)

and with Eq. (52),

(73)

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

10 An equivalent Cartesian formula

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

(74)

Apply the triple-product identity:

(75)

Hence

(76)

Since

(77)

Eq. (39) becomes

(78)

Both forms are mathematically identical:

(79)

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

11 A physical reading of the two terms

The eccentricity vector can be read as a competition between

(80)

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. (36) says these two changes are exactly synchronized:

(81)

Their difference is therefore frozen in inertial space.

That is the deepest intuition behind Eq. (39): 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.

12 The Laplace–Runge–Lenz vector

A closely related conserved quantity is

(82)

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

(83)

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.

13 Worked GPS-like numerical example

Consider the Earth-centered inertial state

(84)

(85)

with

(86)

First compute the specific angular momentum,

(87)

which gives approximately

(88)

with

(89)

The position magnitude is

(90)

so

(91)

Now evaluate

(92)

The result is

(93)

Its magnitude is

(94)

The orbit is therefore a low-eccentricity ellipse.

The specific mechanical energy is

(95)

which corresponds to

(96)

The semilatus rectum obtained dynamically is

(97)

while the orbital-element expression gives

(98)

The agreement provides a useful implementation check.

The angle between e and r is

(99)

which recovers the true anomaly used to generate the state.

Finally, the radial speed is

(100)

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

14 Implementation recipe from a Cartesian state

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

r	= ∥r∥,	(101)
h	= r × v,	(102)
e	= −,	(103)
e	= ∥e∥.	(104)

For e > 0, the unit periapsis direction is

(105)

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.

15 What to remember

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

(106)

Its final form is

(107)

Its magnitude is

(108)

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

(109)

For an ellipse,

(110)

and energy supplies the equivalent relation

(111)

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