## Short thesis description.

he thesis is an application of the Lax pair and the Riemann-Hilbert problem (RHP) techniques to an old non linear boundary partial differential equation (PDE) problem originated in Einstein's theory of General Relativity.

The Lax Pair technique is roughly the following. Consider a pair of the linear differential equations

where the is unknown matrix valued function and are known regular matrix valued functions. If a regular solution exists the we must have

or

We compute both left and right hand sides and obtain

The last expression is a non-linear PDE if and depend on some unknown function of interest. For example, we may have

where the function is the unknown function of interest and the functions and have known functional form. Hence, we may regard a non-linear PDE as a solvability condition for some linear system of differential equations.

How this helps? Suppose we introduce a complex spectral parameter into the matrices and : , . Consequently, the solution would also depend on the parameter . For the particular functions , one regularly can prove that the first non trivial term of asymptotic of the solution depends on the . Hence, if we construct a general solution of the linear PDE system

then we recover a class of solutions of the non-linear PDE .

It is an elementary fact from complex analysis that if a function of a complex argument has known singularities and it is otherwise analytic then the information about the singularities determines the function. The task of recovering of such a function is called the Riemann-Hilbert problem.

The solution exists on the -Riemann surface of the functions , . The spectral analysis of the system reveals necessary singularities of the solution. A researcher may be able to add freedom into the problem by discovering a way to introduce more singularities. This way the problem gives rise to a RHP on some restriction of the Riemann surface of the functions , . The construction of a solution of such RHP is a developed area. The technique is called the Backer-Akhiezer function.

If this all has been performed then there are still two remaining questions. First, the solution is specified in terms of the RHP parameters while the original PDE comes from physics. One has to identify the physical meaning of the parameters. Second, the general procedure outlined above provides no statements about boundary conditions. One has to find a particular structure of the RHP which produces the solution with needed boundary conditions.

In the presented thesis all the steps described above have been performed for the following physical problem. Two black holes are rotating around the common axis of symmetry. The situation is stationary. The goal is to find a relativistic correction to the Newton gravity law that is due to the rotation of the bodies.

The axial and time symmetry of the problem allows reduction of the four dimensional space time to a two dimensional non linear PDE

where the and are unknown functions and some singular boundary conditions are imposed on the axis of symmetry between the rotating black holes as well as some conditions of decay at the -infinity.