Charge-parity switching effects and optimisation of transmon-qubit design parameters

Enhancing the performance of noisy quantum processors requires improving our understanding of error mechanisms and the ways to overcome them. A judicious selection of qubit design parameters, guided by an accurate error model, plays a pivotal role in improving the performance of quantum processors. In this study, we identify optimal ranges for qubit design parameters, grounded in comprehensive noise modeling. To this end, we commence by analyzing a previously unexplored error mechanism that can perturb diabatic two-qubit gates due to charge-parity switches caused by quasiparticles. We show that such charge-parity switching can be the dominant quasiparticle-related error source in a controlled-Z gate between two qubits. Moreover, we also demonstrate that quasiparticle dynamics, resulting in uncontrolled charge-parity switches, induce a residual longitudinal interaction between qubits in a tunable-coupler circuit. Our analysis of optimal design parameters is based on a performance metric for quantum circuit execution that takes into account the fidelity and frequencies of the appearance of both single and two-qubit gates in the circuit. This performance metric together with a detailed noise model enables us to find an optimal range for the qubit design parameters. Substantiating our findings through exact numerical simulations, we establish that fabricating quantum chips within this optimal parameter range not only augments the performance metric but also ensures its continued improvement with the enhancement of individual qubit coherence properties. Conversely, straying from the optimal parameter range can lead to the saturation of the performance metric. Our systematic analysis offers insights and serves as a guiding framework for the development of the next generation of transmon-based quantum processors.


Introduction
While quantum processors continue to progress towards practical use, the errors present in current systems are still the most limiting factor.A dominant error in superconductor-based quantum computers is decoherence.There have been several proposals to mitigate it, either by designing new qubit types [1][2][3] or by further optimizing the existing designs, typically focusing on increasing the coherence times of the circuit [4][5][6] .However, in the latter case, one often encounters trade-offs between different circuit properties.For example in transmon qubits, which have emerged as the most widely-used qubit type in large-scale experiments [7][8][9][10][11] , the suppression of charge noise comes at the cost of low anharmonicity which sets a lower bound on the duration of single-qubit operations 12 .
This illustrates the importance of understanding the different errors affecting the quantum hardware, as well as the fact that an informed design of the circuit parameters must consider a plethora of possible error sources which are not necessarily limited only to the coherence properties of the circuit, but also include leakage errors in single qubit gates due to low anharmonicity, state preparation errors due to finite-temperature heating effects [13][14][15] , as well as the parity-switching error presented in this manuscript.To elaborate further, only taking into account the coherence properties of the transmon and the low anharmonicity, one way to achieve better performance is to increase the transmon anharmonicity (in order to suppress potential leakage errors) while keeping the frequency of the transmon fixed.The latter condition, under the assumption of constant quality factors, ensures the coherence properties of the circuit remain unchanged.However this inadvertently leads into the regime where the transmon charge dispersion becomes more significant.One of the errors that is exponentially more pronounced in the low E J /E C regime is related to the presence of electron-like excitations of the superconducting condensate, referred to as quasiparticles and the charge-parity of the transmon, thus prompting us to analyze these effects further so that a trade-off between the different error sources can be made.

arXiv:2309.17168v4 [quant-ph] 29 Mar 2024
< l a t e x i t s h a 1 _ b a s e 6 4 = " m T 5 D L y V G b 5 k 0 i B 8 p 8 m j P W g z 6 1 V g = " > A A A B 8 n i c b V D L S g N B E J z 1 G e M r 6 t H L Y B C 8 G H a D r 2 P Q i 8 c I 5 g G b N c x O e p M h s z P L z K w Q 1 n y G F w + K e P V r v P k 3 T p I 9 a G J B Q 1 H V T X d X m H C m j e t + O 0 v L K 6 t r 6 4 W N 4 u b W 9 s 5 u a W + / q W W q K D S o 5 F K 1 Q 6 K B M w E N w w y H d q K A x C G H V j i 8 m f i t R 1 C a S X F v R g k E M e k L F j F K j J X 8 p + r D a U c R 0 e f Q L Z X d i j s F X i R e T s o o R 7 1 b + u r 0 J E 1 j E I Z y o r X v u Y k J M q I M o x z G x U 6 q I S F 0 S P r g W y p I D D r I p i e P 8 b F V e j i S y p Y w e K r + n s h I r P U o D m 1 n T M x A z 3 s T 8 T / P T 0 1 0 F W R M J K k B Q W e L o p R j I / H k f 9 x j C q j h I 0 s I V c z e i u m A K E K N T a l o Q / D m X 1 4 k z W r F u 6 i c 3 5 2 V a 9 d 5 H A V 0 i I 7 Q C f L Q J a q h W 1 R H D U S R R M / o F b 0 5 x n l x 3 p 2 P W e u S k 8 8 c o D 9 w P n 8 A 3 a u R A A = = < / l a t e x i t > |2 i < l a t e x i t s h a 1 _ b a s e 6 4 = " N 7 T 5 P o b j c V b 2 k z F I k L J W I I R a R x U = " > A A A B 8 n i c b V D L S g N B E J z 1 G e M r 6 t H L Y B A E I e w G X 8 e g F 4 8 R z A M 2 a 5 i d 9 C Z D Z m e W m V k h r P k M L x 4 U 8 e r X e P N v n C R 7 0 M S C h q K q m + 6 u M O F M G 9 f 9 d p a W V 1 b X 1 g s b x c 2 t 7 Z 3 d 0 t 5 + U 8 t U U W h Q y a V q h 0 Q D Z w I a h h k O 7 U Q B i U M O r X B 4 M / F b j 6 A 0 k + L e j B I I Y t I X L G K U G C v 5 T 9 W H 0 4 4 i o s + h W y q 7 F X c K v E i 8 n J R R j n q 3 9 N X p S Z r G I A z l R G v f c x M T Z E Q Z R j m M i 5 1 U Q 0 L o k P T B t 1 S Q G H S Q T U 8 e 4 2 O r 9 H A k l S 1 h 8 F T 9 P Z G R W O t R H N r O m J i B n v c m 4 n + e n 5 r o K s i Y S F I D g s 4 W R S n H R u L J / 7 j H F F D D R 5 Y Q q p i 9 F d M B U Y Q a m 1 L R h u D N v 7 x I m t W K d 1 E 5 v z s r 1 6 7 z O A r o E B 2 h E + S h S 1 R D t 6 i O G o g i i Z 7 R K 3 p z j P P i v D s f s 9 Y l J 5 8 5 Q H / g f P 4 A 2 p W Q / g = = < / l a t e x i t > |2 + i < l a t e x i t s h a 1 _ b a s e 6 4 = " 2 H 9 6 i N j Y R 7 0 X f B 3 L S / Q / e e z R H s k = " > A A A B 8 n i c b V D L S g N B E O y N r x h f U Y 9 e B o M g C G F X f B 2 D X j x G M A / Y r G F 2 M p s M m Z 1 Z Z m a F s O Y z v H h Q x K t f 4 8 2 / c Z L s Q R M L G o q q b r q 7 w o Q z b V z 3 2 y k s L a + s r h X X S x u b W 9 s 7 5 d 2 9 p p a p I r R B J J e q H W J N O R O 0 Y Z j h t J 0 o i u O Q 0 1 Y 4 v J n 4 r U e q N J P i 3 o w S G s S 4 L 1 j E C D Z W 8 p + 8 h 5 O O w q L P a b d c c a v u F G i R e D m p Q I 5 6 t / z V 6 U m S x l Q Y w r H W v u c m J s i w M o x w O i 5 1 U k 0 T T I a 4 T 3 1 L B Y 6 p D r L p y W N 0 Z J U e i q S y J Q y a q r 8 n M h x r P Y p D 2 x l j M 9 D z 3 k T 8 z / N T E 1 0 F G R N J a q g g s 0 V R y p G R a P I / 6 j F F i e E j S z B R z N 6 K y A A r T I x N q W R D 8 O Z f X i T N 0 6 p 3 U T 2 / O 6 v U r v M 4 i n A A h 3 A M H l x C D W 6 h D g 0 g I O E Z X u H N M c 6 L 8 + 5 8 z F o L T j 6 z D 3 / g f P 4 A 2 Q i Q / Q = = < / l a t e x i t > |1 + i < l a t e x i t s h a 1 _ b a s e 6 Y w S Y 6 W H J 6 + j i O x z 7 J b K X s W b w l 0 k f k 7 K k K P W L X 1 1 e j F N B U p D O d G 6 7 X u J C T K i D K M c x 8 V O q j E h d E j 6 2 L Z U E o E 6 y K Y H j 9 1 j q / T c K F a 2 p H G n 6 u + J j A i t R y K 0 n Y K Y g Z 7 3 J u J / X j s 1 0 V W Q M Z m k B i W d L Y p S 7 p r Y n X z v 9 p h C a v j I E k I V s 7 e 6 d E A U o c Z m V L Q h + P M v L 5 L G a c W / q J z f n Z W r 1 3 k c B T i E I z g B H y 6 h C r d Q g z p Q E P A M r / D m K O f F e X c + Z q 1 L T j 5 z A H / g f P 4 A u G m Q X w = = < / l a t e x i t > |0i < l a t e x i t s h a 1 _ b a s e 6 4 = " j 7 e 4 3 Q a / j A 5 Y F o 6 E b N k t D 2 y Y q h 0 = " > A A A C A H i c b V D L S s N A F J 3 U V 6 2 v q A s X b g a L U J G W p P h a F k V w W c E + o I 1 h M p 2 0 Q y e T M D M R S u j G X 3 H j Q h G 3 f o Y 7 / 8 Z J m o V W D w y c e 8 6 9 3 L n H i x i V y r K + j M L C 4 t L y S n G 1 t L a + s b l l b u + 0 Z R g L T F o 4 Z K H o e k g S R j l p K a o Y 6 U a C o M B j p O O N r 1 K / 8 0 C E p C G / U 5 O I O A E a c u p T j J S W X H P v 2 q 3 f H 1 e 4 O z y C V Z g W 1 a x w z b J V s z L A v 8 T O S R n k a L r m Z 3 8 Q 4 j g g X G G G p O z Z V q S c B A l F M S P T U j + W J E J 4 j I a k p y l H A Z F O k h 0 w h Y d a G U A / F P p x B T P 1 5 0 S C A i k n g a c 7 A 6 R G c t 5 L x f + 8 X q z 8 C y e h P I o V 4 X i 2 y I 8 Z V C F M 0 4 A D K g h W b K I J w o L q v 0 I 8 Q g J h p T M r 6 R D s + Z P / k n a 9 Z p / V T m 9 P y o 3 L P I 4 i 2 A c H o A J s c A 4 a 4 A Y 0 Q Q t g M A V P 4 A W 8 G o / G s / F m v M 9 a C 0 Y + s w t + w f j 4 B o n S k 9 I = < / l a t e x i t > E + 2 (n g ) E 2 (n g ) n V C S p J g J P B 0 U p g z q G e S S w R y X B m o 0 M Q V h S s y v E A y Q R 1 i a 4 k g n B n T 9 5 k b T O q u 5 l 9 a J x X q 5 V Z n E U w Q E 4 A h X g g i t Q A 7 e g D p o A g 0 f w D F 7 B m / V k v V j v 1 s e 0 t G D N e v b B H 1 i f P 6 l f l q g = < / l a t e x i t > Qubit 1 (q 1 ) < l a t e x i t s h a 1 _ b a s e 6 4 = " 3 3 y n z N K C P 2 z q T i B V B H M I s E 9 q k k W L O J I Q h L a n a F e I Q k w t o E V z I h u I s n / y X t W t W 9 q J 4 3 z 8 r 1 S h 5 H E R y A I 1 A B L r g E d X A D G q A F M H g A T + A F v F q P 1 r P 1 Z r 3 P S w t W 3 r M P f s H 6 + A a s c p a q < / l a t e x i t > Qubit 2 (q 2 ) < l a t e x i t s h a 1 _ b a s e 6 4 = " I S V P B d + X L Z J + y q p k e 0 r O X x g N U H 4 = " < l a t e x i t s h a 1 _ b a s e 6 4 = " 3 3 y n z N K The lumped element model of the tunable-coupler circuit used in the implementation of diabatic CPHASE gates, consisting of two computational transmons (dark and light blue) referred to in the following as Qubits 1 and 2 (with indices q 1 and q 2 ), a flux-tunable coupler (green), denoted with the index c, and capacitive couplings between the transmons (black).The readout resonators and drive lines for the implementation of single-qubit gates are not included in the schematics.
Such detrimental quasiparticle-induced effects have, in turn, motivated research for finding mitigation strategies such as normal-metal traps [24][25][26] , band-gap engineering [27][28][29] and improved qubit design 30 .There have also been efforts towards designing new types of superconducting qubits that are expected to be intrinsically robust against quasiparticle tunneling 31 .
Suppression of charge-noise susceptibility of the transmon is achieved by adding a large shunt capacitor in parallel with a Josephson junction of a Cooper-pair box 12 .However, the energy levels of the transmon exhibit a weak 2e-periodic charge dispersion and are, thus, not completely independent of the offset charge on the transmon island.Since the presence of a quasiparticle shifts the island charge by e, the energy spectrum of the transmon can, thus, be divided into two distinct manifolds based on the parity of the number of quasiparticles on the island as shown in Fig. 1a.
Switching of the charge-parity can occur either due to a preexisting quasiparticle tunneling across the Josephson junction on to the transmon island, or due to photon-assisted breaking of a Cooper pair, as depicted in Fig. 1b 16 .The timescale of these stochastic parity-switches is referred to as the parity-switching time and, as we argue in the Results section, it is typically much shorter compared to the time needed to obtain meaningful statistics from the quantum computer.
Since the difference between the two parity manifolds is strongly suppressed, quasiparticle effects in transmons have mainly been analyzed in the context of quasiparticle-induced decoherence.However, the parity-dependent energy splitting of the higher-excited states is much larger compared to that of the first-excited and ground states, as shown in Fig. 1a.Therefore, parity-switching in the second excited state can potentially become a notable source of error for example in diabatic CZ (controlled-Z) or more generally CPHASE (controlled-phase) gates 32 .There are different implementations of this gate scheme, most notably, statically coupled transmons 33 and parametric drives 34 , tunable couplers and parametric drives 35 and the frequency-tunable coupler architecture [35][36][37][38][39][40][41][42][43][44] .Nonetheless, in all of these gate implementations, the second excited state of one transmon is populated during the gate operation, suggesting the gate could be susceptible to the charge-parity of one transmon in the system.
Besides the ability to perform fast and high-fidelity gate operations, one of the main reasons for the introduction of the tunable coupler is the fact that the static residual-ZZ interaction between the qubits can be completely suppressed by tuning the coupler transmon to a specific frequency 36,38 , thus leading to reduced cross-talk and spectator qubit related decoherence 45 .When determining this frequency, it has been shown that the level repulsion between the higher levels of the tunable coupler system are relevant 36,39 , thus implying that quasiparticle dynamics might affect our ability to effectively decouple the qubits.On the contrary, these level repulsions can serve as a means to enhance the ZZ interaction strength and facilitate the implementation of CPHASE gates entirely adiabatically, without populating the higher excited states 46 .Building upon this insight, we anticipate that parity switching may likewise exert a significant influence on the performance of adiabatic gate implementations.
In this paper, we develop an analytical theory of parity switches in a tunable coupler based architecture which is currently one of the most promising platforms for large-scale quantum computing 7,8,10,11 .We demonstrate that the effects of a parity switch in a two-qubit gate can be a relevant source of error, even in the transmon regime.Moreover, we show that this previously unidentified error can, in certain parameter regimes, be the dominant quasiparticle-induced error mechanism during a diabatic two-qubit gate, indicated by a comparison to currently achievable parity switching rates observed in superconducting circuits.Furthermore, we demonstrate that the inherent stochastic nature of parity-switching events limits the ability to suppress any unwanted longitudinal interactions between the qubits coupled through tunable couplers 47 .We find that the magnitude of the unwanted interactions make this effect relevant as coherence times advance into the millisecond regime 48,49 .
Secondly, we introduce a systematic and versatile analysis aimed at identifying the optimal qubit design parameters.This comprehensive analysis involves accurately modeling dominant error sources and utilizing a performance metric tailored specifically for the quantum circuit under consideration.Importantly, we utilize analytical expressions for various infidelity terms within the performance metric, which significantly aids the scalability of our approach.By examining the performance metric's reliance on the transmon qubit design parameters, we delineate an optimal parameter range that maximizes overall performance.We demonstrate how our analysis and findings can offer invaluable guidance in the advancement of transmonbased quantum processors.Furthermore, since the established methodology in this paper is general-purpose, we anticipate its broader applicability across various hardware platforms.

