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]
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)
\sim(10-4{pc
\sim(10-1{pc
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)
\sim(100kpc/100km/s,1011M\odot/1011)
\sim(10Mpc/1000km/s,1014M\odot/1077=m\nu)
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,Here in the N-body system, any individual member,
mi
mj
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.
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
g
\sigma\varphi\equiv
2 | |
\sqrt{\sigma | |
\varphi\varphi |
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
Twice kinetic energy per unit mass of the above uniform sphere is
M\proptor3\proptox3
The average Virial per unit mass can be computed from averaging its local value
r ⋅ (-\nabla\Phi)
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
The radial pressure gradient
The reason for the discrepancy is partly due to centrifugal force
0.2643V0=\sigma\varphi<\sigmar=0.5V0
r=0.8488r0
0.2643V0=\sigma\varphi>\sigmar=0
Now we can verify that
\beta\ne0
\langleV\varphi\rangle\ne0
\partial\varphi\Phi(x,t)=0
\partial\thetan(x,t)=0
{\partial\langleV\varphi\rangle\over\partialt}
-{\langleV\varphi\rangle\overtfric}
{\partialVr\over\partialt}=0
\nablaxmn(x,t)=0
{\partialn(x,t)\over\partialt}=0
\langleV(x,t)\rangle=0
\langleV(x)\rangle<0
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
So the fluid temperature per unit mass, i.e., the 1-dimensional velocity dispersion squared would be
Along the rotational z-axis,
z=\pmz0
z=0
P(0,0)={M0\over4\pi
2 | |
R | |
0 |
z0}{-GM0log[1-(1+R0/z
-2 | |
0) |
]\over2z0}.
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
tfric → infty
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},
~j=~r,~\theta,~\varphi
E=
In summary, in the spherical Jeans eq.,
\overline{r\partialr\Phi}=
2} | |
\overline{v | |
cir |
=\overline{GM\overr}=\overline{\langle
2 | |
v | |
t |
\rangle}
\overline{globalaverage