Designing Kerr interactions using multiple superconducting qubit types in a single circuit

The engineering of Kerr interactions has great potential for quantum information processing applications in multipartite quantum systems and for investigation of many-body physics in a complex cavity-qubit network. We study how coupling multiple different types of superconducting qubits to the same cavity modes can be used to modify the self- and cross-Kerr effects acting on the cavities and demonstrate that this type of architecture could be of significant benefit for quantum technologies. Using both analytical perturbation theory results and numerical simulations, we first show that coupling two superconducting qubits with opposite anharmonicities to a single cavity enables the effective self-Kerr interaction to be diminished, while retaining the number splitting effect that enables control and measurement of the cavity field. We demonstrate that this reduction of the self-Kerr effect can maintain the fidelity of coherent states and generalised Schr\"{o}dinger cat states for much longer than typical coherence times in realistic devices. Next, we find that the cross-Kerr interaction between two cavities can be modified by coupling them both to the same pair of qubit devices. When one of the qubits is tunable in frequency, the strength of entangling interactions between the cavities can be varied on demand, forming the basis for logic operations on the two modes. Finally, we discuss the feasibility of producing an array of cavities and qubits where intermediary and on-site qubits can tune the strength of self- and cross-Kerr interactions across the whole system. This architecture could provide a way to engineer interesting many-body Hamiltonians and a useful platform for quantum simulation in circuit quantum electrodynamics.

Abstract. The engineering of Kerr interactions has great potential for quantum information processing applications in multipartite quantum systems and for investigation of many-body physics in a complex cavity-qubit network. We study how coupling multiple different types of superconducting qubits to the same cavity modes can be used to modify the self-and cross-Kerr effects acting on the cavities and demonstrate that this type of architecture could be of significant benefit for quantum technologies.
Using both analytical perturbation theory results and numerical simulations, we first show that coupling two superconducting qubits with opposite anharmonicities to a single cavity enables the effective self-Kerr interaction to be diminished, while retaining the number splitting effect that enables control and measurement of the cavity field. We demonstrate that this reduction of the self-Kerr effect can maintain the fidelity of coherent states and generalised Schrödinger cat states for much longer than typical coherence times in realistic devices. Next, we find that the cross-Kerr interaction between two cavities can be modified by coupling them both to the same pair of qubit devices. When one of the qubits is tunable in frequency, the strength of entangling interactions between the cavities can be varied on demand, forming the basis for logic operations on the two modes. Finally, we discuss the feasibility of producing an array of cavities and qubits where intermediary and on-site qubits can tune the strength of self-and cross-Kerr interactions across the whole system. This architecture could provide a way to engineer interesting many-body Hamiltonians and a useful platform for quantum simulation in circuit quantum electrodynamics. arXiv:1709.04026v1 [quant-ph] 12 Sep 2017

