openfermion.linalg.erpa_eom_hamiltonian

Evaluate \(\sum_{a,b,c,d}h_{a, b, d, c}<\psi[p^ q, [a^ b^ c d, r^ s]]\psi>\)

openfermion.linalg.erpa_eom_hamiltonian(
    h_ijkl: numpy.ndarray, tpdm: numpy.ndarray, p: int, q: int, r: int, s: int
) -> Union[float, complex]

h_ijkl two-body integral tensors of full reduced Hamiltonian
tpdm two-RDM
p left creation op index
q left annihilation op index
r right creation op index
s right annihilation op index

float or complex number