Parametric conversion and amplification based on three-wave mixing are powerful primitives for efficient quantum operations. For superconducting qubits, such operations can be realized with a quadrupole Josephson junction element, the Josephson Ring Modulator, which behaves as a loss-less three-wave mixer. However, combining multiple quadrupole elements is a difficult task so it would be advantageous to have a three-wave dipole element that could be tessellated for increased power handling and/or information throughput. Here, we present a dipole circuit element with third-order nonlinearity, which implements three-wave mixing. Experimental results for a non-degenerate amplifier based on the proposed third-order nonlinearity are reported.

In quantum devices based on superconductors, Josephson junctions provide a nonlinear interaction between electromagnetic modes, which is purely dispersive. However, because the Josephson potential is an even function of the superconducting phase difference φ, this nonlinearity is, to lowest order, of the form φ4. This Kerr nonlinearity is useful for engineering interactions between modes, but it imparts undesired frequency shifts. This problem is analogous to that generated by χ(3) media in nonlinear optics.1 Such frequency shifts become especially problematic as the number of interacting modes, and therefore frequency crowding, increases. An alternative strategy is to use a minimal φ3 nonlinearity for engineering the same useful interactions between modes while minimizing these unwanted Kerr frequency shifts.

A form of this φ3 nonlinearity has been realized with four Josephson junctions arranged in a ring threaded by DC magnetic flux. This device, called the Josephson ring modulator (JRM),2,3 specifically provides a trilinear Hamiltonian term of the form φXφYφZ between three modes named X, Y, and Z. However, the JRM is a quadrupole element, i.e., it imposes a current and phase relation between four nodes of a circuit. Thus, a question naturally arises, would it be possible to engineer a third-order φ3 nonlinearity in a dipole device similar to the Josephson junction?

An inherent advantage of such a dipole element is that it is more modular and easier to integrate into complex circuits. In particular, such an element can be tessellated, as with a Josephson junction,4 for improved power handling capabilities and information throughput. Furthermore, it can provide a three-wave coupling between a superconducting qubit and an oscillator.

In this letter, we report an experimental realization of a pure φ3 nonlinear dipole element, demonstrating its performance in a non-degenerate parametric amplifier. Because this dispersive element is asymmetric in the transformation of φφ, in contrast with the superconducting quantum interference device (SQUID)5,6 and the superconducting low-inductance undulatory galvanometer (SLUG),7,8 we have named it the Superconducting Nonlinear Asymmetric Inductive eLement (SNAIL). After the completion of our work, we became aware of a theoretical proposal for a similar element, given in Ref. 9.

Figure 1(a) shows the circuit we propose to achieve cubic nonlinearity in a dipole circuit element. The SNAIL consists of a superconducting loop of n large Josephson junctions and a single smaller junction (tunneling energies EJ and αEJ, respectively), which we thread with DC magnetic flux Φext. The proposed circuit has the inductive energy:

USNAIL(φ)=αEJcos(φ)nEJcos(φextφn),
(1)

where φ is the superconducting phase across the small junction, φext=2πΦext/Φ0 is the reduced applied magnetic flux, and Φ0=h/2e is the magnetic flux quantum. Note that Eq. (1) is only a function of a single degree of freedom φ as we eliminated the dynamics due to any intra-array modes and considered only common excitations across the array of n junctions. This reduction is valid when EJEC for each junction where EC=e2/2CJ is the Coulomb charging energy of the junction with capacitance CJ, and when C0CJ/n2 where C0 is the capacitance to ground of each island between junctions.10 We will assume in the following that these conditions are realized.

FIG. 1.

(a) Circuit for the Superconducting Nonlinear Asymmetric Inductive eLement (SNAIL) reduced to one degree of freedom φ. The loop of n = 3 large junctions and one smaller junction (tunneling energies EJ and αEJ, respectively) is threaded with an external DC flux Φext. (b) An example SNAIL potential for α=0.29 and Φext=0.41Φ0 that includes third-order nonlinearity without fourth-order nonlinearity. (c) Colormaps of the (α,Φext) parameter space for the third-order (top) and fourth-order (bottom) nonlinear terms show regions where c30 while c4=0 [white region in bottom panel (c)]. The black-hatched regions in (c) correspond to unwanted hysteretic double well behavior similar to the flux qubit regime. The crosses in (c) mark the set of parameters chosen for (b). This working point is optimal (see text).

