Abstract:
A generalization of the well-known result concerning the survival probability of a critical branching process in random environment $Z_k$ is considered. The triangular array scheme of branching processes in random environment $Z_{k,n}$ that are close to $Z_k$ for large $n$ is studied. The equivalence of the survival probabilities for the processes $Z_{n,n}$ and $Z_n$ is obtained under rather natural assumptions on the closeness of $Z_{k,n}$ and $Z_k$.
Bibliography: 7 titles.
Keywords:random walks, branching processes, random environments.
Branching processes in random environment (BPREs) are a generalization of well-known Galton–Watson branching processes. Unlike the latter, the distribution of the number of offspring in each generation in a BPRE depends on some random factor, called an environment. Like in the case of Galton–Watson branching processes, the subcritical, critical and supercritical cases can be considered for BPREs. There are results concerning the survival probability for all three types of BPREs in the case when the environment consists of independent and identically distributed random elements.
In this paper we consider a critical BPRE. The asymptotic behaviour of the survival probability of the process in the case when the generating function of the distribution for the number of offspring of a single particle is linear fractional was deduced in [1]. With no assumption on the explicit form of the distribution for the number of offspring of a single particle, some results were obtained in [2] and generalized in [3].
A sequence $\Xi=\{\xi_i,\,i \in \mathbb{N}\}$ of independent and identically distributed random elements with values in a measurable space $(Y, \mathcal{G})$ is called a random environment. We denote by $(\Omega, \mathcal{F}, \mathsf{P}^{\Xi})$ the probability space on which a random environment $\Xi$ is defined.
We consider a family of generating functions
$$
\begin{equation*}
\{f_y,\, y \in Y\}.
\end{equation*}
\notag
$$
The process $\mathcal{Z}=\{Z_k,\,k \geqslant 0\}$ is called a branching process in a random environment $\Xi$.
For each fixed $\omega \in \Omega$ we introduce a probability measure $\mathsf{P}_{\omega}$ on the space of sequences $\mathbb{N}_0^{\infty}$, $\mathbb{N}_0 :=\{0\} \cup \mathbb{N}$. For any $G \subset \mathbb{N}_0^{\infty}$ and $V \in \mathcal{F}$ we set
and extend the probability measure $\mathsf{P}$ to $(\mathbb{N}_0^{\infty} \times \Omega, \sigma(2^{\mathbb{N}_0^{\infty}} \times \mathcal{F}))$. Note that
By the definition of $\Xi$ the random variables $X_i$, $i \geqslant 1$, are independent and identically distributed. The sequence of random variables $\{S_k,\,k \in \mathbb{N}\}$ is called the associated random walk.
Assumption 1. $\mathsf{P}$-almost surely, $F_0'(1), F_0''(1) \in (0, \infty)$ and
where $c_{-}=\sum_{k=1}^{\infty} k^{-1} (\mathsf{P}(S_k<0)-1/2)$ and $\Upsilon$ is some positive constant.
Remark 1. Theorem 1 follows from the results obtained in [3] and [4], but this fact requires a separate substantiation, which is presented in § 3.
In this paper we consider the triangular array scheme of branching processes $\{Z_{k,n},\, {k \!\leqslant\! n}\}$ in the random environment $\Xi$. Our goal is to find conditions on $Z_{k,n}$ under which the probabilities $\mathsf{P}(Z_n>0)$ and $\mathsf{P}(Z_{n,n}>0)$ are equivalent as $n \to \infty$.
The paper is structured as follows. In § 2 the main result is presented. In § 3 a theorem on survival probability for a BPRE is proved. In § 4 several assertions used to substantiate the main result are established.
The set of random variables $\{Z_{k,n},\,0 \leqslant k \leqslant n,\,Z_{0,n}=1\}$ is said to be a perturbed branching process in the random environment $\Xi$ (PBPRE).
We state the following assumption on the smallness of perturbations.
Assumption 3. For all $y \in Y$, $0 \leqslant i<n$ and $s \in [0, 1]$,
Theorem 3 (see [4], Theorem 5.1). Relation (1.1) holds under Assumptions 1 and 6.
In the proof of Theorem 5.1 in [4], Assumption C is used only to prove Lemma 5.5. We prove this lemma in the case when Assumptions 1 and 2 hold.
The statement of Lemma 5.5 involves a certain measure $\mathsf{P}^+$. We will not dwell on the details of its construction: it was thoroughly described in [4], § 5.
Proof. In was deduced in the proof of Lemma 5.5 in [4] that by Assumption 1, for any $\delta \in (0, 1/2)$ there exists a set $\Omega''=\Omega''(\delta)$, $\mathsf{P}^+(\Omega'')=1$, such that for any $\omega \in \Omega''$ the inequality
for $x=x_j=\exp\{j^{1/2-\delta_1}\}$. The second term on the right-hand side of (3.10) is finite by Assumption 2. We prove the finiteness of the first term on the right-hand side of (3.10).
If the first factor on the right-hand side of (3.11) is finite, then the first term on the right-hand side of (3.10) is too. Assume the contrary. Then there exists an increasing sequence $\{n_k \in \mathbb{N} \mid k \in \mathbb{N}\}$ such that $h_{n_k}^0>1/n_k$. The sequence $h_j^0$ does not increase, which yields
Since $n_k \geqslant k$ for all $k \in \mathbb{N}$ the relation (3.12) contradicts Assumption 2 in view of Cauchy’s criterion.
By virtue of the convergence of the two series on the right-hand side of (3.10) and the Borel–Cantelli lemma, there exists a set $\Omega'''$, $\mathsf{P}^+(\Omega''')=1$, such that for any $\omega \in \Omega'''$ there exists a positive function $D_2(\omega)$ such that
This relation replaces Lemma 5.8 in [4], which is used there in the proof of Theorem 5.1. The other assertions in [4] involved in the proof of Theorem 5.1, namely Lemmas 5.2 and 5.9 and Theorem 4.6, use only Assumption 1.
for arbitrary natural $k$ and $l$. Using (4.2) and (4.3) we arrive at the required equality (4.1).
The lemma is proved.
To prove Theorem 2 we divide the double sum in the representation from Lemma 3 into the two parts corresponding to the subsets of indices $k \leqslant m$ and $k>m$, respectively. We show that the second part is negligible in comparison with the first for large $m$. For this purpose we need the following estimate.
Lemma 4. Let $\widehat{L}_k$ denote $\min\{S_i + \theta(i) \mid 0 \leqslant i \leqslant k\}$. Then, under Assumption 1, for any $\delta \in (0, 1/2)$ there exists a constant $K$ such that
We estimate the quantity $p_2$. We define recursively a random variable $\nu_i$. Set ${\nu_0=0}$ and assume that $\nu_i$ is defined for $i \geqslant 0$; then
We estimate the first term on the right-hand side of (4.8). We set $\Delta_0 \,{:=}\, 0$ and $\Delta_i :=\nu_i-\nu_{i-1}$ for $i>0$ and note that the $\Delta_i$ are independent, identically distributed random variables for $i>0$. Since the events $\{\rho_k \leqslant r\}$ and $\{\nu_{r+1}>k\}$ coincide for an arbitrary $r$, Markov’s inequality and the subadditivity of the function $g(x) :=x^{1/2-\delta/2}$, $x \geqslant 0$, imply the relation
By virtue of (4.13), Markov’s inequality, and since the $\eta_i$ are nonnegative, the second term on the right-hand side of (4.8) is estimated as follows:
Let $\omega \in Q_{k,n}$. On the strength of the fact that the event $\{Z_{k,n}>0\}$ is embedded in $\{Z_{\widehat{\tau}_k,n}>0\}$ and $Z_{\widehat{\tau}_k,n}$ is integer valued, we have
Owing to Lemma 6, we have reduced the investigation of the asymptotic behaviour of the sequence $\sum_{k=0}^n \sum_{l=1}^{\infty} A_{k,n,l} B_{n-k,l,n}$ to examining the sums $\sum_{k=0}^m \sum_{l=1}^{\infty} A_{k,n,l} B_{n-k,l,n}$. Now we constrain the values of $l$.
Lemma 7. Under Assumptions 1 and 3, for fixed $m$ there exists a sequence of positive numbers $\{\beta_M= \beta_M(m),\,M \in \mathbb{N}\}$ such that $\beta_M \to 0$ as $M \to \infty$ and
Proof. We prove that for each $k \in \mathbb{N}$ there exists a sequence of positive numbers $\{\beta_M^{(k)},\,M \in \mathbb{N}\}$ tending to zero as $M \to \infty$ such that
Assume that the left-hand side of (4.35) does not tend to zero for some $k$. Then there are a positive number $\varepsilon$ and an increasing sequence of natural numbers $\{M_r,\,r \in \mathbb{N}\}$ such that
If $\sup\{n_r \mid r \in \mathbb{N}\}<\infty$, then there exists a natural number that is infinitely often repeated in the sequence $\{n_r,\,r \in \mathbb{N}\}$. We denote this number by $n$. By the definition of $n$ there is an increasing sequence of natural numbers $\widetilde{M}_r$ such that
However, relation (4.37) contradicts the continuity of the probability measure.
If $\sup\{n_r \mid r \in \mathbb{N}\}=\infty$, then there exist increasing sequences of natural numbers $\{\widetilde{n}_r,\,r \in \mathbb{N}\}$ and $\{\widetilde{M}_r,\,r \in \mathbb{N}\}$ such that
for all $r>R$. By virtue of Assumption 3 and Remark 2, the left-hand side of (4.40) tends to $\mathsf{P}(Z_k>M)$ as $r \to \infty$, which yields the estimate $\mathsf{P}(Z_k>M) \geqslant \varepsilon$. However, this contradicts (4.39).
In all above cases we arrive at a contradiction; hence (4.35) is true. By the inequality $B_{n-k,l,n} \leqslant \mathsf{P}(L_{k,n} \geqslant 0)$, Lemma 1 and (4.35) we have
Using Lemmas 6 and 7 we can reduce the problem to the examination of a finite number of combinations $(k, l)$. In this case we can use Assumption 3 to switch from $F_{i-1,n}$ to $F_{i-1}$ for $i \leqslant k$.
Lemma 8. Under Assumptions 1 and 3, for fixed $m$ and for any $\varepsilon>0$ there exist $M=M(\varepsilon)$ and $N=N(\varepsilon)$ such that
It follows from Assumption 3 and Remark 2 that for $k \leqslant m$ $Z_{k,n}$ converges to $Z_k$ in distribution as $n \to \infty$ for all $(\xi_1, \dots, \xi_k)$. By Lebesgue’s dominated convergence theorem we have
It follows from Lemmas 6–8 that the main contribution to the asymptotic behaviour of the survival probability of a PBPRE is made by the sum $\sum_{k=0}^m \sum_{l=1}^M A_{k,l} B_{n-k,l,n}$. It remains to establish a result on the asymptotic behaviour of $B_{n-k,l,n}$. Since the number of terms in the sum under consideration is finite, there is no question of uniform asymptotics with respect to $k$ and $l$.
In what follows we need an expression for survival probability in the case when the environment is fixed.
Lemma 9 (see [4], Propositions 1.3 and 1.4). The survival probability of the PBPRE $Z_{k,n}$ satisfies
To prove the next assertion we must pass to the limit under the condition that some part of the trajectory of the sequence $S_i$, $i \geqslant 0$, is nonnegative. A tool making it possible to take the limit is the measure $\mathsf{P}^+$. The properties of this measure are described in § 3.
Lemma 10. Under Assumptions 1 and 5, for any $k \in \mathbb{N}_0$ there is a set $\Omega'$, $\mathsf{P}^+(\Omega')=1$, such that for all $\omega\in \Omega'$ and $\varepsilon>0$ there exists a parameter ${R=R(\omega, \varepsilon)}$ such that
Proof. Our argument is like in the proof of Lemma 2.
By virtue of Assumption 1, relation (3.7) for $\delta :=\min\{\delta_2, \delta_3\}/2$ holds for some $\Omega''=\Omega''(\delta)$ such that $\mathsf{P}^+(\Omega'')=1$.
To estimate the sum (4.49) we need an estimate for the probability
for $x=x_j=\exp\{j^{1/2-\delta_3}\}$. The convergence of the series (4.52) and the Borel–Cantelli lemma imply the existence of a set $\Omega'''$, $\mathsf{P}^+(\Omega''')=1$, such that for any $\omega \in \Omega'''$ there is a positive function $D_3(\omega)$ such that
Since $\delta<\delta_2$ and $\delta<\delta_3$, the terms of the series on the right-hand side of (4.54) are exponentially small. It follows that for all $\omega \in \Omega'$ and $\varepsilon>0$ there exists $R=R(\omega, \varepsilon)$ such that the left-hand side of inequality (4.49) is at most $\varepsilon$ for $r>R$ and $n>k+r$.
By Lemma 1 and Assumption 4 the second term on the right-hand side of (4.58) admits the estimate
$$
\begin{equation}
l \, \mathsf{P}(\overline{Q}_n \mid L_{k,n} \geqslant 0) \leqslant l \frac{\mathsf{P}(\overline{Q}_n)}{\mathsf{P}(L_{k,n} \geqslant 0)}\to 0, \qquad n \to \infty.
\end{equation}
\tag{4.59}
$$
We estimate the first term on the right in (4.58). Fix $\varepsilon>0$. By virtue of Lemmas 9 and 10 there exists a set $\Omega' \in \Omega, \mathsf{P}^+(\Omega')=1$ such that for any $\omega \in \Omega'$ there exists $R=R(\omega, \varepsilon)$ such that
Since $\varepsilon$ is arbitrary, we conclude that the sequence of random variables ${\mathsf{I}_{Q_n} |\widetilde{\pi}_{k,n}-\widetilde{\pi}_{k,n}^{\,0}|}$ converges to $0$ $\mathsf{P}^+$-almost surely as $n \to \infty$. As this sequence is uniformly bounded, it follows from [4], Lemma 5.2, that
Proof. Fix $\omega \in \Omega$ and $j \in \mathbb{N}_0$. By Assumption 3 and Theorem 3.1.1 in [7] there exist a probability space $(\widehat{\Omega}_{\omega}, \widehat{\mathcal{F}}_{\omega}, \widehat{\mathsf{P}}_{\omega})$ and random variables $\widehat{Y}_{\omega, n}$, $j<n$, and $\widehat{Y}_{\omega}$ defined on it such that $F_{j;\omega}$ is the generating function for $\widehat{Y}_{\omega}$, $F_{j,n;\omega}$ is the generating function for $\widehat{Y}_{\omega,n}$, and $\widehat{Y}_{\omega,n} \to \widehat{Y}_{\omega}$ as $n \to \infty$ $\widehat{\mathsf{P}}_{\omega}$-almost surely. We denote the mean on this space by $\widehat{\mathsf{E}}_{\omega}$. By Fatou’s lemma,
for arbitrary $\omega \in \Omega$ and $n \in \mathbb{N}$. Note that, by virtue of Assumption 4, the sequence $\mathsf{I}_{Q_n}$ converges to $1$ in probability as $n \to \infty$. By Riesz’s theorem, there is a subsequence $\mathsf{I}_{Q_{n_k}}$ that converges to $1$ $\mathsf{P}$-almost surely as $k \to \infty$. It follows that
Lemma 12 implies that Assumption 2 is fulfilled. By Assumptions 1 and 2, Theorem 1 is true for the BPRE $\{Z_n,\,n \geqslant 0\}$, which implies the convergence
The author is grateful to A. V. Shklyaev for his permanent support of this research. The author is grateful to anonymous referees for their comments, which made it possible to improve the presentation considerably.
Bibliography
1.
M. V. Kozlov, “On the asymptotic behavior of the probability of non-extinction for critical branching processes in a random environment”, Theory Probab. Appl., 21:4 (1977), 791–804
2.
J. Geiger and G. Kersting, Theory Probab. Appl., 45:3 (2001), 517–525
3.
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
4.
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.
5.
D. Denisov, A. Sakhanenko and V. Wachtel, “First-passage times for random walks with nonidentically distributed increments”, Ann. Probab., 46:6 (2018), 3313–3350
6.
V. V. Petrov, Sums of independent random variables, Akademie-Verlag, Berlin, 1975, x+348 pp.
7.
A. V. Skorokhod, “Limit theorems for stochastic processes”, Theory Probab. Appl., 1:3 (1956), 261–290
Citation:
V. V. Kharlamov, “Asymptotic behaviour of the survival probability of almost critical branching processes in a random environment”, Mat. Sb., 215:1 (2024), 131–152; Sb. Math., 215:1 (2024), 119–140