sources of program support: MRSEC Toeplitz at Brandeis April
|
|
- Elmer Gardner
- 6 years ago
- Views:
Transcription
1 Landau Quasiparticle Methods for Calculating Eigenvalues and Eigenvectors of Large Matrices Hui Dai, Zachary Geary, and Leo Kadanoff # # LeoP@UChicago.edu Matrices are used to describe the modes of oscillation of physical objects. Generally, N by N matrices describe objects with N modes. We often calculate the eigenvectors of large matrices of this kind by making use of calculations for infinite matrices. This approach follows the spirit of Landau quasiparticle methods and also the kinds of scalingt used in critical phenomena. Propagation of information over the whole range of spatial indices in these eigenvectors is demonstrated by power law and logarithmic N-dependence of the eigenvalues and eigenvectors. The dependences can be interpreted in terms of the transfer of information along the line of objects described by the matrix. sources of program support: MRSEC Toeplitz at Brandeis April
2 Toeplitz at Brandeis April
3 Q: How do I find a worthwhile problem? Toeplitz at Brandeis April
4 Q: How do I find a worthwhile problem? A: Find a subject which will have a lasting impact Toeplitz at Brandeis April
5 Q: How do I find a worthwhile problem? A: Find a subject which will have a lasting impact Q: How can I find one of those? Toeplitz at Brandeis April
6 Q: How do I find a worthwhile problem? A: Find a subject which will have a lasting impact Q: How can I find one of those? A: Find something University simple of Chicago MRSEC, and general which nobody understands Toeplitz at Brandeis April
7 Q: How do I find a worthwhile problem? A: Find a subject which will have a lasting impact Q: How can I find one of those? A: Find something University simple of Chicago MRSEC, and general which nobody understands Q: How can I do that? Toeplitz at Brandeis April
8 Q: How do I find a worthwhile problem? A: Find a subject which will have a lasting impact Q: How can I find one of those? A: Find something University simple of Chicago MRSEC, and general which nobody understands Q: How can I do that? A: That s an interesting problem Toeplitz at Brandeis April
9 what is a Toeplitz matrix? answer: T(j,k)=T(j-k) When does it arise? Toeplitz at Brandeis April
10 Applications of Toeplitz matrices: Oscillation frequencies of molecules in a solid: You can calculate frequencies of oscillations of a regular solid by Fourier transformation A row of N different molecules have frequencies given by the eigenvalues of an N by N Toeplitz matrix A row of N different molecules have frequencies given by the eigenvalues of an N by N Toeplitz matrix Toeplitz at Brandeis April
11 Other applications 2d Ising model: Onsager solution: N by M system Z=det Q N Q is an M by M matrix correlation functions obtained by local changes in Q correlation functions are given as determinants of Toeplitz matrices. Quantum Information in a two-part system: Mutual information is a determinant of Toeplitz matrix average correlations of eigenvalues in ensemble of random matrices: reduced to determinant of Toeplitz matrix Known application of eigenvalues or eigenfunctions: none yet Toeplitz at Brandeis April
12 Other applications 2d Ising model: Onsager solution: N by M system Z=det Q N Q is an M by M matrix correlation functions obtained by local changes in Q correlation functions are given as determinants of Toeplitz matrices. Quantum Information in a two-part system: Mutual information is a determinant of Toeplitz matrix average correlations of eigenvalues in ensemble of random matrices: reduced to determinant of Toeplitz matrix Known application of eigenvalues or eigenfunctions: none yet IT IS GOOD TO COME INTO A SUBJECT EARLY Toeplitz at Brandeis April
13 Toeplitz eigenvalues: formal properties N 1 N 1 k =0 j =0 T(j k)ψ s k (ψ s k ) tr ψ t k k s = ε s ψ s k (ψ s j ) tr T(j k) = ε s (ψ s k ) tr orthonormality = δ s,t (ψ s j ) tr ψ s k = δ j,k if T(j-k) is real. then We arrange eigenstates in order according to the phase of the eigenvalue completeness tr=transpose j -->N-j-1 this transpose symmetry works like a complex conjugation for hermetian matrices {for all Toeplitz matrices) (ψ s j )* = ψ N s 1 j in this case, the transpose symmetry is hermetian conjugation Toeplitz at Brandeis April
14 describe a family of Toeplitz matrices via Fourier Transform: The relation between a and T is a kind of Fourier transform. a is called the symbol. The inverse Fourier transform. is of the form, which mostly makes sense when z is on the unit circle. a(z) = z j j = T(j) We often write z = e ip note that the function a(z) on unit circle describes family of matrices for values of all N Toeplitz at Brandeis April
15 There is a special case of Toeplitz matrix which was used for statistical mechanics, especially critical properties of 2d Ising Model by University McCoy of Chicago MRSEC, and Wu Kadanoff Kadanoff-Ceva special case: Hartwig & Fisher Toeplitz at Brandeis April
16 A Singular Case: The matrices T are singular when the function a is singular on the unit circle. We use the form of Fisher and Hartwig for the symbol a α,β (z)=(2-z-1/z) α (-z) β University of Chicago MRSEC, another IRG way X of writing: a(e -ip ) = (2(1- cos p)) α (-e - ip ) β This form gives for large j-k, β describes jumps in a; α describes zeros T(j-k)~ j-k-β -(1+2 α) Toeplitz at Brandeis April
17 Previous Results Szego looked at det T, which is equivalent to discussing the geometric average of eigenvalues. He proved that, for non-singular a, the large N average has the form (log det T)/N ---> const+ O(1/N) Hartwig and Fisher studied the singular case and found an additional term in the geometric average proportional to (logn)(α 2 β 2 ) / N Harold Widom studied the eigenvalues for large N and proved that in the singular case the distribution of eigenvalues approached the (-infinity,infinity) distribution Lee, Dai, and Bettleheim showed that these eigenvalues each have a correction of the form (logn)f (α,β,z) / N + O(1/ N) Toeplitz at Brandeis April
18 Geometric Strategy Real Problem: on [0, 1,2,...,N-1] Strategy Zero: solve on [-infinity, infinity] (used by previous authors) Strategy One: solve University of Chicago on [0, MRSEC, 1,2, infinity] (this work) Strategy Two: solve on [-infinity,... N-2, N-1] (future work) Strategy Three: patch together One and Two (future work) Toeplitz at Brandeis April
19 Strategy zero replace j,k = {0,1..., N-1} by j,k in set {-infinity, infinity}. Then eigenvalue problem is k = T(j k)ψ k = εψ j This can be solved by simple Fourier transformation. Or simply try a solution of an exponential University form of Chicago MRSEC, ψ j = e ipj p has to be real ε = a(e ip ) Then the eigenvalue is Thus, eigenvalue is determined in (-infinity, infinity) case Following Landau-quasiparticle methods we extend this analysis to include finite-n case Toeplitz at Brandeis April
20 ε = a(e ip ) p has to be real [-, ] example eigenvalues: for z =exp(ip) a(z)=(2-z-1/z) α (-z) β β = 1/ 2 α =1/ 3 Toeplitz at Brandeis April
21 As N--> infinity, finite N eigenvalues converge to (-infinity,infinity) result Eigenvalues for N=40, 100, 200, infinity Calculation performed using Mathematica. Eigenvalues for N=101, 201, 401, 801,1601,infinity * means p is close to π/2 We shall particularly study eigenvectors for these eigenvalues Toeplitz at Brandeis April
22 But first jump ahead to results for eigenvalues: j ε (-infinity, infinity) j ε(0, infinity) 1<β<0 case I: all numbers inside the curve are eigenvalues j ε(0, infinity) 0<β<1 case II: no eigenvalues at all j ε{0, N-1} nonetheless we can use the (0,infinity problem) to get information of the eigenfunctions of the [0,n-1] problem Toeplitz at Brandeis April
23 Example I: Positive Translation : T a(z)=z represent an arbitrary state by eigenvalue equation is T Ψ=εΨ so that ɛψ 0 = 0 ψ j z j+1 = ɛ j=0 j=0 ψ j z j ψ(z) = or j=0 ψ j z j ɛψ j = ψ j 1 so that ψ j+1 = ψ j /ɛ for j=0,1,2 all components of the eigenfunction are zero, so no eigenvalue solution Toeplitz at Brandeis April
24 Example II: Negative Translation : T a(z) = 1/z ψ(z) = ψ j=1 zψ(z) = ψ j z j 1 j=1 eigenvalue equation is T Ψ=εΨ ψ j z j (1/z) z j = z j-1 for j>0 definition (1/z) 1 = 0 or Text so that ψ j+1 = ɛψ j so that for j=0,1,2,... ψ j = ψ 0 ɛ j for j=0,1,2 require ε to be on or inside unit circle Toeplitz at Berkeley January 9,
25 Finite N Results for the eigenvectors Hui Dai, Zachary Geary, Leo Kadanoff The eigenvectors are calculated with the aid of Mathematica Recall: infinite N eigenvectors are of form exp(ipj) with p being real. Hence the magnitude does not vary with j. Large N eigenvectors are roughly of this same form but p has an imaginary part so that the magnitude of the eigenvectors can vary as j varies. In fact, the result looks like Toeplitz at Brandeis April
26 Magnitude of eigenvector finite-n eigenvector roughly varies as exp(ipj) over much of range but p has an imaginary part Empirically Im(p)~ (1+2α) ln N /N drop in eigenfunction over entire range proportional to N -(1+2 α) β = 1/ 2 α =1/ 3 N=101 Toeplitz at Brandeis April
27 Phase of eigenvector increases by approximately (ln i)/i=π/2 as j increases by one unit β = 1/ 2 α =1/ 3 We see that eigenvector is approximately exp(ipj) except at beginning and end Toeplitz at Brandeis April
28 Compare two different N s β = 1/ 2 α =1/ 3 Shape of eigenstate depends only on j/(n-1) except at very beginning and very end. Scale of Eigenstate varies as Toeplitz at Brandeis April
29 Compare two different N s β = 1/ 2 α =1/ 3 Shape of eigenstate depends only on j/(n-1) except at very beginning and very end. Scale of Eigenstate varies as Toeplitz at Brandeis April
30 Scaling Behavior Three Different Regions of Scaling Central Region: eigenfunction depends upon j/(n-1). shape of wave function versus j/(n-1) is independent of N. Four branches visible near j=n-1. beginning points: near j=0. eigenfunction depends upon j. Size of region remains constant as N goes to infinity. ending points: here near j=n-1. eigenfunction depends upon N-j. Size of region remains constant as N goes to infinity. (I think) Toeplitz at Brandeis April
31 Eigenfunction shows N-independent behavior for small j Toeplitz at Brandeis April
32 For larger j, eigenfunction shows an N- independent behavior when it is plotted as a function of j/(n-1) Toeplitz at Brandeis April
33 Next calculate form of eigenvector: Strategy one Instead of considering j,k in {0,1..., N-1} consider j,k in set {0, 1,...,infinity}. Calculate solution in this case by using Weiner-Hopf (complex analytic) technique. existence of solution is determined by winding number = the (signed) number of times a(z) circles the origin when z circles it once in positive sense Case I: Winding number = -1. There is a unique nontrivial eigenfunction. (example a(z)=1/z) Case II: Winding number > -1. There is no nontrivial eigenfunction. (example a(z)=z, winding number = 1) acceptable solution if 0 < α < β < 1 winding number is 1 when 0 < α < β < 1 case I thus no solution Toeplitz at Brandeis April
34 simplified examples: Case I: many eigenvalues T(j k) = δ j k, 1 N 1 k =0 T(j k)ψ k = εψ j a(z)=1/z note that T transfers information to larger index values ψ j +1 = εψ j j=0,1,2...infinity ψ j = Aε j All eigenvalues on or in unit circle are acceptable ψ j = e ipj Im p greater than or equal to zero, therefore, magnitude of ε=1/z=e -ip eigenfunction decreases as j gets larger. Toeplitz at Brandeis April
35 Case II: no solutions T(j k) = δ j k,1 N 1 k =0 T(j k)ψ k εψ 0 = 0 εψ j = ψ j 1 = εψ j a(z)=z j=0,1,2...infinity j=1,2...infinity note that in this case T transfers information to smaller index values ψ j = 0 only trivial solution: eigenfunction is zero Toeplitz at Brandeis April
36 The Calculational method As said earlier, we can calculate the eigenfunction for the case involving the Hartwig-Fisher form of a(z) in the situation in which the spatial indices run from zero to infinity and the values University of Chicago of MRSEC, the parameter are in the right range. In the following, I copy the presentation of our joint work provided by my coauthor, Hui Dai. Toeplitz at Brandeis April
37 Equation can be solved by Fourier transformation Toeplitz at Brandeis April
38 positive and negative powers of z non-negative powers of z negative powers of z Note one equation in two unknowns--solve by using analytic properties Toeplitz at Brandeis April
39 first factorize K(z) = K + (z) /[zk (z)] Kψ + = ψ then separate K + (z)ψ + (z) = zk (z)ψ (z) then solve K + (z)ψ + (z) = zk (z)ψ (z) ψ + = Const / K + = j=0 z j ψ j = const Toeplitz at Brandeis April
40 G + path surrounds all singularities outside unit circle Toeplitz at Brandeis April
41 Pushing the Contour Outward z Start with zero in a(z )-ε (at z =zc) (gives branch cut in exponent) branch cut in a(z ) blue contour at unit circle pole at z =z Convert to branch point, cut, in a(z )-ε branch cut in a(z ) circle at infinity Toeplitz at Brandeis April
42 Result: An evaluation of the eigenfunction terms: = z j ψ j j (z / z c ) j j =0 z c = i e iδ /N e γ(logn )/N a few low powers of z give some N-independent terms for small j term in integrand like (t-1) 2α gives additive term ψ j 1/ j 2α+1 Toeplitz at Brandeis April
43 Weiner Hopf solution Exact Eigenfunction β = 1/ 2 α =1/ 3 Toeplitz at Brandeis April
44 Weiner Hopf solution Exact Eigenfunction β = 1/ 2 α =1/ 3 Toeplitz at Brandeis April
45 Weiner Hopf solution Exact Eigenfunction β = 1/ 2 α =1/ 3 looks the same except for j=100 point Toeplitz at Brandeis April
46 Weiner-Hopf gives excellent solution for j/n<0.9 Plot absolute error in WH approximation vs j/(n-1) Plot relative error in WH approximation vs j/(n-1) Here N=1000, error is of order 10-7 Relative error remains small until large values of j/(n-1) Toeplitz at Brandeis April
47 Eigenfunction is a sum of an exponential and an algebraic term N=100 N=1000 red = (zc) -j green=.2/(j+1) 1+2 α blue=data Toeplitz at Brandeis April
48 terms: (z / z c ) j j =0 explains overall exponential decay, and period approximately 4, z = e iq e γ (log N )/N c a few low powers of z give some N-independent University of Chicago terms MRSEC, for small j explains N-independent bumps for small j term in integrand like (t-1) 2α gives additive term ψ j 1/ j 2α+1 explains behavior of eigenvalue... Toeplitz at Brandeis April
49 Explain behavior of eigenvalue Near j=n-1 end two terms balance out and permit system to reach end point. Exponential decay term ψ N 1 z c N = e γ(logn )/N N γ logn = e = N γ balances against algebraic decay ψ N 1 N (2α+1) giving γ = 2α +1 p = p 0 + i(2α +1)(logN) / N + O(1/ N) and an eigenvalue of the form and a p-value of the form ε = (2 e ip 0 e ip 0 ) α ( e ip 0 ) β [1 i logn N (2α +1)(β + α cot(p 0 / 2))] This is the first determination of the eigenvalue shift for finite N Toeplitz at Brandeis April
50 β=-1/2 α=1/4 theory gives 1+2α=1.5 observation gives offset=exp(-1.16)=.22 β=-1/2 α=1/3 theory gives 1+2α=1.67 observation gives offset=exp(-1.7) =.18 Toeplitz at Brandeis April
51 Conclusions I: What do we have? We have an accurate method of calculating eigenfunctions which works except at the low end of the eigenfunction (j near N-1 for negative β). The method works for both sign of β since the interchange of β and -β simply trades which region of j is the high side. no exponential term for positive β We have a determination of eigenvalue good to (log N)/N. These ideas apparently work across the entire eigenvalue spectrum but have not been fully evaluated near the ends of the spectrum. We do NOT have an understanding of the Hartwig-Fisher result for the determinant. Maybe the result comes from the ends of the spectrum. We expect more progress soon. Toeplitz at Brandeis April
52 j ε (-infinity, infinity) Conclusions II: Eigenvalues j ε(0, infinity) 1<β<0 j ε(0, infinity) 0<β<1 j ε{0, N-1} Toeplitz at Brandeis April
53 Conclusions III: information transfer j ε (-infinity, infinity) j ε(0, infinity) 1<β<0 K must transfer information in both directions information mostly transferred from zero to infinity j ε(0, infinity) 0<β<1 information mostly transferred from infinity to zero j ε{0, N-1} K must transfer information in both directions Toeplitz at Brandeis April
54 References Hui Dai, Zachary Geary, and Leo Kadanoff this work, arxiv: (2009) Toeplitz at Brandeis April
Page 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 informationUniform individual asymptotics for the eigenvalues and eigenvectors of large Toeplitz matrices
Uniform individual asymptotics for the eigenvalues and eigenvectors of large Toeplitz matrices Sergei Grudsky CINVESTAV, Mexico City, Mexico The International Workshop WIENER-HOPF METHOD, TOEPLITZ OPERATORS,
More information22.3. Repeated Eigenvalues and Symmetric Matrices. Introduction. Prerequisites. Learning Outcomes
Repeated Eigenvalues and Symmetric Matrices. Introduction In this Section we further develop the theory of eigenvalues and eigenvectors in two distinct directions. Firstly we look at matrices where one
More informationChapter Two Elements of Linear Algebra
Chapter Two Elements of Linear Algebra Previously, in chapter one, we have considered single first order differential equations involving a single unknown function. In the next chapter we will begin to
More informationCHAPTER 3. Matrix Eigenvalue Problems
A SERIES OF CLASS NOTES FOR 2005-2006 TO INTRODUCE LINEAR AND NONLINEAR PROBLEMS TO ENGINEERS, SCIENTISTS, AND APPLIED MATHEMATICIANS DE CLASS NOTES 3 A COLLECTION OF HANDOUTS ON SYSTEMS OF ORDINARY DIFFERENTIAL
More informationA = 3 B = A 1 1 matrix is the same as a number or scalar, 3 = [3].
Appendix : A Very Brief Linear ALgebra Review Introduction Linear Algebra, also known as matrix theory, is an important element of all branches of mathematics Very often in this course we study the shapes
More informationA VERY BRIEF LINEAR ALGEBRA REVIEW for MAP 5485 Introduction to Mathematical Biophysics Fall 2010
A VERY BRIEF LINEAR ALGEBRA REVIEW for MAP 5485 Introduction to Mathematical Biophysics Fall 00 Introduction Linear Algebra, also known as matrix theory, is an important element of all branches of mathematics
More informationQuantum Physics II (8.05) Fall 2004 Assignment 3
Quantum Physics II (8.5) Fall 24 Assignment 3 Massachusetts Institute of Technology Physics Department Due September 3, 24 September 23, 24 7:pm This week we continue to study the basic principles of quantum
More informationRepeated Eigenvalues and Symmetric Matrices
Repeated Eigenvalues and Symmetric Matrices. Introduction In this Section we further develop the theory of eigenvalues and eigenvectors in two distinct directions. Firstly we look at matrices where one
More informationPhase Transitions and Critical Behavior:
II Phase Transitions and Critical Behavior: A. Phenomenology (ibid., Chapter 10) B. mean field theory (ibid., Chapter 11) C. Failure of MFT D. Phenomenology Again (ibid., Chapter 12) // Windsor Lectures
More informationMAGNETIC BLOCH FUNCTIONS AND VECTOR BUNDLES. TYPICAL DISPERSION LAWS AND THEIR QUANTUM NUMBERS
MAGNETIC BLOCH FUNCTIONS AND VECTOR BUNDLES. TYPICAL DISPERSION LAWS AND THEIR QUANTUM NUMBERS S. P. NOVIKOV I. In previous joint papers by the author and B. A. Dubrovin [1], [2] we computed completely
More informationAPPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.
APPENDIX A Background Mathematics A. Linear Algebra A.. Vector algebra Let x denote the n-dimensional column vector with components 0 x x 2 B C @. A x n Definition 6 (scalar product). The scalar product
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 informationFinal Exam Practice Problems Answers Math 24 Winter 2012
Final Exam Practice Problems Answers Math 4 Winter 0 () The Jordan product of two n n matrices is defined as A B = (AB + BA), where the products inside the parentheses are standard matrix product. Is the
More informationPhysics 202 Laboratory 5. Linear Algebra 1. Laboratory 5. Physics 202 Laboratory
Physics 202 Laboratory 5 Linear Algebra Laboratory 5 Physics 202 Laboratory We close our whirlwind tour of numerical methods by advertising some elements of (numerical) linear algebra. There are three
More informationHilbert Spaces. Hilbert space is a vector space with some extra structure. We start with formal (axiomatic) definition of a vector space.
Hilbert Spaces Hilbert space is a vector space with some extra structure. We start with formal (axiomatic) definition of a vector space. Vector Space. Vector space, ν, over the field of complex numbers,
More informationImage Registration Lecture 2: Vectors and Matrices
Image Registration Lecture 2: Vectors and Matrices Prof. Charlene Tsai Lecture Overview Vectors Matrices Basics Orthogonal matrices Singular Value Decomposition (SVD) 2 1 Preliminary Comments Some of this
More informationIDENTIFICATION AND ANALYSIS OF TIME-VARYING MODAL PARAMETERS
IDENTIFICATION AND ANALYSIS OF TIME-VARYING MODAL PARAMETERS By STEPHEN L. SORLEY A THESIS PRESENTED TO THE GRADUATE SCHOOL OF THE UNIVERSITY OF FLORIDA IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE
More informationMathematical Methods for Engineers and Scientists 1
K.T. Tang Mathematical Methods for Engineers and Scientists 1 Complex Analysis, Determinants and Matrices With 49 Figures and 2 Tables fyj Springer Part I Complex Analysis 1 Complex Numbers 3 1.1 Our Number
More informationUnitary Dynamics and Quantum Circuits
qitd323 Unitary Dynamics and Quantum Circuits Robert B. Griffiths Version of 20 January 2014 Contents 1 Unitary Dynamics 1 1.1 Time development operator T.................................... 1 1.2 Particular
More informationPROBLEMS In each of Problems 1 through 12:
6.5 Impulse Functions 33 which is the formal solution of the given problem. It is also possible to write y in the form 0, t < 5, y = 5 e (t 5/ sin 5 (t 5, t 5. ( The graph of Eq. ( is shown in Figure 6.5.3.
More informationMOMENTS OF HYPERGEOMETRIC HURWITZ ZETA FUNCTIONS
MOMENTS OF HYPERGEOMETRIC HURWITZ ZETA FUNCTIONS ABDUL HASSEN AND HIEU D. NGUYEN Abstract. This paper investigates a generalization the classical Hurwitz zeta function. It is shown that many of the properties
More informationProperties of Linear Transformations from R n to R m
Properties of Linear Transformations from R n to R m MATH 322, Linear Algebra I J. Robert Buchanan Department of Mathematics Spring 2015 Topic Overview Relationship between the properties of a matrix transformation
More informationA matrix is a rectangular array of. objects arranged in rows and columns. The objects are called the entries. is called the size of the matrix, and
Section 5.5. Matrices and Vectors A matrix is a rectangular array of objects arranged in rows and columns. The objects are called the entries. A matrix with m rows and n columns is called an m n matrix.
More informationQuantum Computing Lecture 2. Review of Linear Algebra
Quantum Computing Lecture 2 Review of Linear Algebra Maris Ozols Linear algebra States of a quantum system form a vector space and their transformations are described by linear operators Vector spaces
More informationA matrix is a rectangular array of. objects arranged in rows and columns. The objects are called the entries. is called the size of the matrix, and
Section 5.5. Matrices and Vectors A matrix is a rectangular array of objects arranged in rows and columns. The objects are called the entries. A matrix with m rows and n columns is called an m n matrix.
More informationIntroduction to Instantons. T. Daniel Brennan. Quantum Mechanics. Quantum Field Theory. Effects of Instanton- Matter Interactions.
February 18, 2015 1 2 3 Instantons in Path Integral Formulation of mechanics is based around the propagator: x f e iht / x i In path integral formulation of quantum mechanics we relate the propagator to
More informationLinear Algebra using Dirac Notation: Pt. 2
Linear Algebra using Dirac Notation: Pt. 2 PHYS 476Q - Southern Illinois University February 6, 2018 PHYS 476Q - Southern Illinois University Linear Algebra using Dirac Notation: Pt. 2 February 6, 2018
More informationLinear Algebra- Final Exam Review
Linear Algebra- Final Exam Review. Let A be invertible. Show that, if v, v, v 3 are linearly independent vectors, so are Av, Av, Av 3. NOTE: It should be clear from your answer that you know the definition.
More information2. As we shall see, we choose to write in terms of σ x because ( X ) 2 = σ 2 x.
Section 5.1 Simple One-Dimensional Problems: The Free Particle Page 9 The Free Particle Gaussian Wave Packets The Gaussian wave packet initial state is one of the few states for which both the { x } and
More informationMatrix Representation
Matrix Representation Matrix Rep. Same basics as introduced already. Convenient method of working with vectors. Superposition Complete set of vectors can be used to express any other vector. Complete set
More informationChapter 5. Linear Algebra. A linear (algebraic) equation in. unknowns, x 1, x 2,..., x n, is. an equation of the form
Chapter 5. Linear Algebra A linear (algebraic) equation in n unknowns, x 1, x 2,..., x n, is an equation of the form a 1 x 1 + a 2 x 2 + + a n x n = b where a 1, a 2,..., a n and b are real numbers. 1
More informationDECAY OF SINGLET CONVERSION PROBABILITY IN ONE DIMENSIONAL QUANTUM NETWORKS
DECAY OF SINGLET CONVERSION PROBABILITY IN ONE DIMENSIONAL QUANTUM NETWORKS SCOTT HOTTOVY Abstract. Quantum networks are used to transmit and process information by using the phenomena of quantum mechanics.
More informationMITOCW ocw f99-lec23_300k
MITOCW ocw-18.06-f99-lec23_300k -- and lift-off on differential equations. So, this section is about how to solve a system of first order, first derivative, constant coefficient linear equations. And if
More informationVector Spaces. Vector space, ν, over the field of complex numbers, C, is a set of elements a, b,..., satisfying the following axioms.
Vector Spaces Vector space, ν, over the field of complex numbers, C, is a set of elements a, b,..., satisfying the following axioms. For each two vectors a, b ν there exists a summation procedure: a +
More informationEigenvalues and Eigenvectors
Contents Eigenvalues and Eigenvectors. Basic Concepts. Applications of Eigenvalues and Eigenvectors 8.3 Repeated Eigenvalues and Symmetric Matrices 3.4 Numerical Determination of Eigenvalues and Eigenvectors
More informationEigenvalues and Eigenvectors
5 Eigenvalues and Eigenvectors 5.2 THE CHARACTERISTIC EQUATION DETERMINANATS n n Let A be an matrix, let U be any echelon form obtained from A by row replacements and row interchanges (without scaling),
More informationIntroduction to Mathematical Physics
Introduction to Mathematical Physics Methods and Concepts Second Edition Chun Wa Wong Department of Physics and Astronomy University of California Los Angeles OXFORD UNIVERSITY PRESS Contents 1 Vectors
More informationVector spaces and operators
Vector spaces and operators Sourendu Gupta TIFR, Mumbai, India Quantum Mechanics 1 2013 22 August, 2013 1 Outline 2 Setting up 3 Exploring 4 Keywords and References Quantum states are vectors We saw that
More informationMTH Linear Algebra. Study Guide. Dr. Tony Yee Department of Mathematics and Information Technology The Hong Kong Institute of Education
MTH 3 Linear Algebra Study Guide Dr. Tony Yee Department of Mathematics and Information Technology The Hong Kong Institute of Education June 3, ii Contents Table of Contents iii Matrix Algebra. Real Life
More informationBackground Mathematics (2/2) 1. David Barber
Background Mathematics (2/2) 1 David Barber University College London Modified by Samson Cheung (sccheung@ieee.org) 1 These slides accompany the book Bayesian Reasoning and Machine Learning. The book and
More informationarxiv: v1 [cond-mat.stat-mech] 3 Jun 2018
Ising Models with Holes: Crossover Behavior arxiv:1806.00873v1 [cond-mat.stat-mech] 3 Jun 2018 Helen Au-Yang and Jacques H.H. Perk Department of Physics, Oklahoma State University, 145 Physical Sciences,
More informationDiscrete Time Fourier Transform (DTFT)
Discrete Time Fourier Transform (DTFT) 1 Discrete Time Fourier Transform (DTFT) The DTFT is the Fourier transform of choice for analyzing infinite-length signals and systems Useful for conceptual, pencil-and-paper
More informationAn Inverse Mass Expansion for Entanglement Entropy. Free Massive Scalar Field Theory
in Free Massive Scalar Field Theory NCSR Demokritos National Technical University of Athens based on arxiv:1711.02618 [hep-th] in collaboration with Dimitris Katsinis March 28 2018 Entanglement and Entanglement
More informationarxiv:math/ v1 [math.fa] 1 Jan 2001
arxiv:math/0101008v1 [math.fa] 1 Jan 2001 On the determinant formulas by Borodin, Okounkov, Baik, Deift, and Rains A. Böttcher We give alternative proofs to (block case versions of) some formulas for Toeplitz
More information-state problems and an application to the free particle
-state problems and an application to the free particle Sourendu Gupta TIFR, Mumbai, India Quantum Mechanics 1 2013 3 September, 2013 Outline 1 Outline 2 The Hilbert space 3 A free particle 4 Keywords
More information8.04: Quantum Mechanics Professor Allan Adams. Problem Set 7. Due Tuesday April 9 at 11.00AM
8.04: Quantum Mechanics Professor Allan Adams Massachusetts Institute of Technology Thursday April 4 Problem Set 7 Due Tuesday April 9 at 11.00AM Assigned Reading: E&R 6 all Li. 7 1 9, 8 1 Ga. 4 all, 6
More informationReview of Linear Algebra
Review of Linear Algebra Definitions An m n (read "m by n") matrix, is a rectangular array of entries, where m is the number of rows and n the number of columns. 2 Definitions (Con t) A is square if m=
More informationEigenvalues and Eigenvectors
LECTURE 3 Eigenvalues and Eigenvectors Definition 3.. Let A be an n n matrix. The eigenvalue-eigenvector problem for A is the problem of finding numbers λ and vectors v R 3 such that Av = λv. If λ, v are
More informationCITS2401 Computer Analysis & Visualisation
FACULTY OF ENGINEERING, COMPUTING AND MATHEMATICS CITS2401 Computer Analysis & Visualisation SCHOOL OF COMPUTER SCIENCE AND SOFTWARE ENGINEERING Topic 7 Matrix Algebra Material from MATLAB for Engineers,
More informationOutline 1. Real and complex p orbitals (and for any l > 0 orbital) 2. Dirac Notation :Symbolic vs shorthand Hilbert Space Vectors,
chmy564-19 Fri 18jan19 Outline 1. Real and complex p orbitals (and for any l > 0 orbital) 2. Dirac Notation :Symbolic vs shorthand Hilbert Space Vectors, 3. Theorems vs. Postulates Scalar (inner) prod.
More informationREVIEW REVIEW. Quantum Field Theory II
Quantum Field Theory II PHYS-P 622 Radovan Dermisek, Indiana University Notes based on: M. Srednicki, Quantum Field Theory Chapters: 13, 14, 16-21, 26-28, 51, 52, 61-68, 44, 53, 69-74, 30-32, 84-86, 75,
More informationQuantum Field Theory II
Quantum Field Theory II PHYS-P 622 Radovan Dermisek, Indiana University Notes based on: M. Srednicki, Quantum Field Theory Chapters: 13, 14, 16-21, 26-28, 51, 52, 61-68, 44, 53, 69-74, 30-32, 84-86, 75,
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 informationReview of scalar field theory. Srednicki 5, 9, 10
Review of scalar field theory Srednicki 5, 9, 10 2 The LSZ reduction formula based on S-5 In order to describe scattering experiments we need to construct appropriate initial and final states and calculate
More informationLecture Notes in Linear Algebra
Lecture Notes in Linear Algebra Dr. Abdullah Al-Azemi Mathematics Department Kuwait University February 4, 2017 Contents 1 Linear Equations and Matrices 1 1.2 Matrices............................................
More informationLinear Algebra. Matrices Operations. Consider, for example, a system of equations such as x + 2y z + 4w = 0, 3x 4y + 2z 6w = 0, x 3y 2z + w = 0.
Matrices Operations Linear Algebra Consider, for example, a system of equations such as x + 2y z + 4w = 0, 3x 4y + 2z 6w = 0, x 3y 2z + w = 0 The rectangular array 1 2 1 4 3 4 2 6 1 3 2 1 in which the
More informationChem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Fall Semester 2006 Christopher J. Cramer. Lecture 5, January 27, 2006
Chem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Fall Semester 2006 Christopher J. Cramer Lecture 5, January 27, 2006 Solved Homework (Homework for grading is also due today) We are told
More informationLecture 6. Eigen-analysis
Lecture 6 Eigen-analysis University of British Columbia, Vancouver Yue-Xian Li March 7 6 Definition of eigenvectors and eigenvalues Def: Any n n matrix A defines a LT, A : R n R n A vector v and a scalar
More informationMITOCW watch?v=y6ma-zn4olk
MITOCW watch?v=y6ma-zn4olk PROFESSOR: We have to ask what happens here? This series for h of u doesn't seem to stop. You go a 0, a 2, a 4. Well, it could go on forever. And what would happen if it goes
More informationConsistent Histories. Chapter Chain Operators and Weights
Chapter 10 Consistent Histories 10.1 Chain Operators and Weights The previous chapter showed how the Born rule can be used to assign probabilities to a sample space of histories based upon an initial state
More informationChapter 8. Rigid transformations
Chapter 8. Rigid transformations We are about to start drawing figures in 3D. There are no built-in routines for this purpose in PostScript, and we shall have to start more or less from scratch in extending
More informationSection 7.2: One-to-One, Onto and Inverse Functions
Section 7.2: One-to-One, Onto and Inverse Functions In this section we shall developed the elementary notions of one-to-one, onto and inverse functions, similar to that developed in a basic algebra course.
More informationTOEPLITZ OPERATORS. Toeplitz studied infinite matrices with NW-SE diagonals constant. f e C :
TOEPLITZ OPERATORS EFTON PARK 1. Introduction to Toeplitz Operators Otto Toeplitz lived from 1881-1940 in Goettingen, and it was pretty rough there, so he eventually went to Palestine and eventually contracted
More informationRandom matrices and the Riemann zeros
Random matrices and the Riemann zeros Bruce Bartlett Talk given in Postgraduate Seminar on 5th March 2009 Abstract Random matrices and the Riemann zeta function came together during a chance meeting between
More informationQueens College, CUNY, Department of Computer Science Numerical Methods CSCI 361 / 761 Spring 2018 Instructor: Dr. Sateesh Mane.
Queens College, CUNY, Department of Computer Science Numerical Methods CSCI 361 / 761 Spring 2018 Instructor: Dr. Sateesh Mane c Sateesh R. Mane 2018 8 Lecture 8 8.1 Matrices July 22, 2018 We shall study
More informationQuantum Mechanics for Mathematicians: Energy, Momentum, and the Quantum Free Particle
Quantum Mechanics for Mathematicians: Energy, Momentum, and the Quantum Free Particle Peter Woit Department of Mathematics, Columbia University woit@math.columbia.edu November 28, 2012 We ll now turn to
More informationA group G is a set of discrete elements a, b, x alongwith a group operator 1, which we will denote by, with the following properties:
1 Why Should We Study Group Theory? Group theory can be developed, and was developed, as an abstract mathematical topic. However, we are not mathematicians. We plan to use group theory only as much as
More informationPhysics 342 Lecture 26. Angular Momentum. Lecture 26. Physics 342 Quantum Mechanics I
Physics 342 Lecture 26 Angular Momentum Lecture 26 Physics 342 Quantum Mechanics I Friday, April 2nd, 2010 We know how to obtain the energy of Hydrogen using the Hamiltonian operator but given a particular
More informationMATH 583A REVIEW SESSION #1
MATH 583A REVIEW SESSION #1 BOJAN DURICKOVIC 1. Vector Spaces Very quick review of the basic linear algebra concepts (see any linear algebra textbook): (finite dimensional) vector space (or linear space),
More informationEigenvalues and Eigenvectors
5 Eigenvalues and Eigenvectors 5.2 THE CHARACTERISTIC EQUATION DETERMINANATS nn Let A be an matrix, let U be any echelon form obtained from A by row replacements and row interchanges (without scaling),
More informationGeometric Aspects of Sturm-Liouville Problems I. Structures on Spaces of Boundary Conditions
Geometric Aspects of Sturm-Liouville Problems I. Structures on Spaces of Boundary Conditions QINGKAI KONG HONGYOU WU and ANTON ZETTL Abstract. We consider some geometric aspects of regular Sturm-Liouville
More informationChapter 2. Linear Algebra. rather simple and learning them will eventually allow us to explain the strange results of
Chapter 2 Linear Algebra In this chapter, we study the formal structure that provides the background for quantum mechanics. The basic ideas of the mathematical machinery, linear algebra, are rather simple
More informationReal-Space Renormalisation
Real-Space Renormalisation LMU München June 23, 2009 Real-Space Renormalisation of the One-Dimensional Ising Model Ernst Ising Real-Space Renormalisation of the One-Dimensional Ising Model Ernst Ising
More informationSolving Homogeneous Systems with Sub-matrices
Pure Mathematical Sciences, Vol 7, 218, no 1, 11-18 HIKARI Ltd, wwwm-hikaricom https://doiorg/112988/pms218843 Solving Homogeneous Systems with Sub-matrices Massoud Malek Mathematics, California State
More informationEigenvectors and Hermitian Operators
7 71 Eigenvalues and Eigenvectors Basic Definitions Let L be a linear operator on some given vector space V A scalar λ and a nonzero vector v are referred to, respectively, as an eigenvalue and corresponding
More informationMathematical Introduction
Chapter 1 Mathematical Introduction HW #1: 164, 165, 166, 181, 182, 183, 1811, 1812, 114 11 Linear Vector Spaces: Basics 111 Field A collection F of elements a,b etc (also called numbers or scalars) with
More informationPhysics 218 Quantum Mechanics I Assignment 6
Physics 218 Quantum Mechanics I Assignment 6 Logan A. Morrison February 17, 2016 Problem 1 A non-relativistic beam of particles each with mass, m, and energy, E, which you can treat as a plane wave, is
More informationNPTEL
NPTEL Syllabus Selected Topics in Mathematical Physics - Video course COURSE OUTLINE Analytic functions of a complex variable. Calculus of residues, Linear response; dispersion relations. Analytic continuation
More information8 Square matrices continued: Determinants
8 Square matrices continued: Determinants 8.1 Introduction Determinants give us important information about square matrices, and, as we ll soon see, are essential for the computation of eigenvalues. You
More informationE = φ 1 A The dynamics of a particle with mass m and charge q is determined by the Hamiltonian
Lecture 9 Relevant sections in text: 2.6 Charged particle in an electromagnetic field We now turn to another extremely important example of quantum dynamics. Let us describe a non-relativistic particle
More informationChapter 3. Determinants and Eigenvalues
Chapter 3. Determinants and Eigenvalues 3.1. Determinants With each square matrix we can associate a real number called the determinant of the matrix. Determinants have important applications to the theory
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 informationChapter 2 Ensemble Theory in Statistical Physics: Free Energy Potential
Chapter Ensemble Theory in Statistical Physics: Free Energy Potential Abstract In this chapter, we discuss the basic formalism of statistical physics Also, we consider in detail the concept of the free
More informationQuantum Information Types
qitd181 Quantum Information Types Robert B. Griffiths Version of 6 February 2012 References: R. B. Griffiths, Types of Quantum Information, Phys. Rev. A 76 (2007) 062320; arxiv:0707.3752 Contents 1 Introduction
More informationOn the quantum theory of rotating electrons
Zur Quantentheorie des rotierenden Elektrons Zeit. f. Phys. 8 (98) 85-867. On the quantum theory of rotating electrons By Friedrich Möglich in Berlin-Lichterfelde. (Received on April 98.) Translated by
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 informationChecking Consistency. Chapter Introduction Support of a Consistent Family
Chapter 11 Checking Consistency 11.1 Introduction The conditions which define a consistent family of histories were stated in Ch. 10. The sample space must consist of a collection of mutually orthogonal
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 information4. Linear transformations as a vector space 17
4 Linear transformations as a vector space 17 d) 1 2 0 0 1 2 0 0 1 0 0 0 1 2 3 4 32 Let a linear transformation in R 2 be the reflection in the line = x 2 Find its matrix 33 For each linear transformation
More informationPseudospectra and Nonnormal Dynamical Systems
Pseudospectra and Nonnormal Dynamical Systems Mark Embree and Russell Carden Computational and Applied Mathematics Rice University Houston, Texas ELGERSBURG MARCH 1 Overview of the Course These lectures
More 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 informationEigenvalues in Applications
Eigenvalues in Applications Abstract We look at the role of eigenvalues and eigenvectors in various applications. Specifically, we consider differential equations, Markov chains, population growth, and
More informationDS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra.
DS-GA 1002 Lecture notes 0 Fall 2016 Linear Algebra These notes provide a review of basic concepts in linear algebra. 1 Vector spaces You are no doubt familiar with vectors in R 2 or R 3, i.e. [ ] 1.1
More informationElementary Linear Algebra
Matrices J MUSCAT Elementary Linear Algebra Matrices Definition Dr J Muscat 2002 A matrix is a rectangular array of numbers, arranged in rows and columns a a 2 a 3 a n a 2 a 22 a 23 a 2n A = a m a mn We
More information1 Algebra of State Vectors
J. Rothberg October 6, Introduction to Quantum Mechanics: Part Algebra of State Vectors What does the State Vector mean? A state vector is not a property of a physical system, but rather represents an
More information1 Expansion in orthogonal functions
Notes "orthogonal" January 9 Expansion in orthogonal functions To obtain a more useful form of the Green s function, we ll want to expand in orthogonal functions that are (relatively) easy to integrate.
More informationCh 10.1: Two Point Boundary Value Problems
Ch 10.1: Two Point Boundary Value Problems In many important physical problems there are two or more independent variables, so the corresponding mathematical models involve partial differential equations.
More informationSection 3: Complex numbers
Essentially: Section 3: Complex numbers C (set of complex numbers) up to different notation: the same as R 2 (euclidean plane), (i) Write the real 1 instead of the first elementary unit vector e 1 = (1,
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More information