stochastic geometry has recently been used for the analysis and performance evaluation of wireless  ( ad hoc as well as cellular ) networks ; in this approach , one models node locations as a spatial point process , e.g. , homogeneous poisson point processes , and one computes various network statistics , e.g. , interference , successful transmission probability , coverage  ( or , outage ) probability etc . as spatial averages .
this often leads to tractable performance metrics that are amenable to parametric optimization with respect to network parameters ( node density , protocol parameters , etc . ) .
more precisely , this approach yields spatial averages of the performance metrics for given network parameters ; then the parameters can be chosen to optimize performance .
this approach takes a macroscopic view of the network with the underlying assumption that all nodes in the network have identical statistical characteristics .    in practice , due to randomness and heterogeneity in networks
, nodes need to adapt to local spatial and temporal conditions  ( e.g. , channel conditions and topology ) to reach optimum network wide performance .
for example , nodes in wireless lans adjust their window sizes based on acknowledgment feedback ; in cellular networks nodes are scheduled based on channel conditions and adapt their transmit powers based on the measured sinrs , which in turn depend on the transmit powers set by other nodes . in all such scenarios ,
distributed adaptive algorithms are used to reach a desired network wide operating point e.g. that maximizing some utility . while the behavior of such distributed optimization protocols is often well understood on a given topology , there are usually no analytical characterizations of the statistical properties of the optimal state in large random and heterogeneous networks .
the main aim of this work is to use stochastic geometry to study spatial adaptations of medium access control in aloha that aim at optimizing certain utilities . while we identify a utility for which stochastic geometry can be used to compute the spatial distribution of map and the expected utility , we are far from being able to do so for all types of utilities within the @xmath0-fair class and we discuss the difficulties to be faced .
let us start with a review of the state of the art on aloha .
wireless spectrum is well known to be a precious and scarce shared resource .
medium access control  ( mac ) algorithms are employed to coordinate access to the shared wireless medium .
an efficient mac protocol should ensure high system throughput , and should also distribute the available bandwidth fairly among the competing nodes .
the simplest of the mac protocols , aloha and slotted aloha , with a `` random access '' spirit , were introduced by abramson  @xcite and roberts  @xcite respectively . in these protocols ,
only one node could successfully transmit at a time .
reference  @xcite modeled node locations as spatial point processes , and also modeled channel fadings , interferences and sinr based reception .
this allowed for spatial reuse and multiple simultaneous successful transmissions depending on sinr levels at the corresponding receivers .
all the above protocols prescribe identical attempt probabilities for all the nodes .
reference  @xcite further proposed opportunistic aloha in which nodes transmission attempts are modulated by their channel conditions .    among the initial attempts of map adaptation in aloha , reference
@xcite analyzed protocol model and proposed stochastic approximation based strategies that were based on receiver feedback and were aimed at stabilizing the network .
references  @xcite also optimized nodes attempt probabilities  ( or thresholds ) in order to maximize the spatial density of successful transmissions .
reference  @xcite analyzed both plain and opportunistic aloha in a network where all the nodes communicate to one access point .
they assumed statistically identical rayleigh faded channels with no dependence on geometry  ( i.e. , no path loss components ) .
they demonstrated a paradoxical behavior where plain aloha yields better aggregate throughput than the opportunistic one .
reference  @xcite also studied optimal random access with sinr based reception .
however , they considered constant channel gains .
they developed a centralized algorithm that maximizes the network throughput , and also an algorithm that leads to max - min fair operation .
reference  @xcite modeled network as an undirected graph and studied aloha under the protocol model .
they designed distributed algorithms that are either proportional fair or max - min fair .
reference  @xcite built upon the model of  @xcite , and formulated the channel access problem as a non - cooperative game among users .
they considered throughput and delay as performance metrics and proposed pricing schemes that induce socially optimum behavior at equilibrium .
however , they set time average quantities  ( e.g. , throughput , delay ) as utilities  ( or costs ) , and concentrated on symmetric nash equilibria .
consequently , in their analysis , dependence on local conditions vanishes .    in none of the above aloha models ,
nodes account for both wireless channel randomness and local topology for making their random access decisions , as we do in the present paper .
there is a vast literature on the modeling of csma by stochastic geometry which will not be reviewed in detail here .
the very nature of this mac protocol is adaptive as each node senses the network and acts in order to ensure that certain exclusion rules are satisfied , namely that neighboring nodes do not access the channel simultaneously .
however , csma as such is designed to guarantee a reasonable scheduling , not to optimize any utility of the throughput .
the closest reference to our work is probably  @xcite where the authors study an adaptation of the exclusion range and of the transmit power of a csma node to the location of the closest interferer .
this adaptation aims at maximizing the mean number of nodes transmitting per unit time and space  ( while respecting the above exclusion rules ) .
this mean number is however only a surrogate of the rate . in addition , the adaptation is only w.r.t . to the location of the nearest interferer .
we study spatial adaptation of aloha in ad hoc networks .
the network setting is described in section  [ sec : network - model ] .
we consider quasi - static networks in which mobiles learn the topology , and incorporate this information in their medium access probability  ( map ) selection .
we consider the cases where nodes are benevolent and cooperate in a distributed way to maximize the global network throughput or to reach either a proportional fair or max - min fair sharing of the network resources .
we analyze the case where nodes account only for their closest interferers , for all nodes in a given ball around them , or even all nodes in the network .
section  [ sec : aggr - interf ] is focused on the distributed algorithms that maximize the aggregate throughput or lead to max - min fairness in such networks . in the proportional fair case ,
we show that nodes can compute the optimal maps as solutions to certain fixed point equations . in the maximum throughput case
, the optimal maps are obtained through a gibbs sampling based algorithm . in the max min case ,
the optimal maps are obtained as the solution of a convex optimization problem .
section  [ ss : stocg ] contains the stochastic geometry results .
the model features nodes forming a realization of a homogeneous poisson point process in the euclidean plane .
we compute the map distribution in such a network in the proportional fair case using shot noise field theory . to the best of our knowledge
, this distribution is the first example of successful combination of stochastic geometry and adaptive protocol design aimed at optimizing ceratin utility function within this aloha setting .
we also show that the mean value of the logarithm of the throughput obtained by a typical node can be derived from this distribution . finally , we discuss the difficulties to be faced in order to extend the result to other types of utilities .
the numerical results are gathered in section  [ ss : numerics ] .
the aim of this section is two - fold : 1 ) check the analytical results against simulation and 2 ) quantify the gains brought by adaption within this setting .
we model the ad - hoc wireless network as a set of transmitters and their corresponding receivers , all located in the euclidean plane .
this is often referred to as `` bipole model ''  ( * ? ? ?
* chapter  16 ) .
there are @xmath1 transmitter - receiver pairs communicating over a shared channel .
the transmitters follow the slotted version of the aloha medium access control  ( mac ) protocol  ( see figure  [ fig : bipole - model ] ) .
a transmitter , in each transmission attempt , sends one packet which occupies one slot .
each transmitter uses unit transmission power .
we assume that each node has an infinite backlog of packets to transmit to its receiver .
the euclidean distance between transmitter @xmath2 and receiver @xmath3 is @xmath4 , and the path - loss exponent is @xmath0  ( @xmath5 ) .
we also assume rayleigh faded channels with @xmath6 being the random fading between transmitter @xmath2 and receiver @xmath3 .
moreover , we assume that the random variables @xmath7 are independent and identically distributed with mean @xmath8 . thus all @xmath6s have cumulative distribution function  ( cdf ) @xmath9 with @xmath10 .
all the receivers are also subjected to white gaussian thermal noise with variance @xmath11 , which is also constant across slots .
we assume that a receiver successfully receives the packet of the corresponding transmitter if the received sinr exceeds a threshold @xmath12 .
let @xmath13 be the indicator variable indicating whether transmitter @xmath3 transmits in a slot , and @xmath14 be @xmath3 s medium access probability .
thus @xmath15 when node @xmath3 transmits , the received sinr at the corresponding receiver is @xmath16 then the probability of successful reception @xmath17 can be calculated as follows .
@xmath18 thus @xmath19 where @xmath20 further simplifying , @xmath21 then , the rate or throughput of transmitter @xmath3 is given by @xmath22 .
the thermal noise appears merely as a constant multiplicative factor in the expression for the successful transmission probability  ( see  ) .
moreover , in interference limited networks , the impact of thermal noise is negligible as compared to interference .
we focus on such networks , and thus we ignore the thermal noise factor throughout .
in this section , we analyze adaptations of spatial aloha that maximize aggregate throughput or achieve proportional fairness or max - min fairness .      the throughput maximizing medium access probabilities solve the following optimization problem .
we first argue that the optimum in the above optimization problem is attained at one of the vertices of the hypercube formed by the constraint set . to see this ,
suppose @xmath24^{\mathcal{n}}$ ] is an optimal solution , and @xmath25 for some @xmath26 . clearly , @xmath27
since the partial derivative is independent of @xmath14 , @xmath14 can be set to either @xmath28 or @xmath29 without reducing the value of the objective function .
this proves our claim . in the following we focus only on such extreme solutions .
then the above problem is equivalent to finding an @xmath30 such that @xmath31 if and only if @xmath32 is an optimal solution .
thus we are interested in @xmath33    [ [ an - iterative - solution ] ] an iterative solution + + + + + + + + + + + + + + + + + + + + +    we can pose this problem as a strategic form game with the users as players  @xcite . for each player
its action @xmath34 lies in @xmath35 , and the utility function @xmath36 is given by @xmath37 this is a potential game with the above objective function as the potential function  @xcite .
thus the best response dynamics converges to a nash equilibrium .
this algorithm can be implemented in a distributed fashion if each node @xmath3 knows @xmath38 for all @xmath2 , and also @xmath39 and @xmath40 for all @xmath41 after each iteration . however
, a nash equilibrium can be a suboptimal solution to the above optimization problem . to alleviate this problem
, we propose a gibbs sampler based distributed algorithm , wherein each node @xmath3 chooses action @xmath29 with probability @xmath42 the parameter @xmath43 is called the temperature .
the gibbs sampler dynamics converges to a steady state which is the gibbs distribution associated with the aggregate throughput and the temperature @xmath43  @xcite . in other words
, we are led to the following distribution on the action profiles : @xmath44 where @xmath45 is a normalizing constant .
when @xmath43 goes to @xmath28 in an appropriate way  ( i.e. , as @xmath46 , where @xmath47 is the time ) , the distribution @xmath48 converges to a dirac mass at the action profile @xmath49 with maximum aggregate utility if it is unique .
notice that the aggregate utility @xmath50 equals the aggregate throughput .
thus the action profile @xmath49 is a solution to the original throughput optimization problem .
the first two terms in the utility function @xmath51 can be seen as `` selfish '' part of user @xmath3 , whereas the last summation term is `` altruistic '' part .
the user makes a decision based on whether the `` selfish '' part dominates or viceversa .    in a quasi - static network where topology continuously changes ,
although at a slower time scale , different sets of nodes are likely to be scheduled to transmit under different topologies .
thus , in terms of long term performance , maximum throughput medium access is not grossly unfair .
the proportional fair medium access problem can be formulated as follows .
@xmath52 the objective function can be rewritten as @xmath53 we thus have a convex separable optimization problem .
the partial derivative of the objective function with respect to @xmath14 is @xmath54 which is continuous and decreasing in @xmath14 over @xmath55 $ ] .
we conclude that at optimality , for each @xmath56 , @xmath57 observe that user @xmath3 s optimal attempt probability is independent of the attempt probabilities of other users . in particular , if @xmath58 , user @xmath3 can perform iterations @xmath59 autonomously .
furthermore , @xmath60 clearly , @xmath61 , i.e. , @xmath62 is a contraction .
thus the fixed point iterations converge to the optimal @xmath14 starting from any @xmath63 $ ] .
the characterization of the optimal attempt probabilities reflects the altruistic behavior of users .
more precisely , user @xmath3 s attempt probability is a function of @xmath64 which are measures of @xmath3 s interference to all other users .
in particular , if @xmath65 , i.e. , if @xmath3 s transmission does not cause significant interference to the other users , then @xmath3 transmits in all the slots . unlike the throughput maximization problem
, there is no `` selfish '' component in the decision making rule .    as the target sinr @xmath66 , @xmath67 for all @xmath68 , and the proportional fair attempt probabilities satisfy @xmath69 for all @xmath26 .
this yields @xmath70 for all @xmath26 .
this is expected , because in the limiting case a transmission can succeed if and only if there is no other concurrent transmission .
this is hence aloha without spatial reuse , and it is well known that in this case , the optimal access probability is @xmath71 asymptotically  @xcite .
our analysis in this section follows  @xcite .
the max - min fair medium access problem can be formulated as @xmath72 where constraint functions are defined for all @xmath73^{\mathcal{n}}$ ] .
the following is an equivalent convex optimization problem  ( see  @xcite for details ) : @xmath74    the lagrange function of this problem is given by  @xcite @xmath75 with @xmath76 being the lagrange multipliers .
minimization of the lagrange function  ( which is concave in @xmath77 and @xmath78 ) gives @xmath79 wang and kar  @xcite suggest that the lagrange multipliers be updated using the gradient projection method . more precisely , for all @xmath26 , @xmath80^+ , \label{eqn : aggr - interf - lambda - update}\end{aligned}\ ] ] where @xmath81 is the step size at the @xmath82th iteration .
further more ,  ( * ? ? ?
* theorem  2 ) implies that a solution arbitrary close to an optimal solution can be reached via appropriate choice of step sizes . however , all the users need to exchange variables in order to perform updates .    finally , the _ directed link graph _ corresponding to our network is a directed graph in which each vertex stands for a user  ( i.e. , a transmitter - receiver pair ) in the network .
there is an edge from vertex @xmath3 to vertex @xmath2 in the directed link graph if transmission of user @xmath3 affects the success of transmission of user @xmath2 .
two vertices @xmath3 and @xmath2 are said to be connected if either of the following two conditions hold :    1 .
there is an edge from @xmath3 to @xmath2 or viceversa , 2 .
there are vertices @xmath83 such that @xmath84 and @xmath85 are connected for @xmath86 .
clearly , the directed link graph for our network model is a complete graph ; for any pair of vertices @xmath3 and @xmath2 there is an edge from @xmath3 to @xmath2 and also from @xmath2 to @xmath3 .
in particular , the directed link graph is a single _ strongly connected component
_  @xcite .
thus  ( * ? ? ?
* corollary  1 ) implies that the above optimization also obtains the lexicographic max - min fair medium access probabilities that yield identical rates for all the users .
note that a user needs to know the entire topology , and in a few cases , also needs to communicate with all the nodes to implement the adaptation rules developed in sections  [ sec : aggr - interf - max - thput]-[sec : aggr - interf - max - min - fair ] . in this section , we carry out analysis under the simplifying assumption that the aggregate interference at a receiver is dominated by the transmission from the closest interferer .
this is a reasonable approximation in a moderately dense network , specifically when the path loss attenuations are high . throughout this section
, we use the notation @xmath87 for all @xmath56 . in words , @xmath88 is the strongest interferer of node @xmath3 , and @xmath89 is the set of nodes to which node @xmath3 is the strongest interferer .
we assume that there is always a unique @xmath88 for each @xmath3 .
then , accounting only for the nearest interferer , the approximate probability of successful transmission for node @xmath3 is @xmath90 the analysis of sections  [ sec : aggr - interf ] can be adapted to this simplified scenario .
the throughput maximization problem can now be posed as follows .
@xmath91    as in section  [ sec : aggr - interf - max - thput ] , we can argue that some @xmath92 attains the optimal throughput . again , an equivalent optimization problem is @xmath93 or alternatively , @xmath94    we now formulate a strategic form game among users , with action sets @xmath35 and utility functions given by @xmath95 again a gibbs sampler based algorithm yields the optimal set of the transmitting users . also , user @xmath3 only needs to know the distances of user @xmath88 and all the receivers in @xmath89 and their actions to make its decision .
notice that user @xmath3 must choose @xmath96 if @xmath97 such users can set their actions to @xmath29 , and need not undergo gibbs sampler updates .
[ [ discussion ] ] discussion + + + + + + + + + +    consider a scenario where a node s closest interferer does not transmit ,
i.e. , has zero attempt probability .
nonetheless , this node always has an active closest interferer  ( unless there are no other nodes in the network ) .
a better approximation of the success probabilities , and hence of the throughput , is obtained by always accounting for the closest active interferer . towards this , let us define @xmath98 for all @xmath32 .
we are now faced with the following optimization problem .
@xmath99 we can now define users utility functions as follows .
@xmath100 where @xmath101 .
the analogous distributed algorithm  ( gibbs sampler based ) can again be shown to lead to the optimal solution .
we now aim to solve the following optimization problem . @xmath102
following the discussion in section  [ sec : aggr - interf - prop - fair ] , we obtain @xmath103 again , if @xmath104 , iterations @xmath59 converge to the optimal @xmath14 starting from any @xmath63 $ ] .
if @xmath105 , i.e. , if @xmath3 s transmission does not cause significant interference to the users for whom @xmath3 is closest interferer , then @xmath3 transmits in all the slots .
the same user may not transmit  ( in any slot ) under the maximum throughput objective if @xmath106 and @xmath88 and users in @xmath89 transmit .
we can have explicit formulae for the attempt probabilities in a few special cases .    1 .
suppose @xmath89 is singleton for each @xmath3 .
if @xmath107 , then @xmath108 2 .
_ linear network topology :  _ we now consider a scenario where nodes are placed along a line , say @xmath109 , and are indexed sequentially .
we also assume that for any node @xmath3 the potential interferers are the two immediate neighbors @xmath110 and @xmath111 .
this also amounts to assuming @xmath112 for all @xmath3 .
then , assuming @xmath113 , the proportional fair attempt probability of node @xmath3 satisfies @xmath114 this is quadratic equation in @xmath14 , which on solving gives @xmath115      similar to section  [ sec : aggr - interf - max - min - fair ] , the max - min fair medium access problem is @xmath116 again , the constraint functions are defined for all @xmath73^{\mathcal{n}}$ ] .
let us recall the definition of the directed link graph associated with the network . in this section
, we only account for the interference due the closest interferer .
thus there is an edge from vertex @xmath3 to vertex @xmath2 if and only if @xmath117 .
we assume that the directed link graph is connected  ( i.e. , all the vertices in the graph are connected with each other ) .
if it is not connected , the max - min fair medium access problem on the entire graph decomposes into separate max - min fair medium access problems on each of the connected subgraphs , which can be solved independently .
we now pursue the following convex optimization problem which is equivalent to the above max - min fair medium access optimization problem .
@xmath118 the last set of constraints along with the connected assumption  ( of the directed link graph ) forces @xmath119 to be equal for all @xmath26 .
this confirms equivalence to the initial optimization problem .
now the lagrange function is @xmath120 where @xmath121 , are the lagrange multipliers .
an approach similar to section  [ sec : aggr - interf - max - min - fair ] prescribes the following update rules . for all @xmath122 , @xmath123^+ , \label{eqn : close - interf - lambda - update } \\
\mu_{ij}(n+1 )   = & [ \mu_{ij}(n ) + \beta(n)(\theta_i - \theta_j)]^+ , \label{eqn : close - interf - mu - update}\end{aligned}\ ] ] where @xmath81 is the step size at the @xmath82th iteration as before .
observe that any user @xmath3 can perform updates  - via local information exchange .
more precisely , it needs to communicate only with user @xmath88 and the users in @xmath89 .
the aim of the present paper is primarily of theoretical nature and it is beyond our scope to discuss implementation issues .
let us however stress that the discussed adaptations of the map are implementable .
we will focus on the proportional fair case in view the main focus of the paper .
assume each receiver has a distinctive pilot signal with fixed power . since we assumed a quasi - static network in which the nodes move at a slower time scale , each node
can then learn the distance that separates it from a given receiver by listening to its pilot signal and by performing a time average over  ( so as to smooth out fading ) .
once this data is available for all receivers , a given transmitter can then solve the key fixed point equation that characterizes its optimal map .
distinctive pilot signals can be obtained by a collection of orthogonal codes chosen at random by the receivers . in practice , it is enough for a transmitter to detect the `` dominant '' receivers  ( i.e. , those within a certain distance to it ) , so that the scheme will work in an infinite network with a finite  ( properly tuned ) number of such codes .
we now assume that the transmitting nodes are scattered on the euclidian plane according to a homogeneous poisson point process of intensity @xmath124 . for each transmitter ,
its corresponding receiver is at distance @xmath125 in a random direction .
the traffic and channel models are the same as in section  [ sec : network - model ] . as before the transmitters use slotted aloha to access the channel , and a receiver successfully receives the packet from its transmitter if the received sinr exceeds a threshold @xmath12 .
finally , transmitters adapt their attempt probabilities as described in section  [ sec : aggr - interf ] .
each transmitter is associated with a multi dimensional mark that carries information about the adaptive transmission probability and the transmission status .
let @xmath126 denote a marked poisson point , where    * @xmath127 denotes the poisson point process of intensity @xmath124 , representing the location of transmitters in the euclidean plane .
* @xmath128 denote the marks of the poisson point process @xmath129 , which consist of three components : * * @xmath130 denote the angles from transmitters to receivers .
these angles are i.i.d . and
uniform on @xmath131 $ ] and independent of @xmath132 . we will call them the primary marks . * * @xmath133 denote the maps of the nodes ; @xmath134 is a secondary mark ( i.e. a functionals of @xmath132 and its primary marks , see below ) . *
* @xmath135 are indicator functions that take value one if a given node decides to transmit in a given time slot , and zero otherwise .
clearly , @xmath136 .
in particular , given @xmath134 , @xmath137 is independent of everything else including @xmath138 .
the locations of the receivers will be denoted by @xmath139 with @xmath140 .
it follows from the displacement theorem @xcite that @xmath141 is also a homogeneous poisson point process of intensity @xmath124 .
the above assumptions will be referred to as the poisson model .
we will also consider below a more general case where the above marked point process is just stationary .
let us consider response functions @xmath142 defined for each @xmath143 as follows @xmath144 for all @xmath145 , the shot noise field @xmath146 associated with the above response function and the marked point process @xmath129 is @xmath147 notice that this shot noise is not that representing the interference at @xmath148 .
it rather measures the effect of the presence of a transmitter at @xmath148 on the whole set of receivers .
consider a typical node at the origin , @xmath149 , with marks @xmath150 .
let @xmath151 denote the palm distribution of the stationary marked point process @xmath129  ( * ? ? ?
* chapter  1 ) .
the fixed point equation determining the map of node @xmath152 reads  ( see  ) @xmath153 we have a similar equation for each node and the sequence @xmath133 is readily seen to be a sequence of marks of @xmath132 and its primary marks .
it follows directly from monotonicity arguments that @xmath154 notice that we have not used the specific assumptions on the point process so far .
hence we have the following general connection between the optimal map distribution and the shot noise @xmath155 :    for all stationary marked point processes @xmath156  ( not necessarily poisson ) , for all @xmath157 , @xmath158 with @xmath151 the palm distribution of @xmath129 .
we now use the fact that @xmath129 is an independently marked poisson point process  ( * ? ? ?
* definition  2.1 ) . from slivnyak
s theorem  ( * ? ? ?
* theorem  1.13 ) , @xmath159 for all @xmath160 .
consequently , @xmath161    it follows from  ( * ? ? ?
* proposition  2.6 ) and from the fact that @xmath162 is a homogeneous poisson point process that one can write the laplace transform @xmath163 of the shot noise @xmath164 as @xmath165 where @xmath166 .
[ thm : map - distribution ] under the above poisson assumptions , the attempt probability of the typical node has the distribution @xmath167 with @xmath168 given by ( [ eq : lapl ] ) .
let @xmath169 denote the density of the shot noise field @xmath170 .
then @xmath171 where @xmath172 if @xmath173 and 0 otherwise .
now using parseval s theorem @xmath174 with @xmath175 the fourier transform of the real valued random variable @xmath176 and @xmath177 the complex conjugate of @xmath178 .
the claim follows after substituting @xmath179 and @xmath180 .
[ lem : laplace - transform ] for @xmath181 , the laplace transform @xmath182 can be simplified as @xmath183      this subsection is devoted to the analysis of the mean value of the logarithm of the throughput of the typical node : @xmath184 with @xmath185 defined in  . since we know the cdf @xmath186 of @xmath187 , the first term poses no problem .
the second term can be rewritten as  ( see  ) @xmath188.\ ] ] under the law @xmath151 , the points @xmath189 of @xmath132 form a homogeneous poisson point process of intensity @xmath124 .
however , the marks @xmath190 do not have the law identified in the last section .
in fact , the mark @xmath134 of a point @xmath191 satisfies the following modified fixed point equation : @xmath192 with the convention that @xmath193 if there is no solution in @xmath55 $ ] .
we can use the same argument as above to conclude that @xmath194 iff @xmath195    conditioned on there being two nodes at @xmath28 and @xmath148 , the other points form a homogeneous poisson point process of intensity @xmath124 .
this allows one to prove the following .
[ thm : map - typical - receiver ] under the above poisson assumptions , given that there is a node at 0 and a node at @xmath196 , the attempt probability of the node at @xmath148 has the distribution @xmath197 with @xmath198 and @xmath140 .
due to the circular symmetry , the first integral in the expression of @xmath199 depends on @xmath148 only through @xmath200 .
thus the density of @xmath201 also depends on @xmath200 only , and it will be denoted by @xmath202 when @xmath203 .
the density of @xmath187 identified in the last subsection will be denoted by @xmath186 .
the main result of this section is :    [ thm : mean - utility ] under the above poisson assumptions , in the proportional fair case , the mean utility of a typical node is @xmath204    see appendix  [ proof : thm - mean - utility ] .
a preliminary concern when trying to use euclidean stochastic geometry of the infinite plane in the maximum throughput and max - min fairness cases is that it is not clear whether the associated infinite dimensional optimization problems make sense in the first place . in the proportional fair case
, each node computes its optimizing map in one step as the solution of a fixed point equation that is almost surely well defined  ( in terms of a shot noise ) even in the infinite poisson population case .
unfortunately , this does not extend to the other two cases .
this does not mean that there is no hope at all .
consider for instance the maximum throughput case accounting only for the closest interferer , and further simplify it by measuring interference at the transmitter rather than at the receiver . in other words , consider the same optimization problem as in section  [ sec : close - interf - max - thput ] but with @xmath88 being the closest transmitter to transmitter @xmath3  ( rather than to receiver @xmath3 ) .
then , in the poisson case , the infinite dimensional optimization problem can be shown to reduce to a _
countable collection of finite optimization problems_. this follows from the fact that there are no infinite `` descending chains '' in a homogeneous poisson point process  ( * ? ? ?
* chapter  2 )  ( a _ descending chain _ is a sequence of nodes @xmath205 such that @xmath206 for all @xmath82 ) . as a result
, the poisson point process can be decomposed into a countable collection of finite `` descending trees '' , where each path from the leaves to the root is a descending chain .
the associated optimization is hence well defined and the problem can be reduced to evaluating the solution of the optimization problem of section  [ sec : close - interf - max - thput ] on the typical descending tree .
hence , there is hope to progress on this and on related cases .
this will however not be pursued here and is left for future research .
in this section , we study the proposed adaptive spatial aloha schemes quantitatively . we compute various metrics formulated through the stochastic geometry based analysis , and we also perform simulation .
the simulation not only validates the analytical model , but also illustrates the performance of schemes for which we do not have an analytical characterization .
we used maple and matlab to evaluate the integrals of section [ ss : stocg ] .
the infinite integral that shows up in the expression of the laplace transform ( [ eq : lapl ] ) is handled without truncation by maple and matlab .
the singularity at @xmath207 in the contour integrals ( [ eq : parsev ] ) leveraging parseval s theorem is a false singularity and it is also handled without further work by either maple or matlab .
the matlab code is particularly efficient and is used throughout the analytical evaluations described below .
we consider a two dimensional square plane with side length @xmath208 , and @xmath1 nodes placed independently over the plane according to the uniform distribution ; this corresponds to @xmath209 in the stochastic geometry model .
, the poisson random variable with mean @xmath210 is highly concentrated around its mean .
thus we can use @xmath210 nodes for all the realizations in our simulation . ]
each node has its receiver randomly located on the unit circle around it , again as per the uniform distribution .
thus @xmath211 for all @xmath3 .
we set @xmath181 and @xmath212 . to nullify the edge effect , we take into account only the nodes falling in the @xmath213 square around the center while computing various metrics . while all other parameters remain , we vary @xmath208 and @xmath1 for different simulations . for each parameter
set we calculate the average of the performance metric of interest over @xmath214 independent network realizations .
we validate the analytical expression against the simulation for the case of proportional fair medium access . for illustration , we plot the cumulative distribution function  ( c.d.f . ) of the map in figure  [ fig : map - cdf - analysis ] . here
we set @xmath215 and consider two values of @xmath1 , @xmath216 and @xmath217 , which correspond to @xmath218 and @xmath219 respectively .
the plots show that the stochastic geometry based formula  ( see theorem  [ thm : map - distribution ] ) quite accurately predicts the nodes behavior in simulation .
we also study the distributions @xmath220  ( also referred to as @xmath202 for @xmath203 ) defined in theorem  [ thm : map - typical - receiver ] .
figure  [ fig : map - cdf - receiver ] shows their plots for @xmath219 and two values of @xmath200 , @xmath221 and @xmath222 .
again the plots based on the analytical expression and those based on simulation closely match .
clearly , under @xmath151 , nodes closer to the origin are more likely to be inactive .    , the map distribution of a node at distance @xmath223 from origin , under @xmath151.,width=288,height=240 ]      in this section we illustrate the performance of various adaptive schemes and their benefit over plain aloha .
we compute the performance metrics via simulation and also through analytical expressions whenever the latter are available . in such cases the analytical results and the simulation
validate each other .
first we set @xmath224 and @xmath225 .
we consider the maximum throughput medium access , however , only accounting for the closest interferers . in figure
[ fig : gibbs - closei ] , we show steady state behavior of our gibbs sampling based algorithms ; we have set the temperature @xmath226 . as expected , the improved maximum throughput medium access  ( see the discussion at the end of section  [ sec : close - interf - max - thput ] ) insures a better exclusion behavior . under this scheme ,
a lesser number of nodes transmit , and neighboring nodes are unlikely to transmit simultaneously .
so , this is expected to deliver better aggregate throughput .
now we keep @xmath208 fixed at @xmath227 , but vary @xmath1 from @xmath228 to @xmath229 ; this corresponds to varying @xmath124 from @xmath230 to @xmath231 in the analytical expressions .
we evaluate the aggregate throughputs of various aloha schemes including plain aloha .
the average throughputs are plotted in figure  [ fig : throughput ] . _
although in some of the schemes we derive the attempt probabilities only considering the closest interferers , we always take into account the aggregate interference while calculating the throughput_. when the number of nodes is small , both the throughput maximizing medium access and plain aloha have identical performance ; both prescribe attempt probabilities close to one for all the nodes .
when the number of nodes increases beyond @xmath232 , the throughput maximizing medium access significantly underestimates the interference , and thus its performance deteriorates . on the other hand ,
the aggregate interference based proportional fair scheme significantly outperforms plain aloha in terms of aggregate throughput also .
this benefit is sustained even as the number of nodes increases .
we also notice that the improved version of maximum throughput medium access  ( see the discussion at the end of section  [ sec : close - interf - max - thput ] ) yields best performance among all the schemes , and its performance does not deteriorate until a much higher number of nodes .        in figure
[ fig : log - throughput ] , we plot the logarithms of the aggregate throughputs corresponding to the two proportional fair medium access schemes .
the figure illustrates that as the number of nodes increases , the performance of the closest interferer based medium access worsens in comparison to the performance of the aggregate interference based scheme - this is not visible merely looking at the corresponding aggregate throughputs  ( see figure  [ fig : throughput ] ) ) .        in figure
[ fig : map - cdf ] , we plot the c.d.f . of the map for the proportional fair case .
our objective is to compare the case when nodes account for the aggregate interference with when they account for the closest interferer only .
we set @xmath215 and plot map distributions corresponding to two values of @xmath1 , @xmath233 and @xmath234  ( corresponding to @xmath219 and @xmath218 respectively ) .
as expected , the nodes attempt more aggressively when they account for the closest interferer only . while there is a favored probability interval in the aggregate interference case ,
there are more than one such intervals in the closest interferer case .
also , in the latter case , about @xmath235 of nodes attempt in almost all slots , irrespective of @xmath1 , the total number of nodes .
this can be understood by noticing that the probability that a node is not the closest interferer to any other node is not sensitive to @xmath1 .
we have shown the feasibility of the performance analysis of distributed adaptive protocols that aim at maximizing some global utility in a large random network using stochastic geometry .
more precisely , the most natural distributed adaptation of the medium access probability of aloha that aims at proportional fairness optimization was shown to have a tractable optimal map distribution .
this distribution is obtained from the law of a certain shot noise field that describes the interference created by a typical node to all receivers but his . in the poisson case ,
the distribution of the optimal map is obtained as a non - singular contour integral which is amenable to an efficient evaluation using classical numerical tools .
the network performance at optimum can in turn be deduced from the latter using campbell s formula .
this approach is shown to provide an analytic way of quantifying the gains brought by this proportionally fair adaptive version of aloha compared to plain aloha .
this line of thoughts opens several research directions .
the first one is the extension to other types of fairness , still in the framework of aloha .
we indicate that this is possible at least under certain simplifications of the interference model .
the second and broader question is whether this approach can be extended to mac protocols other than aloha .
an example would be an adaptation of the exclusion radius of csma / ca to the full environment of a node aiming at maximizing some utility of the throughput .
a third general question concerns the evaluation of the `` price of decentralization '' .
when the discussed protocols are suboptimal because of their greedy / distributed nature , is it possible to use stochastic geometry to evaluate the typical discrepancy between the performance of the distributed scheme and the optimal centralized one ?
this work was carried out at laboratory of information , networking and communication sciences  ( lincs ) paris , and was supported by inria - alcatel lucent bell labs joint research center .    10    n.  abramson . :
another alternative for computer communications . in _ proceedings of the fall joint computer conference  ( afips ) _ , montvale , new jersey , november 1970 .
f.  baccelli and b.  blaszczyszyn . , volume 3  issue  3 - 4 of _ foundations and trends in networking_. now publishers , 2009 .
f.  baccelli and b.  blaszczyszyn . ,
volume 4  issue  1 - 2 of _ foundations and trends in networking_. now publishers , 2009 .
f.  baccelli , b.  blaszczyszyn , and p.  muhlethaler .
an aloha protocol for multihop mobile wireless networks .
, 52(2):421436 , february 2006 .
f.  baccelli , b.  blaszczyszyn , and p.  muhlethaler .
stochastic analysis of spatial and opportunistic aloha .
, 27(7):11051119 , september 2009 .
f.  baccelli , j.  li , t.  richardson , s.  shakkottai , s.  subramanian , and x.  wu . on optimizing csma for wide area ad hoc networks . ,
72(1 - 2):3168 , 2012 .
d.  p. bertsekas and r.  g. gallager .
. prentice hall , 1992 .
s.  boyd and l.  vandenberghe . .
cambridge university press , 2004 .
p.  bremaud . .
number  31 in texts in applied mathematics .
springer - verlag , 1999 .
m.  franceschetti and r.  meester . .
cambridge university press , 2008 .
b.  hajek .
stochastic approximation methods for decentralized control of multiaccess communications . ,
31(2):176184 , march 1985 .
m.  k. hanawal , e.  altman , and f.  baccelli .
stochastic geometry based medium access games . in _
ieee infocom _ ,
orlando , florida , usa , march 2012 .
hsu and h .- j .
channel - aware aloha with sinr capture : when is the knowledge of channel not helpful ?
in _ proceedings of ieee pimrc _ , toronto , canada , september 2011 .
mohsenian - rad , v.  w. wong , and r.  schober .
optimal sinr - based random access . in _ proceedings of ieee infocom _ ,
san diego , ca , usa , march 2010 .
d.  monderer and l.  s. shapley .
potential games . , 14(1):124143 , may 1996 .
m.  osborne and a.  rubinstein . . the mit press , 1999 .
l.  g. roberts .
packet system with and without slots and capture .
, 5(2):2842 , april 1975 .    x.  wang and k.  kar .
distributed algorithms for max - min fair rate allocation in aloha networks . in _
42nd annual allerton conference on communication , control , and computing _ , urbana , il , usa , september - october 2004 .    x.  wang , k.  kar , and j .- s .
lexicographic max - min fair rate allocation in random access wireless networks . in _
45th ieee conference on decision and control _ , san diego , ca , usa , december 2006 .
let @xmath236 be a point process with marks and @xmath237 be the point process @xmath238 shifted by @xmath239 . given that @xmath240 , @xmath241 where the mapping @xmath242 associates with @xmath243 and the real number @xmath45 the solution of @xmath244 with the usual convention if there is no solution in @xmath55 $ ] .
it follows from slivnyak s theorem that @xmath245=\\ & & \hspace{-.7 cm } \mathbb{e}^0 \left[\sum_{n\ne 0 } \log \left(1-\frac{h(||x_n - y_0||,\phi^r\setminus \{y_0,y_n\}\circ\theta_{x_n})}{\frac{\vert x_n - y_0\vert^{\alpha}}{t r_0^{\alpha } } + 1}\right ) \right]=\\ & & \hspace{-.7cm}\frac 1 { 2\pi } \int\limits_0^{2\pi } \mathbb{e } \sum_{n } \log \left(1-\frac{h(||x_n - ( r_0,\phi)||,\phi^r \setminus \{y_n\}\circ\theta_{x_n } ) } { \frac{\vert x_n-(r_0,\phi)\vert^{\alpha}}{t r_0^{\alpha } } + 1}\right ) { \rm d}\phi.\end{aligned}\ ] ] it now follows from campbell s formula that @xmath246\\ & & \hspace{-.6 cm } = \int_{\mathbb{r}^2 } \int_v   \log\left(1-\frac{v\overline r_0}{||x-(r_0,\phi)||^\alpha+\overline r_0}\right ) \lambda { \rm d}x\\ & & \hspace{3 cm } \mathbb{p}^0 ( h(||x-(r_0,\phi)||,\phi^r\setminus\{y_0\})= { \rm d}v)\\ & & \hspace{-.6 cm } = \int_{\mathbb{r}^2 } \int_v   \log\left(1-\frac{v\overline r_0}{||x-(r_0,\phi)||^\alpha+\overline r_0}\right ) ) \lambda { \rm d}x \\ & & \hspace{3 cm } \mathbb{p } ( h(||x-(r_0,\phi)||,\phi^r)= { \rm d}v),\end{aligned}\ ] ] where the last relation follows from slivnyak s theorem .


