in recent years the spin - boson model @xcite
@xcite has experienced a strong revival , as it is well suited to describe dissipative and decoherence effects on the dynamics of a two - level system ( tls ) or qubit coupled to a bath .
crucial for the effects of the environment on the dynamics of the tls is the shape of the spectral density of the harmonic bath .
it is common to assume an ohmic spectral density , which is linear in the continuous bath modes . in this work we concentrate on a so - called structured bath ,
for which the spectral density is ohmic at low frequencies but exhibits a lorenztian - shaped peak at a certain frequency @xmath0 .
it has been shown in @xcite that a spin - boson model with such an effective spectral density can be exactly mapped on the model of a tls which is coupled to a single harmonic oscillator ( ho ) of frequency @xmath0 , where the latter feels the influence of an ohmic bath.due to its wide applicability the tls - oscillator system has been object of intense research along the years .
so it reflects for example the physics of single atoms with a large electric dipole moment coupled to the microwave photons of a cavity @xcite , or quantum dots in photonic crystals @xcite .
more recently the model has received quite some attention in the field of quantum computation , where two - level systems are used to implement the two logical states of a qubit .
we will especially focus on the solid - state implementation of such systems . here , two prominent realizations of a qubit - oscillator system are the cooper - pair box ( cpb ) @xcite  @xcite coupled to a transmission line resonantor @xcite  @xcite and the josephson flux qubit @xcite read out by a dc - squid @xcite  @xcite .
inspired by experiments with real atoms interacting with a cavity mode , one speaks for the cpb case of circuit quantum electrodynamics , as now the cpb plays the role of an artificial atom and the waveguide acts as a cavity . from such a setup one expects a huge step towards the realization of a quantum computer , as the transmission line resonator can be used to couple qubits together @xcite , store the information of qubits or to provide non - demolition read - out schemes @xcite . concerning the flux qubit
, the read - out usually happens through a damped dc - squid , which is inductively coupled to the qubit .
however , through the squid enviromental noise is transferred to the qubit leading to decoherence and dissipation within its dynamics .
the effect of this noise on the qubit depends very much on the strength @xmath1 of the coupling between qubit and squid and one faces a conflicting situation . on the one hand one
wants a strong coupling for a good read - out resolution . on the other hand
the coupling should be minimized to keep the negative effects of the environment as small as possible . in @xcite
it has been shown that the qubit - squid system can be described by a spin - boson model with an effective spectral density @xmath2 exhibting a peak at the plasma frequency @xmath0 of the squid . applying the above mentioned mapping an equivalent point of view
is to consider the squid as an lc - circuit coupled to the ohmic bath and model it as a harmonic oscillator .
a detailed description of a nondestructive read - out scheme is e.g. given in @xcite .
the spin - boson model can be formally solved using e.g. real - time path integral methods @xcite . however , in order to get closed - form analytical results , approximations must be invoked .
a quite common one is the so - called weak coupling approximation ( wca ) , which is perturbative in the bath spectral density @xcite .
however , it has been shown that for strong qubit - ho coupling @xmath1 and for small detuning @xmath3 , where @xmath4 is the qubit energy splitting , such an approximation breaks down @xcite , as coherent exchange processes between tls and oscillator are disregared .
for an unbiased qubit the non - interacting blip approximation ( niba ) used in @xcite @xcite circumvents this problem as it is non - perturbative in the coupling @xmath1 and therefore takes correctly into account the influence of the oscillator on the tls .
moreover , it allows an analytic treatment of the dynamics .
however , the niba is known to break down for a biased qubit at low temperatures @xcite .
another approach , which treats the system non - perturbatively in the bath is the flow - equation renormalization method @xcite , where the spin - boson hamiltonian is diagonalized using infinitesimal unitary transformations .
however , whithin this approach analytical solutions are difficult to find .
recently a polaron transformation was used by huang to obtain analytically the population dynamics and confirm the shiba s relation for an unbiased tls @xcite.in the case in which the qubit and the ho are considered as the central quantum system being coupled to an ohmic bath , the numerical , ab - initio quasiadiabatic propagator path - integral ( quapi ) method @xcite is a nice tool as it enables to cover both the resonant regime , where the oscillator frequency is close to the qubit energy splitting , and the dispersive regime with the oscillator being far detuned from the qubit @xcite .
moreover , it can be applied to a biased as well as to an unbiased tls and therefore be used as a testbed for analytical results . for qubits
being operated at the degeneracy point , which means an unbiased tls , very often a rotating wave approximation ( rwa ) is applied @xcite , which is expected to be valid for small detuning and yields as starting point the jaynes - cummings hamiltonian @xcite .
this model was first used to study a two - state atom interacting with a single , close to resonance cavity mode of the electromagnetic field and predicts e.g. the repeated revival and collapse of rabi oscillations within the atomic excitation probability . by condsidering the tls - ho system in the representation of displaced ho states , brito were able to truncate the infinite hilbert space of this system without loosing the effects of the ho on the tls dynamics @xcite .
however , so far none of these works could provide an analytical expression for the dynamics of the dissipative qubit being valid for zero as well as non - zero detuning and for both a biased and unbiased tls . in this work an analytic expression for the dissipative qubit s dynamics which includes the effects of a finite detuning and of a static bias is derived . specifically , starting from the qubit - ho perspective , the eigenvalues and eigenfunctions of the non - dissipative tls - ho system are found approximately using van - vleck perturbation theory up to second order in the coupling @xmath1 .
notice that no rotating wave approximation is required .
dissipation effects are then evaluated by solving a born - markov master equation for the reduced density matrix in the system s eigenbasis .
the structure of the work is as follows .
the dissipative tls - ho hamiltonian and the main dynamical quantities are introduced in . inspired by the work of goorden @xcite , we demonstrate in how the eigenstates and eigenenergies of the non - dissipative hamiltonian can be found approximately using van - vleck perturbation theory @xcite . in this way we can provide an analytical formula for the non - dissipative dynamics whichs takes into account the full hilbert space of the qubit - ho system . after that , we show how for low temperatures ( @xmath5 ) this infinite hilbert space can be truncated and discuss the relevant contributions of the ho to the dynamics . in the influence of the environment
is investigated , by looking at solutions of the bloch - redfield equations .
specifically , analytical expressions fo the tls dynamics are obtained and compared with numerical solutions .
the main physical features of the coupled tls - ho system are discussed in . to illustrate the effects of counter - rotating terms in the hamiltonian of the qubit - ho system , which are neglected performing a rwa
, we compare in our calculations to results obtained from the jaynes - cummings model .
in this section we introduce the hamiltonian for a qubit coupled through a harmonic oscillator to a thermal bath . further , a formula for the population difference between the qubit s two logical states is derived .      to set up the model we consider the hamiltonian of a qubit - ho system , @xmath6 , which is coupled to an environmental bath , @xmath7 , by the interaction hamiltonian @xmath8 , so that the total hamiltonian becomes @xmath9 the hamiltonian , @xmath10 consists of @xmath11 the hamiltonian of the tls / qubit and the harmonic oscillator , and the interaction term @xmath12 the hamiltonian of the tls
is given in the subspace @xmath13 , corresponding to a clockwise or counterclockwise current in the superconducting loop of a three - junction josephson qubit or more generally to the qubit s two logical states . in the case of a superconducting flux - qubit ,
the energy bias @xmath14 can be tuned by an applied external flux , @xmath15 , and is zero at the so - called degeneracy point .
the tunnelling amplitude is described by @xmath16 . for @xmath17
the states @xmath18 and @xmath19 are eigenstates of @xmath20 , whereas at the degeneracy point those eigenstates are a symmetric and antisymmetric superposition of the two logical states .
further , @xmath21 and @xmath22 are the annihilation and creation operator for the ho with frequency @xmath0 , and @xmath1 characterizes the coupling strength .
we also introduce the energy splitting @xmath23 between the groundstate @xmath24 and the excited state @xmath25 of the tls . using the transformation @xmath26 with @xmath27 and @xmath28
, we obtain the hamiltonian of the tls in the this basis : @xmath29 .
the states @xmath19 and @xmath18 become in the energy basis @xmath30 the hamiltonian @xmath31 is diagonal in the eigenbasis @xmath32 with @xmath33 being the occupation number : @xmath34 . for the eigenbasis of the combined hamiltonian @xmath35 we write
@xmath36 following caldeira and leggett @xcite , we model the environmental influences originating from the circuitry surrounding the qubit and the oscillator as a bath of harmonic oscillators being coupled bilinearly to the ho .
thus , the environment is described by @xmath37 and the interaction hamiltonian is @xmath38 the operators @xmath39 and @xmath40 are the creation and destruction operator , respectively , for the @xmath41 bath oscillator , @xmath42 is its frequency and @xmath43 gives the coupling strength .
the whole bath can be described by its spectral density , which we consider to be ohmic : @xmath44 in @xcite it is shown that the above model is equivalent to that of a tls being coupled directly to a harmonic bath including the single oscillator of frequency @xmath0 ; i.e. , a spin - boson model @xcite with a peaked effective spectral density , @xmath45 the relation between @xmath46 and the coupling parameter @xmath1 between the qubit and the ho is @xmath47 @xcite .
this second perspective is suitable for calculating the dynamics of the qubit using a path - integral approach , as it was done for example in @xcite for the case of an unbiased qubit ( @xmath48 ) .
the approach in @xcite , however , being based on the niba @xcite , is not suitable to investigate the low temperature dynamics of a biased tls .
thus , in this manuscript we will consider the tls and the single oscillator as central quantum system and solve the bloch - redfield master equations for the density matrix of this system , which are valid also for the case of a biased tls .
the main goal of this work is to determine the dynamics @xmath49 of the qubit .
that means , we wish to calculate the population difference @xmath50 between the @xmath51 and @xmath52 states of the qubit .
the reduced density matrix of the tls , @xmath53 is found after tracing out the oscillator and bath degrees of freedom from the total density matrix @xmath54 . in turn
@xmath55 is the reduced density matrix of the qubit - ho system .
how to calculate this density matrix will be shown later .
after some algebra , illustrated in more detail in [ apppt ] , we arrive at an expression for @xmath49 , given in terms of diagonal and off - diagonal elements of @xmath56 in the tls - ho eigenbasis @xmath57 .
it reads @xmath58 where @xmath59   + 2 \sin \theta { \langle j { \rm g}|n\rangle } { \langle j { \rm e}|n\rangle } \right\ } \rho_{n n}(t ) \label{p0},\\    \fl p_{nm}(t ) = 2 \sum_j \biggl\ { \cos \theta \biggl[{\langle j { \rm g}|n\rangle } { \langle m|j { \rm g}\rangle } - { \langle j { \rm e}|n\rangle } { \langle m|j { \rm e}\rangle }     \biggr ] \nonumber\\       + \sin \theta \biggl[{\langle j { \rm e}|n\rangle } { \langle m|j { \rm g}\rangle } + { \langle j { \rm e}|m\rangle } { \langle n|j { \rm g}\rangle }     \biggr ] \biggr\ } \re \ { \rho_{nm}(t ) \ } \label{pnm}\end{aligned}\ ] ] with @xmath60 . how to determine the eigenstates of @xmath6 is described in the next section .
in this section we show how to find the eigenvalues of the unperturbed qubit - ho hamiltonian @xmath6 approximately by using van - vleck perturbation theory @xcite .
the idea is to take advantage of the degenerate or doublet structure of the energy spectrum of the uncoupled ( @xmath61 ) tls - ho system near resonance , e.g. at @xmath62 .
then , as long as the perturbation is small compared to the energy separation of the different doublets , the full hamiltonian will exhibit a similar spectrum of bundled energy levels .
the eigenenergies of the uncoupled tls - ho system are immediately found by applying the hamiltonian @xmath63 on the eigenstates in : @xmath64 the dashed lines in show the energy spectrum corresponding to vs. the oscillator frequency @xmath0 for the five lowest eigenstates . except for
the groundstate , @xmath65 , the states @xmath66 and @xmath67 are degenerate in the resonant case @xmath68 .
close to resonance the spectrum exhibits a doublet structure . with the coupling
being switched on , the full hamiltonian @xmath69 reads @xmath70 in the basis @xmath71 . in order to diagonalize the hamiltonian @xmath72 we consider @xmath73 as a small perturbation , which is resonable as long as @xmath74 . applying van - vleck perturbation theory we construct an effective hamiltonian , @xmath75 having the same eigenvalues as @xmath76 but no matrix elements connecting states which are far off from degeneracy . thus , @xmath77 will be block - diagonal with all quasi - degenerate energy levels being in one common block . as in our case
always two states are nearly degenerate , each block of @xmath78 builds a two - by - two matrix .
this can be easily diagonalized in order to determine the eigenstates .
following @xcite we calculate the transformation matrix @xmath79 up to second order in @xmath1 .
the general formulas for both an arbitrary hamiltonian and @xmath76 are given in [ appvanvleck ] .
the only surviving matrix elements of the effective hamiltonian , apart from the ones being of zeroth order in @xmath1 , are @xmath80 and @xmath81 . \label{heff2}\end{aligned}\ ] ] thus , @xmath82 has the matrix structure @xmath83 where the section shown corresponds to the basis states @xmath84 and @xmath85 . from this form
it is easy to calculate the eigenstates and eigenenergies .
the groundstate @xmath86 , which is an eigenstate of @xmath87 , has the eigenenergy @xmath88 the other eigenstates of @xmath87 are , @xmath89 , @xmath90 corresponding to the eigenenergies @xmath91\\                               = \hbar \left [ ( j+\frac{1}{2 } ) \omega + w_1 + w_0 \mp   \frac{1}{2 } \sqrt{\delta_j^2 + 4(j+1 ) |\delta|^2 } \right ] , } \ ] ] @xmath92 by construction these are also eigenenergies of @xmath93 . using the transformation ( [ vanvlecktrafo ] ) we get the eigenvectors of @xmath93 as @xmath94 the energy spectrum of @xmath93 is shown in for the case of an unbiased tls ( @xmath95 ) .
we want to emphasize that our findings are also valid for the more general case @xmath96 . at resonance , where the spectrum for the uncoupled case is degenerated ,
avoided crossings can be seen .
the gap between two formerly degenerated levels for @xmath97 is @xmath98 which is as predicted by the jaynes - cumming model @xcite .    as we will show in ,
the second order correction @xmath99 in , whichleads to a shift in the resonance frequency , is a result of the counter - rotating terms in @xmath93 .
as such it can be interpreted as a bloch - siegert shift @xcite .      with the coupling to the bath being turned off ,
the time evolution of the density matrix of the qubit - ho system is given by @xmath100 and consequently @xmath101 with @xmath102 .
with that becomes @xmath103 where we defined @xmath104 . from we notice that the dynamics of the qubit is characterized by an infinite number of oscillation frequencies rather than showing rabi oscillations with a single distinct frequency .
this is clearly a consequence of the coupling of the ho to the tls .
further we assume that at @xmath105 the qubit starts in the state @xmath51 and that the occupation numbers of the ho are boltzmann distributed , so that @xmath106 where @xmath107 is the partition function of the oscillator and @xmath108 denotes the inverse temperature of the system . in the tls - ho eigenbasis
this becomes @xmath109\\                \times   \left [ \cos \left ( \frac{\theta}{2 } \right ) { \langle j{\rm g}|m\rangle } + \sin \left ( \frac{\theta}{2 } \right ) { \langle j { \rm e}|m\rangle } \right].}\ ] ]      with we found a formula which describes using the approximate eigenenergies and eigenstates in and the non - dissipative dynamics up to second order in @xmath1 , thereby taking into account all oscillator levels .
thus , we still have to deal with an infinite hilbert space .
typically experiments , see e.g. in @xcite , run in a temperature regime for which @xmath110 . considering the exponential function in we assume the higher oscillator levels to be only sparsely populated and the maximum value of the sum in is truncated to @xmath111 . nevertheless , states @xmath112 with @xmath113 still play a role in the dynamics . in fact , due to the van - vleck transformation @xmath114 , for example the state @xmath115\ ] ] yields nonvanishing contributions to the matrix elements @xmath116 and @xmath117 occurring in due to the fact that the energy eigenstates of the coupled tls - ho system are made of linear combinations which involve also these states .
using together with one finds that coefficients @xmath118 with @xmath119 are of higher than second order in @xmath1 .
the same is valid for @xmath120 , @xmath121 , @xmath122 and @xmath123 .
thus , those terms play no role in our calculation of @xmath49 . furthermore
, @xmath124 neglecting also these contributions we find that @xmath125 for @xmath126 . in the end it will be sufficient to concentrate on eigenstates of @xmath93 up to @xmath127 .
this trunctation leaves us with ten possible oscillation frequencies @xmath128 , where @xmath129 and @xmath130 .    as an example
we calculate the dynamics of an unbiased tls ( @xmath95 ) . here
the coefficients @xmath131 , @xmath132 , @xmath133 , @xmath134 and @xmath135 vanish due to symmetry , so that @xmath136 additionally as a benchmark we consider the mostly studied resonant case , where @xmath137 .
in this case we find with the transitions frequencies @xmath138 the dynamical quantity @xmath49 and its fourier transform are shown in .
one clearly sees the influence of the coupling to the ho on the dynamics of the tls . instead of rabi oscillations with a single frequency
, @xmath49 oscillates with six different frequencies , which are in the fourier spectrum symmetrically located around the point @xmath139 . among those frequencies @xmath140 and @xmath141 are dominating .
they correspond to transitions between the first or second energy level of the qubit - ho system and its groundstate and their weight is almost equal . to summarize , one notices that due to the coupling with the oscillator additional frequencies are induced into the qubit dynamics .
theoretically , the number of those frequencies is infinite . at low temperatures , however , transitions between the lower energy levels of the system are clearly dominating . again , for simplicity we have shown here the case of an unbiased tls . for @xmath96
the behaviour is similar only that in the fourier spectrum the weight difference of the two dominating peaks will be more pronounced .
in the preceding section we neglected the influence of the bath on the qubit - ho system . yet , in order to model a realistic situation , we have to pay attention to environmental influences , as they lead to decoherence and dissipation in the dynamics of the qubit , which is harmful for quantum computing application .
thus , we will now consider the full hamiltonian @xmath142 .
as shown in , we need for the calculation of the population difference @xmath49 the density matrix @xmath56 of the qubit - ho system . starting from the liouville equation of motion for the full density matrix @xmath143 of @xmath142 , @xmath144,\ ] ] where the index stands for the interaction picture and following @xcite
, we can provide a born - markov master equation for @xmath56 being in the schrdinger picture and expressed in the basis of the eigenstates of @xmath93 : @xmath145 the free dynamics of the system is given by the first term of the right - hand side in the above equation .
the rate coefficients are defined as @xmath146 x_{nk } x_{lm } \\
-\delta_{ml } \sum_{l^\prime } g(\omega_{l^\prime k } ) n_{l^\prime k } x_{n l^\prime } x_{l^\prime k }   + \delta_{nk } \sum_{k^\prime } g(\omega_{l k^\prime } ) n_{k^\prime
l } x_{l k^\prime } x_{k^\prime m}}\ ] ] with @xmath147 $ ] and @xmath148 . for the derivation of the master equation , besides the born and markov approximations , some more assumptions have been made , which we briefly mention .
first , we consider our system and the bath to be initially ( at @xmath105 ) uncorrelated ; i.e. , @xmath149 with @xmath150 and @xmath151 the partition function of the bath .
further , with the bath consisting of infinite degrees of freedom , we assume the effects of the interaction with the qubit - ho system to dissipate away quickly , such that the bath remains in thermal equilibrium for all times @xmath152 : @xmath153 besides , an initial slip term which occurs due to the sudden coupling of the system to the bath is neglected @xcite . and
as last approximation the lamb - shift of the oscillation frequencies @xmath128 was not taken into account @xcite .      in @xmath154 describes matrix elements of the operator @xmath155 in the qubit - ho eigenbasis . by use of , and those states were expressed in the basis @xmath156 , and we will also calculate the oscillator matrix elements in this basis . for that purpose the operator @xmath157 is defined .
four different situations can be distinguished .
there are matrix elements were neither the oscillator nor the qubit state are changed , namely @xmath158 and @xmath159 with @xmath160 .
we see that those elements are independent of @xmath161 , the occupation number of the oscillator .
next , we look at the case where a single quantum is emitted or absorbed from the oscillator and get @xmath162 and @xmath163 with @xmath164 for a transition within the qubit we have @xmath165 . and finally , if the qubit and the oscillator state are changed simultaneously , one obtains @xmath166 and @xmath167 , @xmath168 comparing the magnitude of the transition terms , we notice that those consisting in changes of the oscillator occupation only are the dominant ones , as they have a part which is of zeroth order in @xmath1 . further , for the case in which the qubit is operated at the degeneracy point @xmath169 and @xmath170 vanish . with those results we can calculate the matrix elements @xmath154 .
they are given in [ apposcillatormatrixelements ] .
like in we assume the system to be operated at low temperatures and thus take as highest qubit - ho state the eigenstate @xmath127 . for determination of @xmath49 the formulas of can be used .
unlike in the non - dissipative case @xmath56 is not given anymore by the simple expression .
rather we have to solve a system of coupled differential equations , namely . to do this analytically we will follow three different approaches and compare them finally to the numerical solution of .
we start by introducing @xmath171 which yields the set of differential equations for @xmath172 @xmath173      as a first approach we make the full secular approximation ; i.e. , we neglect fast rotating terms in and keep only contributions where @xmath174 vanishes . in this way the off - diagonal elements of @xmath175 are decoupled from the diagonal ones so that @xmath176 the equation for the off - diagonal elements is then @xmath177 which becomes with @xmath178 as through the fsa the oscillatory motion of the dynamics is separated from the relaxation one we can divide into two parts , @xmath179 where @xmath180 describes the relaxation and @xmath181 the dephasing parts of the dynamics . with the latter takes the form @xmath182 with the dephasing rates @xmath183 .
expressions for the dephasing coefficients @xmath184 can be found in [ ratecoeff ] and the initial conditions @xmath185 are given by .
the diagonal elements are more difficult to obtain , as one has to solve a system of coupled differential equations , . calculating the corresponding rate coefficients of this system for the five lowest eigenstates ,
we find that there are only eight independent ones , namely @xmath186 , @xmath187 , @xmath188 , @xmath189 , @xmath190 , @xmath191 , @xmath192 and @xmath193 .
they are given by @xmath194 where @xmath161 and @xmath195 adopt the above values .
furthermore , @xmath196 , @xmath197 , @xmath198 and @xmath199 vanish .
the remaining rate coefficients are combinations of the above .
we find that @xmath200 and @xmath201    however , the system is still too complicated to be solved analytically .
thus , we invoke a further approximation : we consider the factor @xmath202 $ ] with @xmath203 in and use that @xmath204 .
it depends strongly on the temperature @xmath205 for which value of @xmath128 this limit is reached approximately .
for the parameters we are working with one usually is in the region where @xmath206 .
thus , we will neglect in the following terms containing the factor @xmath207 .
furthermore , one sees from that @xmath208 and @xmath209 . with
that @xmath210 and @xmath211 can be neglected . using
the matrix of the system ( [ diagonalfsa ] ) becomes @xmath212 the eigenvalues and eigenvectors of this matrix and the associated time evolution of the elements @xmath213 are given in  of [ appdiagonalelements ] .
unlike for the dephasing part , we can not extract a simple analytical expression for the relaxation rate as @xmath214 now consists of a sum of several exponential functions , cf . together with  .
but still we are able to provide an analytical formula for @xmath49 using .      in order to obtain a simple expression for the relaxation part
, we consider the long - time dynamics of the system . in other words , rather than looking at the many relaxation contributions to the populations @xmath213 , we focus on the smallest eigenvalue of the relaxation coefficients , as it will dominate at long times .
further , we consider only the rate matrix associated to the three lowest qubit - ho eigenstates , @xmath215 , @xmath216 and @xmath217 in and obtain with - that @xmath218 here , we have not neglected the rate coefficients containing the term @xmath207 and further took @xmath188 into account despite of being of third order in @xmath1 as such contribution removes the degeneracy between the two lowest eigenvalues at resonance , cf .
the smallest eigenvalue reads @xmath219^{1/2 } \biggr\}.\end{aligned}\ ] ] with the system being detuned this expression can be simplified further , namely @xmath220    in the relaxation rate @xmath221 as obtained from is shown as a function of the oscillator frequency @xmath0 .
clearly , it is maximal at resonance ( @xmath222 ) , whereas it decays with @xmath0 being detuned from the resonance .
this effect has already been predicted by blais @xcite .
as the qubit is not directly coupled to the bath but rather through the oscillator , the latter being detuned filters out the environmental noise at the qubit transition frequency .
additionally , we show the second smallest eigenvalue of . we notice that close to the resonant point ( @xmath222 ) there is an avoided crossing .
finally , we find that @xmath223 where like in @xmath104 . for getting @xmath224 we have in principle to find the steady - state solution of .
here , we just assume for @xmath225 a boltzmann distribution for the qubit - ho system , so that @xmath226 with @xmath227 .
thus , @xmath228
+ 2 \sin \theta { \langle j { \rm g}|n\rangle } { \langle j { \rm e}|n\rangle } \right\ } \rho_{n n}(\infty ) \label{pinf}.\ ] ] the formula for the long - time dynamics is obtained as @xmath229 to get further insight on the dominant frequencies we evaluate the fourier transform of according to @xmath230 yielding @xmath231.\end{aligned}\ ] ]      an improvement to the fsa is to take into account certain non - vanishing contributions of @xmath174 .
we have to keep in mind , that there are quasi - degenerate levels close to resonance . in our case
the first with second energy level and the third with fourth one build a doublet , meaning that they are close together in energy space .
the level spacing is approximately proportional to @xmath1 for the former and @xmath232 for the latter . because of that and as the transitions from level three and four are less probable , we will in the following only consider the first and second level as being almost degenerate . taking this into account in
we arrive for the diagonal elements at @xmath233 a numerical analysis shows that the effect of the last two terms on the right - hand side of the above equation will in the worst case lead to very small wiggles in @xmath213 and play no role in @xmath49 .
thus , we finally write @xmath234 which is the same equation as we got in the fsa approach . however , the off - diagonal contributions @xmath235 , @xmath236 , @xmath237 , @xmath238 , @xmath239 and @xmath240 have to be examined more carefully . from ( [ masterwithansatz ] )
we find that one has to solve the equations @xmath241 with @xmath242,@xmath243 or @xmath244 as solution one gets @xmath245 here , the oscillation frequencies and the decay of the off - diagonal elements are given by @xmath246\ ] ] with @xmath247 ^ 2 + 4 \pi^2 \mathcal{l}_{nm , jk } \mathcal{l}_{jk , nm}}.\ ] ] the amplitudes of the oscillations are given through the coefficients @xmath248}{2 r_{nm , jk}}\ ] ] and @xmath249.\ ] ]    thus , we have again all ingredients to calculate analytically the relaxation and dephasing part of . for the psa we can not provide a simple expression for the dephasing rates as in the fsa , where we had @xmath250 . as one can see from and ,
@xmath251 and @xmath252 are a mixture of contributions decaying with @xmath253 and @xmath254 .
similar to our findings for the relaxation rate , also here the smallest eigenvalue will dominate the dephasing behaviour . from the right graph in
we find that this is @xmath255 .
comparing it with the dephasing rates we got using the fsa , left graph in , we see that for negative detuning ( @xmath256 ) the rate @xmath257 approximates @xmath258 , whereas for positive detuning ( @xmath259 ) this is done by @xmath260 . in the fsa @xmath261 and
@xmath262 correspond to the frequencies @xmath140 and @xmath141 , respectively . in the psa
the frequency @xmath263 is given by @xmath264 for @xmath256 and @xmath265 for @xmath259 .
hence , for negative detuning oscillations with frequency @xmath141 will dominate the dynamics , while those with @xmath140 will almost vanish . for positive detuning
it is the other way round . in and
this behaviour is reflected by the coefficients @xmath266 and @xmath267 . around resonance ( @xmath268 )
the psa tells us by and that the dephasing rates and frequencies are a mixture of @xmath262 and @xmath261 or @xmath140 and @xmath141 , respectively . from the left graph in one
notices further that the fsa rate @xmath261 grows linearly with @xmath0 for positive detuning .
however , as the weight of the corresponding frequency @xmath141 will be almost zero , @xmath261 will give no relevant contribution to @xmath269 in this regime but the dephasing will rather be associated to the fsa rate @xmath262 .
hence , out of resonance the fsa will still fairly well describe the dynamics of @xmath49 . comparing the expressions for @xmath270 and @xmath271 given in [ ratecoeff ] by and with the approximative expressions for the relaxation rate at positive and negative detuning
, we see that for zero bias ( @xmath272 ) the psa dephasing rate is equal to @xmath273 . for a biased system an additional term
is added depending on the spectral density of the bath at @xmath274 .
in we compare the three analytical solutions described above to the numerical solution of the master equation for the case of an unbiased tls being at resonance with the oscillator . concerning both the dynamics of @xmath49 and its fourier spectrum we see a good agreement between the different solutions .
the one being closest to the numerical solution is the psa solution .
we also want to mention that going to stronger damping @xmath275 , the fsa results start to show deviations from the numerical solution . here , one should use the psa only .
however , for the parameter regime used in the following , we will mainly apply due to its simple , analytical form .
having solved the master equation analytically and numerically we can examine the dynamics of the system and its fourier transform for different situations .
first , we will look at a qubit operated at the degeneracy point ( @xmath48 ) being in and out of resonance with the oscillator .
then , we will concentrate on the biased qubit in the same regime of parameters .      for unbiased qubits
we can compare our predictions with the analytical results obtained in @xcite by starting from a spin - boson model with the effective spectral density . in @xcite
a so - called weak damping approximation ( wda ) based on the non - interacting blip approximation ( niba ) is applied .
the wda allows a non - perturbative treatment of the coupling between the tls and ho and hence can reproduce the occurence of two dominating frequencies as expected e.g. from exact quapi calculations @xcite .
the niba , and hence the wda , however , become not reliable for a biased tls .
we find that the overall agreement between our approach and the wda is very good .
however , in the wda solution the frequencies are slightly shifted compared to the ones obtained from our master equation .
this may result from the perturbative expansion we have performed with respect to @xmath1 by applying the van - vleck perturbation theory .
first , we look at the resonant case shown in .    in agreement with previous works @xcite , we find that the dynamics is dominated by two frequencies corresponding to @xmath140 and @xmath141 with separation being approximately @xmath276 .
the weight of the latter is a bit larger .
the reason for the bigger weight is that at resonance ( @xmath222 ) the qubit - ho eigenstate @xmath277 is not a symmetric superposition of the states @xmath278 and @xmath279 unlike it is predicted by the janyes - cummings model ( cf e.g. @xcite ) .
we notice that the two unequal peaks have indeed been experimentally observed in @xcite ( see fig .
4b therein ) . considering the states
@xmath280 and @xmath281 in and , one already sees that for a symmetric superposition of these states we need that @xmath282 vanishes or that @xmath283^{-\frac{1}{2}}$ ] ( cf . ) . besides , in order to get the qubit - ho eigenstates one still has to perform the van - vleck transformation , which adds contributions to @xmath216 and @xmath217 from states corresponding to oscillator levels higher than @xmath111 .
thus , our system behaves for @xmath284 as being negatively detuned , which means that the peak belonging to the higher frequency dominates , as we will show below .
slightly increasing @xmath0 will give a stronger weight to the peak at @xmath140 .
this effect is not very pronounced for the non - dissipative dynamics of the unbiased qubit ( ) , as there the two frequencies are still almost equally weighted .
looking however at the fourier transform of the dissipative dynamics , one notices that the relaxation rate also contributes to the weight of the peaks with a prefactor @xmath285 . as for a negative detuned system
@xmath262 is slightly bigger than @xmath261 , the difference between the two peaks becomes more clear in the dissipative case . for @xmath96 the effect can already be noticed in the non - dissipative case.next , we consider in the case of negative detuning , where @xmath286 .
no matter which approach one is looking at , clearly the frequency @xmath141 is dominating .
furthermore , paying attention to the timescale of the dynamics , one notices that the relaxation time is _ enhanced _ compared to the one we found for the resonant system .
this behaviour was already explained by the formula for the relaxation rate .
again , the numerical and the solution obtained by using the long - time ansatz in agree quite well with each other , whereas the amplitude of the oscillation with frequency @xmath141 is stronger in the wda approach .    also remarkable is the fact that looking at the fourier transform in one sees in the inset already small contributions of the higher oscillator levels .
the transitions corresponding to @xmath287 and @xmath288 give raise to small additional peaks , while the contributions of @xmath289 and @xmath290 are negatively weighted and cause dips .
the wda approach does not show this additional contributions .
they are , however , confirmed by the numerical quapi calculations in @xcite ( see figure 2 therein ) .    in the case of positive detuning ( @xmath291 ) shown in
we find a quite good agreement between all three approaches . also for postive detuning the relaxation time is enhanced compared to the resonant case .
contrary to the negatively detuned situation the additional peaks have vanished .
besides , now the frequency @xmath140 is dominating the dynamics .
this behaviour , namely that for negative detuning @xmath141 and for positive detuning @xmath140 is dominating , was already found in @xcite .
we will briefly explain how one can explain this observation physically .    for this
we look at . for a detuned system ( @xmath292 ) the qubit - ho eigenstates are not symmetric superpositions of the states @xmath293 and @xmath294 .
they rather asymptotically approach the eigenstates of the uncoupled qubit - ho hamiltonian . in
we see that for a negatively detuned system ( line a ) the qubit - ho eigenstate @xmath295 approaches the state @xmath296 , whereas the main contribution to the state @xmath297 will come from the state @xmath298 . from the left diagram in
we see that the state @xmath217 is energetically higher than the state @xmath216 . however , due to the boltzmann distributed occupation of the oscillator , the state @xmath299 will be more populated than the state @xmath300 and consequently also @xmath217 will exhibit a larger population than @xmath216 , as the latter only feels a small contribution from the state @xmath299 .
thus , transitions from @xmath217 to the groundstate are more likely to occur than those from @xmath216 to the groundstate .
this explains the dominance of @xmath141 in and . in this case
the frequency @xmath301 and @xmath302 . as far as not excluded by selection rules , minor peaks from transitions to the levels lying in between can be also seen.for positive detuning ( line c in ) @xmath295 approaches @xmath298 , while @xmath297 is close to @xmath296 . from the right graph in
we see that the state @xmath216 , being lowest in energy apart from the groundstate , is now also more probable to be occupied than @xmath217 . therefore , as confirmed by and , the frequency @xmath140 is dominating whereas @xmath141 is represented only by a small peak in the fourier spectrum
. furthermore , as there are no additional energy levels between the state @xmath216 , which is most probably to be populated , and the ground level , other transitions than those corresponding to @xmath140 or @xmath141 are very unlikely to occur . in the dip corresponding to @xmath303 appears only very faintly .
we will now examine a qubit being operated at finite bias .
we consider the case @xmath304 . for negative
bias - offset the behaviour is analogous .
again three different situations are taken into account : the qubit being in resonance with the oscillator ( @xmath305 ) , negative ( @xmath306 ) and positive ( @xmath307 ) detuning.for the resonant case ( @xmath222 ) depicted in we see a similiar behaviour as for the unbiased qubit . again
two frequencies , @xmath140 and @xmath141 , are dominating the dynamics .
left to the peak at @xmath140 a small dip can be found in the fourier spectrum .
this corresponds to the transition @xmath303 .
for infinite time the dynamics relaxes to an equilibrium value which is nonzero in contrast to the unbiased case .
this can be seen in the fourier spectrum through a relaxation peak at @xmath308 .
the peak arises because of the term @xmath309 in .
the first part of this sum gives rise to the negative shift of this peak .
the reason that for the analytical solution the peak is not as strongly shifted as for the numerical one is technical : in order to plot the delta function in we gave it a finite width , which surpresses the negative contribution of the first term in . like for the unbiased qubit
the highest energy level playing a role for the dynamics is @xmath310 ; i.e. , only the ground and first excited level of the oscillator are of importance .    in the dynamics and its fourier transform for a negatively detuned qubit - oscillator system with @xmath96 are shown .
like for the unbiased case detuning gives raise to longer relaxation times for the qubit .
also in agreement with the unbiased case is the dominance of the frequency @xmath141 .
we see that for small @xmath152 the long - time solution slightly overestimates the maxima of the oscillations and underestimats its minima .
furthermore , we get here the unphysical situation that the maximum of the third oscillation in @xmath49 exceeds the value of one .
the reason for that behaviour is that , by construction , we underestimate with the relaxation at short times . as for certain paramteres
the term @xmath311 in can become negative , it increases too fast towards the equilibrium and gives thus raise to the observed deviations in the short time behaviour . on a longer timescale both graphs agree quite well .    for the case of positive detuning ( @xmath259 ) , which is presented in ,
the upward shift of the dynamics obtained from and compared to the numerical graph of @xmath49 at small times is even stronger . to visualize that it is not a failure of the fsa approach we show in and additionally the analytical fsa solution of and calculated in .
the latter agrees very well with the numerical solution . at long time - scales and for the fourier spectrum
all three approaches match with each other very well.to conclude this paragraph we want to mention that all the results found both for the unbiased and the biased qubit confirm the numerical quapi results in @xcite .
so far we have always considered the qubit for certain values of @xmath14 and finite or zero detuning . in this section ,
we fix the oscillator frequency at @xmath284 .
that means that an unbiased qubit will be at resonance with the oscillator .
changing the bias to positive or negative values will always lead to negative detuning , as @xmath312 .
shows a density plot of the fourier transform of the symmetrized correlation function against the bias of the qubit and the fourier frequency @xmath313 .
we consider this correlation function rather than @xmath49 , as it is symmetric in the bias @xmath14 .
the symmetrized correlation function is defined as follows @xcite : @xmath314 where @xmath315 .
expressed in terms of the population difference @xmath49 this becomes , @xmath316 with @xmath317 and @xmath318 being symmetric and antisymmetric in @xmath14 and @xmath319 .
the fourier transform of @xmath320 is defined as @xmath321 considering now we see that for any bias the spectrum is dominated by two frequencies , namely @xmath140 and @xmath141 .
detuning the system @xmath141 gets more and more important , as we could already observe in the two previous sections for the positively detuned systems .
furthermore , the peaks are shifted to higher frequency values and at @xmath322 the relaxation peak occurs .
we want to compare these results to a circuit qed experiment performed by wallraff @xcite .
there the qubit is realized by a cooper pair box , which is coupled to a superconducting transmission line resonator .
the properties of the system are determined by probing the resonator spectroscopically .
the amplitude of a microwave probe beam transmitted through the resonator is measured versus the probe frequency and the gate charge of the cooper pair box ( see figure 4 in @xcite ) . via the gate charge
the qubit can be detuned in situ from the degeneracy point .
the frequency of the resonator is chosen in such a way that it is in resonance with a qubit being operated at the degeneracy point .
for the resonant case two dominating frequencies , being almost equally weighted and symmetrically positioned around the cavity frequency , are observed .
going away from the degeneracy point the system becomes detuned and the frequency of the cavity dominates .
the behaviour we observe in is similar .
however , as we are looking at the dynamics of the qubit , it corresponds to a spectroscopic measurement on the tls rather than on the oscillator .
as explained above the two lowest excited states of the coupled tls - ho system , namely @xmath216 and @xmath217 , evolve from an almost symmetric superposition of basis states @xmath323 at resonance @xmath324 to the states @xmath325 and @xmath326 ( cf the left graph in ) . for @xmath327 the two peaks of the rabi splitting
are observed . for @xmath96 , which means negative detuning in this case , the peak with the lower frequency corresponding to @xmath140 approaches more and more the frequency @xmath0 of the oscillator , as the state @xmath216 becomes @xmath325 for large detuning and then @xmath328 .
furthermore , the transition peak at @xmath140 gets weaker as also the occupation probability of @xmath216 decreases . at @xmath329
the symmetrized correlation function vanishes at @xmath140 and increases again for higher values of @xmath330 . here , the amplitude @xmath331 in @xmath49 changes its sign .
in contrast the peak at @xmath141 becomes stronger with the detuning and approaches more and more the qubit splitting energy @xmath332 , as @xmath217 approaches @xmath299 and then @xmath333 .
additionally , looking at the logarithmic plot one sees around @xmath334 a peak appearing , which corresponds to the frequency @xmath303 and is forbidden at @xmath95 . for large detuning
it arises from transitions from @xmath299 to @xmath325 and therefore has the value @xmath335 .
the amplitude of this peak is very small compared to the peaks at @xmath140 and @xmath141 and is not resolved in the experiment of wallraff .
van - vleck perturbation theory enabled us to find approximately the eigenstates and eigenenergies of the full hamiltonian of the tls - ho system without performing a rotating - wave approximation . using those eigenstates and eigenenergies in a born - markov master equation we could calculate the dynamics of such a system under the influence of an environmental bath . in the following
we will show how the results change if we neglect counter - rotating terms in the tls - ho hamiltonian for @xmath95 .
for this we rewrite the interaction part in as @xmath336 where we identified with @xmath337 and @xmath338 a rotating and counter - rotating part of @xmath339 , respectively , and introduced the two - level transition operators @xmath340 . neglecting the counter - rotating part @xmath338 in @xmath93 leads to the jaynes - cummings hamiltonian @xmath341
this hamiltonian can be diagonalized exactly and its eigenstate and eigenvalues can for example be found in @xcite . in order to see the effect of not taking into account the counter - rotating terms we diagonalize @xmath342 using van - vleck perturbation theory .
looking at the formula fo the effective hamiltonian in [ appvanvleck ] and keeping in mind that we set @xmath14 to zero , we see that the second order contributions in @xmath1 vanish neglecting @xmath338 ; i.e. , @xmath343 , @xmath344 are zero in . considering further the transformation matrix @xmath79 , equations and in [ appvanvleck ] show that @xmath345 for @xmath346 .
thus , with we find that _ for the jaynes - cummings model @xmath347 is identical to @xmath342 _ and therefore the eigenstates of @xmath347 are simultaneously eigenstates of @xmath342 .
consequently , one can determine from  the eigenstates and eigenenergies of @xmath342 .
the energy of the groundstate @xmath348 is @xmath349 . for the higher states we get @xmath350 corresponding to the eigenenergies
@xmath351                               = \hbar \left [ ( j+\frac{1}{2 } ) \omega   \mp   \frac{1}{2 } \sqrt{\delta_{\rm jc}^2 + 4(j+1 ) |\delta|^2 } \right],\ ] ] with @xmath352 and @xmath353 . comparing these eigenstates and eigenenergies to the ones found for @xmath93 , , and , we see that the counter - rotating terms yield second order corrections in @xmath1 not present in the jaynes - cummings hamiltonian .
these corrections give rise to a very prominent effect concerning the resonance condition between tls and ho .
from @xmath354 we find the tls being in resonance with the oscillator for @xmath97 .
considering @xmath355 in this resonance condition is shifted to @xmath356.\ ] ] this second order correction to the resonance frequency due to counter - rotating terms is known as bloch - siegert shift @xcite . the eigenstates and are always a superposition of two basis states of the unperturbed system .
this is like for the eigenstates and of the effective hamiltonian .
however , in order to find the eigenstates @xmath357 of @xmath358 we had to apply the transformation @xmath359 on the effective eigenstates so that @xmath357 is in the end a superpostion of _ several _ states of the basis @xmath360.to calculate the reduced density matrix of the qubit - ho system described by @xmath342 and taking into account the influence of an environmental oscillator bath we can again use the born - markov master equation .
we just have to use and as basis states and the corresponding eigenenergies . for the population difference
@xmath49 we apply , which becomes for @xmath48 @xmath361   \re \ { \rho_{nm}(t ) \}.\end{aligned}\ ] ]      before doing a qualitative comparison between the results obtained from the original hamiltonian @xmath93 and the simplified form @xmath342 , we want to analyse which transitions between the different eigenstates of @xmath342 yield contributions to @xmath49 . for this
it is helpful to rewrite and . in the following we neglect the upper index jc denoting an eigenstate of @xmath342 . for the state @xmath357
we have three possibilities : the first one corresponds to @xmath362 . in this case
the only non - vanishing component of @xmath215 is @xmath363 .
second , @xmath364 can be an even number , which means , expressed in terms of oscillator quanta , that @xmath365 . then @xmath366 and third for an odd state @xmath367
we find @xmath368 it is quite easy to see that the part @xmath369 in vanishes for any @xmath364 .
that means that the diagonal elements of the reduced density matrix yield no contributions to @xmath49 for @xmath48 and the equilibrium value of the dynamics will be also zero .
not so easy to see are the combinations of @xmath364 and @xmath370 yielding contributions to the off - diagonal part @xmath371 . considering transitions to the groundstate ( @xmath372 ) , we find using that @xmath373 that @xmath374 with ( [ evenstate ] ) and ( [ oddstate ] ) these elements are non - zero only if @xmath375 or @xmath376 . for the more general case with @xmath377 and
both being different from zero we have to investigate products of components like @xmath378 .
for @xmath379 we find that @xmath380 and that @xmath381 so that @xmath382 for the case of transitions between odd states , one gets that @xmath383 and similarly under exchange of @xmath384 and @xmath385 , so that @xmath386 for transitions from an even to an odd state and vice versa we have to pay attention to the fact that @xmath387 which yields the selection rule @xmath388 .
further , @xmath389 yielding @xmath390 . to sum up : @xmath371 is non - zero if one of the three following cases is valid : @xmath391 , @xmath392 , @xmath388 or @xmath390 .
the allowed transitions are shown in figure [ figselecrules ] .
one sees that transitions between almost degenerate levels are forbidden .
this behaviour we have also found in for the non - dissipative dynamics resulting from @xmath93 .      in the following
we compare the numerical solution of the born - markov master equation originating from the jaynes - cummings hamiltonian with the solution using the eigenstates of the full hamiltonian @xmath93 , which we found by applying van - vleck perturbation theory .
the fixed parameters we use are in units of @xmath16 : @xmath95 , @xmath393 , @xmath394 and @xmath395 . the oscillator frequency @xmath0 is varied .
for all the three possible cases ( positive , negative and zero detuning ) one notices from figures [ positivedetuning ]  [ nodetuning ] that the jaynes - cummings approach _ underestimates _ the dephasing rate @xmath396 ( means a larger peak at frequency @xmath140 ) and _ overestimates _ the rate @xmath397 ( smaller peak at @xmath141 ) compared to the approach with the full hamiltonian.the case of positive detuning ( @xmath398 ) is shown in figure [ positivedetuning ] . as here the dynamics is dominated by the frequency @xmath140 , the equilibrium value is reached on a too long time scale using the jaynes - cummings hamiltonian .    on the contrary , for negative detuning ( @xmath286 ) , which is shown in figure [ negativedetuning ] and where @xmath141 dominates
, the equilibrium value is reached too fast within the jaynes - cummings approach .
furthermore , considering the graph of the fourier transform we find that small contributions which come from higher level transitions , and which have already been discussed in , are not caught be the jaynes - cummings approach .    for the resonant case ( @xmath284 ) we find that the jaynes - cummings method predicts @xmath140 to be slightly dominating whereas the approach starting from @xmath6 results in @xmath141 being dominating .
the reason for this discrepancy is that due to the counter - rotating terms we have for @xmath6 no symmetric or antisymmetric superposition of the unperturbed eigenstates at @xmath284 in contrast to the jaynes - cummings model . to conclude this section
we can say that for an unbiased tls - ho system the jaynes - cummings model gives a good insight in the qualitative behaviour of @xmath49 both for a slightly detuned and a non - detuned system .
however , it under- or overestimates dephasing times for the system .
furthermore , we find taking into account counter - rotating terms in @xmath358 that at @xmath222 the dressed eigenstates are not a symmetric or antisymmetric superposition of the uncoupled states . moreover , the effects of transitions between states of different manifolds are neglected .
in conclusion , we discussed the dynamics of a biased and unbiased tls coupled through a harmonic oscillator to an environmental bath described by an ohmic spectral density . in particular , we examined the regime of weak damping and moderate coupling between oscillator and tls . an equivalent description of our system is provided by the spin - boson model with a structured spectral density .
in contrast to many other works in this field , our starting point was not the jaynes - cummings hamiltonian , but a more general one given in , where no initial rotating wave approximation has been applied .
in we provided with a formal expression for the population difference and showed in how the hamiltonian of the coupled qubit - ho system can be diagonalized approximately using van - vleck perturbation theory .
this approach is valid both for the oscillator being in resonance with the qubit ( @xmath222 ) and for finite detuning ( @xmath399 ) . in an analytical expression for the non - dissipative dynamics
was provided up to second order in the qubit - ho coupling @xmath1 taking into account the infinite hilbert space of the system . for low temperatures ( @xmath400 )
we truncated the hilbert space and found that transition processes between the groundstate and the two first excited energy levels of the qubit - ho system dominate the dynamics ( ) . in the influence of the bath was taken into account by solving the born - markov master equation for the density matrix of the qubit - ho system . to do this
analytically we considered two variants of the secular approximation : first , in , the full secular approximation ( fsa ) , where all fast oscillating terms are neglected , and second , in , the partial secular approximation , where attention was paid to the fact that the first two excited energy levels are almost degenerate . using an ansatz for the long time dynamics in , we could provide a general expression for the relaxation and dephasing rates of the qubit , showing that the relaxation time can be enhanced by detuning the oscillator into the off - resonant regime . it was found that all three approaches agree quite well with the numerical solution of the master equation .
the dynamics of both a biased and an unbiased qubit were intesively studied for zero and finite detuning in .
the results agree qualitatively with the numerical findings within the ab - initio quapi approach @xcite .
furthermore , in a good agreement with the results of the weak damping approximation performed in @xcite for a symmetric spin - boson model was found .
besides , we saw that at resonance ( @xmath222 ) the first two excited qubit oscillator states are not a symmetric or antisymmetric superposition of the states @xmath401 and @xmath325 , as predicted by the jaynes - cummings model , and thus could give an explanation for the differently weighted peaks of @xmath140 and @xmath141 in the fourier spectrum of the dynamics .
we further could explain the dominance of frequency @xmath141 in the case of negative detuning @xmath402 and of frequency @xmath140 for positive detuning , respectively .
moreover , we showed that for large negative detuning @xmath141 approaches the energy splitting @xmath403 of the qubit , whereas @xmath140 approximates the oscillator frequency @xmath0 .
this behaviour agrees nicely with spectroscopic experiments performed on a circuit qed architecture @xcite .
in we compared our results for an unbiased system to the ones obtained starting with the jaynes - cummings model .
we visualize the effects of the counter - rotating terms in the bloch - siegert shift of the resonance frequency and in contributions of states with larger oscillator number to the tls dynamics . apart from this
the jaynes - cummings model and our approach agree quite well for the non - dissipative case .
also for the dissipative case at resonance an initial rwa represents a good approximation to our starting hamiltonian and seems to be favourable as it is analytically exactly diagonalizable . for detuned systems , however , we find discrepancies concerning relaxation and dephasing times and a rwa becomes less appropriate to give precise results .
thus , we think that our approach represents an improvement as it is valid in a wider parameter range avoiding an initial rotating wave approximation . to our knowledge it provides for the first time analytical results for the dynamics of an unbiased _ and _ biased qubit coupled to a structured environment being valid both in the resonant _ and _ off - resonant regime .
furthermore , due to the generality of the qubit - oscillator model , we expect our results to be of interest for a wide range of experimental applications .
we acknowledge financial support under dfg program sfb631 .
further we would like to thank g. begemann for his help with some of the figures .
here , we show how to derive given in . in order to trace out the oscillator degrees of freedom we transform @xmath404 into the basis @xmath71 and by using @xmath405 we find @xmath406 performing the trace over the oscillator @xmath407 similarily
, we find for the off - diagonal elements of the reduced density matrix @xmath408 , \end{aligned}\ ] ] @xmath409.\end{aligned}\ ] ] using and in we can express @xmath49 in the energy basis , yielding @xmath410 + \sin \theta \left [ \rho_{\rm red;ge}(t)+\rho_{\rm red;eg}(t ) \right]\\     \fl = \cos(\theta ) \biggl (     \sum_{j=0}^\infty \sum_n \left [ { \langle j{\rm g}|n\rangle}^2 - { \langle j { \rm e}|n\rangle}^2 \right ] \rho_{nn } ( t ) \\      + \sum_{j=0}^\infty
\sum_{\stackrel{n , m}{n \neq m } } \left [ { \langle j { \rm g}|n\rangle } { \langle m|j { \rm g}\rangle } -   { \langle j { \rm g}|n\rangle } { \langle m|j { \rm g}\rangle } \right ] \re \{\rho_{nm}(t ) \ } \biggr ) \\      \fl \phantom{\mathrel{= } } + \sin ( \theta ) \biggl ( 2 \sum_{j=0}^\infty \sum_n { \langle j { \rm e}|n\rangle } \rho_{nn}(t ) { \langle n|j { \rm g}\rangle } \\
+ \sum_{j=0}^\infty \sum_{\stackrel{n , m}{n \neq m } } \left [ { \langle j { \rm e}|n\rangle}{\langle m|j { \rm g}\rangle } + { \langle j { \rm e}|m\rangle}{\langle j { \rm g}|n\rangle}\right ] \re \{\rho_{nm}(t)\}\biggr)}.\ ] ] with we can write @xmath411 where @xmath412   + 2 \sin \theta { \langle j { \rm g}|n\rangle } { \langle j { \rm e}|n\rangle } \right\ } \rho_{n n}(t ) , \\    \fl p_{nm}(t ) = 2 \sum_j \biggl\ { \cos \theta \biggl[{\langle j { \rm g}|n\rangle } { \langle m|j { \rm g}\rangle } - { \langle j { \rm e}|n\rangle } { \langle m|j { \rm e}\rangle }     \biggr ] \nonumber\\       + \sin \theta \biggl[{\langle j { \rm e}|n\rangle } { \langle m|j { \rm g}\rangle } + { \langle j { \rm e}|m\rangle } { \langle n|j { \rm g}\rangle }     \biggr ] \biggr\ } \re \{\rho_{nm}(t ) \}. \end{aligned}\ ] ]
let us consider the hamiltonian @xmath413 consisting of the free hamiltonian @xmath414 and a small perturbation @xmath415 , which is proportional to the coupling constant @xmath1 .
additionally we assume , that the energy levels @xmath416 of @xmath414 are grouped into manifolds , with @xmath46 being the index of the manifold and @xmath417 is used to distinguish between different energy levels within the same manifold .
the energy levels @xmath416 are eigenenergies of @xmath414 : @xmath418 through the perturbation @xmath419 different manifolds are coupled together . as long as the coupling @xmath1 is small , namely that @xmath420 for @xmath421 , also the energy levels of the total hamiltonian @xmath142 are clustered into manifolds . using the transformation @xmath422
, we construct an effective hamiltonian @xmath423 , which acts only within the individual manifolds ; i.e. , @xmath424 for @xmath421 , and has the same eigenvalues as @xmath142 within the manifolds .
we expand @xmath79 and @xmath423 in terms of the small parameter @xmath1 up to second order : @xmath425 for calculating @xmath426 and @xmath427 we use that @xmath428 for @xmath421 and furthermore choose that @xmath79 has no matrix elements within a manifold , namely @xmath429 .
now , one can iteratively calculate @xmath79 and @xmath423 order by order .
for details see @xcite . here
, we give only the results .
for the transformation one has @xmath430 and @xmath431 \nonumber \\
+ \sum_k \frac{1}{e_{j,\beta } - e_{j,\alpha } } \frac{{\langle j,\alpha| } v { |k,\beta\rangle } { \langle k,\beta| } v { |j,\beta\rangle}}{e_{k , \beta}-e_{j , \alpha } } \nonumber
\\          + \sum_k \frac{1}{e_{j,\beta } - e_{j,\alpha } } \frac{{\langle j,\alpha| } v { |k,\alpha\rangle } { \langle k,\alpha| } v { |j,\beta\rangle}}{e_{k , \alpha}-e_{j , \beta } } , \textrm{\quad for \quad } \alpha \neq \beta.\end{aligned}\ ] ] the effective hamiltonian is up to second order @xmath432 + \or(g^3).\nonumber\end{aligned}\ ] ]    in the case of the hamiltonian @xmath76 the first order matrix elements are @xmath433 and for the second order contributions @xmath434 where e. g. @xmath435 these matrix elements change sign under index transposition and all other matrix elements vanish .
finally , we get the transformation up to second order in @xmath1 : @xmath436
here , we give the matrix elements @xmath154 specified in . @xmath437
matrix elements including the groundstate are given separately because of the special shape of @xmath215 : @xmath438 all other matrix elements are zero .
here , we give the rate coefficients of the master equation for the reduced density matrix elements .
they are : @xmath439 @xmath440 @xmath441 @xmath442 @xmath443 @xmath444 @xmath445 @xmath446 @xmath447 @xmath448 @xmath449 @xmath450 @xmath451 @xmath452 @xmath453 @xmath454
the solutions of the fsa master equation for the diagonal elements reads : @xmath455 @xmath456 @xmath457 @xmath458 @xmath459
99 leggett a j 1987 1 weiss u 1999 _ quantum dissipative systems _
2nd edn ( singapore : world scientific ) grifoni m and hnggi p 1998 _ phys .
rep . _ * 304 * 229 garg a , onunchic j n and ambegaokar v 1985 _ j. chem .
phys . _ * 83 * 4491 raimond j m , brune m and haroche s 2001 565 yoshie t , scherer a , hendrickson j , khitrova g , gibbs h m , rupper g , ell c , shcheklin o b and deppe d g 2004 _ nature _ * 432 * 200 reithmaier j p 2004 _ nature _ * 432 * 197 nakamura y , pashkin yu a and tsai j s 1999 _ nature _ * 398 * 786 makhlin y , schn g and shnirman a 2001 357 vion d , aassime a , cottet a , joyez p , pothier h , urbina c , esteve d and devoret m h 2002 _ science _ * 296 * 886 collin e , ithier g , aassime a , joyez p , vion d and esteve d 2004 157005 blais a , huang r - s , wallraff a , girvin s m and schoelkopf r j 2004 a * 69 * 062320 wallraff a , schuster d i , blais a , frunzio l , huang r - s , majer j , kumar s , girvin s m and schoelkopf r j 2004 _ nature _ * 431 * 162 schuster d i , wallraff a , blais a , frunzio l , huang r - s , majer j , girvin s m and schoelkopf r j 2005 123602 wallraff a , schuster d i , blais a , frunzio l , majer j , devoret m h , girvin s m and schoelkopf r j 2005 060501 schuster d i 2007 _ nature _ * 445 * 515 wallraff a , schuster d i , blais a , gambetta j m , schreier j , frunzio l , devoret m h , girvin s m and schoelkopf r j 2007 050501 mooij j e , orlando t p , levitov l , tian l , van der wal c h and lloyd s 1999 _ science _ * 285 * 1036 van der wal c h , ter haar a c j , wilhelm f k , schouten r n , harmans c j p m , orlando t p , lloyd s and mooij j e 2000 _ science _ * 290 * 773 chiorescu i , nakamura y , harmans c j p m and mooij j e 2003 _ science _ * 299 * 1869 chiorescu i , bertet p , semba k , nakamura y , harmans c j p m and mooij j e 2004 _ nature _ * 431 * 159 johansson j , saito s , meno t , nakano h , ueda m , semba k and takayanagi h 2006 127006 blais a , gambetta j , wallraff a , schuster d i , girvin s m , devoret m h and schoelkopf r j 2007 a * 75 * 032329 tian l , lloyd s and orlando t p 2002 b * 65 * 144516 van der wal c h , wilhelm f k , harmans c j p m and mooij j e 2003 _ eur . phys . j. _
b * 31 * 2003 lupacu a , verwijs c j m , schouten r n , harmans c j p m and mooij j e 2004 177006 thorwart m , paladino e and grifoni m 2004 _ chem .
* 296 * 333 wilhelm f k , kleff s and von delft j 2004 _ chem . phys . _
* 296 * 345 goorden m c , thorwart m and grifoni m 2004 267005 goorden m c , thorwart m and grifoni m 2005 _ eur .
phys . j. _
b * 45 * 405 nesi f , grifoni m and paladino e 2007 316 kleff s , kehrein s and von delft j 2003 _ physica e _ * 18 * 343 kleff s , kehrein s and von delft j 2004 b * 70 * 014516 huang p and zheng h 2008 _ preprint _ cond - mat/0707.0725v3 makri n and makarov d e 1995 _ j. chem . phys . _ * 102 * 4600 makri n and makarov d e 1995 _ j. chem .
phys . _ * 102 * 4611 jaynes e t and cummings f w 1963 _ ieee proc . _ * 51 * 90 shore b w and knight p l 1993 _ j. mod .
_ * 40 * 1195 brito f and caldeira a o 2008 _ preprint _ quant - ph/08060992v1 shavitt i and redmon l t 1980 _ j. chem .
* 73 * 5711 cohen - tannoudji d , dupont - roc j and grynberg g 1992 _ atom - photon interactions : basic processes and applications _
( new york / chichester / brisbane / toronto / singapore : john wiley & sons ) caldeira a o and leggett a j 1983 374 blum k 1996 _ density matrix theory and applications _
2nd edn ( new york and london : plenum press ) louisell w h 1973 _ quantum statistical properties of radiation _
( new york / london / sydney/ toronto : john wiley & sons ) bloch f and siegert a 1940 522 vogel w , welsch d g and wallentowitz s 2001 _ quantum optics : an introduction _
2nd edn ( berlin : wiley - vch )