Activated Dynamics of the Cascading 2-GREM Scaling limits & Aging

Size: px
Start display at page:

Download "Activated Dynamics of the Cascading 2-GREM Scaling limits & Aging"

Transcription

1 Activated Dynamics of the Cascading 2-GREM Scaling limits & Aging Véronique Gayrard CNRS, Université d Aix-Marseille Based on joint work with L.R. Fontes, arxiv: Spin Glasses and Related Topics, Banff September 30-October 5, 2018

2 V n = { 1,1} n. 2-level GREM Hamiltonian (H n (σ),σ V n ) centered gaussian process on (Ω,F,P) with where and E(H n (σ)h n (η)) = na(1 ρ n (σ,η)), ρ n (σ,η) = 1 n 1 (min{i σ(i) η(i)} 1) A(x) 1 a 0 p 1 x Alternatively, setting n 1 = pn, n = n 1 +n 2, σ = σ 1 σ 2 V n1 V n2, H n (σ) = H n (1) (σ 1 )+H n (2) (σ 1 σ 2 ) = ang σ (1) 1 + (1 a)ng σ (2) 1 σ 2, where {g (1) σ 1,g (2) σ 1 σ 2 ; σ V n } are i.i.d. N(0,1).

3 STATIC Ground states (extremes of H n ): H n (1) (σ1 1 ) H(1) n (σ1 2 ) H(1) n (σ x 1 1 )... and for each x 1 {1,...,2 n 1} H n (2) (σ x 1 1 σx 11 2 ) H(2) n (σ x 1 1 σx 12 2 ) H(2) n (σ x 1 2 )... 1 σx 1x 2 a > p a < p Free energy for a > p exhibits two discontinuities at the critical inverse temperatures β 1 = β p a < β = 2ln2 < β 2 = β 1 p 1 a

4 Which dynamics to choose? Glauber dynamics: Markov jump process (X n (t),t > 0) on V n with rates e βh n(σ) λn (σ,η) = e βh n(η) λn (η,σ) σ η, and λ n (σ,η) = 0 else. Random Hopping: Metropolis: λ RH n (σ,η) e βh n(σ), σ η. Because H n is hierarchical, i.e. λ Metro n (σ,η) e β[h n(η) H n (σ)] +, σ η. H n (σ) = H (1) n (σ 1 )+H (2) n (σ 1 σ 2 ), σ = σ 1 σ 2 V n1 V n2, so is Metropolis λ Metro n (σ,η) e β[h n(η) H n (σ)] +½σ 1 +e β[h(2) n (η) H n (2) (σ)] +½σ 2. η η We choose λ n (σ,η) e βh n(σ) ½σ 1 η +e βh(2) n (σ) ½σ 2 η.

5 AGING... t 0 t w t w +t time unit = c n µ n : initial distribution c n : time scale C n (t w,t), t w,t > 0: time-time correlation function, e.g., C n (t w,t) = P µn [n 1( ) X n (c n t w ),X n (c n (t w +t)) 1 δ ], 0 < δ < 1...occurs if C (ρ) 1 or lim n C n(t w,t) = C (t/t w ), lim tw,t 0 t/tw=ρ lim n C n(t w,t) = C (ρ). Asymptotically, C n (t w,t) is not time translation invariant ρ

6 Random hopping dynamics of the REM static critical temperature: β = 2log2. time-scale: given β eff > 0, c n (β eff ) = exp { nββ eff β 2β eff ( log ( β 2 eff n/2 ) +log4π )}. start = uniform distribution on V n if β eff β, any distribution else. Stationary (II) Ergodic (IV) Stationary (I) Aging (III)

7 Metropolis dynamics of the REM Stationary (II) Stationary (I) Aging (III) Phase transition lines are not understood: aging is proved for that is, on time scales β > β eff, β eff < β (1 cst. ) logn n, (V.G., PTRF 18) c n (β eff ) < e cst nlogn c n (β ).

8 K-processes Markov processes in random environment on countably infinite state space, N { }, with a single unstable state { }. Simplest example: the uniform K-process (as in the RHD of the REM). when in a stable state x, waits an exponential time of mean value γ(x), x N γ(x) < then jumps to the unstable state, & from there, enters any finite set of stable states with uniform distribution. More general re-entrance mechanisms (not uniform but weighted) or state spaces (hierarchical) needed to describe more involved models or dynamics. Ergodic (infinite volume) dynamics of mean-field spin glasses are modeled in physics through so-called (finite volume) trap models. Two GREM-like trap models were proposed (Bouchaud & Dean, Sasaki & Nemoto). Is one of them the right one?

