Issue 
A&A
Volume 556, August 2013



Article Number  A93  
Number of page(s)  27  
Section  Astrophysical processes  
DOI  https://doi.org/10.1051/00046361/201220607  
Published online  05 August 2013 
Kinetic theory of spatially inhomogeneous stellar systems without collective effects^{⋆}
Laboratoire de Physique Théorique (IRSAMC), CNRS and UPS, Université de
Toulouse,
31062
Toulouse Cedex 4,
France
email:
chavanis@irsamc.upstlse.fr
Received:
21
October
2012
Accepted:
27
March
2013
We review and complete the kinetic theory of spatially inhomogeneous stellar systems when collective effects (dressing of the stars by their polarization cloud) are neglected. We start from the BBGKY hierarchy issued from the Liouville equation and consider an expansion in powers of 1/N in a proper thermodynamic limit. For N → +∞, we obtain the Vlasov equation describing the evolution of collisionless stellar systems like elliptical galaxies. This corresponds to the mean field approximation. At the order 1/N, we obtain a kinetic equation describing the evolution of collisional stellar systems like globular clusters. This corresponds to the weak coupling approximation. This equation coincides with the generalized Landau equation derived from a more abstract projection operator formalism. This equation does not suffer logarithmic divergences at large scales since spatial inhomogeneity is explicitly taken into account. Making a local approximation, and introducing an upper cutoff at the Jeans length, it reduces to the VlasovLandau equation which is the standard kinetic equation of stellar systems. Our approach provides a simple and pedagogical derivation of these important equations from the BBGKY hierarchy which is more rigorous for systems with longrange interactions than the twobody encounters theory. Making an adiabatic approximation, we write the generalized Landau equation in angleaction variables and obtain a Landautype kinetic equation that is valid for fully inhomogeneous stellar systems and is free of divergences at large scales. This equation is less general than the recently derived LenardBalescutype kinetic equation since it neglects collective effects, but it is substantially simpler and could be useful as a first step. We discuss the evolution of the system as a whole and the relaxation of a test star in a bath of field stars. We derive the corresponding FokkerPlanck equation in angleaction variables and provide expressions for the diffusion coefficient and friction force.
Key words: gravitation / methods: analytical / globular clusters: general
Appendices are available in electronic form at http://www.aanda.org
© ESO, 2013
1. Introduction
In its simplest description, a stellar system can be viewed as a collection of N classical point mass stars in Newtonian gravitational interaction (Spitzer 1987; Heggie & Hut 2003; Binney & Tremaine 2008). As understood early on by Hénon (1964), selfgravitating systems experience two successive types of relaxation, a rapid collisionless relaxation towards a quasistationary state (QSS) that is a virialized state in mechanical equilibrium but not in thermodynamical equilibrium, followed by a slow collisional relaxation. One might think that, because of the development of stellar encounters, the system will reach, at sufficiently long times, a statistical equilibrium state described by the MaxwellBoltzmann distribution. However, it is wellknown that unbounded stellar systems cannot be in strict statistical equilibrium^{1} because of the permanent escape of high energy stars (evaporation) and the gravothermal catastrophe (core collapse). Therefore, the statistical mechanics of stellar systems is essentially an outofequilibrium problem which must be approached through kinetic theories.
The first kinetic equation was written by Jeans (1915). Neglecting encounters between stars, he described the dynamical evolution of stellar systems by the collisionless Boltzmann equation coupled to the Poisson equation. This purely mean field description applies to large groups of stars such as elliptical galaxies whose ages are much less than the collisional relaxation time. A similar equation was introduced by Vlasov (1938) in plasma physics to describe the collisionless evolution of a system of electric charges interacting by the Coulomb force. The collisionless Boltzmann equation coupled selfconsistently to the Poisson equation is often called the Vlasov equation^{2}, or the VlasovPoisson system.
The concept of collisionless relaxation was first understood by Hénon (1964) and King (1966). LyndenBell (1967) developed a statistical theory of this process and coined the term “violent relaxation”. In the collisionless regime, the evolution of the star cluster is described by the VlasovPoisson system. Starting from an unsteady or unstable initial condition, the VlasovPoisson system develops a complicated mixing process in phase space. Because the Vlasov equation is timereversible, it never achieves a steady state but develops filaments at smaller and smaller scales. However, the coarsegrained distribution function obtained by locally averaging the finegrained distribution function over the filaments usually achieves a steady state on a few dynamical times. LyndenBell (1967) tried to predict the QSS resulting from violent relaxation by developing a statistical mechanics of the Vlasov equation. He derived a distribution function formally equivalent to the FermiDirac distribution (or to a superposition of FermiDirac distributions). However, when coupled to the Poisson equation, these distributions have an infinite mass. Therefore, the VlasovPoisson system has no statistical equilibrium state (in the sense of LyndenBell). This is clear evidence that violent relaxation is incomplete (LyndenBell 1967). Incomplete relaxation is caused by inefficient mixing and nonergodicity. In general, the fluctuations of the gravitational potential δΦ(r, t) that drive the collisionless relaxation last only for a few dynamical times and die out before the system has mixed efficiently (Tremaine et al. 1986). Understanding the origin of incomplete relaxation, and developing models of incomplete violent relaxation to predict the structure of galaxies, is a very difficult problem (Arad & Johansson 2005). Some models of incomplete violent relaxation have been proposed based on different physical arguments (Bertin & Stiavelli 1984; Stiavelli & Bertin 1987; Hjorth & Madsen 1991; Chavanis et al. 1996; Chavanis 1998; Levin et al. 2008).
On longer timescales, stellar encounters (sometimes referred to as collisions by an abuse of language) must be taken into account. This description is particularly important for small groups of stars such as globular clusters whose ages are of the order of the collisional relaxation time. Chandrasekhar (1942, 1943a,b) developed a kinetic theory of stellar systems to determine the timescale of collisional relaxation and the rate of escape of stars from globular clusters^{3}. To simplify the kinetic theory, he considered an infinite homogeneous system of stars. He started from the general FokkerPlanck equation and determined the diffusion coefficient and the friction force (first and second moments of the velocity increments) by considering the mean effect of a succession of twobody encounters^{4}. This approach was based on Jeans (1929) demonstration that the cumulative effect of the weak deflections resulting from the relatively distant encounters is more important than the effect of occasional large deflections produced by relatively close encounters. However, this approach leads to a logarithmic divergence at large scales that is more difficult to remove in stellar dynamics than in plasma physics because of the absence of Debye shielding for the gravitational force^{5}. Chandrasekhar & von Neumann (1942) developed a completely stochastic formalism of gravitational fluctuations and showed that the fluctuations of the gravitational force are given by the Holtzmark distribution (a particular Lévy law) in which the nearest neighbor plays a prevalent role. From these results, they argued that the logarithmic divergence has to be cutoff at the interparticle distance l (see also Jeans 1929; and Spitzer 1940). However, since the interparticle distance is smaller than the Debye length, the same arguments should also apply in plasma physics, which is not the case. Therefore, the conclusions of Chandrasekhar and von Neumann are usually taken with circumspection. In particular, Cohen et al. (1950) argued that the logarithmic divergence should be cutoff at the Jeans length which gives an estimate of the system’s size. While in neutral plasmas the effective interaction distance is limited to the Debye length, in a selfgravitating system the distance between interacting particles is only limited by the system’s size. Therefore, the Jeans length is the gravitational analogue of the Debye length. These kinetic theories lead to a collisional relaxation time scaling as t_{R} ~ (N/lnN)t_{D}, where t_{D} is the dynamical time and N the number of stars in the system. Chandrasekhar (1949) also developed a Brownian theory of stellar dynamics and showed that, from a qualitative point of view, the results of kinetic theory can be understood very simply in that framework^{6}. In particular, he showed that a dynamical friction is necessary to maintain the MaxwellBoltzmann distribution of statistical equilibrium and that the coefficients of friction and diffusion are related to each other by an Einstein relation which is a manifestation of the fluctuationdissipation theorem. This relation is confirmed by his more precise kinetic theory based on twobody encounters. It is important to emphasize, however, that Chandrasekhar did not derive the kinetic equation for the evolution of the system as a whole. Indeed, he considered the Brownian motion of a test star in a fixed distribution of field stars (bath) and derived the corresponding FokkerPlanck equation. This equation has been used by Chandrasekhar (1943b), Spitzer & Härm (1958), Michie (1963), King (1965), and more recently by Lemou & Chavanis (2010) to study the evaporation of stars from globular clusters in a simple setting.
King (1960) noted that if we were to describe the dynamical evolution of the cluster as a whole, the distribution of the field stars should evolve in time in a selfconsistent manner so that the kinetic equation must be an integrodifferential equation. The kinetic equation obtained by King, from the results of Rosenbluth et al. (1957), is equivalent to the Landau equation, although written in a different form. It is interesting to note, for historical reasons, that none of the previous authors seemed to be aware of the work of Landau (1936) in plasma physics. There is, however, an important difference between stellar dynamics and plasma physics. Neutral plasmas are spatially homogeneous because of electroneutrality and Debye shielding. By contrast, stellar systems are spatially inhomogeneous. The abovementioned kinetic theories developed for an infinite homogeneous system can be applied to an inhomogeneous system only if we make a local approximation. In that case, the collision term is calculated as if the system were spatially homogeneous or as if the collisions could be treated as local. Then, the effect of spatial inhomogeneity is only retained in the advection (Vlasov) term which describes the evolution of the system due to mean field effects. This leads to the VlasovLandau equation which is the standard kinetic equation of stellar dynamics. To our knowledge, this equation has been first written (in a different form), and studied, by Hénon (1961). Hénon also exploited the timescale separation between the dynamical time t_{D} and the relaxation time t_{R} ≫ t_{D} to derive a simplified kinetic equation for f(ϵ,t), where ϵ = v^{2}/2 + Φ(r, t) is the individual energy of a star by unit of mass, called the orbitaveraged FokkerPlanck equation. In this approach, the distribution function f(r, v, t), averaged over a short timescale, is a steady state of the Vlasov equation of the form f(ϵ,t) which slowly evolves in time, on a long timescale, due to the development of collisions (i.e. correlations caused by finite N effects or graininess). Hénon used this equation to obtain a more relevant value for the rate of evaporation from globular clusters, valid for inhomogeneous systems. Cohn (1980) solved the orbitaveraged FokkerPlanck equation numerically to describe the collisional evolution of star clusters. His treatment accounts both for the escape of high energy stars put forward by Spitzer (1940) and for the phenomenon of core collapse resulting from the gravothermal catastrophe discovered by Antonov (1962) and LyndenBell & Wood (1968) on the basis of thermodynamics and statistical mechanics. The local approximation, which is a crucial step in the kinetic theory, is supported by the stochastic approach of Chandrasekhar & von Neumann (1942) showing the preponderance of the nearest neighbor. However, this remains a simplifying assumption which is not easily controllable. In particular, as we have already indicated, the local approximation leads to a logarithmic divergence at large scales that is difficult to remove. This divergence would not have occurred if full account of spatial inhomogeneity had been given from the start.
The effect of spatial inhomogeneity was investigated by Severne & Haggerty (1976), Parisot & Severne (1979), Kandrup (1981), and Chavanis (2008a,b). In particular, Kandrup (1981) derived a generalized Landau equation from the Liouville equation by using projection operator technics. This generalized Landau equation is interesting because it takes into account effects of spatial inhomogeneity which were neglected in previous approaches. Since the finite extension of the system is properly accounted for, there is no divergence at large scales^{7}. Furthermore, this approach clearly shows which approximations are needed in order to recover the traditional Landau equation. Unfortunately, the generalized Landau equation remains extremely complicated for practical applications.
In addition, this equation is still approximate as it neglects collective effects and considers binary collisions between naked particles. As in any weakly coupled system, the particles engaged in collisions are dressed by the polarization clouds caused by their own influence on other particles. Collisions between dressed particles have quantitatively different outcomes than collisions between naked ones. In the case of plasmas, collective effects are responsible for Debye shielding and they are accounted for in the LenardBalescu equation. They eliminate the logarithmic divergence that occurs at large scales in the Landau equation. For selfgravitating systems, they lead to antishielding and are more difficult to analyze because the system is spatially inhomogeneous^{8}. If we consider a finite homogeneous system, and take collective effects into account, one finds a severe divergence of the diffusion coefficient when the size of the domain reaches the Jeans scale (Weinberg 1993). This divergence, which is related to the Jeans instability, does not occur in a stable spatially inhomogeneous stellar system. Some authors like Thorne (1968), Miller (1968), Gilbert (1968, 1970), and Lerche (1971) attempted to take collective effects and spatial inhomogeneity into account. They obtained very complicated kinetic equations that have not found application until now. They managed, however, to show that collective effects are equivalent to increasing the effective mass of the stars, hence diminishing the relaxation time. Since, on the other hand, the effect of spatial inhomogeneity is to increase the relaxation time (Parisot & Severne 1979), the two effects act in opposite directions and may balance each other.
Recently, Heyvaerts (2010) derived from the BBGKY hierarchy issued from the Liouville equation a kinetic equation in angleaction variables that takes both spatial inhomogeneity and collective effects into account. To calculate the collective response, he used FourierLaplace transforms and introduced a biorthogonal basis of pairs of densitypotential functions (Kalnajs 1971a). The kinetic equation derived by Heyvaerts is the counterpart for spatially inhomogeneous selfgravitating systems of the LenardBalescu equation for plasmas. Following his work, we showed that this equation could be obtained equivalently from the Klimontovich equation by making the socalled quasilinear approximation (Chavanis 2012a). We also developed a test particle approach and derived the corresponding FokkerPlanck equation in angleaction variables, taking collective effects into account. This provides general expressions of the diffusion coefficient and friction force for spatially inhomogeneous stellar systems.
In a sense, these equations solve the problem of the kinetic theory of stellar systems since they take into account both spatial inhomogeneity and collective effects. However, the drawback is that they are extremely complicated to solve (in addition of being complicated to derive). In an attempt to reduce the complexity of the problem, we shall derive in this paper a kinetic equation that is valid for spatially inhomogeneous stellar systems but that neglects collective effects. Collective effects may be less crucial in stellar dynamics than in plasma physics. In plasma physics, they must be taken into account in order to remove the divergence at large scales that appears in the Landau equation. In the case of stellar systems, this divergence is removed by the spatial inhomogeneity of the system, not by collective effects. Actually, previous kinetic equations based on the local approximation ignore collective effects and already give satisfactory results. We shall therefore ignore collective effects and derive a kinetic equation (in positionvelocity and angleaction variables) that is the counterpart for spatially inhomogeneous selfgravitating systems of the Landau equation for plasmas. Our approach has three main interests. First, the derivation of this Landautype kinetic equation is considerably simpler than the derivation of the LenardBalescutype kinetic equation, and it can be done in the physical space without having to introduce LaplaceFourier transforms nor biorthogonal basis of pairs of densitypotential functions. This offers a more physical derivation of kinetic equations of stellar systems that may be of interest for astrophysicists. Secondly, our approach is sufficient to remove the largescale divergence that occurs in kinetic theories based on the local approximation. It represents, therefore, a conceptual progress in the kinetic theory of stellar systems. Finally, this equation is simpler than the LenardBalescutype kinetic equation derived by Heyvaerts (2010), and it could be useful as a first step before considering more complicated effects. Its drawback is to ignore collective effects but this may not be crucial as we have explained (as suggested by Weinberg’s work, collective effects become important only when the system is close to instability). This Landautype equation was previously derived for systems with arbitrary longrange interactions^{9} in various dimensions of space (Chavanis 2007, 2008a,b, 2010) but we think that it is important to discuss these results in the specific case of selfgravitating systems with complements and amplification.
The paper is organized as follows. In Sect. 2, we study the dynamical evolution of a spatially inhomogeneous stellar system as a whole. Starting from the BBGKY hierarchy issued from the Liouville equation, and neglecting collective effects, we derive a general kinetic equation valid at the order 1/N in a proper thermodynamic limit. For N → +∞, it reduces to the Vlasov equation. At the order 1/N we recover the generalized Landau equation derived by Kandrup (1981) from a more abstract projection operator formalism. This equation is free of divergence at large scales since spatial inhomogeneity has been properly accounted for. Making a local approximation and introducing an upper cutoff at the Jeans length, we recover the standard VlasovLandau equation which is usually derived from a kinetic theory based on twobody encounters. Our approach provides an alternative derivation of this fundamental equation from the more rigorous Liouville equation. It has therefore some pedagogical interest. In Sect. 3, we study the relaxation of a test star in a steady distribution of field stars. We derive the corresponding FokkerPlanck equation and determine the expressions of the diffusion and friction coefficients. We emphasize the difference between the friction by polarization and the total friction (this difference may have been overlooked in previous works). For a thermal bath, we derive the Einstein relation between the diffusion and friction coefficients and obtain the explicit expression of the diffusion tensor. This returns the standard results obtained from the twobody encounters theory but, again, our presentation is different and offers an alternative derivation of these important results. For that reason, we give a short review of the basic formulae. In Sect. 4, we derive a Landautype kinetic equation written in angleaction variables and discuss its main properties. This equation, which does not make the local approximation, applies to fully inhomogeneous stellar systems and is free of divergence at large scales. We also develop a test particle approach and derive the corresponding FokkerPlanck equation in angleaction variables. Explicit expressions are given for the diffusion tensor and friction force, and they are compared with previous expressions obtained in the literature.
2. Evolution of the system as a whole
2.1. The BBGKY hierarchy
We consider an isolated system of N stars with identical mass m in Newtonian interaction. Their dynamics is fully described by the Hamilton equations (1)This Hamiltonian system conserves the energy E = H, the mass M = Nm, and the angular momentum L = ∑ _{i}mr_{i} × v_{i}. As recalled in the Introduction, stellar systems cannot reach a statistical equilibrium state in a strict sense. In order to understand their evolution, it is necessary to develop a kinetic theory.
We introduce the Nbody distribution function P_{N}(r_{1}, v_{1}, ..., r_{N}, v_{N}, t) giving the probability density of finding at time t the first star with position r_{1} and velocity v_{1}, the second star with position r_{2} and velocity v_{2}, etc. Basically, the evolution of the Nbody distribution function is governed by the Liouville equation (2)where (3)is the gravitational force by unit of mass experienced by the ith star due to its interaction with the other stars. Here, Φ_{d}(r) denotes the exact gravitational potential produced by the discrete distribution of stars and F(j → i) denotes the exact force by unit of mass created by the jth star on the ith star. The Liouville equation (Eq. (2)), which is equivalent to the Hamilton equations (Eq. (1)), contains too much information to be exploitable. In practice, we are only interested in the evolution of the onebody distribution P_{1}(r, v, t).
From the Liouville equation we can construct the complete BBGKY hierarchy for the reduced distribution functions (4)where the notation x stands for (r, v). The generic term of this hierarchy reads (5)This hierarchy of equations is not closed since the equation for the onebody distribution P_{1}(x_{1}, t) involves the twobody distribution P_{2}(x_{1}, x_{2}, t), the equation for the twobody distribution P_{2}(x_{1}, x_{2}, t) involves the threebody distribution P_{3}(x_{1}, x_{2}, x_{3}, t), and so on.
It is convenient to introduce a cluster representation of the distribution functions. Specifically, we can express the reduced distribution P_{j}(x_{1}, ..., x_{j}, t) in terms of products of distribution functions P_{j′ < j}(x_{1}, ..., x_{j′}, t) of lower order plus a correlation function (see e.g. Eqs. (8) and (9) below). Considering the scaling of the terms in each equation of the BBGKY hierarchy, we can see that there exist solutions of the whole BBGKY hierarchy such that the correlation functions scale as 1/N^{j − 1} in the proper thermodynamic limit N → +∞ defined in Appendix A. This implicitly assumes that the initial condition has no correlation, or that the initial correlations respect this scaling^{10}. If this scaling is satisfied, we can consider an expansion of the BBGKY hierarchy in terms of the small parameter 1/N. This is similar to the expansion of the BBGKY hierarchy in plasma physics in terms of the small parameter 1/Λ, where Λ ≫ 1 represents the number of charges in the Debye sphere (Balescu 2000). However, in plasma physics, the system is spatially homogeneous (because of Debye shielding which restricts the range of interaction) while, for stellar systems, spatial inhomogeneity must be taken into account. This brings additional terms in the BBGKY hierarchy that are absent in plasma physics.
2.2. The truncation of the BBGKY hierarchy at the order 1/N
The first two equations of the BBGKY hierarchy are We decompose the two and threebody distributions in the form where is the correlation function of order j. Substituting Eqs. (8) and (9) in Eqs. (6) and (7), and simplifying some terms, we obtain (10)(11)Equations (10) and (11) are exact for all N, but they are not closed. As explained previously, we shall close these equations at the order 1/N in the thermodynamic limit N → +∞. In this limit P_{1} ~ 1, , and . On the other hand, we can introduce dimensionless variables such that r ~ v ~ t ~ m ~ 1 and F(i → j) ~ G ~ 1/N (see Appendix A).
The advection term in the lefthand side (l.h.s.) of Eq. (10) is of order 1, and the collision term C in the righthand side (r.h.s.) is of order 1/N. Let us now consider the terms in Eq. (11) one by one. The first four terms correspond to the Liouville equation. The Liouville operator ℒ = ℒ_{0} + ℒ′ + ℒ_{m.f.} describes the complete twobody problem, including the inertial motion, the interaction between the stars (1,2) and the mean field produced by the other stars. The terms ℒ_{0} and ℒ_{m.f.} are of order 1/N while the term ℒ′ is of order 1/N^{2}. Therefore, the interaction term ℒ′ can a priori^{11} be neglected in the Liouville operator. This corresponds to the weak coupling approximation where only the mean field term ℒ_{0} + ℒ_{m.f.} is retained. The fifth term in Eq. (11) is a source term expressible in terms of the onebody distribution; it is of order 1/N. If we consider only the mean field Liouville operator ℒ_{0} + ℒ_{m.f.} and the source term , as we shall do in this paper, we can obtain a kinetic equation for stellar systems that is the counterpart of the Landau equation in plasma physics. The sixth term is of order 1/N and it corresponds to collective effects (i.e. the dressing of the particles by the polarization cloud). In plasma physics, this term leads to the LenardBalescu equation. It takes into account dynamical screening and regularizes the divergence at large scales that appears in the Landau equation. In the case of stellar systems, there is no largescale divergence because of the spatial inhomogeneity of the system. Therefore, collective effects are less crucial in the kinetic theory of stellar systems than in plasma physics. However, this term has been properly taken into account by Heyvaerts (2010) who obtained a kinetic equation of stellar systems that is the counterpart of the LenardBalescu equation in plasma physics. The last two terms are of the order 1/N^{2} and they will be neglected. In particular, the threebody correlation function , of order 1/N^{2}, can be neglected at the order 1/N. In this way, the hierarchy of equations is closed and a kinetic equation involving only twobody encounters can be obtained.
If we introduce the notations f = NmP_{1} (distribution function) and (twobody correlation function), we get at the order 1/NWe have introduced the mean force (by unit of mass) created on star 1 by all the other stars (14)and the fluctuating force (by unit of mass) created by star 2 on star 1 (15)Equations (12) and (13) are exact at the order 1/N. They form the right basis to develop the kinetic theory of stellar systems at this order of approximation. Since the collision term in the r.h.s. of Eq. (12) is of order 1/N, we expect that the relaxation time of stellar systems scales as ~ Nt_{D} where t_{D} is the dynamical time. As we shall see, the discussion is more complicated because of the presence of logarithmic corrections in the relaxation time and the absence of a strict statistical equilibrium state.
2.3. The limit N → +∞: the Vlasov equation (collisionless regime)
In the limit N → +∞, for any fixed interval of time [0, T] , the correlations between stars can be neglected. Therefore, the mean field approximation becomes exact and the Nbody distribution function factorizes in N onebody distribution functions (16)Substituting this factorization in the Liouville equation, and integrating over x_{2}, x_{3}, ..., x_{N}, we find that the smooth distribution function f(r, v, t) = NmP_{1}(r, v, t) is the solution of the Vlasov equation (17)This equation also results from Eq. (12) if we neglect the correlation function g(1,2) in the r.h.s. and replace N − 1 by N.
The Vlasov equation describes the collisionless evolution of stellar systems for times shorter than the relaxation time ~ Nt_{D}. In practice N ≫ 1 so that the domain of validity of the Vlasov equation is huge (see the end of Sect. 2.8). As recalled in the Introduction, the VlasovPoisson system develops a process of phase mixing and violent relaxation leading to a QSS on a very short timescale, of the order of a few dynamical times t_{D}. Elliptical galaxies are in such QSSs. LyndenBell (1967) developed a statistical mechanics of the Vlasov equation in order to describe this process of violent relaxation and predict the QSS achieved by the system. Unfortunately, the predictions of his statistical theory are limited by the problem of incomplete relaxation. Kinetic theories of violent relaxation, which may account for incomplete relaxation, have been developed by Kadomtsev & Pogutse (1970), Severne & Luwel (1980), and Chavanis (1998, 2008a,b).
2.4. The order O(1/N): the generalized Landau equation (collisional regime)
If we neglect strong collisions and collective effects, the first two equations of the BBGKY hierarchy (Eqs. (12) and (13)) reduce to
(18)(19)The first equation gives the evolution of the onebody distribution function. The l.h.s. corresponds to the (Vlasov) advection term. The r.h.s. takes into account correlations (finite N effects, graininess, discreteness effects) between stars that develop as a result of their interactions. These correlations correspond to encounters (collisions).
Equation (19) may be viewed as a linear first order differential equation in time. It can be symbolically written as (20)where ℒ = ℒ_{0} + ℒ_{m.f.} is a mean field Liouville operator and is a source term expressible in terms of the onebody distribution. This equation may be solved by the method of characteristics. Introducing the Green function (21)constructed with the mean field Liouville operator ℒ, we obtain (22)where we have assumed that no correlation is present initially so that g(x, x_{1}, t = 0) = 0; if correlations are present initially, it can be shown that they are rapidly washed out^{12}. Substituting Eq. (22) in Eq. (18), we obtain (23)In writing this equation, we have adopted a Lagrangian point of view. The coordinates r_{i} appearing after the Green function must be viewed as and . Therefore, in order to evaluate the integral of Eq. (23), we must move the stars following the trajectories determined by the selfconsistent mean field.
The kinetic equation (Eq. (23)) is valid at the order 1/N so it describes the collisional evolution of the system (ignoring collective effects) on a timescale of order Nt_{D}. Equation (23) is a nonMarkovian integrodifferential equation. It takes into account delocalizations in space and time (i.e. spatial inhomogeneity and memory effects). Actually, the Markovian approximation is justified in the N → +∞ limit because the timescale ~ Nt_{D} over which f(r, v, t) changes is long compared to the correlation time τ_{corr} ~ t_{D} over which the integrand in Eq. (23) has significant support^{13}. Therefore, we can compute the correlation function (Eq. (22)) by assuming that the distribution function is frozen at time t. This corresponds to the Bogoliubov ansatz in plasma physics. If we replace f(r, v, t − τ) and f(r_{1}, v_{1}, t − τ) by f(r, v, t) and f(r_{1}, v_{1}, t) in Eq. (23) and extend the time integral to +∞, we obtain (24)Similarly, we can compute the trajectories of the stars by assuming that the mean field is independent on τ and equal to its value at time t so that and .
The structure of the kinetic equation (Eq. (24)) has a clear physical meaning. The l.h.s. corresponds to the Vlasov advection term arising from mean field effects. The r.h.s. can be viewed as a collision operator C_{N} [f] taking finite N effects into account. For N → +∞, it vanishes and we recover the Vlasov equation. For finite N, it describes the cumulative effect of binary collisions between stars. The collision operator is a sum of two terms, a diffusion term and a friction term. The coefficients of diffusion and friction are given by generalized Kubo formulae, i.e. they involve the time integral of the autocorrelation function of the fluctuating force. The kinetic equation (Eq. (24)) bears some resemblance with the FokkerPlanck equation. However, it is more complicated since it is an integrodifferential equation, not a differential equation (see Sect. 3).
Equation (24) may be viewed as a generalized Landau equation. Since the spatial inhomogeneity and the finite extension of the system are properly taken into account, there is no divergence at large scales. There is, however, a logarithmic divergence at small scales which is due to the neglect of the interaction term ℒ′ in the Liouville operator (see Sect. 2.6). At large scales (i.e. for large impact parameters), this term can be neglected and the trajectories of the stars are essentially due to the mean field. Thus, we can make the weak coupling approximation leading to the Landau equation. This approximation describes weak collisions. However, at small scales (i.e. for small impact parameters), we can no longer ignore the interaction term ℒ′ in the Liouville operator and we have to solve the classical twobody problem. This is necessary to correctly describe strong collisions for which the trajectory of the particles deviates strongly from the mean field motion. When the mean field Green function G (constructed with ℒ_{0} + ℒ_{m.f.}) is replaced by the total Green function G′ (constructed with ℒ_{0} + ℒ′ + ℒ_{m.f.}), taking into account the interaction term, the generalized Landau equation (Eq. (24)) is replaced by a more complicated equation which can be viewed as a generalized Boltzmann equation. This equation is free of divergence (since it takes both spatial inhomogeneity and strong collisions into account), but it is unnecessarily complicated because it does not exploit the dominance of weak collisions over rare strong collisions for the gravitational potential. Indeed, a star suffers a large number of weak distant encounters and very few close encounters. A better practical procedure is to use the generalized Landau equation (Eq. (24)) with a cutoff at small scales in order to take into account our inability to describe strong collisions by this approach (see Sect. 2.6).
The generalized Landau equation (Eq. (24)) was derived by Kandrup (1981) from the Liouville equation by using the projection operator formalism. It can also be derived from the Liouville equation by using the BBGKY hierarchy or from the Klimontovich equation by making a quasilinear approximation (Chavanis 2008a,b).
2.5. The VlasovLandau equation
Selfgravitating systems are intrinsically spatially inhomogeneous. However, the collision operator at position r can be simplified by making a local approximation and performing the integrals as if the system were spatially homogeneous with the density ρ = ρ(r). This amounts to replacing f(r_{1}, v_{1}, t) by f(r, v_{1}, t) and by F(i → j) in Eq. (24). This local approximation is motivated by the work of Chandrasekhar & von Neumann (1942) who showed that the distribution of the gravitational force is a Lévy law (called the Holtzmark distribution) dominated by the contribution of the nearest neighbor. With this local approximation, Eq. (24) becomes (25)where we have used F(0 → 1) = −F(1 → 0). The Green function G_{0} corresponds to the free motion of the particles associated with the Liouville operator ℒ_{0}. Using Eqs. (C.1) and (C.2), the foregoing equation can be rewritten as (26)where F(1 → 0, t − τ) is expressed in terms of the Lagrangian coordinates. The integrals over τ and r_{1} can be calculated explicitly (see Appendices C and D). We then find that the evolution of the distribution function is governed by the VlasovLandau equation (27)where we have noted w = v − v_{1}, f = f(r, v, t), f_{1} = f(r, v_{1}, t), and where (2π)^{3}û(k) = −4πG/k^{2} represents the Fourier transform of the gravitational potential. Under this form, we see that the collisional evolution of a stellar system is due to a condition of resonance k·v = k·v′ (with v ≠ v′) encapsulated in the δfunction. This δfunction expresses the conservation of energy.
The VlasovLandau equation may also be written as (see Appendix C) where (30)is the Coulombian logarithm that has to be regularized with appropriate cutoffs (see Sect. 2.6). The r.h.s. of Eq. (28) is the original form of the collision operator given by Landau (1936) for the Coulombian interaction^{14}. It applies to weakly coupled plasmas. We note that the potential of interaction only appears in the constant A which merely determines the relaxation time. The structure of the Landau equation is independent of the potential. The Landau equation was originally derived from the Boltzmann equation in the limit of weak deflections Δv ≪ 1 (Landau 1936)^{15}. In the case of plasmas, the system is spatially homogeneous and the advection term is absent in Eq. (28). In the case of stellar systems, when we make the local approximation, the spatial inhomogeneity of the system is only retained in the advection term of Eq. (28). This is why this kinetic equation is referred to as the VlasovLandau equation. This is the fundamental kinetic equation of stellar systems.
2.6. Heuristic regularization of the divergences
To obtain the VlasovLandau equation (Eq. (28)), we have made a local approximation. This amounts to calculating the collision operator at each point as if the system were spatially homogeneous. As a result of this homogeneity assumption, a logarithmic divergence appears at large scales in the Coulombian logarithm (Eq. (30)). In plasma physics, this divergence is cured by the Debye shielding. A charge is surrounded by a polarization cloud of opposite charges which reduces the range of the interaction. When collective effects are properly taken into account, as in the LenardBalescu equation, no divergence appears at large scales and the Debye length arises naturally. Heuristically, we can use the Landau equation and introduce an upper cutoff at the Debye length λ_{D} ~ (k_{B}T/ne^{2})^{1/2}. For selfgravitating systems, there is no shielding and the divergence is cured by the finite extent of the system. The interaction between two stars is only limited by the size of the system. When spatial inhomogeneity is taken into account, as in the generalized Landau equation (Eq. (24)), no divergence occurs at large scales. Heuristically, we can use the VlasovLandau equation (Eq. (28)) and introduce an upper cutoff at the Jeans length λ_{J} ~ (k_{B}T/Gm^{2}n)^{1/2} which gives an estimate of the system’s size R.
The Coulombian logarithm (Eq. (30)) also diverges at small scales. As explained previously, this is due to the neglect of strong collisions that produce important deflections. Indeed, for collisions with small impact parameter, the mean field approximation is clearly irrelevant and it is necessary to solve the twobody problem exactly (see Appendix H). Accordingly, the smallscale divergence is cured by the proper treatment of strong collisions. Heuristically, we can use the Landau equation and introduce a lower cutoff at the gravitational Landau length (the gravitational analogue of the Landau length in plasma physics) which corresponds to the impact parameter leading to a deflection at 90°.
Introducing a largescale cutoff at the Jeans length λ_{J} and a smallscale cutoff at the Landau length λ_{L}, and noting that , we find that the Coulombian logarithm scales as where is the number of stars in the cluster^{16}.
2.7. Properties of the VlasovLandau equation
The VlasovLandau equation conserves the mass M = ∫f drdv and the energy . It also monotonically increases the Boltzmann entropy S = −∫(f/m)ln (f/m) drdv in the sense that Ṡ ≥ 0 (Htheorem). As a result of the local approximation, the proof of these properties is the same^{17} as for the spatially homogeneous Landau equation (Balescu 2000). Because of these properties, we might expect that a stellar system will relax towards the Boltzmann distribution which maximizes the entropy at fixed mass and energy. However, we know that there is no maximum entropy state for an unbounded selfgravitating system (the Boltzmann distribution has infinite mass). Therefore, the VlasovLandau equation does not relax towards a steady state and the entropy does not reach a stationary value. Actually, the entropy increases permanently as the system evaporates. But since evaporation is a slow process, the system may achieve a QSS that is close to the Boltzmann distribution. A typical quasistationary distribution is the MichieKing model (31)where ϵ = v^{2}/2 + Φ(r) is the energy and j = r × v the angular momentum. This distribution takes into account the escape of high energy stars and the anisotropy of the velocity distribution. It can be derived, under some approximations, from the VlasovLandau equation by using the condition that f = 0 if the energy of the star is larger than the escape energy ϵ_{m} (Michie 1963; King 1965). The MichieKing distribution reduces to the isothermal distribution f ∝ e^{−βmϵ} for low energies. In this sense, we can define a relaxation time for a stellar system. From the VlasovLandau equation (Eq. (28)), we find that the relaxation time scales as (32)where v_{m} is the root mean square (rms) velocity of the stars. Introducing the dynamical time t_{D} ~ λ_{J}/v_{m} ~ R/v_{m}, we obtain the scaling (33)The fact that the ratio between the relaxation time and the dynamical time depends only on the number of stars and scales as N/lnN was first noted by Chandrasekhar (1942).
A simple estimate of the evaporation time was calculated by Ambartsumian (1938) and Spitzer (1940) and gives t_{evap} ≃ 136t_{R}. More precise values have been obtained by studying the evaporation process in an artificially uniform medium (Chandrasekhar 1943b; Spitzer & Härm 1958; Michie 1963; King 1965; Lemou & Chavanis 2010) or in a more realistic inhomogeneous cluster (Hénon 1961). Since t_{evap} ≫ t_{R}, we can consider that the system relaxes towards a steady distribution given by Eq. (31) on a timescale t_{R} and that this distribution slowly evolves on a longer timescale as the stars escape^{18}. The characteristic time in which the system’s stars evaporate is t_{evap}. The evaporation is one reason for the evolution of stellar systems. However, as demonstrated by Antonov (1962) and LyndenBell & Wood (1968), stellar systems may evolve more rapidly because of gravothermal catastrophe. In that case, the MichieKing distribution changes significatively as a result of core collapse. This evolution has been described by Cohn (1980), and it leads ultimately to the formation of a binary star surrounded by a hot halo. A configuration of this type can have arbitrarily large entropy. Cohn (1980) found that the entropy increases permanently during core collapse, confirming that the VlasovLandau equation has no equilibrium state.
Even if the system were confined within a small box so as to prevent both the evaporation and the gravothermal catastrophe, there would be no statistical equilibrium state in a strict sense because there is no global entropy maximum (Antonov 1962). A configuration in which some subset of the particles are tightly bound together (e.g. a binary star), and in which the remaining particles share the energy thereby released, may have an arbitrarily large entropy. However, such configurations, which require strong correlations, are generally reached very slowly (on a timescale much larger than (N/lnN)t_{D}) as a result of encounters involving many particles. To describe these configurations, one would have to take high order correlations into account in the kinetic theory. These configurations may be relevant in systems with a small number of stars (Chabanol et al. 2000), but when N is large the picture is different. On a timescale of the order of (N/lnN)t_{D} the onebody distribution function is expected to reach the Boltzmann distribution which is a local entropy maximum. This state is metastable, but its lifetime is expected to be very large, scaling as e^{N}, so that it is stable in practice (Chavanis 2005, 2006). In this sense, there exist true statistical equilibrium states for selfgravitating systems confined within a small box. However, we may argue that this situation is highly artificial.
Finally, using very different methods based on ergodic theory, Gurzadyan & Savvidy (1986) argue that, because of collective behaviour, the relaxation time for stellar systems scales like intermediate between the dynamical time t_{D} (violent relaxation) and the binary relaxation time t_{R} ~ (N/lnN)t_{D}.
2.8. Dynamical evolution of stellar systems: a short review
Using the kinetic theory, we can identify different phases in the dynamical evolution of stellar systems.
A selfgravitating system initially outofmechanical equilibrium undergoes a process of violent collisionless relaxation towards a virialized state. In this regime, the dynamical evolution of the cluster is described by the VlasovPoisson system. The phenomenology of violent relaxation has been described by Hénon (1964), King (1966), and LyndenBell (1967). Numerical simulations that start from cold and clumpy initial conditions generate a QSS that fits the de Vaucouleurs R^{1/4} law for the surface brightness of elliptical galaxies quite well (van Albada 1982). The inner core is almost isothermal (as predicted by LyndenBell 1967) while the velocity distribution in the envelope is radially anisotropic and the density profile decreases as r^{4}. One success of LyndenBell’s statistical theory of violent relaxation is to explain the isothermal core of elliptical galaxies without recourse to collisions. In contrast, the structure of the halo cannot be explained by LyndenBell’s theory as it results from an incomplete relaxation. Models of incompletely relaxed stellar systems have been elaborated by Bertin & Stiavelli (1984), Stiavelli & Bertin (1987), and Hjorth & Madsen (1991). These theoretical models nicely reproduce the results of observations and numerical simulations (Londrillo et al. 1991; Trenti et al. 2005). In the simulations, the initial condition needs to be sufficiently clumpy and cold to generate enough mixing required for a successful application of the statistical theory of violent relaxation. Numerical simulations starting from homogeneous spheres (see e.g. Roy & Perez 2004; Levin et al. 2008; Joyce et al. 2009) show little angular momentum mixing and lead to different results. In particular, they display a larger amount of mass loss (evaporation) than simulations starting from clumpy initial conditions. Clumps thus help the system to reach a universal final state from a variety of initial conditions, which can explain the similarity of the density profiles observed in elliptical galaxies.
On longer timescales, encounters between stars must be taken into account and the dynamical evolution of the cluster is governed by the VlasovLandauPoisson system. The first stage of the collisional evolution is driven by evaporation. Because of a series of weak encounters, the energy of a star can gradually increase until it reaches the local escape energy; in that case the star leaves the system. Numerical simulations (Spitzer 1987; Binney & Tremaine 2008) show that during this regime the system reaches a QSS that slowly evolves in amplitude as a result of evaporation as the system loses mass and energy. This quasistationary distribution function is close to the MichieKing model (Eq. (31)). The system has a corehalo structure. The core is isothermal while the stars in the outer halo move in predominantly radial orbits. Therefore, the distribution function in the halo is anisotropic. The density follows the isothermal law ρ ~ r^{2} in the central region (with a core of almost uniform density) and decreases as ρ ~ r^{−7/2} in the halo (for an isolated cluster with ϵ_{m} = 0) or tends to zero as ρ ∝ (r_{t} − r)^{3/2} (for a cluster with a finite tidal radius). Because of evaporation, the halo expands while the core shrinks as required by energy conservation. During the evaporation process, the central density increases permanently. At some point in the evolution, the system undergoes an instability related to the Antonov (1962) instability^{19} and the gravothermal catastrophe sets in (LyndenBell & Wood 1968). This instability is due to the negative specific heat of the inner system that evolves by losing energy and thereby growing hotter. The energy lost is transferred outward by stellar encounters. Hence the temperature always decreases outward, and the center continually loses energy, shrinks, and heats up. This leads to core collapse. Mathematically speaking, core collapse would generate a finite time singularity. When the evolution is modeled by the orbitaveraged FokkerPlanck equation, Cohn (1980) finds that the collapse is selfsimilar, that the central density becomes infinite in a finite time, and that the density behaves as ρ ~ r^{2.23}. The invariant profile found by Cohn differs from the MichieKing distribution (for which ρ ~ r^{2}) beyond a radius of about 10r_{core}. Larson (1970) and LyndenBell & Eggleton (1980) find similar results by modeling the evolution of the system by fluid equations. In all cases, the authors find a singular density profile at the collapse time that is integrable at r = 0. This means that “the singularity contains no mass”^{20}. In reality, if we come back to the Nbody system, there is no singularity and core collapse is arrested by the formation of binary stars caused by threebody collisions. These binaries can release sufficient energy to stop the collapse and even drive a reexpansion of the cluster in a postcollapse regime (Inagaki & LyndenBell 1983; Goodman 1984). Then, in principle, a series of gravothermal oscillations should follow (Bettwieser & Sugimoto 1984).
At the present epoch, small groups of stars such as globular clusters (N ~ 10^{5}, t_{D} ~ 10^{5} yr, age ~ 10^{10} yr, t_{R} ~ 10^{10} yr) are in the collisional regime. They are either in QSSs described by the MichieKing model or experiencing core collapse. By contrast, large clusters of stars like elliptical galaxies (N ~ 10^{11}, t_{D} ~ 10^{8} yr, age ~ 10^{10} yr, t_{R} ~ 10^{19} yr) are still in the collisionless regime and their apparent organization is a result of an incomplete violent relaxation.
3. Test star in a thermal bath
3.1. The FokkerPlanck equation
We now consider the relaxation of a test star (i.e. a tagged particle) evolving in a steady distribution of field stars^{21}. Because of the encounters with the field stars, the test star has a stochastic motion. We call P(r, v, t) the probability density of finding the test star at position r with velocity v at time t. The evolution of P(r, v, t) can be obtained from the generalized Landau equation (Eq. (24)) by considering that the distribution function of the field stars is fixed. Therefore, we replace f(r, v, t) by P(r, v, t) and f(r_{1}, v_{1}, t) by f(r_{1}, v_{1}) where f(r_{1}, v_{1}) is the steady distribution of the field stars. This procedure transforms the integrodifferential equation (Eq. (24)) into the differential equation (34)where ⟨ F ⟩ (r) = −∇Φ(r) is the static mean force created by the field stars with density ρ(r_{1}). This equation does not present any divergence at large scales. We can understand the above procedure as follows (see Chavanis 2012b, for more details). Equations (24) and (34) govern the evolution of the distribution function of a test star (described by the coordinates r and v) interacting with field stars (described by the running coordinates r_{1} and v_{1}). In Eq. (24), all the stars are equivalent so the distribution of the field stars f(r_{1}, v_{1}, t) changes with time exactly like the distribution of the test star f(r, v, t). In Eq. (34), the test star and the field stars are not equivalent since the field stars form a bath. The field stars have a steady given distribution f(r_{1}, v_{1}) while the distribution of the test star P(r, v, t) changes with time.
If we make a local approximation and use the VlasovLandau equation (Eq. (26)), we obtain (35)Denoting the advection operator by d/dt, Eq. (35) can be written in the form of a FokkerPlanck equation (36)involving a diffusion tensor (37)and a friction force (38)If we start directly from Eq. (27), which amounts to performing the integrals over τ and r_{1} in the previous expressions, we obtain the FokkerPlanck equation (39)with the diffusion and friction coefficients Using Eq. (40), we can easily establish the identity (42)where D = D^{μμ} = tr(D).
The diffusion tensor D^{μν} results from the fluctuations of the gravitational force caused by the granularities in the distribution of the field stars. It can be derived directly from the formula (see Appendix G) (43)obtained from Eq. (45a). The friction force F_{pol} results from the retroaction of the field stars to the perturbation caused by the test star, as in a polarization process. It can be derived from a linear response theory (Marochnik 1968; Kalnajs 1971b; Kandrup 1983; Chavanis 2008b). It will be called the friction by polarization to distinguish it from the total friction (see below). Equations (35)–(41) have been derived within the local approximation. More general formulae, valid for fully inhomogeneous stellar systems, are given in Kandrup (1983) and Chavanis (2008b). The friction force has also been calculated by Tremaine & Weinberg (1984), Bekenstein & Maoz (1992), Maoz (1993), and Nelson & Tremaine (1999) using different approaches.
Since the diffusion tensor depends on the velocity v of the test star, it is useful to rewrite Eq. (36) in a form that is fully consistent with the general FokkerPlanck equation (44)with (45)By identification, we find that (46)Therefore, when the diffusion coefficient depends on the velocity, the total friction is different from the friction by polarization. Substituting Eqs. (40) and (41) in Eq. (46), and using an integration by parts, we find that the diffusion and friction coefficients may be written as These expressions can be obtained directly from the equations of motion by expanding the trajectories of the stars in powers of 1/N in the limit N → +∞ (Chavanis 2008b). We recall that Eqs. (47) and (48) display a logarithmic divergence at small and large scales that must be regularized by introducing proper cutoffs as explained in Sect. 2.6. In astrophysics, the diffusion and friction coefficients of a star were first calculated by Chandrasekhar (1943a) from a twobody encounters theory (see also Cohen et al. 1950; Gasiorowicz et al. 1956; and Rosenbluth et al. 1957). The expressions obtained by these authors are different from those given above, but they are equivalent (see Sect. 3.5). In plasma physics, the diffusion and friction coefficients of a charge were first calculated by Hubbard (1961a) who took collective effects into account, thereby eliminating the divergence at large scales. When collective effects are neglected, his expressions reduce to Eqs. (47) and (48). On the other hand, strong collisions have been taken into account by Chandrasekhar (1943a) in astrophysics and by Hubbard (1961b) in plasma physics. In that case, there is no divergence at small impact parameters in the diffusion and friction coefficients, and the Landau length appears naturally (see Appendix H).
The two forms of the FokkerPlanck equation (Eqs. (36) and (44)) have their own interest. The expression (Eq. (44)) where the diffusion coefficient is placed after the two derivatives ∂^{2}(DP) involves the total friction force F_{friction} and the expression (Eq. (36)) where the diffusion coefficient is placed between the derivatives ∂D∂P involves the friction by polarization F_{pol}. Astrophysicists are used to the form in Eq. (44). However, it is the form in Eq. (36) that stems from the Landau equation (Eq. (26)). We shall come back to this observation in Sect. 3.5.
From Eqs. (40) and (41), we easily obtain (49)Combining Eq. (46) with Eq. (49), we get (50)Therefore, the friction force F_{friction} is equal to twice the friction by polarization F_{pol} (for a test star with mass m interacting with field stars with mass m_{f}, this factor two is replaced by (m + m_{f})/m; see Appendix F). This explains the difference of factor 2 in the calculations of Chandrasekhar (1943a) who determined F_{friction} and in the calculations of Kalnajs (1971b) and Kandrup (1983) who determined F_{pol}.
3.2. The Einstein relation
In the central region of the system, the distribution of the field stars is close to the MaxwellBoltzmann distribution^{22}(51)where β = 1/(k_{B} T) is the inverse temperature and ρ(r_{1}) ∝ e^{−βmΦ(r1)} is the density given by the Boltzmann law. Therefore, the field stars form a thermal bath. Substituting the identity (52)in Eq. (41), using the δfunction to replace k·v_{1} by k·v, and comparing the resulting expression with Eq. (40), we find that (53)The friction coefficient is given by an Einstein relation expressing the fluctuationdissipation theorem^{23}. We emphasize that the Einstein relation is valid for the friction force by polarization F_{pol}, not for the total friction F_{friction} (we do not have this subtlety in standard Brownian theory where the diffusion coefficient is constant). Using Eq. (43), we can rewrite Eq. (53) in the form (54)This relation is usually called the Kubo formula. More general expressions of the Kubo formula valid for fully inhomogeneous stellar systems are given in Kandrup (1983) and Chavanis (2008b). Using the Einstein relation, the FokkerPlanck equation (Eq. (36)) takes the form (55)where the diffusion coefficient is given by Eq. (40) with Eq. (51). This equation is similar to the Kramers equation in Brownian theory (Kramers 1940) except that the diffusion coefficient is a tensor and that it depends on the velocity of the test star. For an isotropic distribution function (e.g. the Maxwellian), it can be put in the form (56)where D_{ ∥ } and D_{ ⊥ } are the diffusion coefficients in the directions parallel and perpendicular to the velocity of the test star. We note that D ≡ D^{μμ} = D_{ ∥ } + D_{ ⊥ }. Using D^{μν}v^{ν} = D_{ ∥ }v^{μ}, the friction by polarization (Eq. (53)) can be written as (57)The friction is proportional and opposite to the velocity of the test star, and the friction coefficient is given by the Einstein relation ξ = D_{ ∥ }βm. The total friction is (58)This is the Chandrasekhar dynamical friction.
The steady state of the FokkerPlanck equation (Eq. (55)) is the MaxwellBoltzmann distribution (Eq. (51)). Since the FokkerPlanck equation admits an Htheorem for the Boltzmann free energy (Risken 1989), one can prove that P(r, v, t) converges towards the MaxwellBoltzmann distribution (Eq. (51)) for t → +∞. In other words, the test star acquires the distribution of the field stars (thermalization). We recall, however, that for selfgravitating systems, these results cannot hold everywhere in the cluster since the MaxwellBoltzmann distribution does not exist globally.
If we assume that the distribution P(r, v, t) of the test star is isotropic, the FokkerPlanck equation becomes (59)or (60)It can be obtained from an effective Langevin equation (61)where R(t) is a Gaussian white noise satisfying ⟨ R(t) ⟩ = 0 and ⟨ R^{μ}(t)R^{ν}(t′) ⟩ = δ^{μν}δ(t − t′). Since the diffusion coefficient depends on the velocity, the noise is multiplicative. The force acting on the star consists of two parts, a part derivable from a smoothed distribution of matter −∇Φ, and a residual random part attributable to the fluctuations of that distribution. In turn, the random part can be described by a Gaussian white noise multiplied by a velocitydependent factor plus a friction force equal to the friction by polarization F_{pol} = −D_{ ∥ }(v)βmv and a spurious drift F_{spurious} = (1/2)(∂D_{ ∥ }/∂v) attributable to the multiplicative noise (we have used the Stratonovich representation). If we neglect the velocity dependence of the diffusion coefficient, we recover the results of Chandrasekhar (1943a, 1949) obtained in his Brownian theory.
3.3. The diffusion tensor for isothermal systems
When the velocity distribution of the field stars is given by the Maxwellian distribution (Eq. (51)), the diffusion tensor (Eq. (40)) may be calculated as follows. If we introduce the representation (62)for the δfunction in Eq. (40), the diffusion tensor can be rewritten as (63)where is the threedimensional Fourier transform of the velocity distribution. This equation can be directly obtained from Eq. (37) or Eq. (43) (see Appendix G). This shows that the auxiliary integration variable t in Eq. (63) represents the time. For the Maxwellian distribution (Eq. (51)), is a Gaussian. If we perform the integration over t (which is the onedimensional Fourier transform of a Gaussian), we find that the diffusion tensor can be expressed as (64)Alternatively, this expression can be obtained from Eq. (40) by introducing a cartesian system of coordinates for v_{1} with the zaxis taken along the direction of k, and performing the integration. With the notation , Eq. (64) can be rewritten as (65)where (66)We note that the potential of interaction only appears in a multiplicative constant that fixes the relaxation time (see below). Using Eq. (C.7), we get (67)where lnΛ is the Coulombian logarithm (Eq. (30)) and is the mean square velocity of the field stars. The diffusion tensor may be written as (68)where t_{R} is the local relaxation time defined below in Eq. (76). This relation emphasizes the scaling .
Introducing a spherical system of coordinates with the zaxis in the direction of x, we can easily compute the integral in Eq. (66). Then, we can write the normalized diffusion tensor in the form (69)where (70)with (71)The error function is defined by (72)We have the asymptotic behaviours We note that when x → 0. The longitudinal and transverse diffusion coefficients and the friction force are plotted in Fig. 1.
Fig. 1 Normalized diffusion coefficients G_{ ∥ }(x), G_{ ⊥ }(x) and friction force xG_{ ∥ }(x) for a thermal bath. The friction is maximum for x ≃ 0.97, i.e. when the velocity of the test star is approximately equal to 0.792 times the rms velocity of the field stars (its value is (xG_{ ∥ })_{max} = 2.383). 
3.4. The relaxation time
We can use the preceding results to estimate the relaxation time of the velocity distribution of the test particle towards the Maxwellian distribution (thermalization). If we set , the FokkerPlanck equation (Eq. (55)) can be rewritten as (75)where t_{R} is the local relaxation time (76)The prefactor is equal to 0.482 but, of course, this numerical factor may vary depending on the definition of the relaxation time. The relaxation time is inversely proportional to the local density ρ(r). Therefore, the relaxation time is shorter in regions of high density (core) and longer in regions of low density (halo). Introducing the dynamical time t_{D} = λ_{J}/v_{m}, we get (77)We note that the relaxation time of a test particle in a bath is of the same order as the relaxation time of the system as a whole (see Sect. 2.7).
We can also get an estimate of the relaxation time by the following argument (Spitzer 1987). If the diffusion coefficient were constant, the typical velocity of the test star (in one spatial direction) would increase as ⟨ (Δv)^{2} ⟩ /3 ~ 2D_{ ∥ }t. The relaxation time t_{r} is the typical time at which the typical velocity of the test star has reached its equilibrium value so that ⟨ (Δv)^{2} ⟩ (t_{r}) = ⟨ v^{2} ⟩ (+∞). Since D_{ ∥ } depends on v, the description of the diffusion process is more complex. However, the formula resulting from the previous arguments with D_{ ∥ } = D_{ ∥ }(v_{m}) should provide a good estimate of the relaxation time. Using Eq. (67) and comparing with Eq. (76) we obtain t_{r} = K_{3}t_{R}, where . Numerically, K_{3} = 0.13587547....
Finally, we can estimate the relaxation time by where ξ is the friction coefficient. Using the Einstein relation ξ = D_{ ∥ }βm (see Eq. (57)) with D_{ ∥ } = D_{ ∥ }(v_{m}) we find that .
3.5. The Rosenbluth potentials
It is possible to obtain simple expressions of the diffusion and friction coefficients for any isotropic distribution of the bath. If we start from the expression of the VlasovLandau equation (Eq. (28)), we find that the FokkerPlanck equation (Eq. (35)) can be written as (78)where K^{μν} is defined by Eq. (29). Comparing Eq. (78) with Eqs. (36) and (44), the diffusion and friction coefficients are given by Using the identities (81)and (82)the coefficients of diffusion and friction can be rewritten as where (85)are the socalled Rosenbluth potentials (Rosenbluth et al. 1957). In terms of these potentials, the FokkerPlanck equations (Eqs. (36) and (44)) may be rewritten as If the field particles have an isotropic velocity distribution, the Rosenbluth potentials take the particularly simple form (see e.g. Binney & Tremaine 2008) When g = g(v), the diffusion tensor (Eq. (83)) can be put in the form of Eq. (56) with (90)Using Eq. (89), we obtain On the other hand, when h = h(v), the friction force (Eq. (84)) can be written as (93)Using Eq. (88), we get (94)This expression can be obtained directly from Eq. (84) by noting (Binney & Tremaine 2008) that h(v) in Eq. (85) is similar to the gravitational potential Φ(r) produced by a distribution of mass ρ(r), where v plays the role of r and f(v) the role of ρ(r). Therefore, if f(v) is isotropic, Eq. (94) is equivalent to the expression of the gravitational field F = −GM(r)r/r^{3} produced by a spherically symmetric distribution of matter, where M(r) is the mass within the sphere of radius r (Gauss theorem). This formula shows that the friction is due only to field stars with a velocity less than the velocity of the test star. This observation was first made by Chandrasekhar (1943a).
The previous expressions for the diffusion and friction coefficients are valid for any isotropic distribution of the field particles. When f(v) is the Maxwell distribution, we recover the results of Sect. 3.2. If we substitute Eqs. (91)–(93) into Eq. (36) or (44), we get a FokkerPlanck equation describing the evolution of a test particle in a bath with a prescribed distribution f(v)^{24}. Alternatively, if we come back to the original Landau kinetic equation (Eq. (28)) which can be written as (95)assume an isotropic velocity distribution, use D^{μν}v^{ν} = D_{ ∥ }v^{μ}, and substitute the general expressions of the diffusion and friction coefficients (Eqs. (90)–(94)) with f = f(v,t), we obtain the integrodifferential equation (96)or, in more explicit form (97)describing the evolution of the system as a whole. Under this form, Eq. (97) applies to an artificial infinite homogeneous distribution of stars. This equation has been studied by King (1960) in his investigations on the evaporation of globular clusters. Within the local approximation, Eq. (97) also represents the simplification of the collision operator that occurs in the r.h.s. of the VlasovLandau equation (Eq. (28)) when the velocity distribution of the stars is isotropic. In that case, we must restore the space variable and the advection term in Eq. (97). From this equation, implementing an adiabatic approximation, we can derive the orbitaveraged FokkerPlanck equation which was used by Hénon (1961) and Cohn (1980) to study the collisional evolution of globular clusters. It reads (98)where (99)is proportional to the phase space volume available to stars with an energy less than ϵ.
3.6. Comparison with the twobody encounters theory
In the previous sections, we have derived the standard kinetic equations of stellar systems from the the Liouville equation by using the BBGKY hierarchy. In classical textbooks of astrophysics (Spitzer 1987; Heggie & Hut 2003; Binney & Tremaine 2008), these equations are derived in a different manner. One usually starts from the FokkerPlanck equation (Eq. (44)), makes a local approximation, and evaluates the diffusion tensor ⟨ Δv^{μ}Δv^{ν} ⟩ and the friction force ⟨ Δv^{μ} ⟩ by considering the mean effect of a succession of twobody encounters. This twobody encounters theory was pioneered by Chandrasekhar (1942) and further developed by Rosenbluth et al. (1957). This approach leads directly to the expressions in Eqs. (83) and (84) of the diffusion and friction coefficients of a test star in a bath of field stars. These expressions are then substituted in the FokkerPlanck equation (Eq. (44)). Finally, arguing that the field stars and the test star should evolve in the same manner, the FokkerPlanck equation is transformed into the integrodifferential equation (Eq. (97)) describing the evolution of the system as a whole (King 1960). When an adiabatic approximation is implemented (Hénon 1961), one finally obtains the orbitaveraged FokkerPlanck equation (Eq. (98)).
In this paper we have proceeded the other way round. Starting from the Liouville equation, using the BBGKY hierarchy, making a local approximation, and neglecting strong collisions, we have derived the VlasovLandau equation (Eq. (28)) describing the evolution of the system as a whole. Then, making a bath approximation^{25}, we have obtained the FokkerPlanck equation in the form of Eq. (36) with the diffusion and friction coefficients given by Eqs. (79) and (80). This equation can then be transformed into Eq. (44) with Eqs. (83) and (84). Our approach emphasizes the importance of the Landau equation in the kinetic theory of stellar systems, while this equation does not appear in the works of Chandrasekhar (1942), Rosenbluth et al. (1957), King (1960), and Hénon (1961), nor in the standard textbooks of stellar dynamics by Spitzer (1987), Heggie & Hut (2003), and Binney & Tremaine (2008).
Actually, the kinetic equation derived by these authors is equivalent to the Landau equation^{26}, but it is written in a different form. They write the FokkerPlanck equation in the form of Eq. (44) with the diffusion coefficient placed after the second derivative (∂^{2}D) while the Landau equation (Eq. (28)) is related to the FokkerPlanck equation (Eq. (36)) in which the diffusion coefficient is inserted between the first derivatives (∂D∂). This difference is important on a physical point of view for two reasons. First, the Landau equation isolates the friction by polarization F_{pol} while the equation derived by Chandrasekhar (1942) and Rosenbluth et al. (1957) involves the total friction F_{friction}. Secondly, the Landau equation has a nice symmetric structure from which we can immediately deduce all the conservation laws of the system (conservation of mass, energy, impulse, and angular momentum) and the Htheorem for the Boltzmann entropy (Balescu 2000). These properties are less apparent in the equations derived by King (1960) and Hénon (1961) for the evolution of the system as a whole. It is interesting to note that the symmetric structure of the kinetic equation was not realized by early stellar dynamicists while the Landau equation was known long before in plasma physics.
Finally, the approach based on the Liouville equation and on the BBGKY hierarchy is more rigorous than the twobody encounters theory because it relaxes the assumption of locality and does not produce any divergence at large scales^{27}. It leads to the generalized Landau equation (Eq. (24)) that is perfectly wellbehaved at large scales contrary to the VlasovLandau equation (Eq. (28)) in which a largescale cutoff has to be introduced in order to avoid a divergence. This divergence is due to the longrange nature of the gravitational potential which precludes a rigorous application of the twobody encounters theory that is valid for potentials with shortrange interactions. Actually, the twobody encounters theory is marginally applicable to the gravitational force (it only generates a weak logarithmic divergence at large scales) and this is why it is successful in practice. The generalized Landau equation (Eq. (24)) represents a conceptual improvement of the VlasovLandau equation (Eq. (28)) because it goes beyond the local approximation and takes fully into account the spatial inhomogeneity of the system. Unfortunately, this equation is very complicated to be of much practical use. It can, however, be simplified by using angleaction variables as we show in the next section.
4. Kinetic equations with angleaction variables
4.1. Adiabatic approximation
In order to deal with spatially inhomogeneous systems, it is convenient to introduce angleaction variables (Goldstein 1956; Binney & Tremaine 2008). Angleaction variables have been used by many authors in astrophysics in order to solve dynamical stability problems (Kalnajs 1977; Goodman 1988; Weinberg 1991; Pichon & Cannon 1997; Valageas 2006a) or to compute the diffusion and friction coefficients of a test star in a cluster (LyndenBell & Kalnajs 1972; Tremaine & Weinberg 1984; Binney & Lacey 1988; Weinberg 1998; 2001; Nelson & Tremaine 1999; Pichon & Aubert 2006; Valageas 2006b; Chavanis 2007, 2010). By construction, the Hamiltonian H in angle and action variables depends only on the actions J = (J_{1}, ..., J_{d}) that are constants of the motion; the conjugate coordinates w = (w_{1}, ..., w_{d}) are called the angles (see Appendix B). Therefore, any distribution of the form f = f(J) is a steady state of the Vlasov equation. According to the Jeans theorem, this is not the general form of Vlasov steady states. However, if the potential is regular, for all practical purposes, any timeindependent solution of the Vlasov equation may be represented by a distribution of the form f = f(J) (strong Jeans theorem).
We shall assume that the system has reached a QSS described by a distribution f = f(J) as a result of a violent collisionless relaxation involving only mean field effects. Because of finite N effects, the distribution function f slowly evolves in time. Finite N effects are taken into account in the collision operator appearing in the r.h.s. of Eq. (24). Since this term is of order 1/N, the effect of collisions (granularities, finite N effects, correlations) is a very slow process that takes place on a relaxation timescale t_{R} ~ (N/lnN)t_{D} (see Sect. 2.7). Therefore, there is a timescale separation between the dynamical time t_{D} that is the timescale during which the system reaches a steady state of the Vlasov equation through phase mixing and violent collisionless relaxation, and the collisional relaxation time t_{R} that is the timescale during which the system reaches an almost isothermal distribution because of finite N effects.
Because of this timescale separation, the distribution function is stationary on the dynamical timescale. It will evolve through a sequence of QSSs that are steady states of the Vlasov equation, depending only on the actions J, slowly changing in time as a result of the cumulative effect of encounters (finite N effects). Indeed, the system readjusts itself dynamically at each step of the collisional process. The distribution function averaged over a short dynamical timescale can be approximated by (100)Therefore, the distribution function is a function f = f(J, t) of the actions alone that slowly evolves in time under the effect of collisions. This is similar to an adiabatic approximation. The system is approximately in mechanical equilibrium at each stage of the dynamics and the collisions slowly drive it towards an almost isothermal distribution, corresponding to a quasithermodynamical equilibrium state.
4.2. Evolution of the system as a whole: a Landautype equation with angleaction variables
Introducing angleaction variables w and J, the generalized Landau equation (Eq. (24)) becomes^{28}(101)with (102)To obtain Eq. (101), we have averaged Eq. (24) over w (to simplify the expressions, the average ⟨ . ⟩ = (2π)^{3}∫dw is implicit), written the scalar products as Poisson brackets, and used the invariance of the Poisson brackets and of the phase space volume element on a change of canonical variables. Introducing the Fourier transform of the potential with respect to the angles (103)so that (104)we get (105)Substituting this expression in Eq. (101), we obtain (106)With angleaction variables, the equations of motion of a star determined by the mean field take the very simple form (see Appendix B) (107)where Ω(J, t) is the angular frequency of the orbit with action J. As explained previously, we have neglected the variation of the mean field on a timescale of the order of the dynamical time so it is considered to be frozen when we compute the stellar trajectories (adiabatic or Bogoliubov assumption). Substituting these relations in Eq. (106) and making the transformations l → −l and l_{1} → −l_{1} in the second term (friction term), we obtain successively (108)and (109)It is easy to establish that (110)Therefore, the kinetic equation can be rewritten as (111)Integrating over w_{1}, and recalling that this expression has to be averaged over w, we obtain (112)Making the transformation τ → −τ, then (k, k_{1}) → (−k, −k_{1}), and adding the resulting expression to Eq. (112), we get (113)Finally, using the identity (Eq. (C.11)), we obtain the kinetic equation (114)This kinetic equation was previously derived for systems with arbitrary longrange interactions in various dimensions of space (Chavanis 2007, 2010) and it is here specifically applied to stellar systems. Since collective effects are neglected, this kinetic equation can be viewed as a Landautype equation with angleaction variables describing the evolution of spatially inhomogeneous stellar systems. The collisional evolution of these systems is due to a condition of resonance k·Ω(J, t) = k_{1}·Ω(J_{1}, t) (with (k_{1}, J_{1}) ≠ (k, J)) encapsulated in the δfunction. This δfunction expresses the conservation of energy. It can be shown (Chavanis 2007) that the kinetic equation (Eq. (114)) conserves the mass M = ∫f dJ and the energy E = ∫fϵ(J) dJ, and monotonically increases the Boltzmann entropy S = −∫(f/m)ln (f/m) dJ (Htheorem). However, as explained in Sect. 2.7, this equation does not reach a steady state because of the absence of statistical equilibrium for stellar systems^{29}.
4.3. Relaxation of a star in a thermal bath: a FokkerPlanck equation with angleaction variables
Implementing a test particle approach as in Sect. 3, we find that the equation for P(J, t), the probability density of finding the test star with an action J at time t, is (115)The angular frequency Ω(J) is now a static function determined by the distribution f(J) of the field stars. Equation (115) can be written in the form of a FokkerPlanck equation (116)involving a diffusion tensor (117)and a friction by polarization (118)Writing the FokkerPlanck equation in the usual form (119)with (120)we find that the relation between the friction by polarization and the total friction is (121)Substituting Eqs. (117) and (118) in Eq. (121) and using an integration by parts, we find that the diffusion and friction coefficients are given by These expressions can be obtained directly from the Hamiltonian equations of motion by expanding the trajectories of the stars in powers of 1/N in the limit N → +∞ (Valageas 2006a).
Let us assume that the field stars form a thermal bath with the Boltzmann distribution (124)where ϵ(J) is the energy of a star in an orbit with action J. As we have explained before, this distribution is not defined globally for a selfgravitating system. However, it holds approximately for stars with low energy^{30}. Using the identity ∂ϵ/∂J = Ω(J) (see Appendix B), we find that (125)Substituting this relation in Eq. (118), using the δfunction to replace k_{1}·Ω(J_{1}) by k·Ω(J), and comparing the resulting expression with Eq. (117), we finally get (126)which is the appropriate Einstein relation for our problem. For a thermal bath, using Eq. (126), the FokkerPlanck equation (Eq. (116)) can be written as (127)where D^{μν}(J) is given by Eq. (117) with Eq. (124). Recalling that Ω(J) = ∂ϵ/∂J, this equation is similar to the Kramers equation in Brownian theory (Kramers 1940). This is a driftdiffusion equation describing the evolution of the distribution P(J, t) of the test star in an effective potential U_{eff}(J) = ϵ(J) produced by the field stars. For t → +∞, the distribution of the test star relaxes towards the Boltzmann distribution (Eq. (124)). This takes place on a typical relaxation time t_{R} ~ (N/lnN)t_{D}. Again, this is valid only in the part of the cluster where the Boltzmann distribution holds approximately.
4.4. The impact theory and the wave theory
The approach developed in this paper has removed the problematic divergence at large scales that appears in the VlasovLandau equation when one makes a local approximation. However, this approach, which is based on a weak coupling approximation, yields a divergence at small scales because it does not take strong collisions into account. Stellar systems are therefore described by two complementary kinetic theories depending on the value of the impact parameter λ (see Appendix H). The twobody encounters theory, or impact theory, based on the Boltzmann equation or on the FokkerPlanck equation is appropriate to describe strong collisions when λ ~ λ_{L}. In this theory there is no divergence at small scales. However, this theory does not take spatial inhomogeneity and collective interactions between stars into account. As a result, it exhibits a logarithmic divergence at large scales. This implies that the integrand in the Boltzmann equation is valid only for impact parameters sufficiently smaller than the Jeans length (λ ≪ λ_{J}). The wave theory based on the Landau (or LenardBalescu) equation written with angleaction variables is appropriate to describe spatial inhomogeneity and collective effects when λ ~ λ_{J}. In this theory, there is no divergence at large scales. However, this theory does not take into account the curvature of orbits at small impact parameters. As a result, it exhibits a logarithmic divergence at small scales. This implies that the integrand in the Landau (or LenardBalescu) equation with angleaction variables is valid only for wavelengths sufficiently larger than the Landau length (λ ≫ λ_{L}). For N ≫ 1, the range of validity of these two theories greatly overlaps in the intermediate region λ_{L} ≪ λ ≪ λ_{J} corresponding to the domain of validity of the VlasovLandau equation. This equation describes only weak collisions and presents logarithmic divergences both at small and large scales. Combining the impact theory and the wave theory, we can motivate (but not rigorously justify) using the VlasovLandau equation with a smallscale cutoff at the Landau length and a largescale cutoff at the Jeans length. The success of this equation is that the divergence of the integral over the impact parameter (or wavenumber) is only logarithmic so that its value depends only weakly on the choice of the cutoffs. However, a better kinetic equation is the Landau (or LenardBalescu) equation written with angleaction variables with a smallscale cutoff at the Landau length.
5. Conclusion
Starting from the Liouville equation, using a truncation of the BBGKY hierarchy at the order 1/N, and neglecting strong collisions and collective effects, we have derived a kinetic equation (Eq. (24)) in physical space that can be viewed as a generalized Landau equation. This equation was previously derived by Kandrup (1981) using projection operator technics. A nice feature of this equation is that it does not present any divergence at large scales since the spatial inhomogeneity of the system and its finite extent are properly accounted for. When a local approximation is implemented, and a cutoff is introduced heuristically at the Jeans length, we recover the VlasovLandau equation (Eq. (27)) which is the standard equation of stellar dynamics. On the other hand, using angleaction variables, we have derived a Landautype equation (Eq. (114)) for fully inhomogeneous stellar systems. We have also developed a test particle approach and derived the corresponding FokkerPlanck equations (Eqs. (39) and (115)) in positionvelocity space and angleaction space respectively. Explicit expressions have been given for the diffusion and friction coefficients. We have distinguished the friction by polarization from the total friction. A limitation of the approach presented here is that it neglects collective effects^{31}. More general kinetic equations, corresponding to LenardBalescutype equations taking spatial inhomogeneity and collective effects into account, have been derived recently by Heyvaerts (2010) from the Liouville equation and by Chavanis (2012a) from the Klimontovich equation (these approaches based on the BBGKY hierarchy or on the quasilinear approximation are equivalent but the formalism is different). These kinetic equations are more general than those derived in the present paper, but they are also more complicated (to derive and to solve). Therefore, the equations presented in this paper may be useful as a first step.
We have also discussed the differences between the present approach based on the BBGKY hierarchy and the more classical twobody encounters theory (Chandrasekhar 1942). The twobody encounters theory, which is usually adapted to shortrange potentials (Boltzmann 1872; Chapman & Cowling 1939), can take strong collisions into account so it does not yield any divergence at small scales. However, this approach cannot take spatial inhomogeneity into account so it yields a divergence at large scales. This divergence is due to the longrange nature of the gravitational potential. In addition, the twobody encounters theory does not take collective effects into account; these effects are specific to systems with longrange interactions. By contrast, the approach based on the BBGKY hierarchy takes into account the spatial inhomogeneity of the system and collective effects. Therefore, it does not yield any divergence at large scales. However, it fails to take strong collisions into account because of the weak coupling approximation so it yields a divergence at small scales. In a sense, the gravitational potential is intermediate between shortrange and longrange potentials because both strong collisions and collective effects must be taken into account. Therefore, the approaches adapted to shortrange or longrange potentials are both marginally applicable to selfgravitating systems (they yield a logarithmic divergence at large or small scales respectively). This is why the kinetic equations of stellar dynamics can be obtained in different manners (i.e. from the impact theory or from the wave theory) that turn out to be complementary to each other.
In this paper, we have assumed that the system is isolated from the surrounding. As a result, the source of noise is a result of discreteness (finite N) effects internal to the system. The case where the noise is caused by external sources (perturbations on a galaxy, cosmological environment on dark matter halos) is also interesting. It has been considered by several authors such as Weinberg (2001), Ma & Bertschinger (2004), and Pichon & Aubert (2006) who developed appropriate kinetic theories.
Finally, we would like to stress that the kinetic theories discussed in this paper (leading to complicated equations) do not replace the valuable tools (Nbody simulations and FokkerPlanck codes) used by most astronomers. Nevertheless, they complement these numerical technics and complete our understanding of the gravitational Nbody problem. We also hope that the simplified kinetic equations (Eqs. (114), (115), and (127)) derived in this paper (neglecting collective effects) will be amenable to a numerical analysis.
A few years earlier, Landau (1936) had developed a kinetic theory of Coulombian plasmas taking twobody encounters into account. Starting from the Boltzmann (1872) equation, and making a weak deflection approximation, he derived a kinetic equation for the collisional evolution of neutral plasmas. His approach leads to a divergence at large scales that he removed heuristically by introducing a cutoff at the Debye length λ_{D} (Debye & Hückel 1923) which is the size over which the electric field produced by a charge is screened by the cloud of opposite charges. Later, Lenard (1960) and Balescu (1960) developed a more precise kinetic theory taking collective effects into account. They derived a more elaborate kinetic equation, free of divergence at large scales, in which the Debye length appears naturally. This justifies the heuristic procedure of Landau.
The stochastic evolution of a star is primarily due to many small deflections produced by relatively distant encounters. The problem of treating particles undergoing numerous weak deflections was originally encountered in relation to the Brownian motion of large molecules which are thermally agitated by the smaller field molecules. The stochastic character of the many small impulses which act on a suspended particle is usually described by the FokkerPlanck equation. Chandrasekhar (1943c) noted the analogy between stellar dynamics and Brownian theory and employed a FokkerPlanck equation to describe the evolution of the velocity distribution function of stars. Indeed, he argued that a given star undergoes many smallangle (large impact parameter) collisions in a time that is small compared with that in which its position or velocity changes appreciably. Since the Newtonian potential makes the cumulative effect of these small momentum transfers dominant, the stochastic methods of the FokkerPlanck equation should be more appropriate for a stellar system than the Boltzmann approach.
There remains a logarithmic divergence at small scales as in the original Landau (1936) theory due to the neglect of strong collisions. This divergence can be cured heuristically by introducing a cutoff at the Landau length λ_{L} corresponding to the impact parameter leading to a deflection at 90° (Landau 1936). This divergence does not occur in the theory of Chandrasekhar (1942, 1943a,b) which takes strong collisions into account.
In a plasma, since the Coulomb force between electrons is repulsive, each particle of the plasma tends to attract to it particles of opposite charge and to repel particles of like charge, thereby creating a kind of cloud of opposite charge which screens the interaction at the scale of the Debye length. In the gravitational case, since the Newton force is attractive, the situation is considerably different. The test star draws neighboring stars into its vicinity and these add their gravitational force to that of the test star itself. The bare gravitational force of the test star is thus augmented rather than shielded. The polarization acts to increase the effective gravitational mass of a test star.
For a review on the dynamics and thermodynamics of systems with longrange interactions, see Campa et al. (2009).
Actually, the interaction term becomes large at small scales r_{2} − r_{1} → 0, so its effect is not totally negligible (in other words, the expansion in terms of 1/N is not uniformly convergent). In particular, the interaction term ℒ′ must be taken into account in order to describe strong collisions with small impact parameter that lead to large deflections very different from the mean field trajectory corresponding to ℒ_{0} + ℒ_{m.f.}. This is important to regularize the divergence at small scales that appears in the Landau equation. We shall return to this problem in Sect. 2.4.
The term corresponding to g(x, x_{1}, t = 0) in the kinetic equation has been called the “destruction term” by Prigogine and Résibois because it describes the destruction of the effect of the initial correlations (Balescu 2000).
It is sometimes argued that the Markovian approximation is not justified for stellar systems because the force autocorrelation decreases slowly as 1/t (Chandrasekhar 1944). However, this result is only true for an infinite homogeneous system (see Appendix G). For spatially inhomogeneous distributions, the correlation function decreases more rapidly and the Markovian approximation is justified (Severne & Haggerty 1976).
We note that Eq. (25) with G_{0} replaced by the total Green function G′ taking into account the interaction term is equivalent to the Boltzmann equation. Indeed, for spatially homogeneous systems, the Boltzmann equation can be derived from the BBGKY hierarchy (Eqs. (10) and (11)) by keeping the Liouville operator ℒ = ℒ_{0} + ℒ′ describing the twobody problem exactly and the source (Balescu 2000). Therefore, the procedure used by Landau which amounts to expanding the Boltzmann equation in the limit of weak deflections is equivalent to the one presented here that starts from the BBGKY hierarchy and neglects ℒ′ in the Liouville operator.
This distribution is not steady in the sense that the coefficients A, β, and r_{a} slowly vary in time as the stars escape and the system loses mass and energy. However, the distribution keeps the same form during the evaporation process (King 1965). Hénon (1961) previously found a selfsimilar solution of the orbitaveraged FokkerPlanck equation. He showed that the core contracts as the cluster evaporates and that the central density is infinite (the structure of the core resembles the singular isothermal sphere with a density ρ ∝ r^{2}). He argued that the concentration of energy, without concentration of mass, at the center of the system is attributable to the formation of tight binary stars. The invariant profile found by Hénon does not exactly coincide with the MichieKing distribution, but is reasonably close.
The collisional evolution of the system can be measured precisely in terms of the scale escape energy . The onset of instability corresponds to x_{0} ≃ 9.3 which is the value at which the King model becomes thermodynamically unstable (Cohn 1980).
Actually, the number of stars in the core decreases as the collapse proceeds. When <10^{2} stars remain in the core, strong collisions become important. Goodman (1983) showed that the evolution caused by strong collisions is very similar to that caused by numerous weak ones, except for the rate of evaporation which is greatly increased.
In plasma physics, this steady distribution is the Boltzmann distribution of statistical equilibrium which is the steady state of the Landau equation. In the case of stellar systems, there is an intrinsic difficulty since no statistical equilibrium state exists in a strict sense: the Boltzmann distribution has infinite mass, and the VlasovLandau equation has no steady state because of the escape of highenergy stars. However, we have seen that the system can reach a quasisteady distribution (e.g. a MichieKing distribution) and that this distribution changes on an evaporation timescale that is long with respect to the collisional relaxation time. Therefore, we can consider that this distribution is steady on the collisional timescale over which the FokkerPlanck approach applies.
The collisional evolution of a star, under the effect of twobody encounters, can be understood as an interplay between two competing effects. The fluctuations of the gravitational field induce a diffusion in velocity space which tends to increase the speed of the star. This effect is counterbalanced by a friction (dissipation) which results in a systematic deceleration along the direction of motion. As emphasized by Chandrasekhar (1943a, 1949) in his Brownian theory of stellar motion, the Einstein relation (Eq. (53)) guarantees that the Maxwell distribution (Eq. (51)) is a steady state of the FokkerPlanck equation (Eq. (36)).
It is important to emphasize, however, that the twobody encounters theory of Chandrasekhar (1942) can take strong collisions into account. Therefore, there is no divergence at small impact parameters and the gravitational Landau length λ_{L} appears naturally in the Coulombian logarithm lnΛ, contrary to the Landau theory which ignores strong collisions (see Appendix H).
It is moreover applicable with almost no modification to other systems with longrange interactions for which a twobody encounters theory is not justified (Chavanis 2010).
The adiabatic assumption (Eq. (100)) is consistent with the Bogoliubov ansatz of kinetic theory. Since the correlation time of the fluctuations is of the order of the dynamical time or shorter, we can freeze the distribution function at time t to compute the integral over τ (see Sect. 2.4). This distribution function, which is a steady state of the Vlasov equation, defines a set of angleaction variables w and J that we can use to perform the integral over τ. Then, the distribution function f(J, t) evolves with time on a longer timescale according to Eq. (101).
For selfgravitating systems in lower dimensions of space, or for systems with longrange interactions with a smooth potential like the HMF model, a statistical equilibrium state exists. In this case, it can be shown that the Landautype equation (Eq. (114)) relaxes towards the Boltzmann distribution on a timescale Nt_{D} provided there are enough resonances (see Chavanis 2007).
As discussed in Appendix E, collective effects can substantially decrease the relaxation time for systems at the verge of instability.
If we account for shortrange interactions (quantum mechanics, soften potential, hard spheres) there is a new dimensionless parameter in the problem that we shall generically call μ. In that case, the proper thermodynamic limit of selfgravitating systems with shortrange interactions corresponds to N → +∞ in such a way that ϵ, η, and μ are of order unity. Some examples of scalings are given by Chavanis & Rieutord (2003).
The approach of Chandrasekhar takes into account strong collisions with an impact parameter smaller than the Landau length λ_{L} that yield a deflection at an angle larger than 90°. As we have seen, this can suppress the divergence at small scales. However, the Chandrasekhar approach does not take into account the possibility of forming binary stars which correspond to bound states with strong correlation between particles. In other words, Chandrasekhar only considers hyperbolic trajectories and not elliptical ones in the twobody problem. The effect of binary stars in the kinetic theory requires a special treatment (Heggie 1975).
Actually, Landau studies the evolution of the system as a whole while Chandrasekhar studies the relaxation of a test star in a bath of field stars. Landau obtains an integrodifferential equation (the Landau equation) that conserves the energy (microcanonical description) while Chandrasekhar obtains a differential equation (the FokkerPlanck or Kramers equation) that does not conserve the energy and involves a fixed temperature (canonical description). However, these two equations are clearly related as explained in Sect. 3.
References
 Antonov, V. A. 1962, Vest. Leningr. Gos. Univ., 7, 135 [Google Scholar]
 Ambartsumian, V. A. 1938, Ann. Leningrad State Univ., 22, 19 [Google Scholar]
 Arad, I., & Johansson, P. H. 2005, MNRAS, 362, 252 [NASA ADS] [CrossRef] [Google Scholar]
 Balescu, R. 1960, Phys. Fluids, 3, 52 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Balescu, R. 2000, Statistical Dynamics: Matter out of Equilibrium (Imperial College Press) [Google Scholar]
 Bekenstein, J., & Maoz, E. 1992, ApJ, 390, 79 [NASA ADS] [CrossRef] [Google Scholar]
 Bertin, G., & Stiavelli, M. 1984, A&A, 137, 26 [NASA ADS] [Google Scholar]
 Bettwieser, E., & Sugimoto, D. 1984, MNRAS, 208, 493 [NASA ADS] [Google Scholar]
 Binney, J., & Lacey, C. 1988, MNRAS, 230, 597 [NASA ADS] [CrossRef] [Google Scholar]
 Binney, J., & Tremaine, S. 2008, Galactic Dynamics (Princeton Series in Astrophysics) [Google Scholar]
 Boltzmann, L. 1872, Wien, Ber., 66, 275 [Google Scholar]
 Campa, A., Dauxois, T., & Ruffo, S. 2009, Phys. Rep., 480, 57 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Chabanol, M. L., Corson, F., & Pomeau, Y. 2000, Europhys. Lett., 50, 148 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Chandrasekhar, S. 1942, Principles of Stellar Dynamics (University of Chicago press) [Google Scholar]
 Chandrasekhar, S. 1943a, ApJ, 97, 255 [NASA ADS] [CrossRef] [Google Scholar]
 Chandrasekhar, S. 1943b, ApJ, 97, 263 [NASA ADS] [CrossRef] [Google Scholar]
 Chandrasekhar, S. 1943c, Rev. Mod. Phys., 15, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Chandrasekhar, S. 1944, ApJ, 99, 47 [NASA ADS] [CrossRef] [Google Scholar]
 Chandrasekhar, S. 1949, Rev. Mod. Phys., 21, 383 [NASA ADS] [CrossRef] [Google Scholar]
 Chandrasekhar, S., & von Neumann, J. 1942, ApJ, 95, 489 [Google Scholar]
 Chapman, S., & Cowling, T. G. 1939, The Mathematical Theory of Nonuniform Gases (Cambridge Mathematical Library) [Google Scholar]
 Chavanis, P. H. 1998, MNRAS, 300, 981 [NASA ADS] [CrossRef] [Google Scholar]
 Chavanis, P. H. 2005, A&A, 432, 117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Chavanis, P. H. 2006, Int. J. Mod. Phys. B, 20, 3113 [Google Scholar]
 Chavanis, P. H. 2007, Physica A, 377, 469 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Chavanis, P. H. 2008a, Physica A, 387, 787 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Chavanis, P. H. 2008b, Physica A, 387, 1504 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Chavanis, P. H. 2010, J. Stat. Mech., 5, 05019 [CrossRef] [Google Scholar]
 Chavanis, P. H. 2012a, Physica A, 391, 3680 [Google Scholar]
 Chavanis, P. H. 2012b, EPJ Plus, 127, 19 [Google Scholar]
 Chavanis, P. H., & Rieutord, M. 2003, A&A, 412, 1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Chavanis, P. H., Sommeria, J., & Robert, R. 1996, ApJ, 471, 385 [NASA ADS] [CrossRef] [Google Scholar]
 Cohen, R. S., Spitzer, L., & Routly, P. M. 1950, Phys. Rev., 80, 230 [NASA ADS] [CrossRef] [Google Scholar]
 Cohn, H. 1980, ApJ, 242, 765 [NASA ADS] [CrossRef] [Google Scholar]
 Debye, P., & Hückel, E. 1923, Phys. Z., 24, 305 [Google Scholar]
 de Vega, H. J., & Sanchez, N. 2002, Nucl. Phys. B, 625, 409 [NASA ADS] [CrossRef] [Google Scholar]
 Gasiorowicz, S., Neuman, M., & Riddell, J. R. 1956, Phys. Rev., 101, 922 [NASA ADS] [CrossRef] [Google Scholar]
 Gilbert, I. H. 1968, ApJ, 152, 1043 [NASA ADS] [CrossRef] [Google Scholar]
 Gilbert, I. H. 1970, ApJ, 159, 239 [Google Scholar]
 Goldstein, H. 1956, Classical Mechanics (Cambridge MA: AddisonWesley) [Google Scholar]
 Goodman, J. 1983, ApJ, 270, 700 [NASA ADS] [CrossRef] [Google Scholar]
 Goodman, J. 1984, ApJ, 280, 298 [NASA ADS] [CrossRef] [Google Scholar]
 Goodman, J. 1988, ApJ, 329, 612 [NASA ADS] [CrossRef] [Google Scholar]
 Gurzadyan, V. G., & Savvidy, G. K. 1986, A&A, 160, 203 [NASA ADS] [Google Scholar]
 Heggie, D. C. 1975, MNRAS, 173, 729 [NASA ADS] [CrossRef] [Google Scholar]
 Heggie, D., & Hut, P. 2003, The Gravitational MillionBody Problem: A Multidisciplinary Approach to Star Cluster Dynamics (Cambridge University Press) [Google Scholar]
 Hénon, M. 1961, Ann. Astrophys., 24, 369 [Google Scholar]
 Hénon, M. 1964, Ann. Astrophys., 27, 83 [Google Scholar]
 Hénon, M. 1982, A&A, 114, 211 [NASA ADS] [Google Scholar]
 Heyvaerts, J. 2010, MNRAS, 407, 355 [NASA ADS] [CrossRef] [Google Scholar]
 Hjorth, J., & Madsen, J. 1991, MNRAS, 253, 703 [NASA ADS] [Google Scholar]
 Hubbard, J. 1961a, Proc. R. Soc. Lond., 260, 114 [Google Scholar]
 Hubbard, J. 1961b, Proc. R. Soc. Lond., 261, 371 [Google Scholar]
 Inagaki, S., & LyndenBell, D. 1983, MNRAS, 205, 913 [NASA ADS] [Google Scholar]
 Jeans, J. H. 1915, MNRAS, 76, 70 [NASA ADS] [CrossRef] [Google Scholar]
 Jeans, J. H. 1929, Astronomy and Cosmogony (Cambridge Univ. Press) [Google Scholar]
 Joyce, M., Marcos, B., & SylosLabini, F. 2009, PRE, 397, 775 [Google Scholar]
 Kac, M., Uhlenbeck, G. E., & Hemmer, P. C. 1963, J. Math. Phys., 4, 216 [NASA ADS] [CrossRef] [Google Scholar]
 Kadomtsev, B. B., & Pogutse, O. P. 1970, PRL, 25, 1155 [NASA ADS] [CrossRef] [Google Scholar]
 Kalnajs, A. 1971a, ApJ, 166, 275 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Kalnajs, A. 1971b, Astro. Space Sci., 13, 279 [Google Scholar]
 Kalnajs, A. 1977, ApJ, 212, 637 [NASA ADS] [CrossRef] [Google Scholar]
 Kandrup, H. 1981, ApJ, 244, 316 [NASA ADS] [CrossRef] [Google Scholar]
 Kandrup, H. 1983, ApSS, 97, 435 [Google Scholar]
 Katz, J. 2003, Found. Phys., 33, 223 [CrossRef] [MathSciNet] [Google Scholar]
 King, I. R. 1960, AJ, 65, 122 [Google Scholar]
 King, I. R. 1965, AJ, 70, 376 [NASA ADS] [CrossRef] [Google Scholar]
 King, I. R. 1966, AJ, 71, 64 [Google Scholar]
 Kramers, H. A. 1940, Physica, 7, 284 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Landau, L. D. 1936, Phys. Z. Sowj. Union, 10, 154 [Google Scholar]
 Larson, R. B. 1970, MNRAS, 147, 323 [NASA ADS] [Google Scholar]
 Lee, E. P. 1968, ApJ, 151, 687 [NASA ADS] [CrossRef] [Google Scholar]
 Lemou, M., & Chavanis, P. H. 2010, Physica A, 389, 1021 [NASA ADS] [CrossRef] [Google Scholar]
 Lenard, A. 1960, Ann. Phys. (N. Y.), 10, 390 [NASA ADS] [CrossRef] [Google Scholar]
 Lerche, I. 1971, ApJ, 166, 207 [NASA ADS] [CrossRef] [Google Scholar]
 Levin, Y., Pakter, R., & Rizzato, F. 2008, PRE, 78, 021130 [Google Scholar]
 Londrillo, P., Messina, A., & Stiavelli, M. 1991, MNRAS, 250, 54 [NASA ADS] [Google Scholar]
 LyndenBell, D. 1967, MNRAS, 136, 101 [NASA ADS] [CrossRef] [Google Scholar]
 LyndenBell, D., & Eggleton, P. P. 1980, MNRAS, 191, 483 [NASA ADS] [Google Scholar]
 LyndenBell, D., & Kalnajs, A. J. 1972, MNRAS, 157, 1 [NASA ADS] [CrossRef] [Google Scholar]
 LyndenBell, D., & Wood, R. 1968, MNRAS, 138, 495 [NASA ADS] [Google Scholar]
 Ma, C. P., & Bertschinger, E. 2004, ApJ, 612, 28 [NASA ADS] [CrossRef] [Google Scholar]
 Maoz, E. 1993, MNRAS, 263, 75 [NASA ADS] [Google Scholar]
 Marochnik, L. S. 1968, Sov. Astr., 11, 873 [Google Scholar]
 Messer, J., & Spohn, H. 1982, J. Stat. Phys., 29, 561 [NASA ADS] [CrossRef] [Google Scholar]
 Michie, R. 1963, MNRAS, 125, 127 [NASA ADS] [CrossRef] [Google Scholar]
 Miller, R. H. 1966, ApJ, 146, 831 [NASA ADS] [CrossRef] [Google Scholar]
 Monaghan, J. J. 1978, MNRAS, 184, 25 [NASA ADS] [Google Scholar]
 Nelson, R. W., & Tremaine, S. 1999, MNRAS, 306, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Padmanabhan, T. 1990, Phys. Rep., 188, 285 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Parisot, J. P., & Severne, G. 1979, Astro. Space Sci., 61, 121 [NASA ADS] [CrossRef] [Google Scholar]
 Pichon, C., & Aubert, D. 2006, MNRAS, 368, 1657 [NASA ADS] [CrossRef] [Google Scholar]
 Pichon, C., & Cannon, R. C. 1997, MNRAS, 291, 616 [NASA ADS] [CrossRef] [Google Scholar]
 Risken, H. 1989, The FokkerPlanck equation (Springer) [Google Scholar]
 Rosenbluth, M., MacDonald, W., & Judd, D. 1957, Phys. Rev., 107, 1 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Rosseland, S. 1928, MNRAS, 88, 208 [NASA ADS] [Google Scholar]
 Roy, F., & Perez, J. 2004, MNRAS, 348, 62 [NASA ADS] [CrossRef] [Google Scholar]
 Schwarzschild, K. 1924, Seeliger Festschrift, 94 [Google Scholar]
 Severne, G., & Haggerty, M. J. 1976, Astro. Space Sci., 45, 287 [Google Scholar]
 Severne, G., & Luwel, M. 1980, Astro. Space Sci., 72, 293 [Google Scholar]
 Smart, W. M. 1938, Stellar Dynamics (Cambridge Univ. Press) [Google Scholar]
 Spitzer, L. 1940, MNRAS, 100, 396 [NASA ADS] [CrossRef] [Google Scholar]
 Spitzer, L. 1987, Dynamical Evolution of Globular Clusters (Princeton Series in Astrophysics) [Google Scholar]
 Spitzer, L., & Härm, R. 1958, ApJ, 127, 544 [NASA ADS] [CrossRef] [Google Scholar]
 Stiavelli, M., & Bertin, G. 1987, MNRAS, 229, 61 [NASA ADS] [Google Scholar]
 Thorne, R. M. 1968, ApJ, 151, 671 [NASA ADS] [CrossRef] [Google Scholar]
 Tremaine, S., & Weinberg, M. 1984, MNRAS, 209, 729 [NASA ADS] [CrossRef] [Google Scholar]
 Tremaine, S., Hénon, M., & LyndenBell, D. 1986, MNRAS, 219, 285 [NASA ADS] [CrossRef] [Google Scholar]
 Trenti, M., Bertin, G., & van Albada, T. S. 2005, A&A, 433, 57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Valageas, P. 2006a, A&A, 450, 445 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Valageas, P. 2006b, PRE, 74, 6606 [Google Scholar]
 van Albada, T. S. 1982, MNRAS, 201, 939 [NASA ADS] [CrossRef] [Google Scholar]
 Vlasov, A. A. 1938, Zh. Eksp. i Teor. Fiz., 8, 291 [Google Scholar]
 Weinberg, M. D. 1991, ApJ, 368, 66 [NASA ADS] [CrossRef] [Google Scholar]
 Weinberg, M. D. 1993, ApJ, 410, 543 [Google Scholar]
 Weinberg, M. D. 1998, MNRAS, 297, 101 [NASA ADS] [CrossRef] [Google Scholar]
 Weinberg, M. D. 2001, MNRAS, 328, 311 [NASA ADS] [CrossRef] [Google Scholar]
