See also

A Jupyter notebook version of this tutorial can be downloaded here.

Single transmon qubit spectroscopy

Here we will carry out qubit spectroscopy on a single transmon in order to find the \(|0\rangle \rightarrow |1\rangle\) drive frequency.

from __future__ import annotations

import json

import numpy as np
import rich  # noqa:F401
from qcodes.instrument import find_or_create_instrument
from qcodes.parameters import ManualParameter

import quantify_core.data.handling as dh
from qblox_instruments import Cluster, ClusterType
from quantify_core.measurement.control import MeasurementControl
from quantify_core.visualization.pyqt_plotmon import PlotMonitor_pyqt as PlotMonitor
from quantify_scheduler.device_under_test.quantum_device import QuantumDevice
from quantify_scheduler.gettables import ScheduleGettable
from quantify_scheduler.instrument_coordinator import InstrumentCoordinator
from quantify_scheduler.instrument_coordinator.components.qblox import (
    ClusterComponent,
)
from quantify_scheduler.operations.gate_library import Measure, Reset, X
from quantify_scheduler.schedules.schedule import Schedule

Setup

In this section we configure the hardware configuration which specifies the connectivity of our system.

The experiments of this tutorial are meant to be executed with a Qblox Cluster controlling a transmon system. The experiments can also be executed using a dummy Qblox device that is created via an instance of the Cluster class, and is initialized with a dummy configuration. When using a dummy device, the analysis will not work because the experiments will return np.nan values.

Configuration file

This is a template hardware configuration file for a 2-qubit system with a flux-control line which can be used to tune the qubit frequency. We will only work with qubit 0.

The hardware setup is as follows, by cluster slot: - QCM (Slot 2) - Flux line for q0. - QCM-RF (Slot 6) - Drive line for q0 using fixed 80 MHz IF. - QRM-RF (Slot 8) - Readout line for q0 using a fixed LO set at 7.5 GHz.

Note that in the hardware configuration below the mixers are uncorrected, but for high fidelity experiments this should also be done for all the modules.

with open("configs/tuning_transmon_coupled_pair_hardware_config.json") as hw_cfg_json_file:
    hardware_cfg = json.load(hw_cfg_json_file)

# Enter your own dataset directory here!
dh.set_datadir(dh.default_datadir())

Data will be saved in:
/root/quantify-data

Scan For Clusters

We scan for the available devices connected via ethernet using the Plug & Play functionality of the Qblox Instruments package (see Plug & Play for more info).

!qblox-pnp list

No devices found

cluster_ip = None  # To run this tutorial on hardware, fill in the IP address of the cluster here
cluster_name = "cluster0"

Connect to Cluster

We now make a connection with the Cluster.

cluster = find_or_create_instrument(
    Cluster,
    recreate=True,
    name=cluster_name,
    identifier=cluster_ip,
    dummy_cfg=(
        {
            2: ClusterType.CLUSTER_QCM,
            4: ClusterType.CLUSTER_QRM,
            6: ClusterType.CLUSTER_QCM_RF,
            8: ClusterType.CLUSTER_QRM_RF,
        }
        if cluster_ip is None
        else None
    ),
)

Quantum device settings

Here we initialize our QuantumDevice and our qubit parameters, checkout this tutorial for further details.

In short, a QuantumDevice contains device elements where we save our found parameters. Here we are loading a template for 2 qubits, but we will only use qubit 0.

quantum_device = QuantumDevice.from_json_file("devices/transmon_device_2q.json")
qubit = quantum_device.get_element("q0")
quantum_device.hardware_config(hardware_cfg)

Configure measurement control loop

We will use a MeasurementControl object for data acquisition as well as an InstrumentCoordinator for controlling the instruments in our setup.

The PlotMonitor is used for live plotting.

All of these are then associated with the QuantumDevice.