Introduction
Experimental progress in coherent superconducting circuits has resulted in a wide variety of qubit designs, from the prototypical flux [1,2,3,4], phase [5,6,7,8] and charge [9,10,11] qubits to more modern designs such as transmon [12] and fluxonium [13] circuits. The mathematical description and quantum dynamics of qubits when coupled transversely to resonators is provided by circuit quantum electrodynamics (cQED). While these devices are all designed to approximate a two-level system in a superconducting circuit, they all possess many additional levels which can affect their interactions with other elements. In particular, when a qubit is coupled to a linear resonator these extra levels can determine the type of nonlinearity that is induced in that resonator. When a superconducting qubit is far detuned from a resonator mode but is nonetheless strongly coupled to it, the second order effect of this induced nonlinearity is known as the Kerr effect. In superconducting devices the Kerr interaction can be significant on the single-photon level [14], a regime which is difficult to reach in optical atom-cavity systems. While the nonlinearity induced by a qubit, via state dependent shifts, enables information to be encoded and controlled in coherent quantum states of a cavity mode, the Kerr effect also distorts these states and rapidly reduces their fidelity over time [15]. The strength of the Kerr effect is therefore an important property of the resonator mode. Here we study how the induced Kerr effects can be engineered by combining different types of superconducting qubits in a single circuit. The existence of mature circuit designs at different experimental groups suggests that implementing this type of design, which two different superconducting qubits on a single substrate, might now be feasible.
While longitudinal coupling between devices is increasingly being studied in superconducting circuits [16,17,18,19], we will consider only consider only the transverse coupling described by cQED. In this configuration, coupled qubits can give rise to Kerr interactions both on single cavities and between pairs of cavities. The self-Kerr is a nonlinear shifting of a resonator frequency as a function of the number of photons in the mode. A simple quantum system where this can been seen is the quantum Duffing oscillator, with its term proportional to (a † a) 2 in the Hamiltonian, where a is the photon annihilation operator for a resonator mode. In the classical limit, this becomes the quadratic dependence of the refractive index on the electric field strength, sometimes known as self-phase modulation [20,21]. This effect manifests itself at the second order in a series expansion of the Jaynes-Cummings interaction in the dispersive limit [22]. With two resonators, this idea can be extended to the cross-Kerr effect, known as cross-phase modulation in optics, which is indicated by a term proportional to a † ab † b. Cross-Kerr can be realised by coupling two cavities of different frequency via a single qubit device. Both of these effect are also of interest in a variety of analogous quantum systems, such as cavity optomechanics [23,24], where large nonlinearities can feasibly be produced.
In optical systems Kerr effects are typically very weak due to the weak coupling that can be realised between natural atoms and optical modes. By contrast, in cQED the strong couplings in the dispersive regime mean that very strong self-and cross-Kerr effects can be produced, and both nonlinearities have been demonstrated on the level of single photons in recent experiments [25,14,26]. It is known that this large interaction strength [27] is required to implement logic qubit gates using the cross-Kerr interaction [28,29], while strong self-Kerr can be used to generate Schrödinger cat states in cavity [30]. Recent developments in three-dimensional microwave cavities have produced coherence times of nearly 1ms [31], leading to great interest in using the cavity field to store continuous-variable quantum information in Schrödinger cat states [32]. Specific protocols to prepare [33], encode quantum information in [15], stabilise [34], entangle [18] and error correct [35] such cavity states have been demonstrated experimentally. Recent work also discusses how to transmit generalised cat states out of a cavity into a transmission line [36]. In many cases, for example if trying to build a long-lived quantum memory, it is desirable to have very weak nonlinearites so that stored states can be preserved with high fidelity. A method to design different combinations of Kerr nonlinearities, and even vary them dynamically could therefore be extremely useful.
In this paper we demonstrate that, by combining different types of superconducting qubits with different anharmonicities in a single circuit, the strength of both the self-and cross-Kerr interactions can be designed. In addition, the use of tunable qubits allows these interactions to be controlled dynamically to a much larger extent than is possible in circuits based on one type of qubit. We focus on dynamics of quantum states in the single photon regime and well within the coherent timescales where quantum information processing is designed to operate. In Section 2 we start by considering modifying self-Kerr by coupling two different superconducting qubits to a single mode. As coherence times improve, self-Kerr becomes more significant in relation to the cavity decay rate and introduces larger phase distortions to these cavity states, reducing the fidelity of any stored quantum state. In Section 2.1 a time-independent perturbation theory is used to show that the self-Kerr interaction can be passively eliminated from the cavity field while still retaining the number splitting required for controlling and measuring the cavity. In Section 2.2 we show that this result persists when a more realistic model of transmon qubits, the quantum Duffing oscillator with several levels, is considered. Using exact numerical simulation of the Schrödinger equation, we show that the perturbation theory method well describes the behaviour of the full quantum system. In Section 2.3 we demonstrate numerically that this self-Kerr cancellation is also possible in the case where the qubit, in this case fluxonium [13] cannot be modelled as a Duffing oscillator. These results suggest that such a setup could be used to increase the fidelity of quantum memories. A similar cancellation idea has recently been demonstrated using a single multi-level system in cQED [37].
In Section 3, we extend the principle to modifying the cross-Kerr interaction between two cavities. The cross-Kerr effect is of particular interest in quantum optics because of its ability to entangle two optical degrees of freedom and therefore has the potential to be used for quantum information processing. The interaction can be used to perform entangling gates between qubits encoded in superpositions of coherent states or Fock states depending on the computational basis that is being used. For travelling photonic states at optical frequencies it is difficult to achieve sufficient strong nonlinearities to perform the entanglement of any reasonable timescale. However this can be significantly improved in trapped photonic states in microwave circuit QED. If two modes, both starting in the state a i |0 + b i |1 , interact by the cross-Kerr interaction for a time t, the resulting state is (1) The interaction only acts on the final term of the superposition as it has no effect if there is a vacuum state in either mode, and waiting a time t = π/g will therefore change the sign of this term. This realises a controlled-phase gate (e iπ a † a b † b ) to obtain a maximally entangled state between two modes. A similar scheme has recently been theoretically proposed for entangling propagating photons [38]. We consider two cavities which are coupled together by two intermediary qubits with opposite anharmonicities. In this configuration we can realise two different regimes. In one regime the cross-Kerr interaction is entirely cancelled by the presence of a second superconducting qubit providing excellent mutual isolation of the cavities (Section 3.1). In the other regime one of the superconducting qubits is moderately detuned and the cross-Kerr produces a maximally entangled state between the cavities (Section 3.2). Switching between these two configurations through the use of a tunable qubit allows an entangling gate to be performed and then the interaction switched off to maintain this entangled state, forming the basis of a quantum logic gate.
Finally, in Section 4 we consider how this Kerr engineering could scale up to manybody systems. Circuit QED has opened avenues in many-body physics by providing a potential way to engineer an array of cavities and using superconducting qubits to modify the on-site and coupling parameters. SQUID loops have also been studied as a way to induce correlations between quantum resonators through the use of on-site and intermediate devices [39]. This kind of setup is of interest for studying phenomena such as quantum phase transitions [40,41], quench dynamics [42] and for implementation of a quantum simulator [43] in a driven-dissipative system. It is theoretically possible to realise an analogue of the Bose-Hubbard model [44,45] in circuit QED using coupling to qubits to provide the on-site nonlinearity, known as the Jaynes-Cummings lattice [46,47,48]. In other work, the phase diagram of a line of cavities possessing only cross-Kerr interactions has also been solved theoretically [49], and a very recent experiment reports a quantum phase transition in a line of 72 superconducting cavities [50]. We here consider a line of cavities, each with an on-site qubit and coupled by intermediary qubits along the line and show that, in principle, the use of qubits with different anharmonicities in this setup could allow the realisation of novel many-body Hamiltonians. The tunable qubits would also enable the study of quantum quenches and other interesting phenomena [51], where parameters of the system Hamiltonian are changed suddenly, revealing information about many-body systems.

