Spectral analysis of two doubly infinite Jacobi operators

Size: px
Start display at page:

Download "Spectral analysis of two doubly infinite Jacobi operators"

Transcription

1 Spectral analysis of two doubly infinite Jacobi operators František Štampach jointly with Mourad E. H. Ismail Stockholm University Spectral Theory and Applications conference in memory of Boris Pavlov March 13, 2016 František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

2 Contents 1 Introduction 2 Spectral resolution of A 3 Spectral resolution of B František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

3 Two doubly-infinite Jacobi matrices We analyze the spectral properties of Jacobi operators A and B acting on vectors of the standard basis of l 2 Z) as: and Ae n = q n+1 e n 1 + q n e n+1, n Z, where q 0, 1) and α R. Be n = e n 1 + αq n e n + e n+1, n Z, František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

4 Two doubly-infinite Jacobi matrices We analyze the spectral properties of Jacobi operators A and B acting on vectors of the standard basis of l 2 Z) as: and Ae n = q n+1 e n 1 + q n e n+1, n Z, where q 0, 1) and α R. Be n = e n 1 + αq n e n + e n+1, n Z, The semi-infinite parts of A and B has been studied before mainly within the theory of Orthogonal Polynomials: František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

5 Two doubly-infinite Jacobi matrices We analyze the spectral properties of Jacobi operators A and B acting on vectors of the standard basis of l 2 Z) as: and Ae n = q n+1 e n 1 + q n e n+1, n Z, where q 0, 1) and α R. Be n = e n 1 + αq n e n + e n+1, n Z, The semi-infinite parts of A and B has been studied before mainly within the theory of Orthogonal Polynomials: A with n < 0 [Al-Salam, Ismail 83]: Rogers Ramanujan s functions and continued fractions. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

6 Two doubly-infinite Jacobi matrices We analyze the spectral properties of Jacobi operators A and B acting on vectors of the standard basis of l 2 Z) as: and Ae n = q n+1 e n 1 + q n e n+1, n Z, where q 0, 1) and α R. Be n = e n 1 + αq n e n + e n+1, n Z, The semi-infinite parts of A and B has been studied before mainly within the theory of Orthogonal Polynomials: A with n < 0 [Al-Salam, Ismail 83]: Rogers Ramanujan s functions and continued fractions. A with n > 0 [Chen, Ismail 98, FS 16]: Indeterminate moment problem, Nevanlinna functions. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

7 Two doubly-infinite Jacobi matrices We analyze the spectral properties of Jacobi operators A and B acting on vectors of the standard basis of l 2 Z) as: and Ae n = q n+1 e n 1 + q n e n+1, n Z, where q 0, 1) and α R. Be n = e n 1 + αq n e n + e n+1, n Z, The semi-infinite parts of A and B has been studied before mainly within the theory of Orthogonal Polynomials: A with n < 0 [Al-Salam, Ismail 83]: Rogers Ramanujan s functions and continued fractions. A with n > 0 [Chen, Ismail 98, FS 16]: Indeterminate moment problem, Nevanlinna functions. B with n 0 [Ismail, Mulla 87]: q-chebyshev polynomials. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

8 Two doubly-infinite Jacobi matrices We analyze the spectral properties of Jacobi operators A and B acting on vectors of the standard basis of l 2 Z) as: and Ae n = q n+1 e n 1 + q n e n+1, n Z, where q 0, 1) and α R. Be n = e n 1 + αq n e n + e n+1, n Z, The semi-infinite parts of A and B has been studied before mainly within the theory of Orthogonal Polynomials: A with n < 0 [Al-Salam, Ismail 83]: Rogers Ramanujan s functions and continued fractions. A with n > 0 [Chen, Ismail 98, FS 16]: Indeterminate moment problem, Nevanlinna functions. B with n 0 [Ismail, Mulla 87]: q-chebyshev polynomials. The spectrum of any associated semi-infinite Jacobi operator is never known explicitly α 0) but is expressible in terms of zeros of certain special functions.functions František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

9 Special functions 1/2 - basic hypergeometric series Let 0 < q < 1, r, s Z +. Recall the basic hypergeometric function [ ] a1, a r φ 2,... a r s ; q, z b 1, b 2,... b s is defined by the power series a 1 ; q) na 2 ; q) n... a r ; q) n 1) s r+1)n q s r+1)nn 1)/2 z n b n=0 1 ; q) nb 2 ; q) n... b s; q) n q; q) n where z, a 1, a 2,..., a r C, b 1, b 2,..., b s C \ q Z and a; q) 0 = 1, a; q) n = 1 a)1 aq)... 1 aq n 1 ) is the q-pochhammer symbol. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

10 Special functions 1/2 - basic hypergeometric series Let 0 < q < 1, r, s Z +. Recall the basic hypergeometric function [ ] a1, a r φ 2,... a r s ; q, z b 1, b 2,... b s is defined by the power series a 1 ; q) na 2 ; q) n... a r ; q) n 1) s r+1)n q s r+1)nn 1)/2 z n b n=0 1 ; q) nb 2 ; q) n... b s; q) n q; q) n where z, a 1, a 2,..., a r C, b 1, b 2,..., b s C \ q Z and a; q) 0 = 1, a; q) n = 1 a)1 aq)... 1 aq n 1 ) is the q-pochhammer symbol. Here we will need only 0 φ 1 and 1 φ 1. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

11 Special functions 2/2 - theta functions The theta function: 1 θ qz) := z; q) q/z; q) = q nn 1)/2 z) n q; q) n= František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

12 Special functions 2/2 - theta functions The theta function: 1 θ qz) := z; q) q/z; q) = q nn 1)/2 z) n q; q) n= Jacobi s theta functions: ϑ 1 z q) = iq 1/4 e iz q 2 ; q 2 ) θ q 2 e 2iz) ϑ 2 z q) = q 1/4 e iz q 2 ; q 2 ) θ q 2 e 2iz) ϑ 3 z q) = q 2 ; q 2 ) θ q 2 qe 2iz) ϑ 4 z q) = q 2 ; q 2 ) θ q 2 qe 2iz) František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

13 Contents 1 Introduction 2 Spectral resolution of A 3 Spectral resolution of B František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

14 Operators associated with the Jacobi matrix A Operator A with Dom A = spane n n Z} acting as has deficiency indices 1, 1). Ae n = q n+1 e n 1 + q n e n+1 František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

15 Operators associated with the Jacobi matrix A Operator A with Dom A = spane n n Z} acting as Ae n = q n+1 e n 1 + q n e n+1 has deficiency indices 1, 1). Thus, there is a one-parameter family of self-adjoint extensions A t, t R }. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

16 Operators associated with the Jacobi matrix A Operator A with Dom A = spane n n Z} acting as Ae n = q n+1 e n 1 + q n e n+1 has deficiency indices 1, 1). Thus, there is a one-parameter family of self-adjoint extensions A t, t R }. Let D := ψ l 2 Aψ l 2 }. By using the theory of self-adjoint extensions and simple structure of matrix A one gets: František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

17 Operators associated with the Jacobi matrix A Operator A with Dom A = spane n n Z} acting as Ae n = q n+1 e n 1 + q n e n+1 has deficiency indices 1, 1). Thus, there is a one-parameter family of self-adjoint extensions A t, t R }. Let D := ψ l 2 Aψ l 2 }. By using the theory of self-adjoint extensions and simple structure of matrix A one gets: Proposition self-adjoint extensions) For t R }, operators A t, acting as A t ψ = Aψ, with domains Dom A t = ψ D if t R, or are all self-adjoint extensions of A. lim n q n ψ 2n+1 + tψ 2n ) = 0 Dom A = ψ D } lim n q n qψ 2n 1 tψ 2n ) = 0, } lim n q n ψ 2n = 0, František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