Charge-Parity Modelling
In this section, we first present how a parity switch affects a single transmon qubit.The Hamiltonian of an individual transmon, not taking into account potential higher order contributions to the Josephson energy 50 , is given by 12,51 where the operator n represents the dimensionless charge and φ is the superconducting phase operator across the Josephson junction.The variables E C , E J and n g represent the charging energy of an electron (i.e. the energy required to add a single electron of the Cooper-pair to the transmon island), Josephson energy and dimensionless offset charge, respectively.The variables n and φ are related via the canonical commutation relation [ φ , n] = i.Additionally we have included a discrete parity variable P ∈ {−1, +1}, corresponding to the parity of the number of electrons that have tunneled across the junction.The parity term has the same effect as a shift of the offset charge by ∆n g = 1/2.Denoting the eigenenergies of the original transmon Hamiltonian in Eq. 1 as E i with i ∈ {0, 1, 2, 3, ...}, the difference between the energy levels of the different parity states can be asymptotically approximated by 12 where the superscript refers to the parity and the charge dispersion ε m is given by While the exponential suppression of the charge dispersion with the ratio E J /E C is well-known and the main reason for the introduction of the transmon, the formula in Eq. 3 also predicts a significant increase in the charge dispersion of higher excited states.This means that even though the effect of a parity switch may be small in the computational subspace, the effect can be significantly more pronounced if higher-excited states are involved in the operation of two-qubit gates.For example, we find |ε 2 /ε 1 | ∼ 40 for E J /E C ∼ 50.This difference is even more pronounced since certain effects scale with the square of the charge dispersion, as we show in the following.Consequently, the Hamiltonian of the single transmon in Eq. 1 can be approximated in the low-energy manifold and in the asymptotic limit of E J /E C ≫ 1 as where â are bosonic annihilation operators, and hω 22 are the parity-averaged expressions for the transmon (angular) frequency and the anharmonicity in the asymptotic limit, respectively.Here, we have taken into account the fact that the different parities have almost identical parameters, and the small differences between them are taken into account with the two parameters depending on the parity, δ ω(P, n g ) and δ α(P, n g ).The parity P therefore divides the eigenstates of the Hamiltonian in Eq. 1 into two distinct manifolds, as illustrated in Fig. 1a.Due to the rapid scaling of the charge dispersion shown in Eq. 3, we neglect in the following the effect of the parity switching on the first excited state, i.e. we set δ ω = 0, and only focus on the second excited state.Consequently, δ α(P, n g ) = Pε 2 cos(2πn g )/(2h).

The Diabatic CPHASE Gate
][43] .We show the circuit schematics of the tunable-coupler setup in Fig. 1d.Here, the two computational transmons, which we refer to as Qubits 1 and 2 (q 1,2 ), are capacitively coupled with each other and to a third, frequency-tunable, transmon which is referred to as the tunable coupler or simply coupler (c).
The main operation principle of the diabatic CPHASE gate is shown in Fig. 1c.The CPHASE gate is implemented by tuning the frequency of the coupler closer to the frequency of the computational transmons by using a flux pulse.The conditional phase is collected during a Rabi oscillation between the |11⟩ and |02⟩ states of the computational qubits, as illustrated in Fig. 1c.We model the circuit pictured in Fig. 1d with the Hamiltonian 38,39 Since the second excited state of one of the computational transmons is significantly populated during the gate operation, the non-degeneracy of the two parity levels can have a direct effect, and therefore quasiparticle tunneling and photon-assisted pair breaking in the transmon can become a notable source of error in the gate operation.Further details about the numerical modeling and definition of the computational basis in relation to the states of the three constituent transmons are found in the Methods section.

Effective Model
The tunable coupler circuit Hamiltonian in Eq. 5 is difficult to analyze and we must often rely on numerical studies 39 .It is therefore beneficial to introduce an effective Hamiltonian that can approximate the physics of the system.Similar to Refs. 38,52  we introduce the Schrieffer-Wolff transformation as a means to decouple the computational transmon states from the coupler states, and assume that the decoupled coupler remains in the ground state during the gate operations.Unlike the approach taken in Refs.38,52 , where the Hilbert space of the local transmons is truncated to the computational subspace, we also include the |02⟩ state to account for the Rabi oscillation that is used to accumulate the conditional phase.More details of the Schrieffer-Wolff transformation can be found in the Methods section, where we show that the diabatic gate can be modelled with the following effective unitary Here, we have denoted the conditional phase with φ (t) and the population of the |11⟩ state with P 11 .Note that this operator is not necessarily trace-preserving, as part of the population of the |11⟩ state might remain in the |02⟩ state, due to potential 4/27 calibration errors.More explicitly, P 11 is given by where we have defined the qubit-qubit detuning and the Rabi frequency of the |11⟩ ↔ |02⟩ transition as ∆ = ωq 1 − ωq 2 and 11,02 , respectively.The parameters with the tilde denote the perturbed parameters of the original full Hamiltonian from Eq. 5, which were derived by applying the Schrieffer-Wolff transformation (see Methods, more specifically Eqs.32-35).The conditional phase in Eq. 6 is given by Eqs. 7 and 8 can be used to assess the susceptibility of the gate parameters to a small perturbation, such as a parity switch.Due to the larger charge dispersion of the second excited state, and the fact that the second excited state of Qubit 2 (q 2 ) is populated during the gate operation, we can assume that the main contribution of the parity switch is the perturbation of the anharmonicity α q 2 .
By treating the parity-dependent contribution to the anharmonicity δ α from Eq. 4 as a small perturbation, we can obtain the parity-dependent expressions for the conditional phase and |11⟩ population P 11 (t g , P q 2 ) ≈ 1 + 1 2 where the Taylor expansion of the optimal gate parameters for a small perturbation of α q 2 evaluated at the mean (parity averaged) anharmonicity α q 2 from Eq. 4 was employed.At this point, we stress again that δ α(P, n g ) = Pε 2 cos(2πn g )/2.While the above expressions are completely general also in the non-perturbative regime, the relations given in Eqs.7 and 8 can be used to obtain analytical expressions for ∂ φ /∂ α q 2 and ∂ 2 P 11 /∂ α 2 q 2 , which determine the susceptibility of the gate to charge-parity switches.
We realize from Eq. 7 that the implementation of a high-fidelity gate with an arbitrary conditional phase φ 0 requires P 11 (t g ) = 1, otherwise some population remains outside of the computational subspace.Therefore, the following condition for the gate time t g must hold: Ωt g = n • 2π, in which n ∈ N is an integer number.This condition enables us to further simplify the relations for the susceptibility of the conditional phase to a parity switch, and we obtain up to the leading order, and 11,02 In the derivation of Eq. 11 we have neglected the terms proportional to g q 1 c g q 2 c /(ω q 2 − ω c + α q 2 ) 3 and g 2 q 2 c /(ω q 2 − ω c + α q 2 ) 3 and higher orders.Additionally Eq. 12 only contains the lowest order scaling of the result.A list of all the assumptions used in the derivation of Eqs.11 and 12 is included in the Methods section.

Gate Fidelity Limitations
So far we have quantified how a parity switching event can affect the parameters of the gate unitary.In order to describe the gate performance in a quantum circuit, we also need to consider how frequently parity switching events occur.
Experiments of parity switching lifetimes typically show the parity switching time to lie in the broad range of T P ∼ 100 µs − 1 s 51,[53][54][55][56][57] .Even though the parity switching lifetime of transmons might increase in the future, e.g., due to better design and improved shielding, it appears that the parity lifetime may be fundamentally upper bounded by the high-energy quasiparticle burst events, which are observed to happen once every 10 − 50 s [58][59][60][61] .

Kraus Operator Description
Comparing the realistic range of parity lifetimes of superconducting qubits to the duration of a single two-qubit gate t g , which is typically in the range of tens to hundreds of nanoseconds 36,37,41,43,62 , we observe that t g ≪ T P .However, any meaningful application of a quantum computer will include re-running an algorithm, comprised of a large number of non-parallel gates N gates , in order to reduce the statistical uncertainty of the observable being evaluated.Denoting the number of runs by N shots , the uncertainty of any observable generally scales as 1/ √ N shots 63 and therefore N shots must be large.All together, the time needed to execute a full algorithm T alg roughly scales as T alg ∝ N shots • N gates • t g , but it may be realistically even longer, due to the time needed to measure and reinitialize the quantum computer, and any possible pulse schedule compilation of the control electronics 64 .
As an example, in Refs. 7,9,65 sate of the art devices were run for a total of several minutes in order to obtain meaningful results, which implies the following clear separation of timescales The left side of Eq. 13 indicates that the probability for a parity switch occurring during the operation of a single gate is very low, while the right-hand side suggests that a large number of parity switches can occur during an execution of an algorithm.This means that the effect of the charge-parity switch (CPS) on a diabatic CPHASE gate can be described by the following Kraus operators 66 acting on the two-qubit density matrix ρ where Û± are Kraus operators corresponding to the different parity implementations of the two-qubit gate.Eq. 14 can be interpreted as a stochastic application of two different gate operators; by assuming that the target conditional phase is φ 0 , they can be written as with δ φ = ∂ φ /∂ α q 2 ε q 2 2 cos(2πn q 2 g )/h and δ P 11 = ∂ 2 P 11 /∂ α 2 q 2 ε q 2 2 cos(2πn q 2 g )/h 2 /2, which is a result of Eqs. 9 and 10, together with δ α(P, n g ) = Pε 2 cos(2πn g )/(2h).We have additionally assumed that both parities are equally likely, however this assumption can also be easily relaxed.Note that the channel defined in Eq. 14 is not trace preserving when δ P 11 > 0.More specifically tr{CPHASE CPS [ ρ]} ≤ tr{ ρ}, where the equality holds when δ P 11 = 0.In this case the map is both completely positive and trace preserving and the completeness relation of the Kraus operators ∑ i=+,− Û † i Ûi = I holds.We have defined the Kraus operators in Eq. 15 so that there are small errors associated with each parity state.We show in the following that this corresponds to higher average gate fidelities, compared to having one parity state with a perfect fidelity and the second parity state with a larger error.This assumption is, therefore, equivalent to optimally calibrating the gate with respect to the parity switching error.

Gate Fidelity
Using the Kraus operator description in Eq. 15 allows us to make a formal analysis of performance metric of an arbitrary conditional phase gate in the presence of parity switches and resulting phase and leakage errors.
The average gate fidelity F in the presence of leakage can be computed by 67,68 with the leakage parameter Here ÛCPHASE is the unitary operator of an ideal CPHASE gate and d is the dimension of the computational Hilbert space, which in our case is d = 4.The fidelity F of the operation given in Eq. 14 for small perturbations δ φ and δ P 11 is given by the following expression

