In the general case, the CBE has non less than 7 independent variables and cannot be solved. One way to reduce the complexity of the problem is to consider moments of the CBE (e.g. the continuity equation, the Jeans equations, the Virial Theorem), which are defined in the 3D space and can be connected to observations. They allow us to evaluate the general properties of the systems under exam and the properties of the CBE solutions, but they do not provide the full physics of the stellar systems, which is stored in the exact CBE solutions. Exact CBE solutions can be obtained for simple systems.
In this lecture we are interested in deriving some of these solutions as an exercise to understand how observed quantities can be a reflection of the internal dynamics of equilibrium systems. We start by deriving an important theorem of the galaxy dynamics: the Jeans Theorem.
First of all we recall that a function of the phase-space coordinates : is an integral of motion if and only if which means that:
Integral of motions are important to define orbits. In fact, if at any given time we know 5 integral of motion in the 6-D space, each of them defines a 5-D hyper-surface and their intersection defines a 1-D trajectory, the orbit.
Jeans Theorem: Any steady-state solution of the CBE depends on the phase-space coordinates only through integrals of motion in the galactic potential, and any function of the integrals yields a steady-state solution of the CBE. Suppose is a steady-state solution of the CBE:
otherwise written as:
which is means that is an integral of motion.
Conversely, if are integrals of motion and is any function of variables then:
because by definition of integral of motion:
Thus is a solution of the CBE.
We will be mainly interested in the second proposition of the Jeans Theorem in the sense that any distribution function, which is written in terms of integral of motion is assured to be a solution of the CBE. There is a strong version of the Jeans Theorem (that we do not prove here), which states that the distribution function of steady-state galaxy in which almost all orbits are regular with non-resonant frequencies may be presumed to be a function only of three independent isolating integrals.
Spherical systems admit 4 integrals of motion: the energy and angular momentum. A function of these four integrals is a solution of the CBE. If the system is also spherically symmetric, then depends only on the magnitude of ; that is: We remark that, if the potential is provided by the tracer’s mass, then self-consistency holds, and:
or, using spherical symmetry,
For the following calculations we will define a relative potential: where is some constant value, and the relative energy of a star:
In particular is chosen in order to have for , and for .
Obviously, the relative potential of an isolated system satisfies the Poisson’s equation of the form:
with boundary conditions as
The equation of state of an isothermal body of gas is , where is a constant.
The equation of hydrostatic support of an isothermal gas is:
where is the Boltzmann constant, and are the pressure and temperature of the gas, is the mass of the particle, and is the mass of the system within By multiplying the equation above by and differentiating, we obtain:
where we have used
Let us assume now that the stellar system has a distribution function of the form:
Integrating over all velocities, we derive the number density of the system:
thus, the Poisson’s Equation becomes:
By comparing equation (1) with eq. (4) we see that, if we define: the two equations are identical. This means that the structure of an isothermal self-gravitating sphere of gas is identical with the structure of collisionless system of stars whose density in phase space is given by equation (2).
In particular we can find out which number density distribution is given for a isothermal sphere by solving equation (4). If we assume the solutions of the form of a power law , then the left side of the equation:
and the right side:
They are equal if: and which implies:
This is the density of the singular isothermal sphere because the density is infinite at It is possible to derive other solutions, which do not have this kind of problem (and are more physical). A well know solution is the so called King’s law which has the form: where is the central density and: is the King radius.
In equation (4) we define: and such that:
which, using Eq. (3), can be re-written as:
The new density law, which can be derived by numerical integration of equation (6), is characterized by the presence of a “core” in the central regions where the density flattens to constant value At large distance () it goes as the singular isothermal sphere (eq. 5), which in the new coordinates writes:
It can be shown that a good approximation for the solution of equation (6) in the central regions is the modified Hubble law:
The surface density of the isothermal sphere in units of for is:
We remark that, if
is the mass interior to the spherical radius then the circular speed is written as: which, by equation(4), can be also re-written as:
where the term tends to -2 at for the isothermal sphere. Thus the circular speed at large distances from the center is constant to
We finally show how one can derive the distribution function of a spherically symmetric star system starting from the density profile By definition we know that the density profile is the integral of the DF over all velocities. If we write as a function of the energy and we replace the integral over all velocities as the integral over using a change of variable, we obtain:
As far as is a monotonic function of we can consider as a function of and obtain:
This is an Abel integral equation of the kind already found previously in this course, which can be solved to obtain :
which is also called Eddington’s formula from its original inventor, sir Arthur Eddington (1916).
This formula states that one can derive the DF (as a function of the energy only) of a spherical (non rotating) systems just by knowing its density distribution. However, there is no guarantee that the solution of this equation will satisfy the physical requirement of being nowhere negative.
One can show that this is true if and only if:
is an increasing function of
8. Spiral arms
10. Scale relations
14. Galaxy dynamics
19. Galaxy clusters