################################ # 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 = [M[1], M[2], M[3], M[4]] 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 + (M[3] - 2*m2)*x[1]*x[2]*x[3] + (M[4] - 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] + (M[3] - 3*m2)*x[1]*x[3]*x[6] - m2*x[1]*x[4]^2 + (-2*m2 + s)*x[1]*x[4]*x[5] + (M[4] - 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] + (M[3] - 2*m2)*x[2]*x[3]*x[4] + (M[3] - 2*m2)*x[2]*x[3]*x[5] + (M[3] - 2*m2)*x[2]*x[3]*x[6] - m2*x[2]*x[4]^2 + (M[3] - 3*m2)*x[2]*x[4]*x[5] + (M[4] - 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] + (M[4] - 2*m2)*x[3]*x[4]*x[6] - m2*x[3]*x[5]^2 + (M[4] - 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 + (M[4] - 2*m2)*x[4]*x[5]*x[6] - m2*x[5]^2*x[6] - m2*x[5]*x[6]^2 parameters = [M[3], M[4], 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] = M[3] χ[1] = 82 weights[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, 0, 0, -1], [0, -1, -1, 0, 1, 1], [0, -1, -1, 1, 0, 1], [0, -1, -1, 1, 1, 0], [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, -1, -1, 0, 0, -1], [0, 0, 0, 1, 1, 1], [-1, -1, -1, 0, -1, 0], [0, 1, 0, 1, 1, 0], [0, -1, -1, 0, 0, 1], [0, -1, -1, 0, 1, 0], [0, -1, -1, 1, 0, 0], [1, -1, -1, 0, 0, 0], [1, 0, 1, 0, 0, 1], [0, -1, -1, -1, -1, 0], [0, -1, -1, -1, 0, -1], [1, 0, 0, 0, 1, 1], [1, 0, 0, 1, 0, 1], [1, 0, 0, 1, 1, 0], [0, 0, 0, 1, 1, 0], [0, 0, 0, 1, 0, 1], [0, -1, -1, 0, 0, 0], [1, 0, 0, 0, 0, 1], [1, 0, 0, 0, 1, 0]] computed_with[1] = ["PLD_sym", "PLD_num"] ################################ # Component 2 ################################ D[2] = M[3] - 4*m2 χ[2] = 102 weights[2] = [[0, -1, -1, 1, 1, 1], [1, -1, -1, 0, 1, 1], [1, -1, -1, 1, 0, 1], [1, -1, -1, 1, 1, 0], [0, -1, -1, 0, 1, 1], [0, -1, -1, 1, 0, 1], [0, -1, -1, 1, 1, 0], [1, -1, -1, 0, 0, 1], [1, -1, -1, 0, 1, 0], [1, -1, -1, 1, 0, 0], [0, -1, -1, 0, 0, 1], [0, -1, -1, 0, 1, 0], [0, -1, -1, 1, 0, 0], [1, -1, -1, 0, 0, 0], [0, -1, -1, 0, 0, 0]] computed_with[2] = ["PLD_sym", "PLD_num"] ################################ # Component 3 ################################ D[3] = M[3] - 9*m2 χ[3] = 114 weights[3] = [[-1, 0, -1, 0, 0, -1], [0, -1, 0, -1, -1, 0]] computed_with[3] = ["PLD_sym"] ################################ # Component 4 ################################ D[4] = M[3] - m2 χ[4] = 110 weights[4] = [[-1, 0, -1, 0, 0, -1], [-1, 0, -1, 0, 0, -1], [-1, 0, -1, 0, 0, -1], [0, -1, 0, -1, -1, 0], [0, -1, 0, -1, -1, 0], [0, -1, 0, -1, -1, 0]] computed_with[4] = ["PLD_sym"] ################################ # Component 5 ################################ D[5] = M[3]^2 - 2*M[3]*M[4] - 2*M[3]*s + M[4]^2 - 2*M[4]*s + s^2 χ[5] = 70 weights[5] = [[-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[5] = ["PLD_num"] ################################ # Component 6 ################################ D[6] = M[3]^2*m2 - 2*M[3]*M[4]*m2 + M[3]*M[4]*s - 5*M[3]*m2*s + M[4]^2*m2 - 5*M[4]*m2*s + 9*m2^2*s + 4*m2*s^2 χ[6] = 114 weights[6] = [[-1, -1, 0, -1, -1, 0], [-1, 0, -1, 0, -1, -1]] computed_with[6] = ["PLD_num"] ################################ # Component 7 ################################ D[7] = M[3]^2*m2 - 5*M[3]*M[4]*m2 + M[3]*M[4]*s - 2*M[3]*m2*s + 4*M[4]^2*m2 + 9*M[4]*m2^2 - 5*M[4]*m2*s + m2*s^2 χ[7] = 112 weights[7] = [[-1, 0, -1, -1, 0, -1], [0, -1, 0, -1, -1, -1]] computed_with[7] = ["PLD_num"] ################################ # Component 8 ################################ D[8] = M[3]^2*m2 - 5//4*M[3]*M[4]*m2 + 1//4*M[3]*M[4]*s + 9//4*M[3]*m2^2 - 5//4*M[3]*m2*s + 1//4*M[4]^2*m2 - 1//2*M[4]*m2*s + 1//4*m2*s^2 χ[8] = 114 weights[8] = [[-1, -1, -1, -1, 0, 0], [0, -1, -1, 0, -1, -1]] computed_with[8] = ["PLD_num"] ################################ # Component 9 ################################ D[9] = M[3]^4*M[4]^2 - 2*M[3]^4*M[4]*m2 + M[3]^4*m2^2 - M[3]^3*M[4]^3 - M[3]^3*M[4]^2*m2 - M[3]^3*M[4]^2*s + 2*M[3]^3*M[4]*m2^2 + 5*M[3]^3*M[4]*m2*s - 4*M[3]^3*m2^2*s + M[3]^2*M[4]^4 - M[3]^2*M[4]^3*m2 - M[3]^2*M[4]^3*s + 3*M[3]^2*M[4]^2*m2^2 + 3*M[3]^2*M[4]^2*m2*s + M[3]^2*M[4]^2*s^2 - 8*M[3]^2*M[4]*m2^2*s - 4*M[3]^2*M[4]*m2*s^2 + 6*M[3]^2*m2^2*s^2 - 2*M[3]*M[4]^4*m2 + 2*M[3]*M[4]^3*m2^2 + 5*M[3]*M[4]^3*m2*s - 8*M[3]*M[4]^2*m2^2*s - 4*M[3]*M[4]^2*m2*s^2 + 10*M[3]*M[4]*m2^2*s^2 + M[3]*M[4]*m2*s^3 - 4*M[3]*m2^2*s^3 + M[4]^4*m2^2 - 4*M[4]^3*m2^2*s + 6*M[4]^2*m2^2*s^2 - 4*M[4]*m2^2*s^3 + m2^2*s^4 χ[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] = M[3]^4*m2^2 - 2*M[3]^4*m2*s + M[3]^4*s^2 - 4*M[3]^3*M[4]*m2^2 + 5*M[3]^3*M[4]*m2*s - M[3]^3*M[4]*s^2 + 2*M[3]^3*m2^2*s - M[3]^3*m2*s^2 - M[3]^3*s^3 + 6*M[3]^2*M[4]^2*m2^2 - 4*M[3]^2*M[4]^2*m2*s + M[3]^2*M[4]^2*s^2 - 8*M[3]^2*M[4]*m2^2*s + 3*M[3]^2*M[4]*m2*s^2 - M[3]^2*M[4]*s^3 + 3*M[3]^2*m2^2*s^2 - M[3]^2*m2*s^3 + M[3]^2*s^4 - 4*M[3]*M[4]^3*m2^2 + M[3]*M[4]^3*m2*s + 10*M[3]*M[4]^2*m2^2*s - 4*M[3]*M[4]^2*m2*s^2 - 8*M[3]*M[4]*m2^2*s^2 + 5*M[3]*M[4]*m2*s^3 + 2*M[3]*m2^2*s^3 - 2*M[3]*m2*s^4 + M[4]^4*m2^2 - 4*M[4]^3*m2^2*s + 6*M[4]^2*m2^2*s^2 - 4*M[4]*m2^2*s^3 + m2^2*s^4 χ[10] = 114 weights[10] = [[-1, -1, -1, -1, -1, 0], [-1, -1, -1, 0, -1, -1]] computed_with[10] = ["PLD_num"] ################################ # Component 11 ################################ D[11] = M[3]^4*m2^2 - 4*M[3]^3*M[4]*m2^2 + M[3]^3*M[4]*m2*s - 4*M[3]^3*m2^2*s + 6*M[3]^2*M[4]^2*m2^2 - 4*M[3]^2*M[4]^2*m2*s + M[3]^2*M[4]^2*s^2 + 10*M[3]^2*M[4]*m2^2*s - 4*M[3]^2*M[4]*m2*s^2 + 6*M[3]^2*m2^2*s^2 - 4*M[3]*M[4]^3*m2^2 + 5*M[3]*M[4]^3*m2*s - M[3]*M[4]^3*s^2 - 8*M[3]*M[4]^2*m2^2*s + 3*M[3]*M[4]^2*m2*s^2 - M[3]*M[4]^2*s^3 - 8*M[3]*M[4]*m2^2*s^2 + 5*M[3]*M[4]*m2*s^3 - 4*M[3]*m2^2*s^3 + M[4]^4*m2^2 - 2*M[4]^4*m2*s + M[4]^4*s^2 + 2*M[4]^3*m2^2*s - M[4]^3*m2*s^2 - M[4]^3*s^3 + 3*M[4]^2*m2^2*s^2 - M[4]^2*m2*s^3 + M[4]^2*s^4 + 2*M[4]*m2^2*s^3 - 2*M[4]*m2*s^4 + m2^2*s^4 χ[11] = 114 weights[11] = [[-1, -1, 0, -1, -1, -1], [-1, 0, -1, -1, -1, -1]] computed_with[11] = ["PLD_num"] ################################ # Component 12 ################################ D[12] = M[4] χ[12] = 82 weights[12] = [[0, 1, 1, -1, 1, -1], [1, 0, 1, -1, 1, -1], [1, 1, 0, -1, 1, -1], [1, 1, 1, -1, 0, -1], [-1, -1, 0, -1, 0, 0], [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], [0, 0, -1, 0, -1, -1], [1, 1, 0, -1, 0, -1], [-1, -1, 0, -1, 0, -1], [0, 0, 1, 0, 1, 1], [0, 1, 1, 0, 1, 0], [-1, 0, 0, -1, -1, -1], [0, 0, 0, -1, 1, -1], [0, 0, 1, -1, 0, -1], [0, 1, 0, -1, 0, -1], [1, 0, 1, 0, 1, 0], [0, -1, -1, -1, 0, -1], [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, 1, 0, 1, 0], [0, 1, 1, 0, 0, 0], [0, 0, 0, -1, 0, -1], [1, 0, 0, 0, 1, 0], [1, 1, 0, 0, 0, 0]] computed_with[12] = ["PLD_sym", "PLD_num"] ################################ # Component 13 ################################ D[13] = M[4] - 4*m2 χ[13] = 102 weights[13] = [[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, -1]] computed_with[13] = ["PLD_sym", "PLD_num"] ################################ # Component 14 ################################ D[14] = M[4] - 9*m2 χ[14] = 114 weights[14] = [[-1, -1, 0, -1, 0, 0], [0, 0, -1, 0, -1, -1]] computed_with[14] = ["PLD_sym"] ################################ # Component 15 ################################ D[15] = M[4] - m2 χ[15] = 110 weights[15] = [[-1, -1, 0, -1, 0, 0], [-1, -1, 0, -1, 0, 0], [-1, -1, 0, -1, 0, 0], [0, 0, -1, 0, -1, -1], [0, 0, -1, 0, -1, -1], [0, 0, -1, 0, -1, -1]] computed_with[15] = ["PLD_sym"] ################################ # Component 16 ################################ D[16] = m2 χ[16] = 28 weights[16] = [[-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[16] = ["PLD_sym", "PLD_num"] ################################ # Component 17 ################################ D[17] = m2 - 1//4*s χ[17] = 102 weights[17] = [[-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[17] = ["PLD_sym", "PLD_num"] ################################ # Component 18 ################################ D[18] = m2 - 1//9*s χ[18] = 114 weights[18] = [[-1, 0, -1, -1, 0, 0], [0, -1, 0, 0, -1, -1]] computed_with[18] = ["PLD_sym"] ################################ # Component 19 ################################ D[19] = m2 - s χ[19] = 110 weights[19] = [[-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[19] = ["PLD_sym"] ################################ # Component 20 ################################ D[20] = s χ[20] = 82 weights[20] = [[-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, 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]] computed_with[20] = ["PLD_sym", "PLD_num"]