Homotopy approach for scattering amplitude for running QCD coupling

In this paper we proposed the homotopy approach for solving the nonlinear Balitsky-Kovchegov (BK) evolution equation with running QCD coupling. The approach consists of two steps. First, is the analytic solution to the nonlinear evolution equation for the simplified, leading twist kernel. Second, is the iteration procedure that allow us to calculate corrections analytically or semi-numerically. For the leading twist kernel it is shown that the first iteration leads to $\leq 1\%$ accuracy. The $\zeta = -\frac{4N_c}{b_0}Y \ln\Lb \bas (1/Q^2_s(Y))/\bas\Lb r^2\Rb\Rb$ ($r$ is the dipole size, $Q_s$ is the saturation scale) and geometric scaling behaviour of the scattering amplitude are discussed as well as the dependence on the value of the infrared cutoff.