Solving Einstein's Vacuum Equations with Physics-Informed Neural Networks: Boundary Conditions and Domain Decomposition

We investigate the application of Physics-Informed Neural Networks (PINNs) to the numerical solution of Einstein's vacuum field equations for static spacetimes. We first reproduce the Schwarzschild solution and then extend the method to the axisymmetric $q-$metric, a nontrivial exact solution characterized by a mass quadrupole moment. We analyze the influence of boundary conditions, domain decomposition, and equation redundancy on the convergence and stability of the training process. The proposed framework accurately reconstructs the metric functions in the computational domain while maintaining small residual errors. Our results demonstrate that PINNs provide a robust and flexible approach to solving Einstein's equations and offer a promising foundation for investigating gravitational configurations for which exact analytical solutions are unknown.