Showing papers in "arXiv: Probability in 2006"

Multi-Dimensional G-Brownian Motion and Related Stochastic Calculus under G-Expectation

Abstract: We develop a notion of nonlinear expectation --G-expectation-- generated by a nonlinear heat equation with infinitesimal generator G. We first study multi-dimensional G-normal distributions. With this nonlinear distribution we can introduce our G-expectation under which the canonical process is a multi dimensional G-Brownian motion. We then establish the related stochastic calculus, especially stochastic integrals of Ito's type with respect to our G-Brownian motion and derive the related Ito's formula. We have also obtained the existence and uniqueness of stochastic differential equation under our G-expectation.

Processes on Unimodular Random Networks

Abstract: We investigate unimodular random networks. Our motivations include their characterization via reversibility of an associated random walk and their similarities to unimodular quasi-transitive graphs. We extend various theorems concerning random walks, percolation, spanning forests, and amenability from the known context of unimodular quasi-transitive graphs to the more general context of unimodular random networks. We give properties of a trace associated to unimodular random networks with applications to stochastic comparison of continuous-time random walk.

Asymptotic properties of realized power variations and related functionals of semimartingales

Abstract: This paper is concerned with the asymptotic behavior of sums of terms which are a test function f evaluated at successive increments of a discretely sampled semimartingale. Typically the test function is a power function (when the power is 2 we get the realized quadratic variation) . We prove a variety of ``laws of large numbers'', that is convergence in probability of these sums, sometimes after normalization. We also exhibit in many cases the rate of convergence, as well as associated central limit theorems.

Non-local Dirichlet Forms and Symmetric Jump Processes

Abstract: We consider the symmetric non-local Dirichlet form $(E, F)$ given by \[ E (f,f)=\int_{R^d} \int_{R^d} (f(y)-f(x))^2 J(x,y) dx dy \] with $F$ the closure of the set of $C^1$ functions on $R^d$ with compact support with respect to $E_1$, where $E_1 (f, f):=E (f, f)+\int_{R^d} f(x)^2 dx$, and where the jump kernel $J$ satisfies \[ \kappa_1|y-x|^{-d-\alpha} \leq J(x,y) \leq \kappa_2|y-x|^{-d-\beta} \] for $0<\alpha< \beta <2, |x-y|<1$. This assumption allows the corresponding jump process to have jump intensities whose size depends on the position of the process and the direction of the jump. We prove upper and lower estimates on the heat kernel. We construct a strong Markov process corresponding to $(E, F)$. We prove a parabolic Harnack inequality for nonnegative functions that solve the heat equation with respect to $E$. Finally we construct an example where the corresponding harmonic functions need not be continuous.

Empirical graph Laplacian approximation of Laplace–Beltrami operators: Large sample results

Abstract: Let ${M}$ be a compact Riemannian submanifold of ${{\bf R}^m}$ of dimension $\scriptstyle{d}$ and let ${X_1,...,X_n}$ be a sample of i.i.d. points in ${M}$ with uniform distribution. We study the random operators $$ \Delta_{h_n,n}f(p):=\frac{1}{nh_n^{d+2}}\sum_{i=1}^n K(\frac{p-X_i}{h_n})(f(X_i)-f(p)), p\in M $$ where ${K(u):={\frac{1}{(4\pi)^{d/2}}}e^{-\|u\|^2/4}}$ is the Gaussian kernel and ${h_n\to 0}$ as ${n\to\infty.}$ Such operators can be viewed as graph laplacians (for a weighted graph with vertices at data points) and they have been used in the machine learning literature to approximate the Laplace-Beltrami operator of ${M,}$ ${\Delta_Mf}$ (divided by the Riemannian volume of the manifold). We prove several results on a.s. and distributional convergence of the deviations ${\Delta_{h_n,n}f(p)-{\frac{1}{|\mu|}}\Delta_Mf(p)}$ for smooth functions ${f}$ both pointwise and uniformly in ${f}$ and ${p}$ (here ${|\mu|=\mu(M)}$ and ${\mu}$ is the Riemannian volume measure). In particular, we show that for any class ${{\cal F}}$ of three times differentiable functions on ${M}$ with uniformly bounded derivatives $$ \sup_{p\in M}\sup_{f\in F}\Big|\Delta_{h_n,p}f(p)-\frac{1}{|\mu|}\Delta_Mf(p)\Big|= O\Big(\sqrt{\frac{\log(1/h_n)}{nh_n^{d+2}}}\Big) a.s. $$ as soon as $$ nh_n^{d+2}/\log h_n^{-1}\to \infty and nh^{d+4}_n/\log h_n^{-1}\to 0, $$ and also prove asymptotic normality of ${\Delta_{h_n,p}f(p)-{\frac{1}{|\mu|}}\Delta_Mf(p)}$ (functional CLT) for a fixed ${p\in M}$ and uniformly in ${f}.$

