dileptons produced in heavy ion collisions interact only electromagnetically with their hadronic surroundings . especially the early phases of such collisions , when dilepton production rates are largest ,
can therefore be explored with lepton pairs .
some of the initial ideas on photon @xcite and dilepton @xcite @xcite emission have met with considerable interest , and this area has now expanded into a major component of research in ultrarelativistic heavy ion collisions @xcite . on the experimental side , dilepton experiments for hadronic collisions
@xcite@xcite have been complemented only recently by nucleus - nucleus experiments at very high energies @xcite@xcite .
since the spectrum of lepton pairs and of real photons has been proposed as a signature of the quark - gluon plasma @xcite , it is imperative quantitatively to understand the emission mechanisms from the confined sector of qcd , which in this sense constitutes the background .    guided by theoretical estimates
, we can roughly divide the lepton pair production cross section , expressed in terms of its invariant mass @xmath1 , into a number of distinct regions . for @xmath2 ,
the drell - yan contribution dominates , while pairs arising from the semileptonic decay of @xmath3 are also important @xcite . for lepton
pair invariant masses around the light vector meson masses ( @xmath4 , @xmath5 , and @xmath6 ) , two - body @xcite and three - body @xcite reactions dominate .
the latter processes also extend their influence to the region @xmath7 , where the quark - gluon plasma may become visible for sufficiently high initial temperatures @xcite . in this paper
, we concentrate on the regions of low invariant masses , the so - called `` soft '' limit @xmath8 mev . there ,
production via bremsstrahlung is expected to play an important role and many authors have recently concerned themselves with such reactions @xcite@xcite .
since particle production plays an important role in ultrarelativistic heavy ion collisions , we consider only events that arise from microscopic meson - meson dynamics . in high energy
heavy ion experiments , the meson to baryon ratio is such that the absence of baryons in our treatment is not expected to constitute a major hurdle to phenomenology . more specifically , we restrict our studies to those of pion - pion bremsstrahlung .
the production of photons or dileptons necessarily involves 3- or 4-particle final states . except for special cases such as cross sections at fixed @xmath9 @xcite ,
such final states have universally been handled under some approximation .
it is our goal in this paper to go beyond these by presenting an essentially exact formulation for the @xmath10 cross section ( as a function of dilepton invariant mass ) and , in the process , to test how far these approximations may be trusted . in order to achieve these goals , we formulate the pion - pion interaction in terms of a relativistic one boson exchange model ( obe ) .    note that we do not attempt to compare directly to experimental data because the main point at issue is a comparison of theoretical approaches to bremsstrahlung . to do justice to small-@xmath1 data
would require spacetime integration over the fireball region as well as inclusion of dalitz decays and the landau - pomeranchuk effect , a substantial task which will be postponed to future work .
our paper is organized as follows : in the next section , we discuss the set of approximations known under the collective banner of the soft photon aproximation ( spa ) .
next , we derive exact formulae for the bremsstrahlung generation of lepton pairs .
section [ sec : threesp ] is devoted to technical issues having to do with phase space ; sections [ sec : curint ] and [ sec : xsexact ] then derive cross sections based respectively on various spa variants and the exact formulation .
our obe model is presented in section [ sec : fullobe ] ; finally , section [ sec : results ] contains our numerical results in the form of cross sections @xmath11 and rates @xmath12 and a discussion on the issues raised .
the soft photon approximation ( spa ) has been widely used in calculating bremsstrahlung dilepton spectra .
it is based on the early realization @xcite that the cross section for production of low - energy real photons is dominated by the corresponding hadronic amplitude . in practice , this means neglecting the photon momentum in the hadronic matrix element as well as all photon emission from vertices and internal lines ( see section [ sec : spa ] ) .      in order for the spa to be a valid approximation , two conditions @xcite
must be met :    1 .   the photon energy must be much smaller than the energy of any one of the hadrons participating in the scattering , @xmath13 2 .
the process of radiating the photon must be separable from the hadronic interaction process , i.e.  the hadronic and electromagnetic time and distance scales must be sufficiently different to permit separate treatment .
the range of the hadron - hadron interaction is given roughly by @xmath14 , the inverse mass of the exchange boson , while the distance the hadron may propagate with off - shell energy @xmath15 before emitting the photon is about @xmath16 .
the spa is valid only when @xmath17 , or @xmath18    implicit in these equations is , of course , a specific lorentz frame with respect to which the energies are measured .
for the special case where two hadrons of equal mass @xmath19 collide , these conditions can be re - written in their cms frame as @xmath20 in a simple bremsstrahlung experiment , these limits are easily satisfied by selecting only photons or dileptons of low energy in the laboratory frame .
an early dilepton paper by rckl , for example @xcite , calculated spectra for @xmath21 and @xmath22 gev while restricting transverse momentum to below 500 mev .    in the complex multiparticle systems formed in the course of nucleus - nucleus collisions ,
however , the situation is more complicated : there are many binary collisions , and their respective cms frames do not generally coincide either with one another or with the overall nucleus - nucleus cms frame . a photon that is soft in a particular hadron - hadron cms therefore does not have to be soft in the laboratory and vice versa .
( this fact has to be taken into account when addressing issues related to the spectrum of real photons @xcite . )    in a situation where there are many cms frames , it therefore becomes necessary to look at relativistically invariant quantities .
for lepton pairs , the extra degree of freedom provided by a non - vanishing invariant mass is a suitable variable of choice . looking at invariant masses means , however , that @xmath23 is no longer fixed but must vary over its full kinematic range , which for our example of colliding equal - mass hadrons is given by ( see section [ sec : curint ] ) @xmath24 once the upper kinematic limit for @xmath23 approaches or exceeds the bounds set by either eq .
( [ iad ] ) or ( [ iae ] ) , the soft photon approximation clearly fails .    in figure 1
, we show the three functions ( [ iad ] ) and ( [ iae ] ) and ( [ iah ] ) for the case @xmath25 mev , @xmath26 mev and dilepton invariant masses @xmath27 and @xmath28 mev .
it is immediately clear that the assumptions underlying the spa are not fulfilled even for small @xmath1 : the kinematic range accessible to @xmath23 is never much smaller than the limits set by the spa .
the situation becomes even worse for larger @xmath1 .
closely related to the above is a common misconception that the spa must necessarily be valid for small invariant masses .
the confusion stems from the fact that , because @xmath29 a small value for @xmath23 necessarily implies small @xmath1 .
however , the reverse is not true : both @xmath23 and @xmath9 can be large even while adding up to small @xmath1 .
an extreme counterexample is of course a hard real photon with zero invariant mass but large energy .
one can not , therefore , a priori expect spa calculations to coincide with full hadronic calculations even when @xmath1 is small .
section [ sec : results ] provides examples of such deviations ; see especially the ratios between approximations and the exact calculation in figures 12 and 13 which do not approach unity for small @xmath1 .
this point has been emphasized before @xcite .
we are therefore led to conclude that the spa does not generally apply to rates or cross sections measured in terms of the dilepton invariant mass .
correspondingly , doubt must be cast on the validity of the spa in calculations pertaining to such measurements .
that said , one may still ask _ by how much _ spa - based calculations differ from the more complicated true situation : is it a matter of a few percent , or orders of magnitude ?
full calculations such as presented in section [ sec : fullobe ] below are cumbersome and time - consuming , and it is therefore of value to understand quantitatively to what percentage the spa approximation can be trusted .
it is this _ quantitative _ question that we attempt to cast some light on by making a model calculation capable both of `` true '' answers and an spa approximation .
clearly , a single model calculation such as presented below can not answer this quantitative question exhaustively .
nevertheless , its results provide at least an indicator of the reliability of the spa : if , within our simple obe , dileptons can or can not reasonably be described within the spa , it is possible or even likely that the same or similar conclusions would emerge from other models .      rather than being just a single step , the soft photon approximation as implemented in a practical context can consist of a number of stages . starting with the most fundamental ,
these are :    1 .
the virtual photon producing the lepton pair is emitted only by external legs of the hadronic reaction ; emission from internal parts of the hadronic reaction is neglected because the radiation from internal lines forms a sub - leading contribution @xcite .
2 .   the dependence of the hadronic matrix element @xmath30 on the photon momentum @xmath31 is neglected .
3 .   the photon momentum @xmath31 is neglected in the phase space delta functions .
in addition , we shall be considering in this paper the effect of    1 .   approximating the virtual photon current by the real photon current , ( @xmath32 ;
see eqs .
( [ brk ] ) and ( [ brl ] ) ) .
\(a ) and ( b ) have formed the basis for most bremsstrahlung dilepton calculations as they are hard to avoid .
the restriction of phase space ( c ) , on the other hand , can be circumvented fairly easily by inserting later a phase space correction factor ( section [ sec : ruckl ] ) or by avoiding ( c ) altogether using eqs .
( [ brv ] ) , ( [ fcm ] ) and sections [ sec : threesp ] and [ sec : xsexact ] .    the change from virtual to real photon currents in step ( d ) is not directly dependent on ( c ) and can be implemented either jointly or separately .
we shall be investigating both alternatives and determining their impact on the overall dilepton rate . in section
[ sec : fullcross ] and section [ sec : fullobe ] and beyond , we shall refrain from using even the basic assumptions ( a ) and ( b ) and calculate exact cross sections and rates within the same lagrangian .
we wish to find the cross section for a dilepton pair @xmath33 emitted in the semi - elastic scattering of two pions @xmath34 and @xmath35 into @xmath36 and @xmath37 , @xmath38 the lepton pair forms via a virtual photon @xmath39 with four - momentum @xmath40 and mass @xmath41 emitted from either the central blob or external legs of the hadronic reaction @xmath42 the leading contribution will come from photon emission from external pion legs .
first , we follow lichard @xcite in a leading - terms derivation but place more emphasis on the role of phase space ; the full formalism supplanting the leading - terms derivation is presented in section [ sec : fullcross ] .
four diagrams contribute to the leading - term cross section , one of which is shown in figure 2 . in the notation of figure 2
, the total matrix element can be written as products of purely hadronic reactions @xmath30 with currents @xmath43 and a leptonic part @xmath44 @xcite , @xmath45\ ;    l_\mu \;,\end{aligned}\ ] ] where @xmath46 and , in slightly condensed notation , the terms of the pionic current are @xmath47 with @xmath48 the charge of pion @xmath49 in units of the proton charge . note
that @xmath43 contains the propagator for the off - shell pion . making approximation ( b ) ,
one now neglects the dependence on @xmath31 of the four hadronic matrix elements ; for example @xmath50 permitting the usual factorization of @xmath51 into a purely hadronic matrix element and a photon - plus - lepton part : @xmath52 where now the virtual - photon current @xmath53 contains the sum over the four diagrams . since approximations ( a ) and ( b ) are valid only for photons of small @xmath31 and by implication for small @xmath1 , the extra terms in numerator and denominator are usually neglected ( in what we called approximation ( d ) ) to form the real - photon current @xmath54 squaring the total matrix element and summing over lepton polarizations @xmath55 yields , using either current @xcite , @xmath56 \,,\end{aligned}\ ] ] where @xmath57 and the leptonic tensor is @xmath58 at this point , @xmath31 is still defined only as the sum of lepton momenta @xmath59 .    with @xmath60^{1/2}$ ] the incoming flux ,
the unpolarized cross section is an integral over phase space of _ four _ outgoing particles , @xmath61 where we write @xmath62 for short .
this can be reduced to a _
three_-particle phase space integral by changing variables and integrating over redundant degrees of freedom in the dilepton cms , @xmath63 where the positron solid angle @xmath64 must be kept due to the @xmath65 term , and @xmath66 is the lepton mass . the dependence of @xmath65 on @xmath64 is integrated out , @xmath67 d\omega_+ = { 8\pi\over 3}\left(1 + { 2\mu^2\over m^2}\right ) ( -j^2 ) \,,\ ] ] so that , using eqs .
[ brm][brs ] , we can express the total cross section in terms of a leptonic factor @xmath68 and a `` hadronic plus virtual photon '' cross section @xmath69 , @xmath70 with @xmath71 and where @xmath72 is a function of _ three_-particle phase space @xmath73 with the understanding that its derivation is valid under approximations ( a ) and ( b ) only , eq .  ( [ brt ] ) is exact and makes no assumptions with respect to phase space . in particular , we note that @xmath31 is still contained in the delta function constraining the phase space in eq .
( [ brw ] ) .
next , we shall show that eq .
( [ brt ] ) is valid even when approximations ( a ) and ( b ) are not made but that @xmath69 takes on a form different from eq .
( [ brv ] ) .      when calculating dilepton cross sections for all possible diagrams , including internal emission of photons , it is no longer possible to factorize the matrix element into a hadronic part and a current as in eq .
( [ bre ] ) .
instead , we consider , for a given pion - pion reaction , the full matrix element @xmath74 where @xmath19 runs over all contributing diagrams .
the sum of diagrams is gauge invariant , @xmath75 but not necessarily the individual terms . unlike in the previous sections
@xmath76 contains , besides the hadronic interaction , the photon vertex and where applicable the external pion propagator , both of which were previously part of the current .    starting again with the cross section in four - phase space ( [ bro ] ) and transforming from @xmath77 to @xmath78 using eq .
( [ brqq ] ) , the unpolarized dilepton cross section is @xmath79 squaring and summing over lepton polarizations , we obtain in analogy to eq .
( [ brm ] ) @xmath80 where as before @xmath81 . in the dilepton cms , @xmath82 and
the second term in @xmath83 integrates to @xmath84 since @xmath85 in this frame , @xmath86 and hence in the dilepton cms @xmath87 so , using gauge invariance ( [ fcc ] ) , @xmath88 hence @xmath89 and the full cross section becomes @xmath90 which is identical to eq .
( [ brt ] ) , but with the virtual photon cross section reading @xmath91 comparing to eq .
( [ brv ] ) , we have the correspondence @xmath92 @xmath93 .    equations ( [ fcl ] ) and ( [ fcm ] ) , together with use of the 3-phase space variables of section [ sec : threesp ] provide a fully covariant formalism capable of exact treatment of both hadronic and electromagnetic sectors of a given model . within the physical limitations implicit in a given model , it makes no further assumptions or approximations .
in particular , none of the approximations listed in section [ sec : spa ] are needed .      in the context of nuclear collisions , one
is interested more in overall dilepton rates than in fixed-@xmath94 cross sections . not attempting to account fully for either quantum statistical or flow effects here
, we shall make the simplest assumption of a locally thermal boltzmann gas of pions at a ( local ) temperature @xmath95 .
the rate of pair production of invariant mass @xmath1 per unit of four - volume then is @xcite @xmath96 where @xmath97 is the basic three - particle kinematic function @xcite , which for equal masses @xmath98 has the more familiar form @xmath99 , and @xmath100 is the spin degeneracy factor which for pions is unity .
@xmath101 is the modified bessel function .
to complete our preparatory work , we summarize the treatment of three - particle phase space in terms of relativistic invariants .
the relations of this section will be useful first in 3-phase space integration of the currents ( [ brk ] ) and ( [ brl ] ) and later in the exact treatment of bremsstrahlung dilepton production .    three - particle phase space has been studied extensively and we here merely outline the procedure ; for details see ref .  @xcite .
as shown in figure 3 , the kinematics for a reaction @xmath102 can be described in terms of the five invariants @xmath103 , @xmath104 , @xmath105 , @xmath106 , and @xmath107 . for our purposes
, we shall be using @xmath108 and @xmath31 interchangeably as referring to the virtual photon , i.e.  @xmath109 and @xmath110 .
the other momenta refer , as usual , to the incoming and outgoing hadrons .
inverting , one gets @xmath111 the corresponding phase space integral is @xmath112 where the weighting is given in terms of the cayley determinant @xmath113 the cayley determinant is quadratic in any of its arguments ; specifically , we use the form @xmath114 where the kinematic limits on @xmath115 are given by @xmath116^{1/2 } , \end{aligned}\ ] ] written in terms of yet another determinant and the basic four - particle kinematic function @xmath117 the kinematic limits on @xmath118 are @xmath119 the limits on @xmath120 and @xmath121 are given by the chew - low plot : for values of @xmath121 bounded by @xmath122[s + m_1 ^ 2 - ( m_2 + m)^2 ] \nonumber\\ & \pm & { 1\over 2s }      \lambda^{1/2}(s , m_a^2,m_b^2)\ ;      \lambda^{1/2}(s,(m_2{+}m)^2,m_1 ^ 2 ) \,,\end{aligned}\ ] ] the limits on @xmath120 are @xmath123[m_a^2 + m_1 ^ 2 - t_1 ] \nonumber\\     & & \ \ \ \   + { 1\over 2m_a^2 }      \lambda^{1/2}(s , m_a^2,m_b^2)\ ;      \lambda^{1/2}(t_1,m_a^2,m_1 ^ 2 ) \ , , \\
\label{fvkm } s_2 ^ - & = & ( m_2 + m)^2 \,.\end{aligned}\ ] ] although not directly relevant for pions , we include for completeness the case where the hadronic masses @xmath124 and @xmath125 are unequal . in this case , @xmath121 may become larger than @xmath126 given above , and the corresponding additional phase space is delimited by @xmath127 and for this additional domain , @xmath120 is limited by @xmath128 of eq .
( [ fvk ] ) and the corresponding lower branch , given by a minus sign before the @xmath129-factors .
clearly , @xmath121 resembles the 2-phase space invariant @xmath130 , but they are distinct quantities due to the presence of the additional momentum @xmath31 .
integration over pure phase space yields @xmath131 where @xmath132 and we have written @xmath133 for short .
this yields identical results to the pure phase space formula ( [ icg ] ) below , as it should .
we note in passing that the ratio of the kinematic domains of ( 3-space ) @xmath121 to ( 2-space ) @xmath130 is exactly the phase space reduction factor found in ( [ icj ] ) below .
all the elements for our calculation are now in place . in this section ,
we look at three approximations of increasing sophistication .
we proceed first to find explicit expressions for the pion - pion - photon cross section @xmath69 of eq .
( [ brv ] ) based on the factorization into a current and hadronic matrix element of section [ sec : external ] .
these approximations are based on an angular average of @xmath134 in the @xmath135 center - of - momentum frame . using three - space invariants ,
we integrate the current covariantly in section [ sec : xs3sp ] .      in 1976 , rckl proposed @xcite that the dilepton cross section be written in terms of the _ real photon _
cross section : @xmath136 this ansatz , now in common use , was utilized by haglin et al . @xcite and others previously @xcite to derive dilepton production rates . in our notation
, the rckl formula reads @xmath137 which becomes , using the formulas of the previous section , and approximating @xmath138 and @xmath139 , @xmath140 the real - photon cross section @xmath141 can be written down immediately from the virtual - photon cross section @xmath69 of ( [ brv ] ) , @xmath142 neglecting @xmath31 in the delta function of @xmath143 , the 3-phase space integral of eq .
( [ brw ] ) factorizes into the usual 2-phase space and the @xmath9-integration @xmath144 with @xmath145 leading all in all to @xmath146 where @xmath147 is the elastic hadronic cross section . integrating @xmath134 first with respect to the photon phase space with @xmath148 leads to @xmath149 where @xmath150/2\sqrt{s}$ ] is the upper kinematic boundary .
the angular average over @xmath151 was found by haglin et al .  for the real - photon current ( [ brl ] ) in terms of velocities @xmath152 in the @xmath153 ) cms as @xcite @xmath154 \,,\end{aligned}\ ] ] where @xmath155 \ , , \\ d_{ij } & \equiv &         \left [    ( \bbox{\beta}_i -         \bbox{\beta}_j)^2                - ( \bbox{\beta}_i { \times } \bbox{\beta}_j)^2         \right]^{1/2 } \,.\end{aligned}\ ] ] to convert these angular averages into functions of 2-phase space invariants @xmath156
, it is necessary to invoke approximation ( c ) restricting phase space : while @xmath157 in the chosen frame , the presence of @xmath9 means that @xmath158 , i.e.  @xmath159 and @xmath160 are not true cms velocities . ignoring therefore @xmath9 at this point , the velocities and hence the angular averages can be written @xcite in terms of @xmath94 , @xmath130 and @xmath23 , @xmath161 \,,\end{aligned}\ ] ] with @xmath162 as usual .
having expediently eliminated @xmath31 at two points , it is very important in this formulation to correct for phase space , i.e.  to alleviate the effects of approximation ( c ) .
pure phase space integration of @xmath143 ( _ including _ @xmath31 in the delta function ) yields exactly @xmath163 where @xmath164 integration over pure two - particle phase space , on the other hand , yields @xmath165 so @xmath166 , unconstrained by @xmath31 , grossly overestimates @xmath167 , making it necessary to insert _ a posteriori _ the factor @xcite @xmath168 into the @xmath23-integration ( [ rrh ] ) .
the `` rckl approximation '' for dilepton production therefore reads , in its final form , @xmath169 this is the form used in ref .
@xcite .
starting not from the rckl formula but rather from lichard s improved formalism @xcite of section [ sec : external ] , the derivation proceeds along the same lines as the above .
one finds correspondingly for the real - photon current @xmath170 differing from the above by a constant factor @xmath171 and the extra factors entering the exact version ( [ bru ] ) of @xmath68 .    the corresponding virtual - photon current yields the third approximation within the angle - averaged current family
, @xmath172 where the angular average of ( [ brk ] ) is @xmath173 where @xmath174 ) reduces to eq .
( [ hgk ] ) for @xmath175 .
we shall be testing the effect of using the virtual over real - photon currents within the present approximation in section [ sec : results ] below .
we note in passing that eq .
( [ hgh ] ) would be preferable to eq .
( [ hgk ] ) for very small hadron masses @xmath19 , since then the latter is logarithmically sensitive to @xmath176 while the former is regulated by the dilepton invariant mass @xmath1 .
three approximations were needed for the above angle - averaged current cross sections : the rckl formula , converting cms velocities @xmath152 to invariants @xmath94 and @xmath130 leading to eqs .
( [ hge ] ) and ( [ hgf ] ) , and correcting for phase space using eq .
( [ icj ] ) . at the cost of increased computing time ,
all three can be avoided by utilizing the formalism of three - particle phase space invariants .
within 3-phase space , and writing the flux as @xmath177 , the bremsstrahlung cross section of eq.([brv ] ) becomes ) be possible to restore @xmath31 to @xmath30 and then to include it with @xmath178 in the inner integral , but because of missing contact terms this procedure is not gauge invariant and shall not be pursued . ]
@xmath179   \,.\ ] ] this differs from eq .
( [ icklo ] ) in the way the current is treated : by an angular average plus @xmath23-integration for the earlier , by full 3-phase space invariants in the present case .
note that @xmath30 is here a function of @xmath121 rather than @xmath130 .
in all the above approximations , neglecting the dependence on @xmath31 of the squared hadronic matrix element @xmath180 is unavoidable because this had been a precondition to factorizing into electromagnetic and hadronic parts in eq .
( [ bre ] ) .
there have been attempts go beyond this on - shell matrix element approximation by expanding @xmath30 in @xmath31 ( see refs .
@xcite ) , but this procedure involves successively more complicated derivatives without guaranteeing convergence .    using 3-phase space invariants and the exact for the emission of real photons and to order @xmath181 for dilepton emission . ] expressions ( [ fcm ] ) for the pion - pion - photon cross section and ( [ fcl ] ) for the dilepton cross section , such expansions in @xmath31 become superfluous . to calculate explicit dilepton cross sections , we merely insert from our favorite model for pion - pion scattering the complete set of matrix elements , and write the elements of the squared sum of ( [ fcm ] ) in terms of these invariants :
@xmath182 this casting in terms of invariants means that @xmath31 and all consequent off - shell effects in @xmath183 are fully taken care of . can be implemented within 2-phase space when only the _ differential _ cross section @xmath184 need be calculated ; see e.g.  ref.@xcite ; for @xmath185 , this becomes impossible . ]
putting it all together , we obtain @xmath186   \,.\ ] ] while the cross sections may be exact , microscopic models for the pion - pion interaction are far from perfect . in the next section , we make use of a simple model to illustrate the present cross section and the approximations enumerated before .
the main point of this section , however , is to stress that , apart from model - induced ailments , the exact treatment of any model for the reaction @xmath187 is possible .
to compare dilepton cross sections under the exact formula and various approximations introduced so far , we must make a model : little if any experimental data for the differential pion - pion - photon cross section is available . to this purpose , we return to a simple lagrangian with @xmath188 , @xmath4 and @xmath189 exchange used previously @xcite which , when @xmath130-integrated , successfully reproduced experimental data for elastic @xmath190 scattering @xcite : @xmath191 here , we include the full momentum dependence of the @xmath192 propagator which now reads @xcite @xmath193 the @xmath194 field also has a @xmath195-dependent term in the unitary gauge propagator @xcite , @xmath196 \,,\ ] ] but this term contributes to matrix elements only in the form @xmath197 etc .  and is thus zero for equal pion masses .    for @xmath130- and @xmath198-channel exchange of an @xmath199-meson ( @xmath200 ) with momentum @xmath195 by the pions ,
we implement monopole strong form factors @xcite @xmath201 in figure 4 , we show the results of fitting this model to the elastic scattering data as before @xcite but with the different propagator ( [ rtd ] ) .
again , the fit is reasonably good , but the strong dependence on momentum - dependent terms in the @xmath192 propagator shows up in the region above the @xmath192 peak .
best - fit values for the free parameters now are : @xmath202 mev , @xmath203 mev , @xmath204 mev ; @xmath205 , @xmath206 , @xmath207 . throughout
, we take @xmath208 mev .    for the current - based approximations of section [ sec : curint ] , the above lagrangian is all we need ; for the complete cross section of section [ sec : xsexact ] , it must be augmented by internal photon emission pieces
. this can be achieved by minimal substitution , thereby guaranteeing gauge invariance ( for a similar approach , see the `` quantum hadrodynamics '' lagrangians of ref .
@xcite ) . to the hadronic lagrangian ( [ rtb ] ) we therefore add @xmath209 details regarding these pieces and corresponding vertex factors can be found in the appendix .
the @xmath210 interaction results in a modified strong form factor for @xmath130- or @xmath198-channel charged @xmath4 exchange . if @xmath211 and @xmath212 are the usual form factors at the @xmath213 and @xmath214 vertices respectively , with @xmath215 , electromagnetic gauge invariance requires that the form factor for the emission of a photon by the @xmath4 be @xmath216 where @xmath217 and @xmath218 are the respective @xmath4 propagators . for the monopole form ( [ rth ] ) of
@xmath219 this becomes @xcite @xmath220 note that we do not include strict vector meson dominance coupling of the photons to vector mesons .
this is equivalent to the introduction of electromagnetic form factors .
for the specific application we are considering , small @xmath1 s , those would be @xmath221 1 . for previous obe models ,
see e.g.  refs.@xcite and references therein .
to summarize : the differential cross section for dilepton production as a function of invariant mass , @xmath222 , can be found in a number of approximations and an exact way ( within the limitations of the hadronic reaction model ) .
the approximations are : the rckl approximation ( [ ick ] ) , two approximations , eqs .  ( [ icklz ] ) and ( [ icklo ] ) , based on a more careful derivation by lichard , using the real and virtual photon currents respectively , and an approximation integrating the current covariantly in 3-phase space , eq .  ( [ bvp ] ) .
the `` exact '' formula is given by eq .
( [ xxc ] ) .    of these
, the rckl and lichard cross sections can be found easily and quickly within 2-phase space .
the 3-phase space current approximation is somewhat more difficult but still uses 2-phase space hadronic matrix elements . implementing the exact formulation ( [ xxc ] ) involves much more effort : if there are @xmath223 diagrams contributing to a given reaction , the sum in ( [ xxc ] ) contains @xmath224 terms .
since an individual term @xmath225 in ( [ xxb ] ) is itself a dirac sum , only partial factorization of the overall sum is possible .
further complications arise from the presence of imaginary pieces in the matrix elements .    to quantify the differences between the approximations and the exact formulation , we have studied all five distinct pion - pion reactions .
writing @xmath226 as shorthand for the reaction @xmath227 and so on , we have calculated within our obe model cross sections for @xmath226 , @xmath228 , @xmath229 , @xmath230 and @xmath231 .
a total of 36 diagrams contribute at tree level to the first two reactions ( 24 for @xmath39 emitted by external pion lines , 12 by vertices and the exchange meson ) ; the remaining reactions are made up of 16 diagrams each ( 8 external , 8 internal ) .
numerical results were checked by performing the following consistency checks : setting @xmath232 , @xmath233 , @xmath234 , @xmath235 and @xmath236 in the _ hadronic _ part of the exact matrix elements ( while keeping these variables where they enter the `` current '' part ) and using only the external - emission diagrams reproduces exactly the results of the 3-phase space integrated current ( [ bvp ] ) .
gauge invariance in the @xmath188 , @xmath4 and @xmath192 fields provided another sensitive test .
figures 59 show the cross sections for these five reactions as a function of @xmath237 . in every case ,
the left and right panels show @xmath238 for @xmath27 mev and 300 mev respectively .
final - state symmetrization factors were included where appropriate .
initial - state symmetrization was also included for @xmath228 and @xmath230 in order to facilitate use within a thermal pion gas environment .
when the initial - state pions are considered distinguishable , these two cross sections should be multiplied by a factor 2 for these reactions .
because they are identical in structure to their charge - conjugate versions , cross sections for the reactions @xmath231 and @xmath228 were doubled  this , too , should be corrected for when desired .    to prevent overcrowding , only the rckl plus the virtual-@xmath39 current approximations ( [ icklo ] ) and ( [ bvp ] )
are shown , together with the exact calculations .
all these were computed using the same obe model and parameter values specified in section [ sec : fullobe ] above .    in the complicated structures and deviations in figs .
59 , the following points are of interest :    1 .
the @xmath4 peak in the reaction @xmath226 is overestimated by a factor 2.5 and 1.5 by the rckl and lichard approximations respectively for @xmath27 mev ; for @xmath239 mev , the overestimation is much greater ( figure 5 ) .
only the 3-space current approximation does an adequate job here .
similarly , the @xmath192 peak of the exact result lies well below the corresponding approximations .
overestimation ( 1.75.0 ) of the @xmath4 peak also occurs for approximations in @xmath231 .
2 .   in the reaction @xmath231 , the 3-space current approximation ,
on the other hand , lies well below the exact result ( figure 6 ) .
overestimation factors for @xmath238 are of order 0.54 for the reactions @xmath230 and @xmath229 , depending on approximation and cms energy . here , too , the 3-space current approximation generally does better than the others but again tends to underestimate in some parts ( figs .  7 , 8) .
the largest discrepancy between approximations and the exact result occur for the reaction @xmath228 ( figure 9 ) : for the rckl approximation , factors 3 ( for 10 mev ) to 30 ( for 300 mev ) arise , while the lichard approximation yields corresponding overestimation factors of 1.9 and 1420 .
since , however , this reaction contributes only little when contributions from all reactions are added up , this effect is not of much interest in a heavy - ion context . 5 .
there is generally not much difference between using the real - photon ( not shown ) or virtual - photon currents in the lichard approximations ; this is in agreement with the conclusions of ref .
@xcite . for @xmath239 mev
, some differences can be discerned , but both deviate in general more from the exact result than from each other
.    figures 10 and 11 show the dilepton production rates per unit spacetime summed over all seven reactions , calculated using the boltzmann formula ( [ rtk ] ) , for temperatures @xmath240 and 200 mev respectively .
these temperatures probably lie below and above an expected transition to a new quark - gluon phase and so represent extrema in their behavior . corresponding ratios of approximation over exact rates are shown in figures 12 and 13 . the rckl approximation is the worst , as expected ; overestimation is in the range 24 for @xmath240 mev , and 28 for @xmath241 mev .
the lichard approximations overestimate by factors 1.44 , depending on temperature and @xmath94 . here ,
the _ real-@xmath242 _
3-phase space current , the lower of the two dotted lines , does surprisingly well , deviating from the exact result by less than 20% throughout .
comparing this fact , however , with the substantial deviations shown in the @xmath94-differential cross sections of figs .
59 , one must conclude that much of the `` agreement '' seen on the level of rates is due to averaging over differences seen in the @xmath94-dependent cross sections and that such `` agreement '' therefore does not validate the approximations .
we also note the fact that none of the approximations approaches the exact result for small values of @xmath1 : even for the smallest value shown ( @xmath27 mev ) , the discrepancy is still above 40% for the lichard approximations and larger than a factor 2 for the rckl approximation .
figure 14 shows the fractional contributions to the total @xmath243 mev rate for the rckl , virtual-@xmath242 current lichard , and exact cross sections .
as discussed , the @xmath228 contribution virtually falls away in the exact calculation ; by contrast the reaction @xmath231 plus charge conjugate becomes more important , in second place behind the dominant reaction @xmath226 .    finally , we ask where the discrepancies between approximations and exact cross section originate : is this due to the inclusion of diagrams with photon emission from vertices and propagators ?
a partial answer is provided by figure 15 , where the cross sections for reactions @xmath244 and @xmath229 are plotted .
the upper lines correspond , as before , to the rckl and two lichard approximations , while the solid line again represents the exact result .
the lowest dash - dotted line , on the other hand , represents the exact result but excluding all internal diagrams and their cross terms with external ones .
the difference between this lower line and the exact result ( solid line ) therefore represents the contribution of the internal diagrams ; while the difference between the lower dash - dotted line and the upper lines ( approximations ) represents the change in cross section due to inclusion / exclusion of @xmath31 in the _ external_-emission diagrams .
we see that the contribution of internal emission is not all that large , albeit nonnegligible . by far
the most important effect on @xmath238 is the inclusion of the full dependence of @xmath245 on the photon momentum @xmath31 .
in fact , the effect is so large that for the reaction @xmath228 the cross section would become negative if internal diagrams were not included ! in other words , among the list of approximations listed in section [ sec : spa ] , approximation ( b ) is the most far - reaching .    the importance of including @xmath31 fully in the matrix element can be understood qualitatively : as illustrated in figs .
59 for one reaction , the elastic pion - pion cross sections vary considerably within the region of interest . adding or subtracting a sizeable amount of energy - momentum in the form of @xmath31
therefore naturally leads to a substantial change in cross section also .
this change is , however , neglected in all the approximations listed : there are no resonances in the cross section which could provide some amount of cancellation by pion - pion reactions being `` shifted into '' and `` shifted out of '' a resonance peak . ] .
we believe that this ( unavoidable ) neglect is at the heart of the considerable differences between approximations and exact result seen in figs .  59 .    in summary ,
the various approximations for bremsstrahlung from pion - pion collisions can not in general be believed beyond at best a factor 2 or more , depending on invariant mass and temperature .
we have shown in sections [ sec : fullcross ] and [ sec : xsexact ] how 3-phase space can be fully taken into account .
the most important effect , it turns out , is not the inclusion of photon emission from vertices or propagators but the strong dependence of the cross section on the inclusion of the photon momentum @xmath31 into the _ hadronic _
collision part .
since this calculation neglected baryonic degrees of freedom and resonances , and because of the technical difficulties of a further integration over spacetime , our results are not applied directly to experimental data ; this is a matter for the future .
all we can say at this stage is that it is likely that the approximations overestimate the dilepton yield throughout the mass region of interest , and that the degree of overestimation rises both with increasing invariant mass @xmath1 and temperature @xmath95 . since the bremsstrahlung channel competes with various dalitz decays , and since various experiments take much care in applying cuts to eliminate dalitz pairs as far as possible @xcite a quantitative estimate would necessarily have to take such experimental cuts into account also .
in addition , attention should be paid to the issue of the landau - pomeranchuk effect , which is expected to suppress dilepton bremsstrahlung rates @xcite .
our results , taken in conjunction with these calculations , therefore appear to indicate even larger suppression of bremsstrahlung dileptons than had been thought previously .
we thank p.  lichard for enlightening discussions and his careful derivation of dilepton formulas .
thanks to a.  drees , j.  zhang , j.  c.  pan , p.  lipa and b.  buschbeck for useful discussions .
this work was supported in part by the natural sciences and engineering research council of canada , by the qubec fcar fund , by a nato collaborative research grant , by the national science foundation under grant number 9403666 and by the austrian fonds zur frderung der wissenschaftlichen forschung ( fwf ) .
besides the hadronic interaction lagrangian ( [ rtb ] ) with its vertex factors @xmath246 we need the electromagnetic interaction lagrangian linear in @xmath247 , @xmath248 with the pion current @xmath249 the @xmath194 field tensor , @xmath250 and @xmath251 the usual electromagnetic tensor , the individual contributions are @xmath252_3     + { 1\over 2 }       e f_{\mu\nu } [ \bbox{\rho}^\mu \bbox{{\times } \rho}^\nu]_3 \\ & = & ie a_\mu \left [      \rho_{- \nu } (   \partial^\mu\rho_+^\nu                    - \partial^\nu\rho_+^\mu )    - \rho_{+ \nu } (   \partial^\mu\rho_-^\nu                    - \partial^\nu\rho_-^\mu )               \right ]    + ie \partial_\mu a_\nu \left [         \rho_-^\mu   \rho_+^\nu - \rho_-^\nu   \rho_+^\mu                           \right]\ , , \nonumber\\ \label{ggi } { \cal l}_{\pi\pi\sigma\gamma } & = & - 2 e g_\sigma a_\mu j_\pi^\mu \ , ,   \\
\label{ggj } { \cal l}_{\pi\pi\rho\gamma } & = & - e g_\rho a_\mu \left [      \bbox{\pi { \times}}(\bbox{\pi { \times } \rho}^\mu )                      \right]_3 \\ & = & e g_\rho a_\mu     \left [ 2 \pi^+ \pi^- \rho_0^\mu            - \pi^- \pi^0 \rho_+^\mu            - \pi^+ \pi^0 \rho_-^\mu     \right]\ , , \nonumber\\ \label{ggk } { \cal l}_{\pi\pi f\gamma } & = & - e g_f f_{\mu\nu } \left [ a^\mu j_\pi^\nu + a^\nu j_\pi^\mu                       \right]\,.\end{aligned}\ ] ] the corresponding vertex factors are @xmath253      & \mbox{\hspace*{5mm}for }         \rho_\pm^\alpha(p ) \to \rho_\pm^\beta(p^\prime ) a^\mu(q ) \\ % \label{vtppsg } \gamma_{\pi\pi\sigma\gamma } & = &     \pm 2ie g_\sigma ( p+p^\prime)^\mu & \mbox{\hspace*{5mm}for }                               \pi^{\pm}(p ) \to                              \pi^{\pm}(p^\prime )                              \sigma a^\mu \\ % \label{vtpprg } \gamma_{\pi\pi\rho\gamma } & = & \left\ {    \begin{array}{l }      -2ie g_\rho g^{\mu\nu } \\      + ie g_\rho g^{\mu\nu } \\      + ie g_\rho g^{\mu\nu } \\    \end{array } \right . &    \begin{array}{l }      \mbox{\hspace*{4mm}for }       \pi^{\pm}(p ) \to \pi^{\pm}(p^\prime ) \rho_0^\nu a^\mu
\\ %      \mbox{\hspace*{4mm}for }       \pi^{\pm}(p ) \to \pi^0(p^\prime ) \rho_{\pm}^\nu a^\mu \\
%      \mbox{\hspace*{4mm}for }       \pi^{0}(p ) \rho_{\pm}^\nu \to \pi^{\pm}(p^\prime ) a^\mu \\
\end{array } \\
% \label{vtppfg } \gamma_{\pi\pi f\gamma } & = &     \pm ie g_f \left [ g^{\alpha\mu}(p+p^\prime)^\beta                     + g^{\beta \mu}(p+p^\prime)^\alpha \right ]                            & \mbox{\hspace*{5mm}for }                               \pi^{\pm}(p ) \to                              \pi^{\pm}(p^\prime )                              f^{\alpha\beta } a^\mu    \ , .
\end{array}\end{aligned}\ ] ]                    g.  roche , talk presented at the nato advanced study institute , bodrum / turkey , september 1993 , and references therein ; in : proceedings of the nato asi on hot and dense nuclear matter , edited by w.  greiner , h.  stcker and a.  gallmann , nato asi series b335 .
serot and j.d .
walecka , _ the relativistic nuclear many - body problem _
, advances in nuclear physics , vol .
16 , edited by j.  w.  negele and e.  vogt , plenum , new york ( 1986 ) ; b.d .
serot and j.d .
walecka , rep .
phys .  * 55 * , 1855 ( 1992 ) .
figure 1 : breakdown of the soft photon approximation : the spa is valid only when @xmath23 is much smaller than the two solid lines shown .
the region between the dashed lines is the domain of integration for @xmath23 when calculating cross sections as a function of dilepton invariant mass .
figure 2 : one of the four contributing diagrams following from spa factorization of the matrix element into leptonic tensor @xmath83 , photon emission current @xmath255 , and elastic hadronic matrix element @xmath30 .
figure 5 : cross section for the reaction @xmath257 as a function of @xmath237 for fixed dielectron invariant masses @xmath27 and 300 mev .
solid line : exact obe calculation eq .
( [ xxc ] ) .
dash - dotted line : rckl approximation ( [ ick ] )  this is the commonly - used version of the spa .
other approximations are : dashed line : angular - averaged virtual @xmath242 current ( [ icklo ] ) .
dotted line : 3-phase space current ( [ bvp ] ) .
figure 10 : total bremsstrahlung yield of @xmath266 for all seven pion - pion reactions , for a boltzmann gas with temperature @xmath240 mev .
solid line : exact obe calculation , dash - dotted line : rckl approximation , dashed line : angular - averaged virtual-@xmath39 current , dotted line : 3-phase space current , virtual-@xmath39 current ( upper ) real-@xmath242 current ( lower ) .
figure 12 : ratios of spa approximation calculations divided by exact obe rate for @xmath240 mev .
dash - dotted line : rckl / exact , dashed line : ( angular - averaged virtual @xmath242 current)/exact .
upper dotted line : ( 3-phase space virtual @xmath242 current)/exact .
lower dotted line : ( 3-phase space real @xmath242 current)/exact .
figure 14 : fractional contribution to the total dielectron rate at @xmath243 mev as a function of dilepton invariant mass @xmath1 for ( a ) the rckl , ( b ) virtual - photon lichard approximations and ( c ) the exact calculation .
solid line ( upper ) : @xmath226 , dashed line : @xmath231 plus c.c .
, solid line ( lower ) : @xmath230 , dash - dotted line : @xmath229 , dotted line : @xmath228 plus c.c .
figure 15 : contribution of external - emission vs.  external - plus - internal diagrams for the reactions @xmath229 ( left ) and @xmath228 ( right ) , both for @xmath239 mev .
lines are as in figs .  7 and 9 . the new dash - dotted line below the ( solid line ) exact calculation represents contributions arising solely from emission of @xmath39 by an external pion line , but taking @xmath31 into account in @xmath245 , in contrast to the approximations ( upper lines ) which neglected @xmath31 .


