Qedma function fails on simple 12Q circuit

[nonhermitian]

Steps to reproduce the problem

N = 12
qc = QuantumCircuit(N)

qc.x(range(N))
qc.h(range(N))

for kk in range(N // 2, 0, -1):
    qc.ch(kk, kk - 1)
for kk in range(N // 2, N - 1):
    qc.ch(kk, kk + 1)

op = SparsePauliOp('Z'*12)

qesem_function  = catalog.load('qedma/qesem')

qesem_job = qesem_function.run(
    backend_name=backend.name,
    action = "execute",
    instance = backend._instance,
    pubs = pubs,
)

What is the current behavior?

It fails with

{'error': 'An unexpected error occurred. Please try again later. If the issue persists, please contact support.'}

What is the expected behavior?

It should work because it is a standard input that works on other functions like the IBM circuit function

[johannesgreiner]

for reference, the underlying cause of the failure seems to be that the construction of the PUB should allow single observables outside of an array. This example works:

qc = QuantumCircuit(12)

qc.x(range(12))
qc.h(range(12))

for kk in range(12 // 2, 0, -1):
    qc.ch(kk, kk - 1)
for kk in range(12 // 2, 12 - 1):
    qc.ch(kk, kk + 1)

op = SparsePauliOp('Z'*12)

pubs = [(qc, [op])]

options = {
    "default_precision": 0.1,
    "execution_mode": "batch",
}
   
# Run the function
job = qesem_function.run(
   pubs=pubs,
   instance=IQP_INSTANCE,
   backend_name = BACKEND,
   options=options,
)