the abelian sandpile model and rotor - routing model are combinatorial models of various dynamics on graphs that were rediscovered several times in several different communities , ranging from combinatorics to self - organized criticality , arithmetic geometry , and algorithms .
( for example , @xcite ; see @xcite for more details and references . )
the stable configurations in the abelian sandpile model on a graph @xmath1 form a group called the sandpile group or the critical group of @xmath1 ; it is also known as the picard group of the graph , which we will denote @xmath0 in analogy with the picard group of a riemann surface .
it has long been known that the order of @xmath0 is equal to the number of spanning trees of @xmath1 ( see , e.g. , @xcite and references therein ) .
indeed , this is a form of kirchoff s matrix  tree theorem . a natural first question to ask , then ,
is whether there is a canonical bijection between the elements of @xmath0 and the spanning trees of @xmath1 .
what does `` canonical '' mean here ?
at the very least , we should require that the bijection be invariant under automorphisms of @xmath1 , for certainly it should nt depend on what we name the vertices and edges .
but there is clearly no such bijection : many graphs , e.g.  a complete graph or an @xmath2-cycle , do nt even have a distinguished spanning tree that could correspond to the identity element of @xmath0 .    so it is natural to ask instead for the next best thing : _
is there a canonical free , transitive action of the sandpile group on the set of spanning trees ? _ in other words , is there a canonical @xmath0-torsor structure on the set of spanning trees of @xmath1 ? by `` canonical , '' we still mean an action that is invariant under automorphisms of @xmath1 .
such an action would give a bijection between the sandpile group and the set of spanning trees , once we decide on a spanning tree that corresponds to the identity element .
but this is still too much to ask for .
for example , consider the graph @xmath1 with two vertices and @xmath2 parallel edges between them .
the sandpile group @xmath0 is cyclic of order @xmath2 , and each edge of @xmath1 is itself a spanning tree .
so an automorphism - invariant action of @xmath3 amounts to an automorphism - invariant cyclic ordering of the edges , which clearly does not exist .
this example , and others like it , suggest that we might hope for a free transitive action that is canonically defined on @xmath1 once we ve fixed the additional structure of a cyclic ordering of the edges incident to each vertex . a graph together with such a system of cyclic orderings is called a _ ribbon graph _ ( or sometimes `` combinatorial embedding '' ) .
we therefore arrive at the following question .
[ q : intro ] is there a free and transitive action of @xmath0 on the spanning trees of @xmath1 that is invariant under automorphisms preserving a ribbon graph structure on @xmath1 ?
holroyd _ et al . _
@xcite gave one possible way to answer question  [ q : intro ] .
given a ribbon graph _ and _ a choice of a basepoint vertex , they used rotor - routing to define a free and transitive action of @xmath0 on the spanning trees of @xmath1 .
this action is explained in section  2 .
therefore an affirmative answer to question  [ q : intro ] would follow if one could prove that the basepoint in the construction of holroyd _ et al _ is unnecessary .
in other words , we have the following question , which to our knowledge was first asked by ellenberg @xcite : does the rotor - routing action of @xmath0 on the spanning trees of @xmath1 depend on the choice of basepoint ?    in this paper we prove that if @xmath1 is a _ planar _ ribbon graph , then miraculously , the action defined by rotor routing is independent of the basepoint and hence canonical .
furthermore , we show that this characterization is tight : only planar ribbon graphs have this property .
our main theorem is thus :    [ theorem : main ] let @xmath1 be a connected ribbon graph .
the action of the sandpile group @xmath0 on the set @xmath4 of spanning trees is independent of the choice of basepoint vertex if and only if @xmath1 is a planar ribbon graph .
the proof is based on three key ideas .
first , the rotor - routing action of the sandpile group on spanning trees can be partially modeled via _ rotor - routing on unicycles _
( @xcite ) .
this is a related dynamical system with the property that rotor - routing becomes periodic , rather than terminating after finitely many steps .
rotor - routing on unicycles is described in section  [ sec : reversibility ] .
the second main idea is that the independence of the sandpile action on spanning trees can be described in terms of _ reversibility of cycles_. we introduce the notion of reversibility ( previously considered in @xcite only for planar graphs ) , and prove in proposition  [ prop : reversibility2 ] that reversibility is a well - defined property of cycles in a ribbon graph .
we establish a relation between reversibility and basepoint - independence in section  [ sec : spanningtrees ] .
third , reversibility is closely related to whether a cycle separates the surface corresponding to the ribbon graph into two components .
we prove in lemma  [ lemma : almostseparating ] that these conditions are almost equivalent .
moreover , although they are not equivalent for individual cycles , we prove in proposition  [ prop : planar_is_reversible ] that _ all _ cycles are reversible if and only if all cycles are separating , in which case the ribbon graph is planar .
by a _ graph _ we mean a connected and undirected graph , with multiple edges allowed but no self - loops .
we write @xmath5 and @xmath6 for the vertex set and edge set of @xmath1 , respectively . throughout
, @xmath1 denotes a _ ribbon graph _ : a graph together with a cyclic ordering of the edges incident to each vertex .
ribbon graphs are equivalent to graphs embedded into orientable surfaces ( hence their alternate name `` combinatorial embeddings '' ) : from any ribbon graph one can reconstruct a closed orientable surface in which it is naturally embedded , and conversely the local orientation of the surface at a vertex gives the cyclic ordering ( see , e.g. , ( * ? ? ?
* theorem  3.7 ) ) .
we denote by @xmath4 the set of spanning trees of @xmath1 .
when the graph is clear from context , we sometimes write @xmath7 for @xmath4 .    * the rotor - routing model . * a _ rotor configuration _
@xmath8 on a directed graph is an assignment to each vertex @xmath9 of an outgoing edge @xmath10 $ ] based at @xmath9 .
when we speak about rotor - routing or rotor configurations on a ( ribbon ) graph @xmath1 , we will always mean rotor configurations on the underlying eulerian directed graph @xmath11 obtained by replacing each undirected edge of @xmath1 with two oppositely oriented edges .    fix a ribbon graph @xmath1 .
the _ rotor - routing model _ is a deterministic process on pairs @xmath12 consisting of a rotor configuration @xmath8 and a vertex @xmath13 .
we think of @xmath14 as the position of a _ chip _ that will move around the graph as the process is iterated .
one step of the rotor - routing process consists of replacing @xmath12 by a new pair @xmath15 defined as follows .
first , the rotor configuration @xmath16 coincides with @xmath8 at every vertex except at @xmath14 , where @xmath17 $ ] is the edge leaving @xmath14 that follows @xmath18 $ ] in the cyclic order .
second , the new position of the chip @xmath9 is the vertex at the other end of the edge @xmath17 $ ] from @xmath14 . in short ,
the rotor - routing process rotates the rotor @xmath18 $ ] at the current position of the chip to the next position , and then moves the chip along this edge in the new direction of the rotor . in this paper
we will frequently iterate the rotor - routing process until a certain `` stopping condition '' is reached , and reason about the resulting configuration and the steps taken to reach it .    *
the sandpile group . *
the group @xmath19 of _ divisors _ on the graph @xmath1 is the free abelian group of formal @xmath20-linear combinations @xmath21 of the vertices of @xmath1 . in this paper
we will be mostly interested in the subgroup @xmath22 of `` degree-0 '' divisors : @xmath23 if we fix a vertex @xmath24 ( for `` root '' ) , @xmath22 is freely generated by the basis @xmath25 as @xmath26 ranges over the other vertices of @xmath1 .
the _ graph laplacian _ of a graph @xmath1 is an operator @xmath27 closely related to the adjacency matrix of @xmath1 .
with respect to the basis of vertices @xmath26 , the operator @xmath28 corresponds to the symmetric matrix whose diagonal entries are @xmath29 = degree of @xmath26 and whose off - diagonal entries @xmath30 .
the sum @xmath31 of each row and column is 0 , which shows both that the image of @xmath28 is contained in @xmath22 , and that @xmath32 .
the image of @xmath28 is a finite - index subgroup of @xmath22 , and the _ sandpile group _
@xmath0 is the quotient : @xmath33 for any graph @xmath1 the group @xmath0 is a finite abelian group whose order is , up to sign , the determinant @xmath34 of any @xmath35 submatrix of the graph laplacian .
people who study sandpile groups usually prefer to work with the distinguished representatives of each divisor class that are recurrent under a certain dynamical process
. we will not describe the beautiful story of sandpile dynamics here , but that viewpoint does implicitly play a key role in the rotor - routing action that we now define .
* action of @xmath0 on @xmath4 . * given a choice of root @xmath36 , holroyd _ et al . _
@xcite defined a natural action of @xmath22 on @xmath4 , which we now describe .
this action is trivial on any divisor in @xmath37 and thus descends to an action of @xmath0 on @xmath4 .    .
the rotor - routing process in the bottom row shows that @xmath38 and thus @xmath39 so the action of @xmath40 on @xmath41 is independent of our choice of @xmath14 or @xmath9 as basepoint . here
, the white circle indicates the position of the chip , and the rotors are all oriented clockwise relative to the page.,height=192 ]    for a fixed root @xmath24 , given a spanning tree @xmath42 and a divisor @xmath43 , we will write @xmath44 for the image of @xmath41 under the rotor - routing action ( based at @xmath24 ) .
we define the action of @xmath22 by describing the action @xmath45 of the generators @xmath46 for each @xmath47 , as follows . given a spanning tree @xmath41 , orient the edges of @xmath41 towards @xmath24 .
the resulting collection of directed edges defines a rotor configuration on @xmath1 , except that there is no rotor at the root @xmath24 ; we think of this as a rotor configuration @xmath8 on the graph obtained from @xmath11 by removing the edges leaving @xmath24 . place a chip on the vertex @xmath14 , and
iterate the rotor - routing process starting with @xmath12 .
stop when the chip first reaches @xmath24 ( which it necessarily must ( * ? ? ?
* lemma  3.6 ) ) , and call the resulting state @xmath48 .
the rotor configuration @xmath16 is in fact another spanning tree @xmath49 oriented towards @xmath24 , and we define @xmath50 to be the spanning tree @xmath49 .
see figure  [ f : planar ] .
since @xmath22 is an abelian group , the action of different divisors @xmath51 will commute : @xmath52 .
however , we emphasize that the actions corresponding to _ different _ root vertices @xmath24 and @xmath53 may _ not _ commute , so it is important to distinguish between @xmath54 and @xmath55 .
the action of @xmath22 descends to the sandpile group @xmath0 .
in fact , for any choice of root vertex @xmath24 , the resulting action of @xmath0 on @xmath4 is free and transitive ( * ? ? ?
* lemmas 3.17 and 3.19 ) .
( this property has been vastly generalized to abelian networks by bond and levine in their recent preprint @xcite . )
we can now restate ellenberg s question @xcite as : if @xmath56 is the permutation defined by this action , does @xmath57 for all @xmath58 and @xmath59 ?
in other words , is the rotor - routing @xmath0-torsor structure on the set of spanning trees @xmath4 independent of basepoint ?
theorem  [ theorem : main ] completely answers this question .
a _ path _ @xmath60 in the graph @xmath1 is an edge - injective map @xmath61 from the length-@xmath62 path @xmath63 ; it may or may not be vertex - injective .
a _ cycle _ @xmath64 in @xmath1 is the image of a vertex - injective and edge - injective map @xmath65 from the length-@xmath62 cycle @xmath66 ( a circle subdivided into @xmath62 edges ) .
we consider both paths @xmath60 and cycles @xmath64 to be oriented , but we identify maps @xmath67 that differ by cyclic permutation ( in other words , we do not specify a basepoint for cycles ) .
we denote by @xmath68 or @xmath69 the same path or cycle with the opposite orientation .
given an oriented cycle @xmath64 in a ribbon graph and a vertex @xmath70 , the edges incident to @xmath14 naturally fall into three classes : the two edges involved in @xmath64 , the edges ( if any ) on the left of @xmath64 , and the edges ( if any ) on the right of @xmath64 .
the edges to the right of @xmath64 , for example , are the ones that occur after the out - edge at @xmath14 in @xmath64 and before the in - edge at @xmath14 in @xmath64 , in the cyclic order at @xmath14 .
these three classes are disjoint , since our graphs have no self - loops ; however , if an edge @xmath71 has both endpoints @xmath14 and @xmath9 on @xmath64 , it is possible that @xmath71 is on the left of @xmath64 at @xmath14 but on the right of @xmath64 at @xmath9 .    a cycle @xmath64 is _ nonseparating _
if there exists a path @xmath60 in @xmath1 whose endpoints lie on @xmath64 and is disjoint from @xmath64 otherwise , with its first edge on the left of @xmath64 and its last edge on the right of @xmath64 .
such a path @xmath60 is a _ witness
_ that @xmath64 is nonseparating .
the cycle @xmath64 is _ separating _ if no such witness exists ; this is equivalent to saying that @xmath64 is a separating curve on the surface associated to the ribbon graph @xmath1 .
the ribbon graph @xmath1 is _ planar _ if every cycle is separating ( equivalently , if its associated surface is a sphere ) .    * unicycles and rotor - routing . *
a _ unicycle _ @xmath72 consists of a rotor configuration @xmath8 which contains exactly one directed cycle @xmath73 , together with a vertex @xmath26 lying on the cycle @xmath73 .
it is not difficult to see that applying rotor - routing to the configuration @xmath72 preserves these conditions , and thus takes unicycles to unicycles ( * ? ? ?
* lemma 3.4 ) .
in fact , the following is is a key observation that we will use many times :    [ lemma : periodic ] let @xmath72 be a unicycle on a graph with @xmath74 edges . iterating the rotor routing process @xmath75 times starting at @xmath72
, the chip traverses each edge of @xmath1 exactly once in each direction , each rotor makes exactly one full turn , and the stopping state is @xmath76 .
we write @xmath77 if the configuration @xmath15 can be obtained from @xmath12 by iterating the rotor - routing process some positive number of times . by lemma  [ lemma : periodic ]
, the relation @xmath77 is in fact an equivalence relation on unicycles .
the following lemma shows that each equivalence class is naturally in bijection with the directed edges of @xmath1 .
we write @xmath78 if @xmath77 and @xmath79=e$ ] .
[ lemma : therotorknows ] let @xmath12 be a unicycle . for any vertex @xmath9 and
any directed edge @xmath71 based at @xmath9 , there exists a unique rotor configuration @xmath16 such that @xmath80 .
run the rotor - routing process starting at @xmath12 until the next occurrence of @xmath12 . according to lemma  [ lemma : periodic
] , each rotor makes precisely one full rotation in this process .
since the rotor at @xmath9 only advances when the chip is at @xmath9 , this implies that at some intermediate stage @xmath15 we had @xmath79=e$ ] .
the next rotor - routing step then advances the rotor at @xmath9 , so @xmath16 is unique .
the _ first _ time the chip arrives at @xmath9 when rotor - routing is iterated starting with @xmath12 , the current rotor configuration @xmath16 satisfies @xmath79=\rho[y]$ ] , since @xmath9 has not previously been visited .
accordingly we write @xmath81 as an abbreviation for @xmath82}(\sigma , y)$ ] . finally , when @xmath15 has already been defined , we write @xmath83 as shorthand for the rotor routing process beginning with @xmath12 and ending at the first occurrence of @xmath84 .
[ def : reversible ] given a unicycle @xmath72 with directed cycle @xmath85 , let @xmath86 denote the configuration obtained from @xmath8 by reversing the edges of @xmath64 , and keeping all other rotors unchanged .
we say that @xmath64 is _ reversible _ if @xmath87 .    before using this terminology
, we need to verify that it is well - defined : in other words , that reversibility really is a property of the _ cycle _ @xmath64 , and does not depend on the choice of @xmath8 nor the choice of @xmath26 . note that by definition , @xmath64 is reversible if and only if @xmath69 is reversible .
we will prove in proposition  [ prop : reversibility2 ] that reversibility is well - defined , but we first need the following technical result .
this result will also be used , along with lemma  [ lemma : almostseparating ] , to show that reversible cycles are quite close to being separating .
[ prop : lcrc ] let @xmath64 be a directed cycle . if there exists a unicycle @xmath72 with @xmath85 such that @xmath88 , then the vertices @xmath89 can be partitioned into two sets @xmath90 such that :    1 .   if @xmath89 is adjacent to @xmath91 along an edge lying on the left of @xmath64 at @xmath14 , then @xmath92 .
2 .   if @xmath89 is adjacent to @xmath91 along an edge lying on the right of @xmath64 at @xmath14 , then @xmath93 .
if @xmath89 is adjacent to @xmath94 , then @xmath9 and @xmath95 are either both in @xmath96 or both in @xmath97 .
any separating curve @xmath64 separates the surface into two parts @xmath96 and @xmath97 as in this proposition .
the key utility of proposition  [ prop : lcrc ] is that it allows us to define the sets @xmath96 and @xmath97 knowing only that @xmath64 is _ reversible _ , without knowing whether or not @xmath64 is separating .    for any vertex @xmath9 in @xmath1 , let @xmath98 be its degree .
consider the rotor - routing process that takes @xmath72 to @xmath99 .
for any vertex @xmath9 _ not _ lying on @xmath64 , the rotor ends at its initial position @xmath10={\overline}{\rho}[y]$ ] , so by lemma  [ lemma : periodic ] , the vertex @xmath9 was visited either 0 or @xmath98 times .
we define @xmath100 and @xmath101 . since @xmath1 has no isolated vertices ( by our convention , all graphs are connected ) , these two sets are disjoint .
since no directed edge is traversed more than once , in the latter case @xmath93 the chip must arrive at @xmath9 along each of the @xmath98 adjacent edges exactly once and leave @xmath9 along each edge exactly once .
now consider a vertex @xmath14 lying on @xmath64 .
since the cycle @xmath64 has been reversed , the rotor at @xmath14 begins at @xmath18=c[x]$ ] and ends at @xmath102={\overline}{c}[x]$ ] .
therefore the rotor rotates @xmath103 times , where @xmath104 is 1 plus the number of edges at @xmath14 on the _ right _ of @xmath64 , and so the chip visits @xmath14 precisely @xmath103 times .
the first @xmath105 times the chip leaves @xmath14 , it leaves along the @xmath105 edges lying on the right of @xmath64 ; the final time it leaves @xmath14 , it leaves along the edge @xmath102={\overline}{c}[x]$ ] of the reversed cycle @xmath64 ( since this is the final state of the rotor ) .
if @xmath89 is adjacent to @xmath91 along an edge lying on the right of @xmath64 at @xmath14 , the chip will leave @xmath14 along this edge at some point , so @xmath93 . similarly , if @xmath89 is adjacent to @xmath91 along an edge lying on the left of @xmath64 at @xmath14 , the chip does not leave @xmath14 along this edge ; in other words , the chip does not arrive at @xmath9 along this edge .
this implies that @xmath92 , because we noted above a vertex @xmath9 lying in @xmath97 must be visited along every one of its incoming edges , and this would include the edge from @xmath14 . finally , if @xmath89 is adjacent to @xmath94 and the chip traverses this edge in either direction , then both @xmath9 and @xmath95 lie in @xmath97 ; if not , both lie in @xmath96 .
the three conditions above uniquely characterize @xmath96 and @xmath97 as those vertices that can be connected by a path to the left side or right side of @xmath64 respectively ; since this description does not depend on @xmath8 or @xmath26 , the sets @xmath96 and @xmath97 depend only on the cycle @xmath64 .
this seems at first to prove that @xmath64 is separating .
however , it does not rule out the possibility of a witness consisting of a single edge joining the left side of @xmath91 to the right side of @xmath106 .
in fact , nonseparating cycles of this form can indeed be reversible , as we will see in the proof of proposition  [ prop : planar_is_reversible ] .
[ prop : reversibility2 ] let @xmath64 be a directed cycle . if @xmath107 for _ some _ unicycle @xmath108 with @xmath109 , then @xmath110 for _ any _ unicycle @xmath76 with @xmath85 . in other words ,
the reversibility of a cycle @xmath64 , as in definition  [ def : reversible ] , is well - defined : it only depends on @xmath64 , and not on a choice of unicycle .    given a unicycle @xmath72 with @xmath85 ,
define its _ maximal reversal _
@xmath111 to be the first state @xmath111 encountered in the rotor - routing process such that @xmath112 and @xmath113={\overline}{c}[u]$ ] . in other words ,
the chip is on the cycle @xmath64 at the vertex @xmath114 , and the rotor at @xmath114 is about to rotate _
past _ the reversed cycle @xmath69 to the edges on the left of @xmath64 for the first time .
we will prove that as long as @xmath107 for _ some _ @xmath115 and some @xmath116 , the maximal reversal of @xmath72 is @xmath99 .
then we will have @xmath117 as desired .
( we can not use lemma  [ lemma : therotorknows ] to deduce this from @xmath113={\overline}{c}[u]={\overline}{\rho}[u]$ ] , because we do not yet know that @xmath118 . )
assume that @xmath107 for some @xmath115 with @xmath109 , and let @xmath96 and @xmath97 be the sets defined by proposition  [ prop : lcrc ] .
the set of directed edges traversed in the process @xmath119 was described in the proof of that proposition : the chip traverses ( 1 ) every directed edge leaving @xmath93 , ( 2 ) each directed edge leaving @xmath91 on the right side of @xmath64 , and ( 3 ) each edge @xmath120 $ ] leaving @xmath91 .
this set of directed edges only depends on @xmath64 , not on @xmath115 or @xmath121 , so we will denote it by @xmath122 .
for any vertex @xmath123 , let @xmath124 denote the number of edges in @xmath122 leaving @xmath123 .
( specifically , @xmath125 for @xmath91 , @xmath126 for @xmath93 , and @xmath127 for @xmath92 . )
since the chip ends up back at @xmath26 , it leaves each vertex @xmath123 the same number of times that it arrives there , so @xmath124 is also the number of directed edges in @xmath122 _ arriving _ at @xmath123 .
now consider an arbitrary rotor configuration @xmath72 with @xmath85 and @xmath128 .
the key to this lemma is the observation that if @xmath111 is the maximal reversal of @xmath72 , the process @xmath129 only traverses edges in @xmath122 .
assume otherwise , and let @xmath71 be the first directed edge _ not _ in @xmath122 traversed .
since the chip begins on @xmath64 , and no edges connect @xmath97 to @xmath96 , this first edge @xmath71 is based at some @xmath91 . by definition of @xmath122 ,
either @xmath71 lies on the left side of @xmath64 , or @xmath130 $ ] itself . in either case , the rotor at @xmath14 must rotate past @xmath120 $ ] before it can reach the edge @xmath71 .
but this is a contradiction , since by definition @xmath111 is the first time the rotor at @xmath91 is @xmath131 $ ] while the chip is there .    in the process @xmath132
, the rotor at @xmath112 rotates from @xmath133=c[u]$ ] to @xmath113={\overline}{c}[u]$ ] , so the chip leaves @xmath114 precisely @xmath134 times .
however , if @xmath135 , the chip must have arrived at @xmath114 precisely @xmath136 times , since it did not originate there . since the total number of edges in @xmath122 arriving at @xmath114 is only @xmath134 , this is a contradiction .
this shows that @xmath137 , and moreover that the chip has arrived at @xmath26 along each of the @xmath138 edges in @xmath122 directed towards @xmath26 .
one of these is the edge @xmath139 $ ] , where @xmath140 is the next vertex after @xmath26 in the cycle @xmath64 , so we must have @xmath141={\overline}{c}[w]$ ] as well .
but if the rotor at @xmath123 has reversed from @xmath142=c[w]$ ] to @xmath141={\overline}{c}[w]$ ] , the chip must have left @xmath124 times .
to arrive @xmath124 times at @xmath123 , the rotor at the _ next _ vertex after @xmath123 must have been reversed , and so on . by induction , we conclude that we have @xmath17={\overline}{c}[x]={\overline}{\rho}[x]$ ] for all @xmath91 .
in particular , the directed cycle of @xmath16 is @xmath143 .
since the process @xmath144 does not visit @xmath96 , we know that @xmath79=\rho[y]={\overline}{\rho}[y]$ ] for all @xmath92 .
it remains to show that @xmath79=\rho[y]$ ] for @xmath93 . to do this ,
let @xmath145 be the maximal reversal of @xmath146 .
since @xmath143 , the process @xmath147 does not visit @xmath148 , so @xmath149=\sigma[y]$ ] for all @xmath93 .
but the previous paragraph shows that @xmath150=\overline{\sigma}[v]=\rho[v]$ ] . since @xmath151
, lemma  [ lemma : therotorknows ] implies that @xmath152 .
therefore @xmath79=\rho[y]$ ] for all @xmath89 , and @xmath17={\overline}{\rho}[x]$ ] for all @xmath91 , demonstrating that @xmath153 as desired .
[ lemma : almostseparating ] let @xmath64 be a directed cycle .
if @xmath64 is reversible , then any path @xmath60 witnessing that @xmath64 is nonseparating must have length 1 .
conversely , if @xmath64 is a separating cycle , then @xmath64 is reversible
.    let @xmath60 be a path of length @xmath154 , with vertices @xmath155 , with endpoints @xmath156 and @xmath157 .
if @xmath60 is a witness that @xmath64 is nonseparating , the second vertex @xmath158 lies in @xmath96 , and the second - to - last vertex @xmath159 lies in @xmath97 . but this contradicts condition ( 3 ) of proposition  [ prop : lcrc ] , which says that membership in @xmath96 is locally constant along paths in @xmath160 .
therefore any witness @xmath60 must have length 1 .    to prove that any separating cycle is reversible
, we will go through the proof of proposition  [ prop : lcrc ] in reverse ; the resulting argument closely parallels ( * ? ? ?
* corollary 4.11 ) , where this was proved for planar ribbon graphs .
( in fact , a special case of reversibility , for cycles in an @xmath161 planar grid graph , was proved even earlier , in ( * ? ? ?
* proposition i ) . )
if @xmath64 is a separating cycle , we can partition the vertices @xmath89 into @xmath162 , where @xmath163 ( resp .
@xmath164 ) if there exists a path from @xmath9 to @xmath64 ending with an edge on the left side ( resp .
right side ) of @xmath64 .
every @xmath9 can be connected to @xmath64 by some path , since @xmath1 is connected , so @xmath165 . if @xmath166 could be connected to both sides of @xmath64 , splicing these paths together would give a witness that @xmath64 was nonseparating , so @xmath167 .
let @xmath168 be the ribbon graph obtained from @xmath1 by deleting all vertices in @xmath169 , all edges adjacent to @xmath169 , and all edges lying on the right of @xmath64 .
because @xmath64 is separating , no edge on the left side of @xmath64 is removed ( every such edge either joins @xmath91 to @xmath163 , or joins @xmath91 to @xmath106 , in which case it lies on the left side of @xmath64 at both endpoints ) . since no rotor in @xmath170 points along a removed edge
, we can restrict the rotor configuration @xmath8 to @xmath168 as @xmath171 .
it is easy to run the rotor routing process starting at @xmath172 , because at each @xmath91 the two edges of @xmath64 are now adjacent in the cyclic order .
thus at each step , the rotor at one vertex on @xmath64 is reversed , and the chip moves to the previous vertex on @xmath64 . after a number of steps equal to the length of @xmath64
, each rotor on @xmath64 has moved once , the cycle @xmath64 has been reversed , and the chip is back at @xmath26 ; in other words , the process takes @xmath173 .
since rotor - routing on @xmath168 is periodic , continuing gives a rotor - routing process on @xmath168 taking @xmath174 to @xmath172 .
but this segment of the process does not involve any of the deleted edges .
therefore if the rotor - routing process is run in parallel on @xmath1 starting with @xmath99 , the two processes will take precisely the same series of steps .
when the former terminates at @xmath175 , the latter will therefore be at @xmath76 , demonstrating that @xmath176 as desired .
[ prop : planar_is_reversible ] a connected ribbon graph @xmath1 without loops is planar if and only if all cycles on @xmath1 are reversible .
if @xmath1 is planar , the jordan curve theorem implies that every cycle is separating , so lemma  [ lemma : almostseparating ] implies that every cycle is reversible .
for the converse , suppose @xmath1 is not planar , so it has at least one nonseparating cycle , but that every cycle on @xmath1 is reversible .
our first claim is that every nonseparating cycle on @xmath1 has length 2 .
indeed , let @xmath60 be a path witnessing that the cycle @xmath64 is nonseparating , joining the left side of @xmath91 to the right side of @xmath106 .
this path splits @xmath64 into two directed arcs from @xmath14 to @xmath177 . call them @xmath178 and @xmath179 , labeled so that @xmath180 .
let @xmath181 be the cycle @xmath182 , and let @xmath183 be the cycle @xmath184 .
the arc @xmath178 is now a witness that @xmath181 is nonseparating , since @xmath178 lies on the left of @xmath181 at @xmath14 but on the right of @xmath181 at @xmath177 .
similarly @xmath179 witnesses that @xmath183 is nonseparating , since @xmath179 lies on the right of @xmath183 at @xmath14 but on the left of @xmath183 at @xmath177 . since all three cycles @xmath64 , @xmath181 , and @xmath183 are reversible by assumption , lemma  [ lemma : almostseparating ] implies that the paths @xmath60 , @xmath179 and @xmath178 each have length 1 .
in particular , the original cycle @xmath185 has length 2 , as claimed .
let us rename the edge @xmath60 to @xmath186 , and the cycle @xmath64 to @xmath187 , to expose the cyclic symmetry of our notation .
then the cyclic order of these edges at @xmath14 and @xmath177 is the same , namely @xmath188 and @xmath189 respectively , and @xmath190 with indices modulo 3 .
see figure  [ fig : a1a2a3 ] .
let us say that an edge at @xmath14 _ lies between _ @xmath191 and @xmath192 , for @xmath193 , @xmath194 , or @xmath195 , if it comes after @xmath191 and before @xmath192 in the cyclic order of edges at @xmath14 .
similarly , we will say that an edge at @xmath177 lies between @xmath196 and @xmath197 if it comes after @xmath196 and before @xmath197 in the cyclic ordering of edges at @xmath177 .     and
@xmath186 in the proof of proposition  [ prop : planar_is_reversible ] .
the orientation at each vertex is clockwise with respect to the page.,width=144 ]    before moving on , we observe that _ any _ path from @xmath14 to @xmath177 not passing through @xmath14 or @xmath177 except at its endpoints will witness that one of the cycles @xmath198 is nonseparating , so by lemma  [ lemma : almostseparating ] any such path has length 1 . indeed ,
each edge at @xmath14 other than the @xmath191 lies on the right side of exactly one cycle @xmath198 , and on the left side of the other two ; similarly , each edge at @xmath177 lies on the left of exactly one @xmath198 and on the right of the other two .
thus the pigeonhole principle implies that no path from @xmath14 to @xmath177 can lie on the same side of @xmath198 for all three simultaneously .
we may also assume that no edge @xmath71 connects @xmath14 to @xmath177 and lies between @xmath179 and @xmath186 at @xmath14 , after repeating the following reduction step .
suppose @xmath71 is such an edge .
if @xmath71 lies between @xmath199 and @xmath200 at @xmath177 , then replace @xmath179 with @xmath71 .
if @xmath71 lies between @xmath199 and @xmath201 at @xmath177 , then replace @xmath186 with @xmath71 .
if @xmath71 lies between @xmath200 and @xmath201 at @xmath177 , then replace either @xmath179 or @xmath186 with @xmath71 .
in each case , this reduces the number of edges at @xmath14 between @xmath179 and @xmath186 while preserving the cyclic orderings of @xmath188 and @xmath202 at @xmath14 and @xmath177 , so we may repeat this step until there are no edges to @xmath177 lying between @xmath179 and @xmath186 at @xmath14 .
let @xmath12 be any unicycle with @xmath203 ( e.g.  by adding @xmath178 to a spanning tree that uses @xmath179 ; see  [ sec : spanningtrees ] for details ) .
we have @xmath18=a_1 $ ] and @xmath204={\overline}{a}_2 $ ] by definition .
define @xmath16 by @xmath205 .
we had @xmath102=a_2 $ ] and @xmath206={\overline}{a}_1 $ ] , and the latter implies @xmath207={\overline}{a}_1 $ ] .
but we know moreover that @xmath17=a_3 $ ] .
this is because in the rotor routing process @xmath208 , starting from the last time the chip leaves @xmath14 and ending at its arrival at @xmath177 , the chip traces out a path from @xmath14 to @xmath177 .
that path has length 1 by the argument above , which means that the chip arrives for the first time at @xmath177 along some edge from @xmath14 .
but we have assumed that @xmath186 is the first such edge that the rotor at @xmath14 encounters when it starts at @xmath179 .
finally , define @xmath209 by @xmath210 .
we had @xmath211=a_1 $ ] and @xmath212={\overline}{a_3}$ ] , and the former implies @xmath213=a_1 $ ] .
we have assumed that all cycles are reversible , so @xmath214 and @xmath215 . by transitivity
this implies @xmath216 ; since @xmath213=\rho[x]=a_1 $ ] , it must be that @xmath152 .
however , this would mean that the rotor @xmath217 $ ] ends up at @xmath218 , which is impossible : this rotor began at @xmath212={\overline}{a}_3 $ ] , so the first edge to @xmath14 that it encounters can not be @xmath218 ( it would encounter @xmath199 first , if not some other edge ) .
this contradiction completes the proof that not all cycles in @xmath1 can be reversible if @xmath1 has nonseparating cycles .
* spanning trees and unicycles .
* the reason that we studied unicycles so carefully in the previous section is that they are closely related to spanning trees , as we now explain .
let @xmath42 be a spanning tree .
any two vertices @xmath14 and @xmath9 are connected by a unique geodesic path in @xmath41 ; we write @xmath219 for that path , oriented from @xmath14 to @xmath9 . therefore if @xmath71 is an edge of @xmath1 connecting @xmath9 to @xmath14 , the undirected graph @xmath220 contains a unique cycle .
now given a directed edge @xmath71 from @xmath9 to @xmath14 and a spanning tree @xmath41 , define a rotor configuration @xmath221 as follows .
let @xmath222 denote the collection of directed edges obtained from @xmath41 by directing each edge towards the `` root '' @xmath9 .
then the collection of directed edges @xmath223 forms a rotor configuration @xmath221 with unique directed cycle @xmath224 .
we will denote this directed cycle by @xmath225 .
if @xmath41 already contains @xmath71 , then @xmath225 is the length 2 cycle consisting only of @xmath226 ; note that this cycle is trivially reversible in any ribbon graph .
if we took instead the edge @xmath227 from @xmath14 to @xmath9 , the rotor configuration @xmath228 comes from the directed graph @xmath229 , with directed cycle @xmath230
. the following observation will be fundamental for us .      for @xmath232
this follows because @xmath233 for vertices @xmath234 , we have @xmath235=t_y[z]$ ] , since the geodesics @xmath236 and @xmath237 share an initial segment . then @xmath238=\rho_{{\overline}e}(t)[z]$ ] , and @xmath238 = { \overline}{\rho_e(t)}[z]$ ] since @xmath95 is not on the cycle .    given a spanning tree @xmath239 and an edge @xmath71 from @xmath9 to @xmath14 , the rotor configuration @xmath240 is a unicycle since @xmath241 .
conversely , given a unicycle @xmath72 , deleting the directed edge @xmath242 $ ] yields a spanning tree @xmath243 with all edges directed towards the root @xmath26 .
this correspondence between spanning trees and rotor configurations lets us partially model the action of the sandpile group on spanning trees .
[ prop : rotorscompute ] given a spanning tree @xmath41 in a ribbon graph @xmath1 , let @xmath71 be an edge not in @xmath41 , and @xmath53 any vertex on the @xmath41-geodesic between the endpoints @xmath244 of @xmath71 ( including @xmath14 or @xmath9 ) .
let @xmath245 ( so that @xmath246 ) .
then @xmath247 is a unicycle , and if @xmath16 is the rotor configuration defined by @xmath248 , then @xmath249    the rotor configuration @xmath8 is equal to @xmath223 , where @xmath222 is obtained by directing each edge of @xmath41 towards @xmath9 .
the configuration @xmath8 then has a unique cycle , namely @xmath250 , and by assumption @xmath53 is on that cycle , so @xmath247 is a unicycle .
next , by definition , @xmath251 is obtained by beginning with @xmath222 , placing the chip at @xmath53 , and running the rotor - routing process until the chip arrives at @xmath9 for the first time .
the presence of the rotor @xmath10=e$ ] does not affect this process , since @xmath9 is not reached until the last step .
therefore the rotor - routing process @xmath252 takes precisely the same steps as the process @xmath253 .
we conclude that @xmath254 . since @xmath255 by definition , this proves the lemma .    * basepoint - independence and reversible cycles . * before giving our last few lemmas and proving our main theorem , let us give some intuition as to why reversibility should be at all related to basepoint - independence . consider a spanning tree @xmath256 , and let @xmath14 and @xmath9 be adjacent vertices connected by an edge @xmath71 from @xmath14 to @xmath9 .
suppose that we wish to compute @xmath257 for some divisors @xmath51 .
the first step is to compute @xmath258 . by proposition  [ prop : rotorscompute ] , the computation of @xmath259 can be modeled ( at least for some divisors ) by a rotor - routing process ; this process terminates with the rotor configuration @xmath260 , which is the spanning tree @xmath41 oriented towards the basepoint @xmath14 .
the next step is to compute @xmath261 , and to do this we first need to re - orient the tree @xmath41 towards the basepoint @xmath9 .
this does not affect any edges except those on the @xmath41-geodesic between @xmath14 and @xmath9 , which are reversed . if we add the edge @xmath71 to our rotor configuration @xmath260 to obtain @xmath262 , the reversal of the @xmath41-geodesic @xmath219 amounts to the reversal of the cycle @xmath224 in the rotor configuration @xmath221 .
if this cycle is reversible , then we can compute both steps @xmath263 and @xmath264 as part of a _ single rotor - routing process _
, even though these actions involve two _ different _ basepoints .
this basic observation is the heart of the relation between reversibility and basepoint - independence , upon which our main theorem rests .
the next lemma capitalizes on this observation .
[ lemma : ongeodesic ] given a spanning tree @xmath41 in a ribbon graph @xmath1 and an edge @xmath71 from @xmath9 to @xmath14 , let @xmath53 be any vertex on the geodesic @xmath219 ( including @xmath265 or @xmath266 ) , and set @xmath267 .
let @xmath245 ; as in proposition  [ prop : rotorscompute ] , @xmath247 is a unicycle .
define rotor configurations @xmath270 , @xmath16 , and @xmath209 by @xmath271 note that @xmath272={\overline}{\rho}[x]={\overline}{e}$ ] , and similarly @xmath79=\rho[y]=e$ ] implies that @xmath213={\overline}{\sigma}[x]={\overline}{e}$ ] .
in fact , by lemma  [ lemma : ovrho ] and proposition  [ prop : rotorscompute ] , we have      thus we see that equality is exactly the condition @xmath276 .
moreover , since @xmath272=\tau[x]={\overline}{e}$ ] , lemma  [ lemma : therotorknows ] implies that this is equivalent to @xmath277 .
consider the following diagram relating these six rotor configurations : @xmath278\ar_{\txt{\large ? } } @{2{<}~2{>}}[dl]&&(\sigma , y)\ar^{\txt{\large ? } } @{2{<}~2{>}}[dr]\\    ( { \overline}{\rho},s)\ar@{2{<}~2{>}}[dr]&&&&({\overline}{\sigma},y)\ar@{2{<}~2{>}}[dl]\\    & ( \psi , x)\ar_{\txt{\large ? } } @{2{<}~2{>}}[rr]&&(\tau , x ) } \ ] ] three of these edges hold by definition , namely @xmath279 , @xmath280 , and @xmath281 . by the results of section  [ sec : reversibility ] , the upper left edge @xmath282 is equivalent to the cycle @xmath283 being reversible .
similarly the upper right edge is equivalent to @xmath284 being reversible .
finally , we have already shown that the bottom edge is equivalent to the assertion
. therefore parts ( a ) and ( b ) of the lemma each assert that if two of the remaining three edges hold , the third does as well . but
this follows immediately from the transitivity of this equivalence relation .      by lemma  [ lemma : almostseparating ] ,
the hypothesis of lemma  [ lemma : ongeodesic](a ) is always satisfied if @xmath1 is planar , so @xmath286 now applying @xmath287 to both sides , we conclude that @xmath288 .
we emphasize that we are _ not _ yet claiming that @xmath289 , since the condition of corollary  [ corollary : first ] that @xmath53 lies on the @xmath41-geodesic from @xmath14 to @xmath9 need not hold for all trees @xmath41 .
[ lemma : rununtil ] let @xmath290 be a spanning tree of an arbitrary ribbon graph @xmath1 .
let @xmath14 , @xmath95 , and @xmath53 be arbitrary vertices such that @xmath53 is the first vertex on @xmath291 visited by the rotor - routing process starting with @xmath292
. then @xmath293    define @xmath294 by @xmath295 ; by definition , the underlying spanning tree @xmath49 of @xmath294 is equal to @xmath296 .
similarly , define @xmath297 by @xmath298 , with underlying spanning tree @xmath299 .    as we noted in the proof of lemma  [ lemma : ovrho ] , the configurations @xmath260 and
@xmath300 coincide for all vertices @xmath301 , while on this geodesic , the difference between @xmath260 and @xmath300 is that @xmath291 is reversed and thus replaced by @xmath302 .
therefore if the rotor - routing processes are run in parallel starting with @xmath292 and @xmath303 , the same steps will be taken until the chip first arrives at @xmath291 . by our assumption on @xmath53 ,
this occurs when the chip reaches @xmath53 , i.e.  when the second process reaches @xmath304 .
the first process coincides except that the geodesic between @xmath14 and @xmath53 is directed towards @xmath14 instead .
but by the first sentence of this paragraph , this is the directed tree @xmath305 obtained from @xmath49 by directing its edges towards @xmath14 . to sum up
, we have observed that @xmath306 .    by definition
, @xmath307 can be computed by starting with @xmath308 and running the rotor - routing process until the chip reaches @xmath14 .
the above shows that the first process reaches @xmath308 as an intermediate stage ; therefore continuing this process eventually yields @xmath309 , since @xmath298 by definition .
this shows that @xmath310 . expanding out our notation
, this says that @xmath311 applying @xmath312 to both sides yields @xmath313 , as desired .
let @xmath53 be the first vertex on @xmath219 that is reached by the rotor - routing process starting with @xmath292 .
note that @xmath53 is _ a fortiori _
the first vertex visited on @xmath302 or on @xmath315 as well .
we use this to apply lemma  [ lemma : rununtil ] twice , giving @xmath316 write @xmath49 for the tree @xmath296 appearing in .
notice that by definition of @xmath53 , in the rotor routing process that computes @xmath296 , the chip never touches any vertex on @xmath219 until it reaches @xmath53 , so @xmath317 .
so @xmath53 lies on @xmath318 by corollary  [ corollary : first ] applied to @xmath49 ( with the roles of @xmath14 and @xmath9 reversed ) , we have @xmath319 the desired result follows once this equation is combined with , by applying @xmath320 to the left - hand side of and applying @xmath321 to the right - hand side of .
( @xmath322 ) first , assume that @xmath1 is a planar ribbon graph .
our goal is to prove for any @xmath43 that @xmath323 as elements of @xmath324 for arbitrary vertices @xmath14 and @xmath9 .
since @xmath1 is connected , it suffices to prove this when @xmath14 and @xmath9 are adjacent , with the general case following by induction .
fix adjacent vertices @xmath14 and @xmath9 .
the group @xmath22 is generated by divisors of the form @xmath325 as @xmath95 ranges over the vertices of @xmath1 .
but we proved in corollary  [ corollary : fulladjacent ] that @xmath326 when @xmath14 and @xmath9 are adjacent .
it follows that @xmath323 for any @xmath43 .
this completes the proof that if @xmath1 is a planar ribbon graph , the action of @xmath0 on @xmath7 is independent of the basepoint .
( @xmath327 ) conversely , assume that the action of @xmath0 on @xmath7 is independent of the basepoint . choose any edge @xmath71 of @xmath1 with endpoints @xmath14 and @xmath9 . we claim that for any @xmath239 and any @xmath43 , the cycle @xmath225 is reversible if and only if the cycle @xmath328 is reversible .
( we have dropped the subscripted basepoints , since by assumption the action of @xmath0 does not depend on them . )
we first prove the claim in the case when @xmath329 for some @xmath95 .    as in the proof of corollary
[ corollary : fulladjacent ] , let @xmath53 be the first vertex on @xmath219 reached by the rotor - routing process starting with @xmath292 . since we have assumed independence of basepoint , this process @xmath330 computes @xmath331 .
since this process does not reach @xmath219 until the last step , the cycle @xmath225 remains unchanged , i.e. @xmath332 .
we may now apply lemma  [ lemma : ongeodesic ] to @xmath333 , since @xmath53 lies on the cycle @xmath332 .
the condition is always satisfied by our assumption that the action is independent of the basepoint , so applying lemma  [ lemma : ongeodesic](b ) implies that @xmath334 completing the proof of the claim for @xmath329 .
since @xmath22 is generated by such elements , the claim follows .
we now show that every cycle @xmath64 in @xmath1 is reversible . given a cycle @xmath64 ,
let @xmath71 be any edge of @xmath64 .
extend the edge @xmath71 to a spanning tree @xmath41 , and separately extend the path @xmath335 to a spanning tree @xmath49 .
since @xmath0 acts transitively on @xmath7 , there exists @xmath43 for which @xmath336 .
thus the claim above shows that @xmath225 is reversible if and only if @xmath268 is reversible . but
@xmath225 is the trivial cycle @xmath226 , which is always reversible in any ribbon graph , while @xmath268 is our original cycle @xmath64 .
we conclude that every cycle in @xmath1 is reversible . by proposition
[ prop : planar_is_reversible ] , @xmath1 is a planar ribbon graph .
* acknowledgements .
* we are very grateful to jordan ellenberg for calling our attention to the question considered in this paper , and we thank math overflow for providing a venue for the question . the first and third authors would also like to thank the second author for suggesting that the three of us get together to work on this problem , which was great fun !
a.  holroyd , l.  levine , k.  mszros , y.  peres , j.  propp and d.  wilson , chip - firing and rotor - routing on directed graphs , in _ in and out of equilibrium 2 _ , eds . v.  sidoravicius and m.  e.  vares , _ progress in probability _ 60 , birkhuser ( 2008 ) , 331364 . `
arxiv:0801.3306 ` .


