Bounds for (generalised) Lyapunov exponents for deterministic and random products of shears
|
|
- Emil Cain
- 5 years ago
- Views:
Transcription
1 Bounds for (generalised) Lyapunov exponents for deterministic and random products of shears Rob Sturman School of Mathematics University of Leeds Applied & Computational Mathematics seminar, 15 March 2017 University of Edinburgh Joint work with James Springham (Leeds) Jean-Luc Thiffeault (University of Wisconsin-Madison) Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
2 Periodic, deterministic, random Motivation fluid mixing Motivation Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
3 Periodic, deterministic, random Motivation fluid mixing Motivation Basic model is a combination of shears: ( ) ( ) G = & F = Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
4 Periodic, deterministic, random Periodic The Arnold Cat Map compose F and G H = G F (a toral automorphism) ( ) ( ) = ( ) is hyperbolic + Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
5 Periodic, deterministic, random Periodic The Arnold Cat Map compose F and G H = G F (a toral automorphism) ( ) ( ) = Lyapunov exponents ( ) is hyperbolic + 1 lim n n log DHn x v = log(3 + 5)/2 Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
6 Periodic, deterministic, random Periodic The Arnold Cat Map compose F and G H = G F (a toral automorphism) ( ) ( ) = Lyapunov exponents ( ) is hyperbolic + 1 lim n n log DHn x v = log(3 + 5)/2 2 The Cat Map is strong mixing: A, B, lim n µ(hn (A) B) = µ(a)µ(b) Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
7 Periodic, deterministic, random Periodic The Arnold Cat Map compose F and G H = G F (a toral automorphism) ( ) ( ) = Lyapunov exponents ( ) is hyperbolic + 1 lim n n log DHn x v = log(3 + 5)/2 2 The Cat Map is strong mixing: A, B, lim n µ(hn (A) B) = µ(a)µ(b) 3 Exponential decay of correlations (φ H n )ψdµ φdµ ψdµ = O(θ n ) Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
8 Periodic, deterministic, random Deterministic Linked twist maps on the torus Id P = {(x, y) y 1/2}, Q = {(x, y) x 1/2} Id and define F(x, y) = G(x, y) = { (x + 2y, y) if (x, y) P (x, y) { (x, y + 2x) (x, y) if x / P if (x, y) Q if x / Q S And finally H(x, y) = G F(x, y) Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
9 Periodic, deterministic, random Deterministic LTM S Id Id A trajectory always landing in overlap S is uniformly hyperbolic Almost all trajectories land outside S The closer to a boundary, the longer the return to S Arbitrarily long sequences of Fs and Gs But these are rare enough that Lyapunov exponent λ is non-zero No uniform growth = algebraic decay of correlations (polynomial, not exponential, mixing) Computing value of λ is hard Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
10 Periodic, deterministic, random LTM return time partition Deterministic Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
11 Periodic, deterministic, random Random matrix products Random G = ( ) & F = ( ) At each iterate, choose F or G with probability 1/2. Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
12 Periodic, deterministic, random Random Random braiding of particle paths Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
13 Periodic, deterministic, random Tangled magnetic fields Random Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
14 Periodic, deterministic, random Lyapunov exponents Definition of Lyapunov exponents Diffeomorphisms: 1 λ(x, v) = lim n n log D xh n v Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
15 Periodic, deterministic, random Lyapunov exponents Definition of Lyapunov exponents Diffeomorphisms: 1 λ(x, v) = lim n n log D xh n v For products of random matrices M N = N k=1 A i k : 1 λ = lim N N E log M 1 N = lim N N E log X N almost surely, where X N = A in X N 1. Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
16 Periodic, deterministic, random Lyapunov exponents Definition of Lyapunov exponents Diffeomorphisms: 1 λ(x, v) = lim n n log D xh n v For products of random matrices M N = N k=1 A i k : 1 λ = lim N N E log M 1 N = lim N N E log X N almost surely, where X N = A in X N 1. Convergence is given by the celebrated Furstenberg-Kesten theorem (1960). Equivalence of the two definitions is given by the Oseledets multiplicative ergodic theorem (1968) λ is independent of choice of norm Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
17 Periodic, deterministic, random Hard to compute analytically... Lyapunov exponents There is a paucity of exact results concerning Lyapunov exponents for random matrices, as famously lamented by Kingman (1973): Pride of place among the unsolved problems of subadditive ergodic theory must go to the calculation of [the Lyapunov exponent]... In none of the applications described here is there an obvious mechanism for obtaining an exact numerical value, and indeed this usually seems to be a problem of some depth. Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
18 Periodic, deterministic, random Existing work... but easy computationally The submultiplicativity of matrix norms provides the most popular upper bound in the literature : E k = 2 k E {log M 2 k }, (for products drawn from 2 possible matrices) The E k decrease monotonically to λ (for any matrix norm) as k This is easy, if not efficient the number of matrix product calculations required increases exponentially with k Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
19 Periodic, deterministic, random Existing work Many (very good) algorithms and approximations... Mannion, Products of 2 2 random matrices (1993) Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
20 Periodic, deterministic, random Existing work Many (very good) algorithms and approximations... Mannion, Products of 2 2 random matrices (1993) Pollicott, Inventiones mathematicae, 181, 1 (2010): Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
21 Periodic, deterministic, random Existing work Many (very good) algorithms and approximations... Mannion, Products of 2 2 random matrices (1993) Pollicott, Inventiones mathematicae, 181, 1 (2010): Protasov & Jungers, Lower and upper bounds for the largest Lyapunov exponent of matrices (2013) Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
22 Periodic, deterministic, random Our contribution Why is this different? Rigorous, explicit (elementary) bounds (Not algorithmic, no complex functions, no Fourier decomposition) Positive and negative entries in matrices (but see also D. Viswanath, Random Fibonacci sequences and the number Mathematics of Computation of the American Mathematical Society (2000)) Generalised Lyapunov exponents Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
23 Periodic, deterministic, random Generalised Lyapunov exponents Generalised Lyapunov exponents Often want more than knowledge of the infinite time limit of matrix product growth Study deviations from the long-time average Growth rate of the qth moment of the matrix product norm Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
24 Periodic, deterministic, random Generalised Lyapunov exponents Generalised Lyapunov exponents Often want more than knowledge of the infinite time limit of matrix product growth Study deviations from the long-time average Growth rate of the qth moment of the matrix product norm Define the generalized Lyapunov exponents as 1 l(q) = lim N N log E M N q, Note that l (0) is the usual Lyapunov exponent l(1) measures topological entropy (in the sense of growth of material lines) Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
25 Periodic, deterministic, random Generalised Lyapunov exponents Generalised Lyapunov exponents are hard to compute J. Vanneste, Estimating generalized Lyapunov exponents for products of random matrices. Physical Review E 81.3 (2010) For large q, l(q) is controlled by exceedingly rare realisations of the matrix products, and hence it is difficult to estimate reliably using Monte Carlo numerical methods, even with importance sampling. Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
26 Bounding Lyapunov exponents Express the problem differently Group the matrices together W.l.o.g. assume that the first matrix in the product is F Group F s and G s together into J blocks e.g. (FFGGG)(FG)(FGG)(FFFFGG)(FGG)... so the random product is J M J = F a j G b j, a j + b j = n j, j=1 J n j = N, j=1 with 1 a i, b i < n i, so n i 2. Now the a i and b i are the i.i.d. random variables, with identical probability distribution P(x) = 2 x, x 1. Hence, the joint distribution P(a, b) = P(a) P(b) = 2 (a+b). We have the expected values Ea = Eb = 2, so En = 4 Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
27 Bounding Lyapunov exponents Express the problem differently Group the matrices together Recall 1 λ = lim N N E log X N almost surely, where X N = FX N 1 or GX N 1. Then where 1 λ = lim J 4J E log X J X J = F a J G b J X J 1 X J 1 F a J 1G b J 1X J 2 X J 2 F a 1G b 1X 0. X 0 We will compute bounds on each term in the RHS product, by bounding the orientation of the vectors X i, and recalling that the choice of norm is arbitrary. Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
28 Shears Bounding Lyapunov exponents The simple case We take the specific matrices ( ) 1 0 F =, G = 1 1 and let ( ) 1 1, 0 1 ( ) M ab = F a G b 1 b =. a 1 + ab A more general problem takes the matrices F = ( ) 1 0, G = α 1 ( ) ( ) 1 β, M 0 1 ab = F a G b 1 bβ =. aα 1 + aαbβ Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
29 Bounding Lyapunov exponents Invariant cones Simple case (single, positive shears) Take v u/v = 1 F = ( ) 1 0, G = 1 1 ( ) u Lemma The cone C = {0 u/v 1} in tangent space is invariant under M ab = F a G b for all a, b 1. Moreover, this is the minimal such cone. Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
30 Norm bounds I. Bounding Lyapunov exponents Invariant cones Lemma (L -norm) The norm M ab X for a vector X C satisfies the bounds 1 + ab M abx X 1 + a + ab Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
31 Norm bounds I. Bounding Lyapunov exponents Invariant cones Lemma (L -norm) The norm M ab X for a vector X C satisfies the bounds 1 + ab M abx X 1 + a + ab Lemma (L 1 -norm) The norm M ab X 1 for a vector X C satisfies the bounds (a + b + ab) M abx 1 X b + ab Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
32 Norm bounds II. Bounding Lyapunov exponents Invariant cones Lemma (L 2 -norm) The norm M ab X 2 for a vector X C satisfies the bounds M ab X 2 2 X ( 2 + C ab + ) C ab (C ab + 4), where and M ab X 2 2 X 2 2 C ab = (a + b) 2 + a 2 b 2. { (1 + ab) 2 + b 2 min ( (1 + a + ab) 2 + (1 + b) 2) 1 2 Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
33 Putting it together Bounding Lyapunov exponents Invariant cones 1 λ = lim J 4J E log X J 1 = lim J 4J E log F aj G bj X J 1 X J 1 F a 1G b 1X 0 X 0 Using the geometric probability distribution for the a s and b s, we then have, for the L -norm: a,b=1 2 a b log(1 + ab) 4λ a,b=1 and similar expressions for the other two norms. 2 a b log(1 + a + ab) Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
34 Accuracy Accuracy in the simple case Simple case The lowest upper bound (U 2 ) and largest lower bound (L 1 ) differ by about 8%. The true value (from explicit calculation via a standard algorithm) is Norm Lower bound Upper bound L 1 L 1 (1, 1) = U 1 (1, 1) = L 2 L 2 (1, 1) = U 2 (1, 1) = L L (1, 1) = U (1, 1) = Table: Bounds for the maximal Lyapunov exponent for the matrix product in the case α = β = 1. Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
35 Accuracy Positive shears of α and β General shears The invariant cone is now C = {0 u/v 1/α} Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
36 Accuracy Positive shears of α and β General shears The invariant cone is now Theorem C = {0 u/v 1/α} The Lyapunov exponent λ(α, β) for the product M N satisfies max L k(α, β) 4λ(α, β) min U k(α, β) k {1,2, } k {1,2, } where L k (α, β) = U k (α, β) = 2 a b log φ k (a, b, α, β) a,b=1 a,b=1 2 a b log ψ k (a, b, α, β), Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
37 Accuracy General shears φ 1 (a, b, α, β) = 1 + α 1+α (a + bβ + aαbβ) ((1 + aαbβ) 2 + b 2 β 2 ) 1/2 φ 2 (a, b, α, β) = min ( 1 ( 1+α α 2 (1 + a + aαbβ) 2 + (1 + αbβ) 2)) 1/2 { aαbβ φ (a, b, α, β) = min {max(1 + aαbβ, bβ), max(α(1 + a + aαbβ), 1 + ψ 1 (a, b, α, β) = 1 + bβ + aαbβ ( ) 1/2 ψ 2 (a, b, α, β) = 2 + C aαbβ + C aαbβ (C aαbβ + 4) { ψ (a, b, α, β) = (1 + a + aαbβ) α 1 max(1 + aα + aαbβ, 1 + bβ) α < 1 where C aαbβ = (aα + bβ) 2 + (aαbβ) 2. Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
38 Accuracy of the bounds Accuracy General shears λ α = β Linf upper Linf lower L1 upper L1 lower L2 upper L2 lower standard Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
39 Accuracy Asymptotic explicit bounds General shears Corollary The Lyapunov exponent λ(α, β) for the product M N for α, β 1 satisfies K + log αβ 4λ K + log( αβ + 1 ) + 1 log(1 + αβ), αβ 2 where K = 2 a b log ab a,b=1 Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
40 Accuracy of the bounds Accuracy General shears Bounds λ α = β Linf upper L2 upper L1 upper Linf lower L2 lower L1 lower Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
41 Accuracy of the bounds Accuracy General shears Upper bound - Lower bound Linf L1 L2 Combined α = β Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
42 Accuracy Improving the bounds Comparing two independent geometric distributions Lemma When a and b are both drawn from i.i.d. geometric distributions with parameter 1/2, we have P(a = b) = P(a > b) = P(b > a) = 1/3. Proof. We have P(a = b) = P(a = i b = i) = i=1 i=1 2 2i = 1/4 1 1/4 = 1 3. The remaining two equalities follow by symmetry. Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
43 Accuracy Improving the bounds The cone can be made smaller Lemma The cone C = {0 u v 1 α } is mapped into the following cones, in the following cases: 1 When a = b, M ab (C) = {0 u v 1+αβ α(2+αβ) }; 2 when a > b, M ab (C) = {0 u v 1+αβ α(3+2αβ) }; 3 when a < b, M ab (C) = C. and hence expressions like 4λ a,b=1 2 a b log 1 (ψ(a = b) + ψ(a > b) + ψ(a < b)) 3 Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
44 Accuracy Improving the bounds Accuracy of the bounds, improved Upper bound - Lower bound Linf L1 L2 Linf improved L1 improved L2 improved Combined Combined, improved α = β Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
45 Accuracy Improving the bounds Generalised Lyapunov exponents Theorem We have 4l(q, α, β) 4l(q, α, β) max k {1,2, } max k {1,2, } log log with φ k, ψ k and defined as above. a,b=1 a,b=1 2 a b (φ k (a, b, α, β)) q 2 a b (ψ k (a, b, α, β)) q Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
46 Accuracy Generalised Lyapunov exponent Improving the bounds l(q) q Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
47 Accuracy Generalised Lyapunov exponents Generalised Lyapunov exponents When q is a positive integer we can evaluate this easily by expanding the power q. Since ( ) ( ) 2 a b a n b m = 2 a a n 2 b b m a,b=1 a=1 b=1 we require values of the polylogarithm Li q ( 1 2 ), defined by z k Li s (z) = k s. k=1 For integer q = s we have special values 2 a a n = 1, 2, 6, 26, 150, 1082, 9366,... for n = 0, 1, 2, 3, 4, 5, 6,... a=1 and so the L norm, for example, gives... Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
48 Accuracy Generalised Lyapunov exponents Generalised Lyapunov exponents Corollary Generalised Lyapunov exponents in the case α = β = 1 are bounded by: 1 4 log 5 l(1) 1 4 log log 45 l(2) 1 4 log log 797 l(3) 1 4 log log l(4) 1 4 log log l(5) 1 4 log Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
49 Accuracy Switch one of the shears Reversing the twist Take F = ( ) 1 0, G = α 1 ( ) 1 β 0 1 u/v = 1 v Now ( ) M ab = F a G b 1 bβ =. aα 1 aαbβ u is only hyperbolic a, b when αβ > 4. We ll assume α, β > 2 for ease. Lemma The cone C = { 1 u/v 0} in tangent space is invariant under M ab = F a G b for all a, b 1. Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
50 Accuracy of the bounds Accuracy Reversing the twist 0.20 Ũ 1 L Ũ 2 L 2 Ũ L Ũ 1 ˆL Ũ 2 ˆL 2 Ũ ˆL Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
51 Conclusions Conclusions Alternating shears in (trivially) uniformly hyperbolic LTMs are more interesting non-zero Lyapunov exponents We can derive pretty accurate, explicit, upper and lower bounds for Lyapunov exponents for random products of shears of varying strength Generalised Lyapunov exponents can also be bounded (which are also hard to compute) Our results, unusually, include matrices with negative entries Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
52 Conclusions Next?... Linked twist maps also produce a sequence FFGFFGFGGGFFGFGGGF... = J F a j G b j but now the a j, b j are no longer random we don t have a probability distribution. But a block (F a G b ) represents the return map to the hyperbolic overlap regions... and such a block still possess the invariant cone and bounds on the growth of the vector norm... and Kac s lemma states that given a set Q with measure µ(q), the expected return time is 1/µ(Q)... j=1 so this leads us back to the return time partition. Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
53 Next?... Conclusions x 1 = y 1 = 1/2 Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
54 Next?... Conclusions x 1 = y 1 = 1/2 x 1 = y 1 = 3/4 Rob Sturman (University of Leeds) Mixing for products of shears 15 March / 44
Bounds for (generalised) Lyapunov exponents for random products of shears
Bounds for (generalised) Lyapunov exponents for random products of shears Rob Sturman Department of Mathematics University of Leeds LAND Seminar, 13 June 2016 Leeds Joint work with Jean-Luc Thiffeault
More informationConcentration inequalities for linear cocycles and their applications to problems in dynamics and mathematical physics *
Concentration inequalities for linear cocycles and their applications to problems in dynamics and mathematical physics * Silvius Klein Pontifical Catholic University of Rio de Janeiro (PUC-Rio), Brazil
More informationEntropy, mixing, and independence
Entropy, mixing, and independence David Kerr Texas A&M University Joint work with Hanfeng Li Let (X, µ) be a probability space. Two sets A, B X are independent if µ(a B) = µ(a)µ(b). Suppose that we have
More informationTHE MULTIPLICATIVE ERGODIC THEOREM OF OSELEDETS
THE MULTIPLICATIVE ERGODIC THEOREM OF OSELEDETS. STATEMENT Let (X, µ, A) be a probability space, and let T : X X be an ergodic measure-preserving transformation. Given a measurable map A : X GL(d, R),
More information5 Birkhoff s Ergodic Theorem
5 Birkhoff s Ergodic Theorem Birkhoff s Ergodic Theorem extends the validity of Kolmogorov s strong law to the class of stationary sequences of random variables. Stationary sequences occur naturally even
More informationSRB measures for non-uniformly hyperbolic systems
SRB measures for non-uniformly hyperbolic systems Vaughn Climenhaga University of Maryland October 21, 2010 Joint work with Dmitry Dolgopyat and Yakov Pesin 1 and classical results Definition of SRB measure
More informationRobustly transitive diffeomorphisms
Robustly transitive diffeomorphisms Todd Fisher tfisher@math.byu.edu Department of Mathematics, Brigham Young University Summer School, Chengdu, China 2009 Dynamical systems The setting for a dynamical
More informationMET Workshop: Exercises
MET Workshop: Exercises Alex Blumenthal and Anthony Quas May 7, 206 Notation. R d is endowed with the standard inner product (, ) and Euclidean norm. M d d (R) denotes the space of n n real matrices. When
More informationLyapunov Exponents of Linear Cocycles *
Lyapunov Exponents of Linear Cocycles * Continuity via Large Deviations Silvius Klein Norwegian University of Science and Technology (NTNU) *based on joint work with Pedro Duarte, University of Lisbon
More informationLecture Notes Introduction to Ergodic Theory
Lecture Notes Introduction to Ergodic Theory Tiago Pereira Department of Mathematics Imperial College London Our course consists of five introductory lectures on probabilistic aspects of dynamical systems,
More information7.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 informationDYNAMICAL SYSTEMS PROBLEMS. asgor/ (1) Which of the following maps are topologically transitive (minimal,
DYNAMICAL SYSTEMS PROBLEMS http://www.math.uci.edu/ asgor/ (1) Which of the following maps are topologically transitive (minimal, topologically mixing)? identity map on a circle; irrational rotation of
More informationL(C G (x) 0 ) c g (x). Proof. Recall C G (x) = {g G xgx 1 = g} and c g (x) = {X g Ad xx = X}. In general, it is obvious that
ALGEBRAIC GROUPS 61 5. Root systems and semisimple Lie algebras 5.1. Characteristic 0 theory. Assume in this subsection that chark = 0. Let me recall a couple of definitions made earlier: G is called reductive
More informationRational Exponents. Polynomial function of degree n: with leading coefficient,, with maximum number of turning points is given by (n-1)
Useful Fact Sheet Final Exam Interval, Set Builder Notation (a,b) = {x a
More informationAverage laws in analysis
Average laws in analysis Silvius Klein Norwegian University of Science and Technology (NTNU) The law of large numbers: informal statement The theoretical expected value of an experiment is approximated
More informationPeriodic Sinks and Observable Chaos
Periodic Sinks and Observable Chaos Systems of Study: Let M = S 1 R. T a,b,l : M M is a three-parameter family of maps defined by where θ S 1, r R. θ 1 = a+θ +Lsin2πθ +r r 1 = br +blsin2πθ Outline of Contents:
More informationSmooth Ergodic Theory and Nonuniformly Hyperbolic Dynamics
CHAPTER 2 Smooth Ergodic Theory and Nonuniformly Hyperbolic Dynamics Luis Barreira Departamento de Matemática, Instituto Superior Técnico, 1049-001 Lisboa, Portugal E-mail: barreira@math.ist.utl.pt url:
More informationON RIGIDITY OF ALGEBRAIC DYNAMICAL SYSTEMS
ON RIGIDITY OF ALGEBRAIC DYNAMICAL SYSTEMS ALEX CLARK AND ROBBERT FOKKINK Abstract. We study topological rigidity of algebraic dynamical systems. In the first part of this paper we give an algebraic condition
More informationProblems in hyperbolic dynamics
Problems in hyperbolic dynamics Current Trends in Dynamical Systems and the Mathematical Legacy of Rufus Bowen Vancouver july 31st august 4th 2017 Notes by Y. Coudène, S. Crovisier and T. Fisher 1 Zeta
More informationTEST CODE: PMB SYLLABUS
TEST CODE: PMB SYLLABUS Convergence and divergence of sequence and series; Cauchy sequence and completeness; Bolzano-Weierstrass theorem; continuity, uniform continuity, differentiability; directional
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 informationThe Multivariate Gaussian Distribution [DRAFT]
The Multivariate Gaussian Distribution DRAFT David S. Rosenberg Abstract This is a collection of a few key and standard results about multivariate Gaussian distributions. I have not included many proofs,
More informationChaos and Liapunov exponents
PHYS347 INTRODUCTION TO NONLINEAR PHYSICS - 2/22 Chaos and Liapunov exponents Definition of chaos In the lectures we followed Strogatz and defined chaos as aperiodic long-term behaviour in a deterministic
More informationLimits at Infinity. Horizontal Asymptotes. Definition (Limits at Infinity) Horizontal Asymptotes
Limits at Infinity If a function f has a domain that is unbounded, that is, one of the endpoints of its domain is ±, we can determine the long term behavior of the function using a it at infinity. Definition
More informationPhysical measures of discretizations of generic diffeomorphisms
Ergod. Th. & Dynam. Sys. (2018), 38, 1422 1458 doi:10.1017/etds.2016.70 c Cambridge University Press, 2016 Physical measures of discretizations of generic diffeomorphisms PIERRE-ANTOINE GUIHÉNEUF Université
More informationApril 13, We now extend the structure of the horseshoe to more general kinds of invariant. x (v) λ n v.
April 3, 005 - Hyperbolic Sets We now extend the structure of the horseshoe to more general kinds of invariant sets. Let r, and let f D r (M) where M is a Riemannian manifold. A compact f invariant set
More informationBalanced Truncation 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 2004: MODEL REDUCTION Balanced Truncation This lecture introduces balanced truncation for LTI
More informationNUMERICALLY COMPUTING THE LYAPUNOV EXPONENTS OF MATRIX-VALUED COCYCLES
NUMERICALLY COMPUTING THE LYAPUNOV EXPONENTS OF MATRIX-VALUED COCYCLES RODRIGO TREVIÑO This short note is based on a talk I gave at the student dynamical systems seminar about using your computer to figure
More informationMatrix Energy. 1 Graph Energy. Christi DiStefano Gary Davis CSUMS University of Massachusetts at Dartmouth. December 16,
Matrix Energy Christi DiStefano Gary Davis CSUMS University of Massachusetts at Dartmouth December 16, 2009 Abstract We extend Ivan Gutmans idea of graph energy, stemming from theoretical chemistry via
More information8. Limit Laws. lim(f g)(x) = lim f(x) lim g(x), (x) = lim x a f(x) g lim x a g(x)
8. Limit Laws 8.1. Basic Limit Laws. If f and g are two functions and we know the it of each of them at a given point a, then we can easily compute the it at a of their sum, difference, product, constant
More informationChapter III. Stability of Linear Systems
1 Chapter III Stability of Linear Systems 1. Stability and state transition matrix 2. Time-varying (non-autonomous) systems 3. Time-invariant systems 1 STABILITY AND STATE TRANSITION MATRIX 2 In this chapter,
More informationSTAT 7032 Probability Spring Wlodek Bryc
STAT 7032 Probability Spring 2018 Wlodek Bryc Created: Friday, Jan 2, 2014 Revised for Spring 2018 Printed: January 9, 2018 File: Grad-Prob-2018.TEX Department of Mathematical Sciences, University of Cincinnati,
More informationRepresentations of moderate growth Paul Garrett 1. Constructing norms on groups
(December 31, 2004) Representations of moderate growth Paul Garrett Representations of reductive real Lie groups on Banach spaces, and on the smooth vectors in Banach space representations,
More informationMA257: INTRODUCTION TO NUMBER THEORY LECTURE NOTES
MA257: INTRODUCTION TO NUMBER THEORY LECTURE NOTES 2018 57 5. p-adic Numbers 5.1. Motivating examples. We all know that 2 is irrational, so that 2 is not a square in the rational field Q, but that we can
More informationSOME EXAMPLES OF NONUNIFORMLY HYPERBOLIC COCYCLES. L.-S. Young 1
SOME EXAMPLES OF NONUNIFORMLY HYPERBOLIC COCYCLES L.-S. Young Abstract. We consider some very simple examples of SL(2, R)-cocycles and prove that they have positive Lyapunov exponents. These cocycles form
More informationMinimal basis for connected Markov chain over 3 3 K contingency tables with fixed two-dimensional marginals. Satoshi AOKI and Akimichi TAKEMURA
Minimal basis for connected Markov chain over 3 3 K contingency tables with fixed two-dimensional marginals Satoshi AOKI and Akimichi TAKEMURA Graduate School of Information Science and Technology University
More informationProblem: A class of dynamical systems characterized by a fast divergence of the orbits. A paradigmatic example: the Arnold cat.
À È Ê ÇÄÁ Ë ËÌ ÅË Problem: A class of dynamical systems characterized by a fast divergence of the orbits A paradigmatic example: the Arnold cat. The closure of a homoclinic orbit. The shadowing lemma.
More informationLUCK S THEOREM ALEX WRIGHT
LUCK S THEOREM ALEX WRIGHT Warning: These are the authors personal notes for a talk in a learning seminar (October 2015). There may be incorrect or misleading statements. Corrections welcome. 1. Convergence
More informationPeriodic approximation of exceptional Lyapunov exponents for sem
Periodic approximation of exceptional Lyapunov exponents for semi-invertible operator cocycles Department of Mathematics, University of Rijeka, Croatia (joint work with Lucas Backes, Porto Alegre, Brazil)
More informationBraids of entangled particle trajectories
Braids of entangled particle trajectories Jean-Luc Thiffeault Department of Mathematics University of Wisconsin Madison Institute for Mathematics and its Applications University of Minnesota Twin Cities
More informationMATH 423 Linear Algebra II Lecture 3: Subspaces of vector spaces. Review of complex numbers. Vector space over a field.
MATH 423 Linear Algebra II Lecture 3: Subspaces of vector spaces. Review of complex numbers. Vector space over a field. Vector space A vector space is a set V equipped with two operations, addition V V
More informationAN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES
AN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES JOEL A. TROPP Abstract. We present an elementary proof that the spectral radius of a matrix A may be obtained using the formula ρ(a) lim
More information6 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 informationarxiv: v1 [math.co] 7 Jul 2014
Sum-ratio estimates over arbitrary finite fields Oliver Roche-Newton arxiv:1407.1654v1 [math.co] 7 Jul 2014 July 11, 2018 Abstract The aim of this note is to record a proof that the estimate max{ A+A,
More informationA. Pikovsky. Deapartment of Physics and Astronomy, Potsdam University. URL: pikovsky
Lyapunov exponents and all that A. Pikovsky Deapartment of Physics and Astronomy, Potsdam University URL: www.stat.physik.uni-potsdam.de/ pikovsky Alexandr Mikhailovich Lyapunov * 1857 Jaroslavl 1870 1876
More informationNonamenable Products are not Treeable
Version of 30 July 1999 Nonamenable Products are not Treeable by Robin Pemantle and Yuval Peres Abstract. Let X and Y be infinite graphs, such that the automorphism group of X is nonamenable, and the automorphism
More informationThe dynamics of mapping classes on surfaces
The dynamics of mapping classes on surfaces Eriko Hironaka May 16, 2013 1 Introduction to mapping classes and the minimum dilatation problem In this section, we define mapping classes on surfaces, and
More informationLINEAR RECURSIVE SEQUENCES. The numbers in the sequence are called its terms. The general form of a sequence is
LINEAR RECURSIVE SEQUENCES BJORN POONEN 1. Sequences A sequence is an infinite list of numbers, like 1) 1, 2, 4, 8, 16, 32,.... The numbers in the sequence are called its terms. The general form of a sequence
More informationCHILE LECTURES: MULTIPLICATIVE ERGODIC THEOREM, FRIENDS AND APPLICATIONS
CHILE LECTURES: MULTIPLICATIVE ERGODIC THEOREM, FRIENDS AND APPLICATIONS. Motivation Context: (Ω, P) a probability space; σ : Ω Ω is a measure-preserving transformation P(σ A) = P(A) for all measurable
More informationPHY411 Lecture notes Part 5
PHY411 Lecture notes Part 5 Alice Quillen January 27, 2016 Contents 0.1 Introduction.................................... 1 1 Symbolic Dynamics 2 1.1 The Shift map.................................. 3 1.2
More informationGroup construction in geometric C-minimal non-trivial structures.
Group construction in geometric C-minimal non-trivial structures. Françoise Delon, Fares Maalouf January 14, 2013 Abstract We show for some geometric C-minimal structures that they define infinite C-minimal
More informationCombining the cycle index and the Tutte polynomial?
Combining the cycle index and the Tutte polynomial? Peter J. Cameron University of St Andrews Combinatorics Seminar University of Vienna 23 March 2017 Selections Students often meet the following table
More informationSYMBOLIC DYNAMICS FOR HYPERBOLIC SYSTEMS. 1. Introduction (30min) We want to find simple models for uniformly hyperbolic systems, such as for:
SYMBOLIC DYNAMICS FOR HYPERBOLIC SYSTEMS YURI LIMA 1. Introduction (30min) We want to find simple models for uniformly hyperbolic systems, such as for: [ ] 2 1 Hyperbolic toral automorphisms, e.g. f A
More informationTopological methods for stirring and mixing
Topological methods for stirring and mixing Jean-Luc Thiffeault Department of Mathematics University of Wisconsin Madison WASCOM 2011 Brindisi, Italy, 16 June 2011 1 / 30 The Taffy Puller This may not
More informationPOISSON PROCESSES 1. THE LAW OF SMALL NUMBERS
POISSON PROCESSES 1. THE LAW OF SMALL NUMBERS 1.1. The Rutherford-Chadwick-Ellis Experiment. About 90 years ago Ernest Rutherford and his collaborators at the Cavendish Laboratory in Cambridge conducted
More informationFree probability and quantum information
Free probability and quantum information Benoît Collins WPI-AIMR, Tohoku University & University of Ottawa Tokyo, Nov 8, 2013 Overview Overview Plan: 1. Quantum Information theory: the additivity problem
More informationRandom variables. DS GA 1002 Probability and Statistics for Data Science.
Random variables DS GA 1002 Probability and Statistics for Data Science http://www.cims.nyu.edu/~cfgranda/pages/dsga1002_fall17 Carlos Fernandez-Granda Motivation Random variables model numerical quantities
More informationEntropy for zero-temperature limits of Gibbs-equilibrium states for countable-alphabet subshifts of finite type
Entropy for zero-temperature limits of Gibbs-equilibrium states for countable-alphabet subshifts of finite type I. D. Morris August 22, 2006 Abstract Let Σ A be a finitely primitive subshift of finite
More information4. Ergodicity and mixing
4. Ergodicity and mixing 4. Introduction In the previous lecture we defined what is meant by an invariant measure. In this lecture, we define what is meant by an ergodic measure. The primary motivation
More informationTopological Dynamics
Topological Dynamics Probing dynamical systems using loops Jean-Luc Thiffeault Department of Mathematics University of Wisconsin Madison Chaos/Xaoc Anniversary Conference, 26 July 2009 Collaborators: Sarah
More informationSection Properties of Rational Expressions
88 Section. - Properties of Rational Expressions Recall that a rational number is any number that can be written as the ratio of two integers where the integer in the denominator cannot be. Rational Numbers:
More informationTHEODORE VORONOV DIFFERENTIABLE MANIFOLDS. Fall Last updated: November 26, (Under construction.)
4 Vector fields Last updated: November 26, 2009. (Under construction.) 4.1 Tangent vectors as derivations After we have introduced topological notions, we can come back to analysis on manifolds. Let M
More informationON HYPERBOLIC MEASURES AND PERIODIC ORBITS
ON HYPERBOLIC MEASURES AND PERIODIC ORBITS ILIE UGARCOVICI Dedicated to Anatole Katok on the occasion of his 60th birthday Abstract. We prove that if a diffeomorphism on a compact manifold preserves a
More information1 Acceptance-Rejection Method
Copyright c 2016 by Karl Sigman 1 Acceptance-Rejection Method As we already know, finding an explicit formula for F 1 (y), y [0, 1], for the cdf of a rv X we wish to generate, F (x) = P (X x), x R, is
More information1.5 Approximate Identities
38 1 The Fourier Transform on L 1 (R) which are dense subspaces of L p (R). On these domains, P : D P L p (R) and M : D M L p (R). Show, however, that P and M are unbounded even when restricted to these
More informationAbstracts. Furstenberg The Dynamics of Some Arithmetically Generated Sequences
CHAOS AND DISORDER IN MATHEMATICS AND PHYSICS Monday 10:00-11:00 Okounkov Algebraic geometry of random surfaces 11:30-12:30 Furstenberg Dynamics of Arithmetically Generated Sequences 12:30-14:30 lunch
More informationALMOST SURE CONVERGENCE OF RANDOM GOSSIP ALGORITHMS
ALMOST SURE CONVERGENCE OF RANDOM GOSSIP ALGORITHMS Giorgio Picci with T. Taylor, ASU Tempe AZ. Wofgang Runggaldier s Birthday, Brixen July 2007 1 CONSENSUS FOR RANDOM GOSSIP ALGORITHMS Consider a finite
More informationAsymptotics of minimax stochastic programs
Asymptotics of minimax stochastic programs Alexander Shapiro Abstract. We discuss in this paper asymptotics of the sample average approximation (SAA) of the optimal value of a minimax stochastic programming
More informationRod motions Topological ingredients Optimization Data analysis Open issues References. Topological chaos. Jean-Luc Thiffeault
Topological chaos Jean-Luc Thiffeault Department of Mathematics University of Wisconsin Madison Chaos and Complex Systems Seminar Physics Dept., University of Wisconsin, Madison 4 October 2011 1 / 29 The
More informationMA 123 (Calculus I) Lecture 3: September 12, 2017 Section A2. Professor Jennifer Balakrishnan,
What is on today Professor Jennifer Balakrishnan, jbala@bu.edu 1 Techniques for computing limits 1 1.1 Limit laws..................................... 1 1.2 One-sided limits..................................
More informationB553 Lecture 1: Calculus Review
B553 Lecture 1: Calculus Review Kris Hauser January 10, 2012 This course requires a familiarity with basic calculus, some multivariate calculus, linear algebra, and some basic notions of metric topology.
More information(II.D) There is exactly one rref matrix row-equivalent to any A
IID There is exactly one rref matrix row-equivalent to any A Recall that A is row-equivalent to B means A = E M E B where E i are invertible elementary matrices which we think of as row-operations In ID
More informationTHE N-VALUE GAME OVER Z AND R
THE N-VALUE GAME OVER Z AND R YIDA GAO, MATT REDMOND, ZACH STEWARD Abstract. The n-value game is an easily described mathematical diversion with deep underpinnings in dynamical systems analysis. We examine
More informationTopology, pseudo-anosov mappings, and fluid dynamics
Topology, pseudo-anosov mappings, and fluid dynamics Jean-Luc Thiffeault Department of Mathematics University of Wisconsin Madison Institute for Mathematics and its Applications University of Minnesota
More informationHypocoercivity and Sensitivity Analysis in Kinetic Equations and Uncertainty Quantification October 2 nd 5 th
Hypocoercivity and Sensitivity Analysis in Kinetic Equations and Uncertainty Quantification October 2 nd 5 th Department of Mathematics, University of Wisconsin Madison Venue: van Vleck Hall 911 Monday,
More informationTwo viewpoints on measure valued processes
Two viewpoints on measure valued processes Olivier Hénard Université Paris-Est, Cermics Contents 1 The classical framework : from no particle to one particle 2 The lookdown framework : many particles.
More informationInvariant measure for random walks on ergodic environments on a strip
Invariant measure for random walks on ergodic environments on a strip Dmitry Dolgopyat and Ilya Goldsheid Department of Mathematics and Institute of Physical Science & Technology University of Maryland
More informationThe Quantum Hall Conductance: A rigorous proof of quantization
Motivation The Quantum Hall Conductance: A rigorous proof of quantization Spyridon Michalakis Joint work with M. Hastings - Microsoft Research Station Q August 17th, 2010 Spyridon Michalakis (T-4/CNLS
More information18.175: Lecture 2 Extension theorems, random variables, distributions
18.175: Lecture 2 Extension theorems, random variables, distributions Scott Sheffield MIT Outline Extension theorems Characterizing measures on R d Random variables Outline Extension theorems Characterizing
More informationEconomics 204 Fall 2011 Problem Set 2 Suggested Solutions
Economics 24 Fall 211 Problem Set 2 Suggested Solutions 1. Determine whether the following sets are open, closed, both or neither under the topology induced by the usual metric. (Hint: think about limit
More informationNotes 6 : First and second moment methods
Notes 6 : First and second moment methods Math 733-734: Theory of Probability Lecturer: Sebastien Roch References: [Roc, Sections 2.1-2.3]. Recall: THM 6.1 (Markov s inequality) Let X be a non-negative
More informationAN INTRODUCTION TO THE FUNDAMENTAL GROUP
AN INTRODUCTION TO THE FUNDAMENTAL GROUP DAVID RAN Abstract. This paper seeks to introduce the reader to the fundamental group and then show some of its immediate applications by calculating the fundamental
More informationWaiting times, recurrence times, ergodicity and quasiperiodic dynamics
Waiting times, recurrence times, ergodicity and quasiperiodic dynamics Dong Han Kim Department of Mathematics, The University of Suwon, Korea Scuola Normale Superiore, 22 Jan. 2009 Outline Dynamical Systems
More informationM4P52 Manifolds, 2016 Problem Sheet 1
Problem Sheet. Let X and Y be n-dimensional topological manifolds. Prove that the disjoint union X Y is an n-dimensional topological manifold. Is S S 2 a topological manifold? 2. Recall that that the discrete
More informationOn the classification of irrational numbers
arxiv:506.0044v [math.nt] 5 Nov 07 On the classification of irrational numbers José de Jesús Hernández Serda May 05 Abstract In this note we make a comparison between the arithmetic properties of irrational
More informationRandomized Algorithms
Randomized Algorithms Prof. Tapio Elomaa tapio.elomaa@tut.fi Course Basics A new 4 credit unit course Part of Theoretical Computer Science courses at the Department of Mathematics There will be 4 hours
More informationMARKOV PARTITIONS FOR HYPERBOLIC SETS
MARKOV PARTITIONS FOR HYPERBOLIC SETS TODD FISHER, HIMAL RATHNAKUMARA Abstract. We show that if f is a diffeomorphism of a manifold to itself, Λ is a mixing (or transitive) hyperbolic set, and V is a neighborhood
More informationMATH 23b, SPRING 2005 THEORETICAL LINEAR ALGEBRA AND MULTIVARIABLE CALCULUS Midterm (part 1) Solutions March 21, 2005
MATH 23b, SPRING 2005 THEORETICAL LINEAR ALGEBRA AND MULTIVARIABLE CALCULUS Midterm (part 1) Solutions March 21, 2005 1. True or False (22 points, 2 each) T or F Every set in R n is either open or closed
More informationAlgebraic Techniques for Switching Systems
Algebraic Techniques for Switching Systems And applications to Wireless Control Neworks Raphaël Jungers (UCL, Belgium) UCLA May 2014 Outline Joint spectral characteristics Application: WCNs and switching
More informationA genus 2 characterisation of translation surfaces with the lattice property
A genus 2 characterisation of translation surfaces with the lattice property (joint work with Howard Masur) 0-1 Outline 1. Translation surface 2. Translation flows 3. SL(2,R) action 4. Veech surfaces 5.
More informationRUSSELL JAY HENDEL. F i+1.
RUSSELL JAY HENDEL Abstract. Kimberling defines the function κ(n) = n 2 α n nα, and presents conjectures and open problems. We present three main theorems. The theorems provide quick, effectively computable,
More information7 Planar systems of linear ODE
7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution
More informationLOG-LINEAR ODES AND APPLICATIONS TO THE RICCI FLOW FOR HOMOGENEOUS SPACES
LOG-LINEAR ODES AND APPLICATIONS TO THE RICCI FLOW FOR HOMOGENEOUS SPACES MICHAEL KREISEL AND MARY RUSSELL Introduction Background Given a Riemannian manifold (M, g), recall the definition of Ricci curvature
More informationMathematical Statistical Physics Solution Sketch of the math part of the Exam
Prof. B. Paredes, Prof. P. Pickl Summer Term 2016 Lukas Nickel, Carlos Velasco July 22, 2016 Mathematical Statistical Physics Solution Sketch of the math part of the Exam Family Name: First name: Student
More informationA Topological Theory of Stirring
A Topological Theory of Stirring Jean-Luc Thiffeault Department of Mathematics Imperial College London University of Wisconsin, 15 December 2006 Collaborators: Matthew Finn Emmanuelle Gouillart Olivier
More informationCSL361 Problem set 4: Basic linear algebra
CSL361 Problem set 4: Basic linear algebra February 21, 2017 [Note:] If the numerical matrix computations turn out to be tedious, you may use the function rref in Matlab. 1 Row-reduced echelon matrices
More informationa b c d e GOOD LUCK! 3. a b c d e 12. a b c d e 4. a b c d e 13. a b c d e 5. a b c d e 14. a b c d e 6. a b c d e 15. a b c d e
MA Elem. Calculus Fall 07 Exam 07-09- Name: Sec.: Do not remove this answer page you will turn in the entire exam. No books or notes may be used. You may use an ACT-approved calculator during the exam,
More informationRigorous approximation of invariant measures for IFS Joint work April 8, with 2016 S. 1 Galat / 21
Rigorous approximation of invariant measures for IFS Joint work with S. Galatolo e I. Nisoli Maurizio Monge maurizio.monge@im.ufrj.br Universidade Federal do Rio de Janeiro April 8, 2016 Rigorous approximation
More informationAliprantis, Border: Infinite-dimensional Analysis A Hitchhiker s Guide
aliprantis.tex May 10, 2011 Aliprantis, Border: Infinite-dimensional Analysis A Hitchhiker s Guide Notes from [AB2]. 1 Odds and Ends 2 Topology 2.1 Topological spaces Example. (2.2) A semimetric = triangle
More information6.1 Polynomial Functions
6.1 Polynomial Functions Definition. A polynomial function is any function p(x) of the form p(x) = p n x n + p n 1 x n 1 + + p 2 x 2 + p 1 x + p 0 where all of the exponents are non-negative integers and
More information