In this manuscript, we address open questions raised by Dieker and Yakir (2014), who proposed a novel method of estimating (discrete) Pickands constants $\mathcal{H}^\delta_\alpha$ using a family of estimators $\xi^\delta_\alpha(T)$, $T>0$, where $\alpha\in(0,2]$ is the Hurst parameter, and $\delta\geq0$ is the step size of the regular discretization grid. We derive an upper bound for the discretization error $\mathcal{H}_\alpha^0 - \mathcal{H}_\alpha^\delta$, whose rate of convergence agrees with Conjecture 1 of Dieker and Yakir (2014) in the case $\alpha\in(0,1]$ and agrees up to logarithmic terms for $\alpha\in(1,2)$. Moreover, we show that all moments of $\xi_\alpha^\delta(T)$ are uniformly bounded and the bias of the estimator decays no slower than $\exp\{-\mathcal CT^{\alpha}\}$, as T becomes large.

In this paper, we show that, with probability $1$ , a random Beltrami field exhibits chaotic regions that coexist with invariant tori of complicated topologies. The motivation to consider this question, which arises in the study of stationary Euler flows in dimension 3, is V.I. Arnold’s 1965 speculation that a typical Beltrami field exhibits the same complexity as the restriction to an energy hypersurface of a generic Hamiltonian system with two degrees of freedom. The proof hinges on the obtention of asymptotic bounds for the number of horseshoes, zeros and knotted invariant tori and periodic trajectories that a Gaussian random Beltrami field exhibits, which we obtain through a nontrivial extension of the Nazarov–Sodin theory for Gaussian random monochromatic waves and the application of different tools from the theory of dynamical systems, including Kolmogorov–Arnold–Moser (KAM) theory, Melnikov analysis and hyperbolicity. Our results hold both in the case of Beltrami fields on ${\mathbb {R}}^3$ and of high-frequency Beltrami fields on the 3-torus.

Gaussian process regression is widely used to model an unknown function on a continuous domain by interpolating a discrete set of observed design points. We develop a theoretical framework for proving new moderate deviations inequalities on different types of error probabilities that arise in GP regression. Two specific examples of broad interest are the probability of falsely ordering pairs of points (incorrectly estimating one point as being better than another) and the tail probability of the estimation error at an arbitrary point. Our inequalities connect these probabilities to the mesh norm, which measures how well the design points fill the space.

We consider a branching random walk on a d-ary tree of height n ($n \in \mathbb{N}$), in the presence of a hard wall which restricts each value to be positive, where d is a natural number satisfying $d\geqslant2$. We consider the behaviour of Gaussian processes with long-range interactions, for example the discrete Gaussian free field, under the condition that it is positive on a large subset of vertices. We observe a relation with the expected maximum of the processes. We find the probability of the event that the branching random walk is positive at every vertex in the nth generation, and show that the conditional expectation of the Gaussian variable at a typical vertex, under positivity, is less than the expected maximum by order of $\log n$.

We derive a new theoretical lower bound for the expected supremum of drifted fractional Brownian motion with Hurst index $H\in(0,1)$ over a (in)finite time horizon. Extensive simulation experiments indicate that our lower bound outperforms the Monte Carlo estimates based on very dense grids for $H\in(0,\tfrac{1}{2})$. Additionally, we derive the Paley–Wiener–Zygmund representation of a linear fractional Brownian motion in the general case and give an explicit expression for the derivative of the expected supremum at $H=\tfrac{1}{2}$ in the sense of Bisewski, Dȩbicki and Rolski (2021).

We study an optimal reinsurance problem for a diffusion model, in which the drift of the claim follows an Ornstein–Uhlenbeck process. The aim of the insurer is to maximize the expected exponential utility of its terminal wealth. We consider two cases: full information and partial information. Full information occurs when the insurer directly observes the drift; partial information occurs when the insurer observes only its claims. By applying stochastic control and by solving the corresponding Hamilton–Jacobi–Bellman equations, we find the value function and the optimal reinsurance strategy under both full and partial information. We determine a relationship between the value function and reinsurance strategy under full information with the value function and reinsurance strategy under partial information.