def configure_measurement_control_loop(
    device: QuantumDevice, cluster: Cluster, live_plotting: bool = False
) -> tuple[MeasurementControl, InstrumentCoordinator]:
    meas_ctrl = find_or_create_instrument(MeasurementControl, recreate=True, name="meas_ctrl")
    ic = find_or_create_instrument(InstrumentCoordinator, recreate=True, name="ic")

    # Add cluster to instrument coordinator
    ic_cluster = ClusterComponent(cluster)
    ic.add_component(ic_cluster)

    if live_plotting:
        # Associate plot monitor with measurement controller
        plotmon = find_or_create_instrument(PlotMonitor, recreate=False, name="PlotMonitor")
        meas_ctrl.instr_plotmon(plotmon.name)

    # Associate measurement controller and instrument coordinator with the quantum device
    device.instr_measurement_control(meas_ctrl.name)
    device.instr_instrument_coordinator(ic.name)

    return (meas_ctrl, ic)


meas_ctrl, instrument_coordinator = configure_measurement_control_loop(quantum_device, cluster)

Qubit spectroscopy

def two_tone_spec_sched_nco(
    qubit,  # noqa: ANN001
    spec_pulse_frequencies: np.array,
    repetitions: int = 1,
) -> Schedule:
    """
    Generate a batched schedule for performing fast two-tone spectroscopy.

    Using the X gate to perform the frequency sweep on the qubit.

    Parameters
    ----------
    qubit
        qubit that should be used.
    spec_pulse_frequencies
        Sample frequencies for the spectroscopy pulse in Hertz.
    repetitions
        The amount of times the Schedule will be repeated.

    """
    sched = Schedule("two-tone", repetitions)

    for acq_idx, spec_pulse_freq in enumerate(spec_pulse_frequencies):
        sched.add(Reset(qubit.name))
        sched.add(X(qubit.name, freq=spec_pulse_freq))
        sched.add(Measure(qubit.name, acq_index=acq_idx), rel_time=200e-9)
    return sched

freqs = ManualParameter(name="freq", unit="Hz", label="Frequency")
freqs.batched = True

qubit_spec_sched_kwargs = dict(
    qubit=qubit,
    spec_pulse_frequencies=freqs,
)

gettable = ScheduleGettable(
    quantum_device,
    schedule_function=two_tone_spec_sched_nco,
    schedule_kwargs=qubit_spec_sched_kwargs,
    real_imag=False,
    batched=True,
)

meas_ctrl.gettables(gettable)

quantum_device.cfg_sched_repetitions(300)
center = 6.1e9
frequency_setpoints = np.linspace(center - 20e6, center + 20e6, 300)
meas_ctrl.settables(freqs)
meas_ctrl.setpoints(frequency_setpoints)

qs_ds = meas_ctrl.run("Two-tone")
qs_ds

Starting batched measurement...
Iterative settable(s) [outer loop(s)]:
         --- (None) ---
Batched settable(s):
         freq
Batch size limit: 300

