## Abstract

We propose and numerically evaluate a protocol to generate an arbitrary quantum state of the magnetization in a magnet. The protocol involves repeatedly exciting a frequency-tunable superconducting transmon and transferring the excitations to the magnet via a microwave cavity. To avoid decay, the protocol must be much shorter than magnon lifetime. Speeding up the protocol by simply shortening the pulses leads to non-resonant leakage of excitations to higher levels of the transmon accompanied by higher decoherence. We discuss how to correct for such leakages by applying counter pulses to de-excite these higher levels. In our protocol, states with a maximum magnon occupation of up to ∼9 and average magnon number up to ∼4 can be generated with fidelity 0.75.

Export citation and abstract BibTeX RIS

Original content from this work may be used under the terms of the Creative Commons Attribution 4.0 license. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.

Magnets have found commercial applications in magnetic field sensing and information storage, and are a promising building block for long-range information transfer [1, 2] and low-power logic devices [3]. Recently, there has been an interest in bringing such applications into the quantum domain, known as 'quantum magnonics' [4–7], in which the unique properties of the collective excitation of the magnetization, called magnons, are harnessed. Such properties include a high degree of frequency tunability, and low dissipation. The latter is a particular characteristic to the material of choice for quantum magnonics, the ferrimagnetic insulator yttrium iron garnet (YIG) [8], in which evidence for macroscopic (mm-long) coherence lengths were observed [9, 10]. Furthermore, magnons can couple coherently to other systems, such as optical light [11–14], microwaves (MWs) [15–17] and mechanical vibrations [18–20]. In particular, the strong coupling regime between cavity MW photons and magnons, where information can be coherently exchanged between the two excitations at rates 10 MHz much faster than each of their decay rates 1 MHz, can be routinely achieved in experiments. In contrast, the strong coupling regime in the optical domain is notoriously more difficult to achieve, however there are promising developments in this direction [21–26]. Possible applications of such hybrid magnonic quantum systems include magnon-based quantum transducers between light and MWs, as well as transduction involving other degrees of freedom such as phonons [20, 27] and electrons [28–30]. Moreover, platform specific applications such as enhanced magnetic field sensing by using non-classical states of the magnetization can be envisioned [31], potentially reaching the Heisenberg limit [32].

A quantum platform involving magnets requires generation, manipulation, and detection of non-classical states of the magnetization. While single-magnon detection has been demonstrated [33–35] by exploiting the coupling of magnets to superconducting circuits mediated by MW cavities [4, 36–38], experiments so far have probed coherent or thermal magnon states whereas 'true' non-classical states (defined e.g. by a Wigner function with negative regions) remain to be realized. Theoretical proposals in this direction include the generation of cat states of the magnetization using MW cavity photons [39], the generation of single magnon Fock states [40] and cat states [41] via optical means, and the generation of cat states using qubit [42]. In thin films, the generation of entangled pairs of travelling magnons was proposed [43]. All of these proposals are specialized to their respective target states, while there is no known method to deterministically generate an arbitrary magnon state.

In this manuscript we fill in this gap by proposing a protocol to this end, tailored for magnets coupled to a superconducting transmon via a MW cavity, as depicted in figure 1. Our method is based on the protocol proposed in [44], which shows that the deterministic generation of an arbitrary quantum state of a harmonic oscillator can be accomplished by coupling the oscillator to a qubit. A superconducting transmon is a weakly anharmonic oscillator whose first two levels can be treated as a qubit. We show how the protocol of [44] can be modified to take into account the anharmonicity of the superconducting transmon, in order to generate arbitrary quantum magnon states with high fidelity.

The paper is organized as follows. We discuss the system and effective Hamiltonian in section 1. In section 2, we approximate the superconducting transmon by a two-level system and review the protocol from [44] for the cavity mangnonic system. In section 3, we discuss how to correct for errors due to deviations from the two-level approximation, in particular due to interferences from higher energy levels present in the transmon. In section 4, we discuss the fidelity of state generation for paradigmatic examples of non-classical states. We discuss possible improvements to our protocol that can increase the fidelity of generation further in section 5. We review the main findings qualitatively in section 6.

