galaxy clusters are extensively observed in optical , x - ray and radio bands . in the radio band , a cluster can be observed in the the rayleigh - jeans side of the cosmic microwave background spectrum , as a dip in the brightness temperature , due to sunyaev - zeldovich effect ( sz ) ( sunyaev & zeldovich 1972 ; for a comprehensive review , see birkinshaw 1999 ) .
the sz distortion appears as a decrement for wavelengths @xmath4 mm ( frequencies @xmath5 ghz ) and as an increment for wavelengths @xmath6 mm .
sz effect has the advantage that the sz intensity , unlike that of the x - ray , does not suffer from the @xmath7 cosmological dimming . as discussed by numerous authors ( birkinshaw & hughes 1994 ; silverberg et.al .
1997 ) , one can combine the x - ray and radio observations for clusters to determine cosmological parameters .
this has been done in the recent years to determine the hubble constant @xmath0 ( birkinshaw 1999 ) .
the sz signal is , however , weak and difficult to detect .
recent high signal to noise detection have been made over a wide range in wavelengths using single dish observations : at radio wavelengths ( herbig et al , 1995 , hughes & birkinshaw 1998 ) , millimeter wavelengths ( holzapfel et al . 1997 , pointecouteau et al . 1999 ) and submillimeter wavelengths ( komatsu et al . 1999 ) .
interferometric observations have also been carried out to image the sz effect ( jones et al .
1993 , saunders et al .
1999 , reese et al .
1999 , grego et al .
other than estimating @xmath8 combining sz and x - ray data , sz effect alone can also be used to determine the cosmological mass density @xmath9 of the universe ( bartlett & silk 1994 , oukbir & blanchard 1992 , 1997 , blanchard & bartlett 1998 ) . however
, these procedures generally assume the cluster gas is spherical , unclumped and isothermal
. almost all clusters , however , show departures from these simplistic assumptions with some to a large extent .    departures from these simple assumptions can lead to systematic errors in the determination of the different cosmological parameters ( inagaki et al 1995 ) , especially , it was seen that non - isothermality of the cluster can lead to a substantial error in values of the cosmological parameters .
temperature structure in a cluster can be the result of the shape of the gravitational potential ( navarro et al 1997 , makino et al 1998 ) , or it can arise due to the fact that the initial falling gas in the cluster potential is less shock heated than the later falling gas ( evrard 1990 ) .
in fact hydrodynamical simulations of isolated clusters also show a definite temperature structure and can introduce error in the value of @xmath0 when compared to the traditional isothermal @xmath10-models ( yoshikawa 1998 ) .
roettiger et al ( 1997 ) have shown that cluster mergers can result in deviations from both sphericity and isothermality .
observationally , the main handicap arises from the fact that the thermal structure of clusters are hard to measure , and the temperatures generally taken in analyses are the x - ray emission weighted temperature , usually measured over a few core radii .
so , in general an isothermal description of the cluster is taken ( or sometimes a phenomenological temperature model based on the coma cluster : for example see eq 73 of birkinshaw 1999 ) .    in this paper
, we study another important phenomenon that can substantially change the temperature structure , _ viz . _ a cooling flow .
cooling flows in clusters of galaxies ( for an introduction , see fabian et al 1984 ) is a well established fact by now , and it is seen that around @xmath11 of clusters exhibit cooling flows in their core with @xmath12 of them having cooling flows of more than @xmath13 ( markevitch et al 1998 , peres et al 1998 , allen et al 1999 ) . in the largest systems , the mass deposition rate can be as high as @xmath14 ( allen 2000 ) .
the idealised picture of a cooling flow is as following : initially when the cluster forms , the infalling gas is heated from gravitational collapse . with time
this gas cools slowly and a quasi hydrostatic state emerges .
however in the central region , where energy is lost due to radiation faster than elsewhere , an inward ` cooling flow ' initially arises due to the pressure gradient ( fabian 1994 ) .
this can modify the sz decrement and act as a systematic source of error in the determination of the cosmological parameters .
schlickeiser ( 1991 ) has shown that free - free emission from cold gas in the cooling flow can actually lead to an apparent decrease of the sz effect at the centre .
since the central cooling flow region is generally very small , the isothermal @xmath10-model of cluster gas can still be used for the majority of cluster region even for cooling flow cluster , with the extra precaution of excluding the central x - ray spike from the x - ray fit , and a corresponding change made in the fitting of the sz decrement .
this is only possible for nearby clusters , however with well resolved cluster cores .
naively , the change in the central sz decrement @xmath15 can be seen as follows : for a non cooling flow cluster , the central decrement is given by the line of sight integral of the electron pressure through the cluster centre along the full extent of the cluster .
if the cluster has a maximum radius @xmath16 , then the central sz decrement at rj wavelengths can be written as @xmath17 . for a cluster with a cooling flow ,
let us suppose that the electron pressure @xmath18 , drastically falls below a certain radius @xmath19 , which is typically well inside the core of the cluster .
the resulting central decrement is then @xmath20 .
depending on the distance of @xmath19 from the cluster centre [ @xmath21  to  @xmath22,there wiil be a change in the value of @xmath15 by @xmath23
. however , this simplistic view may not be true .
this estimate assumes that the pressure profile remains a @xmath10 profile outside the radius @xmath19 . the pressure profile , however , need not follow the @xmath10 profile once cooling flow starts and it can deviate from it substantially even for radii much larger than @xmath19 . as a matter of fact , there can actually be an increase in the pressure for a large region inside the cooling flow , before a sudden drop inside @xmath19 .
since the usual proecedure for estimating the hubble constant depends on fitting @xmath10 profiles to the sz and x - ray profiles , to estimate @xmath24 , this change in the pressure profile due to cooling flow can distort the estimation of @xmath24 and hence , the value of @xmath8 in a non - trivial way .
we study this effect in detail in later sections .    in this paper
, we have investigated the problem of cooling flow induced change in the temperature and density profile , its effect on the sz effect , and its subsequent effect on the determination of cosmological parameters . in
2 we briefly review the physics of sunyaev zeldovich effect ;  3 is devoted to the physics of cooling flows and discussing the cooling flow solutions ; in  4 we look at the effect of cooling flow solutions on sz effect , both on the determination of @xmath0 and @xmath9 ; we conclude in  5 with a brief comment on how this work differs from other work and the relevance of this paper . for the sz effect
, our notation and approach mainly follows that described in barbosa et al(1996 ) .
the integral of the electron pressure along any line - of - sight through the cluster determines the magnitude of the distortion of the apparent brightness temperature of the cosmic microwave background ( cmb ) due to sz effect .
this is quantified in terms of the compton @xmath25-parameter : @xmath26 where @xmath27 is the boltzman constant , @xmath28 is the gas temperature , @xmath29 is the electron rest mass , @xmath30 is the electron number density , @xmath31 is the velocity of light and @xmath32 is the thomson scattering cross - section .
this occurs through the inverse - compton scattering , by the hot intracluster gas , of the cmb photons propagating through the cluster medium , and the energy transfer in this interaction between hot electrons and cmb photons resulting in a distortion to the cmb spectrum .
the sz surface brightness at a position @xmath33 of the cluster with respect to the mean cmb intensity is given by @xmath34    _ x _ is a dimensionless frequency parameter @xmath35 where @xmath36 is the planck constant , @xmath37 is the observing frquency and @xmath38 is the cmb temperature at the present epoch : @xmath39k .
the function @xmath40 describes the spectral shape of the effect @xmath41.\ ] ] since the total photon number is conserved in the inverse compton scattering process , upscattering of the photons , the spectral dependance gets an unique shape , through a decrement in the brightness temperature at lower frequencies while an increase is observed at higher frequencies .
the sunyaev - zeldovich effect provides a unique observing approach to traditional methods  which use x - ray temperature and x - ray luminosity .
the x - ray studies have a major disadvantage because of ` cosmological dimming ' , the surface brightness of distant x - ray sources falls off as @xmath42 , and for this reason , obtaining samples of clusters at cosmological distances is challenging .
the sz effect has the distinct advantage of being independent of the distance to the cluster .
the sz flux density from a cluster will diminish with distance to the cluster as the square of the angular - size distance ; in contrast to x - ray flux densities from clusters which diminish as the square of the luminosity distance to the cluster .
if the observations resolve clusters , particularly at lower redshifts , the observed sky sz temperature distribution will be sensitive to the thermal electron temperature structure within the clusters ; once again this may be contrasted with x - ray emission images of cluster gas distributions which are mainly sensitive to the gas density distribution .
the method for the determination of hubble constant using sunyaev - zeldovich effect uses two observable quantities : 1 ) @xmath43 of the cmb due to sz effect ; 2 ) the x - ray surface brightness @xmath44 of the cluster .
these can be written as ( r ) = -2_l_min^l_maxk_b t_e _ t n_e dl , s_x(r ) = 1 _ l_min^l_max dl_x dl , where @xmath45 is the distance to the line of sight from the cluster centre , @xmath46 and @xmath47 give the extension of the cluster along the line of sight , @xmath48 is the x - ray emissivity and @xmath49 the line element along the line of sight .
the x - ray emissivity in the frequency band @xmath50 to @xmath51 can be written as = n_e^2 ( t_e;_1 , _ 2 , z ) , where ( t_e;_1 , _ 2 , z ) = 2 ^ 1/2 16e^6 a ( t_e;_1 , _ 2 , z ) , where @xmath52 du .\end{aligned}\ ] ] in the above equations we have assumed primordial abundance of hydrogen and helium and have set @xmath53 , @xmath54 is the electron charge , @xmath55 , @xmath56 , and @xmath57 is the velocity averaged gaunt factor for the ion of charge ze ( kellog 1975 ) .
traditionally , to model the cluster gas distribution one takes the following density and temperature profiles ( cavaliere & fusco - femiano 1978 ) n_e(r ) = n_e0 ^ -3/2 , t_e(r ) = t_iso = constant , where @xmath58 is the central electron density and @xmath24 the core radius of the cluster .
the above expressions are used as an empirical fitting model , and the parameter @xmath59 is regarded as the fitting parameter .
the equation holds for @xmath60 , where @xmath61 is the maximum ` effective ' extension of the cluster .
conventionally , @xmath62 , and then from equations ( 5),(6),(10),(11 ) , we get @xmath63}^ { 1/2 - 3\beta /2 } , \end{aligned}\ ] ] @xmath64}^{1/2 - 3\beta } , \end{aligned}\ ] ] where @xmath65 is the gamma function . since both the central cmb decrement and the x - ray surface brightness are observed , one can then combine equations ( 12 ) and ( 13 ) to estimate the core radius as r_c , est & = & ^2_obs ( 3 - 1/2)(3/2)^2 + & & m_e^2 c^4 ^-1/2 . in the above equation @xmath66
is the angular core radius observed in x - ray , and @xmath67 is the x - ray flux averaged temperature ( obtained from fitting the observed x - ray spectrum to the theoretical spectrum expected from isothermal case ) .
this x - ray emission weighted temperature is given by t_iso .
the point to be noted is that @xmath68 is the virial radius of the cluster and its choice depends on the observer . if the temperature has a spatial structure then the @xmath69 inferred from such a procedure may give different values depending on how much of the cluster is taken in making the above average . it has been seen ( yoshikawa et al 1998 ) , that this can lead to a substantial change in sz effect inferred , and thus to the value @xmath0 .
the angular diameter distance @xmath70 can be approximated for nearby ( @xmath71 ) clusters as d_a = cz .
thus finally we have got an estimate of @xmath72 as @xmath73 .
now , if from other observations we know the cosmological parameters @xmath9 and @xmath74 , then one can estimate the hubble constant    as can be seen from the above equations , the value of @xmath0 depends crucially on the many assumptions of isothermality and @xmath10-model density distribution of the cluster . cooling flow changes both of these and so it can significantly affect the value of the hubble constant .
from x - ray spectra of clusters , it is known that the continuum emission is thermal bremsstrahlung in nature and originates from diffuse intracluster gas with densities @xmath75 and temperatures around @xmath76 k. the gas is usually believed to be in hydrostatic equilibrium . if , however , the density in the inner region is large enough , so that the cooling time is less than the age of the cluster , then there is a ` cooling flow ' ( fabian et al 1984 and references therein ) . of course there would be a flow only when the dynamical time is also shorter than the cooling time scale ( @xmath77 ) .
the basic equations of cooling flow are : & = & -u _ b d _ b dr - 2 u r - _ b , + udu+ 1 d(_b ) & = & - gm_t(r ) , + ud(3 ) - u d_b & = & ( ) , [ eq : flow ] where @xmath78 , @xmath79 is the mean molecular weight , and @xmath80 is the proton mass . for steady flows with constant mass flux ,
this implies @xmath82 for steady flows .
( note that in cooling flows both @xmath83 and @xmath84 are negative .
the subscripts @xmath85 refers to baryons and @xmath86 refers to total i.e baryons + dark matter .
however , we assume the baryonic contribution to the total mass negligible w.r.t to the dark matter contribution . )
@xmath87 describes the distribution of the total mass and depends on the details of dark matter density profiles ( see below ) .
@xmath88 is the cooling function defined so that @xmath89 is the rate of cooling per unit volume .
we use an analytical fit to the optically thin cooling function as given by sarazin & white ( 1987 ) , & & ( ( ) 10 ^ -22 erg cm^3 s^-1 ) = 4.7 + & + & 0.313 t^0.08 + & + & 6.42 t^-0.2 + & + & 0.000439 t^0.35 .
this fit is accurate to within @xmath90% accuracy , for a plasma with solar metallicity , within @xmath91 k. for @xmath92 k , it underestimates cooling by a factor of order unity ( compared to the exact cooling function , as in , e.g. , schmutzler & tscharnuter , 1993 ) , and therefore is a conservative fit to use , as far as the effect of cooling is concerned .    for non - steady flows , we adopt the formalism of white & sarazin ( 1987 ) , where the mass deposition rate , @xmath93 , is characterised by a ` gas - loss efficiency ' parameter @xmath94
. one writes @xmath95 where @xmath96 is the local isobaric cooling rate ( @xmath97 ) .
it has been found that @xmath98 models can reproduce the observed variation of mass flux ( @xmath99 ) ( sarazin & graney 1991 ) .
fabian ( 1994 ) has noted that these models of white & sarazin ( 1987 ) yield good approximations to the emission weighted mean temperature and density profiles for cooling flow clusters .
we also note that rizza et al ( 2000 ) have used the steady flow models of white & sarazin ( 1987 ) to simulate cooling flows .
we first discuss cooling flows with @xmath100constant . with @xmath101
, one can eliminate the density from eq .
[ eq : flow ] to get two differential equations :    & = & u   + ddr & = & 2 [ eq : steady ] these equations have singularities at the sonic radius @xmath19 where @xmath102 .
a necessary condition of singularity is that the numerators of eq .
[ eq : steady ] vanish at the sonic radius
. therefore ( mathews & bregman 1978 )    r_s = ( 3/10_s )    we have used two different dark matter profiles for the cluster .
the first model ( model a ) has been discussed earlier in the literature in the context of cooling flows in cluster ( white & sarazin 1987 ; wise & sarazin 1993 ) with a density profile , _
d=\ {    ll _ o 1+(r / r_core)^2 + _ o , g 1+(r / r_c , g)^2 & + _ o 1+(r / r_core)^2 &    . here
@xmath103 gm @xmath104 and @xmath105 kpc describe the profile of the cluster mass , and @xmath106 gm @xmath104 and @xmath107 kpc describe the profile of the galaxy in the centre of the cluster .
model b does not have the galaxy in the center , and so it is described simply by @xmath108 $ ] .    with these assumptions ,
the solutions for steady cooling flows , @xmath100 constant ) are fully characterized by ( 1 ) the inflow rate , @xmath84 , and ( 2 ) the temperature of the gas @xmath109 at the sonic radius @xmath19 .
obviously , the cooling flow solutions are only valid within the cooling radius @xmath110 where @xmath111 .
we assume a value of @xmath112 gyr for all models .
we assume that outside the cooling radius , gas obeys quasi - hydrostatic equilibrium ( sarazin 1986 ) .
although this means matching the cooling flow solutions with nonzero @xmath83 to @xmath113 solutions outside , in reality the velocity of gas at the cooling radius is very small ( for a @xmath114 m@xmath115 yr@xmath116 with @xmath117 gm @xmath104 , at @xmath118 kpc implies a velocity of @xmath119 km s@xmath116 ) , which is close to the limit of turbulence in the cluster gas ( jaffe 1980 ) , and smaller than the sound velocity ( @xmath120 km s@xmath116 ) .
the velocity of the flow at the cooling radius is , therefore , for all practical purposes , sufficiently small to be matched to the solution of hydrostatic equilibrium outside .
( in this approach , we avoid the time consuming search for the critical value of @xmath84 for which the flow solutions behave isothermally at @xmath121 ( see sulkanen 1993 ) . )    as in the usual assumptions for the interpretation of sz effect , we assume that the gas outside the cooling radius is isothermal , with a constant temperature profile .
the density , therefore , obeys @xmath122^{-3 \beta /2}$ ] , where @xmath123 , and @xmath69 is the temperature of the gas at and outside the cooling radius .    for models with non - zero @xmath94 ( model c : has the same mass profiles as model a )
, the solutions are characterized by @xmath109 and the value of @xmath84 at the cooling radius , @xmath124 .
since a fraction of mass drops out of the flow in this case , the inflow velocity need not rise fast and so it is possible to find completely subsonic solutions .
we numerically solved the flow equations for the parameters listed in table 1 .
the density , temperature and pressure profiles for three cases are presented in figures 1 , 2 & 3 .
we mark the position of @xmath110 in each case , and we mark @xmath19 for the cases of transonic flows ( when @xmath125 constant ) . beyond @xmath110
we match a hydrostatic solution , as explained above , for the respective potentials .
we also present , for comparison , the behaviour if the solutions outside @xmath110 are assumed to extend inwards ( that is , if no cooling flow is assumed ) .
we will postpone the discussion on the effect of these profiles on the sz decrement to a later section , and only discuss the qualitative aspects of the solutions here .
the solution a1 is similar to that presented by wise & sarazin ( 1993 ) ( their figure 1 ; although they chose to characterise the solutions by the temperature at @xmath110 and not @xmath109 as we have done here ) .
it is also similar ( qualitatively ) to the solution for a262 presented by sulkanen ( 1989 ) .
as the latter authors have noted , the effect of having a galactic potential in the center is to have a flatter temperature profile for @xmath126 , than in the case of no galactic potential .
this aspect is clearly seen while comparing our solutions with and without galactic potentials in the center .
our calculations for the case without the central galaxy are admittedly flawed in the very inner regions where the gas density is larger than the dark matter density , which results in an incorrect determination of the graviational potential in the inner region .
however , this happens only inside a region @xmath127 kpc from the centre , and should not influence our final results by a large extent .
a word of explanation for the pressure profiles is in order here . naively speaking
, it would appear that the pressure profile inside the cooling radius should have lower values than the corresponding case of hydrostatic equilibrium .
the fact that it is not always so has been noted in the literature ( e.g. , soker & sarazin 1988 , fig 1 of sulkanen 1989 ) .
the reason for the pressure bump just outside of the sonic radius is that the flow in this inner region is not pressure driven , but rather by gravity ( see also soker & sarazin 1988 ) .
this is why the bump in the profile depends on the presence and absence of the galaxy in the centre .
and this profile leads to the curious result that the presence of cooling flow can lead to the overestimation of the hubble constant as discussed in the next section    the model with mass deposition ( c1 ) is shown in figure 3 .
the local mass flux is found to be almost proportional to the radius , consistent with various observations ( fabian , nulsen & canizares 1984 ; thomas 1987 ) , and , therefore , is probably a realistic model for cooling flow clusters . in this case
the temperature drops gradually all the way through , since the velocity does not rise too fast .
the deposited mass is assumed to impart negligible pressure and the pressure refers only to that of gas taking part in the flow .
in this section we discuss the sz and x - ray profiles of clusters with cooling flows . we compare these with profiles from cluster having gas in hydrostatic equilibrium , and comment on the reliability of measuring hubble constant .
the effect of cooling flow and the subsequent increased bremsstrahlung emmison is seen in the sudden increase in the x - ray flux in the innermost region of the cluster ( figure 5 ) .
the signature of the cooling flow is seen in the form of the central spike in the x - ray profile .
the x - ray profile is only affected slightly by the drop in temperature and it is the dependence on the gas density that holds . the temperature dependance becomes important only near the sonic point .
outside @xmath110 , the x - ray profile is the same as that in the hydrostatic cases .
the sz distortion is proportional to the line of sight integral of the pressure , and the sudden increase of the gas density inside the cooling radius is moderated by the decrease of the gas temperature . as a result
there is a gradual increase in the gas pressure .
near the sonic point the temperature drops drastically by orders of magnitude , and results in sudden decrease in pressure .
however , since this change in pressure becomes acute only within @xmath128 of the core radius , it contributes negligibly to the line - of - sight integral of the gas pressure , and leads to an increase in the sz distortion inside the cooling radius for all models considered ( see figure 4 ) . like the x - ray profiles ,
the sz profiles outside @xmath110 is the same as that for the corresponding hydrosstatic cases .
the sz profiles have been calculated in the rayleigh - jeans limit ( @xmath129 ) where @xmath130 of eqn .
4 goes to @xmath131 .
in general , however , the profiles should be calculated using eqn 2 .
our results below are independent of the observational frequency , since the profiles at different frequencies have similar shapes , with the amplitude of the sz distortion scaled either up or down .
once both profiles are known , one can determine the deviation in the value of the hubble constant using eqns(14 ) & ( 16 ) .
the deviation from the idealistic case can be parametrised as    f_h = h_0,esth_0,true the above formula has been used to determine the deviation of the estimated value of @xmath0 from the actual value , for models listed in table-1 .
the effect of cooling flow on the determination of the cosmological parameters are summarised in table-2 .    to begin with ,
one has to get best fitted values for @xmath24 ( or @xmath132 ) and @xmath10 from different profiles . since , the estimation of the hubble constant depends on the determination of these parameters from the profiles , we look at this issue in more detail .
we must keep in mind that the best fitted value of @xmath24 ( or @xmath132 ) and @xmath10 depends on whether one decides to fit the x - ray or the sz profiles , and the choice can lead to significant differences in the estimated value of @xmath0 .
one of the reasons for the strong dependance on the nature of the profile can be the non - isothermality of the cluster gas .
recent observations indicate that intracluster gas has a temperature structure , see markevitch et al 1998 .
this is because the y - parameter depends on the integral over @xmath28 , while emissivity of thermal bremsstrahlung depends on @xmath133 .
the dependance of the gaunt factor on @xmath28 isindirect and weak .
yoshikawa et al 1998 have shown that gas temperature drop in the central regions ( their fig .
3 ) , should increase both @xmath24 and @xmath10 fitted to @xmath134 , and to alesser extent to @xmath135 , as compared to those compared to @xmath136 .
this discrepancy increases at higher redshifts .
however , in their case , there is little change in the gas density profile .
clumpiness can also give rise to different fits , resulting in an overestimation of the hubble constant ( inagaki etal 1995 ) .
there are two other important points that have to be kept in mind while fitting the profiles .
first , we must remember that we are trying to fit a cluster having a finite profile with the formulae ( eqns 12 & 13 ) for isothermal @xmath10 profiles which is derived assuming the cluster to be of infinite extent .
this can , by itself , lead to an overestimation of @xmath0 ( inagaki et al 1995 ) . thus
to have a good fit one must choose a segment of the profile such that , within that segment , the profiles ( sz or x - ray ) for a finite cluster do not differ much from those of a _ hypothetical _ cluster of infinite extent .
we found that sz and x - ray profiles of clusters start differing from those of infinite size at radii greater than @xmath137 times the core radius .
hence , we have restricted our fitting to radii within @xmath138 .
next , one must also be careful to exclude the region close to the sonic point , so that the x - ray spike is excluded from the fit .
also , the central portion in the sz profile should be avoided as its inclusion can give an _ apparent _ central decrement less than its neighbouring points ( see schlickeiser 1991 ) .
we have fitted the sz and x - ray profiles varying the inner cutoff radius and the results for a representative solution for each class of model are tabulated in tables 3 & 4 .
thus , all fittings were done for profiles extending from @xmath139 to @xmath140 .    as can be seen from table-2 ,
cooling flows can lead to an overestimation of the hubble constant .
however , we must emphasise , that it may not be possible to _ a priori _ estimate the amount of bias introduced in the measurement of the hubble constant due to cooling flows .
there is no simple correlation between the amount of cooling ( i.e @xmath84 ) and the change in the estimated @xmath0 from the actual value .
the total change depends not only on @xmath84 , but also on the position of cooling radius , sonic radius , temperature at the sonic point and the isothermal temperature characterising the hydrostatic cases , with which comparisons are made .
specifically , the fitted values of @xmath141 and @xmath10 for cooling flow models differ from hydrostatic models according to shape of underlying profiles , which is marked by two important features , firstly the central excess of x - ray flux ( or excess decrement of sz flux ) , and secondly the the deviation from the smooth hydrostatic profile inside @xmath110 , the amount of overestimation mainly depends on these factors . for models with a central galaxy potential
, there is always an over - estimation of @xmath0 , which is greater than the models without the central galaxy .
for the realistic cases of models c1 and c2 , where we have a variable @xmath84 with @xmath45 inside the cooling radius , the deviation of estimated hubble constant from its actual value is almost the same .
they are also greater than that of models a and b , having similar mass flow rates .
this may be due to the fact that the maximum deviation in pressure from the hydrostatic cases is more in non steady cases , than in steady flows . also non - steady cases are marked by the absence of the sonic radii and the subsequent drop in temperature .
we note that although the different choice of fitting may change the absolute determination of cosmological parameters , the trend i.e deviation from the correct values remains more or less unaffected .
it is interesting to note that for b - type model ( c1 ) , which include mass deposition in cooling flows , the deviations decrease as one excludes a greater part of the cooling flow region ( table 4 & 5 ) .
the other models show an increase instead . here , we remind ourselves that models with mass deposition i.e c - type models are more realistic ( fabian 1994 ) .
it is possible that the unsually high value of deviation ( table 4 ) and the counter - intuitive trend of increasing deviation with decreasing portion of cooling flow region used for fitting ( table 4 & 5 ) , arise because of the unrealistic modelling of cooling flows .
if we take the model b1 as a realistic one , then table 4 & 5 show that to obtain a value of the hubble constant within an accuracy of @xmath2 , one should have @xmath142 . in most cases , @xmath143 ( fabian 1984 ) .
however , as @xmath110 can not be determined without actually detecting a cooling flow in a cluster , we suggest that _ a significant portion of the profile within @xmath24 should be excluded as a precaution_. the sz and the x - ray profiles for the different models are shown in figures 4 & 5 .
our work on the effect of temperature structure of clusters and its effect on sz decrement differs from other previous work of this nature in following way : this work takes into account the change in density profile as well as the temperature profile since both becomes important in the central region of the cluster .
also , previous authors have looked at the issue of non - isothermality of a cluster at radii greater than the core radius of the cluster , whereas we look at temperature change at regions inside the core radius . for them
the density profile can still be well approximated by a @xmath10 profile , whereas for cooling flow solutions , density profile is vastly different .
further , they have neglected radiative cooling in their work .
we for the first time look at sz effect in presence of radiative cooling , by first solving the cooling flow equations for reasonable and physical solutions .    in summary , we find that the presence of a cooling flow in a cluster can lead to an overestimation of the hubble constant determined from the sunyaev - zeldovich decrement .
we have used different models of cooling flows , with and without mass deposition , and found the deviation in the estimated value of the hubble constant in the case of a cooling flow from that of hydrostatic equilibrium .
we have used the usual procedure of fitting the sz and x - ray profile with a @xmath10 profile to get an estimated value of @xmath24 , and then compared with that for the case of gas in hydrostatic equilibrium in order to estimate the deviation in the hubble constant . for the more realistic models with mass deposition ( varying @xmath1 with radius )
, we found that the deviation decreases with the exclusion of greater portions of the cooling flow region .
quantitatively , we found that for the deviation to be less than @xmath2 , one should exclude a portion of the profile upto @xmath144 . since @xmath110 is difficult to estimate without actually detecting a cooling flow , we have suggested that a significant portion of the profile inside @xmath24 should be excluded , to be safe .
there can be another important implication of the effect of cooling flows . with the upcoming satellite missions ( map & planck ) , we have come to the point
where there are efforts to constrain @xmath145 with surveys of blank sz fields ( bartlett et al 1998 ; bartlett 2000 ) , ultimately giving rise to sz - selected catalog of clusters ( aghanim et al 1997 ) .
this method relies on estimating the number of sz sources brighter than a given threshold flux ( barbosa et al 1996 ) .
the point to be noted is that since these surveys are essentially flux limited in nature , the validity of the analysis to determine @xmath145 depends crucially on the _ one - to - one _ association of flux - limits to corresponding mass limits of clusters . from our analysis above , it seems that it may not be possible to associate a unique cluster mass to a given sz distortion , given the uncertainty due to the presence of cooling flows .
this might lead to contamination in sz cluster catalogs and the inference of @xmath145 .
recently , attempts have been made to constrain @xmath9 from variance measurement of brightness temperature in blank fields ( subrahmanyan et al 1998 ) and comparing them to simulated fields ( majumdar & subrahmanyan , 2000 ) , of cumulative sz distortions from a cosmolgical distribution of clusters .
these results may also be systematically affected due to the presence of clusters having cooling flows .
the estimations made in this paper strictly applies to cases where the image of the sz effect is directly obtained by single dish observations . for interferometric observations , the interferometer samples the fourier transform of the sky brightness rather than the direct image of the sky .
the fourier conjugate variables to the right ascension and declination form the @xmath146 plane in the fourier domain . due to spatial filtering by an interferometer ,
it is necessary that models be fitted directly to the data in the @xmath146 plane , rather that to the image after deconvolution .
we do not forsee drastic change from our inferences in such cases since the result mainly depends on the deviation of the sz and x - ray profile in case of a cooling flow from those in hydrostatic equilibrium .
this , however , should be looked in greater detail in future .
we also note that with the growing number of high quality images of the sz effect with interferometers , which have greater resolution than single dish antennas , the shape parameters of the clusters can be directly determined from the sz dataset rather than from an x - ray image ( grego et al .
2000 ) .    finally , we would like to add , though the calculations presented in this paper were done using the dark matter profile ( eqn .
21 ) , which is `` commonly used '' for calculating cooling flow solutions , it is inconsistent with the dark matter profile ( navarro et al 1997 ) found in numerical simulations .
( for a comparison of mass and gas distribution in clusters having cooling flows with different dark matter profiles , see waxman & miralda - escude , 1995 ) .
moreover , we have neglected the self gravity of the gas .
suto et al ( 1998 ) have calculated the effect of including the self gravity of the gas in determining the gas density profile . to make strong conclusions about the effect of cooling flows in the determination of the hubble constant
, one should take both the points mentioned above into account .
we are grateful to joseph silk for his comments on the manuscript .
we also thank the referee , dr asantha cooray , for numerous suggestions which have helped us in improving the paper .
one of the authors ( sm ) would like to thank the raman research institute , bangalore , for hospitality and for providing computational facilities .
sm acknowledges help recieved from r.sridharan on using idl .
finally , he would also like to express his gratitude to pijushpani bhattacharjee for his constant encouragement .
aghanim , n. , de luca , a. , bouchet , f. r. , gispert , r. & puget , j. l. , 1997 , a & a , 325 , 9 allen , s. w. , fabian , a. c. , johnstone , r. m. , arnaud , k. a. , & nulsen , p. e. j. , 1999 : astro - ph/9910188 allen , s. w. , 2000 : astro - ph/0002506 barbosa , d. , bartlett , j. g. , blanchard , a. & oukbir , j. , 1996 , astro - ph/9607036 barbosa , d. , bartlett , j. g. , blanchard , a. & oukbir , j. , 1996 , a & a , 314 , 13 bartlett , j. g. & silk , j. , 1994 , apj , 423 , 12 bartlett , j. g. , blanchard , a. & barbosa , d. , astro - ph/9808308 bartlett , j. g , . 2000
: astro - ph/0001267 birkinshaw , m. & hughes , j. p. , 1994 ,
apj , 420 , 33 birkinshaw , m. 1999 , phys rep .
, 310 , 97 blanchard , a & bartlett , j. g. , 1998 , a & a , 332 , 49l cavaliere , a. & fusco - femiano , r. , 1978 , a & a , 70 , 677 evrard , a. e. , 1990 , apj , 363 , 349 fabian , a. c. , nulsen , p. e. j. & canizares , c. r. 1984 , nature , 310 , 733 fabian , a. c. 1994 , araa , 32 , 277 grego , l. , carlstrom , j. e. , joy , m. k. , reese , e. d. , holder , j. p. , patel , s. , cooray , a. , & holzapfel , w. l. , 2000 , : astro - ph/0003085 herbig , t. , lawrence , c. r. , redhead , a. c. s. , & gulkis , s. , 1995 , apj , 449 , l5 holzapfel , w. l. , et al .
1997 , apj , 480 , 449 reese , e. d. , 1999 , preprint ( astro - ph/9912010 ) hughes , j. p. , & birkinshaw , m. , 1998 , apj , 501 , 1 inagaki , y. , suginohara , t. & suto , y. , 1995 , pasj , 47 , 411 jaffe , w. 1980 , apj , 241 , 925 jones , m. et al .
1993 , nature , 365 , 320 kellog , e. , baldwin , j. r. & koch , d. , 1975 , apj , 199 , 299 komatsu , e. , kitayama , t. , suto , y. , hattori , m. , kawabe , r. , matsuo , h. , schindler , s. , & yoshikawa , k. , 1999 , apj , 516 , l1 majumdar , s. & subrahmanyan , r. , 2000 , mnras , 312 , 724 makino , n. , sasaki , s. & suto , y. , 1998 , apj , 497 , 555 markevitch , m. , forman , w. r. , sarazin , c. l. & vikhlinin , a. , 1998 , apj , 504 , 27 mathews , w. g. , & bregman , j. n. , 1978 , apj , 224 , 308 navarro , j. f. , frenk , c. s. , & white , s. d. m. , 1997 , apj , 40 , 493 oukbir , j. & blanchard , a. , 1992 , a & a , 262 , 21l oukbir , j. & blanchard , a. , 1997 , a & a , 317 , 1 padmanabhan , t. , 1993 , structure formation in the universe , cambridge university press , cambridge peres , c. b. , fabian , a. c. , edge , a. c. , allen , a. w. , johnstone , r. m. , & white , d. a. , 1998 , mnras , 298 , 416 pointecouteau , e. , giard , m. , benoit , a. , dsert , f. x. , aghanim , n. , coron , n. , lamarre . , j. a. , & delabrouille , j. , 1999 , apj , 519 , l115 press , w. h. & schecter , p. , 1974 ,
apj , 187 , 225 reese , e. d. , et al .
1999 : astro - ph/9912071 rizza , e. , loken , c. , bliton , m. , roettiger , k. , & burns , j. , 2000 , apj , 119 , 21 roettiger , k. , stone , j. m. & mushotzky , r. f. , 1997 , apj , 482 , 588 sarazin , c. l. 1986 , rev .
phys . , 58 , 1 sarazin , c. l. & white , r. e. 1987 , apj , 320 , 32 sarazin , c. l. & graney , c. m. 1991 , apj , 375 , 532 saunders , r. et al .
1999 : astro - ph/9904168 schlickeiser , r. 1991 , a & a , 128 , 23l schmutzler , t. & tscharnuter , w. m. 1993 , a & a , 273 , 318 silverberg , r. f. et al , 1997 , apj , 485 , 22 soker , n. & sarazin , c. l. 1988 , apj , 327 , 66 subrahmanyan , r. , kesteven , m. j. , ekers , r. d. , sinclair , m. & silk , j. 1998 , mnras , 298 , 1189 sulkanen , m. e. , burns , j. o. & norman , m. l. 1989 , apj , 344 , 604 sunyaev , r. a. & zeldovich , ya . b. 1972 , comm .
space phys.,4 , 173 suto , y. , sasaki , s. & makino , y. , 1998 , apj , 509 , 544 thomas , p. , fabian , a. c. & nulsen , p.e.j .
1987 , mnras , 228 , 973 .
waxman , e. & miralda - escude , j. , 1995 , apj , 451 , 451 white , r. e. , & sarazin , c. l. , 1987 , apj , 318 , 629 wise , m. w. & sarazin c. l. , 1993 , apj , 415 , 58


building scalable multi - qubit systems is presently the main challenge towards the implementation of a solid - state quantum information processor  @xcite .
high - fidelity single qubit gates based both on semi-  @xcite and super - conducting technologies are nowadays available .
in particular , for the three basic types of superconducting qubits ( charge , flux and phase ) single - qubit operations with high quality factor have been demonstrated in different laboratories  @xcite .
further improvement has been recently achieved with circuit - qed architectures  @xcite .
multi - qubit systems instead have been proved harder to operate , the main limitation arising from the broadband solid state noise .
the requirements for building an elementary quantum processor are in fact quite demanding on the efficiency of the protocols .
this includes both severe constraints on readout and a sufficient isolation from fluctuations to reduce decoherence effects .
based on experiments on the different josephson junction ( jj ) setups , there is presently a general consensus on the most common noise spectral forms and on the main consequences on systems evolutions .
typically noise is broadband and structured , i. e. the noise spectrum extends to several decades , it is non - monotonic , sometimes a few resonances are present .
noise with @xmath0 spectrum is common to virtually all nanodevices .
its physical origin varies from device to device and depends on the specific material .
different implementations are in fact more sensitive either to charge , flux or to critical current fluctuations with spectral density scaling as the inverse of the frequency .
the presence of slow components in the environment makes the decay of the coherent signal strongly dependent on the experimental protocol being used  @xcite .
measurements protocols requiring numerous repetitions are particularly sensitive to the unstable device calibration due to low - frequency fluctuations .
the leading effect is defocusing of the measured signal , analogous to inhomogeneous broadening in nmr  @xcite .
the intrinsic high - frequency cut - off of @xmath0 noise is hardly detectable , measurements typically extending to @xmath1  hz ( recently charge noise up to @xmath2  mhz has been detected in a set  @xcite ) .
incoherent energy exchanges between system and environment , leading to relaxation and decoherence , occur at typical operating frequencies ( about @xmath2 ghz ) .
indirect measurements of noise spectrum in this frequency range often suggest a white or ohmic behavior  @xcite .
in addition , narrow resonances at selected frequencies ( sometimes resonant with the nanodevice relevant energy scales ) have being observed  @xcite . in certain devices
they originate from the circuitry  @xcite and may eventually be reduced by circuit design .
more often , resonances are signatures of the presence of spurious fluctuators which also show up in the time resolved evolution , unambiguously proving the discrete nature of the noise sources  @xcite .
such fluctuators may severely limit the reliability of nanodevices  @xcite .
explanation of this rich physics is beyond phenomenological theories describing the environment as a set of harmonic oscillators . on the other side , an accurate characterization of the noise sources
might be a priory inefficient , since a huge number of parameters would be required for a microscopic description .
therefore one may follow a different route , consisting in classifying noise sources on the basis of their effects instead of on their nature and to understand , case by case , which is the efficient description of the environment .
the required information may in fact depend on the specific protocol , especially when the environment is long - time correlated .
this program is meaningful for quantum information where relevant time scales are much smaller than the decoherence time .
this means that in the generic favorable situation coupling with the environment has simple effects on the system dynamics . as a difference with condensed matter physics where
long time behavior is emphasized and interesting problems involve entanglement of a many - body system .
here we illustrate a road - map to treat broadband noise which allows to obtain reasonable approximations by systematically including only the relevant information on the environment , out of the huge parametrization needed to specify it .
the multi - stage approach for the different classes of broadband noise has been originally introduced for single qubit gates  @xcite .
the obtained predictions for the decay of the coherent signal are in agreement with observations in various jj implementations and in different protocols  @xcite .
we mention the observed decay of ramsey fringes in charge - phase qubits  @xcite and recent results on flux / phase qubits  @xcite . here
we extend the procedure to complex solid - state architectures . as an illustrative case
, we perform a systematic analysis of the effects and interplay of low and high frequency noise components in a two - qubit gate in a fixed coupling scheme .
such a systematic analysis points out that efficient operations in the solid state require an accurate preliminary characterization of the noise spectral characteristics and tuning appropriately the device working point .
the paper is organized as follows : in section 2 we introduce the general framework of solid - state multi qubit gates and illustrate the characteristics of the most common spectral forms . in section 3
we propose a classification of the noise sources and in section 4 we present a multi - stage approach to deal with broadband noise .
section 5 is dedicated to adiabatic noise . in section 6 the multistage approach
is applied to a universal two - qubit gate .
section 7 summarizes our main findings .
a multi - qubit gate can be modeled by the following hamiltonian ( @xmath3 ) @xmath4 where for each qubit , labeled by @xmath5 , @xmath6 includes set of parameters intrinsic to the device and time dependent classical control fields .
manipulation of tunable fields allows control of the dynamics and design of arbitrary single qubit operations via unitary transformations .
interactions among qubits and the needed additional parametrization to describe the multi - qubit gate are included in @xmath7 .
depending on the design , the control hamiltonian may span a reduced subspace of the qubit liouville space .
in other words , there is a limited number of ports available for control , for tuning , for state preparation and measurement .
noise is also coupled via these ports .
it is rather usual that only one of the control fields allows the required fast addressing of each qubit .
therefore we focus on a model where both the control fields , @xmath8 , and the environment are coupled to a single qubit operator , say @xmath9 , @xmath10 where the polar angles @xmath11 define qubit-@xmath12 working point .
the environment hamiltonian is @xmath13 and @xmath14 is a collective environment variable acting on qubit @xmath12 .
model ( [ h ] ) implies a projection of the device hamiltonian onto the subspace spanned by the two lowest energy eigenstates for each qubit .
this description is valid provided that manipulations with the control fields do not induce leakage to higher energy states of the device  @xcite .
we leave unspecified the nature of the noise sources described by @xmath15 , relevant cases being either discrete or gaussian fluctuations .
in addition , the environment hamiltonian may include correlations among noise sources affecting each qubit . in the spirit of the present analysis
, we assume that the only information at our disposal is the power spectrum of @xmath16 fluctuations @xmath17 where @xmath18 denotes the equilibrium average with respect to @xmath13 .
we remark that the effect of the environment on the system dynamics is completely characterized by @xmath19 only if the environment is composed of harmonic oscillators or if it is weakly coupled to the system and short time correlated .
of course this is not the case if low - energy excitations determine memory effects .
this is the typical situation in the solid state where in general additional statistical information on the environment is required .
we assume that each noise component @xmath16 has broadband spectrum @xmath20 , @xmath21 followed by a white or ohmic flank at frequencies @xmath22 .
low- and high - frequency cutoffs depend on the specific setup
. impurities of various origin responsible for random telegraph fluctuations contribute lorentzian peaks to the power spectrum @xmath23 .
such a spectrum originates from classical fluctuations @xmath24 , where @xmath25 randomly switches between two values @xmath26 with rate @xmath27 .
an ensemble of @xmath28 bistable impurities with a distribution of switching rates @xmath29 , @xmath30 $ ] and average coupling strength @xmath31 , gives rise to @xmath0-spectrum , @xmath32\,\,\omega^{-1}$ ]  @xcite
. selected lorentzian peaks may be visible in the spectrum if individual impurities are strongly coupled , i. e. when @xmath33 . instead , damped coherent fluctuators at frequencies @xmath34 , in the simplest cases contribute to the total spectrum with additional peaks , @xmath35 .
spectroscopic evidence of coherent impurities of frequency @xmath36 close to the qubit larmor frequency has been reported in @xcite .
remarkably , the possibility to exploit spurious quantum two - level systems as qubits  @xcite or for quantum memory operations has been recently demonstrated  @xcite . in these cases a quantum description of the impurity
is required  @xcite .
the above description illuminates that in the solid state we have to deal with broadband and structured noise . in other words
, the noise spectrum extends to several decades , it is non - monotonic , sometimes a few resonances are present .
the various noise sources responsible for the above phenomenology have a qualitative different influence on the system evolution .
this naturally leads to a classification of the noise sources according to the effects produced rather than to their specific nature .
the effects of high - amplitude noise at low frequencies , like @xmath0 noise , vary from protocol to protocol .
this feature is typical of non - markovian baths .
quantum operations necessarily require repetitions of single detections , each leading to boolean answer
. therefore , even `` single shot '' measurements result from numerous repetitions of single runs in an overall process which may last minutes . in the presence of low frequency fluctuations
this leads to unstable device calibration and random clock frequencies in the various repetitions . as a result ,
a de - focused signal is observed , a phenomenon analogous to inhomogeneous broadening in nmr .
re - focusing protocols , like echo or some dynamical decoupling schemes , allow partial recovery of the signal  @xcite .
since the environment is long - time correlated , statistical information beyond the power spectrum may be required to describe its effects .
this is the case for instance in echo protocols  @xcite .
environments with long - time memory belong to the class of adiabatic noise , for which the born - oppenheimer approximation is applicable .
we classify this part of the noise spectrum as `` adiabatic noise '' .
noise at higher frequency , around the system typical scales ( e.g. single qubit larmor frequencies ) , results in incoherent energy exchanges between the quantum device and the environment . in particular , high - frequency noise is responsible for spontaneous decay . for systems relevant for quantum information ,
the system - bath coupling is rather small , in addition high - frequency noise is short - time correlated .
therefore , in simplest cases , effects can be described by a born - markov master equation  @xcite .
for single qubit gates it leads to the relaxation and secular dephasing times , @xmath37 and @xmath38  @xcite .
we classify this part of the noise spectrum as `` quantum noise '' .
finally , resonances in the spectrum unveil the presence of discrete noise sources which severely effect the system performances , in particular reliability of devices .
this is the case when classical impurities are slow enough to induce a visible bistable instability in the system intrinsic frequency .
for instance , single qubit gates in the presence of random telegraph noise ( switching rate @xmath27 ) may display two effective frequencies , @xmath39 and @xmath40 , depending on the impurity state .
their visibility is measured by the ratio @xmath41 .
beatings can be observed in the `` strong coupling '' regime @xmath42  @xcite .
quantum impurities may also entangle with the device .
this additionally leads to a variety of features , like peculiar temperature dependencies of decay rates  @xcite . under these conditions ,
knowledge of the power spectrum is absolutely insufficient and in order to describe these effects the relevant system hilbert space has to be enlarged to include the responsible environmental degrees of freedom .
effects in general depend on the specific protocol and require a microscopic model of the fluctuators .
we classify this part of the noise spectrum as `` strongly - coupled noise '' .
each noise class requires a specific approximation scheme , which is not appropriate for the other classes .
the overall effect results from the interplay of the three classes of noise . in the following section we will illustrate a multi - scale theory to deal with solid state broadband noise .
we are interested to a reduced description of the @xmath43-qubit system , expressed by the reduced density matrix @xmath44 .
it is formally obtained by tracing out environmental degrees of freedom from the total density matrix @xmath45 , which depends on quantum ( q ) , adiabatic ( a ) and strongly coupled ( sc ) variables .
the elimination procedure can be conveniently performed by separating in the interaction hamiltonian , @xmath46 , various noise classes , e.g. by formally writing @xmath47 adiabatic noise is typically correlated on a time scale much longer than the inverse of the natural frequencies @xmath48 , then application of the born - oppenheimer approximation is equivalent to replace @xmath49 with a classical stochastic field @xmath50 .
this approach is valid when the contribution of adiabatic noise to spontaneous decay is negligible , a necessary condition being @xmath51 .
this condition is usually satisfied at short enough times , since @xmath52 is substantially different from zero only at frequencies @xmath53 .
this fact already suggests a route to trace - out different noise classes in the appropriate order .
the total density matrix parametrically depends on the specific realization of the slow random drives @xmath54 and may be written as @xmath55 .
the first step is to trace out quantum noise . in the simplest cases this requires solving a master equation . in a second stage
, the average over all the realizations of the stochastic processes , @xmath56 , is performed .
this leads to a reduced density matrix for the @xmath43-qubit system plus the strongly coupled degrees of freedom .
these have to be traced out in a final stage by solving the heisenberg equations of motion , or by approaches suitable to the specific microscopic hamiltonian or interaction .
for instance , the dynamics may be solved exactly for some special quantum impurity models at pure dephasing , @xmath57 , when impurities are longitudinally coupled to each qubit  @xcite .
the ordered multi - stage elimination procedure can be formally written as @xmath58 \ ,
p[\vec x(t ) ] \ ;   tr_{q } \big [ \ , w^{q , sc}\big(t | \vec x(t)\big ) \ , \big]\right\ } \ , .\ ] ] in the following section we concentrate on the elimination of adiabatic noise
. the ordered procedure will be illustrated in section 6 for a two - qubit gate .
in general , the adiabatic approximation holds true for times short enough to fulfill the necessary condition , @xmath59 . in the peculiar , pure dephasing regime , @xmath57 ,
relaxation processes are forbidden and the adiabatic approximation is exact for any @xmath60 .
in addition , the adiabatic scheme can be applied also in the presence of correlations between processes @xmath61 and @xmath62  @xcite . for the sake of clarity here
we consider adiabatic noise affecting independently each qubit .
moreover , we exclude time - dependent drives in @xmath8 and in @xmath63 . the procedure can be straightforwardly extended for instance to rabi oscillations and other ac - drives  @xcite .
suppose we are able to diagonalize , for each realization @xmath56 of the stochastic processes , the system hamiltonian , @xmath64 and denote @xmath65 an instantaneous eigenstate of ( [ eq : hamiltonian - stochastic ] ) with eigenvalue @xmath66 .
if the system is prepared in a pure state , @xmath67 , in the adiabatic approximation , each component @xmath68 evolves in time according to @xmath69 where @xmath70 is the dynamic phase .
the stray fields @xmath56 affect the time evolution of the gate fidelity via modifications of the phases @xmath71 and of the eigenstates @xmath72 .
in addition , @xmath73 reflects imperfect preparation due to the random field .
this occurs , for instance , when the system is prepared by a reset to a pure state vector @xmath67 followed by an initialization pulse whose effect depends also on the stray field at @xmath74 .
the evolution of the system conditional density matrix can be presented in a compact form @xmath75 \ ; { \mathrm e}^{- i \int_0^t d s \ , \omega_{m p}(s ) }   \ , , \end{aligned}\ ] ] where the instantaneous splittings appearing in the phase factor are @xmath76 and we have introduced the operator @xmath77 = \ , { \,| \ , { m(\vec x_t ) } \,\rangle   \,}{\ , \langle\,{m(\vec x_0 ) } \,|\ , }   \rho(0 ) { \,| \ , { p(\vec x_0 ) } \,\rangle   \ , } { \ , \langle\,{p(\vec x_t ) } \,|\,}\ ] ] which contains information about the preparation and the eigenstates errors . finally , the average of @xmath78 over the stochastic process yields the system density matrix , which can be presented as a path integral @xmath79 \ ,
p[\vec x(s ) ] \ , \rho^n\big(t|\vec x(s)\big )   \ , .\end{aligned}\ ] ] here @xmath80 $ ] contains information both on the stochastic processes and on details of the specific protocol .
it is convenient to split it as follows @xmath81 \,=\ , f[\vec x(s ) ] \ ; p[\vec x(s ) ] \ , , \ ] ] where @xmath82 $ ] is the probability of the realization @xmath83 .
the filter function @xmath84 $ ] describes the specific operation . for most of present day experiments on solid - state qubits
@xmath84=1 $ ] . for an open - loop feedback protocol , which allows initial control of some collective variable of the environment ,
say @xmath85 , we should put @xmath84 \propto \ , \delta(\vec x_0 ) = \ , \pi_i \delta(x_{i0})$ ] .
a critical issue is the identification of @xmath82 $ ] for the specific noise sources , as those displaying @xmath0 power spectrum .
if we sample the stochastic process at times @xmath86 , with @xmath87 and @xmath88 , we can identify @xmath89 = \lim_{m \to \infty }   p_{m+1}(\vec x_m , t;\dots ; \vec x_1 , t_1 ; \vec x_0 , 0 )   \ , ,   \label{sampling}\ ] ] where @xmath90 is a @xmath91 joint probability and we have used the shorthand notation @xmath92 . in the following we will propose a systematic method to select only the relevant statistical information on the stochastic process out of the full characterization included in @xmath82 $ ] .
we would like to remark , at this point , that in the adiabatic treatment the word decoherence is perhaps abused .
decoherence is ultimately due to entanglement of the qubit to a quantum environment  @xcite , whereas classical adiabatic noise produces only de - focusing . however , unless the signal can be _ totally _ re - focused , and this is impossible in practice , the behavior of @xmath93 is the same as for true decoherence  @xcite . in other words , although proper and improper mixed states at the fundamental level are well distinct concepts  @xcite , they are not distinct in the density matrix description .
we also observe that present day `` applied '' research on solid - state coherent nanodevices for quantum information focuses on the `` short - time '' dynamics , since a signal which is almost decayed is useless . in this context , methods as the adiabatic approximation are valuable even if they are not accurate at long time scales and are not valid down to zero temperature .
the longitudinal approximation consists in neglecting modifications of the eigenstates @xmath72 and preparation effects . without loss of generality
, we may assume vanishing average of the stochastic processes after the preparation pulse , @xmath94 .
the longitudinal approximation amounts to put @xmath95 where @xmath96 is an eigenstate of the system hamiltonian @xmath97 .
then @xmath98 \approx   { \,| \ , { m } \,\rangle   \ , } \rho^n_{mp}(0 ) { \ , \langle\,{p } \,|\,}$ ] is a projected element of the initial density matrix and eq.([eq : pathint - adiabatic ] ) simplifies to @xmath99 \ ,
p[\vec x(s ) ] \ ;   { \mathrm e}^{- i \int_0^t \!\!d s \ ,   \omega_{m p}(s ) }   \ , .\end{aligned}\ ] ] the significance of the longitudinal approximation is easily illustrated for a single qubit gate , when the system hamiltonian reduces to @xmath100 in this case , retaining variations of the splitting @xmath39 amounts to consider only fluctuations of the length of the vector @xmath101 .
this would be the only effect if the noise would act longitudinally , @xmath102 .
variations of the eigenstates @xmath103 originate instead from `` transverse '' variations of @xmath101 .
the longitudinal assumption has been performed in ref.@xcite to discuss the effect of gaussian adiabatic environments .
the present approach automatically provides constraints on its validity and shows that whereas errors due to transverse fluctuations are weakly dependent on time , phase errors accumulate .
therefore , transverse noise in the adiabatic approximation has possibly some effect only at very short times , but the phase damping channel eventually prevails .
these considerations strongly depend on the amplitude of the noise .
we checked analytically and with simulations that they hold true for realistic figures of noise as those measured in experiments by zorin et al .
@xcite .    in the longitudinal approximation , diagonal elements of the reduced density matrix in the eigenstate basis , @xmath104
, do not decay .
instead , the decay of the off - diagonal elements results from @xmath105 $ ] , where the complex average phase is given by @xmath106 \ , p[x(s ) ] \ ,   { \mathrm e}^{i \int_0^t \!\!ds   \ , \omega[x(s ) ] } \ , .\ ] ] we notice that the longitudinal approximation may be exact for certain protocols , which therefore are not affected by adiabatic transverse noise .
for instance , this is the case if the system is described by ( [ eq : hamiltonian - stochastic2 ] ) and it is prepared in a state @xmath107 .
if we measure @xmath108 , since @xmath109 do not depend on @xmath110 , the result is not affected by transverse fluctuations and @xmath111 \
, \exp\{{{\cal i } m \,}\ !
\phi(t)\}$ ] .
this the case of the decaying oscillations pattern measured with ramsey interference  @xcite , if the effect of imperfect @xmath112 pulses is negligible .      a standard approximation of the path integrals ( [ eq : pathint - adiabatic],[eq : pathint - adiabatic - longitudinal ] ) consists in neglecting the time dependence in the path , @xmath113 and taking @xmath114=1 $ ] . in this static path approximation ( spa )
the problem reduces to ordinary integrations with @xmath115 . in the single qubit case , for instance , eq.([eq : blur ] )
gives average phase shift @xmath116 where @xmath117 ^ 2 + ( \omega \cos \theta)^2}$ ] .
eq.([eq : blur - static ] ) describes the effect of a distribution of stray energy shifts @xmath118 and corresponds to the rigid lattice breadth contribution to inhomogeneous broadening  @xcite . in experiments with solid state devices
this approximation describes the measurement procedure consisting in signal acquisition and averaging over a large number @xmath119 of repetitions of the protocol , for an overall time @xmath120 ( which may also be minutes in actual experiments ) . due to slow fluctuations of the solid state environment calibration , the initial value , @xmath121 , fluctuates during the repetitions blurring the average signal , independently on the measurement being single - shot or not .
the probability @xmath122 describes the distribution of the random variable obtained by sampling the stochastic process @xmath25 at the initial time of each repetition , i. e. at times @xmath123 , @xmath124 .
if @xmath125 results from many independent random variables of a multimode environment , the central limit theorem applies and @xmath122 is a gaussian distribution with standard deviation @xmath126 @xmath127 with integration limits @xmath128 , and the high - frequency cut - off of the @xmath0 spectrum , @xmath129 . in the spa
the distribution standard deviation , @xmath126 , is the only adiabatic noise characteristic parameter .
if the equilibrium average of the stochastic process vanishes , eq.([eq : blur - static ] ) reduces to @xmath130 ^ 2 + ( \omega \cos \theta)^2 } } \ , .
\label{spa - oneqb}\ ] ] a convenient approximation is obtained by expanding @xmath131 to second order in @xmath125 , which leads to  @xcite @xmath132 the short - times decay of coherent oscillations qualitatively depends on the working point .
in fact , the suppression of the signal , @xmath133 $ ] , turns from a @xmath134 behavior at @xmath135 to a power law , @xmath136^{-1/4}$ ] , at @xmath137 . in these limits
eq.([ref : quadratic ] ) reproduces known results for gaussian @xmath0 environments .
in particular , at @xmath138 we obtain the short - times limit , @xmath139 , of the exact result of ref.@xcite . at @xmath140 the short and intermediate times result of ref .
@xcite is reproduced .
the fact that results of a diagrammatic approach with a quantum environment , as those of ref .
@xcite , can be reproduced and generalized already at the simple spa level makes the semi - classical approach quite promising .
it shows that , at least for not too long times ( but surely longer than times of interest for quantum state processing ) , the quantum nature of the environment may not be relevant for the class of problems which can be treated in the born - oppenheimer approximation .
notice also that the spa itself has surely a wide validity since it does not require information about the _
dynamics _ of the noise sources , provided they are slow  @xcite .
going beyond the spa amounts to sample more accurately the stochastic process @xmath83 in ( [ sampling ] ) .
the first correction to the spa is obtained by parametrizing the random process as follows @xmath141 . inserting this expression in ( [ eq : pathint - adiabatic ] ) and approximating @xmath82 $ ] in ( [ sampling ] ) with the conditional probability @xmath142 $ ] , we obtain @xmath143 the joint probability depends on the statistics of the noise sources .
therefore the first correction to the spa distinguishes discrete and gaussian processes . considering again a single qubit
, the average phase ( [ eq : blur ] ) at the working point @xmath144 for generic gaussian noise becomes @xmath145 t \over \omega } \big ] + { 1 \over 2 }   \ln \big[1
+   i   { \sigma^2 \pi(t)\,t \over 3 \ , \omega } \big ] \ , , \ ] ] where @xmath146 is a transition probability , depending on the stochastic process . for ornstein - uhlenbeck processes it reduces to the result of ref .
the first correction suggests that the spa , in principle valid for @xmath147 , may have a broader validity .
this is illustrated in fig.([figure1 ] ) where the adiabatic approximation ( numerical evaluation of eq.([eq : blur ] ) ) is compared with the exact ( numerical ) dynamics in the presence of @xmath0 noise and with the analytic forms resulting from the spa and its first correction . for @xmath0 noise due to
a set of bistable impurities the spa valid also for @xmath148 , if @xmath149 .
of course the adiabatic approximation is tenable if @xmath150 .    by sampling more accurately the adiabatic process
@xmath56 it is possible to selectively include the statistical information needed for the specific measurement process .
for instance , echo protocols are able to partly re - focus the signal , in other words de - focusing described by the spa is almost canceled by the echo pulse sequence .
the decay of the echo signal is due to the un - canceled dynamics of the low - frequency fluctuations and to quantum noise .
the leading effect of adiabatic noise can be estimated by a proper parametrization of @xmath56 , similar to the one considered in the present paragraph  @xcite .
in the present section we apply the multi - stage approach to a universal two - qubit gate based on a fixed coupling scheme .
capacitive and inductive fixed couplings  @xcite have been used to demonstrate two - qubit logic gates in different jj implementations  @xcite .
entanglement is generated by tuning single - qubit energy spacing to achieve mutual resonance . to this end , during the gate operation at least one qubit has to be moved away from the working point of minimal sensitivity to parameters variations , the `` optimal point '' .
this has so far represented the main drawback of the fixed - coupling scheme for josephson implementations , with the exception of phase qubits  @xcite .
more recent proposals have attempted to solve this problem by introducing tunable coupling schemes  @xcite .
most of them rely on additional circuit elements and gain their tunability from ac - driving  @xcite or from `` adiabatic '' couplers  @xcite .
some of these coupling schemes have been tested in experiments and are potentially scalable  @xcite .
none of these implementations is however totally immune from imperfections .
in general , introducing additional on - chip circuits elements opens new ports to noise . the possibility to employ
`` minimal '' fixed coupling schemes has been recently reconsidered in ref .
@xcite , pointing out the possibility to single out `` optimal coupling '' conditions which ensure reasonable protection from @xmath0 fluctuations .    in order to implement a @xmath151 gate in a fixed coupling scheme
we need two resonant qubits with a transverse coupling , as described by the hamiltonian @xmath152 where @xmath153 is the coupling strength , and @xmath154 the pseudo - spin operators whose eigenstates @xmath155 ( eigenvalues @xmath156 ) are the computational states of qubit @xmath12 .
eigenvalues and eigenvectors of @xmath157 are given by -6pt @xmath158 where @xmath159 , @xmath160 with @xmath161 and we have used the shorthand notation @xmath162 .
the level structure of the coupled system is schematically illustrated in fig.([splittings ] ) . as a result of the diagonal block structure of the hamiltonian ( [ eq : hamiltonian ] ) in the computational space ,
the two - qubit hilbert space is factorized in two subspaces spanned by pairs of eigenvectors .
a system prepared in @xmath163 , freely evolving for a time @xmath164 , yields the entangled state @xmath165 , corresponding to a @xmath151 operation .
the pair of states @xmath166 and @xmath167 span the subspace where the @xmath151 gate is realized , which we name swap - subspace .
the subspace spanned by the pair of states @xmath168 and @xmath169 is termed z - subspace .
eq.([eq : hamiltonian ] ) models , for instance , capacitive coupled charge qubits or inductively coupled flux qubits  @xcite .
typically , these system can be driven by pulses acting along the transverse direction , @xmath9 . the same port introduces noise into the system  @xcite .
this is the case of the charge - phase two - port design of ref .
@xcite . ] . here
we consider noise sources acting transversely with respect to each qubit , i. e. the interaction hamiltonian takes the form @xmath170 each noise component @xmath16 has broadband spectrum @xmath20 , @xmath21 followed by a white flank at frequencies @xmath171 . correlated noise sources acting on both qubits have been addressed in @xcite .
if the system is initialized in the swap subspace , for instance in the state @xmath172 , bitwise readout gives the qubit 1 switching probability @xmath173 , i.e. the probability that it will pass to the state @xmath174 , and the probability @xmath175 of finding the qubit 2 in the initial state @xmath174 .
cyclic anti - correlation of the probabilities signals the formation of the entangled state during the @xmath151 operation . in terms of the two qubit
reduced density matrix in the eigenstate basis the switching probabilities read @xmath176 + \rho_{00}(t ) + \left [ \ , \rho_{33}(t ) - \rho_{00}(t ) \right ] \sin^2 \frac{\varphi}{2 }   \nonumber \\ & & + { \rm re}[\rho_{12}(t ) ] +   { \rm re } [ \rho_{03}(t ) ] \sin \varphi \label{eq : psw1 } \\
\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!p^{(2)}(t )    = \langle - | { \rm tr}_1
\rho(t ) | - \rangle & = & \frac{1}{2 } \left [ \ , \rho_{11}(t ) + \rho_{22}(t ) \right ]   +   \rho_{00}(t )    +   \left [ \ , \rho_{33}(t ) - \rho_{00}(t ) \right ] \sin^2 \frac{\varphi}{2 }   \nonumber \\   & & -   { \rm re}[\rho_{12}(t ) ] + { \rm re}[\rho_{03}(t ) ] \sin \varphi .
\label{eq : psw2}\end{aligned}\ ] ] the matrix elements entering the above probabilities can be evaluated in the multi - stage approach .
we split the interaction as in eq.([int - plit ] ) , @xmath177 .
low - frequency noise is treated in the adiabatic and longitudinal approximation .
in addition we limit the analysis to the spa .
we denote with @xmath178 the eigenvalues of @xmath179 with @xmath180 .    _
first stage : elimination of quantum noise _ quantum noise is traced out by solving the born - markov master equation for the reduced density matrix  @xcite with eigenvalues parametrically dependent on the random fields @xmath181 . in the system eigenstate basis   with @xmath180 ] and performing the secular approximation ( to be self - consistently checked ) it takes the standard form  @xcite : @xmath183 where @xmath184 .
the rates @xmath185 , @xmath186 depend on the real parts of the lesser and greater green s functions , describing emission ( absorption ) rates to ( from ) the quantum reservoirs  @xcite .
in addition , for white quantum noise energy shifts are vanishing and do not appear in eq .
( [ eq : coherences ] ) .
because of the symmetry of ( [ eq : hamiltonian ] ) , dissipative transitions inside each -  swap or z  - subspace are forbidden .
the allowed inelastic energy exchange processes are evidenced in fig.([splittings ] ) .
the only independent emission rates are @xmath187 , @xmath188 .
they read @xmath189 \\
\gamma_{20 } = \frac{1}{8 } \ , ( 1-\sin \varphi ) \ , [ c_{x_1}(\omega_{20 } )   +   c_{x_2}(\omega_{20 } ) ] \ , , \end{array } \label{eq : ratesgeneral}\end{aligned}\ ] ] where the absorption rates , @xmath190 , are related to the spectrum of quantum noise , @xmath191 .
emission rates have the same form with @xmath192 replaced by @xmath193 . for the considered initial condition , @xmath194 , the only non vanishing elements of the reduced density matrix in the eigenbasis of @xmath195 are the populations and the swap coherence , @xmath196 . for independent quantum noise sources acting on the two qubits , the swap decay rate reads @xmath197 = \frac{1}{2 } \ , [ \gamma_1 + \gamma_2 ] $ ] , where @xmath198 , @xmath199 , are relaxation rates of the swap levels . from eq.([eq :
coherences ] ) we get @xmath200 the secular approximation is valid provided that @xmath201 .
this condition is fulfilled , for instance , for white noise levels extrapolated from single charge - phase qubit experiments , as it can be evinced from the slow decay due to quantum noise reported in fig.([figure2 ] ) , gray line .
equations ( [ eq : populations ] ) for the populations do not decouple even in the secular limit .
general solutions are quite lengthy . here
we report the approximate forms in the small temperature limit with respect to the uncoupled qubits splittings , @xmath202 . in this case level @xmath203
remains unpopulated and @xmath204   \quad , \quad   \rho_{11}(t )   \approx   \frac{1}{2 } \
, e^{- \gamma_{10 } t } \quad , \quad \rho_{22}(t ) \approx \frac{1}{2 } \
, e^{- \gamma_{20 } t } \ , .
\label{populations}\end{aligned}\ ] ]    _ second stage : elimination of adiabatic noise _ we now consider the effect of low - frequencies in the adiabatic , longitudinal and static path approximations .
populations are unaffected by adiabatic noise , whereas the swap coherence , @xmath205 , has to be averaged over @xmath206 . here
we disregard the negligible @xmath206 dependence of the rates @xmath207 via @xmath208 . under these simplifying assumption ,
we are left with the following average @xmath209 with @xmath210}$ ] .
the swap splitting can be estimated by treating in perturbation theory @xmath211 , with @xmath180 , with respect to @xmath157 .
this leads to @xmath212 the average in eq .
( [ eq : staticint ] ) can be evaluated in analytic form and gives  @xcite @xmath213   \label{eq : spa - swap}\ ] ] where @xmath214 , and @xmath215 $ ] is the k - bessel function of order zero  @xcite .
we considered equal standard deviations for noise sources acting on both qubits , @xmath216 . inserting ( [ populations ] ) and ( [ eq : spa - swap ] ) in ( [ eq : psw1 ] ) and ( [ eq : psw2 ] ) we obtain the switching probabilities in the multi - stage approach .
out of phase oscillations signals two - qubit states anti - correlations and follows from @xmath217 $ ] , with @xmath218 +   \cos^2\left ( \frac{\varphi}{2}\right ) \approx   - \frac{1}{4 } \
, \cos \varphi \ , \left [ e^{- \gamma_{10}t } + e^{- \gamma_{20}t } \right ] +   g^2/4 $ ] .
the efficiency of the gate results from the interplay of quantum and adiabatic noise .
high - frequency noise levels expected from single qubit experiments  @xcite weakly affect the switching probability , whose decay is mainly due to low - frequency noise , fig .
( [ figure2 ] ) .
remarkably , considerable recovery of short - times oscillation amplitude may be achieved by an optimal choice of the coupling strength , @xmath219  @xcite .
this regime is illustrated in fig.([figure3 ] ) .
in this article we presented a road - map to treat broadband noise typical of solid state nanodevices .
the introduced multi - stage approach allows to obtain reasonable approximations by systematically including only the relevant information on the complex environment , out of the huge parametrization which would be required for a microscopic description .
since the environment is in general long - time correlated , the required information depends on the specific protocol .
the predictions obtained with the present approach are in agreement with observations in various single qubit jj implementations and in different protocols  @xcite .
we extended the procedure to deal with complex solid - state architectures .
this is a required step in order to predict efficiency and possibly appropriately design of nanodevices for quantum information processing . both because of the complexity of architectures and of the unavoidable broadband nature of solid state noise , theoretical tools allowing systematic and controlled approximations are particularly valuable .    as an illustrative case
, we performed a simplified analysis of the effects and interplay of low and high frequency noise components in a two - qubit gate in a fixed coupling scheme .
our results points out that efficient operations in the solid state require an accurate preliminary characterization of the noise spectral characteristics and tuning appropriately the device working point .
we acknowledge support from the eu - eurosqip ( ist-3 - 015708-ip ) .
99 m. nielsen , i. chuang , `` quantum computation and quantum information '' , cambridge univ .
press , 2005 .
t. hyashi _
et al . _ ,
lett . * 91 * , 226804 ( 2003 ) ; j. r. petta _ et al .
lett . * 93 * , 186802 ( 2004 ) ; j.r .
et al . _ ,
science * 309 * , 2180 ( 2005 ) ; j. gorman , d.g .
hasko , and d.a .
williams , phys .
lett . * 95 * , 090502 ( 2005 ) ; f. h. l. koppens _ et al .
_ , nature * 442 * , 766 ( 2006 ) .
y. nakamura _
et al . _ ,
nature * 398 * , 786 ( 1999 ) ; y. yu _ et al .
_ , science * 296 * , 889 ( 2002 ) ; j.m . martinis _
et al . _ ,
lett . * 89 * , 117901 ( 2002 ) ; i. chiorescu _
et al . _ ,
science , * 299 * , 1869 , ( 2003 ) ; t. yamamoto _ et al .
_ , nature * 425 * , 941 ( 2003 ) ; s. saito _ et al .
* 93 * , 037001 ( 2004 ) ; j. johansson
_ et al . _ ,
ibid . * 96 * , 127006 ( 2006 ) ; f. deppe _
et al . _ , nat .
* 4 * , 686 ( 2008 ) ; d. vion _
et al . _ ,
science * 296 * , 886 ( 2002 ) .
j. koch et al .
a * 76 * , 042319 ( 2007 ) ; j. a. schreier _
et al . _ ,
b * 77 * , 180502(r ) ( 2008 ) . d. i. schuster , _ et al .
_ nature * 431 * , 162 ( 2004 ) ; d. i. schuster _
nature * 445 * , 515 ( 2007 ) ; m. gppl _
nature * 454 * , 315 ( 2008 ) ; lev s. bishop _ et al .
_ nature physics * 5 * , 105 ( 2008 ) .
y. nakamura _ et .
al _ , phys .
lett . * 88 * , 047901 ( 2002 ) .
g. falci , e. paladino , r. fazio , in _ quantum phenomena of mesoscopic systems _ ,
b. l. altshuler and v. tognetti eds .
, proc . of the international school `` enrico fermi '' , varenna 2002 , ios press ( 2003 ) , cond - mat/0312550 .
lett . * 94 * , 167002 ( 2005 ) . c. p. slichter _
principles of magnetic resonance _ , springer - verlag , berlin ( 1996 ) s. kafanof _ et al .
b , * 78 * , 125411 ( 2008 ) .
o. astafiev _
lett . * 93 * , 267007 ( 2004 ) .
g. ithier _
et al . _ ,
b * 72 * , 134519 ( 2005 ) . r.w .
et al . _ ,  phys .
lett . * 93 * , 077003 ( 2004 ) ; k.b .
cooper _ et al .
_ , ibid . * 93 * , 180401 ( 2004 ) .
j. eroms _ et al .
. lett . * 89 * , 122516 ( 2006 ) . c. h. van der wal , _ et al .
_ , science 290 , 773777 ( 2000 ) . t. duty _ et al .
_ , phys . rev .
b * 69 * , 140503(r ) ( 2004 ) ; y. m. galperin _
96 * , 097009 ( 2006 ) ; j bergli , y m galperin and b l altshuler , new j. phys . *
11 * , 025002 ( 2009 ) .
f. chiarello private communication 2008 .
experiment based on the setup of s.poletto _ et al .
_ new j. phys . * 11 * , 013009 ( 2008 ) .
r. fazio , g. m. palma , and j. siewert phys .
83 * , 5385 ( 1999 ) . , rev .
* 60 * , 537 ( 1988 ) .
j. claudon _
et al . _ ,
b * 76 * , 024508 ( 2007 ) ; l. tian and r.w .
simmonds , phys .
99 * , 137002 ( 2007 ) ; f. deppe _
b * 76 * , 214503 ( 2007 ) ; z. kim , _
b * 78 * , 144506 ( 2008 ) ; a. lupascu _
_ arxiv:0810.0590 .
a. m. zagoskin _
lett . * 97 * , 077001 ( 2006 ) .
m. neeley , _ et al .
_ nature physics * 4 * , 523 ( 2008 ) .
e. paladino _
b * 77 * , 041303(r)(2008 )
. s. ashab _
_ , new j. phys .
* 8 * , 103 ( 2006 ) ; n. p oxtoby _ et al .
_ new j. phys . * 11 * , 063028 ( 2009 ) . c. cohen - tannoudji , j. dupont - roc and g. grynberg _ atom - photon interactions _ , wiley - interscience ( 1993 ) .
e. paladino _
et al . _ ,
* 88 * , 228304 ( 2002 ) .
correlated noise and cross - talk effects have been addressed in a. darrigo _ et al .
njp * 10 * , 115006 ( 2008 ) . g. falci _
_ in preparation .
w. zurek , physics today * 44 * , 36 ( 1991 ) .
g. m. palma , k.a .
suominen , a. k. ekert , proc .
london a , * 452 * , 567 ( 1996 ) .
d. g. cory _
* 81 * , 2152 ( 1998 ) . c. isham , _ lectures on quantum theory : mathematical and structural foundations _ , imperial college press , london ( 1995 ) .
y. makhlin , a. shnirman , phys .
lett . * 92 * , 178301 ( 2004 ) .
a.b . zorin _
et al . _ ,
b * 53 * , 13682 ( 1996 ) .
( [ ref : quadratic ] ) is also valid for ornstein - uhlenbeck processes , see @xcite . for random telegraph noise
the discrete nature of the process modifies this result see e. paladino _
state phys . ,
* 43 * , 747 ( 2003 ) k. rabenstein , v.a .
sverdlov , d.v .
averin jetp letters * 79 * , 646 ( 2004 ) .
et al . _ ,
nature * 398 * , 305 ( 1999 ) ; j.  q. you _ et al .
lett . * 89 * , 197902 ( 2002 ) . .
a. pashkin _
et al . _ ,
nature * 421 * , 823 ( 2003 ) ; a.  j. berkley _
et al . _ , science * 300 * , 1548 ( 2003 ) ; t.  yamamoto _ et al .
_ , nature * 425 * , 941 ( 2003 ) ; r.  mcdermott _ et al .
_ , science * 307 * , 1299 ( 2005 ) ; j.  b. majer _
lett . * 94 * , 090501 ( 2005 ) ; m. steffen _
et al . _ ,
science * 313 * , 1423 ( 2006 ) ; j.h .
et al . _ ,
nature * 447 * , 836 ( 2007 ) .
m. a. sillanp , j. i. park and r. w. simmonds nature * 449 * , 438 ( 2007 ) .
averin , c. bruder , phys .
* 91 * , 057003 ( 2003 ) ; a. blais _
et al . _ ,
ibid . * 90 * , 127901 ( 2003 ) ; f. plastina , g. falci , phys .
b * 67 * , 224514 ( 2003 ) ; b.plourde _ et al . _ ,
ibid . * 70 * , 140501(r ) ( 2004 ) ; a.o .
niskanen , y. nakamura , j.s .
tsai , ibid .
* 73 * , 094506 ( 2006 ) ; p. bertet , c.  j. harmans , j.e .
mooij , ibid .
* 73 * , 064512 ( 2006 ) ; y - d wang , a. kemp , k. semba , ibid .
* 79 * , 024502 ( 2009 ) . c.  rigetti , a.  blais , and m.  devoret , phys .
* 94 * , 240502 ( 2005 ) ; yu - xi liu _ et al .
_ , ibid . * 96 * , 067003 ( 2006 ) ; g. s. paraoanu , phys . rev .
b 74 , 140504(r ) ( 2006 ) .
t.  v. filippov _
et al . _ , ieee trans .
* 13 * , 1005 ( 2003 ) ; j.  lantz _ et al .
b * 70 * , 140507(r ) ( 2004 ) ; a.  maassen van den brink _
et al . _ , new j. phys . * 7 * , 230 ( 2005 ) .
a. izmalkov _
* 93 * , 037003 ( 2004 ) ; m. grajcar _
_ , prb * 72 * , 02503 ( 2005 ) ; t. hime _
et al . _ ,
science * 314 * , 1427 ( 2006 ) ; a. o. niskanen _
_ , science * 316 * 723 ( 2007 ) ; a. o. niskanen _
et al . _ ,
nature * 447 * , 386 ( 2007 ) ; j. majer _
_ , ibid . * 449 * , 443 ( 2007 ) ; s. h. w. van der ploeg _
et al . _ ,
* 98 * , 057004 ( 2007 ) ; a. fay _
* 100 * , 187003 ( 2008 ) ; a. o. niskanen _
et al . _ ,
b * 77 * , 064505 ( 2008 ) ; k. harrabi _
_ , phys . rev .
b 79 , 020507 ( 2009 ) .
e. paladino , a. mastellone , a. darrigo , g. falci , arxiv:0906.3115 u. weiss , `` quantum dissipative systems '' , third edition , world scientific 2008 . m. abramowitz , i. a. stegun `` handbook of mathematical functions '' , dover ( 1965 ) .
e. paladino _ at al .
_ physica e in press ( 2009 ) .