18 Operators associated with the Jacobi matrix A Operator A with Dom A = spane n n Z} acting as Ae n = q n+1 e n 1 + q n e n+1 has deficiency indices 1, 1). Thus, there is a one-parameter family of self-adjoint extensions A t, t R }. Let D := ψ l 2 Aψ l 2 }. By using the theory of self-adjoint extensions and simple structure of matrix A one gets: Proposition self-adjoint extensions) For t R }, operators A t, acting as A t ψ = Aψ, with domains Dom A t = ψ D if t R, or are all self-adjoint extensions of A. In addition, lim n q n ψ 2n+1 + tψ 2n ) = 0 Dom A = ψ D } lim n q n qψ 2n 1 tψ 2n ) = 0, } lim n q n ψ 2n = 0, σ ca t ) = σessa t ) = 0}, t R }. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

19 Solutions of the eigenvalue equation The q-exponential function: E qz) = n=0 q n2 /4 ) z n = 1 φ 1 0; q 1/2 ; q 1/2, q 1/4 z q; q) n František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

20 Solutions of the eigenvalue equation The q-exponential function: Sequences ψ ±, where E qz) = n=0 q n2 /4 ) z n = 1 φ 1 0; q 1/2 ; q 1/2, q 1/4 z q; q) n ψ n ± := ±i) n q n/2 E q 2 ±ixq n ), are two linearly independent solutions of the difference equation Aψ = xψ. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

21 Solutions of the eigenvalue equation The q-exponential function: Sequences ψ ±, where E qz) = n=0 q n2 /4 ) z n = 1 φ 1 0; q 1/2 ; q 1/2, q 1/4 z q; q) n ψ n ± := ±i) n q n/2 E q 2 ±ixq n ), are two linearly independent solutions of the difference equation Aψ = xψ. By inspection of the asymptotic behavior of ψ ± n, as n ±, one gets: ψ ± l 2 + ), however, ψ ± / l 2 ). František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

22 Solutions of the eigenvalue equation The q-exponential function: Sequences ψ ±, where E qz) = n=0 q n2 /4 ) z n = 1 φ 1 0; q 1/2 ; q 1/2, q 1/4 z q; q) n ψ n ± := ±i) n q n/2 E q 2 ±ixq n ), are two linearly independent solutions of the difference equation Aψ = xψ. By inspection of the asymptotic behavior of ψ ± n, as n ±, one gets: ψ ± l 2 + ), however, ψ ± / l 2 ). Hence, one expects there are non-trivial coefficients a = ax) and b = bx) such that aψ + + bψ l 2 ). František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

23 The l 2 -solution Proposition For all x C \ 0}, the sequence ) ) ϕx) := θ q iq 1/2 x ψ ) x) + θ q iq 1/2 x ψ +) x), is the non-trivial solution of Aφ = xφ which belongs to l 2 Z). In addition, within the space l 2 Z), this solution is given uniquely up to a multiplicative constant. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

24 The l 2 -solution Proposition For all x C \ 0}, the sequence ) ) ϕx) := θ q iq 1/2 x ψ ) x) + θ q iq 1/2 x ψ +) x), is the non-trivial solution of Aφ = xφ which belongs to l 2 Z). In addition, within the space l 2 Z), this solution is given uniquely up to a multiplicative constant. Moreover, ϕ nx) = 1; q) x n q nn 1)/2 0φ 1 ; 0; q 2, q 2n+4 x 2) and František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

25 The l 2 -solution Proposition For all x C \ 0}, the sequence ) ) ϕx) := θ q iq 1/2 x ψ ) x) + θ q iq 1/2 x ψ +) x), is the non-trivial solution of Aφ = xφ which belongs to l 2 Z). In addition, within the space l 2 Z), this solution is given uniquely up to a multiplicative constant. Moreover, and ϕ nx) = 1; q) x n q nn 1)/2 0φ 1 ; 0; q 2, q 2n+4 x 2) q 2 ϕx) 2 ; q 2) 2 l 2 = 4 ) θ q; q 2 2 q 2 z 2) František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

26 The secular equation Theorem secular equation) For t R }, one has spec c A t ) = 0} and spec p A t ) coincides with the set of roots of the secular equation: xθ q 4 q 2 x 2) + tθ q 4 x 2) = 0, for t R, and θ q 4 x 2) = 0, for t =. In addition, all eigenvalues of A t are simple. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

27 The secular equation Theorem secular equation) For t R }, one has spec c A t ) = 0} and spec p A t ) coincides with the set of roots of the secular equation: xθ q 4 q 2 x 2) + tθ q 4 x 2) = 0, for t R, and θ q 4 x 2) = 0, for t =. In addition, all eigenvalues of A t are simple. Corollary: } spec p A 0 ) = ±q 2n+1 n Z } and spec p A ) = ±q 2n n Z. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

28 The secular equation Theorem secular equation) For t R }, one has spec c A t ) = 0} and spec p A t ) coincides with the set of roots of the secular equation: xθ q 4 q 2 x 2) + tθ q 4 x 2) = 0, for t R, and θ q 4 x 2) = 0, for t =. In addition, all eigenvalues of A t are simple. Corollary: } spec p A 0 ) = ±q 2n+1 n Z } and spec p A ) = ±q 2n n Z. How to solve the secular equation in general? František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

29 The secular equation Theorem secular equation) For t R }, one has spec c A t ) = 0} and spec p A t ) coincides with the set of roots of the secular equation: xθ q 4 q 2 x 2) + tθ q 4 x 2) = 0, for t R, and θ q 4 x 2) = 0, for t =. In addition, all eigenvalues of A t are simple. Corollary: } spec p A 0 ) = ±q 2n+1 n Z } and spec p A ) = ±q 2n n Z. How to solve the secular equation in general? Reparametrize t = Φs) and use nice properties of Jacobi s theta functions... František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

30 Reparametrization t = Φs) The function Φs) := iq 1/2 ϑ 4 is q 2 ) ϑ ) 1 is q 2 is real-valued, strictly decreasing on 0, 2 ln q), and maps [0, 2 ln q) onto R }. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

31 Reparametrization t = Φs) The function Φs) := iq 1/2 ϑ 4 is q 2 ) ϑ ) 1 is q 2 is real-valued, strictly decreasing on 0, 2 ln q), and maps [0, 2 ln q) onto R }. For the inverse function, one has Φ 1 dx t) = Cq) Dq) ) t + x 2 q + x 2) where and Cq) = q 1/2 ϑ 2 0 q 2 ) ϑ 3 0 q 2 ) 0 q 2 ) Dq) = qϑ2 3 ϑ 2 ). 2 0 q 2 František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

32 Spectrum fully explicitly Using the reparametrization t = Φs), the secular equation simplifies to θ q 2 e s x ) θ q 2 e s x ) = 0. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

