Source code for qiskit_experiments.library.characterization.resonator_spectroscopy

# This code is part of Qiskit.
#
# (C) Copyright IBM 2021.
#
# This code is licensed under the Apache License, Version 2.0. You may
# obtain a copy of this license in the LICENSE.txt file in the root directory
# of this source tree or at http://www.apache.org/licenses/LICENSE-2.0.
#
# Any modifications or derivative works of this code must retain this
# copyright notice, and modified files need to carry a notice indicating
# that they have been altered from the originals.

"""Spectroscopy experiment class for resonators."""

import warnings
from typing import Iterable, Optional, Sequence, Tuple
import numpy as np

from qiskit import QuantumCircuit
from qiskit import pulse
from qiskit.circuit import Parameter
from qiskit.exceptions import QiskitError
from qiskit.providers import Backend

from qiskit_experiments.framework import BackendData, BackendTiming, Options
from qiskit_experiments.library.characterization.spectroscopy import Spectroscopy
from qiskit_experiments.database_service import Resonator
from .analysis.resonator_spectroscopy_analysis import ResonatorSpectroscopyAnalysis


[docs] class ResonatorSpectroscopy(Spectroscopy): """An experiment to perform frequency spectroscopy of the readout resonator. # section: overview This experiment does spectroscopy on the readout resonator. It applies the following circuit .. parsed-literal:: ┌─┐ q: ┤M├ └╥┘ c: 1/═╩═ 0 where a spectroscopy pulse is attached to the measurement instruction. An initial circuit can be added before the measurement by setting the ``initial_circuit`` experiment option. If set, the experiment applies the following circuit: .. parsed-literal:: ┌─────────────────┐┌─┐ q: ┤ initial_circuit ├┤M├ └─────────────────┘└╥┘ c: 1/════════════════════╩═ 0 Side note 1: when doing readout resonator spectroscopy, each measured IQ point has a frequency dependent phase. Close to the resonance, the IQ points start rotating around in the IQ plan. This effect must be accounted for in the data processing to produce a meaningful signal. The default data processing workflow will therefore reduce the two- dimensional IQ data to one-dimensional data using the magnitude of each IQ point. Side node 2: when running readout resonator spectroscopy in a parallel experiment the user will need to specify the memory slot to use. This can easily be done with the code shown below. .. code:: specs = [] for slot, qubit in enumerate(qubits): specs.append(ResonatorSpectroscopy( physical_qubits=[qubit], backend=backend2, memory_slot=slot, )) exp = ParallelExperiment(specs, backend=backend2) # section: warning Some backends may not have the required functionality to properly support resonator spectroscopy experiments. The experiment may not work or the resulting resonance may not properly reflect the properties of the readout resonator. # section: example The resonator spectroscopy experiment can be run by doing: .. code:: python qubit = 1 spec = ResonatorSpectroscopy([qubit], backend) exp_data = spec.run().block_for_results() exp_data.figure(0) This will measure the resonator attached to qubit 1 and report the resonance frequency as well as the kappa, i.e. the line width, of the resonator. # section: analysis_ref :class:`ResonatorSpectroscopyAnalysis` """ @classmethod def _default_experiment_options(cls) -> Options: """Default option values used for the spectroscopy pulse. All units of the resonator spectroscopy experiment are given in seconds. Experiment Options: amp (float): The amplitude of the spectroscopy pulse. Defaults to 1 and must be between 0 and 1. duration (float): The duration in seconds of the spectroscopy pulse. sigma (float): The standard deviation of the spectroscopy pulse in seconds. width (float): The width of the flat-top part of the GaussianSquare pulse in seconds. initial_circuit (QuantumCircuit): A single-qubit initial circuit to run before the measurement/spectroscopy pulse. The circuit must contain only a single qubit and zero classical bits. If None, no circuit is appended before the measurement. Defaults to None. memory_slot (int): The memory slot that the acquire instruction uses in the pulse schedule. The default value is ``0``. This argument allows the experiment to function in a :class:`.ParallelExperiment`. """ options = super()._default_experiment_options() options.amp = 1 options.duration = 480e-9 options.sigma = 60e-9 options.width = 360e-9 options.initial_circuit = None options.memory_slot = 0 return options
[docs] def set_experiment_options(self, **fields): # Check that the initial circuit is for a single qubit only. if "initial_circuit" in fields: initial_circuit = fields["initial_circuit"] if ( initial_circuit is not None and initial_circuit.num_qubits != 1 or initial_circuit.num_clbits != 0 ): raise QiskitError( "Initial circuit must be for exactly one qubit and zero classical bits. Got " f"{initial_circuit.num_qubits} qubits and {initial_circuit.num_clbits} " "classical bits instead." ) return super().set_experiment_options(**fields)
def __init__( self, physical_qubits: Sequence[int], backend: Optional[Backend] = None, frequencies: Optional[Iterable[float]] = None, absolute: bool = True, **experiment_options, ): """Initialize a resonator spectroscopy experiment. A spectroscopy experiment run by setting the frequency of the readout drive. The parameters of the GaussianSquare spectroscopy pulse can be specified at run-time through the experiment options. Args: physical_qubits: List containing the resonator on which to run readout spectroscopy. backend: Optional, the backend to run the experiment on. frequencies: The frequencies to scan in the experiment, in Hz. The default values range from -20 MHz to 20 MHz in 51 steps. If the ``absolute`` variable is set to True then a center frequency obtained from the backend's defaults is added to each value of this range. absolute: Boolean to specify if the frequencies are absolute or relative to the resonator frequency in the backend. The default value is True. experiment_options: Key word arguments used to set the experiment options. Raises: QiskitError: If no frequencies are given and absolute frequencies are desired and no backend is given or the backend does not have default measurement frequencies. """ with warnings.catch_warnings(): warnings.filterwarnings( "ignore", message="deprecation of Qiskit Pulse", module="qiskit_experiments", category=DeprecationWarning, ) analysis = ResonatorSpectroscopyAnalysis() if frequencies is None: frequencies = np.linspace(-20.0e6, 20.0e6, 51) if absolute: frequencies += self._get_backend_meas_freq( BackendData(backend) if backend is not None else None, physical_qubits[0], ) super().__init__( physical_qubits, frequencies, backend, absolute, analysis, **experiment_options ) @staticmethod def _get_backend_meas_freq(backend_data: Optional[BackendData], qubit: int) -> float: """Get backend meas_freq with error checking""" if backend_data is None: raise QiskitError( "Cannot automatically compute absolute frequencies without a backend." ) if len(backend_data.meas_freqs) < qubit + 1: raise QiskitError( "Cannot retrieve default measurement frequencies from backend. " "Please set frequencies explicitly or set `absolute` to `False`." ) return backend_data.meas_freqs[qubit] @property def _backend_center_frequency(self) -> float: """Returns the center frequency of the experiment. Returns: The center frequency of the experiment. Raises: QiskitError: If the experiment does not have a backend set. """ return self._get_backend_meas_freq(self._backend_data, self.physical_qubits[0]) def _template_circuit(self) -> QuantumCircuit: """Return the template quantum circuit.""" circuit = QuantumCircuit(1, 1) if self.experiment_options.initial_circuit is not None: circuit.append(self.experiment_options.initial_circuit, [0]) circuit.measure(0, 0) return circuit def _schedule(self) -> Tuple[pulse.ScheduleBlock, Parameter]: """Create the spectroscopy schedule.""" timing = BackendTiming(self.backend) if timing.dt is None: raise QiskitError(f"{self.__class__.__name__} requires a backend with a dt value.") duration = timing.round_pulse(time=self.experiment_options.duration) sigma = self.experiment_options.sigma / timing.dt width = self.experiment_options.width / timing.dt qubit = self.physical_qubits[0] freq_param = Parameter("frequency") with pulse.build(backend=self.backend, name="spectroscopy") as schedule: pulse.shift_frequency(freq_param, pulse.MeasureChannel(qubit)) pulse.play( pulse.GaussianSquare( duration=duration, amp=self.experiment_options.amp, sigma=sigma, width=width, ), pulse.MeasureChannel(qubit), ) pulse.acquire(duration, qubit, pulse.MemorySlot(self.experiment_options.memory_slot)) return schedule, freq_param def _metadata(self): """Update metadata with the resonator components.""" metadata = super()._metadata() metadata["device_components"] = list(map(Resonator, self.physical_qubits)) return metadata
[docs] def circuits(self): """Create the circuit for the spectroscopy experiment. The circuits are based on a GaussianSquare pulse and a frequency_shift instruction encapsulated in a measurement instruction. Returns: circuits: The circuits that will run the spectroscopy experiment. """ sched, freq_param = self._schedule() circs = [] for freq in self._frequencies: freq_shift = freq - self._backend_center_frequency if self._absolute else freq freq_shift = np.round(freq_shift, decimals=3) sched_ = sched.assign_parameters({freq_param: freq_shift}, inplace=False) circuit = self._template_circuit() circuit.add_calibration("measure", self.physical_qubits, sched_) self._add_metadata(circuit, freq) circs.append(circuit) return circs