# Boltzmann equation

(Redirected from Boltzmann transport equation)
"BTE" redirects here. For other uses, see BTE (disambiguation).
The place of the Boltzmann kinetic equation on the stairs of model reduction from microscopic dynamics to macroscopic continuum dynamics (Illustration to the content of the book[1]).

The Boltzmann equation or Boltzmann transport equation (BTE) describes the statistical behaviour of a thermodynamic system not in a state of equilibrium, devised by Ludwig Boltzmann in 1872.[2] The classic example of such a system is a fluid with temperature gradients in space causing heat to flow from hotter regions to colder ones, by the random but biased transport of the particles making up that fluid. In the modern literature the term Boltzmann equation is often used in a more general sense, referring to any kinetic equation that describes the change of a macroscopic quantity in a thermodynamic system, such as energy, charge or particle number.

The equation arises not by analyzing the individual positions and momenta of each particle in the fluid but rather by considering a probability distribution for the position and momentum of a typical particle—that is, the probability that the particle occupies a given very small region of space (mathematically written ${\displaystyle d^{3}{\bf {r}}}$, where d means "differential") centered at the position ${\displaystyle {\bf {r}}}$, and has momentum nearly equal to a given momentum vector ${\displaystyle {\bf {p}}}$ (thus occupying a very small region of momentum space ${\displaystyle d^{3}{\bf {p}}}$), at an instant of time.

The Boltzmann equation can be used to determine how physical quantities change, such as heat energy and momentum, when a fluid is in transport. One may also derive other properties characteristic to fluids such as viscosity, thermal conductivity, and electrical conductivity (by treating the charge carriers in a material as a gas).[2] See also convection-diffusion equation.

The equation is a nonlinear integro-differential equation, and the unknown function in the equation is a probability density function in six-dimensional space of a particle position and momentum. The problem of existence and uniqueness of solutions is still not fully resolved, but some recent results are quite promising.[3][4]

## Overview

### The phase space and density function

The set of all possible positions r and momenta p is called the phase space of the system; in other words a set of three coordinates for each position coordinate x, y, z, and three more for each momentum component px, py, pz. The entire space is 6-dimensional: a point in this space is (r, p) = (x, y, z, px, py, pz), and each coordinate is parameterized by time t. The small volume ("differential volume element") is written

${\displaystyle {\text{d}}^{3}\mathbf {r} \,{\text{d}}^{3}\mathbf {p} ={\text{d}}x\,{\text{d}}y\,{\text{d}}z\,{\text{d}}p_{x}\,{\text{d}}p_{y}\,{\text{d}}p_{z}.}$

Since the probability of N molecules which all have r and p within ${\displaystyle d^{3}{\bf {r}}}$ ${\displaystyle d^{3}{\bf {p}}}$ is in question, at the heart of the equation is a quantity f which gives this probability per unit phase-space volume, or probability per unit length cubed per unit momentum cubed, at an instant of time t. This is a probability density function: f(r, p, t), defined so that,

${\displaystyle {\text{d}}N=f(\mathbf {r} ,\mathbf {p} ,t)\,{\text{d}}^{3}\mathbf {r} \,{\text{d}}^{3}\mathbf {p} }$

is the number of molecules which all have positions lying within a volume element ${\displaystyle d^{3}{\bf {r}}}$ about r and momenta lying within a momentum space element ${\displaystyle d^{3}{\bf {p}}}$ about p, at time t.[5] Integrating over a region of position space and momentum space gives the total number of particles which have positions and momenta in that region:

{\displaystyle {\begin{aligned}N&=\int \limits _{\mathrm {positions} }{\text{d}}^{3}\mathbf {r} \int \limits _{\mathrm {momenta} }{\text{d}}^{3}\mathbf {p} \,f(\mathbf {r} ,\mathbf {p} ,t)\\[5pt]&=\iiint \limits _{\mathrm {positions} }\quad \iiint \limits _{\mathrm {momenta} }f(x,y,z,p_{x},p_{y},p_{z},t)\,{\text{d}}x\,{\text{d}}y\,{\text{d}}z\,{\text{d}}p_{x}\,{\text{d}}p_{y}\,{\text{d}}p_{z}\end{aligned}}}

which is a 6-fold integral. While f is associated with a number of particles, the phase space is for one-particle (not all of them, which is usually the case with deterministic many-body systems), since only one r and p is in question. It is not part of the analysis to use r1, p1 for particle 1, r2, p2 for particle 2, etc. up to rN, pN for particle N.

