Dynamics of gravitational collapse in the axisymmetric Einstein-Vlasov system

We numerically investigate the dynamics 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 non-stationary 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.