recent progress in the physics of bose - einstein condensation ( bec ) of atomic gases ( see refs .
@xcite and references therein ) stimulated a further development of theoretical methods used previously for the description of superfluidity , superconductivity , and collective effects in heavy atoms and nuclei ( see refs.@xcite ) .
these methods treat condensates as mean fields , while interactions between particles , according to the bogolubov theory , lead to formation of fluctuations ( see refs.@xcite ) . in the case of a positive scattering length ,
the condensate ( mean field ) is stable , while fluctuations can be treated by using a bogolubov transformation to quasiparticle states .
recent experiments @xcite drew attention to the physics of condensates with a negative scattering length . in this case , a free condensate ( with a homogeneous mean field ) is unstable , collapsing within a finite time @xcite , while a trapped condensate is stable at the low occupations of the experiment @xcite . in the jila experiments on rb @xcite , the effect of feshbach resonance ( see ref .
@xcite ) was used for changing the sign of the scattering length .
this effect allows the formation of a condensate with a positive scattering length , followed by a study of the dynamics of collapse after applying a sudden change of sign of the scattering length .
a formal extension of the bogolubov transformation to the case of negative scattering lengths would lead to complex energies of the quasiparticles , but actually such a transformation does not exist in this case , as was proven by bogolubov ( see ref .
nevertheless , the development of fluctuations can substantially affect the dynamics of an unstable condensate ( see the theoretical analysis @xcite of the jila experiments @xcite ) .
the present paper describes the dynamics of fluctuations by using the parametric approximation employed in ref .
@xcite for a model problem of a two - mode atom - molecule condensate , and in ref.@xcite for a description of quantum effects on curve crossing .
this method is similar to the parametric approximation used in quantum optics ( see ref .
@xcite ) , and in the description of interaction of electromagnetic and matter waves in refs .  @xcite .
it is also closely related to the linearization approximation in quantum soliton theory ( see ref .
@xcite and references therein ) and to the theory of quantum integrals of motion ( see refs .
@xcite and references therein ) .
the fluctuations are treated here by a linearization of the exact quantum equations of motion for the field operators . in the case of a stable bec ,
the resulting linear equations describe the motion of quantum oscillators , corresponding to the bogolubov quasiparticles .
however , in the case of unstable condensates , these equations describe inverted oscillators , leading to a rapid growth of the fluctuations .
this growth remains super - exponential as long as the fluctuations do not reach a macroscopic occupation .
this dynamics offers a more precise picture compared to the exponential growth obtained from a liapunov analysis of the classical mean field .
initially , the dynamics can be described as a quantum process of spontaneous decay of the condensate into the fluctuation vacuum ( see refs.@xcite ) , accompanied by the formation of entangled pairs ( see also refs .
@xcite ) . as in the case of fluctuations in a stable condensate
( see ref .
@xcite ) , the pairs are formed in two - mode squeezed states .
the formally calculated energy of a bogolubov quasiparticle becomes complex also when the energy of the excitation is lower than the energy of a particle in the condensate state .
this case can take place if the condensate itself occupies an excited translational state , as realized recently in the nist experiments on four - wave mixing of atomic condensate waves @xcite .
the loss of condensate atoms in such a situation was studied in ref .
@xcite , using the method of complex scattering lengths , without taking into account second - quantized properties of the atomic scattering .
the results of the present theory , which take into account these properties , are in agreement with ref .
@xcite at short collision times [ see eq.([ndotstat ] ) below ] .
however , at long collision times the effects of bose enhancement lead to a substantial increase of the losses [ see eq .
( [ ndotenh ] ) below ] .
notice that the necessity of accounting for the quantum effects was mentioned in ref .
@xcite itself .
these effects should appear under the conditions of the nist experiments @xcite .
the present method can describe only qualitatively this case of large losses , comparable in magnitude to the initial occupation .
this case may be treated by the incorporation of quantum effects into the slowly - varying - envelope approximation developed in ref .
@xcite .
quantum effects on the loss incurred by collisions of two condensate wavepackets were studied in ref .
@xcite , using an approach similar to the present one . however , the results of the present work disagree with those of ref .
the losses quoted in ref .
@xcite , which have been evaluated only for short collision times , differ also from the results of ref.@xcite .
these disagreements are apparently due to some unjustified neglections in ref .  @xcite ( see discussion in sec .  [ gauss ] ) .    as in the case of fluctuations in a ground - state condensate ,
the scattered particles are formed here in entangled ( squeezed ) states .
therefore the processes considered here may open up possibilities for the formation of an entangled gas remained after the decay of a condensate , or after the collision of wavepackets .
the entangled gas can have many application ( see ref .
@xcite ) .
this paper is organized as follows .
section [ ground ] considers the dynamics of instabilities in a ground - state condensate with a negative scattering length .
the equations of motion for the condensate and the fluctuations are derived in sec .
[ condsep ] by using a projection onto a coherent state .
the equations for the fluctuations are solved in sec .
[ stabunstmod ] . in the dependence on the excitation energy ,
the fluctuations are divided into two classes , of stable and unstable modes , which are respectively associated with regular and inverted oscillators in sec.[reginvosc ] .
the dynamics and statistical properties of the instabilities are considered in sec .
[ instdyn ] , and sec.[conddepl ] describes the dynamics of condensate depletion .
the collision of two matter waves is studied in sec.[twocond ] .
general equations for the condensate and the fluctuations are derived in sec .
[ tcsep ] .
section [ loss ] considers elastic losses in collisions of two matter waves .
the collisions of rectangular and gaussian wavepackets are treated in detail in secs .
[ rect ] and [ gauss ] , respectively .
a system of units in which planck s constant is @xmath0 is used below .
consider a homogeneous gas of bosons of mass @xmath1 , in a finite quantization box of volume @xmath2 , described by the following hamiltonian in the momentum representation ( see refs .
@xcite ) , @xmath3 here @xmath4 is the annihilation operator for bosons in the momentum @xmath5-mode , @xmath6 is the kinetic energy , and @xmath7 describes the elastic scattering . in the model case of a zero - range contact potential
, @xmath8 is a constant , and can be expressed in terms of the elastic scattering length @xmath9 as @xmath10 ( a modification of the contact - potential model was recently proposed in ref .
@xcite which allows to avoid ultraviolet divergences ) .
the hamiltonian ( [ h ] ) leads to the following equations of motion for the annihilation operators : @xmath11 following a general method @xcite let us separate the condensate part by a projection onto a coherent state @xmath12 defined so that @xmath13 describes the ground - state condensate occupying the momentum state @xmath14 .
the annihilation operator can be written as @xmath15 where @xmath16 , the operators @xmath17 represent the non - condensate part ( fluctuations ) having the zero expectation value , @xmath18 and @xmath19 is the conserved total number of bosons .
the exponential factor in eq .
( [ psiphi ] ) shifts the zero energy to the energy of non - excited condensate .
substitution of eq .
( [ psiphi ] ) into eq .
( [ pside ] ) , combined with the property ( [ csiev ] ) , produces equations of motion for @xmath20 and @xmath21 . assuming a small occupation of the fluctuations ,
let us keep in these equations only terms of lowest non - vanishing order in @xmath22 and @xmath23 , obtaining the following equations : @xmath24 where @xmath25 and @xmath26 is the condensate density .
this approximate method , that leads to linear operator equations for the fluctuations , is similar to the parametric approximation used in quantum optics ( see ref .
it has been applied to the study of a molecular bec dissociation in refs .
@xcite ) , and is valid as long as the depletion of the condensate is negligibly small .
an applicability criterion for this method is presented in sec .
[ conddepl ] below [ see eq .
( [ valran ] ) ] .
the present approximation is gapless [ see eq .
( [ lambda ] ) below ] and , therefore , non - conserving ( see ref .
@xcite ) , although the average number of the atoms is conserved [ see eq .
( [ nbos ] ) above ] .
equation ( [ csip ] ) combined with the corresponding equation for @xmath27 provide a set of two linear operator equations ( similar to the set considered in refs .
this set can be solved exactly , using common methods reducing it to an eigenproblem , by assuming the time - independence of @xmath28 and @xmath29 .
this assumption is corroborated by eq .
( [ phi ] ) , which shows that the variation of @xmath30 is of the second order in @xmath21 and therefore can be neglected .
indeed , the validity of the parametric approximation requires a small variation of the condensate density , and the dominant part of the condensate phase is accumulated in the exponential factor in eq .
( [ psiphi ] ) .    as in the case considered in ref .
@xcite , the form of the solution depends on the ratio of the detuning @xmath31 to the coupling strength @xmath29 .
all momentum modes observing @xmath32 are unstable , and the corresponding solutions attain the aperiodic form @xmath33 where @xmath34 the condition ( [ inst ] ) can be obeyed only when @xmath35 , i. e. , with a negative scattering length , and for excitations with an energy @xmath36 .
otherwise , whenever @xmath37 , the modes are stable and the solutions are oscillatory functions of the form @xmath38    in a case of intervening mode ( @xmath39 ) the solution is a linear function of time @xmath40 this case can take place for the spurious ( goldstowne ) mode with a vanishing excitation energy ( @xmath41 , see refs.@xcite ) , or at @xmath42 .
the stable modes can be treated by the bogolubov transformation .
it is therefore instructive to establish a connection between the present and the bogolubov methods .
the equations of motion ( [ csip ] ) can be derived from the following hamiltonian constituent involving two counter - propagating modes , @xmath43 the transformation @xmath44 introduces a new pair of field operators , @xmath45 and @xmath46 , whose equations of motion are decoupled by separating the hamiltonian into @xmath47 , where @xmath48 ( hereafter in this subsection the arguments @xmath5 and @xmath49 are dropped out , for simplicity s sake . )
a system described by the hamiltonian ( [ hj ] ) was solved in ref .
@xcite , and the solutions have a form similar to eq .
( [ solap ] ) , or eq .  ( [ solp ] ) , according to the ratio of @xmath50 to @xmath51 .
let us introduce canonically conjugate operators for the coordinates @xmath52 and the momenta @xmath53 , such that @xmath54 , writing out the hamiltonian ( [ hj ] ) as @xmath55 if @xmath56 , this hamiltonian describes a harmonic oscillator and can be transformed into the standard form representing free bogolubov quasiparticles [ the total transformation , incorporating ( [ csi12 ] ) , being indeed the bogolubov transformation ] .
the motion of the harmonic oscillator is periodic , in agreement with the periodic behavior of the solution ( [ solp ] ) .
a different situation takes place if the condition ( [ inst ] ) holds . in this case one of the coefficients ( preceding @xmath57 or @xmath58 ) is negative , and the hamiltonian ( [ hpq ] ) represents an inverted oscillator .
it can not be transformed to a form representing free quasiparticles , and therefore a bogolubov transformation can not exist , in agreement with ref .
the infinite aperiodic motion of the inverted oscillator is exhibited in eq .
( [ solap ] ) and in the instability of the mode .    in a case of intervening mode ( @xmath39 ) eq .
( [ lambda ] ) gives @xmath59 , one of the coefficients ( preceding @xmath58 or @xmath57 ) vanishes , and the hamiltonian ( [ hpq ] ) represents a motion without restoring forces or an oscillator with infinitely large mass , respectively ( see ref .
such systems demonstrate linear time - dependence of @xmath53 or @xmath52 , respectively , in agreement with the solution ( [ soli ] ) .
the bogolubov transformation does not exist in this case as well . in the case of the goldstowne mode ( @xmath60 )
the hamiltonian ( [ hp ] ) has already the form ( [ hj ] ) , and the transformation ( [ csi12 ] ) is not necessary .
this mode brings only one term @xmath61 into the total hamiltonian , in agreement with ref .
@xcite ( up to the exchange of the coordinate and the momentum ) .
consider the dynamics of fluctuations in an unstable bec , starting from the state of a pure condensate @xmath62 this problem is not only of a speculative interest ; it actually describes the recent jila experiment @xcite , in which the elastic scattering length was varied very fast from the zero value ( of the ideal condensate ) to a negative value ( of an unstable condensate ) .
given the initial conditions ( [ purecond ] ) , eq .  ( [ solap ] ) leads to the following dynamics for the occupation of the unstable modes , @xmath63 and for the two - particle correlations , @xmath64 it is not surprising that this dynamics is similar to the dynamics of spontaneous dissociation of a molecular condensate , studied in ref .
@xcite , as the non - condensate modes can be described in terms of the fields @xmath65 [ see eq .
( [ csi12 ] ) ] by the hamiltonian ( [ hj ] ) , which is equivalent to the one considered in ref .  @xcite . by analogy to the results of ref .
@xcite , the non - condensate atoms are formed in squeezed states involving atom pairs of opposite momenta .
these states are similar to the two - mode squeezed states of fluctuations in a stable bec studied in ref .  @xcite .
the states are perfectly squeezed and have zero variance in the relative number of particles with momenta @xmath5 and @xmath66 , being optimally entangled in the entropic sense .
the uncertainty of the generalized coordinates @xmath67 where the operators @xmath68 ( with @xmath69 ) are defined by eq .
( [ csi12 ] ) , attains the maximal and minimal values @xmath70 at the two orthogonal values of the respective phase angles @xmath71 , @xmath72    a common method for studying bec stability is the liapunov analysis ( see e.  g.  ref .
this method produces equations for small perturbations @xmath73 to the mean field .
these equations are equivalent to eq .
( [ csip ] ) , in which the operators @xmath22 are replaced by @xmath74-number functions @xmath73 . in the case
studied here one obtains for the liapunov exponents the expression @xmath75 given by eq.([lambda ] ) .
this type of analysis attributes an exponential gain to the small perturbation , @xmath76 , when the condition ( [ inst ] ) is obeyed .
this result is in qualitative agreement with the quantum analysis used here , which attributes an even faster ( super - exponential ) gain to the fluctuations [ see eq .
( [ den ] ) ] , as @xmath77 the dynamics of the fluctuations has been studied in ref.@xcite .
it is shown there that the quantum ( squeezed ) state always grows faster than a classical ( coherent ) state .
therefore the liapunov analysis concerning classical fields underestimates the fluctuation growth , becoming correct only when @xmath78 , where the fluctuations are large enough to be considered as classical ones .    a formal application of the bogolubov transformation at @xmath79 ( although incorrect , as had been observed by bogolubov himself ;
see ref.@xcite ) produces quasiparticles with imaginary energies @xmath80 , which correspond to exponential growth .
however , the transformed field operators do not obey the proper bose commutation rules .
the homogeneous boson gas with a negative scattering length considered above always has unstable modes due to the continuous nature of the excitation spectrum .
the situation changes in the case of a trapped gas having a discrete excitation spectrum .
following ref .
@xcite , consider @xmath81 as a function of the condensate density @xmath82 for a fixed value of @xmath83 chosen so that @xmath84 , where @xmath85 is the energy of an excited state of the trapped condensate . at low condensate density @xmath86
is imaginary , corresponding to the well - known stability of a trapped condensate with a negative scattering length at low occupations ( see ref.@xcite ) . at @xmath87
the value of @xmath88 becomes real , and therefore describes an instability of the condensate ( leading to its collapse ) .
although the present approach does not take into account the inhomogeneity of the condensate and the variation of its size , the dependence of @xmath89 on @xmath82 [ see eq .
( [ lambda ] ) ] is in qualitative agreement with the one obtained from the numerical calculations of ref .
@xcite for the monopole mode .
let us consider now the dynamics of condensate depletion . in the limit @xmath90 , eqs .
( [ den ] ) , ( [ cor ] ) , and ( [ phi ] ) lead to the following equation for the condensate density , @xmath91 where @xmath92 ) and stable ( @xmath93 ) modes , respectively , @xmath94 is the boundary between the stable and unstable modes , and the contact potential ( [ ua ] ) is used hereafter in this section .
the statistical weight of intervening mode with @xmath39 is negligible compared to the continuums of the stable and unstable modes . unlike a case of stable condensate , where the unrestricted behavior of the goldstowne mode requires a special treatment ( see ref .
@xcite ) , in the present case the effect of linearly - growing intervening modes is negligible compared to the effect of exponentially - growing unstable modes .
substitution of @xmath95 allows us to express @xmath96 in terms of the weber functions @xmath97 ( see ref .
@xcite ) : @xmath98 where @xmath99 , and @xmath100 is the nonlinear interaction time ( see ref .
approximate expressions for @xmath96 , holding at small and large values of @xmath49 , have the following respective forms , @xmath101 the contribution of the stable modes can be evaluated by substituting @xmath102 , as a result of which @xmath103 its limiting values can be approximately expressed as @xmath104 finally , the rate of condensate depletion can be expressed as @xmath105 @xmath106    at the two distinguishable loss regimes . at the regime of non - stationarity
@xmath107 the depletion is mostly due to an occupation of the stable modes , which carry much more statistical weight than the unstable modes [ @xmath108 , see eqs.([iinsta ] ) and ( [ istaba ] ) ] .
the statistical weight is infinitely large , leading to a singularity in the depletion rate ( [ ndot ] ) , at @xmath109 .
this singularity is absent in eq .
( [ deltan ] ) , and can be eliminated by using of an interatomic interaction with a finite range in place of the contact potential ( [ ua ] ) . at @xmath110
the occupation of the modes with @xmath111 begins to decrease due to their oscillating behavior , leading to decrease of the depletion rate . at the regime of bose - enhancement @xmath112
the depletion is dominated by the exponentially increasing occupation of the unstable modes [ @xmath113 , see eqs.([iinsta ] ) and ( [ istaba ] ) ] .
the approximate expressions ( [ deltan ] ) are compared in fig .
[ fg ] to results of numerical integration of eq .
( [ ndot ] ) .
one should be reminded that the parametric approximation , on which all the analysis above rests , is applicable only as long as @xmath114 . in the case of one excited mode @xcite ,
this approximation is in good agreement with numerical results while @xmath115 . in the case of a dilute condensate ( @xmath116 ,
see ref .
@xcite ) , the validity range @xmath117 includes a region @xmath118 where the exponential bose - enhancement factor can reach a rather large value , of the order of @xmath119 ) predicts small depletion at times @xmath120 ( with @xmath121 , given the dilution parameter @xmath122 is small ) .
therefore , the development of instabilities in bec , including collapse , starts with a delay time @xmath123 after turning on the interatomic interactions , although the collapse itself is a nonlinear phenomenon that can not be treated by the parametric approximation . under the conditions of the @xmath124rb experiment @xcite ( where @xmath125@xmath126 and @xmath127 nm ) , the expected delay is @xmath128ms ( cf .
[ fg ] ) , in agreement with the experimentally observed delay time for the collapse . in this case
the parametric approximation is valid as long as @xmath129 , allowing for values of the exponential bose - enhancement factor in eq .
( [ deltane ] ) exceeding 100 .
it should be noted that the condensate dynamics can be sensitive to initial correlations and the spatial inhomogeneity neglected here .
the approach developed in sec .  [ ground ] above can also be applied to a collision of two free counterpropagating bec waves with momenta @xmath130 ( see ref .
@xcite ) . in this case
the annihilation operator should be represented as @xmath131 where @xmath132 here @xmath133 describe the two condensate fields with momenta @xmath130 , respectively , the coherent states @xmath134 are defined by @xmath135 and the conserved total number of atoms is @xmath136 substitution of eq .
( [ psi2p ] ) into eq.([pside ] ) gives the following equations of motion for the fluctuations , @xmath137 where @xmath138 and @xmath139 are the densities of the condensate waves .
the term @xmath140 \hat{\xi } \left ( { \bf p}-2{\bf p}_{0},t\right )   \nonumber \\ & & + \phi ^{2}_{+}\left ( t\right ) u\left ( |{\bf p}-{\bf p}_{0}|\right )   \hat{\xi } ^{\dag } \left ( -{\bf p}+2{\bf p}_{0},t\right ) \ } \nonumber \\ \times&&\exp\{2i[s_+(t)-s_-(t)]\ }   \nonumber \\
+ \{&&\phi ^{*}_{+}\left ( t\right ) \phi _ { -}\left ( t\right ) \left [   u\left ( 2{\bf p}_{0}\right ) + u\left ( |{\bf p}+{\bf p}_{0}|\right )   \right ] \hat{\xi } \left ( { \bf p}+2{\bf p}_{0},t\right )   \nonumber \\ & & + \phi ^{2}_{-}\left ( t\right ) u\left ( |{\bf p}+{\bf p}_{0}|\right )   \hat{\xi } ^{\dag } \left ( -{\bf p}-2{\bf p}_{0},t\right ) \ } \nonumber \\ \times&&\exp\{2i[s_-(t)-s_+(t)]\}\end{aligned}\ ] ] describes couplings of fluctuations with momenta different by @xmath141 . whenever @xmath142 the disparity of the energies of such fluctuations far exceeds the coupling strengths , and the term @xmath143 can be neglected in eq .
( [ csip2 ] ) .
this allows us also to neglect correlations of such fluctuations in the equations of motion for the condensate fields , which then attain the form @xmath144 equation ( [ csip2 ] ) with @xmath143 neglected has solutions of the form eqs.([solap ] ) , ( [ solp ] ) , or ( [ soli ] ) , according to the ratio of @xmath145 to @xmath146 , in which @xmath51 , @xmath50 , and @xmath75 are respectively replaced by @xmath147 , @xmath146 , and @xmath148    the condition of instability @xmath149 can be obeyed now even for positive values of the scattering length ( and hence of @xmath8 ) , as @xmath146 contains a negative term @xmath150 [ see eq .
( [ d2 ] ) ] . in the case of a contact potential ( [ ua ] )
an instability takes place when @xmath151 or when @xmath152 where @xmath153 whenever condition ( [ highp ] ) is obeyed , the momenta @xmath154 are given approximately by @xmath155 the instability interval can include the region @xmath156 .
it may appear that this fact can lead to the absurd conclusion that a ground - state condensate ( @xmath157 ) with a positive scattering length is unstable .
however , the analysis above is valid only when the condition ( [ highp ] ) is obeyed .
otherwise , the terms neglected in eqs.([csip2 ] ) and ( [ phi2 ] ) have to be taken into account , leading to stability of the condensate at @xmath157 .    as in the case of a ground - state condensate ,
the scattered atoms with momenta @xmath5 and @xmath66 are correlated , forming optimally entangled and optimally squeezed states [ see discussion following eq.([cor ] ) ] .
consider now the elastic - scattering loss of atoms from the two colliding counterpropagating waves bounded only by the normalization box .
substitution of eqs .
( [ den ] ) and ( [ cor ] ) ( with the replacement of @xmath50 , @xmath51 , and @xmath75 by @xmath146 , @xmath147 , and @xmath158 , respectively ) into eq .
( [ phi2 ] ) leads to the following equations for the two condensate densities , @xmath159 where latexmath:[\[\begin{aligned } & & i_{\text{unst}}={1\over m p_0}\int\limits^{p { } _ { + } } _ { p { } _ { -}}p^{2}{\sinh\left ( 2\lambda _ { 2}\left ( p\right ) t\right ) \over   \lambda _ { 2}\left ( p\right ) } d p \nonumber \\ & & i_{-}={1\over m p_0}\int\limits^{p { } _ { -}}_{0}p^{2}{\sin\left ( 2     p\right ) |}d p \label{iunstpm } \\ & & i_{+}={1\over m p_0}\int\limits^{\infty } _ { p { } _ { + } } p^{2}{\sin\left ( 2|\lambda _ { 2}\left ( p\right ) |t\right ) \over    the contributions of the unstable ( @xmath96 ) and stable ( @xmath161 ) modes , respectively .    the substitution of @xmath162 allows us to express @xmath96 in terms of the modified struve function @xmath163 ( see ref .
@xcite ) , as @xmath164 where @xmath99 , and here @xmath165 .
approximate expressions for @xmath96 at small and large values of @xmath166 have the respective forms @xmath167 given the condition @xmath168 , the substitution of @xmath169 produces the estimate @xmath170 the integral @xmath171 can be estimated under the same condition by retaining only the first term in the taylor series for @xmath172 .    given the opposite condition , @xmath173 , the contributions of the stable modes @xmath174 can be evaluated by substituting @xmath175 resulting in @xmath176 where @xmath177 and @xmath178 ( with @xmath179 ) . the total contribution of the stable modes , given by @xmath180 where @xmath181 is the modified bessel function ( see ref.@xcite ) , can be approximately expressed as    @xmath182    combining eqs .
( [ npm ] ) , ( [ iinst2 ] ) , and ( [ ipm ] ) , one obtains the total rate of the condensate loss @xmath183    unlike the case of ground - state condensate , characterized by a single time - scale @xmath184 , an additional time - scale , @xmath185 , appears in the present case , leading to appearance of a third regime of the condensate loss . at short times , @xmath107 ,
the loss is dominated by the occupation of the stable modes due to their higher statistical weight .    at the regime of non - stationarity @xmath168
, the wavefunction of the scattered particles does not have enough time to become stationary , and the estimate ( [ imnonst ] ) gives the loss rate @xmath186 the singularity at @xmath109 has the same nature as in the case of the ground - state condensate with a negative scattering length [ see discussion following eq.([deltan ] ) ] .
at the regime of stationarity , @xmath187 , the loss rate becomes time independent , and given by @xmath188 in full agreement with theory @xcite .    at the regime of bose - enhancement , @xmath112 , the loss
is principally determined by the occupation of unstable modes which increase exponentially due to the effect of bose enhancement . the loss rate
can then be approximately expressed as @xmath189    the approximate expressions ( [ npmdotn ] ) , ( [ ndotstat ] ) , and ( [ ndotenh ] ) are compared in fig .
[ fc ] to results of numerical integration in eqs .
( [ iunstpm ] ) .
the contribution of stable modes decreases at long times as @xmath49 increases [ see eq.([ipmlt ] ) ] . in order to explain the decrease ,
let us consider the asymptote of @xmath161 using , in the spirit of scattering theory @xcite , the substitution @xmath190 in the absence of many - body effects ( @xmath191 ) , one obtains @xmath192 , expressing the regular dispersion law for free particles . also , @xmath193 , @xmath194 , and hence @xmath195 , which also hold for small @xmath49 [ see eq .
( [ ipmst ] ) ]
. therefore the scattering process can be described by a stationary scattering theory at any time @xmath196 ( see ref.@xcite ) .
contrarily , in the case of a finite @xmath8 , many - body effects modify the dispersion law , so that @xmath197 , and hence @xmath198 .
the discussion above is related to the scattering of two condensate waves bounded only by a normalization box .
a more realistic yet simple model is that of a half - collision of rectangular wavepackets . consider a homogeneous condensate occupying a region of length @xmath199 along the @xmath200 axis ( leaving the results independent of size and shape in the @xmath201-plane ) . at the moment
@xmath109 a bragg scattering pulse is applied , creating a wavepacket of density @xmath202 and momentum @xmath141 directed along the @xmath200 axis , while the rest of the atoms , with density @xmath203 , remains in the @xmath14 wavepacket .
losses take place only at the overlapping parts of the two wavepackets .
the overlapping time @xmath204 varies from 0 to @xmath205 , between the leading and trailing edges of the moving wavepacket , where @xmath206 is the total duration of the collision of the two wavepackets .
the average density loss can be expressed as @xmath207 where @xmath208 is defined by eq.([npmdota ] ) .    at the regime of non - stationarity , @xmath209 , or @xmath210 ,
the loss rate is given by eq .
( [ npmdotn ] ) , and @xmath211 at the regime of stationarity , @xmath212 , the loss rate is determined by eq .
( [ ndotstat ] ) , and the loss has a form @xmath213 the same expression for the loss can also be obtained by using the method of complex scattering length @xcite .    at the regime of bose - enhancement ,
@xmath214 , the loss rate eq .
( [ ndotenh ] ) gives the average density loss @xmath215 the approximate expressions ( [ dnpmn ] ) , ( [ dnpms ] ) , and ( [ dnpmbe ] ) are compared in fig .
[ fc ] to results of numerical integration in eqs .
( [ iunstpm ] ) and ( [ dnpma ] ) .    the parametric approximation and , hence all the analysis above ,
are applicable only as long as @xmath216 . using eq .
( [ dnpmbe ] ) , this condition can be rewritten in the form of a limit on the duration of collision consistent with the approximation : @xmath217 for example , in the nist na experiment @xcite ( where @xmath218@xmath126 , @xmath219 nm , @xmath220 m , and @xmath221@xmath222 ; see ref .
@xcite ) , the parametric approximation is applicable while @xmath223 , which allows for values of the exponential bose - enhancement factor in eq .
( [ dnpmbe ] ) of more than a 100 .
this validity range substantially exceeds the range of applicability of the complex scattering - length method , @xmath224 . however , in the four - wave - mixing nist experiment @xcite , @xmath225 , which requires taking into account a considerable depletion of the condensates during the collision .
this task may be performed by combining the present approach with that of the slowly - varying - envelope approximation @xcite .
the consideration of wavepackets of a more realistic shape requires the solution of eq .
( [ csip2 ] ) for time - dependent @xmath146 and @xmath147 .
this solution can be obtained by using a perturbation theory , which is applicable as long as @xmath107 , and does not describe the effects of bose enhancement .
the evaluation is made simpler in the coordinate representation , in which case the linearized equations of motion for the field operators of the fluctuations @xmath226 have the form @xmath227 where @xmath228 is a unit - normalized wavefunction describing particles in the colliding wavepackets , @xmath229 is the total number of particles in the wavepackets , and @xmath230 is the interatomic interaction potential .
the first - order perturbative solution of eq .
( [ csir ] ) can be expressed in the form @xmath231 here the zero - order solution @xmath232 is expressed in terms of free ( plane ) waves @xmath233 and the annihilation operators @xmath22 for the fluctuation in the momentum representation . as the following equations turn out to be independent of the first term in eq .
( [ xi1 ] ) , including the function @xmath234 , it is not necessary to specify that function .
the second term includes the function @xmath235    considering scattering into the vacuum of fluctuations , such that @xmath236 , and therefore @xmath237 one can express the number of the lost atoms in the form @xmath238 where @xmath239 here @xmath240 and @xmath241 is the heavyside step function . as one can see @xmath242 is a time - dependent green function for two free particles .
therefore @xmath243 is nothing else but the second - order term for the diagonal matrix element of the transition amplitude , obtained by time - dependent perturbation theory for a collision of two particles , each in the state @xmath244 .
if the size of the wavepacket obeys @xmath245 , the result can be reduced to a form given by the stationary scattering theory ( see ref .
@xcite ) , which is proportional to the imaginary part of the elastic scattering length , in agreement with ref .  @xcite .
consider two colliding gaussian wavepackets described by the wavefunction @xmath246 where @xmath247 and @xmath248 is the number of particles in each wavepacket .
substituting eq .
( [ phig ] ) into eq .
( [ deltan ] ) one can obtain @xmath249 where @xmath250 is the collision duration .
equation ( [ deltanst ] ) coincide with the expression which can be obtained by using of the complex scattering length method @xcite and neglecting the depletion of the wavepackets .
the recent preprint @xcite , posted during preparation of the present paper , treats the loss of atoms from colliding bec wavepackets by using an approach similar to the one presented here , with solutions of the form of eq .
( [ solap ] ) published previously in ref .
@xcite . however , unjustified neglections of substantial terms in the hamiltonian lead to an expression for @xmath147 which is smaller by a factor of four than the one given here by eq .
( [ g2 ] ) .
this neglection should therefore result in an underestimation of the condensate loss .    the above - mentioned work , ref .
@xcite , considers collision of gaussian wavepackets to first - order perturbation theory only , which do not describe the effects of bose enhancement .
finite analytical expressions of the kind presented here were not derived there , and the results of numerical calculations disagree with results of the complex scattering length method @xcite even after the correcting for the missing factor four in @xmath147 are taken into account .
in contrast , the analytical expression eq .
( [ deltanst ] ) of the present paper , based on the same approximations , is in an agreement with the complex scattering length method .
fluctuations in bec are treated by using a parametric approximation for the exact quantum equations of motion for the field operators .
depending on the excitation energy , the fluctuations are divided into stable modes , which can be represented as bogolubov quasiparticles ( field oscillators ) , and unstable modes , which can be represented as inverted oscillators and demonstrate a rapid growth .
initially the dynamics of unstable modes is a quantum process of spontaneous decay into the vacuum of the fluctuations with a super - exponential rate .
this approach is applied to two physical problems .
the first one involves the dynamics of a bec after a fast change of the scattering length from a non - negative to a negative one . in this case , the experimentally observed @xcite delay of the condensate collapse is described .
the second problem involves the elastic scattering losses on collision of two bec wavepackets . at short collision times
[ see eq.([ndotstat ] ) ] the losses obtained here are in agreement with those of the method of complex scattering length @xcite , with specific reference to the collisions of non - confined waves , of rectangular and of gaussian wavepackets . at longer collision times
[ see eq .
( [ ndotenh ] ) ] quantum effects lead to bose - enhancement of the elastic scattering losses .
the atoms lost from the condensates in both cases form correlated pairs in optimally entangled and optimally squeezed states .
the author is very grateful to p.  s.  julienne and k.  burnett for helpful discussions .
special appreciation is expressed to a.  ben - reuven for invaluable discussions and the help in preparation of the present paper .
this work was partially supported by the national science foundation through a grant for the institute for theoretical atomic and molecular physics at harvard university and the smithsonian astrophysical observatory .
parkins and d.f .
walls , phys .
rep . * 303 * , 1 ( 1998 ) .
f. dalfovo , s. giorgini , l.p .
pitaevskii , and s. stringari , rev .
phys . * 71 * , 463 ( 1999 ) .
e. m. lifshitz and l. p. pitaevskii , _ statistical physics , pt . 2 _ ( pergamon , oxford , 1980 ) . a. l. fetter and j. d. walecka , _ quantum theory of many - body systems _
( mcgraw , san francisco , 1971 ) .
n. n. bogolubov and n. n. bogolubov , jr .
_ an introduction to quantum statistical mechanics _
( gordon and breach , lausanne , 1994 ) . j .-
blaizot and g. ripka , _ quantum theory of finite systems _
( mit press , cambridge , ma , 1986 ) .
n. p. proukakis and k. burnett , j. res .
nist * 101 * , 475 ( 1996 ) .
n. p. proukakis , k. burnett , and h. t. c. stoof , phys .
rev a * 57 * , 1230 ( 1998 ) .
a. griffin , phys .
b * 53 * , 9341 ( 1996 ) .
r. walser , j. williams , j. cooper , and m. holland , phys .
a * 59 * , 3878 ( 1999 ) .
sackett , j.m .
gerton , m. welling and r.g .
hulet , phys .
lett . * 82 * , 876 ( 1999 ) .
j. m. gerton , d. strekalov , i. prodan , and r. g. hulet , nature ( london ) * 408 * , 692 , ( 2000 ) .
j. l. roberts _ et al .
* 86 * , 4211 ( 2001 ) .
e. a. donley _
et al . _ , nature ( london ) * 412 * , 295 ( 2001 ) .
e. timmermans , p. tommasini , m .
hussein , and a. kerman , phys .
* 315 * , 199 ( 1999 ) .
duine and h. t. c. stoof , phys .
lett . * 86 * , 2204 ( 2001 ) .
a. vardi , v. a. yurovsky , and j. r. anglin , cond - mat/0105439 ; phys .
rev a ( in press ) .
v. a. yurovsky , a. ben - reuven , and p. s. julienne , cond - mat/0108372 ; phys .
( submitted ) .
m. o. scully and m. s. zubairy , _ quantum optics _
( university press , cambridge , 1997 ) .
m. g. moore and p. meystre , phys .
a * 59 * , r1754 ( 1999 ) ; c. k. law and n. p. bigelow , phys . rev a * 58 * , 4791 ( 1998 ) .
h. a. haus and c. x. yu , j. opt .
b * 17 * , 618 ( 2000 ) . v. v. dodonov and v. i. manko , in _ invariants and evolution of nonstationary quantum systems _ , proceedings of lebedev physics institute vol . 183 , edited by m. a. markov ( nova science publishers , commack , ny , 1989 ) , p. 3 ; o. v. manko , phys .
a * 228 * , 29 ( 1997 ) .
a. vardi and j. r. anglin , phys .
lett . * 86 * , 568 ( 2001 ) .
j. r. anglin and a. vardi , phys .
rev a * 64 * , 013605 ( 2001 ) .
u. v. poulsen and k. molmer , phys .
rev a * 63 * , 023604 ( 2001 ) .
d. c. roberts , t. gasenzer , and k. burnett , quant - ph/0106060 .
et al . _ , nature ( london ) * 398 * , 218 ( 1999 ) . y. b. band , m. trippenbach , j. p. burke , and p. s. julienne , phys . rev
. lett . * 84 * , 5462 ( 2000 ) .
m. trippenbach , y. b. band , and p. s. julienne , phys .
a * 62 * , 023608 ( 2000 ) .
r. bach , m. trippenbach , and k. rzazewski , cond - mat/0107117 .
m. olshanii and l. pricopenko , cond - mat 0101275 .
m. lewenstein and l. you , phys .
lett . * 77 * 3489 , ( 1996 ) .
r. j. dodd _
et al _ , phys .
a * 54 * , 661 ( 1996 ) .
h. bateman and a. erdely , _ higher transcendental functions , vol . 2 _
( mc graw - hill , new york , 1953 ) .
m. l. goldberger and k. w. watson , _ collision theory _
( wiley , new york , 1964 ) .


comprehension of the phenomenon of _ decoherence _ in open quantum systems has always attracted much attention , in particular as a prerequisite to understand the transition from quantum to classical behavior .
the dissipative two - state or _ spin - boson _ model has been thoroughly studied in wide regions of the parameter space with diverse methods and techniques since the 80 s @xcite .    in the last decade , the subject of decoherence has experienced renaissance following the growing interest in the field of quantum state manipulation and quantum computation @xcite .
any noise source sensitively leads to a narrowing of the quantum coherence domain .
this entails severe limitations for coupled qubits to perform logic quantum operations . for this reason ,
extensive understanding of the decoherence mechanisms is indispensable .    in this work ,
we focus upon a model which is a generalization of the single spin - boson model to the case of two spins which mutually interact via an ising - type coupling and are coupled to independent environments made up by bosons .
the first analysis of this model relying on the influence functional method was given by dub and stamp @xcite .
they obtained results for the dynamics in analytic form in restricted regions of the parameter space by omitting certain classes of path contributions and bath correlations .
several other previous studies on the same or related models relied on the master equation and/or perturbative redfield approach @xcite .
besides the weak - coupling assumption , often the secular approximation @xcite is made , which breaks down however when the spectrum becomes degenerate .
the model allows , for instance , to study decoherence and relaxation of two coupled qubits @xcite , or the influence of a bistable impurity on the qubit dynamics @xcite . the latter may significantly degrade coherence in josephson phase qubits @xcite
other possible application is study of coherence effects in coupled molecular magnets @xcite .    in earlier works ,
the model has been analyzed in the pure dephasing regime both by the feynman - vernon method @xcite and the lindblad approach @xcite .
here we extend the work in ref .
@xcite beyond the pure dephasing regime and include the full dynamics of the qubit .
in particular , we are interested in the competition between decoherence and relaxation to the equilibrium state . here
we focus on the white - noise regime .
we shall derive the exact solution for the reduced density matrix without restriction on the parameters of the model and analyze it in the coherent and incoherent domains and in the crossover regions in between .
the model and relevant quantities of the reduced dynamics are introduced in section 2 .
section 3 deals with the exact formal solution for the reduced dynamics . in section 4 , the path sum is carried out in the white - noise domain without any further approximation , and analytic expressions for the relevant expectation values in laplace space are presented . after an overview of the qualitative features of the dynamics in section 5 , we present in section 6 explicit expressions for decoherence and relaxation in the various parameter regimes ranging from low temperture and/or weak coupling to high temperature and/or strong coupling .
finally , we study in section 7 the influence of a nonlinear spin - boson environment on the second spin in the various limits .
we demonstrate that it behaves in the weak - coupling limit as a bosonic ( linear ) bath with a resonant spectral structure .
we consider two two - state systems which are coupled to each other via an ising - type coupling and to independent bosonic environments . in pseudospin representation , we choose the generalized spin - boson hamiltonian ( we use units where @xmath0 ) @xmath1 in the basis formed by the localized eigen states @xmath2 and @xmath3 of @xmath4 and @xmath5 , respectively , @xmath6 and @xmath7 represent the tunneling couplings between the localized states , and the coupling term @xmath8 acts as a mutual bias energy of strength @xmath9 . the collective bath modes
@xmath10 $ ] ( @xmath11 ) represent fluctuating bias forces .
the hamiltonian is very rich in content and may model diverse physical situations .
it may describe two coupled qubits or a qubit @xmath12 in contact with a complex environment formed by a bistable dissipative impurity @xmath13 .
other possible realizations are coupled molecular magnets of which the low - energy states can be viewed as a spin @xcite .    for the model ( [ ham1 ] ) , all effects of the environments are captured by the power spectrum of the collective bath modes @xmath14 with the spectral density of the coupling @xcite @xmath15 here , the second form represents the ohmic case with a high - frequency cut - off @xmath16 .
alternatively , one may choose that the two spins are coupled to a common bath @xcite .
here we study the effects of independent environments .
this case is realistic in most physical systems of actual interest .
the density matrix of a single spin has four matrix elements , the two populations that we shall label as @xmath17 and @xmath18 , and the two coherences with labels @xmath19 and @xmath20 .
the two - spin density matrix has 16 matrix elements @xmath21 .
we choose for convenience that the first ( second ) index refers to the states @xmath22 ( @xmath23 ) of the @xmath12spin ( @xmath13spin ) .
the matrix elements can be expressed in terms of expectation values of 15 operators , @xmath24 , @xmath25 and @xmath26 ( @xmath27 and @xmath28 ) .
the 4 pure populations may then be written as @xmath29/4 \ ; , \\[3 mm ] \rho_{\rm 1,3}^{}(t ) & = & [ \ , 1 \,+\,\langle\sigma_z^{}\rangle_t \,-\,\langle\tau_z^{}\rangle_t \,-\,\langle\sigma_z^{}\tau_z^{}\rangle_t\,]/4 \ ; , \\[3 mm ] \rho_{\rm 3,1}^{}(t ) & = & [ \ , 1 \,-\,\langle\sigma_z^{}\rangle_t \,+\,\langle\tau_z^{}\rangle_t \,-\,\langle\sigma_z^{}\tau_z^{}\rangle_t \,]/4\ ; , \\[3 mm ] \rho_{\rm 3,3}^{}(t ) & = & [ \ , 1 \,-\,\langle\sigma_z^{}\rangle_t \,-\,\langle\tau_z^{}\rangle_t \,+\,\langle\sigma_z^{}\tau_z^{}\rangle_t \,]/4\;.   \end{array}\ ] ] corresponding expressions hold for the 4 pure coherences and the 8 hybrid states . for instance , we have @xmath30/4 \ ; .\ ] ] here we are predominantly interested in the populations .
throughout we will choose that the reduced system starts out from the initial state @xmath31 while the heat reservoirs are in thermal equilibrium at temperature @xmath32 .    in the absence of the environment
, the hamiltonian @xmath33 can be easily transformed into diagonal form @xmath34 the eigenfrequencies are @xmath35 @xmath36 and they obey the vieta relations @xmath37 \omega_-^2+\omega_+^2 & = & 2 (   \delta_1 ^ 2+\delta_2 ^ 2+v^2 _ { } ) \ ; , \quad\ ; \omega_+^2\omega_-^2 \;=\ ; ( \delta_1 ^ 2+\delta_2 ^ 2+v^2_{})^2_{}- 4\delta_1 ^ 2\delta_2 ^ 2 \ ; .
\end{array}\ ] ]    the liouville equations @xmath38 $ ] ( @xmath39 ) , where the set @xmath40 represents the above 15 operators , yield 15 coupled equations .
these are conveniently solved in laplace space .
for instance , we get @xmath41    \langle \sigma_z\tau_z ( \lambda ) \rangle & = & { \displaystyle \frac{(v^2 + \lambda^2 ) ( v^2 + \delta_1 ^ 2 + \delta_2 ^ 2 + \lambda^2)}{\lambda\ , ( \lambda^2+\omega_{+}^{2})(\lambda^2+\omega_{-}^{2 } ) } } \ ; .   \end{array}\ ] ]    here we are interested in the evolution of the two - spin system without restricting ourselves to weak damping .
therefore , we refrain from employing the perturbative redfield approach .
rather we calculate the reduced dynamics with use of the feynman - vernon influence functional method .
we show that the solution is available in analytic form in the white noise limit for general parameters @xmath42 and @xmath32 .
within the feynman - vernon method , the exact formal expression for the rdm of the two - spin system is the quadruple path integral @xmath43\ , { \mathcal a}^ * [ \sigma']\ , { \mathcal a}[\tau]\ , { \mathcal a}^ * [ \tau ' ] \,{\cal b}[\sigma,\sigma';\tau,\tau ' ] \,\mathcal{f}[\sigma,\sigma';\tau,\tau']\ ] ] with appropriately chosen boundary values for the spin paths . here , each of the paths @xmath44 starts out from the localized state @xmath2 at time zero .
they end up at time @xmath45 in the states @xmath46 , respectively , where @xmath47 .
the functional @xmath48 $ ] is the amplitude for the free spin @xmath12 to follow the path @xmath49 , the functional @xmath50 $ ] represents the coupling of the two spins ( see below ) , and the functional @xmath51 $ ] introduces the environmental influences .    for uncorrelated baths
, we have @xmath52\,=\,\mathcal{f}_1[\xi_1^{},\eta_1^{}]\,\,\mathcal{f}_2[\xi_2^{},\eta_2^{}]$ ] , where @xmath53 = \int_0^t \!\!\!dt ' \int_0^{t'}\!\!\ !
dt '' \left [   \dot{\xi}_\zeta ( t')q_\zeta'(t'-t'')\dot{\xi}_\zeta(t '' )   \,+\,{\rm i}\,\dot{\xi}_\zeta ( t')q_\zeta''(t'-t'')\dot{\eta}_\zeta(t '' ) \right ] \;.\ ] ] here we have introduced symmetric and antisymmetric spin paths , @xmath54\ ; , \hspace{1 cm } \eta_1^{}(t')\;=\ ; \frac{1}{2}[\sigma(t')+\sigma'(t ' ) ] } \;,\\[3 mm ] \xi_2^{}(t ' ) & = & { \displaystyle \frac{1}{2}[\tau(t')-\tau'(t')]\ ; , \hspace{1 cm } \eta_2(^{}t')\;=\;\frac{1}{2}[\tau(t')+\tau'(t ' ) ] } \ ; . \end{array}\ ] ] the correlator @xmath55 is the second integral of the force autocorrelation function @xmath56 ( see eq .
( [ powerspec ] ) ) . in the ohmic scaling limit
, we have @xmath57   + i \pi k_\zeta \ , { \rm sgn}(t)\ ; , \qquad \zeta=1,\,2\ ; .\ ] ] here , @xmath58 is the usual dimensionless ohmic coupling strength for the spin @xmath59 , and @xmath60 is the inverse temperature .    to handle the quadruple path integral ( [ rdm ] ) , we follow the procedure for the single spin - boson problem @xcite and write it as an integral over two paths , one for each spin .
each such path visits the diagonal `` sojourn '' states and the off - diagonal `` blip '' states of the respective spin
. a path which starts and ends in a sojourn state must contain an even number of transitions with amplitude @xmath61 for each flip of spin @xmath59 .
the flips occur at times @xmath62 for spin 1 and at times @xmath63 for spin 2 . upon labeling the sojourn and
blip states with charges @xmath64 ( for spin 1 ) and @xmath65 ( for spin 2 ) , each with values @xmath66 , the paths with @xmath67 and @xmath68 transitions , respectively , may be written as @xmath69\ ; , \\ \xi^{(n_1^{})}_1(t ' ) & = &   \sum_{j=1}^{n_1^ { } } \xi_{1,j}^ { }   \left[\theta(t'-t_{2j-1})-\theta(t'-t_{2j } ) \right]\ ; ,   \\
\eta^{(n_2^{})}_2(t ' ) & = & \sum_{j=0}^{n_2^ { } } \eta_{2,j}^ { } \left[\theta(t'-s_{2j}^{})-\theta(t'-s_{2j+1}^ { } ) \right]\ ; , \\ \xi^{(n_2^{})}_2(t ' )   & = & \sum_{j=1}^{n_2^ { } } \xi_{2,j}^ { }   \left[\theta(t'-s_{2j-1}^{})-\theta(t'-s_{2j}^ { } ) \right]\ ; .
\end{aligned}\ ] ] upon introducing the notation @xmath70 , we may write the bath correlations between the blip pair @xmath71 of spin @xmath59 in the compact form @xmath72 with this , the influence functional for the paths ( [ sojblip ] ) reads @xmath73\ ,   \exp\big[-\sum_{j=2}^{n_\zeta^{}}\sum_{k=1}^{j-1}\xi_{\zeta j}^{}\xi_{\zeta , k}^{}\lambda_{\zeta;j , k}^{}\big ] } \\[4 mm ]   & & \times\ ; { \displaystyle    \,\exp\big[{\rm i } \,\pi k_\zeta^{}\sum_{k=0}^{n_\zeta^{}-1 } \xi_{\zeta , k+1}^{}\eta_{\zeta , k}^ { }     \big ] } \ ; .
\end{array}\ ] ] here the first and second term represent the intrablip and interblip correlations , respectively .
the phase term is specific to the ohmic scaling limit and represents correlations of the sojourns with their subsequent blips .
the sum over all paths now means ( i ) to sum over all possible intermediate sojourn and blip states of the two spins the paths with a given number of transitions can visit , ( ii ) to integrate over the ( for each spin ) time - ordered jumps of these paths , and ( iii ) to sum over the possible number of transitions the two spins can take , @xmath74 & \times & { \displaystyle \sum_{n_2^{}=0}^\infty   \sum _ { \{\eta_{2,j}^{}=\pm 1\ } }   \sum _ { \{\xi_{2,j}^{}=\pm 1\ } }   \int_{t_{\rm i}^{}}^{t_{\rm f}^ { } } \!\ ! { \rm d}s_{2n_2^{}}^ { } \int_{t_{\rm i}^{}}^{s_{2n_2^ { } } } \!\ !   { \rm d}s_{2n_2^{}-1}^ { }   \cdots   \int_{t_{\rm i}^{}}^{s_2^ { } }   \!\!{\rm d}s_1^{}\;\cdots    }    \ ; . \end{array}\ ] ] for coherences , the number of transitions in the respective spin path is odd .
next , we rewrite the double path sum ( [ pathsum ] ) with time - ordering for each spin in terms of a single path over the 16 possible states with time ordering of all the flip times @xmath75 . in this representation ,
the system lingers for some period in a particular state of the double - spin rdm and then it flips with amplitude @xmath76 or @xmath77 to another state . for the longitudinal spin - spin coupling @xmath78 in ( [ ham1 ] ) , the coupling factor @xmath79 in eq .
( [ rdm ] ) is unity when the system dwells on one of the four pure soujourn or one of the four pure blip states , and it is sensitive to the coupling when it stays , say for a period @xmath80 , in one of the eight hybrid states . in detail , we have @xmath81 & = & { \rm e}^{+i v u_j^{}}_{}\ ; , \quad\quad\mbox{for}\quad\quad j \;\in\ ; ( 1,2),\,(2,1),\,(4,3),\,(3,4)\ ; , \\[3 mm ] { \cal b}[\,j\ , ] & = & { \rm e}^{-i v u_j^{}}_{}\ ; , \quad\quad\mbox{for}\quad\quad j \;\in\ ; ( 1,4),\,(4,1),\,(2,3),\,(3,2)\ ; .   \end{array}\ ] ]    combination of the above expressions yields the exact formal solution for the dynamics of the rdm of the two - spin system in the ohmic scaling limit . evidently , because of the nonconvolutive form of the bath correlations in the influence functional ( [ infl ] ) , the path sum can not be performed in analytic form .
alternatively , one may recast the exact formal series expression for the populations in the form of generalized master equations in which the kernels , by definition , are the irreducible components of path segments with diagonal initial and final states @xcite . in the general case ,
the kernels are given by an infinite series in @xmath82 and @xmath83 , with the time integrals in each summand being again in nonconvolutive form .
additional difficulties in performing the path sum ( [ pathsum ] ) arises from the spin - spin coupling ( [ intfac ] ) .
in the white - noise limit of eq .
( [ powerspec ] ) , @xmath84 , where @xmath85 is a scaled thermal energy , the bath correlation function takes the form accounts for deviation of the actual high frequency behavior of @xmath86 from the exponential cut - off form in eq .
( [ specdens ] ) @xcite . ]
@xmath87 \,+\ , \vartheta_\zeta^ { } |t|   \,+\ , { \rm i}\ , \pi k_\zeta \ , { \rm sgn}(t ) \;.\ ] ] this expression emerges directly from eq .
( [ ohmcorr ] ) in the high temperature or long - time limit @xmath88 .
the first term in eq .
( [ wncorr ] ) leads to an adiabatic ( franck - condon - type ) renormalization factor made up by modes in the frequency range @xmath89 .
it is natural to assimilate this term , together with the phase term , into an effective temperature - dependent tunneling matrix element , @xmath90 where @xmath91 is the standard renormalized tunneling matrix element .
all dynamical effects of the environmental coupling are captured by the second term @xmath92 in eq .
( [ wncorr ] ) . from this
we see that the weight of the thermal energy relative to the systems energies @xmath93 and @xmath9 is assessed by the scaled thermal energy @xmath94 .
based on experience from the single spin - boson system we should expect that the form ( [ wncorr ] ) is a remarkably good approximation in the parameter range @xcite @xmath95 and this is corroborated indeed by our study . for a single unbiased spin with ohmic damping @xmath96 , the coherent - incoherent `` phase''-transition is at temperature @xmath97 with @xmath98 @xcite .
therefore we should expect , that , for @xmath99 , the white noise form ( [ wncorr ] ) is valid not only in the incoherent regime but also in a sizeable domain of the coherent regime .
this shall be confirmed subsequently .    for the form ( [ wncorr ] ) of @xmath100
, the interblip correlations cancel out exactly in eq .
( [ blipcor ] ) , @xmath101 . as a result ,
each term of the infinite series for the rdm becomes a convolution .
this makes the path sum accomplishable .
we shall now exemplify , by taking @xmath102 and @xmath103 as examples , that the path sum for the laplace transform of the rdm can be carried out exactly in analytic form .
this is achieved by first calculating the kernels and then summing up the respective geometrical series of these objects .      by definition ,
the kernels represent irreducible path segments which interpolate between pure sojourn states .
irreducibility means that these segments can not be separated into uncorrelated pieces without removing bath correlations .
the analysis gives that every contribution to the kernel @xmath106 of @xmath107 displays initially and finally a transition of spin @xmath12 with any number of even hops of spin @xmath13 at intermediate times , as shown in the diagrams of fig .
there are no other contributions .
for all other irreducible diagrams , one may think of , e.g. , where either the first or the last flip , or both of them , are flips of spin @xmath13 , the respective contributions from the two different final states of spin @xmath13 cancel each other .
sketch of the irreducible kernels ( with amputated legs ) @xmath108 ( a ) , @xmath109 ( b ) , and @xmath110 ( c ) .
the solid and dashed curve represent the correlations of bath 1 and of bath 2 , respectively . ]
it is convenient to write @xmath111 , where @xmath112 is of order @xmath113 .
the sum over all spin states of order @xmath114 yields the expressions @xmath115 @xmath116 & & \!\!\!\!\times \ ; { \displaystyle   \,\frac{1}{2}\,[\ , \cos(\pi k_1^ { } + v\tau_1^ { } + v \tau_3^ { } ) -    \cos(\pi k_1^ { } + v\tau_1^ { } - v \tau_3^{})\ , ] \cos(\pi\,k_2^ { } ) } \ ; , \end{array}\ ] ] @xmath117 here we have taken into account that , for the white - noise form ( [ wncorr ] ) , a correlation of bath @xmath59 stretching over an interval between neighboring hops effectively leads to a shift of the laplace variable in the respective time integral , @xmath118 .
it is convenient to split the kernels into the contributions which are even and odd in the coupling @xmath9 , @xmath119 .
the resulting expressions may be written as @xmath120 & \times & { \displaystyle   \!\!\!\left\ { 1 \;+\ ; \frac{v^2_{}}{(\lambda + \vartheta_1^{})(\lambda + \vartheta_1^{}+\vartheta_2^{})}\ ,        \frac{\bar\delta_{2}^2 } {   ( \lambda + \vartheta_1^{})^2_{}+ v^2 _ { }    +    \frac{\lambda+\vartheta_1^{}}{\lambda + \vartheta_1^ { } + \vartheta_2^{}}\bar\delta_2^{2 } } \right\ }   } \ ; , \end{array}\ ] ] @xmath121 & \times & { \displaystyle   \!\!\!\left\ { 1 \;-\ ; \frac{\lambda+\vartheta_1^{}}{\lambda + \vartheta_1^{}+\vartheta_2^{}}\ ,        \frac{\bar\delta_{2}^2 } {   ( \lambda + \vartheta_1^{})^2_{}+ v^2 _ { }    +    \frac{\lambda+\vartheta_1^{}}{\lambda + \vartheta_1^ { } + \vartheta_2^{}}\bar\delta_2^{2 } } \right\ }   } \ ; . \end{array}\ ] ] paths which visit a pure sojourn state at intermediate times yield reducible contributions .
taking into account all possibilities of such visits yields a geometrical series in the kernel @xmath122 , while @xmath123 occurs only once as initial irreducible contribution .
thus we get for @xmath107 the concise form @xmath124 algebraic manipulation gives @xmath107 in the form of a simple fraction with denominator and numerator in the form of polynomials of degree four .
we get @xmath125 & + & \!\!\!(v^2_{}\vartheta_1^ { } + \bar\delta_2 ^ 2\vartheta_1^ { } + \vartheta_1 ^ 3 + v^2_{}\vartheta_2^ { }   + \vartheta_1 ^ 2\vartheta_2^{})\lambda + v\bar\delta_1^{}\tan(\pi k_1^ { } ) ( \vartheta_1^{}+\vartheta_2^ { } + \lambda )   \ ; ,   \end{array}\ ] ] @xmath126 & + & ( v^2\ , \vartheta_1^ { } + 2\ , \bar\delta_{1}^2 \,\vartheta_1^ { } + \bar\delta_{2}^2 \
, \vartheta_1^ { } +   \vartheta_1 ^ 3 +     v^2 _ { } \,\vartheta_2^ { } + \bar\delta_{1}^2\ , \vartheta_2^ { } + \vartheta_1 ^ 2\ , \vartheta_2^{})\,\lambda   \\[3 mm ] & + & \bar\delta_{1}^2 \,\vartheta_1 ^ 2 + \bar\delta_{1}^2\ , \vartheta_1^{}\ , \vartheta_2^ { } +   \bar\omega_{}^2 \,\bar\delta_{}^2   \ ; .
\end{array}\ ] ] here we have introduced the eigenfrequencies @xmath127 , @xmath128 of the undamped coupled two - spin system .
the bar denotes adiabatic renormalization in eq .
( [ def_eigen1 ] ) according to eq .
( [ adia ] ) .
the pole at @xmath129 in eq .
( [ expecsigma ] ) yields the equilibrium value @xmath130 hence @xmath131 is negligibly small for @xmath132 , while , as @xmath133 , it takes the proper ( white noise ) equilibrium value @xmath134 of the single biased spin boson system .
the dynamical poles @xmath135 ( @xmath136 ) are given by a quartic equation with real coefficients .
upon collecting the various pole contributions , we get in the time domain @xmath137 evidently , the expectation @xmath138 takes similar form , @xmath139 here , the polynomials @xmath140 and @xmath141 follow from the expressions ( [ sz1 ] ) and ( [ sz2 ] ) by interchange of the indices 1 and 2 .
sketch of the irreducible kernels @xmath142 ( left ) and @xmath143 ( right ) .
the intervals are dressed by self - energy contributions of spin @xmath12 ( black circle ) and spin @xmath13 ( black square ) , as sketched in fig .
[ f4 ] . ]
self - energy terms due to spin @xmath12 ( circle ) and spin @xmath13 ( square ) . ]      consider first contributions to the kernel of @xmath145 in which the first and last flip are made by spin @xmath12 , as sketched in fig .  [ f2 ] . in the intervals of the bare diagrams , either spin @xmath12 or spin @xmath13 or both stay offdiagonal .
every interval in which spin @xmath59 dwells in a sojourn state is dressed by selfenergy contributions schematically given for spin @xmath12 ( circle ) and spin @xmath13 ( square ) in fig .
diagram [ f2](a ) yields @xmath146 \,\alpha_1^{}(\lambda )   \ ; , \end{aligned}\ ] ] where the functions @xmath147 @xmath148 are given by @xmath149 the nested diagram in fig .
[ f2](b ) produces the additional factor @xmath150 , @xmath151 higher - order nested diagrams of the type sketched in fig .
[ f2 ] class into a geometrical series .
all these terms are readily summed up to the contribution @xmath152 similarly , we find that diagrams ( a ) and ( b ) in fig .
[ f3 ] yield the expressions @xmath153 \,\alpha_1^{}(\lambda)\alpha_2^{}(\lambda )   \ ; , \ ] ] @xmath154 with all higher - order nested diagrams of this type added , one finds again a geometrical series in @xmath155 , @xmath156 clearly , we must add those contributions resulting from the terms @xmath157 and @xmath158 by interchange of the spins @xmath12 and @xmath13 .
the analysis shows that there are no other contributions .
again , we split the kernel @xmath159 into the parts which are even and odd in the coupling @xmath9 .
we readily get for @xmath160 $ ] the forms @xmath161 @xmath162 these expression represent the entity of irreducible path segments .
next , we observe that the sum of two - spin paths with any number of interim visits of pure sojourn states yields a geometrical series of these objects . in the part which is odd in @xmath9 , the first irreducible path section is again described by the kernel @xmath163 .
thus we get @xmath164 the second form is a simple fraction with polynomials @xmath165 and @xmath166 of sixth order , @xmath167   \\[2 mm ] & & + \ ; v \tan ( \pi k_1^ { } ) \ , \bar\delta_1 ^ 2\,(\vartheta_2^{}+\lambda)\,[\,\bar\delta_1 ^ 2\ , -\
, \bar\delta_2 ^ 2\ , + \ , ( \vartheta_2+\lambda)(\vartheta_1^{}+\vartheta_2^{}+\lambda )   \ , ]       \\[2 mm ]   & & + \;v \tan ( \pi k_2^ { } ) \ , \bar\delta_2 ^ 2\,(\vartheta_1^{}+\lambda)\,[\,\bar\delta_2 ^ 2 \,-\ , \bar\delta_1 ^ 2 \,+\ , ( \vartheta_1+\lambda ) ( \vartheta_1^{}+\vartheta_2^{}+\lambda )   \ , ]       \\[2 mm ]     & & + \ ; \lambda \ ,
v^4_{}(\vartheta_1^{}+\vartheta_2^{}+\lambda )   \,+\,\lambda v^2_{}[\,\bar\delta_2 ^ 2\,(\vartheta_1^{}+\lambda ) +   \bar\delta_1 ^ 2\,(\vartheta_2^{}+\lambda)\ , ] \\[2 mm ] & & + \ ; \lambda \,v^2 _ { }   ( \vartheta_1^{}+\vartheta_2^{}+\lambda ) [ \,\vartheta_1 ^ 2+\vartheta_2 ^ 2 + 2\lambda ( \vartheta_1^{}+\vartheta_2^{}+\lambda)\,]\\[2 mm ] & & + \ ; \lambda \,(\vartheta_1^{}+\lambda)(\vartheta_2^{}+\lambda ) [ \,\bar\delta_1 ^ 2\,(\vartheta_1+\lambda ) + \bar\delta_2 ^ 2\,(\vartheta_2+\lambda)\ , ] \\[2 mm ] & & + \ ; \lambda\,(\vartheta_1^{}+\lambda)^2_{}(\vartheta_2^{}+\lambda)^2_{}(\vartheta_1^{}+\vartheta_2^{}+\lambda ) \ ; , \end{array}\ ] ] @xmath168 & & + \;(\vartheta_1 + \vartheta_2)(2\bar\omega_+^2 + 2 \bar\omega_-^2 + \vartheta_1 ^ 2 + 6\vartheta_1\vartheta_2 + \vartheta_2 ^ 2)\,\lambda^3 _ { } \\[2 mm ] & & + \;[\,\bar\omega_+^2\,\bar\omega_-^2 + \vartheta_1 \vartheta_2(2 \vartheta_1 + \vartheta_2)(\vartheta_1 + 2\vartheta_2 ) + v^2 ( 3\vartheta_1 ^ 2 + 4\vartheta_1\vartheta_2 + 3\vartheta_2 ^ 2 ) \\[2 mm ] & & + \ ; \bar\delta_1 ^ 2(2\vartheta_1 + \vartheta_2)(\vartheta_1 + 3\vartheta_2 ) +     \bar\delta_2 ^ 2(2\vartheta_2 + \vartheta_1)(\vartheta_2 + 3\vartheta_1 ) ] \,\lambda_{}^2\\[2 mm ] & & + \;(\vartheta_1 + \vartheta_2)[\,\bar\omega_+^2\ , \bar\omega_-^2   + \vartheta_1 ^ 2\vartheta_2 ^ 2   + ( v^2 + \bar\delta_2 ^ 2)\vartheta_1 ^ 2 + ( v^2_{}+\bar\delta_1 ^ 2)\vartheta_2 ^ 2 \\[2 mm ] & & + \ ; 3\vartheta_1^{}\vartheta_2^{}(\bar\delta_1 ^ 2+\bar\delta_2 ^ 2 ) \,]\,\lambda   \;+ \ ; v^2(\vartheta_1 + \vartheta_2)(\bar\delta_1 ^ 2\vartheta_1 + \bar\delta_2 ^ 2 \vartheta_2)\\[2 mm ] & & + \ ;   \vartheta_1\vartheta_2 [ \ , ( \bar\delta_1 ^ 2 - \bar\delta_2 ^ 2)^2   + ( \bar\delta_2 ^ 2\vartheta_1 + \bar\delta_1 ^ 2\vartheta_2 ) ( \vartheta_1^{}+\vartheta_2^ { } )   \ , ] \ ; . \end{array}\ ] ]    sketch of the irreducible kernels @xmath169 ( left ) and @xmath170 ( right ) .
again , the intervals are dressed as sketched in fig .
[ f4 ] . ]    the odd powers of the pole function can be removed with a shift .
putting @xmath171 \equiv \tilde d(x)$ ] , we obtain a polynomial with even powers , @xmath172\,x^4 _ { } \\[3 mm ] & & + \;[\,\bar\omega_+^2\bar\omega_-^2\ , + \,(\bar\delta_1 ^ 2+\bar\delta_2 ^ 2 - 2v^2_{})\vartheta_1^ { } \vartheta_2^ { }
-\bar\delta_1 ^ 2\vartheta_1 ^ 2-\bar\delta_2 ^
2\vartheta_2 ^ 2 \\[3 mm ] & & \quad\ ; + \;\frac{1}{16 }   ( \vartheta_1^{}-\vartheta_2^{})^2_{}(3\vartheta_1 ^ 2 + 2\vartheta_1^{}\vartheta_2^{}+3\vartheta_2 ^ 2 ) \,]\,x^2 _ { } \\[3 mm ] & & -\frac{1}{4 } [ \,v^2_{}(\vartheta_1^{}+\vartheta_2^ { } ) \,+\ , ( \bar\delta_1 ^ 2-\bar\delta_2 ^ 2)(\vartheta_2^{}-\vartheta_1^ { } ) \,+\,\frac{1}{4}(\vartheta_1 ^ 2-\vartheta_2 ^ 2)(\vartheta_1^{}-\vartheta_2^ { } ) \,]^2 _ { }   \;. \\[3 mm ] \end{array}\ ] ] thus we have in the time domain @xmath173 where the @xmath174 are the zeros of @xmath175 , and the equilibrium value is @xmath176    the expressions ( [ expecsigma ] )  ( [ expectau ] ) and ( [ expecsigtau ] )  ( [ sigtauz3 ] ) are the main results of this work .
they represent the exact analytical solutions for @xmath107 , @xmath177 , and @xmath145 , in the white - noise limit for general coupling @xmath9 and general effective reservoir couplings @xmath178 and @xmath179 . except for use of the form ( [ wncorr ] )
, no other approximation has been made .
we remark that for all other initial and final states of the rdm we would find the same pole functions ( [ sz2 ] ) and ( [ sigtauz1 ] )
. only the numerator function would be different .
the behaviors of the four dynamical poles of @xmath180 and the six dynamical poles of @xmath181 , and the respective amplitudes are quite multifarious . in this section , we sketch the characteristics for the symmetric system , @xmath182 and @xmath183 .      in the coupling range @xmath184
, there are three crossover temperatures , denoted by @xmath185 , and @xmath186 ( see figs .  [ fig : sz_crosstemp ] and [ fig : sz_re_im ] )
.    in the regime @xmath187 the dynamics is coherent and described by a superposition of two damped oscillations . for @xmath188 ,
the oscillations have different frequency and the same damping rate , and the amplitudes are comparable in magnitude . on the other hand , in the range @xmath189 , they have the same frequency , but different decrement and the amplitude belonging to the larger decrement is negligibly small .    in the temperature regime @xmath190 ,
the dynamics is incoherent . in the regime @xmath191 ,
the four poles are real , and the two smallest rates have largest amplitudes and dominate the relaxation process .    in the so - called kondo regime @xmath192 ,
the dominant pole is real and approaches @xmath193 , and its residuum goes to @xmath194 , as temperature is increased .
the other real pole takes the value @xmath195 , while its residuum drops to zero .
there is also a damped oscillation of which the frequency and rate approach asymptotically @xmath9 and @xmath196 , but the amplitude becomes negligibly small .
the phenomenon that , in the kondo regime for @xmath197 , incoherent relaxation slows down with increasing temperature , is already well - known in the single spin - boson problem @xcite .
[ fig : sz_fourier_time ] shows the transition from coherent to incoherent dynamics as @xmath196 is raised .
one can also see that at high @xmath196 the effective damping decreases with increasing @xmath196 .    for @xmath198 ,
there is only one crossover .
it separates the regime with two complex conjugate poles from the regime with one pair of complex conjugate poles and two real poles . above @xmath199 ,
the relaxation is governed by the kondo pole .      for a symmetric system ,
two poles are cancelled .
the reduced pole equation reads @xmath200 the characteristic behavior of the poles as function of the scaled temperature @xmath196 is shown in fig .
[ fig : sztz_re_im ] . at low @xmath196 ,
all poles contribute to the dynamics . in the kondo regime @xmath201 ,
the leading pole behaves as @xmath202 , and the amplitudes of the other contributions are negligibly small .    :
crossover temperatures as functions of @xmath9 , @xmath203.,width=6,height=4 ]    : real ( a ) and imaginary ( b ) part of @xmath135 , @xmath204 , @xmath203.,title="fig:",width=6,height=4 ] : real ( a ) and imaginary ( b ) part of @xmath135 , @xmath204 , @xmath203.,title="fig:",width=6,height=4 ]    : plots of @xmath205 in the fourier regime ( a ) and time regime ( b ) , @xmath204 , @xmath203 , @xmath206 .
red ( small - dashed ) @xmath207 , green ( dashed ) @xmath208 , blue ( medium - dashed ) @xmath209 , grey ( long - dashed ) @xmath210.,title="fig:",width=6,height=4 ] : plots of @xmath205 in the fourier regime ( a ) and time regime ( b ) , @xmath204 , @xmath203 , @xmath206 .
red ( small - dashed ) @xmath207 , green ( dashed ) @xmath208 , blue ( medium - dashed ) @xmath209 , grey ( long - dashed ) @xmath210.,title="fig:",width=6,height=4 ]    : real and imaginary part of @xmath135 , @xmath204 , @xmath203.,title="fig:",width=6,height=4 ] : real and imaginary part of @xmath135 , @xmath204 , @xmath203.,title="fig:",width=6,height=4 ]
below the first crossover temperature , @xmath211 , the real parts of all poles are varying linearly with @xmath32 . in detail , the poles and amplitudes are as follows .      there is a superposition of two damped oscillations with frequencies @xmath212 \lambda_{3,4}^ { } & = &   { \displaystyle \pm\,\mathrm{i}\,\bar\omega\,-\,\frac{2\,\bar\omega^2 _ { } - \bar\delta^2 _ { } -   \bar\delta_1 ^ 2}{2\,(\bar\omega^2 - \bar\delta^2)}\,\vartheta_1 \,-\ ,   \frac{\bar\delta_2 ^ 2-\bar\delta^2_{}}{2\,(\bar\omega^2_{}- \bar\delta^2_{})}\,\vartheta_2^ { } \,+\,{\cal o } ( \vartheta^2 _ { } ) \ ; , } \end{array}\ ] ] and amplitudes @xmath213      there is a superposition of two damped oscillations and two contributions describing incoherent relaxation towards the equilibrium value @xmath214 , @xmath215 \lambda_5^ { } & = & { \displaystyle -\,\frac{\bar\omega^2 _ { } - \bar\delta _ 1 ^ 2}{\bar\omega^2 _ { } - \bar\delta^2 _ { } } \vartheta_1^ { } \,-\ ,   \frac{\bar\omega^2 _ { } -\bar\delta_2 ^ 2}{\bar\omega^2 _ { } - \bar\delta^2 _ { } } \vartheta_2^ { }   \ ; , \quad \lambda_6^ { } \;=\ ; -\,\frac{\bar\delta_1 ^ 2 - \bar\delta^2_{}}{\bar\omega^2 _ { } - \bar\delta^2 _ { } } \vartheta_1^ { } \,-\ ,   \frac{\bar\delta_2 ^ 2 - \bar\delta^2_{}}{\bar\omega^2 _ { } - \bar\delta^2 _ { } } \vartheta_2^ { }    \ ; ,   } \end{array}\ ] ] where terms of order @xmath216 are disregarded . in leading order ,
the amplitudes are @xmath217      when the coupling @xmath9 and/or the scaled temperatures @xmath218 are large compared to the other frequencies , the amplitudes of three ( five ) pole contributions to @xmath102 ( @xmath103 ) are negligibly small , and among the dynamical poles only the real pole with smallest modulus is relevant . hence the two spins essentially behave as a single spin which relaxes incoherently to the equilibrium state according to @xmath219 \langle\sigma_z^{}\tau_z^{}\rangle_t & = & { \displaystyle \big(\ , 1\,-\,\langle\sigma_z^{}\tau_z^{}\rangle_{\rm eq}\,\big)\ , \rme^{-\gamma_{\sigma\tau}^{}t } _ { } \ , + \ , \langle\sigma_z^{}\tau_z^{}\rangle_{\rm eq}^ { }   } \ ; . \end{array}\ ] ]      the relaxation rate is found from @xmath220 with the form ( [ sz2 ] ) as @xmath221 this reduces in the parameter regime @xmath222 to @xmath223 in this regime , @xmath180 is independent of the coupling @xmath9 and hence independent of the dynamics of the @xmath13-spin . the temperature dependence
@xmath224 distinguishes the so - called kondo regime , in which , for @xmath225 , the relaxation dynamics slows down as temperature is increased .    on the other hand ,
when @xmath226 , the two spins are locked together @xcite , and the effective tunneling matrix element is @xmath227 , as follows from ( [ def_eigen1 ] ) with ( [ def_eigen2 ] ) .
we then get from eq .
( [ relaxrate ] ) @xmath228 \ ; .\ ] ] this yields the limiting expressions @xmath229 the former is the relaxation rate of the biased single spin - boson system at low @xmath178 .
the latter describes kondo - like joint relaxation of the locked spins .      in the incoherent regime , the relaxation rate @xmath230 of the effective single spin receives rate contributions from both the @xmath12- and the @xmath13-spin as if these were independent biased spins in contact with their own heat reservoir .    in the large - coupling limit , @xmath231
, the relaxation rate is found as @xmath232 the individual contributions are single - spin rates in the large - bias regime . in the high temperature limit @xmath233 ,
on the other hand , both rate contributions are kondo - like , @xmath234    consider next the regime @xmath235 , in which spin @xmath12 behaves kondo - like , as in eq .
( [ kondo1 ] ) .
hence the dynamics of the @xmath12-spin is slow compared to that of the @xmath13-spin .
thus we should expect that @xmath103 approaches the dynamics of the biased single spin - boson case as @xmath178 is increased .
taking into account terms of linear order in @xmath236 in the pole equation , the expression ( [ expecsigtau ] ) with ( [ sigtauz2 ] ) and ( [ sigtauz1 ] ) assumes the form @xmath237 indeed , in the limit @xmath238 , this form reduces just to the analytic expression for @xmath177 of the biased single spin - boson system in the white - noise limit @xcite .
let us now view spin @xmath13 with reservoir 2 as an environment for spin @xmath12 .
this complex environment is in general non - gaussian and non - markovian @xcite .
recently , the same model has been studied numerically using a markovian master equation approach @xcite . to proceed , we first note that in the absence of bath 1 , @xmath239 , the pole equation @xmath240 is still of fourth order .
there is no reduction in the general case . in fig .
[ fig : sbe_re_im ] we show plots of the four poles as functions of @xmath179 for a particular set of parameters .     with spin - boson environment : real ( a ) and imaginary ( b ) part of @xmath135 . at low @xmath179
, @xmath180 is a superposition of two damped oscillations . at high @xmath179 ,
there is one damped oscillation and one relevant relaxation contribution .
the parameters are @xmath241 , @xmath242 , @xmath243.,title="fig:",width=6,height=4 ]   with spin - boson environment : real ( a ) and imaginary ( b ) part of @xmath135 . at low @xmath179 , @xmath180 is a superposition of two damped oscillations . at high
@xmath179 , there is one damped oscillation and one relevant relaxation contribution .
the parameters are @xmath241 , @xmath242 , @xmath243.,title="fig:",width=6,height=4 ]      simplification occurs , however , when @xmath179 is very large compared to the other frequencies . in this regime , the kernel ( [ kplus ] ) reduces to the form @xmath244 where @xmath245 is the relaxation rate of spin @xmath13 in the kondo regime . with this high - temperature expression for the kernel ,
the quantity @xmath107 is found to read @xmath246 this expression describes the dynamics of spin @xmath12 coupled to a spin - boson environment , where the latter is in the kondo regime . to leading order in @xmath247 , the poles of the expression ( [ szsb ] )
are ( see fig .
[ fig : sbe_re_im ] at large @xmath179 ) @xmath248 and the amplitudes read @xmath249 the expressions ( [ szsb ] )  ( [ szsb2 ] ) may now be compared with the corresponding ones of a fictive single biased spin - boson system with parameters @xmath82 and @xmath9 in the white noise limit at scaled temperature @xmath250 .
the part that is symmetric in the bias reads @xmath251 we see that with the identification @xmath252 the expressions ( [ szsb ] ) and ( [ szsb3 ] ) are quite similar .
observe , however , that the damping rate of the oscillation is somewhat different because of the term @xmath253 in eq .
( [ szsb3 ] ) instead of @xmath254 in eq .
( [ szsb ] ) @xcite .
most importantly and interestingly in the correspondence , temperature maps on the inverse of it .
we should expect that , in the weak - coupling limit , the spin - boson environment is gaussian and can be represented by a resonant power spectrum of a bath of bosons .
the gaussian approximation of the spin - boson environment is found by matching the power spectrum of the coupling of the @xmath12-spin to the spin - boson environment @xcite ( with normalization as in eq .
( [ powerspec ] ) ) @xmath255\,\big\rangle\ , \rme^{{\rm i}\,\omega t}_{},\ ] ] with that of a harmonic oscillator bath . in the white - noise limit of an unbiased spin ,
the symmetrized equilibrium correlation function @xmath256 coincides with the expectation @xmath257 .
thus we obtain @xmath258 the resulting power spectrum is that of a structured bath of bosons with a resonance of width @xmath179 at frequency @xmath259 , @xmath260 due to the coupling to the spin - boson environment , the spin @xmath12 performs damped oscillation , @xmath261 . upon calculating the decoherence rate in order @xmath262 , i.e. , the so - called one - boson - exchange contribution of the effective boson bath , we obtain @xmath263 } \ ; .\ ] ] the analysis is completed by observing that this form emerges also upon calculating directly @xmath264 from the pole equation @xmath220 with the form ( [ sz2 ] ) .
we have studied the dynamics of a spin or qubit @xmath12 coupled to another spin , which could be , for instance , another qubit , or a bistable impurity , or a measuring device .
we have solved the dynamics exactly for white - noise reservoir couplings , and we have studied the rich behaviors of the dynamics in diverse limits ranging from weak coupling and/or low temperatures to strong coupling and/or high temperature . we have also analyzed the effects of a spin - boson environment on the spin dynamics in the gaussian and non - gaussian domains .    this paper has not attempted to perform applications to already available experiments , instead we have tried to make some general points on complementary regimes and on the crossovers in between .    one possible simple generalization beyond the white noise limit would be to replace the white - noise bath correlations in time intervals in which the laplace variable is irrelevant by the full quantum noise correlation , for instance @xmath265 the advantage would be two - fold : ( i ) the noise integral is known in analytic form for the ohmic correlation function ( [ ohmcorr ] ) @xcite , and ( ii ) with this substitution the algebraic form of pole equation and residua would be left unchanged
.    one further generalization of the hamiltonian ( [ ham1 ] ) is an applied bias acting on one or both of the spins , e.g. , of the form @xmath266 and @xmath267 .
this important extra ingredient could be taken into account exactly in the white noise regime , and would lead to additional shifts of the laplace variable in all blip states of the @xmath12- and @xmath13-spin . then one would end up with expressions of the form ( [ expecsigma ] ) and ( [ expecsigtau ] ) with polynomials ramped up by bias terms .
this extension will be discussed elsewhere .
finally , extension of the analysis of the dynamics to the regime @xmath268 requires to revert to the original expression ( [ powerspec ] ) and compute its effect perturbatively in the one - boson - exchange approximation .
this can be done either with the self - energy method presented in ref .
@xcite or with the redfield approach .
one then find , e.g. , that the actual equilibrium state of @xmath214 for @xmath269 is @xmath270 \ ; .\ ] ] this reduces for @xmath271 to the previous form ( [ sztzeq ] ) found in the white - noise regime .
the corresponding extension of the analysis and of the results given in subsection [ subseclowtemp ] to the domain @xmath272 will be reported elsewhere .
the authors wish to thank g. falci and e. paladino for valuable discussions . financial support by the dfg through sfb / tr 21 is gratefully acknowledged .