## 1. System

The hybrid system we considered is depicted in figure 1: inside a MW cavity, magnons can couple to a superconducting transmon [4, 36]. The effective coupling Hamiltonian can be obtained by tracing out the MW field. We assume a strong coupling regime, such that the dissipation can be ignored while calculating the effective Hamiltonian. In this section, we present the final results while the detailed calculations are provided in appendix

When the magnet is much smaller than the wavelength of the MWs (∼cm), only the total magnetization **
M
** couples to the MWs, so we need to consider only the uniform magnon mode known as the Kittel mode. It can be quantized via a Holstein–Primakoff transformation [45, 46] that, in the limit of , reduces to:

where is the magnon annihilation operator satisfying the bosonic commutation relations , *M*_{s
} is the saturation magnetization, and the zero-point fluctuations of the magnetization are given by:

with *γ* > 0 being the absolute gyromagnetic ratio and *V*_{m
} being the volume of the magnet. Assuming a spherical magnet and an external magnetic field *H*_{app}, the Hamiltonian is with which can be tuned to be in the range of ∼2*π* **×** 5–10 GHz. For a YIG sample with volume m^{3}, kA m^{−1}, and GHz T^{−1}, we have , implying that the deviations from the ground state are small.

A flux-tunable transmon consists of two Josephson junctions in parallel with Josephson energy respectively, forming a SQUID loop along with a capacitor with capacitance *C* also in parallel [47–50]. When the charging energy is such that , the energy levels of the circuit resemble that of a weakly anharmonic oscillator,

where and with the effective Josephson energy being:

Here Φ is the flux inside the SQUID loop and is the superconducting flux quantum. Such a set of energy levels can be labelled using an annihilation operator and its corresponding number operator . The number of quanta exchanged between the two Josephson junctions is given by:

The Hamiltonian becomes:

The frequency can be tuned dynamically via the external flux in SQUID Φ, to be in a typical range of –20 GHz. Tuning via Φ can be achieved at a fast timescale 10 ns [51, 52], but it comes at the cost of increased decoherence proportional to the allowed tunability , as a higher tunability implies larger susceptibility to flux noises [49, 53]. For our case, we assume that is tunable within a range of GHz.

We assume that the magnet and the transmon are both placed inside a MW cavity which can be modelled as a set of harmonic modes with frequencies and annihilation operator . The details of the modes depend on the shape and size of the cavity. The Hamiltonian for the unloaded cavity is . Its coupling to both the magnet and the transmon corresponds to exchange of quanta and can be modelled as:

where the coupling with the magnetic field at the location of the magnet, and where is the electric field at the location of the transmon. In deriving the Hamiltonian, we have used a rotating wave approximation to drop the terms , and . This is justified since for the frequencies we are considering. If for all *r* modes, and , we can define cavity-dressed magnon and transmon modes with annihilation operators and respectively. Ignoring contributions that are quartic or higher order in and ,

Details about this Hamiltonian are provided in appendix *ω*_{m
} and transmon's frequency *ω*_{q
} stay close to their bare values and respectively. A schematic representation of each term of the Hamiltonian is shown in figure 1.

The function gives the anharmonicity of the transmon and gives the magnon-transmon coupling. Below, we use the same characters and for the values of anharmonicity and coupling. Typically –300 MHz increasing with *n* and the coupling *g* can achieve ∼2*π* **×** 20–30 MHz depending on the system parameters.

Restricting the analysis to only the first two levels of the transmon, we get the Jaynes–Cummings Hamiltonian:

where , are the first two excited states of the transmon, and . As mentioned before, the frequency *ω*_{q
} is tunable. When , the magnons and transmon are decoupled besides a negligible dispersive coupling, while at resonance, , the two subsystems can exchange quanta.

## 2. Ideal state generation