Weak approximation of stochastic differential equations and application to derivative pricing

Abstract: The authors present a new simple algorithm to approximate weakly stochastic differential equations in the spirit of [1] and [2]. They apply it to the problem of pricing Asian options under the Heston stochastic volatility model, and compare it with other known methods. It is shown that the combination of the suggested algorithm and quasi-Monte Carlo methods makes computations extremely fast. [1] Shigeo Kusuoka, ``Approximation of Expectation of Diffusion Process and Mathematical Finance,'' Advanced Studies in Pure Mathematics, Proceedings of Final Taniguchi Symposium, Nara 1998 (T. Sunada, ed.), vol. 31 2001, pp. 147--165. [2] Terry Lyons and Nicolas Victoir, ``Cubature on Wiener Space,'' Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences 460 (2004), pp. 169--198.

A Delayed Black and Scholes Formula II

Abstract: This article is a sequel to [A.H.M.P]. In [A.H.M.P], we develop an explicit formula for pricing European options when the underlying stock price follows a non-linear stochastic delay equation with fixed delays in the drift and diffusion terms. In this article, we look at models of the stock price described by stochastic functional differential equations with variable delays. We present a class of examples of stock dynamics with variable delays that permit an explicit form for the option pricing formula. As in [A.H.M.P], the market is complete with no arbitrage. This is achieved through the existence of an equivalent martingale measure. In subsequent work, the authors intend to test the models in [A.H.M.P] and the present article against real market data.

The noise in the circular law and the Gaussian free field

Abstract: Fill an n x n matrix with independent complex Gaussians of variance 1/n. As n approaches infinity, the eigenvalues {z_k} converge to a sum of an H^1-noise on the unit disk and an independent H^{1/2}-noise on the unit circle. More precisely, for C^1 functions of suitable growth, the distribution of sum_{k=1}^n (f(z_k)-E f(z_k)) converges to that of a mean-zero Gaussian with variance given by the sum of the squares of the disk H^1 and the circle H^{1/2} norms of f. Moreover, with p_n the characteristic polynomial, log|p_n|- E log|p_n| tends to the planar Gaussian free field conditioned to be harmonic outside the unit disk. Finally, for polynomial test functions f, we prove that the limiting covariance structure is universal for a class of models including Haar distributed unitary matrices.

Recent advances in invariance principles for stationary sequences

Abstract: In this paper we survey some recent results on the central limit theorem and its weak invariance principle for stationary sequences. We also describe several maximal inequalities that are the main tool for obtaining the invariance principles, and also they have interest in themselves. The classes of dependent random variables considered will be martingale-like sequences, mixing sequences, linear processes, additive functionals of ergodic Markov chains.

Statistical Romberg extrapolation: A new variance reduction method and applications to option pricing

Abstract: We study the approximation of $\mathbb{E}f(X_T)$ by a Monte Carlo algorithm, where $X$ is the solution of a stochastic differential equation and $f$ is a given function. We introduce a new variance reduction method, which can be viewed as a statistical analogue of Romberg extrapolation method. Namely, we use two Euler schemes with steps $\delta$ and $\delta^{\beta},0<\beta<1$. This leads to an algorithm which, for a given level of the statistical error, has a complexity significantly lower than the complexity of the standard Monte Carlo method. We analyze the asymptotic error of this algorithm in the context of general (possibly degenerate) diffusions. In order to find the optimal $\beta$ (which turns out to be $\beta=1/2$), we establish a central limit type theorem, based on a result of Jacod and Protter for the asymptotic distribution of the error in the Euler scheme. We test our method on various examples. In particular, we adapt it to Asian options. In this setting, we have a CLT and, as a by-product, an explicit expansion of the discretization error.

