The Numerical Range of a Matrix

Total Page:16

File Type:pdf, Size:1020Kb

The Numerical Range of a Matrix The Numerical Range of a Matrix Kristin A. Camenga and Patrick X. Rault, joint work with Dan Rossi, Tsvetanka Sendova, and Ilya M. Spitkovsky Houghton College and SUNY Geneseo January 13, 2015 Rault & Camenga 1/13/15 January 13, 2015 1 / 23 Numerical Range Definition (Numerical Range) Let A be an n × n matrix with entries in C. Then the numerical range of A is given by n ∗ n W (A) = fhAx; xi : x 2 C ; kxk = 1g = fx Ax : x 2 C ; kxk = 1g: Rault & Camenga 1/13/15 January 13, 2015 2 / 23 Some Examples 2 1 0 1 3 2 −1 2 2 3 A = 4 0 i :5 5 B = 4 0 −1 2 5 0 0 −i 0 0 −1 Rault & Camenga 1/13/15 January 13, 2015 3 / 23 W (A) also contains the eigenvalues of A. Proof: If λ is an eigenvalue of A, then we pick a corresponding unit eigenvector, x. Then hAx; xi = hλx; xi = λhx; xi = λ. Basic Properties The numerical range of a matrix, W (A) is a compact and convex subset of C. Rault & Camenga 1/13/15 January 13, 2015 4 / 23 Basic Properties The numerical range of a matrix, W (A) is a compact and convex subset of C. W (A) also contains the eigenvalues of A. Proof: If λ is an eigenvalue of A, then we pick a corresponding unit eigenvector, x. Then hAx; xi = hλx; xi = λhx; xi = λ. Rault & Camenga 1/13/15 January 13, 2015 4 / 23 Unitary Matrices A matrix U is unitary if U∗U = I or U−1 = U∗. If U is a unitary n × n matrix, W (U∗AU) = W (A) for any n × n matrix A. We say W (A) is invariant under unitary similarities. If A is unitarily reducible, that is, unitarily similar to the direct sum A1 ⊕ A2, then W (A) is the convex hull of W (A1) [ W (A2). A1 0 A1 ⊕ A2 = 0 A2 Rault & Camenga 1/13/15 January 13, 2015 5 / 23 Unitarily Reducible Example 2 1 0 0 0 3 1 0 2 i 6 1 2 0 0 7 A1 = ; A2 = ; A1 ⊕ A2 = 6 7 1 2 1 3 4 0 0 2 i 5 0 0 1 3 Rault & Camenga 1/13/15 January 13, 2015 6 / 23 Normal Matrices A square matrix A is normal if A∗A = AA∗. A matrix is normal if and only if it is unitarily similar to a diagonal matrix. In this case, A is unitarily reducible to a direct sum of 1 × 1 matrices which are its eigenvalues, so W (A) is the convex hull of its eigenvalues. 2 2 0 0 0 3 6 0 1 − i 0 0 7 A = 6 7 4 0 0 1 + i 0 5 0 0 0 0 Rault & Camenga 1/13/15 January 13, 2015 7 / 23 Characterizing Shapes: 2 × 2 matrices If A is an irreducible 2 × 2 matrix, W (A) is an ellipse with foci at the eigenvalues. 0 3 A = −1 4 Rault & Camenga 1/13/15 January 13, 2015 8 / 23 Characterizing Shapes: 3 × 3 matrices Rault & Camenga 1/13/15 January 13, 2015 9 / 23 Characterizing Shapes: 4 × 4 Doubly Stochastic matrices A doubly stochastic matrix is one whose real number entries are non-negative and each row and column sums to 1. For example: 2 2 2 1 3 5 0 5 5 2 1 2 6 5 0 5 5 7 6 3 2 7 4 0 5 0 5 5 1 2 2 5 5 5 0 Key fact: Every doubly stochastic matrix is unitarily reducible to the direct sum of the matrix [1] with another matrix: ∗ U AU = [1] ⊕ A1; where U is a unitary matrix with real entries. Rault & Camenga 1/13/15 January 13, 2015 10 / 23 Characterizing Shapes: 4 × 4 Doubly Stochastic matrices Given the characterization of shapes of the numerical range for 3 × 3 T matrices and unitary reducibility U AU = [1] ⊕ A1;, there are three possibilities: W (A1) is the convex hull of a point and an ellipse (with the point lying either inside or outside the ellipse); the boundary of W (A1) contains a flat portion, with the rest of it lying on a 4th degree algebraic curve; W (A1) has an ovular shape, bounded by a 6th degree algebraic curve. All possibilities occur in these three categories and we can characterize which occurs. Rault & Camenga 1/13/15 January 13, 2015 11 / 23 4 × 4 D-S matrices with ellipse-based numerical range Theorem Let A be a 4 × 4 doubly stochastic matrix. Then the boundary of W (A) consists of elliptical arcs and line segments if and only if tr A3 − tr(AT A2) µ := tr A − 1 + tr(AT A) − tr A2 is an eigenvalue of A (multiple, if µ = 1). If, in addition, α = tr A−1−3µ > 0; β = (tr A−1−3µ)2−tr(AT A)+1+2(det A)/µ+µ2 > 0; then W (A) also has corner points at µ and 1, and thus four flat portions on the boundary. Otherwise, 1 is the only corner point of W (A); and @W (A) consists of two flat portions and one elliptical arc. Rault & Camenga 1/13/15 January 13, 2015 12 / 23 4 × 4 D-S examples: ellipses and line segments α = 7=24, β = −59=576 α = 43=96, β = −779=9216 α ≈ 0:65, β ≈ 0:31 α ≈ 0:77, β ≈ 0:21 Rault & Camenga 1/13/15 January 13, 2015 13 / 23 4 × 4 D-S Examples: 4th degree curves and flat portions 4th degree curves and 1 with flat 4th degree curve and 1 with flat portion on the left portion on the right Rault & Camenga 1/13/15 January 13, 2015 14 / 23 4 × 4 D-S Examples: Ovular Rault & Camenga 1/13/15 January 13, 2015 15 / 23 Graphing the boundary of the numerical range iφ Let Aφ = Ae . W (A) vs. W (Aφ) 0.4 0.2 0.2 0.1 f -0.4 -0.2 0.2 0.4 -0.4 -0.2 0.2 0.4 -0.1 -0.2 -0.2 -0.4 Rault & Camenga 1/13/15 January 13, 2015 16 / 23 Alternative (Kippenhahn, 1951): Let F (x : y : t) := det(H0x + K0y + It): Then W (A) is the convex hull of the dual curve to F (x : y : t) = 0. Graphing the boundary of the numerical range iφ Let Aφ = Ae . 0.4 W (A) vs. W (Aφ) ∗ ∗ A Aφ+Aφ Aφ−Aφ W f Hφ := 2 ; Kφ := 2i 0.2 W (Hφ) = Re(W (Aφ)) =[ λmin; λmax]. H L lmin W H f Hφv = λmaxv. f -0.3 -0.2 -0.1 0.1 0.2 0.3 λmax = hHφv; vi = RehAv; vi. lmax hAv; vi2 W (A). H L -0.2 Singularity: W (Aφ) has a vertical flat portion ) λmax has multiplicity 2 -0.4 Rault & Camenga 1/13/15 January 13, 2015 17 / 23 Graphing the boundary of the numerical range iφ Let Aφ = Ae . 0.4 W (A) vs. W (Aφ) ∗ ∗ A Aφ+Aφ Aφ−Aφ W f Hφ := 2 ; Kφ := 2i 0.2 W (Hφ) = Re(W (Aφ)) =[ λmin; λmax]. H L lmin W H f Hφv = λmaxv. f -0.3 -0.2 -0.1 0.1 0.2 0.3 λmax = hHφv; vi = RehAv; vi. lmax hAv; vi2 W (A). H L -0.2 Singularity: W (Aφ) has a vertical flat portion ) λmax has multiplicity 2 -0.4 Alternative (Kippenhahn, 1951): Let F (x : y : t) := det(H0x + K0y + It): Then W (A) is the convex hull of the dual curve to F (x : y : t) = 0. Rault & Camenga 1/13/15 January 13, 2015 17 / 23 The Gau-Wu number Definition (Gau-Wu number, 2013) k(A) = max #fx1;:::; xk g 8j;hAxj ;xj i2@W (A) fx1;:::;xk g orthonormal n fx1;:::;xk }⊂C Basic results: 1 ≤ k(A) ≤ n Points on parallel support lines of W (A) come from orthogonal vectors. Rault & Camenga 1/13/15 January 13, 2015 18 / 23 Vertical flat portion `1 ) pair of orthogonal eigenvectors u; v of H0, with hBu; ui; hBv; vi 2 `1 \ @W (B). Let `2 be a parallel support line, and let hAw; wi 2 `2 \ @W (B). Then w ? u; w ? v. Figure: B 2 M ( ) 3 C Thus k(B) = 3. Basic examples Basic results: 1 ≤ k(A) ≤ n k(A) ≥ 2 if n ≥ 2 Figure: A 2 M2(C). k(A) = 2 Rault & Camenga 1/13/15 January 13, 2015 19 / 23 Let `2 be a parallel support line, and let hAw; wi 2 `2 \ @W (B). Then w ? u; w ? v. Thus k(B) = 3. Basic examples Basic results: 1 ≤ k(A) ≤ n k(A) ≥ 2 if n ≥ 2 Figure: A 2 M2(C). k(A) = 2 Vertical flat portion `1 ) pair of orthogonal eigenvectors u; v of H0, with hBu; ui; hBv; vi 2 `1 \ @W (B). Figure: B 2 M3(C) Rault & Camenga 1/13/15 January 13, 2015 19 / 23 Thus k(B) = 3. Basic examples Basic results: 1 ≤ k(A) ≤ n k(A) ≥ 2 if n ≥ 2 Figure: A 2 M2(C). k(A) = 2 Vertical flat portion `1 ) pair of orthogonal eigenvectors u; v of H0, with hBu; ui; hBv; vi 2 `1 \ @W (B). Let `2 be a parallel support line, and let hAw; wi 2 `2 \ @W (B). Then w ? u; w ? v. Figure: B 2 M3(C) Rault & Camenga 1/13/15 January 13, 2015 19 / 23 Basic examples Basic results: 1 ≤ k(A) ≤ n k(A) ≥ 2 if n ≥ 2 Figure: A 2 M2(C).
Recommended publications
  • Some Properties of the Essential Numerical Range on Banach Spaces
    Pure Mathematical Sciences, Vol. 6, 2017, no. 1, 105 - 111 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/pms.2017.7912 Some Properties of the Essential Numerical Range on Banach Spaces L. N. Muhati, J. O. Bonyo and J. O. Agure Department of Pure and Applied Mathematics Maseno University, P.O. Box 333 - 40105, Maseno, Kenya Copyright c 2017 L. N. Muhati, J. O. Bonyo and J. O. Agure. This article is dis- tributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Abstract We study the properties of the essential algebraic numerical range as well as the essential spatial numerical range for Banach space operators. Mathematics Subject Classification: 47A10; 47A05; 47A53 Keywords: Essential algebraic numerical range, Essential spatial numer- ical range, Calkin algebra 1 Introduction Let A be a complex normed algebra with a unit and let A∗ denote its dual space. For a comprehensive theory on normed algebras, we refer to [4]. We define the algebraic numerical range of an element a 2 A by V (a; A) = ff(a): f 2 A∗; f(1) = 1 = kfkg. Let X denote a complex Banach space and L(X) be the Banach algebra of all bounded linear operators acting on X. We de- note the algebraic numerical range of any T 2 L(X) by V (T; L(X)). For T 2 L(X), the spatial numerical range is defined by W (T ) = fhT x; x∗i : x 2 X; x∗ 2 X∗; kxk = 1 = kx∗k = hx; x∗ig.
    [Show full text]
  • Numerical Range of Operators, Numerical Index of Banach Spaces, Lush Spaces, and Slicely Countably Determined Banach Spaces
    Numerical range of operators, numerical index of Banach spaces, lush spaces, and Slicely Countably Determined Banach spaces Miguel Martín University of Granada (Spain) June 5, 2009 Preface This text is not a book and it is not in its nal form. This is going to be used as classroom notes for the talks I am going to give in the Workshop on Geometry of Banach spaces and its Applications (sponsored by Indian National Board for Higher Mathematics) to be held at the Indian Statistical Institute, Bangalore (India), on 1st - 13th June 2009. Miguel Martín Contents Preface 3 Basic notation7 1 Numerical Range of operators. Surjective isometries9 1.1 Introduction.....................................9 1.2 The exponential function. Isometries....................... 15 1.3 Finite-dimensional spaces with innitely many isometries............ 17 1.3.1 The dimension of the Lie algebra..................... 19 1.4 Surjective isometries and duality......................... 20 2 Numerical index of Banach spaces 23 2.1 Introduction..................................... 23 2.2 Computing the numerical index.......................... 23 2.3 Numerical index and duality............................ 28 2.4 Banach spaces with numerical index one..................... 32 2.5 Renorming and numerical index.......................... 34 2.6 Finite-dimensional spaces with numerical index one: asymptotic behavior.. 38 2.7 Relationship to the Daugavet property....................... 39 2.8 Smoothness, convexity and numerical index 1 .................. 43 3 Lush spaces 47 3.1 Examples of lush spaces.............................. 49 5 3.2 Lush renormings.................................. 52 3.3 Some reformulations of lushness.......................... 53 3.4 Lushness is not equivalent to numerical index 1 ................. 55 3.5 Stability results for lushness............................ 56 4 Slicely countably determined Banach spaces 59 4.1 Slicely countably determined sets........................
    [Show full text]
  • Numerical Range and Functional Calculus in Hilbert Space
    Numerical range and functional calculus in Hilbert space Michel Crouzeix Abstract We prove an inequality related to polynomial functions of a square matrix, involving the numerical range of the matrix. We also show extensions valid for bounded and also unbounded operators in Hilbert spaces, which allow the development of a functional calculus. Mathematical subject classification : 47A12 ; 47A25 Keywords : functional calculus, numerical range, spectral sets. 1 Introduction In this paper H denotes a complex Hilbert space equipped with the inner product . , . and h i corresponding norm v = v, v 1/2. We denote its unit sphere by Σ := v H; v = 1 . For k k h i H { ∈ k k } a bounded linear operator A (H) we use the operator norm A and denote by W (A) its ∈ L k k numerical range : A := sup A v , W (A) := Av,v ; v ΣH . k k v ΣH k k {h i ∈ } ∈ Recall that the numerical range is a convex subset of C and that its closure contains the spectrum of A. We keep the same notation in the particular case where H = Cd and A is a d d matrix. × The main aim of this article is to prove that the inequality p(A) 56 sup p(z) (1) k k≤ z W (A) | | ∈ holds for any matrix A Cd,d and any polynomial p : C C (independently of d). ∈ → It is remarkable that the completely bounded version of this inequality is also valid. More precisely, we consider now matrix-valued polynomials P : C Cm,n, i.e., for z C, P (z) = → ∈ (p (z)) is a matrix, with each entry p (.) a polynomial C C.
    [Show full text]
  • Numerical Range of Square Matrices a Study in Spectral Theory
    Numerical Range of Square Matrices A Study in Spectral Theory Department of Mathematics, Linköping University Erik Jonsson LiTH-MAT-EX–2019/03–SE Credits: 16 hp Level: G2 Supervisor: Göran Bergqvist, Department of Mathematics, Linköping University Examiner: Milagros Izquierdo, Department of Mathematics, Linköping University Linköping: June 2019 Abstract In this thesis, we discuss important results for the numerical range of general square matrices. Especially, we examine analytically the numerical range of complex-valued 2 × 2 matrices. Also, we investigate and discuss the Gershgorin region of general square matrices. Lastly, we examine numerically the numerical range and Gershgorin regions for different types of square matrices, both contain the spectrum of the matrix, and compare these regions, using the calculation software Maple. Keywords: numerical range, square matrix, spectrum, Gershgorin regions URL for electronic version: urn.kb.se/resolve?urn=urn:nbn:se:liu:diva-157661 Erik Jonsson, 2019. iii Sammanfattning I denna uppsats diskuterar vi viktiga resultat för numerical range av gene- rella kvadratiska matriser. Speciellt undersöker vi analytiskt numerical range av komplexvärda 2 × 2 matriser. Vi utreder och diskuterar också Gershgorin områden för generella kvadratiska matriser. Slutligen undersöker vi numeriskt numerical range och Gershgorin områden för olika typer av matriser, där båda innehåller matrisens spektrum, och jämför dessa områden genom att använda beräkningsprogrammet Maple. Nyckelord: numerical range, kvadratisk matris, spektrum, Gershgorin områden URL för elektronisk version: urn.kb.se/resolve?urn=urn:nbn:se:liu:diva-157661 Erik Jonsson, 2019. v Acknowledgements I’d like to thank my supervisor Göran Bergqvist, who has shown his great interest for linear algebra and inspired me for writing this thesis.
    [Show full text]
  • Morphophoric Povms, Generalised Qplexes, Ewline and 2-Designs
    Morphophoric POVMs, generalised qplexes, and 2-designs Wojciech S lomczynski´ and Anna Szymusiak Institute of Mathematics, Jagiellonian University,Lojasiewicza 6, 30-348 Krak´ow, Poland We study the class of quantum measurements with the property that the image of the set of quantum states under the measurement map transforming states into prob- ability distributions is similar to this set and call such measurements morphophoric. This leads to the generalisation of the notion of a qplex, where SIC-POVMs are re- placed by the elements of the much larger class of morphophoric POVMs, containing in particular 2-design (rank-1 and equal-trace) POVMs. The intrinsic geometry of a generalised qplex is the same as that of the set of quantum states, so we explore its external geometry, investigating, inter alia, the algebraic and geometric form of the inner (basis) and the outer (primal) polytopes between which the generalised qplex is sandwiched. In particular, we examine generalised qplexes generated by MUB-like 2-design POVMs utilising their graph-theoretical properties. Moreover, we show how to extend the primal equation of QBism designed for SIC-POVMs to the morphophoric case. 1 Introduction and preliminaries Over the last ten years in a series of papers [2,3, 30{32, 49, 68] Fuchs, Schack, Appleby, Stacey, Zhu and others have introduced first the idea of QBism (formerly quantum Bayesianism) and then its probabilistic embodiment - the (Hilbert) qplex, both based on the notion of symmetric informationally complete positive operator-valued measure
    [Show full text]
  • Numerical Ranges and Applications in Quantum Information
    Numerical Ranges and Applications in Quantum Information by Ningping Cao A Thesis presented to The University of Guelph In partial fulfilment of requirements for the degree of Doctor of Philosophy in Mathematics & Statistics Guelph, Ontario, Canada © Ningping Cao, August, 2021 ABSTRACT NUMERICAL RANGES AND APPLICATIONS IN QUANTUM INFORMATION Ningping Cao Advisor: University of Guelph, 2021 Dr. Bei Zeng Dr. David Kribs The numerical range (NR) of a matrix is a concept that first arose in the early 1900’s as part of efforts to build a rigorous mathematical framework for quantum mechanics andother challenges at the time. Quantum information science (QIS) is a new field having risen to prominence over the past two decades, combining quantum physics and information science. In this thesis, we connect NR and its extensions with QIS in several topics, and apply the knowledge to solve related QIS problems. We utilize the insight offered by NRs and apply them to quantum state inference and Hamiltonian reconstruction. We propose a new general deep learning method for quantum state inference in chapter 3, and employ it to two scenarios – maximum entropy estimation of unknown states and ground state reconstruction. The method manifests high fidelity, exceptional efficiency, and noise tolerance on both numerical and experimental data. A new optimization algorithm is presented in chapter 4 for recovering Hamiltonians. It takes in partial measurements from any eigenstates of an unknown Hamiltonian H, then provides an estimation of H. Our algorithm almost perfectly deduces generic and local generic Hamiltonians. Inspired by hybrid classical-quantum error correction (Hybrid QEC), the higher rank (k : p)-matricial range is defined and studied in chapter 5.
    [Show full text]
  • On Invariance of the Numerical Range and Some Classes of Operators In
    ON INVARIANCE OF THE NUMERICAL RANGE AND SOME CLASSES OF OPERATORS IN HILBERT SPACES By Arthur Wanyonyi Wafula A thesis submitted in fulfillment of the requirements for the award of the Degree of Doctor of Philosophy in Pure Mathematics School of Mathematics,University of Nairobi 15 April 2013 G^ o a op h University of NAIROBI Library 0430140 4 0.1 DECLARATION This thesis is my orginal work and has not been presented for a degree in any other University S ignat ure.....D ate . .....A* / ° 5 / 7'tTL. m h I Arthur Wanyonyi Wanambisi Wafula This thesis has been submitted for examination with our approval as University super­ visors; Signature Date 2- Signature P rof. J. M. Khalagai Prof. G.P.Pokhariyal 1 0.2 DEDICATION This work is dedicated to my family ; namely, my wife Anne Atieno and my sons Davies Magudha and Mishael Makali. n 0.3 ACKNOWLEDGEMENT My humble acknowledgement go to the Almighty God Jehovah who endows us withdhe gift of the Human brain that enables us to be creative and imaginative. I would also like to express my sincere gratitude to my Supervisor Prof.J.M.Khalagai for tirelessly guiding me in my research.I also wish to thank Prof.G.P.Pokhariyal for never ceasing to put me on my toes. I sincerely wish to thank the Director of the School of Mathematics DrJ .Were for his invaluable support. Furthermore,! am indebted to my colleagues Mrs.E.Muriuki, Dr.Mile, Dr.I.Kipchichir, Mr R.Ogik and Dr.N.'Owuor who have been a source of inspiration.
    [Show full text]
  • Uncertainty Relations and Joint Numerical Ranges
    Konrad Szyma´nski Uncertainty relations and joint numerical ranges MSc thesis under the supervision of Prof. Karol Zyczkowski˙ arXiv:1707.03464v1 [quant-ph] 11 Jul 2017 MARIANSMOLUCHOWSKIINSTITUTEOFPHYSICS, UNIWERSYTETJAGIELLO NSKI´ KRAKÓW,JULY 2 0 1 7 Contents Quantum states 1 Joint Numerical Range 5 Phase transitions 11 Uncertainty relations 19 Conclusions and open questions 23 Bibliography 25 iii To MS, AS, unknow, DS, JJK, AM, FF, KZ,˙ SM, P, N. Introduction Noncommutativity lies at the heart of quantum theory and provides a rich set of mathematical and physical questions. In this work, I address this topic through the concept of Joint Numerical Range (JNR) — the set of simultaneously attainable expectation values of multiple quantum observables, which in general do not commute. The thesis is divided into several chapters: Quantum states Geometry of the set of quantum states of size d is closely related to JNR, which provides a nice tool to analyze the former set. The problem of determining the intricate structure of this set is known to quickly become hard as the dimensionality grows (approximately as d2). The difficulty stems from the nonlinear constraints put on the set of parameters. In this chapter I briefly introduce the formalism of density matrices, which will prove useful in the later sections. Joint Numerical Range Joint Numerical Range of a collection of k operators, or JNR for short, is an object capturing the notion of simultaneous measurement of averages — expectation values of multiple observables. This is precisely the set of values which are simultaneously attainable for fixed observables (F1, F2,..., Fk) over a given quantum state r 2 Md: if we had a function taking quantum state and returning the tuple of average values: E(r) = (hF1ir , hF1ir ,..., hFkir),(1) Joint Numerical Range L(F1,..., Fk) is precisely E[Md].
    [Show full text]
  • Spectral Inclusion for Unbounded Block Operator Matrices ✩
    Journal of Functional Analysis 256 (2009) 3806–3829 www.elsevier.com/locate/jfa Spectral inclusion for unbounded block operator matrices ✩ Christiane Tretter Mathematisches Institut, Universität Bern, Sidlerstr. 5, 3012 Bern, Switzerland Received 15 September 2008; accepted 23 December 2008 Available online 23 January 2009 Communicated by N. Kalton Abstract In this paper we establish a new analytic enclosure for the spectrum of unbounded linear operators A admitting a block operator matrix representation. For diagonally dominant and off-diagonally dominant block operator matrices, we show that the recently introduced quadratic numerical range W 2(A) contains the eigenvalues of A and that the approximate point spectrum of A is contained in the closure of W 2(A). This provides a new method to enclose the spectrum of unbounded block operator matrices by means of the non-convex set W 2(A). Several examples illustrate that this spectral inclusion may be considerably tighter than the one by the usual numerical range or by perturbation theorems, both in the non-self-adjoint case and in the self-adjoint case. Applications to Dirac operators and to two-channel Hamiltonians are given. © 2009 Elsevier Inc. All rights reserved. Keywords: Spectrum; Unbounded linear operator; Block operator matrix; Numerical range; Quadratic numerical range 1. Introduction The spectra of linear operators play a crucial role in many branches of mathematics and in numerous applications. Analytic information on the spectrum is, in general, hard to obtain and numerical approximations may not be reliable, in particular, if the operator is not self-adjoint or normal. In this paper, we suggest a new analytic enclosure for the spectrum of unbounded linear operators that admit a so-called block operator matrix representation.
    [Show full text]
  • Numerical Range: (In) a Matrix Nutshell
    Numerical range: (in) a matrix nutshell Panayiotis J. Psarrakos∗ Michael J. Tsatsomeros† August 12, 2002 ∗Department of Mathematics, National Technical University, Zografou Campus, Athens 15780, Greece. †Mathematics Department, Washington State University, Pullman, Washington 99164-3113. 1 The definition and a name The numerical range, field of values, Wertovorrat, Hausdorff domain, range of values have been, among others, competing names for the same notion. The first two have outlived the competition but the advantage goes decisively to the first: The numerical range of an operator T in a complex Hilbert space H. This is the term coined by Marshall Stone in 1932 [29] for the subset of the complex plane C given by W (T ) = Tx,x : x H, x = 1 . {h i ∈ k k } In technical terms here, ‘T is a bounded linear operator on a separable Hilbert space H with inner product , and induced norm x = x,x ’. But to keep our lives simple1, let’s confine our h · · i k k h i discussion to finite dimensions and talk aboutp the numerical range of an n n complex matrix A ∗ ×n ∗ (we write A n(C)), where x = √x x is the Euclidean length of x C and x its conjugate transpose; that∈ M is, k k ∈ W (A) = x∗Ax : x Cn, x∗x = 1 . { ∈ } This is a remarkable set for many reasons. It succinctly captures information about A as a transfor- mation, particularly about the eigenvalues and, perhaps more importantly, about the eigenspaces of A. The shape of W (A) (more accurately its boundary) is also related to the combinatorial structure of A, by which we mean the pattern of signs of its entries or the pattern of its nonzero entries, some- times represented by a graph.
    [Show full text]
  • Randomized Methods for Computing Low-Rank Approximations of Matrices
    Randomized methods for computing low-rank approximations of matrices by Nathan P. Halko B.S., University of California, Santa Barbara, 2007 M.S., University of Colorado, Boulder, 2010 A thesis submitted to the Faculty of the Graduate School of the University of Colorado in partial fulfillment of the requirements for the degree of Doctor of Philosophy Department of Applied Mathematics 2012 This thesis entitled: Randomized methods for computing low-rank approximations of matrices written by Nathan P. Halko has been approved for the Department of Applied Mathematics Per-Gunnar Martinsson Keith Julian David M. Bortz Francois G. Meyer Date The final copy of this thesis has been examined by the signatories, and we find that both the content and the form meet acceptable presentation standards of scholarly work in the above mentioned discipline. iii Halko, Nathan P. (Ph. D., Applied Mathematics) Randomized methods for computing low-rank approximations of matrices Thesis directed by Professor Per-Gunnar Martinsson Randomized sampling techniques have recently proved capable of efficiently solving many standard problems in linear algebra, and enabling computations at scales far larger than what was previously possible. The new algorithms are designed from the bottom up to perform well in modern computing environments where the expense of communication is the primary constraint. In extreme cases, the algorithms can even be made to work in a streaming environment where the matrix is not stored at all, and each element can be seen only once. The dissertation describes a set of randomized techniques for rapidly constructing a low-rank ap- proximation to a matrix.
    [Show full text]
  • On Characterization of Posinormal Operators and Their Applications in Quantum Mechanics
    Vol-1, Issue-2 PP. 212-221 ISSN: 2394-5788 ON CHARACTERIZATION OF POSINORMAL OPERATORS AND THEIR APPLICATIONS IN QUANTUM MECHANICS JOHN P. AKELLO C. RWENYO SABASI N. B. OKELO DEPARTMENT OF OMAORO and P. O. SCHOOL OF MATHEMATICS MATHEMATICS MOGOTU AND ACTUARIAL SCIENCE, KIBABII UNIVERSITY DEPARTMENT OF JARAMOGI OGINGA ODINGA COLLEGE MATHEMATICS UNIVERSITY OF SCIENCE KISII UNIVERSITY AND TECHNOLOGY, BONDO-KENYA ABSTRACT Studies on Hilbert space operators have been carried out by several researchers and mathematicians with interesting results obtained. A posinormal operator is one of the operators under concern with interesting characters. Properties of posinormal operators have not been exhausted hence this calls for intense characterization of these operators. Of great concern is the norm inequality involving these operators. In this paper, we give a detailed study of the norm property of posinormal operators when they are sefadjoint. The objectives of the study are to: establish norm inequalities for posinormal operators; characterize posinormal operators and; determine applications of posinormal operators in other fields. The methodology involved the use of tensor products as a technical approach in determining these norm inequalities. The results show that the norm of a posinormal operator is equal to the norm of any normal operator when the operators are selfadjoint. The results obtained are very important in the classification of Hilbert space operators and their applications in other fields like quantum theory. General Terms: Hilbert space, Posinormal operators, Selfadjoint, norm inequality, Applications and Quantum Mechanics. 1. INTRODUCTION The concept of arbitrary vector spaces can be generalized to inner product spaces and complete inner product spaces called Hilbert Spaces.
    [Show full text]