![\"image0\"](\"figures/Banner.jpeg\")
\n",
"\n",
""
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "adf7949a",
"metadata": {
"execution": {
"iopub.execute_input": "2025-03-12T18:03:26.665907Z",
"iopub.status.busy": "2025-03-12T18:03:26.664972Z",
"iopub.status.idle": "2025-03-12T18:03:28.469542Z",
"shell.execute_reply": "2025-03-12T18:03:28.468194Z"
},
"tags": [
"imports",
"header_0"
]
},
"outputs": [],
"source": [
"from __future__ import annotations\n",
"\n",
"from pathlib import Path\n",
"\n",
"import numpy as np\n",
"import rich # noqa:F401\n",
"from qcodes.parameters import ManualParameter\n",
"from simulated_data import get_simulated_res_spec_data\n",
"\n",
"import quantify_core.data.handling as dh\n",
"from quantify_core.analysis.base_analysis import Basic2DAnalysis\n",
"from quantify_core.analysis.spectroscopy_analysis import ResonatorSpectroscopyAnalysis\n",
"from quantify_core.data.handling import load_dataset\n",
"from quantify_scheduler import QuantumDevice, Schedule, ScheduleGettable\n",
"from quantify_scheduler.device_under_test.transmon_element import BasicTransmonElement\n",
"from quantify_scheduler.enums import BinMode\n",
"from quantify_scheduler.operations import IdlePulse, Measure\n",
"from quantify_scheduler.operations.acquisition_library import SSBIntegrationComplex\n",
"from quantify_scheduler.operations.pulse_library import (\n",
" SetClockFrequency,\n",
" SquarePulse,\n",
")\n",
"from quantify_scheduler.resources import ClockResource\n",
"\n",
"from utils import display_dict, initialize_hardware, run_schedule, show_connectivity # noqa:F401"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "3a01f27c",
"metadata": {
"execution": {
"iopub.execute_input": "2025-03-12T18:03:28.474776Z",
"iopub.status.busy": "2025-03-12T18:03:28.473250Z",
"iopub.status.idle": "2025-03-12T18:03:28.478940Z",
"shell.execute_reply": "2025-03-12T18:03:28.478030Z"
},
"tags": [
"config",
"header_1"
]
},
"outputs": [],
"source": [
"hw_config_path = \"configs/tuning_transmon_coupled_pair_hardware_config.json\""
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "570b1b7f",
"metadata": {
"execution": {
"iopub.execute_input": "2025-03-12T18:03:28.481033Z",
"iopub.status.busy": "2025-03-12T18:03:28.480871Z",
"iopub.status.idle": "2025-03-12T18:03:28.486164Z",
"shell.execute_reply": "2025-03-12T18:03:28.485223Z"
},
"tags": [
"header_2"
]
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Data will be saved in:\n",
"/root/quantify-data\n"
]
}
],
"source": [
"# Enter your own dataset directory here!\n",
"dh.set_datadir(dh.default_datadir())"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "ef1712d5",
"metadata": {
"execution": {
"iopub.execute_input": "2025-03-12T18:03:28.488581Z",
"iopub.status.busy": "2025-03-12T18:03:28.488063Z",
"iopub.status.idle": "2025-03-12T18:03:29.085420Z",
"shell.execute_reply": "2025-03-12T18:03:29.083904Z"
},
"tags": [
"header_3"
]
},
"outputs": [],
"source": [
"quantum_device = QuantumDevice(\"transmon_device\")\n",
"quantum_device.hardware_config.load_from_json_file(hw_config_path)\n",
"qubit = BasicTransmonElement(\"q0\")\n",
"qubit.measure.acq_channel(0)\n",
"quantum_device.add_element(qubit)\n",
"cluster_ip = None\n",
"meas_ctrl, inst_coord, cluster = initialize_hardware(quantum_device, ip=cluster_ip)"
]
},
{
"cell_type": "markdown",
"id": "bb0ae42b",
"metadata": {},
"source": [
"# Resonator Spectroscopy\n",
"Here we go through resonator discovery and punchout spectroscopy for identifying the resonator and measuring it's resonant frequency."
]
},
{
"cell_type": "markdown",
"id": "1d621a21",
"metadata": {
"tags": [
"no_demo"
]
},
"source": [
"## Setup\n",
"In this section we configure the hardware configuration which specifies the connectivity of our system.\n",
"\n",
"The experiments of this tutorial are meant to be executed with a Qblox Cluster controlling a transmon system.\n",
"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.\n",
"When using a dummy device, the analysis will not work because the experiments will return `np.nan` values.\n",
"\n",
"### Configuration file\n",
"\n",
"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.\n",
"\n",
"The hardware connectivity is as follows, by cluster slot:\n",
"\n",
"- **QCM** (Slot 2)\n",
" - $\\text{O}^{1}$: Flux line for `q0`.\n",
" - $\\text{O}^{2}$: Flux line for `q1`.\n",
"\n",
"- **QCM-RF** (Slot 6)\n",
" - $\\text{O}^{1}$: Drive line for `q0` using fixed 80 MHz IF.\n",
" - $\\text{O}^{2}$: Drive line for `q1` using fixed 80 MHz IF.\n",
"\n",
"- **QRM-RF** (Slot 8)\n",
" - $\\text{O}^{1}$ and $\\text{I}^{1}$: Shared readout line for `q0`/`q1` using a fixed LO set at 7.5 GHz.\n",
"\n",
"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."
]
},
{
"cell_type": "markdown",
"id": "a1e7864e",
"metadata": {
"tags": [
"no_demo"
]
},
"source": [
"### Quantum device settings\n",
"Here we initialize our `QuantumDevice` and our qubit parameters, checkout this [tutorial](https://quantify-os.org/docs/quantify-scheduler/tutorials/Operations%20and%20Qubits.html) for further details.\n",
"\n",
"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."
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "d06a025b",
"metadata": {
"execution": {
"iopub.execute_input": "2025-03-12T18:03:29.089191Z",
"iopub.status.busy": "2025-03-12T18:03:29.088806Z",
"iopub.status.idle": "2025-03-12T18:03:29.094216Z",
"shell.execute_reply": "2025-03-12T18:03:29.093261Z"
}
},
"outputs": [],
"source": [
"def set_readout_attenuation_hardware_config(attenuation_dB: int):\n",
" hwcfg = quantum_device.hardware_config()\n",
" output_att = hwcfg[\"hardware_options\"][\"output_att\"]\n",
" output_att[f\"{qubit.ports.readout()}-{qubit.name}.ro\"] = attenuation_dB\n",
" quantum_device.hardware_config(hwcfg)\n",
"\n",
"\n",
"set_readout_attenuation_hardware_config(32)\n",
"qubit.measure.acq_delay(452e-9)"
]
},
{
"cell_type": "markdown",
"id": "33c13669",
"metadata": {
"title": "BODY"
},
"source": [
"## Resonator Spectroscopy"
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "2276596a",
"metadata": {
"execution": {
"iopub.execute_input": "2025-03-12T18:03:29.096444Z",
"iopub.status.busy": "2025-03-12T18:03:29.096278Z",
"iopub.status.idle": "2025-03-12T18:03:29.103858Z",
"shell.execute_reply": "2025-03-12T18:03:29.102920Z"
},
"tags": [
"hide_input"
]
},
"outputs": [],
"source": [
"# This schedule can also be imported rom quantify_scheduler.schedules\n",
"\n",
"\n",
"def heterodyne_spec_sched_nco(\n",
" pulse_amp: float,\n",
" pulse_duration: float,\n",
" frequencies: np.ndarray,\n",
" acquisition_delay: float,\n",
" integration_time: float,\n",
" port: str,\n",
" clock: str,\n",
" init_duration: float = 10e-6,\n",
" repetitions: int = 1,\n",
" port_out: str | None = None,\n",
") -> Schedule:\n",
" \"\"\"\n",
" Generate a batched schedule for performing fast heterodyne spectroscopy.\n",
"\n",
" Using the :class:`~quantify_scheduler.operations.pulse_library.SetClockFrequency`\n",
" operation for doing an NCO sweep.\n",
"\n",
" Parameters\n",
" ----------\n",
" pulse_amp\n",
" Amplitude of the spectroscopy pulse in Volt.\n",
" pulse_duration\n",
" Duration of the spectroscopy pulse in seconds.\n",
" frequencies\n",
" Sample frequencies for the spectroscopy pulse in Hertz.\n",
" acquisition_delay\n",
" Start of the data acquisition with respect to the start of the spectroscopy\n",
" pulse in seconds.\n",
" integration_time\n",
" Integration time of the data acquisition in seconds.\n",
" port\n",
" Location on the device where the acquisition is performed.\n",
" clock\n",
" Reference clock used to track the spectroscopy frequency.\n",
" init_duration\n",
" The relaxation time or dead time.\n",
" repetitions\n",
" The amount of times the Schedule will be repeated.\n",
" port_out\n",
" Output port on the device where the pulse should be applied. If `None`, then use\n",
" the same as `port`.\n",
"\n",
" \"\"\"\n",
" sched = Schedule(\"Fast heterodyne spectroscopy (NCO sweep)\", repetitions)\n",
" sched.add_resource(ClockResource(name=clock, freq=frequencies.flat[0]))\n",
"\n",
" if port_out is None:\n",
" port_out = port\n",
"\n",
" for acq_idx, freq in enumerate(frequencies):\n",
" sched.add(IdlePulse(duration=init_duration), label=f\"buffer {acq_idx}\")\n",
"\n",
" sched.add(\n",
" SetClockFrequency(clock=clock, clock_freq_new=freq),\n",
" label=f\"set_freq {acq_idx} ({clock} {freq:e} Hz)\",\n",
" )\n",
"\n",
" spec_pulse = sched.add(\n",
" SquarePulse(\n",
" duration=pulse_duration,\n",
" amp=pulse_amp,\n",
" port=port_out,\n",
" clock=clock,\n",
" ),\n",
" label=f\"spec_pulse {acq_idx})\",\n",
" )\n",
"\n",
" sched.add(\n",
" SSBIntegrationComplex(\n",
" duration=integration_time,\n",
" port=port,\n",
" clock=clock,\n",
" acq_index=acq_idx,\n",
" acq_channel=0,\n",
" bin_mode=BinMode.AVERAGE,\n",
" ),\n",
" ref_op=spec_pulse,\n",
" ref_pt=\"start\",\n",
" rel_time=acquisition_delay,\n",
" label=f\"acquisition {acq_idx})\",\n",
" )\n",
"\n",
" return sched"
]
},
{
"cell_type": "code",
"execution_count": 7,
"id": "777242b2",
"metadata": {
"execution": {
"iopub.execute_input": "2025-03-12T18:03:29.106001Z",
"iopub.status.busy": "2025-03-12T18:03:29.105818Z",
"iopub.status.idle": "2025-03-12T18:03:29.112182Z",
"shell.execute_reply": "2025-03-12T18:03:29.111251Z"
}
},
"outputs": [],
"source": [
"freqs = ManualParameter(name=\"freq\", unit=\"Hz\", label=\"Frequency\")\n",
"freqs.batched = True\n",
"freqs.batch_size = 100\n",
"center = 7.364e9\n",
"frequency_setpoints = np.linspace(center - 10e6, center + 10e6, 300)\n",
"\n",
"spec_sched_kwargs = dict(\n",
" pulse_amp=1 / 6,\n",
" pulse_duration=2e-6,\n",
" frequencies=freqs,\n",
" acquisition_delay=196e-9,\n",
" integration_time=2e-6,\n",
" init_duration=10e-6,\n",
" port=qubit.ports.readout(),\n",
" clock=qubit.name + \".ro\",\n",
")\n",
"gettable = ScheduleGettable(\n",
" quantum_device,\n",
" schedule_function=heterodyne_spec_sched_nco,\n",
" schedule_kwargs=spec_sched_kwargs,\n",
" real_imag=False,\n",
" batched=True,\n",
")\n",
"\n",
"\n",
"def simulated_get():\n",
" gettable.old_get()\n",
" return get_simulated_res_spec_data(f=frequency_setpoints, fr=center)\n",
"\n",
"\n",
"if cluster_ip is None:\n",
" gettable.old_get = gettable.get\n",
" gettable.get = simulated_get\n",
"\n",
"meas_ctrl.gettables(gettable)"
]
},
{
"cell_type": "code",
"execution_count": 8,
"id": "87515de7",
"metadata": {
"execution": {
"iopub.execute_input": "2025-03-12T18:03:29.114587Z",
"iopub.status.busy": "2025-03-12T18:03:29.114260Z",
"iopub.status.idle": "2025-03-12T18:03:29.405835Z",
"shell.execute_reply": "2025-03-12T18:03:29.404958Z"
}
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Starting batched measurement...\n",
"Iterative settable(s) [outer loop(s)]:\n",
"\t --- (None) --- \n",
"Batched settable(s):\n",
"\t freq \n",
"Batch size limit: 100\n",
"\n"
]
},
{
"data": {
"application/vnd.jupyter.widget-view+json": {
"model_id": "d5cd863cadba4b818f2149898055e135",
"version_major": 2,
"version_minor": 0
},
"text/plain": [
"Completed: 0%| [ elapsed time: 00:00 | time left: ? ] it"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/html": [
"<xarray.Dataset> Size: 7kB\n",
"Dimensions: (dim_0: 300)\n",
"Coordinates:\n",
" x0 (dim_0) float64 2kB 7.354e+09 7.354e+09 ... 7.374e+09 7.374e+09\n",
"Dimensions without coordinates: dim_0\n",
"Data variables:\n",
" y0 (dim_0) float64 2kB 13.17 13.44 13.25 13.23 ... 14.97 15.03 14.88\n",
" y1 (dim_0) float64 2kB -134.9 -130.7 -124.1 ... 7.038 12.15 17.46\n",
"Attributes:\n",
" tuid: 20250312-180329-116-9bd3d0\n",
" name: resonator spectroscopy\n",
" grid_2d: False\n",
" grid_2d_uniformly_spaced: False\n",
" 1d_2_settables_uniformly_spaced: False<xarray.Dataset> Size: 128kB\n",
"Dimensions: (dim_0: 4000)\n",
"Coordinates:\n",
" x0 (dim_0) float64 32kB 7.357e+09 7.357e+09 ... 7.367e+09 7.367e+09\n",
" x1 (dim_0) float64 32kB 0.01 0.01 0.01 0.01 0.01 ... 0.5 0.5 0.5 0.5\n",
"Dimensions without coordinates: dim_0\n",
"Data variables:\n",
" y0 (dim_0) float64 32kB nan nan nan nan nan ... nan nan nan nan nan\n",
" y1 (dim_0) float64 32kB nan nan nan nan nan ... nan nan nan nan nan\n",
"Attributes:\n",
" tuid: 20250312-180333-013-d332fe\n",
" name: resonator punchout\n",
" grid_2d: True\n",
" grid_2d_uniformly_spaced: True\n",
" 1d_2_settables_uniformly_spaced: False\n",
" xlen: 400\n",
" ylen: 10| setting | \nvalue | \n
|---|---|
| instrument_type | \nCluster | \n
| sequence_to_file | \nFalse | \n
| ref | \ninternal | \n
| setting | \nvalue | \n
|---|---|
| auto_lo_cal | \non_lo_interm_freq_change | \n
| auto_sideband_cal | \non_interm_freq_change | \n
| setting | \nvalue | \n
|---|---|
| interm_freq | \n80000000.0 | \n
| setting | \nvalue | \n
|---|---|
| lo_freq | \n7.500000e+09 | \n
| setting | \nvalue | \n
|---|---|
| interm_freq | \n80000000.0 | \n
| setting | \nvalue | \n
|---|---|
| config_type | \nquantify_scheduler.backends.qblox_backend.Qblo... | \n
| setting | \nvalue | \n
|---|---|
| auto_lo_cal | \non_lo_interm_freq_change | \n
| auto_sideband_cal | \non_interm_freq_change | \n
| setting | \nvalue | \n
|---|---|
| num_bins | \n400 | \n
| index | \n0 | \n
| setting | \nvalue | \n
|---|---|
| offset_ch0_path_I | \nNone | \n
| offset_ch0_path_Q | \nNone | \n
| offset_ch1_path_I | \nNone | \n
| offset_ch1_path_Q | \nNone | \n
| in0_gain | \nNone | \n
| in1_gain | \nNone | \n
| distortion_corrections | \n[{'bt': {'coeffs': None, 'config': 'QbloxFilte... | \n
| out0_lo_freq_cal_type_default | \nLoCalEnum.ON_LO_INTERM_FREQ_CHANGE | \n
| out1_lo_freq_cal_type_default | \nLoCalEnum.OFF | \n
| lo0_freq | \n7500000000.0 | \n
| lo1_freq | \nNone | \n
| out0_att | \n32 | \n
| out1_att | \nNone | \n
| in0_att | \nNone | \n
| setting | \nvalue | \n
|---|---|
| auto_lo_cal | \non_lo_interm_freq_change | \n
| auto_sideband_cal | \non_interm_freq_change | \n
| setting | \nvalue | \n
|---|---|
| auto_lo_cal | \non_lo_interm_freq_change | \n
| auto_sideband_cal | \non_interm_freq_change | \n
| setting | \nvalue | \n
|---|---|
| repetitions | \n80 | \n
| setting | \nvalue | \n
|---|---|
| acq_channel | \n0 | \n
| acq_indices | \n[0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13,... | \n
| thresholded_trigger_count | \nNone | \n
| setting | \nvalue | \n
|---|---|
| acq_protocol | \nSSBIntegrationComplex | \n
| bin_mode | \naverage | \n
| acq_return_type | \n<class 'complex'> | \n
| repetitions | \n80 | \n
| setting | \nvalue | \n
|---|---|
| instrument_type | \nQRM_RF | \n
| setting | \nvalue | \n
|---|---|
| reference_source | \ninternal | \n
| sync_on_external_trigger | \nNone | \n
| setting | \nvalue | \n
|---|---|
| instrument_type | \nQCM_RF | \n
| setting | \nvalue | \n
|---|---|
| sync_en | \nTrue | \n
| channel_name | \ncomplex_output_0 | \n
| channel_name_measure | \n[complex_input_0] | \n
| connected_output_indices | \n(0, 1) | \n
| connected_input_indices | \n(0, 1) | \n
| seq_fn | \nNone | \n
| thresholded_acq_trigger_write_en | \nNone | \n
| thresholded_acq_trigger_write_address | \nNone | \n
| thresholded_acq_trigger_write_invert | \nFalse | \n
| nco_en | \nTrue | \n
| init_offset_awg_path_I | \n0.0 | \n
| init_offset_awg_path_Q | \n0.0 | \n
| init_gain_awg_path_I | \n1.0 | \n
| init_gain_awg_path_Q | \n1.0 | \n
| modulation_freq | \n-135993568.443528 | \n
| mixer_corr_phase_offset_degree | \n0.0 | \n
| mixer_corr_gain_ratio | \n1.0 | \n
| auto_sideband_cal | \nSidebandCalEnum.ON_INTERM_FREQ_CHANGE | \n
| integration_length_acq | \n1000 | \n
| thresholded_acq_threshold | \nNone | \n
| thresholded_acq_rotation | \nNone | \n
| ttl_acq_input_select | \nNone | \n
| ttl_acq_threshold | \nNone | \n
| ttl_acq_auto_bin_incr_en | \nNone | \n
| setting | \nvalue | \n
|---|---|
| graph | \n[[cluster0.module6.complex_output_0, q0:mw], [... | \n
| setting | \nvalue | \n
|---|---|
| instrument_type | \nQCM | \n
| setting | \nvalue | \n
|---|---|
| lo_freq | \n7.500000e+09 | \n
| setting | \nvalue | \n
|---|---|
| q0:mw-q0.01 | \n10 | \n
| q1:mw-q1.01 | \n10 | \n
| q0:res-q0.ro | \n32 | \n
| q1:res-q1.ro | \n60 | \n