################################ # Diagram information ################################ name = "npltrb" edges = [[1, 2], [2, 3], [3, 5], [4, 5], [1, 5], [2, 4]] nodes = [1, 1, 3, 4] internal_masses = [m2, m2, m2, m2, m2, m2] external_masses = [M2, M2, M2, M2] U = x[1]*x[2] + x[1]*x[3] + x[1]*x[4] + x[1]*x[6] + x[2]*x[4] + x[2]*x[5] + x[2]*x[6] + x[3]*x[4] + x[3]*x[5] + x[3]*x[6] + x[4]*x[5] + x[5]*x[6] F = -m2*x[1]^2*x[2] - m2*x[1]^2*x[3] - m2*x[1]^2*x[4] - m2*x[1]^2*x[6] - m2*x[1]*x[2]^2 + (M2 - 2*m2)*x[1]*x[2]*x[3] + (M2 - 3*m2)*x[1]*x[2]*x[4] + (-2*m2 + s)*x[1]*x[2]*x[5] - 3*m2*x[1]*x[2]*x[6] - m2*x[1]*x[3]^2 + (-3*m2 + s)*x[1]*x[3]*x[4] + (-2*m2 + s)*x[1]*x[3]*x[5] + (M2 - 3*m2)*x[1]*x[3]*x[6] - m2*x[1]*x[4]^2 + (-2*m2 + s)*x[1]*x[4]*x[5] + (M2 - 2*m2)*x[1]*x[4]*x[6] + (-2*m2 + s)*x[1]*x[5]*x[6] - m2*x[1]*x[6]^2 - m2*x[2]^2*x[4] - m2*x[2]^2*x[5] - m2*x[2]^2*x[6] + (M2 - 2*m2)*x[2]*x[3]*x[4] + (M2 - 2*m2)*x[2]*x[3]*x[5] + (M2 - 2*m2)*x[2]*x[3]*x[6] - m2*x[2]*x[4]^2 + (M2 - 3*m2)*x[2]*x[4]*x[5] + (M2 - 2*m2)*x[2]*x[4]*x[6] - m2*x[2]*x[5]^2 + (-3*m2 + s)*x[2]*x[5]*x[6] - m2*x[2]*x[6]^2 - m2*x[3]^2*x[4] - m2*x[3]^2*x[5] - m2*x[3]^2*x[6] - m2*x[3]*x[4]^2 - 3*m2*x[3]*x[4]*x[5] + (M2 - 2*m2)*x[3]*x[4]*x[6] - m2*x[3]*x[5]^2 + (M2 - 3*m2)*x[3]*x[5]*x[6] - m2*x[3]*x[6]^2 - m2*x[4]^2*x[5] - m2*x[4]*x[5]^2 + (M2 - 2*m2)*x[4]*x[5]*x[6] - m2*x[5]^2*x[6] - m2*x[5]*x[6]^2 parameters = [M2, m2, s] variables = [x[1], x[2], x[3], x[4], x[5], x[6]] χ_generic = 116 f_vector = [36, 120, 169, 125, 51, 11] ################################ # Component 1 ################################ D[1] = M2 χ[1] = 52 weights[1] = [[0, -1, -1, 1, 1, 1], [0, 1, 1, -1, 1, -1], [1, -1, -1, 0, 1, 1], [1, -1, -1, 1, 0, 1], [1, -1, -1, 1, 1, 0], [1, 0, 1, -1, 1, -1], [1, 1, 0, -1, 1, -1], [1, 1, 1, -1, 0, -1], [-1, -1, 0, -1, 0, 0], [-1, 0, -1, 0, 0, -1], [0, -1, -1, 0, 1, 1], [0, -1, -1, 1, 0, 1], [0, -1, -1, 1, 1, 0], [0, 0, 1, -1, 1, -1], [0, 1, 0, -1, 1, -1], [0, 1, 1, -1, 0, -1], [0, -1, 0, -1, -1, 0], [1, -1, -1, 0, 0, 1], [1, -1, -1, 0, 1, 0], [1, -1, -1, 1, 0, 0], [1, 0, 0, -1, 1, -1], [1, 0, 1, -1, 0, -1], [0, 0, -1, 0, -1, -1], [1, 1, 0, -1, 0, -1], [-1, -1, -1, 0, 0, -1], [0, 0, 0, 1, 1, 1], [-1, -1, -1, 0, -1, 0], [-1, -1, 0, -1, 0, -1], [0, 0, 1, 0, 1, 1], [0, 1, 0, 1, 1, 0], [0, 1, 1, 0, 1, 0], [-1, 0, 0, -1, -1, -1], [0, -1, -1, 0, 0, 1], [0, -1, -1, 0, 1, 0], [0, -1, -1, 1, 0, 0], [0, 0, 0, -1, 1, -1], [0, 0, 1, -1, 0, -1], [0, 1, 0, -1, 0, -1], [1, -1, -1, 0, 0, 0], [1, 0, 1, 0, 0, 1], [0, -1, -1, -1, -1, 0], [1, 0, 1, 0, 1, 0], [1, 0, 0, 0, 1, 1], [1, 0, 0, 1, 0, 1], [1, 0, 0, 1, 1, 0], [1, 0, 0, -1, 0, -1], [0, 0, -1, -1, -1, -1], [1, 1, 0, 0, 1, 0], [1, 1, 0, 1, 0, 0], [1, 1, 1, 0, 0, 0], [0, 0, 0, 1, 1, 0], [0, 0, 0, 1, 0, 1], [0, 0, 1, 0, 1, 0], [0, 1, 1, 0, 0, 0], [0, -1, -1, 0, 0, 0], [0, 0, 0, -1, 0, -1], [1, 0, 0, 0, 0, 1], [1, 0, 0, 0, 1, 0], [1, 1, 0, 0, 0, 0]] computed_with[1] = ["PLD_sym", "PLD_num"] ################################ # Component 2 ################################ D[2] = M2 - 1//4*s χ[2] = 70 weights[2] = [[-1, -1, -1, -1, 0, 0], [-1, -1, 0, -1, -1, 0], [-1, 0, -1, -1, 0, -1], [-1, 0, -1, 0, -1, -1], [0, -1, 0, -1, -1, -1], [0, -1, -1, 0, -1, -1], [-1, -1, -1, -1, 0, -1], [0, 0, 0, 0, 1, 1], [-1, -1, -1, -1, -1, 0], [-1, -1, -1, 0, -1, -1], [-1, -1, 0, -1, -1, -1], [0, 0, 1, 0, 0, 1], [0, 1, 0, 0, 1, 0], [-1, 0, -1, -1, -1, -1], [0, 1, 0, 1, 0, 0], [1, 0, 1, 0, 0, 0], [0, -1, -1, -1, -1, -1], [1, 0, 0, 1, 0, 0], [-1, -1, -1, -1, -1, -1], [0, 0, 0, 0, 1, 0], [0, 0, 0, 0, 0, 1], [0, 0, 0, 1, 0, 0], [0, 0, 1, 0, 0, 0], [0, 1, 0, 0, 0, 0], [1, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0]] computed_with[2] = ["PLD_num"] ################################ # Component 3 ################################ D[3] = M2 - 4*m2 χ[3] = 89 weights[3] = [[0, -1, -1, 1, 1, 1], [0, 1, 1, -1, 1, -1], [1, -1, -1, 0, 1, 1], [1, -1, -1, 1, 0, 1], [1, -1, -1, 1, 1, 0], [1, 0, 1, -1, 1, -1], [1, 1, 0, -1, 1, -1], [1, 1, 1, -1, 0, -1], [0, -1, -1, 0, 1, 1], [0, -1, -1, 1, 0, 1], [0, -1, -1, 1, 1, 0], [0, 0, 1, -1, 1, -1], [0, 1, 0, -1, 1, -1], [0, 1, 1, -1, 0, -1], [1, -1, -1, 0, 0, 1], [1, -1, -1, 0, 1, 0], [1, -1, -1, 1, 0, 0], [1, 0, 0, -1, 1, -1], [1, 0, 1, -1, 0, -1], [1, 1, 0, -1, 0, -1], [0, -1, -1, 0, 0, 1], [0, -1, -1, 0, 1, 0], [0, -1, -1, 1, 0, 0], [0, 0, 0, -1, 1, -1], [0, 0, 1, -1, 0, -1], [0, 1, 0, -1, 0, -1], [1, -1, -1, 0, 0, 0], [0, -1, -1, -1, 0, -1], [1, 0, 0, -1, 0, -1], [0, -1, -1, 0, 0, 0], [0, 0, 0, -1, 0, -1]] computed_with[3] = ["PLD_sym", "PLD_num"] ################################ # Component 4 ################################ D[4] = M2 - 9*m2 χ[4] = 112 weights[4] = [[-1, -1, 0, -1, 0, 0], [-1, 0, -1, 0, 0, -1], [0, -1, 0, -1, -1, 0], [0, 0, -1, 0, -1, -1]] computed_with[4] = ["PLD_sym"] ################################ # Component 5 ################################ D[5] = M2 - m2 χ[5] = 92 weights[5] = [[-1, -1, 0, -1, 0, 0], [-1, 0, -1, 0, 0, -1], [0, -1, 0, -1, -1, 0], [0, 0, -1, 0, -1, -1], [-1, -1, 0, -1, -1, 0], [-1, 0, -1, 0, -1, -1]] computed_with[5] = ["PLD_sym", "PLD_num"] ################################ # Component 6 ################################ D[6] = M2 - s χ[6] = 102 weights[6] = [[-1, -1, -1, -1, 0, -1], [-1, -1, -1, -1, -1, 0], [-1, -1, -1, 0, -1, -1], [-1, -1, 0, -1, -1, -1], [-1, 0, -1, -1, -1, -1], [0, -1, -1, -1, -1, -1], [-1, -1, -1, -1, -1, -1]] computed_with[6] = ["PLD_num"] ################################ # Component 7 ################################ D[7] = M2^2 - 10*M2*m2 + 9*m2^2 + 4*m2*s χ[7] = 114 weights[7] = [[-1, -1, 0, -1, -1, 0], [-1, 0, -1, 0, -1, -1]] computed_with[7] = ["PLD_num"] ################################ # Component 8 ################################ D[8] = M2^2*s + 9*M2*m2^2 - 7*M2*m2*s + m2*s^2 χ[8] = 112 weights[8] = [[-1, -1, -1, -1, 0, 0], [-1, 0, -1, -1, 0, -1], [0, -1, 0, -1, -1, -1], [0, -1, -1, 0, -1, -1]] computed_with[8] = ["PLD_num"] ################################ # Component 9 ################################ D[9] = M2^4 - 6*M2^3*m2 + 9*M2^2*m2^2 + M2^2*m2*s - 6*M2*m2^2*s + m2^2*s^2 χ[9] = 114 weights[9] = [[-1, -1, -1, -1, 0, -1], [0, -1, -1, -1, -1, -1]] computed_with[9] = ["PLD_num"] ################################ # Component 10 ################################ D[10] = m2 χ[10] = 28 weights[10] = [[-1, 0, 1, 1, 0, 1], [-1, 1, 0, 1, 0, 1], [-1, 1, 1, 0, 0, 1], [-1, 1, 1, 1, 0, 0], [0, -1, 0, 1, 1, 1], [0, 0, -1, 1, 1, 1], [0, 1, 1, -1, 1, 0], [0, 1, 1, 0, 1, -1], [1, -1, 0, 0, 1, 1], [1, -1, 0, 1, 0, 1], [1, -1, 0, 1, 1, 0], [1, 0, 1, -1, 1, 0], [0, 0, 1, 1, -1, 1], [1, 0, 1, 0, 1, -1], [1, 0, -1, 0, 1, 1], [1, 0, -1, 1, 0, 1], [1, 0, -1, 1, 1, 0], [1, 1, 0, -1, 1, 0], [0, 1, 0, 1, -1, 1], [1, 1, 0, 0, 1, -1], [1, 1, 1, -1, 0, 0], [0, 1, 1, 0, -1, 1], [0, 1, 1, 1, -1, 0], [1, 1, 1, 0, 0, -1], [-1, 0, 0, 1, 0, 1], [-1, 0, 1, 0, 0, 1], [-1, 0, 1, 1, 0, 0], [-1, -1, 0, 0, 0, 0], [-1, 1, 0, 0, 0, 1], [-1, 1, 0, 1, 0, 0], [-1, 0, -1, 0, 0, 0], [-1, 1, 1, 0, 0, 0], [-1, 0, 0, -1, 0, 0], [-1, 0, 0, 0, 0, -1], [0, -1, 0, 0, 1, 1], [0, -1, 0, 1, 0, 1], [0, -1, 0, 1, 1, 0], [0, 0, -1, 0, 1, 1], [0, 0, -1, 1, 0, 1], [0, 0, -1, 1, 1, 0], [0, 0, 1, -1, 1, 0], [0, 1, 0, -1, 1, 0], [0, 1, 1, -1, 0, 0], [0, 0, 1, 0, 1, -1], [0, 1, 0, 0, 1, -1], [0, 1, 1, 0, 0, -1], [1, -1, 0, 0, 0, 1], [1, -1, 0, 0, 1, 0], [0, -1, 0, -1, 0, 0], [1, -1, 0, 1, 0, 0], [0, -1, 0, 0, -1, 0], [0, -1, 0, 0, 0, -1], [1, 0, 0, -1, 1, 0], [1, 0, 1, -1, 0, 0], [0, 0, 0, 1, -1, 1], [0, 0, 1, 0, -1, 1], [0, 0, 1, 1, -1, 0], [1, 0, 0, 0, 1, -1], [1, 0, 1, 0, 0, -1], [1, 0, -1, 0, 0, 1], [1, 0, -1, 0, 1, 0], [0, 0, -1, -1, 0, 0], [1, 0, -1, 1, 0, 0], [0, 0, -1, 0, -1, 0], [0, 0, -1, 0, 0, -1], [1, 1, 0, -1, 0, 0], [0, 1, 0, 0, -1, 1], [0, 1, 0, 1, -1, 0], [1, 1, 0, 0, 0, -1], [0, 0, 0, -1, -1, 0], [0, 1, 1, 0, -1, 0], [0, 0, 0, 0, -1, -1], [-1, 0, 0, 0, 0, 1], [-1, 0, 0, 1, 0, 0], [-1, 0, 1, 0, 0, 0], [-1, -1, 0, 0, 0, -1], [-1, 1, 0, 0, 0, 0], [0, -1, 0, 0, 0, 1], [0, -1, 0, 0, 1, 0], [0, -1, 0, 1, 0, 0], [0, 0, -1, 0, 0, 1], [0, 0, -1, 0, 1, 0], [0, 0, -1, 1, 0, 0], [0, 0, 0, -1, 1, 0], [0, 0, 1, -1, 0, 0], [0, 1, 0, -1, 0, 0], [0, 0, 0, 0, 1, -1], [0, 0, 1, 0, 0, -1], [0, 1, 0, 0, 0, -1], [1, -1, 0, 0, 0, 0], [1, 0, 0, -1, 0, 0], [0, 0, 0, 0, -1, 1], [0, 0, 0, 1, -1, 0], [0, 0, 1, 0, -1, 0], [1, 0, 0, 0, 0, -1], [1, 0, -1, 0, 0, 0], [0, 0, -1, -1, -1, 0], [0, 1, 0, 0, -1, 0], [-1, 0, 0, 0, 0, 0], [0, -1, 0, 0, 0, 0], [0, 0, -1, 0, 0, 0], [0, 0, 0, -1, 0, 0], [0, 0, 0, 0, 0, -1], [0, 0, 0, 0, -1, 0], [-1, -1, -1, -1, -1, -1]] computed_with[10] = ["PLD_sym", "PLD_num"] ################################ # Component 11 ################################ D[11] = m2 - 1//4*s χ[11] = 101 weights[11] = [[-1, 0, 1, 1, -1, 1], [-1, 1, 0, 1, -1, 1], [-1, 1, 1, 0, -1, 1], [-1, 1, 1, 1, -1, 0], [-1, 0, 0, 1, -1, 1], [-1, 0, 1, 0, -1, 1], [-1, 0, 1, 1, -1, 0], [-1, 1, 0, 0, -1, 1], [-1, 1, 0, 1, -1, 0], [-1, 1, 1, 0, -1, 0], [-1, 0, 0, 0, -1, 1], [-1, 0, 0, 1, -1, 0], [-1, 0, 1, 0, -1, 0], [-1, 1, 0, 0, -1, 0], [-1, 0, 0, 0, -1, 0]] computed_with[11] = ["PLD_sym", "PLD_num"] ################################ # Component 12 ################################ D[12] = m2 - 1//9*s χ[12] = 114 weights[12] = [[-1, 0, -1, -1, 0, 0], [0, -1, 0, 0, -1, -1]] computed_with[12] = ["PLD_sym"] ################################ # Component 13 ################################ D[13] = m2 - s χ[13] = 110 weights[13] = [[-1, 0, -1, -1, 0, 0], [-1, 0, -1, -1, 0, 0], [-1, 0, -1, -1, 0, 0], [0, -1, 0, 0, -1, -1], [0, -1, 0, 0, -1, -1], [0, -1, 0, 0, -1, -1]] computed_with[13] = ["PLD_sym"] ################################ # Component 14 ################################ D[14] = s χ[14] = 30 weights[14] = [[-1, 0, 1, 1, -1, 1], [-1, 1, 0, 1, -1, 1], [-1, 1, 1, 0, -1, 1], [-1, 1, 1, 1, -1, 0], [-1, 0, 0, 1, -1, 1], [-1, 0, 1, 0, -1, 1], [-1, 0, 1, 1, -1, 0], [-1, 0, -1, -1, 0, 0], [-1, 1, 0, 0, -1, 1], [-1, 1, 0, 1, -1, 0], [-1, 1, 1, 0, -1, 0], [0, -1, 0, 0, -1, -1], [-1, 0, 0, 0, -1, 1], [-1, 0, 0, 1, -1, 0], [-1, -1, -1, 0, -1, 0], [-1, 0, 1, 0, -1, 0], [-1, -1, 0, 0, -1, -1], [0, 0, 1, 1, 0, 1], [-1, 1, 0, 0, -1, 0], [0, 1, 0, 0, 1, 1], [-1, 0, -1, -1, -1, 0], [0, 1, 0, 1, 0, 1], [-1, 0, 0, -1, -1, -1], [0, 1, 1, 0, 0, 1], [0, 1, 1, 1, 0, 0], [1, 0, 1, 1, 0, 0], [-1, -1, -1, -1, 0, -1], [-1, 0, 0, 0, -1, 0], [0, 0, 0, 1, 0, 1], [0, 0, 1, 1, 0, 0], [0, 1, 0, 0, 0, 1], [0, 1, 1, 0, 0, 0], [0, -1, -1, -1, -1, -1], [-1, -1, -1, -1, -1, -1], [0, 0, 0, 0, 1, 0], [1, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0]] computed_with[14] = ["PLD_sym", "PLD_num"]