Modifying cavity self-Kerr for quantum memory
The self-Kerr induced by a qubit on a linear resonator can be derived from diagonalising the Jaynes-Cummings (JC) Hamiltonian where a is the cavity annihilation operator, σ ± are the qubit raising and lowering operators, ω c is the cavity frequency, ω q is the qubit frequency and g is the cavityqubit coupling. The JC model assumes that the rotating wave approximation (RWA) [52] holds and that interactions that create excitations simultaneously in the cavity and transmon are therefore very unlikely. This means that the model is only valid when |ω c − ω q | |ω c + ω q |. Outside of this parameter range, the Rabi model is required [53], giving rise to interesting new physics [54]. We discuss the breakdown of the RWA in the context of our results in Appendix Appendix A.4. The eigenvalues of the JC model can be solved for exactly by diagonalising the Hamiltonian in the |n, ↑ , |n + 1, ↓ basis, with adjacent blocks uncoupled because the total number of excitations conserved by the interaction. The Hamiltonian can be written in 2 × 2 matrix blocks where n+1 is the total number of excitations in the system (0 ≤ n ≤ N ). Separately, we also have the ground state |0, ↓ and the state |N, ↑ which has the maximum number of excitations permitted in the chosen basis. As the blocks are not coupled, we can diagonalise them separately, giving the eigenvalues where ∆ = ω c − ω q is the detuning between the cavity and qubit. The eigenstates of the system are |n, + = cos θ n 2 |n + 1, ↓ + sin θ n 2 |n, ↑ , |n, − = − sin θ n 2 |n + 1, ↓ + cos θ n 2 |n, ↑ , with θ n = arctan(2g √ n + 1/∆). In the dispersive regime ∆ g, we can expand the square root term of Eqn. 4 in the small parameter g/∆ and assume that the qubit always remains in its ground state. This gives use the approximate eigenvalues Designing Kerr interactions using multiple superconducting qubit types in a single circuit6 The nonlinear n 2 term describes the cavity self-Kerr, the sign of which is completely determined by the sign of the detuning. This dispersive energy spectrum can also found by approximately diagonalising the JC Hamiltonian by application of an appropriate unitary transformation [22,55] or by calculating a perturbation expansion up the fourth order in the interaction [56].

Eliminating cavity self-Kerr with a pair of two-level qubits
A single qubit induces a self-Kerr interaction proportional to g 4 /∆ 3 on the cavity, which can be modified by changing either the coupling or the detuning between the qubit and the cavity. Detuning the qubit very far from the cavity, however, turns off the number splitting required for control and read out of the cavity field. We therefore add a second qubit and investigate whether this can be used instead to modify the nonlinearity of the cavity. The Hamiltonian for the one-cavity and two-qubit (1C2Q) system is given by where g 1 and g 2 are the couplings to the two qubits, with their respective detunings ∆ 1 = ω c − ω q1 and ∆ 2 = ω c − ω q2 . If ∆ 1 = ∆ 2 this is known as the two-atom Dicke model [57]. We can also diagonalise this system block-wise in 4 × 4 blocks given by for ∆ + = ∆ 1 + ∆ 2 and ∆ − = ∆ 1 − ∆ 2 . The analytical form of the eigenvalues of this block is very complicated in general, but if we consider the case g 1 = g 2 = g, and ∆ 1 = −∆ 2 = ∆, then the eigenvalues are (n + 1)ω c |n, ↓↓ or |n, ↑↑ , (n + 1)ω c ± 4g 2 (n + 3 2 ) + ∆ 2 |n, ↑↓ or |n, ↓↑ .
As for the single-qubit case in Equations (5) and (6), the eigenvectors in the dispersive regime are close to being the bare eigenstates. The first eigenvalue applies to the states where the qubits are both in the same state. This means that, in the dispersive regime, if both qubits remain in the ground state the cavity self-Kerr is completely eliminated. The cavity state feels no nonlinearity, but the number splitting of the qubit transition frequency remains, which allows the qubit to still be used to control and measure the cavity. This result is applicable to the flux qubits used in a recent experiment, which have strong coupling and very large anharmonicities and can therefore be considered as close to a true two-level system [58].

