An axisymmetric collapse of non-rotating gravitational waves is numerically investigated in the subcritical regime where no black holes form but where curvature attains a maximum and decreases, following the dispersion of the initial wave packet. We focus on a curvature invariant with dimensions of length, and find that near the threshold for black hole formation it reaches a maximum along concentric rings of finite radius around the axis. In this regime the maximal value of the invariant exhibits a power-law scaling with the approximate exponent 0.38, as a function of a parametric distance from the threshold. In addition, the variation of the curvature in the critical limit is accompanied by increasing amount of echos, with nearly equal temporal and spatial periods. The scaling and the echoing patterns, and the corresponding constants, are independent of the initial data and coordinate choices.

We describe the first axisymmetric numerical code based on the generalized harmonic formulation of the Einstein equations which is regular at the axis. We test the code by investigating gravitational collapse of distributions of complex scalar field in a Kaluza-Klein spacetime. One of the key issues of the harmonic formulation is the choice of the gauge source functions, and we conclude that a damped wave gauge is remarkably robust in this case. Our preliminary study indicates that evolution of regular initial data leads to formation both of black holes with spherical and cylindrical horizon topologies. Intriguingly, we find evidence that near threshold for black hole formation the number of outcomes proliferates. Specifically, the collapsing matter splits into individual pulses, two of which travel in the opposite directions along the compact dimension and one which is ejected radially from the axis. Depending on the initial conditions, a curvature singularity develops inside the pulses.

We investigate the effects of pair creation on the internal geometry of a black hole, which forms during the gravitational collapse of a charged massless scalar field. Classically, strong central Schwarzschild-like singularity forms, and a null, weak, mass-inflation singularity arises along the Cauchy horizon, in such a collapse. We consider here the discharge, due to pair creation, below the event horizon and its influence on the {\it dynamical formation} of the Cauchy horizon. Within the framework of a simple model we are able to trace numerically the collapse. We find that a part of the Cauchy horizon is replaced by the strong space-like central singularity. This fraction depends on the value of the critical electric field, $E_{\rm cr}$, for the pair creation.