Deterministic entanglement of photons in two superconducting microwave resonators

Deterministic entanglement of photons
in two superconducting microwave resonators

H. Wang Department of Physics, University of California, Santa Barbara, CA 93106, USA Department of Physics and Zhejiang California International NanoSystems Institute,
Zhejiang University, Zhejiang 310027, China
   Matteo Mariantoni Department of Physics, University of California, Santa Barbara, CA 93106, USA    Radoslaw C. Bialczak Department of Physics, University of California, Santa Barbara, CA 93106, USA    M. Lenander Department of Physics, University of California, Santa Barbara, CA 93106, USA    Erik Lucero Department of Physics, University of California, Santa Barbara, CA 93106, USA    M. Neeley Department of Physics, University of California, Santa Barbara, CA 93106, USA    A. O’Connell Department of Physics, University of California, Santa Barbara, CA 93106, USA    D. Sank Department of Physics, University of California, Santa Barbara, CA 93106, USA    M. Weides Department of Physics, University of California, Santa Barbara, CA 93106, USA    J. Wenner Department of Physics, University of California, Santa Barbara, CA 93106, USA    T. Yamamoto Department of Physics, University of California, Santa Barbara, CA 93106, USA Green Innovation Research Laboratories, NEC Corporation, Tsukuba, Ibaraki 305-8501, Japan    Y. Yin Department of Physics, University of California, Santa Barbara, CA 93106, USA    J. Zhao Department of Physics, University of California, Santa Barbara, CA 93106, USA    John M. Martinis Department of Physics, University of California, Santa Barbara, CA 93106, USA    A. N. Cleland Department of Physics, University of California, Santa Barbara, CA 93106, USA
September 7, 2019

Quantum entanglement, one of the defining features of quantum mechanics, has been demonstrated in a variety of nonlinear spin-like systems. Quantum entanglement in linear systems has proven significantly more challenging, as the intrinsic energy level degeneracy associated with linearity makes quantum control more difficult. Here we demonstrate the quantum entanglement of photon states in two independent linear microwave resonators, creating -photon NOON states as a benchmark demonstration. We use a superconducting quantum circuit that includes Josephson qubits to control and measure the two resonators, and we completely characterize the entangled states with bipartite Wigner tomography. These results demonstrate a significant advance in the quantum control of linear resonators in superconducting circuits.

Quantum superposition and entanglement have been demonstrated experimentally using spin-like physical systems ranging from atoms to electronic circuits haroche:2006 ; blatt:2008 ; hanson:2007 ; neumann:2008 ; clarke:2008 ; neeley:2010 ; dicarlo:2010 . These systems all display strong nonlinearity, and are used because this nonlinearity allows straightforward quantum control by classical means. The quantum control of linear systems, exemplified by the harmonic oscillator, is by contrast more difficult, and has only been achieved using nonlinear intermediaries: Atoms haroche:2006 ; mabuchi:2002 to control optical cavities, ions to control ion motion leibfried:2003 ; jost:2009 , and superconducting qubits to control photons in microwave resonators schoelkopf:2008 ; hofheinz:2008 ; hofheinz:2009 ; wang:2009b . Quantum entanglement of cavity photons still presents a significant challenge: Experiments have demonstrated maximally-entangled photons in different polarization modes of the same cavity rauschenbeutel:2001 ; papp:2009 , but the entanglement of photons in two physically distinct cavities mariantoni:2008 ; merkel:2010 ; strauch:2010 has proven more elusive.

Figure 1: (Color online) (a) Device circuit schematic. Coupling resonator is connected to and using 20 MHz coupling capacitors, while storage resonator () is connected to () through a 17 MHz coupling capacitor. (b) Qubit spectroscopy, showing probability (color bar) vs. microwave frequency and flux bias for each qubit. Avoided-level crossings near  GHz (dash-dot lines) are due to the coupling resonator and near 6.3 GHz (dashed lines) due to each qubit’s storage resonator. Lower right panel shows magnified view of circled area, upper right panel shows three qubit levels.

Here we show the deterministic generation of entangled photon states in two spatially-separated microwave resonators, achieved by manipulating the photon states with a pair of superconducting phase qubits. We use as a benchmark the generation of NOON states boto:2000 ; nagata:2007 ; shalm:2009 ; jones:2009 ; afek:2010 , comprising a total of photons in the two resonators ( and ), entangled in the quantum state