FIG. 1.

(a) Circuit for the Superconducting Nonlinear Asymmetric Inductive eLement (SNAIL) reduced to one degree of freedom φ. The loop of n = 3 large junctions and one smaller junction (tunneling energies EJ and αEJ, respectively) is threaded with an external DC flux Φext. (b) An example SNAIL potential for α=0.29 and Φext=0.41Φ0 that includes third-order nonlinearity without fourth-order nonlinearity. (c) Colormaps of the (α,Φext) parameter space for the third-order (top) and fourth-order (bottom) nonlinear terms show regions where c30 while c4=0 [white region in bottom panel (c)]. The black-hatched regions in (c) correspond to unwanted hysteretic double well behavior similar to the flux qubit regime. The crosses in (c) mark the set of parameters chosen for (b). This working point is optimal (see text).

Close modal

With α0.8 and Φext0.5Φ0, this circuit is well known as the flux qubit,11,12 which has a double-well potential. However, here we propose using a different parameter set to create a potential with a single minimum. Moreover, we can adjust the potential to cancel the fourth-order (Kerr) term while keeping a substantial cubic term for a particular choice of α and Φext.

Through numerical minimization of Eq. (1), we analyze the SNAIL's mixing capabilities by Taylor expanding about the minimum φmin to obtain the effective potential for φ̃=φφmin:

Ueff(φ̃)/EJ=c2φ̃2+c3φ̃3+c4φ̃4+,
(2)

where (c2, c3, c4, ) are numerically determinable coefficients whose specific values depend on n, α, and Φext. Our desire for a pure cubic nonlinearity without any Kerr translates to a Ueff with nonzero c3, but c4=0. Note that the case of n = 1, corresponding to the SQUID, always gives c3=0 since the potential about φmin is always a pure cosine irrespective of the values of α and Φext. Interestingly, α and Φext are not enough to break the φ̃φ̃ symmetry in this case. Additionally, in the n1 limit, the array behaves as a linear inductance and the potential approaches the fluxonium qubit/RF-SQUID regime.13,14 Replacing the array with a geometric inductance in an RF-SQUID configuration also achieves the desired nonlinearity.9 Experimentally, we choose the smallest n2 that is easy to fabricate, which is n = 3 for our Dolan bridge process. We specialize all further analysis to the case n = 3 as depicted in Fig. 1(a), but the results are easily extendable to different values of n.

For the SNAIL to achieve Kerr-free three-wave mixing, we must optimize on the parameter space of (α, Φext). The top (bottom) panel of Fig. 1(c) shows c3 (c4) as a function of the (α, Φext) parameter space. Focusing on the top of Fig. 1(c), we want to maximize c3 while avoiding any hysteretic double-well behavior, marked by the black-hatched region. This is done by restricting ourselves to α0.5. The Kerr-free region of parameter space corresponds to c4=0 [the white region in the bottom panel of Fig. 1(c)]. As such, we choose an optimal set of parameters α=0.29 and Φext=0.41Φ0, marked by the crosses in Fig. 1(c), which maximizes c3 while c4=0 and the potential has a single minimum. We plot one period of the USNAIL potential [Eq. (1)] for this optimal set of parameters in Fig. 1(b).

In order to experimentally verify our model for the SNAIL as a cubic nonlinearity, we integrate it into a design for a non-degenerate three-wave mixing amplifier. Such an amplifier consists of two spatially and spectrally separate modes that are coupled by a pump applied at the sum of their respective frequencies. Achieving amplification with this pump condition is an explicit signature of a cubic nonlinearity.1 The design [Fig. 2(a)] is based on that of the Josephson Ring Modulator (JRM),2,3 a quadrupole element with four identical Josephson junctions in a superconducting ring. By threading the ring with an external magnetic flux, the JRM achieves the three-wave coupling φXφYφZ between modes X, Y, and Z. This coupling allows the device to act as a parametric amplifier and frequency converter.

