the electron transport via resonant tunnelling in mesoscopic systems has been the subject of extensive theoretical research due to recent development in fabrication of small electronic devices and their potential applications .
some interest has been focused on the transport properties of a quantum dot ( qd ) under the influence of external time - dependent fields .
new effects have been observed and theoretically described , e.g. , photon - assisted tunnelling through small quantum dots , @xcite , photon - electron pumps @xcite and others . in most theoretical investigations
a qd placed between two leads was considered ( e.g. , refs .
@xcite ) and the current flowing through a qd under periodic modulation of the qd electronic structure or periodic ( nonperiodic ) modulation of the tunnelling barriers and electron energy levels of both ( left and right ) leads was calculated .
one of the important problems of the mesoscopic physics is the interference of the charge carries .
this interference appears when two ( or more ) transmission channels for electron tunnelling exist .
such possibility exists , e.g. in the electron transport through a qd embedded in a ring in the aharonov - bohm geometry and much theoretical interest has been paid to description of the phase coherence in this and related systems , e.g. refs . @xcite .
another experimental situation in which the interference may occur can be realized in the scanning tunnelling microscope ( stm ) .
the recent experimental and theoretical studies with a low - temperature stm of a single atom deposited on a metallic surface showed the asymmetric fano resonances in the tunnelling spectra , e.g. refs .
in the stm measurements tip probes the transmission of electrons either through the adsorbed atom or directly from the surface .
the transport of electrons through both channels leads to an asymmetric shape of the conductance curves which is typical behavior for the fano resonance resulting from constructive and destructive interference processes .
the quantum interference can be also observed in the mesoscopic system with multiple energy levels @xcite .
a model which incorporates a weak direct nonresonant transmission through a qd as well as the resonant tunneling channel was also considered in ref .
@xcite in the context of the large value of the transmission phase found in the experiment for the kondo regime of a qd @xcite .
a number of works has been devoted to the problem of the electron transport in the multiterminal qd systems and here we mention only a few of them . in ref .
@xcite the conductance of the n - lead system was considered showing that the kondo resonance at equilibrium is split into n-1 peaks . in ref .
@xcite an explicit form for the transmission coefficient in the electron transport through a qd connected with three leads was found .
the electron transport and shot noise in a multi - terminal coupled qd system in which each lead was disturbed by classical microwave fields were studied in ref .
multiterminal qd systems or magnetic junctions were also intensively investigated in context of the spin - dependent transport , e.g. refs .
a three - terminal qd system was studied in ref .
@xcite to measure of the nonequilibrium qd density of states ( splitting of the kondo resonance peak ) .
the cross correlations of the currents and the differential conductance of the qd coupled with three leads described by the infinite - u anderson hamiltonian were considered in ref .
the general formulation of the time - dependent spin - polarized transport in a system consisting of the resonant tunnelling structure coupled with several magnetic terminals was considered by zhao _
et al_. @xcite and as an application of this formalism the electron transport in a system with two terminals under an ac external field was investigated .        in all papers
mentioned above and in the studies relating to the electron transport through a qd with the additional ( bridge ) nonresonant transmission channels , the time - dependent external fields were not applied and the considered systems were driven out of the equilibrium only by means of a dc voltage bias ( see , however , ref .
@xcite )    in this paper , we address the issue of a qd coupled with three leads with additional , nonresonant coupling between leads driven out of equilibrium by means of a dc voltage bias and time - dependent external fields .
the qd is connected with three metal leads and one of these leads , say the left ( l ) lead , is coupled with the two remaining leads , say the first ( @xmath0 ) and the second ( @xmath1 ) right leads .
the possible stm experimental setup corresponding to our model system is presented in fig .
1 . in literature
, different theoretical approaches have been developed to treat the time - dependent electron transport in the mesoscopic systems .
the most popular nonequilibrium green s function method depends on the two time arguments and for time - dependent problems it is a rather difficult task to calculate them without any approximations ( e.g. beyond the wide - band limit ) .
therefore , in our treatment of the time - dependent problems we use the evolution operator which , as a rule , essentially depends on one time argument ( e.g. refs .
such an approach is especially well suited for the problem with time - dependent coupling between the qd and leads .
the outline of the paper is as follows . in sec .
ii we start with the model and method for the derivation of the qd charge and currents . in sec .
iii we present the results for the time - averaged currents and their derivatives with respect to the qd energy level position ( or equivalently , with respect to the gate voltage ) obtained for different time - dependence of the parameters characterizing the considered system .
we consider also the transient current characteristics in the case of the rectangular - pulse modulations imposed on the qd - lead barriers .
the last section presents the conclusions and in the appendix we give the short derivation of the evolution operator matrix elements needed in the qd charge and current calculations .
we consider a qd coupled through the tunneling barriers @xmath2 @xmath3 to three metal leads .
one of these leads , say the left lead ( l ) is coupled additionally with the remaining two leads , say the first and second right leads ( @xmath0 , @xmath1 ) by the tunneling barriers @xmath4 . in the following we will denote the wave vectors associated with the left lead by the letter @xmath5 and the wave vectors corresponding to the first and second right leads by the letters @xmath6 and @xmath7 , respectively .
the chemical potentials @xmath8 of the three metal leads may not be equal , and their difference is not necessarily small .
we write the hamiltonian of the considered system in the form @xmath9 , where @xmath10 @xmath11 the operators @xmath12 , @xmath13 are the annihilation ( creation ) operators of the electrons in the corresponding leads and the dot , respectively . for simplicity
the dot is characterised only by a single level @xmath14 and the intra - dot electron - electron coulomb interaction is neglected .
the neglect of coulomb interaction is quite reasonable in some systems and , as we are going to concentrate on the investigations of the influence of the third lead ( in comaprison with the qd - two leads system ) and the additional tunneling channels between the leads on the time - dependent transport , then in the first step ignoring the coulomb interaction should be justified .
we consider our mesoscopic system in the presence of external microwave ( mw ) fields which are applied to the dot and three leads . in most theoretical treatments of photon - assisted electron tunneling
it is assumed that the driving field equals the applied external field .
however , the situation is more complicated and the internal potential can be different from the applied potential @xcite .
one of the consequences will be , e.g. the asymmetry between the corresponding left and right sidebands @xcite .
the main feature of the time - dependent transport remains , however , unchanged and in our treatment as usual we assume that in the adiabatic approximation the energy levels of the leads and qd are driven with the frequency @xmath15 and the amplitudes @xmath16 ( @xmath17 ) , @xmath18 and read @xmath19 , @xmath20 , respectively .    the time - evolution of the qd charge and the current flowing in the system can be described in terms of the time - evolution operator @xmath21 given by the equation of motion ( in the interaction representation ) : @xmath22 with @xmath23 and @xmath24 where @xmath25 denotes the time ordering and the units such that @xmath26 have been chosen .
here we have assumed that the interactions between the qd and leads , as well as between the left and right leads are switched on in the distant past @xmath27 .
the qd charge is calculated according to the formula ( cf .
@xcite ) : @xmath28 here @xmath29 and @xmath30 denote the matrix elements of the evolution operator calculated within the basis functions containing the single - particle functions @xmath31 , @xmath32 , @xmath33 and @xmath34 corresponding to three leads and qd , respectively . @xmath35 and
@xmath36 represent the initial filling of the corresponding single - particle states .
the tunneling current flowing , e.g. from the left lead , @xmath37 , can be obtained using the time derivative of the total number of electrons in the left lead , @xmath38 ( in the following we put @xmath39 ) , where @xmath40 in the following only the matrix elements of the evolution operator present in eqs .
( 4 ) and ( 5 ) are required and they can be obtained solving the corresponding sets of coupled differential equations constructed according to eq .
( 3 ) with @xmath41 written as follows : @xmath42 where @xmath43 and @xmath44 correspond to @xmath31 , @xmath45 , @xmath33 or @xmath34 , respectively .    as the example ,
the matrix element @xmath29 required for the calculation of the first term of the qd charge , eq .
4 , can be obtained solving the following set of coupled differential equations @xmath46 the total number of coupled equations in this case is equal to @xmath47 , @xmath48 being the number of discrete wave vectors @xmath5 , @xmath6 and @xmath7 taken to perform the corresponding summation over the wave vectors . usually , the number @xmath48 equal to 100 - 200 is sufficient to achieve the desired accuracy of the calculations .
we have solved numerically this and other similar sets of the coupled differential equations needed in calculations of all matrix elements of the evolution operator present in eqs .
we have used this method for the special case of time - dependent couplings of the qd with leads and the couplings of the left lead with two right leads .
the set of eqs .
7 - 9 in the case of vanishing over - dot couplings between the left and right leads is greatly simplified and gives , e.g. for @xmath29 : @xmath49 where @xmath50 and @xmath51 is the fourier transform of the @xmath52-th lead density of states and @xmath53 .
similar simplifications occur in the calculations of other matrix elements of @xmath21 required in the formulas for @xmath54 and @xmath55 .
however , for the nonvanishing couplings @xmath56 and @xmath57 ( over - dot bridge between the left and right leads ) one has to solve the starting set of eqs .
much more analytical calculations can be done in the case of constant values of the tunneling matrix elements present in our model . in this case
the general equation satisfied by @xmath29 is derived in the appendix and under the wide band limit ( wbl ) approximation , e.g. @xcite this equation takes the simple form @xmath58 here @xmath59 , @xmath60 , @xmath61 being the bandwidth of the lead energy band ( @xmath62 ) and @xmath63 . in the appendix
we give the derivations of all functions needed for calculation of the qd charge and currents .
we assumed the simplified assumption that all tunneling matrix elements are independent of the wave vectors .
the interactions between the qd and leads are assumed to be equal between themselves and denoted by @xmath64 and the interactions between the left and two rights leads ( i.e. @xmath65 and @xmath66 ) corresponding to the over - dot tunneling channels are also equal one with another and denoted by @xmath67 .
it is easy to show that the first term of the general formula for the qd charge , eq .
4 , together with the solution of eq .
( 12 ) , @xmath68 , tends to zero for @xmath69 as @xmath70 .
the next terms of the qd charge formula can be calculated using the functions @xmath71 , @xmath72 and @xmath73 , eqs .
43,47 , being the solutions of the corresponding differential equations , eqs .
finally , the time - averaged qd charge is given by @xmath74 where @xmath75 @xmath76 where @xmath77 denotes the time - averaging and @xmath78 denotes the fermi function of the @xmath52-th ( @xmath79 lead .
noticing , that @xmath80 in order to calculate the current @xmath37 the functions @xmath81 , @xmath82 , @xmath83 and @xmath84 are required and they are given in the appendix in eqs .
37 , 44 , 49 .
after lengthly but straightforward calculations we obtain for the time averaged current leaving the left lead the following formula :    @xmath85 \,,\end{aligned}\ ] ]    where @xmath86 @xmath87 and @xmath88 denotes the bessel function .
the corresponding formula for the time - averaged current @xmath89 leaving @xmath90-lead , @xmath91 , can not be written in such symmetrical form as in eq .
18 , because the @xmath90-lead is coupled with @xmath92-lead only . for @xmath93
we have : @xmath94 \,,\end{aligned}\ ] ] where @xmath95 , @xmath96 , @xmath97 and @xmath98 for @xmath99 .
note , that the integrals present in the formula for the qd charge and currents , eqs .
17,18,21 , can be easily performed analytically and final algebraic expressions can be obtained . especially simple and transparent form can be given for the conductance @xmath100 , @xmath101 .
for example , @xmath102 reads as : @xmath103 where @xmath104^{-1 } \,,\end{aligned}\ ] ] @xmath105 analyzing eq .
22 one can find the origin of the asymmetric line shapes in differential conductance resulting from the interference of resonant and nonresonant tunneling paths .
for the case of @xmath106 we observe the lorentzian resonances localized at @xmath107 .
the amplitudes of these resonances are determined by the @xmath108-th order bessel functions calculated for the argument @xmath109 . for the case of nonvanishing @xmath110 ,
the resulting curve is a superposition of the lorentzian - like resonances and asymmetric parts weighed by the factors @xmath111 and @xmath112 , respectively .
the lorentzian - like resonance is centered at @xmath113 with the peak width at half maximum ( fwhm ) equal to @xmath114 and the maximum value equal to @xmath115 .
the asymmetric part of the differential conductance corresponding to the @xmath116th sideband is also centered in the same point with the distance between its maximum and minimum equal to @xmath117 and the absolute values of these extrema are equal to @xmath118 . for comparison , in the case of the qd
coupled with two leads the corresponding lorentzian - like part of the differential conductance corresponding to the @xmath116th sideband is centered at @xmath119 , with fwhm equal to @xmath120 and the maximum value equal to @xmath121 . knowing the explicit expressions for the currents one can check the following relations between different elements of the conductance matrix @xmath122 , e.g. @xcite .
current conservation implies @xmath123 , @xmath124 . on the other hand , @xmath125 only for @xmath126 . for other relations between the amplitudes @xmath18 and @xmath127 , @xmath128 , @xmath129 we have :    @xmath130    for @xmath131 , and    @xmath132    for @xmath133 .
we consider the qd coupled with three , say the left , the first right and second right metal leads with the additional over - dot ( bridge ) couplings between the left and right leads .
we present the results for the time - averaged currents and derivatives of the average current with respect to the qd energy level in the presence of external microwave fields which are applied to the dot and three leads , respectively .
the time - dependent currents are also calculated in the case when the periodic rectangular - pulse external fields are applied to each qd - lead barrier .
we have taken for @xmath110 the values comparable with @xmath134 and estimated @xmath134 ( assuming its @xmath5-independence , @xmath135 ) using the relation @xmath136 , where @xmath137 is the @xmath138-lead bandwidth and @xmath139 ( @xmath140 , @xmath141 was assumed ) . in our calculations we assumed @xmath142 ,
all energies are given in @xmath143 units , time in @xmath144 units , the current , its derivatives and frequency are given in @xmath145 and @xmath146 units , respectively .
firstly , we consider the time - averaged currents in the presence of the time - varying ( harmonic case ) external fields . here
we give the explicit formula for the averaged current , @xmath147 , performing the corresponding integrals in the general formula given in eq .
( 18 ) , and compare it with the current @xmath148 flowing in the system of the qd coupled with two leads only : @xmath149 \nonumber\\ & + & { \gamma \over \pi}{x(1 - 2x^2 ) \over ( 1 + 2x^2)^3 } \nonumber\\   & \times & \left .
\sum_k \left [ j_k^2   \left ( { \delta_d-\delta_l \over \omega } \right ) \ln(g_l)- j_k^2 \left ( { \delta_d-\delta_i \over \omega } \right ) \ln(g_i ) \right ] \right\},\end{aligned}\ ] ] where @xmath150 , @xmath151 and @xmath152 , whereas for @xmath148 we have : @xmath153 \nonumber\\ & + & { \gamma \over \pi}{x(1-x^2 ) \over ( 1+x^2)^3 } \nonumber\\   & \times & \sum_k \left [ j_k^2   \left ( { \delta_d-\delta_l \over \omega } \right ) \ln(g_l^{(2)})- j_k^2 \left ( { \delta_d-\delta_i \over \omega } \right ) \ln(g_r^{(2 ) } ) \right],\end{aligned}\ ] ] where @xmath154 and @xmath155 , @xmath156 . note , that @xmath147 consists of the two terms and each term is similar in its structure to the current flowing in the qd - two leads ( qd-2l ) system , @xmath148
. however , due to the interference of the charge carriers propagating along the different ways the arguments of the _ arctangens _ and _ logarytmic _ functions are different and the individual terms in @xmath147 and @xmath148 are weighted by different @xmath157-dependent factors .     against @xmath14 in the system of a qd coupled with three ( solid line ) or two ( broken line ) leads . the microwave field is applied to the qd with @xmath158 and @xmath159 , @xmath160 , @xmath161 , @xmath162 , @xmath163 . ]    in figure [ fig2 ] we compare the averaged values of the current flowing from the left lead in the systems in which the qd is coupled with three or two leads , the solid and broken curves , respectively .
the external microwave field is applied only to the qd and dc bias between the left and right leads is small in comparison with @xmath15 , @xmath18 and @xmath143 .
the coupling @xmath110 is assumed to be zero .
in such a case the sidebands on the current curve are clearly visible . the width of the corresponding peaks is smaller for the case of the qd coupled with two leads .
analyzing the expressions for @xmath164 and @xmath148 one can obtain ( for @xmath165 and @xmath161 ) the subsequent peaks in @xmath147 ( as functions of @xmath14 ) in the form @xmath166 with the fwhm equal @xmath167 . for comparison , in the case of the qd
coupled with two leads the corresponding peaks are described by the functions @xmath168 with the fwhm equal @xmath169
.     against @xmath14 in the system of a qd coupled with three ( thick lines ) or two ( thin lines ) leads .
the solid , broken and dotted curves correspond to @xmath170 and @xmath171 , respectively .
@xmath172 , @xmath173 , @xmath174 , @xmath159 . ]    next , in fig .
[ fig3 ] we analyze the influence of the over - dot coupling between the left and right leads . again , as in fig .
[ fig2 ] , we consider an external microwave field applied only to the qd with @xmath172 and @xmath175 .
the thin ( thick ) lines correspond to the qd coupled with two ( three ) leads .
we show @xmath147 for three values of the inter - leads coupling strength represented by the parameter @xmath176 .
for @xmath177 we have well defined sideband peaks as previously shown in fig .
next , we show the results for @xmath178 and @xmath179 . note , that the expression for @xmath147 ( and @xmath148 ) consists of three terms .
the first term depends on the difference @xmath180 and does not depend on the qd energy level @xmath14 .
the second term corresponds to the loretzian - type contribution to a given sideband ( it disappears for @xmath181 ) and the last term corresponds to an asymmetric contribution ( it disappears for @xmath182 ) , respectively .
this last term influences the sideband shape only for nonvanishing over - dot coupling between leads and is the most prominent sign of the interference effects . for @xmath183
this term disappears and the resulting sidebands have a form of symmetric dips due to the negative value of the coefficient @xmath184 in eq .
27 . in this case
the nonresonant tunneling channels modify the dips center position and its fwhm in comparison with the position and fwhm of sidebands presented for @xmath177 .
note , that for a qd coupled with two leads , the corresponding sidebands ( the thin dotted line ) are not fully symmetric curves as in this case the last term of eq .
28 disappears for @xmath185 and not for @xmath186 ( c.f .
@xcite ) . for @xmath181
the corresponding sidebands are fully asymmetric curves as in this case the second term in eq .
27 ( which introduces asymmetry ) disappears .
again , for a qd coupled with two leads the corresponding sidebands ( the thin broken curve ) are described by not fully asymmetric curves as the second term of eq . 28 , which gives a symmetric contribution to sidebands , disappears for @xmath187 and not for @xmath188 .    , @xmath189 and @xmath190 ( solid , broken and dotted lines , respectively ) against @xmath14 in the system of a qd coupled with three leads .
the upper ( lower ) panel corresponds to @xmath177 ( @xmath191 ) . @xmath172 and other parameters as in fig .
the values of @xmath192 for @xmath177 ( upper panel ) have been multiplied by a factor of 20 for illustrating purposes . ]    in the next fig .
[ fig4 ] we show the currents flowing in the qd - three leads ( qd-3l ) system for @xmath177 and @xmath191 . for vanishing values of the over - dot coupling the current @xmath147
is characterized by a sequence of the symmetric peaks , but the current @xmath190 is a superposition of the asymmetric structures placed in the points where the symmetric sidebands occur on @xmath147 curve . analyzing the current @xmath190 according to eq .
29 we have for the parameters in fig . [ fig4 ] :    @xmath193,\end{aligned}\ ] ]    where @xmath194 .
one can see that each sideband is the sum of the peak ( two first terms in eq .
29 ) and the dip ( the last two terms in eq .
29 ) resulting in the asymmetric structure shown in fig .
[ fig4 ] .     and @xmath15 . the upper ( middle ) panel shows @xmath147 ( @xmath190 ) for @xmath177 and the lower panel gives @xmath147 for @xmath191 .
@xmath172 , @xmath173 , @xmath195 . ]    to learn more about the influence of the third electrode and additional over - dot coupling between leads we present in fig .
[ fig5 ] the currents @xmath147 for @xmath177 and @xmath191 and @xmath190 for @xmath177 . for vanishing coupling between leads ( @xmath177 )
the current @xmath147 exhibits a well known sideband structure for @xmath196 and for small frequencies , @xmath197 , the two broad maxima at @xmath198 are present . at the same time , the current @xmath190 exhibits the asymmetric structures centered on the @xmath199 plane at the points where photon sidebands occur on the @xmath147 curves .
these asymmetric structures on the @xmath190 curve exist also at @xmath197 . on the other hand , similar structure of the @xmath164 ( the lower panel of fig .
[ fig5 ] ) is obtained for @xmath191 , i.e. we observe a number of asymmetric resonances separated by the photon energy for @xmath200 .
notice the similarity of both pictures , i.e. @xmath147 calculated for @xmath191 and @xmath190 for @xmath177 , respectively .
note , however , the different scale for these currents .     against @xmath14 for different values of the over - dot coupling strength between the left and right leads .
the broken , dotted , thin solid and thick solid lines correspond to @xmath201 and @xmath171 , respectively . @xmath202 and @xmath203 , @xmath204 , @xmath205 , @xmath206 , @xmath159 . ]    in fig .
[ fig6 ] we present @xmath147 for different over - dot coupling assuming a strong asymmetry of the applied microwave field ( ac potential is applied only to the @xmath0-lead in the qd-3l system ) .
the additional @xmath1-lead is characterized by the chemical potential @xmath207 . for better demonstration of the influence of the ac potential on @xmath147 we moved down each curve by the constant value @xmath202 ,
for @xmath177 we observe a shoulder on the left side of the main peak and a small negative current for the positive values of @xmath14 .
this picture is similar to the known results ( for @xmath177 ) obtained experimentally and theoretically in the qd-2l systems , e.g. refs .
@xcite . with the increasing over - dot coupling @xmath110
the height of the main resonant peak decreases and disappears at all for @xmath191 but the enhancement of the current for negative values of @xmath14 ( for @xmath208 ) is more and more expressive . at the same time , for all values of @xmath157 we observe a negative current for the small positive @xmath14 with a greater absolute value for stronger coupling between leads . for greater values of @xmath157 the shape of the @xmath147 curve
is changed dramatically and for @xmath186 , @xmath147 becomes nearly reversed in comparison with that calculated for @xmath177 .     against @xmath14 for the case of the microwave field
applied only to the qd with @xmath209 ( thin line ) or @xmath158 ( thick line ) - upper panel .
the lower panel corresponds to the case of the microwave field applied to the qd and @xmath1 lead with @xmath210 and @xmath211 or @xmath212 ( thin or thick lines , respectively ) .
the broken lines show the results when the microwave field applied to the qd and @xmath1 lead are out of phase ( with the phase difference of @xmath213 ) , @xmath205 , @xmath206 , @xmath159 . ]    in fig .
[ fig7 ] we analyze the influence of the third lead ( here named as @xmath1 ) on the current @xmath147 when the external microwave field is applied to this lead and to the qd with @xmath210 . for comparison ,
we add in the upper panel the results for @xmath147 obtained for the case when this additional third lead is not irradiated by the microwave field . in this case , as before , see e.g. fig .
[ fig2 ] , we observe typical sidebands on the current curves ( the difference between the lead chemical potentials is small in comparison with the amplitude @xmath18 ) .
however , after including the third lead irradiated by the external microwave field the dependence of the current @xmath147 on the gate voltage ( or equivalently on the qd energy level position ) is quite different - compare the thin or thick solid lines of the upper and lower panels . for smaller values of @xmath214 and @xmath129 ,
the averaged current @xmath215 is very similar to the corresponding current @xmath216 obtained by applying the external microwave field only to one lead ( see the broken line in fig .
these curves are , however , related between themselves by a relation @xmath217 .
now we can observe a small negative current at small negative values of @xmath14 and some enhancement of the current on the right side of the main peak .
similarly , the significant differences between the corresponding currents are observed also for greater values of the amplitudes @xmath214 and @xmath129 ( compare the thick solid lines in the upper and lower panels of fig .
[ fig7 ] ) .
note , that very similar behavior of the current @xmath147 as the function of the gate voltage is observed if we compare the case when the microwave field is applied only to one lead and the case when the microwave field is applied simultaneously to the qd and @xmath218lead but with the phase difference of @xmath213 - compare the thin broken curve in fig .
[ fig7 ] with the broken curve in fig .
[ fig6 ] .     against @xmath14 in the system of a qd coupled with three ( thick lines ) or two ( thin lines ) leads .
the upper ( lower ) panel corresponds to @xmath177 ( @xmath191 ) . in a qd - three leads system @xmath219 , @xmath220 and @xmath221 and @xmath222 ( solid , broken and dotted curves , respectively ) . in a qd
coupled with two leads @xmath219 , @xmath223 and @xmath222 ( solid , broken and dotted curves , respectively ) .
the values of solid curves have been multiplied by a factor of 10 for the illustrating purposes .
@xmath159 , @xmath172 , @xmath173 . ]
next , we consider the qd-2l system for different values of the chemical potential of the right lead ( @xmath221 and 10 ) and show the results of including the next lead ( @xmath224lead ) to the system with @xmath220 . in fig .
[ fig8 ] we present the results obtained for @xmath147 in the case when the microwave field is applied only to the qd . for the case @xmath174 and @xmath177 the sidebands are very clearly visible and the corresponding peaks are lower and broader for the qd-3l system as we discussed before ( the upper panel ) . for the non - zero over - dot coupling between the left and right leads the sideband peaks get asymmetric forms and for @xmath191 they become fully asymmetric ( the lower panel ) .
the current for the qd-2l system is also composed of a number of asymmetric components although now these forms are fully asymmetric for @xmath225 as we know from the earlier discussion . for greater values of @xmath226 ,
the corresponding current flowing in the qd-2l system achieves ( for vanishing , as well as for nonzero coupling between leads ) greater negative values and its dependence on the qd energy level position is well - marked in comparison to the results characterizing the qd-2l system .     against @xmath14 for the qd-2l systems : @xmath227 - thin solid lines and @xmath219 , @xmath228 - thick solid lines and the qd-3l system : @xmath229 , @xmath230 - broken lines .
the upper ( lower ) part corresponds to @xmath177 ( @xmath191 ) and @xmath159 , @xmath231 , @xmath173 .
the dotted lines correspond to the sum of the currents flowing in two qd-2l systems . ]    in order to emphasize the influence of the additional lead on the currents we show in fig .
[ fig9 ] the current @xmath164 for qd-2l and qd-3l systems calculated for the parameters for which the corresponding curves are relatively simple .
we assumed the small amplitude of the qd energy level oscillations , @xmath231 , and @xmath159 , for which ( for @xmath177 ) only the central peak corresponding to elastic tunneling is visible on the @xmath147 curves .
we show @xmath164 obtained for two different qd-2l systems which can be viewed as components of more complicated qd-3l system .
we observe that due to the interference effects , the current @xmath147 flowing in the qd-3l system is not simply a sum of currents ( dotted lines ) flowing in the corresponding qd-2l systems .
the difference between this sum and the current corresponding to the qd-3l system is relatively large and exists independently of the coupling between leads .     against @xmath14 and @xmath15 for @xmath177 .
we compare the results obtained for the qd coupled with two leads - the upper ( middle ) panel - for @xmath232 , @xmath233 , @xmath234 , @xmath235 , @xmath236 ( @xmath232 , @xmath237 , @xmath234 , @xmath235 , @xmath238 ) with the results obtained for the qd coupled with three leads - the lower panel - @xmath234 , @xmath239 , @xmath240 , @xmath241 , @xmath232 , @xmath242 , @xmath243 , @xmath159 . ]    in order to present more information about the differences between the electron transport in the qd-3l and qd-2l systems we display in fig .
[ fig10 ] the derivative @xmath244 as a function of @xmath14 and @xmath15 .
the lowest panel corresponds to the qd-3l system and two other panels correspond to the qd-2l systems .
these qd-2l systems are characterized by such parameters that combined together give us the considered qd-3l system .
one can see , that the considered characteristics of the electron transport in the qd-3l system are not simply the algebraic sum of the corresponding curves of both qd-2l systems .
in all three cases shown in fig [ fig10 ] , the position of the corresponding minima and maxima ( along the @xmath245axis ) can be identified with the values of the leads chemical potentials .
however , the corresponding structures are less clear in the case of the qd-3l system in comparison with those for the qd-2l models .    , @xmath0 and @xmath1 .
the @xmath92-lead is coupled with the qd only - the left panels and with the qd and two other leads , @xmath246 - the right panels .
the couplings between the qd and @xmath247 leads are changed periodically .
the upper , middle and lower panels correspond to @xmath248 and -3 , respectively .
@xmath249 . the thin , thick and broken curves correspond to @xmath250 and @xmath251 currents , respectively . ]    in the next step of our investigations of the electron transport in the qd-3l systems we consider the time - dependent currents flowing in response to the time - dependent barriers between the qd and leads or in response to suddenly removed ( or inserted ) connection of the qd with one of the leads . as a first we assume a rectangular - pulse modulation applied to the qd-@xmath0 lead and qd-@xmath1 lead barriers .
we assume that these modulations are with a phase difference of @xmath213 . in the first ( second ) half - cycle , @xmath252 ( @xmath253 ) and the qd
is coupled only to the @xmath1 lead ( @xmath0 lead ) .
in addition , the qd is coupled to the next , say @xmath254lead , with a constant value @xmath255 . in the following
we consider the time - dependent currents @xmath256 , @xmath257 and @xmath258 for the three specific conditions : @xmath259 , @xmath260 and @xmath261 .
in addition , we assume @xmath262 , @xmath263 and take for the period of the considered barrier modulation @xmath264 . in a such case
we integrate numerically the corresponding set of the differential equations for the matrix elements of the evolution operator and in the next step calculate the currents according to the formula @xmath265 , where @xmath266 or @xmath1 and @xmath267 is given in eq .
5 ( or similar to it ) .
we checked that the qd charge hardly depends on the additional @xmath110 couplings .
although the qd charge is almost insensitive to the additional over - dot couplings the currents demonstrate such dependence .
especially visible are the differences for the case when the chemical potential @xmath268 of the third electrode @xmath92 lies between the chemical potentials of two other leads , see fig .
[ fig11 ] b , e . for other values of @xmath268 ,
the influence of the over - dot tunneling channels for the parameters considered here is smaller .
note , that after abrupt changing of the coupling between the qd and @xmath0 , or @xmath1 leads the currents @xmath269 and @xmath251 rapidly change too , and after a short time reach the steady values .
the qd coupled with three leads could be considered as the three - state system .
we observe that for some values of the lead chemical potentials the currents change their values periodically e.g. from zero to the positive value ( see @xmath251 in fig .
[ fig11]b ) , from zero to the negative value ( see @xmath270 in fig . [
fig11]b ) or from the negative to the positive value ( see @xmath215 in fig .
[ fig11]b).the additional couplings between leads modify the values of the currents but the qualitative picture remains the same .
note , that in the first moment after abrupt changing of the coupling between the qd and @xmath0 or @xmath1 lead we have @xmath271 as in this case @xmath272 ( not shown here ) .
after some delay time the qd charge stabilizes around its equilibrium value , the currents tend to constant values and their sum is equal to zero .    , @xmath0 and @xmath1 leads in the system shown in fig .
the thin solid , thick solid and broken curves correspond to different couplings of the @xmath92 lead to the other elements of the system : ( @xmath273 , @xmath274 ) , ( @xmath273 , @xmath275 , @xmath276 ) and ( @xmath277 ) , respectively . ]    in the last step we consider the response of the currents to the abrupt inclusion into the qd-2l system of the third electrode ( in our case , @xmath92-lead ) .
the results are presented in fig .
[ fig12 ] together with the schematic view of the subsequent tunnel connections between the qd and the three leads .
we show the time - dependent currents @xmath37 , @xmath278 and @xmath258 corresponding to the three different ways of inclusion of the @xmath92-lead .
we assumed the chemical potentials @xmath279 , @xmath260 and @xmath280 . consider the current @xmath251 flowing from the @xmath1-lead characterized by the highest chemical potential @xmath226 . before adding to the system of the @xmath92-lead the current @xmath251 has a constant value and flows from @xmath1-lead through the qd energy level to @xmath0-lead . when the @xmath92-lead is included into the system ( the tunneling coupling @xmath281 changes abruptly at @xmath282 from zero to nonzero value ) the current @xmath251 is almost unchanged - its value decreases slightly without any transients at short times after the time @xmath282 .
next , we consider the case when the @xmath92-lead is abruptly connected simultaneously with the qd and @xmath0-lead .
now the current @xmath251 decreases significantly during the short time after the moment of inclusion of @xmath92-lead and settles to its constant value .
note , that @xmath251 decreases despite the additional charge transfer channel between @xmath1 and @xmath92 leads ( through the qd energy level ) .
the destructive interference appears in this case as we have two transmission channels for tunneling electrons between @xmath1 and @xmath0 leads .
however , the constructive interference is visible if we consider the next case when the @xmath92-lead is abruptly coupled with the qd-2l system assuming nonzero values of @xmath281 , @xmath283 and @xmath284 .
now , we have one additional charge transfer channel ( from @xmath1-lead ) in comparison with the former case .
the current @xmath251 rapidly increases with some fluctuations and after the time @xmath285 decreases to the constant value greater than in qd-2l system .
similar analysis can be made considering the currents @xmath286 and @xmath270 although the transient current changes are more visible now at short times after the abrupt inclusion of additional electron tunneling channels .
the above discussion concerns the specific values of the lead chemical potentials and the position of the qd energy level .
nevertheless , the similar qualitative conclusions can be made also for other values characterizing the considered system .
we have studied the time - dependent tunneling transport through the qd coupled with three metal leads using the evolution operator technique .
the time - dependent qd charge and currents were determined in terms of the appropriate evolution operator matrix elements . applying the wide band limit to the integrodifferential equations satisfied by the evolution operator matrix elements we were able to give the analytical expressions for the time - averaged currents and differential conductance .
we considered the external harmonic microwave fields applied to different parts of the considered system , as well as , the rectangular - pulse modulation imposed on different qd - leads barriers .
in addition , we have studied also the time dependence of the currents due to abrupt inclusion into the qd - two leads system of the third electrode .
we have considered also the effect of the additional couplings between leads ( we coupled one of the leads with the other two leads ) on the conductance and current flowing in the system .
our main results can be summarized as follows : + - for the vanishing nonresonant tunneling path , @xmath106 , and for the parameters for which the photon - assisted sidebands are clearly visible on the @xmath147 curve the subsequent sideband peaks have the lorentz - type form with the fwhm equal @xmath167 in comparison with @xmath169 for the qd - two leads system . for the increasing value of @xmath110 the form of the sidebands transforms from the lorentz - type to the fully asymmetric form for @xmath287 . for greater @xmath157 the form of the sidebands changes and gains again the lorentz - type shape for @xmath183 . for the qd
coupled with two leads the corresponding values of @xmath157 are equal to @xmath288 and @xmath289 , respectively .
+ - for the vanishing @xmath110 the differential conductance curve , e.g. @xmath290 , possesses the sidebands of the lorentz - type localized at @xmath291 . for @xmath292 these sidebands
are described by the superposition of two parts , the lorentz - type and the asymmetric one centered at @xmath293 weighed by the factors @xmath294 and @xmath295 , respectively . for the qd - two leads system
the corresponding @xmath116th sidebands are centered at @xmath296 and their symmetric and asymmetric parts are weighted by the factors @xmath297 and @xmath298 , respectively .
+ - the symmetry properties of the sidebands corresponding to the current flowing from the given lead depend on the position of the chemical potential of this lead in comparison with the chemical potentials of the other two leads . taking , for example @xmath106 and @xmath226 localized in the middle between @xmath268 and @xmath299 one can observe on the @xmath164 curve the sidebands of nearly regular ( lorentz - type ) forms while the sidebands on the @xmath190 curve have asymmetric structures .
however , in the presence of the nonresonant tunneling path the sidebands on the @xmath164 curve change their form and for @xmath181 they have a fully asymmetric shape . on the other hand ,
the sidebands on the @xmath190 curve have for @xmath181 a nearly lorentz - like shape .
+ - especially large interference effects can be observed if we compare the current flowing in the qd-3l system with the sum of currents flowing in the two qd-2l systems which can be viewed as the components of the considered more complicated qd-3l system .
the difference between them is relatively large independently of the over - dot coupling between leads ( see fig .
[ fig9 ] ) + - in the case of strong asymmetry of the applied external field ( @xmath300 ) we observe for @xmath177 a shoulder on the left side of the main resonant peak on the @xmath147 curve vs. the gate voltage . with the increasing over - dot coupling between @xmath92 and
@xmath301leads the main resonant peak disappears and transforms in a dip for strong coupling @xmath110 . at the same time
a shoulder of this curve increases with the increasing @xmath110 .
+ - let us consider the time dependence of currents flowing in response to the time - dependent barriers between the qd and two leads ( we assume a constant coupling of the qd with the third lead ) .
for the assumed rectangular - pulse modulation applied to the qd-@xmath0 lead and qd-@xmath1 lead barriers one can consider the qd-3l system as a three - state one .
for example , the currents @xmath302 , @xmath278 and @xmath258 change their values periodically between zero and positive , positive and zero and negative and positive values , respectively ( see fig .
the additional couplings between the leads introduce only small quantitative changes .
+ * acknowledgements : * the work of one of us ( rt ) has been partially supported by the kbn grant no .
pbz - min-008/p03/2003 .      in this section
we present the derivations of the general equations satisfied by the required functions @xmath29 , @xmath303 , @xmath304 and @xmath305 , @xmath306 , @xmath307 needed for the calculations of the qd charge @xmath54 and the currents flowing in the considered system . in the next step , using the wbl approximation we simplify these equations and give the analytical solutions for them . let us begin from the derivation of the integro - differential equation satisfied by @xmath29 . writing down the formal solution of eq . ( 11 ) @xmath308 and inserting them to the formal solution of eq .
( 10 ) @xmath309 one can obtain after straightforward calculations the function @xmath81 expressed in the terms of @xmath29 , only @xmath310 where @xmath311 we remember , that the wave vectors @xmath312 and @xmath7 correspond to the left lead and the first and second right leads , respectively . inserting into eq .
9 the expressions for the functions @xmath313 , @xmath314 and @xmath315 , eqs .
( 17 , 18 ) , we can write the integro - differential equation for @xmath29 in the form @xmath316 u_{dd}(t_1,t_0 ) \nonumber\\ & + & \sum^\infty_{j=1}(-i)^{2j-1 } \int\limits^t_{t_0}\int\limits^{t_1}_{t_0 }   \ldots \int\limits^{t_{2j-1}}_{t_0 } dt_1 \ldots dt_{2j } \sum_{\vec k_1,\vec k_2 , \ldots,\vec k_j } [ r_{d\vec k_1}(t , t_1 ) r_{\vec k_1\vec k_2 } ( t_2,t_3 )   \ldots\nonumber\\ & & r_{\vec k_{j-1}\vec k_j}(t_{2j-2},t_{2j-1 } ) \tilde v_{\vec k_jd}(t_{2j } ) \nonumber\\ & + &   \tilde v_{d\vec k_1}(t ) r_{\vec k_1\vec k_2}(t_1,t_2)\ldots r_{\vec k_jd } ( t_{2j-1},t_{2j } ) ] u_{dd}(t_{2j},t_0 ) \nonumber\\ & + & \sum^\infty_{j=1}(-i)^{2j } \int\limits^t_{t_0}\int\limits^{t_1}_{t_0 }   \ldots \int\limits^{t_{2j}}_{t_0 } dt_1\ldots dt_{2j+1 } \sum_{\vec k_1\vec k_2   \ldots \vec k_j } [ r_{d\vec k_1}(t , t_1 ) r_{\vec k_1\vec k_2 } ( t_2,t_3 ) \ldots \nonumber\\ & & r_{\vec k_{j-1}\vec k_j } ( t_{2j-2},t_{2j-1 } ) \ , r_{\vec k_jd}(t_{2j},t_{2j+1 } ) \\ & + & \sum_{\vec k_{j + 1 } } \tilde v_{d\vec k_1}(t ) \ , r_{\vec k_1\vec k_2}(t_1,t_2 ) \ldots r_{\vec k_j,\vec k_{j+1}}(t_{2j-1},t_{2j } ) \tilde v_{\vec k_{j+1}d } ( t_{2j+1 } ) ] u_{dd}(t_{2j+1},t_0 ) \ , .
\nonumber % 10 ( a.5)\end{aligned}\ ] ] this rather untractable general equation can be greatly simplified using the wbl approximation .
assuming @xmath317 , @xmath318 , the multidimensional time integrations and summations over the wave vectors can be performed giving in result @xmath319 with the solution @xmath320 here @xmath321 and @xmath60 , @xmath61 being the bandwidth of the lead energy band ( @xmath62 ) . the function @xmath81 given in eq .
( 19 ) can be reduced within the wbl approximation to the form @xmath322 where @xmath323 in order to calculate @xmath71 we write down accordingly with eq .
( 3 ) the corresponding set of coupled differential equations for the functions @xmath324 , @xmath82 , @xmath325 and @xmath326 .
the subsequent steps of the calculations are similar to those performed in the derivation of eq .
( 34 ) . inserting the formal solutions for @xmath82 , @xmath327 @xmath328 \ , , \nonumber\\ \\ % 13     ( a.11 ) & & u_{\vec p\vec k}(t , t_0 ) =   - i\int\limits^t_{t_0 }   dt_1 \left[\tilde v_{\vec pd}(t_1 ) u_{d\vec k } ( t_1,t_0 ) +    \sum_{\vec k_1 } \tilde v_{\vec pk_1}(t_1 ) u_{\vec k_1\vec k } ( t_1,t_0)\right ] \
, , % 14   ( a.12)\end{aligned}\ ] ] into the differential equation satisfied by @xmath71 one obtains the derivative @xmath329 expressed in terms of @xmath71 and @xmath82 . on the other hand , on the basis of eqs .
( 33 , 34 ) the function @xmath82 can be represented in the form containing only @xmath330 .
finally , one obtains @xmath331 \nonumber\\ & + & ( -i)^3 \int\limits^t_{t_0 } \int\limits^{t_1}_{t_0 }
dt_1 dt_2 \left[\sum_{\vec k_1,\vec q_1 }   \tilde v_{d\vec k_1}(t ) \tilde v_{\vec k_1\vec q_1}(t_1 ) \tilde v_{\vec q_1\vec k }   ( t_2 ) + \sum_{\vec k_1,\vec r_1 } \ldots   \right ] \nonumber \\ & + & ( -i)^4 \int\limits^t_{t_0 } \int\limits^{t_1}_{t_0 } \int\limits^{t_2}_{t_0 }   dt_1 dt_2 dt_3 \left[\sum_{\vec q_1,\vec k_1,\vec q_2 }   \tilde v_{d\vec q_1}(t ) \tilde v_{\vec q_1\vec k_1}(t_1 ) \tilde v_{\vec k_1\vec q_2 }   ( t_2 ) \tilde v_{\vec q_2\vec k}(t_3 ) \right .
+ \nonumber\\ & & \hspace{3cm}\left .
\sum_{\vec r_1,\vec k_1\vec q_2 } \ldots + \sum_{\vec q_1\vec k_1 \vec r_2 } \ldots + \sum_{\vec r_1\vec k_1\vec r_2 } \ldots \right ] + \ldots \nonumber\\ & + & ( -i)^2 \int\limits^t_{t_0 } dt_1 \left[\sum_{\vec k_1 } \tilde v_{\vec d\vec k_1}(t )   \tilde v_{\vec k_1d } ( t_1 )   + \sum_{\vec q_1 } \tilde v_{d\vec
q_1}(t )    \tilde v_{\vec q_1d}(t_1)\right .   \nonumber\\ & & \left .
\hspace{3 cm } + \sum_{\vec r_1}\tilde v_{d\vec r_1}(t ) \tilde v_{\vec r_1d}(t_1)\right ]    u_{d\vec k}(t_1,t_0 )   \nonumber \\ & + & ( -i)^3 \int\limits^t_{t_0 } \int\limits^{t_1}_{t_0 } dt_1 dt_2    \left[\sum_{\vec k_1,\vec q_1 } \tilde v_{d\vec k_1}(t ) \tilde v_{\vec k_1\vec q_1}(t_1 ) \tilde v_{\vec q_1d}(t_2 ) \right.\nonumber\\ & & \hspace{3cm}\left .
+ \sum_{\vec k_1,\vec r_1 } \ldots + \sum_{\vec q_1\vec k_1 }     \ldots + \sum_{\vec r_1\vec k_1 } \ldots \right ] u_{d\vec k}(t_2,t_0 ) + \ldots % 15   ( a.13)\end{aligned}\ ] ] in the wbl approximation this equation reduces to the form @xmath332 with the solution @xmath333 where @xmath334 . taking into account eqs .
( 33 , 34 ) the function @xmath82 can be written in terms of @xmath335 @xmath336    in order to calculate @xmath337 we first write down the coupled set of equations ( on the basis of eq .
( 3 ) ) satisfied by the functions @xmath72 , @xmath338 , @xmath339 and @xmath83 . inserting the formal solutions for the functions @xmath338 , @xmath339 and @xmath83 into the differential equation for the function @xmath72 one obtains the derivative @xmath340 expressed in terms of @xmath341 , @xmath83 , @xmath338 and @xmath339 .
inserting again into this equation the formal solutions for the functions @xmath338 , @xmath339 and @xmath83 , and repeating this process again and again , one obtains : @xmath342\nonumber\\ & + & ( -i)^4 \int\limits^t_{t_0 } \int\limits^{t_1}_{t_0 } \int\limits^{t_2}_{t_0 }   dt_1 dt_2 dt_3 \left[\sum_{\vec k_1\vec q_1\vec k_2 } \tilde v_{d\vec k_1}(t )   \tilde v_{\vec k_1\vec q_1}(t_1 ) \tilde v_{\vec q_1\vec k_2}(t_2 )   \tilde v_{\vec k_2\vec q } ( t_3 )   \right .
\nonumber \\ & & \left .
\hspace{3 cm } + \sum_{\vec k_1\vec r_1\vec k_2}\ldots   + \ldots \right ] +   \ldots   \nonumber\\ % % % % % % % % % % % % % & + & ( -i)^2 \int\limits^t_{t_0 } dt_1 \sum_{\vec p=\vec k_1,\vec q_1,\vec r_1 }    \tilde v_{d\vec p}(t ) \tilde v_{\vec p d}(t_1 )     u_{d\vec q } ( t_1,t_0 )   \nonumber \\ & + & ( -i)^3 \int\limits^t_{t_0 } \int\limits^{t_1}_{t_0 } dt_1 dt_2    \left[\sum_{\vec q_1,\vec k_1 } \tilde v_{d\vec q_1}(t )   \tilde v_{\vec q_1\vec k_1}(t_1 )    \tilde v_{\vec k_1d}(t_2 ) \right.\nonumber\\ & + & \left .
\sum_{\vec r_1\vec k_1 } \ldots + \sum_{\vec k_1\vec q_1 }     \ldots + \sum_{\vec k_1\vec r_1 } \ldots \right ] u_{d\vec q}(t_2,t_0 ) + \ldots % 18   ( a.19)\end{aligned}\ ] ] in the wbl approximation this equation is reduced to the following form : @xmath343 where @xmath344 and has the solution as follows : @xmath345 the function @xmath73 is identical with @xmath72 .
to calculate , e.g. the current @xmath37 , we still need the functions @xmath346 and @xmath347 .
the function @xmath83 can be obtained solving the set of the coupled differential equations for the functions @xmath72 , @xmath348 , @xmath339 and @xmath83 . writing down the formal solution for @xmath83 and inserting into it , in the first step , the formal solutions for @xmath349 and @xmath339 and , in the second step , the formal solutions for @xmath72 and @xmath83 , and so on , one obtains @xmath350
u_{d\vec q}(t_1,t_0 ) \nonumber\\ & + & ( -i)^3 \int\limits^t_{t_0 } \int\limits^{t_1}_{t_0}\int\limits^{t_2}_{t_0 } dt_1 dt_2    dt_3 \left[\sum_{\vec q_1,\vec k_1 }   \tilde v_{\vec kq_1}(t_1 ) \tilde v_{\vec q_1\vec k_1}(t_2 ) \tilde v_{\vec
k_1d}(t_3 ) + \sum_{\vec r_1\vec k_1 } \ldots \right ] u_{d\vec q } ( t_3,t_0 ) \nonumber \\ & + & \ldots + ( -i ) \int\limits^t_{t_0 } dt_2 \tilde v_{\vec k\vec q } ( t_1,t_0 )   + ( -i)^3 \int\limits^t_{t_0 } \int\limits^{t_1}_{t_0 } \int\limits^{t_2}_{t_0 }    dt_1 dt_2 d_3\nonumber \\ & & \hspace{3 cm } \left[\sum_{\vec q_1,\vec k_1}\tilde v_{\vec kq_1}(t_1 ) \tilde v_{\vec q_1\vec k_1}(t_2 )   \tilde v_{\vec k_1d}(t_3 ) + \sum_{\vec r_1\vec k_1 } \ldots\right ] + \ldots % 26   ( a.24)\end{aligned}\ ] ] under the wbl approximation this equation reduces to the form : @xmath351 where @xmath341 is given in eq . (
the function @xmath84 has the identical form as @xmath352 .
m. tsukada and n. shima , in : `` dynamical processes and ordering on solid surfaces '' , eds .
a. yoshimori and m. tsukada ( springer , berlin 1995 ) , p. 34 ; s. tsuneyuki , n. shima , m. tsukada , surf .
* 186 * 26 ( 1987 )