H g ( E ) the corresponding expansion to the order ξ 2 it is given by: For the magnetic combination H g ( E ) the expansion to the order ξ 2 reads Modelling H g ( E $^{)}$( x, ξ, t = 0 ) On fig 1 we show the results of the numerical computation of G ( E ) 0 contribution into the electric combination of gluon GPDs H g ( E ) at t = 0 for several different values of ξ . As the numerical input for forward gluon distributions we used MRST LO fit [36] for nucleon parton distributions at Q 2 = 1 GeV $^{2}$. 14 H g ( E $^{)}$( x, ξ, t ) = 5 12 xG ( E ) 0 ( x, t ) − 1 6 x 2 ∂ ∂x G ( E ) 0 ( x, t ) + 1 8 ∫ 1 x dy ( x y ) 3 2 G ( E ) 0 ( y, t ) + ξ 2 [ 1 32 ( − 7 3 x + x ) G ( E ) 0 ( x, t ) + 1 24 (1 + x $^{2}$) ∂ ∂x G ( E ) 0 ( x, t ) + 1 24 x (1 − x $^{2}$) ∂ 2 ∂x $_{2}$G ( E ) 0 ( x, t ) + ∫ 1 x dyG ( E ) 0 ( y, t ) ( 5 128 ( x y ) 1 2 + 3 128 ( x y ) 3 2 + 1 y 2 ( 3 128 ( y x ) 1 2 + 5 128 ( x y ) 1 2 )) + 7 48 xG ( E ) 2 ( x, t ) − 1 24 x 2 ∂ ∂x G ( E ) 2 ( x, t )+ ∫ 1 x dyG ( E ) 2 ( y, t ) ( 35 256 ( y x ) 1 2 + 5 128 ( x y ) 1 2 + 3 256 ( x y ) 3 2 )] + O ( ξ $^{4}$) . (38) H g ( M $^{)}$( x, ξ, t ) = 1 8 xG ( M ) 0 ( x, t ) − 1 4 x 2 ∂ ∂x G ( M ) 0 ( x, t ) + 1 6 x 3 ∂ 2 ∂x $_{2}$G ( M ) 0 ( x, t ) − 1 16 ∫ 1 x dy ( x y ) 3 2 G ( M ) 0 ( y, t ) + ξ 2 [ 1 16 x G ( M ) 0 ( x, t ) + 1 32 (1 − 5 x $^{2}$) ∂ ∂x G ( M ) 0 ( x, t ) − 1 8 x (1 − x $^{2}$) ∂ 2 ∂x $_{2}$G ( M ) 0 ( x, t ) − 1 24 x $^{2}$(1 − x $^{2}$) ∂ 3 ∂x $_{3}$G ( M ) 0 ( x, t ) − ∫ 1 x dyG ( M ) 0 ( y, t ) ( 15 256 ( x y ) 1 2 + 15 256 ( x y ) 3 2 + 1 y 2 ( 3 256 ( y x ) 1 2 + 15 256 ( x y ) 1 2 )) − 5 48 xG ( M ) 2 ( x, t ) − 1 48 x 2 ∂ ∂x G ( M ) 2 ( x, t ) + 1 24 x 3 ∂ 2 ∂x $_{2}$G ( M ) 2 ( x, t ) − ∫ 1 x dyG ( M ) 2 ( y, t ) ( 35 512 ( y x ) 1 2 + 15 256 ( x y ) 1 2 + 15 512 ( x y ) 3 2 )] + O ( ξ $^{4}$) . (39)