FIG. 2.

(a) A circuit design for a non-degenerate three-wave mixing amplifier which uses the SNAIL and is based on the design of the JRM. The elements in the dashed blue box represent a coupled microstrip line whose even mode electrical length is λS/4. The coupled microstrip is open at DC, and so, unlike the JRM, there is no exterior loop in this design. Depictions of the X, Y, and Z eigenmodes of the device are shown by the diamonds on the right. The + and − signs refer to the relative phase of the voltage across the four SNAILs. In (b), we show an optical image of the four SNAIL mixing structure, where the inset shows an individual SNAIL with n = 3 large junctions. The arrows highlight the four junctions. For discussion of the black veil around the junctions, see text.

FIG. 2.

(a) A circuit design for a non-degenerate three-wave mixing amplifier which uses the SNAIL and is based on the design of the JRM. The elements in the dashed blue box represent a coupled microstrip line whose even mode electrical length is λS/4. The coupled microstrip is open at DC, and so, unlike the JRM, there is no exterior loop in this design. Depictions of the X, Y, and Z eigenmodes of the device are shown by the diamonds on the right. The + and − signs refer to the relative phase of the voltage across the four SNAILs. In (b), we show an optical image of the four SNAIL mixing structure, where the inset shows an individual SNAIL with n = 3 large junctions. The arrows highlight the four junctions. For discussion of the black veil around the junctions, see text.

Close modal

To achieve the same coupling term, we arrange four identical SNAILs into a ring. Like the JRM, this Wheatstone bridge of SNAILs cancels all cubic terms in the Hamiltonian of the device except the desired φXφYφZ. This cancellation assumes that all four SNAILs are identical. Note that inverting the polarity of one SNAIL would destroy this cancellation, since Ueff [Eq. (2)] depends on the sign of φ̃. While the JRM relies on the flux through the exterior ring to achieve a cubic potential, in our design each SNAIL individually achieves φ̃3, and so we can eliminate the exterior DC loop in the ring. This is advantageous as it prevents an extra DC persistent current from biasing the SNAILs and allows for future tessellation without introducing a prohibitively large DC loop. Moreover, the SNAIL allows flexibility for other applications currently unfeasible with the JRM, such as endowing a single mode with a cubic nonlinearity. Further applications may also prefer the expression of terms such as a φX2φZ for degenerate three-wave mixing,15 achievable by simply inverting the polarity of two opposite SNAILs in the ring. In this way, the SNAIL allows the circuit designer to separate the size of the loop for flux biasing from the design choice of which third-order terms are expressed and which are canceled.

For enhancement of the nonlinearity and maximal isolation, we designed a resonant embedding structure using microstrip resonators,16,17 which couple to the X, Y, and Z modes of the SNAIL ring [Fig. 2(a)]. The highest frequency mode, called the idler, couples to the Y mode of the ring via two λI/4 sections of standard microstrip transmission line (red cylinders), where λI is the wavelength of the idler mode. For the X mode of the ring, called the signal, on one side we use similar coupling through a λS/4 section of standard microstrip transmission line (blue cylinder). On the other side, however, we use an impedance-matched coupled microstrip transmission line (blue-boxed elements), which breaks the exterior loop at DC but allows RF coupling to φX. Gap capacitors couple the signal and idler modes to 50 Ω transmission lines to set the quality factor Q150 for this device. A small coupling capacitor weakly couples our pump port to the Z mode of the ring to allow driving of parametric processes for three-wave mixing.

