Galerkin approximation of holomorphic eigenvalue problems: weak T-coercivity and T-compatibility
Martin Halla

TL;DR
This paper develops a new framework using weak T-coercivity and T-compatibility concepts to analyze Galerkin approximations of holomorphic eigenvalue problems, ensuring their regularity and convergence.
Contribution
It introduces the concepts of weak T-coercivity and T-compatibility, providing a general technique for proving regularity of Galerkin approximations in complex eigenvalue problems.
Findings
Framework improves previous results on eigenvalue approximation.
Applicable to a wide range of holomorphic Fredholm operator problems.
Ensures convergence and regularity of Galerkin methods in non-coercive settings.
Abstract
We consider Galerkin approximations of holomorphic Fredholm operator eigenvalue problems for which the operator values don't have the structure "coercive+compact". In this case the regularity (in sense of [O. Karma, Numer. Funct. Anal. Optim. 17 (1996)]) of Galerkin approximations is not unconditionally satisfied and the question of convergence is delicate. We report a technique to prove regularity of approximations which is applicable to a wide range of eigenvalue problems. In particular, we introduce the concepts of weak T-coercivity and T-compatibility and prove that for weakly T-coercive operators, T-compatibility of Galerkin approximations implies their regularity. Our framework immediately improves the results of [T. Hohage, L. Nannen, BIT 55(1) (2015)], is immediately applicable to analyze approximations of eigenvalue problems related to [A.-S. Bonnet-Ben Dhia, C. Carvalho, P.…
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Galerkin approximation of holomorphic eigenvalue problems:
weak T-coercivity and T-compatibility
Martin Halla
(Date: August 12th, 2019.)
Abstract.
We consider Galerkin approximations of holomorphic Fredholm operator eigenvalue problems for which the operator values don’t have the structure “coercive+compact”. In this case the regularity (in sense of [O. Karma, Numer. Funct. Anal. Optim. 17 (1996)]) of Galerkin approximations is not unconditionally satisfied and the question of convergence is delicate. We report a technique to prove regularity of approximations which is applicable to a wide range of eigenvalue problems. In particular, we introduce the concepts of weak T-coercivity and T-compatibility and prove that for weakly T-coercive operators, T-compatibility of Galerkin approximations implies their regularity.
Our framework immediately improves the results of [T. Hohage, L. Nannen, BIT 55(1) (2015)], is immediately applicable to analyze approximations of eigenvalue problems related to [A.-S. Bonnet-Ben Dhia, C. Carvalho, P. Ciarlet, Num. Math. 138(4) (2018)] and is already applied in [G. Unger, preprint (2017)].
Key words and phrases:
holomorphic eigenvalue problem, non-linear eigenvalue problem, approximation analysis, T-coercivity.
2010 Mathematics Subject Classification:
47J10, 65H17, 65N25.
The analysis of approximations for holomorphic Fredholm operator eigenvalue problems has a long history [15], [24], [25], [17], [18], [19] and is usually performed in the framework of discrete approximation schemes [21] and regular approximations of operator functions [14], [1]. In this framework a complete convergence analysis and asymptotic error estimates for eigenvalues are given by Karma in [18], [19]. If the discrete approximation scheme is chosen as a Galerkin scheme, then the assumptions of [18], [19] reduce to a single non-trivial assumption: the regular approximation property (see Definition 1.5 for the meaning of regularity). If the operator values are of the form “coercive+compact”, the regularity of Galerkin approximations is unconditionally satisfied. However, if the operator values are not of this kind the question of spectrally converging approximations is very delicate. This can already be observed for linear eigenvalue problems, see e.g. [3], [2]. Though it is little known how to prove regularity of approximations. In Theorem 1.8 we report a new condition on the Galerkin spaces to ensure the regularity of Galerkin approximations such that [18], [19] can be applied. This condition is stronger than the classical regularity condition. However, it suffices for a wide variety of applications. On the side, we report in Lemma 2.6 new asymptotic error estimates on eigenspaces for regular Galerkin approximations (which are not provided by [18], [19]). The latter is an improvement of Unger [22, Theorem 4.3.7]. We combine our approach with the results of [18], [19] in Proposition 2.7 and Corollary 2.8.
As preparation for the forthcoming concept of weakly T-coercive operators (operator functions) we remind the reader how Fredholmness of operators is usually established. In the case of coercive operators Fredholmness is trivial. The same holds for weakly coercive operators , i.e. is a compact perturbation of a coercive operator. Else wise we may construct an isomorphism such that is weakly coercive ( denotes the adjoint operator of ), which yields the Fredholmness of . The name “T-coercivity” originates from Bonnet-Ben Dhia, Ciarlet, Zwölf [5]. The notion was introduced to analyze differential operators with sign-changing coefficients in the principal part which occur e.g. in the modeling of meta materials. The technique is also applied in the analysis of interior transmission eigenvalue problems, see e.g. [9], [8]. Although as far as we know, the concept goes back to a remark by Buffa [6] (wherein ) on non-coercive operators with applications to Maxwell equations. For an operator to be (weakly) -coercive means that is already (weakly) coercive. However, in eigenvalue problems the operator values will be in general not bijective (precisely at the eigenvalues). Thus the nomenclature of T-coercivity is not meaningful for our purposes and we will rely on the term weak T-coercivity. In general the Galerkin spaces will not be -invariant and hence one cannot reproduce the above analysis on the approximation level. An invariance condition is indeed not necessary, but can be relaxed. We will make precise in which sense the Galerkin spaces have to interact with the operator to ensure regularity. It will turn out that the existence of bounded linear operators from the Galerkin spaces to themselves such that
[TABLE]
with
[TABLE]
is sufficient. We call this property “-compatibility”. The norm (2) was termed “discrete norm” by Descloux, Nassif and Rappaz [10], [11] wherein it was used in a different but familiar context. In our context it was already employed by Hohage and Nannen [16] for the analysis of perfectly matched layer and Hardy space infinite element methods in cylindrical waveguides; and also by Bonnet-Ben Dhia, Ciarlet and Carvalho [7], [4] for the analysis of finite element methods for equations which involve meta materials. Both works [16], [4] prove weak T-coercivity and T-compatibility. Thus our results can directly be applied to improve the results of [16] and to establish convergence results for approximations of the eigenvalue problems related to [4]. Note that the negative material parameters in meta materials are e.g. of the kind with being the eigenvalue parameter. Hence such eigenvalue problems are indeed non-linear.
However, the original motivation for this article was to provide a framework for the convergence analysis of boundary element discretizations of boundary integral formulations of Maxwell eigenvalue problems and is already applied by Unger [23]. Although the Maxwell eigenvalue problem is of linear nature, its formulation as boundary integral equation becomes non-linear due to the dependency of the fundamental solution on the frequency.
The remainder of this article is structured as follows. In Section 1 we introduce the notion of weak T-coercivity and T-compatibility. In Theorem 1.8 we prove that T-compatibility implies regularity. In Section 2 we report in Lemma 2.6 an approximation result on eigenspaces for regular Galerkin approximations of holomorphic Fredholm operator eigenvalue problems. We merge our results with the results of Karma [18], [19] in Proposition 2.7 and Corollary 2.8.
1. Weak T-coercivity and T-compatibility
Let be a Hilbert space with scalar body and scalar product and associated norm . Let be the space of bounded linear operators from to with operator norm for . For we denote its adjoint operator by , i.e. for all . For a closed subspace let be the space of bounded linear operators from to with norm for and denote the orthogonal projection from to . Henceforth we assume that is a sequence of closed subspaces of such that converges point-wise to the identity, i.e. for each .
Definition 1.1**.**
Let and be bijective. The operator is called
- (1)
coercive, if , 2. (2)
weakly coercive, if there exists a compact operator such that is coercive, 3. (3)
-coercive if is coercive, 4. (4)
weakly -coercive if is weakly coercive.
Due to the Lemma of Lax-Milgram every coercive operator is invertible. Every weakly -coercive operator is Fredholm with index zero. For a (weakly) coercive operator it is true that the Galerkin approximations inherit the (weak) coercivity, while for (weakly) -coercive operators it is in general wrong.
We note that if is weakly coercive, then is so too. Vice-versa, if is weakly coercive, then so is . Hence we could alternatively define to be (weakly) right -coercive, if is (weakly) coercive. However, we stick to the former variant because it is more convenient.
For an operator or or a sum of such we define the “discrete norm”
[TABLE]
Definition 1.2**.**
Consider and . We say that converges to in discrete norm, if
[TABLE]
We define in the following what we mean by -compatible approximations of weakly -coercive operators.
Definition 1.3**.**
Let be weakly -coercive. Then we call the sequence of Galerkin approximations -compatible, if is a sequence of index zero Fredholm operators and there exists a sequence of index zero Fredholm operators such that converges to in discrete norm: .
Definition 1.4**.**
A sequence is said to be compact, if for every subsequence exists in turn a converging subsubsequence.
Definition 1.5**.**
A sequence is called regular, if for every bounded sequence the compactness of already implies the compactness of .
Next we briefly elaborate on the notion of regularity for readers who are totally unfamiliar with this concept. Regularity of Galerkin approximations is a meaningful generalization of stability and well suited for the approximation analysis of eigenvalue problems. Consider for example bijective and its Galerkin approximation . In this case regularity of implies stability: Assume that is not stable. Thus there exists with for each such that . If is regular, there exists a subsequence and such that . It follows . Since is bijective, it follows which is a contradiction to .
On the other hand, consider a holomorphic Fredholm operator function with non-empty resolvent set and sequences of eigenvalues with normalized eigenelements of the Galerkin approximation (i.e. ) such that (see Section 2 for definitions and details). If is regular for each , then is indeed an eigenvalue of (i.e. there occurs no spectral pollution): Due to the continuity of with respect to , implies . If is regular, there exists a subsequence and such that . It follows and , i.e. is an eigenvalue of with normalized eigenelement .
Our next goal is to prove in Theorem 1.8 that -compatible Galerkin approximations of weakly -coercive operators are regular. In preparation we formulate the next two lemmata.
Lemma 1.6**.**
Let and be a sequence of operators with and . Then there exist a constant and an index such that
[TABLE]
for all . If is bijective and is Fredholm with index zero for each , then there exist a constant and an index such that is also bijective for all and
[TABLE]
Proof.
Let . With the triangle inequality we deduce
[TABLE]
and hence
[TABLE]
Since the right hand side of the previous inequality is bounded. Similar, with the inverse triangle inequality we deduce
[TABLE]
and hence
[TABLE]
It hold and . Thus let be such that and for all . It follows
[TABLE]
for all . For the last claim let be such that for all . Again with the inverse triangle inequality and
[TABLE]
it follows
[TABLE]
for all . We deduce that is injective. Since is Fredholm with index zero its bijectivity follows. The norm estimate holds due to . ∎
Lemma 1.7**.**
Let be weakly -coercive and be compact such that is coercive. Let be a -compatible Galerkin approximation of . Then there exist and , such that is invertible and
[TABLE]
for all .
Proof.
Let be large enough such that is bijective (see Lemma 1.6). We compute
[TABLE]
with coercivity constant
[TABLE]
Since is uniformly bounded from above and below (see Lemma 1.6) and converges to in discrete norm by assumption, it follows the existence of and such that
[TABLE]
for all . Hence is injective. Since is Fredholm with index zero and is compact, is Fredholm with index zero too. Thus is bijective. The norm estimate follows now from
[TABLE]
for any bijective . ∎
Theorem 1.8**.**
Let be weakly -coercive and
[TABLE]
be a -compatible Galerkin approximation. Then is regular.
Proof.
Without loss of generality let be a bounded sequence, and be such that . Let be compact such that is coercive. Let and . Since is compact and is bounded, there exist a subsequence and such that . It follows
[TABLE]
Due to Lemma 1.7 there exist and , such that for all operator is invertible and . For we compute
[TABLE]
The first term on the right hand side of the latter inequality converges to zero, as previously discussed. The second and third term converge to zero, because converges point-wise to the identity. Hence
[TABLE]
∎
2. Holomorphic eigenvalue problems
We refer the reader to [13] and [20, Appendix] for theory on holomorphic (Fredholm) operator functions. Let be an open, connected and non-empty subset of . Let be an operator function. An operator function is called holomorphic, if it is complex differentiable. An operator function is called Fredholm, if is Fredholm for each . We denote the resolvent set and spectrum of an operator function as
[TABLE]
For an operator function we denote by the operator function defined by for each and by A^{-1}(\cdot)\colon\rho\big{(}A(\cdot)\big{)}\to L(X) the operator function defined by for each \lambda\in\rho\big{(}A(\cdot)\big{)}. Note that for a holomorphic operator function the operator function defined by is holomorphic as well. For a holomorphic operator function denote by the derivative of . It is well known (see e.g. [12, Theorem 8.2]) that for a holomorphic Fredholm operator function such that is bijective for at least one , the spectrum \sigma\big{(}A(\cdot)\big{)} is discrete, has no accumulation points in and every \lambda\in\sigma\big{(}A(\cdot)\big{)} is an eigenvalue. That is, there exists such that . In this case we call an eigenelement. An ordered collection of elements in is called a Jordan chain at if is an eigenelement corresponding to and if
[TABLE]
The elements of a Jordan chain are called generalized eigenelements and the closed linear hull of all generalized eigenelements of at is called the generalized eigenspace for at . For an eigenelement we denote by the maximal length of a Jordan chain at beginning with and
[TABLE]
The maximal length of a Jordan chain is always finite, see e.g. [20, Lemma A.8.3]. Next we generalize Definitions 1.1, 1.3, 1.5 and Theorem 1.8 to operator functions.
Definition 2.1**.**
Let be operator functions and \rho\big{(}T(\cdot)\big{)}=\Lambda. is (weakly) (-)coercive, if is (weakly) (-)coercive for each .
Definition 2.2**.**
Let be weakly -coercive. Then we call the sequence of Galerkin approximations -compatible, if is compatible for each .
Definition 2.3**.**
Let be an operator function. The sequence of Galerkin approximations is regular, if is regular for each
Theorem 2.4**.**
Let be weakly -coercive and
[TABLE]
be a -compatible Galerkin approximation. Then is regular.
Proof.
Follows from Theorem 2.4. ∎
Next we prepare to apply [18], [19].
Lemma 2.5**.**
Let be a holomorphic Fredholm operator function and let be a sequence of closed subspaces of with orthogonal projections onto , such that converges point-wise to the identity. Then the Galerkin scheme \big{(}P_{n}A(\cdot)|_{X_{n}})_{n\in\mathbb{N}} is a discrete approximation scheme in the sense of [18].
Proof.
For a Galerkin scheme it holds with the notation of [18]
[TABLE]
Assumptions a1)-a4) of [18] follow all from the point-wise convergence of . ∎
Next we generalize Theorem 4.3.7 of [22].
Lemma 2.6**.**
Let be open, be a Hilbert space and be the space of bounded linear operators from to . Let be a holomorphic Fredholm operator function with non-empty resolvent set and be a sequence of closed subspaces of with orthogonal projections onto , such that converges point-wise to the identity, i.e. for all . Let be the Galerkin approximation of defined by for each . Let the assumptions of [18, Theorem 2, Theorem 3] and [19, Theorem 2, Theorem 3] be satisfied. Let be a compact set with rectifiable boundary \partial\tilde{\Lambda}\subset\rho\big{(}A(\cdot)\big{)} and \tilde{\Lambda}\cap\sigma\big{(}A(\cdot)\big{)}=\{\lambda_{0}\}. Then there exist and such that for all
[TABLE]
for all \lambda_{n}\in\sigma\big{(}A_{n}(\cdot)\big{)}\cap\tilde{\Lambda} and all with .
Proof.
We proceed as in [22]: Theorem 4.3.7 of [22] requires a special form of the operator function . However its proof uses this assumption only to apply Lemma 4.2.1 of [22]. Hence we need to establish the result of [22, Lemma 4.2.1] without the assumption on the special form of . However, the result of [22, Lemma 4.2.1] already follows from [18, Theorem 2 ii)].
∎
Next we apply [18], [19] and Lemma 2.6.
Proposition 2.7**.**
Let be open, connected and non-empty, be a Hilbert space and be the space of bounded linear operators from to itself. Let be a holomorphic Fredholm operator function with non-empty resolvent set \rho\big{(}A(\cdot)\big{)}\neq\emptyset. Let be a sequence of closed subspaces of with orthogonal projections onto , such that converges point-wise to the identity, i.e. for each . Let be the Galerkin approximation of defined by for each . Assume that is Fredholm with index zero for each and . Assume that is a regular approximation of (see Definition 2.3). Then the following results hold.
- i)
For every eigenvalue of exists a sequence converging to with being an eigenvalue of for almost all . 2. ii)
Let be a sequence of normalized eigenpairs of , i.e.
[TABLE]
and , so that , then
- a)
* is an eigenvalue of ,* 2. b)
* is a compact sequence and its cluster points are normalized eigenelements of .* 3. iii)
For every compact the sequence is stable on , i.e. there exist and such that for all and all . 4. iv)
For every compact with rectifiable boundary \partial\tilde{\Lambda}\subset\rho\big{(}A(\cdot)\big{)} exists an index such that
[TABLE]
for all , whereby denotes the generalized eigenspace of an operator function at .
Let be a compact set with rectifiable boundary \partial\tilde{\Lambda}\subset\rho\big{(}A(\cdot)\big{)}, \tilde{\Lambda}\cap\sigma\big{(}A(\cdot)\big{)}=\{\lambda_{0}\} and
[TABLE]
whereby denotes the complex conjugate of and the adjoint operator function of defined by for each . Then there exist and such that for all
- v)
[TABLE]
for all \lambda_{n}\in\sigma\big{(}A_{n}(\cdot)\big{)}\cap\tilde{\Lambda}, whereby denotes the maximal length of a Jordan chain of at the eigenvalue , 2. vi)
[TABLE]
whereby is the weighted mean of all the eigenvalues of in
[TABLE] 3. vii)
[TABLE]
for all \lambda_{n}\in\sigma\big{(}A_{n}(\cdot)\big{)}\cap\tilde{\Lambda} and all with .
Proof.
The first three claims follow with [18, Theorem 2], if we can proof that the required assumptions are satisfied. First of all a Galerkin scheme is a discrete approximation scheme due to Lemma 2.5. The operator function are holomorphic by assumption. It follows that is also holomorphic. and are index zero Fredholm operator functions by assumption. Assumption b1 \rho\big{(}A(\cdot)\big{)}\neq\emptyset is also an assumption of this theorem. Assumption b2 follows from Lemma 1.7 (at least for sufficiently large ). Assumption b3 follows from . Assumption b4 follows from the point-wise convergence of the projections . Assumption b5 is also an assumption of this theorem.
The fourth claim follows with [18, Theorem 3], if we can proof the required assumption (R). We can chose as injection, i.e. . Hence . Since ii) follows from the point-wise convergence of the projections .
The fifth and sixth claim follow with [19, Theorem 2, Theorem 3], if we can proof their required assumptions. Assumption a1-a4 are canonical satisfied by Galerkin schemes. We already proved that Assumptions b1-b5 are satisfied. We can chose . For [19, Theorem 3] we can chose the same as before.
For the proof of the seventh claim we refer to Lemma 2.6. ∎
Finally we combine Theorem 2.4 and Proposition 2.7.
Corollary 2.8**.**
Let be open, connected and non-empty, be a Hilbert space and be the space of bounded linear operators from to . Let be a holomorphic weakly -coercive operator function (see Definition 2.1) with non-empty resolvent set \rho\big{(}A(\cdot)\big{)}\neq\emptyset. Let be a sequence of closed subspaces of with orthogonal projections onto , such that converges point-wise to the identity, i.e. for each . Let be the Galerkin approximation of defined by for each . Assume that is -compatible (see Definition 2.2). Then results i)-vii) of Proposition 2.7 hold.
Proof.
Since is weakly -coercive, it is Fredholm with index zero. Since is -compatible, it is Fredholm with index zero and regular. ∎
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] P. M. Anselone and M. L. Treuden, Regular operator approximation theory , Pacific J. Math. 120 (1985), no. 2, 257–268. MR 810769 (87a:47017)
- 2[2] Douglas N. Arnold, Richard S. Falk, and Ragnar Winther, Finite element exterior calculus: from Hodge theory to numerical stability , Bull. Amer. Math. Soc. (N.S.) 47 (2010), no. 2, 281–354. MR 2594630 (2011 f:58005)
- 3[3] Daniele Boffi, Franco Brezzi, and Lucia Gastaldi, On the problem of spurious eigenvalues in the approximation of linear elliptic problems in mixed form , Math. Comp. 69 (2000), no. 229, 121–140. MR 1642801 (2000 i:65175)
- 4[4] Anne-Sophie Bonnet-Ben Dhia, Camille Carvalho, and Patrick Ciarlet, Mesh requirements for the finite element approximation of problems with sign-changing coefficients , Numerische Mathematik 138 (2018), no. 4, 801–838.
- 5[5] Anne-Sophie Bonnet-Ben Dhia, Patrick Ciarlet, and Carlo Maria Zwölf, Time harmonic wave diffraction problems in materials with sign-shifting coefficients , Journal of Computational and Applied Mathematics 234(6) (2010), 1912–1919.
- 6[6] Annalisa Buffa, Remarks on the discretization of some noncoercive operator with applications to heterogeneous maxwell equations , SIAM Journal on Numerical Analysis 43 (2005), no. 1, 1–18.
- 7[7] Camille Carvalho, Etude mathématique et numérique de structures plasmonique avec des coins. , Ph.D. thesis, 12 2015.
- 8[8] Lucas Chesnel, Interior transmission eigenvalue problem for maxwell’s equations: the T-coercivity as an alternative approach , Inverse problems 28 (2012), 065005.
