NIEP
🚧 Under Development

Site Notice: This research hub and its interactive spectral visualization tools are actively under development. Features, 3D dataset galleries, and solvers are being continuously expanded and refined.

Research & Visual Computing

The Nonnegative Inverse Eigenvalue Problem

An open-source, interactive collection of computational tools and geometric visualizations designed to explore, understand, and advance research on matrix spectra and nonnegativity.

Interactive Visualizations & Solvers

Select a tool to begin interactive spectral analysis

Karpelevič Region Viewer

Visualize the exact region $\mathcal{K}_n$ of all complex numbers that can arise as eigenvalues of an $n \times n$ stochastic matrix. Features live Farey sequence root calculation, Ito boundary polynomial solver, point containment testing, and dynamic power orbit visualization ($z^1, z^2, \dots, z^k$).

✓ Arbitrary order $n$ selection ($n = 2$ to $16$)
✓ Real-time containment test: $z \in \mathcal{K}_n$
✓ Power-up animation with trajectory trail
✓ High-DPI interactive canvas with pan & zoom

Trace Polytope Gallery

Explore interactive 2D polygons and 3D polyhedra capturing necessary and sufficient trace polytope conditions for the Symmetric (SNIEP) and Real (RNIEP) inverse eigenvalue problems for dimensions n = 4, 5, and 6. Filter by matrix class, dimension, and vertex signature.

✓ Filter by SNIEP, RNIEP, and matrix dimension
✓ Interactive Plotly 3D rotatable geometry
✓ Complete vertex listings & facet coordinates
✓ Searchable polytope index

Numerical Spectra Realizer

Given a target spectrum, determine whether it can be realized by a Nonnegative, Stochastic, Symmetric, or Stochastic Symmetric matrix. Configure numerical tolerances, run Orsi alternating projection solvers with real Schur factorizations, and inspect the realizing or best-approximation matrix.

✓ 4 target classes: Nonnegative, Stochastic, Symmetric, Doubly Stochastic
✓ Configurable tolerance (ε = 10−2 to 10−8) & best candidate fallback
✓ Exact Suleimanova, Soules, & circulant direct realizations
✓ Interactive heatmap, matrix export (LaTeX, NumPy, MATLAB), & plot

Academic Profile & Research

Curated academic website of Benjamin J. Clark, Ph.D. in Mathematics (WSU) and Software Engineer at Google. Browse research papers on matrix-preserving polynomials, conference presentations, student lecture notes, and curriculum vitae.

✓ 3 Journal Articles & 3 arXiv Preprints
✓ WSU Course teaching & downloadable lecture notes
✓ 8 Conference and seminar presentation slides
✓ Academic & software engineering CV
Theory & Research Foundations

About the Nonnegative Inverse Eigenvalue Problem

A comprehensive overview of spectral nonnegativity, major theoretical breakthroughs, and foundational literature.

📐 Core Formulation & Subproblems

The mathematical framework governing entrywise nonnegative realizations and their subclasses.

01

The NIEP Formulation

Given a multiset of $n$ complex numbers $\sigma = \{\lambda_1, \lambda_2, \dots, \lambda_n\}$, does there exist an entrywise nonnegative matrix $A \ge 0$ of order $n$ whose spectrum equals $\sigma$?

$$A = (a_{ij}) \in \mathbb{R}^{n \times n}, \quad a_{ij} \ge 0, \quad \det(\lambda I - A) = \prod_{i=1}^n (\lambda - \lambda_i)$$

Despite over a century of study since the roots of Perron-Frobenius theory, a general constructive characterization for arbitrary dimension $n \ge 5$ remains one of the premier unsolved problems in combinatorial matrix theory.

02

Subproblems & Variants

Several foundational variants isolate specific algebraic and topological structures:

  • SNIEP (Symmetric): Requires $A = A^T$, restricting all eigenvalues to $\mathbb{R}$.
  • RNIEP (Real): Requires $A \ge 0$ with real spectrum (not necessarily symmetric).
  • Stochastic NIEP: Requires row sums equal to $1$, constraining spectral radius $\rho(A) = 1$.
  • Nonzero Spectrum Problem: Realization by an $N \times N$ matrix with $N \ge n$, padding with $N - n$ zero eigenvalues.
03

Necessary Conditions

Any realizable spectrum $\sigma$ must satisfy universal invariants:

$$\begin{aligned} \mathbf{1.\; Reality:} & \quad \sigma = \bar{\sigma} \quad (\text{closed under complex conjugation}) \\[2pt] \mathbf{2.\; Perron\; Root:} & \quad \rho = \max_{1 \le i \le n} |\lambda_i| \in \sigma \\[2pt] \mathbf{3.\; Power\; Sum\; Traces:} & \quad s_k = \sum_{i=1}^n \lambda_i^k = \operatorname{tr}(A^k) \ge 0 \quad (\forall k \ge 1) \\[2pt] \mathbf{4.\; Loewy\text{-}London:} & \quad s_k^m \le n^{m-1} s_{km} \quad (\forall k, m \ge 1) \end{aligned}$$

While necessary, these trace conditions alone are not sufficient for dimensions $n \ge 5$.

🏆 Landmark Results & Milestones

Chronological developments establishing the boundaries, characterizations, and solvability of the NIEP.

Foundation 1907 · 1912

