In this note, we formulate a ‘one-sided’ version of Wormald’s differential equation method. In the standard ‘two-sided’ method, one is given a family of random variables that evolve over time and which satisfy some conditions, including a tight estimate of the expected change in each variable over one-time step. These estimates for the expected one-step changes suggest that the variables ought to be close to the solution of a certain system of differential equations, and the standard method concludes that this is indeed the case. We give a result for the case where instead of a tight estimate for each variable’s expected one-step change, we have only an upper bound. Our proof is very simple and is flexible enough that if we instead assume tight estimates on the variables, then we recover the conclusion of the standard differential equation method.

We introduce and study a game-theoretic model to understand the spread of an epidemic in a homogeneous population. A discrete-time stochastic process is considered where, in each epoch, first, a randomly chosen agent updates their action trying to maximize a proposed utility function, and then agents who have viral exposures beyond their immunity get infected. Our main results discuss asymptotic limiting distributions of both the cardinality of the subset of infected agents and the action profile, considered under various values of two parameters (initial action and immunity profile). We also show that the theoretical distributions are almost always achieved in the first few epochs.

We prove a large deviation principle for the slow-fast rough differential equations (RDEs) under the controlled rough path (RP) framework. The driver RPs are lifted from the mixed fractional Brownian motion (FBM) with Hurst parameter $H\in (1/3,1/2)$. Our approach is based on the continuity of the solution mapping and the variational framework for mixed FBM. By utilizing the variational representation, our problem is transformed into a qualitative property of the controlled system. In particular, the fast RDE coincides with Itô stochastic differential equation (SDE) almost surely, which possesses a unique invariant probability measure with frozen slow component. We then demonstrate the weak convergence of the controlled slow component by averaging with respect to the invariant measure of the fast equation and exploiting the continuity of the solution mapping.

Let $\{\omega _n\}_{n\geq 1}$ be a sequence of independent and identically distributed random variables on a probability space $(\Omega , \mathcal {F}, \mathbb {P})$, each uniformly distributed on the unit circle $\mathbb {T}$, and let $\ell _n=cn^{-\tau }$ for some $c>0$ and $0<\tau <1$. Let $I_{n}=(\omega _n,\omega _n+\ell _n)$ be the random interval with left endpoint $\omega _n$ and length $\ell _n$. We study the asymptotic property of the covering time $N_n(x)=\sharp \{1\leq k\leq n: x\in I_k\}$ for each $x\in \mathbb {T}$. We prove the quenched central limit theorem for the covering time, that is, $\mathbb {P}$-almost surely,

$$ \begin{align*}\frac{N_n(x)-\mathbb{E}_{\mathbb{P}}(N_n(x))}{\sqrt{\sum_{k=1}^n \ell_k(1-\ell_k)}}\end{align*} $$

converges in law to the standard normal distribution.

We study the last exit time that a spectrally negative Lévy process is below zero until it reaches a positive level b, denoted by $g_{\tau_b^+}$. We generalize the results of the infinite-horizon last exit time explored by Chiu and Yin (2005) by incorporating a random horizon $\tau_b^+$, which represents the first passage time above b. We derive an explicit expression for the joint Laplace transform of $g_{\tau_b^+}$ and $\tau_b^+$ by utilizing a hybrid observation scheme approach proposed by Li, Willmot, and Wong (2018). We further study the optimal prediction of $g_{\tau_b^+}$ in the $L_1$ sense, and find that the optimal stopping time is the first passage time above a level $y_b^{\ast}$, with an explicit characterization of the stopping boundary $y_b^{\ast}$. As examples, Brownian motion with drift and the Cramér–Lundberg model with exponential jumps are considered.

We consider stochastic differential equations (SDEs) driven by a fractional Brownian motion with a drift coefficient that is allowed to be arbitrarily close to criticality in a scaling sense. We develop a comprehensive solution theory that includes strong existence, path-by-path uniqueness, existence of a solution flow of diffeomorphisms, Malliavin differentiability and $\rho $-irregularity. As a consequence, we can also treat McKean-Vlasov, transport and continuity equations.

