gfactors#

gfactors = {'PPR': (0, 0, sqrt(5)/(5*sqrt(2*J_i + 1)*Abs(2*J_i - 1))), 'PQQ': (0, sqrt(3)*(J_i - 1)/(6*J_i*(2*J_i - 1)*sqrt(2*J_i + 1)), -sqrt(5)*(J_i + 1)/(10*J_i*sqrt(2*J_i + 1)*Abs(2*J_i - 1))), 'PQR': (0, -sqrt(3)*sqrt((J_i - 1)*(2*J_i - 1))*sqrt(J_i + 1)/(6*J_i*(2*J_i - 1)*(2*J_i + 1)), -sqrt(5)*sqrt(J_i - 1)*sqrt(J_i + 1)/(10*J_i*sqrt(2*J_i - 1)*(2*J_i + 1))), 'PRP': (1/(3*(2*J_i + 1)**(3/2)), -sqrt(3)*(J_i + 1)/(6*J_i*(2*J_i + 1)**(3/2)), sqrt(5)*(J_i + 1)*(2*J_i + 3)/(30*J_i*(2*J_i - 1)*(2*J_i + 1)**(3/2))), 'PRQ': (-1/(3*(2*J_i + 1)**(3/2)), sqrt(3)/(6*J_i*(2*J_i + 1)**(3/2)), sqrt(5)*(2*J_i + 3)/(30*J_i*(2*J_i + 1)**(3/2))), 'PRR': (1/(3*(2*J_i + 1)**(3/2)), sqrt(3)/(6*(2*J_i + 1)**(3/2)), sqrt(5)/(30*(2*J_i + 1)**(3/2))), 'QPQ': (0, -sqrt(3)*sqrt(J_i - 1)*sqrt(J_i + 1)/(6*J_i*sqrt(2*J_i - 1)*(2*J_i + 1)), -sqrt(5)*sqrt(J_i - 1)*sqrt(J_i + 1)/(10*J_i*sqrt(2*J_i - 1)*(2*J_i + 1))), 'QPR': (0, sqrt(3)*(J_i + 1)/(6*J_i*(2*J_i + 1)**(3/2)), -sqrt(5)*Abs(J_i - 1)/(10*J_i*(2*J_i + 1)**(3/2))), 'QQP': (-1/(3*(2*J_i + 1)**(3/2)), sqrt(3)/(6*J_i*(2*J_i + 1)**(3/2)), sqrt(5)*(2*J_i + 3)/(30*J_i*(2*J_i + 1)**(3/2))), 'QQQ': (1/(3*(2*J_i + 1)**(3/2)), -sqrt(3)/(6*J_i*(J_i + 1)*(2*J_i + 1)**(3/2)), sqrt(5)*(2*J_i - 1)*(2*J_i + 3)/(30*J_i*(J_i + 1)*(2*J_i + 1)**(3/2))), 'QQR': (-1/(3*(2*J_i + 1)**(3/2)), -sqrt(3)/((2*J_i + 1)**(3/2)*(6*J_i + 6)), sqrt(5)*Abs(2*J_i - 1)/((2*J_i + 1)**(3/2)*(30*J_i + 30))), 'QRP': (0, sqrt(3)*J_i/((2*J_i + 1)**(3/2)*(6*J_i + 6)), -sqrt(5)*(J_i + 2)/((2*J_i + 1)**(3/2)*(10*J_i + 10))), 'QRQ': (0, -sqrt(3)*sqrt(J_i)*sqrt(J_i + 2)/((2*J_i + 1)*sqrt(2*J_i + 3)*(6*J_i + 6)), -sqrt(5)*sqrt(J_i)*sqrt(J_i + 2)/((2*J_i + 1)*sqrt(2*J_i + 3)*(10*J_i + 10))), 'RPP': (1/(3*(2*J_i + 1)**(3/2)), sqrt(3)/(6*(2*J_i + 1)**(3/2)), sqrt(5)/(30*(2*J_i + 1)**(3/2))), 'RPQ': (-1/(3*(2*J_i + 1)**(3/2)), -sqrt(3)/((2*J_i + 1)**(3/2)*(6*J_i + 6)), sqrt(5)*Abs(2*J_i - 1)/((2*J_i + 1)**(3/2)*(30*J_i + 30))), 'RPR': (1/(3*(2*J_i + 1)**(3/2)), -sqrt(3)*J_i/((2*J_i + 1)**(3/2)*(6*J_i + 6)), sqrt(5)*J_i*(2*J_i - 1)/((2*J_i + 1)**(3/2)*(2*J_i + 3)*(30*J_i + 30))), 'RQP': (0, -sqrt(3)*sqrt(J_i)*sqrt(J_i + 2)/((2*J_i + 1)*sqrt(2*J_i + 3)*(6*J_i + 6)), -sqrt(5)*sqrt(J_i)*sqrt(J_i + 2)/((2*J_i + 1)*sqrt(2*J_i + 3)*(10*J_i + 10))), 'RQQ': (0, sqrt(3)*(J_i + 2)/(sqrt(2*J_i + 1)*(2*J_i + 3)*(6*J_i + 6)), -sqrt(5)*J_i/(sqrt(2*J_i + 1)*(2*J_i + 3)*(10*J_i + 10))), 'RRP': (0, 0, sqrt(5)/(5*sqrt(2*J_i + 1)*(2*J_i + 3)))}#

Isotropic version of angular momentum-dependent combinations of three Wigner-6j factors involved in evaluating four-fold dipole interaction operator.