6/27
Here, we observe that the effect on the fidelity is of second order in the charge dispersion, due to the coherent nature of a conditional phase error.However, we have shown in Eq. 38 of the Methods section that the infidelity of a series of N gates is given by 1 − F N ≈ N(1 − F ).This indicates that the error scales linearly with the number of gates, as is typical of incoherent errors, but quadratically in the error parameter δ φ , as is expected from a coherent error 69,70 .
In deriving Eq. 38, one can observe that the numerical prefactors in front of both terms (here 3/80 and 1/32) increase, and therefore the average gate fidelity decreases, if the gate is not calibrated in such a way that the error is equally distributed between both parity states.For this reason, we have chosen the Kraus operators as given in Eq. 15.
We conclude from Eq. 17, together with Eqs.11 and 12, that the dominant effect of the parity switch event is the shift in the conditional phase, rather than leakage, since q 2 in the perturbative regime considered in this manuscript.Moreover the magnitude of the shift in conditional phase is given by Eq. 9 and Eq.11.Further simplifying Eq. 38 therefore results in Numerical simulations Here, we compare the above results with the numerically exact treatment of full the Hamiltonian in Eq. 5.The parity switch effect is taken into account in the numerical experiments by considering the transmon Hamiltonian from Eq. 4. In other words, the system is simulated for all 2 3 = 8 possible parity states and the results are averaged accordingly.According to Eq. 2, the magnitude of the effect also depends on the offset voltage of Qubit 2, which is typically not known.Thus, we assume for simplicity that n q 2 g = 0.Alternatively, as long as all the effects remain second order in the charge dispersion, one can also define an average charge dispersion as εq

Parity-Switch Induced Diabatic Gate Errors
To achieve high-fidelity gate simulations, it is crucial to carefully select the Hamiltonian parameters in Eq. 5. We have described how these parameters were chosen so that high-fidelity gates are possible with arbitrary ratios of E J /E C of Qubit 2 in Table 1 of the Methods section.We also note that in general, we use parameters that closely resemble those in the implementation presented in Ref. 36 .We examine the effect of a parity-switch on the fidelity of a CZ gate and compare the perturbative analytical results to a full simulation of the Hamiltonian in Eq. 5.The pulse shape used in the simulations (defined in Eq. 27) is parameterized by its amplitude A, ramping up time τ b and plateau duration τ c .The full duration of the flux pulse is therefore given by T = τ c + 2τ b .The conditional phase of the gate in the simulations is obtained from propagating the state |ψ(t = 0)⟩ = 1 2 (1, 1, 1, 1) T and extracting the conditional phase of the |11⟩ state.Fig. 2a compares the analytical results in Eqs. 9 and 11 to the numerically obtained values and confirms that the two approaches agree up to a good accuracy.Since the magnitude of the charge dispersion in the numerical analysis is the same as in the analytical treatment, the error stems completely from the approximations made in evaluating ∂ φ /∂ α q 2 .Furthermore, Fig. 2a also clearly shows that the parity-switching-induced shift in the conditional phase scales exponentially with the ratio of E J /E C .This is due to the scaling of the charge dispersion in Eq. 3.
In Fig. 2b, we show numerical data of the full gate fidelity after the parity switch and, as a comparison, the corresponding analytical result calculated using Eq.18.The average gate fidelity in the numerical example is obtained from propagating a number of input states, reconstructing the effective superoperator of the gate from these simulations and subsequently using Eq.16 to obtain the average gate fidelity of the gate.We observe that the infidelity of the numerical simulation flattens for E J /E C ≳ 75, which is due to other errors in the gate implementation, such as leakage transitions during the pulse ramping up and down, which was confirmed by monitoring the population of the computational states during gate operation (see Supplementary Information for more details).In this region (in grey) the effect of a parity switch is not seen since it is too small compared to other errors.On the other hand at lower values of E J /E C , the numerical results overlap with the parity switching error predicted by the shift in the conditional phase.While clearly demonstrating the magnitude of the error, this result also shows that the leakage error contribution in Eq. 17 is negligible in the perturbative regime.The gate durations in Fig. 2 are typically 45 ns ≲ t g ≲ 60 ns, with τ c ∼ t g .
Additionally, Fig. 2b compares the magnitude of the error due to quasiparticle related decoherence to the parity switch induced error described in this work.Since the quasiparticle induced characteristic decay times T 1 and T φ depend on a large number of parameters and there are two possible quasiparticle generating mechanisms (see Fig. 1b) 16 , we only provide an upper bound based on the parity switching time observed in the references cited in the caption of Fig. 2.This upper bound is determined by noting that in the computational subspace 55  56 .More specifically, references [1] and [2] have reported the values T P = 1.25 ms and T P = 2.5 ms respectively.In reference [3], parity-switching times between 1 ms and 1.5 s were demonstrated.For this reference, we have used the value of T P = 20 ms, which approximately corresponds to the median parity-switching time of all the samples.The grey region is the region in which errors due to unwanted transitions during the gate operation are more prominent and the effect of a parity switch is negligible.It therefore represents the lower bound of the infidelity, which is largely independent of the ratio E J /E C .
where the rates Γ +− i j represent the transition rates between states |i ± ⟩ → | j ∓ ⟩ in different parity manifolds of a single transmon.Since amplitude damping noise has a larger effect on the fidelity compared to pure dephasing (see Table 2), we furthermore assume the worst-case scenario in which each quasiparticle-induced parity switching event results in a T 1 decay, so that 1/T qp 1 = Γ +− 01 + Γ +− 10 ≈ 2/T P .This expression for T qp 1 , together with experimentally measured values of T P , provides an approximate upper bound for the magnitude of the effect of the decoherence.Even though each green line in Fig. 2b uses a constant measured T P , independent of E J /E C , the plotted infidelity contribution is not constant due to the varying gate duration We observe from Fig. 2b that depending on the ratio of E J /E C and the parity lifetimes T P , the contribution of parity switching to the infidelity of the two-qubit gate system can dominate the contribution from quasiparticle-induced decoherence.Note that the two-qubit gates are also typically the noisiest building blocks of a quantum algorithm 7,8,11,69 .We have also provided additional numerical results in the Supplementary Information, showing that the effect of a parity switch on the leakage becomes the main contribution to the infidelity at shorter gate times.
In order to compare the effect of the charge-parity switches to 1/ f -charge noise, which was previously thought to dominate the low E J /E C regime, we included a comparison of the infidelity of a two-qubit gate due to 1/ f -type charge noise in the Supplementary Information, since the contribution of both of these error sources scales with the charge dispersion of the transmon.However, unlike the charge dispersion of the transmon which is maximized when cos(2πn g ) = 1, the low-frequency charge noise decoherence rate is maximal at sin(2πn g ) = 1 .This means that these two errors are mutually exclusive, i.e. if we were hypothetically able to tune n g to the value where the charge dispersion of the transmon is equal to zero, that point corresponds to the maximal decay rate due to low-frequency charge noise 12 .Nonetheless, in the comparison we have assumed in both the charge-parity switching analysis and the low-frequency charge noise analysis that we are at the noise hotspot, i.e. the value of n g where the effects are maximal, thus slightly overestimating both errors.We observe that at lower E J /E C the charge-parity switching error is dominant and vice versa.The crossover between the infidelity due to the charge-noise-induced dephasing in the computational subspace and the parity switching exhibits a crossover at E J /E C ∼ 80, which corresponds to a gate infidelity of 1 − F ∼ 10 −8 , i.e. the charge-parity switching error is dominant for E J /E C ≲ 80.However, due to the increased charge dispersion in the higher-excited levels of the transmon and the utilization of the second excited state, we have additionally analyzed the effect of the second excited state charge-noise-induced dephasing which was found to be much larger compared to the effects in the computational subspace.In this case, the crossover between the charge-parity error, which was again dominant at lower E J /E C , was found to occur at E J /E C ∼ 65 and infidelities on the order of 1 − F ∼ 10 −6 .We note here that our analysis overestimates the effect of the charge noise in a realistic scenario by neglecting the time correlations in the noise and should be treated as an upper bound of the magnitude of the effect.These results demonstrate that the parity-switching error is the dominant charge dispersion related error of the two-qubit gate in the regime of E J /E C ≲ 65.

