This work is concerned with the study of singular limits for the Vlasov-Poisson system in the case of massless electrons (VPME), which is a kinetic system modelling the ions in a plasma. Our objective is threefold: first, we provide a mean field derivation of the VPME system in dimensions d = 2, 3 from a system of N extended charges. Secondly, we prove a rigorous quasineutral limit for initial data that are perturbations of analytic data, deriving the Kinetic Isothermal Euler (KIE) system from the VPME system in dimensions d = 2, 3. Lastly, we combine these two singular limits in order to show how to obtain the KIE system from an underlying particle system. scaling, the VPME system becomesIn real plasmas, the Debye length is typically very small. In this case, the plasma is called quasineutral, in reference to the fact that the plasma appears to be neutral overall at the observation scale. In the physics literature, quasineutrality is often included in the very definition of plasma. It is therefore interesting to consider the limit in which ε tends to zero. The formal limit of (1.3) is the kinetic isothermal Euler system:(1.4) A related equation, derived in the quasineutral limit from a version of (1.3) in which the electron density e U is approximated by the linearisation 1 + U , was named Vlasov-Dirac-Benney by Bardos and studied in [2] and [4]. In the classical case, where E is replaced by a pressure term which is a Lagrange multiplier corresponding to the constraint ρ = 1, Bossy-Fontbona-Jabin-Jabir [6] showed local-in-time existence of analytic solutions in the one-dimensional case. Global existence of weak solutions is not known for (1.4).The rigorous justification of the quasineutral limit is a non-trivial and subtle problem. The limit has a direct correspondence to a long-time limit for the Vlasov-Poisson system, and is therefore vulnerable to known instability mechanisms inherent to the physical system under consideration. In fact, for the classical system (1.1), it was shown by Hauray and Han-Kwan in [19] that the quasineutral limit is false in general if the initial data are assumed to have only Sobolev regularity.Rigorous results on the quasineutral limit go back to the works of Brenier-Grenier [9] and Grenier [15] for the classical system (1.1). A result of particular relevance for our purposes is the work of Grenier [16], proving the limit for the classical system assuming uniformly analytic data. The works of Han-Kwan-Iacobelli [20,21] extended this result to data that are very small, but possibly rough, perturbations of the uniformly analytic case, in dimension 1, 2 and 3.In the massless electrons case, Han-Kwan-Iacobelli [21] showed a rigorous limit in dimension one, again for rough perturbations of analytic data, while Han-Kwan-Rousset [22] consider Penrose-stable data with sufficiently high Sobolev regularity. In this work, we extend the results of [21] to higher dimensions by showing a rigorous quasineutral limit for the VPME system (1.3) in dimension 2 and 3, for data that are very small, but...