Vlasov–Poisson as an initial-value problem

Use a one-dimensional spatial model with one retained velocity coordinate. Let F describe number density per unit velocity, so its velocity integral is the physical three-dimensional number density. This reduced distribution can be obtained by integrating a full distribution over the two perpendicular velocities in an unmagnetized longitudinal problem. It is not normalized to one unless divided by local density. Its units differ from those of the original six-dimensional phase-space density, even though both are often called f.

The first velocity moment gives mean flow, while the central second moment gives parallel pressure after multiplication by mass. Subtracting the flow matters: a narrow beam moving rapidly can have large bulk kinetic energy but small pressure. Higher moments measure asymmetry and heat flow that are absent from density and pressure alone. A fluid model retaining only low moments therefore needs a closure; merely taking moments of Vlasov does not determine the unretained moments.

For electrons with charge minus e and an immobile uniform ion background, the electrostatic acceleration is minus eE divided by electron mass. The density perturbation sets the spatial derivative of E through Poisson's equation. On a periodic domain, integrating that equation requires global neutrality. An arbitrary nonneutral initial distribution is incompatible with a periodic electric field unless the model explicitly changes its background or includes a different boundary condition.

Poisson's equation determines the field only up to a spatially uniform component in this formulation. A common periodic electrostatic convention chooses zero mean E and fixes the additive constant of potential separately. The initial distribution must also remain nonnegative. A formal small perturbation that becomes more negative than the equilibrium tail somewhere is not an admissible particle density. For open boundaries, specify incoming distribution data consistently with the selected electrostatic boundary conditions rather than importing periodic assumptions.

Following a characteristic in position and velocity makes the distribution constant. The characteristic itself depends on the field generated by the whole distribution, so particles do not simply free stream unless that field is absent. The phase-space flow is incompressible: deformation can stretch and fold a distribution without erasing its fine structure. A numerical method that smooths those folds may change entropy even while its density plot looks accurate.

With suitable periodic or decaying boundaries, total particle number and the sum of electron kinetic and electrostatic field energy are conserved. Field energy can decrease while particle kinetic energy increases, so decay of an electric oscillation does not establish collisional dissipation. Retaining the evolving velocity structure is central to interpreting collisionless damping and later echoes. This conservation statement belongs to the ideal continuum model; finite resolution, boundary losses, or added collisions require their own accounting.