Variational Principles for Nonlinear Eigenvalue Problems
|
|
- Richard Rudolf Carroll
- 6 years ago
- Views:
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 Heinrich Voss voss@tu-harburg.de Hamburg University of Technology Institute of Mathematics TUHH Heinrich Voss Nonlinear eigenvalue problems Eigenvalue problems
More informationComputable 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 informationA 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 informationIterative 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 informationDedicated 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 informationA 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 informationON 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 informationSolving 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 informationAn 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 informationAvailable 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 informationSPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS
SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS TSOGTGEREL GANTUMUR Abstract. After establishing discrete spectra for a large class of elliptic operators, we present some fundamental spectral properties
More informationEfficient 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 informationEigenvalue 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 informationSpectral 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 informationAnalysis 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 informationSPECTRAL 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 informationLinear 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 informationSTATIONARY 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 informationAsymptotic 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 informationMathematical 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 informationSecond 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 information08a. 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 informationYour 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 informationLECTURE # 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 informationMAT 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 informationThe 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 informationThe 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 informationQUALIFYING 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 informationITERATIVE 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 informationMatrix 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 informationACM/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 informationSolving 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 informationEcon 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 informationMa/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 informationChap 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 informationFinal 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 informationAA242B: 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 informationAppendix 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 informationLECTURE 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 informationStudy 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 informationLINEAR 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 informationHilbert Spaces. Hilbert space is a vector space with some extra structure. We start with formal (axiomatic) definition of a vector space.
Hilbert Spaces Hilbert space is a vector space with some extra structure. We start with formal (axiomatic) definition of a vector space. Vector Space. Vector space, ν, over the field of complex numbers,
More informationI 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 informationLINEAR 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 informationSolving 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 informationLinear Algebra: Matrix Eigenvalue Problems
CHAPTER8 Linear Algebra: Matrix Eigenvalue Problems Chapter 8 p1 A matrix eigenvalue problem considers the vector equation (1) Ax = λx. 8.0 Linear Algebra: Matrix Eigenvalue Problems Here A is a given
More informationOrdinary 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 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 informationIr 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 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 informationIn 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 informationSPECTRAL 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 informationMATH 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 informationIntroduction 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 informationThe quantum state as a vector
The quantum state as a vector February 6, 27 Wave mechanics In our review of the development of wave mechanics, we have established several basic properties of the quantum description of nature:. A particle
More informationLecture 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 informationHomework 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 informationReduction 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 informationMath 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 informationKeeping σ 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 informationRecall 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 informationAMS526: 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 informationSemi-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 informationFunctional 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 informationOn 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 informationNonlinear 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 informationMath 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 information3.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 informationCONTENTS. 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 informationarxiv: 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 information5 Compact linear operators
5 Compact linear operators One of the most important results of Linear Algebra is that for every selfadjoint linear map A on a finite-dimensional space, there exists a basis consisting of eigenvectors.
More informationIntroduction 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 informationAbstract 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 information1. 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 informationNumerical 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 informationIMPORTANT 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 informationTranspose & 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 informationDS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra.
DS-GA 1002 Lecture notes 0 Fall 2016 Linear Algebra These notes provide a review of basic concepts in linear algebra. 1 Vector spaces You are no doubt familiar with vectors in R 2 or R 3, i.e. [ ] 1.1
More informationMath 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 informationWOMP 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 informationQuadratic 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 informationLecture 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 informationFurther 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 informationVectors 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 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 informationChapter 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 informationTranspose & 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 94, Fall 214) by Baisheng Yan Department of Mathematics Michigan State University yan@math.msu.edu Contents Chapter 1.
More informationLecture 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 information1 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 informationMTH 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 informationFurther 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 informationCONVERGENCE 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 informationABSTRACT. 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 informationRecitation 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 informationIMPORTANT 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 informationTHE 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 informationSTABILITY. 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 informationExact 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 informationRNDr. 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