It is assumed the particles in the system are identical (so each has an identical mass m). For a mixture of more than one chemical species, one distribution is needed for each, see below.

### Principal statement

The general equation can then be written:[6]

${\displaystyle {\frac {\partial f}{\partial t}}=\left({\frac {\partial f}{\partial t}}\right)_{\mathrm {force} }+\left({\frac {\partial f}{\partial t}}\right)_{\mathrm {diff} }+\left({\frac {\partial f}{\partial t}}\right)_{\mathrm {coll} }}$

where the "force" term corresponds to the forces exerted on the particles by an external influence (not by the particles themselves), the "diff" term represents the diffusion of particles, and "coll" is the collision term – accounting for the forces acting between particles in collisions. Expressions for each term on the right side are provided below.[6]

Note that some authors use the particle velocity v instead of momentum p; they are related in the definition of momentum by p = mv.

## The force and diffusion terms

Consider particles described by f, each experiencing an external force F not due to other particles (see the collision term for the latter treatment).

Suppose at time t some number of particles all have position r within element ${\displaystyle d^{3}{\bf {r}}}$ and momentum p within ${\displaystyle d^{3}{\bf {p}}}$. If a force F instantly acts on each particle, then at time t + Δt their position will be r + Δr = r + pΔt/m and momentum p + Δp = p + FΔt. Then, in the absence of collisions, f must satisfy

${\displaystyle f\left(\mathbf {r} +{\frac {\mathbf {p} }{m}}\,\Delta t,\mathbf {p} +\mathbf {F} \,\Delta t,t+\Delta t\right)\,d^{3}\mathbf {r} \,d^{3}\mathbf {p} =f(\mathbf {r} ,\mathbf {p} ,t)\,d^{3}\mathbf {r} \,d^{3}\mathbf {p} }$

Note that we have used the fact that the phase space volume element ${\displaystyle d^{3}{\bf {r}}}$ ${\displaystyle d^{3}{\bf {p}}}$ is constant, which can be shown using Hamilton's equations (see the discussion under Liouville's theorem). However, since collisions do occur, the particle density in the phase-space volume ${\displaystyle d^{3}{\bf {r}}}$ '${\displaystyle d^{3}{\bf {p}}}$ changes, so

{\displaystyle {\begin{aligned}dN_{\mathrm {coll} }&=\left({\frac {\partial f}{\partial t}}\right)_{\mathrm {coll} }\Delta td^{3}\mathbf {r} d^{3}\mathbf {p} \\[5pt]&=f\left(\mathbf {r} +{\frac {\mathbf {p} }{m}}\Delta t,\mathbf {p} +\mathbf {F} \Delta t,t+\Delta t\right)d^{3}\mathbf {r} d^{3}\mathbf {p} -f(\mathbf {r} ,\mathbf {p} ,t)\,d^{3}\mathbf {r} \,d^{3}\mathbf {p} \\[5pt]&=\Delta f\,d^{3}\mathbf {r} \,d^{3}\mathbf {p} \end{aligned}}}

(1)

where Δf is the total change in f. Dividing (1) by ${\displaystyle d^{3}{\bf {r}}}$ ${\displaystyle d^{3}{\bf {p}}}$ Δt and taking the limits Δt → 0 and Δf → 0, we have

${\displaystyle {\frac {df}{dt}}=\left({\frac {\partial f}{\partial t}}\right)_{\mathrm {coll} }}$

(2)

The total differential of f is:

{\displaystyle {\begin{aligned}df&={\frac {\partial f}{\partial t}}\,dt+\left({\frac {\partial f}{\partial x}}\,dx+{\frac {\partial f}{\partial y}}\,dy+{\frac {\partial f}{\partial z}}\,dz\right)+\left({\frac {\partial f}{\partial p_{x}}}\,dp_{x}+{\frac {\partial f}{\partial p_{y}}}\,dp_{y}+{\frac {\partial f}{\partial p_{z}}}\,dp_{z}\right)\\[5pt]&={\frac {\partial f}{\partial t}}dt+\nabla f\cdot d\mathbf {r} +{\frac {\partial f}{\partial \mathbf {p} }}\cdot d\mathbf {p} \\[5pt]&={\frac {\partial f}{\partial t}}dt+\nabla f\cdot {\frac {\mathbf {p} dt}{m}}+{\frac {\partial f}{\partial \mathbf {p} }}\cdot \mathbf {F} \,dt\end{aligned}}}