Elements of Stochastic Calculus via Regularisation

Abstract: This paper first summarizes the foundations of stochastic calculus via regularization and constructs through this procedure Ito and Stratonovich integrals. In the second part, a survey and new results are presented in relation with finite quadratic variation processes, Dirichlet and weak Dirichlet processes.

Thresholds for virus spread on networks

Abstract: We study how the spread of computer viruses, worms, and other self-replicating malware is affected by the logical topology of the network over which they propagate. We consider a model in which each host can be in one of 3 possible states - susceptible, infected or removed (cured, and no longer susceptible to infection). We characterise how the size of the population that eventually becomes infected depends on the network topology. Specifically, we show that if the ratio of cure to infection rates is larger than the spectral radius of the graph, and the initial infected population is small, then the final infected population is also small in a sense that can be made precise. Conversely, if this ratio is smaller than the spectral radius, then we show in some graph models of practical interest (including power law random graphs) that the final infected population is large. These results yield insights into what the critical parameters are in determining virus spread in networks.

Critical Percolation Exploration Path and SLE(6): a Proof of Convergence

Abstract: It was argued by Schramm and Smirnov that the critical site percolation exploration path on the triangular lattice converges in distribution to the trace of chordal SLE(6). We provide here a detailed proof, which relies on Smirnov's theorem that crossing probabilities have a conformally invariant scaling limit (given by Cardy's formula). The version of convergence to SLE(6) that we prove suffices for the Smirnov-Werner derivation of certain critical percolation crossing exponents and for our analysis of the critical percolation full scaling limit as a process of continuum nonsimple loops.

Probabilistic approach for granular media equations in the non uniformly convex case

Abstract: We use here a particle system to prove a convergence result as well as a deviation inequality for solutions of granular media equation when the confinement potential and the interaction potential are no more uniformly convex. Proof is straightforward, simplifying deeply proofs of Carrillo-McCann-Villani \cite{CMV,CMV2} and completing results of Malrieu \cite{malrieu03} in the uniformly convex case. It relies on an uniform propagation of chaos property and a direct control in Wasserstein distance of solutions starting with different initial measures. The deviation inequality is obtained via a $T\_1$ transportation cost inequality replacing the logarithmic Sobolev inequality which is no more clearly dimension free.

Analysis of the Rosenblatt process

Abstract: We analyze {\em the Rosenblatt process} which is a selfsimilar process with stationary increments and which appears as limit in the so-called {\em Non Central Limit Theorem} (Dobrushin and Major (1979), Taqqu (1979)). This process is non-Gaussian and it lives in the second Wiener chaos. We give its representation as a Wiener-It\^o multiple integral with respect to the Brownian motion on a finite interval and we develop a stochastic calculus with respect to it by using both pathwise type calculus and Malliavin calculus.

Dynamics of trap models

Abstract: These notes cover one of the topics of the class given in the Les Houches Summer School ``Mathematical statistical physics'' in July 2005. The lectures tried to give a summary of the recent mathematical results about the long-time behaviour of dynamics of (mean-field) spin-glasses and other disordered media. We have chosen here to restrict the scope of these notes to the dynamics of trap models only, but to cover this topic in somewhat more depth.

Spatial birth and death processes as solutions of stochastic equations

Abstract: Spatial birth and death processes are obtained as solutions of a sys- tem of stochastic equations. The processes are required to be locally nite, but may involve an innite population over the full (noncompact) type space. Con- ditions are given for existence and uniqueness of such solutions, and for temporal and spatial ergodicity. For birth and death processes with constant death rate, a sub-criticality condition on the birth rate implies that the process is ergodic and converges exponentially fast to the stationary distribution.

Stein's method for concentration inequalities

Abstract: We introduce a version of Stein's method for proving concentration and moment inequalities in problems with dependence. Simple illustrative examples from combinatorics, physics, and mathematical statistics are provided.

Continuous-time GARCH processes

