Variational Principles for Nonlinear Eigenvalue Problems

Size: px
Start display at page:

Download "Variational Principles for Nonlinear Eigenvalue Problems"

Transcription

1 Variational Principles for Nonlinear Eigenvalue Problems Heinrich Voss Hamburg University of Technology Institute of Mathematics TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

2 Nonlinear eigenvalue problem The term nonlinear eigenvalue problem is not used in a unique way in the literature TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

3 Nonlinear eigenvalue problem The term nonlinear eigenvalue problem is not used in a unique way in the literature On the one hand: Parameter dependent nonlinear (with respect to the state variable u) operator equations are discussed concerning positivity of solutions multiplicity of solution T (λ, u) = 0 dependence of solutions on the parameter; bifurcation (change of ) stability of solutions TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

4 In this presentation For λ D C let T (λ) be a linear self-adjoint and bounded operator on a Hilbert space H. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

5 In this presentation For λ D C let T (λ) be a linear self-adjoint and bounded operator on a Hilbert space H. Find λ D and x 0 such that T (λ)x = 0. (1) Then λ is called an eigenvalue of T ( ), and x a corresponding eigenelement. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

6 In this presentation For λ D C let T (λ) be a linear self-adjoint and bounded operator on a Hilbert space H. Find λ D and x 0 such that T (λ)x = 0. (1) Then λ is called an eigenvalue of T ( ), and x a corresponding eigenelement. Nonlinear eigenproblems of this type arise in dynamic/stability analysis of structures and in fluid mechanics TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

7 In this presentation For λ D C let T (λ) be a linear self-adjoint and bounded operator on a Hilbert space H. Find λ D and x 0 such that T (λ)x = 0. (1) Then λ is called an eigenvalue of T ( ), and x a corresponding eigenelement. Nonlinear eigenproblems of this type arise in dynamic/stability analysis of structures and in fluid mechanics electronic behavior of semiconductor hetero-structures TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

8 In this presentation For λ D C let T (λ) be a linear self-adjoint and bounded operator on a Hilbert space H. Find λ D and x 0 such that T (λ)x = 0. (1) Then λ is called an eigenvalue of T ( ), and x a corresponding eigenelement. Nonlinear eigenproblems of this type arise in dynamic/stability analysis of structures and in fluid mechanics electronic behavior of semiconductor hetero-structures vibration of fluid-solid structures TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

9 In this presentation For λ D C let T (λ) be a linear self-adjoint and bounded operator on a Hilbert space H. Find λ D and x 0 such that T (λ)x = 0. (1) Then λ is called an eigenvalue of T ( ), and x a corresponding eigenelement. Nonlinear eigenproblems of this type arise in dynamic/stability analysis of structures and in fluid mechanics electronic behavior of semiconductor hetero-structures vibration of fluid-solid structures vibration of sandwich plates TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

10 In this presentation For λ D C let T (λ) be a linear self-adjoint and bounded operator on a Hilbert space H. Find λ D and x 0 such that T (λ)x = 0. (1) Then λ is called an eigenvalue of T ( ), and x a corresponding eigenelement. Nonlinear eigenproblems of this type arise in dynamic/stability analysis of structures and in fluid mechanics electronic behavior of semiconductor hetero-structures vibration of fluid-solid structures vibration of sandwich plates accelerator design TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

11 In this presentation For λ D C let T (λ) be a linear self-adjoint and bounded operator on a Hilbert space H. Find λ D and x 0 such that T (λ)x = 0. (1) Then λ is called an eigenvalue of T ( ), and x a corresponding eigenelement. Nonlinear eigenproblems of this type arise in dynamic/stability analysis of structures and in fluid mechanics electronic behavior of semiconductor hetero-structures vibration of fluid-solid structures vibration of sandwich plates accelerator design vibro-acoustics of piezoelectric/poroelastic structures TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

12 In this presentation For λ D C let T (λ) be a linear self-adjoint and bounded operator on a Hilbert space H. Find λ D and x 0 such that T (λ)x = 0. (1) Then λ is called an eigenvalue of T ( ), and x a corresponding eigenelement. Nonlinear eigenproblems of this type arise in dynamic/stability analysis of structures and in fluid mechanics electronic behavior of semiconductor hetero-structures vibration of fluid-solid structures vibration of sandwich plates accelerator design vibro-acoustics of piezoelectric/poroelastic structures nonlinear integrated optics TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

13 In this presentation For λ D C let T (λ) be a linear self-adjoint and bounded operator on a Hilbert space H. Find λ D and x 0 such that T (λ)x = 0. (1) Then λ is called an eigenvalue of T ( ), and x a corresponding eigenelement. Nonlinear eigenproblems of this type arise in dynamic/stability analysis of structures and in fluid mechanics electronic behavior of semiconductor hetero-structures vibration of fluid-solid structures vibration of sandwich plates accelerator design vibro-acoustics of piezoelectric/poroelastic structures nonlinear integrated optics regularization of total least squares problems TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

14 In this presentation For λ D C let T (λ) be a linear self-adjoint and bounded operator on a Hilbert space H. Find λ D and x 0 such that T (λ)x = 0. (1) Then λ is called an eigenvalue of T ( ), and x a corresponding eigenelement. Nonlinear eigenproblems of this type arise in dynamic/stability analysis of structures and in fluid mechanics electronic behavior of semiconductor hetero-structures vibration of fluid-solid structures vibration of sandwich plates accelerator design vibro-acoustics of piezoelectric/poroelastic structures nonlinear integrated optics regularization of total least squares problems stability of delay differential equations TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

15 Outline 1 Two examples 2 Variational characterization of eigenvalue problems Overdamped Problems Nonoverdamped problems 3 An unsymmetric linear eigenproblem 4 Concluding remarks TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

16 Outline Two examples 1 Two examples 2 Variational characterization of eigenvalue problems Overdamped Problems Nonoverdamped problems 3 An unsymmetric linear eigenproblem 4 Concluding remarks TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

17 Two examples Example 1: Electronic behavior of quantum dots Semiconductor nanostructures have attracted tremendous interest in the past few years because of their special physical properties and their potential for applications in micro and optoelectronic devices. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

18 Two examples Example 1: Electronic behavior of quantum dots Semiconductor nanostructures have attracted tremendous interest in the past few years because of their special physical properties and their potential for applications in micro and optoelectronic devices. In such nanostructures, the free carriers are confined to a small region of space by potential barriers, and if the size of this region is less than the electron wavelength, the electronic states become quantized at discrete energy levels. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

19 Two examples Example 1: Electronic behavior of quantum dots Semiconductor nanostructures have attracted tremendous interest in the past few years because of their special physical properties and their potential for applications in micro and optoelectronic devices. In such nanostructures, the free carriers are confined to a small region of space by potential barriers, and if the size of this region is less than the electron wavelength, the electronic states become quantized at discrete energy levels. The ultimate limit of low dimensional structures is the quantum dot, in which the carriers are confined in all three directions, thus reducing the degrees of freedom to zero. Therefore, a quantum dot can be thought of as an artificial atom. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

20 Problem Two examples Determine relevant energy states (i.e. eigenvalues) and corresponding wave functions (i.e. eigenfunctions) of a three-dimensional quantum dot embedded in a matrix. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

21 Problem ct. Two examples Governing equation: Schrödinger equation ( 2 ) 2m(x, E) Φ + V (x)φ = EΦ, x Ω q Ω m where is the reduced Planck constant, m(x, E) is the electron effective mass, and V (x) is the confinement potential. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

22 Problem ct. Two examples Governing equation: Schrödinger equation ( 2 ) 2m(x, E) Φ + V (x)φ = EΦ, x Ω q Ω m where is the reduced Planck constant, m(x, E) is the electron effective mass, and V (x) is the confinement potential. m and V are discontinous across the heterojunction. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

23 Problem ct. Two examples Governing equation: Schrödinger equation ( 2 ) 2m(x, E) Φ + V (x)φ = EΦ, x Ω q Ω m where is the reduced Planck constant, m(x, E) is the electron effective mass, and V (x) is the confinement potential. m and V are discontinous across the heterojunction. Boundary and interface conditions Φ = 0 on outer boundary of matrix Ω m 1 Φ BenDaniel Duke condition m m n = 1 Φ Ωm m q n Ωq on interface TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

24 Variational form Two examples Find E R and Φ H 1 0 (Ω), Φ 0, Ω := Ω q Ω m, such that a(φ, Ψ; E) := Φ Ψ dx + m q (x, E) 2 Ω q Ω m 1 Φ Ψ dx m m (x, E) + V q (x)φψ dx + V m (x)φψ dx Ω q Ω m = E ΦΨ dx =: Eb(Φ, Ψ) for every Ψ H0 1 (Ω) Ω TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