with photons in resonator and zero in , superposed with the state with the occupation numbers reversed. Such a state has the same degree of entanglement as the Bell state, but with excitations. We also generate MOON states, in which, e.g., resonator has or zero quanta, entangled with resonator with zero or quanta.We fully characterize the two-resonator photon states using bipartite Wigner tomography, which represents a non-trivial extension of single-cavity Wigner tomography haroche:2006 ; leibfried:2003 ; hofheinz:2008 ; hofheinz:2009 ; wang:2009b , and allows us to distinguish entanglement from an incoherent ensemble.

To accomplish this goal, we developed a new quantum circuit comprising two superconducting phase qubits martinis:2009 and three microwave resonators. A sketch of the circuit topology is shown in Figs. 1(a). The circuit includes a coupling resonator , connected to both qubits, and two state storage resonators and , each coupled to one qubit. The resonator frequencies are all different which allows us to frequency-select the qubit-resonator interactions.

The basic method for generating two-resonator entangled states, illustrated in Fig. 2, is to excite and then entangle the two qubits using the coupling resonator. We can swap the resulting Bell state to the two storage resonators, creating an NOON state . If we want to generate higher photon states, we instead selectively excite each qubit to its next higher energy level (Fig. 1(b)), generating the state , thus using the qubits as “qutrits” you:2007 ; neeley:2009 . The required microwave excitation is selective, due to the anharmonicity of the qubits. The qubit excitation is then swapped to each storage resonator through the qubit transition, creating a four-fold entangled state , where the first two letters indicate the qubit states, and the second two numbers the storage resonator states. We then re-excite the qubits to their states, and again swap the excitation to the resonators, generating . This process can be repeated until the entangled state has photons. In the final step, each qubit’s transition is brought on resonance with the corresponding storage resonator, swapping the last excitation and leaving the system in , an -photon NOON state.

Figure 2: (Color online) NOON state preparation sequence. Resonators are represented by dashed lines, and the qubit and transitions by dark and light color solid lines, respectively. (1): is excited to and half-swapped to , generating the Bell state . (2): Coupling resonator swapped to , generating . (3a): NOON state generated by fully swapping each qubit to its storage resonator. For higher states: (4): Qubits excited to . (5): One photon swapped into storage resonators, generating . Steps (4) and (5) are repeated times, generating . (6): Final photon transfer generates (7) (or higher ) NOON state.

We analyze the final resonator state using the qubits as probes. The simplest analysis uses a coincidence measurement: We bring the qubits into resonance with their corresponding storage resonators for an interaction time , after which both qubits are measured simultaneously. The preparation and measurement sequence is repeated times, yielding the joint-qubit state probabilities, , , , and , where is the probability of measuring the first qubit in its ground state with the second qubit in its excited state, and so on. We then vary the interaction time , capturing the evolution of these probabilities. If a resonator has photons, the th photon will swap between the qubit and resonator at a rate scaling as , while for more complicated states, the interaction is a sum of components oscillating at their respective frequencies, weighted by the photon occupation probabilities hofheinz:2009 ; wang:2009b .

Figure 3: (Color online) Qubit coincidence probability measurements for , , and NOON states, and for a mixed state. (blue) and (red) oscillate with interaction time at a rate , while (brown) always remains small. For the mixed state, behavior is nearly identical to the NOON state. Lines are fits to the data. Statistical errors, from the measured probability spread of , are shown only for PRLsupplementary . Horizontal and vertical axes are same for all plots.

For resonators entangled in a NOON state, a joint measurement should correspond to either photons in one resonator and zero in the other, or to the reverse situation; the measurement of the qubits “collapses” the system onto one or the other outcome. Thus in one measurement at most one of the qubits will be in the excited state. When averaged over many measurements, the maximum probability of measuring a particular qubit in is , while the probability of measuring both qubits in should be zero. Therefore we expect that and will oscillate between and , will be zero, and should equal .

Figure 4: (Color online) NOON and mixed state density matrix amplitudes, reconstructed in the photon-number basis from bipartite Wigner tomography; states are labeled where is the photon number in resonator and that in . Bar heights and colors represent matrix element amplitudes. The dominant amplitudes for all three states are in the expected locations, although the off-diagonal elements decrease with , due to the finite qubit dephasing time. Errors for the density matrix elements are not shown but are small PRLsupplementary .

Coincidence measurements are shown in Fig. 3 for NOON states up to , and are consistent with these expectations. The oscillations seen in these measurements are however insufficient proof of resonator entanglement, as an incoherent mixed state can give the same results. To demonstrate this, we have controllably generated a synthetic ensemble comprising a population of states and states PRLsupplementary . Coincidence measurements (Fig. 3) of this synthetic mixed state generate outcomes identical to those of the NOON state.

