def circuit(weight):
init_state = np.flip(np.array([1,1,0,0,0]))
qml.BasisState(init_state, wires=wires)
DoubleExcitationUnitary(weight, wires=pphh)
def circuit(weight):
init_state = np.flip(np.array([1, 1, 0, 0, 0]))
qml.BasisState(init_state, wires=range(N))
DoubleExcitationUnitary(weight, wires1=wires1, wires2=wires2)
return [qml.expval(qml.PauliZ(w)) for w in range(N)]