Perturbation Theory , B =
|
|
- Shawn Giles Young
- 5 years ago
- Views:
Transcription
1 Appendix D Perturbation Theory D. Simple Examples Let A = 3, B = 4 6. (D.) The eigenvalues of A are, and 3, where λ = has multiplicity. To find the eigenvalue of the perturbed matrix A + ǫb corresponding to the unperturbed eigenvalue λ =, we use non-degenerate perturbation theory, so λ = λ + ǫλ + O ( ǫ ) (D.) = + ǫ x T Bx x T x + O ( ǫ ) (D.3) = + ǫ [ ] 4 6 [ ] + O ( ǫ ) (D.4) = + 4ǫ + O ( ǫ ), (D.5) noting that A = A T, so x, the eigenvector of AT with eigenvalue, is the same as x, the eigenvector of A with eigenvalue. We can do likewise for λ = 3 to find λ = 3 + 6ǫ + O ( ǫ ). (D.6) 3
2 For λ =, we must use degenerate perturbation theory, which says that Bc = λ Mc (D.7) where so equation D.7 becomes [ B ij v T i Bv j M ij v T i v j, (D.8) ] c = λ [ hence λ is i or i, so the perturbed eigenvalues for λ = are ] c, (D.9) λ = + ǫ { i i + O ( ǫ ) (D.) In summary, the eigenvalues of A + ǫb are ± iǫ λ = + 4ǫ + O ( ǫ ) (D.) 3 + 6ǫ These are compared to numerically calculated eigenvalues in figure D.. If instead, B = 4, (D.) 6 then perturbation theory says that for ǫ small, the eigenvalues of A + ǫb are ± ǫ λ = + 4ǫ + O ( ǫ ) (D.3) 3 + 6ǫ These are compared to numerically calculated eigenvalues in figure D.. D. Wilkinson s Matrix We now choose to look at an example where perturbation theory is very poor at approximating the perturbed eigenvalues. Following Wilkinson [3], we let 9 8 A = (D.4) 4
3 Real Component Real Component of First Eigenvalue Imag. Component Imag. Component of st Eigenvalue...5. Real Component of Second Eigenvalue. Real Component Imag. Component of nd Eigenvalue. Imag. Component..5. Real Component Real Component of Third Eigenvalue Imag. Component of 3rd Eigenvalue. Imag. Component..5. Real Component of Fourth Eigenvalue 3. Real Component Imag. Component of 4th Eigenvalue. Imag. Component..5. calculated by MATLAB predicted by perturbation theory Figure D.: The real (first column) and imaginary (second column) components of the four eigenvalues of A + ǫb, plotted versus ǫ, with A and B defined by equation D.. 5
4 Eigenvalue Eigenvalue Eigenvalue Eigenvalue. First Eigenvalue Second Eigenvalue Third Eigenvalue Fourth Eigenvalue.5. calculated by MATLAB predicted by perturbation theory Figure D.: The four eigenvalues (all real) of A+ǫB, plotted versus ǫ, where B is defined by equation D.. 6
5 This has eigenvalues,,..., and the eigenvector corresponding to eigenvalue λ is ( ) λ ( λ)!. ( )! (D.5) }. λ zeroes The eigenvector of A T corresponding to eigenvalue λ is D.. Let Eigenvalues. ()!.. () λ (λ )! } λ zeroes B =..... (D.6) (D.7) Then perturbation theory says that the eigenvalues of A+ǫB are given, to order ǫ, by λ n = λ n + ǫλ n = n + ǫ x T Bx x T x 9 ( ) n = n + ǫ ( n)!(n )! (D.8) (D.9) (D.) or, if we generalise A and B to be m by m matrices, λ n = n + ǫ mm ( ) m n (m n)!(n )! (D.) 7
6 Real Component of Eigenvalue Real Component of Tenth Eigenvalue 5 5 predicted range of validity of perturbation theory calculated by MATLAB predicted by perturbation theory Figure D.3: The th eigenvalue in the m = case, as a function of ǫ, as computed numerically by MATLAB. The eigenvalue predicted by perturbation theory is also shown. The predicted range gives the number η for which we must have ǫ η in order for perturbation theory to be valid, as in equation D.. where this is valid for m m ǫ. (D.) (m n)!(n )! For example, when m = 5, perturbation theory predicts λ 8 = 8 + ǫ (D.3) so we expect this to be valid for ǫ 8. The eigenvalues in the m = case, as a function of ǫ, were computed numerically in MATLAB and are shown in figures D.4 D.9. Also, the th eigenvalues in the m = case, as a function of ǫ, is compared to that predicted by perturbation theory. This is shown in figure D.3. Error in Computation Since the Wilkinson matrix has such poorly behaved eigenvalues, we expect that numerical computations of its perturbed eigenvalues will be difficult. To illustrate this, we compute the eigenvalues of the perturbed Wilkinson matrix for different values of m and ǫ using two different programs, MATLAB and Mathematica. The calculations were performed to 5 digits. The result is shown 8
7 Real Component of Eigenvalues Real Component of Eigenvalues Imaginary Component of Eigenvalues Eigenvalues 5 5 Real Component of Eigenvalues Imaginary Component of Eigenvalues Imaginary Component of Eigenvalues predicted range of validity of perturbation theory calculated by MATLAB Figure D.4: The eigenvalues in the m = case, as a function of ǫ, as computed numerically by MATLAB. The predicted range gives the number η, as a function of λ, for which we must have ǫ η in order for perturbation theory to be valid, as in equation D.. 9
8 Eigenvalues Imaginary Component 5 Real Component 3 Figure D.5: A 3D representation of the eigenvalues in the m = case, as a function of ǫ, as computed by MATLAB.
9 Size of Matrix Difference Between Mathematica and MATLAB Eigenvalues Region Where Eigenvalues Become All Real e+ e e 4 e 6 e 8 e e e 4 Figure D.6: The difference between the computed eigenvalues of A + ǫb according to MATLAB and Mathematica, as a function of m and ǫ. The dashed lines indicate the region in which the eigenvalues, according to Mathematica, become entirely real.
10 Size of Matrix e+ e e 4 e 6 e 8 e 5 e e Difference Between Mathematica and Perturbation Theory Eigenvalues Region Where Eigenvalues Become All Real Figure D.7: The difference between the computed eigenvalues of A + ǫb according to Mathematica and the predicted eigenvalues from perturbation theory, as a function of m and ǫ. The dashed lines indicate the region in which the eigenvalues, according to Mathematica, become entirely real. This also corresponds to the approximate boundary of the region of validity of perturbation theory.
11 Size of Matrix e+ e e 4 e 6 e 8 e 5 e e Difference Between MATLAB and Perturbation Theory Eigenvalues Region Where Eigenvalues Become All Real Figure D.8: The difference between the computed eigenvalues of A + ǫb according to MATLAB and the predicted eigenvalues from perturbation theory, as a function of m and ǫ. The dashed lines indicate the region in which the eigenvalues, according to Mathematica, become entirely real. This also corresponds to the approximate boundary of the region of validity of perturbation theory. 3
12 in figure D.6. More precisely, the plot is of the maximum of the absolute value of the complex difference between the eigenvalues as computed by MATLAB and Mathematica. To determine which is most accurate, we also plot each as compared with perturbation theory (figures D.7 and D.8). Observe that, within the region of validity of perturbation theory, Mathematica agrees very well with perturbation theory. We assume then that the error depicted in figure D.7 is due to errors from perturbation theory, not Mathematica. D.. Sherman-Woodbury-Morrison Formula If we wish to solve the related problem Ax + ǫbx = b, with A and B as above, then perturbation theory predicts x = ( A A ǫba ) b. (D.4) Since for our choice of B, we can write B = uv T, then we can apply the Sherman- Woodbury-Morrison formula which says that the exact solution is x = (A A ǫba ) + ǫv T A b, (D.5) u so we expect perturbation theory to be valid for ǫv T A u. Now 9!! ( ) 8!! ( ) 7! 9!! ( )! 8! 9! ( ) 7! 9! A 7! = 8!... so v T A u = ( )9! for m =. We can generalise this for any m to obtain!! (D.6) (D.7) We then expect perturbation theory to be valid for v T A u = ( m)m. (D.8) m! ǫ ( v T A u ) (D.9) ( ) ( m) m ǫ. m! (D.3) For example, when m = 5, this means that perturbation theory is valid for ǫ 9. 4
13 Size of Perturbation Imaginary Component 5 pseudospectrum eigenvalues of perturbed matrix Real Component 3 Figure D.9: A 3D version of the pseudospectrum for the unperturbed matrix (the surface), in the m = case, in comparison with the eigenvalues of the perturbed matrix (the lines), as computed by MATLAB. The lines in this figure are identical to those in figure D.5. D..3 Pseudospectrum Following Embree and Trefethen [7], we define the pseudospectrum of A to be Λ ǫ (A) = {z C : z Λ(A + E) for some E with E ǫ} (D.3) where Λ(A + E) is the set of eigenvalues of the matrix (A + E) and is a matrix norm induced by a vector norm. We choose to be the -norm. Since ǫb = ǫ, then the eigenvalues of the perturbed matrix A + ǫb are elements of Λ ǫ (A). The pseudospectrum of A in the m = case is shown in figures D.9 and D.. Note that the pseudospectrum of A agrees very well with the eigenvalues of A + ǫb as computed by MATLAB, although the eigenvalues of A + ǫb are all strictly within the boundary predicted by the pseudospectrum. This is especially visible in figure D.. This suggests that there exist some matrices C with C = such that the eigenvalues of A + ǫc are slightly further away from the unperturbed eigenvalues than the eigenvalues of A + ǫb. 5
14 5 4 Imaginary Component Real Component eigenvalues of unperturbed matrix eigenvalues of perturbed matrix contour plot of pseudospectrum 3 Figure D.: The pseudospectrum for the unperturbed matrix, in the m = case. 6
15 D.3 Domain Perturbations Consider now the first order corrections to the eigenvalues of a differential equation where the domain has been perturbed by some perturbation of order ǫ. We look in specific at the problem φ + λφ = in Ω ǫ (D.3a) φ = on δω ǫ (D.3b) where and δω ǫ : R = + ǫf (θ) (D.33) f (θ) = n= a n e inθ, a n = a n. (D.34) We will now find the first order corrections to the eigenvalues by two different methods, and show that they produce the same result. D.3. General Method For the first method, we start by writing φ = φ + ǫφ + ǫ φ +... λ = λ + ǫλ + ǫ λ +... (D.35) (D.36) and equating terms of equal order from equation D.3a, we obtain that O () : φ + λ φ = (D.37) O (ǫ) : φ + λ φ = λ φ (D.38) L(φ) = λ φ. (D.39) Performing a Taylor expansion of the boundary conditions, = φ( + ǫf (θ),θ) (D.4) = φ(,θ) + ǫf (θ) φ,r (,θ) + ǫ f (θ) φ,rr (,θ) and equating terms of equal order, (D.4) O () : = φ (,θ) O (ǫ) : = φ (,θ) + f (θ) φ,r (,θ) (D.4) (D.43) 7
16 First Eigenvalue The unperturbed eigenfunction of the first eigenvalue is φ = J (j r) (D.44) where j mk is the kth positive zero of the Bessel function J m. Now we know that (φ, L(φ )) (φ, L(φ )) = (φ φ,r φ φ,r ) ds (D.45) where (f,g) D D fg dx. Using equations D.37, D.39, D.4 and D.43, equation D.45 becomes λ (φ,φ ) = f (θ)(φ,r ) ds so λ = = = j D D f (θ)(φ,r) ds D (φ ) dx ( an e inθ) (J (j ) j ) dθ θ=π θ= n= θ=π dθ r= θ= r= r (J (j r)) dθ = (J (j )) j θ=π ( a + ( θ= n= an e inθ + a n e inθ)) dθ π (J (j )) ( θ=π a + (Re (a n ) cos (nθ) + θ= n= Im (a n ) sin (nθ)) ) dθ (D.46) (D.47) (D.48) (D.49) (D.5) (D.5) = j a. (D.5) The first eigenvalue is then λ = j ǫj a + O ( ǫ ). (D.53) Second and Third Eigenvalues The second and third eigenvalues of the unperturbed system are the same, so we have φ = c v + c v (D.54) where c and c are arbitrary constants and v = J (j r) cos (θ), v = J (j r) sin (θ). (D.55) 8
17 As before, we can use (v, L(φ )) (φ, L(v )) = D (v φ,r φ v,r ) ds, (D.56) which, upon applying equations D.37, D.39, D.4 and D.43, becomes λ (v,φ ) = f (θ) φ,r v,r ds (D.57) D λ c (v,v ) λ c (v,v ) = c f v,r ds + c f v,r v,r ds. (D.58) D We can do likewise, with v in equation D.56 replaced with v to obtain λ c (v,v ) λ c (v,v ) = c f v,r v,r ds + c f v,r ds. (D.59) D We can express these two equations together as [ ][ ] f v,r f v,r v,r c f v,r v,r f v,r c D D = λ [ (v,v ) (v,v ) (v,v ) (v,v ) ] [ c (D.6) Using the definition of v and v from equation D.55 and the result that r= r= r (J (j r)) dr = J (j ) J (j ) this simplifies to [ ] [ f v,r f v,r v,r c f v,r v,r f v,r c ] c ]. = (J (j )), (D.6) = π (J (j )) λ [ c c ]. (D.6) If we re-write this as [ A B B C ] [ c c ] [ c = Q c ] (D.63) then the eigenvalues are given by Q = A + C ± A + C AC + 4B, (D.64) 9
18 where A = = D θ=π θ= f (θ) (v,r (,θ)) ds (D.65) ( a n e )(J inθ (j ) j cos (θ)) dθ (D.66) = (J (j )) j = (J (j )) j n= θ=π θ= θ=π θ= ( ( a + a + ( an e inθ + a n e inθ)) cos (θ) dθ (D.67) n= (Re (a n ) cos (nθ) n= +Im (a n ) sin(nθ)) ) cos (θ) dθ (D.68) = (J (j )) jπ (Re(a ) + a ) (D.69) A + ξ. (D.7) We likewise obtain that C = f (θ) (v,r (,θ)) ds (D.7) D ( θ=π = a n e )(J inθ (j ) j sin(θ)) dθ (D.7) θ= n= = (J (j )) jπ ( Re (a ) + a ) (D.73) = A + ξ (D.74) and B = = D θ=π θ= f (θ)v,r (,θ) v,r (,θ) ds (D.75) ( a n e )(J inθ (j ) j ) sin(θ)cos (θ) dθ (D.76) n= = (J (j )) j π Im(a ), (D.77) so that Q = A + C ± A + C AC + 4B = ξ ± A + B. (D.78) (D.79)
19 From the definition of Q, we obtain that λ π (J (j )) = Q (D.8) π λ (J (j )) = (J (j )) ja π ± (J (j )) jπ a (D.8) λ = j (± a a ) (D.8) So the second eigenvalue is and the third is D.3. λ = j + ǫj ( a a ) + O ( ǫ ) λ = j + ǫj ( a a ) + O ( ǫ ). Method of Assumed Solution (D.83) (D.84) For this method, we follow Wolf and Keller [3]. We start by assuming that the solution is of the form φ(r,θ,ǫ) = A n (ǫ) J n (kr) e inθ (D.85) where n= A n (ǫ) = δ n m α n + ǫβ n + ǫ γ n +... R (θ,ǫ) = + ǫ n a n e inθ + ǫ n b n e inθ + O ( ǫ 3) (D.86) (D.87) λ k = k + ǫk + ǫ k +... Applying boundary conditions, φ(r (θ,ǫ),θ,ǫ) = A n (ǫ) J n (kr (θ,ǫ))e inθ =. n (D.88) (D.89) Expanding this in a Taylor series, and writing kr = k + (k k ) + k (R ), we obtain A n {J n (k ) + J n (k ) (k k + k (R )) + n } J n (k ) (k k + k (R )) +... e inθ =. (D.9) Applying equations D.86 D.88 and equating terms of equal order, we obtain that O () : ( α m e imθ + α m e imθ) J m (k ) = O (ǫ) : ( ( α m J me imθ + α m J me imθ) k + k a l e )+ ilθ l β n J n e inθ =. n (D.9) (D.9)
20 From equation D.9, we obtain that k = j mp, so m =, p = corresponds to the first eigenvalue and m =, p = will give the second and third eigenvalues. First Eigenvalue Under the condition that m =, equation D.9 becomes ( ) = α J k + k a l e ilθ + β n J n e inθ. n Equating terms that have θ =, we obtain l (D.93) = α J (j ) (k + k a ) + β J (j ) (D.94) = α J (j ) (k + k a ) (D.95) = k + k a, (D.96) so and k k = a (D.97) λ = k which agrees with equation D.53. Second and Third Eigenvalues (D.98) = k + ǫk k + O ( ǫ ) (D.99) = j ǫj a + O ( ǫ ), (D.) Again following Wolf and Keller [3], we equate the coefficients of e imθ in equation D.9 to obtain and = α m J m (k ) k a m + α m J m (k ) (k + k a ) + β m J m (k ) (D.) = α m J m (k ) k a m + α m J m (k ) (k + k a ) (D.) = α m k a m + α m (k + k a ) (D.3) k k = α ma m α m a. (D.4) Alternately, equating the coefficients of e imθ in equation D.9, we obtain Equating equations D.4 and D.5, we obtain that α m a m α m eiω = α m a m α m k k = α ma m α m a. (D.5) e iω (D.6)
21 where so that ( ) αm a m ω = arg α m ( ) αm a m = arg α m e iω = e iω e iω = ±. (D.7) (D.8) (D.9) (D.) It then follows that k = ± α m a m k α m a (D.) = ± a m a, (D.) so the second and third eigenvalues are λ = j + ǫj (± a a ) + O ( ǫ ) (D.3) which agrees with equations D.83 and D.84. D.3.3 Result Similarly, Wolf and Keller [3] compute the second order correction to the first eigenvalue. Then under the condition that the area of the domain does not change under the perturbation (which causes a =, among other things), they obtain that Aλ II = πj π/ j j [ Aλ I πj ] / + j J (j +O [ Aλ I πj ) /J (j ) ]. (D.4) where A is the area of the domain, and λ I and λ II are the first and second eigenvalues, respectively. Normalising A to unity, we obtain the plot shown in figure D.. 3
22 8 Second Eigenvalue (λ II ) First Eigenvalue (λ I ) Figure D.: The first and second eigenvalues for perturbations to a circle. When there are no perturbations, λ I is at a minimum. 4
A = 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 informationSpin Dynamics Basic Theory Operators. Richard Green SBD Research Group Department of Chemistry
Spin Dynamics Basic Theory Operators Richard Green SBD Research Group Department of Chemistry Objective of this session Introduce you to operators used in quantum mechanics Achieve this by looking at:
More informationMath 337, Summer 2010 Assignment 5
Math 337, Summer Assignment 5 Dr. T Hillen, University of Alberta Exercise.. Consider Laplace s equation r r r u + u r r θ = in a semi-circular disk of radius a centered at the origin with boundary conditions
More informationPerturbation Theory 1
Perturbation Theory 1 1 Expansion of Complete System Let s take a look of an expansion for the function in terms of the complete system : (1) In general, this expansion is possible for any complete set.
More informationPart IA. Vectors and Matrices. Year
Part IA Vectors and Matrices Year 2018 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2018 Paper 1, Section I 1C Vectors and Matrices For z, w C define the principal value of z w. State de Moivre s
More information( ) = 9φ 1, ( ) = 4φ 2.
Chemistry 46 Dr Jean M Standard Homework Problem Set 6 Solutions The Hermitian operator A ˆ is associated with the physical observable A Two of the eigenfunctions of A ˆ are and These eigenfunctions are
More informationFINAL EXAM MATH303 Theory of Ordinary Differential Equations. Spring dx dt = x + 3y dy dt = x y.
FINAL EXAM MATH0 Theory of Ordinary Differential Equations There are 5 problems on 2 pages. Spring 2009. 25 points Consider the linear plane autonomous system x + y x y. Find a fundamental matrix of the
More information) z r θ ( ) ( ) ( ) = then. Complete Solutions to Examination Questions Complete Solutions to Examination Questions 10.
Complete Solutions to Examination Questions 0 Complete Solutions to Examination Questions 0. (i We need to determine + given + j, j: + + j + j (ii The product ( ( + j6 + 6 j 8 + j is given by ( + j( j
More informationFinal: Solutions Math 118A, Fall 2013
Final: Solutions Math 118A, Fall 2013 1. [20 pts] For each of the following PDEs for u(x, y), give their order and say if they are nonlinear or linear. If they are linear, say if they are homogeneous or
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 informationMATH325 - QUANTUM MECHANICS - SOLUTION SHEET 11
MATH35 - QUANTUM MECHANICS - SOLUTION SHEET. The Hamiltonian for a particle of mass m moving in three dimensions under the influence of a three-dimensional harmonic oscillator potential is Ĥ = h m + mω
More informationExtra Problems for Math 2050 Linear Algebra I
Extra Problems for Math 5 Linear Algebra I Find the vector AB and illustrate with a picture if A = (,) and B = (,4) Find B, given A = (,4) and [ AB = A = (,4) and [ AB = 8 If possible, express x = 7 as
More informationAMS 212A Applied Mathematical Methods I Appendices of Lecture 06 Copyright by Hongyun Wang, UCSC. ( ) cos2
AMS 22A Applied Mathematical Methods I Appendices of Lecture 06 Copyright by Hongyun Wang UCSC Appendix A: Proof of Lemma Lemma : Let (x ) be the solution of x ( r( x)+ q( x) )sin 2 + ( a) 0 < cos2 where
More informationNonlinear Autonomous Systems of Differential
Chapter 4 Nonlinear Autonomous Systems of Differential Equations 4.0 The Phase Plane: Linear Systems 4.0.1 Introduction Consider a system of the form x = A(x), (4.0.1) where A is independent of t. Such
More informationTEST CODE: MIII (Objective type) 2010 SYLLABUS
TEST CODE: MIII (Objective type) 200 SYLLABUS Algebra Permutations and combinations. Binomial theorem. Theory of equations. Inequalities. Complex numbers and De Moivre s theorem. Elementary set theory.
More informationExamples and MatLab. Vector and Matrix Material. Matrix Addition R = A + B. Matrix Equality A = B. Matrix Multiplication R = A * B.
Vector and Matrix Material Examples and MatLab Matrix = Rectangular array of numbers, complex numbers or functions If r rows & c columns, r*c elements, r*c is order of matrix. A is n * m matrix Square
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 informationMTH4101 CALCULUS II REVISION NOTES. 1. COMPLEX NUMBERS (Thomas Appendix 7 + lecture notes) ax 2 + bx + c = 0. x = b ± b 2 4ac 2a. i = 1.
MTH4101 CALCULUS II REVISION NOTES 1. COMPLEX NUMBERS (Thomas Appendix 7 + lecture notes) 1.1 Introduction Types of numbers (natural, integers, rationals, reals) The need to solve quadratic equations:
More informationQuantum Physics III (8.06) Spring 2008 Final Exam Solutions
Quantum Physics III (8.6) Spring 8 Final Exam Solutions May 19, 8 1. Short answer questions (35 points) (a) ( points) α 4 mc (b) ( points) µ B B, where µ B = e m (c) (3 points) In the variational ansatz,
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science : Dynamic Systems Spring 2011
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science 6.24: Dynamic Systems Spring 20 Homework 9 Solutions Exercise 2. We can use additive perturbation model with
More informationPhysics 115C Homework 2
Physics 5C Homework Problem Our full Hamiltonian is H = p m + mω x +βx 4 = H +H where the unperturbed Hamiltonian is our usual and the perturbation is H = p m + mω x H = βx 4 Assuming β is small, the perturbation
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 informationREPRESENTATION THEORY WEEK 7
REPRESENTATION THEORY WEEK 7 1. Characters of L k and S n A character of an irreducible representation of L k is a polynomial function constant on every conjugacy class. Since the set of diagonalizable
More information(Again, this quantity is the correlation function of the two spins.) With z chosen along ˆn 1, this quantity is easily computed (exercise):
Lecture 30 Relevant sections in text: 3.9, 5.1 Bell s theorem (cont.) Assuming suitable hidden variables coupled with an assumption of locality to determine the spin observables with certainty we found
More informationMath Theory of Partial Differential Equations Lecture 3-2: Spectral properties of the Laplacian. Bessel functions.
Math 412-501 Theory of Partial ifferential Equations Lecture 3-2: Spectral properties of the Laplacian. Bessel functions. Eigenvalue problem: 2 φ + λφ = 0 in, ( αφ + β φ ) n = 0, where α, β are piecewise
More informationMath 489AB Exercises for Chapter 1 Fall Section 1.0
Math 489AB Exercises for Chapter 1 Fall 2008 Section 1.0 1.0.2 We want to maximize x T Ax subject to the condition x T x = 1. We use the method of Lagrange multipliers. Let f(x) = x T Ax and g(x) = x T
More informationMath 304 Answers to Selected Problems
Math Answers to Selected Problems Section 6.. Find the general solution to each of the following systems. a y y + y y y + y e y y y y y + y f y y + y y y + 6y y y + y Answer: a This is a system of the
More informationSecond quantization: where quantization and particles come from?
110 Phys460.nb 7 Second quantization: where quantization and particles come from? 7.1. Lagrangian mechanics and canonical quantization Q: How do we quantize a general system? 7.1.1.Lagrangian Lagrangian
More informationDifferentiation of functions of covariance
Differentiation of log X May 5, 2005 1 Differentiation of functions of covariance matrices or: Why you can forget you ever read this Richard Turner Covariance matrices are symmetric, but we often conveniently
More informationThis ODE arises in many physical systems that we shall investigate. + ( + 1)u = 0. (λ + s)x λ + s + ( + 1) a λ. (s + 1)(s + 2) a 0
Legendre equation This ODE arises in many physical systems that we shall investigate We choose We then have Substitution gives ( x 2 ) d 2 u du 2x 2 dx dx + ( + )u u x s a λ x λ a du dx λ a λ (λ + s)x
More informationA NOVEL METHOD TO DERIVE EXPLICIT KLT KERNEL FOR AR(1) PROCESS. Mustafa U. Torun and Ali N. Akansu
A NOVEL METHOD TO DERIVE EXPLICIT KLT KERNEL FOR AR() PROCESS Mustafa U. Torun and Ali N. Akansu Department of Electrical and Computer Engineering New Jersey Institute of Technology University Heights,
More informationMechanics Physics 151
Mechanics Phsics 151 Lecture 8 Rigid Bod Motion (Chapter 4) What We Did Last Time! Discussed scattering problem! Foundation for all experimental phsics! Defined and calculated cross sections! Differential
More informationTEST CODE: MMA (Objective type) 2015 SYLLABUS
TEST CODE: MMA (Objective type) 2015 SYLLABUS Analytical Reasoning Algebra Arithmetic, geometric and harmonic progression. Continued fractions. Elementary combinatorics: Permutations and combinations,
More information21 Laplace s Equation and Harmonic Functions
2 Laplace s Equation and Harmonic Functions 2. Introductory Remarks on the Laplacian operator Given a domain Ω R d, then 2 u = div(grad u) = in Ω () is Laplace s equation defined in Ω. If d = 2, in cartesian
More information4.3 Classical model at low temperature: an example Continuum limit and path integral The two-point function: perturbative
Contents Gaussian integrals.......................... Generating function.......................2 Gaussian expectation values. Wick s theorem........... 2.3 Perturbed gaussian measure. Connected contributions.......
More information3 + 4i 2 + 3i. 3 4i Fig 1b
The introduction of complex numbers in the 16th century was a natural step in a sequence of extensions of the positive integers, starting with the introduction of negative numbers (to solve equations of
More informationReview of Fundamental Equations Supplementary notes on Section 1.2 and 1.3
Review of Fundamental Equations Supplementary notes on Section. and.3 Introduction of the velocity potential: irrotational motion: ω = u = identity in the vector analysis: ϕ u = ϕ Basic conservation principles:
More informationSpiral Wave Chimeras
Spiral Wave Chimeras Carlo R. Laing IIMS, Massey University, Auckland, New Zealand Erik Martens MPI of Dynamics and Self-Organization, Göttingen, Germany Steve Strogatz Dept of Maths, Cornell University,
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra for Computational and Data Sciences)
AMS526: Numerical Analysis I (Numerical Linear Algebra for Computational and Data Sciences) Lecture 1: Course Overview; Matrix Multiplication Xiangmin Jiao Stony Brook University Xiangmin Jiao Numerical
More information24 Solving planar heat and wave equations in polar coordinates
24 Solving planar heat and wave equations in polar coordinates Now that all the preparations are done, I can return to solving the planar heat and wave equations in domains with rotational symmetry. 24.1
More informationInvertibility and stability. Irreducibly diagonally dominant. Invertibility and stability, stronger result. Reducible matrices
Geršgorin circles Lecture 8: Outline Chapter 6 + Appendix D: Location and perturbation of eigenvalues Some other results on perturbed eigenvalue problems Chapter 8: Nonnegative matrices Geršgorin s Thm:
More informationRAJASTHAN P.E.T. MATHS-1995
RAJASTHAN P.E.T. MATHS-1995 1. The equation of the normal to the circle x2 + y2 = a2 at point (x y ) will be : (1) x y - xy = (2) xx - yy = (3) x y + xy = (4) xx + yy = 2. Equation of the bisector of the
More informationFunctional Analysis Review
Outline 9.520: Statistical Learning Theory and Applications February 8, 2010 Outline 1 2 3 4 Vector Space Outline A vector space is a set V with binary operations +: V V V and : R V V such that for all
More informationPH.D. PRELIMINARY EXAMINATION MATHEMATICS
UNIVERSITY OF CALIFORNIA, BERKELEY Dept. of Civil and Environmental Engineering FALL SEMESTER 2014 Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem
More informationAPPENDIX 2.1 LINE AND SURFACE INTEGRALS
2 APPENDIX 2. LINE AND URFACE INTEGRAL Consider a path connecting points (a) and (b) as shown in Fig. A.2.. Assume that a vector field A(r) exists in the space in which the path is situated. Then the line
More informationSingular Value Decomposition (SVD)
School of Computing National University of Singapore CS CS524 Theoretical Foundations of Multimedia More Linear Algebra Singular Value Decomposition (SVD) The highpoint of linear algebra Gilbert Strang
More informationMath 108b: Notes on the Spectral Theorem
Math 108b: Notes on the Spectral Theorem From section 6.3, we know that every linear operator T on a finite dimensional inner product space V has an adjoint. (T is defined as the unique linear operator
More informationNOTES ON LINEAR ALGEBRA CLASS HANDOUT
NOTES ON LINEAR ALGEBRA CLASS HANDOUT ANTHONY S. MAIDA CONTENTS 1. Introduction 2 2. Basis Vectors 2 3. Linear Transformations 2 3.1. Example: Rotation Transformation 3 4. Matrix Multiplication and Function
More informationBOUNDARY PROBLEMS IN HIGHER DIMENSIONS. kt = X X = λ, and the series solutions have the form (for λ n 0):
BOUNDARY PROBLEMS IN HIGHER DIMENSIONS Time-space separation To solve the wave and diffusion equations u tt = c u or u t = k u in a bounded domain D with one of the three classical BCs (Dirichlet, Neumann,
More informationQualification Exam: Mathematical Methods
Qualification Exam: Mathematical Methods Name:, QEID#41534189: August, 218 Qualification Exam QEID#41534189 2 1 Mathematical Methods I Problem 1. ID:MM-1-2 Solve the differential equation dy + y = sin
More information11 a 12 a 21 a 11 a 22 a 12 a 21. (C.11) A = The determinant of a product of two matrices is given by AB = A B 1 1 = (C.13) and similarly.
C PROPERTIES OF MATRICES 697 to whether the permutation i 1 i 2 i N is even or odd, respectively Note that I =1 Thus, for a 2 2 matrix, the determinant takes the form A = a 11 a 12 = a a 21 a 11 a 22 a
More informationSpecial Symmetrical Mtrix 1
Special Symmetrical Mtrix Eigenvalues and eigenvectors for Special Symmetrical Matrix We shall consider the following special symmetrical matrix with n rows and n columns denoted y S: 0 0 0 0 S = 0 0 0
More information+ i. cos(t) + 2 sin(t) + c 2.
MATH HOMEWORK #7 PART A SOLUTIONS Problem 7.6.. Consider the system x = 5 x. a Express the general solution of the given system of equations in terms of realvalued functions. b Draw a direction field,
More information(8.51) ẋ = A(λ)x + F(x, λ), where λ lr, the matrix A(λ) and function F(x, λ) are C k -functions with k 1,
2.8.7. Poincaré-Andronov-Hopf Bifurcation. In the previous section, we have given a rather detailed method for determining the periodic orbits of a two dimensional system which is the perturbation of a
More informationChapter 6 Higher Dimensional Linear Systems
hapter 6 Higher Dimensional Linear Systems Linear systems in R n : Observations: X x. x n X X a a n..... a n a nn ; Let T the the coordinate change matrix such that T T is in canonical form ;then Y T X
More informationQuantum Physics II (8.05) Fall 2002 Assignment 11
Quantum Physics II (8.05) Fall 00 Assignment 11 Readings Most of the reading needed for this problem set was already given on Problem Set 9. The new readings are: Phase shifts are discussed in Cohen-Tannoudji
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 informationLAPLACE EQUATION. = 2 is call the Laplacian or del-square operator. In two dimensions, the expression of in rectangular and polar coordinates are
LAPLACE EQUATION If a diffusion or wave problem is stationary (time independent), the pde reduces to the Laplace equation u = u =, an archetype of second order elliptic pde. = 2 is call the Laplacian or
More informationQuasi-Hermitian operators
Nuclear Physics Institute, Řež Faculty of Nuclear Sciences and Physical Engineering, Prague Stará Lesná 2008 Outline 1 Quasi-Hermiticity dim H = N dim H = 2 Pseudo-Hermiticity Antilinear symmetry J-self-adjointness
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 information3/31/ Product of Inertia. Sample Problem Sample Problem 10.6 (continue)
/1/01 10.6 Product of Inertia Product of Inertia: I xy = xy da When the x axis, the y axis, or both are an axis of symmetry, the product of inertia is zero. Parallel axis theorem for products of inertia:
More informationforms Christopher Engström November 14, 2014 MAA704: Matrix factorization and canonical forms Matrix properties Matrix factorization Canonical forms
Christopher Engström November 14, 2014 Hermitian LU QR echelon Contents of todays lecture Some interesting / useful / important of matrices Hermitian LU QR echelon Rewriting a as a product of several matrices.
More informationx x2 2 + x3 3 x4 3. Use the divided-difference method to find a polynomial of least degree that fits the values shown: (b)
Numerical Methods - PROBLEMS. The Taylor series, about the origin, for log( + x) is x x2 2 + x3 3 x4 4 + Find an upper bound on the magnitude of the truncation error on the interval x.5 when log( + x)
More information(x 1, y 1 ) = (x 2, y 2 ) if and only if x 1 = x 2 and y 1 = y 2.
1. Complex numbers A complex number z is defined as an ordered pair z = (x, y), where x and y are a pair of real numbers. In usual notation, we write z = x + iy, where i is a symbol. The operations of
More information18.03 The Heat Equation
18.03 The Heat Equation Jeremy Orloff 1 Introduction The goal in this note is to derive the heat equation from Newton s law of cooling. We can get most of the way using just mathematics. The final piece
More informationDavid Eigen. MA112 Final Paper. May 10, 2002
David Eigen MA112 Fina Paper May 1, 22 The Schrodinger equation describes the position of an eectron as a wave. The wave function Ψ(t, x is interpreted as a probabiity density for the position of the eectron.
More informationSolutions for MAS277 Problems
Solutions for MAS77 Problems Solutions for Chapter problems: Inner product spaces 1. No. The only part that fails, however, is the condition that f, f should mean that f. If f is a non-zero function that
More informationProblem Set 8 Mar 5, 2004 Due Mar 10, 2004 ACM 95b/100b 3pm at Firestone 303 E. Sterl Phinney (2 pts) Include grading section number
Problem Set 8 Mar 5, 24 Due Mar 1, 24 ACM 95b/1b 3pm at Firestone 33 E. Sterl Phinney (2 pts) Include grading section number Useful Readings: For Green s functions, see class notes and refs on PS7 (esp
More informationReview Questions REVIEW QUESTIONS 71
REVIEW QUESTIONS 71 MATLAB, is [42]. For a comprehensive treatment of error analysis and perturbation theory for linear systems and many other problems in linear algebra, see [126, 241]. An overview of
More informationTime-Independent Perturbation Theory
4 Phys46.nb Time-Independent Perturbation Theory.. Overview... General question Assuming that we have a Hamiltonian, H = H + λ H (.) where λ is a very small real number. The eigenstates of the Hamiltonian
More informationMATH 220: INNER PRODUCT SPACES, SYMMETRIC OPERATORS, ORTHOGONALITY
MATH 22: INNER PRODUCT SPACES, SYMMETRIC OPERATORS, ORTHOGONALITY When discussing separation of variables, we noted that at the last step we need to express the inhomogeneous initial or boundary data as
More informationConformal maps. Lent 2019 COMPLEX METHODS G. Taylor. A star means optional and not necessarily harder.
Lent 29 COMPLEX METHODS G. Taylor A star means optional and not necessarily harder. Conformal maps. (i) Let f(z) = az + b, with ad bc. Where in C is f conformal? cz + d (ii) Let f(z) = z +. What are the
More informationBefore you begin read these instructions carefully:
NATURAL SCIENCES TRIPOS Part IB & II (General Friday, 30 May, 2014 9:00 am to 12:00 pm MATHEMATICS (2 Before you begin read these instructions carefully: You may submit answers to no more than six questions.
More informationENGI Linear Approximation (2) Page Linear Approximation to a System of Non-Linear ODEs (2)
ENGI 940 4.06 - Linear Approximation () Page 4. 4.06 Linear Approximation to a System of Non-Linear ODEs () From sections 4.0 and 4.0, the non-linear system dx dy = x = P( x, y), = y = Q( x, y) () with
More informationMie Theory, Cross-Sections, Efficiency Factors and Phase Function
Appendix C Mie Theory, Cross-Sections, Efficiency Factors and Phase Function C.1 Mie Theory The concept behind the Mie Theory is explained in section.6.1. However, the mathematical derivation of the theory
More informationMaxwell Equations Waves Helmholtz Wave Equation Dirac Delta Sources Power Matrices Complex Analysis Optimization. EM & Math - Basics
EM & Math - S. R. Zinka srinivasa_zinka@daiict.ac.in October 16, 2014 Outline 1 Maxwell Equations 2 Waves 3 Helmholtz Wave Equation 4 Dirac Delta Sources 5 Power 6 Matrices 7 Complex Analysis 8 Optimization
More informationChapter 14 Three Ways of Treating a Linear Delay Differential Equation
Chapter 14 Three Ways of Treating a Linear Delay Differential Equation Si Mohamed Sah and Richard H. Rand Abstract This work concerns the occurrence of Hopf bifurcations in delay differential equations
More informationDepartment of Mathematics
INDIAN INSTITUTE OF TECHNOLOGY, BOMBAY Department of Mathematics MA 04 - Complex Analysis & PDE s Solutions to Tutorial No.13 Q. 1 (T) Assuming that term-wise differentiation is permissible, show that
More informationMath 46, Applied Math (Spring 2008): Final
Math 46, Applied Math (Spring 2008): Final 3 hours, 80 points total, 9 questions, roughly in syllabus order (apart from short answers) 1. [16 points. Note part c, worth 7 points, is independent of the
More information2017 YEAR 5 PROMOTION EXAMINATION MATHEMATICS 9758
RAFFLES INSTITUTION 07 YEAR 5 PROMOTION EXAMINATION MATHEMATICS 9758 September/October 07 Total Marks: 00 3 hours Additional materials: Answer Paper List of Formulae (MF6) READ THESE INSTRUCTIONS FIRST
More informationMATA Y, Tutorial #11, Lucas Ashbury-Bridgwood
MATA 018Y, Tutorial #11, Lucas Ashbury-Bridgwood 1 today Q5 pick-up + solution A11 Q5 pick-up 0:100:15 3 Q5 solution 0:150:35 A11.1 A11 #3e 0:300:50 0 Find the eigenvalues A and compute det A, det A 3
More information8.05, Quantum Physics II, Fall 2013 TEST Wednesday October 23, 12:30-2:00pm You have 90 minutes.
8.05, Quantum Physics II, Fall 03 TEST Wednesday October 3, :30-:00pm You have 90 minutes. Answer all problems in the white books provided. Write YOUR NAME and YOUR SECTION on your white books). There
More informationThe Higgins-Selkov oscillator
The Higgins-Selkov oscillator May 14, 2014 Here I analyse the long-time behaviour of the Higgins-Selkov oscillator. The system is ẋ = k 0 k 1 xy 2, (1 ẏ = k 1 xy 2 k 2 y. (2 The unknowns x and y, being
More informationPartial Differential Equations
M3M3 Partial Differential Equations Solutions to problem sheet 3/4 1* (i) Show that the second order linear differential operators L and M, defined in some domain Ω R n, and given by Mφ = Lφ = j=1 j=1
More informationQuantum dynamics with non-hermitian PT -symmetric operators: Models
Hauptseminar Theoretische Physik 01 Quantum dynamics with non-hermitian PT -symmetric operators: Models Mario Schwartz 13.06.01 Mario Schwartz PT -Symmetric Operators: Models 1 / 36 Overview Hauptseminar
More informationAdvanced Computational Fluid Dynamics AA215A Lecture 2 Approximation Theory. Antony Jameson
Advanced Computational Fluid Dynamics AA5A Lecture Approximation Theory Antony Jameson Winter Quarter, 6, Stanford, CA Last revised on January 7, 6 Contents Approximation Theory. Least Squares Approximation
More informationSymmetry and Properties of Crystals (MSE638) Stress and Strain Tensor
Symmetry and Properties of Crystals (MSE638) Stress and Strain Tensor Somnath Bhowmick Materials Science and Engineering, IIT Kanpur April 6, 2018 Tensile test and Hooke s Law Upto certain strain (0.75),
More informationEXERCISES ON DETERMINANTS, EIGENVALUES AND EIGENVECTORS. 1. Determinants
EXERCISES ON DETERMINANTS, EIGENVALUES AND EIGENVECTORS. Determinants Ex... Let A = 0 4 4 2 0 and B = 0 3 0. (a) Compute 0 0 0 0 A. (b) Compute det(2a 2 B), det(4a + B), det(2(a 3 B 2 )). 0 t Ex..2. For
More informationBTL What is the value of m if the vector is solenoidal. BTL What is the value of a, b, c if the vector may be irrotational.
VALLIAMMAI ENGINEERING OLLEGE SRM NAGAR, KATTANDKULATHUR Department of Mathematics MA65 - MATHEMATIS II QUESTION BANK - 6 UNIT - I VETOR ALULUS Part - A. Find, if at (, -, ). BTL-. Find the Directional
More informationK th -order linear difference equation (Normal Form)
The standard form of the general K th -order difference equation is x n + 1 = F(x n, x n - 1, x n - 2,..., x n - K + 1, P) where P is a vector of parameters. However, linear equations are normally expressed
More informationMath October 20, 2006 Exam 2: Solutions
Math 412-51 October 2, 26 Exam 2: Solutions Problem 1 (5 pts) Solve the heat equation in a rectangle < x < π, < y < π, subject to the initial condition and the boundary conditions u t = 2 u x + 2 u 2 y
More informationPoint estimates for Green s matrix to boundary value problems for second order elliptic systems in a polyhedral cone
Maz ya, V. G., Roßmann, J.: Estimates for Green s matrix 1 ZAMM Z. angew. Math. Mech. 00 2004 0, 1 30 Maz ya, V. G.; Roßmann, J. Point estimates for Green s matrix to boundary value problems for second
More informationPlatzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Introduction to PDEs and Numerical Methods Lecture 6: Numerical solution of the heat equation with FD method: method of lines, Euler
More informationMAT01A1: Complex Numbers (Appendix H)
MAT01A1: Complex Numbers (Appendix H) Dr Craig 14 February 2018 Announcements: e-quiz 1 is live. Deadline is Wed 21 Feb at 23h59. e-quiz 2 (App. A, D, E, H) opens tonight at 19h00. Deadline is Thu 22 Feb
More informationSolutions for Chapter 3
Solutions for Chapter Solutions for exercises in section 0 a X b x, y 6, and z 0 a Neither b Sew symmetric c Symmetric d Neither The zero matrix trivially satisfies all conditions, and it is the only possible
More informationSOUTH AFRICAN TERTIARY MATHEMATICS OLYMPIAD
SOUTH AFRICAN TERTIARY MATHEMATICS OLYMPIAD. Determine the following value: 7 August 6 Solutions π + π. Solution: Since π
More informationDegenerate Perturbation Theory
Physics G6037 Professor Christ 12/05/2014 Degenerate Perturbation Theory The treatment of degenerate perturbation theory presented in class is written out here in detail. 1 General framework and strategy
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 informationMath 46, Applied Math (Spring 2009): Final
Math 46, Applied Math (Spring 2009): Final 3 hours, 80 points total, 9 questions worth varying numbers of points 1. [8 points] Find an approximate solution to the following initial-value problem which
More information