User guide

Introduction

quantify-scheduler is a python module for writing quantum programs featuring a hybrid gate-pulse control model with explicit timing control. It extends the circuit model from quantum information processing by adding a pulse-level representation to operations defined at the gate-level, and the ability to specify timing constraints between operations. Thus, a user is able to mix gate- and pulse-level operations in a quantum circuit.

In quantify-scheduler, both a quantum circuit consisting of gates and measurements and a timed sequence of control pulses are described as a Schedule . The Schedule contains information on when operations should be performed. When adding operations to a schedule, one does not need to specify how to represent this Operation on all (both gate and pulse) abstraction levels. Instead, this information can be added later during Compilation. This allows the user to effortlessly mix the gate- and pulse-level descriptions as is required for many experiments. We support similar flexibility in the timing constraints, one can either explicitly specify the timing using ScheduleBase.schedulables, or rely on the compilation which will use the duration of operations to schedule them back-to-back.

Creating a schedule

The most convenient way to interact with a Schedule is through the quantify_scheduler API. In the following example, we will create a function to generate a Schedule for a Bell experiment and visualize one instance of such a circuit.

# import the Schedule class and some basic operations.
from quantify_scheduler import Schedule
from quantify_scheduler.operations.gate_library import Reset, Measure, CZ, Rxy, X90

def bell_schedule(angles, q0:str, q1:str, repetitions: int):

    for acq_index, angle in enumerate(angles):

        sched = Schedule(f"Bell experiment on {q0}-{q1}")

        sched.add(Reset(q0, q1))  # initialize the qubits
        sched.add(X90(qubit=q0))
        # Here we use a timing constraint to explicitly schedule the second gate to start
        # simultaneously with the first gate.
        sched.add(X90(qubit=q1), ref_pt="start", rel_time=0)
        sched.add(CZ(qC=q0, qT=q1))
        sched.add(Rxy(theta=angle, phi=0, qubit=q0) )
        sched.add(Measure(q0, acq_index=acq_index))  # denote where to store the data
        sched.add(Measure(q1, acq_index=acq_index), ref_pt="start")

    return sched


sched = bell_schedule(
    angles=[45.0],
    q0="q0",
    q1="q1",
    repetitions=1024)
# visualize the circuit
f, ax = sched.plot_circuit_diagram()

_images/a807fbd61eb6b260ab8b1b20e5046b192710db169ccac86d4d81c0768b30e845.png

Tip

Creating schedule generating functions is a convenient design pattern when creating measurement code. See the section on execution for an example of how this is used in practice.

Concepts and terminology

quantify-scheduler can be understood by understanding the following concepts.

  • Schedules describe when an operation needs to be applied.

  • Operations describe what needs to be done.

  • Resources describe where an operation should be applied.

  • Compilation: between different abstraction layers for execution on physical hardware.

The following table shows an overview of the different concepts and how these are represented at the quantum-circuit layer and quantum-device layer.

Table 1 Overview of concepts and their representation at different levels of abstraction.

Concept

Quantum-circuit layer

Quantum-device layer

When

Schedule

What

Operation

Gates and Measurements

Pulses and acquisition protocols

Where

Resource

Qubits

Ports and clocks

Quantum-circuit layer

The Quantum-circuit description is an idealized mathematical description of a schedule.

Gates and measurements

In this description operations are quantum gates that act on idealized qubits as part of a quantum circuit. Operations can be represented by (idealized) unitaries acting on qubits. The gate_library contains common operations (including the measurement operation) described at the quantum-circuit level.

The Measure is a special operation that represents a measurement on a qubit. In addition to the qubit it acts on, one also needs to specify where to store the data.

Qubits

At the gate-level description, operations are applied to qubits. Qubits are represented by strings corresponding to the name of a qubit (e.g., q0, q1, A1, QL, qubit_1, etc.). Valid qubits are strings that appear in the device configuration file used when compiling the schedule.

Visualization

A Schedule containing operations can be visualized using a circuit diagram by calling its method plot_circuit_diagram().

Alternatively, one can plot the waveforms in schedules using plot_pulse_diagram() (the default plotting backend is matplotlib, but it is possible to use plotly by adding the argument plot_backend='plotly'):

