arxiv: v3 [math.pr] 7 May 2012

Similar documents
A Scaling Analysis of a Cat and Mouse Markov Chain

Ergodic Theorems. Samy Tindel. Purdue University. Probability Theory 2 - MA 539. Taken from Probability: Theory and examples by R.

Functional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals

arxiv: v2 [math.pr] 4 Sep 2017

P i [B k ] = lim. n=1 p(n) ii <. n=1. V i :=

Lecture 7. µ(x)f(x). When µ is a probability measure, we say µ is a stationary distribution.

Lecture 7. We can regard (p(i, j)) as defining a (maybe infinite) matrix P. Then a basic fact is

Markov processes and queueing networks

Lecture 20: Reversible Processes and Queues

Chapter 7. Markov chain background. 7.1 Finite state space

Other properties of M M 1

Stochastic process. X, a series of random variables indexed by t

Necessary and sufficient conditions for strong R-positivity

Brownian Motion. 1 Definition Brownian Motion Wiener measure... 3

Queueing Networks and Insensitivity

Statistics 150: Spring 2007

Selected Exercises on Expectations and Some Probability Inequalities

Strong approximation for additive functionals of geometrically ergodic Markov chains

9 Brownian Motion: Construction

Stability and Rare Events in Stochastic Models Sergey Foss Heriot-Watt University, Edinburgh and Institute of Mathematics, Novosibirsk

Applied Stochastic Processes

arxiv: v1 [math.pr] 6 Jan 2014

Lecture 5. If we interpret the index n 0 as time, then a Markov chain simply requires that the future depends only on the present and not on the past.

Mean-field dual of cooperative reproduction

Positive and null recurrent-branching Process

Lecture 12. F o s, (1.1) F t := s>t

MS&E 321 Spring Stochastic Systems June 1, 2013 Prof. Peter W. Glynn Page 1 of 10

Lecture 11: Introduction to Markov Chains. Copyright G. Caire (Sample Lectures) 321

STOCHASTIC PROCESSES Basic notions

Markov Chains and Stochastic Sampling

MATH 56A: STOCHASTIC PROCESSES CHAPTER 2

Modern Discrete Probability Spectral Techniques

Let (Ω, F) be a measureable space. A filtration in discrete time is a sequence of. F s F t

Stochastic Process (ENPC) Monday, 22nd of January 2018 (2h30)

Lectures on Stochastic Stability. Sergey FOSS. Heriot-Watt University. Lecture 4. Coupling and Harris Processes

Stochastic Processes (Week 6)

Uniformly Uniformly-ergodic Markov chains and BSDEs

The strictly 1/2-stable example

Properties of an infinite dimensional EDS system : the Muller s ratchet

Brownian motion. Samy Tindel. Purdue University. Probability Theory 2 - MA 539

Stochastic flows associated to coalescent processes

Lecture Notes 7 Random Processes. Markov Processes Markov Chains. Random Processes

SMSTC (2007/08) Probability.

Markov Chains CK eqns Classes Hitting times Rec./trans. Strong Markov Stat. distr. Reversibility * Markov Chains

T. Liggett Mathematics 171 Final Exam June 8, 2011

Universal examples. Chapter The Bernoulli process

Logarithmic scaling of planar random walk s local times

Interlude: Practice Final

Regular Variation and Extreme Events for Stochastic Processes

4 Sums of Independent Random Variables

Theory and Applications of Stochastic Systems Lecture Exponential Martingale for Random Walk

LIMITS FOR QUEUES AS THE WAITING ROOM GROWS. Bell Communications Research AT&T Bell Laboratories Red Bank, NJ Murray Hill, NJ 07974

Brownian Motion. Chapter Definition of Brownian motion

Lecture 17 Brownian motion as a Markov process

PERCOLATION IN SELFISH SOCIAL NETWORKS

Latent voter model on random regular graphs

Irreducibility. Irreducible. every state can be reached from every other state For any i,j, exist an m 0, such that. Absorbing state: p jj =1

Lecture 9. d N(0, 1). Now we fix n and think of a SRW on [0,1]. We take the k th step at time k n. and our increments are ± 1

Itô s excursion theory and random trees

IEOR 8100: Topics in OR: Asymptotic Methods in Queueing Theory. Fall 2009, Professor Whitt. Class Lecture Notes: Wednesday, September 9.

Markov Chains. X(t) is a Markov Process if, for arbitrary times t 1 < t 2 <... < t k < t k+1. If X(t) is discrete-valued. If X(t) is continuous-valued

SUMMARY OF RESULTS ON PATH SPACES AND CONVERGENCE IN DISTRIBUTION FOR STOCHASTIC PROCESSES

1. Stochastic Processes and filtrations

Practical conditions on Markov chains for weak convergence of tail empirical processes

A Type of Shannon-McMillan Approximation Theorems for Second-Order Nonhomogeneous Markov Chains Indexed by a Double Rooted Tree

Simultaneous drift conditions for Adaptive Markov Chain Monte Carlo algorithms

Classification of Countable State Markov Chains

Lecture 21 Representations of Martingales

1 Continuous-time chains, finite state space

Recap. Probability, stochastic processes, Markov chains. ELEC-C7210 Modeling and analysis of communication networks

MATH 564/STAT 555 Applied Stochastic Processes Homework 2, September 18, 2015 Due September 30, 2015

Lecture 6. 2 Recurrence/transience, harmonic functions and martingales

ON METEORS, EARTHWORMS AND WIMPS

MAT 570 REAL ANALYSIS LECTURE NOTES. Contents. 1. Sets Functions Countability Axiom of choice Equivalence relations 9

Weak quenched limiting distributions of a one-dimensional random walk in a random environment

CDA6530: Performance Models of Computers and Networks. Chapter 3: Review of Practical Stochastic Processes

MS&E 321 Spring Stochastic Systems June 1, 2013 Prof. Peter W. Glynn Page 1 of 10. x n+1 = f(x n ),

Figure 10.1: Recording when the event E occurs

Math 6810 (Probability) Fall Lecture notes

BRANCHING PROCESSES 1. GALTON-WATSON PROCESSES

2. Transience and Recurrence

12 Markov chains The Markov property

Markov Chains Handout for Stat 110

Example: physical systems. If the state space. Example: speech recognition. Context can be. Example: epidemics. Suppose each infected

INTRODUCTION TO MARKOV CHAINS AND MARKOV CHAIN MIXING

Continuous time Markov chains

I forgot to mention last time: in the Ito formula for two standard processes, putting

Central limit theorems for ergodic continuous-time Markov chains with applications to single birth processes

Markov Chains. Chapter Existence and notation. B 2 B(S) and every n 0,

Examples of Countable State Markov Chains Thursday, October 16, :12 PM

2. The averaging process

18.175: Lecture 8 Weak laws and moment-generating/characteristic functions

3. The Voter Model. David Aldous. June 20, 2012

Stochastic integration. P.J.C. Spreij

TAUBERIAN THEOREM FOR HARMONIC MEAN OF STIELTJES TRANSFORMS AND ITS APPLICATIONS TO LINEAR DIFFUSIONS

Asymptotics for posterior hazards

Weighted Sums of Orthogonal Polynomials Related to Birth-Death Processes with Killing

FUNCTIONAL CENTRAL LIMIT (RWRE)

Invariance Principle for Variable Speed Random Walks on Trees

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science

Transcription:

