These were calculated in "Finite Volume at Two-loops in Chiral Perturbation Theory" J. Bijnens, T. Rössler LU TP 14-38, arXiv:1411.6384, JHEP 1501 (2015) 034 Expressions in terms of physical masses and physical Fpi (infinite volume) +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ Pion mass: massp4V = + AbV(mp2)*Fpi^-2 * ( - 1/2*mp2 ) + AbV(me2)*Fpi^-2 * ( 1/6*mp2 ); massp6LV = + AbV(mp2)*F^-4 * ( 40*L8r*mp2^2 + 80*L6r*mp2^2 - 24*L5r*mp2^2 - 48*L4r *mp2^2 + 28*L3r*mp2^2 + 32*L2r*mp2^2 + 56*L1r*mp2^2 ) + AbV(mk2)*F^-4 * ( 32*L8r*mp2*mk2 + 64*L6r*mp2*mk2 - 16*L5r*mp2*mk2 - 64*L4r*mp2*mk2 + 20*L3r*mp2*mk2 + 16*L2r*mp2*mk2 + 64*L1r*mp2*mk2 ) + AbV(me2)*F^-4 * ( 8*L8r*mp2^2 - 64/3*L7r*mp2*mk2 + 64/3*L7r*mp2^2 + 64/3*L6r*mp2*mk2 - 16/3*L6r*mp2^2 - 32/9*L5r*mp2*mk2 + 8/9*L5r*mp2^2 - 64/3*L4r*mp2*mk2 + 16/3*L4r*mp2^2 + 16/3*L3r*mp2*mk2 - 4/3*L3r* mp2^2 + 16/3*L2r*mp2*mk2 - 4/3*L2r*mp2^2 + 64/3*L1r*mp2*mk2 - 16/3* L1r*mp2^2 ) + A23bV(mp2)*F^-4 * ( - 12*L3r*mp2 - 48*L2r*mp2 - 24*L1r*mp2 ) + A23bV(mk2)*F^-4 * ( - 12*L3r*mp2 - 48*L2r*mp2 ) + A23bV(me2)*F^-4 * ( - 4*L3r*mp2 - 12*L2r*mp2 ); massp6RV = + AbV(mp2)*F^-4 * ( pi16*mp2*mk2 + pi16*mp2^2 - 2*Ab(mp2)*mp2 - Ab(mk2) *mp2 + 1/4*B0b(mp2,mp2)*mp2^2 + 1/12*B0b(me2,me2)*mp2^2 ) + AbV(mp2)^2*F^-4 * ( - 3/8*mp2 ) + AbV(mp2)*AbV(mk2)*F^-4 * ( - 1/2*mp2 ) + AbV(mp2)*AbV(me2)*F^-4 * ( - 1/12*mp2 ) + AbV(mp2)*B0bV(mp2,mp2)*F^-4 * ( 1/4*mp2^2 ) + AbV(mp2)*B0bV(me2,me2)*F^-4 * ( 1/12*mp2^2 ) + AbV(mk2)*F^-4 * ( pi16*mp2*mk2 - 1/2*Ab(mp2)*mp2 - 1/2*Ab(mk2)*mp2 - 1/2*Ab(me2)*mp2 - 2/9*B0b(me2,me2)*mp2*mk2 ) + AbV(mk2)^2*F^-4 * ( - 1/4*mp2 ) + AbV(mk2)*AbV(me2)*F^-4 * ( - 1/2*mp2 ) + AbV(mk2)*B0bV(me2,me2)*F^-4 * ( - 2/9*mp2*mk2 ) + AbV(me2)*F^-4 * ( 1/3*Ab(mp2)*mp2 - 1/3*Ab(mk2)*mp2 - 1/12*B0b(mp2, mp2)*mp2^2 + 4/27*B0b(me2,me2)*mp2*mk2 - 7/108*B0b(me2,me2)*mp2^2 ) + AbV(me2)^2*F^-4 * ( 1/72*mp2 ) + AbV(me2)*B0bV(mp2,mp2)*F^-4 * ( - 1/12*mp2^2 ) + AbV(me2)*B0bV(me2,me2)*F^-4 * ( 4/27*mp2*mk2 - 7/108*mp2^2 ) + HHbV(mp2,mp2,mp2,mp2)*F^-4 * ( 5/6*mp2^2 ) + HHbV(mp2,mk2,mk2,mp2)*F^-4 * ( mp2*mk2 - 5/8*mp2^2 ) + HHbV(mp2,me2,me2,mp2)*F^-4 * ( 1/18*mp2^2 ) + HHbV(mk2,mk2,me2,mp2)*F^-4 * ( 1/2*mp2*mk2 + 1/24*mp2^2 ) + HH1bV(mp2,mk2,mk2,mp2)*F^-4 * ( mp2^2 ) + HH1bV(me2,mk2,mk2,mp2)*F^-4 * ( - mp2^2 ) + HH21bV(mp2,mp2,mp2,mp2)*F^-4 * ( 3*mp2^2 ) + HH21bV(mp2,mk2,mk2,mp2)*F^-4 * ( - 3/8*mp2^2 ) + HH21bV(mk2,mp2,mk2,mp2)*F^-4 * ( 3*mp2^2 ) + HH21bV(me2,mk2,mk2,mp2)*F^-4 * ( 9/8*mp2^2 ) + HH27bV(mp2,mp2,mp2,mp2)*F^-4 * ( - 3*mp2 ) + HH27bV(mp2,mk2,mk2,mp2)*F^-4 * ( 3/8*mp2 ) + HH27bV(mk2,mp2,mk2,mp2)*F^-4 * ( - 3*mp2 ) + HH27bV(me2,mk2,mk2,mp2)*F^-4 * ( - 9/8*mp2 ); ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ Kaon mass massp4V = + AbV(me2)*Fpi^-2 * ( - 1/4*me2 - 1/12*mp2 ); massp6LV = + AbV(mp2)*F^-4 * ( 24*L8r*mp2*mk2 + 48*L6r*mp2*mk2 - 12*L5r*mp2*mk2 - 48*L4r*mp2*mk2 + 15*L3r*mp2*mk2 + 12*L2r*mp2*mk2 + 48*L1r*mp2*mk2 ) + AbV(mk2)*F^-4 * ( 48*L8r*mk2^2 + 96*L6r*mk2^2 - 24*L5r*mk2^2 - 64*L4r *mk2^2 + 30*L3r*mk2^2 + 36*L2r*mk2^2 + 72*L1r*mk2^2 ) + AbV(me2)*F^-4 * ( 64/3*L8r*mk2^2 - 56/3*L8r*mp2*mk2 + 16/3*L8r*mp2^2 + 64/3*L7r*mk2^2 - 32*L7r*mp2*mk2 + 32/3*L7r*mp2^2 + 64/3*L6r*mk2^2 - 16/3*L6r*mp2*mk2 - 64/9*L5r*mk2^2 + 4/3*L5r*mp2*mk2 - 8/9*L5r* mp2^2 - 64/3*L4r*mk2^2 + 16/3*L4r*mp2*mk2 + 28/3*L3r*mk2^2 - 7/3*L3r* mp2*mk2 + 16/3*L2r*mk2^2 - 4/3*L2r*mp2*mk2 + 64/3*L1r*mk2^2 - 16/3* L1r*mp2*mk2 ) + A23bV(mp2)*F^-4 * ( - 9*L3r*mk2 - 36*L2r*mk2 ) + A23bV(mk2)*F^-4 * ( - 18*L3r*mk2 - 60*L2r*mk2 - 24*L1r*mk2 ) + A23bV(me2)*F^-4 * ( - L3r*mk2 - 12*L2r*mk2 ); massp6RV = + AbV(mp2)*F^-4 * ( 3/4*pi16*mk2^2 - 3/16*Ab(mp2)*mk2 - 3/4*Ab(mk2)*mk2 - 3/16*Ab(me2)*mk2 - 1/6*B0b(me2,me2)*mp2*mk2 ) + AbV(mp2)^2*F^-4 * ( - 3/32*mk2 ) + AbV(mp2)*AbV(mk2)*F^-4 * ( - 3/4*mk2 ) + AbV(mp2)*AbV(me2)*F^-4 * ( - 3/16*mk2 ) + AbV(mp2)*B0bV(me2,me2)*F^-4 * ( - 1/6*mp2*mk2 ) + AbV(mk2)*F^-4 * ( 3/4*pi16*mk2^2 + 3/4*pi16*mp2*mk2 - 3/4*Ab(mp2)*mk2 - 3/2*Ab(mk2)*mk2 + 4/9*B0b(me2,me2)*mk2^2 ) + AbV(mk2)^2*F^-4 * ( - 3/4*mk2 ) + AbV(mk2)*B0bV(me2,me2)*F^-4 * ( 4/9*mk2^2 ) + AbV(me2)*F^-4 * ( 1/2*pi16*mk2^2 + 1/4*pi16*mp2*mk2 - 41/48*Ab(mp2)* mk2 + 1/12*Ab(mp2)*mp2 - 2/3*Ab(mk2)*mk2 + 19/48*Ab(me2)*mk2 - 1/12* Ab(me2)*mp2 - 8/27*B0b(me2,me2)*mk2^2 + 7/54*B0b(me2,me2)*mp2*mk2 ) + AbV(me2)^2*F^-4 * ( 25/288*mk2 ) + AbV(me2)*B0bV(me2,me2)*F^-4 * ( - 8/27*mk2^2 + 7/54*mp2*mk2 ) + HHbV(mp2,mp2,mk2,mk2)*F^-4 * ( - 15/32*mk2^2 + 3/4*mp2*mk2 ) + HHbV(mp2,mk2,me2,mk2)*F^-4 * ( 13/16*mk2^2 ) + HHbV(mk2,mk2,mk2,mk2)*F^-4 * ( 3/4*mk2^2 ) + HHbV(mk2,me2,me2,mk2)*F^-4 * ( 181/288*mk2^2 ) + HH1bV(mk2,mp2,mp2,mk2)*F^-4 * ( 3/4*mk2^2 ) + HH1bV(mk2,mp2,me2,mk2)*F^-4 * ( - 3/2*mk2^2 ) + HH1bV(mk2,me2,me2,mk2)*F^-4 * ( - 5/4*mk2^2 ) + HH21bV(mp2,mp2,mk2,mk2)*F^-4 * ( 9/4*mk2^2 ) + HH21bV(mk2,mp2,mp2,mk2)*F^-4 * ( - 9/32*mk2^2 ) + HH21bV(mk2,mp2,me2,mk2)*F^-4 * ( 27/16*mk2^2 ) + HH21bV(mk2,mk2,mk2,mk2)*F^-4 * ( 9/4*mk2^2 ) + HH21bV(mk2,me2,me2,mk2)*F^-4 * ( 27/32*mk2^2 ) + HH27bV(mp2,mp2,mk2,mk2)*F^-4 * ( - 9/4*mk2 ) + HH27bV(mk2,mp2,mp2,mk2)*F^-4 * ( 9/32*mk2 ) + HH27bV(mk2,mp2,me2,mk2)*F^-4 * ( - 27/16*mk2 ) + HH27bV(mk2,mk2,mk2,mk2)*F^-4 * ( - 9/4*mk2 ) + HH27bV(mk2,me2,me2,mk2)*F^-4 * ( - 27/32*mk2 ); +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ Eta mass: massp4V = + AbV(mp2)*Fpi^-2 * ( 1/2*mp2 ) + AbV(mk2)*Fpi^-2 * ( - me2 - 1/3*mp2 ) + AbV(me2)*Fpi^-2 * ( 8/9*mk2 - 7/18*mp2 ); massp6LV = + AbV(mp2)*F^-4 * ( 24*L8r*mp2^2 - 64*L7r*mp2*mk2 + 64*L7r*mp2^2 + 64* L6r*mp2*mk2 - 16*L6r*mp2^2 - 32/3*L5r*mp2*mk2 + 8/3*L5r*mp2^2 - 64* L4r*mp2*mk2 + 16*L4r*mp2^2 + 16*L3r*mp2*mk2 - 4*L3r*mp2^2 + 16*L2r* mp2*mk2 - 4*L2r*mp2^2 + 64*L1r*mp2*mk2 - 16*L1r*mp2^2 ) + AbV(mk2)*F^-4 * ( 256/3*L8r*mk2^2 - 224/3*L8r*mp2*mk2 + 64/3*L8r* mp2^2 + 256/3*L7r*mk2^2 - 128*L7r*mp2*mk2 + 128/3*L7r*mp2^2 + 256/3* L6r*mk2^2 - 64/3*L6r*mp2*mk2 - 256/9*L5r*mk2^2 + 16/3*L5r*mp2*mk2 - 32/9*L5r*mp2^2 - 256/3*L4r*mk2^2 + 64/3*L4r*mp2*mk2 + 112/3*L3r*mk2^2 - 28/3*L3r*mp2*mk2 + 64/3*L2r*mk2^2 - 16/3*L2r*mp2*mk2 + 256/3*L1r* mk2^2 - 64/3*L1r*mp2*mk2 ) + AbV(me2)*F^-4 * ( 896/9*L8r*mk2^2 - 1024/9*L8r*mp2*mk2 + 344/9*L8r* mp2^2 + 1024/9*L7r*mk2^2 - 1664/9*L7r*mp2*mk2 + 640/9*L7r*mp2^2 + 256/ 3*L6r*mk2^2 - 128/3*L6r*mp2*mk2 + 16/3*L6r*mp2^2 - 832/27*L5r*mk2^2 + 896/27*L5r*mp2*mk2 - 280/27*L5r*mp2^2 - 256/9*L4r*mk2^2 + 128/9* L4r*mp2*mk2 - 16/9*L4r*mp2^2 + 64/3*L3r*mk2^2 - 32/3*L3r*mp2*mk2 + 4/ 3*L3r*mp2^2 + 128/3*L2r*mk2^2 - 64/3*L2r*mp2*mk2 + 8/3*L2r*mp2^2 + 128/3*L1r*mk2^2 - 64/3*L1r*mp2*mk2 + 8/3*L1r*mp2^2 ) + A23bV(mp2)*F^-4 * ( - 16*L3r*mk2 + 4*L3r*mp2 - 48*L2r*mk2 + 12*L2r* mp2 ) + A23bV(mk2)*F^-4 * ( - 16/3*L3r*mk2 + 4/3*L3r*mp2 - 64*L2r*mk2 + 16* L2r*mp2 ) + A23bV(me2)*F^-4 * ( - 16*L3r*mk2 + 4*L3r*mp2 - 32*L2r*mk2 + 8*L2r* mp2 - 32*L1r*mk2 + 8*L1r*mp2 ); massp6RV = + AbV(mp2)*F^-4 * ( Ab(mp2)*mp2 - Ab(mk2)*mp2 - 1/4*B0b(mp2,mp2)*mp2^2 + 4/9*B0b(me2,me2)*mp2*mk2 - 7/36*B0b(me2,me2)*mp2^2 ) + AbV(mp2)^2*F^-4 * ( - 1/8*mp2 ) + AbV(mp2)*AbV(mk2)*F^-4 * ( - 3/2*mp2 ) + AbV(mp2)*AbV(me2)*F^-4 * ( 1/12*mp2 ) + AbV(mp2)*B0bV(mp2,mp2)*F^-4 * ( - 1/4*mp2^2 ) + AbV(mp2)*B0bV(me2,me2)*F^-4 * ( 4/9*mp2*mk2 - 7/36*mp2^2 ) + AbV(mk2)*F^-4 * ( 8/3*pi16*mk2^2 + 2/3*pi16*mp2*mk2 - 1/3*pi16*mp2^2 - 8/3*Ab(mp2)*mk2 - 7/6*Ab(mp2)*mp2 - 2/3*Ab(mk2)*mk2 + 3/2*Ab(mk2)* mp2 - 8/3*Ab(me2)*mk2 + 7/6*Ab(me2)*mp2 - 32/27*B0b(me2,me2)*mk2^2 + 14/27*B0b(me2,me2)*mp2*mk2 ) + AbV(mk2)^2*F^-4 * ( mk2 + 3/4*mp2 ) + AbV(mk2)*AbV(me2)*F^-4 * ( - 32/9*mk2 + 3/2*mp2 ) + AbV(mk2)*B0bV(me2,me2)*F^-4 * ( - 32/27*mk2^2 + 14/27*mp2*mk2 ) + AbV(me2)*F^-4 * ( 16/9*Ab(mp2)*mk2 - 8/9*Ab(mp2)*mp2 - 8/3*Ab(mk2)* mk2 + 10/9*Ab(mk2)*mp2 + 8/9*Ab(me2)*mk2 - 2/9*Ab(me2)*mp2 + 1/12* B0b(mp2,mp2)*mp2^2 + 4/9*B0b(mk2,mk2)*mk2^2 + 64/81*B0b(me2,me2)* mk2^2 - 56/81*B0b(me2,me2)*mp2*mk2 + 49/324*B0b(me2,me2)*mp2^2 ) + AbV(me2)^2*F^-4 * ( 8/27*mk2 - 31/216*mp2 ) + AbV(me2)*B0bV(mp2,mp2)*F^-4 * ( 1/12*mp2^2 ) + AbV(me2)*B0bV(mk2,mk2)*F^-4 * ( 4/9*mk2^2 ) + AbV(me2)*B0bV(me2,me2)*F^-4 * ( 64/81*mk2^2 - 56/81*mp2*mk2 + 49/324* mp2^2 ) + HHbV(mp2,mp2,me2,me2)*F^-4 * ( 1/6*mp2^2 ) + HHbV(mp2,mk2,mk2,me2)*F^-4 * ( 3/2*mp2*mk2 + 1/8*mp2^2 ) + HHbV(mk2,mk2,me2,me2)*F^-4 * ( 50/9*mk2^2 - 11/3*mp2*mk2 + 5/8*mp2^2 ) + HHbV(me2,me2,me2,me2)*F^-4 * ( 128/243*mk2^2 - 112/243*mp2*mk2 + 49/ 486*mp2^2 ) + HH1bV(mp2,mk2,mk2,me2)*F^-4 * ( - 4*mp2*mk2 + mp2^2 ) + HH1bV(me2,mk2,mk2,me2)*F^-4 * ( - 32/3*mk2^2 + 20/3*mp2*mk2 - mp2^2 ) + HH21bV(mp2,mk2,mk2,me2)*F^-4 * ( 6*mk2^2 - 3*mp2*mk2 + 3/8*mp2^2 ) + HH21bV(me2,mk2,mk2,me2)*F^-4 * ( 6*mk2^2 - 3*mp2*mk2 + 3/8*mp2^2 ) + HH27bV(mp2,mk2,mk2,me2)*F^-4 * ( - 9/2*mk2 + 9/8*mp2 ) + HH27bV(me2,mk2,mk2,me2)*F^-4 * ( - 9/2*mk2 + 9/8*mp2 ); ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ Pion decay constant: decayp4V = + AbV(mp2)*Fpi^-1 * ( 1 ) + AbV(mk2)*Fpi^-1 * ( 1/2 ); decayp6LV = + AbV(mp2)*F^-3 * ( 6*L5r*mp2 + 12*L4r*mp2 - 14*L3r*mp2 - 16*L2r*mp2 - 28*L1r*mp2 ) + AbV(mk2)*F^-3 * ( 4*L5r*mp2 + 16*L4r*mk2 - 10*L3r*mk2 - 8*L2r*mk2 - 32*L1r*mk2 ) + AbV(me2)*F^-3 * ( 2/3*L5r*mp2 + 16/3*L4r*mk2 - 4/3*L4r*mp2 - 8/3*L3r* mk2 + 2/3*L3r*mp2 - 8/3*L2r*mk2 + 2/3*L2r*mp2 - 32/3*L1r*mk2 + 8/3* L1r*mp2 ) + A23bV(mp2)*F^-3 * ( 6*L3r + 24*L2r + 12*L1r ) + A23bV(mk2)*F^-3 * ( 6*L3r + 24*L2r ) + A23bV(me2)*F^-3 * ( 2*L3r + 6*L2r ); decayp6RV = + AbV(mp2)*F^-3 * ( - 1/2*pi16*mk2 - 1/4*pi16*mp2 + 1/2*Ab(mp2) + 1/2* Ab(mk2) ) + AbV(mp2)*B0bV(mp2,mp2)*F^-3 * ( - 1/2*mp2 ) + AbV(mk2)*F^-3 * ( - 1/2*pi16*mk2 - 1/8*pi16*mp2 + 1/2*Ab(mp2) + 1/4* Ab(mk2) ) + AbV(me2)*F^-3 * ( 1/6*pi16*mk2 - 1/6*pi16*mp2 + 1/6*Ab(mp2) - 1/6*Ab( mk2) ) + AbV(me2)*B0bV(mp2,mp2)*F^-3 * ( 1/6*mp2 ) + AbV(me2)*B0bV(mk2,mk2)*F^-3 * ( - 1/6*mk2 ) + HHbV(mp2,mp2,mp2,mp2)*F^-3 * ( - 1/2*mp2 ) + HHbV(mp2,mk2,mk2,mp2)*F^-3 * ( - 1/2*mk2 + 1/16*mp2 ) + HHbV(mk2,mk2,me2,mp2)*F^-3 * ( - 1/4*mk2 + 1/16*mp2 ) + HH27bV(mp2,mp2,mp2,mp2)*F^-3 * ( 3/2 ) + HH27bV(mp2,mk2,mk2,mp2)*F^-3 * ( - 3/16 ) + HH27bV(mk2,mp2,mk2,mp2)*F^-3 * ( 3/2 ) + HH27bV(me2,mk2,mk2,mp2)*F^-3 * ( 9/16 ) + HHdV(mp2,mp2,mp2,mp2)*F^-3 * ( 5/12*mp2^2 ) + HHdV(mp2,mk2,mk2,mp2)*F^-3 * ( 1/2*mp2*mk2 - 5/16*mp2^2 ) + HHdV(mp2,me2,me2,mp2)*F^-3 * ( 1/36*mp2^2 ) + HHdV(mk2,mk2,me2,mp2)*F^-3 * ( 1/4*mp2*mk2 + 1/48*mp2^2 ) + HH1dV(mp2,mk2,mk2,mp2)*F^-3 * ( 1/2*mp2^2 ) + HH1dV(me2,mk2,mk2,mp2)*F^-3 * ( - 1/2*mp2^2 ) + HH21dV(mp2,mp2,mp2,mp2)*F^-3 * ( 3/2*mp2^2 ) + HH21dV(mp2,mk2,mk2,mp2)*F^-3 * ( - 3/16*mp2^2 ) + HH21dV(mk2,mp2,mk2,mp2)*F^-3 * ( 3/2*mp2^2 ) + HH21dV(me2,mk2,mk2,mp2)*F^-3 * ( 9/16*mp2^2 ) + HH27dV(mp2,mp2,mp2,mp2)*F^-3 * ( - 3/2*mp2 ) + HH27dV(mp2,mk2,mk2,mp2)*F^-3 * ( 3/16*mp2 ) + HH27dV(mk2,mp2,mk2,mp2)*F^-3 * ( - 3/2*mp2 ) + HH27dV(me2,mk2,mk2,mp2)*F^-3 * ( - 9/16*mp2 ); +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ Kaon decay constant: decayp4V = + AbV(mp2)*Fpi^-1 * ( 3/8 ) + AbV(mk2)*Fpi^-1 * ( 3/4 ) + AbV(me2)*Fpi^-1 * ( 3/8 ); decayp6LV = + AbV(mp2)*F^-3 * ( 3/2*L5r*mk2 + 3/2*L5r*mp2 + 12*L4r*mp2 - 15/2*L3r* mp2 - 6*L2r*mp2 - 24*L1r*mp2 ) + AbV(mk2)*F^-3 * ( 3*L5r*mk2 + 3*L5r*mp2 + 16*L4r*mk2 - 15*L3r*mk2 - 18*L2r*mk2 - 36*L1r*mk2 ) + AbV(me2)*F^-3 * ( 1/6*L5r*mk2 + 3/2*L5r*mp2 + 16/3*L4r*mk2 - 4/3*L4r* mp2 - 14/3*L3r*mk2 + 7/6*L3r*mp2 - 8/3*L2r*mk2 + 2/3*L2r*mp2 - 32/3* L1r*mk2 + 8/3*L1r*mp2 ) + A23bV(mp2)*F^-3 * ( 9/2*L3r + 18*L2r ) + A23bV(mk2)*F^-3 * ( 9*L3r + 30*L2r + 12*L1r ) + A23bV(me2)*F^-3 * ( 1/2*L3r + 6*L2r ); decayp6RV = + AbV(mp2)*F^-3 * ( - 15/32*pi16*mk2 + 3/16*pi16*mp2 - 3/64*Ab(mp2) + 9/32*Ab(mk2) + 9/64*Ab(me2) + 3/16*B0b(me2,me2)*mp2 ) + AbV(mp2)^2*F^-3 * ( - 15/128 ) + AbV(mp2)*AbV(mk2)*F^-3 * ( 3/32 ) + AbV(mp2)*AbV(me2)*F^-3 * ( 9/64 ) + AbV(mp2)*B0bV(mp2,mp2)*F^-3 * ( - 3/16*mp2 ) + AbV(mp2)*B0bV(me2,me2)*F^-3 * ( 3/16*mp2 ) + AbV(mk2)*F^-3 * ( - 9/16*pi16*mk2 - 3/8*pi16*mp2 + 27/32*Ab(mp2) + 9/ 16*Ab(mk2) - 9/32*Ab(me2) - 1/2*B0b(me2,me2)*mk2 ) + AbV(mk2)^2*F^-3 * ( 3/32 ) + AbV(mk2)*AbV(me2)*F^-3 * ( - 9/32 ) + AbV(mk2)*B0bV(me2,me2)*F^-3 * ( - 1/2*mk2 ) + AbV(me2)*F^-3 * ( - 3/32*pi16*mk2 - 3/16*pi16*mp2 + 37/64*Ab(mp2) - 11/32*Ab(mk2) + 9/64*Ab(me2) + 1/3*B0b(me2,me2)*mk2 - 7/48*B0b(me2, me2)*mp2 ) + AbV(me2)^2*F^-3 * ( 9/128 ) + AbV(me2)*B0bV(mp2,mp2)*F^-3 * ( 1/16*mp2 ) + AbV(me2)*B0bV(mk2,mk2)*F^-3 * ( - 1/4*mk2 ) + AbV(me2)*B0bV(me2,me2)*F^-3 * ( 1/3*mk2 - 7/48*mp2 ) + HHbV(mp2,mp2,mk2,mk2)*F^-3 * ( 3/64*mk2 - 3/8*mp2 ) + HHbV(mp2,mk2,me2,mk2)*F^-3 * ( - 9/32*mk2 ) + HHbV(mk2,mk2,mk2,mk2)*F^-3 * ( - 3/8*mk2 ) + HHbV(mk2,me2,me2,mk2)*F^-3 * ( - 9/64*mk2 ) + HH27bV(mp2,mp2,mk2,mk2)*F^-3 * ( 9/8 ) + HH27bV(mk2,mp2,mp2,mk2)*F^-3 * ( - 9/64 ) + HH27bV(mk2,mp2,me2,mk2)*F^-3 * ( 27/32 ) + HH27bV(mk2,mk2,mk2,mk2)*F^-3 * ( 9/8 ) + HH27bV(mk2,me2,me2,mk2)*F^-3 * ( 27/64 ) + HHdV(mp2,mp2,mk2,mk2)*F^-3 * ( - 15/64*mk2^2 + 3/8*mp2*mk2 ) + HHdV(mp2,mk2,me2,mk2)*F^-3 * ( 13/32*mk2^2 ) + HHdV(mk2,mk2,mk2,mk2)*F^-3 * ( 3/8*mk2^2 ) + HHdV(mk2,me2,me2,mk2)*F^-3 * ( 181/576*mk2^2 ) + HH1dV(mk2,mp2,mp2,mk2)*F^-3 * ( 3/8*mk2^2 ) + HH1dV(mk2,mp2,me2,mk2)*F^-3 * ( - 3/4*mk2^2 ) + HH1dV(mk2,me2,me2,mk2)*F^-3 * ( - 5/8*mk2^2 ) + HH21dV(mp2,mp2,mk2,mk2)*F^-3 * ( 9/8*mk2^2 ) + HH21dV(mk2,mp2,mp2,mk2)*F^-3 * ( - 9/64*mk2^2 ) + HH21dV(mk2,mp2,me2,mk2)*F^-3 * ( 27/32*mk2^2 ) + HH21dV(mk2,mk2,mk2,mk2)*F^-3 * ( 9/8*mk2^2 ) + HH21dV(mk2,me2,me2,mk2)*F^-3 * ( 27/64*mk2^2 ) + HH27dV(mp2,mp2,mk2,mk2)*F^-3 * ( - 9/8*mk2 ) + HH27dV(mk2,mp2,mp2,mk2)*F^-3 * ( 9/64*mk2 ) + HH27dV(mk2,mp2,me2,mk2)*F^-3 * ( - 27/32*mk2 ) + HH27dV(mk2,mk2,mk2,mk2)*F^-3 * ( - 9/8*mk2 ) + HH27dV(mk2,me2,me2,mk2)*F^-3 * ( - 27/64*mk2 ); ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ Eta decay constant: decayp4V = + AbV(mk2)*Fpi^-1 * ( 3/2 ); decayp6LV = + AbV(mp2)*F^-3 * ( 2*L5r*mp2 + 12*L4r*mp2 - 6*L3r*mp2 - 6*L2r*mp2 - 24 *L1r*mp2 ) + AbV(mk2)*F^-3 * ( 8/3*L5r*mk2 + 4*L5r*mp2 + 16*L4r*mk2 - 14*L3r*mk2 - 8*L2r*mk2 - 32*L1r*mk2 ) + AbV(me2)*F^-3 * ( 32/9*L5r*mk2 - 14/9*L5r*mp2 + 16/3*L4r*mk2 - 4/3* L4r*mp2 - 8*L3r*mk2 + 2*L3r*mp2 - 16*L2r*mk2 + 4*L2r*mp2 - 16*L1r*mk2 + 4*L1r*mp2 ) + A23bV(mp2)*F^-3 * ( 6*L3r + 18*L2r ) + A23bV(mk2)*F^-3 * ( 2*L3r + 24*L2r ) + A23bV(me2)*F^-3 * ( 6*L3r + 12*L2r + 12*L1r ); decayp6RV = + AbV(mk2)*F^-3 * ( - 3/2*pi16*mk2 - 3/8*pi16*mp2 + 3/2*Ab(mp2) + 3/4* Ab(mk2) ) + AbV(me2)*F^-3 * ( 1/2*pi16*mk2 - 1/2*Ab(mk2) ) + AbV(me2)*B0bV(mk2,mk2)*F^-3 * ( - 1/2*mk2 ) + HHbV(mp2,mk2,mk2,me2)*F^-3 * ( - 9/16*mp2 ) + HHbV(mk2,mk2,me2,me2)*F^-3 * ( - 3/4*mk2 + 3/16*mp2 ) + HH27bV(mp2,mk2,mk2,me2)*F^-3 * ( 27/16 ) + HH27bV(me2,mk2,mk2,me2)*F^-3 * ( 27/16 ) + HHdV(mp2,mp2,me2,me2)*F^-3 * ( 1/12*mp2^2 ) + HHdV(mp2,mk2,mk2,me2)*F^-3 * ( 3/4*mp2*mk2 + 1/16*mp2^2 ) + HHdV(mk2,mk2,me2,me2)*F^-3 * ( 25/9*mk2^2 - 11/6*mp2*mk2 + 5/16*mp2^2 ) + HHdV(me2,me2,me2,me2)*F^-3 * ( 64/243*mk2^2 - 56/243*mp2*mk2 + 49/972 *mp2^2 ) + HH1dV(mp2,mk2,mk2,me2)*F^-3 * ( - 2*mp2*mk2 + 1/2*mp2^2 ) + HH1dV(me2,mk2,mk2,me2)*F^-3 * ( - 16/3*mk2^2 + 10/3*mp2*mk2 - 1/2* mp2^2 ) + HH21dV(mp2,mk2,mk2,me2)*F^-3 * ( 3*mk2^2 - 3/2*mp2*mk2 + 3/16*mp2^2 ) + HH21dV(me2,mk2,mk2,me2)*F^-3 * ( 3*mk2^2 - 3/2*mp2*mk2 + 3/16*mp2^2 ) + HH27dV(mp2,mk2,mk2,me2)*F^-3 * ( - 9/4*mk2 + 9/16*mp2 ) + HH27dV(me2,mk2,mk2,me2)*F^-3 * ( - 9/4*mk2 + 9/16*mp2 );