Dynamical upper bounds in quantum mechanics. Laurent Marin. Helsinki, November 2008

Size: px
Start display at page:

Download "Dynamical upper bounds in quantum mechanics. Laurent Marin. Helsinki, November 2008"

Transcription

1 Dynamical upper bounds in quantum mechanics Laurent Marin Helsinki, November 2008 ynamical upper bounds in quantum mechanics (Helsinki, November 2008) Laurent Marin 1 / 30

2 Let H be a discrete self-adjoint Schrödinger operator on separable Hilbert space l 2 (Z) into himself. [Hψ](n) = ψ(n + 1) + ψ(n 1) + V (n)ψ(n) V is a function Z R, the so-called potential. ynamical upper bounds in quantum mechanics (Helsinki, November 2008) Laurent Marin 2 / 30

3 In this work, we consider only Sturmian potential : V (n) = ( (n + 1)β nβ )V V is a positive constant called strength of potential or couplage. β is an irrationnal number. Link to physics : Modeling electron/phonon spectrum in a quasicrystal. Only 2 values possible for the potential. Quasiperiodical potential : no period but any sequence repeats infinitly many times ynamical upper bounds in quantum mechanics (Helsinki, November 2008) Laurent Marin 3 / 30

4 Potential construction : By graphic means : line with β slope, iterated rotation by angle β. By word concatenation : The potential coincide with sturmian words, W 1 = V W 0 = 0 W k+1 = W a k+1 k W k 1, k 0. ynamical upper bounds in quantum mechanics (Helsinki, November 2008) Laurent Marin 4 / 30

5 Unitary dynamical evolution : ψ(t) = e ith ψ(0). In this case, ψ is known to spread out with time over the canonical basis of l 2 (Z). Dynamical field consist to quantify how fast it spreads. ynamical upper bounds in quantum mechanics (Helsinki, November 2008) Laurent Marin 5 / 30

6 The probability for the system to be in n over the time T (in average) : a(n, T ) = 2 T 0 e 2t/T e ith ψ(0), δ n 2 dt. We denote the time average outside probabilities P(N, T ) = a(n, T ), n >N ynamical upper bounds in quantum mechanics (Helsinki, November 2008) Laurent Marin 6 / 30

7 For p > 0, one defines moments of order p X p ψ(0) = n n p a(n, T ), And their growth exponents : β ψ(0) (p) = lim inf T log X p ψ(0) p log T log X p β + ψ(0) ψ(0) (p) = lim sup T p log T. ynamical upper bounds in quantum mechanics (Helsinki, November 2008) Laurent Marin 7 / 30

8 For all α [0, + ], S (α) = lim inf T log P(T α 2, T ) log T and S + (α) = lim sup T log P(T α 2, T ) log T The following critical exponents are particular of interest : α ± l = sup{α 0 : S ± (α) = 0}, α ± u = sup{α 0 : S ± (α) < }. ynamical upper bounds in quantum mechanics (Helsinki, November 2008) Laurent Marin 8 / 30

9 α l ± is the (lower and upper) rates of propagation of the essential part of the wavepacket and α u ± as the rates of propagation of the fastest part. In particular, if γ > α u + then P(T γ, T ) goes to 0 faster than any inverse power of T. They verify 0 α ± l α ± u. lim p 0 β± ψ (p) = α± l lim p β± ψ (p) = α± u. ynamical upper bounds in quantum mechanics (Helsinki, November 2008) Laurent Marin 9 / 30

10 (B. Simon) For this model, α ± u 1. Until 2005, it was an open question, if one can improve this so-called balistic bound. In this work, we improve upper bounds for the fast part of the wavepacket under assumption on the irrational number. ynamical upper bounds in quantum mechanics (Helsinki, November 2008) Laurent Marin 10 / 30

11 q 1 = 0, q 0 = 1, q k+1 = a k+1 q k + q k 1 Let β be an irrational number and H β with a sturmian log q potential associated to β. V > 20. If D = lim sup k k k is finite then α u + 2D log ( ). V 8 3 ynamical upper bounds in quantum mechanics (Helsinki, November 2008) Laurent Marin 11 / 30

12 ψ(n + 1) + ψ(n 1) + V (n)ψ(n) = zψ(n) (1) with z C and ψ a non-zero vector. One can reformulate this in terms of transfer matrix. ( ) ( ) ψ(n + 1) ψ(1) = F (n, z) ψ(n) ψ(0) T (n, z)...t (1, z) n 1, F (n, z) = Id n = 0, [T (n, z)] 1...[T (0, z)] 1 n 1. ( ) z V (m) 1 T (m, z) = 1 0 ynamical upper bounds in quantum mechanics (Helsinki, November 2008) Laurent Marin 12 / 30

13 The transfer matrix play a key role in all the theory. Of course, in eigenvalue/vector study, and spectrum study. But also in dynamical field, with this estimate : (Damanik, Tcheremtchansev) P(N, T ) exp( cn) ( K + T 3 K max 1 k N F (k, E + i 2) 1 T ) de, ynamical upper bounds in quantum mechanics (Helsinki, November 2008) Laurent Marin 13 / 30

14 p 1 = 1, p 0 = 0, q 1 = 0, q 0 = 1, p k+1 = a k+1 p k + p k 1 We denote q k+1 = a k+1 q k + q k 1 M k (z) = F (q k, z) M k+1 (z) = M k 1 (z)m k (z) a k+1 ynamical upper bounds in quantum mechanics (Helsinki, November 2008) Laurent Marin 14 / 30

15 Let t k,p be the trace of the matrix M k 1 M p k. The evolution along the p index is given by and consequently, t k,p+1 = t k+1,0 t k,p t k,p 1, t k,p+1 = S p (t k+1,0 )t k,1 S p 1 (t k+1,0 )t k,0 S p denotes the p th Tchebychev polynomial of the second kind The evolution along the k index is related to the p-evolution by t k+1,0 = t k,ak+1, ynamical upper bounds in quantum mechanics (Helsinki, November 2008) Laurent Marin 15 / 30