33 Spectrum fully explicitly Using the reparametrization t = Φs), the secular equation simplifies to θ q 2 e s x ) θ q 2 e s x ) = 0. Theorem Let t R }, then where spec p A t ) = e s q 2Z e s q 2Z dx s = Cq) Dq) ) t + x 2 q + x 2). In addition, the family of corresponding eigenvectors ϕ±e s q 2N ) N Z}, where forms an orthogonal basis of l 2 Z). ϕ nx) = 1; q) x n q nn 1)/2 0φ 1 ; 0; q 2, q 2n+4 x 2), František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

34 Contents 1 Introduction 2 Spectral resolution of A 3 Spectral resolution of B František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

35 The discrete Schrödinger operator B The second Jacobi matrix B determines the unique operator B = U + U + αv where U is the forward shift operator and V is the self-adjoint diagonal operator: Ue n = e n+1 and Ve n = q n e n, n Z. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

36 The discrete Schrödinger operator B The second Jacobi matrix B determines the unique operator B = U + U + αv where U is the forward shift operator and V is the self-adjoint diagonal operator: Ue n = e n+1 and Ve n = q n e n, n Z. Proposition essential spectrum) The operator B is self-adjoint and one has σessb) = [ 2, 2]. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

37 Solutions of the eigenvalue equation Using properties of the Hahn-Exton q-bessel functions one verifies: f nz) := 1) n α n q 2 1 nn+1) z 1 α 1 q n+1 ; q ) 1 φ 1 0; z 1 α 1 q n+1 ; q, zα 1 q n+1) and g nz) := z n zαq 1 n ; q ) 1 φ 1 0; zαq 1 n ; q, qz 2). are two solutions of the equation for all α, z 0. Bψ = z + z 1 )ψ. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

38 Solutions of the eigenvalue equation Using properties of the Hahn-Exton q-bessel functions one verifies: f nz) := 1) n α n q 2 1 nn+1) z 1 α 1 q n+1 ; q ) 1 φ 1 0; z 1 α 1 q n+1 ; q, zα 1 q n+1) and g nz) := z n zαq 1 n ; q ) 1 φ 1 0; zαq 1 n ; q, qz 2). are two solutions of the equation for all α, z 0. Note that z z + z 1 : Bψ = z + z 1 )ψ. z 0 < z < 1} C \ [ 2, 2], outside σessa) ) e iθ θ [0, π]} [ 2, 2], inside σessa) ) František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

39 Solutions of the eigenvalue equation Using properties of the Hahn-Exton q-bessel functions one verifies: f nz) := 1) n α n q 2 1 nn+1) z 1 α 1 q n+1 ; q ) 1 φ 1 0; z 1 α 1 q n+1 ; q, zα 1 q n+1) and g nz) := z n zαq 1 n ; q ) 1 φ 1 0; zαq 1 n ; q, qz 2). are two solutions of the equation for all α, z 0. Note that z z + z 1 : Bψ = z + z 1 )ψ. z 0 < z < 1} C \ [ 2, 2], outside σessa) ) e iθ θ [0, π]} [ 2, 2], inside σessa) ) The solutions f z) and gz) are linearly independent iff z / α 1 q Z 0} since W f, g) = z 1 θ q αz). František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

40 The point spectrum Detailed asymptotic analysis of solutions f and g yields: If 0 < z < 1 and z / α 1 q Z 0}, then l 2 + ) / l 2 + ) f z) / l 2 gz) ) l 2 ). František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

41 The point spectrum Detailed asymptotic analysis of solutions f and g yields: If 0 < z < 1 and z / α 1 q Z 0}, then l 2 + ) / l 2 + ) f z) / l 2 gz) ) l 2 ) If z = 1 the asymptotic behavior of solutions is very different and, in the end, it implies that: For α R and x [ 2, 2], there is no non-trivial solution of Bψ = xψ belonging to l 2 Z).. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

42 The point spectrum Detailed asymptotic analysis of solutions f and g yields: If 0 < z < 1 and z / α 1 q Z 0}, then l 2 + ) / l 2 + ) f z) / l 2 gz) ) l 2 ) If z = 1 the asymptotic behavior of solutions is very different and, in the end, it implies that: For α R and x [ 2, 2], there is no non-trivial solution of Bψ = xψ belonging to l 2 Z). Theorem point spectrum) If α 0, then } σb) \ [ 2, 2] = σ pb) = α 1 q m + αq m m > log q α and all points from this set are simple eigenvalues of B.. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

43 The point spectrum Detailed asymptotic analysis of solutions f and g yields: If 0 < z < 1 and z / α 1 q Z 0}, then l 2 + ) / l 2 + ) f z) / l 2 gz) ) l 2 ) If z = 1 the asymptotic behavior of solutions is very different and, in the end, it implies that: For α R and x [ 2, 2], there is no non-trivial solution of Bψ = xψ belonging to l 2 Z). Theorem point spectrum) If α 0, then } σb) \ [ 2, 2] = σ pb) = α 1 q m + αq m m > log q α and all points from this set are simple eigenvalues of B. Further, eigenvectors v m corresponding to eigenvalues α 1 q m + αq m, can be chosen as v m = v m,j } j=, with v m,j = f j α 1 q m) = 1) j α j q 1 2 jj+1) q m+j+1 ; q) 1 φ 1 0; q m+j+1 ; q, α 2 q m+j+1).. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

44 The point spectrum Detailed asymptotic analysis of solutions f and g yields: If 0 < z < 1 and z / α 1 q Z 0}, then l 2 + ) / l 2 + ) f z) / l 2 gz) ) l 2 ) If z = 1 the asymptotic behavior of solutions is very different and, in the end, it implies that: For α R and x [ 2, 2], there is no non-trivial solution of Bψ = xψ belonging to l 2 Z). Theorem point spectrum) If α 0, then } σb) \ [ 2, 2] = σ pb) = α 1 q m + αq m m > log q α and all points from this set are simple eigenvalues of B. Further, eigenvectors v m corresponding to eigenvalues α 1 q m + αq m, can be chosen as v m = v m,j } j=, with v m,j = f j α 1 q m) = 1) j α j q 1 2 jj+1) q m+j+1 ; q) 1 φ 1 0; q m+j+1 ; q, α 2 q m+j+1).. In addition, v m l 2 Z) = α m q mm+1)/2 1 α 2 q 2m q; q), m > log q α. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

45 The absolutely continuous part of the spectral measure Let us denote E k,l ) := e k, E B )e l, k, l Z, where E B stands for the spectral measure of B. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

46 The absolutely continuous part of the spectral measure Let us denote E k,l ) := e k, E B )e l, k, l Z, where E B stands for the spectral measure of B. To determine the spectral measure in the essential spectrum we use the formula E k,l a, b)) = lim lim 1 b δ Gk,l x + iɛ) G k,l x iɛ) ) dx, δ 0+ ɛ 0+ 2πi a+δ where G k,l z) := e k, B z) 1 e l = 1 g k z)f l z), k l, W f, g) g l z)f k z), k l, František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

47 The absolutely continuous part of the spectral measure Let us denote E k,l ) := e k, E B )e l, k, l Z, where E B stands for the spectral measure of B. To determine the spectral measure in the essential spectrum we use the formula E k,l a, b)) = lim lim 1 b δ Gk,l x + iɛ) G k,l x iɛ) ) dx, δ 0+ ɛ 0+ 2πi a+δ where G k,l z) := e k, B z) 1 e l = 1 g k z)f l z), k l, W f, g) g l z)f k z), k l, Proposition Let α 0 and 2 a < b 2. Then for any k, l Z, it holds E k,l [a, b]) = 1 2π φa f l φ b e iφ) f k e iφ) e 2iφ ; q ) 2 αe iφ, qα 1 e iφ ; q ) dφ where φ a = arccos a/2) and φ b = arccos b/2). Consequently, σ acb) = [ 2, 2]. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

