arxiv: v1 [math-ph] 4 Mar 2008
|
|
- Evan Fletcher
- 6 years ago
- Views:
Transcription
1 arxiv: v1 [math-ph] 4 Mar 2008 On Spectral Analysis of an Operator Pencil with Applications to Wave Propagation in Periodic and Frequency Dependent Materials Christian Engström and Markus Richter July 2, 2008 Abstract We study wave propagation in periodic and frequency dependent materials. The approach in this paper leads to spectral analysis of a quadratic operator pencil where the spectrum parameter relates to the quasimomentum and the frequency is a parameter. We show that the underlying operator has a discrete spectrum, where the eigenvalues are symmetrically placed with respect to the real and imaginary axis. Moreover, we derive bounds on the eigenvalues that is tight in low-frequency applications. Numerical calculations in two dimensions illustrates the method. 1 Introduction In this paper, we study wave propagation in periodic and frequency dependent materials. The underlying operator has a band-structure, which by appropriate choice of the geometrical structure has gaps. Such structures have numerous applications in optics, photonics and microwave engineering. Other important applications are low-frequency wave propagation and homogenization of periodic composite materials. Electromagnetic wave propagation is in general frequency dependent, since dispersive effects such as resonances or relaxation processes are present in the material. The frequency dependence of the material is in most cases disregarded in spectral analysis. This is justified in some cases where the frequency dependence is weak, but the frequency dispersion can in general not be ignored. We study wave propagation in a dispersive media that for a frequency interval is characterized by a frequency ω and space x dependent real-valued permittivity ǫ(x, ω). The waves propagate in a non-magnetic material with the permittivity ǫ(x 1, x 2, ω) independent of the third coordinate x 3. Let E and H denote the electric and magnetic waves, respectively. The x 3 independent electromagnetic wave (E, H) is decomposed into transverse electric (TE) polarized waves (E 1, E 2, 0, 0, 0, H 3 ) and transverse magnetic (TM) polarized waves (0, 0, E 3, H 1, H 2, 0). This decomposition reduces the spectral problem for the 1
2 Maxwell operator to one scalar equation for H 3 and one equation for E 3. We consider TM polarized waves, which gives a Helmholtz type of equation in E 3, but the analysis also applies to the TE case. The approach in this paper leads to spectral analysis of a quadratic operator pencil where the spectrum parameter relates to the quasimomentum (the Floquet-Bloch wave vector) and the frequency is a parameter. We show that the underlying operator has a discrete spectrum, where the eigenvalues are symmetrically placed with respect to the real and imaginary axis (Hamiltonian symmetry). Moreover, we derive bounds on the eigenvalues that is tight in lowfrequency applications. Finally, a finite element method is used to discretize the operator pencil and the isotropic implicitly restarted Arnoldi algorithm (eigs from MATLAB) is used to compute eigenvalues of the resulting large sparse matrix pencil. 2 Bloch solutions and spectrum problems Let Γ denote the lattice 2πZ 2 and ǫ a measurable function with ǫ(x + γ, ω) = ǫ(x, ω), γ Γ (1) for almost all x = (x 1, x 2 ) R 2. A function with the property (1) is called Γ-periodic. The basic period cell of the lattice Γ is Ω = ( π, π] 2 and Ω = 4π 2 denotes the size of the basic cell. We assume that ǫ is in L (R 2 ) for each ω [ω 0, ω 1 ] and that ǫ is invertible, positive and bounded 0 < c 0 ǫ(x, ω) c 1, (2) for almost all x R 2 and ω [ω 0, ω 1 ]. Let = R 2 /Γ denote the torus in two dimensions and, the scalar product in R 2. Denote by C ( ) the space of smooth complex-valued Γ-periodic functions. The space L 2 ( ) is defined to be the completion of C ( ) in the L 2 -norm and H 1 ( ) is the completion of C ( ) in the H 1 -norm. Let denote the gradient with respect to x R 2 and let denote the divergence. The transverse magnetic polarized waves (0, 0, E 3, H 1, H 2, 0) can be written in terms of the electric field E alone. The equation v = ω 2 ǫ(x, ω)v (3) models a time-harmonic monochromatic electromagnetic wave with frequency ω and electric polarization E(x) = (0, 0, v(x)). A Bloch solution of (3) is a solution on the form v(x) = e i k,x u(x), (4) where u is a Γ-periodic function and k C 2 is the quasimomentum (or the Floquet-Bloch wave vector) [11]. Since (e i k,x u(x)) = e i k,x ( + ik)u(x) the Bloch solutions of (3) are Γ-periodic solutions of ( + ik) ( + ik)u(x) = ω 2 ǫ(x, ω)u. (5) 2
3 The frequency ω as a multi-valued function of the quasimomentum k is called the dispersion relation and the graph of the dispersion relation defines the Bloch variety [11]. Definition 1 The complex Bloch variety is the set B = {(k, ω) C 2 C the problem (3)has a non zerobloch solution} (6) In the case of smooth coefficients, analytic properties of the Bloch variety B are discussed in Kuchment [11]. The Bloch solutions of (3) can be determined from spectrum problems. In the non-dispersive case ǫ = ǫ(x), Bloch solutions with k R 2 are given by eigenvectors of a selfadjoint operator. In the dispersive case ǫ = ǫ(x, ω), the spectrum problem is non-linear, which complicates the analysis. We will in the dispersive case consider solutions (k, ω) C 2 R and use the real frequency ω as a parameter. The two cases are discussed below. 2.1 The frequency independent case The symbol L 2 (R 2, ǫdx) is used for the weighted L 2 -space with the scalar product (u, v) ǫ = u(x)v(x)ǫ(x)dx. (7) R 2 We introduce in L 2 (R 2, ǫdx) the quadratic form q(v) = v 2 dx, v H 1 (R 2, ǫdx). (8) R 2 The quadratic form q is semi-bounded from below and closed. The form q defines in L 2 (R 2, ǫdx) the associated selfadjoint operator A(x, )v := ǫ 1 (x) v. (9) The spectrum of elliptic operators with periodic coefficient has been studied in detail by many authors including Odeh and Keller [14], Reed and Simon [15], and Figotin and Kuchment [4]. The spectrum σ(a) has a band structure and the existence of band-gaps for certain geometrical structures has been proved [4, 5, 6]. From the Floquet-Bloch theory [15, 11] follows the spectrum of A in the infinite periodic medium can be obtained from the k - dependent family of Γ-periodic solutions (5). We give below the well-known construction of the band structure. Define a family of quadratic forms depending on k R 2 by q k (u) = ( + ik)u 2 dx, u H 1 (, ǫdx), (10) where k is the quasimomentum. The quadratic form q k defines for every real k a unique selfadjoint operator A(k) in L 2 (, ǫdx) that we write as A(k)u := A(x, + ik)u = ǫ 1 (x)( + ik) ( + ik)u. (11) 3
4 The spectral parameter ω 2 represent the square of the time frequency ω of the wave. Since each operator A(k) is selfadjoint in L 2 (, ǫdx) and has a compact resolvent, the spectrum of A(k) consists of eigenvalues 0 ω 2 1(k) ω 2 2(k)... The main spectrum result is that the spectrum of A, that is defined for functions on R 2, is the union σ(a) = S n, S n = [min k ωn(k), 2 maxωn(k)]. 2 (12) k The intervals S n are called the stability zones of the operator A and the complement of σ(a) is called the instability zone. The successive segments S n and S n+1 may overlap, but we have a gap in the spectrum σ(a) if the intersection S n S n+1 is empty for some n. That is, if ω 2 belongs to the instability zone for all k R 2. The spectrum σ(a) can also be described in terms of the Bloch variety (6). The real Bloch variety is the intersection B real = B R 3 (13) and the spectrum of A is the projection of the real Bloch variety onto the ω- axis. 3 The frequency dependent case We assumed above that the permittivity ǫ is independent of the frequency ω. The spectrum σ(a) coincides in the frequency independent case with the range of the dispersion relation ω 2 (k). The permittivity can be constant in a frequency interval, but only free space is independent of the frequency. If we exclude free space, the Kramers-Kronig relations [2, 7] implies that the imaginary part of ǫ(ω) cannot be zero for all ω. We focus on the dispersive case ǫ = ǫ(x, ω), where ǫ is real-valued in the frequency interval [ω 0, ω 1 ]. That is, we assume that there is no absorption line in the frequency interval under consideration. This is a reasonable model for many dialectic materials (insulators) for a range of frequencies in the micro-wave regime and in the optical regime [7]. The permittivity has certain holomorphic properties as a function of ω, which make it possible to analyze ω(k) with the theory for holomorphic operatorvalued functions. A more tempting alternative is to study the multi-valued vector function k(ω), which independent of ǫ leads to a quadratic polynomial operator pencil. Below, we study the quadratic pencil in the case of a real-valued permittivity function. 3.1 The quadratic operator pencil We presented in section 2.1 results for the shifted operator (11) in the nondispersive case, where the spectral parameter is related to the time frequency 4
5 ω. We consider below the dispersive case and relate the spectrum parameter to the quasimomentum k. We reduce the spectrum problem in k to a problem in the complex amplitude of the vector k. Let k = λθ, where λ C and θ on the unit sphere S 1. We write in some cases the complex spectrum parameter λ on the form λ = λ r + λ i i, (14) where λ r and λ i are real numbers. Define a family of spectrum problems as ( + iλθ) ( + iλθ)u(x) = ω 2 ǫ(x, ω)u, (15) where u and ǫ are Γ periodic functions. This spectrum problem in (λ, u) corresponds to the problem A(k)u = ω 2 u in (11), but we consider k(ω) and not ω(k). The family of spectrum problems (15) is quadratic in the spectral parameter λ, with parameter ω. We write below the classical problem (15) on a weak form. Define the sesquilinear forms a 0 : H 1 ( ) H 1 ( ) C a 1 : L 2 ( ) H 1 ( ) C a 2 : L 2 ( ) L 2 ( ) C (16) where a 0 (u, v) = u, v ω 2 ǫuvdx, T 2 a 1 (u, v) = 2i u θ, v dx, T 2 a 2 (u, v) = uvdx. (17) The weak solutions of (15) are defined as eigenvectors u of H 1 ( )\{0} and eigenvalues λ C which satisfy for all v H 1 ( ). λ 2 a 2 (u, v) + λa 1 (u, v) + a 0 (u, v) = 0, (18) Lemma 1 The sesquilinear forms a 0, a 1 and a 2 are bounded for all elements u, v H 1 ( ). Proof: The sesquilinear forms a 0 and a 2 are trivially bounded and the boundedness of a 1 follows from the inequality ab δ a 2 /2 + b 2 /(2δ), δ > 0. Let (, ) H1 ( ) denote the scalar product in H 1 ( ). The spectrum problem (18) can be written λ 2 (A 2 u, v) H 1 ( ) + λ(a 1 u, v) H 1 ( ) + (A 0 u, v) H 1 ( ) = 0, (19) 5
6 where the operators A 0, A 1 and A 2 are defined by the bounded sesquilinear forms (A 0 u, v) H 1 ( ) = a 0 (u, v), (A 1 u, v) H 1 ( ) = a 1 (u, v), (A 2 u, v) H 1 ( ) = a 2 (u, v). (20) According to Riesz theorem are these operators on H 1 ( ) unique, linear and bounded [8]. The equality (19) that holds for all v H 1 ( ) is equivalent to the operator equation A 0 u + λa 1 u + λ 2 A 2 u = 0 (21) in H 1 ( ). Let L(H 1 ( )) denote the set of all bounded linear operators on H 1 ( ). We introduce the operator pencil K : C L(H 1 ( )) K(λ) = A 0 + λa 1 + λ 2 A 2, λ C. (22) The quadratic eigenvalue problem is then: Find λ C and a non-zero u H 1 ( ) such that K(λ)u = 0. (23) The adjoint operator pencil is K (λ) = A 0 + λa 1 + λ2 A 2, λ C (24) and the pencil K(λ) is said to be selfadjoint if A 0, A 1 and A 2 are selfadjoint. The lemma below shows that the pencil K(λ) is selfadjoint and hence closed. Lemma 2 The operators A 0, A 1 and A 2 are self-adjoint Proof: The lemma follows from the equalities (A 0 u, v) H1 ( ) = a 0 (u, v) = (A 0 v, u) H1 ( ) = (v, A 0 u) H 1 ( ) = (A 0 u, v) H 1 ( ) (A 1 u, v) H 1 ( ) = a 1 (u, v) = (A 1 v, u) H 1 ( ) = (v, A 1 u) H 1 ( ) = (A 1u, v) H 1 ( ), (A 2 u, v) H 1 ( ) = a 2 (u, v) = (A 2 v, u) H 1 ( ) = (v, A 2 u) H 1 ( ) = (A 2 u, v) H 1 ( ). 3.2 Spectral properties We consider in this section spectral properties of the operator pencil (22). Fredholm theory is used to prove that the spectrum of K is discrete, where the eigenvalues are symmetrically placed with respect to the real and imaginary axis. Moreover, we prove that the pencil K(λ) cannot be reduced to a monic bounded pencil, that is, to a pencil on the form K(λ) = B 0 + λb 1 + λ 2 I, λ C. (25) Nonmonic pencils has a more complex spectrum structure [12]. Lemma 3 The operator A 2 has an unbounded inverse. 6
7 Proof: Since A 2 > 0 we have [1] A 2 u 2 H 1 ( ) A 2 H 1 ( )(A 2 u, u) H 1 ( ) = c u 2 dx. (26) The operator A 2 has a bounded inverse if and only if there exist a positive constant c such that A 2 u H 1 ( ) c u H 1 ( ) for all u H 1 ( ). The inequality (26) implies that such constant c cannot exist (take for instance u(x) = e i n,x ), but A 2 has an inverse since the null space N(A 2 ) is trivial. A bounded operator T is Fredholm if the range R(T) is closed and if the null space N(T), and the complement of the range R(T) are finite dimensional. The Fredholm index of a Fredholm operator T is defined as indt = nul T def T, (27) where the nullity nul T is the dimension of the null space and the deficiency def T is the dimension of the complement to the range R(T). The resolvent set ρ is the set of all λ C, such that the operator T(λ) is boundedly invertible and the spectrum is defined by the complement σ(t) = C\ρ(T). (28) The essential spectrum σ ess (T) of T is the set of all λ σ(t) such that T(λ) is not a Fredholm operator. The complement of the essential spectrum σ p (T) = σ(t)\σ ess (T) (29) is called the discrete spectrum. The discrete spectrum consists of eigenvalues of finite geometrical multiplicity, nul T <. We shall show that the spectrum of K(λ) is discrete. The proof consists of two steps. We prove that the operator pencil K(λ) is a sum of compact operators and of a Fredholm operator. Lemma 4 The operators A 1 and A 2 are compact in H 1 ( ). Proof: We show that A 1 is compact. Let {u n } H 1 ( ) be a given bounded sequence. Then since H 1 ( ) is compactly embedded into L 2 ( ) [17], the sequence {u n } has a convergent subsequence {u m } in L 2 ( ). Let u m and u m denote two element in the subsequence {u m }. From the boundedness of the Hermitian form a 1 follows A 1 u m A 1 u m 2 H 1 ( ) = a 1(u m u m, A 1 u m A 1 u m ) C A 1 u m A 1 u m H 1 ( ) u m u m L 2 ( ). (30) Then since {u m } is a Cauchy sequence in L 2 ( ), the sequence {A 1 u m } is a Cauchy sequence in H 1 ( ), which converge since H 1 ( ) is complete. Using the inequality u L 2 ( ) C u H 1 ( ), the compactness of A 2 is proved in the same way. 7
8 Theorem 1 The operator pencil K(λ) is a Fredholm operator with index zero. Proof: The operator A 0 in (22) can be written on the form where A (1) 0 and A (2) (A 0 u, v) H 1 ( ) = (A (1) 0 u, v) H 1 ( ) + (A (2) 0 u, v) H 1 ( ), (31) (A (1) 0 u, v) H 1 ( ) = ω 2 0 are defined by the bounded sesquilinear forms ǫuvdx, (A (2) 0 u, v) H 1 ( ) = u, v dx. (32) The operator A (1) 0 is selfadjoint and compact since ǫ L ( ) is real-valued and H 1 ( ) is compactly embedded into L 2 ( ). Represent u, v H 1 ( ) with its Fourier series: u(x) = n Z 2 û n e i n,x, v(x) = m Z 2 ˆv m e i m,x. (33) Since the gradient of u and of v are in L 2 ( ) the Plancherel s identity [19] gives (A (2) 0 u, v) H 1 ( ) = dx n Z 2 inû n e i n,x, = Ω n Z 2 n 2 û nˆv n, m Z 2 imˆv me i m,x (34) which implies that the dimension of the null space is one. The dimension of the orthogonal complement of the range is also one, since R(T) = N(T) for self-adjoint operators. That is, A (2) 0 is a Fredholm operator with index (27) zero. The self-adjoint operator pencil K(λ) is a sum of A (2) 0 and compact operators, which is a Fredholm operator with the index of A (2) 0 [8] Hamiltonian symmetry We show in this section that the eigenvalues λ of (18) are symmetrically placed with respect to the real and imaginary axes. Lemma 2 implies that the spectrum of K(λ) is symmetric with respect to R, since K(λ) is invertible if and only if K (λ) is invertible. Lemma 5 The operators A 0, A 1 and A 2 have the properties A 0 u = A 0 u, A 1 u = A 1 u and A 2 u = A 2 u. Proof: The lemma follows from the equalities (A 0 u, v) H1 ( ) = (A 0 u, v) H1 ( ) = a 0(u, v) = (A 0 u, v) H1 ( ), (A 1 u, v) H1 ( ) = (A 1 u, v) H1 ( ) = a 1(u, v) = (A 1 u, v) H1 ( ), (A 2 u, v) H 1 ( ) = (A 2 u, v) H 1 ( ) = a 2(u, v) = (A 2 u, v) H 1 ( ). The theorem below states the Hamiltonian structure of the spectrum σ(k(λ)). 8
9 Theorem 2 The spectrum of the quadratic operator pencil K(λ) consists of eigenvalues of finite multiplicity, λ C, that have the Hamiltonian structure (λ, λ, λ, λ). Proof: Theorem (1) states that K(λ) is a Fredholm operator, which implies that the spectrum σ(k(λ)) consists of eigenvalues of finite multiplicity. Lemma 5 gives the relation K(λ)u = A 0 u + λa 1 u + λ 2 A 2 u = A 0 u λa 1 u + λ 2 A 2 u = K( λ)u, (35) which implies that λ is an eigenvalue whenever λ is an eigenvalue. The two relations K ( λ) = A 0 λa 1 + λ 2 A 2 = K( λ), (36) follows from lemma 2. K (λ) = A 0 + λa 1 + λ 2 A 2 = K(λ), (37) 3.3 Band gaps and spectral estimates The complex Bloch variety of the pencil K(λ) is the set B(K) = {(k, ω) C 2 R K(λ)u = 0 has a non zerosolution}, (38) where the ω-dependence in K(λ) is a real parameter. The real Bloch variety of the pencil K(λ) is the intersection B real (K) = B(K) R 3 (39) We state below the condition for a band-gap in terms the Bloch variety. Definition 2 The frequency ω is in a band gap if the projection of the real Bloch variety onto the ω-axis is empty. Solutions on the form (k, ω) R 2 R exist if the frequency ω is not in a band gap. That is, propagating Bloch waves exist for the given frequency ω and we say that ω belongs to the stability zone. Notice that we let k be real in the frequency dependent case and defined a band gap in (12) as an empty intersection of two successive bands [min ωn 2 (k), max k k ω 2 n (k)] [min k ω2 n+1 (k), max k ωn+1 2 (k)] = {} (40) or equally in terms of the projection of the real Bloch variety (13) onto the ω-axis. We define the dual lattice Γ = {q R 2 γ, q 2πZ, γ Γ} (41) where Γ = 2πZ 2. The basic cell (the Brillouin zone), of the dual lattice Γ = Z 2 is ( Ω = 1 2 2], 1 2. (42) 9
10 The real part of all eigenvalues λ to (15) is Ω - periodic, since the Bloch solutions (4) are periodic in the real part of k = θλ. We introduce the set B Ω = {(λ r θ, λ i θ, ω) Ω R 2 R K(λ)u = 0 has a non zerosolution}, (43) which is a subset of the Bloch variety (38). This set corresponds to solutions with the real part of the wave vector k = θλ in the Brillouin zone Ω and a real time frequency ω. The frequency ω is in a band-gap if λ i 0 for all eigenvalues of K(λ), with θλ r Ω. A gap in the spectrum of K(λ) is also called a polarization gap since we only consider the polarization E(x) = (0, 0, u(x)). When ω = 0, the lemma below shows that the eigenvalues can be calculated explicitly. We use the lemma in the numerical calculation in section 5. Lemma 6 Let ω = 0 and θ = (1, 0) or θ = (0, 1). Then all quasimomentum k C 2 can be written on the form k = Z 2 + iz 2. (44) Proof: Represent u, v H 1 ( ) with its Fourier series: u(x) = n Z 2 û n e i n,x, v(x) = m Z 2 ˆv m e i m,x. (45) The equation (K(λ)u, v) H 1 ( ) = 0 gives n Z 2 (n + λθ) 2 û nˆv n = 0. (46) The coefficients are zero whenever û n 0. This condition (n + θλ) 2 = 0 can be written as { n 2 + 2λ r θ, n + λ 2 r λ2 i = 0, λ i θ, n + λ r λ i = 0. (47) Assume λ r = 0. Then follows λ i = ± n and θ, n = 0, which implies λ i Z when θ = (1, 0) or θ = (0, 1). If λ i = 0 and λ r θ Ω can the first equation be estimated as n 2 + 2λ r θ, n + λ 2 r λ 2 i n 2 n + λ 2 r. (48) Then follows λ r = 0, since n 2 n 0. If λ r 0 and λ i 0 can the eigenvalues be written on the form λ r = θ, n, λ i = ± n 2 θ, n 2. (49) The Hamiltoninan structure implies that λ r = θ, n, λ i = ± n 2 θ, n 2 (50) 10
11 also are eigenvalues. For θ = (1, 0) or θ = (0, 1) are the quasimomentum k then integer points in C 2. The frequency ω = 0 is an eigenvalue of A 0 = A 0 (ω) and the corresponding eigenvector is constant. An eigenvalue ω of the operator-valued function A 0 (ω) corresponds to a solution of the classical eigenvalue problem u = ω 2 ǫ(x, ω)u, (51) where u and ǫ are Γ-periodic. Assume that (u, ω) = (ũ, ω) is a solution of A 0 (ω)u = 0. The pencil K(λ) can then be written on the form K(λ)ũ = λ(a 1 ũ + λa 2 ũ). (52) Zero is always an eigenvalue of the pencil and a nonzero eigenvalue λ is a solution of A 1 ũ + λa 2 ũ = 0, (53) where ũ N(A 0 ( ω)). All eigenvalues of (52) are on the form λ = 2 θ, n since an expansion in Fourier series gives n Z 2 (2 θ, n + λ) ˆũ n 2 = 0. (54) The eigenvalues λ are all real and the number of eigenvalues at ω = ω depends on the dimension of the null space N(A 0 ( ω)) Estimates of the numerical range Below we present estimates of the spectrum that are tight when ω is close to zero. We use the numerical range of an operator-valued function that is defined below. Definition 3 The numerical range of K(λ) is defined as W(K) = {λ C ; (K(λ)u, u) H1 ( ) = 0 for some u H 1 ( )\{0}} (55) The spectrum of K(λ) is contained in the closure of W(K) [12]. The values in the numerical range of K(λ) are symmetrically placed with respect to the real and imaginary axes. This Hamiltonian symmetry is proved in the same way as the Hamiltonian structure of the spectrum in section The space of all periodic functions in H 1 ( ) with meanvalue zero is denoted V = {u H 1 ( ); udx = 0}. (56) Lemma 7 If λ C with λ r 0 and λ i 0. Then u V is a necessary condition for a non-trivial solution of (K(λ)u, u) H 1 ( ) = 0. (57) 11
12 Proof: Represent u H 1 ( ) with its Fourier series: u(x) = n Z 2 û n e i n,x. (58) When substituting the Fourier series into the eigenvalue equation (K(λ)u, u) H1 ( ) = 0, the equation for the Fourier coefficients becomes Ω n Z 2 (n + λθ) 2 û n 2 = ω 2 ǫ(x, ω) n Z 2 û n e i n,x 2 dx. (59) Since ǫ is real, the imaginary part of the right side of the equation is zero and we obtain the condition n Z 2 ( n, θ λ i + λ r λ i ) û n 2 = 0. (60) From the assumption λ r 0, λ i 0 and condition (60) follows û 0 = 0. We give below an upper bound and a lower bound on the real part of λ 2 that is independent of the geometrical structure. Lemma 8 Let λ r θ Ω. Then all λ W(K) have the property λ 2 r λ2 i ω2 ess sup ǫ(, ω). (61) Proof: Represent u H 1 ( ) with its Fourier series: u(x) = n Z 2 û n e i n,x. (62) A complex λ is in W(K) if Ω n Z 2 (n + λθ) 2 û n 2 = ω 2 ǫ(x, ω) u 2 dx. (63) The real part of the terms in the Fourier expansion is R{(n + λθ) 2 } = n, n + 2 n, λ r θ + λ 2 r λ 2 i. (64) Since λ r θ Ω the n-dependent terms can be estimated as n, n + 2 n, λ r θ n 2 n 0, (65) which implies } R { Ω n Z2 (n + λθ) 2 û n 2 ω 2 ǫ(x, ω) u 2 dx Ω [(λ 2 r λ 2 i) ω 2 ess sup ǫ(, ω)] u L 2 ( ). (66) 12
13 The left hand side of the inequality is zero when λ is in W(K), but the right hand side is according to the inequality above negative if (61) is not satisfied. From lemma 8, we obtain the upper bound λ r λ 2 i + ω2 ess supǫ(, ω), (67) but we have no restriction on λ i. We also derive lower bounds on the real part of λ 2. Lemma 9 Let λ r θ Ω. Then all λ W(K) with λ r 0 and λ i 0 have the property λ 2 r λ 2 i ω 2 ess inf ǫ(, ω) a(ω), (68) where a(ω) is non-negative. Proof: The terms in the Fourier expansion (59) can be estimated as R(n + λθ) 2 n 2 + n + λ 2 r λ 2 i (69) which implies that the Fourier series in (59) satisfies the inequality R(n + λθ) 2 û n 2 (λ 2 r λ2 i ) û n 2 + ( n 2 + n ) û n 2. (70) n Z 2 n Z 2 n Z 2 From the inequality above and the positiveness of the permittivity ǫ follows } R { Ω n Z2 (n + λθ) 2 û n 2 ω 2 ǫ(x, ω) u 2 dx Ω [(λ 2 r λ 2 i ) ω 2 ess inf ǫ(, ω)] û n 2 + Ω (71) ( n 2 + n ) û n 2 n Z 2 n Z 2 The left side of the inequality is zero when λ is in W(K), but the right side is according to the inequality above negative if (68) is not satisfied, with a(ω) = Ω n Z 2 ( n 2 + n ) û n 2. (72) The estimate (68), with a = 0, trivially holds for a constant eigenvector u. From lemma (6) follows that a constant function is an eigenvector when ω = 0 and λ = 0 is the only real solution with λ r θ Ω. The eigenspace of K(0) at ω = 0 consists only of constant functions as (K(0)u, v) H1 ( ) = u, v dx (73) has a one-dimensional null space. The value on a(ω) depends on the geometry, but we expect b to be close to zero for low frequencies ω. 13
14 We obtain tighter bounds on the real part of λ 2 if we include geometrical information. Assume that the material ǫ consists of N components N ǫ(x, ω) = ǫ j (ω)χ j (x) (74) j=1 where ǫ j is positive and χ j is the characteristic function of component j. We give below an upper bound on the real part of λ 2 in terms of the volume fraction of the N components ǫ j. Lemma 10 Assume that the material is a N component composite (74). An element λ W(λ) with λ r θ Ω satisfies the inequality N λ 2 r λ 2 i ω 2 f j ǫ j, (75) j=1 where f j is the volume fraction of component ǫ j. Proof: The integral of ǫ u 2 over is N N ǫ j χ j u 2 dx = ǫ j u Ω 2 dx. (76) j Ω j=1 The real part of (n + λθ) 2 satisfies R{(n + λθ) 2 } λ 2 r λ 2 i which implies that (59) can be estimated as } R { Ω n Z2 (n + λθ) 2 û n 2 ω 2 ǫ(x, ω) u 2 dx j=1 Ω [(λ 2 r λ2 i ) N ω2 ǫ j u Ω 2 dx. j j=1 (77) The integral over Ω j satisfies u 2 2 dx = û n e i n,x dx Ω j û n 2. (78) Ω j Ω j n Z 2 n Z 2 From (K(λu, u) = 0 and the inequalities (77) and (78) follows N 0 Ω (λ2 r λ2 i ) Ω j ω2 Ω ǫ j. (79) The bound (75) is optimal in the sense that we can construct a geometry that satisfies the equality. Let ǫ = ǫ(x 1, ω) denote the permittivity of a material that is constant in x 2. Take θ = (0, 1) and u = 1. Then j=1 π 2π(λ 2 r λ2 i ) = ω2 ǫ(x 1, ω)dx 1, (80) which implies (75) with equality in the N component case. π 14
15 4 Numerical treatment In order to compute approximate eigenvalues λ C that satisfy a 0 (u, v) + λa 1 (u, v) + λ 2 a 2 (u, v) = 0 (81) for some u H 1 ( ) \ {0} and all v H 1 ( ), we discretize the problem by a Galerkin ansatz. We use a finite element method with quadratic elements to construct N real-valued basis functions φ 1,...,φ N that span an N-dimensional subspace V H 1 ( ). The periodic boundary conditions inferred by the topology of were incorporated directly in the finite element space. After assembling the finite element matrices we obtain the matrix pencil M(λ) : C C N N M(λ) = M 0 + λm 1 + λ 2 M 2, λ C. (82) The matrices M 0, M 1 and M 2 are the descrete versions of (81), where the entries are (M 0 ) mn = φ n, φ m ω 2 ǫ(ω)φ n φ m dx T 2 (M 1 ) mn = i φ n θ, φ m dx i φ n, θ φ m dx (83) T 2 (M 2 ) mn = φ n φ m dx for all i, j = 1,...,N. The matrix M 0 is real and symmetric, the matrix M 1 is purely imaginary and Hermitian, and the matrix M 2 is real, symmetric and positive definite. All three matrices are typically large and sparse. Since we prefer to work with real matrices only, we define the matrices B 0 = M 0, B 1 = im 1, B 2 = M 2, and the matrix pencil G : C C N N Obviously, G(µ) = B 0 + µb 1 + µ 2 B 2, µ C. (84) G(µ) = M( iµ), (85) i.e., the eigenvalues µ of G are related to the eigenvalues λ of M by λ = iµ. The pencil G is called a gyroscopic pencil, since B 0 is real and symmetric, B 1 is real and skew-symmetric, and B 2 is real, symmetric and positive definite. Pencils of this form naturally occur in the analysis of gyroscopic systems [18]. There is an on-going research effort in devising efficient numerical methods for these pencils [13, 10]. In this work we used a simple linearization of the quadratic eigenvalue problem in combination with a Krylov space method [16] to compute a few number of eigenvalues with the lowest modulus. There are several ways to linearize a quadratic eigenvalue problems [18]. The one we used was proposed in [13] it leads to an eigenvalue problem of the form Wx = µx, (86) 15
16 Case I Case II Gap open close open close Table 1: Case I: Constant permittivity from Dobson ǫ 1 = 1 and ǫ 2 = 8.9. Case II: The frequency dependent model (87). where W R 2N 2N is a Hamiltonian matrix, thus preserving the structure of the quadratic spectrum. The MATLAB routine eigs can be applied to this problem to obtain the sought approximations. 5 Numerical experiments We use the numerical algorithm in section 4 to compute band gaps of a rectangular array of dielectric cylinders in air. The permittivity of air is ǫ 1 = 1 for all frequencies under consideration. The permittivity of the dielectric cylinders is generally a complicated function of the frequency. It is common in practice to use least-squares minimization to determine a rational function approximation from frequency-domain data. The choice of material model is not important for the algorithm, since the frequency dependence only is a parameter. We will consider the simple Lorentz model ǫ 2 (ω) = a + b c ω 2, (87) where a = 4.5, b = 7.12, and c = 2. The constants in (87) are chosen such that the mean value (ǫ(0) + ǫ(0.8))/2 is 8.9. We assume that this model is valid for all frequencies in the range ω [0, 0.8]. This assumption make it possible to compare the frequency dependent model with the calculations for the frequency independent case in Dobson [3]. Dobson uses the cylinder radii r = 0.378d, where d is the lattice constant 2π. One can show that the band structure can be calculated reducing the angles to 0 θ π/4, which is called the irreducible Brillouin zone [9]. Figure 1 shows the eigenvalues λ(ω, θ) with smallest imaginary part as θ is varied over the irreducible Brillouin zone and ω [0, 0.8]. The band gap frequencies in Table 1 match well with Dobson s results [3]. A frequency ω is in a band gap, in the sense of definition 2, when all eigenvalues have a non-zero imaginary part. The band gaps in Table 1 corresponds in Figure 1 to the tubes. The frequency dependent model gives a monotone permittivity function with ǫ(0) = 8.06 and ǫ(0.8) = Observe that the difference between the mean value 8.9 and the maximum and minimum only is 9.4 percent. Figure 1 shows the imaginary part of λ(ω, θ) and Table 1 gives the band gap frequencies for the frequency dependent model (87). 16
17 Im (ka/π) 0 Im (ka/π) θ/π ω/2π θ/π ω/2π Figure 1: Imaginary parts of λ. Left: Material parameters ǫ 1 = 1 and ǫ 2 = 8.9. Right: ǫ 1 = 1 and ǫ 2 according to the Lorentz model (87). 6 Acknowledgments The authors express their gratitude to Thomas Gauss, Melanie Reimers, and Arne Schneck from research training group 1294 for helpful discussions and acknowledges the support of the German Research Foundation (RTG 1294). References [1] M. S. Birman and M. Z. Solomjak, Spectral Theory of Self-Adjoint Operators in Hilbert Space, D. Reidel Publishing, Dordrecht, [2] M. Cessenat, Mathematical Methods in Electromagnetism, Series on Advances in Mathematics for Applied Sciences Vol. 41, World Scientific Publisher, Singapore, [3] D. C. Dobson, An efficient band structure calculations in 2d photonic crystals, J. Comput. Phys., 149 (1999), pp [4] A. Figotin and P. Kuchment, Band-gap structure of spectra of periodic dielectric and acoustic media. i scalar model, SIAM J. Appl. Math., 56 (1996), pp [5], Band-gap structure of spectra of periodic dielectric and acoustic media. ii two-dimensional photonic crystals, SIAM J. Appl. Math., 56 (1996), pp [6], Spectral properties of classical waves in high-contrast periodic media, SIAM J. Appl. Math., 58 (1998), pp [7] J. D. Jackson, Classical Electrodynamics, John Wiley & Sons, New York, third ed.,
18 [8] T. Kato, Perturbation Theory for Linear Operators, Springer-Verlag, Berlin, [9] C. Kittel, Introduction to Solid State Physics, John Wiley & Sons, New York, 6 ed., [10] D. Kressner, Numerical Methods for General and Structured Eigenvalue Problems, Springer-Verlag, Berlin, [11] P. Kuchment, Floquet Theory for Partial Differential Equations, Birkhäuser, Basel, [12] A. S. Markus, Introduction to the Spectral Theory of Polynomial Operator Pencils, Transl. Math. Monogr. 71, AMS, Providence, [13] V. Mehrmann and D. Watkins, Structure-preserving methods for computing eigenpairs of large sparse skew-hamiltonian/hamiltonian pencils., SIAM J. Sci. Comput., 22 (2001), pp [14] F. Odeh and J. B. Keller, Partial differential equations with periodic coefficients and Bloch waves, J. Math. Phys., 5 (1964), pp [15] M. Reed and B. Simon, Methods of modern mathematical physics, vol. IV: Analysis of Operators, Academic Press, New York, [16] D. C. Sorensen, Implicit application of polynomial filters in a k-step arnoldi method, SIAM J. Matrix Anal. Appl., 13 (1992), pp [17] M. E. Taylor, Pseudodifferential Operators, Princeton University Press, Princeton, [18] F. Tisseur and K. Meerbergen, The quadratic eigenvalue problem., SIAM Rev., 43 (2001), pp [19] K. Yosida, Functional Analysis, Springer-Verlag, Berlin, DEPARTMENT OF MATHEMATICS, KARLSRUHE UNIVERSITY, GERMANY address: christian.engstroem@math.uni-karlsruhe.de, markus.richter@math.uni-karlsruhe.de 18
arxiv: v2 [math-ph] 1 Jul 2008
ON THE SPECTRUM OF AN OPERATOR PENCIL WITH APPLICATIONS TO WAVE PROPAGATION IN PERIODIC AND FREQUENCY DEPENDENT MATERIALS CHRISTIAN ENGSTRÖM AND MARKUS RICHTER Abstract. We study wave propagation in periodic
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 informationPerturbation Theory for Self-Adjoint Operators in Krein spaces
Perturbation Theory for Self-Adjoint Operators in Krein spaces Carsten Trunk Institut für Mathematik, Technische Universität Ilmenau, Postfach 10 05 65, 98684 Ilmenau, Germany E-mail: carsten.trunk@tu-ilmenau.de
More informationSELF-ADJOINTNESS OF SCHRÖDINGER-TYPE OPERATORS WITH SINGULAR POTENTIALS ON MANIFOLDS OF BOUNDED GEOMETRY
Electronic Journal of Differential Equations, Vol. 2003(2003), No.??, pp. 1 8. ISSN: 1072-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) SELF-ADJOINTNESS
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 informationTHE FORM SUM AND THE FRIEDRICHS EXTENSION OF SCHRÖDINGER-TYPE OPERATORS ON RIEMANNIAN MANIFOLDS
THE FORM SUM AND THE FRIEDRICHS EXTENSION OF SCHRÖDINGER-TYPE OPERATORS ON RIEMANNIAN MANIFOLDS OGNJEN MILATOVIC Abstract. We consider H V = M +V, where (M, g) is a Riemannian manifold (not necessarily
More informationA computer-assisted Band-Gap Proof for 3D Photonic Crystals
A computer-assisted Band-Gap Proof for 3D Photonic Crystals Henning Behnke Vu Hoang 2 Michael Plum 2 ) TU Clausthal, Institute for Mathematics, 38678 Clausthal, Germany 2) Univ. of Karlsruhe (TH), Faculty
More informationON THE EXISTENCE OF TRANSMISSION EIGENVALUES. Andreas Kirsch1
Manuscript submitted to AIMS Journals Volume 3, Number 2, May 29 Website: http://aimsciences.org pp. 1 XX ON THE EXISTENCE OF TRANSMISSION EIGENVALUES Andreas Kirsch1 University of Karlsruhe epartment
More informationTitle: Localized self-adjointness of Schrödinger-type operators on Riemannian manifolds. Proposed running head: Schrödinger-type operators on
Title: Localized self-adjointness of Schrödinger-type operators on Riemannian manifolds. Proposed running head: Schrödinger-type operators on manifolds. Author: Ognjen Milatovic Department Address: Department
More informationDiscreteness of Transmission Eigenvalues via Upper Triangular Compact Operators
Discreteness of Transmission Eigenvalues via Upper Triangular Compact Operators John Sylvester Department of Mathematics University of Washington Seattle, Washington 98195 U.S.A. June 3, 2011 This research
More informationQUASI-UNIFORMLY POSITIVE OPERATORS IN KREIN SPACE. Denitizable operators in Krein spaces have spectral properties similar to those
QUASI-UNIFORMLY POSITIVE OPERATORS IN KREIN SPACE BRANKO CURGUS and BRANKO NAJMAN Denitizable operators in Krein spaces have spectral properties similar to those of selfadjoint operators in Hilbert spaces.
More informationA perturbation-based reduced basis method for parametric eigenvalue problems
A perturbation-based reduced basis method for parametric eigenvalue problems Eric Cancès 1, Virginie Ehrlacher 1, David Gontier 2, Damiano Lombardi 3 and Antoine Levitt 1 1. CERMICS, Ecole des Ponts ParisTech
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 quadratic eigenvalue problem (QEP) is to find scalars λ and nonzero vectors u satisfying
I.2 Quadratic Eigenvalue Problems 1 Introduction The quadratic eigenvalue problem QEP is to find scalars λ and nonzero vectors u satisfying where Qλx = 0, 1.1 Qλ = λ 2 M + λd + K, M, D and K are given
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 informationBASIC MATRIX PERTURBATION THEORY
BASIC MATRIX PERTURBATION THEORY BENJAMIN TEXIER Abstract. In this expository note, we give the proofs of several results in finitedimensional matrix perturbation theory: continuity of the spectrum, regularity
More informationE. YABLONOVITCH photonic crystals by using level set methods
Appl. Phys. B 81, 235 244 (2005) DOI: 10.1007/s00340-005-1877-3 Applied Physics B Lasers and Optics C.Y. KAO 1, Maximizing band gaps in two-dimensional S. OSHER 2 E. YABLONOVITCH photonic crystals by using
More information1 Math 241A-B Homework Problem List for F2015 and W2016
1 Math 241A-B Homework Problem List for F2015 W2016 1.1 Homework 1. Due Wednesday, October 7, 2015 Notation 1.1 Let U be any set, g be a positive function on U, Y be a normed space. For any f : U Y let
More informationRiesz bases of Floquet modes in semi-infinite periodic waveguides and implications
Riesz bases of Floquet modes in semi-infinite periodic waveguides and implications Thorsten Hohage joint work with Sofiane Soussi Institut für Numerische und Angewandte Mathematik Georg-August-Universität
More informationThe Residual Spectrum and the Continuous Spectrum of Upper Triangular Operator Matrices
Filomat 28:1 (2014, 65 71 DOI 10.2298/FIL1401065H Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat The Residual Spectrum and the
More informationApplied Mathematics Letters
Applied Mathematics Letters 25 (2012) 1926 1930 Contents lists available at SciVerse ScienceDirect Applied Mathematics Letters journal homepage: wwwelseviercom/locate/aml Bragg structure and the first
More informationON THE SPECTRUM OF NARROW NEUMANN WAVEGUIDE WITH PERIODICALLY DISTRIBUTED δ TRAPS
ON THE SPECTRUM OF NARROW NEUMANN WAVEGUIDE WITH PERIODICALLY DISTRIBUTED δ TRAPS PAVEL EXNER 1 AND ANDRII KHRABUSTOVSKYI 2 Abstract. We analyze a family of singular Schrödinger operators describing a
More informationMath 302 Outcome Statements Winter 2013
Math 302 Outcome Statements Winter 2013 1 Rectangular Space Coordinates; Vectors in the Three-Dimensional Space (a) Cartesian coordinates of a point (b) sphere (c) symmetry about a point, a line, and a
More informationClasses of Linear Operators Vol. I
Classes of Linear Operators Vol. I Israel Gohberg Seymour Goldberg Marinus A. Kaashoek Birkhäuser Verlag Basel Boston Berlin TABLE OF CONTENTS VOLUME I Preface Table of Contents of Volume I Table of Contents
More informationSHARP BOUNDARY TRACE INEQUALITIES. 1. Introduction
SHARP BOUNDARY TRACE INEQUALITIES GILES AUCHMUTY Abstract. This paper describes sharp inequalities for the trace of Sobolev functions on the boundary of a bounded region R N. The inequalities bound (semi-)norms
More informationPerturbation theory for eigenvalues of Hermitian pencils. Christian Mehl Institut für Mathematik TU Berlin, Germany. 9th Elgersburg Workshop
Perturbation theory for eigenvalues of Hermitian pencils Christian Mehl Institut für Mathematik TU Berlin, Germany 9th Elgersburg Workshop Elgersburg, March 3, 2014 joint work with Shreemayee Bora, Michael
More informationNotation. 0,1,2,, 1 with addition and multiplication modulo
Notation Q,, The set of all natural numbers 1,2,3, The set of all integers The set of all rational numbers The set of all real numbers The group of permutations of distinct symbols 0,1,2,,1 with addition
More informationIrina Pchelintseva A FIRST-ORDER SPECTRAL PHASE TRANSITION IN A CLASS OF PERIODICALLY MODULATED HERMITIAN JACOBI MATRICES
Opuscula Mathematica Vol. 8 No. 008 Irina Pchelintseva A FIRST-ORDER SPECTRAL PHASE TRANSITION IN A CLASS OF PERIODICALLY MODULATED HERMITIAN JACOBI MATRICES Abstract. We consider self-adjoint unbounded
More informationNEW RESULTS ON TRANSMISSION EIGENVALUES. Fioralba Cakoni. Drossos Gintides
Inverse Problems and Imaging Volume 0, No. 0, 0, 0 Web site: http://www.aimsciences.org NEW RESULTS ON TRANSMISSION EIGENVALUES Fioralba Cakoni epartment of Mathematical Sciences University of elaware
More informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
More informationPeriodic Schrödinger operators with δ -potentials
Networ meeting April 23-28, TU Graz Periodic Schrödinger operators with δ -potentials Andrii Khrabustovsyi Institute for Analysis, Karlsruhe Institute of Technology, Germany CRC 1173 Wave phenomena: analysis
More informationReal Analysis Notes. Thomas Goller
Real Analysis Notes Thomas Goller September 4, 2011 Contents 1 Abstract Measure Spaces 2 1.1 Basic Definitions........................... 2 1.2 Measurable Functions........................ 2 1.3 Integration..............................
More informationEIGENFUNCTIONS OF DIRAC OPERATORS AT THE THRESHOLD ENERGIES
EIGENFUNCTIONS OF DIRAC OPERATORS AT THE THRESHOLD ENERGIES TOMIO UMEDA Abstract. We show that the eigenspaces of the Dirac operator H = α (D A(x)) + mβ at the threshold energies ±m are coincide with the
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 informationLECTURE 5 APPLICATIONS OF BDIE METHOD: ACOUSTIC SCATTERING BY INHOMOGENEOUS ANISOTROPIC OBSTACLES DAVID NATROSHVILI
LECTURE 5 APPLICATIONS OF BDIE METHOD: ACOUSTIC SCATTERING BY INHOMOGENEOUS ANISOTROPIC OBSTACLES DAVID NATROSHVILI Georgian Technical University Tbilisi, GEORGIA 0-0 1. Formulation of the corresponding
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 informationElements of Positive Definite Kernel and Reproducing Kernel Hilbert Space
Elements of Positive Definite Kernel and Reproducing Kernel Hilbert Space Statistical Inference with Reproducing Kernel Hilbert Space Kenji Fukumizu Institute of Statistical Mathematics, ROIS Department
More informationOperators with Closed Range, Pseudo-Inverses, and Perturbation of Frames for a Subspace
Canad. Math. Bull. Vol. 42 (1), 1999 pp. 37 45 Operators with Closed Range, Pseudo-Inverses, and Perturbation of Frames for a Subspace Ole Christensen Abstract. Recent work of Ding and Huang shows that
More information16 1 Basic Facts from Functional Analysis and Banach Lattices
16 1 Basic Facts from Functional Analysis and Banach Lattices 1.2.3 Banach Steinhaus Theorem Another fundamental theorem of functional analysis is the Banach Steinhaus theorem, or the Uniform Boundedness
More informationSpectrum (functional analysis) - Wikipedia, the free encyclopedia
1 of 6 18/03/2013 19:45 Spectrum (functional analysis) From Wikipedia, the free encyclopedia In functional analysis, the concept of the spectrum of a bounded operator is a generalisation of the concept
More informationIntroduction to Spectral Theory
P.D. Hislop I.M. Sigal Introduction to Spectral Theory With Applications to Schrodinger Operators Springer Introduction and Overview 1 1 The Spectrum of Linear Operators and Hilbert Spaces 9 1.1 TheSpectrum
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 informationMath Introduction to Numerical Analysis - Class Notes. Fernando Guevara Vasquez. Version Date: January 17, 2012.
Math 5620 - Introduction to Numerical Analysis - Class Notes Fernando Guevara Vasquez Version 1990. Date: January 17, 2012. 3 Contents 1. Disclaimer 4 Chapter 1. Iterative methods for solving linear systems
More informationWaveguide Propagation Modes and Quadratic Eigenvalue Problems
Waveguide Propagation Modes and Quadratic Eigenvalue Problems Prof. André Nicolet, Institut Fresnel, UMR CNRS 6133, Université Paul Cézanne, Marseille, France. Prof. Christophe Geuzaine, Mathematics Department,
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 informationThe value of a problem is not so much coming up with the answer as in the ideas and attempted ideas it forces on the would be solver I.N.
Math 410 Homework Problems In the following pages you will find all of the homework problems for the semester. Homework should be written out neatly and stapled and turned in at the beginning of class
More informationSpectral theory for compact operators on Banach spaces
68 Chapter 9 Spectral theory for compact operators on Banach spaces Recall that a subset S of a metric space X is precompact if its closure is compact, or equivalently every sequence contains a Cauchy
More informationChapter 7: Bounded Operators in Hilbert Spaces
Chapter 7: Bounded Operators in Hilbert Spaces I-Liang Chern Department of Applied Mathematics National Chiao Tung University and Department of Mathematics National Taiwan University Fall, 2013 1 / 84
More informationReview and problem list for Applied Math I
Review and problem list for Applied Math I (This is a first version of a serious review sheet; it may contain errors and it certainly omits a number of topic which were covered in the course. Let me know
More informationHyperbolicity of Renormalization for C r Unimodal Maps
Hyperbolicity of Renormalization for C r Unimodal Maps A guided tour Edson de Faria Department of Mathematics IME-USP December 13, 2010 (joint work with W. de Melo and A. Pinto) Ann. of Math. 164(2006),
More informationMATH 31BH Homework 1 Solutions
MATH 3BH Homework Solutions January 0, 04 Problem.5. (a) (x, y)-plane in R 3 is closed and not open. To see that this plane is not open, notice that any ball around the origin (0, 0, 0) will contain points
More informationInequalities of Babuška-Aziz and Friedrichs-Velte for differential forms
Inequalities of Babuška-Aziz and Friedrichs-Velte for differential forms Martin Costabel Abstract. For sufficiently smooth bounded plane domains, the equivalence between the inequalities of Babuška Aziz
More informationA Spectral Characterization of Closed Range Operators 1
A Spectral Characterization of Closed Range Operators 1 M.THAMBAN NAIR (IIT Madras) 1 Closed Range Operators Operator equations of the form T x = y, where T : X Y is a linear operator between normed linear
More informationModelling in photonic crystal structures
Modelling in photonic crystal structures Kersten Schmidt MATHEON Nachwuchsgruppe Multiscale Modelling and Scientific Computing with PDEs in collaboration with Dirk Klindworth (MATHEON, TU Berlin) Holger
More informationAN EFFICIENT METHOD FOR BAND STRUCTURE CALCULATIONS IN 2D PHOTONIC CRYSTALS
AN EFFICIENT METHOD FOR BAND STRUCTURE CALCULATIONS IN 2D PHOTONIC CRYSTALS DAVID C. DOBSON Abstract. An efficient method for band structure calculations in dielectric photonic crystals is presented. The
More informationThe Dirichlet-to-Neumann operator
Lecture 8 The Dirichlet-to-Neumann operator The Dirichlet-to-Neumann operator plays an important role in the theory of inverse problems. In fact, from measurements of electrical currents at the surface
More informationCHAPTER 4 ELECTROMAGNETIC WAVES IN CYLINDRICAL SYSTEMS
CHAPTER 4 ELECTROMAGNETIC WAVES IN CYLINDRICAL SYSTEMS The vector Helmholtz equations satisfied by the phasor) electric and magnetic fields are where. In low-loss media and for a high frequency, i.e.,
More informationAPPROXIMATION OF MOORE-PENROSE INVERSE OF A CLOSED OPERATOR BY A SEQUENCE OF FINITE RANK OUTER INVERSES
APPROXIMATION OF MOORE-PENROSE INVERSE OF A CLOSED OPERATOR BY A SEQUENCE OF FINITE RANK OUTER INVERSES S. H. KULKARNI AND G. RAMESH Abstract. Let T be a densely defined closed linear operator between
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 informationSpectral Theory, with an Introduction to Operator Means. William L. Green
Spectral Theory, with an Introduction to Operator Means William L. Green January 30, 2008 Contents Introduction............................... 1 Hilbert Space.............................. 4 Linear Maps
More informationOn compact operators
On compact operators Alen Aleanderian * Abstract In this basic note, we consider some basic properties of compact operators. We also consider the spectrum of compact operators on Hilbert spaces. A basic
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 informationAbsence of bound states for waveguides in 2D periodic structures
Absence of bound states for waveguides in 2D periodic structures Maria Radosz Rice University (Joint work with Vu Hoang) Mathematical and Numerical Modeling in Optics Minneapolis, December 13 2016 1 /
More informationBOUNDARY VALUE PROBLEMS IN KREĬN SPACES. Branko Ćurgus Western Washington University, USA
GLASNIK MATEMATIČKI Vol. 35(55(2000, 45 58 BOUNDARY VALUE PROBLEMS IN KREĬN SPACES Branko Ćurgus Western Washington University, USA Dedicated to the memory of Branko Najman. Abstract. Three abstract boundary
More informationChapter 1. Preliminaries. The purpose of this chapter is to provide some basic background information. Linear Space. Hilbert Space.
Chapter 1 Preliminaries The purpose of this chapter is to provide some basic background information. Linear Space Hilbert Space Basic Principles 1 2 Preliminaries Linear Space The notion of linear space
More informationarxiv: v1 [math.gr] 8 Nov 2008
SUBSPACES OF 7 7 SKEW-SYMMETRIC MATRICES RELATED TO THE GROUP G 2 arxiv:0811.1298v1 [math.gr] 8 Nov 2008 ROD GOW Abstract. Let K be a field of characteristic different from 2 and let C be an octonion algebra
More informationINF-SUP CONDITION FOR OPERATOR EQUATIONS
INF-SUP CONDITION FOR OPERATOR EQUATIONS LONG CHEN We study the well-posedness of the operator equation (1) T u = f. where T is a linear and bounded operator between two linear vector spaces. We give equivalent
More informationReview of Some Concepts from Linear Algebra: Part 2
Review of Some Concepts from Linear Algebra: Part 2 Department of Mathematics Boise State University January 16, 2019 Math 566 Linear Algebra Review: Part 2 January 16, 2019 1 / 22 Vector spaces A set
More informationSectorial Forms and m-sectorial Operators
Technische Universität Berlin SEMINARARBEIT ZUM FACH FUNKTIONALANALYSIS Sectorial Forms and m-sectorial Operators Misagheh Khanalizadeh, Berlin, den 21.10.2013 Contents 1 Bounded Coercive Sesquilinear
More informationSpectral analysis for a class of linear pencils arising in transport theory
Spectral analysis for a class of linear pencils arising in transport theory Petru A. Cojuhari AGH University of Science and Technology Kraków, Poland Vienna, December 17 20, 2016 Petru A. Cojuhari Spectral
More informationElementary linear algebra
Chapter 1 Elementary linear algebra 1.1 Vector spaces Vector spaces owe their importance to the fact that so many models arising in the solutions of specific problems turn out to be vector spaces. The
More informationIterative Solution of a Matrix Riccati Equation Arising in Stochastic Control
Iterative Solution of a Matrix Riccati Equation Arising in Stochastic Control Chun-Hua Guo Dedicated to Peter Lancaster on the occasion of his 70th birthday We consider iterative methods for finding the
More informationOctober 25, 2013 INNER PRODUCT SPACES
October 25, 2013 INNER PRODUCT SPACES RODICA D. COSTIN Contents 1. Inner product 2 1.1. Inner product 2 1.2. Inner product spaces 4 2. Orthogonal bases 5 2.1. Existence of an orthogonal basis 7 2.2. Orthogonal
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 informationOn the Ritz values of normal matrices
On the Ritz values of normal matrices Zvonimir Bujanović Faculty of Science Department of Mathematics University of Zagreb June 13, 2011 ApplMath11 7th Conference on Applied Mathematics and Scientific
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 informationDecomposition of Riesz frames and wavelets into a finite union of linearly independent sets
Decomposition of Riesz frames and wavelets into a finite union of linearly independent sets Ole Christensen, Alexander M. Lindner Abstract We characterize Riesz frames and prove that every Riesz frame
More informationElectromagnetic fields and waves
Electromagnetic fields and waves Maxwell s rainbow Outline Maxwell s equations Plane waves Pulses and group velocity Polarization of light Transmission and reflection at an interface Macroscopic Maxwell
More informationPartial Differential Equations and the Finite Element Method
Partial Differential Equations and the Finite Element Method Pavel Solin The University of Texas at El Paso Academy of Sciences ofthe Czech Republic iwiley- INTERSCIENCE A JOHN WILEY & SONS, INC, PUBLICATION
More informationAn introduction to some aspects of functional analysis
An introduction to some aspects of functional analysis Stephen Semmes Rice University Abstract These informal notes deal with some very basic objects in functional analysis, including norms and seminorms
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 informationAN EXAMPLE OF SPECTRAL PHASE TRANSITION PHENOMENON IN A CLASS OF JACOBI MATRICES WITH PERIODICALLY MODULATED WEIGHTS
AN EXAMPLE OF SPECTRAL PHASE TRANSITION PHENOMENON IN A CLASS OF JACOBI MATRICES WITH PERIODICALLY MODULATED WEIGHTS SERGEY SIMONOV Abstract We consider self-adjoint unbounded Jacobi matrices with diagonal
More informationComputations of Critical Groups at a Degenerate Critical Point for Strongly Indefinite Functionals
Journal of Mathematical Analysis and Applications 256, 462 477 (2001) doi:10.1006/jmaa.2000.7292, available online at http://www.idealibrary.com on Computations of Critical Groups at a Degenerate Critical
More informationarxiv: v2 [math.fa] 17 May 2016
ESTIMATES ON SINGULAR VALUES OF FUNCTIONS OF PERTURBED OPERATORS arxiv:1605.03931v2 [math.fa] 17 May 2016 QINBO LIU DEPARTMENT OF MATHEMATICS, MICHIGAN STATE UNIVERSITY EAST LANSING, MI 48824, USA Abstract.
More informationWave equation on manifolds and finite speed of propagation
Wave equation on manifolds and finite speed of propagation Ethan Y. Jaffe Let M be a Riemannian manifold (without boundary), and let be the (negative of) the Laplace-Beltrami operator. In this note, we
More informationEXISTENCE OF GUIDED MODES ON PERIODIC SLABS
SUBMITTED FOR: PROCEEDINGS OF THE FIFTH INTERNATIONAL CONFERENCE ON DYNAMICAL SYSTEMS AND DIFFERENTIAL EQUATIONS June 16 19, 2004, Pomona, CA, USA pp. 1 8 EXISTENCE OF GUIDED MODES ON PERIODIC SLABS Stephen
More informationKernel Method: Data Analysis with Positive Definite Kernels
Kernel Method: Data Analysis with Positive Definite Kernels 2. Positive Definite Kernel and Reproducing Kernel Hilbert Space Kenji Fukumizu The Institute of Statistical Mathematics. Graduate University
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 information6 Inner Product Spaces
Lectures 16,17,18 6 Inner Product Spaces 6.1 Basic Definition Parallelogram law, the ability to measure angle between two vectors and in particular, the concept of perpendicularity make the euclidean space
More informationhere, this space is in fact infinite-dimensional, so t σ ess. Exercise Let T B(H) be a self-adjoint operator on an infinitedimensional
15. Perturbations by compact operators In this chapter, we study the stability (or lack thereof) of various spectral properties under small perturbations. Here s the type of situation we have in mind:
More informationKeywords: Bloch-Floquet Theorem, Periodic Media, Photonic Crystals, Plasmonic Crystals, Electronic Crystals, Phononic Crystals
On the Bloch Theorem and Orthogonality Relations Sina Khorasani School of Electrical Engineering Sharif University of Technology P. O. Box 11365-9363 Tehran, Iran Email: khorasani@sina.sharif.edu Abstract
More informationFoundations of Matrix Analysis
1 Foundations of Matrix Analysis In this chapter we recall the basic elements of linear algebra which will be employed in the remainder of the text For most of the proofs as well as for the details, the
More informationFinite Element Analysis of Acoustic Scattering
Frank Ihlenburg Finite Element Analysis of Acoustic Scattering With 88 Illustrations Springer Contents Preface vii 1 The Governing Equations of Time-Harmonic Wave Propagation, 1 1.1 Acoustic Waves 1 1.1.1
More informationMath 307 Learning Goals. March 23, 2010
Math 307 Learning Goals March 23, 2010 Course Description The course presents core concepts of linear algebra by focusing on applications in Science and Engineering. Examples of applications from recent
More informationMATH 240 Spring, Chapter 1: Linear Equations and Matrices
MATH 240 Spring, 2006 Chapter Summaries for Kolman / Hill, Elementary Linear Algebra, 8th Ed. Sections 1.1 1.6, 2.1 2.2, 3.2 3.8, 4.3 4.5, 5.1 5.3, 5.5, 6.1 6.5, 7.1 7.2, 7.4 DEFINITIONS Chapter 1: Linear
More informationCOMPACT OPERATORS. 1. Definitions
COMPACT OPERATORS. Definitions S:defi An operator M : X Y, X, Y Banach, is compact if M(B X (0, )) is relatively compact, i.e. it has compact closure. We denote { E:kk (.) K(X, Y ) = M L(X, Y ), M compact
More informationRemark By definition, an eigenvector must be a nonzero vector, but eigenvalue could be zero.
Sec 6 Eigenvalues and Eigenvectors Definition An eigenvector of an n n matrix A is a nonzero vector x such that A x λ x for some scalar λ A scalar λ is called an eigenvalue of A if there is a nontrivial
More informationOn Riesz-Fischer sequences and lower frame bounds
On Riesz-Fischer sequences and lower frame bounds P. Casazza, O. Christensen, S. Li, A. Lindner Abstract We investigate the consequences of the lower frame condition and the lower Riesz basis condition
More information285K Homework #1. Sangchul Lee. April 28, 2017
285K Homework #1 Sangchul Lee April 28, 2017 Problem 1. Suppose that X is a Banach space with respect to two norms: 1 and 2. Prove that if there is c (0, such that x 1 c x 2 for each x X, then there is
More information