From continuous to the discrete Fourier transform: classical and q
|
|
- Buddy Sharp
- 6 years ago
- Views:
Transcription
1 From continuous to the discrete Fourier transform: classical and quantum aspects 1 Institut de Mathématiques de Bourgogne, Dijon, France 2 Sankt-Petersburg University of Aerospace Instrumentation (SUAI), St-Petersburg, Russia. s: matveev@u-bourgogne.fr, vladimir.matveev9@gmail.com lecture delivered at the international conference "Recent Advances in quantum Integrable Systems" ( RAQIS 14 ) 1-5 September 2014, Dijon, France September 5-th, 2014
2 1 Some general features of FT and DFT 2 The connection of the Gauss sum with the trace of the DFT matrix 3 The complete systems of the eigenvectors of Φ(n) and the real orthogonal matrices O(n) diagonalizing DFT 4 5 Poisson summation formula and intertwining relation between FT and DFT 6 7 Generating eigenvectors of DFT from absolutely convergent series. 8
3 We explain and generalise certain results by Mehta (1987) relating the eigenfunctions of quantum harmonic oscillator with the eigenvectors of the discrete Fourier transform (DFT). Using far more general construction based on a kind of Poisson summation, we associate the eigenvectors of the DFT with the sums of any absolutely convergent series. This construction in particular provides the overcomplete bases of the DFT eigenvectors by means of Jacobian theta functions or ν-theta functions and leads in a quite natural way to some new identities and addition theorems in the theory of special and q-special functions.
4 References 1.M.L. Mehta "Eigenvalues and Eigenvectors of the Finite Fourier Transform" J.Math.Phys. 28 (4) p (1987) 2.V.B. Matveev "Intertwininig relations between the Fourier transform and discrete Fourier transform, the related functional densities and beyond" Inverse Problems 17, (2001) 3. M Ruzzi Jacobi theta functions and discretee Fourier transforms, J.Math.Phys. 47, (2006) 4. D.V. Chudnovsky, G.V.Chudnovsky T.Morgan "Variance of signals and Their Finite Fourier transforms" in vol. D.Chudnovsky and G. Chudnovsky eds. additive number theory : Festschrift in honour of M.B. Nathanson, Springer Science +Buisness Media,LCC 2010 p R. Malekar and H.Bhate "Discrete Fourier Transform and Jacobi function identities" J.Math. Phys (2010) 6. R.Malekar and H.Bhate " Discrete Fourier transform and Riemann identities for theta functions" Applied Mathematical Letters (2012) 7. M.Atakishiyeva, N. Atakishiev and T. Koornwinder "On a q-extension of Mehta s eigenvectors of the finite Fourier transform for q a root of unity. " arxiv : v2 [math.ca] 21 May M.Atakishieva and N.Atakishiev "On integral and finite Fourier transforms of q-hermite polynomials" (2009 ) 9. D. Mumford "Tata lectures on theta" v.2
5 Spectral theorem for the operator roots of unity Suppose U is a bounded operator in a Hilbert space, U = U n = I, q := e 2iπ n, n Z +. n 1 m=0 q m P m, P m = 1 n 1 q jm U j. n n 1 P j = I, Pj 2 = I, UP j = q j P j, P j P m = 0, m j. j=0 Obviously U k = j=1 q jk P j. j=0
6 Spectral multiplicities m j of the eigenvalue q j of U In a finite dimensional case we have: m j := Tr P j = 1 n n 1 k=0 n 1 f (U) := f (q j )P j. j=1 q kj Tr U k.
7 Fourier Transform (FT) 1 f (k) = F(f )(k) := 2π R f (x)e ikx dx Inversion formula: f (x) = 1 2π R f (k)e ikx dk. Since (F 2 f )(x) = f ( x) F 4 = I, (F 3 f )(x) = (F 1 f )(x) = f ( x), (1) where I is an identity operator.
8 Spectral decomposition of F F = P 0 + ip 1 P 2 ip 3, P j = ( i) jk F k k=0 Therefore any continuous absolutely integrable function f (x) generates the 4 eigenfunctions f j of F: F(f j ) = i j f j, j = 0, 1, 2, 3 f j (x) = f (x) + ( i) j f (x) + ( i) 2j f ( x) + ( i) 3j f ( x). This allows to generate the eigenfunctions of F quite different from the most known ψ n = e x2 2 H n (x), which are also eigenfunctions of the quantum harmonic oscillator H.
9 Example: f := exp( αx 2 + itx) with α > 0, t C F(f )(x) = 1 e (x t)2 4α 2α, F 2 (f )(x) = f ( x), F 3 (f )(x) = e (x+t)2 4α. Hence we can generate the 2-parametric family of eigenfunctions f j (x, α, t) = P j f of the Fourier operator:
10 f 0 = e αx2 +itx + e αx 2 itx + 1 e (x t)2 4α 2α e (x+t)2 4α, ( = 2 e αx2 cos tx + 1 e x2 +t 2 ) 4α cosh xt, 2α 2α f 1 = e αx2 +itx e αx 2 itx i 1 e (x t)2 4α e (x+t)2 4α 2α ( = 2i e αx2 sin tx + 1 e x2 +t 2 ) 4α sinh xt, 2α 2α f 2 = e αx2 +itx + e αx 2 itx 1 e (x t)2 4α + e (x+t)2 4α 2α ( = 2 e αx2 cos tx 1 e x2 +t 2 ) 4α cosh xt, 2α 2α f 3 = e αx2 +itx e αx 2 itx + i 1 e (x t)2 4α ie (x+t)2 4α 2α ( = 2i e αx2 sin tx + 1 e x2 +t 2 ) 4α sinh xt. 2α 2α
11 Since F is linear operator, the linear combinations of the functions above, corresponding to the different values of parameters α and t, are also producing the eigenfunctions. We can even integrate them with respect to t and α with some densities ρ(α, t), with compact support, vanishing for α ɛ, ɛ > 0. In this way we can generate eigenfunctions of F depending on any number of the functional parameters ρ i (α, t).
12 Creation and annihilation operators a := d d x + x, a := d d x + x, F a = ia F, F a = ia For any eigenfunction f S of F, corresponding to the eigenvalue i k, we have : F(a j f ) = i k+j a j f, F((a ) j f ) = i k j (a ) j f.
13 Schur matrix and DFT Φ jk (n) = 1 n q jk, j, k = 0,..., n 1, q = e 2iπ n, Φjk = Φ kj, Φ 1 = Φ. The DFT (discrete Fourier transform ) f f describes the action of Φ on the finite sequences,(vectors), f = (f 0,..., f n 1 ) : fk = 1 n 1 f j e 2πijk n, fk±n = f k n j=0
14 Inversion formula for the DFT and Parceval identity f j = 1 n 1 fk e 2πijk n n k=0 or f = Φ 1 f. n 1 n 1 f j 2 = f j 2. j=0 j=0
15 Square of the DFT matrix Φ 2 = In other words, Φ 2 jm = δ n j,m = Φ 2 mj, (Φ2 f ) j = f n j, Φ 4 = I, Φ 3 = Φ 1 = Φ = Φ. Tr Φ 2 = 3 + ( 1)n 2
16 The quantity G = n Tr Φ = is the famous Gauss sum: 1 n G(n) = Tr F = 1 n 1 n s=0 q s2 = 1 + ( i)n. 1 i
17 Spectral decomposition of the DFT operator Φ = 3 i j p j, p j = 1 ( i) jk Φ k, 4 j=0 k 3 p j = I, p j p k = p j δ jk. j=0 The spectral projectors p j are real symmetric matrices!
18 Generating eigenvectors of DFT from any finite sequence It follows from the spectral decomposition above, that any finite sequence sequence f j, j = 0,..., n 1 generates the eigenvector of Φ : Φv = i k v via a formula v j (k) = f j + ( i) k fj + ( 1) k f n j + ( i) 3k f j. Of course for the sequences f j extended periodically we can replace f n j by f j.
19 Spectral multiplicities of Φ(n) n = 4m + 2 m(1) = m + 1 m(i) = m m( 1) = m + 1 m( i) = m, n = 4m m(1) = m + 1 m(i) = m m( 1) = m m( i) = m 1 n = 4m + 1 m(1) = m + 1 m(i) = m m( 1) = m m( i) = m n = 4m + 3 m(1) = m + 1 m(i) = m + 1 m( 1) = m + 1 m( i) = m
20 Characteristic polynomial P(z) and determinant of Φ k=3 P(z) = (z i k ) m k, det F = k=0 k=3 k=1 i km k = i n(n 1) 2 +(n 1) 2 = i 3n2 5n 2 +1,
21 The explicit formulas for spectral projectors p j given above provide the complete set of eigenvectors of Φ: we have to take m 0 first columns of p 0, m 1 first nontrivial, (since first columns of p 1 and p 3 are formed by zeros), columns of p 1, m 2 first columns of p 2 and m 3 first nontrivial columns of p 3 to get it. The completeness follows from the relation p j = I. The eigenvector v(k, m), equal to the m-th column of p k, has the components v j (m, k) (equal to the matrix elements p mj = p jm of p k ), : 4v j (k, m) = δ jm + ( 1) k δ n j,m + ( i) k qjm 3k q jm + ( i) n n
22 Φ(2) := 1 2 ( ),, F 2 = I, p 1 = p 3 = 0. Orthogonal matrix O 2 reducing F to the form diag[1, 1] = O 1 FO is O 2 =
23 n=4 The DFT matrix Φ is Φ = i 1 i i 1 i Here p 3 = 0, since m( i) = m 3 = 0. O 4 : O 1 4 Φ O 4 = diag[1, 1, 1, i] O 4 =
24 Generalized Poisson summation formula Let E be the space of the continuous functions absolutely integrable together with their Fourier transform and f E. Poisson summation formula: ab r Z f (b(ar + x)) = 2π m Z f (2πm(ab) 1 )e 2πimx a.
25 In order to create the DFT structure in the RHS of the Poisson formula, we specify a, b and x and write m as follows : 2π a = n, x = j, b = j = 0,..., n 1; n r Z m = nr + k, r Z, k = 0,..., n 1. The Poisson formula now takes a form: ( ) [ 2π f (nr + j) = 1 n 1 ( )] 2π f (nr + k) q jk. n n n k=0 r Z
26 Consider a following mapping M : M : E C n, M(f ) := v, v j = r Z f ( ) 2π (nr + j) ; j = 0,..., n 1. n we can interpret Poisson formula as an intertwining relation between FT and DFT. Namely the following statement holds.
27 Theorem The mapping M is intertwining F and Φ : Φ M = M F. The mapping M is surjective, it maps any eigenvector of F to some eigenvector of Φ. Any continuous, absolutely integrable eigenfunction ψ of F generates the eigenvector M(ψ) of Φ, corresponding to the same eigenvalue : F(ψ) = i k ψ Φ M(ψ) = i k M(ψ), k = 0, 1, 2, 3. In particular case, when f is any eigenfunction of H, this connection was found by Mehta (1987)
28 Let m Z g m be any absolutely convergent series and g(x) := m Z g m q mx,, q := exp 2iπ n, f j := r Z g rn+j. Obviously f j = f j+n. n 1 j=0 m Z qjl g nm+j = k Z g k q kl = g(l), l = 0,..., n 1; m Z g nm+j = 1 n n 1 l=0 g(l) q lj, q = e 2πi n.
29 Of course the second formula is immediate consequence of the first obtained by applying the inversion formula for DFT to its LHS. To prove the first formula we set r = nm + j, which implies that j = r nm, q jl = q rl lnm = q rl. Now it is clear that a double sum in the LHS of can be written as a single sum which is exactly the same as the RHS which completes the proof.
30 Theorem Let p Z g p is any absolutely converging series. Then vector v(k) with the components v j (k), v j (k) = m Z (g nm+j +( 1) k g nm j )+ ( i)k n m Z (g m +( 1) k g m )e 2πimj n, is an eigenvector of the DFT : Φv = i k v.
31 We can compare now the explicit expressions for the eigenvectors of Φ coming from spectral theorem with those obtained from the formula above. Assuming that g n are the functions depending on some parameters, we can expect that for the functions represented by the RHS of we can get in this way some functional identities. We will discuss here only the identities for generalized theta functions (ν-theta functions) although there are many other applications.
32 ν-theta functions The function θ(x, τ, ν), which we will shortly call ν-theta function, is defined by the formula θ(x, τ, ν) = m= e πiτm2ν +2πimx, ν Z +, Im τ > 0
33 It is obvious that θ(x, τ, ν) is an entire function of x satisfying the relation θ(x + 1, τ, ν) = θ(x, τ, ν). θ(x, τ, ν) reduces to the usual theta function θ(x, τ) when ν = 1.
34 ν-theta function satisfies the following PDE 2 (2π) 2ν 1 ( 1) ν θ τ = i 2ν θ x 2ν.
35 ν-theta function θ a,b (x, τ, ν) with characterstics [a, b] is introduced by the formulas θ a,b (x, τ, ν) = n Z exp[πiτ(n + a) 2ν + 2πi(n + a)(x + b)]. This function also gives the solution to a same parabolic PDE.
36 In the case ν = 1, θ a,b (x, τ, ν) reduces to the usual theta function with characteristics θ a,b (x, τ), simply connected with θ(x, τ): θ(x, τ) = θ(x + aτ + b, τ) exp[πia 2 τ + 2πia(x + b)].
37 Suppose that ( πiτm 2ν g m = exp n 2ν + 2πimx ). n We obviously have the formula g mn+j = θ j,0(x, τ, ν) n m Z
38 Proposition For any τ with Imτ > 0 vector v(τ, x, k, ν) with the components v j, v j (x, τ, k) = θ j,0(x, τ, ν) + ( 1)k θ j,0(x, τ, ν) ( n ) n + 1 n [( i) k θ j+x n, τ, ν + ( i) 3k θ n 2ν is the eigenvector of the DFT : Φ(n)v(k) = i k v(k). ( x j n, τ n 2ν, ν )]. This result already implies interesting functional relations between the ν-theta functions.
39 Let e k are the eigenvectors of Φ(n) : Φ(n)e k = i k e k. We obviously have the relation: Φ(n)(e 0 + e 1 + e 2 + e 3 ) = e 0 + ie 1 e 2 ie 3 Applying this relation to the eigenvectors of Φ(n) given by the series above we get the following statement:
40 Proposition In the case ν = 1 the result formulated below is well known and describes the connection between two bases in the space of entire functions with prescribed quasi periodicity properties (Mumford, Tata lectures on Thetas). The ν-theta functions ( x + k θ n, τ n 2ν, ν ), θ j/n,0 (x, τ, ν), k, j = 0,..., n 1 are connected by the discrete Fourier transform : θ j/n,0 (x, τ, ν) = 1 n θ ( x+k n n 1, τ n 2ν, ν ) = n 1 k=0 q jk θ ( x+k n, τ, ν ), n 2ν j=0 qjk θ j/n,0 (x, τ, ν).
41 In particular, taking j = 0 k = 0 in the formulas above we we obtain that n 1 θ(x, τ, ν) = 1 n k=0 θ ( x+k n, τ, ν ), n 2ν θ ( x n, τ, ν ) = n 1 n 2ν j=0 θ j/n,0(x, τ, ν)
42 Unitarity of the discrete Fourier transform implies now the following relations between the ν-theta functions (Parseval identity): n 1 ( x + k θ n k=0 ), τ 2 n 2ν, ν n 1 = n k=0 θ k/n,0 (x, τ, ν) 2
43 In the simplest case ν = 1, n = 2 we have two eigenvalues of Φ(2),- λ 1 = 1 and λ 2 = 1, and two eigenvectors v1, v( 1) The generaal formula expressing the components of these eigenvectors by means of Jacobi theta functions reads, (j = 0, 1): v j (1) = θ j/2,0 (x, τ)+θ j/2,0 (x, τ)+ 1 2 [θ ( x + j 2, τ 4 v j ( 1) = θ j/2,0 (x, τ)+θ j/2,0 (x, τ) 1 ( x + j [θ 2 2, τ 4 ) ( x j + θ ) + θ 2, τ 4 ( x j 2, τ 4 )], )]
44 Since Φ(2)(v(1) + v(2)) = v(1) v(2) the formulas above imply two identities: ( x θ(x, τ)+θ 1/2,0 = θ 2 4), τ, θ(x, τ) θ 1/2,0 (x, τ) = θ ( x + 1 2, τ ). 4
45 Suppose that Φv(x, τ) = v(x, τ). v(x + τ/2, τ) and v(y + 1, τ) are again eigenvectors, corresponding to the same eigenvalue λ 1 = 1. Therefore det [v(x + τ/2, τ), v(y + 1, τ)] = 0. This gives a quadratic relation between theta functions allowing to prove extended and classical quartic Riemann identities.
46 THANK YOU FOR YOUR ATTENTION!
1 Assignment 1: Nonlinear dynamics (due September
Assignment : Nonlinear dynamics (due September 4, 28). Consider the ordinary differential equation du/dt = cos(u). Sketch the equilibria and indicate by arrows the increase or decrease of the solutions.
More informationTHE SELBERG TRACE FORMULA OF COMPACT RIEMANN SURFACES
THE SELBERG TRACE FORMULA OF COMPACT RIEMANN SURFACES IGOR PROKHORENKOV 1. Introduction to the Selberg Trace Formula This is a talk about the paper H. P. McKean: Selberg s Trace Formula as applied to a
More informationDunkl operators and Clifford algebras II
Translation operator for the Clifford Research Group Department of Mathematical Analysis Ghent University Hong Kong, March, 2011 Translation operator for the Hermite polynomials Translation operator for
More informationPhysics 137A Quantum Mechanics Fall 2012 Midterm II - Solutions
Physics 37A Quantum Mechanics Fall 0 Midterm II - Solutions These are the solutions to the exam given to Lecture Problem [5 points] Consider a particle with mass m charge q in a simple harmonic oscillator
More informationFinite temperature form factors in the free Majorana theory
Finite temperature form factors in the free Majorana theory Benjamin Doyon Rudolf Peierls Centre for Theoretical Physics, Oxford University, UK supported by EPSRC postdoctoral fellowship hep-th/0506105
More informationConnection Formula for Heine s Hypergeometric Function with q = 1
Connection Formula for Heine s Hypergeometric Function with q = 1 Ryu SASAKI Department of Physics, Shinshu University based on arxiv:1411.307[math-ph], J. Phys. A in 48 (015) 11504, with S. Odake nd Numazu
More informationA Vector-Space Approach for Stochastic Finite Element Analysis
A Vector-Space Approach for Stochastic Finite Element Analysis S Adhikari 1 1 Swansea University, UK CST2010: Valencia, Spain Adhikari (Swansea) Vector-Space Approach for SFEM 14-17 September, 2010 1 /
More information1 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 information2.4.8 Heisenberg Algebra, Fock Space and Harmonic Oscillator
.4. SPECTRAL DECOMPOSITION 63 Let P +, P, P 1,..., P p be the corresponding orthogonal complimentary system of projections, that is, P + + P + p P i = I. i=1 Then there exists a corresponding system of
More informationMathematical Methods wk 2: Linear Operators
John Magorrian, magog@thphysoxacuk These are work-in-progress notes for the second-year course on mathematical methods The most up-to-date version is available from http://www-thphysphysicsoxacuk/people/johnmagorrian/mm
More informationLecture 15 Review of Matrix Theory III. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore
Lecture 15 Review of Matrix Theory III Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Matrix An m n matrix is a rectangular or square array of
More informationShort-time Fourier transform for quaternionic signals
Short-time Fourier transform for quaternionic signals Joint work with Y. Fu and U. Kähler P. Cerejeiras Departamento de Matemática Universidade de Aveiro pceres@ua.pt New Trends and Directions in Harmonic
More informationBaker-Akhiezer functions and configurations of hyperplanes
Baker-Akhiezer functions and configurations of hyperplanes Alexander Veselov, Loughborough University ENIGMA conference on Geometry and Integrability, Obergurgl, December 2008 Plan BA function related
More informationLecture 6 Quantum Mechanical Systems and Measurements
Lecture 6 Quantum Mechanical Systems and Measurements Today s Program: 1. Simple Harmonic Oscillator (SHO). Principle of spectral decomposition. 3. Predicting the results of measurements, fourth postulate
More informationSpectral properties of quantum graphs via pseudo orbits
Spectral properties of quantum graphs via pseudo orbits 1, Ram Band 2, Tori Hudgins 1, Christopher Joyner 3, Mark Sepanski 1 1 Baylor, 2 Technion, 3 Queen Mary University of London Cardiff 12/6/17 Outline
More informationStrauss PDEs 2e: Section Exercise 4 Page 1 of 6
Strauss PDEs 2e: Section 5.3 - Exercise 4 Page of 6 Exercise 4 Consider the problem u t = ku xx for < x < l, with the boundary conditions u(, t) = U, u x (l, t) =, and the initial condition u(x, ) =, where
More informationChapter 2 Heisenberg s Matrix Mechanics
Chapter 2 Heisenberg s Matrix Mechanics Abstract The quantum selection rule and its generalizations are capable of predicting energies of the stationary orbits; however they should be obtained in a more
More informationSymmetries for fun and profit
Symmetries for fun and profit Sourendu Gupta TIFR Graduate School Quantum Mechanics 1 August 28, 2008 Sourendu Gupta (TIFR Graduate School) Symmetries for fun and profit QM I 1 / 20 Outline 1 The isotropic
More informationx 3y 2z = 6 1.2) 2x 4y 3z = 8 3x + 6y + 8z = 5 x + 3y 2z + 5t = 4 1.5) 2x + 8y z + 9t = 9 3x + 5y 12z + 17t = 7
Linear Algebra and its Applications-Lab 1 1) Use Gaussian elimination to solve the following systems x 1 + x 2 2x 3 + 4x 4 = 5 1.1) 2x 1 + 2x 2 3x 3 + x 4 = 3 3x 1 + 3x 2 4x 3 2x 4 = 1 x + y + 2z = 4 1.4)
More information2. The Schrödinger equation for one-particle problems. 5. Atoms and the periodic table of chemical elements
1 Historical introduction The Schrödinger equation for one-particle problems 3 Mathematical tools for quantum chemistry 4 The postulates of quantum mechanics 5 Atoms and the periodic table of chemical
More informationREPRESENTATION THEORY WEEK 7
REPRESENTATION THEORY WEEK 7 1. Characters of L k and S n A character of an irreducible representation of L k is a polynomial function constant on every conjugacy class. Since the set of diagonalizable
More informationWave packet decompositions adapted to (non-self-adjoint) operators
1 / 42 Wave packet decompositions adapted to (non-self-adjoint) operators Joe Viola Laboratoire de Mathématiques Jean Leray Université de Nantes Joseph.Viola@univ-nantes.fr 26 March 2018 2 / 42 Outline
More informationPHY 396 K. Solutions for problems 1 and 2 of set #5.
PHY 396 K. Solutions for problems 1 and of set #5. Problem 1a: The conjugacy relations  k,  k,, Ê k, Ê k, follow from hermiticity of the Âx and Êx quantum fields and from the third eq. 6 for the polarization
More informationFORMULA SHEET FOR QUIZ 2 Exam Date: November 8, 2017
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Electromagnetism II November 5, 207 Prof. Alan Guth FORMULA SHEET FOR QUIZ 2 Exam Date: November 8, 207 A few items below are marked
More informationRepresentation theory and quantum mechanics tutorial Spin and the hydrogen atom
Representation theory and quantum mechanics tutorial Spin and the hydrogen atom Justin Campbell August 3, 2017 1 Representations of SU 2 and SO 3 (R) 1.1 The following observation is long overdue. Proposition
More informationarxiv:solv-int/ v1 4 Sep 1995
Hidden symmetry of the quantum Calogero-Moser system Vadim B. Kuznetsov Institute of Mathematical Modelling 1,2,3 Technical University of Denmark, DK-2800 Lyngby, Denmark arxiv:solv-int/9509001v1 4 Sep
More informationThe Density Matrix for the Ground State of 1-d Impenetrable Bosons in a Harmonic Trap
The Density Matrix for the Ground State of 1-d Impenetrable Bosons in a Harmonic Trap Institute of Fundamental Sciences Massey University New Zealand 29 August 2017 A. A. Kapaev Memorial Workshop Michigan
More information1. 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 informationMP 472 Quantum Information and Computation
MP 472 Quantum Information and Computation http://www.thphys.may.ie/staff/jvala/mp472.htm Outline Open quantum systems The density operator ensemble of quantum states general properties the reduced density
More informationQMI PRELIM Problem 1. All problems have the same point value. If a problem is divided in parts, each part has equal value. Show all your work.
QMI PRELIM 013 All problems have the same point value. If a problem is divided in parts, each part has equal value. Show all your work. Problem 1 L = r p, p = i h ( ) (a) Show that L z = i h y x ; (cyclic
More informationPhysics 342 Lecture 26. Angular Momentum. Lecture 26. Physics 342 Quantum Mechanics I
Physics 342 Lecture 26 Angular Momentum Lecture 26 Physics 342 Quantum Mechanics I Friday, April 2nd, 2010 We know how to obtain the energy of Hydrogen using the Hamiltonian operator but given a particular
More informationLecture II Search Method for Lax Pairs of Nonlinear Partial Differential Equations
Lecture II Search Method for Lax Pairs of Nonlinear Partial Differential Equations Usama Al Khawaja, Department of Physics, UAE University, 24 Jan. 2012 First International Winter School on Quantum Gases
More informationMathematical Preliminaries
Mathematical Preliminaries Introductory Course on Multiphysics Modelling TOMAZ G. ZIELIŃKI bluebox.ippt.pan.pl/ tzielins/ Table of Contents Vectors, tensors, and index notation. Generalization of the concept
More informationMatrix Theory. A.Holst, V.Ufnarovski
Matrix Theory AHolst, VUfnarovski 55 HINTS AND ANSWERS 9 55 Hints and answers There are two different approaches In the first one write A as a block of rows and note that in B = E ij A all rows different
More informationFunctional 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 informationGeneralized Burgers equations and Miura Map in nonabelian ring. nonabelian rings as integrable systems.
Generalized Burgers equations and Miura Map in nonabelian rings as integrable systems. Sergey Leble Gdansk University of Technology 05.07.2015 Table of contents 1 Introduction: general remarks 2 Remainders
More informationTalk 3: Examples: Star exponential
Akira Yoshioka 2 7 February 2013 Shinshu Akira Yoshioka 3 2 Dept. Math. Tokyo University of Science 1. Star exponential 1 We recall the definition of star exponentials. 2 We show some construction of star
More informationUniversality of distribution functions in random matrix theory Arno Kuijlaars Katholieke Universiteit Leuven, Belgium
Universality of distribution functions in random matrix theory Arno Kuijlaars Katholieke Universiteit Leuven, Belgium SEA 06@MIT, Workshop on Stochastic Eigen-Analysis and its Applications, MIT, Cambridge,
More informationlinearly indepedent eigenvectors as the multiplicity of the root, but in general there may be no more than one. For further discussion, assume matrice
3. Eigenvalues and Eigenvectors, Spectral Representation 3.. Eigenvalues and Eigenvectors A vector ' is eigenvector of a matrix K, if K' is parallel to ' and ' 6, i.e., K' k' k is the eigenvalue. If is
More informationProducts of random matrices
PHYSICAL REVIEW E 66, 664 Products of random matrices A. D. Jackson* and B. Lautrup The Niels Bohr Institute, Copenhagen, Denmark P. Johansen The Institute of Computer Science, University of Copenhagen,
More information11 a 12 a 21 a 11 a 22 a 12 a 21. (C.11) A = The determinant of a product of two matrices is given by AB = A B 1 1 = (C.13) and similarly.
C PROPERTIES OF MATRICES 697 to whether the permutation i 1 i 2 i N is even or odd, respectively Note that I =1 Thus, for a 2 2 matrix, the determinant takes the form A = a 11 a 12 = a a 21 a 11 a 22 a
More informationConformal maps. Lent 2019 COMPLEX METHODS G. Taylor. A star means optional and not necessarily harder.
Lent 29 COMPLEX METHODS G. Taylor A star means optional and not necessarily harder. Conformal maps. (i) Let f(z) = az + b, with ad bc. Where in C is f conformal? cz + d (ii) Let f(z) = z +. What are the
More informationPartial Differential Equations
M3M3 Partial Differential Equations Solutions to problem sheet 3/4 1* (i) Show that the second order linear differential operators L and M, defined in some domain Ω R n, and given by Mφ = Lφ = j=1 j=1
More information(super)algebras and its applications
for Lie (super)algebras and its applications Institute of Nuclear Physics, Moscow State University, 119 992 Moscow, Russia The International School Advanced Methods of Modern Theoretical Physics: Integrable
More informationQuantum Field Theory III
Quantum Field Theory III Prof. Erick Weinberg January 19, 2011 1 Lecture 1 1.1 Structure We will start with a bit of group theory, and we will talk about spontaneous symmetry broken. Then we will talk
More information4.3 Classical model at low temperature: an example Continuum limit and path integral The two-point function: perturbative
Contents Gaussian integrals.......................... Generating function.......................2 Gaussian expectation values. Wick s theorem........... 2.3 Perturbed gaussian measure. Connected contributions.......
More information2-Dimensional Density Matrices
A Transformational Property of -Dimensional Density Matrices Nicholas Wheeler January 03 Introduction. This note derives from recent work related to decoherence models described by Zurek and Albrecht.
More informationRandom matrices and the Riemann zeros
Random matrices and the Riemann zeros Bruce Bartlett Talk given in Postgraduate Seminar on 5th March 2009 Abstract Random matrices and the Riemann zeta function came together during a chance meeting between
More informationLecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3.12 in Boas)
Lecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3 in Boas) As suggested in Lecture 8 the formalism of eigenvalues/eigenvectors has many applications in physics, especially in
More informationThe elliptic sinh-gordon equation in the half plane
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 8 25), 63 73 Research Article The elliptic sinh-gordon equation in the half plane Guenbo Hwang Department of Mathematics, Daegu University, Gyeongsan
More informationAsymptotics of Hermite polynomials
Asymptotics of Hermite polynomials Michael Lindsey We motivate the study of the asymptotics of Hermite polynomials via their appearance in the analysis of the Gaussian Unitary Ensemble (GUE). Following
More information(V.C) Complex eigenvalues and other topics
V.C Complex eigenvalues and other topics matrix Let s take a closer look at the characteristic polynomial for a 2 2 namely A = a c f A λ = detλi A = λ aλ d bc b d = λ 2 a + dλ + ad bc = λ 2 traλ + det
More informationSlepian Functions Why, What, How?
Slepian Functions Why, What, How? Lecture 1: Basics of Fourier and Fourier-Legendre Transforms Lecture : Spatiospectral Concentration Problem Cartesian & Spherical Cases Lecture 3: Slepian Functions and
More informationQuantum Physics II (8.05) Fall 2004 Assignment 3
Quantum Physics II (8.5) Fall 24 Assignment 3 Massachusetts Institute of Technology Physics Department Due September 3, 24 September 23, 24 7:pm This week we continue to study the basic principles of quantum
More informationMP463 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 informationProblem set 2 The central limit theorem.
Problem set 2 The central limit theorem. Math 22a September 6, 204 Due Sept. 23 The purpose of this problem set is to walk through the proof of the central limit theorem of probability theory. Roughly
More informationQuantum 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 informationEXERCISES IN MODULAR FORMS I (MATH 726) (2) Prove that a lattice L is integral if and only if its Gram matrix has integer coefficients.
EXERCISES IN MODULAR FORMS I (MATH 726) EYAL GOREN, MCGILL UNIVERSITY, FALL 2007 (1) We define a (full) lattice L in R n to be a discrete subgroup of R n that contains a basis for R n. Prove that L is
More informationRepresentation of objects. Representation of objects. Transformations. Hierarchy of spaces. Dimensionality reduction
Representation of objects Representation of objects For searching/mining, first need to represent the objects Images/videos: MPEG features Graphs: Matrix Text document: tf/idf vector, bag of words Kernels
More informationThe spectral zeta function
The spectral zeta function Bernd Ammann June 4, 215 Abstract In this talk we introduce spectral zeta functions. The spectral zeta function of the Laplace-Beltrami operator was already introduced by Minakshisundaram
More informationLinear Algebra: Matrix Eigenvalue Problems
CHAPTER8 Linear Algebra: Matrix Eigenvalue Problems Chapter 8 p1 A matrix eigenvalue problem considers the vector equation (1) Ax = λx. 8.0 Linear Algebra: Matrix Eigenvalue Problems Here A is a given
More informationTopics in Harmonic Analysis Lecture 1: The Fourier transform
Topics in Harmonic Analysis Lecture 1: The Fourier transform Po-Lam Yung The Chinese University of Hong Kong Outline Fourier series on T: L 2 theory Convolutions The Dirichlet and Fejer kernels Pointwise
More informationMathematical Physics Homework 10
Georgia Institute of Technology Mathematical Physics Homework Conner Herndon November, 5 Several types of orthogonal polynomials frequently occur in various physics problems. For instance, Hermite polynomials
More informationEigenvalue PDFs. Peter Forrester, M&S, University of Melbourne
Outline Eigenvalue PDFs Peter Forrester, M&S, University of Melbourne Hermitian matrices with real, complex or real quaternion elements Circular ensembles and classical groups Products of random matrices
More information2.1 Calculation of the ground state energy via path integral
Chapter 2 Instantons in Quantum Mechanics Before describing the instantons and their effects in non-abelian gauge theories, it is instructive to become familiar with the relevant ideas in the context of
More informationPhysics 221A Fall 1996 Notes 16 Bloch s Theorem and Band Structure in One Dimension
Physics 221A Fall 1996 Notes 16 Bloch s Theorem and Band Structure in One Dimension In these notes we examine Bloch s theorem and band structure in problems with periodic potentials, as a part of our survey
More informationLucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche
Scuola di Dottorato THE WAVE EQUATION Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche Lucio Demeio - DIISM wave equation 1 / 44 1 The Vibrating String Equation 2 Second
More information2 Quantization of the Electromagnetic Field
2 Quantization of the Electromagnetic Field 2.1 Basics Starting point of the quantization of the electromagnetic field are Maxwell s equations in the vacuum (source free): where B = µ 0 H, D = ε 0 E, µ
More informationFinite-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 informationPART IV Spectral Methods
PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,
More informationStatistical Interpretation
Physics 342 Lecture 15 Statistical Interpretation Lecture 15 Physics 342 Quantum Mechanics I Friday, February 29th, 2008 Quantum mechanics is a theory of probability densities given that we now have an
More informationDEN: Linear algebra numerical view (GEM: Gauss elimination method for reducing a full rank matrix to upper-triangular
form) Given: matrix C = (c i,j ) n,m i,j=1 ODE and num math: Linear algebra (N) [lectures] c phabala 2016 DEN: Linear algebra numerical view (GEM: Gauss elimination method for reducing a full rank matrix
More informationEigenmodes for coupled harmonic vibrations. Algebraic Method for Harmonic Oscillator.
PHYS208 spring 2008 Eigenmodes for coupled harmonic vibrations. Algebraic Method for Harmonic Oscillator. 07.02.2008 Adapted from the text Light - Atom Interaction PHYS261 autumn 2007 Go to list of topics
More informationCollocation based high dimensional model representation for stochastic partial differential equations
Collocation based high dimensional model representation for stochastic partial differential equations S Adhikari 1 1 Swansea University, UK ECCM 2010: IV European Conference on Computational Mechanics,
More informationGeometric control and dynamical systems
Université de Nice - Sophia Antipolis & Institut Universitaire de France 9th AIMS International Conference on Dynamical Systems, Differential Equations and Applications Control of an inverted pendulum
More informationSymmetry and Properties of Crystals (MSE638) Stress and Strain Tensor
Symmetry and Properties of Crystals (MSE638) Stress and Strain Tensor Somnath Bhowmick Materials Science and Engineering, IIT Kanpur April 6, 2018 Tensile test and Hooke s Law Upto certain strain (0.75),
More informationConsider 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 informationThe Kernel Trick. Robert M. Haralick. Computer Science, Graduate Center City University of New York
The Kernel Trick Robert M. Haralick Computer Science, Graduate Center City University of New York Outline SVM Classification < (x 1, c 1 ),..., (x Z, c Z ) > is the training data c 1,..., c Z { 1, 1} specifies
More informationEigenvalues and Eigenvectors
/88 Chia-Ping Chen Department of Computer Science and Engineering National Sun Yat-sen University Linear Algebra Eigenvalue Problem /88 Eigenvalue Equation By definition, the eigenvalue equation for matrix
More informationPROBABILITY DISTRIBUTIONS AND ORTHOGONAL POLYNOMIALS ASSOCIATED WITH THE ONE-PARAMETER FIBONACCI GROUP. 1. Introduction
Communications on Stochastic Analysis Vol. 10, No. 1 (2016) 117-130 Serials Publications www.serialspublications.com PROBABILITY DISTRIBUTIONS AND ORTHOGONAL POLYNOMIALS ASSOCIATED WITH THE ONE-PARAMETER
More informationIntroduction to relativistic quantum mechanics
Introduction to relativistic quantum mechanics. Tensor notation In this book, we will most often use so-called natural units, which means that we have set c = and =. Furthermore, a general 4-vector will
More informationQuantum mechanics in one hour
Chapter 2 Quantum mechanics in one hour 2.1 Introduction The purpose of this chapter is to refresh your knowledge of quantum mechanics and to establish notation. Depending on your background you might
More informationQuantum chaos on graphs
Baylor University Graduate seminar 6th November 07 Outline 1 What is quantum? 2 Everything you always wanted to know about quantum but were afraid to ask. 3 The trace formula. 4 The of Bohigas, Giannoni
More informationLecture 12: Detailed balance and Eigenfunction methods
Miranda Holmes-Cerfon Applied Stochastic Analysis, Spring 2015 Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility),
More informationHilbert Space Problems
Hilbert Space Problems Prescribed books for problems. ) Hilbert Spaces, Wavelets, Generalized Functions and Modern Quantum Mechanics by Willi-Hans Steeb Kluwer Academic Publishers, 998 ISBN -7923-523-9
More informationClass notes: Approximation
Class notes: Approximation Introduction Vector spaces, linear independence, subspace The goal of Numerical Analysis is to compute approximations We want to approximate eg numbers in R or C vectors in R
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More informationPLEASE LET ME KNOW IF YOU FIND TYPOS (send to
Teoretisk Fysik KTH Advanced QM (SI2380), Lecture 2 (Summary of concepts) 1 PLEASE LET ME KNOW IF YOU FIND TYPOS (send email to langmann@kth.se) The laws of QM 1. I now discuss the laws of QM and their
More informationTopics in Fourier analysis - Lecture 2.
Topics in Fourier analysis - Lecture 2. Akos Magyar 1 Infinite Fourier series. In this section we develop the basic theory of Fourier series of periodic functions of one variable, but only to the extent
More informationA complete criterion for convex-gaussian states detection
A complete criterion for convex-gaussian states detection Anna Vershynina Institute for Quantum Information, RWTH Aachen, Germany joint work with B. Terhal NSF/CBMS conference Quantum Spin Systems The
More informationPH.D. PRELIMINARY EXAMINATION MATHEMATICS
UNIVERSITY OF CALIFORNIA, BERKELEY SPRING SEMESTER 207 Dept. of Civil and Environmental Engineering Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem
More informationGeneral Exam Part II, Fall 1998 Quantum Mechanics Solutions
General Exam Part II, Fall 1998 Quantum Mechanics Solutions Leo C. Stein Problem 1 Consider a particle of charge q and mass m confined to the x-y plane and subject to a harmonic oscillator potential V
More informationarxiv: v1 [hep-ph] 30 Dec 2015
June 3, 8 Derivation of functional equations for Feynman integrals from algebraic relations arxiv:5.94v [hep-ph] 3 Dec 5 O.V. Tarasov II. Institut für Theoretische Physik, Universität Hamburg, Luruper
More informationFrom the mesoscopic to microscopic scale in random matrix theory
From the mesoscopic to microscopic scale in random matrix theory (fixed energy universality for random spectra) With L. Erdős, H.-T. Yau, J. Yin Introduction A spacially confined quantum mechanical system
More informationAlgebra Workshops 10 and 11
Algebra Workshops 1 and 11 Suggestion: For Workshop 1 please do questions 2,3 and 14. For the other questions, it s best to wait till the material is covered in lectures. Bilinear and Quadratic Forms on
More informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
More informationLarge Deviations, Linear Statistics, and Scaling Limits for Mahler Ensemble of Complex Random Polynomials
Large Deviations, Linear Statistics, and Scaling Limits for Mahler Ensemble of Complex Random Polynomials Maxim L. Yattselev joint work with Christopher D. Sinclair International Conference on Approximation
More informationOPSF, Random Matrices and Riemann-Hilbert problems
OPSF, Random Matrices and Riemann-Hilbert problems School on Orthogonal Polynomials in Approximation Theory and Mathematical Physics, ICMAT 23 27 October, 2017 Plan of the course lecture 1: Orthogonal
More informationPartial differential equations (ACM30220)
(ACM3. A pot on a stove has a handle of length that can be modelled as a rod with diffusion constant D. The equation for the temperature in the rod is u t Du xx < x
More informationNon-perturbative effects in ABJM theory
October 9, 2015 Non-perturbative effects in ABJM theory 1 / 43 Outline 1. Non-perturbative effects 1.1 General aspects 1.2 Non-perturbative aspects of string theory 1.3 Non-perturbative effects in M-theory
More information