25 Two examples Electron effective mass The dependence of m(x, E) on E can be derived from the eight-band k p analysis and effective mass theory. Projecting the 8 8 Hamiltonian onto the conduction band results in the single Hamiltonian eigenvalue problem with { mq (E), x Ω m(x, E) = q m m (E), x Ω m ( 1 m j (E) = P2 j E + g j V j 1 E + g j V j + δ j ), j {m, q} where m j is the electron effective mass, V j the confinement potential, P j the momentum, g j the main energy gap, and δ j the spin-orbit splitting in the jth region. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

26 Two examples Electron effective mass The dependence of m(x, E) on E can be derived from the eight-band k p analysis and effective mass theory. Projecting the 8 8 Hamiltonian onto the conduction band results in the single Hamiltonian eigenvalue problem with { mq (E), x Ω m(x, E) = q m m (E), x Ω m ( 1 m j (E) = P2 j E + g j V j 1 E + g j V j + δ j ), j {m, q} where m j is the electron effective mass, V j the confinement potential, P j the momentum, g j the main energy gap, and δ j the spin-orbit splitting in the jth region. Other types of effective mass (taking into account the effect of strain, e.g.) appear in the literature. They are all rational functions of E where 1/m(x, E) is monotonically decreasing with respect to E, and that s all we need. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

27 Properties Two examples a(,, E) bilinear, symmetric, bounded, H0 1 (Ω) elliptic for E 0 b(, ) bilinear, positive definite, bounded, completely continuous By the Lax Milgram lemma the variational eigenproblem is equivalent to T (E)Φ = 0 where T (E) : H0 1 (Ω) H0 1 (Ω), E 0, is a family of self-adjoint and bounded operators. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

28 Two examples Example 2: Free vibrations of fluid-solid structures can be modelled in terms of solid displacement and fluid pressure and one obtains the classical form of an eigenproblem div [σ(u)] + ω 2 ρ s u = 0 in Ω s, p + ω2 c 2 p = 0 in Ω f, σ(u) n pn = 0 on Γ I, p n + ω 2 ρ f u n = 0 on Γ I, u = 0 on Γ D, p n = 0 on Γ N, u: solid displacement p: fluid pressure λ = ω 2 : eigenparameter σ(u): linearized stress tensor ρ s, ρ f : densities of solid and fluid Γ I Ω f Ω s Interface conditions: equilibrium of accelerations and of force densities. Γ O TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

29 Two examples Variational and operator form Find λ := ω 2 C and (u, p) HΓ 1 D (Ω s ) 3 H 1 (Ω f ) such that a s (v, u) + c(v, p) = λb s (v, u) and a f (q, p) = λ( c(u, q) + b f (q, p)). for every (v, q) HΓ 1 D (Ω s ) 3 H 1 (Ω f ). TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

30 Two examples Variational and operator form Find λ := ω 2 C and (u, p) H 1 Γ D (Ω s ) 3 H 1 (Ω f ) such that a s (v, u) + c(v, p) = λb s (v, u) and for every (v, q) H 1 Γ D (Ω s ) 3 H 1 (Ω f ). a f (q, p) = λ( c(u, q) + b f (q, p)). which (using the Lax-Milgram Lemma) can be transformed into a linear (but not self-adjoint) eigenvalue problem K s u + Cp = λm s u K f p = λ( C u + M f p) (2a) (2b) where K s : H 1 Γ D (Ω s ) 3 H 1 Γ D (Ω s ) 3 is self-adjoint, elliptic, bounded,... TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

31 Two examples Rational form of fluid-solid eigenproblem Let 0 < σ 1 σ 2... denote the eigenvalues of the decoupled eigenproblem K s u = σm s u and denote by u 1, u 2,... corresponding orthonormal eigenfunctions. Then the spectral theorem yields (K s λm s ) 1 u = n=1 1 σ n λ u n, u u n. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

32 Two examples Rational form of fluid-solid eigenproblem Let 0 < σ 1 σ 2... denote the eigenvalues of the decoupled eigenproblem K s u = σm s u and denote by u 1, u 2,... corresponding orthonormal eigenfunctions. Then the spectral theorem yields (K s λm s ) 1 u = n=1 1 σ n λ u n, u u n. If λ is not contained in the spectrum of the decoupled solid eigenproblem, then λ is an eigenvalue of the coupled fluid-solid problem if and only if it is an eigenvalue of the rational eigenvalue problem T (λ)p := K f p + λm f p + n=1 λ σ n λ C np, C n p := u n, Cp C u n. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

33 Two examples Rational form of fluid-solid eigenproblem Let 0 < σ 1 σ 2... denote the eigenvalues of the decoupled eigenproblem K s u = σm s u and denote by u 1, u 2,... corresponding orthonormal eigenfunctions. Then the spectral theorem yields (K s λm s ) 1 u = n=1 1 σ n λ u n, u u n. If λ is not contained in the spectrum of the decoupled solid eigenproblem, then λ is an eigenvalue of the coupled fluid-solid problem if and only if it is an eigenvalue of the rational eigenvalue problem T (λ)p := K f p + λm f p + n=1 λ σ n λ C np, C n p := u n, Cp C u n. T (λ) : H 1 (Ω f ) H 1 (Ω f ) is self-adjoint and bounded. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

34 Two examples Quadratic form of fluid-solid eigenproblem Another self-adjoint form of the fluid-solid eigenproblem is obtained if the second equation in (2) is multiplied by ω and p is substituted by p =: ωw. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

35 Two examples Quadratic form of fluid-solid eigenproblem Another self-adjoint form of the fluid-solid eigenproblem is obtained if the second equation in (2) is multiplied by ω and p is substituted by p =: ωw. Then problem (2) is equivalent to the quadratic eigenvalue problem (( ) ( ) ( )) ( ) Ks O O C + ω O C ω 2 Ms O u = 0. O O M f w K f TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

36 Outline Variational characterization of eigenvalue problems 1 Two examples 2 Variational characterization of eigenvalue problems Overdamped Problems Nonoverdamped problems 3 An unsymmetric linear eigenproblem 4 Concluding remarks TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

37 Outline Variational characterization of eigenvalue problems Overdamped Problems 1 Two examples 2 Variational characterization of eigenvalue problems Overdamped Problems Nonoverdamped problems 3 An unsymmetric linear eigenproblem 4 Concluding remarks TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

38 Overdamped Problems Variational characterization for linear eigenproblems Let A : H H a bounded linear and self-adjoint operator in a Hilbert space H. Then those eigenvalues λ 1 λ 2... above the essential spectrum of A (if there are any) can be characterized by three fundamental variational principles, TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

39 Overdamped Problems Variational characterization for linear eigenproblems Let A : H H a bounded linear and self-adjoint operator in a Hilbert space H. Then those eigenvalues λ 1 λ 2... above the essential spectrum of A (if there are any) can be characterized by three fundamental variational principles,by Rayleigh s principle λ n = max{r(x) : x, x i = 0, i = 1,..., n 1}, R(x) := Ax, x x, x. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

40 Overdamped Problems Variational characterization for linear eigenproblems Let A : H H a bounded linear and self-adjoint operator in a Hilbert space H. Then those eigenvalues λ 1 λ 2... above the essential spectrum of A (if there are any) can be characterized by three fundamental variational principles,by Rayleigh s principle λ n = max{r(x) : x, x i = 0, i = 1,..., n 1}, R(x) := the maxmin characterization due to Poincaré λ n = max dim V =n min R(x). x V, x 0 Ax, x x, x. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

41 Overdamped Problems Variational characterization for linear eigenproblems Let A : H H a bounded linear and self-adjoint operator in a Hilbert space H. Then those eigenvalues λ 1 λ 2... above the essential spectrum of A (if there are any) can be characterized by three fundamental variational principles,by Rayleigh s principle λ n = max{r(x) : x, x i = 0, i = 1,..., n 1}, R(x) := the maxmin characterization due to Poincaré λ n = max dim V =n min R(x). x V, x 0 the minmax characterization due to Courant, Fischer, and Weyl. λ n = min dim V =n 1 max R(x). x V,x 0 Ax, x x, x. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

42 Variational characterizations are very powerful tools when studying self-adjoint linear operators on a Hilbert space H. Bounds for eigenvalues, comparison theorems, interlacing results and monotonicity of eigenvalues can be proved easily with these characterizations, to name just a few. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52 Variational characterization of eigenvalue problems Overdamped Problems Variational characterization for linear eigenproblems Let A : H H a bounded linear and self-adjoint operator in a Hilbert space H. Then those eigenvalues λ 1 λ 2... above the essential spectrum of A (if there are any) can be characterized by three fundamental variational principles,by Rayleigh s principle λ n = max{r(x) : x, x i = 0, i = 1,..., n 1}, R(x) := the maxmin characterization due to Poincaré λ n = max dim V =n min R(x). x V, x 0 the minmax characterization due to Courant, Fischer, and Weyl. λ n = min dim V =n 1 max R(x). x V,x 0 Ax, x x, x.

43 Rayleigh functional Overdamped Problems Let f : { J H R (λ, x) T (λ)x, x be continuous, and assume that for every fixed x H, x 0, the real equation has at most one solution in J. f (λ, x) = 0 (3) TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

44 Rayleigh functional Overdamped Problems Let f : { J H R (λ, x) T (λ)x, x be continuous, and assume that for every fixed x H, x 0, the real equation has at most one solution in J. f (λ, x) = 0 (3) Then equation (3) implicitly defines a functional p on some subset D of H \ {0} which we call the Rayleigh functional, and which is exactly the Rayleigh quotient in case of a linear eigenproblem T (λ) = λi A. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

45 Rayleigh functional Overdamped Problems Let f : { J H R (λ, x) T (λ)x, x be continuous, and assume that for every fixed x H, x 0, the real equation has at most one solution in J. f (λ, x) = 0 (3) Then equation (3) implicitly defines a functional p on some subset D of H \ {0} which we call the Rayleigh functional, and which is exactly the Rayleigh quotient in case of a linear eigenproblem T (λ) = λi A. We further assume that for every x D, x 0 and λ J, λ p(x) it holds that f (λ, x)(λ p(x)) > 0 which generalizes the definiteness of the operator B for the generalized linear eigenproblem T (λ) := λb A. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

46 Overdamped problems Overdamped Problems If the Rayleigh functional p is defined on the entire space H \ {0} then the eigenproblem T (λ)x = 0 is called overdamped. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

47 Overdamped problems Overdamped Problems If the Rayleigh functional p is defined on the entire space H \ {0} then the eigenproblem T (λ)x = 0 is called overdamped. This notation is motivated by the finite dimensional quadratic eigenvalue problem T (λ)x = λ 2 Mx + λαcx + Kx = 0 where M, C and K are symmetric and positive definite matrices. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

48 Overdamped problems Overdamped Problems If the Rayleigh functional p is defined on the entire space H \ {0} then the eigenproblem T (λ)x = 0 is called overdamped. This notation is motivated by the finite dimensional quadratic eigenvalue problem T (λ)x = λ 2 Mx + λαcx + Kx = 0 where M, C and K are symmetric and positive definite matrices. α = 0 all eigenvalues on imaginary axis TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

49 Overdamped problems Overdamped Problems If the Rayleigh functional p is defined on the entire space H \ {0} then the eigenproblem T (λ)x = 0 is called overdamped. This notation is motivated by the finite dimensional quadratic eigenvalue problem T (λ)x = λ 2 Mx + λαcx + Kx = 0 where M, C and K are symmetric and positive definite matrices. α = 0 increase α all eigenvalues on imaginary axis eigenvalues go into left half plane as conjugate complex pairs TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

50 Overdamped problems Overdamped Problems If the Rayleigh functional p is defined on the entire space H \ {0} then the eigenproblem T (λ)x = 0 is called overdamped. This notation is motivated by the finite dimensional quadratic eigenvalue problem T (λ)x = λ 2 Mx + λαcx + Kx = 0 where M, C and K are symmetric and positive definite matrices. α = 0 increase α increase α all eigenvalues on imaginary axis eigenvalues go into left half plane as conjugate complex pairs complex pairs reach real axis, run in opposite directions TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

51 Overdamped problems Overdamped Problems If the Rayleigh functional p is defined on the entire space H \ {0} then the eigenproblem T (λ)x = 0 is called overdamped. This notation is motivated by the finite dimensional quadratic eigenvalue problem T (λ)x = λ 2 Mx + λαcx + Kx = 0 where M, C and K are symmetric and positive definite matrices. α = 0 increase α increase α increase α all eigenvalues on imaginary axis eigenvalues go into left half plane as conjugate complex pairs complex pairs reach real axis, run in opposite directions all eigenvalues on the negative real axis TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

52 Overdamped problems Overdamped Problems If the Rayleigh functional p is defined on the entire space H \ {0} then the eigenproblem T (λ)x = 0 is called overdamped. This notation is motivated by the finite dimensional quadratic eigenvalue problem T (λ)x = λ 2 Mx + λαcx + Kx = 0 where M, C and K are symmetric and positive definite matrices. α = 0 increase α increase α increase α increase α all eigenvalues on imaginary axis eigenvalues go into left half plane as conjugate complex pairs complex pairs reach real axis, run in opposite directions all eigenvalues on the negative real axis all eigenvalues going to the left are smaller than all eigenvalues going to the right system is overdamped TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

53 Overdamped Problems Quadratic overdamped problems For quadratic overdamped systems the two solutions p ± (x) = 1 ( ) α Cx, x ± α 2 2 Cx, x 2 4 Mx, x Kx, x / Mx, x. of the quadratic equation T (λ)x, x = λ 2 Mx, x + λα Cx, x + Kx, x = 0 (4) are real, and they satisfy sup x 0 p (x) < inf x 0 p + (x). TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

54 Overdamped Problems Quadratic overdamped problems For quadratic overdamped systems the two solutions p ± (x) = 1 ( ) α Cx, x ± α 2 2 Cx, x 2 4 Mx, x Kx, x / Mx, x. of the quadratic equation T (λ)x, x = λ 2 Mx, x + λα Cx, x + Kx, x = 0 (4) are real, and they satisfy sup x 0 p (x) < inf x 0 p + (x). Hence, equation (4) defines two Rayleigh functionals p and p + corresponding to the intervals J := (, inf p +(x)) and J + := (sup p (x), ). x 0 x 0 TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

55 Rayleigh s principle Overdamped Problems For general (not necessarily quadratic) overdamped problems Hadeler (1967 for the finite dimensional case, and 1968 for dim H = ) generalized Rayleigh s principle proving that the eigenvectors are orthogonal with respect to the generalized scalar product T (p(x)) T (p(y)) x, y, if p(x) p(y) [x, y] := p(x) p(y) (5) T (p(x))x, y, if p(x) = p(y) which is symmetric, definite and homogeneous, but in general is not bilinear. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

56 Overdamped Problems Rayleigh s principle (Hadeler 1967,1968) Let T (λ) : H H, λ J be a family of linear self-adjoint and bounded operators such that (1) is over-damped, and assume that for λ J there exists ν(λ) > 0 such that T (λ) + ν(λ)i is compact. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

57 Overdamped Problems Rayleigh s principle (Hadeler 1967,1968) Let T (λ) : H H, λ J be a family of linear self-adjoint and bounded operators such that (1) is over-damped, and assume that for λ J there exists ν(λ) > 0 such that T (λ) + ν(λ)i is compact. Let T ( ) be continuously differentiable and suppose that T (λ)x, x λ > 0 for ever x 0. λ=p(x) TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

58 Overdamped Problems Rayleigh s principle (Hadeler 1967,1968) Let T (λ) : H H, λ J be a family of linear self-adjoint and bounded operators such that (1) is over-damped, and assume that for λ J there exists ν(λ) > 0 such that T (λ) + ν(λ)i is compact. Let T ( ) be continuously differentiable and suppose that T (λ)x, x λ > 0 for ever x 0. λ=p(x) Then problem T (λ)x = 0 has at most a countable set of eigenvalues in J which we assume to be ordered by magnitude λ 1 λ 2..., where each eigenvalue is counted according to its multiplicity. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

59 Overdamped Problems Rayleigh s principle (Hadeler 1967,1968) Let T (λ) : H H, λ J be a family of linear self-adjoint and bounded operators such that (1) is over-damped, and assume that for λ J there exists ν(λ) > 0 such that T (λ) + ν(λ)i is compact. Let T ( ) be continuously differentiable and suppose that T (λ)x, x λ > 0 for ever x 0. λ=p(x) Then problem T (λ)x = 0 has at most a countable set of eigenvalues in J which we assume to be ordered by magnitude λ 1 λ 2..., where each eigenvalue is counted according to its multiplicity. The corresponding eigenvectors x 1, x 2,... can be chosen orthonormally with respect to the generalized scalar product (5), and the eigenvalues can be determined recurrently by λ n = min{p(x) : [x, x i ] = 0, i = 1,..., n 1, x 0}. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

60 Minmax principle Overdamped Problems Poincaré s maxmin characterization was first generalized by Duffin (1955) to overdamped quadratic eigenproblems of finite dimension, and for more general overdamped problems of finite dimension it was proved by Rogers (1964). TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

61 Minmax principle Overdamped Problems Poincaré s maxmin characterization was first generalized by Duffin (1955) to overdamped quadratic eigenproblems of finite dimension, and for more general overdamped problems of finite dimension it was proved by Rogers (1964). Infinite dimensional eigenvalue problems were studied by Turner (1967), Langer (1968), and Weinberger (1969) who proved generalizations of both, the maxmin characterization of Poincaré and of the minmax characterization of Courant, Fischer and Weyl for quadratic (and by Turner (1968) for polynomial) overdamped problems. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

62 Minmax principle Overdamped Problems Poincaré s maxmin characterization was first generalized by Duffin (1955) to overdamped quadratic eigenproblems of finite dimension, and for more general overdamped problems of finite dimension it was proved by Rogers (1964). Infinite dimensional eigenvalue problems were studied by Turner (1967), Langer (1968), and Weinberger (1969) who proved generalizations of both, the maxmin characterization of Poincaré and of the minmax characterization of Courant, Fischer and Weyl for quadratic (and by Turner (1968) for polynomial) overdamped problems. The corresponding generalizations for general overdamped problems of infinite dimension were derived by Hadeler (1968). Similar results (weakening the compactness or smoothness requirements) are contained in Rogers (1968), Werner (1971), Abramov (1973), Hadeler (1975), Markus (1985), Maksudov & Gasanov (1992), and Hasanov (2002). TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

63 Overdamped Problems Minmax and Maxmin principle (Hadeler 1968) Let T (λ) : H H, λ J be a family of linear self-adjoint and bounded operators such that (1) is over-damped, and assume that for λ J there exists ν(λ) > 0 such that T (λ) + ν(λ)i is compact. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

64 Overdamped Problems Minmax and Maxmin principle (Hadeler 1968) Let T (λ) : H H, λ J be a family of linear self-adjoint and bounded operators such that (1) is over-damped, and assume that for λ J there exists ν(λ) > 0 such that T (λ) + ν(λ)i is compact. Let T ( ) be continuously differentiable and suppose that T (λ)x, x λ > 0 for ever x 0. λ=p(x) TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

65 Overdamped Problems Minmax and Maxmin principle (Hadeler 1968) Let T (λ) : H H, λ J be a family of linear self-adjoint and bounded operators such that (1) is over-damped, and assume that for λ J there exists ν(λ) > 0 such that T (λ) + ν(λ)i is compact. Let T ( ) be continuously differentiable and suppose that T (λ)x, x λ > 0 for ever x 0. λ=p(x) Let the eigenvalues λ n of T (λ)x = 0 be numbered in non-decreasing order regarding their multiplicities. Then they can be characterized by the following two variational principles λ n = min dim V =n = max dim V =n 1 max p(x) x V, x 0 min p(x). x V, x 0 TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

66 Overdamped Problems Example 1: Quantum dot problem For the quantum dot problem the family of operators T (λ) satisfies the general conditions of the variational characterizations. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

67 Overdamped Problems Example 1: Quantum dot problem For the quantum dot problem the family of operators T (λ) satisfies the general conditions of the variational characterizations. For x 0 it holds that a(x, x, 0) > 0, b(x, x) > 0 and λ a(x, x, λ) is monotonically decreasing for λ 0. Hence, f (λ, x) = T (λ)x, x = λb(x, x) a(x, x, λ) = 0 has exactly one positive solution p(x), and f (λ, x) λ > 0. λ=p(x) TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

68 Overdamped Problems Example 1: Quantum dot problem For the quantum dot problem the family of operators T (λ) satisfies the general conditions of the variational characterizations. For x 0 it holds that a(x, x, 0) > 0, b(x, x) > 0 and λ a(x, x, λ) is monotonically decreasing for λ 0. Hence, f (λ, x) = T (λ)x, x = λb(x, x) a(x, x, λ) = 0 has exactly one positive solution p(x), and f (λ, x) λ > 0. λ=p(x) Thus, the quantum dot problem has a countable number of non-negative eigenvalues which allow for all three variational characterizations. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

69 Overdamped Problems Example 2: Fluid-solid vibration For the rational eigenproblem governing free vibrations the family of operators T (λ) satisfies the general conditions of the variational characterizations in every interval J n := (σ n 1, σ n ). TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

70 Overdamped Problems Example 2: Fluid-solid vibration For the rational eigenproblem governing free vibrations the family of operators T (λ) satisfies the general conditions of the variational characterizations in every interval J n := (σ n 1, σ n ). f (λ, p) := T (λ)p, p = K f p, p + λ M f p, p + n=1 λ σ n λ u n, Cp 2 is monotonically increasing, such that f (λ, p) = 0 has at most one solution in J n, but the Rayleigh functional is not defined on the entire space. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

71 Outline Variational characterization of eigenvalue problems Nonoverdamped problems 1 Two examples 2 Variational characterization of eigenvalue problems Overdamped Problems Nonoverdamped problems 3 An unsymmetric linear eigenproblem 4 Concluding remarks TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

72 Non-overdamped problems Nonoverdamped problems For nonoverdamped eigenproblems (i.e. D(p) H \ {0}) the natural ordering to call the smallest eigenvalue the first one, the second smallest the second one, etc., is not appropriate. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

73 Non-overdamped problems Nonoverdamped problems For nonoverdamped eigenproblems (i.e. D(p) H \ {0}) the natural ordering to call the smallest eigenvalue the first one, the second smallest the second one, etc., is not appropriate. This is obvious if we make a linear eigenvalue T (λ)x := (λi A)x = 0 nonlinear by restricting it to an interval J which does not contain the smallest eigenvalue of A. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

74 Non-overdamped problems Nonoverdamped problems For nonoverdamped eigenproblems (i.e. D(p) H \ {0}) the natural ordering to call the smallest eigenvalue the first one, the second smallest the second one, etc., is not appropriate. This is obvious if we make a linear eigenvalue T (λ)x := (λi A)x = 0 nonlinear by restricting it to an interval J which does not contain the smallest eigenvalue of A. Then all conditions are satisfied, p is the restriction of the Rayleigh quotient R A to D(p) := {x 0 : R A (x) J}, and inf x D(p) p(x) will not be an eigenvalue. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

75 Enumeration of eigenvalues Nonoverdamped problems λ J is an eigenvalue of T ( ) if and only if µ = 0 is an eigenvalue of the linear problem T (λ)y = µy. The key idea is to orientate the number of λ on the location on the eigenvalue µ = 0 in the spectrum of the linear operator T (λ). TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

76 Enumeration of eigenvalues Nonoverdamped problems λ J is an eigenvalue of T ( ) if and only if µ = 0 is an eigenvalue of the linear problem T (λ)y = µy. The key idea is to orientate the number of λ on the location on the eigenvalue µ = 0 in the spectrum of the linear operator T (λ). We assume that for every λ J it holds that the supremum of the essential spectrum of T (λ) is negative ( for instance: there exists ν(λ) > 0 such that T (λ) + ν(λ)i is compact). TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

77 Enumeration of eigenvalues Nonoverdamped problems λ J is an eigenvalue of T ( ) if and only if µ = 0 is an eigenvalue of the linear problem T (λ)y = µy. The key idea is to orientate the number of λ on the location on the eigenvalue µ = 0 in the spectrum of the linear operator T (λ). We assume that for every λ J it holds that the supremum of the essential spectrum of T (λ) is negative ( for instance: there exists ν(λ) > 0 such that T (λ) + ν(λ)i is compact). If λ J is an eigenvalue of T ( ) then there exists n N such that 0 = max dimv =n min x V, x 0 T (λ)x, x. x, x TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

78 Enumeration of eigenvalues Nonoverdamped problems λ J is an eigenvalue of T ( ) if and only if µ = 0 is an eigenvalue of the linear problem T (λ)y = µy. The key idea is to orientate the number of λ on the location on the eigenvalue µ = 0 in the spectrum of the linear operator T (λ). We assume that for every λ J it holds that the supremum of the essential spectrum of T (λ) is negative ( for instance: there exists ν(λ) > 0 such that T (λ) + ν(λ)i is compact). If λ J is an eigenvalue of T ( ) then there exists n N such that 0 = max dimv =n min x V, x 0 T (λ)x, x. x, x In this case we assign n to the eigenvalue λ of problem T (λ)x = 0 as its number and call λ an n-th eigenvalue of T ( ). TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

79 Nonoverdamped problems Minmax characterization (V.&B.Werner 1982, V. 2010) Let T (λ), λ J be a family of linear self-adjoint and bounded operators on a Hilbert space H depending continuously on a parameter λ J where J is an open real (not necessarily bounded) interval. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

80 Nonoverdamped problems Minmax characterization (V.&B.Werner 1982, V. 2010) Let T (λ), λ J be a family of linear self-adjoint and bounded operators on a Hilbert space H depending continuously on a parameter λ J where J is an open real (not necessarily bounded) interval. Assume that for x 0 f (λ, x) := T (λ)x, x = 0 has at most one solution p(x) J, and let D be the domain of the Rayleigh functional p, (λ p(x))f (λ, x) > 0 for every x D and λ p(x), and that the supremum of the essential spectrum of T (λ) is negative for every λ J. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

81 Nonoverdamped problems Minmax characterization (V.&B.Werner 1982, V. 2010) Let T (λ), λ J be a family of linear self-adjoint and bounded operators on a Hilbert space H depending continuously on a parameter λ J where J is an open real (not necessarily bounded) interval. Assume that for x 0 f (λ, x) := T (λ)x, x = 0 has at most one solution p(x) J, and let D be the domain of the Rayleigh functional p, (λ p(x))f (λ, x) > 0 for every x D and λ p(x), and that the supremum of the essential spectrum of T (λ) is negative for every λ J. Then the nonlinear eigenvalue problem T (λ)x = 0 has at most a countable set of eigenvalues in J, and it holds that: If λ n J is an n-th eigenvalue then λ n = min dim V =n, V D sup x D V p(x). (6) TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

82 Nonoverdamped problems Minmax characterization (V.&B.Werner 1982, V. 2010) Let T (λ), λ J be a family of linear self-adjoint and bounded operators on a Hilbert space H depending continuously on a parameter λ J where J is an open real (not necessarily bounded) interval. Assume that for x 0 f (λ, x) := T (λ)x, x = 0 has at most one solution p(x) J, and let D be the domain of the Rayleigh functional p, (λ p(x))f (λ, x) > 0 for every x D and λ p(x), and that the supremum of the essential spectrum of T (λ) is negative for every λ J. Then the nonlinear eigenvalue problem T (λ)x = 0 has at most a countable set of eigenvalues in J, and it holds that: If λ n J is an n-th eigenvalue then If conversely λ n = λ n = min dim V =n, V D inf dim V =n, V D sup x D V sup x D V p(x). (6) p(x) J then λ n is an n-th eigenvalue of T (λ)x = 0 and (6) holds. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

83 Sketch of proof Variational characterization of eigenvalue problems Nonoverdamped problems Step 1 (technical): Let λ J, and assume that V is a finite dimensional subspace of H such that V D. Then it holds that < < λ = sup p(x) min T (λ)x, x = > x V D(p) x V 0 (7) > TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

84 Sketch of proof Variational characterization of eigenvalue problems Nonoverdamped problems Step 1 (technical): Let λ J, and assume that V is a finite dimensional subspace of H such that V D. Then it holds that < < λ = sup p(x) min T (λ)x, x = > x V D(p) x V 0 (7) > Step 2: If λ n is an n-th eigenvalue, then µ n (λ n ) = 0, and µ n (λ n ) = max dim V =n min T (λ n)x, x = x V, x =1 min T (λ n )x, x. x V, x =1 TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

85 Sketch of proof Variational characterization of eigenvalue problems Nonoverdamped problems Step 1 (technical): Let λ J, and assume that V is a finite dimensional subspace of H such that V D. Then it holds that < < λ = sup p(x) min T (λ)x, x = > x V D(p) x V 0 (7) > Step 2: If λ n is an n-th eigenvalue, then µ n (λ n ) = 0, and µ n (λ n ) = max dim V =n min T (λ n)x, x = x V, x =1 min T (λ n )x, x. x V, x =1 Hence, min x V, x =1 T (λ n )x, x 0 for every V with dim V = n, and (7) implies sup x V D p(x) λ n = sup x V D p(x). Hence, λ n is a minmax value of p. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

86 Nonoverdamped problems Theorem (minmax for extreme eigenvalues) Assume that the conditions of the minmax characterization hold and that inf p(x) J. x D TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

87 Nonoverdamped problems Theorem (minmax for extreme eigenvalues) Assume that the conditions of the minmax characterization hold and that inf p(x) J. x D If λ n J for some n N then every V H j with V D(p) and λ j = sup x V D(p) p(x) is contained in D {0}, and the characterization (6) can be replaced by λ j = min dim V = j V D {0} max p(v), j = 1,..., n. (8) v V, x 0 TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

88 Nonoverdamped problems Maxmin characterization (V. 2003) Assume that the conditions of the minmax characterization are satisfied. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

89 Nonoverdamped problems Maxmin characterization (V. 2003) Assume that the conditions of the minmax characterization are satisfied. If there is an n-th eigenvalue λ n J of T (λ)x = 0, then λ n = max V H n 1 V D inf p(v), v V D and the maximum is attained by W := span{u 1,..., u n 1 } where u j denotes an eigenvector corresponding to the j-largest eigenvalue µ j (λ n ) of T (λ n ). TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

90 Nonoverdamped problems Maxmin characterization (V. 2003) Assume that the conditions of the minmax characterization are satisfied. If there is an n-th eigenvalue λ n J of T (λ)x = 0, then λ n = max V H n 1 V D inf p(v), v V D and the maximum is attained by W := span{u 1,..., u n 1 } where u j denotes an eigenvector corresponding to the j-largest eigenvalue µ j (λ n ) of T (λ n ). Proof takes advantage of the following Lemma: Let λ J, and let V be a finite dimensional subspace of H such that V D. Then it holds that λ < = > inf x V D(p) p(x) max x V, x =1 T (λ)x, x < = > 0 TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

91 Nonoverdamped problems Rayleigh s principle for nonoverdamped problems Under the conditions of the minmax characterization assume that J contains n 1 eigenvalues λ 1 λ n (where λ i is an ith eigenvalue) with corresponding [, ] orthogonal eigenvectors x 1,..., x n. If there exists x D with [x i, x] = 0 for i = 1,..., n then J contains an (n + 1)th eigenvalue, and λ n+1 = inf{p(x) : [x j, x] = 0, i = 1,..., n}. (9) TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

92 Nonoverdamped problems Rayleigh s principle for nonoverdamped problems Under the conditions of the minmax characterization assume that J contains n 1 eigenvalues λ 1 λ n (where λ i is an ith eigenvalue) with corresponding [, ] orthogonal eigenvectors x 1,..., x n. If there exists x D with [x i, x] = 0 for i = 1,..., n then J contains an (n + 1)th eigenvalue, and λ n+1 = inf{p(x) : [x j, x] = 0, i = 1,..., n}. (9) The generalized scalar product (5) has to be modified according to T (p(x)) T (p(y)) x, y, if p(x) p(y) [x, y] := p(x) p(y) x, y, if p(x) = p(y) if T is not differentiable at p(x) = p(y). TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

93 Sketch of proof Variational characterization of eigenvalue problems Nonoverdamped problems Let {y k } be a minimizing sequence of (9) such that y k = 1, [x j, y k ] = 0, j = 1,..., n, p(y k ) λ n+1. TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

94 Sketch of proof Variational characterization of eigenvalue problems Nonoverdamped problems Let {y k } be a minimizing sequence of (9) such that y k = 1, [x j, y k ] = 0, j = 1,..., n, p(y k ) λ n+1. There exist unique n ỹ k = y k + c kj x j such that T (λ n+1 )x j, ỹ k = 0, j = 1,..., n, and it holds that y k ỹ k 0 and p(ỹ k ) λ n+1. j=1 TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

95 Sketch of proof Variational characterization of eigenvalue problems Nonoverdamped problems Let {y k } be a minimizing sequence of (9) such that y k = 1, [x j, y k ] = 0, j = 1,..., n, p(y k ) λ n+1. There exist unique n ỹ k = y k + c kj x j such that T (λ n+1 )x j, ỹ k = 0, j = 1,..., n, and it holds that y k ỹ k 0 and p(ỹ k ) λ n+1. j=1 For V k := span{x 1,..., x n, ỹ k } it can be shown that lim sup{p(z) : z V k} = λ n+1, k and it follows from the minmax characterization that λ n+1 is an eigenvalue of T ( ). TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

96 Further literature Nonoverdamped problems Minmax characterizations for non-overdamped nonlinear eigenvalue problems were proved (independently from our work) by Barston (1974) for some extreme eigenvalues of finite dimensional quadratic eigenproblems TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

97 Further literature Nonoverdamped problems Minmax characterizations for non-overdamped nonlinear eigenvalue problems were proved (independently from our work) by Barston (1974) for some extreme eigenvalues of finite dimensional quadratic eigenproblems and for the infinite dimensional case by Griniv & Meln ik (1996): T (λ) = A(λ) I, A(λ) compact TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

98 Further literature Nonoverdamped problems Minmax characterizations for non-overdamped nonlinear eigenvalue problems were proved (independently from our work) by Barston (1974) for some extreme eigenvalues of finite dimensional quadratic eigenproblems and for the infinite dimensional case by Griniv & Meln ik (1996): T (λ) = A(λ) I, A(λ) compact Binding, Eschwé & H. Langer (2000) TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

99 Further literature Nonoverdamped problems Minmax characterizations for non-overdamped nonlinear eigenvalue problems were proved (independently from our work) by Barston (1974) for some extreme eigenvalues of finite dimensional quadratic eigenproblems and for the infinite dimensional case by Griniv & Meln ik (1996): T (λ) = A(λ) I, A(λ) compact Binding, Eschwé & H. Langer (2000) Eschwé & M. Langer (2004) unbounded operators TUHH Heinrich Voss Variational Principles for Nonlinear EVPs Lausanne, August, / 52

CHAPTER 5 : NONLINEAR EIGENVALUE PROBLEMS

CHAPTER 5 : NONLINEAR EIGENVALUE PROBLEMS CHAPTER 5 : NONLINEAR EIGENVALUE PROBLEMS Heinrich Voss voss@tu-harburg.de Hamburg University of Technology Institute of Mathematics TUHH Heinrich Voss Nonlinear eigenvalue problems Eigenvalue problems

More information

Computable error bounds for nonlinear eigenvalue problems allowing for a minmax characterization

Computable error bounds for nonlinear eigenvalue problems allowing for a minmax characterization Computable error bounds for nonlinear eigenvalue problems allowing for a minmax characterization Heinrich Voss voss@tuhh.de Joint work with Kemal Yildiztekin Hamburg University of Technology Institute

More information

A MAXMIN PRINCIPLE FOR NONLINEAR EIGENVALUE PROBLEMS WITH APPLICATION TO A RATIONAL SPECTRAL PROBLEM IN FLUID SOLID VIBRATION

A MAXMIN PRINCIPLE FOR NONLINEAR EIGENVALUE PROBLEMS WITH APPLICATION TO A RATIONAL SPECTRAL PROBLEM IN FLUID SOLID VIBRATION A MAXMIN PRINCIPLE FOR NONLINEAR EIGENVALUE PROBLEMS WITH APPLICATION TO A RATIONAL SPECTRAL PROBLEM IN FLUID SOLID VIBRATION HEINRICH VOSS Abstract. In this paper we prove a maxmin principle for nonlinear

More information

Iterative projection methods for sparse nonlinear eigenvalue problems

Iterative projection methods for sparse nonlinear eigenvalue problems Iterative projection methods for sparse nonlinear eigenvalue problems Heinrich Voss voss@tu-harburg.de Hamburg University of Technology Institute of Mathematics TUHH Heinrich Voss Iterative projection

More information

Dedicated to Ivo Marek on the occasion of his 75th birthday. T(λ)x = 0 (1.1)

Dedicated to Ivo Marek on the occasion of his 75th birthday. T(λ)x = 0 (1.1) A MINMAX PRINCIPLE FOR NONLINEAR EIGENPROBLEMS DEPENDING CONTINUOUSLY ON THE EIGENPARAMETER HEINRICH VOSS Key words. nonlinear eigenvalue problem, variational characterization, minmax principle, variational

More information

A Jacobi Davidson-type projection method for nonlinear eigenvalue problems

A Jacobi Davidson-type projection method for nonlinear eigenvalue problems A Jacobi Davidson-type projection method for nonlinear eigenvalue problems Timo Betce and Heinrich Voss Technical University of Hamburg-Harburg, Department of Mathematics, Schwarzenbergstrasse 95, D-21073

More information

ON SYLVESTER S LAW OF INERTIA FOR NONLINEAR EIGENVALUE PROBLEMS

ON SYLVESTER S LAW OF INERTIA FOR NONLINEAR EIGENVALUE PROBLEMS ON SYLVESTER S LAW OF INERTIA FOR NONLINEAR EIGENVALUE PROBLEMS ALEKSANDRA KOSTIĆ AND HEINRICH VOSS Key words. eigenvalue, variational characterization, principle, Sylvester s law of inertia AMS subject

More information

Solving Regularized Total Least Squares Problems

Solving Regularized Total Least Squares Problems Solving Regularized Total Least Squares Problems Heinrich Voss voss@tu-harburg.de Hamburg University of Technology Institute of Numerical Simulation Joint work with Jörg Lampe TUHH Heinrich Voss Total

More information

An Arnoldi Method for Nonlinear Symmetric Eigenvalue Problems

An Arnoldi Method for Nonlinear Symmetric Eigenvalue Problems An Arnoldi Method for Nonlinear Symmetric Eigenvalue Problems H. Voss 1 Introduction In this paper we consider the nonlinear eigenvalue problem T (λ)x = 0 (1) where T (λ) R n n is a family of symmetric

More information

Available online at ScienceDirect. Procedia Engineering 100 (2015 ) 56 63

Available online at   ScienceDirect. Procedia Engineering 100 (2015 ) 56 63 Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 100 (2015 ) 56 63 25th DAAAM International Symposium on Intelligent Manufacturing and Automation, DAAAM 2014 Definite Quadratic

More information

SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS

SPECTRAL 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 information

Efficient Methods For Nonlinear Eigenvalue Problems. Diploma Thesis

Efficient Methods For Nonlinear Eigenvalue Problems. Diploma Thesis Efficient Methods For Nonlinear Eigenvalue Problems Diploma Thesis Timo Betcke Technical University of Hamburg-Harburg Department of Mathematics (Prof. Dr. H. Voß) August 2002 Abstract During the last

More information

Eigenvalue Problems CHAPTER 1 : PRELIMINARIES

Eigenvalue Problems CHAPTER 1 : PRELIMINARIES Eigenvalue Problems CHAPTER 1 : PRELIMINARIES Heinrich Voss voss@tu-harburg.de Hamburg University of Technology Institute of Mathematics TUHH Heinrich Voss Preliminaries Eigenvalue problems 2012 1 / 14

More information

Spectral Properties of Elliptic Operators In previous work we have replaced the strong version of an elliptic boundary value problem

Spectral Properties of Elliptic Operators In previous work we have replaced the strong version of an elliptic boundary value problem Spectral Properties of Elliptic Operators In previous work we have replaced the strong version of an elliptic boundary value problem L u x f x BC u x g x with the weak problem find u V such that B u,v

More information

Analysis Preliminary Exam Workshop: Hilbert Spaces

Analysis Preliminary Exam Workshop: Hilbert Spaces Analysis Preliminary Exam Workshop: Hilbert Spaces 1. Hilbert spaces A Hilbert space H is a complete real or complex inner product space. Consider complex Hilbert spaces for definiteness. If (, ) : H H

More information

SPECTRAL THEOREM FOR COMPACT SELF-ADJOINT OPERATORS

SPECTRAL THEOREM FOR COMPACT SELF-ADJOINT OPERATORS SPECTRAL THEOREM FOR COMPACT SELF-ADJOINT OPERATORS G. RAMESH Contents Introduction 1 1. Bounded Operators 1 1.3. Examples 3 2. Compact Operators 5 2.1. Properties 6 3. The Spectral Theorem 9 3.3. Self-adjoint

More information

Linear algebra and applications to graphs Part 1

Linear algebra and applications to graphs Part 1 Linear algebra and applications to graphs Part 1 Written up by Mikhail Belkin and Moon Duchin Instructor: Laszlo Babai June 17, 2001 1 Basic Linear Algebra Exercise 1.1 Let V and W be linear subspaces

More information

STATIONARY SCHRÖDINGER EQUATIONS GOVERNING ELECTRONIC STATES OF QUANTUM DOTS IN THE PRESENCE OF SPIN ORBIT SPLITTING

STATIONARY SCHRÖDINGER EQUATIONS GOVERNING ELECTRONIC STATES OF QUANTUM DOTS IN THE PRESENCE OF SPIN ORBIT SPLITTING STATIONARY SCHRÖDINGER EQUATIONS GOVERNING ELECTRONIC STATES OF QUANTUM DOTS IN THE PRESENCE OF SPIN ORBIT SPLITTING MARTA BETCKE AND HEINRICH VOSS Abstract. In this work we derive a pair of nonlinear

More information

Asymptotic distribution of eigenvalues of Laplace operator

Asymptotic distribution of eigenvalues of Laplace operator Asymptotic distribution of eigenvalues of Laplace operator 23.8.2013 Topics We will talk about: the number of eigenvalues of Laplace operator smaller than some λ as a function of λ asymptotic behaviour

More information

Mathematical foundations - linear algebra

Mathematical foundations - linear algebra Mathematical foundations - linear algebra Andrea Passerini passerini@disi.unitn.it Machine Learning Vector space Definition (over reals) A set X is called a vector space over IR if addition and scalar

More information

Second Order Elliptic PDE

Second Order Elliptic PDE Second Order Elliptic PDE T. Muthukumar tmk@iitk.ac.in December 16, 2014 Contents 1 A Quick Introduction to PDE 1 2 Classification of Second Order PDE 3 3 Linear Second Order Elliptic Operators 4 4 Periodic

More information

08a. Operators on Hilbert spaces. 1. Boundedness, continuity, operator norms

08a. Operators on Hilbert spaces. 1. Boundedness, continuity, operator norms (February 24, 2017) 08a. Operators on Hilbert spaces Paul Garrett garrett@math.umn.edu http://www.math.umn.edu/ garrett/ [This document is http://www.math.umn.edu/ garrett/m/real/notes 2016-17/08a-ops

More information

Your first day at work MATH 806 (Fall 2015)

Your first day at work MATH 806 (Fall 2015) Your first day at work MATH 806 (Fall 2015) 1. Let X be a set (with no particular algebraic structure). A function d : X X R is called a metric on X (and then X is called a metric space) when d satisfies

More information

LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)

LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,

More information

MAT 578 FUNCTIONAL ANALYSIS EXERCISES

MAT 578 FUNCTIONAL ANALYSIS EXERCISES MAT 578 FUNCTIONAL ANALYSIS EXERCISES JOHN QUIGG Exercise 1. Prove that if A is bounded in a topological vector space, then for every neighborhood V of 0 there exists c > 0 such that tv A for all t > c.

More information

The following definition is fundamental.

The following definition is fundamental. 1. Some Basics from Linear Algebra With these notes, I will try and clarify certain topics that I only quickly mention in class. First and foremost, I will assume that you are familiar with many basic

More information

The converse is clear, since

The converse is clear, since 14. The minimal polynomial For an example of a matrix which cannot be diagonalised, consider the matrix ( ) 0 1 A =. 0 0 The characteristic polynomial is λ 2 = 0 so that the only eigenvalue is λ = 0. The

More information

QUALIFYING EXAMINATION Harvard University Department of Mathematics Tuesday September 21, 2004 (Day 1)

QUALIFYING EXAMINATION Harvard University Department of Mathematics Tuesday September 21, 2004 (Day 1) QUALIFYING EXAMINATION Harvard University Department of Mathematics Tuesday September 21, 2004 (Day 1) Each of the six questions is worth 10 points. 1) Let H be a (real or complex) Hilbert space. We say

More information

ITERATIVE PROJECTION METHODS FOR SPARSE LINEAR SYSTEMS AND EIGENPROBLEMS CHAPTER 3 : SEMI-ITERATIVE METHODS

ITERATIVE PROJECTION METHODS FOR SPARSE LINEAR SYSTEMS AND EIGENPROBLEMS CHAPTER 3 : SEMI-ITERATIVE METHODS ITERATIVE PROJECTION METHODS FOR SPARSE LINEAR SYSTEMS AND EIGENPROBLEMS CHAPTER 3 : SEMI-ITERATIVE METHODS Heinrich Voss voss@tu-harburg.de Hamburg University of Technology Institute of Numerical Simulation

More information

Matrix Theory, Math6304 Lecture Notes from March 22, 2016 taken by Kazem Safari

Matrix Theory, Math6304 Lecture Notes from March 22, 2016 taken by Kazem Safari Matrix Theory, Math6304 Lecture Notes from March 22, 2016 taken by Kazem Safari 1.1 Applications of Courant-Fisher and min- or -min Last time: Courant -Fishert -min or min- for eigenvalues Warm-up: Sums

More information

ACM/CMS 107 Linear Analysis & Applications Fall 2017 Assignment 2: PDEs and Finite Element Methods Due: 7th November 2017

ACM/CMS 107 Linear Analysis & Applications Fall 2017 Assignment 2: PDEs and Finite Element Methods Due: 7th November 2017 ACM/CMS 17 Linear Analysis & Applications Fall 217 Assignment 2: PDEs and Finite Element Methods Due: 7th November 217 For this assignment the following MATLAB code will be required: Introduction http://wwwmdunloporg/cms17/assignment2zip

More information

Solving Polynomial Eigenproblems by Linearization

Solving Polynomial Eigenproblems by Linearization Solving Polynomial Eigenproblems by Linearization Nick Higham School of Mathematics University of Manchester higham@ma.man.ac.uk http://www.ma.man.ac.uk/~higham/ Joint work with D. Steven Mackey and Françoise

More information

Econ Slides from Lecture 7

Econ Slides from Lecture 7 Econ 205 Sobel Econ 205 - Slides from Lecture 7 Joel Sobel August 31, 2010 Linear Algebra: Main Theory A linear combination of a collection of vectors {x 1,..., x k } is a vector of the form k λ ix i for

More information

Ma/CS 6b Class 23: Eigenvalues in Regular Graphs

Ma/CS 6b Class 23: Eigenvalues in Regular Graphs Ma/CS 6b Class 3: Eigenvalues in Regular Graphs By Adam Sheffer Recall: The Spectrum of a Graph Consider a graph G = V, E and let A be the adjacency matrix of G. The eigenvalues of G are the eigenvalues

More information

Chap 3. Linear Algebra

Chap 3. Linear Algebra Chap 3. Linear Algebra Outlines 1. Introduction 2. Basis, Representation, and Orthonormalization 3. Linear Algebraic Equations 4. Similarity Transformation 5. Diagonal Form and Jordan Form 6. Functions

More information

Final Review Written by Victoria Kala SH 6432u Office Hours R 12:30 1:30pm Last Updated 11/30/2015

Final Review Written by Victoria Kala SH 6432u Office Hours R 12:30 1:30pm Last Updated 11/30/2015 Final Review Written by Victoria Kala vtkala@mathucsbedu SH 6432u Office Hours R 12:30 1:30pm Last Updated 11/30/2015 Summary This review contains notes on sections 44 47, 51 53, 61, 62, 65 For your final,

More information

AA242B: MECHANICAL VIBRATIONS

AA242B: MECHANICAL VIBRATIONS AA242B: MECHANICAL VIBRATIONS 1 / 17 AA242B: MECHANICAL VIBRATIONS Solution Methods for the Generalized Eigenvalue Problem These slides are based on the recommended textbook: M. Géradin and D. Rixen, Mechanical

More information

Appendix A Functional Analysis

Appendix A Functional Analysis Appendix A Functional Analysis A.1 Metric Spaces, Banach Spaces, and Hilbert Spaces Definition A.1. Metric space. Let X be a set. A map d : X X R is called metric on X if for all x,y,z X it is i) d(x,y)

More information

LECTURE 7. k=1 (, v k)u k. Moreover r

LECTURE 7. k=1 (, v k)u k. Moreover r LECTURE 7 Finite rank operators Definition. T is said to be of rank r (r < ) if dim T(H) = r. The class of operators of rank r is denoted by K r and K := r K r. Theorem 1. T K r iff T K r. Proof. Let T

More information

Study Guide for Linear Algebra Exam 2

Study Guide for Linear Algebra Exam 2 Study Guide for Linear Algebra Exam 2 Term Vector Space Definition A Vector Space is a nonempty set V of objects, on which are defined two operations, called addition and multiplication by scalars (real

More information

LINEAR ALGEBRA REVIEW

LINEAR ALGEBRA REVIEW LINEAR ALGEBRA REVIEW JC Stuff you should know for the exam. 1. Basics on vector spaces (1) F n is the set of all n-tuples (a 1,... a n ) with a i F. It forms a VS with the operations of + and scalar multiplication

More information

Hilbert 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. 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 information

I teach myself... Hilbert spaces

I teach myself... Hilbert spaces I teach myself... Hilbert spaces by F.J.Sayas, for MATH 806 November 4, 2015 This document will be growing with the semester. Every in red is for you to justify. Even if we start with the basic definition

More information

LINEAR ALGEBRA SUMMARY SHEET.

LINEAR ALGEBRA SUMMARY SHEET. LINEAR ALGEBRA SUMMARY SHEET RADON ROSBOROUGH https://intuitiveexplanationscom/linear-algebra-summary-sheet/ This document is a concise collection of many of the important theorems of linear algebra, organized

More information

Solving a rational eigenvalue problem in fluidstructure

Solving a rational eigenvalue problem in fluidstructure Solving a rational eigenvalue problem in fluidstructure interaction H. VOSS Section of Mathematics, TU Germany Hambu~g-Harburg, D-21071 Hamburg, Abstract In this paper we consider a rational eigenvalue

More information

Linear Algebra: Matrix Eigenvalue Problems

Linear 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 information

Ordinary Differential Equations II

Ordinary Differential Equations II Ordinary Differential Equations II February 23 2017 Separation of variables Wave eq. (PDE) 2 u t (t, x) = 2 u 2 c2 (t, x), x2 c > 0 constant. Describes small vibrations in a homogeneous string. u(t, x)

More information

Math 108b: Notes on the Spectral Theorem

Math 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 information

Ir O D = D = ( ) Section 2.6 Example 1. (Bottom of page 119) dim(v ) = dim(l(v, W )) = dim(v ) dim(f ) = dim(v )

Ir O D = D = ( ) Section 2.6 Example 1. (Bottom of page 119) dim(v ) = dim(l(v, W )) = dim(v ) dim(f ) = dim(v ) Section 3.2 Theorem 3.6. Let A be an m n matrix of rank r. Then r m, r n, and, by means of a finite number of elementary row and column operations, A can be transformed into the matrix ( ) Ir O D = 1 O

More information

Functional Analysis Review

Functional 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 information

In particular, if A is a square matrix and λ is one of its eigenvalues, then we can find a non-zero column vector X with

In particular, if A is a square matrix and λ is one of its eigenvalues, then we can find a non-zero column vector X with Appendix: Matrix Estimates and the Perron-Frobenius Theorem. This Appendix will first present some well known estimates. For any m n matrix A = [a ij ] over the real or complex numbers, it will be convenient

More information

SPECTRAL THEOREM FOR SYMMETRIC OPERATORS WITH COMPACT RESOLVENT

SPECTRAL THEOREM FOR SYMMETRIC OPERATORS WITH COMPACT RESOLVENT SPECTRAL THEOREM FOR SYMMETRIC OPERATORS WITH COMPACT RESOLVENT Abstract. These are the letcure notes prepared for the workshop on Functional Analysis and Operator Algebras to be held at NIT-Karnataka,

More information

MATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators.

MATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators. MATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators. Adjoint operator and adjoint matrix Given a linear operator L on an inner product space V, the adjoint of L is a transformation

More information

Introduction to Quantum Spin Systems

Introduction to Quantum Spin Systems 1 Introduction to Quantum Spin Systems Lecture 2 Sven Bachmann (standing in for Bruno Nachtergaele) Mathematics, UC Davis MAT290-25, CRN 30216, Winter 2011, 01/10/11 2 Basic Setup For concreteness, consider

More information

The quantum state as a vector

The 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 information

Lecture Notes on PDEs

Lecture Notes on PDEs Lecture Notes on PDEs Alberto Bressan February 26, 2012 1 Elliptic equations Let IR n be a bounded open set Given measurable functions a ij, b i, c : IR, consider the linear, second order differential

More information

Homework 11 Solutions. Math 110, Fall 2013.

Homework 11 Solutions. Math 110, Fall 2013. Homework 11 Solutions Math 110, Fall 2013 1 a) Suppose that T were self-adjoint Then, the Spectral Theorem tells us that there would exist an orthonormal basis of P 2 (R), (p 1, p 2, p 3 ), consisting

More information

Reduction of nonlinear eigenproblems with JD

Reduction of nonlinear eigenproblems with JD Reduction of nonlinear eigenproblems with JD Henk van der Vorst H.A.vanderVorst@math.uu.nl Mathematical Institute Utrecht University July 13, 2005, SIAM Annual New Orleans p.1/16 Outline Polynomial Eigenproblems

More information

Math Linear Algebra Final Exam Review Sheet

Math Linear Algebra Final Exam Review Sheet Math 15-1 Linear Algebra Final Exam Review Sheet Vector Operations Vector addition is a component-wise operation. Two vectors v and w may be added together as long as they contain the same number n of

More information

Keeping σ fixed for several steps, iterating on µ and neglecting the remainder in the Lagrange interpolation one obtains. θ = λ j λ j 1 λ j σ, (2.

Keeping σ fixed for several steps, iterating on µ and neglecting the remainder in the Lagrange interpolation one obtains. θ = λ j λ j 1 λ j σ, (2. RATIONAL KRYLOV FOR NONLINEAR EIGENPROBLEMS, AN ITERATIVE PROJECTION METHOD ELIAS JARLEBRING AND HEINRICH VOSS Key words. nonlinear eigenvalue problem, rational Krylov, Arnoldi, projection method AMS subject

More information

Recall that any inner product space V has an associated norm defined by

Recall that any inner product space V has an associated norm defined by Hilbert Spaces Recall that any inner product space V has an associated norm defined by v = v v. Thus an inner product space can be viewed as a special kind of normed vector space. In particular every inner

More information

AMS526: Numerical Analysis I (Numerical Linear Algebra)

AMS526: Numerical Analysis I (Numerical Linear Algebra) AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 16: Rayleigh Quotient Iteration Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 10 Solving Eigenvalue Problems All

More information

Semi-Classical Dynamics Using Hagedorn Wavepackets

Semi-Classical Dynamics Using Hagedorn Wavepackets Semi-Classical Dynamics Using Hagedorn Wavepackets Leila Taghizadeh June 6, 2013 Leila Taghizadeh () Semi-Classical Dynamics Using Hagedorn Wavepackets June 6, 2013 1 / 56 Outline 1 The Schrödinger Equation

More information

Functional Analysis. Martin Brokate. 1 Normed Spaces 2. 2 Hilbert Spaces The Principle of Uniform Boundedness 32

Functional Analysis. Martin Brokate. 1 Normed Spaces 2. 2 Hilbert Spaces The Principle of Uniform Boundedness 32 Functional Analysis Martin Brokate Contents 1 Normed Spaces 2 2 Hilbert Spaces 2 3 The Principle of Uniform Boundedness 32 4 Extension, Reflexivity, Separation 37 5 Compact subsets of C and L p 46 6 Weak

More information

On the spectrum of the Laplacian

On the spectrum of the Laplacian On the spectrum of the Laplacian S. Kesavan Institute of Mathematical Sciences, CIT Campus, Taramani, Chennai - 600113. kesh@imsc.res.in July 1, 2013 S. Kesavan (IMSc) Spectrum of the Laplacian July 1,

More information

Nonlinear Eigenvalue Problems: An Introduction

Nonlinear Eigenvalue Problems: An Introduction Nonlinear Eigenvalue Problems: An Introduction Cedric Effenberger Seminar for Applied Mathematics ETH Zurich Pro*Doc Workshop Disentis, August 18 21, 2010 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction

More information

Math 4A Notes. Written by Victoria Kala Last updated June 11, 2017

Math 4A Notes. Written by Victoria Kala Last updated June 11, 2017 Math 4A Notes Written by Victoria Kala vtkala@math.ucsb.edu Last updated June 11, 2017 Systems of Linear Equations A linear equation is an equation that can be written in the form a 1 x 1 + a 2 x 2 +...

More information

3.024 Electrical, Optical, and Magnetic Properties of Materials Spring 2012 Recitation 3 Notes

3.024 Electrical, Optical, and Magnetic Properties of Materials Spring 2012 Recitation 3 Notes 3.024 Electrical, Optical, and Magnetic Properties of Materials Spring 2012 Outline 1. Schr dinger: Eigenfunction Problems & Operator Properties 2. Piecewise Function/Continuity Review -Scattering from

More information

CONTENTS. 4 Hausdorff Measure Introduction The Cantor Set Rectifiable Curves Cantor Set-Like Objects...

CONTENTS. 4 Hausdorff Measure Introduction The Cantor Set Rectifiable Curves Cantor Set-Like Objects... Contents 1 Functional Analysis 1 1.1 Hilbert Spaces................................... 1 1.1.1 Spectral Theorem............................. 4 1.2 Normed Vector Spaces.............................. 7 1.2.1

More information

arxiv: v2 [math.oc] 31 Jul 2017

arxiv: v2 [math.oc] 31 Jul 2017 An unconstrained framework for eigenvalue problems Yunho Kim arxiv:6.09707v [math.oc] 3 Jul 07 May 0, 08 Abstract In this paper, we propose an unconstrained framework for eigenvalue problems in both discrete

More information

5 Compact linear operators

5 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 information

Introduction to Electronic Structure Theory

Introduction to Electronic Structure Theory Introduction to Electronic Structure Theory C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology June 2002 Last Revised: June 2003 1 Introduction The purpose of these

More information

Abstract Vector Spaces and Concrete Examples

Abstract Vector Spaces and Concrete Examples LECTURE 18 Abstract Vector Spaces and Concrete Examples Our discussion of linear algebra so far has been devoted to discussing the relations between systems of linear equations, matrices, and vectors.

More information

1. Select the unique answer (choice) for each problem. Write only the answer.

1. Select the unique answer (choice) for each problem. Write only the answer. MATH 5 Practice Problem Set Spring 7. Select the unique answer (choice) for each problem. Write only the answer. () Determine all the values of a for which the system has infinitely many solutions: x +

More information

Numerical Methods - Numerical Linear Algebra

Numerical Methods - Numerical Linear Algebra Numerical Methods - Numerical Linear Algebra Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Numerical Linear Algebra I 2013 1 / 62 Outline 1 Motivation 2 Solving Linear

More information

IMPORTANT DEFINITIONS AND THEOREMS REFERENCE SHEET

IMPORTANT DEFINITIONS AND THEOREMS REFERENCE SHEET IMPORTANT DEFINITIONS AND THEOREMS REFERENCE SHEET This is a (not quite comprehensive) list of definitions and theorems given in Math 1553. Pay particular attention to the ones in red. Study Tip For each

More information

Transpose & Dot Product

Transpose & Dot Product Transpose & Dot Product Def: The transpose of an m n matrix A is the n m matrix A T whose columns are the rows of A. So: The columns of A T are the rows of A. The rows of A T are the columns of A. Example:

More information

DS-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 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 information

Math 102, Winter Final Exam Review. Chapter 1. Matrices and Gaussian Elimination

Math 102, Winter Final Exam Review. Chapter 1. Matrices and Gaussian Elimination Math 0, Winter 07 Final Exam Review Chapter. Matrices and Gaussian Elimination { x + x =,. Different forms of a system of linear equations. Example: The x + 4x = 4. [ ] [ ] [ ] vector form (or the column

More information

WOMP 2001: LINEAR ALGEBRA. 1. Vector spaces

WOMP 2001: LINEAR ALGEBRA. 1. Vector spaces WOMP 2001: LINEAR ALGEBRA DAN GROSSMAN Reference Roman, S Advanced Linear Algebra, GTM #135 (Not very good) Let k be a field, eg, R, Q, C, F q, K(t), 1 Vector spaces Definition A vector space over k is

More information

Quadratic forms. Here. Thus symmetric matrices are diagonalizable, and the diagonalization can be performed by means of an orthogonal matrix.

Quadratic forms. Here. Thus symmetric matrices are diagonalizable, and the diagonalization can be performed by means of an orthogonal matrix. Quadratic forms 1. Symmetric matrices An n n matrix (a ij ) n ij=1 with entries on R is called symmetric if A T, that is, if a ij = a ji for all 1 i, j n. We denote by S n (R) the set of all n n symmetric

More information

Lecture 4 Eigenvalue problems

Lecture 4 Eigenvalue problems Lecture 4 Eigenvalue problems Weinan E 1,2 and Tiejun Li 2 1 Department of Mathematics, Princeton University, weinan@princeton.edu 2 School of Mathematical Sciences, Peking University, tieli@pku.edu.cn

More information

Further Mathematical Methods (Linear Algebra)

Further Mathematical Methods (Linear Algebra) Further Mathematical Methods (Linear Algebra) Solutions For The 2 Examination Question (a) For a non-empty subset W of V to be a subspace of V we require that for all vectors x y W and all scalars α R:

More information

Vectors To begin, let us describe an element of the state space as a point with numerical coordinates, that is x 1. x 2. x =

Vectors To begin, let us describe an element of the state space as a point with numerical coordinates, that is x 1. x 2. x = Linear Algebra Review Vectors To begin, let us describe an element of the state space as a point with numerical coordinates, that is x 1 x x = 2. x n Vectors of up to three dimensions are easy to diagram.

More information

APPENDIX 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.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 information

Chapter 2 Formulas and Definitions:

Chapter 2 Formulas and Definitions: Chapter 2 Formulas and Definitions: (from 2.1) Definition of Polynomial Function: Let n be a nonnegative integer and let a n,a n 1,...,a 2,a 1,a 0 be real numbers with a n 0. The function given by f (x)

More information

Transpose & Dot Product

Transpose & Dot Product Transpose & Dot Product Def: The transpose of an m n matrix A is the n m matrix A T whose columns are the rows of A. So: The columns of A T are the rows of A. The rows of A T are the columns of A. Example:

More information

************************************* Applied Analysis I - (Advanced PDE I) (Math 940, Fall 2014) Baisheng Yan

************************************* Applied Analysis I - (Advanced PDE I) (Math 940, Fall 2014) Baisheng Yan ************************************* Applied Analysis I - (Advanced PDE I) (Math 94, Fall 214) by Baisheng Yan Department of Mathematics Michigan State University yan@math.msu.edu Contents Chapter 1.

More information

Lecture 8 : Eigenvalues and Eigenvectors

Lecture 8 : Eigenvalues and Eigenvectors CPS290: Algorithmic Foundations of Data Science February 24, 2017 Lecture 8 : Eigenvalues and Eigenvectors Lecturer: Kamesh Munagala Scribe: Kamesh Munagala Hermitian Matrices It is simpler to begin with

More information

1 The postulates of quantum mechanics

1 The postulates of quantum mechanics 1 The postulates of quantum mechanics The postulates of quantum mechanics were derived after a long process of trial and error. These postulates provide a connection between the physical world and the

More information

MTH 503: Functional Analysis

MTH 503: Functional Analysis MTH 53: Functional Analysis Semester 1, 215-216 Dr. Prahlad Vaidyanathan Contents I. Normed Linear Spaces 4 1. Review of Linear Algebra........................... 4 2. Definition and Examples...........................

More information

Further Mathematical Methods (Linear Algebra)

Further Mathematical Methods (Linear Algebra) Further Mathematical Methods (Linear Algebra) Solutions For The Examination Question (a) To be an inner product on the real vector space V a function x y which maps vectors x y V to R must be such that:

More information

CONVERGENCE THEORY. G. ALLAIRE CMAP, Ecole Polytechnique. 1. Maximum principle. 2. Oscillating test function. 3. Two-scale convergence

CONVERGENCE THEORY. G. ALLAIRE CMAP, Ecole Polytechnique. 1. Maximum principle. 2. Oscillating test function. 3. Two-scale convergence 1 CONVERGENCE THEOR G. ALLAIRE CMAP, Ecole Polytechnique 1. Maximum principle 2. Oscillating test function 3. Two-scale convergence 4. Application to homogenization 5. General theory H-convergence) 6.

More information

ABSTRACT. Name: Daniil R. Sarkissian Department: Mathematical Sciences

ABSTRACT. Name: Daniil R. Sarkissian Department: Mathematical Sciences ABSTRACT Name: Daniil R. Sarkissian Department: Mathematical Sciences Title: Theory and Computations of Partial Eigenvalue and Eigenstructure Assignment Problems in Matrix Second-Order and Distributed-Parameter

More information

Recitation 1 (Sep. 15, 2017)

Recitation 1 (Sep. 15, 2017) Lecture 1 8.321 Quantum Theory I, Fall 2017 1 Recitation 1 (Sep. 15, 2017) 1.1 Simultaneous Diagonalization In the last lecture, we discussed the situations in which two operators can be simultaneously

More information

IMPORTANT DEFINITIONS AND THEOREMS REFERENCE SHEET

IMPORTANT DEFINITIONS AND THEOREMS REFERENCE SHEET IMPORTANT DEFINITIONS AND THEOREMS REFERENCE SHEET This is a (not quite comprehensive) list of definitions and theorems given in Math 1553. Pay particular attention to the ones in red. Study Tip For each

More information

THE SINGULAR VALUE DECOMPOSITION AND LOW RANK APPROXIMATION

THE SINGULAR VALUE DECOMPOSITION AND LOW RANK APPROXIMATION THE SINGULAR VALUE DECOMPOSITION AND LOW RANK APPROXIMATION MANTAS MAŽEIKA Abstract. The purpose of this paper is to present a largely self-contained proof of the singular value decomposition (SVD), and

More information

STABILITY. Phase portraits and local stability

STABILITY. Phase portraits and local stability MAS271 Methods for differential equations Dr. R. Jain STABILITY Phase portraits and local stability We are interested in system of ordinary differential equations of the form ẋ = f(x, y), ẏ = g(x, y),

More information

Exact Envelope Function Theory Band Structure of Semiconductor Heterostructure

Exact Envelope Function Theory Band Structure of Semiconductor Heterostructure University of Southern Denmark Mads Clausen Institute Exact Envelope Function Theory Band Structure of Semiconductor Heterostructure Daniele Barettin daniele@mci.sdu.dk Summary Introduction k p Model Exact

More information

RNDr. Petr Tomiczek CSc.

RNDr. Petr Tomiczek CSc. HABILITAČNÍ PRÁCE RNDr. Petr Tomiczek CSc. Plzeň 26 Nonlinear differential equation of second order RNDr. Petr Tomiczek CSc. Department of Mathematics, University of West Bohemia 5 Contents 1 Introduction

More information