Superconducting circuit protected by two-Cooper-pair tunneling

TL;DR

This paper introduces a novel superconducting qubit design that leverages two-Cooper-pair tunneling to create a protected quantum state, potentially enhancing qubit coherence by reducing environmental sensitivity.

Contribution

The work proposes a new superconducting circuit element enabling parity-based qubit encoding with improved protection against decoherence.

Findings

Ground state manifold is nearly degenerate and indexed by Cooper pair parity

Large capacitance shunting enhances qubit protection

Potential for increased qubit coherence and stability

Abstract

We present a protected superconducting qubit based on an effective circuit element that only allows pairs of Cooper pairs to tunnel. These dynamics give rise to a nearly degenerate ground state manifold indexed by the parity of tunneled Cooper pairs. We show that, when the circuit element is shunted by a large capacitance, this manifold can be used as a logical qubit that we expect to be insensitive to multiple relaxation and dephasing mechanisms.

Figures

Tables

Table 1. Table 1: Circuit parameters used for numerical simulations.

${ϵ}_J/h$	${ϵ}_C/h$	${ϵ}_L/h$	$x$
15 GHz	2 GHz	1 GHz	0.02

Table 2. Table 2: Expected relaxation times T 1 subscript 𝑇 1 T_{1} and pure dephasing times T ϕ subscript 𝑇 ϕ T_{\upphi} at φ ext = π subscript 𝜑 ext 𝜋 \varphi_{\text{ext}}=\pi through various channels. The operators 𝒪 𝒪 \mathcal{O} coupling to the bath and the noisy parameters λ 𝜆 \lambda are listed in addition to the associated quality factors and spectral density amplitudes, respectively. Relaxation and dephasing times are shown for varying cases of inductive disorder: δ L = 0.0 , 0.3 , 0.6 , 0.9 subscript 𝛿 𝐿 0.0 0.3 0.6 0.9 \delta_{L}=0.0,0.3,0.6,0.9 (see also Fig.  5 ). The operators N i subscript 𝑁 𝑖 N_{i} depict the number of Cooper pairs having tunneled across the two Josephson junctions, i.e.  N i = n ± 1 2 ​ ( N − η ) subscript 𝑁 𝑖 plus-or-minus 𝑛 1 2 𝑁 𝜂 N_{i}=n\pm\frac{1}{2}(N-\eta) Entries that read “ ∞ \infty ” represent numerical infinity. *These values are reported for clarity but have no bearing on the estimates shown (see Secs.  III.3 , III.4 ).

Loss channel	$𝒪$	Quality factor	$T_1$ (ms)
			${\delta }_L=0.0$	${\delta }_L=0.3$	${\delta }_L=0.6$	${\delta }_L=0.9$
Capacitive	$2e{N}_i$	$Q_{\text{cap}}\sim {\mathrm{1\times 10}}^6\text{<cite>Wang <em>et al.</em> (<a>2015</a>)</cite>}$	$\mathrm{780\hspace{0.17em}000}$	$\mathrm{17\hspace{0.17em}000}$	$\mathrm{1\hspace{0.17em}000}$	$18$
Inductive	${\varphi }_0{\varphi }_i$	$Q_{\text{ind}}\sim {\mathrm{500\times 10}}^6\text{<cite>Pop <em>et al.</em> (<a>2014</a>)</cite>}$	$0.61$	$0.79$	$1.1$	$1.4$
Purcell	$2e\eta$	$Q_{\text{cap}}\sim {\mathrm{1\times 10}}^6\text{<cite>Wang <em>et al.</em> (<a>2015</a>)</cite>}$	$\infty$	$\mathrm{2\hspace{0.17em}500}$	$380$	$470$
Quasiparticle	$2{\varphi }_0\sin \frac{{\phi }_i}{2}$	$1/x_{\text{qp}}^{\ast }\sim {\mathrm{0.3\times 10}}^6\text{<cite>Pop <em>et al.</em> (<a>2014</a>)</cite>}$	$\infty$	$\infty$	$\infty$	$\infty$
Dephasing channel	$\lambda$	Spectral density amplitude	$T_{\varphi }$ (ms)
			${\delta }_L=0.0$	${\delta }_L=0.3$	${\delta }_L=0.6$	${\delta }_L=0.9$
Charge	$N_{\text{g}}$	${\sqrt{A_{N_{\text{g}}}}}^{\ast }\sim {\mathrm{1\times 10}}^{-4}\text{<cite>Zorin <em>et al.</em> (<a>1996</a>)</cite>}$	$0.0037$	$0.15$	$74$	$3.3\times {10}^6$
Flux	${\phi }_{\text{ext}}$	$\sqrt{A_{{\phi }_{\text{ext}}}}/2\pi \sim {\mathrm{3\times 10}}^{-6}\text{<cite>Quintana <em>et al.</em> (<a>2017</a>)</cite>}$	$0.022$	$0.13$	$0.67$	$1.8$
Shot	$n_{\text{p}}$	$n_{\text{th}}/Q_{\text{cap}}\sim {\mathrm{1\times 10}}^{-7}\text{<cite>Wang <em>et al.</em> (<a>2015</a>)</cite>}$	$4.6$	$4.8$	$5.3$	$8.7$
Critical current	${ϵ}_J$	$\sqrt{A_{{ϵ}_J}}/{ϵ}_J\sim {\mathrm{5\times 10}}^{-7}\text{<cite>van Harlingen <em>et al.</em> (<a>2004</a>)</cite>}$	$210$	$40$	$8.2$	$2.9$

Equations

$$-\frac{1}{2}{E}_J\sum _{N=-\infty }^{\infty }\left(|N⟩⟨N+2|+|N+2⟩⟨N|\right)=-E_J\cos 2\phi ,$$

$$H=4E_C(N-N_{\text{g}}{)}^2-E_J\cos 2\phi ,$$

$$\Delta E\approx 16E_C\sqrt{\frac{2}{\pi }}{\left(\frac{2E_J}{E_C}\right)}^{3/4}{\text{e}}^{-\sqrt{2E_J/E_C}}\cos (\pi N_{\text{g}}),$$

$$\varphi ={\phi }_1+{\phi }_2\phi =\frac{1}{2}({\phi }_1-{\phi }_2)\theta =\frac{1}{2}({\varphi }_1-{\varphi }_2),$$

$$H$$

$$+{ϵ}_L\left[\frac{1}{4}(\varphi -{\phi }_{\text{ext}}{)}^2+{\theta }^2\right]-2{ϵ}_J\cos \phi \cos \frac{\varphi }{2},$$

$$\varphi =\frac{1}{1+z}\left(2\left|\phi -2\pi \text{round}\frac{\phi }{2\pi }\right|+z{\phi }_{\text{ext}}\right),$$

$$H_{\text{eff}}$$

$$-\frac{16}{3\pi }{ϵ}_L(\pi -{\varphi }_{\text{ext}})\cos \phi -{ϵ}_J\left(1-\frac{5}{4}z\right)\cos 2\phi$$