A more complete resonator measurement, that can resolve entangled from mixed states, uses bipartite Wigner tomography, a nontrivial extension of single resonator tomography leibfried:2003 ; hofheinz:2009 ; wang:2009b . This involves injecting a coherent Gaussian microwave pulse into each of the storage resonators, with controlled amplitude and phase, displacing the resonator states in phase space. The resonators are then simultaneously measured with a joint probability measurement, now as a function of the amplitude and phase of the coherent pulses. From the complete set of measurements, the two-resonator density matrix can be calculated PRLsupplementary .

In Fig. 4 we display the amplitudes of the density matrices measured for resonator NOON states up to , as well as for the mixed state. While there are non-idealities, the desired non-zero matrix elements are clearly apparent for the NOON states, while for the mixed state, the density matrix has only zero-valued off-diagonal elements. The state preparation fidelities, , are found to be (), (), and (). For , the most probable entanglement of formation audenaert:2001 is , while for and , and 0.28, respectively; for the mixed state, the is zero. We also calculate the negativity , where are the eigenvalues of the partial transpose of the density matrix shalm:2009 ; peres:1996 , and indicates entanglement. The negativities are found to be (), (), and (); for the mixed state we find zero with an upper bound of 0.001. The decrease of these values with photon number is compatible with expectations: The state preparation requires phase-coherence of the four-element entangled states for most of the preparation sequence, which is limited by the qubit coherence time PRLsupplementary . Other than this technical limitation, the deterministic generation is completely scalable to large .

A hallmark of NOON states is their rapid phase evolution boto:2000 ; nagata:2007 ; jones:2009 ; afek:2010 , which can be verified by Wigner tomography using two distinct methods. For the state, after entangling the qubits in a Bell state, we wait for a variable time and then swap the state into the storage resonators. The density matrices measured at three different delay times are shown in Fig. 5(a). The phases of the off-diagonal elements rotate with time due to the qubit-resonator frequency difference, as in Fig. 5(b), showing the expected linear dependence.

This phase-measurement method suffers from the short qubit dephasing time. A second method is to change the phase reference for the coherent pulses used in the Wigner tomography, avoiding storage of the state in the qubit. We add an additional phase to the pulses applied to resonator only. The resulting density matrices show the expected rotation of the off-diagonal elements. In Fig. 5(c) we plot the off-diagonal phase angle for different ; the state evolves three times faster than the state, as expected.

Figure 5: (Color online) Phase sensitivity of NOON states up to . (a) Density matrix for NOON state at different times after qubit entanglement. Each element is represented by an arrow, with orientation determined by the phase angle (scale on bottom left). Off-diagonal elements rotate with delay time (black:  ns, blue:  ns, red:  ns), due to frequency difference between the qubit operating point and resonator. Dephasing causes decrease in amplitude of off-diagonal elements. (b) Phase angle of the upper left non-zero off-diagonal element in (a) versus delay time. Line is fit to a rotation rate of  MHz, which corresponds well to set frequency difference. Error bars indicate maximum phase uncertainty. (c) Rotation angle of off-diagonal element versus controlled, additional phase angle used in coherent state displacement. Lines are fits (slopes indicated on plot), consistent with expected phase sensitivity. Confidence bounds are as in (b), with uncertainty increasing with due to increased phase sensitivity.

We also used tomography to measure the NOON state decay PRLsupplementary . We find that the off-diagonal elements decrease at approximately the same rate as the diagonal elements, with a decay time s, consistent with a Markovian environment wang:2009b .

We note that the generation sequence allows an additional flexibility: We can add different numbers of quanta to the resonators, thus generating MOON states, . An example with and is detailed in PRLsupplementary .

The capabilities we have demonstrated here, generating complex entangled photon states in two resonators, hold promise for new quantum architectures in which superconducting resonators play a more central role. The protocol can be extended in a natural way to entangle larger numbers of resonators, allowing, e.g., the direct generation of resonator GHZ and W states leibfried:2003 ; blatt:2008 ; nagata:2007 ; neumann:2008 ; papp:2009 ; jones:2009 . The longer coherence times achievable in superconducting resonators will be of direct utility in performing more complex quantum algorithms, furthering the capabilities of superconducting quantum circuit architectures.

