the experiments on resistivity of high mobility two - dimensional electron gas ( 2deg ) in presence of a relatively weak magnetic field and microwave radiation led to a discovery of striking zero - resistance states ( zrs ) induced by a microwave field by mani _
_ @xcite and zudov _ et al .
other experimental groups also found the microwave induced zrs in various 2deg samples ( see e.g. @xcite ) .
a similar behavior of resistivity is also observed for electrons on a surface of liquid helium in presence of magnetic and microwave fields @xcite .
these experimental results obtained with different systems stress the generic nature of zrs .
various theoretical explications for this striking phenomenon have been proposed during the decade after the first experiments @xcite .
an overview of experimental and theoretical results is give in the recent review @xcite .
in our opinion the most intriguing feature of zrs is their almost periodic structure as a function of the ratio @xmath0 between the microwave frequency @xmath1 and cyclotron frequency @xmath2 .
( in the following we are using units with electron charge @xmath3 and mass @xmath4 equal to unity ) . indeed , a hamiltonian of electron in a magnetic field is equivalent to an oscillator , it has a magneto - plasmon resonance at @xmath5 but in a linear oscillator there are no matrix elements at @xmath6 and hence a relatively weak microwave field is not expected to affect electron dynamics and resistivity properties of transport . of course
, one can argue that impurities can generate harmonics being resonant at high @xmath7 but zrs is observed only in high mobility samples and thus the density of impurities is expected to be rather low .
it is also important to note that zrs appears at high landau levels @xmath8 so that a semiclassical analysis of the phenomenon seems to be rather relevant .    in this work
we develop the theoretical approach proposed in @xcite .
this approach argues that impurities produce only smooth potential variations inside a bulk of a sample so that zrs at high @xmath9 appear from the orbits moving along sharp sample boundaries .
it is shown @xcite that collisions with boundaries naturally generate high harmonics and that a moderate microwave field gives stabilization of edge channel transport of electrons in a vicinity of @xmath10 producing at these @xmath9 a resistance going to zero with increasing microwave power .
this theory is based on classical dynamics of electrons along a sharp edge .
the treatment of relaxation processes is modeled in a phenomenological way by a dissipative term in the newton equations .
additional noise term in the dynamical equations takes into account thermal fluctuations .
the dissipation leads to synchronization of cyclotron phase with a phase of microwave field producing stabilization of edge transport along the edges in a vicinity of resonant @xmath11 values .
thus , according to the edge stabilization theory @xcite the zrs phase is related to a universal synchronization phenomenon which is a well established concept in nonlinear sciences @xcite .    while the description of edge transport stabilization @xcite captures a number of important features observed in zrs experiments it assumes that the contribution of bulk orbits in transport is negligibly small .
this assumption is justified for smooth potential variations inside the bulk of a sample .
however , a presence of isolated small scale scatterers inside the bulk combined with a smooth potential component can significantly affect the transport properties of electrons ( see e.g. @xcite ) .
also the majority of theoretical explanations of zrs phenomenon considers only a contribution of scattering in a bulk @xcite .
thus it is necessary to analyze how a scattering on a single impurity is affected by a combined action of magnetic and microwave fields . in this work
we perform such an analysis modeling impurity by a rigid circular disk of finite radius .
we show that the dynamics in a vicinity of disk has significant similarities with dynamics of orbits along a sharp edge leading to appearance of zrs type features in a resistivity dependence on @xmath9 .
the paper is composed as follows : in section ii we discuss the dynamics in edge vicinity , in section iii we analyze scattering on a single disk , in section iv we study scattering on many disks when their density is low , here we determine the resistivity dependence on @xmath9 and other system parameters , physical scales of zrs effects are analyzed in section v , effects of two microwave driving fields and other theory predictions are considered in section vi , discussion of the results is given in section vii .
we study various models which we list here for a reader convenience : wall model described by the newton equations ( [ eq1 ] ) , ( [ eq2 ] ) with microwave field polarization perpendicular to the wall ( model ( w1 ) equivalent to model ( 1 ) in @xcite ) ; the chirikov standard map description ( [ eq3 ] ) of the wall model dynamics called model ( w2 ) ( equivalent to model ( 2 ) in @xcite at parameter @xmath12 ) ; the single disk model with radial microwave field called model ( dr1 ) ; the chirikov standard map description ( [ eq3 ] ) of model ( dr1 ) called model ( dr2 ) ( here @xmath13 in ( [ eq3 ] ) , @xmath14 ) ; the model of a single disk in a linearly polarized microwave field and static electric field called model ( d1 ) ; the model of transport in a system with many disks called model ( d2 ) which extends the model ( d1 ) ; extension of model ( d2 ) with disk roughness and dissipation in space called model ( d3 ) ; the wall model ( w2 ) extended to two microwave fields is called model ( w3 ) .
we remind first the approach developed in @xcite . here , the classical electron dynamics is considered in a proximity of the fermi surface and in a vicinity of sample edge modeled as a specular wall .
the motion is described by newton equations @xmath15 where a dimensionless vector @xmath16 describes microwave driving field @xmath17 . here
an electron velocity @xmath18 is measured in units of fermi velocity @xmath19 and @xmath20 describes a relaxation processes to the fermi surface .
we also use the dimensionless amplitude of velocity oscillations induced by a microwave field @xmath21 . as in @xcite , in the following we use units with @xmath22 .
the last two terms @xmath23 and @xmath24 in ( [ eq1 ] ) account for elastic collisions with the wall and small angle scattering .
disorder scattering is modeled as random rotations of @xmath25 by small angles in the interval @xmath26 with poissonian distribution over time interval @xmath27 .
the amplitude of noise is assumed to be relatively small so that the mean free path @xmath28 is much larger than the cyclotron radius @xmath29 .
we note that the dissipative term is also known as a gaussian thermostat @xcite or as a landau - stuart dissipation @xcite . the dynamical evolution described by eq .
( [ eq1 ] ) is simulated numerically using the runge - kutta method .
following @xcite we call this system model ( w1 ) ( equivalent to model ( 1 ) in @xcite ) .
we note that for typical experimental zrs parameters we have : electron density @xmath30 , effective electron mass @xmath31 , microwave frequency @xmath32 , fermi energy @xmath33 , corresponding to @xmath34 , with fermi velocity @xmath35 . at such a frequency
the cyclotron resonance @xmath36 takes place at @xmath37 with the cyclotron radius @xmath38 . at such a magnetic field
we have the energy spacing between landau levels @xmath39 corresponding to a landau level @xmath40 . for a microwave field strength @xmath41
we have the parameter @xmath42 . with these physical values of system parameters
we can always recover the physical quantities from our dimensionless units with @xmath43 .
examples of orbits running along the edge of specular wall are given in @xcite ( see fig.1 there ) .
a microwave field creates resonances between the microwave frequency @xmath1 and a frequency of nonlinear oscillations of orbits colliding with the wall . due to a specular nature of this collisions
the electron motion has high harmonics of cyclotron frequency that leads to appearance of resonances around @xmath44 ( there is an additional shift of approximate value @xmath45 to @xmath11 values due to a finite width of nonlinear resonance ) .    to characterize the dynamical motion it is useful to construct the poincar section following the standard methods of nonlinear systems @xcite .
we consider the hamiltonian case at @xmath46 in absence of noise .
also we choose a linear polarized microwave field being perpendicular to the wall which is going along @xmath47-axis ( same geometry as in @xcite ) . in this case
the generalized momentum @xmath48 is an integral of motion since there are no potential forces acting on electron along the wall ( here we use the landau gauge with a vector potential @xmath49 ) .
the momentum @xmath50 determines a distance @xmath51 between a cyclotron center and the wall , which also remains constant in time .
the hamiltonian of the system has the form : @xmath52 where @xmath53 is the wall potential being zero or infinity for @xmath54 or @xmath55 .
thus , we have here a so called case of one and half degrees of freedom ( due to periodic time dependence of hamiltonian on time ) and the poincar section has continuous invariant curves in the integrable regions of phase space @xcite .    ) for @xmath56 ( left column ) and @xmath57 ( right column ) and different amplitudes of microwave field @xmath58 ( from bottom to top ) . here , the integral @xmath59 , trajectories start from wall with fixed @xmath60 . data for model ( w1 ) at @xmath46 , @xmath61 .
]    the poincar sections for ( [ eq1 ] ) , ( [ eq2 ] ) at @xmath62 and various amplitudes of microwave field @xmath63 are shown in fig .
it shows a velocity @xmath64 at moments of collision with the wall at @xmath65 as a function of microwave phase @xmath66 at these moments of time .
all orbits initially start at the wall edge @xmath65 with the initial velocity @xmath67 .
the value of @xmath68 is the integral of motion .
however , the kinetic energy of electron @xmath69 varies with time .
we see that at a small @xmath70 the main part of the phase space is covered by invariant curves corresponding to integrable dynamics .
however , a presence of chaotic component with scattered points is also visible in a vicinity of separatrix of resonances , especially at large @xmath71 .
the points at @xmath72 close to zero correspond to orbits only slightly touching the wall , while the orbits at @xmath73 have a large cyclotron radius and collide with the wall almost perpendicularly .
there are also sliding orbits which have the center of cyclotron orbit inside the wall @xmath74 but we do not discuss they here .
indeed , the orbits , which only slightly touch the wall ( @xmath75 ) , play the most important role for transport since the scattering angles in the bulk are small for high mobility samples and an exchange between bulk and edge goes via such type of dominant orbits @xcite .
we note that the section of fig .
[ fig1 ] at @xmath57 , @xmath76 is in a good agreement with those shown in fig.1b of @xcite .
however , here we have single invariant curves while in @xcite the curves have a certain finite width .
this happens due to the fact that in @xcite the poincar section was done with trajectories having different values of the integral @xmath68 that gave some broadening of invariant curves . for a fixed integral value
we have no overlap between invariant curves as it is well seen in fig .  [ fig1 ] here .    the phase space in fig .
[ fig1 ] has a characteristic resonance at a certain @xmath77 value which position depends on @xmath9 @xcite
. an approximate description of the electron dynamics and phase space structure can be obtained on a basis of the chirikov standard map @xcite,@xcite . in this description developed in @xcite an electron velocity
has an oscillating component @xmath78 ( assuming that @xmath79 ) and a collision with the wall gives a change of modulus of @xmath64 by @xmath80 ( like a collision with a moving wall ) . for small collision angles
the time between collisions is @xmath81 .
indeed , @xmath82 is the cyclotron period .
however , the time between collisions is slightly smaller by an amount @xmath83 : at @xmath84 an electron moves in an effective triangular well created by the lorentz force and like for a stone thrown against a gravitational field this gives the above reduction of @xmath85 ( formally this expression for @xmath85 is valid for sliding orbits but for orbits slightly touching the wall we have the same @xmath85 but with minus that gives the correction @xmath86 ) .
the same result can be obtained via semiclassical quantization of edge states developed in @xcite .
it also can be found from a geometric overlap between the wall and cyclotron circle .
this yields an approximate dynamics description in terms of the chirikov standard map @xcite : @xmath87 with the chaos parameter @xmath88 .
usually we are in the integrable regime with @xmath89 due to small values of @xmath63 used in experiments .
a developed chaos appears at @xmath90 @xcite . here
bars mark the new values of variables going from one collision to a next one , @xmath64 is the velocity component perpendicular to the wall , @xmath91 is the microwave phase at the moment of collision .
here we introduced a dimensionless parameter @xmath92 which is equal to @xmath12 for the case of the wall model @xmath93 considered here .
however , we will show that for the dynamics around disk with a radial field in model ( dr1 ) we have the same map ( [ eq3 ] ) with @xmath94 .
due to that it is convenient to write all formula with @xmath92 .
we note that a similar map ( [ eq3 ] ) describes also a particle dynamics in a one - dimensional triangular well and a monochromatic field @xcite .
the term @xmath95 in ( [ eq3 ] ) describes dissipation and noise .
the later gives fluctuations of velocity @xmath64 at each iteration ( @xmath96 ; corresponding to random rotation of velocity vector in ( [ eq1 ] ) ) .
damping from electron - phonon and electron - electron collisions contribute to @xmath97 .
the poincar sections of this map are in a good agreement with those obtained from the hamiltonian dynamics as it is seen in fig .
[ fig1 ] here and fig.1 in @xcite .
following @xcite we call this system model ( w2 ) ( equivalent to model ( 2 ) in @xcite ) .
a phase shift of @xmath98 by @xmath99 does not affect the dynamics and due to that the phase space structure changes periodically with integer values of @xmath9 .
indeed , the position of the main resonance corresponds to a change of phase by an integer number of @xmath99 values @xmath100 that gives the position of resonance at @xmath101 where @xmath4 is the nearest integer of @xmath102 and @xmath103 is the fractional part of @xmath9 . due to this relation
we have the different resonance position for @xmath56 and @xmath104 being in agreement with the data of fig .
[ fig1 ] at small values of @xmath63 when nonlinear corrections are small ( we have here @xmath12 ) .
thus at @xmath57 we have the resonance position at @xmath105 in agreement with fig .
[ fig1 ] ( right bottom panel ) . for @xmath106
we have @xmath107 and at @xmath56 the resonance position moves to negative value @xmath108 .
thus , at @xmath109 the resonance separatrix easily moves particles out from the edge at @xmath110 where they escape to the bulk due to noise .
in contrast at @xmath111 particles move along separatrix closer to the edge being then captured inside the resonance which gives synchronization of cyclotron phase with the microwave phase .
this mechanism stabilizes the transmission along the edge .     of electrons as a function of their dimensionless cyclotron center position @xmath112 between two walls and the frequency ratio @xmath0 .
the distance between specular walls is @xmath113 .
the amplitude of microwave field is @xmath114 with polarization parallel ( left panel ) and perpendicular ( right panel ) to the walls ( see text for more details ) . here
@xmath115 , @xmath116 , @xmath27 , the variation of @xmath0 is obtained by changing magnetic field ( @xmath117 ) keeping @xmath118 ; 100 electrons are simulated at each @xmath9 up to time @xmath119 .
density is proportional to color changing from zero ( black ) to maximal density ( white ) .
data for model ( w1 ) . ]    in @xcite it is shown that the orbits started in edge vicinity are strongly affected by a microwave field that leads to zrs type oscillations of transmission along the edge and longitudinal resistivity @xmath120 . the zrs structure appears both in the frame of dynamics described by ( [ eq1 ] ) ( model ( w1 ) ) and map description ( [ eq3 ] ) ( model ( w2 ) ) .
the physical mechanism is based on synchronization of a cyclotron phase with a phase of microwave driving that leads to stabilization of electron transport along the edge .
an extensive amount of numerical data has been presented in @xcite and we think there is no need to add more . here , we simply want to illustrate that even those orbits which start in the bulk are affected by this synchronization effect .
for that we take a band of two walls with a band width between them being @xmath121 . initially
100 trajectories are distributed randomly in a bulk part between walls when a cyclotron radius is not touching the walls @xmath122 .
their dynamics is followed during the run time @xmath119 according to eq .
( [ eq1 ] ) and a density distribution @xmath123 averaged in a time interval @xmath124 is obtained for a range of @xmath125 ( 261 values of @xmath9 are taken homogeneously in this interval ) .
the value of @xmath126 approximately corresponds to a distance propagation along the wall of @xmath127 at typical values @xmath128 , @xmath129 .
this is comparable with a usual sample size used in experiments @xcite .
similar values of @xmath130 were used in @xcite .     as a function of @xmath112 for microwave polarization parallel ( left panels ) and perpendicular ( right panels ) to the walls .
here we have no microwave at top panels , @xmath131 , @xmath132 at middle panels , @xmath131 , @xmath133 at bottom panels . in all panels we have noise amplitude @xmath116 as in fig .
[ fig2 ] , dashed curves in bottom panels are obtained with @xmath134 .
simulations are done with 500 trajectories , other parameters are the same as in fig .
data for model ( w1 ) . ]    the dependence of density @xmath135 on @xmath51 and @xmath9 are shown in fig .
[ fig2 ] for two polarizations of microwave field .
the data show that orbits from a bulk can be captured in edge vicinity for a long time giving an increase of density in a vicinity of edge .
this capture is significant around resonance values @xmath136 .
this is confirmed by a direct comparison of density profiles in fig .
[ fig3 ] at @xmath137 and @xmath138 . in the later case
we have a large density peak due to trajectories trapped in a resonance ( see fig .  [ fig1 ] ) where they are synchronized with a microwave field .
an increase of noise amplitude @xmath139 gives a significant reduction of the amplitude of these resonant peaks ( fig .
[ fig3 ] bottom panels ) .
the increase of density is more pronounced for polarization perpendicular to the wall in agreement with data shown in fig.2 of @xcite .
we also performed numerical simulations using eq .
( [ eq1 ] ) with a smooth wall modeled by a potential @xmath140 . for large values @xmath141
( e.g. @xmath142 ) we find the poincar sections to be rather similar to those shown in fig .
[ fig1 ] that gives a similar structure of electron density as in figs .
[ fig2],[fig3 ] .
a finite wall rigidity can produce a certain shift of optimal capture conditions appearing as a result of additional correction to a cyclotron period due to a part of orbit inside the wall .
the data presented in this section show that electrons from the bulk part of the sample can be captured for a long time in edge vicinity thereby increasing the electron density near the edge .
this effect is very similar to the accumulation of electrons on the edges of the electron cloud under zrs conditions that was reported for surface electrons on helium in @xcite .
however we have to emphasize that the confinement potential for surface electrons is very different from the hard wall potential assumed in our simulations , as a consequence our results can not be applied directly to this case .
it is possible that the formation of ballistic channels on the edge of the sample combined with the redistribution of the electrons density can effectively short the bulk contribution and induce directly a vanishing @xmath120 .
however , it is also important to understand how a scattering on impurities inside the bulk is affected by a microwave radiation .
we study this question in next sections .
it should be noted that resistivity properties of a regular lattice of disk antidots in 2deg had been studied experimentally @xcite and theoretically @xcite . but
effects of microwave field were not considered till present .    in our studies we model an impurity as a rigid disk of fixed radius @xmath143 keeping @xmath144 and changing @xmath117 . in a magnetic field
a cyclotron radius moves in a free space only due to a static _ dc-_electric field @xmath145 .
we fix the direction of @xmath145 along @xmath146axis and measure its strength by a dimensionless parameter @xmath147 . even in absence of a microwave field a motion in a vicinity of disk in crossed static electric and magnetic fields of moderate strength is not so simple .
the studies presented in @xcite and @xcite show that dynamics in disk vicinity is described by a symplectic disk map which is rather similar to the map ( [ eq3 ] ) .
it is characterized by a chaos parameter @xmath148 where @xmath149 is the drift velocity ; @xmath150 gives an amplitude of change of radial velocity at collision .
orbits from a vicinity of disk can escape for @xmath151 @xcite .     at @xmath152 ( from top to bottom ) ; right column panels : @xmath153 at @xmath76 . here
the integral of orbital momentum is @xmath154 , trajectories start from disk with fixed tangent velocity component @xmath155 .
data for model ( dr1 ) . ]
we start our analysis from the construction of the poincar section in presence of microwave field at zero static field .
to have a case with one and half degrees of freedom we start from a model case when a microwave field is directed only along radius from a disk center .
the dynamics is described by eq .
( [ eq1 ] ) with a dimensionless microwave amplitude @xmath63 .
the dynamical evolution is obtained numerically by the runge - kutta method . at first
we consider a case without dissipation and noise .
due to radial force direction the orbital momentum is an additional integral of motion ( as @xmath68 for the wall case ) and thus we have again 3/2 degrees of freedom .
we call this disk model with radial microwave field as model ( dr1 ) .
the poincar sections at the moments of collisions with disk are shown in fig .
[ fig4 ] ,  [ fig5 ] for model ( dr1 ) . here , @xmath156 is the radial component of electron velocity and @xmath98 is a microwave phase both taken at the moment of collision with disk .
we see that the phase space structure remains approximately the same when @xmath9 is increased by unity ( compare @xmath157 panels in fig .  [ fig4 ] ) .
this happens for orbits only slightly touching the disk ( small @xmath156 ) since the microwave phase change during a cyclotron period is shifted by an integer amount of @xmath99 ( in a first approximation at @xmath158 ) .
the similarity between the wall and disk cases is directly seen from fig .  [ fig5 ] as well as periodicity with @xmath159 .
in fact in the case of disk with a radial field the dynamics can be also described by the chirikov standard map ( [ eq3 ] ) where @xmath64 should be understood as a radial velocity @xmath156 at the moment of collision .
the second equation has the same form since the change of the phase between two collisions is given by the same equation but with the parameter @xmath94 .
this expression for @xmath92 is obtained from the geometry of slightly intersecting circles of radius @xmath160 for disk and radius @xmath161 for cyclotron orbit ( the angle segment of cyclotron circle is @xmath162 ) . for @xmath163
this expression naturally reproduces the wall case while at @xmath158 we have the correction term proportional to @xmath164 going to zero that also well corresponds to the geometry of two disks .
after such modification of @xmath92 we find that the resonance positions @xmath165 are proportional to @xmath92 .
thus the model ( dr1 ) reduced to the map description ( [ eq3 ] ) at @xmath166 is called model ( dr2 ) .    the expression for @xmath167 works
rather well . indeed , for @xmath168 in fig .
[ fig5 ] we obtain @xmath169 for model ( w2 ) and @xmath170 for model ( dr2 ) .
these values are in a good agreement with numerical values @xmath171 for model ( w2 ) and @xmath172 for model ( dr1 ) . in the later case
the agreement is less accurate due to a larger size of nonlinear resonance . the comparison of poincar sections given by the chirikov standard map ( [ eq3 ] ) and the dynamics from newton equations , shown in fig .
[ fig5 ] , confirms the validity of map description .    according to the well established results for the chirikov standard map @xcite we find for models ( w1 ) , ( w2 ) and ( dr1 ) , ( dr2 ) the width of separatrix @xmath173 and the corresponding resonance energy width @xmath174 : @xmath175 where @xmath176 is the resonance shift produced by a resonance half width @xmath177 .
this relation shows that for the disk case this energy is increased by a factor @xmath92 compared to the wall case . in majority of our numerical simulations
we have @xmath178 .
thus a radial field models ( dr1 ) , ( dr2 ) represent a useful approximation to understand the properties of dynamics in a disk vicinity but a real situation corresponds to a linear microwave polarization and the poincar section analysis should be modified to understand the dynamics in this case .     and @xmath168 ( top , middle ) and @xmath179 ( bottom ) .
top panels are obtained from the chirikov standard map ( [ eq3 ] ) at @xmath12 ( left top panel ) , corresponding to the wall model ( w2 ) , and at @xmath180 ( right top panel ) , corresponding to the disk model ( dr2 ) .
middle and bottom panels are obtained from solution of newton equations ( [ eq1 ] ) for wall ( left ) and disk ( right ) .
other parameters are as in fig .
[ fig1 ] and fig .
[ fig4 ] , @xmath64 and @xmath156 are expressed in units of @xmath19 .
there is no dissipation and no noise . ]
due to that we start to analyze the scattering problem on a disk in presence of weak static field @xmath181 and microwave field @xmath63 using eq .
( [ eq1 ] ) . for the scattering problem
we find more simple to have dissipation to work only at the time moments of electron collisions with disk : at such time moments the radial component of electron velocity is reduced by a factor @xmath182 , the reduction is done only if the kinetic energy of electron is larger than the fermi energy .
such a dissipation can be induced by phonon excitations inside the antidot disk .
we fix geometry directing _ dc-_field along @xmath146axis and microwave along @xmath183axis .
the noise is modeled in the same way as above in eq .
( [ eq1 ] ) .
we call this system disk model ( d1 ) .    , @xmath57 .
_ top right panel _ : temporary captured path at @xmath152 , @xmath57 . _
bottom left panel _ : path captured forever at @xmath152 , @xmath57 .
_ bottom right panel _
: no capture at @xmath152 , @xmath106 .
the trajectory part colliding with disk is shown by black curve , its part before and after collisions is shown in gray .
the red ( light gray ) points and curves show the trajectory of cyclotron center . here
the dissipation parameter is @xmath184 ; the static electric field is directed along @xmath146axis and @xmath185 ; microwave field is directed along @xmath183axis .
there is no noise here .
coordinates @xmath186 are expressed in units of @xmath160 .
data for model ( d1 ) . ]    examples of electron cyclotron trajectories scattering on disk are shown in fig .
[ fig6 ] . in absence of microwave field
a trajectory escapes from disk rather rapidly .
a similar situation appears at @xmath106 and microwave field with @xmath152 .
in contrast for @xmath57 and @xmath152 a trajectory can be captured for a long time or even forever depending on initial impact parameter .    for some impact parameters
a trajectory can be captured for a very long time @xmath187 , in certain cases in absence of noise we have @xmath188 . at such long capture times
the collisions with disk become synchronized with the phase @xmath98 of microwave field at the moment of collisions . this
is directly illustrated in fig .
[ fig7 ] where we show the angle @xmath189 of a collision point on disk , counted from @xmath146axis , in dependence on @xmath98 . indeed , the dependence @xmath189 on @xmath98 forms a smooth curve corresponding to synchronization of two phases .
at the same time the radial velocity at collisions @xmath156 moves along some smooth invariant curve @xmath190 in the phase space @xmath191 . however , to make a correct comparison with the radial field models ( dr1 ) , ( dr2 ) we should take into account that the cyclotron circle rotates around disk so that we should draw the poincar section in the rotational phase @xmath192 . in this representation
we see the appearance of the resonance ( see right column of fig .  [ fig7 ] ) that is similar to those seen in figs .
[ fig4 ] ,  [ fig5 ] for the radial field models .     and microwave phase @xmath98 .
_ left column _ : dependence of angle @xmath189 of collision point on disk , counted from @xmath146axis , and radial velocity @xmath156 , taken at collision , on microwave phase @xmath98 .
_ right column _ : same as in left column with @xmath192 .
here @xmath193 , @xmath152 , @xmath194 , @xmath184 , @xmath195 , there is no noise ; points are shown for times @xmath196 , the capture time of this orbit is @xmath197 .
data for model ( d1 ) . ]    a more direct correspondence between radial field models ( dr1 ) , ( dr2 ) and the model ( d1 ) with a linearly polarized microwave field is well seen from the poincar sections shown in the rotation frame of phase @xmath192 in fig .
[ fig8 ] . in this frame
we see directly the resonance at @xmath198 being very similar to the wall case and the radial field model . however , the positions of resonance at @xmath199 are different from those in fig .
[ fig5 ] . of course in the rotation frame
the orbital momentum is only approximately conserved that gives a broadening of invariant curves in fig .
[ fig8 ] .
with @xmath192 for @xmath168 , @xmath200 ( left top ) ; @xmath193 , @xmath200 ( right top ) ; @xmath201 , @xmath76 ( left middle ) ; @xmath193 , @xmath76 ( right middle ) ; @xmath201 , @xmath152 ( left bottom ) ; @xmath193 , @xmath152 ( right bottom ) ; @xmath156 is expressed in units of @xmath19 .
data for model ( d1 ) , no noise no dissipation . ]
we explain this as follows . for the linear polarized field of model ( d1 )
the radial component of microwave field is proportional to @xmath202 where we kept only slow frequency component of radial field ( the neglected term with @xmath203 gives resonant values @xmath204 ) .
the radial field @xmath205 gives kicks to the radial velocity component at collisions with disk similar to the case of model ( dr2 ) described by eq .
( [ eq3 ] ) : @xmath206 , @xmath207 . here
we use the radial field component phase @xmath208 at a moment of collision with disk ( the tangent component does not change @xmath156 and can be neglected ) .
the phase variation @xmath209 has the first term @xmath210 being the same as for the radial field model ( dr2 ) , and an additional term related to rotation around disk with @xmath211 which comes from geometry .
indeed , the segment angles of intersections of circles @xmath160 and @xmath161 are : for disk radius @xmath160 it is @xmath212 and for cyclotron radius @xmath161 it is @xmath213 .
thus , their ratio is @xmath214 in agreement with the geometrical scaling .
this result can be obtained from the expression for @xmath215 by interchange of two disks that gives the above expression for @xmath216 ( at @xmath217 both shifts @xmath218 and @xmath219 are equal ) .    thus again the dynamical description is reduced down to the chirikov standard map with slightly modified parameters giving us for the model ( d1 ) the chaos parameter @xmath220 being usually smaller than unity , resonance position @xmath167 , resonance width @xmath173 and the resonance energy width @xmath174 : @xmath221 where @xmath176 is a shift of resonance produced by a finite separatrix half width @xmath222 . for our numerical simulations we have @xmath178 with @xmath223 , @xmath224 and @xmath225 .    at @xmath226 eq .
( [ eq5 ] ) gives the values @xmath227 at @xmath168 while the numerical data of fig .
[ fig8 ] give @xmath228 , and we have at @xmath193 the theory value @xmath229 being in a good agreement with the numerical value @xmath230 of fig .  [ fig8 ] .
for @xmath201 we have the resonance position at @xmath231 corresponding to the bulk and thus the resonance is absent .
the resonance width in fig .
[ fig8 ] at @xmath193 , @xmath200 can be estimated as @xmath232 that is in a satisfactory agreement with the theoretical value @xmath233 from ( [ eq5 ] ) .
we remind that in model ( d1 ) we have only approximate conservation of orbital momentum that gives a broadening of invariant curves and makes determination of the resonance width less accurate . in spite of this broadening
we see that the resonance description by the chirikov standard map works rather well .    in fig .
[ fig9 ] we show the dependence of average capture time @xmath187 on @xmath9 in the model ( d1 ) .
the averaging is done over @xmath234 trajectories scattered on disk at all such impact parameters that cyclotron orbit can touch the disk . here
, we show the ratio @xmath235 where @xmath236 is an average capture time in absence of microwave . according to our numerical data we have approximate dependence @xmath237 corresponding to a period of nonlinear oscillations in a disk map description discussed in @xcite .
on @xmath9 at @xmath152 shown at various static fields : @xmath238 ( top left ) , @xmath239 ( top right ) , @xmath240 ( bottom left ) , @xmath241 ( bottom right ) . here
@xmath184 , there is no noise .
data for model ( d1 ) . ]    the data of fig .
[ fig9 ] show a clear periodic dependence of capture time @xmath187 on @xmath9 corresponding to the periodicity variation of poincar section with @xmath9 ( see figs .  [ fig4 ] ,  [ fig5 ] .  [ fig8 ] ) .
this structure is especially well visible at weak static fields . with an increase of @xmath181 this structure
is suppressed .
indeed , at large @xmath181 even without microwave field the trajectories can escape from disk as it discussed in @xcite and microwave field does not affect the scattering in this regime .
the distributions of capture times are shown in fig .
[ fig10 ] .
we clearly see that at resonant values of @xmath9 a microwave field leads to appearance of long capture times . for example
, we have the probability to be captured for @xmath242 being @xmath243 at @xmath193 while at @xmath106 we have @xmath244(left panel in fig .
[ fig10 ] ) ; and we have @xmath245 at @xmath246 while at @xmath247 we have @xmath248 ( right panel in fig .
[ fig10 ] ) .
these data confirm much stronger capture at certain resonant values of @xmath9 .     of capture times @xmath187 for parameters of fig .
[ fig9 ] . _ left panel : _
@xmath106 ( dashed curve for minimum of capturing probability ) and @xmath193 ( full curve at maximum of capturing probability ) at @xmath194 . _ right panel : _ similar cases at @xmath247 ( dashed curve ) and @xmath246 ( full curve ) at @xmath249 .
data are obtained with @xmath250 trajectories started at different impact parameters and running up to time @xmath251 . here
@xmath187 is expressed in units of @xmath252 .
data for model ( d1 ) . ]    the data of figs .
[ fig9 ] ,  [ fig10 ] show that the scattering process on disk is strongly modified by a microwave field . however , to determine the conductivity properties of a sample we need to know what is an average displacement @xmath253 along static field after a scattering on a single disk .
indeed , in our model a dissipation is present only during collisions with disk while in a free space between disks the dynamics is integrable and hamiltonian .
hence during such a free space motion there is no displacement along the static field ( the dissipative part of conductivity or resistivity appears only due to dissipation on disk ) . the dependence of @xmath253 on @xmath9 is shown in fig .  [ fig11 ] . in absence of microwave field at @xmath254
we find @xmath255 that corresponds to a simple estimate @xmath256 .
the numerical data show that @xmath253 is practically independent of @xmath181 and that @xmath257 in absence of dissipation at @xmath258 . in presence of microwave field
we see that the displacement along the static field has strong periodic oscillations with @xmath9 .
the striking feature of fig .
[ fig11 ] is the appearance of windows of zero displacement @xmath259 at resonance values @xmath260 .
we discuss how this scattering on a single disk modifies the resistance of a sample with large number of disks in next section .     along static field after scattering on a single disk shown as a function of @xmath9 at various amplitudes of static and microwave fields : @xmath261 , @xmath262 ( gray dashed curve ) ; @xmath194 , @xmath262 ( black dashed curve ) ; @xmath261 , @xmath152 ( gray full curve ) ; @xmath194 , @xmath152 ( black full curve ) .
the data are obtained by averaging over @xmath263 scattered trajectories with random impact parameters ; here @xmath184 , noise amplitude @xmath264 .
data for model ( d1 ) . ]
to determine a resistance of a sample with many disks we use the following scattering disk model .
the scattering on a single disk in a static electric field @xmath181 is computed as it is described in the previous section with a random impact parameter inside the collision cross section @xmath265 .
after that a trajectory evolves along @xmath266-axis according to the exact solution of hamiltonian eq .
( [ eq1 ] ) ( no dissipation and no noise ) up to a collision with next disk which is taken randomly on a distance between @xmath267 and @xmath268 where @xmath269 is a mean free path along @xmath183axis and @xmath270 is a two - dimensional density of disks ( of course @xmath271 ) . in a vicinity of disk
the dynamical evolution is obtained by runge - kutta solution of dynamical equations as it was the case in previous section .
we use low disk density with @xmath272 .
the collision with disk is done with a random impact parameter in the @xmath146axis of disk vicinity : the impact parameter is taken randomly in the interval @xmath273 $ ] around disk center .
noise acts only when a center of cyclotron radius of trajectory is on a distance @xmath274 from disk center so that a collision with disk is possible . after scattering on a disk
a free propagation follows up to next collision with disk .
along such a trajectory we compute the average displacement @xmath275 and @xmath276 after a time interval @xmath277 . in this way
the number of collisions with disks is @xmath278 and a total displacement in @xmath146axis is @xmath279 where @xmath253 is an average displacement on one disk discussed in previous section .
we compute the global displacements @xmath280 on a time interval @xmath281 averaging data over 200 trajectories .
we call this system model ( d2 ) .     and @xmath282 on magnetic field @xmath283 ( resistivity is expressed in arbitrary numerical units ) .
blue and gray curves show respectively @xmath120 and @xmath282 in absence of microwave field .
black curve with points show @xmath120 dependence of @xmath284 at microwave field @xmath152 ( here @xmath284 is expressed in units of @xmath285 ) . here
@xmath194 , @xmath184 , noise amplitude @xmath264 , @xmath27 .
resonant values @xmath11 are shown by arrows ; bars show statistical errors for @xmath120 .
data are obtained by averaging over 200 trajectories propagating up to time @xmath286 .
data for model ( d2 ) . ]
then the current components are equal to @xmath287 , @xmath288 and conductivity components are @xmath289 , @xmath290 ( the current is computed per one electron ) .
we work in the regime of weak _
dc-_field where @xmath291 scales linearly with @xmath145 .
the current @xmath292 is determined by the drift velocity @xmath293 . since the mean free path is large compared to disk size @xmath294 we have an approximate relation @xmath295 , @xmath296 . as in 2deg experiments
@xcite we have in our simulations @xmath297 ( see fig .  [ fig12 ] ) .
the resistivity is obtained by the usual inversion of conductivity tensor with @xmath298 , @xmath299 .
the dependence of @xmath120 , @xmath282 , expressed in arbitrary numerical units , on magnetic field @xmath300 is shown in fig .
[ fig12 ] .    in absence of microwave field
we find @xmath301 and @xmath302 similar to experiments @xcite .
for small noise amplitude ( e.g. @xmath264 ) we have @xmath120 growing linearly with @xmath284 ( see fig .  [ fig12 ] ) but at larger amplitudes ( e.g. @xmath303 ) its increase with @xmath284 becomes practically flat showing only @xmath304 increase in a give range of @xmath284 variation .
( left panel ) and @xmath282 ( right panel ) as function of @xmath0 at various noise amplitudes @xmath264 ( black curve ) , @xmath305 ( blue / dark curve ) , @xmath306 ( red / gray curve ) .
here @xmath152 , @xmath194 , other parameters are as in fig .
[ fig12 ] .
curves are drown though numerical points obtained with a step @xmath307 .
data for model ( d2 ) . ]    in presence of microwave field the dependence of @xmath120 on @xmath284 is characterized by periodic oscillations with minimal @xmath120 values being close to zero at resonant values of @xmath308 well visible in fig .
[ fig12 ] . the dependence of @xmath120 , @xmath282 rescaled to their values @xmath309 ,
@xmath310 in absence of microwave field are shown in fig .
[ fig13 ] at various amplitudes of noise and fixed @xmath63 , and in fig .  [ fig14 ] at various @xmath63 and fixed noise amplitude @xmath139 .
we see that increase on noise leads to an increase of minimal values of @xmath120 at resonant values @xmath11 . in a similar way
a decrease of microwave power leads to increase of minimal values of @xmath120 at @xmath11 . at the same time
the hall resistance @xmath282 is only weakly affected by microwave radiation as it also happens in zrs experiments .     at various microwave amplitudes @xmath200 ( red / gray curve ) ,
@xmath306 ( blue / dark curve ) , @xmath311 ( black curve ) ; amplitude of noise is fixed at @xmath264 .
data for model ( d2 ) . ]
these results are in a qualitative agreement with the zrs experiments . on the basis of our numerical studies
we attribute the appearance of approximately zero resistance at @xmath11 values in our bulk model of disk scatterers to long capture times of orbits in disk vicinity at these @xmath11 values ( see fig .  [ fig9 ] ) . during this time @xmath187 noise
gives fluctuations of collisional phase @xmath189 and due to that a cyclotron circle escapes from disk practically at random displacement @xmath253 that after averaging gives average @xmath257 .
since resistivity is determined by the average value of @xmath253 this leads to appearance of zrs .
we note that this mechanism is different from the one of edge transport stabilization discussed here and in @xcite . however , both mechanisms are related to a long capture times near edge or near disk that happens due to synchronization of cyclotron phase with microwave field phase and capture inside the nonlinear resonance .
of trajectories at the moment of collisions with disks for parameters of fig .
[ fig14 ] at @xmath152 : @xmath168 ( top left ) , @xmath193 ( top right ) , @xmath201 ( bottom left ) , @xmath312 ( bottom right ) . here
@xmath192 where @xmath91 is a microwave phase at the moments of collisions with disk and @xmath189 is the angle on disk at collision moment , counted from @xmath47-axis ( same as in figs .
[ fig7 ] ,  [ fig8 ] ) .
data are obtained from @xmath313 trajectories iterated up to time @xmath286 .
density of points is shown by color with black at zero and white at maximum density .
the average number of collisions per disk per trajectory is @xmath314 , @xmath315 , @xmath316 , @xmath317 respectively for @xmath168 , @xmath318 , @xmath319 , @xmath320 .
data are obtained for model ( d2 ) . ]    to illustrate the capture inside the resonance we present the distributions of trajectories from fig .
[ fig14 ] shown in the phase space plane @xmath321 at the moments of collisions with disks in fig .
[ fig15 ] .
this is similar to the poincar sections of fig .
[ fig8 ] however , now we consider the real case of diffusion and scattering on many disks in the model ( d2 ) with noise and dissipation .
we see that for @xmath193 orbits are captured in a vicinity of the center of nonlinear resonance at @xmath322 well seen in fig .
for @xmath168 we have a density maximum located at smaller values of @xmath167 and @xmath322 even if there is a certain shift of @xmath167 produced by a significant resonance width at @xmath152 . at @xmath201
we have a density maximum at @xmath323 corresponding to an unstable fixed point of separatrix .
the total number of collision points @xmath324 in this case is by a factor @xmath325 smaller than in the case of stable fixed point at @xmath193 .
a similar situation is seen in the case of wall model ( w1 ) ( see fig.1d , f in @xcite ) even if there the ratio between number of captured points was significantly larger .
the results of fig .
[ fig15 ] show that in the zrs phase the collisions with disk indeed create synchronization of cyclotron and microwave phases and capture of trajectories inside the nonlinear resonance .    however , there are also some distinctions between bulk disk model ( d2 ) and experimental observations .
the first one is that there are minima for @xmath326 but there are no peaks which are well visible around integer @xmath9 values in zrs experiments @xcite,@xcite and numerical simulations of transport along the edge @xcite .
the second one is appearance of small negative values of @xmath120 at @xmath11 values .
on @xmath9 in the disk model with dissipation at disk collisions at rate @xmath327 and dissipation in disk vicinity with rate @xmath328 ; a disk roughness gives additional angle rotations with amplitude @xmath329 ( see text ) ; the amplitude of noise in disk vicinity is @xmath330 .
here we have @xmath152 , @xmath194 ; 51 numerical points in @xmath9 are connected by lines to adapt an eye .
data are obtained by averaging over 100 trajectories propagating up to time @xmath286 .
data for model ( d3 ) . ]
we attribute the absence of peaks to a specific dissipation mechanism which takes place only at disk collisions . it is rather convenient to run long trajectories using exact solution for free propagation between disks .
indeed , in this scheme there is no dissipation during this free space propagation and thus these parts of trajectories have no displacement along static field .
we also tested a dissipation model with additional @xmath331 for @xmath332 and @xmath333 for @xmath334 .
this dissipation works only in a disk vicinity when the distance between disk center and cyclotron center is smaller than @xmath335 .
the dissipation @xmath336 on disk remains unchanged .
we also added a certain roughness of disk surface modeled as an additional random angle rotation of velocity vector in the range @xmath337 , done at the moment of collision with disk .
we call this system disk model ( d3 ) .
the results for the resistivity ratio @xmath326 are shown in fig .
[ fig16 ] .
they show an appearance of clear peaks of @xmath326 in presence of such additional dissipation in vicinity of integer @xmath9 .
there is also a small shift of minima from integer plus @xmath45 to integer plus @xmath338 .
the second point of distinctions from zrs experiments is a small negative value of @xmath120 at resonant @xmath9 values .
it is relatively small for disk model ( d2 ) ( see figs .
[ fig13 ] ,  [ fig14 ] ) and it becomes more pronounced for disk model ( d3 ) ( see fig .
[ fig16 ] ) .
it is possible that a scattering on disk in presence of dissipation , noise , static and microwave field gives a negative displacement @xmath253 which generates such negative @xmath120 values .
we expect that in the limit of static field going to zero this effect disappears .
indeed , the negative values become smaller at smaller @xmath181 according to data of fig .
[ fig11 ] but unfortunately the small @xmath181 limit is also very difficult to investigate numerically .
we consider that at this stage of the theory the presence of negatives values for @xmath120 does not constitute a critical disagreement .
the escape parameters for electrons that have been captured on an impurity for a time long enough to make many rotations around it , are likely to strongly depend on the model for the electron impurity interaction and further theoretical work on a more microscopic model is needed . in general
a zero average displacement along the field direction seems natural for a smooth distribution of trapping times with a characteristic time scale much larger than the rotation time around the impurity ( this assumption does not seem to hold for our model , see for example the sharp features on fig .
finally in section vii we propose a slightly different mechanism by which the combination of trapping on impurities investigated here and electron - electron interactions can lead to zrs .
the zrs experiments @xcite show that the resistance @xmath120 in the zrs minima scales according to arrehenius law @xmath339 with a certain energy scale dependent on a strength of microwave field . in typical experimental conditions
one finds very large @xmath340 at @xmath341 ( see e.g. fig.3 in @xcite ) .
these data also indicate the dependence @xmath342 .
this energy scale @xmath343 is very large being only by a factor @xmath344 smaller than the fermi energy @xmath345 . at the same time
the amplitude of microwave field is rather weak corresponding to @xmath346 at field of @xmath347 or ten times larger at @xmath348 ( unfortunately it is not known what is an amplitude of microwave field acting on an electron ) .
as in @xcite we argue that the arrehenius scale is determined by the energy resonance width ( [ eq4 ] ) with @xmath349 .
indeed , the resonance forms an energy barrier for a particle trapped inside the resonance by dissipative effects being analogous to a wash - board potential .
an escape from this potential well requires to overcome the energy @xmath350 leading to the arrehenius law for @xmath120 dependence on temperature . assuming the case of the wall with @xmath12 we obtain at @xmath351 the activation temperature @xmath352 being is a satisfactory agreement with the experimental observations . the theoretical relation ( [ eq4 ] ) also reproduces the experimental dependence @xmath353 at @xmath12 . in his relation @xmath354
being confirmed by the numerical simulations presented in @xcite .
this dependence is in a satisfactory agreement with the power dependence found in experiments @xcite . in other samples
one finds that the dependence @xmath355 works in a better way .
we think that higher terms in a nonlinear resonance can be responsible for scaling @xmath355 being different from the relation ( [ eq4 ] ) .
also a finite rigidity of the wall or disk scatterers can be responsible for appearance of higher power of @xmath63 .
the energy scale @xmath350 on disks is enhanced by a factor @xmath94 for the case of radial field ( [ eq4 ] ) .
however , we showed that for a linear polarization the scale @xmath350 is given by eq.([eq5 ] ) and thus there is no enhancement at large @xmath92 .
indeed , we performed direct simulations at parameters of fig .
[ fig14 ] with the reduced value of disk radius by a factor 2 .
the numerical data give approximately the same traces @xmath326 vs. @xmath9 at @xmath356 without visible signs of deeper minima at small @xmath63 .
this confirms the theoretical expressions ( [ eq5 ] ) . in any case , for small values @xmath158 one should analyze the quantum scattering problem which is significantly more complex compared to the classical case .
we may assume that in a quantum case one should replace @xmath160 by a magnetic length @xmath357 .
in such a case we are getting @xmath358 that gives @xmath359 at @xmath360 and microwave amplitude @xmath361 .
however , in this case we obtain the scale @xmath362 being practically independent of @xmath9 which differs from experimental data . in any case in experiments
the size of impurities is small compared to @xmath161 and a quantum treatment is required to reproduce the correct picture for @xmath120 dependence on parameters in the zrs phase .
another point is related to the positions of zrs minima on @xmath9 axis .
we remind that that for the wall model the resonance is located at @xmath363 ( [ eq4 ] ) and that the separatrix width is @xmath364 .
the capture of trajectories from the bulk is most efficient when a half width of separatrix touch the border of bulk at @xmath365 with @xmath366 that gives the expression @xmath367 for the wall case . at @xmath368 , @xmath193
this gives @xmath369 being in a good agreement with the numerical data @xmath370 for @xmath120 dependence on @xmath9 ( see figs.2,3 in @xcite with a visible tendency of @xmath103 growth with @xmath9 ) .
for the data presented here in fig .
[ fig14 ] for the disk case at @xmath152 , @xmath178 , @xmath193 we obtain from ( [ eq5 ] ) @xmath371 that is slightly less than the numerical value @xmath372 for minima location .
we attribute this difference to an approximate nature of expression for the resonance width at relatively strong microwave fields .
we also note that in experiments an additional contribution to the value of @xmath103 can appear due to a finite rigidity of disk and wall potentials .
the theoretical models presented here and in @xcite reproduce main experimental features of zrs experiments @xcite,@xcite
. however , it would be useful to have some additional theoretical predictions which can be tested experimentally .
a certain characteristic feature of both wall and disk models is appearance of nonlinear resonance .
for example , according to the wall model ( w2 ) described by eq .
( [ eq3 ] ) the dynamics inside the resonance is very similar to dynamics of a pendulum .
the frequency of phase oscillations inside the resonance is @xmath373 @xcite . here
the frequency is expressed in number of map iterations and since the time between collisions is approximately @xmath82 we obtain the physical frequency @xmath374 of these resonant oscillations being @xmath375 . at @xmath376
this frequency is significantly smaller than the driving microwave frequency .
the dynamics inside the resonance should be very sensitive to perturbations at frequency @xmath377 that gives : @xmath378    to check this theoretical expectation we study numerically the effect of additional microwave driving with dimensional amplitude @xmath379 ( @xmath380 ) and frequency @xmath381 .
we use the wall model ( w2 ) based on the chirikov standard map described here and in @xcite . as for the additional driving frequency
we have @xmath382 where @xmath383 is the field strength of microwave frequency @xmath381 , we assume that both main and additional fields are collinear and perpendicular to the wall . in presence of second frequency the map ( [ eq3 ] ) takes the form : @xmath384 the only modification appears in the first equation since now the change of velocity at collision depends on both fields ; the second equation remains the same as in ( [ eq3 ] ) .
as in the model ( w2 ) the term @xmath385 describes the effects of dissipation with rate @xmath97 and noise with amplitude @xmath386 of random velocity angle rotations .
we call this system model ( w3 ) .     on frequency ratio @xmath387 in model ( w3 ) described by the map ( [ eq7 ] ) .
here , the test driving at frequency @xmath381 has fixed amplitude @xmath388 ; the main microwave driving is located in the zrs phase at @xmath389 , and its amplitude takes values @xmath390 ( blue curve ) , @xmath391 ( green curve ) , @xmath311 ( magenta curve ) , @xmath392 ( red curve ) ( these curves follow from top to bottom at @xmath393 ) .
the values of @xmath120 are computed at fixed @xmath388 and corresponding @xmath63 ; the values of @xmath394 are computed at @xmath395 and @xmath76 .
the data are obtained at noise amplitude @xmath396 and dissipation @xmath397 ; averaging is done over 2000 trajectories for 5000 iterations of map ( [ eq7 ] ) .
_ right panel : _ dependence of peak position @xmath387 on main microwave driving amplitude @xmath63 obtained from data of left panel at @xmath388 and additional data at @xmath398 ( blue points ) , the theory dependence ( [ eq6 ] ) at @xmath399 is shown by the straight red line .
data for model ( w3 ) .
, title="fig : " ]   on frequency ratio @xmath387 in model ( w3 ) described by the map ( [ eq7 ] ) . here , the test driving at frequency @xmath381 has fixed amplitude @xmath388 ; the main microwave driving is located in the zrs phase at @xmath389 , and its amplitude takes values @xmath390 ( blue curve ) , @xmath391 ( green curve ) , @xmath311 ( magenta curve ) , @xmath392 ( red curve ) ( these curves follow from top to bottom at @xmath393 ) .
the values of @xmath120 are computed at fixed @xmath388 and corresponding @xmath63 ; the values of @xmath394 are computed at @xmath395 and @xmath76 .
the data are obtained at noise amplitude @xmath396 and dissipation @xmath397 ; averaging is done over 2000 trajectories for 5000 iterations of map ( [ eq7 ] ) . _ right panel : _ dependence of peak position @xmath387 on main microwave driving amplitude @xmath63 obtained from data of left panel at @xmath388 and additional data at @xmath398 ( blue points ) , the theory dependence ( [ eq6 ] ) at @xmath399 is shown by the straight red line .
data for model ( w3 ) .
, title="fig : " ]    in the model ( w3 ) the resistance @xmath120 is computed numerically in the same way as in the model ( w2 ) described in @xcite : the displacement along the edge between collisions is @xmath400 ; it determines the total displacement @xmath253 along the edge during the total computation time @xmath401 ; then @xmath402 where @xmath403 is an effective diffusion rate along the edge . to see the effect of additional weak test driving @xmath379 at frequency @xmath381 we
place the system in the zrs phase at @xmath389 and measure the variation of rescaled resistance @xmath404 . here
@xmath120 is the resistance in presence of both microwave fields @xmath63 and @xmath379 while @xmath394 is the resistance at @xmath395 and a certain fixed @xmath63 .
the dependence of @xmath404 on the frequency ratio is shown in fig .
[ fig17 ] at left panel .
the main feature of this data is appearance of peak at low frequency ratio @xmath405 . in the range
@xmath406 the testing field @xmath379 is nonresonant and does not affect @xmath120 however at @xmath405 it becomes resonant to the pendulum oscillations in the wall vicinity and hence strongly modifies @xmath120 value .
the dependence of this resonance ratio @xmath387 on amplitude of main driving field @xmath63 is shown in right panel of fig .
[ fig17 ] .
the numerical data are in a good agreement with the above theoretical expression ( [ eq6 ] ) .
the theoretical dependence ( [ eq6 ] ) allows to check the synchronization theory of edge state stabilization .
it also allows to measure the strength of main microwave driving force acting on an electron that still remains an experimental challenge .
the experimental testing of relation ( [ eq6 ] ) requires to work with good zrs samples which have every low resistance in zrs minima since this makes the effect of testing field @xmath381 to be more visible .
we note that the recent experiments in a low frequency regime @xmath407 @xcite demonstrate that @xmath120 is sensitive to low frequency driving .
the expression ( [ eq6 ] ) is written for the case when @xmath120 is mainly determined by a transport along edges .
if the dominant contribution is given by bulk disk scatterers then a certain numerical coefficient @xmath408 should be introduced in the right part of the expression . according to the data of fig .
[ fig8 ] and eqs .
( [ eq5 ] ) we estimate @xmath409 ( the separatrix width is smaller for the disk case compared to the wall case at the same @xmath63 ) .
another interesting experimental possibility of our theory verification is to take a hall bar of high mobility 2deg sample and put on it antidots with regular or disordered distribution ( it is important to have no direct collision - less path for a cyclotron radius in crossed _ dc-_electric and magnetic fields ) with a low density of antidot disks @xmath410 ( as in our numerical studies ) so that an average distance between antidots is larger than the cyclotron radius @xmath161 .
the regular antidot lattices have been already realized experimentally @xcite .
the effect of microwave field on electron transport in a regular lattice has been studied in the frame of ratchet transport in asymmetric lattices @xcite .
even a case of symmetric circular antidots has been studied in @xcite but the lattice was regular and no special attention was paid to analysis of resistivity at zrs resonant regime with @xmath136 .
we think that the experimental condition of @xcite can be relatively easy modified to observe the zrs effect on disk scatterers discussed here .
above we presented theoretical and numerical results which in our opinion explain the appearance of microwave induced zrs in high mobility samples .
the synchronization theory of zrs proposed in @xcite and extended here is based on a clear physical picture : high harmonics @xmath411 are generated by collisions with sharp edge boundary or isolated impurities which are modeled here by specular disks .
the zrs phases appear in a vicinity of resonant values @xmath412 . at these @xmath11 values
the cyclotron phase of electron motion becomes synchronized with the microwave phase due to dissipative processes present in the system .
for trajectories at the edge vicinity this synchronization gives stabilization of propagation along edge channels that creates an exponential drop of resistivity contribution of these channels with decreasing amplitude of thermal noise and increasing amplitude of microwave field .
the contribution to resistivity from trajectories in the bulk is analyzed in the frame of scattering on many well separated disk impurities . here again the synchronization of cyclotron phase with the microwave phase takes place approximately at the same resonant @xmath11 values . at these @xmath11 values
the synchronization leads to long time capture of trajectories in disk vicinity . during this long time
an initial cyclotron phase is washed out by noisy fluctuations and many rotations around disk and thus an electron escapes from a disk with an average zero displacement along the applied _
dc-_field even if dynamics in disk vicinity is dissipative .
this provides the main mechanism of suppression of dissipative resistivity contribution from isolated impurities in the bulk . as a result
the contribution of bulk to dissipative conductivity @xmath413 is suppressed , as it was assumed in @xcite , and the main contribution to current is given by electron propagation along edge states stabilized by a microwave field .
as we showed above the resonance width or resonance energy scale @xmath350 are approximately the same for the disk and wall cases ( see eqs .
( [ eq4 ] ,  [ eq5 ] ) ) .
we note that for the disk case the energy @xmath350 is not sensitive to the disk radius as soon as it is significantly smaller than the cyclotron radius .
thus we expect that at @xmath11 values the conductivity @xmath413 in the bulk is suppressed by a microwave field and at these fields the current is flowing essentially along stabilized edge states . in the case of corbino geometry
we have radial conductivity @xmath414 which is determined by the bulk scattering and now the minima of @xmath414 are located at @xmath11 values ( see e.g figs .  [ fig12 ] , [ fig13],[fig14 ] where @xmath415 ) .
the zrs experiments performed in the corbino geometry give minima of @xmath414 at these @xmath11 values ( see e.g. @xcite ) being in agreement with the synchronization theory .
it is interesting to note that the nonlinear dynamics in vicinity of edge and disk impurity is well described by the chirikov standard map @xcite .
the map description explains the location of resonances at integer values of @xmath9 with an additional shift @xmath370 produced by a finite separatrix width of nonlinear resonance .
a finite rigidity of wall or disk potential can give a modification of this shift @xmath103 .
our results show that the zrs phases at @xmath11 appear only at weak noise corresponding to high mobility samples .
strong noise destroys synchronization and trajectories are no more captured at edge or disk vicinity .
we also note that internal sample potentials with significant gradients act like a strong local _ dc-_field which destroys stability regions around disk impurities or near edge . due to
that the zrs effect exists only in high mobility samples .
the resistance at zrs minima drops significantly with the growth of microwave field strength since it increases the amplitude of nonlinear resonance which captures the synchronized trajectories .
the synchronization theory of zrs is based on classical dynamics of noninteracting electrons .
it is possible that electron - electron interaction effects can also suppress the contribution to resistivity from neutral short range range scatterers ( interface roughness , adatoms , ... ) .
indeed , long capture times can increase the electron density around these short ranged impurities transforming them into long range charged scatters that the other electrons can circumvent by adiabatically following the long range component of the disorder potential thereby avoiding a scattering event .
however , the theoretical description of this short - ranged impurity cloaking mechanism for the zrs effect remains a serious challenge .
another important step remains the development of a quantum synchronization theory for zrs . even if in experiments the landau level is relatively high @xmath416 , there are only about ten landau levels inside a nonlinear resonance @xcite and quantum effects should play a significant role .
the general theoretical studies show that the phenomenon of quantum synchronization persists at small effective values of planck constant @xmath417 but it becomes destroyed by quantum fluctuations at certain large values of @xmath417 @xcite .    the importance of quantum zrs theory is also related to a short range nature of the impurities considered here , typically on a scale of a few nanometers or even less . we have modeled these impurities by disks with a radius that was only several times ( in fact @xmath9 times ) smaller than the cyclotron radius which is not so close to microscopic reality .
we could argue that in the quantum case a nanometer sized impurity would act effectively as an impurity of a size of quantum magnetic length @xmath418 .
this gives a ratio @xmath419 which is comparable with the one used in our simulations with @xmath420 but of course a quantum treatment of scattering on nanometer size impurity in crossed electric and magnetic and also microwave fields remains a theoretical challenge .
we note that such type of scattering can be efficiently analyzed by tools of quantum chaotic scattering ( see e.g. @xcite ) and we expect that these tools will allow to make a progress in the quantum theory development of striking zrs phenomenon .
we hope that the synchronization theory of microwave induced zrs phenomenon described here can be tested in further zrs experiments .
this work was supported in part by anr france pnano project nanoterra ; ovz was partially supported by the ministry of education and science of russian federation .
mani , j.h .
smet , k. von klitzing , v.  narayanamurti , w.b .
johnson , and v.  umansky , nature * 420 * , 646 ( 2002 ) .
zudov , r.r .
du , l.n .
pfeiffer , and k.w .
west , phys .
lett . * 90 * , 046807 ( 2003 ) .
dorozhkin , jetp lett . * 77 * , 577 ( 2003 ) j. h. smet , b. gorshunov , c. jiang , l.  pfeiffer , k.  west , v.  umansky , m.  dressel , r.  meisels , f.  kuchar , and k. von klitzing , phys .
lett . * 95 * , 116804 ( 2005 ) a.a .
bykov , a.k .
bakarov , d.r .
islamov , and a.i .
toropov , jetp lett . *
84 * , 391 ( 2006 ) d.  konstantinov , and k.  kono , phys .
lett . * 105 * , 226801 ( 2010 ) d.konstantinov , a.d.chepelianskii , and k.kono , jour .
japan * 81 * , 093601 ( 2012 ) i.a .
dmitriev , a.d .
mirlin , d.g .
polyakov , and m.a .
zudov , rev .
mod . phys . * 84 * , 1709 ( 2012 ) a.d .
chepelianskii , and d.l .
shepelyansky , phys .
b * 80 * , 241308(r ) ( 2009 ) a.  pikovsky , m.  rosenblum , and j.  kurths , _ synchronization : a universal concept in nonlinear sciences _ , cambridge univ .
press , cambridge ( 2001 ) a.d .
chepelianskii , _ non linear transport in drift - diffusion equations under magnetic field _ , + arxiv:1110.2033[cond - mat.mes - hall ] ( 2011 ) w.  g.  hoover , _ time reversibility , computer simulation , and chaos _ , world sci . , singapore ( 1999 ) b.v .
chirikov , phys . rep . * 52 * , 263 ( 1979 ) a.j .
lichtenberg , and m.a .
lieberman , _ regular and chaotic dynamics _ , springer , berlin ( 1992 ) .
b.chirikov , and d.shepelyansky , scholarpedia * 3(3 ) * , 3550 ( 2008 ) y.  avishai , and g.  montambaux , eur .
j. b * 79 * , 215 ( 2011 ) f.  benvenuto , g.  casati , i.  guarneri , and d.  l.  shepelyansky , z. phys .
b cond . mat . *
84 * , 159 ( 1991 ) d.  weiss , m.  l.  roukes , a.  menschig , p.  grambow , k. von klitzing , g.  weimann , phys .
* 66 * , 2790 ( 1991 ) g.  m.  gusev , v.  t.  dolgopolov , z.  d.  kvon , a.  a.  shashkin , v.  m.  kudryashov , l.  v.  litvin and yu .
nastaushev , pisma zh .
. fiz . * 54 * , 369 ( 1991 ) r.  fleischmann , t.  geisel , and r.  ketzmerick , phys .
lett . * 68 * , 1367 ( 1992 ) t.  geisel , r.  ketzmerick and o.  schedletzky , phys .
rev . lett .
* 69 * , 1680 ( 1992 ) n.  berglund , a.  hansen , e.  h.  hauge , and j.  piasecki , phys .
lett . * 77 * , 2149 ( 1996 ) a.d .
chepelianskii , and d.l .
shepelyansky , phys .
b * 63 * , 165310 ( 2001 ) a.d .
chepelianskii , j.  laidet , i.  farrer , h.e .
beere , d.a .
ritchie , and h.  bouchiat , phys .
b * 86 * , 205108 ( 2012 ) s.  sassine , yu .  krupko , j .- c .
portal , z.d .
kvon , r.  murali , k.p .
martin , g.  hill , and a.d .
wieck , phys .
b * 78 * , 045431 ( 2008 ) c.l .
yang , m.a .
zudov , t.a .
knuuttila , r.r .
du , l.n .
pfeiffer , and k.w .
west , phys .
* 91 * , 096803 ( 2003 ) a.a .
bykov , jetp lett . * 87 * , 551 ( 2008 ) o.v .
zhirov , and d.l .
shepelyansky , eur .
j. d * 38 * , 375 ( 2006 ) j.  wiersig , and m.  hentschel , phys .
lett . * 100 * , 033901 ( 2008 ) a.  eberspcher , j.  main , and g.  wunner , phys .
e * 82 * , 046201 ( 2010 )