from quantify_scheduler.operations.pulse_library import SquarePulse, RampPulse
from quantify_scheduler.compilation import determine_absolute_timing

schedule = Schedule("waveforms")
schedule.add(SquarePulse(amp=0.2, duration=4e-6, port="P"))
schedule.add(RampPulse(amp=-0.1, offset=.2, duration=6e-6, port="P"))
schedule.add(SquarePulse(amp=0.1, duration=4e-6, port="Q"), ref_pt='start')
determine_absolute_timing(schedule)

_ = schedule.plot_pulse_diagram(sampling_rate=20e6)

_images/737dadc28df781c3d3107f8f21abae8ce02f75489668260d15aa277b552653eb.png

Summary

  • Gates are described by unitaries.

  • Gates are applied to qubits.

  • Measurements are applied to qubits.

  • Qubits are represented by strings.

Quantum-device layer

The quantum-device layer describes waveforms and acquisition protocols applied to a device. These waveforms can be used to implement the idealized operations expressed on the quantum-circuit layer, or can be used without specifying a corresponding representation at the quantum-circuit layer.

Pulses and acquisition operations

The pulse-level description typically contains parameterization information, such as amplitudes, durations and so forth required to synthesize the waveform on control hardware. pulse_library module contains a collection of commonly used pulses.

Measurements are decomposed into pulses and acquisition operations. Similarly to pulse operations, acquisition operations contain their timing information and correspondent acquisition protocols. acquisition_library module contains a collection of commonly used acquisition operations.

Ports and clocks

To specify where an operation is applied, the quantum-device layer description needs to specify both the location in physical space as well as in frequency space.

For many systems, it is possible to associate a qubit with an element or location on a device that a signal can be applied to. We call such a location on a device a port. Like qubits, ports are represented as strings (e.g., P0, feedline_in, q0:mw_drive, etc.). In the last example, a port is associated with a qubit by including the qubit name at the beginning of the port name (separated by a colon :).

Associating a qubit can be useful when visualizing a schedule and or keeping configuration files readable. It is, however, not required to associate a port with a single qubit. This keeps matters simple when ports are associated with multiple qubits or with non-qubit elements such as tunable couplers.

Besides the physical location on a device, a pulse is typically applied at a certain frequency and with a phase. These two parameters are stored in a ClockResource. Each ClockResource also has a name to be easily identified. The name should identify the purpose of the clock resource, not the value of the frequency. By storing the frequency and phase in a clock, we can adjust the frequency of a transition, but refer to it with the same name.

Similar to ports, clocks can be associated with qubits by including the qubit name in the clock name (again, this is not required). If the frequency of a clock is set to 0 (zero), the pulse is applied at baseband and is assumed to be real-valued.

Fig. 1 shows how the resources (qubit, port and clock) map to a physical device.

_images/Device_ports_clocks.svg

Fig. 1 Resources are used to indicate where operations are applied. (a) Ports (purple) indicate a location on a device. By prefixing the name of a qubit in a port name (separated by a colon :) a port can be associated with a qubit (red), but this is not required. (b) Clocks (blue) denote the frequency and phase of a signal. They can be set to track the phase of a known transition. By prefixing the name of a qubit in a clock name (separated by a colon :) a clock can be associated with a qubit (red), but this is not required. Device image from Dickel (2018) .

Acquisition protocols

When we define an acquisition, we must specify how to acquire data and how to process it to provide a meaningful result. For example, typical readout of a superconducting qubit state will consist of the following steps:

  1. Sending a spectroscopy pulse to a port of a device, that is connected to a readout resonator of the qubit.

  2. Acquiring a raw voltage trace of the reflected or transmitted signal using microwave digitizer equipment.

  3. Taking a weighted sum of the returned signal to obtain a single complex number.

  4. Assign most likely readout outcome (0 or 1) based on a pre-calibrated threshold.

Description of these processing steps is called an acquisition protocol. To define an acquisition protocol, the developer must define two things: an algorithm to process the data and the schema of a data array returned by an instrument coordinator component (ICC). ICC’s job is to retrieve data from the instrument, process it and return it in a required format. If the readout equipment supports hardware acceleration for part of the processing steps, its ICC can and should utilize it.

