Lecture 4: Gravitational Potential and N-Body Equations II.

Physics 141/241

Orbits in Spherical Potentials

Consider the motion of a star in a spherically-symmetric potential, MATH. The orbit of the star remains in a plane perpendicular to the angular momentum vector, and it's natural to adopt a polar coordinate system; call the coordinates MATH and $\varphi $. The system has $n=2$ degrees of freedom, so the phase space has $4$ dimensions.

The equations of motion can be derived by starting with the Lagrangian $L=K-V$,

(1) MATH,

where MATH and MATH. In what follows, we will use a choice of unit mass and $m$ will not show in the equations:

(2) MATH.

Differentiating with respect to MATH and MATH yields the momenta conjugate to $R$ and $\varphi $,

(3) MATH,

(4) MATH,

where $v_{R}$ and $v_{\varphi }$ are velocities in the radial and azimuthal directions. The Hamiltonian (energy) may now be expressed as a function of the coordinates and conjugate momenta:

(5) MATH.

The equations of motions are

(6) MATH,

(7) MATH,

(8) MATH,

(9) MATH.

Here $dJ/dt=0$ because the conjugate coordinate phi does not appear in $H$; $\varphi $ is called a cyclic coordinate.

The system has two integrals of motion. One, of course, is the total energy $E$, numerically equal to the value of $H$. The other is the angular momentum $J$. These quantities are given by

(10) MATH,

(11) MATH.

Each of these integrals of motion defines a hypersurface in phase space, and the orbit must remain in the intersection of these hypersurfaces. This can be visualized by ignoring the $\varphi $ coordinate and drawing surfaces of constant $E$ and J in the three-dimensional space MATH. Surfaces of constant $E$ are figures of revolution about the $R$ axis, while surfaces of constant $J$ are hyperbolas in the $(r,v_{\varphi })$ plane. The intersection of these surfaces is a closed curve, and an orbiting star travels around this curve.

For an orbit of a given $J$, the system may be reduced to one degree of freedom by defining the effective potential,

(12) MATH;

the corresponding equations of motion are then just

(13) MATH $v_{R}$,

(14) MATH $-d\Psi /dR$.

Because $\Psi (R)$ diverges as $R->0$, the star is energetically prohibited from coming too close to the origin, and shuttles back and forth between turning points $R_{\min }$ and $R_{\max }$.

In addition to its periodic radial motion described by $Eqs.(6-9)$, a star also executes a periodic azimuthal motion as it orbits the center of the potential. If the radial and azimuthal periods are incommensurate, as is usually the case, the resulting orbit never returns to its starting point in phase space; in coordinate space such an orbit is a rosette. The Keplerian potential is a very special case in which the radial and azimuthal periods of all bound orbits are equal. The only other potential in which all orbits are closed is the harmonic potential generated by a uniform sphere; here the radial period is half the azimuthal one and all bound orbits are ellipses centered on the bottom of the potential well. Thus in the Keplerian case all stars advance in azimuth by MATH between successive pericenters, while in the harmonic case they advance by MATH. Galaxies typically have mass distributions intermediate between these extreme cases, so most orbits in spherical galaxies are rosettes advancing by MATH between pericenters.

N-body Equations of Motion

Any system in which stellar collisions are rare may be idealized as a collection of $N$ point-sized bodies, each with mass $m_{i}$, position MATH, and velocity $\vec{v}_{i}$. The hamiltonian for such a system is

(15) MATH,

where $H$ depends on all body positions and velocities, the first sum runs over all $N$ bodies, the second runs over all pairs of bodies (twice, hence the factor of $1/2$), and $G$ is the gravitational constant. Then the equations of motion are

(16) MATH,

(17) MATH,

where the sum runs over all bodies except body $i$.

The famous Aarseth code in Fortran (Aarseth Fortran code) calculates the general motion of the N-body system in computer simulations.

N-body systems obey several basic conservation laws. A symmetry of the Hamiltonian is a transformation which leaves the physical system unchanged. For example, translation in time, $t->t+\Delta t$, is a symmetry of Eq. 15 because $H$ is not an explicit function of time; consequently the total system energy $E=K+U=H$ is conserved. Likewise, symmetry with respect to translation in space MATH, implies conservation of total linear momentum, and symmetry with respect to rotation gives rise to conservation of total angular momentum.

Virial Parameters

Another general result shown by manipulating Eqs. 16,17 is the scalar virial theorem which states that for a system in equilibrium,

(18) MATH

where $K$ and $U$ are the total kinetic and potential energy, respectively, and the angle-brackets indicate time-averages. Since $E=K+U$, the time-averaged kinetic and potential energies are related to the conserved total energy by

(19) MATH MATH .

The total mass $M$ and total energy $E$ of an N-body system thus define characteristic velocity and length scales

(20) MATH,

(21) MATH.

These are sometimes known as the virial velocity and radius, respectively.

The quantity $t_{c}=R/V$ is an estimate of the time a typical star takes to cross the system. This timescale may be expressed in several different ways; for example, in terms of the total mass $M$ and energy $E$, it is

(22) MATH.

Note that $M$ and $E$ are conserved, so $t_{c}$ is a constant even for systems which are far from dynamical equilibrium. In such cases $t_{c}$ approximates the time-scale over which the system evolves toward equilibrium.

Another expression for $t_{c}$ follows from the substitution $V^{2}=GM/R$ valid for systems near equilibrium:

(23) MATH.

Here the quantity $M/R^{3}$, which has units of density, appears. In systems with galaxy-like density profiles, the virial radius is approximately proportional to the half-mass radius: $R=2.5R_{h}$. Using this relationship, it follows that

(23) MATH,

where $\rho _{h}$ is the mean density within $R_{h}$. Since the crossing time is just supposed to indicate a typical time-scale for orbital motion, it is usual to drop the numerical constant, and define

(24) MATH.