  Minkowski space is shown to be globally stable as a solution to the massive
Einstein--Vlasov system. The proof is based on a harmonic gauge in which the
equations reduce to a system of quasilinear wave equations for the metric,
satisfying the weak null condition, coupled to a transport equation for the
Vlasov particle distribution function. Central to the proof is a collection of
vector fields used to control the particle distribution function, a function of
both spacetime and momentum variables. The vector fields are derived using a
general procedure, are adapted to the geometry of the solution and reduce to
the generators of the symmetries of Minkowski space when restricted to acting
on spacetime functions. Moreover, when specialising to the case of vacuum, the
proof provides a simplification of previous stability works.