The Acquisition data format chapter briefly describes the implication of acquisition protocols and binning mode on a format of the data returned by ICC. For a detailed description of all acquisition protocols defined in quantify-scheduler and supported by at least some of the backends consult the acquisition protocols reference section.

Acquisition channel and acquisition index

quantify-scheduler identifies each acquisition in a resulting dataset with two integer numbers: acquisition channel and acquisition index.

Acquisition channel is a stream of acquisition data that corresponds to a single device element measured with the same acquisition protocol. Each acquisition within the same acquisition channel must return data of a uniform size and structure, that is described by the definition of an acquisition protocol. The order number of an acquisition within the acquisition channel during a schedule run is called an acquisition index.

On a quantum circuit layer they roughly correspond to a qubit being measured and a number of the measurement of a given qubit in a schedule. However, if a qubit (or, more precisely, a device element) can be measured using more than one acquisition protocol, that will require defining several acquisition channels for it. When you are specifying a schedule on circuit level, acquisition channel is supposed to be configured within a submodule of a device element, that corresponds to a given measurement (i.e. measure() submodule of a BasicTransmonElement}). Acquisition index should be specified at the instantiation of a measurement operation.

Summary

  • Pulses are described as parameterized waveforms.

  • Pulses are applied to ports at a frequency specified by a clock.

  • Ports and clocks are represented by strings.

  • Acquisition protocols describe the processing steps to perform on an acquired signal in order to interpret it.

  • Acquisition channel and acquisition index describe how to find an acquisition result in a dataset.

Compilation

Different compilation steps are required to go from a high-level description of a schedule to something that can be executed on hardware. The scheduler supports multiple compilation steps, the most important ones are the step from the quantum-circuit layer (gates) to the quantum-device layer (pulses), and the one from the quantum-device layer into instructions suitable for execution on physical hardware. The compilation is performed by a QuantifyCompiler, which is configured through the CompilationConfig:

class CompilationConfig(**data)[source]

Base class for a compilation config.

Subclassing is generally required to create useful compilation configs, here extra fields can be defined.

backend: Type[QuantifyCompiler][source]

A reference string to the QuantifyCompiler class used in the compilation.

device_compilation_config: Union[DeviceCompilationConfig, Dict, None][source]

The DeviceCompilationConfig used in the compilation from the quantum-circuit layer to the quantum-device layer.

hardware_compilation_config: Optional[HardwareCompilationConfig][source]

The HardwareCompilationConfig used in the compilation from the quantum-device layer to the control-hardware layer.

The compilers are described in detail in Compilers. This is schematically shown in Fig. 2.

_images/compilation_overview.svg

Fig. 2 A schematic overview of the different abstraction layers and the compilation process. Both a quantum circuit, consisting of gates and measurements of qubits, and timed sequences of control pulses are represented as a Schedule . The information specified in the device configuration is used during compilation to add information on how to represent Operation s specified at the quantum-circuit level as control pulses. The information in the hardware description, hardware options, and connectivity is then used to compile the control pulses into instructions suitable for hardware execution.

In the first compilation step, pulse information is added to all operations that are not valid pulses (see Operation.valid_pulse) based on the information specified in the Device compilation configuration.

A second compilation step takes the schedule at the pulse level and translates this for use on a hardware backend. This compilation step is performed using a hardware dependent compiler and uses the information specified in the Hardware compilation configuration.

Note

We use the term “device” to refer to the physical object(s) on the receiving end of the control pulses, e.g. a thin-film chip inside a dilution refrigerator.

And we employ the term “hardware” to refer to the instruments (electronics) that are involved in the pulse generations / signal digitization.

Device compilation configuration

The device compilation configuration is used to compile from the quantum-circuit layer to the quantum-device layer. This datastructure is auto-generated by the QuantumDevice using the parameters stored in the DeviceElements.

class DeviceCompilationConfig(**data)[source]

Information required to compile a schedule to the quantum-device layer.

From a point of view of Compilation this information is needed to convert a schedule defined on a quantum-circuit layer to a schedule defined on a quantum-device layer.

Hardware compilation configuration

