Matrix powers and linear recurrences

Nonnegative powers of rational matrices

A published 2024 theorem computes the complete set of exponents at which a rational matrix has no negative entry. It decides both whether one positive power works and whether every sufficiently large power works. Those are different questions, and the distinction matters mathematically and algorithmically.

FPRD-MN-D01 · definitions

Three questions with different quantifiers

Let M∈Qd×dM\in\mathbb Q^{d\times d} and read Mn≥0M^n\ge0 entrywise. The three standard questions are

PN(M): ∃n≥1Mn≥0,ENN(M): ∃N ∀n≥NMn≥0,AN(M): ∀n≥0Mn≥0. \begin{aligned} \mathsf{PN}(M)&:\ \exists n\ge1\quad M^n\ge0,\\ \mathsf{ENN}(M)&:\ \exists N\ \forall n\ge N\quad M^n\ge0,\\ \mathsf{AN}(M)&:\ \forall n\ge0\quad M^n\ge0. \end{aligned}

The first asks for one positive exponent; the second asks for an entire tail. If the first allowed n=0n=0, every matrix would satisfy it because M0=I≥0M^0=I\ge0. Keeping the exponent domain explicit prevents the existential, eventual, and everywhere variants from being conflated.

FPRD-MN-T01 · published theorem

The good exponents are effectively periodic

Theorem 7 (D'Costa–Ouaknine–Worrell). For every rational square matrix MM, the setS(M)={n∈N:Mn≥0}S(M)=\{n\in\mathbb N:M^n\ge0\}is effectively computable and ultimately periodic.

The proof first chooses a modulus that separates root-of-unity collisions among eigenvalues. Each residue class then gives nondegenerate linear recurrence sequences for the matrix entries. Here dominated means that an entry recurrence has a unique characteristic root of greatest modulus. If every entry sequence in a residue is dominated, their eventual signs are computable and can be intersected.

The load-bearing step is global. Ifn0≡ℓ(modL)n_0\equiv \ell\pmod L andMn0≥0M^{n_0}\ge0, then the powersMn0(1+Lt)M^{n_0(1+Lt)} remain in the same residue and have nonnegative matrix realizations. Applied entry by entry, Berstel's theorem (or Perron–Frobenius theory) forces those nondegenerate subsequences, and hence the original residue sequences, to be dominated. Therefore a single nondominated entry makes the entire residue empty.

This one-matrix theorem is unconditional. Schanuel's conjecture is used for the paper's broader theorem about commuting matrix semigroups, not for this result.

FPRD-MN-C01 · corollary

Both natural readings are decidable

An effective ultimately periodic set has a finite exceptional prefix and a periodic tail. Consequently,

  • power nonnegativity asks whether S(M)S(M) contains a positive exponent;
  • eventual nonnegativity asks whether S(M)S(M) is cofinite.

Both tests are decidable from the published construction. This is an application of an existing theorem, not a claim that FPRD Lab independently closed the decision problem.

FPRD-MN-C02 · counterexample

One good power does not give a good tail

For the scalar matrix M=[−1]M=[-1], every even power is [1][1] and every odd power is [−1][-1]. Thus one positive nonnegative power exists, but the powers are not eventually nonnegative.

More generally, if Mk≥0M^k\ge0, then every multiple Mrk=(Mk)rM^{rk}=(M^k)^r is nonnegative. This gives an arithmetic progression, not necessarily a cofinite tail.

FPRD-MN-R01 · reduction

All entries can be interleaved into one recurrence

Each sequence uij(n)=(Mn)iju_{ij}(n)=(M^n)_{ij}is an LRS of order at most dd. Number the s=d2s=d^2 entries byr=0,…,s−1r=0,\ldots,s-1 and definevsn+r=ur(n).v_{sn+r}=u_r(n).

If Ur(x)=Nr(x)/D(x)U_r(x)=N_r(x)/D(x), thenV(x)=∑r=0s−1xrUr(xs)=∑rxrNr(xs)D(xs). V(x)=\sum_{r=0}^{s-1}x^rU_r(x^s) =\frac{\sum_r x^rN_r(x^s)}{D(x^s)}. Hence vv has order at mostd3d^3, and its ultimate nonnegativity is equivalent to eventual nonnegativity of the whole matrix sequence.

This reduction is valid but not decisive: general Ultimate Positivity for one LRS remains open. The matrix theorem succeeds by using structure shared across all entries.

FPRD-MN-F01 · counterfamily

A polynomial cutoff misses valid instances

Take a signed cycle of odd prime length ppwith one negative edge. Its first nonnegative power is2p2p. For block lengthsp1,…,pkp_1,\ldots,p_k, the block diagonal sum first becomes nonnegative at

2lcm⁡(p1,…,pk)=2∏i=1kpi. 2\operatorname{lcm}(p_1,\ldots,p_k)=2\prod_{i=1}^k p_i.

Choosing the first odd primes gives dimensiond=∑ipid=\sum_i p_i and a least good exponent exp⁡(Θ(dlog⁡d))\exp(\Theta(\sqrt{d\log d})). Thus no dimension-polynomial forward search is sound.

The exponent still has a short binary description. This family rules out one algorithmic strategy; it is not a decision-problem complexity lower bound.

Independent finite checks

A dependency-free exact checker verified 3,564 Cayley–Hamilton entry identities and 4,212 interleaved values for all 81 two-by-two matrices over {−1,0,1}\{-1,0,1\}. It also enumerated 4,282 signed permutation matrices through dimension five and reproduced the first-hit value2p2p for signed cycles of lengths 3, 5, 7, and 11. These computations test the reductions and hostile examples; the published spectral proof establishes the unrestricted theorem. Read the output.

Sources

The decision theorem is Julian D'Costa, Joël Ouaknine, and James Worrell, Nonnegativity Problems for Matrix Semigroups, LIPIcs 289 (STACS 2024), article 27, pages 27:1–27:16, DOI 10.4230/LIPIcs.STACS.2024.27. Theorem 7 supplies effective ultimate periodicity; Remark 8 isolates the positive-power test.

The standard eventual quantifiers and the earlier reduction are from S. Akshay, Supratik Chakraborty, and Debtanu Pal, On Eventual Non-negativity and Positivity for the Weighted Sum of Powers of Matrices, IJCAR 2022, LNCS 13385, pages 671–690, DOI 10.1007/978-3-031-10769-6_39, arXiv:2205.09190. Definition 1 uses a universal tail; it does not support a one-exponent reading.

Problem source. Eventual non-negativity of matrices, contributed by S. Akshay. The source is cited here for provenance; it is not a research category.