Theory Wiki · Subproblem Analysis

Symmetric Nonnegative Inverse Eigenvalue Problem

The study of real spectra realizable by symmetric entrywise nonnegative matrices ($A = A^T \ge 0$), orthogonal eigenspaces, and convex trace polyhedra.

1. Definition & Orthogonal Geometry

Enforcing matrix symmetry imposes severe geometric constraints on eigenvector frames.

The Symmetric Nonnegative Inverse Eigenvalue Problem (SNIEP) asks: given a multiset of $n$ real numbers $\sigma = \{\lambda_1, \lambda_2, \dots, \lambda_n\} \subset \mathbb{R}$, does there exist an entrywise nonnegative, symmetric matrix $A \in \mathbb{R}^{n \times n}$ whose spectrum is $\sigma$?

Definition: The SNIEP

Formal Statement

Find a matrix $A = (a_{ij}) \in \mathbb{R}^{n \times n}$ such that:

$$\begin{aligned} A &= A^T \quad (\text{symmetry}) \\ a_{ij} &\ge 0 \quad (\text{entrywise nonnegativity}) \\ \sigma(A) &= \{\lambda_1, \lambda_2, \dots, \lambda_n\} \end{aligned}$$

By the Spectral Theorem for real symmetric matrices, any realization $A$ is orthogonally diagonalizable:

$$A = Q \Lambda Q^T, \quad Q Q^T = I_n, \quad \Lambda = \operatorname{diag}(\lambda_1, \dots, \lambda_n)$$

This means the entrywise nonnegativity condition $a_{ij} \ge 0$ translates into a system of quadratic inequalities on the rows of an orthogonal matrix $Q = (q_{ik})$:

$$a_{ij} = \sum_{k=1}^n \lambda_k q_{ik} q_{jk} \ge 0 \quad (\forall 1 \le i, j \le n)$$

Because the rows of $Q$ must form an orthonormal basis of $\mathbb{R}^n$ ($\sum_k q_{ik}^2 = 1$ and $\sum_k q_{ik} q_{jk} = 0$ for $i \ne j$), the eigenvectors cannot tilt freely. This rigid orthogonality makes the SNIEP strictly harder to satisfy than general real realizations when dimension $n \ge 5$.

2. Necessary Invariants for SNIEP

Moments, trace Frobenius norms, and spectral partitioning.

1

Real Spectrum

$$\sigma \subset \mathbb{R}$$

All eigenvalues of a real symmetric matrix are guaranteed to be real.

2

Perron Root

$$\lambda_1 = \max_{i} |\lambda_i| \in \sigma$$

The maximal eigenvalue in absolute value is positive and belongs to $\sigma$.

3

Frobenius Trace Norm

$$s_{2k} = \|A^k\|_F^2 \ge 0$$

Even power sum traces equal the squared Frobenius norm of the symmetric power $A^k$.

4

Fiedler Partitioning

$$\lambda_1 + \lambda_n \ge 0$$

For irreducible $A$, the sum of extreme eigenvalues satisfies positivity bounds.

3. Sufficient Conditions & Realization Constructions

Techniques to constructively synthesize symmetric nonnegative realizations.

Fiedler's Suleĭmanova Theorem (1974)

Exact Sufficiency

A spectrum $\sigma = \{\lambda_1, \lambda_2, \dots, \lambda_n\}$ is called a Suleĭmanova spectrum if exactly one eigenvalue is positive and all others are non-positive:

$$\lambda_1 > 0 \ge \lambda_2 \ge \lambda_3 \ge \dots \ge \lambda_n$$

Miroslav Fiedler proved that if $\sigma$ is a Suleĭmanova spectrum, then $\sigma \in \operatorname{SNIEP}(n)$ if and only if the trace is nonnegative: $\sum_{i=1}^n \lambda_i \ge 0$. Furthermore, Fiedler provided an explicit recursive orthogonal similarity construction to build $A = A^T \ge 0$.

Soules Matrices & Step Bases (Soules, 1983)

