Stellar dynamics explained

Stellar dynamics is the branch of astrophysics which describes in a statistical way the collective motions of stars subject to their mutual gravity. The essential difference from celestial mechanics is that the number of body

N\gg10.

Typical galaxies have upwards of millions of macroscopic gravitating bodies and countless number of neutrinos and perhaps other dark microscopic bodies. Also each star contributes more or less equally to the total gravitational field, whereas in celestial mechanics the pull of a massive body dominates any satellite orbits.[1]

Connection with fluid dynamics

Stellar dynamics also has connections to the field of plasma physics.[2] The two fields underwent significant development during a similar time period in the early 20th century, and both borrow mathematical formalism originally developed in the field of fluid mechanics.

In accretion disks and stellar surfaces, the dense plasma or gas particles collide very frequently, and collisions result in equipartition and perhaps viscosity under magnetic field. We see various sizes for accretion disks and stellar atmosphere, both made of enormous number of microscopic particle mass,

(L/V,M/N)

\sim(10-8

55
pc/500km/s,1M
\odot/10

=mp)

at stellar surfaces,

\sim(10-4{pc

}/10,0.1M_/10^\sim m_) around Sun-like stars or km-sized stellar black holes,

\sim(10-1{pc

}/100,10M_/10^\sim m_) around million solar mass black holes (about AU-sized) in centres of galaxies.

The system crossing time scale is long in stellar dynamics, where it is handy to note that

1000pc/1km/s=1000Myr=HubbleTime/14.

The long timescale means that, unlike gas particles in accretion disks, stars in galaxy disks very rarely see a collision in their stellar lifetime. However, galaxies collide occasionally in galaxy clusters, and stars have close encounters occasionally in star clusters.

As a rule of thumb, the typical scales concerned (see the Upper Portion of P.C.Budassi's Logarithmic Map of the Universe) are

(L/V,M/N)

\sim(10pc/10km/s,1000M\odot/1000)

for M13 Star Cluster,

\sim(100kpc/100km/s,1011M\odot/1011)

for M31 Disk Galaxy,

\sim(10Mpc/1000km/s,1014M\odot/1077=m\nu)

for neutrinos in the Bullet Clusters, which is a merging system of N = 1000 galaxies.

Connection with Kepler problem and 3-body problem

At a superficial level, all of stellar dynamics might be formulated as an N-body problemby Newton's second law, where the equation of motion (EOM) for internal interactions of an isolated stellar system of N members can be written down as,m_i\frac = \sum_^N \frac. Here in the N-body system, any individual member,

mi

is influenced by the gravitational potentials of the remaining

mj

members.

In practice, except for in the highest performance computer simulations, it is not feasible to calculate rigorously the future of a large N system this way. Also this EOM gives very little intuition. Historically, the methods utilised in stellar dynamics originated from the fields of both classical mechanics and statistical mechanics. In essence, the fundamental problem of stellar dynamics is the N-body problem, where the N members refer to the members of a given stellar system. Given the large number of objects in a stellar system, stellar dynamics can address both the global, statistical properties of many orbits as well as the specific data on the positions and velocities of individual orbits.

Concept of a gravitational potential field

Stellar dynamics involves determining the gravitational potential of a substantial number of stars. The stars can be modeled as point masses whose orbits are determined by the combined interactions with each other. Typically, these point masses represent stars in a variety of clusters or galaxies, such as a Galaxy cluster, or a Globular cluster. Without getting a system's gravitational potential by adding all of the point-mass potentials in the system at every second, stellar dynamicists develop potential models that can accurately model the system while remaining computationally inexpensive.[3] The gravitational potential,

\Phi

, of a system is related to the acceleration and the gravitational field,

g

by:\frac }=\mathbf =-\nabla _\Phi (\mathbf),~~\Phi (\mathbf _)=-\sum _^.

Likewise the rms velocity in the rotation direction is computed by a weighted mean as follows, E.g.,\begin\langle\mathbf_\varphi^2\rangle(|\mathbf|)&\equiv \\& = \\& = 0.25V_0^2 = 0.5 \langle V_t^2 \rangle \\& = \\& =\langle\mathbf_\theta^2\rangle(|\mathbf|), \\\end

Here

\langle

2
V
t

\rangle=\langle

2
V
\theta

+

2\rangle
V
\varphi
2.
=0.5V
0