Abstract: A family of continuous-time generalized autoregressive conditionally heteroscedastic processes, generalizing the $\operatorname {COGARCH}(1,1)$ process of Kl\"{u}ppelberg, Lindner and Maller [J. Appl. Probab. 41 (2004) 601--622], is introduced and studied. The resulting $\operatorname {COGARCH}(p,q)$ processes, $q\ge p\ge 1$, exhibit many of the characteristic features of observed financial time series, while their corresponding volatility and squared increment processes display a broader range of autocorrelation structures than those of the $\operatorname {COGARCH}(1,1)$ process. We establish sufficient conditions for the existence of a strictly stationary nonnegative solution of the equations for the volatility process and, under conditions which ensure the finiteness of the required moments, determine the autocorrelation functions of both the volatility and the squared increment processes. The volatility process is found to have the autocorrelation function of a continuous-time autoregressive moving average process.

Scheduling control for queueing systems with many servers: asymptotic optimality in heavy traffic

Abstract: A multiclass queueing system is considered, with heterogeneous service stations, each consisting of many servers with identical capabilities. An optimal control problem is formulated, where the control corresponds to scheduling and routing, and the cost is a cumulative discounted functional of the system's state. We examine two versions of the problem: ``nonpreemptive,'' where service is uninterruptible, and ``preemptive,'' where service to a customer can be interrupted and then resumed, possibly at a different station. We study the problem in the asymptotic heavy traffic regime proposed by Halfin and Whitt, in which the arrival rates and the number of servers at each station grow without bound. The two versions of the problem are not, in general, asymptotically equivalent in this regime, with the preemptive version showing an asymptotic behavior that is, in a sense, much simpler. Under appropriate assumptions on the structure of the system we show: (i) The value function for the preemptive problem converges to $V$, the value of a related diffusion control problem. (ii) The two versions of the problem are asymptotically equivalent, and in particular nonpreemptive policies can be constructed that asymptotically achieve the value $V$. The construction of these policies is based on a Hamilton--Jacobi--Bellman equation associated with $V$.

Concentration inequalities for dependent Random variables via the martingale method

Abstract: The martingale method is used to establish concentration inequalities for a class of dependent random sequences on a countable state space, with the constants in the inequalities expressed in terms of certain mixing coefficients. Along the way, bounds are obtained on martingale differences associated with the random sequences, which may be of independent interest. As applications of the main result, concentration inequalities are also derived for inhomogeneous Markov chains and hidden Markov chains, and an extremal property associated with their martingale difference bounds is established. This work complements and generalizes certain concentration inequalities obtained by Marton and Samson, while also providing different proofs of some known results.

Error expansion for the discretization of Backward Stochastic Differential Equations

Abstract: We study the error induced by the time discretization of a decoupled forward-backward stochastic differential equations $(X,Y,Z)$. The forward component $X$ is the solution of a Brownian stochastic differential equation and is approximated by a Euler scheme $X^N$ with $N$ time steps. The backward component is approximated by a backward scheme. Firstly, we prove that the errors $(Y^N-Y,Z^N-Z)$ measured in the strong $L\_p$-sense ($p \geq 1$) are of order $N^{-1/2}$ (this generalizes the results by Zhang 2004). Secondly, an error expansion is derived: surprisingly, the first term is proportional to $X^N-X$ while residual terms are of order $N^{-1}$.

Exponential Approximation by Stein's Method and Spectral Graph Theory

Abstract: General Berry-Esseen bounds are developed for the exponential distri- bution using Stein's method and a new concentration inequality approach. As an application, a sharp error term is obtained for Hora's result that the spectrum of the Johnson graph has an exponential limit.

Partially Reflected Brownian Motion: A Stochastic Approach to Transport Phenomena

Abstract: Transport phenomena are ubiquitous in nature and known to be important for various scientific domains. Examples can be found in physics, electrochemistry, heterogeneous catalysis, physiology, etc. To obtain new information about diffusive or Laplacian transport towards a semi-permeable or resistive interface, one can study the random trajectories of diffusing particles modeled, in a first approximation, by the partially refle cted Brownian motion. This stochastic process turns out to be a convenient mathematical foundation for discrete, semi-continuous and continuous theoretical descriptions of diffusive transport. This paper presents an overview of these topics with a special emphasis on the close relation between stochastic processes with partial reflections and Laplacian transport phenomena. We give selected examples of these phenomena followed by a brief introduction to the partially reflected Brownian motion and related probabilistic topics (e.g., local t ime process and spread harmonic measure). A particular attention is paid to the use of the Dirichlet-to-Neumann operator. Some practical consequences and further perspectives are discussed.