Orthogonal Basis

George W. Soules introduced a powerful class of orthogonal matrices $R \in O(n)$ constructed from binary tree decompositions. If $R$ is a Soules matrix (its first column is strictly positive $r_1 > 0$, and subsequent columns possess specific step sign patterns), then for any spectrum $\sigma = (\lambda_1, \dots, \lambda_n)$ with $\lambda_1 \ge \dots \ge \lambda_n$ satisfying monotonic step bounds:

$$A = R \operatorname{diag}(\lambda_1, \dots, \lambda_n) R^T \ge 0$$

Soules bases remain the premier practical method for generating continuous families of symmetric nonnegative matrices.

4. The Symmetry Gap: $\operatorname{SNIEP}(n) \subsetneq \operatorname{RNIEP}(n)$

Why orthogonal eigenvectors fail where general non-symmetric realizations succeed.

For $n \le 4$, the Symmetric and Real NIEP coincide ($\operatorname{SNIEP}(n) = \operatorname{RNIEP}(n)$). In 1996, Charles R. Johnson, Thomas J. Laffey, and Raphael Loewy proved that the two problems diverge (first demonstrated at $n = 6$), and subsequent results by Meehan (1998), Laffey & Meehan (1999), and Hartwig & Loewy established $\operatorname{SNIEP}(n) \subsetneq \operatorname{RNIEP}(n)$ for every $n \ge 5$.

Order-5 Separating Spectrum

Hartwig–Loewy · Laffey–Meehan (1999)

In the trace-zero benchmark family $\sigma(t) = \{3+t, 3-t, -2, -2, -2\}$, taking $t = \tfrac{1}{2}$ and scaling by $2$ yields:

$$\sigma = \{7, 5, -4, -4, -4\}$$

Its power sums satisfy $s_1 = 0$, $s_2 = 122$, $s_3 = 276$, $s_4 = 3794$, $s_5 = 16860 > 0$, and $4s_4 - s_2^2 = 292 > 0$, and $\sigma$ is realized by an explicit $5 \times 5$ nonnegative integer matrix (see RNIEP). Yet no symmetric $5 \times 5$ nonnegative matrix can realize $\sigma$:

$$\{7, 5, -4, -4, -4\} \in \operatorname{RNIEP}(5) \setminus \operatorname{SNIEP}(5)$$

Mechanism: Any symmetric $A = A^T \ge 0$ of order $5$ with $\operatorname{tr}(A) = 0$ has zero diagonal ($a_{ii} = 0$), which forces its ordered eigenvalues $\lambda_1 \ge \dots \ge \lambda_5$ to obey $\lambda_2 + \lambda_5 \le 0$ (McDonald & Neumann 2000; Spector 2009; hence $\sigma(t) \in \operatorname{SNIEP}(5) \iff t \ge 1$). For $\{7, 5, -4, -4, -4\}$, $\lambda_2 + \lambda_5 = 5 + (-4) = 1 > 0$, ruling out any symmetric nonnegative realization.

5. Convex Geometry & Trace Polytopes

Visualizing the feasible spectral region as convex polyhedral slices.

While the full moment image in power-sum coordinates $(s_1, s_2, \dots, s_k)$ is a curved semialgebraic region, fixing an eigenvector frame (such as a Soules basis or character-table Perron similarity) or slicing along normalized spectral families produces convex polyhedral cross-sections and semialgebraic trace projections.

Because $\operatorname{SNIEP}(n) \subsetneq \operatorname{RNIEP}(n)$ for $n \ge 5$, the realizable symmetric region sits strictly inside the RNIEP region, giving researchers an interactive geometric window into the exact boundary where symmetry fails.

Interactive 3D Trace Polytope Gallery

Explore interactive 2D and 3D WebGL projections comparing order-5 SNIEP vs. RNIEP trace polytopes and realizability regions.

Open Polytope Gallery →

6. Primary Literature & Citations

Foundational papers from the NIEP Comprehensive Bibliography

7. See Also & Related Articles