(3)

where ∇ is the gradient operator, · is the dot product,

${\displaystyle {\frac {\partial f}{\partial \mathbf {p} }}=\mathbf {\hat {e}} _{x}{\frac {\partial f}{\partial p_{x}}}+\mathbf {\hat {e}} _{y}{\frac {\partial f}{\partial p_{y}}}+\mathbf {\hat {e}} _{z}{\frac {\partial f}{\partial p_{z}}}=\nabla _{\mathbf {p} }f}$

is a shorthand for the momentum analogue of ∇, and êx, êy, êz are Cartesian unit vectors.

### Final statement

Dividing (3) by dt and substituting into (2) gives:

${\displaystyle {\frac {\partial f}{\partial t}}+{\frac {\mathbf {p} }{m}}\cdot \nabla f+\mathbf {F} \cdot {\frac {\partial f}{\partial \mathbf {p} }}=\left({\frac {\partial f}{\partial t}}\right)_{\mathrm {coll} }}$

In this context, F(r, t) is the force field acting on the particles in the fluid, and m is the mass of the particles. The term on the right hand side is added to describe the effect of collisions between particles; if it is zero then the particles do not collide. The collisionless Boltzmann equation is often called the Vlasov equation.

This equation is more useful than the principal one above, yet still incomplete, since f cannot be solved unless the collision term in f is known. This term cannot be found as easily or generally as the others – it is a statistical term representing the particle collisions, and requires knowledge of the statistics the particles obey, like the Maxwell–Boltzmann, Fermi–Dirac or Bose–Einstein distributions.

## The collision term (Stosszahlansatz) and molecular chaos

A key insight applied by Boltzmann was to determine the collision term resulting solely from two-body collisions between particles that are assumed to be uncorrelated prior to the collision. This assumption was referred to by Boltzmann as the "Stosszahlansatz", and is also known as the "molecular chaos assumption". Under this assumption the collision term can be written as a momentum-space integral over the product of one-particle distribution functions:[2]

${\displaystyle \left({\frac {\partial f}{\partial t}}\right)_{\mathrm {coll} }=\iint gI(g,\Omega )[f(\mathbf {p'} _{A},t)f(\mathbf {p'} _{B},t)-f(\mathbf {p} _{A},t)f(\mathbf {p} _{B},t)]\,d\Omega \,d^{3}\mathbf {p} _{A}\,d^{3}\mathbf {p} _{B}.}$

where pA and pB are the momenta of any two particles (labeled as A and B for convenience) before a collision, p′A and p′B are the momenta after the collision,

${\displaystyle g=|\mathbf {p} _{B}-\mathbf {p} _{A}|=|\mathbf {p'} _{B}-\mathbf {p'} _{A}|}$

is the magnitude of the relative momenta (see relative velocity for more on this concept), and I(g, Ω) is the differential cross section of the collision, in which the relative momenta of the colliding particles turns through an angle θ into the element of the solid angle dΩ, due to the collision.

## Simplifications to the Boltzmann equation collision term

Since much of the challenge in solving the Boltzmann equation originates with the complex collision term, attempts have been made to 'model' and simplify the collision term. The best known model equation is due to Bhatnagar, Gross and Krook.[7] The assumption in the BGK approximation is that the effect of molecular collisions is to force a non-equilibrium distribution function at a point in physical space back to a Maxwellian equilibrium distribution function and that the rate at which this occurs is proportional to the molecular collision frequency. The Boltzmann equation is therefore modified to the BGK form,

${\displaystyle {\frac {\partial f}{\partial t}}+{\frac {\mathbf {p} }{m}}\cdot \nabla f+\mathbf {F} \cdot {\frac {\partial f}{\partial \mathbf {p} }}=\nu (f_{0}-f)}$

where ${\displaystyle \nu }$ is the molecular collision frequency and ${\displaystyle f_{0}}$ is the local Maxwellian distribution function given the gas temperature at this point in space.

## General equation (for a mixture)

For a mixture of chemical species labelled by indices i = 1,2,3...,n the equation for species i is:[2]