Acknowledgments. We thank K. Audenaert for assistance in entanglement-of-formation calculations. This work was supported by IARPA under ARO award W911NF-04-1-0204. M.M. acknowledges support from an Elings Postdoctoral Fellowship. H.W. acknowledges partial support by the Fundamental Research Funds for the Central Universities in China (Program No. 1A3000*172210301[26]). Devices were made at the UC Santa Barbara Nanofabrication Facility, a part of the NSF-funded National Nanotechnology Infrastructure Network.


  • (1) S. Haroche and J. M. Raimond, Exploring the Quantum - Atoms, Cavities and Photons (Oxford University Press, Oxford UK, 2006).
  • (2) R. Blatt and D. Wineland, Nature (London) 453, 1008 (2008).
  • (3) R. Hanson, L. P. Kouwenhoven, J. R. Petta, S. Tarucha, and L. M. K. Vandersypen, Rev. Mod. Phys. 79, 1217 (2007).
  • (4) P. Neumann et al., Science 320, 1326 (2008).
  • (5) J. Clarke and F. K. Wilhelm, Nature (London) 453, 1031 (2008).
  • (6) M. Neeley et al., Nature (London) 467, 570 (2010).
  • (7) L. DiCarlo et al., Nature (London) 467, 574 (2010).
  • (8) H. Mabuchi and A. C. Doherty, Science 298, 1372 (2002).
  • (9) D. Leibfried, R. Blatt, C. Monroe, and D. J. Wineland, Rev. Mod. Phys. 75, 281 (2003).
  • (10) J. D. Jost et al., Nature (London) 459, 683 (2009).
  • (11) R. J. Schoelkopf and S. M. Girvin, Nature (London) 451, 664 (2008).
  • (12) M. Hofheinz et al., Nature (London) 454, 310 (2008).
  • (13) M. Hofheinz et al., Nature (London) 459, 546 (2009).
  • (14) H. Wang et al., Phys. Rev. Lett. 103, 200404 (2009).
  • (15) A. Rauschenbeutel et al., Phys. Rev. A 64, 050301 (2001).
  • (16) S. B. Papp et al., Science 324, 764 (2009).
  • (17) M. Mariantoni et al., Phys. Rev. B 78, 104508 (2008).
  • (18) S. T. Merkel and F. W. Wilhelm, New J. Phys. 12, 093036 (2010).
  • (19) F. W. Strauch, K. Jacobs, and R. W. Simmonds, Phys.Rev. Lett. 105, 050501 (2010).
  • (20) A. N. Boto et al., Phys. Rev. Lett. 85, 2733 (2000).
  • (21) T. Nagata et al., Science 316, 726 (2007).
  • (22) L. K. Shalm, R. B. A. Adamson, and A. M. Steinberg, Nature (London) 457, 67 (2009).
  • (23) J. A. Jones et al., Science 324, 1166 (2009).
  • (24) I. Afek, O. Ambar, and Y. Silberberg, Science 328, 879 (2010).
  • (25) J. M. Martinis, Quantum Information Processing 8, 81 (2009).
  • (26) J. Q. You, Y.-X. Liu, C. P. Sun, and F. Nori, Phys. Rev. B 75, 104516 (2007).
  • (27) M. Neeley et al., Science 325, 722 (2009).
  • (28) K. Audenaert, F. Verstraete, and B. DeMoor, Phys. Rev. A 64, 052304 (2001).
  • (29) A. Peres, Phys. Rev. Lett. 77, 1413 (1996).
  • (30) See supplementary below for additional discussion and analysis.

Supplementary Material for

“Deterministic entanglement of photons in two superconducting microwave resonators”

I Materials and Methods

The device fabrication is similar to that published previously s_wang:2008 . The half-wavelength superconducting coplanar waveguide resonators are made of rhenium deposited on a -axis single-crystal sapphire substrate, with a m-wide center signal trace and 10 m gaps to the ground plane metallization on either side of the center trace. We place a single lithographed shorting strap connecting the two ground planes at the midpoint of each resonator to improve the quality of the grounding. This point is a voltage node for the fundamental half-wave resonant mode, so that there is minimal additional dielectric loss from the shorting strap’s underlying amorphous Si insulating film s_oconnell:2008 .

In the circuit layout, the coupling resonator is designed to have a higher resonance frequency than the two state storage resonators and . This prevents the qubit frequencies from having to cross the resonator frequency during NOON state amplification. The two storage resonators and are designed with slightly different resonance frequencies, to avoid possible interference between the resonators. The full frequency span in the design was chosen to be about MHz, within the dynamic range of our custom microwave electronics. The two superconducting phase qubits and coplanar waveguide resonators are fabricated together, using our standard multi-layer process s_martinis:2009 . We use interdigitated coupling capacitors between the qubits and the resonators, calculated to each have a coupling capacitance of 1.9 fF. The actual coupling strengths vary slightly with resonator frequency; the detailed component parameters are listed in Table S1.