The hardware compilation configuration is used to compile from the quantum-device layer to the control-hardware layer. Currently, this datastructure is parsed from a user-defined dict that should be passed to the quantum_device.hardware_config parameter.

class HardwareCompilationConfig(**data)[source]

Information required to compile a schedule to the control-hardware layer.

From a point of view of Compilation this information is needed to convert a schedule defined on a quantum-device layer to compiled instructions that can be executed on the control hardware.

backend: Callable[[Schedule, Any], Schedule][source]

A . separated string specifying the location of the compilation backend this configuration is intended for.

connectivity: Union[Connectivity, Dict][source]

Datastructure representing how ports on the quantum device are connected to ports on the control hardware.

hardware_description: Dict[str, Union[ClusterDescription, PulsarQCMDescription, PulsarQRMDescription, ZIHDAWG4Description, ZIHDAWG8Description, ZIUHFQADescription, LocalOscillatorDescription]][source]

Datastructure describing the control hardware instruments in the setup and their high-level settings.

hardware_options: Optional[HardwareOptions][source]

The HardwareOptions used in the compilation from the quantum-device layer to the control-hardware layer.

Hardware Description

The HardwareDescription datastructure specifies a control hardware instrument in the setup, along with its instrument-specific settings. There is a specific HardwareDescription datastructure for each of the currently supported instrument types, which are discriminated through the hardware_type and instrument_type fields.

Hardware Options

The HardwareOptions datastructure contains the settings used in compiling from the quantum-device layer to a set of instructions for the control hardware. Most hardware options are structured as Dict[str, HardwareOption], where the keys are the port-clock combinations on which these settings should be applied.

class HardwareOptions(**data)[source]

Datastructure containing the hardware options for each port-clock combination.

distortion_corrections: Optional[Dict[str, DistortionCorrection]][source]

Dictionary containing the distortion corrections (values) that should be applied to waveforms on a certain port-clock combination (keys).

latency_corrections: Optional[Dict[str, LatencyCorrection]][source]

Dictionary containing the latency corrections (values) that should be applied to operations on a certain port-clock combination (keys).

mixer_corrections: Optional[Dict[str, MixerCorrections]][source]

Dictionary containing the mixer corrections (values) that should be used for signals on a certain port-clock combination (keys).

modulation_frequencies: Optional[Dict[str, ModulationFrequencies]][source]

Dictionary containing the modulation frequencies (values) that should be used for signals on a certain port-clock combination (keys).

power_scaling: Optional[Dict[str, PowerScaling]][source]

Dictionary containing the gain/attenuation settings (values) that should be used for signals on a certain port-clock combination (keys).

Connectivity

The Connectivity datastructure describes how ports on the quantum device are connected to ports on the control hardware.

Note

The Connectivity datastructure is currently under development. Information on the connectivity between port-clock combinations on the quantum device and ports on the control hardware is currently included in the old-style hardware configuration file, which should be included in the "connectivity" field of the HardwareCompilationConfig.

Execution

Warning

This section describes functionality that is not fully implemented yet. The documentation describes the intended design and may change as the functionality is added.

Different kinds of instruments

In order to execute a schedule, one needs both physical instruments to execute the compiled instructions as well as a way to manage the calibration parameters used to compile the schedule. Although one could use manually written configuration files and send the compiled files directly to the hardware, the Quantify framework provides different kinds of Instruments to control the experiments and the management of the configuration files (Fig. 3).

_images/instruments_overview.svg

Fig. 3 A schematic overview of the different kinds of instruments present in an experiment. Physical instruments are QCoDeS drivers that are directly responsible for executing commands on the control hardware. On top of the physical instruments is a hardware abstraction layer, that provides a hardware agnostic interface to execute compiled schedules. The instruments responsible for experiment control are treated to be as stateless as possible 1 . The knowledge about the system that is required to generate the configuration files is described by the QuantumDevice and DeviceElements. Several utility instruments are used to control the flow of the experiments.

Physical instruments

QCoDeS instrument drivers are used to represent the physical hardware. For the purpose of quantify-scheduler, these instruments are treated as stateless, the desired configurations for an experiment being described by the compiled instructions. Because the instruments correspond to physical hardware, there is a significant overhead in querying and configuring these parameters. As such, the state of the instruments in the software is intended to track the state of the physical hardware to facilitate lazy configuration and logging purposes.