Eliminating self-Kerr with two Duffing-like qubits
In practice, many superconducting qubits are not ideal two-level systems, but are instead relatively weakly anharmonic oscillators. The transmon [12], for example, can be modelled as a quantum Duffing oscillator [22] with a constant anharmonicity when only the lowest few levels are relevant. With many qubit levels the exact diagonalisation used above becomes increasingly difficult, but we can still use perturbation theory to derive expressions for the eigenenergies in the dispersive limit. The Hamiltonian for a cavity coupled to a single transmon is where the transmon field has annihilation operator b and nonlinearity χ. To calculate perturbation theory results that include the self-Kerr, we need to calculate terms up to the fourth order in the interaction Hamiltonian. The presence of the third transmon level adds a correction to the energies associated with two excitations moving from the cavity into the qubit (full details given in Appendix A), giving the modified cavity eigenenergies to fourth order in g/∆. We can see that, with a finitely anharmonic qubit, the sign of the nonlinearity induced on the cavity is not fully determined by ∆, but also depends on the sign of χ [59]. Two qubit levels are therefore only sufficient to predict even the sign of the interaction when |χ| |∆|, which is not generally satisfied in the dispersive regime of cQED. While χ is always negative in transmon devices, flux qubits can behave like a Duffing oscillator in the absence of external flux and produce energy levels with a positive anharmonicity of a few hundred megahertz [3].
Next, we consider coupling a second Duffing-like qubit to the same cavity. The total Hamiltonian is equivalent to and the eigenenergies are To cancel the cavity self-Kerr with the additional qubit levels included, we now need to choose ∆ 1 = −∆ 2 and χ 1 = −χ 2 . We therefore need to use two different types of qubit in a single circuit with opposite nonlinearities, for example transmons and the flux qubits discussed above. We must also ensure that when designing the circuit we choose the parameters such that g i |ω c − ω i | |ω c + ω i | so that we remain within the dispersive regime and RWA.
We can numerically solve the Schrödinger equation to see whether this result holds when all orders of H 1C2T are included. The action of the Kerr effect on a coherent state is to distort the state, transiently producing various coherent state superpositions before completing a full revival of the initial state [14]. This gives rise to a periodic collapse and revival of the cavity amplitude | a | when starting with a coherent initial condition. In Figure 2 we show that our two-qubit setup appears to cancel not just the cavity self-Kerr, but all orders of the cavity nonlinearity when the detunings and anharmonicities are equal and opposite. This occurs in this idealised case because the energy spectrum is completely symmetrical in the frame of the cavity. In practise, the extent to which the interaction can be reduced may be limited by small non-RWA effects, as discussed in Appendix A.4. We neglect decoherence processes as self-Kerr only has significant effects when it is many times greater than the dissipation and the effects we are looking at therefore occur on much shorter timescales than the dissipation. We see that a single transmon causes the periodic revivals expected, while after the addition of the second qubit an initial coherent state undergoes almost trivial evolution. Some high frequency oscillations appear, caused by the fact that the bare cavity eigenstates are not eigenvalues of the full coupled system. This demonstrates that different types of qubit used together in the same circuit can be used to achieve modifications of the Kerr effect that are not possible using transmons only, in this case with potential to improve the quality of quantum memory.
MHz, χ 2 = 300 MHz, and g 1 = g 2 = 100 MHz. With a single transmon, the amplitude oscillates due to the cavity self-Kerr. This is completely eliminated by the addition of the second device and the amplitude remains constant, apparently indicating that all orders of the anharmonicity are cancelled. The small reduction in amplitude and the faster oscillations are due to the fact that the eigenstates of the total system are not pure cavity eigenstates.