In this section, we review the protocol to generate an arbitrary state of the magnetization via its coupling to a two-level system. This is discussed in [44] theoretically for a general harmonic oscillator. The protocol consists of: (a) exciting the qubit; (b) partially transferring the excitation to the oscillator and (c) repeating the process multiple times to construct the Fock-state superposition that describes the target state. In this setup, the transmon can be excited via a time-dependent drive, while the transmon and the magnons can be put in and out of resonance by tuning the transmon frequency (via an applied flux), effectively changing the coupling rate (see figure 2). The procedure was experimentally implemented to generate quantum states of MW photons coupled to a superconducting transmon [52].

Download figure:

Standard image High-resolution imageIf we start from the qubit's state while the magnon is in the ground state , switching the coupling on for a period and in the absence of dissipation evolves the state as:

where is *n*-magnon state. Thus, we can create an arbitrary superposition of and . If we repeat this process *N* times, we can, in principle, generate an arbitrary magnon state containing up to a maximum of *N* magnons. The protocol then consists of a series of operations in which first the transmon is excited, and then the qubit's frequency is changed to bring it in resonance with the magnons for a time of the order .

To make the process concrete, consider the Hamiltonian restricted to the first two levels of the transmon , where the static part is:

and the dynamic part is:

here and in what follows, the superscript (2) indicates the approximation of the transmon as a two-level system. In this expression, is the externally induced change in the transmon's frequency, and is the transmon drive amplitude. We assume that the qubit and the magnons are far detuned, , implying that magnons and qubit are decoupled at . We review the main steps of the derivation here leaving the detailed calculations to appendices

A non-zero can excite a qubit. The time evolution operator for and reads:

where with *t*_{g
} a pre-determined gate time. This applies a Bloch rotation of magnitude 2*θ* around the axis .

After the qubit is excited, we want to (partially) transfer the quanta to magnons. If the coupling is switched on for a time ( where the detuning is ) and then switched off for a time , (), keeping the qubit driving off , then the time evolution operator can be written as where the phases evolve as:

and the quanta are exchanged via:

with the total number operator .

A given target state containing a maximum of *N*_{m
} magnons, can be systematically achieved by adding magnons one by one as described above. For that, we need a set of time evolution operators and (see figure 2) such that:

and each has a maximum of *j* excitations. Explicitly, for there are a maximum of and *j* magnons in and respectively. The coefficients are found by reversing the problem in order to systematically remove magnons one by one from . Inductively, we assume that is found. Explicitly, we find to ensure achieved via:

Next, we find to ensure giving:

Then, we find and the induction continues until all quanta are removed. The procedure to inductively obtain the protocol's parameters is schematically shown in figure 3.

Download figure:

Standard image High-resolution image## 3. Realistic state generation

The protocol discussed in the previous section is suitable to generate states with only a small number of excitations. States with higher number occupation require longer protocols and hence are spoiled by dissipation. The protocol duration can be decreased via two ways. Firstly, the gate time to excite the qubit can be decreased. Nevertheless, for a transmon qubit, a shorter gate time, where *α* is the level-dependent anharmonicity (see below equation (8)), implies a broader frequency spectrum allowing for significant non-resonant transitions to higher transmon levels. Secondly, a stronger coupling between magnons and the transmon speeds up the excitation transfer step. However, a very strong coupling, where *g* is the level-dependent magnon-transmon coupling (see below equation (8)), implies a leakage of magnon excitations into higher transmon states. Such interferences limit the protocol time and eventually the largest size of the harmonic oscillator state that can be created. In this section, we show that transitions to higher levels can be incorporated by adding extra pulses for cancelling out such transition amplitudes. The extra pulses are analogous to the ones used for correcting qubit gates [54, 55] and can significantly decrease the protocol time allowing for states containing a large magnon number.

For clarity in presentation, we consider the case of and:

while the calculations for the general case are presented in appendices

and the dynamic part is:

For a target magnon state , we want to find a set of operations and satisfying (see figure 2):

where refers to the magnon state and transmon state . As before, starting from , we inductively find where the transfer operator removes a magnon, i.e. ensures:

and the gate operator removes transmon excitations, i.e. ensures:

for .

### 3.1. Transfer of excitation