Hardware abstraction layer

Because different physical instruments have different interfaces, a hardware abstraction layer serves to provide a uniform interface. This hardware abstraction layer is implemented as the InstrumentCoordinator to which individual InstrumentCoordinatorComponents are added that provide the uniform interface to the individual instruments.

The quantum device and the device elements

The knowledge of the system is described by the QuantumDevice and DeviceElements. The QuantumDevice directly represents the device under test (DUT) and contains a description of the connectivity to the control hardware as well as parameters specifying quantities like cross talk, attenuation and calibrated cable-delays. The QuantumDevice also contains references to individual DeviceElements, representations of elements on a device (e.g, a transmon qubit) containing the (calibrated) control-pulse parameters.

Because the QuantumDevice and the DeviceElements are an Instrument, the parameters used to generate the configuration files can be easily managed and are stored in the snapshot containing the experiment’s metadata.

Experiment flow

To use schedules in an experimental setting, in which the parameters used for compilation as well as the schedules themselves routinely change, we provide a framework for performing experiments making use of the concepts of quantify-core. Central in this framework are the schedule quantify_scheduler.gettables that can be used by the quantify_core.measurement.control.MeasurementControl and are responsible for the experiment flow.

This flow is schematically shown in Fig. 4.

_images/experiments_control_flow.svg

Fig. 4 A schematic overview of the experiments control flow.

Let us consider the example of an experiment used to measure the coherence time \(T_1\). In this experiment, a \(\pi\) pulse is used to excite the qubit, which is left to idle for a time \(\tau\) before it is measured. This experiment is then repeated for different \(\tau\) and averaged.

In terms of settables and gettables to use with the quantify_core.measurement.control.MeasurementControl, the settable in this experiment is the delay time \(\tau\), and the gettable is the execution of the schedule.

We represent the settable as a qcodes.instrument.parameter.ManualParameter:

from qcodes.instrument.parameter import ManualParameter

tau = ManualParameter("tau", label=r"Delay time", initial_value=0, unit="s")

To execute the schedule with the right parameters, the ScheduleGettable needs to have a reference to a template function that generates the schedule, the appropriate keyword arguments for that function, and a reference to the QuantumDevice to generate the required configuration files.

For the \(T_1\) experiment, quantify-scheduler provides a schedule generating function as part of the quantify_scheduler.schedules.timedomain_schedules: the quantify_scheduler.schedules.timedomain_schedules.t1_sched().

from quantify_scheduler.schedules.timedomain_schedules import t1_sched
schedule_function = t1_sched

Inspecting the quantify_scheduler.schedules.timedomain_schedules.t1_sched(), we find that we need to provide the times \(\tau\), the name of the qubit, and the number of times we want to repeat the schedule. Rather than specifying the values of the delay times, we pass the parameter tau.

qubit_name = "q0"
sched_kwargs = {
    "times": tau,
    "qubit": qubit_name,
    "repetitions": 1024 # could also be a parameter
}

The ScheduleGettable is set up to evaluate the value of these parameter on every call of ScheduleGettable.get. This flexibility allows the user to create template schedules that can then be measured by varying any of it’s input parameters using the quantify_core.measurement.control.MeasurementControl.

Similar to how the schedule keyword arguments are evaluated for every call to ScheduleGettable.get, the device config and hardware config files are re-generated from the QuantumDevice for every iteration. This ensures that if a calibration parameter is changed on the QuantumDevice, the compilation will be affected as expected.

from quantify_scheduler.device_under_test.quantum_device import QuantumDevice
device = QuantumDevice(name="quantum_sample")

These ingredients can then be combined to perform the experiment:

from quantify_core.measurement import MeasurementControl
meas_ctrl = MeasurementControl("meas_ctrl")
t1_gettable = ScheduleGettable(
    device=device,
    schedule_function=schedule_function,
    schedule_kwargs=sched_kwargs
)

meas_ctrl.settables(tau)
meas_ctrl.setpoints(times)
meas_ctrl.gettables(t1_gettable)
label = f"T1 experiment {qubit_name}"
dataset = meas_ctrl.run(label)

