qiskit.ignis.mitigation.expval_meas_mitigator_circuits¶

expval_meas_mitigator_circuits
(num_qubits, method='CTMP', labels=None)[source]¶ Generate measurement error mitigator circuits and metadata.
Use the
ExpvalMeasMitigatorFitter
class to fit the execution results to construct a calibrated expectation value measurement error mitigator. Parameters
num_qubits (
int
) – the number of qubits to calibrate.method (
Optional
[str
]) – the mitigation method'complete'
,'tensored'
, or'CTMP'
.labels (
Optional
[List
[str
]]) – Optional, custom labels to run for calibration. If None the method will determine the default label values.
 Returns
 (circuits, metadata) the measurement error characterization
circuits, and metadata for the fitter.
 Return type
tuple
 Mitigation Method:
The
'complete'
method will generate all \(2^n\) computational basis states measurement circuits and fitting will return aCompleteExpvalMeasMitigator
. This method should only be used for small numbers of qubits.The
'tensored'
method will generate two input state circuits of the all 0 and all 1 states on number of qubits unless custom labels are specified. Ftting will return aTensoredExpvalMeasMitigator
. This method assumes measurement errors are uncorrelated between qubits.The
'CTMP'
method will generate input state circuits, unless custom labels are specified. The default input states must obey the following cirterion: for every pair of qubits, projection of the input states on the two qubits contains all four possible assignments to the qubits (00, 01, 10, 11). For n<7, these would be the all 1 state and the \(n\) states with a single qubit in the 1 state and all others in the 0 state (also the all 0 state, if n<3). For n>=7, these would be the all 0 state, the all 1 state, and 2*ceil(log2(n)) states resulting from the following procedure: For each qubits, write its index in binary form, horizontically. For example: with 8 qubits 0, 1, 2,…, 7, for qubit 7 we write: 1 1 1 And for all 8 qubits, we obtain ceil(log2(n))=3 lines: 00001111 00110011 01010101 One can see that the every column is the binary form of the column number. Then write again, the same lines, negated: 11110000 11001100 10101010 The all 0 and all 1 states guarantee that each pair of qubits has input states with projections 00 and 11. The other lines guarantee the projections 01 and 10 (since the qubits are different, when written in binary form, there must be a digit in which they differ). Fitting will return aCTMPExpvalMeasMitigator
.
Example
The following example shows calibrating a 5qubit expectation value measurement error mitigator using the
'tensored'
method.from qiskit import execute from qiskit.test.mock import FakeVigo import qiskit.ignis.mitigation as mit backend = FakeVigo() num_qubits = backend.configuration().num_qubits # Generate calibration circuits circuits, metadata = mit.expval_meas_mitigator_circuits( num_qubits, method='tensored') result = execute(circuits, backend, shots=8192).result() # Fit mitigator mitigator = mit.ExpvalMeasMitigatorFitter(result, metadata).fit() # Plot fitted Nqubit assignment matrix mitigator.plot_assignment_matrix()
<matplotlib.axes._subplots.AxesSubplot at 0x7f6ddff405b0>
The following shows how to use the above mitigator to apply measurement error mitigation to expectation value computations
from qiskit import QuantumCircuit # Test Circuit with expectation value 1. qc = QuantumCircuit(num_qubits) qc.x(range(num_qubits)) qc.measure_all() # Execute shots = 8192 seed_simulator = 1999 result = execute(qc, backend, shots=8192, seed_simulator=1999).result() counts = result.get_counts(0) # Expectation value of Z^N without mitigation expval_nomit, error_nomit = mit.expectation_value(counts) print('Expval (no mitigation): {:.2f} ± {:.2f}'.format( expval_nomit, error_nomit)) # Expectation value of Z^N with mitigation expval_mit, error_mit = mit.expectation_value(counts, meas_mitigator=mitigator) print('Expval (with mitigation): {:.2f} ± {:.2f}'.format( expval_mit, error_mit))
Expval (no mitigation): 0.66 ± 0.01 Expval (with mitigation): 1.02 ± 0.01