The Annals of Applied Probability 2012, Vol. 22, No. 2, 792 826 DOI: 10.1214/11-AAP785 c Institute of Mathematical Statistics, 2012 arxiv:0905.2259v3 [math.pr] 7 May 2012 A SCALING ANALYSIS OF A CAT AND MOUSE MARKOV CHAIN By Nelly Litvak 1 and Philippe Robert University of Twente and INRIA Paris Rocquencourt If C n is a Markov chain on a discrete state space S, a Markov chain C n,m n on the product space S S, the cat and mouse Markov chain, is constructed. The first coordinate of this Markov chain behaves like the original Markov chain and the second component changes only when both coordinates are equal. The asymptotic properties of this Markov chain are investigated. A representation of its invariant measure is, in particular, obtained. When the state space is infinite it is shown that this Markov chain is in fact null recurrent if the initial Markov chain C n is positive recurrent and reversible. In this context, the scaling properties of the location of the second component, the mouse, are investigated in various situations: simple random walks in Z and Z 2 reflected a simple random walk in N and also in a continuous time setting. For several of these processes, a time scaling with rapid growth gives an interesting asymptotic behavior related to limiting results for occupation times and rare events of Markov processes. 1. Introduction. The PageRank algorithm of Google, as designed by Brin and Page [10] in 1998, describes the web as an oriented graph S whose nodes are the web pages and the html links between these web pages, the links of the graph. In this representation, the importance of a page is defined as its weight for the stationary distribution of the associated random walk on the graph. Several off-line algorithms can be used to estimate this equilibrium distribution on such a huge state space, they basically use numerical procedures matrix-vector multiplications. See Berkhin [4], for example. Several on-line algorithms that update the ranking scores while exploring Received October 2010; revised April 2011. 1 Supported by The Netherlands Organisation for Scientific Research NWO under Meervoud Grant 632.002.401. AMS 2000 subject classifications. 60J10, 90B18. Key words and phrases. Cat and mouse Markov chains, scaling of null recurrent Markov chains. This is an electronic reprint of the original article published by the Institute of Mathematical Statistics in The Annals of Applied Probability, 2012, Vol. 22, No. 2, 792 826. This reprint differs from the original in pagination and typographic detail. 1

2 N. LITVAK AND P. ROBERT the graph have been recently proposed to avoid some of the shortcomings of off-line algorithms, in particular, in terms of computational complexity. The starting point of this paper is an algorithm designed by Abiteboul et al. [1] to compute the stationary distribution of a finite recurrent Markov chain. In this setting, to each node of the graph is associated a number, the cash of the node. The algorithm works as follows: at a given time, the node x with the largest value V x of cash is visited, V x is set to 0 and the value of the cash of each of its d x neighbors is incremented by V x /d x. Another possible strategy to update cash variables is as follows: a random walker updates the values of the cash at the nodes of its random path in the graph. This policy is referred to as the Markovian variant. Both strategies have the advantage of simplifying the data structures necessary to manage the algorithm. It turns out that the asymptotic distribution, in terms of the number of steps of the algorithm, of the vector of the cash variables gives an accurate estimation of the equilibrium distribution; see Abiteboul et al. [1] for the complete description of the procedure to get the invariant distribution. See also Litvak and Robert [23]. The present paper does not address the problem of estimating the accuracy of these algorithms, it analyzes the asymptotic properties of a simple Markov chain which appears naturally in this context. Cat and mouse Markov chain. IthasbeenshowninLitvakandRobert[23] that, for the Markovian variant of the algorithm, the distribution of the vector of the cash variables can be represented with the conditional distributions of a Markov chain C n,m n on the discrete state space S S. The sequence C n, representing the location of the cat, is a Markov chain with transition matrix P =px,y associated to the random walk on the graph S. The second coordinate, the location of the mouse, M n has the following dynamic: If M n C n, then M n+1 =M n, If M n =C n, then, conditionally on M n, the random variable M n+1 has distribution pm n,y,y S and is independent of C n+1. This can be summarized as follows: the cat moves according to the transition matrix P =px,y and the mouse stays idle unless the cat is at the same site, in which case the mouse also moves independently according to P = px,y. The terminology cat and mouse problem is also used in a somewhat different way in game theory, the cat playing the role of the adversary. See Coppersmith et al. [11] and references therein. The asymptotic properties of this interesting Markov chain C n,m n for a number of transition matrices P are the subject of this paper. In particular, the asymptotic behavior of the location mouse M n is investigated. The

A SCALING ANALYSIS OF A CAT AND MOUSE MARKOV CHAIN 3 distribution of M n plays an important role in the algorithm designed by Abiteboul et al. [1]; see Litvak and Robert [23] for further details. It should be noted that M n is not, in general, a Markov chain. Outline of the paper. Section 2 analyzes the recurrence properties of the Markov chain C n,m n when the Markov chain C n is recurrent. A representation of the invariant measure of C n,m n in terms of the reversed process of C n is given. Since the mouse moves only when the cat arrives at its location, it may seem quite likely that the mouse will spend most of the time at nodes which are unlikely for the cat. It is shown that this is indeed the case when the state space is finite and if the Markov chain C n is reversible but not in general. When the state space is infiniteandif themarkov chain C n is reversible, it turns out that the Markov chain C n,m n is in fact null recurrent. A precise description of the asymptotic behavior of the sequence M n is done via a scaling in time and space for several classes of simple models. Interestingly, the scalings used are quite diverse, as it will be seen. They are either related to asymptotics of rare events of ergodic Markov chains or to limiting results for occupation times of recurrent random walks: 1 Symmetric simple random walks. The cases of symmetric simple random walks in Z d with d=1 and 2 are analyzed in Section 3. Note that for d 3 the Markov chain C n is transient so that in this case the location of the mouse does not change with probability 1 after some random time: In the one-dimensional case, d=1, if M 0 =C 0 =0, on the linear time scale t nt, as n gets large, it is shown that the location of the mouse is of the order of 4 n. More precisely, the limit in distribution of the process M nt / 4 n,t 0 is a Brownian motion B 1 t taken at the local time at 0 of another independent Brownian motion B 2 t. See Theorem 2 below. This result can be roughly described as follows. Under this linear time scale the location of the cat, a simple symmetrical random walk, is of the order of n by Donsker s theorem. It turns out that it will encounter n times the mouse. Since the mouse moves only when it encounters the cat and that it also follows the sample path of a simple random walk, after n steps its order of magnitude will be therefore of the order of 4 n. When d=2, on the linear time scale t nt, the location of the mouse is of the order of logn. More precisely, the finite marginals of the rescaled processes M expnt / n,t 0 converge to the corresponding finite marginals of a Brownian motion in R 2 on a time scale which

