Note

This is the documentation for the current state of the development branch of Qiskit Experiments. The documentation or APIs here can change prior to being released.

Quantum Volume

Quantum Volume (QV) is a single-number metric that can be measured using a concrete protocol on near-term quantum computers of modest size. The QV method quantifies the largest random circuit of equal width and depth that the computer successfully implements. Quantum computing systems with high-fidelity operations, high connectivity, large calibrated gate sets, and circuit rewriting toolchains are expected to have higher quantum volumes. See the Qiskit Textbook for an explanation on the QV method, which is described in Refs. [1] [2].

The Quantum Volume is determined by the largest successful circuit depth $d_{max}$, and equals to $2^{d_{max}}$. In this experiment, we generate QV circuits using the qiskit.circuit.library.QuantumVolume class on $d$ qubits, which contain $d$ layers, where each layer consists of random 2-qubit unitary gates from $SU(4)$, followed by a random permutation on the $d$ qubit. Then these circuits run on the quantum backend and on an ideal simulator (either qiskit_aer.AerSimulator or qiskit.quantum_info.Statevector).

A depth $d$ QV circuit is successful if it has mean heavy-output probability > 2/3 with confidence level > 0.977 (corresponding to z_value = 2), and at least 100 trials have been ran.

Note

This tutorial requires the qiskit-aer and qiskit-ibm-runtime packages to run simulations. You can install them with python -m pip install qiskit-aer qiskit-ibm-runtime.

from qiskit_experiments.framework import BatchExperiment
from qiskit_experiments.library import QuantumVolume

# For simulation
from qiskit_aer import AerSimulator
from qiskit_ibm_runtime.fake_provider import FakeSydneyV2

backend = AerSimulator.from_backend(FakeSydneyV2())

QV experiment

To run the QV experiment we need need to provide the following QV parameters, in order to generate the QV circuits and run them on a backend and on an ideal simulator:

• qubits: The number of qubits or list of physical qubits for the experiment.

• trials: The number of trials to run the quantum volume circuit (the default is 100).

• seed: Seed or generator object for random number generation. If None then default_rng will be used.

• simulation_backend: The simulator backend to use to generate the expected results. the simulator must have a save_probabilities method. If None, AerSimulator will be used (in case AerSimulator is not installed, Statevector will be used).

Note: In some cases, 100 trials are not enough to obtain a QV greater than 1 for the specified number of qubits. In this case, adding more trials may reduce the error bars to allow passing the threshold.

The analysis results of the QV Experiment are:

• The mean heavy-output probabilities (HOP) and standard deviation

• The calculated quantum volume, which will be None if the experiment does not pass the threshold

Extra data included in the analysis results includes

• The heavy HOPs for each trial

• Confidence level (should be greater than 0.977)

• The number of trials and depth of the QV circuits

• Whether the QV circuit was successful

qubits = tuple(range(4)) # Can use specific qubits. for example [2, 4, 7, 10]

qv_exp = QuantumVolume(qubits, seed=42)
# Transpile options like optimization_level affect only the real device run and not the simulation run
# Run options affect both simulation and real device runs
qv_exp.set_transpile_options(optimization_level=3)

# Run experiment
expdata = qv_exp.run(backend).block_for_results()

# View result data
display(expdata.figure(0))

display(expdata.analysis_results(dataframe=True))

	name	experiment	components	value	quality	backend	run_time	HOPs	two_sigma	depth	trials	success	confidence
