### abstract

- A new approach to global QCD analysis is developed. The main ingredients are two QCD-based evolution equations. The first one is the Balitsky-Kovchegov non-linear equation, which sums higher twists while preserving unitarity. The second equation is linear and it is responsible for the correct short distance behavior of the theory, namely it includes the DGLAP kernel. Our approach allows for extrapolation of the parton distributions to the very high energies available at the LHC as well as very low photon virtualities, \(Q^2\ll 1{\mathrm{GeV^2}}\). All existing low x data on the F2 structure function are reproduced using two fitting parameters, the other parameters were taken as constants. The result is \(\chi^2/df = 1\). Analyzing the parameter \(\lambda\equiv \partial\ln F_2/\partial(\ln 1/x)\) at very low x and Q2 well below \(1 {\mathrm{GeV^2}}\) we find \(\lambda \simeq 0.08 \)-0.1. This is a result which agrees with the “soft pomeron” intercept without involving soft physics.