4 N. LITVAK AND P. ROBERT is an independent discontinuous stochastic process with independent and nonhomogeneous increments. 2 Reflected simple random walk. Section 4 investigates the reflected simple random walk on the integers. A jump of size +1 resp., 1 occurs with probability p [resp., 1 p] and the quantity ρ=p/1 p is assumed to be strictly less than 1 so that the Markov chain C n is ergodic. Ifthelocation ofthemouseisfaraway fromtheorigin,that is,m 0 =n with n large and the cat is at equilibrium, a standard result shows that it takes a duration of time of the order of ρ n for the cat to hit the mouse. This suggests an exponential time scale t ρ n t to study the evolution of the successive locations of the mouse. For this time scale it is shownthat thelocation of themouseis still of theorderof n as long as t<w where W is some nonintegrable random variable. At time t=w on the exponential time scale, the mouse has hit 0 and after that time the process M tρ n /n oscillates between 0 and above 1/2 on every nonempty time interval. 3 Continuous time random walks. Section 5 introduces the cat and mouse process for continuous time Markov processes. In particular, a discrete Ornstein Uhlenbeck process, the M/M/ queue, is analyzed. This is a birth and death process whose birth rates are constant and the death rate at n N is proportional to n. When M 0 =n, contrary to the case of the reflected random walk, there does not seem to exist a time scale for which a nontrivial functional theorem holds for the corresponding rescaled process. Instead, it is possible to describe the asymptotic behavior of the location of the mouse after the pth visit of the cat. It has a multiplicative representation of the form nf 1 F 2 F p where F p are i.i.d. random variables on [0, 1]. The examples analyzed are quite specific. They are, however, sufficiently representative of the different situations for the dynamic of the mouse: 1 One considers the case when an integer valued Markov chain C n is ergodic and the initial location of the mouse is far away from 0. The correct time scale to investigate the evolution of the location of the mouse is given by the duration of time for the occurrence of a rare event for the original Markov chain. When the cat hits the mouse at this level, before returning to the neighborhood of 0, it changes the location of the mouse by an additive resp., multiplicative step in the case of the reflected random walk resp., M/M/ queue. 2 For null recurrent homogeneous random walks, the distribution of the duration of times between two visits of the cat to the mouse do not depend on the location of the mouse but it is nonintegrable. The main problem is therefore to get a functional renewal theorem associated to an i.i.d. sequence T n of nonnegative random variables such that ET 1 =+.