a94885b6	mean_HOP	QuantumVolume	[Q0, Q1, Q2, Q3]	0.73+/-0.04	bad	aer_simulator_from(fake_sydney)	None	[0.662109375, 0.6796875, 0.7568359375, 0.69824...	0.088805	4	100	None	None
1c1afd27	quantum_volume	QuantumVolume	[Q0, Q1, Q2, Q3]	1	bad	aer_simulator_from(fake_sydney)	None	None	None	4	100	False	0.922702

Adding trials

Adding more trials may reduce the error bars to allow passing the threshold (beside the error bars - QV experiment must have at least 100 trials to be considered successful). In case you want to add less than 100 additional trials, you can modify the amount of trials added before re-running the experiment.

qv_exp.set_experiment_options(trials=60)
expdata2 = qv_exp.run(backend, analysis=None).block_for_results()
expdata2.add_data(expdata.data())
qv_exp.analysis.run(expdata2).block_for_results()

# View result data
display(expdata2.figure(0))
display(expdata2.analysis_results(dataframe=True))

	name	experiment	components	value	quality	backend	run_time	HOPs	two_sigma	depth	trials	success	confidence
5c4fdb6a	mean_HOP	QuantumVolume	[Q0, Q1, Q2, Q3]	0.730+/-0.035	bad	aer_simulator_from(fake_sydney)	None	[0.6337890625, 0.7080078125, 0.7529296875, 0.7...	0.070173	4	160	None	None
94186e8e	quantum_volume	QuantumVolume	[Q0, Q1, Q2, Q3]	1	bad	aer_simulator_from(fake_sydney)	None	None	None	4	160	False	0.9651

Calculating Quantum Volume using a batch experiment

Run the QV experiment with an increasing number of qubits to check what is the maximum Quantum Volume for the specific device. To reach the real system’s Quantum Volume, one must run more trials and additional enhancements might be required (See Ref. [2] for details).

exps = [QuantumVolume(tuple(range(i)), trials=200) for i in range(3, 6)]

batch_exp = BatchExperiment(exps)
batch_exp.set_transpile_options(optimization_level=3)

# Run
batch_expdata = batch_exp.run(backend).block_for_results()

Extracting the maximum Quantum Volume.

qv_values = batch_expdata.analysis_results("quantum_volume", dataframe=True).value.tolist()

print(f"Max quantum volume is: {max(qv_values)}")

Max quantum volume is: 32

for i in range(batch_exp.num_experiments):
    print(f"\nComponent experiment {i}")
    display(batch_expdata.figure(i))
display(batch_expdata.analysis_results(dataframe=True))

Component experiment 0

Component experiment 1

Component experiment 2

	name	experiment	components	value	quality	backend	run_time	HOPs	two_sigma	depth	trials	success	confidence
3b36a9e4	mean_HOP	QuantumVolume	[Q0, Q1, Q2]	0.795+/-0.029	good	aer_simulator_from(fake_sydney)	None	[0.8466796875, 0.765625, 0.7001953125, 0.83593...	0.05714	3	200	None	None
df2f407c	quantum_volume	QuantumVolume	[Q0, Q1, Q2]	8	good	aer_simulator_from(fake_sydney)	None	None	None	3	200	True	0.999996
00045234	mean_HOP	QuantumVolume	[Q0, Q1, Q2, Q3]	0.729+/-0.031	bad	aer_simulator_from(fake_sydney)	None	[0.7119140625, 0.6689453125, 0.646484375, 0.73...	0.062833	4	200	None	None
afd149e7	quantum_volume	QuantumVolume	[Q0, Q1, Q2, Q3]	1	bad	aer_simulator_from(fake_sydney)	None	None	None	4	200	False	0.976984
0833e575	mean_HOP	QuantumVolume	[Q0, Q1, Q2, Q3, Q4]	0.732+/-0.031	good	aer_simulator_from(fake_sydney)	None	[0.708984375, 0.703125, 0.7587890625, 0.712890...	0.062614	5	200	None	None
96d97eaa	quantum_volume	QuantumVolume	[Q0, Q1, Q2, Q3, Q4]	32	good	aer_simulator_from(fake_sydney)	None	None	None	5	200	True	0.98202

References

[1]

Andrew W. Cross, Lev S. Bishop, Sarah Sheldon, Paul D. Nation, and Jay M. Gambetta, Validating quantum computers using randomized model circuits, Phys. Rev. A 100, 032328 (2019). https://arxiv.org/pdf/1811.12926

[2] (1,2)

Petar Jurcevic et. al. Demonstration of quantum volume 64 on a superconducting quantum computing system, https://arxiv.org/pdf/2008.08571

See also

• API documentation: quantum_volume

• Qiskit Textbook: Measuring Quantum Volume