We discuss how to find satisfying equation (23). Here, we write down the salient features of the derivation in a self-contained manner while the calculations are provided in detail in appendix

We start with the same ansatz as in the simpler two-level case, i.e. for and = 0 for , along with no transmon driving , for which:

This ansatz can be simplified by noticing that the number of excitations are conserved during the evolution, i.e. . An expression for restricted to the finite-dimensional space of can be found by numerically diagonalizing the two Hamiltonians, and , again restricted to this space. Finally, we get a condition on of the form:

where can be found in terms of the eigenvectors and eigenvalues of the two Hamiltonians (see equation (B.8)). This equation can be numerically solved for .

The above ansatz yields a solution in most cases, nevertheless, it can fail for some states, for example, when . Intuitively, when the coupling is on, the excitation from transfers to all but with different oscillation frequencies . Unless these oscillation frequencies are commensurate, the amplitude of the wave-function in will stay non-zero. In cases where this ansatz fails, we can apply a second operation,

where we choose as the values found from a two-level approximation (see equation (18)), and are found by applying the above optimization to . In principle, this process can be continued until a solution is found. However, in our simulations, we never had to go beyond a second step.

### 3.2. Transmon gates

We discuss the main steps of deriving that satisfies equation (24), while the detailed calculations are provided in appendix

As noted before, the application of ensures that has no excitations above the level *j*. We systematically apply gates transferring amplitudes from higher excitation states to lower ones, i.e. (suppressing the index *j*):

where we want to ensure (see condition equation (24)) , and so on. First, we discuss how to implement .

We try the ansatz and where and *φ*_{app} are real, giving:

We can separate the phases by:

where we ignore the interaction between magnons and the transmon, since , and:

For , where *t*_{g
} is the gate time, the solution is the same as that in the two-level approximation, choosing and satisfying:

where . For concreteness, we choose:

The approximate time-dependence is given by:

where .

However, a drive amplitude with a general form can induce transitions to higher transmon levels, e.g. transitions between and . This is problematic as a significant transition amplitude implies excitation number *j* + 2, while we want to have no excitations above *j*. Thus, we consider , where are higher frequency components that cancel out higher order transitions, making .

The error in the transmon gate is given by:

where,

cannot be found analytically. To linear order in the perturbation, we ignore the corrections due to non-commutativity, and impose:

This in turn implies the following set of equations that the drive amplitude must satisfy:

The first condition is a no-bias condition ensuring that . The second condition, with , ensures that the transitions cancel out at the end of the gate. Similarly, the transitions are cancelled by the third condition. In our simulations, we find that for typical parameters, the population beyond level 2 is negligible, so we can ignore the conditions for .

The above conditions leave sufficient freedom in choosing . In what follows, we choose and:

where the parameters *ε*_{µ
} are governed by an underdetermined set of linear equations that can be found from equations (38)–(40). Given , we can find *φ*_{app} by minimizing .

We point out that although we choose a linear approximation, the corrections are not always small. However, in our simulations we find that the error matrix is dominated by extraneous phases for all transmon levels. One might think that these can be removed by adding higher order corrections, however, we observed that such higher order corrections typically shifts the errors from to instead of removing them. Thus, we account for these extra phases by modifying the next steps of the induction.

If is significant (we take the threshold as 0.1 in our simulations), a second gate between and , say , needs to be applied. The correction method is analogous to above. The levels are found to be approximately unoccupied, so we do not need to apply more gates. The time required for is either or depending on whether a second gate is required. If we want to treat the transmon as a two-level system, we would need gate times 7 (found empirically) to avoid leakage, which is even longer than applying two gates as we do here.

## 4. Results

To evaluate the above algorithm, we numerically solve the time evolution that implements the protocol for the density matrix, using QuTiP [56, 57], satisfying:

where the Hamiltonian is given by equations (20) and (21). Assuming zero temperature, we take the Lindblad operators pertaining to magnon decay , transmon decay and transmon dephasing . We did not include magnon dephasing as it has not yet been experimentally observed, although there are theoretical reasons for pure dephasing of magnons [58]. We choose typical parameters GHz, GHz, , , MHz, and coupling MHz. For a given target state , the algorithm described in the previous sections gives the signals and , along with the protocol time . Starting from the ground state, , we numerically solve the above Lindblad equation with the given signals, Δ and , and find the reduced density matrix for the magnon, :

