A paraxial model of an unmagnetized, collisionless plasma plume expanding into vacuum is presented. Electrons are treated kinetically, relying on the adiabatic invariance of their radial action integral for the integration of Vlasov's equation, whereas ions are treated as a cold species. The quasi-2D plasma density; self-consistent electric potential; electron pressure, temperature, and heat fluxes are analyzed. In particular, the model yields the collisionless cooling of electrons, which differs from the Boltzmann relation and the simple polytropic laws usually employed in fluid and hybrid PIC/fluid plume codes.This article presents a quasi-2D kinetic model of a plasma plume and uses it to characterize several features of the electron expansion. The model assumes a collisionless, unmagnetized, quasineutral, steady-state plasma, providing a good description of the plume far-region. The integration of the electron model relies on the first-order