48 Summary Theorem If α 0, then and σ essb) = σ acb) = [ 2, 2], } σ pb) = σ d B) = α 1 q m + αq m m > log q α. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

49 Summary Theorem If α 0, then and σ essb) = σ acb) = [ 2, 2], } σ pb) = σ d B) = α 1 q m + αq m m > log q α. In addition, for M R a Borel set, we have E k,l M) = 1 f l e iφ) f k e iφ) e 2iφ ; q ) 2 2π 2 cos φ [ 2,2] A αe iφ, qα 1 e iφ ; q ) dφ α 2 q 2m) α 2m q mm+1) f l α 1 q m) f k α 1 q m). q; q) 2 m> log α α 1 q m +αq m M František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

50 One consequence for special functions Recall the Hanhn-Exton or third Jackson s) q-bessel function is defined as J νz; q) = z ν qν+1 ; q) q; q) 1φ 1 0; q ν+1 ; q, qz 2). František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

51 One consequence for special functions Recall the Hanhn-Exton or third Jackson s) q-bessel function is defined as J νz; q) = z ν qν+1 ; q) q; q) 1φ 1 0; q ν+1 ; q, qz 2). Elements of the eigenvectors v m are expressible in terms of Hanhn-Exton q-bessel function and the formula for the norm of the eigenvectors yields n Z Jn 2 z; q) = 1, z < 1. 1 z2 František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

52 One consequence for special functions Recall the Hanhn-Exton or third Jackson s) q-bessel function is defined as J νz; q) = z ν qν+1 ; q) q; q) 1φ 1 0; q ν+1 ; q, qz 2). Elements of the eigenvectors v m are expressible in terms of Hanhn-Exton q-bessel function and the formula for the norm of the eigenvectors yields n Z Jn 2 z; q) = 1, z < 1. 1 z2 This formula seems to be new Really?) and it generalizes the well-known summation formula for the Bessel functions of the first kind: Jn 2 z) = 1, z < 1. n Z František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

53 The end Thank you! František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20

SPECTRAL THEOREM FOR COMPACT SELF-ADJOINT OPERATORS

SPECTRAL THEOREM FOR COMPACT SELF-ADJOINT OPERATORS SPECTRAL THEOREM FOR COMPACT SELF-ADJOINT OPERATORS G. RAMESH Contents Introduction 1 1. Bounded Operators 1 1.3. Examples 3 2. Compact Operators 5 2.1. Properties 6 3. The Spectral Theorem 9 3.3. Self-adjoint

More information

j=1 u 1jv 1j. 1/ 2 Lemma 1. An orthogonal set of vectors must be linearly independent.

