Probability (math.PR)

  • PDF
    Fyodorov, Hiary & Keating established an intriguing connection between the maxima of log-correlated processes and the ones of the Riemann zeta function on a short interval of the critical line. In particular, they suggest that the analogue of the free energy of the Riemann zeta function is identical to the one of the Random Energy Model in spin glasses. In this paper, the connection between spin glasses and the Riemann zeta function is explored further. We study a random model of the Riemann zeta function and show that its two-overlap distribution corresponds to the one of a one-step replica symmetry breaking (1-RSB) spin glass. This provides evidence that the local maxima of the zeta function are strongly clustered.
  • PDF
    We access the edge of Gaussian beta ensembles with one spike by analyzing high powers of the associated tridiagonal matrix models. In the classical cases beta=1, 2, 4, this corresponds to studying the fluctuations of the largest eigenvalues of additive rank one perturbations of the GOE/GUE/GSE random matrices. In the infinite-dimensional limit, we arrive at a one-parameter family of random Feynman-Kac type semigroups, which features the stochastic Airy semigroup of Gorin and Shkolnikov [13] as an extreme case. Our analysis also provides Feynman-Kac formulas for the spiked stochastic Airy operators, introduced by Bloemendal and Virag [6]. The Feynman-Kac formulas involve functionals of a reflected Brownian motion and its local times, thus, allowing to study the limiting operators by tools of stochastic analysis. We derive a first result in this direction by obtaining a new distributional identity for a reflected Brownian bridge conditioned on its local time at zero.
  • PDF
    This work concerns the Ornstein-Uhlenbeck type process associated to a positive self-similar Markov process $(X(t))_{t\geq 0}$ which drifts to $\infty$, namely $U(t):= {\rm e}^{-t}X({\rm e}^t-1)$. We point out that $U$ is always a (topologically) recurrent Markov process and identify its invariant measure in terms of the law of the exponential functional $\hat I := \int_0^\infty \exp(\hat\xi_s) {\rm d}s$, where $\hat\xi$ is the dual of the real-valued Lévy process $\xi$ related to $X$ by the Lamperti transformation. This invariant measure is infinite (i.e. $U$ is null-recurrent) if and only if $\xi_1\not \in L^1(\mathbb{P})$. In that case, we determine the family of Lévy processes $\xi$ for which $U$ fulfills the conclusions of the Darling-Kac theorem. Our approach relies crucially on another generalized Ornstein-Uhlenbeck process that can be associated to the Lévy process $\xi$, namely $V(t) := \exp(\xi_t)\left(\int_0^t \exp(-\xi_s){\rm d}s +V(0)\right)$, and properties of time-substitutions based on additive functionals.
  • PDF
    We analyze the metastable states near criticality of the bootstrap percolation on Galton-Watson trees. We find that, depending on the exact choice of the offspring distribution, it is possible to have several distinct metastable states, with varying scaling of their duration while approaching criticality.
  • PDF
    Let $\boldsymbol{X}(t)=(X_1(t),\ldots,X_d(t))$ be a Gaussian vector process and $g(t)$ be a continuous function. The asymptotics of distribution of $\left\|\boldsymbol{X}(t)\right\|_p$, the $L^p$ norm for Gaussian finite-dimensional vector, have been investigated in numerous literatures. In this contribution we are concerned with the exact tail asymptotics of $\left\|\boldsymbol{X}(t)\right\|^c_p,\ c>0, $ with trend $g(t)$ over $[0,T]$. Both scenarios that $\boldsymbol{X}(t)$ is locally stationary and non-stationary are considered. Important examples include $\sum_{i=1}^d \left|X_i(t)\right|+g(t)$ and chi-square processes with trend, i.e., $\sum_{i=1}^d X_i^2(t)+g(t)$. These results are of interest in applications in engineering, insurance and statistics, etc.
  • PDF
    This work focuses on recurrence and ergodicity of switching diffusions consisting of continuous and discrete components, in which the discrete component takes values in a countably infinite set and the rates of switching at current time depend on the value of the continuous component over an interval including certain past history. Sufficient conditions for recurrence and ergodicity are given. Moreover, the relationship between systems of partial differential equations and recurrence when the switching is past-independent is established under suitable conditions.
  • PDF
    For a general class of large non-Hermitian random block matrices $\mathbf{X}$ we prove that there are no eigenvalues away from a deterministic set with very high probability. This set is obtained from the Dyson equation of the Hermitization of $\mathbf{X}$ as the self-consistent approximation of the pseudospectrum.
  • PDF
    The convex hull $P_{n}$ of a Gaussian sample $X_{1},...,X_{n}$ in $R^{d}$ is a Gaussian polytope. We prove that the expected number of facets $E f_{d-1} (P_n)$ is monotonically increasing in $n$. Furthermore we prove this for random polytopes generated by uniformly distributed points in a $d$-dimensional ball.
  • PDF
    We present an elementary approach to concentration of disordered Hamiltonians. Assuming differentiability of the limiting free energy $F$ with respect to the inverse temperature $\beta$, we show that the Hamiltonian concentrates around the energy level $F'(\beta)$ under the free energy and Gibbs average.
  • PDF
    We consider uniform random permutations in proper substitution-closed classes and study their limiting behavior in the sense of permutons. The limit depends on the generating series of the simple permutations in the class. Under a mild sufficient condition, the limit is an elementary one-parameter deformation of the limit of uniform separable permutations, previously identified as the Brownian separable permuton. This limiting object is therefore in some sense universal. We identify two other regimes with different limiting objects. The first one is degenerate; the second one is nontrivial and related to stable trees. These results are obtained thanks to a characterization of permuton convergence by convergence of expected pattern densities. The limit of expected pattern densities is then computed by using the substitution tree encoding of permutations and performing singularity analysis on the tree series.
  • PDF
    The famous results of Komlós, Major and Tusnády (see [15] and [17]) state that it is possible to approximate almost surely the partial sums of size n of i.i.d. centered random variables in L p (p > 2) by a Wiener process with an error term of order o(n 1/p). Very recently, Berkes, Liu and Wu [3] extended this famous result to partial sums associated with functions of an i.i.d. sequence, provided a condition on a functional dependence measure in L p is satisfied. In this paper, we adapt the method of Berkes, Liu and Wu to partial sums of functions of random iterates. Taking advantage of the Markovian setting, we shall give new dependent conditions, expressed in terms of a natural coupling (in L $\infty$ or in L 1), under which the strong approximation result holds with rate o(n 1/p). As we shall see our conditions are well adapted to a large variety of models, including left random walks on GL d (R), contracting iterated random functions, autoregressive Lipschitz processes, and some ergodic Markov chains. We also provide some examples showing that our L 1-coupling condition is in some sense optimal.
  • PDF
    We study statistical model checking of continuous-time stochastic hybrid systems. The challenge in applying statistical model checking to these systems is that one cannot simulate such systems exactly. We employ the multilevel Monte Carlo method (MLMC) and work on a sequence of discrete-time stochastic processes whose executions approximate and converge weakly to that of the original continuous-time stochastic hybrid system with respect to satisfaction of the property of interest. With focus on bounded-horizon reachability, we recast the model checking problem as the computation of the distribution of the exit time, which is in turn formulated as the expectation of an indicator function. This latter computation involves estimating discontinuous functionals, which reduces the bound on the convergence rate of the Monte Carlo algorithm. We propose a smoothing step with tunable precision and formally quantify the error of the MLMC approach in the mean-square sense, which is composed of smoothing error, bias, and variance. We formulate a general adaptive algorithm which balances these error terms. Finally, we describe an application of our technique to verify a model of thermostatically controlled loads.
  • PDF
    In this paper, we consider isotropic and stationary max-stable, inverse max-stable and max-mixture processes $X=(X(s))\_{s\in\bR^2}$ and the damage function $\cD\_X^{\nu}= |X|^\nu$ with $0<\nu<1/2$. We study the quantitative behavior of a risk measure which is the variance of the average of $\cD\_X^{\nu}$ over a region $\mathcal{A}\subset \bR^2$. This kind of risk measure has already been introduced and studied for \verosome max-stable processes in \citekoch2015spatial. %\textcolorredIn this study, we generalised this risk measure to be applicable for several models: asymptotic dependence represented by max-stable, asymptotic independence represented by inverse max-stable and mixing between of them. We evaluated the proposed risk measure by a simulation study.
  • PDF
    In this paper, we are concerned with the stochastic SIS (susceptible-infected-susceptible) and SIR (susceptible-infected-recovered) models on high-dimensional lattices with random edge weights, where a susceptible vertex is infected by an infectious neighbor at rate proportional to the weight on the edge connecting them. All the edge weights are assumed to be i.i.d.. Our main result gives mean field limits for survival probabilities of the two models as the dimension grows to infinity, which extends the main conclusion given in \citeXue2017 for classic stochastic SIS model.
  • PDF
    We consider a non-attractive three state contact process on $\mathbb Z$ and prove that there exists a regime of survival as well as a regime of extinction. In more detail, the process can be regarded as an infection process in a dynamic environment, where non-infected sites are either healthy or passive. Infected sites can recover only if they have a healthy site nearby, whereas non-infected sites may become infected only if there is no healthy and at least one infected site nearby. The transition probabilities are governed by a global parameter $q$: for large $q$, the infection dies out, and for small enough $q$, we observe its survival. The result is obtained by a coupling to a discrete time Markov chain, using its drift properties in the respective regimes.
  • PDF
    The generalized Dickman distribution ${\cal D}_\theta$ with parameter $\theta>0$ is the unique distribution for a non-negative random variable $W$ such that $W=_d W^*$, where \begineqnarray W^*=_d U^1/\theta(W+1) \qquad (1) \endeqnarray with $U \sim {\cal U}[0,1]$ independent of $W$, and $=_d$ denoting equality in distribution. Members of this family appear in number theory, stochastic geometry, perpetuities and the study of algorithms. We obtain bounds in Wasserstein type distances between ${\cal D}_\theta$ and \begineqnarray W_n= \frac1n \sum_i=1^n Y_k B_k \qquad (2) \endeqnarray where $B_1,\ldots,B_n, Y_1, \ldots, Y_n$ are independent with $B_k \sim {\rm Ber}(1/k), E[Y_k]=k, {\rm Var}(Y_k)=\sigma_k^2$ and provide an application to the minimal directed spanning tree in $\mathbb{R}^2$, and also obtain such bounds when the Bernoulli variables in $(2)$ are replaced by Poissons. In addition, we broaden the class of generalized Dickman distributions by studying the fixed points of the transformation \begineqnarray* s(W^*)=_d U^1/\thetas(W+1) \endeqnarray* generalizing $(1)$, that allows the introduction of non-identity utility functions $s(\cdot)$ into Vervaat perpetuities. We obtain distributional bounds for recursive methods that can be used to simulate from this family, with special attention paid to the generalized Dickman case when $s(x)=x$.
  • PDF
    We consider sums of increments given by a functional of a stationary Markov chain. Letting $T$ be the first return time of the partial sums process to $(-\infty,0]$, under general assumptions we give determine the asymptotic behavior of the survival probability, $\mathbb{P}(T\ge t)\sim Ct^{-1/2}$ for an explicit constant $C$. Our analysis is based on a novel connection between the survival probability and the running maximum of the time-reversed process, and relies on a functional central limit theorem for Markov chains. Our result extends the classic theorem of Sparre Anderson on sums of mean zero and independent increments to the case of correlated increments. As applications, we recover known clustering results for the 3-color cyclic cellular automaton and the Greenberg-Hastings model in one dimension, and we prove a new clustering result for the 3-color firefly cellular automaton.
  • PDF
    Let $X_1,\ldots,X_n$ be a standard normal sample in $\mathbb R^d$. We compute exactly the expected volume of the Gaussian polytope $\mathrm{conv}[X_1,\ldots,X_n]$, the symmetric Gaussian polytope $\mathrm{conv}[\pm X_1,\ldots,\pm X_n]$, and the Gaussian zonotope $[0,X_1]+\ldots+[0,X_n]$ by exploiting their connection to the regular simplex, the regular crosspolytope, and the cube with the aid of Tsirelson's formula. The expected volumes of these random polytopes are given by essentially the same expressions as the intrinsic volumes and external angles of the regular polytopes. For all these quantities, we obtain asymptotic formulae which are more precise than the results which were known before. More generally, we determine the expected volumes of some heteroscedastic random polytopes including $ \mathrm{conv}[l_1X_1,\ldots,l_nX_n] $ and $ \mathrm{conv}[\pm l_1 X_1,\ldots, \pm l_n X_n], $ where $l_1,\ldots,l_n\geq 0$ are parameters, and the intrinsic volumes of the corresponding deterministic polytopes. Finally, we relate the $k$-th intrinsic volume of the regular simplex $S^{n-1}$ to the expected maximum of independent standard Gaussian random variables $\xi_1,\ldots,\xi_n$ given that the maximum has multiplicity $k$. Namely, we show that $$ V_k(S^n-1) = \frac (2\pi)^\frac k2 k! ⋅\lim_\varepsilon↓0 \varepsilon^1-k \mathbb E [\max{\xi_1,\ldots,\xi_n} 1_{\xi_(n) - \xi_(n-k+1)≤\varepsilon}], $$ where $\xi_{(1)} \leq \ldots \leq \xi_{(n)}$ denote the order statistics. A similar result holds for the crosspolytope if we replace $\xi_1,\ldots,\xi_n$ by their absolute values.
  • PDF
    In this paper we consider a large class of symmetric Markov processes $X=(X_t)_{t\ge0}$ on $\R^d$ generated by non-local Dirichlet forms, which include jump processes with small jumps of $\alpha$-stable-like type and with large jumps of super-exponential decay. Let $D\subset \R^d$ be an open (not necessarily bounded and connected) set, and $X^D=(X_t^D)_{t\ge0}$ be the killed process of $X$ on exiting $D$. We obtain explicit criterion for the compactness and the intrinsic ultracontractivity of the Dirichlet Markov semigroup $(P^{D}_t)_{t\ge0}$ of $X^D$. When $D$ is a horn-shaped region, we further obtain two-sided estimates of ground state in terms of jumping kernel of $X$ and the reference function of the horn-shaped region $D$.
  • PDF
    This work deals with the numerical approximation of backward stochastic differential equations (BSDEs). We propose a new algorithm which is based on the regression-later approach and the least squares Monte Carlo method. We give some conditions under which our numerical algorithm convergences and solve two practical experiments to illustrate its performance.
  • PDF
    We prove a martingale-coboundary representation for random fields with a completely commuting filtration. For random variables in L2 we present a necessary and sufficient condition which is a generalization of Heyde's condition for one dimensional processes from 1975. For Lp spaces with 2 ≤p < ∞we give a necessary and sufficient condition which extends Volny's result from 1993 to random fields and improves condition of El Machkouri and Giraudo from 2016 (arXiv:1410.3062). In application, new weak invariance principle and estimates of large deviations are found.
  • PDF
    We define a class of so-called thinnable ideals $\mathcal{I}$ on the positive integers which includes several well-known examples, e.g., the collection of sets with zero asymptotic density, sets with zero logarithmic density, and some summable ideals. It is shown that, given a sequence $(x_n)$ taking values in a topological space and a thinnable ideal $\mathcal{I}$, the set of cluster points of $(x_n)$ is equal to the set of cluster points of almost all its subsequences. Then, we obtain a characterization of ideal convergence, which improves the main result in [Trans. Amer. Math. Soc. \textbf347 (1995), 1811--1819].
  • PDF
    The goal of these notes is to fill some gaps in the literature about random walks in the Cauchy domain of attraction, which has been in many cases left aside because of its additional technical difficulties. We prove here several results in that case: a Fuk-Nagaev inequality and a local version of it ; a large deviation theorem ; two types of local large deviation theorems. We also derive two important applications of these results: a sharp estimate of the tail of the first ladder epochs, and renewal theorems -- extending standard renewal theorems to the case of random walks. Most of our techniques carry through to the case of random walks in the domain of attraction of an $\alpha$-stable law with $\alpha \in(0,2)$, so we also present results in that case, since many of them seem to be missing in the literature.
  • PDF
    The paper studies the higher-order absolute differences taken from progressive terms of time-homogenous binary Markov chains. Two theorems presented are the limiting theorems for these differences, when their order $k$ converges to infinity. Theorems 1 and 2 assert that there exist some infinite subsets $E$ of natural series such that $k$th order differences of every such chain converge to the equi-distributed random binary process as $k$ growth to infinity remaining on $E$. The chains are classified into two types and $E$ depend only on the type of a given chain. Two kinds of discrete capacities for subsets of natural series are defined, and in their terms such sets $E$ are described.
  • PDF
    We consider a Cox process with Poisson shot noise intensity which has been widely applied in insurance, finance, queue theory, statistic, and many other fields. Cox process is flexible because its intensity depends on not only the time but also a stochastic process, it can be considered as a two step randomization procedure. In this paper, we study precise deviations for shot noise Cox process using the recent mod-$\phi$ convergence method.

Recent comments

Mile Gu Nov 20 2015 05:04 UTC

Good question! There shouldn't be any contradiction with the correspondence principle. The reason here is that the quantum models are built to simulate the output behaviour of macroscopic, classical systems, and are not necessarily macroscopic themselves. When we compare quantum and classical comple

...(continued)
hong Nov 20 2015 00:40 UTC

Interesting results. But, just wondering, does it contradict to the correspondence principle?

Richard Kueng Jul 28 2015 07:01 UTC

fyi: our quantum implications are presented in Subsection 2.2 (pp 7-9).

Zoltán Zimborás May 28 2014 04:42 UTC

It's a bit funny to look at a formally verified proof of the CLT :), here it is online:
https://github.com/avigad/isabelle.

Piotr Migdał Apr 18 2014 18:43 UTC

A podcast summarizing this paper, by Geoff Engelstein: [The Dice Tower # 351 - Dealing with the Mockers (43:55 - 50:36)](http://dicetower.coolstuffinc.com/tdt-351-dealing-with-the-mockers), and [an alternative link on the BoardGameGeek](http://boardgamegeek.com/boardgamepodcastepisode/117163/tdt-351

...(continued)