NIEP
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 are identical: any real spectrum realizable non-symmetrically is also realizable symmetrically. However, in 1996, Charles Johnson, Thomas Laffey, and Raphael Loewy proved that for all $n \ge 5$, the SNIEP is strictly smaller than the RNIEP.

The Johnson-Laffey-Loewy Spectrum

Separation Landmark (1996)

Consider the multiset of order 5:

$$\sigma = \{4, 2, 2, -4, -4\}$$

Notice that $\sum \lambda_i = 4 + 2 + 2 - 4 - 4 = 0$ and all power sum traces are positive. It can be realized by an asymmetric nonnegative matrix $A \ge 0$ (hence $\sigma \in \operatorname{RNIEP}(5)$). However, Johnson, Laffey, and Loewy proved that no symmetric $5 \times 5$ nonnegative matrix can have spectrum $\sigma$:

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

Mechanism: If $A = A^T \ge 0$ realized $\sigma$, the trace of $A$ is zero, forcing all diagonal entries to be zero: $a_{ii} = 0$. Since $A$ is symmetric, this forces specific zero sub-block patterns that constrain the 2-dimensional eigenspace corresponding to $-4$, leading to an algebraic contradiction.

5. Convex Geometry & Trace Polytopes

Visualizing the feasible spectral region as convex polyhedral slices.

When candidate spectra are parameterized by power sum coordinates $(s_1, s_2, \dots, s_k)$, the set of realizable symmetric spectra forms a convex polyhedron called the SNIEP Trace Polytope.

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

Interactive 3D Trace Polytope Gallery

Explore interactive 3D WebGL projections comparing SNIEP vs. RNIEP trace polytopes for dimensions $n = 4, 5, 6$.

Open Polytope Gallery

6. See Also & Related Articles

← Return to NIEP Research Hub