j=1 u 1jv 1j. 1/ 2 Lemma 1. An orthogonal set of vectors must be linearly independent. Lecture Notes: Orthogonal and Symmetric Matrices Yufei Tao Department of Computer Science and Engineering Chinese University of Hong Kong taoyf@cse.cuhk.edu.hk Orthogonal Matrix Definition. Let u = [u

More information

On rational approximation of algebraic functions. Julius Borcea. Rikard Bøgvad & Boris Shapiro

On rational approximation of algebraic functions. Julius Borcea. Rikard Bøgvad & Boris Shapiro On rational approximation of algebraic functions http://arxiv.org/abs/math.ca/0409353 Julius Borcea joint work with Rikard Bøgvad & Boris Shapiro 1. Padé approximation: short overview 2. A scheme of rational

More information

The following definition is fundamental.

The following definition is fundamental. 1. Some Basics from Linear Algebra With these notes, I will try and clarify certain topics that I only quickly mention in class. First and foremost, I will assume that you are familiar with many basic

More information

The roots of The Third Jackson q-bessel Function

The roots of The Third Jackson q-bessel Function The roots of The Third Jacson q-bessel Function L. D. Abreu Departamento de Matemática Universidade de Coimbra Coimbra, Portugal J. Bustoz Department of Mathematics, Arizona State University Tempe, Arizona

More information

Kaczmarz algorithm in Hilbert space

Kaczmarz algorithm in Hilbert space STUDIA MATHEMATICA 169 (2) (2005) Kaczmarz algorithm in Hilbert space by Rainis Haller (Tartu) and Ryszard Szwarc (Wrocław) Abstract The aim of the Kaczmarz algorithm is to reconstruct an element in a

More information

MP463 QUANTUM MECHANICS

MP463 QUANTUM MECHANICS MP463 QUANTUM MECHANICS Introduction Quantum theory of angular momentum Quantum theory of a particle in a central potential - Hydrogen atom - Three-dimensional isotropic harmonic oscillator (a model of

More information

BIRTH-DEATH PROCESSES AND q-continued FRACTIONS

BIRTH-DEATH PROCESSES AND q-continued FRACTIONS BIRTH-DEATH PROCESSES AND q-continued FRACTIONS TONY FENG, RACHEL KIRSCH, ELISE MCCALL, AND MATT WAGE Abstract. In the 997 paper of Parthasarathy, Lenin, Schoutens, and Van Assche [6], the authors study

More information

Orthogonal Polynomial Ensembles

Orthogonal Polynomial Ensembles Chater 11 Orthogonal Polynomial Ensembles 11.1 Orthogonal Polynomials of Scalar rgument Let wx) be a weight function on a real interval, or the unit circle, or generally on some curve in the comlex lane.

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

Moreover this binary operation satisfies the following properties

Moreover this binary operation satisfies the following properties Contents 1 Algebraic structures 1 1.1 Group........................................... 1 1.1.1 Definitions and examples............................. 1 1.1.2 Subgroup.....................................

More information

Jacobi-Angelesco multiple orthogonal polynomials on an r-star

Jacobi-Angelesco multiple orthogonal polynomials on an r-star M. Leurs Jacobi-Angelesco m.o.p. 1/19 Jacobi-Angelesco multiple orthogonal polynomials on an r-star Marjolein Leurs, (joint work with Walter Van Assche) Conference on Orthogonal Polynomials and Holomorphic

More information

Spectral Analysis of Matrices - An Introduction for Engineers

Spectral Analysis of Matrices - An Introduction for Engineers Spectral Analysis of Matrices - An Introduction for Engineers Timo Weidl Department of Mathematics Royal Institute of Technology S-144 Stockholm weidl@mathkthse January 13, 26 Plan of the Talk 1 A practical

More information

MATH 23a, FALL 2002 THEORETICAL LINEAR ALGEBRA AND MULTIVARIABLE CALCULUS Solutions to Final Exam (in-class portion) January 22, 2003

MATH 23a, FALL 2002 THEORETICAL LINEAR ALGEBRA AND MULTIVARIABLE CALCULUS Solutions to Final Exam (in-class portion) January 22, 2003 MATH 23a, FALL 2002 THEORETICAL LINEAR ALGEBRA AND MULTIVARIABLE CALCULUS Solutions to Final Exam (in-class portion) January 22, 2003 1. True or False (28 points, 2 each) T or F If V is a vector space

More information

Operator algebras related to bounded positive operator. operator with simple spectrum

Operator algebras related to bounded positive operator. operator with simple spectrum Operator algebras related to bounded positive operator with simple spectrum Institute of Mathematics, University of Bia lystok joint work with A. Odzijewicz, A. Tereszkiewicz 13 June 2007 Varna Coherent

More information

Mathematical foundations - linear algebra

Mathematical foundations - linear algebra Mathematical foundations - linear algebra Andrea Passerini passerini@disi.unitn.it Machine Learning Vector space Definition (over reals) A set X is called a vector space over IR if addition and scalar

More information

Functional Analysis Review

Functional Analysis Review Outline 9.520: Statistical Learning Theory and Applications February 8, 2010 Outline 1 2 3 4 Vector Space Outline A vector space is a set V with binary operations +: V V V and : R V V such that for all

More information

Pseudospectra and Nonnormal Dynamical Systems

Pseudospectra and Nonnormal Dynamical Systems Pseudospectra and Nonnormal Dynamical Systems Mark Embree and Russell Carden Computational and Applied Mathematics Rice University Houston, Texas ELGERSBURG MARCH 1 Overview of the Course These lectures

More information

SPECTRAL THEOREM FOR SYMMETRIC OPERATORS WITH COMPACT RESOLVENT

SPECTRAL THEOREM FOR SYMMETRIC OPERATORS WITH COMPACT RESOLVENT SPECTRAL THEOREM FOR SYMMETRIC OPERATORS WITH COMPACT RESOLVENT Abstract. These are the letcure notes prepared for the workshop on Functional Analysis and Operator Algebras to be held at NIT-Karnataka,

More information

EIGENVALUES AND EIGENVECTORS 3

EIGENVALUES AND EIGENVECTORS 3 EIGENVALUES AND EIGENVECTORS 3 1. Motivation 1.1. Diagonal matrices. Perhaps the simplest type of linear transformations are those whose matrix is diagonal (in some basis). Consider for example the matrices

More information

Matrix orthogonal polynomials and group representations

Matrix orthogonal polynomials and group representations Matrix orthogonal polynomials and group representations Inés Pacharoni Universidad Nacional de Córdoba, Argentina (Joint work with Juan Tirao) Discrete Systems and Special Functions Cambridge, UK. July

More information

JUST THE MATHS UNIT NUMBER 9.9. MATRICES 9 (Modal & spectral matrices) A.J.Hobson

JUST THE MATHS UNIT NUMBER 9.9. MATRICES 9 (Modal & spectral matrices) A.J.Hobson JUST THE MATHS UNIT NUMBER 9.9 MATRICES 9 (Modal & spectral matrices) by A.J.Hobson 9.9. Assumptions and definitions 9.9.2 Diagonalisation of a matrix 9.9.3 Exercises 9.9.4 Answers to exercises UNIT 9.9

More information

Matrix Lie groups. and their Lie algebras. Mahmood Alaghmandan. A project in fulfillment of the requirement for the Lie algebra course

Matrix Lie groups. and their Lie algebras. Mahmood Alaghmandan. A project in fulfillment of the requirement for the Lie algebra course Matrix Lie groups and their Lie algebras Mahmood Alaghmandan A project in fulfillment of the requirement for the Lie algebra course Department of Mathematics and Statistics University of Saskatchewan March

More information

13 Maximum Modulus Principle

13 Maximum Modulus Principle 3 Maximum Modulus Principle Theorem 3. (maximum modulus principle). If f is non-constant and analytic on an open connected set Ω, then there is no point z 0 Ω such that f(z) f(z 0 ) for all z Ω. Remark

More information

1. Matrix multiplication and Pauli Matrices: Pauli matrices are the 2 2 matrices. 1 0 i 0. 0 i

1. Matrix multiplication and Pauli Matrices: Pauli matrices are the 2 2 matrices. 1 0 i 0. 0 i Problems in basic linear algebra Science Academies Lecture Workshop at PSGRK College Coimbatore, June 22-24, 2016 Govind S. Krishnaswami, Chennai Mathematical Institute http://www.cmi.ac.in/~govind/teaching,

More information

MATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators.

MATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators. MATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators. Adjoint operator and adjoint matrix Given a linear operator L on an inner product space V, the adjoint of L is a transformation

More information

Asymptotic distribution of eigenvalues of Laplace operator

Asymptotic distribution of eigenvalues of Laplace operator Asymptotic distribution of eigenvalues of Laplace operator 23.8.2013 Topics We will talk about: the number of eigenvalues of Laplace operator smaller than some λ as a function of λ asymptotic behaviour

More information

Math Ordinary Differential Equations

Math Ordinary Differential Equations Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x

More information

Review problems for MA 54, Fall 2004.

Review problems for MA 54, Fall 2004. Review problems for MA 54, Fall 2004. Below are the review problems for the final. They are mostly homework problems, or very similar. If you are comfortable doing these problems, you should be fine on

More information

MOMENTS OF HYPERGEOMETRIC HURWITZ ZETA FUNCTIONS

MOMENTS OF HYPERGEOMETRIC HURWITZ ZETA FUNCTIONS MOMENTS OF HYPERGEOMETRIC HURWITZ ZETA FUNCTIONS ABDUL HASSEN AND HIEU D. NGUYEN Abstract. This paper investigates a generalization the classical Hurwitz zeta function. It is shown that many of the properties

More information

MATH 304 Linear Algebra Lecture 20: The Gram-Schmidt process (continued). Eigenvalues and eigenvectors.

MATH 304 Linear Algebra Lecture 20: The Gram-Schmidt process (continued). Eigenvalues and eigenvectors. MATH 304 Linear Algebra Lecture 20: The Gram-Schmidt process (continued). Eigenvalues and eigenvectors. Orthogonal sets Let V be a vector space with an inner product. Definition. Nonzero vectors v 1,v

More information

Eigenvalues and Eigenvectors

Eigenvalues and Eigenvectors Sec. 6.1 Eigenvalues and Eigenvectors Linear transformations L : V V that go from a vector space to itself are often called linear operators. Many linear operators can be understood geometrically by identifying

More information

Symmetric and self-adjoint matrices

Symmetric and self-adjoint matrices Symmetric and self-adjoint matrices A matrix A in M n (F) is called symmetric if A T = A, ie A ij = A ji for each i, j; and self-adjoint if A = A, ie A ij = A ji or each i, j Note for A in M n (R) that

More information

Analysis IV : Assignment 3 Solutions John Toth, Winter ,...). In particular for every fixed m N the sequence (u (n)

Analysis IV : Assignment 3 Solutions John Toth, Winter ,...). In particular for every fixed m N the sequence (u (n) Analysis IV : Assignment 3 Solutions John Toth, Winter 203 Exercise (l 2 (Z), 2 ) is a complete and separable Hilbert space. Proof Let {u (n) } n N be a Cauchy sequence. Say u (n) = (..., u 2, (n) u (n),

More information

DETERMINANT IDENTITIES FOR THETA FUNCTIONS

DETERMINANT IDENTITIES FOR THETA FUNCTIONS DETERMINANT IDENTITIES FOR THETA FUNCTIONS SHAUN COOPER AND PEE CHOON TOH Abstract. Two proofs of a theta function identity of R. W. Gosper and R. Schroeppel are given. A cubic analogue is presented, and

More information

Foundations of Matrix Analysis

Foundations of Matrix Analysis 1 Foundations of Matrix Analysis In this chapter we recall the basic elements of linear algebra which will be employed in the remainder of the text For most of the proofs as well as for the details, the

More information

Quadratic forms. Here. Thus symmetric matrices are diagonalizable, and the diagonalization can be performed by means of an orthogonal matrix.

Quadratic forms. Here. Thus symmetric matrices are diagonalizable, and the diagonalization can be performed by means of an orthogonal matrix. Quadratic forms 1. Symmetric matrices An n n matrix (a ij ) n ij=1 with entries on R is called symmetric if A T, that is, if a ij = a ji for all 1 i, j n. We denote by S n (R) the set of all n n symmetric

More information

Krein s formula and perturbation theory

Krein s formula and perturbation theory ISSN: 40-567 Krein s formula and perturbation theory P.Kurasov S.T.Kuroda Research Reports in athematics Number 6, 2000 Department of athematics Stockholm University Electronic versions of this document

More information

A Right Inverse Of The Askey-Wilson Operator arxiv:math/ v1 [math.ca] 5 Oct 1993 B. Malcolm Brown and Mourad E. H. Ismail

A Right Inverse Of The Askey-Wilson Operator arxiv:math/ v1 [math.ca] 5 Oct 1993 B. Malcolm Brown and Mourad E. H. Ismail A Right Inverse Of The Askey-Wilson Operator arxiv:math/9310219v1 [math.ca] 5 Oct 1993 B. Malcolm Brown and Mourad E. H. Ismail Abstract We establish an integral representation of a right inverse of the

More information

ACM 104. Homework Set 4 Solutions February 14, 2001

ACM 104. Homework Set 4 Solutions February 14, 2001 ACM 04 Homework Set 4 Solutions February 4, 00 Franklin Chapter, Problem 4, page 55 Suppose that we feel that some observations are more important or reliable than others Redefine the function to be minimized

More information

Fredholm Theory. April 25, 2018

Fredholm Theory. April 25, 2018 Fredholm Theory April 25, 208 Roughly speaking, Fredholm theory consists of the study of operators of the form I + A where A is compact. From this point on, we will also refer to I + A as Fredholm operators.

More information

1 Distributions (due January 22, 2009)

1 Distributions (due January 22, 2009) Distributions (due January 22, 29). The distribution derivative of the locally integrable function ln( x ) is the principal value distribution /x. We know that, φ = lim φ(x) dx. x ɛ x Show that x, φ =

More information

Definition and basic properties of heat kernels I, An introduction

Definition and basic properties of heat kernels I, An introduction Definition and basic properties of heat kernels I, An introduction Zhiqin Lu, Department of Mathematics, UC Irvine, Irvine CA 92697 April 23, 2010 In this lecture, we will answer the following questions:

More information

Trace Class Operators and Lidskii s Theorem

Trace Class Operators and Lidskii s Theorem Trace Class Operators and Lidskii s Theorem Tom Phelan Semester 2 2009 1 Introduction The purpose of this paper is to provide the reader with a self-contained derivation of the celebrated Lidskii Trace

More information

Basic Calculus Review

Basic Calculus Review Basic Calculus Review Lorenzo Rosasco ISML Mod. 2 - Machine Learning Vector Spaces Functionals and Operators (Matrices) Vector Space A vector space is a set V with binary operations +: V V V and : R V

More information

Problem 1A. Find the volume of the solid given by x 2 + z 2 1, y 2 + z 2 1. (Hint: 1. Solution: The volume is 1. Problem 2A.

Problem 1A. Find the volume of the solid given by x 2 + z 2 1, y 2 + z 2 1. (Hint: 1. Solution: The volume is 1. Problem 2A. Problem 1A Find the volume of the solid given by x 2 + z 2 1, y 2 + z 2 1 (Hint: 1 1 (something)dz) Solution: The volume is 1 1 4xydz where x = y = 1 z 2 This integral has value 16/3 Problem 2A Let f(x)

More information

RIESZ BASES AND UNCONDITIONAL BASES

RIESZ BASES AND UNCONDITIONAL BASES In this paper we give a brief introduction to adjoint operators on Hilbert spaces and a characterization of the dual space of a Hilbert space. We then introduce the notion of a Riesz basis and give some

More information

Real Variables # 10 : Hilbert Spaces II

Real Variables # 10 : Hilbert Spaces II randon ehring Real Variables # 0 : Hilbert Spaces II Exercise 20 For any sequence {f n } in H with f n = for all n, there exists f H and a subsequence {f nk } such that for all g H, one has lim (f n k,

More information

Chapter 1. Preliminaries. The purpose of this chapter is to provide some basic background information. Linear Space. Hilbert Space.

Chapter 1. Preliminaries. The purpose of this chapter is to provide some basic background information. Linear Space. Hilbert Space. Chapter 1 Preliminaries The purpose of this chapter is to provide some basic background information. Linear Space Hilbert Space Basic Principles 1 2 Preliminaries Linear Space The notion of linear space

More information

1 Fourier transform as unitary equivalence

1 Fourier transform as unitary equivalence Tel Aviv University, 009 Intro to functional analysis 1 1 Fourier transform as unitary equivalence 1a Introduction..................... 1 1b Exponential map................... 1c Exponential map as an

More information

Dmitry Pelinovsky Department of Mathematics, McMaster University, Canada

Dmitry Pelinovsky Department of Mathematics, McMaster University, Canada Spectrum of the linearized NLS problem Dmitry Pelinovsky Department of Mathematics, McMaster University, Canada Collaboration: Marina Chugunova (McMaster, Canada) Scipio Cuccagna (Modena and Reggio Emilia,

More information

ARTICLE IN PRESS. JID:YAIMA AID:6155 /FLA [m1l; v1.231; Prn:21/02/2018; 9:28] P.1 (1-56) Advances in Mathematics ( )

ARTICLE IN PRESS. JID:YAIMA AID:6155 /FLA [m1l; v1.231; Prn:21/02/2018; 9:28] P.1 (1-56) Advances in Mathematics ( ) JID:YAIMA AID:655 /FLA [ml; v.3; Prn:/0/08; 9:8] P. -56 Advances in Mathematics Contents lists available at ScienceDirect Advances in Mathematics 5 www.elsevier.com/locate/aim 5 Nevanlinna theory of the

More information

An identity of Andrews and the Askey-Wilson integral

An identity of Andrews and the Askey-Wilson integral Ramanujan J DOI 0.007/s39-008-922-4 An identity of Andrews and the Askey-Wilson integral Zhi-Guo Liu Received: 6 July 2007 / Accepted: 7 January 2008 Springer Science+Business Media, LLC 2008 Abstract

More information

Heisenberg s inequality

Heisenberg s inequality Heisenberg s inequality Michael Walter 1 Introduction 1 Definition. Let P be a probability measure on R. Then its mean or expectation is defined as E(P) := x dp(x) Its variance is the expected square deviation

More information

Central Limit Theorems for linear statistics for Biorthogonal Ensembles

Central Limit Theorems for linear statistics for Biorthogonal Ensembles Central Limit Theorems for linear statistics for Biorthogonal Ensembles Maurice Duits, Stockholm University Based on joint work with Jonathan Breuer (HUJI) Princeton, April 2, 2014 M. Duits (SU) CLT s

More information

MORE NOTES FOR MATH 823, FALL 2007

MORE NOTES FOR MATH 823, FALL 2007 MORE NOTES FOR MATH 83, FALL 007 Prop 1.1 Prop 1. Lemma 1.3 1. The Siegel upper half space 1.1. The Siegel upper half space and its Bergman kernel. The Siegel upper half space is the domain { U n+1 z C

More information

Math Homework 8 (selected problems)

Math Homework 8 (selected problems) Math 102 - Homework 8 (selected problems) David Lipshutz Problem 1. (Strang, 5.5: #14) In the list below, which classes of matrices contain A and which contain B? 1 1 1 1 A 0 0 1 0 0 0 0 1 and B 1 1 1

More information

Eigenvalues and Eigenvectors

Eigenvalues and Eigenvectors Eigenvalues and Eigenvectors Definition 0 Let A R n n be an n n real matrix A number λ R is a real eigenvalue of A if there exists a nonzero vector v R n such that A v = λ v The vector v is called an eigenvector

More information

Formalism of quantum mechanics

Formalism of quantum mechanics Formalism of quantum mechanics August 22, 204 Contents Introduction 2 Classical physics 3 Quantum physics 3 4 Harmonic oscillator 0 Introduction Quite a bit of the serious mathematical theory of self-adjoint

More information

Integer Partitions With Even Parts Below Odd Parts and the Mock Theta Functions

Integer Partitions With Even Parts Below Odd Parts and the Mock Theta Functions Integer Partitions With Even Parts Below Odd Parts and the Mock Theta Functions by George E. Andrews Key Words: Partitions, mock theta functions, crank AMS Classification Numbers: P84, P83, P8, 33D5 Abstract

More information

5 Compact linear operators

5 Compact linear operators 5 Compact linear operators One of the most important results of Linear Algebra is that for every selfadjoint linear map A on a finite-dimensional space, there exists a basis consisting of eigenvectors.

More information

Born-Oppenheimer Approximation in a singular system Prague. O.T.Turgut Bogazici Univ.

Born-Oppenheimer Approximation in a singular system Prague. O.T.Turgut Bogazici Univ. Born-Oppenheimer Approximation in a singular system Prague O.T.Turgut Bogazici Univ. June 5, 206 How can we use Born-Oppenheimer in a singular system (if we can)? What is the usual approach: Let us consider

More information

Real Analysis Problems

Real Analysis Problems Real Analysis Problems Cristian E. Gutiérrez September 14, 29 1 1 CONTINUITY 1 Continuity Problem 1.1 Let r n be the sequence of rational numbers and Prove that f(x) = 1. f is continuous on the irrationals.

More information

Stable Process. 2. Multivariate Stable Distributions. July, 2006

Stable Process. 2. Multivariate Stable Distributions. July, 2006 Stable Process 2. Multivariate Stable Distributions July, 2006 1. Stable random vectors. 2. Characteristic functions. 3. Strictly stable and symmetric stable random vectors. 4. Sub-Gaussian random vectors.

More information

Further Mathematical Methods (Linear Algebra)

Further Mathematical Methods (Linear Algebra) Further Mathematical Methods (Linear Algebra) Solutions For The Examination Question (a) To be an inner product on the real vector space V a function x y which maps vectors x y V to R must be such that:

More information

Physics 505 Homework No. 1 Solutions S1-1

Physics 505 Homework No. 1 Solutions S1-1 Physics 505 Homework No s S- Some Preliminaries Assume A and B are Hermitian operators (a) Show that (AB) B A dx φ ABψ dx (A φ) Bψ dx (B (A φ)) ψ dx (B A φ) ψ End (b) Show that AB [A, B]/2+{A, B}/2 where

More information

Lecture I: Asymptotics for large GUE random matrices

Lecture I: Asymptotics for large GUE random matrices Lecture I: Asymptotics for large GUE random matrices Steen Thorbjørnsen, University of Aarhus andom Matrices Definition. Let (Ω, F, P) be a probability space and let n be a positive integer. Then a random

More information

Bessel s Equation. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Fall Department of Mathematics

Bessel s Equation. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Fall Department of Mathematics Bessel s Equation MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018 Background Bessel s equation of order ν has the form where ν is a constant. x 2 y + xy

More information

Some Combinatorial and Analytical Identities

Some Combinatorial and Analytical Identities Some Combinatorial and Analytical Identities Mourad EH Ismail Dennis Stanton November 6, 200 Abstract We give new proofs and explain the origin of several combinatorial identities of Fu and Lascoux, Dilcher,

More information

Computationally, diagonal matrices are the easiest to work with. With this idea in mind, we introduce similarity:

Computationally, diagonal matrices are the easiest to work with. With this idea in mind, we introduce similarity: Diagonalization We have seen that diagonal and triangular matrices are much easier to work with than are most matrices For example, determinants and eigenvalues are easy to compute, and multiplication

More information

Chapter 2 Spectra of Finite Graphs

Chapter 2 Spectra of Finite Graphs Chapter 2 Spectra of Finite Graphs 2.1 Characteristic Polynomials Let G = (V, E) be a finite graph on n = V vertices. Numbering the vertices, we write down its adjacency matrix in an explicit form of n

More information

Eigenvalues and Eigenvectors A =

Eigenvalues and Eigenvectors A = Eigenvalues and Eigenvectors Definition 0 Let A R n n be an n n real matrix A number λ R is a real eigenvalue of A if there exists a nonzero vector v R n such that A v = λ v The vector v is called an eigenvector

More information

Spectral Theory of Orthogonal Polynomials

Spectral Theory of Orthogonal Polynomials Spectral Theory of Orthogonal Polynomials Barry Simon IBM Professor of Mathematics and Theoretical Physics California Institute of Technology Pasadena, CA, U.S.A. Lecture 1: Introduction and Overview Spectral

More information

Linear ODEs. Existence of solutions to linear IVPs. Resolvent matrix. Autonomous linear systems

Linear ODEs. Existence of solutions to linear IVPs. Resolvent matrix. Autonomous linear systems Linear ODEs p. 1 Linear ODEs Existence of solutions to linear IVPs Resolvent matrix Autonomous linear systems Linear ODEs Definition (Linear ODE) A linear ODE is a differential equation taking the form

More information

Bergman shift operators for Jordan domains

Bergman shift operators for Jordan domains [ 1 ] University of Cyprus Bergman shift operators for Jordan domains Nikos Stylianopoulos University of Cyprus Num An 2012 Fifth Conference in Numerical Analysis University of Ioannina September 2012

More information

The L p -dissipativity of first order partial differential operators

The L p -dissipativity of first order partial differential operators The L p -dissipativity of first order partial differential operators A. Cialdea V. Maz ya n Memory of Vladimir. Smirnov Abstract. We find necessary and sufficient conditions for the L p -dissipativity

More information

Iterative Methods for Ax=b

Iterative Methods for Ax=b 1 FUNDAMENTALS 1 Iterative Methods for Ax=b 1 Fundamentals consider the solution of the set of simultaneous equations Ax = b where A is a square matrix, n n and b is a right hand vector. We write the iterative

More information

Calculus C (ordinary differential equations)

Calculus C (ordinary differential equations) Calculus C (ordinary differential equations) Lesson 9: Matrix exponential of a symmetric matrix Coefficient matrices with a full set of eigenvectors Solving linear ODE s by power series Solutions to linear

More information

LINEAR ALGEBRA MICHAEL PENKAVA

LINEAR ALGEBRA MICHAEL PENKAVA LINEAR ALGEBRA MICHAEL PENKAVA 1. Linear Maps Definition 1.1. If V and W are vector spaces over the same field K, then a map λ : V W is called a linear map if it satisfies the two conditions below: (1)

More information

Scattering matrices and Weyl functions

Scattering matrices and Weyl functions Scattering matrices and Weyl functions Jussi Behrndt a, Mark M. Malamud b, Hagen Neidhardt c April 6, 2006 a)technische Universität Berlin, Institut für Mathematik, Straße des 7. Juni 36, D 0623 Berlin,

More information

Quantum Computing Lecture 2. Review of Linear Algebra

Quantum Computing Lecture 2. Review of Linear Algebra Quantum Computing Lecture 2 Review of Linear Algebra Maris Ozols Linear algebra States of a quantum system form a vector space and their transformations are described by linear operators Vector spaces

More information

In this section again we shall assume that the matrix A is m m, real and symmetric.

In this section again we shall assume that the matrix A is m m, real and symmetric. 84 3. The QR algorithm without shifts See Chapter 28 of the textbook In this section again we shall assume that the matrix A is m m, real and symmetric. 3.1. Simultaneous Iterations algorithm Suppose we

More information

Remarks on Some Basic Hypergeometric Series

Remarks on Some Basic Hypergeometric Series Remarks on Some Basic Hypergeometric Series Changgui Zhang Many results in Mathematical Analysis seem to come from some obvious computations. For a few years, we have been interested in the analytic theory

More information

Commutative Banach algebras 79

Commutative Banach algebras 79 8. Commutative Banach algebras In this chapter, we analyze commutative Banach algebras in greater detail. So we always assume that xy = yx for all x, y A here. Definition 8.1. Let A be a (commutative)

More information

MATH 20F: LINEAR ALGEBRA LECTURE B00 (T. KEMP)

MATH 20F: LINEAR ALGEBRA LECTURE B00 (T. KEMP) MATH 20F: LINEAR ALGEBRA LECTURE B00 (T KEMP) Definition 01 If T (x) = Ax is a linear transformation from R n to R m then Nul (T ) = {x R n : T (x) = 0} = Nul (A) Ran (T ) = {Ax R m : x R n } = {b R m

More information

Chapter 3. Vector spaces

Chapter 3. Vector spaces Chapter 3. Vector spaces Lecture notes for MA1111 P. Karageorgis pete@maths.tcd.ie 1/22 Linear combinations Suppose that v 1,v 2,...,v n and v are vectors in R m. Definition 3.1 Linear combination We say

More information

Eigenvalues, Eigenvectors, and an Intro to PCA

Eigenvalues, Eigenvectors, and an Intro to PCA Eigenvalues, Eigenvectors, and an Intro to PCA Eigenvalues, Eigenvectors, and an Intro to PCA Changing Basis We ve talked so far about re-writing our data using a new set of variables, or a new basis.

More information

Fredholm determinant with the confluent hypergeometric kernel

Fredholm determinant with the confluent hypergeometric kernel Fredholm determinant with the confluent hypergeometric kernel J. Vasylevska joint work with I. Krasovsky Brunel University Dec 19, 2008 Problem Statement Let K s be the integral operator acting on L 2

More information

OWN ZEROS LUIS DANIEL ABREU

OWN ZEROS LUIS DANIEL ABREU Pré-Publicações do Departamento de Matemática Universidade de Coimbra Preprint Number 4 3 FUNCTIONS q-orthogonal WITH RESPECT TO THEIR OWN ZEROS LUIS DANIEL ABREU Abstract: In [4], G. H. Hardy proved that,

More information

. The following is a 3 3 orthogonal matrix: 2/3 1/3 2/3 2/3 2/3 1/3 1/3 2/3 2/3

. The following is a 3 3 orthogonal matrix: 2/3 1/3 2/3 2/3 2/3 1/3 1/3 2/3 2/3 Lecture Notes: Orthogonal and Symmetric Matrices Yufei Tao Department of Computer Science and Engineering Chinese University of Hong Kong taoyf@cse.cuhk.edu.hk Orthogonal Matrix Definition. An n n matrix

More information

Explicit inverse of a tridiagonal (p, r) Toeplitz matrix

Explicit inverse of a tridiagonal (p, r) Toeplitz matrix Explicit inverse of a tridiagonal (p, r Toeplitz matrix A.M. Encinas, M.J. Jiménez Dpt. Matemàtiques, Universitat Politècnica de Catalunya - BarcelonaTech Abstract Tridiagonal matrices appears in many

More information

Problems Session. Nikos Stylianopoulos University of Cyprus

Problems Session. Nikos Stylianopoulos University of Cyprus [ 1 ] University of Cyprus Problems Session Nikos Stylianopoulos University of Cyprus Hausdorff Geometry Of Polynomials And Polynomial Sequences Institut Mittag-Leffler Djursholm, Sweden May-June 2018

More information

Review and problem list for Applied Math I

Review and problem list for Applied Math I Review and problem list for Applied Math I (This is a first version of a serious review sheet; it may contain errors and it certainly omits a number of topic which were covered in the course. Let me know

More information

Electromagnetism HW 1 math review

Electromagnetism HW 1 math review Electromagnetism HW math review Problems -5 due Mon 7th Sep, 6- due Mon 4th Sep Exercise. The Levi-Civita symbol, ɛ ijk, also known as the completely antisymmetric rank-3 tensor, has the following properties:

More information

Lecture 10: Solving the Time-Independent Schrödinger Equation. 1 Stationary States 1. 2 Solving for Energy Eigenstates 3

Lecture 10: Solving the Time-Independent Schrödinger Equation. 1 Stationary States 1. 2 Solving for Energy Eigenstates 3 Contents Lecture 1: Solving the Time-Independent Schrödinger Equation B. Zwiebach March 14, 16 1 Stationary States 1 Solving for Energy Eigenstates 3 3 Free particle on a circle. 6 1 Stationary States

More information

DRIFT OF SPECTRALLY STABLE SHIFTED STATES ON STAR GRAPHS

DRIFT OF SPECTRALLY STABLE SHIFTED STATES ON STAR GRAPHS DRIFT OF SPECTRALLY STABLE SHIFTED STATES ON STAR GRAPHS ADILBEK KAIRZHAN, DMITRY E. PELINOVSKY, AND ROY H. GOODMAN Abstract. When the coefficients of the cubic terms match the coefficients in the boundary

More information

On Tricomi and Hermite-Tricomi Matrix Functions of Complex Variable

On Tricomi and Hermite-Tricomi Matrix Functions of Complex Variable Communications in Mathematics and Applications Volume (0), Numbers -3, pp. 97 09 RGN Publications http://www.rgnpublications.com On Tricomi and Hermite-Tricomi Matrix Functions of Complex Variable A. Shehata

More information

Applications of Diagonalization

Applications of Diagonalization Applications of Diagonalization Hsiu-Hau Lin hsiuhau@phys.nthu.edu.tw Apr 2, 200 The notes cover applications of matrix diagonalization Boas 3.2. Quadratic curves Consider the quadratic curve, 5x 2 4xy

More information

ECE 275A Homework #3 Solutions

ECE 275A Homework #3 Solutions ECE 75A Homework #3 Solutions. Proof of (a). Obviously Ax = 0 y, Ax = 0 for all y. To show sufficiency, note that if y, Ax = 0 for all y, then it must certainly be true for the particular value of y =

More information

Lecture 7: Positive Semidefinite Matrices

Lecture 7: Positive Semidefinite Matrices Lecture 7: Positive Semidefinite Matrices Rajat Mittal IIT Kanpur The main aim of this lecture note is to prepare your background for semidefinite programming. We have already seen some linear algebra.

More information