Resolving quanta of collective spin excitations in a millimeter-sized ferromagnet

+ See all authors and affiliations

Science Advances  05 Jul 2017:
Vol. 3, no. 7, e1603150
DOI: 10.1126/sciadv.1603150


Combining different physical systems in hybrid quantum circuits opens up novel possibilities for quantum technologies. In quantum magnonics, quanta of collective excitation modes in a ferromagnet, called magnons, interact coherently with qubits to access quantum phenomena of magnonics. We use this architecture to probe the quanta of collective spin excitations in a millimeter-sized ferromagnetic crystal. More specifically, we resolve magnon number states through spectroscopic measurements of a superconducting qubit with the hybrid system in the strong dispersive regime. This enables us to detect a change in the magnetic moment of the ferromagnet equivalent to a single spin flipped among more than 1019 spins. Our demonstration highlights the strength of hybrid quantum systems to provide powerful tools for quantum sensing and quantum information processing.

  • Hybrid Quantum Systems
  • quantum magnonics
  • superconducting qubits
  • ferromagnetic insulators
  • magnons


Engineering interactions between photons and quanta of excitations in atomic and solid-state systems are central to the development of quantum technologies. Cavity and circuit quantum electrodynamics have enabled the realization of many gedanken experiments in quantum optics (13), as well as offering a promising platform for quantum computing (47). Ideas from these fields have been successfully transposed to other architectures such as optomechanical systems (8). In these systems, phonons in mechanical modes can interact with photons of both microwave and optical domains, offering a promising platform for transducing quantum information between microwave-only quantum systems, such as superconducting qubits, and photons in optical fibers (810).

Recently, a novel approach enabling bidirectional conversion between microwave and optical photons using their mutual interaction with macroscopic collective spin excitation modes in ferromagnetic insulators has been explored (1114). Combined with the coherent interaction between the magnons in these magnetostatic modes and superconducting qubits (15, 16), this demonstration opens up the possibility of transducing quantum information between a superconducting quantum processor and photons in optical fibers. However, a key ingredient for using quantum magnonics systems as quantum transducers is the capability of encoding arbitrary qubit states in magnon nonclassical states.


The coherent interaction between magnons and a superconducting qubit was previously demonstrated through the observation of a magnon-vacuum Rabi splitting of the qubit (15). Here, we go further and explore the off-resonant, dispersive regime of quantum magnonics, a promising regime to probe and control magnon states using the superconducting qubit. In circuit quantum electrodynamics, the strong dispersive interaction between a qubit and a microwave cavity has been used as a quantum sensor of the electromagnetic field of the cavity (17, 18) and to create and manipulate nonclassical states of microwave light (19, 20). Here, we demonstrate the ability to reach the strong dispersive regime in quantum magnonics by probing magnons at the level of single quanta in a millimeter-sized ferromagnet using spectroscopic measurements of the qubit. Our demonstration opens up the possibility of measuring, creating, and manipulating macroscopic nonclassical magnon states in ferromagnetic systems.

Our hybrid system consists of a superconducting qubit and a single-crystalline yttrium iron garnet (YIG) sphere inside a three-dimensional microwave cavity (Fig. 1A). The transmon-type superconducting qubit has a resonant frequency of 7.9905 GHz (21). A pair of permanent magnets and a coil are used to apply a magnetic field B0 to the YIG sphere, making it a single-domain ferromagnet. The electric dipole of the qubit and the magnetic dipole of the ferromagnet couple to the electric and magnetic fields of the cavity modes, respectively. The YIG sphere is placed near the antinode of the magnetic field of the transverse electric TE102 cavity mode, so that the cavity field is nearly uniform throughout the 0.5-mm sphere. This makes the uniformly precessing mode, or Kittel mode, the most dominantly coupled magnetostatic mode. As detailed in section S2, the interaction strength gm-c between the TE102 cavity mode, or coupler mode, and the Kittel mode reaches the strong coupling regime (2226), where gm-c/2π = 22.5 MHz is much larger than both the cavity (κc/2π = 2.08 MHz) and magnon (γm/2π = 1.3 MHz) linewidths.

Fig. 1 Hybrid system and qubit-magnon coherent interaction.