${\displaystyle {\frac {\partial f_{i}}{\partial t}}+{\frac {\mathbf {p} _{i}}{m_{i}}}\cdot \nabla f_{i}+\mathbf {F} \cdot {\frac {\partial f_{i}}{\partial \mathbf {p} _{i}}}=\left({\frac {\partial f_{i}}{\partial t}}\right)_{\mathrm {coll} }}$

where fi = fi(r, pi, t), and the collision term is

${\displaystyle \left({\frac {\partial f_{i}}{\partial t}}\right)_{\mathrm {coll} }=\sum _{j=1}^{n}\iint g_{ij}I_{ij}(g_{ij},\Omega )[f'_{i}f'_{j}-f_{i}f_{j}]\,d\Omega \,d^{3}\mathbf {p'} .}$

where f′ = f′(p′i, t), the magnitude of the relative momenta is

${\displaystyle g_{ij}=|\mathbf {p} _{i}-\mathbf {p} _{j}|=|\mathbf {p'} _{i}-\mathbf {p'} _{j}|}$

and Iij is the differential cross-section as before, between particles i and j. The integration is over the momentum components in the integrand (which are labelled i and j). The sum of integrals describes the entry and exit of particles of species i in or out of the phase space element.

## Applications and extensions

### Conservation equations

The Boltzmann equation can be used to derive the fluid dynamic conservation laws for mass, charge, momentum and energy.[8]:p 163 For a fluid consisting of only one kind of particle, the number density n is given by:

${\displaystyle n=\int f\,d^{3}p}$

The average value of any function A is:

${\displaystyle \langle A\rangle ={\frac {1}{n}}\int Af\,d^{3}p}$

Since the conservation equations involve tensors, the Einstein summation convention will be used where repeated indices in a product indicate summation over those indices. Thus ${\displaystyle \mathbf {x} \rightarrow x_{i}}$ and ${\displaystyle \mathbf {p} \rightarrow p_{i}=mw_{i}}$ where ${\displaystyle w_{i}}$ is the particle velocity vector. Define ${\displaystyle g(p_{i})}$ as some function of momentum ${\displaystyle p_{i}}$ only, which is conserved in a collision. Assume also that the force ${\displaystyle F_{i}}$ is a function of position only, and that f is zero for ${\displaystyle p_{i}\rightarrow \pm \infty }$. Multiplying the Boltzmann equation by g and integrating over momentum yields four terms which, using integration by parts, can be expressed as:

${\displaystyle \int g{\frac {\partial f}{\partial t}}\,d^{3}p={\frac {\partial }{\partial t}}(n\langle g\rangle )}$
${\displaystyle \int {\frac {p_{j}g}{m}}{\frac {\partial f}{\partial x_{j}}}\,d^{3}p={\frac {1}{m}}{\frac {\partial }{\partial x_{j}}}(n\langle gp_{j}\rangle )}$
${\displaystyle \int gF_{j}{\frac {\partial f}{\partial p_{j}}}\,d^{3}p=-nF_{j}\left\langle {\frac {\partial g}{\partial p_{j}}}\right\rangle }$
${\displaystyle \int g\left({\frac {\partial f}{\partial t}}\right)_{\mathrm {coll} }\,d^{3}p=0}$

where the last term is zero since g is conserved in a collision. Letting ${\displaystyle g=m}$, the mass of the particle, the integrated Boltzmann equation becomes the conservation of mass equation:[8]:pp 12,168

${\displaystyle {\frac {\partial }{\partial t}}\rho +{\frac {\partial }{\partial x_{j}}}(\rho V_{j})=0}$

where ${\displaystyle \rho =mn}$ is the mass density and ${\displaystyle V_{i}=\langle w_{i}\rangle }$ is the average fluid velocity.

Letting ${\displaystyle g=mw_{i}}$, the momentum of the particle, the integrated Boltzmann equation becomes the conservation of momentum equation:[8]:pp 15,169

${\displaystyle {\frac {\partial }{\partial t}}(\rho V_{i})+{\frac {\partial }{\partial x_{j}}}(\rho V_{i}V_{j}+P_{ij})-nF_{i}=0}$

where ${\displaystyle P_{ij}=\rho \langle (w_{i}-V_{i})(w_{j}-V_{j})\rangle }$ is the pressure tensor. (The viscous stress tensor plus the hydrostatic pressure.)

Letting ${\displaystyle g={\tfrac {1}{2}}mw_{i}w_{i}}$, the kinetic energy of the particle, the integrated Boltzmann equation becomes the conservation of energy equation:[8]:pp 19,169

${\displaystyle {\frac {\partial }{\partial t}}(u+{\tfrac {1}{2}}\rho V_{i}V_{i})+{\frac {\partial }{\partial x_{j}}}(uV_{j}+{\tfrac {1}{2}}\rho V_{i}V_{i}V_{j}+J_{qj}+P_{ij}V_{i})-nF_{i}V_{i}=0}$

where ${\displaystyle u={\tfrac {1}{2}}\rho \langle (w_{i}-V_{i})(w_{i}-V_{i})\rangle }$ is the kinetic thermal energy density and ${\displaystyle J_{qi}={\tfrac {1}{2}}\rho \langle (w_{i}-V_{i})(w_{k}-V_{k})(w_{k}-V_{k})\rangle }$ is the heat flux vector.

### Hamiltonian mechanics

In Hamiltonian mechanics, the Boltzmann equation is often written more generally as

${\displaystyle {\hat {\mathbf {L} }}[f]=\mathbf {C} [f],\,}$

where L is the Liouville operator describing the evolution of a phase space volume and C is the collision operator. The non-relativistic form of L is

${\displaystyle {\hat {\mathbf {L} }}_{\mathrm {NR} }={\frac {\partial }{\partial t}}+{\frac {\mathbf {p} }{m}}\cdot \nabla +\mathbf {F} \cdot {\frac {\partial }{\partial \mathbf {p} }}\,.}$

### Quantum theory and violation of particle number

It is possible to write down relativistic Boltzmann equations for relativistic quantum systems in which the number of particles is not conserved in collisions. This has several applications in physical cosmology,[9] including the formation of the light elements in big bang nucleosynthesis, the production of dark matter and baryogenesis. It is not a priori clear that the state of a quantum system can be characterized by a classical phase space density f. However, for a wide class of applications a well-defined generalization of f exists which is the solution of an effective Boltzmann equation that can be derived from first principles of quantum field theory.[10]

### General relativity and astronomy

The Boltzmann equation is of use in galactic dynamics. A galaxy, under certain assumptions, may be approximated as a continuous fluid; its mass distribution is then represented by f; in galaxies, physical collisions between the stars are very rare, and the effect of gravitational collisions can be neglected for times far longer than the age of the universe.

Its generalization in general relativity is

${\displaystyle {\hat {\mathbf {L} }}_{\mathrm {GR} }=p^{\alpha }{\frac {\partial }{\partial x^{\alpha }}}-\Gamma ^{\alpha }{}_{\beta \gamma }p^{\beta }p^{\gamma }{\frac {\partial }{\partial p^{\alpha }}},}$

where Γαβγ is the Christoffel symbol of the second kind (this assumes there are no external forces, so that particles move along geodesics in the absence of collisions), with the important subtlety that the density is a function in mixed contravariant-covariant (xi, pi) phase space as opposed to fully contravariant (xi, pi) phase space.[11][12]

In physical cosmology, the study of processes in the early universe often requires to take into account the effects of quantum mechanics and general relativity.[9] In the very dense medium formed by the primordial plasma after the big bang, particles are continuously created and annihilated. In such an environment quantum coherence and the spatial extension of the wavefunction can affect the dynamics, making it questionable whether the classical phase space distribution f that appears in the Boltzmann equation is suitable to describe the system. In many cases it is, however, possible to derive an effective Boltzmann equation for a generalized distribution function from first principles of quantum field theory.[10] This includes the formation of the light elements in big bang nucleosynthesis, the production of dark matter and baryogenesis.

## Solutions to the Boltzmann equation

It was proven only in 2010 that exact solutions to the Boltzmann equation are always mathematically well-behaved. This means that if a system obeying the Boltzmann equation is perturbed, then it will return to equilibrium, rather than diverging to infinity or behaving otherwise.[13][14] However, this existence proof is not helpful for solving the equation in realistic scenarios. Indeed, such statements only tell us whether the solution subject to specified conditions exist, but not how to find them. In practice, numerical methods are used to find approximate solutions to the various forms of the Boltzmann equation with applications ranging from hypersonic aerodynamics in rarefied gas flows[15] to plasma flows.[16]

## Notes

1. ^ A.N. Gorban, I.V. Karlin,Invariant manifolds for physical and chemical kinetics, Lecture Notes in Physics, January 2005.
2. ^ a b c d Encyclopaedia of Physics (2nd Edition), R. G. Lerner, G. L. Trigg, VHC publishers, 1991, ISBN (Verlagsgesellschaft) 3-527-26954-1, ISBN (VHC Inc.) 0-89573-752-3
3. ^ DiPerna, R. J.; Lions, P.-L. (1989). "On the Cauchy problem for Boltzmann equations: global existence and weak stability". Ann. of Math. (2). 130 (2): 321–366. doi:10.2307/1971423.
4. ^ Philip T. Gressman & Robert M. Strain (2010). "Global classical solutions of the Boltzmann equation with long-range interactions". Proceedings of the National Academy of Sciences. 107 (13): 5744–5749. arXiv:. Bibcode:2010PNAS..107.5744G. doi:10.1073/pnas.1001185107. PMC . PMID 20231489.
5. ^ Huang, Kerson (1987). Statistical Mechanics (Second ed.). New York: Wiley. p. 53. ISBN 0-471-81518-7.
6. ^ a b McGraw Hill Encyclopaedia of Physics (2nd Edition), C.B. Parker, 1994, ISBN 0-07-051400-3
7. ^ Bhatnagar, P. L.; Gross, E. P.; Krook, M. (1954-05-01). "A Model for Collision Processes in Gases. I. Small Amplitude Processes in Charged and Neutral One-Component Systems". Physical Review. 94 (3): 511–525. Bibcode:1954PhRv...94..511B. doi:10.1103/PhysRev.94.511.
8. ^ a b c d de Groot, S.R.; Mazur, P. (1984). Non-Equilibrium Thermodynamics. New York: Dover Publications Inc. ISBN 0-486-64741-2. Retrieved 2013-01-31.
9. ^ a b Edward Kolb & Michael Turner (1990). The Early Universe. Westview Press. ISBN 9780201626742.
10. ^ a b M. Drewes; C. Weniger; S. Mendizabal (8 January 2013). "The Boltzmann equation from quantum field theory". Phys. Lett. B. 718 (3): 1119–1124. arXiv:. Bibcode:2013PhLB..718.1119D. doi:10.1016/j.physletb.2012.11.046.
11. ^ Debbasch, Fabrice; Willem van Leeuwen (2009). "General relativistic Boltzmann equation I: Covariant treatment". Physica A. 388 (7): 1079–1104. Bibcode:2009PhyA..388.1079D. doi:10.1016/j.physa.2008.12.023.
12. ^ Debbasch, Fabrice; Willem van Leeuwen (2009). "General relativistic Boltzmann equation II: Manifestly covariant treatment". Physica A. 388 (9): 1818–34. Bibcode:2009PhyA..388.1818D. doi:10.1016/j.physa.2009.01.009.
13. ^ Philip T. Gressman & Robert M. Strain (2010). "Global classical solutions of the Boltzmann equation with long-range interactions". PNAS. 107: 5744. Bibcode:2010PNAS..107.5744G. doi:10.1073/pnas.1001185107. PMC . PMID 20231489.
14. ^ "Mathematicians Solve 140-Year-Old Boltzmann Equation". news.upenn.edu. Retrieved 2016-01-20.
15. ^ Evans, Ben; Morgan, Ken; Hassan, Oubay (2011-03-01). "A discontinuous finite element solution of the Boltzmann kinetic equation in collisionless and BGK forms for macroscopic gas flows". Applied Mathematical Modelling. 35 (3): 996–1015. doi:10.1016/j.apm.2010.07.027.
16. ^ Pareschi, L.; Russo, G. (2000-01-01). "Numerical Solution of the Boltzmann Equation I: Spectrally Accurate Approximation of the Collision Operator". SIAM Journal on Numerical Analysis. 37 (4): 1217–1245. doi:10.1137/S0036142998343300. ISSN 0036-1429.

## References

• Harris, Stewart (1971). An introduction to the theory of the Boltzmann equation. Dover Books. p. 221. ISBN 978-0-486-43831-3.. Very inexpensive introduction to the modern framework (starting from a formal deduction from Liouville and BBGKY) in which the Boltzmann equation is placed. Most of statistical mechanics textbooks like Huang even treat the topic with the original arguments used by Boltzmann. For what concerns the deduction of the equation, these books rather use a heuristic explanation hiding the range of validity and the characteristic assumptions that distinguish Boltzmann's from other transport equations like Fokker–Planck or Landau equations.