A SCALING ANALYSIS OF A CAT AND MOUSE MARKOV CHAIN 5 More precisely, if Nt= i 1 1 {T1 + +T i t}, one has to find φn such that the sequence of processes Nnt/φn,t 0 converges as n goes to infinity. When the tail distribution of T 1 has a polynomial decay, several technical results are available. See Garsia and Lamperti [12], for example. This assumption is nevertheless not valid for the two-dimensional case. In any case, it turns out that the best way especially for d=2 to get such results is to formulate the problem in terms of occupation times of Markov processes for which several limit theorems are available. This is the key of the results in Section 3. The fact that for all the examples considered jumps occur on the nearest neighbors does not change this qualitative behavior. Under more general conditions analogous results should hold. Additionally, this simple setting has the advantage of providing explicit expressions for most of the constants involved. 2. The cat and mouse Markov chain. In this section we consider a general transition matrix P = px,y,x,y S on a discrete state space S. Throughout the paper, it is assumed that P is aperiodic, irreducible without loops, that is, px,x=0 for all x S and with an invariant measure π. Note that it is not assumed that π has a finite mass. The sequence C n will denote a Markov chain with transition matrix P =px,y. It will represent the sequence of nodes which are sequentially updated by the random walker. The transition matrix of the reversed Markov chain C n is denoted by and, for y S, one defines p x,y= πy πx py,x H y =inf{n>0:c n =y} and H y =inf{n>0:c n =y}. The Markov chain C n,m n on S S referred to as the cat and mouse Markov chain is introduced. Its transition matrix Q=q, is defined as follows: for x, y, z S, { q[x,y,z,y]=px,z, if x y; 1 q[y,y,z,w]=py,zpy,w. The process C n [resp., M n ] will be defined as the position of the cat resp., the mouse. Note that the position C n of the cat is indeed a Markov chain with transition matrix P =p,. The position of the mouse M n changes only when the cat is at the same position. In this case, starting

6 N. LITVAK AND P. ROBERT from x S they both move independently according to the stochastic vector px,. Since thetransition matrix of C n is assumed to beirreducibleand aperiodic, it is not difficult to check that the Markov chain C n,m n is aperiodic and visits with probability 1 all the elements of the diagonal of S S. In particular, there is only one irreducible component. Note that C n,m n itself is not necessarily irreducible on S S, as the following example shows: take S ={0,1,2,3} and the transition matrix p0,1 =p2,3=p3,1=1 and p1,2=1/2=p1,0; in this case the element 0,3 cannot be reached starting from 1, 1. Theorem 1 Recurrence. The Markov chain C n,m n on S S with transition matrix Q defined by relation 1 is recurrent: the measure ν defined as 2 H y νx,y=πxe x n=1 pc n,y, x,y S, is invariant. Its marginal on the second coordinate is given by, for y S, ν 2 y def. = x Sνx,y=E π pc 0,yH y, and it is equal to π on the diagonal, νx,x=πx for x S. In particular, with probability 1, the elements of S S for which ν is nonzero are visited infinitely often and ν is, up to a multiplicative coefficient, the unique invariant measure. The recurrence property is not surprising: the positive recurrence property of the Markov chain C n shows that cat and mouse meet infinitely often with probability one. The common location at these instants is a Markov chain with transition matrix P and therefore recurrent. Note that the total mass of ν, νs S= y SE π pc 0,yH y can be infinite when S is countable. See Kemeny et al. [20] for an introduction on recurrence properties of discrete countable Markov chains. The measure ν 2 on S is related to the location of the mouse under the invariant measure ν. Proof of Theorem 1. FromtheergodicityofC n itisclearthatνx,y is finite for x, y S. One has first to check that ν satisfies the equations of invariant measure for the Markov chain C n,m n, 3 νx,y= z y νz,ypz,x+ z νz,zpz,xpz,y, x,y S.

For x, y S, 4 and 5 A SCALING ANALYSIS OF A CAT AND MOUSE MARKOV CHAIN 7 z yνz,ypz,x= z y H y πxp x,ze z H y =πxe x νz, zpz, xpz, y z S n=2 pc n,y n=1 pc n,y = H z 1 πxp x,zpz,ye z pcn,z. z S n=0 The classical renewal argument for the invariant distribution π of the Markov chain Cn, and any bounded function f on S, gives that H 1 z 1 E π f= E z Hz E z fcn ; see Theorem 3.2, page 12, of Asmussen [3], for example. In particular, we have πz=1/e z Hz, and H z 1 E z pcn,z =E z HzE π pc0,z= x S πxpx,z πz n=0 6 = πz πz =1. n=0 Substituting the last identity into 5, we obtain z Sνz,zpz,xpz,y= πxp x,zpz,y z S 7 =πxe x pc1,y. Relations 3 5 and 7 show that ν is indeed an invariant distribution. At the same time, from 6 one gets the identity νx,x=πx for x S. x Sνx,y= The second marginal is given by, for y S, πxe x pct,y1 {H y t} t 1 x S = t 1E π pc t,y1 {H y t}

8 N. LITVAK AND P. ROBERT = πxp x,z 1 p z 1,z 2 p z t 1,z t pz t,y x Sz 1,...,z t 1 yz t S = pz 1,xpz 2,z 1 pz t,z t 1 πz t pz t,y x Sz 1,...,z t 1 yz t S = t 1E π pc 0,y1 {Hy t}=e π pc 0,yH y, and the theorem is proved. The representation 2 of the invariant measure can be obtained formally through an iteration of the equilibrium equations 3. Since the first coordinate of C n,m n is a Markov chain with transition matrix P and ν is the invariant measure for C n,m n, the first marginal of ν is thus equal to απ for some α>0, that is, νx,y=απx, x S. y The constant α is in fact the total mass of ν. In particular, from 2, one gets that the quantity hx def. = y S H y E x n=1 pc n,y, x S, is independent of x S and equal to α. Note that the parameter α can be infinite. Proposition 1 Location of the mouse in the reversible case. If C n is a reversible Markov chain, with the definitions of the above theorem, for y S, the relation ν 2 y=1 πy holds. If the state space S is countable, the Markov chain C n,m n is then null recurrent. Proof. For y S, by reversibility, ν 2 y=e π pc 0,yH y = x πxpx,ye x H y = x πypy,xe x H y =πye y H y 1 =1 πy. The proposition is proved. Corollary 1 Finite state space. If the state space S is finite with cardinality N, then C n,m n converges in distribution to C,M such

A SCALING ANALYSIS OF A CAT AND MOUSE MARKOV CHAIN 9 that 8 with H y PC =x,m =y=α 1 πxe x n=1 α= y SE π pc 0,yH y pc n,y, x,y S, in particular, PC = M = x = α 1 πx. If the Markov chain C n is reversible, then PM =y= 1 πy N 1. Tetali [29] showed, via linear algebra, that if C n is a general recurrent Markov chain, then 9 E π pc 0,yH y N 1. y S See also Aldous and Fill [2]. It follows that the value α =N 1 obtained for reversible chains is the maximal possible value of α. The constant α 1 is the probability that the cat and mouse are at the same location. In the reversible case, Corollary 1 implies the intuitive fact that the less likely a siteis forthecat, themorelikely it is forthemouse. Thisis, however, false in general. Consider a Markov chain whose state space S consists of r cycles with respective sizes m 1,...,m r with one common node 0, r S ={0} {k,i:1 i m k }, k=1 and with the following transitions: for 1 k r and 2 i m k, pk,i,k,i 1=1, pk,1,0=1 and p0,k,m k = 1 r. Define m=m 1 +m 2 + +m r. It is easy to see that π0= r and πy= 1 m+r m+r, y S {0}. One gets that for the location of the mouse, for y S, { πym mk +1, if y=k,m ν 2 y=e π pc 0,yH y = k,1 k r, πy, otherwise. Observe that for any y distinct from 0 and k,m k, we have π0>πy and ν 2 0>ν 2 y; the probability to find a mouse in 0 is larger than in y. Note that in this example one easily obtains c=1/r. 3. Random walks in Z and Z 2. In this section the asymptotic behavior of the mouse when the cat follows a recurrent random walk in Z and Z 2 is

10 N. LITVAK AND P. ROBERT analyzed. The jumps of the cat are uniformly distributed on the neighbors of the current location. 3.1. One-dimensional random walk. The transition matrix P of this random walk is given by px,x+1= 1 2 =px,x 1, x Z. Decomposition into cycles. If the cat and the mouse start at the same location, they stay together a random duration of time G which is geometrically distributed with parameter 1/2. Once they are at different locations for the first time, they are at distance 2 so that the duration of time T 2 until they meet again has the same distribution as the hitting time of 0 by the random walk which starts at 2. The process C n,m n,0 n G+T 2 is definedas acycle. Thesamplepath of themarkov chain C n,m n canthus be decomposed into a sequence of cycles. It should be noted that, during a cycle, the mouse moves only during the period with duration G. Since one investigates the asymptotic properties of the sample paths of M n on the linear time scale t nt for n large, to get limit theorems one should thus estimate the number of cycles that occur in a time interval [0,nt]. For this purpose, we compare the cycles of the cat and mouse process to the cycles of a simple symmetric random walk, which are the time intervals between two successive visits to zero by the process C n. Observe that a cycle of C n is equal to 1+T 1, where T 1 is the time needed to reach zero starting from 1. Further, T 2 is the sum of two independent random variables distributed as T 1. Hence, one guesses that on the linear time scale t nt the number of cycles on [0,nt] for C n,m n is asymptotically equivalent to 1/2 of the number of cycles on [0,nt] for C n, as n. It is well known that the latter number is of the order n. Then the mouse makes order of n steps of a simple symmetric random walk, and thus its location must be of the order 4 n. To make this argument precise, we first prove technical Lemma 1, which says that only o n of C n -cycles can befitted into the time interval of the order n. Next, Lemma 2 proves that the number of cycles of length T 2 +2 on [0,nt], scaled by n, converges to 1/2 of the local time of a Brownian motion, analogously to the corresponding result for the number of cycles of a simple symmetric random walk [22]. Finally, the main limiting result for the location of the mouse is given by Theorem 2. Lemma 1. For any x, ε>0 and K >0, k+ ε 1 n n lim P n + inf 0 k x n i=k 1+T 1,i K =0,

A SCALING ANALYSIS OF A CAT AND MOUSE MARKOV CHAIN 11 where T 1,i are i.i.d. random variables with the same distribution as the first hitting time of 0 of C n, T 1 =inf{n>0:c n =0 with C 0 =1}. Proof. If E is an exponential random variable with parameter 1 independent of the sequence T 1,i, by using the fact that, for u 0,1, Eu T 1 =1 1 u 2 /u, then for n 2, ε n 1 logp 1+T 1,i E n Denote by i=0 m n = inf 0 k x n the above relation gives x n k+ ε 1 Pm n E P n hence, k=0 + n=2 n i=k = ε n log1 1 e 2/ n ε 4 n. k+ ε 1 n n i=k 1+T 1,i E Pm n E<+ 1+T 1,i x n +1e ε 4 n, and, consequently, with probability 1, there exists N 0 such that, for any n N 0, we have m n >E. Since PE K>0, the lemma is proved. Lemma 2. Let, for n 1, T 2,i i.i.d. random variables with the same distribution as T 2 =inf{k>0:c k =0 with C 0 =2} and u n = + l=1 1 { l k=1 2+T 2,k<n}, then the process u tn / n converges in distribution to L B t/2, where L B t is the local time process at time t 0 of a standard Brownian motion. Proof. The variable T 2 can be written as a sum T 1 +T 1 of independent randomvariables T 1 and T 1 havingthesamedistributionas T 1 definedinthe above lemma. For k 1, the variable T 2,k can be written as T 1,2k 1 +T 1,2k. Clearly, 1 + 1 2 { l k=1 1+T 1,k<n} 1 2 u n 1 + 1 2 { l k=1 1+T 1,k<n}. l=1 l=1

12 N. LITVAK AND P. ROBERT Furthermore, + n 1 1 { l k=1 1+T 1,k<n},n 1 dist. = r n def. = 1 {Cl =0},n 1, l=1 where C n is the symmetric simple random walk. A classical result by Knight [22] see also Borodin [8] and Perkins [25] gives that the process r nt / n converges in distribution to L B t as n gets large. The lemma is proved. l=1 The main result of this section can now be stated. Theorem 2 Scaling of the location of the mouse. If C 0,M 0 N 2, the convergence in distribution 1 dist. 4 M n nt,t 0 = B 1 L B2 t,t 0 lim n + holds, where B 1 t and B 2 t are independent standard Brownian motions on R and L B2 t is the local time process of B 2 t at 0. The location of the mouse at time T is therefore of the order of 4 T as T gets large. The limiting process can be expressed as a Brownian motion sloweddownbytheprocessofthelocaltimeat 0ofanindependentBrownian motion. The quantity L B2 T can be interpreted as the scaled duration of time the cat and the mouse spend together. Proof of Theorem 2. Without loss of generality, one can assume that C 0 =M 0. A coupling argument is used. Take: i.i.d. geometric random variables G i such that PG 1 p=1/2 p 1 for p 1; Ck aandcb j,k,j 1,i.i.d.independentsymmetricrandomwalksstarting from 0; and assume that all these random variables are independent. One denotes, for m=1, 2 and j 1, Define T b m,j =inf{k 0:C b j,k =m}. { C a C k,m k = k,ck a, 0 k<g 1, CG a 1 2I 1 +I 1 C1,k G b 1,CG a 1, G 1 k τ 1, with I 1 =CG a 1 CG a 1 1, τ 1 =G 1 +T2,1 b. It is not difficult to check that [C k,m k,0 k τ 1 ]

A SCALING ANALYSIS OF A CAT AND MOUSE MARKOV CHAIN 13 has the same distribution as the cat and mouse Markov chain during a cycle as defined above. Define t 0 =0 and t i =t i 1 +τ i, s 0 =0 and s i =s i 1 +G i. The i+1th cycle is defined as Ck t a i +s i,ck t a i +s i, t i k<t i +G i+1, C k,m k = Cs a i+1 2I i+1 +I i+1 Ci+1,k t b i G i+1,cs a i+1, t i +G i+1 k t i+1, with I i+1 =Cs a i+1 Cs a i+1 1 and τ i+1 =G i+1 +T2,i+1 b. Thesequence C n,m n has the same distribution as the Markov chain with transition matrix Q defined by relation 1. With this representation, the location M n of the mouse at time n is given by Cκ a n, where κ n is the number of steps the mouse has made up to time n, formally defined as [ + i 1 ] [ def. + i κ n = G l +n t i 1 1 {ti 1 n t i 1 +G i }+ G l ]1 {ti 1 +G i <n<t i }, i=1 l=1 in particular, 10 ν n l=1 ν n+1 G l κ n with ν n defined as the number of cycles of the cat and mouse process up to time n: { l+1 ν n =inf{l:t l+1 >n}=inf l: G k +T2,k }. b >n Define l=1 k=1 G l i=1 { l+1 ν n =inf l: 2+T2,k }, b >n k=1 then, for δ>0, on the event {ν n >ν n +δ n}, ν n+δ n n k=1 ν n+δ n n+ Hence, ν n+1 ν 2+T2,k b n+δ [G k +T2,k b ]+ k=ν n+2 k=1 ν n+1 2+T2,k b G k 2. ν n+δ n k=ν n+2 k=1 n k=ν n+2 ν n+1 2+T2,k b G k 2; k=1 l=1 ν n+1 2+T2,k b G k 2 k=1

14 N. LITVAK AND P. ROBERT since T1,k b 2+Tb 2,k, the relation { 11 {ν n ν n >δ n} { inf 1 l ν n inf 1 l ν n l+ δ n k=l l+ δ n k=l } ν n+1 T1,k b G k 2,ν n >ν n k=1 T b 1,k sup 1 l ν n k=1 } l+1 G k 2 holds. Since EG 1 =2, Donsker s theorem gives the following convergence in distribution: tk +1 1 dist. G k 2,0 t 1 = varg 1 Wt,0 t 1, K lim K + k=1 where Wt is a standard Brownian motion, and, therefore, 12 lim K + For t>0, define 1 sup K 1 l K k=1 n s,0 s t def. = By relation 11 one gets that, for 0 s t, { 13 { n s>δ} { inf 1 l ν ns l+1 G k 2 dist. = varg 1 sup Wt. 0 t 1 1 n ν ns ν ns,0 s t. l+ δ n k=l l+ δ n inf 1 l ν nt k=l T b 1,k sup 1 l ν ns k=1 T b 1,k sup 1 l ν nt k=1 } l+1 G k 2 } l+1 G k 2. Letting ε>0, by Lemma 2 and relation 12, there exist some x 0 >0 and n 0 such that if n n 0, then, respectively, l+1 14 Pν nt x 0 n ε and P sup G k 2 x 0 n ε. By using relation 13, { } n s>δ sup 0 s t {ν nt x 0 n} 1 l x 0 n k=1

A SCALING ANALYSIS OF A CAT AND MOUSE MARKOV CHAIN 15 { l+ δ n inf 1 l ν nt k=l {ν nt x 0 n} { T b 1,k sup 1 l ν nt k=1 l+ δ n inf T1,k b sup 1 l x 0 n k=l } l+1 G k 2,ν nt <x 0 n 1 l x 0 n k=1 } l+1 G k 2. WithasimilardecompositionwiththepartialsumsofG k 2,relations14 give the inequality, for n n 0, k+ δ 1 n P sup n s>δ 2ε+P inf n T1,i b x 0. 0 s t 1 k x 0 n By Lemma 1, the left-hand side is thus arbitrarily small if n is sufficiently large. In a similar way the same results holds for the variable sup n s: 0 s t. The variable sup n s :0 s t converges therefore in distribution to 0. Consequently, by using relation 10 and the law of large numbers, the same property holds for sup 0 s t 1 n κ ns 2ν ns. Donsker stheoremgivesthatthesequenceofprocessesc a ns / 4 n,0 s t converges in distribution to B 1 s,0 s t. In particular, for ε and δ>0, there exists some n 0 such that if n n 0, then 1 P 4 C a n nu Ca nv ε; δ sup 0 u,v t, u v δ see Billingsley [6], for example. Since M n =Cκ a n for any n 1, the processes 1 1 4 M n ns,0 s t and 4 C a n 2ν ns,0 s t have therefore the same asymptotic behavior for the convergence in distribution. Since, by construction C a k and ν n are independent, with Skorohod s representation theorem, one can assume that, on an appropriate probability space with two independent Brownian motions B 1 s and B 2 s, the convergences i=k lim n + Ca ns / 4 n,0 s t=b 1 s,0 s t, lim n + ν ns / n=l B2 s/2,0 s t hold almost surely for the norm of the supremum. This concludes the proof of the theorem.

16 N. LITVAK AND P. ROBERT 3.2. Random walk in the plane. The transition matrix P of this random walk is given by, for x Z 2, px,x+1,0=px,x 1,0=px,x+0,1=px,x 0,1= 1 4. Decomposition into cycles. In the one-dimensional case, when the cat and the mouse start at the same location, when they are separated for the first time, they are at distance 2, so that the next meeting time has the same distribution as the hitting time of 0 for the simple random walk when it starts at 2. For d = 2, because of the geometry, the situation is more complicated. When the cat and the mouse are separated for the first time, there are several possibilities for the patterns of their respective locations and not only one as for d=1. A finite Markov chain has to be introduced that describes the relative position of the mouse with respect to the cat. Let e 1 =1,0, e 1 = e 1, e 2 =0,1, e 2 = e 2 and the set of unit vectors of Z 2 is denoted by E = {e 1,e 1,e 2,e 2 }. Clearly, when the cat and the mouse are at the same site, they stay together a geometric number of steps whose mean is 4/3. When they are just separated, up to a translation, a symmetry or a rotation, if the mouse is at e 1, the cat will be at e 2, e 2 or e 1 with probability 1/3. The next time the cat will meet the mouse corresponds to one of the instants of visit to E by the sequence C n. If one considers only these visits, then, up to a translation, it is not difficult to see that the position of the cat and of the mouse is a Markov chain with transition matrix Q R defined below. Definition 1. Let e 1 =1,0, e 1 = e 1, e 2 =0,1, e 2 = e 2 and the set of unit vectors of Z 2 is denoted by E ={e 1,e 1,e 2,e 2 }. If C n is a random walk in the plane, R n denotes the sequence in E such that R n is the sequence of unit vectors visited by C n and 15 r ef def. = PR 1 =f R 0 =e, e,f E. A transition matrix Q R on E 2 is defined as follows: for e, f, g E, Q R e,g,f,g=r ef, e g, 16 Q R e,e,e, e=1/3, Q R e,e,e,e=q R e,e,e, e=1/3, with the convention that e, e are the unit vectors orthogonal to e, µ R denotes the invariant probability distribution associated to Q R and D E is the diagonal of E 2. A characterization of the matrix R is as follows. Let τ + =infn>0:c n E and τ =infn 0:C n E,

A SCALING ANALYSIS OF A CAT AND MOUSE MARKOV CHAIN 17 then clearly r ef =PC τ + =f C 0 =e. For x Z 2, define φx=pc τ =e 1 C 0 =x. By symmetry, it is easily seen that the coefficients of R can be determined by φ. For x / E, by looking at the state of the Markov chain at time 1, one gets the relation φx def. = φx+e 1 +φx+e 1 +φx+e 2 +φx+e 2 4φx=0 and φe i = 0 if i { 1,2, 2} and φe 1 = 1. In other words, φ is the solution of a discrete Dirichlet problem: it is a harmonic function for the discretelaplacianonz 2 withfixedvaluesone.classically, thereisaunique solution to the Dirichlet problem; see Norris [24], for example. An explicit expression of φ is, apparently, not available. Theorem 3. If C 0,M 0 N 2, the convergence in distribution of finite marginals 1 dist. lim n M e n + nt,t 0 = WZt holds, with 16µR D E Zt= L B T t, 3π where µ R is the probability distribution on E 2 introduced in Definition 1, the process Wt=W 1 t,w 2 t is a two-dimensional Brownian motion and: L B t the local time at 0 of a standard Brownian motion Bt on R independent of Wt. For t 0, T t =inf{s 0:Bs=t}. Proof. The proof follows the same lines as before: a convenient construction of the process to decouple the time scale of the visits of the cat and the motion of the mouse. The arguments which are similar to the ones used in the proof of the one-dimensional case are not repeated. Let R n,s n be the Markov chain with transition matrix Q R that describes the relative positions of the cat and the mouse at the instances of visits of C n,m n to E E up to rotation, symmetry and translation. For N visits to the set E E, the proportion of time the cat and the mouse will have met is given by 1 N N 1 {Rl =S l }; l=1 this quantity converges almost surely to µ R D E.

18 N. LITVAK AND P. ROBERT Now one has to estimate the number of visits of the cat to the set E. Kasahara [18] see also Bingham [7] and Kasahara [17] gives that, for the convergence in distribution of the finite marginals, the following convergence holds: lim n + 1 n e nt i=0 1 {Ci E} dist. 4 = π L BT t. The rest of the proof follows the same lines as in the proof of Theorem 2. Remark. Tanaka s Formula see Rogers and Williams [27] gives the relation Tt LT t =t sgnbsdbs, 0 where sgnx = 1 if x < 0 and +1 otherwise. Since the process T t has independent increments and that the T t s are stopping times, one gets that LT t has also independent increments. With the function t T t being discontinuous, the limiting process WZt is also discontinuous. This is related to the fact that the convergence of processes in the theorem is minimal: it is only for the convergence in distribution of finite marginals. For t 0, the distribution of LT t is an exponential distribution with mean 2t; see Borodin and Salminen [9], for example. The characteristic function of 16µR D E LT t W 1 3π at ξ C such that Reξ=0 can be easily obtained as Ee iξw 1[Zt] = α2 0 α 2 0 +ξ2 t with α 0 = 3π 4 µ R D E. With a simple inversion, one gets that the density of this random variable is a bilateral exponential distribution given by α 0 2 t exp α 0 y, y R. t The characteristic function can be also represented as + Ee iξw1[zt] = α2 0 α 2 =exp e iξu 1Πt,udu 0 +ξt with Πt,u= e α 0 u / t, u R. u Πt,udu is in fact the associated Lévy measure of the nonhomogeneous process with independentincrements W 1 Zt. SeeChapter 5of Gikhman and Skorohod [13].

A SCALING ANALYSIS OF A CAT AND MOUSE MARKOV CHAIN 19 4. The reflected random walk. In this section the cat follows a simple ergodic random walk on the integers with a reflection at 0; an asymptotic analysis of the evolution of the sample paths of the mouse is carried out. Despite being a quite simple example, it exhibits already an interesting scaling behavior. Let P denote the transition matrix of the simple reflected random walk on N, px,x+1=p, x 0, 17 px,x 1=1 p, x 0, p0,0=1 p. It is assumed that p 0,1/2 so that the corresponding Markov chain is positive recurrent and reversible and its invariant probability distribution is a geometric random variable with parameter ρ def. = p/1 p. In this case, one can check that the measure ν on N 2 defined in Theorem 1 is given by νx,y=ρ x 1 ρ, 0 x<y 1, νy 1,y=ρ y 1 1 ρ1 p, νy,y=ρ y 1 ρ, νy+1,y=ρ y+1 1 ρp, νx,y=ρ x 1 ρ, x>y+1. The following proposition describes the scaling for the dynamics of the cat. Proposition 2. If, for n 1, T n =inf{k >0:C k =n}, then, as n goes to infinity, the random variable T n /E 0 T n converges in distribution to an exponentially distributed random variable with parameter 1 and lim E 0T n ρ n = 1+ρ n + 1 ρ 2 with ρ=p/1 p. If C 0 =n, then T 0 /n converges almost surely to 1+ρ/1 ρ. Proof. The first convergence result is standard; see Keilson [19] for closely related results. Note that the Markov chain C n has the same distribution as the embedded Markov chain of the M/M/1 queue with arrival rate p and service rate q. The first part of the proposition is therefore a discrete analogue of the convergence result of Proposition 5.11 of Robert [26]. If C 0 =n and define by induction τ n =0 and, for 0 i n, τ i =inf{k 0:C k+τi+1 =i}, hence, τ n + +τ i is the first time when the cat crosses level i. The strong Markov property gives that the τ i,0 i n 1 are i.i.d. A standard cal-

20 N. LITVAK AND P. ROBERT culation see Grimmett and Stirzaker [15], e.g. gives that Eu τ 1 = 1 1 4pqu, 0 u 1, 2pu hence, Eτ 0 =1+ρ/1 ρ. Since T 0 =τ n 1 + +τ 0, the last part of the proposition is therefore a consequence of the law of large numbers. Additive jumps. An intuitive picture of the main phenomenon is as follows. It is assumed that the mouse is at level n for some n large. If the cat starts at 0, according to the above proposition, it will take of the order of ρ n steps to reach the mouse. The cat and the mouse will then interact for a short amount of time until the cat returns in the neighborhood of 0, leaving the mouse at some new location M. Note that, because n is large, the reflection condition does not play a role for the dynamics of the mouse at this level and by spatial homogeneity outside 0, one has that M =n+m where M is some random variable whose distribution is independent of n. Hence, when the cat has returned to the mouse k times after hitting 0 and then went back to 0 again, the location of the mouse can be represented as n+m 1 + +M k, where M i are i.i.d. with the same distribution as M. Roughly speaking, on the exponential time scale t ρ n t, it will be seen that the successive locations of the mouse can be represented with the random walk associated to M with a negative drift, that is, EM <0. The section is organized as follows: one investigates the properties of the random variable M and the rest of the section is devoted to the proof of the functional limit theorem. The main ingredient is also a decomposition of the sample path of C n,m n into cycles. A cycle starts and ends with the cat at 0 and the mouse is visited at least once by the cat during the cycle. Free process. Let C n,m n be the cat and mouse Markov chain associated to the simple random walk on Z without reflection the free process: p x,x+1=p=1 p x,x 1 x Z. Proposition 3. If C 0,M 0 =0,0, then the asymptotic location of the mouse for the free process M =lim n M n is such that, for u C such that u =1, 18 in particular, Furthermore, the relation 19 Eu M = EM = 1 ρ ρ1 ρu 2 ρ 2 u 2 +1+ρu 1, 1 and E =1. ρ M 1 E sup <+ n 0 ρ M n

A SCALING ANALYSIS OF A CAT AND MOUSE MARKOV CHAIN 21 holds. If S k is the random walk associated to a sequence of i.i.d. random variables with the same distribution as M and E i are i.i.d. exponential random variables with parameter 1 + ρ/1 ρ 2, then the random variable W defined by 20 W = + k=0 ρ S k E k is almost surely finite with infinite expectation. Proof. Let τ =inf{n 1:C n <M n}, then, by looking at the different cases, one has { 1, if M1 =1,C 1 = 1, M τ = 1+M τ, if M 1 =1,C 1 =1, 1+M τ, if M 1 = 1, where M τ is an independent r.v. with the same distribution as M τ. Hence, for u C such that u =1, one gets that Eu M τ= 1 p 1 u +p2 u Eu M τ+p1 pu holds. Since M τ C τ =2, after time τ, the cat and the mouse meet again with probability ρ 2. Consequently, M M τ,i, 1+G dist. = where M τ,i are i.i.d. random variables with the same distribution as M τ and G is an independent geometrically distributed random variable with parameter ρ 2. This identity gives directly the expression 18 for the characteristic function of M and also the relation EM = 1/ρ. Recall that the mouse can move one step up only when it is at the same location as the cat, hence, one gets the upper bound i=1 supm n U def. = 1+sup n 0 n 0 andthefactthat U 1hasthesamedistributionastheinvariantdistribution of the reflected random walk C n, that is, a geometric distribution with parameter ρ gives directly inequality 19. Let N =N t be a Poisson process with rate 1 ρ 2 /1+ρ, then one can check the following identity for the distributions: 21 W dist. = + 0 ρ S N t dt. C n

22 N. LITVAK AND P. ROBERT By the law of large numbers, S Nt /t converges almost surely to 1 + ρ/[1 ρ 2 ρ]. One gets therefore that W is almost surely finite. From 18, one gets u Eu M can be analytically extended to the interval 1+ρ 1 ρ1+3ρ 2ρ 2 <u< 1+ρ+ 1 ρ1+3ρ 2ρ 2 in particular, for u=1/ρ and its value is Eρ M =1. This gives by 20 and Fubini s theorem that EW=+. Note that Eρ M =1 implies that the exponential moment Eu M of the random variable M is finite for u in the interval [1,1/ρ]. Exponential functionals. The representation 21 shows that the variable W is an exponential functional of a compound Poisson process. See Yor [30]. It can be seen as the invariant distribution of the auto-regressive process X n defined as X n+1 def. = ρ An X n +E n, n 0. The distributions of these random variables are investigated in Guillemin et al. [16] when A n are nonnegative. See also Bertoin and Yor [5]. The above proposition shows that W has a heavy-tailed distribution. As it will be seen in the scaling result below, this has a qualitative impact on the asymptotic behavior of the location of the mouse. See Goldie [14] for an analysis of the asymptotic behavior of tail distributions of these random variables. A scaling for the location of the mouse. Therest of the section is devoted to the analysis of the location of the mouse when it is initially far away from the location of the cat. Define and, for k 1, 22 s 1 =inf{l 0:C l =M l } and t 1 =inf{l s 1 :C l =0} s k+1 =inf{l t k :C l =M l } and t k+1 =inf{l s k+1 :C l =0}. Proposition 2 suggests an exponential time scale for a convenient scaling of the location of the mouse. When the mouse is initially at n and the cat at the origin, it takes the duration s 1 of the order of ρ n so that the cat reaches this level. Just after that time, the two processes behave like the free process on Z analyzed above, hence, when the cat returns to the origin at time t 1, the mouse is at position n+m. Note that on the extremely fast exponential time scale t ρ n t, the finite time that the cat and mouse spend together is vanishing, and so is the time needed for the cat to reach zero from n+m linear in n by the second statement of Proposition 2.

A SCALING ANALYSIS OF A CAT AND MOUSE MARKOV CHAIN 23 Hence, on the exponential time scale, s 1 is a finite exponential random variable, and s 2 is distributed as a sum of two i.i.d. copies of s 1. The following proposition presents a precise formulation of this description, in particular, a proof of the corresponding scaling results. For the sake of simplicity, and because of the topological intricacies of convergence in distribution, in a first step the convergence result is restricted on the time interval [0,s 2 ], that is, on the two first cycles. Theorem 4 below gives the full statement of the scaling result. Proposition 4. If M 0 =n 1 and C 0 =0, then, as n goes to infinity, the random variable M t1 n,ρ n t 1 converges in distribution to M,E 1 and the process M tρ n 1 {0 t<ρ n s 2 } n converges in distribution for the Skorohod topology to the process 1 {t<e1 +ρ M E 2 }, where the distribution of M is as defined in Proposition 3, and it is independent of E 1 and E 2, two independent exponential random variables with parameter 1+ρ/1 ρ 2. Proof. For T >0, D[0,T],R denotes the space of cadlag functions, that is, of right continuous functions with left limits, and d 0 is the metric on this space defined by, for x, y D[0,T],R, d 0 x,y= inf ϕ H [ sup 0 s<t<t ϕt ϕs log t s + sup 0 s<t ] xϕs ys, where H is the set of nondecreasing functions ϕ such that ϕ0 = 0 and ϕt=t. See Billingsley [6]. An upper index n is added on the variables s 1, s 2, t 1 to stress the dependence on n. Take three independent Markov chains Ck a, Cb k and Cc k with transition matrix P such that C0 a=cc 0 =0, Cb 0 =n and, for i=a, b, c, Tp i denotes the hitting time of p 0 for Ck i. Since C k,m k,s n 1 k tn 1 has the same distribution as n+c k,n+m k,0 k <Tb 0, by the strong Markov property, the sequence M k,k s n 2 has the same distribution as N k,0 k Tn a +T0 b+tc n, where n, k T n a, n+m k T, n a Ta n k Ta n +Tb 0, 23 N k = n+m, T a T0 b n +T0 b k Ta n +T0 b+tc n+m. T 0 b

24 N. LITVAK AND P. ROBERT Here Ck b n,m k,0 k Tb 0 is a sequence with the same distribution as the free process with initial starting point 0,0 and killed at the hitting time of n by the first coordinate. Additionally, it is independent of the Markov chains Ck a and Cc k. In particular, the random variable M t 1 n, the jump of the mouse from its initial position when the cat hits 0, has the same distribution as M. Since T b T0 b 0 converges almost surely to infinity, M T0 b is converging in distribution to M. Proposition 2 and the independence of Ck a and Cc k show that the sequences ρ n Tn a and ρ n Tn c converge in distribution to two independent exponential random variables E 1 and E 2 with parameter 1+ρ/1 ρ 2. By using Skorohod s Representation theorem,see Billingsley [6] up to a change of probability space, it can be assumed that these convergences hold for the almost sure convergence. By representation 23, the rescaled process M tρ n /n1 {0 t<ρ n s 2 }, t T has the same distribution as 1, t<ρ n Tn, a x n t def. = 1+ 1 n M ρ n t Tn, a ρn Tn a t<ρ n Tn a +T0 b, 1+ 1 n M, ρ n T a T0 b n +Tb 0 t<ρn Tn a+tb 0 +Tc n+m T 0, b 0, t ρ n Tn a +T0 b+tc n+m T 0, b for t T. Proposition 2 shows that T0 b /n converges almost surely to 1 ρ/ 1+ρ so that ρ n Tn a +T0 b converges to E 1 and, for n 1, ρ n T c n+m T b 0 =ρ M T b 0ρ n+m T b 0T c n+m T b 0 ρ M E 2, almost surely as n goes to infinity. Additionally, one has also almost surely. Define 1 lim n + n sup M k =0, k 0 x =1 {t<t E1 +ρ M E 2 }, where a b=mina,b for a, b R. Time change. For n 1 and t>0, define u n resp., v n as the minimum resp., maximum of t ρ n T a n +Tb 0 +Tc n+m T b 0 and t E 1 +ρ M E2, and ϕ n s= v n u n s, 0 s u n, v n +s u n T v n T u n, u n <s T.

A SCALING ANALYSIS OF A CAT AND MOUSE MARKOV CHAIN 25 Notingthat ϕ n H,byusingthisfunctioninthedefinitionofthedistanced 0 on D[0,T],R to have an upper bound of dx n,x and with the above convergence results, one gets that, almost surely, the sequence dx n,x converges to 0. The proposition is proved. Theorem 4 Scaling for the location of the mouse. If M 0 =n, C 0 =0, then the process M tρ n 1 {t<ρ n t n} n converges in distribution for the Skorohod topology to the process 1 {t<w}, where W is the random variable defined by 20. If H 0 is the hitting time of 0 by M n, H 0 =inf{s 0:M s =0}, then, as n goes to infinity, ρ n H 0 converges in distribution to W. Proof. In the same way as in the proof of Proposition 4, it can be proved that for p 1, the random vector [M tk n,ρ n t k,1 k p] converges in distribution to the vector k 1 S k, ρ S i E i and, for k 0, the convergence in distribution M tρ 24 lim n 1 n + n {0 t<ρ n t k } =1 {t<e1 +ρ S 1E 2 + +ρ S k 1E k } i=0 holds for the Skorohod topology. Let φ:[0,1] R + be defined by φs=eρ sm, then φ0 =φ1=1 and φ 0<0, since φ is strictly convex then for all s<1, φs<1. If C 0 =M 0 =n, and the samplepath of M k n,k 0 follows the sample path of areflected randomwalk startingat 0, wehave, inparticular, that the supremum of its successive values is integrable. By Proposition 3, as n goes to infinity, M t1 n is converging in distribution to M. Lebesgue s theorem gives therefore that the averages are also converging, hence, since EM is negative, there exists N 0 such that if n N 0, 25 E n,n M t1 def. = EM t1 M 0 =C 0 =n n+ 1 2 EM =n 1 2ρ. Note that t 1 has the same distribution as T 0 in Proposition 2 when C 0 =n. Proposition 2 now implies that there exists K 0 0 so that, for n N 0, 26 ρ n/2 E 0,n t 1 K 0.

26 N. LITVAK AND P. ROBERT The identity E1/ρ M =1 implies that Eρ M /2 <1, and inequality 19 and Lebesgue s theorem imply that one can choose 0<δ<1 and N 0, so that 27 Eρ n Mt 1 /2 δ holds for n N 0. Let ν =inf{k 1:M tk N 0 } and, for k 1, G k the σ-field generated by the random variables C j,m j for j t k. Because of inequality 25, one can check that the sequence M tk ν + 1 2ρ k ν,k 0 is a super-martingale with respect to the filtration G k, hence, EM tk ν + 1 2ρ Ek ν EM 0=n. Since the location of the mouse is nonnegative, by letting k go to infinity, one gets that Eν 2ρn. In particular, ν is almost surely a finite random variable. Intuitively, t ν is the time when the mouse reaches the area below a finite boundary N 0. Our goal now is to prove that the sequence ρ n t ν converges in distribution to W. For p 1 and on the event {ν p}, 28 ρ n t ν t p 1/2 = ν 1 1/2 ν 1 ρ n t k+1 t k ρ n t k+1 t k. k=p For k p, inequality 26 and the strong Markov property give that the relation k=p ρ Mt k /2 E[ t k+1 t k G k ]=ρ Mt k /2 E 0,Mtk [ t 1 ] K 0 holds on the event {ν >k} {M tk >N 0 }. One gets therefore that E ρ n t k+1 t k 1 {k<ν} =Eρ n Mt k /2 1 {k<ν} ρ Mt k /2 E[ t k+1 t k G k ] K 0 Eρ n Mt k /2 1 {k<ν} holds, and, with inequality 27 and again the strong Markov property, Eρ n Mt k /2 1 {k<ν} =Eρ k 1 j=0 Mt j+1 Mt j /2 1 {k<ν} δeρ k 2 j=0 Mt j+1 Mt j /2 1 {k 1<ν} δ k. Relation 28 gives therefore that E ρ n t ν t p K 0δ p 1 δ.

For ξ 0, 29 A SCALING ANALYSIS OF A CAT AND MOUSE MARKOV CHAIN 27 Ee ξρn t ν Ee ξρn t p E1 e ξρn t ν t p + +Pν <p = + 0 K 0δ p 1 δ ξe ξu Pρ n t ν t p udu+pν <p + 0 ξ u e ξu du+pν <p by using Markov s inequality for the random variable ρ n t ν t p. Since ρ n t p converges in distribution to E 0 + ρ S 1 E 1 + + ρ Sp E p, one can prove that, for ε>0, by choosing a fixed p sufficiently large and that if n is large enough, then the Laplace transforms at ξ 0 of the random variables ρ n t ν and W are at a distance less than ε. At time t ν the location M tν of the mouse is x N 0 and the cat is at 0. Since the sites visited by M n are a Markov chain with transition matrix px,y, with probability 1, thenumber R of jumpsfor themouseto reach 0 is finite. By recurrence of C n, almost surely, the cat will meet the mouse R times in a finite time. Consequently, if H 0 is the time when the mouse hits 0 for the first time, then by the strong Markov property, the difference H 0 t ν is almost surely a finite random variable. The convergence in distribution of ρ n H 0 to W is therefore proved. Nonconvergence of scaled process after W. Theorem 4 could suggest that the convergence holds for a whole time axis, that is, M tρ lim n,t 0 =1 n + n {t<w},t 0 for the Skorohod topology. That is, after time W the rescaled process stays at 0 like for fluid limits of stable stochastic systems. However, it turns out that this convergence does not hold at all for the following intuitive and nonrigorous reason. Each time the cat meets the mouse at x large, the location of the mouse is at x+m when the cat returns to 0, where M is the random variable defined in Proposition 3. In this way, after the kth visit of the cat, the mouse is at the kth position of a random walk associated to M starting at x. Since E1/ρ M =1, Kingman s result see Kingman [21] implies that the hitting time of δn, with 0<δ<1, by this random walk started at 0 is of the order of ρ δn. For each of the steps of the random walk, the cat needs also of the order of ρ δn units of time. Hence, the mouse reaches the level δn in order of ρ 2δn steps, and this happens on any finite interval [s,t] on the time scale t ρ n t only if δ 1/2. Thus, it is very likely that the next relation holds: lim P n + sup s u t M uρ n n = 1 2 =1.