The Density Matrix for the Ground State of 1-d Impenetrable Bosons in a Harmonic Trap
|
|
- Marsha Margaret Sims
- 6 years ago
- Views:
Transcription
1 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 Center for Applied Interdisciplinary Mathematics August 25-29, 2017
2 The Plan 1-d Impenetrable Bosons
3 The Plan 1-d Impenetrable Bosons Density Matrix and Bose Condensation
4 The Plan 1-d Impenetrable Bosons Density Matrix and Bose Condensation Homogeneous and Periodic Case
5 The Plan 1-d Impenetrable Bosons Density Matrix and Bose Condensation Homogeneous and Periodic Case Recurrence Relations
6 The Plan 1-d Impenetrable Bosons Density Matrix and Bose Condensation Homogeneous and Periodic Case Recurrence Relations Orthogonal Polynomial System & Semi-Classical OPS Class
7 The Plan 1-d Impenetrable Bosons Density Matrix and Bose Condensation Homogeneous and Periodic Case Recurrence Relations Orthogonal Polynomial System & Semi-Classical OPS Class Lax pairs
8 The Plan 1-d Impenetrable Bosons Density Matrix and Bose Condensation Homogeneous and Periodic Case Recurrence Relations Orthogonal Polynomial System & Semi-Classical OPS Class Lax pairs L(1, 1, 3; 2) member of degenerations from 2-variable Garnier system
9 1-d Impenetrable Bosons Hamiltonian N bosons x j (, ) N 2 H = x 2 i=1 i N + x 2 j + 2c δ(x i x j ), c > 0 j=1 1 i<j N Ground State eigenfunction Ψ 0 (x 1,..., x N ) satisfies subject to boundary conditions 1. Ψ 0 xi =x j +0 = Ψ 0 xi =x j 0 ( Ψ0 2. Ψ ) 0 x i x j xi =x j +0 HΨ 0 = E 0 Ψ 0 ( Ψ0 Ψ ) 0 = 2c Ψ x i x j 0 xi =x j ±0 xi =x j 0 Bethe Ansatz solution without harmonic trap [Lieb+Liniger 1963] Impenetrable Bose wavefunction c [Girardeau 1960]. The wave function for impenetrable bosons must vanish at coincident coordinates, Ψ 0 (x 1,..., x i,..., x j,..., x N ) = 0 for x i = x j, (i j), This means that for any fixed ordering of the particles x 1 < x 2 <... < x N there is no distinction between impenetrable bosons and free fermions.
10 Ground State Wavefunction Schrödinger Operator N 2 H R = x 2 j=1 j N + x 2 j j=1 Normalised single particle eigenstates {φ k (x)} k=0,1,... φ k (x) = 2 k c k e x2 /2 H k (x), (c k ) 2 = π 1/2 2 k k! H k (x) denotes the Hermite polynomial of degree k. Ground state wave function Ψ 0 as a Slater determinant of distinct single-particle states 1 Ψ 0 (x 1,..., x N ) = (N!) 1/2 det[φ k(x j )] j=1,...,n k=0,...,n 1 where the factor of (N!) 1/2 is included so that dx 1 dx N Ψ 0 (x 1,..., x N ) 2 = 1.
11 Density Matrix Density Matrix is defined as ρ N (x; y) N dx 2... dx N Ψ 0 (x, x 2,..., x N )Ψ 0 (y, x 2,..., x N ), R R so that dx ρ N (x; x) = N R Using this one must solve the eigenvalue problem ρ N (x; y)ψ k (y) dy = λ k ψ k (x), k Z 0. R Note λ 0 > λ 1 >... > λ N 1 > 0, N 1 λ j = N j=0
12 GUE Averages Recall ρ N (x; y) = N dx 2... dx N Ψ 0 (x, x 2,..., x N )Ψ 0 (y, x 2,..., x N ) R R Vandermonde determinant formula with p j (x) a monic polynomial of degree j, det[p j 1 (x k )] j,k=1,...,n = det[x j 1 ] k j,k=1,...,n = (x k x j ) 1 j<k N Product formula for Ψ 0 Ψ 0 (x 1,..., x N ) = 1 C N N j=1 e x2 j /2 1 j<k N N 1 x k x j, (C N ) 2 = N! (c l ) 2. This coincides precisely with the multivariate probability distribution function (p.d.f.) for a particular class of random matrices l=0 Ψ 0 2 = E GUEN namely, GUE N, the Gaussian unitary ensemble of N N complex Hermitian matrices, Density matrix ρ N+1 (x; y) can be written as an average over this eigenvalue p.d.f ρ N+1 (x; y) = 1 N c 2 e x2 /2 y 2 /2 E GUEN x x l y x l N l=1
13 Density Matrix and Bose Condensation Criteria for Bose-Einstein condensation (BEC) are - Onsager-Penrose 1956 criteria for BEC λ 0 (N) N = e O(1) BEC o(1) no BEC as N Off Diagonal Long Range Order criteria of Yang 1962 ρ (x) x C > 0 NB: already in thermodynamic limit N
14 Thermodynamic limit of Homogeneous System [Jimbo Miwa Mori Sato 1980] Density Matrix ρ (x; 0) is now identified with a τ-function of the fifth Painlevé system ρ (x; 0) = ρ 0 exp ( πρ 0 x σ V (t) dt ), 0 t where ρ 0 denotes the bulk density. σ V satisfies the non-linear equation subject to the x 0 boundary condition (xσ V )2 + 4(xσ V σ V 1) ( xσ V σ V + (σ V )2) = 0 σ V (x) x 0 x2 3 + x3 3π + O(x4 ). [Lenard 1972] Eigenvalues of the density matrix are Fourier coefficients L λ k = ρ N (x; 0)e 2πikx/L dx. 0 The N asymptotic expansion ( ) 1/2 π ρ N (x; 0) ρ 0 A, A = G4 (3/2) N sin(πρ 0 x/n) 2π where G(x) denotes the Barnes G-function, valid for x/n fixed. Therefore a specific c computable from the above. λ 0 c N
15 Finite N Periodic System [FFGW 2003] Finite system of N + 1 particles on a ring of circumference L, x [0, L]. Wavefunction 1 Ψ 0 (x 1,..., x N ) = L N/2 (N!) 1/2 2 sin π(x k x j ) L 1 j<k N Density Matrix ρ N+1 (x; 0) ρ N+1 (x; 0) = 1 N L E ( πx U(N) 2 sin L θ ) l ( ) θl 2 sin 2 2 l=1 Density Matrix can be identified with a member of a τ-function sequence in the sixth Painlevé system with Okamoto parameters b = ( 1 2 N, N, 1 2 N, N). Thermodynamic limit made by substituting x = e 2it/N and taking N with N/L fixed.
16 Determinants The general Heine identity N! det [ Determinant formula with matrix elements ] dx w(x)x j+k 2 = R j,k=1,...,n dx 1 R dx N R N w(x l ) l=1 ρ N+1 (x; y) = 1 c 2 e x2 /2 y 2 /2 det[a j,k (x; y)] j,k=1,...,n N a j,k (x; y) = 2 j k+2 c j 1 c k 1 To simplify further, we note 1 j<k N dt x t y t H j 1 (t)h k 1 (t)e t2 (x t)(y t), t [x, y] x t y t = (x t)(y t), t [x, y] x j x k 2
17 Orthogonal Polynomial System (OPS) Orthogonal polynomial system {p n (t)} defined by the orthogonality relations n=0 dt w(t)p n (t)t m 0 0 m < n = R h n m = n A consequence of the orthogonality relation is the three term recurrence relation a n+1 p n+1 (t) = (t b n )p n (t) a n p n 1 (t), n 1 Central objects in our probabilistic model are the Hankel determinants n := det[µ j+k 2 ] j,k=1,...,n, n 1, 0 := 1, in terms of the moments of the weight {µ n } n=0,1,..., µ n := dt w(t)t n R Alternatively n = 1 n! dt 1... R dt n R result in orthogonal polynomial theory n w(t l ) l=1 1 j<k n a 2 n = n+1 n 1 2, n 1, n (t k t j ) 2, n 1
18 Spectral Derivative/1st Lax pair [Bauldry 1990, Bonan and Clark 1990] Orthogonal polynomials satisfy a system of coupled first order linear differential equations with respect to t ( d/dt) Wp n = (Ω n V)p n a n Θ n p n 1, n 1, with the coefficient functions V(t), W(t), Θ n (t), Ω n (t). Compatibility Condition Coefficient functions satisfy the recurrence relations Ω n+1 + Ω n = (t b n )Θ n, n 0, (Ω n+1 Ω n )(t b n ) = W + a 2 n+1 Θ n+1 a 2 nθ n 1, n 0
19 Semi-Classical OPS Class [Maroni, > ] Semi-Classical weights defined by 1 d 2V(t) w(t) = w(t) dt W(t) with V(t) and W(t) irreducible polynomials in t. Deformed Hermite Orthogonal Polynomials with weight w(t; x, y) = t x t y e t2, t R, = [1 2χ (x,y) (t)](t x)(t y)e t2 where 1 if t I χ I (t) = 0 otherwise
20 Spectral Data [Magnus 1995,...] Numerator and denominator polynomials W(t) = (t x)(t y), 2V(t) = 2t 3 + 2(x + y)t 2 + 2t(1 xy) x y N.B. we have finite singularities at t = x, y. Coefficient functions are and Θ n (t) = 2t 2 + 2(x + y b n )t Ω n (t) = t 3 + (x + y)t 2 + (n + 1 xy 2a 2 n)t + 2n (x + y)b n + 2(1 xy) 2[a 2 n+1 + a2 n + b 2 n], n 0 n 1 + b i + (x + y)(2a 2 n n 1 2 ) 2a2 n(b n + b n 1 ), n 1 i=0
21 Recurrence Relations: Discrete Garnier Equations Proposition Density matrix is given by (N + 1)! ρ N+1 (x; y) = C 2 e x2 /2 y 2 /2 N, N+1 and N+1 = a 2 2 N N, N 1 subject to 0 = 1, 1 = µ 0. Coupled recurrences for the polynomial coefficients are n b i + (n + 1 xy)b n b 3 n + (x + y)[b 2 n n 1] i=0 + a 2 n+1 [x + y 2b n b n+1 ] a 2 n[x + y 2b n b n 1 ] = 0 and a 2 n+1 [2n + 5 2xy + 2(x + y)(b n+1 + b n ) 2(a 2 n+2 + a2 n+1 + b2 n+1 + b2 n + b n+1 b n )] a 2 n[2n + 1 2xy + 2(x + y)(b n + b n 1 ) 2(a 2 n + a 2 n 1 + b2 n + b 2 n 1 + b nb n 1 )] + (b n x)(b n y) = 0
22 Initial Data and Moments The initial data for b n and a 2 n are given by where for x < y b 1 = 0, b 0 = µ 1 µ 0, a 2 0 = 0, a2 1 = µ 0µ 2 µ 2 1 µ 2 0 µ 0 = π( xy)[1 + erf (x) erf (y)] + ye x2 xe y2 µ 1 = 1 2 π(x + y)[1 + erf (x) erf (y)] e x 2 + e y2 µ 2 = π( xy)[1 + erf (x) erf (y)] + (y 1 2 x)e x2 (x 1 2 y)e y2
23 Classical Solutions to L(1, 1, 3; 2) System Deformed Hermite Orthogonal Polynomials with weight w(t; x, y) = (t x) µ (t y) ν e t2, t Γ where Γ {, x, y, }. Numerator and denominator polynomials W(t) = (t x)(t y), 2V(t) = 2t(t x)(t y) + ν(t x) + µ(t y) Coefficient functions are and Θ n (t) = 2(t x)(t y) 2b n (t x y) + θ n, n 0 Ω n (t) = t(t x)(t y) (2a 2 n n)(t x y) ν(t x) µ(t y) + ω n, n 1
24 Discrete Garnier Equations Proposition Coupled first order recurrences n n + 1 for system {a 2 n, b n, θ n, ω n } θ n + 2b n x = [ω n + (2a 2 n n)x][ω n + (2a 2 n n)x ν(x y)] a 2 n(θ n 1 + 2b n 1 x) θ n + 2b n y = [ω n + (2a 2 n n)y][ω n + (2a 2 n n)y + µ(x y)] a 2 n(θ n 1 + 2b n 1 y) a 2 n+1 + a2 n = b 2 n (2n µ + ν θ n) ω n+1 + ω n = (x + y)θ n + b n (2xy θ n ) µx νy
25 2 2 matrix system Spectral Derivatives Y n (t; x, y) = ( pn ) ɛ n /w p n 1 ɛ n 1 /w where the associated functions are ɛ n (s; x, y) = dt w(t) p n(t) Γ s t, Proposition s Γ Assume Reµ > 0, Reν > 0 and x y. The spectral derivatives satisfy { t Y n = 2t ( ) ( ) a 0 1 n + A x 1 0 t x + A } y Y n AY n t y where the residue matrices are ω n + (2a 2 n n)y x y a n(θ n + 2b n y) x y A x (t; x, y; µ, ν) a n (θ n 1 + 2b n 1 y) x y ω n + (2a 2 n n)y x y µ, (1) and A y = A x x y,µ ν. Properties: TrA x = µ, det A x = 0
26 Two-variable Garnier Systems L(1, 1, 1, 1, 1; 2) κ 0, κ 1, κ, θ 1, θ 2 L(1, 1, 1, 2; 2) κ 0, κ 1, κ, θ 2 L(1, 1, 3; 2) κ 0, κ 1, κ L(1, 2, 2; 2) κ 0, κ 1, κ L(1, 4; 2) κ 0, κ L(2, 3; 2) κ 0, κ L(5; 2) α L( 9 2 ; 2)
27 Deformation Derivatives Proposition Assume Reµ > 0, Reν > 0 and x y. The deformation derivatives with respect to x or y satisfy { 1 x Y 1 n = 2 x y { 1 y Y 1 n = 2 x y ( ) θn 2b n y 0 0 θ n 1 + 2b n 1 y ( ) θn + 2b n x 0 0 θ n 1 2b n 1 x A } x Y n B x Y n t x A } y Y n B y Y n t y
28 N Asymptotics Special diagonal case: x = y = 0 yields b n = 0 and with a 2 n+1 [ 2n + 5 2a 2 n+2 [ ] n+1] 2a2 = a 2 n 2n + 1 2a 2 n 2a 2 n 1 a 2 0 = 0, a2 1 = 3 2 Solution?
29 N Asymptotics Special diagonal case: x = y = 0 yields b n = 0 and with a 2 n+1 [ 2n + 5 2a 2 n+2 [ ] n+1] 2a2 = a 2 n 2n + 1 2a 2 n 2a 2 n 1 a 2 0 = 0, a2 1 = 3 2 Solution? a 2 n = n ( 1)n
30 N Asymptotics Special diagonal case: x = y = 0 yields b n = 0 and with a 2 n+1 [ 2n + 5 2a 2 n+2 [ ] n+1] 2a2 = a 2 n 2n + 1 2a 2 n 2a 2 n 1 a 2 0 = 0, a2 1 = 3 2 Solution? or for general y = x 0 a 2 n = n ( 1)n
31 N Asymptotics Special diagonal case: x = y = 0 yields b n = 0 and with a 2 n+1 [ 2n + 5 2a 2 n+2 [ ] n+1] 2a2 = a 2 n 2n + 1 2a 2 n 2a 2 n 1 Solution? or for general y = x 0 a 2 0 = 0, a2 1 = 3 2 a 2 n = n ( 1)n [ a 2 n = n (n + 2)H 2 n+1 (n + 1)H ] [ n+2h n n H 2 n 1 (n 1)H ] nh n 2 2 [ ] 2 (n + 1)H 2 n n H n+1 H n 1
32 Where to from here? Can the density matrix eigenvalues (or just the lowest one) be accessed directly and analytically using these ideas? the theory is easily generalisable to calculate other quantities, higher correlations etc, take the thermodynamic limit, finite temperature extensions may be possible within the structures, the situation for finite interparticle repulsion c < is more challenging, certainly other non-integrable potentials can be handled with a little more effort with the resulting difference and differential systems being non-integrable extensions to 2 dimensions would require a major task at the level of the basic mathematical theory
33 Selected References M. D. Girardeau. Relationship between systems of impenetrable bosons and fermions in one dimension. J. Math. Phys., 1: , E. H. Lieb and W. Liniger. Exact analysis of an interacting Bose gas. I. The general solution and the ground state. Phys. Rev. (2), 130: , A. Lenard. Momentum distribution in the ground state of the one-dimensional system of impenetrable bosons. J. Math. Phys., 5(7): , H. G. Vaidya and C. A. Tracy. One particle reduced density matrix of impenetrable bosons in one dimension at zero temperature. J. Math. Phys., 20(11): , M. Jimbo, T. Miwa, Y. Môri, and M. Sato. Density matrix of an impenetrable Bose gas and the fifth Painlevé transcendent. Phys. D, 1(1):80-158, H. Kimura. The degeneration of the two-dimensional Garnier system and the polynomial Hamiltonian structure. Ann. Mat. Pura Appl. 155:(4), 25-74, A. P. Magnus. Painlevé-type differential equations for the recurrence coefficients of semi-classical orthogonal polynomials. J. Comput. Appl. Math., 57(1-2): , P. J. Forrester, N. E. Frankel, T. M. Garoni, and N. S. Witte. Painlevé transcendent evaluations of finite system density matrices for 1d impenetrable bosons. Comm. Math. Phys., 238(1-2): , H. Kawamuko. On the Garnier system of half-integer type in two variables. Funkcial. Ekvac., 52(2): , 2009.
Triangular matrices and biorthogonal ensembles
/26 Triangular matrices and biorthogonal ensembles Dimitris Cheliotis Department of Mathematics University of Athens UK Easter Probability Meeting April 8, 206 2/26 Special densities on R n Example. n
More informationNear extreme eigenvalues and the first gap of Hermitian random matrices
Near extreme eigenvalues and the first gap of Hermitian random matrices Grégory Schehr LPTMS, CNRS-Université Paris-Sud XI Sydney Random Matrix Theory Workshop 13-16 January 2014 Anthony Perret, G. S.,
More informationMultiple orthogonal polynomials. Bessel weights
for modified Bessel weights KU Leuven, Belgium Madison WI, December 7, 2013 Classical orthogonal polynomials The (very) classical orthogonal polynomials are those of Jacobi, Laguerre and Hermite. Classical
More informationPainlevé Transcendents and the Hankel Determinants Generated by a Discontinuous Gaussian Weight
Painlevé Transcendents and the Hankel Determinants Generated by a Discontinuous Gaussian Weight arxiv:1803.10085v [math-ph] May 018 Chao Min and Yang Chen November 9, 018 Abstract This paper studies the
More informationOrthogonal 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 informationTemperature Correlation Functions in the XXO Heisenberg Chain
CONGRESSO NAZIONALE DI FISICA DELLA MATERIA Brescia, 13-16 June, 1994 Temperature Correlation Functions in the XXO Heisenberg Chain F. Colomo 1, A.G. Izergin 2,3, V.E. Korepin 4, V. Tognetti 1,5 1 I.N.F.N.,
More informationGeometric Dyson Brownian motion and May Wigner stability
Geometric Dyson Brownian motion and May Wigner stability Jesper R. Ipsen University of Melbourne Summer school: Randomness in Physics and Mathematics 2 Bielefeld 2016 joint work with Henning Schomerus
More informationarxiv:cond-mat/ v1 [cond-mat.soft] 3 Feb 2000
Many-body solitons in a one-dimensional condensate of hard core bosons M. D. Girardeau 1,2 and E. M. Wright 2 1 Institute of Theoretical Science, University of Oregon, Eugene, OR 97403 2 Optical Sciences
More informationThis ODE arises in many physical systems that we shall investigate. + ( + 1)u = 0. (λ + s)x λ + s + ( + 1) a λ. (s + 1)(s + 2) a 0
Legendre equation This ODE arises in many physical systems that we shall investigate We choose We then have Substitution gives ( x 2 ) d 2 u du 2x 2 dx dx + ( + )u u x s a λ x λ a du dx λ a λ (λ + s)x
More informationOrthogonal Polynomials, Perturbed Hankel Determinants. and. Random Matrix Models
Orthogonal Polynomials, Perturbed Hankel Determinants and Random Matrix Models A thesis presented for the degree of Doctor of Philosophy of Imperial College London and the Diploma of Imperial College by
More informationLecture 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 informationPainlevé VI and Hankel determinants for the generalized Jacobi weight
arxiv:98.558v2 [math.ca] 3 Nov 29 Painlevé VI and Hankel determinants for the generalized Jacobi weight D. Dai Department of Mathematics, City University of Hong Kong, Tat Chee Avenue, Kowloon, Hong Kong
More informationLecture 4.6: Some special orthogonal functions
Lecture 4.6: Some special orthogonal functions Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 4340, Advanced Engineering Mathematics
More informationBEC in one dimension
BEC in one dimension Tilmann John 11. Juni 2013 Outline 1 one-dimensional BEC 2 theoretical description Tonks-Girardeau gas Interaction exact solution (Lieb and Liniger) 3 experimental realization 4 conclusion
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 informationNONINTEGER FLUXES, DOLBEAULT COMPLEXES, AND SUPERSYMMETRIC QUANTUM MECHANICS. based on [ ] and [ ] Hannover, August 1, 2011
NONINTEGER FLUXES, DOLBEAULT COMPLEXES, AND SUPERSYMMETRIC QUANTUM MECHANICS based on [1104.3986] and [1105.3935] Hannover, August 1, 2011 WHAT IS WRONG WITH NONINTEGER FLUX? Quantization of Dirac monopole
More informationRANDOM MATRIX THEORY AND TOEPLITZ DETERMINANTS
RANDOM MATRIX THEORY AND TOEPLITZ DETERMINANTS David García-García May 13, 2016 Faculdade de Ciências da Universidade de Lisboa OVERVIEW Random Matrix Theory Introduction Matrix ensembles A sample computation:
More informationFinal Exam May 4, 2016
1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.
More informationMarkov operators, classical orthogonal polynomial ensembles, and random matrices
Markov operators, classical orthogonal polynomial ensembles, and random matrices M. Ledoux, Institut de Mathématiques de Toulouse, France 5ecm Amsterdam, July 2008 recent study of random matrix and random
More informationQuantum Mechanics Solutions. λ i λ j v j v j v i v i.
Quantum Mechanics Solutions 1. (a) If H has an orthonormal basis consisting of the eigenvectors { v i } of A with eigenvalues λ i C, then A can be written in terms of its spectral decomposition as A =
More informationDeterminantal point processes and random matrix theory in a nutshell
Determinantal point processes and random matrix theory in a nutshell part II Manuela Girotti based on M. Girotti s PhD thesis, A. Kuijlaars notes from Les Houches Winter School 202 and B. Eynard s notes
More informationAlmost Product Evaluation of Hankel Determinants
Almost Product Evaluation of Hankel Determinants Ömer Eğecioğlu Department of Computer Science University of California, Santa Barbara CA 93106 omer@csucsbedu Timothy Redmond Stanford Medical Informatics,
More informationarxiv:hep-th/ v1 14 Oct 1992
ITD 92/93 11 Level-Spacing Distributions and the Airy Kernel Craig A. Tracy Department of Mathematics and Institute of Theoretical Dynamics, University of California, Davis, CA 95616, USA arxiv:hep-th/9210074v1
More informationLEGENDRE POLYNOMIALS AND APPLICATIONS. We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.
LEGENDRE POLYNOMIALS AND APPLICATIONS We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.. Legendre equation: series solutions The Legendre equation is
More informationLecture 6. Four postulates of quantum mechanics. The eigenvalue equation. Momentum and energy operators. Dirac delta function. Expectation values
Lecture 6 Four postulates of quantum mechanics The eigenvalue equation Momentum and energy operators Dirac delta function Expectation values Objectives Learn about eigenvalue equations and operators. Learn
More informationPainlevé IV and degenerate Gaussian Unitary Ensembles
Painlevé IV and degenerate Gaussian Unitary Ensembles arxiv:math-ph/0606064v1 8 Jun 006 Yang Chen Department of Mathematics Imperial College London 180 Queen s Gates London SW7 BZ UK M. V. Feigin Department
More informationTwo special equations: Bessel s and Legendre s equations. p Fourier-Bessel and Fourier-Legendre series. p
LECTURE 1 Table of Contents Two special equations: Bessel s and Legendre s equations. p. 259-268. Fourier-Bessel and Fourier-Legendre series. p. 453-460. Boundary value problems in other coordinate system.
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 informationSimple Harmonic Oscillator
Classical harmonic oscillator Linear force acting on a particle (Hooke s law): F =!kx From Newton s law: F = ma = m d x dt =!kx " d x dt + # x = 0, # = k / m Position and momentum solutions oscillate in
More informationMaximal height of non-intersecting Brownian motions
Maximal height of non-intersecting Brownian motions G. Schehr Laboratoire de Physique Théorique et Modèles Statistiques CNRS-Université Paris Sud-XI, Orsay Collaborators: A. Comtet (LPTMS, Orsay) P. J.
More informationPage 404. Lecture 22: Simple Harmonic Oscillator: Energy Basis Date Given: 2008/11/19 Date Revised: 2008/11/19
Page 404 Lecture : Simple Harmonic Oscillator: Energy Basis Date Given: 008/11/19 Date Revised: 008/11/19 Coordinate Basis Section 6. The One-Dimensional Simple Harmonic Oscillator: Coordinate Basis Page
More informationPainlevé equations and orthogonal polynomials
KU Leuven, Belgium Kapaev workshop, Ann Arbor MI, 28 August 2017 Contents Painlevé equations (discrete and continuous) appear at various places in the theory of orthogonal polynomials: Discrete Painlevé
More informationProblem 1: Spin 1 2. particles (10 points)
Problem 1: Spin 1 particles 1 points 1 Consider a system made up of spin 1/ particles. If one measures the spin of the particles, one can only measure spin up or spin down. The general spin state of a
More informationDifferential Equations for Dyson Processes
Differential Equations for Dyson Processes Joint work with Harold Widom I. Overview We call Dyson process any invariant process on ensembles of matrices in which the entries undergo diffusion. Dyson Brownian
More informationDeterminantal point processes and random matrix theory in a nutshell
Determinantal point processes and random matrix theory in a nutshell part I Manuela Girotti based on M. Girotti s PhD thesis and A. Kuijlaars notes from Les Houches Winter School 202 Contents Point Processes
More informationA Probability Review
A Probability Review Outline: A probability review Shorthand notation: RV stands for random variable EE 527, Detection and Estimation Theory, # 0b 1 A Probability Review Reading: Go over handouts 2 5 in
More informationBulk scaling limits, open questions
Bulk scaling limits, open questions Based on: Continuum limits of random matrices and the Brownian carousel B. Valkó, B. Virág. Inventiones (2009). Eigenvalue statistics for CMV matrices: from Poisson
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 informationStochastic nonlinear Schrödinger equations and modulation of solitary waves
Stochastic nonlinear Schrödinger equations and modulation of solitary waves A. de Bouard CMAP, Ecole Polytechnique, France joint work with R. Fukuizumi (Sendai, Japan) Deterministic and stochastic front
More informationLecture 5: Harmonic oscillator, Morse Oscillator, 1D Rigid Rotor
Lecture 5: Harmonic oscillator, Morse Oscillator, 1D Rigid Rotor It turns out that the boundary condition of the wavefunction going to zero at infinity is sufficient to quantize the value of energy that
More informationTESTING ADS/CFT. John H. Schwarz STRINGS 2003
TESTING ADS/CFT John H. Schwarz STRINGS 2003 July 6, 2003 1 INTRODUCTION During the past few years 1 Blau et al. constructed a maximally supersymmetric plane-wave background of type IIB string theory as
More informationThe Quantum Heisenberg Ferromagnet
The Quantum Heisenberg Ferromagnet Soon after Schrödinger discovered the wave equation of quantum mechanics, Heisenberg and Dirac developed the first successful quantum theory of ferromagnetism W. Heisenberg,
More informationExtreme eigenvalue fluctutations for GUE
Extreme eigenvalue fluctutations for GUE C. Donati-Martin 204 Program Women and Mathematics, IAS Introduction andom matrices were introduced in multivariate statistics, in the thirties by Wishart [Wis]
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 informationIntroduction to Orthogonal Polynomials: Definition and basic properties
Introduction to Orthogonal Polynomials: Definition and basic properties Prof. Dr. Mama Foupouagnigni African Institute for Mathematical Sciences, Limbe, Cameroon and Department of Mathematics, Higher Teachers
More informationLecture 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 informationMath 266, Midterm Exam 1
Math 266, Midterm Exam 1 February 19th 2016 Name: Ground Rules: 1. Calculator is NOT allowed. 2. Show your work for every problem unless otherwise stated (partial credits are available). 3. You may use
More informationSecond Quantization Method for Bosons
Second Quantization Method for Bosons Hartree-Fock-based methods cannot describe the effects of the classical image potential (cf. fig. 1) because HF is a mean-field theory. DFF-LDA is not able either
More informationSection 10: Many Particle Quantum Mechanics Solutions
Physics 143a: Quantum Mechanics I Section 10: Many Particle Quantum Mechanics Solutions Spring 015, Harvard Here is a summary of the most important points from this week (with a few of my own tidbits),
More informationExponential tail inequalities for eigenvalues of random matrices
Exponential tail inequalities for eigenvalues of random matrices M. Ledoux Institut de Mathématiques de Toulouse, France exponential tail inequalities classical theme in probability and statistics quantify
More informationThe Multivariate Normal Distribution. In this case according to our theorem
The Multivariate Normal Distribution Defn: Z R 1 N(0, 1) iff f Z (z) = 1 2π e z2 /2. Defn: Z R p MV N p (0, I) if and only if Z = (Z 1,..., Z p ) T with the Z i independent and each Z i N(0, 1). In this
More informationPhysics 250 Green s functions for ordinary differential equations
Physics 25 Green s functions for ordinary differential equations Peter Young November 25, 27 Homogeneous Equations We have already discussed second order linear homogeneous differential equations, which
More informationMultifractality in simple systems. Eugene Bogomolny
Multifractality in simple systems Eugene Bogomolny Univ. Paris-Sud, Laboratoire de Physique Théorique et Modèles Statistiques, Orsay, France In collaboration with Yasar Atas Outlook 1 Multifractality Critical
More informationKet space as a vector space over the complex numbers
Ket space as a vector space over the complex numbers kets ϕ> and complex numbers α with two operations Addition of two kets ϕ 1 >+ ϕ 2 > is also a ket ϕ 3 > Multiplication with complex numbers α ϕ 1 >
More informationRandom matrix pencils and level crossings
Albeverio Fest October 1, 2018 Topics to discuss Basic level crossing problem 1 Basic level crossing problem 2 3 Main references Basic level crossing problem (i) B. Shapiro, M. Tater, On spectral asymptotics
More informationORTHOGONAL POLYNOMIALS
ORTHOGONAL POLYNOMIALS 1. PRELUDE: THE VAN DER MONDE DETERMINANT The link between random matrix theory and the classical theory of orthogonal polynomials is van der Monde s determinant: 1 1 1 (1) n :=
More informationBefore you begin read these instructions carefully.
MATHEMATICAL TRIPOS Part IB Tuesday, 31 May, 2016 9:00 am to 12:00 pm PAPER 1 Before you begin read these instructions carefully. Each question in Section II carries twice the number of marks of each question
More informationOn level crossing in deterministic and random matrix pencils
On level crossing in deterministic and random matrix pencils May 3, 2018 Topics to discuss 1 Basic level crossing problem 2 3 4 Main references (i) B. Shapiro, M. Tater, On spectral asymptotics of quasi-exactly
More informationarxiv:cond-mat/ v1 [cond-mat.str-el] 15 Jul 2005
Correlation functions of one-dimensional Bose-Fermi mixtures Holger Frahm and Guillaume Palacios Institut für Theoretische Physik, Universität Hannover, Appelstr. 2, 30167 Hannover, Germany (Dated: July
More informationOn the concentration of eigenvalues of random symmetric matrices
On the concentration of eigenvalues of random symmetric matrices Noga Alon Michael Krivelevich Van H. Vu April 23, 2012 Abstract It is shown that for every 1 s n, the probability that the s-th largest
More informationTwo-variable Wilson polynomials and the generic superintegrable system on the 3-sphere
Two-variable Wilson polynomials and the generic superintegrable system on the 3-sphere Willard Miller, [Joint with E.G. Kalnins (Waikato) and Sarah Post (CRM)] University of Minnesota Special Functions
More informationChapter 29. Quantum Chaos
Chapter 29 Quantum Chaos What happens to a Hamiltonian system that for classical mechanics is chaotic when we include a nonzero h? There is no problem in principle to answering this question: given a classical
More informationNumerical Methods for Random Matrices
Numerical Methods for Random Matrices MIT 18.95 IAP Lecture Series Per-Olof Persson (persson@mit.edu) January 23, 26 1.9.8.7 Random Matrix Eigenvalue Distribution.7.6.5 β=1 β=2 β=4 Probability.6.5.4.4
More informationQuantum Mechanics: Vibration and Rotation of Molecules
Quantum Mechanics: Vibration and Rotation of Molecules 8th April 2008 I. 1-Dimensional Classical Harmonic Oscillator The classical picture for motion under a harmonic potential (mass attached to spring
More informationPHY 407 QUANTUM MECHANICS Fall 05 Problem set 1 Due Sep
Problem set 1 Due Sep 15 2005 1. Let V be the set of all complex valued functions of a real variable θ, that are periodic with period 2π. That is u(θ + 2π) = u(θ), for all u V. (1) (i) Show that this V
More informationChemistry 432 Problem Set 4 Spring 2018 Solutions
Chemistry 4 Problem Set 4 Spring 18 Solutions 1. V I II III a b c A one-dimensional particle of mass m is confined to move under the influence of the potential x a V V (x) = a < x b b x c elsewhere and
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 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 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 informationMath 308 Final Exam Practice Problems
Math 308 Final Exam Practice Problems This review should not be used as your sole source for preparation for the exam You should also re-work all examples given in lecture and all suggested homework problems
More informationLecture #8: Quantum Mechanical Harmonic Oscillator
5.61 Fall, 013 Lecture #8 Page 1 Last time Lecture #8: Quantum Mechanical Harmonic Oscillator Classical Mechanical Harmonic Oscillator * V(x) = 1 kx (leading term in power series expansion of most V(x)
More informationDouble contour integral formulas for the sum of GUE and one matrix model
Double contour integral formulas for the sum of GUE and one matrix model Based on arxiv:1608.05870 with Tom Claeys, Arno Kuijlaars, and Karl Liechty Dong Wang National University of Singapore Workshop
More informationNumerical Evaluation of Standard Distributions in Random Matrix Theory
Numerical Evaluation of Standard Distributions in Random Matrix Theory A Review of Folkmar Bornemann s MATLAB Package and Paper Matt Redmond Department of Mathematics Massachusetts Institute of Technology
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 informationCOMPLEX HERMITE POLYNOMIALS: FROM THE SEMI-CIRCULAR LAW TO THE CIRCULAR LAW
Serials Publications www.serialspublications.com OMPLEX HERMITE POLYOMIALS: FROM THE SEMI-IRULAR LAW TO THE IRULAR LAW MIHEL LEDOUX Abstract. We study asymptotics of orthogonal polynomial measures of the
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 informationEach problem is worth 34 points. 1. Harmonic Oscillator Consider the Hamiltonian for a simple harmonic oscillator. 2ml 2 0. d 2
Physics 443 Prelim # with solutions March 7, 8 Each problem is worth 34 points.. Harmonic Oscillator Consider the Hamiltonian for a simple harmonic oscillator H p m + mω x (a Use dimensional analysis to
More informationChemistry 532 Problem Set 7 Spring 2012 Solutions
Chemistry 53 Problem Set 7 Spring 01 Solutions 1. The study of the time-independent Schrödinger equation for a one-dimensional particle subject to the potential function leads to the differential equation
More informationBessel Functions Michael Taylor. Lecture Notes for Math 524
Bessel Functions Michael Taylor Lecture Notes for Math 54 Contents 1. Introduction. Conversion to first order systems 3. The Bessel functions J ν 4. The Bessel functions Y ν 5. Relations between J ν and
More informationMATH3383. Quantum Mechanics. Appendix D: Hermite Equation; Orthogonal Polynomials
MATH3383. Quantum Mechanics. Appendix D: Hermite Equation; Orthogonal Polynomials. Hermite Equation In the study of the eigenvalue problem of the Hamiltonian for the quantum harmonic oscillator we have
More informationIntroduction to Integrability
Introduction to Integrability Problem Sets ETH Zurich, HS16 Prof. N. Beisert, A. Garus c 2016 Niklas Beisert, ETH Zurich This document as well as its parts is protected by copyright. This work is licensed
More informationEEL 5544 Noise in Linear Systems Lecture 30. X (s) = E [ e sx] f X (x)e sx dx. Moments can be found from the Laplace transform as
L30-1 EEL 5544 Noise in Linear Systems Lecture 30 OTHER TRANSFORMS For a continuous, nonnegative RV X, the Laplace transform of X is X (s) = E [ e sx] = 0 f X (x)e sx dx. For a nonnegative RV, the Laplace
More informationAtkins & de Paula: Atkins Physical Chemistry 9e Checklist of key ideas. Chapter 8: Quantum Theory: Techniques and Applications
Atkins & de Paula: Atkins Physical Chemistry 9e Checklist of key ideas Chapter 8: Quantum Theory: Techniques and Applications TRANSLATIONAL MOTION wavefunction of free particle, ψ k = Ae ikx + Be ikx,
More informationModels with fundamental length 1 and the finite-dimensional simulations of Big Bang
Models with fundamental length 1 and the finite-dimensional simulations of Big Bang Miloslav Znojil Nuclear Physics Institute ASCR, 250 68 Řež, Czech Republic 1 talk in Dresden (June 22, 2011, 10.50-11.35
More informationBrief review of Quantum Mechanics (QM)
Brief review of Quantum Mechanics (QM) Note: This is a collection of several formulae and facts that we will use throughout the course. It is by no means a complete discussion of QM, nor will I attempt
More informationDerivation of quantum hydrodynamic equations with Fermi-Dirac and Bose-Einstein statistics
Derivation of quantum hydrodynamic equations with Fermi-Dirac and Bose-Einstein statistics Luigi Barletti (Università di Firenze) Carlo Cintolesi (Università di Trieste) 6th MMKT Porto Ercole, june 9th
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 informationThe nature of superfluidity in the cold atomic unitary Fermi gas
The nature of superfluidity in the cold atomic unitary Fermi gas Introduction Yoram Alhassid (Yale University) Finite-temperature auxiliary-field Monte Carlo (AFMC) method The trapped unitary Fermi gas
More informationPhysics 215 Quantum Mechanics 1 Assignment 5
Physics 15 Quantum Mechanics 1 Assignment 5 Logan A. Morrison February 10, 016 Problem 1 A particle of mass m is confined to a one-dimensional region 0 x a. At t 0 its normalized wave function is 8 πx
More informationMathematical Tripos Part IB Michaelmas Term Example Sheet 1. Values of some physical constants are given on the supplementary sheet
Mathematical Tripos Part IB Michaelmas Term 2015 Quantum Mechanics Dr. J.M. Evans Example Sheet 1 Values of some physical constants are given on the supplementary sheet 1. Whenasampleofpotassiumisilluminatedwithlightofwavelength3
More informationFourier series on triangles and tetrahedra
Fourier series on triangles and tetrahedra Daan Huybrechs K.U.Leuven Cambridge, Sep 13, 2010 Outline Introduction The unit interval Triangles, tetrahedra and beyond Outline Introduction The topic of this
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 informationSusceptibility of the Two-Dimensional Ising Model
Susceptibility of the Two-Dimensional Ising Model Craig A. Tracy UC Davis September 2014 1. Definition of 2D Ising Model Outline Outline 1. Definition of 2D Ising Model 2. Why is the nearest neighbor zero-field
More informationA Detailed Look at a Discrete Randomw Walk with Spatially Dependent Moments and Its Continuum Limit
A Detailed Look at a Discrete Randomw Walk with Spatially Dependent Moments and Its Continuum Limit David Vener Department of Mathematics, MIT May 5, 3 Introduction In 8.366, we discussed the relationship
More informationSymmetry properties, density profiles and momentum distribution of multicomponent mixtures of strongly interacting 1D Fermi gases
Symmetry properties, density profiles and momentum distribution of multicomponent mixtures of strongly interacting 1D Fermi gases Anna Minguzzi LPMMC Université Grenoble Alpes and CNRS Multicomponent 1D
More informationBose-Einstein condensates under rotation: The structures within
Bose-Einstein condensates under rotation: The structures within Peter Mason Centre de Mathématiques Appliquées, Ecole Polytechnique Paris, France In collaboration with Amandine Aftalion and Thierry Jolicoeur:
More informationCentral 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 informationDegree Master of Science in Mathematical Modelling and Scientific Computing Mathematical Methods I Thursday, 12th January 2012, 9:30 a.m.- 11:30 a.m.
Degree Master of Science in Mathematical Modelling and Scientific Computing Mathematical Methods I Thursday, 12th January 2012, 9:30 a.m.- 11:30 a.m. Candidates should submit answers to a maximum of four
More informationFredholm 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 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 information