Using standard Dolan bridge techniques for aluminum Josephson junctions on silicon, we fabricated a ring of four SNAILs according to the above design [Fig. 2(b)]. The inset of Fig. 2(b) shows an electron micrograph of a single SNAIL element with n = 3 large junctions (critical currents I0=7.1μA) and a single smaller junction (I0=2.0μA). The ratio of junction sizes gave α=0.29 as desired. For the array of large junctions, we utilized a technique developed for the original fluxonium arrays13,14 that allows for creation of three equivalent junctions using only two Dolan bridges. Characteristics of the standard microstrip embedding structure were similar to past work.17 The black veil surrounding the junctions in the micrograph results from undesired leftover resist in those locations during deposition.18 

For experimental characterization, we mounted the device at the base stage of a helium dilution refrigerator with access to both the signal and idler ports. A uniform external magnetic flux Φext was applied to each SNAIL by a magnet coil mounted under the sample. Figure 3(a) shows the resonant frequency fS (fI) of the signal (idler) mode as a function of Φext, determined through reflection measurements off of the signal (idler) port (example shown in the inset). Sweeping Φext in the opposite direction resulted in the same mode frequencies, implying that our device was not hysteretic. The dashed black line corresponds to the calculated optimal flux Φext=0.41Φ0 for this sample with α=0.29.

FIG. 3.

(a) Resonant frequencies of the Signal (red) and Idler (blue) modes as a function of external applied flux. The inset shows the phase of the reflection coefficient at Φext=0.41Φ0 (black dashed line), used to fit resonant frequency at each flux. Biasing at the dashed line in (a) and applying a strong microwave pump tone at fP=fS+fI, the top panel of (b) shows the reflection gain off the Idler port. Gain increased and bandwidth decreased as a function of pump power (indicated by color). As a proxy for noise temperature, the bottom panel of (b) shows the Noise Visibility Ratio (NVR=Pon/Poff). The NVR of our device was comparable to that of other quantum-limited amplifiers measured in the same system, suggesting that the noise performance of the SNAIL is similar to that of the nearly quantum-limited JRM.

FIG. 3.

(a) Resonant frequencies of the Signal (red) and Idler (blue) modes as a function of external applied flux. The inset shows the phase of the reflection coefficient at Φext=0.41Φ0 (black dashed line), used to fit resonant frequency at each flux. Biasing at the dashed line in (a) and applying a strong microwave pump tone at fP=fS+fI, the top panel of (b) shows the reflection gain off the Idler port. Gain increased and bandwidth decreased as a function of pump power (indicated by color). As a proxy for noise temperature, the bottom panel of (b) shows the Noise Visibility Ratio (NVR=Pon/Poff). The NVR of our device was comparable to that of other quantum-limited amplifiers measured in the same system, suggesting that the noise performance of the SNAIL is similar to that of the nearly quantum-limited JRM.

Close modal

Biasing at this optimal flux point, we applied a strong microwave pump tone on the pump port at frequency fP=fS+fI to drive the parametric amplification process. Achieving parametric amplification with this frequency condition provides the explicit signature of three-wave mixing and hence the presence of a φ3 nonlinearity. The top panel of Fig. 3(b) shows the reflection gain off of the signal port. Increasing the pump power (color) increased the gain and decreased the bandwidth in accordance with the gain-bandwidth product, standard for cavity-based parametric amplifiers.19 We also observed a slight shift in frequency for the maximum gain as we increased the pump power. This shift is similar in magnitude to that observed with a JRM. Further study of this shift and its relation to the Kerr nonlinearity20 are left for future work. In the bottom panel of Fig. 3(b), we plot the Noise Visibility Ratio, a proxy for noise temperature, which is the ratio between the noise power spectral density with the pump on and the pump off (NVR=Pon/Poff). The NVR of our device was comparable to that of JRM-based quantum-limited amplifiers measured in the same system,2 showing that the SNAIL can potentially be used as a quantum-limited amplifier as well.

In future work, the dipole nature of the SNAIL enables the replacement of a single SNAIL with an array of SNAILs [Fig. 4(a)] for increased dynamic range or bandwidth. This strategy is directly analogous to the tessellation of DC-SQUIDs in the JPA4 or the tessellation of transmission line unit cells in the Josephson traveling wave parametric amplifier,21 both four-wave mixing amplifiers. The SNAIL allows this tessellation for three-wave mixing amplifiers.9 