spatially periodic non - centrosymmetric systems are able to transport nonequilibrium particles in the absence of an average macroscopic force resulting in the ratchet effect  @xcite .
ratchet effects whose prerequisites are simultaneous breaking of both thermal equilibrium and spatial inversion symmetry can be realized in a great variety of forms and in particular , as electric transport in semiconductor systems  @xcite .
recent experiments demonstrated that terahertz radiation induced ratchet effects , can be efficiently excited in semiconductor quantum wells ( qw )  @xcite and graphene  @xcite with lateral superlattice structures .
these experiments allowed one to explore basic physics of the ratchet effects in low dimensional electron systems , provide information on non - equilibrium transport in such systems and demonstrate that ratchet effects can be applied for room temperature terahertz radiation detection  @xcite .
the latter , besides high sensitivity and short response times , offer new functionality being a good candidate for all - electric detection of the radiation polarization state including radiation helicity  @xcite .    here
we report on the observation and study of magnetic quantum ratchet effect in ( cd , mn)te/(cd , mg)te diluted magnetic heterostructures and cdte / cdmgte qws superimposed with lateral asymmetric superlattices . applying magnetic field @xmath2 along the growth direction
we observe that the ratchet current exhibits sign - alternating @xmath1-periodic oscillations with amplitudes immensely larger than the ratchet signal at zero magnetic field .
the results are analyzed in terms of the theory of magnetic ratchet effects in qw structures with a lateral asymmetric periodic potential  @xcite .
we show that the photocurrent generation is based on the combined action of a spatially periodic in - plane potential and the spatially modulated light due to the near field effects of radiation diffraction .
corresponding theoretical analysis describes the experimental results .
the paper is organized as follows .
first , we present our samples and results of magneto - transport characterization ( sec .
[ sample ] ) , as well as briefly describe the experimental technique ( sec .
[ experi ] ) . in sec .
[ erstes ] we discuss the results on the photocurrents generated at zero magnetic field .
section  [ magnetic ] contains the main experimental results on the magnetic quantum ratchet effect .
it has three subsections describing the basic results , the laser beam scan across the structure and the effect of the electrostatic potential acting on the electron gas . in the following sec .  [ discussion ]
we present the theory , calculate dependencies of the photocurrent on magnetic field , and compare them with the experimental data . finally , in sec .
[ summary ] , we summarize the results and discuss prospectives of future experimental studies of the magnetic quantum ratchet effect .
experiments are carried out on ( cd , mn)te / cdmgte and cdte / cdmgte single qw structures grown by molecular beam epitaxy on ( 001)-oriented gaas substrates  @xcite . during the growth the fluxes of cd , te , mg and mn have been supplied from elemental sources while iodine flux has been obtained from a zni@xmath3 source .
the schematic layout of the ( cd , mn)te / cdmgte layer structure and sketch of the qw are presented in fig .
[ structure ] .
the thick buffer ( @xmath4  6  @xmath5 m ) , consisting of cdte and cd@xmath6mg@xmath7te layers and cdte / cd@xmath6mg@xmath7te short period superlattices , have been grown to reduce the number of dislocations resulting from the strong lattice mismatch between gaas and cd@xmath8mg@xmath9te .
the well width in both kinds of structures is 9.7  nm and the cd@xmath6mg@xmath7te alloy serves as a barrier material .
the composition of barriers has been determined from photoluminescence spectra . in order to obtain a two - dimensional electron gas
the structures have been modulation doped by iodine donors incorporated into the 5  nm thick region of the top barrier at the distance of either 10 or 15  nm away from the qw , in ( cd , mn)te / cdmgte and cdte / cdmgte structures , respectively .
doped regions have been overgrown by an undoped cap layer with thickness of either 50  nm , for ( cd , mn)te / cdmgte qws , or 75  nm , for cdte / cdmgte qws . the density of the two - dimensional electron gas , @xmath10 , and the electron mobility , @xmath11 , determined at liquid helium temperature ( 4.2 k ) for ( cd , mn)te / cdmgte structures are @xmath12 @xmath13 and @xmath14  @xmath15/vs , and for cdte / cdmgte structures @xmath16  @xmath13 and @xmath17  @xmath18/vs .
( cd , mn)te / cdmgte structures contain a single qw made of ( cd , mn)te digital magnetic alloy  @xcite , while cdte / cdmgte structures contain a single qw made of non - magnetic cdte . in the former structures two evenly spaced cd@xmath19mn@xmath20te thin layers
were inserted during the qw growth , see fig .
[ structure](a ) .
these two layers were three - monolayer - thick ( ml ) and separated from each other and from the barriers by 8 ml thick cdte ( 1 ml is 0.324 nm thick ) .
incorporation of mn@xmath21 ions carrying localized spin @xmath22 into the qw region leads to a strong enhancement of the effective @xmath23 factor of band carriers , and hence to an enhanced zeeman splitting .
this has been shown in previous magneto - photoluminescence ( pl ) studies performed on the ( cd , mn)te / cdmgte samples made from the same wafer as used in the current study .
magneto - pl showed strong red - shift of the pl line with increasing magnetic field  @xcite . from fitting of the modified brillouin function  @xcite to the field dependence of the pl line position the effective average concentration of mn in the digital alloy has been estimated to be @xmath24 = 0.015 .
the samples have also been characterized by electrical transport measurements . at low temperatures pronounced shubnikov - de - haas ( sdh ) oscillations and
very well resolved quantum hall plateaus have been observed .
two characteristic dependencies measured in hall bar and van - der - pauw geometries in ( cd , mn)te structures are shown in figs .  [ transport](a ) and  [ transport](b ) , respectively . while the amplitude of the oscillations increases in cdte and ( cd , mn)te structures at rather high temperatures with raising magnetic field ,
the amplitude of the oscillations in ( cd , mn)te at low temperatures is a more complex function of the magnetic field @xmath25 and shows a beating like pattern .
the latter is clearly seen in fig .  [ transport](a ) for @xmath26 k ( @xmath27  t ) and in panel ( b ) for 4.2  k ( @xmath28  t ) .
such a behavior in diluted magnetic semiconductor ( dms ) structures is well known to be caused by the exchange interaction of electrons with mn@xmath21 ions resulting in the exchange - enhanced zeeman splitting  @xcite .
mn@xmath20te layers inserted during qw growth .
( b ) cross - section of the dual - grating gate superlattice formed by metal fingers deposited on top of the qw structure .
structure composition is given on the right side of the sketch .
the supercell of the grating gate fingers consists of metal stripes having two different widths @xmath29 and @xmath30 separated by spacings @xmath31 and @xmath32 .
this asymmetric supercell is repeated @xmath33 times to create a superlattice with period @xmath34 , see table  [ tab1 ] .
( c ) photograph of the sample together with schematic top view of the dual - grating gate superlattice .
thin red and thick blue lines sketch the top gates having different thicknesses and spacings . ]      for the ratchet effect experiments we fabricated a dual grating top gate ( dgg ) superlattices on top of the qw structures .
a sketch of the gate fingers and a corresponding optical micrograph are shown in figs .  [ structure](b ) and  [ structure](c ) , respectively .
the grating - gate supercell consists of two metal stripes having widths @xmath29 and @xmath30 and spacings @xmath31 and @xmath32 .
the widths of the thin stripes were either @xmath35 or 1.85  @xmath11 m , that of the wide stripes @xmath36  @xmath11 m , see tab .
spacings in most of the structures were @xmath37  @xmath11 m and @xmath38  @xmath11 m . to clarify the influence of the superlattice parameters on the ratchet effect , we additionally fabricated one dgg structure with substantially different spacings @xmath39  @xmath11 m and @xmath40  @xmath11 m between the stripes . the supercell is repeated @xmath33 times to produce an asymmetric superlattice with the period @xmath41 , see refs .
the two subgrating gates , each formed by stripes of identical width , can be biased independently . for this purpose
all thin ( top gate g1 ) and thick ( top gate g2 ) grating stripes have been connected by additional gold stripes , see thick horizontal lines in fig .
[ structure](c ) .
the benefit of this geometry is that the periodic lateral electrostatic potential of the dgg structure can be varied in a controllable way .
several @xmath42 and @xmath43mm@xmath44 size samples of the same batch were prepared .
the dgg structure has been fabricated on one half of the samples so that the other half of the sample surface remains unpatterned serving as reference part .
it has a total area of @xmath45 where @xmath46 is the overlap length of two gate lattices being of the order of hundreds of micrometer .
.sample parameters shown in detail in figs .
[ structure](b ) and  [ structure](c ) [ cols="^,^,^,^,^,^,^",options="header " , ]     . [ tab1 ]    in all samples , besides the sample # 4 , the gate fingers have been made by electron beam lithography and subsequent deposition of 25 or 30  nm thick gold films .
the dgg structure of sample # 4 has been fabricated by depositing @xmath47 nm thick dy and @xmath48 nm thick au films .
pre - magnetizating the hard magnetic dy - based superlattice enables realization of a magnetic ratchet which features an inhomogeneous magnetic field , recently suggested in  @xcite .    for the photocurrent and magneto - transport measurements
several pairs of ohmic contacts have been prepared , see inset in figs.[transport](b ) and  [ gatedepb0 ] .
contact pads were placed in a way that the photo - induced currents can be measured perpendicularly to the metal fingers ( @xmath49 , contacts 1 and 4 ) and parallel to them ( @xmath50 , contacts 3 and 5 ) .
two additional contacts ( 2 and 6 ) were used for detecting the photocurrent signals from the unpatterned area as a reference .
magneto - transport data for the ( cd , mn)te dgg structure at liquid helium temperature are shown in fig .
[ transport](b ) .
note that the application of gate voltages to individual gates does not visibly influence the period of @xmath1-oscillations of the longitudinal resistance , @xmath51 ( not shown ) .
this is due to the fact that the area of dgg fingers is a very small fraction of the whole sample area .
moreover , the variation of bias voltage applied to gate 1 ( @xmath52 ) also does not substantially affect the value of @xmath51 at a fixed magnetic field .
the only detected deviation from the value of @xmath51 at zero gate voltages is an increase of the resistance by @xmath53  % for @xmath54  v. the gate voltage dependencies of the longitudinal resistance measured at one of the maxima of the @xmath51 oscillations are shown in the inset in fig .
[ transport](b ) .
the dependencies are obtained at fixed magnetic field by sweeping the potential applied to one gate and holding the other one at zero bias .
the observed behavior of @xmath51 will be addressed below in the discussion of photocurrent data .
note that while the overall characteristic of the data obtained at different cooldowns is the same the onset of the resistance increase can be shifted on the gate voltage scale by @xmath55  v. this is ascribed to cooldown dependent charge trapping in the insulator .     and 4.2  k in ( cd , mn)te structures .
( a ) unpatterned hall bar sample .
( b ) square shaped samples with a dgg superlattice .
for this measurement a van - der - pauw geometry has been used .
inset in ( b ) shows the dependence of the longitudinal resistance , @xmath51 , on the gate voltage @xmath56 ( @xmath57 ) measured for @xmath58 ( @xmath59 ) .
measurements are presented for @xmath60 t corresponding to a maximum of the @xmath51-oscillations .
bottom inset in ( b ) sketches the dgg sample . ]
for thz excitation we applied a continuous wave ( _ cw _ ) molecular optically pumped laser  @xcite .
the laser operated at the frequency @xmath61 thz ( photon energy @xmath62 mev , wavelength @xmath63 m ) .
the incident power about 30  mw was modulated at about @xmath47 or 625  hz by an optical chopper .
the radiation at normal incidence was focused onto a spot of about @xmath64  mm diameter .
the spatial beam distribution had an almost gaussian profile , measured by a pyroelectric camera  @xcite .
taking into account the size of the superlattice we obtain that the power irradiating the structure is @xmath65mw .
the radiation intensity and electric field strength on the sample are @xmath66mw/@xmath18 and @xmath67v / cm , respectively .
note that almost all experimental data , except the scan of the beam across the sample , are obtained for radiation focused on the dgg superlattice .
the structures were placed in a temperature variable optical magneto - cryostat .
the photocurrents were studied in the temperature range between 2 and 9  k. the laser radiation was linearly polarized along the @xmath68-axis .
the photoresponse is measured by the voltage drop @xmath69 across a load resistor @xmath70  ohm  @xmath71 using standard lock - in technique . here
@xmath72 is the sample resistance .
the benefit of using of the small value of @xmath73 is that the detected signal is unaffected by the sample resistance variation and is just proportional to the electric current generated by the thz radiation .
the current is calculated via @xmath74 .
on the gate voltage @xmath52 ( @xmath57 ) applied to the stripes at zero potential of the other gate @xmath75 ( @xmath59 ) .
the polarization independent photocurrent contribution @xmath76 is extracted from the photocurrent polarization dependencies ( not shown ) .
the data are obtained for ( cd , mn)te qw sample # 1 at zero magnetic field . here and in the further plots , except for the one showing results on scan across the sample ,
the data are presented for the laser beam focused on the dgg superlattice .
insets show experimental geometries . here and in the further plots ,
the data are obtained for linearly polarized radiation with the electric field vector @xmath77 directed perpendicularly to the metal gate fingers . ]
we begin with briefly introducing the results obtained at zero magnetic field .
while the paper is devoted to the observation of magneto - ratchet effect these results are important for further analysis of the data on magnetic field induced photocurrents .    illuminating the dgg superlattice at normal incidence
we have observed a photocurrent , whose magnitude and direction are sensitive to the electrostatic potentials applied to the first ( g1 ) and the second ( g2 ) gate sublattices .
the photoresponse is detected perpendicularly to the gate stripes as well as along them .
the fact that the photoresponce is generated by normally incident radiation provides a first indication that it is formed due to the presence of the superlattice .
indeed it is well known that in unpatterned ( 001)-oriented qw structures a photocurrent signal is only detectable under oblique incidence  @xcite .
the photocurrent consists of a contribution independent of the radiation polarization and a substantially smaller one dependent on the orientation of electric field vector @xmath78 of the linearly polarized radiation in respect to the orientation of the gate stripes .
this has been demonstrated by measuring the photocurrent as a function the orientation of the @xmath79-field vector in respect to the @xmath68-axis ( not shown ) .
the same polarization behavior has been previously observed for the ratchet effect in lateral gaas - based structures  @xcite as well as in dgg structures fabricated on top of inalas / ingaas / inalas / inp high electron mobility transistors ( hemt )  @xcite and graphene  @xcite .    for unbiased gates , a non - zero built - in electrostatic potential
is formed due to the presence of metal stripes in the qw vicinity .
a direct evidence for the ratchet effect comes from the variation of the photocurrent by application of the bias voltages to the individual gates .
indeed the ratchet effect is known to be proportional to the averaged product @xmath80 of the derivative of the coordinate dependent electrostatic potential @xmath81 and the distribution of the electric near - field @xmath82  @xcite . in the following ,
we call the parameter @xmath83 the _ lateral asymmetry_. the value of @xmath83 may change sign depending on @xmath81 .
consequently a variation of individual gate voltages should result in a change of the ratchet current including reversal of its direction , see the discussion in sec .
[ discussion ] .
exactly this behavior has been observed in the experiment . in order to tune the lateral asymmetry
, we applied different bias voltages @xmath84 to the grating gates .
figure  [ gatedepb0 ] demonstrates the influence of the gate voltage variation on the amplitude and sign of the polarization independent photocurrent @xmath85  @xcite .
holding one of the gates at zero bias and varying the gate voltage on the other one we could controllably change the lateral asymmetry .
the figure reveals that inversion of this lateral asymmetry obtained either by change of polarity of voltage applied to one gate or by interchange of gate voltage polarities applied to narrow and wide gates results in a change of sign of the photocurrent .
this observation agrees well with the signature of ratchet currents : @xmath86 .
note that , while the overall characteristics of the data obtained at different cooldowns is the same , the photocurrent sign inversions can be shifted on the gate voltage scale by @xmath55  v. similarly to the results addressed in sec .
[ imp ] this is ascribed to cooldown dependent charge trapping in the insulator .
we also note that for large negative gate voltages all contributions vanish , which can be attributed to closing the transistor channel beneath the metal stripes of the dgg superlattice .
figure  [ gatedepb0 ] also reveals that at large positive gate voltages the photocurrent saturates due to the effect of high electron densities .    to summarize , experiments at zero magnetic field provide a consistent picture demonstrating that the photocurrents are caused by the ratchet effect .
they are ( i )  generated due to the lateral asymmetry , ( ii )  characterized by specific polarization dependencies for directions along and across the metal stripes , and ( iii )  change the direction upon reversing the lateral asymmetry .
these results are in full agreement with the theory of ratchet effects excited by the polarized thz electric field in asymmetric lateral superlattices , discussed below in sec .  [ discussion ] .
the overall behavior of the photocurrent is also in qualitative agreement with that of the electronic ratchet effects observed in semiconductor qw structures and graphene with a lateral superlattice  @xcite .
as a function of the magnetic field @xmath2 measured in ( cd , mn)te qw sample  # 1 for two combinations of the gate voltages .
solid blue line show the data for the gate voltages @xmath87 and @xmath88  v and the dashed red line for @xmath89  v and @xmath90 , the setups denoted as @xmath91  v and @xmath92  v. the vertical dashed line indicates zero magnetic field .
the inset illustrates the experimental setup . ]
now we turn to the main part of the paper devoted to magnetic field induced ratchet effects .
figure  [ jxjy4k ] shows the normalized photocurrent excited by linearly polarized radiation as a function of magnetic field applied perpendicularly to the qw plane .
the principal observation is that with raising magnetic field the photocurrent drastically increases and , at high magnetic fields , exhibits sign - alternating @xmath1-periodic oscillations .
moreover , fig .
[ jxjy4k ] demonstrates that the reversal of the lateral asymmetry results in an inversion of the magneto - photocurrent direction providing a first clear indication that it is driven by the ratchet effect .
measurements of the magneto - photocurrent as a function of radiation power show that it scales linearly with the radiation power , i.e. quadratically with the radiation electric field , @xmath93 , see inset in fig .
[ intensity ] .
note that the positions of oscillation maxima / minima are independent of the radiation power , fig .
[ intensity ] .
similar results , particularly magneto - oscillations with current amplitudes much larger than the photocurrent at @xmath94 , are obtained for all studied samples . studying various samples at different experimental conditions we found that the current in the @xmath95-direction is much smaller than @xmath49 .
data from the ( cd , mn)te qw sample  # 4 with ferromagnetic dy gate are shown in fig .
the period and magnitude of the @xmath1-oscillations correlate with transport behavior and are very close to those shown for ( cd , mn)te qw samples with gold dgg structure presented for sample  # 1 , fig .
[ jxjy4k ] .     measured at different levels of power ranging from 0.8 to 4.3  mw in ( cd , mn)te qw sample # 1 .
inset shows radiation power dependence of the photocurrent amplitude @xmath96 for @xmath60  t. solid line is a linear fit . ]
measured as a function of magnetic field @xmath2 in ( cd , mn)te qw sample  # 4 with the dy grating . ]
figure  [ cdte ] shows another example of magneto - photocurrent measured for nonmagnetic cdte qw sample  # 5 with gold dgg . while oscillations are clearly detected the signal is superimposed with a substantial background current . because of the background the photocurrent , in contrast to samples # 1 and # 4 , does nt change sign with varying magnetic field .
the background photocurrent has also been obtained for ( cd , mn)te qw samples # 2 and # 3 .
we attribute an appearance of the background to imperfections of our large - size superlattices which are present in some structures as confirmed by optical microscope images  @xcite .
the imperfections locally reduce symmetry and give rise to ratchet unrelated magneto - gyrotropic photogalvanic currents in the `` bulk '' of the qw , see e.g. ref .
@xcite .     as a function of magnetic field @xmath2 in cdte - based qw sample  # 5 at two gate voltage sequences : @xmath97
v@xmath98 and @xmath99@xmath1002.5  v. panels ( b ) and ( c ) show dependencies of the photocurrent magnitude on the gate voltage @xmath52 ( @xmath57 ) obtained for zero potential on the other gate @xmath57 ( @xmath52 ) .
the data in panels ( b ) and ( c ) are obtained for the photocurrent maxima / minima at @xmath101 and 4.9  t , respectively . ]      to prove that the magneto - ratchet effect stems from the irradiation of the superlattice we scanned the laser spot across the sample along the @xmath68-direction .
the photocurrent @xmath49 was measured for @xmath102 t and the gate voltage combination @xmath1030.5  v , i.e. at a maximum of one of the magneto - oscillations , fig .
[ jxjy4k ] .
the photocurrent generated by linearly polarized radiation with @xmath104 as a function of the radiation spot position and the corresponding experimental geometry are shown in fig .
the photocurrent reaches its maximum for the laser spot centered at the superlattice and rapidly decays with the spot shifting away . for the beam spot touching the sample edge
the current increases again .
this is caused by an edge photocurrent  @xcite known for graphene  @xcite and semiconductor qws  @xcite .
comparison with the laser beam spatial distribution measured by a pyroelectric camera ( dashed curve in fig .  [ scan ] ) shows that the change of the photocurrent for the scans across the dgg structure only slightly deviates from the gaussian intensity profile .
note that the dgg area is smaller than the beam spot .
these findings unambiguously demonstrate that the photocurrent is caused by irradiation of the superlattice .     measured in ( cd , mn)te qw sample  # 1 at a maximum of the photocurrent oscillation at @xmath102  t and @xmath87
, @xmath88 v. inset shows the setup with the laser spot scanned along the @xmath68-axis .
positions @xmath105 and @xmath106 correspond to the borders of the superlattice whereas @xmath107 denotes the sample edge position .
the dashed line represents the laser beam spatial distribution which is measured by a pyroelectric camera and scaled to the current maximum . ]      in order to explore the role of the lateral asymmetry , a prerequisite for the ratchet photocurrents , we have systematically studied magneto - photocurrents for different combinations of the gate voltages .
the results are shown for different samples in figs .
[ gatedepb0 ] ,  [ jxjy4k ] ,  [ cdte ] and  [ asymmetry ] .
( @xmath57 ) obtained for zero potential on the other gate @xmath57 ( @xmath52 ) .
the data are obtained for ( cd , mn)te qw sample  # 1 at @xmath60  t and @xmath108  t , panels ( a ) and ( b ) , respectively .
the inset shows the magnetic field dependence of the photocurrent measured for different gate voltage combinations .
vertical dashed lines demonstrate that the maxima / minima positions of the magneto - photocurrent are not shifted for all gate voltage combinations .
the numbers in the inset correspond to : 1 @xmath109 @xmath110  v ; 2 @xmath109 @xmath111  v ; 3 @xmath109 @xmath112  v ; 4 @xmath109 @xmath113  v ; 5 @xmath109 @xmath114  v and 6 @xmath109 @xmath115  v. ]    while the magnitude and the sign of the oscillations are strongly affected by the variation of the gate potentials the maxima / minima positions remain unchanged . to show that this conclusion is valid for the whole range of gate potentials used in our measurements we studied the photocurrent for a large set of different magnitudes and polarities of @xmath52 and @xmath116 . these results are presented in the inset in fig .
[ asymmetry ] for @xmath2 ranging from 2.7 to 4.1  t.    using the fact that maxima positions do not shift upon variation of @xmath2 we fixed the magnetic field and controllably varied the lateral asymmetry by changing the potential applied to one of the gates and holding the other one at zero bias .
figures  [ asymmetry](a ) and  [ asymmetry](b ) show the results obtained for ( cd , mn)te qw sample  # 1 at magnetic field strengths 3.9 and @xmath108  t , respectively .
the figures reveal that the polarity of photocurrent changes with the sign of the asymmetry parameter @xmath83 , eq .  .
this is primarily seen by comparing the curves obtained by variation of the potentials applied either to the thin ( dashed lines ) or thick ( solid lines ) gate stripes .
the observed stronger change of the signal for @xmath117  v corresponds to an increase of the longitudinal resistance detected in magneto - transport measurements , see inset in fig .
[ transport](b ) .
it should be noted that , due to the abovementioned built - in potential , the ratchet current can be detected even for zero gate voltages at both gates .
the same behavior upon variation of the lateral asymmetry has been observed for dy - based dgg sample  # 4 . studying the photocurrent as a function of the gate voltage in cdte qw sample  # 5 and in ( cd , mn)te samples  # 2 and  # 3 we obtained , that , despite the background
, the photocurrent magnitude is also controlled by the lateral asymmetry .
corresponding results obtained for cdte qw dgg structure  # 5 at @xmath1184.9  t are shown in figs .
[ cdte](b ) and  ( c ) .     as
a function of @xmath2 measured in ( cd , mn)te qw dgg sample  # 1 for three temperatures and two gate voltage combinations . ]
now we turn to the analysis of the @xmath1-oscillations of the photocurrent which are obviously related to the sdh oscillations .
comparison with the magneto - transport data reveals that the period of the @xmath1-oscillations is indeed equal to that of the longitudinal magneto - resistance oscillations .
this result is confirmed for all samples and temperatures used in our experiments .
in the cdte samples and ( cd , mn)te sample at relatively high temperatures the amplitude of oscillations grows with the magnetic field , see e.g. fig .
[ temperatur ] .
by contrast , in diluted magnetic semiconductor ( cd , mn)te qw dgg samples at liquid helium temperature @xmath119 the rise of the magneto - photocurrent amplitude at low fields is followed by its reduction at higher fields , see fig .
[ temperatur](b ) . at the lower temperature @xmath120
k this reduction , now detected at even lower magnetic fields , is followed by a further increase of the photocurrent magnitude at higher @xmath2 , see data in fig .
[ temperatur](c ) .
similar `` beats '' of the magneto - oscillation amplitudes are detected for the longitudinal resistance , fig .
[ transport ] , in agreement with the presence of the exchange enhanced zeeman splitting in dms materials  @xcite .
the enhanced zeeman splitting is caused by the exchange interaction between electrons and mn@xmath21 ions and described by the modified brillouin function .
the spin splitting first grows linearly with @xmath2 and then saturates .
the temperature increase leads to a decrease of exchange interaction and shifts the saturation to higher fields .
all these features are detected in both transport and magneto - photocurrent experiments .
moreover , due to the enhanced spin splitting the oscillations measured at low temperatures and high magnetic fields become spin - split which results in the beats .
in low - dimensional semiconductor structures with lateral superlattices , a @xmath0 electric current is generated due to the action of an electromagnetic wave s @xmath121 electric field  @xcite .
the effect of the superlattice is twofold : it generates a one - dimensional periodic electrostatic potential @xmath81 acting upon the 2d carriers with @xmath68 being the superlattice principal axis , and causes a periodic spatial modulation of the thz electric field due to the near field diffraction .
figure  [ field ] shows coordinate dependence of the thz electric near - field @xmath122 calculated for the ( cd , mn)te qws based dgg structure for two combinations of the gate voltages .
the electric field distribution caused by the near - field diffraction is calculated for radiation with frequency @xmath123  thz applying a self - consistent electromagnetic approach based on the integral equation method described in detail in ref .
figure  [ absorption ] shows the calculated thz absorption spectrum of the structure .
it is seen that the plasmonic resonances in the qw are excited at frequencies well below the operation frequency of 2.54  thz .
therefore in the experimental frequency range the absorption follows the drude law , see the inset in fig
.  [ absorption ] . as a result , the near - field distribution in the qw shown in fig .
[ field ] is caused almost solely by the metal dual - grating gate and not by the plasma oscillations .
hence one can neglect the dependence of the near field on the voltage applied to the dgg .
moreover , since the radiation frequency is much higher than that of plasmonic resonances we consider only electronic ratchet mechanism below .
normalized to the electric field of the incident thz wave .
shaded regions correspond to the positions of dgg metal fingers .
results are obtained for parameters of ( cd , mn)te qw sample # 1 with dielectric constant @xmath124 at the frequency 2.54  thz .
the momentum relaxation time @xmath125 is obtained from the electron mobility specified in sec .  [ sample ] . ]    . ]     on the gate voltages @xmath52 ( dashed line ) and @xmath57 ( solid ) fitted to experimental curves shown in fig .
[ gatedepb0 ] ( @xmath46 is defined in fig .  [ structure ] ) .
the lines are obtained after eq .
( [ j_zero_b ] ) . to be specific
, we used for calculations @xmath77 in accordance with the experimental conditions .
for the calculations we used the ( cd , mn)te dgg qw structure parameters discussed in sec .
[ sample ] , electric near - field @xmath122 shown in fig .
[ field ] , lattice temperature @xmath126  k , frequency @xmath61  thz and effective mass @xmath127 .
the only fitting parameter used here is the relaxation time @xmath128  ns .
inset shows the calculated frequency dependence of the photocurrent @xmath49 together with the experimental result for @xmath61  thz ( dot ) . ]
the ratchet current is generated by a combined action on carriers of both the radiation near field , transmitted through the grating , and the periodic static potential @xmath81  @xcite .
a crucial condition for @xmath0 electric current generation is that the lateral superlattice is asymmetric .
applying the voltages @xmath129 and @xmath130 in a controllable way , one can vary the potential @xmath81 and reverse the lateral asymmetry @xmath83 .
consequently , the ratchet photocurrent can also reverse its sign as observed in experiment , see fig .
[ gatedepb0 ] .
the ratchet current has components both perpendicular and parallel to the metal fingers and consists , in general , of polarization - independent , linear- and circular - ratchet contributions .
all these individual ratchet currents have been observed and described in various semiconductor - based structures  @xcite as well as in graphene with lateral periodical top gates  @xcite .
the mechanism leading to photocurrent formation can be illustrated on the polarization - independent photocurrent caused by the seebeck ratchet ( thermoratchet ) effect .
this contribution dominates at zero magnetic field ( sec .
[ erstes ] , fig .
[ gatedepb0 ] ) and , as we show below , is most relevant to the observed magnetic ratchet current .
the seebeck ratchet effect is caused by inhomogeneous heating of two - dimensional electron gas and subsequent relaxation of the electron temperature . due to the near - field space modulation , the field @xmath82 acting upon electrons depends on the coordinate @xmath68 , fig .
[ field ] .
the field heats the electron gas changing the effective electron temperature locally to @xmath131 .
here we assume that , within a short time , the electron distribution becomes a quasi equilibrium fermi - dirac distribution with the temperature @xmath132 .
the space - modulated temperature correction is defined by the energy balance and can be found from the following equation  @xcite @xmath133 here @xmath134 is the boltzmann constant , @xmath135 is the temperature relaxation time , @xmath125 is the momentum relaxation time , @xmath136 is the electron charge , and @xmath137 is the electron effective mass .
we take into account that , in the present experiment @xmath138 , where @xmath139 is the cyclotron frequency .    as a result of the inhomogeneous heating , electrons diffuse from warmer to colder regions , and form a nonequilibrium density profile given by @xmath140 where @xmath33 is the average electron density , @xmath141 is the fermi energy . the ratchet current at zero magnetic field can be represented as a drift current of the electrons in the electric field @xmath142 of the space modulated electrostatic potential @xmath81 . from ohm s law @xmath143 with @xmath144 we have @xmath145 where the bar denotes averaging over the superlattice period .
finally , the ratchet current density at zero magnetic field is given by  @xcite @xmath146 where @xmath83 is given by eq .  .
figure  [ popovfateev ] shows the photocurrent calculated after eq .   together with the data of fig .
[ gatedepb0 ] .
the result of fig .
[ popovfateev ] is obtained for the near - field @xmath122 shown in fig .
[ field ] .
parameters used for the calculations correspond to the sample  # 1 .
the only adjustable parameter is the temperature relaxation time taken as @xmath128  ns . this value agrees with the temperature relaxation time of few nanoseconds reported for a degenerate electron gas in gaas - based heterojunctions  @xcite .
figure  [ popovfateev ] shows qualitative and reasonably good quantitative agreement between theory and experiment .
the above analysis suggests that , in the absence of magnetic field , the ratchet current is rather weak for degenerate electrons because , according to eq .
, the relative concentration correction @xmath147 is much smaller than the temperature correction @xmath148 . below
we show that the situation changes drastically in the presence of a quantizing magnetic field @xmath149 .
the magnetic field induced ratchet currents in the directions perpendicular , @xmath150 , and parallel , @xmath151 , to the stripes of the superlattice are given by  @xcite @xmath152 ^ 2 }          { \sinh{z}-z\cosh{z } \over \sinh^2{z } } \delta,\ ] ] @xmath153 here the magneto - oscillations are described by the factor @xmath154 ( @xmath155 ) @xmath156 while the dependence on @xmath157 describes the oscillation suppression with increasing temperature , eqs .
- are derived in the first order in @xmath154 . as compared with the theory developed in ref .
@xcite we take into account the spin subband splitting @xmath158 due to the zeeman effect . for the dms based ratchet structures the exchange enhanced zeeman splitting is given by  @xcite @xmath159 where @xmath160 is the electron land factor in the absence of magnetic impurities , @xmath161 is bohr s magneton , @xmath162 is the mn @xmath23-factor , @xmath163 is the mn - spin system temperature , @xmath164 and @xmath165 account for the mn - mn antiferromagnetic interaction , @xmath166 is the modified brillouin function , and @xmath167 is the exchange integral .
equations  ,   predict that the ratchet current components @xmath150 and @xmath168 have opposite parity upon magnetic field reversal due to the photocurrent deflection by the lorentz force  @xcite .
the ratchet current becomes an oscillating function of the magnetic field as the landau levels move through the fermi level and the photocurrent @xmath1-oscillation period is the same as for sdh oscillations .
equation  ( [ j ] ) explains the reason for the ratchet current oscillations with giantly enhanced magnitude as compared to the ratchet current at zero magnetic field .
the crucial issue is that the ratchet current in a quantized magnetic field arises due to heating - induced corrections to the conductivity rather than a variation of the electron density , and the conductivity in this regime is extremely sensitive to electron temperature variations . indeed , even a weak change of the electron temperature near the dingle temperature results in exponentially strong changes of the conductivity  @xcite .
therefore , the sensitivity of the conductivity to electron gas heating results in a huge enhancement of the ratchet current .
the ratio @xmath169 is governed by a factor @xmath170 and can reach about two orders of magnitude  @xcite .
it should be stressed that this enhancement occurs in moderate fields where the resistance exhibits weak sdh oscillations .
heating damps the conductivity oscillations , therefore , depending on the sign of the correction , heating can either increase or decrease the conductivity because ( @xmath171 ) changes sign as a function of @xmath2 .
as a result , the heating - induced ratchet current is an oscillating function of the field with a zero mean value .
and @xmath151 calculated after eqs .
,   for ( cd , mn)te qw sample at @xmath126  k. the structure parameters are relevant to the experiment presented in fig .  [ jxjy4k ] : @xmath172  mev , @xmath173 and @xmath174 . for other parameters we used the literature values for _
n_-(cd , mn)te : @xmath175 , @xmath176  mev , see ref .
@xcite , and @xmath127 , see ref .
for the calculations we also used @xmath177 and literature data of @xmath178  k and @xmath179 for @xmath180 , see ref .
solid and dashed curves show the results for two asymmetries of equal magnitudes but opposite signs of @xmath83 .
the oscillation amplitude is much higher than the ratchet current @xmath181 at zero magnetic field ( the latter is not seen at the chosen scale ) . ]    the calculated magnetic field dependence of the ratchet currents @xmath150 and @xmath151 are shown in fig .
[ theoryjxjy ] .
the results are presented for two gate voltage combinations with equal magnitude but opposite sign of @xmath83 .
equations  -([jy ] ) reveal the most important features of the magnetic ratchet effect : the current is proportional to @xmath182 , oscillates around zero with the same period as the longitudinal magnetoresistance , the amplitude of the current oscillations is immensely larger than the ratchet current at zero magnetic field , and the sign of the photocurrent reverses by changing the sign of @xmath183 .
all these features have been observed experimentally , see figs .
[ jxjy4k]-[cdte ] , [ asymmetry ] , and [ temperatur ] .
equations  ( [ j])- also describe the more complicated oscillating ratchet current behavior observed in the dms structures at low temperatures , see figs .  [ temperatur ] and [ theorytemp ] for experiment and theory , respectively . indeed , due to the temperature dependence of the brillouin function at low temperatures the zeeman splitting becomes comparable to the landau levels separation .
calculated curves for three values of the lattice temperatures and , consequently , different contributions of the zeeman effect , are shown in fig .
[ theorytemp ] .
the calculations show , as in experiment , beats in the oscillations of @xmath184 .
they stem from the interplay of zeeman and landau splittings and become more pronounced with temperature decreasing .
comparison of figs .
[ temperatur ] and [ theorytemp ] demonstrates that eq .
( [ j ] ) describes the qualitative behavior of the photocurrent quite well .
we note that the beats in @xmath185 are substantially influenced by the mn spin system temperature , which , because of electron gas heating , can be higher than the lattice one  @xcite .
this increase due to the strong dependence of the brillouin function on @xmath163 , results in a reduced exchange part of @xmath158 in eq .
( [ deltaz ] ) .
this shifts the onset of the exchange enhanced zeeman splitting to higher magnetic fields .
allowance for this effect should further improve the agreement between experimental data and calculations .
calculated for ( cd , mn)te qw sample at three different temperatures .
curves are obtained for the same values of @xmath186 , @xmath11 , @xmath24 , @xmath23 , @xmath187 , @xmath137 , @xmath188 , @xmath165 and @xmath164 as used in fig . 14 .
solid and dashed curves show , similarly to fig .
[ theoryjxjy ] , the results for two opposite lateral asymmetries . ]
to summarize , we have experimentally demonstrated and theoretically explained the magnetic quantum ratchet effect , i.e. , the ratchet effect in quantizing magnetic fields . to generate the ratchet current we use an asymmetric interdigitated top gate superlattice on top of the investigated quantum wells .
electric currents driven by terahertz electric field exhibit sign - changing magneto - oscillations with an amplitude giantly enhanced as compared to the photocurrent at zero magnetic field . the current amplitude and direction can be controllably changed by the variation of voltages applied to individual gate sublattices . the effect is observed in structures with various qw compositions and superlattice parameters .
the photocurrent generation mechanism can be well described in terms of semiclassical theory of magnetic ratchet effects .
the observed effect is driven by the periodic modulation of the electron temperature caused by near field diffraction .
the theory of the seebeck ratchet effect in the presence of quantizing magnetic field shows that the ratchet current follows the oscillations of the longitudinal resistance . in the dms structures at low temperatures
the investigated effect is strongly influenced by the exchange enhanced zeeman splitting caused by the exchange interaction of electrons with mn@xmath21 ions .
finally we note , that while here we have dealt with a homogeneous magnetic field , the observed ratchet photocurrent in superlattices made of hard magnetic material ( dy ) opens a possibility to study the ratchet effect driven by an inhomogeneous periodic magnetic field recently suggested in ref .
this kind of experiments would apply remnant magnetization of dysprosium together with the enhanced magnetic properties of dms qws and is a future task .
we thank i. a. dmitriev and d. r. yakovlev for helpful discussions . the support from the dfg priority program sfb  689 , the volkswagen stiftung program and rfbr ( projects 15 - 02 - 02989 , 16 - 02 - 01037 , 16 - 02 - 00375 , and 17 - 52 - 53063 ) is gratefully acknowledged .
the work of dvf was supported by the grant from the president of the russian federation ( no .
mk-5447.2016.2 ) . the research in poland
was partially supported by the national science centre ( poland ) through grant no .
dec-2012/06/a / st3/00247 and by the foundation for polish science through the ira programme financed by eu within sg op programme .
f. jlicher , a. ajdari , and j. prost , rev .
phys . * 69 * , 1269 ( 1997 ) .
p. reimann , phys . rep . * 361 * , 57 ( 2002 ) .
h. linke ( ed . ) , _ ratchets and brownian motors : basics , experiments and applications _ , special issue appl .
phys . a : mater .
sci . process .
a * 75 * , 167 ( 2002 ) .
p. hnggi and f. marchesoni , rev .
phys . * 81 * , 387 ( 2009 ) .                          v. v. popov , d. v. fateev , e. l. ivchenko , and s. d. ganichev , phys .
b * 91 * , 235436 ( 2015 ) .
p. olbrich , j. karch , e. l. ivchenko , j. kamann , b. mrz , m. fehrenbacher , d. weiss , and s. d. ganichev , phys .
b * 83 * , 165320 ( 2011 ) .
p. olbrich , j. kamann , m. knig , j. munzert , l. tutsch , j. eroms , d. weiss , m .- h .
liu , l. e. golub , e. l. ivchenko , v. v. popov , d. v. fateev , k. v. mashinsky , f. fromm , th .
seyller , and s. d. ganichev , phys .
b * 93 * , 075422 ( 2016 ) .
t. watanabe , s. a. boubanga - tombet , y. tanimoto , d. fateev , v. popov , d. coquillat , w. knap , y. m. meziani , yuye wang , h. minamide , h. ito , and t. otsuji , ieee sensors j. * 3 * , 89 ( 2013 ) .
y. kurita , g. ducournau , d. coquillat , a. satou1 , k. kobayashi , s. boubanga tombet , y. m. meziani , v. v. popov , w. knap , t. suemitsu , and t. otsuji , appl .
104 * , 251114 ( 2014 ) .
s. a. boubanga - tombet , y. tanimoto , a. satou , t. suemitsu , y. wang , h. minamide , h. ito , d. v. fateev , v. v. popov , and t. otsuji , appl .
. lett . * 104 * , 262104 ( 2014 ) .
v. v. popov , d. v. fateev , t. otsuji , y. m. meziani , d. coquillat , and w. knap , appl .
. lett . * 99 * , 243504 ( 2011 ) .
p. faltermeier , p. olbrich , w. probst , l. schell , t. watanabe , s. a. boubanga - tombet , t. otsuji , and s. d. ganichev , j. appl . phys . * 118 * , 084301 ( 2015 )
. s. d. ganichev , w. weber , j. kiermaier , s. n. danilov , d. schuh , w. wegscheider , ch .
gerl , d. bougeard , g. abstreiter and w. prettl , j. appl .
physics * 103 * , 114504 ( 2008 ) .
s. n. danilov , b. wittmann , p. olbrich , w. eder , w. prettl , l. e. golub , e. v.beregulin , z. d. kvon , n. n. mikhailov , s. a. dvoretsky , v. a. shalygin , n. q. vinh , a. f. g. van der meer , b. murdin , and s. d. ganichev , j. appl . phys . * 105 * , 013106 ( 2009 ) .                p. olbrich , c. zoth , p. lutz , c. drexler , v. v. belkov , ya .
v. terentev , s. a. tarasenko , a. n. semenov , s. v. ivanov , d. r. yakovlev , t. wojtowicz , u. wurstbauer , d. schuh , and s. d.ganichev , phys .
b * 86 * , 085310 ( 2012 ) .      j. a. gaj , r. planel , and g. fishman , solid state commun . * 29 * , 435 ( 1979 ) .
j. k. furdyna , j. appl . phys . * 64 * , r29 ( 1988 ) .
t.  dietl , in _
handbook on semiconductors _
, vol.3b , ed . t.s .
moss ( north - holland , amsterdam , 1994 ) .
j. karch , p. olbrich , m.  schmalzbauer , c. zoth , c.  brinsteiner , m. fehrenbacher , u.  wurstbauer , m.m .
glazov , s. a. tarasenko , e. l. ivchenko , d. weiss , j.  eroms , and s. d.  ganichev , phys .
lett . * 97 * , 227402 ( 2010 ) .
j. karch , c. drexler , p. olbrich , m. fehrenbacher , m. hirmer , m. m. glazov , s. a. tarasenko , e. l. ivchenko , b. birkner , j. eroms , d. weiss , r. yakimova , s. lara - avila , s. kubatkin , m. ostler , t. seyller , and s. d. ganichev , phys .
* 107 * , 276601 ( 2011 ) .
m. m.glazov and s. d.ganichev , physics reports * 535 * , 101 ( 2014 ) .
m.dantscher , d.a.kozlov , m.t.scherr , s.gebert , j.brenfnger , m.v.durnev , s.a.tarasenko , v.v.belkov , n.n.mikhailov , s.a.dvoretsky , z.d.kvon , d.weiss , and s.d.ganichev , arxiv cond - mat : 1612.08854 ( 2016 ) .
k .- m.dantscher , d. a. kozlov , p. olbrich , c. zoth , p.faltermeier , m.lindner , g. v.budkin , s. a. tarasenko , v. v.belkov , z. d.kvon , n. n. mikhailov , s. a. dvoretsky , d.weiss , b. jenichen , and s. d.ganichev , phys .
b * 92 * , 165314 ( 2015 ) .