/usr/local/lib/python3.9/site-packages/quantify_scheduler/backends/types/qblox.py:1220: ValidationWarning: Setting `auto_lo_cal=on_lo_interm_freq_change` will overwrite settings `dc_offset_i=0.0` and `dc_offset_q=0.0`. To suppress this warning, do not set either `dc_offset_i` or `dc_offset_q` for this port-clock.
  warnings.warn(
/usr/local/lib/python3.9/site-packages/quantify_scheduler/backends/types/qblox.py:1235: ValidationWarning: Setting `auto_sideband_cal=on_interm_freq_change` will overwrite settings `amp_ratio=1.0` and `phase_error=0.0`. To suppress this warning, do not set either `amp_ratio` or `phase_error` for this port-clock.
  warnings.warn(

<xarray.Dataset> Size: 7kB
Dimensions:  (dim_0: 300)
Coordinates:
    x0       (dim_0) float64 2kB 6.08e+09 6.08e+09 ... 6.12e+09 6.12e+09
Dimensions without coordinates: dim_0
Data variables:
    y0       (dim_0) float64 2kB nan nan nan nan nan nan ... nan nan nan nan nan
    y1       (dim_0) float64 2kB nan nan nan nan nan nan ... nan nan nan nan nan
Attributes:
    tuid:                             20240902-014501-705-e85dbc
    name:                             Two-tone
    grid_2d:                          False
    grid_2d_uniformly_spaced:         False
    1d_2_settables_uniformly_spaced:  False

xarray.Dataset

• Dimensions:
  • dim_0: 300
• Coordinates: (1)
  • x0
    (dim_0)
    float64
    6.08e+09 6.08e+09 ... 6.12e+09

    name :

    freq

    long_name :

    Frequency

    units :

    Hz

    batched :

    True

    array([6.08000000e+09, 6.08013378e+09, 6.08026756e+09, 6.08040134e+09,
           6.08053512e+09, 6.08066890e+09, 6.08080268e+09, 6.08093645e+09,
           6.08107023e+09, 6.08120401e+09, 6.08133779e+09, 6.08147157e+09,
           6.08160535e+09, 6.08173913e+09, 6.08187291e+09, 6.08200669e+09,
           6.08214047e+09, 6.08227425e+09, 6.08240803e+09, 6.08254181e+09,
           6.08267559e+09, 6.08280936e+09, 6.08294314e+09, 6.08307692e+09,
           6.08321070e+09, 6.08334448e+09, 6.08347826e+09, 6.08361204e+09,
           6.08374582e+09, 6.08387960e+09, 6.08401338e+09, 6.08414716e+09,
           6.08428094e+09, 6.08441472e+09, 6.08454849e+09, 6.08468227e+09,
           6.08481605e+09, 6.08494983e+09, 6.08508361e+09, 6.08521739e+09,
           6.08535117e+09, 6.08548495e+09, 6.08561873e+09, 6.08575251e+09,
           6.08588629e+09, 6.08602007e+09, 6.08615385e+09, 6.08628763e+09,
           6.08642140e+09, 6.08655518e+09, 6.08668896e+09, 6.08682274e+09,
           6.08695652e+09, 6.08709030e+09, 6.08722408e+09, 6.08735786e+09,
           6.08749164e+09, 6.08762542e+09, 6.08775920e+09, 6.08789298e+09,
           6.08802676e+09, 6.08816054e+09, 6.08829431e+09, 6.08842809e+09,
           6.08856187e+09, 6.08869565e+09, 6.08882943e+09, 6.08896321e+09,
           6.08909699e+09, 6.08923077e+09, 6.08936455e+09, 6.08949833e+09,
           6.08963211e+09, 6.08976589e+09, 6.08989967e+09, 6.09003344e+09,
           6.09016722e+09, 6.09030100e+09, 6.09043478e+09, 6.09056856e+09,
    ...
           6.10943144e+09, 6.10956522e+09, 6.10969900e+09, 6.10983278e+09,
           6.10996656e+09, 6.11010033e+09, 6.11023411e+09, 6.11036789e+09,
           6.11050167e+09, 6.11063545e+09, 6.11076923e+09, 6.11090301e+09,
           6.11103679e+09, 6.11117057e+09, 6.11130435e+09, 6.11143813e+09,
           6.11157191e+09, 6.11170569e+09, 6.11183946e+09, 6.11197324e+09,
           6.11210702e+09, 6.11224080e+09, 6.11237458e+09, 6.11250836e+09,
           6.11264214e+09, 6.11277592e+09, 6.11290970e+09, 6.11304348e+09,
           6.11317726e+09, 6.11331104e+09, 6.11344482e+09, 6.11357860e+09,
           6.11371237e+09, 6.11384615e+09, 6.11397993e+09, 6.11411371e+09,
           6.11424749e+09, 6.11438127e+09, 6.11451505e+09, 6.11464883e+09,
           6.11478261e+09, 6.11491639e+09, 6.11505017e+09, 6.11518395e+09,
           6.11531773e+09, 6.11545151e+09, 6.11558528e+09, 6.11571906e+09,
           6.11585284e+09, 6.11598662e+09, 6.11612040e+09, 6.11625418e+09,
           6.11638796e+09, 6.11652174e+09, 6.11665552e+09, 6.11678930e+09,
           6.11692308e+09, 6.11705686e+09, 6.11719064e+09, 6.11732441e+09,
           6.11745819e+09, 6.11759197e+09, 6.11772575e+09, 6.11785953e+09,
           6.11799331e+09, 6.11812709e+09, 6.11826087e+09, 6.11839465e+09,
           6.11852843e+09, 6.11866221e+09, 6.11879599e+09, 6.11892977e+09,
           6.11906355e+09, 6.11919732e+09, 6.11933110e+09, 6.11946488e+09,
           6.11959866e+09, 6.11973244e+09, 6.11986622e+09, 6.12000000e+09])

• Data variables: (2)
  • y0
    (dim_0)
    float64
    nan nan nan nan ... nan nan nan nan

    name :

    magn

    long_name :

    Magnitude0

    units :

    V

    batched :

    True

    batch_size :

    1024

    array([nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan])

  • y1
    (dim_0)
    float64
    nan nan nan nan ... nan nan nan nan

    name :

    phase

    long_name :

    Phase0

    units :

    deg

    batched :

    True

    batch_size :

    1024

    array([nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,
           nan])

• Indexes: (0)
• Attributes: (5)

  tuid :

  20240902-014501-705-e85dbc

  name :

  Two-tone

  grid_2d :

  False

  grid_2d_uniformly_spaced :

  False

  1d_2_settables_uniformly_spaced :

  False

qubit.clock_freqs.f01(6.1e9)

rich.print(quantum_device.hardware_config())

{
    'config_type': 'quantify_scheduler.backends.qblox_backend.QbloxHardwareCompilationConfig',
    'hardware_description': {
        'cluster0': {
            'instrument_type': 'Cluster',
            'modules': {
                '6': {'instrument_type': 'QCM_RF'},
                '2': {'instrument_type': 'QCM'},
                '8': {'instrument_type': 'QRM_RF'}
            },
            'sequence_to_file': False,
            'ref': 'internal'
        }
    },
    'hardware_options': {
        'output_att': {'q0:mw-q0.01': 10, 'q1:mw-q1.01': 10, 'q0:res-q0.ro': 60, 'q1:res-q1.ro': 60},
        'mixer_corrections': {
            'q0:mw-q0.01': {
                'auto_lo_cal': 'on_lo_interm_freq_change',
                'auto_sideband_cal': 'on_interm_freq_change',
                'dc_offset_i': None,
                'dc_offset_q': None,
                'amp_ratio': 1.0,
                'phase_error': 0.0
            },
            'q1:mw-q1.01': {
                'auto_lo_cal': 'on_lo_interm_freq_change',
                'auto_sideband_cal': 'on_interm_freq_change',
                'dc_offset_i': None,
                'dc_offset_q': None,
                'amp_ratio': 1.0,
                'phase_error': 0.0
            },
            'q0:res-q0.ro': {
                'auto_lo_cal': 'on_lo_interm_freq_change',
                'auto_sideband_cal': 'on_interm_freq_change',
                'dc_offset_i': None,
                'dc_offset_q': None,
                'amp_ratio': 1.0,
                'phase_error': 0.0
            },
            'q1:res-q1.ro': {
                'auto_lo_cal': 'on_lo_interm_freq_change',
                'auto_sideband_cal': 'on_interm_freq_change',
                'dc_offset_i': None,
                'dc_offset_q': None,
                'amp_ratio': 1.0,
                'phase_error': 0.0
            }
        },
        'modulation_frequencies': {
            'q0:mw-q0.01': {'interm_freq': 80000000.0},
            'q1:mw-q1.01': {'interm_freq': 80000000.0},
            'q0:res-q0.ro': {'lo_freq': 7500000000.0},
            'q1:res-q1.ro': {'lo_freq': 7500000000.0}
        }
    },
    'connectivity': {
        'graph': [
            ['cluster0.module6.complex_output_0', 'q0:mw'],
            ['cluster0.module6.complex_output_1', 'q1:mw'],
            ['cluster0.module2.real_output_0', 'q0:fl'],
            ['cluster0.module2.real_output_1', 'q1:fl'],
            ['cluster0.module8.complex_output_0', 'q0:res'],
            ['cluster0.module8.complex_output_0', 'q1:res']
        ]
    }
}

quantum_device.to_json_file("devices/")

'devices/device_2q_2024-09-02_01-45-02_UTC.json'