Haydock s recursive solution of self-adjoint problems. Discrete spectrum
|
|
- Curtis Higgins
- 5 years ago
- Views:
Transcription
1 Haydock s recursive solution of self-adjoint problems. Discrete spectrum Alexander Moroz Wave-scattering.com wavescattering@yahoo.com January 3, 2015 Alexander Moroz (WS) Recursive solution January 3, / 37
2 Overview Haydock s recursive solution Alexander Moroz (WS) Recursive solution January 3, / 37
3 Overview Haydock s recursive solution Novel numerical method Alexander Moroz (WS) Recursive solution January 3, / 37
4 Overview Haydock s recursive solution Novel numerical method Rabi model Alexander Moroz (WS) Recursive solution January 3, / 37
5 Overview Haydock s recursive solution Novel numerical method Rabi model Exactly solvable models Alexander Moroz (WS) Recursive solution January 3, / 37
6 Overview Haydock s recursive solution Novel numerical method Rabi model Exactly solvable models Quasi-exactly solvable (QES) models Alexander Moroz (WS) Recursive solution January 3, / 37
7 Overview Haydock s recursive solution Novel numerical method Rabi model Exactly solvable models Quasi-exactly solvable (QES) models Conclusions Alexander Moroz (WS) Recursive solution January 3, / 37
8 The message I m attempting to transmit Tridiagonality is fundamental Alexander Moroz (WS) Recursive solution January 3, / 37
9 The message I m attempting to transmit Tridiagonality is fundamental Three-term recurrence relations are indispensable Alexander Moroz (WS) Recursive solution January 3, / 37
10 The message I m attempting to transmit Tridiagonality is fundamental Three-term recurrence relations are indispensable Orthogonal polynomial systems are unavoidable Alexander Moroz (WS) Recursive solution January 3, / 37
11 The message I m attempting to transmit Tridiagonality is fundamental Three-term recurrence relations are indispensable Orthogonal polynomial systems are unavoidable an efficient numerical method for energy levels Alexander Moroz (WS) Recursive solution January 3, / 37
12 The message I m attempting to transmit Tridiagonality is fundamental Three-term recurrence relations are indispensable Orthogonal polynomial systems are unavoidable an efficient numerical method for energy levels characterization of integrable and QES models in terms of intrinsic polynomial properties Alexander Moroz (WS) Recursive solution January 3, / 37
13 The message I m attempting to transmit Tridiagonality is fundamental Three-term recurrence relations are indispensable Orthogonal polynomial systems are unavoidable an efficient numerical method for energy levels characterization of integrable and QES models in terms of intrinsic polynomial properties confine integrable spectra to four basic classes Alexander Moroz (WS) Recursive solution January 3, / 37
14 The message I m attempting to transmit Tridiagonality is fundamental Three-term recurrence relations are indispensable Orthogonal polynomial systems are unavoidable an efficient numerical method for energy levels characterization of integrable and QES models in terms of intrinsic polynomial properties confine integrable spectra to four basic classes characterization of chaos? Alexander Moroz (WS) Recursive solution January 3, / 37
15 Trivial and always valid statement For any linear self-adjoint Ĥ there always exists an orthonormal basis {e n } n=0 such that eigenstates E of Ĥ Ĥ E = E E can be expanded as E = p n (E)e n n=0 Alexander Moroz (WS) Recursive solution January 3, / 37
16 Haydock s recursive solution (H1) the expansion coefficients {p n (E)} are polynomials with degree p n = n, orthonormal with respect to the density of states (DOS), n 0 (E), p n (E)p m (E)n 0 (E)dE = δ nm Alexander Moroz (WS) Recursive solution January 3, / 37
17 Haydock s recursive solution (H1) the expansion coefficients {p n (E)} are polynomials with degree p n = n, orthonormal with respect to the density of states (DOS), n 0 (E), p n (E)p m (E)n 0 (E)dE = δ nm (H2) energy eigenstates E are the generating function of the orthogonal polynomials Alexander Moroz (WS) Recursive solution January 3, / 37
18 Haydock s recursive solution (H1) the expansion coefficients {p n (E)} are polynomials with degree p n = n, orthonormal with respect to the density of states (DOS), n 0 (E), p n (E)p m (E)n 0 (E)dE = δ nm (H2) energy eigenstates E are the generating function of the orthogonal polynomials (H3) the orthogonality of the energy eigenstates, E E = δ EE yields a dual orthogonality relation n 0 (E) p n (E)p n (E ) = δ EE n=0 where E and E are both eigenvalues Alexander Moroz (WS) Recursive solution January 3, / 37
19 Tridiagonality is generic There always exists an orthonormal basis {e n } n=0 self-adjoint operator takes on a tridiagonal form such that a given He n = a n e n + b n+1 e n+1 + b n e n 1 (1) with real recurrence coefficients {a n } and {b n }, where b n 0, n 0. Alexander Moroz (WS) Recursive solution January 3, / 37
20 Proof of He n = a n e n + b n+1 e n+1 + b n e n 1 n = 0: find a 0, b 1, e 1 He 0 = a 0 e 0 + b 1 e 1 a 0 = e 0, He 0 b 2 1 = (H a 0)e 0, (H a 0 )e 0 e 1 = (H a 0 )e 0 /b 1 Alexander Moroz (WS) Recursive solution January 3, / 37
21 Proof of He n = a n e n + b n+1 e n+1 + b n e n 1 n = 0: find a 0, b 1, e 1 He 0 = a 0 e 0 + b 1 e 1 a 0 = e 0, He 0 b 2 1 = (H a 0)e 0, (H a 0 )e 0 e 1 = (H a 0 )e 0 /b 1 Induction step to determine a n, b n+1, e n+1 e n 1, He n = e n, He n 1 b n a n = e n, He n b n+1 e n+1 = (H a n )e n b n e n 1 b 2 n+1 = (H a n)e n b n e n 1, (H a n )e n b n e n 1 e n+1 = [(H a n )e n b n e n 1 ]/b n+1 By construction e n+1 normalized to one and orthogonal to e n and e n 1. Alexander Moroz (WS) Recursive solution January 3, / 37
22 Final check: e n+1 orthogonal to e n 2,...,e 0 Tridiagonality of He n, self-adjointness, and normalizability b n+1 e m, e n+1 = e m, He n = He m, e n Alexander Moroz (WS) Recursive solution January 3, / 37
23 Final check: e n+1 orthogonal to e n 2,...,e 0 Tridiagonality of He n, self-adjointness, and normalizability b n+1 e m, e n+1 = e m, He n = He m, e n Tridiagonality for m < n 1 He m is a linear combination of e m 1, e m, and e m+1, all of which have zero overlap with e n if m + 1 < n. Alexander Moroz (WS) Recursive solution January 3, / 37
24 Relation to orthogonal polynomials I Wave function expansion E = p n (E)e n n=0 Alexander Moroz (WS) Recursive solution January 3, / 37
25 Relation to orthogonal polynomials I Wave function expansion E = p n (E)e n n=0 Three-term recurrence relation (TTRR) Expansion coefficients satisfy Ep n (E) = a n p n (E) + b n+1 p n+1 (E) + b n p n 1 (E) (2) with b k 0 and an initial condition p 0 = 1 and p 1 = 0. {p n (E)} are by the very definition orthogonal polynomials. Alexander Moroz (WS) Recursive solution January 3, / 37
26 Favard s theorem (Theorem I-4.4 of Chihara s book) For given arbitrary sequences of complex numbers {c n } and {λ n } in the TTRR xp n = P n+1 + c n P n + λ n P n 1 there always exists a moment functional, that is a linear functional L acting in the space of (complex) monic polynomials C[E], such that the polynomials P n defined by the TTRR are orthogonal under L: L(P k P l ) = 0 k l N The functional L is unique if we impose the normalization condition L(P 0 ) = L(1) = µ 0, where µ 0 is a chosen positive constant. Alexander Moroz (WS) Recursive solution January 3, / 37
27 Relation to orthogonal polynomials II Shorthand for basis e n = p n (H)e 0 Alexander Moroz (WS) Recursive solution January 3, / 37
28 Relation to orthogonal polynomials II Shorthand for basis e n = p n (H)e 0 Orthogonality relations e m, e n = p m (H)e 0, p n (H)e 0 = δ mn Alexander Moroz (WS) Recursive solution January 3, / 37
29 Relation to orthogonal polynomials III Change to the orthonormal basis of eigenstates e 0 = k ω k ψ(e k ) Alexander Moroz (WS) Recursive solution January 3, / 37
30 Relation to orthogonal polynomials III Change to the orthonormal basis of eigenstates e 0 = k ω k ψ(e k ) Orthogonal polynomials of a discrete variable e m, e n = k ω k 2 p m (E k )p n (E k ) = δ mn Alexander Moroz (WS) Recursive solution January 3, / 37
31 Relation to orthogonal polynomials III Change to the orthonormal basis of eigenstates e 0 = k ω k ψ(e k ) Orthogonal polynomials of a discrete variable e m, e n = k ω k 2 p m (E k )p n (E k ) = δ mn Weight function is the local density of function (DOS) n(e) = k ω k 2 δ(e E k ) Alexander Moroz (WS) Recursive solution January 3, / 37
32 Weight function more rigorously (pp of Chihara s book) lim n x nl = ξ l σ lim j ξ j Alexander Moroz (WS) Recursive solution January 3, / 37
33 Weight function more rigorously (pp of Chihara s book) lim n x nl = ξ l σ lim j ξ j (a) ξ l = σ = (l 1) Alexander Moroz (WS) Recursive solution January 3, / 37
34 Weight function more rigorously (pp of Chihara s book) lim n x nl = ξ l σ lim j ξ j (a) ξ l = σ = (l 1) (b) < ξ 1 < ξ 2 <... < ξ l = σ for some l 1 Alexander Moroz (WS) Recursive solution January 3, / 37
35 Weight function more rigorously (pp of Chihara s book) lim n x nl = ξ l σ lim j ξ j (a) ξ l = σ = (l 1) (b) < ξ 1 < ξ 2 <... < ξ l = σ for some l 1 (c) < ξ 1 < ξ 2 <... < ξ l <... < σ = Alexander Moroz (WS) Recursive solution January 3, / 37
36 Weight function is the set of limit points of flows of zeros Alexander Moroz (WS) Recursive solution January 3, / 37
37 Numerical recipe Choose N c N 0 and determine the first N 0 zeros x Ncl, l N 0, of P Nc (x). Usually a good starting point is to take N c N Because P Nc (x) has N c simple zeros, any omission of a zero could be easily identified. Alexander Moroz (WS) Recursive solution January 3, / 37
38 Numerical recipe Choose N c N 0 and determine the first N 0 zeros x Ncl, l N 0, of P Nc (x). Usually a good starting point is to take N c N Because P Nc (x) has N c simple zeros, any omission of a zero could be easily identified. Gradually increase the cut-off value of N c. The latter is what drives the incessant flows of polynomial zeros x Ncl, wherein each flow is characterized by the parameter l. Alexander Moroz (WS) Recursive solution January 3, / 37
39 Numerical recipe Choose N c N 0 and determine the first N 0 zeros x Ncl, l N 0, of P Nc (x). Usually a good starting point is to take N c N Because P Nc (x) has N c simple zeros, any omission of a zero could be easily identified. Gradually increase the cut-off value of N c. The latter is what drives the incessant flows of polynomial zeros x Ncl, wherein each flow is characterized by the parameter l. Monitor convergence of the respective flows induced by the very first n zeros of P Nc (x). Each flow is a monotonically decreasing sequence having necessary a fixed limit point. Terminate your calculations when the N 0 -th zero of P Nc (x) converged to ξ N0 within predetermined accuracy. Then as a rule all other flows x Ncl with l < N 0 have converged, too. Alexander Moroz (WS) Recursive solution January 3, / 37
40 Rabi model Hamiltonian ĤR Ĥ R = ω1â â + gσ 1 (â + â) + µσ 3 ω 0... TLS resonance frequency, ω... cavity mode frequency µ = ω 0 /2, [â, â ] = 1, g... a coupling constant 1... the unit matrix σ j... the Pauli matrices Alexander Moroz (WS) Recursive solution January 3, / 37
41 Rabi model Hamiltonian ĤR Ĥ R = ω1â â + gσ 1 (â + â) + µσ 3 ω 0... TLS resonance frequency, ω... cavity mode frequency µ = ω 0 /2, [â, â ] = 1, g... a coupling constant 1... the unit matrix σ j... the Pauli matrices In what follows: κ = g/ω, = µ/ω, = 1 Alexander Moroz (WS) Recursive solution January 3, / 37
42 Convergence toward the 1000th eigenvalue of the Rabi model in positive parity eigenspace n, ,1 1E-3 1E-5 1E-7 1E-9 1E-11 κ=0.2, =0.4 κ=1.0, =0.7 1E (n-1000)/ ɛ 999 = , , for (κ, ) = (0.2, 0.4), (1, 0.7), (5, 0.4), respectively Alexander Moroz (WS) Recursive solution January 3, / 37
43 An approximation of the spectrum ɛ l = l κ 2 of displaced harmonic oscillator by the zeros x nl nl - l 10 0,1 1E-3 1E-5 1E ,1 1E-3 1E-5 1E ,1 1E-3 1E-5 1E l κ=2 n=100 n=500 n=1000 Alexander Moroz (WS) Recursive solution January 3, / 37
44 I m looking for somebody to analyze level statistics. My input: Thousands of energy levels per parity subspace. Alexander Moroz (WS) Recursive solution January 3, / 37
45 Comparison with Braak s PRL 107, (2011) A transcendental function is constructed [ G ± (ζ) = K n (ζ, κ) 1 ] κ n ζ n n=0 and one determines eigenvalues as zeros of G ± (ζ). Alexander Moroz (WS) Recursive solution January 3, / 37
46 Comparison with Braak s PRL 107, (2011) A transcendental function is constructed [ G ± (ζ) = K n (ζ, κ) 1 ] κ n ζ n n=0 and one determines eigenvalues as zeros of G ± (ζ). The expansion coefficients K n are determined by a TTRR φ n+1 f n(ζ) (n + 1) φ n + 1 n + 1 φ n 1 = 0 where f n (ζ) = 2κ + 1 ) (n ζ 2. 2κ n ζ Alexander Moroz (WS) Recursive solution January 3, / 37
47 Shortcomings of Braak s approach Relies heavily on unproven hypotheses that G ± (ζ) takes on a zero value between its subsequent poles at ζ = n and ζ = n + 1 once - by implicitly presuming that at one of the poles G ± (ζ) goes to + and at the neighboring pole goes to, with a monotonic behavior from + to between the poles; Alexander Moroz (WS) Recursive solution January 3, / 37
48 Shortcomings of Braak s approach Relies heavily on unproven hypotheses that G ± (ζ) takes on a zero value between its subsequent poles at ζ = n and ζ = n + 1 once - by implicitly presuming that at one of the poles G ± (ζ) goes to + and at the neighboring pole goes to, with a monotonic behavior from + to between the poles; twice - implicitly presuming that G ± (ζ) goes to one of ± at both subsequent poles, and in between the poles it has rather a featureless behavior, e.g., similar to a cord hanging on two posts; Alexander Moroz (WS) Recursive solution January 3, / 37
49 Shortcomings of Braak s approach Relies heavily on unproven hypotheses that G ± (ζ) takes on a zero value between its subsequent poles at ζ = n and ζ = n + 1 once - by implicitly presuming that at one of the poles G ± (ζ) goes to + and at the neighboring pole goes to, with a monotonic behavior from + to between the poles; twice - implicitly presuming that G ± (ζ) goes to one of ± at both subsequent poles, and in between the poles it has rather a featureless behavior, e.g., similar to a cord hanging on two posts; none - occurs under the similar circumstances as described in the previous item, if the cord is too short, e.g., it does not stretch sufficiently up or down so as to cross the abscissa. Alexander Moroz (WS) Recursive solution January 3, / 37
50 Shortcomings of Braak s approach Relies heavily on unproven hypotheses that G ± (ζ) takes on a zero value between its subsequent poles at ζ = n and ζ = n + 1 once - by implicitly presuming that at one of the poles G ± (ζ) goes to + and at the neighboring pole goes to, with a monotonic behavior from + to between the poles; twice - implicitly presuming that G ± (ζ) goes to one of ± at both subsequent poles, and in between the poles it has rather a featureless behavior, e.g., similar to a cord hanging on two posts; none - occurs under the similar circumstances as described in the previous item, if the cord is too short, e.g., it does not stretch sufficiently up or down so as to cross the abscissa. Numerically no advantage over the classical Schweber s solution [Ann. Phys. (N.Y.) 41, 205 (1967)]. Alexander Moroz (WS) Recursive solution January 3, / 37
51 Comparison with Braak s PRL 107, (2011) - II The TTRR for the coefficients K n is not a fundamental one, i.e. not of Haydock s type, because of f n (ζ) is not a linear function of energy. Alexander Moroz (WS) Recursive solution January 3, / 37
52 Comparison with Braak s PRL 107, (2011) - II The TTRR for the coefficients K n is not a fundamental one, i.e. not of Haydock s type, because of f n (ζ) is not a linear function of energy. Any relation with the orthogonal polynomials {p n } is eluded. Alexander Moroz (WS) Recursive solution January 3, / 37
53 Comparison with Braak s PRL 107, (2011) - II The TTRR for the coefficients K n is not a fundamental one, i.e. not of Haydock s type, because of f n (ζ) is not a linear function of energy. Any relation with the orthogonal polynomials {p n } is eluded. Instead at looking straight at the zeros of p N of sufficiently large degree, a complicated composite object, G ± (ζ), is formed. Alexander Moroz (WS) Recursive solution January 3, / 37
54 Is Rabi model integrable or solvable? Alexander Moroz (WS) Recursive solution January 3, / 37
55 Structure relation and bispectral condition Divided-difference operator or discrete derivative D x f (x) = f (ι +(x)) f (ι (x)) ι + (x) ι (x) Alexander Moroz (WS) Recursive solution January 3, / 37
56 Structure relation and bispectral condition Divided-difference operator or discrete derivative Structure relation D x f (x) = f (ι +(x)) f (ι (x)) ι + (x) ι (x) D x p n (x) = B n (x)p n (x) + A n (x)p n 1 (x) where D x : Π n [x] Π n 1 [x], Π n [x] being the linear space of polynomials in x over C with degree at most n Z 0 Alexander Moroz (WS) Recursive solution January 3, / 37
57 Consequences of a structure relation if there is one structure relation there is another - simply apply the fundamental TTRR Alexander Moroz (WS) Recursive solution January 3, / 37
58 Consequences of a structure relation if there is one structure relation there is another - simply apply the fundamental TTRR a pair of mutually adjoint raising and lowering ladder operators Alexander Moroz (WS) Recursive solution January 3, / 37
59 Consequences of a structure relation if there is one structure relation there is another - simply apply the fundamental TTRR a pair of mutually adjoint raising and lowering ladder operators orthogonal polynomials satisfy in general a second-order difference equation = bispectral property Alexander Moroz (WS) Recursive solution January 3, / 37
60 Consequences of a structure relation if there is one structure relation there is another - simply apply the fundamental TTRR a pair of mutually adjoint raising and lowering ladder operators orthogonal polynomials satisfy in general a second-order difference equation = bispectral property allows one to introduce a discrete analogue of the Bethe Ansatz equations Alexander Moroz (WS) Recursive solution January 3, / 37
61 Consequences of a structure relation if there is one structure relation there is another - simply apply the fundamental TTRR a pair of mutually adjoint raising and lowering ladder operators orthogonal polynomials satisfy in general a second-order difference equation = bispectral property allows one to introduce a discrete analogue of the Bethe Ansatz equations the polynomials are hypergeometric of Askey-Wilson type Alexander Moroz (WS) Recursive solution January 3, / 37
62 D x requires four primary classes of special non-uniform lattices 1 the linear lattice Alexander Moroz (WS) Recursive solution January 3, / 37
63 D x requires four primary classes of special non-uniform lattices 1 the linear lattice 2 the linear q-lattice Alexander Moroz (WS) Recursive solution January 3, / 37
64 D x requires four primary classes of special non-uniform lattices 1 the linear lattice 2 the linear q-lattice 3 the quadratic lattice Alexander Moroz (WS) Recursive solution January 3, / 37
65 D x requires four primary classes of special non-uniform lattices 1 the linear lattice 2 the linear q-lattice 3 the quadratic lattice 4 the q-quadratic lattice Alexander Moroz (WS) Recursive solution January 3, / 37
66 D x requires four primary classes of special non-uniform lattices 1 the linear lattice 2 the linear q-lattice 3 the quadratic lattice 4 the q-quadratic lattice Either Λ = {x x = u 2 n 2 + u 1 n + u 0, n N} Alexander Moroz (WS) Recursive solution January 3, / 37
67 D x requires four primary classes of special non-uniform lattices 1 the linear lattice 2 the linear q-lattice 3 the quadratic lattice 4 the q-quadratic lattice Either or Λ = {x x = u 2 n 2 + u 1 n + u 0, n N} Λ = {x x = u 2 q n + u 1 q n + u 0, n N} where u j are real constants and the real parameter q satisfies 0 < q < 1 Alexander Moroz (WS) Recursive solution January 3, / 37
68 Is Rabi model integrable or solvable? Alexander Moroz (WS) Recursive solution January 3, / 37
69 Characterization of tridiagonal He n 1 = a n 1 e n 1 + b n e n + b n 1 e n 2 (A) b n = 0 for some n = N > 0 Alexander Moroz (WS) Recursive solution January 3, / 37
70 Characterization of tridiagonal He n 1 = a n 1 e n 1 + b n e n + b n 1 e n 2 (A) b n = 0 for some n = N > 0 (B) b n 0 for all n 0 Alexander Moroz (WS) Recursive solution January 3, / 37
71 He n 1 = a n 1 e n 1 + b n e n + b n 1 e n 2 with b N = 0 for some N > 0 p N (E) enters the TTRR first for n = N. One has p n+n (E) = p N (E)Q n (E) (n 0) i.e. p N +1 (E) is proportional to p N (E). The quotient polynomials Q n (E) form a new orthogonal sequence of polynomials. Alexander Moroz (WS) Recursive solution January 3, / 37
72 He n 1 = a n 1 e n 1 + b n e n + b n 1 e n 2 with b N = 0 for some N > 0 p N (E) enters the TTRR first for n = N. One has p n+n (E) = p N (E)Q n (E) (n 0) i.e. p N +1 (E) is proportional to p N (E). The quotient polynomials Q n (E) form a new orthogonal sequence of polynomials. The polynomials {p n (E)} N 1 n=0 decouple from {p n(e)} n=n. Alexander Moroz (WS) Recursive solution January 3, / 37
73 He n 1 = a n 1 e n 1 + b n e n + b n 1 e n 2 with b N = 0 for some N > 0 p N (E) enters the TTRR first for n = N. One has p n+n (E) = p N (E)Q n (E) (n 0) i.e. p N +1 (E) is proportional to p N (E). The quotient polynomials Q n (E) form a new orthogonal sequence of polynomials. The polynomials {p n (E)} N n=0 1 decouple from {p n(e)} n=n. The operator Ĥ has necessary a finite dimensional invariant subspace V N = {e n } N n=0 1. Alexander Moroz (WS) Recursive solution January 3, / 37
74 Quasi-exact solvabality of a linear differential operator Ĥ Ĥ preserves a finite dimensional subspace V N of a Hilbert space L 2 (S) (S R), ĤV N V N, dim V N = N <, on which the operator Ĥ is naturally defined Alexander Moroz (WS) Recursive solution January 3, / 37
75 Quasi-exact solvabality of a linear differential operator Ĥ Ĥ preserves a finite dimensional subspace V N of a Hilbert space L 2 (S) (S R), ĤV N V N, dim V N = N <, on which the operator Ĥ is naturally defined The basis of V N admits an explicit analytic form V N = span {e 0 (x),..., e N 1 (x)} Alexander Moroz (WS) Recursive solution January 3, / 37
76 Basic properties of finite OPS defined by xp n = P n+1 + c n P n + λ n P n 1 If λ k > 0 and c k are real for k = 0, 1,..., N 1 then: Alexander Moroz (WS) Recursive solution January 3, / 37
77 Basic properties of finite OPS defined by xp n = P n+1 + c n P n + λ n P n 1 If λ k > 0 and c k are real for k = 0, 1,..., N 1 then: the polynomials of the resulting finite orthogonal polynomial sequence {P k } N k=1 have real and simple zeros Alexander Moroz (WS) Recursive solution January 3, / 37
78 Basic properties of finite OPS defined by xp n = P n+1 + c n P n + λ n P n 1 If λ k > 0 and c k are real for k = 0, 1,..., N 1 then: the polynomials of the resulting finite orthogonal polynomial sequence {P k } N k=1 have real and simple zeros the zeros of {P k } N k=1 interlace (the Sturm property) Alexander Moroz (WS) Recursive solution January 3, / 37
79 Basic properties of finite OPS defined by xp n = P n+1 + c n P n + λ n P n 1 If λ k > 0 and c k are real for k = 0, 1,..., N 1 then: the polynomials of the resulting finite orthogonal polynomial sequence {P k } N k=1 have real and simple zeros the zeros of {P k } N k=1 interlace (the Sturm property) ν P (E) = N 1 k=0 w k θ(e E k ), where θ(x) is Heaviside s step function, E k are the N zeros of P N, and the coefficients w k can be found by solving N 1 l=0 P k (E l ) w l = δ k0, k = 0, 1,..., N 1. Alexander Moroz (WS) Recursive solution January 3, / 37
80 Basic properties of finite OPS defined by xp n = P n+1 + c n P n + λ n P n 1 If λ k > 0 and c k are real for k = 0, 1,..., N 1 then: the polynomials of the resulting finite orthogonal polynomial sequence {P k } N k=1 have real and simple zeros the zeros of {P k } N k=1 interlace (the Sturm property) ν P (E) = N 1 k=0 w k θ(e E k ), where θ(x) is Heaviside s step function, E k are the N zeros of P N, and the coefficients w k can be found by solving N 1 l=0 w k > 0 for 0 k < N P k (E l ) w l = δ k0, k = 0, 1,..., N 1. Alexander Moroz (WS) Recursive solution January 3, / 37
81 Basic properties of finite OPS defined by xp n = P n+1 + c n P n + λ n P n 1 If λ k > 0 and c k are real for k = 0, 1,..., N 1 then: the polynomials of the resulting finite orthogonal polynomial sequence {P k } N k=1 have real and simple zeros the zeros of {P k } N k=1 interlace (the Sturm property) ν P (E) = N 1 k=0 w k θ(e E k ), where θ(x) is Heaviside s step function, E k are the N zeros of P N, and the coefficients w k can be found by solving N 1 l=0 w k > 0 for 0 k < N dν P (E) is unique P k (E l ) w l = δ k0, k = 0, 1,..., N 1. Alexander Moroz (WS) Recursive solution January 3, / 37
82 Normalization paradox Orthogonality L[π(x)p n (x)] = 0 for every polynomial π(x) of degree k < n. Hence the norm of p n (E) for n N vanishes, i.e. 0 1 b N +1 L[E N p N (E)] Paradox One cannot satisfy (H1), i.e. that {p n (E)} are orthonormal with respect to the density of states (DOS), n 0 (E), p n (E)p m (E)n 0 (E)dE = δ nm Alexander Moroz (WS) Recursive solution January 3, / 37
83 Resolution Assume that the respective spectra obtained in V N and L 2 (S)\V N form two disjoint intervals on the real axis, S = S qes S, sup S qes < inf S. Alexander Moroz (WS) Recursive solution January 3, / 37
84 Resolution Assume that the respective spectra obtained in V N and L 2 (S)\V N form two disjoint intervals on the real axis, S = S qes S, sup S qes < inf S. On S qes is n 0 (E)dE represented by the discrete measure dν P (E) Alexander Moroz (WS) Recursive solution January 3, / 37
85 Resolution Assume that the respective spectra obtained in V N and L 2 (S)\V N form two disjoint intervals on the real axis, S = S qes S, sup S qes < inf S. On S qes is n 0 (E)dE represented by the discrete measure dν P (E) On S is n 0 (E)dE represented by dν PQ (E) = 1 p 2 N (E) dν Q(E) (E S ) Alexander Moroz (WS) Recursive solution January 3, / 37
86 Conclusions tridiagonality is fundamental Alexander Moroz (WS) Recursive solution January 3, / 37
87 Conclusions tridiagonality is fundamental TTRR are indispensable Alexander Moroz (WS) Recursive solution January 3, / 37
88 Conclusions tridiagonality is fundamental TTRR are indispensable OPS are unavoidable Alexander Moroz (WS) Recursive solution January 3, / 37
89 Conclusions tridiagonality is fundamental TTRR are indispensable OPS are unavoidable an efficient numerical method for energy levels Alexander Moroz (WS) Recursive solution January 3, / 37
90 Conclusions tridiagonality is fundamental TTRR are indispensable OPS are unavoidable an efficient numerical method for energy levels characterization of integrable and QES models in terms of intrinsic polynomial properties Alexander Moroz (WS) Recursive solution January 3, / 37
91 Conclusions tridiagonality is fundamental TTRR are indispensable OPS are unavoidable an efficient numerical method for energy levels characterization of integrable and QES models in terms of intrinsic polynomial properties confine integrable spectra to four basic classes Alexander Moroz (WS) Recursive solution January 3, / 37
92 Conclusions tridiagonality is fundamental TTRR are indispensable OPS are unavoidable an efficient numerical method for energy levels characterization of integrable and QES models in terms of intrinsic polynomial properties confine integrable spectra to four basic classes characterization of chaos? Alexander Moroz (WS) Recursive solution January 3, / 37
93 R. Haydock, p. 217 of his review: The generality of this result suggests a new way of looking at quantum mechanics. Since any quantum system can be transformed into a chain model, we need only investigate such chain models in order to see the varieties of quantum phenomena that are possible. To understand a particular physical system, we need only find the chain model appropriate to it. Alexander Moroz (WS) Recursive solution January 3, / 37
94 The End Alexander Moroz (WS) Recursive solution January 3, / 37
95 Further reading This talk has been based on excerpts from my arxiv: ; EPL 100, (2012); AP 338, (2013); AP 340, (2014); J. Phys. A: Math. Theor. 47(49), (2014); AP 351, (2014) Visit me at F77 source codes and numerical data files for the Rabi model can be freely downloaded from My all F77 scattering codes are available at For list of my unfinished projects see Alexander Moroz (WS) Recursive solution January 3, / 37
The 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 informationQuantum wires, orthogonal polynomials and Diophantine approximation
Quantum wires, orthogonal polynomials and Diophantine approximation Introduction Quantum Mechanics (QM) is a linear theory Main idea behind Quantum Information (QI): use the superposition principle of
More informationwhich implies that we can take solutions which are simultaneous eigen functions of
Module 1 : Quantum Mechanics Chapter 6 : Quantum mechanics in 3-D Quantum mechanics in 3-D For most physical systems, the dynamics is in 3-D. The solutions to the general 3-d problem are quite complicated,
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 informationLecture 4 Quantum mechanics in more than one-dimension
Lecture 4 Quantum mechanics in more than one-dimension Background Previously, we have addressed quantum mechanics of 1d systems and explored bound and unbound (scattering) states. Although general concepts
More informationBloch Wilson Hamiltonian and a Generalization of the Gell-Mann Low Theorem 1
Bloch Wilson Hamiltonian and a Generalization of the Gell-Mann Low Theorem 1 Axel Weber 2 arxiv:hep-th/9911198v1 25 Nov 1999 Instituto de Física y Matemáticas, Universidad Michoacana de San Nicolás de
More informationIsotropic harmonic oscillator
Isotropic harmonic oscillator 1 Isotropic harmonic oscillator The hamiltonian of the isotropic harmonic oscillator is H = h m + 1 mω r (1) = [ h d m dρ + 1 ] m ω ρ, () ρ=x,y,z a sum of three one-dimensional
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 information(a) [15] Use perturbation theory to compute the eigenvalues and eigenvectors of H. Compute all terms up to fourth-order in Ω. The bare eigenstates are
PHYS85 Quantum Mechanics II, Spring 00 HOMEWORK ASSIGNMENT 6: Solutions Topics covered: Time-independent perturbation theory.. [0 Two-Level System: Consider the system described by H = δs z + ΩS x, with
More information1 Mathematical preliminaries
1 Mathematical preliminaries The mathematical language of quantum mechanics is that of vector spaces and linear algebra. In this preliminary section, we will collect the various definitions and mathematical
More informationPhysics 221A Fall 2017 Notes 27 The Variational Method
Copyright c 2018 by Robert G. Littlejohn Physics 221A Fall 2017 Notes 27 The Variational Method 1. Introduction Very few realistic problems in quantum mechanics are exactly solvable, so approximation methods
More informationLecture 4 Quantum mechanics in more than one-dimension
Lecture 4 Quantum mechanics in more than one-dimension Background Previously, we have addressed quantum mechanics of 1d systems and explored bound and unbound (scattering) states. Although general concepts
More informationErik A. van Doorn Department of Applied Mathematics University of Twente Enschede, The Netherlands
REPRESENTATIONS FOR THE DECAY PARAMETER OF A BIRTH-DEATH PROCESS Erik A. van Doorn Department of Applied Mathematics University of Twente Enschede, The Netherlands University of Leeds 30 October 2014 Acknowledgment:
More informationNormal form for the non linear Schrödinger equation
Normal form for the non linear Schrödinger equation joint work with Claudio Procesi and Nguyen Bich Van Universita di Roma La Sapienza S. Etienne de Tinee 4-9 Feb. 2013 Nonlinear Schrödinger equation Consider
More informationREFINMENTS OF HOMOLOGY. Homological spectral package.
REFINMENTS OF HOMOLOGY (homological spectral package). Dan Burghelea Department of Mathematics Ohio State University, Columbus, OH In analogy with linear algebra over C we propose REFINEMENTS of the simplest
More informationThe Simple Harmonic Oscillator
The Simple Harmonic Oscillator Asaf Pe er 1 November 4, 215 This part of the course is based on Refs [1] [3] 1 Introduction We return now to the study of a 1-d stationary problem: that of the simple harmonic
More informationMagnetic wells in dimension three
Magnetic wells in dimension three Yuri A. Kordyukov joint with Bernard Helffer & Nicolas Raymond & San Vũ Ngọc Magnetic Fields and Semiclassical Analysis Rennes, May 21, 2015 Yuri A. Kordyukov (Ufa) Magnetic
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 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 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 information3 Symmetry Protected Topological Phase
Physics 3b Lecture 16 Caltech, 05/30/18 3 Symmetry Protected Topological Phase 3.1 Breakdown of noninteracting SPT phases with interaction Building on our previous discussion of the Majorana chain and
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 informationProblem Set 6: Solutions Math 201A: Fall a n x n,
Problem Set 6: Solutions Math 201A: Fall 2016 Problem 1. Is (x n ) n=0 a Schauder basis of C([0, 1])? No. If f(x) = a n x n, n=0 where the series converges uniformly on [0, 1], then f has a power series
More informationNANOSCALE SCIENCE & TECHNOLOGY
. NANOSCALE SCIENCE & TECHNOLOGY V Two-Level Quantum Systems (Qubits) Lecture notes 5 5. Qubit description Quantum bit (qubit) is an elementary unit of a quantum computer. Similar to classical computers,
More informationIntroduction to orthogonal polynomials. Michael Anshelevich
Introduction to orthogonal polynomials Michael Anshelevich November 6, 2003 µ = probability measure on R with finite moments m n (µ) = R xn dµ(x)
More informationLecture 3 Dynamics 29
Lecture 3 Dynamics 29 30 LECTURE 3. DYNAMICS 3.1 Introduction Having described the states and the observables of a quantum system, we shall now introduce the rules that determine their time evolution.
More informationMATH 31 - ADDITIONAL PRACTICE PROBLEMS FOR FINAL
MATH 3 - ADDITIONAL PRACTICE PROBLEMS FOR FINAL MAIN TOPICS FOR THE FINAL EXAM:. Vectors. Dot product. Cross product. Geometric applications. 2. Row reduction. Null space, column space, row space, left
More informationThe 3 dimensional Schrödinger Equation
Chapter 6 The 3 dimensional Schrödinger Equation 6.1 Angular Momentum To study how angular momentum is represented in quantum mechanics we start by reviewing the classical vector of orbital angular momentum
More informationIntroduction to the Mathematics of the XY -Spin Chain
Introduction to the Mathematics of the XY -Spin Chain Günter Stolz June 9, 2014 Abstract In the following we present an introduction to the mathematical theory of the XY spin chain. The importance of this
More informationIntroduction to Ergodic Theory
Introduction to Ergodic Theory Marius Lemm May 20, 2010 Contents 1 Ergodic stochastic processes 2 1.1 Canonical probability space.......................... 2 1.2 Shift operators.................................
More informationName: Final Exam MATH 3320
Name: Final Exam MATH 3320 Directions: Make sure to show all necessary work to receive full credit. If you need extra space please use the back of the sheet with appropriate labeling. (1) State the following
More informationLecture 10. The Dirac equation. WS2010/11: Introduction to Nuclear and Particle Physics
Lecture 10 The Dirac equation WS2010/11: Introduction to Nuclear and Particle Physics The Dirac equation The Dirac equation is a relativistic quantum mechanical wave equation formulated by British physicist
More informationSolutions: Problem Set 3 Math 201B, Winter 2007
Solutions: Problem Set 3 Math 201B, Winter 2007 Problem 1. Prove that an infinite-dimensional Hilbert space is a separable metric space if and only if it has a countable orthonormal basis. Solution. If
More informationOn some tridiagonal k Toeplitz matrices: algebraic and analytical aspects. Applications
On some tridiagonal k Toeplitz matrices: algebraic and analytical aspects Applications R Álvarez-Nodarse1,, J Petronilho 3, and NR Quintero 4, 1 Dpto Análisis Matemático, Universidad de Sevilla, Apdo 1160,
More informationMultiple Meixner polynomials and non-hermitian oscillator Hamiltonians arxiv: v2 [math.ca] 8 Nov 2013
Multiple Meixner polynomials and non-hermitian oscillator Hamiltonians arxiv:1310.0982v2 [math.ca] 8 Nov 2013 F Ndayiragie 1 and W Van Assche 2 1 Département de Mathématiques, Université du Burundi, Campus
More informationMassachusetts Institute of Technology Physics Department
Massachusetts Institute of Technology Physics Department Physics 8.32 Fall 2006 Quantum Theory I October 9, 2006 Assignment 6 Due October 20, 2006 Announcements There will be a makeup lecture on Friday,
More informationProblem 1: A 3-D Spherical Well(10 Points)
Problem : A 3-D Spherical Well( Points) For this problem, consider a particle of mass m in a three-dimensional spherical potential well, V (r), given as, V = r a/2 V = W r > a/2. with W >. All of the following
More informationarxiv: v1 [quant-ph] 8 Sep 2010
Few-Body Systems, (8) Few- Body Systems c by Springer-Verlag 8 Printed in Austria arxiv:9.48v [quant-ph] 8 Sep Two-boson Correlations in Various One-dimensional Traps A. Okopińska, P. Kościk Institute
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 informationPhysics 221A Fall 1996 Notes 12 Orbital Angular Momentum and Spherical Harmonics
Physics 221A Fall 1996 Notes 12 Orbital Angular Momentum and Spherical Harmonics We now consider the spatial degrees of freedom of a particle moving in 3-dimensional space, which of course is an important
More information5 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 informationTheoretical design of a readout system for the Flux Qubit-Resonator Rabi Model in the ultrastrong coupling regime
Theoretical design of a readout system for the Flux Qubit-Resonator Rabi Model in the ultrastrong coupling regime Ceren Burçak Dağ Supervisors: Dr. Pol Forn-Díaz and Assoc. Prof. Christopher Wilson Institute
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 informationSpectral 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(Dynamical) quantum typicality: What is it and what are its physical and computational implications?
(Dynamical) : What is it and what are its physical and computational implications? Jochen Gemmer University of Osnabrück, Kassel, May 13th, 214 Outline Thermal relaxation in closed quantum systems? Typicality
More informationTime Independent Perturbation Theory Contd.
Time Independent Perturbation Theory Contd. A summary of the machinery for the Perturbation theory: H = H o + H p ; H 0 n >= E n n >; H Ψ n >= E n Ψ n > E n = E n + E n ; E n = < n H p n > + < m H p n
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 informationAlex Eremenko, Andrei Gabrielov. Purdue University. eremenko. agabriel
SPECTRAL LOCI OF STURM LIOUVILLE OPERATORS WITH POLYNOMIAL POTENTIALS Alex Eremenko, Andrei Gabrielov Purdue University www.math.purdue.edu/ eremenko www.math.purdue.edu/ agabriel Kharkov, August 2012
More informationQuantum Mechanics Solutions
Quantum Mechanics Solutions (a (i f A and B are Hermitian, since (AB = B A = BA, operator AB is Hermitian if and only if A and B commute So, we know that [A,B] = 0, which means that the Hilbert space H
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 informationNewton s Method and Localization
Newton s Method and Localization Workshop on Analytical Aspects of Mathematical Physics John Imbrie May 30, 2013 Overview Diagonalizing the Hamiltonian is a goal in quantum theory. I would like to discuss
More informationPreliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012
Instructions Preliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012 The exam consists of four problems, each having multiple parts. You should attempt to solve all four problems. 1.
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 informationMATH 56A: STOCHASTIC PROCESSES CHAPTER 2
MATH 56A: STOCHASTIC PROCESSES CHAPTER 2 2. Countable Markov Chains I started Chapter 2 which talks about Markov chains with a countably infinite number of states. I did my favorite example which is on
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 informationOverview of Topological Cluster-State Quantum Computation on 2D Cluster-State
Overview of Topological Cluster-State Quantum Computation on 2D Cluster-State based on High-threshold universal quantum computation on the surface code -Austin G. Fowler, Ashley M. Stephens, and Peter
More informationQUANTUM THEORY OF LIGHT EECS 638/PHYS 542/AP609 FINAL EXAMINATION
Instructor: Professor S.C. Rand Date: April 5 001 Duration:.5 hours QUANTUM THEORY OF LIGHT EECS 638/PHYS 54/AP609 FINAL EXAMINATION PLEASE read over the entire examination before you start. DO ALL QUESTIONS
More informationHarmonic Oscillator Eigenvalues and Eigenfunctions
Chemistry 46 Fall 217 Dr. Jean M. Standard October 4, 217 Harmonic Oscillator Eigenvalues and Eigenfunctions The Quantum Mechanical Harmonic Oscillator The quantum mechanical harmonic oscillator in one
More informationThe Phase Space in Quantum Field Theory
The Phase Space in Quantum Field Theory How Small? How Large? Martin Porrmann II. Institut für Theoretische Physik Universität Hamburg Seminar Quantum Field Theory and Mathematical Physics April 13, 2005
More informationBi-Orthogonal Approach to Non-Hermitian Hamiltonians with the Oscillator Spectrum: Generalized Coherent States for Nonlinear Algebras
Bi-Orthogonal Approach to Non-Hermitian Hamiltonians with the Oscillator Spectrum: Generalized Coherent States for Nonlinear Algebras arxiv:1706.04684v [quant-ph] 3 Oct 017 Oscar Rosas-Ortiz and Kevin
More informationLecture on: Numerical sparse linear algebra and interpolation spaces. June 3, 2014
Lecture on: Numerical sparse linear algebra and interpolation spaces June 3, 2014 Finite dimensional Hilbert spaces and IR N 2 / 38 (, ) : H H IR scalar product and u H = (u, u) u H norm. Finite dimensional
More informationApplied Physics 150a: Homework #3
Applied Physics 150a: Homework #3 (Dated: November 13, 2014) Due: Thursday, November 20th, anytime before midnight. There will be an INBOX outside my office in Watson (Rm. 266/268). 1. (10 points) The
More informationFinding eigenvalues for matrices acting on subspaces
Finding eigenvalues for matrices acting on subspaces Jakeniah Christiansen Department of Mathematics and Statistics Calvin College Grand Rapids, MI 49546 Faculty advisor: Prof Todd Kapitula Department
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 informationLecture #1. Review. Postulates of quantum mechanics (1-3) Postulate 1
L1.P1 Lecture #1 Review Postulates of quantum mechanics (1-3) Postulate 1 The state of a system at any instant of time may be represented by a wave function which is continuous and differentiable. Specifically,
More informationHomework 4, 5, 6 Solutions. > 0, and so a n 0 = n + 1 n = ( n+1 n)( n+1+ n) 1 if n is odd 1/n if n is even diverges.
2..2(a) lim a n = 0. Homework 4, 5, 6 Solutions Proof. Let ɛ > 0. Then for n n = 2+ 2ɛ we have 2n 3 4+ ɛ 3 > ɛ > 0, so 0 < 2n 3 < ɛ, and thus a n 0 = 2n 3 < ɛ. 2..2(g) lim ( n + n) = 0. Proof. Let ɛ >
More informationDegenerate Perturbation Theory. 1 General framework and strategy
Physics G6037 Professor Christ 12/22/2015 Degenerate Perturbation Theory The treatment of degenerate perturbation theory presented in class is written out here in detail. The appendix presents the underlying
More informationCentral density. Consider nuclear charge density. Frois & Papanicolas, Ann. Rev. Nucl. Part. Sci. 37, 133 (1987) QMPT 540
Central density Consider nuclear charge density Frois & Papanicolas, Ann. Rev. Nucl. Part. Sci. 37, 133 (1987) Central density (A/Z* charge density) about the same for nuclei heavier than 16 O, corresponding
More information1 Math 241A-B Homework Problem List for F2015 and W2016
1 Math 241A-B Homework Problem List for F2015 W2016 1.1 Homework 1. Due Wednesday, October 7, 2015 Notation 1.1 Let U be any set, g be a positive function on U, Y be a normed space. For any f : U Y let
More informationA Fourier-Bessel Expansion for Solving Radial Schrödinger Equation in Two Dimensions
A Fourier-Bessel Expansion for Solving Radial Schrödinger Equation in Two Dimensions H. TAŞELI and A. ZAFER Department of Mathematics, Middle East Technical University, 06531 Ankara, Turkey ABSTRACT The
More informationA classification of sharp tridiagonal pairs. Tatsuro Ito, Kazumasa Nomura, Paul Terwilliger
Tatsuro Ito Kazumasa Nomura Paul Terwilliger Overview This talk concerns a linear algebraic object called a tridiagonal pair. We will describe its features such as the eigenvalues, dual eigenvalues, shape,
More informationThe infinite square well in a reformulation of quantum mechanics without potential function
The infinite square well in a reformulation of quantum mechanics without potential function A.D. Alhaidari (a), T.J. Taiwo (b) (a) Saudi Center for Theoretical Physics, P.O Box 32741, Jeddah 21438, Saudi
More informationLinear Algebra. Min Yan
Linear Algebra Min Yan January 2, 2018 2 Contents 1 Vector Space 7 1.1 Definition................................. 7 1.1.1 Axioms of Vector Space..................... 7 1.1.2 Consequence of Axiom......................
More informationSPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS
SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS TSOGTGEREL GANTUMUR Abstract. After establishing discrete spectra for a large class of elliptic operators, we present some fundamental spectral properties
More informationThe quantum state as a vector
The quantum state as a vector February 6, 27 Wave mechanics In our review of the development of wave mechanics, we have established several basic properties of the quantum description of nature:. A particle
More informationContents. Preface for the Instructor. Preface for the Student. xvii. Acknowledgments. 1 Vector Spaces 1 1.A R n and C n 2
Contents Preface for the Instructor xi Preface for the Student xv Acknowledgments xvii 1 Vector Spaces 1 1.A R n and C n 2 Complex Numbers 2 Lists 5 F n 6 Digression on Fields 10 Exercises 1.A 11 1.B Definition
More informationQuantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours.
Quantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours. There are 10 problems, totalling 180 points. Do all problems. Answer all problems in the white books provided.
More informationSolutions to chapter 4 problems
Chapter 9 Solutions to chapter 4 problems Solution to Exercise 47 For example, the x component of the angular momentum is defined as ˆL x ŷˆp z ẑ ˆp y The position and momentum observables are Hermitian;
More informationPart V. 17 Introduction: What are measures and why measurable sets. Lebesgue Integration Theory
Part V 7 Introduction: What are measures and why measurable sets Lebesgue Integration Theory Definition 7. (Preliminary). A measure on a set is a function :2 [ ] such that. () = 2. If { } = is a finite
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 informationProblem Set No. 3: Canonical Quantization Due Date: Wednesday October 19, 2018, 5:00 pm. 1 Spin waves in a quantum Heisenberg antiferromagnet
Physics 58, Fall Semester 018 Professor Eduardo Fradkin Problem Set No. 3: Canonical Quantization Due Date: Wednesday October 19, 018, 5:00 pm 1 Spin waves in a quantum Heisenberg antiferromagnet In this
More informationC. Show your answer in part B agrees with your answer in part A in the limit that the constant c 0.
Problem #1 A. A projectile of mass m is shot vertically in the gravitational field. Its initial velocity is v o. Assuming there is no air resistance, how high does m go? B. Now assume the projectile is
More informationNonlinear Integral Equation Formulation of Orthogonal Polynomials
Nonlinear Integral Equation Formulation of Orthogonal Polynomials Eli Ben-Naim Theory Division, Los Alamos National Laboratory with: Carl Bender (Washington University, St. Louis) C.M. Bender and E. Ben-Naim,
More informationPhysics 4022 Notes on Density Matrices
Physics 40 Notes on Density Matrices Definition: For a system in a definite normalized state ψ > the density matrix ρ is ρ = ψ >< ψ 1) From Eq 1 it is obvious that in the basis defined by ψ > and other
More informationLecture 11: Clifford algebras
Lecture 11: Clifford algebras In this lecture we introduce Clifford algebras, which will play an important role in the rest of the class. The link with K-theory is the Atiyah-Bott-Shapiro construction
More informationStructured Krylov Subspace Methods for Eigenproblems with Spectral Symmetries
Structured Krylov Subspace Methods for Eigenproblems with Spectral Symmetries Fakultät für Mathematik TU Chemnitz, Germany Peter Benner benner@mathematik.tu-chemnitz.de joint work with Heike Faßbender
More informationin terms of the classical frequency, ω = , puts the classical Hamiltonian in the form H = p2 2m + mω2 x 2
One of the most important problems in quantum mechanics is the simple harmonic oscillator, in part because its properties are directly applicable to field theory. The treatment in Dirac notation is particularly
More informationA new approach to quantum metrics. Nik Weaver. (joint work with Greg Kuperberg, in progress)
A new approach to quantum metrics Nik Weaver (joint work with Greg Kuperberg, in progress) Definition. A dual operator system is a linear subspace V of B(H) such that I V A V implies A V V is weak*-closed.
More informationWeighted Sums of Orthogonal Polynomials Related to Birth-Death Processes with Killing
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 8, Number 2, pp. 401 412 (2013) http://campus.mst.edu/adsa Weighted Sums of Orthogonal Polynomials Related to Birth-Death Processes
More informationSpectral Theorem for Self-adjoint Linear Operators
Notes for the undergraduate lecture by David Adams. (These are the notes I would write if I was teaching a course on this topic. I have included more material than I will cover in the 45 minute lecture;
More information5.4 Given the basis e 1, e 2 write the matrices that represent the unitary transformations corresponding to the following changes of basis:
5 Representations 5.3 Given a three-dimensional Hilbert space, consider the two observables ξ and η that, with respect to the basis 1, 2, 3, arerepresentedby the matrices: ξ ξ 1 0 0 0 ξ 1 0 0 0 ξ 3, ξ
More informationCHM 532 Notes on Creation and Annihilation Operators
CHM 53 Notes on Creation an Annihilation Operators These notes provie the etails concerning the solution to the quantum harmonic oscillator problem using the algebraic metho iscusse in class. The operators
More informationA PRIMER ON SESQUILINEAR FORMS
A PRIMER ON SESQUILINEAR FORMS BRIAN OSSERMAN This is an alternative presentation of most of the material from 8., 8.2, 8.3, 8.4, 8.5 and 8.8 of Artin s book. Any terminology (such as sesquilinear form
More informationHighest-weight Theory: Verma Modules
Highest-weight Theory: Verma Modules Math G4344, Spring 2012 We will now turn to the problem of classifying and constructing all finitedimensional representations of a complex semi-simple Lie algebra (or,
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 informationMATH 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 informationarxiv:quant-ph/ v1 9 Oct 2006
Adiabatic Elimination in a Lambda System E. Brion, L.H. Pedersen, and K. Mølmer Lundbeck Foundation Theoretical Center for Quantum System Research, Department of Physics and Astronomy, University of Aarhus,
More informationNormalization integrals of orthogonal Heun functions
Normalization integrals of orthogonal Heun functions Peter A. Becker a) Center for Earth Observing and Space Research, Institute for Computational Sciences and Informatics, and Department of Physics and
More informationIncompatibility Paradoxes
Chapter 22 Incompatibility Paradoxes 22.1 Simultaneous Values There is never any difficulty in supposing that a classical mechanical system possesses, at a particular instant of time, precise values of
More informationNotes on Quantum Mechanics
Notes on Quantum Mechanics Kevin S. Huang Contents 1 The Wave Function 1 1.1 The Schrodinger Equation............................ 1 1. Probability.................................... 1.3 Normalization...................................
More information