Finally, we calculate the fidelity of state generation which can be interpreted as the probability of finding the magnetization in the target state .

We compare the Wigner function of the target state and the achieved state for some special cases. The Wigner function is defined as:

The Wigner function can be interpreted as follows. Define a symmetric product as an average over all permutations of *p*
's and *q*
's. For example,

In a given state , the average of the above operator is given by [59, 60]:

We can convert this into expectation values of the magnetization operators and via the substitution (note that we are using the dressed operator instead of the bare magnetization operator as the distinction between them is small) getting:

where . The above equation suggests an interpretation of as a quasi-classical probability distribution. However, note that is not always positive and a negative implies a lack of a probabilistic interpretation considered to be a signature of non-classicality. The unit , defined in equation (2), is given by:

where *V*_{m
} is the volume of the magnet. This is a very small angle but it can be measured via the magnet's coupling to the transmon [34, 35].

In figure 4, we consider the case of Fock states and their superposition with the ground state where . The fidelity decreases with increasing *n*. This is caused by two factors. The dominant factor is an increase in the protocol time with the number of magnons leading to higher decay. Secondly, each gate introduces a small amount of error, including residual transfer to higher levels, that compounds with the number of gates, linearly increasing with *n*. It appears that Fock states are more difficult to generate than the superposition states, which could be due to a larger average number of magnons in , compared to , leading to a faster decay. For and , we also show Wigner functions of the target state, and the achieved state . The essential feature of are the azimuthal fringes of order *n*, which is well captured for *n* = 6 and weakly reproduced for *n* = 10. The essential feature of are the *n* circular fringes, that are destroyed by dephasing.

Download figure:

Standard image High-resolution imageA superposition of two semi-classical states, known as a cat state, is useful for error correction protocols because of their insensitivity to particle loss noise [61], particularly as carriers of information [62–64] and in metrology [65–68]. In case of a magnet, we can define semi-classical states as the minimum-uncertainty states with isotropic quantum fluctuations, given by coherent states [69]:

for an arbitrary complex number *ξ*. A coherent state has average and can be interpreted as the quantum state 'closest' to the classical state with magnetization and . An 'odd' cat state given by has two semi-classical components separated by the magnetization ; the same holds for an 'even' cat state . Figure 5 shows the fidelity of generation as a function of *ξ*. The Wigner function for an ideal cat state has two blobs at the expected peaks along with fringes in the centre. The presence of the fringes signify the quantum coherence between the two semi-classical components. In contrast, the Wigner function for a classical distribution, would still have the blobs but no fringes.

Download figure:

Standard image High-resolution imageIn the achieved states, we find such fringes for cat sizes up to 4 until which the fidelities are not too low. At cat size 5, we see that the achieved state are closer to a classical probability distribution (they have 0.85 fidelity with the corresponding classical states ). Differences in the odd and cat states for low sizes can be explained on the basis that odd cats have larger number of magnons, , than even cats, , an effect that is negligible for *ξ* > 2. Fidelity for generation of any other state can be found using the publicly available codes [70].

## 5. Improvements

Here, we focussed on adapting and characterizing the protocol for a system with a magnet and a transmon. There are several methods by which the protocol can be improved, as we now elaborate.

### 5.1. Pre-processing

We considered only a specific type of operations, where the transmon is excited and the excitation is transferred to the magnons. With the help of the MW cavity, one can also apply other operations, in particular, a displacement given by caused by a coherent excitation of the magnet by MWs. This can be used to decrease the time of the protocol as one can choose an appropriate *β*, such that has fewer magnons than . For example, if , then:

implying that has fewer average magnons than . A better *β* can be found numerically that decreases the highest amplitude in to within some given tolerance. Then, we run the protocol to create and apply a MW pulse to displace keeping a high fidelity. Pre-processing with more complicated analog signals can decrease the duration further [42].

