A computer-implemented method for determining a control sequence for performing a sequence of
qubit interactions on a plurality of qubits on a
quantum device to simulate a
quantum many-body Hamiltonian. The method includes determining a sequence of two-
qubit interactions based on decomposing a multi-
qubit interaction term (120) into a sequence of three interaction terms, where a first interaction term (111) includes a primal operator O, a second interaction term (112) includes an ancillary operator H, and a third interaction term (113) includes a negative of the primal operator, where H and O are each a
tensor product of at least two Pauli matrices, and where the first interaction term, the second interaction term, and the third interaction term produced by the
decomposition include fewer qubits than the original multi-qubit interaction term. The method further includes iterative
decomposition of the interaction terms in terms of the primal and ancillary operators until the multi-qubit interaction term (110) is decomposed into a sequence of two-qubit interaction terms (120). The present invention further relates to computer programs,
quantum circuits,
quantum devices, and methods for solving problems involving many-body interactions.