Let (Y, Z) denote the solution to a forward-backward stochastic differential equation (FBSDE). If one constructs a random walk $B^n$ from the underlying Brownian motion B by Skorokhod embedding, one can show $L_2$-convergence of the corresponding solutions $(Y^n,Z^n)$ to $(Y, Z).$ We estimate the rate of convergence based on smoothness properties, especially for a terminal condition function in $C^{2,\alpha}$. The proof relies on an approximative representation of $Z^n$ and uses the concept of discretized Malliavin calculus. Moreover, we use growth and smoothness properties of the partial differential equation associated to the FBSDE, as well as of the finite difference equations associated to the approximating stochastic equations. We derive these properties by probabilistic methods.

In a classical, continuous-time, optimal stopping problem, the agent chooses the best time to stop a stochastic process in order to maximise the expected discounted return. The agent can choose when to stop, and if at any moment they decide to stop, stopping occurs immediately with probability one. However, in many settings this is an idealistic oversimplification. Following Strack and Viefers we consider a modification of the problem in which stopping occurs at a rate which depends on the relative values of stopping and continuing: there are several different solutions depending on how the value of continuing is calculated. Initially we consider the case where stopping opportunities are constrained to be event times of an independent Poisson process. Motivated by the limiting case as the rate of the Poisson process increases to infinity, we also propose a continuous-time formulation of the problem where stopping can occur at any instant.

We establish a fundamental property of bivariate Pareto records for independent observations uniformly distributed in the unit square. We prove that the asymptotic conditional distribution of the number of records broken by an observation given that the observation sets a record is Geometric with parameter 1/2.

We characterize the second order subexponentiality of an infinitely divisible distribution on the real line under an exponential moment assumption. We investigate the asymptotic behaviour of the difference between the tails of an infinitely divisible distribution and its Lévy measure. Moreover, we study the second order asymptotic behaviour of the tail of the $t$th convolution power of an infinitely divisible distribution. The density version for a self-decomposable distribution on the real line without an exponential moment assumption is also given. Finally, the regularly varying case for a self-decomposable distribution on the half line is discussed.

We prove an extension of Pisier’s inequality (1986) with a dimension-independent constant for vector-valued functions whose target spaces satisfy a relaxation of the UMD property.

The ruin probability in the classical Brownian risk model can be explicitly calculated for both finite and infinite time horizon. This is not the case for the simultaneous ruin probability in the two-dimensional Brownian risk model. Relying on asymptotic theory, we derive in this contribution approximations for both simultaneous ruin probability and simultaneous ruin time for the two-dimensional Brownian risk model when the initial capital increases to infinity.

The limit behavior of partial sums for short range dependent stationary sequences (with summable autocovariances) and for long range dependent sequences (with autocovariances summing up to infinity) differs in various aspects. We prove central limit theorems for partial sums of subordinated linear processes of arbitrary power rank which are at the border of short and long range dependence.

The equilibrium properties of allocation algorithms for networks with a large number of nodes with finite capacity are investigated. Every node receives a flow of requests. When a request arrives at a saturated node, i.e. a node whose capacity is fully utilized, an allocation algorithm may attempt to reallocate the request to a non-saturated node. For the algorithms considered, the reallocation comes at a price: either extra capacity is required in the system, or the processing time of a reallocated request is increased. The paper analyzes the properties of the equilibrium points of the associated asymptotic dynamical system when the number of nodes gets large. At this occasion the classical model of Gibbens, Hunt, and Kelly (1990) in this domain is revisited. The absence of known Lyapunov functions for the corresponding dynamical system significantly complicates the analysis. Several techniques are used. Analytic and scaling methods are used to identify the equilibrium points. We identify the subset of parameters for which the limiting stochastic model of these networks has multiple equilibrium points. Probabilistic approaches are used to prove the stability of some of them. A criterion of exponential stability with the spectral gap of the associated linear operator of equilibrium points is also obtained.

