Showing papers in "arXiv: Probability in 2008"

The conditioned reconstructed process

Abstract: We investigate a neutral model for speciation and extinction, the constant rate birth-death process. The process is conditioned to have $n$ extant species today, we look at the tree distribution of the reconstructed trees-- i.e. the trees without the extinct species. Whereas the tree shape distribution is well-known and actually the same as under the pure birth process, no analytic results for the speciation times were known. We provide the distribution for the speciation times and calculate the expectations analytically. This characterizes the reconstructed trees completely. We will show how the results can be used to date phylogenies.

Function spaces and capacity related to a Sublinear Expectation: application to G-Brownian Motion Pathes

Abstract: In this paper we give some basic and important properties of several typical Banach spaces of functions of $G$-Brownian motion pathes induced by a sublinear expectation--G-expectation. Many results can be also applied to more general situations. A generalized version of Kolmogorov's criterion for continuous modification of a stochastic process is also obtained. The results can be applied to continuous time dynamic and coherent risk measures in finance in particular for path-dependence risky positions under situations of volatility model uncertainty.

Liouville Quantum Gravity and KPZ

Abstract: Consider a bounded planar domain D, an instance h of the Gaussian free field on D (with Dirichlet energy normalized by 1/(2\pi)), and a constant 0 < gamma < 2. The Liouville quantum gravity measure on D is the weak limit as epsilon tends to 0 of the measures \epsilon^{\gamma^2/2} e^{\gamma h_\epsilon(z)}dz, where dz is Lebesgue measure on D and h_\epsilon(z) denotes the mean value of h on the circle of radius epsilon centered at z. Given a random (or deterministic) subset X of D one can define the scaling dimension of X using either Lebesgue measure or this random measure. We derive a general quadratic relation between these two dimensions, which we view as a probabilistic formulation of the KPZ relation from conformal field theory. We also present a boundary analog of KPZ (for subsets of the boundary of D). We discuss the connection between discrete and continuum quantum gravity and provide a framework for understanding Euclidean scaling exponents via quantum gravity.

Random matrices: Universality of ESDs and the circular law

Abstract: Given an $n \times n$ complex matrix $A$, let $$\mu_{A}(x,y):= \frac{1}{n} |\{1\le i \le n, \Re \lambda_i \le x, \Im \lambda_i \le y\}|$$ be the empirical spectral distribution (ESD) of its eigenvalues $\lambda_i \in \BBC, i=1, ... n$. We consider the limiting distribution (both in probability and in the almost sure convergence sense) of the normalized ESD $\mu_{\frac{1}{\sqrt{n}} A_n}$ of a random matrix $A_n = (a_{ij})_{1 \leq i,j \leq n}$ where the random variables $a_{ij} - \E(a_{ij})$ are iid copies of a fixed random variable $x$ with unit variance. We prove a \emph{universality principle} for such ensembles, namely that the limit distribution in question is {\it independent} of the actual choice of $x$. In particular, in order to compute this distribution, one can assume that $x$ is real of complex gaussian. As a related result, we show how laws for this ESD follow from laws for the \emph{singular} value distribution of $\frac{1}{\sqrt{n}} A_n - zI$ for complex $z$. As a corollary we establish the Circular Law conjecture (in both strong and weak forms), that asserts that $\mu_{\frac{1}{\sqrt{n}} A_n}$ converges to the uniform measure on the unit disk when the $a_{ij}$ have zero mean.

On optimality of the barrier strategy in de Finetti's dividend problem for spectrally negative L\'{e}vy processes

Abstract: We consider the classical optimal dividend control problem which was proposed by de Finetti [Trans. XVth Internat. Congress Actuaries 2 (1957) 433--443]. Recently Avram, Palmowski and Pistorius [Ann. Appl. Probab. 17 (2007) 156--180] studied the case when the risk process is modeled by a general spectrally negative Levy process. We draw upon their results and give sufficient conditions under which the optimal strategy is of barrier type, thereby helping to explain the fact that this particular strategy is not optimal in general. As a consequence, we are able to extend considerably the class of processes for which the barrier strategy proves to be optimal.

Multivariate normal approximation using exchangeable pairs

Abstract: Since the introduction of Stein's method in the early 1970s, much research has been done in extending and strengthening it; however, there does not exist a version of Stein's original method of exchangeable pairs for multivariate normal approximation. The aim of this article is to fill this void. We present three abstract normal approximation theorems using exchangeable pairs in multivariate contexts, one for situations in which the underlying symmetries are discrete, and real and complex versions of a theorem for situations involving continuous symmetry groups. Our main applications are proofs of the approximate normality of rank k projections of Haar measure on the orthogonal and unitary groups, when k = o(n).

