Lowest order is one
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

Everything normalized to F_0,
i.e. lowest order is 1

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

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

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

       + ee^2*KE1r * ( 6*Q(4562) )

       + ee^2*KE2r * ( 6*Q(4562) )

       + ee^2*KE5r * ( 2*Q(1)^2 + 2*Q(2)^2 )

       + ee^2*KE6r * ( 2*Q(1)^2 + 2*Q(2)^2 )

       + ee^2*KE12r * ( 2*Q(12)^2 )

       + ee^2*KE18r * ( 4*Q(1)*Q(2) )

       + ee^2*KE19r * ( 2*Q(1)*Q(2) )

       + Ab(me14)*F^-2 * ( 1/4 )

       + Ab(me15)*F^-2 * ( 1/4 )

       + Ab(me16)*F^-2 * ( 1/4 )

       + Ab(me24)*F^-2 * ( 1/4 )

       + Ab(me25)*F^-2 * ( 1/4 )

       + Ab(me26)*F^-2 * ( 1/4 )

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

       + B1b(mg2,m11,m11)*ee^2 * (  - Q(12)^2 )

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

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

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

       + ee^2*KE1r * ( 6*Q(4562) )

       + ee^2*KE2r * ( 6*Q(4562) )

       + ee^2*KE5r * ( 2*Q(122) )

       + ee^2*KE6r * ( 2*Q(122) )

       + ee^2*KE12r * ( 2*Q(12)^2 )

       + ee^2*KE18r * ( 4*Q(1)*Q(2) )

       + ee^2*KE19r * ( 2*Q(1)*Q(2) )

       + Ab(me14)*F^-2 * ( 1/4 )

       + Ab(me15)*F^-2 * ( 1/4 )

       + Ab(me16)*F^-2 * ( 1/4 )

       + Ab(me24)*F^-2 * ( 1/4 )

       + Ab(me25)*F^-2 * ( 1/4 )

       + Ab(me26)*F^-2 * ( 1/4 )

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

       + B1b(mg2,m11,m11)*ee^2 * (  - Q(12)^2 )

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


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

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

       + ee^2*KE1r * ( 6*Q(4562) )

       + ee^2*KE2r * ( 6*Q(4562) )

       + ee^2*KE5r * ( 2*Q(122) )

       + ee^2*KE6r * ( 2*Q(122) )

       + ee^2*KE12r * ( 2*Q(12)^2 )

       + ee^2*KE18r * ( 4*Q(1)*Q(2) )

       + ee^2*KE19r * ( 2*Q(1)*Q(2) )

       + Ab(me14)*F^-2 * ( 1/4 )

       + Ab(me15)*F^-2 * ( 1/4 )

       + Ab(me16)*F^-2 * ( 1/4 )

       + Ab(me24)*F^-2 * ( 1/4 )

       + Ab(me25)*F^-2 * ( 1/4 )

       + Ab(me26)*F^-2 * ( 1/4 )

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

       + B1b(mg2,m11,m11)*ee^2 * (  - Q(12)^2 )

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

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

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

       + ee^2*KE1r * ( 6*Q(4562) )

       + ee^2*KE2r * ( 6*Q(4562) )

       + ee^2*KE5r * ( 2*Q(132) )

       + ee^2*KE6r * ( 2*Q(132) )

       + ee^2*KE12r * ( 2*Q(13)^2 )

       + ee^2*KE18r * ( 4*Q(1)*Q(3) )

       + ee^2*KE19r * ( 2*Q(1)*Q(3) )

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

       + Ab(me14)*F^-2 * ( 1/4 )

       + Ab(me15)*F^-2 * ( 1/4 )

       + Ab(me16)*F^-2 * ( 1/4 )

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

       + Ab(me34)*F^-2 * ( 1/4 )

       + Ab(me35)*F^-2 * ( 1/4 )

       + Ab(me36)*F^-2 * ( 1/4 )

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

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

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

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

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

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

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

       + ee^2*KE1r * ( 6*Q(4562) )

       + ee^2*KE2r * ( 6*Q(4562) )

       + ee^2*KE5r * ( 2*Q(132) )

       + ee^2*KE6r * ( 2*Q(132) )

       + ee^2*KE12r * ( 2*Q(13)^2 )

       + ee^2*KE18r * ( 4*Q(1)*Q(3) )

       + ee^2*KE19r * ( 2*Q(1)*Q(3) )

       + Ab(me2)*F^-2 * (  - 1/12*RV(me2,1,3) )

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

       + Ab(me14)*F^-2 * ( 1/4 )

       + Ab(me15)*F^-2 * ( 1/4 )

       + Ab(me16)*F^-2 * ( 1/4 )

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

       + Ab(me34)*F^-2 * ( 1/4 )

       + Ab(me35)*F^-2 * ( 1/4 )

       + Ab(me36)*F^-2 * ( 1/4 )

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

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

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

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

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


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

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

       + ee^2*KE1r * ( 6*Q(4562) )

       + ee^2*KE2r * ( 6*Q(4562) )

       + ee^2*KE5r * ( 2*Q(132) )

       + ee^2*KE6r * ( 2*Q(132) )

       + ee^2*KE12r * ( 2*Q(13)^2 )

       + ee^2*KE18r * ( 4*Q(1)*Q(3) )

       + ee^2*KE19r * ( 2*Q(1)*Q(3) )

       + Ab(mp2)*F^-2 * (  - 1/12*RV(mp2,me2,1,3) )

       + Ab(me2)*F^-2 * (  - 1/12*RV(me2,mp2,1,3) )

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

       + Ab(me14)*F^-2 * ( 1/4 )

       + Ab(me15)*F^-2 * ( 1/4 )

       + Ab(me16)*F^-2 * ( 1/4 )

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

       + Ab(me34)*F^-2 * ( 1/4 )

       + Ab(me35)*F^-2 * ( 1/4 )

       + Ab(me36)*F^-2 * ( 1/4 )

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

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

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

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

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