$$\text{∙}/\text{∘}=\{\begin{array}{ll}\text{leavevmode}to10.09pt\text{vbox}to6.86pt\text{pgfpicture}\text{makeatletter}\text{lower}-0.2ptto0.0pt\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{definecolor}pgfstrokecolorrgb\mathrm{0,0,0}\text{pgfsys@color@rgb@stroke}000\text{pgfsys@invoke}\text{pgfsys@color@rgb@fill}000\text{pgfsys@invoke}\text{pgfsys@setlinewidth}0.4pt\text{pgfsys@invoke}\text{nullfont}to0.0pt\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{pgfsys@transformcm}-0.86620.49973-0.49973-0.86621.25948pt6.22845pt\text{pgfsys@invoke}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{pgfsys@setlinewidth}0.4pt\text{pgfsys@invoke}\text{pgfsys@moveto}6.45831pt3.22916pt\text{pgfsys@moveto}0.86113pt6.45831pt\text{pgfsys@stroke}\text{pgfsys@invoke}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{pgfsys@setlinewidth}0.4pt\text{pgfsys@invoke}\text{pgfsys@moveto}0.0pt0.0pt\text{pgfsys@curveto}-1.78343pt0.0pt-3.22916pt1.44572pt-3.22916pt3.22916pt\text{pgfsys@curveto}-3.22916pt4.58667pt-2.3801pt5.79929pt-1.10445pt6.26358pt\text{pgfsys@stroke}\text{pgfsys@invoke}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{pgfsys@setlinewidth}0.4pt\text{pgfsys@invoke}\text{pgfsys@moveto}0.0pt0.0pt\text{pgfsys@curveto}1.78343pt0.0pt3.22916pt1.44572pt3.22916pt3.22916pt\text{pgfsys@curveto}3.22916pt4.58667pt2.3801pt5.79929pt1.10445pt6.26358pt\text{pgfsys@stroke}\text{pgfsys@invoke}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{hss}\text{pgfsys@discardpath}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{hss}\text{lxSVG@closescope}\text{endpgfpicture}/\text{leavevmode}to10.09pt\text{vbox}to6.86pt\text{pgfpicture}\text{makeatletter}\text{lower}-0.2ptto0.0pt\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{definecolor}pgfstrokecolorrgb\mathrm{0,0,0}\text{pgfsys@color@rgb@stroke}000\text{pgfsys@invoke}\text{pgfsys@color@rgb@fill}000\text{pgfsys@invoke}\text{pgfsys@setlinewidth}0.4pt\text{pgfsys@invoke}\text{nullfont}to0.0pt\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{pgfsys@transformcm}0.86620.49973-0.499730.8662-1.25948pt6.22845pt\text{pgfsys@invoke}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{pgfsys@setlinewidth}0.4pt\text{pgfsys@invoke}\text{pgfsys@moveto}-6.45831pt3.22916pt\text{pgfsys@moveto}-0.86113pt6.45831pt\text{pgfsys@stroke}\text{pgfsys@invoke}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{pgfsys@setlinewidth}0.4pt\text{pgfsys@invoke}\text{pgfsys@moveto}0.0pt0.0pt\text{pgfsys@curveto}-1.78343pt0.0pt-3.22916pt1.44572pt-3.22916pt3.22916pt\text{pgfsys@curveto}-3.22916pt4.58667pt-2.3801pt5.79929pt-1.10445pt6.26358pt\text{pgfsys@stroke}\text{pgfsys@invoke}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{pgfsys@setlinewidth}0.4pt\text{pgfsys@invoke}\text{pgfsys@moveto}0.0pt0.0pt\text{pgfsys@curveto}1.78343pt0.0pt3.22916pt1.44572pt3.22916pt3.22916pt\text{pgfsys@curveto}3.22916pt4.58667pt2.3801pt5.79929pt1.10445pt6.26358pt\text{pgfsys@stroke}\text{pgfsys@invoke}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{hss}\text{pgfsys@discardpath}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{hss}\text{lxSVG@closescope}\text{endpgfpicture},& \text{for }{\phi }_{\text{ext}}\mathrm{mod}2\pi <\pi \\ -/+,& \text{for }{\phi }_{\text{ext}}\mathrm{mod}2\pi =\pi \\ \text{leavevmode}to10.09pt\text{vbox}to6.86pt\text{pgfpicture}\text{makeatletter}\text{lower}-0.2ptto0.0pt\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{definecolor}pgfstrokecolorrgb\mathrm{0,0,0}\text{pgfsys@color@rgb@stroke}000\text{pgfsys@invoke}\text{pgfsys@color@rgb@fill}000\text{pgfsys@invoke}\text{pgfsys@setlinewidth}0.4pt\text{pgfsys@invoke}\text{nullfont}to0.0pt\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{pgfsys@transformcm}0.86620.49973-0.499730.8662-1.25948pt6.22845pt\text{pgfsys@invoke}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{pgfsys@setlinewidth}0.4pt\text{pgfsys@invoke}\text{pgfsys@moveto}-6.45831pt3.22916pt\text{pgfsys@moveto}-0.86113pt6.45831pt\text{pgfsys@stroke}\text{pgfsys@invoke}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{pgfsys@setlinewidth}0.4pt\text{pgfsys@invoke}\text{pgfsys@moveto}0.0pt0.0pt\text{pgfsys@curveto}-1.78343pt0.0pt-3.22916pt1.44572pt-3.22916pt3.22916pt\text{pgfsys@curveto}-3.22916pt4.58667pt-2.3801pt5.79929pt-1.10445pt6.26358pt\text{pgfsys@stroke}\text{pgfsys@invoke}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{pgfsys@setlinewidth}0.4pt\text{pgfsys@invoke}\text{pgfsys@moveto}0.0pt0.0pt\text{pgfsys@curveto}1.78343pt0.0pt3.22916pt1.44572pt3.22916pt3.22916pt\text{pgfsys@curveto}3.22916pt4.58667pt2.3801pt5.79929pt1.10445pt6.26358pt\text{pgfsys@stroke}\text{pgfsys@invoke}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{hss}\text{pgfsys@discardpath}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{hss}\text{lxSVG@closescope}\text{endpgfpicture}/\text{leavevmode}to10.09pt\text{vbox}to6.86pt\text{pgfpicture}\text{makeatletter}\text{lower}-0.2ptto0.0pt\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{definecolor}pgfstrokecolorrgb\mathrm{0,0,0}\text{pgfsys@color@rgb@stroke}000\text{pgfsys@invoke}\text{pgfsys@color@rgb@fill}000\text{pgfsys@invoke}\text{pgfsys@setlinewidth}0.4pt\text{pgfsys@invoke}\text{nullfont}to0.0pt\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{pgfsys@transformcm}-0.86620.49973-0.49973-0.86621.25948pt6.22845pt\text{pgfsys@invoke}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{pgfsys@setlinewidth}0.4pt\text{pgfsys@invoke}\text{pgfsys@moveto}6.45831pt3.22916pt\text{pgfsys@moveto}0.86113pt6.45831pt\text{pgfsys@stroke}\text{pgfsys@invoke}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{pgfsys@setlinewidth}0.4pt\text{pgfsys@invoke}\text{pgfsys@moveto}0.0pt0.0pt\text{pgfsys@curveto}-1.78343pt0.0pt-3.22916pt1.44572pt-3.22916pt3.22916pt\text{pgfsys@curveto}-3.22916pt4.58667pt-2.3801pt5.79929pt-1.10445pt6.26358pt\text{pgfsys@stroke}\text{pgfsys@invoke}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{pgfsys@beginscope}\text{pgfsys@invoke}\text{pgfsys@setlinewidth}0.4pt\text{pgfsys@invoke}\text{pgfsys@moveto}0.0pt0.0pt\text{pgfsys@curveto}1.78343pt0.0pt3.22916pt1.44572pt3.22916pt3.22916pt\text{pgfsys@curveto}3.22916pt4.58667pt2.3801pt5.79929pt1.10445pt6.26358pt\text{pgfsys@stroke}\text{pgfsys@invoke}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{hss}\text{pgfsys@discardpath}\text{pgfsys@invoke}\text{lxSVG@closescope}\text{pgfsys@endscope}\text{hss}\text{lxSVG@closescope}\text{endpgfpicture},& \text{for }{\phi }_{\text{ext}}\mathrm{mod}2\pi >\pi \end{array},$$

$$H_{\text{int}}=\frac{C_{\text{g}}}{C_{\text{shunt}}+C_{\text{g}}}(2e\eta )V$$

$$H_{\text{int}}=\frac{L_{\text{s}}}{2L}\left({\varphi }_0\varphi \right)I$$

$$|{\mathcal{O}}_{\psi }{|}^2\equiv \frac{|⟨\psi |\mathcal{O}|0\text{∘}⟩{|}^2}{⟨0\text{∘}|{\mathcal{O}}^{†}\mathcal{O}|0\text{∘}⟩}$$

$$H^{\prime }=2{ϵ}_J{\delta }_J\sin \phi \sin \frac{\varphi }{2}.$$

$$H_{\text{eff}}^{\prime }=-\frac{16}{3\pi }{ϵ}_J{\delta }_J\left(\sin \phi -\frac{1}{5}\sin 3\phi \right),$$

$$H^{\prime }=-8{ϵ}_C\frac{{\delta }_C}{1-{\delta }_C^2}n(N-N_{\text{g}}-\eta ),$$

$$H^{\prime }={ϵ}_L\frac{{\delta }_L}{1-{\delta }_L^2}(\varphi -{\phi }_{\text{ext}})\theta .$$

$$\frac{1}{T_1}=\frac{1}{{\hslash }^2}|⟨0+|\mathcal{O}|0-⟩{|}^2\left[S_{\mathcal{EE}}(\Delta \omega )+S_{\mathcal{EE}}(-\Delta \omega )\right].$$

$$S_{VV}(\omega )+S_{VV}(-\omega )=\frac{2\hslash }{C_J{\mathcal{Q}}_{\text{cap}}(\omega )}\coth \frac{\hslash |\omega |}{2k_{\text{B}}T}.$$

$$S_{II}(\omega )+S_{II}(-\omega )=\frac{2\hslash }{L_i{\mathcal{Q}}_{\text{ind}}(\omega )}\coth \frac{\hslash |\omega |}{2k_{\text{B}}T}.$$

$$S_{\text{qp}}(\omega )+S_{\text{qp}}(-\omega )=2\hslash \omega \text{Re}{Y}_{\text{qp}}(\omega )\coth \frac{\hslash \omega }{2k_{\text{B}}T},$$

$$\frac{1}{T_{\text{upphi}}}=\frac{\pi }{(2\text{e}{)}^2}ϵ/\hslash ,$$

$$\frac{1}{T_{\upphi}}=A_{\varphi_{\text{ext}}}\bigg{|}\frac{\partial^{2}\Delta E}{\partial\varphi_{\text{ext}}^{2}}\bigg{|},$$

$$\frac{1}{T_{\text{upphi}}}=n_{\text{th}}\kappa \frac{{\chi }^2}{{\chi }^2+{\kappa }^2},$$

$$\frac{1}{T_{\text{upphi}}}=\sqrt{A_{{ϵ}_J}}\left|\frac{\partial \Delta E}{\partial {ϵ}_J}\right|,$$

$$4E_C{\left(i\frac{\partial }{\partial \phi }+N_{\text{g}}\right)}^2\psi -E_J\cos 2\phi \psi =E\psi .$$

$$g^{\prime \prime }+\left(\frac{E_J}{4E_C}\cos 2\phi +\frac{E}{4E_C}\right)g=0,$$

$$E_k(N_{\text{g}})=4E_C{a}_{N_{\text{g}}+\ell (k,N_{\text{g}})}\left(-\frac{E_J}{8E_C}\right),$$

$$E_k(N_{\text{g}})\approx E_k(1/2)-\frac{1}{2}{ϵ}_k\cos (\pi N_{\text{g}})$$

$${ϵ}_k\approx (-1{)}^k\frac{4^{k+2}}{k!}{E}_C\sqrt{\frac{2}{\pi }}{\left(\frac{2E_J}{E_C}\right)}^{\frac{2k+3}{4}}{\text{e}}^{-\sqrt{2E_J/E_C}},$$

$$\displaystyle\frac{\hbar^{2}}{8\epsilon_{C}}\biggl{[}2\ddot{\varphi}+\frac{1}{x}(\ddot{\varphi}+\ddot{\theta})\biggr{]}$$

Full text

I Introduction

I.1 Motivation

Superconducting circuits are widely recognized as a powerful potential platform for quantum computation and now stand at the frontier of quantum error correction Devoret and Schoelkopf (2013). Future progress will likely stem from two complementary strategies: (i) active error correction characterized by measurement-based Kelly et al. (2015); Ristè et al. (2015); Ofek et al. (2016) and autonomous Murch et al. (2012); Shankar et al. (2013); Leghtas et al. (2015); Puri et al. (2017) stabilization, and (ii) passive error correction characterized by protected qubits (Douçot and Ioffe, 2012, and references therein). We address strategy (ii) in this article by designing an experimentally accessible protected qubit.

The transmon qubit Koch et al. (2007) has proven to be a remarkably successful prototype of a protected qubit. Its circuit (depicted in Fig. 1a) contains a Josephson junction whose potential energy is $U=-E_{J}\cos\varphi$, with $E_{J}$ being the tunneling energy and $\varphi$ being the superconducting phase across the junction. The large shunt capacitance imposes a large effective mass for the analogous “particle in a box,” confining the low-energy wavefunctions near $\varphi=0$, the only minimum for $\varphi\in(-\pi,\pi)$. This confinement suppresses the susceptibility of the qubit to offset charge noise and renders the energy spectrum approximately harmonic with level spacing much smaller than $E_{J}$ (see Fig. 1b). On the other hand, circuit elements with degenerate phase states that only allow tunneling of pairs of Cooper pairs, meaning their potential energy is $U=-E_{J}\cos 2\varphi$, have been developed in recent years as a building block for topologically protected qubits Douçot and Vidal (2002); Gladchenko et al. (2009). In this article, we propose a transmon-like qubit with additional protection from environmental noise by combining the large shunt capacitance of the transmon with such a $\cos 2\varphi$ circuit element (the cross-hatched box in Fig. 1c). In this case, the wavefunctions are localized near $\varphi=0,\pi$ (see Fig. 1d), resulting in a nearly degenerate harmonic level arrangement. While the detrimental effects of offset charge noise are similarly suppressed in this circuit, sensitivity of the qubit to other decoherence mechanisms is also reduced, owing to the conservation of Cooper pair number parity.

In this article, we introduce a few-body transmon-type qubit where the charge carriers are exclusively pairs of Cooper pairs. Our central result is that there exists an experimentally attainable parameter regime for which conservative predictions of relaxation and dephasing times exceed $1\,$ms, i.e. an order of magnitude higher than those of typical transmons, given the same environmental noise Rigetti et al. (2012); Wang et al. (2019). In the remainder of Sec. I, we describe a toy model for the protected qubit. We proceed by analytically and numerically examining the Hamiltonian for the full superconducting circuit in Sec. II. Our attention in Sec. III then turns to properties of the ground state manifold, which we envision using as a protected qubit. A brief discussion about the concept of protection and examples of protected qubits, as well as our perspectives on readout and control, follows in Sec. IV. Finally, we summarize our results in Sec. V.

I.2 $\cos 2\varphi$ element

We first examine the advantages of the ideal circuit in Fig. 1c as a protected qubit. This circuit can be viewed as a Josephson-junction-like element (the cross-hatched box) shunted by a capacitance. Pairs of Cooper pairs are the only charge excitations permitted to tunnel through this element Douçot and Vidal (2002). In the Cooper pair number basis, the potential energy assumes the form

[TABLE]

where $E_{J}$ is the effective tunneling energy of the process. This expression follows from the conjugacy relation $[\varphi,N]=i$, where $N$ is the number of Cooper pairs that have tunneled. The invariance of the potential under translations in $\varphi$ by multiples of $\pi$ implies that half-fluxons are able to traverse the element.

The shunt capacitance and other charging effects produce a quadratic kinetic energy, yielding the Hamiltonian

[TABLE]

where $E_{C}$ is the charging energy and $N_{\text{g}}$ is the offset charge. This offset charge has been introduced due to the periodicity of the Hamiltonian in $\varphi$, which reflects the presence of a superconducting island in the circuit (as colored in Fig. 1c).

Since the circuit element only allows pairs of Cooper pairs to tunnel, the parity of the number of Cooper pairs that have tunneled is preserved under the action of the Hamiltonian. This property leads to two nearly degenerate ground states $|+\rangle$ and $|-\rangle$, which only consist of even and odd Cooper pair number states, respectively Gladchenko et al. (2009). Since these states have no overlap in charge space (equivalently, they have opposite periodicity in phase space—see Fig. 1d), we have $\langle-|\mathcal{O}|+\rangle\approx 0$ for any sufficiently local operator $\mathcal{O}$. Furthermore, the states 111This notation is chosen for consistency with Sec. II, where there is a correspondence to persistent current handedness. $|\leavevmode\hbox to13.32pt{\vbox to8.26pt{\pgfpicture\makeatletter\hbox{\hskip 6.65831pt\lower-1.59998pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ } {}{{}}{} {}{} {}{{}{}}{}{}{}{{}}{{}}{{}{}}{{}{}}{{{{}{}{{}} }}{{}}}{{{{}{}{{}} }}{{}}}{{{{}{}{{}} }}{{}}}{{{{}{}{{}} }}{{}}\pgfsys@beginscope\pgfsys@invoke{ } {\pgfsys@beginscope\pgfsys@invoke{ } {{}} {{{ {\pgfsys@beginscope\pgfsys@setlinewidth{0.32pt}\pgfsys@setdash{}{0.0pt}\pgfsys@roundcap\pgfsys@roundjoin{} {}{}{} {}{}{} \pgfsys@moveto{-1.19998pt}{1.59998pt}\pgfsys@curveto{-1.09998pt}{0.99998pt}{0.0pt}{0.09999pt}{0.29999pt}{0.0pt}\pgfsys@curveto{0.0pt}{-0.09999pt}{-1.09998pt}{-0.99998pt}{-1.19998pt}{-1.59998pt}\pgfsys@stroke\pgfsys@endscope}} }{{{}}{{{}}{\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@transformcm{0.8662}{0.49973}{-0.49973}{0.8662}{-1.25948pt}{6.22845pt}\pgfsys@invoke{ }\pgfsys@invoke{ \lxSVG@closescope }\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}}{{}}}} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}{{}{}}{{}{}}{{{{}{}{{}} }}{{}} {{{}}} } \pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }{}\pgfsys@moveto{-6.45831pt}{3.22916pt}\pgfsys@moveto{-0.86113pt}{6.45831pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{}{{}}{} {}{} {}{{}{}}{}{}{}{{}}{{}}{{}{}}{{}{}}{{{{}{}{{}} }}{{}}}{{{{}{}{{}} }}{{}}}{{{{}{}{{}} }}{{}}}{{{{}{}{{}} }}{{}}\pgfsys@beginscope\pgfsys@invoke{ } {\pgfsys@beginscope\pgfsys@invoke{ } {{}} {{}{{{}}{{{}}{\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@transformcm{-0.8662}{0.49973}{-0.49973}{-0.8662}{1.25948pt}{6.22845pt}\pgfsys@invoke{ }\pgfsys@invoke{ \lxSVG@closescope }\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}}{{}}}} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}{{}{}}{{}{}}{{{{}{}{{}} }}{{}} {{{}}} } \pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }{}\pgfsys@moveto{6.45831pt}{3.22916pt}\pgfsys@moveto{0.86113pt}{6.45831pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{}{{}}{}{}{{{{}{}{}{}}}{{}{}{}{}}}{} {} {}{}\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }{}\pgfsys@moveto{0.0pt}{0.0pt}\pgfsys@curveto{-1.78343pt}{0.0pt}{-3.22916pt}{1.44572pt}{-3.22916pt}{3.22916pt}\pgfsys@curveto{-3.22916pt}{4.58667pt}{-2.3801pt}{5.79929pt}{-1.10445pt}{6.26358pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{}{{}}{}{}{{{{}{}{}{}}} {{}{}{}{}}}{} {} {}{}\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }{}\pgfsys@moveto{0.0pt}{0.0pt}\pgfsys@curveto{1.78343pt}{0.0pt}{3.22916pt}{1.44572pt}{3.22916pt}{3.22916pt}\pgfsys@curveto{3.22916pt}{4.58667pt}{2.3801pt}{5.79929pt}{1.10445pt}{6.26358pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope\hss}}\lxSVG@closescope\endpgfpicture}}\rangle=\frac{1}{\sqrt{2}}(|+\rangle\pm|-\rangle)$ are respectively localized near $\varphi=0,\pi$ (see Fig. 1d). Because these states have suppressed overlap in phase space for large $E_{J}/E_{C}$ (i.e. they are roughly inversely periodic in charge space), we have $\langle\hskip 0.43057pt\leavevmode\hbox to10.09pt{\vbox to6.86pt{\pgfpicture\makeatletter\hbox{\hskip 3.42915pt\lower-0.2pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ } {}{{}}{} {}{} {}{{}{}}{}{}{}{{}}{{}}{{}{}}{{}{}}{{{{}{}{{}} }}{{}}}{{{{}{}{{}} }}{{}}}{{{{}{}{{}} }}{{}}}{{{{}{}{{}} }}{{}}\pgfsys@beginscope\pgfsys@invoke{ } {\pgfsys@beginscope\pgfsys@invoke{ } {{}} {{}{{{}}{{{}}{\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@transformcm{-0.8662}{0.49973}{-0.49973}{-0.8662}{1.25948pt}{6.22845pt}\pgfsys@invoke{ }\pgfsys@invoke{ \lxSVG@closescope }\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}}{{}}}} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}{{}{}}{{}{}}{{{{}{}{{}} }}{{}} {{{}}} } \pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }{}\pgfsys@moveto{6.45831pt}{3.22916pt}\pgfsys@moveto{0.86113pt}{6.45831pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{}{{}}{}{}{{{{}{}{}{}}}{{}{}{}{}}}{} {} {}{}\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }{}\pgfsys@moveto{0.0pt}{0.0pt}\pgfsys@curveto{-1.78343pt}{0.0pt}{-3.22916pt}{1.44572pt}{-3.22916pt}{3.22916pt}\pgfsys@curveto{-3.22916pt}{4.58667pt}{-2.3801pt}{5.79929pt}{-1.10445pt}{6.26358pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{}{{}}{}{}{{{{}{}{}{}}} {{}{}{}{}}}{} {} {}{}\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }{}\pgfsys@moveto{0.0pt}{0.0pt}\pgfsys@curveto{1.78343pt}{0.0pt}{3.22916pt}{1.44572pt}{3.22916pt}{3.22916pt}\pgfsys@curveto{3.22916pt}{4.58667pt}{2.3801pt}{5.79929pt}{1.10445pt}{6.26358pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope\hss}}\lxSVG@closescope\endpgfpicture}}|\mathcal{O}|\leavevmode\hbox to10.09pt{\vbox to6.86pt{\pgfpicture\makeatletter\hbox{\hskip 6.65831pt\lower-0.2pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ } {}{{}}{} {}{} {}{{}{}}{}{}{}{{}}{{}}{{}{}}{{}{}}{{{{}{}{{}} }}{{}}}{{{{}{}{{}} }}{{}}}{{{{}{}{{}} }}{{}}}{{{{}{}{{}} }}{{}}\pgfsys@beginscope\pgfsys@invoke{ } {\pgfsys@beginscope\pgfsys@invoke{ } {{}} {{}{{{}}{{{}}{\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@transformcm{0.8662}{0.49973}{-0.49973}{0.8662}{-1.25948pt}{6.22845pt}\pgfsys@invoke{ }\pgfsys@invoke{ \lxSVG@closescope }\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}}{{}}}} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}{{}{}}{{}{}}{{{{}{}{{}} }}{{}} {{{}}} } \pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }{}\pgfsys@moveto{-6.45831pt}{3.22916pt}\pgfsys@moveto{-0.86113pt}{6.45831pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{}{{}}{}{}{{{{}{}{}{}}}{{}{}{}{}}}{} {} {}{}\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }{}\pgfsys@moveto{0.0pt}{0.0pt}\pgfsys@curveto{-1.78343pt}{0.0pt}{-3.22916pt}{1.44572pt}{-3.22916pt}{3.22916pt}\pgfsys@curveto{-3.22916pt}{4.58667pt}{-2.3801pt}{5.79929pt}{-1.10445pt}{6.26358pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{}{{}}{}{}{{{{}{}{}{}}} {{}{}{}{}}}{} {} {}{}\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }{}\pgfsys@moveto{0.0pt}{0.0pt}\pgfsys@curveto{1.78343pt}{0.0pt}{3.22916pt}{1.44572pt}{3.22916pt}{3.22916pt}\pgfsys@curveto{3.22916pt}{4.58667pt}{2.3801pt}{5.79929pt}{1.10445pt}{6.26358pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope\hss}}\lxSVG@closescope\endpgfpicture}}\hskip 0.86108pt\rangle\approx 0$ for similarly local $\mathcal{O}$. These are precisely the conditions for simultaneously suppressing spurious transitions and phase changes between the states Ioffe and Feigel’man (2002) [resembling a Gottesman-Kitaev-Preskill (GKP) encoding on a circle Gottesman et al. (2001)].

