arxiv: v2 [math.fa] 24 Apr 2017

Transcription

1 ENCLOSURE OF THE NUMERICAL RANGE OF A CLASS OF NON-SELFADJOINT RATIONAL OPERATOR FUNCTIONS arxiv: v2 [math.fa] 24 Apr 217 CHRISTIAN ENGSTRÖM AND AXEL TORSHAGE Abstract. In this paper we introduce an enclosure of the numerical range of a class of rational operator functions. In contrast to the numerical range the presented enclosure can be computed exactly in the infinite dimensional case as well as in the finite dimensional case. Moreover, the new enclosure is minimal given only the numerical ranges of the operator coefficients and many characteristics of the numerical range can be obtained by investigating the enclosure. We introduce a pseudonumerical range and study an enclosure of this set. This enclosure provides a computable upper bound of the norm of the resolvent. 1. Introduction The spectral properties of operator functions play an important role in mathematical analysis and in many applications [Tre8, Lif89, APT2]. A classic enclosure of the spectrum is the closure of the numerical range [Mar88]. Furthermore, the norm of the resolvent in a point ω is under some conditions bounded by a quantity that depend on the distance from ω to the numerical range [MM1]. Knowledge of the numerical range is also important in perturbation theory and in several other branches of operator theory [Kat95]. However, in most cases it is not possible to analytically determine the numerical range, not even in the finite dimensional case. The geometric properties of the numerical range of matrix polynomials and rational matrix functions have been studied extensively [LR94, AMP2] and it is possible to numerically approximate the shape of the numerical range of matrix polynomials [CNP2]. However, as matrix functions generated by a discretization of a differential equation are very large, the available algorithms are very time consuming. Furthermore, the methods developed for matrix problems are not applicable in the infinite dimensional case. In this paper we introduce an enclosure of the numerical range of a class of rational operator functions whose values are linear operators in a Hilbert space H. Importantly, this new enclosure is applicable in the infinite dimensional case as well as in the finite dimensional case. Let A and B be selfadjoint operators in H, where B is non-zero and bounded. We consider rational operator functions of the form (1.1) T pωq : A ω 2 ω 2 c idω ω 2 B, dom T pωq dom A, ω P Cztδ`, δ u, where c and d are non-negative real numbers, and δ are the poles of the coefficient of B. If d, then the operator function ω 2 ÞÑ T pω 2 q is selfadjoint. This function 1991 Mathematics Subject Classification. 47J1, 47A56, 47A12. Key words and phrases. Non-linear spectral problem, numerical range, psedospectra, resolvent estimate. 1

2 2 CHRISTIAN ENGSTRÖM AND AXEL TORSHAGE has been studied extensively [FM91, AL95, ALL1, EL4]. In the case B ě the rational function is the first Schur complement of a selfadjoint block operator matrix [AL95, Tre8]. The non-selfadjoint case, d ą (as well as the case d ), has applications in electromagnetic field theory and cover important applications in optics [Zio3, EKE12, ELT17]. The presented enclosure of the numerical range is minimal, given only the numerical ranges of A and of B, and we will show that this enclosure can be computed exactly. Resolvent estimates and pseudospectra can be used to investigate quantitative properties of non-normal operators and operator functions [Dav7, TE5]. In particular, estimates of the resolvent of bounded analytic operator functions were considered in [MM1]. To derive a computable estimate for (1.1), we introduce a pseudonumerical range and study an enclosure of this set. The derived enclosure of the pseudonumerical range provides a computable upper bound of the norm of the resolvent in the complement of the new enclosure of the numerical range. This enclosure of the pseudospectra can be used to understand how the resolvent behaves outside the enclosure of the numerical range. Moreover, the enclosure of the pseudospectra shows where the resolvent potentially is large and can in the finite dimensional case be combined with a numerical estimate of the pseudospectra [TE5]. The the paper is organized as follows: In Section 2, we present the enclosure of the numerical range, the theoretical framework used in the paper, and conditions for determining if ω P C belong to the enclosure. Our main results are Theorem 2.9 and the algorithm in Proposition 2.19, which can be used to determine the enclosure of the numerical range. In Section 3, properties of the boundary of the enclosure are analyzed in detail. Our main results are conditions for the existence of a strip in the complement of the numerical range given in Proposition 3.15 and in Proposition Moreover, Proposition 3.16 and Proposition 3.31 provide important properties of the strip. In Section 4, the ɛ-pseudonumerical range is introduced and we determine an enclosure of this set. Our main results are Theorem 4.3, which shows how the boundary of the enclosure of the pseudospectra can be determined and Corollary 4.6 gives an estimate of the resolvent of (1.1). Throughout this paper, we use the following notation. Let ω R and ω I denote the real and imaginary parts of ω, respectively. If M is a subset of an Euclidean space, then BM denotes the boundary of M. Further, we denote by? the principal square root. 2. Enclosure of the numerical range In this section we derive an enclosure of the numerical range of the operator function (1.1). Define for a non-negative real number c and a positive d the constants c (2.1) θ : c d2 4, δ : θ i d 2. Note that the operator-valued function (1.1) is defined for ω P C : Cztδ`, δ u, where δ` and δ are the poles of T and the domain is independent of ω P C. For u P dom T ztu we define the functionals α u pau, uq{pu, uq, β u pbu, uq{pu, uq,

3 NUMERICAL RANGE OF NONSELFADJOINT RATIONAL OPERATOR FUNCTIONS 3 and (2.2) t pαu,β uqpωq : pt pωqu, uq pu, uq α u ω 2 ω 2 c idω ω 2 β u, ω P C. The numerical range of T is by definition ď W pt q : tω P C : t pαu,β uqpωq u. updom Aztu For convenience we will in some cases not explicitly write the dependence of u in the functionals α u and β u. To simplify the investigation of W pt q, we define the polynomial (2.3) p pα,βq pωq : t pα,βq pωqpc idω ω 2 q pα ω 2 qpc idω ω 2 q βω 2. For fixed values on the constants c and d we order the roots (2.4) r n : R ˆ R Ñ C, n 1,..., 4, of p α,β such that they are continuous functions of pα, βq P R 2. The numerical range of T can then be written as 4ď ď (2.5) W pt q r n pα u, β u q. n 1 updom Aztu From (2.5) it is apparent that W pt q consists of at most four components, i.e., W pt q is a union of at most four maximal connected subsets of W pt q. Let R : R Y t 8u denote the extended line of real numbers and denote by C : C Y t8u the Riemann sphere. We extend the functions r n, n 1,..., 4 to r n : R ˆ R Ñ C such that the extension coincides with the limit values. For a given set X Ă RˆR let W X pt q Ă C denote the set 4ď (2.6) W X pt q : r n pxq, r n pxq : ď r n pα, βq. n 1 pα,βqpx The roots r n of (2.3) are given by particular pairs pα u, β u q in (2.7) Ω : W paq ˆ W pbq Ă R ˆ R. Hence, by taking X Ω in (2.6), we get the enclosure 4ď (2.8) W Ω pt q r n pωq Ą W pt q. n 1 Moreover, from this definition it follows that W Ω pt q is the minimal set that encloses W pt q given only W paq and W pbq. Lemma 2.1. The polynomial p pα,βq defined in (2.3) has in the limits α Ñ 8 the roots δ`, δ, and 8, where 8 is a double root. Proof. Define p 2 pωq : pω δ`qpω δ q, then the roots of p α,β coincide with those of p pα,βq pωq (2.9) p 2 pωq ` ω 2 β p 2pωq. α α The poles δ` and δ are roots of p 2 and (2.9) is for large α a small perturbation of p 2. Then, since the roots of a polynomial depend continuously on its coefficients,

4 4 CHRISTIAN ENGSTRÖM AND AXEL TORSHAGE δ` and δ are roots in the limits α Ñ 8. There can be no other finite roots in the limit since the perturbation of p 2 is arbitrary small. Proposition 2.2. The enclosure W Ω pt q as defined in (2.8) has the following properties: i) W Ω pt q is symmetric with respect to the imaginary axis. ii) P W Ω pt q if and only if P W paq or c. iii) δ` P W Ω pt q if and only if W paq is unbounded or P W pbq or c. iv) δ P W Ω pt q if and only if W paq is unbounded or P W pbq. v) 8 P W Ω pt q if and only if W paq is unbounded. Proof. i) The polynomial p pα,βq piωq has real coefficients. Hence, the symmetry follows from the complex conjugate root theorem. ii) Follows directly from (2.3) and (2.8). iii) c implies δ` and δ` P W Ω pt q then follows from ii). The number p pα,βq pδ`q βδ 2` is zero for β, which implies δ` P W Ω pt q if P W pbq. If W paq is unbounded the statement follows directly from Lemma 2.1. Suppose none of the above holds, then p pα,βq pδ`q βδ 2`, and since W paq is bounded, p pα,βq pωq in a neighborhood of δ`. The proof of iv) is similar to iii) with the difference δ for c. v) is immediate from Lemma 2.1. Corollary 2.3. Let W Ω pt q denote the enclosure (2.8) and take ω P t, δ`, δ, 8u. Then ω P W Ω pt q if and only if r n pα, βq ω for some n P t1, 2, 3, 4u and pα, βq P BΩ. Proof. Similar to Proposition 2.2. The following propositions provide simple tests for ω P W Ω pt q. Proposition 2.4. Let W Ω pt q denote the enclosure (2.8) and assume that ω is a point on the imaginary axis with ω iω I P irzt, δ`, δ u. Then ω P W Ω pt q if and only if at least one of following conditions hold: (2.1) ω 2 I ωi 2 c ` dω I ` ωi 2 inf W pbq P W paq, ω 2 I ωi 2 c ` dω I ` ωi 2 sup W pbq P W paq, c ` dω I ` ωi 2 ωi 2 pωi 2 ` inf W paqq P W pbq. Proof. By definition iω I P W Ω pt q X irzt, δ`, δ u if and only if there exists a pα, βq P Ω such that (2.11) α ωi 2 c ` dω I ` ωi 2 β. Thus α is a non-constant real linear function in β. Since pα, βq P Ω and β belongs to a bounded set, r n pα 1, β 1 q iω I for some pair pα 1, β 1 q P BΩ. Equation (2.11) has two solutions unless the pair is a corner of Ω. Hence it is enough to investigate three of the line segments on BΩ to determine if iω I P W Ω pt q. The converse holds trivially. Let D denote the open disk! ˇ (2.12) D : ω : ˇω ` i c ˇ ă c ) Ă C. d d ω 2 I

