We prove the stability and convergence of the high order discontinuous Galerkin scheme to spherically symmetric Einstein-scalar equations for a class of large initial data that ensures the formation of a black hole. Having chosen the Bondi coordinate system, we achieve L2 stability and obtain the optimal error estimates.

We construct a non-polynomial local discontinuous Galerkin (LDG) scheme for the prescribed mean curvature equation to approximate boundary gradient blow-up solutions and obtain error estimates.

In this paper, we show the numerical solution for spherically symmetric SU(2) EinsteinYang-Mills (EYM) equations. We show the existence of entropy weak solution for EYM.

In this letter, we find the first dynamically stable non-singular solution spherically symmetric SU(2) Einstein-Yang-Mills equation. This solutions is regular at r=0 and asymptotically flat. Since the Yang-Mills field strength decay exponentially, the Einstein-Yang-Mills particle perhaps can be used to be a candidate for dark matter.

We present stable solution of static spherically symmetric Einstein-Yang-Mills equations with the SU(2) gauge group. This solution is asymptotically flat and regular at r = 0 and with nontrivial Yang-Mills(YM) connection. With quantized values of the Arnowitt-Deser-Misner (ADM) mass, the solutions asymptotically approach the Schwarzschild solution and have zero global YM charges. Numerical evidences suggest that this solution is both linearly and nonlinearly stable and has a ring of generic curvature singularities along the horizon. An effective counterexample to the no-hair conjecture is provided by this stable solution. Moreover, the stable black hole solution suggests that the coupling of gauge field to gravity in early Universe will generate a new type of black holes. Their stability means that these might be a possible new source of primordial black holes left over from the early Universe and serves as a possible new candidate for dark matter.

The critical collapse of a scalar field is a threshold solution of black hole formation, in which a naked singularity arises. We study here the curvature strength of this singularity using a numerical ansatz. The behavior of the Jacobi volume forms is examined along a non-spacelike geodesic in the limit of approach to the singularity. These are seen to be vanishing, thus showing that all physical objects will be crushed to zero size near the singularity. Consequently, although the critical collapse is considered to be a fine-tuning problem, the naked singularity forming is gravitationally strong. This implies that the spacetime cannot be extended beyond the singularity, thus making the singularity genuine and physically interesting. These results imply that the nature of critical collapse may need to be examined and explored further.