Protocols¶
Generates a 

Generates a 

Generates a 

Generates a 

Generates a 

Generates a 

Generates a 

Generates a 

Generates a 

Computes the quantum capacity (QCAP) bound of the given circuit. 

Randomly compiles the given circuit into many new random circuits which implement the same algorithm. 
Make CB¶

trueq.
make_cb
(cycle, n_random_cycles, n_circuits=30, n_decays=20, targeted_errors=None, twirling_group=None, propagate_correction=False)¶ Generates a
CircuitCollection
to estimate the process fidelity of a cycle using cycle benchmarking (CB). See also the CB guide.import trueq as tq # generate a circuit collection to run CB on the 0th qubit, with 30 circuits, # for each length in [4, 20] and each with 3 randomly chosen Pauli decay strings circuits = tq.make_cb({0: tq.Gate.x}, [4, 20], 30, 3) # draw the first circuit circuits[0].draw()
 Parameters
cycle (
Cycle
dict
) – The cycle to benchmark or an argument to instantiate the cycle.n_random_cycles (
Iterable
) – A list of positive integers specifying how many random cycles will be generated during the protocol, e.g.[4, 20]
.n_circuits (
int
) – The number of circuits for each number of random cycles.n_decays (
int
) – An integer specifying the total number of randomly chosen Pauli decay strings used to measure the process infidelity or the probability of each error. Warning: Setting this value lower thanmin(20, 4 ** n_qubits  1)
may result in a biased estimate of the process fidelity, and setting this value lower thanmin(40, 4 ** n_qubits  1)
may result in a biased estimate of the probability for nonidentity errors.targeted_errors (
Iterable
) – A list of Pauli strings, e.g.["ZIZIZ", "XYXYX"]
that specify which errors to target. By default, the PauliI * n_qubits
is used to estimate the probability of no error, i.e. the process fidelity.twirling_group (
Iterable
) – The twirling group to use in this protocol. This can be one of"P"
(Pauli),"C"
(Clifford),"SU"
(special unitary), or"I"
(Identity). The twirling group will be the tensor product of the specified group on the keys of the cycle to be benchmarked. A mixture of groups on different labels can be specified by an iterable of groups and labels, for example,[["P", 0], ["C", 1, 2], ["C", 3]]
, where["C", 1, 2]
specifies the twoqubit Clifford group acting on the labels[1, 2]
and["C", 3]
specifies the onequbit Clifford group acting on the label3
. Any unspecified elements oflabels
will be set to"P"
. Warning:"I"
is provided as a convenience method for expert users to identify coherent errors and/or crosstalk even though an exponential decay is not expected.propagate_correction (
bool
) – Whether to propagate correction gates to the end of the circuit or compile them in to neighbouring cycles. Warning: can result in arbitrary twoqubit gates at the end of the circuit!
 Returns
A collection of CB circuits.
 Return type
Crosstalk Diagnostics¶

trueq.
make_crosstalk_diagnostics
(labels, n_random_cycles, n_circuits=30, subsets=None, include_xrb=True, twirling_group=None)¶ Generates a
CircuitCollection
to measure the discrepancy between gate quality while applying gates simultaneously versus while applying gates in an isolated way to each individual system. See also the CTD guide.This discrepancy is assessed by running SRB simultaneously on the specified systems, as well as running it on each system in turn. Optionally, and
True
by default, the coherence of crosstalk errors is assessed using XRB. Therefore, this function is equivalent to concatenating the circuits from multiple calls tomake_srb()
(and optionallymake_xrb()
) with different label configurations.import trueq as tq # generate a circuit collection to run crosstalk diagnostics in single qubit # mode for qubits 5, 6, 7, and 8 circuits = tq.make_crosstalk_diagnostics([5, 6, 7, 8], [4, 100]) # reduce the number of circuits to perform by excluding XRB circuits # this means we cannot distinguish between coherent errors (due to static # crosstalk) and incoherent errors (due to fluctuating crosstalk) circuits = tq.make_crosstalk_diagnostics( [5, 6, 7, 8], [4, 100], include_xrb=False ) # we can also twirl some pairs of systems using entangling gates circuits = tq.make_crosstalk_diagnostics( [5, [6, 7], 8], [4, 100], include_xrb=False ) # using the subsets option, we can customize exactly which subsets of the full # simultaneous twirl are performed in isolation. Here, the default value would # have resulted in subsets [[5],[6],[[7, 8]]] circuits = tq.make_crosstalk_diagnostics( [5, 6, [7, 8]], [4, 100], include_xrb=False, subsets=[[5], [5, 6], [[7, 8]]] )
 Parameters
labels (
Iterable
) – A list of which sets of system labels are to be twirled together in each circuit, e.g.[3, [1, 2], 4]
.n_random_cycles (
Iterable
) – A list of positive integers specifying how many random cycles will be generated during the protocol, e.g.[6, 20]
.n_circuits (
int
) – The number of circuits for at eachn_random_cycles
for each protocol.subsets (
Iterable
) – A list of subsets of the given labels to be used. The default value ofNone
results in each subset being a member of the given labels. See the last example above.include_xrb (
bool
) – Whether to include XRB circuits in the output.twirling_group (
Iterable
) – The twirling group to use in this protocol. This can be either"C"
(Clifford) or"SU"
(special unitary). The twirling group will be the tensor product of the specified group on the labels to be benchmarked. A mixture of groups on different labels can be specified by an iterable of groups and labels such as[["SU", 0], ["C", 1, 2], ["C", 3]]
, where["C", 1, 2]
specifies the twoqubit Clifford group acting on the labels[1, 2]
and["C", 3]
specifies the onequbit Clifford group acting on the label3
. Any unspecified labels will be set to"C"
.
 Returns