16 The behavior of the sequence of trace t k,0 is interesting in dynamics but also for the spectrum of the operator H. σ = {E R, s.t.{t k,0 (E)} k bounded} We can be more precise and tell with the next criteria when this sequence is bounded or not. ynamical upper bounds in quantum mechanics (Helsinki, November 2008) Laurent Marin 16 / 30

17 A necessary and sufficient condition that {t k,0 (E)} k be unbounded is that t N+1,0 (E) > 2, t N,0 (E) > 2, t N 1,0 (E) 2 for some N 0. This N is unique. Denote G k = G k 1 + a k G k 2, G 0 = 1, G 1 = 1. We have t k,0 (E) e cg k N k > N. ynamical upper bounds in quantum mechanics (Helsinki, November 2008) Laurent Marin 17 / 30

18 Denote σ k,p = {E R, t k,p (E) 2}. (i) the set σ k,p is made of pq k + q k 1 distinct intervals, (ii) the σ k,0 coincide with the spectrum of periodic operator H k replacing β with p k /q k in definition) (iii) σ σ k+1,0 σ k,0 Moreover, σ the spectrum of H verify σ = k=n (σ k+1,0 σ k,0 ) ynamical upper bounds in quantum mechanics (Helsinki, November 2008) Laurent Marin 18 / 30

19 For a given k, we call type I gap : a band of σ k,1 included in a band of σ k,0 and therefore in a gap of σ k+1,0, type II band : a band of σ k+1,0 included in a band of σ k, 1 and in a gap of σ k,0, type III band : a band of σ k+1,0 included in a band of σ k,0 and in a gap of σ k,1. ynamical upper bounds in quantum mechanics (Helsinki, November 2008) Laurent Marin 19 / 30

20 (Raymond) At a given level k, (i) a type I gap contains an unique type II band of σ k+2,0. (ii) a type II band contains (a k+1 + 1) bands of type I of σ k+1,1. They are alternated with (a k+1 ) type III bands of σ k+2,0 (iii) a type III band contains (a k+1 ) bands of type I of σ k+1,1. They are alternated with (a k+1 1) type III bands of σ k+2,0 ynamical upper bounds in quantum mechanics (Helsinki, November 2008) Laurent Marin 20 / 30

21 (Liu, Wen) Each band length has a lower bound with coefficient of Q n and an upper bound with coefficient of P n. with c 1 = 3 V 8 Q n = with c 2 = 1 V +5. P n = 0 ca n c 1 /a n 0 c 1 /a n c 1 /a n 0 c 1 /a n 0 c a n c 2 (a n + 2) 3 0 c 2 (a n + 2) 3 c 2 (a n + 2) 3 0 c 2 (a n + 2) 3 ynamical upper bounds in quantum mechanics (Helsinki, November 2008) Laurent Marin 21 / 30

22 By now, we define the periodic approximants spectrum not only in R but in C. σ k,0 = {z C : t k,0 (z) 2} If k 3, and V > 20 then there exist constants c,d > 0 such that q k 1 j=1 q B(x (j) k, r k 1 k) σ k,0 B(x (j) k, R k) j=1 where {x (j) k } 1 j q k 1 are the zeros of t k,0, r k = c. min band length and R k = d.max band length. ynamical upper bounds in quantum mechanics (Helsinki, November 2008) Laurent Marin 22 / 30

23 As x (j) k are real, we have σ k,0 {z C : Imz < R k } {z C : Imz < dq γ(v ) k }. for a suitable γ(v ). This implies σ σ k,0 σ k+1,0 {z C : Imz < dq γ(v ) k }. (2) We should have R k < dq γ(v ) k chosen by taking : so a suitable γ can be γ(v ) lim sup k log c 1. k 2 log q k ynamical upper bounds in quantum mechanics (Helsinki, November 2008) Laurent Marin 23 / 30

24 For T > 1, denote by k(t ) the unique integer with and let q γ(v ) k(t ) 1 d T qγ(v ) k(t ) d N(T ) = q k(t )+ k(t ). N(T ) C ν T 1 γ(v ) T ν. (3) ynamical upper bounds in quantum mechanics (Helsinki, November 2008) Laurent Marin 24 / 30