Product-form stationary distributions for deficiency zero chemical reaction networks

Abstract: We consider stochastically modeled chemical reaction systems with mass-action kinetics and prove that a product-form stationary distribution exists for each closed, irreducible subset of the state space if an analogous deterministically modeled system with mass-action kinetics admits a complex balanced equilibrium. Feinberg's deficiency zero theorem then implies that such a distribution exists so long as the corresponding chemical network is weakly reversible and has a deficiency of zero. The main parameter of the stationary distribution for the stochastically modeled system is a complex balanced equilibrium value for the corresponding deterministically modeled system. We also generalize our main result to some non-mass-action kinetics.

A New Central Limit Theorem under Sublinear Expectations

Abstract: We describe a new framework of a sublinear expectation space and the related notions and results of distributions, independence. A new notion of G-distributions is introduced which generalizes our G-normal-distribution in the sense that mean-uncertainty can be also described. W present our new result of central limit theorem under sublinear expectation. This theorem can be also regarded as a generalization of the law of large number in the case of mean-uncertainty.

Large deviations for infinite dimensional stochastic dynamical systems

Abstract: The large deviations analysis of solutions to stochastic differential equations and related processes is often based on approximation. The construction and justification of the approximations can be onerous, especially in the case where the process state is infinite dimensional. In this paper we show how such approximations can be avoided for a variety of infinite dimensional models driven by some form of Brownian noise. The approach is based on a variational representation for functionals of Brownian motion. Proofs of large deviations properties are reduced to demonstrating basic qualitative properties (existence, uniqueness and tightness) of certain perturbations of the original process.

Ising models on locally tree-like graphs

Abstract: We consider ferromagnetic Ising models on graphs that converge locally to trees. Examples include random regular graphs with bounded degree and uniformly random graphs with bounded average degree. We prove that the "cavity" prediction for the limiting free energy per spin is correct for any positive temperature and external field. Further, local marginals can be approximated by iterating a set of mean field (cavity) equations. Both results are achieved by proving the local convergence of the Boltzmann distribution on the original graph to the Boltzmann distribution on the appropriate infinite random tree.

Information Percolation with Equilibrium Search Dynamics

Abstract: We solve for the equilibrium dynamics of information sharing in a large population. Each agent is endowed with signals regarding the likely outcome of a random variable of common concern. Individuals choose the effort with which they search for others from whom they can gather additional information. When two agents meet, they share their information. The information gathered is further shared at subsequent meetings, and so on. Equilibria exist in which agents search maximally until they acquire sufficient information precision, and then minimally. A tax whose proceeds are used to subsidize the costs of search improves information sharing and can in some cases increase welfare. On the other hand, endowing agents with public signals reduces information sharing and can in some cases decrease welfare.

Stochastic equations of non-negative processes with jumps

Abstract: We study stochastic equations of non-negative processes with jumps. The existence and uniqueness of strong solutions are established under Lipschitz and non-Lipschitz conditions. The comparison property of two solutions are proved under suitable conditions. The results are applied to stochastic equations driven by one-sided Levy processes and those of continuous state branching processes with immigration.

Special invited paper. Large deviations

Abstract: This paper is based on Wald Lectures given at the annual meeting of the IMS in Minneapolis during August 2005. It is a survey of the theory of large deviations.

Rough evolution equations

Abstract: We generalize Lyons' rough paths theory in order to give a pathwise meaning to some nonlinear infinite-dimensional evolution equation associated to an analytic semigroup and driven by an irregular noise. As an illustration, we discuss a class of linear and nonlinear 1d SPDEs driven by a space--time Gaussian noise with singular space covariance and Brownian time dependence.

Brownian motion on the Sierpinski carpet

Abstract: We prove that, up to scalar multiples, there exists only one local regular Dirichlet form on a generalized Sierpinski carpet that is invariant with respect to the local symmetries of the carpet. Consequently for each such fractal the law of Brownian motion is uniquely determined and the Laplacian is well defined.

Stein's method and normal approximation of Poisson functionals

Abstract: We combine Stein's method with a version of Malliavin calculus on the Poisson space. As a result, we obtain explicit Berry-Ess\'een bounds in Central Limit Theorems (CLTs) involving multiple Wiener-It\^o integrals with respect to a general Poisson measure. We provide several applications to CLTs related to Ornstein-Uhlenbeck L\'evy processes.

The Alexander-Orbach conjecture holds in high dimensions

Abstract: We examine the incipient infinite cluster (IIC) of critical percolation in regimes where mean-field behavior has been established, namely when the dimension d is large enough or when d>6 and the lattice is sufficiently spread out. We find that random walk on the IIC exhibits anomalous diffusion with the spectral dimension d_s=4/3, that is, p_t(x,x)= t^{-2/3+o(1)}. This establishes a conjecture of Alexander and Orbach. En route we calculate the one-arm exponent with respect to the intrinsic distance.

