superconducting transmission line resonators can be used as quantum data buses , quantum memories , and single microwave photon detectors  @xcite .
they usually work  in the microwave regime and can also be used as quantum nodes in so - called quantum networks  @xcite .
it is well known that the entanglement is one of the most important resources for quantum information processing  @xcite , and microwave photons play a critical role in quantum state control for solid state quantum devices . therefore , engineering arbitrarily entangled microwave photon states is a very fundamental issue for both solid state quantum information processing and quantum optics  @xcite on superconducting quantum chips .
usually , nonclassical photon states of a single - mode cavity field are generated through the interaction between the cavity field and the two - level atom .
the methods of generating nonclassical photon states can be classified into two ways .
one is to engineer appropriate hamiltonians in different evolution durations by tuning experimental parameters when the target state is being generated  @xcite .
the other one is to obtain the target state via appropriately designed measurements  @xcite .
the former one is deterministic , while the latter one is probabilistic and usually has a low probability to succeed .
if the nonclassical state is generated using natural atomic systems , the latter method is usually more practical since most of parameters are not possible or not easy to be tuned .
however , in artificial atomic systems , the former method is more appropriate because system parameters can be artificially controlled .
for example , superconducting quantum circuits ( sqcs )  @xcite provide us a very convenient way to deterministically engineer nonclassical states of a single - mode microwave field by varying the system parameters  @xcite .
the method of deterministically generating entangled photon states using atomic systems can be tracked to that of generating entangled phonon states of two vibrational modes  @xcite , in a trapped ion interacting with laser fields , by using different sideband transitions .
however , the number of steps in such a method  @xcite exponentially depends on the maximum phonon numbers .
a few proposals were put forward to overcome the exponential dependence of the phonon number by introducing auxiliary atomic energy levels  @xcite , using phonon number dependent interactions  @xcite , or employing multiphonon transitions of high phonon numbers  @xcite .
these methods have successfully reduced the number of steps into quadratic polynomials of the maximum phonon numbers .
the generation of entangled microwave photon states of two modes using superconducting qubit has been studied  @xcite , where a classically driven superconducting qubit with time - dependent frequency is coupled to two microwave fields in two separated cavities .
the interaction hamiltonian between the superconducting qubit and the cavity fields of two modes is described by the jaynes - cummings model .
therefore , there is only single photon transition in each step .
however , the photon - number - dependent stark effects  @xcite induced by the qubit - field coupling make it possible to independently implement operations for photon states .
thus , the number of steps also quadratically depends on the maximum photon number .
it has been shown that the superconducting qubit and the cavity field can have both transverse and longitudinal couplings when the inversion symmetry of the qubit potential energy is broken  @xcite .
the longitudinal coupling can induce multiphoton transitions  @xcite in different sidebands as in trapped ions  @xcite and thus arbitrary photon states of a single - mode cavity field can be more conveniently engineered  @xcite . motivated by studies  @xcite , we study a method to generate entangled microwave photon states in two separated cavities coupled by a superconducting qubit using multiphoton transitions .
we first show that the longitudinal couplings can induce two - mode multiphoton processes similar to those in trapped ions  @xcite , and then study an efficient way to generate superposed two - mode photon states .
the paper is organized as below . in sec .
[ sec : inter - mode ] , an effective hamiltonian , similar to that of trapped ions with two vibrational modes  @xcite , is derived , and then different sideband transitions are discussed . in sec .
[ sec : algorithm ] , a new algorithm is introduced to generate arbitrary superpositions of two - mode photon states . in sec .
[ sec : parsel ] , we discuss how to choose parameters to obtain a high fidelity of the target state . in sec .
[ sec : simresults ] , we numerically study the effects of both imperfect control pulses and the environment on the generated target state .
[ sec : discussions ] , the advantages and experimental feasibility of our method are discussed . finally , we summarize our results in sec .  [
sec : conclusions ] .
and the second one has the frequency @xmath0 .
the coupling strength is @xmath1 ( @xmath2 ) between the qubit and the first ( second ) cavity field .
the qubit is driven by a classical field ( in the middle with the black color ) with the frequency @xmath3 and rabi frequency @xmath4.,scaledwidth=46.0% ]
as schematically shown in fig .
[ fig : model ] , we study a system where a superconducting qubit ( sq ) , modeled as a two level system , is coupled to two single - mode microwave fields in two separated cavities and driven by a classical field .
the system hamiltonian can be given by @xmath5 here , @xmath6 and @xmath7 are the free hamiltonians of the sq and the cavity fields , respectively .
moreover , @xmath8 is the interaction hamiltonian between the sq and cavity fields , and @xmath9 is the interaction hamiltonian between the sq and the classical field . in the qubit basis ,
the qubit hamiltonian is given by @xmath10 with @xmath11 and @xmath12 @xmath13 .
the parameter @xmath14 is the qubit frequency .
the kets @xmath15 and @xmath16 denote the ground and excited states of the qubit , respectively .
the free hamiltonian of two cavity fields is given by @xmath17 where @xmath18 ( @xmath19 is the annihilation ( creation ) operator of the @xmath20th cavity field with its frequency @xmath21 and @xmath22 .
the interaction hamiltonian between the qubit and two cavity fields is @xmath23 where @xmath24 is the coupling strength between the @xmath20th cavity field and the qubit , and @xmath25 is a parameter which depends on the inversion symmetry of the qubit potential energy .
similarly , the interaction hamiltonian between the qubit and classical field is given by @xmath26 where @xmath4 is the coupling strength ( or rabi frequency ) between the qubit and the driving field .
the parameters @xmath3 and @xmath27 are the driving frequency and driving phase , respectively .    in eqs .
( [ eq : hg_tilde ] ) and ( [ eq : hd_tilde ] ) , when the qubit potential energy possesses inversion symmetry , i.e. , @xmath28 , there are only transverse couplings between the qubit and cavity fields  @xcite . if the rotating wave approximation is further made and there is no driving ( @xmath29 ) , eq .
( [ eq:1 ] ) is reduced to extensively studied jaynes - cummings model  @xcite .
when the qubit potential energy possesses a broken inversion symmetry  @xcite , i.e. , @xmath30 , there are both transverse and longitudinal couplings between the qubit and microwave fields .
the broken inversion symmetry of the qubit potential energy can be achieved when the bias charge for the charge qubit or the bias flux for the flux qubit is tuned off the optimal point  @xcite .
but for the phase qubit , the inversion symmetry of the potential energy is always broken @xcite . here
, we will study a general method and not specify a particular qubit .
we now change the qubit basis into the current basis of the flux qubit or the charge basis of the charge qubit .
this is equivalent to diagonalizing the operator @xmath31 . in the new basis , the hamiltonian in eq .
( [ eq:1 ] ) becomes @xmath32 here , the hamiltonians @xmath33 , @xmath34 , and @xmath35 are given by @xmath36 with @xmath37 , and @xmath38 @xmath39 .
hereafter , the parameters @xmath40 and @xmath41 are called transverse and longitudinal frequencies of the qubit , respectively . the kets @xmath42 * * and * * @xmath43 * * are persistent current states of the flux qubit or charge states of the charge qubit . here , * * @xmath44 * * is the rotation operator along the @xmath45-axis , with * * @xmath46 * * @xmath47 . the parameter @xmath48 results in longitudinal couplings between the qubit and microwave fields in eq .
( [ eq:1 ] ) .
below , we will show that @xmath48 can induce two - mode multiphoton processes in the qubit , and then use these multiphton processes to generate arbitrary superpositions of two - mode photon states .      to see how the multiphoton processes can be induced by the longitudinal coupling when @xmath48 , we now apply a unitary transformation @xmath49   , \label{eq : d}\ ] ] to the hamiltonian in eq .
( [ eq : h ] ) .
then , we obtain an effective hamiltonian @xmath50   .\nonumber\end{aligned}\ ] ] it is clear that @xmath51 is the displacement operator  @xcite of two - mode cavity fields .
the displacement quantity is @xmath52 for the @xmath20th cavity field .
hereafter , we will call the picture after the operator @xmath51 as the displacement picture .
the ratios @xmath53 are called the lamb - dicke parameters in analogy to trapped ions  @xcite .    to understand the classical - field - assisted multiphoton transitions of two cavity fields in the qubit
, we apply to eq .
( [ eq : heff ] ) a time - dependent unitary transformation @xmath54   , \label{eq : ud}\ ] ] with @xmath55 .
then , another effective hamiltonian @xmath56   , \label{eq : heffd}\ ] ] can be derived , with the time - dependent term @xmath57   .\ ] ] here , @xmath58 is the bessel function of the first kind .
equation  ( [ eq : heffd ] * ) * shows that multiphoton transitions with different modes can be controlled by the classical field as in trapped ions  @xcite .    in the interaction picture with the free hamiltonian @xmath59 , equation  ( [ eq : heffd ] ) becomes @xmath60 where @xmath61 is the coupling strength between the qubit and cavity field with each different transition process , and its algebraic form is@xmath62 t+in\phi\right\ } \nonumber\\ &   \times\exp\left (   -{\sum_{l}\frac{\eta_{l}^{2}}{2}}\right )   j_{n}\left ( x\right )   \frac{\left (   -1\right )   ^{n_{1}+n_{2}}\eta_{1}^{m_{1}+n_{1}}\eta_{2}^{m_{2}+n_{2}}}{m_{1}!n_{1}!m_{2}!n_{2}!}. \label{eq : jnab}\ ] ] equation  ( [ eq : hint ] )
describes the classical - field - assisted two - mode multiphoton processes as in trapped ions  @xcite .
the magnitude of @xmath63 depends on @xmath64 , @xmath65 , and @xmath66 .
we find @xmath67 where the properties of @xmath68 have been studied in ref .
@xcite . the specific expression * * of @xmath69
is given by @xmath70   t+in\phi\right\ }   .\nonumber\end{aligned}\ ] ] similarly to eq .
( [ eq : jnab_abs ] ) , the magnitude of @xmath69 can be rewritten as @xmath71 it is clear that @xmath72 is independent of the reduced driving strength @xmath65 . from eqs .
( [ eq : jnab_abs])-([eq : jn_abs ] ) , we know that both @xmath73 and @xmath74 can be changed by adjusting @xmath65 and @xmath66 in a similar way . by introducing new variables @xmath75 , we expand eq .
( [ eq : hint ] ) in the fock state basis , and then have @xmath76 with @xmath77 and @xmath78 . here , @xmath79 denotes the ladder operator of the @xmath20th cavity field .
the time - dependent transition element @xmath80 is given by @xmath81 with @xmath82 , @xmath83 replaced by @xmath84 , @xmath85 respectively . the complex transition amplitude @xmath86 and detuning @xmath87 are respectively @xmath88 the parameter @xmath89
is given by @xmath90 with @xmath91{cc}1 , & k<0\\ 0 , & k\geq0 \end{array } \right .   , \\
l_{n}^{\left (   k\right )   } \left (   z\right )    &   = \sum\limits_{l=0}^{n}\frac{\left (   n+k\right )   ! } { \left (   n - l\right )   ! } \frac{\left (   -1\right ) ^{l}}{\left (   l+k\right )   ! } \frac{z^{l}}{l!}.\end{aligned}\ ] ] here , @xmath92 is the generalized laguerre polynomials .
it is clear that the classical - field - assisted multiphoton transitions can be derived from eq .
( [ eq : hint2 ] ) using different frequency - matching conditions .
we now give detailed discussions on how to engineer two - mode multiphoton processes by tuning the driving field .
let us assume that the driving field is tuned to satisfy the resonant condition @xmath93 then eq .
( [ eq : hint2 ] ) can be reduced to an effective hamiltonian @xmath94 when unwanted terms are neglected , that is , @xmath95 the time evolution operator governed by the hamiltonian in eq .
( [ eq : h_sideband1 ] ) is given by @xmath96 recall that @xmath97 and @xmath98 as defined previously . here
the new parameters used in eq .
( [ eq : u_sideband1 ] ) are respectively @xmath99 as shown in eqs .
( [ eq : h_sideband1 ] ) and ( [ eq : u_sideband1 ] ) , @xmath100 photons in the @xmath20th resonator can be either created if @xmath101 or annihilated if @xmath102 while the qubit is flipped up .
similarly , @xmath100 photons in the @xmath20th resonator can be either created if @xmath102 or annihilated if @xmath101 while the qubit is flipped down .
thus different sideband excitations can be constructed , depending on the values of @xmath103 and @xmath104 .    because the hamiltonian derived in eq .
( [ eq : hint ] ) is similar to that of the trapped ions  @xcite , the algorithm using two - mode multiphonon processes in trapped ions can be directly applied into our model , and different superpositions of two - mode photons can be generated .
as a special case , two - mode fock states of high photon numbers can in principle be more efficiently generated with just two steps as single - mode fock states of high photon numbers  @xcite .
however , we here design a new algorithm via different sideband transitions of low photon numbers by tuning the driving field with properly selecting the parameters @xmath105 , @xmath64 , @xmath21 , and @xmath66 .
the detailed discussions of parameter selection will be given in sec .
[ sec : parsel ] .
let us first study a universal algorithm for generating arbitrary two - mode microwave photon states using sideband transitions with the following four hamiltonians @xmath106 , @xmath107 , @xmath108 , and @xmath109 . here , for the compact of notations
, we have used @xmath110 to represent @xmath111 with @xmath112 .
for instance , @xmath106 is actually @xmath94 with @xmath113 , @xmath114 , and @xmath115 . for different @xmath116 , @xmath103 , and @xmath104 , the interaction hamiltonian @xmath94 and its time evolution operator @xmath117 have already been given in sec .
[ sec : eng_inter - mode ] .
below , we will first study how to generate the target state by choosing pulse durations , frequencies , and phases of the driving fields at each generation step with different sideband excitations , and then we will apply our algorithm to the generation of noon states and discuss particular properties of the algorithm .
we note that the state generation in our algorithm is studied in the displacement picture with the unitary transformation as shown in eq .
( [ eq : d ] ) .
the arbitrary quantum states , we expect to be generated , is written as@xmath118 where @xmath119 means that the first and second cavities contain @xmath120 and @xmath83 photons , respectively , and @xmath121 means that the qubit is in the ground state .
besides , @xmath122 and @xmath123 mean the maximum photon number and the probability amplitude on the state @xmath124 , respectively .
we assume that the system is initially in the state @xmath125    we suppose the target state @xmath126 can be generated by alternately switching on and off the two - mode transitions @xmath106 , @xmath107 , @xmath108 , and @xmath109 . with the designed time evolution operators
, the state generation procedure can be represented by,@xmath127   \left\vert \psi _ { 0}\right\rangle \nonumber\\ =   &   \left [   \bar{u}_{f}^{\dag}\left (   t_{f}\right )   u_{\bar{1}}^{p_{f}}\left (   t_{f}\right )   \bar{u}_{f}\left (   0\right )   \right ]   \cdots\nonumber\\ &   \left [   \bar{u}_{1}^{\dag}\left (   t_{1}\right )   u_{\bar{1}}^{p_{1}}\left ( t_{1}\right )   \bar{u}_{1}\left (   0\right )   \right ]   \left\vert 0,0\right\rangle \left\vert g\right\rangle , \label{eq : evolution}\ ] ] where @xmath128 denotes the transition type for the @xmath129th step , and @xmath130 is the time duration for the @xmath129th step .
the time evolution operator @xmath131 is given by @xmath132 \nonumber\\ &   \times\exp\left [   ix\frac{\sigma_{z}}{2}\sin\left (   \tilde{\omega}_{\nu } t_{\nu}+\phi_{\nu}\right )   \right ]   .
\label{eq : ubar}\ ] ] as discussed above , the transitions of different types can be achieved by changing the frequency @xmath3 of the driving field , which is denoted by @xmath133 for the @xmath129th step .
the phase of the driving field for the @xmath129th step is denoted by @xmath134 .
we can express eq .
( [ eq : evolution ] ) in another equivalent form of iteration , @xmath135 with @xmath136 and @xmath137 given in eqs .
( [ eq : psi_0 ] ) and ( [ eq : psi_f ] ) , respectively .
the ket @xmath138 is the state after the @xmath129th step .
we note that the subscript @xmath139 of @xmath137 in eq .
( [ eq : evolution ] ) denotes the number of the final step .
equation  ( [ eq : iteration ] ) means that the initial state is restored from the target state by a composition of sideband transitions with proper time durations , frequencies and phases of driving fields .
it is a recursion algorithm .    .
the @xmath82 and @xmath83-axis respectively denote the photon number of the first and second mode .
the two - mode photon state is denoted by @xmath140 .
the qubit state is represented by the @xmath141-axis with @xmath142 or @xmath143 respectively denoting the ground state @xmath144 or excited state @xmath145 .
the state component @xmath146 is represented by a block at the location @xmath147 .
if a state component is occupied , we color the corresponding block with red ; otherwise , the block is left uncolored .
the arrows respectively represent the 10 , 01 , 11 , and 00  transitions with transition types labeled aside them .
the solid arrow indicates a desired population transfer from the starting state to the end state , while the dashed arrow indicates the inevitable irrelevant oscillation when the desired population transfer is implemented .
( a ) schematic diagram for transferring the populations on states @xmath148 , @xmath149 , and @xmath150 to the state @xmath151 .
this is achieved by consecutively using 01 , 10 , 01 , and 10  transitions .
( b ) schematic diagram for transferring populations on states @xmath152 and @xmath151 to the state @xmath153 to the state @xmath154 .
this is achieved by using a 10  transition .
( d ) schematic diagram for transferring the  population on the state @xmath155    without loss of generality , we use the maximum photon number @xmath156 as an example to show our algorithm .
the more general case with arbitrary @xmath122 is given in appendix .
[ append : univ_alg ] .
the detailed steps for generating the target state@xmath157 with @xmath156 using our recursion algorithm are described as the following four procedures .    * procedure ( i)*. as schematically shown in fig .
[ fig : general](a ) , from the final state @xmath137 , we first transfer the populations in the state space spanned by @xmath158 to the state @xmath151 .
this procedure consists of four steps as schematically shown in below @xmath159{0\bar{1}}\left\vert 0,1\right\rangle \left\vert e\right\rangle \xrightarrow[f-1]{\bar{1}0}\left\vert
1,1\right\rangle \left\vert g\right\rangle \xrightarrow[f-2]{0\bar{1}}\nonumber\\ &   \left\vert 1,0\right\rangle \left\vert e\right\rangle \xleftarrow[f-3]{\bar{1}0}\left\vert 2,0\right\rangle \left\vert g\right\rangle .\label{eq : proc_i}\ ] ] in eq .
( [ eq : proc_i ] ) , the transition type and step number is labeled respectively above and below the arrow .
the arrow points to the direction of the population transfer . in step @xmath129 , the population transfer is accomplished by properly tuning @xmath133 , @xmath130 , and @xmath134 .
after this procedure , we obtain the state @xmath160 which only has popupaltions in the space @xmath161 .
we note that two additional oscillations @xmath162 will also occur inevitably when the population transfer from the state @xmath163 to the state @xmath164 is implemented .
but they do not cause population leakage outside the original space and no extra steps should be taken for them .
we call these oscillations as irrelevant oscillations . for this procedure ,
these irrelevant oscillations are schematically shown by dashed arrows in fig .
[ fig : general](a ) .
the irrelevant oscillation can also occur in the following procedures and is shown by dashed arrows .
* procedure ( ii)*. as schematically shown in fig .
[ fig : general](b ) , starting from the state @xmath160 , we need to transfer the populations in the state space spanned by @xmath165 to the state @xmath166 .
this procedure consists of following two steps @xmath167{1\bar{1}}\left\vert 1,0\right\rangle \left\vert e\right\rangle \xrightarrow[f-5]{00}\left\vert 1,0\right\rangle \left\vert g\right\rangle .\ ] ] after this procedure , we obtain the state @xmath168 , which only has popupaltions in the space @xmath169 .    * procedure ( iii)*. this procedure is similar to procedure ( i ) . as schematically shown in fig .
[ fig : general](c ) , starting from the state @xmath168 , here we need to transfer the population on the state @xmath166 to the state @xmath154 .
this procedure consists of only one step as below @xmath170{\bar{1}0}\left\vert 0,0\right\rangle \left\vert e\right\rangle .\ ] ] after this procedure , we obtain the state @xmath171 , which only has popupaltions in the space @xmath172 .
* procedure ( iv)*. this procedure is similar to the procedure ( ii ) . as schematically shown in fig .
[ fig : general](d ) , starting from the state @xmath171 , we need to transfer the population on the state @xmath154 to the state @xmath173 .
this procedure only consists of one step as below @xmath174{00}\left\vert 0,0\right\rangle \left\vert g\right\rangle .\ ] ] we thus obtain the state @xmath175 .
therefore , the target @xmath126 can be generated from the initial state @xmath176 using inverse processes from the procedure ( iv ) to the procedure ( i ) .
we note @xmath177 .
therefore , the generation of the target state with @xmath178 needs @xmath179 steps .
our algorithm takes a quadratic number of steps while an exponential one is required in ref .
let us now analyze the reason .
our algorithm employs four interaction hamiltonians , @xmath107 , @xmath108 , @xmath109 , and @xmath106 given in eq .
( [ eq : h_sideband1 ] ) . however , four interaction hamiltonians @xmath180 , @xmath181 , @xmath182 , and @xmath183are employed in ref .
the former three interaction hamiltonians between our algorithm and those in ref .
@xcite are qualitatively identical since they convert the same number of bosons for either mode when the two - level system is excited . however , the last ones show fundamental difference between our algorithm and that in ref .
@xcite , because ours creates one boson ( photon ) of one mode but annihilate one boson ( photon ) of the other when the two - level system is excited .
but in ref .
@xcite , one boson for both modes can be simultaneously created when the two - level system is excited .
this difference is critical for us to design an algorithm which can keep track of the populations with a constant total boson ( photon ) number .
therefore , there is no population leakage outside the original space .
however , the algorithm in ref .
@xcite has population leakage . obviously , if the last interaction in ref .
@xcite is changed to @xmath184 , a theoretically equivalent algorithm to ours can also be developed . in this sense
, our algorithm can be regarded as the improved version of that in ref .
@xcite .
let us now study how to choose the pulse duration @xmath130 , the frequency @xmath133 and phase @xmath134 of the driving field to generate a target state in the @xmath129th step for different types of transitions .
we suppose that the population transfer is taken as following @xmath185{p_{\nu}}\left\vert n_{1}+k_{1},n_{2}+k_{2}\right\rangle \left\vert e\right\rangle , \ ] ] in the @xmath129th step , where the transition type @xmath186 should be switched on based on the previous discussions . thus the driving frequency is taken as @xmath187 from the resonant condition in eq .
( [ eq : res_con ] ) . by introducing the notations@xmath188 then from  eq .
( [ eq : iteration ] ) , we need to solve the equation,@xmath189 we thus have the explicit solution for the pulse duration @xmath130 as @xmath190 the phase of the driving field is determined by @xmath191 here , the notation @xmath192 is the value of @xmath193 for the @xmath129th step , which is given in eq .
( [ eq : phink ] ) and depends on @xmath134 .
still recall @xmath194 with @xmath195 .    similarly , if the population transfer is taken as @xmath196{p_{\nu}}\left\vert n_{1},n_{2}\right\rangle \left\vert g\right\rangle .\ ] ] in the @xmath129th step .
the explicit solution for @xmath130 is then@xmath197 and the phase of the driving field is determined by @xmath198    according to the target state , the time duration , frequency and phase of the driving field for each step can be calculated using above equations .
for example , if the @xmath199 " transition is used in the @xmath200rd step , then we use eq .
( [ eq : t_nu_e ] ) and eq .
( [ eq : phi_nu_e ] ) to obtain @xmath201 and @xmath202 by setting @xmath203 .
as an example , we now apply our algorithm to the generation of the noon state , i.e. , the target state is @xmath204 the recursion algorithm restoring @xmath126 to the vacuum state @xmath173 is schematically shown in fig .
[ fig : noon ] for the maximum photon number @xmath156 . in fig .
[ fig : noon](a ) , we can find that all the populations in the hilbert space spanned by @xmath205 can be transferred to the state @xmath151 by consecutively using transitions 01 `` , 10 '' , 01 `` , and 10 '' , i.e. , @xmath159{0\bar{1}}\left\vert 0,1\right\rangle \left\vert e\right\rangle \xrightarrow[f-1]{\bar{1}0}\left\vert 1,1\right\rangle \left\vert g\right\rangle \xrightarrow[f-2]{0\bar{1}}\nonumber\\ &   \left\vert 1,0\right\rangle \left\vert e\right\rangle \xleftarrow[f-3]{\bar{1}0}\left\vert 2,0\right\rangle \left\vert g\right\rangle .\end{aligned}\ ] ] after this procedure , as schematically shown in fig .
[ fig : noon](b ) , all the populations on the state latexmath:[$|1,0\rangle    @xmath173 by consecutively using transitions 00 `` , 10 '' , and 00 " , i.e. , @xmath207{00}\left\vert 1,0\right\rangle \left\vert g\right\rangle \xrightarrow[f-5]{\bar{1}0}\nonumber\\ &   \left\vert 0,0\right\rangle \left\vert e\right\rangle \xrightarrow[f-6]{00}\left\vert 0,0\right\rangle \left\vert g\right\rangle .\end{aligned}\ ] ] the total step number is thus @xmath208 for generating the noon state @xmath209 .
more generally , given an arbitrary @xmath122 , the total step number for generating the noon state in eq .
( [ eq : noon ] ) is @xmath210 the step number for generating noon state has been greatly reduced in comparison with that for generating an arbitrary state in eq .
( [ eq : f_general ] ) .
obviously , the noon state can be generated without using the 11  transition .
if we assume the lamb - dicke parameter @xmath211 , which is usually the case even in the ultrastrong regime in superconducting circuit qed systems  @xcite . from eq .
( [ eq : rabi ] ) , we have the rabi frequencies @xmath212 , @xmath213 , @xmath214 , and @xmath215 .
thus , the transition 11  generally takes more time among the four types of transitions employed by us .
therefore , our algorithm may show a better efficiency for generating noon sates than generating arbitrary entangled states .
this is especially true when the maximal photon number @xmath122 is higher and the lamb - dicke parameter @xmath66 is smaller .    .
( a ) schematic diagram for transferring the population in the space latexmath:[$\{|0,2\rangle    @xmath151 .
this is achieved by consecutively using 01 , 10 , 01 ,  and 10  transitions .
( b ) schematic diagram for transferring the population on the state @xmath151 to the state @xmath173 .
this is achieved by consecutively using 00 , 10 ,  and 00  transitions .
, scaledwidth=48.0% ]
in all of the above studies , we make an approximation that all unwanted terms have been neglected
. however , these neglected terms will affect the fidelity of the prepared target state .
let us now discuss how to minimize the effect of these unwanted terms in eq .
( [ eq : hint2 ] ) on the target state by choosing appropriate parameters . in principle , the effects of these unwanted terms can be perfectly removed by pulse calibration techniques . here
, we study a method to minimize the effect of these unwanted terms by choosing the parameters when the pulse calibration can not be used .    in our algorithm ,
we have used four interactions @xmath107 , @xmath108 , @xmath106 , and @xmath109 , all of them are constructed by the terms with the bessel function @xmath217 in eq .
( [ eq : hint2 ] ) . here ,
in the subscript of the bessel function , we also use @xmath218 to denote @xmath219 if @xmath220 .
we hope to suppress all the terms with the bessel functions @xmath221 for @xmath222 .
we focus on the case * * @xmath223 considering possible experimental conditions . in this case , only lower order bessel functions @xmath224 , @xmath225 , and @xmath226 play significant roles .
thus , we need only to find proper parameters such that the effect of the terms with @xmath224 , @xmath227 , and @xmath226 are negligibly small .
our idea is to make those terms nonresonant by properly choosing the parameters @xmath64 and @xmath105 of the qubit , and frequencies @xmath228 and @xmath0 of two microwave modes .
that is , we assume that the frequency of the @xmath20th cavity mode satisfies @xmath229 where @xmath230 is a positive integer and @xmath231 is the greatest common divisor of @xmath228 and @xmath0 . assuming that the @xmath232transition is switched on , i.e. , the transition detuning @xmath233 , then from eq .
( [ eq : detuning ] ) , the frequency @xmath3 of the driving field must satisfy the condition @xmath234 from eqs .
( [ eq : detuning ] ) ,  ( [ eq : oil ] ) ,  and  ( [ eq : omega_wave ] ) , the detuning of the term with @xmath235 , @xmath236 , @xmath237 is then given by @xmath238 thus the terms with the bessel function @xmath227 will have the detuning @xmath239 we expect that the terms with @xmath227 are nonresonant .
thus , the relation that @xmath240 must hold . a simple but sufficient condition
is @xmath241 where @xmath242 is an integer . similarly , for the terms with @xmath243 , @xmath244 , and @xmath245 ,
the sufficient conditions can be given by @xmath246 the conditions in eqs .
( [ eq : con1])-([eq : con_2 ] ) can be summarized as @xmath247 we can also assume that the longitudinal frequency @xmath105 of the qubit is @xmath248 where @xmath249 is the integer part , and @xmath250is the fraction part . to meet eq .
( [ eq : con_omegaz ] ) , there should be @xmath251    the nonresonant terms with @xmath252 , @xmath253 , and @xmath254 still have effect on the desired time evolution .
these effects can be further eliminated by decreasing the stark shifts caused by the terms with @xmath255 in eq .
( [ eq : rabi ] ) . the ideal case is@xmath256 or equivalently,@xmath257 for @xmath258 , @xmath259 , @xmath260 , @xmath261 , and @xmath262 , where the constraint condition for @xmath120 and @xmath263 denotes the working space of our algorithm .
equation  ( [ eq : con_rabi0 ] ) means that the stark shifts should be negligibly smaller than the rabi frequencies for state generation . considering that @xmath122 is the maximum photon number of the target state , and using eq .
( [ eq : dnk1k2p ] ) and eq .
( [ eq : ozform ] ) , we can obtain * * @xmath264 here , @xmath265 means @xmath65 rounded down and @xmath266 means @xmath65 rounded up .
we thus reduce eq .
( [ eq : con_rabi ] ) to @xmath267 a condition much stronger than eq .
( [ eq : con_rabi ] ) . if eqs .
( [ eq : con_rabi2 - 1])-([eq : con_rabi2 - 4 ] ) are fulfilled , the nonresonant terms can in principle be suppressed . we know from eq .
( [ eq : rabi ] ) that eqs .
( [ eq : con_rabi2 - 1])-([eq : con_rabi2 - 4 ] ) can be satisfied if , for example , the parameter @xmath64 of the qubit is tuned sufficiently small , assuming that the reduced driving frequency @xmath65 and lamb - dicke parameters @xmath66 have been appropriately chosen .    beside the terms with bessel functions @xmath268 where @xmath269 , there are also unwanted terms with the bessel function @xmath217 , which , however , also satisfy the resonant condition @xmath270 here , we have used eqs .
( [ eq : omega_wave ] ) and ( [ eq : dnk1k2p ] ) to obtain eq .
( [ eq : unrester ] ) .
the lamb - dicke parameters satisfy the condition @xmath271 for circuit qed systems even in the ultrastrong regime  @xcite . from eq .
( [ eq : oil ] ) , we know that @xmath272 and @xmath273 are coprime numbers .
we can further make @xmath272 ( or @xmath273 ) sufficiently large .
thus the unwanted resonant terms will possess large @xmath274 ( or @xmath275 ) . in this way
, the effects of these terms will be suppressed due to the exponential decrease via the term @xmath276 in eq .
( [ eq : rabi ] ) .
the condition , that the term @xmath277 is negligibly small , can be summarized as that @xmath272 and @xmath273 should satisfy @xmath278    we now summarize the condition that minimizes the effects of unwanted terms .
the parameter @xmath105 of the qubit should satisfy eq .
( [ eq : ozform ] ) and eq .
( [ eq : con_r ] ) .
however , the parameter @xmath64 of the qubit is mainly constrained by current experiments .
for example , typical values of @xmath279 are in the range @xmath280 ghz .
the frequencies of the cavity modes @xmath21 should satisfy eq .
( [ eq : oil ] ) and eq .
( [ eq : con_li ] ) .
the values of the reduced driving frequency @xmath55 and lamb - dick parameter @xmath53 should satisfy eq .
( [ eq : con_rabi ] ) or stronger conditions eqs .
( [ eq : con_rabi2 - 1])-([eq : con_rabi2 - 4 ] ) .
appropriate values of @xmath65 and @xmath66 can be obtained via numerical simulations , which will be discussed below in sec .
[ sec : unter ] .
we now further numerically simulate the effect of the unwanted terms on the generation of target states by using examples of generating the following two target states @xmath281 for some given parameters .
it is obvious that @xmath282 is an entangled state where every state component is evenly occupied .
we thus call @xmath282 the evenly - populated state .
the state @xmath283 is a two - photon noon state  @xcite .
both @xmath282 and @xmath283 possess a maximum photon number @xmath156 .
the fidelities for generating these two states @xmath282 and @xmath283 are defined as @xmath284 here , @xmath285 and @xmath286 are respectively the actually generated states via the total hamiltonian in eq .
( [ eq : h ] ) .
we now determine the detailed experimental parameters . from eq .
( [ eq : ozform ] ) and eq .
( [ eq : con_r ] ) , we set @xmath287 @xmath288 , @xmath289 , and @xmath290 ghz , which corresponds to @xmath291 ghz . from eq .
( [ eq : con_rabi ] ) or eqs .
( [ eq : con_rabi2 - 1])-([eq : con_rabi2 - 4 ] ) , the parameter @xmath279 should be made smaller , e.g. , we set @xmath292 ghz .
the frequency of the @xmath20th cavity , i.e. , @xmath21 , is determined by eq .
( [ eq : oil ] ) and ( [ eq : con_li ] ) .
since the microwave fields are usually of several gigahertz , here we set @xmath293 , and @xmath294 , thus yielding @xmath295 ghz and @xmath296 ghz .
the lamb - dicke parameters for the first and second cavity modes are set to be identical , i.e. , @xmath297    we vary the lamb - dick parameter @xmath298 and the reduced driving frequency @xmath55 to simulate the effect of the unwanted terms on the fidelity of the expected target states in eqs .
( [ eq : psig ] ) and ( [ eq : psinoon ] ) .
the pulses are taken according to the calculation of chapter iii b. that is , in the @xmath129th step , we use a sinusoidal driving with the driving frequency @xmath133 .
since the sinusoidal driving lasts for a duration @xmath130 , the driving field can be considered as square - windowed sinusoidal signal and thus , strictly speaking , is not delta - shaped in the spectrum .
the simulation results for generating target states in eq .
( [ eq : psig ] ) and eq .
( [ eq : psinoon ] ) are listed in table  [ tab : psig ] and table  [ tab : psinoon ] , respectively .
we can easily find that larger reduced driving strengths @xmath65 and lamb - dick parameters @xmath298 can usually make the fidelity higher . for the evenly - populated state @xmath299 , the largest fidelity 0.939 can be obtained at @xmath300 and @xmath301 . however , for the noon state @xmath302 , the largest fidelity 0.92 can be obtained at @xmath303 and @xmath304 .
[ c]c|c|cccccc & + & 0.2 & 0.3714 & 0.4571 & 0.5429 & 0.6286 & 0.7143 + & 0.3 & 0.115 & 0.393 & 0.519 & 0.608 & 0.739 & 0.675 + & 0.7857 & 0.589 & 0.817 & 0.872 & 0.851 & 0.911 & 0.852 + @xmath65 & 1.0286 & 0.615 & 0.85 & 0.876 & 0.906 & 0.886 & 0.857 + & 1.2714 & 0.724 & 0.872 & 0.886 & 0.906 & 0.852 & 0.878 + & 1.7571 & 0.821 & 0.939 & 0.899 & 0.859 & 0.825 & 0.837 + & 2 & 0.867 & 0.915 & 0.838 & 0.876 & 0.887 & 0.859 +    [ tab : psig ]    [ c]c|c|cccccc & + & 0.2 & 0.3714 & 0.4571 & 0.5429 & 0.6286 & 0.7143 + & 0.3 & 0.108 & 0.34 & 0.403 & 0.395 & 0.461 & 0.687 + & 0.7857 & 0.675 & 0.815 & 0.833 & 0.862 & 0.829 & 0.868 + @xmath65 & 1.0286 & 0.78 & 0.857 & 0.846 & 0.867 & 0.883 & 0.813 + & 1.5143 & 0.871 & 0.877 & 0.873 & 0.877 & 0.787 & 0.806 + & 1.7571 & 0.844 & 0.918 & 0.889 & 0.902 & 0.862 & 0.819 + & 2 & 0.876 & 0.909 & 0.832 & 0.92 & 0.853 & 0.806 +    [ tab : psinoon ]
in * * the above , we only discuss the effect of unwanted terms on the generation of target states .
we now study the effect of dissipation on the fidelities of target states by numerical simulation for given parameters .
when the environmental effect is included , the dynamical evolution of the sqc can be described by the master equation @xmath305   + \mathcal{d}\left [   \sqrt{\gamma_{eg}}\tilde{\sigma}_{ge}\right ]   \rho+\mathcal{d}\left [   \sqrt{\gamma_{ee}}\tilde{\sigma}_{ee}\right ]   \rho\nonumber\\ &   + \mathcal{d}\left [   \sqrt{\gamma_{gg}}\tilde{\sigma}_{gg}\right ] \rho+\mathcal{d}\left [   \sqrt{\kappa_{1}}a_{1}\right ]   \rho+\mathcal{d}\left [ \sqrt{\kappa_{2}}a_{2}\right ]   \rho , \label{eq : me_schr}\ ] ] where @xmath306 and @xmath307 are the reduced density operator and the hamiltonian of the whole system , respectively . the total hamiltonian has been given in eq .
( [ eq : h ] ) .
the compact notation @xmath308 \rho=\left (   2c\rho c^{\dag}-c^{\dag}c\rho-\rho c^{\dag}c\right )   /2 $ ] represents the lindblad - type dissipation .
we have noted that @xmath309 is the basis of @xmath310 , but the qubit dissipation is determined by the qubit basis @xmath311 .
the ground ( @xmath312 ) and excited ( @xmath313 ) states of the qubit are given by the eigenstates of eq .
( [ eq : hq ] ) .
if we define @xmath314 with @xmath315 and @xmath316 , and also define @xmath317 with @xmath318 and @xmath319 .
we can easily verify @xmath320 where @xmath321 , and @xmath322 . in eq .
( [ eq : me_schr ] ) , @xmath323 is the pure - relaxation rate from the qubit excited state to the ground state . besides , @xmath324 and @xmath325 are the pure - dephasing rates originating from disturbed qubit eigenstates
. the decay rates of the first and the second cavity fields are denoted by @xmath326 and @xmath327 , respectively .    using parameters in sec .
[ sec : unter ] and taking the reduced driving * * strength @xmath328 and lamb - dicke parameter @xmath301 from table  [ tab : psig ] and table  [ tab : psinoon ] , we find that the highest fidelity @xmath329 is achieved for generating the evenly - populated state @xmath299 in eq .
( [ eq : psig ] ) , and a high fidelity @xmath330 is also reached for generating the noon state @xmath331 in eq .
( [ eq : psinoon ] ) .
we now assume that the decay rates in eq .
( [ eq : me_schr ] ) are taken as @xmath332 , @xmath333 mhz , and @xmath334 mhz .
we assume that the density operators @xmath335 and @xmath336 are the actually generated states for the target states @xmath299 and @xmath331 .
then the fidelities can be redefined as @xmath337 we perform numerical simulations using the above parameters and obtain @xmath338 @xmath339 and @xmath340 .
the total time for generating @xmath299 is @xmath341 and that for generating @xmath331 is @xmath342 .
both @xmath343 and @xmath344 are too small to induce significant decoherence at the decay rates specified by us .
thus , the fidelity losses induced by dissipation are fairly small , which are @xmath345 for the evenly - populated state @xmath346 and @xmath347 for the noon state @xmath331 .
we now discuss the advantages and disadvantages between our methods and the previous ones  @xcite for generating arbitrarily entangled states of two microwave fields or two vibrational modes .
[ c]c|c|c|c|c & * pop . leak . * & * no . at . lev .
* & * st . eff . * & * mult . proc . * + ref .
@xcite & yes & 2 & no & l. pn .
@xcite & no & 3 & no & l. pn .
@xcite & no & 3 & no & h. pn .
@xcite & no & 2 & yes & none + ref .
@xcite & no & 2 & no & h. pn .
@xcite & no & 2 & yes & none + our proposal & no & 2 & no & l.
+    [ tab : compareuniversal ]    the brief comparison between these methods is listed in table .
[ tab : compareuniversal ] . in detail ,
@xcite provided an algorithm to generate arbitrarily entangled states of two vibrational modes .
but due to population leakage outside the original space , it takes an exponential complexity of the number of steps .
the succeeding proposals  @xcite overcome the exponential drawback in several ways : ( 1 ) a third atomic level is used to shield oscillations that cause population leakage  @xcite .
but the disadvantage is that higher energy levels of systems usually have larger decay rates , which inevitably reduce the fidelities of the target states .
( 2 ) boson - number - dependent stark effects are used to realize independent operations of particular states  @xcite .
but the disadvantage is that the detunings of nonresonant terms are usually less by one order of the coupling strengths between the two - level system and boson modes .
this means that the rabi frequencies are smaller , and the longer generation time is required .
( 3 ) multiphoton processes of high photon number are used to shield oscillations that cause population leakage or reduce the number of steps  @xcite .
but the disadvantage is that if the coupling strengths between the atom and cavity fields are not high enough , then the rabi frequencies become small , especially for states with high photon numbers , which obviously indicates longer generation time .    besides the advantage that there is no population leakage , our method has also the following advantages compared with previous ones .
@xcite : ( 1 ) it only uses the two energy levels of the qubit .
thus , the fidelities of the target states should be higher because there is no other auxiliary energy levels .
( 2 ) the detunings of the nonresonant terms are in the order of the resonator frequencies .
they are usually bigger than the coupling strengths between the qubit and resonator modes .
thus the rabi frequency can be made bigger than those using boson - number - dependent stark effects .
( 3 ) we use multiphoton processes of low photon number , i.e. , one photon at most is converted for either mode .
thus the rabi frequency can be bigger than those using multiphoton processes of higher photon number , especially when the coupling strengths between the qubit and cavity modes are not very big . of course
, stronger couplings will further enhance the rabi frequencies and hence reduce the generation time .
we point out that the real supercoducting qubit circuits are mutilevel systems , the information leakage to higher levels is not avoidable .
however , the leakage can be neglected when the transition frequency between the first excited state and the second excited state is much larger than the qubit frequency .
for example , in the flux qubit circuits , due to its large anharmonicity of energy levels , the information leakage is negligibly small .
however , for the transmon and phase qubit , the anharmonicity is very weak .
thus , the pulse should be carefully calibrated to avoid information leakage to higher levels .
the pulse calibration can be done as in ref .
@xcite .
we now compare the differences between our algorithm and other ones for generating noon states . ref .
@xcite uses mutliphoton processes to generate noon states . in superconducting systems
, this means a low generation efficiency if the lamb - dicke parameter is not sufficiently big .
@xcite uses synchronization technology to generate noon states , but the time duration for synchronization between two steps can be quite long and there exists inevitably information leakage . ref .
@xcite and its experimental realization  @xcite use two phase qubits with three active energy levels to generate noon states of two cavity modes .
the experimental setup is complex and the high energy levels of qubits will reduce the decoherence time . * * ref .
@xcite uses photon number - dependent stark effects to achieve independent operations .
thus the rabi frequency is smaller than the qubit - cavity coupling strengths .
@xcite requires that two qubits be initially prepared in a bell state and finally get decoupled from the qubits and cavity fields .
@xcite uses one qubit but still needs one additional level to shield unwanted resonances .
more recently , ref .
@xcite uses one qubit of four levels which resonantly interacts with two resonators simultaneously to speed up the generation process of noon states .    when applied to generating noon states , our algorithm has new features besides the common advantages for generating arbitrary two - mode photon states : ( 1 ) only carrier processes  @xcite and one - photon processes are used . in this case
, even though the coupling strengths between the qubit and cavity modes are small , large rabi frequencies can still be obtained .
( 2 ) the number of steps is reduced to linear dependence on the maximum photon number .
these advantages indicate less generation time and thus guarantee a higher efficiency than preceding methods .
now we discuss the experimental feasibility of our scheme .
table  [ tab : psig ] and table  [ tab : psinoon ] show that without pulse calibration , higher fidelities can be achieved at bigger * * lamb - dicke parameters @xmath298 and reduced driving frequencies @xmath65 .
these values are already in the ultrastrong regime . ref .
@xcite has reported ultrastrong couplings between three resonator modes and a flux qubit , where the lamb - dicke parameter @xmath298 can reach as high as @xmath348 . in the ultrastrong regime
, rabi frequencies can be made to approach the magnitude of @xmath64 , which usually ranges from @xmath259 to @xmath349 ghz .
the decay rates of the qubit and cavity fields are usually in the magnitude of megahertz .
thus the dissipation has small effect on the fidelities of target states . for singe - mode microwave fields , fock states with up to six photons  @xcite and fock state superpositions  @xcite
have been experimentally demonstrated using phase qubits .
the noon state up to @xmath200 photons has also been experimentally reported  @xcite .
we thus hope that our proposal is also experimentally feasible in the near future .
in summary , we have proposed an approach to generate arbitrary superpositions of photon states of two microwave fields in two separated cavities .
our method mainly depends on the coexistence of transverse and longitudinal couplings between the qubit and cavity fields . employing the longitudinal couplings
, we derive a hamiltonian which is similar to that of trapped ions interacting with two vibrational modes  @xcite . using four simple interaction hamiltonians derived from the longitudinal coupling , we design the state generation algorithm . our algorithm can be regarded as the improved version of that in  @xcite when the transverse and longitudinal couplings coexist in circuit qed systems .
but it has remedied the drawback that the number of steps exponentially depends on the maximal photon number , which is replaced by a quadratic dependence .
compared with previous ones with quadratic complexity , our algorithm does not require atomic energy levels higher than two  @xcite , boson - number - dependent stark effects  @xcite , or multiboson processes of high boson numbers  @xcite .
when applied to the generation of noon states , whose engineering has been extensively studied  @xcite , our algorithm needs only carrier and one - photon sideband transitions .
meanwhile , the number of steps only linearly depends on the maximum photon numbers .
in fact , these properties for generating noon states can be generalized to any states with a constant total photon number of both modes .
we have also discussed how to avoid the effect of unwanted terms on the generation of target state .
our numerical results show that fidelities above @xmath350 can be reached in the ultrastrong regime for the two - photon evenly - populated state and noon state when the environmental effect is neglected
. the generation time can be very short , in which case , the environment has small effect on fidelities of the target states .
we here note that due to the similarity of two - mode interaction hamiltonians , the algorithm using two - mode multi - phonon processes in ref .
@xcite can be directly applied into our model .
thus , two - mode fock states with high photon numbers can be generated with just two steps as one - mode fock states in  @xcite .
we have noted that our method for generating noon states is similar to a recent algorithm simplified from the one which employs stark effects to generate arbitrary entangled states  @xcite .
yxl is supported by the national basic research program of china grant no .
2014cb921401 , the nsfc grants no .  61025022 , and no .
the basic principle of our algorithm is state space reduction .
the state space of the target state in eq .
( [ eq : psi_f ] ) can be denoted  as @xmath351 we will implement @xmath352 procedures , with each procedure containing some steps .    in the 1st procedure , we aim to clear the populations in the subspace @xmath353 of @xmath354 .
this can be achieved via alternatively switching @xmath122 10  transitions and @xmath122 01  transitions , i.e. , @xmath355{0\bar{1}}\left\vert 0,n_{\max}-1\right\rangle \left\vert e\right\rangle \xrightarrow[f-1]{\bar{1}0}\left\vert 1,n_{\max}-1\right\rangle \left\vert g\right\rangle \cdots\nonumber\\ &   \xrightarrow[f-2n_{\mathrm{max}}+2]{0\bar{1}}\left\vert n_{\max } -1,0\right\rangle \left\vert e\right\rangle \xleftarrow[f-2n_{\mathrm{max}}+1]{\bar{1}0}\left\vert n_{\max},0\right\rangle \left\vert g\right\rangle .\end{aligned}\ ] ] to transfer the populations in the subspace @xmath356 of @xmath354 to the state @xmath357 . thus , @xmath354 is reduced to the state space @xmath358 where@xmath359      in the @xmath361th procedure , we aim to clear the populations in the subspace @xmath362 of @xmath363 .
this can be achieved via alternatively switching @xmath364 11  transitions and @xmath364 00  transitions , i.e.,@xmath365{1\bar{1}}\left\vert 1,n_{\max}-\mu-1\right\rangle \left\vert e\right\rangle \nonumber\\ &   \xrightarrow[f - n_{2\mu}-1]{00}\left\vert 1,n_{\max}-\mu-1\right\rangle \left\vert g\right\rangle \cdots\nonumber\\ &   \xrightarrow[f - n_{2\mu}-2(n_{\mathrm{max}}-\mu)+2]{1\bar{1}}\left\vert n_{\max}-\mu,0\right\rangle \left\vert e\right\rangle \nonumber\\ &   \xrightarrow[f - n_{2\mu}-2(n_{\mathrm{max}}-\mu)+1]{00}\left\vert n_{\max } -\mu,0\right\rangle \left\vert g\right\rangle .\end{aligned}\ ] ] with @xmath366 , to transfer the populations in the subspace @xmath367 of @xmath363 to the state @xmath368 . here
thus , @xmath363 is reduced to the state space @xmath369 where @xmath370    in the @xmath371th procedure , we aim to clear the populations in the subspace @xmath372 of @xmath373 .
this can be achieved via alternatively switching @xmath364 10  transitions and @xmath374 01  transitions , i.e. , @xmath375{\bar{1}0}\left\vert 1,n_{\max}-\mu-1\right\rangle \left\vert g\right\rangle \nonumber\\ &   \xrightarrow[f - n_{2\mu+1}-1]{0\bar{1}}\left\vert 1,n_{\max}-\mu -2\right\rangle \left\vert e\right\rangle \cdots\nonumber\\ &   \xrightarrow[f - n_{2\mu+1}-2(n_{\mathrm{max}}-\mu)+3]{\bar{1}0}\left\vert n_{\max}-\mu-1,0\right\rangle \left\vert e\right\rangle \nonumber\\ &   \xleftarrow[f - n_{2\mu+1}-2(n_{\mathrm{max}}-\mu)+2]{0\bar{1}}\left\vert n_{\max}-\mu,0\right\rangle \left\vert g\right\rangle .\end{aligned}\ ] ] with @xmath376 , to transfer the populations in the subspace @xmath377  of @xmath373 to the state @xmath378 .
thus , @xmath373 is reduced to the state space @xmath379 where@xmath380      in the ( @xmath352)th procedure , we aim to clear the populations in the subspace @xmath381 of @xmath382 .
this can be achieved via switching one 00  transitions , i.e.,@xmath383{00}\left\vert 0,0\right\rangle \left\vert g\right\rangle .\ ] ] to transfer the populations in the subspace @xmath384 of @xmath382 to the state @xmath385 .
thus , @xmath382 is reduced to the state space @xmath386 where @xmath387 is namely the intial state space for @xmath388 in eq .
( [ eq : psi_0 ] ) .
base on the above discussion , we now calculate the number of steps to generate the target state in eq .
( [ eq : psi_f ] ) with an arbitrary maximum photon number @xmath122 . for alternatively switching 11
`` and 00 '' transitions , we need @xmath389 steps , given by @xmath390 for alternatively switching 10 `` and 01 '' transitions , we need @xmath391 steps , given by @xmath392 therefore , the total step number @xmath139 is @xmath393 for generating the target state in eq .
( [ eq : psi_f ] ) .
h. wang , m. mariantoni , r. c. bialczak , m. lenander , e. lucero , m. neeley , a. d. oconnell , d. sank , m. weides , j. wenner , t. yamamoto , y. yin , j. zhao , j. m. martinis , and a. n. cleland , phys .
106 * , 060401 ( 2011 ) .