Bifurcating solutions of Rotating-Bowen-York initial data with a positive cosmological constant

A generalization of the Bowen–York initial data to the case with the positive cosmological constant is investigated numerically.  We follow the construction presented recently by Bizoń, Pletka and Simon, and solve numerically the corresponding Lichnerowicz equation on a compactified domain S1 × S2 . We find and describe new solutions, bifurcating from those discovered by Bizoń et al . We provide numerical arguments suggesting the absence of additional branches of solutions.