This paper addresses the asymptotic analysis of sojourn functionals of spatiotemporal Gaussian random fields with long-range dependence (LRD) in time, also known as long memory. Specifically, reduction theorems are derived for local functionals of nonlinear transformation of such fields, with Hermite rank $m\geq 1,$ under general covariance structures. These results are proven to hold, in particular, for a family of nonseparable covariance structures belonging to the Gneiting class. For $m=2,$ under separability of the spatiotemporal covariance function in space and time, the properly normalized Minkowski functional, involving the modulus of a Gaussian random field, converges in distribution to the Rosenblatt-type limiting distribution for a suitable range of values of the long-memory parameter.

This paper investigates spatial data on the unit sphere. Traditionally, isotropic Gaussian random fields are considered as the underlying mathematical model of the cosmic microwave background (CMB) data. We discuss the generalized multifractional Brownian motion and its pointwise Hölder exponent on the sphere. The multifractional approach is used to investigate the CMB data from the Planck mission. These data consist of CMB radiation measurements at narrow angles of the sky sphere. The results obtained suggest that the estimated Hölder exponents for different CMB regions do change from location to location. Therefore, the CMB temperature intensities are multifractional. The methodology developed is used to suggest two approaches for detecting regions with anomalies in the cleaned CMB maps.

This paper studies the joint tail asymptotics of extrema of the multi-dimensional Gaussian process over random intervals defined as $P(u)\;:\!=\; \mathbb{P}\{\cap_{i=1}^n (\sup_{t\in[0,\mathcal{T}_i]} ( X_{i}(t) +c_i t )>a_i u )\}$ , $u\rightarrow\infty$ , where $X_i(t)$ , $t\ge0$ , $i=1,2,\ldots,n$ , are independent centered Gaussian processes with stationary increments, $\boldsymbol{\mathcal{T}}=(\mathcal{T}_1, \ldots, \mathcal{T}_n)$ is a regularly varying random vector with positive components, which is independent of the Gaussian processes, and $c_i\in \mathbb{R}$ , $a_i>0$ , $i=1,2,\ldots,n$ . Our result shows that the structure of the asymptotics of P(u) is determined by the signs of the drifts $c_i$ . We also discuss a relevant multi-dimensional regenerative model and derive the corresponding ruin probability.

We consider a continuous Gaussian random field living on a compact set $T\subset \mathbb{R}^{d}$ . We are interested in designing an asymptotically efficient estimator of the probability that the integral of the exponential of the Gaussian process over T exceeds a large threshold u. We propose an Asmussen–Kroese conditional Monte Carlo type estimator and discuss its asymptotic properties according to the assumptions on the first and second moments of the Gaussian random field. We also provide a simulation study to illustrate its effectiveness and compare its performance with the importance sampling type estimator of Liu and Xu (2014a).

Latouche and Nguyen (2015b) constructed a sequence of stochastic fluid processes and showed that it converges weakly to a Markov-modulated Brownian motion (MMBM). Here, we construct a different sequence of stochastic fluid processes and show that it converges strongly to an MMBM. To the best of our knowledge, this is the first result on strong convergence to a Markov-modulated Brownian motion. Besides implying weak convergence, such a strong approximation constitutes a powerful tool for developing deep results for sophisticated models. Additionally, we prove that the rate of this almost sure convergence is $o(n^{-1/2} \log n)$ . When reduced to the special case of standard Brownian motion, our convergence rate is an improvement over that obtained by a different approximation in Gorostiza and Griego (1980), which is $o(n^{-1/2}(\log n)^{5/2})$ .

We study shot noise processes with cluster arrivals, in which entities in each cluster may experience random delays (possibly correlated), and noises within each cluster may be correlated. We prove functional limit theorems for the process in the large-intensity asymptotic regime, where the arrival rate gets large while the shot shape function, cluster sizes, delays, and noises are unscaled. In the functional central limit theorem, the limit process is a continuous Gaussian process (assuming the arrival process satisfies a functional central limit theorem with a Brownian motion limit). We discuss the impact of the dependence among the random delays and among the noises within each cluster using several examples of dependent structures. We also study infinite-server queues with cluster/batch arrivals where customers in each batch may experience random delays before receiving service, with similar dependence structures.