More concretely, the ground state splitting obeys

[TABLE]

for large $E_{J}/E_{C}$ (see App. A). The two ground state energies oscillate out of phase with one another in $N_{\text{g}}$. Moreover, this shows that the splitting, as well as the charge dispersion, is exponentially suppressed in $E_{J}/E_{C}$. Thus, the role of the shunt capacitance is to decrease the charging energy $E_{C}$ and hence mitigate offset charge noise, much like in the transmon qubit Koch et al. (2007).

II Superconducting circuit

II.1 Hamiltonian

We now detail the superconducting circuit necessary to practically implement the sought-after $\cos 2\varphi$ Josephson element. This circuit, depicted in Fig. 2a, is composed of two identical arms, each containing a Josephson junction in series with superinductances Manucharyan et al. (2009); Masluk et al. (2012), arranged in parallel Kitaev (2006); Brooks et al. (2013) and shunted by a large capacitance. These superinductances are split in half and placed on either side of the respective Josephson junctions to avoid large capacitances to ground. Kirchhoff’s current law allows us to treat the phases across the superinductances in series as equal. When the external magnetic flux threading the inductive loop reaches half of a flux quantum, i.e. when $\varphi_{\text{ext}}=\pi$, the loop approximates the cross-hatched box in Fig. 1c. At this particular bias point, a Cooper pair can only tunnel through one of the Josephson junctions if it is accompanied by another Cooper pair tunneling through the other junction (in either direction). Conversely, a fluxon traversing a single Josephson junction corresponds to a half-fluxon traversing the whole element.