9 Back to the cascading 2-GREM (a > p) Set λ n (σ,η) e β(h(1) n (σ 1 )+H n (2) (σ 1 σ 2 )) ½σ 1 +e βh(2) n (σ) ½σ 2 η η c i n max σ i V ni e βh(i) n (σ) e nββ ap..., i = 1, e nββ (1 a)(1 p)..., i = 2. Two sources of dynamical phase transitions on extreme scales Static: β 1 = β p a β = 2ln2 β 2 = β 1 p 1 a Fine tuning: Given ζ n < b n n 2β 2 2 let β n (a,p,ζ n ) be the solution in β of In explicit form: We say that β 1 = β 1 n (a,p,ζ n ) is 2 n 2/c 1 n = e ζn. β n (a,p,ζ n ) β FT( 1 ζ n bn ), β FT β 2 1 p pa. above at below fine tuning temperature if ζ n + ζ n ζ ζ n.

10 Assume that a > p,β > β 2 and β FT > β 2. β 1 β β 2 β FT Theorem (Scaling limits).the rescaled process X n (c n t) has three possible scaling limits, Y = Y 1 Y 2, above (c n = c 2 n), below (c n = 2 n 2c 1 nc 2 n) and at fine tuning (c n = e ζ c 2 n), all of them ergodic processes expressed through K-processes. Convergence is in law w.r.t. the environment. Choose as time-time correlation function C(t w,t) = 1P(A 1 A 2 )+pp(a 1 A c 2 )+(1 p)p(ac 1 A 2) where for i = 1,2, A i = A i (t w,t) = {Y i does not jump between times t w and t w +t}. Theorem (Aging). ( ) Asl 1 α2 1+ρ, above FT, ( ) ( ) lim C(t w,t) = pasl 1 α1 α tw,t ρ +(1 p)aslα2 1 1+ρ, at FT, ( ) t/tw ρ pasl 1 α1, below FT, 1+ρ where α i = β i /β (0,1) and Asl is the p.d.f. of the arcsine law.

11 γn(x 1 1 ) = e βh(1) n (σ x 1 1 ) /c 1 n x 1, γn(x 2 1 x 2 ) = e βh(2) n (σ x 1 1 σx 2 2 ) /c 2 n Set S = S 1 S 2, S i = {1,...,M i }, i = 1,2. For x = x 1 x 2,y = y 1 y 2 S, [ x (1 λ 1 n )+λ x 1 n νn(x 1 1 ) ] 1 Proba(x y) = (1+o(1)) M2, if x 1 = y 1, λ x 1 n νn 1(y 1) M 1, else, 2 where λ x 1 n λ x 1 = 1, above FT, 1 1+e ζ M 2 γ 1 (x 1 ) and ν 1 n is the random measure on S 1 ν 1 n (x 1), at FT, 0, below FT, γ1 (x 1 ) z S γ 1 (z), 1 γ1 (x 1 )λ x 1 z S γ 1 (z)λ 1 1 M, 1 above FT, z, at FT, below FT.

12 None of the GREM-like trap models of theoretical physics correctly predicted the dynamics of the GREM at extreme (i.e. ergodic) time-scales.

13 THANK YOU FOR YOUR ATTENTION!

Lecture 20: Reversible Processes and Queues