Cube root fluctuations for the corner growth model associated to the exclusion process

Abstract: We study the last-passage growth model on the planar integer lattice with exponential weights. With boundary conditions that represent the equilibrium exclusion process as seen from a particle right after its jump we prove that the variance of the last-passage time in a characteristic direction is of order t^{2/3}. With more general boundary conditions that include the rarefaction fan case we show that the last-passage time fluctuations are still of order t^{1/3}, and also that the transversal fluctuations of the maximal path have order t^{2/3}. We adapt and then build on a recent study of Hammersley's process by Cator and Groeneboom, and also utilize the competition interface introduced by Ferrari, Martin and Pimentel. The arguments are entirely probabilistic, and no use is made of the combinatorics of Young tableaux or methods of asymptotic analysis.

The mixing time of the giant component of a random graph

Abstract: We show that the total variation mixing time of the simple random walk on the giant component of supercritical Erdos-Renyi graphs is log^2 n. This statement was only recently proved, independently, by Fountoulakis and Reed. Our proof follows from a structure result for these graphs which is interesting in its own right. We show that these graphs are "decorated expanders" - an expander glued to graphs whose size has constant expectation and exponential tail, and such that each vertex in the expander is glued to no more than a constant number of decorations.

Stochastic Differential Equations: A Wiener Chaos Approach

Abstract: A new method is described for constructing a generalized solution for stochastic differential equations. The method is based on the Cameron-Martin version of the Wiener Chaos expansion and provides a unified framework for the study of ordinary and partial differential equations driven by finite- or infinite-dimensional noise with either adapted or anticipating input. Existence, uniqueness, regularity, and probabilistic representation of this Wiener Chaos solution is established for a large class of equations. A number of examples are presented to illustrate the general constructions. A detailed analysis is presented for the various forms of the passive scalar equation and for the first-order It\^{o} stochastic partial differential equation. Applications to nonlinear filtering if diffusion processes and to the stochastic Navier-Stokes equation are also discussed.

The configurational measure on mutually avoiding SLE paths

Abstract: We define multiple chordal SLEs in a simply connected domain by considering a natural configurational measure on paths. We show how to construct these measures so that they are conformally covariant and satisfy certain boundary perturbation and Markov properties, as well as a cascade relation. As an example of our construction, we derive the scaling limit of Fomin's identity in the case of two paths directly; that is, we prove that the probability that an SLE(2) and a Brownian excursion do not intersect can be given in terms of the determinant of the excursion hitting matrix. Finally, we define the lambda-SAW, a one-parameter family of measures on self-avoiding walks on Z^2.

Overshoots and undershoots of L\'{e}vy processes

Abstract: We obtain a new fluctuation identity for a general Levy process giving a quintuple law describing the time of first passage, the time of the last maximum before first passage, the overshoot, the undershoot and the undershoot of the last maximum. With the help of this identity, we revisit the results of Kluppelberg, Kyprianou and Maller [Ann. Appl. Probab. 14 (2004) 1766--1801] concerning asymptotic overshoot distribution of a particular class of Levy processes with semi-heavy tails and refine some of their main conclusions. In particular, we explain how different types of first passage contribute to the form of the asymptotic overshoot distribution established in the aforementioned paper. Applications in insurance mathematics are noted with emphasis on the case that the underlying Levy process is spectrally one sided.

Deviation bounds for additive functionals of Markov process

Abstract: In this paper we derive non asymptotic deviation bounds for $$¶_ u (|\frac 1t \int_0^t V(X_s) ds - \int V d\mu | \geq R)$$ where $X$ is a $\mu$ stationary and ergodic Markov process and $V$ is some $\mu$ integrable function. These bounds are obtained under various moments assumptions for $V$, and various regularity assumptions for $\mu$. Regularity means here that $\mu$ may satisfy various functional inequalities (F-Sobolev, generalized Poincare etc...).