Black hole initial data in Gauss-Bonnet gravity: Momentarily static case

We study the method for generating the initial data of black hole systems in Gauss-Bonnet gravity. The initial data are assumed to be momentarily static and conformally flat. Although the equation for the conformal factor is highly nonlinear, it is successfully solved by numerical relaxation for one-black-hole and two-black-hole systems. The common apparent horizon is studied in the two-black-hole initial data, and the result suggests that the Penrose inequalities are satisfied in this system. This is the first step for simulating black hole collisions in higher-curvature theories.