component (GHz) (GHz) (GHz) (ns) (ns)
- -
- -
- -
- or
Coupling strength
Table S1: Resonator and qubit parameters. The qubit parameters are quoted for when the qubits are off-resonance (see below). The qubit nonlinearity is the frequency difference between the and transitions. The phase coherence time is obtained using a Ramsey interference experiment, which yields the Ramsey time , from which we calculate . measured for resonators similar to those used here s_wang:2009b satisfies ; we assume the same applies here. for the qubits decreases with increasing length of the pulse sequence due to the nature of the phase noise. For most state generation sequences used in this experiment, of typical length to ns, the qubit is in the range of to ns. Coupling strengths correspond to the measured splitting in frequency units, with appearing in the Hamiltonian shown in Eq. (S1).

Ii Generation Sequence Tune-Up

The time required for each qubit-resonator -SWAP is calibrated separately. The swap times obtained from these calibrations scale correctly as with the number of photons in the resonator s_hofheinz:2009 and also depend on the state of the qubit. Examples of the swap calibrations for a one-photon swap between qubit and resonator are shown in Fig. S1, for both the and transitions. The swap time for the transition is approximately times that for the transition. This scaling is as expected, as the multi-level phase qubit can be well-approximated as a weakly nonlinear harmonic oscillator for the energy levels confined by the qubit’s metastable potential well. The scaling confirms that we can use harmonic-oscillator-like raising and lowering operators for the three-level qubit, as was done in the numerical simulations (see below).

Figure S1: (Color online) One-photon swaps between qubit and resonator . (a) Pulse sequence for measuring swaps between the qubit transition and resonator . When off-resonance, the qubit is excited by a pair of pulses, with the first pulse taking the qubit from to and the second pulse making the transition from to . The qubit transition is then tuned close to resonance with resonator and the qubit and resonator left to interact for a time . The triangular pulse at the end is used to measure the probability of the qubit being in the state at the end of the sequence. (b) Left: Qubit swaps with resonator . The qubit is prepared as in (a), but without the second pulse . The qubit probability (color bar) is plotted as a function of interaction time and frequency tuning . Oscillations in probability are due to the qubit excitation swapping with the resonator. Right: Qubit swaps with resonator ; the preparation of the qubit is as in (a) and the plot shows the -state probability as a function of interaction time and qubit flux bias . Red dashed lines indicate the on-resonance values used in the experiment. The swap frequency for the on-resonance transition is times that for the on-resonance transition.

Figure S2 shows the detailed pulse sequence used to generate and measure the NOON state. Sequence steps in general are calibrated and checked separately to maximize preparation fidelity, when possible. For example, we first optimize the qubit Bell state preparation at the end of step I in the preparation sequence detailed in Fig. S2. The fidelity for the Bell state is above , with entanglement of formation 0.59, which agrees well with numerical simulations, performed using a pure dephasing time ns for the qubits (see below). State tomography of the qubits is also as expected: Fig. S2(b) shows the density matrices extracted from coupled qubit tomography, measured at different times during the NOON state preparation. We note that at the end of the sequence, both qubits should return to their ground states. Experimentally we observe small populations in the excited states due to decoherence and pulse imperfections. The exact qubit state after the NOON state generation is measured and is used as the initial state for the qubits when performing Wigner tomography on the storage resonators (see below).

Figure S2: (Color online) NOON state preparation. (a) Sequence to generate and measure NOON state. Qubit transition frequencies (dark-color solid lines) and (light solid lines) are tuned to control coupling of each qubit to resonators (dashed lines); in red, in blue. Vertical lines divide sequence into steps; circled numbers indicate points where qubit tomography was performed. Step I: A pulse excites to , followed by a entangling and resonator . A full -SWAP between and disentangles the resonator and generates the qubit Bell state . For the NOON state this entanglement is transferred to the storage resonators (not shown). For , pulses excite each qubit to yield . Both qubits’ transitions are then tuned to the respective storage resonators, and an -SWAP transfers one photon to the resonators. For , this is repeated to transfer a second photon (not shown). In step III, both qubits’ transitions are tuned to the resonators and an -SWAP transfers the remaining excitation to the resonators, completing the generation, leaving both qubits in their ground states. Step IV: State analysis. The storage resonator () is displaced in phase space by a Gaussian pulse (), followed by an on-resonance interaction with each qubit’s transition for a time , followed by joint-qubit state readout. Many repetitions of this sequence are combined to generate joint-qubit state probabilities , , , and versus . The joint-qubit state probabilities can be used to obtain the two-resonator density matrix (see text). (b) Qubit density matrices extracted from tomography at circled points in (a). (1) Qubit Bell state. (2) Qubit state after transferring first photon into resonators. This is a state with qubits and resonators entangled, so the density matrix is as expected. (3) Qubit state after NOON state generation. Both qubits are very nearly in their ground states, and thus can be used to measure the resonators.