Parity-Switching Effects on the Adiabatic Interaction
The tunable coupler architecture analyzed thus far can also be used to implement an adiabatic CPHASE gate.In this case, the level repulsion between the states of the double-excitation manifold is used to obtain a phase shift of the |11⟩ state, which accumulates the conditional phase with the rate where hω i j are the eigenenergies of the eigenstates which form the computational basis, as described in the Methods section.In this gate implementation the population ideally never leaves the computational subspace, however we will show that, also in the adiabatic implementation, the gate will suffer from charge-parity switching errors of a similar magnitude as in the diabatic case.Additionally, one of the main benefits of the tunable coupler architecture is the ability to also suppress any residual interactions between the qubits during idling 38,39,47 , with unwanted ZZ coupling strengths demonstrated to be below 1 kHz 36 .In this section, we also show that uncontrolled parity-switching sets a lower bound on the minimum achievable unwanted ZZ interaction in such systems.
In the presence of parity-switches in all three transmons of the tunable coupler system, the ZZ coupling rate of each of the distinct parity configurations will slightly differ.Since these differences are proportional to the charge dispersion, and therefore small, we can describe the parity-dependent ZZ coupling of this system with a first-order Taylor expansion, so that with P i ∈ {−1, +1} ∀i.While we have still assumed ε 2 ≫ ε 1 , we have also acknowledged that the first-excited-state charge dispersion and the derivatives ∂ ζ ZZ /∂ ω i can, in certain cases, be significantly larger compared to ∂ ζ ZZ /∂ α i .ζ ZZ therefore corresponds to the ZZ coupling strength of the system with parity-averaged parameters.However, in order to evaluate the derivatives in Eq.21, or the general value of the coupling rate ζ ZZ with fixed Hamiltonian parameters (without including parity effects), it is necessary to go to fourth-order perturbation theory in the coupling strengths, as was done in Refs. 36,39 While the complete fourth-order expression is impractical for obtaining any analytical insight, by further assuming the hierarchy of the system parameters Σ ic ≫ ∆ ic ≫ g ic ≫ g 12 , where Σ i j = ω i + ω j and ∆ i j = ω i − ω j , for i, j ∈ {q 1 , c, q 2 }, the cumbersome expressions can be significantly simplified into 39 ζ ZZ ≈ 2 (α q 1 + α q 2 ) g2 01,10 − 2ν g01,10 (2α where ν = g q 1 c g q 2 c /(2∆ q 1 c ∆ q 2 c ) ∼ 10 −3 .Focusing on the idling regime, where ζ 0 ZZ = 0, the difference between the derivatives obtained from Eq. 22 and the numerically exact result is shown in Fig. 3a.We observe qualitatively good agreement in the vicinity of the two coupler idling frequencies of the system (vertical black dotted lines).
Since the system has eight uncontrolled, rapidly (compared to experimental timescales, as in Eq. 13) switching parity states, we define a parity-averaged idling interaction strength: This definition is the parity-averaged mean-square of the idling ZZ strength defined in Eq. 21.Parity-averaged coupling strengths are shown in Fig. 3b, where we observe a good agreement between the analytical result and the numerical data.The possible values of the quantity ζ ZZ (P q 1 , P c , P q 2 ) for the third data point [also shown in panel a] are visualized in the bar graph.
The bar graph demonstrates that Qubit 2 (q 2 ) is the dominant contributor to the perturbation, with small corrections due to the parity of Qubit 1 (q 1 ).This can be seen by noting that the residual coupling strength is dependent mostly on the parity of Qubit 2 (q 2 ).We have extended this analysis to larger values of ∆ q 1 q 2 in the Supplementary Information, to show that the parity-switching induced ZZ coupling strength is not always dependent on the parity of a single transmon in the system.The parity of the coupler in this set of parameters is largely irrelevant due to the high ratio of E J c /E C c ≈ 250, but this is not always the case.For example, considering the gate implementation from Ref. 41 , the ratio of E J c /E C c for the coupler is much smaller and therefore, the coupler parity in this implementation has a much larger effect on the strength of the residual ZZ coupling.The above results therefore present a fundamental limit on the magnitude of unwanted interactions that can be achieved in the tunable coupler architecture.However, the overall effect on an algorithm is more complex as it depends on the duration of the execution (since the interaction is always "on"), which in turn depends on the coherence times.We estimate that this error becomes relevant if the coherence times are of the order of 1/ ⟨ ζ 2 ZZ ⟩ CPS , since a significant unwanted conditional phase is accumulated if the algorithm duration ∆T (which is determined by the coherence times) is long enough such that . For E J /E C ∼ 50, the coherence time (or algorithm duration) needed to observe the parity-switchinginduced residual coupling strength is on the order of 1 ms.Note that current fabrication processes are indeed already approaching this value 48,49 .
Moving on the adiabatic CPHASE gate, we have observed in Fig. 3a,b that the state contributing most to the unwanted ZZ interaction in the idling regime is the second excited state of the high-frequency qubit, provided that ω q 1 − ω q 2 ≃ α q 2 .This is 10/27 due to the hybridization of the computational |11⟩ state with the local transmon state |0 q 1 0 c 2 q 2 ⟩, which was also confirmed by the additional results included in the Supplementary Information.More importantly, this means that according to the same reasoning as in the diabatic gate case, the adiabatic gate can be described by the same pair of Kraus operators as defined in Eq. 15, and since we are assuming adiabaticity, we can automatically set δ P 11 = 0. Again, applying the same reasoning as in the diabatic case, the fidelity of the adiabatic gate is given by In the second step we have again assumed a flat pulse with duration t g and neglected the ramping up and down phases of the flux pulse.With this simplification, the conditional phase is given by φ 0 = ζ 0 ZZ t g .Contrasting the expressions in Eq. 24 and Eq. 18, we can see that if the quantity ∂ ζ ZZ /∂ α q 2 • φ 0 /|ζ 0 ZZ | is on the order of ∼ 10 ns, the error of the adiabatic gate is similar to the diabatic case (displayed in Fig. 2).More specifically, the conditional phase shift due to a charge-parity switch is given by δ φ d ≈ t g /2 • ε q 2 2 /h in the diabatic case (according to Eq. 11) and by δ φ /h in the adiabatic case (ignoring the factor cos(2πn q 2 g )).3c, where we can see that in the vicinity of the horizontal line where ω q 1 − ω q 2 = α q 2 , the hybridization with the second excited state results in an increased sensitivity to the charge-parity switches.However, the exact value of ∂ ζ ZZ /∂ α q 2 • φ 0 /|ζ 0 ZZ | is very dependent on the value of the qubit-qubit detuning, but generally the regime where ω q 1 − ω q 2 < α q 2 is more favorable.In the regime where ω q 1 − ω q 2 ≈ α q 2 , the observed effect of a charge-parity switch on the fidelity (which is proportional to the square of the plotted values according to Eq. 24) can be an order of magnitude larger than in the diabatic case.Conversely, in the regime where ω q 1 − ω q 2 < α q 2 , the infidelity is up to two orders of magnitude smaller.More specifically, at E J /E C = 50 (corresponding to ε q 2 2 ≈ 250 h kHz), the maximal fidelity achievable in the range of the plot in Fig. 3c is F ∼ 1 − 10 −7 , and the minimal fidelity is F ∼ 1 − 10 −2 .

Optimal Qubit Parameters
Having established the magnitude of the parity-induced error on a two-qubit gate, we have shown that this error can be mitigated by increasing the E J /E C ratio of the transmon.However, there are other error sources present in such architectures 40 , and while increasing the E J /E C ratio will suppress the parity-switching errors, it may also increase the contribution of other possible error sources.Therefore, in order to find better parameters for future transmon-based quantum computers, we must evaluate the contributions of all errors affecting the system.In particular, we estimate optimal regions for the qubit parameters E J and E C , where the errors contributing to the gate and state preparation infidelities are minimized.
5. Errors in the two-qubit gate operation due to parity switch effects that are analyzed in this manuscript.
In Methods, we show how the above error sources scale with the transmon Hamiltonian parameters E J and E C .We have not included any errors related to the control and calibration of the individual gates, as such errors do not explicitly depend on the qubit parameters and their inclusion, therefore, would not significantly alter the presented results.Similarly, measurement errors are present, but have no explicit dependence on E J and E C .For simplicity, we consider the same tunable coupler system as in the previous section, but arranged in a square grid, as pictured in Fig. 4a.In order to avoid frequency crowding issues with such a connectivity 82 , the qubits in the array are divided into low and high-frequency transmons.Furthermore, to prevent next-nearest-neighbor frequency collisions, even qubits assigned to the high-or low-frequency groups must possess non-identical frequencies, as demonstrated in Ref 83 .Nonetheless, the frequency spread within each group is considerably smaller compared to the disparity between the means of these two groups.While our analysis focuses on a single qubit pair, extending the results to encompass a parameter distribution is straightforward.Additionally, we will define a broader region of optimal parameters spanning several hundred MHz in frequency, ensuring that all transmons within the square grid can fit within this designated range.
q q Y E J z l l 1 d J p 1 F 3 r u q X D x e 1 5 m 0 R R x l O 4 B T O w Y F r a M I 9 t K A N B J 7 g G V 7 h z Z p a L 9 a 7 9 b F o L V n F z D H 8 g f X 5 A x q + k l U = < / l a t e x i t > ⇡ 2

< l a t e x i t s h a 1 _ b a s e 6 4 = " Z n v B d M L j C d 5 R N T t 3 E X P Y h Y C G t Z Y = " >
q q Y E J z l l 1 d J p 1 F 3 r u q X D x e 1 5 m 0 R R x l O 4 B T O w Y F r a M I 9 t K A N B J 7 g G V 7 h z Z p a L 9 a 7 9 b F o L V n F z D H 8 g f X 5 A x q + k l U = < / l a t e x i t > ⇡ 2 < l a t e x i t s h a 1 _ b a s e 6 4 = " s r X F r C w s r p j a m     Values of E J and E C of Qubit 2 (q 2 ) from different experimental implementations of the tunable coupler transmon architecture.Both parameters are extracted from the reported ω q 2 and α q 2 , and therefore the points are only an approximation 50 37 and [6] Google Quantum AI (2022) 8 .The black lines represent constant ratios of E J /E C and the red lines correspond to contours of constant qubit frequency ω q 2 .c Schematic representation of the circuit used to infer the weights in Eqs. 25 and 26, with the state preparation pictured on the left (in purple), single-qubit π and π/2-rotations in green and blue respectively, and the two-qubit CZ gate in orange.N is the number of times the pictured circuit (without the state preparation) is repeated before measurement and therefore an integer determined by the reference coherence time.In our case, we consider N = T 1,0 10(t TQG +2t SQG ) .d The function 1 − P defined in Eq. 26 plotted for different values of the second (higher frequency) transmon E J and E C .We consider a single-qubit gate implemented with a Gaussian DRAG pulse with a duration of t SQG = 16 ns, a two-qubit gate duration of t TQG = 50 ns and three different reference coherence times indicated on top of each panel.For all three cases, the reference T φ ,0 = T 1,0 at E J = 12 • h GHz and E C = 0.2 • h GHz for Qubit 2 (q 2 ) with the parameters of Qubit 1 (q 1 ) given in Table 1.All the parameters (transmon parameters and decay times) are scaled accordingly to different values of E J and E C , and for each qubit individually, as described in Tables 1 and 2. The striped green area marks the region in the plot with the lowest values of 1 − P, defined by the 10th percentile of the plotted values.The darker green contour is obtained with a density matrix simulation of the circuit from panel c with the same errors, but instead of evaluating the function 1 − P, it is obtained by minimizing the infidelity of the state before measurement.The lighter green dotted contour is the optimal region obtained using the advanced noise model from Table 2.
In the following, we consider two noise models.It has been demonstrated that the gate infidelity due to decoherence is, up to first order, independent of the unitary dynamics 84 .Therefore, in the basic noise model, presented in Table 2, both T 1 and T 2 decay are modeled by applying the corresponding Kraus operators after the gate operation.We further assume that the detuning between the computational transmons is consistently ω q 1 − ω q 2 ≈ α q 2 , and the anharmonicity of the qubits is similar.Given these conditions, we parameterize the entire system in terms of E J and E C of the higher-frequency transmon.Examples of parameter values, obtained from experimental demonstrations of CZ gates, are plotted in Fig. 4b.
In the advanced noise model (see Table 2), we account for the population leaving the computational subspace during the operation of the diabatic CZ gate 85 .Additionally, in the advanced noise model we also consider the qubits to be detuned by α q 1 /2 ≈ α q 2 /2 during the execution of single-qubit gates (in order to prevent driving unwanted leakage transitions) 86 .The low-frequency qubit is then tuned so that ω q 1 − ω q 2 = α q 2 only during the two-qubit gate operation, thus increasing the flux-noise susceptibility during a CZ gate, according to Eq. 47.
In order to quantify the performance of an algorithm execution with a specified pair of parameters E J and E C in mind, we define a performance metric P which we will then maximize.We further define P as one minus a weighted sum of the  2 (top) and the value of the performance metric P as a function of coherence times with assuming T φ ,0 = T 1,0 at the reference point (bottom).
Here we considered the same reference point as used in Fig. 4d. a A well-designed system, where the parameters E J and E C are optimally adapted so that P is maximized for each value of the coherence time on the x-axis.b P at various coherence times while keeping E J and E C fixed, corresponding to the values marked as point [5] in Fig. 4b.infidelity contributions of all the relevant errors listed in Table 2.This sum can be written as where the summation runs across all infidelity contributions, or more explicitly for the simple algorithm pictured in Fig. 4c and the basic noise model from Table 2 1 The first term in Eq. 25 and first line in Eq. 26 correspond to the errors of the two-qubit gate, the second term and line to single-qubit gate errors and the last to the state preparation error.The performance metric for the advanced noise model can be derived from Table 2.A similar fidelity approximation was defined in Refs. 7,87 We have introduced additional weights in the sum, in order to account for the relative number of single and two-qubit gates, and also to correctly take into account the fact that the error in state preparation occurs only once, while the gate error is significantly amplified after a number of applications of the operation.As the determination of the position of maximal P within the (E C , E J ) landscape, using which we find the optimal range for parameters, depends exclusively on the relative values of the weights w i, j , we proceed by assuming that all error terms associated with the two-qubit gate are assigned weights of w TQG = 1.As pictured in Fig. 4c, we analyze a circuit where we perform four single-qubit gates per each two-qubit gate, with half those single-qubit gates being π rotations which are more susceptible to leakage.Note that this ratio of single to two-qubit gates arises naturally with the introduction of randomized compiling into the algorithm 69 .The weights w SQG,T 1,φ for the decoherence induced infidelity during a single-qubit gate therefore have a value of 2 (since a π and π/2 rotation are applied) per qubit, while the leakage error w SQG,leak.= 1 per each qubit.Lastly, the thermal state preparation error is weighted so that w SP = 10(t TQG + 2t SQG )/T 1,0 , since this quantity is the inverse of the approximate number of single and two-qubit gates that we can perform within a specified coherence time T 1,0 at the reference point.Therefore, the performance metric P approximates the fidelity of the simple circuit displayed in Fig. 4c.
The values of 1 − P from Eq. 26 for certain parameters are plotted in Fig. 4d.Due to the fact that the maximum of P is an optimal solution only for the circuit pictured in Fig. 4c, we have plotted an optimal region of parameters defined by the 10th percentile of the points with the lowest value of 1 − P.This region is defined so that the parameters of the system can be further fine-tuned within the optimal area.This optimal region is also contrasted with the optimal region obtained with a density matrix simulation of the same circuit, and the good agreement between the curves shows that P is a valid performance metric.Figure 6.General quantum processor parameter optimization procedure.The procedure is composed of two parts: the preparatory phase and the actual optimization.The optimization loop is composed of first obtaining the necessary data to evaluate the performance metric from Eq. 25 (step 3), after which we can determine whether the value of P is good enough or whether the optimization should continue.In the latter case, the performance metric is extrapolated to different values of the quantum processor parameters (from step 1.b) we wish to tune as plotted in Fig. 4d and a better parameter set is identified in step 5, based on the results from step 4. Afterwards, a new device is fabricated with the improved parameter set and the cycle repeats.
We also observe good alignment between the optimal regions of the two distinct considered noise models, with the second model anticipated to offer greater accuracy.The discrepancy is larger at lower coherence times, since both models differ mainly in the decoherence modeling of the diabatic CZ gate.Further numerical tests have shown that the more accurate infidelity formulas are the main reason for the difference 85 , rather than the increased flux-sensitivity of Qubit 1 during two-qubit gates, or the differing transmon idling frequency configurations.The color scale corresponding to the 1 − P has four distinct less-favorable regions, corresponding to the error sources in Table 2.For small E C and, consequently, low anharmonicity α q 2 , the overall fidelities are low due to relatively large errors caused by leakage during single-qubit gates.For small E J and large E C , E J /E C has a low value and therefore the system experiences relatively large errors arising from parity switching events.If both E J and E C are large, the coherence times are short and, therefore, the dominating source of infidelity.Moving perpendicularly to the constant frequency contours towards lower values of E C and E J , results in lower frequencies and, thus increased errors caused by thermal excitations.
On the other hand, when considering near-future transmons with coherence times close to 0.5 ms 48,49 in Fig. 4d, we observe that the optimal region is shifted towards larger values of E C , and is even limited by the parity-switching effects in the bottom right corner, thus demonstrating the importance of this effect in future QPU design, when two-qubit gate infidelities surpass the limit of 10 −3 infidelity.
By comparing the data in Fig. 4b and d, we observe that the implementation by Xu, H., et al. (2021) 37 is close to the parity-switching induced error region.Since the pure dephasing times of the system were not reported we are not able to assess the relative contribution of the parity switching error to this gate implementation.However, using the results shown in Eqs.11 and 17, we observe that the parity-switching-error induced infidelity for the parameters is estimated to be 1 − F ≈ 2 • 10 −3 .As a comparison, the reported Γ 1 decay rates (measured in the idling configuration) in the reference were 1/Γ q 1 1 = 20.8µs and 1/Γ q 2 1 = 28.8µs 37 .Together with the reported gate time, these decay times correspond to an infidelity contribution of approximately 6 • 10 −4 , as given by the expression in Table 2.These results show that for this specific implementation of the diabatic CZ gate, the described effects of the parity switches are comparable, and possibly even greater than the T 1 decay induced infidelity.
The shifts of the optimal region with increasing coherence times seen in Fig. 4d mean that the qubit parameters must be adapted to the currently achievable coherence times.Fig. 5a shows how the performance metric P increases with the achievable coherence properties of the system, provided that the parameters E J and E C are adapted to the coherence times.On the other hand, Fig. 5b displays the value of the function P without changing E J and E C .Although the performance metric P on panels a and b is initially largely coherence limited, the assumption of fixed parameters E J and E C in the simulations of Fig. 5b shows that simply increasing the coherence times does not necessarily give better fidelities if the effects of other sources of error are not taken into account.In this case, the parity switching error becomes the dominant error source, as seen from the upper panel in Fig. 5b.
It is also important to mention what other aspects of successfully operating a transmon tunable coupler based quantum 14/27 computer were not included in the presented analysis.
• We are neglecting any possible cross-talk effects between next-nearest neighbors.
• Idling errors were discussed in this work, but not included in this analysis, as their contribution depends significantly on the specific algorithm being implemented.Taking also this effect into account, the area with low E J /E C is even less favorable.
• Other sources of decoherence are expected to have smaller contributions and therefore do not significantly affect the findings presented here.
• The TLS environment is random, meaning that the scaling shown in Table 2 only holds for the average of a large number of qubits.Additionally, the T 1 decay rate is also heavily influenced by the design, i.e., the exact geometry of the capacitor pads of the transmon, meaning that the presented results are only valid for a comparison of qubits with similar designs.
• Recently, two-qubit gates with tunable couplers idled below the qubit frequencies have been demonstrated 35,41 .In such implementations, the presence of thermal coupler excitations limit the gate performance at lower qubit frequencies in a more complex manner.Our analysis in the low frequency regime therefore only holds for the original tunable coupler proposals, where the coupler is idled above the computational transmons.
While the results in Figs.4d and 5a,b demonstrate how the knowledge about the error contributions can be used to guide the design of superconducting circuits, the scheme can easily be generalized to any type of quantum processor and any set of error sources.The only requirement is that the infidelity contribution of each error source depending on the tunable parameters of the processor can be evaluated either analytically or numerically.Especially the latter is typically feasible for errors which are local, i.e. do not affect a large number of components.
The schematic in Fig. 6 shows how the optimization procedure can be generalized to different quantum processor parameters and errors.The preparatory phase includes first identifying the parameters of the quantum processor we wish to optimize as well as the error sources whose contributions depend on these parameters (step 1).In the transmon example this corresponds to E J and E C of the high-frequency qubit and the list of error sources in Table 2.In the second step an analysis of the infidelity contributions of the selected error sources with respect to the tunable parameters is performed, as described in Table 2.At this point, analytical formulas, such as the one in Eq. 18, are invaluable.However, in certain cases numerical simulations are needed to obtain these relations.Generally speaking, numerical simulations are feasible in the case when the error source is local to a small number of qubits (i.e. the size of the Hilbert system for the simulation is not too large).
Once the preparatory phase is complete, the optimization step of the processor can commence by first obtaining the scaling coefficients needed to evaluate the infidelity contributions of each error (step 3 in Fig. 6), and hence evaluate the performance metric P.This step can include measurements of the average coherence times of the qubits, their effective temperatures, etc.In some cases, such as for the parity switching error described in this manuscript, no additional measurements are required.Once the performance metric is evaluated one can decide whether to proceed with the optimization based on a predefined condition which can be chosen arbitrarily.As an example, in Fig. 4d, the convergence condition is satisfied if the current QPU parameters are within the optimal region defined by the top 10th percentile.If the value of the performance metric was found to be insufficient, a new set of parameters is chosen based on the the landscape of the performance metric (step 5), as shown in Fig. 4d.

Discussion
We have provided a novel framework for the optimization of circuit parameters that can be used to guide the future design of transmon-based quantum computers.Our findings reveal the presence of a distinct global performance peak within the E J and E C parameter space, which has not been identified before.Moreover, our optimization procedure can be straightforwardly extended to more error sources, provided that the scaling of the error as a function of the system parameters is known analytically, or the infidelity contribution can be evaluated numerically.The latter is typically realistic as long as the error is sufficiently local, i.e., it depends only on the parameters of a handful of transmons at most.While we have based our analysis on transmon qubits connected via tunable couplers, the same principles, albeit with different error sources, can be applied to different types of qubits 2 or co-design chips 88 .Additionally, more parameters than just E J and E C can be optimized, e.g., also the gate durations can be considered as free parameters since they are realistically easy to adjust in experiments.The limiting factor here are the errors for which the analytical behavior is unknown and numerical interpolation in a large parameter space is too demanding.
We have chosen to present a range of optimal parameters, rather than prescribing a single optimal value for E J and E C since additional improvements of the noise models are expected to improve the accuracy and reliability of the procedure.An extension of the framework in the tunable coupler example would include also optimizing the E J and E C of the coupler transmon, as well as the couplings between the three transmons of the system.We are not aware of any analytical results linking the gate infidelity with these parameters, however it is known that the coupler coherence affects the effective qubit decoherence due to the hybridization of the states 39,40 .Correctly taking into account the effects of different coupling strengths necessitates a good understanding of the effect of the state hybridization on single-qubit rotations 52,86 and the two-qubit gate duration (Eq.11).We believe that the inclusion of pulse-level numerical modeling, as presented in Ref. 40 , can also take into account such effects.Another straightforward improvement of our framework would be to consider a set of (possibly random) parameters for E J and E C in the evaluation of the performance metric in order to better mimic the spread in transmon frequencies required to avoid unwanted resonances in a more realistic lattice 82 .This would also allow us to better estimate the effects of flux noise.
It is important to acknowledge the potential influence of other noise mechanisms, not explicitly addressed in this study, on the optimal design parameters.Specifically for superconducting circuits, the list of error sources in Table 2 can be expanded to include more contributions to the circuit coherence times, such as losses due to charge-coupling to an impedance, current fluctuations in the flux-bias lines, etc.All of these contributions to the coherence times can be evaluated using open-source software such as scqubits 89 .We believe that even taking into account cross-talk effects is possible, using similar approaches as in Ref. 87 .Nonetheless, the results in Fig. 4d show that the optimal parameter regimes presented are robust to small perturbations in the idling configuration, and higher-order corrections to the noise models.Our findings, as depicted in Fig. 4d, reveal that a two-fold variation in the reference coherence time only marginally adjusts the optimal parameter domain.This observation underscores the robustness of our results, suggesting that the presence of additional, potentially sub-leading noise mechanisms omitted from our simulations are unlikely to precipitate a drastic alteration in the presented outcomes.
Additionally, we have established how parity switching affects the commonly implemented tunable-coupler mediated diabatic and adiabatic CZ gate in a transmon based quantum computer, both analytically and numerically.We have shown that the parity switching error can be the main quasiparticle-related error source of the two-qubit gate.Moreover, we have demonstrated that the experimental implementation of the gate presented in Ref. 37 may have a comparable, if not larger, contribution of parity-switching errors compared to all T 1 decay mechanisms.
While the tunable-coupler-based diabatic and adiabatic CPHASE gates are more relevant due to its implementation in leading large-scale experiments 7,8,10,11 , we believe that the effects described in this manuscript should be considered in any current or future transmon-based quantum gate which utilizes higher-excited states, such as the gate schemes proposed in Refs. 34,35,90 .Accordingly, we do not believe that the parity-switching effects play a major role in implementing iSWAP-like interactions in the tunable coupler architecture, since the states involved in the interaction are not significantly hybridized with the second excited state of any transmon.Furthermore, the effects on parity-switches on more general fSIM gates 91 are expected to be of a similar magnitude as in the adiabatic CZ case presented in Fig. 3c.
One of the primary anticipated advantages of incorporating tunable couplers into the system is the potential for on-demand complete suppression of ZZ-type interactions among the qubits.However, our research in this context has revealed that the stochastic nature of parity switches imposes constraints on this proposition, practically establishing a lower bound on the achievable minimum ZZ coupling strength.The magnitude of this "always on" interaction should be an important consideration when running longer algorithms.More specifically, this effect becomes relevant if the algorithm is long enough to accumulate a considerable conditional phase due to the unwanted coupling strengths shown in Fig. 3b.Since current coherence times are approaching the 1 ms limit 48,49 , the residual idling strength can become relevant if the described effects are not taken into account in the design of the transmon parameters.
One way of mitigating the parity-switching effects would be to attempt to tune the offset charge n g to the point where both parity manifolds are degenerate.However, such a solution is not practical, since the environmental charge noise would result in a drift of the offset charge n g as was demonstrated in Refs. 92,93 We further note that cos(2πn g ) = 0 is the low-frequency charge noise hotspot 12 in which the qubit frequency is maximally sensitive to the fluctuations of the offset charge.Therefore, the qubit is expected to have lower coherence times at this particular value for n g .

Modelling the Tunable Coupler Circuit and Diabatic CPHASE Gate
The dependence of the transmon parameters on the parity and offset charge is not explicitly shown in Eq. 5 in order to simplify the notation.However, we note that the whole system has 2 3 = 8 distinct parity states.We emphasize also that the couplings between the transmons in Eq. 5, g i j = β i j √ ω i ω j also depend on the frequencies, meaning that while g q 1 q 2 is constant, g q 1,2 c is time dependent.The dimensionless prefactors β i j depend on the coupling capacitances, as well as self-capacitances of the transmons in the lumped-element circuit model 38 .
The computational basis in this scheme is formed by the eigenstates of the Hamiltonian from Eq. 5 in the idling configuration (defined below), rather than the local (uncoupled) transmon states 39 .Since the couplings act only as a perturbation to the uncoupled states, we identify the full Hamiltonian eigenstates corresponding to the uncoupled states.More specifically, the computational state |i j⟩, i, j ∈ {0, 1} is the eigenstate |ψ⟩ of the Hamiltonian in Eq. 5 with the maximal overlap |⟨ψ|i q 1 0 c j q 2 ⟩|.This notation is employed throughout this manuscript and we typically omit the subscript indices q 1,2 and c.The kets with three indices (e.g.|i jk⟩) always denote the local (uncoupled) Fock states of the three-transmon system.The kets with only two indices (e.g.|i j⟩) are used to denote the eigenstates of the whole system that are closest to the local (uncoupled) state |i0 j⟩.We denote with ω i j the angular frequency of the computational state |i j⟩.These eigenenergies of these states are also used to computed the ZZ interaction strength ζ ZZ , as defined in Eq. 20.
However, since the coupler frequency is tunable, it is possible to find one or two frequencies ω c for which ζ ZZ = 0 36,38,41 , provided that the qubit-qubit detuning ω q 1 − ω q 2 ∈ [α q 2 , −α q 1 ].These special frequencies are referred to as the coupler idling frequencies and denoted with ω idle c .A variety of pulse shapes can be used to implement the gate effectively and without inducing too many unwanted transitions 36,39,62 .Even though our analytical analysis makes minimal assumptions about the pulse shape, we need to choose a specific form for the numerical simulations.In our simulations we use the flattop Gaussian pulse described by the formula 36,46 The flattop Gaussian is obtained by a convolution of a step function with duration τ c and a Gaussian with parameter σ .The reasoning behind this choice is that the convolution of the flattop pulse with a Gaussian strongly suppresses the spectral component of the flattop pulse at higher frequencies, thus reducing the probability for unwanted transitions.An additional rise time of τ b is also introduced, which we fix to τ b = 2 √ 2σ .Since any gate must have a finite duration, we introduce a cut-off at time T = 2τ b + τ c .The constant C is then introduced to correct for the discontinuity at the beginning and the end of a pulse with finite duration.The coupler frequency is varied accordingly ω c (t) = ω idle c − A f (t).More details about the numerical simulations are given in Ref. 40 .The parity of the whole system is then switched from one state to the other according to the charge dispersion given by Eqs. 2 and 3.
Since in our analytical approach we approximate the realistic pulse shape shown in Eq. 27 with a square pulse (i.e., we neglect any dynamics during the ramping up and down of the flux pulse), we have to be careful when comparing the results to those obtained with the numerical data.Particularly, the pulse duration T in the numerics is in general different than the gate duration t g we have defined in our analytical derivation, since the full duration T also includes the time needed to ramp the pulse up and down.When we compare our analytical results with numerical simulations, we first extract the number n of Rabi oscillations from the simulated data by monitoring the population of the |002⟩ state during the gate operation.After n is determined, the analytical gate time is adjusted such that t g = nπ/Ω.The effective gate time obtained in this manner does not differ significantly from the duration of the flat part of the pulse, typically less than 5 ns.

The Schrieffer-Wolff Transformation
The Schrieffer-Wolff (SW) transformation used in this manuscript was first introduced in Ref. 94 and a similar transformation has been applied to the computational subspace of the two-qubit system in Refs. 38,52 The aim of the SW transformation in our case is not to diagonalize the system, but rather to decouple the coupler states from the computational transmons, thus enabling us to only study the reduced system, i.e. we want to simplify the full Hamiltonian into a more tractable reduced model, containing only the relevant states (i.e. the states which have a significant population).
The reason for this is the fact that the computational basis of the system is defined by the eigenstates of the system which have the maximum overlap with states of the form |i q 1 0 c j q 2 ⟩, i.e. with the coupler always in the ground state.Any excitations of the coupler therefore lead to errors, so in order to analyze the ideal gate dynamics, we constrain ourselves only to the energy levels of the computational subspace and the second excited state used in the Rabi oscillation during the CPHASE gate operation.
In general, a SW transformation is obtained by noting that any unitary operator can be written as Û = e Ŝ = 1 where Ŝ is anti-hermitian, Ŝ = − Ŝ † .Consequently a unitary transformation of an arbitrary Hamiltonian Ĥ can be expanded in terms of Ŝ as As is typical in perturbation theory, we introduce the index α for bookkeeping purposes, and split the full Hamiltonian into a diagonal part and two off-diagonal perturbations, so that Ĥ = Ĥ0 + α V1 + α 2 V2 .Additionally we rewrite the operator Ŝ as a first order operator Ŝ → α Ŝ, since if the perturbation is small, Û should be close to identity.In our case, Ĥ is given in Eq. 5 and the first-order perturbation V1 = − ∑ i=q 1 ,q 2 hg ic ( â † i − âi )( â † c − âc ) corresponds to the capacitive couplings of the two transmons (Qubits 1 and 2, q 1,2 ) to the coupler (c), the direct coupling between the qubits V2 = −hg q 1 q 2 ( â † q 1 − âq 1 )( â † q 2 − âq 2 ) is a second order perturbation, while Ĥ0 is a sum of the three independent anharmonic oscillator Hamiltonians.This hierarchy is chosen due to the fact that in all practical scenarios g q 1 q 2 ≪ g q 1 c , g q 2 c 38 .
By plugging the ansaetze into Eq.28, and grouping the terms with the same order of α, we obtain Looking at the first-order term, it is natural to choose Ŝ such that Ŝ, Ĥ0 = − V1 , i.e. so that we cancel any couplings to the coupler states up to lowest order.However in order to do so and account for the couplings to the higher state correctly, we generalize the transformation from Refs. 38,52 = ∑ n q i ,n c ∈{0,1} (n q i + 1)(n c + 1) We have additionally defined ∆ q i c = ω q i − ω c , Σ q i c = ω q i + ω c and the operators πn,m k = |n⟩⟨m|, acting in the Hilbert space of k ∈ {q 1 , c, q 2 }.
Since we have assumed the coupler remains in the ground state at all times, the effective Hamiltonian is defined on the set of states {|0 q 1 0 c 0 q 2 ⟩, |0 q 1 0 c 1 q 2 ⟩, |1 q 1 0 c 0 q 2 ⟩, |1 q 1 0 c 1 q 2 ⟩, |0 q 1 0 c 2 q 2 ⟩}.Additionally, we neglect any couplings outside of the subspace of interest, however the resulting effective Hamiltonian still contains terms coupling the levels |000⟩ ↔ |101⟩ and |000⟩ ↔ |002⟩.These couplings are neglected in the rotating-wave approximation, as these transitions do not conserve the total occupation number.
By additionally setting the energy of the ground state to zero we arrive at the effective subspace Hamiltonian from Eq. 36.The perturbative parameter values are given by g01,10 = g q 1 q 2 + g q 1 c g q 2 c 2 1 g11,02 = √ 2 g q 1 q 2 + g q 1 c g q 2 c 2 1

Deriving the Effective Diabatic Gate Unitary
We consider here the truncated Hilbert space spanned by the states {|00⟩, |01⟩, |10⟩, |11⟩, |02⟩}.In our analytical considerations, we thus neglect the possibility of leakage, which is a good assumption for high-fidelity gates.We emphasize that the leakage effects are included in our numerical simulations which have been made with the full Hamiltonian defined in Eq. 5.
We additionally neglect the counter-rotating terms within the truncated Hilbert space and, after applying the Schrieffer-Wolff trnasofrmation to the Hamiltonian in Eq. 5, we obtain where the parameters with the tilde denote the perturbed parameters of the original full Hamiltonian from Eq. 5 (see Eqs. 32-35)  up to second order in g q i c /(ω q i − ω c ).Without any loss of generality, we have also assumed that the second excited state of the second qubit (q 2 ) is used to perform the Rabi oscillation.The pulse shape (see Eq. 27) used to implement a diabatic CPHASE gate can typically be divided into three distinct stages: (i) fast sweep of the coupler frequency to an operation value close to the resonance with the qubits; (ii) long constant-frequency plateau at the operation frequency; (iii) fast sweep back to idling frequency of the coupler.Since the plateau is typically longer compared to the fast sweeps 36 , and the effective coupling strengths g01,10 and g11,02 are significantly smaller if the coupler is further detuned, we assume that the plateau stage is the only contributor to the dynamics and the effective Hamiltonian from Eq. 36 is constant in time.We show later that this approximation is valid for realistic gate durations close to 40 ns and longer.If this approximation cannot be made, the form of the effective Hamiltonian remains unchanged.However, the analytical results derived in the remainder of this section become more complex due to the time dependence of the effective Hamiltonian.
The effective Hamiltonian in Eq. 36 is block-diagonal.Using matrix exponentiation, one can readily obtain an effective unitary time-evolution operator.Assuming that the computational transmons are detuned by approximately one anharmonicity, i.e. ωq 1 − ωq 2 ≈ αq 2 , the amplitude of the Rabi oscillation between the single-excitation states |01⟩ and |10⟩ is lower than that of the transition between |11⟩ and |02⟩ because | ωq 1 − ωq 2 | ≫ g01,10 .Thus, we can neglect the interaction between the single-excitation states.
Since we are interested in the operator acting on the computational subspace of the system, we further truncate the subspace by excluding the non-computational |02⟩ state.After also accounting for the single-qubit phases, which is typically done via virtual Z-rotations 95 , we obtain the effective time-evolution operator Û(t) from Eq. 6. Eq. 7 for P 11 is derived by recognizing the block-diagonal structure of the effective Hamiltonian in Eq. 36.This simplifies the analysis, reducing it to a standard two-level Rabi oscillation 63 .Eq. 8 is obtained by subtracting the phases of the single-excitation states from the phase of the |11⟩ state.

List of Assumptions
Here we summarize all of the approximations involved in deriving Eqs.11 and 12 in the main text.These expressions are valid under the following assumptions: 1.The initial assumptions used to derive the effective Hamiltonian in Eq. 36 are valid, meaning that: (a) The second order perturbation theory used for the Schrieffer-Wolff transformation is valid, i.e. g 2 q 2 c /(ω q 2 − ω c ) 2 ≪ 1.
(b) The gate has low leakage outside of the considered subspace, spanned by the computational states and second excited state of the relevant qubit.
(c) The rotating wave approximation for the counter-rotating coupling terms is justified.
3. The coupler frequency in the interaction regime is relatively constant.If this is not the case, the time dependency of the perturbative parameters must be taken into account.
4. The gate has low leakage outside of the computational subspace.
5. The simplified formula in Eq. 11 additionally neglects the terms proportional to g q 1 c g q 2 c /(∆ q 2 c + α q 2 ) 3 and g 2 q 2 c /(∆ q 2 c + α q 2 ) 3 , and smaller.6.We have assumed that the main contributor to the perturbation is the Qubit 2 (q 2 ) whose second excited state is populated during the gate operation.However, if the system is designed in such a way that the charge dispersion of any of the other two transmons is significantly larger, their effects might not be negligible anymore.
Note that both assumptions about the leakage are automatically fulfilled if the gate has a high-fidelity.

Fidelity Scaling
Scaling of the fidelity with the number of gates, i.e. computing the average gate fidelity of sequential application of N gates, can be found by using Eq.14 to first define the map corresponding to a series of gates where we have used the fact that the operators in Eq. 15 are diagonal and therefore commute with each other and also Û † + = Û− .Using the effective Kraus operators from Eq. 37, and only considering the conditional phase error, we combine this with the 19/27 fidelity definition from Eq. 16, and obtain This result indicates that calibrating the gate such that the error is evenly split between the two parities not only increases the single-gate fidelity but also leads to a more generous scaling of the infidelity ∝ N, compared to the purely coherent error case for which the error scales as ∝ N 2 69 .

Gate parameters for high-fidelity simulations
Finding good gate parameters, both for the Hamiltonian as well as for the pulse for high-fidelity simulations, is not a trivial task.
Here, we discuss how to find optimal Hamiltonian parameters at different qubit frequencies and anharamonicities.In general, the only prerequisite is that the qubits are detuned by approximately one anharmonicity, which can be seen from Eq. 36.
As seen from Table 1, we keep some of the parameters in the simulation fixed, while others depend on α q 2 , which is varied.These parameters are based on the experimental values from Ref. 36 and the coupling coefficients from the Hamiltonian in Eq. 5 are related to β i j via the following relation g i j = β i j √ ω i ω j , as described in the main text.The small perturbation to α q 1 is there in order to ensure that the energy levels are significantly non-degenerate for perturbation theory to apply.This is also completely realistic as a typical fabrication procedure results in seemingly random deviations from the designed values.More specifically, the gates simulated in Fig. 2(a,b) are obtained by varying α q 2 with values α q 2 ∈ [195, 230, 250, 270, 300] h MHz.Additionally, in order to certify that the analyzed effect is not limited to the choice of parameters presented above, the second data point of Fig. 2(a,b) at E J /E C ≈ 50, was generated in the same way as in Table 1, but with the change to ω q 2 = 5.1 • 2πGHz.The coupler frequency in the idling configuration ω idle c , i.e. before and after performing a gate, is determined by diagonalizing the Hamiltonian to fulfill the condition ζ ZZ = 0 from Eq. 20.The pulse parameters from Eq. 27 are obtained by numerically optimizing the fidelity of the gate, with fixed σ = 5 ns and τ b = 2 √ 2σ .Typical values of the amplitude are A ∼ 1 − 1.2 • 2π GHz and τ c ∼ 60 ns.

Scaling of Transmon Error Sources
In this section we derive the scaling of the noise parameters shown in Table 2 with the transmon energies E J and E C .The noise is modeled by appending the appropriate noise channel after the gate unitary, and the calculated infidelity is thus independent of the unitary dynamics, as was shown in Ref. 84 .

Charge noise T 1
The T 1 decay time of transmon devices is believed to be currently limited by the presence of a discrete number of environmental two-level systems (TLSs) which couple to the qubit via their electrical dipole [71][72][73][74][75][76][77] .Therefore assuming that the charge noise in the system is weak enough, the interaction Hamiltonian between a transmon and an environmental TLS can be derived from Eq. 1, by replacing n g → n g + δ ng , as in Refs. 72,77  so that In the above equation, we have used the asymptotic expression for the number operator n, derived already in Ref. 12 .The operator δ ng is defined in the Hilbert space of the TLS, and is related to the parameters of the TLS -more specifically, its electrical dipole 75 .More importantly, δ ng does not explicitly depend on E J and E C .Since the majority of the T 1 experiments on transmons display exponential decays 96 , the dynamics can be captured by the Lindblad equation.As a second order approximation is assumed in the derivation of the Lindblad equation, the resulting decay rates are found to be proportional to the square of the coupling coefficient of the transmon to the TLS environment 66 .In our case this translates to

20/27
Here we have omitted the noise spectrum of the environment, since current models assume a flat noise spectrum without any dependence on the qubit frequency 75 .

Flux noise T φ
While it is not strictly necessary for the computational transmons to be flux-tunable, it is often desired as flux-tunability additionally enables the implementation of an iSWAP gate with the same architecture 36,38 .Having slow flux tunability is also desirable in order to avoid resonances with TLSs in the environment which can severely limit the T 1 decay time 76,77 , as well as helping with the issue of frequency crowding 82 .However, unlike the random noise spectrum of a TLS environment producing T 1 dynamics, the noise spectrum of magneticflux noise is typically observed to have a 1/ f α frequency dependence 71,[78][79][80] , with α ∼ 1.The large noise-spectrum amplitude at lower frequencies means that the long-time correlation results in non-Markovian dynamics 97 .Assuming that the noise is slow enough, in order for the adiabatic approximation to hold, the interaction Hamiltonian due to a slowly fluctuating magnetic environment δ Φ is given by 12 Ĥq−flux = ∂ Ĥ Similarly to the T 1 -decay scenario, we make the assumption that the environment operator δ Φ remains independent of the transmon parameters.In reality, the magnitude of this operator is contingent upon the inner product between the magnetic dipole operator of the spins and the surface vector of the SQUID loop.The form of the flux dispersion ∂ ω/∂ Φ of a split-junction transmon, with Josephson energies E J 1 and E J 2 , is determined by the relations 12 E J (Φ) = E JΣ cos(πΦ) 1 + d 2 tan 2 (πΦ), with E JΣ = E J 1 + E J 2 and d = |E J 1 − E J 2 |/(E J 1 + E J 2 ).We assume that d ≲ 1, since values of d ≈ 1 result in overall longer coherence times, and the computational transmons need to be tuned over a smaller frequency range, typically on the order of ∼ 100 MHz 36 .In this regime, the flux dispersion can be approximated as and the decay rate scales with the coupling coefficient 98,99 , Ideally, the system is designed so that both the idling frequency of the computational qubits (i.e. the frequency at which single-qubit rotations are performed), as well as the computational transmon frequency during the two-qubit gate operation are positioned at the first-order flux-insensitive sweet-spot.In that case, it can be seen from Eq. 47, that the flux noise decay rate is identical in both the idling configuration as during the gate operation.This is true for both computational qubits.While we have explicitly addressed only the first-order flux dispersion ∂ ω/∂ Φ, it is evident that the aforementioned conclusion remains valid even if ∂ ω/∂ Φ = 0.In such cases, we must consider the second order dispersion ∂ 2 ω/∂ Φ 2 .A simplified scaling relationship of Γ φ with E C and E J can be derived under the condition that the applied flux Φ is independent of E J and E C , or equivalently, when we are interested in the flux-averaged decay rate Γφ = 1 0 dΦ, Γφ (Φ).This yields the relation: which holds for all values of d.
The averaging over the flux is relevant, in a more realistic scenario, where the frequencies of the system may be calibrated in order to avoid resonances with individual TLSs and other elements on the chip 87 .Due to fabrication inaccuracies 83 and the stochastic nature of the TLSs, it is impossible to deterministically predict the optimal frequencies in an array of transmons, which prompts us to consider the flux-averaged decay rate scaling in Eq. 49.
Realistic gate implementations frequently involve fast flux tuning of one of the computational transmons to enhance gate fidelity 36,91 .This flux tuning influences the flux dispersion of the qubit and consequently impacts the decay rate during the two-qubit gate implementation.While predicting the exact detunings of all qubits proves challenging, in a more realistic noise model, we incorporate the flux dependence of the dephasing rate from Eq. 49 to better model the dephasing during the operation of one of the transmons in the two-qubit gate.Typically, tuning the low-frequency transmon (q 1 ) proves more advantageous, as the second excited state of Qubit 2 is inherently more susceptible to decoherence.
In the advanced noise model, we therefore consider that the frequency of Qubit 1 satisfies ω q 1 − ω q 2 = α q 2 /2 ≈ α q 1 /2 for all values of E J and E C during single-qubit gate operation.This will reduce any potential leakage due to driving the |10⟩ ↔ |02⟩ transition when instead of |10⟩ ↔ |11⟩ 86 .More specifically we consider Qubit 1 to be weakly tunable (d = 0.9) and we assume its frequency is shifted so that the two-qubit gate resonance condition is satisfied, i.e. ω q 1 − ω q 2 = α q 2 during the two-qubit gate operation.Otherwise the transmons are assumed to be idling 10 MHz away from the sweet-spot.This information is sufficient to compute the decay rate of Qubit 1 at a different flux bias, by using Eq.47.We believe that pinpointing more specific detuning values would require pulse-level simulation.

Single-Qubit Gate Leakage
The low anharmonicity of the transmon is a limiting factor in the operation of single-qubit gates 40,81 as fast operations drive a part of the population from the computational subspace into the second excited state.A straightforward and effective scheme for mitigating this effect while implementing single-qubit gates known as Derivative Removal by Adiabatic Gate (DRAG) was presented in Ref. 81 .
Since typical pulse amplitudes used to perform single-qubit operations are typically much lower compared to the qubit frequency, and the drive is resonant with the qubit, we assume that the rotating-wave approximation is accurate.In the frame rotating with the qubit frequency, the effective Hamiltonian depends on the pulse parameters, anharmonicity, and the detuning between the drive and qubit frequencies.This indicates that the amount of leakage does not explicitly depend on the qubit frequency.
Thus, the average gate fidelity, similar to the one defined in Eq. 16, but in this case for a single-qubit rotation around the x or y-axes, depends only on the qubit anharmonicity (or charging energy), since α ≃ −E C in the transmon limit.This means that even though analytical results are not available, the relationship can be determined numerically and interpolated.This relationship, albeit with different parameters, has already been plotted in Ref. 81 , and generally follows a dependence of P leak.∝ E −γ C , 5 ≲ γ ≲ 6, with higher exponents observed at lower E C .The independence of the single-qubit gate infidelity of the qubit frequency (within the transmon regime) was also verified numerically.
The single-qubit gate parameters assumed in Fig. 4d and Fig 5 were a DRAG Gaussian pulse, with a σ = 4 ns and a total duration t SQG = 4σ .The amplitudes of both DRAG components are numerically optimized before interpolating the dependence of P leak.on α, which was used to generate Fig. 4d.The pulse drive frequency is assumed to be resonant with the qubit.More details are available in Ref. 40 .We note here, that most of the SQG infidelity is due to leakage, rather than phase errors.

Thermal Excitation Error
While other error sources affect the gate performance, the thermal-excitation error considered here only affects the state preparation.The average gate fidelity, which is defined as the fidelity averaged over Haar random distributed input states, is therefore not applicable, since this error only affects one input state.We therefore replace the average gate infidelity 67 with a state infidelity.This can also be qualitatively thought of as replacing the Haar random distribution with a delta-like distribution with a peak at the |0⟩ state.By modeling the thermal excitation as a bit-flip channel E [ ρ] = (1 − P |1⟩ ) ρ + P |1⟩ σx ρ σx it is straightforward to see that the state infidelity 63 where we have additionally assumed that the temperature is low enough such that the population of the higher-excited states is negligible.
Tables Table 1.Table of Hamiltonian parameters (Eq.5) with which high-fidelity CZ gates are possible.
4.8 • h GHz ω q 1 ω q 2 + α q 2 + 10 • h MHz β q 1 c 0.015 β q 2 c 0.015 β q 1 q 2 0.001 4 = " M 6 E H L w P 5 Y c f M B 1 4 o f 7 X a 1 b n 8 e 1 8 = " > A A A B 8 n i c b V D L S g N B E O y N r x h f U Y 9 e B o P g x b A r v o 5 B L x 4 j m A d s 1 j A 7 m U 2 G z M 4 s M 7 N C W P M Z X j w o 4 t W v 8 e b f O E n 2 o I k F D U V V N 9 1 d Y c K Z N q 7 7 7 R S W l l d W 1 4 r r p Y 3 N r e 2 d 8 u 5 e U 8 t U E d o g k kv V D r G m n A n a M M x w 2 k 4 U x X H I a S s c 3 k z 8 1 i N V m k l x b 0 Y J D W L c F y x i B B s r + U / e w 0 l H Y d H n t F u u u F V 3 C r R I v J x U I E e 9 W / 7 q 9 C R J Y y o M 4 V h r 3 3 M T E 2 R Y G U Y 4 H Z c 6 q a Y J J k P c p 7 6 l A s d U B 9 n 0 5 D E 6 s k o P R V L Z E g Z N 1 d 8 T G Y 6 1 H s W h 7 Y y x G e h 5 b y L + 5 / m p i a 6 C j I k k N V S Q 2 a I o 5 c h I N P k f 9 Z i i x P C R J Z g o Z m 9 F Z I A V J s a m V L I h e P M v L 5 L m a d W 7 q J 7 f n V V q 1 3 k c R T i A Q z g G D y 6 h B r d Q h w Y Q k P A M r / D m G O f F e X c + Z q 0 F J 5 / Z h z 9 w P n 8 A 3 B 6 Q / w = = < / l a t e x i t > |1 i < l a t e x i t s h a 1 _ b a s e 6 4 = " v n U I s Z Y g p 1 H 6 l S N k f b Q n h 7 a K B 9 E = " > A A A B 8 H i c b V D L S g N B E O z 1 G e M r 6 t H L Y h A 8 h V 3 x d Q x 6 8 R j B P C R Z w u y k N x k y M 7 v M z A p h z V d 4 8 a C I V z / H m 3 / j J N m D J h Y 0 F F X dd H e F C W f a e N 6 3 s 7 S 8 s r q 2 X t g o b m 5 t 7 + y W 9 v Y b O k 4 V x T q N e a x a I d H I m c S 6 Y Y l a t e x i t s h a 1 _ b a s e 6 4 = " Z d C b w O 5 0 u + 1 2 t s 1 e 7 I d YO P 4 d R k k = " > A A A C B X i c b V D L S g M x F M 3 U V 6 2 v U Z e 6 C B a h g p Q Z 8 b U s u H H Z g n 1 A Z x g y a a Y N T T J j k h H K 0 I 0 b f 8 W N C 0 X c + g / u / B s z b R f a e i D h c M 6 9 3 H t P m D C q t O N 8 W 4 W l 5 Z X V t e J 6 a W N z a 3 v H 3 t 1 r q T i V m D R x z G L Z C Z E i j A r S 1 F Q z 0 k k k Q T x k p B 0 O b 3 K / / U C k o r G 4 0 6 O E + B z 1 B Y 0 o R t p I g X 3 o c a Q H k m e N N K Q a e q f Q H e d / 5 T 5 w T w K 7 7 F S d C e A i c W e k D G a o B / a X 1 4 t x y o n Q m C G l u q 6 T a D 9 D U l P M y L j k p Y o k C A 9 R n 3 Q N F Y g T 5 W e T K 8 b w 2 C g 9 G M X S P K H h R P 3 d k S G u 1 I i H p j L f W c 1 7 u f i f 1 0 1 1 d O 1 5 f b b + V Q P O n q U = " > A A A C B X i c b V D L S g M x F M 3 U V 6 2 v U Z e 6 C B a h g p S Z 4 m t Z c O O y B f u A z j B k 0 r Q N T T J j k h H K 0 I 0 b f 8 W N C 0 X c + g / u / B s z 7 S y 0 e i D h c M 6 9 3 H t P G D O q t O N 8 W Y W l 5 Z X V t e J 6 a W N z a 3 v H 3 t 1 r q y i R m L R w x C L Z D Z E i j A r S 0 l Q z 0 o 0 l Q T x k p B O O r z O / c 0 + k o p G 4 1 Z O Y + B w N B R 1 Q j L S R A v v Q 4 0 i P J E + b S U g 1 9 E 5 h b Z r 9 l b u g d h L Y Z a f q z A D / E j c n Z Z C j E d i f X j / C C S d C Y 4 a U 6r l O r P 0 U S U 0 x I 9 O S l y g S I z x G Q 9 I z V C B O l J / O r p j C Y 6 P 0 4 S C S 5 g k N Z + r P j h R x p S Y 8 N J X Z z m r R y 8 T / v F 6 i B 1 d + S k W c a C L w f N A g Y

Figure 3 .
Figure 3.Comparison of the analytical formulas with numerical results obtained via exact diagonalization in the idling regime.a Largest derivatives in Eq.21 obtained numerically (solid lines) and from Eq. 22 (filled circles) using the parameters from Table1and α q 2 = −270 • h MHz.The x-axis represents the coupler frequency, with the two idling frequencies denoted with black dashed vertical lines.Additional plots of the first-order derivatives at different qubit detunings are available in the Supplementary Information.b Comparison of numerical and analytical results for the root-mean-squared coupling strength at the idling point (defined in Eq. 23), with the bar graph displaying the numerical values of the coupling for the 8 different parity states for the highlighted data point.The parameters of the simulation are chosen identically as in Fig.2and are listed in the Methods section, more specifically in Table1.The parameters for the bar plot correspond to the parameters in panel a. c The susceptibility of an adiabatic CZ (φ 0 = π) gate to a parity switch (defined in Eq. 24) at different coupler frequencies and qubit detunings, measured in GHz (left axis) and in units of α q 2 (right axis), using the same parameters as in panel a.The black contours represent different values of ζ 0 ZZ , as indicated on the plot.The white dashed line indicates the detuning used in panel a.

11 / 27 <
l a t e x i t s h a 1 _ b a s e 6 4 = " r 2 C 9 k / D p P C c S R W a U m u d Z S T g s M R Q = " > A A A B + H i c b V D L S s N A F L 2 p r 1 o f j b p 0 E y y C q 5 K I r 2 X R j c s W 7 A P a U C b T a Tt 0 Z h L m I d T Q L 3 H j Q h G 3 f o o 7 / 8 Z J m 4 W 2 H h g 4 n H M v 9 8 y J E k a V 9 v 1 v p 7 C 2 v r G 5 V d w u 7 e z u 7 Z f d g 8 O W i o 3 E p I l j F s t O h B R h V J C m p p q R T i I J 4 h E j 7 W h y l / n t R y I V j c W D n i Y k 5 G g k 6 J B i p K 3 U d 8 s 9 j v R Y 8 r R h I q r V r O 9 W / K o / h 7 d K g p x U I E e 9 7 3 7 1 B j E 2 n A i N G V K q G / i J D l M k N c W M z E o 9 o 0 i C 8 A S N S N d S g T h R Y T o P P v N O r T L w h r G 0 T 2 h v r v 7 e S B F X a s o j O 5 n F V M t e J v 7 n d Y 0 e 3 o Q p F Y n R R O D F o a F h n o 6 9 r A V v Q C X B m k 0 t Q V h S m 9 X D Y y Q R 1r a r k i 0 h W P 7 y K m m d V 4 O r 6 m X j o l K 7 z e s o w j G c w B k E c A 0 1 u I c 6 N A G D g W d 4 h T f n y X l x 3 p 2 P x W j B y X e O 4 A + c z x 9 1 Q Z O g < / l a t e x i t > Qubits < l a t e x i t s h a 1 _ b a s e 6 4 = " P u s n j p f 7 4 l d d u O A O T O B o R X 3 7 N F Q = " > A A A B + n i c b V D L S s N A F L 3 x W e s r 1 a W b Y B F c l U R 8 L Y v d u K x g H 9 C G M p l O 2 6 H z C D M T p c R + i h s X i r j 1 S 9 z 5 N 0 7 a L L T 1 w M D h 3 H O 5 Z 0 4 U M 6 q N 7 3 8 7 K 6 t r 6 x u b h a 3 i 9 s 7 u 3 r 5 b O m h q m S h M G l g y q d o R 0 o R R Q R q G G k b a s S K I R 4 y 0 o n E t m 7 c e i N J U i n s z i U n I 0 V D Q A c X I W K n n l r o c m Z H i a U 0 m M b P G a c 8 t + x V / B m + Z B D k p Q 4 5 6 z / 3 q 9 i V O O B E G M 6 R 1 J / B j E 6 Z I G Y o Z m R a 7 i S Y x w m M 0 J B 1 L B e J E h + k s + t Q 7 s U r f G 0 h l n z D e T P 2 9 k S K u 9 Y R H 1 p k F 1 Y u z T P x v 1 k n M 4 D p M q Y g T Q w S e H x o k z D P S y 3 r w + l Q R b N j E E o Q V t V k 9 P E I K Y W N L K N o S g s U v L 5 P m W S W 4 r F z c n Z e r N 3 k d B T i C Y z i F A K 6 g C r d Q h w Z g e I R n e I U 3 5 8 l 5 c d 6 d j 7 l 1 x c l 3 D u E P n M 8 f D y u U i Q = = < / l a t e x i t > Couplers < l a t e x i t s h a 1 _ b a s e 6 4 = " M y r 3 K c m i y / f K 9 X Q u u l c h F k n R y g M = " > A A A B 6 H i c b V D L S g N B E O y N r x h f U Y 9 e B o P g K e y K r 2 P W 4 g y r U g Q H C M 7 z C m / P o v D j v z s e 8 N e d k M 4 f w B 8 7 n D 6 n R j N w = < / l a t e x i t > N < l a t e x i t s h a 1 _ b a s e 6 4 = " D E L t b W N a r n R N L O x I H N y K K Q w D 8 m s = " > A A A B 6 n i c b V D L S g N B E O y N r x h f U Y 9 e B o P g K e y K r 2 P Q i 8 e I 5 g H J E m Y n s 8 m Q 2 d l l p l c I S z 7 B i w d F v P p F 3 v w b J 8 k e N L G g o a j q p r s r S K Q w 6 L r f T m F l d W 1 9 o 7 h Z 2 t r e 2 d 0 r 7 x 8 0 TZ x q x h s s l r F u B 9 R w K R R v o E D J 2 4 n m N A o k b w W j 2 6 n f e u L a i F g 9 4 j j h f k Q H S o S C U b T S Q z c R v X L F r b o z k G X i 5 a Q C O e q 9 8 l e 3 H 7 M 0 4 g q Z p M Z 0 P D d B P 6 M a B Z N 8 U u q m h i e U j e i A d y x V N O L G z 2 a n T s i J V f o k j L U t h W S m / p 7 I a G T M O A p s Z 0 R x a B a 9 q f i f 1 0 k x v P Y z o Z I U u W L z R W E q C c Z k + j f p C 8 0 Z y r E l l G l h b y V s S D V l a N M p 2 R C 8 x Z e X S fO s 6 l 1 W L + 7 P K 7 W b P I 4 i H M E x n I I H V 1 C D O 6 h D A x g M 4 B l e 4 c 2 R z o v z 7 n z M W w t O P n M I f + B 8 / g B S / Y 3 X < / l a t e x i t > ⇡ < l a t e x i t s h a 1 _ b a s e 6 4 = " D E L t b W N a r n R N L O x I H N y K K Q w D 8 m s = " > A A A B 6 n i c b V D L S g N B E O y N r x h f U Y 9 e B o P g K e y K r 2 P Q i 8 e I 5 g H J E m Y n s 8 m Q 2 d l l p l c I S z 7 B i w d F v P p F 3 v w b J 8 k e N L G g o a j q p r s r S K Q w 6 L r f T m F l d W 1 9 o 7 h Z 2 t r e 2 d 0 r 7 x 8 0 T Z x q x h s s l r F u B 9 R w K R R v o E D J 2 4 n m N A o k b w W j 2 6 n f e u L a i F g 9 4 j j h f k Q H S o S C U b T S Q z c R v X L F r b o z k G X i 5 a Q C O e q 9 8 l e 3 H 7 M 0 4 g q Z p M Z 0 P D d B P 6 M a B Z N 8 U u q m h i e U j e i A d yx V N O L G z 2 a n T s i J V f o k j L U t h W S m / p 7 I a G T M O A p s Z 0 R x a B a 9 q f i f 1 0 k x v P Y z o Z I U u W L z R W E q C c Z k + j f p C 8 0 Z y r E l l G l h b y V s S D V l a N M p 2 R C 8 x Z e X S fO s 6 l 1 W L + 7 P K 7 W b P I 4 i H M E x n I I H V 1 C D O 6 h D A x g M 4 B l e 4 c 2 R z o v z 7 n z M W w t O P n M I f + B 8 / g B S / Y 3 X < / l a t e x i t > ⇡ < l a t e x i t s h a 1 _ b a s e 6 4 = " Z n v B d M L j C d 5 R N T t 3 E X P Y h Y C G t Z Y = " > A A A B 9 H i c b V D L S s N A F L 2 p r 1 p f V Z d u B o v g q i T F 1 7 L o x m U F + 4 A m l M l 0 0 g 6 d T O L M p F B C v s O N C 0 X c + j H u / B u n a R b a e u D C 4 Z x 7 u f c e P + Z M a d v + t k p r 6 x u b W + X t y s 7 u 3 v 5 B 9 f C o o 6 J E E t o m E Y 9 k z 8 e K c i Z o W z P N a S + W F I c + p 1 1 / c j f 3 u 1 A A A B 9 H i c b V D L S s N A F L 2 p r 1 p f V Z d u B o v g q i T F 1 7 L o x m U F + 4 A m l M l 0 0 g 6 d T O L M p F B C v s O N C 0 X c + j H u / B u n a R b a e u D C 4 Z x 7 u f c e P + Z M a d v + t k p r 6 x u b W + X t y s 7 u 3 v 5 B 9 f C o o 6 J E E t o m E Y 9 k z 8 e K c i Z o W z P N a S + W F I c + p 1 1 / c j f 3 u 1 1 p A I I F n e I U 3 6 8 l 6 s d 6 t j / l o y S p 2 D u A P r M 8 f j V q S Z w = = < / l a t e x i t > t e x i t s h a 1 _ b a s e 6 4 = " W 9 f T K / Z + t a M S 0 / f v S 8 0 i + 1 g E N F 4 = " > A A A C E 3 i c b V D L S s N A F J 3 4 r P U V d e l m s A g i p S S i 1 Y 1 Q d O O y Q l / Q h D C Z T t q h M 0 m Y m Q g l 5 B / c + C t u X C j i 1 o 0 7 / 8 Z J m 4 W 2 H h g 4 c 8 6 9 3 H u P H z M q l W V 9 G 0 v L K 6 t r 6 6 W N 8 u b W 9 s 6 u u b 6 8 0 r g p 4 i i B Q 3 A E T o A N L k E D 3 I E m a A M M H s E z e A V v x p P x Y r w b H 7 P S J a P o O Q B / Y H z + A D E i m + s = < / l a t e x i t > T 1,0 = T ,0 = 50 µs < l a t e x i t s h a 1 _ b a s e 6 4 = " d g O Z u x B Y 2 O p A c u f c Z 7 u b C R z x l z w h x l 9 x 4 0 I R t y 7 c + T d O 2 i y 0 9 c D A m X P u 5 d 5 7 / J h R q S z r 2 y i t r K 6 t b 5 Q 3 K 1 v b O 7 t 7 5 v 5 B V 0 a J w K S D I x a J v o 8 k Y T Q k H U U V I / 1 Y E M R 9 R n r + 5 D b 3 e w 9 E S B q F b T W N i c v R K K Q B x U h p y T P P 2 l 5 q 1 6 w M X k P N n H h M 5 x / b s q B T g w 5 P H I 7 U W P B U Z p 5 Z t e r W D H C Z 2 A W p g g I t z / x y h h F O O A k V Z k j K g W 3 F y k 2 R U B Q z k l W c R J I Y 4 Q k a k Y G m I e J E u u n s q A y e a G U I g 0 j o F y o 4 U 3 9 3 p I h L O e W + r s w 3 l I t e L v 7 n D R I V X L k p D e N E k R D P B w U J g y q C e U J w S A X r H i g 7 B X j x 5 m X T P 6 3 a j 3 r i / q D Z v i j j K 4 A g c g 1 N g g 0 v Q B H e g B T o A g 0 f w D F 7 B m / F k v B j v x s e 8 t G Q U P Y f g D 4 z P H 6 N G n C E = < / l a t e x i t > T 1,0 = T ,0 = 100 µs < l a t e x i t s h a 1 _ b a s e 6 4 = " I G p I o F O k P d 8 d u c 1 a g q F 0 / n k r g 8I = " > A A A C F H i c b V D L S s N A F J 3 4 r P U V d e l m s A i C p S S i 1 Y 1 Q d O O y Q l / Q h D C Z T t q h M 0 m Y m Q g l 5 C P c + C t u X C j i 1 o U 7 / 8 Z J m 4 W 2 H h g 4 c 8 6 9 3 H u P H z M q l W V 9 G 0 v L K 6 t r 6 6 W N 8 u b W 9 s 6 u u b f f k V E i M G n j i E W i 5 y N J G A 1 J W 1 H F S C 8 W B H G f k a 4 / v s 3 9 7 g M R k k Z h S 0 1 i 4 n I 0 D G l A M V J a 8 s z T l p f a V S u D 1 1 A z J x 7 R 2 e f C s q B T h Q 5 P H I 7 U S P B U Z p 5 Z s W r W F H C R 2 A W p g A J N z / x y B h F O O A k V Z k j K v m 3 F y k 2 R U B Q z k p W dR J I Y 4 T E a k r 6 m I e J E u u n 0 q A w e a 2 U A g 0 j o F y o 4 V X 9 3 p I h L O e G + r s w 3 l P N e L v 7 n 9 R M V X L k p D e N E k R D P B g U J g y q

Figure 4 .
Figure 4. Determining optimal transmon parameters.a Schematic representation of the square grid architecture, with high-(dark blue) and low-frequency (light blue) transmons, connected via tunable couplers (green).bValues of E J and E C of Qubit 2 (q 2 ) from different experimental implementations of the tunable coupler transmon architecture.Both parameters are extracted from the reported ω q 2 and α q 2 , and therefore the points are only an approximation50 .The annotations refer to the following references [1] Collodo, M. C., et al. (2020) 46 , [2] Xu, Y., et al. (2020) 44 [3] Sung, Y., et al. (2021) 36 , [4] Wu, Y., et al. (2021) 11 , [5] Xu, H., et al. (2021)37 and[6] Google Quantum AI (2022)8 .The black lines represent constant ratios of E J /E C and the red lines correspond to contours of constant qubit frequency ω q 2 .c Schematic representation of the circuit used to infer the weights in Eqs. 25 and 26, with the state preparation pictured on the left (in purple), single-qubit π and π/2-rotations in green and blue respectively, and the two-qubit CZ gate in orange.N is the number of times the pictured circuit (without the state preparation) is repeated before measurement and therefore an integer determined by the reference coherence time.In our case, we consider N = Figure 4. Determining optimal transmon parameters.a Schematic representation of the square grid architecture, with high-(dark blue) and low-frequency (light blue) transmons, connected via tunable couplers (green).bValues of E J and E C of Qubit 2 (q 2 ) from different experimental implementations of the tunable coupler transmon architecture.Both parameters are extracted from the reported ω q 2 and α q 2 , and therefore the points are only an approximation50 .The annotations refer to the following references [1] Collodo, M. C., et al. (2020) 46 , [2] Xu, Y., et al. (2020) 44 [3] Sung, Y., et al. (2021) 36 , [4] Wu, Y., et al. (2021) 11 , [5] Xu, H., et al. (2021)37 and[6] Google Quantum AI (2022)8 .The black lines represent constant ratios of E J /E C and the red lines correspond to contours of constant qubit frequency ω q 2 .c Schematic representation of the circuit used to infer the weights in Eqs. 25 and 26, with the state preparation pictured on the left (in purple), single-qubit π and π/2-rotations in green and blue respectively, and the two-qubit CZ gate in orange.N is the number of times the pictured circuit (without the state preparation) is repeated before measurement and therefore an integer determined by the reference coherence time.In our case, we consider N =

Figure 5 .
Figure 5.The relative contributions to the value of 1 − P (in Eq. 25) from the error sources listed in Table2(top) and the value of the performance metric P as a function of coherence times with assuming T φ ,0 = T 1,0 at the reference point (bottom).Here we considered the same reference point as used in Fig.4d. a A well-designed system, where the parameters E J and E C are optimally adapted so that P is maximized for each value of the coherence time on the x-axis.b P at various coherence times while keeping E J and E C fixed, corresponding to the values marked as point[5] in Fig.4b.
Illustration of two parity switching mechanisms.The vertical axis represents the energy, with the left and right regions corresponding to the two sides of the Josephson junction with the middle region corresponding to the insulator.The light grey area corresponds to the density of states of a BCS superconductor on both sides of the junction.Orange: A high energy photon breaks a Cooper pair (dark blue) thus generating two quasiparticles (light blue), with one quasiparticle tunneling across the junction.Red: A preexisting quasiparticle tunnels across the junction.c The energy level diagram of the states involved in the operation of a diabatic CPHASE gate.The Rabi oscillation between the levels |11⟩ ↔ |02⟩ is marked with green arrows, however the larger charge dispersion of the second excited state means that both parity states of the |02⟩ level cannot be considered degenerate anymore.d e x i t > |01i Figure1.A schematic representation of the effect of a parity switch on a single transmon and on a two-qubit gate.a The energy diagram of a single offset-charge-sensitive transmon with E J /E C = 10 with two distinct parity manifolds, marked with + and -.While the ground state |0⟩ also comprises of two distinct parity levels, the difference between them is not visible.b

b
Figure 2.Effects of a parity switch on the fidelity of a CZ gate.aAnalyticalpredictions (red circles) of the conditional-phase difference of the CZ gate between the two parities.The theory predictions from Eqs. 11 and 9 are compared to the phase extracted from a full numerical simulation of the gate (blue crosses), with the Hamiltonian parameters from Table1.The inset displays the relative absolute error between the analytical predictions and the numerical results.The x-axis represents the E J /E C ratio (bottom) and second excited state charge dispersion (top) of the qubit with the higher frequency (Qubit 2).The inset displays the relative error of the analytical prediction.b Gate infidelity after a parity switch on Qubit 2, where the theory predictions (red circles) only take into account the effect on the conditional phase, according to Eq. 18.The green lines are estimates of the upper bound for the infidelity contribution of quasiparticle induced decoherence on the two-qubit gate system, with the parity lifetimes T P measured by different references: [1] Risté, D., et al. (2013) 54 , [2] Diamond, S., et al. (2022) 53 and [3] Kurter, C., et al. (2022)

Table 2 .
Scaling analysis for the most relevant error sources in an architecture of flux-tunable transmons, with two noise models.