Dynamical systems perturbed by big jumps: meta stability and non local search
|
|
- Laureen Newton
- 5 years ago
- Views:
Transcription
1 Dynamical systems perturbed by big jumps: meta stability and non local search Ilya Pavlyukevich Humboldt Universität zu Berlin / Universität Heidelberg Institut für Mathematik Workshop III: Transport Systems in Geography, Geosciences, and Networks May 5 9, 28, IPAM, UCLA Supported by the DFG research projekt Stochastic Dynamics of Climate States
2 . Motivation Greenland ice-core data allows to reconstruct Earth s climate up to 2. years before present. International projects: GRIP (328 m), GISP2 ( m), NGRIP ( m)
3 2. Paleo proxy data. Dansgaard Oeschger events Paleo data proxies: oxygen isotopes, dust, volcanic markers etc. Global climate during the last glacial ( 2 b. p.) has experienced at least 2 abrupt and large amplitude shifts (Dansgaard Oeschger events) T Age (thousands of years before present) rapid warming by 5 C within at most a few decades plateau phase with slow cooling lasting several centuries rapid drop to cold stadial conditions Simulations: Ganopolsky&Rahmstorf, Potsdam Institute for Climate Impact Research
4 3. Paleo proxy data: detailed look. The calcium (Ca) signal from the GRIP ice core: about 8, data points from kyr to 9 kyr before present. Typical interjump time: 2 years, mean waiting time 47 years What triggers the transitions? Langevin equation for climate dynamics dx(t) = U (X(t)) dt + ε dl(t) U double-well potential, wells correspond to the climate states. P. Ditlevsen (Geophys. Res. Lett. 999): histogram analysis of the noise. Noise L has α stable component with α.75.
5 4. The random perturbation L. Lévy processes L is a Lévy process: independent increments L(s) L(t) L(s) B is a Brownian motion: independent increments stationary increments L(t) L(s) d = L(t s) stationary increments L() = stochastically continuous (no fixed jumps) B() = B(t) B(s) N (, t s) right continuous with left limits continuous paths Description via Fourier transform: ϕ t (λ) = Ee iλl(t) L(t + s) = L(t) + (L(t + s) L(t)) d = L(t) + L (s) ϕ t+s (λ) = ϕ s (λ) ϕ t (λ) ϕ t (λ) = Ee iλl(t) = e tφ(λ)
6 Poisson process, intensity c > Ee iλl(t) = e iλk e ct(ct)k = e c(eiλ ) k! k= = exp 5. Examples [ t ] (e iλy )cδ (dy) Compound Poisson [ process, c > ], jumps σ Ee iλl(t) = exp t (e iλy )cσ(dy) General case, Lévy Hinchin Formula Ee iλl(t) [ ( = exp t d 2 λ2 + iµλ + )].5 (e iλy iλy I{ y })ν(dy) d, µ R, ν({}) =, (y 2 )ν(dy) <
7 6. α stable Lévy processes (Lévy flights) Jump measure ν(y) = y +α, α (, 2). ( Ee iλl(t) = exp t (e iλy iλy I{ y }) dy y +α Markov process with a generator Af(x) = [f(x + y) f(x) yf (x)i{ y }] Af(x) = ( ) α/2 f(x) L is a purely jump process. ) dy y +α, fractional Laplacian Control of jumps: J(t) = L(t) L(t ) the jump size at time t >. ( dy N(t, A) = {s (, t] : J(s) A} Poisson t For example: {s (, t] : J(s) } Poisson ( t A y +α ) dy y +α )
8 7. α stable Lévy processes II α =.75 α =.75 Length of sample paths: finite for α (, ) and infinite for α [, 2). Self-similarity: (L(ct)) t = (c /α L(t)) t. Hurst parameter H = α 2. Heavy tails: P(L(t) u) c uα, u. Explicit form of the Fourier transform Ee iλl(t) = e c(α) λ αt, α = α = 2 Cauchy process Brownian motion π + x 2 e x2 2 2π
9 8. Object of study. Simple system with Lévy perturbation Small noise (ε ) asymptotics of solutions of SDE X ε x(t) = x U (X ε x(s)) ds + εl(t), ε. L α stable (symmetric) Lévy process (maybe + Brownian motion) Regular n well potential U(x) (smooth, U (m i ), U (s i ), U (x) > x +δ s sn, x ± ) x Meta stable behaviour Transitions between the wells m m n Regular one well potential (smooth, U () > ) U(x) Exit time: σ x (ε) = inf{t : X ε x(t) / [ b, a]} a, b < (b = ) b a x P. Imkeller & I. Pavlyukevich Stoch. Proc. Appl. 6, 26; J. Phys. A: Math. Gen. 39, 26; ESAIM: P&S, 28
10 9. What is known Freidlin, Wentzell (Random Perturbations of Dynamical Systems, 979): Gaussian perturbations: ˆX ε x(t) = x U ( ˆX ε x(s)) ds + εw (t), ε. Perturbation: εw standard Brownian motion of small amplitude (in R d ). Locally infinitely divisible perturbations leading to Gaussian εl ε with jump part [ ] A ε ν(x, dy) f(x) = f(x + εy) f(x) εy, f(x) R d \{} ε Lévy measure ν has all exponential moments. For example, εl ε (t) = επ /ε (t) t, π standard Poisson process. Effects of heavy tails Godovanchuk (Theor. Prob. Appl. 26, 982): Large deviations for Markov processes with heavy jumps. Samorodnitsky, Grigoriou (Stoch. Proc. Appl. 5, 69 97, 23): Tails of solutions of SDEs driven by Lévy processes with power tails.
11 . First exit: Gaussian vs. Lévy U(x) U( b)>u(a) h=u(a) U(x) b ˆX ε x(t) = x x a U ( ˆX ε x(s)) ds + εw (t) Large deviations (Freidlin Wentzell): b X ε x(t) = x x a U (X ε x(s)) ds + εl(t) P(e (2h δ)/ε2 < ˆσ x (ε) < e (2h+δ)/ε2 ) Mean exit time (Kramers, Day, Bovier): ε ( ) π 2h Eˆσ x (ε) U (a) U () exp ε 2 Exponential exit (Day, Bovier) ( ) ˆσx (ε) P Eˆσ x (ε) > u exp ( u) P( ε α δ < σ x (ε) < ε α+δ ) P Eσ x (ε) α [ ε α a α + ] b α ( ) σx (ε) Eσ x (ε) > u exp ( u)
12 . Transitions For small > denote B i = {y : y m i } and τ i x(ε) = inf{t : X ε x(t) j i B j } τ s U(x) s n x m n Theorem. For x B i the following holds as ε : m P(ε α τ i x(ε) > u) e q iu, Eτ i x(ε) ε α q i, q i = P(X ε (τ i x(ε)) B j ) q ij q i, q ij = R\(s i,s i ) (s j,s j ) dy y m i +α, dy y m i +α, i j.
13 Theorem Meta stability Let x (s j, s j ) and t >. Then X ε x ( ) t ε α d Y mj (t), where Y is a Markov chain on {m,..., m n } with a generator Q = (q ij ), q ii = q i. Remark. Y has the invariant measure π(dy) = π j >, n π j δ mj (dy), j= where Q T π =.
14 3. Meta stable behaviour. Gaussian case ˆX ε x(t) = x U ( ˆX ε x(s)) ds + εw (t) Different life times for each well: Eτ i x(ε) exp(v i /ε 2 ) n critical exponents = λ < λ < < λ n such that for Λ (λ i, λ i+ ) the process ˆX ( x ε te Λ/ε2) converges to µ(λ, x) Double well potential: critical time scale: λ ε 2π U () U (m 2 ) exp(2h/ε2 ) Convergence of fin. dim. distrib.: (Kipnis, Newman; Mathieu) ˆX ε (λ ε t) Ŷ (t) h < H m U(x) m 2 H h x Generator of Ŷ ( ), and Ŷ () = { m, if x <, m 2, if x >.
15 4. Doubts and questions X ε x(t) = x U (X ε x(s)) ds + εl(t), ε. Why α stable (non-gaussian, heavy-tail) noise? Why small noise?
16 5. A natural system with small noise Is there a Gaussian system with a priori small noise? Simulated annealing (Kirkpatrick et al., Geman&Geman, Černy, 985): Ẑ(t) = z σ(t) temperature, σ(t) as t. U (Ẑ(s)) ds + σ(s)dw (s), ANNEALING: gradual cooling of steel (copper, glass etc.) in order to induce softenness, to relieve internal stresses, to refine the crystalline structure. I. Pavlyukevich: Stoch. Proc. Appl., 28; J. Phys. A: Math. and Theor., 27
17 6. Gaussian diffusion: long time behaviour Stochastic optimisation: look for a global minimum m of U. ˆX(t) = x U ( ˆX(s)) ds + εw (t) Generator A ε f = ε2 2 f U f.5 Invariant measure: A εµ = µ ε (dx) = c(ε)e 2U(x)/ε2.5 dx µ ε δ m, ε The spectrum { Λ ε k, Ψε k } k of A ε is discrete. Λ ε =, Ψ ε (x) = Λ ε exp( Θ/ε 2 ) SPECTRAL GAP (Friedman, Day, Bovier et al.) E x f( ˆX ε (t)) = f, L 2 (µ ε )(x) + e Λε t f, Ψ ε L 2 (µ ε )Ψ ε (x) + Law( ˆX ε (t)) µ ε, t
18 Ẑ z (t) = z INTUITION: 7. Gaussian simulated annealing U (Ẑz(s)) ds + Ẑ(t) ˆX σ(t) (t) FOURIER EXPANSION : σ(s)dw (s), σ(t) = ( ) /2 θ. ln(λ + t) E,z f(ẑ(t)) f, L 2 (µ σ(t) ) + e Λσ(t) t f, Ψ σ(t) L 2 (µ σ(t) )Ψ σ(t) + CONVERGENCE: µ σ(t) δ m, Λ σ(t) t e Θ/σ(t)2 t = t (t + λ) Θ/θ {, θ > Θ,, θ < Θ, t +, θ > Θ Z(t) m. This works, see Chiang&Hwang&Sheu, Holley&Kusuoka&Stroock and others HOWEVER: very slow cooling, local search, unknown critical rate Θ...
19 8. Modifications: big jumps Szu and Hartley (Phys. Lett. A 22, 987): fast simulated annealing Discrete schema with Cauchy jumps ξ k X kh = X (k )h + U (X (k )h )h + ξ k (k )h U = U(X kh ) U(X (k )h ) Metropolis acceptance probability: Accept always if U Accept with probability exp ( U kh) if U > Non-local search Fast algebraic cooling Convergence??? Probably works??? Applications to image recognition
20 9. Simulated annealing with α stable processes Analogously to Gaussian simulated annealing: Z(t) = z U (Z(s)) ds + dl(s) (λ + s) θ, λ >, θ >. Convergence of Z(t) as t? Time-homogeneous case: (s j,s j ) ( ) t X ε ε α Y (t) Y Markov chain on {m,..., m n } with generator Q = (q ij ) dy q ij = y m i +α, q dy ii = q i = y m i +α R\(s i,s i ) Law(Y (t)) π, Q π = Z(t) X /(λ+t)θ (t) Y ( t ) (t + λ) αθ { π, αθ <, no convergence, αθ >.
21 Slow cooling: convergence Theorem 3. [slow cooling] Let αθ <. Then Eτ i (λ) λαθ q i, P(Z(τ i (λ)) B j ) q ij q i, λ, τ s U(x) s n m n x Z(t) π, t. m
22 2. Fast cooling: trapping Theorem 4. [trapping] Let αθ >. Then P(τ i (λ) < ) = O(λ αθ ), λ, Eτ i (λ) =
23 22. Gaussian vs. α stable simulated annealing Simulated annealing with α stable process: allows to determine the measure π, i.e. the sizes of potential wells the search is non-local polynomially fast cooling BUT: how to find the global minimum of U?
24 23. Simulated annealing with stable-like processes In progress: I.P., J. Comp. Phys., 27 V (t) = v U (V (s )) ds + dh(v (s ), s) (λ + s) θ, λ >, θ > H a stable-like process with the jump measure dy ν x (dy) = y +α(x) U(x) A α(x) s s n x a s s n x E m n E m n m α(x) = { A, U(x) E, trapping, a, U(x) > E, exit, m < a < A < 2, aθ < < Aθ.
25 24. An example of non-local search The lowest minmum: U(4.9, 9.9) =.46 The second { lowest minimum: U( 9.7,.) =.85.8, if U(x), α(x) = θ =.75.., if U(x) >, Look for a local minimum: U(m) V k = V k U(V k )h ξh k (V k ) (λ+(k )h) θ k 2 6, λ = 4, V [ 2, 2] 2, h =. Success ratio: 96 out of.
26 25. Numerical example in R Stable-like simulated annealing a =., A =.8, θ =.75 λ N first k = 5 4 k = 5 5 N k k N k k Gaussian SDE θ = 3 λ N first k = 5 4 k = 5 5 N k k N k k Fast simulated annealing λ N first k = 5 4 k = 5 5 N k k N k k Gaussian simulated annealing θ = 3 λ N first k = 5 4 k = 5 5 N k k N k k
27 26. Numerical example in R 4. Shekel s function S,4 (y) = i= c i + y a i 2, y R4. Stable-like simulated annealing a =.2, A =.9, θ =.6 λ N first k = 5 4 k = 5 5 N k k N k k Gaussian SDE θ = λ N first k = 5 4 k = 5 5 N k k N k k Fast simulated annealing λ N first k = 5 4 k = 5 5 N k k N k k Gaussian simulated annealing θ = λ N first k = 5 4 k = 5 5 N k k N k k
28 27. Numerical example in R 6. Hartman s function 4 6 H 4,6 (y) = c i exp a ij (y i p ij ) 2,, y R 6. i= j= Stable-like simulated annealing a =.5, A =.9, θ =.6 λ N first k = 5 4 k = 5 5 N k k N k k Gaussian SDE θ = 6 λ N first k = 5 4 k = 5 5 N k k N k k Fast simulated annealing λ N first k = 5 4 k = 5 5 N k k N k k Gaussian simulated annealing θ = 6 λ N first k = 5 4 k = 5 5 N k k N k k
29 Pictures Page Slides/greenland.gif Page 2 From: Ganopolski, A. and Rahmstorf, S., Abrupt glacial climate changes due to stochastic resonance. Physical Review Letters 88(3), 385+, 22. Page 3 From: Ditlevsen, P. D., Observation of α stable noise induced millenial climate changes from an ice record, Geophysical Research Letters 26(), , 999. Page Page 8 From: H. Szu, Automated fault recognition by image correlation neutral network technique. IEEE Transactions on Industrial Electronics, 4(2):97 28, 993. References. Imkeller, P. and Pavlyukevich, I. First exit times of SDEs driven by stable Lévy processes, Stochastic Processes and their Applications 6 (4), 6 642, Imkeller, P. and Pavlyukevich, I. Metastable behaviour of small noise Lévy driven diffusions, ESAIM: Probability and Statistics, Imkeller, P. and Pavlyukevich, I. Lévy flights: transitions and meta stability, Journal of Physics A: Mathematical and General 39, L237 L246, Pavlyukevich, I. Simulated annealing for Lévy driven jump diffusions, Stochastic Processes and their Applications 8, 7 5, Pavlyukevich, I. Cooling down Lévy flights, Journal of Physics A: Mathematical and Theoretical 4, , Pavlyukevich, I. Lévy flights, non-local search and simulated annealing, Journal of Computational Physics 226 (2), , 27 Ilya Pavlyukevich, Humboldt Universität zu Berlin, Institut für Mathematik, Rudower Chaussee 25, 2489 Berlin, Germany pavljuke@math.hu-berlin.de pavljuke 28
Dynamics of reaction-diffusion equations with small Lévy noise
Dynamics of reaction-diffusion equations with small Lévy noise A. Debussche J. Gairing C. Hein M. Högele P. Imkeller I. Pavlyukevich ENS Cachan, U de los Andes Bogota, HU Berlin, U Jena Banff, January
More informationResidence-time distributions as a measure for stochastic resonance
W e ie rstra ß -In stitu t fü r A n g e w a n d te A n a ly sis u n d S to ch a stik Period of Concentration: Stochastic Climate Models MPI Mathematics in the Sciences, Leipzig, 23 May 1 June 2005 Barbara
More informationSystems Driven by Alpha-Stable Noises
Engineering Mechanics:A Force for the 21 st Century Proceedings of the 12 th Engineering Mechanics Conference La Jolla, California, May 17-20, 1998 H. Murakami and J. E. Luco (Editors) @ASCE, Reston, VA,
More informationRandom Dynamical Systems with Non-Gaussian Fluctuations Some Advances & Perspectives
Random Dynamical Systems with Non-Gaussian Fluctuations Some Advances & Perspectives Jinqiao Duan Illinois Institute of Technology Department of Applied Mathematics Chicago, IL 60616, USA www.iit.edu/~duan
More informationDynamical systems with Gaussian and Levy noise: analytical and stochastic approaches
Dynamical systems with Gaussian and Levy noise: analytical and stochastic approaches Noise is often considered as some disturbing component of the system. In particular physical situations, noise becomes
More information2012 NCTS Workshop on Dynamical Systems
Barbara Gentz gentz@math.uni-bielefeld.de http://www.math.uni-bielefeld.de/ gentz 2012 NCTS Workshop on Dynamical Systems National Center for Theoretical Sciences, National Tsing-Hua University Hsinchu,
More informationFirst exit times for Lévy-driven diffusions with exponentially light jumps
First exit times for Lévy-driven diffusions with exponentially light jumps P. Imkeller, I. Pavlyukevich and T. Wetzel Institut für Mathematik Humboldt-Universität zu Berlin Rudower Chaussee 25 12489 Berlin
More informationStochastic resonance: non-robust and robust tuning notions
Stochastic resonance: non-robust and robust tuning notions S. Herrmann 1 P. Imkeller I. Pavlyukevich 3 14th February 003 Abstract We consider a dynamical system describing the diffusive motion of a particle
More informationPower Variation of α-stable Lévy-Processes
Power Variation of α-stable Lévy-Processes and application to paleoclimatic modelling Jan Gairing C. Hein, P. Imkeller Department of Mathematics Humboldt Universität zu Berlin Weather Derivatives and Risk
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationGARCH processes continuous counterparts (Part 2)
GARCH processes continuous counterparts (Part 2) Alexander Lindner Centre of Mathematical Sciences Technical University of Munich D 85747 Garching Germany lindner@ma.tum.de http://www-m1.ma.tum.de/m4/pers/lindner/
More informationMetastability for the Ginzburg Landau equation with space time white noise
Barbara Gentz gentz@math.uni-bielefeld.de http://www.math.uni-bielefeld.de/ gentz Metastability for the Ginzburg Landau equation with space time white noise Barbara Gentz University of Bielefeld, Germany
More informationOrnstein-Uhlenbeck processes for geophysical data analysis
Ornstein-Uhlenbeck processes for geophysical data analysis Semere Habtemicael Department of Mathematics North Dakota State University May 19, 2014 Outline 1 Introduction 2 Model 3 Characteristic function
More informationMulti-Factor Lévy Models I: Symmetric alpha-stable (SαS) Lévy Processes
Multi-Factor Lévy Models I: Symmetric alpha-stable (SαS) Lévy Processes Anatoliy Swishchuk Department of Mathematics and Statistics University of Calgary Calgary, Alberta, Canada Lunch at the Lab Talk
More information{σ x >t}p x. (σ x >t)=e at.
3.11. EXERCISES 121 3.11 Exercises Exercise 3.1 Consider the Ornstein Uhlenbeck process in example 3.1.7(B). Show that the defined process is a Markov process which converges in distribution to an N(0,σ
More informationStochastic differential equations in neuroscience
Stochastic differential equations in neuroscience Nils Berglund MAPMO, Orléans (CNRS, UMR 6628) http://www.univ-orleans.fr/mapmo/membres/berglund/ Barbara Gentz, Universität Bielefeld Damien Landon, MAPMO-Orléans
More informationStochastic Differential Equations.
Chapter 3 Stochastic Differential Equations. 3.1 Existence and Uniqueness. One of the ways of constructing a Diffusion process is to solve the stochastic differential equation dx(t) = σ(t, x(t)) dβ(t)
More informationLocal vs. Nonlocal Diffusions A Tale of Two Laplacians
Local vs. Nonlocal Diffusions A Tale of Two Laplacians Jinqiao Duan Dept of Applied Mathematics Illinois Institute of Technology Chicago duan@iit.edu Outline 1 Einstein & Wiener: The Local diffusion 2
More informationStochastic (intermittent) Spikes and Strong Noise Limit of SDEs.
Stochastic (intermittent) Spikes and Strong Noise Limit of SDEs. D. Bernard in collaboration with M. Bauer and (in part) A. Tilloy. IPAM-UCLA, Jan 2017. 1 Strong Noise Limit of (some) SDEs Stochastic differential
More informationRandom and Deterministic perturbations of dynamical systems. Leonid Koralov
Random and Deterministic perturbations of dynamical systems Leonid Koralov - M. Freidlin, L. Koralov Metastability for Nonlinear Random Perturbations of Dynamical Systems, Stochastic Processes and Applications
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 15. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationAn adaptive numerical scheme for fractional differential equations with explosions
An adaptive numerical scheme for fractional differential equations with explosions Johanna Garzón Departamento de Matemáticas, Universidad Nacional de Colombia Seminario de procesos estocásticos Jointly
More informationfor all f satisfying E[ f(x) ] <.
. Let (Ω, F, P ) be a probability space and D be a sub-σ-algebra of F. An (H, H)-valued random variable X is independent of D if and only if P ({X Γ} D) = P {X Γ}P (D) for all Γ H and D D. Prove that if
More informationThe concentration of a drug in blood. Exponential decay. Different realizations. Exponential decay with noise. dc(t) dt.
The concentration of a drug in blood Exponential decay C12 concentration 2 4 6 8 1 C12 concentration 2 4 6 8 1 dc(t) dt = µc(t) C(t) = C()e µt 2 4 6 8 1 12 time in minutes 2 4 6 8 1 12 time in minutes
More informationOn Optimal Stopping Problems with Power Function of Lévy Processes
On Optimal Stopping Problems with Power Function of Lévy Processes Budhi Arta Surya Department of Mathematics University of Utrecht 31 August 2006 This talk is based on the joint paper with A.E. Kyprianou:
More informationStationary distributions of non Gaussian Ornstein Uhlenbeck processes for beam halos
N Cufaro Petroni: CYCLOTRONS 2007 Giardini Naxos, 1 5 October, 2007 1 Stationary distributions of non Gaussian Ornstein Uhlenbeck processes for beam halos CYCLOTRONS 2007 Giardini Naxos, 1 5 October Nicola
More informationIntroduction to Machine Learning CMU-10701
Introduction to Machine Learning CMU-10701 Markov Chain Monte Carlo Methods Barnabás Póczos & Aarti Singh Contents Markov Chain Monte Carlo Methods Goal & Motivation Sampling Rejection Importance Markov
More informationLeast Squares Estimators for Stochastic Differential Equations Driven by Small Lévy Noises
Least Squares Estimators for Stochastic Differential Equations Driven by Small Lévy Noises Hongwei Long* Department of Mathematical Sciences, Florida Atlantic University, Boca Raton Florida 33431-991,
More informationNumerical Methods with Lévy Processes
Numerical Methods with Lévy Processes 1 Objective: i) Find models of asset returns, etc ii) Get numbers out of them. Why? VaR and risk management Valuing and hedging derivatives Why not? Usual assumption:
More informationFunctional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals
Functional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals Noèlia Viles Cuadros BCAM- Basque Center of Applied Mathematics with Prof. Enrico
More informationJump-type Levy Processes
Jump-type Levy Processes Ernst Eberlein Handbook of Financial Time Series Outline Table of contents Probabilistic Structure of Levy Processes Levy process Levy-Ito decomposition Jump part Probabilistic
More informationGARCH processes probabilistic properties (Part 1)
GARCH processes probabilistic properties (Part 1) Alexander Lindner Centre of Mathematical Sciences Technical University of Munich D 85747 Garching Germany lindner@ma.tum.de http://www-m1.ma.tum.de/m4/pers/lindner/
More informationEvgeny Spodarev WIAS, Berlin. Limit theorems for excursion sets of stationary random fields
Evgeny Spodarev 23.01.2013 WIAS, Berlin Limit theorems for excursion sets of stationary random fields page 2 LT for excursion sets of stationary random fields Overview 23.01.2013 Overview Motivation Excursion
More informationDynamics of the evolving Bolthausen-Sznitman coalescent. by Jason Schweinsberg University of California at San Diego.
Dynamics of the evolving Bolthausen-Sznitman coalescent by Jason Schweinsberg University of California at San Diego Outline of Talk 1. The Moran model and Kingman s coalescent 2. The evolving Kingman s
More informationHigher order weak approximations of stochastic differential equations with and without jumps
Higher order weak approximations of stochastic differential equations with and without jumps Hideyuki TANAKA Graduate School of Science and Engineering, Ritsumeikan University Rough Path Analysis and Related
More informationPSEUDO RESONANCE INDUCED QUASI-PERIODIC BEHAVIOR IN STOCHASTIC THRESHOLD DYNAMICS
January 4, 8:56 WSPC/INSTRUCTION FILE Ditlevsen Braunfinal Stochastics and Dynamics c World Scientific Publishing Company PSEUDO RESONANCE INDUCED QUASI-PERIODIC BEHAVIOR IN STOCHASTIC THRESHOLD DYNAMICS
More informationSelf-intersection local time for Gaussian processes
Self-intersection local Olga Izyumtseva, olaizyumtseva@yahoo.com (in collaboration with Andrey Dorogovtsev, adoro@imath.kiev.ua) Department of theory of random processes Institute of mathematics Ukrainian
More informationIntroduction to asymptotic techniques for stochastic systems with multiple time-scales
Introduction to asymptotic techniques for stochastic systems with multiple time-scales Eric Vanden-Eijnden Courant Institute Motivating examples Consider the ODE {Ẋ = Y 3 + sin(πt) + cos( 2πt) X() = x
More informationp 1 ( Y p dp) 1/p ( X p dp) 1 1 p
Doob s inequality Let X(t) be a right continuous submartingale with respect to F(t), t 1 P(sup s t X(s) λ) 1 λ {sup s t X(s) λ} X + (t)dp 2 For 1 < p
More informationAnomalous Lévy diffusion: From the flight of an albatross to optical lattices. Eric Lutz Abteilung für Quantenphysik, Universität Ulm
Anomalous Lévy diffusion: From the flight of an albatross to optical lattices Eric Lutz Abteilung für Quantenphysik, Universität Ulm Outline 1 Lévy distributions Broad distributions Central limit theorem
More informationSmall-time asymptotics of stopped Lévy bridges and simulation schemes with controlled bias
Small-time asymptotics of stopped Lévy bridges and simulation schemes with controlled bias José E. Figueroa-López 1 1 Department of Statistics Purdue University Seoul National University & Ajou University
More informationDefinition: Lévy Process. Lectures on Lévy Processes and Stochastic Calculus, Braunschweig, Lecture 2: Lévy Processes. Theorem
Definition: Lévy Process Lectures on Lévy Processes and Stochastic Calculus, Braunschweig, Lecture 2: Lévy Processes David Applebaum Probability and Statistics Department, University of Sheffield, UK July
More informationExercises. T 2T. e ita φ(t)dt.
Exercises. Set #. Construct an example of a sequence of probability measures P n on R which converge weakly to a probability measure P but so that the first moments m,n = xdp n do not converge to m = xdp.
More informationLévy Processes and Infinitely Divisible Measures in the Dual of afebruary Nuclear2017 Space 1 / 32
Lévy Processes and Infinitely Divisible Measures in the Dual of a Nuclear Space David Applebaum School of Mathematics and Statistics, University of Sheffield, UK Talk at "Workshop on Infinite Dimensional
More informationHomework #6 : final examination Due on March 22nd : individual work
Université de ennes Année 28-29 Master 2ème Mathématiques Modèles stochastiques continus ou à sauts Homework #6 : final examination Due on March 22nd : individual work Exercise Warm-up : behaviour of characteristic
More informationPerformance Evaluation of Generalized Polynomial Chaos
Performance Evaluation of Generalized Polynomial Chaos Dongbin Xiu, Didier Lucor, C.-H. Su, and George Em Karniadakis 1 Division of Applied Mathematics, Brown University, Providence, RI 02912, USA, gk@dam.brown.edu
More informationPoisson Jumps in Credit Risk Modeling: a Partial Integro-differential Equation Formulation
Poisson Jumps in Credit Risk Modeling: a Partial Integro-differential Equation Formulation Jingyi Zhu Department of Mathematics University of Utah zhu@math.utah.edu Collaborator: Marco Avellaneda (Courant
More informationBranching Processes II: Convergence of critical branching to Feller s CSB
Chapter 4 Branching Processes II: Convergence of critical branching to Feller s CSB Figure 4.1: Feller 4.1 Birth and Death Processes 4.1.1 Linear birth and death processes Branching processes can be studied
More informationConvex Optimization CMU-10725
Convex Optimization CMU-10725 Simulated Annealing Barnabás Póczos & Ryan Tibshirani Andrey Markov Markov Chains 2 Markov Chains Markov chain: Homogen Markov chain: 3 Markov Chains Assume that the state
More informationOn Stochastic Adaptive Control & its Applications. Bozenna Pasik-Duncan University of Kansas, USA
On Stochastic Adaptive Control & its Applications Bozenna Pasik-Duncan University of Kansas, USA ASEAS Workshop, AFOSR, 23-24 March, 2009 1. Motivation: Work in the 1970's 2. Adaptive Control of Continuous
More informationModel Reduction and Stochastic Resonance
Model eduction and Stochastic esonance P. Imkeller I. Pavlyukevich 4th June Abstract We provide a mathematical underpinning of the physically widely known phenomenon of stochastic resonance, i.e. the optimal
More information1. Stochastic Processes and filtrations
1. Stochastic Processes and 1. Stoch. pr., A stochastic process (X t ) t T is a collection of random variables on (Ω, F) with values in a measurable space (S, S), i.e., for all t, In our case X t : Ω S
More informationMetastability for interacting particles in double-well potentials and Allen Cahn SPDEs
Nils Berglund nils.berglund@univ-orleans.fr http://www.univ-orleans.fr/mapmo/membres/berglund/ SPA 2017 Contributed Session: Interacting particle systems Metastability for interacting particles in double-well
More informationStability of Stochastic Differential Equations
Lyapunov stability theory for ODEs s Stability of Stochastic Differential Equations Part 1: Introduction Department of Mathematics and Statistics University of Strathclyde Glasgow, G1 1XH December 2010
More informationSmall ball probabilities and metric entropy
Small ball probabilities and metric entropy Frank Aurzada, TU Berlin Sydney, February 2012 MCQMC Outline 1 Small ball probabilities vs. metric entropy 2 Connection to other questions 3 Recent results for
More informationChange of Measure formula and the Hellinger Distance of two Lévy Processes
Change of Measure formula and the Hellinger Distance of two Lévy Processes Erika Hausenblas University of Salzburg, Austria Change of Measure formula and the Hellinger Distance of two Lévy Processes p.1
More informationA Concise Course on Stochastic Partial Differential Equations
A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original
More informationMonte Carlo methods for sampling-based Stochastic Optimization
Monte Carlo methods for sampling-based Stochastic Optimization Gersende FORT LTCI CNRS & Telecom ParisTech Paris, France Joint works with B. Jourdain, T. Lelièvre, G. Stoltz from ENPC and E. Kuhn from
More informationCTRW Limits: Governing Equations and Fractal Dimensions
CTRW Limits: Governing Equations and Fractal Dimensions Erkan Nane DEPARTMENT OF MATHEMATICS AND STATISTICS AUBURN UNIVERSITY August 19-23, 2013 Joint work with Z-Q. Chen, M. D Ovidio, M.M. Meerschaert,
More informationMetastability in a class of parabolic SPDEs
Metastability in a class of parabolic SPDEs Nils Berglund MAPMO, Université d Orléans CNRS, UMR 6628 et Fédération Denis Poisson www.univ-orleans.fr/mapmo/membres/berglund Collaborators: Florent Barret,
More informationHardy-Stein identity and Square functions
Hardy-Stein identity and Square functions Daesung Kim (joint work with Rodrigo Bañuelos) Department of Mathematics Purdue University March 28, 217 Daesung Kim (Purdue) Hardy-Stein identity UIUC 217 1 /
More informationAnomalous diffusion in biology: fractional Brownian motion, Lévy flights
Anomalous diffusion in biology: fractional Brownian motion, Lévy flights Jan Korbel Faculty of Nuclear Sciences and Physical Engineering, CTU, Prague Minisymposium on fundamental aspects behind cell systems
More informationLower Tail Probabilities and Related Problems
Lower Tail Probabilities and Related Problems Qi-Man Shao National University of Singapore and University of Oregon qmshao@darkwing.uoregon.edu . Lower Tail Probabilities Let {X t, t T } be a real valued
More informationSpectral functions of subordinated Brownian motion
Spectral functions of subordinated Brownian motion M.A. Fahrenwaldt 12 1 Institut für Mathematische Stochastik Leibniz Universität Hannover, Germany 2 EBZ Business School, Bochum, Germany Berlin, 23 October
More informationSimulation of diffusion. processes with discontinuous coefficients. Antoine Lejay Projet TOSCA, INRIA Nancy Grand-Est, Institut Élie Cartan
Simulation of diffusion. processes with discontinuous coefficients Antoine Lejay Projet TOSCA, INRIA Nancy Grand-Est, Institut Élie Cartan From collaborations with Pierre Étoré and Miguel Martinez . Divergence
More informationRegular Variation and Extreme Events for Stochastic Processes
1 Regular Variation and Extreme Events for Stochastic Processes FILIP LINDSKOG Royal Institute of Technology, Stockholm 2005 based on joint work with Henrik Hult www.math.kth.se/ lindskog 2 Extremes for
More informationA Short Introduction to Diffusion Processes and Ito Calculus
A Short Introduction to Diffusion Processes and Ito Calculus Cédric Archambeau University College, London Center for Computational Statistics and Machine Learning c.archambeau@cs.ucl.ac.uk January 24,
More informationSome Properties of NSFDEs
Chenggui Yuan (Swansea University) Some Properties of NSFDEs 1 / 41 Some Properties of NSFDEs Chenggui Yuan Swansea University Chenggui Yuan (Swansea University) Some Properties of NSFDEs 2 / 41 Outline
More informationLangevin Methods. Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 10 D Mainz Germany
Langevin Methods Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 1 D 55128 Mainz Germany Motivation Original idea: Fast and slow degrees of freedom Example: Brownian motion Replace
More informationMathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( )
Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio (2014-2015) Etienne Tanré - Olivier Faugeras INRIA - Team Tosca November 26th, 2014 E. Tanré (INRIA - Team Tosca) Mathematical
More informationI forgot to mention last time: in the Ito formula for two standard processes, putting
I forgot to mention last time: in the Ito formula for two standard processes, putting dx t = a t dt + b t db t dy t = α t dt + β t db t, and taking f(x, y = xy, one has f x = y, f y = x, and f xx = f yy
More informationHandbook of Stochastic Methods
C. W. Gardiner Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sciences Third Edition With 30 Figures Springer Contents 1. A Historical Introduction 1 1.1 Motivation I 1.2 Some Historical
More informationUCSD ECE250 Handout #27 Prof. Young-Han Kim Friday, June 8, Practice Final Examination (Winter 2017)
UCSD ECE250 Handout #27 Prof. Young-Han Kim Friday, June 8, 208 Practice Final Examination (Winter 207) There are 6 problems, each problem with multiple parts. Your answer should be as clear and readable
More informationBrownian Motion. 1 Definition Brownian Motion Wiener measure... 3
Brownian Motion Contents 1 Definition 2 1.1 Brownian Motion................................. 2 1.2 Wiener measure.................................. 3 2 Construction 4 2.1 Gaussian process.................................
More informationPartial Differential Equations with Applications to Finance Seminar 1: Proving and applying Dynkin s formula
Partial Differential Equations with Applications to Finance Seminar 1: Proving and applying Dynkin s formula Group 4: Bertan Yilmaz, Richard Oti-Aboagye and Di Liu May, 15 Chapter 1 Proving Dynkin s formula
More informationThe multidimensional Ito Integral and the multidimensional Ito Formula. Eric Mu ller June 1, 2015 Seminar on Stochastic Geometry and its applications
The multidimensional Ito Integral and the multidimensional Ito Formula Eric Mu ller June 1, 215 Seminar on Stochastic Geometry and its applications page 2 Seminar on Stochastic Geometry and its applications
More informationHarnack Inequalities and Applications for Stochastic Equations
p. 1/32 Harnack Inequalities and Applications for Stochastic Equations PhD Thesis Defense Shun-Xiang Ouyang Under the Supervision of Prof. Michael Röckner & Prof. Feng-Yu Wang March 6, 29 p. 2/32 Outline
More informationErnesto Mordecki. Talk presented at the. Finnish Mathematical Society
EXACT RUIN PROBABILITIES FOR A CLASS Of LÉVY PROCESSES Ernesto Mordecki http://www.cmat.edu.uy/ mordecki Montevideo, Uruguay visiting Åbo Akademi, Turku Talk presented at the Finnish Mathematical Society
More informationSimulation methods for stochastic models in chemistry
Simulation methods for stochastic models in chemistry David F. Anderson anderson@math.wisc.edu Department of Mathematics University of Wisconsin - Madison SIAM: Barcelona June 4th, 21 Overview 1. Notation
More informationCOMPUTER-AIDED MODELING AND SIMULATION OF ELECTRICAL CIRCUITS WITH α-stable NOISE
APPLICATIONES MATHEMATICAE 23,1(1995), pp. 83 93 A. WERON(Wroc law) COMPUTER-AIDED MODELING AND SIMULATION OF ELECTRICAL CIRCUITS WITH α-stable NOISE Abstract. The aim of this paper is to demonstrate how
More informationMoment Properties of Distributions Used in Stochastic Financial Models
Moment Properties of Distributions Used in Stochastic Financial Models Jordan Stoyanov Newcastle University (UK) & University of Ljubljana (Slovenia) e-mail: stoyanovj@gmail.com ETH Zürich, Department
More informationHypoelliptic multiscale Langevin diffusions and Slow fast stochastic reaction diffusion equations.
Hypoelliptic multiscale Langevin diffusions and Slow fast stochastic reaction diffusion equations. Wenqing Hu. 1 (Joint works with Michael Salins 2 and Konstantinos Spiliopoulos 3.) 1. Department of Mathematics
More informationA random perturbation approach to some stochastic approximation algorithms in optimization.
A random perturbation approach to some stochastic approximation algorithms in optimization. Wenqing Hu. 1 (Presentation based on joint works with Chris Junchi Li 2, Weijie Su 3, Haoyi Xiong 4.) 1. Department
More informationContinuous dependence estimates for the ergodic problem with an application to homogenization
Continuous dependence estimates for the ergodic problem with an application to homogenization Claudio Marchi Bayreuth, September 12 th, 2013 C. Marchi (Università di Padova) Continuous dependence Bayreuth,
More informationBrownian intersection local times
Weierstrass Institute for Applied Analysis and Stochastics Brownian intersection local times Wolfgang König (TU Berlin and WIAS Berlin) [K.&M., Ann. Probab. 2002], [K.&M., Trans. Amer. Math. Soc. 2006]
More informationPotential Theory on Wiener space revisited
Potential Theory on Wiener space revisited Michael Röckner (University of Bielefeld) Joint work with Aurel Cornea 1 and Lucian Beznea (Rumanian Academy, Bukarest) CRC 701 and BiBoS-Preprint 1 Aurel tragically
More informationA Lévy-Fokker-Planck equation: entropies and convergence to equilibrium
1/ 22 A Lévy-Fokker-Planck equation: entropies and convergence to equilibrium I. Gentil CEREMADE, Université Paris-Dauphine International Conference on stochastic Analysis and Applications Hammamet, Tunisia,
More informationExponential functionals of Lévy processes
Exponential functionals of Lévy processes Víctor Rivero Centro de Investigación en Matemáticas, México. 1/ 28 Outline of the talk Introduction Exponential functionals of spectrally positive Lévy processes
More informationDensities for the Navier Stokes equations with noise
Densities for the Navier Stokes equations with noise Marco Romito Università di Pisa Universitat de Barcelona March 25, 2015 Summary 1 Introduction & motivations 2 Malliavin calculus 3 Besov bounds 4 Other
More informationOperational Risk and Pareto Lévy Copulas
Operational Risk and Pareto Lévy Copulas Claudia Klüppelberg Technische Universität München email: cklu@ma.tum.de http://www-m4.ma.tum.de References: - Böcker, K. and Klüppelberg, C. (25) Operational VaR
More informationLocal time path integrals and their application to Lévy random walks
Local time path integrals and their application to Lévy random walks Václav Zatloukal (www.zatlovac.eu) Faculty of Nuclear Sciences and Physical Engineering Czech Technical University in Prague talk given
More informationSelf-similar Markov processes
Self-similar Markov processes Andreas E. Kyprianou 1 1 Unversity of Bath Which are our favourite stochastic processes? Markov chains Diffusions Brownian motion Cts-time Markov processes with jumps Lévy
More informationON ADDITIVE TIME-CHANGES OF FELLER PROCESSES. 1. Introduction
ON ADDITIVE TIME-CHANGES OF FELLER PROCESSES ALEKSANDAR MIJATOVIĆ AND MARTIJN PISTORIUS Abstract. In this note we generalise the Phillips theorem [1] on the subordination of Feller processes by Lévy subordinators
More informationStochastic Analysis. Prof. Dr. Andreas Eberle
Stochastic Analysis Prof. Dr. Andreas Eberle March 13, 212 Contents Contents 2 1 Lévy processes and Poisson point processes 6 1.1 Lévy processes.............................. 7 Characteristic exponents.........................
More informationELEMENTS OF PROBABILITY THEORY
ELEMENTS OF PROBABILITY THEORY Elements of Probability Theory A collection of subsets of a set Ω is called a σ algebra if it contains Ω and is closed under the operations of taking complements and countable
More informationLatent voter model on random regular graphs
Latent voter model on random regular graphs Shirshendu Chatterjee Cornell University (visiting Duke U.) Work in progress with Rick Durrett April 25, 2011 Outline Definition of voter model and duality with
More informationModelling Interactions Between Weather and Climate
Modelling Interactions Between Weather and Climate p. 1/33 Modelling Interactions Between Weather and Climate Adam Monahan & Joel Culina monahana@uvic.ca & culinaj@uvic.ca School of Earth and Ocean Sciences,
More informationHomogenization with stochastic differential equations
Homogenization with stochastic differential equations Scott Hottovy shottovy@math.arizona.edu University of Arizona Program in Applied Mathematics October 12, 2011 Modeling with SDE Use SDE to model system
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationMemory and hypoellipticity in neuronal models
Memory and hypoellipticity in neuronal models S. Ditlevsen R. Höpfner E. Löcherbach M. Thieullen Banff, 2017 What this talk is about : What is the effect of memory in probabilistic models for neurons?
More information