We consider the symmetric and antisymmetric combinations of superconducting phase coordinates

[TABLE]

and their conjugate charges $\{n,N,\eta\}$ in numbers of Cooper pairs. Note that the prefactor in the definition of $\varphi$ is chosen to bring the coordinate into agreement with the phase drop across the inductive loop in the limit that $\theta$ vanishes. The Hamiltonian reads

[TABLE]

where $\epsilon_{C}$ and $\epsilon_{J}$ are the single junction charging and tunneling energies, $2\epsilon_{L}$ is the inductive energy of each superinductance, $x\equiv C_{J}/C_{\text{shunt}}$ is the ratio of the junction capacitance to the shunt capacitance, and $N_{\text{g}}$ is the offset charge on the superconducting island (see Fig. 2a).

From this expression, it is clear that this circuit has three strongly coupled modes. The $\phi$ mode is flux dependent and is strongly and nonlinearly coupled, via the Josephson junctions, to the $\varphi$ mode. The $\varphi$ mode is offset charge dependent and strongly but linearly capacitively coupled to the $\theta$ mode. Our analysis and the effects observed in the remainder of this article require the parameter regime $\epsilon_{L}\ll\epsilon_{J}$, $\epsilon_{C}\lesssim\epsilon_{J}$, and $x\ll 1$ 222Specifically, the theory in Sec. II.2 breaks down at $\epsilon_{L}\sim\epsilon_{J}$ and the protection in Sec. III breaks down for $x\gtrsim 0.1$.. In particular, the parameters chosen for numerical simulations are listed in Tab. 1 and are similar to those of recent fluxonium devices Manucharyan et al. (2009); Earnest et al. (2018).

II.2 Semiclassical theory

In order to gain insight into the structure of the Hamiltonian in Eq. 2, we briefly revisit its potential energy $U$ in the $\varphi_{1}\varphi_{2}$-plane, which is plotted in Fig. 2b for $\theta=0$. The cosine terms in the potential form a two-dimensional “egg carton” of wells. The minimum of the quadratic term in the potential occurs at $\varphi_{1}+\varphi_{2}=\varphi_{\text{ext}}$, which generally falls between adjacent diagonal ridges of cosine wells. At the special value of $\varphi_{\text{ext}}=\pi$, these two ridges are degenerate. Near this value of the external flux, we consider the path of the system between neighboring potential minima by numerically solving for the three-dimensional instanton Matveev et al. (2002) trajectory (see App. B).

We then constrain the system to this maximally probable tunneling path. For the parameters mentioned above, this path is well-described by

[TABLE]

where $z\equiv\epsilon_{L}/\epsilon_{J}$ is a small parameter 333The main discrepancy between the approximate path in Eq. 3 and the calculated path in Fig. 2b is the sharp bending near the minima, which arises due to coupling to the $\theta$ variable.. Plugging this expression into the Hamiltonian in Eq. 2, approximating the resulting one-dimensional potential by the first few terms in its Fourier series, and Taylor expanding about $z=0$ yields $H\approx H_{\text{eff}}$ with

[TABLE]

to leading order. Here, $\phi_{\text{ext}}=\left|\varphi_{\text{ext}}-4\pi\,\text{round}\,\frac{\varphi_{\text{ext}}}{4\pi}\right|$ and we have discarded terms higher than the second harmonic or $O(\epsilon_{L})$ (see App. B for additional terms). This treatment exposes the “$\cos 2\varphi$ nature” of the potential at $\varphi_{\text{ext}}=\pi$, where the $\cos\varphi$ term vanishes. By comparison to Eq. 1, we see that the added complication is that the $\varphi$ mode is strongly coupled to the $\theta$ mode. The resulting hybridization is a central ingredient to understanding properties of the system beyond the ground state manifold (see Sec. II.4).

We comment that this approximation neglects quantum fluctuations that are perpendicular to the path in Eq. 3. This is consequently a semiclassical approximation: we have minimized the energy of the system with respect to the dynamical coordinate orthogonal to the trajectory Matveev et al. (2002). Moreover, from Fig. 2b, it is clear that the approximation we have made is that fluxons traverse a single Josephson junction at a time.

II.3 Energy spectrum

From numerical diagonalization of Eq. 2 (see App. C), we obtain the dependence of the energy levels on external flux as shown in Fig. 3a Smith et al. (2016). At $\varphi_{\text{ext}}=\pi$ (the dashed line in Fig. 3a), the spectrum resembles a doubled harmonic oscillator with energy $\sqrt{16x\epsilon_{L}\epsilon_{C}}$. Once $\varphi_{\text{ext}}$ deviates from $\pi$, half of the energy levels increase in energy linearly with slope $\sim\frac{32}{3}\epsilon_{L}$ 444This term is simply the potential energy at $\varphi=\pi$ relative to that at $\varphi=0$ (with $\phi_{\text{ext}}=0$) found from Eq. 4.. The other half of the energy levels form a flux-independent harmonic ladder.

We can understand this level structure as that of two emergent modes (see also App. D). The first mode is flux-dependent and its excitations correspond to the number of magnetic flux vortices, or fluxons, enclosed by the inductive loop in Fig. 2a. In turn, the number of enclosed fluxons identically maps onto the magnitude and chirality of the circulating persistent current in the inductive loop. The second mode is flux-independent and its excitations correspond to quantized charge density oscillations, or plasmons, across the inductive loop/shunt capacitance in Fig. 2a. Each plasmon involves the two superinductances (energy $2\epsilon_{L}$ in parallel) and the shunt capacitance (energy $x\epsilon_{C}$), and hence has energy $\sqrt{16x\epsilon_{L}\epsilon_{C}}$. Hereafter, we refer to these modes as the “fluxon mode” and the “plasmon mode,” respectively 555We comment that the Aharonov-Bohm effect manifests itself in this circuit through the fluxon mode, in which charge carriers circle a pinned flux. Conversely, the Aharonov-Casher effect involves the plasmon mode, in which flux vortices circle a pinned charge.. Additionally, we assign the labels |m\,\text{\raisebox{1.72218pt}{\scriptscriptstyle\bullet}}\,\rangle to the lowest-energy states, where $m$ denotes the number of plasmons and  $\scriptscriptstyle\bullet$  (or  $\scriptscriptstyle\circ$ ) denotes the presence (or absence) of a fluxon excitation relative to the ground state 666Here we restrict the Hilbert space to the two fluxon states with lowest energy.. Note that the fluxon placeholder indices are defined by

[TABLE]

where $|m\pm\rangle=|m\rangle\otimes\frac{1}{\sqrt{2}}(|\leavevmode\hbox to10.09pt{\vbox to6.86pt{\pgfpicture\makeatletter\hbox{\hskip 6.65831pt\lower-0.2pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ } {}{{}}{} {}{} {}{{}{}}{}{}{}{{}}{{}}{{}{}}{{}{}}{{{{}{}{{}} }}{{}}}{{{{}{}{{}} }}{{}}}{{{{}{}{{}} }}{{}}}{{{{}{}{{}} }}{{}}\pgfsys@beginscope\pgfsys@invoke{ } {\pgfsys@beginscope\pgfsys@invoke{ } {{}} {{}{{{}}{{{}}{\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@transformcm{0.8662}{0.49973}{-0.49973}{0.8662}{-1.25948pt}{6.22845pt}\pgfsys@invoke{ }\pgfsys@invoke{ \lxSVG@closescope }\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}}{{}}}} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}{{}{}}{{}{}}{{{{}{}{{}} }}{{}} {{{}}} } \pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }{}\pgfsys@moveto{-6.45831pt}{3.22916pt}\pgfsys@moveto{-0.86113pt}{6.45831pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{}{{}}{}{}{{{{}{}{}{}}}{{}{}{}{}}}{} {} {}{}\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }{}\pgfsys@moveto{0.0pt}{0.0pt}\pgfsys@curveto{-1.78343pt}{0.0pt}{-3.22916pt}{1.44572pt}{-3.22916pt}{3.22916pt}\pgfsys@curveto{-3.22916pt}{4.58667pt}{-2.3801pt}{5.79929pt}{-1.10445pt}{6.26358pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{}{{}}{}{}{{{{}{}{}{}}} {{}{}{}{}}}{} {} {}{}\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }{}\pgfsys@moveto{0.0pt}{0.0pt}\pgfsys@curveto{1.78343pt}{0.0pt}{3.22916pt}{1.44572pt}{3.22916pt}{3.22916pt}\pgfsys@curveto{3.22916pt}{4.58667pt}{2.3801pt}{5.79929pt}{1.10445pt}{6.26358pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope\hss}}\lxSVG@closescope\endpgfpicture}}\hskip 0.86108pt\rangle\pm|\hskip 0.43057pt\leavevmode\hbox to10.09pt{\vbox to6.86pt{\pgfpicture\makeatletter\hbox{\hskip 3.42915pt\lower-0.2pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ } {}{{}}{} {}{} {}{{}{}}{}{}{}{{}}{{}}{{}{}}{{}{}}{{{{}{}{{}} }}{{}}}{{{{}{}{{}} }}{{}}}{{{{}{}{{}} }}{{}}}{{{{}{}{{}} }}{{}}\pgfsys@beginscope\pgfsys@invoke{ } {\pgfsys@beginscope\pgfsys@invoke{ } {{}} {{}{{{}}{{{}}{\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@transformcm{-0.8662}{0.49973}{-0.49973}{-0.8662}{1.25948pt}{6.22845pt}\pgfsys@invoke{ }\pgfsys@invoke{ \lxSVG@closescope }\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}}{{}}}} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}{{}{}}{{}{}}{{{{}{}{{}} }}{{}} {{{}}} } \pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }{}\pgfsys@moveto{6.45831pt}{3.22916pt}\pgfsys@moveto{0.86113pt}{6.45831pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{}{{}}{}{}{{{{}{}{}{}}}{{}{}{}{}}}{} {} {}{}\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }{}\pgfsys@moveto{0.0pt}{0.0pt}\pgfsys@curveto{-1.78343pt}{0.0pt}{-3.22916pt}{1.44572pt}{-3.22916pt}{3.22916pt}\pgfsys@curveto{-3.22916pt}{4.58667pt}{-2.3801pt}{5.79929pt}{-1.10445pt}{6.26358pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{}{{}}{}{}{{{{}{}{}{}}} {{}{}{}{}}}{} {} {}{}\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }{}\pgfsys@moveto{0.0pt}{0.0pt}\pgfsys@curveto{1.78343pt}{0.0pt}{3.22916pt}{1.44572pt}{3.22916pt}{3.22916pt}\pgfsys@curveto{3.22916pt}{4.58667pt}{2.3801pt}{5.79929pt}{1.10445pt}{6.26358pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope\hss}}\lxSVG@closescope\endpgfpicture}}\rangle)$ (as in Sec. I.2) and the index represents the persistent current direction. This labeling serves the purpose of assigning quantum numbers consistently for all external flux values (see the colors in Fig. 3a 777Note that the coloring changes for only the odd-numbered plasmon states, due to the fact that the eigenstate with even overall parity has lower energy for each plasmon state, and that the overall parity includes both a plasmon and a fluxon component.), except the particular case $\varphi_{\text{ext}}\!\!\!\mod{2\pi}=0$ that we do not focus on.