25 ( K P(N(T ), T ) exp( cn(t )) + T 3 K exp( cn(t )) + T 3 e 2cG k(t ) max 1 q n N(T ) From this bound, we see that P(N(T ), T ) goes to 0 faster than any inverse power of T. As N(T ) C ν T 1 γ(v ) T ν. M n(e we have α + u 1 γ(v ) + ν ynamical upper bounds in quantum mechanics (Helsinki, November 2008) Laurent Marin 25 / 30

26 Let β be an irrational number and H β with a sturmian log q potential associated to β. V > 20. If D = lim sup k k k is finite then α u + 2D log ( ). V 8 3 Moreover, for an irrational with continued fraction expansion containing no 1, the dynamical upper bound becomes α u + D log ( ). V 8 3 ynamical upper bounds in quantum mechanics (Helsinki, November 2008) Laurent Marin 26 / 30

27 It is clear that taking V large enough, one can obtain a non trivial bound that is smaller than 1. The condition D finite is true for almost all irrational. Moreover, (Khintchin) For almost all β with respect to Lebesgue measure, D = lim sup k log q k k = D K = π2 12 log 2. ynamical upper bounds in quantum mechanics (Helsinki, November 2008) Laurent Marin 27 / 30

28 There exist an irrational number β with D = + such that for any V > 20. α + u = 1. Proof : Sturmian potentials of operator H β and H βn same first p n+1 + q n+1 values. have the ynamical upper bounds in quantum mechanics (Helsinki, November 2008) Laurent Marin 28 / 30

29 (Last) Let H 1 = + V 1 and H 2 = + V 2 on l 2 (Z), and such that V 1 (k), V 2 (k) < C for all k Z and C a constant. Let T > 0 and ε > 0 be fixed so there exist L(T, ε), δ > 0 such that V 1 (k) V 2 (k) < δ for all k < L, then X 2 H 1,ψ(0) X 2 H 2,ψ(0) < ε. ynamical upper bounds in quantum mechanics (Helsinki, November 2008) Laurent Marin 29 / 30

30 As H βn is a periodic potential operator, one has X 2 β n T > C n T 2 choose T n big enough such that C n > 1 log T n. One can then choose a n+1 such that L(T n, ε) p n+1 + q n+1. Inductively, we have a sequence T n going to infinity and an irrational number β with X 2 β Tn > T 2 n log T n ε > T 2 δ n, δ > 0 α + u β δ 1 (2) > 1 δ, δ > 0. ynamical upper bounds in quantum mechanics (Helsinki, November 2008) Laurent Marin 30 / 30

Dynamics of SL(2, R)-Cocycles and Applications to Spectral Theory

Dynamics of SL(2, R)-Cocycles and Applications to Spectral Theory Dynamics of SL(2, R)-Cocycles and Applications to Spectral Theory David Damanik Rice University Global Dynamics Beyond Uniform Hyperbolicity Peking University Beijing International Center of Mathematics

More information

Root Finding For NonLinear Equations Bisection Method

Root Finding For NonLinear Equations Bisection Method Root Finding For NonLinear Equations Bisection Method P. Sam Johnson November 17, 2014 P. Sam Johnson (NITK) Root Finding For NonLinear Equations Bisection MethodNovember 17, 2014 1 / 26 Introduction The

More information

THE FRACTAL DIMENSION OF THE SPECTRUM OF THE FIBONACCI HAMILTONIAN

THE FRACTAL DIMENSION OF THE SPECTRUM OF THE FIBONACCI HAMILTONIAN THE FRACTAL DIMENSION OF THE SPECTRUM OF THE FIBONACCI HAMILTONIAN DAVID DAMANIK, MARK EMBREE, ANTON GORODETSKI, AND SERGUEI TCHEREMCHANTSEV Abstract. We study the spectrum of the Fibonacci Hamiltonian

More information

Palindromic complexity of infinite words associated with simple Parry numbers

Palindromic complexity of infinite words associated with simple Parry numbers Palindromic complexity of infinite words associated with simple Parry numbers Petr Ambrož (1)(2) Christiane Frougny (2)(3) Zuzana Masáková (1) Edita Pelantová (1) March 22, 2006 (1) Doppler Institute for

More information

1.2. Direction Fields: Graphical Representation of the ODE and its Solution Let us consider a first order differential equation of the form dy

1.2. Direction Fields: Graphical Representation of the ODE and its Solution Let us consider a first order differential equation of the form dy .. Direction Fields: Graphical Representation of the ODE and its Solution Let us consider a first order differential equation of the form dy = f(x, y). In this section we aim to understand the solution

More information

RESOLVENT METHODS FOR QUANTUM WALKS WITH AN APPLICATION TO A THUE-MORSE QUANTUM WALK

RESOLVENT METHODS FOR QUANTUM WALKS WITH AN APPLICATION TO A THUE-MORSE QUANTUM WALK RESOLVENT METHODS FOR QUANTUM WALKS WITH AN APPLICATION TO A THUE-MORSE QUANTUM WALK JAKE FILLMAN Abstract. In this expository work, we discuss spatially inhomogeneous quantum walks in one dimension and

More information

Finite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product

Finite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )

More information

Approximation exponents for algebraic functions in positive characteristic

Approximation exponents for algebraic functions in positive characteristic ACTA ARITHMETICA LX.4 (1992) Approximation exponents for algebraic functions in positive characteristic by Bernard de Mathan (Talence) In this paper, we study rational approximations for algebraic functions

More information

Transfer matrices and transport for Schrödinger operators

Transfer matrices and transport for Schrödinger operators Transfer matrices and transport for Schrödinger operators François Germinet UMR 8524 CNRS, UFR de Mathématiques, Université de Lille 1, F-59655 Villeneuve d Ascq Cédex, France E-mail: germinet@agat.univ-lille1.fr

More information

Semigroup Growth Bounds

Semigroup Growth Bounds Semigroup Growth Bounds First Meeting on Asymptotics of Operator Semigroups E.B. Davies King s College London Oxford, September 2009 E.B. Davies (KCL) Semigroup Growth Bounds Oxford, September 2009 1 /

More information

Quasiperiodic Schrodinger operators: sharp arithmetic spectral transitions and universal hierarchical structure of eigenfunctions

Quasiperiodic Schrodinger operators: sharp arithmetic spectral transitions and universal hierarchical structure of eigenfunctions Quasiperiodic Schrodinger operators: sharp arithmetic spectral transitions and universal hierarchical structure of eigenfunctions S. Jitomirskaya Atlanta, October 10, 2016 Almost Mathieu operators (H λ,α,θ

More information

Heat kernels of some Schrödinger operators

Heat kernels of some Schrödinger operators Heat kernels of some Schrödinger operators Alexander Grigor yan Tsinghua University 28 September 2016 Consider an elliptic Schrödinger operator H = Δ + Φ, where Δ = n 2 i=1 is the Laplace operator in R

More information

A Smooth Operator, Operated Correctly

A Smooth Operator, Operated Correctly Clark Department of Mathematics Bard College at Simon s Rock March 6, 2014 Abstract By looking at familiar smooth functions in new ways, we can make sense of matrix-valued infinite series and operator-valued

More information

Recitation: 10 11/06/03

Recitation: 10 11/06/03 Recitation: 10 11/06/03 Ensembles and Relation to T.D. It is possible to expand both sides of the equation with F = kt lnq Q = e βe i If we expand both sides of this equation, we apparently obtain: i F

More information

Combinatorics On Sturmian Words. Amy Glen. Major Review Seminar. February 27, 2004 DISCIPLINE OF PURE MATHEMATICS

Combinatorics On Sturmian Words. Amy Glen. Major Review Seminar. February 27, 2004 DISCIPLINE OF PURE MATHEMATICS Combinatorics On Sturmian Words Amy Glen Major Review Seminar February 27, 2004 DISCIPLINE OF PURE MATHEMATICS OUTLINE Finite and Infinite Words Sturmian Words Mechanical Words and Cutting Sequences Infinite

More information

Notes on uniform convergence

Notes on uniform convergence Notes on uniform convergence Erik Wahlén erik.wahlen@math.lu.se January 17, 2012 1 Numerical sequences We begin by recalling some properties of numerical sequences. By a numerical sequence we simply mean

More information

The Sommerfeld Polynomial Method: Harmonic Oscillator Example

The Sommerfeld Polynomial Method: Harmonic Oscillator Example Chemistry 460 Fall 2017 Dr. Jean M. Standard October 2, 2017 The Sommerfeld Polynomial Method: Harmonic Oscillator Example Scaling the Harmonic Oscillator Equation Recall the basic definitions of the harmonic

More information

MATRIX INTEGRALS AND MAP ENUMERATION 2

MATRIX INTEGRALS AND MAP ENUMERATION 2 MATRIX ITEGRALS AD MAP EUMERATIO 2 IVA CORWI Abstract. We prove the generating function formula for one face maps and for plane diagrams using techniques from Random Matrix Theory and orthogonal polynomials.

More information

FUNCTIONAL ANALYSIS LECTURE NOTES: WEAK AND WEAK* CONVERGENCE

FUNCTIONAL ANALYSIS LECTURE NOTES: WEAK AND WEAK* CONVERGENCE FUNCTIONAL ANALYSIS LECTURE NOTES: WEAK AND WEAK* CONVERGENCE CHRISTOPHER HEIL 1. Weak and Weak* Convergence of Vectors Definition 1.1. Let X be a normed linear space, and let x n, x X. a. We say that

More information

DAVID DAMANIK AND DANIEL LENZ. Dedicated to Barry Simon on the occasion of his 60th birthday.

DAVID DAMANIK AND DANIEL LENZ. Dedicated to Barry Simon on the occasion of his 60th birthday. UNIFORM SZEGŐ COCYCLES OVER STRICTLY ERGODIC SUBSHIFTS arxiv:math/052033v [math.sp] Dec 2005 DAVID DAMANIK AND DANIEL LENZ Dedicated to Barry Simon on the occasion of his 60th birthday. Abstract. We consider

More information

Quantum Mechanics for Mathematicians: Energy, Momentum, and the Quantum Free Particle

Quantum Mechanics for Mathematicians: Energy, Momentum, and the Quantum Free Particle Quantum Mechanics for Mathematicians: Energy, Momentum, and the Quantum Free Particle Peter Woit Department of Mathematics, Columbia University woit@math.columbia.edu November 28, 2012 We ll now turn to

More information

Page 684. Lecture 40: Coordinate Transformations: Time Transformations Date Revised: 2009/02/02 Date Given: 2009/02/02

Page 684. Lecture 40: Coordinate Transformations: Time Transformations Date Revised: 2009/02/02 Date Given: 2009/02/02 Page 684 Lecture 40: Coordinate Transformations: Time Transformations Date Revised: 2009/02/02 Date Given: 2009/02/02 Time Transformations Section 12.5 Symmetries: Time Transformations Page 685 Time Translation

More information

= 1 2 x (x 1) + 1 {x} (1 {x}). [t] dt = 1 x (x 1) + O (1), [t] dt = 1 2 x2 + O (x), (where the error is not now zero when x is an integer.

= 1 2 x (x 1) + 1 {x} (1 {x}). [t] dt = 1 x (x 1) + O (1), [t] dt = 1 2 x2 + O (x), (where the error is not now zero when x is an integer. Problem Sheet,. i) Draw the graphs for [] and {}. ii) Show that for α R, α+ α [t] dt = α and α+ α {t} dt =. Hint Split these integrals at the integer which must lie in any interval of length, such as [α,

More information

Consider a system of n ODEs. parameter t, periodic with period T, Let Φ t to be the fundamental matrix of this system, satisfying the following:

Consider a system of n ODEs. parameter t, periodic with period T, Let Φ t to be the fundamental matrix of this system, satisfying the following: Consider a system of n ODEs d ψ t = A t ψ t, dt ψ: R M n 1 C where A: R M n n C is a continuous, (n n) matrix-valued function of real parameter t, periodic with period T, t R k Z A t + kt = A t. Let Φ

More information

Kolmogorov Complexity and Diophantine Approximation

Kolmogorov Complexity and Diophantine Approximation Kolmogorov Complexity and Diophantine Approximation Jan Reimann Institut für Informatik Universität Heidelberg Kolmogorov Complexity and Diophantine Approximation p. 1/24 Algorithmic Information Theory

More information

INITIAL POWERS OF STURMIAN SEQUENCES

INITIAL POWERS OF STURMIAN SEQUENCES INITIAL POWERS OF STURMIAN SEQUENCES VALÉRIE BERTHÉ, CHARLES HOLTON, AND LUCA Q. ZAMBONI Abstract. We investigate powers of prefixes in Sturmian sequences. We obtain an explicit formula for ice(ω), the

More information

11. Recursion Method: Concepts

11. Recursion Method: Concepts University of Rhode Island DigitalCommons@URI Nonequilibrium Statistical Physics Physics Course Materials 16 11. Recursion Method: Concepts Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative

More information

1 Math 241A-B Homework Problem List for F2015 and W2016

1 Math 241A-B Homework Problem List for F2015 and W2016 1 Math 241A-B Homework Problem List for F2015 W2016 1.1 Homework 1. Due Wednesday, October 7, 2015 Notation 1.1 Let U be any set, g be a positive function on U, Y be a normed space. For any f : U Y let

More information

arxiv: v1 [math.sp] 25 Jun 2012

arxiv: v1 [math.sp] 25 Jun 2012 THE DENSITY OF STATES MEASURE OF THE WEAKLY COUPLED FIBONACCI HAMILTONIAN DAVID DAMANIK AND ANTON GORODETSKI arxiv:06.60v [math.sp] Jun 0 Abstract. We consider the density of states measure of the Fibonacci

More information

Fine scales of decay rates of operator semigroups

Fine scales of decay rates of operator semigroups Operator Semigroups meet everything else Herrnhut, 5 June 2013 A damped wave equation 2 u u + 2a(x) u t2 t = 0 (t > 0, x Ω) u(x, t) = 0 (t > 0, x Ω) u(, 0) = u 0 H0 1 (Ω), u t (, 0) = u 1 L 2 (Ω). Here,

More information

Asymptotic Spectral Properties of the Schrödinger Operator with Thue-Morse Potential

Asymptotic Spectral Properties of the Schrödinger Operator with Thue-Morse Potential Asymptotic Spectral Properties of the Schrödinger Operator with Thue-Morse Potential William Clark Rachael Kline Michaela Stone Cornell SMI Advisors: May Mei and Andrew Zemke August 2, 2013 Introduction

More information

Dirichlet uniformly well-approximated numbers

Dirichlet uniformly well-approximated numbers Dirichlet uniformly well-approximated numbers Lingmin LIAO (joint work with Dong Han Kim) Université Paris-Est Créteil (University Paris 12) Approximation and numeration Université Paris Diderot-Paris

More information

Euler s Method and Logistic Growth (BC Only)

Euler s Method and Logistic Growth (BC Only) Euler s Method Students should be able to: Approximate numerical solutions of differential equations using Euler s method without a calculator. Recognize the method as a recursion formula extension of

More information

1D-Quantum Systems. a review ( ) Sponsoring. Jean BELLISSARD. NSF grant No

1D-Quantum Systems. a review ( ) Sponsoring. Jean BELLISSARD. NSF grant No 1D-Quantum Systems Sponsoring (1980-1993) a review NSF grant No. 0901514 Jean BELLISSARD Georgia Institute of Technology, Atlanta School of Mathematics & School of Physics e-mail: jeanbel@math.gatech.edu

More information

dynamical Diophantine approximation

dynamical Diophantine approximation Dioph. Appro. Dynamical Dioph. Appro. in dynamical Diophantine approximation WANG Bao-Wei Huazhong University of Science and Technology Joint with Zhang Guo-Hua Central China Normal University 24-28 July

More information

Report on some recent advances in Diophantine approximation

Report on some recent advances in Diophantine approximation Report on some recent advances in Diophantine approximation Michel Waldschmidt Université Pierre et Marie Curie Paris 6, UMR 7586 IMJ Institut de Mathématiques de Jussieu, 4, Place Jussieu, Paris, F 75005

More information

6 Linear Equation. 6.1 Equation with constant coefficients

6 Linear Equation. 6.1 Equation with constant coefficients 6 Linear Equation 6.1 Equation with constant coefficients Consider the equation ẋ = Ax, x R n. This equating has n independent solutions. If the eigenvalues are distinct then the solutions are c k e λ

More information

General Power Series

General Power Series General Power Series James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University March 29, 2018 Outline Power Series Consequences With all these preliminaries

More information

Assignment 4. u n+1 n(n + 1) i(i + 1) = n n (n + 1)(n + 2) n(n + 2) + 1 = (n + 1)(n + 2) 2 n + 1. u n (n + 1)(n + 2) n(n + 1) = n

Assignment 4. u n+1 n(n + 1) i(i + 1) = n n (n + 1)(n + 2) n(n + 2) + 1 = (n + 1)(n + 2) 2 n + 1. u n (n + 1)(n + 2) n(n + 1) = n Assignment 4 Arfken 5..2 We have the sum Note that the first 4 partial sums are n n(n + ) s 2, s 2 2 3, s 3 3 4, s 4 4 5 so we guess that s n n/(n + ). Proving this by induction, we see it is true for

More information

RESULTS ON THE BETA FUNCTION. B(x, y) = for x, y > 0, see for example Sneddon [4]. It then follows that

RESULTS ON THE BETA FUNCTION. B(x, y) = for x, y > 0, see for example Sneddon [4]. It then follows that SARAJEVO JOURNAL OF MATHEMATICS Vol.9 (21) (2013), 101 108 DOI: 10.5644/SJM.09.1.09 RESULTS ON THE BETA FUNCTION BRIAN FISHER 1. Introduction The Beta function is usually defined by B(x, y) = 1 0 t x 1

More information

Factorizations of the Fibonacci Infinite Word

Factorizations of the Fibonacci Infinite Word 2 3 47 6 23 Journal of Integer Sequences, Vol. 8 (205), Article 5.9.3 Factorizations of the Fibonacci Infinite Word Gabriele Fici Dipartimento di Matematica e Informatica Università di Palermo Via Archirafi

More information

QUARTIC POWER SERIES IN F 3 ((T 1 )) WITH BOUNDED PARTIAL QUOTIENTS. Alain Lasjaunias

QUARTIC POWER SERIES IN F 3 ((T 1 )) WITH BOUNDED PARTIAL QUOTIENTS. Alain Lasjaunias QUARTIC POWER SERIES IN F 3 ((T 1 )) WITH BOUNDED PARTIAL QUOTIENTS Alain Lasjaunias 1991 Mathematics Subject Classification: 11J61, 11J70. 1. Introduction. We are concerned with diophantine approximation

More information

On approximation of real, complex, and p-adic numbers by algebraic numbers of bounded degree. by K. I. Tsishchanka

On approximation of real, complex, and p-adic numbers by algebraic numbers of bounded degree. by K. I. Tsishchanka On approximation of real, complex, and p-adic numbers by algebraic numbers of bounded degree by K. I. Tsishchanka I. On approximation by rational numbers Theorem 1 (Dirichlet, 1842). For any real irrational

More information

A Glimpse of Quantum Computation

A Glimpse of Quantum Computation A Glimpse of Quantum Computation Zhengfeng Ji (UTS:QSI) QCSS 2018, UTS 1. 1 Introduction What is quantum computation? Where does the power come from? Superposition Incompatible states can coexist Transformation

More information

THE DENSITY OF STATES MEASURE OF THE WEAKLY COUPLED FIBONACCI HAMILTONIAN

THE DENSITY OF STATES MEASURE OF THE WEAKLY COUPLED FIBONACCI HAMILTONIAN THE DENSITY OF STATES MEASURE OF THE WEAKLY COUPLED FIBONACCI HAMILTONIAN DAVID DAMANIK AND ANTON GORODETSKI Abstract. We consider the density of states measure of the Fibonacci Hamiltonian and show that,

More information

Real Analysis Math 131AH Rudin, Chapter #1. Dominique Abdi

Real Analysis Math 131AH Rudin, Chapter #1. Dominique Abdi Real Analysis Math 3AH Rudin, Chapter # Dominique Abdi.. If r is rational (r 0) and x is irrational, prove that r + x and rx are irrational. Solution. Assume the contrary, that r+x and rx are rational.

More information

Gaussian Random Fields

Gaussian Random Fields Gaussian Random Fields Mini-Course by Prof. Voijkan Jaksic Vincent Larochelle, Alexandre Tomberg May 9, 009 Review Defnition.. Let, F, P ) be a probability space. Random variables {X,..., X n } are called

More information

Uniform Diophantine approximation related to b-ary and β-expansions

Uniform Diophantine approximation related to b-ary and β-expansions Uniform Diophantine approximation related to b-ary and β-expansions Lingmin LIAO (joint work with Yann Bugeaud) Université Paris-Est Créteil (University Paris 12) Universidade Federal de Bahia April 8th

More information

Quantum Mechanics I Physics 5701

Quantum Mechanics I Physics 5701 Quantum Mechanics I Physics 5701 Z. E. Meziani 02/23//2017 Physics 5701 Lecture Outline 1 General Formulation of Quantum Mechanics 2 Measurement of physical quantities and observables 3 Representations

More information

Math Solutions to homework 5

Math Solutions to homework 5 Math 75 - Solutions to homework 5 Cédric De Groote November 9, 207 Problem (7. in the book): Let {e n } be a complete orthonormal sequence in a Hilbert space H and let λ n C for n N. Show that there is

More information

Eigenvalues of Robin Laplacians on infinite sectors and application to polygons

Eigenvalues of Robin Laplacians on infinite sectors and application to polygons Eigenvalues of Robin Laplacians on infinite sectors and application to polygons Magda Khalile (joint work with Konstantin Pankrashkin) Université Paris Sud /25 Robin Laplacians on infinite sectors 1 /

More information

Title. Author(s)Arai, Asao. CitationLetters in Mathematical Physics, 85(1): Issue Date Doc URL. Rights. Type.

Title. Author(s)Arai, Asao. CitationLetters in Mathematical Physics, 85(1): Issue Date Doc URL. Rights. Type. Title On the Uniqueness of Weak Weyl Representations of th Author(s)Arai, Asao CitationLetters in Mathematical Physics, 85(1): 15-25 Issue Date 2008-07 Doc URL http://hdl.handle.net/2115/38135 Rights The

More information

Aperiodic Hamiltonians

Aperiodic Hamiltonians Spectrum Approximation for Aperiodic Hamiltonians Sponsoring Jean BELLISSARD Westfälische Wilhelms-Universität, Münster Department of Mathematics Georgia Institute of Technology, Atlanta School of Mathematics

More information

Brownian Motion. 1 Definition Brownian Motion Wiener measure... 3

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

to 1D Periodic Approximants Aperiodic Hamiltonians Sponsoring Jean BELLISSARD CRC 701, Bielefeld, Germany

to 1D Periodic Approximants Aperiodic Hamiltonians Sponsoring Jean BELLISSARD CRC 701, Bielefeld, Germany Sponsoring Periodic Approximants to 1D Aperiodic Hamiltonians CRC 701, Bielefeld, Germany Jean BELLISSARD Georgia Institute of Technology, Atlanta School of Mathematics & School of Physics e-mail: jeanbel@math.gatech.edu

More information

Introduction to Quantum Spin Systems

Introduction to Quantum Spin Systems 1 Introduction to Quantum Spin Systems Lecture 2 Sven Bachmann (standing in for Bruno Nachtergaele) Mathematics, UC Davis MAT290-25, CRN 30216, Winter 2011, 01/10/11 2 Basic Setup For concreteness, consider

More information

Chapter 3: Baire category and open mapping theorems

Chapter 3: Baire category and open mapping theorems MA3421 2016 17 Chapter 3: Baire category and open mapping theorems A number of the major results rely on completeness via the Baire category theorem. 3.1 The Baire category theorem 3.1.1 Definition. A

More information

On lower bounds for integration of multivariate permutation-invariant functions

On lower bounds for integration of multivariate permutation-invariant functions On lower bounds for of multivariate permutation-invariant functions Markus Weimar Philipps-University Marburg Oberwolfach October 2013 Research supported by Deutsche Forschungsgemeinschaft DFG (DA 360/19-1)

More information

The geometrical semantics of algebraic quantum mechanics

The geometrical semantics of algebraic quantum mechanics The geometrical semantics of algebraic quantum mechanics B. Zilber University of Oxford November 29, 2016 B.Zilber, The semantics of the canonical commutation relations arxiv.org/abs/1604.07745 The geometrical

More information

Fréchet differentiability of the norm of L p -spaces associated with arbitrary von Neumann algebras

Fréchet differentiability of the norm of L p -spaces associated with arbitrary von Neumann algebras Fréchet differentiability of the norm of L p -spaces associated with arbitrary von Neumann algebras (Part I) Fedor Sukochev (joint work with D. Potapov, A. Tomskova and D. Zanin) University of NSW, AUSTRALIA

More information

MAT 135B Midterm 1 Solutions

MAT 135B Midterm 1 Solutions MAT 35B Midterm Solutions Last Name (PRINT): First Name (PRINT): Student ID #: Section: Instructions:. Do not open your test until you are told to begin. 2. Use a pen to print your name in the spaces above.

More information

called the homomorphism induced by the inductive limit. One verifies that the diagram

called the homomorphism induced by the inductive limit. One verifies that the diagram Inductive limits of C -algebras 51 sequences {a n } suc tat a n, and a n 0. If A = A i for all i I, ten A i = C b (I,A) and i I A i = C 0 (I,A). i I 1.10 Inductive limits of C -algebras Definition 1.10.1

More information

Sobolev Spaces. Chapter Hölder spaces

Sobolev Spaces. Chapter Hölder spaces Chapter 2 Sobolev Spaces Sobolev spaces turn out often to be the proper setting in which to apply ideas of functional analysis to get information concerning partial differential equations. Here, we collect

More information

Fall TMA4145 Linear Methods. Solutions to exercise set 9. 1 Let X be a Hilbert space and T a bounded linear operator on X.

Fall TMA4145 Linear Methods. Solutions to exercise set 9. 1 Let X be a Hilbert space and T a bounded linear operator on X. TMA445 Linear Methods Fall 26 Norwegian University of Science and Technology Department of Mathematical Sciences Solutions to exercise set 9 Let X be a Hilbert space and T a bounded linear operator on

More information

Continued Fractions New and Old Results

Continued Fractions New and Old Results Continued Fractions New and Old Results Jeffrey Shallit School of Computer Science University of Waterloo Waterloo, Ontario N2L 3G1 Canada shallit@cs.uwaterloo.ca https://www.cs.uwaterloo.ca/~shallit Joint

More information

Infinite Continued Fractions

Infinite Continued Fractions Infinite Continued Fractions 8-5-200 The value of an infinite continued fraction [a 0 ; a, a 2, ] is lim c k, where c k is the k-th convergent k If [a 0 ; a, a 2, ] is an infinite continued fraction with

More information

8 Periodic Linear Di erential Equations - Floquet Theory

8 Periodic Linear Di erential Equations - Floquet Theory 8 Periodic Linear Di erential Equations - Floquet Theory The general theory of time varying linear di erential equations _x(t) = A(t)x(t) is still amazingly incomplete. Only for certain classes of functions

More information

Variational construction of periodic and connecting orbits in the planar Sitnikov problem. Mitsuru Shibayama(Kyoto University)

Variational construction of periodic and connecting orbits in the planar Sitnikov problem. Mitsuru Shibayama(Kyoto University) Variational construction of periodic and connecting orbits in the planar Sitnikov problem Mitsuru Shibayama(Kyoto University) 1 The three-body problem Consider the planar three-body problem which is governed

More information

Physics 106a, Caltech 13 November, Lecture 13: Action, Hamilton-Jacobi Theory. Action-Angle Variables

Physics 106a, Caltech 13 November, Lecture 13: Action, Hamilton-Jacobi Theory. Action-Angle Variables Physics 06a, Caltech 3 November, 08 Lecture 3: Action, Hamilton-Jacobi Theory Starred sections are advanced topics for interest and future reference. The unstarred material will not be tested on the final

More information

Real Variables: Solutions to Homework 9

Real Variables: Solutions to Homework 9 Real Variables: Solutions to Homework 9 Theodore D Drivas November, 20 xercise 0 Chapter 8, # : For complex-valued, measurable f, f = f + if 2 with f i real-valued and measurable, we have f = f + i f 2

More information

MATH 162. Midterm 2 ANSWERS November 18, 2005

MATH 162. Midterm 2 ANSWERS November 18, 2005 MATH 62 Midterm 2 ANSWERS November 8, 2005. (0 points) Does the following integral converge or diverge? To get full credit, you must justify your answer. 3x 2 x 3 + 4x 2 + 2x + 4 dx You may not be able

More information

2. As we shall see, we choose to write in terms of σ x because ( X ) 2 = σ 2 x.

2. As we shall see, we choose to write in terms of σ x because ( X ) 2 = σ 2 x. Section 5.1 Simple One-Dimensional Problems: The Free Particle Page 9 The Free Particle Gaussian Wave Packets The Gaussian wave packet initial state is one of the few states for which both the { x } and

More information

Quantum mechanics. Chapter The quantum mechanical formalism

Quantum mechanics. Chapter The quantum mechanical formalism Chapter 5 Quantum mechanics 5.1 The quantum mechanical formalism The realisation, by Heisenberg, that the position and momentum of a quantum mechanical particle cannot be measured simultaneously renders

More information

Quantum Field Theory II

Quantum Field Theory II Quantum Field Theory II T. Nguyen PHY 391 Independent Study Term Paper Prof. S.G. Rajeev University of Rochester April 2, 218 1 Introduction The purpose of this independent study is to familiarize ourselves

More information

arxiv: v1 [math.sp] 14 Jan 2010

arxiv: v1 [math.sp] 14 Jan 2010 SPECTRAL AND QUANTUM DYNAMICAL PROPERTIES OF THE WEAKLY COUPLED FIBONACCI HAMILTONIAN DAID DAMANIK AND ANTON GORODETSKI arxiv:1001.2552v1 [math.sp] 14 Jan 2010 Abstract. We consider the spectrum of the

More information

Introduction to the GRAPE Algorithm

Introduction to the GRAPE Algorithm June 8, 2010 Reference: J. Mag. Res. 172, 296 (2005) Journal of Magnetic Resonance 172 (2005) 296 305 www.elsevier.com/locate/jmr Optimal control of coupled spin dynamics: design of NMR pulse sequences

More information

Criterions on periodic feedback stabilization for some evolution equations

Criterions on periodic feedback stabilization for some evolution equations Criterions on periodic feedback stabilization for some evolution equations School of Mathematics and Statistics, Wuhan University, P. R. China (Joint work with Yashan Xu, Fudan University) Toulouse, June,

More information

MATH 117 LECTURE NOTES

MATH 117 LECTURE NOTES MATH 117 LECTURE NOTES XIN ZHOU Abstract. This is the set of lecture notes for Math 117 during Fall quarter of 2017 at UC Santa Barbara. The lectures follow closely the textbook [1]. Contents 1. The set

More information

Solutions Final Exam May. 14, 2014

Solutions Final Exam May. 14, 2014 Solutions Final Exam May. 14, 2014 1. (a) (10 points) State the formal definition of a Cauchy sequence of real numbers. A sequence, {a n } n N, of real numbers, is Cauchy if and only if for every ɛ > 0,

More information

Errata Applied Analysis

Errata Applied Analysis Errata Applied Analysis p. 9: line 2 from the bottom: 2 instead of 2. p. 10: Last sentence should read: The lim sup of a sequence whose terms are bounded from above is finite or, and the lim inf of a sequence

More information

Numerical Simulation of Spin Dynamics

Numerical Simulation of Spin Dynamics Numerical Simulation of Spin Dynamics Marie Kubinova MATH 789R: Advanced Numerical Linear Algebra Methods with Applications November 18, 2014 Introduction Discretization in time Computing the subpropagators

More information

285K Homework #1. Sangchul Lee. April 28, 2017

285K Homework #1. Sangchul Lee. April 28, 2017 285K Homework #1 Sangchul Lee April 28, 2017 Problem 1. Suppose that X is a Banach space with respect to two norms: 1 and 2. Prove that if there is c (0, such that x 1 c x 2 for each x X, then there is

More information

TOEPLITZ OPERATORS. Toeplitz studied infinite matrices with NW-SE diagonals constant. f e C :

TOEPLITZ OPERATORS. Toeplitz studied infinite matrices with NW-SE diagonals constant. f e C : TOEPLITZ OPERATORS EFTON PARK 1. Introduction to Toeplitz Operators Otto Toeplitz lived from 1881-1940 in Goettingen, and it was pretty rough there, so he eventually went to Palestine and eventually contracted

More information

Applied Analysis (APPM 5440): Final exam 1:30pm 4:00pm, Dec. 14, Closed books.

Applied Analysis (APPM 5440): Final exam 1:30pm 4:00pm, Dec. 14, Closed books. Applied Analysis APPM 44: Final exam 1:3pm 4:pm, Dec. 14, 29. Closed books. Problem 1: 2p Set I = [, 1]. Prove that there is a continuous function u on I such that 1 ux 1 x sin ut 2 dt = cosx, x I. Define

More information

Functional Analysis Exercise Class

Functional Analysis Exercise Class Functional Analysis Exercise Class Week 2 November 6 November Deadline to hand in the homeworks: your exercise class on week 9 November 13 November Exercises (1) Let X be the following space of piecewise

More information

Generalized string equations for Hurwitz numbers

Generalized string equations for Hurwitz numbers Generalized string equations for Hurwitz numbers Kanehisa Takasaki December 17, 2010 1. Hurwitz numbers of Riemann sphere 2. Generating functions of Hurwitz numbers 3. Fermionic representation of tau functions

More information

Substitutions and symbolic dynamical systems, Lecture 7: The Rudin-Shapiro, Fibonacci and Chacon sequences

Substitutions and symbolic dynamical systems, Lecture 7: The Rudin-Shapiro, Fibonacci and Chacon sequences Substitutions and symbolic dynamical systems, Lecture 7: The Rudin-Shapiro, Fibonacci and Chacon sequences 18 mars 2015 The Rudin-Shapiro sequence The Rudin-Shapiro sequence is the fixed point u beginning

More information

Introduction to orthogonal polynomials. Michael Anshelevich

Introduction to orthogonal polynomials. Michael Anshelevich Introduction to orthogonal polynomials Michael Anshelevich November 6, 2003 µ = probability measure on R with finite moments m n (µ) = R xn dµ(x)

More information

Counting patterns in digital expansions

Counting patterns in digital expansions Counting patterns in digital expansions Manfred G. Madritsch Department for Analysis and Computational Number Theory Graz University of Technology madritsch@math.tugraz.at Séminaire ERNEST Dynamique, Arithmétique,

More information

Inverse function theorems

Inverse function theorems Inverse function theorems Ivar Ekeland CEREMADE, Université Paris-Dauphine Ravello, May 2011 Ivar Ekeland (CEREMADE, Université Paris-Dauphine) Inverse function theorems Ravello, May 2011 1 / 18 Ekeland

More information

10.1. The spectrum of an operator. Lemma If A < 1 then I A is invertible with bounded inverse

10.1. The spectrum of an operator. Lemma If A < 1 then I A is invertible with bounded inverse 10. Spectral theory For operators on finite dimensional vectors spaces, we can often find a basis of eigenvectors (which we use to diagonalize the matrix). If the operator is symmetric, this is always

More information

Calculus in Gauss Space

Calculus in Gauss Space Calculus in Gauss Space 1. The Gradient Operator The -dimensional Lebesgue space is the measurable space (E (E )) where E =[0 1) or E = R endowed with the Lebesgue measure, and the calculus of functions

More information

Manual for SOA Exam MLC.

Manual for SOA Exam MLC. Chapter 10. Poisson processes. Section 10.5. Nonhomogenous Poisson processes Extract from: Arcones Fall 2009 Edition, available at http://www.actexmadriver.com/ 1/14 Nonhomogenous Poisson processes Definition

More information

FURTHER NUMBER THEORY

FURTHER NUMBER THEORY FURTHER NUMBER THEORY SANJU VELANI. Dirichlet s Theorem and Continued Fractions Recall a fundamental theorem in the theory of Diophantine approximation. Theorem. (Dirichlet 842). For any real number x

More information

Quantum Searching. Robert-Jan Slager and Thomas Beuman. 24 november 2009

Quantum Searching. Robert-Jan Slager and Thomas Beuman. 24 november 2009 Quantum Searching Robert-Jan Slager and Thomas Beuman 24 november 2009 1 Introduction Quantum computers promise a significant speed-up over classical computers, since calculations can be done simultaneously.

More information

Lecture 4: Numerical solution of ordinary differential equations

Lecture 4: Numerical solution of ordinary differential equations Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor

More information

Continued Fractions New and Old Results

Continued Fractions New and Old Results Continued Fractions New and Old Results Jeffrey Shallit School of Computer Science University of Waterloo Waterloo, Ontario N2L 3G1 Canada shallit@cs.uwaterloo.ca https://www.cs.uwaterloo.ca/~shallit Joint

More information

RICE UNIVERSITY. Positive Lyapunov Exponent for Ergodic Schrödinger Operators. Helge Krüger

RICE UNIVERSITY. Positive Lyapunov Exponent for Ergodic Schrödinger Operators. Helge Krüger RICE UNIVERSITY Positive Lyapunov Exponent for Ergodic Schrödinger Operators by Helge Krüger A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE Doctor of Philosophy Approved,

More information

Introduction to Ergodic Theory

Introduction to Ergodic Theory Introduction to Ergodic Theory Marius Lemm May 20, 2010 Contents 1 Ergodic stochastic processes 2 1.1 Canonical probability space.......................... 2 1.2 Shift operators.................................

More information

The reference [Ho17] refers to the course lecture notes by Ilkka Holopainen.

The reference [Ho17] refers to the course lecture notes by Ilkka Holopainen. Department of Mathematics and Statistics Real Analysis I, Fall 207 Solutions to Exercise 6 (6 pages) riikka.schroderus at helsinki.fi Note. The course can be passed by an exam. The first possible exam

More information