Fractal-Dimensional Properties of Subordinators

TL;DR

This paper investigates the fractal properties of subordinators, establishing a central limit theorem for the box-counting covering number, and introduces a new process to better understand the subordinator's dimension in relation to its Lévy measure.

Contribution

It provides a CLT for the box-counting number of subordinators and introduces a new process simplifying the analysis of their fractal dimensions.

Findings

Established a CLT for the number of boxes needed to cover a subordinator's range.

Proposed a new process by truncating jumps, facilitating easier analysis of fractal properties.

Connected the box-counting dimension to the Lévy measure of the subordinator.

Abstract

This work looks at the box-counting dimension of sets related to subordinators (non-decreasing L\'evy processes). It was recently shown in [Savov, 2014] that almost surely $\lim_{\delta\to0}U(\delta)N(t,\delta) = t$, where $N(t,\delta)$ is the minimal number of boxes of size at most $\delta$ needed to cover a subordinator's range up to time $t$, and $U(\delta)$ is the subordinator's renewal function. Our main result is a central limit theorem (CLT) for $N(t,\delta)$, complementing and refining work in [Savov, 2014]. Box-counting dimension is defined in terms of $N(t,\delta)$, but for subordinators we prove that it can also be defined using a new process obtained by shortening the original subordinator's jumps of size greater than $\delta$. This new process can be manipulated with remarkable ease in comparison to $N(t,\delta)$, and allows better understanding of the box-counting dimension of a subordinator's range in terms of its L\'evy measure, improving upon [Corollary 1, Savov, 2014]. Further, we shall prove corresponding CLT and almost sure convergence results for the new process.

Equations

$$\Phi (\lambda )=\text{d}\lambda +{\int }_0^{\infty }(1-e^{-\lambda x})\Pi (dx),$$

$$\underset{\delta \to 0}{\mathrm{liminf}}\frac{I(2\delta )}{I(\delta )}>1,$$

$$\frac{N(t,\delta )-ta(\delta )}{t^{\frac{1}{2}}b(\delta )}\stackrel{d}{\to }\mathcal{N}(0,1),$$

$$\underset{\delta \to 0}{\lim }\frac{L(t,\delta )}{\mu (\delta )}=t,$$

$$\frac{L(t,\delta )-t\mu (\delta )}{t^{\frac{1}{2}}v(\delta )\frac{}{}}\stackrel{d}{\to }\mathcal{N}(0,1)$$

$$\mathbb{E}[L(t,\delta )]=t\mu (\delta ),\text{Var}(L(t,\delta ))=tv(\delta ).$$

$$\underset{\delta \to 0}{\lim }\frac{U(\delta {)}^{\frac{7}{3}}}{{\sigma }_{\delta }^2}=0.$$

$$\left|\mathbb{P}\left(\frac{Y_1-\mathbb{E}[Y]+\cdots +Y_n-\mathbb{E}[Y]}{\text{Var}(Y{)}^{\frac{1}{2}}\sqrt{n}}\ge x\right)-\mathbb{P}(Z\ge x)\right|\le \frac{c\mathbb{E}[|Y-\mathbb{E}[Y]{|}^3]}{\text{Var}(Y{)}^{\frac{3}{2}}\sqrt{n}}.$$

$$\underset{\delta \to 0+}{\lim }\left|\mathbb{P}\left(\frac{N(1,\delta )-a(\delta )}{b(\delta )}\le x\right)-\mathbb{P}\left(Z\le x\right)\right|=0.$$

$$\mathbb{P}\left(\frac{N(1,\delta )-a(\delta )}{b(\delta )}\le x\right)=\mathbb{P}\left(N(1,\delta )\le a(\delta )+xb(\delta )\right),$$

$$\begin{split}\left(\ref{b}\right)=\mathbb{P}\big{(}N(1,\delta)\leq n\big{)}&=\mathbb{P}\left(T_{\delta}^{(n)}\geq 1\right)\hskip 5.69046pt=\hskip 5.69046pt\mathbb{P}\left(\hskip 6.48714pt\sum_{i=1}^{n}T_{\delta,i}\geq 1\hskip 6.48714pt\right)\\ &=\mathbb{P}\left(\hskip 11.38092pt\sum_{i=1}^{n}\big{(}T_{\delta,i}-U(\delta)\big{)}\geq 1-nU(\delta)\hskip 2.27626pt\right)\\ &=\mathbb{P}\left(\frac{\sum_{i=1}^{n}\big{(}T_{\delta,i}-U(\delta)\big{)}}{\sqrt{n\sigma_{\delta}^{2}}}\geq\frac{1-nU(\delta)}{\sqrt{n\sigma_{\delta}^{2}}}\right).\end{split}$$

$$n=\lceil a(\delta )+xb(\delta )\rceil \sim a(\delta )+xb(\delta )=a(\delta )+o(a(\delta ))\sim a(\delta )=U(\delta {)}^{-1}.$$

$$-x^{\prime }:$$

$$=\frac{1-1-xb(\delta )U(\delta )}{(a(\delta )+xb(\delta ){)}^{\frac{1}{2}}{\sigma }_{\delta }}\sim \frac{-xb(\delta )U(\delta )}{U(\delta {)}^{-\frac{1}{2}}{\sigma }_{\delta }}=\frac{-xb(\delta )U(\delta {)}^{\frac{3}{2}}}{{\sigma }_{\delta }}=-x.$$

$$\begin{array}{rl}\left|\mathbb{P}\left(\frac{N(1,\delta )-a(\delta )}{b(\delta )}\le x\right)-\mathbb{P}\left(Z\le x\right)\right|\le \left|\mathbb{P}\left(Z\ge -x^{\prime }\right)-\mathbb{P}\left(Z\ge -x\right)\right|& \\ +\left|\mathbb{P}\left(\frac{1}{\sqrt{n{\sigma }_{\delta }^2}}\sum _{i=1}^n\left(T_{\delta ,i}-U(\delta )\right)\ge -x^{\prime }\right)-\mathbb{P}\left(Z\ge -x^{\prime }\right)\right|& \\ =\left|\mathbb{P}\left(\frac{1}{\sqrt{n{\sigma }_{\delta }^2}}\sum _{i=1}^n\left(T_{\delta ,i}-U(\delta )\right)\ge -x^{\prime }\right)-\mathbb{P}\left(Z\ge -x^{\prime }\right)\right|& +o(1).\end{array}$$

$$\left(\text{ref}laststep\right)\le C\frac{\mathbb{E}[|T_{\delta }-U(\delta ){|}^3]}{{\sigma }_{\delta }^3{n}^{\frac{1}{2}}}+o(1)\sim C\frac{U(\delta {)}^{\frac{1}{2}}\mathbb{E}[|T_{\delta }-U(\delta ){|}^3]}{{\sigma }_{\delta }^3}.$$

$$\left(\text{ref}laststep\right)\le 8C\frac{U(\delta {)}^{\frac{1}{2}}U(\delta {)}^3}{{\sigma }_{\delta }^3}=8C{\left(\frac{U(\delta {)}^{\frac{7}{3}}}{{\sigma }_{\delta }^2}\right)}^{\frac{2}{3}}.$$

$$\underset{\delta \to 0+}{\mathrm{limsup}}\frac{\mathbb{E}[T_{\delta }^m]}{U(\delta {)}^m}<\infty .$$

$$\frac{\mathbb{E}[T_{\delta }^m]}{U(\delta {)}^m}=\mathbb{E}\left[{\left(\frac{T_{\delta }}{U(\delta )}\right)}^m\right]={\int }_0^{\infty }m{y}^{m-1}\mathbb{P}\left(\frac{T_{\delta }}{U(\delta )}>y\right)dy.$$

$$\frac{\mathbb{E}[T_{\delta }^m]}{U(\delta {)}^m}={\int }_0^{\infty }m{y}^{m-1}\mathbb{P}(X_{yU(\delta )}\le \delta )dy={\int }_0^{\infty }m{y}^{m-1}\mathbb{P}(e^{-\frac{1}{\delta }{X}_{yU(\delta )}}\ge e^{-1})dy.$$

$$\frac{\mathbb{E}[T_{\delta }^m]}{U(\delta {)}^m}\le {\int }_0^{\infty }m{y}^{m-1}{e}^{1-yU(\delta )\Phi (1/\delta )}dy\le {\int }_0^{\infty }m{y}^{m-1}{e}^{1-cy}dy,$$

$$\underset{\delta \to 0}{\mathrm{liminf}}\frac{I(\delta )}{I(2^{-n}\delta )}>\frac{1}{\eta }.$$

$$\underset{\delta \to 0}{\mathrm{liminf}}\frac{I(\delta )}{I(\delta /2)}=\underset{\delta \to 0}{\mathrm{liminf}}\frac{{\int }_0^{\delta }\overline{\Pi }(y)(dy)}{{\int }_0^{\delta /2}\overline{\Pi }(y)dy}=B.$$

$$\underset{\delta \to 0+}{\mathrm{liminf}}\frac{I(\delta )}{I(2^{-n}\delta )}=\underset{\delta \to 0+}{\mathrm{liminf}}\left(\frac{I(\delta )}{I(2^{-1}\delta )}\frac{I(2^{-1}\delta )}{I(2^{-2}\delta )}\cdots \frac{I(2^{-(n-1)}\delta )}{I(2^{-n}\delta )}\right)$$

$$\ge$$

$$\underset{\delta \to 0}{\mathrm{limsup}}\delta {\lambda }_{\delta }<\infty ,$$

$$\eta I(\delta )+I(2^{-n}\delta )$$

$$=e^{-{\lambda }_{\delta }{2}^{-n}\delta }\delta \overline{\Pi }(\delta )+e^{-{\lambda }_{\delta }{2}^{-n}\delta }{\int }_0^{\delta }x\Pi (dx)+{\int }_0^{2^{-n}\delta }x\Pi (dx).$$

$$\left(\text{ref}ttrr\right)\ge e^{-{\lambda }_{\delta }{2}^{-n}\delta }\delta \overline{\Pi }(\delta )+{\int }_{2^{-n}\delta }^{\delta }{e}^{-{\lambda }_{\delta }{2}^{-n}\delta }x\Pi (dx)+{\int }_0^{2^{-n}\delta }{e}^{-{\lambda }_{\delta }x}x\Pi (dx).$$

$$\left(\text{ref}ttrr\right)$$

$$=e^{-{\lambda }_{\delta }{2}^{-n}\delta }\delta \overline{\Pi }(\delta )+{\int }_0^{\delta }{e}^{-{\lambda }_{\delta }x}x\Pi (dx)\ge g({\lambda }_{\delta })=\frac{\delta }{(1+\alpha )U(\delta )}\ge \frac{I(\delta )}{(1+\alpha )K},$$

$${\sigma }_{\delta }^2=\text{Var}(T_{\delta })\ge \text{Var}(T_{\delta };|T_{\delta }-U(\delta )|\ge \alpha U(\delta ))$$

Full text

\newpagestyle

main\sethead[0][][0 \sectiontitle] \hdrtitle0 \sethead[0][ADAM BARKER][] FRACTAL-DIMENSIONAL PROPERTIES OF SUBORDINATORS0 \setfoot[][][]

FRACTAL-DIMENSIONAL PROPERTIES OF SUBORDINATORS

ADAM BARKER

Abstract

This work looks at the box-counting dimension of sets related to subordinators (non-decreasing Lévy processes). It was recently shown in [24] that almost surely $\lim_{\delta\rightarrow 0}U(\delta)N(t,\delta)=t$, where $N(t,\delta)$ is the minimal number of boxes of size at most $\delta$ needed to cover a subordinator’s range up to time $t$, and $U(\delta)$ is the subordinator’s renewal function. Our main result is a central limit theorem (CLT) for $N(t,\delta)$, complementing and refining work in [24].

Box-counting dimension is defined in terms of $N(t,\delta)$, but for subordinators we prove that it can also be defined using a new process obtained by shortening the original subordinator’s jumps of size greater than $\delta$. This new process can be manipulated with remarkable ease in comparison to $N(t,\delta)$, and allows better understanding of the box-counting dimension of a subordinator’s range in terms of its Lévy measure, improving upon [24, Corollary 1]. Further, we shall prove corresponding CLT and almost sure convergence results for the new process.

Introduction & Background

We shall mostly study the minimal number, $N(t,\delta)$, of intervals of length at most $\delta$ needed to cover the range $\{X_{s}:0\leq s\leq t\}$ of a subordinator $(X_{s})_{s\geq 0}$. The main result in this paper is a central limit theorem for $N(t,\delta)$, complementing the almost sure convergence result $\lim_{\delta\rightarrow 0}U(\delta)N(\delta,t)=t$, almost surely, where $U(\delta)$ denotes the renewal function of the subordinator, see [24, Theorem 1.1].

Prior to the results in [24], most works on box-counting dimension focused only on finding the value of $\lim_{\delta\rightarrow 0}\log(N(t,\delta))/\log(1/\delta)$, which defines the box-counting dimension. However, working with $N(t,\delta)$ itself allows precise understanding of its fluctuations around its mean, inaccessible at the log scale.

We will introduce an alternative “box-counting scheme” to $N(t,\delta)$, which allows us to understand the dimension of the range in terms of the Lévy measure, complementing results formulated in terms of the renewal function.

The fractal dimensional study of sets such as the range or graph of Lévy processes, and especially subordinators, has a very rich history. There are many works which study the box-counting, Hausdorff, and packing dimensions of sets related to Lévy processes [4, 6, 9, 11, 12, 13, 16, 17, 18, 19, 20, 24, 25, 26].

A Lévy process is a stochastic process in $\mathbb{R}^{d}$ which has stationary, independent increments, and starts at the origin. A subordinator $X:=(X_{t})_{t\geq 0}$ is a non-decreasing real-valued Lévy process. The Laplace exponent $\Phi$ of a subordinator $X$ is defined by the relation $e^{-\Phi(\lambda)}=\mathbb{E}[e^{-\lambda X_{1}}]$ for $\lambda\geq 0$. By the Lévy Khintchine formula [1, p72], $\Phi$ can always be expressed as

[TABLE]

where d is the linear drift, and $\Pi$ is the Lévy measure, which determines the size and intensity of the jumps (discontinuities) of $X$, and satisfies the condition $\int_{0}^{\infty}(1\wedge x)\Pi(dx)<\infty$. The renewal function is the expected first passage time above $\delta$, $U(\delta):=\mathbb{E}[T_{\delta}]$, where $T_{\delta}:=\int_{0}^{\infty}\mathbbm{1}_{\{X_{t}\leq\delta\}}dt$.

If the Lévy measure is infinite, then infinitesimally small jumps occur at an infinite rate, almost surely. We will not study processes with finite Lévy measure, as they have only finitely many jumps, and hence no fractal structure.

The box-counting dimension of a set in $\mathbb{R}^{d}$ is $\lim_{\delta\rightarrow 0}\log(N(\delta))/\log(1/\delta)$, where $N(\delta)$ is the minimal number of $d$-dimensional boxes of side length $\delta$ needed to cover the set. The limsup and liminf respectively define the upper and lower box-counting dimensions. For further background reading, we refer to [1, 2] for subordinators, [7, 21, 23] for Lévy processes, and[9, 26] for fractals.

The paper is structured as follows: Section 2 outlines the statements of all of the main results; Section 3 contains the proof of the CLT result for $N(t,\delta)$ and the lemmas required for this proof; Section 4 contains the proofs of all of the main results on the new process $L(t,\delta)$; Section 5 extends this work to the graph of a subordinator, and considers the special case of a subordinator with regularly varying Laplace exponent.

Main Results

A Central Limit Theorem for $N(t,\delta)$

Expanding upon Bertoin’s result [2, Theorem 5.1], the following almost sure limiting behaviour of $N(t,\delta)$ was determined by Savov [24, Theorem 1.1].

Theorem 2.1** (Savov, 2014).**

If a subordinator has infinite Lévy measure or a non-zero drift, then for all $t>0$, $\lim_{\delta\rightarrow 0+}U(\delta)N(t,\delta)=t$ almost surely.

We will complement and refine this work with a CLT on $N(t,\delta)$. When the subordinator has no drift, we require a mild condition on the Lévy measure:

[TABLE]

where $I(u):=\int_{0}^{u}\overline{\Pi}(x)dx$, and $\overline{\Pi}(x):=\Pi((x,\infty))$.

Remark 2.2**.**

Condition $\left(\ref{intcond}\right)$ has many equivalent formulations, see [1, Ex. III.7] and [3, Section 2.1]. We emphasise that $\left(\ref{intcond}\right)$ is far less restrictive than regular variation (or even $\mathcal{O}$-regular variation) of the Laplace exponent, and appears naturally in the context of the law of the iterated logarithm (see e.g. [1, p87]).

Theorem 2.3**.**

For every driftless subordinator with Lévy measure satisfying $\left(\ref{intcond}\right)$, for any $t>0$, $N(t,\delta)$ satisfies the following central limit theorem:

[TABLE]

as $\delta\rightarrow 0$, where $a(\delta):=U(\delta)^{-1}$, and $b(\delta):=U(\delta)^{-\frac{3}{2}}\textup{Var}(T_{\delta})^{\frac{1}{2}}$.

An Alternative Box-Counting Scheme, $L(t,\delta)$

Definition 2.4**.**

The process of $\delta$-shortened jumps, $\tilde{X}^{\delta}:=(\tilde{X}_{t}^{\delta})_{t\geq 0}$, is obtained by shortening all jumps of $X$ of size larger than $\delta$ to instead have size $\delta$. That is, $\tilde{X}^{\delta}$ is the subordinator with Laplace exponent $\tilde{\Phi}^{\delta}(u)=\textup{d}u+\int_{0}^{\delta}(1-e^{-ux})\tilde{\Pi}^{\delta}(dx)$ and Lévy measure $\tilde{\Pi}^{\delta}(dx)=\Pi(dx)\mathbbm{1}_{\{x<\delta\}}+\overline{\Pi}(\delta)\Delta_{\delta}$, where $\Delta_{\delta}$ denotes a unit point mass at $\delta$, and $\Pi$ is the Lévy measure of $X$.

Definition 2.5**.**

For $\delta,t>0$, $L(t,\delta)$ is defined by $L(t,\delta):=\frac{1}{\delta}\tilde{X}_{t}^{\delta}$.

We will see in Theorem 2.7 that $L(t,\delta)$ can replace $N(t,\delta)$ in the definition $\lim_{\delta\rightarrow 0}\log(N(t,\delta))/\log(1/\delta)$ of the box-counting dimension of the range of $X$. Then we will prove almost sure convergence and CLT results for $L(t,\delta)$.

Remark 2.6**.**

The log scale at which box-counting dimension is defined allows flexibility among functions to be taken in place of the optimal count. In particular, there is freedom between functions related by $f\asymp g$ asymptotically, where the notation means that there exist positive constants $A,B$ such that $Af(x)\leq g(x)\leq Bf(x)$ for all $x$. For more details, we refer to [9, p42].

Theorem 2.7**.**

For all $\delta,t>0$, for every subordinator, $N(t,\delta)\asymp L(t,\delta)$. In particular, by Remark 2.6, $L(t,\delta)$ can be used to define the box-counting dimension of the range, i.e. $\lim_{\delta\rightarrow 0}\log(N(t,\delta))/\log(1/\delta)=\lim_{\delta\rightarrow 0}\log(L(t,\delta))/\log(1/\delta)$.

Theorem 2.8**.**

For every subordinator with infinite Lévy measure, for all $t>0$,

[TABLE]

almost surely, where $\mu(\delta):=\frac{1}{\delta}(\textup{d}+I(\delta))$, and $I(\delta)=\int_{0}^{\delta}\overline{\Pi}(y)dy$.

Remark 2.9**.**

It can be deduced from [2, Prop 1.4] that $U(\delta)^{-1}\asymp\frac{1}{\delta}(I(\delta)+\ \text{{d}})$, for any subordinator. Theorems 2.1, 2.7 and 2.8 allow us to understand this relationship in terms of geometric properties of subordinators.

Theorem 2.10**.**

For every subordinator with infinite Lévy measure, for all $t>0$,

[TABLE]

as $\delta\rightarrow 0$, where $\mu(\delta)=\frac{1}{\delta}(\textup{d}+I(\delta))$, and $v(\delta):=\frac{1}{\delta}\left[\int_{0}^{\infty}(x\wedge\delta)^{2}\Pi(dx)\right]^{\frac{1}{2}}.$

Remark 2.11**.**

Applying Remark 2.4, the Lévy Khintchine formula $\left(\ref{lk}\right)$, and the fact that for any integrable function $f$, $\int_{0}^{\delta}f(x)\ \tilde{\Pi}^{\delta}(dx)=\int_{0}^{\infty}f(x\wedge\delta)\ \Pi(dx)$, it follows that for all $\delta,t>0$, the mean and variance of $L(t,\delta)$ are given by

[TABLE]

Computing the moments of $L(t,\delta)$ is remarkably simple in comparison to the moments of $N(t,\delta)$, which are not well known. This is a key benefit of using $L(t,\delta)$ to study the box-counting dimension of the range of a subordinator.

Proof of Theorem 2.3

A Sufficient Condition for Theorem 2.3

We will first work towards proving the following sufficient condition:

Lemma 3.1**.**

For every subordinator with infinite Lévy measure, a sufficient condition for the convergence in distribution $\left(\ref{thm}\right)$, with $\sigma_{\delta}^{2}:=\textup{Var}(T_{\delta})$, is

[TABLE]

The proof of Lemma 3.1 relies upon the Berry-Esseen Theorem, a very useful result for proving central limit theorem results as it provides the speed of convergence, which is stated here in Lemma 3.2. See [10, p542] for more details.

Lemma 3.2**.**

(Berry-Esseen Theorem) Let $Z\sim\mathcal{N}(0,1)$. There exists a finite constant $c>0$ such that for every collection of iid random variables $(Y_{k})_{k\in\mathbb{N}}$ with the same distribution as $Y$, where $Y$ has finite mean, finite absolute third moment, and finite non-zero variance, for all $n\in\mathbb{N}$ and $x\in\mathbb{R}$,

[TABLE]

For brevity, we will only provide calculations for $t=1$. The proofs for different values of $t$ are essentially the same. Recall the definitions $a(\delta):=U(\delta)^{-1}$, $\sigma_{\delta}^{2}:=\textup{Var}(T_{\delta})$, and $b(\delta):=U(\delta)^{-\frac{3}{2}}\sigma_{\delta}$. We shall aim to prove that for all $x\in\mathbb{R}$,

[TABLE]

For each $\delta>0$, $\left(\ref{berry}\right)$ provides an upper bound, and then under condition (2), we can prove that this bound converges to zero as $\delta\rightarrow 0$.

Proof of Lemma 3.1.

Let $T_{\delta}^{(k)}$ denote the $k$th time at which $N(t,\delta)$ increases, and let $T_{\delta,k}$, $k\in\mathbb{N}$, denote iid copies of $T_{\delta}^{(1)}$. By the strong Markov property, $T_{\delta}^{(k)}$ and $\sum_{i=1}^{k}T_{\delta,i}$ have the same distribution. Then, with $n:=\lceil a(\delta)+xb(\delta)\rceil$, where $\lceil\cdot\rceil$ denotes the ceiling function,

[TABLE]

and since $N(1,\delta)$ only takes integer values, using the fact that $T_{\delta}^{(n)}$ has the same distribution as the sum of $n$ iid copies of $T_{\delta}^{(1)}$, it follows that

[TABLE]

It follows from Lemma 3.3 that $\sigma_{\delta}^{2}\leq\mathbb{E}[T_{\delta}^{2}]\leq cU(\delta)^{2}$, which then implies that $b(\delta)=o(a(\delta))$ as $\delta\rightarrow 0$. Then, as $\delta\to 0$, the asymptotic behaviour of $n$ is

[TABLE]

It follows, with $x^{\prime}$ depending on $x$ and $\delta$, that as $\delta\rightarrow 0$,

[TABLE]

Now, by the triangle inequality and symmetry of the normal distribution, combining $(\ref{probexpre})$ and $(\ref{xx})$, it follows that as $\delta\to 0$, for any $x\in\mathbb{R}$,

[TABLE]

Recall that we wish to show that $\left(\ref{laststep}\right)$ converges to zero. By the Berry-Esseen Theorem and the fact that $n\sim U(\delta)^{-1}$, it follows that as $\delta\to 0$,

[TABLE]

Applying the triangle inequality, then Lemma 3.3 with $m=2$ and $m=3$ to $\mathbb{E}[|T_{\delta}-U(\delta)|^{3}]$, it follows that

[TABLE]

Therefore if the condition $\left(\ref{lemmm}\right)$ as in the statement of Lemma 3.1 holds, then the desired convergence in distribution $\left(\ref{thm}\right)$ follows, as required.

∎

Lemma 3.3**.**

For every subordinator with infinite Lévy measure, for all $m\geq 1$,

[TABLE]

Proof of Lemma 3.3.

First, by the moments and tails lemma (see [15, p26]),

[TABLE]

By the definition of $T_{\delta}$, it follows that $X_{u}\geq\delta$ if and only if $T_{\delta}\leq u$, and then

[TABLE]

Now, applying Markov’s inequality, the definition $\mathbb{E}[e^{-\lambda X_{t}}]=e^{-t\Phi(\lambda)}$, and the fact that $U(\delta)\Phi(1/\delta)\geq c$ for some constant $c$ (see [2, Prop 1.4]),

[TABLE]

which is finite and independent of $\delta$. Therefore the $\limsup$ is finite, as required. ∎

Proof of Theorem 2.3

Theorem 2.3 is proven by a contradiction, using Lemma 3.4 to show that the sufficient condition in Lemma 3.6 holds.

Lemma 3.4**.**

Recall the definition $I(\delta):=\int_{0}^{\delta}\overline{\Pi}(x)dx$. The condition $\left(\ref{intcond}\right)$ implies that for each $\eta\in(0,1)$, there exists a sufficiently large integer $n$ such that

[TABLE]

Proof of Lemma 3.4.

The integral condition $\left(\ref{intcond}\right)$ imposes that for some $B>1$,

[TABLE]

Then, by effectively replacing $1/2$ with $2^{-n}$ (so $1/2$ is replaced by a smaller constant), we can replace $B$ with $B^{n}$, which can be made arbitrarily large by choice of $n$. This follows by splitting up the fraction,

[TABLE]

where we simply take $n$ sufficiently large that $B^{n}>1/\eta$.

∎

Using Lemma 3.4 for a contradiction is the step in the proof of Theorem 2.3 which requires the condition $\left(\ref{intcond}\right)$. In order to prove Theorem 2.3, we require the notation introduced in Definition 3.5. We refer to [14, p93] for more details.

Definition 3.5**.**

Recalling from Remark 2.4 that the process $\tilde{X}^{\delta}$ has Laplace exponent $\tilde{\Phi}^{\delta}(u)=\textup{d}u+\int_{0}^{\delta}(1-e^{-ux})\Pi(dx)+(1-e^{-u\delta})\overline{\Pi}(\delta)$, we define:

(i) $g(u):=\frac{d}{du}\tilde{\Phi}^{\delta}(u)=\textup{d}+\int_{0}^{\delta}xe^{-ux}\tilde{\Pi}^{\delta}(dx)$,

(ii) $R(u):=\tilde{\Phi}^{\delta}(u)-ug(u)=\int_{0}^{\delta}\left(1-e^{-ux}(1+ux)\right)\tilde{\Pi}^{\delta}(dx)$,

(iii) $\lambda_{\delta}$ denotes the unique solution to $g(\lambda_{\delta})=x_{\delta}$, for $\textup{d}<x_{\delta}<\textup{d}+\int_{0}^{\delta}x\tilde{\Pi}^{\delta}(dx)$.

One can ignore the drift d in Definition 3.5, since $\textup{d}=0$ throughout Section 3. The proof of Theorem 2.3 now requires the following lemma:

Lemma 3.6**.**

For $\alpha>0$, $t=(1+\alpha)U(\delta)$, and $g(\lambda_{\delta})=x_{\delta}=\delta/t$, if

[TABLE]

then the desired convergence in distribution $\left(\ref{thm}\right)$, as in Theorem 2.3, holds.

Proof of Theorem 2.3.

Assume for a contradiction that there exists a sequence $(\delta_{m})_{m\geq 1}$ converging to zero, such that $\lim_{m\to\infty}\lambda_{\delta_{m}}\delta_{m}=\infty$. That is to say, assume that the sufficient condition in Lemma 3.6 doesn’t hold. For brevity, we omit the dependence of $\delta_{m}$ on $m$. Hence for all fixed $\eta,n>0$, $\eta\geq e^{-\lambda_{\delta}2^{-n}\delta}$ for all small enough $\delta>0$. By Fubini’s theorem, $I(\delta)=\int_{0}^{\delta}\overline{\Pi}(x)dx=\int_{0}^{\delta}x\tilde{\Pi}^{\delta}(dx)$, so

[TABLE]

Removing part of the first integral and noting $1\geq e^{-\lambda_{\delta}x}$ for all $x>0$,

[TABLE]

Now, $e^{-\lambda_{\delta}2^{-n}\delta}\geq e^{-\lambda_{\delta}x}$ for $x\geq 2^{-n}\delta$. So for $g(\lambda_{\delta})=x_{\delta}=\frac{\delta}{(1+\alpha)U(\delta)}$, where $\alpha>0$ is fixed and chosen sufficiently large that $x_{\delta}<\int_{0}^{\delta}x\tilde{\Pi}^{\delta}(dx)$ for all $\delta$ (this is possible by the relation $U(\delta)^{-1}\asymp I(\delta)/\delta$, see [2, Prop 1.4]),

[TABLE]

where the last two inequalities respectively follow by Definition 2.4, Definition 3.5 $(i)$ with $\textup{d}=0$, and the relation $U(\delta)^{-1}\asymp I(\delta)/\delta$, see [1, p74]. So for a constant $K>0$, for all sufficiently small $\delta>0$, we have shown $\eta I(\delta)+I(2^{-n}\delta)\geq\frac{I(\delta)}{(1+\alpha)K}$.

Taking $\eta>0$ small enough that $\frac{1}{(1+\alpha)K}\geq 2\eta$, it follows that $I(2^{-n}\delta)\geq\eta I(\delta)$, and hence $I(\delta)/I(2^{-n}\delta)\leq 1/\eta$. But in Lemma 3.4 we showed that for each fixed $\eta>0$, there is sufficiently large $n$ such that $\liminf_{\delta\to 0}I(\delta)/I(2^{-n}\delta)>1/\eta$, which is a contradiction, so the sufficient condition as in Lemma 3.6 must hold.

∎

Remark 3.7**.**

For a driftless subordinator, Theorem 2.3 holds under the same condition $\left(\ref{intcond}\right)$ applied to the function $H(y):=\int_{0}^{y}x\Pi(dx)$ rather than the integrated tail function $I$. The integrated tail $I(y)=H(y)+y\overline{\Pi}(y)$ depends on the large jumps of $X$ since $\overline{\Pi}(x)=\Pi((x,\infty))$, but $H$ does not depend on the large jumps, so these conditions are substantially different.

With only minor changes, the argument as in the proof of Theorem 2.3 works with $H$ in place of $I$. Under condition $\left(\ref{intcond}\right)$ for $H$ in place of $I$, one can prove that Lemma 3.4 holds with $H$ in place of $I$. Then we assume for a contradiction that there exists a sequence $(\delta_{m})_{m\geq 1}$ converging to zero, such that $\lim_{m\to\infty}\lambda_{\delta_{m}}\delta_{m}=\infty$. But then as in the proof of Theorem 2.3, one can deduce that $\eta H(\delta)+H(2^{-n}\delta)\geq\frac{1}{(1+\alpha)K^{\prime}}H(\delta)$, which contradicts the analogous Lemma 3.4 result with $H$ in place of $I$. **

Remark 3.8**.**

Theorem 2.3 can also be proven for subordinators with a drift $\textup{d}>0$, under a stronger regularity condition. For $Y_{t}:=X_{t}-\textup{d}t$, define $\Phi_{Y}$ as the Laplace exponent of $Y$. The convergence in distribution $\left(\ref{thm}\right)$ holds whenever $\limsup_{x\to 0}x^{-5/6}\Phi_{Y}(x)<\infty$. This is proven using Remark 3.10, the inequality $\mathbb{P}(Y_{t}<a)\geq 1-Cth(a)$ for all Lévy processes (see [22, p954] for details), and the asymptotic expansion of $U(\delta)$ as in [8, Theorem 4].**

Proofs of Lemmas 3.9, 3.12, 3.6

Lemmas 3.9, 3.12, and 3.6 give sufficient conditions for Theorem 2.3 to hold. The proofs for these lemmas are facilitated by Lemma 3.11, which was proven in 1987 by Jain and Pruitt [14, p94]. Recall that $\tilde{X}^{\delta}$ denotes the process with $\delta$-shortened jumps, as defined in Definition 2.4.

Lemma 3.9**.**

The convergence in distribution $\left(\ref{thm}\right)$ as in Theorem 2.3 holds if for some $\alpha\in(0,1]$, $\liminf_{\delta\rightarrow 0}\left[\mathbb{P}\left(\tilde{X}^{\delta}_{(1+\alpha)U(\delta)}\leq\delta\right)+\mathbb{P}\left(\tilde{X}^{\delta}_{(1-\alpha)U(\delta)}\geq\delta\right)\right]>0$.

Proof of Lemma 3.9.

For all $\alpha>0$, recalling that $\mathbb{E}[T_{\delta}]=U(\delta)$,

[TABLE]

[TABLE]

For the desired convergence in distribution $\left(\ref{thm}\right)$ to hold, it is sufficient by Lemma 3.1 to show that $\lim_{\delta\rightarrow 0}U(\delta)^{\frac{7}{3}}/\sigma_{\delta}^{2}=0$. Now,

[TABLE]

Note that $T_{\delta}\geq t$ if and only if $\tilde{X}^{\delta}_{t}\leq\delta$ since jumps of size larger than $\delta$ do not occur in either case, and so $X_{t}=\tilde{X}^{\delta}_{t}$ when $T_{\delta}\geq t$. It follows that $\left(\ref{thm}\right)$ holds if

[TABLE]

∎

Remark 3.10**.**

The condition in Lemma 3.9 is not optimal. If for $\varepsilon\in\left(0,\frac{1}{6}\right)$, $\lim_{\delta\rightarrow 0}U(\delta)^{2\varepsilon-\frac{1}{3}}\left[\mathbb{P}\left(\tilde{X}^{\delta}_{U(\delta)+U(\delta)^{1+\varepsilon}}\leq\delta\right)+\mathbb{P}\left(\tilde{X}^{\delta}_{U(\delta)-U(\delta)^{1+\varepsilon}}\geq\delta\right)\right]=\infty,$ then the convergence in distribution $\left(\ref{thm}\right)$ follows too. This stronger condition does not lead to any more generality than the condition $\left(\ref{intcond}\right)$ for driftless subordinators.

Lemma 3.11** (Jain, Pruitt [14, Lemma 5.2]).**

There exists $c>0$ such that for every $\varepsilon>0$, $t\geq 0$ and $x_{\delta}>0$ satisfying $\textup{d}=g(\infty)<x_{\delta}<g(0)=\textup{d}+\int_{0}^{\delta}x\tilde{\Pi}^{\delta}(dx)$,

[TABLE]

Lemma 3.12**.**

For $\alpha>0$, $t=(1+\alpha)U(\delta)$, and $g(\lambda_{\delta})=x_{\delta}=\delta/t$, if

[TABLE]

then the desired convergence in distribution $\left(\ref{thm}\right)$, as in Theorem 2.3, holds.

Proof of Lemma 3.12.

Applying the inequality (16) from Lemma 3.11,

[TABLE]

Now, letting $\limsup_{\delta\rightarrow 0}tR(\lambda_{\delta})<\infty$, we will consider two separate cases:

(i) If $\liminf_{\delta\rightarrow 0}tR(\lambda_{\delta})=\beta>0$, then by choice of $\varepsilon>0$ such that $\frac{1+\varepsilon}{\varepsilon^{2}}=\frac{\beta}{2c}$, the lower bound in (17) is larger than a positive constant as $\delta\rightarrow 0$.

(ii) If $\liminf_{\delta\rightarrow 0}tR(\lambda_{\delta})=0$, then imposing $\varepsilon=2c/(tR(\lambda_{\delta}))$, the lower bound in (17) is again larger than a positive constant as $\delta\rightarrow 0$. The desired convergence in distribution $\left(\ref{thm}\right)$ then follows in each case by Lemma 3.9.

∎

Proof of Lemma 3.6.

Noting that $1-e^{-y}(1+y)\leq y$ for all $y>0$,

[TABLE]

[TABLE]

Then by the relation $U(\delta)I(\delta)\leq C\delta$ for a constant $C$ (see [2, Prop 1.4]),

[TABLE]

So we can conclude that if $\limsup_{\delta\rightarrow 0}\delta\lambda_{\delta}<\infty$, then the desired convergence in distribution $\left(\ref{thm}\right)$ follows by Lemma 3.12.

∎

Proofs of Results on $L(t,\delta)$

Firstly, we prove Theorem 2.7, which confirms that $L(t,\delta)$ can replace $N(t,\delta)$ in the definition of the box-counting dimension of the range. This is done by showing that $L(t,\delta)\asymp N(t,\delta)$, which is known to be sufficient by Remark 2.6.

Proof of Theorem 2.7.

The jumps of the original subordinator $X$ and the process with shortened jumps $\tilde{X}^{\delta}$ are all the same size, other than jumps bigger than size $\delta$. The optimal number of intervals to cover the range, $N(X,t,\delta)$, always increases by $1$ at each jump bigger than size $\delta$, regardless of its size, so it follows that $N(X,t,\delta)=N(\tilde{X}^{\delta},t,\delta)$, with the obvious notation.

Instead of counting the number $N(X,t,\delta)$ of boxes needed to cover the range of $X$, consider those needed for the range of the subordinator $X^{(0,\delta)}$ with Lévy measure $\Pi(dx)\mathbbm{1}_{\{x<\delta\}}$ (so all jumps of size larger than $\delta$ are removed), and adding $Y_{t}^{\delta}$, which counts the number of jumps larger than size $\delta$ of $X$. It follows that $N(X,t,\delta)\leq N(X^{(0,\delta)},t,\delta)+Y_{t}^{\delta}\leq 2N(X,t,\delta)$.

Consider $M(X^{(0,\delta)},t,\delta)$, the number of intervals in a lattice of side length $\delta$ to intersect with the range of $X^{(0,\delta)}$. It is easy to show that $N(t,\delta)\asymp M(t,\delta)$ (see [9, p42]). Also, $M(X^{(0,\delta)},t,\delta)=\lceil\frac{1}{\delta}X_{t}^{(0,\delta)}\rceil$, since $X^{(0,\delta)}$ has no jumps of size larger than $\delta$. Now, $\frac{1}{\delta}X_{t}^{(0,\delta)}\asymp\lceil\frac{1}{\delta}X_{t}^{(0,\delta)}\rceil$ for small enough $\delta$, and hence

[TABLE]

[TABLE]

By Remark 2.6, $\lim_{\delta\rightarrow 0}\frac{\log(L(t,\delta))}{\log(1/\delta)}=\lim_{\delta\rightarrow 0}\frac{\log(N(t,\delta))}{\log(1/\delta)}$, and hence $L(t,\delta)$ can be used to define the box-counting dimension of the range of any subordinator.

∎

Next we will prove the CLT result for $L(t,\delta)$, working with $t=1$ for brevity. The proof is essentially the same for other values of $t>0$. We will show convergence of the Laplace transform of $\frac{1}{v(\delta)}(L(1,\delta)-\mu(\delta))$ to that of the standard normal distribution. Recall that $Z\sim\mathcal{N}(0,1)$ has Laplace transform $\mathbb{E}[e^{-\lambda Z}]=e^{\lambda^{2}/2}$.

Proof of Theorem 2.10.

By Remark 2.4 and $\left(\ref{lk}\right)$, $\delta L(t,\delta)=\tilde{X}_{t}^{\delta}$ is a subordinator with Laplace exponent $\tilde{\Phi}^{\delta}$, and it follows that for any $\lambda\geq 0$,

[TABLE]

Recalling the definition $\mu(\delta)=\frac{1}{\delta}(\textup{d}+I(\delta))$, where $I(\delta):=\int_{0}^{\delta}x\tilde{\Pi}^{\delta}(dx)$, and writing $\tilde{\Phi}^{\delta}$ in the Lévy Khintchine representation as in $\left(\ref{lk}\right)$, it follows that

[TABLE]

[TABLE]

Then applying the fact that $\frac{y^{2}}{2}-\frac{y^{3}}{6}\leq y-1+e^{-y}\leq\frac{y^{2}}{2}$ for all $y>0$,

[TABLE]

By the definition of $v(\delta)$, it follows that $v(\delta)^{2}=\frac{1}{\delta^{2}}\int_{0}^{\delta}x^{2}\tilde{\Pi}^{\delta}(dx)$, and so

[TABLE]

It is then sufficient, in order to show that $\left(\ref{terms}\right)$ converges to $\frac{\lambda^{2}}{2}$, to prove that

[TABLE]

Again by the definition of $v(\delta)$, for $\left(\ref{lcltcond}\right)$ to hold we require both

[TABLE]

Squaring the expression in $\left(\ref{int1}\right)$, since $x\leq\delta$ within each integral, it follows that

[TABLE]

By the binomial expansion, $(a+b)^{3}\geq 3a^{2}b$ for $a,b>0$, and then

[TABLE]

since the Lévy measure is infinite. For $\left(\ref{int2}\right)$, simply observe that

[TABLE]

∎

Next we will prove the almost sure convergence result for $L(t,\delta)$. If there is a drift and the Lévy measure is finite, then the result is trivial. So we need only consider cases with infinite Lévy measure, and begin with the zero drift case. Using a Borel-Cantelli argument (see [15, p32] for details), we shall prove that $\liminf_{\delta\to 0}L(t,\delta)/\mu(\delta)=\limsup_{\delta\to 0}L(t,\delta)/\mu(\delta)=t$ almost surely.

First, we will prove the almost sure convergence to $t$ along a subsequence $\delta_{n}$ converging to zero. Then, by monotonicity of $\mu(\delta)$ and $L(t,\delta)$, we will deduce that for all $\delta$ between $\delta_{n}$ and $\delta_{n+1}$, $L(t,\delta)/\mu(\delta)$ also tends to $t$ as $\delta_{n}\to 0$.

Proof of Theorem 2.8.

For all $\varepsilon>0$, by Chebyshev’s inequality and Remark 2.11,

[TABLE]

[TABLE]

Recall that $\mu(\delta)=\int_{0}^{\infty}\frac{1}{\delta}(x\wedge\delta)\ \Pi(dx)$, so since $\frac{1}{\delta}(x\wedge\delta)$ is non-decreasing as $\delta$ decreases, it follows that $\mu(\delta)$ is non-decreasing as $\delta$ decreases. Now, $\lim_{\delta\to 0}\mu(\delta)=\infty$, and $\mu$ is continuous, so it follows that for any fixed $r\in(0,1)$ there is a decreasing sequence $\delta_{n}$ such that $\mu(\delta_{n})=r^{-n}$ for each $n$. Then $\left(\ref{summ}\right)$ is finite, so by the Borel-Cantelli lemma, $\lim_{n\to\infty}L(t,\delta_{n}/\mu(\delta_{n})=t$ almost surely.

When there is no drift, $L(t,\delta)$ is given by changing the original subordinator’s jump sizes from $y$ to $\frac{1}{\delta}(y\wedge\delta)$. By monotonicity of this map, it follows that for a fixed sample path of the original subordinator, each individual jump of the process $L(t,\delta_{n+1})$ is at least as big as the corresponding jump of the process $L(t,\delta_{n})$. So $L(t,\delta)$ is non-decreasing as $\delta$ decreases, and so for all $\delta_{n+1}\leq\delta\leq\delta_{n}$,

[TABLE]

Then by our choice of the subsequence $\delta_{n}$, it follows that for all $\delta_{n+1}\leq\delta\leq\delta_{n}$,

[TABLE]

and since $\lim_{n\to\infty}L(t,\delta_{n})/\mu(\delta_{n})=t$, it follows that

[TABLE]

Taking limits as $r\to 1$, it follows that $\lim_{\delta\to 0}L(t,\delta)/\mu(\delta)=t$ almost surely.

For a process with a positive drift $\textup{d}>0$ and infinite Lévy measure, denote the scaling term obtained by removing the drift as $\hat{\mu}(\delta):=\mu(\delta)-\textup{d}/\delta$. Then the above Borel-Cantelli argument for $\hat{\mu}$ yields the almost sure limit along a subsequence $\hat{\delta}_{n}$ as in $\left(\ref{summ}\right)$. Then since the functions $\mu(\delta)$ and $L(t,\delta)$ are again monotone in $\delta$ when there is a drift, the argument applies as in $\left(\ref{lll}\right)$.

∎

Remark 4.1**.**

Theorem 2.8 is formulated in terms of the characteristics of the subordinator (i.e. the drift and Lévy measure). For $N(t,\delta)$, the almost sure behaviour in Theorem 2.1 is formulated in terms of the renewal function, and in order to write this in terms of the characteristics, the expression is more complicated than for $L(t,\delta)$. For details, see [24, Corollary 1] and [8, Prop 1], the latter of which is very powerful for understanding the asymptotics of $U(\delta)$ for subordinators with a positive drift, significantly improving upon results in [5]. **

Extensions and Special Cases

Extensions: Box-Counting Dimension of the Graph

The graph of a subordinator $X$ up to time $t$ is the set $\{(s,X_{s}):0\leq s\leq t\}$. The box-counting dimensions of the range and graph are closely related. This is evident when we consider the mesh box counting schemes $M_{G}(t,\delta)$, $M_{R}(t,\delta)$, denoting graph and range respectively. The mesh box-counting scheme counts the number of boxes in a lattice of side length $\delta$ to intersect with a set.

Remark 5.1**.**

For every subordinator with infinite Lévy measure or a positive drift, $M_{G}(t,\delta)=\left\lfloor t/\delta\right\rfloor+M_{R}(t,\delta),$ where $\lfloor\cdot\rfloor$ denotes the floor function. Indeed, $M_{R}(t,\delta)$ increases by 1 if and only if $M_{G}(t,\delta)$ increases by 1 and the new box for the graph lies directly above the previous box. For each integer $n$, $M_{G}(t,\delta)$ also increases at time $n\delta$, the new box directly to the right of the previous box.

Remark 5.2**.**

It follows that the graph of every subordinator $X$ has the same box-counting dimension as the range of $X^{\prime}_{t}:=t+X_{t}$, the original process plus a unit drift.

Proposition 5.3**.**

For every subordinator with drift d $>0$, the box-counting dimensions of the range and graph agree almost surely.

Proof of Proposition 5.3.

Letting $T_{(\delta,\infty)}$ denote the first passage time of the subordinator above $\delta$, consider an optimal covering of the graph with squares of side length $\delta$ as follows:

Starting with $[0,\delta]\times[0,\delta]$, at time $T_{1}:=\min(T_{(\delta,\infty)},\delta)$, add a new box $[T_{1},T_{1}+\delta]\times[X_{T_{1}},X_{T_{1}}+\delta]$, and so on. Denote the number of these boxes by $N_{G}(t,\delta)$, and write $N_{R}(t,\delta)$ as the optimal number of boxes needed to cover the range.

If d $\geq 1$, then we have $T_{1}=T_{(\delta,\infty)}$ because $X_{\delta}\geq\text{d}\delta$. It follows that each time $N_{G}(t,\delta)$ increases by $1$, so does $N_{R}(t,\delta)$, and vice versa, so $N_{G}(t,\delta)=N_{R}(t,\delta)$, and the box-counting dimension of the range and graph are equal when $\textup{d}\geq 1$.

For d $\in(0,1)$, a similar argument applies with a covering of $\frac{\delta}{\textup{d}}\times\delta$ rectangles rather than $\delta\times\delta$ squares. Starting with $[0,\frac{\delta}{\textup{d}}]\times[0,\delta]$, at time $T_{1}$, add a new box $[T_{1},T_{1}+\frac{\delta}{\textup{d}}]\times[X_{T_{1}},X_{T_{1}}+\delta]$, and so on. The number of these boxes is again $N_{R}(\delta,t)$, since $X_{\frac{\delta}{\textup{d}}}\geq\delta$. By Remark 2.6 , this covering of rectangles can still be used to define the box-counting dimension of the range, since for $k:=\left\lceil\frac{1}{\textup{d}}\right\rceil$, with $N_{G}(t,\delta)$ and $N_{G}^{\prime}(t,\delta)$ as the number of squares and of rectangles respectively,

[TABLE]

∎

Remark 5.4**.**

The box-counting dimension of the graph of every subordinator is 1 almost surely, since subordinators have bounded variation (BV) almost surely. The same is true for the graph of all BV functions/processes, including in particular every Lévy process without a Gaussian component, whose Lévy measure satisfies $\int(1\wedge|x|)\Pi(dx)<\infty$. By Proposition 5.3, the box-counting dimension of the range of every subordinator with drift $\textup{d}>0$ is 1 almost surely. **

Special Cases: Regular Variation of the Laplace Exponent

Corollary 5.5 is analogous to [24, Corollary 2], with $L(t,\delta)$ in place of $N(t,\delta)$. This allows very fine comparisons, not visible at the log-scale, to be made between subordinators whose Laplace exponents are regularly varying with the same index.

Corollary 5.5**.**

Consider a subordinator whose Laplace exponent is regularly varying at infinity, such that $\Phi(\lambda)\sim\lambda^{\alpha}F(\lambda)$ for $\alpha\in(0,1)$, where $F(\cdot)$ is a slowly varying function. Then almost surely as $\delta\to 0$, for all $t>0$,

[TABLE]

Proof of Corollary 5.5.

Note that $\textup{d}=0$, i.e. there is no drift, when the Laplace exponent is regularly varying of index $\alpha\in(0,1)$. By Theorem 2.8, as $\delta\to 0$,

[TABLE]

Since $\Phi$ is regularly varying at [math], as $x\to 0$, $\overline{\Pi}(x)\sim\Phi(\frac{1}{x})/\Gamma(1-\alpha)$ (see [1, p75]). Then by Karamata’s Theorem (see [3, Prop. 1.5.8]), almost surely as $\delta\to 0$,

[TABLE]

∎

Corollary 5.6 strengthens the result of Theorem 2.7 when the Laplace exponent $\Phi$ is regularly varying. The result can not be strengthened in general, as the relationship between $\mu(\delta)$ and $U(\delta)^{-1}$ is “$\asymp$” rather than “$\sim$” (see[2, Prop. 1.4]).

Corollary 5.6**.**

For a subordinator with Laplace exponent $\Phi$ regularly varying at infinity with index $\alpha\in(0,1)$, for all $t>0$, almost surely as $\delta\to 0$,

[TABLE]

Corollary 5.6 follows immediately from Corollary 5.5 and [24, Corollary 2], which says that when the Laplace exponent $\Phi$ is regularly varying at infinity, such that $\Phi(\lambda)\sim\lambda^{\alpha}F(\lambda)$ for $\alpha\in(0,1)$, where $F(\cdot)$ is a slowly varying function, for all $t>0$, almost surely as $\delta\to 0$,

[TABLE]

Remark 5.7**.**

For $\alpha\in(0,1)$, $\Gamma(2-\alpha)\Gamma(1+\alpha)$ takes values between $\pi/4$ and $1$. So $L(t,\delta)$ and $N(t,\delta)$ are closely related when the Laplace exponent is regularly varying, but as $\delta\to 0$, $L(t,\delta)$ grows to infinity slightly faster than $N(t,\delta)$.

Acknowledgements

Many thanks to Mladen Savov for guiding the author towards this interesting topic, and for numerous helpful discussions related to this work. Further thanks to Ron Doney for his feedback on an early draft of this paper, and thanks to an anonymous referee for their valuable comments on this work.