II.4 Wavefunctions

Calculated charge wavefunctions $\langle N|\psi\rangle$ and phase wavefunctions $\langle\varphi,\phi|\psi\rangle$, obtained from numerical diagonalization of Eq. 2, are shown in Figs. 3b, 3c for the four lowest-energy eigenstates at $\varphi_{\text{ext}}=\pi$. Roughly, the phase wavefunctions are computed by projection of the $\theta$ coordinate and a Fourier transform to the $\varphi\phi$-plane, while the charge wavefunctions are computed by projection of the $\theta$ coordinate and constraint to the trajectory in Eq. 3 (see App. C for details).

The charge wavefunctions are grid states with Fock-state envelopes Gottesman et al. (2001); Bell et al. (2014). For fluxon excitation index $+/-$, these grid states are superpositions of even/odd Cooper pair number states 888The correspondence between Cooper pair number parity and the symmetry of the phase wavefunction can be seen from the $\pi$-periodicity/$\pi$-antiperiodicity of $|\!+\!/-\rangle$ and $\langle N|\psi\rangle=\frac{1}{2\pi}\int_{0}^{2\pi}\text{d}\varphi\,\text{e}^{iN\varphi}\langle\varphi|\psi\rangle$ Devoret (1997).. Additionally, $m$ corresponds to the order of the Fock state envelope. Note that a logical qubit encoded in $|0+\rangle$ and $|0-\rangle$ is protected from spurious transitions except those mediated by operators that flip Cooper pair parity.

On the other hand, the phase wavefunctions are approximately Fock states localized within the potential energy wells (see Sec. II.2) Earnest et al. (2018). The fluxon index $+/-$ denotes whether the state $|m\pm\rangle$ is a symmetric (bonding) or antisymmetric (antibonding) superposition of states localized within opposite ridges of potential wells. These ridges correspond to persistent currents of opposite chirality, and hence also to the absence/presence of a fluxon in the inductive loop of the circuit Dempster et al. (2014). In this picture, $m$ refers to the Fock order of the localized states. Finally, operators that flip Cooper pair parity correspond to odd functions of $\phi$ or functions of $\varphi$ with period an odd division of $2\pi$, which can be seen from Fig. 3c to mediate the transition $|m+\rangle\leftrightarrow|m-\rangle$

III Qubit

In Sec. II, we analyzed the multi-mode Hamiltonian describing the superconducting circuit in Fig. 2a. Numerical diagonalization of this Hamiltonian showed the emergence of a linear plasmon mode and a nonlinear fluxon mode 999The number of plasmon and fluxon modes form the rows and columns of a crude superconducting circuit periodic table. Charge-based qubits such as the transmon have one plasmon mode and no fluxon modes. Flux-based qubits such as the fluxonium have one fluxon mode and no plasmon modes. This circuit has two plasmon modes (one of which is neglected at low frequencies) and one fluxon mode. The larger resulting Hilbert space is the resource used to achieve protection.. In this section, we consider the properties of the logical qubit formed by $\{|0+\rangle,|0-\rangle\}$, the two lowest-energy eigenstates at $\varphi_{\text{ext}}=\pi$, which generalizes to \{|0\,\text{\raisebox{1.72218pt}{\scriptscriptstyle\circ}}\,\rangle,|0\,\text{\raisebox{1.72218pt}{\scriptscriptstyle\bullet}}\,\rangle\} away from $\varphi_{\text{ext}}=\pi$.

III.1 Matrix elements

To better elucidate which types of operators can and cannot induce transitions between the two states of the qubit, we examine the relevant matrix elements corresponding to capacitive and inductive coupling. This discussion is particularly relevant to understanding the expected dominant loss mechanisms (Sec. III.3) and designing a measurement and control apparatus that does not directly couple to the qubit (Sec. IV.3).

For capacitive coupling, a generic voltage $V$ couples to the superconducting island of the circuit in Fig. 2a via a gate capacitance $C_{\text{g}}$ and will append the term

[TABLE]

to the Hamiltonian in Eq. 2, in addition to dressing the shunt capacitance. This voltage may be a degree of freedom of another mode in the embedding circuit, a noise source, or an ac drive. We therefore see that the susceptibility of undergoing a transition from the ground state, due to capacitive coupling to the qubit island, is directly related to the matrix element \langle\psi|\eta|0\,\text{\raisebox{1.72218pt}{\scriptscriptstyle\circ}}\,\rangle.

For inductive coupling, a generic current $I$ couples to the circuit via a small inductance $L_{\text{s}}$ shared with the inductive loop, which adds the term

[TABLE]

to the Hamiltonian in Eq. 2. Here, $\phi_{0}=\hbar/2e$ is the reduced magnetic flux quantum and $L$ is the superinductance in each arm of the qubit (i.e. $\epsilon_{L}=\phi_{0}^{2}/L$). Like the voltage source, this current may represent an internal or environmental degree of freedom. We see that the susceptibility of undergoing a transition from the ground state, due to inductive coupling to the inductive loop, is related to the matrix element \langle\psi|\phi|0\,\text{\raisebox{1.72218pt}{\scriptscriptstyle\circ}}\,\rangle.

Limiting the Hilbert space to the six lowest-energy eigenstates $\{|m\pm\rangle:m=0,1,2\}$, we numerically compute the normalized matrix elements

[TABLE]

from the ground state |0\,\text{\raisebox{1.72218pt}{\scriptscriptstyle\circ}}\,\rangle for the operators $\mathcal{O}=\eta,\phi$. Results are plotted in Fig. 4. Note that $\sum_{\psi}|\mathcal{O}_{\psi}|^{2}=1$ and $|\mathcal{O}_{\psi}|^{2}>0$, so we may reasonably consider these as transition probabilities via $\mathcal{O}$.

We see from Fig. 4a that transitions mediated by capacitive coupling to the qubit island are only allowed from |0\,\text{\raisebox{1.72218pt}{\scriptscriptstyle\circ}}\,\rangle to |1\,\text{\raisebox{1.72218pt}{\scriptscriptstyle\circ}}\,\rangle. These selection rules result from the decoupling of the even and odd Cooper pair number parity manifolds (see App. D). Most importantly, transitions between qubit states are forbidden, meaning capacitive coupling offers a promising ingredient for indirect qubit measurement and control (see Sec. IV.3). Conversely, inductive coupling to the inductive loop of the qubit permits transitions between |0\,\text{\raisebox{1.72218pt}{\scriptscriptstyle\circ}}\,\rangle and |0\,\text{\raisebox{1.72218pt}{\scriptscriptstyle\bullet}}\,\rangle in the vicinity of $\varphi_{\text{ext}}=\pi$, as shown in Fig. 4b. This effect arises because the operator $\phi$ induces transitions between the Cooper pair parity manifolds, as can be seen from the Fourier series for Eq. 3. As a consequence, we expect that relaxation of the qubit will be primarily due to inductive loss in the superinductances (see Sec. III.3).

III.2 Disorder

A highly symmetric superconducting circuit is usually fragile in view of unavoidable fabrication imperfections Dempster et al. (2014). The symmetry of our circuit involving the two inductive arms in Fig. 2a may be broken in three parameters: the Josephson energies of the junctions, the capacitances of the junctions, or the superinductances. To analyze these effects, we numerically diagonalize Eq. 2 and examine the energy splitting $\Delta E$ 101010The energy splitting $\Delta E$ is evaluated at $N_{\text{g}}=0$. as well as the charge dispersion $\epsilon=\max_{N_{\text{g}}}\Delta E-\min_{N_{\text{g}}}\Delta E$ of the $\{|0+\rangle,|0-\rangle\}$ manifold at $\varphi_{\text{ext}}=\pi$. A dimensionless quantity $\delta\in[0,1)$ is introduced to parameterize the extent of asymmetry in all three cases, and the $\delta$-dependence of the energies $\Delta E$ and $\epsilon$ is studied.

III.2.1 Disorder in $\epsilon_{J}$

We model disorder in the Josephson energies of the junctions by allowing the values of $\epsilon_{J}$ to deviate. We therefore set the left and right junction tunneling energies to $(1\pm\delta_{J})\epsilon_{J}$, respectively, where $\delta_{J}$ is the aforementioned asymmetry parameter. The Hamiltonian in Eq. 2 is perturbed by the term

[TABLE]

