The charge breeder process make a natural representation of the test particle approach, which is an importanttopic of plasma physics: in detail, the exotic particles injected are the 'test particles', which propagates insidethe ECRIS plasma, which act as background: our particle motion description (theoretical and with numerical programs) includes full 3D details of collisions (affecting beam stopping and diffusion) and ECRIS field (solenoids, hexapoles, and ambipolar field models).At the same time, the progressive ionisation of particles (or recombination) is considered. The statistical properties of collisions are represented with an ensemble of particles. The numerical simulation confirm a three phase model of charge breeding: a) injection and slowing down in a typical length (also analytically estimated);b) diffusion mainly towards a central plasma core, due to ambipolar potential (this phase is longer, so mostionization happens here); c) random collision causes extraction and beam formation. Application includes: 1) optimisation of the charge breeder conditions; 2) validation of estimate of the background condition; 3) extension and comparison to other physics regimes. On the other side, the ionisation dynamic equilibrium that attains in phase b) can be studied with analytical model giving precise bounds on charge breeder efficiency, comparable with EBIS case. Therefore charge breeder studies leads to a better understanding of ECRIS physics. LNL-INFN is also involved in many research and projects on acceleration of exotic nuclear beams.