Improved quantum memory using fluxonium qubits
To demonstrate that this principle applies to other qubit devices, even those that cannot be modelled as a Duffing oscillator, we also show simulations using a fluxonium qubit [13] to cancel the cavity self-Kerr, as was proposed in an earlier paper [60]. The transmon device provides a negative nonlinearity and can be used to control the cavity, while we use the fluxonium to provide an opposite, positive nonlinearity. While our perturbation theory results do not apply fully to this system, as fluxonium has nonzero matrix element for all qubit transitions, rather than just adjacent levels, we see that significant cancellation can still be achieved.
The Hamiltonian for the cavity-transmon-fluxonium system is where ω c is the cavity frequency, ω t is the transmon fundamental frequency, χ is the transmon nonlinearity, ω f,j is the energy of the j-th level of the fluxonium, λ t,jk and λ f,jk are the coupling strengths of the j → k qubit transitions to the cavity mode, a is the cavity annihilation operator, b is the transmon annihilation operator and |k t,f are the eigenstates of the qubits. The fluxonium energy levels and the relative magnitudes of the coupling constants are determined by the device parameters: the Josephson, charging and inductive energies and the flux through the device. In Figure 3 we show how a coherent state with initial amplitude α = 2 evolves both with and without the additional fluxonium. We see that when only the transmon is present the cavity state experiences a full revival after approximately 100 µs, while adding the fluxonium means that the cavity amplitude is held almost constant for this full period, at least an order of magnitude reduction of the self-Kerr. We also show Wigner functions at various times during the evolution along with their fidelity with the initial condition, which is calculated by We see that the self-Kerr in the absence of the fluxonium distorts the state dramatically, producing Schrödinger cat states and other coherent state superpositions before returning to the initial coherent state. When the fluxonium is included, the fidelity remains above 0.90 for 50 µs, before the state starts to become slightly squeezed. This is likely to be due to the presence of higher-order linearities, which are not perfectly cancelled by the fluxonium. In this case, the ratio of the interaction strengths with and without the fluxonium is approximately 10, rather than the theoretical infinite limit given with two Duffing-like qubits.
To store quantum information we can encode it in a superposition of coherent states in the cavity field. Instead of the typical Schrödinger cat states, we may want to use a space of logical computation states which all have positive parity where N + α is a normalisation factor. These logical qubit states are of interest because the loss of a single photon from the state moves the state into an orthogonal subspace to the computational space. This allows easier reconstruction of the state and makes them more suitable for the storage of quantum information [35]. Again, we want to see how these states evolve with and without the fluxonium present. In Fig. 4, we plot the fidelity as a function of time with initial states with a = 0, 0.1, 0.2, 0.3, 0.4, 0.5, noting that the a > 0.5 states are just rotations of the a < 0.5 case. This time fidelities are given by For these superpositions, collapses and revivals occur much more rapidly than for coherent states. While this means that less time must be waited for a full revival, it also means that this time must be calculated very accurately if it is to be corrected for. Introducing the fluxonium means that all states are preserved with greater than 0.50 fidelity for 100 µs, and are preserved with above 90% fidelity for longer than a full revival with just the transmon present. These results show that, by using multiple types of superconducting qubit, we can eliminate the distortions caused by self-Kerr [35] when using the cavity as a quantum memory. While these are deterministic, it may not be practical to wait for the next full revival before performing an operation on the cavity state, for example to switch the state back into a superconducting qubit by a qc-map operation [15]. There is therefore a benefit associated with reducing the cavity self-Kerr as far as possible compared with the decoherence time by design, while maintaining the number splitting that enables the cavity field to be controlled via a coupled qubit. This scheme can be implemented passively rather than requiring relatively complex controls to implement cancellation, such as using specific phase gates [61].

Cross-Kerr engineering for entangling gates between two cavities
Having demonstrated that the cavity self-Kerr can be modified by the use of multiple qubit types, we turn to performing a similar engineering of the cross-Kerr interaction between two cavities. We start by considering two cavities which are both coupled to a single two-level atom with the Hamiltonian where the two cavities have annihilation operators a i , frequencies ω i and coupling to the qubit g i . We can also estimate the cross-Kerr effect using the fourth-order perturbation theory. The approximate eigenvalues in the dispersive limit are For the remainder of this section we will be less concerned with the dispersive shifts of the cavity, and will extract the quadratic correction terms Designing Kerr interactions using multiple superconducting qubit types in a single circuit13 which describe the self-and cross-Kerr interactions. As in the case of self-Kerr on a single mode, considering two qubit levels is not sufficient to describe the qualitative behaviour of the system for weakly anharmonic devices. Adding in a third transmon level produces the new expressions where again χ is the transmon anharmonicity. We see that, as with the self-Kerr the sign of the coupling between the two cavity fields is now determined by the sign of the nonlinearity if the qubit is either above or below both cavities.

Modifying the cross-Kerr interaction with two Duffing-like qubits
We now consider how we can modify the cross-Kerr between the modes by adding a second device, also coupled to both cavities. Specifically, we are interested in whether we can completely switch off the cross-Kerr in some range of parameters. While a single flux-tunable qubit offers some degree of tunability of the frequency, and therefore the cross-Kerr, the range over which the device can be tuned is relatively small. In this configuration, switching off the interaction entirely is difficult as the interaction reduces proportional to ∆ −3 . In an extended system where we want to perform many such operations in turn, even a factor of ten reduction in the cross-Kerr will produce a significant loss of entanglement between the two cavities in the time it takes to perform a small number of gates elsewhere in the system. A way to engineer greater tunability o the interaction is therefore desirable. The Hamiltonian for the two-cavity-two-qubit system is Similar to the self-Kerr calculations in the previous section, the self-and cross-Kerr terms due to both qubits behave additively in the perturbation expansion. The coefficients for this system are X 2C2D = 2g 2 11 g 2 21 χ 1 (∆ 11 + ∆ 21 ) ∆ 2 11 ∆ 2 21 (∆ 11 + ∆ 21 + χ 1 ) where we have defined detunings ∆ ij = ω c,i − ω q,j and couplings g ij between the ith cavity and the j-th qubit. We can now see that, similar to the self-Kerr, we can choose the system parameters so that there is no cross-Kerr, and without making any of the energy levels of the system degenerate. By simplifying the parameters such as g 11 = g 22 , g 12 = g 21 , χ 1 = χ = −χ 2 , ∆ 11 = −∆ 22 , ∆ 12 = −∆ 21 , the self-and cross-Kerr terms are equivalent to