(A) Schematic illustration of a ferromagnetic YIG sphere and a superconducting transmon qubit inside a three-dimensional microwave cavity. A magnetic field B0 is applied to the YIG sphere using permanent magnets and a coil. The magnetostatic mode in which spins uniformly precess in the ferromagnetic sphere, or the Kittel mode, couples to the magnetic field of the cavity modes. The qubit and the Kittel mode interact through virtual excitations in the cavity modes at a rate gq-m. (B) The spectrum of the qubit is measured by probing the change of the reflection coefficient Re(Δr) of a microwave excitation resonant, with the probe mode at frequency ωp as a function of the spectroscopy frequency ωs and the coil current I, changing the magnetic field at the ferromagnet. The avoided crossing indicates a coherent interaction between the qubit and the Kittel mode. Vertical dashed lines indicate that I = −4.25 mA, where the qubit and the Kittel mode are hybridized (Fig. 1D), and that I = −5.02 mA, where the qubit-magnon interaction is in the dispersive regime (Figs. 2 to 4). (C) Magnon-vacuum Rabi splitting of the qubit on resonance, with the Kittel mode at I = −4.25 mA. From the fit, we extract the qubit-magnon coupling rate gq-m/2π = 7.79 MHz.

The coupling of the superconducting qubit and the Kittel mode to the same cavity modes creates an effective interaction between these two macroscopic systems (15, 16). To verify that this interaction is coherent, we perform spectroscopic measurements of the qubit with the hybrid system in the quantum regime at a temperature of 10 mK. Although the qubit-magnon coupling is mostly provided by the TE102 cavity mode at 8.4563 GHz (coupler mode), we use the dispersive interaction of the qubit with the TE103 cavity mode at 10.4492 GHz (probe mode) to read out the qubit. This scheme avoids measurement-induced dephasing caused by photon number fluctuations in the coupler mode (16, 17). The change in the reflection coefficient r of a probe microwave tone, resonant with the probe mode and containing less than one photon on average, is measured while exciting the qubit with a spectroscopic microwave tone at frequency ωs. Figure 1B shows the real part of the change of the reflection coefficient, Δr, measured for different currents I in the coil, changing the frequency ωmB0 of the magnons in the Kittel mode. The avoided crossing indicates the coherent interaction between the qubit and the Kittel mode (15). The qubit-magnon coupling strength gq-m of 7.79 MHz, obtained from the magnon-vacuum Rabi splitting of the qubit (Fig. 1C), is much larger than both the power-broadened qubit linewidth γq/2π = 1.74 MHz and the magnon linewidth γm/2π= 1.3 MHz.

We now investigate the dispersive regime of our hybrid system, where the detuning between the bare qubit and Kittel mode frequencies, |ωqbare − ωmbare|, is much larger than gq-m. The exchange of quanta of excitations between the qubit and the Kittel mode, through virtual photons in the coupler cavity mode, is then highly suppressed. The dispersive part of the qubit-magnon Hamiltonian is then given byEmbedded Image(1)where Embedded Image, with |g(e)〉 the ground (excited) state of the transmon qubit, Embedded Image (Embedded Image) is the magnon creation (annihilation) operator, and χq-m is the qubit-magnon dispersive shift (4, 16). This dispersive interaction makes the qubit and magnon frequencies dependent on the state of the other system. More precisely, the qubit frequency ωq(nm) depends on the magnon number state |nm = {0, 1, 2, …}〉, and the magnon frequency ωmi depends on the transmon state |i = {g, e, f, …}〉. As illustrated in Fig. 2A, the strong dispersive regime, where |2χq-m| > max[γqm], enables the observation of magnon number states |nm〉 via magnon number–dependent ac Stark shift of the qubit frequency (17, 18).

Fig. 2 Dispersive qubit-magnon interaction.

