Paper data;

These circuits are to reproduce the results of Fig. 4f (\Delta \mathcal{C}^{(4)}) of the paper "Practical Quantum Advantage in Computing Local Observables".

The files are meant to be used in the following way. Imagine one had some `cirq` backend which could run on 65 qubits. Then, the results are obtained as follows:

import numpy as np

import qsimcirq

import cirq

import pickle

simulator = qsimcirq.QSimSimulator()

C4_circuits = pickle.load(open('circuits_cycle_23_1020_iswaps', 'rb'))

C4_mc_circuits = pickle.load(open('circuits_cycle_23_1020_iswaps_mc_list', 'rb'))

msmt_qubit = cirq.GridQubit(10, 9)

C4 = np.array([simulator.simulate_expectation_values(C[:-1], observables=[cirq.Z.on(msmt_qubit)])[0].real for C in C4_circuits])

C4_mc = np.array([simulator.simulate_expectation_values(C[:-1], observables=[cirq.Z.on(msmt_qubit)])[0].real for MClist in C4_mc_circuits for C in MClist]).reshape((20, 50)).mean(axis=0)

answer = C4 - C4_mc