FIG. 4.

(a) Replacing a SNAIL with an array of SNAILs will allow for increased dynamic range or bandwidth, similar to the tessellation of DC-SQUIDs in the JPA. The dipolar nature of the SNAIL allows this tessellation for three-wave mixing amplifiers. (b) An example of a circuit in which the SNAIL is used as a tunable coupling device between two modes (red and blue boxes).

FIG. 4.

(a) Replacing a SNAIL with an array of SNAILs will allow for increased dynamic range or bandwidth, similar to the tessellation of DC-SQUIDs in the JPA. The dipolar nature of the SNAIL allows this tessellation for three-wave mixing amplifiers. (b) An example of a circuit in which the SNAIL is used as a tunable coupling device between two modes (red and blue boxes).

Close modal

Another relevant application of the SNAIL is for Kerr-free tunable coupling between two cavities [Fig. 4(b)]. By pumping at the difference frequency between the red and blue modes, the three-wave mixing capabilities of the SNAIL generate the effective interaction geff(t)(abar+arab), where ar and ab are the annihilation operators for the red and blue modes respectively. Such tunable couplings have been used as tools for quantum optics in the microwave regime;22,23 however, previous schemes always required that the red and blue modes inherit self-Kerr and cross-Kerr terms from the nonlinearity.24 In principle, using a SNAIL with c4=0 allows one to generate the same effective interaction without the harmful Kerr terms.

In conclusion, we have introduced a dipole realization of a pure φ3 nonlinearity, the SNAIL, which performs three-wave mixing in the microwave regime of superconducting quantum circuits. We have also demonstrated the SNAIL's three-wave mixing capabilities by integrating it into a non-degenerate three-wave mixing amplifier, whose performance matched that of similar quantum-limited amplifiers based on the JRM. However, the SNAIL has the advantage of supporting parametric three-wave processes without introducing extra Kerr-induced frequency shifts. Left for future work is the demonstration of this feature. Furthermore, the dipole nature of the SNAIL allows for the integration of cubic nonlinearity wherever the circuit designer deems necessary. Specifically, this enables the construction of three-wave mixing circuits that share the work between many junctions to increase bandwidth and input power handling capabilities, without the addition of prohibitively large DC loops. Thus, the SNAIL will likely be at the heart of the next generation of robust quantum-limited parametric amplifiers and converters.

We acknowledge fruitful discussions with Clarke Smith and Luigi Frunzio. We also acknowledge the Yale Quantum Institute. Facilities use was supported by YINQE, the Yale SEAS cleanroom, and NSF MRSEC DMR 1119826. This research was supported by ARO under Grant Nos. W911NF-14-1-0011 and W911NF-16-10349, AFOSR under Grant No. FA9550-15-1-0029, and MURI-ONR Grant No. N00014-16-1-2270.