Iii Numerical Simulations

Numerical simulations were performed using the model Hamiltonian


where is the Hamiltonian of the qubit , and ( and ) are the raising and lowering operators for the 3-level qubit (resonator ), is the coupling strength between qubit and resonator , with a sum over all possible qubit-resonator combinations, such that and is the time-dependent, two-tone ( and ) microwave drive on qubit .

The 3-level qubit Hamiltonian was approximated as


where for simplicity we assumed a constant nonlinearity MHz, so that MHz. We approximated the multi-level qubit and by the raising and lower operators for the lowest three levels of a harmonic oscillator, as discussed above.

Decoherence was approximated using the Lindblad master equation taking into account the Markovian environment s_lindblad:1976 , where two characteristic decay times, the energy relaxation time and the pure dephasing time , were used for each resonator and qubit.

The simulations do not directly account for the non-Markovian character of the phase noise in the qubits. To account for this, we used a sequence-time dependent for each qubit, as obtained from Ramsey interference measurements. We used ns for 50 ns-long sequences and ns for 100 ns-long sequences. The resulting simulations agree reasonably well with the experimental measurements.

Iv Bipartite Wigner Tomography

iv.1 Displacement Pulses

The bipartite Wigner tomography is an extension of a method described elsewhere s_leibfried:1996 ; s_hofheinz:2009 ; s_wang:2009b . We displace each resonator with Gaussian pulses and (resonator and , respectively), with variable phase and amplitude. The values of and are distributed over several concentric circles in the complex plane, centered on the origin, where the distribution of values varies approximately with the size of the NOON state. The radii of the circles run through the set , in square-root of photon number units s_hofheinz:2009 . The pulse values are evenly distributed on each circle, with complex values , where ranges from 1 to and is an integer ranging from 1 (), 5 or 6 (), up to to (), depending on the number of photons in the NOON state.

We use every possible combination of values of and distributed over the circles of the same radius for tomography, i.e., for each value of , we use all values with the same amplitude as . The total number of displacement pulse combinations is thus quite large and increases with the photon number in the NOON state, typically involving of order a few hundred pulses. The displacement pulses can be expressed as


iv.2 Photon Populations

For an initial joint-resonator density matrix , the displacement pulses shift the density matrix to


By bringing both qubits (initially in their states) on resonance with the resonators, the joint number states contained in two resonators, i.e., the diagonal elements of , can be read out, as each diagonal element swaps with the qubits at a different rate s_hofheinz:2009 , resulting in a distinct time-dependence for the probabilities , , , and . These can be numerically simulated using the device parameters from Table S1.

As displayed in Fig. S2(b), there is a small non-zero occupation of the excited state of each qubit after the state generation sequence, due to decoherence and pulse imperfections. We use the measured qubit state after the state generation sequence as the qubit initial condition when numerically simulating the tomography. Using these simulations, we obtain the time dependence for the each of the probabilities , , , and corresponding to different combinations of photon number (Fock) states in the two storage resonators. Examples of these probability traces are shown in Fig. S3(a) for some selected initial states. The time-dependent traces for these probabilities, for the set of Fock states , are then used to decompose the experimentally-measured time traces, which yields the probability distribution for the Fock number states contained in the storage resonators. This thus yields the diagonal elements of the experimentally-measured displaced density matrix .

We obtain the diagonal elements of by doing a least-squares fit of the time-dependent probabilities, corrected for measurement fidelity. We use the MATLAB packages YALMIP and SeDuMi for the fitting. The number of fitting parameters is the number of diagonal elements, directly determined by the maximum photon number state contained in the resonators, plus the number of photon quanta added by the displacement pulses and . Fits are done with constraints and to return meaningful probability values. Examples of these fits are shown in Fig. S3, for the NOON state.