See Fig. 5 for a plot of $\Delta E$ and $\epsilon$ as a function of $\delta_{J}$. The important feature in these plots is that the charge dispersion decreases exponentially while the splitting increases exponentially with $\delta_{J}$ 111111The extrapolation of the charge dispersion as a function of the disorder beyond $\delta=0.6$ is due to numerical instabilities inherent to low-energy quantities paired with disappearance of an efficient diagonalization basis for large disorder. This is particularly obvious in the vicinity of $\delta=1$, at which point the number of modes of the circuit reduces from three to two.. These features arise from the effective Hamiltonian in Eq. 4 being accompanied by $2\pi$-periodic terms in the presence of disorder. In this case of disorder in $\epsilon_{J}$, the approximations in Sec. II.2 lead to $H^{\prime}\approx H_{\text{eff}}^{\prime}$ with

[TABLE]

which evidently permits the tunneling of single Cooper pairs across the element. The resulting qubit retains characteristics of the symmetric circuit as well as the asymmetric/transmon-like circuit.

III.2.2 Disorder in $\epsilon_{C}$

Analogous to Sec. III.2.1, we set the left and right junction charging energies to $\epsilon_{C}/(1\pm\delta_{C})$, respectively. Aside from dressing the charging energy, the Hamiltonian in Eq. 2 inherits the term

[TABLE]

whose effect on the qubit manifold is plotted in Fig. 5.

This form for capacitive disorder is assumed due to the fact that, if $\delta_{J}=\delta_{C}\equiv\delta_{A}$, then the product $\epsilon_{J}\epsilon_{C}$ is kept constant for both junctions under the effects of disorder. This corresponds to the physical case where the junction plasma frequencies are fixed by oxidation, but their areas differ due to fabrication imperfections. Here, the junction areas are $(1\pm\delta_{A})A$ because $A\propto\sqrt{\epsilon_{J}/\epsilon_{C}}$. The consequences of area disorder are plotted in Fig. 5.

III.2.3 Disorder in $\epsilon_{L}$

Following the same procedure, we set the left and right superinductance inductive energies to $\epsilon_{L}/(1\pm\delta_{L})$, respectively. This form is taken in order to fix the total linear inductance in the loop. Aside from dressing the inductive energy, the Hamiltonian in Eq. 2 is perturbed by

[TABLE]

Note in Fig. 5 that the charge dispersion and energy splitting follow the same general trend for inductive disorder as for the other three. The key difference is that the charge dispersion decreases more quickly than for any other form of disorder. Oppositely, the splitting is initially the same as for area disorder, but the slope decreases in $\delta_{L}$. We conclude that disorder allows us to engineer a circuit with a sufficiently non-degenerate ground state manifold whose charge dispersion is largely suppressed. For reasons that will become clear in Sec. III.4, these features are extremely valuable for designing a qubit that is protected from dephasing.

III.3 Relaxation

We model loss due to an arbitrary channel using Fermi’s Golden Rule. This gives us the relaxation rate, including both emission and absorption, of

[TABLE]

In this expression, $\mathcal{O}$ is the operator of the circuit coupling to a noisy bath variable $\mathcal{E}(t)$, the noise spectral density of which is given by $S_{\mathcal{EE}}(\omega)$. Note that $\Delta E=\hbar\,\Delta\omega$. In this subsection, we calculate the relaxation rates for the expected dominant loss mechanisms of the qubit based on numerical diagonalization of the full Hamiltonian in Eq. 2. We perform this calculation for various degrees of inductive disorder $\delta_{L}$ (see Fig. 5).

We consider four possible loss channels for the qubit: capacitive loss, inductive loss, Purcell loss, and quasiparticle tunneling Pop et al. (2014). Capacitive loss involves dielectric dissipation in the Josephson junction capacitances. In this case, we have \mathcal{O}=2eN_{i}=2e\big{[}n\pm\frac{1}{2}(N-\eta)\big{]} for the charge across the $i$-th junction and $\mathcal{E}=V$ for a bath voltage with

[TABLE]

Here, $C_{J}=e^{2}/2\epsilon_{C}$ is the junction capacitance, $T$ is the temperature, and $\mathcal{Q}_{\text{cap}}(\omega)$ is a frequency-dependent quality factor Braginsky et al. (1987) with nominal value $Q_{\text{cap}}\sim 1\times 10^{6}$ Wang et al. (2015) (see App. E for more details).

Depending on the specific implementation, inductive loss may occur within the superinductances via quasiparticle tunneling Pop et al. (2014). This situation can be modeled by taking $\mathcal{O}=\phi_{0}\phi_{i}$ for the flux across the $i$-th superinductance and $\mathcal{E}=I$ for a bath current with

[TABLE]

In this expression, $L_{i}=(1\pm\delta_{L})L$ is the $i$-th superinductance and $\mathcal{Q}_{\text{ind}}(\omega)$ is a frequency-dependent quality factor with nominal value $Q_{\text{ind}}\sim 500\times 10^{6}$ Pop et al. (2014). As shown in App. E, this frequency dependence is included to extrapolate to small qubit transition frequencies.

Loss due to the Purcell effect should mainly arise from coupling of the qubit to the plasmon mode, which we model as dielectric loss in the shunt capacitance. As such, we have $\mathcal{O}=2e\eta$ and $\mathcal{E}=V$ for a bath voltage with the same noise spectral density as Eq. 7, but with $C_{J}$ replaced by $C_{\text{shunt}}$.

Finally, quasiparticle tunneling is expected to occur across either Josephson junction and contribute to qubit relaxation Pop et al. (2014). In this case, we have $\mathcal{O}=2\phi_{0}\sin\frac{\varphi_{i}}{2}$ for the $i$-th junction and the noise spectral density

[TABLE]

with $\text{Re}\,Y_{\text{qp}}(\omega)$ being the dissipative part of the Josephson junction admittance Catelani et al. (2011) (see App. E for the explicit form). Note that this admittance scales linearly in both $\epsilon_{J}$ and the quasiparticle density $x_{\text{qp}}$ defined relative to that of Cooper pairs.

The calculated relaxation rates and the corresponding components of Eq. 6 are shown in Tab. 2 for all four loss channels at $\varphi_{\text{ext}}=\pi$. For both capacitive and inductive loss, we note that measurements have not yet been performed using qubits with extremely small energy splittings, and the calculation presented here relies on an extrapolation to such frequencies. A key feature is the complete absence of quasiparticle loss, as in the fluxonium qubit Pop et al. (2014). Additionally, we see that the asymmetric qubit is marginally less susceptible to inductive loss than the symmetric qubit. This improvement comes at the cost of the susceptibility to capacitive and Purcell loss. However, we emphasize that the lifetimes shown in Tab. 2 are conservative estimates that demonstrate $T_{1}\gtrsim 1$ ms, which is at least competitive with state-of-the-art qubit implementations Rigetti et al. (2012); Barends et al. (2013); Pop et al. (2014); Bell et al. (2014); Yan et al. (2016); Groszkowski et al. (2018).

III.4 Pure dephasing

In this subsection, we examine the dependence of the qubit transition energy $\Delta E$ on various system parameters $\lambda$ corresponding to different dephasing mechanisms. We estimate the dephasing times due to offset charge noise, external flux noise, photon shot noise in the plasmon mode, and critical current noise. The noise spectral densities are all assumed to be $1/f$ and hence given by $S_{\lambda\lambda}(\omega)=2\pi A_{\lambda}/|\omega|$, where $\sqrt{A_{\lambda}}$ is the amplitude, with the sole exception of shot noise whose spectral density is assumed to be Lorentzian.

III.4.1 Charge noise

In the case of perfect symmetry, the charge dispersion $\epsilon$ is identically mapped to $\Delta E$. As a consequence, resilience to dephasing from offset charge noise demands a high degree of degeneracy, making experimental implementation difficult. Fortunately, this issue can be avoided by introducing inductive disorder into the system (see Sec. III.2.3). In the slow-varying limit for charge noise, the pure dephasing rate is bounded by

[TABLE]

where e is Euler’s number Koch et al. (2007). In Tab. 2, we see that this yields a strict bound on the decoherence time in the case of perfect symmetry, which is greatly alleviated in the presence of inductive disorder. Note that this estimate does not explicitly involve $\sqrt{A_{N_{\text{g}}}}$ because the slow-varying limit has been taken where the offset charge assumes a random value for each measurement, but does not fluctuate within a given measurement.

III.4.2 Flux noise

At $\varphi_{\text{ext}}=\pi$, the qubit manifold is only susceptible to external flux noise to second order. However, one can easily obtain 121212For example, this can done using the double-well Hamiltonian $H=\frac{1}{2}(\Delta E)\sigma_{x}+\frac{1}{2}\pi\epsilon_{L}(\varphi_{\text{ext}}-\pi)\sigma_{z}$ for $\varphi_{\text{ext}}$ in the vicinity of $\pi$, where the two spin states correspond to $\phi=0,\pi$. One then has the flux-dependent splitting $\delta E=\sqrt{(\Delta E)^{2}+\pi^{2}\epsilon_{L}^{2}(\varphi_{\text{ext}}-\pi)^{2}}$ and hence $\partial^{2}(\delta E)/\partial\varphi_{\text{ext}}^{2}=(\pi\epsilon_{L})^{2}/\Delta E$ (at $\varphi_{\text{ext}}=\pi$). the relation $\partial^{2}\Delta E/\partial\varphi_{\text{ext}}^{2}\propto 1/\Delta E$ (at $\varphi_{\text{ext}}=\pi$), which shows that insensitivity to flux-noise-induced dephasing requires sufficiently weak degeneracy, at odds with the requirement for resilience to charge noise. Once more, introducing inductive disorder avoids this issue. Neglecting components that depend on details of the implementation, the pure dephasing rate Cottet (2002) is bounded by

[TABLE]

where $\sqrt{A_{\varphi_{\text{ext}}}}/2\pi\sim 3\times 10^{-6}$ is the amplitude of the noise spectral density Quintana et al. (2017); Kou et al. (2017) for $\varphi_{\text{ext}}$. As in the case of charge noise, this places a stringent limit on the decoherence time of the qubit in the absence of inductive disorder, as we see in Tab. 2.

III.4.3 Shot noise

We expect the dominant contribution to dephasing due to thermal photons to arise from coupling to the plasmon mode. This mode has an angular frequency $\omega_{\text{p}}\approx\sqrt{16x\epsilon_{C}\epsilon_{L}}/\hbar$ and an occupation $n_{\text{p}}$ with mean given by $n_{\text{th}}=1/(\text{e}^{\hbar\omega_{\text{p}}/k_{\text{B}}T}-1)$. To this end, we place a bound on the pure dephasing time $T_{\upphi}$ using the expression Rigetti et al. (2012)

[TABLE]

