We consider a single server queue that has a threshold to change its arrival process and service speed by its queue length, which is referred to as a two-level GI/G/1 queue. This model is motivated by an energy saving problem for a single server queue whose arrival process and service speed are controlled. To obtain its performance in tractable form, we study the limit of the stationary distribution of the queue length in this two-level queue under scaling in heavy traffic. Except for a special case, this limit corresponds to its diffusion approximation. It is shown that this limiting distribution is truncated exponential (or uniform if the drift is null) below the threshold level and exponential above it under suitably chosen system parameters and generally distributed interarrival times and workloads brought by customers. This result is proved under a mild limitation on arrival parameters using the so-called basic adjoint relationship (BAR) approach studied in Braverman, Dai, and Miyazawa (2017, 2024) and Miyazawa (2017, 2024). We also intuitively discuss about a diffusion process corresponding to the limit of the stationary distribution under scaling.

The payoff in the Chow–Robbins coin-tossing game is the proportion of heads when you stop. Stopping to maximize expectation was addressed by Chow and Robbins (1965), who proved there exist integers ${k_n}$ such that it is optimal to stop at n tosses when heads minus tails is ${k_n}$. Finding ${k_n}$ was unsolved except for finitely many cases by computer. We prove an $o(n^{-1/4})$ estimate of the stopping boundary of Dvoretsky (1967), which then proves ${k_n} = \left\lceil {\alpha \sqrt n \,\, - 1/2\,\, + \,\,\frac{{\left( { - 2\zeta (\! -1/2)} \right)\sqrt \alpha }}{{\sqrt \pi }}{n^{ - 1/4}}} \right\rceil $ except for n in a set of density asymptotic to 0, at a power law rate. Here, $\alpha$ is the Shepp–Walker constant from the Brownian motion analog, and $\zeta$ is Riemann’s zeta function. An $n^{ - 1/4}$ dependence was conjectured by Christensen and Fischer (2022). Our proof uses moments involving Catalan and Shapiro Catalan triangle numbers which appear in a tree resulting from backward induction, and a generalized backward induction principle. It was motivated by an idea of Häggström and Wästlund (2013) to use backward induction of upper and lower Value bounds from a horizon, which they used numerically to settle a few cases. Christensen and Fischer, with much better bounds, settled many more cases. We use Skorohod’s embedding to get simple upper and lower bounds from the Brownian analog; our upper bound is the one found by Christensen and Fischer in another way. We use them first for yet many more examples and a conjecture, then algebraically in the tree, with feedback to get much sharper Value bounds near the border, and analytic results. Also, we give a formula that gives the exact optimal stop rule for all n up to about a third of a billion; it uses the analytic result plus terms arrived at empirically.

Let X be the sum of a diffusion process and a Lévy jump process, and for any integer $n\ge 1$ let $\phi_n$ be a function defined on $\mathbb{R}^2$ and taking values in $\mathbb{R}$, with adequate properties. We study the convergence of functionals of the type

\begin{align*} \Delta_n\sum_{i=1}^{[t/\Delta_n]} \phi_n\!\left(\alpha_n\left[\frac{X_{(i-1)\Delta_n}}{\alpha_n}\right], \:\frac{\alpha_n}{\sqrt{\Delta_n}}\left(\left[\frac{X_{i\Delta_n}}{\alpha_n}\right]-\left[\frac{X_{(i-1)\Delta_n}}{\alpha_n}\right]\right)\right),\end{align*}

$(\Delta_n)$

$(\alpha_n)$

$n\to +\infty$

$\frac{\alpha_n}{\sqrt{\Delta_n}}$

$[0,+\infty)$]

This paper concerns an insurance firm’s surplus process observed at renewal inspection times, with a focus on assessing the probability of the surplus level dropping below zero. For various types of inter-inspection time distributions, an explicit expression for the corresponding transform is given. In addition, Cramér–Lundberg-type asymptotics are established. Also, an importance sampling-based Monte Carlo algorithm is proposed, and is shown to be logarithmically efficient.

The Wright–Fisher model, originating in Wright (1931) is one of the canonical probabilistic models used in mathematical population genetics to study how genetic type frequencies evolve in time. In this paper we bound the rate of convergence of the stationary distribution for a finite population Wright–Fisher Markov chain with parent-independent mutation to the Dirichlet distribution. Our result improves the rate of convergence established in Gan et al. (2017) from $\mathrm{O}(1/\sqrt{N})$ to $\mathrm{O}(1/N)$. The results are derived using Stein’s method, in particular, the prelimit generator comparison method.

We study an optimal inventory control problem under a reflected jump–diffusion netflow process with state-dependent jumps, in which the intensity of the jump process can depend on the inventory level. We examine the well-posedness of the associated integro-differential Hamilton–Jacobi–Bellman (ID-HJB) equation with Neumann boundary condition in the classical sense. To achieve this, we first establish the existence of viscosity solutions to the ID-HJB equation of an auxiliary control problem with a compact policy space, which is proved to be equivalent to the primal problem. We reformulate the ID-HJB equation as a Neumann HJB equation with the (non-local) integral term expressed in terms of the value function of the auxiliary problem and prove the existence of a unique classical solution to the Neumann HJB equation. Then, the well-posedness of the primal ID-HJB equation follows from the unique classical solution of the Neumann HJB equation and the existence of viscosity solutions to the auxiliary ID-HJB equation. Based on this classical solution, we characterize the optimal (admissible) inventory control strategy and show the verification result for the primal control problem.