### 5.2. High fidelity gates

In our analysis, we find that non-ideal transmon gates are responsible for a significant part of fidelity loss which increases linearly with the number of gates necessary to generate the target state. We took a simple case of correcting at first order perturbation similar to [54, 55]. Simply adding higher order terms in a Dyson or a Magnus expansion does not necessarily improve gate fidelities. When we correct for second order term in Magnus expansion, we find that we can decrease the transitions and significantly, but at the expense of a high transition probability^{
3
}
. This gets corrected only at the third order in Magnus expansion, at which point the transition might play a role. However, there is a rich literature on speeding up gates to ∼2 with low loss of fidelity, such as DRAG [71, 72] or optimal control [73–76], that can not only speed up but also improve the accuracy of gates.

### 5.3. Combining steps

We followed a semi-analytical approach where we systematically add an excitation one by one. In principle, some steps can be combined, e.g. exciting the transmon and transferring the excitation to magnon can occur simultaneously. This becomes infeasible to handle analytically but relevant pulses can, in principle, be found using optimal control theory. A similar case has been discussed in the two-level approximation [77] for a general harmonic oscillator states with 2–3 quanta. Although this approach should result in the least possible errors, it is computationally very expensive and scaling such an approach to high magnon numbers [74, 78] would still require a good initial guess for time-dependent functions describing the transmon driving and the transmon frequency.

## 6. Conclusions

We proposed and numerically demonstrated a protocol to generate an arbitrary quantum state of the magnetization using its coupling to a superconducting transmon, which is a weakly anharmonic oscillator. The protocol involves repeatedly applying transmon gates and transferring an excitation to the magnons. As longer protocol times lead to higher dissipation, we can improve fidelity of generation by reducing the gate time, *t*_{g
}, and the transfer time, *t*_{on}, in the notation of sections 2 and 3. Fast gate pulses with duration ∼2 (*A* is defined in equation (19)) would imply leakage of excitation into higher transmon levels that can be reduced by modifying the input pulses [51, 71, 73, 74]. We use correction pulses by removing errors at first order in a Magnus expansion, similar to [54, 55], bringing the gate time down to . Correcting for second order term in the expansion does not seem to decrease errors, but rather transfer the population to higher levels. The transfer time, *t*_{on}, can be reduced with a higher magnon-transmon coupling. As the magnon-MW coupling can reach very high values, it is likely to be experimentally feasible to have high magnon-transmon couplings. However, a coupling causes leakage of magnons into higher transmon levels leading to loss of fidelities, so we choose . For such couplings, the leakage is small and can be incorporated by renormalizing the duration of the coupling, .

The main source of errors is the magnon dissipation putting a limit on the size of the states. We find that the protocol faithfully produces states up to a largest occupation of ∼9 magnons and an average of ∼4 magnons. We note that our protocol creates states in dressed magnons, given by instead of bare magnons . As discussed in appendix

With a setup to generate an arbitrary magnetization, there still lies a question of how to experimentally measure it, which we leave to future work. The figures can be generated via the publicly available codes [70].

## Acknowledgments

The authors thank H Ribeiro for useful discussions. We acknowledge funding from the Max Planck Society and from the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) through Project-ID 429529648-TRR 306 QuCoLiMa ('Quantum Cooperativity of Light and Matter').

## Data availability statement

The codes to generate the plots in the paper are available at the following URL/DOI: https://github.com/sancharsharma/Magnon-State-Generation.

## Appendix A.: Derivation of effective Hamiltonian

In this appendix, we derive the effective Hamiltonian for a magnet coupled to a transmon indirectly via a MW cavity.

## A.1. Microwaves

The MW cavity is modelled by the classical Hamiltonian density:

It contains several modes labelled by *r* giving the quantization:

and similar holds for **
E
**. Here,

*H*_{r }and are respectively the mode profile and the annihilation operator of mode

*r*. With appropriate normalization, this gives the Hamiltonian,

where is the resonance frequency of mode *r*.

## A.2. Magnet

