We numerically investigate the dynamcis near black hole formation of solutions to the Einstein-Vlasov system in axisymmetry. Our results are obtained using a particle-in-cell and finite difference code based on the (2+1)+1 formulation of the Einstein field equations in axisymmetry. Solutions are launched from generic type initial data and exhibit type I critical behaviour. In particular we find lifetime scaling in solutions containing black holes, and support that the critical solutions are stationary. Our results contain examples of solutions that form black holes, perform damped oscillations, and appear to disperse. We prove that complete dispersal of the solution implies that it has nonpositive binding energy.