where $\chi$ is the dispersive shift of the qubit on the plasmon mode and $\kappa=\omega_{\text{p}}/\mathcal{Q}_{\text{cap}}(\omega_{\text{p}})$ is the linewidth of the plasmon mode. Note that, as in Sec. III.3, we use the frequency-dependent dielectric quality factor $\mathcal{Q}_{\text{cap}}(\omega)$ with nominal value $Q_{\text{cap}}\sim 1\times 10^{6}$ Wang et al. (2015). In all cases, this form of noise is not expected to limit the decoherence time of the qubit (see Tab. 2).

III.4.4 Critical current noise

The final dephasing mechanism we investigate is critical current noise, or fluctuations in the Josephson energy $\epsilon_{J}$. As in the case of flux noise, we discard factors that are sensitive to experimental details and take the dephasing rate to be bounded by Cottet (2002)

[TABLE]

where $\sqrt{A_{\epsilon_{J}}}\sim 5\times 10^{-7}\epsilon_{J}$ is the amplitude of the noise spectral density van Harlingen et al. (2004) for $\epsilon_{J}$. As in the case of photon shot noise, we do not expect critical current noise to place a strict bound on the decoherence time (see Tab. 2).

At this point, it is clear that inductive asymmetry constitutes a necessary ingredient in the protection of this qubit. In light of Eq. 5, we attribute this to the resulting hybridization of the $\phi$ and $\theta$ modes, which are not directly coupled in the symmetric case (see Eq. 2). The fluxon transition between the qubit states inherits some character of the plasmon transition, thereby breaking the correspondence between $\epsilon$ and $\Delta E$ and reducing the flux matrix element $\langle 0\!-\!|\phi|0+\rangle$.

IV Discussion

IV.1 Protection

We have proposed a superconducting circuit whose ground state manifold can be made protected, at the Hamiltonian level, from multiple common sources of noise. Here we describe some of the other strategies for achieving protection. The simplest manifestations involve improving quality factors associated with coupling to different thermal baths Martinis et al. (2005); Houck et al. (2008); Barends et al. (2013) and reducing noise spectral densities for different Hamiltonian parameters Paik et al. (2011); Rigetti et al. (2012); Kumar et al. (2016). At a higher level, a popular approach to suppressing qubit relaxation has been to localize wavefunctions in disparate regions of phase space to lessen transition matrix elements Pop et al. (2014); Bell et al. (2014); Earnest et al. (2018); Lin et al. (2018). On the other hand, delocalization of the same wavefunctions has been shown to mitigate dephasing effects by reducing qubit sensitivity to Hamiltonian parameters van der Wal et al. (2000); Vion et al. (2002); Koch et al. (2007); Manucharyan et al. (2009); Steffen et al. (2010); Yan et al. (2016). Superconducting circuits with multiple degrees of freedom whose qubit wavefunctions are both localized and delocalized combine both of these approaches. In these circuits, quantum information is diffused among constituent local degrees of freedom, providing protection from local perturbations. The many-body limit promises topological protection Douçot and Vidal (2002); Ioffe and Feigel’man (2002); Douçot et al. (2003); Kitaev (2006), in which global operators are necessary to manipulate logical qubits, but the few-body case described here offers an experimentally realistic approximation Cottet et al. (2002); Gladchenko et al. (2009); Brooks et al. (2013); Bell et al. (2014); Dempster et al. (2014); Kou et al. (2017); Groszkowski et al. (2018). In the following subsection, we outline the key differences between our circuit and similar proposals.

IV.2 Comparison to other protected qubits

IV.2.1 Heavy fluxonium

The heavy fluxonium Earnest et al. (2018); Lin et al. (2018), in which a conventional fluxonium qubit is shunted by a large capacitance, shares many similar features with the circuit in Fig. 2a. In fact, our circuit roughly reduces to that of the conventional fluxonium in the extremely asymmetric limit of $\delta_{L}\rightarrow 1$, and the eigenstates of both circuits are similarly represented in terms of persistent currents. The two main differences are that our circuit contains an additional small Josephson junction and that the shunt capacitance is placed across the junction along with a subset of the array junctions. An analogy can be drawn between $\delta_{L}$ in our circuit and the shunt capacitance in the heavy fluxonium, because they both suppress the qubit transition frequency. Physically, the qubit eigenstates differ when these analogous parameters fall to small values. While the eigenstates collapse to Cooper pair number parity states in our circuit, this is not true for the heavy fluxonium because the variable $\varphi$ is not compact.

IV.2.2 Rhombus qubit

The circuit in Fig. 2a bears resemblance to the rhombus Blatter et al. (2001), with two central differences. First, in the rhombus, the superinductances are replaced by single Josephson junctions, or arrays of a few larger junctions Douçot and Ioffe (2012); Bell et al. (2014). This changes the parameter regime of the circuit from $\epsilon_{L}\ll\epsilon_{J}$ to $\epsilon_{L}\sim\epsilon_{J}$, decreasing the amplitude of the $\cos 2\varphi$ term in the Hamiltonian. Second, the shunt capacitance is replaced by a gate capacitance to a voltage source Bell et al. (2014). This is akin to substituting an electrostatic gate for a shunt capacitance to obtain the Cooper pair box Vion et al. (2002) from the transmon Koch et al. (2007); overall suppression of the charge dispersion is traded for the ability to bias the circuit at its charge sweet spot.

Finally, the rhombus by itself is not designed to be a protected qubit. Rather, when multiple rhombi are arranged into a one-dimensional chain (or a two-dimensional fabric), the ground states are eigenstates of a nonlocal operator, which provides topological protection Douçot and Vidal (2002); Ioffe and Feigel’man (2002); Douçot et al. (2003); Douçot and Ioffe (2012). On one hand, our qubit does not require such scaling to achieve protection. On the other hand, the protection we predict is inherently susceptible to local perturbations (see Secs. III.3, III.4).

IV.2.3 0-$\pi$ qubit

The 0-$\pi$ qubit Brooks et al. (2013) also has a similar superconducting circuit to that in Fig. 2a, but with three essential distinctions. First, the pairs of superinductances on each arm of the circuit are combined, altering the embedding capacitance matrix, in particular the capacitances to ground. Second, there is the addition of a second large capacitance shunting the inductive loop between its two horizontally oriented nodes. When this second capacitance is precisely $C_{\text{shunt}}$, this permits the exact decoupling of the $\theta$ mode from the $\varphi$ mode in Eq. 2. Third, the 0-$\pi$ qubit is operated in a parameter regime where $\epsilon_{L}\sim 0.01\epsilon_{J}$, as opposed to our circuit, where $\epsilon_{L}\sim 0.1\epsilon_{J}$ Dempster et al. (2014). This additional order of magnitude in the superinductance may prove marginally less accessible experimentally Manucharyan et al. (2009); Masluk et al. (2012).

Notably, the inductive loop in the 0-$\pi$ qubit is threaded with $\varphi_{\text{ext}}\approx 0$ at its working point instead of $\varphi_{\text{ext}}=\pi$. This leads to a substantial change in the physics of the ground state manifold; $|0+\rangle$ and $|0-\rangle$ are approximately localized in distinct potential wells Dempster et al. (2014). In our case, these states are approximately the symmetric and antisymmetric superpositions of the localized wavefunctions, which are themselves localized in distinct Cooper pair number parities.

IV.3 Readout and control

Protected qubits face the serious obstacle of realizing state manipulation and measurement while remaining sufficiently isolated from their environments to preserve their coherence. We envision performing readout and control using an ancillary mode structure, which enables cascaded dispersive readout and Raman indirect transitions, as outlined below.

For readout, we aim to exploit the sizable native dispersive shift of the qubit on the plasmon mode, $\chi/2\pi\sim-20\,\text{MHz}$ (see Sec. III.4.3). Unfortunately, the small anharmonicity, which is at most of order $10\,\text{MHz}$, of the plasmon mode makes readout of the plasmon mode using a linear mode (e.g. a microwave cavity) difficult. As a remedy, we propose introducing an ancillary anharmonic mode by which to measure the plasmon state. In this scheme, a dispersive interaction will be mediated by the ancillary mode between the readout and plasmon modes, thereby enabling dispersive measurement with two readout tones Kirchmair et al. (2013).

The above-mentioned ancillary anharmonic mode will also be useful for control of the protected qubit. As shown in Sec. III.1, direct transitions mediated by capacitive coupling to the qubit superconducting island are strictly forbidden in the symmetric case, and weakly forbidden in the asymmetric case. Note that this coupling method is chosen to minimize the contributions to decoherence resulting from the introduction of the readout/control circuit. For manipulation, we therefore propose transitioning through the $|1+\rangle$ or $|1-\rangle$ state, as in a conventional $\Lambda$ system. An additional complication is that, even in the presence of inductive asymmetry, the joint fluxon-plasmon transitions $|0+\rangle\leftrightarrow|1-\rangle$ and $|0-\rangle\leftrightarrow|1+\rangle$ are forbidden.

Using an ancillary mode based on the Superconducting Nonlinear Asymmetric Inductive eLement (SNAIL) Frattini et al. (2017), it is possible to engineer selection rules where the joint fluxon-plasmon transitions are weakly permitted while the qubit transition is left unaffected, as has recently been demonstrated in a fluxonium qubit Vool et al. (2018). This stems from the feature of the SNAIL that allows, at specific flux biases, third-order nonlinearities without fourth-order nonlinearities. If we parameterize the broken selection rule using $y\in[0,1)$, the matrix element of $|0+\rangle\leftrightarrow|1-\rangle$ relative to $|0+\rangle\leftrightarrow|1+\rangle$, then the gate speed via sequential direct transitions and a stimulated Raman transition is estimated to be $(1+y^{-1})t_{\mathrm{gate}}$ and $2\Delta y^{-1}t_{\mathrm{gate}}^{2}$, respectively Steck (2019). Here, $\Delta\sim 100\,\text{MHz}$ is the detuning of the Raman drive from the intermediate state and $t_{\mathrm{gate}}\sim 10\,\text{ns}$ is the $\pi$-pulse time Chow et al. (2010). This also introduces an additional Purcell loss channel $|0-\rangle\rightarrow|1+\rangle\rightarrow|0+\rangle$ (and vice versa), which corresponds to the relaxation time $T_{1}\sim y^{-2}\mathcal{Q}_{\text{cap}}(\omega_{\mathrm{p}})/\omega_{\mathrm{p}}$. We therefore expect gate speeds to be slowed by a factor of $y^{-1}<10$ while lifetimes are kept in the millisecond range.

V Conclusion

In summary, we have designed a few-body superconducting circuit in which the charge carriers are well-approximated by pairs of Cooper pairs at a particular bias point. The Josephson tunneling element that supports these charge carriers is characterized by a $\cos 2\varphi$ term in the Hamiltonian, whose emergence we have shown analytically. Our numerical simulations supplement these arguments and demonstrate protection against a variety of common relaxation and dephasing sources. We find that this protection is substantially enhanced in the presence of disorder. Finally, we compared our circuit to similar proposals and offered our perspectives on readout and control.