(A) Schematic illustration of the hybrid system in the strong dispersive regime. A microwave excitation at frequency ωmw is used to create a magnon coherent state in the Kittel mode. The excitation is detuned from the magnon frequency, with the qubit in the ground state, ωmg, by Δmw = ωmg − ωmw. In the strong dispersive regime, magnon number states |nm〉 (of probability distribution Embedded Image) are mapped into the qubit spectrum as peaks at frequencies Embedded Image, separated by 2χq-m + Δmw and with a spectral weight closely related to Embedded Image. (B) Measurement of the qubit spectrum for a coil current I = −5.02 mA as a function of the Kittel mode excitation frequency ωmw and the spectroscopy frequency ωs. The excitation frequency producing the maximum magnon-induced ac Stark shift of the qubit from ωq (horizontal dashed line) yields an estimation of ωmg/2π ≈ 7.95 GHz (vertical dashed line). The Kittel mode spectrum, measured via its dispersive interaction with the probe mode, appears as a faint vertical line at ~7.95 GHz. The signature corresponding to the two-photon transition involving both the spectroscopy and the excitation photons and exciting both the qubit and a magnon (fig. S1) is indicated by the diagonal dashed line given by Embedded Image, calculated with χq-m/2π = 1.5 MHz at ωmg/2π = 7.95 GHz.

The qubit-magnon dispersive regime is investigated through spectroscopic measurements of the qubit while exciting the Kittel mode at frequency ωmw, detuned by Δmw = ωmg − ωmw from the dressed magnon frequency ωmg. The measurement of the qubit spectrum while sweeping ωmw for a coil current of −5.02 mA and a Kittel mode excitation power Pmw of 7.9 fW is shown in Fig. 2B. Near resonant excitation Δmw ~ 0, the qubit is ac Stark–shifted by the magnon occupancy in the Kittel mode, a signature of the qubit-magnon dispersive interaction similar to the qubit-photon counterpart in circuit quantum electrodynamics experiments (17, 27). The positive magnon-induced ac Stark shift shows that χq-m > 0, and the excitation frequency producing the maximum shift indicates that ωmg/2π ≈ 7.95 GHz. Both these features are consistent with the hybrid system being in the straddling regime (fig. S1) (21, 28). Notably, the signature corresponding to the two-photon transition, from |g, nm = 0〉 to |e, nm = 1〉, involving both the spectroscopy and the excitation photons and exciting both the qubit and a magnon (fig. S1), is also visible at Embedded Image.

We now focus on resolving the magnon number states through measurements with the excitation frequency close to resonance with the Kittel mode (Δmw≪γm). In the qubit spectra shown in Fig. 3, the excitation frequency is fixed at 7.95 GHz, close to resonance with ωmg for I = −5.02 mA (Fig. 2B). The microwave excitation creates a magnon coherent state in the Kittel mode. When coherently driving the Kittel mode, we observe peaks in the qubit spectrum at frequencies higher than the zero-magnon peak. As shown next, these peaks correspond to different numbers of magnons in the Kittel mode.

Fig. 3 Resolving magnon number states.

(A) Measurements of the qubit spectrum at different Kittel mode excitation powers Pmw for a coil current of −5.02 mA and Kittel mode excitation frequency of 7.95 GHz. Black lines show fits of the data to the spectrum of a qubit dispersively coupled to a harmonic oscillator. Color-coded shaded areas show components of the spectrum corresponding to different photon number states |np〉 of the probe cavity mode and magnon number states |nm〉. Vertical offsets are shown by horizontal dashed lines. (B) Measured qubit spectra as a function of Pmw. For clarity, after subtracting a power-dependent offset, Re(Δr) is normalized relative to its maximum for each drive power. For both (A) and (B), vertical dashed lines indicate the frequencies of the qubit |g〉 ↔ |e〉 transitions corresponding to the first four magnon number states, neglecting a power-dependent ac Stark shift, which is small relative to the dispersive shift per magnon for this range of Pmw.