Multivariate normal approximation using Stein's method and Malliavin calculus

Abstract: We combine Stein's method with Malliavin calculus in order to obtain explicit bounds in the multidimensional normal approximation (in the Wasserstein distance) of functionals of Gaussian fields. Our results generalize and refine the main findings by Peccati and Tudor (2005), Nualart and Ortiz-Latorre (2007), Peccati (2007) and Nourdin and Peccati (2007b, 2008); in particular, they apply to approximations by means of Gaussian vectors with an arbitrary, positive definite covariance matrix. Among several examples, we provide an application to a functional version of the Breuer-Major CLT for fields subordinated to a fractional Brownian motion.

On percolation in random graphs with given vertex degrees

Abstract: We study the random graph obtained by random deletion of vertices or edges from a random graph with given vertex degrees. A simple trick of exploding vertices instead of deleting them, enables us to derive results from known results for random graphs with given vertex degrees. This is used to study existence of giant component and existence of k-core. As a variation of the latter, we study also bootstrap percolation in random regular graphs. We obtain both simple new proofs of known results and new results. An interesting feature is that for some degree sequences, there are several or even infinitely many phase transitions for the k-core.

Large Deviations of Vector-valued Martingales in 2-Smooth Normed Spaces

Abstract: In this paper, we derive exponential bounds on probabilities of large deviations for “light tail” martingales taking values in finitedimensional normed spaces Our primary emphasis is on the case where the bounds are dimension-independent or nearly so We demonstrate that this is the case when the norm on the space can be approximated, within an absolute constant factor, by a norm which is differentiable on the unit sphere with a Lipschitz continuous gradient We also present various examples of spaces possessing the latter property

A decomposition of the bifractional Brownian motion and some applications

Abstract: In this paper we show a decomposition of the bifractional Brownian motion with parameters H,K into the sum of a fractional Brownian motion with Hurst parameter HK plus a stochastic process with absolutely continuous trajectories. Some applications of this decomposition are discussed.

On ergodicity of some Markov processes

Abstract: We formulate a criterion for the existence and uniqueness of an invariant measure for a Markov process taking values in a Polish phase space. In addition, weak-$^*$ ergodicity, that is, the weak convergence of the ergodic averages of the laws of the process starting from any initial distribution, is established. The principal assumptions are the existence of a lower bound for the ergodic averages of the transition probability function and its local uniform continuity. The latter is called the e-property. The general result is applied to solutions of some stochastic evolution equations in Hilbert spaces. As an example, we consider an evolution equation whose solution describes the Lagrangian observations of the velocity field in the passive tracer model. The weak-$^*$ mean ergodicity of the corresponding invariant measure is used to derive the law of large numbers for the trajectory of a tracer.

Limit theorems for additive functionals of a Markov chain

Abstract: Consider a Markov chain $\{X_n\}_{n\ge 0}$ with an ergodic probability measure $\pi$. Let $\Psi$ a function on the state space of the chain, with $\alpha$-tails with respect to $\pi$, $\alpha\in (0,2)$. We find sufficient conditions on the probability transition to prove convergence in law of $N^{1/\alpha}\sum_n^N \Psi(X_n)$ to a $\alpha$-stable law. ``Martingale approximation'' approach and ``coupling'' approach give two different sets of conditions. We extend these results to continuous time Markov jump processes $X_t$, whose skeleton chain satisfies our assumptions. If waiting time between jumps has finite expectation, we prove convergence of $N^{-1/\alpha}\int_0^{Nt} V(X_s) ds$ to a stable process. In the case of waiting times with infinite average, we prove convergence to a Mittag-Leffler process.

Stein's method and exact Berry--Esseen asymptotics for functionals of Gaussian fields

Abstract: We show how to detect optimal Berry--Esseen bounds in the normal approximation of functionals of Gaussian fields. Our techniques are based on a combination of Malliavin calculus, Stein's method and the method of moments and cumulants, and provide de facto local (one-term) Edgeworth expansions. The findings of the present paper represent a further refinement of the main results proven in Nourdin and Peccati [Probab. Theory Related Fields 145 (2009) 75--118]. Among several examples, we discuss three crucial applications: (i) to Toeplitz quadratic functionals of continuous-time stationary processes (extending results by Ginovyan [Probab. Theory Related Fields 100 (1994) 395--406] and Ginovyan and Sahakyan [Probab. Theory Related Fields 138 (2007) 551--579]); (ii) to ``exploding'' quadratic functionals of a Brownian sheet; and (iii) to a continuous-time version of the Breuer--Major CLT for functionals of a fractional Brownian motion.

