Abstract:
Let $\{S_n,\,n\geqslant 0\}$ be a random walk with increments that belong (without centering) to the domain of attraction of an $alpha$-stable law, that is, there exists a process $\{Y_t,\,t\geqslant 0\}$ such that $S_{nt}/a_{n}$$\Rightarrow$$Y_t$, $t\geqslant 0$, as $n\to\infty$ for some scaling constants $a_n$. Assuming that $S_{0}=o(a_n)$ and $S_n\leqslant \varphi (n)=o(a_n)$, we prove several conditional limit theorems for the distribution of the random variable $S_{n-m}$ given that $m=o(n)$ and $\min_{0\leqslant k\leqslant n}S_k\geqslant 0$. These theorems supplement the assertions established by Caravenna and Chaumont in 2013. Our results are used to study the population size of a critical branching process evolving in an unfavourable environment.
Bibliography: 28 titles.
Keywords:
random walks, stable law, conditional limit theorems, branching processes, unfavourable random environment.
Ministry of Science and Technology (MOST) of China
G20221740071
The research of V. A. Vatutin was performed at the Steklov International Mathematical Center and supported by the Ministry of Science and Higher Education of the Russian Federation (agreement no. 075-15-2022-265) and also by the Ministry of Science and Technology of the People's Republic of China (grant no. G20221740071). The research of C. Dong was supported by the Ministry of Science and Technology of the People's Republic of China (grant no. G20221740071).
The research of E. E. Dyakonova was performed at the Steklov International Mathematical Center and supported by the Ministry of Science and Higher Education of the Russian Federation (agreement no. 075-15-2022-265).
The investigations of the conditional distributions of random walks conditioned to stay positive or nonnegative have a long history (see, for example, [4], [6]–[12], [14], [17], [18], [20] and [25]). The interest in various invariance principles for conditional processes, to which most of these works are devoted, is due not only to the natural development of the theory of random walks, but also to the widespread use of such results in the theory of branching processes evolving in constant and random environments alike (see [2], [19] and [26]), in statistical physics, in particular, in investigations of random polymers (see [13]), and also in other areas.
In this paper we also analyze the properties of random walks conditioned to stay nonnegative. Our research was motivated by the results due to Caravenna and Chaumont [11]. They showed that if the distribution of the step of a random walk belongs (without centering) to the domain of attraction of a stable law, then, provided that this walk stays nonnegative on an interval of length $n$, the distribution of the appropriately normalized excursion of length $n$ generated by this random walk converges to the distribution of the excursion of a stable Lévy process conditioned to stay nonnegative on the interval $[0,1]$. We supplement this result by a study of the behaviour of the appropriately rescaled trajectory of the random walk in a left neighbourhood of the $n-m$ excursion, where $m=o(n)$. It turns out that three different distributions occur in the limit, depending on the rate of convergence of $m/n$ to zero. These three regimes lead to three different rates of growth for the number of particles in critical branching processes evolving in an unfavourable random environment and conditioned to survive by a distant moment $n$. Our assertions for critical branching processes in random environment supplement the corresponding theorems established in [26], [27] and [28], where the distributions of the number of particles at times $m=o(n)$ and $m=[nt]$, $0<t\leqslant 1$, were considered under the condition that the processes survive by time $n$.
The paper is structured as follows.
In § 2 we introduce the main concepts, describe our basic conditions imposed on random walks, recall some known local limit theorems, and prove auxiliary assertions for random walks conditioned to stay nonnegative.
In § 3 several conditional limit theorems are presented for excursions of a lattice random walk conditioned to stay nonnegative.
In § 4 we prove the corresponding conditional limit theorems for excursions of a random walk whose step has an absolutely continuous distribution.
In § 5 we establish a conditional limit theorem for an almost surely convergent sequence of random variables.
In § 6 we apply our results for random walks conditioned to stay nonnegative to the description of the distribution of the population size of a critical branching process evolving in an unfavourable environment and conditioned to survive by a distant moment of time.
§ 2. Notation and assumptions
For the unification of presentation and convenience of references we adhere mainly to the notation and assumptions introduced in [11].
In what follows we let $C_{1},C_{2},\dots$ denote some absolute constants, which do not necessarily coincide in distinct formulae.
We introduce the sets $\mathbb{N}\mathbf:=\{1,2,\dots\} $ and $\mathbb{N}_{0}:=\mathbb{N} \cup \{0\}$. Given two positive sequences $\{c_{n},\,n\in \mathbb{N}\} $ and $\{d_{n},\,n\in\mathbb{N}\}$, we write, as usual,
Recall that a positive sequence $\{c_{n},\,n\in \mathbb{N}\} $ (or a real function $c(x))$ is called regularly varying at infinity with index $\gamma \in \mathbb{R}$, which is denoted by $c_{n}\in R_{\gamma}$ or $c(x)\in R_{\gamma}$, respectively, if $c_{n}\sim n^{\gamma}l(n) $ ($c(x)\sim x^{\gamma}l(x)$, respectively), where $l(x)$ is a slowly varying function, that is, a positive function such that $l(cx)/l(x)\to 1$ as ${x\to \infty}$ for any fixed $c>0$.
be a subset of the space $\mathbb{R}^{2}$. Given a pair $(\alpha,\beta)\in \mathcal{A}$ and a random variable $X$, we write $X\in \mathcal{D}(\alpha,\beta)$ if the distribution of $X$ belongs (without centering) to the domain of attraction of a stable law with density $g(x)=g_{\alpha,\beta }(x)$, $x\in (-\infty,+\infty)$, and the characteristic function is
with independent identically distributed increments. Throughout this section we assume that the random walk $(\mathcal{S}=\{S_{n},\,n\in \mathbb{N}_{0}\},\,\mathbb{P})$ satisfies the following conditions.
Condition A1. The increments $X_{n}$, $n=1,2,\dots$, of $\mathcal{S}$ belong to $\mathcal D(\alpha,\beta)$.
This means, in particular, that there is an increasing sequence of positive numbers
as $n\to \infty $; the symbol $\Longrightarrow $ denotes convergence in distribution in the space $D[0,+\infty)$ with the Skorokhod topology, and the process $(\mathcal{Y}=\{Y_t,\,t\geqslant 0\},\,\mathbf{P})$ is strictly stable, that is, has stationary independent increments and marginal distributions described by the characteristic functions
Let $\Omega^{\mathrm{RW}}:=\mathbb{R}^{\mathbb{N}_{0}}$ be the space of right-continuous functions with jumps at integer points on the nonnegative half-axis, and let $\Omega :=D([0,\infty),\mathbb{R}) $ be the space of right-continuous real functions on $[0,\infty)$ that have finite limits from the left. The space $\Omega $ endowed with the Skorokhod topology turns to a Polish space with the corresponding $\sigma$-algebra of Borel sets. We also set $\Omega_{N}^{\mathrm{RW}}:=\mathbb{R}^{\{0,1,\dots,N\}}$ and $\Omega_t:=D([0,t],\mathbb{R})$ for $t\in (0,\infty)$.
We let $\mathbb{P}_x$ denote the law of the random walk starting at the point $x\in \mathbb{R}$, that is, the probability law on $\Omega^{\mathrm{RW}}$ generated by the random walk $\mathcal{S}+x$, where $\mathbb{P}$ is the law of a random walk $\mathcal{S}$ starting at zero; given a process $\mathcal{Y}=\{Y_t,\,t\geqslant 0\} $ with distribution $\mathbf{P}$ on $\Omega $, we let $\mathbf{P}_a$ denote the distribution of the shifted process $\mathcal{Y}+a$ for $a\in \mathbb{R}$.
Condition A2. One of the following constraints holds.
$\bullet$ (The $(h;c)$-lattice case.) The measure $\mathbb{P}$ of the increment $X_1$ of a random walk is concentrated on the lattice $c+hZ$, where the step $h>0$ is chosen to be maximum possible (that is, the distribution of $X_1$ is not concentrated on any lattice $c_{0}+h_{0}Z$ for $h_{0}>h$ and $c_{0}\in R$). Note that we can take $c\in [0,h)$.
$\bullet$ (The absolutely continuous case.) The measure $\mathbb{P}$ of the increment $X_1$ of a random walk is absolutely continuous with respect to the Lebesgue measure on $\mathbb{R}$, and there exists $n\in \mathbb N$ such that the density $f_{n}(x):=\mathbb{P}(S_{n}{\in}dx){/}dx$ of the random variable $S_{n}$ is bounded (and therefore $f_{n}(x){\in}L^{\infty}$).
The law of the bridge of length $N\in \mathbb{N}$ of a positive random walk starting at the point $x\in [0,\infty)$ and ending at $y\in [0,\infty)$ is the probability law on the set $\Omega_{N}^{\mathrm{RW}}$ defined by
where $\mathcal{K}(N-1):=\{s_{1}>0,\dots,s_{N-1}>0\} $ is the domain of integration and $f(\,\cdot\,)=f_{1}(\,\cdot\,)$ is the density of the random walk increment $X_1$.
For $t\in [0,\infty)$ and $a,b\in [0,\infty)$ we let $\mathbf{P}_{a,b}^{\uparrow,t}$ denote the law on the set $\Omega_{t}$ corresponding to the bridge of length $t$ of a nonnegative Lévy process starting at $a$ and ending at $b$, which can informally be defined by
(see § 6.1 in [11], where it is explained in detail that such a law is well defined).
In the case when $\alpha =2$ and $\rho =1/2$, that is, when the process $\mathcal{Y}$ is a standard Brownian motion, the law $\mathbf{P}_{a,b}^{\uparrow,1}(\,\cdot\,) $ specifies the distribution of the Brownian excursion.
For $N\in \mathbb{N}$ we define the rescaling map $\phi_{N}\colon \Omega _{N}^{\mathrm{RW}}\to \Omega_{1}$ by
where $\{a_{N},\,N\in \mathbb{N}\} $ is the scaling sequence appearing in Condition A1.
Using this definition, we let $\mathbb{P}_{x,y}^{\uparrow,N}\,{\circ}\, \phi_{N}^{-1}$ denote the law induced on the set $\Omega_{1}:=D([0,1],\mathbb{R}) $ via $\mathbb{P}_{x,y}^{\uparrow,N}$ and $\phi_{N}$.
The following important theorem was established in [11].
Theorem 1. Let $a,b\in [0,\infty)$, and let $\{x_{N},\,N\in \mathbb{N}\} $ and $\{y_{N},\,N\in \mathbb{N}\} $ be two nonnegative sequences such that $x_{N}/a_{N}\to a$ and $y_{N}/a_{N}\to b$ (in the $(h;c)$-lattice case, we assume that $(y_{N}-x_{N}) \in Nc+h\mathbb{Z}$ for all $N\in \mathbb{N}$). If Conditions A1 and A2 hold, then
for any $z>0$. Therefore, $S_{N-m}/a_{N}\to 0$ in probability as $N\to \infty $; thus, Theorem 1 gives little information about the distribution of the random variable $S_{N-m}$ in this situation.
In this section we indicate a centering and a rescaling of the random variable $S_{N-m}$ that ensure the convergence of the transformed random variable to a proper nondegenerate distribution. It turns out that the form of the limit distribution depends substantially on the rate of change of the parameter $y_{N}$. In addition, to include also the absolutely continuous case into consideration, we will show that
where $A(z)$ is a proper nondegenerate distribution, the form of which is different in the cases when $y_{N}/a_{m}\to 0$, $y_{N}/a_{m}\to T\in (0,\infty)$ and $y_{N}/a_{m}\to \infty $.
To prove the results stated we need new notation, in which the infimum over an empty set is infinite.
We set $S_{0}:=0$ and $\tau _{0}^{\pm}:=0$; for $k\geqslant 1$ we let
Furthermore, if Condition A1 holds, then $\lim_{n\to\infty}\mathbf{P}(S_n>0)=\rho\in (0,1)$; therefore, the random variables $\tau_1^+$ and $\tau_1^-$ are proper and (see, for example, [21] or [23]) we have
Note that by [25], relations (15) and (31), and [11], relation (3.18), there are positive constants $\widehat{C}$, $\mathcal{C}^{+}$, and $\mathcal{C}^{-}$ such that
We will repeatedly use this fact referring to the results obtained in [14] and [25].
In what follows we need various functions of the form $T(x_{N}$, $y_{N})$ depending on the real parameters $x_{N}$ and $y_{N}$.
Given a positive sequence $a_{n}$, $n=1,2,\dots$, for brevity we will write
“$T(x_{N},y_{N})=o(1)$ or $T(x_{N},y_{N})\sim c>0$ uniformly with respect to $x_{N}=o(a_{N})$, $y_{N}=o(a_{m})$ as $\min(m,N)\to \infty$”,
to indicate that
“$T(x_{N},y_{N})=o(1)$ or $T(x_{N},y_{N})\sim c>0$ uniformly with respect to $x_{N}\in (0,\delta_{N}a_{N}]$, $y_{N}\in (0,\delta_{N}a_{m}]$ for any positive sequence $\delta_{N}\to 0$ as $\min(m,N)\to \infty$”.
To simplify the presentation below we recall several local limit theorems established in [11], § 4, and valid in the lattice case.
Let $g(\,\cdot\,)$ be the density of a random variable $Y_{1}$ and $g^{+}(\,\cdot\,)$ be the density of the meander of a Lévy process $\mathcal{Y}$ at time $t=1$ (see [8]):
(see Lemma 4 in [12], where the properties of this density were described in detail). Similarly, let $g^{-}(\,\cdot\,)$ be the density of the meander of the Levy process $-\mathcal{Y}$ at time $t=1$. Finally, for $a,b\in [0,\infty)$ we introduce the function
Lemma 1. Assume that Conditions A1 and A2 (in the $(h;c)$-lattice case) are satisfied. Then for $x,y\geqslant 0$, where $(y-x) \in nc+h\mathbb{Z}$, the following relations hold as ${n\to \infty}$:
(1) uniformly with respect to $x=o(a_{n})$ and $y\geqslant 0$,
Note that (2.12) is a consequence of Liggett’s invariance principle for bridges (see [20]) and Gnedenko’s local theorem (see [16]).
The renewal functions introduced above were defined in terms of characteristics of the weak ladder moments of random walks. It is worth noting that both the strict ladder moments $\{\widehat{\tau}_{k}^{\pm},\,k\geqslant 0\} $ and the quantities $\{\widehat{H}_{k}^{\pm},\,k\geqslant 0\} $ defined by the relations $\widehat{\tau}_{0}^{\pm}:=0$, $\widehat{H}_{0}^{\pm}:=0$ and
for $k\geqslant 1$ are also often used to study the properties of random walks and branching processes in a random environment (see, for example, [1], [3], [14], [25] and [27]). The sequences $\{\widehat{H}_{k}^{\pm},\,k\geqslant 0\} $ generate the renewal functions
In the absolutely continuous case results similar to Lemma 1 were also deduced in [11].
Lemma 2. Assume that Conditions A1 and A2 (in the absolutely continuous case) are satisfied. Then the following relations hold for $x,y\geqslant 0$ as $n\to \infty$:
(1) uniformly with respect to $x=o(a_{n})$ and $y\geqslant 0$,
Relation (2.16) was not explicitly indicated in [11]. However, it is straightforward to derive it from (2.15) and (2.2) using symmetry, or, more precisely, by considering the random walk $-\mathcal{S}$ instead of $\mathcal{S}$, swapping $x$ and $y$, and replacing each symbol ‘$+$’ by ‘$-$’.
Relation (2.18) is a consequence of Liggett’s invariance principle for bridges (see [20]) and the local theorem proved by Stone in [24].
Remark 1. Note that the value of the constant $C^{\ast}$ is not universal in the following sense. If we replace the scaling sequence by $\overline{a}_{n}=ca_{n}$ for some $c>0$, then
Furthermore, let $U(w):=C^{\ast \ast}w^{\alpha (1-\rho)},w\geqslant 0$. It follows from the proof of Theorem 1.1 in [10] (see the reasoning between (3.11) and (3.12) and the definition (3.1) there) that
Considering the distribution of the meander of the process $-\mathcal{Y}$ at $t=1$ and noting that the positivity parameter of the process $-\mathcal{Y}$ is $1-\rho $, we can easily derive from this fact using symmetry considerations that
The second and third assertions in the lemma can be deduced using summation and integration, respectively.
The lemma is proved.
§ 3. Limit theorems for the lattice case
Following [11], we consider below only the $(1;0)$-lattice case of Condition A2, that is, we assume that the distribution of $X_{1}$ is concentrated on the set $\mathbb{Z}$ and is aperiodic. The general $(h;c)$-lattice case requires only more cumbersome notation, while the proofs need no new ideas.
We aim at studying the asymptotic behaviour of the probabilities
under the condition that $\max (x_{N},y_{N})/a_{N}\to 0$ and $m=o(N)$ as $N\to \infty $.
We consider three cases separately: $y_{N}/a_{m}\to 0$, $y_{N}/a_{m}\to T\in (0,\infty)$ and $y_{N}/a_{m}\to \infty $.
3.1. The case when $y_{N}/a_{m}\to 0$
Lemma 5. Assume that Conditions A1 and A2 are satisfied and the random variable $X_{1}$ has a $(1;0)$-lattice distribution. If $\min (m,N)\to \infty $ so that $m=o(N)$, then
uniformly with respect to $k$, $\varepsilon a_{m}\leqslant k\leqslant za_{m}$, $y_{N}=o(a_{m})$, and $x_{N}=o(a_{N})$. Using (2.6) we infer the relation
as $N\to \infty $ uniformly with respect to $y_{N}=o(a_{m})$ and $x_{N}=o(a_{N})$. Since $\varepsilon >0$ can be chosen to be arbitrarily small, based on (3.6) and (3.2) we deduce that
as $N\to\infty$ uniformly with respect to $x_{N}=o(a_{N})$ and $y_{N}=o(a_{m})$.
The lemma is proved.
Corollary 1. Assume that Conditions A1 and A2 are satisfied and the random variable $X_{1}$ has a $(1;0)$-lattice distribution. If $\min (m,N)\to \infty $ so that $m=o(N)$, then
for any $z\in (0,\infty)$ uniformly with respect to $x_{N}=o(a_{N})$ and $y_{N}=o(a_{m})$.
Using summation, we derive the required assertion from (3.9), (3.2) and (3.3).
3.2. The case when $y_{N}\asymp a_{m}$
Let $0<t_{0}<t_{1}<\infty $ be fixed positive numbers.
Lemma 6. Assume that Conditions A1 and A2 are satisfied and the random variable $X_{1}$ has a $(1;0)$-lattice distribution. If $\min (m,N)\to \infty $ so that $m=o(N)$, then
as $m\to \infty $ uniformly with respect to $\varepsilon a_{m}\leqslant k\leqslant za_{m}$ and $j\in [ t_{0}a_{m},t_{1}a_{m}]$. In view of (3.7) and (3.10) it follows that
Unlike §§ 3.1 and 3.2, where the distribution of the random variable $S_{N-m}$ was considered, here we look at the distribution of the difference $S_{N-m}-S_{N}$.
Lemma 7. Assume that Conditions A1 and A2 are satisfied and the random variable $X_{1}$ has a $(1;0)$-lattice distribution. If $\min (m,N)\to \infty $ and $a_{m}=o(y_{N})$, then
according to the invariance principle for random walks such that the distributions of their increments belong (without centering) to the domain of attraction of a stable law. Note that $b_{N-m}\sim b_{N}$ if $m=o(N)$ as $N\to \infty $; therefore,
as $N\to \infty $. Since we can choose the parameter $M$ to be arbitrarily large and $\varepsilon (M)>0$ in (3.16) to be arbitrarily small, it is true that
uniformly with respect to $x_{N}=o(a_{N})$ and $y_{N}=o(a_{N})$. In view of (2.1) this implies the assertion of the lemma.
Corollary 3. Assume that Conditions A1 and A2 are satisfied and the random variable $X_{1}$ has a $(1;0)$-lattice distribution. If $m\to \infty $ and $a_{m}=o(y_{N})$, then
Proof. The proofs of assertions (1)–(3) repeat conceptually the arguments used for the corresponding assertions in Corollaries 1–3. Therefore, we only verify (4.1).
as $N\to \infty $ uniformly with respect to $x_{N}=o(a_{N})$ and $y_{N}\sim Ta_{m}=o(a_{N})$.
Relation (4.1) is proved, and the proof of Lemma 8 is complete.
§ 5. Limit theorem for almost surely convergent sequences
In this section we need a conditional measure $\mathbb{P}_x^{\uparrow}(\,\cdot\,)$ generating a random walk starting at $x$ and conditioned to stay nonnegative on the whole time axis (see, for example, [4], [10] and [25]). This new measure is defined as follows: for $x\geqslant 0$, any $N\in \mathbb{N}$ and each measurable set $B$ in the $\sigma $-algebra generated by the random variables $S_{1},\dots,S_{N}$ we set
Theorem 2. Assume that Conditions A1 and A2 are satisfied and $H_{1},H_{2},\dots$ is a uniformly bounded sequence of random variables adapted to a filtration $\widetilde{\mathcal{F}} = \{\widetilde{\mathcal{F}}_{k},\,k\in\mathbb{N}\}$ that converges $\mathbb{P}_{0}^{\uparrow}$-almost surely to a random variable $H_{\infty}$ as ${n\to \infty}$. Assume that the parameter $m=m(n)$ tends to infinity as $n\to \infty $ so that $m=o(n)$. Then the following hold:
where $\mathcal{S}'=\{S_{n}',\,n=0,1,2,\dots\} $ is a probabilistic copy of the random walk $\mathcal{S}$ that is independent of the set $\{S_{j},\,j=0,1,\dots,k\} $. Owing to Corollaries 1 and 2 and relation (2.21), there are constants $C$, $C_{1}$, $C_{2}$ and $C_{3}$ such that
First letting $n$ and then $k$ tend to infinity and applying the dominated convergence theorem, we see that the right-hand side of the previous relation tends to zero for each $\lambda >1$.
Based on this result, we arrive at the chain of equalities
Without loss of generality let $H_{\infty}>0$ and $A_{1}(z\lambda^{-1/\alpha}) \, \mathbb{E}_{0}^{\uparrow}[H_{\infty}] \leqslant 1$; in view of this assumption, we conclude that
Combining the above estimates proves (5.1) in the $(1;0)$-lattice case.
To prove (5.2) for the $(1;0)$-lattice case we must repeat almost verbatim the arguments used in the verification of (5.1), replacing the reference to Corollary 1 by a reference to Corollary 2.
To verify (5.3) in the $(1;0)$-lattice case we repeat almost verbatim the arguments used to prove assertion (5.1), replacing the reference to Corollary 1 by a reference to Corollary 3.
To prove assertion (5.1) in the absolutely continuous case we replace $\sum_{j=0}^{\varphi (n)}V^{+}(j)$ in the above reasoning by $\displaystyle\int_{0}^{\varphi (n)}V^{+}(w)\,dw$ and use Lemma 8. Similar arguments also make it possible to verify (5.2) and (5.3) in the case when the random variable $X_{1}$ has an absolutely continuous distribution.
In this section we use the results obtained for random walks to study the population size of critical branching processes evolving in an unfavourable random environment. To describe the model under consideration and state the problems to be solved in detail, we introduce the space $\mathcal{M}$ $=\{\mathfrak{f}\} $ of all probability measures on the set $\mathbb{N}_{0}$. To simplify our notation we identify a measure $\mathfrak{f}=\{\mathfrak{f}(\{0\}),\mathfrak{f}(\{ 1\}),\dots\} \in \mathcal{M}$ with the corresponding probability generating function
and do not distinguish between $\mathfrak{f}$ and $f$. The space $\mathcal{M}$ $=\{\mathfrak{f}\} =\{f\} $ endowed with a metric induced by the distance in variation between distributions turns into a Polish space. Let
be a sequence of independent probabilistic copies of the random variable $F$. The infinite sequence $\mathcal{E}=\{F_{n},\,n\in \mathbb{N}\} $ is called a random environment.
A sequence of nonnegative integer-valued random variables $\mathcal{Z}=\{Z_{n},\,n\in \mathbb{N}_{0}\} $ defined on some probability space $(\Omega,\mathcal{F},\mathbb{P})$ is called a branching process in a random environment (BPRE) if $Z_{0}$ is independent of $\mathcal{E}$ and, under the condition that $\mathcal{E}$ is fixed, the process $\mathcal{Z}$ is a Markov chain:
for all $n\in \mathbb{N}$, $z_{n-1}\in \mathbb{N}_{0}$, and $f_{1},f_{2},\ldots\in \mathcal{M}$, where $\xi_{n1},\xi_{n2},\dots $ is a sequence of independent identically distributed random variables with distribution $f_{n}$. Thus, $Z_{n-1}$ is the size of the $(n-1)$st generation of the branching process and $f_{n}$ is the distribution of the number of offspring of each particle in the $(n-1)$st generation.
It is known (see [1], Theorem 1.1 and Corollary 1.2) that if Conditions B1 and B2 are satisfied, then there are a constant $\theta \in (0,\infty)$ and a sequence $l(1),l(2),\dots$ slowly varying at infinity such that
for any $t\in [0,1] $ and $x\geqslant 0$, where $\mathcal{Y}^{+}=\{Y_{t}^{+},\,0\leqslant t\leqslant1\} $ denotes the meander of a strictly stable process $\mathcal{Y}$ of index $\alpha$.
Thus, if a BPRE is critical, then the random variables $\log Z_{[nt]}$, $t\in (0,1]$, under the condition $Z_{n}>0$ and $S_{n}$ (the value of the associated random walk ensuring the survival of the population by the moment $n$) grow as $a_{n}$ multiplied by random factors.
These results were supplemented in [28] by an analysis of the distributions of the appropriately rescaled random variables $\log Z_{[nt]}$, $t\in (0,1]$, considered under the condition that $Z_{n}>0$ and $S_{n}\leqslant \varphi (n)$, where $\varphi (n)\to \infty $ as $n\to \infty$ so that $\varphi (n)=o(a_{n})$. In this case, in view of (6.1), the event $\{S_{n}\leqslant \varphi(n)\} $ can be interpreted as unfavourable for the evolution of the critical branching process in random environment conditioned to survive by a distant moment of time.
We introduce the notation
$$
\begin{equation*}
A_{\mathrm{u.s}}:=\{Z_{n}>0\text{ for all }n>0\}
\end{equation*}
\notag
$$
for the event meaning that the process survives throughout its existence and set
Note that, due to [27], Theorem 1, if Conditions B1 and B2 are satisfied, the distribution of the increments of the associated random walk is absolutely continuous, and $\varphi(n)=o(a_{n})$, then
A similar asymptotics with $\displaystyle\int_{0}^{\varphi (n)}V^{+}(w)\,dw$ replaced by $\sum_{j=0}^{\varphi (n)}V^{+}(j)$ takes place for the $(1;0)$-lattice case.
Along with the asymptotic behaviour of the survival probability of a critical BPRE, the following theorem describes the rate of growth of the population size on the logarithmic scale.
Theorem 3 (see [28]). Assume that Conditions B1 and B2 are satisfied. If $\varphi (n)\to \infty $ as $n\to \infty $ so that $\varphi (n)=o(a_{n})$, then
for any $t\in (0,1) $ and $x\in [0,\infty)$, where $\mathcal{Y}^{++}=\{Y_{t}^{++},\,0\leqslant t\leqslant1\} $ is the excursion of a strictly stable process $\mathcal{Y}$ of index $\alpha $.
There is a substantial difference between the scaling orders in (6.4) and (6.5), which shows that there must be a phase transition in the rate of growth of the population size in the case when we consider the process inside an interval $[n-m,n)$, where $m=o(n)$ as $m\to \infty $. Such a transition actually takes place, and the aim of this section is to prove that three distinct limit distributions occur when the random variables $\log Z_{n-m}$ are rescaled by $a_{m}$. The form of these distributions depends on which of the following three conditions holds as $\min (n-m,m)\to \infty $: the function $\varphi(n)$ is of order $o(m)$; the function $\varphi (n)$ is proportional to $m$; $m=o(\varphi (n))$.
Our main result is as follows. Recall that in the lattice case, for the simplicity of presentation we agree to restrict ourselves to the case of a $(1;0)$-lattice.
Theorem 4. Assume that Conditions B1 and B2 are satisfied, $\min (m,n)\to \infty $ and $m=o(n)$. Then:
(1) if $\varphi (n)\to \infty $ as $n\to \infty $ so that $\varphi (n)=o(a_{m})$, then
Proof. For integers $0\leqslant r\leqslant n$, we consider the rescaled population size process $\mathcal{X}^{r,n}=\{X_{t}^{r,n},\,0\leqslant t\leqslant 1\} $ defined by
It was shown in [28], Theorem 1, that if $r_{1},r_{2},\dots$ is a sequence of natural numbers such that $r_{n}\leqslant n$ and $r_{n}\to \infty $ and $\varphi (n)\to \infty $ as $n\to \infty $ so that $\varphi (n)=o(a_{n})$, then
for any continuity point $x\in (0,\infty)$ of the distribution of the random variable $W$. Using (6.7), we conclude that for any $\varepsilon >0$ there exists $M=M(m,n)$ such that
and study separately the asymptotic behaviour of the terms on the right-hand side of this equality as $\min(m,n)\to \infty $, $m=o(n)$. In view of (6.9) it suffices to analyze the quantity
From this point until the end of the proof we assume that the distribution of the random variable $X_{1}$ is absolutely continuous. To prove the required assertions in the case when $X_{1}$ has a lattice distribution, it is necessary to replace $\displaystyle\int $ by $\sum $ throughout what follows.
$\mathbb{P}^{\uparrow}$-almost surely as $n-j\to \infty $. In addition, $H_{\infty}(k) >0$ $\mathbb{P}^{\uparrow}$-almost surely by Proposition 3.1 in [1]. Furthermore, for $q\in (-\sqrt{\varphi (n)},0]$ we have
Now we consider separately the cases when $\varphi (n)=o(a_{m})$, $\varphi (n)\sim Ta_{m}$ and $a_{m}=o(\varphi (n))$.
1) Assume that $\varphi (n)=o(a_{m})$. In this case, setting $H_{n}:=\mathbb{P}({Z_{n-j}\!>\!0\mid \mathcal{E}, Z_{0}\!=\!k})$ in (5.1) and recalling (6.12), for each $q\in [-\sqrt{\varphi (n)},0] $ we obtain
3) Finally, we consider the case when $m=m(n)\to \infty $ and $a_{m}=o(\varphi(n))$ as ${n\to \infty}$. We introduce the notation $S_{n-m,n}:=S_{n-m}-S_{n}$ and write the equality
for $z\in (-\infty,\infty)$. It is straightforward to verify that if we replace $S_{n-m}$ by $S_{n-m,n}$ and $S_{n-m-j}$ by $S_{n-m-j,n-j}$ in all relations between (6.10) and (6.11), then all estimates between (6.10) and (6.11) remain valid. Thus, we have
The authors are grateful to the referee, whose comments made it possible to improve the presentation of the results in the paper.
Bibliography
1.
V. I. Afanasyev, J. Geiger, G. Kersting and V. A. Vatutin, “Criticality for branching processes in random environment”, Ann. Probab., 33:2 (2005), 645–673
2.
V. I. Afanasyev, “Invariance principle for a critical Galton-Watson process attaining a high level”, Theory Probab. Appl., 55:4 (2011), 559–574
3.
V. I. Afanasyev, C. Böinghoff, G. Kersting and V. A. Vatutin, “Limit theorems for weakly subcritical branching processes in random environment”, J. Theoret. Probab., 25:3 (2012), 703–732
4.
J. Bertoin and R. A. Doney, “On conditioning a random walk to stay nonnegative”, Ann. Probab., 22:4 (1994), 2152–2167
5.
N. H. Bingham, C. M. Goldie and J. L. Teugels, Regular variation, Encyclopedia Math. Appl., 27, Cambridge Univ. Press, Cambridge, 1987, xx+491 pp.
6.
E. Bolthausen, “On a functional central limit theorem for random walks conditioned to stay positive”, Ann. Probab., 4:3 (1976), 480–485
7.
A. Bryn-Jones and R. A. Doney, “A functional limit theorem for random walk conditioned to stay non-negative”, J. London Math. Soc. (2), 74:1 (2006), 244–258
8.
L. Chaumont, “Excursion normalisée, méandre at pont pour les processus de Lévy stables”, Bull. Sci. Math., 121:5 (1997), 377–403
9.
F. Caravenna, “A local limit theorem for random walks conditioned to stay positive”, Probab. Theory Related Fields, 133:4 (2005), 508–530
10.
F. Caravenna and L. Chaumont, “Invariance principles for random walks conditioned to stay positive”, Ann. Inst. Henri Poincaré Probab. Stat., 44:1 (2008), 170–190
11.
F. Caravenna and L. Chaumont, “An invariance principle for random walk bridges conditioned to stay positive”, Electron. J. Probab., 18 (2013), 60, 32 pp.
12.
L. Chaumont and R. A. Doney, “Invariance principles for local times at the maximum of random walks and Lévy processes”, Ann. Probab., 38:4 (2010), 1368–1389
13.
F. den Hollander, Random polymers, École d'Été de Probabilités de Saint-Flour XXXVII – 2007, Lecture Notes in Math., 1974, Springer-Verlag, Berlin, 2009, xiv+258 pp.
14.
R. A. Doney, “Local behaviour of first passage probabilities”, Probab. Theory Related Fields, 152:3–4 (2012), 559–588
15.
W. Feller, An introduction to probability theory and its applications, v. 2, John Wiley & Sons, Inc., New York–London–Sydney, 1966, xviii+626 pp.
16.
B. V. Gnedenko and A. N. Kolmogorov, Limit distributions for sums of independent random variables, Addison-Wesley Publishing Company, Inc., Cambridge, MA, 1954, ix+264 pp.
17.
D. L. Iglehart, “Functional central limit theorems for random walks conditioned to stay positive”, Ann. Probab., 2:2 (1974), 608–619
18.
W. D. Kaigh, “An invariance principle for random walk conditioned by a late return to zero”, Ann. Probab., 4:1 (1976), 115–121
19.
G. Kersting and V. Vatutin, Discrete time branching processes in random environment, Math. Stat. Ser., John Wiley & Sons, London; ISTE, Hoboken, NJ, 2017, xiv+286 pp.
20.
T. M. Liggett, “An invariance principle for conditioned sums of independent random variables”, J. Math. Mech., 18:6 (1968), 559–570
21.
B. A. Rogozin, “The distribution of the first ladder moment and height and fluctuation of a random walk”, Theory Probab. Appl., 16:4 (1971), 575–595
22.
E. Seneta, Regularly varying functions, Lecture Notes in Math., 508, Springer-Verlag, Berlin–New York, 1976, v+112 pp.
23.
Ya. G. Sinai, “On the distribution of the first positive sum for a sequence of independent random variables”, Theory Probab. Appl., 2:1 (1957), 122–129
24.
C. Stone, “A local limit theorem for nonlattice multi-dimensional distribution functions”, Ann. Math. Statist., 36:2 (1965), 546–551
25.
V. A. Vatutin and V. Wachtel, “Local probabilities for random walks conditioned to stay positive”, Probab. Theory Related Fields, 143:1–2 (2009), 177–217
26.
V. Vatutin and E. Dyakonova, “Path to survival for the critical branching processes in a random environment”, J. Appl. Probab., 54:2 (2017), 588–602
27.
V. A. Vatutin and E. E. Dyakonova, Critical branching processes evolving in an unfavorable random environment, 2022, 15 pp., arXiv: 2209.13611
28.
V. A. Vatutin and E. E. Dyakonova, “Population size of a critical branching process evolving in an unfavorable environment”, Theory Probab. Appl., 68:3 (2023), 411–430
Citation:
V. A. Vatutin, C. Dong, E. E. Dyakonova, “Random walks conditioned to stay nonnegative and branching processes in an unfavourable environment”, Sb. Math., 214:11 (2023), 1501–1533
This publication is cited in the following 3 articles:
Vladimir A. Vatutin, Elena E. Dyakonova, “Branching processes under nonstandard conditions”, Stoch. Qual. Control, 39:1 (2024), 1–1
V. A. Vatutin, E. E. Dyakonova, “On the prospective minimum of the random walk conditioned to stay nonnegative”, Discrete Math. Appl., 34:6 (2024), 337–362
V. A. Vatutin, C. Dong, E. E. Dyakonova, “Some functionals for random walks and critical branching processes in an extremely unfavourable random environment”, Sb. Math., 215:10 (2024), 1321–1350