Likewise \langle\mathbf_r^2\rangle(\mathbf) = = \left(\sqrt \right)^2.

So the pressure tensor or dispersion tensor is \begin\sigma^2_(\mathbf)= & \\ =&\langle\mathbf_i\mathbf_j\rangle- \langle\mathbf_i\rangle\langle\mathbf_j\rangle \\= & \begin \left[1-({r \over r_0})^2\right]\left(\right)^2 & 0 & 0 \\ 0 & \left(\right)^2 & 0 \\ 0 & 0 & \left[1- ({8 \over 3 \pi})^2\right]\left(\right)^2 \end \endwith zero off-diagonal terms because of the symmetric velocity distribution. Note while there is no Dark Matter in producing the previous flat rotation curve, the price is shown by the reduction factor

{8\over3\pi}=0.8488

in the random velocity spread in the azimuthal direction. Among the diagonal dispersion tensor moments,

\sigma\theta\equiv

2
\sqrt{\sigma
\theta\theta
} = 0.5V_0 is the biggest among the three at all radii, and

\sigma\varphi\equiv

2
\sqrt{\sigma
\varphi\varphi
} \ge \sigma_r \equiv \sqrt only near the edge between

0.8488r0\ler\ler0

.

The larger tangential kinetic energy than that of radial motion seen in the diagonal dispersions is often phrased by an anisotropy parameter \beta(r) \equiv 1 - = 1 - = - \le 0; a positive anisotropy would have meant that radial motion dominated, and a negative anisotropy means that tangential motion dominates (as in this uniform sphere).

A worked example of Virial Theorem

Twice kinetic energy per unit mass of the above uniform sphere is

\begin & = \overline \equiv \langle \overline \rangle \\& = M_0^ \int_0^ \langle V_\theta^2+V_\varphi^2+V_r^2 \rangle dM \\& = M_0^ \int_0^1 \left(++\right) d(x^3 M_0) = 0.6 V_0^2, ~~ x \equiv = \left(\right)^,\endwhich balances the potential energy per unit mass of the uniform sphere, inside which

M\proptor3\proptox3

.

The average Virial per unit mass can be computed from averaging its local value

r(-\nabla\Phi)

, which yields\begin & = \overline\\&=M_0^ \int_0^ \mathbf \cdot (\rho d\mathbf^3) = -M_0^ \int_0^ dM \\& = -M_0^ \int_^ = - = -0.6 V_0^2,\end as required by the Virial Theorem. For this self-gravitating sphere, we can also verify that the virial per unit mass equals the averages of half of the potential \begin & = \overline \\& = M_0^ \int_^ d (M_0 x^3) \\& = = .\endHence we have verified the validity of Virial Theorem for a uniform sphere under self-gravity, i.e., the gravity due to the mass density of the stars is also the gravity that stars move in self-consistently; no additional dark matter halo contributes to its potential, for example.

A worked example of Jeans Equation in a uniform sphere

Jeans Equation is a relation on how the pressure gradient of a system should be balancing the potential gradient for an equilibrium galaxy. In our uniform sphere, the potential gradient or gravity is \nabla \Phi = = \ge 0, ~~\Omega = .

The radial pressure gradient - = - - = + 0 \ge 0.

The reason for the discrepancy is partly due to centrifugal force = > 0, and partly due to anisotropic pressure \begin &=0.25 \Omega^2 r \ge 0\\ & = 0.25\Omega^2 r - = \pm, \endso

0.2643V0=\sigma\varphi<\sigmar=0.5V0

at the very centre, but the two balance at radius

r=0.8488r0

, and then reverse to

0.2643V0=\sigma\varphi>\sigmar=0

at the very edge.

Now we can verify that\begin & = (-\sum_ V_i \partial_i \langle V_r\rangle) - \cancel - \nabla_r \Phi + \sum_\\&= - 0 - + \left[-{d (\rho \sigma_r^2) \over \rho dr} + {\sigma_\theta^2 + \sigma_\varphi^2 -2\sigma_r^2 \over r} \right] \\& = -(\Omega^2 r) + \\& \left[{\Omega^2 r \over 2} + {(0.5V_0)^2 + (0.2643V_0)^2- 2 \times 0.25 \Omega^2 (r_0^2-r^2) \over r} \right] \\& = 0.\end Here the 1st line above is essentially the Jeans equation in the r-direction, which reduces to the 2nd line, the Jeans equation in an anisotropic (aka

\beta\ne0

) rotational (aka

\langleV\varphi\rangle\ne0

) axisymmetric (

\partial\varphi\Phi(x,t)=0

) sphere (aka

\partial\thetan(x,t)=0

) after much coordinate manipulations of the dispersion tensor; similar equation of motion can be obtained for the two tangential direction, e.g.,