Let $G$ be a group. The notion of linear sofic approximations of $G$ over an arbitrary field $F$ was introduced and systematically studied by Arzhantseva and Păunescu [3]. Inspired by one of the results of [3], we introduce and study the invariant $\kappa _F(G)$ that captures the quality of linear sofic approximations of $G$ over $F$. In this work, we show that when $F$ has characteristic zero and $G$ is linear sofic over $F$, then $\kappa _F(G)$ takes values in the interval $[1/2,1]$ and $1/2$ cannot be replaced by any larger value. Further, we show that under the same conditions, $\kappa _F(G)=1$ when $G$ is torsion-free. These results answer a question posed by Arzhantseva and Păunescu [3] for fields of characteristic zero. One of the new ingredients of our proofs is an effective non-concentration estimates for random walks on finitely generated abelian groups, which may be of independent interest.

In this paper, the model of bisexual branching processes affected by viral infectivity and with random control functions in independent and identically distributed (i.i.d.) random environments is established and the Markov property is given firstly. Then the relations of the probability generating functions of this model are studied, and some sufficient conditions for process extinction under common mating functions are presented. Finally, the limiting behaviors of the considered model after proper normalization, such as the sufficient conditions for the convergence in L1 and L2 and almost everywhere convergence, are investigated under the condition that the random control functions are super additive.

We consider the extremes of the logarithm of the characteristic polynomial of matrices from the C$\beta $E ensemble. We prove convergence in distribution of the centered maxima (of the real and imaginary parts) toward the sum of a Gumbel variable and another independent variable, which we characterize as the total mass of a ‘derivative martingale’. We also provide a description of the landscape near extrema points.

We study a version of the stochastic control problem of minimizing the sum of running and controlling costs, where control opportunities are restricted to independent Poisson arrival times. Under a general setting driven by a general Lévy process, we show the optimality of a periodic barrier strategy, which moves the process upward to the barrier whenever it is observed to be below it. The convergence of the optimal solutions to those in the continuous-observation case is also shown.

Let $B^{H}$ be a d-dimensional fractional Brownian motion with Hurst index $H\in(0,1)$, $f\,:\,[0,1]\longrightarrow\mathbb{R}^{d}$ a Borel function, and $E\subset[0,1]$, $F\subset\mathbb{R}^{d}$ are given Borel sets. The focus of this paper is on hitting probabilities of the non-centered Gaussian process $B^{H}+f$. It aims to highlight how each component f, E and F is involved in determining the upper and lower bounds of $\mathbb{P}\{(B^H+f)(E)\cap F\neq \emptyset \}$. When F is a singleton and f is a general measurable drift, some new estimates are obtained for the last probability by means of suitable Hausdorff measure and capacity of the graph $Gr_E(f)$. As application we deal with the issue of polarity of points for $(B^H+f)\vert_E$ (the restriction of $B^H+f$ to the subset $E\subset (0,\infty)$).

Early investigation of Pólya urns considered drawing balls one at a time. In the last two decades, several authors have considered multiple drawing in each step, but mostly for schemes involving two colors. In this manuscript, we consider multiple drawing from urns of balls of multiple colors, formulating asymptotic theory for specific urn classes and addressing more applications. The class we consider is affine and tenable, built around a ‘core’ square matrix. We examine cases where the urn is irreducible and demonstrate its relationship to matrix irreducibility for its core matrix, with examples provided. An index for the drawing schema is derived from the eigenvalues of the core. We identify three regimes: small, critical, and large index. In the small-index regime, we find an asymptotic Gaussian law. In the critical-index regime, we also find an asymptotic Gaussian law, albeit with a difference in the scale factor, which involves logarithmic terms. In both of these regimes, we have explicit forms for the structure of the mean and the covariance matrix of the composition vector (both exact and asymptotic). In all three regimes we have strong laws.

We provide explicit small-time formulae for the at-the-money implied volatility, skew, and curvature in a large class of models, including rough volatility models and their multi-factor versions. Our general setup encompasses both European options on a stock and VIX options, thereby providing new insights on their joint calibration. The tools used are essentially based on Malliavin calculus for Gaussian processes. We develop a detailed theoretical and numerical analysis of the two-factor rough Bergomi model and provide insights on the interplay between the different parameters for joint SPX–VIX smile calibration.

We solve the non-discounted, finite-horizon optimal stopping problem of a Gauss–Markov bridge by using a time-space transformation approach. The associated optimal stopping boundary is proved to be Lipschitz continuous on any closed interval that excludes the horizon, and it is characterized by the unique solution of an integral equation. A Picard iteration algorithm is discussed and implemented to exemplify the numerical computation and geometry of the optimal stopping boundary for some illustrative cases.