Entangling gates with photons using flux-tunable qubits in two cavities
The effect of the cross-Kerr interaction is to create entanglement between the two cavities. When self-Kerr effect and other nonlinearites are present, additional phases are introduced which change what gate is applied but do not change the degree of entanglement present in the system. A useful way to see this is by using the purity [62] of the state and its individual parts. For a maximally entangled state of two qubits, the full two-qubit system will be pure, whereas tracing out either subsystem will produce a single qubit which is in a completely mixed state.
To perform entangling gates, we now imagine two distinct regimes. In the first, the frequency of the second qubit is selected so that the cross-Kerr is cancelled and the entanglement operation is switched off. In the second the qubit is tuned further from the cavities, X becomes non-zero and entanglement occurs. If this interaction can be switched on for the correct duration, then a maximally entangling operation can be performed. Switching between the two states is achieved by changing the qubit frequency, as can be achieved using an external flux in some transmon devices [63]. To test our perturbation theory results, we fully simulate the evolution of the systems under Figure 6. Plots of purities P (ρ 12 ) in red and P (ρ 1 ) in blue showing the effect of switching the cross-Kerr between mode off and on. P (ρ 2 ) is identical to P (ρ 1 ). (a) If we tune the second qubit according to the values predicted to the perturbation theory calculations ∆ 12 = −1.2 GHz, ∆ 21 = −1 GHz, then the cavities remain in a product state, shown by both P (ρ 12 ) and P (ρ 1 ) remaining close to 1. Using a single qubit which is detuned by a further 1 GHz from both cavities (shown in grey) does not turn off the entangling operation as effectively. (b) When the second qubit is detuned so that ∆ 12 = −2.2 GHz, ∆ 21 = −2 GHz then maximal entanglement is achieved after 30 µs, almost as quickly as would be achieved if the second device was not present. Other system parameters are ∆ 11 = 1 GHz, ∆ 21 = 1200 GHz, g 11 = 80 MHz, g 12 = 87.6 MHz, g 21 = 87.6 MHz, g 22 = 80 MHz, χ 1 = −300 MHz, χ 2 = 300 MHz.
H 2C2T , starting with both cavities in intial state given by a superposition of a vacuum and a photon state (|0 + |1 )/ √ 2. We show the results of these simulations in Figure  6, where we plot the evolution of the system in both the on and off configurations. We plot the purity of the two-qubit state with both cavities traced out P (ρ 12 ) = tr (ρ 12 ) and the purity of the individual qubit states with the remainder of the system traced out P (ρ i ) = tr i (ρ 12 ).
When the qubits are arranged to cancel the cross-Kerr, the system remains in a product state over the full time interval, as shown by the fact that both the two-qubit and single-qubit states have a purity very close to one. The rapid oscillations in the purity are caused by the entanglement of the cavities and qubits that occurs even in the dispersive regime. When one qubit is detuned 1 GHz further from the cavities, we see that the purity of the single qubit states oscillate between 1 and 0.5, with a period of 60 µs. This means that if the interaction is switched on for 30 µs then the two cavity states become maximally entangled. We compare this to the reduction in the cross-Kerr that can be achieved by detuning a single qubit by 1 GHz, giving an on-off ratio for the interaction of 8. We see that the second qubit allows us to switch off the interaction much more effectively, in principle an infinite on-off ratio, without significantly slowing the entanglement in the 'on' configuration. In practise, the range over which a qubit can be tuned is limited and so it is difficult to reduce the cross-Kerr using a single qubit further then shown in Figure 6. We therefore have significant advantages over using a single tunable qubit to implement gates.
As with the self-Kerr interaction, it is likely that the present of non-RWA terms will ultimately restrict how much the cross-Kerr can be reduced but, significantly, the controllability of the cavity field via a qubit is maintained, as there is always a qubit close enough that we remain in the dispersive regime. To perform specific logic operations on the two cavities, it may also be possible to improve this scheme by adding extra qubits to each cavity which cancel the self-Kerr. This ability to perform a cross-Kerr based gate between two modes becomes even more valuable when extended to an array of cavities, where it may enable quantum computation when combined with a single qubit rotation implemented via one of the qubits.