5 NUMERICAL RANGE OF NONSELFADJOINT RATIONAL OPERATOR FUNCTIONS 5 Lemma 2.5. Let W Ω pt q denote the enclosure (2.8) and denote by D the disk (2.12), then it holds that BD X W Ω pt q Ă t, δ`, δ, 2ic{du. Proof. Assume that ω P BDztδ`, δ uxw Ω pt q, then the imaginary part of t pα,βq pωq in (2.2) is 2ω R ω I, which is zero only for ω P t, 2ic{du. Proposition 2.6. Let W Ω pt q and D be the enclosure (2.8) and the disk (2.12), respectively. Take ω P CzpiRYtδ`, δ, 8uq. Then ω P W Ω pt q if and only if ω R BD and 2 2ω I `p ω (2.13) ˆβpωq : R ` ωi 2 ` dω I ` cq 2 ` ωr 2 p2ω I ` dq 2 d ω 2 P W pbq, ` 2cω I and (2.14) ˆαpωq : p2ω I ` dq ω 4 d ω 2 P W paq. ` 2cω I Proof. Assume that ω P W Ω pt q for some ω R ir Y tδ`, δ u, then the real and imaginary parts of the equality t pα,βq pωq give the following linear system of equations: (2.15) 2ω I `p ω 2 R ` ω 2 I ` dω I ` cq 2 ` ω 2 Rp2ω I ` dq 2 pd ω 2 ` 2cω I qβ, (2.16) p2ω I ` dq ω 4 pd ω 2 ` 2cω I qα. The expression d ω 2 `2cω I is only zero for ω P BD. Hence (2.13) and (2.14) follows from (2.15), (2.16), and Lemma 2.5. Define the sets Π (2.17) β : tω P CzD : ω I ď u, Π α : tω P CzD : ω I ě d{2u Y tω P D : ω I ď d{2u. Corollary 2.7 presents several general properties of the enclosure W Ω pt q. In particular iii) iv) show that Π β and Π α determine the the sign of ˆβpωq and of ˆαpωq, where ˆβ and ˆα are defined in Proposition 2.6. Corollary 2.7. Let W Ω pt q denote the enclosure (2.8) and denote by D the disk (2.12). Let ˆα and ˆβ be the functions defined in (2.13) and in (2.14), respectively. Let Π β and Π α denote the sets (2.17). Then the following properties hold: i) If ω R ir, then ω R W Ω pt q provided that ω I is large enough. ii) For sequences tω n u P W Ω pt q, with ωr n Ñ 8, n Ñ 8, it holds that ˆβpω n q 2ωI n pωn R q2 {d and ωi n Oppωn R q 2 q. iii) If ω R ir Y BD then ˆβpωq ě if and only if ω P Π β. iv) If ω R ir Y BD then ˆαpωq ě if and only if ω P Π α. Proof. i) The value ˆβpωq gets arbitrary large as ω I Ñ 8 but W pbq is bounded. ii) Assume tω n u P W Ω pt q, ωr n Ñ 8, then i) implies that ωn I is bounded and ˆβpω n q 2ωI n pωn R q2 {d from (2.15). Hence, the boundedness of ˆβpω n q yields that ωi n Oppωn R q 2 q. iii) and iv) follow by straightforward calculations. Lemma 2.8. The functions r n in (2.4) have the following properties: i) For given ω P CzpiR Y tδ`, δ, 8uq there is a unique pair pα, βq P R 2 such that r n pα, βq ω for some n P t1, 2, 3, 4u. Further, if r m pα, βq ω, m n, then α c, β d 2 {4, and ω ac d 2 {16 id{4.

6 6 CHRISTIAN ENGSTRÖM AND AXEL TORSHAGE ii) For given ω P irztδ`, δ u and β P R, the unique α P R such that r n pα, βq ω for some n P t1, 2, 3, 4u is α ωi 2 c ` dω I ` ωi 2 β. iii) For given ω P irztu and α P R, the unique β P R such that r n pα, βq ω for some n P t1, 2, 3, 4u is ˆ β pc ` dω I ` ωiq 2 1 ` 1 ωi 2 α. Proof. i) Proposition 2.6 yields that ω P W Ω pt q if and only if pˆαpωq, ˆβpωqq P Ω. Thus r m pα, βq ω is only possible for pα, βq pˆαpωq, ˆβpωqq. Assume ω r n pα, βq r m pα, βq, n m. Then, ω is also a double root since ω R ir and roots of p pα,βq are symmetric with respect to the imaginary axis. The result is then obtained using an ansatz with these two double roots. ii) iii) Follows trivially from the definition of p α,β. In Theorem 2.9, we show that the enclosure of the numerical range W Ω pt q is closely related to the set 4ď (2.18) W BΩ pt q r n pbωq. Theorem 2.9. Let W Ω pt q denote the enclosure (2.8) and let W BΩ pt q denote the set (2.18). Then the following equalities hold: i) W Ω pt q X ir W BΩ pt q X ir. ii) BW Ω pt qzir W BΩ pt qzir. n 1 Proof. i) The inclusion W Ω pt q X ir Ą W BΩ pt q X ir is clear from (2.8) and (2.18). Let ω P W Ω pt q X ir, then the result follows from Corollary 2.3 and Proposition 2.4. ii) Assume δ` P W Ω pt qzir, then δ` P BW Ω pt q follows from Lemma 2.5 and Corollary 2.3 implies δ` P W BΩ pt q. The proof for δ is similar to the proof for δ` and for 8 the result follows directly. Apart from ir Y tδ`, δ, 8u, Lemma 2.8 i) yields that r n : R ˆ R Ñ C» R 2 is injective. Then Br n pωqzir r n pbωqzir is a consequence of the invariance of domain theorem [Bro12]. Hence, BW Ω pt qzir Ă W BΩ pt qzir and ii) follows from Lemma 2.8 i). Corollary 2.1. The boundary of W Ω pt qzir is W BΩ pt qzir. Definition Let N : H for d ă 2? c and N : rδ, δ`s for d ě 2? c. Define the sets τ 1 : tinf W paq, inf W pbqu Y tsup W paq, sup W pbqu, τ 2 : tinf W paq, sup W pbqu Y tsup W paq, inf W pbqu, Ť4 Ť4 R 1 : n 1 r npτ 1 q X irzn, R 2 : n 1 r npτ 2 q X N. Let m : R 1 9YR 2 Ñ N denote a counting function, where for iµ P R j we set mpiµq : n 1 ω 2 I 4ÿ #tτ P τ j : r n pτq iµu. Due to continuity Y n r n p 8, βq tδ, i8u and Y n r n p8, βq tδ, 8u.

