Source author record

Mark Embree

Mark Embree appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA math.SP math-ph math.MP Numerical Analysis cond-mat.mtrl-sci cond-mat.supr-con math.CV math.DS math.OC quant-ph

Catalog footprint

What is connected

8works

11topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2026arXiv

A parametric Keldysh decomposition

Contour integral algorithms seek to compute a small number of eigenvalues located within a bounded region of the complex plane. These methods can be applied to both linear and nonlinear matrix eigenvalue problems. In the latter case, the foundation of these methods comes from the Keldysh decomposition, which breaks the nonlinear matrix-valued function into two parts: a rational function whose poles match the desired eigenvalues, and a remainder term that is analytic within the target region. Under contour integration this analytic part vanishes (via Cauchy's theorem), leaving only the component containing the desired eigenvalues. We propose an extension of the Keldysh decomposition for matrix-valued functions that depend analytically on an additional parameter. We establish key properties of this parametric Keldysh decomposition, and introduce an algorithm for solving parametric nonlinear eigenvalue problems that is based upon it.

preprint2022arXiv

How Descriptive are GMRES Convergence Bounds?

GMRES is a popular Krylov subspace method for solving linear systems of equations involving a general non-Hermitian coefficient matrix. The conventional bounds on GMRES convergence involve polynomial approximation problems in the complex plane. Three popular approaches pose this approximation problem on the spectrum, the field of values, or pseudospectra of the coefficient matrix. We analyze and compare these bounds, illustrating with six examples the success and failure of each. When the matrix departs from normality due only to a low-dimensional invariant subspace, we discuss how these bounds can be adapted to exploit this structure. Since the Arnoldi process that underpins GMRES provides approximations to the pseudospectra, one can estimate the GMRES convergence bounds as an iteration proceeds.

preprint2021arXiv

Mathematics of magic angles in a model of twisted bilayer graphene

We provide a mathematical account of the recent Physical Reviews Letter by Tarnopolsky--Kruchkov--Vishwanath. The new contributions are a spectral characterization of magic angles, its accurate numerical implementation and an exponential estimate on the squeezing of all bands as the angle decreases. Pseudospectral phenomena due to the non-hermitian nature of operators appearing in the model play a crucial role in our analysis.

preprint2020arXiv

Contour Integral Methods for Nonlinear Eigenvalue Problems: A Systems Theoretic Approach

Contour integral methods for nonlinear eigenvalue problems seek to compute a subset of the spectrum in a bounded region of the complex plane. We briefly survey this class of algorithms, establishing a relationship to system realization techniques in control theory. This connection motivates a new general framework for contour integral methods (for linear and nonlinear eigenvalue problems), building on recent developments in multi-point rational interpolation of dynamical systems. These new techniques, which replace the usual Hankel matrices with Loewner matrix pencils, incorporate general interpolation schemes and permit ready recovery of eigenvectors. Because the main computations (the solution of linear systems associated with contour integration) are identical for these Loewner methods and the traditional Hankel approach, a variety of new eigenvalue approximations can be explored with modest additional work. Numerical examples illustrate the potential of this approach. We also discuss how the concept of filter functions can be employed in this new framework, and show how contour methods enable a data-driven modal truncation method for model reduction.

preprint2015arXiv

Fast singular value decay for Lyapunov solutions with nonnormal coefficients

Lyapunov equations with low-rank right-hand sides often have solutions whose singular values decay rapidly, enabling iterative methods that produce low-rank approximate solutions. All previously known bounds on this decay involve quantities that depend quadratically on the departure of the coefficient matrix from normality: these bounds suggest that the larger the departure from normality, the slower the singular values will decay. We show this is only true up to a threshold, beyond which a larger departure from normality can actually correspond to faster decay of singular values: if the singular values decay slowly, the numerical range cannot extend far into the right-half plane.

preprint2014arXiv

Spectral Approximation for Quasiperiodic Jacobi Operators

Quasiperiodic Jacobi operators arise as mathematical models of quasicrystals and in more general studies of structures exhibiting aperiodic order. The spectra of these self-adjoint operators can be quite exotic, such as Cantor sets, and their fine properties yield insight into associated dynamical systems. Quasiperiodic operators can be approximated by periodic ones, the spectra of which can be computed via two finite dimensional eigenvalue problems. Since long periods are necessary to get detailed approximations, both computational efficiency and numerical accuracy become a concern. We describe a simple method for numerically computing the spectrum of a period-$K$ Jacobi operator in $O(K^2)$ operations, and use it to investigate the spectra of Schrödinger operators with Fibonacci, period doubling, and Thue-Morse potentials.

preprint2014arXiv

Spectral Properties of Schrödinger Operators Arising in the Study of Quasicrystals

We survey results that have been obtained for self-adjoint operators, and especially Schrödinger operators, associated with mathematical models of quasicrystals. After presenting general results that hold in arbitrary dimensions, we focus our attention on the one-dimensional case, and in particular on several key examples. The most prominent of these is the Fibonacci Hamiltonian, for which much is known by now and to which an entire section is devoted here. Other examples that are discussed in detail are given by the more general class of Schrödinger operators with Sturmian potentials. We put some emphasis on the methods that have been introduced quite recently in the study of these operators, many of them coming from hyperbolic dynamics. We conclude with a multitude of numerical calculations that illustrate the validity of the known rigorous results and suggest conjectures for further exploration.

preprint2007arXiv

The Fractal Dimension of the Spectrum of the Fibonacci Hamiltonian

We study the spectrum of the Fibonacci Hamiltonian and prove upper and lower bounds for its fractal dimension in the large coupling regime. These bounds show that as $λ\to \infty$, $\dim (σ(H_λ)) \cdot \log λ$ converges to an explicit constant ($\approx 0.88137$). We also discuss consequences of these results for the rate of propagation of a wavepacket that evolves according to Schrödinger dynamics generated by the Fibonacci Hamiltonian.