We introduce a fractional Klein-Kramers equation which describes sub-ballistic superdiffusion in phase space in the presence of a space-dependent external force field. This equation defines the differential Lévy walk model whose solution is shown to be non-negative. In the velocity coordinate, the probability density relaxes in Mittag-Leffler fashion towards the Maxwell distribution whereas in the space coordinate, no stationary solution exists and the temporal evolution of moments exhibits a competition between Brownian and anomalous contributions. Classically, Brownian stochastic transport processes in the phase space spanned by velocity v and coordinate x are described by the deterministic Klein-Kramers equation (KKE). In the low and high friction limits, the KKE reduces to the Rayleigh equation which describes the relaxation of the velocity probability density function (pdf) towards the Maxwell distribution, and the FokkerPlanck-Smoluchowski equation controlling the temporal approach of the Gibbs-Boltzmann equilibrium, respectively [1, 2]. The KKE is therefore a fundamental equation in non-equilibrium systems dynamics.Brownian transport is characterised through the Gaussian pdf and the linear time dependence x 2 (t) = 2Kt of the mean squared displacement in the force-free diffusion limit, its universality being guaranteed by the central limit theorem [2]. In a broad variety of systems, however, it has been found that correlations in space or time give birth to anomalous transport whose pdf is non-Gaussian and/or whose mean squared displacement is non-linear in time [3]. These systems include charge carrier transport in amorphous semiconductors [4], tracer dispersion in convection rolls and rotating flows [5,6] The systems we are interested in fall into the broad class whose force-free diffusion behaviour is characterised through the power-law form x 2 (t) ∝ t κ , which separates into subdiffusion (0 < κ < 1) and superdiffusion (κ > 1). The continuous time random walk (CTRW) model has proved to be a well-suited framework which accounts for such anomalous diffusion for the entire spectrum of κ [14]. Especially in the sub-ballistic superdiffusive domain 1 < κ < 2, Lévy walks which couple long flight times with a time cost have been a successful tool, e.g., in fluid dynamics [6,15]. The space-time coupling of Lévy walks leads to finite moments and they are therefore fundamentally different from Lévy flights which exhibit a diverging variance [16].In the presence of external force fields, the CTRW approach is less flexible. It has been realised that fractional equations constitute a tailor-made framework to formulate the underlying dynamics equations in coordinate and phase space; see, for instance, [16,17,18,19,20] and references therein. In the subdiffusive domain 0 < κ < 1, a fractional KKE was derived from CTRW models [21] and from the Chapman-Kolmogorov equation [22]. A similar consistent generalisation to systems in the regime 1 < κ < 2 is still outstanding. In this note, we propose the fractional KKE