We consider Markov processes that alternate continuous motions and jumps in a general locally compact Polish space. Starting from a mechanistic construction, a first contribution of this article is to provide conditions on the dynamics so that the associated transition kernel forms a Feller semigroup, and to deduce the corresponding infinitesimal generator. As a second contribution, we investigate the ergodic properties in the special case where the jumps consist of births and deaths, a situation observed in several applications including epidemiology, ecology, and microbiology. Based on a coupling argument, we obtain conditions for convergence to a stationary measure with a geometric rate of convergence. Throughout the article, we illustrate our results using general examples of systems of interacting particles in $\mathbb{R}^d$ with births and deaths. We show that in some cases the stationary measure can be made explicit and corresponds to a Gibbs measure on a compact subset of $\mathbb{R}^d$. Our examples include in particular Gibbs measures associated to repulsive Lennard-Jones potentials and to Riesz potentials.

Persistent Betti numbers are a major tool in persistent homology, a subfield of topological data analysis. Many tools in persistent homology rely on the properties of persistent Betti numbers considered as a two-dimensional stochastic process $ (r,s) \mapsto n^{-1/2} (\beta^{r,s}_q ( \mathcal{K}(n^{1/d} \mathcal{X}_n))-\mathbb{E}[\beta^{r,s}_q ( \mathcal{K}( n^{1/d} \mathcal{X}_n))])$. So far, pointwise limit theorems have been established in various settings. In particular, the pointwise asymptotic normality of (persistent) Betti numbers has been established for stationary Poisson processes and binomial processes with constant intensity function in the so-called critical (or thermodynamic) regime; see Yogeshwaran et al. (Prob. Theory Relat. Fields 167, 2017) and Hiraoka et al. (Ann. Appl. Prob. 28, 2018).

In this contribution, we derive a strong stabilization property (in the spirit of Penrose and Yukich, Ann. Appl. Prob. 11, 2001) of persistent Betti numbers, and we generalize the existing results on their asymptotic normality to the multivariate case and to a broader class of underlying Poisson and binomial processes. Most importantly, we show that multivariate asymptotic normality holds for all pairs (r, s), $0\le r\le s<\infty$, and that it is not affected by percolation effects in the underlying random geometric graph.

We present a closed-form solution to a discounted optimal stopping zero-sum game in a model based on a generalised geometric Brownian motion with coefficients depending on its running maximum and minimum processes. The optimal stopping times forming a Nash equilibrium are shown to be the first times at which the original process hits certain boundaries depending on the running values of the associated maximum and minimum processes. The proof is based on the reduction of the original game to the equivalent coupled free-boundary problem and the solution of the latter problem by means of the smooth-fit and normal-reflection conditions. We show that the optimal stopping boundaries are partially determined as either unique solutions to the appropriate system of arithmetic equations or unique solutions to the appropriate first-order nonlinear ordinary differential equations. The results obtained are related to the valuation of the perpetual lookback game options with floating strikes in the appropriate diffusion-type extension of the Black–Merton–Scholes model.

This paper studies a novel Brownian functional defined as the supremum of a weighted average of the running Brownian range and its running reversal from extrema on the unit interval. We derive the Laplace transform for the squared reciprocal of this functional, which leads to explicit moment expressions that are new to the literature. We show that the proposed Brownian functional can be used to estimate the spot volatility of financial returns based on high-frequency price observations.

We derive some key extremal features for stationary kth-order Markov chains that can be used to understand how the process moves between an extreme state and the body of the process. The chains are studied given that there is an exceedance of a threshold, as the threshold tends to the upper endpoint of the distribution. Unlike previous studies with $k>1$, we consider processes where standard limit theory describes each extreme event as a single observation without any information about the transition to and from the body of the distribution. Our work uses different asymptotic theory which results in non-degenerate limit laws for such processes. We study the extremal properties of the initial distribution and the transition probability kernel of the Markov chain under weak assumptions for broad classes of extremal dependence structures that cover both asymptotically dependent and asymptotically independent Markov chains. For chains with $k>1$, the transition of the chain away from the exceedance involves novel functions of the k previous states, in comparison to just the single value, when $k=1$. This leads to an increase in the complexity of determining the form of this class of functions, their properties, and the method of their derivation in applications. We find that it is possible to derive an affine normalization, dependent on the threshold excess, such that non-degenerate limiting behaviour of the process, in the neighbourhood of the threshold excess, is assured for all lags. We find that these normalization functions have an attractive structure that has parallels to the Yule–Walker equations. Furthermore, the limiting process is always linear in the innovations. We illustrate the results with the study of kth-order stationary Markov chains with exponential margins based on widely studied families of copula dependence structures.