1. Introduction
The search for self-gravitating fluid instabilities is an interesting subject for the determination of the structure formation of interstellar gas clouds. The first analysis of the self-gravitating fluid instabilities from the hydrodynamic equations coupled with the Newtonian Poisson equation was performed by Jeans [
1] in 1902. In this work, he obtained a wavelength cutoff (Jeans wavelength) where perturbations with wavelengths smaller than the Jeans wavelength propagate as harmonic waves in time, while perturbations with wavelengths larger than the Jeans wavelength grow or decay in time. The gravitational collapse of self-gravitating interstellar gas clouds, which is connected with the exponential growth in time of the mass density perturbations, is known as the Jeans instability [
2,
3,
4]. From a physical point of view, we may say that a collapse of a mass density inhomogeneity happens whenever the outwards pressure force of the self-gravitating gas cloud is smaller than the inwards gravitational force.
There are two main methods to treat the problem of self-gravitating fluid instabilities, one by using the hydrodynamic equations coupled with the Newtonian Poisson equation and another to consider the collisionless Boltzmann equation coupled with the Newtonian Poisson equation (see, for example, [
3,
4]).
The self-gravitating fluid instabilities were also considered in the
gravity theory [
5,
6,
7], where a modified dispersion relation implies new unstable modes.
Recently, the self-gravitating fluid instabilities were examined within the framework of the hydrodynamic equations and Poisson equations, which follow from the post-Newtonian theory (see [
8,
9]) in [
10,
11,
12]. The same problem was analysed within the framework of a post-Newtonian collisonless Boltzmann equation (see [
13,
14,
15]) in [
16,
17]. In these works, it was shown that the mass which is necessary for an overdensity to begin the gravitational collapse in the post-Newtonian theory is smaller than the one given by the Newtonian theory.
In the present work, we reexamine the self-gravitating fluid instabilities from the post-Newtonian collisionless Boltzmann equation coupled with the post-Newtonian Poisson equations in order to determine the difference between the Newtonian and post-Newtonian theories in the growth of the instabilities dictated by the mass density contrast. It is shown that the growth of the mass density contrast in the post-Newtonian theory is more accentuated than that of the Newtonian theory. The methodology used here differs from the one given in [
16] since the amplitudes of the gravitational potentials and one-particle distribution function perturbations are considered to be time-dependent. Furthermore, the time-dependent amplitude of the one-particle distribution function is considered a linear combination of the summational invariants from the post-Newtonian kinetic theory.
The structure of the paper is the following: in
Section 2, we introduce the post-Newtonian expressions for the Boltzmann and Poisson equations and the energy–momentum tensor defined in terms of the one-particle distribution function. The perturbations of the one-particle distribution function and gravitational potentials from a stationary equilibrium background is the subject of
Section 3. In
Section 4, the perturbations are represented by plane waves with small time-dependent amplitudes, and from the system of differential equations for the amplitudes, a time-evolution differential equation for the mass density contrast is found. The analysis of the solutions of the time-evolution differential equation for the mass density contrast is the topic of
Section 5. The summary and conclusions of the work are given in
Section 6 and
Appendix A close the work.
2. Post-Newtonian Boltzmann Equation
The statistical description of a self-gravitating gas can be described by the Boltzmann equation, which governs the space-time evolution of the one-particle distribution function
defined in the phase space spanned by the spatial coordinates
and velocity
of the particles. Here, we are interested in analysing the collisionless post-Newtonian Boltzmann equation (see [
9] and the reference therein), which can be written in the first post-Newtonian approximation as [
16]
Above
is the three-velocity component of the particle four-velocity
in first post-Newtonian approximation [
8]
Furthermore, in (
1),
U and
are scalar gravitational potentials, while
is a vector gravitational potential. They are defined in the first post-Newtonian approximation of the metric tensor
components by [
8]
The gravitational potentials
U,
and
satisfy Poisson equations, which are given in terms of the energy–momentum tensor split in orders of
denoted by
, namely
The energy–momentum tensor is defined in kinetic theory of gases by the one-particle distribution function [
18]
where the first post-Newtonian approximation of the invariant integration element is given by [
19]
3. Perturbations of the One-Particle Distribution Function and Gravitational Potentials
We shall write the one-particle distribution function and the gravitational potentials as sums of background and perturbed terms, where the background terms correspond to an equilibrium state.
For a relativistic gas, the equilibrium state of the one-particle distribution function is given by the Maxwell–Jüttner distribution function [
18]. For a stationary equilibrium background, where the hydrodynamic velocity vanishes the first post-Newtonian approximation of the Maxwell–Jüttner distribution function
is [
19]
In the above equation, is the non-relativistic Maxwellian distribution function, which is a function of the particle rest mass m, mass density , gas particle velocity , and the dispersion velocity .
By denoting the background terms of the gravitational potentials by the subscript zero and the perturbed terms of the one-particle distribution function and gravitational potentials by the subscript 1, we write
From the substitution of the representations (
9)–(
12) into the Boltzmann Equation (
1), we obtain an equation for the space-time evolution of the perturbed one-particle distribution function [
16]
Above, the background Maxwell–Jüttner distribution function is denoted by
. Here, we note that the background equation derived from (
1) is identically satisfied when
,
and
.
The conditions of vanishing background potential gravitational gradients do not satisfy the Poisson equations, and as usual, we have to assume “Jeans swindle” (see [
4]), where the Poisson Equations (
4) and (
5) are valid only for the perturbed distribution function and gravitational potentials. Hence, from (
4) and (
5) together with the (
6) and the representations (
9)–(
12) yield
4. Representation of Perturbations as Plane Waves
To go further in the analysis of the instabilities, we represent the perturbations as plane waves of wave number vector
and time-dependent small amplitudes
,
,
and
, namely
From the insertion of the plane wave representations (
17) and (
18) into the perturbed Boltzmann Equation (
13), we obtain
while the Poisson Equations (
14)–(
16) with the plane wave representations (
17) and (
18) become
Now, we follow the methodology developed in [
9] and consider the perturbed one-particle distribution function
proportional to a sum of summational invariants, which are quantities that are conserved in the binary collision of the particles. In the relativistic kinetic theory of gases, the summational invariants are the rest mass of a particle
m and the momentum four-vector
. Here, we write the perturbed one-particle distribution function as a linear combination of the summational invariants
, where
and
are unknowns that do not depend on the momentum four-vector
.
The post-Newtonian approximation of the linear combination of the summational invariants
obtained from (
2) and (
3) and by considering terms up to the order
reads
where new unknowns
,
and
were introduced defined by
From (
23), we can identify the summational invariants in the post-Newtonian approximation, namely
By neglecting the terms of order , the above summational invariants reduce to those of the non-relativistic theory: 1, and .
The representation of the perturbed one-particle distribution function follows from the product of the background Maxwell–Jüttner distribution function
and the linear combination of the summational invariants (
23):
The dependence of the perturbed gravitational potentials
,
and
on the unknown amplitudes
,
and
of the perturbed one-particle distribution function follows from the insertion of (
27) into (
20)–(
22) and the integration of the resulting equations, yielding
In the above equations, we introduced the modulus of the wave number vector
and the Jeans wavelength
. Note that (
29) follows from the scalar product of the vector equation for
by
.
Now, we introduce the perturbed distribution function (
27) into the perturbed Boltzmann Equation (
19) and multiply the resulting equation by each of the summational invariants given in (
26). From the integration of the resulting equations, by taking into account the invariant element of integration (
7), the following system of differential equations emerges:
If we combine (
31) and (
32), we obtain the following differential equation:
Hence, we can express the unknown in terms of the unknown and without loss of generality we can choose the constant .
Here, we introduce the mass density contrast
, which refers to a parameter that indicates where a local increase in the matter density happens. The density contrast is defined by the ratio of the perturbed and unperturbed mass densities, namely,
An evolution equation for the mass density contrast can be obtained from the differentiation of (
31) with respect to time by using the relationships (
34) and (
35) and the elimination of
from (
33) and of the gravitational potentials
,
and
by considering (
28), (
29) and (
30), respectively. After some arrangements and by introducing the dimensionless time
, the wavelength of the perturbation
and the Jeans wavelength
, we obtain
5. The Time Evolution of the Mass Density Contrast
The solution of the differential Equation (
36) for the mass density contrast reads
where
and
are arbitrary constants.
Let us analyse the numerical solutions of the mass density contrast differential Equation (
36) for the cases where the wavelength of the perturbation is bigger than the Jeans wavelength
and smaller
.
The mass density contrast differential Equation (
36) was solved for the initial conditions
and
(say), and the numerical solutions are plotted in
Figure 1 and
Figure 2 for different values of the ratio between the dispersion velocity and the light speed, namely,
, which corresponds to the Newtonian theory and
and
to the post-Newtonian theory. Furthermore, in the evaluation of (
36), it was supposed that
, which can be justified by the virial theorem, where the Newtonian gravitational potential can be approximated with the square of the dispersion velocity.
In
Figure 1, the time evolution of the mass density contrast
is shown for the case where the wavelength of the perturbation is bigger than the Jeans wavelength (
) and by considering different values of the ratio between the dispersion velocity and the light speed
. We observe from this figure that in this case, the mass density contrast grows in dimensionless time
. We may also infer from this figure that the growth of the instability in the post-Newtonian theory differs from the one of the Newtonian theory, the former growing more rapidly than the latter. Furthermore, the increase in the ratio
implies a more accentuated growth of the mass density contrast with the dimensionless time. As was pointed out previously, the so-called Jeans instability is connected with the gravitational collapse of self-gravitating interstellar gas clouds, which are associated with the mass density perturbations, which grow exponentially in time [
2,
3,
4].
The time evolution of the mass density is plotted in
Figure 2 for the case where the perturbation wavelength is smaller than the Jeans wavelength (
), where different values of the ratio between the dispersion velocity and the light speed
were considered. We conclude from this figure that the mass density contrast propagates as harmonic waves in time. Note that there is no remarkable difference between the Newtonian and post-Newtonian solutions for the values adopted here.
6. Conclusions
In this work, the self-gravitating fluid instabilities were analysed within the framework of the post-Newtonian Boltzmann equation. The post-Newtonian Boltzmann equation is coupled with three Poison equations for the post-Newtonian gravitational potentials, which are functions of the energy–momentum tensor and determined from the knowledge of the post-Newtonian Maxwell–Jüttner distribution function. The one-particle distribution function and the gravitational potentials were perturbed from their background states represented by plane waves characterised by a wave number vector and time-dependent small amplitudes. The time-dependent amplitude of the one-particle distribution function was supposed to be a linear combination of the summational invariants of the post-Newtonian kinetic theory. From the coupled system of differential equations for the time-dependent amplitudes of the one-particle distribution function and gravitational potentials, an evolution equation for the mass density contrast was obtained. It was shown that for perturbation wavelengths smaller than the Jeans wavelength, the mass density contrast propagates as harmonic waves in time. For perturbation wavelengths greater than the Jeans wavelength, the mass density contrast grows in time, and the instability growth in the post-Newtonian theory is more accentuated than that of the Newtonian theory.
This work improves the former works on Jeans instability, where it was shown that in the post-Newtonian theory, the mass necessary for an overdensity to begin the gravitational collapse is smaller than in the Newtonian theory. Here, we have shown that in the post-Newtonian theory, the instabilities grow more rapidly than those in the Newtonian theory.