In paper:
chi(i) = mii
Q(4562) = \overline Q_2
Q(ij) = q_ij
Bb?(m1,m2,1,m3) are the primed functions in the paper
chi(456) is \overline\chi_1
The RC and RV notation has additional indices with sea quarks, remove those
  to get the notation of the paper
RD and RS correspond to R^z
Q(132) = Q(1)**2+Q(3)**2
Q(122) = Q(1)**2+Q(2)**2

P...xij means D_val=i, d_sea=j

lowest order is

P2x1221x11 = chi(1)+2 ZE F^2 (Q(1)-Q(2))^2;
P2x1221x21 = (chi(1)+chi(3)/2+2 ZE F^2 (Q(1)-Q(3))^2;

  P4x1221x11 =
       + F^-2*L4r * (  - 24*chi(1)*chi(4) )

       + F^-2*L5r * (  - 8*chi(1)^2 )

       + F^-2*L6r * ( 48*chi(1)*chi(4) )

       + F^-2*L8r * ( 16*chi(1)^2 )

       + ee^2*ZE*L4r * (  - 48*chi(4)*Q(12)^2 )

       + ee^2*ZE*L5r * (  - 16*chi(1)*Q(12)^2 )

       + ee^2*KE1r * (  - 12*chi(1)*Q(4562) )

       + ee^2*KE2r * (  - 12*chi(1)*Q(4562) )

       + ee^2*KE5r * (  - 4*chi(1)*Q(122) )

       + ee^2*KE6r * (  - 4*chi(1)*Q(122) )

       + ee^2*KE7r * ( 12*chi(1)*Q(4562) )

       + ee^2*KE8r * ( 12*chi(1)*Q(4562) + 12*chi(4)*Q(12)^2 )

       + ee^2*KE9r * ( 4*chi(1)*Q(1)^2 + 4*chi(1)*Q(2)^2 )

       + ee^2*KE10r * ( 4*chi(1)*Q(1)^2 + 4*chi(1)*Q(2)^2 + 8*chi(1)*Q(12)^2 )

       + ee^2*KE11r * ( 8*chi(1)*Q(12)^2 )

       + ee^2*KE18r * (  - 8*chi(1)*Q(1)*Q(2) )

       + ee^2*KE19r * (  - 4*chi(1)*Q(1)*Q(2) )

       + Ab(m11)*F^-2 * (  - 1/3*chi(1) )

       + Ab(m11)*ee^2 * ( Q(12)^2 )

       + Ab(m14)*ee^2*ZE * ( 6*Q(12)^2 )

       + Bb(mg2,m11,m11)*ee^2 * ( 4*chi(1)*Q(12)^2 )

       + Bb(m11,m11,0)*F^-2 * (  - 1/3*chi(1)*RD(1,4) )

       + B1b(mg2,m11,m11)*ee^2 * (  - 4*chi(1)*Q(12)^2 );

  P4x1221x12 =
       + F^-2*L4r * (  - 24*chi(1)*chi(456) )

       + F^-2*L5r * (  - 8*chi(1)^2 )

       + F^-2*L6r * ( 48*chi(1)*chi(456) )

       + F^-2*L8r * ( 16*chi(1)^2 )

       + ee^2*ZE*L4r * (  - 48*chi(456)*Q(12)^2 )

       + ee^2*ZE*L5r * (  - 16*chi(1)*Q(12)^2 )

       + ee^2*KE1r * (  - 12*chi(1)*Q(4562) )

       + ee^2*KE2r * (  - 12*chi(1)*Q(4562) )

       + ee^2*KE5r * (  - 4*chi(1)*Q(122) )

       + ee^2*KE6r * (  - 4*chi(1)*Q(122) )

       + ee^2*KE7r * ( 12*chi(1)*Q(4562) )

       + ee^2*KE8r * ( 12*chi(1)*Q(4562) + 12*chi(456)*Q(12)^2 )

       + ee^2*KE9r * ( 4*chi(1)*Q(1)^2 + 4*chi(1)*Q(2)^2 )

       + ee^2*KE10r * ( 4*chi(1)*Q(1)^2 + 4*chi(1)*Q(2)^2 + 8*chi(1)*Q(12)^2 )

       + ee^2*KE11r * ( 8*chi(1)*Q(12)^2 )

       + ee^2*KE18r * (  - 8*chi(1)*Q(1)*Q(2) )

       + ee^2*KE19r * (  - 4*chi(1)*Q(1)*Q(2) )

       + Ab(me2)*F^-2 * (  - 1/3*chi(1)*RS(me2,4,6,1,1) )

       + Ab(m11)*F^-2 * (  - 1/3*chi(1)*RC(1,4,6,me2,me2) )

       + Ab(m11)*ee^2 * ( Q(12)^2 )

       + Ab(m14)*ee^2*ZE * ( 4*Q(12)^2 )

       + Ab(m16)*ee^2*ZE * ( 2*Q(12)^2 )

       + Bb(mg2,m11,m11)*ee^2 * ( 4*chi(1)*Q(12)^2 )

       + Bb(m11,m11,0)*F^-2 * (  - 1/3*chi(1)*RD(1,4,6,me2) )

       + B1b(mg2,m11,m11)*ee^2 * (  - 4*chi(1)*Q(12)^2 );

   P4x1221x13 =
       + F^-2*L4r * (  - 24*chi(1)*chi(456) )

       + F^-2*L5r * (  - 8*chi(1)^2 )

       + F^-2*L6r * ( 48*chi(1)*chi(456) )

       + F^-2*L8r * ( 16*chi(1)^2 )

       + ee^2*ZE*L4r * (  - 48*chi(456)*Q(12)^2 )

       + ee^2*ZE*L5r * (  - 16*chi(1)*Q(12)^2 )

       + ee^2*KE1r * (  - 12*chi(1)*Q(4562) )

       + ee^2*KE2r * (  - 12*chi(1)*Q(4562) )

       + ee^2*KE5r * (  - 4*chi(1)*Q(122) )

       + ee^2*KE6r * (  - 4*chi(1)*Q(122) )

       + ee^2*KE7r * ( 12*chi(1)*Q(4562) )

       + ee^2*KE8r * ( 12*chi(1)*Q(4562) + 12*chi(456)*Q(12)^2 )

       + ee^2*KE9r * ( 4*chi(1)*Q(1)^2 + 4*chi(1)*Q(2)^2 )

       + ee^2*KE10r * ( 4*chi(1)*Q(1)^2 + 4*chi(1)*Q(2)^2 + 8*chi(1)*Q(12)^2 )

       + ee^2*KE11r * ( 8*chi(1)*Q(12)^2 )

       + ee^2*KE18r * (  - 8*chi(1)*Q(1)*Q(2) )

       + ee^2*KE19r * (  - 4*chi(1)*Q(1)*Q(2) )

       + Ab(mp2)*F^-2 * (  - 1/3*chi(1)*RS(mp2,4,5,6,1,1,me2) )

       + Ab(me2)*F^-2 * (  - 1/3*chi(1)*RS(me2,4,5,6,1,1,mp2) )

       + Ab(m11)*F^-2 * (  - 1/3*chi(1)*RC(1,4,5,6,mp2,me2) )

       + Ab(m11)*ee^2 * ( Q(12)^2 )

       + Ab(m14)*ee^2*ZE * ( 2*Q(12)^2 )

       + Ab(m15)*ee^2*ZE * ( 2*Q(12)^2 )

       + Ab(m16)*ee^2*ZE * ( 2*Q(12)^2 )

       + Bb(mg2,m11,m11)*ee^2 * ( 4*chi(1)*Q(12)^2 )

       + Bb(m11,m11,0)*F^-2 * (  - 1/3*chi(1)*RD(1,4,5,6,mp2,me2) )

       + B1b(mg2,m11,m11)*ee^2 * (  - 4*chi(1)*Q(12)^2 );

   P4x1331x21 =
       + F^-2*L4r * (  - 24*chi(4)*m13 )

       + F^-2*L5r * (  - 8*m13^2 )

       + F^-2*L6r * ( 48*chi(4)*m13 )

       + F^-2*L8r * ( 16*m13^2 )

       + ee^2*ZE*L4r * (  - 48*chi(4)*Q(13)^2 )

       + ee^2*ZE*L5r * (  - 16*Q(13)^2*m13 )

       + ee^2*KE1r * (  - 12*Q(4562)*m13 )

       + ee^2*KE2r * (  - 12*Q(4562)*m13 )

       + ee^2*KE5r * (  - 4*Q(132)*m13 )

       + ee^2*KE6r * (  - 4*Q(132)*m13 )

       + ee^2*KE7r * ( 12*Q(4562)*m13 )

       + ee^2*KE8r * ( 12*chi(4)*Q(13)^2 + 12*Q(4562)*m13 )

       + ee^2*KE9r * ( 4*chi(1)*Q(1)^2 + 4*chi(3)*Q(3)^2 )

       + ee^2*KE10r * ( 4*chi(1)*Q(1)^2 + 4*chi(3)*Q(3)^2 + 8*Q(13)^2*m13 )

       + ee^2*KE11r * ( 8*Q(13)^2*m13 )

       + ee^2*KE18r * (  - 8*Q(1)*Q(3)*m13 )

       + ee^2*KE19r * (  - 4*Q(1)*Q(3)*m13 )

       + Ab(m11)*F^-2 * (  - 1/3*RS(1,4,3)*m13 )

       + Ab(m13)*ee^2 * ( Q(13)^2 )

       + Ab(m14)*ee^2*ZE * ( 6*Q(1)*Q(13) )

       + Ab(m33)*F^-2 * (  - 1/3*RS(3,4,1)*m13 )

       + Ab(m34)*ee^2*ZE * (  - 6*Q(1)*Q(13) + 6*Q(13)^2 )

       + Bb(mg2,m13,m13)*ee^2 * ( 4*Q(13)^2*m13 )

       + B1b(mg2,m13,m13)*ee^2 * (  - 4*Q(13)^2*m13 );

   P4x1331x22 =
       + F^-2*L4r * (  - 24*chi(456)*m13 )

       + F^-2*L5r * (  - 8*m13^2 )

       + F^-2*L6r * ( 48*chi(456)*m13 )

       + F^-2*L8r * ( 16*m13^2 )

       + ee^2*ZE*L4r * (  - 48*chi(456)*Q(13)^2 )

       + ee^2*ZE*L5r * (  - 16*Q(13)^2*m13 )

       + ee^2*KE1r * (  - 12*Q(4562)*m13 )

       + ee^2*KE2r * (  - 12*Q(4562)*m13 )

       + ee^2*KE5r * (  - 4*Q(132)*m13 )

       + ee^2*KE6r * (  - 4*Q(132)*m13 )

       + ee^2*KE7r * ( 12*Q(4562)*m13 )

       + ee^2*KE8r * ( 12*chi(456)*Q(13)^2 + 12*Q(4562)*m13 )

       + ee^2*KE9r * ( 4*chi(1)*Q(1)^2 + 4*chi(3)*Q(3)^2 )

       + ee^2*KE10r * ( 4*chi(1)*Q(1)^2 + 4*chi(3)*Q(3)^2 + 8*Q(13)^2*m13 )

       + ee^2*KE11r * ( 8*Q(13)^2*m13 )

       + ee^2*KE18r * (  - 8*Q(1)*Q(3)*m13 )

       + ee^2*KE19r * (  - 4*Q(1)*Q(3)*m13 )

       + Ab(me2)*F^-2 * (  - 1/3*RS(me2,4,6,1,3)*m13 )

       + Ab(m11)*F^-2 * (  - 1/3*RS(1,4,6,3,me2)*m13 )

       + Ab(m13)*ee^2 * ( Q(13)^2 )

       + Ab(m14)*ee^2*ZE * ( 4*Q(1)*Q(13) - 2*Q(4)*Q(13) - 2*Q(5)*Q(13) )

       + Ab(m16)*ee^2*ZE * ( 2*Q(1)*Q(13) - 2*Q(6)*Q(13) )

       + Ab(m33)*F^-2 * (  - 1/3*RS(3,4,6,1,me2)*m13 )

       + Ab(m34)*ee^2*ZE * (  - 4*Q(1)*Q(13) + 2*Q(4)*Q(13) + 2*Q(5)*Q(13) + 4
         *Q(13)^2 )

       + Ab(m36)*ee^2*ZE * (  - 2*Q(1)*Q(13) + 2*Q(6)*Q(13) + 2*Q(13)^2 )

       + Bb(mg2,m13,m13)*ee^2 * ( 4*Q(13)^2*m13 )

       + B1b(mg2,m13,m13)*ee^2 * (  - 4*Q(13)^2*m13 );


   P4x1331x23 =
       + F^-2*L4r * (  - 24*chi(456)*m13 )

       + F^-2*L5r * (  - 8*m13^2 )

       + F^-2*L6r * ( 48*chi(456)*m13 )

       + F^-2*L8r * ( 16*m13^2 )

       + ee^2*ZE*L4r * (  - 48*chi(456)*Q(13)^2 )

       + ee^2*ZE*L5r * (  - 16*Q(13)^2*m13 )

       + ee^2*KE1r * (  - 12*Q(4562)*m13 )

       + ee^2*KE2r * (  - 12*Q(4562)*m13 )

       + ee^2*KE5r * (  - 4*Q(132)*m13 )

       + ee^2*KE6r * (  - 4*Q(132)*m13 )

       + ee^2*KE7r * ( 12*Q(4562)*m13 )

       + ee^2*KE8r * ( 12*chi(456)*Q(13)^2 + 12*Q(4562)*m13 )

       + ee^2*KE9r * ( 4*chi(1)*Q(1)^2 + 4*chi(3)*Q(3)^2 )

       + ee^2*KE10r * ( 4*chi(1)*Q(1)^2 + 4*chi(3)*Q(3)^2 + 8*Q(13)^2*m13 )

       + ee^2*KE11r * ( 8*Q(13)^2*m13 )

       + ee^2*KE18r * (  - 8*Q(1)*Q(3)*m13 )

       + ee^2*KE19r * (  - 4*Q(1)*Q(3)*m13 )

       + Ab(mp2)*F^-2 * (  - 1/3*RS(mp2,4,5,6,1,3,me2)*m13 )

       + Ab(me2)*F^-2 * (  - 1/3*RS(me2,4,5,6,1,3,mp2)*m13 )

       + Ab(m11)*F^-2 * (  - 1/3*RS(1,4,5,6,3,mp2,me2)*m13 )

       + Ab(m13)*ee^2 * ( Q(13)^2 )

       + Ab(m14)*ee^2*ZE * ( 2*Q(1)*Q(13) - 2*Q(4)*Q(13) )

       + Ab(m15)*ee^2*ZE * ( 2*Q(1)*Q(13) - 2*Q(5)*Q(13) )

       + Ab(m16)*ee^2*ZE * ( 2*Q(1)*Q(13) - 2*Q(6)*Q(13) )

       + Ab(m33)*F^-2 * (  - 1/3*RS(3,4,5,6,1,mp2,me2)*m13 )

       + Ab(m34)*ee^2*ZE * (  - 2*Q(1)*Q(13) + 2*Q(4)*Q(13) + 2*Q(13)^2 )

       + Ab(m35)*ee^2*ZE * (  - 2*Q(1)*Q(13) + 2*Q(5)*Q(13) + 2*Q(13)^2 )

       + Ab(m36)*ee^2*ZE * (  - 2*Q(1)*Q(13) + 2*Q(6)*Q(13) + 2*Q(13)^2 )

       + Bb(mg2,m13,m13)*ee^2 * ( 4*Q(13)^2*m13 )

       + B1b(mg2,m13,m13)*ee^2 * (  - 4*Q(13)^2*m13 );