As a final remark, we comment that engineering a circuit whose potential energy is dominated by a $\cos 2\varphi$ term opens the door to more exotic designs with potential energies of the form $\cos\mu\varphi$, with $\mu\in\mathbb{N}$, which could be obtained by introducing additional loops in the circuit. These could be tremendously valuable for quantum simulation, realizing nearly degenerate ground state manifolds with greater multiplicities, or performing degeneracy-preserving measurements of photon number parity Cohen et al. (2017).

Acknowledgements.

We thank J. Cohen, S. M. Girvin, L. I. Glazman, M. Mirrahimi, R. J. Schoelkopf, K. Serniak, and S. Shankar for their insights. This work was supported by the Army Research Office Grant No. W911NF-14-1-0011 and W911NF-18-1-0212. W.C.S. was supported by DoDMURI award No. FP057123-C. Facilities use was supported by the Yale Institute of Nanoscience and Quantum Engineering under National Science Foundation Grant No. MRS1119826.

Appendix A Mathieu equation

The time-independent Schrödinger equation for Eq. 1 assumes the form

[TABLE]

Using the trial function $g\equiv\text{e}^{-iN_{\text{g}}\varphi}\psi$, this equation becomes

[TABLE]

which is the Mathieu equation Meixner et al. (1980); Koch et al. (2007). The $k$-th eigenenergy is given by

[TABLE]

where $a$ denotes the Mathieu characteristic value and $\ell$ is an integer that sorts the eigenvalues. The crucial difference between this solution and that for the Cooper pair box is that the characteristic exponent $2(N_{\text{g}}+\ell)$ is replaced by $N_{\text{g}}+\ell$ Koch et al. (2007). This gives rise to the nearly degenerate ground states in the Hamiltonian in Eq. 1. In particular, for $E_{J}\gg E_{C}$, we can use the asymptotic form of the Mathieu characteristic values to find

[TABLE]

with $\epsilon_{k}$ being the charge dispersion of the $k$-th level. Importantly, for even values of $k$, we have $E_{k}=E_{k+1}\approx(\sqrt{32E_{J}E_{C}}-2E_{C})k$ (neglecting constants and nonlinear terms) and $\epsilon_{k}=-\epsilon_{k+1}$. Additionally, the expression for the ground state splitting in Sec. I.2 is justified by

[TABLE]

and the fact that $\Delta E=E_{1}-E_{0}=\epsilon_{0}\cos(\pi N_{\text{g}})$.

Appendix B Instanton treatment

The average tunneling trajectory between potential minima for the Hamiltonian in Eq. 2 can be found by examining an equivalent classical problem. This is done by inverting the potential and solving the classical equations of motion for a particle that starts at one maxima and ends at the other, each with zero kinetic energy Matveev et al. (2002). Consequently, this classical trajectory requires an infinite amount of time. The corresponding differential equations,

[TABLE]

are solved numerically. Carrying out the procedure in Sec. II.2, but retaining terms through the fourth harmonic and $O(z^{2})$, results in the extended version of the effective Hamiltonian

[TABLE]

which clearly reduces to Eq. 4. At $\varphi_{\text{ext}}=\pi$, we see that the dominant correction term is of the form $\cos 4\varphi$. Moreover, we observe that the odd Fourier terms exactly vanish for the symmetric circuit at this flux bias, while the coefficients for the even Fourier terms decay roughly in powers of $z$.

Appendix C Numerical diagonalization

For numerical diagonalization of either Eq. 2 or Eq. 4, we employ a charge basis for the $\varphi$ mode in order to efficiently capture the dependence of the Hamiltonian on $N_{\text{g}}$ Devoret (1997). This is due to the exact periodicity of the Hamiltonian in $\varphi$, which inhibits the use of a harmonic oscillator basis to capture the multi-well physics apparent in Fig. 2b. For the $\theta$ and $\phi$ modes, we use harmonic oscillator bases. Explicitly, the full Hamiltonian in Eq. 2 can be written as

[TABLE]

after discarding constants, where $a^{\dagger}/a$ and $b^{\dagger}/b$ are bosonic creation/annihilation operators for the $\phi$ and $\theta$ modes, respectively. Additionally, the zero point fluctuation amplitudes \phi_{\text{zpf}}=\big{(}\frac{8\epsilon_{C}}{\epsilon_{L}}\big{)}^{1/4} and \eta_{\text{zpf}}=\frac{1}{2}\big{(}\frac{\epsilon_{L}}{x\epsilon_{C}}\big{)}^{1/4} have been introduced for convenience. The diagonalization basis is then

[TABLE]

with $|N\rangle$ being the Cooper pair number eigenstates of the operator $N$ and $|p\rangle/|q\rangle$ being the photon number eigenstates of $a^{\dagger}a/b^{\dagger}b$. For numerical convergence, the dimensions $N_{0}\gtrsim 7$, $p_{0}\gtrsim 7$, and $q_{0}\gtrsim 30$ are used, depending on the desired truncation accuracy. The remaining operators are written as

[TABLE]

in this basis, while the final cosine term in Eq. 8 has an exact matrix expression Smith et al. (2016), which we do not repeat here. Although the above discussion applies to the case of perfect symmetry, the discussion in Sec. III.2 allows a straightforward extension to the disordered case. Finally, we comment that, despite the high dimensionality of the matrix in Eq. 8, diagonalization can be made efficient by taking advantage of its sparsity.

Now we detail the method used to compute the wavefunctions $\langle N|\psi\rangle$ and $\langle\varphi,\phi|\psi\rangle$ (up to an overall phase) for $|\psi\rangle=|m\pm\rangle$ plotted in Figs. 3b, 3c. First, the coordinate space wavefunction $\langle\varphi,\phi,\theta|\psi\rangle$ is constructed using $\langle\varphi|N\rangle=\text{e}^{-iN\varphi}$ Devoret (1997) and the expressions for $\langle\phi|p\rangle$ and $\langle\theta|q\rangle$ as harmonic oscillator wavefunctions:

[TABLE]

where the complex coefficients $c_{Npq}$ are obtained by diagonalization (and consistent choice of overall phase). Then, the $\theta$-dependence is removed by projection onto $\theta=0$ and normalization. This amounts to computing

[TABLE]

which is plotted in Fig. 3c. Next, the $\phi$-dependence is removed by invoking Eq. 3. In this step, we project $\phi$ onto $\phi(\varphi)$ and normalize to obtain

[TABLE]

Finally, the Fourier transform of this function yields the charge wavefunctions

[TABLE]

which are plotted in Fig. 3b.

Appendix D Parity sectors

To explain the energy level structure at $\varphi_{\text{ext}}=\pi$, we factor Eq. 4 into the form Aasen et al. (2016)

[TABLE]

The parity sectors are governed by the Hamiltonians

[TABLE]

In the above, $k_{+}=0$ and $k_{-}=1$ while $\tilde{\varphi}$ and $\tilde{N}$ should be viewed as conjugate operators corresponding to pairs of Cooper pairs. Note that this expression is an exact reformulation of Eq. 4. In the limit that $\epsilon_{J}\gg\epsilon_{C}$, we may Taylor expand the potential about $\tilde{\varphi}=0$ and retain the first few terms. This step discards the effects of the offset charge, rendering $H_{+}$ and $H_{-}$ identical. We then decompose the linear part of the Hamiltonian into normal modes, which yields

[TABLE]

to leading order in $z$. The corresponding transformation is given by

[TABLE]

Eq. 9 reveals two weakly anharmonic modes: the plasmon mode at the low frequency of $\sqrt{16x\epsilon_{L}\epsilon_{C}}/h$ and a junction self-resonant mode at the high frequency of $\sqrt{8\epsilon_{J}\epsilon_{C}}/h$. These modes are coupled by a quartic nonlinearity, which has the primary effect of inducing a Kerr shift on the junction self-resonant mode. At frequencies lower than $\sqrt{8\epsilon_{J}\epsilon_{C}}/h$, the energy level structure is that of a two-fold degenerate harmonic oscillator, in agreement with the simulated energy spectrum at $\varphi_{\text{ext}}=\pi$ (see the dashed line in Fig. 3a).

Appendix E Coherence estimates

To expand on Secs. III.3 and III.4, where relaxation pure dephasing times were calculated for four values of $\delta_{L}$, we plot in Fig. 6 the results for the entire span of $\delta_{L}$. In addition, in Fig. 6c we also plot the expected relaxation, pure dephasing, and decoherence times combining the effects of all considered channels. Here, the overall decoherence time is defined using the relation $1/T_{2}=1/2T_{1}+1/T_{\upphi}$. We also note that, for the estimate of the dephasing rate due to charge noise, we do not use the extrapolation shown in Fig. 5a for the cases of disorder in $\epsilon_{J}$, $\epsilon_{C}$, and $\epsilon_{A}$. Numerical instabilities are not encounted for disorder in $\epsilon_{L}$ until $\delta_{L}\approx 0.95$ (which is outside the plotted range for both Fig. 5a and Fig. 6b).

In the remainder of this appendix, we describe the assumed frequency dependence of the three quality factors used for the relaxation time estimates in Sec. III.3. First, for the dielectric quality factor, we use the form Braginsky et al. (1987); Pop et al. (2014)

[TABLE]

so that the nominal value $Q_{\text{cap}}\sim 1\times 10^{6}$ corresponds to measurements performed at the resonant frequency of $6\,\text{GHz}$ Wang et al. (2015).

Second, for quasiparticle loss, the dissipative part of the admittance for a Josephson junction with tunneling energy $\epsilon_{J}$ is given Catelani et al. (2011) to be

[TABLE]

where $R_{\text{K}}=h/e^{2}$ is the resistance quantum, $\Delta$ is the superconducting gap, and $K_{0}$ is the modified Bessel function of the second kind. This expression holds under the assumption that the quasiparticle bath is in thermal equilibrium at temperatures $T\ll\Delta/k_{\text{B}}$, which may not be correct, depending on the implementation Serniak et al. (2018). As we have seen, the matrix elements for quasiparticle loss vanish for this circuit, so this admittance has little bearing on our relaxation time estimates.

On the other hand, inductive loss has been suggested to occur due to quasiparticle tunneling across Josephson junctions that compose superinductances Pop et al. (2014). If we expect the frequency-dependence of the dissipative part of the admittance to agree between quasiparticle and inductive loss, then we arrive at

[TABLE]

so that the nominal value $Q_{\text{ind}}\sim 500\times 10^{6}$ corresponds to measurements performed at the resonant frequency of $0.5\,\text{GHz}$.