In this paper, we solve exit problems for a one-sided Markov additive process (MAP) which is exponentially killed with a bivariate killing intensity $\omega(\cdot,\cdot)$ dependent on the present level of the process and the current state of the environment. Moreover, we analyze the respective resolvents. All identities are expressed in terms of new generalizations of classical scale matrices for MAPs. We also remark on a number of applications of the obtained identities to (controlled) insurance risk processes. In particular, we show that our results can be applied to the Omega model, where bankruptcy takes place at rate $\omega(\cdot,\cdot)$ when the surplus process becomes negative. Finally, we consider Markov-modulated Brownian motion (MMBM) as a special case and present analytical and numerical results for a particular choice of piecewise intensity function $\omega(\cdot,\cdot)$ .

The Gaussian polytope $\mathcal P_{n,d}$ is the convex hull of n independent standard normally distributed points in $\mathbb{R}^d$ . We derive explicit expressions for the probability that $\mathcal P_{n,d}$ contains a fixed point $x\in\mathbb{R}^d$ as a function of the Euclidean norm of x, and the probability that $\mathcal P_{n,d}$ contains the point $\sigma X$ , where $\sigma\geq 0$ is constant and X is a standard normal vector independent of $\mathcal P_{n,d}$ . As a by-product, we also compute the expected number of k-faces and the expected volume of $\mathcal P_{n,d}$ , thus recovering the results of Affentranger and Schneider (Discr. and Comput. Geometry, 1992) and Efron (Biometrika, 1965), respectively. All formulas are in terms of the volumes of regular spherical simplices, which, in turn, can be expressed through the standard normal distribution function $\Phi(z)$ and its complex version $\Phi(iz)$ . The main tool used in the proofs is the conic version of the Crofton formula.

We consider the Voronoi diagram generated by n independent and identically distributed $\mathbb{R}^{d}$ -valued random variables with an arbitrary underlying probability density function f on $\mathbb{R}^{d}$ , and analyze the asymptotic behaviors of certain geometric properties, such as the measure, of the Voronoi cells as n tends to infinity. We adapt the methods used by Devroye et al. (2017) to conduct a study of the asymptotic properties of two types of Voronoi cells: (1) Voronoi cells that have a fixed nucleus; (2) Voronoi cells that contain a fixed point. We show that the geometric properties of both types of cells resemble those in the case when the Voronoi diagram is generated by a homogeneous Poisson point process. Additionally, for the second type of Voronoi cells, we determine the limiting distribution, which is universal in all choices of f, of the re-scaled measure of the cells.

Both sequential Monte Carlo (SMC) methods (a.k.a. ‘particle filters’) and sequential Markov chain Monte Carlo (sequential MCMC) methods constitute classes of algorithms which can be used to approximate expectations with respect to (a sequence of) probability distributions and their normalising constants. While SMC methods sample particles conditionally independently at each time step, sequential MCMC methods sample particles according to a Markov chain Monte Carlo (MCMC) kernel. Introduced over twenty years ago in [6], sequential MCMC methods have attracted renewed interest recently as they empirically outperform SMC methods in some applications. We establish an $\mathbb{L}_r$ -inequality (which implies a strong law of large numbers) and a central limit theorem for sequential MCMC methods and provide conditions under which errors can be controlled uniformly in time. In the context of state-space models, we also provide conditions under which sequential MCMC methods can indeed outperform standard SMC methods in terms of asymptotic variance of the corresponding Monte Carlo estimators.

This paper investigates the random horizon optimal stopping problem for measure-valued piecewise deterministic Markov processes (PDMPs). This is motivated by population dynamics applications, when one wants to monitor some characteristics of the individuals in a small population. The population and its individual characteristics can be represented by a point measure. We first define a PDMP on a space of locally finite measures. Then we define a sequence of random horizon optimal stopping problems for such processes. We prove that the value function of the problems can be obtained by iterating some dynamic programming operator. Finally we prove via a simple counter-example that controlling the whole population is not equivalent to controlling a random lineage.

