<ph f="cmex"> </ph><ph f="cmbx">Random walks in a Dirichlet environment</ph>

### Christophe Sabot

Laboratoire de Probabilites de Paris 6, 4 place Jussieu, 75252 Paris Cedex 05 E-mail address : enriquez@ccr.jussieu.fr Unite de Mathematiques Pures et Appliquees, 46, Allee d'Italie F-69364 Lyon Cedex 07 E-mail address : christophe.sabot@umpa.ens-lyon.fr
Abstract: This paper states a law of large numbers for a random walk in a random iid environment on ${\mathbb{Z}}^{d}$  , where the environment follows some Dirichlet distribution. Moreover, we give explicit bounds for the asymptotic velocity of the process and also an asymptotic expansion of this velocity at low disorder.

1 Introduction

After a first breakthrough of Kalikow [3, giving a transience criterion for non-reversible multidimensional Random Walks in Random Environment, Sznitman and Zerner proved, several years later, a law of large numbers in [7, followed by a central limit theorem proved by Sznitman in [6. A generalization to the case of mixing environments was proved afterwards by Comets and Zeitouni in [1(we refer to [8for an overview of the subject). Despite these progresses, many important questions, concerning recurrence or explicit criteria for a ballistic behavior, remain largely open.
Among random walks in random environment, random walks in an iid Dirichlet environment take a special place, since their annealed law coincides with the law of some transition reinforced random walk having an affine reinforcement (see [2). These reinforced walks are defined as follows.
At time 0, we attribute, in a translation invariant way, a weight to each oriented edge of ${\mathbb{Z}}^{d}$  , and each time the walk crosses an edge, the weight of this edge is increased by one. Finally, the walk is a nearest neighbour walk, which chooses, at each time, an outgoing edge with a probability which is proportionnal to its weight.
In this paper, we state a law of large numbers for such random walks, under a simple and explicit condition on the weights. Moreover, we give explicit bounds for the asymptotic velocity of these walks and also an asymptotic expansion of this velocity at low disorder. Low disorder corresponds, in the random environment model, to the case where the law of the transition probabilities is concentrated around its mean value, and, in its reinforcement interpretation, to the case where the initial weights of the transitions are large, so that these weights are not significantly affected during the life of the walk (at least, if the walk is transient).
Let us precise that these walks do not enter the class of walks considered in [3, [4and in several other works, asking the law of the environment to satisfy a uniform ellipticity condition.
This ellipticity hypothesis is usually used in two ways :
_ in the definition of Kalikow's auxiliary Markov chain which involves the expection of the Green function of the walk killed when exiting a given set. The uniform ellipticity is then a comfortable assumption for checking the integrability of this Green function.
_ in the estimates of the drift of Kalikow's auxiliary Markov chain, the ellipticity condition often plays a key role. We overcome this difficulty by using an integration by part formula.
In section 2, we give the definition of random walks in Dirichlet environment, remind their connection with transition reinforced random walks and we present our main results. In section 3, we present an integration by part formula that will be the key analytic tool in the proof of our results. Indeed, in section 4, it is shown how one can take advantage of the special form of the law of the environment, in order to estimate, using the formula of Section 3, the drift of the killed Kalikow's auxiliary walk.
In section 5, we study the integrability of the Green function of the walk which ensures the existence of the original (non killed) Kalikow's auxiliary walk and finish the proof of our first result by applying the law of large numbers of Sznitman and Zerner [7. In section 6, we follow the scheme of [4to get bounds for the asymptotic velocity of the walk, and deduce an expansion of the asymptotic velocity at low disorder.

2 Definitions and statement of the results

We denote by ${T}_{2d}:=\left\{\left({x}_{1},...,{x}_{2d}\right)\in \right]0,1{\right]}^{2d},s.t.,{\sum }_{i=1}^{2d}{x}_{i}=1\right\}$  , and by $\left({e}_{i}{\right)}_{1\le i\le 2d}$  the family of unitary vectors of ${\mathbb{Z}}^{d}$  , defined as follows: $\left({e}_{i}{\right)}_{1\le i\le d}$  is the canonical basis of ${\mathbb{R}}^{d}$  , and for all $j\in \left\{d+1,...,2d\right\}$  ${e}_{j}=-{e}_{j-d}$  .
For all $\stackrel{⃗}{\alpha }:=\left({\alpha }_{1},...,{\alpha }_{2d}\right)\in \right]0,+\infty {\left[}^{2d}$  , we denote by ${\lambda }^{\stackrel{⃗}{\alpha }}$  the Dirichlet probability measure on ${T}_{2d}$  with parameters $\left({\alpha }_{1},...,{\alpha }_{2d}\right)$  i.e. the measure on ${T}_{2d}$  :
$\frac{\Gamma \left({\alpha }_{1}+...+{\alpha }_{2d}\right)}{\Gamma \left({\alpha }_{1}\right)...\Gamma \left({\alpha }_{2d}\right)}{x}_{1}^{{\alpha }_{1}-1}...{x}_{2d}^{{\alpha }_{2d}-1}d{x}_{1}...d{x}_{2d-1}.$  For a unit vector $e$  of ${\mathbb{Z}}^{d}$  , we will sometimes write, for reading conveniences, ${\alpha }_{e}$  for the weight ${\alpha }_{i}$  where $i$  is such that ${e}_{i}=e$  .
Let us now introduce random walks in an iid Dirichlet environment on ${\mathbb{Z}}^{d}$  .
We define an environment as an element $\omega =\left(\omega \left(x\right){\right)}_{x\in {\mathbb{Z}}^{d}}$  where at any vertex $x$  , $\omega \left(x\right):=\left(\omega \left(x,x+{e}_{1}\right),...,\omega \left(x,x+{e}_{2d}\right)\right)$  belongs to ${T}_{2d}$  . We set $\mu :={\otimes }_{x\in {\mathbb{Z}}^{d}}{\lambda }^{\stackrel{⃗}{\alpha }}$  , so that $\mu$  is a probability measure on the environments such that $\left(\omega \left(x\right){\right)}_{x\in {\mathbb{Z}}^{d}}$  are independent random variables of law ${\lambda }^{\stackrel{⃗}{\alpha }}$  .
We denote by ${P}^{\omega }$  the law of the Markov chain in the environment $\omega$  starting at $0$  defined by:
$\forall x\in {\mathbb{Z}}^{d},\forall k\in \mathbb{N},\forall i=1,...,2d,{P}^{\omega }\left({X}_{k+1}=x+{e}_{i}|{X}_{k}=x\right)=\omega \left(x,x+{e}_{i}\right).$  The law of the random walk in random environment (or the so-called annealed measure) is the probability measure ${P}^{\mu }=\int {P}^{\omega }d\mu \left(\omega \right)$  .
In [2, we show that random walks in iid environment have the law of some reinforced random walk. The following proposition states that the case of a Dirichlet environment corresponds to a quite natural law of reinforcement:
Proposition 1. The measure ${P}^{\mu }$  satisfies that ${P}^{\mu }$  -almost everywhere, ${P}^{\mu }\left({X}_{n+1}=x+{e}_{i}|\sigma \left({X}_{k},k\le n\right)\right)=\frac{{\alpha }_{i}+{N}_{i}\left(n,{X}_{n}\right)}{{\sum }_{k=1}^{r}{\alpha }_{k}+{N}_{k}\left(n,{X}_{n}\right)}$  where $\stackrel{⃗}{N}\left(n,x\right)=\left({N}_{i}\left(n,x\right){\right)}_{1\le i\le 2d}$  and ${N}_{i}\left(n,x\right)={\sum }_{l=0}^{n-1}{1}_{\left\{{X}_{l+1}-{X}_{l}={e}_{i},{X}_{l}=x\right\}}$  .
We refer the reader to [2for the proof.

2.1 Bounds of the asymptotic velocity

We can now state our first result:
Theorem 1. Let $\stackrel{⃗}{\alpha }:=\left({\alpha }_{1},...,{\alpha }_{2d}\right)\in \right]0,+\infty {\left[}^{2d}$  , and $\mu ={\otimes }_{x\in {\mathbb{Z}}^{d}}{\lambda }^{\stackrel{⃗}{\alpha }}$  a probability measure on the environment. Let us assume that there exists $i\in \left\{1,...,2d\right\}$  such that ${\alpha }_{{e}_{i}}>1+{\alpha }_{-{e}_{i}}$  .
The process ${X}_{n}$  is transient under ${P}^{\mu }$  , and $\exists v\in {\mathbb{R}}^{d}\\left\{0\right\},\text{such that}{P}^{\mu }\left(\frac{{X}_{n}}{n}{\to }_{n\to \infty }v\right)=1.$  Moreover, for all $i\in \left\{1,...,d\right\}$  , $\frac{1}{\left({\sum }_{k=1}^{2d}{\alpha }_{k}\right)-1}\left(min\left(0,{\alpha }_{{e}_{i}}-1\right)-{\alpha }_{-{e}_{i}}\right)\le v.{e}_{i}\le \frac{1}{\left({\sum }_{k=1}^{2d}{\alpha }_{k}\right)-1}\left({\alpha }_{{e}_{i}}-min\left(0,{\alpha }_{-{e}_{i}}-1\right)\right)$
Remark 1: The assumption on ${\alpha }_{i}$  ensures that the set ${}^{d}{\prod }_{i=1}\left[min\left(0,{\alpha }_{{e}_{i}}-1\right)-{\alpha }_{-{e}_{i}},{\alpha }_{{e}_{i}}-min\left(0,{\alpha }_{-{e}_{i}}-1\right)\right]$  does not contain 0. It is a key ingredient in the check of Kalikow's transience condition.
Remark 2: When the ${\alpha }_{i}$  's are large, $v$  becomes close to the vector $\frac{1}{{\sum }_{k=1}^{2d}{\alpha }_{k}}{\sum }_{i=1}^{d}\left({\alpha }_{{e}_{i}}-{\alpha }_{-{e}_{i}}\right){e}_{i}$  .
This is not surprising if one thinks at the corresponding reinforced walk: the initial weights of the transitions are large enough so that they are not significantly affected during the life of the walk, and the law of the walk becomes close to the law of the Markov chain with probability transition $\frac{{\alpha }_{i}}{{\sum }_{k=1}^{2d}{\alpha }_{k}}$  in the direction ${e}_{i}$  .
Remark 3: In dimension 1, the condition of theorem 1 is actually optimal. Indeed from [5, we know that the asymptotic velocity is not null if and only if either ${E}^{\mu }\left[\frac{\omega \left(0,{e}_{1}\right)}{\omega \left(0,-{e}_{1}\right)}\right]>1$  or ${E}^{\mu }\left[\frac{\omega \left(0,-{e}_{1}\right)}{\omega \left(0,{e}_{1}\right)}\right]>1$  , which corresponds exactly to ${\alpha }_{{e}_{1}}>1+{\alpha }_{-{e}_{1}}$  or ${\alpha }_{-{e}_{1}}>1+{\alpha }_{{e}_{1}}$  .

2.2 Expansion of the velocity in the limit of large parameters

We turn now to the second result of the paper, which gives the asymptotic velocity of the walk in the limit of large parameters ${\alpha }_{k}$  . Let us remind that, in the limit of large parameters ${\alpha }_{k}$  , the environment is concentrated around its mean value.
Let us fix some notations. We consider some fixed transition probabilities $\left({m}_{i}\right)\in {T}_{2d},$  and a parameter $\gamma >0$  (aimed to tend to $\infty$  ). We consider the weights ${\alpha }_{k}=\gamma {m}_{k},$  so that the expectation of the transition probability, ${E}_{\mu }\left(\omega \left(x,x+{e}_{i}\right)\right),$  is independent of $\gamma$  and equal to ${m}_{i}$  .
The mean environment $\left({m}_{i}\right)$  defines the transition probabilities of an homogeneous walk on ${\mathbb{Z}}^{d}$  , which is ballistic with asymptotic velocity ${d}_{m}={\sum }_{k=1}^{2d}{m}_{k}{e}_{k},$  when the mean drift ${d}_{m}$  is not null. We denote by ${G}^{m}$  its Green function.
The following result gives an estimate in $O\left(\frac{1}{{\gamma }^{2}}\right)$  of the asymptotic velocity (in section 6 we give explicit bounds for this estimate).
Theorem 2. Assume ${d}_{m}\ne 0$  .
For $\gamma$  large enough, Theorem 1 applies, i.e. there exists $v\ne 0$  such that ${lim}_{n\to \infty }\frac{{X}_{n}}{n}=v,{P}^{\mu }a.s.$  .
Moreover, when $\gamma$  is large, we have the following expansion for $v$  :
$v={d}_{m}-\frac{{d}_{m}}{\gamma }\left({G}^{m}\left(0,0\right)-1\right)+O\left(\frac{1}{{\gamma }^{2}}\right).$
Remark 1: Surprisingly, the second order of the expansion is colinear to the mean drift ${d}_{m}$  . We see that $\left({G}^{m}\left(0,0\right)-1\right)>0$  , which means that there is a slowdown effect, since the second order term is directed in the opposite direction to the mean drift.
Remark 2: In [4, the second author gave an expansion of the asymptotic velocity in the case of a uniformly elliptic environment. In this work, several of the estimates relied strongly on the ellipticity condition, so that the proofs of [4have here to be modified. Nevertheless, if we apply the formula of [4to this case, we get the same expansion (many simplifications occur due to the particular expression of the covariance matrix). It is not surprising that the speed-up effect obtained in some cases of [4is not observed in the case of a Dirichlet environment. The example of [4, section 2, was based, indeed, on some correlation between the transition probabilities in orthogonal directions. Here, there is a kind of independence of the transition probabilities in each direction, in the following sense: under $\mu ={\otimes }_{x\in {\mathbb{Z}}^{d}}{\lambda }^{\stackrel{⃗}{\alpha }}$  , the law of $\omega \left(z,z+{e}_{i}\right)$  is independent of the law of $\left(\frac{\omega \left(z,z+{e}_{k}\right)}{1-\omega \left(z,z+{e}_{i}\right)}{\right)}_{k\ne i}$  .
Remark 3: The Green function ${G}^{m}\left(0,0\right)$  has the following explicit Fourier expression ${G}^{m}\left(0,0\right)=\frac{1}{\left(2\pi {\right)}^{d}}{\int }_{\left[0,2\pi {\right]}^{d}}\frac{1}{1-2{\sum }_{i=1}^{d}\sqrt{{m}_{{e}_{i}}{m}_{-{e}_{i}}}cos\left({\theta }_{i}\right)}d{\theta }_{1}\cdots d{\theta }_{d}.$  (we refer to Step 2 of the proof of Proposition 3).

3 An integration by part formula

In this section, we present an integration by part formula on ${T}_{2d}$  that will appear to be the key analytic tool in the estimation of the drift of Kalikow's auxiliary walk.
Lemma 1. For all $\stackrel{⃗}{\alpha }\in \right]0,+\infty {\left[}^{2d}$  , and all differentiable function $f$  on ${\mathbb{R}}^{2d}$  , ${\int }_{{T}_{2d}}fd{\lambda }^{\stackrel{⃗}{\alpha }}=\frac{{\alpha }_{1}+...+{\alpha }_{d}}{{\alpha }_{1}}{\int }_{{T}_{2d}}{x}_{1}.fd{\lambda }^{\stackrel{⃗}{\alpha }}+\frac{1}{{\alpha }_{1}}{\int }_{{T}_{2d}}{x}_{1}.\left(\left({\sum }_{k=1}^{2d}{x}_{k}\frac{\partial f}{\partial {x}_{k}}\right)-\frac{\partial f}{\partial {x}_{1}}\right)d{\lambda }^{\stackrel{⃗}{\alpha }}.$
Proof: We use the well known identity between the Dirichlet law ${\lambda }^{\stackrel{⃗}{\alpha }}$  and the law of the vector $\left(\frac{{Z}_{1}}{{\sum }_{i=1}^{2d}{Z}_{i}},...,\frac{{Z}_{2d}}{{\sum }_{i=1}^{2d}{Z}_{i}}\right)$  where the random variables ${Z}_{i}$  are independent variables following the gamma distribution $\Gamma \left({\alpha }_{i},1\right)$  of density $\frac{1}{\Gamma \left({\alpha }_{i}\right)}{z}^{{\alpha }_{i}-1}{e}^{-z}$  on ${\mathbb{R}}_{+}$  .
This identity implies ${\int }_{{T}_{2d}}fd{\lambda }^{\stackrel{⃗}{\alpha }}=$  $\frac{1}{\Gamma \left({\alpha }_{1}\right)...\Gamma \left({\alpha }_{2d}\right)}{\int }_{{\mathbb{R}}_{+}^{2d}}f\left(\frac{{z}_{1}}{{\sum }_{i=1}^{2d}{z}_{i}},...,\frac{{z}_{2d}}{{\sum }_{i=1}^{2d}{z}_{i}}\right){e}^{-{\sum }_{i=1}^{2d}{z}_{i}}{z}_{1}^{{\alpha }_{1}-1}...{z}_{2d}^{{\alpha }_{2d}-1}d{z}_{1}...d{z}_{2d}.$  Integrating by part with respect to ${z}_{1}$  , we get ${\int }_{{T}_{2d}}fd{\lambda }^{\stackrel{⃗}{\alpha }}=\frac{1}{\Gamma \left({\alpha }_{1}+1\right)...\Gamma \left({\alpha }_{2d}\right)}{\int }_{{\mathbb{R}}_{+}^{2d}}\left(\stackrel{~}{f}-\frac{\partial \stackrel{~}{f}}{\partial {z}_{1}}\right){e}^{-{\sum }_{i=1}^{2d}{z}_{i}}{z}_{1}^{{\alpha }_{1}}...{z}_{2d}^{{\alpha }_{2d}-1}d{z}_{1}...d{z}_{2d}$  where $\stackrel{~}{f}\left({z}_{1},...,{z}_{2d}\right):=f\left(\frac{{z}_{1}}{{\sum }_{i=1}^{2d}{z}_{i}},...,\frac{{z}_{2d}}{{\sum }_{i=1}^{2d}{z}_{i}}\right)$  .
Now, we decompose this last integral into the $\stackrel{~}{f}$  -part and the $\frac{\partial \stackrel{~}{f}}{\partial {z}_{1}}$  -part.
Using, in the reverse sense, the “Gamma” interpretation of the Dirichlet law ${\lambda }^{\left({\alpha }_{1}+1,{\alpha }_{2},...,{\alpha }_{2d}\right)}$  , the $\stackrel{~}{f}$  -part becomes $\frac{\Gamma \left({\alpha }_{1}+...+{\alpha }_{2d}+1\right)}{\Gamma \left({\alpha }_{1}+1\right)...\Gamma \left({\alpha }_{2d}\right)}{\int }_{{T}_{2d}}f.{x}_{1}^{{\alpha }_{1}}...{x}_{2d}^{{\alpha }_{2d}-1}d{x}_{1}...d{x}_{2d-1}=\frac{{\alpha }_{1}+...+{\alpha }_{2d}}{{\alpha }_{1}}{\int }_{{T}_{2d}}{x}_{1}.fd{\lambda }^{\stackrel{⃗}{\alpha }}.$  Now, the $\frac{\partial \stackrel{~}{f}}{\partial {z}_{1}}$  -part writes $-\frac{1}{{\alpha }_{1}\Gamma \left({\alpha }_{1}\right)...\Gamma \left({\alpha }_{2d}\right)}{\int }_{{\mathbb{R}}_{+}^{d}}\left({z}_{1}.\frac{\partial \stackrel{~}{f}}{\partial {z}_{1}}\right){e}^{-{\sum }_{i=1}^{2d}{z}_{i}}{z}_{1}^{{\alpha }_{1}-1}...{z}_{2d}^{{\alpha }_{2d}-1}d{z}_{1}...d{z}_{2d}$  and
 $\begin{array}{cc}{z}_{1}.\frac{\partial \stackrel{~}{f}}{\partial {z}_{1}}& =\left(\frac{{z}_{1}}{{\sum }_{i=1}^{2d}{z}_{i}}-\frac{{z}_{1}^{2}}{\left({\sum }_{i=1}^{2d}{z}_{i}{\right)}^{2}}\right){\stackrel{~}{f}}_{1}-\frac{{z}_{1}{z}_{2}}{\left({\sum }_{i=1}^{2d}{z}_{i}{\right)}^{2}}{\stackrel{~}{f}}_{2}-...-\frac{{z}_{1}{z}_{2d}}{\left({\sum }_{i=1}^{2d}{z}_{i}{\right)}^{2}}{\stackrel{~}{f}}_{d}\end{array}$
 $\begin{array}{cc}& =\left(\frac{{z}_{1}}{{\sum }_{i=1}^{2d}{z}_{i}}\right)\left({\stackrel{~}{f}}_{1}-{\sum }_{k=1}^{2d}\left(\frac{{z}_{k}}{{\sum }_{i=1}^{2d}{z}_{i}}\right){\stackrel{~}{f}}_{k}\right)\end{array}$
where ${\stackrel{~}{f}}_{k}\left({z}_{1},..,{z}_{2d}\right)=\frac{\partial f}{\partial {x}_{k}}\left(\frac{{z}_{1}}{{\sum }_{i=1}^{2d}{z}_{i}},...,\frac{{z}_{2d}}{{\sum }_{i=1}^{2d}{z}_{i}}\right)$  . The “Gamma” interpretation of the Dirichlet law ${\lambda }^{\left({\alpha }_{1},...,{\alpha }_{2d}\right)}$  (used for the third time) allows to conclude.

4 Kalikow's auxiliary walk

We remind here the generalization of Kalikow's auxiliary walk (see [3) which was already presented in [4.
Let $U$  be a connected subset of ${\mathbb{Z}}^{d}$  , and $\delta \in \right]0,1\right]$  . We denote by $\partial U$  the boundary set of $U$  , i.e.
$\partial U:=\left\{z\in {\mathbb{Z}}^{d}\U,\exists x\in U,|z-x|=1\right\}$  .
For all $z\in U$  and ${z}^{\prime }\in U\cup \partial U$  , and for all environment $\omega$  , we introduce the Green function of the random walk under the environment $\omega$  killed at rate $\delta$  and at the boundary of $U$  :
${G}_{U,\delta }^{\omega }\left(z,{z}^{\prime }\right)={E}_{z}^{\omega }\left({\sum }_{k=0}^{{T}_{U}}{\delta }^{k}{1}_{{X}_{k}={z}^{\prime }}\right)$  where ${T}_{U}=inf\left\{k,{X}_{k}\in {\mathbb{Z}}^{d}\U\right\}$  .
In the sequel, we will drop the subscript $\delta$  when $\delta =1$  , and we will write ${G}_{U}^{\omega }\left(z,{z}^{\prime }\right)$  instead of ${G}_{U,1}^{\omega }\left(z,{z}^{\prime }\right)$  .
We introduce now the generalized Kalikow's transition probabilities (originally, Kalikow's transition probabilities were introduced in the case $\delta =1$  ):
${\stackrel{^}{\omega }}_{U,\delta ,{z}_{0}}\left(z,z+{e}_{i}\right)=\frac{{E}_{\mu }\left[{G}_{U,\delta }^{\omega }\left({z}_{0},z\right)\omega \left(z,z+{e}_{i}\right)\right]}{{E}_{\mu }\left[{G}_{U,\delta }^{\omega }\left({z}_{0},z\right)\right]}.$  In order to give bounds for these transition probabilities, we will be lead to apply the integration by part formula of the previous section to the functions ${G}_{U,\delta }^{\omega }\left(x,y\right)$  , viewed as functions of the environment $\omega$  .
For this purpose, we need the following lemma which gives the expression of the derivatives of these functions:
Lemma 2. For all connected subset $U$  of ${\mathbb{Z}}^{d}$  , for all ${x}_{1},{x}_{2},{x}_{3},{x}_{4}\in U$  , and for all $\delta \in \right]0,1\left[$  , $\frac{\partial {G}_{U,\delta }^{\omega }\left({x}_{1},{x}_{4}\right)}{\partial \left(\omega \left({x}_{2},{x}_{3}\right)\right)}=\delta {G}_{U,\delta }^{\omega }\left({x}_{1},{x}_{2}\right){G}_{U,\delta }^{\omega }\left({x}_{3},{x}_{4}\right)$
Proof: Let us define the transition matrix ${\Omega }^{U}:=\left({\omega }^{U}\left(x,y\right){\right)}_{\left(x,y\right)\in \left(U\cup \partial U{\right)}^{2}}$  , by ${\omega }^{U}\left(x,y\right)=\omega \left(x,y\right)$  if $x\in U$  , and ${\omega }^{U}\left(x,y\right)=0$  if $x\in \partial U$  . We have ${G}_{U,\delta }^{\omega }\left({x}_{1},{x}_{4}\right)={\sum }_{n\ge 0}{\delta }^{n}\left({\Omega }^{U}{\right)}_{\left({x}_{1},{x}_{4}\right)}^{n}$  and $\frac{\partial \left({\Omega }^{U}{\right)}_{\left({x}_{1},{x}_{4}\right)}^{n}}{\partial \left(\omega \left({x}_{2},{x}_{3}\right)\right)}={\sum }_{{k}_{1}+{k}_{2}=n-1}\left({\Omega }^{U}{\right)}_{\left({x}_{1},{x}_{2}\right)}^{{k}_{1}}\left({\Omega }^{U}{\right)}_{\left({x}_{3},{x}_{4}\right)}^{{k}_{2}}$  so that, taking the derivatives term by term in the sum defining ${G}_{U,\delta }^{\omega }\left({x}_{1},{x}_{4}\right)$  , we obtain the result. We turn now to the estimation of the transition probabilities:
Proposition 2. For all connected subset $U$  of ${\mathbb{Z}}^{d}$  , for all ${z}_{0},z\in U$  , for all $\delta \in \right]0,1\left[$  and all $i=1,...,2d$  , $\bullet$  if $\left({\sum }_{k=1}^{2d}{\alpha }_{k}\right)>1$  , then $min\left(0,\frac{{\alpha }_{i}-1}{\left({\sum }_{k=1}^{2d}{\alpha }_{k}\right)-1}\right)\le {\stackrel{^}{\omega }}_{U,\delta ,{z}_{0}}\left(z,z+{e}_{i}\right)\le \frac{{\alpha }_{i}}{\left({\sum }_{k=1}^{2d}{\alpha }_{k}\right)-1}$  $\bullet$  if $\left({\sum }_{k=1}^{2d}{\alpha }_{k}\right)<1$  , then $0\le {\stackrel{^}{\omega }}_{U,\delta ,{z}_{0}}\left(z,z+{e}_{i}\right)\le \frac{1-{\alpha }_{i}}{1-\left({\sum }_{k=1}^{2d}{\alpha }_{k}\right)}$
Proof: For the clarity of notations we give the proof for $i=1$  .
Lemma 2 yields $\frac{\partial {G}_{U,\delta }^{\omega }\left({z}_{0},z\right)}{\partial \left(\omega \left(z,z+{e}_{i}\right)\right)}=\delta {G}_{U,\delta }^{\omega }\left({z}_{0},z\right){G}_{U,\delta }^{\omega }\left(z+{e}_{i},z\right).$  We now apply Lemma 1 with $f={G}_{U,\delta }^{\omega }\left({z}_{0},z\right)$  , viewed as a function of the only variables ${x}_{i}:=\omega \left(z,z+{e}_{i}\right)$  for $i=1,...,2d$  , and we get ${E}_{\mu }\left[{G}_{U,\delta }^{\omega }\left({z}_{0},z\right)\right]=\frac{{\alpha }_{1}+...+{\alpha }_{2d}}{{\alpha }_{1}}{E}_{\mu }\left[{G}_{U,\delta }^{\omega }\left({z}_{0},z\right)\omega \left(z,z+{e}_{1}\right)\right]$  $+\frac{1}{{\alpha }_{1}}{E}_{\mu }\left[\omega \left(z,z+{e}_{1}\right).{G}_{U,\delta }^{\omega }\left({z}_{0},z\right)\left(\delta {\sum }_{k=1}^{2d}\omega \left(z,z+{e}_{k}\right){G}_{U,\delta }^{\omega }\left(z+{e}_{k},z\right)-\delta {G}_{U,\delta }^{\omega }\left(z+{e}_{1},z\right)\right)\right]\left(1\right)$  We recall that $\delta {\sum }_{k=1}^{2d}\omega \left(z,z+{e}_{k}\right){G}_{U,\delta }^{\omega }\left(z+{e}_{k},z\right)={G}_{U,\delta }^{\omega }\left(z,z\right)-1$  so that the second term in the right side of (1) writes $\frac{1}{{\alpha }_{1}}{E}_{\mu }\left[\omega \left(z,z+{e}_{1}\right).{G}_{U,\delta }^{\omega }\left({z}_{0},z\right)\left({G}_{U,\delta }^{\omega }\left(z,z\right)-1-\delta {G}_{U,\delta }^{\omega }\left(z+{e}_{1},z\right)\right)\right]$  so that we get
 $\begin{array}{cc}{E}_{\mu }\left[{G}_{U,\delta }^{\omega }\left({z}_{0},z\right)\right]=& \frac{{\alpha }_{1}+...+{\alpha }_{2d}}{{\alpha }_{1}}{E}_{\mu }\left[{G}_{U,\delta }^{\omega }\left({z}_{0},z\right)\omega \left(z,z+{e}_{1}\right)\right]\end{array}$
 $\begin{array}{cc}& +\frac{1}{{\alpha }_{1}}{E}_{\mu }\left[\omega \left(z,z+{e}_{1}\right).{G}_{U,\delta }^{\omega }\left({z}_{0},z\right)\left({G}_{U,\delta }^{\omega }\left(z,z\right)-1-\delta {G}_{U,\delta }^{\omega }\left(z+{e}_{1},z\right)\right)\right]\end{array}$
and for the ratio ${\stackrel{^}{\omega }}_{U,\delta ,{z}_{0}}\left(z,z+{e}_{1}\right)=\frac{{E}_{\mu }\left[\omega \left(z,z+{e}_{1}\right){G}_{U,\delta }^{\omega }\left({z}_{0},z\right)\right]}{{E}_{\mu }\left[{G}_{U,\delta }^{\omega }\left({z}_{0},z\right)\right]}$  , ${\stackrel{^}{\omega }}_{U,\delta ,{z}_{0}}\left(z,z+{e}_{1}\right)=\frac{{\alpha }_{1}}{\left({\sum }_{k=1}^{2d}{\alpha }_{k}\right)-1}$  $+\frac{1}{\left({\sum }_{k=1}^{2d}{\alpha }_{k}\right)-1}\frac{{E}_{\mu }\left[\omega \left(z,z+{e}_{1}\right).{G}_{U,\delta }^{\omega }\left({z}_{0},z\right)\left({G}_{U,\delta }^{\omega }\left(z,z\right)-\delta {G}_{U,\delta }^{\omega }\left(z+{e}_{1},z\right)\right)\right]}{{E}_{\mu }\left[{G}_{U,\delta }^{\omega }\left({z}_{0},z\right)\right]}\left(2\right)$  But, ${\sum }_{k=1}^{2d}\omega \left(z,z+{e}_{k}\right)\left({G}_{U,\delta }^{\omega }\left(z,z\right)-\delta {G}_{U,\delta }^{\omega }\left(z+{e}_{k},z\right)\right)=1$  and therefore, for all $k=1,...,2d$  , $0\le {G}_{U,\delta }^{\omega }\left(z,z\right)-\delta {G}_{U,\delta }^{\omega }\left(z+{e}_{k},z\right)\le \frac{1}{\omega \left(z,z+{e}_{1}\right)}.$  These inequalities allow to bound the ratio in the second term of the right side of (2), between 0 and 1, and this finishes the proof.

5 Proof of Theorem 1

We gather now all the ingredients of the proof of Theorem 1. We want to apply Sznitman and Zerner's law of large numbers [7. From a careful reading of the proof of this law of large numbers, we can see that the only conditions that need to be fullfilled, are the integrability of the Green function ${G}_{U}^{\omega }\left({z}_{0},{z}_{0}\right)$  for all bounded $U$  , and Kalikow's condition.
The integrability of the Green function is proved in the following lemma:
Lemma 3. If there exists $i\in \left\{1,...,2d\right\}$  , such that ${\alpha }_{i}>1$  , then for all connected subset $U$  of ${\mathbb{Z}}^{d}$  and all ${z}_{0}\in U$  , ${E}_{\mu }\left[{G}_{U}^{\omega }\left({z}_{0},{z}_{0}\right)\right]$  is finite.
Proof: For the clarity of notations, we suppose that ${\alpha }_{1}>1$  .
Define now by $N$  the least integer such that ${z}_{0}+N{e}_{1}$  belongs to $\partial U$  .
We have the following lower bound for the probability $P\left(\omega ,{z}_{0},U\right)$  to reach $\partial U$  from ${z}_{0}$  without returning to ${z}_{0}$  :
$P\left(\omega ,{z}_{0},U\right){\ge }^{N-1}{\prod }_{k=0}\omega \left({z}_{0}+k{e}_{1},{z}_{0}+\left(k+1\right){e}_{1}\right).$  The number of returns to ${z}_{0}$  before hitting $\partial U$  , being a geometric variable whose parameter is precisely $P\left(\omega ,{z}_{0},U\right)$  , its expectation ${G}_{U}^{\omega }\left({z}_{0},{z}_{0}\right)$  is equal to $\frac{1}{P\left(\omega ,{z}_{0},U\right)}$  .
We are now lead to examine the integrability of ${E}_{\mu }\left[{\left({\prod }_{k=0}^{N-1}\omega \left({z}_{0}+k{e}_{1},{z}_{0}+\left(k+1\right){e}_{1}\right)\right)}^{-1}\right]$  which is equal to $\left({\int }_{{T}_{2d}}\frac{1}{{x}_{1}}d{\lambda }^{\stackrel{⃗}{\alpha }}{\right)}^{N}$  which is finite since ${\alpha }_{1}>1$  . We now have to check Kalikow's condition.
We notice first that, under the assumption of Theorem 1, Lemma 3 applies and Kalikow's auxiliary walk is well defined. Then, the monotone convergence theorem allows to make $\delta$  converge to 1 in the inequalities of Proposition 2.
We then deduce that the drift of Kalikow's walk belongs to ${\frac{1}{\left({\sum }_{k=1}^{2d}{\alpha }_{k}\right)-1}}^{d}{\prod }_{i=1}\left[min\left(0,{\alpha }_{{e}_{i}}-1\right)-{\alpha }_{-{e}_{i}},{\alpha }_{{e}_{i}}-min\left(0,{\alpha }_{-{e}_{i}}-1\right)\right]$  which does not contain 0, under the assumption of Theorem 1. This proves Kalikow's transience condition.
In order to estimate the asymptotic velocity of the process, we apply directly Proposition 3.2 of [4which makes the link between $v$  and the drift of Kalikow's walk.
Remark: in Lemma 3, we only got a sufficient condition for the integrability of the Green function to hold. A better result about this question would not have ameliorated the statement of Theorem 1 as far as our check of Kalikow's condition requires a stronger assumption.

6 Proof of Theorem 2

Theorem 2 of section 2 is actually a consequence of a more precise result, where the $“O”$  in the expansion is replaced by an explicit upper bound.
Let us fix some notations: we set $\gamma ={\sum }_{i=1}^{2d}{\alpha }_{i},$  and ${m}_{i}={m}_{{e}_{i}}=\frac{{\alpha }_{i}}{\gamma }={E}^{{\lambda }^{\left(\alpha \right)}}\left(\omega \left({e}_{i}\right)\right).$  When $\gamma$  is large, the environment $\left(\omega \left(x,{e}_{i}\right)\right)$  tends to concentrate around its mean $\left({m}_{i}\right)$  , what can be seen from the expression of the correlations $Co{v}_{\mu }\left(\omega \left(x,x+{e}_{i}\right),\omega \left(x,x+{e}_{j}\right)\right)=\left\{\begin{array}{c}-\frac{{m}_{i}{m}_{j}}{\gamma +1},ifi\ne j\\ \frac{{m}_{i}\left(1-{\sum }_{k\ne i}{m}_{k}\right)}{\gamma +1},ifi=j,\end{array}$  The mean environment $\left({m}_{i}\right)$  defines the transition probabilities of an homogeneous walk on ${\mathbb{Z}}^{d}$  , and we define ${k}_{m}=2{\sum }_{i=1}^{d}\sqrt{{m}_{{e}_{i}}{m}_{-{e}_{i}}},$  so that $1-{k}_{m}={\sum }_{i=1}^{d}\left(\sqrt{{m}_{{e}_{i}}}-\sqrt{{m}_{-{e}_{i}}}{\right)}^{2},$  measures the non-symmetry of the walk. When ${k}_{m}<1$  , this walk is ballistic with asymptotic velocity ${d}_{m}={\sum }_{i=1}^{2d}{m}_{i}{e}_{i},$  and we denote by ${G}^{m}\left(\cdot ,\cdot \right)$  its Green function. Let us define ${\eta }_{m}=\frac{{max}_{i}\sqrt{\frac{{m}_{{e}_{i}}}{{m}_{-{e}_{i}}}}}{1-{k}_{m}}.$
Proposition 3. Assume we are in the condition of application of Theorem 1, and that $\frac{2d}{\gamma }{\eta }_{m}\le 1,$  then we have the following estimate $|v-{d}_{m}\left(1-\frac{1}{\gamma -1}\left({G}^{m}\left(0,0\right)-1\right)\right)|\le 16{\left(\frac{d}{\gamma }\right)}^{2}\frac{{\eta }_{m}^{2}}{1-\frac{2d}{\gamma }{\eta }_{m}}.$
Proof: Considering the domain $U={\mathbb{Z}}^{d}$  , a killing parameter $\delta <1$  and ${z}_{0}=0$  , we get from formula (2) ${\stackrel{^}{\omega }}_{\delta }\left(z,z+{e}_{i}\right)={m}_{i}+\frac{{m}_{i}}{\gamma -1}-\frac{1}{\gamma -1}\frac{{E}_{\mu }\left[{G}_{\delta }^{\omega }\left(0,z\right)\omega \left(z,z+{e}_{i}\right)\left({G}_{\delta }^{\omega }\left(z,z\right)-\delta {G}_{\delta }^{\omega }\left(z+{e}_{i},z\right)\right)\right]}{{E}_{\mu }\left[{G}_{\delta }\left(0,z\right)\right]}.$  In the sequel, we will sometimes forget the superscript $\omega$  in ${G}_{\delta }^{\omega }$  , when there will be no ambiguity.
Let us introduce a new probability on the environments $\stackrel{~}{\mu }\left(d\omega \right)$  given by $\stackrel{~}{\mu }\left(d\omega \right)=\frac{{G}_{\delta }^{\omega }\left(0,z\right)}{{E}_{\mu }\left({G}_{\delta }^{\omega }\left(0,z\right)\right)}\mu \left(d\omega \right).$  We see that $\frac{{E}_{\mu }\left[{G}_{\delta }\left(0,z\right)\omega \left(z,z+{e}_{i}\right)\left({G}_{\delta }\left(z,z\right)-\delta {G}_{\delta }\left(z+{e}_{i},z\right)\right)\right]}{{E}_{\mu }\left[{G}_{\delta }\left(0,z\right)\right]}={E}_{\stackrel{~}{\mu }}\left[\left({G}_{\delta }\left(z,z\right)-\delta {G}_{\delta }\left(z+{e}_{i},z\right)\right)\omega \left(z,z+{e}_{i}\right)\right].$  We proceed as in [4, and apply Kalikow's formula (cf. the generalized version in [4, Proposition 3.1) to the measure $\stackrel{~}{\mu }$  .
It means that we have ${E}_{\stackrel{~}{\mu }}\left[{G}_{\delta }^{\omega }\left(z,z\right)\omega \left(z,z+{e}_{i}\right)\right]={G}_{\delta }^{{\stackrel{~}{\omega }}^{z}}\left(z,z\right){\stackrel{~}{\omega }}^{z}\left(z,z+{e}_{i}\right),$  where ${\stackrel{~}{\omega }}^{z}$  is the auxiliary transition probability given by ${\stackrel{~}{\omega }}^{z}\left(y,y+{e}_{j}\right)=\frac{{E}_{\stackrel{~}{\mu }}\left[{G}_{\delta }^{\omega }\left(z,y\right)\omega \left(y,y+{e}_{j}\right)\right]}{{E}_{\stackrel{~}{\mu }}\left[{G}_{\delta }^{\omega }\left(z,y\right)\right]}.$  Similarly, ${E}_{\stackrel{~}{\mu }}\left[{G}_{\delta }^{\omega }\left(z+{e}_{i},z\right)\omega \left(z,z+{e}_{i}\right)\right]={G}_{\delta }^{{\stackrel{~}{\omega }}^{z+{e}_{i}}}\left(z+{e}_{i},z\right){\stackrel{~}{\omega }}^{z+{e}_{i}}\left(z,z+{e}_{i}\right),$  where ${\stackrel{~}{\omega }}^{z+{e}_{i}}$  is the auxiliary transition probability given by ${\stackrel{~}{\omega }}^{z+{e}_{i}}\left(y,y+{e}_{j}\right)=\frac{{E}_{\stackrel{~}{\mu }}\left[{G}_{\delta }^{\omega }\left(z+{e}_{i},y\right)\omega \left(y,y+{e}_{j}\right)\right]}{{E}_{\stackrel{~}{\mu }}\left[{G}_{\delta }^{\omega }\left(z+{e}_{i},y\right)\right]}.$  Step 1: We want to estimate the transition probabilities ${\stackrel{~}{\omega }}^{z}$  and ${\stackrel{~}{\omega }}^{z+{e}_{i}}$  .
Lemma 2 yields $\left(\frac{\partial }{\partial \omega \left(y,y+{e}_{k}\right)}-\frac{\partial }{\partial \omega \left(y,y+{e}_{j}\right)}\right){G}_{\delta }^{\omega }\left(\cdot ,z\right)=\delta {G}_{\delta }^{\omega }\left(\cdot ,y\right)\left({G}_{\delta }^{\omega }\left(y+{e}_{k},z\right)-{G}_{\delta }^{\omega }\left(y+{e}_{j},z\right)\right),$  moreover ${\sum }_{k=1}^{2d}\omega \left(y,y+{e}_{k}\right)\left({G}_{\delta }^{\omega }\left(y,z\right)-\delta {G}_{\delta }^{\omega }\left(y+{e}_{k},z\right)\right)={1}_{y=z}.$  Using the integration by part formula given in Lemma 1, we get
 $\begin{array}{cc}{m}_{{e}_{j}}{E}_{\mu }\left[{G}_{\delta }\left(0,z\right){G}_{\delta }\left(z,y\right)\right]=& {E}_{\mu }\left[{G}_{\delta }\left(0,z\right){G}_{\delta }\left(z,y\right)\omega \left(y,y+{e}_{j}\right)\right]\end{array}$
 $\begin{array}{cc}& \end{array}$
 $\begin{array}{cc}& +\frac{1}{\gamma }{E}_{\mu }\left[{G}_{\delta }\left(0,y\right)\left({G}_{\delta }\left(y,z\right)-\delta {G}_{\delta }\left(y+{e}_{j},z\right)-{1}_{y=z}\right){G}_{\delta }\left(z,y\right)\omega \left(y,y+{e}_{j}\right)\right]\end{array}$
 $\begin{array}{cc}& \end{array}$
 $\begin{array}{cc}& +\frac{1}{\gamma }{E}_{\mu }\left[{G}_{\delta }\left(0,z\right){G}_{\delta }\left(z,y\right)\left({G}_{\delta }\left(y,y\right)-\delta {G}_{\delta }\left(y+{e}_{j},y\right)-1\right)\omega \left(y,y+{e}_{j}\right)\right]\end{array}$
But we have $0\le \omega \left(y,y+{e}_{j}\right)\left({G}_{\delta }\left(y,y\right)-\delta {G}_{\delta }\left(y+{e}_{j},y\right)\right)\le 1,\left(3\right)$  and if $y\ne z$  $|{G}_{\delta }\left(0,y\right)\omega \left(y,y+{e}_{j}\right)\left({G}_{\delta }\left(y,z\right)-\delta {G}_{\delta }\left(y+{e}_{j},z\right)\right)|\le \left(2d-1\right){G}_{\delta }\left(0,z\right).\left(4\right)$  Indeed, for all $k=1,\cdots ,2d$  , we have ${G}_{\delta }^{\omega }\left(y+{e}_{k},z\right)\ge {E}_{y+{e}_{k}}^{\omega }\left[{\delta }^{{T}_{y}}\right]{G}_{\delta }^{\omega }\left(y,z\right),$  where ${T}_{y}$  is the hitting time of $y$  (equal to infinity if the random walk never hits $y$  ).
Since $\frac{1}{1-\delta {\sum }_{k}\omega \left(y,y+{e}_{k}\right){E}_{y+{e}_{k}}^{\omega }\left[{\delta }^{{T}_{y}}\right]}={G}_{\delta }^{\omega }\left(y,y\right),$  we get $\omega \left(y,y+{e}_{k}\right)\left({G}_{\delta }^{\omega }\left(y,z\right)-\delta {G}_{\delta }^{\omega }\left(y+{e}_{k},z\right)\right)\le \frac{{G}_{\delta }^{\omega }\left(y,z\right)}{{G}_{\delta }^{\omega }\left(y,y\right)}.$  But, we also have ${\sum }_{k=1}^{2d}\omega \left(y,y+{e}_{k}\right)\left({G}_{\delta }^{\omega }\left(y,z\right)-\delta {G}_{\delta }^{\omega }\left(y+{e}_{k},z\right)\right)=0,ify\ne z,$  so that we have $|\omega \left(y,y+{e}_{j}\right)\left({G}_{\delta }^{\omega }\left(y,z\right)-\delta {G}_{\delta }^{\omega }\left(y+{e}_{j},z\right)\right)|\le \left(2d-1\right)\frac{{G}_{\delta }^{\omega }\left(y,z\right)}{{G}_{\delta }^{\omega }\left(y,y\right)},$  which immediately implies the estimate (4).
The inequalities (3) and (4) imply that $|{m}_{{e}_{j}}{E}_{\mu }\left[{G}_{\delta }\left(0,z\right){G}_{\delta }\left(z,y\right)\right]-{E}_{\mu }\left[{G}_{\delta }\left(0,z\right){G}_{\delta }\left(z,y\right)\omega \left(y,y+{e}_{j}\right)\right]|\le \frac{2d}{\gamma }{E}_{\mu }\left[{G}_{\delta }\left(0,z\right){G}_{\delta }\left(z,y\right)\right].$  This gives the following estimate for ${\stackrel{~}{\omega }}^{z}$  $|{m}_{{e}_{j}}-{\stackrel{~}{\omega }}^{z}\left(y,y+{e}_{j}\right)|\le \frac{2d}{\gamma }.$  The same procedure gives the same estimate for ${\stackrel{~}{\omega }}^{z+{e}_{i}}$  .
Hence, we see that ${E}_{\stackrel{~}{\mu }}\left[{G}_{\delta }^{\omega }\left(z,z\right)\omega \left(z,z+{e}_{i}\right)\right]={G}_{\delta }^{m+\Delta m}\left(z,z\right)\left({m}_{i}+\Delta m\left(z,z+{e}_{i}\right)\right),$  where $\Delta m\left(z,z+{e}_{i}\right)$  is a correction to the homogeneous transition probability $\left({m}_{i}\right)$  uniformly bounded by $|\Delta m|\le \frac{2d}{\gamma }.$  The same reasoning holds for ${E}_{\stackrel{~}{\mu }}\left[{G}_{\delta }^{\omega }\left(z+{e}_{i},z\right)\omega \left(z,z+{e}_{i}\right)\right]={G}_{\delta }^{m+\Delta m}\left(z+{e}_{i},z\right)\left({m}_{i}+\Delta m\left(z,z+{e}_{i}\right)\right),$  even if the correction term $\Delta m$  is not the same.
Step 2: We compare now the Green function ${G}_{\delta }^{m+\Delta m}$  with ${G}_{\delta }^{m}$  .
This is done in [4, but we reproduce the main lines of the proof, since we want to obtain explicit bounds. We first introduce the symmetrizing function ${\phi }^{m}\left(z\right){=}^{d}{\prod }_{i=1}{\sqrt{\frac{{m}_{{e}_{i}}}{{m}_{-{e}_{i}}}}}^{{z}_{i}}.$  The Green function ${G}_{\delta }^{m}$  is transformed into ${G}_{\delta }^{m}={M}_{\phi }^{-1}{G}_{\delta {k}_{m}}^{s}{M}_{\phi },\left(5\right)$  where ${M}_{\phi }$  is the operator of multiplication by $\phi$  , and ${G}_{\delta {k}_{m}}^{s}$  is the Green function of the symmetric random walk with transition probability ${s}_{{e}_{i}}={s}_{-{e}_{i}}=\frac{\sqrt{{m}_{{e}_{i}}{m}_{-{e}_{i}}}}{2{\sum }_{k=1}^{2d}\sqrt{{m}_{{e}_{k}}{m}_{-{e}_{k}}}},i=1,\dots ,d,$  with killing rates $\delta {k}_{m}$  where ${k}_{m}=2{\sum }_{k=1}^{2d}\sqrt{{m}_{{e}_{k}}{m}_{-{e}_{k}}}.$  We refer to Step 2 of the proof of Lemma 4.3 in [4for precisions about this fact.
Hence, we see that ${\sum }_{y}{G}_{\delta k}^{s}\left(0,y\right)\le \frac{1}{1-\delta {k}_{m}},$  which means that $\parallel {G}_{\delta k}^{s}{\parallel }_{\infty }\le \frac{1}{1-\delta {k}_{m}}.$  We also have ${G}_{\delta }^{m+\Delta m}-{G}_{\delta }^{m}=-{G}_{\delta }^{m}\left(I-\left(I-\delta \Delta {P}_{m}{G}_{\delta }^{m}{\right)}^{-1}\right),$  where $\Delta {P}_{m}$  is the matrix $\left(\Delta {P}_{m}{\right)}_{x,x+{e}_{i}}=\Delta m\left(x,x+{e}_{i}\right)$  (and null anywhere else).
Thus, we get ${G}_{\delta }^{m+\Delta m}-{G}_{\delta }^{m}=\delta {M}_{\phi }^{-1}{G}_{\delta {k}_{m}}^{s}{M}_{\phi }\Delta {P}_{m}{M}_{\phi }^{-1}{G}_{\delta {k}_{m}}^{s}{\left(I-\delta {M}_{\phi }\Delta {P}_{m}{M}_{\phi }^{-1}{G}_{\delta {k}_{m}}^{s}\right)}^{-1}{M}_{\phi },$  but $\parallel {M}_{\phi }\Delta {P}_{m}{M}_{\phi }^{-1}{\parallel }_{\infty }\le \left({max}_{i}\phi \left({e}_{i}\right)\right)\frac{2d}{\gamma }$  so that we get $\left({G}_{\delta }^{m+\Delta m}-{G}_{\delta }^{m}\right)\left(x,y\right)\le \phi \left(y-x\right)\frac{2d}{\gamma }\frac{1}{1-{k}_{m}}{\eta }_{m}\frac{1}{1-\frac{2d}{\gamma }{\eta }_{m}},$  and ${G}_{\delta }^{m+\Delta m}\left(x,y\right)\le \phi \left(y-x\right)\frac{1}{1-{k}_{m}}\frac{1}{1-\frac{2d}{\gamma }{\eta }_{m}}.$  This implies that, for all $i\in \left\{1,...,2d\right\}$  , $|{\stackrel{^}{\omega }}_{\delta }\left(z,z+{e}_{i}\right)-{m}_{i}+\frac{{m}_{i}}{\gamma -1}\left({G}_{\delta }^{m}\left(0,0\right)-{G}_{\delta }^{m}\left({e}_{i},0\right)-1\right)|$
 $\begin{array}{c}\le \frac{1}{\gamma -1}\left(2\frac{2d}{\gamma }{\eta }_{m}^{2}\frac{1}{1-\frac{2d}{\gamma }{\eta }_{m}}+2\frac{2d}{\gamma }{\eta }_{m}\frac{1}{1-\frac{2d}{\gamma }{\eta }_{m}}\right)\end{array}$
 $\begin{array}{c}\le 8\frac{d}{{\gamma }^{2}}\frac{{\eta }_{m}^{2}}{1-\frac{2d}{\gamma }{\eta }_{m}}\left(6\right)\end{array}$
(we used here ${\eta }_{m}\ge 1$  ).
The sum ${\sum }_{i=1}^{2d}\left({m}_{i}-\frac{{m}_{i}}{\gamma -1}\left({G}_{\delta }^{m}\left(0,0\right)-{G}_{\delta }^{m}\left({e}_{i},0\right)-1\right)\right).{e}_{i}$  tends to ${d}_{m}\left(1-\frac{1}{\gamma -1}\left({G}^{m}\left(0,0\right)-1\right)\right)$  when $\delta$  tends to 1. Indeed, the sum ${\sum }_{i=1}^{2d}{m}_{i}{G}^{m}\left({e}_{i},0\right).{e}_{i}$  cancels, due to the fact that for each $i\in \left\{1,...,d\right\}$  ${m}_{{e}_{i}}{G}^{m}\left({e}_{i},0\right)$  and ${m}_{-{e}_{i}}{G}^{m}\left(-{e}_{i},0\right)$  are both equal to the common value $\sqrt{{m}_{{e}_{i}}{m}_{-{e}_{i}}}{G}_{k}^{s}\left({e}_{i},0\right)$  (cf formula (5)).
The triangular inequality combined with the $2d$  inequalities (6) gives that, for all $z$  , ${limsup}_{\delta \to 1}\parallel {d}_{{\stackrel{^}{\omega }}_{\delta }}\left(z\right)-{d}_{m}\left(1-\frac{1}{\gamma -1}\left({G}^{m}\left(0,0\right)-1\right)\right)\parallel \le 16{\left(\frac{d}{\gamma }\right)}^{2}\frac{{\eta }_{m}^{2}}{1-\frac{2d}{\gamma }{\eta }_{m}},$  where ${d}_{{\stackrel{^}{\omega }}_{\delta }}\left(z\right)={\sum }_{k=1}^{2d}\stackrel{^}{\omega }\left(z,z+{e}_{k}\right){e}_{k}$  is the local drift of the transition probability ${\stackrel{^}{\omega }}_{\delta }$  . Proposition 3.2 of [4allows to conclude. References

1. Comets, F., Zeitouni, O., A law of large numbers for random walks in random mixing environments. Ann. Probab. 32 (2004), no. 1B, 880–914.
2. Enriquez, N., Sabot, C., Edge oriented reinforced random walks and RWRE. C. R. Math. Acad. Sci. Paris 335 (2002), no. 11, 941–946.
3. Kalikow, S., Generalized random walk in a random environment. Ann. Probab. 9 (1981), no. 5, 753–768.
4. Sabot, C., Ballistic random walks in random environment at low disorder. Ann. Probab. 32 (2004), no. 4, 2996–3023.
5. Solomon, F., Random walks in a random environment. Ann. Probability 3 (1975), 1–31.
6. Sznitman, A-S., Slowdown estimates and central limit theorem for random walks in random environment. J. Eur. Math. Soc. (JEMS) 2 (2000), no. 2, 93–143.
7. Sznitman, A-S., Zerner, M., A law of large numbers for random walks in random environment. Ann. Probab. 27 (1999), no. 4, 1851–1869.
8. Zeitouni, O., Random walks in random environment. Lectures on probability theory and statistics, 189–312, Lecture Notes in Math., 1837, Springer, Berlin, 2004.

Laboratoire de Probabilites de Paris 6, 4 place Jussieu, 75252 Paris Cedex 05 E-mail address : enriquez@ccr.jussieu.fr Unite de Mathematiques Pures et Appliquees, 46, Allee d'Italie F-69364 Lyon Cedex 07 E-mail address : christophe.sabot@umpa.ens-lyon.fr