1.
R. W.
Boyd
,
Nonlinear Optics
, 3rd ed. (
Elsevier Inc
.,
2008
).
2.
N.
Bergeal
,
F.
Schackert
,
M.
Metcalfe
,
R.
Vijay
,
V. E.
Manucharyan
,
L.
Frunzio
,
D. E.
Prober
,
R. J.
Schoelkopf
,
S. M.
Girvin
, and
M. H.
Devoret
,
Nature
465
,
64
(
2010
).
3.
N.
Bergeal
,
R.
Vijay
,
V. E.
Manucharyan
,
I.
Siddiqi
,
R. J.
Schoelkopf
,
S. M.
Girvin
, and
M. H.
Devoret
,
Nat. Phys.
6
,
296
(
2010
).
4.
M. A.
Castellanos-Beltran
,
K. D.
Irwin
,
G. C.
Hilton
,
L. R.
Vale
, and
K. W.
Lehnert
,
Nat. Phys.
4
,
929
(
2008
).
5.
J. E.
Zimmerman
and
A. H.
Silver
,
Phys. Rev.
141
,
367
(
1966
).
6.
The SQUID Handbook
, edited by
J.
Clarke
and
A. I.
Braginski
(
Wiley
,
Weinheim
,
2004
), Vol. 1.
7.
J.
Clarke
,
Philos. Mag.
13
,
115
(
1966
).
8.
D.
Hover
,
Y.-F.
Chen
,
G. J.
Ribeill
,
S.
Zhu
,
S.
Sendelbach
, and
R.
McDermott
,
Appl. Phys. Lett.
100
,
063503
(
2012
).
9.
A.
Zorin
,
Phys. Rev. Appl.
6
,
034006
(
2016
).
10.
K. A.
Matveev
,
A. I.
Larkin
, and
L. I.
Glazman
,
Phys. Rev. Lett.
89
,
096802
(
2002
).
11.
J. E.
Mooij
,
T. P.
Orlando
,
L.
Levitov
,
L.
Tian
,
C. H.
van der Wal
, and
S.
Lloyd
,
Science
285
,
1036
(
1999
).
12.
C. H.
van der Wal
,
A. C. J.
ter Haar
,
F. K.
Wilhelm
,
R. N.
Schouten
,
C. J. P. M.
Harmans
,
T. P.
Orlando
,
S.
Lloyd
, and
J. E.
Mooij
,
Science
290
,
773
(
2000
).
13.
V. E.
Manucharyan
,
J.
Koch
,
L. I.
Glazman
, and
M. H.
Devoret
,
Science
326
,
113
(
2009
).
14.
N. A.
Masluk
, “
Reducing the losses of the fluxonium artificial atom
,” Ph.D. thesis,
Yale University
,
2012
.
15.
T.
Yamamoto
,
K.
Inomata
,
M.
Watanabe
,
K.
Matsuba
,
T.
Miyazaki
,
W. D.
Oliver
,
Y.
Nakamura
, and
J. S.
Tsai
,
Appl. Phys. Lett.
93
,
042510
(
2008
).
16.
B.
Abdo
,
A.
Kamal
, and
M.
Devoret
,
Phys. Rev. B
87
,
014508
(
2013
).
17.
F.
Schackert
, “
A practical quantum-limited parametric amplifier based on the josephson ring modulator
,” Ph.D. thesis,
Yale University
,
2013
.
18.

In principle, this undesired material places an upper bound on the internal quality factor of our device, but this bound is significantly higher than the external quality factor. This black veil is likely to be the result of the resist dose associated with our large junctions.

19.
A. A.
Clerk
,
M. H.
Devoret
,
S. M.
Girvin
,
F.
Marquardt
, and
R. J.
Schoelkopf
,
Rev. Mod. Phys.
82
,
1155
(
2010
).
20.
G.
Liu
,
T.-C.
Chien
,
X.
Cao
,
O.
Lanes
,
E.
Alpern
,
D.
Pekker
, and
M.
Hatridge
, preprint arXiv:1703.04425 [cond-mat, physics:quant-ph] (
2017
).
21.
C.
Macklin
,
K.
O'Brien
,
D.
Hover
,
M. E.
Schwartz
,
V.
Bolkhovsky
,
X.
Zhang
,
W. D.
Oliver
, and
I.
Siddiqi
,
Science
350
,
307
(
2015
).
22.
E.
Flurin
,
N.
Roch
,
J.
Pillet
,
F.
Mallet
, and
B.
Huard
,
Phys. Rev. Lett.
114
,
090503
(
2015
).
23.
K. M.
Sliwa
,
M.
Hatridge
,
A.
Narla
,
S.
Shankar
,
L.
Frunzio
,
R. J.
Schoelkopf
, and
M. H.
Devoret
,
Phys. Rev. X
5
,
041020
(
2015
).
24.
W.
Pfaff
,
C. J.
Axline
,
L. D.
Burkhart
,
U.
Vool
,
P.
Reinhold
,
L.
Frunzio
,
L.
Jiang
,
M. H.
Devoret
, and
R. J.
Schoelkopf
, preprint arXiv:1612.05238 [cond-mat, physics:quant-ph] (
2016
).