If the magnet is much smaller than MW wavelengths (∼cm), then only the total magnetization **
M
** couples to the MWs. Under a DC field along with magnitude

*H*

_{app}, its classical Hamiltonian density is:

where **
H
** is the magnetic field due to the cavity. The magnetization is quantized by the Holstein-Primakoff transformation [45, 46], given by a form:

where and *β* are to be found. We know that a spin in an applied field undergoes a Larmor precession with frequency where *γ* is the absolute value of the gyromagnetic ratio. This implies that the Kittel mode should have the resonance frequency, i.e. after quantization:

This implies where (to be found) is the saturation magnetization reduced by zero point fluctuations and is the magnetization reduction because of one (delocalized) spin flip. Using , we find the relations:

This is solved explicitly by:

Typically , so we can resort to the leading order corrections,

For , we can simply substitute:

This gives the magnon Hamiltonian (coupled to the cavity):

where and:

We assumed rotating wave approximation which is valid when .

## A.3. Transmon

We consider a flux-tunable transmon that consists of two Josephson junctions in parallel forming a SQUID loop. Its Hamiltonian is [48, 49]:

Here, is the number of quanta exchanged between the two junctions, is the number of quanta excited by MWs, and are the fluxes through each junction. Under weak coupling, we have . Flux quantization constraints where *m* is an arbitrary integer, Φ is the flux through the SQUID loop, and is the superconducting flux quantum. Defining ,

where the phase is:

and the effective junction energy is:

In the transmon limit , the fluctuations in are small, so we can expand the cosine and get an approximately quadratic Hamiltonian. This suggests the quantization,

that reduces the Hamiltonian to:

ignoring higher order terms. Here and .

## A.4. Dressed magnons and transmons

The dressed states are found by diagonalizing the total Hamiltonian where:

and the interaction:

However, exact diagonalization is not possible analytically, so we resort to a perturbative treatment, in particular a Schrieffer–Wolff (SW) transformation [79]. The SW transformation perturbatively gives the rotation operator such that is diagonal to a given order. We find to linear order in coupling constants *G* that gives the rotated Hamiltonian,

To remove first order corrections, we choose:

where the transmon susceptibility is:

ensuring . Note that we need that implies a limit on the number of excitations s.t. for all *r*,

are satisfied. This gives the rotated Hamiltonian,

As there are four terms in both and , a total of 16 commutators are required. If we ignore all non-resonant terms, such as , and use the Hermiticity of we are left with:

For commuting with all photon operators,

This leads to

where

This procedure yields the rotated Hamiltonian. To return back to the original one, we define the dressed operators: , , and . Up to linear order in *G*, the dressed magnon operator is:

and the dressed transmon operator is:

This gives the final Hamiltonian:

Here, is the effective MW Hamiltonian that is not relevant for our discussion. The renormalized magnon's frequency is Stark shifted by the MWs:

Similarly, the transmon's frequency is:

The Stark shift of higher transmon levels renormalizes the anharmonicity to:

Finally, the coupling of each transmon transition to the magnons is given by:

## Appendix B.: Transfer of quanta

In this appendix, we show how to obtain the periods required for transferring quanta from transmon to magnons. In what follows, denotes a state with *m* magnons and transmon at *q*th level. Given a state having a maximum of *j* excitations, i.e. for and a tunable Hamiltonian:

we want to find s.t. , where the time-evolution operator is:

with Texp being the time-ordered exponential operator. Here, and as derived in the previous section.

First, consider the ansatz with , s.t. the magnons and transmons are resonant for a time-period *t*_{on}, effectively switching on the coupling. Adding an off-period of *t*_{off} afterwards, we get and we want to find to satisfy the above condition. We can derive the time-evolution unitary explicitly using where the total number of excitations are given by:

Thus, we divide the wave-function as:

During the off-time, we have the equations of motion, (using the Schrödinger equation ):

This is a set of linear differential equations which can be solved by diagonalizing the dynamical matrix. For the case *N* = *j*, we can write:

where the vector , and the dynamical matrix is diagonalized by . Similar analysis with works for the on-time,

We want to ensure . We find by putting and then, the condition becomes:

where the vector . This can be numerically solved for .

In the special case of a two-level approximation (ignoring higher levels) and for sufficiently large detuning , the above condition reduces to:

This gives:

In some special cases, equation (B.8) will have no solution. Then, we try the ansatz with a pulse of width , followed by an off-time and a second pulse of an undetermined time *t*_{on},

where . After another off-period of time *t*_{off}, the optimization condition becomes (following the same analysis as before):

which can be solved for . This process can be repeated indefinitely, however, in practice, we did not observe the need to go beyond a second iteration.

## Appendix C.: Transmon gates

We now describe how the transmon gates can be implemented in our scheme. In what follows, denotes a state with *m* magnons and transmon at *q*th level. Given an initial state with a maximum of *j* magnons and excitations, i.e. for and , and a tunable Hamiltonian:

we want to find s.t. , where the time-evolution operator:

with Texp being the time ordered exponential operator and *t*_{g
} being a pre-determined gate time.

Note that we have ignored the coupling between magnons and transmons assuming large enough detuning . Then, the evolution is given by:

Here,

An ideal gate between and , with angular parameters is:

barring an unimportant global phase . The parameters for this gate can be found by ensuring , where giving:

and,

To implement this, we try the ansatz with real to be found. We separate out the phases by the transformation,

where

and satisfies:

where by definition .

First, consider the two-level approximation, where we want to solve:

where is the 'bare' solution to be chosen. The unitary becomes:

where

For such a case, we have:

This resembles the ideal gate except for the extra phases due to the free evolution of the states. For this two-level approximation, we can ensure , where by choosing:

Notice that we can choose any satisfying . For concreteness, we choose a sinusoidal:

## C.1. Including all levels

Transition between and are non-resonant by the frequency *α*_{2}, so the above two-level approximation works when [assuming that the anharmonicity increases with levels, i.e. ]. For a general *t*_{g
}, we modify the pulse and define the error as . The error satisfies:

where the error Hamiltonian is:

This can be divided as:

with the reduced Hamiltonians defined as:

We perform a Magnus expansion to find [54, 55]:

so, to ensure up to first order, we require which translates to:

We used . The last condition ensures that there are no transitions, and we ignored all higher transitions as they are approximately unoccupied throughout the protocol.

## C.2. Satisfying the conditions

The conditions allow for a lot of freedom in choosing . Here, we take an ansatz of a sum of sinusoids, ensuring that the pulse starts and ends at 0,

where the set of indices *ν* can be freely chosen. Equation (C.20) reads:

Equation (C.22) is:

Equation (C.21) is:

To simplify the integration, we define a real integral:

giving

For , this gives:

If we choose the gate time s.t. for some integer *N*, we find that the odd and even frequency components can be separated, and in particular, we need only the odd components to satisfy all the above conditions.

For simulation purposes, we choose the specific case of , the bare transmon, and choose for some natural number *N* giving (for some arbitrary set of indices *H*),

This is a set of linear equations which can be solved to give *ε*_{µ
}. To find an expression for , we use the fact that it satisfies the same recursion relation as that of Bessel functions:

with the initial conditions,

Then, an expression for can be found recursively.

## C.3. Phase

Given , we find *φ*_{app} to minimize . We have where can be decomposed as in equation (C.8)]. Then,

As the leakage into higher levels is suppressed via an appropriate choice of , we expect the 2x2 sub-block of for the subspace spanned by to be approximately unitary. For concreteness, we find the unitary matrix closest to this sub-block, say , in 2-norm via using the singular value decomposition . This gives:

for some constants . Then,

giving the amplitude:

This is minimized at:

## C.4. Gates between higher levels

Typically, there would be leakage from magnons into higher levels of transmon, which needs to be corrected via applying gates between the higher levels. Consider the case when we want to apply a gate between and with angular parameters . We try for a pulse shape:

Again, we have the 'bare' pulse which we take as half-sines. To remove the errors at first order in the Magnus expansion, we find that here satisfies the same set of conditions equations (C.20)–(C.22). Thus, the same calculations as above apply.