Online material
Appendix A: The thermodynamic limit and the weak coupling approximation
In this Appendix, we discuss the proper thermodynamic limit of selfgravitating systems and justify the weak coupling approximation.
We call v_{m} the root mean square velocity of the stars and R the system’s size. The kinetic energy and the potential energy U ~ N^{2}Gm^{2}/R in the Hamiltonian (Eq. (1)) are comparable provided that (this fundamental scaling may also be obtained from the virial theorem). As a result, the energy scales as and the kinetic temperature, defined by , scales as k_{B}T ~ GNm^{2}/R ~ E/N. Inversely, these relations may be used to define R and v_{m} as a function of the energy E (conserved quantity in the microcanonical ensemble) or as a function of the temperature T (fixed quantity in the canonical ensemble).
The proper thermodynamic limit of a selfgravitating system corresponds to N → +∞ in such a way that the normalized energy ϵ = ER/(GN^{2}m^{2}) and the normalized temperature η = βGNm^{2}/R are of order unity. Of course, the usual thermodynamic limit N,V → +∞ with N/V ~ 1 is not applicable to selfgravitating systems since these systems are spatially inhomogeneous.
From R and v_{m}, we define the dynamical time . If we rescale the distances by R, the velocities by v_{m}, and the times by t_{D}, we find that the equations of the BBGKY hierarchy depend on a single dimensionless parameter η/N. This is the equivalent of the plasma parameter in plasma physics. It is usually argued that the correlation functions scale as . Therefore, when N → +∞, we can expand the equations of the BBGKY hierarchy in powers of 1/N ≪ 1. This corresponds to a weak coupling approximation (see below). We note that the thermodynamic limit is welldefined for the outofequilibrium problem although there is no statistical equilibrium state (i.e. the density of state and the partition function diverge).
By a suitable normalization of the parameters, we can take R ~ v_{m} ~ t_{D} ~ m ~ 1. In that case, we must impose G ~ 1/N. This is the Kac prescription (Kac et al. 1963; Messer & Spohn 1982). With this normalization, E ~ N, S ~ N and T ~ 1. The energy and the entropy are extensive but they remain fundamentally nonadditive (Campa et al. 2009). The temperature is intensive. This normalization is very convenient since the length, velocity and time scales are of order unity. Furthermore, since the coupling constant G scales as 1/N, this immediately shows that a regime of weak coupling holds when N ≫ 1.
Other normalizations of the parameters are possible. For example, Gilbert (1968) considers the limit N → +∞ with R ~ v_{m} ~ t_{D} ~ G ~ 1 and m ~ 1/N. In that case, E ~ 1, S ~ N, and T ~ 1/N. On the other hand, de Vega & Sanchez (2002) define the thermodynamic limit as N → +∞ with m ~ G ~ v_{m} ~ 1 and R ~ N in order to have E ~ N, S ~ N, and T ~ 1. This normalization is natural since m and G should not depend on N. However, in that case, the dynamical time t_{D} = R/v_{m} ~ N diverges with the number of particles. Therefore, this normalization is not very convenient to develop a kinetic theory of stellar systems^{32}. If we impose G ~ 1, E ~ N, T ~ 1 and t_{D} ~ 1, we get R ~ N^{1/5} and m ~ N^{−2/5} (this corresponds to ρ ~ 1 since ). We stress that all these scalings are equally valid. The important thing is that ϵ and η are O(1) when N → +∞. One should choose the most convenient scaling which, in our opinion, is the Kac scaling^{33}.
We can also present the preceding results in the following manner, by analogy with plasma physics. A fundamental length scale in selfgravitating systems is the Jeans length λ_{J} = (4πGβmρ)^{−1/2}. This is the counterpart of the Debye length λ_{D} = (4πe^{2}βρ/m)^{−1/2} in plasma physics. Since , we find that λ_{J} ~ R. Therefore, the Jeans length is of the order of the system’s size. On the other hand, a fundamental time scale in selfgravitating systems is provided by the gravitational pulsation ω_{G} = (4πGρ)^{1/2}. This is the counterpart of the plasma pulsation ω_{P} = (4πe^{2}ρ/m^{2})^{1/2} in plasma physics. We can define the dynamical time by . If we rescale the distances by λ_{J} and the times by t_{D}, we find that the equations of the BBGKY hierarchy depend on a single dimensionless parameter , where gives the number of stars in the Jeans sphere. This is the counterpart of the number of electrons in the Debye sphere in plasma physics. When Λ ≫ 1, we can expand the equations of the BBGKY hierarchy in terms of the small parameter 1/Λ.
Finally, we show that 1/Λ may be interpreted as a coupling parameter. The coupling parameter Γ is defined as the ratio of the interaction strength at the mean interparticle distance Gm^{2}n^{1/3} (resp. e) to the thermal energy k_{B}T. This leads to (resp. ). If we define the coupling parameter g as the ratio of the interaction strength at the Jeans (resp. Debye) length Gm^{2}/λ_{J} (resp. e^{2}/λ_{D}) to the thermal energy k_{B}T, we get g = 1/Λ (or g = E_{pot}/E_{kin} = (Gm^{2}/R)/k_{B}T = η/N with the initial variables). Therefore, the expansion of the BBGKY hierarchy in terms of the coupling parameter Γ or g is equivalent to an expansion in terms of the inverse of the number of particles in the Jeans sphere (resp. Debye sphere ). The weak coupling approximation is therefore justified when Λ ≫ 1.
Appendix B: Angleaction variables
In Sect. 4, we have explained that during its collisional evolution a stellar system passes by a succession of QSSs that are steady states of the Vlasov equation slowly changing under the effect of close encounters (finite N effects). The slowly varying distribution function f(r, v) determines a potential Φ(r) and a oneparticle Hamiltonian ϵ = v^{2}/2 + Φ(r) that we assume to be integrable. Therefore, it is possible to use angleaction variables constructed with this Hamiltonian (Goldstein 1956; Binney & Tremaine 2008). This construction is done adiabatically, i.e. the distribution function, and the angleaction variables, slowly change in time.
A particle with coordinates (r, v) in phase space is described equivalently by the angleaction variables (w, J). The Hamiltonian equations for the conjugate variables (r, v) are (B.1)In terms of the variables (r, v), the dynamics is complicated because the potential explicitly appears in the second equation. Therefore, this equation dv/dt = −∇Φ cannot be easily integrated except if Φ = 0, i.e. for a spatially homogeneous system. In that case, the velocity v is constant and the unperturbed equations of motion reduce to r = vt + r_{0}, i.e. to a rectilinear motion at constant velocity. Now, the angleaction variables are constructed so that the Hamiltonian does not depend on the angles w. Therefore, the Hamiltonian equations for the conjugate variables (w, J) are (B.2)where Ω(J) is the angular frequency of the orbit with action J. From these equations, we find that J is a constant and that w = Ω(J)t + w_{0}. Therefore, the equations of motion are very simple in these variables. They extend naturally the trajectories at constant velocity for spatially homogeneous systems. This is why this choice of variables is relevant to develop the kinetic theory. Of course, even if the description of the motion becomes simple in these variables, the complexity of the problem has not completely disappeared. It is now embodied in the relation between the position and momentum variables and the angle and action variables which can be very complicated.
Appendix C: Calculation of K^{μν}
In this Appendix, we compute the tensor K^{μν} that appears in the VlasovLandau equation (Eq. (28)). Within the local approximation, we can proceed as if the system were spatially homogeneous. In that case, the mean field force vanishes, ⟨ F ⟩ = 0, and the unperturbed equations of motion (i.e. for N → +∞) reduce to corresponding to a rectilinear motion at constant velocity. The collision term in the kinetic equation (Eq. (26)) can be written as (C.3)with (C.4)The force by unit of mass created by particle 1 on particle 0 is given by (C.5)where u(r − r′) = −G/r − r′ is the gravitational potential. The Fourier transform and the inverse Fourier transform of the potential are defined by (C.6)For the gravitational interaction: (C.7)Substituting Eq. (C.6b) in Eq. (C.5), and writing explicitly the Lagrangian coordinates, we get (C.8)Using the equations of motion (C.1) and (C.2), and introducing the notations x = r − r_{1} and w = v − v_{1}, we obtain (C.9)Therefore, (C.10)Using the identity (C.11)and integrating over x and k′, we find that (C.12)Performing the transformation τ → −τ, then k → −k, and adding the resulting expression to Eq. (C.12), we get (C.13)Using the identity (C.11), we finally obtain (C.14)which leads to Eq. (27).
Introducing a spherical system of coordinates in which the z axis is taken in the direction of w, we find that (C.15)Using k_{x} = ksinθcosφ, k_{y} = ksinθsinφ and k_{z} = kcosθ, it is easy to see that only K_{xx}, K_{yy} and K_{zz} can be nonzero. The other components of the matrix K^{μν} vanish by symmetry. Furthermore, (C.16)Using the identity , we get (C.17)With the change of variables s = cosθ, we obtain (C.18)so that, finally, (C.19)On the other hand, (C.20)In conclusion, we obtain (C.21)with (C.22)Using Eq. (C.7), this leads to Eqs. (28)–(30).
Appendix D: Another derivation of the Landau equation
For an infinite homogeneous system, the distribution function and the twobody correlation function can be written as f(0) = f(v, t) and g(0,1) = g(r − r_{1},v, v_{1}, t). In that case, Eqs. (18) and (19) become where we have defined x = r − r_{1} and w = v − v_{1}. Using the Bogoliubov ansatz, we shall treat the distribution function f as a constant, and determine the asymptotic value g(x, v, v_{1}, +∞) of the correlation function. Introducing the Fourier transforms of the potential of interaction and of the correlation function, Eq. (D.1) may be replaced by (D.3)where we have used the reality condition ĝ(−k) = ĝ(k)^{ ∗ }. On the other hand, taking the LaplaceFourier transform of Eq. (D.2) and assuming that no correlation is present initially (if there are initial correlations, their effect becomes rapidly negligible), we get (D.4)Taking the inverse Laplace transform of Eq. (D.4), and using the residue theorem, we find that the asymptotic value t → +∞ of the correlation function, determined by the pole ω = 0, is (D.5)Using the Plemelj formula (D.6)we get (D.7)Substituting Eq. (D.7) in Eq. (D.3), we obtain the Landau equation (Eq. (27)).
Appendix E: LenardBalescu equation for homogeneous stellar systems
If we assume that the system is spatially homogeneous (or make the local approximation), and take collective effects into account, the VlasovLandau Eq. (27) is replaced by the VlasovLenardBalescu equation (E.1)where ϵ(k, ω) is the dielectric function (E.2)The Landau equation is recovered by taking ϵ(k, k·v)^{2} = 1. The LenardBalescu equation generalizes the Landau equation by replacing the bare potential of interaction û(k) by a dressed potential of interaction (E.3)The dielectric function in the denominator takes the dressing of the particles by their polarization cloud into account. In plasma physics, this term corresponds to a screening of the interactions. The LenardBalescu equation accounts for dynamical screening since the velocity v of the particles explicitly appears in the effective potential. However, for Coulombian interactions, it is a good approximation to neglect the deformation of the polarization cloud due to the motion of the particles and use the static results on screening (Debye & Hückel 1923). This amounts to replacing the dynamic dielectric function ϵ(k, k·v) by the static dielectric function ϵ(k, 0). In this approximation, û_{dressed}(k, k·v) is replaced by the DebyeHückel potential corresponding to u_{DH}(x) = (e^{2}/m^{2})e^{−kDr}/r in physical space. If we make the same approximation for stellar systems, we find that û_{dressed}(k, k·v) is replaced by (E.4)corresponding to u_{DH}(x) = −Gcos(k_{J}r)/r in physical space. In this approximation, the VlasovLenardBalescu equation (Eq. (E.1)) takes the same form as the VlasovLandau equation (Eq. (28)) except that A is now given by A = 2πmG^{2}Q with (E.5)where R is the system’s size. We see that Q diverges algebraically, as (λ_{J} − R)^{1}, when R → λ_{J} instead of yielding a finite Coulombian logarithm lnΛ ~ lnN when collective effects are neglected^{34}. This naive approach shows that collective effects (which account for antishielding) tend to increase the diffusion coefficient and consequently tend to reduce the relaxation time. This is the conclusion reached by Weinberg (1993) with a more precise approach. However, his approach is not fully satisfactory since the system is assumed to be spatially homogeneous and the ordinary LenardBalescu equation is used. In that case, the divergence when R → λ_{J} is a manifestation of the Jeans instability that a spatially homogeneous selfgravitating system experiences when the size of the perturbation overcomes the Jeans length. For inhomogeneous systems, the Jeans instability is suppressed so the results of Weinberg should be used with caution. Heyvaerts (2010) and Chavanis (2012a) derived a more satisfactory LenardBalescu equation that is valid for spatially inhomogeneous stable selfgravitating systems. This equation does not present any divergence at large scales. However, this equation is complicated and it is difficult to measure the importance of collective effects. The approach of Weinberg (1993), and the arguments given in this Appendix, suggest that collective effects reduce the relaxation time of inhomogeneous stellar systems. This reduction should be particularly strong for a system close to instability because of the enhancement of fluctuations (Monaghan 1978).
These considerations show that it is not possible to make a local approximation and simultaneously take collective effects into account. The usual procedure is to ignore collective effects, make a local approximation, and introduce a largescale cutoff at the Jeans length. This is the procedure that is usually followed in stellar dynamics (Binney & Tremaine 2008). The only rigorous manner to take collective effects into account is to use the LenardBalescu equation written with angleaction variables.
Appendix F: Multispecies systems
It is straightforward to generalize the kinetic theory of stellar systems for several species of stars. The VlasovLandau equation (Eq. (27)) is replaced by (F.1)where f^{a}(r, v, t) is the distribution function of species a normalized such that , and the sum ∑ _{b} runs over all species. We can use this equation to give a new interpretation of the test particle approach developed in Sect. 3. We make three assumptions: (i) We assume that the system is composed of two types of stars, the test stars with mass m and the field stars with mass m_{f}; (ii) we assume that the number of test stars is much lower than the number of field stars; and (iii) we assume that the field stars are in a steady distribution f(r, v). Because of assumption (ii), the collisions between the field stars and the test stars do not alter the distribution of the field stars so that the field stars remain in their steady state. The collisions of the test stars among themselves are also negligible, so they only evolve as a result of collisions with the field stars. Therefore, if we call P(r, v, t) the distribution function of the test stars (to have notations similar to those of Sect. 3 with, however, a different interpretation), its evolution is given by the FokkerPlanck equation obtained from Eq. (F.1) yielding (F.2)The diffusion and friction coefficients are given by We recall that the diffusion coefficient is due to the fluctuations of the gravitational force produced by the field stars, while the friction by polarization is due to the perturbation on the distribution of the field stars caused by the test stars. This explains the occurrence of the masses m_{f} and m in Eqs. (F.3) and (F.4), respectively.
Using Eq. (46) and noting that (F.5)we get (F.6)If we assume furthermore that m ≫ m_{f}, we find that F_{friction} ≃ F_{pol}. However, in general, the friction force is different from the friction by polarization. The other results of Sect. 3 can be easily generalized to multispecies systems.
If the field stars have an isothermal distribution, then from Eq. (F.4), and the FokkerPlanck equation (Eq. (F.2)) reduces to (F.7)where D^{μν}(v) is given by Eq. (F.3). Using the results of Sect. 3, we get (F.8)where is the velocity dispersion of the field stars and ρ_{f} = n_{f}m_{f} their density. On the other hand, lnΛ = ln (λ_{J}/λ_{L}), where is the Jeans length and is the Landau length. This yields . The total friction is . If the distribution of the field stars is f ∝ e^{−βmfv2/2}, the equilibrium distribution of the test stars is P ∝ e^{−βmv2/2} ∝ f^{m/mf}. When m ≫ m_{f}, the evolution is dominated by frictional effects; when m ≪ m_{f} it is dominated by diffusion. The relaxation time scales like where is the r.m.s. velocity of the test stars. Therefore, , where ξ is the friction coefficient associated to the friction by polarization. The true friction coefficient is ξ_{ ∗ } = (1 + m_{f}/m)ξ. We get , where we have used .
Appendix G: Temporal correlation tensor of the gravitational force
The diffusion of the stars is caused by the fluctuations of the gravitational force. For an infinite homogeneous system (or in the local approximation), the diffusion tensor can be derived from Eq. (43) that is wellknown in Brownian theory. This expression involves the temporal autocorrelation tensor of the gravitational force experienced by a star. It can be written as (G.1)We first compute the tensor (G.2)Proceeding as in Appendix C, we find (G.3)Introducing a system of spherical coordinates with the zaxis in the direction of w, and using Eq. (C.7), we obtain after some calculations (G.4)According to Eq. (C.4), we have (G.5)Therefore, Eqs. (G.4) and (G.5) lead to Eq. (29) with lnΛ replaced by (G.6)where t_{min} and t_{max} are appropriate cutoffs. The upper cutoff should be identified with the dynamical time t_{D}. On the other hand, the divergence at short times is due to the inadequacy of our assumption of straightline trajectories to describe very close encounters. If we take t_{min} = λ_{L}/v_{m} and t_{max} = λ_{J}/v_{m} ~ t_{D}, we find that lnΛ′ = lnΛ. In Appendix C, we have calculated K^{μν} by integrating first over time then over space. This yields the Coulombian logarithm (Eq. (30)). Here, we have integrated first over space then over time. This yields the Coulombian logarithm (Eq. (G.6)). As discussed by Lee (1968), these two approaches are essentially equivalent. We remark, however, that the calculations of Sect. C can be performed for arbitrary potentials, while the calculations of this section explicitly use the specific form of the gravitational potential.
According to Eqs. (G.1), (G.2), and (G.4), the force autocorrelation function can be written as (G.7)In particular (G.8)The τ^{1} decay of the autocorrelation function of the gravitational force was first derived by Chandrasekhar (1944) with a different method. This result has also been obtained, and discussed, by Cohen et al. (1950) and Lee (1968). According to Eqs. (43) and (G.7), the diffusion tensor is given by^{35}(G.9)This returns Eq. (79) with lnΛ′ instead of lnΛ.
When f(v_{1}) is the Maxwell distribution (Eq. (51)), we can compute the force autocorrelation tensor from Eq. (G.7) by using the Rosenbluth potentials as in Sect. 3.5. Alternatively, combining Eqs. (G.1) and (G.3), we have (G.10)where is the threedimensional Fourier transform of the distribution function. For the Maxwell distribution (Eq. (51)), we obtain (G.11)Introducing a system of spherical coordinates with the zaxis in the direction of v, and using Eq. (C.7), we obtain after some calculations (G.12)where G^{μν}(x) is defined in Sect. 3.3. In particular, (G.13)Integrating Eq. (G.12) over time, we recover the expression (Eq. (67)) of the diffusion tensor for a Maxwellian distribution with lnΛ′ instead of lnΛ.
Finally, the autocorrelation tensor of the gravitational field at two different points (at the same time) is (G.14)
Appendix H: The different kinetic equations
The standard kinetic equations (Boltzmann, FokkerPlanck, Vlasov, Landau, and LenardBalescu) and their generalizations can be derived from the BBGKY hierarchy. In this Appendix, we show the connection between these different equations, and discuss their domains of validity, without entering into technical details.
The first two equations of the BBGKY hierarchy may be written symbolically as where f is the onebody distribution function, g the twobody correlation function, and h the threebody correlation function. In the first equation, is the Vlasov operator taking into account the free motion of the particles and the advection by the mean field . On the other hand, the collision term C [g] describes the effect of twobody correlations on the evolution of the distribution function. In the second equation, ℒ = ℒ_{0} + ℒ′ + ℒ_{m.f.} is a twobody Liouville operator where ℒ_{0} describes the free motion of the particles, ℒ′ describes the exact twobody interaction, and ℒ_{m.f.} takes into account the effect of the mean field in the twobody problem. The term describes collective effects and the term describes threebody correlations. Finally, is a source term depending on the onebody distribution function.
As explained in Sect. 2.2, for N ≫ 1 we can expand the equations of the BBGKY hierarchy in terms of the small parameter 1/N. For N → +∞, the encounters are negligible and we get the Vlasov equation (Jeans 1915; Vlasov 1938). This corresponds to the mean field approximation. At the order 1/N, we can neglect threebody correlations () and strong collisions (ℒ′ = 0) that are of order 1/N^{2}. This corresponds to the weak coupling approximation. In fact, since the gravitational potential is singular at r = 0, the twobody correlation function g(r_{1}, r_{2}) becomes large when r_{1} − r_{2} → 0 due to the effect of strong collisions. As a result, it is necessary to take the term ℒ′g into account at small scales. When threebody correlations are neglected (), Eqs. (H.1) and (H.2) are closed. However, it does not appear possible to solve these equations explicitly without further approximation. The usual strategy is to solve these equations for small, intermediate, and large impact parameters, and then connect these limits.
If we neglect strong collisions (ℒ′ = 0), collective effects (), and make a local approximation (ℒ_{m.f.} = 0), we get the Landau equation (Eqs. (28) and (29)) with (H.3)This factor presents a logarithmic divergence at small and large scales. The Landau equation may be derived in different manners that are actually equivalent to the above procedure. Landau (1936) obtained his equation by starting from the Boltzmann equation and using a weak deflection approximation Δv ≪ 1. In his calculations, he approximated the trajectories of the particles by straight lines even for collisions with small impact parameters. This is equivalent to starting from the FokkerPlanck equation (Eq. (44)) and calculating the first and second moments of the velocity increments (Eq. (45)) resulting from a succession of binary encounters by using a straight line approximation.
If we neglect collective effects (), make a local approximation (ℒ_{m.f.} = 0), but take strong collisions into account (ℒ′ ≠ 0), we get the Boltzmann equation (see Balescu 2000). This equation does not present any divergence at small scales. Since the system is dominated by weak encounters, we can expand the Boltzmann equation for weak deflections Δv ≪ 1 while taking into account the effect of strong collisions. This is equivalent to the treatment of Chandrasekhar (1942, 1943a,b) who started from the FokkerPlanck equation (Eq. (44)) and calculated the first and second moments of the velocity increments (Eq. (45)) resulting from a succession of binary collisions by taking strong collisions into account. In the calculation of ⟨ Δv^{μ} ⟩ and ⟨ Δv^{μ}Δv^{ν} ⟩ , he used the exact trajectory of the stars (i.e. he solved the twobody problem exactly) and took into account the strong deflections due to collisions with small impact parameters. In the dominant approximation lnN ≫ 1, his approach, completed by Rosenbluth et al. (1957), leads to the FokkerPlanck equation (Eq. (44)) with the expressions of the diffusion tensor and friction force (Eqs. (83) and (84)) expressed in terms of the Rosenbluth potentials (Eq. (85)) with (H.4)As shown in Sect. 3.5, the resulting FokkerPlanck equation can be transformed into the Landau equation (Eqs. (28) and (29)). As in the Landau approach, the factor calculated in Eq. (H.4) presents a logarithmic divergence at large scales. However, contrary to the Landau approach, it does not present a logarithmic divergence at small scales since the effect of strong collisions is taken into account explicitly in the Chandrasekhar approach^{36}. As a result, the gravitational Landau length appears naturally in the calculations of Chandrasekhar. This establishes that the relevant smallscale cutoff in Eq. (H.3) is the Landau length.
The Landau and the Chandrasekhar kinetic theories make the same assumption: binary encounters and an expansion in powers of the momentum transfer. They actually differ in the order in which these are introduced. Landau starts from the Boltzmann equation and considers a weak deflection approximation while Chandrasekhar directly starts from the FokkerPlanck equation, but calculates the coefficients of diffusion and friction with the binarycollision picture. At that level, their theories are equivalent^{37} since the FokkerPlanck equation can be precisely obtained from the Boltzmann equation in the limit of weak deflections. The crucial difference is that Landau makes a weak coupling assumption and ignores strong collisions (i.e. the bending of the trajectories) while Chandrasekhar takes them into account. On the other hand, both theories ignore spatial inhomogeneity and collective effects.
If we neglect strong collisions (ℒ′ = 0) and collective effects (), but take spatial inhomogeneity (ℒ_{m.f.} ≠ 0) into account we get the generalized Landau equation (Eqs. (24) and (114)) derived in this paper (see also Kandrup 1981; and Chavanis 2008a,b). This equation presents a logarithmic divergence at small scales since strong collisions are neglected, but not at large scales since the finite extent of the system is taken into account. This suggests that the relevant largescale cutoff in Eqs. (H.3) and (H.4) is the Jeans length.
If we neglect strong collisions (ℒ′ = 0) but take collective effects () and spatial inhomogeneity (ℒ_{m.f.} ≠ 0) into account, we get the generalized LenardBalescu equation derived by Heyvaerts (2010) and Chavanis (2012a).
The ordinary Landau equation () is intermediate between the Boltzmann equation and the generalized LenardBalescu (and generalized Landau) equation. It describes the effect of weak collisions but ignores strong collisions, spatial inhomogeneity, and collective effects. It can be obtained from the generalized LenardBalescu equation (ℒ′ = 0) by making a local approximation (ℒ_{m.f.} = 0) and neglecting collective effects (), or from the Boltzmann equation () by considering the limit of small deflections (ℒ′ = 0).
Actually, these kinetic equations describe the effect of collisions at different scales (the scale λ may be interpreted as the impact parameter). For λ ~ λ_{L} (small impact parameters), the collisions are strong and we must solve the twobody problem exactly. For λ ~ l (intermediate impact parameters), the collisions are weak and we can make a weak coupling approximation. For λ ~ λ_{J} (large impact parameters), we must take spatial inhomogeneity and collective effects into account. The Boltzmann equation is valid for λ ≪ λ_{J}. It describes strong
collisions (λ ~ λ_{L}) and weak collisions (λ ~ l). The generalized Landau and generalized LenardBalescu equations are valid for λ ≫ λ_{L}. They describe weak collisions (λ ~ l), spatial inhomogeneity and collective effects (λ ~ λ_{J}). The ordinary Landau equation is valid for λ_{L} ≪ λ ≪ λ_{J}. It describes weak collisions. When we go beyond the domains of validity of these equations, divergences occur and appropriate cutoffs must be introduced.
Connecting these different limits, we find that the best description of stellar systems is provided by the generalized LenardBalescu equation with a smallscale cutoff at the Landau length. If we neglect collective effects, we get the generalized Landau equation (Eqs. (24) and (114)) with a smallscale cutoff at the Landau length. Finally, if we make a local approximation and neglect collective effects, the VlasovLandau equation (Eqs. (28) and (29)) with a smallscale cutoff at the Landau length and a largescale cutoff at the Jeans length provides a relevant description of stellar systems and has the advantage of the simplicity.
All Figures
Fig. 1 Normalized diffusion coefficients G_{ ∥ }(x), G_{ ⊥ }(x) and friction force xG_{ ∥ }(x) for a thermal bath. The friction is maximum for x ≃ 0.97, i.e. when the velocity of the test star is approximately equal to 0.792 times the rms velocity of the field stars (its value is (xG_{ ∥ })_{max} = 2.383). 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.