This note corrects an error in the definition of the rate function in Jacquier, Pakkanen, and Stone (2018) and slightly simplifies some proofs.

We introduce a definition of long range dependence of random processes and fields on an (unbounded) index space $T\subseteq \mathbb{R}^d$ in terms of integrability of the covariance of indicators that a random function exceeds any given level. This definition is specifically designed to cover the case of random functions with infinite variance. We show the value of this new definition and its connection to limit theorems via some examples including subordinated Gaussian as well as random volatility fields and time series.

We provide upper and lower bounds for the mean $\mathscr{M}(H)$ of $\sup_{t\geq 0} \{B_H(t) - t\}$ , with $B_H(\!\cdot\!)$ a zero-mean, variance-normalized version of fractional Brownian motion with Hurst parameter $H\in(0,1)$ . We find bounds in (semi-) closed form, distinguishing between $H\in(0,\frac{1}{2}]$ and $H\in[\frac{1}{2},1)$ , where in the former regime a numerical procedure is presented that drastically reduces the upper bound. For $H\in(0,\frac{1}{2}]$ , the ratio between the upper and lower bound is bounded, whereas for $H\in[\frac{1}{2},1)$ the derived upper and lower bound have a strongly similar shape. We also derive a new upper bound for the mean of $\sup_{t\in[0,1]} B_H(t)$ , $H\in(0,\frac{1}{2}]$ , which is tight around $H=\frac{1}{2}$ .

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.

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.

We determine the order of magnitude of $\mathbb{E}|\sum _{n\leqslant x}f(n)|^{2q}$, where $f(n)$ is a Steinhaus or Rademacher random multiplicative function, and $0\leqslant q\leqslant 1$. In the Steinhaus case, this is equivalent to determining the order of $\lim _{T\rightarrow \infty }\frac{1}{T}\int _{0}^{T}|\sum _{n\leqslant x}n^{-it}|^{2q}\,dt$.

In particular, we find that $\mathbb{E}|\sum _{n\leqslant x}f(n)|\asymp \sqrt{x}/(\log \log x)^{1/4}$. This proves a conjecture of Helson that one should have better than squareroot cancellation in the first moment and disproves counter-conjectures of various other authors. We deduce some consequences for the distribution and large deviations of $\sum _{n\leqslant x}f(n)$.

The proofs develop a connection between $\mathbb{E}|\sum _{n\leqslant x}f(n)|^{2q}$ and the $q$th moment of a critical, approximately Gaussian, multiplicative chaos and then establish the required estimates for that. We include some general introductory discussion about critical multiplicative chaos to help readers unfamiliar with that area.

Motivated by the classical De Bruijn's identity for the additive Gaussian noise channel, in this paper we consider a generalized setting where the channel is modelled via stochastic differential equations driven by fractional Brownian motion with Hurst parameter H ∈ (0, 1). We derive generalized De Bruijn's identity for Shannon entropy and Kullback–Leibler divergence by means of Itô's formula, and present two applications. In the first application we demonstrate its equivalence with Stein's identity for Gaussian distributions, while in the second application, we show that for H ∈ (0, 1/2], the entropy power is concave in time while for H ∈ (1/2, 1) it is convex in time when the initial distribution is Gaussian. Compared with the classical case of H = 1/2, the time parameter plays an interesting and significant role in the analysis of these quantities.

In this paper we introduce the multivariate Brownian semistationary (BSS) process and study the joint asymptotic behaviour of its realised covariation using in-fill asymptotics. First, we present a central limit theorem for general multivariate Gaussian processes with stationary increments, which are not necessarily semimartingales. Then, we show weak laws of large numbers, central limit theorems, and feasible results for BSS processes. An explicit example based on the so-called gamma kernels is also provided.