since the foundation of quantum physics there has been huge interest in the nontrivial transition from the quantum to the classical regime . an important observation in this context is the ehrenfest theorem stating that the dynamics of quantum mechanical expectation values is determined by the classical equations of motion @xcite .
going beyond expectation values , an _
ehrenfest time _
scale has been identified as the time below which the quantum time evolution is well approximated by the corresponding classical dynamics @xcite .
the ehrenfest time is the time it takes for the chaotic classical dynamics to stretch an initially narrow wave packet to some relevant classical length scale such as the system size @xmath5 .
since the stretching is exponential in classically chaotic systems , one has @xmath6 with the lyapunov exponent @xmath7 of the classical dynamics , the initial spread @xmath8 of the wave packet and the fermi momentum @xmath9 .    in recent years
, there has been much interest in determining the influence of the ehrenfest time on stationary transport quantities such as the conductance @xcite , its variance @xcite and its behavior under decoherence @xcite , shot noise @xcite and higher moments of the current @xcite , and on time dependent quantities such as the spectral form factor @xcite , the survival probability @xcite and the fidelity @xcite .
most of these papers used the trajectory - based semiclassical approach to transport @xcite which currently is the method of choice for investigating ehrenfest - time dependences of quantum observables @xcite .
the leading order quantum correction to the conductance was found in refs .
@xcite to decay exponentially with the ehrenfest time . qualitatively speaking this can be understood by noting first that this contribution originates from loop diagrams [ see fig .
[ fig : orbt](b ) ] , and as such depends on the return probability , and second that , neglecting system - dependent nongeneric processes , this return probability contribution essentially vanishes for times shorter than the ehrenfest time .
such an intuitive interpretation of the ehrenfest - time dependence does not always work : though inherently of nonclassical nature , the leading order contribution to the conductance variance of systems ideally coupled to external leads turns out to be independent of the ehrenfest time @xcite .
( gray boxes ) .
a classical trajectory traversing the system is shown by a solid ( blue ) line .
the second dashed ( red ) line on the left indicates a trajectory backreflected at the barrier .
this occurs with probability @xmath10 . ]    in this paper , we determine semiclassically the ehrenfest - time dependence of the variance @xmath11 of the conductance for a chaotic system coupled to external leads via nonideal contacts modeled by tunnel barriers of transparency @xmath12 .
the situation is depicted in fig .
the presence of tunnel barriers has the dramatic effect that @xmath11 increases or decreases with @xmath0 , depending on the value of @xmath3 . for @xmath4
, we even observe an enhancement of the variance above the universal value in the presence of time reversal symmetry for equal lead widths of @xmath13 @xcite upon increasing @xmath0 .
this is very surprising , given the quantal nature of the conductance fluctuations . in the range @xmath14 we find a reduction of @xmath11 as @xmath0 increases that is strongest around @xmath15 , and recover the @xmath0-independent behavior of @xmath16 at @xmath17 observed in refs .
@xcite . the precise dependence on @xmath3 is depicted in fig .
[ condfluc ] .    .
the function @xmath18 , defined in eqs .
( [ li22.9 ] ) and ( [ li23 ] ) , measures the deviation from the universal rmt variance . in this graph
the parameter @xmath19 is plotted as a function of the tunnel barrier transparency @xmath3 for time - reversal symmetric systems .
the solid red line gives the analytical prediction , eq .
( [ li23 ] ) , and the blue dots are results of numerical simulations .
error bars indicate the standard deviation over the ensemble of calculated data . , scaledwidth=47.0% ]    to obtain these results we first , in sec
.  [ semiapprox ] , introduce the semiclassical approximation to the conductance variance .
then , in sec
.  [ part1 ] , we analytically calculate the @xmath3-dependence of the diagrammatic contributions to the variance , to leading order in the total number of open channels and linear in the ehrenfest time .
we list all relevant diagrams and calculate their contributions . given their number , we also identify the most relevant ones and specify the range in @xmath3 where they are particularly important . in sec .
[ part2 ] we present numerical results that confirm our analytical results and finally conclude in sec .  [ concs ] .
within the scattering approach to quantum transport @xcite , the energy averaged conductance @xmath20 ( in units of @xmath21 ) can be expressed in terms of the transmission @xmath22 as @xmath23 with @xmath24 denoting an average over an energy window that is classically small but quantum mechanically large .
this leads to the following expression for the variance @xmath25 ^ 2\right\rangle-\left\langle{\rm tr}\left(\textbf{tt}^\dagger\right)\right\rangle^{2}.\ ] ] the scattering matrix elements are related by the fisher - lee relation @xcite to the projection of the green function onto the transverse directions in the leads .
performing the projection to leading order in @xmath26 and approximating the green function semiclassically , one obtains @xmath27 with the heisenberg time @xmath28 , the time conjugate to the mean level spacing . here
the sum is over the scattering trajectories @xmath29 which connect channel @xmath30 in the entrance ( or say left ) lead and channel @xmath31 in the exit ( or right ) lead in fig .
the summands contain rapidly oscillating phases depending on the classical actions @xmath32 of the considered classical trajectories , and classical stability prefactors @xmath33 whose precise form is given for example in @xcite .
inserting eq .
( [ in6 ] ) into ( [ in4.9 ] ) we obtain the semiclassical expression for the conductance , @xmath34 using ( [ in6 ] ) in ( [ in5 ] ) yields the semiclassical expression its variance , @xmath35 with the channel sums in eqs .
( [ in6.9 ] ) and ( [ in7 ] ) running over all open lead channels ( @xmath36 in the left and @xmath37 in the right lead ) . if we consider contributions in the first term in eq .
( [ in7 ] ) where @xmath29 and @xmath38 form a correlated pair ( with self - encounters ) and @xmath39 and @xmath40 form a separate correlated pair , we simply recreate the second term .
we can thus remove the second term in the above equation by removing such pairs from the semiclassical treatment of the first term . in terms of trajectories
we then obtain @xmath41 where the trajectories @xmath42 go from channel @xmath30 in the entrance lead to channel @xmath31 in the exit lead .
likewise trajectories @xmath43 go from channel @xmath44 to channel @xmath45 . because we have removed terms from correlated trajectories where @xmath46 and @xmath47 [ this restriction is denoted by the double bracket in eq .
( [ condvartraj ] ) ] we are left with quadruplets where all four trajectories interact through encounters .    before performing the energy average the approximations for @xmath20 and @xmath16 in eqs .
( [ in6.9 ] ) and ( [ condvartraj ] ) are rapidly fluctuating as a function of energy for @xmath48 .
thus only contributions from very similar trajectories survive the average .
the classical contribution to eq .
( [ in6.9 ] ) results from equal trajectories @xmath49 , the so called diagonal approximation @xcite , for an illustration see fig .
[ fig : orbt](a ) .
( b).,scaledwidth=45.0% ]    here eq .
( [ in6.9 ] ) yields @xmath50}=\frac{1}{\th}\sum_{a , b } \sum_{\gamma(a\to b)}\left|a_\gamma\right|^2.\ ] ] from here on , superscripts refer to the corresponding figure . to evaluate the remaining @xmath29-summation the sum rule for open systems @xcite that transforms the sum over orbits into an integral over their durations
is applied @xmath51 here @xmath1 is the dwell time of the cavity , i.e.  the typical time a classical particle remains inside the chaotic system .
this is given by @xmath52 with @xmath53 .
this finally yields @xmath54}\approx\frac{n_{\rm l}n_{\rm r}}{n_{\rm l}+n_{\rm r}}.\ ] ] quantum corrections to this result are obtained from pairs of slightly different trajectories .
the considered trajectories are almost identical differing only in how they are connected within self encounters , for the pair considered in @xcite , see fig .
[ fig : orbt](b ) . here
the orbits possess close self encounters with one orbit crossing and the other anticrossing there leading to a different direction of traversal of the closed loop . considering this pair in eq .
( [ in6.9 ] ) leads to the leading order quantum correction to the conductance . to determine its contribution the action difference between the partner trajectories and the number of crossings needs to be determined .
the calculation is done here within the phase - space approach , in the context of the conductance it was first performed in ref .
we will follow the latter approach throughout this paper .
there a poincar surface of section is considered inside the encounter region and the difference along the stable and unstable directions of the piercing points of the two stretches , @xmath55 and @xmath56 , respectively , is used to characterize an encounter . in terms of these coordinates the action difference for the orbit pair in fig .
[ fig : orbt ] is given by @xmath57 @xcite .
the weight @xmath58 that additionally depends on the duration of the orbit @xmath59 measuring the number of encounters is obtained from the ergodicity of the flow as @xcite @xmath60 where @xmath61 is the duration of the encounter . in general , the action difference @xmath62 and the weight @xmath63 depend on the trajectory configuration considered . in total
we obtain for the quantum correction @xmath64 resulting from the diagram in fig .
[ fig : orbt](b ) @xmath65}&\approx&\int_{-c}^c\!\!{\rm d}s{\rm d}u \!\!\int_{2t_{\rm enc}}^\infty\!\!{\rm d}t \,w({s},{u}){\rm e}^{\frac{i\delta s}{\hbar } } { \rm e}^{-\frac{\left(t - t_{\rm enc}\right)}{\tau_d}}\nonumber\\&\approx&-\frac{n_{\rm l}n_{\rm r}}{\left(n_{\rm l}+n_{\rm r}\right)^2}{\rm e}^{-\te/\td}\end{aligned}\ ] ] with @xmath66 . in the first line additionally
the survival probability correction during the encounter @xcite is taken into account . during the encounter
the stretches are so close that the orbit escapes either during the first stretch or does not escape at all leading to the enhanced survival probability in eq .
( [ quantumco ] ) .
the @xmath67-integrals in eq .
( [ quantumco ] ) are performed as described in @xcite , by substituting @xmath68 and @xmath69 .
the @xmath70-integral then essentially cancels the @xmath71 in the denominator and the @xmath72-integral yields ( after a partial integration ) the contribution @xmath73 .    to treat orbits differing at several places , ref .
@xcite introduces the splitting of the orbit into encounters and links . during the encounters the orbits are close to themselves but the orbit and its partner are differently connected . due to the exponential separation of neighboring trajectories in the chaotic case
these last essentially an ehrenfest time , as will become clear from the calculation below .
the stretches are connected by the links , where ` links ' denote the long parts of the trajectory where the trajectory and its partner are essentially identical ( up to time reversal symmetry ) . with a suitable change of variables in the calculation , in the rmt limit @xmath74
, one can treat different encounters as distinct and separate the semiclassical contribution into a product of contributions over each of the links and encounters .
the total contribution can therefore be obtained by diagrammatic rules @xcite .    away from this limit , and for the ehrenfest - time dependence , the encounters may start to overlap and for the conductance variance the trajectories can be seen to meet and surround periodic orbits trapped inside the system , i.e.  they have encounters with periodic orbits , and must be treated as part of a continuous family @xcite .    when we include tunnel barriers ( as in fig .
[ sys ] ) , three main changes occur that were originally described in @xcite :    1 .
the particles enter and leave the cavity with the probability @xmath3 , which leads to a factor @xmath75 for each trajectory pair .
2 .   while without tunnel barriers every particle that hits the lead leaves the system , now only the ratio @xmath3 of the particles hitting the lead leaves the cavity . for the links ,
the effective dwell time is therefore @xmath76 and the dwell time in the exponential in ( [ in8 ] ) should be replaced by this effective dwell time .
however , if trajectory stretches are correlated , as they are during encounters , then the whole configuration is lost if just one stretch of the encounter leaves the system .
this happens with a probability @xmath77 for @xmath78 correlated stretches .
the dwell time in such a situation is therefore replaced by @xmath79 $ ] .
additional encounter diagrams become possible ; for an example which contributes to the energy averaged conductance @xmath20 , see fig .
[ fig : orbt](c ) . in this case
one encounter stretch can be moved into the lead forcing the other to be backreflected at the opening .
note that configurations where both stretches are backreflected at the opening are already taken into account by the modified dwell time explained above .
although the effective dwell times are altered by the tunnel barriers , the action difference and weight functions are unaffected so that in the rmt limit ( @xmath74 ) contributions can still be obtained by diagrammatic rules : the contribution of each link is now given by @xmath80 .
the stretches of an encounter of @xmath78 orbits yield @xmath81 . for the ehrenfest - time dependence ,
however , these changes render the calculation of the contributions to the conductance variance much more difficult compared to @xmath82 : due to the discontinuous form of the effective dwell time , the contributions from diagrams with a different number of surroundings of trapped periodic orbits need to be split and treated separately .
also the possibility for encounter stretches to be backreflected at the tunnel barriers increases the number of diagrams considerably .
here we calculate the leading order in @xmath83 contributions to the variance of the conductance for non - zero ehrenfest time in the presence of tunnel barriers .
we show all the relevant diagrams and calculate their contributions .
the results given here are valid in the unitary case ; results for the orthogonal case can be obtained by multiplying the total by a factor of @xmath84 .
we start with the contributions important in the rmt limit which allow us to recover the rmt result .
the corresponding rmt calculation was performed in @xcite by brouwer and beenakker .
first we consider two 2-encounters ( encounters involving two trajectory stretches ) in a row , see fig .
[ ana1](a ) .
this diagram also occurred in the @xmath82 treatment of @xcite , but we will explain how for @xmath85 other diagrams with backreflected stretches can be derived from this one .        when the encounters are inside the cavity and not touching the tunnel barriers , as depicted in fig .
[ ana1](a ) , we start with eq .
( [ condvartraj ] ) and use the sum rule ( [ in8 ] ) , the action difference @xmath86 ( where the subscripts refer to the different encounters ) and the weight function @xcite @xmath87 containing the durations of the two trajectories indicated by solid lines in fig .
[ ana1](a ) , @xmath88 and @xmath89 , respectively and the durations of the two encounters of the trajectories @xcite @xmath90 with a classical constant @xmath44 again of the order one . from eqs .  ( [ in7 ] ) and ( [ in8 ] ) we write @xmath91}&= & \frac{\gamma^4n_{\rm l}^2n_{\rm r}^2}{\th^2}\prod_{i=1}^6 \int_0^\infty { \rm d}t_i\,{\rm e}^{-\gamma t_i/\td } \int_{-c}^c{\rm d}{\boldsymbol s}{\rm d}{\boldsymbol u}\frac{1}{\omega^2}\nonumber\\ & & \times\frac{{\rm e}^{\frac{{\rm i}}{\hbar}{\boldsymbol s}{\boldsymbol u } } } { t_{{\rm enc},1}t_{{\rm enc},2}}{\rm e}^{-p_2\left(t_{{\rm enc},1}+t_{{\rm enc},2}\right)/\td } , \end{aligned}\ ] ] with the superscript again referring to the corresponding figure : fig .
[ ana1 ] . as explained in sec .
[ semiapprox ] , the trajectory quadruplet leads to the overall factor @xmath92 while the sum over possible channels provides @xmath93 .
the six links have an effective dwell time of @xmath76 , while each 2-encounter experiences the dwell time @xmath94 as explained at the end of the last section .
the @xmath95-integrals are performed , as described after eq .
( [ quantumco ] ) , by substituting @xmath96 and @xmath97 .
each @xmath98-integral yields ( after a partial integration ) the contribution @xmath99 as already obtained for the conductance for @xmath82 .
finally the @xmath100-integrals yield @xmath101 } = \frac{n_{\rm l}^2n_{\rm r}^2}{\left(n_{\rm l}+n_{\rm r}\right)^4 } \frac{p_2 ^ 2}{\gamma^2}{\rm e}^{-2p_2\te/\td},\ ] ] which generalizes the result for @xmath82 from @xcite .
the representation of eq .
( [ li0.2 ] ) , and the integrals to arrive at ( [ li0.3 ] ) , also nicely illustrate how the diagrammatic rules introduced above arise in this context .
we now turn to the new diagrams that arise due to the tunnel barriers where some of the links are shrunk until an encounter touches a barrier .
first we consider the case in which just one link connecting the encounter to the opening is removed and the corresponding encounter stretch now starts in the opening ; as in the example in fig .
[ fig : orbt](b ) . as the stretches during an encounter lie very close to each other ,
the other encounter stretch has to be backreflected at the opening ( so only one link is lost ) .
this contribution can therefore only exist for @xmath85 .
the changes in the analytical calculation mainly affect the weight function , see also @xcite . if we shrink a link on trajectory 1 then , compared to eq .
( [ li0 ] ) , the orbit of duration @xmath88 now only involves two links so the factor @xmath102 is replaced by @xmath103 . moreover , for the encounter which touches the barrier , @xmath104 is replaced by an integration variable @xmath105 that is integrated from zero to the @xmath104 defined in ( [ li0.1 ] ) .
this variable measures the length of the encounter that remains inside the system , i.e.  which has not yet been moved into the lead .
performing again the steps described after ( [ li0.1 ] ) yields an expression similar to ( [ li0.2 ] ) but with @xmath106 instead of @xmath107 link factors and @xmath104 in the exponential in the second line replaced by an integration variable @xmath105 that is integrated from zero to @xmath104 . because half of the encounter is backreflected at the tunnel barrier , additionally this contribution
is multiplied by @xmath108 .
the contribution from the lower limit of the @xmath105-integral is zero @xcite , leading to @xmath109 } = -\frac{4n_{\rm l}^2n_{\rm r}^2}{\left(n_{\rm l}+n_{\rm r}\right)^4 } \frac{(1-\gamma)p_2}{\gamma}{\rm e}^{-2p_2\te/\td},\ ] ] with the ` @xmath110 ' denoting that one link was removed .
the prefactor 4 is due to the fact that there are four such links we can remove .
analogously , the contribution where two links that connect two _ different _ encounters to the opening are removed is @xmath111}= \frac{4n_{\rm l}^2n_{\rm
r}^2}{\left(n_{\rm l}+n_{\rm r}\right)^4 } ( 1-\gamma)^2{\rm e}^{-2p_2\te/\td},\ ] ] where there are again four possibilities of picking two links to remove .    additionally we can , and this is a possibility which also exists for @xmath82 , remove both links connecting the _ same _ encounter to the opening .
this means that the trajectories tunnel straight into the encounter . as the encountering orbits are so close together this means that @xmath112 or @xmath113 so that there is only one channel summation in the lead where the encounter touches .
also we have two links fewer and one integral over the part of the encounter that remains inside the system , yielding @xmath114}= -\frac{\left(n_{\rm l}^2n_{\rm r}+n_{\rm l}n_{\rm r}^2\right)}{\left(n_{\rm l}+n_{\rm r}\right)^3 } p_2{\rm e}^{-2p_2\te/\td},\ ] ] where the additional ` ( s ) ' in the superscript indicates that two links were removed at the same encounter .
we can further remove one link from the other encounter to obtain @xmath115}= \frac{2\left(n_{\rm l}^2n_{\rm r}+n_{\rm l}n_{\rm r}^2\right)}{\left(n_{\rm l}+n_{\rm r}\right)^3 } \gamma(1-\gamma){\rm e}^{-2p_2\te/\td}.\ ] ] finally , when all four links connecting the encounter to the leads are removed , we have @xmath116}= \frac{n_{\rm l}n_{\rm r}}{\left(n_{\rm l}+n_{\rm r}\right)^2}\gamma^2 { \rm e}^{-2p_2\te/\td}.\ ] ] with time reversal symmetry , however , we can also have fig .
[ ana1](b ) where , because channels @xmath30 and @xmath44 are in the left lead and channels @xmath31 and @xmath45 in the right lead , we can not shrink more than two links simultaneously .
similarly , we can not remove two links from the same encounter .
since @xmath117 , the contributions ( [ li0.6][li0.7 ] ) actually cancel so that the diagram in fig .
[ ana1](b ) and the ones obtained by cutting links provide the same contribution as the diagram in fig .
[ ana1](a ) .
time reversal symmetry then still gives a factor 2 in this case , while for all the following cases , diagrams related by time reversal symmetry provide the same contributions directly .
next we treat the diagrams with two 2-encounters which lie along a trapped periodic orbit ; depicted in fig .
although the links connect the encounter stretches in a different way compared to fig .
[ ana1 ] these diagrams again contain two 2-encounters and the same number of links .
the possibilities for shrinking links are identical as for fig .
[ ana1](b ) and so each diagram provides the same contributions as above and the same total contribution as the configuration in fig .
[ ana1](a ) .
they are multiplied , however , by a factor @xmath84 since the orbits here have two possibilities to traverse the enclosed periodic orbits ( schematically we can also reflect the diagrams horizontally ) .
the last relevant diagram type to obtain the rmt result is shown in fig .
this diagram contains one 3-encounter with a periodic orbit .
its contribution is calculated in an analogous way to eq .
( [ li0.2 ] ) .
here we have one link fewer , one encounter instead of two and involving three rather than two orbit stretches .
therefore we have @xmath118}= -\frac{2n_{\rm l}^2n_{\rm r}^2}{\left(n_{\rm l}+n_{\rm r}\right)^4 } \frac{p_3}{\gamma}{\rm e}^{-p_3\te/\td}.\ ] ] alternatively this result can be obtained from the contribution of a 3-encounter , @xmath119 , in the appendix of @xcite by including the correct dwell times as well as the contributions from the links .    also in this case encounter stretches can be shrunk and removed : first removing one link connecting the encounter to the opening yields @xmath120}=
\frac{8n_{\rm l}^2n_{\rm r}^2}{\left(n_{\rm l}+n_{\rm r}\right)^4 } ( 1-\gamma)^2{\rm e}^{-p_3\te/\td}\ ] ] and second removing both links : @xmath121}= \frac{2\left(n_{\rm l}n_{\rm r}^2+n_{\rm l}^2n_{\rm r}\right)}{\left(n_{\rm l}+n_{\rm r}\right)^3 } \gamma(1-\gamma){\rm e}^{-p_3\te/\td}.\ ] ] having calculated all the contributions in the rmt limit , we can obtain the rmt result by setting @xmath122 . when summing the the results in eqs .
( [ li0.3][li0.91 ] ) , we obtain @xmath123 , \end{aligned}\ ] ] which agrees with the rmt prediction in @xcite . in appendix
[ diffprobs ] , we use these diagrams to obtain the rmt result when each channel has a different tunneling probability .
having gone through all the diagrams that contribute at zero ehrenfest time , we now turn to those diagrams whose contribution vanishes at zero ehrenfest time . for these contributions , the periodic orbits in figs .
[ ana2 ] and [ ana3 ] become important and we now view those diagrams as trajectories which have an encounter with the periodic orbit , rather than with each other . for example , in fig .
[ ana3 ] we could start with the solid trajectory which passes from channel @xmath44 to @xmath45 and the dashed trajectory from @xmath30 to @xmath31 and build the rest of the diagram from those starting points and the periodic orbit .
both of those trajectories encounter the periodic orbit once . in the semiclassical treatment of fig .
[ ana3 ] above , it was implicitly assumed that these encounters occur at the same point along the periodic orbit .
the resulting 3-encounter can therefore be considered as an ` aligned ' 3-encounter , but for the further ehrenfest time dependence we also need to consider the situation where the two encounters with the periodic orbit occur at different points along the periodic orbit but still overlap . in this case
we have a ` non - aligned ' 3-encounter , while when the encounters no longer overlap we return to the two separate 2-encounters of fig .
[ ana2 ] .
we now derive the ehrenfest - time dependence of a non - aligned 3-encounter , whose base trajectories are depicted in fig .
[ ana4](a ) . in order to obtain the complete trajectory quadruplet that contributes in ( [ condvartraj ] ) we first need to include an additional traversal of the periodic orbit with one of the base trajectories to obtain the original trajectory pair [ which has positive action in ( [ condvartraj ] ) ] . the partner pair ( with negative action )
is then created by including the extra traversal of the periodic orbit with the other base trajectory . in this way
we recover a diagram like fig .
[ ana3](a ) from fig .  [ ana4](a ) and a small action difference ( the action of the periodic orbit itself cancels ) . after writing the contribution of a non - aligned 3-encounter in an analogous way as for the aligned 3-encounter , as explained in eqs .
( [ li0][li0.2 ] ) , as product of link and encounter contributions , it can be evaluated by making use of the results for @xmath124 in the appendix of @xcite . of course , with the tunnel barriers , the dwell times must be modified compared to @xcite : the link dwell time is modified by the factor @xmath125 , while the parts of the encounter where only one trajectory is correlated to the periodic orbit  these contributions are called fringes in @xcite  have the factor @xmath126 since we have two stretches close to each other .
likewise , when both trajectories are correlated with the periodic orbit we have three stretches in total and the corresponding factor @xmath127 . with these corrections ,
the contribution is @xmath128}&=&\frac{4n_{\rm l}^2n_{\rm r}^2}{\left(n_{\rm l}+n_{\rm r}\right)^4}\frac{p_2 ^ 2 } { \gamma\left(2p_2-p_3\right)}\nonumber\\ & & \times\left({\rm e}^{-p_3\te/\td}-{\rm e}^{-2p_2\te/\td}\right ) . \end{aligned}\ ] ] in order to proceed to the additional diagrams that arise from touching the tunnel barriers , we first reconsider this contribution in detail along the lines of @xcite . as for eqs .
( [ li0][li0.2 ] ) , we start with the weight function for two base trajectories encountering a periodic orbit , as explained in @xcite @xmath129 here @xmath130 and @xmath131 are the durations of the links of the base trajectories that connect the periodic orbit encounters to the lead , while @xmath132 is the period of the periodic orbit .
the integral over the period corresponds to summing over all periodic orbits which can be encountered .
the encounter times are as in ( [ li0.1 ] ) , but using the stable and unstable distances between the encounter stretches and the periodic orbit itself .
finally , @xmath105 measures the time difference between the midpoints of the two encounter stretches in fig .
[ ana4 ] , and the @xmath105-integral covers the different arrangements of the stretches relative to each other .
the limits of the two last integrals in ( [ li1.1 ] ) are not given as they depend on the specific configuration considered below .
the first two integrals in ( [ li1.1 ] ) can again be transformed , together with the integrals from the sum rule ( [ in8 ] ) , into a product of link and encounter contributions . the general expression for the contributions to the variance from diagrams containing an enclosed periodic orbit that does not touch
the lead becomes @xmath133 like the limits of the @xmath132- and @xmath105-integrals , the function @xmath134 determining the survival probability along the periodic orbit ( and the encountering trajectory stretches ) will be specified for each contribution separately below .
for example , for the non - aligned 3-encounter , we can arrange the different alignments in terms of the durations of the two encounters .
we let @xmath135 denote the longer encounter and @xmath136 the shorter . in fig .
[ ana4](a ) , we do not yet allow the encounter stretches to surround the periodic orbit so we impose the restriction @xmath137 .
then we can first consider the case that the shorter encounter lies inside the longer .
we will also refer to this case later as a ` generalized ' 3-encounter .
the time difference between the midpoints of the encounters therefore satisfies @xmath138 , while the survival probability function in ( [ li1.2 ] ) is given by @xmath139 the different terms simply correspond to the parts of the periodic orbit which are followed by one , two and three trajectory stretches respectively .
as this survival probability is independent of @xmath105 , the @xmath105-integral simply yields @xmath140 .
performing the remaining integrals , we obtain @xmath141 of course , the range of @xmath105 includes 0 , the case in which the encounters are perfectly aligned .
( [ li1.755 ] ) therefore includes the contribution of the 3-encounter in fig .
[ ana3](a ) , and for just the non - aligned 3-encounter we subtract ( [ li0.8 ] ) to give @xmath142 now we consider the case in which the shorter encounter is no longer fully inside the longer , but where the two encounters still overlap .
this can be further separated according to the total length of the encounters , @xmath143 , compared to the periodic orbit .
when @xmath144 , we have the range @xmath145 .\ ] ] and the survival probability @xmath146 where the terms are now expressed as corrections due to additional correlated stretches . integrating the exponential depending on this survival probability in ( [ li1.2 ] ) with respect to @xmath105 yields @xmath147/\td}\right .
\\ & & \left.-{\rm e}^{-\left[p_3t_{{\rm enc , min}}+p_2\left(t_{{\rm enc , max}}-t_{{\rm enc , min}}\right)+\gamma\left(\tau_p - t_{{\rm enc , max}}\right)\right]/\tau_d}\right ) .\nonumber\end{aligned}\ ] ] performing the remaining integrals in ( [ li1.2 ] ) leads to @xmath148 \nonumber\end{aligned}\ ] ] if the encounters are longer than the periodic orbit , @xmath149 , but we still do not allow them to overlap at both ends , we have the restriction @xmath150 , \ ] ] while the survival probability remains as in ( [ li1.3 ] ) . the @xmath105-integral yields @xmath151/\td } \right . \\ & & \left.-{\rm e}^{-\left[p_3t_{{\rm enc , min}}+p_2 \left(t_{{\rm enc , max}}-t_{{\rm enc , min}}\right ) + \gamma\left(\tau_p - t_{{\rm enc , max}}\right)\right]/\td}\right ) , \nonumber \end{aligned}\ ] ] which finally leads to @xmath152\nonumber\\ & & -\frac{1}{\gamma\left(p_3 - 2p_2+\gamma\right ) } \left[\left(p_3-p_2\right)p_2{\rm e}^{-p_3\te/\td}\right.\nonumber\\ & & \left.\left .- p_2\left(p_3-p_2+\gamma\right){\rm e}^{-\left(p_3+\gamma\right ) \te/\td}\right ] \right\rbrace .\end{aligned}\ ] ] as can be easily checked , the sum of the contributions of the non - aligned 3-encounter in eqs .
( [ li1.756],[li1.6],[li1.7 ] ) equals the contribution calculated directly in ( [ li1 ] ) .
the reason why we have discussed this more complicated route , is that we can use it to easily calculate the contributions when we start to shrink links and allow the encounter to touch the tunnel barriers .
for example , if we remove one link , then by performing analogous steps as explained before ( [ li0.4 ] ) , we obtain @xmath153}&= & \frac{8n_{\rm l}^2n_{\rm r}^2}{\left(n_{\rm l}^2+n_{\rm r}^2\right)^4 } \frac{p_2\left(1-\gamma\right ) } { \left(2p_2-p_3\right)}\left[1+\left(1-\gamma\right)\right ] \nonumber \\ & & \times\left({\rm e}^{-2p_2\te/\td}-{\rm e}^{-p_3\te/\td}\right ) ,   \end{aligned}\ ] ] where in the square brackets , the @xmath154 results from configurations where the enclosed periodic orbit touches the tunnel barrier where only one encounter stretch is correlated with the orbit , and the @xmath155 from configurations where the periodic orbit touches the tunnel barrier while both encounter stretches are correlated with it . likewise ,
if two links connecting _ different _ stretches to the lead are removed , we obtain the contribution @xmath156}&= & \frac{4n_{\rm l}^2n_{\rm r}^2}{\left(n_{\rm l}^2+n_{\rm r}^2\right)^4 } \frac{\gamma\left(1-\gamma \right)^2\left(2-\gamma\right)^2 } { \left(2p_2-p_3\right ) } \nonumber \\ & & \times\left({\rm e}^{-p_3\te/\td}-{\rm e}^{-2p_2\te/\td}\right ) .
\end{aligned}\ ] ]    next we turn to the configuration in fig .
[ ana4](b ) with encounter stretches overlapping at both ends .
we still have @xmath149 , but a different restriction on @xmath105 : @xmath157 .\ ] ] the survival probability is again independent of @xmath105 , and given by @xmath158 so that the @xmath105-integral yields @xmath159 . performing the remaining integrals we obtain @xmath160}&= & \frac{2n_{\rm l}^2n_{\rm
r}^2}{\left(n_{\rm l}+n_{\rm r}\right)^4 } \left[\frac{p_2 ^ 2}{\left(2p_2-p_3\right)^2}\right .
\nonumber\\ & & \times \left({\rm e}^{-2p_2\te/\td}-{\rm e}^{-p_3\te/\td}\right ) \nonumber\\ & & \left .
+ \frac{p_2\left(p_3-p_2\right)\te}{\left(2p_2-p_3\right)\td } { \rm e}^{-p_3\te/\td } \right].\end{aligned}\ ] ] again the encounter stretches can be brought to the lead by first removing one link @xmath161}&= & \frac{8n_{\rm l}^2n_{\rm r}^2}{\left(n_{\rm l}+n_{\rm r}\right)^4 } \left[\frac{p_2\gamma(1-\gamma)^2}{\left(2p_2-p_3\right)^2}\right.\nonumber\\ & & \times \left({\rm e}^{-p_3\te/\td}-{\rm e}^{-2p_2\te/\td}\right ) \nonumber \\ & & \left.-\frac{p_3\gamma(1-\gamma)^2\te}{2\left(2p_2-p_3\right)\td } { \rm e}^{-p_3\te/\td } \right ] , \end{aligned}\ ] ] and second by removing two links connecting _
different _ stretches to the opening @xmath162}&= & \frac{8n_{\rm
l}^2n_{\rm r}^2}{\left(n_{\rm l}+n_{\rm r}\right)^4 } \left[\frac { \gamma^2(1-\gamma)^4}{\left(2p_2-p_3\right)^2}\right.\nonumber\\ & & \times\left({\rm e}^{-2p_2\te/\td}-{\rm e}^{-p_3\te/\td}\right ) \nonumber\\ & & \left.+\frac{\gamma^2(1-\gamma)^4\te}{\left(2p_2-p_3\right)\td } { \rm e}^{-p_3\te/\td}\right ] .
\end{aligned}\ ] ]      as the durations of the encounters are of the order of the ehrenfest time , for vanishing ehrenfest time we only considered the situation where the encounter stretches could be partially reflected from the tunnel barriers in _
one _ lead .
however , for increasing ehrenfest time configurations where encounter stretches are partially reflected from the tunnel barriers in _ both _ leads , i.e.  where they touch the opening at both ends , become important . if one of the 2-encounters in fig .
[ ana1 ] were to be partially reflected at both ends then one of the links between the two encounters would need to tunnel through the barrier and exit the system so that the rest of the diagram could not be completed . in figs .
[ ana2 ] and [ ana3 ] , however , as long as the trajectory stretches which follow the periodic orbit are reflected at the tunnel barriers and remain in the system , we can allow the other links to tunnel through the barrier and exit the system .
for the base trajectories in fig .
[ ana4](a ) this means we can allow both links of one of the base trajectories , on either side of the _ same _ encounter , to shrink into the start and end lead .
we start with the configuration in fig .
[ ana2 ] and for the calculation we note that the contributions of the different orbit parts in ( [ condvartraj ] ) are multiplicative @xcite .
we can therefore reconnect the orbits in such a way that they split into parts whose contributions have previously been calculated .
first there is the remaining base trajectory , with a 2-encounter with the enclosed periodic orbit and two links connecting it to the opening .
this contributes the factor @xmath163 the rest of the diagram involves the periodic orbit itself and the encountering stretch that tunnels through to start and end in the leads .
the contribution is calculated by summing over all enclosed periodic orbits , using the sum rule ( [ in8 ] ) .
remember , that when allowing the encounter to move into the lead an additional time integral occurred measuring the part of the stretch that lies still inside the system .
now there are two time integrations representing the amount of the encounter which is cut short in each of the two leads . the first time integral cancels @xmath104 , while the second essentially yields a factor @xmath164 $ ] . for the details of the calculation we refer to @xcite .
this part of the diagram then contributes @xmath165 so that this configuration of fig .
[ ana2 ] altogether yields @xmath166}&= & -\frac{4n_{\rm l}^2n_{\rm r}^2}{\left(n_{\rm l}+n_{\rm r}\right)^4}(1-\gamma)^2 { \rm e}^{-p_2\te/\td}\nonumber\\ & & \times\left(1-{\rm e}^{-p_2\te/\td}\right),\end{aligned}\ ] ] where the ` ( s ) ' in the superscript denotes that the two links were removed along the same base trajectory from the same encounter stretch . one factor 2 in the last equation derives again from the mirror symmetry of this configuration explained after eq .
( [ li0.7 ] ) , the other from the fact that each of the two stretches can touch the opening at both ends .
furthermore , one of the links of the other base trajectory may also be shrunk so that the other encounter tunnels into the lead at one end .
this contribution is @xmath167}&= & \frac{8n_{\rm l}^2n_{\rm r}^2}{\left(n_{\rm l}+n_{\rm r}\right)^4}\frac{\gamma(1-\gamma)^3}{p_2 } { \rm e}^{-p_2\te/\td}\nonumber\\ & & \times\left(1-{\rm e}^{-p_2\te/\td}\right).\end{aligned}\ ] ] removing all 4 links , so that both encounter stretches on the base trajectories tunnel into the leads at both ends , is also possible .
the contribution is simply given by the square of ( [ li1.13 ] ) .
we can repeat this treatment for the aligned 3-encounter of fig .
[ ana3 ] . if one encounter stretch tunnels into the leads at both ends we obtain @xmath168}&= & \frac{8n_{\rm l}^2n_{\rm r}^2}{\left(n_{\rm l}+n_{\rm r}\right)^4}\frac{\gamma(1-\gamma)^4}{p_3 } \left(1-{\rm e}^{-p_3\te/\td}\right),\nonumber\\\end{aligned}\ ] ] where , because of the alignment and proximity of the two encounter stretches , the periodic orbit and the other encounter stretch must be backreflected at the tunnel barriers . allowing the ends of the second encounter stretch to progressively tunnel through into the leads as well
, we have , if it tunnels into the lead at one end @xmath169}&= & \frac{4\left(n_{\rm l}n_{\rm
r}^2+n_{\rm l}^2n_{\rm r}\right)}{\left(n_{\rm l}+n_{\rm r}\right)^3 } \frac{\gamma^2(1-\gamma)^3}{p_3}\nonumber\\ & & \times\left(1-{\rm e}^{-p_3\te/\td}\right)\end{aligned}\ ] ] and at both ends @xmath170}= \frac{2n_{\rm l}n_{\rm r}}{\left(n_{\rm l}+n_{\rm r}\right)^2}\frac{\gamma^3(1-\gamma)^2}{p_3 } \left(1-{\rm e}^{-p_3\te/\td}\right).\nonumber\\\end{aligned}\ ] ] such configurations can be also considered for the diagrams in fig.[ana4 ] and their contributions can be calculated analogously . to simplify the results , however
, we will later perform an expansion of the contributions in powers of the ehrenfest time and only retain terms up to linear order . as the results for the diagrams of fig .
[ ana4 ] are all of higher order in @xmath0 than the linear one , we will not focus on their explicit form here .      a further effect , and one that actually causes the independence of the conductance variance of the ehrenfest time for @xmath82 , are correlations during fringes near periodic orbits .
encounter fringes refer to regions where the two base trajectories which encounter the periodic orbit are correlated with each other , but are no longer correlated with the periodic orbit itself ; see fig .
[ ana5 ] for a schematic depiction .
as the two base trajectories leave the periodic orbit correlated with each other , we can consider the encounters with the periodic orbit to be aligned .
the encounter of both trajectories with the periodic orbit has length @xmath171 while the fringes , during which the base trajectories are correlated have lengths @xmath172 and @xmath173 where the subscripts refer to the fact that the stable and unstable distances between the base trajectories must be small for them to remain correlated . here
, we treat the case in which the encounter length is shorter than the periodic orbit . by this
we generalize the calculation of @xcite to @xmath174 .
there , it was shown that the contribution derived from ( [ li1.2 ] ) can be expressed as @xmath175}&= & \frac{2n_{\rm l}^2n_{\rm r}^2c^2\lambda}{\left(\pi\hbar\right)^2 \left(n_{\rm l}+n_{\rm r}\right)^4}\int_{1-b}^{b-1}{\rm d}s'{\rm d}u ' \int_0 ^ 1{\rm d}u\,u\nonumber\\ & & \times\int_{t_{\rm enc}}^\infty { \rm d}\tau_p { \rm e}^{-\left[\left(p_3-\gamma\right){t_{\rm enc}}+\gamma\tau_p\right]/\td } \nonumber\\ & & \times\left(\left|\frac{s'}{b-1}\right|^{\frac{p_2}{\lambda\td}}-1\right ) \left(\left|\frac{u'}{b-1}\right|^{\frac{p_2}{\lambda\td}}-1\right)\nonumber\\ & & \times\cos\left[\frac{c^2u\left(s'-u'\right)}{\hbar}\right]\end{aligned}\ ] ] where @xmath31 is again a classical constant of order unity .
the stable and unstable coordinates are defined as in @xcite : as correlation effects away from the periodic orbits get only important when the encountering orbits approach and leave this at classically close points one first characterize the correlation of all three orbit pieces ( i.e.  the periodic orbit and the two encountering orbits ) by the coordinate @xmath176 .
furthermore the difference between the stable and the unstable coordinates of the encountering orbits at the beginning and the end of the encounter is denoted by @xmath177 and @xmath178 .
these coordinates characterize the time the encountering orbits remain correlated before and after they are correlated with the periodic orbit , respectively . in ( [ li1.19 ] )
, we included the correct dwell times for the different parts of the trajectory : @xmath76 for the isolated links , @xmath94 for the fringes and @xmath179 for the encounter stretches correlated with the periodic orbit .
following the steps in @xcite , this contribution evaluates to @xmath180}&= & \frac{2n_{\rm l}^2n_{\rm r}^2}{\left(n_{\rm l}+n_{\rm
r}\right)^4 } \frac{p_2 ^ 2}{\gamma\left(2p_2-p_3\right)}\nonumber\\ & & \times\left({\rm e}^{-p_3\te/\td}-{\rm e}^{-2p_2\te/\td}\right).\end{aligned}\ ] ] without tunnel barriers , as shown in @xcite , further contributions can be obtained when the fringes start in the leads .
similarly , with tunnel barriers , we can have the fringes tunnel into the leads .
this leads to an additional time integral over the duration of the fringe that remains inside the system .
when one fringe tunnels into the lead , two links are removed and the contribution becomes @xmath181}&= & -\frac{2\left(n_{\rm l}^2n_{\rm r}+n_{\rm l}n_{\rm r}^2\right)}{\left(n_{\rm l}+n_{\rm r}\right)^3 } \frac{\gamma p_2}{\left(2p_2-p_3\right)}\nonumber\\ & & \times\left({\rm e}^{-p_3\te/\td}-{\rm e}^{-2p_2\te/\td}\right).\end{aligned}\ ] ] the ` ( s ) ' in the superscript refers to the fact that the links were removed from the same side of the encounter . when both fringes , and all four links are removed , we have @xmath182}&= & \frac{2n_{\rm l}n_{\rm r}}{\left(n_{\rm l}+n_{\rm r}\right)^2 } \frac{\gamma^3}{\left(2p_2-p_3\right)}\nonumber\\ & & \times\left({\rm e}^{-p_3\te/\td}-{\rm e}^{-2p_2\te/\td}\right).\end{aligned}\ ] ] with the tunnel barriers though , new possibilities also arise . for example , when a fringe tunnels into the lead , one of the trajectories could be backscattered , as for the 2-encounters previously . similarly to ( [ li0.4 ] ) , we obtain here @xmath183}&= & -\frac{8n_{\rm l}^2n_{\rm r}^2}{\left(n_{\rm l}+n_{\rm r}\right)^4 } \frac{(1-\gamma)p_2}{\left(2p_2-p_3\right)}\nonumber\\ & & \times\left({\rm e}^{-p_3\te/\td}-{\rm e}^{-2p_2\te/\td}\right).\end{aligned}\ ] ] furthermore we can move additionally one fringe into the lead yielding @xmath184}&= & \frac{4\left(n_{\rm l}n_{\rm r}^2+n_{\rm l}^2n_{\rm r}\right)}{\left(n_{\rm l}+n_{\rm r}\right)^3 } \frac{(1-\gamma)\gamma^2}{\left(2p_2-p_3\right)}\nonumber\\ & & \times\left({\rm e}^{-p_3\te/\td}-{\rm e}^{-2p_2\te/\td}\right).\end{aligned}\ ] ] if both fringes are partially backscattered at the leads , we have the contribution @xmath185}&= & \frac{8n_{\rm l}^2n_{\rm r}^2}{\left(n_{\rm l}+n_{\rm r}\right)^4}\frac{\gamma(1-\gamma)^2 } { \left(2p_2-p_3\right)}\nonumber\\ & & \times\left({\rm e}^{-p_3\te/\td}-{\rm e}^{-2p_2\te/\td}\right).\end{aligned}\ ] ]    next we can allow the periodic orbit itself , rather than the fringes , to touch the lead and be reflected from the tunnel barriers . the encounter stretches may then both tunnel into the lead , or one may be reflected .
we recall that the case in which both are reflected is already included in the survival probability of the encounter .
one fringe is therefore removed from the systems , while the fringe that remains could then be inside the system , tunnel into the lead or be partially backreflected at the opening .
if the periodic orbit touches both leads during the encounter , we return to the situation in the previous subsection , so here we only consider the cases in which at least some of the second fringe remains .
the removal of one fringe entirely implies the following changes in ( [ li1.19 ] ) : with only one fringe time left , we only have one of the factors in the third line of ( [ li1.19 ] ) . also the encounter time is replaced by an integration variable ( with a range from @xmath186 to @xmath187 ) which accounts for the amount of the encounter left after the periodic orbit is reflected from the tunnel barriers . performing , with these changes , again the steps in @xcite analogous to those following ( [ li1.19 ] ) above
, we obtain the contributions for the various different cases .
first , when the remaining fringe lies inside the system , we could have the encounter ( which touches the lead ) partially reflected @xmath188}&= & -\frac{4n_{\rm l}^2n_{\rm r}^2}{\left(n_{\rm l}+n_{\rm r}\right)^4}\frac{(1-\gamma)^2p_2}{p_3}\nonumber\\ & & \times\left(1-{\rm e}^{-p_3\te/\td}\right){\rm e}^{-p_2\te/\td},\end{aligned}\ ] ] where ` -1f-1l ' refers to the fact that one fringe and one link were removed .
likewise , both encounter stretches could tunnel directly into the lead so that two links are removed on the same side of the encounter @xmath189}&= & -\frac{\left(n_{\rm l}n_{\rm r}^2+n_{\rm l}^2n_{\rm r}\right)}{\left(n_{\rm l}+n_{\rm r}\right)^3 } \frac{(1-\gamma)p_2\gamma}{p_3}\nonumber\\ & & \times\left(1-{\rm e}^{-p_3\te/\td}\right){\rm e}^{-p_2\te/\td}.\end{aligned}\ ] ] the remaining fringe could also be paritally reflected from the tunnel barriers . for the encounter stretches
have the same two possibilities as above .
this leads to @xmath190}&= & \frac{8n_{\rm l}^2n_{\rm
r}^2}{\left(n_{\rm l}+n_{\rm r}\right)^4}\frac{(1-\gamma)^3\gamma}{p_3}\nonumber\\ & & \times\left(1-{\rm e}^{-p_3\te/\td}\right){\rm e}^{-p_2\te/\td } , \end{aligned}\ ] ] and @xmath191}&= & \frac{2\left(n_{\rm l}^2n_{\rm r}+n_{\rm l}n_{\rm r}^2\right)}{\left(n_{\rm l}+n_{\rm r}\right)^3 } \frac{(1-\gamma)^2\gamma^2}{p_3}\nonumber\\ & & \times\left(1-{\rm e}^{-p_3\te/\td}\right){\rm e}^{-p_2\te/\td}.\end{aligned}\ ] ] finally the remaining fringe can tunnel directly into the lead , giving the contributions @xmath192}&= & \frac{2\left(n_{\rm l}^2n_{\rm r}+n_{\rm l}n_{\rm r}^2\right)}{\left(n_{\rm l}+n_{\rm r}\right)^3 } \frac{(1-\gamma)^2\gamma^2}{p_3}\nonumber\\ & & \times\left(1-{\rm e}^{-p_3\te/\td}\right){\rm e}^{-p_2\te/\td } , \end{aligned}\ ] ] and @xmath193}&= & \frac{2n_{\rm l}n_{\rm r}}{\left(n_{\rm l}+n_{\rm r}\right)^2 } \frac{(1-\gamma)^2\gamma^3}{p_3}\nonumber\\ & & \times\left(1-{\rm e}^{-p_3\te/\td}\right){\rm e}^{-p_2\te/\td}.\end{aligned}\ ] ]      having treated all encounter configurations with periodic orbits where the encounter stretches are shorter than the enclosed periodic orbit we now turn to the corresponding contributions where the encounter stretches are longer than the enclosed periodic orbit . in this context only diagrams with stretches that both possess the same number of traversals around the enclosed periodic orbit yield a contribution , for a justification see @xcite . in this context
we consider that each encounter stretch has @xmath194 full windings around the enclosed periodic orbit .
we take here the two encounter times @xmath195 , @xmath196 and the corresponding primitive times @xmath197 .
these primitive times are again shorter than the enclosed periodic orbit and we can now also consider the different cases depicted in figs .
[ ana2][ana5 ] , and how they can be arranged .
we start with the case in which the two primitive encounter times , @xmath198 , do not overlap .
when @xmath199 , the corresponding diagram is depicted in fig .
[ ana2 ] . to explain the calculation we begin with eq .
( [ li1.2 ] ) for @xmath199 and perform the integrals over the links @xmath200 , and over @xmath105 where the latter integral leads to the factor @xmath201 : @xmath202}_{k=0}&= & \frac{n_{\rm
l}^2n_{\rm r}^2}{\left(n_{\rm l}+n_{\rm r}\right)^{4 } } \int_{-c}^c { \rm d}\boldsymbol{s}{\rm d}\boldsymbol{u } \frac{\th^2}{\omega^2}\frac{{\rm e}^{\frac{{\rm i}}{\hbar}\boldsymbol{s}\boldsymbol{u } } } { t_{{\rm enc},1 } t_{{\rm enc},2 } } \nonumber \\ & & \times\int_{t_{{\rm encs}}}^\infty   { \rm d}\tau_p\left(\tau_p - t_{{\rm encs}}\right ) \nonumber \\ & & \times { \rm e}^{-\left[\left(p_2-\gamma\right)t_{{\rm encs } } + \gamma \tau_p\right]/\td } .\end{aligned}\ ] ] to turn to @xmath203 , @xmath104 is then replaced only inside the integrand of the @xmath132-integral by @xmath198 , these times are shorter than @xmath132 by definition .
the @xmath104 before the @xmath132-integral compensate the overcounting of equivalent positions of the poincar surface of sections by the @xmath204-integrals and thus do not need to be altered .
the possibilities for placing these stretches around the periodic orbit are the same as they were for @xmath104 for @xmath199 .
this allows us to treat @xmath205 in essentially the same way we treated @xmath199 before . as @xmath206
, we reexpress the primitive encounter times in terms of @xmath104 and @xmath132 . the limits of the @xmath132-integration are also altered for @xmath203 : as the primitive encounters do not ( yet ) overlap , we have the condition @xmath207 so that the lower limit is now @xmath208 .
since we remove all @xmath194 complete windings of the periodic orbit from both encounter times , the shorter encounter must be @xmath209 so that the upper limit is @xmath210 .
finally @xmath211 is replaced by @xmath212 ( where in this context @xmath3 is defined as @xmath213 ) since during the primitive encounter stretches we have a total of @xmath214 correlated stretches while elsewhere we have @xmath215 .
all these replacements finally yield @xmath216}_{{\rm mt}}&= & \frac{n_{\rm l}^2n_{\rm r}^2}{\left(n_{\rm
l}+n_{\rm r}\right)^{4 } } \int_{-c}^c { \rm d}\boldsymbol{s}{\rm d}\boldsymbol{u } \frac{\th^2}{\omega^2}\frac{{\rm e}^{\frac{{\rm i}}{\hbar}\boldsymbol{s}\boldsymbol{u } } } { t_{{\rm enc},1 } t_{{\rm enc},2 } } \\ & & \times\int_{\frac{t_{{\rm encs}}}{2k+1}}^{\frac{t_{{\rm enc , min}}}{k } }   { \rm d}\tau_p\left[(2k+1)\tau_p - t_{{\rm encs}}\right]\nonumber\\ & & \times{\rm e}^{-\left[\left(p_{2k+2}-p_{2k+1}\right ) \left(t_{{\rm encs}}-2k\tau_p\right)+p_{2k+1}\tau_p\right]/\td } , \nonumber\end{aligned}\ ] ] where the additional index ` mt ' indicates that this contribution results from multiple traversals @xmath205 .
performing now the remaining integrals in the same way as for example in eqs .
( [ li1.4][li1.7 ] ) we obtain this contribution . as the full expression
is rather involved we just give here terms up to linear order in @xmath0 : @xmath217}_{\rm mt } & = & -\frac{2n_{\rm l}^2n_{\rm r}^2}{\left(n_{\rm l}+n_{\rm r}\right)^4}\sum_{k=1}^\infty \frac{\left(p_{2k+2}-p_{2k+1}\right)\te}{k^2(2k+1)\td } \nonumber \\ & & { } + { \mathcal o}\left(\te^2\right ) .\end{aligned}\ ] ] we also omitted here terms that are also nonzero for @xmath218 as they all cancel with the contributions below .
also in this case one link can be removed leading to a configuration where one encounter stretch touches the opening .
this implies that in ( [ li1.23.5 ] ) one factor @xmath219 is removed , the duration of the encounter stretch which touches the tunnel barrier is replaced by an integration variable integrated from @xmath186 to @xmath104 , and additionally a factor @xmath220 arises to account for the probability that all the @xmath215 orbital parts surrounding the periodic orbit are backreflected when they hit the lead .
this contribution then yields @xmath221}_{\rm mt}&= & \frac{4n_{\rm
l}^2n_{\rm r}^2}{\left(n_{\rm l}+n_{\rm r}\right)^4}\sum_{k=1}^\infty \frac{\left(p_{2k+2}-p_{2k+1}\right)\te}{k^2(2k+1)\td}\nonumber\\ & & { } + { \mathcal o}\left(\te^2\right ) . \end{aligned}\ ] ] shrinking two links connecting two _ different _ encounter stretches to the opening , we obtain ( following the same steps as just described also for the other @xmath104 ) no contribution which is linear in @xmath222 .
shrinking two links connecting the _ same _ encounter we need again to introduce in ( [ li1.23.5 ] ) two additional time integrals as described before ( [ li1.13 ] ) .
this also finally yields zero contribution linear in @xmath0 , because the terms resulting from the two limits of the @xmath132-integral cancel .
the same also holds for contributions following .
similar contributions are also obtained in the other cases . for a generalized 3-encounter introduced around eq .
( [ li1.75 ] ) we obtain the integrals given in ( [ ap2 ] ) which provide the contribution @xmath223}_{\rm mt}&= & \frac{2n_{\rm l}^2n_{\rm r}^2}{\left(n_{\rm l}+n_{\rm r}\right)^4}\sum_{k=1}^\infty \frac{\te}{k(k+1)\td } \\ & & \times \left(2p_{2k+2}-p_{2k+1}-p_{2k+3}\right ) + { \mathcal o}\left(\te^2\right ) , \nonumber\end{aligned}\ ] ] when the encounters lie inside the system and @xmath224}_{\rm mt}&= & -\frac{4n_{\rm l}^2n_{\rm r}^2}{\left(n_{\rm l}+n_{\rm r}\right)^4}\sum_{k=1}^\infty \frac{\te}{k(k+1)\td } \\ & & \times\left(2p_{2k+2}-p_{2k+1}-p_{2k+3}\right ) + { \mathcal o}\left(\te^2\right ) , \nonumber \end{aligned}\ ] ] when one link is removed .
additionally we can shrink both links connecting the 3-encounter to the opening on one side .
this amounts to removing two links , and introducing an additional integral over the part of the encounter that remains inside the system , in ( [ ap2 ] ) , yielding @xmath225}_{\rm mt}&= & \frac{2\left(n_{\rm l}n_{\rm r}^2+n_{\rm l}^2n_{\rm r}\right)}{\left(n_{\rm l}+n_{\rm r}\right)^3 } \sum_{k=1}^\infty\frac{\te}{k(k+1)\td } \\ & & \times\left(2p_{2k+2}-p_{2k+1}-p_{2k+3}\right ) + { \mathcal o}\left(\tau_e^2\right ) .
\nonumber\end{aligned}\ ] ] again zero contribution is obtained when more links are removed .    for a non - aligned 3-encounter , when the encounter lies inside the system we obtain from eqs .
( [ ap31],[ap32 ] ) @xmath226}_{\rm mt}&= & \frac{4n_{\rm l}^2n_{\rm r}^2}{\left(n_{\rm
l}+n_{\rm r}\right)^4}\sum_{k=1}^\infty \left[\frac{\left(p_{2k+3}-p_{2k+2}\right)\te}{(k+1)(2k+1)\td}-\right.\nonumber\\ & & \left.-\frac{\left(p_{2k+2}-p_{2k+1}\right)\te}{k(2k+1)\td}\right ] + { \mathcal o}\left(\te^2\right ) \end{aligned}\ ] ] and when one link is removed : @xmath227}_{\rm mt}&= & -\frac{8n_{\rm l}^2n_{\rm r}^2}{\left(n_{\rm l}+n_{\rm r}\right)^4}\sum_{k=1}^\infty \left[\frac{\left(p_{2k+3}-p_{2k+2}\right)\te}{(k+1)(2k+1)\td}\right.\nonumber\\ & & + \frac{\left(2p_{2k+2}-p_{2k+1}-p_{2k+3}\right)\te}{2k(k+1)\td}\nonumber\\ & & \left.-\frac{\left(p_{2k+2}-p_{2k+1}\right)\te}{k(2k+1)\td}\right ] + { \mathcal o}\left(\te^2\right).\end{aligned}\ ] ] finally we treat the configuration in which the encounter stretches overlap at both ends . when @xmath199 this configuration is depicted in fig .
[ ana4](b ) . with multiple traversals of the periodic orbit , the contribution given in eq .
( [ ap4 ] ) provides , when the stretches lie inside the system , @xmath228}_{\rm mt}&= & -\frac{2n_{\rm l}^2n_{\rm r}^2}{\left(n_{\rm l}+n_{\rm r}\right)^4}\sum_{k=1}^\infty \frac{\left(p_{2k+3}-p_{2k+2}\right)\te}{(k+1)^2(2k+1)\td}\nonumber\\ & & + { \mathcal o}\left(\te^2\right)\end{aligned}\ ] ] and when one link is removed @xmath229}_{\rm mt}&= & \frac{4n_{\rm l}^2n_{\rm
r}^2}{\left(n_{\rm l}+n_{\rm r}\right)^4}\sum_{k=1}^\infty \frac{\left(p_{2k+3}-p_{2k+2}\right)\te}{(k+1)^2(2k+1)\td}\nonumber\\ & & + { \mathcal o}\left(\te^2\right).\end{aligned}\ ] ] no further contributions are obtained when taking into account fringes as in fig .
[ ana5 ] . this can be checked from eq .
( [ li1.19 ] ) by adjusting the limits of the @xmath132-integral appropriately .
then , the terms from the two limits of the @xmath132-integral cancel to linear order in @xmath0 .      having calculated all the contributions we can now sum them to obtain the overall contributions to the conductance variance , to leading order in the number of open channels in the leads .
first we can check that all the terms after sec .
[ discrete ] are zero for @xmath122 , so that the rmt result ( [ li22.1 ] ) is preserved .
the case that is especially interesting for the comparison with the numerics is the contribution to the conductance variance when @xmath230 , for arbitrary @xmath3 and of linear order in @xmath0 .
we therefore consider for @xmath231 , in the presence of time reversal symmetry , @xmath232+{\alpha}\frac{\te}{\td}.\end{aligned}\ ] ] we obtain @xmath233 \end{aligned}\ ] ] by summing the contributions linear in @xmath222 from eqs .
( [ li0.3]-[li0.91 ] , [ li1 ] , [ li1.76]-[li1.8 ] , [ li1.10]-[li1.12 ] , [ li1.14]-[li1.16 ] , [ li1.20]-[li1.2297 ] , [ li1.24]-[li1.32 ] ) . here
@xmath234 denotes the polylogarithmic function , which arises when performing the @xmath194-summations above .
the function @xmath235 is shown as the solid ( red ) line in fig .
[ condfluc ] . having derived the complete set of contributions , we now discuss the main terms that lead to the two extrema in fig .
[ condfluc ] .
first we note that the contributions from multiply - traversed periodic orbits are quite small .
when @xmath82 they provide no contribution linear in @xmath0 , as can be checked from eqs .
( [ li1.24]-[li1.32 ] ) . for most of the contributions for @xmath85
the index @xmath194 of the @xmath236 s in the sum is much larger than one . in this case
the differences of the @xmath236 s in eqs .
( [ li1.24]-[li1.32 ] ) tend to one and thus the contributions from these equations tend to zero .
then the result is almost equal to the vanishing value for @xmath82 .    in fig .
[ condfluc ] , the maximum lies in a region of small @xmath3 . in this regime ,
most of the diagrams approximately cancel each other as can be seen by expanding their contributions in a taylor series for small @xmath3 : the ones from diagrams containing encounters inside the system [ e.g.  the contribution ( [ li0.3 ] ) ] are canceled by the contributions obtained when one and two links are shrunk [ e.g.  the contributions ( [ li0.4][li0.7 ] ) ] . however , this cancellation does not hold for the 3-encounter with the periodic orbit ( see fig .
[ ana3 ] ) . in this case
we have two possibilities : first that the encounter lies inside the system or touches the lead at one end [ see eqs .
( [ li0.8][li0.91 ] ) ] . in total
this leads to an ehrenfest - time dependent factor of @xmath237 .
additionally , we also need to take into account the contributions from encounters touching the leads on both sides leading to an increasing contribution proportional to @xmath238 with a larger prefactor than the first one [ see eqs .
( [ li1.18][li1.16 ] ) ] . together
these cause the first peak in fig .
[ condfluc ] .    in the case of the dip
the factor @xmath108 , i.e.  the probability for the particle to be backreflected at the opening , is quite small .
this implies that we get the main contributions from diagrams with encounters inside the system ( i.e.  the ones also obtained for @xmath82 ) depicted in figs .
[ ana1][ana5 ] .
these together yield a negative contribution .
we numerically confirm our main analytical prediction , eqs .
( [ li22.9 ] ) and ( [ li23 ] ) .
the model we use is the open kicked rotator with time - dependent hamiltonian  @xcite @xmath239 with @xmath240 the free flight time .
depending on the kicking strength @xmath241 the dynamics changes from integrable for @xmath242 to fully chaotic for @xmath243 . in the latter regime ,
local exponential instability is characterized by the lyapunov exponent @xmath244     as a function of @xmath245 for various values of @xmath3 .
the straight lines are linear fits constrained to go to @xmath246 at @xmath122
. their slope give the parameter @xmath19 defined in eqs .
( [ li22.9 ] ) and ( [ li23]).,scaledwidth=47.0% ]    the quantities @xmath247 and @xmath248 are introduced to break the hamiltonian s two symmetries @xcite , and investigate different symmetry classes with ( @xmath247 or @xmath249 ) and without ( @xmath250 ) time reversal symmetry .
the ehrenfest time in this system is determined , within a constant of little relevance , by @xmath251 with @xmath252 the size of the hilbert space , determined by the quantization of the hamiltonian via discretization of the coordinates as , e.g. @xmath253 , @xmath254 . here
@xmath44 is a system - dependent constant of order one that is of classical origin , and as such does not depend on @xmath252 . a quantum representation of the hamiltonian ( [ num0 ] )
is then provided by the unitary @xmath255 floquet matrix @xmath256 , giving the time evolution for one iteration of the map defined by @xmath257 in a time interval @xmath258 . for our specific choice of the kicked rotator ,
the floquet operator has matrix elements @xmath259 } \\ & \times & \sum_m e^{2 \pi i m(l - l')/m }     e^{-(imk/4\pi ) \cos(2\pi ( m+m_0)/m ) } \nonumber\end{aligned}\ ] ] with @xmath260 and @xmath261 .
transport can finally be investigated once absorbing phase - space strips are introduced to model contacts to leads .
this is achieved by means of projection operators @xmath134 . in our case of tunnel - coupled leads ,
the latter are @xmath262 diagonal matrices with entries @xmath263 , assuming that the system is coupled to all @xmath264 channels with the same transparency @xmath265 .
the scattering matrix is finally defined as  @xcite @xmath266^{-1 } u p^t \ , .\ ] ] the kicked rotator model is particularly performant to investigate the semiclassical limit with @xmath267 , as it allows for a rather large variation of the system size and hence the ehrenfest time .
previous such investigations of the open kicked rotator are presented e.g. in refs .
@xcite , and we refer the reader to these works for further details of the model .
the variance of the conductance is calculated by varying the quasienergies @xmath268 at given lead position .
averages of @xmath16 are further performed with different lead positions .
we first show in fig .
[ condflucraw ] the behavior of @xmath11 as a function of @xmath269 .
we clearly see @xmath3-dependent behaviors , as @xmath16 increases for small @xmath3 and decreases for larger @xmath3 until it becomes independent of @xmath0 for @xmath82 .
we extract the linear slope of the @xmath11 vs. @xmath269 curve via a linear fit , setting @xmath122 where the curves intersect the rmt universal value @xmath270
@xcite . further dividing this slope by @xmath3 gives us @xmath19 .
there is still an uncertainty in @xmath19 due to an uncertainty in the precise value of the lyapunov exponent  the latter has been found numerically to deviate in open systems from its exact value ( [ num1 ] ) , which has been attributed to finite - time effects  and because of the arbitrary constant @xmath44 of order one in the definition of @xmath0 .
we remove this uncertainty by forcing numerical and analytical data to agree for @xmath271 .
once this is done , there is no free parameter left .
we compare the so obtained numerical values for @xmath19 with the analytical prediction of eq .
( [ li23 ] ) in fig .
[ condfluc ] .
the excellent agreement between numerical data and the analytical curve fully confirms our theory .
we also checked , but do not show , that the linear @xmath0-dependence of @xmath272 is halved when time - reversal symmetry is broken .
in this article we determined the dependence of the conductance variance of a chaotic cavity with tunnel barriers on the tunneling probability @xmath3 .
in particular we find an ehrenfest - time dependence for the general case @xmath273 of non - perfect coupling .
we focused on the contribution which is linear in the ehrenfest time and at leading order in the total number of open channels @xmath83 .
we predict a nonmonotonous sinusoidal behavior with one maximum at @xmath274 and a ( roughly 50% ) deeper minimum at @xmath275 .
this analytical result was derived semiclassically by systematically considering a rather large number of possible configurations , but the general behavior derives from a much smaller set .
we discussed the main contributions which lead to each of the two peaks to provide a better intuitive understanding of the overall structure .
finally , we compared these analytical predictions with numerical simulations performed for the kicked rotor ( in the chaotic regime ) .
there was good agreement between the two , with the analytical curve within the error bars of all the numerically determined datapoints .
although , even at leading order , we treated a large number of possible semiclassical diagrams , they could be constructed , as we showed in this article , in a controlled way from the diagrams that exist when the tunnel barriers are absent .
this calculation could therefore possibly be extended to higher order contributions in @xmath276 and @xmath0 , but a more natural next step would be to extend the ehrenfest - time dependence of the ( leading order in @xmath83 ) results for the full counting statistics of transport moments @xcite to include tunnel barriers in the leads and to account for the extra diagrammatic possibilities that then arise .    on the numerical side ,
we emphasize that the agreement observed in this article is the first observed for the conductance properties of systems with tunnel barriers . neither the results for the conductance nor for shot noise obtained in @xcite could yet be confirmed by numerical simulations .
further numerical investigation in this area would therefore be highly desirable .
finally it would be interesting to check the observed phenomena experimentally .
although ehrenfest - time effects in antidot superlattices were observed more than ten years ago @xcite , none of the other ehrenfest - time dependencies predicted for chaotic systems could be checked so far experimentally .
this has been mainly due to the small range through which the system size can be varied in quantum dots .
however , with the additional parameter of the tunneling probability and in view of the expected double peak structure , the prediction in this article could be amenable to experimental verification .
this work was funded by the deutsche forschungsgemeinschaft ( research unit for 760 ) ( dw and kr ) and by the alexander von humboldt foundation ( jk ) .
pj has been partially supported by the nsf under grant no dmr-0706319 and by the swiss center of excellence manep .
it is a great pleasure to thank c.  petitjean for discussions and for drawing our attention to this problem .
furthermore we thank r.  whitney for discussions at the early stage of this project .
here we give the necessary integrals for obtaining the contributions from diagrams containing enclosed periodic orbits that are surrounded @xmath194-times by each encounter stretch .
in these expressions the link and @xmath105-integrals in eq .
( [ li1.2 ] ) have already been performed , but not yet the @xmath132- and the @xmath277-integrals .
these expressions are especially useful because they not only allow one to obtain the contributions from encounters inside the system but ( via integrations with respect to the corresponding encounter times ) also the contributions when encounters touch the openings .    the corresponding expression obtained for two independent 2-encounters
was already given in the main text in eq .
( [ li1.23.5 ] ) . in the case of a generalized 3-encounter
we obtain @xmath278}_{\rm mt}&= & \frac{2n_{\rm l}^2n_{\rm r}^2}{\left(n_{\rm l}+n_{\rm r}\right)^4 } \int_{-c}^c{\rm d}\boldsymbol{s}{\rm d}\boldsymbol{u } \int_{\frac{t_{{\rm enc , max}}}{k+1}}^{\frac{t_{{\rm enc , min}}}{k}}{\rm d}\tau_p\frac{\th^2}{\omega^2}\nonumber\\ & & \times\frac{{\rm e}^{\frac{{\rm i}}{\hbar}\boldsymbol{s}\boldsymbol{u } } } { t_{{\rm enc},1}t_{{\rm enc},2}}\left(t_{\rm enc , max}-t_{\rm enc , min}\right)\nonumber\\ & & \times{\rm e}^{-\left(p_{2k+3}-p_{2k+2}\right)\left(t_{\rm enc , min}-k\tau_p\right)/\td}\\ & & \times{\rm e}^{-\left[\left(p_{2k+2}-p_{2k+1}\right)\left(t_{\rm enc , max}-k\tau_p\right)+p_{2k+1}\tau_p\right]/\td}\nonumber.\end{aligned}\ ] ] for a non - aligned 3-encounter , similar to ( [ li1.4],[li1.5 ] ) , it is given first for @xmath279 by @xmath280}_{\rm mt(1)}&= & \frac{2n_{\rm l}^2n_{\rm r}^2}{\left(n_{\rm l}+n_{\rm r}\right)^4}\frac{2\td}{\left(p_{2k+3}+p_{2k+1}-2p_{2k+2}\right)}\nonumber \\ & & \times \int_{-c}^c{\rm d}\boldsymbol{s}{\rm d}\boldsymbol{u } \frac{\th^2}{\omega^2 } \frac{{\rm e}^{\frac{{\rm i}}{\hbar}\boldsymbol{s}\boldsymbol{u } } } { t_{{\rm enc},1}t_{{\rm enc},2 } } \int_{\frac{t_{{\rm encs}}}{2k+1}}^{\frac{t_{{\rm enc , min}}}{k } } { \rm d}\tau_p\nonumber\\ & & \times{\rm e}^{-p_{2k+1}\tau_p/\tau_d}\nonumber\\&&\times\left[{\rm e}^{-\left(p_{2k+2}-p_{2k+1}\right ) \left(t_{{\rm encs}}-2k\tau_p\right)/\td}\right.\nonumber\\ & & { } -{\rm e}^{-\left(p_{2k+3}-p_{2k+2}\right ) \left(t_{\rm enc , min}-k\tau_p\right)/\td } \nonumber\\ & & \left.\times{\rm e}^{-\left(p_{2k+2}-p_{2k+1}\right ) \left(t_{\rm enc , max}-k\tau_p\right)/\td}\right ] , \end{aligned}\ ] ] and second for @xmath281 by @xmath282}_{\rm mt(2)}&= & \frac{2n_{\rm l}^2n_{\rm
r}^2}{\left(n_{\rm l}+n_{\rm r}\right)^4}\frac{2\td}{\left(p_{2k+3}+p_{2k+1}-2p_{2k+2}\right)}\nonumber \\ & & \times \int_{-c}^c{\rm d}\boldsymbol{s}{\rm d}\boldsymbol{u } \frac{\th^2}{\omega^2 } \frac{{\rm e}^{\frac{{\rm i}}{\hbar}\boldsymbol{s}\boldsymbol{u } } } { t_{{\rm enc},1}t_{{\rm enc},2 } } \int_{\frac{t_{{\rm enc , max}}}{k+1}}^{\frac{t_{{\rm encs}}}{2k+1 } } { \rm d}\tau_p \nonumber\\ & & \times\left[{\rm e}^{-\left ( 2p_{2k+2}-p_{2k+3}\right)\tau_p/\td}\right.\nonumber\\ & & \times{\rm e}^{-\left(p_{2k+3}-p_{2k+2}\right ) \left(t_{{\rm encs}}-2k\tau_p\right)/\td } \nonumber\\ & & { } -{\rm e}^{-\left(p_{2k+3}-p_{2k+2}\right ) \left(t_{\rm enc , min}-k\tau_p\right)/\td } \nonumber \\ & & \times{\rm e}^{-\left(p_{2k+2}-p_{2k+1}\right ) \left(t_{\rm enc , max}-k\tau_p\right)/\td } \nonumber \\ & & \left.\times{\rm e}^{-p_{2k+1}\tau_p/\td}\right].\end{aligned}\ ] ] in the case of the two encounters overlapping at both ends the contribution is obtained to be @xmath283}_{\rm mt}&=&\frac{2n_{\rm l}^2n_{\rm r}^2}{\left(n_{\rm l}+n_{\rm r}\right)^4 } \int_{-c}^c{\rm d}\boldsymbol{s}{\rm d}\boldsymbol{u } \int_{\frac{t_{{\rm enc , max}}}{k+1}}^{\frac{t_{{\rm encs}}}{2k+1 } } { \rm d}\tau_p \frac{\th^2}{\omega^2 } \nonumber\\ & & \times\frac{{\rm e}^{\frac{{\rm i}}{\hbar}\boldsymbol{s}\boldsymbol{u } } } { t_{{\rm encs}}}\left[t_{{\rm encs}}-\left(2k+1\right)\tau_p\right]\nonumber\\ & & \times{\rm e}^{-\left(p_{2k+3}-p_{2k+2}\right ) \left(t_{{\rm encs}}-2k\tau_p\right)/\td}\nonumber\\ & & \times{\rm e}^{-\left(2p_{2k+2}-p_{2k+3}\right)\tau_p/\td}.\end{aligned}\ ] ]
here we generalize eq .
( [ li22.1 ] ) to the case of different @xmath284 for the different lead modes .
factors like @xmath285 are then replaced by a sum over the @xmath284 with respect to the @xmath83 open channels .
an analogous replacement is made for the @xmath211 for @xmath286 . in order to keep the notation compact
we introduce @xmath287 for example , the contribution from the diagram in fig .
[ ana1](a ) , previously given in ( [ li0.3 ] ) , becomes @xmath288 } = \frac{g_{(1)}^{2}g_{(2)}^{2}g_{2}^{2}}{g_{1}^{6}}.\ ] ] where the ehrenfest - time dependence is the same as in ( [ li0.3 ] ) although we set @xmath122 here .
the first two terms in the numerator result from the channel summations .
the third term replaces the factor @xmath289 in ( [ li0.3 ] ) and the denominator replaces the previous @xmath290 .    to keep the contributions from configurations where one link connecting the encounter to the opening is removed in a compact form , we define @xmath291 considering again the corresponding contribution from fig.[ana1](a ) , previously given in ( [ li0.4 ] )
, we have @xmath292 } = -\frac{2\left[h_{(1),2}g_{(2)}+g_{(1)}h_{(2),2}\right]g_{2 } g_{(1)}g_{(2)}}{g_{1}^{5}}.\ ] ] the first term derives from when the link connecting the encounter to lead 1 is shrunk while the second term corresponds to when the link connecting an encounter to lead 2 is removed . compared to ( [ ap6 ] ) , the factor @xmath293 takes into account that the particle is entering the system in a certain channel and returns to the same channel after traversing a link .    removing two links from different encounters
likewise leads to @xmath294 } = \frac{4g_{(1)}h_{(1),2}h_{(2),2}g_{(2)}}{g_{1}^{4}}.\ ] ] for the diagram in fig .
[ ana1](a ) we can also remove two links from the same encounter , and for this we define @xmath295 so that we obtain @xmath296}\!&=&\ ! -\frac{\left[i_{(1),1}g_{(2)}^{2}+g_{(1)}^{2}i_{(2),1}\right]g_{2}}{g_{1}^{4 } } , \\
\label{ap8.3 } { \rm var}\,g^{[{\ref{ana1}(a)-3{\rm l}}]}\!&=&\ !
\frac{2\left[g_{(1)}h_{(1),2}i_{(2),1}+i_{(1),1}h_{(2),2}g_{(2)}\right]}{g_{1}^{3 } } , \nonumber\\ \\
\label{ap8.4 } { \rm var}\,g^{[{\ref{ana1}(a)-4{\rm l}}]}\!&=&\ !
\frac{i_{(1),1}i_{(2),1}}{g_{1}^{2 } } .\end{aligned}\ ] ] for the diagram in fig .
[ ana1](b ) we can not shrink two links attached to the same encounter ( or more than two links ) and because of the way the encounters are arranged when we shrink two links we obtain @xmath297 } & = & \frac{2g_{(1)}h_{(1),2}h_{(2),2}g_{(2)}}{g_{1}^{4 } } \\ & & { } + \frac{h_{(1),2}^{2}g_{(2)}^{2}+g_{(1)}^{2}h_{(2),2}^{2}}{g_{1}^{4 } } .
\nonumber\end{aligned}\ ] ] using the definitions above , this is equal to ( [ ap8.1 ] ) plus the results from ( [ ap8.2][ap8.4 ] ) so that time reversal symmetry still leads simply to a factor of 2 .    in an analogous manner , we can show that each diagram in fig .
[ ana2 ] gives twice the contribution of the diagram in fig .
[ ana1](b ) while for the diagrams in fig .
[ ana3 ] we obtain @xmath298 } & = & -\frac{2g_{(1)}^{2}g_{(2)}^{2}g_{3}}{g_{1}^{5}},\\ \label{ap8.7 } { \rm var}\,g^{[{\ref{ana3}-1{\rm l } } ] } & = &   \frac{4\left[g_{(1)}h_{(1),3}g_{(2)}^{2}+g_{(1)}^{2}h_{(2),3}g_{(2)}\right ] } { g_{1}^{4}},\nonumber\\   \\ \label{ap8.8 } { \rm var}\,g^{[{\ref{ana3}-2{\rm l } } ] } & = &   \frac{2\left[i_{(1),2}g_{(2)}^{2}+g_{(1)}^{2}i_{(2),2}\right]}{g_{1}^{3 } } .\end{aligned}\ ] ] summing all these contributions , we obtain the rmt result @xcite for the leading order in @xmath83 contribution to the conductance variance @xmath299 , \end{aligned}\ ] ] for systems without time reversal symmetry , and twice this result for those with . here
the notation introduced in @xcite @xmath300 was used .
here we calculate the shot noise power .
it can be written as @xcite @xmath301\right\rangle - \left\langle { \rm tr}\left[\left(\mathbf{tt}^{\dagger}\right)^2\right]\right\rangle.\ ] ] the average conductance @xmath20 has previously been calculated @xcite , and the first few terms for systems with time reversal symmetry can be written as @xmath302 using the notation of appendix [ diffprobs ] .
the second two terms derive from the diagrams in fig .
[ fig : orbt ] so that without time reversal symmetry only the first term in ( [ gdiffprobeq ] ) remains .
this is included as @xmath303 for systems without time reversal symmetry and @xmath304 for those with time reversal symmetry . the function @xmath305 is given semiclassically in terms of 4 trajectories @xmath306 where the main difference from the conductance variance ( [ condvartraj ] ) is that the trajectories connect different channels and the sum is unrestricted .
we start with the diagram from fig .
[ trajectoryquad ] , where since the partner trajectories cross over in the encounter they automatically travel from and to the correct channels and all four channels are unrestricted .
this means that this structure is now lower order ( in inverse channel number ) than it was for the conductance variance ( it was there contained by taking into account the diagram in fig .
[ ana1](a ) and removing two links connected to the same encounter ) . for the diagram in fig .
[ trajectoryquad ] , we therefore have @xmath307}&=&-\frac{g_{(1)}^{2}{g_{(2)}}^{2}g_2}{g_{1}^4 } \nonumber \\ & & { } + \frac{2\left[g_{(1)}h_{(1),2}g_{(2)}^2+g_{(1)}^2h_{(2),2}g_{(2)}\right]}{g_{1}^3 } \nonumber \\ & & { } + \frac{g_{(1)}^{2}i_{(2),1}+i_{(1),1}g_{(2)}^{2}}{g_{1}^{2}}.\end{aligned}\ ] ] this contribution , along with its ehrenfest - time dependence , was previously calculated in @xcite .        moving to the next order term , which has recently been calculated using rmt @xcite , we again simply need to look at the corresponding structures and consider the possible ways of shrinking links . without time reversal symmetry
there are no possible structures at this order , so we therefore consider systems with time reversal symmetry .
the structures which contribute are depicted in fig .  2 of @xcite .
there is also , in fig .
2(d ) there , a quadruplet involving two independent pairs , one of which is simply a diagonal pair , while the other involves a single 2-encounter , as in fig .
[ fig : orbt ] of this paper . because the start ( or end ) channels must coincide @xmath112 ( or @xmath113 )
there is a subtlety when we consider shrinking one link .
if we shrink the link on the left ( or right ) side of the 2-encounter so that it moves into the lead where the diagonal pair emanate ( or terminate ) we actually have the same case as when the 3-encounter from fig .
2(c ) ( of @xcite ) moves into the lead .
as such , we include this case there , leaving a contribution of @xmath308}&=&-\frac{2\left[g_{(1)}^{2}i_{(2),1}+i_{(1),1}g_{(2)}^{2}\right]g_{2 } } { g_{1}^{4 } } \\ & & { } + \frac{2\left[g_{(1)}h_{(1),2}i_{(2),1}+i_{(1),1}h_{(2),2}g_{(2)}\right ] } { g_{1}^{3}}. \nonumber\end{aligned}\ ] ] the diagram in fig .
2(a ) of @xcite is similar to the diagram in fig .
[ ana1](b ) here and we obtain @xmath309}&=&\frac{4g_{(1)}^{2}g_{(2)}^{2}g_{2}^{2}}{g_{1}^{6 } }   + \frac{8 g_{(1)}h_{(1),2}h_{(2),2}g_{(2)}}{g_{1}^{4 } } \nonumber \\ & & { } -\frac{8\left[g_{(1)}h_{(1),2}g_{(2)}^{2}+g_{(1)}^{2}h_{(2),2}g_{(2)}\right]g_{2 } } { g_{1}^{5 } } \nonumber \\ & & { } + \frac{2\left[g_{(1)}^{2}h_{(2),2}^{2}+h_{(1),2}^{2}g_{(2)}^{2}\right ] } { g_{1}^{4 } } .
\end{aligned}\ ] ] for the diagram in fig .
2(b ) of @xcite we obtain twice the result of fig .
[ ana2](a ) here [ or four times the result of fig .
[ ana1](b ) ] for the conductance variance . for the diagram in fig .
2(c ) of @xcite we likewise obtain twice the result of fig .
[ ana3](a ) here for the conductance variance . with all the next to leading order shot noise contributions
, we can combine them and indeed find the same result as in @xcite .
ehrenfest , z.  physik * 45 * , 455 ( 1927 ) . b.  v.  chirikov , f.  m.  izrailev , and d.  l.  shepelyansky , sov .
rev.sect .
c * 2 * , 209 ( 1981 ) .
i.  l.  aleiner , a.  i.  larkin , phys .
b * 54 * , 14423 ( 1996 ) .
o.  yevtushenko , g.  ltjering , d.  weiss , and k.  richter , phys .
rev.lett .  * 84 * , 542 ( 2000 ) .
i.  adagideli , phys .
b * 68 * , 233308 ( 2003 ) .
p.  w.  brouwer and s.  rahav , phys .
b * 74 * , 075322 ( 2006 ) .
whitney and ph .
jacquod , phys .
lett . * 94 * , 116801 ( 2005 ) ; ph .
jacquod and r.  s.  whitney , phys .
b * 73 * , 195115 ( 2006 ) .
d.  waltner and j.  kuipers , phys .
e * 82 * , 066205 ( 2010 ) .
jacquod and e.  v.  sukhorukov , phys .
lett .  * 92 * , 116801 ( 2004 ) .
j.  tworzydo , a.  tajic , h.  schomerus , p.  w.  brouwer , and c.  w.  j.  beenakker , phys .
lett .  * 93 * , 186806 ( 2004 ) .
a.  altland , p.  w.  brouwer , and c.  tian , phys .  rev .  lett .  * 99 * , 036804 ( 2007 ) .
c.  petitjean , p.  jacquod , and r.  s.  whitney , jetp lett .
* 86 * , 736 ( 2007 ) ; r.  s.  whitney , ph .
jacquod and c.  petitjean , phys .
b * 77 * , 045315 ( 2008 ) .
o.  agam , i.  l.  aleiner , and a.  i.  larkin , phys .
lett .  * 85 * , 3153 ( 2000 ) .
r.  s.  whitney and ph .
jacquod , phys .
lett .  * 96 * , 206804 ( 2006 ) .
p.  braun , s.  heusler , s.  mller , and f.  haake , j.  phys .
a * 39 * , l159 ( 2006 ) .
g.  berkolaiko , j.  m.  harrison , and m.  novaes , j.  phys .
a * 41 * , 365102 ( 2008 ) .
d.  waltner , j.  kuipers , and k.  richter , phys .
b * 83 * , 195315 ( 2011 ) .
i.  l.  aleiner and a.  i.  larkin , phys .
e * 55 * , r1243 ( 1997 ) . c.  tian and a.  i.  larkin , phys .
b * 70 * , 035305 ( 2004 ) .
p.  w.  brouwer , s.  rahav , and c.  tian , phys .
e * 74 * , 066208 ( 2006 ) .
d.  waltner , m.  gutirrez , a.  goussev , and k.  richter , phys .  rev
, 174101 ( 2008 ) .
m.  gutirrez , d.  waltner , j.  kuipers , and k.  richter , phys .  rev .
e * 79 * , 046212 ( 2009 ) .
b.  gutkin , d.  waltner , m.  gutirrez , j.  kuipers , and k.  richter , phys.rev .
e * 81 * , 036222 ( 2010 ) .
h.  u.  baranger , r.  a.  jalabert , and a.  d.  stone , chaos * 3 * , 665 ( 1993 ) . k.  richter and m.  sieber , phys .
lett .  * 89 * , 206801 ( 2002 ) .
k.  richter , _ semiclassical theory of mesoscopic quantum systems _
, springer , berlin , heidelberg ( 2000 ) .
p.  w.  brouwer and s.  rahav , phys .
b * 74 * , 085313 ( 2006 ) .
p.  w.  brouwer and c.  w.  j.  beenakker , j.  math .
phys .  * 37 * , 4904 ( 1996 ) .
m.  bttiker , phys .
lett .  * 57 * , 1761 ( 1986 ) .
y.  imry , in _ directions in condensed matter physics _ ,
g.  grinstein and g.  mazenko ( eds . ) ( singapore : world scientific ) ( 1986 ) .
d.  s.  fisher and p.  a.  lee , phys .
b * 23 * , 6851 ( 1981 ) .
s.  heusler , s.  mller , p.  braun , and f.  haake , phys .
lett .   *
96 * , 066804 ( 2006 ) ; s.  mller , s.  heusler , p.  braun , and f.  haake , new j.  phys .  * 9 * , 1 ( 2007 ) .
s.  mller , s.  heusler , p.  braun , f.  haake , and a.  altland , phys .
rev.lett . , 014103 ( 2004 ) ; phys .
e * 72 * , 046207 ( 2005 ) .
r.  s.  whitney , phys .
b * 75 * , 235404 ( 2007 ) . j.  tworzydo , a.  tajic , h.  schomerus , and c.  w.  j.  beenakker , phys.rev .
b * 68 * , 115313 ( 2003 ) .
h.  schomerus and ph .
jacquod , j. phys .
a * 38 * , 10663 ( 2005 ) .
f.  m.  izrailev , phys .
rep .  * 196 * , 299 ( 1990 ) .
y.  v.  fyodorov and h .- j .
sommers , jetp letters * 72 * , 422 ( 2000 ) ; r.  o.  vallejos and a.  m.ozorio de almeida , ann .
phys .  * 278 * , 86 ( 1999 ) .
m.  bttiker , phys .
* 65 * , 2901 ( 1990 ) .
j.  kuipers , j.  phys .
a * 42 * , 425101 ( 2009 ) .
j.  g.  g.  s.  ramos , a.  l.  r.  barbosa , and a.  m.  s.  macdo , phys.rev .
b * 78 * , 235305 ( 2008 ) .