recent years have seen an increased interest in many - particle spin dynamics , also termed _ spintronics
_  @xcite .
most of this research have been devoted to condensed matter systems , where the electron spin constitutes the spin degree of freedom  @xcite .
however , related or equivalent phenomena may occur by considering the polarization of light  @xcite or internal zeeman levels of cold atoms exposed to optical light fields  @xcite or spatially varying magnetic fields  @xcite .    in many of the spintronics schemes developed in the condensed matter community , a spin - orbit ( so ) coupling plays a crucial role for the dynamics .
the particular form of the so coupling can be of either rashba  @xcite or dresselhaus  @xcite type , both being frequently analyzed in terms of an effective gauge force  @xcite . in other words , one can interpret the so couplings as giving rise to an effective spin - dependent lorentz - type force .
thus , the so coupling can serve as a spin - filter or a stern - gerlach apparatus . as the manufacturing techniques improves , undesired effects deriving from crystal impurities can be greatly suppressed and several successful experiments demonstrating the so coupling and its outings have been reported recently  @xcite .
ultracold atoms in optical lattices  @xcite possess another possibility to implement effective so couplings , and the study of it within a clean environment . while in condensed matter settings where the coupling is an intrinsic property of the systems , for cold atoms it has to be externally created via spatially varying laser fields .
ever since the pioneering experiments on bose - einstein condensation ( bec )  @xcite , the field has witnessed a tremendous progress both experimentally and theoretically  @xcite .
the dispersive dipole coupling between light and atoms enables the possibilities to study matter waves in periodic potentials and the experimental realization of model systems discussed within theoretical solid state physics  @xcite .
today , both the weakly interacting atomic gases as well as more strongly interacting ones can be realized experimentally .
furthermore , utilizing the internal atomic hyperfine levels enable studies of the interplay between internal and external degrees of freedom  @xcite .
one interesting feature is the realization of artificial gauge potentials by means of coupling of internal atomic zeeman levels  @xcite .
the technique of light induced gauge potentials allows for non - abelian behavior to be analyzed and the so couplings are one such example  @xcite .    in this paper , we investigate the dynamics of a spinor bec subjected both to a two - dimensional weakly tilted optical square lattice as well to a set of spatially varying light fields giving rise to an effective so coupling between the internal atomic states .
thereby , its dynamics will be driven by two mechanisms , namely ; ( i ) bloch oscillations ( bos ) arising due to the tilting of the lattice , and ( ii ) an intrinsic spin hall effect ( she ) deriving from the effective so coupling . the two phenomena may be analyzed in terms of an extrinsic force along the longitudinal direction of the lattice tilt and an intrinsic transverse force .
being orthogonal to each other , these forces render a complex correlated motion of the bec in the lattice potential .
furthermore , the particular form of the so coupling renders dirac cones characterized by linear dispersions reminiscent of relativistic electrons . as an outcome of these
, the evolution possesses as well a trembling _ zitterbewegung _ motion  @xcite .
normally , such effect damps out very fast @xcite , but we show that in a lattice situations as in the present paper it may indeed be long - lived and hence more easily observable . moreover , the optical lattice can serve as trapping potential for the atoms in the bo regime and no confining harmonic trap is needed which would typically prevent buildup of spin currents  @xcite .    while so couplings arise as an intrinsic property in many crystals , in cold atom systems it has to be implemented externally .
this , at the one hand may be a source for experimental errors , but at the other it gives a handle for the so coupling ; its strength and form can be tuned experimentally .
the same holds for the optical lattice , the lattice strength and geometry are both control parameters , which has made it possible to observe bos for becs in optical lattices  @xcite . due to the large atomic mass compared to the mass of the electron ,
zitterbewegung is expectedly more clearly manifested in atomic systems .
the first experiments detecting zitterbewegung in atomic systems were recently carried out  @xcite .
a different but related setup can be found in ref .
@xcite , where a spinor condensate accelerating in a periodic magnetic field was studied .
we note that in ref .
@xcite , the spatially dependent magnetic field gives rise to an effective gauge force .
however , bloch oscillations were not considered , and moreover , in the present work the periodic potential is external and independent of the so coupling .
the outline of the paper is as follows . in sec .
[ sec2 ] , we introduce the model system and outline how effective so couplings are rendered in such systems . to start with
we analyze the free case in which the problem is analytically solvable and demonstrate the she and the zitterbewegung effect .
applying the acceleration theorem we also solve the problem under constant acceleration .
the dispersions , once an optical lattice has been included , are solved numerically displaying the dirac cones .
section  [ sec3 ] presents the main analysis .
the full numerical results are presented , where we demonstrate the presence of a she and how it is modified by the longitudinal oscillatory motion . in particular , it is shown that the hall current will also oscillate , however with an increasing oscillation amplitude .
on top of the oscillations deriving from the motion in a periodic potential , we find a persisting zitterbewegung effect . finally , we conclude with a summary in sec .
rashba  @xcite and dresselhaus  @xcite so couplings appear in a variety of solids , ranging from semiconductors  @xcite , thin magnetic films  @xcite , to graphene  @xcite .
the so couplings bring about a degeneracy in the energy dispersions . around these degeneracies ,
the dispersions are linear , reminiscent to that of relativistic particles , and therefore they are called _
dirac cones_. the relativistic properties have been confirmed experimentally in graphene either by _ angle resolved photoemission spectroscopy _ or by considering the conductivity and cyclotron motion of electrons  @xcite .
the so coupling and its corresponding degeneracy induces as well a non - trivial effective gauge structure , giving rise to phenomena such as the intrinsic spin and anomalous hall effects  @xcite . in the field of cold atoms ,
the first light induced effective magnetic fields were demonstrated in ref .
@xcite , and signatures of non - abelian structures were recently observed , see  @xcite and reference therein .
bos with ultracold atoms have already been realized experimentally by numerous groups  @xcite .
we proceed by presenting a scheme which realizes the desired so coupling in an atomic tripod system .      the idea of effective gauge potentials arising from adiabatic evolution dates back to the seminal works by mead on vibrations of molecules  @xcite .
berry  @xcite generalized this finding by demonstrating the existence of an effective gauge structure in the context a quantum adiabatic evolution .
wilczek and zee @xcite further generalized berry s result and found an effective non - abelian gauge potential accompanying adiabatic evolution of degenerate energy eigenstates .
dum and olshanii @xcite showed that the effective gauge potentials can be generated in systems of ultracold atoms subjected to spatially dependent light fields that couple internal atomic states . following an approach for creating non - abelian gauge potentials in time - dependent systems  @xcite , ruseckas _ et al . _
@xcite outlined a scheme how to achieve effective non - abelian gauge potentials by considering an atomic tripod system in spatially dependent light fields .
the coupled four - level system possesses two degenerate dark states that , within certain parameter regimes , can be decoupled from the remaining two bright states via standard adiabatic elimination .
the flexibility of this tripod system , deriving from different laser configurations , allows for a variety of phenomena to be studied , such as , e.g. , relativistic properties @xcite , she s @xcite , novel phases of the bec @xcite , and topological defects as effective magnetic monopoles @xcite or vortices  @xcite .        in the present paper
we consider a tripod setup similar to the one described in refs .
the details of atomic and laser configurations rendering the gauge potentials can be seen in fig .
the internal dynamics , is driven by three lasers that couple the three lower atomic states @xmath0 , @xmath1 , and @xmath2 to the excited state @xmath3 . parameterizing the laser fields as @xmath4 , @xmath5 , and @xmath6 , such that @xmath7 , the two degenerate dark states of the corresponding hamiltonian can be written as @xmath8 e^{-is_{13 } } |1\rangle \\
\\   & & + \left [ e^{-i\pi/4 } \cos \ ! \theta \sin \ !
\varphi - e^{i\pi/4 } \cos \varphi \right ] e^{-is_{23 } } |2\rangle \\ \\   & & - e^{-i\pi/4 } \sin\!\theta |3\rangle \big ) , \\ \\
\displaystyle{\frac{1}{\sqrt{2 } } } \big ( \left [ e^{-i\pi/4 } \cos \ ! \theta \cos \ ! \varphi - e^{i\pi/4 } \sin\varphi \right ] e^{-is_{13 } } |1\rangle \\ \\   & & + \left [ e^{-i\pi/4 } \cos \ ! \theta \sin \ ! \varphi + e^{i\pi/4 } \cos\varphi \right ] e^{-is_{23 } } |2\rangle \\ \\   & & - e^{-i\pi/4 } \sin \ !
\theta |3\rangle \big ) , \\ \\
\end{array}\ ] ] where @xmath9 . provided the field amplitude @xmath10 is much larger than the atom - light detuning @xmath11 ( as depicted in fig .
[ fig4 ] ) , the two dark states can be safely eliminated from the two bright states , resulting in an effective two - level system  @xcite . to derive an effective so coupling between the two internal dark states @xmath12 and @xmath13 , we choose @xmath14 , @xmath15 , @xmath16 , and @xmath17 $ ] .
this can be achieved by coupling the two ground states @xmath0 and @xmath1 to the excited state @xmath3 with two lasers intersecting with an angle @xmath18 and propagating in the @xmath19-plane .
these lasers correspond to a standing wave in the @xmath20-direction and a plane wave in the @xmath21-direction .
a third laser , couples the states @xmath2 and @xmath3 , and is a plane wave in the @xmath22-plane with an angle @xmath23 of propagation relative to the @xmath21-axis . in terms of the angles @xmath18 and @xmath23 , and the wave numbers @xmath24 , @xmath25 , we have @xmath26 with @xmath27 . performing the adiabatic elimination , taking into account that @xmath28 and the fact that we consider ultracold atoms so that excitations due to atomic motion can be neglected , the effective single particle two - level hamiltonian becomes in dimensionless variables @xcite @xmath29 where @xmath30 is an effective zeeman splitting , @xmath31 , @xmath32 .
the pauli matrices @xmath33 , @xmath34 , operate in the space spanned by the dark states , i.e. , @xmath35 , @xmath36 , and @xmath37 .
the @xmath38-dependence of the last two terms originates from non - adiabatic couplings  @xcite .
we note that these do not include @xmath39 and as a consequence the @xmath40-dependence can be fully factored out . in eq .
( [ ham4 ] ) , we have used dimensionless parameters .
explicitly , when we further on introduce a lattice potential whose wave number is @xmath41 , we let the characteristic length scale be @xmath42 and we let the characteristic energy scale be the recoil energy @xmath43 . consequently
, time @xmath44 is measured in units of @xmath45 .
the superscript @xmath46 indicates the free hamiltonian in absence of a spatially dependent optical lattice potentials .
the zeeman shift @xmath30 splits the degeneracy of the two internal states and together with the so coupling it gives rise to a transverse current , an intrinsic anomalous hall effect .
we will , however , focus on the intrinsic she obtained when @xmath47 . in this case ,
the transverse particle current vanishes but , on the other hand , a transverse spin current persists .
the @xmath47 condition can be achieved by internal stark shift of the individual bare atomic levels  @xcite .
thus , in the following we will assume @xmath47 , and furthermore we will also consider the symmetric so coupling characterized by @xmath48 .
it should be pointed out , however , that even by lifting these constrains an interplay between transverse and longitudinal atomic motion will still occur .
one important prerequisite for successful application of the current scheme is that the system can be prepared in the dark subspace spanned by @xmath49 .
this issue has been discussed in  @xcite .
the idea is to slowly turn on the external lasers @xmath50 , @xmath51 , such that the state follows adiabatically states in the dark subspace , preparing a superposition of @xmath49 but no bright state components .
we note that we may rewrite the hamiltonian in eq .
( [ ham4 ] ) ( @xmath47 ) as @xmath52 with the vector potential @xmath53 and the scalar potential @xmath54 these potentials constitute vector and scalar parts of the effective gauge potential .
although being constant in space , the vector potential gives rise to a non - trivial effective gauge field due to its non - abelian structure , @xmath55\neq0 $ ] .
explicitly , this effective field is given by the antisymmetric field tensor @xmath56 , \,\ , \ k , l = x , y,\ ] ] resulting in @xmath57 and @xmath58 otherwise . as a result , the spin - orbital motion of the electrons can be understood in terms of an effective lorentz - type force , akin to that found by wong  @xcite for a color - charged classical particle moving in a non - abelian gauge field .
the corresponding force operator may be derived from heisenberg equations of motion according to  @xcite @xmath59}\\ \\   & = & \displaystyle{v_{so}^2 ( \hat{\bf{p } } \times { \bf e}_z ) \hat{\sigma}_z } , \end{array}\ ] ] where it is understood that the operators on the right - hand side are all given in the heisenberg representation and the prefactor @xmath60 in front of the acceleration is the mass in the chosen units .
note that the force in the @xmath21-direction is proportional to @xmath61 and the force in the @xmath20-direction is proportional to @xmath62 .
this demonstrates the transverse character of @xmath63 and its equivalence to a lorentz or coriolis force . the crucial aspect of eq .
( [ lforce ] ) is , however , that the force is state dependent : atoms in dark state @xmath12 feel a force with opposite sign compared to those in dark state @xmath13 .
therefore , the motion of the particles in the @xmath19-plane may separate the two spin states implying that a non - zero spin current is generated . in the following
, we shall call the @xmath20-direction longitudinal and the @xmath21-direction transverse .
the effect of the above lorentz - type gauge force can be examined by computing how the initial ( @xmath64 ) state @xmath65 .
we assume that the internal and momentum states of the atom are initially uncorrelated and that the initial momentum states are gaussian wave packets of the form @xmath66{4\pi\delta_p^2 } } e^{-\frac{(p_x - p_0)^2+p_y^2}{2\delta_p^2}},\end{aligned}\ ] ] where @xmath67 sets the momentum uncertainty and @xmath68 is the average momentum along the @xmath20-direction .
this choice corresponds to the pure internal atomic state @xmath69 .
given this initial state , the time evolved state is readily obtained as @xmath70 with @xmath71 \nonumber \\   & & \times \psi(p_x , p_y ) e^{-i(p_x^2+p_y^2)t } , \nonumber \\
\psi_2(p_x , p_y , t ) & = & -\frac{1}{\sqrt{2}}\left [ \cos(\lambda t ) + i\sin(\lambda t)e^{i\theta}\right ] \nonumber \\   & & \times \psi(p_x , p_y ) e^{-i(p_x^2+p_y^2)t } .\end{aligned}\ ] ] here , we have used the energy dispersions @xmath72 and the energy eigenstates @xmath73 of the hamiltonian in eq .
( [ ham4 ] ) , which read @xmath74 where @xmath75 and @xmath76 .    in fig .
[ fig1 ] we display the dimensionless averages @xmath77 with @xmath78 .
the factor @xmath79 derives from normalization of the individual wave packet components .
since the initial velocity is along the longitudinal @xmath20-direction , the deviation from @xmath80 is a direct consequence of the state - dependent lorentz - type gauge force @xmath63 .
note that @xmath81 giving a vanishing transverse particle current .
this symmetry is broken by keeping the zeeman term @xmath82 in the hamiltonian of eq .
( [ ham4 ] ) .
the particle current induced by such a @xmath83-term goes under the name intrinsic anomalous hall effect .
after @xmath84 , the separation of the two internal dark states has reached an asymptotic value . until @xmath84 ,
the transverse spin current shows an oscillating behavior deriving from the same mechanism as the one underlying relativistic zitterbewegung in the dirac equation  @xcite . in the dirac equation
the trembling motion describes coherent coupling between positive and negative energy solutions , while in the present situation it describes population swapping between the internal states . in solids ,
the amplitude of these oscillations become extremely small due to the lightness of the electron , while in atomic or ionic systems the heavier masses imply larger amplitudes and possible detection  @xcite . despite heavier masses , the rapid damping of the trembling motion as depicted in fig .
1 is a limitation also in atomic / ionic zitterbwegung experiments  @xcite . in the next subsection
we will argue how this decay may be prevented by studying the effect in the presence of an optical lattice .
( green curve correspond to the dark state @xmath12 and the red curve correspond to the dark state @xmath13 ) as defined in eq .
( [ analyexp ] ) .
the dimensionless parameters in eqs .
( [ ham4 ] ) and ( [ gaussin ] ) are @xmath85 ( @xmath86 ) , @xmath47 , @xmath87 , and @xmath88 .
the initial state is thus a gaussian with a non - zero average momentum along the @xmath20-direction.,width=302 ]    in order to generate a current in the system one applies a constant longitudinal force @xmath89 , e.g. , induced by gravitation in atomic systems . in the particular case of solids ,
the force derives from a voltage across the solid , and the electrons reach a steady state motion in which the applied force is balanced by scattering against crystal impurities . on the other hand ,
a system of cold atoms is , to a good approximation , free from impurities and one would therefore achieve constant acceleration . assuming that the force is weak and that the total time of evolution is relatively short
, we can perform an adiabatic approximation , consisting in replacing @xmath90 by @xmath91 in eq .
( [ timesol ] ) . hence , the momentum grows linearly in time in accordance with classical motion .
how this comes about can be easily understood by considering the time evolution of the hamiltonian @xmath92 for small forces @xmath93 and times @xmath44 . in this regime
, the time - evolution operator can be approximated accordingly @xmath94 so that the first exponential gives the coupled time - evolution while the second simply boosts the @xmath20-momentum by @xmath95 . for longer time periods ,
the large velocity rendered by the acceleration will typically cause non - adiabatic excitations , which we neglect for the moment .
similarly , for large values of @xmath96 , the factorization of the exponentials is no longer justified .
the resulting averages @xmath97 calculated within this approximation are depicted in fig .
[ fig2 ] . as before
, we encounter a build - up of a transverse spin current but in this case the current changes sign after some time and finally there is no net spin separation . by comparison with a full numerical integration of the hamiltonian
we find that the imposed adiabatic approximation breaks down after relatively short times .
nevertheless , numerically we have verified that the dynamics beyond this approximation goes as well asymptotically to zero for large times .
we believe that the origin of this result is a complicated interplay between the transverse and longitudinal forces , together with the coupling between the two internal states . for spinless particles ,
the two forces render a spirally transverse motion and hence a net transverse current . for particles with spin - dependent lorentz forces @xmath63
one would expect a net transverse spin current , but the coupled nature of the two spin states seems to destroy such spin current in the case of constant acceleration .
, but with a constant weak force of magnitude @xmath93 in the @xmath20-direction added to the hamiltonian ( [ ham4 ] ) and assuming the validity of the adiabatic approximation .
the dimensionless parameters are @xmath85 , @xmath47 , @xmath98 , @xmath88 , and @xmath99.,width=302 ]      in the following section we numerically investigate the dynamics in the presence of an optical lattice .
insight into the evolution behavior is gained by considering the spectrum of the hamiltonian including a lattice potential while disregarding the linear force potential .
adding to the so hamiltonian in eq .
( [ ham4 ] ) a spatially oscillating periodic potential , yields @xmath100 the hamiltonian is periodic and consequently possesses a band - gap spectrum @xmath101 with associated bloch eigenstates @xmath102 , characterized by a discrete band - index @xmath103 and quasi - momentum numbers @xmath104 restricted to the first brillouin zone , i.e. , @xmath105 and @xmath106 . a detailed study of bloch theory in the presence of so couplings can be found in ref .  @xcite .
the first two bands @xmath107 and @xmath108 , obtained from numerical diagonalization of the hamiltonian , are shown in fig .
a similar band structure has been obtained in ref .
@xcite for cold atoms in a square optical lattice .
we directly see the effect of the so coupling in terms of dirac cones at @xmath109 , where @xmath110 .
it is further seen that the ground state is not the @xmath111 bloch state , as for a regular square lattice lacking so couplings .
this fact is at the heart of the _ jahn - teller effect _
@xcite , where the dirac cone ( or _ conical intersection _ as it is usually referred to when it appears in position space rather than in momentum space ) induces a symmetry breaking that lowers the ground state energy .
furthermore , we note that the two lowest bands are not gapped , while there is a gap between the second and third band . indeed , due to the so coupled two - level structure of the hamiltonian , the bands come in pairs where each pair possesses dirac cones for quasi - momenta @xmath109 .
the effective number of dirac cones is four ( easily seen by shifting the brillouin zone by @xmath60 in both directions ) , compared to two dirac cones encountered in honeycomb lattices , like in graphene .
we notice that recently it was shown how more complex lattice potentials may render single dirac cones within the first brillouin zone  @xcite .     for a square lattice with rashba spin - orbit coupling .
in this example , @xmath85 and @xmath112 . , width=302 ]    in the adiabatic regime , a non - zero extrinsic force @xmath113 is treated by using the same approximation as above , yielding @xmath114 in the presence of a periodic potential , this adiabatic assumption has become known as the _ acceleration theorem _
it implies that no population is transferred between the individual bands @xmath115 . as a direct consequence ,
an initial state with well localized quasi - momentum within one single band will show an oscillatory motion rather than a constantly accelerating one , once the extrinsic force @xmath113 has been turned on .
these are the well studied bloch oscillations  @xcite , which have been experimentally verified in a variety of systems , such as cold atomic gases  @xcite , becs  @xcite , semiconductor superlattices  @xcite , waveguide arrays  @xcite , and photonic crystals  @xcite .
bos are clearly an outcome of the periodic potential and its characteristic energy spectrum possessing gaps of forbidden energies at the center and edges of the brillouin zone . as the amplitude of the potential @xmath116 vanishes , the gaps close and the acceleration theorem breaks down .
in particular , for weak lattice amplitudes , the gap between the first two energy bands scales as @xmath117 . by linearizing the dispersions around the band gap and replacing @xmath118 by @xmath119 ,
one obtains a realization of the celebrated _ landau - zener model _
@xcite . here
, @xmath120 is the initial quasi - momentum in the @xmath20-direction .
this model is analytically solvable , giving the probability of population transfer between the two states of the corresponding bands .
it has been shown that the decay of bos can be correctly estimated using the landau - zener theory  @xcite .    in the present model ,
the situation becomes qualitatively different . at the dirac points
, the gap vanishes and the acceleration theorem can not be imposed . however , from fig .
[ fig3 ] it is seen that the dirac cones appear between the first two bands , and not between the second and third band .
hence , assuming a state initially prepared on the lowest band one may expect predominant population transfer between the @xmath121 and @xmath122 bands and not to higher bands .
thereby , further acceleration is hindered by the gap between the second and third band and spin - dependent bos may still be encountered . in other words ,
a modified acceleration theorem involving pair of bands could result .
in fact , as already mentioned , the acceleration theorem is a specific example of the more general adiabatic theorem , and in molecular and chemical physics it has long been known that adiabaticity breaks down in the vicinity of a conical intersection  @xcite . in this context , it should be kept in mind that adiabaticity is a widespread concept in physics , and the general definition concerns systems whose dynamics is governed by an explicitly time - dependent hamiltonian .
adiabaticity , as it appears in the acceleration theorem , is in this sense more reminiscent of the _ born - oppenheimer approximation _
@xcite , applicable to time - independent hamiltonians describing light and heavy quantal degrees of freedom corresponding here to the spin and the spatial motion , respectively , of the particles .
our numerical analysis will be restricted to a system consisting of weakly interacting atoms , for which a mean - field description is expected to capture the phenomena we are interested in . in the @xmath123-wave scattering regime , we thereby consider dynamics rendered by a spinor gross - pitaevskii equation  @xcite    @xmath124 \\ \\ & & + \left[\begin{array}{cc } g_{11}n_1(x , y , t)+g_{12}n_2(x , y , t ) & 0 \\ 0 & g_{12}n_1(x , y , t)+g_{22}n_2(x , y , t ) \end{array}\right ] \bigg\}\psi(x , y , t ) .
\end{array}\ ] ]    here , @xmath125= \left[\begin{array}{c}\psi_1(x , y , t)\\ \psi_2(x , y , t)\end{array}\right]\ ] ] is a two - component spinor wave function of the dark states , @xmath126 , @xmath93 the external longitudinal force , and @xmath127 , @xmath128 and @xmath129 are the effective @xmath123-wave scattering amplitudes related to the scattering lengths @xmath130 as @xmath131 , where @xmath132 is the width of the single site wave function in the @xmath40-direction .
this definition of @xmath133 implies that the wave function has been normalized as @xmath134 . in the basis of dark states ,
the scattering lengths are effective ones different from the regular hyperfine scattering lengths
. nevertheless , the ratio between the two scattering amplitudes should approximately be the same regardless of basis , and , in particular , for a spinor condensate such as that of @xmath135 one has @xmath136  @xcite .
we will restrict the analysis to repulsive interactions , @xmath137 , and we will furthermore assume @xmath138 @xcite .    in the absence of a lattice ,
the so coupling gives rise to a non - abelian mead - berry curvature @xcite , which results in a transverse splitting of spin components , the intrinsic she .
however , in figs .
[ fig1 ] and [ fig2 ] the motion along the longitudinal @xmath20-direction was either constant or monotonously accelerating , while for a weakly tilted lattice , where bos dominate the longitudinal motion , it is neither clear how the transverse spin current nor how the trembling zitterbewegung motion will be manifested .
for instance , will there be sufficient time for a transverse spin current to be established within a single bloch period ? in ref .
@xcite , the intrinsic anomalous hall effect was studied in a single particle system confined by an harmonic potential and hence the motion was as well oscillatory as in the case of bos .
the transverse current appeared as a rotation on top of the harmonic oscillatory motion , i.e. , in polar coordinates @xmath139 the gauge potential causes a non - zero current in the @xmath140-direction .
such current would indeed be zero for an so coupling possessing an abelian structure rather than non - abelian , e.g. for couplings on the form @xmath141 with @xmath142 .
it is clear that a corresponding rotational current can not be strictly encountered in the present system since exact polar symmetry is broken by the tilting of the lattice and by the lattice in itself .
the time evolution of the system is solved by employing the split - operator approach  @xcite .
the method relies on factorizing the time - evolution operator into a spatial part and a momentum part valid in the limit of infinitely small time - step propagation .
one convenient feature of this method is that the wave function is in principle obtained at any instant of time in both position @xmath143 and momentum @xmath144 space .
furthermore , for a system bounded from below , the approach is capable of giving an approximate ground state of the system by simply propagating it sufficiently long in imaginary time .    ,
@xmath145 , and @xmath146 for all plots , @xmath147 in the upper plots ( a ) and ( b ) , and @xmath85 in the lower plots ( c ) and ( d ) .
thus , the lower plots correspond to the dispersions depicted in fig .
, width=302 ]    our numerical considerations will focus on propagation of an initial ( approximated ) ground state for the untilted ( @xmath148 ) lattice .
that is , we consider the ground state of a spin - orbit coupled bec in an optical lattice , and then we introduce a quantum quench by suddenly switch on a linear force in the @xmath20-direction .
hence , we begin by propagating an ansatz function in imaginary time with zero extrinsic force @xmath148 .
when convergence of the imaginary time evolution has been reached , we switch to real time propagation , including the extrinsic force by adding the potential term @xmath149 to the hamiltonian .
as predicted by the acceleration theorem , for short times the weak forces @xmath93 will induce a linear increase in the quasi - momentum @xmath118 .
in contrast to bos in a regular square lattice , here the so coupling implies that also the quasi - momentum @xmath150 is indirectly affected by the force @xmath93 . at the same time
, the build up of motion in the @xmath21-direction induces an intrinsic force @xmath151 along the longitudinal @xmath20-direction of bos .
in addition , the linear dispersions induce a trembling zitterbewegung motion .
it is thereby likely to expect a complex motion arising from the interplay between the two forces related to she and bos and the outcomes from linear dispersions .
figure [ fig5 ] displays the absolute values of the initial ground state wave function achieved via imaginary time evolution for @xmath148 .
the momentum wave functions depicted in ( a ) and ( c ) demonstrate the jahn - teller effect as the ground state populates non - zero momentum states  @xcite .
the lower two plots ( c ) and ( d ) correspond to the parameters of the two lowest bloch bands shown in fig .
the fact that the ground state obtained from this numerical method is not symmetric , as seen from the plots , is believed to have a two - fold cause .
first , it seems to derive from the fact that the imaginary time evolution is taken only up to a point where the state evolution has approximately converged , and , secondly , from the fact that even small fluctuations can result in an asymmetric ground state , due the four - fold degeneracy of the ground state of our hamiltonian .
in fact , the complete convergence of the ground state has not been reached since the position space wave functions depicted in ( b ) and ( d ) are localized in contradiction to the true ground state having a bloch structure , and furthermore since the momentum wave functions have a finite width .
this , however , should not be seen as a shortcoming since normally one would cool the atomic gas in the presence of a confining potential which would result in a localized gas .
in addition , if the propagated state extends over the full lattice one would run into troubles when quantities such as the mean spatial position is calculated numerically .
note that the momentum wave functions of ( a ) and ( c ) are different from the quasi - momentum as is evident in the multiple bragg components separated by @xmath152 .
in other words , even though the momentum wave functions extend beyond the first brillouin zone , it populates , to a very good approximation , only the lowest band .
we also note that the position space ground state wave functions , shown in ( b ) and ( d ) , predominantly populate the individual lattice sites .
hence , only judging from figs .
[ fig5 ] ( b ) and ( d ) it would be difficult to predict the presence of a so coupling and the dirac cones that derive from it , which instead manifest themselves in momentum space in terms of a jahn - teller effect .
( a ) , average longitudinal position @xmath153 ( b ) , and average intrinsic force @xmath151 in the transverse direction ( c ) .
red lines correspond to the dark state @xmath12 while green lines to @xmath13 . for @xmath153 the two lines are identical because of symmetry reasons .
the dimensionless parameters are the same as in fig .
[ fig5 ] ( a ) and ( b ) , but now with a non - zero extrinsic force , @xmath99 . ,
width=302 ]    in order to analyze the interplay between the transverse and longitudinal forces , @xmath154 and @xmath113 , we calculate @xmath155 with @xmath78 corresponding to averages in terms of the two dark states @xmath12 and @xmath13 , respectively , @xmath156 , and the last quantity is the intrinsic force in the @xmath21-direction . we find that @xmath157 ( up to numerical accuracy ) and therefore plot only one of the two components .
our numerical results are presented in figs .
[ fig6 ] and [ fig7 ] .
the difference between the two figures lies in the scattering amplitude @xmath158 . in fig .
[ fig7 ] the effective scattering amplitude @xmath158 is already as large as 10 recoil energies @xmath159 and we can thereby conclude that the spin current is rather insensitive to the value of @xmath158 . it is generally known that strong non - linearity can qualitatively change system properties . in becs
it leads to phenomena such as soliton and vortex creation .
normally the greatest impact of the non - linearity is seen around curve crossings  @xcite . beyond a critical value of the non - linearity ( in our case @xmath160 ) , in one dimension swallowtail loops are formed in the vicinity of the curve crossings characterizing additional stationary non - linear states , e.g. solitons .
up to this critical value the dynamics seems rather insensitive at least on a mean - field level .
however , formation of these loops beyond the critical value implies breakdown of the acceleration theorem  @xcite . to date
, it seems that analysis of these effects have been restricted to one dimensional situations and it is therefore not fully understood how non - linearity affects the structure of the dirac cones appearing in two dimensions .
nonetheless , we have numerically found that the general structure of figs .
[ fig6 ] and [ fig7 ] drastically change for very large values of @xmath158 .
all present bo experiments , however , are performed in the weakly interacting regime where it is known that loops are absent , and we thereby restrict our analysis to these experimentally more relevant parameter ranges .
we see from both fig .
[ fig6 ] and [ fig7 ] that the @xmath20 and @xmath21 motions are indeed correlated ; when the bloch oscillation reaches its turning points , the spin separation in the @xmath21-direction is roughly maximal and as the velocity in the longitudinal @xmath20-direction changes sign , the spin separation decreases due to the reversed sign of the intrinsic force .
when the longitudinal velocity is the greatest ( equivalent to maxima of the intrinsic force @xmath151 ) , the spin current is reversed .
more precisely , as the atoms begin to accelerate in the tilted lattice the quasi - momentum @xmath118 increases linearly in time and a spin current is built up in the transverse direction . however , after half a bloch period the quasi - momentum reaches the brillouin boundary and the @xmath20-motion is stopped and then reversed .
consequently , also the spin current is reversed and one encounters a decreasing spin current in the @xmath21-direction .
finally , the spin current has dropped to zero and builds up in the opposite direction instead .
hence , the spin current follows the same oscillatory behavior as the bos .
one may argue that a large spin separation should be hindered by the alternating sign of the velocity in the @xmath20-direction . nonetheless , from our numerical results depicted in figs .
[ fig6 ] and [ fig7 ] , we note that the spin separation in the @xmath21-direction is relatively large .
indeed , @xmath161 exceeds several site - spacings and furthermore this separation amplitude increases in time .
the fact that neither the oscillatory spin nor the oscillatory longitudinal current show a clean perfect harmonic pattern derives from , at the one hand possible non - adiabatic effects and the breakdown of the acceleration theorem , and at the other that the transverse current induces an intrinsic force @xmath162 in the longitudinal direction which adds to the extrinsic force @xmath113 .
the full coupled dynamics is therefore rather complex , but nonetheless , our figures make clear the overall in - phase oscillations between the two currents .    on top of the oscillations of @xmath163 , there is a slight drift toward negative values of @xmath20 .
this is an outcome of non - adiabaticity and breakdown of the acceleration theorem .
however , this drift is rather small and the dynamics is therefore predominantly adiabatic .
this drift can be suppressed by increasing the lattice amplitude @xmath116 on the cost that the amplitude of longitudinal oscillations decreases .
it would not , however , imply a diminishing she since it is related to the longitudinal momentum and not position
.    a closer look at the two figures reveals that in the transverse motions there is an additional oscillating motion different from the spin current being in phase with the bos .
this is believed to be zitterbewegung clearly seen in the freely evolving wavepackets of fig .
contrary to that figure , here these oscillations do not collapse suggesting that it might be preferable to detect the zitterbewegung effect within a tilted lattice , compared to the case of freely propagating atoms .
the survival of zitterbewegung should as well derive from the alternating motion due to the optical lattice .
figures [ fig6 ] and [ fig7 ] have been calculated for an approximated initial ground state of the untilted lattice .
as already pointed out , experimentally the condensate might be prepared in the presence of an external trap , and then the optical lattice as well as the coupling of the internal states are turned on in such a way that the condensate is not in a true ground state of the coupled system . to analyze the effect of such excitations
, we have redone the simulations of figs .
[ fig6 ] and [ fig7 ] for an initial gaussian momentum state of a form depicted in eq .
( [ gaussin ] ) with @xmath98 and @xmath88 .
the results are very similar to the ones of figs .
[ fig6 ] and [ fig7 ] , and we therefore do not present them here .
we have also investigated the effect of a shallow harmonic potential in the transverse direction during the evolution in the tilted lattice .
again , a transverse spin current is obtained , however this time the amplitude of transverse oscillations is decreased due to the confining potential .
this reduction of spin separation would indeed be present in any experiment utilizing a confining trapping potential .
the idea of employing an optical lattice potential as a trapping potential circumvent such shorting . in order to achieve trapping in this case
one must stay within the bo regime , i.e. adiabatic . as figs .
[ fig6 ] and [ fig7 ] demonstrates , even though the atoms oscillate a relatively large spin current can build up during time .    in relation to the measurability of the spin separation
, we may note that even though @xmath164 is considerably larger than the lattice spacing , it does not imply that the two wave functions @xmath165 and @xmath166 do not spatially overlap .
in fact , the overlap is relatively large and fluorescent measurement of a single experimental run would therefore not contain sufficient information to reveal the spin current .
however , this is hardly never the case in true experimental realizations and instead an ensemble of measurements should be considered , where the ensemble average gives @xmath164 . since
this value is presumably larger than the lattice spacing , it should indeed be measurable  @xcite . in the cold atom community ,
time - of - flight measurements are often employed , in which it is the atomic momentum density that is being measured . if the separation of the momentum wave functions is non - zero , the position space wave functions @xmath165 and @xmath166 will set apart further during the time - of - flight .
if this separation is faster than the intrinsic wave function broadening , a finite time - of - flight measurement may enhance the detection efficiency .    , but with an increased atom - atom interaction
; @xmath167 . , width=302 ]
in this paper , we have considered an ultracold atomic system possessing three well - studied properties of condensed matter theory , bloch oscillations , intrinsic spin hall currents , and zitterbewegung.our motivation for this study comes from the fact that bloch oscillations have been verified in the same type of setups @xcite , light induced gauge potentials have been demonstrated @xcite , and zitterbewegung have been detected in ion - trap experiments  @xcite .    the spin - hall effect ( she ) and the bloch oscillations ( bo ) in our system are driven by two perpendicular forces : the extrinsic force appearing due to tilting of the lattice potential and the intrinsic force originating from the effective gauge potential caused by the spin - orbit ( so ) coupling .
we demonstrated that the interplay between these two forces gives rise to a correlated motion in the @xmath19-plane , and despite the fact that the bos are expected to prevent a build up of large transverse spin currents , we argued that the she should indeed be detectable .
we utilized system parameters that should be achievable in current experiments , e.g. , the amplitude of the lattice potential is five recoil energies compared to the experimentally utilized ones of ref .
@xcite ranging from one to 14 recoil energies , and the so coupling two recoil energies compared to 8 recoil energies in ref .
@xcite . for a potential amplitude as large as @xmath168 , the dispersions are almost flat and longitudinal
motion is almost frozen .
the momentum along the longitudinal direction do , however , oscillate within the brillouin size and , interestingly and somewhat surprisingly , one would therefore see an oscillating transverse spin hall current despite lack of longitudinal currents .
we furthermore showed that the effects are insensitive to experimentally relevant atom - atom interaction strengths .
entering a strongly interacting gas seemed to affect our results and we argue that this is likely to derive from a non - linearity induced breakdown of adiabaticity around the dirac points .
an additional interesting aspect being discussed is the analogue of relativistic zitterbewegung .
our numerical results indicate that it may be preferable to consider an optical lattice configuration when studying this effect since it persists over longer time periods .
the optical lattice convey as well the possibility to avoid confining trapping potentials preventing large spin currents .
we should also point out that the system presented in this paper allows for many extensions and different phenomena to be studied .
for example , vortex generation in lattices  @xcite , bose - einstein condensation into non - zero momentum states  @xcite , topologically ordered states of matter  @xcite , or more strongly correlated systems , by means of second quantization of the mean - field hamiltonian , may find realization in atoms possessing an internal tripod structure and trapped in optical lattices .
the flexibility permits as well for various extensions , for example time - dependent lattices or so couplings . a recent paper considered zitterbewegung of ultracold atoms for a driven so coupling  @xcite .
the same idea can be utilized in order to increase the transverse spin current in our model .
m. onoda , s. murakami , and n. nagaosa , phys .
lett . * 93 * , 083901 ( 2004 ) ; k. yu .
bliokh and yu .
p. bliokh , phys .
e * 70 * , 026605 ( 2004 ) ; k. yu .
bliokh and yu .
p. bliokh , phys .
a * 333 * , 181 ( 2004 ) ; k. yu .
bliokh and yu .
p. bliokh , nat .
photonics * 2 * , 748 ( 2008 ) . o. hosten and p. g. kwiat , science * 319 * , 787 ( 2008 ) .
a. m. dudarev , r. b. diener , i. carusotto , and q. niu , phys . rev
. lett . * 92 * , 153005 ( 2004 ) ; x .- j .
liu , x. liu , l. c. kwek , and c. h. oh , phys . rev .
lett . * 98 * , 026602 ( 2007 ) ; j. y. vaishnav , j. ruseckas , c. w. clark , and g. juzelunas , phys .
lett . * 101 * , 265302 ( 2008 ) ; x .- j .
liu , m. f. borunda , x. liu , and j. sinova , phys .
lett . * 102 * , 046402 ( 2009 ) .
m. h. anderson , j. r. enscher , m. r. matthews , c. e. wieman , and e. a. cornell , science * 269 * , 198 ( 1995 ) ; k. b. davies , m. o. mewes , m. r. andrews , n. j. vandruten , d. s. durfee , d. m. kurn , and w. ketterle , phys .
* 75 * , 3969 ( 1995 ) .      c. j. myatt , e. a. burt , r. w. ghrist , e. a. cornell , and c. e. wieman , phys .
lett . * 78 * , 586 ( 1997 ) ; j. stenger , s. inouye , d. m. stamper - kurn , h .- j .
miesner , a. p. chikkatur , and w. ketterle , nature * 396 * , 345 ( 1998 ) ; d. m. stamper - kurn , h. j. miesner , a. p. chikkatur , s. inouye , j. stenger , and w. ketterle , phys . rev . lett . * 83 * , 661 ( 1999 ) ; o. mandel , m. greiner , a. widera , t. rom , t. w. hnsch , and i. bloch , phys .
. lett . * 91 * , 010407 ( 2003 ) ; l. e. sadler , j. m. higbie , s. r. leslie , m. vengalottore , and m. s. stamper - kurn , nature * 443 * , 312 ( 2006 ) ; j. larson and j .- p .
martikainen , phys .
rev . a * 80 * , 033605 ( 2009 ) .
lin , r. l. compton , a. r. perry , w. d. phillips , j. v. porto , and i. b. spielman , phys .
lett . * 102 * , 130401 ( 2009 ) ; y .- j .
lin , r. l. compton , k. jimnez - garcia , j. v. porto , and i. b. spielman , nature * 462 * , 628 ( 2009 ) .
it is worth pointing out that even if a confining trap is absent in the actual bo experiments , it is present during the condensation process and as it is turned off before releasing the condensate in the tilted lattice it implies , in the local density approximation , that the chemical potential varies within the extent of the atomic cloud which renders a source for decoherence , see  @xcite .
o. morsch , j. h. mller , m. cristiani , d. ciampini , and e. arimondo , phys .
* 87 * , 140402 ( 2001 ) ; m. cristiani , o. morsch , j. h. mller , d. ciampini , e. arimondo , phys .
a * 65 * , 063612 ( 2002 ) ; j. h. denschlag , j. e. simsarian , h. haffner , c. mckenzie , a. browaeys , d. cho , k. helmerson , s. l. rolston , and w. d. phillips , j. phys .
b : at . mol .
* 35 * , 3095 ( 2002 ) ; m. gustavsson , e. haller , m. j. mark , j. g. danzl , g. rojas - kopeinig , and h .- c .
ngerl , phys .
lett . * 100 * , 080404 ( 2008 ) .
t. ando , j. phys .
japan * 69 * , 1757 ( 2000 ) ; c. l. kane and e. j. male , phys
* 95 * , 226801 ( 2005 ) ; h. min , j. e. hill , n. a. sinitsyn , b. r. sahu , l. kleinman , and a. h. macdonald , phys .
b * 74 * , 165310 ( 2006 ) .                                        v. y. demikhovskii and d. v. khmiisky , j. exp
* 83 * , 340 ( 2006 ) ; s. smirnov , d. bercioux , and m. grifoni , europhys . lett . * 80 * , 27003 ( 2007 ) ; d. v. khmiisky , phys . rev .
b * 79 * , 205401 ( 2009 ) .
k. leo , p. h. bolivar , f. bruggemann , r. schwedler , and k. kohler , solid .
. commun . *
84 * , 943 ( 1992 ) ; j. feldmann , k. leo , j. shah , d. a. b. miller , j. e. cunningham , t. meier , g. vonplessen , a. schulze , p. thomas , and s. schmittrink , phys .
b * 46 * , 7252 ( 1992 ) .
r. morandotti , u. peschel , j. s. aitchison , h. s. eisenberg , and k. silberberg , phys .
lett . * 83 * , 4756 ( 1999 ) ; t. pertsch , p. dannberg , w. elflein , a. brauer , and f. lederer , phys .
lett . * 83 * , 4752 ( 1999 ) .
t.  gericke , p.  wrtz , d.  reitz , t.  langen , and h.  ott , nature physics * 4 * , 949 ( 2008 ) ; n.  gemelke , x.  zhang , ch .- l .  hung , and ch .
chin , nature * 460 * , 995 ( 2009 ) ; w.  s.  bakr , j.  i.  gillen , a.  peng , s.  flling , and m.  greiner , nature * 462 * , 74 ( 2009 ) ; j. f. sherson , c. weitenberg , m. endres , m. cheneau , i. bloch , and s. kuhr , nature * 18 * , august ( 2010 ) .