Figure S3: (Color online) Joint-state probabilities and photon number distribution analysis for NOON state. (a) Joint-qubit state probabilities from numerical simulations, for a few selected initial resonator states in and , as indicated. Lines are a guide to the eye. Device parameters are as in Table S1. (b) Left panel, measured joint-qubit state probabilities for null displacement pulses . Lines are fits to the data. The fit joint-photon probability distributions are displayed in the right panel. The most significant occupations are in and , as expected for this state. (c) Left panel, measured joint-qubit state probabilities for the displaced state with and (square-root of photon number units). Lines are fits to the data. The fit joint-photon probability distributions , shown in the right panel, are significantly different from those shown in (b) due to the displacement. The matrix elements for the higher photon number states are omitted for viewing, although these may not be negligible for larger displacement pulses.

iv.3 Joint-Resonator Density Matrices

With the diagonal elements of measured for a set of displacements , can be solved for by inverting Eq. (S4) through a linear least-squares fit, while restricting the density matrix to be Hermitian. The resulting density matrices may have small negative eigenvalues due to noise. We use the MATLAB packages YALMIP and SeDuMi to perform semi-definite convex optimization programming, allowing us to find a physical density matrix that is closest to the actual matrix.

When solving for of a NOON state, we restrict the dimension of to an matrix, even though the dimension of can be significantly larger than this due to the displacement pulses. We zero-pad the elements in that have photon indices larger than . This approach is validated by the coincidence measurements (see main text), as we do not detect any frequency components for number states above prior to injecting a displacement pulse.

V NOON State Decay Dynamics

The Wigner tomography allows us to study the decoherence dynamics of the bipartite system s_wang:2009b . The experimental results, compared with numerical simulations, are shown in Fig. S4, with relevant elements in Table S2. We note that the time evolution of the off-diagonal elements in the two-resonator density matrix, which represent inter-resonator coherence, is different from the evolution of the corresponding off-diagonal elements for a single resonator, which represent intra-resonator coherence s_wang:2009b .

Figure S4: (Color online) Decay of the NOON state. (a) Density matrix amplitudes for the NOON state at selected delay times. Bar heights and colors represent the matrix element amplitude, scale on right. The off-diagonal element amplitudes decrease at a rate close to that of the diagonal counterparts. Values for the relevant elements are tabulated in Table S2. (b) Left: Experimental data, showing amplitude of the off-diagonal element as marked by a green ellipse in the first panel in (a), plotted versus post-preparation time. Line is a fit yielding a single decay time , consistent with a Markovian dissipative environment s_walls:1985 . Right: Numerical simulation of the same density matrix element as the experiment, with the line an exponential fit. The simulation uses the measured density matrix as the initial condition and the resonator value listed in Table S1. The simulated decay time is in good agreement with experiment.
16 0.36(1) 0.386(9) 0.328-0.037i 0.009
250 0.31(1) 0.32(1) 0.31+0.01i 0.01
500 0.32(1) 0.31(1) 0.285-0.068i 0.009
750 0.27(1) 0.27(1) 0.259+0.017i 0.008
1000 0.27(1) 0.28(1) 0.261+0.010i 0.008
1250 0.23(1) 0.24(1) 0.215-0.037i 0.008
1500 0.24(1) 0.23(1) 0.229-0.033i 0.008
2000 0.20(1) 0.20(1) 0.198+0.001i 0.006
2500 0.17(1) 0.17(2) 0.155-0.005i 0.008
3000 0.14(2) 0.15(2) 0.13+0.04i 0.01
4000 0.08(2) 0.09(4) 0.06+0.03i 0.01
Table S2: Density matrix elements as a function of delay time for the NOON state shown in Fig. S4. Uncertainties for the diagonal elements are in parentheses, while the magnitude of the uncertainty for the off-diagonal terms are given in the last column.

v.1 Error Analysis

Statistical errors in the qubit probability measurements as well as uncertainty in the amplitude calibration for the displacement pulses and are used to estimate the uncertainty in the amplitude and phase of each density matrix element. These errors are found to be small, in part because the constraints on the analysis filters unrealistic values. Instead we find that slow phase drifts in the electronics, perhaps dominated by ambient temperature fluctuations, give the main phase uncertainties, especially during long measurements. Evaluating a single density matrix usually takes a relatively short time during which these drifts are minimal. However, measuring a series of density matrices such as Fig. 5 in the main paper, takes a much longer time, typically to hours, allowing for more significant drifts. These mostly affect the phases of the density matrix elements, rather than the amplitudes.

v.2 Validation