7 NUMERICAL RANGE OF NONSELFADJOINT RATIONAL OPERATOR FUNCTIONS 7 Proposition Let W BΩ pt q denote (2.18), and let R 1, R 2, τ 1, τ 2, and m be defined as in Definition Assume that c ą, then iµ P W BΩ pt q is an endpoint of a line segment of W BΩ pt q X ir if and only if iµ P R 1 9YR 2 and mpiµq is odd. Further, if iµ is an isolated point of W BΩ pt q X ir, then iµ P R 1 9YR 2 and mpiµq is even. Proof. The result is first shown for iµ R t, δ`, δ, i8u. Assume iµ R N Yt, i8u is an endpoint of a line segment or an isolated point of W BΩ pt q X ir. Then α ` µ 2 ` µ 2 c ` dµ ` µ 2 β, µ for some pα, βq P BΩ. Assume that pα, βq R τ 1, then since 2 c`dµ`µ ą it follows by 2 similar arguments as given in Proposition 2.4 that there exist a pair pα 1, β 1 q P ΩzBΩ such that α 1 ` µ 2 ` µ2 c`dµ`µ β 1. Then, Lemma 2.8 ii) iii) gives a contradiction. 2 Hence pα, βq P τ 1 and iµ P R 1. Assume that iµ is an isolated point, then from the symmetry of the roots with respect to the imaginary axis it follows that mpiµq is even. Assume that iµ is an endpoint of a line segment. From the injectivity proven in Lemma 2.8 ii) iii) follows then that exactly one root must be on the line segment. Thus from the roots symmetry with respect to the imaginary axis it follows that mpiµq is odd. If iµ P N zt, δ`, δ u a similar argument proves the claim for R 2. For the converse, assume iµ P R 1 and mpiµq odd. Then since 2 c`dµ`µ µ ą it follows that iµ is 2 the root for an unique pair pα, βq P τ 1. Assume that iµ is not the endpoint of a line segment, then since it is not an isolated point it is an inner point of a line segment in W BΩ pt q X ir. From injectivity proven in Lemma 2.8 ii) iii), symmetry with respect to the imaginary axis, and that mpiµq is odd, it follows that for pα 1, β 1 q P BΩ sufficiently close to pα, βq there is exactly one simple root on the imaginary axis that is in the vicinity of iµ. Take points iµ 1, iµ 2 in the vicinity of iµ such that, µ µ 1 ă µ ă µ 2 and 2 i c`dµ i`µ ą. Then there is some pα 2 1, β 1 q, pα 2, β 2 q P BΩ such i that µ 2 1 α 1 ` µ 2 1 ` c ` dµ 1 ` µ 2 β 1, 1 µ 2 i c`dµ i`µ 2 i µ 2 2 α 2 ` µ 2 2 ` c ` dµ 2 ` µ 2 β 2. 2 Since ą there is a line of solutions pα, βq intersecting BΩ twice. Hence we can assume that α 1 α 2 α. By continuity there must exist a β 3 between β 1 and β 2 such that pα, β 3 q has the root iµ. But β 3 β which contradicts Lemma 2.8 iii). The proof for iµ P R 2 and mpiµq odd is similar. Assume iµ P t, δ`, δ u, then the result is shown by investigating each case for iµ P R 1 9YR 2 and when iµ is an endpoint of a line segment of W BΩ pt q X ir. Remark If c the point µ is always a solution to (2.3) and similar results as in Proposition 2.12 can for this case be obtained from the reduced cubic polynomial. Proposition Let W BΩ pt q denote (2.18), and let R 1, R 2, and m be defined as in Definition Then W BΩ pt q X ir is obtained from R 1 9YR 2 by the following algorithm: (1) Set I : tiµ P R 1 9YR 2 : mpiµq is oddu and enumerate µ P I increasingly µ 1 ă µ 2 ă....

8 8 CHRISTIAN ENGSTRÖM AND AXEL TORSHAGE (2) Add an interval between iµ j, iµ j`1 in I if j is odd. (3) Set W BΩ pt q X ir I Y pr 1 9YR 2 q. Proof. From Proposition 2.12 it follows that step 1 defines I as the set of endpoints of line segments of R 1 9YR 2, where iµ 1 is the minimal imaginary part of a line segment. Then iµ 2 must be the endpoint of that segment, which is the point with maximum imaginary part. Doing this iteratively gives that for all odd j, iµ j is the minimal imaginary part of a line segment and for even j, iµ j is the maximal imaginary part of a line segment. Hence, step 2 sets I to W BΩ pt q X ir apart from isolated points. These points are added in step 3. The following lemma and Lemma 2.8 implies that W BΩ pt qzir has a finite number of points where more than one curve component intersect. Lemma The roots of the polynomial p α,β defined by (2.3) have the following properties: i) Fix α P Rztu, then p pα, q has a multiple root for at most 4 values β P R. ii) Fix β P Rztu, then p p,βq has a multiple root for at most 5 values α P R. iii) p p,βq has a double root at and the roots aβ ` c d 2 {4 id{2. iv) p pα,q has the roots?α and ac d 2 {4 id{2. Proof. If α or both β and d 2? c, then the discriminant ppα,βq is zero and p pα,βq has a double root. For all other cases, we conclude from definition that ppα,βq is a fifth-degree polynomial in α and a fourth-degree polynomial in β. Definition Two disjoint sets Γ 1, Γ 2 P C are neighbors if BΓ 1 X BΓ 2 contains at least one curve segment. The algorithm presented in Proposition 2.17 is described in Figure 1. Proposition Let W Ω pt q denote the enclosure (2.8) and let W BΩ pt q denote (2.18). Then W Ω pt qzir W BΩ pt qzir if W paq or W pbq is constant. Otherwise W Ω pt qzir is obtained from W BΩ pt qzir by the following algorithm: (1) Let O be the component of CzpW BΩ pt qzirq containing values of ω with arbitrarily large imaginary parts. (2) Let I Ă CzpW BΩ pt qzirq be the union of all components that are neighbors of O. (3) Let O Ă CzpW BΩ pt qzirq be the union of all components that are neighbors of I. (4) If I Y O CzpW BΩ pt qzirq, go to step 2. (5) Set W Ω pt qzir I Y W BΩ pt qzir. Proof. If W paq or W pbq are constant the result follows by definition. Corollary 2.7 shows that only one component of CzpW BΩ pt qzirq contains values of ω with arbitrarily large imaginary parts. Hence the initial set O Ă CzW Ω pt qzir in Step 1 is well-defined. From Lemma 2.8 and Lemma 2.15 it follows that W BΩ pt qzir has a finite number of points with more than one curve component intersecting it and by definition W BΩ pt q has at most 4 components. Hence, CzpW BΩ pt qzirq consists of a finite number of components, which implies that the algorithm will terminate after a finite number of steps. Corollary 2.1 yields that the set W BΩ pt qzir is the boundary of a closed set and if two components of CzpW BΩ pt qzirq are neighbors, one is a

9 NUMERICAL RANGE OF NONSELFADJOINT RATIONAL OPERATOR FUNCTIONS Figure 1. Visualization of how WΩ pt qzir is obtained from WBΩ pt q using the algorithm in Proposition Red and blue denotes points given by α P r, 8s and α P r 8, q, respectively. subset of WΩ pt qzir, and one is a subset of CzWΩ pt qzir. Thus the algorithm gives the sets I Ă WΩ pt qzir and O Ă CzWΩ pt qzir, and the Proposition follows therefore from the termination criteria. Remark In graph theory the algorithm in Proposition 2.17 is related to the 2-colorability of the dual graph of WBΩ pt qzir, [SK86, Theorem 2-3]. 3. Analysis of the enclosure of the numerical range In this section, the boundary of the enclosure is analyzed in detail. We derive conditions for the existence of a strip in the complement of the numerical range and prove properties of that strip.? The function T pλq : T p λq is analytic in the upper half-plane C` and ImpT pλqu, uq ě, for λ P C`, if and only if inf W pbq ě. Hence, T is a Nevanlinna function if and only if inf W pbq ě. Since operator functions with applications in physics often are Nevanlinna functions [AL95, ELT17], we analyze in this section the enclosure WΩ pt qzir under the assumption inf W pbq ě. However, the analysis when inf W pbq is allowed to be negative is similar. Let ω1, ω2, ω3, ω4 be the roots of pα,β as defined in (2.3) and define (3.1) it1 : ω1 ` ω2, it2 : ω3 ` ω4, v1 : ω1 ω2, v2 : ω3 ω4. From the relation between the coefficients and roots of a polynomial it follows that (3.2) t1 ` t2 t 1 t 2 ` v1 ` v2 t 1 v2 ` t 2 v1 v 1 v2 d, α ` β ` c, αd, αc. It is of interest to see when pα,β has purely imaginary solutions and the following result shows that it depends on the sign of α.

10 1 CHRISTIAN ENGSTRÖM AND AXEL TORSHAGE Lemma 3.1. Let p α,β be defined as in (2.3). Then, the following statements hold for the roots on the imaginary axis: i) If α ă, then p α,β has at least two roots of the form iµ, µ P R, where µ ą for exactly one root and µ ď (µ ă if c ą ) for at least one root. ii) If α ą, then all roots of p α,β of the form iµ, µ P R, satisfies µ ď (µ ă if c ą ). If d ă 2? c there are no purely imaginary root and if d ě 2? β ` c there are at least two purely imaginary roots. Proof. i) If β, the result follows from Lemma Assume β ą and that µ is a root of the real function ˆp α,β defined by ˆp α,β pµq : p pα,βq piµq. Then, (3.3) ˆp α,β pµq pµ 2 ` αqpµ 2 ` dµ ` cq ` βµ 2, where ˆp α,β is positive and of even order. For α ă it follows that ˆp α,β pq ď (with equality if and only if c ) and thus there is a positive and a non-positive root (negative if c ą ). There can be no other roots µ ą since ˆp 1 α,βpq αd ă and ˆp 1 α,β is convex on r, 8q. Hence, ˆp α,βpµq for exactly one value µ ą. ii) If α ą then ˆp α,β pµq ą whenever µ ą or d ă 2? c. Assume α ą, d ě 2? β ` c, then ˆp α,β p d{2q ď, which implies that ˆp α,β has at least two real roots. The set W BΩ pt q is given by the values on the rectangle BΩ. Hence, there are two types of curves that are interesting to analyze. In subsection 3.1, β P W pbq is fixed and in subsection 3.2, α P W paq is fixed Variation of the numerical range W paq. The set W BΩ pt q defined in (2.18) was in Proposition 2.17 used to determine the enclosure W Ω pt q. In this section, we will describe the subset of W BΩ pt q obtained by fixing β and varying α P W paq in greater detail. To this end we consider the set (3.4) W Rˆtβu pt q 4ď r n pr ˆ tβuq, n 1 defined according to (2.6). Note that for ω P CzpiR Y tδ`, δ uq, the point ω is in W Rˆtβu pt q if and only if β ˆβpωq, where ˆβpωq is defined in (2.13). The set W Rˆtβu pt q can in the variable α P W paq be parametrized into a union of four curves. For β, the set W Rˆtβu pt q is completely characterized by Lemma 2.15 and we will therefore assume β ą in the rest of Section 3.1. Figure 2 illustrates possible behaviors of W Rˆtβu pt q. Proposition 3.2. Let W Rˆtβu pt q denote the set (3.4) and take β, c ą. Then iµ P W Rˆtβu pt qzir for µ P R if and only if µ or µ is a real solution to ˆβ (3.5) q β pµq : µ 4 ` 2dµ 3 ` p2c ` d 2 qµ 2 ` d 2 ` 2c µ ` c pβ ` cq. The statement above holds also if c, with the exception that zero is not in the set W Rˆtβu pt qzir.

11 NUMERICAL RANGE OF NONSELFADJOINT RATIONAL OPERATOR FUNCTIONS (a) (b) (c) Figure 2. Examples of the set W Rˆtβu pt q Ă C, where red and blue denotes W r,8sˆtβu pt q and W r 8,qˆtβu pt q, respectively. Proof. Assume that iµ for some µ P R is a root of p pα,βq of order greater than one. Hence, iµ has to be at least a double root and we set t 1 2µ, v 1 µ 2 in (3.1). The system (3.2) can then be written as (3.6) 2µ ` t 2 d, 2µt 2 ` µ 2 ` v 2 α ` β ` c, 2µv 2 ` µ 2 t 2 αd, µ 2 v 2 αc. Solving (3.6) shows that µ is a solution if and only if µ or µ is a root of (3.5). Assume iµ P W Rˆtβu pt qzir for some µ P R. Then, by the symmetry with respect to the imaginary axis, iµ is a double root of p pα,βq for some α. Hence µ or µ is a real solution to (3.5). For the converse, assume that one of the poles is purely imaginary and that iµ P tδ`, δ u is a root of q β. Then it follows from (3.5) that c P t, d 2 {4u. Furthermore, for c d 2 {4, µ d{2 is a solution to (3.5). For c d 2 {4, δ` δ id{2 is in the limit α Ñ 8 a root of p pα,βq, but p pα,βq do not have a purely imaginary root for any α P R`, which implies iµ P W Rˆtβu pt qzir. For c is µ a solution to (3.5), and thus a double root of p pα,βq for some α P R. Moreover, zero is a root of p pα,βq for all α P R. Hence, the symmetry with respect to the imaginary axis implies that zero only can belong to the set W Rˆtβu pt qzir if for some α, zero is a triple root of p pα,βq. An ansatz with a triple root implies β, which yields a contradiction. Now assume that µ, or µ is a real solution to (3.5), and iµ R tδ`, δ u. Then iµ is a double root of p pα,βq for some α P R and Lemma 2.8 ii) yields iµ P W Rˆtβu pt qzir. Let qβ denote the discriminant of q β. Corollary 3.3. Let p α,β and q β denote the polynomials (2.3) and (3.5), respectively. Then p α,β has a root iµ, µ P Rztu of multiplicity n ą 1 for some α P R if and only if µ is a root of q β of multiplicity n 1. Proof. The assumptions on the coefficients imply that q β can not have an quadruple root. A straightforward calculation shows that p p,βq has a quadruple root iµ for some α if and only if µ is a triple root of q β. Assume that p p,βq has a triple root iµ for some α. From (3.1) and (3.2) follows that this assumption is equivalent to qβ. Hence, q β has a multiple root and the multiplicity must be two. For

12 12 CHRISTIAN ENGSTRÖM AND AXEL TORSHAGE n 2, we showed that the multiplicity of a root of q β can not be larger than one. Then the result follows from Proposition 3.2. The multiplicity of a real root µ of q β determine the number of segments of W Rˆtβu pt qzir intersecting iµ, (if c there is no intersection in zero). For convenience, we set in Lemma 3.4 some constants to 8. These constants are used in Proposition 3.5. Lemma 3.4. Let q β be the polynomial (3.5) and qβ its discriminant. Then the following properties hold: i) For β ă 4c, qβ has an unique non-negative solution d 1 P p, 2? cq. Set d 2 d 3 8. ii) For 4c ď β ă 8c, qβ has the three non-negative solutions d 1 P p, 2? cq, d 2 P p2? c, 2? βs, and d 3 P r2? β, 8q. iii) For β ě 8c ą, qβ has two non-negative solutions d 1 P p, 2? cq, d 2 P p2? c, 2? βs. Set d 3 8. iv) For c, qβ has two non-negative solutions d 1, d 2 27β{32. Set d 3 8. v) The polynomial q β has zero real roots if d ă d 1 and four real roots if d 2 ď d ď d 3. In all other cases q β has two real roots. Proof. Let ˆd : d 2 {4 and consider fp ˆdq : qβ as a polynomial in ˆd, where each positive root will correspond to exactly one positive solution d. By the definition of the discriminant, we obtain fp ˆdq ˆ27 (3.7) 32β 2 pβ 8cq ˆd 3 32 β2 6βc 24c 2 ˆd 2 3c 2 p5β ` 8cq ˆd ` 8c 3 pβ ` cq. iv) For c the roots are, 27β{32 and the result follows. If c ą, the existence of a root ˆd 1 P p, cq follows from fpq ą, fpcq ă. The discriminant of f is f β 12 c 3 pβ 4cq 3. i) Assume β ă 4c, then f ă and thus f has only one real root. ii) Assume 4c ď β ă 8c, then f ě and f is a cubic polynomial. It can be seen that fpβq ě and fp ˆdq Ñ 8, ˆd Ñ 8. Hence there is one root ˆd2 in pc, βs and one root ˆd 3 in rβ, 8q. iii) Assume β ě 8c. Then fpβq ą, thus there is a root ˆd 2 in pc, βs. In the special case β 8c, f is a quadratic polynomial and thus there are no more roots. Otherwise β ą 8c and then the last root will be negative. v) The sign of f will be negative if and only if d 1 ă d ă d 2 or d ą d 3 and thus in these cases q β has two roots. In all other cases it will either have zero or four roots. When d, q β has no roots and by continuity q β has no roots for d ă d 1. For d 2 ď d ď d 3 it holds that d ą 2? c and q β p d{2 ad 2 {4 cq ă. Then, since the highest order term of q β is positive it must have at least one root and thus four roots. Figures 3.(a) and 3.(c) depict W Rˆtβu pt q for d d 1 and d d 2 d 3, respectively. In Figure 3.(c) the set W Rˆtβu pt q intersects the imaginary axis three times at 1. This can only happen when d 2? β 4? c and µ d{4, which implies that all sets W Rˆtβu pt q with this property are linear scalings of the case presented in Figure 3.(c).

13 NUMERICAL RANGE OF NONSELFADJOINT RATIONAL OPERATOR FUNCTIONS (a) (b) (c) Figure 3. Examples of the set W Rˆtβu pt q Ă C, where red and blue denotes W r,8sˆtβu pt q and W r 8,qˆtβu pt q, respectively. In (a) there is a double root of the function q β in (3.5). In (b) there are two distinct roots. In (c) there is one distinct root and one triple root. Proposition 3.5. Let d k P R, k 1, 2, 3 denote the constants defined in Lemma 3.4 and set I 1 : p d{2 a2d 2 8cq{2, d{2q, l 3 : p d{2 ad 2 4c{2, d{2q, I 2 : p 8, d{2 a2d 2 8c{2q, I 4 : p 8, d{2 ad 2 4c{2q. For c ą, the set W Rˆtβu pt qzir defined as in (3.4) intersects the imaginary axis at zero and in the following points, counting multiplicity: i) For d ă 2? c, there are no intersections if d ă d 1. If d ě d 1 there are two intersections in the interval p 8, d{2q. ii) For d 2? c, the two intersections are d{2 and pd ` 3? 4βdq{2. iii) For 2? c ă d ă 2? β ` c there is one intersection in I 2. Additionally, if d 2 ď d ď d 3 there are three intersections in p d{2, q, and if d ě d 3 one intersection in p d{2, q. iv) For d 2? β ` c, d{2 is an intersection and there is one intersection in the interval I 2. If d ě d 3 there are two further intersections in the interval p d{2, q, and if d ă d 3 there are no further intersections. v) For d ą 2? β ` c, there is one intersection in I 4, and one in p 8, d{2qzi 4. If d ě d 3 there are two further intersections in the interval p d{2, q, and if d ă d 3 there are no further intersections. Proof. Assume c ą, then the intersections will coincide with roots of q β. i) From Lemma 3.4 iv) follows that there are no intersections when d P p, d 1 q and two intersections for d P rd 1, 2? cq. Since q β pµq ą for µ ě d{2, d P rd 1, 2? cq the two intersections are in p 8, d{2q. ii) Follows from straight forward computations. iii) The value q β pµq is negative on p 8, d{2qzi 2, thus there are no intersections in p 8, d{2qzi 2. The function q β is convex on I 2 and q β p 8q ą. Hence, there is one intersection in I 2. All other intersections are in p d{2, q and the number of intersections is given by Lemma 3.4 iv). iv) A straight forward computation show that qp d{2q and the remaining statements follows as in iii).

14 14 CHRISTIAN ENGSTRÖM AND AXEL TORSHAGE v) We have qp d{2q ą, q 1 pµq ą for µ P I 3, and qp pd `?d 2 4cq{2q ă, thus one root of q is in I 3. The function q β is convex on I 4 and q β p 8q ą. Hence, there is one intersection in I 4. All other intersections are in p d{2, q and the number of intersections is given by Lemma 3.4 v). Remark 3.6. The statements of Proposition 3.5 hold also in the case c except that there is no intersection in. Proposition 3.7. Let W Rˆtβu pt qzir be defined as in (3.4). Then the point ω is in W Rˆtβu pt qzir if and only if ω P t, 8u (ω 8 if c ) or ω I and one of the following four equations hold: c b (3.8) ω R P ωi Pω 2 I Q ωi ωi 2, where (3.9) P ωi : c d2 2 dω I βd, 4ω I Q ωi : βc ` c 2 ` 2cdω I ` 4cωI 2. Proof. Due to the symmetry with respect to the imaginary axis we can choose t 1 2ω I in (3.1). The result then follows from straight forward computations, where (3.2) is used. The system (3.2) can be solved for a given α by computing the roots of a fourth order polynomial. However, Proposition 3.7 shows that if ω I is known, ω R can be computed independently of α. Hence, for ω P W Rˆtβu pt qzir this allows us to regard ω R as a multivalued function in ω I. By a horizontal strip S Ă C we mean an open set of the form (3.1) S tω P C : s ă ω I ă s 1 u, where s, s 1 P R and s ă s 1. Definition 3.8. For a closed set Γ Ă C, a horizontal strip S Ă CzΓ as defined in (3.1) is said to be maximal with respect to Γ if (3.11) Γ X pr ` is q Ø, Γ X pr ` is 1 q Ø. The set Γ X pr ` is q is called the local minimum points and Γ X pr ` is 1 q is called the local maximum points. Definition 3.9. For a closed set Γ Ă tω P C : ω I ă s P Ru assume that there are n maximal horizontal strips, S 1,..., S n, with respect to Γ. Further let p min and p max denote the points in Γ with the smallest respectively largest imaginary values in Γ. Define the set nď (3.12) M : Γ X S i Y pp min Y p max q, where the points in M will be called the extreme points of Γ. i 1 In the following a strip is always assumed to be horizontal and maximal with respect to a given set. Figure 2.(c) depicts W Rˆtβu pt qzir for a case with a strip S as defined in (3.11). Note that the point in W Rˆtβu pt qzir with largest imaginary part is always zero.

15 NUMERICAL RANGE OF NONSELFADJOINT RATIONAL OPERATOR FUNCTIONS 15 Corollary 3.1. The smallest imaginary part for a point in W Rˆtβu pt qzir as defined in (3.4) is less than Impδ q. Proof. The claim follows immediately since there exists an ω with ω I Impδ q, ω R satisfying (3.8). Lemma Let W Rˆtβu pt qzir denote the set in (3.4) and let P ωi and Q ωi denote the expressions in (3.9). A point ω P W Rˆtβu pt qzir is an extreme point in the sense of Definition 3.9 if and only if ω I is a distinct root of fpω I q : ωi 2 pp ω 2 I Q ωi q. The roots of f are (3.13) i) ω I d ad b 2 4β, ii) ω I d d 1 ` 4β 4c d A double root of f is only possible if d 2? β ă 4? c, where ω ac d 2 {16 id{4. Assume that ω I is a double root of f. Then, ω is a point where more than one curve component intersects W Rˆtβu pt qzir. Proof. By simple computations it follows that the roots of f are (3.13). Lemma 2.8 i) and Proposition 3.7 imply that a double root µ of f exists if and only if µ d{4?β{2 ă? c. Then, Proposition 3.7 yields that the corresponding points on W Rˆtβu pt q are ω ac d{16 id{4. Pω 2 I Q ωi is non-negative in a neighborhood of d{4, thus it is not an extreme point but a point where more than one curve component intersects. Assume that ω P W Rˆtβu pt qzir is an extreme point and that fpω I q. Then since one of the equations (3.8) hold it follows that Pω 2 I Q ωi ą. Likewise P ωi ap ω 2 I Q ωi ωi 2 ą, since ω R. The function P ωi ap ω 2 I Q ωi ωi 2 is continuous b in ω I. Hence, there exists an open interval I ωi containing ω I such that P µ Pµ 2 Q µ µ 2 ą for each µ P I ωi. Then (3.8) holds for any point in the interval, which contradicts that ω is an extreme point. Hence it follows that ω I is a distinct root of f. Now suppose ω I is a distinct root of f. Then for every open interval I ωi containing ω I, there exists an µ P I ωi such that fpµq ă, and it is thus an extreme point by (3.8). Figure 2.(c) shows a case where the local extreme points as well as the points with smallest imaginary part are not on the imaginary axis. Lemma Let W Rˆtβu pt qzir and P ωi, Q ωi be defined as in (3.4) and in (3.9), respectively. A point iµ P ir, where µ b P Rztu is an extreme point (3.12) to W Rˆtβu pt qzir if and only if P µ ` Pµ 2 Q µ µ 2 and µ is a distinct intersection of the imaginary axis. Proof. From Proposition 3.2 follows that each intersection of the imaginary axis is equivalent to a root of q β as defined in (3.5). Assume iµ P W Rˆtβu pt qzir is an extreme point. Then by Proposition 3.7, Pµ 2 Q µ ě and one of the equations b P µ Pµ 2 Q µ µ 2 hold. If Pµ 2 Q µ, then by continuity iµ is for some α a quadruple root of the polynomial p pα,βq defined in (2.3). Hence, by Proposition 3.7 the root can not be an extreme point. Furthermore, Corollary 3.3 implies that µ is a triple root of q β, thus not distinct. Assume that P 2 µ Q µ ą and

16 16 CHRISTIAN ENGSTRÖM AND AXEL TORSHAGE b P µ ` Pµ 2 Q µ µ 2, then it follows from Lemma 3.11 that iµ P W Rˆtβu pt qzir b is not an extreme point. Hence, we have shown that P µ ` Pµ 2 Q µ µ 2 and Pµ 2 Q µ ą. Assume µ is not a distinct root of q β, then at least two segments b of W Rˆtβu pt qzir intersect iµ. Since P µ Pµ 2 Q µ µ 2 ă, Proposition 3.7 implies that in some interval containing µ there is for a given ω I at most two solutions ω. Combining these results shows that iµ bis not an extreme point and the intersection must then be distinct. Assume P µ ` Pµ 2 Q µ µ 2 and that µ is a distinct root of q β, then there is only one segment of W Rˆtβu pt qzir intersecting b iµ. Furthermore, Pµ 2 Q µ ą and thus P µ Pµ 2 Q µ µ 2 ă. Proposition 3.7 implies that iµ is an extreme point. Proposition Let W Rˆtβu pt q denote the set (3.4). Then every real number is the real part of some point in the set W Rˆtβu pt q. Proof. For ω R ac d 2 {4 note that δ`, δ P W Rˆtβu pt q. In all other cases, equation (2.13) has a solution ω I for given β ˆβpωq and ω R. Then α : ˆαpωq is uniquely given by (2.14). Hence, p pα,βq pωq in (2.3) has for fixed ω R and β a solution for some ω P C, and α P R. Lemma Let W Rˆtβu pt qzir be defined as in (3.4). Then, for each bounded component γ Ă W Rˆtβu pt qzir and α P R either one root of the polynomial p α,β in (2.3) belongs to γ, or one root of p α,β can be written as iµ, for µ P J : rminpγzirq I, maxpγzirq I s. Proof. If the bounded component of W Rˆtβu pt qzir does not intersect ir it contains by continuity a root for all α P R. If the bounded component γ Ă W Rˆtβu pt qzir intersects with ir the curve is closed with an even number of intersections of the imaginary axis (counting multiplicity in (3.5)). Assume that there is an α P R such that no root is on γ, and none of the roots are purely imaginary with imaginary part in J. Then, by definition (3.4), r n pα, βq is not on γ Y ij for n 1, 2, 3, 4. From the continuity of the roots it follows that there exists an α 1 such that r n pα 1, βq P γ Y ij and for a sufficiently small ɛ, the roots r n pα 1 `ɛ, βq, n 1, 2, 3, 4 are not on γ YiJ. Since this can only happen on the imaginary axis it follows that (2.3) for α 1 has a purely imaginary multiple root iµ on γ. Thus for some ordering of the roots r 1 pα 1, βq r 2 pα 1, βq iµ. Since r 1 pα 1 ` ɛ, βq R γ and r 2 pα 1 ` ɛ, βq R γ it follows by continuity that they are imaginary. From Lemma 2.8 ii) follows that one of the roots has a larger imaginary part and one has a smaller imaginary part than µ. Then if both are outside J, it follows that J will consist of only one point. Hence µ is a double root of (3.5) and from Corollary 3.3, iµ is a triple root of p pα,βq. Hence r 3 pα 1, βq iµ and by injectivity (Lemma 2.8) one root will belong to γ for α 1 ` ɛ, which is a contradiction. Definition Denote the k real roots of q β defined in (3.5) by (3.14) µ 1 ď µ 2 ď... ď µ k, where k might be zero.

17 NUMERICAL RANGE OF NONSELFADJOINT RATIONAL OPERATOR FUNCTIONS 17 Proposition Let W Rˆtβu pt qzir be defined as in (3.4) and let d 2 denote the constant defined in Lemma 3.4. Then, there is a unique maximal strip S with respect to W Rˆtβu pt qzir, as in Definition 3.8, if and only if d ą minp2? β, d 2 q. If d ď minp2? β, d 2 q, there is no such strip. Proof. By definition (3.4) the imaginary part of W Rˆtβu pt q is bounded and there is a strip maximal with respect to W Rˆtβu pt qzir if and only if there are at least three extreme points with different imaginary parts (disregarding the point ). If d ą 2? β it follows from Proposition 3.7 that ω R W Rˆtβu pt qzir for ω I d{4 and Corollary 3.1 implies that d{4 is not the least imaginary part. Hence, there exists a strip that is maximal with respect to W Rˆtβu pt qzir containing the set R id{4. If 2? β ě d ą d 2, it follows that (3.5) has four distinct real roots (if d 3 2? β it follows from Lemma 3.4 that d 2 2? β and this contradicts 2? β ě d ą d 2 ). Thus, it follows from Corollary 3.3 that there are five intersections of W Rˆtβu pt qzir with the imaginary axis. Assume that the point with smallest imaginary part is not on the imaginary axis. Then it follows from Lemma 3.11 and d 2 ą 2? c that the least imaginary part is larger than d{2, which contradicts Corollary 3.1. Hence the point with the smallest imaginary b part is on the imaginary axis and thus given by the root µ 1, where P µ1 ` Pµ 2 1 Q µ1 µ 2 1 by Lemma Assume that there is no strip maximal with respect to W Rˆtβu pt qzir. Then Lemma 3.12 implies b P µ ` Pµ 2 Q µ µ 2 ą, with Pµ 2 Q µ ě for all µ P pµ 1, q. b Define for µ P pµ 1, q the function fpµq : P µ Pµ 2 Q µ µ 2. Then Proposition 3.7 implies fpµ i q for i 2, 3, 4. Take i P t2, 3, 4u and assume that fpµq is either positive or negative in an open punctured interval around µ i. Then it follows from Corollary 3.3 that µ i is not a distinct root of q β. Hence fpµq alternates signs between the roots. Proposition 3.7 implies that there must be two bounded components of W Rˆtβu pt qzir. Since d ă 2? β ` c, Proposition 3.5 implies that µ i ą d{2 for i 2, 3, 4. Then from Lemma 3.14 it follows that both poles are larger than d{2, which gives a contradiction and a strip maximal with respect to W Rˆtβu pt qzir must therefore exist. Assume that there are at least two strips maximal with respect to W Rˆtβu pt qzir, then it must be at least three components of W Rˆtβu pt qzir, one is unbounded and two are bounded. Lemma 3.14 implies that the bounded components will both enclose a pole. This yields that the poles are imaginary and thus each bounded component intersects the imaginary axis twice and we have four real roots of (3.5). Since the imaginary parts of these roots approach as the real parts approach 8, the points and 8 will be in the same component with no other intersections of the imaginary axis. This means that for α ě there are two roots in the unbounded component. By Lemma 3.14 there are always at least one root on or enclosed by a bounded component. Hence, for α ě there is only one root in each bounded component. Thus due to symmetry the roots in the bounded components are imaginary for all α ě. Hence, if ω belongs to a bounded component of W Rˆtβu pt qzir then ω I ď d{2 and thus all the solutions µ to (3.5) satisfy µ ď d{2, which contradicts Proposition 3.5.

18 18 CHRISTIAN ENGSTRÖM AND AXEL TORSHAGE Assume d ď minp2? β, d 2 q and that there is a strip maximal with respect to W Rˆtβu pt qzir. Then W Rˆtβu pt qzir has at least three extreme points. It follows that (3.5) has at most two distinct roots, µ 1, µ 2, and thus at most two extreme point on ir. Hence, there has to be an extreme point in CziR. This implies that d ă 2? c since otherwise none of the constants in (3.13) can be real and negative. Then since d ă 2? c and d ď minp2? β, d 2 q, Lemma 3.11 and Proposition 3.5 yield that all possible extreme points have imaginary parts smaller than d{2, which is the imaginary part of the poles. Hence by Lemma 3.14 there is no strip maximal with respect to W Rˆtβu pt qzir. Figure 4 depicts W Rˆtβu pt q and illustrates the claim of Proposition The following proposition gives a detailed description of the strip, S, maximal with respect to W Rˆtβu pt qzir. If there is a strip maximal with respect to W Rˆtβu pt qzir, we let s denote the local minimum and s 1 dentote the local maximum, as defined in Definition 3.8. Moreover, M denotes the set of extreme points (3.12) to W Rˆtβu pt qzir. Proposition Let W Rˆtβu pt qzir be defined as in (3.4). Given the ordering of the roots µ i in Definition 3.15, the following properties hold: i) If β ă 4c there is a unique strip maximal with respect to W Rˆtβu pt qzir if and only if d ą 2? β. If d ă β`4c 2? then the local maximum points are not on the imaginary c axis and s p d ad 2 4βq{4. If d ě β`4c 2? the local maximum c point is iµ 2. The local minimum points are not on the imaginary axis and s 1 p d ` ad 2 4βq{4. ii) If β ě 4c ą there is a unique strip maximal with respect to W Rˆtβu pt qzir if and only if d 2 ă d. The local maximum point is iµ 2. If d ď β`4c 2? then the local minimum point is iµ c 3. If d ą β`4c 2? the c local a minimum points are not on the imaginary axis and s 1 p d ` d2 4βq{4. Proof. We will first show that s µ 2 if the local maximum point is on the imaginary axis and s p d ad 2 4βq{4 if the local maximum is not on the imaginary axis. If a strip maximal with respect to W Rˆtβu pt qzir exists then it follows from Corollary 2.7 that there is exactly one unbounded component and one bounded component of W Rˆtβu pt qzir. By continuity the bounded component intersects the imaginary axis an even number of times. If the local maximum is on the imaginary axis it must thus be the largest root of the bounded component. Hence, the local maximum is the root µ 2 in (3.15) if the bounded component intersects the imaginary axis two times and the root µ 4 in (3.15) if there are four intersections. Assume there are four intersections with the imaginary axis. This leads to a contradiction by arguments analogous to the proof of the uniqueness of the strip maximal with respect to W Rˆtβu pt qzir in Proposition Thus, the local maximum is iµ 2 and similarly it follows that the local minimum is iµ 3. Assume that the local maximum is not on the imaginary axis. Then fps q s 2 pp 2 s Q s q from Lemma 3.11 and s satisfies one of (3.13), i) or ii). Since it is

19 NUMERICAL RANGE OF NONSELFADJOINT RATIONAL OPERATOR FUNCTIONS 19 a local maximum, fps ɛq ą and fps ` ɛq ă for sufficiently small ɛ ą, which implies s p d ` d a 1 ` 4β{p4c d 2 qq{4 or s p d ad 2 4βq{4. Assume that s p d ` d a 1 ` 4β{p4c d 2 qq{4, then d 2 ą 4β ` 4c since otherwise s is not negative. Then fps q : P s ap s 2 Q s s 2 ă and we have a contradiction to Proposition 3.7. Hence, s p d ad 2 4βq{4. The proof for the local minimum points is similar. i) Assume β ă 4c, then by Lemma 3.4 it follows that d 2 ą 2? β. Hence by Proposition 3.16 there is a gap if and only if d ą 2? β. The point ω is a local maximum not on the imaginary axis if and only if ω I p d ad 2 4βq{4 and P ωi ωi 2 ď. The condition P ω I ωi 2 β`4c ď holds if and only if d ě 2?. For the c local minimum the same idea is used. ii) Assume β ě 4c then by Lemma 3.4, d 2 ď 2? β. Hence Proposition 3.16 implies that there is a gap if and only if d ą d 2. The point ω is a local maximum not on the imaginary axis if and only if ω I p d ad 2 4βq{4 and P ωi ωi 2 ą, which never holds. For the local minimum the same idea is used for ω I p d ` ad 2 4βq{4 and then it follows that P ωi ωi 2 β`4c ą holds if and only if c ą and d ą 2? c. Proposition Let W Rˆtβu pt qzir be defined as in (3.4). If c set ˆd : and if c ą let ˆd be the unique solution of ˆβ2 (3.15) c ˆd 3 ` βc 3c2 ˆd 2 ` c 2 p2β ` 3cq 16 ˆd c 3 pβ ` cq, that satisfies ă ˆd ă c. Then if d ă 2a ˆd, the extreme points as in Definition 3.9 with smallest imaginary part of W Rˆtβu pt qzir are not on the imaginary axis and the imaginary part of the points are (3.16) 1 4 d d c 1 ` 4β 4c d 2 If d ě 2a ˆd the point with smallest imaginary part is iµ1 as defined in (3.14). Proof. From Proposition 3.5, Lemma 3.11, and Corollary 3.1 it follows that ω with imaginary part (3.16) and iµ 1 are the only possible points that can have smallest imaginary part. Moreover if the points with smallest imaginary part are not on the imaginary axis then d ă 2? c since otherwise ω I ă Impδ q does not hold. It thus follows that if c the point with smallest imaginary part is iµ. Assume c ą, then ω is a point with smallest imaginary part not on the imaginary axis if and only if d ă 2? c and by Proposition 3.7 follows P ωi ωi 2 ą. This holds if and only if d ă 2a ˆd, where ˆd is the unique solution to (3.15) satisfying ă ˆd ă c. In the remaining case the point with smallest imaginary part must be on the imaginary axis and thus iµ 1. In Figure 4.(b) the local minimum is clearly not on the imaginary axis but by increasing d, Figure 4.(c) is obtained, where the local maximum point is on the imaginary axis..

20 2 CHRISTIAN ENGSTRÖM AND AXEL TORSHAGE (a) (b) (c) Figure 4. Examples of the set W Rˆtβu pt q Ă C, where red and blue denotes W r,8sˆtβu pt q and W r 8,qˆtβu pt q, respectively. The figure describes how W Rˆtβu pt q changes with increasing d. In (a) there is no strip maximal with respect to W Rˆtβu pt qzir. In (b) such a strip exists and the point with smallest imaginary part is moving to the imaginary axis. In (c) the local maximum is moved to the imaginary axis. The operator function T defined in (1.1) depends on d and in the following we study the dependence of W Rˆtβu pt qzir on the parameter d. Moreover, we consider properties of a strip maximal with respect to W Rˆtβu pt qzir as in Definition 3.8. Proposition Let W Rˆtβu pt qzir be defined as in (3.4) and consider the extreme points of W Rˆtβu pt qzir depending on d. Then the extreme points of W Rˆtβu pt qzir in Definition 3.9 are continuous in d and the extreme points with smallest imaginary parts are decreasing with d. Furthermore, if a strip maximal with respect to W Rˆtβu pt qzir exists, then s is strictly decreasing and s 1 is strictly increasing with respect to d. Proof. Proposition 3.17, Proposition 3.18, and Proposition 3.7 yield the continuity of the extreme points. Hence, it is enough to show the results for the extreme points on and off the imaginary axis separately. For each extreme point on W Rˆtβu pt qzir, the result follows directly from Propositions 3.17 and All other extreme points can be written in the form ω iµ, µ P R where µ is a solution to (3.5) and the smallest imaginary part is then µ 1 as defined in (3.14). Let qβ d pµq denote the polynomial (3.5) for a given d ą and take ɛ ą. Then ˆ µ 2 ` dµ ` c ` β 4 (3.17) q d`ɛ β pµq q d βpµq ` 2µ ɛ ` ɛ 2 µ 2. Let µ d i for i 1, 2, 3, 4 denote the real roots of qd β ordered non-decreasingly. If the point with the smallest imaginary part is on the imaginary axis then it is given by µ d 1 and Proposition 3.5 implies pµ d 1q 2 ` dµ d 1 ` c ` β{4 ą. Thus q d`ɛ β pµ d 1q ă for ɛ ą small enough and (3.17) then implies that the smallest imaginary part of W Rˆtβu pt qpdqzir is decreasing in d. Assume that it exists a strip S that is maximal with respect to W Rˆtβu pt qzir and that the local maximum (minimum) point is (is 1 ) is on the imaginary axis. Then Proposition 3.17 implies s µ d 2 (s 1 µ d 3) and we conclude that q d β pµq ą for µ P S. In particular the maximum (minimum) is decreasing (increasing) if and

21 NUMERICAL RANGE OF NONSELFADJOINT RATIONAL OPERATOR FUNCTIONS 21 only if q d`ɛ β pµ d 2q ą pq d`ɛ β pµ d 3q ą q for ɛ ą small enough. In general, q d`ɛ β pµq ą for ɛ ą small enough if and only if µ 2 ` dµ ` c ` β{4 ď, which is equivalent to «c c ff d2 β d2 β µ P d 2 4 c, d 2 ` 4 c : I d. What remains to show is that all local extreme points µ d 2, µ d 3 are in I d. It can be seen that if for some d ą it holds that µ d 2, µ d 3 P I d, then µ d 2, µ d 3 P I d for all d ě d. If β ă 4c then it follows from Proposition 3.17 i) that only the local maximum can be on the imaginary axis. Moreover, the condition d ě d : pβ ` 4cq{p2? cq holds, which together with µ d 2?c P I d yields the result. If β ą 4c, there is a strip maximal with respect to W Rˆtβu pt qzir if and only if d ą d 2. Furthermore µ d2 2 µd2 3 holds. Thus it is enough to show µd2 2 P Id2 to prove the claim. Since β ą 4c, it follows from Lemma 3.4 that d 2 ď 2? β ă d. Hence, for d d 2 ě ɛ ą and d d 2 ` ɛ the set W Rˆtβu pt qzir has both the local minimum and maximum on the the imaginary axis. there is a strip with respect to W Rˆtβu pt qzir it follows from µ d2 2 µd2 3 However, under the assumption that that either µ d 2 is decreasing in the vicinity of d 2 or µ d 3 is in the vicinity of d 2 increasing in d. Hence either µ d 2 P I d or µ d 3 P I d. Since this holds for arbitrarily small ɛ ą, µ d2 2 µd2 3 P. Id2 For β 4c, the result follows immediately since the roots of (3.5) are continuous in β and Proposition Variation of the numerical range W pbq. In this section, we describe the subset of W BΩ pt q obtained when fixing α and varying β P W pbq. Let R` : r, 8s and consider the set 4ď (3.18) W tαuˆr`pt q r n pα, R`q, defined according to (2.6). Remark 3.2. In the definition of W tαuˆr`pt q, we set W pbq r, 8s since r, 8s is the smallest closed interval containing W pbq for all bounded B. The limit of the roots r n pα, βq are and 8 id{2 as β Ñ 8. These points are in W tαuˆr`pt q but for α not in W tαuˆw pbq pt q, for any bounded B. For α this is completely solved in Lemma 2.15, and we assume therefore that α P Rztu. The set W tαuˆr`pt q can in the variable β P R` be parametrized into a union of four curves. For ω P CzpiR Y tδ`, δ uq, ω P W Rˆtβu pt q is equivalent to ω P Π β as defined in (2.17), and α ˆαpωq in equation (2.14). The following results for this curve are similar to the results for W Rˆtβu pt q, but the behavior will greatly depend on the sign of α, as can be seen in Proposition 2.7 and in Figure 5. n 1 Proposition Let W tαuˆr`pt q denote the set (3.18) and take α, c, ν P Rzt 2c{du. Then iν P W tαuˆr`pt qzir if and only if α ą ν, or ν is a solution

22 22 CHRISTIAN ENGSTRÖM AND AXEL TORSHAGE (a) (b) (c) Figure 5. Visualization of how the set W tαuˆr`pt q depends on the sign of α, where α ą in panel (a), α in (b), and α ă in (c). to (3.19) q α pνq : ν 4 ` d 2 ν3 αd ν αc, ν ď 2c 2 d, where the inequality is strict if α ă c. For c, the statements of the proposition hold but zero is in that case not in the set W tαuˆr`pt qzir. Proof. The proof follows the same steps as in Proposition 3.2, with the additional condition β ě, which by Corollary 2.7 iii) simplifies to ν ď 2c d or ν on the imaginary axis. Corollary Let p α,β and q α denote the polynomials (2.3) and (3.19), respectively. Then p α,β has a root iν, ν P Rztu of multiplicity n ą 1 for some β P r, 8q if and only if ν ď 2c{d and ν is a root of q α of multiplicity n 1. Proof. Similar to Corollary 3.3, with the additional condition β ě. Lemma Let q α be the polynomial in (3.19) and let qα be its discriminant. Let d denote the largest real d solving qα. If c and α ă, then d, otherwise d is the unique positive solution. The following properties hold for q α and for d : i) If α ą then d P r4 a maxpα, cq, 8q and the polynomial q α has four real roots if d ě d, and if d ă d it has two real roots. ii) If α ă the polynomial q α has two real roots if d ě d, and if d ă d it has no real roots. If c ą then d P p, 2? cs. Proof. Let ˆd : d 2 {4 and study fp ˆdq : qα as a polynomial in ˆd. By definition, the discriminant of q α is (3.2) fp ˆdq 4α 3 ˆd3 α 2 p27α 2 6αc ` 27c 2 q ˆd 2 ` 192α 3 c 2 ˆd 256α 3 c 3. For c the roots are, 27α{4 and the result follows. If c ą the discriminant of f is f α 9 c 3 pα cq 4 pα ` cq 2. i) Assume α ą then f ď with equality only if c α. In that case ˆd 4c is a triple root of f, else f ă and then there is one real root. Hence, there is in each case exactly one real solution to fp ˆdq. Denote this solution ˆd and define d 2a ˆd. Then d is the unique positive solution to qa. Further since fp4αq, fp4cq ď and fp ˆdq Ñ 8, ˆd Ñ 8 it follows that the unique positive solution ˆd P r4 maxpα, cq, 8q. If ˆd ă ˆd then fp ˆdq ă and thus there are two real roots of

23 NUMERICAL RANGE OF NONSELFADJOINT RATIONAL OPERATOR FUNCTIONS 23 q α. If ˆd ě ˆd then fp ˆdq ď and since q α pq ď there is always at least one root and thus it must be four. ii) Assume α ă, then there are always three real roots since f ě. From fp 8cq ď, fpq ą and fp ˆdq Ñ 8, ˆd Ñ 8, follows that exactly one root is positive. Since fpcq ď the unique solution ˆd is in p, cs. Assume d ă d, then by continuity there is no root since q α has no real root for d and qα ą for d ă d. Assume ˆd ě ˆd, then fp ˆdq ď and thus there are two roots of q α. Proposition Let q α be the polynomial (3.19) and let d denote the largest real d solving qα. Then the set W tαuˆr`pt qzir, defined as in (3.18), intersects the imaginary axis in the following points, counting multiplicity: i) If α ą and c ą, there is one intersection in. There are no further intersections if d ă 2? c, one more intersection if 2? c ď d ă d, and three more intersections if d ě d. If α ą and c there is no intersections if d ă d and two intersections if d ě d. ii) If c ď α ă, there is no intersection if d ă 2? c, and one intersection if d ě 2? c. iii) If α ă c, there is no intersection if d ă d, two intersections if d ď d ă 2? c, and one intersection if d ě 2? c. Proof. The intersections of W tαuˆr pt qzir coincide with the roots of the polynomial (3.19) in Proposition In each case Lemma 3.23 is used to obtain the number of real roots. i) Assume α ą and c ą, then by Proposition 3.21, there is a simple intersection in. For d ă 2? c Proposition 3.21 implies that there are no solutions to (3.19). For d 2? c, the point 2c{d d{2 is a simple root of q α and by Proposition 3.21, there is a simple intersection of the imaginary axis. For ν ă d{2, q α pνq ą and hence there are no more intersections of the imaginary axis. For d ą 2? c, q α p 2c{dq ă, thus there is either one or three solutions of (3.19). If 2? c ă d ă d, there is one intersection. Assume d ě d. Then q α pνq ă for ν P p 2c{d, q and q α is convex for ν ě. Hence, (3.19) has three solutions. Now assume c, then Proposition 3.21 shows that zero does not give an intersection and thus there is one less intersection of the imaginary axis in this case. ii), For ν ď 2c{d the derivative qα 1 pνq is non-positive with equality only if α c. Thus there is at most one solution to (3.19). If 2? c ą d then q α p 2c{dq ą and there is no solution. For 2? c ď d, q α p 2c{dq ď and it follows from Proposition 3.21 and α ě c that (3.19) has one solution. Assume 2? c d and α c, then 2c{d is a double root of (3.5) and by Corollary 3.22 a triple root of p pα,q. From the injectivity stated in Lemma 2.8 iii) it follows that there is only one intersection of the imaginary axis at this point. The result then follows from Proposition iii) If d ă d then there are no real solutions to (3.19). If d ď d ă 2? c, then q α pνq ą for ν ą d{2. For d{2 ě ν ą 2c{d we have ˆ 2c q α pνq ě d ` d 2 ˆ d3 8 αd 2 α ˆc d2 4 d2 4 ˆc d2 ą. 4 Thus, the two roots satisfy ν ď 2c{d. Assume d 2? c, then by Proposition 3.19 the value 2c{d is not an intersection of the imaginary axis. However, there is an intersection at the distinct root 3a αd{2, which shows that we have one root. For

24 24 CHRISTIAN ENGSTRÖM AND AXEL TORSHAGE (a) (b) (c) Figure 6. Examples of the set W tαuˆr`pt q Ă C with different numbers of intersections of the imaginary axis. In the panels (a) and pbq are α ą and in panel (c) is α ă. d ą 2? c ą it follows that q α p 2c{dq ă and it can be seen that (3.19) has one solution. The multiplicity of a real root ν ă 2c{d of q α in (3.19) will determine the number of segments of W tαuˆr`pt qzir intersecting iν. However, this will in general not hold if 2c{d is a root of q α. This case is addressed in Proposition Figure 6.(b), shows an example where 2c{d 2 is an intersection of the imaginary axis, while in Figure 6.(c), 2c{d 2 is not an intersection despite being a root of q α. Cases with different numbers of intersections are illustrated in Figure 5 and in Figure 6. Proposition Assume c ą and let W tαuˆr`pt qzir and Π β denote the sets defined in (3.18) and in (2.17), respectively. Then ω P W tαuˆr`pt qzir if and only if ω 8 or one of the following conditions hold: i) d ď 2? c and ω δ, where the inequality is strict if α ă c. ii) ω P Π β ztr id{2u and one of the below four identities hold: c b (3.21) ω R R ωi Rω 2 I S ωi ωi 2, where (3.22) R ωi αd 2pd`2ω I q, S ω I 2αcω I d`2ω I. If c then R W tαuˆr`pt qzir but all other statements hold as for the case c ą. Proof. Similar to Proposition 3.7 subject to β ě. Remark Since in the limit β Ñ 8 the roots r n pα, βq approach and 8 id{2, it is convenient to assume the imaginary part in infinity is d{2. Using this convention, for α ą the largest imaginary part of W tαuˆr`pt qzir is and the smallest imaginary part is d{2. For α ă the largest imaginary part of W tαuˆr`pt qzir is d{2. Lemma Let W tαuˆr`pt qzir denote the set (3.18) and let R ωi and S ωi denote the expressions in (3.22). A point ω P W tαuˆr`pt qzir, with ω I R t, d{2u

25 NUMERICAL RANGE OF NONSELFADJOINT RATIONAL OPERATOR FUNCTIONS 25 is an extreme point in the sense of Definition 3.9 if and only if it is a distinct root to gpω I q : pd ` 2ω I q 2 pr 2 ω I S ωi q. The roots of g are (3.23) ω I d da 1 α c. 4 A double root of g is only possible if α c ą d 2 {16, where ω ac d 2 {16 id{4. Assume ω I is a double root of f. Then, ω is a point where more than one curve component in W tαuˆr`pt qzir intersect. Proof. The proof is similar to that of Lemma Figure 6.(b) show an example where the set W tαuˆr`pt q contains points with more than one curve component intersecting it in W tαuˆr`pt qzir. Lemma Let W tαuˆr`pt qzir and let R ωi, S ωi denote the expressions in (3.18) and in (3.22), respectively. A point iν P ir with ν P Rzt d{2, u is an extreme point to W tαuˆr`pt qzir in the sense of Definition 3.9 if and only if R ν ` ar 2 ν S ν ν 2 and ν is a distinct intersection of the imaginary axis. Proof. The condition R ν ` ar ν 2 S ν ν 2 implies R ν ą. Hence ν ą d{2 for α ą and ν ă d{2 for α ă. Then it follows from Proposition 2.7 iv) that the points (3.23) are obtained for positive β. The rest of the proof is similar to the proof of Lemma Lemma Let W tαuˆr`pt qzir be defined as in (3.18) and let γ P W tαuˆr`pt qzir be a bounded component such that Czγ consists of more than one component. Then for each β P r, 8q, either one root of the polynomial p α,β in (2.3) belongs to γ or one root of p α,β can be written as iµ, for µ P J : rminpγzirq I, maxpγzirq I s. Proof. The proof is similar to the proof of Lemma 3.29 with the exception that for a bounded component γ the set Czγ does not necessarily consists of more than one component. Proposition 3.3. Let W tαuˆr`pt qzir be defined as in (3.18) and let d be the constant in Lemma Then, if α ą and either α ă c or d ą d there is a unique strip with respect to W tαuˆr`pt qzir. In all other cases, there is no such strip. Proof. Similar to the proof of Proposition 3.16, using Proposition 3.25 and Lemmata 3.27, 3.28 and Definition Denote the k real solutions to (3.19) by (3.24) ν 1 ď ν 2 ď... ď ν k, where k might be zero. Proposition Let W tαuˆr`pt qzir be defined as in (3.18) and ν i as in Definition The extreme points of W tαuˆr`pt qzir in Definition 3.9 have the following properties: i) If α ě c there is a unique maximal strip with respect to W tαuˆr`pt qzir if and only if d ą d. The local maximum point is iν 1 and the local minimum is iν 2.

26 26 CHRISTIAN ENGSTRÖM AND AXEL TORSHAGE ii) If c ą α ą there is a unique maximal strip with respect to W tαuˆr`pt qzir. If (3.25) a) holds then the local maximum points are not on the imaginary axis and s p d d a 1 α{cq{4. If (3.25) a) does not hold, then the local maximum point is iν 1. If (3.25) b) holds then the local minimum points are not on the imaginary axis and s 1 p d ` d a 1 α{cq{4. If (3.25) b) does not hold, then the local minimum point is iν 2. iii) If α ă, the points in ω P W tαuˆr`pt qzir with smallest imaginary parts are not on the imaginary axis if (3.25) a) holds and ω I p d d a 1 ` α{cq{4. If (3.25) a) does not hold, then the point with the smallest imaginary part is iν 1. If c the point with smallest imaginary part is always iν 1. 4? c (3.25) a) d ă b 1 ` a1, b) d ă α c Proof. The proof is similar to Proposition ? c b 1 a1 α c Proposition Let W tαuˆr`pt qzir be defined as in (3.18) and let S denote a strip defined in Definition 3.8. The extreme points of W tαuˆr`pt qzir in Definition 3.9 are continuous in d and have the following properties: i) If α ą c and a maximal strip with respect to W tαuˆr`pt qzir exists, then the local minimum of W tαuˆr`pt qzir is increasing in d, and the local maximum is decreasing in d. ii) If c α ą and a maximal strip with respect to W tαuˆr`pt qzir exists, then the local minimum of W tαuˆr`pt qzir is i? c, and the local maximum is id{4 i a d 2 {16 c. iii) If c ą α ą then the local maximum and local minimum of W tαuˆr`pt qzir are decreasing in d. iv) For α ă the smallest imaginary part of W tαuˆr`pt qzir is decreasing in d. Proof. The continuity in d of the extreme points follows from Propositions 3.32 and i) iii) Assume that a maximal strip with respect to W tαuˆr`pt qzir exists. If an extreme point is located not on the imaginary axis then α ă c by Lemma 3.27 and then the result follows from (3.23). Thus it is sufficient to show the result for the roots of (3.19) and for the case α c it is easy to verify that ν 1 d{4 ad 2 {16 c and ν 2?c, which proves ii). Hence, assume that α c. Similarly as in Proposition 3.19, observe the polynomial (3.19), where d is shifted by some small ɛ ą is (3.26) qα d`ɛ pνq : qαpνq d ` ν 2 pν2 αqɛ. For α ą it follows from Proposition 3.3 that d ě d since otherwise the local minimum and maximum are not on the imaginary axis, and thus there are three roots. The local maximum is iν 1 and the local minimum is iν 2 by Proposition The local maximum (minimum) is decreasing (increasing) if and only if qα d`ɛ pν1 d q ă pqα d`ɛ pν2 d q ă q, for ɛ ą sufficiently small. In general, if qαpνq d, then

27 NUMERICAL RANGE OF NONSELFADJOINT RATIONAL OPERATOR FUNCTIONS 27 qα d`ɛ pνq ą for ɛ ą small enough if and only if pν 2 αq ą, and thus ν ă?α. Hence, what needs to be shown is whether ν1 d, ν2 d ă?α. Since q α p?αq αpα cq, there is an even number of roots in p 8,?αs if α ą c and an odd number of roots p 8,?αs if α ă c. The polynomial qαpνq d is concave only on the interval r d{4, s, hence there can at most be two roots in this interval. From Lemma 3.23 and that d ě d it follows that d{4 ď?α and thus there has to be at least one root in the interval p 8, d{4s. Hence, if α ą c there are two roots in p 8,?αs, thus ν1 d, ν2 d ă?α, which proves i). If α ă c it follows that qαp?αq 1 and qαp d{4q 1 are positive and since q α is concave on r d{4,?αs there is at most one root in that interval. On p 8, d{4s, q α is convex and since q α p d{4q ă it follows that q α has exactly one root in this interval. Since the number of roots is odd in p 8,?αs it follows that there is only one root in the interval, this proves iii). iv) The result follows from Lemma 3.27 if the point with smallest imaginary part is not on the imaginary axis. Assume α ă and that the point with smallest imaginary part is on the imaginary axis. Then the point with smallest imaginary part will be iν1 d, which implies ppν1 d q 2 αq ą and thus the smallest imaginary part is decreasing, due to (3.26). 4. Resolvent estimates In this section, the ɛ-pseudonumerical range is introduced and we determine an enclosure of this set. We show how the boundary of the new enclosure of the pseudospectra can be determined and give an estimate of the resolvent of (1.1). For given ɛ ą the ɛ-pseudospectrum σ ɛ pt q is the union of σpt q and the set of all ω P C such that }T 1 pωq} ą ɛ 1. An equivalent condition for ω P σ ɛ pt q is that there exists a function u P dom T, }u} 1 for which }T pωqu} ă ɛ. Such u is called an approximate eigenvector or ɛ-pseudomode [Dav7, p. 255]. To be able to see how well-behaved T 1 is close to W Ω pt q, as defined in (2.8), we will for ω P C Cztδ`, δ u make an upper estimate of the resolvent for the rational function (1.1). Definition 4.1. For an operator function T define the ɛ-pseudonumerical range as the set W ɛ pt q : W pt q Y tω P CzW pt q : Du P dom T, pt pωqu, uq {}u} 2 ă ɛu. From definition 4.1 it is clear that W ɛ pt q Ą W pt q. Moreover, the inequality }T pωqu}}u} ě pt pωqu, uq for u P dom T yields that W ɛ pt qzσpt q Ą σ ɛ pt qzσpt q. In the following, T is the rational operator function defined in (1.1). Similar to the enclosure of W pt q in (2.8), we define an enclosure of W ɛ pt q Ă C as (4.1) W ɛ ΩpT q : W Ω pt q Y tω P CzW Ω pt q : Dpα, βq P Ω, t pα,βq pωq ă ɛu, where t pα,βq pωq is given in (2.2). Remark 4.2. For ω P tδ`, δ, 8u, ω P W Ω pt q if and only if ω P WΩ ɛ pt q for all ɛ ą. Figure 7 illustrates W Ω pt q and WΩ ɛ pt q in two cases. Note in particular that the distance between a point ω P BWΩ ɛ pt q and BW ΩpT q is not constant. Define the set, (4.2) γ ɛ pt q : BΓ ɛ pt qzw Ω pt q, Γ ɛ pt q : tω P CzW BΩ pt q : Dpα, βq P BΩ, t α,β pωq ă ɛu.

28 28 CHRISTIAN ENGSTRÖM AND AXEL TORSHAGE (a) (b) Figure 7. Examples of the set W Ω pt q in red and blue and the set WΩ 1 pt q in grey with c 6, d 4. In (a) W paq r1, 8s, W pbq r, 11s and in (b) W paq r 8, 8s, W pbq r4, 11s. Theorem 4.3. Assume that ɛ ą and define BW ɛ Ω pt q and γɛ pt q as in (4.1) and in (4.2), respectively. Then BW ɛ Ω pt q γɛ pt q. Proof. Assume that ω R tδ`, δ, 8u. If ω P BW ɛ Ω pt q, then min pα,βqpω t pα,βq pωq ɛ by continuity. Conversely assume min pα,βqpω t pα,βq pωq ɛ and let pα, β q be such that t pα,β qpωq ɛ. Since ɛ ą and t pα,β q is non-constant and holomorphic, the minimum modulus principle states that for each neighborhood N of ω, there is an ω 1 P N such that t pα,β qpω 1 q ă ɛ. Hence, min pα,βqpω t pα,βq pω 1 q ă ɛ and thus ω P BW ɛ Ω pt q is equivalent to min pα,βqpω t pα,βq pωq ɛ. Assume ω R W Ω pt q, then min pα,βqpω t pα,βq pωq ą by definition (4.1). Since t pα,βq pωq is linear in α and β, and nonzero it follows that t pα,βq pωq attains its minimum for pα, βq P BΩ. Hence, min t pα,βqpωq min t pα,βqpωq ɛ, pα,βqpbω pα,βqpω and thus ω P BWΩ ɛ pt q. It then follows that BW Ω ɛ pt qztδ`, δ, 8u γ ɛ pt qztδ`, δ, 8u. By definition W pbq tu, which implies that the points tδ`, δ, 8u are not isolated. Then, since BWΩ ɛ pt q and γɛ pt q are closed sets, BWΩ ɛ pt q γɛ pt q follows by continuity. Remark 4.4. From the minimum modulus principle it follows that WΩ ɛ pt q has no new components compared with W Ω pt q. However, components disconnected in W Ω pt q might be connected in WΩ ɛ pt q; see Figure 7.(a). Define the complex constants (4.3) κ : ω 2 c idω ω 2, λ : ω2. From Definition 2.2 it follows that the absolute value of t pα,βq can be written in the form (4.4) t pα,βq pωq ˆ κi ˆα ˆλI. 1 κ R β λ R

29 NUMERICAL RANGE OF NONSELFADJOINT RATIONAL OPERATOR FUNCTIONS 29 The following proposition is proven by a technique similar to the active set algorithm for constrained linear least squares problems presented in [NW6, Chapter 16.5]. Proposition 4.5. Assume ω P C, and let WΩ ɛ pt q, κ, and λ be defined as in (4.1) and (4.3), respectively. Define the constant ɛ as follows: If κ R define c ɛ : min βκ I ` λ I 2 ` min α λ R 2. βpw pbq αpw paq If κ R, define the constants β inf and β sup as the values in W pbq closest to κiλ I κ R pinf W paq λ R q κ 2 and κiλ I κ R psup W paq λ R q κ 2, respectively. Define the interval inf W paq λr B :, sup W paq λ j R. κ R κ R If B X W pbq H set Ω 1 : tpinf W paq, β inf q, psup W paq, β sup qu. If B X W pbq H set Ω 1 : tpinf W paq, β inf q, psup W paq, β sup q, pα op, β op qu, where β op denotes the value in B X W pbq closest to λ I {κ I and α op : β op κ R ` λ R. If κ I the choice of β op P B X W pbq is arbitrary. Define ɛ : Then ω P W ɛ Ω pt q if and only if ɛ ą ɛ. a min βκi ` λ I 2 ` α κ R β λ R 2. pα,βqpω 1 Proof. Note that ω P WΩ ɛ pt q if and only if ɛ ą min pα,βqpω t pα,βq pωq. This is in the sense of (4.4) a constrained linear least squares optimization problem. If κ R the result is trivial. Otherwise, the minimizing value α of this linear problem is either on one of the endpoints of W paq or an α that makes the second equation solvable in β P W pbq. Computing the optimal β P W pbq in all of these cases gives three possible pairs in Ω that minimizes (4.4) and the result follows. Corollary 4.6. Let W ɛ Ω pt q be defined as in (4.1) and assume that ω P CzW ΩpT q. Then, }T 1 pωq} ď 1{ɛ, where ɛ is given by Proposition 4.5. Proof. Assume that ω P CzpW Ω pt qytδ`, δ uq, then from Proposition 4.5 it follows that min pα,βq t pα,βq pωq ɛ ą. Hence, the result follows from }T pωqu}{}u} ě pt pωqu, uq {}u} 2 ě min pα,βq t pα,βq pωq. Figure 8 illustrates the upper bound of }T 1 pωq} in Corollary 4.6 for arbitrary operators with W paq r 32, 4s, W pbq r, 4s, where a solid line depicts W BΩ pt qzir. For ω P W Ω pt q, we chose particular matrices A and B with the given numerical ranges and compute }T 1 pωq} numerically. Note that on the imaginary axis a simpler result than Proposition 4.5 holds. Proposition 4.7. Let W ɛ Ω pt q be defined as in (4.1). Let Aɛ be an arbitrary selfadjoint operator with W pa ɛ q rinf W paq ɛ, sup W paq ` ɛs and define T ɛ pωq : A ɛ ω 2 ω 2 c idω ω 2 B, dom T ɛ pωq dom A ɛ, ω P C. Let W Ω pt ɛ q be defined as in (2.8). Then W ɛ Ω pt q X ir W ΩpT ɛ q X ir.

30 3 CHRISTIAN ENGSTRÖM AND AXEL TORSHAGE Figure 8. Example of the upper bound of }T 1 pωq} from Corollary 4.6 for ω P CzW Ω pt q with W paq r 32, 4s, W pbq r, 4s, c 4, and d 4. For ω P W Ω pt q the resolvent is computed numerically for particular matrices A and B. Proof. Assume ω P irztδ`, δ u, then as in the proof of Theorem 4.3 we have ω P WΩ ɛ pt q if and only if min pα,βqpω t pα,βq pωq ď ɛ. Hence, ω P WΩ ɛ pt q if and only if (4.5) α ` ωi 2 ` c ` dω I ` ωi 2 β e, for some pα, βq P Ω and for some e with e ď ɛ. Assume ω P irztδ`, δ u satisfies (4.5), then α e P R, α e P W pa ɛ q and thus ω P W Ω pt ɛ q. The converse is obvious. Hence, WΩ ɛ pt q X irztδ`, δ u W Ω pt ɛ q X irztδ`, δ u, which by continuity yields WΩ ɛ pt q X ir W ΩpT ɛ q X ir. The set WΩ ɛ pt q X ir can then be obtained by the algorithm given in Proposition Acknowledgements. The authors gratefully acknowledge the support of the Swedish Research Council under Grant No ω 2 I References [AL95] [ALL1] V. M. Adamjan and H. Langer. Spectral properties of a class of rational operator valued functions. J. Operator Theory, 33(2): , V. Adamjan, H. Langer, and M. Langer. A spectral theory for a λ-rational Sturm- Liouville problem. J. Differential Equations, 171(2): , 21.