Perron-Frobenius Spectral Theory

O. Perron · G. Frobenius

Established that any nonnegative matrix $A \ge 0$ has a real, maximal eigenvalue $\rho(A) \ge |\lambda|$ with an associated nonnegative eigenvector. For irreducible matrices, $\rho(A)$ is algebraically simple, strictly positive, and any peripheral eigenvalues on the circle $|z| = \rho(A)$ are distributed with exact cyclic root-of-unity symmetry.

Sufficient Condition 1949

Suleĭmanova's Theorem

H. Suleĭmanova

Proved the earliest major sufficiency result for real spectra: if $\sigma$ contains a single positive eigenvalue $\lambda_0 > 0$ and all other $\lambda_i \le 0$, then $\sigma$ is realizable by a nonnegative matrix if and only if the trace condition $\sum_i \lambda_i \ge 0$ holds. This became the foundation for companion matrix constructions and perturbed realizations.

Exact Geometry 1951

Karpelevič's Boundary Characterization

F. I. Karpelevič

Completely determined the compact region $\mathcal{K}_n \subset \mathbb{C}$ of all possible eigenvalues of $n \times n$ stochastic matrices. The boundary consists of curvilinear arcs connecting roots of unity corresponding to Farey fractions $p/q \le 1$ with $q \le n$. In 1997, T. Ito reformulated these boundary curves as explicit minimal algebraic polynomial equations.

Polynomial Cone 1978 · 1979

Loewy-London Inequalities & Order 3 Resolution

R. Loewy · D. London

Introduced higher trace inequalities $s_k^m \le n^{m-1} s_{km}$, ruling out infinite candidate spectra that satisfy ordinary trace positivity. They posed the fundamental problem of characterizing the cone $\mathcal{P}_n$ of polynomials preserving nonnegative matrices, completely solving the NIEP for $n = 3$ and trace-zero $n = 4$.

General Solution 1981

Complete Solution for Dimension 3

C. R. Johnson

Provided a constructive proof completely solving the general NIEP for $n = 3$ with arbitrary trace. Proved that for $n = 3$, a spectrum is realizable if and only if it satisfies reality, the Perron-Frobenius condition, nonnegativity of trace, and the second-order Loewy-London trace inequality.

Symbolic Dynamics 1991

Boyle-Handelman Spectral Realization Theorem

M. Boyle · D. Handelman

Resolved the nonzero spectrum problem via shift equivalence in topological Markov chains. Proved that a candidate nonzero spectrum is the nonzero spectrum of some nonnegative matrix of possibly larger size $N \ge n$ if and only if it satisfies the Perron condition and $\operatorname{tr}(A^k) \ge 0$ with strict positivity for some power.

Dimension 4 1998 · 1999

Full Resolution of Dimension 4

H. Meehan · T. J. Laffey

Conclusively resolved the general NIEP for $n = 4$ including positive trace cases. Demonstrated that for order 4, the known Loewy-London inequalities, trace nonnegativity, and Perron conditions constitute necessary and sufficient conditions, constructing companion matrix factorizations for all boundary cases.

📚 Key Literature & Citations

Foundational and modern research papers shaping the spectral theory of nonnegative matrices.

Classical 1907

Zur Theorie der Matrices

Oskar Perron

Mathematische Annalen, Vol. 64, pp. 248–263

Classical 1912

Über Matrizen aus nicht negativen Elementen

Georg Frobenius

Sitzungsberichte der Königlich Preussischen Akademie der Wissenschaften, pp. 456–477

Sufficiency 1949

Stochastic matrices with real characteristic values

H. Suleĭmanova

Doklady Akademii Nauk SSSR, Vol. 66, pp. 343–345

Spectral Boundary 1951

On the characteristic roots of matrices with non-negative elements

F. I. Karpelevič

Izvestiya Rossiiskoi Akademii Nauk. Seriya Matematicheskaya, Vol. 15, No. 4, pp. 361–383

Trace Inequalities 1978

A note on the nonnegative inverse eigenvalue problem

Raphael Loewy · David London

Linear and Multilinear Algebra, Vol. 6, No. 1, pp. 83–90

Order 3 Solution 1981

Row stochastic matrices similar to doubly stochastic matrices

Charles R. Johnson

Linear and Multilinear Algebra, Vol. 10, No. 2, pp. 113–130

Shift Equivalence 1991

The spectra of nonnegative matrices via symbolic dynamics

Mike Boyle · David Handelman

Annals of Mathematics, Vol. 133, No. 2, pp. 249–316

Algebraic Boundary 1997

A new characterization of the Karpelevič region for stochastic matrices

Takehiro Ito

Linear Algebra and its Applications, Vol. 255, pp. 321–346

Order 4 Solution 1999

A result on the symmetric nonnegative inverse eigenvalue problem

Thomas J. Laffey · Helena Meehan

Linear Algebra and its Applications, Vol. 302–303, pp. 295–313

Realization Theory 2007

Nonnegative realization of spectra having negative real parts

Thomas J. Laffey · Helena Šmigoc

Linear Algebra and its Applications, Vol. 420, No. 2–3, pp. 672–689

Spectroids 2016

Perron spectroids and the boundary of the Karpelevič region

Charles R. Johnson · Pietro Paparella

Linear Algebra and its Applications, Vol. 507, pp. 215–229