The bipartite Wigner tomography was validated by several consistency checks. (1) The density matrix is as expected for a range of different states, including the highly entangled NOON states, the energy eigenstates (Fig. S5(a)), the separable (product) state (Fig. S5(b)) and the un-entangled mixed states (see below). (2) The NOON states display the expected phase sensitivity as a function of photon number , as shown in Fig. 5 in the main paper and Table S3. (3) The NOON state fidelity and entanglement of formation agree reasonably well with numerical simulations. (4) The time-dependence of the density matrix elements, showing uniform exponential decay of all elements, is as expected and agrees with numerical simulations, as shown in Fig. S4 and Table S2. (5) The calculated negativities are significantly above zero for the NOON states and precisely zero (within the measurement error) for the unentangled mixed state (see below).

Figure S5: (Color online) Benchmark resonator states used to validate the bipartite Wigner tomography. (a) Photon-number representation of the two-resonator density matrix (error bars not shown), for the energy eigenstate , generated by pumping two photons into resonator while leaving in its ground state. The magnitude of each matrix element is represented by the full length of the corresponding arrow and the phase angle determined by the direction of the arrow in the complex plane (scale on bottom left). Density matrix is as expected from state preparation. (b) Two-resonator density matrix (error bars not shown) for the test resonator state , generated by creating the state in each qubit, followed by a complete -SWAP transfer to each storage resonator. Non-idealities are likely due to inaccuracies in the preparation pulses.
Fig. 5(a) black 0.327 0.460 -0.229+0.239i
blue 0.318 0.464 0.204+0.200i
red 0.301 0.439 0.196-0.156i
Table S3: (Color online) Density matrix elements for the data shown in Figs. 5(a) of the main paper, coded by color (uncertainties not shown). The phase uncertainty in the off-diagonal elements is dominated by the slow phase drift of our electronics (see text).

We note that bipartite Wigner tomography can measure any matrix element with a relatively high accuracy, as we can displace the system by an arbitrary amount in phase space. Even for relatively small off-diagonal elements, such as the desired off-diagonal term in the NOON state, the tomography can unambiguously evaluate this element and measure its sensitivity to external phase perturbations.

v.3 Ensemble of Mixed States

We use a synthetic ensemble of mixed states to illustrate the hazards involved in relying purely on coincidence measurements for demonstrating NOON-state entanglement. The ensemble comprises a mixture of 50% and 50% states, i.e., an ensemble with the same probability of being measured in either of the states forming the NOON state, but without any entanglement. This is done by generating the pure state and measuring the time-dependent joint probabilities , , , and for this state. We then generate the other component of the ensemble, the pure state , and repeat the generation and measurement procedure. Each value of involves 300 repeats of the preparation and measurement sequence for each of the pure states. We then combine the measurement results with equal weights, creating the joint probabilities for the synthetic ensemble; these data are shown in the main paper. The tomographic analysis yielding the density matrices is done in the usual way. The outcome of the ensemble measurements are shown in the main text, with the joint probabilities evolving in a way indistinguishable from the entangled NOON states, but the density matrix for the ensemble revealing a complete lack of entanglement, as witnessed by the negligible values for the off-diagonal elements.

Vi MOON State

The NOON-state generation protocol can be simply generalized to generate MOON states, with different photon numbers in the two entangled resonators. The generation is similar to the NOON state sequence shown in Fig. S2. We assume : After generating the Bell entanglement between two qubits at the end of step I (Fig. S2(a)), we repeat step II times, yielding the four-fold entangled state . The photon amplification and transfer process (step II) is then applied times, but only to qubit and resonator , yielding the state . The final qubit excitations are then transferred in step III, resulting in the MOON state , with the qubits disentangled from the resonators. A MOON state generated in this fashion, with and , is shown in Fig. S6.

Figure S6: (Color online) Measured density matrix (errors not shown) of the MOON state with and . The generation sequence is as described in the text. The state fidelity is , lower than our typical NOON state fidelity, due to technical issues with the particular device used for this experiment; the entanglement of formation is 0.16. The negativity is , indicating a statistically significant entanglement.


  • (S1) H. Wang et al., Phys. Rev. Lett. 101, 240401 (2008).
  • (S2) A. D. O’Connell et al., Appl. Phys. Lett. 92, 112903 (2008).
  • (S3) J. M. Martinis, Quantum Information Processing 8, 81 (2009).
  • (S4) H. Wang et al., Phys. Rev. Lett. 103, 200404 (2009).
  • (S5) M. Hofheinz et al., Nature 459, 546 (2009).
  • (S6) G. Lindblad, Commun. Math. Phys. 48, 119 (1976).
  • (S7) D. Leibfried et al., Phys. Rev. Lett. 77, 4281 (1996).
  • (S8) D. F. Walls and G. J. Milburn, Phys. Rev. A 31, 2403 (1985).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description