Convexity, translation invariance and subadditivity for $g$-expectations and related risk measures

Abstract: Under the continuous assumption on the generator $g$, Briand et al. [Electron. Comm. Probab. 5 (2000) 101--117] showed some connections between $g$ and the conditional $g$-expectation $({\mathcal{E}}_g[\cdot|{\mathcal{F}}_t])_{t\in[0,T]}$ and Rosazza Gianin [Insurance: Math. Econ. 39 (2006) 19--34] showed some connections between $g$ and the corresponding dynamic risk measure $(\rho^g_t)_{t\in[0,T]}$. In this paper we prove that, without the additional continuous assumption on $g$, a $g$-expectation ${\mathcal{E}}_g$ satisfies translation invariance if and only if $g$ is independent of $y$, and ${\mathcal{E}}_g$ satisfies convexity (resp. subadditivity) if and only if $g$ is independent of $y$ and $g$ is convex (resp. subadditive) with respect to $z$. By these conclusions we deduce that the static risk measure $\rho^g$ induced by a $g$-expectation ${\mathcal{E}}_g$ is a convex (resp. coherent) risk measure if and only if $g$ is independent of $y$ and $g$ is convex (resp. sublinear) with respect to $z$. Our results extend the results in Briand et al. [Electron. Comm. Probab. 5 (2000) 101--117] and Rosazza Gianin [Insurance: Math. Econ. 39 (2006) 19--34] on these subjects.

Polynomial processes and their applications to mathematical Finance

Abstract: We introduce a class of Markov processes, called $m$-polynomial, for which the calculation of (mixed) moments up to order $m$ only requires the computation of matrix exponentials. This class contains affine processes, processes with quadratic diffusion coefficients, as well as L\'evy-driven SDEs with affine vector fields. Thus, many popular models such as exponential L\'evy models or affine models are covered by this setting. The applications range from statistical GMM estimation procedures to new techniques for option pricing and hedging. For instance, the efficient and easy computation of moments can be used for variance reduction techniques in Monte Carlo methods.

Efficient rare-event simulation for the maximum of heavy-tailed random walks

Abstract: Let $(X_n:n\geq 0)$ be a sequence of iid rv's with negative mean Set $S_0=0$ and define $S_n=X_1+ +X_n$ We propose an importance sampling algorithm to estimate the tail of $M=\max \{S_n:n\geq 0\}$ that is strongly efficient for both light and heavy-tailed increment distributions Moreover, in the case of heavy-tailed increments and under additional technical assumptions, our estimator can be shown to have asymptotically vanishing relative variance in the sense that its coefficient of variation vanishes as the tail parameter increases A key feature of our algorithm is that it is state-dependent In the presence of light tails, our procedure leads to Siegmund's (1979) algorithm The rigorous analysis of efficiency requires new Lyapunov-type inequalities that can be useful in the study of more general importance sampling algorithms

Matrix representation of the stationary measure for the multispecies TASEP

Abstract: In this work we construct the stationary measure of the N species totally asymmetric simple exclusion process in a matrix product formulation. We make the connection between the matrix product formulation and the queueing theory picture of Ferrari and Martin. In particular, in the standard representation, the matrices act on the space of queue lengths. For N >2 the matrices in fact become tensor products of elements of quadratic algebras. This enables us to give a purely algebraic proof of the stationary measure which we present for N =3.

Time discretization and Markovian iteration for coupled FBSDEs

Abstract: In this paper we lay the foundation for a numerical algorithm to simulate high-dimensional coupled FBSDEs under weak coupling or monotonicity conditions. In particular, we prove convergence of a time discretization and a Markovian iteration. The iteration differs from standard Picard iterations for FBSDEs in that the dimension of the underlying Markovian process does not increase with the number of iterations. This feature seems to be indispensable for an efficient iterative scheme from a numerical point of view. We finally suggest a fully explicit numerical algorithm and present some numerical examples with up to 10-dimensional state space.

A modified lookdown construction for the Xi-Fleming-Viot process with mutation and populations with recurrent bottlenecks

Abstract: Letbe a finite measure on the unit interval. A �-Fleming-Viot process is a probability measure valued Markov process which is dual to a coalescent with multiple collisions (�-coalescent) in analogy to the duality known for the classical Fleming-Viot process and Kingman's coalescent, whereis the Dirac measure in 0. We explicitly construct a dual process of the coalescent with simultaneous multi- ple collisions (�-coalescent) with mutation, the �-Fleming-Viot process with muta- tion, and provide a representation based on the empirical measure of an exchange- able particle system along the lines of Donnelly and Kurtz (1999). We establish pathwise convergence of the approximating systems to the limiting �-Fleming-Viot process with mutation. An alternative construction of the semigroup based on the