Kerr engineering in cavity arrays
Having demonstrated the principle of Kerr engineering and two applications in smaller systems, we now study a line of N superconducting cavities which are coupled together by intermediary qubits with different nonlinearities. This setup is shown in Figure 7. Each cavity also possesses an on-site qubit, with every device modelled as a Duffing oscillator. The total Hamiltonian for this extended system is where the b operators correspond to the intermediary qubits, while the c operators are for on-site qubits. Each cavity is coupled to three qubits -two hopping qubits to the left and right, with couplings g i and h i , and one in the cavity with coupling f i . Note that the very end qubits are coupled to a single cavity, so we can set g 1 = h N = 0 for the case of open boundary condition or leave them finite for realising the case of periodic boundary conditions. We also define ∆ ij = ω i − ω j the detuning between the ith and jth cavities, Γ ij = ω i − Ω j the coupling between ith cavity and jth qubit and Ξ i = ω i −Ω i the detuning between cavity i and its on-site qubit. Once again, we apply time-independent perturbation theory to find the self-Kerr on each site S i and the cross-Kerr between any two cavities x ij . These are also derived fully in Appendix A. The self-Kerr is simply the sum of the individual qubits contributions and the cross-Kerr between adjacent cavities in the chain is while non-adjacent cavities experience no direct interaction up to fourth order in the coupling. We can see from these expressions that by combining different types of qubits (different signs of χ i , η i ) it is possible to arrange the parameters of the chain such that many combinations of self-and cross-Kerr coefficients can be realised and in particular explore the limit of no self-Kerr S i = 0, X i,i+1 = 0. Importantly, it is also possible to choose them such that the dispersive approximation and RWA are obeyed. The existence of the on-site qubits, whose f i parameter does not feature in the expressions for the cross-Kerr, allows the cross-Kerr to be tuned initially by the intermediary qubits and then the f i used to achieve the desire self-Kerr. This allows us to realise an effective Hamiltonian where ω i are rescaled cavity frequencies, S i are the on-site self-Kerr strengths and X i,i+1 is the cross-Kerr between the i-th and i + 1-th cavities. In practice these results only converge if no two cavities or qubits have the same frequency, so selecting the parameters so that the entire system behaves dispersively is non-trivial. However if we inspect two cavities i, i + n then as n increases their respective resonance frequencies can become closer without affecting the fourth order (Kerr) interaction as they only interact in higher orders. In most realistic designs there would be a significant and random fabrication uncertainties but with the advent of tunable qubits and tunable cavities, it is possible to compensate small discrepancies in the nonlinearities by in-situ controlling Γ and Ξ. This ability to sweep the parameters means that cancellation may be possible even though measuring the bare parameters in an extended is extremely challenging. In addition, parameters of the system may drift or fluctuate in time, necessitating detailed experimental investigation. In principle, the potential of combining different types of qubits can also be extended 2D lattices, although this presents even greater complexity.

Conclusion
We have shown the benefit of combining different types of superconducting qubit to engineer Kerr nonlinearities in systems of cavities. Using perturbation theory and numerical simulations of the exact model we have shown that coupling two types of qubit, with opposite anharmonicities, to a single cavity can enable a complete cancellation of the cavity self-Kerr. When one of the qubits is tunable then the Kerr interaction becomes also tunable. Despite the cancellation of the self-Kerr interaction, the cavity field can still be controlled by the number splitting of the qubit due to the coupling to the cavity. We then showed that a similar result can be achieved for cross-Kerr between modes, demonstrating that configurations exist where the cross-Kerr can be switched on and off, enabling a maximally entangling operation to implemented between the two cavities if frequency-tunable qubits are used. This is a significant improvement over the efficacy of coupling the cavity to a single tunable qubit which be reduces the interaction strength polynomially and therefore cannot remove it entirely without leaving the dispersive regime. Finally, we extended this to a line of superconducting cavities and showed that it could be possible to engineer a chain of cavities with arbitrary on-site and interaction Kerr-coefficients. The on-site qubit for each cavity then allows additional tunability of a class of many-body Hamiltonians and can be used to study phase transitions and quenches out of equilibrium. Implementing cross-Kerr mediated gates in a line of superconducting cavities could also provide a platform for quantum computing.