A collection of circuits to diagnose crosstalk.
 Return type
Make IRB¶

trueq.
make_irb
(cycle, n_random_cycles, n_circuits=30, twirling_group=None, propagate_correction=False)¶ Generates a
CircuitCollection
to estimate the process fidelity of specific gates using (simultaneous) interleaved randomized benchmarking (IRB). See also the IRB guide.import trueq as tq # generate a circuit collection to run single qubit IRB on an X gate acting on # qubit 0, with 30 random circuits for each circuit length in [5, 40, 60, 100] circuits = tq.make_irb({0: tq.Gate.x}, [5, 40, 60, 100], 30) # next, generate circuits to run IRB on a cycle with an X gate acting on qubit 0 # and a CZ gate on qubits (2, 3), with 20 random circuits at each circuit length circuits = tq.make_irb({0: tq.Gate.x, (2, 3): tq.Gate.cz}, [2, 100], 20) # finding the Pauli matrix compiled into the first circuit print(circuits[0].key.compiled_pauli) # draw the first circuit circuits[0].draw()
XZI
0 2 3 Key: cycle: Cycle((0,): Gate.x, (2, 3): Gate.cz, immutable=True) twirl: (('C', 0), ('C', 2, 3)) protocol: IRB compiled_pauli: XZI n_random_cycles: 2 Labels: (0,) Name: Gate.cliff22 Aliases: Gate.cliff22 Generators: Y: 69.28 X: 69.28 Z: 69.28 0.50 0.50j 0.50 0.50j 0.50 0.50j 0.50 0.50j 22 Labels: (2, 3) Name: Gate(IX, IY, ...) Locally Equivalent: iSWAP Generators: IZ: 124.64 XX: 124.64 XY: 124.64 ZI: 60.00 ZX: 60.00 IX: 60.00 YZ: 55.36 ... 0.50 0.50 0.50j 0.50j 0.50 0.50 0.50j 0.50j 0.50 0.50 0.50j 0.50j 0.50 0.50 0.50j 0.50j Imm Labels: (0,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (2, 3) Name: Gate.cz Aliases: Gate.cz Locally Equivalent: CNOT Generators: IZ: 90.00 ZI: 90.00 ZZ: 90.00 1.00 1.00 1.00 1.00 CZ CZ Labels: (0,) Name: Gate.h Aliases: Gate.h Gate.cliff12 Generators: Z: 127.28 X: 127.28 0.71j 0.71j 0.71j 0.71j H Labels: (2, 3) Name: Gate(IX, IY, ...) Locally Equivalent: CNOT Generators: YI: 90.00 IX: 63.64 YX: 63.64 YY: 63.64 IY: 63.64 IZ: 45.00 YZ: 45.00 0.50j 0.50 0.50 0.50j 0.50j 0.50 0.50 0.50j 0.50 0.50j 0.50j 0.50 0.50 0.50j 0.50j 0.50 Imm Labels: (0,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (2, 3) Name: Gate.cz Aliases: Gate.cz Locally Equivalent: CNOT Generators: IZ: 90.00 ZI: 90.00 ZZ: 90.00 1.00 1.00 1.00 1.00 CZ CZ Labels: (0,) Name: Gate.cliff4 Aliases: Gate.cliff4 Generators: X: 90.00 0.71j 0.71 0.71 0.71j 4 Labels: (2, 3) Name: Gate(IX, IY, ...) Locally Equivalent: CNOT Generators: YZ: 80.50 XI: 80.50 YY: 80.50 ZI: 80.50 YX: 80.50 IZ: 49.55 XY: 49.55 ... 0.50j 0.50 0.50j 0.50 0.50 0.50j 0.50 0.50j 0.50 0.50j 0.50 0.50j 0.50j 0.50 0.50j 0.50 Imm Labels: (0,) Name: Meas() M Labels: (2,) Name: Meas() M Labels: (3,) Name: Meas() M Note
We invert each random sequence up to an independent random Pauli matrix (accessible as shown in the example above) to diagnose errors to that are missed by always returning to the initial state (see arXiv:1901.00535).
Warning
The estimate of the process fidelity obtained by taking the ratio of the process fidelities from IRB and SRB is subject to a large systematic uncertainty that is typically not reported. This systematic uncertainty arises because of the many ways in which the noise in the twirling group can combine with the noise in the interleaved gates. Both the standard interleaved estimate and the systematic uncertainty are automatically computed by
fit
if the circuit collection contains SRB circuits. The systematic uncertainty can be reduced if the circuit collection also contains XRB circuits. Parameters
cycle (
Cycle
dict
) – The cycle to benchmark or an argument to instantiate the cycle.n_random_cycles (
Iterable
) – A list of positive integers specifying how many random cycles will be generated during the protocol, e.g.[6, 20]
.n_circuits (
int
) – The number of circuits for each number of random cycles.twirling_group (
Iterable
) – The twirling group to use in this protocol. This can be either"C"
(Clifford) or"SU"
(special unitary). The twirling group will be the tensor product of the specified group on the labels of the cycle to be benchmarked. A mixture of groups on different keys can be specified by an iterable of groups and labels such as[["P", 0], ["C", 1, 2], ["C", 3]]
, where["C", 1, 2]
specifies the twoqubit Clifford group acting on the labels[1, 2]
and["C", 3]
specifies the onequbit Clifford group acting on the label3
. Any unspecified keys of the cycle to be benchmarked will be set to"C"
.propagate_correction (
bool
) – Whether to propagate correction gates to the end of the circuit or compile them in to neighbouring cycles. Warning: can result in arbitrary twoqubit gates at the end of the circuit!
 Returns
A collection of IRB circuits.
 Return type
Make KNR¶

trueq.
make_knr
(cycle, n_random_cycles, n_circuits=30, n_bodies=1, twirling_group=None, propagate_correction=False)¶ Generates a
CircuitCollection
to estimate the probabilities of all errors acting on all sets of subsystems targeted by a combination of k gates using kbody noise reconstruction (KNR). See also the KNR guide.For example, if the input cycle is
{0: tq.Gate.id, (1, 3): tq.Gate.cnot, 2: tq.Gate.x}
andn_bodies=2
, then data will be present to reconstruct any Pauli error rate on the subsystems[0, 1, 3]
,[0, 2]
,[1, 2, 3]
(and any subsystems thereof, e.g.[1, 2]
but not[0, 1, 2]
).import trueq as tq # generate a circuit collection to reconsturct twobody marginal # error distributions # on a 4qubit device using 30 circuits for each length in [6, 20] # and for each Pauli subspace in a set of log2(4) subspaces circuits = tq.make_knr({j: tq.Gate.x for j in range(4)}, [6, 20], 30) # draw the first circuit circuits[0].draw()
0 1 2 3 Key: cycle: Cycle((0,): Gate.x, (1,): Gate.x, (2,): Gate.x, (3,): Gate.x, immutable=True) twirl: (('P', 0), ('P', 1), ('P', 2), ('P', 3)) n_bodies: 1 protocol: KNR compiled_pauli: ZZXY n_random_cycles: 6 measurement_basis: XXXZ Labels: (0,) Name: Gate.h Aliases: Gate.h Gate.cliff12 Generators: Z: 127.28 X: 127.28 0.71 0.71 0.71 0.71 H Labels: (1,) Name: Gate.cliff7 Aliases: Gate.cliff7 Generators: Y: 90.00 0.71 0.71 0.71 0.71 7 Labels: (2,) Name: Gate.cliff13 Aliases: Gate.cliff13 Generators: Z: 127.28 X: 127.28 0.71 0.71 0.71 0.71 13 Labels: (3,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Imm Labels: (0,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (1,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (2,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (3,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (0,) Name: Gate.y Aliases: Gate.y Gate.cliff2 Generators: Y: 180.00 1.00 1.00 Y Labels: (1,) Name: Gate.id Aliases: Gate.id Gate.i Gate.cliff0 Locally Equivalent: Identity Generators: I: 0.00 1.00 1.00 ID Labels: (2,) Name: Gate.z Aliases: Gate.z Gate.cliff3 Generators: Z: 180.00 1.00 1.00 Z Labels: (3,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Imm Labels: (0,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (1,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (2,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (3,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (0,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (1,) Name: Gate.z Aliases: Gate.z Gate.cliff3 Generators: Z: 180.00 1.00 1.00 Z Labels: (2,) Name: Gate.z Aliases: Gate.z Gate.cliff3 Generators: Z: 180.00 1.00 1.00 Z Labels: (3,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Imm Labels: (0,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (1,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (2,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (3,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (0,) Name: Gate.y Aliases: Gate.y Gate.cliff2 Generators: Y: 180.00 1.00 1.00 Y Labels: (1,) Name: Gate.y Aliases: Gate.y Gate.cliff2 Generators: Y: 180.00 1.00 1.00 Y Labels: (2,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (3,) Name: Gate.id Aliases: Gate.id Gate.i Gate.cliff0 Locally Equivalent: Identity Generators: I: 0.00 1.00 1.00 ID Imm Labels: (0,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (1,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (2,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (3,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (0,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (1,) Name: Gate.z Aliases: Gate.z Gate.cliff3 Generators: Z: 180.00 1.00 1.00 Z Labels: (2,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (3,) Name: Gate.id Aliases: Gate.id Gate.i Gate.cliff0 Locally Equivalent: Identity Generators: I: 0.00 1.00 1.00 ID Imm Labels: (0,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (1,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (2,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (3,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (0,) Name: Gate.id Aliases: Gate.id Gate.i Gate.cliff0 Locally Equivalent: Identity Generators: I: 0.00 1.00 1.00 ID Labels: (1,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (2,) Name: Gate.z Aliases: Gate.z Gate.cliff3 Generators: Z: 180.00 1.00 1.00 Z Labels: (3,) Name: Gate.y Aliases: Gate.y Gate.cliff2 Generators: Y: 180.00 1.00 1.00 Y Imm Labels: (0,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (1,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (2,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (3,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (0,) Name: Gate.cliff7 Aliases: Gate.cliff7 Generators: Y: 90.00 0.71 0.71 0.71 0.71 7 Labels: (1,) Name: Gate.cliff6 Aliases: Gate.cliff6 Generators: Y: 270.00 0.71 0.71 0.71 0.71 6 Labels: (2,) Name: Gate.h Aliases: Gate.h Gate.cliff12 Generators: Z: 127.28 X: 127.28 0.71 0.71 0.71 0.71 H Labels: (3,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Imm Labels: (0,) Name: Meas() M Labels: (1,) Name: Meas() M Labels: (2,) Name: Meas() M Labels: (3,) Name: Meas() M Note
The number of returned circuits scales mildly with the number of total systems but exponentially with the value of
n_bodies
. Parameters
cycle (
Cycle
dict
) – The cycle to benchmark or an argument to instantiate the cycle.n_random_cycles (
Iterable
) – A list of positive integers specifying how many random cycles will be generated during the protocol, e.g.[6, 20]
.n_circuits (
int
) – The number of circuits for each number of random cycles.n_bodies (
int
) – A number of gate bodies to reconstruct marginal probability distributions for.twirling_group (
Iterable
) – The twirling group to use in this protocol. This can be one of"P"
(Pauli),"C"
(Clifford),"SU"
(special unitary), or"I"
(Identity). The twirling group will be the tensor product of the specified group on the keys of the cycle to be benchmarked. A mixture of groups on different labels can be specified by an iterable of groups and labels, for example,[["P", 0], ["C", 1, 2], ["C", 3]]
, where["C", 1, 2]
specifies the twoqubit Clifford group acting on the labels[1, 2]
and["C", 3]
specifies the onequbit Clifford group acting on the label3
. Any unspecified elements oflabels
will be set to"P"
. Warning:"I"
is provided as a convenience method for expert users to identify coherent errors and/or crosstalk even though an exponential decay is not expected.propagate_correction (
bool
) – Whether to propagate correction gates to the end of the circuit or compile them in to neighbouring cycles. Warning: can result in arbitrary twoqubit gates at the end of the circuit!
 Returns
A collection of KNR circuits.
 Return type
Make RCAL¶

trueq.
make_rcal
(labels, stagger=False, independent=True, batch=None)¶ Generates a
CircuitCollection
to measure the readout errors on the provided system labels. These circuits contain \(I\) and \(X\) gates. See the RCAL guide for more information. Parameters
labels (
Iterable
) – A list of the system labels whose readout errors are to be estimated.stagger (
bool
) – Whether \(X\) gates should appear in separate immutable cycles. If the gates are much worse when applied in parallel, then staggering them will reduce systematic errors, provided that the gate duration multiplied by the number of qubits is much shorter than \(T_1\).independent (
bool
) – Whether it can be assumed that the readout error is independent over qubits. Under this assumption, only two circuits are required, whereas \(2^n\) circuits are required when this assumption does not hold. Note that if you measure the readout calibration matrix in both independent and nonindependent modes, any observed nonindependence could potentially be due to gate crosstalk errors.batch (
int
) – A unique identifier for these calibration circuits. Other circuits which are to be calibrated based on these circuits’ results in particular should set theirtrueq.Circuit.key
batch
attribute to the same value.
 Returns
A circuit collection to measure the readout errors on the provided system labels.
 Return type
Make SC¶

trueq.
make_sc
(cycle, n_random_cycles, n_circuits=30, pauli_decays=None, twirling_group=None, propagate_correction=False)¶ Generates a
CircuitCollection
to optimize the fidelity of a cycle using stochastic calibration (SC). See also the SC guide.import trueq as tq # generate a circuit collection to run SC on the 0th qubit, with 30 circuits, # for each length in [4, 20] and each with a Pauli decay string "Z" circuits = tq.make_sc({0: tq.Gate.x}, [4, 20], 30, "Z") # draw the first circuit circuits[0].draw()
0 Key: cycle: Cycle((0,): Gate.x, immutable=True) twirl: (('P', 0),) protocol: SC analyze_decays: ('Z',) compiled_pauli: Y n_random_cycles: 4 measurement_basis: Z Labels: (0,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Imm Labels: (0,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (0,) Name: Gate.z Aliases: Gate.z Gate.cliff3 Generators: Z: 180.00 1.00j 1.00j Z Imm Labels: (0,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (0,) Name: Gate.id Aliases: Gate.id Gate.i Gate.cliff0 Locally Equivalent: Identity Generators: I: 0.00 1.00 1.00 ID Imm Labels: (0,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (0,) Name: Gate.y Aliases: Gate.y Gate.cliff2 Generators: Y: 180.00 1.00j 1.00j Y Imm Labels: (0,) Name: Gate.x Aliases: Gate.x Gate.cliff1 Generators: X: 180.00 1.00 1.00 X Labels: (0,) Name: Gate.y Aliases: Gate.y Gate.cliff2 Generators: Y: 180.00 1.00j 1.00j Y Imm Labels: (0,) Name: Meas() M  Parameters
cycle (
Cycle
dict
) – The cycle to benchmark or an argument to instantiate the cycle.n_random_cycles (
Iterable
) – A list of positive integers specifying how many random cycles will be generated during the protocol, e.g.[6, 20]
.n_circuits (
int
) – The number of circuits for each number of random cycles.pauli_decays (
Iterable
) – A list of Pauli decay strings, e.g.["ZIZIZ", "XYXYX"]
that specify which elements of the diagonalized channel should be estimated. These should be chosen to anticommute with the Hamiltonian terms of a known noise source to be optimized.twirling_group (
Iterable
) – The twirling group to use in this protocol. This can be one of"P"
(Pauli),"C"
(Clifford),"SU"
(special unitary), or"I"
(Identity). The twirling group will be the tensor product of the specified group on the keys of the cycle to be benchmarked. A mixture of groups on different labels can be specified by an iterable of groups and labels, for example,[["P", 0], ["C", 1, 2], ["C", 3]]
, where["C", 1, 2]
specifies the twoqubit Clifford group acting on the labels[1, 2]
and["C", 3]
specifies the onequbit Clifford group acting on the label3
. Any unspecified elements oflabels
will be set to"P"
. Warning:"I"
is provided as a convenience method for expert users to identify coherent errors and/or crosstalk even though an exponential decay is not expected.propagate_correction (
bool
) – Whether to propagate correction gates to the end of the circuit or compile them in to neighbouring cycles. Warning: can result in arbitrary twoqubit gates at the end of the circuit!
 Returns
A collection of SC circuits.
 Return type
Make SRB¶

trueq.
make_srb
(labels, n_random_cycles, n_circuits=30, twirling_group=None)¶ Generates a
CircuitCollection
to estimate the process fidelity of random gates from a group using (simultaneous) streamlined randomized benchmarking (SRB). See also the SRB guide.import trueq as tq # generate a circuit collection to run single qubit SRB on qubit 5 circuits = tq.make_srb([5], [4, 200], 30) # generate a circuit collection to run simultaneous onequbit SRB on qubits # 2 and 7 and twoqubit SRB on the pairs [0, 1] and [4, 9] circuits = tq.make_srb([[0, 1], 2, 7, [4, 9]], [4, 200], 30) # finding the Pauli matrix compiled into the first circuit print(circuits[0].key.compiled_pauli) # draw the first circuit circuits[0].draw()
XXXYIZ
0 1 2 4 7 9 Key: twirl: (('C', 2), ('C', 7), ('C', 0, 1), ('C', 4, 9)) protocol: SRB compiled_pauli: XXXYIZ n_random_cycles: 4 Labels: (0, 1) Name: Gate(IX, IZ, ...) Locally Equivalent: CNOT Generators: ZI: 63.64 ZZ: 63.64 YY: 63.64 YI: 63.64 ZX: 63.64 ZY: 63.64 YZ: 63.64 ... 0.50 0.50j 0.50 0.50j 0.50 0.50j 0.50 0.50j 0.50 0.50j 0.50 0.50j 0.50 0.50j 0.50 0.50j Labels: (2,) Name: Gate.cliff22 Aliases: Gate.cliff22 Generators: Y: 69.28 X: 69.28 Z: 69.28 0.50 0.50j 0.50 0.50j 0.50 0.50j 0.50 0.50j 22 Labels: (4, 9) Name: Gate(IX, IY, ...) Locally Equivalent: iSWAP Generators: YZ: 63.64 ZI: 63.64 IX: 63.64 ZZ: 63.64 IY: 63.64 XY: 63.64 YI: 63.64 ... 0.50 0.50j 0.50 0.50j 0.50 0.50j 0.50 0.50j 0.50 0.50j 0.50 0.50j 0.50 0.50j 0.50 0.50j Labels: (7,) Name: Gate.cliff15 Aliases: Gate.cliff15 Generators: Y: 127.28 Z: 127.28 0.71 0.71j 0.71j 0.71 15 Imm Labels: (0, 1) Name: Gate(IX, IY, ...) Locally Equivalent: iSWAP Generators: IX: 60.00 YY: 60.00 YI: 60.00 ZX: 60.00 XZ: 60.00 YX: 60.00 XI: 34.64 ... 0.50 0.50j 0.50 0.50j 0.50j 0.50 0.50j 0.50 0.50j 0.50 0.50j 0.50 0.50 0.50j 0.50 0.50j Labels: (2,) Name: Gate.cliff6 Aliases: Gate.cliff6 Generators: Y: 90.00 0.71 0.71 0.71 0.71 6 Labels: (4, 9) Name: Gate(IX, IY, ...) Locally Equivalent: iSWAP Generators: YZ: 120.00 YI: 120.00 IZ: 120.00 ZY: 60.00 XX: 60.00 YY: 60.00 ZZ: 60.00 ... 0.50 0.50j 0.50j 0.50 0.50j 0.50 0.50 0.50j 0.50j 0.50 0.50 0.50j 0.50 0.50j 0.50j 0.50 Labels: (7,) Name: Gate.z Aliases: Gate.z Gate.cliff3 Generators: Z: 180.00 1.00 1.00 Z Imm Labels: (0, 1) Name: Gate(IY, XI, ...) Locally Equivalent: CNOT Generators: YI: 103.92 ZI: 103.92 XY: 103.92 IY: 90.00 YY: 34.64 ZY: 34.64 XI: 34.64 0.71j 0.71 0.71j 0.71 0.71 0.71j 0.71 0.71j Labels: (2,) Name: Gate.cliff16 Aliases: Gate.cliff16 Generators: Z: 69.28 X: 69.28 Y: 69.28 0.50 0.50j 0.50 0.50j 0.50 0.50j 0.50 0.50j 16 Labels: (4, 9) Name: Gate(IX, IY, ...) Locally Equivalent: SWAP Generators: YZ: 60.00 ZX: 60.00 YX: 60.00 XZ: 60.00 ZY: 60.00 XY: 60.00 YI: 34.64 ... 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j Labels: (7,) Name: Gate.cliff18 Aliases: Gate.cliff18 Generators: Z: 69.28 X: 69.28 Y: 69.28 0.50 0.50j 0.50 0.50j 0.50 0.50j 0.50 0.50j 18 Imm Labels: (0, 1) Name: Gate(IX, IY, ...) Locally Equivalent: CNOT Generators: XI: 80.50 YI: 80.50 ZZ: 80.50 ZX: 80.50 ZY: 80.50 IZ: 49.55 IX: 49.55 ... 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j Labels: (2,) Name: Gate.y Aliases: Gate.y Gate.cliff2 Generators: Y: 180.00 1.00 1.00 Y Labels: (4, 9) Name: Gate(IX, IY, ...) Locally Equivalent: CNOT Generators: ZI: 61.25 YY: 61.25 XX: 61.25 IX: 61.25 IY: 61.25 YZ: 37.85 XZ: 37.85 ... 0.50j 0.50 0.50 0.50j 0.50j 0.50 0.50 0.50j 0.50j 0.50 0.50 0.50j 0.50j 0.50 0.50 0.50j Labels: (7,) Name: Gate.z Aliases: Gate.z Gate.cliff3 Generators: Z: 180.00 1.00 1.00 Z Imm Labels: (0, 1) Name: Gate(IX, IY, ...) Locally Equivalent: CNOT Generators: XI: 80.50 YY: 80.50 YZ: 80.50 ZI: 80.50 YX: 80.50 XX: 49.55 ZY: 49.55 ... 0.50 0.50j 0.50 0.50j 0.50 0.50j 0.50 0.50j 0.50 0.50j 0.50 0.50j 0.50 0.50j 0.50 0.50j Labels: (2,) Name: Gate.cliff13 Aliases: Gate.cliff13 Generators: X: 127.28 Z: 127.28 0.71 0.71 0.71 0.71 13 Labels: (4, 9) Name: Gate(IY, XY, ...) Locally Equivalent: CNOT Generators: ZY: 127.28 XY: 127.28 YZ: 69.28 YX: 69.28 IY: 69.28 0.71 0.71 0.71 0.71 0.71 0.71 0.71 0.71 Labels: (7,) Name: Gate.cliff10 Aliases: Gate.cliff10 Generators: X: 127.28 Y: 127.28 0.71 0.71j 0.71 0.71j 10 Imm Labels: (0,) Name: Meas() M Labels: (1,) Name: Meas() M Labels: (2,) Name: Meas() M Labels: (4,) Name: Meas() M Labels: (7,) Name: Meas() M Labels: (9,) Name: Meas() M Note
We invert each random sequence up to an independent random Pauli matrix (accessible as shown in the example above) to diagnose errors to that are missed by always returning to the initial state (see arXiv:1901.00535).
 Parameters
labels (
Iterable
) – A list of which sets of system labels are to be twirled together in each circuit, e.g.[3, [1, 2], 4]
.n_random_cycles (
Iterable
) – A list of positive integers specifying how many random cycles will be generated during the protocol, e.g.[6, 20]
.n_circuits (
int
) – The number of circuits for each number of random cycles.twirling_group (
Iterable
) – The twirling group to use in this protocol. This can be either"C"
(Clifford) or"SU"
(special unitary). The twirling group will be the tensor product of the specified group on the labels to be benchmarked. A mixture of groups on different labels can be specified by an iterable of groups and labels such as[["SU", 0], ["C", 1, 2], ["C", 3]]
, where["C", 1, 2]
specifies the twoqubit Clifford group acting on the labels[1, 2]
and["C", 3]
specifies the onequbit Clifford group acting on the label3
. Any unspecified labels will be set to"C"
.
 Returns
A collection of SRB circuits.
 Return type
Make XRB¶

trueq.
make_xrb
(labels, n_random_cycles, n_circuits=30, twirling_group=None)¶ Generates a
CircuitCollection
to study how coherent the errors in random gates are using (simultaneous) extended randomized benchmarking (XRB). See also the XRB guide.import trueq as tq # generate circuit collections to run XRB on a single qubit (qubit 0), with # 30 random circuits for each length in [4, 50, 500] circuits = tq.make_xrb([0], [4, 50, 500], 30) # generate circuit collections to run twoqubit XRB on qubit pair [5, 6], with # 50 random circuits for each length in [3, 150] circuits = tq.make_xrb([[5, 6]], [3, 150], 50) # generate a circuit collection to run simultaneous onequbit XRB on [5, 6, 7] # with 30 random circuits for each length in [5, 10, 100] circuits = tq.make_xrb([5, 6, 7], [5, 10, 100], 30) # generate a circuit collection to run simultaneous onequbit XRB on qubit 5, # and twoqubit XRB on qubits [1, 2], with 15 random circuits for each circuit # length in [4, 30] circuits = tq.make_xrb([[1, 2], 5], [4, 30], 15) # draw the first circuit circuits[0].draw()
1 2 5 Key: twirl: (('C', 5), ('C', 1, 2)) protocol: XRB seq_label: 4767 n_random_cycles: 4 measurement_basis: XYY Imm Labels: (1, 2) Name: Gate(IX, IY, ...) Locally Equivalent: iSWAP Generators: ZY: 60.00 XI: 60.00 YZ: 60.00 XZ: 60.00 XX: 60.00 IZ: 60.00 IY: 34.64 ... 0.50 0.50j 0.50 0.50j 0.50 0.50j 0.50 0.50j 0.50 0.50j 0.50 0.50j 0.50 0.50j 0.50 0.50j Labels: (5,) Name: Gate.cliff10 Aliases: Gate.cliff10 Generators: X: 127.28 Y: 127.28 0.71 0.71j 0.71 0.71j 10 Imm Labels: (1, 2) Name: Gate(XI, YX, ...) Locally Equivalent: CNOT Generators: YX: 69.28 ZX: 69.28 XI: 69.28 0.50j 0.50 0.50 0.50j 0.50 0.50j 0.50j 0.50 0.50 0.50j 0.50j 0.50 0.50j 0.50 0.50 0.50j Labels: (5,) Name: Gate.id Aliases: Gate.id Gate.i Gate.cliff0 Locally Equivalent: Identity Generators: I: 0.00 1.00 1.00 ID Imm Labels: (1, 2) Name: Gate(IX, XI, ...) Locally Equivalent: iSWAP Generators: XX: 90.00 IX: 63.64 YZ: 63.64 XI: 63.64 ZY: 63.64 ZZ: 45.00 YY: 45.00 0.50 0.50j 0.50j 0.50 0.50 0.50j 0.50j 0.50 0.50 0.50j 0.50j 0.50 0.50 0.50j 0.50j 0.50 Labels: (5,) Name: Gate.z Aliases: Gate.z Gate.cliff3 Generators: Z: 180.00 1.00 1.00 Z Imm Labels: (1, 2) Name: Gate(IX, XY, ...) Locally Equivalent: SWAP Generators: ZX: 90.00 XZ: 63.64 YZ: 63.64 YY: 63.64 XY: 63.64 IX: 45.00 ZI: 45.00 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j 0.35 0.35j Labels: (5,) Name: Gate.cliff10 Aliases: Gate.cliff10 Generators: X: 127.28 Y: 127.28 0.71 0.71j 0.71 0.71j 10 Imm Labels: (1,) Name: Meas() M Labels: (2,) Name: Meas() M Labels: (5,) Name: Meas() M Note
If the circuit collection also contains SRB circuits with the same twirl group, then calling
fit
will also return an estimatee_U
of how much the process fidelity of the twirling groups can be improved by correcting static calibration errors. Parameters
labels (
Iterable
) – A list of which sets of system labels are to be twirled together in each circuit, e.g.[3, [1, 2], 4]
.n_random_cycles (
Iterable
) – A list of positive integers specifying how many random cycles will be generated during the protocol, e.g.[6, 20]
.n_circuits (
int
) – The number of circuits for each number of random cycles.twirling_group (
Iterable
) – The twirling group to use in this protocol. This can be either"C"
(Clifford) or"SU"
(special unitary). The twirling group will be the tensor product of the specified group on the labels to be benchmarked. A mixture of groups on different labels can be specified by an iterable of groups and labels such as[["SU", 0], ["C", 1, 2], ["C", 3]]
, where["C", 1, 2]
specifies the twoqubit Clifford group acting on the labels[1, 2]
and["C", 3]
specifies the onequbit Clifford group acting on the label3
. Any unspecified labels will be set to"C"
.
 Returns
A collection of XRB circuits.
 Return type
QCAP (Quantum Capacity)¶

trueq.
make_qcap
(cycles, n_random_cycles, n_circuits=30, n_decays=20)¶ Generates a
CircuitCollection
to measure the quantum capacity (QCAP) bound of any circuit whose immutable cycles are contained incycles
. See also the QCAP guide.The bound can be obtained by calling
qcap_bound()
on the QCAP circuits once they have been performed and their results have been populated. Note that this function outputs cycle benchmarking circuits, andfit()
will report the individual infidelities of the given cycles.import trueq as tq # make a circuit to assess circuit = tq.Circuit([ {0: tq.Gate.y}, {(0, 1): tq.Gate.cnot}, {0: tq.Gate.h, 1: tq.Gate.h}, {(1, 2): tq.Gate.cnot, (3, 4): tq.Gate.cnot}, {0: tq.Gate.h, 1: tq.Gate.h}, {(1, 2): tq.Gate.cnot, (3, 4): tq.Gate.cnot}, {0: tq.Gate.h, 1: tq.Gate.h}, {(1, 3): tq.Gate.cnot} ]) # generate a circuit collection to measure the QCAP circuits = tq.make_qcap(circuit, [0, 16]) # run the circuits on your hardware or simulator tq.Simulator().add_overrotation(0.05).run(circuits) # generate the QCAP bound tq.qcap_bound(circuit, circuits)
TrueQ formatting will not be loaded without trusting this notebook or rerunning the affected cells. Notebooks can be marked as trusted by clicking "File > Trust Notebook".QCAPQuantum CapacityPaulis (0, 1, 2, 3, 4) Key:

 labels: (0, 1, 2, 3, 4)
 protocol: QCAP
 twirl: (('P', 0), ('P', 1), ('P', 2), ('P', 3), ('P', 4))
${e}_{IU}$The inferred upper bound on the process infidelity of a circuit if it were run under RC.1.1e01 (3.0e03) 0.11159151234687037, 0.003047013470311924Note
This function is a convenience wrapper of
make_cb()
. If, for example, different values ofn_random_cycles
are desired for each cycle, callmake_cb()
manually for each cycle placing all generated circuits in the same collection. However, it is important to set thetwirling_group
option to((P, i) for i in labels)
wherelabels
are all the qubit labels involved in the circuit(s) of interest. Parameters
cycles (
Iterable
Circuit
) – An iterable ofCycle
s (note that aCircuit
is an iterable of cycles). Only the immutable cycles will be considered, and duplicate immutable cycles will be ignored.n_random_cycles (
Iterable
) – A list of positive integers specifying how many random cycles will be generated during the protocol, e.g.[6, 20]
.n_circuits (
int
) – The number of circuits for each random cycle.n_decays (
int
) – An integer specifying the total number of randomly chosen pauli decay strings used to measure the process infidelity. Warning: Setting this value lower thanmin(20, 4 ** n_qubits  1)
may result in a biased estimate.
 Returns
A collection of QCAP circuits.
 Return type

trueq.
qcap_bound
(circuit, fit_or_circuits)¶ Computes the quantum capacity (QCAP) bound of the given circuit. This scalar quantity is a bound on the circuit performance when the circuit is executed using
randomly_compile()
. See also the QCAP guide. Parameters
circuit (
Circuit
) – The circuit you wish to know the QCAP bound of. This will usually be the original circuit that was given tomake_qcap()
, but can be any circuit whoseimmutable
cycles are a subset of the immutable cycles present in the CB circuits or fits offit_or_circuits
.fit_or_circuits – A circuit collection output by
make_qcap()
populated with results, or an estimate collection that resulted from callingfit()
on such a circuit collection. For efficiency, use the latter format if computing a bound on multiple circuits.
 Returns
A bound on the circut performance.
 Return type
Randomly Compile¶

trueq.
randomly_compile
(circuit, n_compilations=30, twirling_group=None, compile_paulis=False, compress=True)¶ Randomly compiles the given circuit into many new random circuits which implement the same algorithm. Random gates are inserted adjacent to gates in the provided circuit, chosen from the specified twirling group (the Pauli group is used by default). See also the RC guide.
Cycles in the input circuit that are immutable will not be chnaged, while other cycles may have random gates compiled into them.
import trueq as tq #Define a circuit which applies an X gate on the 0th qubit, a controlledZ #gate on qubits 0 and 1, and another X gate on the 0th qubit. cycle1 = tq.Cycle({0: tq.Gate.x}) cycle2 = tq.Cycle({(0, 1): tq.Gate.cz}); circuit = tq.Circuit((cycle1, cycle2, cycle1)) circuit.measure_all() #Run randomized compiling on the circuit compiled_circuit = tq.randomly_compile(circuit) # draw the first circuit compiled_circuit[0].draw()
0 1 Key: twirl: (('P', 0), ('P', 1)) protocol: RC n_random_cycles: 1 Labels: (0,) Name: Gate.y Aliases: Gate.y Gate.cliff2 Generators: Y: 180.00 1.00 1.00 Y Labels: (1,) Name: Gate.id Aliases: Gate.id Gate.i Gate.cliff0 Locally Equivalent: Identity Generators: I: 0.00 1.00 1.00 ID Imm Labels: (0, 1) Name: Gate.cz Aliases: Gate.cz Locally Equivalent: CNOT Generators: IZ: 90.00 ZI: 90.00 ZZ: 90.00 1.00 1.00 1.00 1.00 CZ CZ Labels: (0,) Name: Gate.y Aliases: Gate.y Gate.cliff2 Generators: Y: 180.00 1.00 1.00 Y Labels: (1,) Name: Gate.id Aliases: Gate.id Gate.i Gate.cliff0 Locally Equivalent: Identity Generators: I: 0.00 1.00 1.00 ID Imm Labels: (0,) Name: Meas() M Labels: (1,) Name: Meas() M Note
Running randomly_compile() on a circuit several times is likely to return different compiled circuits due to the random nature of the algorithm. The circuit returned will always contain the same number of cycles as the input circuit, and will be logically equivalent, up to a global phase.
Warning
Running randomly_compile() on circuits containing gates acting on multiple systems with a nonClifford gate will typically result in an error because twirling operations cannot be corrected locally.
 Parameters
circuit (
Circuit
) – A circuit to randomly compile.n_compilations (
int
) – The number of random compilations of the input circuit. Each instance will appear as a new circuit in the returned circuit collection.twirling_group (
Iterable
) – The twirling group to use. This can be one of"P"
(Pauli),"C"
(Clifford),"SU"
(special unitary), or"I"
(Identity). The twirling group will be the tensor product of the specified group on thelabels
of the circuit to be randomly compiled. A mixture of groups on different labels can be specified by an iterable of groups and labels, for example,[["P", 0], ["C", 1, 2], ["C", 3]]
, where["C", 1, 2]
specifies the twoqubit Clifford group acting on the labels[1, 2]
and["C", 3]
specifies the onequbit Clifford group acting on the label3
. Any unspecified elements oflabels
will be set to"P"
.compile_paulis (
bool
) – Whether or not to compile a random Pauli gate onto a qubit in the cycle preceding a measurement operation for eachMeas
operation encountered in the circuit. Which Paulis were compiled in are stored as a string in thekey
, where the order is first defined by cycle index, and then by sorted labels of each cycle. (default isFalse
).compress (
bool
) – Whether or not to compress the easy cycles (default isTrue
).
 Returns
A circuit collection containing randomly compiled versions of the circuit.
 Return type
py:class:~trueq.CircuitCollection