Lecture 20: Reversible Processes and Queues Lecture 20: Reversible Processes and Queues 1 Examples of reversible processes 11 Birth-death processes We define two non-negative sequences birth and death rates denoted by {λ n : n N 0 } and {µ n : n

More information

Where Probability Meets Combinatorics and Statistical Mechanics

Where Probability Meets Combinatorics and Statistical Mechanics Where Probability Meets Combinatorics and Statistical Mechanics Daniel Ueltschi Department of Mathematics, University of Warwick MASDOC, 15 March 2011 Collaboration with V. Betz, N.M. Ercolani, C. Goldschmidt,

More information

{σ x >t}p x. (σ x >t)=e at.

{σ 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 information

Lecture 21: Convergence of transformations and generating a random variable

Lecture 21: Convergence of transformations and generating a random variable Lecture 21: Convergence of transformations and generating a random variable If Z n converges to Z in some sense, we often need to check whether h(z n ) converges to h(z ) in the same sense. Continuous

More information

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS

PROBABILITY: 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 information

Uniqueness of random gradient states

Uniqueness of random gradient states p. 1 Uniqueness of random gradient states Codina Cotar (Based on joint work with Christof Külske) p. 2 Model Interface transition region that separates different phases p. 2 Model Interface transition

More information

Metastability for the Ising model on a random graph

Metastability for the Ising model on a random graph Metastability for the Ising model on a random graph Joint project with Sander Dommers, Frank den Hollander, Francesca Nardi Outline A little about metastability Some key questions A little about the, and

More information

III. Quantum ergodicity on graphs, perspectives

III. Quantum ergodicity on graphs, perspectives III. Quantum ergodicity on graphs, perspectives Nalini Anantharaman Université de Strasbourg 24 août 2016 Yesterday we focussed on the case of large regular (discrete) graphs. Let G = (V, E) be a (q +

More information

Gradient Gibbs measures with disorder

Gradient Gibbs measures with disorder Gradient Gibbs measures with disorder Codina Cotar University College London April 16, 2015, Providence Partly based on joint works with Christof Külske Outline 1 The model 2 Questions 3 Known results

More information

A. Bovier () Branching Brownian motion: extremal process and ergodic theorems

A. Bovier () Branching Brownian motion: extremal process and ergodic theorems Branching Brownian motion: extremal process and ergodic theorems Anton Bovier with Louis-Pierre Arguin and Nicola Kistler RCS&SM, Venezia, 06.05.2013 Plan 1 BBM 2 Maximum of BBM 3 The Lalley-Sellke conjecture

More information

Random dynamical systems with microstructure

Random dynamical systems with microstructure Random dynamical systems with microstructure Renato Feres Washington University, St. Louis ESI, July 2011 1/37 Acknowledgements Different aspects of this work are joint with: Tim Chumley (Central limit

More information

Zig-Zag Monte Carlo. Delft University of Technology. Joris Bierkens February 7, 2017

Zig-Zag Monte Carlo. Delft University of Technology. Joris Bierkens February 7, 2017 Zig-Zag Monte Carlo Delft University of Technology Joris Bierkens February 7, 2017 Joris Bierkens (TU Delft) Zig-Zag Monte Carlo February 7, 2017 1 / 33 Acknowledgements Collaborators Andrew Duncan Paul

More information

Aging, rejuvenation and memory: the example of spin glasses

Aging, rejuvenation and memory: the example of spin glasses Aging, rejuvenation and memory: the example of spin glasses F. Bert, J.-P. Bouchaud, V. Dupuis, J. Hammann, D. Hérisson, F. Ladieu, M. Ocio, D. Parker, E. Vincent Service de Physique de l Etat Condensé,

More information

ELEC546 Review of Information Theory

ELEC546 Review of Information Theory ELEC546 Review of Information Theory Vincent Lau 1/1/004 1 Review of Information Theory Entropy: Measure of uncertainty of a random variable X. The entropy of X, H(X), is given by: If X is a discrete random

More information

1 Continuous-time chains, finite state space

1 Continuous-time chains, finite state space Université Paris Diderot 208 Markov chains Exercises 3 Continuous-time chains, finite state space Exercise Consider a continuous-time taking values in {, 2, 3}, with generator 2 2. 2 2 0. Draw the diagramm

More information

Phase transitions and finite-size scaling

Phase transitions and finite-size scaling Phase transitions and finite-size scaling Critical slowing down and cluster methods. Theory of phase transitions/ RNG Finite-size scaling Detailed treatment: Lectures on Phase Transitions and the Renormalization

More information

Chapter 2: Fundamentals of Statistics Lecture 15: Models and statistics

Chapter 2: Fundamentals of Statistics Lecture 15: Models and statistics Chapter 2: Fundamentals of Statistics Lecture 15: Models and statistics Data from one or a series of random experiments are collected. Planning experiments and collecting data (not discussed here). Analysis:

More information

7.1 Coupling from the Past

7.1 Coupling from the Past Georgia Tech Fall 2006 Markov Chain Monte Carlo Methods Lecture 7: September 12, 2006 Coupling from the Past Eric Vigoda 7.1 Coupling from the Past 7.1.1 Introduction We saw in the last lecture how Markov

More information

Monte Carlo methods for sampling-based Stochastic Optimization

Monte 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 information

A mathematical model for Mott s Variable Range Hopping

A mathematical model for Mott s Variable Range Hopping Santiago November 23-28, 2009 1 A mathematical model for Mott s Variable Range Hopping Jean BELLISSARD 1 2 Georgia Institute of Technology, Atlanta, http://www.math.gatech.edu/ jeanbel/ Collaborations:

More information

STA205 Probability: Week 8 R. Wolpert

STA205 Probability: Week 8 R. Wolpert INFINITE COIN-TOSS AND THE LAWS OF LARGE NUMBERS The traditional interpretation of the probability of an event E is its asymptotic frequency: the limit as n of the fraction of n repeated, similar, and

More information

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS

PROBABILITY: 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 information

6 Markov Chain Monte Carlo (MCMC)

6 Markov Chain Monte Carlo (MCMC) 6 Markov Chain Monte Carlo (MCMC) The underlying idea in MCMC is to replace the iid samples of basic MC methods, with dependent samples from an ergodic Markov chain, whose limiting (stationary) distribution

More information

Lattice Gaussian Sampling with Markov Chain Monte Carlo (MCMC)

Lattice Gaussian Sampling with Markov Chain Monte Carlo (MCMC) Lattice Gaussian Sampling with Markov Chain Monte Carlo (MCMC) Cong Ling Imperial College London aaa joint work with Zheng Wang (Huawei Technologies Shanghai) aaa September 20, 2016 Cong Ling (ICL) MCMC

More information

Markov processes and queueing networks

Markov processes and queueing networks Inria September 22, 2015 Outline Poisson processes Markov jump processes Some queueing networks The Poisson distribution (Siméon-Denis Poisson, 1781-1840) { } e λ λ n n! As prevalent as Gaussian distribution

More information

Gärtner-Ellis Theorem and applications.

Gärtner-Ellis Theorem and applications. Gärtner-Ellis Theorem and applications. Elena Kosygina July 25, 208 In this lecture we turn to the non-i.i.d. case and discuss Gärtner-Ellis theorem. As an application, we study Curie-Weiss model with

More information

Stochastic domination in space-time for the supercritical contact process

Stochastic domination in space-time for the supercritical contact process Stochastic domination in space-time for the supercritical contact process Stein Andreas Bethuelsen joint work with Rob van den Berg (CWI and VU Amsterdam) Workshop on Genealogies of Interacting Particle

More information

Scaling limit of the two-dimensional dynamic Ising-Kac model

Scaling limit of the two-dimensional dynamic Ising-Kac model Scaling limit of the two-dimensional dynamic Ising-Kac model Jean-Christophe Mourrat Hendrik Weber Mathematics Institute University of Warwick Statistical Mechanics Seminar, Warwick A famous result Theorem

More information

Practical unbiased Monte Carlo for Uncertainty Quantification

Practical unbiased Monte Carlo for Uncertainty Quantification Practical unbiased Monte Carlo for Uncertainty Quantification Sergios Agapiou Department of Statistics, University of Warwick MiR@W day: Uncertainty in Complex Computer Models, 2nd February 2015, University

More information

Stationary distributions of non Gaussian Ornstein Uhlenbeck processes for beam halos

Stationary 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 information

Weak Ergodicity Breaking WCHAOS 2011

Weak Ergodicity Breaking WCHAOS 2011 Weak Ergodicity Breaking Eli Barkai Bar-Ilan University Bel, Burov, Korabel, Margolin, Rebenshtok WCHAOS 211 Outline Single molecule experiments exhibit weak ergodicity breaking. Blinking quantum dots,

More information

SOLVABLE VARIATIONAL PROBLEMS IN N STATISTICAL MECHANICS

SOLVABLE VARIATIONAL PROBLEMS IN N STATISTICAL MECHANICS SOLVABLE VARIATIONAL PROBLEMS IN NON EQUILIBRIUM STATISTICAL MECHANICS University of L Aquila October 2013 Tullio Levi Civita Lecture 2013 Coauthors Lorenzo Bertini Alberto De Sole Alessandra Faggionato

More information

Hopping transport in disordered solids

Hopping transport in disordered solids Hopping transport in disordered solids Dominique Spehner Institut Fourier, Grenoble, France Workshop on Quantum Transport, Lexington, March 17, 2006 p. 1 Outline of the talk Hopping transport Models for

More information

The Ising model Summary of L12

The Ising model Summary of L12 The Ising model Summary of L2 Aim: Study connections between macroscopic phenomena and the underlying microscopic world for a ferromagnet. How: Study the simplest possible model of a ferromagnet containing

More information

Spectra of Large Random Stochastic Matrices & Relaxation in Complex Systems

Spectra of Large Random Stochastic Matrices & Relaxation in Complex Systems Spectra of Large Random Stochastic Matrices & Relaxation in Complex Systems Reimer Kühn Disordered Systems Group Department of Mathematics, King s College London Random Graphs and Random Processes, KCL

More information

Propagation of Chaos and Poisson Hypothesis

Propagation of Chaos and Poisson Hypothesis arxiv:1610.0849v1 [math.pr] 6 Oct 016 Propagation of Chaos and Poisson Hypothesis Serge Pirogov, Alexander Rybko, Senya Shlosman, and Alexander Vladimirov Aix Marseille Univ, Université de Toulon, CNRS,

More information

Asymptotic series in quantum mechanics: anharmonic oscillator

Asymptotic series in quantum mechanics: anharmonic oscillator Asymptotic series in quantum mechanics: anharmonic oscillator Facultat de Física, Universitat de Barcelona, Diagonal 645, 0808 Barcelona, Spain Thesis supervisors: Dr Bartomeu Fiol Núñez, Dr Alejandro

More information

A STUDY ON CONDENSATION IN ZERO RANGE PROCESSES

A STUDY ON CONDENSATION IN ZERO RANGE PROCESSES J. KSIAM Vol.22, No.3, 37 54, 208 http://dx.doi.org/0.294/jksiam.208.22.37 A STUDY ON CONDENSATION IN ZERO RANGE PROCESSES CHEOL-UNG PARK AND INTAE JEON 2 XINAPSE, SOUTH KOREA E-mail address: cheol.ung.park0@gmail.com

More information

Ferromagnets and the classical Heisenberg model. Kay Kirkpatrick, UIUC

Ferromagnets and the classical Heisenberg model. Kay Kirkpatrick, UIUC Ferromagnets and the classical Heisenberg model Kay Kirkpatrick, UIUC Ferromagnets and the classical Heisenberg model: asymptotics for a mean-field phase transition Kay Kirkpatrick, Urbana-Champaign June

More information

Domino shuffling and TASEP. Xufan Zhang. Brown Discrete Math Seminar, March 2018

Domino shuffling and TASEP. Xufan Zhang. Brown Discrete Math Seminar, March 2018 Domino shuffling and TASEP Xufan Zhang Brown Discrete Math Seminar, March 2018 Table of Contents 1 Domino tilings in an Aztec diamond 2 Domino shuffling 3 Totally asymmetric simple exclusion process(tasep)

More information

Lecture 6: Markov Chain Monte Carlo

Lecture 6: Markov Chain Monte Carlo Lecture 6: Markov Chain Monte Carlo D. Jason Koskinen koskinen@nbi.ku.dk Photo by Howard Jackman University of Copenhagen Advanced Methods in Applied Statistics Feb - Apr 2016 Niels Bohr Institute 2 Outline

More information

Weak quenched limiting distributions of a one-dimensional random walk in a random environment

Weak quenched limiting distributions of a one-dimensional random walk in a random environment Weak quenched limiting distributions of a one-dimensional random walk in a random environment Jonathon Peterson Cornell University Department of Mathematics Joint work with Gennady Samorodnitsky September

More information

Strong approximation for additive functionals of geometrically ergodic Markov chains

Strong approximation for additive functionals of geometrically ergodic Markov chains Strong approximation for additive functionals of geometrically ergodic Markov chains Florence Merlevède Joint work with E. Rio Université Paris-Est-Marne-La-Vallée (UPEM) Cincinnati Symposium on Probability

More information

STAT STOCHASTIC PROCESSES. Contents

STAT STOCHASTIC PROCESSES. Contents STAT 3911 - STOCHASTIC PROCESSES ANDREW TULLOCH Contents 1. Stochastic Processes 2 2. Classification of states 2 3. Limit theorems for Markov chains 4 4. First step analysis 5 5. Branching processes 5

More information

Bayesian data analysis in practice: Three simple examples

Bayesian data analysis in practice: Three simple examples Bayesian data analysis in practice: Three simple examples Martin P. Tingley Introduction These notes cover three examples I presented at Climatea on 5 October 0. Matlab code is available by request to

More information

RENORMALIZATION OF DYSON S VECTOR-VALUED HIERARCHICAL MODEL AT LOW TEMPERATURES

RENORMALIZATION OF DYSON S VECTOR-VALUED HIERARCHICAL MODEL AT LOW TEMPERATURES RENORMALIZATION OF DYSON S VECTOR-VALUED HIERARCHICAL MODEL AT LOW TEMPERATURES P. M. Bleher (1) and P. Major (2) (1) Keldysh Institute of Applied Mathematics of the Soviet Academy of Sciences Moscow (2)

More information

Cluster Expansion in the Canonical Ensemble

Cluster Expansion in the Canonical Ensemble Università di Roma Tre September 2, 2011 (joint work with Dimitrios Tsagkarogiannis) Background Motivation: phase transitions for a system of particles in the continuum interacting via a Kac potential,

More information

ECE276A: Sensing & Estimation in Robotics Lecture 10: Gaussian Mixture and Particle Filtering

ECE276A: Sensing & Estimation in Robotics Lecture 10: Gaussian Mixture and Particle Filtering ECE276A: Sensing & Estimation in Robotics Lecture 10: Gaussian Mixture and Particle Filtering Lecturer: Nikolay Atanasov: natanasov@ucsd.edu Teaching Assistants: Siwei Guo: s9guo@eng.ucsd.edu Anwesan Pal:

More information

Estimation of arrival and service rates for M/M/c queue system

Estimation of arrival and service rates for M/M/c queue system Estimation of arrival and service rates for M/M/c queue system Katarína Starinská starinskak@gmail.com Charles University Faculty of Mathematics and Physics Department of Probability and Mathematical Statistics

More information

Nonparametric regression with martingale increment errors

Nonparametric regression with martingale increment errors S. Gaïffas (LSTA - Paris 6) joint work with S. Delattre (LPMA - Paris 7) work in progress Motivations Some facts: Theoretical study of statistical algorithms requires stationary and ergodicity. Concentration

More information

In Memory of Wenbo V Li s Contributions

In Memory of Wenbo V Li s Contributions In Memory of Wenbo V Li s Contributions Qi-Man Shao The Chinese University of Hong Kong qmshao@cuhk.edu.hk The research is partially supported by Hong Kong RGC GRF 403513 Outline Lower tail probabilities

More information

Minimum average value-at-risk for finite horizon semi-markov decision processes

Minimum average value-at-risk for finite horizon semi-markov decision processes 12th workshop on Markov processes and related topics Minimum average value-at-risk for finite horizon semi-markov decision processes Xianping Guo (with Y.H. HUANG) Sun Yat-Sen University Email: mcsgxp@mail.sysu.edu.cn

More information

Lecture Characterization of Infinitely Divisible Distributions

Lecture Characterization of Infinitely Divisible Distributions Lecture 10 1 Characterization of Infinitely Divisible Distributions We have shown that a distribution µ is infinitely divisible if and only if it is the weak limit of S n := X n,1 + + X n,n for a uniformly

More information

Mixing Times and Hitting Times

Mixing Times and Hitting Times Mixing Times and Hitting Times David Aldous January 12, 2010 Old Results and Old Puzzles Levin-Peres-Wilmer give some history of the emergence of this Mixing Times topic in the early 1980s, and we nowadays

More information

Estimation theory. Parametric estimation. Properties of estimators. Minimum variance estimator. Cramer-Rao bound. Maximum likelihood estimators

Estimation theory. Parametric estimation. Properties of estimators. Minimum variance estimator. Cramer-Rao bound. Maximum likelihood estimators Estimation theory Parametric estimation Properties of estimators Minimum variance estimator Cramer-Rao bound Maximum likelihood estimators Confidence intervals Bayesian estimation 1 Random Variables Let

More information

Introduction to Markov Chains, Queuing Theory, and Network Performance

Introduction to Markov Chains, Queuing Theory, and Network Performance Introduction to Markov Chains, Queuing Theory, and Network Performance Marceau Coupechoux Telecom ParisTech, departement Informatique et Réseaux marceau.coupechoux@telecom-paristech.fr IT.2403 Modélisation

More information

Large sample distribution for fully functional periodicity tests

Large sample distribution for fully functional periodicity tests Large sample distribution for fully functional periodicity tests Siegfried Hörmann Institute for Statistics Graz University of Technology Based on joint work with Piotr Kokoszka (Colorado State) and Gilles

More information

Gradient interfaces with and without disorder

Gradient interfaces with and without disorder Gradient interfaces with and without disorder Codina Cotar University College London September 09, 2014, Toronto Outline 1 Physics motivation Example 1: Elasticity Recap-Gaussian Measure Example 2: Effective

More information

Weak Ergodicity Breaking

Weak Ergodicity Breaking Weak Ergodicity Breaking Eli Barkai Bar-Ilan University 215 Ergodicity Ergodicity: time averages = ensemble averages. x = lim t t x(τ)dτ t. x = xp eq (x)dx. Ergodicity out of equilibrium δ 2 (, t) = t

More information

NUMERICAL METHODS FOR QUANTUM IMPURITY MODELS

NUMERICAL METHODS FOR QUANTUM IMPURITY MODELS NUMERICAL METODS FOR QUANTUM IMPURITY MODELS http://www.staff.science.uu.nl/~mitch003/nrg.html March 2015 Andrew Mitchell, Utrecht University Quantum impurity problems Part 1: Quantum impurity problems

More information

Variational principle for generalized Gibbs measures

Variational principle for generalized Gibbs measures Variational principle for generalized Gibbs measures Külske, C.; Le Ny, A.; Redig, F.H.J. Published: 0/0/2002 Document Version Publisher s PDF, also known as Version of Record (includes final page, issue

More information

Biased random walks on subcritical percolation clusters with an infinite open path

Biased random walks on subcritical percolation clusters with an infinite open path Biased random walks on subcritical percolation clusters with an infinite open path Application for Transfer to DPhil Student Status Dissertation Annika Heckel March 3, 202 Abstract We consider subcritical

More information

Topic 2 Limits and Continuity c and d) Continuity Handout Notes Assigned Problems: Intro book pg 73, 1-3 and 6-8

Topic 2 Limits and Continuity c and d) Continuity Handout Notes Assigned Problems: Intro book pg 73, 1-3 and 6-8 c&d. Continuity Handout. Page 1 of 5 Topic Limits and Continuity c and d) Continuity Handout Notes Assigned Problems: Intro book pg 73, 1-3 and 6-8 Recall Limits and Function Values: We have already studied

More information

AGING OF THE METROPOLIS DYNAMICS ON THE RANDOM ENERGY MODEL

AGING OF THE METROPOLIS DYNAMICS ON THE RANDOM ENERGY MODEL AGIG OF THE METROPOLIS DYAMICS O THE RADOM EERGY MODEL JIŘÍ ČERÝ AD TOBIAS WASSMER Abstract. We study the Metropolis dynamics of the simplest mean-field spin glass model, the Random Energy Model. We show

More information

(Nearly) Scale invariant fluid dynamics for the dilute Fermi gas in two and three dimensions. Thomas Schaefer North Carolina State University

(Nearly) Scale invariant fluid dynamics for the dilute Fermi gas in two and three dimensions. Thomas Schaefer North Carolina State University (Nearly) Scale invariant fluid dynamics for the dilute Fermi gas in two and three dimensions Thomas Schaefer North Carolina State University Outline I. Conformal hydrodynamics II. Observations (3d) III.

More information

Lower Tail Probabilities and Normal Comparison Inequalities. In Memory of Wenbo V. Li s Contributions

Lower Tail Probabilities and Normal Comparison Inequalities. In Memory of Wenbo V. Li s Contributions Lower Tail Probabilities and Normal Comparison Inequalities In Memory of Wenbo V. Li s Contributions Qi-Man Shao The Chinese University of Hong Kong Lower Tail Probabilities and Normal Comparison Inequalities

More information

ELEMENTS OF PROBABILITY THEORY

ELEMENTS 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 information

Lecture Notes 3 Convergence (Chapter 5)

Lecture Notes 3 Convergence (Chapter 5) Lecture Notes 3 Convergence (Chapter 5) 1 Convergence of Random Variables Let X 1, X 2,... be a sequence of random variables and let X be another random variable. Let F n denote the cdf of X n and let

More information

Perfect simulation algorithm of a trajectory under a Feynman-Kac law

Perfect simulation algorithm of a trajectory under a Feynman-Kac law Perfect simulation algorithm of a trajectory under a Feynman-Kac law Data Assimilation, 24th-28th September 212, Oxford-Man Institute C. Andrieu, N. Chopin, A. Doucet, S. Rubenthaler University of Bristol,

More information

Practical Bayesian Optimization of Machine Learning. Learning Algorithms

Practical Bayesian Optimization of Machine Learning. Learning Algorithms Practical Bayesian Optimization of Machine Learning Algorithms CS 294 University of California, Berkeley Tuesday, April 20, 2016 Motivation Machine Learning Algorithms (MLA s) have hyperparameters that

More information

Stat 150 Practice Final Spring 2015

Stat 150 Practice Final Spring 2015 Stat 50 Practice Final Spring 205 Instructor: Allan Sly Name: SID: There are 8 questions. Attempt all questions and show your working - solutions without explanation will not receive full credit. Answer

More information

Lecture 4: Introduction to stochastic processes and stochastic calculus

Lecture 4: Introduction to stochastic processes and stochastic calculus Lecture 4: Introduction to stochastic processes and stochastic calculus Cédric Archambeau Centre for Computational Statistics and Machine Learning Department of Computer Science University College London

More information

The Sherrington-Kirkpatrick model: an overview

The Sherrington-Kirkpatrick model: an overview The Sherrington-Kirkpatrick model: an overview Dmitry Panchenko arxiv:1211.1094v1 [math-ph] 6 ov 2012 ovember 7, 2012 Abstract The goal of this paper is to review some of the main ideas that emerged from

More information

Mean field simulation for Monte Carlo integration. Part II : Feynman-Kac models. P. Del Moral

Mean field simulation for Monte Carlo integration. Part II : Feynman-Kac models. P. Del Moral Mean field simulation for Monte Carlo integration Part II : Feynman-Kac models P. Del Moral INRIA Bordeaux & Inst. Maths. Bordeaux & CMAP Polytechnique Lectures, INLN CNRS & Nice Sophia Antipolis Univ.

More information

Statistics, Data Analysis, and Simulation SS 2013

Statistics, Data Analysis, and Simulation SS 2013 Statistics, Data Analysis, and Simulation SS 2013 08.128.730 Statistik, Datenanalyse und Simulation Dr. Michael O. Distler Mainz, May 21, 2013 3. Parameter estimation 1 Consistency:

More information

Physics 212: Statistical mechanics II Lecture XI

Physics 212: Statistical mechanics II Lecture XI Physics 212: Statistical mechanics II Lecture XI The main result of the last lecture was a calculation of the averaged magnetization in mean-field theory in Fourier space when the spin at the origin is

More information

Consistency of the maximum likelihood estimator for general hidden Markov models

Consistency of the maximum likelihood estimator for general hidden Markov models Consistency of the maximum likelihood estimator for general hidden Markov models Jimmy Olsson Centre for Mathematical Sciences Lund University Nordstat 2012 Umeå, Sweden Collaborators Hidden Markov models

More information

Nonparametric Indepedence Tests: Space Partitioning and Kernel Approaches

Nonparametric Indepedence Tests: Space Partitioning and Kernel Approaches Nonparametric Indepedence Tests: Space Partitioning and Kernel Approaches Arthur Gretton 1 and László Györfi 2 1. Gatsby Computational Neuroscience Unit London, UK 2. Budapest University of Technology

More information

Molecular Theories of Linear Viscoelasticity THE ROUSE MODEL (P. 1)

Molecular Theories of Linear Viscoelasticity THE ROUSE MODEL (P. 1) THE ROUSE ODEL (P. 1) odel polymer dynamics by a system of N + 1 beads connected by N springs. Figure 1: apping the Polymer Chain onto a Chain of Beads Connected by Springs. ROUSE SCALING Recall that a

More information

Brownian Motion. Chapter Definition of Brownian motion

Brownian Motion. Chapter Definition of Brownian motion Chapter 5 Brownian Motion Brownian motion originated as a model proposed by Robert Brown in 1828 for the phenomenon of continual swarming motion of pollen grains suspended in water. In 1900, Bachelier

More information

1 Probability and Random Variables

1 Probability and Random Variables 1 Probability and Random Variables The models that you have seen thus far are deterministic models. For any time t, there is a unique solution X(t). On the other hand, stochastic models will result in

More information

Asymptotics for posterior hazards

Asymptotics for posterior hazards Asymptotics for posterior hazards Pierpaolo De Blasi University of Turin 10th August 2007, BNR Workshop, Isaac Newton Intitute, Cambridge, UK Joint work with Giovanni Peccati (Université Paris VI) and

More information

Anomalous 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 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 information

Exercises with solutions (Set D)

Exercises with solutions (Set D) Exercises with solutions Set D. A fair die is rolled at the same time as a fair coin is tossed. Let A be the number on the upper surface of the die and let B describe the outcome of the coin toss, where

More information

Question Points Score Total: 70

Question Points Score Total: 70 The University of British Columbia Final Examination - April 204 Mathematics 303 Dr. D. Brydges Time: 2.5 hours Last Name First Signature Student Number Special Instructions: Closed book exam, no calculators.

More information

Numerical Methods with Lévy Processes

Numerical 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 information

Winter 2019 Math 106 Topics in Applied Mathematics. Lecture 9: Markov Chain Monte Carlo

Winter 2019 Math 106 Topics in Applied Mathematics. Lecture 9: Markov Chain Monte Carlo Winter 2019 Math 106 Topics in Applied Mathematics Data-driven Uncertainty Quantification Yoonsang Lee (yoonsang.lee@dartmouth.edu) Lecture 9: Markov Chain Monte Carlo 9.1 Markov Chain A Markov Chain Monte

More information

Introduction to Stochastic processes

Introduction to Stochastic processes Università di Pavia Introduction to Stochastic processes Eduardo Rossi Stochastic Process Stochastic Process: A stochastic process is an ordered sequence of random variables defined on a probability space

More information

Logarithmic scaling of planar random walk s local times

Logarithmic scaling of planar random walk s local times Logarithmic scaling of planar random walk s local times Péter Nándori * and Zeyu Shen ** * Department of Mathematics, University of Maryland ** Courant Institute, New York University October 9, 2015 Abstract

More information

XVI International Congress on Mathematical Physics

XVI International Congress on Mathematical Physics Aug 2009 XVI International Congress on Mathematical Physics Underlying geometry: finite graph G=(V,E ). Set of possible configurations: V (assignments of +/- spins to the vertices). Probability of a configuration

More information

Introduction to Information Entropy Adapted from Papoulis (1991)

Introduction to Information Entropy Adapted from Papoulis (1991) Introduction to Information Entropy Adapted from Papoulis (1991) Federico Lombardo Papoulis, A., Probability, Random Variables and Stochastic Processes, 3rd edition, McGraw ill, 1991. 1 1. INTRODUCTION

More information

Week 1 Quantitative Analysis of Financial Markets Distributions A

Week 1 Quantitative Analysis of Financial Markets Distributions A Week 1 Quantitative Analysis of Financial Markets Distributions A Christopher Ting http://www.mysmu.edu/faculty/christophert/ Christopher Ting : christopherting@smu.edu.sg : 6828 0364 : LKCSB 5036 October

More information

Complexity of the quantum adiabatic algorithm

Complexity of the quantum adiabatic algorithm Complexity of the quantum adiabatic algorithm Peter Young e-mail:peter@physics.ucsc.edu Collaborators: S. Knysh and V. N. Smelyanskiy Colloquium at Princeton, September 24, 2009 p.1 Introduction What is

More information

Mixing times via super-fast coupling

Mixing times via super-fast coupling s via super-fast coupling Yevgeniy Kovchegov Department of Mathematics Oregon State University Outline 1 About 2 Definition Coupling Result. About Outline 1 About 2 Definition Coupling Result. About Coauthor

More information

L n = l n (π n ) = length of a longest increasing subsequence of π n.

L n = l n (π n ) = length of a longest increasing subsequence of π n. Longest increasing subsequences π n : permutation of 1,2,...,n. L n = l n (π n ) = length of a longest increasing subsequence of π n. Example: π n = (π n (1),..., π n (n)) = (7, 2, 8, 1, 3, 4, 10, 6, 9,

More information

On Markov Chain Monte Carlo

On Markov Chain Monte Carlo MCMC 0 On Markov Chain Monte Carlo Yevgeniy Kovchegov Oregon State University MCMC 1 Metropolis-Hastings algorithm. Goal: simulating an Ω-valued random variable distributed according to a given probability

More information

Markov operators, classical orthogonal polynomial ensembles, and random matrices

Markov operators, classical orthogonal polynomial ensembles, and random matrices Markov operators, classical orthogonal polynomial ensembles, and random matrices M. Ledoux, Institut de Mathématiques de Toulouse, France 5ecm Amsterdam, July 2008 recent study of random matrix and random

More information

This exam contains 6 questions. The questions are of equal weight. Print your name at the top of this page in the upper right hand corner.

This exam contains 6 questions. The questions are of equal weight. Print your name at the top of this page in the upper right hand corner. GROUND RULES: This exam contains 6 questions. The questions are of equal weight. Print your name at the top of this page in the upper right hand corner. This exam is closed book and closed notes. Show

More information

Filtrations, Markov Processes and Martingales. Lectures on Lévy Processes and Stochastic Calculus, Braunschweig, Lecture 3: The Lévy-Itô Decomposition

Filtrations, Markov Processes and Martingales. Lectures on Lévy Processes and Stochastic Calculus, Braunschweig, Lecture 3: The Lévy-Itô Decomposition Filtrations, Markov Processes and Martingales Lectures on Lévy Processes and Stochastic Calculus, Braunschweig, Lecture 3: The Lévy-Itô Decomposition David pplebaum Probability and Statistics Department,

More information