manipulating quantum systems allows one to accomplish tasks far beyond the reach of classical devices .
quantum information is paradigmatic in this sense : quantum computers will substantially outperform classical machines for several problems  @xcite . though significant progress has been made recently , severe difficulties still remain ,
amongst which decoherence is certainly the most important .
large systems consisting of many qubits must be prepared in fragile quantum states , which are rapidly destroyed by their unavoidable coupling to the environment .
measurement - based feedback and coherent feedback are possible routes towards the preparation , protection and stabilization of such states . for coherent feedback strategy ,
the controller is also a quantum system coupled to the original one ( see  @xcite and the references therein ) .
this paper is devoted to measurement - based feedback where the controller and the control input are classical objects  @xcite .
the results presented here are directly inspired by a recent experiment  @xcite demonstrating that such a quantum feedback scheme achieves the on - demand preparation and stabilization of non - classical states of a microwave field .
following  @xcite and relying on continuous - time lyapunov techniques exploited in  @xcite , an initial measurement - based feedback was proposed in  @xcite .
this feedback scheme stabilizes photon - number states ( fock states ) of a microwave field ( see e.g. @xcite for a physical description of such cavity quantum electrodynamics ( cqed ) systems ) .
the controller consists of a quantum filter that estimates the state of the field from discrete - time measurements performed by probe atoms , and secondly a stabilizing state - feedback that relies on lyapunov techniques .
the discrete - time behavior is crucial for a possible real - time implementation of such controllers .
closed - loop simulations reported in  @xcite have been confirmed by the stability analysis performed in  @xcite . in the experimental implementation
@xcite , the state - feedback has been improved by considering a strict lyapunov function : this ensures better convergence properties by avoiding the passage by high photon numbers during the transient regime .
the goal of this paper is to present , for a class of discrete - time quantum systems , the mathematical methods underlying such improved lyapunov design .
our main result is given in theorem  [ thm : obscontrollermeas ] where closed - loop convergence is proved in presence of delays and measurement imperfections .
the state - feedback scheme may be applied to generic discrete - time finite - dimensional quantum systems using _ controlled _ measurements in order to deterministically prepare and stabilize the system at some pre - specified target state .
the dynamics of these systems may be expressed in terms of a ( classical ) nonlinear controlled markov chain , whose state space consists of the set of density matrices on some hilbert space .
the jumps in this markov chain are induced by quantum measurements and the associated jump probabilities are state - dependent .
these systems are subject to a discrete - time sequence of positive operator valued measurements ( povms @xcite ) and we use these povms to stabilize the system at the target state . by controlled measurements ,
we mean that at each time - step the chosen povm is not fixed but is a function of some classical control signal @xmath0 , similar to  @xcite .
however , we assume that when the control @xmath0 is zero , the chosen povm performs a quantum non - demolition ( qnd ) measurement  @xcite for some orthonormal basis that includes the target state . the feedback - law is based on a lyapunov function that is a linear combination of a set of martingales corresponding to the open - loop qnd measurements .
this lyapunov function determines a `` distance '' between the target state and the current state .
the parameters of this lyapunov function are given by inverting metzler matrices characterizing the impact of the control input on the kraus operators defining the markov processes and povms .
the ( graph theoretic ) properties of the metzler matrices are used to construct families of open - loop supermartingales that become strict supermartingales in closed - loop .
this fact provides directly the convergence to the target state without using the invariance principle .
a common problem that occurs in quantum feedback control is that of delays between the measurement process and the control process  @xcite . in this paper we demonstrate , using a predictive quantum filter , that the proposed scheme works even in the presence of delays .
convergence analysis is done for perfect and imperfect measurements . for imperfect measurements ,
the dynamics of the system are governed by a nonlinear markov chain given in  @xcite .    in both the perfect and imperfect measurement situations we prove a robustness property of the feedback algorithm : the convergence of the closed - loop system
is ensured even when the feedback law is based on the state of a quantum filter that is not initialized correctly .
this robustness property , is similar in spirit to the separation principle proven in  @xcite .
we use the fact that the state space is a convex set and the target state , being a pure state , is an extreme point of this convex state space and therefore can not be expressed as a convex combination of any other states .
one then uses the linearity of the conditional expectation to prove the robustness property .
our result is only valid for target quantum states that are pure states .
the paper is organized as follows . in section
[ sec : openloop ] , we describe the finite dimensional markov model together with the main modeling assumptions in the case of perfect measurements and study the open - loop behavior ( theorem  [ thm : first ] ) which can be seen as a non - deterministic protocol for preparing a finite number of isolated and orthogonal quantum states . in section  [ sec : closedloop ] , we present the main ideas underlying the construction of these control - lyapunov functions @xmath1 : a metzler matrix attached to the second derivative of the measurement operators and a technical lemma assuming this metzler matrix is irreducible .
finally , theorem  [ thm : main ] describes the stabilizing state feedback derived from @xmath1 .
the same analysis is done for the case of imperfect measurements in section  [ sec : imperfect ] . for the state estimations used in the feedback scheme we propose a brief discussion on the quantum filters and prove a rather general robustness property for perfect measurements in section  [ sec : closedloop ] with theorem  [ thm : obscontroller ] and for imperfect ones in section  [ sec : imperfect ] with theorem  [ thm : obscontrollermeas ] .
section  [ sec : photonbox ] is devoted to the experimental implementation that has been done at laboratoire kastler - brossel of ecole normale suprieure de paris .
closed - loop simulations and experimental data complementary to those reported in  @xcite are presented .
we consider a finite dimensional quantum system ( the underlying hilbert space @xmath2 is of dimension @xmath3 ) being measured through a generalized measurement procedure at discrete - time intervals .
the dynamics are described by a nonlinear controlled markov chain . here
, we suppose perfect measurements and no decoherence .
the system state is described by a density operator @xmath4 belonging to @xmath5 , the set of non - negative hermitian matrices of trace one : @xmath6 to each measurement outcome @xmath7 , @xmath8 being the numbers of possible outcomes , is attached the kraus operator @xmath9 depending on @xmath10 and on a scalar control input @xmath11 . for each @xmath0 ,
@xmath12 satisfy the constraint @xmath13 , the identity matrix .
the kraus map @xmath14 is defined by @xmath15 the random evolution of the state @xmath16 at time - step @xmath17 is modeled through the following dynamics : @xmath18 where @xmath19 is the control at step @xmath17 , subject to a delay of @xmath20 steps .
this delay is usually due to delays in the measurement process that can also be seen as delays in the control process .
@xmath21 is a random variable taking values @xmath10 in @xmath22 with probability @xmath23 . for each @xmath10 , @xmath24 is defined when @xmath25 .
we now state some assumptions that we will be using in the remainder of this paper .
[ assum : imp ] for @xmath26 , all @xmath27 are diagonal in the same orthonormal basis @xmath28 : @xmath29 with @xmath30    [ assum : impt ] for all @xmath31 in @xmath32 there  exists @xmath33 such that @xmath34    [ assum : imo ] all @xmath35 are @xmath36 functions of @xmath0 .    assumption  [ assum : imp ] means that when @xmath26 the measurements are quantum non demolition ( qnd ) measurements over the states @xmath37 : when @xmath38 , any @xmath39 ( orthogonal projector on the basis vector @xmath40 ) is a fixed point of  . since @xmath41 , we have @xmath42 for all @xmath43 according to assumption  [ assum : imp ]
. assumption  [ assum : impt ] means that there exists a @xmath10 such that the statistics when @xmath44 for obtaining the measurement result @xmath10 are different for the fixed points @xmath45 and @xmath46 .
this follows by noting that @xmath47 for @xmath43 .
assumption  [ assum : imo ] is a technical assumption we will use in our proofs .      when the control input vanishes ( @xmath48 ) , the dynamics are simply given by @xmath49 where @xmath21 is a random variable with values in @xmath22 .
the probability @xmath50 to have @xmath51 depends on @xmath52 : @xmath53 .
[ thm : first ] consider a markov process @xmath52 obeying the dynamics of   with an initial condition @xmath54 in @xmath5 .
then with probability one , @xmath52 converges to one of the @xmath55 states @xmath56 with @xmath43 and the probability of convergence towards the state @xmath56 is given by @xmath57    this theorem is already proved in  @xcite and also in  @xcite in a slightly different formulation . a direct proof is based on the following lyapunov function : @xmath58 since @xmath59 is a martingale for any @xmath60 and since @xmath61
in subsection  [ ssec : overview ] , we give an overview of the control method and then in subsections  [ ssec : weights ] and  [ ssec : stab ] , we prove the main results . finally in subsection
[ ssec : qfilter ] , we prove a robustness principle that explains how we can ensure convergence even if the initial state @xmath54 is unknown .
theorem  [ thm : first ] shows that the open - loop dynamics are stable in the sense that in each realization @xmath52 converges ( non - deterministically ) to one of the pure states @xmath56 with probability @xmath62 .
the control goal is to make this convergence deterministic toward a chosen @xmath63 playing the role of controller set point .
we build on the ideas in  @xcite to design a controller that is based on a strict lyapunov function for the target state . in this paper
we assume arbitrary controlled kraus operators @xmath64 that can not be decomposed into qnd measurement operators @xmath65 followed by a unique controlled unitary operator @xmath66 with @xmath67 as assumed in  @xcite where @xmath68 . it can be argued that the control we are proposing is non - hamiltonian control  @xcite , as the control parameter @xmath0 is not necessarily a parameter in the interaction hamiltonian and
could indeed be any parameter of an auxiliary system such as the measurement device .    to convey the main ideas involved in the control design , we begin with the case where there are no delays ( @xmath69 ) and we assume the initial state is known .
we wish to use the open - loop supermartingales to design a lyapunov function for the closed - loop system . by an open - loop supermartingale we mean any function @xmath70 such that @xmath71 for all @xmath72 .
the lyapunov function underlying theorem  [ thm : first ] demonstrates how we can construct such open - loop supermartingales .
now , at each time - step @xmath17 , the feedback signal @xmath73 is chosen by minimizing this supermartingale @xmath74 knowing the state @xmath52 : @xmath75}{\text{argmin } } \bigg\ { { { \text{\normalsize$\mathbb e$}}\left(v(\rho_{k+1 } ) |\rho_k = \rho , u_k = u\right)}\bigg\}.\ ] ] here @xmath76 is some small positive number that needs to be determined .
because @xmath77 $ ] and the control @xmath73 is chosen to minimize @xmath74 at each step , we directly have that @xmath74 is a _ closed - loop supermartingale _ , i.e. , @xmath78 for all @xmath72 .
if this supermartingale @xmath74 is bounded from below then we can directly apply the convergence theorem  [ thm : app ] in the appendix to prove that @xmath52 converges with probability one to the set @xmath79    what remains to be done is to choose an appropriate open - loop supermartingale @xmath74 so that the set @xmath80 is restricted to the target state @xmath81 .
the lyapunov function @xmath82 defined in   does not discriminate between the different basis vectors .
we therefore use the lyapunov function @xmath83 where @xmath84 and @xmath85 is a small positive constant .
the weights @xmath86 are strictly positive numbers except for @xmath87 .
this function @xmath88 is clearly a concave function of the open - loop martingales @xmath59 and therefore is an open - loop supermartingale .
moreover , the weights @xmath86 can be used to quantify the distance of the state @xmath4 from the target state @xmath89 ( c.f .
figure  [ fig : distance ] below which shows how @xmath86 are chosen for the experimental setting ) .
a state @xmath4 is in the set @xmath80 if and only if for all @xmath90 $ ] , we have @xmath91 also from the fact that @xmath88 is an open - loop supermartingale , we have for all @xmath72 @xmath92 we prove in lemma  [ lem : first ] below that given any @xmath93 , we can always choose the weights @xmath94 so that @xmath88 satisfies the following property : @xmath95 , @xmath96 has a strict local minimum at @xmath26 for @xmath97 and strict local maxima at @xmath98 for @xmath99 .
this combined with equation   then ensures that for any @xmath99 , there is some @xmath90 $ ] such that @xmath100 therefore using equation  , we know that @xmath56 is in the limit set @xmath80 if and only if @xmath97 .
this idea can easily be extended to the situations where the delay @xmath101 is non zero .
take @xmath102 as state at step @xmath17 : denote by @xmath103 this state where @xmath104 stands for the control input @xmath0 delayed @xmath105 steps .
then the state form of the delayed dynamics   is governed by the following markov chain @xmath106 the goal is to design a feedback law @xmath107 that globally stabilizes this markov chain @xmath108 towards a chosen target state @xmath109 for some @xmath110 . in theorem
[ thm : main ] , we show how to design a feedback relying on the control lyapunov function @xmath111 .
the idea is to use a predictive filter to estimate the state of the system @xmath101 time - steps later .
finally , we address the situation where the initial state of the system is not fully known but only estimated by @xmath112 .
we show under some assumptions on the initial condition that , the feedback law based on the state @xmath113 of the miss - initialized filter still ensures the convergence of @xmath113 as well as the well - initialized conditional state @xmath52 towards @xmath114 .
this demonstrates how the control algorithm is robust to uncertainties in the initialization of the estimated state of the quantum system .
the construction of the control lyapunov function relies on two lemmas .
[ lem : r ] consider the @xmath115 matrix @xmath116 defined by @xmath117 when @xmath118 , the non - negative @xmath119 is a right stochastic matrix .
for @xmath31 , @xmath120 .
thus @xmath116 is a metzler matrix  .
let us prove that the sum of each row vanishes .
this results from identity @xmath121 . deriving twice
versus @xmath0 the relation @xmath122 yields @xmath123 since for @xmath26 , @xmath124 the above sum corresponds to @xmath125 .
therefore , the diagonal elements of @xmath116 are non - positive . if @xmath118 , then @xmath126 and the matrix @xmath119 is well defined with non - negative entries .
since the sum of each row of @xmath116 vanished , the sum of each row of @xmath127 is equal to @xmath128 .
thus @xmath127 is a right stochastic matrix .    to the metzler matrix @xmath116 defined in lemma  [ lem : r ] , we associate its directed graph denoted by @xmath129 this graph admits @xmath130 vertices labeled by @xmath131 . to each strictly positive off - diagonal element of the matrix @xmath116 , say , on the @xmath132th row and the @xmath133th column we associate an edge from vertex @xmath132 towards vertex @xmath134    [ lem : first ] assume the directed graph @xmath135 of the matrix @xmath116 defined in lemma  [ lem : r ] is strongly connected , i.e. , for any @xmath136 , @xmath137 , there exists a chain of @xmath105 distinct elements @xmath138 of @xmath139 such that @xmath140 , @xmath141 and for any @xmath142 , @xmath143 .
take @xmath144 .
then , there exist @xmath145 strictly positive real numbers @xmath146 , @xmath147 , such that    * for any reals @xmath148 , @xmath147 , there exists a unique vector @xmath149 of @xmath150 with @xmath151 such that @xmath152 where @xmath153 is the vector of @xmath150 of components @xmath148 for @xmath147 and @xmath154 ; if additionally @xmath155 for all @xmath147 , then @xmath156 for all @xmath147 . * for any vector @xmath157 , solution of @xmath158 , the function @xmath159 satisfies @xmath160    since the directed graph @xmath135 coincides with the directed graph of the right stochastic matrix @xmath127 defined in lemma  [ lem : r ] , @xmath127 is irreducible .
since it is a right stochastic matrix , its spectral radius is equal to @xmath128 . by perron - frobenius theorem for non - negative irreducible matrices , this spectral radius ,
i.e. , @xmath128 , is also an eigen - value of @xmath127 and of @xmath161 , with multiplicity one and associated to eigen - vectors having strictly positive entries : the right eigen - vector ( @xmath162 ) is obviously @xmath163 ; the left eigen - vector @xmath164 ( @xmath165 ) can be chosen such that @xmath166 .
consequently , the rank of @xmath116 is @xmath145 with @xmath167 and @xmath168 where @xmath169 is the hyper - plane orthogonal to @xmath170 . since @xmath171 , @xmath172 , exists @xmath173 such that @xmath174 . since @xmath167
, there is a unique @xmath173 solution of @xmath175 such that @xmath87 .
the fact that @xmath176 when @xmath177 for @xmath178 , comes from elementary manipulations of @xmath179 showing that @xmath180 .
@xmath181 , set @xmath182 set @xmath183 then @xmath184 is given by @xmath185 and @xmath186 therefore @xmath187 and @xmath188      the main result of this section is expressed through the following theorem .
[ thm : main ] consider the markov chain   with assumptions  [ assum : imp ] ,  [ assum : impt ] and  [ assum : imo ] .
take @xmath63 and assume that the directed graph @xmath135 associated to the metzler matrix @xmath116 of lemma  [ lem : r ] is strongly connected .
take @xmath189 , @xmath190 the solution of @xmath152 with @xmath87 , @xmath177 for @xmath191 , @xmath192 ( see lemma  [ lem : first ] ) and @xmath193 take @xmath194 and consider the following feedback law @xmath195}{\text{argmin}}\big({{\text{\normalsize$\mathbb e$}}\left(w_\epsilon(\chi_{k+1})|\chi_k , u_k=\xi\right)}\big ) : = f(\chi_k)\ ] ] where @xmath196 then there exists @xmath197 such that , for all @xmath1980,u^*]$ ] and @xmath1990,\min_{n\neq \bar n }
\left(\frac{\lambda_{n}}{r_{n , n}}\right)\right]$ ] , the closed - loop markov chain of state @xmath200 with the feedback law of theorem  [ thm : main ] converges almost surely towards @xmath201 for any initial condition @xmath202^\tau$ ] .    for the sake of simplicity ,
first we demonstrate this theorem for @xmath203 and thus for @xmath204 .
we then explain how this proof may be extended to arbitrary @xmath205 . here
, @xmath206 is given by @xmath207 that can also be presented as @xmath208 with @xmath209 : @xmath210}{\text{min}}\sum_{\mu } p_{\mu,\rho_k}^{\beta_{1,k}}v_\epsilon({{\mathbb k}}^{\xi}({{\mathbb m}}_\mu^{\beta_{1,k}}(\rho_k ) ) )      .\end{gathered}\ ] ] since @xmath211 for any @xmath212 , we have @xmath213 where @xmath214 is defined in  .
thus , @xmath215 is a supermartingale . according to theorem  [ thm : app ] ,
the @xmath216-limit set @xmath80 of @xmath108 is included in @xmath217 ~|~ q_1(\rho,\beta_1)=0,~q_2(\rho,\beta_1)=0\right\}.\ ] ] we show in 3 steps that @xmath80 is restricted to the target state . + * step 1 .
* for all @xmath218 , there exists @xmath219 such that @xmath220 ~|~ { \left < n|\rho|n\right > } \geq 1-\delta\bigg\}.$ ] + * proof of step 1 . *
step 1 follows from @xmath221 for all @xmath222 and technical lemma  [ lem : q1 ] of the appendix . +
* there exist @xmath223 and @xmath224 such that @xmath225 ~|~ { \left<\bar{n}|\rho|\bar{n}\right > } \geq 1-\delta'\}.$ ] + * proof of step 2 .
* if @xmath226 , the minimum of the function @xmath227\ni \xi \mapsto f_{\rho,\beta_1}(\xi)=\sum_{\mu } p_{\mu,\rho}^{\beta_1}v_\epsilon({{\mathbb k}}^{\xi}({{\mathbb m}}_\mu^{\beta_1}(\rho ) ) ) $ ] is attained at @xmath228 . by construction of @xmath88 , we know that , for @xmath229 , @xmath39 , @xmath230 if @xmath178 , and @xmath231 the dependence of @xmath232 versus @xmath4 and @xmath233 is continuous from assumption  [ assum : imo ] .
therefore , @xmath234 such that for all @xmath235 in @xmath236 ~|~ { \left < n|\rho|n\right > } \geq 1-\delta'\}$ ] we have @xmath237 therefore , such @xmath235 and @xmath228 , can not minimize @xmath238 and we have @xmath239 ~|~ { \left < n|\rho|n\right > } \geq 1-\delta'\ } = \emptyset .
$ ] therefore we get , @xmath240 ~|~ { \left<\bar{n}|\rho|\bar{n}\right > } \geq 1-\delta'\},\ ] ] by setting @xmath241 in step 1 , for @xmath76 small enough . + * step 3 .
* @xmath76 can be chosen small enough so that @xmath242 . + * proof of step 3 . * by construction of @xmath88 , we have for @xmath243 and @xmath244 , @xmath245 by continuity of @xmath232 with respect to @xmath4 and @xmath233 , we can choose @xmath246 and @xmath76 small enough such that for all @xmath247 $ ] and @xmath4 satisfying @xmath248 we have @xmath249 this implies that on the domain @xmath250 ^ 2 ~|~ { \left < n|\rho|n\right > } \geq 1-\bar{\delta}\},$ ] @xmath251 is a uniformly strongly convex function of @xmath252 $ ] .
thus the argument of its minimum over @xmath252 $ ] is a continuous function of @xmath4 and @xmath233 .
now we choose @xmath253 , where @xmath254 is as in step 2 .
take a convergent subsequence of @xmath255 ( that we still denote by @xmath108 for simplicity sakes ) .
its limit @xmath256 belongs to @xmath80 and also to @xmath257 ~|~ { \left < n|\rho|n\right > } \geq 1-\bar{\delta}\}.$ ] therefore , @xmath258}{\text{argmin } } f_{\rho_k,\beta_{1,k}}(\xi)$ ] converges to @xmath259}{\text{argmin } } f_{\rho_\infty,\beta_{1,\infty}}(\xi)=0 $ ] since @xmath260 . because @xmath261 , @xmath262 tends almost surely towards @xmath263 .
we know that @xmath264 and @xmath265 tend also almost surely to @xmath263 ( theorem  [ thm : app ] ) and by uniform ( in @xmath4 ) continuity with respect to @xmath262 , @xmath266 and @xmath267 tend almost surely to @xmath263 .
since @xmath268 implies @xmath269 , @xmath52 converges almost surely towards @xmath89 .
this completes the proof for @xmath270 . + * extension to @xmath205 .
* for @xmath271 and @xmath103 , the proof is very similar .
we still have   with @xmath214 and @xmath272 given by formulae analogous to   and  : @xmath262 is replaced by @xmath273 ; @xmath274 is replaced by @xmath275 ; @xmath276 is replaced by @xmath277 since kraus maps are trace preserving .
the analogue of technical lemma  [ lem : q1 ] , whose proof relies on similar continuity and compactness arguments , has now the following statement : _
@xmath278 such that @xmath279^\tau $ ] satisfying @xmath280 , @xmath281 such that @xmath282 .
_   + the last part of the proof showing that , for @xmath283 small enough , the control @xmath73 is a continuous function of @xmath108 when @xmath108 is in the neighborhood of the @xmath216-limit set @xmath284^\tau~|~q_1(\chi)=q_2(\chi)=0\}$ ] remains almost the same .
when the measurement process is fully efficient ( i.e. , the detectors are ideal with detection efficiency equal to one and they observe all the quantum jumps ) and the jump model   admits no error , the markov process   represents a natural choice for estimating the hidden state @xmath285 indeed , the estimate @xmath286 of @xmath4 satisfies the following recursive dynamics @xmath287 where the measurement outcome @xmath21 is driven by  . in practice ,
the control @xmath73 defined in theorem  [ thm : main ] could only depend on this estimation @xmath113 replacing @xmath52 in @xmath108 . if @xmath288 , then @xmath289 and theorem  [ thm : main ] ensures convergence towards the target state . otherwise , the following result guaranties the convergence of such observer / controller scheme when @xmath290 .
[ thm : obscontroller ] consider the recursive equation   and assume that the assumptions of theorem  [ thm : main ] are satisfied .
for each measurement outcome @xmath21 given by  , consider the estimation @xmath291 given by   with an initial condition @xmath112 .
set @xmath292 where @xmath293 is given in theorem  [ thm : main ]
. then there exists @xmath197 such that , for all @xmath1980,u^*]$ ] and @xmath1990,\min_{n\neq \bar n } \left(\frac{\lambda_{n}}{r_{n , n}}\right)\right]$ ] , @xmath52 and @xmath113 converge almost surely towards the target state @xmath89 as soon as @xmath290 .
set @xmath294 .
the main idea is that @xmath295 ( where we take the expectation over all jump realizations ) depends linearly on @xmath54 even though we apply the feedback control .
since @xmath19 is a function of @xmath296 , @xmath297 is given by @xmath298 , where @xmath299 .
the linearity of @xmath295 with respect to @xmath54 is thus verified .
now we apply the assumption @xmath290 and therefore one can find a constant @xmath3000,1[$ ] and a well - defined density matrix @xmath301 in @xmath302 such that @xmath303 applying the dominated convergence theorem : @xmath304 by the linearity of @xmath305 with respect to @xmath54 , @xmath306 is given by @xmath307 as @xmath308 and @xmath309 , we necessarily have that both of them converge to @xmath128 since @xmath310 : @xmath311 this implies the almost sure convergence of @xmath52 towards the pure state @xmath312 .
we have also @xmath313 depending linearly on @xmath54 .
thus @xmath314 is given by @xmath315 using the fact that when @xmath316 , @xmath317 for all @xmath17 , we conclude similarly that @xmath113 converges almost surely towards @xmath312 even if @xmath54 and @xmath112 do not coincide .
we now consider the feedback control problem in the presence of classical measurement imperfections with the possibility of detection errors .
this model is a direct generalization of the ones used in  @xcite ( see also e.g. ,  @xcite for an introduction to quantum filtering ) .
the imperfections in the measurement process are described by a classical probabilistic model relying on a left stochastic matrix @xmath318 , @xmath319 and @xmath33 : @xmath320 and for any @xmath10 , @xmath321 .
the integer @xmath322 corresponds to the number of imperfect outcomes and @xmath323 is the probability of having the imperfect outcome @xmath324 knowing the perfect one @xmath10 .
set @xmath325 since @xmath52 follows  , @xmath326 is also governed by a recursive equation  @xcite : @xmath327 where for each @xmath324 , @xmath328 is the super - operator defined by @xmath329 with @xmath330 , and where @xmath331 is a random variable taking values @xmath324 in @xmath332 with probability @xmath333 since @xmath323 is a left stochastic matrix , @xmath334 which is precisely the kraus map   associated with the markov process @xmath52 . by assumption  [ assum : imp ] , each pure state @xmath56 , @xmath131 remains a fixed point of the markov process   when @xmath48 .
we now consider the dynamics of the filter state @xmath335 in the presence of delays in the feedback control .
similar to the case with perfect measurements , let @xmath336 be the filter state at step @xmath17 , where @xmath337 is the feedback control at time - step @xmath17 delayed @xmath338 steps .
then the delay dynamics are determined by the following markov chain @xmath339 instead of assumption  [ assum : impt ] we now assume    [ assum : hatimpt ] for all @xmath31 in @xmath32 there exists @xmath319 such that @xmath340    assumption  [ assum : hatimpt ] means that there exists a @xmath324 such that the statistics when @xmath48 for obtaining the measurement result @xmath324 are different for the fixed points @xmath45 and @xmath46 .
this follows by noting that @xmath341 for @xmath43 .
we now state the analogue of theorem  [ thm : main ] in the case of imperfect measurements .
[ thm : mainmeas ] consider the markov chain   with assumptions  [ assum : imp ] , [ assum : imo ] and  [ assum : hatimpt ] .
take @xmath63 .
assume that the directed graph of metzler matrix @xmath116 of lemma  [ lem : r ] is strongly connected .
take @xmath189 , @xmath190 solution of @xmath152 with @xmath87 , @xmath177 for @xmath191 , @xmath192 ( see lemma  [ lem : first ] ) and set @xmath342 take @xmath194 and consider the following feedback law @xmath343}{\text{argmin}}\big({{\text{\normalsize$\mathbb e$}}\left(w_\epsilon({{\widehat{\chi}}}_{k+1})|{{\widehat{\chi}}}_k , u_k=\xi\right)}\big ) = f({{\widehat{\chi}}}_k ) , $ ] where @xmath344 + then there exists @xmath197 such that , for all @xmath1980,u^*]$ ] and @xmath1990,\min_{n\neq \bar n } \left(\frac{\lambda_{n}}{r_{n , n}}\right)\right]$ ] , the closed - loop markov chain of state @xmath345 with the feedback law of theorem  [ thm : mainmeas ] converges almost surely towards @xmath201 for any initial condition @xmath346^\tau$ ] .
the proof of this theorem is almost identical to that of theorem  [ thm : main ] with @xmath347 replaced by @xmath348 and @xmath349 replaced by @xmath350 and we do not give the details of this proof .
+ for estimating the hidden state @xmath335 needed for the feedback design of theorem  [ thm : mainmeas ] , let us consider the estimate @xmath351 given by @xmath352 where @xmath331 corresponds to the imperfect outcome detected at step @xmath17 .
such @xmath331 is correlated to the perfect and hidden outcome @xmath21 of   through the classical stochastic process attached to  @xmath318 : for each @xmath21 , @xmath331 is a random variable to be equal to @xmath319 with probability @xmath353 . in practice ,
the control @xmath73 defined in theorem  [ thm : mainmeas ] could only depend on this estimation @xmath354 replacing @xmath326 in @xmath355 .
the following result guaranties the convergence of the feedback scheme when @xmath356 .
[ thm : obscontrollermeas ] consider the recursive equation   and take assumptions of theorem  [ thm : mainmeas ] .
consider the estimation @xmath357 given by   with an initial condition @xmath358 .
set @xmath359 where @xmath293 is given by theorem  [ thm : mainmeas ]
. then there exists @xmath197 such that , for all @xmath1980,u^*]$ ] and @xmath1990,\min_{n\neq \bar n } \left(\frac{\lambda_{n}}{r_{n , n}}\right)\right]$ ] , @xmath52 and @xmath354 converge almost surely towards the target state @xmath89 as soon as @xmath356 .
let us first prove that @xmath326 defined by   converges almost surely towards @xmath89 . since @xmath360 , we have @xmath361 .
thus , there exist @xmath362 and @xmath3000,1[$ ] , such that , @xmath363 .
similarly to the proof of theorem  [ thm : obscontroller ] , @xmath364 depends linearly on @xmath365 : @xmath366 is given by @xmath367 . moreover by theorem  [ thm : mainmeas ]
, @xmath366 converges almost surely towards  @xmath128 , thus @xmath368 and @xmath369 converge also towards @xmath128 .
consequently @xmath326 converges almost surely towards @xmath89 .
since @xmath326 is the conditional expectation of @xmath52 knowing the past imperfect outcomes and control inputs and since its limit @xmath89 is a pure state , @xmath52 converges necessarily towards the same pure state almost surely .
convergence of @xmath354 relies on similar arguments exploiting the linearity of @xmath370 versus @xmath365 .
in this section , we give the explicit expression of the feedback controller which has been experimentally tested in laboratoire kastler - brossel ( lkb ) at ecole normal suprieure ( ens ) de paris .
we briefly summarize how the control design elucidated in this paper is applied to the lkb experiment .
we refer the interested reader to  @xcite for more details .
this feedback controller has been obtained by the lyapunov design discussed in previous sections .
[ fig : setup ]    figure  [ fig : setup ] is a sketch of the experimental setup considered in  @xcite .
two - level atoms , with state space spanned by the vectors @xmath371 and @xmath372 ( excited and ground state respectively ) , act as qubits .
the atoms , emitted from box b , interact individually with the electromagnetic microwave field stored in the cavity  @xmath373 and get entangled with it .
the cavity is placed in between two additional cavities , @xmath374 and @xmath375 , where the atomic state can be manipulated at will .
together , they form an atomic interferometer similar to those used in atomic clocks . the atomic state , @xmath371 or @xmath372 , is eventually measured in the detection device  @xmath376 . using the outcomes of those measurements , a quantum filter , run by a real - time computer , estimates the density matrix @xmath4 of the microwave field trapped in the cavity  @xmath373 .
it also uses all the available knowledge on the experimental setup , especially the limited detection efficiency @xmath377 ( percentage of atoms that are actually detected ) , the detection errors @xmath378 and @xmath379 ( @xmath380 is the probability to detect an atom in the wrong state ) .
a controller eventually calculates ( state feedback ) the amplitude @xmath381 of the classical microwave field that is injected into @xmath373 by an external source @xmath382 so as to bring the cavity field closer to the target state .
the succession of the atomic detection , the state estimation and the microwave injection define one iteration of our feedback loop .
it corresponds to a sampling time of around @xmath383s . for a more detailed description see  @xcite .
the state to be stabilized is the cavity state .
the underlying hilbert space is @xmath384 which is assumed to be finite - dimensional ( truncations to @xmath385 photons ) .
it admits @xmath386 as an ortho - normal basis .
each basis vector @xmath387 corresponds to a pure state , called photon - number state ( fock state ) , with precisely @xmath388 photons in the cavity , @xmath389 with notations of subsection  [ ssec : model ] , @xmath390 and the index @xmath391 .
the three main processes that drive the evolution of our system are decoherence , injection , and measurement .
their action onto the system s state can be described by three super - operators @xmath392 , @xmath393 and @xmath394 , respectively .
we refer to  @xcite for more details .
all operators are expressed in the fock - basis @xmath395 truncated to @xmath385 photons .    after taking into account our full knowledge about the experiment , we finally get the following state estimate at step @xmath396 : @xmath397    with the following short descriptions for the super - operators @xmath392 , @xmath393 and @xmath398    * the decoherence manifests itself through spontaneous loss or capture of a photon to or from the environment which is described as follows : @xmath399 where @xmath400 , @xmath401 and @xmath402 with @xmath403 @xmath404 and @xmath405 photon annihilation and creation operators ( @xmath406 and @xmath407 ) and with @xmath408 the photon number operator ( @xmath409 ) . besides , @xmath410 is the mean number of photons in the cavity mode at thermal equilibrium with its environment . *
the evolution of the state @xmath4 after the control injection is modeled through @xmath411 with @xmath412 in reality , the control at step @xmath17 , @xmath413 is subject to a delay of @xmath20 steps which corresponds to the number of flying atoms between the cavity @xmath373 and the detector @xmath376 .
* in the real experiment , the atom source is probabilistic and is characterized by a truncated poisson probability distribution @xmath414 to have @xmath415 atom(s ) in a sample ( we neglect events with more than 2 atoms ) .
this expands the set of the possible detection outcomes to @xmath416 values @xmath417 , related to the following measurement operators , @xmath418 @xmath419 , @xmath420 , + @xmath421 , @xmath422 and @xmath423 where @xmath424 @xmath425 and @xmath426 are physical parameters .
+ the real measurement process is not perfect : the detection efficiency is limited to @xmath427 and the state detection errors are non - zero ( @xmath428 .
these imperfections are taken into account by considering the left stochastic matrix @xmath323 which is given in  @xcite .
consequently , the optimal state estimate after measurement outcome @xmath429 gets the following form : @xmath430 the measurement operators @xmath431 are diagonal in the fock basis @xmath432 , illustrating their quantum non - demolition nature with respect to the photon number operator and thus fulfilling assumption  [ assum : imp ] .
besides , assumption  [ assum : hatimpt ] can also be fulfilled by a proper choice of the experimental parameters @xmath433 and @xmath426
.      for @xmath434 ( no cavity decoherence ) the markov model of density matrix @xmath335 associated to the filter   is exactly of the form   with @xmath435 being fixed - points in open - loop .
similarly , the underlying markov process of the true cavity state @xmath4 , which is unobservable in practice because of detection errors and delays , admits the same fixed points in open - loop . with parameters given in subsection  [ ssec : exp ] ( except @xmath436 ) , these markov processes satisfy assumptions  [ assum : imp]-[assum : impt]-[assum : imo ] for @xmath4 and assumptions  [ assum : imp]-[assum : imo]-[assum : hatimpt ] for @xmath335 .
moreover the metzler matrix @xmath116 of lemma  [ lem : first ] is irreducible .
consequently the assumptions of theorem  [ thm : mainmeas ] are satisfied .
the feedback law proposed in theorem  [ thm : obscontrollermeas ] and relying on the filter state @xmath351 will stabilize globally the unobservable state @xmath4 towards the target photon - number state @xmath437 .
numerous closed - loop simulations show that taking @xmath438 in the feedback law does not destroy stability and does not affect the convergence rates .
this explains why in the simulations and experiments , we set @xmath438 despite the fact that theorem  [ thm : obscontrollermeas ] guaranties convergence only for arbitrary small but strictly positive @xmath85 .
for @xmath439 positive and small , the @xmath440s are no more fixed - points for @xmath4 and @xmath335 in open - loop .
let us detail how to adapt the feedback scheme of theorem  [ thm : obscontrollermeas ] .
at each step of an ideal experiment the control @xmath73 minimizes the lyapunov function @xmath441 ( @xmath85 is set to zero ) calculated for state @xmath442 . in our real experiment
however , we also take into account decoherence and @xmath101 flying not - yet - detected samples and therefore choose @xmath73 to minimize @xmath443 for @xmath444 with @xmath351 given by  . here
we have introduced the kraus map of the real experiment @xmath445 .
the control @xmath0 minimizing @xmath446 is approximated by @xmath259}{\operatorname{argmax}}\left(a_1 \xi + a_2 \xi^2/2\right)$ ] with @xmath447\rho)}$ ] and @xmath448,{\text{\bf{a}}}^\dag\!-\!{\text{\bf{a}}}]\rho)}$ ] , where @xmath449 is the diagonal operator @xmath450 .
the coefficients @xmath451 are computed using lemma  [ lem : first ] where , for @xmath452 , @xmath453 are chosen negative and with a decreasing modulus versus @xmath388 in order to compensate cavity decay .
for @xmath454 , we have compared in simulations different setting and selected the profile displayed in figure  [ fig : distance ] .     used in the feedback for the simulation of figure  [ fig : simulation ] and for the experiment of figure  [ fig : experimental].,height=226 ]      closed - loop simulation of figure  [ fig : simulation ] shows a typical monte - carlo trajectory of the feedback loop aiming to stabilize the 3-photon state @xmath455 .
the experimental parameters used in the simulations are the following : @xmath456 , @xmath457 , @xmath458 , @xmath459 , @xmath460 , @xmath461 , @xmath462 , @xmath463 , @xmath464 , and @xmath465 .
for the feedback , @xmath451 are given in figure  [ fig : distance ] , @xmath438 and @xmath466 .
the initial states @xmath54 and @xmath358 take the following form : @xmath467 .    .
( a ) detection results .
( b ) control lyapunov function .
( c ) control input .
( d ) estimated photon number probabilities : @xmath468 in red , @xmath469 in thick green , and @xmath470 in blue .
( e ) cavity photon number probabilities : @xmath471 in red , @xmath472 in thick green , and @xmath473 in blue .
, scaledwidth=50.0% ]    the results of the experimental implementation of the feedback scheme are presented in figure  [ fig : experimental ] .
figure  [ fig : experimental](e ) shows that the average fidelity of the target state is about @xmath474 . besides , the asymmetry between the distributions for @xmath475 and @xmath476 indicates the presence of quantum jumps occurring preferentially downwards ( @xmath477 .
contrarily to the simulations of figure  [ fig : simulation ] , the cavity photon number probabilities relying on @xmath4 are not accessible in the experimental data of figure  [ fig : experimental ] since we do not have access to the detection errors and to the cavity decoherence jumps  @xcite .
nevertheless , green curves in simulations of figures  [ fig : simulation](d ) and  [ fig : simulation](e ) indicate that when @xmath469 exceeds @xmath478 , @xmath4 coincides , with high probability , with @xmath437 .    .
( a ) detection results .
( b ) evolution of the control function .
( c ) control injection .
( d ) estimated photon number probabilities : @xmath468 in red , @xmath469 in thick green , and @xmath470 in blue .
( e ) estimated photon number probabilities averaged over 4000 experimental trajectories of the feedback loop.,scaledwidth=50.0% ]
we have proposed a lyapunov design for state - feedback stabilization of a discrete - time finite - dimensional quantum system with qnd measurements .
extensions of this design are possible in different directions such as    * replacing the continuous and one - dimensional input @xmath0 by a multi - dimensional one @xmath479 ; * assuming that @xmath0 belongs to a finite set of discrete values ; * taking an infinite dimensional state space as in  @xcite where the truncation to finite photon numbers is removed ; * considering continuous - time systems similar to the ones investigated in  @xcite ; * ensuring convergence towards a sub - space  @xcite instead of a pure - state and thus achieving a goal similar to error correction code as already proposed in  @xcite .
* acknowledgements : * the authors thank michel brune , serge haroche and jean - michel raimond for enlightening discussions and advices .
10 c.  ahn , a.  c. doherty , and a.  j. landahl .
continuous quantum error correction via quantum feedback control .
, 65:042301 , mar 2002 .
h.  amini , m.  mirrahimi , and p.  rouchon .
stabilization of a delayed quantum system : the photon box case - study .
, 57(8):19181930 , 2012 .
h.  amini , p.  rouchon , and m.  mirrahimi . design of strict control - lyapunov functions for quantum systems with qnd measurements . in _
50th ieee conference on decision and control and european control conference _ , pages 81938198 , 2011 .
michel bauer and denis bernard . convergence of repeated quantum nondemolition measurements and wave - function collapse .
, 84:044103 , oct 2011 .
s. bolognani and f. ticozzi . engineering stable discrete - time quantum dynamics via a canonical qr decomposition .
, 55(12):27212734 , 2010 .
l.  bouten , r.  van handel , and m. r.  james .
an introduction to quantum filtering .
46(6):21992241 , 2007 .    l.  bouten and r.  van handel .
, chapter on the separation principle of quantum control .
world scientific , 2008 .
( see also : arxiv : math - ph/0511021 )    l.  bouten , r.  van handel , and m. r.  james . a discrete invitation to quantum filtering and feedback control .
, 51:239316 , 2009 .    m.  brune , j.  bernu , c.  guerlin , s.  deleglise , c.  sayrin , s.  gleyzes , s.  kuhr , i.  dotsenko , j - m .
raimond , s.  haroche .
process tomography of field damping and measurement of fock state lifetimes by quantum nondemolition photon counting in a cavity .
, 101 , 240402 , 2008 .
i.  dotsenko , m.  mirrahimi , m.  brune , s.  haroche , j .-
raimond , and p.  rouchon .
quantum feedback by discrete quantum non - demolition measurements : towards on - demand generation of photon - number states .
, 80 : 013805 - 013813 , 2009 .    c. w.  gardiner and p.  zoller .
quantum noise . , 2000 .
deterministic and nondestructively verifiable preparation of photon number states . , 97(073601 ) , 2006 .
j.  gough and m.r .
james . the series product and its application to quantum feedforward and feedback networks .
, 54(11):2530 2544 , 2009 .    c.  guerlin , j.  bernu , s.  delglise , c.  sayrin , s.  gleyzes , s.  kuhr , m.  brune , j - m .  raimond and s.  haroche
. progressive field - state collapse and quantum non - demolition photon counting .
, 448 : 889894 , 2007 .
s.  haroche and j.m .
oxford university press , 2006 .
james and j.e . gough .
quantum dissipative systems and feedback control design by interconnection .
, 55(8):1806 1821 , 2010 .
k. kashima and n. yamamoto .
control of quantum systems despite feedback delay .
, 54(4):876 881 , 2009 .
kushner . .
holt , rinehart and wilson , inc . , 1971 .
m.  mirrahimi and r.  van handel . stabilizing feedback controls for quantum systems .
, 46(2):445467 , 2007 .
nielsen and i.l . chuang . .
, 2000 .
k.  nishio , k.  kashima , and j.  imura .
effects of time delay in feedback control of linear quantum systems .
, 79:062105 , jun 2009 .
r.  romano , d.  dalessandro .
incoherent control and entanglement for two - dimensional coupled systems .
, 73:022323 , 2006 .
r.  romano , d.  dalessandro .
environment - mediated control of a quantum system .
, 97:080402 , 2006 .    c.  sayrin .
2011 .    c.  sayrin , i.  dotsenko , x.  zhou , b.  peaudecerf , th .
rybarczyk , s.  gleyzes , p.  rouchon , m.  mirrahimi , h.  amini , m.  brune , j.m .
raimond , and s.  haroche .
real - time quantum feedback prepares and stabilizes photon number states .
, 477:7377 , 2011 .
doi:10.1038/nature10376 .
a.  somaraju , i.  dotsenko , c.  sayrin , and p.  rouchon .
design and stability of discrete - time quantum filters with measurement imperfections .
, 50845089 , 2011 .
a.  somaraju , m.  mirrahimi , and p.  rouchon .
approximate stabilization of an infinite dimensional quantum stochastic system .
, 25:1350001,2013 .
h.  wang , m.  hofheinz , m.  ansmann , r. c.  bialczak , e.  lucero , m.  neeley , a. d.  oconnell , d.  sank , j.  wenner , a. n.  cleland , and john m.  martinis .
measurement of the decay of fock states in a superconducting quantum circuit .
, 101 : 240401 , 2008 .
h.m . wiseman and a.c .
optimal unravellings for feedback control in linear quantum systems . , 94(7):070405 , 2005 .
h.m . wiseman and g.j .
milburn . . , 2009
the following theorem is just an application of theorem @xmath128 in  ( * ? ? ?
[ thm : app ] let @xmath480 be a markov chain on the compact state space @xmath481 suppose , there exists a continuous function @xmath482 satisfying @xmath483 where @xmath484 is a non - negative continuous function of @xmath485 then the @xmath216-limit set @xmath486 ( in the sense of almost sure convergence ) of @xmath480 is contained by the following set @xmath487      the proof is the following . for all @xmath493 , condition
@xmath494 implies that @xmath495 . taking the sum over all @xmath496
, we get @xmath497 for all @xmath60 and @xmath10 .
since @xmath498 and @xmath499 are @xmath36-function of @xmath233 , these relations read @xmath500 where @xmath501 stands for @xmath502 and the scalar functions @xmath503 depend continuously on @xmath4 and @xmath504 .
+ let us finish the proof by contradiction .
assume that for all @xmath489 , there exists @xmath505 $ ] satisfying @xmath506 , such that @xmath507 , @xmath508 .
take @xmath373 tending towards @xmath509 . since @xmath510 and @xmath511 remain in a compact set , we can assume , up to some extraction process , that @xmath510 and @xmath511 converge towards @xmath512 and @xmath513 in @xmath5 and @xmath514 $ ] .
since @xmath515 , we have @xmath516 .
since latexmath:[\[\begin{gathered } \label{eq : rhocn }    we have by continuity for @xmath373 tending to @xmath509 and for all @xmath60 and @xmath10 : @xmath500 .
thus there exists @xmath518 such that @xmath519 ( see the proof of theorem  [ thm : first ] ) . since @xmath520 ,
for @xmath373 large enough , @xmath521 and thus @xmath522 taking @xmath523 , by assumption  [ assum : impt ] , there exists @xmath10 such that @xmath524 . replacing   in   yields : @xmath525 thus , there exists @xmath526 , such that for @xmath523 and @xmath373 large enough @xmath527 but @xmath528 .
this is in contradiction with @xmath529 as soon as @xmath530 .