Return whether [op_a, [op_b, op_a_prime]] is trivially zero.

Assumes all the operators are FermionOperators from the dual basis Hamiltonian, broken into the form i^j^ i j + c_i(i^ i) + c_j(j^ j) or i^ j + j^ i, where i and j are modes and c is a constant. For the full dual basis Hamiltonian, i^ i and j^ j can have distinct coefficients c_i and c_j: for jellium they are necessarily the same. If this is the case, jellium_only should be set to True.

The operators are determined by the indices they act on and by whether they are hopping operators (i^ j + j^ i) or number operators (i^ j^ i j + c_i(i^ i) + c_j(j^ j)). a, b, and a_prime are shorthands for alpha, beta, and alpha_prime.

indices_alpha (set): The indices term_alpha acts on. indices_beta (set): The indices term_beta acts on. indices_alpha_prime (set): The indices term_alpha_prime acts on. is_hopping_operator_alpha (bool): Whether term_alpha is a hopping operator. is_hopping_operator_beta (bool): Whether term_beta is a hopping operator. is_hopping_operator_alpha_prime (bool): Whether term_alpha_prime is a hopping operator. jellium_only (bool): Whether the terms are only from the jellium Hamiltonian, i.e. if c_i = c for all number operators i^ i or if it depends on i.

Whether or not the double commutator is trivially zero.