To fit the data of Fig. 3, we used an analytical model of the spectrum of a qubit dispersively coupled to a harmonic oscillator (17). The asymmetric qubit line shape at Pmw = 0 is well reproduced by including in the fit the photonic contribution to the qubit line shape from the dispersive interaction between the qubit and the probe mode (section S4). The fitting parameters for each excitation power are the occupancy of the Kittel mode Embedded Image (where Embedded Image is the projector to the qubit ground state), the qubit-magnon dispersive shift (χq-m), and the excitation detuning (Δmw) (Fig. 3A). More information on the theory and the fitting procedure can be found in sections S3 to S5. We find a detuning Δmw of −0.38 MHz, indicating a bare magnon frequency ωmbare of 7.9515 GHz. The condition for the dispersive regime is therefore respected with a detuning |ωqbare − ωmbare| of 89 MHz, much larger than the qubit-magnon coupling strength. The qubit-magnon dispersive shift χq-m is found to be 1.5 ± 0.1 MHz, in good agreement with the theoretical value of 1.27 MHz (fig. S1). Resolving magnon number states demonstrates that we have reached the strong dispersive regime of quantum magnonics, with the dispersive shift per magnon 2χq-m being larger than both the power-broadened qubit linewidth γq/2π = 0.78 MHz and the magnon linewidth γm/2π = 1.3 MHz.

The average number of magnons Embedded Image in the Kittel mode extracted from the fit of the data is shown in Fig. 4A. At the lowest excitation power of 79 aW, we are able to resolve 0.026 ± 0.012 magnons in the Kittel mode of the ferromagnetic sphere containing about 1.4 × 1018 spins. We are therefore able to resolve a change in the magnetic moment of the ferromagnet equivalent to a single spin flipped among ~5 × 1019 spins. Within our resolution of about 0.01 magnons, no thermal occupancy of the Kittel mode is observed, indicating that the millimeter-sized ferromagnet is indeed in its magnon-vacuum state. Furthermore, the nonlinearity of Embedded Image against the excitation power can be explained by a Kerr nonlinearity in the Kittel mode caused by the anharmonicity of the transmon (3, 29). We find a Kerr coefficient of −0.20Embedded Image MHz in good agreement with the expected value of −0.12 MHz (fig. S1).

Fig. 4 Magnon occupancy and probability distribution.

(A) Magnon occupancy Embedded Image as a function of the excitation power Pmw. Dashed black line and solid orange line show fits to a driven linear and nonlinear Kittel mode, respectively. Inset: Deviations Embedded Image from the linear fit, indicating a significant magnon Kerr nonlinearity. (B) Probability Embedded Image of the first four magnon number states as a function of Pmw. Poisson distributions are shown as solid lines. Inset: Embedded Image for Pmw = 3.1 fW. For both (A) and (B), error bars correspond to 95% confidence intervals.

Finally, we estimate the probability Embedded Image of having nm magnons in the Kittel mode withEmbedded Image(2)where Embedded Image is the qubit spectrum in the analytical model, to which data are fitted, and Embedded Image is its component associated with the magnon number state |nm〉. For 2χq-m≫γm, the probability distribution calculated with Eq. 2 falls back to the Poisson distribution expected for a driven harmonic oscillator (section S5) (17). The probability distributions Embedded Image of the first four magnon number states, shown in Fig. 4B, indicate small deviations from Poisson distributions. This is expected because the qubit-magnon dispersive shift is only slightly larger than the magnon linewidth in our hybrid system (section S5). Nevertheless, our ability to map the probability distribution of magnon number states to the spectrum of a qubit provides a novel tool for investigating quantum states in magnetostatic modes.


Looking forward, the strong dispersive interaction between magnons and a superconducting qubit demonstrated here should enable the encoding of the qubit into a superposition of magnon coherent states in a magnetostatic mode (19, 20). However, to implement this encoding protocol, the qubit-magnon system needs to be deeper into the strong dispersive regime, either by increasing the qubit-magnon coupling strength or by decreasing the magnon linewidth in the quantum regime (23). Together with the recently demonstrated bidirectional conversion between microwave and optical photons in YIG (12), this could pave the way to the transfer of quantum states between superconducting qubits and photons in optical fibers. Combining two very promising candidates for both stationary and flying qubits, such a breakthrough would be an important step toward the realization of a superconducting qubit–based quantum network. Furthermore, the ability to count the number of magnons in a millimeter-sized ferromagnetic insulator in the quantum regime from zero up to a few magnons could be used to study microscopic mechanisms of collective spin excitations, such as decay, scattering, and coupling to a bath of two-level systems. Finally, the demonstrated architecture of quantum magnonics could also be used in applications in spintronics and spin-based quantum information processing (30).