Acknowledgements
where λ is a small parameter compared with the terms in the unperturbed Hamiltonian H 0 , we can write the first four orders of energy corrections in terms of the unperturbed eigenstates |k and eigenvalues E (0) k [56]. We define V nm = n|V |m and E n m = E In the systems we are looking at, with only linear couplings between cavities and qubits in V , applying V to any eigenstate of the unperturbed system will take the system to an orthogonal state, so V nn = 0. If we are only considering systems without any closed loops, it is also impossible to have terms such as V nk 2 V k 2 k 1 V k 1 n , with odd numbers of matrix elements. This is because an odd number of applications of the linear coupling term cannot map you back to the original state. This greatly simplifies the perturbed energies to Appendix A.1. Cross-Kerr interaction with two-level qubits One case we apply perturbation theory to is two cavities that are both coupled to a single qubit. The Hamiltonian in this case is The second order correction to the cavity energy levels includes only terms from paths excitations can take through the system by applying the interaction part of the Hamiltonian twice and that return to the initial configuration. The only such paths are one excitation for either cavity hopping to the qubit and back again, so the correction is given by The fourth order correction includes all paths made up of four applications of the Hamiltonian, where every intermediary state is distinct from the initial (and final) state. The first two terms come from an excitation hopping from one cavity to the other and back again, with opposite signs because they occur in opposite directions.
E 2C2Q(4) n 1 ,n 2 ,g = g 2 1 g 2 2 n 1 (n 2 + 1) ∆ 2 1 ∆ 12 − g 2 1 g 2 2 n 2 (n 1 + 1) We can see both cross-Kerr and self-Kerr terms in the fourth order correction. Without considering additional transmon levels, there exists a configuration where both the selfand cross-Kerr can in principle be eliminated by an appropriate choice of the parameters.

Appendix A.2. Cross-Kerr with transmon levels
If instead of a qubit we use a transmon modelled as a Duffing oscillator, then the Hamiltonian for a single cavity and qubit is There are no additional corrections to the cavity energy levels, with the qubit in the ground state, at second order because it is not possible to access the higher qubit levels and return to the initial state with two applications of the interaction Hamiltonian, so At fourth order, however, two excitations can hop into the transmon, giving us the fourth order correction Now in the two cavity-one transmon case the Hamiltonian is Again, the second order correction is unchanged by the extra transmon levels but there are new terms at the fourth order correction. These are associated with two photons from one qubit going into the qubit and back again, and also terms where one from each cavity goes into the qubit and then return to the cavities, in four different orders. The fourth order correction is therefore E 2C1T (4) n 1 ,n 2 ,g = g 2 1 g 2 2 n 1 (n 2 + 1) ∆ 2 1 ∆ 12 − g 2 1 g 2 2 n 2 (n 1 + 1) We can repeat this with two Duffing-like devices to gives the self-and cross-Kerr coefficients (A.13) X 2C2D = 2g 2 11 g 2 21 χ 1 (∆ 11 + ∆ 21 ) ∆ 2 11 ∆ 2 21 (∆ 11 + ∆ 21 + χ 1 ) + 2g 2 12 g 2 22 χ 2 (∆ 12 + ∆ 22 ) ∆ 2 12 ∆ 2 22 (∆ 12 + ∆ 22 + χ 2 ) . (A.14) When considering for than two qubit levels, it is no longer possible to eliminate selfand cross-Kerr simultaneously.

Appendix A.3. Array of cavities
Finally, we consider a complete line of cavities, each connected by intermediary qubits, and with each cavity possessing its own on-site qubit. The total Hamiltonian using only two-level qubits is i Ω j 2 σ z j + The second order correction is as for the simpler cases above, simply summing the contributions due to each qubit Designing Kerr interactions using multiple superconducting qubit types in a single circuit22 For the fourth order contribution, there are terms associated with: excitations in each cavity hopping to the cavities to the left and right; two excitations in the same cavity hopping onto different qubits; and two excitations from different points in the line hopping onto different qubits. The full correction is n i (n i − 1) n 2 i .

(A.17)
As before we can now extract the self-Kerr from the new energies n 2 i , (A. 18) and we see that the total nonlinearity is given by adding the separate induced nonlinearities due to the individual qubits. The cross-Kerr between adjacent cavities is only affected by the qubit coupled between them The effect of the new non-RWA terms is to add a further correction to the energy shifts, which to second order in the interaction are E 1C1T (2) n,g = g 2 n ∆ + g 2 (n + 1) ω c + ω q . (A.24) At fourth order, these new interaction terms add a large number of additional perturbation terms E 1C1T (4) n,g = − g 4 n 2 ∆ 3 + 2 g 4 n(n − 1) ∆ 2 (2∆ + χ) − g 4 n(n − 1) ∆ 2 (2ω c + χ) + 2g 4 n 2 ∆ 2 (2ω q + χ) + 2g 4 (n + 1) 2 (2ω q + χ)(ω c + ω q ) 2 + 2g 4 (n + 1)(n + 2) (ω c + ω q ) 2 (2ω c + 2ω q + χ) + g 4 (n + 1)(n + 2) 2ω c (ω c + ω q ) 2 − g 4 (n + 1)(n + 2) 2ω c ∆(ω c + ω q ) (A. 25) In the regime ∆ |ω c + ω q |, it must also be the case that individually ω c , ω q ∆. This means that all the new terms that are not seen in the JC model are very small in this limit. For typical experiments, resonators and qubits in the 5-10 GHz range are used, meaning that with with detunings of the order of 1 GHz these terms will be at least a