arxiv: v1 [quant-ph] 2 Jan 2018

Transcription

1 Nonlnear Party Readout wth a Mcrowave Photodetector M. Schöndorf 1 and F. K. Wlhelm 1 1 Theoretcal Physcs, Saarland Unversty, 6613 Saarbrücken, Germany arxv: v1 [quant-ph Jan 018 Robust hgh-fdelty party measurment s an mportant operaton n many applcatons of quantum computng. In ths work we show how n a crcut-qed archtecture, one can measure party n a sngle shot at very hgh contrast by takng advantage of the nonlnear behavor of a strongly drven mcrowave cavty coupled to one or multple qubts. We work n a nonlnear dspersve regme treated n an exact dspersve transformaton. We show that approprate tunng of expermental parameters leads to very hgh contrast n the cavty and therefore to a hgh effcency party readout wth a mcrowave photon counter or another ampltude detector. These tunng condtons are based on nonlnearty and are hence more robust than prevously descrbed lnear tunng schemes. In the frst part of the paper we show n detal how to acheve ths for two qubt party measurements and extend ths to N qubts n the second part of the paper. We also study the QNDness of the protocol. I. INTRODUCTION The platform of crcut quantum electrodynamcs cqed s a promsng canddate for realzng quantum computng n crcuts n a scalable archtecture [1 7. In ths feld, superconductng crcuts are used to realze qubts. The two lowest levels of the nonlnear energy spectrum play the role of the two qubt states. Wavegudes and mcrowave cavtes allow for control and couplng of superconductng qubts [5, 8. Another crucal pont s readout. Currently, readout n superconductng crcuts s mostly realzed usng homodyne feld ampltude detecton [4, Ths scheme requres addtonal devces such as parametrc amplfers to measure the feld ampltudes[10, 1. Whle amplfers are readly avalable, they requre space-consumng mcrowave perpherals such as crculators [13. In [14 we presented a scheme to readout the state of a qubt by couplng t dspersvely to a drven mcrowave cavty and measure f the cavty s brght or dark usng a mcrowave photon counter. It s also possble to measure mult-qubt party states wth ths setup [15. A challenge les n the lmted senstvty of these detectors. Effects lke back reflecton of ncomng photons and wrong rate calbraton lead to photon loss n the counter, such that one needs a relatvely hgh number of photons to actually get a count [16. Ths can be n conflct wth the applcablty of the dspersve approxmaton [17. A way out to ncrease contrast at lmted senstvty s to boost the sgnal and use the nonlnearty of the drven cavty, smlar to how t has been done n the sngle qubt case. In [18 Bossenault et al. studed a M-level system dspersvely coupled to a mcrowave cavty. They showed numercally that gong to hgher drve strengths where n > n crt leads to a nonlnear behavor of the system dynamcs resultng n a huge enhancement of the cavty occupaton that can be used to dstngush the two logcal qubt states. At the same tme Bshop et al. [19 studed the same system wth just two energy levels ncluded. Ths amounts to a bnary pre-measurement of the qubt state. They used the exact dspersve transformaton [0 n a sem-classcal regme to descrbe ths phenomenon mathematcally. The nonlnear effects were also demonstrated n experment for readout of a two level system [19, 1. Here we study ths transton to a nonlnear response of the cavty usng the exact dspersve transformaton and extend t to multple qubts coupled to the transmsson lne whle takng nto account M energy levels of the system representng the qubt. We show how the exact dspersve transformaton s performed for the general case of M energy levels and N qubts and derve an analytcal expresson for the steady state photon occupaton of the cavty dependng on the N-qubt state. Our results match, n the one-qubt case, the numercal results of [18. Analogous to ther results the equatons lead to a strong enhancement n the cavty occupaton dependng on the qubt state. Ths state dependence can only be seen when we nclude hgher energy levels, snce they lead to asymmetrc frequency shfts of the effectve cavty frequency. We furthermore use stablty analyss to derve an expresson for the crtcal drve strength at whch one can observe the strong enhancement n the cavty occupaton. An mportant observaton s that besdes the qubt state, the poston of ths transton also depends on the detunng of the drve frequency and the bare cavty frequency. We show that wth ths dependency one can tune the system such that t s possble to perform any arbtrary two qubt measurement n the logcal bass, ncludng party readout. Whle one drve frequency s enough to perform party measurements for two qubts, we show that one needs N/ dfferent drve frequences to extend the party readout scheme to N qubts. The advantage of the strongly drven regme s the hgh contrast of about 10 5 photons between the dfferent states, such that even a photon detector wth very low effcency can be used to perform the scheme we present n ths paper. Whle t seems possble to use homodyne detecton at frst, here the problem s that an arbtrary detunng between drve and cavty frequency s not possble, because the drve s at the same tme used for readout, whch also causes a phase senstvty of the readout

2 we do not have when we use a mcrowave photon counter nstead. Another pont whch s crucal n ths strong drven regme s the back acton of the hgh cavty occupaton on the qubt state. Snce we want to perform quantum non demolton measurements QND to use the scheme for nstance for quantum error correcton [, 3, the post-measurement qubt state should be the correspondng party egenstate. We look at the effect of decoherence and relaxaton of the qubt n ths regme and show that all the appearng rates of the decoherence channels n the new frame general dspersve frame are of the order of the ncoherent rates n the lab frame. Ths s mportant, to show that ncoherent processes are not orders of magntude greater n the frame we work n. Addtonally we have to study the effect of photon leakage of the cavty on the decoherence of the qubts. Thspapersorganzedasfollows. InSec. IIwelookat the two qubt case ncludng three energy levels per qubt. We present the system of nterest, perform the transformaton and calculate the photon ampltude as well as the correspondng jump postons. Wth these results we present a two qubt party measurement scheme. In Sec. 3 we expand the approach of Sec. II to the general case of N qubts wth M energy levels. We gve the general form of the exact dspersve transformaton and solve t n the same manner as for the two qubt case. Addtonally we present dfferent possble applcatons for readout, lke ground state testng and mult-qubt party measurement. In Sec. IV we take a look at the ONDness of the protocol. In Sec. V we present our concluson. II. TWO QUBIT CASE A. System and Hamltonan Herewelookattwoqubtscoupledtoastrongly,classcally drven mcrowave cavty. Snce most of the current experments use Transmon qubts, whch have a weak anharmoncty, we wll take nto account 3 energy levels nstead of only the two lowest qubt states. In Sec. III where we generalze the whole calculaton, we expand ths to the case of a general M level system. The cavty s addtonally coupled to a mcrowave photon detector, whch s used to dstngush between a brght and dark cavty wthout detectng the phase [14, 15. The setup s shown n Fg 1, where the photon detecton s performed by the Josephson photomultpler [4, 5, but n prncple there are no restrctons on the type of photon detector. The bare qubt and cavty Hamltonan Ĥ0 s gven by Ĥ 0 = ω c â â+ =0 ω 1 ˆΠ 1 + =0 ω ˆΠ. 1 In ths expresson â and â denote the bosonc anhlaton and creaton operator for a cavty mode of frequency ω c, Mcrowave Drve > 1> 0> ω ω 1 10 Mcrowave Cavty Q 1 Q FIG. 1. System schematc. Qubts are coupled to a drven mcrowave cavty. The exstence of photons n the cavty s read out by a mcrowave photon counter. respectvely and ω j s the correspondng frequency of the energy level j, where the upper ndex stands for j the j-th qubt here j = 1,. The operators ˆΠ = j are the projecton operators on the -th qubt energy level of the j-th qubt. To smplfy the calculaton we set = 1. For later applyng the exact dspersve transformaton, we want to rewrte the Hamltonan usng the ˆσ z operators of the two dmensonal subspaces JPM σ j j j z, = ˆΠ 1 + ˆΠ. The result s a Hamltonan that hghlghts transtons where ω j Ĥ 0 = ω c â â+ =1 ω 1 ˆσ 1 z, + =1 are the transformed frequences ω j 1 = 4ωj 10 +ωj 1 3 ω j = ωj 10 +4ωj 1 3 ω ˆσ z,, 3 4, 5 wth ω k j = ω k ω k j. How to perform ths transformaton n general s shown n App. A. Note that for two dfferent qubts t s ω 1 ω. The nteracton between the cavty and the two qubts n the RWA s gven by a Jaynes-Cummngs term [6 for every allowed transton Ĥ nt = Ĥnt,QB1 +Ĥnt,QB 6 = g 1 Î 1 +, + g Î +,, 7 =1 =1 wth couplng strength g j of the partcular transton, nteracton operator Î j j ±, = â ˆσ ±âσ j 8 and the creaton and annhlaton operator of the respectve qubt transton ˆσ j = 1 j 9 ˆσ j = 1 j. 10

3 3 In the case of a Transmon qubt, the couplng matrx elements g j between the correspondng energy levels only depend on that of the 0 1 transton [7 g j = g j To get a full descrpton of the system, we also have to take nto account the classcal cavty drve whch s represented by the Hamltonan [8 Ĥ d = ǫ âe ω dt +â e ω dt, 1 wth drve strength ǫ and drve frequency ω d. Combnng all terms we end up wth the full Hamltonan Ĥ = Ĥ0 +Ĥnt }{{} +Ĥd. 13 Ĥsys To use the setup shown n Fg. 1 for readout we work n the strong dspersve regme. The dspersve regme allows to reach a QND measurement by avodng Rab oscllatons; the strong-dspersve regme allows to resolve all spectral lnes. Addtonally we assume the bad cavty regme, such that we get a herarchy of system parameter constrants, whch can be satsfed n most experments γ 1,γ Φ κ gj ω j, 1 ω g j ω c, 14 c where κ denotes the cavty decay rate and γ 1, γ Φ the qubt decay and dephasng rate, respectvely. As usual, these ncoherent rates need to be smaller than those nduced by the measurement n order to fathfully detect thequbt, elsetwoulddecaybeforethequbtsdetected. B. Exact dspersve transformaton In the low photon number regme n < n crt wth n crt = ω 10 ω c /4g1, one can use the lnear dspersve approxmaton to dagonalze 13. However, we want to gotoregmeswheren n crt and aperturbatveapproxmaton n n/n crt fals to converge. Therefore we use a dfferent approach, the exact dspersve transformaton, whch was ntroduced n [0 and also has been appled n crcut QED [17, 19, 7. The exact dspersve transformaton for two qubts has the parametrc form ˆD = exp =1 Λ 1 1 ˆN Î1, =1 Λ ˆN Î,, 15 j where the Λ ˆN s are functons of ˆNj = â j â+ ˆΠ. Ths operator denotes the exctaton number of the cavty plus the -th energy level. Snce we are n the strong dspersve regme the ˆN s are approxmately good quan- j tum numbers and therefore Λ can be seen as a scalar when performng the transformaton. Before we apply the transformaton we calculate some mportant commutators. It s easy to show that [Îj,1,Ĥ0 = [Îj,,Ĥ0 = ω j 10 ω c }{{} j 1 ω j 0 ω c } {{ } j Î j +,1 16 Î j +,. 17 To smplfy the notaton we ntroduce the followng nested commutator [9 ad A B [A,B ad n AB ad A ad n 1 A B 18 Wth these commutators, we can apply the transformaton on Ĥsys, usng Baker Campbell Hausdorffs formula Ĥsys D = ˆD Ĥ sys ˆD k +1g + j = Ĥ0 + k +1! + k=0 k=0 1 Λj 1 k +1g + j k +1! Λj ad k Λ Î j,1 ad k Λ Î j, Î+,1 Î+,. 19 To get expresson 19 we used some propertes and relatons of the appearng nested commutators that we prove n Appendx D. A more detaled verson of ths calculaton s shown n Appendx C. Here we dsregarded drect two photon transton terms for nstance terms proportonal to â σ 1 σ, snce the probabltes for such transtons are much less than the one photon processes due to the weak anharmoncty of the Transmon potental selecton rules. It s possble to calculate a closed form of the appearng commutators whch reads ad Λ Î k, ad Λ Î k+1, Î+, Î+, = 4 k Λ k N k Î+, = 4 k Λ k+1 N k+1 ˆσ z,, 0 We put 0 nto 19 and end up wth the followng expresson for the transformed system Hamltonan Ĥ D sys = Ĥ0 + [ =1 N q f j f j 1 j,g j,λ j,n j Î+,1 j,g j,λ j,n j ˆσ j z,, 1

4 4 wth f 1 sn Λ N +g cos Λ N N f g sn Λ N + { } 1 cos Λ Nq. 3 j N 4N To obtan a dagonal system Hamltonan we have to choose Λ j such that the Îj +, contrbuton s zero. Settng equal to zero we fnd the followng choce: N j arctan Λ j = λ j N j. 4 wth λ j = g j / j. Fnally we put ths expresson for Λ j nto 19 and end up wth the dagonal system Hamltonan Ĥ D sys = Ĥ0 =1 j 1 1+4λ j N j ˆσ j z,, 5 Ths expresson s exact up to the non party conservng terms we gnored n 19. At ths pont we only moved Ĥsys nto the dspersve frame, but to descrbe the whole setup we addtonally have to transform the drve Hamltonan Ĥd. Snce we are nterested n the regme n n crt, the dve Hamltonan stays n ts orgnal form by gnorng terms of the order n 1/ and λ [19 Ĥ D d Ĥd. 6 For further calculatons t s more convenent to work wth a tme ndependent Hamltonan. Snce Ĥd stll ncludes a tme dependence, we go nto the frame rotatng wth the drve frequency Û = e ˆnω dt. In ths frame the drve s tme ndependent Ĥd = ǫâ + â and the system Hamltonan just ncorporates an addtonal frequency shft n the bare cavty part wth δ c = ω c ω d. Ĥ 0 = δ c â â+ =1 ω j ˆσ j z,. 7 C. Photon ampltude and nstablty The nterestng value whch s crucal for the usage of the setup n Fg. 1 for readout s the cavty occupaton, whch depends on the correspondng state of the qubt. Snce 5 s dagonal, t s relatvely easy to obtan the steady state soluton of the photon ampltude. As mentoned n IIB, we assumethat the qubt occupaton number s constant durng the dynamcs of the system, whch s satsfed because of the dagonal structure of 5 and the strong detunng between cavty and qubt. Therefore the ˆσ j z,s are constant, whch smplfes the followng calculaton sgnfcantly. As a startng pont we use the Louvllan equaton to obtan an equaton of moton for the annhlaton operator of the cavty mode. Addtonally we nclude an ncoherent channel descrbed by the Lndblad operator ˆL κ = κâ [8, whch represents photon loss n the cavty wth rate κ. The adjont master equaton [30 leads to an equaton of moton for the feld operator â n the Hesenberg pcture κ â = [ĤD sys +Ĥd,â 8 â. Puttng n the expressons for ĤD sys and Ĥd, we get â = δ c =1 g j λ j 1+4λ j ˆσ j N j z, κ â ǫ. 9 Conjugaton of 9 leads to the equaton of moton for â. We are nterested n the cavty occupaton n the post rngup state. Usually the steady state descrbes the state reached at t, but for t wth γ 1 t,γ Φ t 1 the qubt state would be completely destroyed. However, for κt 1 the system s n a pseudo steady state, where the behavor s well descrbed by the steady state solutons. Ths s the reason why we work n the bad cavty lmt, such that for ths tme t pseudo we stll meet the condton γ 1 t,γ Φ t 1. Settng â = â = 0 and solvng both equatons for â and â we end up wth an expresson for the photon occupaton n the steady state n = â â = ǫ [δ c χn q + κ 4 wth nonlnear cavty frequency shft χn q = g j λ j 1+4λ j. 30 σ j N j z,. 31 Note that the frequency shft tself depends on the qubt state, snce t ncludes σ j z, such that the photon ampltude depends on the qubt state as well. Another crucal pont s that the N s nclude the photon number n n the cavty, such that 30 represents a transcendental equaton. We can solve the equaton teratvely and the results for some specfc parameters are shown n Fg.. There are three regmes whch can be dstngushed. For low drve strengths we see a lnear response of the cavty up to a crtcal drve strength ǫ 1 and photon number n 1. Ths corresponds to the regon where the system s descrbed by the lnear dspersve approxmaton. After that the ampltude shows a nonlnear behavor bstable

5 5 regon resultng n a strong enhancement of the photon occupaton. Gong to even hgher drve strengths yelds another crtcal pont ǫ,n, where the response of the cavty returns back to a lnear behavor. The specfc values of ǫ 1 and ǫ depend heavly on the state of the qubt. The effectve cavty frequency on the other hand starts at a specfc value whch corresponds to the usual Stark shft and rapdly goes over to the bare cavty frequency n between the regon ǫ 1 < ǫ < ǫ. In the next chapter we wll see that ths nonlnear behavor results from a bfurcaton of the transcendental equaton 30. In the lmt n 0, the expresson for the frequency shft s lm χn q = n 0 =1 g j λ j σ j z,. 3 If we couple one qubt to the cavty and only take the two qubt states nto account, we observe the lnear χ- shft: χ = ±g 1/ 1 see [4. Thus even though the whole calculaton was performed under the assumpton n n crt westll getthe correctexpressonsforsmallvaluesof n, such that we can assume that ourequatons alsoworks well n ths regme. Another pont worth to menton here s that n the case where only two levels are ncluded the frequency shft s completely symmetrc, such that the response when drvng at the bare cavty frequency would be ndependent of the state of the qubt. Ths shows that the pure exstence of hgher levels nfluence the system dynamcs, they do not have to be occuped at all. D. Stablty anlayss As mentoned before, the observed strong nonlnear behavor of the photon ampltude s caused by a bfurcaton of 30. In between the two lnear regmes see Fg. equaton 30 posses a bstable area wth two attractors. Tunng the drve strength through the frst bfurcaton pont, whch appears at ǫ 1, leads to a decson of the cavty dynamcs between the two attractors. Whch of both attractors actually descrbe the cavty state depends on the hstory of the system. In prncple, small fluctuatons nduced by envronment-asssted processes can drve transtons between the two attractors. However, as we see n Fg. the dfference n ampltude between these s about 10 6 photons and these envronmental fluctuatons are assumed to be rather small. Therefore the system tends to stay n the attractor t chooses when ts drven through the frst bfurcaton pont,.e., the transton tme s exponentally long. In ths secton we want to calculate the two crtcal ponts that restrct the bstable area usng stablty analyss smlar to Drummond et al. [31. It s dffcult to use the full expresson 30 for stablty analyss snce n appears n a square root n the denomnator. As we see n Fg. the transton happens around a cavty occupaton cavty occupaton lnear dspersve regme bstable regme logε/mhz cavty frequency shft [MHz logε/mhz FIG.. Photon ampltude top and effectve cavty frequency bottom dependng on the drve strength ǫ for the four dfferent qubt states. The parameters here are ω 10,ω 1,g 1 1 /π = 4.97,4.071,0.1 GHz, ω 10,ω 1,g 1 /π = 4.094,3.868,0.1 GHz, ω c = GHz and δ c = 0. For every state exsts a specfc drve strength ǫ crt, where the frequency rapdly jumps back to the bare cavty frequency and one observes a strong enhancement n the cavty occupaton. of about 10 photons. Therefore t s a good approxmaton to only keep terms up to g 4/ 3, snce for n 10 we stll meet the condton n g 4/ 3 1. Expandng the square root appearng n 5 up to that order we can derve an equaton of moton for the feld ampltudes n the same manner as n the prevous secton [ α t α = [ [ ǫ κ +hn α ǫ [ κ hn α, 33 wth â = α, â = α and photon number n = α α. The functon h depends on the photon number n and s gven by hn = δ c + =1 j λ j λ j4 N j σ j z, 34

6 6 Wth hn the steady state condton ṅ = 0 can be wrtten as ǫ κ = n 4 +h n. 35 Now we assume small fluctuatons αt around the steady state soluton logǫ,an/MHz logǫ crt,plot /MHz TABLE I. Comparson between the analytcal value of ǫ calculated wth 44 and the actual vlaue of ǫ crt n FIG.. We see an almost perfect agreement. αt = α 0 + αt 36 and get a lnearzed equaton for the fluctuaton [ [ α α t α = A α, 37 wth A = n hn n + +hn κ α 0 n hn α 0 hn n hn n n + +hn κ 38 The stablty of equaton 30 s then controlled by the Hurwtz crtera TrA > 0 39 DetA > If these two crtera are fulflled, the egenvalues of the equaton are stable. Therefore the bstablty can only occur f one of the two equatons 40 changes sgn. Snce TrA = κ and we assume to have a cavty decay κ > 0 only the second Hurwtz crteron ndcates an nstablty. The bstable regon s restrcted by the two crtcal ponts that fulfll the condton DetA = 0, whch leads to the followng expresson for the photon number at the crtcal ponts: n 1/ = ω ω 3 4 κ, 41 6χ where we adopted the notaton of [31 by defnng the parameters ω δ c + χ =1 =1 j λ j j λ j4 λ j4 N j QB, σ j z, 4 σ j z, 43 wth N j ˆΠj QB, =. To get the drve strengths ǫ 1 and ǫ correspondng to the two bfurcaton ponts, we have to put the expresson 41 nto the equaton for ǫ 35 κ ǫ 1/ = n 1/ 4 +h n. 44. It s obvous that the photon numbers resultng from equaton 41 has to be postve, such that a bfurcaton only occurs f ω > 3κ /4 and χ ω < 0. The frst of these nequaltes shows that we do not observe a nonlnear behavor f the leakage rate of the cavty s to hgh. On the other hand the second nequalty leads to the fact that the cavty has to be n the blue detuned regme wth respect to the qubt frequences. The second condton also ndcates the borders dotted lne n Fg. 3. Note that ω as well as χ depend on the state of the two qubt subset, such that ǫ 1 and ǫ depend on t as well, whch explans the dfferent poston of the transton n Fg. Whatweknowuptonowsthatthetranstonbetween the two attractors whch we call low and hgh ampltude attractor n the followng occurs at some value n the bstable area ǫ crt ǫ 1,ǫ. The dynamcs of the ampltude depend on the hstory of the system. Startng at a drve strength smaller than the frst bfurcaton pont ǫ 1 and slowly tune t up to hgher drve strengths ams the system to stay n low ampltude attractor untl t reaches the second bfurcaton pont ǫ, where ths attractor no longer exsts and t rapdly jumps nto the hgh ampltude attractor. On the other hand startng at hgher drve strengths than the second bfurcaton pont ǫ leads to a behavor the other way round. The system stays n the hgh ampltude attractor untl t reaches the bfurcaton pont ǫ 1 and then rapdly jumps nto the low ampltude attractor, snce the hgh ampltude attractor does not exst for ǫ < ǫ 1. Therefore the dynamcs of the system depend on how the tunng of the parameter ǫ s performed. The values for ǫ for the parametersn Fg. are gven n Tab. I. Comparng them to the actual values values of ǫ crt n Fg. we see an almost perfect concdence of ǫ wth ǫ crt for all states, whch s due to the fact that we started wth a small photon number when we solved equaton 30 teratvely. In a real experment where one starts wth small drve strength and slowly tunes up the drve strength the system tends to stay n the low ampltudesolutonforeverystateaslongaspossble, hencethe transton n ths case can be assumed to be very closed to ǫ. Therefore we wll assume ǫ to be the actual transton pont of the ampltude n the followng, snce ths s the more reasonable method n experment.

7 7 0 logε crt /MHz party measurement δ [MHz c FIG. 3. Crtcal pont ǫ crt = ǫ dependng on the cavty drve detunng δ c for the same parameters as n Fg.. The dotted lnes ndcate the value of δ c, where the condton ωχ < 0 s no longer fulflled such that no bfurcaton occurs f we ncrease δ c to hgher negatve values and the system stays n the low ampltude attractor for all drve strengths. For a party measurement we have to choose δ c such that we are n the area where the transton to the hgh ampltude attractor occurs frst for 01 and 10 here δ c 0.0. In ths regon 11 does not show a bfurcaton and stays n the low ampltude attractor, such that for the correspondng crtcal drve strength ǫ crt the photon number for 11 s about 0.5 see nset plot. Addtonally we have to drve wth ǫ 01 < ǫ < ǫ 00, such that we do not reach the bfurcaton pont for 00 whch assures that we also have a small photon ampltude f the system s n 00. E. Applcaton to qubt readout and Qubt party measurement In the prevous sectons we assumed that we drve the cavty on resonance δ c = 0 whch corresponds to homodyne detecton schemes, because the drve pulse s also used for state readout [3. Therefore t s not possble to detune drve and cavty frequency arbtrarly. Snce the bfurcaton ponts are fxed by the qubt parameters n ths case, the orderng of the transton dependng on the qubt states s also fxed, e.g. ǫ crt for the 00 state wll always be larger than the other ones. Because of ths, the usage for readout s lmted. If we are nterested n party readout, there s no possblty to dstngush even from odd party states n the case δ c = 0, snce the ǫ crt s of the odd party states le n between the ǫ crt s of the even party states. Currently readout of superconductng qubts s n most cases performed by homodyne detecton schemes [8 11 where the drve and cavty frequency are n resonance or slghtly detuned heterodyne detecton. Therefore the detunng between the drve and the cavty frequency s somehow fxed. However, t s also possble to use a mcrowave photon counter e.g. the JPM for readout. Usng a mcrowave photon counter for the readout process of the cavty gves the possblty to arbtrarly detune the drve from the bare cavty frequency δ c 0. The dependence of the crtcal drve strength ǫ crt on the detunng δ c s shown n Fg. 3. We see that ǫ crt decreases f we go to hgher negatve values of δ c up to a pont, where the condton ωχ < 0 s no longer satsfed and we no longer have a bfurcaton of 30, hence the system stays n the low ampltude attractor lnear regme over the full range of ǫ see dotted vertcal lnes n Fg. 3. To perform two qubt party measurements wth ths set up, one has to drve the system wth a detunng δ c n between the pont where 11 goes over to a stable behavor and the pont where ths happens for 10 crcled area n Fg. 3. In ths regon ǫ crt for the odd party states s smaller than for the 00 state. The 11 state on the other hand stays n the low ampltude soluton, hence the photon number n the cavty f the system s n the 11 state s about 0.5 see Fg. 3 for the correspondng ǫ. All n all wth ths tunng, the cavty s n a low ampltude state f the qubt s n an even party state and vce versa. When we take a closerlook and comparethe frequency shfts n Fg. and Fg. 3, we see that the value of the detunng that gves the border between stable and unstable behavor n Fg. 3 for the respectve state, matches almost perfectly wth the correspondng bare χ-shft at ǫ = 0. We wll see n Sec III that ths behavor can also be observed for more than two qubts. Therefore we can gve an analytc expresson for the optmal drvng pont, f one wants to perform party measurements. The bare χ-shft s gven by 3. The optmal drvng pont les n between the stablty border of 11 and the one of 10 hence s gven by: ω d,opt = ω c + χ 10 +χ 11, 45 where χ j denotes the bare ch shft f the qubts are n the state j. The regme that can be used for party readout hence s bounded by area between two dotted vertcal lnes 11 and 10 n Fg. 3 ω c +χ 11 < ω d < ω c +χ Choosng a detunng n ths regme, whch s about 7 MHz broad for the parameters n Fg. 3, leads to the rght postons of the bfurcaton ponts to perform two qubt party measurements. There s another crucal pont one has to take care of when performng the measurement. We have to drve wth the rght frequency and the rght ntensty ǫ at the same tme. In a real experment, one would start wth a low drve strength and tune the drve strength up nto the regme ǫ 01 crt < ǫ < ǫ00 crt, hold the drve strength n ths regme for κt 1 and then brng the JPM n resonance wth the cavty to read out f the cavty s brght or dark, whch corresponds to odd or even party respectvely. If we tune the system n ths way we have a photon occupaton of about 10 5 photons f the qubts are n an odd party state and about 1 10 f they are n an even party state. Note that there can be a small dfference n photon number between 11 and 00 but compared to the huge

8 8 contrast between the odd and even states, ths does not sgnfcantly nfluence the measurement. By tunng δ c rght one can perform any possble two qubt measurement n the logcal bass, so ths scheme s not restrcted to party measurement. Performng e.g. projectve state measurement just needs to drve the system such that the correspondng ǫ crt s the lowest one. For some of the states we are n the regon where the dynamcs do not show a bfurcaton, but ths agan s no problem because of the same reason as n the party measurement scheme; the dynamcs stay n the low ampltude attractor, hence the photon number s low for the correspondng ǫ. Note that n the case of two dentcal qubts the two odd party states would show exactly the same behavor, whch means ǫ 01 crt = ǫ10 crt. However, n real experments t s often the case that the qubts have dfferent parameters, snce t s hard to produce two completely dentcal qubts. Therefore we assumed slghtly dfferent qubt parameters n Fg 3, and we see that f the parameters not varyng to much, party measurement can stll be performed even f the qubts are not completely ndstngushable. For two qubt party measurements we onlyneed onedrvefrequency. We wll seen Sec III, that ths scheme can be expanded to N qubt party measurements but wth a need of N/ drve frequences, where denotes the floor functon that maps N to the next smaller nteger. III. N QUBIT CASE A. General formulaton and photon ampltude In ths secton we expand our result of Sec. II to N qubts coupled to the readout cavty and we take nto account M energy levels. The bare qubt and cavty Hamltonan Ĥ0 of ths general case has the form Ĥ 0 = ω c â â+ N =1 ω j ˆΠ j, 47 where we used the same notaton as n the prevous secton, but the upper lmts of the two appearng sums are gven by the number of qubts N and the M energy levels taken nto account. Agan we set = 1 for smplcty. For the two qubt case t was not dffcult to operators defned n Eq., snce we just had to solve an equaton system wth two varables. Here we need a general transformaton rule to get the correspondng Hamltonan ncludng rewrte the Hamltonan usng the σ j z, only σ j z, operatorsn the bare qubt part. How to obtan ths transformaton s shown n Appendx A. After ths transformaton we can wrte the bare Hamltonan as Ĥ 0 = ω c â â+ N =1 ω j ˆσ j z,, 48 Mcrowave Drve M-1> 1> 0>... Q 1 Mcrowave Cavty... FIG. 4. System schematc of the general case. N Qubts are coupled to a drven mcrowave cavty. We take nto account M levels. wth the transformaton rule Q N JPM ω j = A 1,k ω k,k 1, 49 where the matrx elements of A are gven by see Appendx A A 1,k = { k M M 1 k M k M k M The nteracton under the RWA leads to a Jaynes- Cummngs term for every possble qubt transton summed up over all qubts Ĥ nt = N =1 g j Î j +,, 51 where the defnton of Îj, s smlar to 8. Snce we are stll assumng Transmon qubts, the couplng matrx elements of the respectve qubt depends on the couplng rate of the correspondng 0 to 1 transton n the same manner as before see Eq. 11. The drve Hamltonan Ĥd does not change n the N qubt case and s therefore gven by 1. The exact dspersve transformaton n the general case s gven by N ˆD = exp Λ j j ˆN Îj,, 5 =1 wth Λ j ˆN defned n 4. Snce [ˆσ z,,ˆσ z,+1 0, the defnton of the s dfferent for hgher levels j = j j 1 ωj j ωj M j, = 1, = M 1, else. ω j 1 + ωj Applyng ths transformaton on the system Hamltonan s a lttle more dffcult than n the prevous secton

9 9 cavty occupaton exctatons 1 exctaton exctatons 3 exctatons 4 exctatons cavty frequency shft [MHz logε/mhz a logε/mhz b FIG. 5. Photon ampltude a and effectve cavty frequency b dependng on the drve strength ǫ. The plot s for four dentcal qubts ncludng 10 energy levels. The parameters here are ω 10,ω 1,g 1 j /π = 4.97,4.071,0.1 GHz, ω c = GHz and δ c = 0. For every set of states wth the same number of exctatons there exsts a specfc drve strength, where the frequency rapdly jumps back to the bare cavty frequency and one observes a strong enhancement n the cavty occupaton. and s done n Appendx B. However, the resultng Hamltonan has a smlar form, except that the appearng sums go to N and M 1 respectvely Ĥ D sys = Ĥ0 N =1 j 1 1+4λ j N j ˆσ j z,. 54 where we agan moved to the frame rotatng wth the drve frequency ω d, such that δ c nstead of ω c appears n Ĥ 0. Lke before we gnored non party conservng transtons, and the drve Hamltonan remans n ts usual form under the transformaton up to order n 1/ and λ. Wth ths Hamltonan we can agan derve equatons of moton for the feld ampltudes and solve the equaton for the photon occupaton n the cavty for the steady state resultng n n = ǫ [δ c χn q + κ 4 55 whch s the same expresson as n secton two see Eq. 30 but the appearng χ shft has more contrbutng terms χn q = N =1 g j λ j 1+4λ j σ j N j z,. 56 The analogy of the expresson we found here and the ones n the prevous secton ndcate that t s very lkely that we also observe a nonlnear behavor comparable to the two qubt case. To see that ths s ndeed the case, we can perform a stablty analyss n the same manner as before. Dong so we get the same expresson for the photon numbers at the two bfurcaton ponts. n 1/ = ω ω 3 4 κ, 57 6χ wth the parameters ω δ c + χ N N =1 =1 j λ j j λ j4 λ j4 N j QB, σ j z, 58 σ j z,. 59 The condtons that the ampltude shows a nstable behavor are the same as for the two qubt case, but wth the changed parameters 58 and 59. Snce we have more qubts and energy levels here, we get a dfferent poston for the bfurcaton ponts for every qubt state assumng that we have slghtly dfferent parameters for every sngle qubt. However, we are only nterested n the occupaton of the two lowest energy levels of the qubts, snce they realze the two mathematcal qubt states needed for quantum computaton. Agan t s the exstence of hgher levels that nfluence the whole system, they do not have to be occuped. Just ncludng the two qubt levels would lead to symmetrc shfts and the values for ǫ crt would no longer be dfferent for all states. The results for the four qubt case s shown n Fg. 5. We see that the system behaves the same as n the two qubt case. The photon ampltude shows a huge enhancementatthesecondbfurcatonpontǫ crt 44. Analytcal expressons for the two bfurcaton drve strengths can be

10 10 obtaned by puttng the general expresson for the bfurcaton photon numbers 57 nto the expresson for the drve strength 35. In Fg. 5 we ncluded 10 energy levels, whch s more than the hghest number of levels that are relevant n practce [18. Otherwse, we show n App. C that n our case only the next lowest level whch s not occuped matters, so n our case M = 3 yelds the correct results. B. Mult-qubt party Measurements In Sec. IIE we have shown how to perform two qubt party measurements n the nonlnear regme usng our setup. Now we want to show that the same can be done for N qubts, we just need more than one drve frequency. To perform an N qubt party measurement we need N/ dfferent drve frequences. In the N qubt case the bfurcaton drve strength ǫ crt depends on the detunng between the drve and the bare cavty frequency as well. Thereforewe can agan take ths as an advantage to tune ǫ crt of the dfferent states such that they ft for party readout. The values of ǫ crt dependng on δ c are shown n Fg. 6, for the case of four dentcal qubts. Hence only qubts wth dfferent exctaton number can be dstngushed. We see that n the four qubt case n between the borders that restrct the nstablty condton les the border of the even party state Therefore t would not be possble to get ǫ crt for the two odd states smaller than for 0011 at the same tme, whch gves rse to the need of two dfferent drve frequences. One between the nstablty borders of 0111 and 1111, whch means ω c +χ 1111 < ω 1 D < ω c +χ 0111 and the other one n between the nstablty borders of 0001 and 0011 : ω c +χ 0011 < ω D < ω c +χ 0001 as shown n Fg. 6. By comparng the frequency shfts n Fg. 5 and Fg. 6 we agan see a concdence wth the nstablty borders and the bare χ-shfts of the respectve states whch explans borders of the two drve frequences. The optmal drve frequences les n the mddle of the respectve regme and are therefore gven by ω 1 D,opt = ω c χ χ 1111 ω 60 D,opt = ω c χ χ 0011, 61 where χ jkl denotes the bare χ-shft of jkl and can be calculated wth 56. When we drve the system wth these two frequences, there exsts a regon of the drve strength max[ǫ 0111 crt ω1 D,ǫ0111 crt ω D < ǫ mn[ǫ 0001 crt ω 1 D,ǫ0000 crt ω D > ǫ 6 where the dynamcs of the two odd party states are descrbed by the hgh ampltude attractor and the dynamcs 0 logε crt /MHz party measurement δ c [MHz FIG. 6. Value of ǫ crt = ǫ dependng on the detunng between drve and cavty frequency for four dentcal qbuts wth the same parameters as n Fg. 5. The dotted lnes ndcate the border, where the nstable behavor dsappears for larger detunngs and the dynamcs are just descrbed by the low ampltude attractor. To measure party we need to drve the system wth two dfferent frequences. One such that δ c les n between the dotted lne of 1111 and 0111 and the other one such that δ c les n between the dotted lne correspondng to 0011 and of the even states by the low ampltude attractor ether snceǫssmallerthen the correspondngǫ crt, orthe state no longer fulflls the condton of nstablty. The calbraton of the experment can be performed n thesamemannerasnthe onedrvefrequencycase. Frst tune the two drves of the system to the rght frequences and then turn up the drve strength nto the regme 6, hold the drve strengths constant for t 1/κ and after that brng the JPM nto resonance to read out the state of the cavty. We can expand ths measurement scheme to N qubts. In ths case we need N/ dfferent drve frequences. Snce we could show that the nstablty borders and the bare cavty shfts are almost dentcal for the two as well as for the four qubt case, we can follow that ths also holds for the N qubt case. The respectve drve frequences have to be n between all nstablty borders of odd and even states as n the two and four qubt case. Let { Ψ } be the subset of odd party states and { Φ j } the subset of even party states of a N qubt system, where and j denote the number of exctatons respectvely. Wth ths notaton the optmal drve strengths are gven by ω D,opt = ω c χ Ψ +χ Φ +1, 63 where the appearng χ-shfts are agan the bare χ-shfts of the correspondng states. Note that n the case of an even number of qubts, t s = 1,...,N/ and j = 1,...,N/+1 and n the case of an odd number of qubts = 1,...,N +1/ and j = 1,...,N +1/.

11 11 IV. IS THE MEASUREMENT PROTOCOL QND? In ths secton we want to take a closer look at the QNDness of our measurement protocol. One crucal pont here s to calculate the transformed ncoherent rates n the new exact dspersve frame and show that they do not ncrease n a much faster way than the orgnal rates, especally do not scale proportonal to the photon number n. Another mportant ncoherent process s the leakage o photons out of the cavty whch s a key pont for the protocol, whch leads to dephasng between superpostons of equal party states. Agan t s to check f ths dephasng rate s n the range of the ntrnsc ncoherent rates of the qubts or f t destroys the QND character of the measurement. A. Transformaton of ncoherent channels for n 1 In ths secton we take a look at qubt ntrnsc ncoherent effects, such as dephasng and relaxaton and how the correspondng Lndblad operators behave under the exact dspersve transformaton. The consequence of such effects n the regme of small photon numbers s well studed [17, 1, usng the Polaron transformaton, therefore we want to see how the system behaves n the regme n whch we are nterested,.e. n 1. Snce we want to perform non-demolton measurements our setup, t s mportant that the appearng dephasng and relaxaton rates are not scalng wth n or some other parameter that s huge n our regme of nterest. The optmal case s f such terms do not appear n the calculatons, or at least lead to untary perturbatons of the qubts, such that we can dmnsh them usng optmal control methods on the correspondng qubt. We want to focus here on the case of the frst secton, where we studed two three level systems coupled to a transmsson lne, agan n the strong dspersve regme. The general case could also be calculated, but t s not necessary here snce all mportant effects that appear n the N qubt case wth M energy levels taken nto account wll also appear n ths easer system. For example leakageto the fourth level wll be less probablethan leakage to the thrd one, such that we get an upper bound for all hgher leakage processes. Furthermore possble nteractons between the qubts nduced by the hgh photon number n the cavty wll also appear n ths smaller system, f they are present, whch s the reason why we nclude the second qubt and not just concentrate on one. As mentoned, we are n the regme n 1: N q = N QB +n n, 64 such that Λ ˆN q acts lke a scalar on qubt operators. Addtonally we assume the semclasscal lmt, such that â α â α. Frst we want to study relaxaton of the frst qubt wth rate γ 1. The correspondng Lndblad operator s ˆL D 1 = γ 1ˆσ z D. Therefore we have to calculate the transformaton of ˆσ 1 : ˆσ 1 = ˆD ˆσ 1 ˆD 65 = e Λ1nI1, +ΛnI, ˆσ 1 e Λ1nI1, ΛNqI, 66 1 = Λ k k! 1nad k I,1 σ 1 +Λ k nad k I, σ 1. k=0 67 The expressons for the frst commutator can be calculated n a closed form ad k 1 I,1 ˆσ 1 = k 1 n k 1 1 kˆσ z,1 68 { ad k ˆσ 1 k = 0 I,1ˆσ 1 = k 1 n k 1 kˆσ x,1 k wth σ x, = Puttng 68 and 69 nto 67 yelds ˆσ 1 D = ˆσ 1 1 arctan λ 1 n ˆσ 1+4λ x,1 λ 1 n 1 n 1+4λ 1 n σ z, k=0 1 k! Λk nad k I,ˆσ 1 71 In the same manner we get a closed expresson for the second operator. Frst we calculate two orders frst order: [I,,ˆσ 1 = α 0 7 second order: α [I,, 0 = α ˆσ. 73 In second order there appears the commutator between I, and ˆσ, whch can be calculated wth the formula 68 and 69, such that we get for k 1 ad k 1 I, = 1 k n k 1 k Λ k 1 nˆσ z, 74 ad k I, = 1 k k 3 n k Λ k nσ x,. 75 Puttng 74 and 75 and the expresson for the frst and second order nto 71, we get the full expresson for ˆσ D 1 whch yelds the followng expresson for the correspondng Lndblad operator: L D 1 = γ 1ˆσ 1 + γ 1 arctan λ 1 n ˆσ 1+4λ x,1 + γ 1 1 n 8 +γ 1 λ 1 n 1+4λ 1 nˆσ z,1 + γ 1 4 arctan 4 λ n ˆσ 1+4λ x, n nλ arctan λ n ˆσ 1+4λ z, 1 n + γ 1 arctanλ nˆσ 1ˆσ γ 1 8 arctan λ nˆσ 76

12 1 We observe as a frst term the orgnal relaxaton wth rate γ 1 and sx addtonal terms. The frst one s a dephasng n the ˆσ x,1 bass, whch leads to a mxng of the ground and frst excted state. The correspondng rate has an arctan n the numerator and a term proportonal to n n the denomnator, whch leads to a rate much smaller than γ 1 for n 1 and at least smaller than γ 1 for every value of n. The same argument holds for the second term, whch leads to a dephasng n the σ x, bass. Here the rate s 1/4 of the rate of the second one and therefore ths term also has no crucal effect on the qubt. The thrd term leads to a dephasng between the ground and frst excted state. The rate here s smaller than the orgnal relaxaton rate γ 1 for all values of n. The fourth term leads to a dephasng between the frst excted state and the second excted state. Here we observe a rate whch wll also be n the order of the orgnal γ 1, snce we have n n the numeratoras well as the denomnator. The ffth and sxth term leads to a nteractons between the second excted state and the ground and frst excted state, respectvely. Anyways the rates here are agan at the order of the orgnal relaxaton rate γ 1. All n all we have shown that relaxaton n ths system leads to addtonal ncoherent effects, but all are happenng wth a rate smaller or comparablewth γ 1. Snce we need condton 14 to perform any quantum mechancal operatons wth the system, nothng crucal happens here. Now we take a look at dephasng between the ground and frst excted state. The correspondng Lndblad operator can be wrtten as L D ϕ = γ ϕσz,1 D. The calculaton can be performed n the same way as for ˆσ 1 D, where we need the followng commutators: ad k 1 I,1 ˆσ z,1 = 1 k 1 k 1 n k 1 σ x,1 77 ad k I,1ˆσ z, = 1 k k n kˆσ z,1 78 ad k 1 I, ˆσ z,1 = 1 k k 1 n k 1ˆσ x, 79 ad k I,ˆσ z,1 = 1 k k 1 n kˆσ z,, 80 whch yelds the followng expresson for the dephasng operator n the new frame: L D ϕ = γ ϕˆσ z,1 arctan λ 1 n γ ϕ ˆσ 1+4λ z,1 γ ϕ arctan λ n ˆσ 1 n 1+4λ z, n λ 1 n γ ϕ 1+4λ nˆσ λ n x,1 γ ϕ nˆσ 1 1+4λ x,. 81 Agan we have the orgnal dephasng term appearng n the Lndblad operator wth rate γ 1. The frst two addtonal terms lead to dephasng between the ground state and the frst excted state and the frst excted state and the second one, respectvely. The rates are extremely small n the regmen 1 such that we can neglect them. The last two terms lead to dephasng n the ˆσ x,1 and ˆσ x, bass, but wth a rate at least smaller than the orgnal dephasng rate γ ϕ. All n all we have shown that we do not have any relevant ncoherent processes wth rates hgher than the relaxaton and dephasng rate of the qubt, whch gves the possblty to perform QND measurements n ths regme. B. Dephasng due to photon leakage In the party measurement protocol we assumed a cavty decay rate κ, whch s greater than the ntrnsc ncoherent rates of the qubts. Ths assumpton s mportant to reach the steady state and measure before the qubt states decay. Snce states wth the same party can lead to dfferent cavty frequency shfts see Fg., the photons leakng out of the cavty carry qubt nformaton. Ths leakage leads to an effectve dephasngsee e.g. [36 of superpostons of party states. We want to study ths process and calculate the respectve dephasng rate. Note that only even party states cause dfferent shfts n the cavty, snce odd party states are n the hgh ampltude attractor, where the frequency s exactly the bare cavty frequency for all states see see Fg.. Therefore we only study the dephasng between even party states. In the measurement protocol, the drve strength s chosen such that the system stays n the low ampltude attractor for even party states. Here the behavor s stll lnearsee, hence we can approxmate the Hamltonan 5 Ĥ Ĥ0 + =1 g j j â âˆσ j z, +Ĥd, 8 wherewe assumedthat wearealreadyn the framerotatng wth the drve frequency. The ncoherent evoluton of the densty matrx s descrbed by the Lndbladan master equaton ρ = [H,ρ+κD[âρ + γ 1, D[ˆσ j + =1 =1 γ Φ, D[ˆσ j 83 z,, whch ncludes three ncoherent processes, photon loss of the cavty wth rate κ, relaxaton wth rate γ 1, and dephasng wth rate γ Φ,. We assume the ntrnsc ncoherent rates to be equal for all qubts. The qubt-cavty densty matrx of an equal superposton of even party states can be wrtten down as for smplcty we label 00 = 0 and 11 = 1 ˆρ = ˆρ c ˆρ c ˆρ c ˆρ c , 84 where ˆρ c j descrbe the feld part of the densty matrx. Puttng ths densty matrx expresson nto 83 we get equatons of moton for the densty matrx elements ˆρ j.

13 13 These can be solved usng the postve-p representaton leadng to the tme evoluton of the densty matrx ˆρt = 1 c j t j α t α j t 85,j=0 For the detaled calculaton see Appendx D. The nduced qubt dephasng s descrbed by the parameter c 10 t = a 10 t α 1 t α 0 t, 86 a sgnfcant decoherence effect on the superposton of equal party states whch leaves the measurement QND. For the N qubt case the result s smlar, but snce there are more than one drve strengths needed, we also have more dephasng channels. However, they are all n the range of a few khz, hence even addng all of them up does not lead to a sgnfcant dephasng. Note that one could as well use the Polaron Transformaton n the manner of [35 to calculate the repsectve dephasng rate wth and a 10 = a 10 0e [ γ+ ωt e [ 4χ1 χ/ t 0 α1t α 0 t dt 87 α 1 t = α s 1 +e[ κ +χ1 χ+δct α 1 0 α s 1 88 α 0 t = α s 0 +e [ κ χ1+δc α 0 0 α s The steady state values of the feld operators are gven by α s ǫ 1 = κ/+χ 1 χ +δ c 90 α s ǫ 0 =. κ/ χ 1 +δ c 91 The results are smlar to [36, but get an addtonal χ contrbuton from the thrd energy level. Ths leads to a no longer symmetrc dependence of the dephasng rate on the detunng δ c, whch can be explaned by the asymmetrc frequency shft of the cavty n the two level case the shft s ±χ, hence symmetrc. Note that we assumed dentcal qubts n the dervaton such that the appearng lnear shfts read χ = g /. The measurement should be performed n a pseudo steady state, when system dynamcs are almost zero. Therefore we assume the lmt κt 1, where the photon leakage nduced dephasng rate can be wrtten as Γ Φ = 4χ 1 χ /Im{α s 1 αs 0 }. 9 Puttng all together we fnally get the followng expresson for the dephasng rate nduced by photon leakage of the cavty Γ Φ = 4κǫ χ κ 4 +δ c +δ cχ +4χ 1 4χ 1χ +κ χ 93 As we see n Fg. 3, for the two qubt party measurement protocol the detunng s δ c 0 MHz and the correspondng drve strength s ǫ 10 MHz. These parameters lead to an effectve dephasng of the qbut wth rate Γ Φ 9 khz. Snce the cavty decay s assumed to be n the range of a few MHz, we can assume that a pseudo steady state s reached before the cavty photon loss has V. CONCLUSION In concluson we have derved a mathematcal descrpton of N superconductng qubts coupled dspersvely to a mcrowave cavty by generalzng the exact dspersve transformaton. We have obtaned that our system of nterest shows a nonlnear behavor for a crtcal drve strength, that results n a huge enhancement 10 5 photons of the photon occupaton n the mcrowave cavty. Ths crtcal drve strength depends on the qubt state and can therefore be used for hgh effcency state readout. Furthermore we have shown ths the state dependent crtcal drve strength can be vared by the detunng between cavty and drve. Due to ths dependence t s possble to perform varous hgh effcency measurements ncludng mult qubt party readout, usng a mcrowave photon counter to measure the cavty occupaton and we have shown how to tune the system to realze these measurements. We gave expressons for the drve frequences to perform mutl qubt party measurements, where one needs N dfferent drve frequences to measure the party of N coupled qubts. Addtonally we studed the effect of relaxaton and dephasng n the hgh occupaton regme and have shown that the appearng ncoherent rates are smaller or equal to the orgnal rates. Also the photon leakage rate of the cavty does not lead to a fast decay of qubt coherence. There are some other ncoherent processes that could be consdered lke broadenng due to photon number varatons of equal party states, but they are all assumed to be not as sgnfcant as the studed process. Ths makes the presented protocol to a canddate for hgh contrast QND party readout. VI. ACKNOWLEDGMENTS We thank Luke Gova, Caleb Howngton, Lukas Thes and Brtton Plourde for frutful dscusons. Ths research s supported by U.S. Army Research Offce Grant No. W911NF