Wigner tomography
Wigner tomography is the process of extracting the Wigner function (a quasi-probabilistic distribution) of the quantum state of a cavity. The Wigner function is related to the density matrix of a state, therefore the tomography procedure will allow for the construction of the density matrix and a full reconstruction of the quantum state. Through a standard procedure one can encode the state of a qubit in a superposition of coherent states in a cavity. Therefore, using Wigner tomography it'll be possible to extract the full density matrix of the qubit.
The Wigner function is defined as follows , where is the photon parity operator and is the same parameter in the coherent states and represents a complex vector in the IQ plane.
Using a qubit coupled to the cavity it's straightforward to extract the photon parity of the cavity form a repeated measurement of the qubit through an additional readout resonator. The parity is related to the qubit state as such: m where and are the probabilities of finding the qubit in the excited and ground state respectively, which can be extracted with repeated measurement.
Notice: The example describes the tomography process assuming the cavity was encoded prior.
#
ConfigThe configuration consists of 4 quantum elements:
cavity_I
andcavity_Q
define single input elements and are the I and Q components of the cavity that we'll perform the tomography on.qubit
is the qubit that's coupled to the cavityrr
is the readout resonator that's coupled to the qubit and used to read its state
Each element has its own IF and LO frequencies, and connection ports. Next, for each element we define the relevant operation and pulse:
- For
cavity_I
andcavity_Q
we define thedisplace_pulse
, which will be the real and imaginary parts of the displace pulse. These were separated due to a needed 2d parameter sweep over the amplitudes of the pulses for the tomography. - For the
qubit
we define thex_pi/2_pulse
which is simply a rotation around the x axis - For the
rr
we define thereadout_pulse
- the pulse used for measuring the resonator.
The waveforms used for the displace_pulse
and x_pi/2_pulse
are Gaussians with different parameters.
Generally to displace a cavity one needs to apply a pulse such that it integrates to the desired .
#
ProgramWe first calculate the revival time of the qubit coupled to the cavity. Then, we decide of the range we want to sample for constructing the Wigner function, and the spacing. Once we defined the required parameters, we proceed to the QUA program.
We first define the QUA fixed variable for the amplitude scaling required to shift the cavity by the desired
We than create 2 QUA for_
loops to iterate over the points of the IQ grid. The inner-most for_
loop is for repeated
measurement of the same point in the IQ plane.
Then, in each cycle we perform the tomography procedure:
- We align both cavity components in order to be played simultaneously. We displace the I and Q components by the real and
imaginary parts of , respectively, this is done using realtime amplitude modulation, by multiplying the pulse
with the function
amp(x)
, wherex
is the scaling parameter. - Next, we align the cavity with the qubit to ensure the pulses meant for the qubit start after reaching the desired coherent state for the cavity.
On the qubit we apply a
x_pi/2
operation to bring it to the equator, wait for the revival time, and eventually apply a secondx_pi/2
operation to project the qubit to the excited or ground state. - Finally, we measure the state using the readout resonator and demodulate the reflected signals to get the qubits state on the IQ plane which can then determine its state.
#
Post processingHaving the I,Q results of repeated measurement of the qubit for different we can extract the parity of the cavity at each point by counting the excited and ground state measurements. We can display the results using a heatmap which represents the IQ plane, with the axes being the real and imaginary parts of .