Motivated by the investigation of probability distributions with finite variance but heavy tails, we study infinitely divisible laws whose Lévy measure is characterized by a radial component of geometric (tempered) stable type. We closely investigate the univariate case: characteristic exponents and cumulants are calculated, as well as spectral densities; absolute continuity relations are shown, and short- and long-time scaling limits of the associated Lévy processes analyzed. Finally, we derive some properties of the involved probability density functions.

We analyse a Markovian SIR epidemic model where individuals either recover naturally or are diagnosed, leading to isolation and potential contact tracing. Our focus is on digital contact tracing via a tracing app, considering both its standalone use and its combination with manual tracing. We prove that as the population size n grows large, the epidemic process converges to a limiting process, which, unlike with typical epidemic models, is not a branching process due to dependencies created by contact tracing. However, by grouping to-be-traced individuals into macro-individuals, we derive a multi-type branching process interpretation, allowing computation of the reproduction number R. This is then converted to an individual reproduction number $R^\mathrm{(ind)}$, which, in contrast to R, decays monotonically with the fraction of app-users, while both share the same threshold at 1. Finally, we compare digital (only) contact tracing and manual (only) contact tracing, proving that the critical fraction of app-users, $\pi_{\mathrm{c}}$, required for $R=1$ is higher than the critical fraction manually contact-traced, $p_{\mathrm{c}}$, for manual tracing.

We study a variant of the classical Markovian logistic SIS epidemic model on a complete graph, which has the additional feature that healthy individuals can become infected without contacting an infected member of the population. This additional ‘self-infection’ is used to model situations where there is an unknown source of infection or an external disease reservoir, such as an animal carrier population. In contrast to the classical logistic SIS epidemic model, the version with self-infection has a non-degenerate stationary distribution, and we derive precise asymptotics for the time to converge to stationarity (mixing time) as the population size becomes large. It turns out that the chain exhibits the cutoff phenomenon, which is a sharp transition in time from one to zero of the total variation distance to stationarity. We obtain the exact leading constant for the cutoff time and show that the window size is of constant (optimal) order. While this result is interesting in its own right, an additional contribution of this work is that the proof illustrates a recently formalised methodology of Barbour, Brightwell and Luczak (2022), ‘Long-term concentration of measure and cut-off’, Stochastic Processes and their Applications 152, 378–423, which can be used to show cutoff via a combination of concentration-of-measure inequalities for the trajectory of the chain and coupling techniques.

In this paper, we study asymptotic behaviors of a subcritical branching Brownian motion with drift $-\rho$, killed upon exiting $(0, \infty)$, and offspring distribution $\{p_k{:}\; k\ge 0\}$. Let $\widetilde{\zeta}^{-\rho}$ be the extinction time of this subcritical branching killed Brownian motion, $\widetilde{M}_t^{-\rho}$ the maximal position of all the particles alive at time t and $\widetilde{M}^{-\rho}:\!=\max_{t\ge 0}\widetilde{M}_t^{-\rho}$ the all-time maximal position. Let $\mathbb{P}_x$ be the law of this subcritical branching killed Brownian motion when the initial particle is located at $x\in (0,\infty)$. Under the assumption $\sum_{k=1}^\infty k ({\log}\; k) p_k <\infty$, we establish the decay rates of $\mathbb{P}_x(\widetilde{\zeta}^{-\rho}>t)$ and $\mathbb{P}_x(\widetilde{M}^{-\rho}>y)$ as t and y respectively tend to $\infty$. We also establish the decay rate of $\mathbb{P}_x(\widetilde{M}_t^{-\rho}> z(t,\rho))$ as $t\to\infty$, where $z(t,\rho)=\sqrt{t}z-\rho t$ for $\rho\leq 0$ and $z(t,\rho)=z$ for $\rho>0$. As a consequence, we obtain a Yaglom-type limit theorem.

Competing and complementary risk (CCR) problems are often modelled using a class of distributions of the maximum, or minimum, of a random number of independent and identically distributed random variables, called the CCR class of distributions. While CCR distributions generally do not have an easy-to-calculate density or probability mass function, two special cases, namely the Poisson–exponential and exponential–geometric distributions, can easily be calculated. Hence, it is of interest to approximate CCR distributions with these simpler distributions. In this paper, we develop Stein’s method for the CCR class of distributions to provide a general comparison method for bounding the distance between two CCR distributions, and we contrast this approach with bounds obtained using a Lindeberg argument. We detail the comparisons for Poisson–exponential, and exponential–geometric distributions.

In this paper, we study the asymptotic behavior of the generalized Zagreb indices of the classical Erdős–Rényi (ER) random graph G(n, p), as $n\to\infty$. For any integer $k\ge1$, we first give an expression for the kth-order generalized Zagreb index in terms of the number of star graphs of various sizes in any simple graph. The explicit formulas for the first two moments of the generalized Zagreb indices of an ER random graph are then obtained from this expression. Based on the asymptotic normality of the numbers of star graphs of various sizes, several joint limit laws are established for a finite number of generalized Zagreb indices with a phase transition for p in different regimes. Finally, we provide a necessary and sufficient condition for any single generalized Zagreb index of G(n, p) to be asymptotic normal.