and the resulting dataset can be analyzed using

# from quantify_core.analysis.t1_analysis import T1Analysis
# analysis = T1Analysis(label=label).run()

Acquisition data format

quantify-scheduler has two primary interfaces for retrieving acquisition results: using InstrumentCoordinator and using quantify_scheduler.gettables.ScheduleGettable.

Each acquisition and measurement operation in the schedule has associated acquisition channel and acquisition index. If you specifiy a schedule using raw acquisition operations (for example, using SSBIntegrationComplex), use acq_channel and acq_index arguments of the operation to specify them. Optionally you may also specify the requested binning mode in bin_mode (it is AVERAGE by default):

schedule.add(
    SSBIntegrationComplex(
        t0=0,
        duration=100e-9,
        port="q0:res",
        clock="q0.ro",
        acq_channel=3,
        acq_index=1,
        bin_mode=BinMode.AVERAGE
    )
)

For a selected acquisition channel acquisition protocol and binning mode must be the same, otherwise compilation will fail.

When circuit-to-device compilation machinery is used, acq_channel should be specified in the DeviceElement being measured, for example in the DispersiveMeasurement submodule of BasicTransmonElement. acq_index and bin_mode are still specified as input parameters to Measure (or another specialised measurement operation supported by the device element).

Retrieve acquisitions through InstrumentCoordinator

InstrumentCoordinator.retrieve_acquisition() method returns an xarray.Dataset:

  • Each xarray.DataArray in the dataset corresponds to one acq_channel.

  • Exact structure of a data array is defined by an acquisition protocol and bin_mode, that are associated with this acquisition channel.

For example, if a schedule contains two qubits (each one has its own acquisition channel, say, 0 and 1), the first of which has been measured three times and second twice using SSBIntegrationComplex acquisition protocol in BinMode.APPEND, the resulting dataset will have the form:


<xarray.Dataset>
Dimensions:  (repetition: 5, acq_index_0: 3, acq_index_1: 2)
Dimensions without coordinates: repetition, acq_index_0, acq_index_1
Data variables:
    0        (repetition, acq_index_0) complex128 (0.948+0.402j) ... (0.185+0...
    1        (repetition, acq_index_1) complex128 (0.007+0.763j) ... (0.292+0...

Definitions of acquisition protocols and correspondent data format can be found in the acquisition protocols reference section. Note that acquisition protocols define the meaning of each dimension of a data array, but do not guarantee a consistent naming of the dimensions in a dataset. Instead, the exact names of dimensions should be retrieved dynamically during the processing of the dataset.

Retrieve acquisition through ScheduleGettable

ScheduleGettable proxies the instrument-coordinator format to a format that can be used with MeasurementControl. Effectively it flattens the data arrays retrieved from the instrument coordinator, splits complex numbers into either real and imaginary parts (if real_imag is set to True) or absolute value and phase (if real_imag is False) and pads the data with nans to fit it into a single array. Data that corresponds to acquisition channel number \(n\) will end up in items number \(2n\) and \(2n+1\) of that array. For example, if real_imag is set to True in the ScheduleGettable, the dataset above will be converted to:


(
    array([0.948, 0.3  , 0.765, 0.944, 0.35 , 0.417, 0.851, 0.351, 0.713,
       0.982, 0.531, 0.57 , 0.487, 0.215, 0.185]),
    array([0.402, 0.046, 0.695, 0.57 , 0.23 , 0.82 , 0.765, 0.679, 0.093,
       0.048, 0.344, 0.529, 0.093, 0.163, 0.499]),
    array([0.007, 0.494, 0.271, 0.9  , 0.277, 0.178, 0.482, 0.036, 0.317,
       0.292]),
    array([0.763, 0.83 , 0.412, 0.089, 0.779, 0.274, 0.071, 0.23 , 0.166,
       0.035])
)

Footnotes


1

quantify-scheduler treats physical instruments as stateless in the sense that the compiled instructions contain all information that specifies the execution of a schedule. However, for performance reasons, it is important to not reconfigure all parameters of all instruments whenever a new schedule is executed. The parameters (state) of the instruments are used to track the state of physical instruments to allow lazy configuration as well as ensure metadata containing the current settings is stored correctly.