Let V be an n-set, and let X be a random variable taking values in the power-set of V. Suppose we are given a sequence of random coupons $X_1, X_2, \ldots $ , where the $X_i$ are independent random variables with distribution given by X. The covering time T is the smallest integer $t\geq 0$ such that $\bigcup_{i=1}^t X_i=V$ . The distribution of T is important in many applications in combinatorial probability, and has been extensively studied. However the literature has focused almost exclusively on the case where X is assumed to be symmetric and/or uniform in some way.

In this paper we study the covering time for much more general random variables X; we give general criteria for T being sharply concentrated around its mean, precise tools to estimate that mean, as well as examples where T fails to be concentrated and when structural properties in the distribution of X allow for a very different behaviour of T relative to the symmetric/uniform case.

Large samples from a light-tailed distribution often have a well-defined shape. This paper examines the implications of the assumption that there is a limit shape. We show that the limit shape determines the upper quantiles for a large class of random variables. These variables may be described loosely as continuous homogeneous functionals of the underlying random vector. They play an important role in evaluating risk in a multivariate setting. The paper also looks at various coefficients of tail dependence and at the distribution of the scaled sample points for large samples. The paper assumes convergence in probability rather than almost sure convergence. This results in an elegant theory. In particular, there is a simple characterization of domains of attraction.

We introduce a unified framework for solving first passage times of time-homogeneous diffusion processes. Using potential theory and perturbation theory, we are able to deduce closed-form truncated probability densities, as asymptotics or approximations to the original first passage time densities, for single-side level crossing problems. The framework is applicable to diffusion processes with continuous drift functions; in particular, for bounded drift functions, we show that the perturbation series converges. In the present paper, we demonstrate examples of applying our framework to the Ornstein–Uhlenbeck, Bessel, exponential-Shiryaev, and hypergeometric diffusion processes (the latter two being previously studied by Dassios and Li (2018) and Borodin (2009), respectively). The purpose of this paper is to provide a fast and accurate approach to estimating first passage time densities of various diffusion processes.

For a continuous-time random walk X = {Xt, t ⩾ 0} (in general non-Markov), we study the asymptotic behaviour, as t → ∞, of the normalized additive functional $c_t\int _0^{t} f(X_s)\,{\rm d}s$, t⩾ 0. Similarly to the Markov situation, assuming that the distribution of jumps of X belongs to the domain of attraction to α-stable law with α > 1, we establish the convergence to the local time at zero of an α-stable Lévy motion. We further study a situation where X is delayed by a random environment given by the Poisson shot-noise potential: $\Lambda (x,\gamma )= {\rm e}^{-\sum _{y\in \gamma } \phi (x-y)},$ where $\phi \colon \mathbb R\to [0,\infty )$ is a bounded function decaying sufficiently fast, and γ is a homogeneous Poisson point process, independent of X. We find that in this case the weak limit has both ‘quenched’ component depending on Λ, and a component, where Λ is ‘averaged’.

We combine the rough Heston model and the CIR (Cox–Ingersoll–Ross) interest rate together to form a rough Heston-CIR model, so that both the rough behaviour of the volatility and the stochastic nature of the interest rate can be captured. Despite the convoluted structure and non-Markovian property of this model, it still admits a semi-analytical pricing formula for European options, the implementation of which involves solving a fractional Riccati equation. The rough Heston-CIR model is more general, taking both the rough Heston model and the Heston-CIR model as special cases. The influence of rough volatility and stochastic interest rate is shown to be significant through numerical experiments.

For a zero-mean, unit-variance stationary univariate Gaussian process we derive the probability that a record at the time n, say $X_n$ , takes place, and derive its distribution function. We study the joint distribution of the arrival time process of records and the distribution of the increments between records. We compute the expected number of records. We also consider two consecutive and non-consecutive records, one at time j and one at time n, and we derive the probability that the joint records $(X_j,X_n)$ occur, as well as their distribution function. The probability that the records $X_n$ and $(X_j,X_n)$ take place and the arrival time of the nth record are independent of the marginal distribution function, provided that it is continuous. These results actually hold for a strictly stationary process with Gaussian copulas.