{\partial\langleV\varphi\rangle\over\partialt}

, which are useful in modelling ocean currents on the rotating earth surface or angular momentum transfer in accretion disks, where the frictional term

-{\langleV\varphi\rangle\overtfric}

is important. The fact that the l.h.s.

{\partialVr\over\partialt}=0

means that the force is balanced on the r.h.s. for this uniform (aka

\nablaxmn(x,t)=0

) spherical model of a galaxy (cluster) to stay in a steady state (aka time-independent equilibrium

{\partialn(x,t)\over\partialt}=0

everywhere) statically (aka with zero flow

\langleV(x,t)\rangle=0

everywhere). Note systems like accretion disk can have a steady net radial inflow

\langleV(x)\rangle<0

everywhere at all time.

A worked example of Jeans equation in a thick disk

Consider again the thick disk potential in the above example. If the density is that of a gas fluid, then the pressure would be zero at the boundary

z=\pmz0

. To find the peak of the pressure, we note that P(R,z) = \int^_z \partial_z\Phi \rho(R) dz = \rho(R) [\Phi(R,z_0) - \Phi(R,z)].

So the fluid temperature per unit mass, i.e., the 1-dimensional velocity dispersion squared would be \sigma^2(R,z) =, ~~ |z| \le z_0

\sigma^2= \log, ~~ Q(z) \equiv R_0 + z_0 + z +\sqrt.

Along the rotational z-axis, \sigma^2(0,z) = \log \sigma(0,z) = \sqrt \sqrt, which is clearly the highest at the centre and zero at the boundaries

z=\pmz0

. Both the pressure and the dispersion peak at the midplane

z=0

. In fact the hottest and densest point is the centre, where

P(0,0)={M0\over4\pi

2
R
0

z0}{-GM0log[1-(1+R0/z

-2
0)

]\over2z0}.

A recap on worked examples on Jeans Eq., Virial and Phase space density

Having looking at the a few applications of Poisson Eq. and Phase space density and especially the Jeans equation, we can extract a general theme, again using the Spherical cow approach.

Jeans equation links gravity with pressure gradient, it is a generalisation of the Eq. of Motion for single particles. While Jeans equation can be solved in disk systems, the most user-friendly version of the Jeans eq. is the spherical anisotropic version for a static

\langle{vj}\rangle=0

frictionless system

tfricinfty

, hence the local velocity speed
2
\sigma
j

(r)=

2}\rangle(r)
\langle{v
j

-

2(r)}
\underbrace{\langle{v
=0

={\int\limitsinftydvrdv\thetadv\varphi

=0
({v}
j}

)2fp\over\int\limitsinftydvrdv\thetadv\varphifp},

everywhere for each of the three directions

~j=~r,~\theta,~\varphi

.One can project the phase space into these moments, which is easily if in a highly spherical system, which admits conservations of energy

E=

and angular momentum J. The boundary of the system sets the integration range of the velocity bound in the system.

In summary, in the spherical Jeans eq., \begin = & \\= & -+, \\= & -, ~~\text \\= &, ~~\text, \langle v_t^2 \rangle \equiv \langle\rangle + \langle\rangle \end which matches the expectation from the Virial theorem

\overline{r\partialr\Phi}=

2}
\overline{v
cir

=\overline{GM\overr}=\overline{\langle

2
v
t

\rangle}

, or in other words, the

\overline{globalaverage

} kinetic energy of an equilibrium equals the average kinetic energy on circular orbits with purely transverse motion.

See also

Further reading

Notes and References

  1. Book: Murdin, Paul . Encyclopedia of Astronomy and Astrophysics . Nature Publishing Group. 2001 . 978-0750304405 . 1. Stellar Dynamics.
  2. https://cds.cern.ch/record/1053485/files/p37.pdf
  3. Book: Galactic Dynamics. Binney. James. Tremaine. Scott . Princeton University Press . 2008 . 978-0-691-13027-9 . Princeton. 35, 63, 65, 698.