Figure S2 shows the instruments and components used in the experiment. Microwave powers Pr, Ps, and Pmw were calibrated using the input of the cavity as the reference point. At that reference point, the reflection coefficient r is in unity when |ωr − ω10p|≫κ10p, where ωr is the readout frequency and ω10p and κ10p are the resonant frequency and the linewidth of the TE10p cavity mode, respectively, with p = 1,2,3…. By taking into account the attenuation in cables outside and inside the dilution refrigerator, the total attenuations between the microwave sources and the input of the cavity are approximately 81, 122, and 121 dB for the readout, spectroscopy, and Kittel mode microwave excitations, respectively.

The yoke, coil, cavity, and YIG sphere of the hybrid system used in the paper were the same as in the study by Tabuchi et al. (15), whereas the transmon qubit was a different one. The oxygen-free copper microwave cavity had dimensions of 24 × 3 × 53 mm3. An SMA (subminiature version A) connector connected to the cavity was used to measure the reflection coefficient r. A pair of disc-shaped neodymium permanent magnets, with a diameter of 10 mm and a thickness of 1 mm each, were placed at the ends of a magnetic yoke made of pure iron. The magnets produced a static field B0 ≈ 0.29 T in the 4-mm gap between them. The magnetic field can be additionally tuned by a current I in a 104-turn superconducting coil. The field-to-current conversion ratio is approximately 1.7 mT/mA. A YIG sphere glued to an aluminum oxide rod along the 〈110〉 crystal axis was mounted in the cavity at the center of the gap between the magnets. The static field was applied in parallel with the 〈100〉 crystal axis. A transmon-type superconducting qubit, consisting of two large-area aluminum pads and a single Josephson junction (Al/Al2O3/Al), was fabricated on a silicon substrate and was mounted inside the cavity. The qubit and the YIG sphere were separated by 35 mm in the horizontal direction. A double-layer magnetic shield made of aluminum and pure iron covered half of the cavity to protect the qubit from the stray magnetic field of the magnet.


Supplementary material for this article is available at

section S1. Hamiltonian of the hybrid system

section S2. Cavity-magnon coupling

section S3. Qubit spectrum in the dispersive regime

section S4. Qubit spectroscopy—Magnon vacuum state

section S5. Qubit spectroscopy—Magnon coherent state

section S6. Magnon Kerr nonlinearity

fig. S1. Qubit-magnon hybrid system.

fig. S2. Experimental setup.

fig. S3. Cavity-magnon coupling.

fig. S4. Power broadening of the qubit spectrum.

fig. S5. Dispersive qubit-magnon interaction.

fig. S6. Probability distributions.

fig. S7. Magnon Kerr nonlinearity.

fig. S8. Effect of the finite Kerr nonlinearity on the magnon probability distribution.

table S1. Parameters of the hybrid system.

table S2. Comparison between experimental and theoretical values.

table S3. Linewidths of the hybrid system.

table S4. Experimental parameters of the measurements.

References (3134)

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.


Acknowledgments: We acknowledge K. Usami, A. Blais, J. Bourassa, and M. Pioro-Ladrière for useful discussions and P.-M. Billangeon for fabricating the transmon qubit. Funding: This work was partly supported by the Project for Developing Innovation System of MEXT (Ministry of Education, Culture, Sports, Science and Technology), JSPS (Japan Society for the Promotion of Science) KAKENHI (Grants-in-Aid for Scientific Research) (grant nos. 26600071 and 26220601), NICT (National Institute of Information and Communications Technology), JST (Japan Science and Technology Agency) ERATO (Exploratory Research for Advanced Technology) (grant no. JPMJER1601), the Murata Science Foundation, the Research Foundation for Opto-Science and Technology, the Tateisi Science and Technology Foundation, the JSPS Summer Program, and the Natural Sciences and Engineering Research Council of Canada. Author contributions: The measurements were carried out by D.L.-Q. and S.I., with input from all authors. D.L.-Q., Y.T., and Y.N. analyzed the data and wrote the manuscript with input from all authors. Y.T., S.I., and Y.N. conceived the hybrid system with input from A.N., T.I., and R.Y. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.
View Abstract

Navigate This Article